Search for: All records

The baculum, a bone in the penis of many mammal species, shows an astonishing level of morphological divergence between species. Despite hundreds of years of interest, biologists have been unable to directly test its function. The goal of the current study is to uncover molecular details that could allow selective disruption of the baculum while allowing normal sexual differentiation and skeletal development. We compare patterns of androgen receptor binding and single cell gene expression in the developing penis, forelimbs and hindlimbs of mice. We identified chondrocytes in all three tissue types, but those from the developing penis show several unique features, including a population of chondrocytes that express bothRunt-related transcription factor 2(Runx2) andAndrogen receptor(Ar). By combining aRunx2-Cre allele with a floxedArallele in mice, we selectively knocked out androgen signaling in late chondrocytes, resulting in a range of defects in baculum morphology. Males with the most disrupted bacula were unable to copulate, and their bacula appears to be disconnected from the corpus cavernosum muscle. Our study provides insights into the diversity of molecular mechanisms leading to bone and offers the first opportunity to directly test the function of the baculum. 

Abstract A fundamental goal in evolutionary biology and population genetics is to understand how selection shapes the fate of new mutations. Here, we test the null hypothesis that insertion–deletion (indel) events in protein-coding regions occur randomly with respect to secondary structures. We identified indels across 11,444 sequence alignments in mouse, rat, human, chimp, and dog genomes and then quantified their overlap with four different types of secondary structure—alpha helices, beta strands, protein bends, and protein turns—predicted by deep-learning methods of AlphaFold2. Indels overlapped secondary structures 54% as much as expected and were especially underrepresented over beta strands, which tend to form internal, stable regions of proteins. In contrast, indels were enriched by 155% over regions without any predicted secondary structures. These skews were stronger in the rodent lineages compared to the primate lineages, consistent with population genetic theory predicting that natural selection will be more efficient in species with larger effective population sizes. Nonsynonymous substitutions were also less common in regions of protein secondary structure, although not as strongly reduced as in indels. In a complementary analysis of thousands of human genomes, we showed that indels overlapping secondary structure segregated at significantly lower frequency than indels outside of secondary structure. Taken together, our study shows that indels are selected against if they overlap secondary structure, presumably because they disrupt the tertiary structure and function of a protein.