Enzyme-coding genes as molecular clocks: The molecular evolution of animal alpha-amylases

Summary We constructed a cDNA library for the beetle,Tribolium castaneum. This library was screened using a cloned amylase gene fromDrosophila melanogaster as a molecular probe. Beetle amylase cDNA clones were isolated from this bank, and the nucleotide sequence was obtained for a cDNA clone with a coding capacity for 228 amino acids. Both the nucleotide sequence and predicted amino acid sequence were compared to our recent results forD. melanogaster alpha-amylases, along with published sequences for other alpha-amylases. The results show that animal alpha-amylases are highly conserved over their entire length. A borader comparison, which includes plant and microbial alpha-amylase sequences, indicates that parts of the gene are conserved between prokaryotes, plants, and animals. We discuss the potential importance of this and other enzyme-coding genes for the construction of molecular phylogenies and for the study of the general question of molecular clocks in evolution.     

Molecular clocks in reptiles: life history influences rate of molecular evolution

Life history has been implicated as a determinant of variation in rate of molecular evolution amongst vertebrate species because of a negative correlation between body size and substitution rate for many molecular data sets. Both the generality and the cause of the negative body size trend have been debated, and the validity of key studies has been questioned (particularly concerning the failure to account for phylogenetic bias). In this study, a comparative method has been used to test for an association between a range of life-history variables-such as body size, age at maturity, and clutch size-and DNA substitution rate for three genes (NADH4, cytochrome b, and c-mos). A negative relationship between body size and rate of molecular evolution was found for phylogenetically independent pairs of reptile species spanning turtles, lizards, snakes, crocodile, and tuatara. Although this study was limited by the number of comparisons for which both sequence and life-history data were available, the results suggest that a negative body size trend in rate of molecular evolution may be a general feature of reptile molecular evolution, consistent with similar studies of mammals and birds. This observation has important implications for uncovering the mechanisms of molecular evolution and warns against assuming that related lineages will share the same substitution rate (a local molecular clock) in order to date evolutionary divergences from DNA sequences.     

Flying rocks and flying clocks: disparity in fossil and molecular dates for birds

Major disparities are recognized between molecular divergence dates and fossil ages for critical nodes in the Tree of Life, but broad patterns and underlying drivers remain elusive. We harvested 458 molecular age estimates for the stem and crown divergences of 67 avian clades to explore empirical patterns between these alternate sources of temporal information. These divergence estimates were, on average, over twice the age of the oldest fossil in these clades. Mitochondrial studies yielded older ages than nuclear studies for the vast majority of clades. Unexpectedly, disparity between molecular estimates and the fossil record was higher for divergences within major clades (crown divergences) than divergences between major clades (stem divergences). Comparisons of dates from studies classed by analytical methods revealed few significant differences. Because true divergence ages can never be known with certainty, our study does not answer the question of whether fossil gaps or molecular dating error account for a greater proportion of observed disparity. However, empirical patterns observed here suggest systemic overestimates for shallow nodes in existing molecular divergence dates for birds. We discuss underlying biases that may drive these patterns.     

Dissimilar rates in molecular evolution

In this work we present an evolutionary tree based on the differences in the physico-chemical properties involved in amino acid substitutions, instead of considering, for its construction, only the number of changes between species. Phylogenetic trees were constructed from the differences in bulkiness, refractivity index, hydrophobicity, polarity and optical rotation of 9 vertebrate calcitonins. A correlation of the form y = a xb was found between the number of changes (x) and the differences in any given physico-chemical property (y). This correlation implies that the evolutionary time can not be evaluated directly from the number of changes between species.     

Dissimilar rates in molecular evolution

In this work we present an evolutionary tree based on the differences in the physico-chemical properties involved in amino acid substitutions, instead of considering, for its construction, only the number of changes between species. Phylogenetic trees were constructed from the differences in bulkiness, refractivity index, hydrophobicity, polarity and optical rotation of 9 vertebrate calcitonins. A correlation of the form y=a x^b was found between the number of changes (x) and the differences in any given physico-chemical property (y). This correlation implies that the evolutionary time can not be evaluated directly from the number of changes between species.     

A mitogenomic timescale for birds detects variable phylogenetic rates of molecular evolution and refutes the standard molecular clock

Current understanding of the diversification of birds is hindered by their incomplete fossil record and uncertainty in phylogenetic relationships and phylogenetic rates of molecular evolution. Here we performed the first comprehensive analysis of mitogenomic data of 48 vertebrates, including 35 birds, to derive a Bayesian timescale for avian evolution and to estimate rates of DNA evolution. Our approach used multiple fossil time constraints scattered throughout the phylogenetic tree and accounts for uncertainties in time constraints, branch lengths, and heterogeneity of rates of DNA evolution. We estimated that the major vertebrate lineages originated in the Permian; the 95% credible intervals of our estimated ages of the origin of archosaurs (258 MYA), the amniote-amphibian split (356 MYA), and the archosaur-lizard divergence (278 MYA) bracket estimates from the fossil record. The origin of modern orders of birds was estimated to have occurred throughout the Cretaceous beginning about 139 MYA, arguing against a cataclysmic extinction of lineages at the Cretaceous/Tertiary boundary. We identified fossils that are useful as time constraints within vertebrates. Our timescale reveals that rates of molecular evolution vary across genes and among taxa through time, thereby refuting the widely used mitogenomic or cytochrome b molecular clock in birds. Moreover, the 5-Myr divergence time assumed between 2 genera of geese (Branta and Anser) to originally calibrate the standard mitochondrial clock rate of 0.01 substitutions per site per lineage per Myr (s/s/l/Myr) in birds was shown to be underestimated by about 9.5 Myr. Phylogenetic rates in birds vary between 0.0009 and 0.012 s/s/l/Myr, indicating that many phylogenetic splits among avian taxa also have been underestimated and need to be revised. We found no support for the hypothesis that the molecular clock in birds "ticks" according to a constant rate of substitution per unit of mass-specific metabolic energy rather than per unit of time, as recently suggested. Our analysis advances knowledge of rates of DNA evolution across birds and other vertebrates and will, therefore, aid comparative biology studies that seek to infer the origin and timing of major adaptive shifts in vertebrates.     

Time-dependent rates of molecular evolution

For over half a century, it has been known that the rate of morphological evolution appears to vary with the time frame of measurement. Rates of microevolutionary change, measured between successive generations, were found to be far higher than rates of macroevolutionary change inferred from the fossil record. More recently, it has been suggested that rates of molecular evolution are also time dependent, with the estimated rate depending on the timescale of measurement. This followed surprising observations that estimates of mutation rates, obtained in studies of pedigrees and laboratory mutation-accumulation lines, exceeded long-term substitution rates by an order of magnitude or more. Although a range of studies have provided evidence for such a pattern, the hypothesis remains relatively contentious. Furthermore, there is ongoing discussion about the factors that can cause molecular rate estimates to be dependent on time. Here we present an overview of our current understanding of time-dependent rates. We provide a summary of the evidence for time-dependent rates in animals, bacteria and viruses. We review the various biological and methodological factors that can cause rates to be time dependent, including the effects of natural selection, calibration errors, model misspecification and other artefacts. We also describe the challenges in calibrating estimates of molecular rates, particularly on the intermediate timescales that are critical for an accurate characterization of time-dependent rates. This has important consequences for the use of molecular-clock methods to estimate timescales of recent evolutionary events.     

The molecular evolution of pituitary growth hormone prolactin and placental lactogen: A protein family showing variable rates of evolution

Summary Pituitary growth hormone and prolactin, together with the homologous placental hormones, comprise a family of related protein hormones. Complete or partial amino acid sequences of seven mammalian growth hormones, six mammalian prolactins and one placental lactogen are available, and have been compared. A phylogenetic tree has been constructed which describes the relationships within the family. At least two gene duplications have occurred during the evolution of these proteins. Rates of evolution in the family have been quite variable, the overall rate of evolution having been apparently fairly slow, but having increased markedly on several occasions, most notably in the evolution of human (and, on the basis of immunological relationships, probably other primate) growth hormones and rat (and possibly other rodent) prolactins.     

Recognition and molecular evolution in bacteria

One principal function of biological molecules in bacteria is to recognize other molecules. This allows cells to assemble for regulated enzymatic catalysis and the integration of biochemical pathways. Recognition is also an essential and a specific property in base pairing of DNA in the double helix. Therefore, recognition events must have been central to early self-assembly of primitive genetic material, genomes, cells, genetic recombination and especially in enzyme-substrate-product recognition events. Molecular recognition events are examined with an emphasis on their central role in early prokaryotic evolution.     

High rates of molecular evolution in hantaviruses

Hantaviruses are rodent-borne Bunyaviruses that infect the Arvicolinae, Murinae, and Sigmodontinae subfamilies of Muridae. The rate of molecular evolution in the hantaviruses has been previously estimated at approximately 10(-7) nucleotide substitutions per site, per year (substitutions/site/year), based on the assumption of codivergence and hence shared divergence times with their rodent hosts. If substantiated, this would make the hantaviruses among the slowest evolving of all RNA viruses. However, as hantaviruses replicate with an RNA-dependent RNA polymerase, with error rates in the region of one mutation per genome replication, this low rate of nucleotide substitution is anomalous. Here, we use a Bayesian coalescent approach to estimate the rate of nucleotide substitution from serially sampled gene sequence data for hantaviruses known to infect each of the 3 rodent subfamilies: Araraquara virus (Sigmodontinae), Dobrava virus (Murinae), Puumala virus (Arvicolinae), and Tula virus (Arvicolinae). Our results reveal that hantaviruses exhibit short-term substitution rates of 10(-2) to 10(-4) substitutions/site/year and so are within the range exhibited by other RNA viruses. The disparity between this substitution rate and that estimated assuming rodent-hantavirus codivergence suggests that the codivergence hypothesis may need to be reevaluated.     

Molecular clocks: four decades of evolution

During the past four decades, the molecular-clock hypothesis has provided an invaluable tool for building evolutionary timescales, and has served as a null model for testing evolutionary and mutation rates in different species. Molecular clocks have also influenced the development of theories of molecular evolution. As DNA-sequencing technologies have progressed, the use of molecular clocks has increased, with a profound effect on our understanding of the temporal diversification of species and genomes.     

Bats, clocks, and rocks: diversification patterns in Chiroptera

Identifying nonrandom clade diversification is a critical first step toward understanding the evolutionary processes underlying any radiation and how best to preserve future phylogenetic diversity. However, differences in diversification rates have not been quantitatively assessed for the majority of groups because of the lack of necessary analytical tools (e.g., complete species-level phylogenies, estimates of divergence times, and robust statistics which incorporate phylogenetic uncertainty and test appropriate null models of clade growth). Here, for the first time, we investigate diversification rate heterogeneity in one of the largest groups studied thus far, the bats (Mammalia: Chiroptera). We use a recent, robust statistical approach (whole-tree likelihood-based relative rate tests) on complete dated species-level supertree phylogenies. As has been demonstrated previously for most other groups, among-lineage diversification rate within bats has not been constant. However, we show that bat diversification is more heterogeneous than in other mammalian clades thus far studied. The whole-tree likelihood-based relative rates tests suggest that clades within the families Phyllostomidae and Molossidae underwent a number of significant changes in relative diversification rate. There is also some evidence for rate shifts within Pteropodidae, Emballonuridae, Rhinolophidae, Hipposideridae, and Vespertilionidae, but the significance of these shifts depends on polytomy resolution within each family. Diversification rate in bats has also not been constant, with the largest diversification rate shifts occurring 30-50 million years ago, a time overlapping with the greatest number of shifts in flowering plant diversification rates.     

Variable evolutionary rates in the molecular evolution of mammalian growth hormones

In mammals pituitary growth hormone (GH) shows a slow basal rate of evolution (0.22 ± 0.03 × 10^–9 substitutions/amino acid site/year) which appears to have increased by at least 25–50-fold on two occasions, during the evolution of primates (to at least 10.8 ± 1.3 X 10^–9 substitutions/amino acid site/year) and artiodactyl ruminants (to at least 5.6 ± 1.3 X 10^–9 substitutions/amino acid site/year). That these rate increases are real, and not due to inadvertent comparison of nonorthologous genes, was established by showing that features of the GH gene sequences that are not expressed as mature hormone do not show corresponding changes in evolutionary rate. Thus, analysis of nonsynonymous substitutions in the coding sequence for the mature protein confirmed the rate increases seen in the primate and ruminant GHs, but analysis of nonsynonymous substitutions in the signal peptide sequence, synonymous substitutions in the coding sequence for signal peptide or mature protein, and 5 and 3 untranslated sequences showed no statistically significant changes in evolutionary rate. Evidence that the increases in evolutionary rate are probably due to positive selection is provided by the observation that in the cases of both ruminant and primate GHs the periods of rapid evolution were followed by a return to a slow rate similar to the basal rate seen in other mammalian GHs. Differences between the biological properties of GHs have been identified which may relate to these periods of rapid adaptive molecular evolution. On the basis of sequence data currently available (but excluding rodent GHs which show an intermediate rate, the basis of which is not clear) for most (90%) of evolutionary time mammalian GHs have been in the slow phase of evolution, with possibly most of the few amino acid substitutions that have occurred being neutral in nature. But most (80%) of the amino acid substitutions that have been introduced into GH during the course of mammalian evolution have been accepted during the rapid phases and were adaptive in nature.     

The structure of cytochromec and the rates of molecular evolution

Summary The x-ray structure analysis of ferricytochromec shows the reasons for the evolutionary conservatism of hydrophobic and aromatic side chains, lysines, and glycines, which had been observed from comparisons of amino acid sequences from over 30 species. It also shows that the negative character of one portion of the molecular surface is conserved, even though individual acidic side chains are not, and that positive charges are localized around two hydrophobic channels leading from the interior to the surface.The reason for the unusual evolutionary conservation of surface features in cytochromesc is probably the interaction of the molecule with two other large macromolecular complexes, its reductase and oxidase. This conservation of surface structure also explains the relatively slow rate of change of cytochromec sequences in comparison with the globins and enzymes of similar size.The rate of evolution of a protein is the rate of occurrence of mutations in the genome modified by the probability that a random change in amino acid sequence will be tolerable in a functioning protein. The observed rates of change in fibrinopeptides, the globins, cytochromec, and several enzymes are interpreted in terms of the proteins' biological roles.Contribution No. 4114 from the Norman W. Church Laboratory of Chemical Biology, California Institute of Technology, Pasadena, California 91109.     

Protein evolution in deep sea bacteria: an analysis of amino acids substitution rates

Background  

Abyssal microorganisms have evolved particular features that enable them to grow in their extreme habitat. Genes belonging to specific functional categories are known to be particularly susceptible to high-pressure; therefore, they should show some evidence of positive selection. To verify this hypothesis we computed the amino acid substitution rates between two deep-sea microorganisms, Photobacterium profundum SS9 and Shewanella benthica KT99, and their respective shallow water relatives.     

Molecular clocks reduce plasmid loss rates: the R1 case

Plasmids control their replication so that the replication frequency per plasmid copy responds to the number of plasmid copies per cell. High sensitivity amplification in replication response to copy number deviations generally reduces variation in copy numbers between different single cells, thereby reducing the plasmid loss rate in a cell population. However, experiments show that plasmid R1 has a gradual, insensitive replication control predicting considerable copy number variation between single cells. The critical step in R1 copy number control is regulation of synthesis of a rate-limiting cis-acting replication protein, RepA. De novo synthesis of a large number of RepA molecules is required for replication, suggesting that copy number control is exercised at multiple steps. In this theoretical kinetic study we analyse R1 multistep copy number control and show that it results in the insensitive replication response found experimentally but that it at the same time effectively prohibits the existence of only one plasmid copy in a dividing cell. In combination with the partition system of R1, this can lead to very high segregational stability. The R1 control mechanism is compared to the different multistep copy number control of plasmid ColE1 that is based on conventional sensitivity amplification. This implies that while copy number control for ColE1 efficiently corrects for fluctuations that have already occurred, R1 copy number control prevents their emergence in cells that by chance start their cycle with only one plasmid copy. We also discuss how regular, clock-like, behaviour of single plasmid copies becomes hidden in experiments probing collective properties of a population of plasmid copies because the individual copies are out of phase. The model is formulated using master equations, taking a stochastic approach to regulation, but the mathematical formalism is kept to a minimum and the model is simplified to its bare essence. This simplicity makes it possible to extend the analysis to other replicons with similar design principles.     

Application of molecular clocks in ornithology revisited

Molecular clocks have seen many applications in ornithology, but many applications are uncritical. In this commentary, I point out logical inconsistencies in many uses of clocks in avian molecular systematics. I call for greater rigor in application of molecular clocks – clocks should only be used when clocklike behavior has been tested and confirmed, and when appropriate calibrations are available. Authors and reviewers should insist on such rigor to assure that systematics is indeed scientific, and not just storytelling.     

Darwinian aspects of molecular evolution at sublinear propagation rates

Mean fitness is non-decreasing in the symmetry sector of the frequency trajectory followed in competitive replication at sublinear propagation rates (parabolic time course). This sector contains the pairwise symmetric distribution of species frequencies and its neighboring states, and represents at least half the possible states of an evolving sublinear system. States in the non-symmetry sector produce a negative rate of change in mean fitness. The heterogeneous steady state attained in a finite sublinear system is destabilized by formation of a variant with above-threshold fitness. Evolution in the post-steady-state interval elevates the fitness threshold for coexistence. Contrary to the proposition that ‘parabolic growth invariably results in the survival of all competing species’, only species with sufficient fitness to avoid subthreshold fitness survive. An erratum to this article is available at .     

Inferring differences in the distribution of reaction rates across conditions

Elucidating changes in the distribution of reaction rates in metabolic pathways under different conditions is a central challenge in systems biology. Here we present a method for inferring regulation mechanisms responsible for changes in the distribution of reaction rates across conditions from correlations in time-resolved data. A reversal of correlations between conditions reveals information about regulation mechanisms. With the use of a small in silico hypothetical network, based on only the topology and directionality of a known pathway, several regulation scenarios can be formulated. Confronting these scenarios with experimental data results in a short list of possible pathway regulation mechanisms associated with the reversal of correlations between conditions. This procedure allows for the formulation of regulation scenarios without detailed prior knowledge of kinetics and for the inference of reaction rate changes without rate information. The method was applied to experimental time-resolved metabolomics data from multiple short-term perturbation-response experiments in S. cerevisiae across aerobic and anaerobic conditions. The method's output was validated against a detailed kinetic model of glycolysis in S. cerevisiae, which showed that the method can indeed infer the correct regulation scenario.