An evolutionary space-time model with varying among-site dependencies

Inferring Indel Parameters using a Simulation-based Approach

In this study, we present a novel methodology to infer indel parameters from multiple sequence alignments (MSAs) based on simulations. Our algorithm searches for the set of evolutionary parameters describing indel dynamics which best fits a given input MSA. In each step of the search, we use parametric bootstraps and the Mahalanobis distance to estimate how well a proposed set of parameters fits input data. Using simulations, we demonstrate that our methodology can accurately infer the indel parameters for a large variety of plausible settings. Moreover, using our methodology, we show that indel parameters substantially vary between three genomic data sets: Mammals, bacteria, and retroviruses. Finally, we demonstrate how our methodology can be used to simulate MSAs based on indel parameters inferred from real data sets.

Spatial autocorrelation of amino Acid replacement rates in the vasopressin receptor family

Evolutionary rates of sites can be independent of one another or correlated in some fashion. Significant spatial autocorrelation was observed for site amino acid replacement rates in vasopressin receptor family proteins (VPRs). Spatial autocorrelation of rates is the propensity of residues to lie near other residues of similar rate in the folded protein structure. Optimal correlation occurred at a distance suggesting that residues in contact had correlated rates. As another way to study the same phenomenon, VPR was partitioned into >40 x 10 A(3) contiguous spatial clusters for amino acid replacement rate estimation. Partitioning was done without preconception of functional regions of the protein and with a random partition control. Cluster rates exhibited an overdispersed distribution suggesting that rates were not randomly distributed in the spatial partitions. In tests, cluster partitioning improved maximum likelihood and Bayesian likelihood models for VPR evolution. Spatial clusters with outlier rates, or lineage-specific clusters differing in rate, proved to contain VPR features likely to be under selection. Thus the spatial autocorrelation observed is probably not just a statistical finding, but likely has an evolutionary basis in protein function.

Models of coding sequence evolution

Probabilistic models of sequence evolution are in widespread use in phylogenetics and molecular sequence evolution. These models have become increasingly sophisticated and combined with statistical model comparison techniques have helped to shed light on how genes and proteins evolve. Models of codon evolution have been particularly useful, because, in addition to providing a significant improvement in model realism for protein-coding sequences, codon models can also be designed to test hypotheses about the selective pressures that shape the evolution of the sequences. Such models typically assume a phylogeny and can be used to identify sites or lineages that have evolved adaptively. Recently some of the key assumptions that underlie phylogenetic tests of selection have been questioned, such as the assumption that the rate of synonymous changes is constant across sites or that a single phylogenetic tree can be assumed at all sites for recombining sequences. While some of these issues have been addressed through the development of novel methods, others remain as caveats that need to be considered on a case-by-case basis. Here, we outline the theory of codon models and their application to the detection of positive selection. We review some of the more recent developments that have improved their power and utility, laying a foundation for further advances in the modeling of coding sequence evolution.

Investigating protein-coding sequence evolution with probabilistic codon substitution models

This review is motivated by the true explosion in the number of recent studies both developing and ameliorating probabilistic models of codon evolution. Traditionally parametric, the first codon models focused on estimating the effects of selective pressure on the protein via an explicit parameter in the maximum likelihood framework. Likelihood ratio tests of nested codon models armed the biologists with powerful tools, which provided unambiguous evidence for positive selection in real data. This, in turn, triggered a new wave of methodological developments. The new generation of models views the codon evolution process in a more sophisticated way, relaxing several mathematical assumptions. These models make a greater use of physicochemical amino acid properties, genetic code machinery, and the large amounts of data from the public domain. The overview of the most recent advances on modeling codon evolution is presented here, and a wide range of their applications to real data is discussed. On the downside, availability of a large variety of models, each accounting for various biological factors, increases the margin for misinterpretation; the biological meaning of certain parameters may vary among models, and model selection procedures also deserve greater attention. Solid understanding of the modeling assumptions and their applicability is essential for successful statistical data analysis.

Characterizing positive and negative selection and their phylogenetic effects

Total evidence and the use of large datasets to overcome uncertainty are the state of the art in systematic analysis. This assumes that the only true phylogenetic signal is ancestry and that functional, structural, and other factors will not add an alternative signal. Using gene families, where individual codon positions were sorted into bins based upon average-pairwise dN/dS ratio, we show that standard, common phylogenetic methods that were designed for stochastic, neutral, site-independent processes, generate less robust phylogenetic signal for bins with strong negative or positive selection. This was true for phylogenetic reconstruction with parsimony, distance, and likelihood methods. Further, we present a case for the potential existence of systematic functional or structural signal that competes with ancestral signal. For the example of positive selection, we simulate the evolution of sequences through three dimensional lattice constructs with folding constraint and changing binding functionality and show that total evidence for these lattice genes presents trees with functional signal, but that the neutral synonymous sites in these genes show the true ancestral signal. In this case, sequence convergence is promoted by functional convergence.

The quest for natural selection in the age of comparative genomics

Continued genome sequencing has fueled progress in statistical methods for understanding the action of natural selection at the molecular level. This article reviews various statistical techniques (and their applicability) for detecting adaptation events and the functional divergence of proteins. As large-scale automated studies become more frequent, they provide a useful resource for generating biological null hypotheses for further experimental and statistical testing. Furthermore, they shed light on typical patterns of lineage-specific evolution of organisms, on the functional and structural evolution of protein families and on the interplay between the two. More complex models are being developed to better reflect the underlying biological and chemical processes and to complement simpler statistical models. Linking molecular processes to their statistical signatures in genomes can be demanding, and the proper application of statistical models is discussed.

Towards realistic codon models: among site variability and dependency of synonymous and non-synonymous rates

Codon evolutionary models are widely used to infer the selection forces acting on a protein. The non-synonymous to synonymous rate ratio (denoted by Ka/Ks) is used to infer specific positions that are under purifying or positive selection. Current evolutionary models usually assume that only the non-synonymous rates vary among sites while the synonymous substitution rates are constant. This assumption ignores the possibility of selection forces acting at the DNA or mRNA levels. Towards a more realistic description of sequence evolution, we present a model that accounts for among-site-variation of both synonymous and non-synonymous substitution rates. Furthermore, we alleviate the widespread assumption that positions evolve independently of each other. Thus, possible sources of bias caused by random fluctuations in either the synonymous or non-synonymous rate estimations at a single site is removed. Our model is based on two hidden Markov models that operate on the spatial dimension: one describes the dependency between adjacent non-synonymous rates while the other describes the dependency between adjacent synonymous rates. The presented model is applied to study the selection pressure across the HIV-1 genome. The new model better describes the evolution of all HIV-1 genes, as compared to current codon models. Using both simulations and real data analyses, we illustrate that accounting for synonymous rate variability and dependency greatly increases the accuracy of Ka/Ks estimation and in particular of positively selected sites. Finally, we discuss the applicability of the developed model to infer the selection forces in regulatory and overlapping regions of the HIV-1 genome.

Evolution after gene duplication: models, mechanisms, sequences, systems, and organisms

Gene duplication is postulated to have played a major role in the evolution of biological novelty. Here, gene duplication is examined across levels of biological organization in an attempt to create a unified picture of the mechanistic process by which gene duplication can have played a role in generating biodiversity. Neofunctionalization and subfunctionalization have been proposed as important processes driving the retention of duplicate genes. These models have foundations in population genetic theory, which is now being refined by explicit consideration of the structural constraints placed upon genes encoding proteins through physical chemistry. Further, such models can be examined in the context of comparative genomics, where an integration of gene-level evolution and species-level evolution allows an assessment of the frequency of duplication and the fate of duplicate genes. This process, of course, is dependent upon the biochemical role that duplicated genes play in biological systems, which is in turn dependent upon the mechanism of duplication: whole genome duplication involving a co-duplication of interacting partners vs. single gene duplication. Lastly, the role that these processes may have played in driving speciation is examined.

Neighboring-site effects of amino acid mutation

Although the context dependence of nucleotide mutation has been supported by accumulating theoretical and experimental evidence, whether this effect can be extended to amino acid mutation remains obscure. As the amino acid doublets (20 x 20) are much more diverse than their nucleotide counterparts (4 x 4), any attempt to address the neighboring-site effects of amino acid mutation was frustrated by deficient amino acid mutation data. Based on the recently revealed 599,745 mutation sites in 45,137 orthologous proteins, we provide solid evidence for the first time to support the existence of neighboring-site effects in amino acid mutation, which is significantly important to improving the prevalent protein-evolution models.

Using evolutionary information and ancestral sequences to understand the sequence-function relationship in GLP-1 agonists

Glucagon-like peptide-1 (GLP-1) is an incretin hormone with therapeutic potential for type 2 diabetes. A variety of GLP-1 sequences are known from amphibian species, and some of these have been tested here and found to be able to bind and activate the human GLP-1 receptor. While little difference was observed for the in vitro potency for the human GLP-1 receptor, larger differences were found in the enzymatic stability of these peptides. Two peptides showed increased enzymatic stability, and they group together phylogenetically, though they originate from Amphibia and Reptilia. We have used ancestral sequence reconstruction to analyze the evolution of these GLP-1 molecules, including the synthesis of new peptides. We find that the increased stability could not be observed in the resurrected peptides from the common ancestor of frogs, even though they maintain the ability to activate the human GLP-1 receptor. Another method, using residue mapping on evolutionary branches yielded peptides that had maintained potency towards the receptor and also showed increased stability. This represents a new approach using evolutionary data in protein engineering.