Domain Expansion and Functional Diversification in Vertebrate Reproductive Proteins.

The rapid evolution of fertilization proteins has generated remarkable diversity in molecular structure and function. Glycoproteins of vertebrate egg coats contain multiple zona pellucida (ZP)-N domains (1-6 copies) that facilitate multiple reproductive functions, including species-specific sperm recognition. In this report, we integrate phylogenetics and machine learning to investigate how ZP-N domains diversify in structure and function. The most C-terminal ZP-N domain of each paralog is associated with another domain type (ZP-C), which together form a "ZP module." All modular ZP-N domains are phylogenetically distinct from nonmodular or free ZP-N domains. Machine learning-based classification identifies eight residues that form a stabilizing network in modular ZP-N domains that is absent in free domains. Positive selection is identified in some free ZP-N domains. Our findings support that strong purifying selection has conserved an essential structural core in modular ZP-N domains, with the relaxation of this structural constraint allowing free N-terminal domains to functionally diversify.

Molecular mechanisms and evolution of fertilization proteins.

Sexual reproduction involves a cascade of molecular interactions between the sperm and the egg culminating in cell-cell fusion. Vital steps mediating fertilization include chemoattraction of the sperm to the egg, induction of the sperm acrosome reaction, dissolution of the egg coat, and sperm-egg plasma membrane binding and fusion. Despite decades of research, only a handful of interacting gamete recognition proteins (GRPs) have been identified across taxa mediating each of these steps, most notably in abalone, sea urchins, and mammals. This review outlines and compares notable GRP pairs mediating sperm-egg recognition in these three significant model systems and discusses the molecular basis of species-specific fertilization driven by GRP function. In addition, we explore the evolutionary theory behind the rapid diversification of GRPs between species. In particular, we focus on how the coevolution between interacting sperm and egg proteins may contribute to the formation of boundaries to hybridization. Finally, we discuss how pairing structural information with evolutionary insights can improve our understanding of mechanisms of fertilization and their origins.

Proteomics support the threespine stickleback egg coat as a protective oocyte envelope.

Ancestrally marine threespine stickleback fish (Gasterosteus aculeatus) have undergone an adaptive radiation into freshwater environments throughout the Northern Hemisphere, creating an excellent model system for studying molecular adaptation and speciation. Ecological and behavioral factors have been suggested to underlie stickleback reproductive isolation and incipient speciation, but reproductive proteins mediating gamete recognition during fertilization have so far remained unexplored. To begin to investigate the contribution of reproductive proteins to stickleback reproductive isolation, we have characterized the stickleback egg coat proteome. We find that stickleback egg coats are comprised of homologs to the zona pellucida (ZP) proteins ZP1 and ZP3, as in other teleost fish. Our molecular evolutionary analyses indicate that across teleosts, ZP3 but not ZP1 has experienced positive Darwinian selection. Mammalian ZP3 is also rapidly evolving, and surprisingly some residues under selection in stickleback and mammalian ZP3 directly align. Despite broad homology, however, we find differences between mammalian and stickleback ZP proteins with respect to glycosylation, disulfide bonding, and sites of synthesis. Taken together, the changes we observe in stickleback ZP protein architecture suggest that the egg coats of stickleback fish, and perhaps fish more generally, have evolved to fulfill a more protective functional role than their mammalian counterparts.

Solution structure of sperm lysin yields novel insights into molecular dynamics of rapid protein evolution.

Protein evolution is driven by the sum of different physiochemical and genetic processes that usually results in strong purifying selection to maintain biochemical functions. However, proteins that are part of systems under arms race dynamics often evolve at unparalleled rates that can produce atypical biochemical properties. In the marine mollusk abalone, lysin and vitelline envelope receptor for lysin (VERL) are a pair of rapidly coevolving proteins that are essential for species-specific interactions between sperm and egg. Despite extensive biochemical characterization of lysin-including crystal structures of multiple orthologs-it was unclear how sites under positive selection may facilitate recognition of VERL. Using a combination of targeted mutagenesis and multidimensional NMR, we present a high-definition solution structure of sperm lysin from red abalone (Haliotis rufescens). Unapparent from the crystallography data, multiple NMR-based analyses conducted in solution reveal clustering of the N and C termini to form a nexus of 13 positively selected sites that constitute a VERL binding interface. Evolutionary rate was found to be a significant predictor of backbone flexibility, which may be critical for lysin bioactivity and/or accelerated evolution. Flexible, rapidly evolving segments that constitute the VERL binding interface were also the most distorted regions of the crystal structure relative to what was observed in solution. While lysin has been the subject of extensive biochemical and evolutionary analyses for more than 30 years, this study highlights the enhanced insights gained from applying NMR approaches to rapidly evolving proteins.

A conversation with Mariana Wolfner, newly inducted member of the National Academy of Sciences.

Molecular Reproduction and Development is delighted to announce that editorial board member Mariana F. Wolfner has been elected to the National Academy of Sciences. Here, Dr Wolfner is interviewed by two of her former postdocs. She discusses her path to studying reproduction and her career as a researcher and mentor.

A Diverse Lipid Antigen-Specific TCR Repertoire Is Clonally Expanded during Active Tuberculosis.

Human T cells that recognize lipid Ags presented by highly conserved CD1 proteins often express semi-invariant TCRs, but the true diversity of lipid Ag-specific TCRs remains unknown. We use CD1b tetramers and high-throughput immunosequencing to analyze thousands of TCRs from ex vivo-sorted or in vitro-expanded T cells specific for the mycobacterial lipid Ag, glucose monomycolate. Our results reveal a surprisingly diverse repertoire resulting from editing of germline-encoded gene rearrangements analogous to MHC-restricted TCRs. We used a distance-based metric (TCRDist) to show how this diverse TCR repertoire builds upon previously reported conserved motifs by including subject-specific TCRs. In a South African cohort, we show that TCRDist can identify clonal expansion of diverse glucose monomycolate-specific TCRs and accurately distinguish patients with active tuberculosis from control subjects. These data suggest that similar mechanisms govern the selection and expansion of peptide and lipid Ag-specific T cells despite the nonpolymorphic nature of CD1.

Conservation of molecular and cellular phenotypes of invariant NKT cells between humans and non-human primates.

Invariant NKT (iNKT) cells in both humans and non-human primates are activated by the glycolipid antigen, α-galactosylceramide (α-GalCer). However, the extent to which the molecular mechanisms of antigen recognition and in vivo phenotypes of iNKT cells are conserved among primate species has not been determined. Using an evolutionary genetic approach, we found a lack of diversifying selection in CD1 genes over 45 million years of evolution, which stands in stark contrast to the history of the MHC system for presenting peptide antigens to T cells. The invariant T cell receptor (TCR)-α chain was strictly conserved across all seven primate clades. Invariant NKT cells from rhesus macaques (Macaca mulatta) bind human CD1D-α-GalCer tetramer and are activated by α-GalCer-loaded human CD1D transfectants. The dominant TCR-ß chain cloned from a rhesus-derived iNKT cell line is nearly identical to that found in the human iNKT TCR, and transduction of the rhesus iNKT TCR into human Jurkat cells show that it is sufficient for binding human CD1D-α-GalCer tetramer. Finally, we used a 20-color flow cytometry panel to probe tissue phenotypes of iNKT cells in a cohort of rhesus macaques. We discovered several tissue-resident iNKT populations that have not been previously described in non-human primates but are known in humans, such as TCR-γδ iNKTs. These data reveal a diversity of iNKT cell phenotypes despite convergent evolution of the genes required for lipid antigen presentation and recognition in humans and non-human primates.

Quantitative evolutionary proteomics of seminal fluid from primates with different mating systems.

BACKGROUND: Genomic data from various organisms have been used to study how sexual selection has shaped genetic diversity in reproductive proteins, and in particular, to elucidate how mating systems may have influenced evolution at the molecular and phenotypic levels. However, large-scale proteomic data including protein identifications and abundances are only now entering the field of evolutionary and comparative genomics. Variation in both protein sequence and expression level may play important roles in the evolution of sexual traits and behaviors. RESULTS: Here, we broadly analyze the components of seminal fluid from primates with diverse mating systems ranging from monogamous to polygynous, and include genomics, proteomics, phylogenetic and quantitative characters into our framework. Our analyses show that seminal fluid proteins are undergoing rapid evolution and some of these quickly evolving proteins may be influenced by sexual selection. Through evolutionary analyses and protein abundance differences, we identified 84 genes whose evolutionary rates or expression levels were correlated with mating system and other sexual characters. We found that many proteins differ in abundance between monogamous and polygynous primate mating systems. Many of these proteins are enriched in the copulatory plug pathway, which suggests that post-zygotic selective barriers are important regardless of mating system type. CONCLUSIONS: This work is the first to comprehensively compare seminal fluid proteins between human and non-human primates using high-throughput proteomics. Our findings highlight the impact of mating system variation on seminal fluid protein evolution and abundance.

Duplicate abalone egg coat proteins bind sperm lysin similarly, but evolve oppositely, consistent with molecular mimicry at fertilization.

Sperm and egg proteins constitute a remarkable paradigm in evolutionary biology: despite their fundamental role in mediating fertilization (suggesting stasis), some of these molecules are among the most rapidly evolving ones known, and their divergence can lead to reproductive isolation. Because of strong selection to maintain function among interbreeding individuals, interacting fertilization proteins should also exhibit a strong signal of correlated divergence among closely related species. We use evidence of such molecular co-evolution to target biochemical studies of fertilization in North Pacific abalone (Haliotis spp.), a model system of reproductive protein evolution. We test the evolutionary rates (d(N)/d(S)) of abalone sperm lysin and two duplicated egg coat proteins (VERL and VEZP14), and find a signal of co-evolution specific to ZP-N, a putative sperm binding motif previously identified by homology modeling. Positively selected residues in VERL and VEZP14 occur on the same face of the structural model, suggesting a common mode of interaction with sperm lysin. We test this computational prediction biochemically, confirming that the ZP-N motif is sufficient to bind lysin and that the affinities of VERL and VEZP14 are comparable. However, we also find that on phylogenetic lineages where lysin and VERL evolve rapidly, VEZP14 evolves slowly, and vice versa. We describe a model of sexual conflict that can recreate this pattern of anti-correlated evolution by assuming that VEZP14 acts as a VERL mimic, reducing the intensity of sexual conflict and slowing the co-evolution of lysin and VERL.

Selection on plant male function genes identifies candidates for reproductive isolation of yellow monkeyflowers.

Understanding the genetic basis of reproductive isolation promises insight into speciation and the origins of biological diversity. While progress has been made in identifying genes underlying barriers to reproduction that function after fertilization (post-zygotic isolation), we know much less about earlier acting pre-zygotic barriers. Of particular interest are barriers involved in mating and fertilization that can evolve extremely rapidly under sexual selection, suggesting they may play a prominent role in the initial stages of reproductive isolation. A significant challenge to the field of speciation genetics is developing new approaches for identification of candidate genes underlying these barriers, particularly among non-traditional model systems. We employ powerful proteomic and genomic strategies to study the genetic basis of conspecific pollen precedence, an important component of pre-zygotic reproductive isolation among yellow monkeyflowers (Mimulus spp.) resulting from male pollen competition. We use isotopic labeling in combination with shotgun proteomics to identify more than 2,000 male function (pollen tube) proteins within maternal reproductive structures (styles) of M. guttatus flowers where pollen competition occurs. We then sequence array-captured pollen tube exomes from a large outcrossing population of M. guttatus, and identify those genes with evidence of selective sweeps or balancing selection consistent with their role in pollen competition. We also test for evidence of positive selection on these genes more broadly across yellow monkeyflowers, because a signal of adaptive divergence is a common feature of genes causing reproductive isolation. Together the molecular evolution studies identify 159 pollen tube proteins that are candidate genes for conspecific pollen precedence. Our work demonstrates how powerful proteomic and genomic tools can be readily adapted to non-traditional model systems, allowing for genome-wide screens towards the goal of identifying the molecular basis of genetically complex traits.

Evolution of the egg: new findings and challenges.

The evolution of the egg is dynamic, and eggs have numerous species-specific properties across vertebrates and invertebrates. Interestingly, although the structure and function of the egg have remained relatively conserved over time, some constituents of the egg's extracellular barriers are undergoing rapid evolution. In this article, we review current ideas regarding sperm-egg interactions, discuss genetic approaches used to elucidate egg gene functions, and highlight the interesting differences that have evolved across taxa. We suggest that the rapid evolution of egg components and the mechanisms behind sperm-egg interactions are integrally connected, and delve in depth into each component of the egg's extracellular matrices. Finally, we discuss the promising future of reproductive research and how high-throughput genomics and proteomics have the potential to revolutionize the field and provide new evidence that will challenge previously held views about the fertilization process.

Positive selection within a diatom species acts on putative protein interactions and transcriptional regulation.

Diatoms are the most species-rich group of microalgae, and their contribution to marine primary production is important on a global scale. Diatoms can form dense blooms through rapid asexual reproduction; mutations acquired and propagated during blooms likely provide the genetic, and thus phenotypic, variability upon which natural selection may act. Positive selection was tested using genome and transcriptome-wide pair-wise comparisons of homologs in three genera of diatoms (Pseudo-nitzschia, Ditylum, and Thalassiosira) that represent decreasing phylogenetic distances. The signal of positive selection was greatest between two strains of Thalassiosira pseudonana. Further testing among seven strains of T. pseudonana yielded 809 candidate genes of positive selection, which are 7% of the protein-coding genes. Orphan genes and genes encoding protein-binding domains and transcriptional regulators were enriched within the set of positively selected genes relative to the genome as a whole. Positively selected genes were linked to the potential selective pressures of nutrient limitation and sea surface temperature based on analysis of gene expression profiles and identification of positively selected genes in subsets of strains from locations with similar environmental conditions. The identification of positively selected genes presents an opportunity to test new hypotheses in natural populations and the laboratory that integrate selected genotypes in T. pseudonana with their associated phenotypes and selective forces.

Reproduction and immunity-driven natural selection in the human WFDC locus.

The whey acidic protein (WAP) four-disulfide core domain (WFDC) locus located on human chromosome 20q13 spans 19 genes with WAP and/or Kunitz domains. These genes participate in antimicrobial, immune, and tissue homoeostasis activities. Neighboring SEMG genes encode seminal proteins Semenogelin 1 and 2 (SEMG1 and SEMG2). WFDC and SEMG genes have a strikingly high rate of amino acid replacement (dN/dS), indicative of responses to adaptive pressures during vertebrate evolution. To better understand the selection pressures acting on WFDC genes in human populations, we resequenced 18 genes and 54 noncoding segments in 71 European (CEU), African (YRI), and Asian (CHB + JPT) individuals. Overall, we identified 484 single-nucleotide polymorphisms (SNPs), including 65 coding variants (of which 49 are nonsynonymous differences). Using classic neutrality tests, we confirmed the signature of short-term balancing selection on WFDC8 in Europeans and a signature of positive selection spanning genes PI3, SEMG1, SEMG2, and SLPI. Associated with the latter signal, we identified an unusually homogeneous-derived 100-kb haplotype with a frequency of 88% in Asian populations. A putative candidate variant targeted by selection is Thr56Ser in SEMG1, which may alter the proteolytic profile of SEMG1 and antimicrobial activities of semen. All the well-characterized genes residing in the WDFC locus encode proteins that appear to have a role in immunity and/or fertility, two processes that are often associated with adaptive evolution. This study provides further evidence that the WFDC and SEMG loci have been under strong adaptive pressure within the short timescale of modern humans.

Trans genomic capture and sequencing of primate exomes reveals new targets of positive selection.

Comparison of protein-coding DNA sequences from diverse primates can provide insight into these species' evolutionary history and uncover the molecular basis for their phenotypic differences. Currently, the number of available primate reference genomes limits these genome-wide comparisons. Here we use targeted capture methods designed for human to sequence the protein-coding regions, or exomes, of four non-human primate species (three Old World monkeys and one New World monkey). Despite average sequence divergence of up to 4% from the human sequence probes, we are able to capture ~96% of coding sequences. Using a combination of mapping and assembly techniques, we generated high-quality full-length coding sequences for each species. Both the number of nucleotide differences and the distribution of insertion and deletion (indel) lengths indicate that the quality of the assembled sequences is very high and exceeds that of most reference genomes. Using this expanded set of primate coding sequences, we performed a genome-wide scan for genes experiencing positive selection and identified a novel class of adaptively evolving genes involved in the conversion of epithelial cells in skin, hair, and nails to keratin. Interestingly, the genes we identify under positive selection also exhibit significantly increased allele frequency differences among human populations, suggesting that they play a role in both recent and long-term adaptation. We also identify several genes that have been lost on specific primate lineages, which illustrate the broad utility of this data set for other evolutionary analyses. These results demonstrate the power of second-generation sequencing in comparative genomics and greatly expand the repertoire of available primate coding sequences.

Detecting coevolution in mammalian sperm-egg fusion proteins.

Interactions between sperm and egg proteins can occur physically between gamete surface-binding proteins, and genetically between gamete proteins that work in complementary pathways in which they may not physically interact. Physically interacting sperm-egg proteins have been functionally identified in only a few species, and none have been verified within mammals. Candidate genes on both the sperm and egg surfaces exist, but gene deletion studies do not support functional interactions between these sperm-egg proteins; interacting sperm-egg proteins thus remain elusive. Cooperative gamete proteins undergo rapid evolution, and it is predicted that these sperm-egg proteins will also have correlated evolutionary rates due to compensatory changes on both the sperm and egg. To explore potential physical and genetic interactions in sperm-egg proteins, we sequenced four candidate genes from diverse primate species, and used regression and likelihood methods to test for signatures of coevolution between sperm-egg gene pairs. With both methods, we found that the egg protein CD9 coevolves with the sperm protein IZUMO1, suggesting a physical or genetic interaction occurs between them. With regression analysis, we found that CD9 and CRISP2 have correlated rates of evolution, and with likelihood analysis, that CD9 and CRISP1 have correlated rates. This suggests that the different tests may reflect different levels of interaction, be it physical or genetic. Coevolution tests thus provide an exploratory method for detecting potentially interacting sperm-egg protein pairs.

Detecting coevolution through allelic association between physically unlinked loci.

Coevolving interacting genes undergo complementary mutations to maintain their interaction. Distinct combinations of alleles in coevolving genes interact differently, conferring varying degrees of fitness. If this fitness differential is adequately large, the resulting selection for allele matching could maintain allelic association, even between physically unlinked loci. Allelic association is often observed in a population with the use of gametic linkage disequilibrium. However, because the coevolving genes are not necessarily in physical linkage, this is not an appropriate measure of coevolution-induced allelic association. Instead, we propose using both composite linkage disequilibrium (CLD) and a measure of association between genotypes, which we call genotype association (GA). Using a simple selective model, we simulated loci and calculated power for tests of CLD and GA, showing that the tests can detect the allelic association expected under realistic selective pressure. We apply CLD and GA tests to the polymorphic, physically unlinked, and putatively coevolving human gamete-recognition genes ZP3 and ZP3R. We observe unusual allelic association, not attributable to population structure, between ZP3 and ZP3R. This study shows that selection for allele matching can drive allelic association between unlinked loci in a contemporary human population, and that selection can be detected with the use of CLD and GA tests. The observation of this selection is surprising, but reasonable in the highly selected system of fertilization. If confirmed, this sort of selection provides an exception to the paradigm of chromosomal independent assortment.

Mass spectrometry and next-generation sequencing reveal an abundant and rapidly evolving abalone sperm protein.

Abalone, a broadcast spawning marine mollusk, is an important model for molecular interactions and positive selection in fertilization, but the focus has previously been on only two sperm proteins, lysin and sp18. We used genomic and proteomic techniques to bring new insights to this model by characterizing the testis transcriptome and sperm proteome of the Red abalone Haliotis rufescens. One pair of homologous, testis-specific proteins contains a secretion signal and is small, abundant, and associated with the acrosome. Comparative analysis revealed that homologs are extremely divergent between species, and show strong evidence for positive selection. The acrosomal localization and rapid evolution of these proteins indicates that they play an important role in fertilization, and could be involved in the species-specificity of sperm-egg interactions in abalone. Our genomic and proteomic characterization of abalone fertilization resulted in the identification of interesting, novel peptides that have eluded detection in this important model system for 20 years.

The molecular basis of sex: linking yeast to human.

Species-specific recognition between egg and sperm, a crucial event that marks the beginning of fertilization in multicellular organisms, mirrors the binding between haploid cells of opposite mating type in unicellular eukaryotes such as yeast. However, as implied by the lack of sequence similarity between sperm-binding regions of invertebrate and vertebrate egg coat proteins, these interactions are thought to rely on completely different molecular entities. Here, we argue that these recognition systems are, in fact, related: despite being separated by 0.6-1 billion years of evolution, functionally essential domains of a mollusc sperm receptor and a yeast mating protein adopt the same 3D fold as egg zona pellucida proteins mediating the binding between gametes in humans.

Proteomics enhances evolutionary and functional analysis of reproductive proteins.

Reproductive proteins maintain species-specific barriers to fertilization, affect the outcome of sperm competition, mediate reproductive conflicts between the sexes, and potentially contribute to the formation of new species. However, the specific proteins and molecular mechanisms that underlie these processes are understood in only a handful of cases. Advances in genomic and proteomic technologies enable the identification of large suites of reproductive proteins, making it possible to dissect reproductive phenotypes at the molecular level. We first review these technological advances and describe how reproductive proteins are identified in diverse animal taxa. We then discuss the dynamic evolution of reproductive proteins and the potential selective forces that act on them. Finally, we describe molecular and genomic tools for functional analysis and detail how evolutionary data may be used to make predictions about interactions among reproductive proteins.