The molecular evolution of sperm zonadhesin

Based on pioneering work of Hardy and Garbers, zonadhesin has become one of the best studied sperm ligands in boreoeutherian mammals, both from a biochemical and evolutionary perspective. Zonadhesin is a mosaic-type protein that localizes to the apical head of spermatozoa. In pig, cattle, rabbit and primates, zonadhesin precursor essentially consists of two or three MAM (meprin/A5 antigen/mu receptor tyrosine phosphatase) domains, one mucin-like domain, one incomplete and four complete D domains (homologous to vWFD). Mouse zonadhesin is distinguished from this general pattern by 20 extra partial D3 domains. While concerted evolution drives the divergence of the mucin-like domain in the ortholog comparison, MAM and D domains mainly diverge under the influence of drift and positive selection, both in the paralog and ortholog comparison. As can be seen particularly well within a putative binding region in the most C-terminal MAM domain, positive selection not only causes amino acid exchanges, but also promotes changes in the pattern of predicted posttranslational modification. Moving window and correlation analyses of sequence evolution and sexual body dimorphism further suggest that sexual selection, especially sperm competition, drives zonadhesin divergence. However, considering its zona pellucida avidity, female cryptic choice might as well contribute to zonadhesin evolution. Despite the general tendency for divergence of zonadhesin, conservation by negative selection dominates the evolution of most codon sites. In accordance, the distribution of EGF (epidermal growth factor)-like motifs, DP-doublets, single cysteines and CGLC motifs suggests a wide conservation of processing, folding and oligomerization of zonadhesin in pig, rabbit and primates.     

Molecular population genetics of the gene encoding the human fertilization protein zonadhesin reveals rapid adaptive evolution

A hallmark of positive selection (adaptive evolution) in protein-coding regions is a d(N)/d(S) ratio >1, where d(N) is the number of nonsynonymous substitutions/nonsynonymous sites and d(S) is the number of synonymous substitutions/synonymous sites. Zonadhesin is a male reproductive protein localized on the sperm head, comprising many domains known to be involved in cell-cell interaction or cell adhesion. Previous studies have shown that VWD domains (homologous to the D domains of the von Willebrand factor) are involved directly in binding to the female zona pellucida (ZP) in a species-specific manner. In this study, we sequenced 47 coding exons in 12 primate species and, by using maximum-likelihood methods to determine sites under positive selection, we show that VWD2, membrane/A5 antigen mu receptor, and mucin-like domains in zonadhesin are rapidly evolving and, thus, may be involved in binding to the ZP in a species-specific manner in primates. In addition, polymorphism data from 48 human individuals revealed significant polymorphism-to-divergence heterogeneity and a significant departure from equilibrium-neutral expectations in the frequency spectrum, suggesting balancing selection and positive selection occurring in zonadhesin (ZAN) within human populations. Finally, we observe adaptive evolution in haplotypes segregating for a frameshift mutation that was previously thought to indicate that ZAN was a potential pseudogene.     

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

Genes with a role in fertilization show a common pattern of rapid evolution. The role played by positive selection versus lack of selective constraints has been more difficult to establish. One problem arises from attempts to detect selection in an overall gene sequence analysis. I have analyzed the pattern of molecular evolution of fertilin, a gene coding for a heterodimeric sperm protein belonging to the ADAM (A disintegrin and A metalloprotease) gene family. A nonsynonymous to synonymous rate ratio (d(N)/d(S)) analysis for different protein domains of fertilin alpha and fertilin beta showed d(N)/d(S) < 1, suggesting that purifying selection has shaped fertilin's evolution. However, an analysis of the distribution of single positively selected codon sites using phylogentic analysis by maximum likelihood (PAML) showed sites within adhesion domains (disintegrin and cysteine-rich) of fertilin beta evolving under positive selection. The region 3' to the EGF-like domain of fertilin alpha, where the transmembrane and cytoplasmic tail regions are supposed to be localized, showed higher d(N) and d(S) than any other fertilin alpha region. However, it was not possible to identify positively selected codon sites due to ambiguous alignments of the carboxy-end region (ClustalX vs. DiAlign2). When this region was excluded from the PAML analysis, most single positively selected codon sites were concentrated within adhesion domains (cysteine-rich and EGF-like). The use of an ancestral sequence prior to a recent duplication event of fertilin alpha among non-Hominidae primates (Macaca, Papio, and Saguinus) revealed that the duplication is partially responsible for masking the detection of positively selected sites within the disintegrin domain. Finally, most ADAM genes with a potential role in sperm maturation and/or fertilization showed significantly higher d(N) estimates than other ADAM genes.     

Prediction of protein modification sites of pyrrolidone carboxylic acid using mRMR feature selection and analysis

Pyrrolidone carboxylic acid (PCA) is formed during a common post-translational modification (PTM) of extracellular and multi-pass membrane proteins. In this study, we developed a new predictor to predict the modification sites of PCA based on maximum relevance minimum redundancy (mRMR) and incremental feature selection (IFS). We incorporated 727 features that belonged to 7 kinds of protein properties to predict the modification sites, including sequence conservation, residual disorder, amino acid factor, secondary structure and solvent accessibility, gain/loss of amino acid during evolution, propensity of amino acid to be conserved at protein-protein interface and protein surface, and deviation of side chain carbon atom number. Among these 727 features, 244 features were selected by mRMR and IFS as the optimized features for the prediction, with which the prediction model achieved a maximum of MCC of 0.7812. Feature analysis showed that all feature types contributed to the modification process. Further site-specific feature analysis showed that the features derived from PCA's surrounding sites contributed more to the determination of PCA sites than other sites. The detailed feature analysis in this paper might provide important clues for understanding the mechanism of the PCA formation and guide relevant experimental validations.     

Adaptive evolution of fertilization proteins within a genus: variation in ZP2 and ZP3 in deer mice (Peromyscus)

Rapid evolution of reproductive proteins has been documented in a wide variety of taxa. In internally fertilized species, knowledge about the evolutionary dynamics of these proteins between closely related taxa is primarily limited to accessory gland proteins in the semen of Drosophila. Investigation of additional taxa and functional classes of proteins is necessary in order to determine if there is a general pattern of adaptive evolution of reproductive proteins between recently diverged species. We performed an evolutionary analysis of 2 egg coat proteins, ZP2 and ZP3, in 15 species of deer mice (genus Peromyscus). Both of these proteins are involved in egg-sperm binding, a critical step in maintaining species-specific fertilization. Here, we show that Zp2 and Zp3 gene trees are not consistent with trees based on nonreproductive genes, Mc1r and Lcat, where species formed monophyletic clades. In fact, for both of the reproductive genes, intraspecific amino acid variation was extensive and alleles were sometimes shared across species. We document positive selection acting on ZP2 and ZP3 and identify specific amino acid sites that are likely targets of selection using both maximum likelihood approaches and patterns of parallel amino acid change. In ZP3, positively selected sites are clustered in and around the region implicated in sperm binding in Mus, suggesting changes may impact egg-sperm binding and fertilization potential. Finally, we identify lineages with significantly elevated rates of amino acid substitution using a Bayesian mapping approach. These findings demonstrate that the pattern of adaptive reproductive protein evolution found at higher taxonomic levels can be documented between closely related mammalian species, where reproductive isolation has evolved recently.     

Positive selection drives rapid evolution of certain amino acid residues in an evolutionarily highly conserved interferon-inducible antiviral protein of fishes

Viperin, an evolutionarily highly conserved interferon-inducible multifunctional protein, has previously been reported to exhibit antiviral activity against a wide range of DNA and RNA viruses. Utilizing the complete nucleotide coding sequence data of fish viperin antiviral genes, and employing the maximum likelihood-based codon substitution models, the present study reports the pervasive role of positive selection in the evolution of viperin antiviral protein in fishes. The overall rate of nonsynonymous (dN) to synonymous (dS) substitutions (dN/dS) for the three functional domains of viperin (N-terminal, central domain and C-terminal) were 1.1, 0.12, and 0.24, respectively. Codon-by-codon substitution analyses have revealed that while most of the positively selected sites were located at the N-terminal amphipathic α-helix domain, few amino acid residues at the C-terminal domain were under positive selection. However, none of the sites in the central domain were under positive selection. These results indicate that, although viperin is evolutionarily highly conserved, the three functional domains experienced differential selection pressures. Taken together with the results of previous studies, the present study suggests that the persistent antagonistic nature of surrounding infectious viral pathogens might be the likely cause for such adaptive evolutionary changes of certain amino acids in fish viperin antiviral protein.     

Tandem Repetitive D Domains of the Sperm Ligand Zonadhesin Evolve Faster in the Paralogue Than in the Orthologue Comparison

Gene duplication is regarded as an important evolutionary mechanism creating genetic and phenotypic novelty. At the same time, the evolutionary mechanisms following gene duplication have been a subject of much debate. Here we analyze the sequence evolution of zonadhesin, a mammalian sperm ligand that binds to the oocyte zona pellucida in a species-specific manner. In pig, rabbit, and primates, precursor zonadhesin comprises, among others, one partial and four complete tandem repetitive D domains. The mouse precursor is distinguished by 20 additional partial D3 domains consisting of 120 amino acids each. This gene structure allows sequence comparison in both paralogues and orthologues. Detailed sequence analysis reveals that D domains evolve faster across paralogues than orthologues. Moreover, at the codon level, partial D3 paralogues of mouse show evidence of positive selection, whereas the corresponding orthologues do not. Individual posttranslational motif patterns and positive selection point to neofunctionalization of partial D3 paralogues of mouse, rather than subfunctionalization. However, as we found additional evidence for homogenization by partial gene conversion, sequence evolution of partial D3 paralogues of mouse might be better described as a combination of divergent and convergent evolution. So far, the divergence at the codon level has outbalanced the convergence at the level of smaller fragments. The probable driving force behind the evolutionary patterns observed is sexual selection. We finally discuss whether the functional determination influences the evolutionary regime acting on sperm ligands and egg receptors, respectively. [Reviewing Editor: Dr. Yves Van de Peer]     

The Wilhelmine E. Key 1999 Invitational lecture. Predicting the evolution of human influenza A

We studied the evolution of the HA1 domain of the H3 hemagglutinin gene from human influenza virus type A. The phylogeny of these genes showed a single dominant lineage persisting over time. We tested the hypothesis that the progenitors of this single evolutionarily successful lineage were viruses carrying mutations at codons at which prior mutations had helped the virus to avoid human immune surveillance. We found evidence that eighteen hemagglutinin codons appeared to have been under positive selection to change the amino acid they encoded in the past. Retrospective tests show that viral lineages undergoing the greatest number of mutations in the positively selected codons were the progenitors of future H3 lineages in nine of eleven recent influenza seasons. Codons under positive selection were associated with antibody combining sites A or B or the sialic acid receptor binding site. However, not all codons in these sites had predictive value. Monitoring new H3 isolates for additional changes in positively selected codons might help identify the most fit extant viral strains that arise during antigenic drift.     

Accelerated evolution and loss of a domain of the sperm-egg-binding protein SED1 in ancestral primates

Proteins involved in sperm-egg binding have been shown to evolve rapidly in several groups of invertebrates and vertebrates. Mammalian SED1 (secreted protein containing N-terminal Notch-like type II epidermal growth factor (EGF) repeats and C-terminal discoidin/F5/8 C domains) is a recently identified sperm surface protein that binds the egg zona pellucida and facilitates sperm-egg adhesion. SED1-null male mice are subfertile. Here we examine the SED1 gene from 11 mammalian species and provide evidence that it underwent accelerated evolution in ancestral primates, most likely driven by positive selection. Specifically, the intensity of the positive selection across various protein domains of SED1 was heterogeneous. Although one of the 2 Notch-like EGF domains, which mediate protein-protein binding, was lost in primate SED1, the second EGF domain evolved under strong positive selection favoring polar to nonpolar amino acid replacements. By contrast, the 2 discoidin/F5/8 type C domains, which are involved in protein-cell membrane binding, do not show definite signs of positive selection. The structural modification and occurrence of directional selection in ancestral primates but not any other lineage suggest that the function of SED1 may have changed during primate evolution. These results reveal a different evolutionary pattern of SED1 from that of many other sperm-egg-binding proteins, which often show diversifying selection occurring in multiple lineages.     

Identification of low density lipoprotein receptor binding domains of human apolipoprotein B-100: a proposed consensus LDL receptor binding sequence of apoB-100

The human liver apoB-100 gene cloned in the lambda gt-11 expression vector expresses fusion proteins reacting with apoB antibodies. A fusion protein induced from a apoB-lambda gt-11 clone reacted with apoB-100 monoclonal antibodies known to block the binding of LDL to the LDL receptor. The fusion protein contains an amino acid sequence domain enriched in positively charged residues which is complementary to the negatively charged amino acids present in the consensus LDL receptor binding domain. This sequence of apoB-100 is proposed as a binding domain for the interaction with the LDL receptor. Comparison of derived amino acid sequences from the entire structure of apoB-100 molecule revealed several similar domains enriched in positively charged amino acids. A consensus sequence of the potential LDL binding domain was identified which contained positively charged amino acids at positions 1, 5 and 8 and a loop of 8-11 amino acids followed by two adjacent positively charged amino acids. These results are interpreted as indicating that there are several potential LDL receptor binding domains in apoB-100.     

Modulating the proteinase inhibitory profile of a plant cystatin by single mutations at positively selected amino acid sites

Cysteine proteinase inhibitors of the cystatin superfamily have several important functions in plants, including the inhibition of exogenous cysteine proteinases during herbivory or infection. Here we used a maximum-likelihood approach to assess whether plant cystatins, like other proteins implicated in host-pest interactions, have been subject to positive selection during the course of their evolution. Several amino acid sites were identified as being positively selected in cystatins from either Poaceae (monocots) and Solanaceae (dicots). These hypervariable sites were located at strategic positions on the protein: on each side of the conserved glycine residues in the N-terminal trunk, within the first and second inhibitory loops entering the active site of target enzymes, and surrounding the larfav motif, a sequence of unknown function conserved among plant cystatins. Supporting the assumption that positively selected, hypervariable sites are indicative of amino acid sites implicated in functional diversity, mutants of the 8th cystatin unit of tomato multicystatin including alternative residues at positively selected sites in the N-terminal trunk exhibited highly variable affinities for the cysteine proteases papain, cathepsin B and cathepsin H. Overall, these observations support the hypothesis that plant cystatins have been under selective pressure to evolve in response to predatory challenges by herbivorous enemies. They also indicate the potential of site-directed mutagenesis at positively selected sites for the generation of cystatins with improved binding properties.     

Evolution of Chitin-Binding Proteins in Invertebrates

Analysis of a group of invertebrate proteins, including chitinases and peritrophic matrix proteins, reveals the presence of chitin-binding domains that share significant amino acid sequence similarity. The data suggest that these domains evolved from a common ancestor which may be a protein containing a single chitin-binding domain. The duplication and transposition of this chitin-binding domain may have contributed to the functional diversification of chitin-binding proteins. Sequence comparisons indicated that invertebrate and plant chitin binding domains do not share significant amino acid sequence similarity, suggesting that they are not coancestral. However, both the invertebrate and the plant chitin-binding domains are cysteine-rich and have several highly conserved aromatic residues. In plants, cysteines have been elucidated in maintaining protein folding and aromatic amino acids in interacting with saccharides [Wright HT, Sanddrasegaram G, Wright CS (1991) J Mol Evol 33:283–294]. It is likely that these residues perform similar functions in invertebrates. We propose that the invertebrate and the plant chitin-binding domains share similar mechanisms for folding and saccharide binding and that they evolved by convergent evolution. Furthermore, we propose that the disulfide bonds and aromatic residues are hallmarks for saccharide-binding proteins. Received: 2 March 1998 / Accepted: 17 July 1998     

Positive selection for gains of N-linked glycosylation sites in hemagglutinin during evolution of H3N2 human influenza A virus

The number of N-linked glycosylation sites in the globular head of hemagglutinin (HA) has increased during evolution of H3N2 human influenza A virus. Here natural selection operating on the gains of N-linked glycosylation sites was examined by using the single-site analysis and the single-substitution analysis. In the single-site analysis, positive selection was not inferred at the amino acid sites where the substitutions generating N-linked glycosylation sites were observed, but was detected at antigenic sites. In contrast, in the single-substitution analysis, positive selection was detected for the amino acid substitutions generating N-linked glycosylation sites. The single-site analysis and the single-substitution analysis appeared to be suitable for detecting recurrent and episodic natural selection, respectively. The gains of N-linked glycosylation sites were likely to be positively selected for the function of shielding antigenic sites from immune responses. At the antigenic sites, positive selection appeared to have operated not only on the radical substitution but also on the conservative substitution in terms of the charge of amino acids, suggesting that the antigenic drift is not a by-product of the evolution of receptor binding avidity in HA of human H3N2 virus.     

Properties of the DNA-binding domain of the simian virus 40 large T antigen.

T antigen (Tag) from simian virus 40 binds specifically to two distinct sites in the viral origin of replication and to single-stranded DNA. Analysis of the protein domain responsible for these activities revealed the following. (i) The C-terminal boundary of the origin-specific and single-strand-specific DNA-binding domain is at or near amino acid 246; furthermore, the maximum of these DNA-binding activities coincides with a narrow C-terminal boundary, spanning 4 amino acids (246 to 249) and declines sharply in proteins with C termini which differ by a few (4 to 10) amino acids; (ii) a polypeptide spanning residues 132 to 246 of Tag is an independent domain responsible for origin-specific DNA binding and presumably for single-stranded DNA binding; and (iii) a comparison of identical N-terminal fragments of Tag purified from mammalian and bacterial cells revealed differential specificity and levels of activity between the two sources of protein. A role for posttranslational modification (phosphorylation) in controlling the DNA-binding activity of Tag is discussed.     

Heterogeneous processing and zona pellucida binding activity of pig zonadhesin

Zonadhesin is a mosaic protein in sperm membrane fractions that binds directly and in a species-specific manner to the extracellular matrix (zona pellucida) of the oocyte. The active form of pig zonadhesin from capacitated, epididymal spermatozoa comprises two covalently associated polypeptide chains of M(r) 105,000 (p105) and M(r) 45,000 (p45). Here we report detection and characterization of multiple zonadhesin isoforms in freshly ejaculated cells. Antibodies to the predicted von Willebrand D0-D1, D1, and D3 domains of pig zonadhesin recognized p105, p45, and additional M(r) 60,000-90,000 polypeptides in particulate fractions of uncapacitated cells. Although the p105/45 form constituted a minority of all zonadhesin forms in sperm membrane fractions, it was the predominant form capable of binding to the pig zona pellucida. Zonadhesin-binding sites were distributed over the entire zona pellucida. Anion exchange chromatography resolved active, p105/45 zonadhesin from the p60-90 inactive forms. Without disulfide bond reduction some zonadhesin was M(r) > or = 300,000, including M(r) 300,000 and 900,000 proteins comprising in part multimers of p105/45. The multimeric forms did not bind the zona pellucida as avidly as did the p105/45 monomer. Expressed D1 and D3 domain fragments containing the CG(L/V)CG sequence motif spontaneously formed multimers at -246 mV E(h) in vitro. Double Cys --> Ser mutants of the D1 fragment formed multimers with the same apparent kinetics as the wild type protein. Zonadhesin localized to the apical head of pig spermatozoa. We conclude that a heterogeneous combination of specific proteolysis and intermolecular disulfide bond formation in the sperm head generates multiple forms of zonadhesin with differing avidities for the zona pellucida.     

An electrostatic switch controls palmitoylation of the large conductance voltage- and calcium-activated potassium (BK) channel

Protein palmitoylation is a major dynamic posttranslational regulator of protein function. However, mechanisms that control palmitoylation are poorly understood. In many proteins, palmitoylation occurs at cysteine residues juxtaposed to membrane-anchoring domains such as transmembrane helices, sites of irreversible lipid modification, or hydrophobic and/or polybasic domains. In particular, polybasic domains represent an attractive mechanism to dynamically control protein palmitoylation, as the function of these domains can be dramatically influenced by protein phosphorylation. Here we demonstrate that a polybasic domain immediately upstream of palmitoylated cysteine residues within an alternatively spliced insert in the C terminus of the large conductance calcium- and voltage-activated potassium channel is an important determinant of channel palmitoylation and function. Mutation of basic amino acids to acidic residues within the polybasic domain results in inhibition of channel palmitoylation and a significant right-shift in channel half maximal voltage for activation. Importantly, protein kinase A-dependent phosphorylation of a single serine residue within the core of the polybasic domain, which results in channel inhibition, also reduces channel palmitoylation. These data demonstrate the key role of the polybasic domain in controlling stress-regulated exon palmitoylation and suggests that phosphorylation controls the domain by acting as an electrostatic switch.     

Lineage-specific differences in evolutionary mode in a salamander courtship pheromone

Functionally equivalent genes may evolve heterogeneously across closely related taxa as a consequence of lineage-specific selective pressures. Such disparate evolutionary modes are especially prevalent in genes that encode postcopulatory reproductive proteins, presumably as a result of sexual selection. We might therefore expect genes that mediate reproduction prior to insemination to evolve in a similar manner. Plethodontid receptivity factor (PRF), a proteinaceous salamander pheromone produced by the male, increases female receptivity during courtship interactions. To test for lineage-specific differences in PRF's evolution, we intensively sampled PRF genes across the eastern Plethodon phylogeny (27 spp.; 34 populations) to compare gene diversification, rates of evolution, modes of selection, and types of amino acid substitution. Our analyses indicate that PRF evolutionary dynamics vary considerably from lineage to lineage. Underlying this heterogeneity, however, are two well-defined transitions in evolutionary mode. The first mode is representative of a typical protein profile, wherein neutral divergence and purifying selection are the dominant features. The second mode is characterized by incessant, cyclical evolution driven by positive selection. In this mode, the positively selected sites are bound by a limited assortment of acceptable amino acids that appear to evolve independently of other sites, resulting in a tremendous number of unique PRF alleles. Several of these selected sites are implicated in receptor binding. These sites are apparently involved in a molecular tango in which the male signal and female receptors coevolve within a confined molecular space. PRF's lineage-specific evolutionary dynamics, in combination with evidence of a molecular tango, highlight the molecular action of sexual selection on a chemical signal that is used during courtship.     

Aminoacyl-tRNA synthetase classes and groups in prokaryotes

Knowledge on the evolution of aminoacyl-tRNA synthetases is crucial to studies on the origins of life. The relationships between the different aminoacyl-tRNA synthetase specificities in prokaryotic organisms are studied in this work. We reconstructed the ancestor sequences and the phylogenetic relationships utilizing the Maximum Likelihood method. The results suggest that in class I the evolution of the N-terminal segment was strongly influenced by the amino acid hydropathy in both domains of prokaryotes. The results for the C-terminal segments of class I were different in the two domains, indicating that its evolution was strongly influenced by the specific types of tRNA modification in each domain. The class II groups in Archaea were more heterogeneous with respect to the hydropathy of amino acids, indicating the interference of other influences. In bacteria, the configuration was also complex but the overall consensual division in two groups was maintained, group IIa forming a single branch with the five hydroapathetic amino acid specificities and group IIb containing the specificities for the moderately hydrophobic together with the hydrophilic amino acids. It is indicated that the aminoacyl-tRNA synthetase in both domains were subjected to different selective forces in diverse parts of the proteins, resulting in complex phylogenetic patterns.     

Peptide binding domains determined through chemical modification of the side-chain functional groups.

A clear understanding of the specific secondary structure and binding domain resulting from the interactions of proteins and peptides with lipid surfaces will provide insight into the specific functions of biologically active molecules. We have shown in earlier studies that the stationary phases used in reverse-phase high-performance liquid chromatography represent a model artificial lipid surface for the study of induced conformational states of peptides on lipid interaction. We have now used reverse-phase high-performance liquid chromatography to determine the binding domains of peptides and, by extension, of proteins to a lipid surface. This approach consists of performing chemical modifications of specific amino acid side-chain functionalities after the interaction of the peptides with the reverse-phase high-performance liquid chromatography C18 groups. The susceptibility to oxidation was also studied after binding of the same peptides to liposomes. Oxidation of a single methionine residue "walked" through an amphipathic alpha-helical 18-mer peptide was selected to illustrate this approach. The extent of oxidation was found to be clearly dictated by the accessibility of the methionine residue to the aqueous mobile phase. The binding domain found for the peptide in its lipid-induced conformational state was unequivocally the entire hydrophobic face of the amphipathic alpha-helix.