Understanding Neutral Genomic Molecular Clocks

The molecular clock hypothesis is a central concept in molecular evolution and has inspired much research into why evolutionary rates vary between and within genomes. In the age of modern comparative genomics, understanding the neutral genomic molecular clock occupies a critical place. It has been demonstrated that molecular clocks run differently between closely related species, and generation time is an important determinant of lineage specific molecular clocks. Moreover, it has been repeatedly shown that regional molecular clocks vary even within a genome, which should be taken into account when measuring evolutionary constraint of specific genomic regions. With the availability of a large amount of genomic sequence data, new insights into the patterns and causes of variation in molecular clocks are emerging. In particular, factors such as nucleotide composition, molecular origins of mutations, weak selection and recombination rates are important determinants of neutral genomic molecular clocks.     

Molecular‐clock methods for estimating evolutionary rates and timescales

The molecular clock presents a means of estimating evolutionary rates and timescales using genetic data. These estimates can lead to important insights into evolutionary processes and mechanisms, as well as providing a framework for further biological analyses. To deal with rate variation among genes and among lineages, a diverse range of molecular‐clock methods have been developed. These methods have been implemented in various software packages and differ in their statistical properties, ability to handle different models of rate variation, capacity to incorporate various forms of calibrating information and tractability for analysing large data sets. Choosing a suitable molecular‐clock model can be a challenging exercise, but a number of model‐selection techniques are available. In this review, we describe the different forms of evolutionary rate heterogeneity and explain how they can be accommodated in molecular‐clock analyses. We provide an outline of the various clock methods and models that are available, including the strict clock, local clocks, discrete clocks and relaxed clocks. Techniques for calibration and clock‐model selection are also described, along with methods for handling multilocus data sets. We conclude our review with some comments about the future of molecular clocks.     

Molecular dating in the genomic era

The comparison of DNA and protein sequences of extant species might be informative for reconstructing the chronology of evolutionary events on Earth. A phylogenetic tree inferred from molecular data directly depicts the evolutionary affinities of species and indirectly allows estimating the age of their origin and diversification. Molecular dating is achieved by assuming the molecular clock hypothesis, i.e., that the rate of change of nucleotide and amino acid sequences is on average constant over geological time. If paleontological calibrations are available, then absolute divergence times of species can be estimated. However, three major difficulties potentially hamper molecular dating : (1) a limited sample of genes and organisms, (2) a limited number of fossil references, and (3) pervasive variations of molecular evolutionary rates among genomes and species. To circumvent these problems, different solutions have been recently proposed. Larger data sets are built with more genes and more species sampled through the mining of an increasing number of genomes. Moreover, independent key fossils are identified to calibrate molecular clocks, and the uncertainty on their age is integrated in subsequent analyses. Finally, models of molecular rate variations are constructed, and incorporated in the so-called relaxed molecular clock approaches. As an illustration of these improvements, we mention that the debated age of the animal (bilaterian metazoans) diversification may have occurred between 642-761 million years ago (Mya), roughly 100 Ma before the Cambrian explosion. Among mammals, the initial diversification of major placental groups may have taken place around 100 Mya, well before the Cretaceous/Tertiary boundary marking the extinction of dinosaurs.     

The active evolutionary lives of echinoderm larvae

Echinoderms represent a researchable subset of a dynamic larval evolutionary cosmos. Evolution of echinoderm larvae has taken place over widely varying time scales from the origins of larvae of living classes in the early Palaeozoic, approximately 500 million years ago, to recent, rapid and large-scale changes that have occurred within living genera within a span of less than a million years to a few million years. It is these recent evolutionary events that offer a window into processes of larval evolution operating at a micro-evolutionary level of evolution of discrete developmental mechanisms. We review the evolution of the diverse larval forms of living echinoderms to outline the origins of echinoderm larval forms, their diversity among living echinoderms, molecular clocks and rates of larval evolution, and finally current studies on the roles of developmental regulatory mechanisms in the rapid and radical evolutionary changes observed between closely related congeneric species.     

Time zones: a comparative genetics of circadian clocks

The circadian clock is a widespread cellular mechanism that underlies diverse rhythmic functions in organisms from bacteria and fungi, to plants and animals. Intense genetic analysis during recent years has uncovered many of the components and molecular mechanisms comprising these clocks. Although autoregulatory genetic networks are a consistent feature in the design of all clocks, the weight of evidence favours their independent evolutionary origins in different kingdoms.     

Beyond Arabidopsis: The circadian clock in non-model plant species

Circadian clocks allow plants to temporally coordinate many aspects of their biology with the diurnal cycle derived from the rotation of Earth on its axis. Although there is a rich history of the study of clocks in many plant species, in recent years much progress in elucidating the architecture and function of the plant clock has emerged from studies of the model plant, Arabidopsis thaliana. There is considerable interest in extending this knowledge of the circadian clock into diverse plant species in order to address its role in topics as varied as agricultural productivity and the responses of individual species and plant communities to global climate change and environmental degradation. The analysis of circadian clocks in the green lineage provides insight into evolutionary processes in plants and throughout the eukaryotes.     

Molecular classification of living organisms

Recent studies in molecular evolution have generated strong conflicts in opinion as to how world living organisms should be classified. The traditional classification of life into five kingdoms has been challenged by the molecular analysis carried out mostly on rRNA sequences, which supported the division of the extant living organisms into three major groups: Archaebacteria, Eubacteria, and Eukaryota. As to the problem of placing the root of the tree of life, the analysis carried out on a few genes has provided discrepant results. In order to measure the genetic distances between species, we have carried out an evolutionary analysis of the glutamine synthetase genes, which previously have been revealed to be good molecular clocks, and of the small and large rRNA genes. All data demonstrate that archaebacteria are more closely related to eubacteria than to eukaryota, thus supporting the classical division of living organisms into two main superkingdoms, Prokaryota and Eukaryota.Abbreviations Mya million years ago - GS glutamine synthetase - Isu large subunit - ssu small subunit - SMC Stationary Markov clock Correspondence to: G. Pesole     

Evolutionary rates of mitochondrial genomes correspond to diversification rates and to contemporary species richness in birds and reptiles

Rates of biological diversification should ultimately correspond to rates of genome evolution. Recent studies have compared diversification rates with phylogenetic branch lengths, but incomplete phylogenies hamper such analyses for many taxa. Herein, we use pairwise comparisons of confamilial sauropsid (bird and reptile) mitochondrial DNA (mtDNA) genome sequences to estimate substitution rates. These molecular evolutionary rates are considered in light of the age and species richness of each taxonomic family, using a random-walk speciation–extinction process to estimate rates of diversification. We find the molecular clock ticks at disparate rates in different families and at different genes. For example, evolutionary rates are relatively fast in snakes and lizards, intermediate in crocodilians and slow in turtles and birds. There was also rate variation across genes, where non-synonymous substitution rates were fastest at ATP8 and slowest at CO3. Family-by-gene interactions were significant, indicating that local clocks vary substantially among sauropsids. Most importantly, we find evidence that mitochondrial genome evolutionary rates are positively correlated with speciation rates and with contemporary species richness. Nuclear sequences are poorly represented among reptiles, but the correlation between rates of molecular evolution and species diversification also extends to 18 avian nuclear genes we tested. Thus, the nuclear data buttress our mtDNA findings.     

On the molecular evolutionary clock

Summary The conceptual framework surrounding the origin of the molecular evolutionary clock and circumstances of this origin are described. In regard to the quest for the best available molecular clocks, a return to protein clocks is conditionally recommended. On the basis of recent data and certain considerations, it is pointed out that the realm of neutrality in evolution is probably less extensive than is now commonly thought, in the three distinct senses of the term neutrality—neutrality as nonfunctionality of mutations, neutrality as equifunctionality of mutations, and neutrality as a mode of fixation of mutations. The possibility is raised that complex sets of interacting components forming a system that is bounded with respect to its environment may quite generally display an intrinsic trend to a quasi-clockwise evolutionary behavior.     

Local molecular clocks in three nuclear genes: divergence times for rodents and other mammals and incompatibility among fossil calibrations

Reconstructing the chronology of mammalian evolution is a debated issue between molecule- and fossil-based inferences. A methodological limitation of molecules is the evolutionary rate variation among lineages, precluding the application of the global molecular clock. We considered 2422 first and second codon positions of the combined ADRA2B, IRBP, and vWF nuclear genes for a well-documented set of placentals including an extensive sampling of rodents. Using seven independent calibration points and a maximum-likelihood framework, we evaluated whether molecular and paleontological estimates of mammalian divergence dates may be reconciled by the local molecular clocks approach, allowing local constancy of substitution rates with variations at larger phylogenetic scales. To handle the difficulty of choosing among all possible rate assignments for various lineages, local molecular clocks were based on the results of branch-length and two-cluster tests. Extensive lineage-specific variation of evolutionary rates was detected, even among rodents. Cross-calibrations indicated some incompatibilities between divergence dates based on different paleontological references. To decrease the impact of a single calibration point, estimates derived from independent calibrations displaying only slight reciprocal incompatibility were averaged. The divergence dates inferred for the split between mice and rats (approximately 13-19 Myr) was younger than previously published molecular estimates. The most recent common ancestors of rodents, primates and rodents, boreoeutherians, and placentals were estimated to be, respectively, approximately 60, 70, 75, and 78 Myr old. Global clocks, local clocks, and quartet dating analyses suggested a Late Cretaceous origin of the crown placental clades followed by a Tertiary radiation of some placental orders like rodents.     

A distance-based least-square method for dating speciation events

Distance-based phylogenetic methods are widely used in biomedical research. However, there has been little development of rigorous statistical methods and software for dating speciation and gene duplication events by using evolutionary distances. Here we present a simple, fast and accurate dating method based on the least-squares (LS) method that has already been widely used in molecular phylogenetic reconstruction. Dating methods with a global clock or two different local clocks are presented. Single or multiple fossil calibration points can be used, and multiple data sets can be integrated in a combined analysis. Variation of the estimated divergence time is estimated by resampling methods such as bootstrapping or jackknifing. Application of the method to dating the divergence time among seven ape species or among 35 mammalian species including major mammalian orders shows that the estimated divergence time with the LS criterion is nearly identical to those obtained by the likelihood method or Bayesian inference.     

Evolution of genes and taxa: a primer

The rapidly growing fields of molecular evolution and systematics have much to offer to molecular biology, but like any field have their own repertoire of terms and concepts. Homology, for example, is a central theme in evolutionary biology whose definition is complex and often controversial. Homology extends to multigene families, where the distinction between orthology and paralogy is key. Nucleotide sequence alignment is also a homology issue, and is a key stage in any evolutionary analysis of sequence data. Models based on our understanding of the processes of nucleotide substitution are used both in the estimation of the number of evolutionary changes between aligned sequences and in phylogeny reconstruction from sequence data. The three common methods of phylogeny reconstruction – parsimony, distance and maximum likelihood – differ in their use of these models. All three face similar problems in finding optimal – and reliable – solutions among the vast number of possible trees. Moreover, even optimal trees for a given gene may not reflect the relationships of the organisms from which the gene was sampled. Knowledge of how genes evolve and at what rate is critical for understanding gene function across species or within gene families. The Neutral Theory of Molecular Evolution serves as the null model of molecular evolution and plays a central role in data analysis. Three areas in which the Neutral Theory plays a vital role are: interpreting ratios of nonsynonymous to synonymous nucleotide substitutions, assessing the reliability of molecular clocks, and providing a foundation for molecular population genetics.     

Nuclear rDNA-based molecular clock of the evolution of triatominae (Hemiptera: reduviidae), vectors of Chagas disease

The evolutionary history and times of divergence of triatomine bug lineages are estimated from molecular clocks inferred from nucleotide sequences of the small subunit SSU (18S) and the second internal transcribed spacer (ITS-2) of the nuclear ribosomal DNA of these reduviids. The 18S rDNA molecular clock rate in Triatominae, and Prosorrhynchan Hemiptera in general, appears to be of 1.8% per 100 million years (my). The ITS-2 molecular clock rate in Triatominae is estimated to be around 0.4-1% per 1 my, indicating that ITS-2 evolves 23-55 times faster than 18S rDNA. Inferred chronological data about the evolution of Triatominae fit well with current hypotheses on their evolutionary histories, but suggest reconsideration of the current taxonomy of North American species complexes.     

Improved integration time estimation of endogenous retroviruses with phylogenetic data

Background

Endogenous retroviruses (ERVs) are genetic fossils of ancient retroviral integrations that remain in the genome of many organisms. Most loci are rendered non-functional by mutations, but several intact retroviral genes are known in mammalian genomes. Some have been adopted by the host species, while the beneficial roles of others remain unclear. Besides the obvious possible immunogenic impact from transcribing intact viral genes, endogenous retroviruses have also become an interesting and useful tool to study phylogenetic relationships. The determination of the integration time of these viruses has been based upon the assumption that both 5′ and 3′ Long Terminal Repeats (LTRs) sequences are identical at the time of integration, but evolve separately afterwards. Similar approaches have been using either a constant evolutionary rate or a range of rates for these viral loci, and only single species data. Here we show the advantages of using different approaches.

Results

We show that there are strong advantages in using multiple species data and state-of-the-art phylogenetic analysis. We incorporate both simple phylogenetic information and Monte Carlo Markov Chain (MCMC) methods to date the integrations of these viruses based on a relaxed molecular clock approach over a Bayesian phylogeny model and applied them to several selected ERV sequences in primates. These methods treat each ERV locus as having a distinct evolutionary rate for each LTR, and make use of consensual speciation time intervals between primates to calibrate the relaxed molecular clocks.

Conclusions

The use of a fixed rate produces results that vary considerably with ERV family and the actual evolutionary rate of the sequence, and should be avoided whenever multi-species phylogenetic data are available. For genome-wide studies, the simple phylogenetic approach constitutes a better alternative, while still being computationally feasible.     

Molecular clocks and evolutionary relationships: Possible distortions due to horizontal gene flow

Summary This paper discusses recent evidence suggesting that genetic information from one species occasionally transfers to another remotely related species. Besides addressing the issue of whether or not the molecular data are consistent with a wide-spread influence of horizontal gene transfer, the paper shows that horizontal gene flow would not necessarily preclude a linear molecular clock or change the rate of molecular evolution (assuming the neutral allele theory). A pervasive influence of horizontal gene transfer is more than just consistent with the data of molecular evolution, it also provides a unique explanation for a number of possibly conflicting phylogenies and contradictory clocks. This phenomenon might explain why some protein clocks are linear while the superoxide dismutase clock is not, how the molecular data on the phylogeny of apes and Australian song birds are not necessarily in conflict with those based on morphology, and, finally, why the mycoplasmas have an accelerated molecular clock.     

Number of ancestral human species: a molecular perspective

Despite the remarkable developments in molecular biology over the past three decades, anthropological genetics has had only limited impact on systematics in human evolution. Genetics offers the opportunity to objectively test taxonomies based on morphology and may be used to supplement conventional approaches to hominid systematics. Our analyses, examining chromosomes and 46 estimates of genetic distance, indicate there may have been only around 4 species on the direct line to modern humans and 5 species in total. This contrasts with current taxonomies recognising up to 23 species.

The genetic proximity of humans and chimpanzees has been used to suggest these species are congeneric. Our analysis of genetic distances between them is consistent with this proposal. It is time that chimpanzees, living humans and all fossil humans be classified in Homo. The creation of new genera can no longer be a solution to the complexities of fossil morphologies. Published genetic distances between common chimpanzees and bonobos, along with evidence for interbreeding, suggest they should be assigned to a single species.

The short distance between humans and chimpanzees also places a strict limit on the number of possible evolutionary side branches that might be recognised on the human lineage. All fossil taxa were genetically very close to each other and likely to have been below congeneric genetic distances seen for many mammals.

Our estimates of genetic divergence suggest that periods of around 2 million years are required to produce sufficient genetic distance to represent speciation. Therefore, Neanderthals and so-called H. erectus were genetically so close to contemporary H. sapiens they were unlikely to have been separate species. Thus, it is likely there was only one species of human (H. sapiens) for most of the last 2 million years. We estimate the divergence time of H. sapiensfrom 16 genetic distances to be around 1.7 Ma which is consistent with evidence for the earliest migration out of Africa. These findings call into question the mitochondrial «African Eve» hypothesis based on a far more recent origin for H. sapiens and show that humans did not go through a bottleneck in their recent evolutionary history.

Given the large offset in evolutionary rates of molecules and morphology seen in human evolution, Homo species are likely to be characterised by high levels of morphological variation and low levels of genetic variability. Thus, molecular data suggest the limits for intraspecific morphological variation used by many palaeoanthropologists have been set too low. The role of phenotypic plasticity has been greatly underestimated in human evolution.

We call into question the use of mtDNA for studies of human evolution. This DNA is under strong selection, which violates the assumption of selective neutrality. This issue should be addressed by geneticists, including a reassessment of its use for molecular clocks. There is a need for greater cooperation between palaeoanthropologists and anthropological geneticists to better understand human evolution and to bring palaeoanthropology into the mainstream of evolutionary biology.     

Human and ape molecular clocks and constraints on paleontological hypotheses

Although the relationships of the living hominoid primates (humans and apes) are well known, the relationships of the fossil species, times of divergence of both living and fossil species, and the biogeographic history of hominoids are not well established. Divergence times of living species, estimated from molecular clocks, have the potential to constrain hypotheses of the relationships of fossil species. In this study, new DNA sequences from nine protein-coding nuclear genes in great apes are added to existing datasets to increase the precision of molecular time estimates bearing on the evolutionary history of apes and humans. The divergence of Old World monkeys and hominoids at the Oligocene-Miocene boundary (approximately 23 million years ago) provides the best primate calibration point and yields a time and 95% confidence interval of 5.4 +/- 1.1 million years ago (36 nuclear genes) for the human-chimpanzee divergence. Older splitting events are estimated as 6.4 +/- 1.5 million years ago (gorilla, 31 genes), 11.3 +/- 1.3 million years ago (orangutan, 33 genes), and 14.9 +/- 2.0 million years ago (gibbon, 27 genes). Based on these molecular constraints, we find that several proposed phylogenies of fossil hominoid taxa are unlikely to be correct.     

Calibration of galliform molecular clocks using multiple fossils and genetic partitions

For more than a century, members of the traditional avian order Galliformes (i.e., pheasants, partridges, junglefowl, and relatives) have been among the most intensively studied birds, but still a comprehensive timeframe for their evolutionary history is lacking. Thanks to a number of recent cladistic interpretations for several galliform fossils, candidates now exist that can potentially be used as accurate internal calibrations for molecular clocks. Here, we describe a molecular timescale for Galliformes based on cytochrome b and ND2 using nine mostly internal fossil-based anchorpoints. Beyond application of calibrations spanning the entire evolutionary history of Galliformes, care was taken to investigate the effects of calibration choice, substitution saturation, and rate heterogeneity among lineages on divergence time estimation. Results show broad consistency in time estimation with five out of the nine total calibrations. Our divergence time estimates, based on these anchorpoints, indicate that the early history of Galliformes took place in the Cretaceous, including the origin of the basal-most megapode and perhaps cracid lineages, but that the remaining morphological diversification likely started in the earliest Tertiary. The multi-calibration/multi-genetic partition approach used here highlights the importance of understanding the genetic saturation, variation, and rate constancy spectra for the accurate calculation of divergence times by use of molecular clocks.     

Primate phylogeny, evolutionary rate variations, and divergence times: a contribution from the nuclear gene IRBP

The first third (ca. 1200 bp) of exon 1 of the nuclear gene encoding the interstitial retinoid-binding protein (IRBP) has been sequenced for 12 representative primates belonging to Lemuriformes, Lorisiformes, Tarsiiformes, Platyrrhini, and Catarrhini, and combined with available data (13 other primates, 11 nonprimate placentals, and 2 marsupials). Phylogenetic analyses using maximum likelihood on nucleotides and amino acids robustly support the monophyly of primates, Strepsirrhini, Lemuriformes, Lorisiformes, Anthropoidea, Catarrhini, and Platyrrhini. It is interesting to note that 1) Tarsiidae grouped with Anthropoidea, and the support for this node depends on the molecular characters considered; 2) Cheirogaleidae grouped within Lemuriformes; and 3) Daubentonia was the sister group of all other Lemuriformes. Study of the IRBP evolutionary rate shows a high heterogeneity within placentals and also within primates. Maximum likelihood local molecular clocks were assigned to three clades displaying significantly contrasted evolutionary rates. Paenungulata were shown to evolve 2.5-3 times faster than Perissodactyla and Lemuriformes. Six independent calibration points were used to estimate splitting ages of the main primate clades, and their compatibility was evaluated. Divergence ages were obtained for the following crown groups: 13.8-14.2 MY for Lorisiformes, 26.5-27.2 MY for Lemuroidea, 39.6-40.7 MY for Lemuriformes, 45.4-46.7 MY for Strepsirrhini, and 56.7-58.4 MY for Haplorrhini. The incompatibility between some paleontological and molecular estimates may reflect the incompleteness of the placental fossil record, and/or indicate that the variable IRBP evolutionary rates are not fully accommodated by local molecular clocks.     

Explosive radiations and the reliability of molecular clocks: island endemic radiations as a test case

The reliability of molecular clocks has been questioned for several key evolutionary radiations on the basis that the clock might run fast in explosive radiations. Molecular date estimates for the radiations of metazoan phyla (the Cambrian explosion) and modern orders of mammals and birds are in many cases twice as old as the palaeontological evidence would suggest. Could some aspect of explosive radiations speed the molecular clock, making molecular date estimates too old? Here we use 19 independent instances of recent explosive radiations of island endemic taxa as a model system for testing the proposed influence of rapid adaptive radiation on the rate of molecular evolution. These radiations are often characterized by many of the potential mechanisms for fast rates in explosive radiations--such as small population size, elevated speciation rate, rapid rate of morphological change, release from previous ecological constraints, and adaptation to new niches--and represent a wide variety of species, islands, and genes. However, we find no evidence of a consistent increase in rates in island taxa compared to their mainland relatives, and therefore find no support for the hypothesis that the molecular clock runs fast in explosive radiations.