Rate Of Protein Sequence Evolution Research Articles

Genome assembly of primitive cultivated potato Solanum stenotomum provides insights into potato evolution.

Genetic diversity is the raw material for germplasm enhancement. Landraces and wild species relatives of potato, which contain a rich gene pool of valuable agronomic traits, can provide insights into the genetic diversity behind the adaptability of the common potato. The diploid plant, Solanum stenotomum (Sst), is believed to have an ancestral relationship with modern potato cultivars and be a potential source of resistance against disease. Sequencing of the Sst genome generated an assembly of 852.85 Mb (N50 scaffold size, 3.7 Mb). Pseudomolecule construction anchored 788.75 Mb of the assembly onto 12 pseudochromosomes, with an anchor rate of 92.4%. Genome annotation yielded 41,914 high-confidence protein-coding gene models and comparative analyses with closely related Solanaceae species identified 358 Sst-specific gene families, 885 gene families with expansion along the Sst lineage, and 149 genes experiencing accelerated rates of protein sequence evolution in Sst, the functions of which were mainly associated with defense responses, particularly against bacterial and fungal infection. Insights into the Sst genome and the genomic variation of cultivated potato taxa are valuable in elaborating the impact of potato evolution in early landrace diploid and facilitate modern potato breeding.

Evolutionary dynamics of sex-biased genes expressed in cricket brains and gonads.

Sex-biased gene expression, particularly sex-biased expression in the gonad, has been linked to rates of protein sequence evolution (nonsynonymous to synonymous substitutions, dN/dS) in animals. However, in insects, sex-biased expression studies remain centred on a few holometabolous species. Moreover, other major tissue types such as the brain remain underexplored. Here, we studied sex-biased gene expression and protein evolution in a hemimetabolous insect, the cricket Gryllus bimaculatus. We generated novel male and female RNA-seq data for two sexual tissue types, the gonad and somatic reproductive system, and for two core components of the nervous system, the brain and ventral nerve cord. From a genome-wide analysis, we report several core findings. Firstly, testis-biased genes had accelerated evolution, as compared to ovary-biased and unbiased genes, which was associated with positive selection events. Secondly, although sex-biased brain genes were much less common than for the gonad, they exhibited a striking tendency for rapid protein sequence evolution, an effect that was stronger for the female than male brain. Further, some sex-biased brain genes were linked to sexual functions and mating behaviours, which we suggest may have accelerated their evolution via sexual selection. Thirdly, a tendency for narrow cross-tissue expression breadth, suggesting low pleiotropy, was observed for sex-biased brain genes, suggesting relaxed purifying selection, which we speculate may allow enhanced freedom to evolve adaptive protein functional changes. The findings of rapid evolution of testis-biased genes and male and female-biased brain genes are discussed with respect to pleiotropy, positive selection and the mating biology of this cricket.

Parallels between natural selection in the cold-adapted crop-wild relative Tripsacum dactyloides and artificial selection in temperate adapted maize.

Artificial selection has produced varieties of domesticated maize that thrive in temperate climates around the world. However, the direct progenitor of maize, teosinte, is indigenous only to a relatively small range of tropical and subtropical latitudes and grows poorly or not at all outside of this region. Tripsacum, a sister genus to maize and teosinte, is naturally endemic to the majority of areas in the western hemisphere where maize is cultivated. A full-length reference transcriptome for Tripsacum dactyloides generated using long-read Iso-Seq data was used to characterize independent adaptation to temperate climates in this clade. Genes related to phospholipid biosynthesis, a critical component of cold acclimation in other cold-adapted plant lineages, were enriched among those genes experiencing more rapid rates of protein sequence evolution in T.dactyloides. In contrast with previous studies of parallel selection, we find that there is a significant overlap between the genes that were targets of artificial selection during the adaptation of maize to temperate climates and those that were targets of natural selection in temperate-adapted T.dactyloides. Genes related to growth, development, response to stimulus, signaling, and organelles were enriched in the set of genes identified as both targets of natural and artificial selection.

Increased rates of protein evolution and asymmetric deceleration after the whole-genome duplication in yeasts

BackgroundWhole-genome duplication (WGD) events have shaped the genomes of eukaryotic organisms. Relaxed selection after duplication along with inherent functional constraints are thought to determine the fate of the paralogs and, ultimately, the evolution of gene function. Here, we investigated the rate of protein evolution (as measured by dN/dS ratios) before and after the WGD in the hemiascomycete yeasts, and the way in which changes in such rates relate to molecular and biological function.ResultsFor most groups of orthologous genes (81%) we observed a change in the rates of evolution after genome duplication. Genes with atypically-low dN/dS ratio before the WGD were prone to increase their rates of evolution after duplication. Importantly, the paralogs were often different in their rates of evolution after the WGD (50% cases), however, this was more consistent with an asymmetric deceleration in the protein-evolution rates, rather than an asymmetric increase of the initial rates. Functional-category analysis showed that regulatory proteins such as protein kinases and transcription factors were enriched in genes that increase their rates of evolution after the WGD. While changes in the rate of protein-sequence evolution were associated to protein abundance, content of disordered regions, and contribution to fitness, these features were an attribute of specific functional classes.ConclusionsOur results indicate that strong purifying selection in ancestral pre-duplication sequences is a strong predictor of increased rates after the duplication in yeasts and that asymmetry in evolution rate is established during the deceleration phase. In addition, changes in the rates at which paralogous sequences evolve before and after WGD are different for specific protein functions; increased rates of protein evolution after duplication occur preferentially in specific protein functions.

Determinants of the rate of protein sequence evolution.

The rate and mechanism of protein sequence evolution have been central questions in evolutionary biology since the 1960s. Although the rate of protein sequence evolution depends primarily on the level of functional constraint, exactly what determines functional constraint has remained unclear. The increasing availability of genomic data has enabled much needed empirical examinations on the nature of functional constraint. These studies found that the evolutionary rate of a protein is predominantly influenced by its expression level rather than functional importance. A combination of theoretical and empirical analyses has identified multiple mechanisms behind these observations and demonstrated a prominent role in protein evolution of selection against errors in molecular and cellular processes.

Local packing density is the main structural determinant of the rate of protein sequence evolution at site level.

Functional and biophysical constraints result in site-dependent patterns of protein sequence variability. It is commonly assumed that the key structural determinant of site-specific rates of evolution is the Relative Solvent Accessibility (RSA). However, a recent study found that amino acid substitution rates correlate better with two Local Packing Density (LPD) measures, the Weighted Contact Number (WCN) and the Contact Number (CN), than with RSA. This work aims at a more thorough assessment. To this end, in addition to substitution rates, we considered four other sequence variability scores, four measures of solvent accessibility (SA), and other CN measures. We compared all properties for each protein of a structurally and functionally diverse representative dataset of monomeric enzymes. We show that the best sequence variability measures take into account phylogenetic tree topology. More importantly, we show that both LPD measures (WCN and CN) correlate better than all of the SA measures, regardless of the sequence variability score used. Moreover, the independent contribution of the best LPD measure is approximately four times larger than that of the best SA measure. This study strongly supports the conclusion that a site's packing density rather than its solvent accessibility is the main structural determinant of its rate of evolution.

Cytonuclear Interactions and Relaxed Selection Accelerate Sequence Evolution in Organelle Ribosomes

Many mitochondrial and plastid protein complexes contain subunits that are encoded in different genomes. In animals, nuclear-encoded mitochondrial proteins often exhibit rapid sequence evolution, which has been hypothesized to result from selection for mutations that compensate for changes in interacting subunits encoded in mutation-prone animal mitochondrial DNA. To test this hypothesis, we analyzed nuclear genes encoding cytosolic and organelle ribosomal proteins in flowering plants. The model angiosperm genus Arabidopsis exhibits low organelle mutation rates, typical of most plants. Nevertheless, we found that (nuclear-encoded) subunits of organelle ribosomes in Arabidopsis have higher amino acid sequence polymorphism and divergence than their counterparts in cytosolic ribosomes, suggesting that organelle ribosomes experience relaxed functional constraint. However, the observed difference between organelle and cytosolic ribosomes was smaller than in animals and could be partially attributed to rapid evolution in N-terminal organelle-targeting peptides that are not involved in ribosome function. To test the role of organelle mutation more directly, we used transcriptomic data from an angiosperm genus (Silene) with highly variable rates of organelle genome evolution. We found that Silene species with unusually fast-evolving mitochondrial and plastid DNA exhibited increased amino acid sequence divergence in ribosomal proteins targeted to the organelles but not in those that function in cytosolic ribosomes. Overall, these findings support the hypothesis that rapid organelle genome evolution has selected for compensatory mutations in nuclear-encoded proteins. We conclude that coevolution between interacting subunits encoded in different genomic compartments within the eukaryotic cell is an important determinant of variation in rates of protein sequence evolution.

The Expansion and Functional Diversification of the Mammalian Ribonuclease A Superfamily Epitomizes the Efficiency of Multigene Families at Generating Biological Novelty

The ribonuclease (RNase) A superfamily is a vertebrate-specific gene family. Because of a massive expansion that occurred during the early mammalian evolution, extant mammals in general have much more RNase genes than nonmammalian vertebrates. Mammalian RNases have been associated with diverse physiological functions including digestion, cytotoxicity, angiogenesis, male reproduction, and host defense. However, it is still uncertain when their expansion occurred and how a wide array of functions arose during their evolution. To answer these questions, we generate a compendium of all RNase genes identified in 20 complete mammalian genomes including the platypus, Ornithorhynchus anatinus. Using this, we delineate 13 ancient RNase gene lineages that arose before the divergence between the monotreme and the other mammals (∼220 Ma). These 13 ancient gene lineages are differentially retained in the 20 mammals, and the rate of protein sequence evolution is highly variable among them, which suggest that they have undergone extensive functional diversification. In addition, we identify 22 episodes of recent expansion of RNase genes, many of which have signatures of adaptive functional differentiation. Exemplifying this, bursts of gene duplication occurred for the RNase1, RNase4, and RNase5 genes of the little brown bat (Myotis lucifugus), which might have contributed to the species’ effective defense against heavier pathogen loads caused by its communal roosting behavior. Our study illustrates how host-defense systems can generate new functions efficiently by employing a multigene family, which is crucial for a host organism to adapt to its ever-changing pathogen environment.

The impact of trans-regulation on the evolutionary rates of metazoan proteins.

Transcription factor (TF) and microRNA (miRNA) are two crucial trans-regulatory factors that coordinately control gene expression. Understanding the impacts of these two factors on the rate of protein sequence evolution is of great importance in evolutionary biology. While many biological factors associated with evolutionary rate variations have been studied, evolutionary analysis of simultaneously accounting for TF and miRNA regulations across metazoans is still uninvestigated. Here, we provide a series of statistical analyses to assess the influences of TF and miRNA regulations on evolutionary rates across metazoans (human, mouse and fruit fly). Our results reveal that the negative correlations between trans-regulation and evolutionary rates hold well across metazoans, but the strength of TF regulation as a rate indicator becomes weak when the other confounding factors that may affect evolutionary rates are controlled. We show that miRNA regulation tends to be a more essential indicator of evolutionary rates than TF regulation, and the combination of TF and miRNA regulations has a significant dependent effect on protein evolutionary rates. We also show that trans-regulation (especially miRNA regulation) is much more important in human/mouse than in fruit fly in determining protein evolutionary rates, suggesting a considerable variation in rate determinants between vertebrates and invertebrates.

Functional and evolutionary analysis of the genome of an obligate fungal symbiont.

Nutritional symbionts of insects include some of the most bizarre genomes studied to date, with extremely reduced size, biased base composition, and limited metabolic abilities. A monophyletic group of aphids within the subfamily Cerataphidinae have lost the bacterial symbiont common to all other Aphididae (Buchnera aphidicola), which have been replaced by a eukaryotic one, the yeast-like symbiont (YLS). As symbionts are expected to experience reduced effective population size (Ne) and largely clonal life cycles, we used this system as a model to test the hypothesis that chronically high levels of genetic drift will result in an increase in size of a eukaryotic genome. We sequenced the genome of the YLS of the aphid Cerataphis brasiliensis and observed elevated rates of protein sequence evolution and intron proliferation in YLS orthologs relative to those of its closest-sequenced relative, consistent with predictions. A moderate amount of repetitive DNA was found along with evidence of directed mutation to prevent proliferation of repetitive elements. Despite increased intron numbers, the overall genome structure appears not to have undergone massive expansion and is around 25 Mb in size. Compared with Buchnera, the YLS appears to have a much broader metabolic repertoire, though many gene families have been reduced in the YLS relative to related fungi. The patterns observed in the YLS genome suggest that its symbiotic lifestyle is permissive to intron proliferation and accelerated sequence evolution, though other factors appear to limit its overall genome expansion.

Differential requirements for mRNA folding partially explain why highly expressed proteins evolve slowly

The cause of the tremendous among-protein variation in the rate of sequence evolution is a central subject of molecular evolution. Expression level has been identified as a leading determinant of this variation among genes encoded in the same genome, but the underlying mechanisms are not fully understood. We here propose and demonstrate that a requirement for stronger folding of more abundant mRNAs results in slower evolution of more highly expressed genes and proteins. Specifically, we show that: (i) the higher the expression level of a gene, the greater the selective pressure for its mRNA to fold; (ii) random mutations are more likely to decrease mRNA folding when occurring in highly expressed genes than in lowly expressed genes; and (iii) amino acid substitution rate is negatively correlated with mRNA folding strength, with or without the control of expression level. Furthermore, synonymous (d(S)) and nonsynonymous (d(N)) nucleotide substitution rates are both negatively correlated with mRNA folding strength. However, counterintuitively, d(S) and d(N) are differentially constrained by selection for mRNA folding, resulting in a significant correlation between mRNA folding strength and d(N)/d(S), even when gene expression level is controlled. The direction and magnitude of this correlation is determined primarily by the G+C frequency at third codon positions. Together, these findings explain why highly expressed genes evolve slowly, demonstrate a major role of natural selection at the mRNA level in constraining protein evolution, and reveal a previously unrecognized and unexpected form of nonprotein-level selection that impacts d(N)/d(S).

Evolution at Two Levels in Fire Ants: The Relationship between Patterns of Gene Expression and Protein Sequence Evolution

Variation in protein sequence and gene expression each contribute to phenotypic diversity, and may be subject to similar selective pressures. Eusocial insects are particularly useful for investigating the evolutionary link between protein sequence and condition-dependent patterns of gene expression because gene expression plays a central role in determining differences between eusocial insect sexes and castes. We investigated the relationship between protein coding sequence evolution and gene expression patterns in the fire ants Solenopsis invicta, S. richteri, and their hybrids to gain greater insight into how selection jointly operates on gene expression and coding sequence. We found that genes with high expression variability within castes and sexes were frequently differentially expressed between castes and sexes, as well as between species and hybrids. These results indicate that genes showing high variation in expression in one context also tend to show high variation in expression in other contexts. Our analyses further revealed that variation in both intra- and interspecific gene expression was positively associated with rate of protein sequence evolution in Solenopsis. This suggests that selective constraints on a gene operate both at the level of protein sequence and at the level of gene expression regulation. Overall, our study provides one of the strongest demonstrations that selective constraints mediate both protein sequence evolution and gene expression variability across different biological contexts and timescales.

Highlight--Tracking Adaptation's Role: Do Larger Populations Evolve Faster?

Highlight--Tracking Adaptation's Role: Do Larger Populations Evolve Faster?

Protein coalitions in a core mammalian biochemical network linked by rapidly evolving proteins

BackgroundCellular ATP levels are generated by glucose-stimulated mitochondrial metabolism and determine metabolic responses, such as glucose-stimulated insulin secretion (GSIS) from the β-cells of pancreatic islets. We describe an analysis of the evolutionary processes affecting the core enzymes involved in glucose-stimulated insulin secretion in mammals. The proteins involved in this system belong to ancient enzymatic pathways: glycolysis, the TCA cycle and oxidative phosphorylation.ResultsWe identify two sets of proteins, or protein coalitions, in this group of 77 enzymes with distinct evolutionary patterns. Members of the glycolysis, TCA cycle, metabolite transport, pyruvate and NADH shuttles have low rates of protein sequence evolution, as inferred from a human-mouse comparison, and relatively high rates of evolutionary gene duplication. Respiratory chain and glutathione pathway proteins evolve faster, exhibiting lower rates of gene duplication. A small number of proteins in the system evolve significantly faster than co-pathway members and may serve as rapidly evolving adapters, linking groups of co-evolving genes.ConclusionsOur results provide insights into the evolution of the involved proteins. We find evidence for two coalitions of proteins and the role of co-adaptation in protein evolution is identified and could be used in future research within a functional context.

Measuring the evolutionary rate of protein–protein interaction

Despite our extensive knowledge about the rate of protein sequence evolution for thousands of genes in hundreds of species, the corresponding rate of protein function evolution is virtually unknown, especially at the genomic scale. This lack of knowledge is primarily because of the huge diversity in protein function and the consequent difficulty in gauging and comparing rates of protein function evolution. Nevertheless, most proteins function through interacting with other proteins, and protein-protein interaction (PPI) can be tested by standard assays. Thus, the rate of protein function evolution may be measured by the rate of PPI evolution. Here, we experimentally examine 87 potential interactions between Kluyveromyces waltii proteins, whose one to one orthologs in the related budding yeast Saccharomyces cerevisiae have been reported to interact. Combining our results with available data from other eukaryotes, we estimate that the evolutionary rate of protein interaction is (2.6 ± 1.6) × 10(-10) per PPI per year, which is three orders of magnitude lower than the rate of protein sequence evolution measured by the number of amino acid substitutions per protein per year. The extremely slow evolution of protein molecular function may account for the remarkable conservation of life at molecular and cellular levels and allow for studying the mechanistic basis of human disease in much simpler organisms.

Impact of translational error‐induced and error‐free misfolding on the rate of protein evolution

What determines the rate of protein evolution is a fundamental question in biology. Recent genomic studies revealed a surprisingly strong anticorrelation between the expression level of a protein and its rate of sequence evolution. This observation is currently explained by the translational robustness hypothesis in which the toxicity of translational error-induced protein misfolding selects for higher translational robustness of more abundant proteins, which constrains sequence evolution. However, the impact of error-free protein misfolding has not been evaluated. We estimate that a non-negligible fraction of misfolded proteins are error free and demonstrate by a molecular-level evolutionary simulation that selection against protein misfolding results in a greater reduction of error-free misfolding than error-induced misfolding. Thus, an overarching protein-misfolding-avoidance hypothesis that includes both sources of misfolding is superior to the translational robustness hypothesis. We show that misfolding-minimizing amino acids are preferentially used in highly abundant yeast proteins and that these residues are evolutionarily more conserved than other residues of the same proteins. These findings provide unambiguous support to the role of protein-misfolding-avoidance in determining the rate of protein sequence evolution.

Molecular evolution of mammalian incretin hormone genes

Incretin hormones are encoded by two different genes in the human genome, the proglucagon ( GCG) and the glucose-dependent insulinotropic polypeptide ( GIP) genes. To better understand the evolution of incretin hormones, as well as the potential for the evolution of species-specific functions for these peptides, we have identified and characterized the genes for these hormones from the genomes of 35 mammalian species, as well as from the genomes of a few non-mammalian vertebrates. Both proglucagon and GIP were found to be single-copy genes in mammals, and exist in stable genomic neighborhoods with conserved flanking gene order. The exon–intron structure of the genes has been conserved within mammals, although variation in the rate of protein sequence evolution for the peptide hormones was observed. Glucagon and GLP-1 sequences are largely invariant among mammals, except the glucagon sequences from hystricomorph rodents. Previous work has shown that the change in glucagon sequences in hystricomorph rodents is associated with the evolution of species-specific functions. GLP-2 sequences have evolved most rapidly, while GIP evolved at an intermediate rate, although both show punctuated rates with a common very rapid phase of evolution on the early mammalian lineage. These observations suggest that GIP and GLP-2 are more likely to have evolved species-specific functions in subgroups of mammals.

Multifunctionality dominantly determines the rate of human housekeeping and tissue specific interacting protein evolution

Elucidation of the determinants of the rate of protein sequence evolution is one of the great challenges in evolutionary biology. It has been proposed that housekeeping genes are evolutionarily slower than tissue specific genes. In the present communication, we have examined different determinants that influence the evolutionary rate variation in human housekeeping and tissue specific proteins present in protein–protein interaction network. Studies on yeast proteome, revealed a predominant role of protein connectivity in determining the rate of protein evolution. However, in human, we did not observe any significant influence of protein connectivity on its evolutionary rate. Rather, a significant impact of the proportion of protein's interacting length (amount of protein interface involved in interaction with its partners), expression level and multifunctionality has been observed in determining the rate of protein evolution. We also observed that multi interface proteins are evolutionarily conserved between housekeeping and tissue specific genes and it has been found that the average number of biological processes they associated in these two sets of genes is similar. Moreover, single interface proteins in housekeeping genes evolve more slowly as compared to tissue specific genes owing to their involvement in different number of biological processes. Partial correlation analysis suggests that the relative importance of three individual factors in determining the evolutionary rate variation between housekeeping and tissue specific proteins is in the order of protein multifunctionality > protein expression level > interacting protein length.

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.

Comparable contributions of structural-functional constraints and expression level to the rate of protein sequence evolution

BackgroundProteins show a broad range of evolutionary rates. Understanding the factors that are responsible for the characteristic rate of evolution of a given protein arguably is one of the major goals of evolutionary biology. A long-standing general assumption used to be that the evolution rate is, primarily, determined by the specific functional constraints that affect the given protein. These constrains were traditionally thought to depend both on the specific features of the protein's structure and its biological role. The advent of systems biology brought about new types of data, such as expression level and protein-protein interactions, and unexpectedly, a variety of correlations between protein evolution rate and these variables have been observed. The strongest connections by far were repeatedly seen between protein sequence evolution rate and the expression level of the respective gene. It has been hypothesized that this link is due to the selection for the robustness of the protein structure to mistranslation-induced misfolding that is particularly important for highly expressed proteins and is the dominant determinant of the sequence evolution rate.ResultsThis work is an attempt to assess the relative contributions of protein domain structure and function, on the one hand, and expression level on the other hand, to the rate of sequence evolution. To this end, we performed a genome-wide analysis of the effect of the fusion of a pair of domains in multidomain proteins on the difference in the domain-specific evolutionary rates. The mistranslation-induced misfolding hypothesis would predict that, within multidomain proteins, fused domains, on average, should evolve at substantially closer rates than the same domains in different proteins because, within a mutlidomain protein, all domains are translated at the same rate. We performed a comprehensive comparison of the evolutionary rates of mammalian and plant protein domains that are either joined in multidomain proteins or contained in distinct proteins. Substantial homogenization of evolutionary rates in multidomain proteins was, indeed, observed in both animals and plants, although highly significant differences between domain-specific rates remained. The contributions of the translation rate, as determined by the effect of the fusion of a pair of domains within a multidomain protein, and intrinsic, domain-specific structural-functional constraints appear to be comparable in magnitude.ConclusionFusion of domains in a multidomain protein results in substantial homogenization of the domain-specific evolutionary rates but significant differences between domain-specific evolution rates remain. Thus, the rate of translation and intrinsic structural-functional constraints both exert sizable and comparable effects on sequence evolution.ReviewersThis article was reviewed by Sergei Maslov, Dennis Vitkup, Claus Wilke (nominated by Orly Alter), and Allan Drummond (nominated by Joel Bader). For the full reviews, please go to the Reviewers' Reports section.