Evolution of protein inhibitors of serine proteinases: Positive Darwinian selection or compositional effects?

Summary In at least two instances involving serine proteinase inhibitors it has been shown that functionally important sites evolve faster and exhibit more interspecific variability than functionally neutral sites. Because these phenomena are difficult to reconcile with the neutral theory of molecular evolution, it has been suggested that the accelerated rate of amino acid substitution at the reactive sites is brought about by positive Darwinian selection. We show that differences in the amino acid composition in the different regions of proteinase inhibitors can account for the differences in the rates of amino acid substitution. By using an index of protein mutability [D. Graur (1985) J Mol Evol 2253–62], we show that the amino acid composition of the reactive center in the ovomucoids andSpi-2 gene products is such that, regardless of function, they are expected to evolve more rapidly than any other polypeptide for which the rate of substitution is known. In addition, the reactive region in theSpi-2 proteins is shown to be free of compositional constraint. Positive Darwinian selection need not be invoked at the present time in these cases.     

Testing for Spatial Clustering of Amino Acid Replacements Within Protein Tertiary Structure

Widely used models of protein evolution ignore protein structure. Therefore, these models do not predict spatial clustering of amino acid replacements with respect to tertiary structure. One formal and biologically implausible possibility is that there is no tendency for amino acid replacements to be spatially clustered during evolution. An alternative to this is that amino acid replacements are spatially clustered and this spatial clustering can be fully explained by a tendency for similar rates of amino acid replacement at sites that are nearby in protein tertiary structure. A third possibility is that the amount of clustering exceeds that which can be explained solely on the basis of independently evolving protein sites with spatially clustered replacement rates. We introduce two simple and not very parametric hypothesis tests that help distinguish these three possibilities. We then apply these tests to 273 homologous protein families. The null hypothesis of no spatial clustering is rejected for 102 of 273 families. The explanation of spatially clustered rates but independent change among sites is rejected for 43 families. These findings need to be reconciled with the common practice of basing evolutionary inferences on models that assume independent change among sites. [Reviewing Editior: Dr. David Pollock]     

Likelihood Analysis of the Chalcone Synthase Genes Suggests the Role of Positive Selection in Morning Glories (<Emphasis Type="Italic">Ipomoea</Emphasis>)

Chalcone synthase (CHS) is a key enzyme in the biosynthesis of flavonoides, which are important for the pigmentation of flowers and act as attractants to pollinators. Genes encoding CHS constitute a multigene family in which the copy number varies among plant species and functional divergence appears to have occurred repeatedly. In morning glories (Ipomoea), five functional CHS genes (A–E) have been described. Phylogenetic analysis of the Ipomoea CHS gene family revealed that CHS A, B, and C experienced accelerated rates of amino acid substitution relative to CHS D and E. To examine whether the CHS genes of the morning glories underwent adaptive evolution, maximum-likelihood models of codon substitution were used to analyze the functional sequences in the Ipomoea CHS gene family. These models used the nonsynonymous/synonymous rate ratio ( = d_N/d_S) as an indicator of selective pressure and allowed the ratio to vary among lineages or sites. Likelihood ratio test suggested significant variation in selection pressure among amino acid sites, with a small proportion of them detected to be under positive selection along the branches ancestral to CHS A, B, and C. Positive Darwinian selection appears to have promoted the divergence of subfamily ABC and subfamily DE and is at least partially responsible for a rate increase following gene duplication.     

Molecular Evolution of a Small Gene Family of Wound Inducible Kunitz Trypsin Inhibitors in Populus

Maximum likelihood models of codon substitutions were used to analyze the molecular evolution of a Kunitz trypsin inhibitor (KTI) gene family in Populus and Salix. The methods support previous assertions that the KTI genes comprise a rapidly evolving gene family. Models that allow for codon specific estimates of the ratio of nonsynonymous to synonymous substitutions (ω) among sites detect positive Darwinian selection at several sites in the KTI protein. In addition, branch-specific maximum likelihood models show that there is significant heterogeneity in ω among branches of the KTI phylogeny. In particular, ω is substantially higher following duplication than speciation. There is also evidence for significant rate heterogeneity following gene duplication, suggesting different evolutionary rates in newly arisen gene duplicates. The results indicate uneven evolutionary rates both between sites in the KTI protein and among different lineages in the KTI phylogeny, which is incompatible with a neutral model of sequence evolution. [Reviewing Editor: Dr. Willie J. Swanson]     

Proteome Composition in Plasmodium falciparum: Higher Usage of GC-Rich Nonsynonymous Codons in Highly Expressed Genes

The parasite Plasmodium falciparum, responsible for the most deadly form of human malaria, is one of the extremely AT-rich genomes sequenced so far and known to possess many atypical characteristics. Using multivariate statistical approaches, the present study analyzes the amino acid usage pattern in 5038 annotated protein-coding sequences in P. falciparum clone 3D7. The amino acid composition of individual proteins, though dominated by the directional mutational pressure, exhibits wide variation across the proteome. The Asn content, expression level, mean molecular weight, hydropathy, and aromaticity are found to be the major sources of variation in amino acid usage. At all stages of development, frequencies of residues encoded by GC-rich codons such as Gly, Ala, Arg, and Pro increase significantly in the products of the highly expressed genes. Investigation of nucleotide substitution patterns in P. falciparum and other Plasmodium species reveals that the nonsynonymous sites of highly expressed genes are more conserved than those of the lowly expressed ones, though for synonymous sites, the reverse is true. The highly expressed genes are, therefore, expected to be closer to their putative ancestral state in amino acid composition, and a plausible reason for their sequences being GC-rich at nonsynonymous codon positions could be that their ancestral state was less AT-biased. Negative correlation of the expression level of proteins with respective molecular weights supports the notion that P. falciparum, in spite of its intracellular parasitic lifestyle, follows the principle of cost minimization. [Reviewing Editor : Dr. Richard Kliman]     

Parallel amino acid replacements in the rhodopsins of the rockfishes (Sebastes spp.) associated with shifts in habitat depth

Among various groups of fishes, a shift in peak wavelength sensitivity has been correlated with changes in their photic environments. The genus Sebastes is a radiation of marine fish species that inhabit a wide range of depths from intertidal to over 600 m. We examined 32 species of Sebastes for evidence of adaptive amino acid substitution at the rhodopsin gene. Fourteen amino acid positions were variable among these species. Maximum likelihood analyses identify several of these to be targets of positive selection. None of these correspond to previously identified critical amino acid sites, yet they may in fact be functionally important. The occurrence of independent parallel changes at certain amino acid positions reinforces this idea. Reconstruction of habitat depths of ancestral nodes in the phylogeny suggests that shallow habitats have been colonized independently in different lineages. The evolution of rhodopsin appears to be associated with changes in depth, with accelerated evolution in lineages that have had large changes in depth.     

Structural and functional differences of Litopenaeus vannamei crustins

Penaeid crustins were described in Litopenaeus vannamei and L. setiferus as proteins belonging to an antibacterial peptide family with similar sequences but different sizes. Six crustin-coding clones were isolated from a cDNA library from L. vannamei hemocytes, sequenced and compared. Two different isoforms (named I and P) were found, based on two nucleotide differences that produce one change in amino acid sequence (Ile/Pro). Other single differences in nucleotide sequences were also noted, but they did not change the translated product. The mRNA steady state levels of crustin I, but not of crustin P, were down regulated by Vibrio alginolyticus inoculation. Thus, the differences among penaeid crustins seem to be associated with one amino acid substitution, which affects their expression after bacterial inoculation. By structural similarity, shrimp crustins seem to belong to an antibacterial WAP-domain containing protein family.     

Amino acid composition and the evolutionary rates of protein-coding genes

Summary Based on the rates of amino acid substitution for 60 mammalian genes of 50 codons or more, it is shown that the rate of amino acid substitution of a protein is correlated with its amino acid composition. In particular, the content of glycine residues is negatively correlated with the rate of amino acid substitution, and this content alone explains about 38% of the total variation in amino acid substitution rates among different protein families. The propensity of a polypeptide to evolve fast or slowly may be predicted from an index or indices of protein mutability directly derivable from the amino acid composition. The propensity of an amino acid to remain conserved during evolutionary times depends not so much on its being features prominently in active sites, but on its stability index, defined as the mean chemical distance [R. Grantham (1974) Science 185862–864] between the amino acid and its mutational derivatives produced by single-nucleotide substitutions. Functional constraints related to active and binding sites of proteins play only a minor role in determining the overall rate of amino acid substitution. The importance of amino acid composition in determining rates of substitution is illustrated with examples involving cytochrome c, cytochrome b₅,ras-related genes, the calmodulin protein family, and fibrinopeptides.     

Selection and Methionine Accumulation in the Fat Body Protein 2 Gene (FBP2), a Duplicate of the Drosophila Alcohol Dehydrogenase (ADH) Gene

The Drosophila fat body protein 2 gene (Fbp2) is an ancient duplication of the alcohol dehydrogenase gene (Adh) which encodes a protein that differs substantially from ADH in its methionine content. In D. melanogaster, there is one methionine in ADH, while there are 51 (20% of all amino acids) in FBP2. Methionine is involved in 46% of amino acid replacements when Fbp2 DNA sequences are compared between D. melanogaster and D. pseudoobscura. Methionine accumulation does not affect conserved residues of the ADH-ADH^r-FBP2 multigene family. The multigene family has evolved by replacement of mildly hydrophobic amino acids by methionine with no apparent reversion. Its short-term evolution was compared between two Drosophila species, while its long-term evolution was compared between two genera belonging respectively to acalyptrate and calyptrate Diptera, Drosophila and Sarcophaga. The pattern of nucleotide substitution was consistent with an independent accumulation of methionines at the Fbp2 locus in each lineage. Under a steady-state model, the rate of methionine accumulation was constant in the lineage leading to Drosophila, and was twice as fast as that in the calyptrate lineage. Substitution rates were consistent with a slight positive selective advantage for each methionine change in about one-half of amino acid sites in Drosophila. This shows that selection can potentially account for a large proportion of amino acid replacements in the molecular evolution of proteins. Received: 12 December 1994 / Accepted: 15 April 1996     

Asymmetric evolution of duplicate genes encoding the CCAAT-binding factor NF-Y in plant genomes

NF-Y is a ubiquitous CCAAT-binding factor composed of NF-YA, NF-YB and NF-YC. Multiple genes encoding NF-Y subunits have been identified in plant genomes. It remains unclear whether the duplicate genes underwent different evolutionary patterns. Likelihood-ratio tests were used to examine whether the amino acid substitution rates are the same between duplicate genes. The influences of selection on evolution were evaluated by comparing the conservative and radical amino acid substitution rates, as well as maximum-likelihood analysis. Some NF-YB and NF-YC duplicates showed significant evidence of asymmetric evolution but not the NF-YA duplicates. Most amino acid replacements in the NF-YB and NF-YC duplicates result in changes in hydropathy, polar requirement and polarity. The physicochemical changes in the sequences of NF-YB seem to be coupled to asymmetric divergence in gene function. Plant NF-Y genes have evolved in different patterns. Relaxed selective constraints following gene duplication are most likely responsible for the unequal evolutionary rates and distinct divergence patterns of duplicate NF-Y genes. Positive selection may have promoted amino acid hydropathy changes in the NF-YC duplicates.     

Comparative analysis of vertebrate <Emphasis Type="Italic">PEPT1</Emphasis> and <Emphasis Type="Italic">PEPT2</Emphasis> genes

The plasma membrane transport proteins belong to SoLute Carrier 15 (SLC15) family and two members of this family have been characterized extensively in higher vertebrates, namely PEPT1 and PEPT2. Despite many efforts have made to define a pharmacophore model for efficient binding and transporting of substrates, there is not a comprehensive study performed to elucidate the evolutionary mechanisms among the SLC15 family members and to statistically evaluate sequence conservation and functional divergence between members. In this study, we compared and contrasted the rates and patterns of molecular evolution of 2 PEPT genes. Phylogenetic tree assembly with all available vertebrate PEPTs suggests that the PEPTs originated by duplications and diverged from a common protein at the base of the eukaryotic tree. Topological structure demonstrates both members share the similar hydrophobic domains (TMDs), which have been constrained by purifying selection. Although both genes show qualitatively similar patterns, their rates of evolution differ significantly due to an increased rate of synonymous substitutions in the structural domains in one copy, suggesting substantial differences in functional constraint on each gene. Site-specific profiles were established by posterior probability analysis revealing significantly divergent regions mainly locate at the hydrophobic region between predicted transmembrane domains 9 and 10 of the proteins. Thus, these results provide the evidence that several amino acid residues with reduced selective constraints are largely responsible for functional divergence between the paralogous PEPTs. These findings may provide a starting point for further experimental verifications.     

Assessing Substitution Variation Across Sites in Grass Chloroplast DNA

We assess the similarity of base substitution processes, described by empirically derived 4 × 4 matrices, using chi-square homogeneity tests. Such significance analyses allow us to assess variation in sequence evolution across sites and we apply them to matrices derived from noncoding sites in different contexts in grass chloroplast DNA. We show that there is statistically significant variation in rates and patterns of mutation among noncoding sites in different contexts and then demonstrate a similar and significant influence of context on substitutions at fourfold degenerate sites of coding regions from grass chloroplast DNA. These results show that context has the same general effect on substitution bias in coding and noncoding DNA: the A+T content of flanking bases is correlated with rate of substitution, transition bias, and GC → AT pressure, while the number of flanking pyrimidines on a single strand is correlated with a mutational bias, or skew, toward pyrimidines. Despite the similarity in general trends, however, when we compare coding and noncoding matrices we find that there is a statistically significant difference between them even when we control for context. Most noticeably, fourfold degenerate sites in coding sequences are undergoing substitution at a higher rate and there are also significant differences in the relationship between pyrimidines skew and the number of flanking pyrimidines. Possible reasons for the differences between coding and noncoding sites are discussed. Furthermore, our analysis illustrates a simple statistical way for comparing substitution processes across sites allowing us to better study variation in evolutionary processes across a genome. [Reviewing Editor: Dr. Martin Kreitman]     

The Rate of Cu,Zn Superoxide Dismutase Evolution

The rate of amino acid replacement in Cu, Zn SOD greatly departs from the expectations of the molecular clock. We examine 27 Cu, Zn SOD sequences available and conclude that: (I) the SOD enzymes from different mammal families differ from each other by roughly the same number of replacements, which is consistent with a simultaneous mammalian radiation; (2) over the most recent 60 million years (MY) the rate of SOD evolution is fairly high (15aa/100aa/100MYR) and may be considered constant; (3) the rate of accumulation of amino acid replacements since the divergence of fungi. plants and animals to the present is inconstant along different branches of the evolutionary tree; moreover it steadily decreases with time, to the same extent in all lineages; (4) some comparisons exhibit divergences that are in any case greater than expected from a Poisson process on the assumption of a molecular clock; (5) plant chloroplast enzymes display fewer differences from each other than cytoplasmic ones; (6) bacteriocuprein (from Photobacterium leiognathi), fluke and human extracellular SOD are all three extremely remotely related to one another and to the SOD of other eukaryotes.

The process of consistent decline of the rate of evolution of Cu. Zn SOD can be described by a number of mathematical functions. We explore simple models that assume constant rates and might be applicable to other proteins or genes that apparently evolve at disparate rates.     

Effective population size and the rate and pattern of nucleotide substitutions

Both the overall rate of nucleotide substitution and the relative proportions of synonymous and non-synonymous substitutions are predicted to vary between species that differ in effective population size (N_e). Our understanding of the genetic processes underlying these lineage-specific differences in molecular evolution is still developing. Empirical analyses indicate that variation in substitution rates and patterns caused by differences in N_e is often substantial, however, and must be accounted for in analyses of molecular evolution.     

Molecular characterization and expression analysis of a gene encoding mannose-binding lectin from bulb of Zephyranthes grandiflora

In this paper we report on the molecular cloning, sequencing and partially characterisation of a lectin from bulb of the Chinese medicinal plant Zephyranthes grandiflora. The full-length cDNA of Z. grandiflora bulb lectin (ZGBL) consisted of 986 bp and contained a 576 bp ORF encoding a 191 amino acid protein. Bioinformatics analysis results clearly indicate that ZGBL belongs to the monocot mannose-binding lectin family, which contains 3 putative mannose-binding sites per subunit. RT-PCR analysis results indicate that ZGBL is constitutively expressed in all the tested tissue types including root, bulb, leaf and flower. Interestingly, ZGBL is more closely related to the Orchidaceae rather than the Amaryllidaceae family on molecular evolution.     

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.     

High‐performance variants of plant diacylglycerol acyltransferase 1 generated by directed evolution provide insights into structure function

Diacylglycerol acyltransferase 1 (DGAT1) catalyzes the acyl‐CoA‐dependent biosynthesis of triacylglycerol, the predominant component of seed oil. In some oil crops, including Brassica napus, the level of DGAT1 activity can have a substantial effect on triacylglycerol production. Structure–function insights into DGAT1, however, remain limited because of the lack of a three‐dimensional detailed structure for this membrane‐bound enzyme. In this study, the amino acid residues governing B. napus DGAT1 (BnaDGAT1) activity were investigated via directed evolution, targeted mutagenesis, in vitro enzymatic assay, topological analysis, and transient expression of cDNA encoding selected enzyme variants in Nicotiana benthamiana. Directed evolution revealed that numerous amino acid residues were associated with increased BnaDGAT1 activity, and 67% of these residues were conserved among plant DGAT1s. The identified amino acid residue substitution sites occur throughout the BnaDGAT1 polypeptide, with 89% of the substitutions located outside the putative substrate binding or active sites. In addition, cDNAs encoding variants I447F or L441P were transiently overexpressed in N. benthamiana leaves, resulting in 33.2 or 70.5% higher triacylglycerol content, respectively, compared with native BnaDGAT1. Overall, the results provide novel insights into amino acid residues underlying plant DGAT1 function and performance‐enhanced BnaDGAT1 variants for increasing vegetable oil production.     

The Relationship Between the Rate of Molecular Evolution and the Rate of Genome Rearrangement in Animal Mitochondrial Genomes

Evolution of mitochondrial genes is far from clock-like. The substitution rate varies considerably between species, and there are many species that have a significantly increased rate with respect to their close relatives. There is also considerable variation among species in the rate of gene order rearrangement. Using a set of 55 complete arthropod mitochondrial genomes, we estimate the evolutionary distance from the common ancestor to each species using protein sequences, tRNA sequences, and breakpoint distances (a measure of the degree of genome rearrangement). All these distance measures are correlated. We use relative rate tests to compare pairs of related species in several animal phyla. In the majority of cases, the species with the more highly rearranged genome also has a significantly higher rate of sequence evolution. Species with higher amino acid substitution rates in mitochondria also have more variable amino acid composition in response to mutation pressure. We discuss the possible causes of variation in rates of sequence evolution and gene rearrangement among species and the possible reasons for the observed correlation between the two rates. [Reviewing Editor: Dr. David Pollock]     

Estimation of DNA Sequence Context-dependent Mutation Rates Using Primate Genomic Sequences

It is understood that DNA and amino acid substitution rates are highly sequence context-dependent, e.g., C --> T substitutions in vertebrates may occur much more frequently at CpG sites and that cysteine substitution rates may depend on support of the context for participation in a disulfide bond. Furthermore, many applications rely on quantitative models of nucleotide or amino acid substitution, including phylogenetic inference and identification of amino acid sequence positions involved in functional specificity. We describe quantification of the context dependence of nucleotide substitution rates using baboon, chimpanzee, and human genomic sequence data generated by the NISC Comparative Sequencing Program. Relative mutation rates are reported for the 96 classes of mutations of the form 5' alphabetagamma 3' --> 5' alphadeltagamma 3', where alpha, beta, gamma, and delta are nucleotides and beta not equal delta, based on maximum likelihood calculations. Our results confirm that C --> T substitutions are enhanced at CpG sites compared with other transitions, relatively independent of the identity of the preceding nucleotide. While, as expected, transitions generally occur more frequently than transversions, we find that the most frequent transversions involve the C at CpG sites (CpG transversions) and that their rate is comparable to the rate of transitions at non-CpG sites. A four-class model of the rates of context-dependent evolution of primate DNA sequences, CpG transitions > non-CpG transitions approximately CpG transversions > non-CpG transversions, captures qualitative features of the mutation spectrum. We find that despite qualitative similarity of mutation rates among different genomic regions, there are statistically significant differences.     

The divergence and positive selection of the plant‐specific BURP‐containing protein family

BURP domain‐containing proteins belong to a plant‐specific protein family and have diverse roles in plant development and stress responses. However, our understanding about the genetic divergence patterns and evolutionary rates of these proteins remain inadequate. In this study, 15 plant genomes were explored to elucidate the genetic origins, divergence, and functions of these proteins. One hundred and twenty‐five BURP protein‐encoding genes were identified from four main plant lineages, including 13 higher plant species. The absence of BURP family genes in unicellular and multicellular algae suggests that this family (1) appeared when plants shifted from relatively stable aquatic environments to land, where conditions are more variable and stressful, and (2) is critical in the adaptation of plants to adverse environments. Promoter analysis revealed that several responsive elements to plant hormones and external environment stresses are concentrated in the promoter region of BURP protein‐encoding genes. This finding confirms that these genes influence plant stress responses. Several segmentally and tandem‐duplicated gene pairs were identified from eight plant species. Thus, in general, BURP domain‐containing genes have been subject to strong positive selection, even though these genes have conformed to different expansion models in different species. Our study also detected certain critical amino acid sites that may have contributed to functional divergence among groups or subgroups. Unexpectedly, all of the critical amino acid residues of functional divergence and positive selection were exclusively located in the C‐terminal region of the BURP domain. In conclusion, our results contribute novel insights into the genetic divergence patterns and evolutionary rates of BURP proteins.