Uncovering the mutation-fixation correlation in short lineages

Background  

We recently reported a highly unexpected positive correlation between the fixation probability of nonsynonymous mutations (estimated by ω) and neutral mutation rate (estimated by K _s) in mammalian lineages. However, this positive correlation was observed for lineages with relatively long divergence time such as the human-mouse lineage, and was not found for very short lineages such as the human-chimpanzee lineage. It was previously unclear how to interpret this discrepancy. It may indicate that the positive correlation between ω and K _s in long lineages is a false finding. Alternatively, it may reflect a biologically meaningful difference between various lineages. Finally, the lack of positive correlation in short lineages may be the result of methodological artifacts.     

A survey of nuclear ribosomal internal transcribed spacer substitution rates across angiosperms: an approximate molecular clock with life history effects

Background  

A full understanding of the patterns and processes of biological diversification requires the dating of evolutionary events, yet the fossil record is inadequate for most lineages under study. Alternatively, a molecular clock approach, in which DNA or amino acid substitution rates are calibrated with fossils or geological/climatic events, can provide indirect estimates of clade ages and diversification rates. The utility of this approach depends on the rate constancy of molecular evolution at a genetic locus across time and across lineages. Although the nuclear ribosomal internal transcribed spacer region (nrITS) is increasingly being used to infer clade ages in plants, little is known about the sources or magnitude of variation in its substitution rate. Here, we systematically review the literature to assess substitution rate variation in nrITS among angiosperms, and we evaluate possible correlates of the variation.     

Correlated evolution of nucleotide substitution rates and allelic variation in Mhc-DRB lineages of primates

Background  

The major histocompatibility complex (MHC) is a key model of genetic polymorphism. Selection pressure by pathogens or other microevolutionary forces may result in a high rate of non-synonymous substitutions at the codons specifying the contact residues of the antigen binding sites (ABS), and the maintenance of extreme MHC allelic variation at the population/species level. Therefore, selection forces favouring MHC variability for any reason should cause a correlated evolution between substitution rates and allelic polymorphism. To investigate this prediction, we characterised nucleotide substitution rates and allelic polymorphism (i.e. the number of alleles detected in relation to the number of animals screened) of several Mhc class II DRB lineages in 46 primate species, and tested for a correlation between them.     

The influence of body size and net diversification rate on molecular evolution during the radiation of animal phyla

Background  

Molecular clock dates, which place the origin of animal phyla deep in the Precambrian, have been used to reject the hypothesis of a rapid evolutionary radiation of animal phyla supported by the fossil record. One possible explanation of the discrepancy is the potential for fast substitution rates early in the metazoan radiation. However, concerted rate variation, occurring simultaneously in multiple lineages, cannot be detected by "clock tests", and so another way to explore such variation is to look for correlated changes between rates and other biological factors. Here we investigate two possible causes of fast early rates: change in average body size or diversification rate of deep metazoan lineages.     

Genomic divergences among cattle,dog and human estimated from large-scale alignments of genomic sequences

Background

Approximately 11 Mb of finished high quality genomic sequences were sampled from cattle, dog and human to estimate genomic divergences and their regional variation among these lineages.

Results

Optimal three-way multi-species global sequence alignments for 84 cattle clones or loci (each >50 kb of genomic sequence) were constructed using the human and dog genome assemblies as references. Genomic divergences and substitution rates were examined for each clone and for various sequence classes under different functional constraints. Analysis of these alignments revealed that the overall genomic divergences are relatively constant (0.32–0.37 change/site) for pairwise comparisons among cattle, dog and human; however substitution rates vary across genomic regions and among different sequence classes. A neutral mutation rate (2.0–2.2 × 10(-9) change/site/year) was derived from ancestral repetitive sequences, whereas the substitution rate in coding sequences (1.1 × 10(-9) change/site/year) was approximately half of the overall rate (1.9–2.0 × 10(-9) change/site/year). Relative rate tests also indicated that cattle have a significantly faster rate of substitution as compared to dog and that this difference is about 6%.

Conclusion

This analysis provides a large-scale and unbiased assessment of genomic divergences and regional variation of substitution rates among cattle, dog and human. It is expected that these data will serve as a baseline for future mammalian molecular evolution studies.

     

INFERRING THE RATES OF BRANCHING AND EXTINCTION FROM MOLECULAR PHYLOGENIES

Molecular techniques provide ancestral phylogenies of extant taxa with estimated branching times. Here we studied the pattern of ancestral phylogeny of extant taxa produced by branching (or cladogenesis) and extinction of taxa, assuming branching processes with time-dependent rates. (1) If the branching rate b and extinction rate c are constant, the semilog plot of the number of ancestral lineages over time is not a straight line but is curvilinear, with increasing slope toward the end, implying that ancestral phylogeny shows apparent increase in the branching rate near the present. The estimate of b and c based on nonlinear fitting is examined by computer simulation. The estimate of branching rate can be usable for a large phylogeny if b is greater than c, but the estimate of extinction rate c is unreliable because of large bias and variance. (2) Gradual decrease in the slope of the semilog plot of the number of ancestral lineages over time, as was observed in a phylogeny of bird families based on DNA hybridization data, can be explained equally well by either the decreasing branching rate or the increasing extinction rate. Infinitely many pairs of branching and extinction rates as functions of time can produce the same ancestral phylogeny. (3) An explosive branching event in the past would appear as a quick increase in the number of ancestral lineages. In contrast, mass extinction occurring in a brief period, if not accompanied by an increase in branching rate, does not produce any rapid change in the number of ancestral lineages at the time. (4) The condition in which the number of ancestral lineages of extant species changes in parallel with the actual number of species in the past is derived.     

Rooting human parechovirus evolution in time

Background  

The Picornaviridae family contains a number of important pathogenic viruses, among which the recently reclassified human parechoviruses (HPeVs). These viruses are widespread and can be grouped in several types. Understanding the evolutionary history of HPeV could answer questions such as how long the circulating lineages last shared a common ancestor and how the evolution of this viral species is shaped by its population dynamics. Using both strict and relaxed clock Bayesian phylogenetics we investigated 1) the substitutions rates of the structural P1 and capsid VP1 regions and 2) evolutionary timescale of currently circulating HPeV lineages.     

Post-duplication charge evolution of phosphoglucose isomerases in teleost fishes through weak selection on many amino acid sites

Background  

The partitioning of ancestral functions among duplicated genes by neutral evolution, or subfunctionalization, has been considered the primary process for the evolution of novel proteins (neofunctionalization). Nonetheless, how a subfunctionalized protein can evolve into a more adaptive protein is poorly understood, mainly due to the limitations of current analytical methods, which can detect only strong selection for amino acid substitutions involved in adaptive molecular evolution. In this study, we employed a comparative evolutionary approach to this question, focusing on differences in the structural properties of a protein, specifically the electric charge, encoded by fish-specific duplicated phosphoglucose isomerase (Pgi) genes.     

Birth-death prior on phylogeny and speed dating

Background  

In recent years there has been a trend of leaving the strict molecular clock in order to infer dating of speciations and other evolutionary events. Explicit modeling of substitution rates and divergence times makes formulation of informative prior distributions for branch lengths possible. Models with birth-death priors on tree branching and auto-correlated or iid substitution rates among lineages have been proposed, enabling simultaneous inference of substitution rates and divergence times. This problem has, however, mainly been analysed in the Markov chain Monte Carlo (MCMC) framework, an approach requiring computation times of hours or days when applied to large phylogenies.     

Adaptive evolution of chloroplast genome structure inferred using a parametric bootstrap approach

Background  

Genome rearrangements influence gene order and configuration of gene clusters in all genomes. Most land plant chloroplast DNAs (cpDNAs) share a highly conserved gene content and with notable exceptions, a largely co-linear gene order. Conserved gene orders may reflect a slow intrinsic rate of neutral chromosomal rearrangements, or selective constraint. It is unknown to what extent observed changes in gene order are random or adaptive. We investigate the influence of natural selection on gene order in association with increased rate of chromosomal rearrangement. We use a novel parametric bootstrap approach to test if directional selection is responsible for the clustering of functionally related genes observed in the highly rearranged chloroplast genome of the unicellular green alga Chlamydomonas reinhardtii, relative to ancestral chloroplast genomes.     

Mutational meltdown or controlled chain reaction: The dynamics of rapid plastome evolution in the hyperdiversity of Poaceae

The study of genomic structural evolution associated with accelerated evolutionary rates that result in avoidance of meltdown and increase biodiversity is becoming ever more possible as the number of available plastomes increases. To more comprehensively analyze rate heterogeneity among monocots and within Poaceae, we sequenced plastomes from four Poaceae species, combined them with publicly available data from ~200 plastomes, and conducted comparative analyses to quantify the pattern of rate heterogeneity between different lineages, functional groups, and periods of evolutionary time. We compared structural differences across the Poaceae to quantify how changes in plastome size correspond to different genomic subunits and the evolution of IR–SC junction boundaries. The substitution rates among ancestral Poaceae were inferred to be exceptionally rapid compared to other monocots but slowed after divergence into extant lineages, which could not be sufficiently explained by positive selection. As such, rapid rates in the ancestral lineage leading to Poaceae might be more closely linked to large-scale structural changes like the loss of ycf1 and ycf2. The total increase in plastome size across Poaceae was positively correlated with the total length of intergenic spacers, tandem repeats, and dispersed repeats as well as large single copy, and inverted repeats (IRs). The continuous evolution of IR–SC junction boundaries was asynchronous with sizes of total genome and subunits across Poaceae. Future work is needed to better understand what factors in ancestral Poaceae evolved to harness such rapid rates of plastome evolution, avoid a mutational meltdown, and escape the stagnation of strong purifying selection as well as if these factors could be utilized to synthetically control rates.     

Inferring clocks when lacking rocks: the variable rates of molecular evolution in bacteria

Background  

Because bacteria do not have a robust fossil record, attempts to infer the timing of events in their evolutionary history requires comparisons of molecular sequences. This use of molecular clocks is based on the assumptions that substitution rates for homologous genes or sites are fairly constant through time and across taxa. Violation of these conditions can lead to erroneous inferences and result in estimates that are off by orders of magnitude. In this study, we examine the consistency of substitution rates among a set of conserved genes in diverse bacterial lineages, and address the questions regarding the validity of molecular dating.     

Non-homogeneous models of sequence evolution in the Bio++ suite of libraries and programs

Background  

Accurately modeling the sequence substitution process is required for the correct estimation of evolutionary parameters, be they phylogenetic relationships, substitution rates or ancestral states; it is also crucial to simulate realistic data sets. Such simulation procedures are needed to estimate the null-distribution of complex statistics, an approach referred to as parametric bootstrapping, and are also used to test the quality of phylogenetic reconstruction programs. It has often been observed that homologous sequences can vary widely in their nucleotide or amino-acid compositions, revealing that sequence evolution has changed importantly among lineages, and may therefore be most appropriately approached through non-homogeneous models. Several programs implementing such models have been developed, but they are limited in their possibilities: only a few particular models are available for likelihood optimization, and data sets cannot be easily generated using the resulting estimated parameters.     

Natural Selection During Functional Divergence to LMP7 and Proteasome Subunit X (PSMB5) Following Gene Duplication

The LMP7 and PSMB5 genes were created through an ancient gene duplication event of their ancestral locus. These proteins contain an active site of proteolysis, and LMP7 replaces PSMB5 as a component of the 20S proteasome after stimulation of cells by interferon-. Replacement of PSMB5 by LMP7 changes the profile of the products of 20S proteasome processing, predisposing digested peptides for transport to and display by the immune system. The purpose of this study is to investigate evolutionary forces influencing functional divergence between LMP7 and PSMB5 following duplication. Levels of synonymous and nonsynonymous substitution rates are estimated to infer differences in levels of natural selection. Estimates of substitution rates indicate that natural selection elevated rates of nonsynonymous substitution in LMP7 following gene duplication, whereas PSMB5 experienced an increase in substitution rate that was not likely due to diversifying natural selection following duplication. Following initial divergence, nearly neutral mutations have dominated gene evolution in both lineages. The LMP7 gene locus provides a rare example of a protein with specialized function arising from duplication and divergence of a housekeeping protein by way of natural selection.Reviewing Editor: Dr. Rasmus Nielsen     

γ-MYN: a new algorithm for estimating Ka and Ks with consideration of variable substitution rates

Background  

Over the past two decades, there have been several approximate methods that adopt different mutation models and used for estimating nonsynonymous and synonymous substitution rates (Ka and Ks) based on protein-coding sequences across species or even different evolutionary lineages. Among them, MYN method (a Modified version of Yang-Nielsen method) considers three major dynamic features of evolving DNA sequences–bias in transition/transversion rate, nucleotide frequency, and unequal transitional substitution but leaves out another important feature: unequal substitution rates among different sites or nucleotide positions.     

Longitudinal survey of <Emphasis Type="Italic">Staphylococcus aureus</Emphasis> in cystic fibrosis patients using a multiple-locus variable-number of tandem-repeats analysis method

Background  

Staphylococcus aureus infection in patients with cystic fibrosis (CF) is frequent and may be due to colonization by a few pathogenic lineages. Systematic genotyping of all isolates, methicillin-susceptible S. aureus (MSSA) as well as methicillin-resistant S. aureus (MRSA) is necessary to identify such lineages and follow their evolution in patients. Multiple-locus variable-number tandem repeat analysis (MLVA/VNTR) was used to survey S. aureus clinical isolates in a French paediatric CF centre.     

Extensive variation in synonymous substitution rates in mitochondrial genes of seed plants

Background  

It has long been known that rates of synonymous substitutions are unusually low in mitochondrial genes of flowering and other land plants. Although two dramatic exceptions to this pattern have recently been reported, it is unclear how often major increases in substitution rates occur during plant mitochondrial evolution and what the overall magnitude of substitution rate variation is across plants.     

Unraveling the biochemistry and provenance of pupylation: a prokaryotic analog of ubiquitination

   

Recently Mycobacterium tuberculosis was shown to possess a novel protein modification, in which a small protein Pup is conjugated to the epsilon-amino groups of lysines in target proteins. Analogous to ubiquitin modification in eukaryotes, this remarkable modification recruits proteins for degradation via archaeal-type proteasomes found in mycobacteria and allied actinobacteria. While a mycobacterial protein named PafA was found to be required for this conjugation reaction, its biochemical mechanism has not been elucidated. Using sensitive sequence profile comparison methods we establish that the PafA family proteins are related to the γ-glutamyl-cysteine synthetase and glutamine synthetase. Hence, we predict that PafA is the Pup ligase, which catalyzes the ATP-dependent ligation of the terminal γ-carboxylate of glutamate to lysines, similar to the above enzymes. We further discovered that an ortholog of the eukaryotic PAC2 (e.g. cg2106) is often present in the vicinity of the actinobacterial Pup-proteasome gene neighborhoods and is likely to represent the ancestral proteasomal chaperone. Pup-conjugation is sporadically present outside the actinobacteria in certain lineages, such as verrucomicrobia, nitrospirae, deltaproteobacteria and planctomycetes, and in the latter two lineages it might modify membrane proteins.     

Whole-gene positive selection, elevated synonymous substitution rates, duplication, and indel evolution of the chloroplast clpP1 gene

Background

Synonymous DNA substitution rates in the plant chloroplast genome are generally relatively slow and lineage dependent. Non-synonymous rates are usually even slower due to purifying selection acting on the genes. Positive selection is expected to speed up non-synonymous substitution rates, whereas synonymous rates are expected to be unaffected. Until recently, positive selection has seldom been observed in chloroplast genes, and large-scale structural rearrangements leading to gene duplications are hitherto supposed to be rare.

Methodology/Principle Findings

We found high substitution rates in the exons of the plastid clpP1 gene in Oenothera (the Evening Primrose family) and three separate lineages in the tribe Sileneae (Caryophyllaceae, the Carnation family). Introns have been lost in some of the lineages, but where present, the intron sequences have substitution rates similar to those found in other introns of their genomes. The elevated substitution rates of clpP1 are associated with statistically significant whole-gene positive selection in three branches of the phylogeny. In two of the lineages we found multiple copies of the gene. Neighboring genes present in the duplicated fragments do not show signs of elevated substitution rates or positive selection. Although non-synonymous substitutions account for most of the increase in substitution rates, synonymous rates are also markedly elevated in some lineages. Whereas plant clpP1 genes experiencing negative (purifying) selection are characterized by having very conserved lengths, genes under positive selection often have large insertions of more or less repetitive amino acid sequence motifs.

Conclusions/Significance

We found positive selection of the clpP1 gene in various plant lineages to correlated with repeated duplication of the clpP1 gene and surrounding regions, repetitive amino acid sequences, and increase in synonymous substitution rates. The present study sheds light on the controversial issue of whether negative or positive selection is to be expected after gene duplications by providing evidence for the latter alternative. The observed increase in synonymous substitution rates in some of the lineages indicates that the detection of positive selection may be obscured under such circumstances. Future studies are required to explore the functional significance of the large inserted repeated amino acid motifs, as well as the possibility that synonymous substitution rates may be affected by positive selection.     

Recent evolution of alternative reproductive modes in the 'living fossil' <Emphasis Type="Italic">Triops cancriformis</Emphasis>

Background  

The Notostraca is a small but ancient crustacean order with a contrasting combination of a conservative morphology and a wide range of reproductive modes. The tadpole shrimp Triops cancriformis, includes bisexual – the putatively ancestral state -, androdioecious and hermaphrodite populations. As hermaphroditism and androdioecy confer a colonisation advantage, we expect the postglacial colonisation of northern Europe to have been effected by lineages with such reproductive modes. Therefore, N European populations should be composed of closely related lineages reflecting a recent range expansion. In contrast, glacial refugia in the south should contain bisexual populations with high haplotype diversity and more population structuring. To test these hypotheses, we analysed the geographic distribution of reproductive modes based on new and published sex ratio data. In addition, we investigated the European phylogeography of T. cancriformis by sequencing over a 1000 bp of mitochondrial DNA (mtDNA) in individuals from a large sample of populations of the three recognised subspecies.