Sequence evolution of the sperm ligand zonadhesin correlates negatively with body weight dimorphism in primates

Gametic selection favours polyandry and selfing

Competition among pollen or sperm (gametic selection) can cause evolution. Mating systems shape the intensity of gametic selection by determining the competitors involved, which can in turn cause the mating system itself to evolve. We model the bidirectional relationship between gametic selection and mating systems, focusing on variation in female mating frequency (monandry-polyandry) and self-fertilisation (selfing-outcrossing). First, we find that monandry and selfing both reduce the efficiency of gametic selection in removing deleterious alleles. This means that selfing can increase mutation load, in contrast to cases without gametic selection where selfing purges deleterious mutations and decreases mutation load. Second, we explore how mating systems evolve via their effect on gametic selection. By manipulating gametic selection, polyandry can evolve to increase the fitness of the offspring produced. However, this indirect advantage of post-copulatory sexual selection is weak and is likely to be overwhelmed by any direct fitness effects of mating systems. Nevertheless, gametic selection can be potentially decisive for selfing evolution because it significantly reduces inbreeding depression, which favours selfing. Thus, the presence of gametic selection could be a key factor driving selfing evolution.

Evolutionary analysis of genes coding for Cysteine-RIch Secretory Proteins (CRISPs) in mammals

BACKGROUND

Cysteine-RIch Secretory Proteins (CRISP) are expressed in the reproductive tract of mammalian males and are involved in fertilization and related processes. Due to their important role in sperm performance and sperm-egg interaction, these genes are likely to be exposed to strong selective pressures, including postcopulatory sexual selection and/or male-female coevolution. We here perform a comparative evolutionary analysis of Crisp genes in mammals. Currently, the nomenclature of CRISP genes is confusing, as a consequence of discrepancies between assignments of orthologs, particularly due to numbering of CRISP genes. This may generate problems when performing comparative evolutionary analyses of mammalian clades and species. To avoid such problems, we first carried out a study of possible orthologous relationships and putative origins of the known CRISP gene sequences. Furthermore, and with the aim to facilitate analyses, we here propose a different nomenclature for CRISP genes (EVAC1-4, "EVolutionarily-analyzed CRISP") to be used in an evolutionary context.

RESULTS

We found differing selective pressures among Crisp genes. CRISP1/4 (EVAC1) and CRISP2 (EVAC2) orthologs are found across eutherian mammals and seem to be conserved in general, but show signs of positive selection in primate CRISP1/4 (EVAC1). Rodent Crisp1 (Evac3a) seems to evolve under a comparatively more relaxed constraint with positive selection on codon sites. Finally, murine Crisp3 (Evac4), which appears to be specific to the genus Mus, shows signs of possible positive selection. We further provide evidence for sexual selection on the sequence of one of these genes (Crisp1/4) that, unlike others, is thought to be exclusively expressed in male reproductive tissues.

CONCLUSIONS

We found differing selective pressures among CRISP genes and sexual selection as a contributing factor in CRISP1/4 gene sequence evolution. Our evolutionary analysis of this unique set of genes contributes to a better understanding of Crisp function in particular and the influence of sexual selection on reproductive mechanisms in general.

Genes Encoding Mammalian Oviductal Proteins Involved in Fertilization are Subjected to Gene Death and Positive Selection

Oviductal proteins play an important role in mammalian fertilization, as proteins from seminal fluid. However, in contrast with the latter, their phylogenetic evolution has been poorly studied. Our objective was to study in 16 mammals the evolution of 16 genes that encode oviductal proteins involved in at least one of the following steps: (1) sperm–oviduct interaction, (2) acrosome reaction, and/or (3) sperm–zona pellucida interaction. Most genes were present in all studied mammals. However, some genes were lost along the evolution of mammals and found as pseudogenes: annexin A5 (ANXA5) and deleted in malignant brain tumor 1 (DMBT1) in tarsier; oviductin (OVGP1) in megabat; and probably progestagen-associated endometrial protein (PAEP) in tarsier, mouse, rat, rabbit, dolphin, and megabat; prostaglandin D2 synthase (PTGDS) in microbat; and plasminogen (PLG) in megabat. Four genes [ANXA1, ANXA4, ANXA5, and heat shock 70 kDa protein 5 (HSPA5)] showed branch-site positive selection, whereas for seven genes [ANXA2, lactotransferrin (LTF), OVGP1, PLG, S100 calcium-binding protein A11 (S100A11), Sperm adhesion molecule 1 (SPAM1), and osteopontin (SPP1)] branch-site model and model-site positive selection were observed. These results strongly suggest that genes encoding oviductal proteins that are known to be important for gamete fertilization are subjected to positive selection during evolution, as numerous genes encoding proteins from mammalian seminal fluid. This suggests that such a rapid evolution may have as a consequence that two isolated populations become separate species more rapidly.

Comparative Sperm Proteomics in Mouse Species with Divergent Mating Systems

Sexual selection is the pervasive force underlying the dramatic divergence of sperm form and function. Although it has been demonstrated that testis gene expression evolves rapidly, exploration of the proteomic basis of sperm diversity is in its infancy. We have employed a whole-cell proteomics approach to characterize sperm divergence among closely related Mus species that experience different sperm competition regimes and exhibit pronounced variation in sperm energetics, motility and fertilization capacity. Interspecific comparisons revealed significant abundance differences amongst proteins involved in fertilization capacity, including those that govern sperm-zona pellucida interactions, axoneme components and metabolic proteins. Ancestral reconstruction of relative testis size suggests that the reduction of zona pellucida binding proteins and heavy-chain dyneins was associated with a relaxation in sperm competition in the M. musculus lineage. Additionally, the decreased reliance on ATP derived from glycolysis in high sperm competition species was reflected in abundance decreases in glycolytic proteins of the principle piece in M. spretus and M. spicilegus. Comparison of protein abundance and stage-specific testis expression revealed a significant correlation during spermatid development when dynamic morphological changes occur. Proteins underlying sperm diversification were also more likely to be subject to translational repression, suggesting that sperm composition is influenced by the evolution of translation control mechanisms. The identification of functionally coherent classes of proteins relating to sperm competition highlights the utility of evolutionary proteomic analyses and reveals that both intensified and relaxed sperm competition can have a pronounced impact on the molecular composition of the male gamete.

Selective constraints on protamine 2 in primates and rodents

BACKGROUND

Protamines are sperm nuclear proteins with a crucial role in chromatin condensation. Their function is strongly linked to sperm head morphology and male fertility. Protamines appear to be affected by a complex pattern of selective constraints. Previous studies showed that sexual selection affects protamine coding sequence and expression in rodents. Here we analyze selective constraints and post-copulatory sexual selection acting on protamine 2 (Prm2) gene sequences of 53 species of primates and rodents. We focused on possible differences in selective constraints between these two clades and on the two functional domains of PRM2 (cleaved- and mature-PRM2). We also assessed if and how changes in Prm2 coding sequence may affect sperm head dimensions.

RESULTS

The domain of Prm2 that is cleaved off during binding to DNA (cleaved-Prm2) was found to be under purifying selection in both clades, whereas the domain that remains bound to DNA (mature-Prm2) was found to be positively selected in primates and under relaxed constraint in rodents. Changes in cleaved-Prm2 coding sequence are significantly correlated to sperm head width and elongation in rodents. Contrary to expectations, a significant effect of sexual selection was not found on either domain or clade.

CONCLUSIONS

Mature-PRM2 may be free to evolve under less constraint due to the existence of PRM1 as a more conserved and functionally redundant copy. The cleaved-PRM2 domain seems to play an important role in sperm head shaping. However, sexual selection on its sequence may be difficult to detect until it is identified which sperm head phenotype (shape and size) confers advantages for sperm performance in different mammalian clades.

Sexual selection and the adaptive evolution of PKDREJ protein in primates and rodents

PKDREJ is a testis-specific protein thought to be located on the sperm surface. Functional studies in the mouse revealed that loss of PKDREJ has effects on sperm transport and the ability to undergo an induced acrosome reaction. Thus, PKDREJ has been considered a potential target of post-copulatory sexual selection in the form of sperm competition. Proteins involved in reproductive processes often show accelerated evolution. In many cases, this rapid divergence is promoted by positive selection which may be driven, at least in part, by post-copulatory sexual selection. We analysed the evolution of the PKDREJ protein in primates and rodents and assessed whether PKDREJ divergence is associated with testes mass relative to body mass, which is a reliable proxy of sperm competition levels. Evidence of an association between the evolutionary rate of the PKDREJ gene and testes mass relative to body mass was not found in primates. Among rodents, evidence of positive selection was detected in the Pkdrej gene in the family Cricetidae but not in Muridae. We then assessed whether Pkdrej divergence is associated with episodes of sperm competition in these families. We detected a positive significant correlation between the evolutionary rates of Pkdrej and testes mass relative to body mass in cricetids. These findings constitute the first evidence of post-copulatory sexual selection influencing the evolution of a protein that participates in the mechanisms regulating sperm transport and the acrosome reaction, strongly suggesting that positive selection may act on these fertilization steps, leading to advantages in situations of sperm competition.

Evolution of protamine genes and changes in sperm head phenotype in rodents

Little is known about the genetic basis of evolutionary changes in sperm phenotype. Postcopulatory sexual selection is associated with differences in protamine gene sequences and promoters and is a powerful force acting on sperm form and function, although links between protamine evolution and sperm phenotype are scarce. Protamines are involved in sperm chromatin condensation, and protamine deficiency negatively affects sperm morphology and male fertility, thus suggesting that they are important for sperm design and function. We examined changes in protamine genes and sperm phenotype in rodents to understand the role of sexual selection on protamine evolution and sperm design. We performed a genotype-phenotype association study using root-to-tip dN/dS (nonsynonymous/synonymous substitutions rate ratio) to account for evolutionary rates and phylogenetic generalized least squares analyses to compare genetic and morphometric data. Evolutionary rates of protamine 1 and the protamine 2 domain cleaved off during chromatin condensation correlated with head size and elongation. Protamine 1 exhibited restricted positive selection on some functional sites, which seemed sufficient to preserve its role in head design. The cleaved-protamine 2, whose relaxation is halted by sexual selection, seems to ensure small, elongated heads that would make sperm more competitive. No association existed between mature-protamine 2 and head phenotype, suggesting little involvement during chromatin condensation and a likely role maintaining the condensed state. Our results suggest that evolutionary changes in protamines could be related to complex developmental modifications in the sperm head. This represents an important step toward understanding the role of changes in gene coding sequences in the divergence of germ cell phenotype.

Mating systems and protein-protein interactions determine evolutionary rates of primate sperm proteins

To assess the relative impact of functional constraint and post-mating sexual selection on sequence evolution of reproductive proteins, we examined 169 primate sperm proteins. In order to recognize potential genome-wide trends, we additionally analysed a sample of altogether 318 non-reproductive (brain and postsynaptic) proteins. Based on cDNAs of eight primate species (Anthropoidea), we observed that pre-mating sperm proteins engaged in sperm composition and assembly show significantly lower incidence of site-specific positive selection and overall lower non-synonymous to synonymous substitution rates (d_N/d_S) across sites as compared with post-mating sperm proteins involved in capacitation, hyperactivation, acrosome reaction and fertilization. Moreover, database screening revealed overall more intracellular protein interaction partners in pre-mating than in post-mating sperm proteins. Finally, post-mating sperm proteins evolved at significantly higher evolutionary rates than pre-mating sperm and non-reproductive proteins on the branches to multi-male breeding species, while no such increase was observed on the branches to unimale and monogamous species. We conclude that less protein–protein interactions of post-mating sperm proteins account for lowered functional constraint, allowing for stronger impact of post-mating sexual selection, while the opposite holds true for pre-mating sperm proteins. This pattern is particularly strong in multi-male breeding species showing high female promiscuity.

Evolutionary Modeling of Genotype-Phenotype Associations, and Application to Primate Coding and Non-coding mtDNA Rate Variation

Variation in substitution rates across a phylogeny can be indicative of shifts in the evolutionary dynamics of a protein or non-protein coding regions. One way to understand these signals is to seek the phenotypic correlates of rate variation. Here, we extended a previously published likelihood method designed to detect evolutionary associations between genotypic evolutionary rate and phenotype over a phylogeny. In simulation with two discrete categories of phenotype, the method has a low false-positive rate and detects greater than 80% of true-positives with a tree length of three or greater and a three-fold or greater change in substitution rate given the phenotype. In addition, we successfully extend the test from two to four phenotype categories and evaluated its performance. We then applied the method to two major hypotheses for rate variation in the mitochondrial genome of primates-longevity and generation time as well as body mass which is correlated with many aspects of life history-using three categories of phenotype through discretization of continuous values. Similar to previous results for mammals, we find that the majority of mitochondrial protein-coding genes show associations consistent with the longevity and body mass predictions and that the predominant signal of association comes from the third codon position. We also found a significant association between maximum lifespan and the evolutionary rate of the control region of the mtDNA. In contrast, 24 protein-coding genes from the nuclear genome do not show a consistent pattern of association, which is inconsistent with the generation time hypothesis. These results show the extended method can robustly identify genotype-phenotype associations up to at least four phenotypic categories, and demonstrate the successful application of the method to study factors affecting neutral evolutionary rate in protein-coding and non-coding loci.

Comparative population genomics of the ejaculate in humans and the great apes

The rapid molecular evolution of reproductive genes is nearly ubiquitous across animals, yet the selective forces and functional targets underlying this divergence remain poorly understood. Humans and closely related species of great apes show strongly divergent mating systems, providing a powerful system to investigate the influence of sperm competition on the evolution of reproductive genes. This is complemented by detailed information on male reproductive biology and unparalleled genomic resources in humans. Here, we have used custom microarrays to capture and sequence 285 genes encoding proteins present in the ejaculate as well as 101 randomly selected control genes in 21 gorillas, 20 chimpanzees, 20 bonobos, and 20 humans. In total, we have generated >25× average genomic coverage per individual for over 1 million target base pairs. Our analyses indicate high levels of evolutionary constraint across much of the ejaculate combined with more rapid evolution of genes involved in immune defense and proteolysis. We do not find evidence for appreciably more positive selection along the lineage leading to bonobos and chimpanzees, although this would be predicted given more intense sperm competition in these species. Rather, the extent of positive and negative selection depended more on the effective population sizes of the species. Thus, general patterns of male reproductive protein evolution among apes and humans depend strongly on gene function but not on inferred differences in the intensity of sperm competition among extant species.

Decoupling of rapid and adaptive evolution among seminal fluid proteins in heliconius butterflies with divergent mating systems

Reproductive proteins often diverge rapidly between species. This pattern is frequently attributed to postmating sexual selection. Heliconius butterflies offer a good opportunity to examine this hypothesis by contrasting patterns of reproductive protein evolution between clades with divergent mating systems. Pupal-mating Heliconius females typically mate only once, limiting opportunity for postmating sexual selection. In contrast, adult-mating females remate throughout life. Reproductive protein evolution is therefore predicted to be slower and show little evidence of positive selection in the pupal-mating clade. We examined this prediction by sequencing 18 seminal fluid protein genes from a dozen Heliconius species and a related outgroup. Two proteins exhibited dN/dS > 1, implicating positive selection in the rapid evolution of at least a few Heliconius seminal fluid proteins. However, contrary to predictions, the average evolutionary rate of seminal fluid proteins was greater among pupal-mating Heliconius. Based on these results, we suggest that positive selection and relaxed constraint can generate conflicting patterns of reproductive protein evolution between mating systems. As predicted, some loci may show elevated evolutionary rates in promiscuous taxa relative to monandrous taxa resulting from adaptations to postmating sexual selection. However, when monandry is derived (as in Heliconius), the opposite pattern may result from relaxed selective constraints.

The molecular evolution of animal reproductive tract proteins: what have we learned from mating-system comparisons?

Postcopulatory sexual selection is thought to drive the rapid evolution of reproductive tract genes in many animals. Recently, a number of studies have sought to test this hypothesis by examining the effects of mating system variation on the evolutionary rates of reproductive tract genes. Perhaps surprisingly, there is relatively little evidence that reproductive proteins evolve more rapidly in species subject to strong postcopulatory sexual selection. This emerging trend may suggest that other processes, such as host-pathogen interactions, are the main engines of rapid reproductive gene evolution. I suggest that such a conclusion is as yet unwarranted; instead, I propose that more rigorous analytical techniques, as well as multigene and population-based approaches, are required for a full understanding of the consequences of mating system variation for the evolution of reproductive tract genes.

Testing the effects of mating system variation on rates of molecular evolution in primates

Post-copulatory sexual selection has been proposed to drive the rapid evolution of reproductive proteins, and, more recently,to increase genome-wide mutation rates. Comparisons of rates of molecular evolution between lineages with different levels of female multiple mating represent a promising, but under-utilized, approach for testing the effects of sperm competition on sequence evolution. Here, I use comparisons between primate species with divergent mating systems to examine the effects of sperm competition on reproductive protein evolution, as well as on sex-averaged mutation rates. Rates of nonsynonymous substitution are higher for testis-specific genes along the chimpanzee lineage in comparison to the human lineage, consistent with expectations. However, the data reported here do not allow firm conclusions concerning the effects of mating system on genome-wide mutation rates, with different results obtained from different species pairs. Ultimately, comparative studies encompassing a range of mating systems and other life history traits will be required to make broad generalizations concerning the genomic effects of sperm competition.

Genotype-phenotype associations: substitution models to detect evolutionary associations between phenotypic variables and genotypic evolutionary rate

Motivation: Mapping between genotype and phenotype is one of the primary goals of evolutionary genetics but one that has received little attention at the interspecies level. Recent developments in phylogenetics and statistical modelling have typically been used to examine molecular and phenotypic evolution separately. We have used this background to develop phylogenetic substitution models to test for associations between evolutionary rate of genotype and phenotype. We do this by creating hybrid rate matrices between genotype and phenotype.

Results: Simulation results show our models to be accurate in detecting genotype–phenotype associations and robust for various factors that typically affect maximum likelihood methods, such as number of taxa, level of relevant signal, proportion of sites affected and length of evolutionary divergence. Further, simulations show that our method is robust to homogeneity assumptions. We apply the models to datasets of male reproductive system genes in relation to mating systems of primates. We show that evolution of semenogelin II is significantly associated with mating systems whereas two negative control genes (cytochrome b and peptidase inhibitor 3) show no significant association. This provides the first hybrid substitution model of which we are aware to directly test the association between genotype and phenotype using a phylogenetic framework.

Availability: Perl and HYPHY scripts are available upon request from the authors.

Contact:to252@cam.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Sexual selection drives weak positive selection in protamine genes and high promoter divergence, enhancing sperm competitiveness

Phenotypic adaptations may be the result of changes in gene structure or gene regulation, but little is known about the evolution of gene expression. In addition, it is unclear whether the same selective forces may operate at both levels simultaneously. Reproductive proteins evolve rapidly, but the underlying selective forces promoting such rapid changes are still a matter of debate. In particular, the role of sexual selection in driving positive selection among reproductive proteins remains controversial, whereas its potential influence on changes in promoter regions has not been explored. Protamines are responsible for maintaining DNA in a compacted form in chromosomes in sperm and the available evidence suggests that they evolve rapidly. Because protamines condense DNA within the sperm nucleus, they influence sperm head shape. Here, we examine the influence of sperm competition upon protamine 1 and protamine 2 genes and their promoters, by comparing closely related species of Mus that differ in relative testes size, a reliable indicator of levels of sperm competition. We find evidence of positive selection in the protamine 2 gene in the species with the highest inferred levels of sperm competition. In addition, sperm competition levels across all species are strongly associated with high divergence in protamine 2 promoters that, in turn, are associated with sperm swimming speed. We suggest that changes in protamine 2 promoters are likely to enhance sperm swimming speed by making sperm heads more hydrodynamic. Such phenotypic changes are adaptive because sperm swimming speed may be a major determinant of fertilization success under sperm competition. Thus, when species have diverged recently, few changes in gene-coding sequences are found, while high divergence in promoters seems to be associated with the intensity of sexual selection.

EST analysis of male accessory glands from Heliconius butterflies with divergent mating systems

Background

Heliconius butterflies possess a remarkable diversity of phenotypes, physiologies, and behaviors that has long distinguished this genus as a focal taxon in ecological and evolutionary research. Recently Heliconius has also emerged as a model system for using genomic methods to investigate the causes and consequences of biological diversity. One notable aspect of Heliconius diversity is a dichotomy in mating systems which provides an unusual opportunity to investigate the relationship between sexual selection and the evolution of reproductive proteins. As a first step in pursuing this research, we report the generation and analysis of expressed sequence tags (ESTs) from the male accessory gland of H. erato and H. melpomene, species representative of the two mating systems present in the genus Heliconius.

Results

We successfully sequenced 933 ESTs clustering into 371 unigenes from H. erato and 1033 ESTs clustering into 340 unigenes from H. melpomene. Results from the two species were very similar. Approximately a third of the unigenes showed no significant BLAST similarity (E-value <10^-5) to sequences in GenBank's non-redundant databases, indicating that a large proportion of novel genes are expressed in Heliconius male accessory glands. In both species only a third of accessory gland unigenes were also found among genes expressed in wing tissue. About 25% of unigenes from both species encoded secreted proteins. This includes three groups of highly abundant unigenes encoding repetitive proteins considered to be candidate seminal fluid proteins; proteins encoded by one of these groups were detected in H. erato spermatophores.

Conclusion

This collection of ESTs will serve as the foundation for the future identification and evolutionary analysis of male reproductive proteins in Heliconius butterflies. These data also represent a significant advance in the rapidly growing collection of genomic resources available in Heliconius butterflies. As such, they substantially enhance this taxon as a model system for investigating questions of ecological, phenotypic, and genomic diversity.

Adaptive evolution in an avian reproductive protein: ZP3

Proteins involved in reproduction appear to be evolving adaptively across taxa. This rapid evolution is thought to be the result of forces involved in sexual selection. One of the most often suggested of these forces is sexual conflict involving sperm competition and polyspermy avoidance. Bird species offer a unique opportunity to test this hypothesis since the avian egg coat tolerates physiological polyspermy, or the penetration of multiple sperm during fertilization, without negative effects on later development. Despite this, and the extensive amount of data gathered on sexual selection in birds, there are limited studies on the patterns of evolution of avian reproductive proteins. Here we present an analysis of the pattern of evolution of Zona Pellucida 3 (ZP3), a protein present on the avian egg coat. We found that, across several galliform and a single anseriform species, ZP3 appears to be diverging by positive adaptive evolution. In an exploratory analysis of portions of the gene in Callipepla californica we also found evidence of a selective sweep at the putative sperm binding region of the protein. In sum, ZP3 in birds, like reproductive proteins in other species, appears to be adaptively evolving. This result suggests that polyspermy avoidance is not sufficient to explain positive Darwinian selection in reproductive proteins across taxonomic groups. Clearly, the inclusion of bird species in the study of reproductive proteins across taxa promises to add greatly to the discussion of the factors driving the widespread phenomenon of adaptive evolution in reproductive proteins.