Amino Acid Sequence Evolution Research Articles

Expression of defensin genes across house fly (Musca domestica) life history gives insight into immune system subfunctionalization.

Animals encounter diverse microbial communities throughout their lifetime, which exert varying selection pressures. Antimicrobial peptides (AMPs), which lyse or inhibit microbial growth, are a first line of defense against some of these microbes. Here we examine how developmental variation in microbial exposure has affected the evolution of expression and amino acid sequences of Defensins (an ancient class of AMPs) in the house fly (Musca domestica). The house fly is a well-suited model for this work because it trophically associates with varying microbial communities throughout its life history and its genome contains expanded families of AMPs, including Defensins. We identified two subsets of house fly Defensins: one expressed in larvae or pupae, and the other expressed in adults. The amino acid sequences of these two Defensin subsets form distinct monophyletic clades, and they are located in separate gene clusters in the genome. The adult-expressed Defensins evolve faster than larval/pupal Defensins, consistent with different selection pressures across developmental stages. Our results therefore suggest that varied microbial communities encountered across life history can shape the evolutionary trajectories of immune genes.

Intrahost evolution of the HIV-2 capsid correlates with progression to AIDS.

HIV-2 infection will progress to AIDS in most patients without treatment, albeit at approximately half the rate of HIV-1 infection. HIV-2 capsid (p26) amino acid polymorphisms are associated with lower viral loads and enhanced processing of T cell epitopes, which may lead to protective Gag-specific T cell responses common in slower progressors. Lower virus evolutionary rates, and positive selection on conserved residues in HIV-2 env have been associated with slower progression to AIDS. In this study we analysed 369 heterochronous HIV-2 p26 sequences from 12 participants with a median age of 30 years at enrolment. CD4% change over time was used to stratify participants into relative faster and slower progressor groups. We analysed p26 sequence diversity evolution, measured site-specific selection pressures and evolutionary rates, and determined if these evolutionary parameters were associated with progression status. Faster progressors had lower CD4% and faster CD4% decline rates. Median pairwise sequence diversity was higher in faster progressors (5.7x10-3 versus 1.4x10-3 base substitutions per site, P<0.001). p26 evolved under negative selection in both groups (dN/dS=0.12). Median virus evolutionary rates were higher in faster than slower progressors - synonymous rates: 4.6x10-3 vs. 2.3x10-3; and nonsynonymous rates: 6.9x10-4 vs. 2.7x10-4 substitutions/site/year, respectively. Virus evolutionary rates correlated negatively with CD4% change rates (ρ = -0.8, P=0.02), but not CD4% level. The signature amino acid at p26 positions 6, 12 and 119 differed between faster (6A, 12I, 119A) and slower (6G, 12V, 119P) progressors. These amino acid positions clustered near to the TRIM5α/p26 hexamer interface surface. p26 evolutionary rates were associated with progression to AIDS and were mostly driven by synonymous substitutions. Nonsynonymous evolutionary rates were an order of magnitude lower than synonymous rates, with limited amino acid sequence evolution over time within hosts. These results indicate HIV-2 p26 may be an attractive therapeutic target.

Epistasis Creates Invariant Sites and Modulates the Rate of Molecular Evolution.

Invariant sites are a common feature of amino acid sequence evolution. The presence of invariant sites is frequently attributed to the need to preserve function through site-specific conservation of amino acid residues. Amino acid substitution models without a provision for invariant sites often fit the data significantly worse than those that allow for an excess of invariant sites beyond those predicted by models that only incorporate rate variation among sites (e.g., a Gamma distribution). An alternative is epistasis between sites to preserve residue interactions that can create invariant sites. Through computer-simulated sequence evolution, we evaluated the relative effects of site-specific preferences and site-site couplings in the generation of invariant sites and the modulation of the rate of molecular evolution. In an analysis of ten major families of protein domains with diverse sequence and functional properties, we find that the negative selection imposed by epistasis creates many more invariant sites than site-specific residue preferences alone. Further, epistasis plays an increasingly larger role in creating invariant sites over longer evolutionary periods. Epistasis also dictates rates of domain evolution over time by exerting significant additional purifying selection to preserve site couplings. These patterns illuminate the mechanistic role of epistasis in the processes underlying observed site invariance and evolutionary rates.

Targeted modulation of protein liquid-liquid phase separation by evolution of amino-acid sequence.

Rationally and efficiently modifying the amino-acid sequence of proteins to control their ability to undergo liquid-liquid phase separation (LLPS) on demand is not only highly desirable, but can also help to elucidate which protein features are important for LLPS. Here, we propose a computational method that couples a genetic algorithm to a sequence-dependent coarse-grained protein model to evolve the amino-acid sequences of phase-separating intrinsically disordered protein regions (IDRs), and purposely enhance or inhibit their capacity to phase-separate. We validate the predicted critical solution temperatures of the mutated sequences with ABSINTH, a more accurate all-atom model. We apply the algorithm to the phase-separating IDRs of three naturally occurring proteins, namely FUS, hnRNPA1 and LAF1, as prototypes of regions that exist in cells and undergo homotypic LLPS driven by different types of intermolecular interaction, and we find that the evolution of amino-acid sequences towards enhanced LLPS is driven in these three cases, among other factors, by an increase in the average size of the amino acids. However, the direction of change in the molecular driving forces that enhance LLPS (such as hydrophobicity, aromaticity and charge) depends on the initial amino-acid sequence. Finally, we show that the evolution of amino-acid sequences to modulate LLPS is strongly coupled to the make-up of the medium (e.g. the presence or absence of RNA), which may have significant implications for our understanding of phase separation within the many-component mixtures of biological systems.

Sequence variation and immunogenicity of the Mycoplasma genitalium MgpB and MgpC adherence proteins during persistent infection of men with non-gonococcal urethritis.

Mycoplasma genitalium is a sexually transmitted bacterial pathogen that infects men and women. Antigenic variation of MgpB and MgpC, the immunodominant adherence proteins of M. genitalium, is thought to contribute to immune evasion and chronic infection. We investigated the evolution of mgpB and mgpC sequences in men with non-gonococcal urethritis persistently infected with M. genitalium, including two men with anti-M. genitalium antibodies at enrollment and two that developed antibodies during follow-up. Each of the four patients was persistently infected with a different strain type and each patient produced antibodies targeting MgpB and MgpC. Amino acid sequence evolution in the variable regions of MgpB and MgpC occurred in all four patients with changes observed in single and multiple variable regions over time. Using the available crystal structure of MgpC of the G37 type strain we found that predicted conformational B cell epitopes localize predominantly to the variable region of MgpC, amino acids that changed during patient infection lie in these epitopes, and variant amino acids are in close proximity to the conserved sialic acid binding pocket. These findings support the hypothesis that sequence variation functions to avoid specific antibodies thereby contributing to persistence in the genital tract.

Epistatic contributions promote the unification of incompatible models of neutral molecular evolution

We introduce a model of amino acid sequence evolution that accounts for the statistical behavior of real sequences induced by epistatic interactions. We base the model dynamics on parameters derived from multiple sequence alignments analyzed by using direct coupling analysis methodology. Known statistical properties such as overdispersion, heterotachy, and gamma-distributed rate-across-sites are shown to be emergent properties of this model while being consistent with neutral evolution theory, thereby unifying observations from previously disjointed evolutionary models of sequences. The relationship between site restriction and heterotachy is characterized by tracking the effective alphabet dynamics of sites. We also observe an evolutionary Stokes shift in the fitness of sequences that have undergone evolution under our simulation. By analyzing the structural information of some proteins, we corroborate that the strongest Stokes shifts derive from sites that physically interact in networks near biochemically important regions. Perspectives on the implementation of our model in the context of the molecular clock are discussed.

A Machine Learning Method for Detecting Autocorrelation of Evolutionary Rates in Large Phylogenies

New species arise from pre-existing species and inherit similar genomes and environments. This predicts greater similarity of the tempo of molecular evolution between direct ancestors and descendants, resulting in autocorrelation of evolutionary rates in the tree of life. Surprisingly, molecular sequence data have not confirmed this expectation, possibly because available methods lack the power to detect autocorrelated rates. Here, we present a machine learning method, CorrTest, to detect the presence of rate autocorrelation in large phylogenies. CorrTest is computationally efficient and performs better than the available state-of-the-art method. Application of CorrTest reveals extensive rate autocorrelation in DNA and amino acid sequence evolution of mammals, birds, insects, metazoans, plants, fungi, parasitic protozoans, and prokaryotes. Therefore, rate autocorrelation is a common phenomenon throughout the tree of life. These findings suggest concordance between molecular and nonmolecular evolutionary patterns, and they will foster unbiased and precise dating of the tree of life.

Dissecting the Genomic Diversification of Late Embryogenesis Abundant (LEA) Protein Gene Families in Plants.

Late embryogenesis abundant (LEA) proteins include eight multigene families that are expressed in response to water loss during seed maturation and in vegetative tissues of desiccation tolerant species. To elucidate LEA proteins evolution and diversification, we performed a comprehensive synteny and phylogenetic analyses of the eight gene families across 60 complete plant genomes. Our integrated comparative genomic approach revealed that synteny conservation and diversification contributed to LEA family expansion and functional diversification in plants. We provide examples that: 1) the genomic diversification of the Dehydrin family contributed to differential evolution of amino acid sequences, protein biochemical properties, and gene expression patterns, and led to the appearance of a novel functional motif in angiosperms; 2) ancient genomic diversification contributed to the evolution of distinct intrinsically disordered regions of LEA_1 proteins; 3) recurrent tandem-duplications contributed to the large expansion of LEA_2; and 4) dynamic synteny diversification played a role on the evolution of LEA_4 and its function on plant desiccation tolerance. Taken together, these results show that multiple evolutionary mechanisms have not only led to genomic diversification but also to structural and functional plasticity among LEA proteins which have jointly contributed to the adaptation of plants to water-limiting environments.

A phylogenomic framework and timescale for comparative studies of tunicates

BackgroundTunicates are the closest relatives of vertebrates and are widely used as models to study the evolutionary developmental biology of chordates. Their phylogeny, however, remains poorly understood, and to date, only the 18S rRNA nuclear gene and mitogenomes have been used to delineate the major groups of tunicates. To resolve their evolutionary relationships and provide a first estimate of their divergence times, we used a transcriptomic approach to build a phylogenomic dataset including all major tunicate lineages, consisting of 258 evolutionarily conserved orthologous genes from representative species.ResultsPhylogenetic analyses using site-heterogeneous CAT mixture models of amino acid sequence evolution resulted in a strongly supported tree topology resolving the relationships among four major tunicate clades: (1) Appendicularia, (2) Thaliacea + Phlebobranchia + Aplousobranchia, (3) Molgulidae, and (4) Styelidae + Pyuridae. Notably, the morphologically derived Thaliacea are confirmed as the sister group of the clade uniting Phlebobranchia + Aplousobranchia within which the precise position of the model ascidian genus Ciona remains uncertain. Relaxed molecular clock analyses accommodating the accelerated evolutionary rate of tunicates reveal ancient diversification (~ 450–350 million years ago) among the major groups and allow one to compare their evolutionary age with respect to the major vertebrate model lineages.ConclusionsOur study represents the most comprehensive phylogenomic dataset for the main tunicate lineages. It offers a reference phylogenetic framework and first tentative timescale for tunicates, allowing a direct comparison with vertebrate model species in comparative genomics and evolutionary developmental biology studies.

Natural vs. random protein sequences: Discovering combinatorics properties on amino acid words

Casual mutations and natural selection have driven the evolution of protein amino acid sequences that we observe at present in nature. The question about which is the dominant force of proteins evolution is still lacking of an unambiguous answer. Casual mutations tend to randomize protein sequences while, in order to have the correct functionality, one expects that selection mechanisms impose rigid constraints on amino acid sequences. Moreover, one also has to consider that the space of all possible amino acid sequences is so astonishingly large that it could be reasonable to have a well tuned amino acid sequence indistinguishable from a random one.In order to study the possibility to discriminate between random and natural amino acid sequences, we introduce different measures of association between pairs of amino acids in a sequence, and apply them to a dataset of 1047 natural protein sequences and 10,470 random sequences, carefully generated in order to preserve the relative length and amino acid distribution of the natural proteins. We analyze the multidimensional measures with machine learning techniques and show that, to a reasonable extent, natural protein sequences can be differentiated from random ones.

Trends in substitution models of molecular evolution.

Substitution models of evolution describe the process of genetic variation through fixed mutations and constitute the basis of the evolutionary analysis at the molecular level. Almost 40 years after the development of first substitution models, highly sophisticated, and data-specific substitution models continue emerging with the aim of better mimicking real evolutionary processes. Here I describe current trends in substitution models of DNA, codon and amino acid sequence evolution, including advantages and pitfalls of the most popular models. The perspective concludes that despite the large number of currently available substitution models, further research is required for more realistic modeling, especially for DNA coding and amino acid data. Additionally, the development of more accurate complex models should be coupled with new implementations and improvements of methods and frameworks for substitution model selection and downstream evolutionary analysis.

Pyvolve: A Flexible Python Module for Simulating Sequences along Phylogenies.

We introduce Pyvolve, a flexible Python module for simulating genetic data along a phylogeny using continuous-time Markov models of sequence evolution. Easily incorporated into Python bioinformatics pipelines, Pyvolve can simulate sequences according to most standard models of nucleotide, amino-acid, and codon sequence evolution. All model parameters are fully customizable. Users can additionally specify custom evolutionary models, with custom rate matrices and/or states to evolve. This flexibility makes Pyvolve a convenient framework not only for simulating sequences under a wide variety of conditions, but also for developing and testing new evolutionary models. Pyvolve is an open-source project under a FreeBSD license, and it is available for download, along with a detailed user-manual and example scripts, from http://github.com/sjspielman/pyvolve.

Integrative analysis of young genes, positively selected genes and lncRNAs in the development of Drosophila melanogaster.

BackgroundYoung genes and genes under positive selection commonly contribute to adaptive phenotypic evolution. Early developmental stages are very important for establishing phenotypes, which might be helpful for studying the evolutionary patterns of these rapidly evolving genes.ResultsHere, we performed a weighted gene co-expression network analysis to identify modules of co-expressed genes at different stages of Drosophila melanogaster development. We found that young genes, including duplicated, orphan, and young lncRNA genes, are significantly enriched among modules associated with specific developmental stages. In addition, genes undergoing rapid amino acid sequence evolution driven by positive selection showed a similar proportion of essentiality with other genes, and enrichment in modules for specific developmental stages.ConclusionsOur integrative analysis revealed important roles for the origin of new genes and rapid amino acid sequence evolution in development that may account for specific phenotype evolution in Drosophila melanogaster.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-014-0241-9) contains supplementary material, which is available to authorized users.

Protein subcellular relocalization of duplicated genes in Arabidopsis.

Gene duplications during eukaroytic evolution, by successive rounds of polyploidy and by smaller scale duplications, have provided an enormous reservoir of new genes for the evolution of new functions. Preservation of many duplicated genes can be ascribed to changes in sequences, expression patterns, and functions. Protein subcellular relocalization (protein targeting to a new location within the cell) is another way that duplicated genes can diverge. We studied subcellular relocalization of gene pairs duplicated during the evolution of the Brassicaceae including gene pairs from the alpha whole genome duplication that occurred at the base of the family. We analyzed experimental localization data from green fluorescent protein experiments for 128 duplicate pairs in Arabidopsis thaliana, revealing 19 pairs with subcellular relocalization. Many more of the duplicate pairs with relocalization than with the same localization showed an accelerated rate of amino acid sequence evolution in one duplicate, and one gene showed evidence for positive selection. We studied six duplicate gene pairs in more detail. We used gene family analysis with several pairs to infer which gene shows relocalization. We identified potential sequence mutations through comparative analysis that likely result in relocalization of two duplicated gene products. We show that four cases of relocalization have new expression patterns, compared with orthologs in outgroup species, including two with novel expression in pollen. This study provides insights into subcellular relocalization of evolutionarily recent gene duplicates and features of genes whose products have been relocalized.

A recurring syndrome of accelerated plastid genome evolution in the angiosperm tribe Sileneae (Caryophyllaceae)

In flowering plants, plastid genomes are generally conserved, exhibiting slower rates of sequence evolution than the nucleus and little or no change in structural organization. However, accelerated plastid genome evolution has occurred in scattered angiosperm lineages. For example, some species within the genus Silene have experienced a suite of recent changes to their plastid genomes, including inversions, shifts in inverted repeat boundaries, large indels, intron losses, and rapid rates of amino acid sequence evolution in a subset of protein genes, with the most extreme divergence occurring in the protease gene clpP. To investigate the relationship between the rates of sequence and structural evolution, we sequenced complete plastid genomes from three species (Silene conoidea, S. paradoxa, and Lychnis chalcedonica), representing independent lineages within the tribe Sileneae that were previously shown to have accelerated rates of clpP evolution. We found a high degree of parallel evolution. Elevated rates of amino acid substitution have occurred repeatedly in the same subset of plastid genes and have been accompanied by a recurring pattern of structural change, including cases of identical inversions and intron loss. This “syndrome” of changes was not observed in the closely related outgroup Agrostemma githago or in the more slowly evolving Silene species that were sequenced previously. Although no single mechanism has yet been identified to explain the correlated suite of changes in plastid genome sequence and structure that has occurred repeatedly in angiosperm evolution, we discuss a possible mixture of adaptive and non-adaptive forces that may be responsible.

A New Formulation of Protein Evolutionary Models that Account for Structural Constraints

Despite the importance of a thermodynamically stable structure with a conserved fold for protein function, almost all evolutionary models neglect site-site correlations that arise from physical interactions between neighboring amino acid sites. This is mainly due to the difficulty in formulating a computationally tractable model since rate matrices can no longer be used. Here, we introduce a general framework, based on factor graphs, for constructing probabilistic models of protein evolution with site interdependence. Conveniently, efficient approximate inference algorithms, such as Belief Propagation, can be used to calculate likelihoods for these models. We fit an amino acid substitution model of this type that accounts for both solvent accessibility and site-site correlations. Comparisons of the new model with rate matrix models and alternative structure-dependent models demonstrate that it better fits the sequence data. We also examine evolution within a family of homohexameric enzymes and find that site-site correlations between most contacting subunits contribute to a higher likelihood. In addition, we show that the new substitution model has a similar mathematical form to the one introduced in Rodrigue et al. (Rodrigue N, Lartillot N, Bryant D, Philippe H. 2005. Site interdependence attributed to tertiary structure in amino acid sequence evolution. Gene 347:207-217), although with different parameter interpretations and values. We also perform a statistical analysis of the effects of amino acids at neighboring sites on substitution probabilities and find a significant perturbation of most probabilities, further supporting the significant role of site-site interactions in protein evolution and motivating the development of new evolutionary models similar to the one described here. Finally, we discuss possible extensions and applications of the new substitution model.

Evolution of the salivary apyrases of blood-feeding arthropods

Phylogenetic analyses of three families of arthropod apyrases were used to reconstruct the evolutionary relationships of salivary-expressed apyrases, which have an anti-coagulant function in blood-feeding arthropods. Members of the 5′nucleotidase family were recruited for salivary expression in blood-feeding species at least five separate times in the history of arthropods, while members of the Cimex-type apyrase family have been recruited at least twice. In spite of these independent events of recruitment for salivary function, neither of these families showed evidence of convergent amino acid sequence evolution in salivary-expressed members. On the contrary, in the 5′-nucleotide family, salivary-expressed proteins conserved ancestral amino acid residues to a significantly greater extent than related proteins without salivary function, implying parallel evolution by conservation of ancestral characters. This unusual pattern of sequence evolution suggests the hypothesis that purifying selection favoring conservation of ancestral residues is particularly strong in salivary-expressed members of the 5′-nucleotidase family of arthropods because of constraints arising from expression within the vertebrate host.

Molecular evolution of the pigmentation gene melanocortin-1 receptor in rodents

Adaptive variation in the melanocortin 1-receptor gene (MC1R), a key locus in melanogenesis, has been identified in some species of rodents. However, in others, MC1R has no causative role in pigmentation phenotypes despite their coat color variation. In this study, we characterized the rates and patterns of MC1R nucleotide and amino acid sequence evolution and, particularly, selective pressures in the separated domains of the protein using a comparative analysis of 43 species representing three major lineages of rodents with variable coat colors. We found high amino acid variation (44% of sites) throughout the protein. Most substitutions were observed in extracellular and transmembrane domains; the intracellular segment was conserved across species. Pairwise non-synonymous substitutions did not vary significantly in different domains among the rodent lineages - i.e., variation was not associated with phylogenetic distance. Phylogeny-based likelihood analysis suggested that purifying selection has mostly shaped the evolutionary course of MC1R. However, a high proportion of sites (27%) were under relaxation of functional constraints (ω = 0.38), and four sites (3, 14, 26, and 251) clearly evolved under positive selection (ω ≅ 2.9). Thus, our data indicate a high proportion of sites evolving under relaxed evolutionary constraints, which might indicate the evolvability of the system in the generation of adaptive changes in specific taxa in rodent lineages.

REvolver: Modeling Sequence Evolution under Domain Constraints

Simulating the change of protein sequences over time in a biologically realistic way is fundamental for a broad range of studies with a focus on evolution. It is, thus, problematic that typically simulators evolve individual sites of a sequence identically and independently. More realistic simulations are possible; however, they are often prohibited by limited knowledge concerning site-specific evolutionary constraints or functional dependencies between amino acids. As a consequence, a protein's functional and structural characteristics are rapidly lost in the course of simulated evolution. Here, we present REvolver (www.cibiv.at/software/revolver), a program that simulates protein sequence alteration such that evolutionarily stable sequence characteristics, like functional domains, are maintained. For this purpose, REvolver recruits profile hidden Markov models (pHMMs) for parameterizing site-specific models of sequence evolution in an automated fashion. pHMMs derived from alignments of homologous proteins or protein domains capture information regarding which sequence sites remained conserved over time and where in a sequence insertions or deletions are more likely to occur. Thus, they describe constraints on the evolutionary process acting on these sequences. To demonstrate the performance of REvolver as well as its applicability in large-scale simulation studies, we evolved the entire human proteome up to 1.5 expected substitutions per site. Simultaneously, we analyzed the preservation of Pfam and SMART domains in the simulated sequences over time. REvolver preserved 92% of the Pfam domains originally present in the human sequences. This value drops to 15% when traditional models of amino acid sequence evolution are used. Thus, REvolver represents a significant advance toward a realistic simulation of protein sequence evolution on a proteome-wide scale. Further, REvolver facilitates the simulation of a protein family with a user-defined domain architecture at the root.

The Population Genomics of a Fast Evolver: High Levels of Diversity, Functional Constraint, and Molecular Adaptation in the Tunicate Ciona intestinalis

Phylogenomics has revealed the existence of fast-evolving animal phyla in which the amino acid substitution rate, averaged across many proteins, is consistently higher than in other lineages. The reasons for such differences in proteome-wide evolutionary rates are still unknown, largely because only a handful of species offer within-species genomic data from which molecular evolutionary processes can be deduced. In this study, we use next-generation sequencing technologies and individual whole-transcriptome sequencing to gather extensive polymorphism sequence data sets from Ciona intestinalis. Ciona is probably the best-characterized member of the fast-evolving Urochordata group (tunicates), which was recently identified as the sister group of the slow-evolving vertebrates. We introduce and validate a maximum-likelihood framework for single-nucleotide polymorphism and genotype calling, based on high-throughput short-read typing. We report that the C. intestinalis proteome is characterized by a high level of within-species diversity, efficient purifying selection, and a substantial percentage of adaptive amino acid substitutions. We conclude that the increased rate of amino acid sequence evolution in tunicates, when compared with vertebrates, is the consequence of both a 2–6 times higher per-year mutation rate and prevalent adaptive evolution.