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.

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.

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.

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.

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.

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.

Structural complexity and molecular heterogeneity of a butterfly ejaculate reflect a complex history of selection.

Male ejaculates are often structurally complex, and this complexity is likely to influence key reproductive interactions between males and females. However, despite its potential evolutionary significance, the molecular underpinnings of ejaculate structural complexity have received little empirical attention. To address this knowledge gap, we sought to understand the biochemical and functional properties of the structurally complex ejaculates of Pieris rapae butterflies. Males in this species produce large ejaculates called spermatophores composed of an outer envelope, an inner matrix, and a bolus of sperm. Females are thought to benefit from the nutrition contained in the soluble inner matrix through increases in longevity and fecundity. However, the indigestible outer envelope of the spermatophore delays female remating, allowing males to monopolize paternity for longer. Here, we show that these two nonsperm-containing spermatophore regions, the inner matrix and the outer envelope, differ in their protein composition and functional properties. We also reveal how these divergent protein mixtures are separately stored in the male reproductive tract and sequentially transferred to the female reproductive tract during spermatophore assembly. Intriguingly, we discovered large quantities of female-derived proteases in both spermatophore regions shortly after mating, which may contribute to spermatophore digestion and hence, female control over remating rate. Finally, we report evidence of past selection on these spermatophore proteins and female proteases, indicating a complex evolutionary history. Our findings illustrate how structural complexity of ejaculates may allow functionally and/or spatially associated suites of proteins to respond rapidly to divergent selective pressures, such as sexual conflict or reproductive cooperation.

Evolutionary signatures amongst disease genes permit novel methods for gene prioritization and construction of informative gene-based networks.

Genes involved in the same function tend to have similar evolutionary histories, in that their rates of evolution covary over time. This coevolutionary signature, termed Evolutionary Rate Covariation (ERC), is calculated using only gene sequences from a set of closely related species and has demonstrated potential as a computational tool for inferring functional relationships between genes. To further define applications of ERC, we first established that roughly 55% of genetic diseases posses an ERC signature between their contributing genes. At a false discovery rate of 5% we report 40 such diseases including cancers, developmental disorders and mitochondrial diseases. Given these coevolutionary signatures between disease genes, we then assessed ERC's ability to prioritize known disease genes out of a list of unrelated candidates. We found that in the presence of an ERC signature, the true disease gene is effectively prioritized to the top 6% of candidates on average. We then apply this strategy to a melanoma-associated region on chromosome 1 and identify MCL1 as a potential causative gene. Furthermore, to gain global insight into disease mechanisms, we used ERC to predict molecular connections between 310 nominally distinct diseases. The resulting "disease map" network associates several diseases with related pathogenic mechanisms and unveils many novel relationships between clinically distinct diseases, such as between Hirschsprung's disease and melanoma. Taken together, these results demonstrate the utility of molecular evolution as a gene discovery platform and show that evolutionary signatures can be used to build informative gene-based networks.

Overlapping Patterns of Rapid Evolution in the Nucleic Acid Sensors cGAS and OAS1 Suggest a Common Mechanism of Pathogen Antagonism and Escape.

A diverse subset of pattern recognition receptors (PRRs) detects pathogen-associated nucleic acids to initiate crucial innate immune responses in host organisms. Reflecting their importance for host defense, pathogens encode various countermeasures to evade or inhibit these immune effectors. PRRs directly engaged by pathogen inhibitors often evolve under recurrent bouts of positive selection that have been described as molecular 'arms races.' Cyclic GMP-AMP synthase (cGAS) was recently identified as a key PRR. Upon binding cytoplasmic double-stranded DNA (dsDNA) from various viruses, cGAS generates the small nucleotide secondary messenger cGAMP to signal activation of innate defenses. Here we report an evolutionary history of cGAS with recurrent positive selection in the primate lineage. Recent studies indicate a high degree of structural similarity between cGAS and 2'-5'-oligoadenylate synthase 1 (OAS1), a PRR that detects double-stranded RNA (dsRNA), despite low sequence identity between the respective genes. We present comprehensive comparative evolutionary analysis of cGAS and OAS1 primate sequences and observe positive selection at nucleic acid binding interfaces and distributed throughout both genes. Our data revealed homologous regions with strong signatures of positive selection, suggesting common mechanisms employed by unknown pathogen encoded inhibitors and similar modes of evasion from antagonism. Our analysis of cGAS diversification also identified alternately spliced forms missing multiple sites under positive selection. Further analysis of selection on the OAS family in primates, which comprises OAS1, OAS2, OAS3 and OASL, suggests a hypothesis where gene duplications and domain fusion events result in paralogs that provide another means of escaping pathogen inhibitors. Together our comparative evolutionary analysis of cGAS and OAS provides new insights into distinct mechanisms by which key molecular sentinels of the innate immune system have adapted to circumvent viral-encoded inhibitors.

The Budding Yeast Ubiquitin Protease Ubp7 Is a Novel Component Involved in S Phase Progression.

DNA damage must be repaired in an accurate and timely fashion to preserve genome stability. Cellular mechanisms preventing genome instability are crucial to human health because genome instability is considered a hallmark of cancer. Collectively referred to as the DNA damage response, conserved pathways ensure proper DNA damage recognition and repair. The function of numerous DNA damage response components is fine-tuned by posttranslational modifications, including ubiquitination. This not only involves the enzyme cascade responsible for conjugating ubiquitin to substrates but also requires enzymes that mediate directed removal of ubiquitin. Deubiquitinases remove ubiquitin from substrates to prevent degradation or to mediate signaling functions. The Saccharomyces cerevisiae deubiquitinase Ubp7 has been characterized previously as an endocytic factor. However, here we identify Ubp7 as a novel factor affecting S phase progression after hydroxyurea treatment and demonstrate an evolutionary and genetic interaction of Ubp7 with DNA damage repair pathways of homologous recombination and nucleotide excision repair. We find that deletion of UBP7 sensitizes cells to hydroxyurea and cisplatin and demonstrate that factors that stabilize replication forks are critical under these conditions. Furthermore, ubp7Δ cells exhibit an S phase progression defect upon checkpoint activation by hydroxyurea treatment. ubp7Δ mutants are epistatic to factors involved in histone maintenance and modification, and we find that a subset of Ubp7 is chromatin-associated. In summary, our results suggest that Ubp7 contributes to S phase progression by affecting the chromatin state at replication forks, and we propose histone H2B ubiquitination as a potential substrate of Ubp7.

Hundreds of Genes Experienced Convergent Shifts in Selective Pressure in Marine Mammals.

Mammal species have made the transition to the marine environment several times, and their lineages represent one of the classical examples of convergent evolution in morphological and physiological traits. Nevertheless, the genetic mechanisms of their phenotypic transition are poorly understood, and investigations into convergence at the molecular level have been inconclusive. While past studies have searched for convergent changes at specific amino acid sites, we propose an alternative strategy to identify those genes that experienced convergent changes in their selective pressures, visible as changes in evolutionary rate specifically in the marine lineages. We present evidence of widespread convergence at the gene level by identifying parallel shifts in evolutionary rate during three independent episodes of mammalian adaptation to the marine environment. Hundreds of genes accelerated their evolutionary rates in all three marine mammal lineages during their transition to aquatic life. These marine-accelerated genes are highly enriched for pathways that control recognized functional adaptations in marine mammals, including muscle physiology, lipid-metabolism, sensory systems, and skin and connective tissue. The accelerations resulted from both adaptive evolution as seen in skin and lung genes, and loss of function as in gustatory and olfactory genes. In regard to sensory systems, this finding provides further evidence that reduced senses of taste and smell are ubiquitous in marine mammals. Our analysis demonstrates the feasibility of identifying genes underlying convergent organism-level characteristics on a genome-wide scale and without prior knowledge of adaptations, and provides a powerful approach for investigating the physiological functions of mammalian genes.

Evolutionary rate covariation identifies new members of a protein network required for Drosophila melanogaster female post-mating responses.

Seminal fluid proteins transferred from males to females during copulation are required for full fertility and can exert dramatic effects on female physiology and behavior. In Drosophila melanogaster, the seminal protein sex peptide (SP) affects mated females by increasing egg production and decreasing receptivity to courtship. These behavioral changes persist for several days because SP binds to sperm that are stored in the female. SP is then gradually released, allowing it to interact with its female-expressed receptor. The binding of SP to sperm requires five additional seminal proteins, which act together in a network. Hundreds of uncharacterized male and female proteins have been identified in this species, but individually screening each protein for network function would present a logistical challenge. To prioritize the screening of these proteins for involvement in the SP network, we used a comparative genomic method to identify candidate proteins whose evolutionary rates across the Drosophila phylogeny co-vary with those of the SP network proteins. Subsequent functional testing of 18 co-varying candidates by RNA interference identified three male seminal proteins and three female reproductive tract proteins that are each required for the long-term persistence of SP responses in females. Molecular genetic analysis showed the three new male proteins are required for the transfer of other network proteins to females and for SP to become bound to sperm that are stored in mated females. The three female proteins, in contrast, act downstream of SP binding and sperm storage. These findings expand the number of seminal proteins required for SP's actions in the female and show that multiple female proteins are necessary for the SP response. Furthermore, our functional analyses demonstrate that evolutionary rate covariation is a valuable predictive tool for identifying candidate members of interacting protein networks.

Digestive organ in the female reproductive tract borrows genes from multiple organ systems to adopt critical functions.

Persistent adaptive challenges are often met with the evolution of novel physiological traits. Although there are specific examples of single genes providing new physiological functions, studies on the origin of complex organ functions are lacking. One such derived set of complex functions is found in the Lepidopteran bursa copulatrix, an organ within the female reproductive tract that digests nutrients from the male ejaculate or spermatophore. Here, we characterized bursa physiology and the evolutionary mechanisms by which it was equipped with digestive and absorptive functionality. By studying the transcriptome of the bursa and eight other tissues, we revealed a suite of highly expressed and secreted gene products providing the bursa with a combination of stomach-like traits for mechanical and enzymatic digestion of the male spermatophore. By subsequently placing these bursa genes in an evolutionary framework, we found that the vast majority of their novel digestive functions were co-opted by borrowing genes that continue to be expressed in nonreproductive tissues. However, a number of bursa-specific genes have also arisen, some of which represent unique gene families restricted to Lepidoptera and may provide novel bursa-specific functions. This pattern of promiscuous gene borrowing and relatively infrequent evolution of tissue-specific duplicates stands in contrast to studies of the evolution of novelty via single gene co-option. Our results suggest that the evolution of complex organ-level phenotypes may often be enabled (and subsequently constrained) by changes in tissue specificity that allow expression of existing genes in novel contexts, such as reproduction. The extent to which the selective pressures encountered in these novel roles require resolution via duplication and sub/neofunctionalization is likely to be determined by the need for specialized reproductive functionality. Thus, complex physiological phenotypes such as that found in the bursa offer important opportunities for understanding the relative role of pleiotropy and specialization in adaptive evolution.

ERC analysis: web-based inference of gene function via evolutionary rate covariation.

UNLABELLED: The recent explosion of comparative genomics data presents an unprecedented opportunity to construct gene networks via the evolutionary rate covariation (ERC) signature. ERC is used to identify genes that experienced similar evolutionary histories, and thereby draws functional associations between them. The ERC Analysis website allows researchers to exploit genome-wide datasets to infer novel genes in any biological function and to explore deep evolutionary connections between distinct pathways and complexes. The website provides five analytical methods, graphical output, statistical support and access to an increasing number of taxonomic groups. AVAILABILITY AND IMPLEMENTATION: Analyses and data at http://csb.pitt.edu/erc_analysis/ CONTACT: nclark@pitt.edu.

Evolutionary rate covariation reveals shared functionality and coexpression of genes.

Evolutionary rate covariation (ERC) is a phylogenetic signature that reflects the covariation of a pair of proteins over evolutionary time. ERC is typically elevated between interacting proteins and so is a promising signature to characterize molecular and functional interactions across the genome. ERC is often assumed to result from compensatory changes at interaction interfaces (i.e., intermolecular coevolution); however, its origin is still unclear and is likely to be complex. Here, we determine the biological factors responsible for ERC in a proteome-wide data set of 4459 proteins in 18 budding yeast species. We show that direct physical interaction is not required to produce ERC, because we observe strong correlations between noninteracting but cofunctional enzymes. We also demonstrate that ERC is uniformly distributed along the protein primary sequence, suggesting that intermolecular coevolution is not generally responsible for ERC between physically interacting proteins. Using multivariate analysis, we show that a pair of proteins is likely to exhibit ERC if they share a biological function or if their expression levels coevolve between species. Thus, ERC indicates shared function and coexpression of protein pairs and not necessarily coevolution between sites, as has been assumed in previous studies. This full interpretation of ERC now provides us with a powerful tool to assign uncharacterized proteins to functional groups and to determine the interconnectedness between entire genetic pathways.

Dynamic digestive physiology of a female reproductive organ in a polyandrous butterfly.

Reproductive traits experience high levels of selection because of their direct ties to fitness, often resulting in rapid adaptive evolution. Much of the work in this area has focused on male reproductive traits. However, a more comprehensive understanding of female reproductive adaptations and their relationship to male characters is crucial to uncover the relative roles of sexual cooperation and conflict in driving co-evolutionary dynamics between the sexes. We focus on the physiology of a complex female reproductive adaptation in butterflies and moths: a stomach-like organ in the female reproductive tract called the bursa copulatrix that digests the male ejaculate (spermatophore). Little is known about how the bursa digests the spermatophore. We characterized bursa proteolytic capacity in relation to female state in the polyandrous butterfly Pieris rapae. We found that the virgin bursa exhibits extremely high levels of proteolytic activity. Furthermore, in virgin females, bursal proteolytic capacity increases with time since eclosion and ambient temperature, but is not sensitive to the pre-mating social environment. Post copulation, bursal proteolytic activity decreases rapidly before rebounding toward the end of a mating cycle, suggesting active female regulation of proteolysis and/or potential quenching of proteolysis by male ejaculate constituents. Using transcriptomic and proteomic approaches, we report identities for nine proteases actively transcribed by bursal tissue and/or expressed in the bursal lumen that may contribute to observed bursal proteolysis. We discuss how these dynamic physiological characteristics may function as female adaptations resulting from sexual conflict over female remating rate in this polyandrous butterfly.

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.

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.

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.

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.