Vocal learning-associated convergent evolution in mammalian proteins and regulatory elements.

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

Evolutionary rate covariation is a reliable predictor of co-functional interactions but not necessarily physical interactions.

Co-functional proteins tend to have rates of evolution that covary over time. This correlation between evolutionary rates can be measured over the branches of a phylogenetic tree through methods such as evolutionary rate covariation (ERC), and then used to construct gene networks by the identification of proteins with functional interactions. The cause of this correlation has been hypothesized to result from both compensatory coevolution at physical interfaces and nonphysical forces such as shared changes in selective pressure. This study explores whether coevolution due to compensatory mutations has a measurable effect on the ERC signal. We examined the difference in ERC signal between physically interacting protein domains within complexes compared to domains of the same proteins that do not physically interact. We found no generalizable relationship between physical interaction and high ERC, although a few complexes ranked physical interactions higher than nonphysical interactions. Therefore, we conclude that coevolution due to physical interaction is weak, but present in the signal captured by ERC, and we hypothesize that the stronger signal instead comes from selective pressures on the protein as a whole and maintenance of the general function.

Molecular evolution of male reproduction across species with highly divergent sperm morphology in diverse murine rodents.

Sperm competition can drive rapid evolution of male reproductive traits, but it remains unclear how variation in sperm competition intensity shapes phenotypic and molecular diversity across clades. Old World mice and rats (subfamily Murinae) comprise a rapid radiation and exhibit incredible diversity in sperm morphology and production. We combined phenotype and sequence data to model the evolution of reproductive traits and genes across 78 murine species. We identified several shifts towards smaller relative testes mass, a trait reflective of reduced sperm competition. Several sperm traits were associated with relative testes mass, suggesting that mating system evolution likely selects for convergent traits related to sperm competitive ability. Molecular evolutionary rates of spermatogenesis proteins also correlated with relative testes mass, but in an unexpected direction. We predicted that sperm competition would result in rapid divergence among species with large relative testes mass, but instead found that many spermatogenesis genes evolve more rapidly in species with smaller relative testes mass due to relaxed purifying selection. While some reproductive genes evolved under positive selection, relaxed selection played a greater role underlying rapid evolution in small testes species. Our work demonstrates that sexual selection can impose strong purifying selection shaping the evolution of male reproduction.

Gadusol is a maternally provided sunscreen that protects fish embryos from DNA damage.

Exposure to ultraviolet radiation (UVR) is harmful to living cells, leading organisms to evolve protective mechanisms against UVR-induced cellular damage and stress.1,2 UVR, particularly UVB (280-320 nm), can damage proteins and DNA, leading to errors during DNA repair and replication. Excessive UVR can induce cellular death. Aquatic organisms face risk of UV exposure as biologically harmful levels of UVB can penetrate >10 m in clear water.3 While melanin is the only known sunscreen in vertebrates, it often emerges late in embryonic development, rendering embryos of many species vulnerable during the earlier stages. Algae and microbes produce a class of sunscreening compounds known as mycosporine-like amino acids (MAAs).4 Fish eggs contain a similar compound called gadusol, whose role as a sunscreen has yet to be tested despite its discovery over 40 years ago.5 The recent finding that many vertebrate genomes contain a biosynthetic pathway for gadusol suggests that many fish may produce and use this molecule as a sunscreen.6 We generated a gadusol-deficient mutant zebrafish to investigate the role of gadusol in protecting fish embryos and larvae from UVR. Our results demonstrate that maternally provided gadusol is the primary sunscreen in embryonic and larval development, while melanin provides modest secondary protection. The gadusol biosynthetic pathway is retained in the vast majority of teleost genomes but is repeatedly lost in species whose young are no longer exposed to UVR. Our data demonstrate that gadusol is a maternally provided sunscreen that is critical for early-life survival in the most species-rich branch of the vertebrate phylogeny.

Reduction of Paraoxonase Expression Followed by Inactivation across Independent Semiaquatic Mammals Suggests Stepwise Path to Pseudogenization.

Convergent adaptation to the same environment by multiple lineages frequently involves rapid evolutionary change at the same genes, implicating these genes as important for environmental adaptation. Such adaptive molecular changes may yield either change or loss of protein function; loss of function can eliminate newly deleterious proteins or reduce energy necessary for protein production. We previously found a striking case of recurrent pseudogenization of the Paraoxonase 1 (Pon1) gene among aquatic mammal lineages-Pon1 became a pseudogene with genetic lesions, such as stop codons and frameshifts, at least four times independently in aquatic and semiaquatic mammals. Here, we assess the landscape and pace of pseudogenization by studying Pon1 sequences, expression levels, and enzymatic activity across four aquatic and semiaquatic mammal lineages: pinnipeds, cetaceans, otters, and beavers. We observe in beavers and pinnipeds an unexpected reduction in expression of Pon3, a paralog with similar expression patterns but different substrate preferences. Ultimately, in all lineages with aquatic/semiaquatic members, we find that preceding any coding-level pseudogenization events in Pon1, there is a drastic decrease in expression, followed by relaxed selection, thus allowing accumulation of disrupting mutations. The recurrent loss of Pon1 function in aquatic/semiaquatic lineages is consistent with a benefit to Pon1 functional loss in aquatic environments. Accordingly, we examine diving and dietary traits across pinniped species as potential driving forces of Pon1 functional loss. We find that loss is best associated with diving activity and likely results from changes in selective pressures associated with hypoxia and hypoxia-induced inflammation.

Gadusol is a maternally provided sunscreen that protects fish embryos from DNA damage.

Ultraviolet radiation (UVR) and its deleterious effects on living cells selects for UVR-protective mechanisms. Organisms across the tree of life evolved a variety of natural sunscreens to prevent UVR-induced cellular damage and stress. However, in vertebrates, only melanin is known to act as a sunscreen. Here we demonstrate that gadusol, a transparent compound discovered over 40 years ago in fish eggs, is a maternally provided sunscreen required for survival of embryonic and larval zebrafish exposed to UVR. Mutating an enzyme involved in gadusol biosynthesis increases the formation of cyclobutane pyrimidine dimers, a hallmark of UVB-induced DNA damage. Compared to the contributions of melanin and the chorion, gadusol is the primary sunscreening mechanism in embryonic and larval fish. The gadusol biosynthetic pathway is retained in the vast majority of teleost genomes but is repeatedly lost in species whose young are no longer exposed to UVR. Our data demonstrate that gadusol is a maternally provided sunscreen that is critical for early-life survival in the most species-rich branch of the vertebrate phylogeny.

Highly Dynamic Gene Family Evolution Suggests Changing Roles for PON Genes Within Metazoa.

Change in gene family size has been shown to facilitate adaptation to different selective pressures. This includes gene duplication to increase dosage or diversification of enzymatic substrates and gene deletion due to relaxed selection. We recently found that the PON1 gene, an enzyme with arylesterase and lactonase activity, was lost repeatedly in different aquatic mammalian lineages, suggesting that the PON gene family is responsive to environmental change. We further investigated if these fluctuations in gene family size were restricted to mammals and approximately when this gene family was expanded within mammals. Using 112 metazoan protein models, we explored the evolutionary history of the PON family to characterize the dynamic evolution of this gene family. We found that there have been multiple, independent expansion events in tardigrades, cephalochordates, and echinoderms. In addition, there have been partial gene loss events in monotremes and sea cucumbers and what appears to be complete loss in arthropods, urochordates, platyhelminths, ctenophores, and placozoans. In addition, we show the mammalian expansion to three PON paralogs occurred in the ancestor of all mammals after the divergence of sauropsida but before the divergence of monotremes from therians. We also provide evidence of a novel PON expansion within the brushtail possum. In the face of repeated expansions and deletions in the context of changing environments, we suggest a range of selective pressures, including pathogen infection and mitigation of oxidative damage, are likely influencing the diversification of this dynamic gene family across metazoa.

Multiple 9-1-1 complexes promote homolog synapsis, DSB repair, and ATR signaling during mammalian meiosis.

DNA damage response mechanisms have meiotic roles that ensure successful gamete formation. While completion of meiotic double-strand break (DSB) repair requires the canonical RAD9A-RAD1-HUS1 (9A-1-1) complex, mammalian meiocytes also express RAD9A and HUS1 paralogs, RAD9B and HUS1B, predicted to form alternative 9-1-1 complexes. The RAD1 subunit is shared by all predicted 9-1-1 complexes and localizes to meiotic chromosomes even in the absence of HUS1 and RAD9A. Here, we report that testis-specific disruption of RAD1 in mice resulted in impaired DSB repair, germ cell depletion, and infertility. Unlike Hus1 or Rad9a disruption, Rad1 loss in meiocytes also caused severe defects in homolog synapsis, impaired phosphorylation of ATR targets such as H2AX, CHK1, and HORMAD2, and compromised meiotic sex chromosome inactivation. Together, these results establish critical roles for both canonical and alternative 9-1-1 complexes in meiotic ATR activation and successful prophase I completion.

Activation by cleavage of the epithelial Na+ channel α and γ subunits independently coevolved with the vertebrate terrestrial migration.

Vertebrates evolved mechanisms for sodium conservation and gas exchange in conjunction with migration from aquatic to terrestrial habitats. Epithelial Na+ channel (ENaC) function is critical to systems responsible for extracellular fluid homeostasis and gas exchange. ENaC is activated by cleavage at multiple specific extracellular polybasic sites, releasing inhibitory tracts from the channel's α and γ subunits. We found that proximal and distal polybasic tracts in ENaC subunits coevolved, consistent with the dual cleavage requirement for activation observed in mammals. Polybasic tract pairs evolved with the terrestrial migration and the appearance of lungs, coincident with the ENaC activator aldosterone, and appeared independently in the α and γ subunits. In summary, sites within ENaC for protease activation developed in vertebrates when renal Na+ conservation and alveolar gas exchange were required for terrestrial survival.

Evolutionary, proteomic, and experimental investigations suggest the extracellular matrix of cumulus cells mediates fertilization outcomes†.

Studies of fertilization biology often focus on sperm and egg interactions. However, before gametes interact, mammalian sperm must pass through the cumulus layer; in mice, this consists of several thousand cells tightly glued together with hyaluronic acid and other proteins. To better understand the role of cumulus cells and their extracellular matrix, we perform proteomic experiments on cumulus oophorus complexes (COCs) in house mice (Mus musculus), producing over 24,000 mass spectra to identify 711 proteins. Seven proteins known to stabilize hyaluronic acid and the extracellular matrix were especially abundant (using spectral counts as an indirect proxy for abundance). Through comparative evolutionary analyses, we show that three of these evolve rapidly, a classic signature of genes that influence fertilization rate. Some of the selected sites overlap regions of the protein known to impact function. In a follow-up experiment, we compared COCs from females raised in two different social environments. Female mice raised in the presence of multiple males produced COCs that were smaller and more resistant to dissociation by hyaluronidase compared to females raised in the presence of a single male, consistent with a previous study that demonstrated such females produced COCs that were more resistant to fertilization. Although cumulus cells are often thought of as enhancers of fertilization, our evolutionary, proteomic, and experimental investigations implicate their extracellular matrix as a potential mediator of fertilization outcomes.

Experimental exchange of paralogous domains in the MLH family provides evidence of sub-functionalization after gene duplication.

Baker's yeast contains a large number of duplicated genes; some function redundantly, whereas others have more specialized roles. We used the MLH family of DNA mismatch repair (MMR) proteins as a model to better understand the steps that lead to gene specialization following a gene duplication event. We focused on two highly conserved yeast MLH proteins, Pms1 and Mlh3, with Pms1 having a major role in the repair of misincorporation events during DNA replication and Mlh3 acting to resolve recombination intermediates in meiosis to form crossovers. The baker's yeast Mlh3 and Pms1 proteins are significantly diverged (19% overall identity), suggesting that an extensive number of evolutionary steps, some major, others involving subtle refinements, took place to diversify the MLH proteins. Using phylogenetic and molecular approaches, we provide evidence that all three domains (N-terminal ATP binding, linker, C-terminal endonuclease/MLH interaction) in the MLH protein family are critical for conferring pathway specificity. Importantly, mlh3 alleles in the ATP binding and endonuclease domains improved MMR functions in strains lacking the Pms1 protein and did not disrupt Mlh3 meiotic functions. This ability for mlh3 alleles to complement the loss of Pms1 suggests that an ancestral Pms1/Mlh3 protein was capable of performing both MMR and crossover functions. Our strategy for analyzing MLH pathway specificity provides an approach to understand how paralogs have evolved to support distinct cellular processes.

A Drosophila screen identifies NKCC1 as a modifier of NGLY1 deficiency.

N-Glycanase 1 (NGLY1) is a cytoplasmic deglycosylating enzyme. Loss-of-function mutations in the NGLY1 gene cause NGLY1 deficiency, which is characterized by developmental delay, seizures, and a lack of sweat and tears. To model the phenotypic variability observed among patients, we crossed a Drosophila model of NGLY1 deficiency onto a panel of genetically diverse strains. The resulting progeny showed a phenotypic spectrum from 0 to 100% lethality. Association analysis on the lethality phenotype, as well as an evolutionary rate covariation analysis, generated lists of modifying genes, providing insight into NGLY1 function and disease. The top association hit was Ncc69 (human NKCC1/2), a conserved ion transporter. Analyses in NGLY1-/- mouse cells demonstrated that NKCC1 has an altered average molecular weight and reduced function. The misregulation of this ion transporter may explain the observed defects in secretory epithelium function in NGLY1 deficiency patients.

MCM8IP activates the MCM8-9 helicase to promote DNA synthesis and homologous recombination upon DNA damage.

Homologous recombination (HR) mediates the error-free repair of DNA double-strand breaks to maintain genomic stability. Here we characterize C17orf53/MCM8IP, an OB-fold containing protein that binds ssDNA, as a DNA repair factor involved in HR. MCM8IP-deficient cells exhibit HR defects, especially in long-tract gene conversion, occurring downstream of RAD51 loading, consistent with a role for MCM8IP in HR-dependent DNA synthesis. Moreover, loss of MCM8IP confers cellular sensitivity to crosslinking agents and PARP inhibition. Importantly, we report that MCM8IP directly associates with MCM8-9, a helicase complex mutated in primary ovarian insufficiency, and RPA1. We additionally show that the interactions of MCM8IP with MCM8-9 and RPA facilitate HR and promote replication fork progression and cellular viability in response to treatment with crosslinking agents. Mechanistically, MCM8IP stimulates the helicase activity of MCM8-9. Collectively, our work identifies MCM8IP as a key regulator of MCM8-9-dependent DNA synthesis during DNA recombination and replication.

Pan-mammalian analysis of molecular constraints underlying extended lifespan.

Although lifespan in mammals varies over 100-fold, the precise evolutionary mechanisms underlying variation in longevity remain unknown. Species-specific genetic changes have been observed in long-lived species including the naked mole-rat, bats, and the bowhead whale, but these adaptations do not generalize to other mammals. We present a novel method to identify associations between rates of protein evolution and continuous phenotypes across the entire mammalian phylogeny. Unlike previous analyses that focused on individual species, we treat absolute and relative longevity as quantitative traits and demonstrate that these lifespan traits affect the evolutionary constraint on hundreds of genes. Specifically, we find that genes related to cell cycle, DNA repair, cell death, the IGF1 pathway, and immunity are under increased evolutionary constraint in large and long-lived mammals. For mammals exceptionally long-lived for their body size, we find increased constraint in inflammation, DNA repair, and NFKB-related pathways. Strikingly, these pathways have considerable overlap with those that have been previously reported to have potentially adaptive changes in single-species studies, and thus would be expected to show decreased constraint in our analysis. This unexpected finding of increased constraint in many longevity-associated pathways underscores the power of our quantitative approach to detect patterns that generalize across the mammalian phylogeny.

Characterization of Female Reproductive Proteases in a Butterfly from Functional and Evolutionary Perspectives.

Molecules that mediate reproductive interactions are some of the most rapidly evolving traits. Researchers have often suggested that this is due to coevolution at key physiological interfaces. However, very few of these interfaces are well understood at the functional level. One such interface is the digestion of the spermatophore in Lepidoptera. Female Lepidoptera have a specialized reproductive organ called the bursa copulatrix that receives and processes the male spermatophore, a complex proteinaceous ejaculate. In the cabbage white butterfly, Pieris rapae, the bursa secretes a mixture of proteases hypothesized to digest the spermatophore. However, these proteases remain biochemically uncharacterized. Using a zymogram approach, we identified six proteases in bursal extracts at sufficiently high concentrations to characterize their in vitro activity. We assessed the modes of action of these bursal enzymes by quantifying their activity following exposure to diagnostic protease inhibitors. A serine protease-specific inhibitor failed to reduce bursal protease digestion of casein. However, a cysteine protease-specific inhibitor did decrease the activity of some proteases. To explore the possible molecular mechanisms responsible for these responses, we created protease homology models. The models mirrored the results of our in vitro experiments, indicating that protease homology models may offer insight into underlying functional mechanisms. Whether the observed bursal protease resistance to known inhibitors is important in the context of spermatophore digestion remains to be tested. However, our results suggest the exciting possibility that bursal protease specificity may have evolved in response to interactions with various proteins and inhibitors present within the female tract during the reproductive process.

Evolution-based screening enables genome-wide prioritization and discovery of DNA repair genes.

DNA repair is critical for genome stability and is maintained through conserved pathways. Traditional genome-wide mammalian screens are both expensive and laborious. However, computational approaches circumvent these limitations and are a powerful tool to identify new DNA repair factors. By analyzing the evolutionary relationships between genes in the major DNA repair pathways, we uncovered functional relationships between individual genes and identified partners. Here we ranked 17,487 mammalian genes for coevolution with 6 distinct DNA repair pathways. Direct comparison to genetic screens for homologous recombination or Fanconi anemia factors indicates that our evolution-based screen is comparable, if not superior, to traditional screening approaches. Demonstrating the utility of our strategy, we identify a role for the DNA damage-induced apoptosis suppressor (DDIAS) gene in double-strand break repair based on its coevolution with homologous recombination. DDIAS knockdown results in DNA double-strand breaks, indicated by ATM kinase activation and 53BP1 foci induction. Additionally, DDIAS-depleted cells are deficient for homologous recombination. Our results reveal that evolutionary analysis is a powerful tool to uncover novel factors and functional relationships in DNA repair.

RERconverge: an R package for associating evolutionary rates with convergent traits.

MOTIVATION: When different lineages of organisms independently adapt to similar environments, selection often acts repeatedly upon the same genes, leading to signatures of convergent evolutionary rate shifts at these genes. With the increasing availability of genome sequences for organisms displaying a variety of convergent traits, the ability to identify genes with such convergent rate signatures would enable new insights into the molecular basis of these traits. RESULTS: Here we present the R package RERconverge, which tests for association between relative evolutionary rates of genes and the evolution of traits across a phylogeny. RERconverge can perform associations with binary and continuous traits, and it contains tools for visualization and enrichment analyses of association results. AVAILABILITY AND IMPLEMENTATION: RERconverge source code, documentation and a detailed usage walk-through are freely available at https://github.com/nclark-lab/RERconverge. Datasets for mammals, Drosophila and yeast are available at https://bit.ly/2J2QBnj. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Robust Method for Detecting Convergent Shifts in Evolutionary Rates.

Identifying genomic elements underlying phenotypic adaptations is an important problem in evolutionary biology. Comparative analyses learning from convergent evolution of traits are gaining momentum in accurately detecting such elements. We previously developed a method for predicting phenotypic associations of genetic elements by contrasting patterns of sequence evolution in species showing a phenotype with those that do not. Using this method, we successfully demonstrated convergent evolutionary rate shifts in genetic elements associated with two phenotypic adaptations, namely the independent subterranean and marine transitions of terrestrial mammalian lineages. Our original method calculates gene-specific rates of evolution on branches of phylogenetic trees using linear regression. These rates represent the extent of sequence divergence on a branch after removing the expected divergence on the branch due to background factors. The rates calculated using this regression analysis exhibit an important statistical limitation, namely heteroscedasticity. We observe that the rates on branches that are longer on average show higher variance, and describe how this problem adversely affects the confidence with which we can make inferences about rate shifts. Using a combination of data transformation and weighted regression, we have developed an updated method that corrects this heteroscedasticity in the rates. We additionally illustrate the improved performance offered by the updated method at robust detection of convergent rate shifts in phylogenetic trees of protein-coding genes across mammals, as well as using simulated tree data sets. Overall, we present an important extension to our evolutionary-rates-based method that performs more robustly and consistently at detecting convergent shifts in evolutionary rates.

Evolutionary rate covariation analysis of E-cadherin identifies Raskol as a regulator of cell adhesion and actin dynamics in Drosophila.

The adherens junction couples the actin cytoskeletons of neighboring cells to provide the foundation for multicellular organization. The core of the adherens junction is the cadherin-catenin complex that arose early in the evolution of multicellularity to link actin to intercellular adhesions. Over time, evolutionary pressures have shaped the signaling and mechanical functions of the adherens junction to meet specific developmental and physiological demands. Evolutionary rate covariation (ERC) identifies proteins with correlated fluctuations in evolutionary rate that can reflect shared selective pressures and functions. Here we use ERC to identify proteins with evolutionary histories similar to the Drosophila E-cadherin (DE-cad) ortholog. Core adherens junction components α-catenin and p120-catenin displayed positive ERC correlations with DE-cad, indicating that they evolved under similar selective pressures during evolution between Drosophila species. Further analysis of the DE-cad ERC profile revealed a collection of proteins not previously associated with DE-cad function or cadherin-mediated adhesion. We then analyzed the function of a subset of ERC-identified candidates by RNAi during border cell (BC) migration and identified novel genes that function to regulate DE-cad. Among these, we found that the gene CG42684, which encodes a putative GTPase activating protein (GAP), regulates BC migration and adhesion. We named CG42684 raskol ("to split" in Russian) and show that it regulates DE-cad levels and actin protrusions in BCs. We propose that Raskol functions with DE-cad to restrict Ras/Rho signaling and help guide BC migration. Our results demonstrate that a coordinated selective pressure has shaped the adherens junction and this can be leveraged to identify novel components of the complexes and signaling pathways that regulate cadherin-mediated adhesion.

Ancient convergent losses of Paraoxonase 1 yield potential risks for modern marine mammals.

Mammals diversified by colonizing drastically different environments, with each transition yielding numerous molecular changes, including losses of protein function. Though not initially deleterious, these losses could subsequently carry deleterious pleiotropic consequences. We have used phylogenetic methods to identify convergent functional losses across independent marine mammal lineages. In one extreme case, Paraoxonase 1 (PON1) accrued lesions in all marine lineages, while remaining intact in all terrestrial mammals. These lesions coincide with PON1 enzymatic activity loss in marine species' blood plasma. This convergent loss is likely explained by parallel shifts in marine ancestors' lipid metabolism and/or bloodstream oxidative environment affecting PON1's role in fatty acid oxidation. PON1 loss also eliminates marine mammals' main defense against neurotoxicity from specific man-made organophosphorus compounds, implying potential risks in modern environments.