Evidence of positive selection at codon sites localized in the C-terminal peptide of ORC6

Origin recognition complex 6 (Orc6) plays a central role in the initiation of DNA replication in all eukaryotic systems. The exact contribution of Orc6 to replication initiation has yet to be elucidated. Here, we analyzed the evolutionary dynamics of Orc6 in 15 vertebrates. Positive selection was detected in the region of exon 6 of the Orc6 gene. Site tests revealed a proportion of codon sites that displayed evidence of positive selection (ω > 1) within the coding sequences of the vertebrate Orc6 gene. Seven positively selected amino acid sites were identified and three were located in exon6. These results suggest that amino acid residues present in the middle region of the protein are more selectively constrained, whereas amino acid residues in the C-terminal peptide of the protein evolve at a faster rate, possibly because of heightened selective pressure during the course of evolution.     

Bayesian analysis suggests that most amino acid replacements in Drosophila are driven by positive selection

One of the principal goals of population genetics is to understand the processes by which genetic variation within species (polymorphism) becomes converted into genetic differences between species (divergence). In this transformation, selective neutrality, near neutrality, and positive selection may each play a role, differing from one gene to the next. Synonymous nucleotide sites are often used as a uniform standard of comparison across genes on the grounds that synonymous sites are subject to relatively weak selective constraints and so may, to a first approximation, be regarded as neutral. Synonymous sites are also interdigitated with nonsynonymous sites and so are affected equally by genomic context and demographic factors. Hence a comparison of levels of polymorphism and divergence between synonymous sites and amino acid replacement sites in a gene is potentially informative about the magnitude of selective forces associated with amino acid replacements. We have analyzed 56 genes in which polymorphism data from D. simulans are compared with divergence from a reference strain of D. melanogaster. The framework of the analysis is Bayesian and assumes that the distribution of selective effects (Malthusian fitnesses) is Gaussian with a mean that differs for each gene. In such a model, the average scaled selection intensity (gamma = N(e)s) of amino acid replacements eligible to become polymorphic or fixed is -7.31, and the standard deviation of selective effects within each locus is 6.79 (assuming homoscedasticity across loci). For newly arising mutations of this type that occur in autosomal or X-linked genes, the average proportion of beneficial mutations is 19.7%. Among the amino acid polymorphisms in the sample, the expected average proportion of beneficial mutations is 47.7%, and among amino acid replacements that become fixed the average proportion of beneficial mutations is 94.3%. The average scaled selection intensity of fixed mutations is +5.1. The presence of positive selection is pervasive with the single exception of kl-5, a Y-linked fertility gene. We find no evidence that a significant fraction of fixed amino acid replacements is neutral or nearly neutral or that positive selection drives amino acid replacements at only a subset of the loci. These results are model dependent and we discuss possible modifications of the model that might allow more neutral and nearly neutral amino acid replacements to be fixed.     

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     

Evaluation of whether accelerated protein evolution in chordates has occurred before, after, or simultaneously with gene duplication

Gene duplication and loss are predicted to be at least of the order of the substitution rate and are key contributors to the development of novel gene function and overall genome evolution. Although it has been established that proteins evolve more rapidly after gene duplication, we were interested in testing to what extent this reflects causation or association. Therefore, we investigated the rate of evolution prior to gene duplication in chordates. Two patterns emerged; firstly, branches, which are both preceded by a duplication and followed by a duplication, display an elevated rate of amino acid replacement. This is reflected in the ratio of nonsynonymous to synonymous substitution (mean nonsynonymous to synonymous nucleotide substitution rate ratio [Ka:Ks]) of 0.44 compared with branches preceded by and followed by a speciation (mean Ka:Ks of 0.23). The observed patterns suggest that there can be simultaneous alteration in the selection pressures on both gene duplication and amino acid replacement, which may be consistent with co-occurring increases in positive selection, or alternatively with concurrent relaxation of purifying selection. The pattern is largely, but perhaps not completely, explained by the existence of certain families that have elevated rates of both gene duplication and amino acid replacement. Secondly, we observed accelerated amino acid replacement prior to duplication (mean Ka:Ks for postspeciation preduplication branches was 0.27). In some cases, this could reflect adaptive changes in protein function precipitating a gene duplication event. In conclusion, the circumstances surrounding the birth of new proteins may frequently involve a simultaneous change in selection pressures on both gene-copy number and amino acid replacement. More precise modeling of the relative importance of preduplication, postduplication, and simultaneous amino acid replacement will require larger and denser genomic data sets from multiple species, allowing simultaneous estimation of lineage-specific fluctuations in mutation rates and adaptive constraints.     

The phylogeny of human globin genes investigated by the maximum parsimony method

Gene phylogenetic trees were constructed by the maximum parsimony method for various sets of ninety six globin chain amino acid sequences spanning plant and animal kingdoms. The method, executed by several computer programs, constructed ancestor and descendant globin messengers on tree topologies which required the least number of nucleotide replacements to account for the evolution of the globins. The human myoglobin-hemoglobin divergence was traced to a gene duplication which occurred either in the first vertebrates or earlier yet in the common ancestor of chordates and annelids, the alpha-beta divergence to a gene duplication in the common ancestor of teleosts and tetrapods, the gamma divergence from typical beta chains to a gene duplication in basal therian mammals, and the delta separation from beta to a duplication in the basal catarrhine primates. Evidence was provided by the globin phylogenies for the hominoid affinities of the gibbon and the close phyletic relationship of the African apes to man. Over the period of teleos-tetrapod divergence the globin messengers evolved at an average rate of 18.5 nucleotide replacements per 100 codons per 10⁸ years, a faster rate than most previous estimates. Very fast and very slow rates were encountered in different globin lineages and at different stages of descent, reducing the effectiveness of globins as molecular clocks. Rates increased with gene duplication and decreased after selection discovered useful specializations in the products of genes which had previously been freer to accept mutations. The early eutherian radiation was characterized by rapid rates of globin evolution, but the later hominoid radiation by extremely slow rates. This pattern was related to more complicated grades of internal organization evolving in human ancestors. The types of nucleotide replacements in the globin messengers over the long course of globin evolution did not seem indicative of any special mutational mechanisms.     

Selection for functional uniformity of tuf duplicates in gamma-proteobacteria

Having an extra copy of a gene is thought to provide some functional redundancy, which results in a higher rate of evolution in duplicated genes. In this article, we estimate the impact of gene duplication on the selection of tuf paralogs, and we find that in the absence of gene conversion, tuf paralogs have evolved significantly slower than when gene conversion has been a factor in their evolution. Thus, tuf gene copies evolve under a selective pressure that ensures their functional uniformity, and gene conversion reduces selection against amino acid substitutions that affect the function of the encoded protein, EF-Tu.     

The early adaptive evolution of calmodulin

Interaction between gene duplication and natural selection in molecular evolution was investigated utilizing a phylogenetic tree constructed by the parsimony procedure from amino acid sequences of 50 calmodulin- family protein members. The 50 sequences, belonging to seven protein lineages related by gene duplication (calmodulin itself, troponin-C, alkali and regulatory light chains of myosin, parvalbumin, intestinal calcium-binding protein, and glial S-100 phenylalanine-rich protein), came from a wide range of eukaryotic taxa and yielded a denser tree (more branch points within each lineage) than in earlier studies. Evidence obtained from the reconstructed pattern of base substitutions and deletions in these ancestral loci suggests that, during the early history of the family, selection acted as a transforming force on expressed genes among the duplicates to encode molecular sites with new or modified functions. In later stages of descent, however, selection was a conserving force that preserved the structures of many coadapted functional sites. Each branch of the family was found to have a unique average tempo of evolutionary change, apparently regulated through functional constraints. Proteins whose functions dictate multiple interaction with several other macromolecules evolved more slowly than those which display fewer protein-protein and protein-ion interactions, e.g., calmodulin and next troponin-C evolved at the slowest average rates, whereas parvalbumin evolved at the fastest. The history of all lineages, however, appears to be characterized by rapid rates of evolutionary change in earlier periods, followed by slower rates in more recent periods. A particularly sharp contrast between such fast and slow rates is found in the evolution of calmodulin, whose rate of change in earlier eukaryotes was manyfold faster than the average rate over the past 1 billion years. In fact, the amino acid replacements in the nascent calmodulin lineage occurred at residue positions that in extant metazoans are largely invariable, lending further support to the Darwinian hypothesis that natural selection is both a creative and a conserving force in molecular evolution.      

The evolution of lysozyme and alpha-lactalbumin

From the analysis of phylogenetic trees constructed from the amino acid sequences and metal-binding properties of various lysozymes c and alpha-lactalbumins, it was found that before the divergence of the lineages of birds and mammals, calcium-binding lysozyme diverged from non-calcium-binding lysozyme. alpha-Lactalbumin evolved from the calcium-binding lysozyme along the mammalian lineage after the divergence of birds and mammals. Rapid evolution took place, not in the process of acquisition of the activity of alpha-lactalbumin, but after the loss of lysozyme activity, due to the change in the distribution of selective pressure on each amino acid site. A general process for the change in function of a protein during evolution is suggested to be as follows: after duplication of the gene, one of their protein products acquires a new function, besides that already present; the old function is eventually lost.     

Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages

The nonsynonymous (amino acid-altering) to synonymous (silent) substitution rate ratio (omega = d(N)/d(S)) provides a measure of natural selection at the protein level, with omega = 1, >1, and <1, indicating neutral evolution, purifying selection, and positive selection, respectively. Previous studies that used this measure to detect positive selection have often taken an approach of pairwise comparison, estimating substitution rates by averaging over all sites in the protein. As most amino acids in a functional protein are under structural and functional constraints and adaptive evolution probably affects only a few sites at a few time points, this approach of averaging rates over sites and over time has little power. Previously, we developed codon-based substitution models that allow the omega ratio to vary either among lineages or among sites. In this paper we extend previous models to allow the omega ratio to vary both among sites and among lineages and implement the new models in the likelihood framework. These models may be useful for identifying positive selection along prespecified lineages that affects only a few sites in the protein. We apply those branch-site models as well as previous branch- and site-specific models to three data sets: the lysozyme genes from primates, the tumor suppressor BRCA1 genes from primates, and the phytochrome (PHY) gene family in angiosperms. Positive selection is detected in the lysozyme and BRCA genes by both the new and the old models. However, only the new models detected positive selection acting on lineages after gene duplication in the PHY gene family. Additional tests on several data sets suggest that the new models may be useful in detecting positive selection after gene duplication in gene family evolution.     

Electric charge divergence in proteins: insights into the evolution of their three-dimensional properties

Previous studies of protein evolution have identified important mutations in various proteins that affect a small number of residues, but dramatically alter protein function. However, the evolutionary process underlying the three-dimensional protein properties, which are determined by a much larger number of residues, remains unclear. Based on a comparative evolutionary analysis of teleost phosphoglucose isomerases (PGIs; EC 5.3.1.9), we previously demonstrated that the relatively weak selection on many amino acid sites has played an important role in the evolution of protein electric charge as a model of three-dimensional protein properties. To ascertain the generality of this finding, we sought further evidence of this type of protein evolution. For this purpose, we analyzed the vertebrate isoforms of fructose-1,6-bisphosphate aldolase (ALD; EC 4.1.2.13), for which electric charges are known to have diverged after gene duplication. The results showed that the divergence in electric charge between the ALD isoforms was also driven by weak selection on many amino acid sites, as in PGI, confirming the generality of earlier findings. To obtain further insights, ALD and PGI were compared to the proteins pancreatic ribonuclease (EC 3.1.27.5) and triose-phosphate isomerase (EC 5.3.1.1), for which electric charges likely evolved through a well-defined mode of molecular evolution; namely, strong selection on specific amino acid sites. Comparison of the number and composition of amino acids on the protein surface suggested that the absolute number of evolutionarily changeable amino acids in a protein affects the strength of selection pressure acting on individual amino acid sites.     

Intrasubtype Reassortments Cause Adaptive Amino Acid Replacements in H3N2 Influenza Genes

Reassortments and point mutations are two major contributors to diversity of Influenza A virus; however, the link between these two processes is unclear. It has been suggested that reassortments provoke a temporary increase in the rate of amino acid changes as the viral proteins adapt to new genetic environment, but this phenomenon has not been studied systematically. Here, we use a phylogenetic approach to infer the reassortment events between the 8 segments of influenza A H3N2 virus since its emergence in humans in 1968. We then study the amino acid replacements that occurred in genes encoded in each segment subsequent to reassortments. In five out of eight genes (NA, M1, HA, PB1 and NS1), the reassortment events led to a transient increase in the rate of amino acid replacements on the descendant phylogenetic branches. In NA and HA, the replacements following reassortments were enriched with parallel and/or reversing replacements; in contrast, the replacements at sites responsible for differences between antigenic clusters (in HA) and at sites under positive selection (in NA) were underrepresented among them. Post-reassortment adaptive walks contribute to adaptive evolution in Influenza A: in NA, an average reassortment event causes at least 2.1 amino acid replacements in a reassorted gene, with, on average, 0.43 amino acid replacements per evolving post-reassortment lineage; and at least ∼9% of all amino acid replacements are provoked by reassortments.     

Testing the neutral fixation of hetero-oligomerism in the archaeal chaperonin CCT

The evolutionary transition from homo-oligomerism to hetero-oligomerism in multimeric proteins and its contribution to function innovation and organism complexity remain to be investigated. Here, we undertake the challenge of contributing to this theoretical ground by investigating the hetero-oligomerism in the molecular chaperonin cytosolic chaperonin containing tailless complex polypeptide 1 (CCT) from archaea. CCT is amenable to this study because, in contrast to eukaryotic CCTs where sub-functionalization after gene duplication has been taken to completion, archaeal CCTs present no evidence for subunit functional specialization. Our analyses yield additional information to previous reports on archaeal CCT paralogy by identifying new duplication events. Analyses of selective constraints show that amino acid sites from 1 subunit have fixed slightly deleterious mutations at inter-subunit interfaces after gene duplication. These mutations have been followed by compensatory mutations in nearby regions of the same subunit and in the interface contact regions of its paralogous subunit. The strong selective constraints in these regions after speciation support the evolutionary entrapment of CCTs as hetero-oligomers. In addition, our results unveil different evolutionary dynamics depending on the degree of CCT hetero-oligomerism. Archaeal CCT protein complexes comprising 3 distinct classes of subunits present 2 evolutionary processes. First, slightly deleterious and compensatory mutations were fixed neutrally at inter-subunit regions. Second, sub-functionalization may have occurred at substrate-binding and adenosine triphosphate-binding regions after the 2nd gene duplication event took place. CCTs with 2 distinct types of subunits did not present evidence of sub-functionalization. Our results provide the 1st in silico evidence for the neutral fixation of hetero-oligomerism in archaeal CCTs and provide information on the evolution of hetero-oligomerism toward sub-functionalization in archaeal CCTs.     

Compensatory gene amplification restores fitness after inter‐species gene replacements

Genes introduced by gene replacements and other types of horizontal gene transfer (HGT) represent a significant presence in many archaeal and eubacterial genomes. Most alien genes are likely to be neutral or deleterious upon arrival and their long‐term persistence may require a mechanism that improves their selective contribution. To examine the fate of inter‐species gene replacements, we exchanged three native S. typhimurium genes encoding ribosomal proteins with orthologues from various other microbes. The results show that replacement of each of these three genes reduces fitness to such an extent that it would provide an effective barrier against inter‐species gene replacements in eubacterial populations. However, these fitness defects could be partially ameliorated by gene amplification that augmented the dosage of the heterologous proteins. This suggests that suboptimal expression is a common fitness constraint for inter‐species gene replacements, with fitness costs conferred by either a lower expression level of the alien protein compared with the native protein or a requirement for an increased amount of the alien protein to maintain proper function. Our findings can explain the observation that duplicated genes are over‐represented among horizontally transferred genes, and suggest a potential coupling between compensatory gene amplification after HGT and the evolution of new genes.     

Toxin-resistant sodium channels: parallel adaptive evolution across a complete gene family

Approximately 75% of vertebrate proteins belong to protein families encoded by multiple evolutionarily related genes, a pattern that emerged as a result of gene and genome duplications over the course of vertebrate evolution. In families of genes with similar or related functions, adaptation to a strong selective agent should involve multiple adaptive changes across the entire gene family. However, we know of no evolutionary studies that have explicitly addressed this point. Here, we show how 4 taxonomically diverse species of pufferfishes (Tetraodontidae) each evolved resistance to the guanidinium toxins tetrodotoxin (TTX) and saxitoxin (STX) via parallel amino acid replacements across all 8 sodium channels present in teleost fish genomes. This resulted in diverse suites of coexisting sodium channel types that all confer varying degrees of toxin resistance, yet show remarkable convergence among genes and phylogenetically diverse species. Using site-directed mutagenesis and expression of a vertebrate sodium channel, we also demonstrate that resistance to TTX/STX is enhanced up to 15-fold by single, frequently observed replacements at 2 sites that have not previously been implicated in toxin binding but show similar or identical replacements in pufferfishes and in distantly related vertebrate and nonvertebrate animals. This study presents an example of natural selection acting upon a complete gene family, repeatedly arriving at a diverse but limited number of adaptive changes within the same genome. To be maximally informative, we suggest that future studies of molecular adaptation should consider all functionally similar paralogs of the affected gene family.     

Parallel evolution in the major haemoglobin genes of eight species of Andean waterfowl

Theory predicts that parallel evolution should be common when the number of beneficial mutations is limited by selective constraints on protein structure. However, confirmation is scarce in natural populations. Here we studied the major haemoglobin genes of eight Andean duck lineages and compared them to 115 other waterfowl species, including the bar-headed goose ( Anser indicus ) and Abyssinian blue-winged goose ( Cyanochen cyanopterus ), two additional species living at high altitude. One to five amino acid replacements were significantly overrepresented or derived in each highland population, and parallel substitutions were more common than in simulated sequences evolved under a neutral model. Two substitutions evolved in parallel in the αA subunit of two (Ala-α8) and five (Thr-α77) taxa, and five identical βA subunit substitutions were observed in two (Ser-β4, Glu-β94, Met-β133) or three (Ser-β13, Ser-β116) taxa. Substitutions at adjacent sites within the same functional protein region were also observed. Five such replacements were in exterior, solvent-accessible positions on the A helix and AB corner of the αA subunit. Five others were in close proximity to inositolpentaphosphate binding sites, and two pairs of independent replacements occurred at two different α¹β¹ intersubunit contacts. More than half of the substitutions in highland lineages resulted in the acquisition of serine or threonine (18 gains vs. 2 losses), both of which possess a hydroxyl group that can hydrogen bond to a variety of polar substrates. The patterns of parallel evolution observed in these waterfowl suggest that adaptation to high-altitude hypoxia has resulted from selection on unique but overlapping sets of one to five amino acid substitutions in each lineage.     

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.     

Evolved neomycin phosphotransferase from an isolate of Klebsiella pneumoniae

A new aminoglycoside resistance gene (aphA1-IAB) confers high-level resistance to neomycin. The sequence of aphA1-IAB is closely related to aphA1 found in the transposons Tn4352, Tn903 and Tn602. For example, aphA1-IAB differs from aphA1-903 at five nucleotides that result in four amino acid replacements. The enzyme encoded by aphA1-IAB has a significantly higher turnover number with neomycin, kanamycin and G418 as substrates than does the aphA1-903 enzyme. A parsimonious phylogenetic tree suggests that aphA1-IAB evolved from an ancestral form that is closely related or identical to the aphA1 found in Tn903. The excess of replacement substitutions over silent substitutions in aphA1-IAB, as well as its convergence toward aphA3 from Staphylococcus aureus, is indicative of selective evolution. Our hypothesis to explain these results is that aphA1-IAB evolved under the selective pressure of neomycin use in relatively recent times.     

Evidence for the adaptive evolution of the carbon fixation gene rbcL during diversification in temperature tolerance of a clade of hot spring cyanobacteria

Determining the molecular basis of enzyme adaptation is central to understanding the evolution of environmental tolerance but is complicated by the fact that not all amino acid differences between ecologically divergent taxa are adaptive. Analysing patterns of nucleotide sequence evolution can potentially guide the investigation of protein adaptation by identifying candidate codon sites on which diversifying selection has been operating. Here, I test whether there is evidence for molecular adaptation of the carbon fixation gene rbcL for a clade of hot spring cyanobacteria in the genus Synechococcus that has diverged in thermotolerance. Amino acid replacements during Synechococcus radiation have resulted in an increase in the number of hydrophobic residues in the RbcLs of more thermotolerant strains. A similar increase in hydrophobicity has been observed for many thermostable proteins. Maximum likelihood models which allow for heterogeneity among codon sites in the ratio of nonsynonymous to synonymous nucleotide substitutions estimated a class of amino acid sites as a target of positive selection. Depending on the model, a single amino acid site that interacts with a flexible element involved in the opening and closing of the active site was estimated with either low or moderate support to be a member of this class. Site-directed mutagenesis approaches are being explored in order to directly test its adaptive significance.