Nathan L. Clark

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 SPs 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.

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.

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.

Evolutionary and Functional Analysis of the Invariant SWIM Domain in the Conserved Shu2/SWS1 Protein Family from Saccharomyces cerevisiae to Homo sapiens

The Saccharomyces cerevisiae Shu2 protein is an important regulator of Rad51, which promotes homologous recombination (HR). Shu2 functions in the Shu complex with Shu1 and the Rad51 paralogs Csm2 and Psy3. Shu2 belongs to the SWS1 protein family, which is characterized by its SWIM domain (CXC...Xn...CXH), a zinc-binding motif. In humans, SWS1 interacts with the Rad51 paralog SWSAP1. Using genetic and evolutionary analyses, we examined the role of the Shu complex in mitotic and meiotic processes across eukaryotic lineages. We provide evidence that the SWS1 protein family contains orthologous genes in early-branching eukaryote lineages (e.g., Giardia lamblia), as well as in multicellular eukaryotes including Caenorhabditis elegans and Drosophila melanogaster. Using sequence analysis, we expanded the SWIM domain to include an invariant alanine three residues after the terminal CXH motif (CXC…Xn…CXHXXA). We found that the SWIM domain is conserved in all eukaryotic orthologs, and accordingly, in vivo disruption of the invariant residues within the canonical SWIM domain inhibits DNA damage tolerance in yeast and protein-protein interactions in yeast and humans. Furthermore, using evolutionary analyses, we found that yeast and Drosophila Shu2 exhibit strong coevolutionary signatures with meiotic proteins, and in yeast, its disruption leads to decreased meiotic progeny. Together our data indicate that the SWS1 family is an ancient and highly conserved eukaryotic regulator of meiotic and mitotic HR.

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.

Genetic and phenotypic influences on copulatory plug survival in mice

Across a diversity of animals, male seminal fluid coagulates upon ejaculation to form a hardened structure known as a copulatory plug. Previous studies suggest that copulatory plugs evolved as a mechanism for males to impede remating by females, but detailed investigations into the time course over which plugs survive in the female’s reproductive tract are lacking. Here, we cross males from eight inbred strains to females from two inbred strains of house mice (Mus musculus domesticus). Plug survival was significantly affected by male genotype. Against intuition, plug survival time was negatively correlated with plug size: long-lasting plugs were small and relatively more susceptible to proteolysis. Plug size was associated with divergence in major protein composition of seminal vesicle fluid, suggesting that changes in gene expression may play an important role in plug dynamics. In contrast, we found no correlation to genetic variation in the protein-coding regions of five genes thought to be important in copulatory plug formation (Tgm4, Svs1, Svs2, Svs4 and Svs5). Our study demonstrates a complex relationship between copulatory plug characteristics and survival. We discuss several models to explain unexpected variation in plug phenotypes.

Evolutionary Rate Covariation in Meiotic Proteins Results from Fluctuating Evolutionary Pressure in Yeasts and Mammals

Evolutionary rates of functionally related proteins tend to change in parallel over evolutionary time. Such evolutionary rate covariation (ERC) is a sequence-based signature of coevolution and a potentially useful signature to infer functional relationships between proteins. One major hypothesis to explain ERC is that fluctuations in evolutionary pressure acting on entire pathways cause parallel rate changes for functionally related proteins. To explore this hypothesis we analyzed ERC within DNA mismatch repair (MMR) and meiosis proteins over phylogenies of 18 yeast species and 22 mammalian species. We identified a strong signature of ERC between eight yeast proteins involved in meiotic crossing over, which seems to have resulted from relaxation of constraint specifically in Candida glabrata. These and other meiotic proteins in C. glabrata showed marked rate acceleration, likely due to its apparently clonal reproductive strategy and the resulting infrequent use of meiotic proteins. This correlation between change of reproductive mode and change in constraint supports an evolutionary pressure origin for ERC. Moreover, we present evidence for similar relaxations of constraint in additional pathogenic yeast species. Mammalian MMR and meiosis proteins also showed statistically significant ERC; however, there was not strong ERC between crossover proteins, as observed in yeasts. Rather, mammals exhibited ERC in different pathways, such as piRNA-mediated defense against transposable elements. Overall, if fluctuation in evolutionary pressure is responsible for ERC, it could reveal functional relationships within entire protein pathways, regardless of whether they physically interact or not, so long as there was variation in constraint on that pathway.

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.

Subterranean mammals show convergent regression in ocular genes and enhancers, along with adaptation to tunneling

The underground environment imposes unique demands on life that have led subterranean species to evolve specialized traits, many of which evolved convergently. We studied convergence in evolutionary rate in subterranean mammals in order to associate phenotypic evolution with specific genetic regions. We identified a strong excess of vision- and skin-related genes that changed at accelerated rates in the subterranean environment due to relaxed constraint and adaptive evolution. We also demonstrate that ocular-specific transcriptional enhancers were convergently accelerated, whereas enhancers active outside the eye were not. Furthermore, several uncharacterized genes and regulatory sequences demonstrated convergence and thus constitute novel candidate sequences for congenital ocular disorders. The strong evidence of convergence in these species indicates that evolution in this environment is recurrent and predictable and can be used to gain insights into phenotype–genotype relationships.

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

Significance Male ejaculates exhibit remarkable diversity, including variation in their spatial and temporal molecular composition. This complexity suggests that ejaculates provide functions far beyond the delivery of sperm. Here, we investigated the molecular and functional specificity of the butterfly spermatophore, a structurally complex ejaculate. We discovered that its two distinct parts originate from separate regions of the male reproductive tract, are transferred sequentially during mating, and seem to be the result of a complex evolutionary history. We also highlight a large and previously unrecognized female contribution to the spermatophore, which calls into question traditional characterizations of females as passive recipients of these male ejaculates. 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.