Amino Acid Sequence Evolution Research Articles

Amphiphilic α-Helical Potential: A Putative Folding Motif Adding Few Constraints to Protein Evolution

Evidence from a number of studies indicates that protein folding is dictated not only by factors stabilizing the native state, but also by potentially independent factors that create folding pathways. How natural selection might cope simultaneously with two independent factors was addressed in this study within the framework of the "Lim-model" of protein folding, which postulates that the early stages of folding of all globular proteins, regardless of their native structure, are directed at least in part by potential to form amphiphilic α-helices. For this purpose, the amphiphilic α-helical potential in randomly ordered amino acid sequences and the conservation in phylogeny of amphiphilic α-helical potential within various proteins were assessed. These analyses revealed that amphiphilic α-helical potential is a common occurrence in random sequences, and that the presence of amphiphilic α-helical potential is present but not conserved in phylogeny within a given protein. The results suggest that the rapid formation of molten globules and the variable behavior of those globules depending on the protein may be a fundamental property of polymers of naturally occurring amino acids more so than a trait that must be derived or maintained by natural selection. Further, the results point toward the utility of randomly occurring process in protein function and evolution, and suggest that the formation of efficient pathways that determine early processes in protein folding, unlike the formation of stable, native protein structure, does not present a substantial hurdle during the evolution of amino acid sequences.

Dramatic Change in Function and Expression Pattern of a Gene Duplicated by Polyploidy Created a Paternal Effect Gene in the Brassicaceae

New gene formation by polyploidy has been an ongoing process during the evolution of various eukaryotes that has contributed greatly to the large number of genes in their genomes. After duplication, some genes that are retained can acquire new functions or expression patterns, or subdivide their functions or expression patterns between duplicates. Here, we show that SHORT SUSPENSOR (SSP) and Brassinosteroid Kinase 1 (BSK1) are paralogs duplicated by a polyploidy event that occurred in the Brassicaceae family about 23 Ma. SSP is involved in paternal control of zygote elongation in Arabidopsis thaliana by transcription in the sperm cells of pollen and then translation in the zygote, whereas BSK1 is involved in brassinosteroid signal transduction. Comparative analysis of expression in 63 different organs and developmental stages revealed that BSK1 and SSP have opposite expression patterns in pollen compared with all other parts of the plant. We determined that BSK1 retains the ancestral expression pattern and function. Thus, SSP has diverged in function after duplication from a component of the brassinosteroid signaling pathway to a paternal regulator of the timing of zygote elongation. The ancestral function of SSP was lost by deletions in the kinase domain. Our sequence rate analysis revealed that SSP but not BSK1 has experienced a greatly accelerated rate of amino acid sequence changes and relaxation of purifying selection. In addition, SSP has been duplicated to create a new gene (SSP-like1) with a completely different expression pattern, a shorter coding sequence that has lost a critical functional domain, and a greatly accelerated rate of amino acid sequence evolution along with evidence for positive selection, together indicative of neofunctionalization. This study illustrates two dramatic examples of neofunctionalization following gene duplication by complete changes in expression pattern and function. In addition, our findings indicate that paternal control of zygote elongation by SSP is an evolutionarily recent innovation in the Brassicaceae family.

Molecular systematics: A synthesis of the common methods and the state of knowledge

The comparative and evolutionary analysis of molecular data has allowed researchers to tackle biological questions that have long remained unresolved. The evolution of DNA and amino acid sequences can now be modeled accurately enough that the information conveyed can be used to reconstruct the past. The methods to infer phylogeny (the pattern of historical relationships among lineages of organisms and/or sequences) range from the simplest, based on parsimony, to more sophisticated and highly parametric ones based on likelihood and Bayesian approaches. In general, molecular systematics provides a powerful statistical framework for hypothesis testing and the estimation of evolutionary processes, including the estimation of divergence times among taxa. The field of molecular systematics has experienced a revolution in recent years, and, although there are still methodological problems and pitfalls, it has become an essential tool for the study of evolutionary patterns and processes at different levels of biological organization. This review aims to present a brief synthesis of the approaches and methodologies that are most widely used in the field of molecular systematics today, as well as indications of future trends and state-of-the-art approaches.

On Theoretical Models of Gene Expression Evolution with Random Genetic Drift and Natural Selection

BackgroundThe relative contributions of natural selection and random genetic drift are a major source of debate in the study of gene expression evolution, which is hypothesized to serve as a bridge from molecular to phenotypic evolution. It has been suggested that the conflict between views is caused by the lack of a definite model of the neutral hypothesis, which can describe the long-run behavior of evolutionary change in mRNA abundance. Therefore previous studies have used inadequate analogies with the neutral prediction of other phenomena, such as amino acid or nucleotide sequence evolution, as the null hypothesis of their statistical inference.Methodology/Principal FindingsIn this study, we introduced two novel theoretical models, one based on neutral drift and the other assuming natural selection, by focusing on a common property of the distribution of mRNA abundance among a variety of eukaryotic cells, which reflects the result of long-term evolution. Our results demonstrated that (1) our models can reproduce two independently found phenomena simultaneously: the time development of gene expression divergence and Zipf's law of the transcriptome; (2) cytological constraints can be explicitly formulated to describe long-term evolution; (3) the model assuming that natural selection optimized relative mRNA abundance was more consistent with previously published observations than the model of optimized absolute mRNA abundances.Conclusions/SignificanceThe models introduced in this study give a formulation of evolutionary change in the mRNA abundance of each gene as a stochastic process, on the basis of previously published observations. This model provides a foundation for interpreting observed data in studies of gene expression evolution, including identifying an adequate time scale for discriminating the effect of natural selection from that of random genetic drift of selectively neutral variations.

The molecular evolution of HIV‐1 protease simulated at atomic detail

Progress in understanding protein folding allows to simulate, with atomic detail, the evolution of amino-acid sequences folding to a given native conformation. A particularly attractive example is the HIV-1 protease, main target of therapies to fight AIDS, which under drug pressure is able to develop resistance within few months from the starting of therapy. By comparing the results of simulations of the evolution of the protease with the corresponding proteomic data, one can approximately determine the value of the associated evolution pressure under which the enzyme has become and, as a consequence, map out the energy landscape in sequence space of the HIV-1 protease. It is found that there are several families of sequences folding to the native conformations of the enzyme. Each of these families are characterized by different sets of highly conserved ("hot") amino acids which play a critical role in the folding and stability of the protease. There are two main possibilities for the virus to move from one family to a different one: (a) in a single generation, through the concerted mutations of the hot amino acids, a highly unlikely event, (b) through a folding path (if it exists), again a very improbable event. In fact, the number of generations needed by the virus to change stepwise its sequence from one family to another is astronomically large. These results point to the "hot" segments of the protease as promising targets for a nonconventional inhibition strategy, likely not to create resistance.

Reduced Effectiveness of Selection Caused by a Lack of Recombination

Genetic recombination associated with sexual reproduction is expected to have important consequences for the effectiveness of natural selection. These effects may be evident within genomes, in the form of contrasting patterns of molecular variation and evolution in regions with different levels of recombination. Previous work reveals patterns that are consistent with a benefit of recombination for adaptation at the level of protein sequence: both positive selection for adaptive variants and purifying selection against deleterious ones appear to be compromised in regions of low recombination [1-11]. Here, we re-examine these patterns by using polymorphism and divergence data from the Drosophila dot chromosome, which has a long history of reduced recombination. To avoid confounding selection and demographic effects, we collected these data from a species with an apparently stable demographic history, Drosophila americana. We find that D. americana dot loci show several signatures of ineffective purifying and positive selection, including an increase in the rate of protein evolution, an increase in protein polymorphism, and a reduction in the proportion of amino acid substitutions attributable to positive selection.

Trends in Prokaryotic Evolution Revealed by Comparison of Closely Related Bacterial and Archaeal Genomes

In order to explore microevolutionary trends in bacteria and archaea, we constructed a data set of 41 alignable tight genome clusters (ATGCs). We show that the ratio of the medians of nonsynonymous to synonymous substitution rates (dN/dS) that is used as a measure of the purifying selection pressure on protein sequences is a stable characteristic of the ATGCs. In agreement with previous findings, parasitic bacteria, notwithstanding the sometimes dramatic genome shrinkage caused by gene loss, are typically subjected to relatively weak purifying selection, presumably owing to relatively small effective population sizes and frequent bottlenecks. However, no evidence of genome streamlining caused by strong selective pressure was found in any of the ATGCs. On the contrary, a significant positive correlation between the genome size, as well as gene size, and selective pressure was observed, although a variety of free-living prokaryotes with very close selective pressures span nearly the entire range of genome sizes. In addition, we examined the connections between the sequence evolution rate and other genomic features. Although gene order changes much faster than protein sequences during the evolution of prokaryotes, a strong positive correlation was observed between the "rearrangement distance" and the amino acid distance, suggesting that at least some of the events leading to genome rearrangement are subjected to the same type of selective constraints as the evolution of amino acid sequences.

Modelling the evolution of protein coding sequences sampled from Measurably Evolving Populations.

Models of nucleotide or amino acid sequence evolution that implement homogeneous and stationary Markov processes of substitutions are mathematically convenient but are unlikely to represent the true complexity of evolution. With the large amounts of data that next generation sequencing promises, appropriate models of evolution are important, particularly when data are collected from ancient and sub-fossil remains, where changes in evolutionary parameters are the norm and not the exception. In this paper, we describe a new codon-based model of evolution that applies to Measurably Evolving Populations (MEPs). A MEP is defined as a population from which it is possible to detect a statistically significant accumulation of substitutions when sequences are obtained at different times. The new model of codon evolution permits changes to the substitution process, including changes to the intensity of selection and the proportions of sites undergoing different selective pressures. In our serial model of codon evolution, changes in the selective regime occur simultaneously across all lineages. Different regions of the protein may also evolve under distinct selective patterns. We illustrate the application of the new model to a dataset of HIV-1 sequences obtained from an infected individual before and after the commencement of antiretroviral therapy.

Mutation screening of the COCH gene in familial and sporadic patients with late onset nonsyndromic sensorineural hearing loss among Chinese population

To investigate the mutational of the coagulation factor C homology (COCH) gene related to autosomal dominant sensorineural nonsyndromic hearing loss (DFNA) with late onset in Chinese population. Peripheral blood samples were collected from he members of 26 DFNA families, members of 19 small DFNA families with un recognized inheritance pattern, and 22 sporadic patients with sensorineural nonsyndromic late onset hearing loss, the hearing loss of all of which occurred during the age range 10 - 40, and 100 normal controls. From different parts of China, these subjects underwent questionnaire survey too. Genomic DNA was isolated, COCH mutation was screened by PCR and sequencing, and restriction endonuclease analysis was used to detect the mutation sites of the COCH gene. The conservation in evolution of the target amino acid sequences was analyzed using CluatalX1.82 software. DNA sequencing of coding regions and exon/intron boundaries of COCH 2 - 12 exons identified a heterozygous G-to-A substitution at position 1625 in exon 12 in a large DFNA family, leading to a C542Y substitution, and a heterozygous T-to-C substitution at position 1535 in exon 12 in a small family, leading to a M512T substitutions. Both the residues of Cys542 and M512 were conserved across human, mouse, chicken, and zebrafish. These mutations were not detected in the 100 control subjects. The C542Y and the M512T mutations cause hearing loss in Chinese DFNA families.

The nematode-arthropod clade revisited: phylogenomic analyses from ribosomal protein genes misled by shared evolutionary biases.

Phylogenetic analysis of major groups of Metazoa using genomic data tends to recover the sister relationships of arthropods and chordates, contradicting the proposed Ecdysozoa (the molting animals), which group the arthropods together with nematodes and relatives. Ribosomal protein genes have been a major data source in phylogenomic studies because they are readily detected as Expressed Sequence Tags (ESTs) due to their high transcription rates. Here we address the debate about the recovery of Ecdysozoa in genomic data by building a new matrix of carefully curated EST and genome sequences for 25 ribosomal protein genes of the small subunit, with focus on new insect sequences in addition to the Diptera sequences generally used to represent the arthropods. Individually, each ribosomal protein gene showed low phylogenetic signal, but in simultaneous analysis strong support emerged for many expected groups, with support increasing linearly with increased gene number. In agreement with most studies of metazoan relationships from genomic data, our analyses contradicted the Ecdysozoa (the putative sister relationship of arthropods and nematodes), and instead supported the affinity of arthropods with chordates. In addition, relationships among holometabolan insects resulted in an unlikely basal position for Diptera. To test for biases in the data that might produce an erroneous arthropod-chordate affinity we simulated sequence data on tree topologies with the alternative arthropod-nematode sister relationships, applying a model of amino acid sequence evolution estimated from the real data. Tree searches on these simulated data still revealed an arthropod-chordate grouping, i.e., the topologies used to simulate the data were not recovered correctly. This suggests that the arthropod-chordate relationships may be obtained erroneously also from the real data even if the alternative topology (Ecdysozoa) represents the true phylogeny. Whereas denser taxon sampling in the future may recover the Ecdysozoa, our analyses demonstrate that recent phylogenomic studies may be affected by as yet unspecified biases in amino acid sequence composition in the model organisms with available genomic data.

Evolution of Amino-Acid Sequences and Codon Usage on the Drosophila miranda Neo-Sex Chromosomes

We have studied patterns of DNA sequence variation and evolution for 22 genes located on the neo-X and neo-Y chromosomes of Drosophila miranda. As found previously, nucleotide site diversity is greatly reduced on the neo-Y chromosome, with a severely distorted frequency spectrum. There is also an accelerated rate of amino-acid sequence evolution on the neo-Y chromosome. Comparisons of nonsynonymous and silent variation and divergence suggest that amino-acid sequences on the neo-X chromosome are subject to purifying selection, whereas this is much weaker on the neo-Y. The same applies to synonymous variants affecting codon usage. There is also an indication of a recent relaxation of selection on synonymous mutations for genes on other chromosomes. Genes that are weakly expressed on the neo-Y chromosome appear to have a faster rate of accumulation of both nonsynonymous and unpreferred synonymous mutations than genes with high levels of expression, although the rate of accumulation when both types of mutation are pooled is higher for the neo-Y chromosome than for the neo-X chromosome even for highly expressed genes.

Site interdependence attributed to tertiary structure in amino acid sequence evolution

Standard likelihood-based frameworks in phylogenetics consider the process of evolution of a sequence site by site. Assuming that sites evolve independently greatly simplifies the required calculations. However, this simplification is known to be incorrect in many cases. Here, a computational method that allows for general dependence between sites of a sequence is investigated. Using this method, measures acting as sequence fitness proxies can be considered over a phylogenetic tree. In this work, a set of statistically derived amino acid pairwise potentials, developed in the context of protein threading, is used to account for what we call the structural fitness of a sequence. We describe a model combining statistical potentials with an empirical amino acid substitution matrix. We propose such a combination as a useful way of capturing the complexity of protein evolution. Finally, we outline features of the model using three datasets and show the approach's sensitivity to different tree topologies.

A novel use of equilibrium frequencies in models of sequence evolution.

Current mathematical models of amino acid sequence evolution are often applied in variants that match their expected amino acid frequencies to those observed in a data set under analysis. This has been achieved by setting the instantaneous rate of replacement of a residue i by another residue j proportional to the observed frequency of the resulting residue j. We describe a more general method that maintains the match between expected and observed frequencies but permits replacement rates to be proportional to the frequencies of both the replaced and resulting residues, raised to powers other than 1. Analysis of a database of amino acid alignments shows that the description of the evolutionary process in a majority (approximately 70% of 182 alignments) is significantly improved by use of the new method, and a variety of analyses indicate that parameter estimation with the new method is well-behaved. Improved evolutionary models increase our understanding of the process of molecular evolution and are often expected to lead to improved phylogenetic inferences, and so it seems justified to consider our new variants of existing standard models when performing evolutionary analyses of amino acid sequences. Similar methods can be used with nucleotide substitution models, but we have not found these to give corresponding significant improvements to our ability to describe the processes of nucleotide sequence evolution.

Constant relative rate of protein evolution and detection of functional diversification among bacterial, archaeal and eukaryotic proteins.

Detection of changes in a protein's evolutionary rate may reveal cases of change in that protein's function. We developed and implemented a simple relative rates test in an attempt to assess the rate constancy of protein evolution and to detect cases of functional diversification between orthologous proteins. The test was performed on clusters of orthologous protein sequences from complete bacterial genomes (Chlamydia trachomatis, C. muridarum and Chlamydophila pneumoniae), complete archaeal genomes (Pyrococcus horikoshii, P. abyssi and P. furiosus) and partially sequenced mammalian genomes (human, mouse and rat). Amino-acid sequence evolution rates are significantly correlated on different branches of phylogenetic trees representing the great majority of analyzed orthologous protein sets from all three domains of life. However, approximately 1% of the proteins from each group of species deviates from this pattern and instead shows variation that is consistent with an acceleration of the rate of amino-acid substitution, which may be due to functional diversification. Most of the putative functionally diversified proteins from all three species groups are predicted to function at the periphery of the cells and mediate their interaction with the environment. Relative rates of protein evolution are remarkably constant for the three species groups analyzed here. Deviations from this rate constancy are probably due to changes in selective constraints associated with diversification between orthologs. Functional diversification between orthologs is thought to be a relatively rare event. However, the resolution afforded by the test designed specifically for genomic-scale datasets allowed us to identify numerous cases of possible functional diversification between orthologous proteins.

What determines the rate of sequence evolution?

High rates of amino-acid sequence evolution have sometimes been considered to be diagnostic for genes undergoing adaptive change. However, two recent studies have shown that rapid evolution of amino-acid sequence can also be congruent with neutrality.

Using protein structural information in evolutionary inference: transmembrane proteins.

We present a model of amino acid sequence evolution based on a hidden Markov model that extends to transmembrane proteins previous methods that incorporate protein structural information into phylogenetics. Our model aims to give a better understanding of processes of molecular evolution and to extract structural information from multiple alignments of transmembrane sequences and use such information to improve phylogenetic analyses. This should be of value in phylogenetic studies of transmembrane proteins: for example, mitochondrial proteins have acquired a special importance in phylogenetics and are mostly transmembrane proteins. The improvement in fit to example data sets of our new model relative to less complex models of amino acid sequence evolution is statistically tested. To further illustrate the potential utility of our method, phylogeny estimation is performed on primate CCR5 receptor sequences, sequences of l and m subunits of the light reaction center in purple bacteria, guinea pig sequences with respect to lagomorph and rodent sequences of calcitonin receptor and K-substance receptor, and cetacean sequences of cytochrome b.

Phylogenetic utility of the nuclear gene arginine decarboxylase: an example from Brassicaceae.

Arginine decarboxylase (ADC) is an important enzyme in the production of putrescine and polyamines in plants. It is encoded by a single or low-copy nuclear gene that lacks introns in sequences studied to date. The rate of Adc amino acid sequence evolution is similar to that of ndhF for the angiosperm family studied. Highly conserved regions provide several target sites for PCR priming and sequencing and aid in nucleotide and amino acid sequence alignment across a range of taxonomic levels, while a variable region provides an increased number of potentially informative characters relative to ndhF for the taxa surveyed. The utility of the Adc gene in plant molecular systematic studies is demonstrated by analysis of its partial nucleotide sequences obtained from 13 representatives of Brassicaceae and 3 outgroup taxa, 2 from the mustard oil clade (order Capparales) and 1 from the related order Malvales. Two copies of the Adc gene, Adc1 and Adc2, are found in all members of the Brassicaceae studied to data except the basal genus Aethionema. The resulting Adc gene tree provides robust phylogenetic data regarding relationships within the complex mustard family, as well as independent support for proposed tribal realignments based on other molecular data sets such as those from chloroplast DNA.

Evolution of DNA or amino acid sequences with dependent sites.

A framework is outlined to study the evolution of DNA or amino acid sequences, if sequence sites do not evolve independently. The units of evolution are nonoverlapping subsequences of length l. Each subsequence evolves independently of the others, but within a subsequence the sequences show a Markov order one dependency. We describe an algorithm to mimic the evolution of such sequences. The influence of dependencies between sites on distance estimates and the reliability of tree reconstruction methods is investigated. We show that an inappropriate model of sequence evolution in the tree reconstruction process will lead to a nonempty Felsenstein zone. Finally, we describe a method to infer l from sequence data. Examples from the evolution of DNA sequences as well as from amino acids are given.

Correlated mutations and residue contacts in proteins

The maintenance of protein function and structure constrains the evolution of amino acid sequences. This fact can be exploited to interpret correlated mutations observed in a sequence family as an indication of probable physical contact in three dimensions. Here we present a simple and general method to analyze correlations in mutational behavior between different positions in a multiple sequence alignment. We then use these correlations to predict contact maps for each of 11 protein families and compare the result with the contacts determined by crystallography. For the most strongly correlated residue pairs predicted to be in contact, the prediction accuracy ranges from 37 to 68% and the improvement ratio relative to a random prediction from 1.4 to 5.1. Predicted contact maps can be used as input for the calculation of protein tertiary structure, either from sequence information alone or in combination with experimental information.

Molecular evolutionary analyses of the small and large subunits of ribulose-1,5-bisphosphate carboxylase/oxygenase

The molecular evolution of the amino acid sequences of the mature small and large subunits of ribulose-1,5-bisphosphate carboxylase/oxygense (Rubisco) was determined. The dataset for each subunit consisted of sequences from 39 different taxa of which 22 are represented with sequence information for both subunits. Phylogenetic trees were reconstructed using distance matrix, parsimony and simultaneous alignment and phylogeny methods. For the small subunit, the latter two methods produced similar trees that differed from the topology of the distance matrix tree. For the large subunit, each of the three tree-building methods yielded a distinct tree. Except for the distance matrix small subunit tree, the tree-building methods produced topologies for the small and large subunit sequences from the nonflowering plant taxa that, for the most part, agree with current taxonomic schemes. With the full datasets, the lack of consistency both among the various trees and with conventional taxonomic relationships was most evident with the Rubisco sequences from angiosperms. It is unlikely that current tree-building methods will be able to reconstruct an unambiguous molecular evolution of either of the Rubisco subunits. Molecular trees, regardless of methodology, showed similar topologies for the small and large subunits from the 22 taxa from which both subunits have been sequenced, indicating that the subunits have changed to the same extent over time. In this case, similar trees were formed because only 4 of the 22 taxa were from dicots. Key words: ribulose-1,5-bisphosphate carboxylase/oxygenase, amino acid sequence, molecular evolution, phyletic trees.