Time Clustered Sampling Can Inflate the Inferred Substitution Rate in Foot-And-Mouth Disease Virus Analyses

With the emergence of analytical software for the inference of viral evolution, a number of studies have focused on estimating important parameters such as the substitution rate and the time to the most recent common ancestor (t _MRCA) for rapidly evolving viruses. Coupled with an increasing abundance of sequence data sampled under widely different schemes, an effort to keep results consistent and comparable is needed. This study emphasizes commonly disregarded problems in the inference of evolutionary rates in viral sequence data when sampling is unevenly distributed on a temporal scale through a study of the foot-and-mouth (FMD) disease virus serotypes SAT 1 and SAT 2. Our study shows that clustered temporal sampling in phylogenetic analyses of FMD viruses will strongly bias the inferences of substitution rates and t _MRCA because the inferred rates in such data sets reflect a rate closer to the mutation rate rather than the substitution rate. Estimating evolutionary parameters from viral sequences should be performed with due consideration of the differences in short-term and longer-term evolutionary processes occurring within sets of temporally sampled viruses, and studies should carefully consider how samples are combined.     

The phylogeny of trypsin-related serine proteases and their zymogens. New methods for the investigation of distant evolutionary relationships.

The sequence of all presently known trypsin-related serine proteases and their zymogens of animal and bacterial origin were optimally aligned on the basis of three different scoring schemes for amino acid comparisons. Sequence homology was found to extend into the activation peptides. The gaps resulting from the alignment of the sequences of the active enzymes formed the basis for a new procedure based on position and number of gaps, which allowed the correct topology of the evolutionary relationship of thrombin and the pancreatic enzymes trypsin, chymotrypsin and elastase to be determined. The procedure was applied in an analogous manner to changes in disulfide bridges as well as to a selected set of amino acid positions.Evolutionary distances between proteins were estimated by minimum, base differences as well as according to the stochastic model of evolution. These distances were used successfully to find the best topology of evolutionary relationships. The fact that the branch lengths in evolutionary trees were less affected by the number of sequences considered when evolutionary distances between contemporary sequences were measured in minimum base differences than when measured according to the stochastic model of evolution, suggested in our specific case, that minimum base differences yielded estimates of evolutionary distance closer to reality than the stochastic model of evolution.All these techniques combined yielded the following picture for the evolution of the four protease families. Prothrombin and the zymogens of the pancreatic serine proteases had a common ancestor with tryptic specificity. After the initial divergence, the gene for trypsinogen duplicated. Evidence was found that the duplicated gene underwent drastic changes for a short period of time to become eventually the common ancestor of chymotrypsin and elastase. The phylogenetic tree elaborated for these enzyme families and the methods introduced to determine its topology, should readily allow determination of the attachment site of branches leading to newly sequenced serine proteases, provided their amino acid sequence can be aligned fairly unambiguously. In addition, the consequences of the alignment of the different serine proteases for the relationship of zymogen to enzyme are discussed.     

Rapid Evolution of Recombinant Saccharomyces cerevisiae for Xylose Fermentation through Formation of Extra-chromosomal Circular DNA

Circular DNA elements are involved in genome plasticity, particularly of tandem repeats. However, amplifications of DNA segments in Saccharomyces cerevisiae reported so far involve pre-existing repetitive sequences such as ribosomal DNA, Ty elements and Long Terminal Repeats (LTRs). Here, we report the generation of an eccDNA, (extrachromosomal circular DNA element) in a region without any repetitive sequences during an adaptive evolution experiment. We performed whole genome sequence comparison between an efficient D-xylose fermenting yeast strain developed by metabolic and evolutionary engineering, and its parent industrial strain. We found that the heterologous gene XylA that had been inserted close to an ARS sequence in the parent strain has been amplified about 9 fold in both alleles of the chromosomal locus of the evolved strain compared to its parent. Analysis of the amplification process during the adaptive evolution revealed formation of a XylA-carrying eccDNA, pXI2-6, followed by chromosomal integration in tandem arrays over the course of the evolutionary adaptation. Formation of the eccDNA occurred in the absence of any repetitive DNA elements, probably using a micro-homology sequence of 8 nucleotides flanking the amplified sequence. We isolated the pXI2-6 eccDNA from an intermediate strain of the evolutionary adaptation process, sequenced it completely and showed that it confers high xylose fermentation capacity when it is transferred to a new strain. In this way, we have provided clear evidence that gene amplification can occur through generation of eccDNA without the presence of flanking repetitive sequences and can serve as a rapid means of adaptation to selection pressure.     

Convergent evolution of domain architectures (is rare)

MOTIVATION: In this paper, we shall examine the evolution of domain architectures across 62 genomes of known phylogeny including all kingdoms of life. We look in particular at the possibility of convergent evolution, with a view to determining the extent to which the architectures observed in the genomes are due to functional necessity or evolutionary descent. We used domains of known structure, because from this and other information we know their evolutionary relationships. We use a range of methods including phylogenetic grouping, sequence similarity/alignment, mutation rates and comparative genomics to approach this difficult problem from several angles. RESULTS: Although we do not claim an exhaustive analysis, we conclude that between 0.4 and 4% of sequences are involved in convergent evolution of domain architectures, and expect the actual number to be close to the lower bound. We also made two incidental observations, albeit on a small sample: the events leading to convergent evolution appear to be random with no functional or structural preferences, and changes in the number of tandem repeat domains occur more readily than changes which alter the domain composition. CONCLUSION: The principal conclusion is that the observed domain architectures of the sequences in the genomes are driven by evolutionary descent rather than functional necessity. CONTACT: gough@supfam.org.     

Adaptive Evolution and the Birth of CTCF Binding Sites in the Drosophila Genome

Changes in the physical interaction between cis-regulatory DNA sequences and proteins drive the evolution of gene expression. However, it has proven difficult to accurately quantify evolutionary rates of such binding change or to estimate the relative effects of selection and drift in shaping the binding evolution. Here we examine the genome-wide binding of CTCF in four species of Drosophila separated by between ∼2.5 and 25 million years. CTCF is a highly conserved protein known to be associated with insulator sequences in the genomes of human and Drosophila. Although the binding preference for CTCF is highly conserved, we find that CTCF binding itself is highly evolutionarily dynamic and has adaptively evolved. Between species, binding divergence increased linearly with evolutionary distance, and CTCF binding profiles are diverging rapidly at the rate of 2.22% per million years (Myr). At least 89 new CTCF binding sites have originated in the Drosophila melanogaster genome since the most recent common ancestor with Drosophila simulans. Comparing these data to genome sequence data from 37 different strains of Drosophila melanogaster, we detected signatures of selection in both newly gained and evolutionarily conserved binding sites. Newly evolved CTCF binding sites show a significantly stronger signature for positive selection than older sites. Comparative gene expression profiling revealed that expression divergence of genes adjacent to CTCF binding site is significantly associated with the gain and loss of CTCF binding. Further, the birth of new genes is associated with the birth of new CTCF binding sites. Our data indicate that binding of Drosophila CTCF protein has evolved under natural selection, and CTCF binding evolution has shaped both the evolution of gene expression and genome evolution during the birth of new genes.     

K-mer natural vector and its application to the phylogenetic analysis of genetic sequences

Based on the well-known k-mer model, we propose a k-mer natural vector model for representing a genetic sequence based on the numbers and distributions of k-mers in the sequence. We show that there exists a one-to-one correspondence between a genetic sequence and its associated k-mer natural vector. The k-mer natural vector method can be easily and quickly used to perform phylogenetic analysis of genetic sequences without requiring evolutionary models or human intervention. Whole or partial genomes can be handled more effective with our proposed method. It is applied to the phylogenetic analysis of genetic sequences, and the obtaining results fully demonstrate that the k-mer natural vector method is a very powerful tool for analysing and annotating genetic sequences and determining evolutionary relationships both in terms of accuracy and efficiency.     

Convergent origins and rapid evolution of spliced leader trans-splicing in Metazoa: Insights from the Ctenophora and Hydrozoa

Replacement of mRNA 5′ UTR sequences by short sequences trans-spliced from specialized, noncoding, spliced leader (SL) RNAs is an enigmatic phenomenon, occurring in a set of distantly related animal groups including urochordates, nematodes, flatworms, and hydra, as well as in Euglenozoa and dinoflagellates. Whether SL trans-splicing has a common evolutionary origin and biological function among different organisms remains unclear. We have undertaken a systematic identification of SL exons in cDNA sequence data sets from non-bilaterian metazoan species and their closest unicellular relatives. SL exons were identified in ctenophores and in hydrozoan cnidarians, but not in other cnidarians, placozoans, or sponges, or in animal unicellular relatives. Mapping of SL absence/presence obtained from this and previous studies onto current phylogenetic trees favors an evolutionary scenario involving multiple origins for SLs during eumetazoan evolution rather than loss from a common ancestor. In both ctenophore and hydrozoan species, multiple SL sequences were identified, showing high sequence diversity. Detailed analysis of a large data set generated for the hydrozoan Clytia hemisphaerica revealed trans-splicing of given mRNAs by multiple alternative SLs. No evidence was found for a common identity of trans-spliced mRNAs between different hydrozoans. One feature found specifically to characterize SL-spliced mRNAs in hydrozoans, however, was a marked adenosine enrichment immediately 3′ of the SL acceptor splice site. Our findings of high sequence divergence and apparently indiscriminate use of SLs in hydrozoans, along with recent findings in other taxa, indicate that SL genes have evolved rapidly in parallel in diverse animal groups, with constraint on SL exon sequence evolution being apparently rare.     

Inferring pattern and process: maximum-likelihood implementation of a nonhomogeneous model of DNA sequence evolution for phylogenetic analysis

A nonhomogeneous, nonstationary stochastic model of DNA sequence evolution allowing varying equilibrium G + C contents among lineages is devised in order to deal with sequences of unequal base compositions. A maximum-likelihood implementation of this model for phylogenetic analyses allows handling of a reasonable number of sequences. The relevance of the model and the accuracy of parameter estimates are theoretically and empirically assessed, using real or simulated data sets. Overall, a significant amount of information about past evolutionary modes can be extracted from DNA sequences, suggesting that process (rates of distinct kinds of nucleotide substitutions) and pattern (the evolutionary tree) can be simultaneously inferred. G + C contents at ancestral nodes are quite accurately estimated. The new method appears to be useful for phylogenetic reconstruction when base composition varies among compared sequences. It may also be suitable for molecular evolution studies.      

A viral sampling design for testing the molecular clock and for estimating evolutionary rates and divergence times

MOTIVATION: The high pace of viral sequence change means that variation in the times at which sequences are sampled can have a profound effect both on the ability to detect trends over time in evolutionary rates and on the power to reject the Molecular Clock Hypothesis (MCH). Trends in viral evolutionary rates are of particular interest because their detection may allow connections to be established between a patient's treatment or condition and the process of evolution. Variation in sequence isolation times also impacts the uncertainty associated with estimates of divergence times and evolutionary rates. Variation in isolation times can be intentionally adjusted to increase the power of hypothesis tests and to reduce the uncertainty of evolutionary parameter estimates, but this fact has received little previous attention. RESULTS: We provide approximations for the power to reject the MCH when the alternative is that rates change in a linear fashion over time and when the alternative is that rates differ randomly among branches. In addition, we approximate the standard deviation of estimated evolutionary rates and divergence times. We illustrate how these approximations can be exploited to determine which viral sample to sequence when samples representing different dates are available.     

The economics of managing evolution

Humans are altering biological systems at unprecedented rates, and these alterations often have longer-term evolutionary impacts. Most obvious is the spread of resistance to pesticides and antibiotics. There are a wide variety of management strategies available to slow this evolution, and there are many reasons for using them. In this paper, we focus on the economic aspects of evolution management and ask: When is it economically beneficial for an individual decision-maker to invest in evolution management? We derive a simple dimensionless inequality showing that it is cost-effective to manage evolution when the percentage increase in the effective life span of the biological resource that management generates is larger than the percentage increase in annual profit that could be obtained by not managing evolution. We show how this inequality can be used to determine optimal investment choices for single decision-makers, to determine Nash equilibrium investment choices for multiple interacting decision-makers, and to examine how these equilibrium choices respond to regulatory interventions aimed at stimulating investment in evolution management. Our results are illustrated with examples involving Bacillus thuringiensis (Bt) crops and antibiotic use in fish farming.

Humans are altering biological systems at unprecedented rates and these alterations often have longer-term evolutionary impacts, such as the spread of resistance to pesticides and antibiotics. In this study, a simple mathematical criterion is derived determining when it is economically beneficial to invest in strategies for controlling and managing evolutionary change.     

Citrullination Was Introduced into Animals by Horizontal Gene Transfer from Cyanobacteria

Protein posttranslational modifications add great sophistication to biological systems. Citrullination, a key regulatory mechanism in human physiology and pathophysiology, is enigmatic from an evolutionary perspective. Although the citrullinating enzymes peptidylarginine deiminases (PADIs) are ubiquitous across vertebrates, they are absent from yeast, worms, and flies. Based on this distribution PADIs were proposed to have been horizontally transferred, but this has been contested. Here, we map the evolutionary trajectory of PADIs into the animal lineage. We present strong phylogenetic support for a clade encompassing animal and cyanobacterial PADIs that excludes fungal and other bacterial homologs. The animal and cyanobacterial PADI proteins share functionally relevant primary and tertiary synapomorphic sequences that are distinct from a second PADI type present in fungi and actinobacteria. Molecular clock calculations and sequence divergence analyses using the fossil record estimate the last common ancestor of the cyanobacterial and animal PADIs to be less than 1 billion years old. Additionally, under an assumption of vertical descent, PADI sequence change during this evolutionary time frame is anachronistically low, even when compared with products of likely endosymbiont gene transfer, mitochondrial proteins, and some of the most highly conserved sequences in life. The consilience of evidence indicates that PADIs were introduced from cyanobacteria into animals by horizontal gene transfer (HGT). The ancestral cyanobacterial PADI is enzymatically active and can citrullinate eukaryotic proteins, suggesting that the PADI HGT event introduced a new catalytic capability into the regulatory repertoire of animals. This study reveals the unusual evolution of a pleiotropic protein modification.     

A test for nucleotide sequence homology

Two macromolecular sequences which have evolved from a common ancestor sequence will tend to include a large number of elements unaffected by replacement mutations in both sequences, as long as the evolutionary rate is not too high or the divergence time is not too great. The positions of corresponding elements may have changed in either daughter sequence due to deletion/insertion mutations involving other sequence elements, but their order can be expected to be the same in both sequences. These sets of correspondences, called matches, may be computed by a recursive algorithm which incorporates constraints on the number of deletion/insertion mutations hypothesized to have occurred. A test is developed which computes the significance of each deletion/insertion hypothesized, based on Monte-Carlo sampling of random sequences with the same base composition as the experimental sequences being tested. Applying the test to 5 S RNAs confirms the relation of Escherichia coli and KB carcinoma 5 S RNAs and establishes the previously undetected homology between Pseudomonas fluorescens and KB 5 S RNAs.     

Effects of Mitochondrial DNA Rate Variation on Reconstruction of Pleistocene Demographic History in a Social Avian Species,Pomatostomus superciliosus

Mitochondrial sequence data is often used to reconstruct the demographic history of Pleistocene populations in an effort to understand how species have responded to past climate change events. However, departures from neutral equilibrium conditions can confound evolutionary inference in species with structured populations or those that have experienced periods of population expansion or decline. Selection can affect patterns of mitochondrial DNA variation and variable mutation rates among mitochondrial genes can compromise inferences drawn from single markers. We investigated the contribution of these factors to patterns of mitochondrial variation and estimates of time to most recent common ancestor (T_MRCA) for two clades in a co-operatively breeding avian species, the white-browed babbler Pomatostomus superciliosus. Both the protein-coding ND3 gene and hypervariable domain I control region sequences showed departures from neutral expectations within the superciliosus clade, and a two-fold difference in T_MRCA estimates. Bayesian phylogenetic analysis provided evidence of departure from a strict clock model of molecular evolution in domain I, leading to an over-estimation of T_MRCA for the superciliosus clade at this marker. Our results suggest mitochondrial studies that attempt to reconstruct Pleistocene demographic histories should rigorously evaluate data for departures from neutral equilibrium expectations, including variation in evolutionary rates across multiple markers. Failure to do so can lead to serious errors in the estimation of evolutionary parameters and subsequent demographic inferences concerning the role of climate as a driver of evolutionary change. These effects may be especially pronounced in species with complex social structures occupying heterogeneous environments. We propose that environmentally driven differences in social structure may explain observed differences in evolutionary rate of domain I sequences, resulting from longer than expected retention times for matriarchal lineages in the superciliosus clade.     

Timing and evolution of the most recent common ancestor of the Korean clade HIV subtype B based on <Emphasis Type="Italic">Nef</Emphasis> and <Emphasis Type="Italic">Vif</Emphasis> sequences

Molecular phylogenetic studies of the HIV-1 isolated from Koreans have suggested the presence of the so-called “Korean clade”, which can be defined as a cluster free of foreign isolates. The Korean clade accounts for more than 60% of Korean isolates and exerts characteristic amino acid sequences. Thus, it is merited to estimate when this Korean clade first emerged in order to understand the evolutionary pattern of the Korean clade. We analyzed and reconstructed the most recent common ancestor (MRCA) sequences from nef (n=229) and vif (n=179) Korean clade sequences. Linear regression analyses of sequence divergence estimates were plotted against sampling years to infer the year in which there was zero divergence from the MRCA sequences. MRCA sequences suggested the Korean clade was first emerged around 1984, before the first detection of HIV-1 in Korea in 1985. Further studies on synonymous and nonsynonymous substitution rates suggested positive selection event for the Korean clade, while other subtype B had undergone negative to neutral evolution.     

Cloning and characterization of cytoplasmic carbonic anhydrase from gills of four Antarctic fish: insights into the evolution of fish carbonic anhydrase and cold adaptation

Although carbonic anhydrase is a ubiquitous enzyme involved in a variety of physiological processes, the information on its evolution and cold adaptation among Antarctic fish is still limited: the only Antarctic fish carbonic anhydrase characterized up-to-date is from Chionodraco hamatus, a member of the Channichthyidae family. In this work, we characterized orthologous genes within two other fish families: Nototheniidae (Trematomus eulepidotus, Trematomus lepidorhinus, Trematomus bernacchii) and Bathydraconidae (Cygnodraco mawsoni). The cDNAs of epithelial gill carbonic anhydrases were cloned and sequenced. Both coding and deduced amino acid sequences were used in phylogenetic analyses. The group of enzymes preferentially expressed in fish erythrocytes (CAIIb) represented the most conserved variant. This result suggests that, although the two variants derived from the same ancestor, CAIIc genes have a more complex evolutionary history than CAIIb. The peculiar distribution of Antarctic CAs among fish CAIIcs suggests that the CAIIc gene appeared at different times through independent duplication events, even after the speciation that led to the differentiation of Antarctic fish families. Using the new CA sequences, we built homology models to trace the expected consequences of sequence variability at the protein structure level. From these analyses, we inferred that sequence variability in Antarctic fish CAs affect important physicochemical properties of these proteins and consequentially influence their reactivity. Furthermore, we searched and tested the validity of various potential molecular trademarks for cold adaptation: significant features that can be related to cold adaptation in fish CAs include reduction of positively charged solvent accessible surfaces and an increased flexibility of N-terminal and C-terminal regions.     

Sequence length bounds for resolving a deep phylogenetic divergence

In evolutionary biology, genetic sequences carry with them a trace of the underlying tree that describes their evolution from a common ancestral sequence. The question of how many sequence sites are required to recover this evolutionary relationship accurately depends on the model of sequence evolution, the substitution rate, divergence times and the method used to infer phylogenetic history. A particularly challenging problem for phylogenetic methods arises when a rapid divergence event occurred in the distant past. We analyse an idealised form of this problem in which the terminal edges of a symmetric four-taxon tree are some factor (λ) times the length of the interior edge. We determine an order λ² lower bound on the growth rate for the sequence length required to resolve the tree (independent of any particular branch length). We also show that this rate of sequence length growth can be achieved by existing methods (including the simple ‘maximum parsimony’ method), and compare these order λ² bounds with an order λ growth rate for a model that describes low-homoplasy evolution. In the final section, we provide a generic bound on the sequence length requirement for a more general class of Markov processes.