Positive selection within sperm-egg adhesion domains of fertilin: an ADAM gene with a potential role in fertilization

Sperm long non-coding RNAs as markers for ram fertility

It is critical in sheep farming to accurately estimate ram fertility for maintaining reproductive effectiveness and for production profitability. However, there is currently a lack of reliable biomarkers to estimate semen quality and ram fertility, which is hindering advances in animal science and technology. The objective of this study was to uncover long non-coding RNAs (lncRNAs) in sperm from rams with distinct fertility phenotypes. Mature rams were allocated into two groups: high and low fertility (HF; n = 31; 94.5 ± 2.8%, LF; n = 25; 83.1 ± 5.73%; P = 0.028) according to the pregnancy rates sired by the rams (average pregnancy rate; 89.4 ± 7.2%). Total RNAs were isolated from sperm of the highest- and lowest-fertility rams (n = 4, pregnancy rate; 99.2 ± 1.6%, and 73.6 ± 4.4%, respectively) followed by next-generation sequencing of the transcripts. We uncovered 11,209 lncRNAs from the sperm of rams with HF and LF. In comparison to each other, there were 93 differentially expressed (DE) lncRNAs in sperm from the two distinct fertility phenotypes. Of these, 141 mRNAs were upregulated and 134 were downregulated between HF and LF, respectively. Genes commonly enriched for 9 + 2 motile cilium and sperm flagellum were ABHD2, AK1, CABS1, ROPN1, SEPTIN2, SLIRP, and TEKT3. Moreover, CABS1, CCDC39, CFAP97D1, ROPN1, SLIRP, TEKT3, and TTC12 were commonly enriched in flagellated sperm motility and sperm motility. Differentially expressed mRNAs were enriched in the top 16 KEGG pathways. Targets of the differentially expressed lncRNAs elucidate functions in cis and trans manner using the genetic context of the lncRNA locus, and lncRNA sequences revealed 471 mRNAs targets of 10 lncRNAs. This study illustrates the existence of potential lncRNA biomarkers that can be implemented in analyzing the quality of ram sperm and determining the sperm fertility and is used in breeding soundness exams for precision livestock farming to ensure food security on a global scale.

De novo protein identification in mammalian sperm using in situ cryoelectron tomography and AlphaFold2 docking

To understand the molecular mechanisms of cellular pathways, contemporary workflows typically require multiple techniques to identify proteins, track their localization, and determine their structures in vitro. Here, we combined cellular cryoelectron tomography (cryo-ET) and AlphaFold2 modeling to address these questions and understand how mammalian sperm are built in situ. Our cellular cryo-ET and subtomogram averaging provided 6.0-Å reconstructions of axonemal microtubule structures. The well-resolved tertiary structures allowed us to unbiasedly match sperm-specific densities with 21,615 AlphaFold2-predicted protein models of the mouse proteome. We identified Tektin 5, CCDC105, and SPACA9 as novel microtubule-associated proteins. These proteins form an extensive interaction network crosslinking the lumen of axonemal doublet microtubules, suggesting their roles in modulating the mechanical properties of the filaments. Indeed, Tekt5 -/- sperm possess more deformed flagella with 180° bends. Together, our studies presented a cellular visual proteomics workflow and shed light on the in vivo functions of Tektin 5.

RNA-sequencing attest increased sperm motility in bovine spermatozoa treated with ethanolic extract of Putranjiva roxburghii

In the course of time, scientific communities have a growing interest in understanding ethano medicines. The Putranjiva roxburghii, a native plant of the Indian Subcontinent is described as a "Child amulet tree" in Ayurveda. Based on the fact that this herbal medicine has an indispensable component of integrative medicine, the present study was planned to assess the effect of ethanolic dried extract of Putranjiva seeds on the motility of X and Y-bearing bovine spermatozoa. The in-vitro effect of seed extract diluted in S-TALP medium on bull semen has been evaluated by Computer Assisted Semen Analysis (CASA) shows a marked increase in the motility of spermatozoa. Motile and non-motile spermatozoa have been separated by glass wool column from the control as well as treated group. The X and Y-bearing sperm quantification have been carried out by droplet digital polymerase chain reaction (ddPCR). The extract didn't exert any differential effect on the motility and viability of X and Y chromosome-bearing spermatozoa. The transcriptome profiling (RNA-Seq) identified 93 differentially expressed genes between the extract treated and control group. It unveils the up-regulation of CATSPER, AKAP3, SPAG, ADAM1B, ADAM2 and ADAM32 genes that are involved in increasing sperm motility. Transcriptome profile also unveil the expression of ZAR1, CYP17A1, APPL2, HOXB4 and SP9 genes involved with embryonic development processes in Putranjiva extract-treated motile spermatozoa. The results envisaged the medicinal value of Putranjiva herb on increased fertility due to combinatory effect like increased sperm motility and favourableness on embryogenesis. The study ruled out the possibility of herbs having any biased effect on the selection of either male or female-bearing spermatozoa in the bull.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03452-4.

Evolution of Fertilization-Related Genes Provides Insights Into Reproductive Health in Natural Ascrotal Mammals

Cryptorchidism is the failure of one or both testes to descend into the bottom of the scrotum. This common congenital malformation in humans and domestic animals is the best characterized risk factor for abnormal sperm functions and infertility. However, current treatment approaches for cryptorchidism do not ensure paternity in all cases. Some lineages of mammals (such as elephants and cetaceans) have natural ascrotal testes (i.e., undescended or incompletely descended testes) and normal sperm motility and fertility, providing an opportunity to understand the genetic basis of cryptorchidism. In this study, we showed that genes associated with sperm motility and competition/fertility in ascrotal mammals experienced frequent, strong selective pressure. The fixation of specific amino acids and positive selection in ascrotal mammals could affect the physicochemical properties and functions of fertilization-related proteins. In a comparison between mammals with undescended testes and incompletely descended testes, discrepancies in genes showing evidence for adaptive evolution and in functional enrichment suggested that multiple molecular mechanisms contribute to the maintenance of fertility in the challenging testicular environment. Our findings revealed substantial heterogeneity in the divergence of fertilization-related genes between natural scrotal and ascrotal mammals and provide insight into molecular mechanisms underlying normal sperm motility and competition in natural ascrotal mammals. We provide a detailed theoretical basis for understanding the pathology of cryptorchidism from a molecular evolutionary perspective. This study may contribute to the establishment of diagnostic and therapeutic targets for sperm motility and fertility disorders due to congenital cryptorchidism in humans and domestic animals.

The Importance of Gene Duplication and Domain Repeat Expansion for the Function and Evolution of Fertilization Proteins

The process of gene duplication followed by gene loss or evolution of new functions has been studied extensively, yet the role gene duplication plays in the function and evolution of fertilization proteins is underappreciated. Gene duplication is observed in many fertilization protein families including Izumo, DCST, ZP, and the TFP superfamily. Molecules mediating fertilization are part of larger gene families expressed in a variety of tissues, but gene duplication followed by structural modifications has often facilitated their cooption into a fertilization function. Repeat expansions of functional domains within a gene also provide opportunities for the evolution of novel fertilization protein. ZP proteins with domain repeat expansions are linked to species-specificity in fertilization and TFP proteins that experienced domain duplications were coopted into a novel sperm function. This review outlines the importance of gene duplications and repeat domain expansions in the evolution of fertilization proteins.

Mouse protein coding diversity: What's left to discover?

For over a century, mice have been used to model human disease, leading to many fundamental discoveries about mammalian biology and the development of new therapies. Mouse genetics research has been further catalysed by a plethora of genomic resources developed in the last 20 years, including the genome sequence of C57BL/6J and more recently the first draft reference genomes for 16 additional laboratory strains. Collectively, the comparison of these genomes highlights the extreme diversity that exists at loci associated with the immune system, pathogen response, and key sensory functions, which form the foundation for dissecting phenotypic traits in vivo. We review the current status of the mouse genome across the diversity of the mouse lineage and discuss the value of mice to understanding human disease.

Insight into bovine (Bos indicus) spermatozoal whole transcriptome profile

Mature spermatozoa harbor both coding and non-coding type of RNAs which regulates spermatogenesis, fertilization and early development. Characterization of bovine sperm transcriptome can provide more insight into the molecular mechanisms involved in these processes. Here, we have analyzed whole transcriptome profile of Bos indicus spermatozoa to access the global RNA expression. RNA-Seq analysis identified 14,306 genes expressed with FPKM >0, while 405 genes expressed when threshold increased to FPKM >5. Functional annotations showed that sperm transcripts were associated with molecular processes (translation, ribosomal small and large subunit assembly) and cellular components (cytosolic small and large ribosomal subunit and membranes) related to known sperm functions at fertilization and spermatogenesis. The RNA-Seq data was validated using droplet digital PCR where both highly abundant gene viz. RN7SL1 and less abundant gene viz. ZFP280B were validated. This study may provide future directions in reproductive biology of Bos indicus.

Occurrence and functional significance of the transcriptome in bovine (Bos taurus) spermatozoa

Mammalian spermatozoa deliver various classes of RNAs to the oocyte during fertilization, and many of them may regulate fertility. The objective of the present study was to determine the composition and abundance of spermatozoal transcripts in fresh bull semen. The entire transcriptome of the spermatozoa from bulls (n = 3) was sequenced using two different platforms (Ion Proton and Illumina) to identify the maximum number of genes present in the spermatozoa. The bovine spermatozoa contained transcripts for 13,833 genes (transcripts per million, TPM > 10). Both intact and fragmented transcripts were found. These spermatozoal transcripts were associated with various stages of spermatogenesis, spermatozoal function, fertilization, and embryo development. The presence of intact transcripts of pregnancy-associated glycoproteins (PAGs) in the spermatozoa suggest a possible influence of sperm transcripts beyond early embryonic development. The specific regions (exon, intron, and exon-intron) of the particular spermatozoal transcripts might help regulate fertilization. This study demonstrates that the use of two different RNA-seq platforms provides a comprehensive profile of bovine spermatozoal RNA. Spermatozoal RNA profiling may be useful as a non-invasive method to delineate possible causes of male infertility and to predict fertility in a manner that is more effective than the conventional methods.

Evolutionary study of the isoflavonoid pathway based on multiple copies analysis in soybean

BACKGROUND

Previous studies suggest that the metabolic pathway structure influences the selection and evolution rates of involved genes. However, most of these studies have exclusively considered a single gene copy encoding each enzyme in the metabolic pathway. Considering multiple-copy encoding enzymes could provide direct evidence of gene evolution and duplication patterns in metabolic pathways. We conducted a detailed analysis of the phylogeny, synteny, evolutionary rate and selection pressure of the genes in the isoflavonoid metabolic pathway of soybeans.

RESULTS

The results revealed that 1) only the phenylalanine ammonia-lyase (PAL) gene family most upstream from the pathway preserved all of the ancient and recent segmental duplications and maintained a strongly conserved synteny among these duplicated segments; gene families encoding branch-point enzymes with higher pleiotropy tended to retain more types of duplication; and genes encoding chalcone reductase (CHR) and isoflavone synthase (IFS) specific for legumes retained only recent segmental duplications; 2) downstream genes evolved faster than upstream genes and were subject to positive selection or relaxed selection constraints; 3) gene members encoding enzymes with high pleiotropy at the branching points were more likely to have undergone evolutionary differentiation, which may correspond to their functional divergences.

CONCLUSIONS

We reconciled our results with existing controversies and proposed that gene copies at branch points with higher connectivity might be under stronger selective constraints and that the gene copies controlling metabolic flux allocation underwent positive selection. Our analyses demonstrated that the structure and function of a metabolic pathway shapes gene duplication and the evolutionary constraints of constituent enzymes.

Positive selection in the adhesion domain of Mus sperm Adam genes through gene duplications and function-driven gene complex formations

BACKGROUND

Sperm and testes-expressed Adam genes have been shown to undergo bouts of positive selection in mammals. Despite the pervasiveness of positive selection signals, it is unclear what has driven such selective bouts. The fact that only sperm surface Adam genes show signals of positive selection within their adhesion domain has led to speculation that selection might be driven by species-specific adaptations to fertilization or sperm competition. Alternatively, duplications and neofunctionalization of Adam sperm surface genes, particularly as it is now understood in rodents, might have contributed to an acceleration of evolutionary rates and possibly adaptive diversification.

RESULTS

Here we sequenced and conducted tests of selection within the adhesion domain of sixteen known sperm-surface Adam genes among five species of the Mus genus. We find evidence of positive selection associated with all six Adam genes known to interact to form functional complexes on Mus sperm. A subset of these complex-forming sperm genes also displayed accelerated branch evolution with Adam5 evolving under positive selection. In contrast to our previous findings in primates, selective bouts within Mus sperm Adams showed no associations to proxies of sperm competition. Expanded phylogenetic analysis including sequence data from other placental mammals allowed us to uncover ancient and recent episodes of adaptive evolution.

CONCLUSIONS

The prevailing signals of rapid divergence and positive selection detected within the adhesion domain of interacting sperm Adams is driven by duplications and potential neofunctionalizations that are in some cases ancient (Adams 2, 3 and 5) or more recent (Adams 1b, 4b and 6).

Transcriptional profiles of mating-responsive genes from testes and male accessory glands of the Mediterranean fruit fly, Ceratitis capitata

BACKGROUND

Insect seminal fluid is a complex mixture of proteins, carbohydrates and lipids, produced in the male reproductive tract. This seminal fluid is transferred together with the spermatozoa during mating and induces post-mating changes in the female. Molecular characterization of seminal fluid proteins in the Mediterranean fruit fly, Ceratitis capitata, is limited, although studies suggest that some of these proteins are biologically active.

METHODOLOGY/PRINCIPAL FINDINGS

We report on the functional annotation of 5914 high quality expressed sequence tags (ESTs) from the testes and male accessory glands, to identify transcripts encoding putative secreted peptides that might elicit post-mating responses in females. The ESTs were assembled into 3344 contigs, of which over 33% produced no hits against the nr database, and thus may represent novel or rapidly evolving sequences. Extraction of the coding sequences resulted in a total of 3371 putative peptides. The annotated dataset is available as a hyperlinked spreadsheet. Four hundred peptides were identified with putative secretory activity, including odorant binding proteins, protease inhibitor domain-containing peptides, antigen 5 proteins, mucins, and immunity-related sequences. Quantitative RT-PCR-based analyses of a subset of putative secretory protein-encoding transcripts from accessory glands indicated changes in their abundance after one or more copulations when compared to virgin males of the same age. These changes in abundance, particularly evident after the third mating, may be related to the requirement to replenish proteins to be transferred to the female.

CONCLUSIONS/SIGNIFICANCE

We have developed the first large-scale dataset for novel studies on functions and processes associated with the reproductive biology of Ceratitis capitata. The identified genes may help study genome evolution, in light of the high adaptive potential of the medfly. In addition, studies of male recovery dynamics in terms of accessory gland gene expression profiles and correlated remating inhibition mechanisms may permit the improvement of pest management approaches.

Phylogenetic and molecular evolution of the ADAM (A Disintegrin And Metalloprotease) gene family from Xenopus tropicalis, to Mus musculus, Rattus norvegicus, and Homo sapiens

ADAM (a disintegrin and metalloprotease) genes have been identified in various tissues and species, and recently associated with several important human diseases such as tumor and asthma. Although various biological processes have been known for the ADAM family in different species including fertilization, neurogenesis, infection and inflammation, little is known about its detailed phylogenetic and molecular evolutionary history. In this study, the ADAMs of Xenopus (Silurana) tropicalis, Mus musculus, Rattus norvegicus, and Homo sapiens were collected and analyzed by using the Bayesian analysis and gene synteny analysis to establish a comprehensive phylogenetic relationship and evolutionary drive of this gene family. It was found that there were more ADAMs in the two rodents than in the amphibian, suggesting an expansion of the ADAM gene family during the early evolution of mammals. All ADAMs from this expansion were retained in both the rodents, but other duplication events occurred subsequently in the two rodents, respectively, leading to the classification of rodent ADAMs as classes I, II and III. Moreover, these duplicated ADAM genes in the rodents were found to be driven by positive selection, which might be the major force to retain them in the genome. Importantly, it was also found that orthologs of ADAM3 and 5 have been lost in humans. These results not only provide valuable information of the evolution of ADAM genes, but may also help in understanding the role of ADAM genes in the pathobiology of relevant diseases.

Sexual selection and the molecular evolution of ADAM proteins

Rapid evolution has been identified for many reproductive genes and recent studies have combined phylogenetic tests and information on species mating systems to test sexual selection. Here we examined the molecular evolution of the ADAM gene family, a diverse group of 35 proteins capable of adhesion to and cleavage of other proteins, using sequence data from 25 mammalian genes. Out of the 25 genes analyzed, all those expressed in male reproductive tissue showed evidence of positive selection. Positively selected amino acids within the protein adhesion domain were only found in sperm surface ADAM proteins (ADAMs 1, 2, 3, 4, and 32) suggesting selection driven by male x female interactions. We tested heterogeneity in rates of evolution of the adhesion domain of ADAM proteins by using sequence data from Hominidae and macaques. The use of the branch and branch-site models (PAML) showed evidence of higher d (N)/d (S) and/or positive selection linked to branches experiencing high postmating selective pressures (chimpanzee and macaque) for Adams 2, 18, and 23. Moreover, we found consistent higher proportion of nonsynonymous relative to synonymous and noncoding sequence substitutions in chimpanzee and/or macaque only for Adams 2, 18, and 23. Our results suggest that lineage-specific sexual selection bouts might have driven the evolution of the adhesion sperm protein surface domains of ADAMs 2 and 18 in primates. Adams 2 and 18 are localized in chromosome 8 of primates and adjacent to each other, so their evolution might have also been influenced by their common genome localization.

Extreme positive selection on a new highly-expressed larval glycoprotein (LGP) gene in Galaxias fishes (Osmeriformes: Galaxiidae)

We describe the intron-exon structure and DNA/protein sequences of a new larval glycoprotein (LGP) gene from nine species of galaxiid fish. The gene has a distant similarity to Danio THP (Tamm-Horsfall urinary glycoprotein; uromodulin) and cichlid SPP120 (seminal plasma glycoprotein) due to conserved features of its zona pellucida (ZP) domain, including eight highly conserved cysteines and a consensus furin cleavage site. Using a combination of 454 sequencing of cDNA and exon-primed intron-spanning sequencing of genomic DNA, we obtained full sequences of the coding region (996 bp) and its intervening sequences (1,459 bp). LGP shows an exceptionally strong signal of positive selection over the entire coding region, as evidenced by d(N)/d(S) values >1. Across nine species of Galaxias, 87/332 (26%) amino acid residues are variable, compared with 9/386 (2%) for mitochondrial cytochrome b (cytb) in the same group of species. Across 36 interspecific pairwise comparisons, genetic distances are in all cases larger for coding region than for introns, by a factor of 2.4-fold on average. Reading frame, gene structure, splice sites, and many ZP motifs are conserved across all species. Together with the fact that the gene is expressed in all species, these results argue clearly against the possibility of a pseudogene. We show by 454 sequencing and quantitative polymerase chain reaction that the transcript is abundant (ca. 0.5%) in newly hatched larvae and appears to be almost absent from a range of adult tissues. We postulate that the strong Darwinian evolution exhibited by this protein may reflect some type of immunoprotection at this vulnerable larval stage.

Sperm proteomics reveals intensified selection on mouse sperm membrane and acrosome genes

Spermatozoa are a focal point for the impact of sexual selection due to sperm competition and sperm-female interactions in a wide range of sexually reproducing organisms. In-depth molecular investigation of the ramifications of these selective regimes has been limited due to a lack of information concerning the molecular composition of sperm. In this study, we utilize three previously published proteomic data sets in conjunction with our whole mouse sperm proteomic analysis to delineate cellular regions of sperm most impacted by positive selection. Interspecific analysis reveals robust evolutionary acceleration of sperm cell membrane genes (which include genes encoding acrosomal and sperm cell surface proteins) relative to other sperm genes, and evidence for positive selection in approximately 22% of sperm cell membrane components was obtained using maximum likelihood models. The selective forces driving the accelerated evolution of these membrane proteins may occur at a number of locations during sperm development, maturation, and transit through the female reproductive tract where the sperm cell membrane and eventually the acrosome are exposed to the extracellular milieu and available for direct cell-cell interactions.

Identification, RNAi knockdown, and functional analysis of an ejaculate protein that mediates a postmating, prezygotic phenotype in a cricket

Postmating, prezygotic phenotypes, especially those that underlie reproductive isolation between closely related species, have been a central focus of evolutionary biologists over the past two decades. Such phenotypes are thought to evolve rapidly and be nearly ubiquitous among sexually reproducing eukaryotes where females mate with multiple partners. Because these phenotypes represent interplay between the male ejaculate and female reproductive tract, they are fertile ground for reproductive senescence – as ejaculate composition and female physiology typically change over an individual's life span. Although these phenotypes and their resulting dynamics are important, we have little understanding of the proteins that mediate these phenotypes, particularly for species groups where postmating, prezygotic traits are the primary mechanism of reproductive isolation. Here, we utilize proteomics, RNAi, mating experiments, and the Allonemobius socius complex of crickets, whose members are primarily isolated from one another by postmating, prezygotic phenotypes (including the ability of a male to induce a female to lay eggs), to demonstrate that one of the most abundant ejaculate proteins (a male accessory gland-biased protein similar to a trypsin-like serine protease) decreases in abundance over a male's reproductive lifetime and mediates the induction of egg-laying in females. These findings represent one of the first studies to identify a protein that plays a role in mediating both a postmating, prezygotic isolation pathway and reproductive senescence.

Adaptive evolution of gamete-recognition proteins in birds

Gamete-recognition proteins have been shown to evolve by positive selection in diverse organism groups, such as marine invertebrates and mammals, although underlying evolutionary mechanisms driving this rapid divergence are poorly understood. However, several hypotheses have been put forward to explain the observed pattern, including different forms of sexual conflict and sperm competition. Because female gametes require more energy to produce than male gametes, female organisms suffer more when fertilisation goes wrong. One process that results in a failed mammalian fertilisation is polyspermy, when >1 sperm fertilises the egg. However in birds, there is no such sexual conflict because multiple sperm typically bind and fuse with the egg. If sexual conflict driven by polyspermy avoidance is important for the evolution of gamete-recognition proteins in vertebrates, we expect to find positive selection in the genes to be less pronounced in birds. We therefore sequenced six genes (ZP1, ZP2, ZP4, ZPAX, CD9, and Acrosin) encoding gamete-recognition proteins in several bird species to test for positive selection. For comparison, we also analysed ortologous sequences in a set of mammalian species. We found no major differences in the occurrence of adaptive evolution and the strength of selection between bird and mammal orthologs. From this we conclude that polyspermy avoidance does not act as the main underlying evolutionary force shaping the rate of evolution in these genes. We discuss other possible processes that could explain positive selection of gamete-recognition proteins in birds and mammals, such as hybridisation avoidance, cryptic female choice, and postcopulatory sperm competition.

Evidence of adaptive evolution of accessory gland proteins in closely related species of the Drosophila repleta group

Accessory gland proteins (Acps) are part of the seminal fluid of Drosophila species. These proteins have important reproductive functions, being responsible for the proper functioning of several steps of the fertilization process. Acps also contribute indirectly for the reproductive success of males by modulating female behavior. Evidence that Acps participate in sperm competition and sexual conflict includes findings that, on average, Acps have fast evolutionary rates, suggestive of adaptive evolution. This is especially true in species of the Drosophila repleta group. Nevertheless, only in a few occasions have robust statistical tests been used to determine whether observed evolutionary rates are in fact due to positive selection on amino acid substitutions between related species. Here we apply maximum likelihood tests for positive selection on 14 Acps of the D. repleta group. To increase statistical robustness, we use at least 8 sequences, all belonging to species of the Drosophila mulleri complex, for each gene analyzed. We found significant evidence of adaptive evolution for 10 of the tested genes. Among these, the ones with a conserved protein domain had positively selected sites within the functional region of the sequence. We also detected one instance of lineage-specific adaptive evolution in a clade formed by 2 sister species.

Diversifying selection and host adaptation in two endosymbiont genomes

BACKGROUND

The endosymbiont Wolbachia pipientis infects a broad range of arthropod and filarial nematode hosts. These diverse associations form an attractive model for understanding host:symbiont coevolution. Wolbachia's ubiquity and ability to dramatically alter host reproductive biology also form the foundation of research strategies aimed at controlling insect pests and vector-borne disease. The Wolbachia strains that infect nematodes are phylogenetically distinct, strictly vertically transmitted, and required by their hosts for growth and reproduction. Insects in contrast form more fluid associations with Wolbachia. In these taxa, host populations are most often polymorphic for infection, horizontal transmission occurs between distantly related hosts, and direct fitness effects on hosts are mild. Despite extensive interest in the Wolbachia system for many years, relatively little is known about the molecular mechanisms that mediate its varied interactions with different hosts. We have compared the genomes of the Wolbachia that infect Drosophila melanogaster, wMel and the nematode Brugia malayi, wBm to that of an outgroup Anaplasma marginale to identify genes that have experienced diversifying selection in the Wolbachia lineages. The goal of the study was to identify likely molecular mechanisms of the symbiosis and to understand the nature of the diverse association across different hosts.

RESULTS

The prevalence of selection was far greater in wMel than wBm. Genes contributing to DNA metabolism, cofactor biosynthesis, and secretion were positively selected in both lineages. In wMel there was a greater emphasis on DNA repair, cell division, protein stability, and cell envelope synthesis.

CONCLUSION

Secretion pathways and outer surface protein encoding genes are highly affected by selection in keeping with host:parasite theory. If evidence of selection on various cofactor molecules reflects possible provisioning, then both insect as well as nematode Wolbachia may be providing substances to hosts. Selection on cell envelope synthesis, DNA replication and repair machinery, heat shock, and two component switching suggest strategies insect Wolbachia may employ to cope with diverse host and intra-host environments.

The evolution of the vertebrate metzincins; insights from Ciona intestinalis and Danio rerio

BACKGROUND

The metzincins are a large gene superfamily of proteases characterized by the presence of a zinc protease domain, and include the ADAM, ADAMTS, BMP1/TLL, meprin and MMP genes. Metzincins are involved in the proteolysis of a wide variety of proteins, including those of the extracellular matrix. The metzincin gene superfamily comprises eighty proteins in the human genome and ninety-three in the mouse. When and how the level of complexity apparent in the vertebrate metzincin gene superfamily arose has not been determined in detail. Here we present a comprehensive analysis of vertebrate metzincins using genes from both Ciona intestinalis and Danio rerio to provide new insights into the complex evolution of this gene superfamily.

RESULTS

We have identified 19 metzincin genes in the ciona genome and 83 in the zebrafish genome. Phylogenetic analyses reveal that the expansion of the metzincin gene superfamily in vertebrates has occurred predominantly by the simple duplication of pre-existing genes rather than by the appearance and subsequent expansion of new metzincin subtypes (the only example of which is the meprin gene family). Despite the number of zebrafish metzincin genes being relatively similar to that of tetrapods (e.g. man and mouse), the pattern of gene retention and loss within these lineages is markedly different. In addition, we have studied the evolution of the related TIMP gene family and identify a single ciona and four zebrafish TIMP genes.

CONCLUSION

The complexity seen in the vertebrate metzincin gene families was mainly acquired during vertebrate evolution. The metzincin gene repertoire in protostomes and invertebrate deuterostomes has remained relatively stable. The expanded metzincin gene repertoire of extant tetrapods, such as man, has resulted largely from duplication events associated with early vertebrate evolution, prior to the sarcopterygian-actinopterygian split. The teleost repertoire of metzincin genes in part parallels that of tetrapods but has been significantly modified, perhaps as a consequence of a teleost-specific duplication event.

Rapid evolution of the trophoblast kunitz domain proteins (TKDPs)-a multigene family in ruminant ungulates

The trophoblast Kunitz domain proteins (TKDPs) are products of the outer cells (trophoblasts) of the placenta of cattle, sheep, and related species. Most are expressed abundantly for only a few days during the time at which the ruminant conceptus is first establishing intimate contacts with the uterine lining. The TKDPs are secretory proteins that possess a carboxyl-terminal peptidase inhibitory domain related to the Kunitz family of serine peptidase inhibitors. On the amino-terminal end are one or more highly unusual regions that are unique to the TKDP genes and have no apparent similarity to any other known sequences. The TKDPs are a rather divergent family that exhibits a good deal of variation among the members. To better understand the reason for such variation, the rates of synonymous (dS) and nonsynonymous (dN), as well as radical (pNR) and conservative (pNC), substitutions were assessed. Phylogenetic trees revealed that the Kunitz domains represented three related groups, whereas the amino-terminal domains formed four groupings. Pairwise comparisons between Kunitz and amino-terminal domain groups demonstrated that dN was consistently greater than dS. In addition, nonsynonymous substitutions in the Kunitz domains tended to be radical (changing charge or polarity), while those in the amino-terminal domains exhibited neither a preponderance of conservative nor radical substitution rates. In summary, the rapid evolution of the TKDPs, coupled with their restricted temporal expression during development, likely reflects the establishment of protein-protein interactions that have evolved to serve the unusual synepitheliochorial placenta of ruminant ungulates.

Positive selection in the evolution of cancer

We hypothesize that forms of antagonistic coevolution have forged strong links between positive selection at the molecular level and increased cancer risk. By this hypothesis, evolutionary conflict between males and females, mothers and foetuses, hosts and parasites, and other parties with divergent fitness interests has led to rapid evolution of genetic systems involved in control over fertilization and cellular resources. The genes involved in such systems promote cancer risk as a secondary effect of their roles in antagonistic coevolution, which generates evolutionary disequilibrium and maladaptation. Evidence from two sources: (1) studies on specific genes, including SPANX cancer/testis antigen genes, several Y-linked genes, the pem homebox gene, centromeric histone genes, the breast cancer gene BRCA1, the angiogenesis gene ANG, cadherin genes, cytochrome P450 genes, and viral oncogenes; and (2) large-scale database studies of selection on different functional categories of genes, supports our hypothesis. These results have important implications for understanding the evolutionary underpinnings of cancer and the dynamics of antagonistically-coevolving molecular systems.

Identification of a positively evolving putative binding region with increased variability in posttranslational motifs in zonadhesin MAM domain 2

Positive selection has been shown to be pervasive in sex-related proteins of many metazoan taxa. However, we are only beginning to understand molecular evolutionary processes on the lineage to humans. To elucidate the evolution of proteins involved in human reproduction, we studied the sequence evolution of MAM domains of the sperm-ligand zonadhesin in respect to single amino acid sites, solvent accessibility, and posttranslational modification. GenBank-data were supplemented by new cDNA-sequences of a representative non-human primate panel. Solvent accessibility predictions identified a probably exposed fragment of 30 amino acids belonging to MAM domain 2 (i.e., MAM domain 3 in mouse). The fragment is characterized by significantly increased rate of positively selected amino acid sites and exhibits high variability in predicted posttranslational modification, and, thus, might represent a binding region in the mature protein. At the same time, there is a significant coincidence of positively selected amino acid sites and non-conserved posttranslational motifs. We conclude that the binding specificity of zonadhesin MAM domains, especially of the presumed epitope, is achieved by positive selection at the level of single amino acid sites and posttranslational modifications, respectively.

Rod and cone opsin families differ in spectral tuning domains but not signal transducing domains as judged by saturated evolutionary trace analysis

The visual receptor of rods and cones is a covalent complex of the apoprotein, opsin, and the light-sensitive chromophore, 11-cis-retinal. This pigment must fulfill many functions including photoactivation, spectral tuning, signal transmission, inactivation, and chromophore regeneration. Rod and cone photoreceptors employ distinct families of opsins. Although it is well known that these opsin families provide unique ranges in spectral sensitivity, it is unclear whether the families have additional functional differences. In this study, we use evolutionary trace (ET) analysis of 188 vertebrate opsin sequences to identify functionally important sites in each opsin family. We demonstrate the following results. (1) The available vertebrate opsin sequences produce a definitive description of all five vertebrate opsin families. This is the first demonstration of sequence saturation prior to ET analysis, which we term saturated ET (SET). (2) The cone opsin classes have class-specific sites compared to the rod opsin class. These sites reside in the transmembrane region and tune the spectral sensitivity of each opsin class to its characteristic wavelength range. (3) The cytoplasmic loops, primarily responsible for signal transmission and inactivation, are essentially invariant in rod versus cone opsins. This indicates that the electrophysiological differences between rod and cone photoreceptors cannot be ascribed to differences in the protein interaction regions of the opsins. SET shows that chromophore binding and regeneration are the only aspects of opsin structure likely to have functionally significant differences between rods and cones, whereas excitatory and adaptational properties of the opsin families appear to be functionally invariant.

Adaptive Molecular Evolution in the Opsin Genes of Rapidly Speciating Cichlid Species

Cichlid fish inhabit a diverse range of environments that vary in the spectral content of light available for vision. These differences should result in adaptive selective pressure on the genes involved in visual sensitivity, the opsin genes. This study examines the evidence for differential adaptive molecular evolution in East African cichlid opsin genes due to gross differences in environmental light conditions. First, we characterize the selective regime experienced by cichlid opsin genes using a likelihood ratio test format, comparing likelihood models with different constraints on the relative rates of amino acid substitution, across sites. Second, we compare turbid and clear lineages to determine if there is evidence of differences in relative rates of substitution. Third, we present evidence of functional diversification and its relationship to the photic environment among cichlid opsin genes. We report statistical evidence of positive selection in all cichlid opsin genes, except short wavelength-sensitive 1 and short wavelength-sensitive 2b. In all genes predicted to be under positive selection, except short wavelength-sensitive 2a, we find differences in selective pressure between turbid and clear lineages. Potential spectral tuning sites are variable among all cichlid opsin genes; however, patterns of substitution consistent with photic environment-driven evolution of opsin genes are observed only for short wavelength-sensitive 1 opsin genes. This study identifies a number of promising candidate-tuning sites for future study by site-directed mutagenesis. This work also begins to demonstrate the molecular evolutionary dynamics of cichlid visual sensitivity and its relationship to the photic environment.

Adaptive evolution of the insulin gene in caviomorph rodents

Insulin is a conservative molecule among mammals, maintaining both its structure and function. Rodents that belong to the Suborder Hystricognathi represent an exception, having a very divergent molecule with unusual physiological properties. In this work, we analyzed the evolutionary pattern of the insulin gene in caviomorph rodents (South American hystricomorph rodents). We found that these rodents have higher rates of nonsynonymous:synonymous substitutions (d(N)/d(S)) than nonhystricomorph rodents and that values are heterogeneous inside the group. We estimated codons under positive selection, specifically the second binding site (A13 and B17) and others related with hexamerization (B18, B20, and B22). In the monomer structure, all selected sites formed a single patch around the second binding site. In the hexamer structure, these amino acids were grouped into three major patches. In this structure, contacts between B chains involved all selected sites (except B18), and between faces in the center of the molecule, all contacts were among selected sites. While there is no clear hypothesis regarding the cause of this drastic change, experimental evidence does show that this group of rodents has some peculiarities in growth function, and, whether coincidental or not, these changes appeared together with important changes in life-history traits.

Positive selection at reproductive ADAM genes with potential intercellular binding activity

Many genes with a role in reproduction, including those implicated in fertilization and spermatogenesis, have been shown to evolve at a faster rate relative to genes associated with other functions and tissues. These survey studies usually group a wide variety of genes with different characteristics and evolutionary histories as reproductive genes based on their site of expression or function. We have examined the molecular evolution of the ADAM (a disintegrin and metalloprotease) gene family, a structurally and functionally diverse group of genes expressed in reproductive and somatic tissue to test whether a variety of protein characteristics such as phylogenetic clusters, tissue of expression, and proteolytic and adhesive function can group fast evolving ADAM genes. We found that all genes were evolving under purifying selection (d(N)/d(S) < 1), although reproductive ADAMs, including those implicated in fertilization and spermatogenesis, evolved at the fastest rate. Genes with a role in binding to cell receptors in endogenous tissue appear to be evolving under purifying selection, regardless of the tissue of expression. In contrast, positive selection of codon sites in the disintegrin/cysteine-rich adhesion domains was detected exclusively in ADAMs 2 and 32, two genes expressed in the testis with a potential role in sperm-egg adhesion. Positive selection was detected in the transmembrane/cytosolic tail region of ADAM genes expressed in a variety of tissues.

Natural selection and evolution of streptococcal virulence genes involved in tissue-specific adaptations

The molecular mechanisms underlying niche adaptation in bacteria are not fully understood. Primary infection by the pathogen group A streptococcus (GAS) takes place at either the throat or the skin of its human host, and GAS strains differ in tissue site preference. Many skin-tropic strains bind host plasminogen via the plasminogen-binding group A streptococcal M protein (PAM) present on the cell surface; inactivation of genes encoding either PAM or streptokinase (a plasminogen activator) leads to loss of virulence at the skin. Unlike PAM, which is present in only a subset of GAS strains, the gene encoding streptokinase (ska) is present in all GAS isolates. In this study, the evolution of the virulence genes known to be involved in skin infection was examined. Most genetic diversity within ska genes was localized to a region encoding the plasminogen-docking domain (beta-domain). The gene encoding PAM displayed strong linkage disequilibrium (P << 0.01) with a distinct phylogenetic cluster of the ska beta-domain-encoding region. Yet, ska alleles of distant taxa showed a history of intragenic recombination, and high intrinsic levels of recombination were found among GAS strains having different tissue tropisms. The data suggest that tissue-specific adaptations arise from epistatic coselection of bacterial virulence genes. Additional analysis of ska genes showed that approximately 4% of the codons underwent strong diversifying selection. Horizontal acquisition of one ska lineage from a commensal Streptococcus donor species was also evident. Together, the data suggest that new phenotypes can be acquired through interspecies recombination between orthologous genes, while constrained functions can be preserved; in this way, orthologous genes may provide a rich and ready source for new phenotypes and thereby play a facilitating role in the emergence of new niche adaptations in bacteria.

Rapid evolution of a pollen-specific oleosin-like gene family from Arabidopsis thaliana and closely related species

It has been shown in a variety of species that genes expressed in reproductive tissues evolve rapidly, which often appears to be the result of positive Darwinian selection. We investigated the evolution of a family of seven pollen-specific oleosin-like proteins (or oleopollenins) in Arabidopsis thaliana and two closely related species. More than 30 kb of a genomic region that harbors the complete, tandemly repeated oleopollenin cluster were sequenced from Arabidopsis lyrata ssp. lyrata, and Boechera drummondii. A phylogenetic analysis of the complete gene cluster from these three species and from Brassica oleracea confirmed its rapid evolution resulting from gene duplication and gene loss events, numerous amino acid substitutions, and insertions/deletions in the coding sequence. Independent duplications were inferred in the lineages leading to Arabidopsis and to Brassica, and gene loss was inferred in the lineage leading to B. drummondii. Comparisons of the ratio of nonsynonymous (d(N)) and synonymous (d(S)) divergence revealed that the oleopollenins are among the most rapidly evolving proteins currently known from Arabidopsis and that they may evolve under positive Darwinian selection. Reverse transcriptase polymerase chain reaction analysis demonstrated the expression of oleopollenins in flowers of the outcrossing A. lyrata, the selfing B. drummondii, and the apomictic Boechera holboellii, suggesting that oleopollenins play an important role in species with different breeding systems. These results are consistent with a putative function in species recognition, but further analyses of protein function and sequence variation in species with different breeding systems are necessary to reveal the underlying causes for the rapid evolution of oleopollenins.

Shall we dance or shall we fight? Using DNA sequence data to untangle controversies surrounding sexual selection

Population and evolutionary genetics studies have largely benefitted from advances in DNA manipulation and sequencing, as well as DNA data analysis techniques. Molecular evolution studies of male reproductive genes show a pattern of rapid evolution shaped, in some cases, by an adaptive selective process. Despite the large body of data on male reproductive genes, the female side of the story has remained unexplored. The few cases of female egg receptors analyzed also show rapid evolution. However, to disentangle between competing hypotheses on how selection operates on male × female molecular interaction leading to fertilization, we need to find male and female molecules that are partners in fertilization. A conflict model of sexual selection (similar to a host-parasite model) assumes a male-driven system where females are being forced under suboptimal conditions. This predicts that the amount of divergence at a female receptor depends on the amount of divergence among the male reproductive proteins that it binds (i.e., males are leading). Under a classical model of runaway sexual selection, female protein receptors might be the key to the rapid molecular changes observed in male reproductive proteins and higher divergence should be expected among female receptors than among their respective male binding proteins.Key words: Reproductive genes, DNA sequence data, sexual selection, coadaptation, conflict.

Adaptive evolution of genes and gene families

The huge influx of genomic sequence and new statistical methods is making the discovery of genes subjected to adaptive evolution increasingly common. The use of comparative genomics to identify adaptive evolution is resulting in predictions of functionally important genes and gene regions. However, the selective pressure driving the adaptive evolution of most genes remains mysterious.