Integration of Bayesian molecular clock methods and fossil-based soft bounds reveals early Cenozoic origin of African lacertid lizards
© Hipsley et al; licensee BioMed Central Ltd. 2009
Received: 20 November 2008
Accepted: 01 July 2009
Published: 01 July 2009
Although current molecular clock methods offer greater flexibility in modelling evolutionary events, calibration of the clock with dates from the fossil record is still problematic for many groups. Here we implement several new approaches in molecular dating to estimate the evolutionary ages of Lacertidae, an Old World family of lizards with a poor fossil record and uncertain phylogeny. Four different models of rate variation are tested in a new program for Bayesian phylogenetic analysis called TreeTime, based on a combination of mitochondrial and nuclear gene sequences. We incorporate paleontological uncertainty into divergence estimates by expressing multiple calibration dates as a range of probabilistic distributions. We also test the reliability of our proposed calibrations by exploring effects of individual priors on posterior estimates.
According to the most reliable model, as indicated by Bayes factor comparison, modern lacertids arose shortly after the K/T transition and entered Africa about 45 million years ago, with the majority of their African radiation occurring in the Eocene and Oligocene. Our findings indicate much earlier origins for these clades than previously reported, and we discuss our results in light of paleogeographic trends during the Cenozoic.
This study represents the first attempt to estimate evolutionary ages of a specific group of reptiles exhibiting uncertain phylogenetic relationships, molecular rate variation and a poor fossil record. Our results emphasize the sensitivity of molecular divergence dates to fossil calibrations, and support the use of combined molecular data sets and multiple, well-spaced dates from the fossil record as minimum node constraints. The bioinformatics program used here, TreeTime, is publicly available, and we recommend its use for molecular dating of taxa faced with similar challenges.
The molecular clock  has become an increasingly common tool among biologists for dating the origins of species or lineages using genetic sequence data. This is commonly done by measuring the genetic distance between two species and, assuming mutations occur at a constant rate, applying an external calibration to convert those distances into geological time (for a review of molecular clock methods, see ). Despite its widespread application, calibration of the clock using independent evidence, typically derived from the fossil record, is still problematic for many groups. While proper calibration dates for major evolutionary events like the mammal-bird or fish-tetrapod split are under constant debate [3–6], less attention has been given to smaller, less inclusive clades, which often have poor fossil records. Such discussions would be particularly useful for evolutionary biologists since it is often these clades that are the subject of more detailed investigations, e.g. in the context of biogeography or diversification. In the present study, we investigate how molecular divergences can be estimated in the absence of a good fossil record, and how fossil calibrations should be applied in such a case.
Here we combine a number of new approaches in molecular dating to assign evolutionary ages to the Old World lizard family Lacertidae (Squamata). Lacertidae, with about 280 species , is the dominant reptile group in Europe and a substantial component of the squamate reptile diversity in Africa. The family is divided into two subfamilies, the Gallotiinae and Lacertinae, with the latter group composed of two monophyletic clades, the mainly Palearctic Lacertini and Eremiadini of Africa (see  for a review of lacertid systematics). Compared to their Eurasian sister taxa, the African radiation shows extraordinary taxic diversity in desert habitats, while mesic-adapted genera in Africa are relatively species-poor. This disparity in species richness is surprising given that desert lacertids are considered evolutionarily younger and therefore have had less time to speciate than their mesic sister taxa, suggesting increased speciation rates in xeric habitat. Testing this hypothesis, however, has been difficult due to our lack of knowledge on the relative antiquity of desert clades. Terrestrial squamate fossils from the Cenozoic of Africa are rare , and this lack of fossil material has seriously hampered our ability to date the main lacertid lineages – a key step towards uncovering the ecological and evolutionary factors shaping their unique biogeographic patterns.
Based on previous molecular clock estimates, lacertids may have entered Africa (and at the same time split from the Palearctic clade) after Eurasia contacted Africa in the Neogene, some 17–19 million years ago (Mya) . Fluctuating climatic conditions and aridification during that time may have promoted speciation in African lacertids through ecological displacement. According to Arnold , competitive interactions among species in mesic habitats forced subordinate taxa into drier, heterogeneous areas, resulting in niche divergence and diversification. However, a fundamental problem with this hypothesis is that other molecular divergence studies [e.g. [10, 11]] have estimated a much older age for crown (and African) lacertids, pushing their origin far into the early Cenozoic. Furthermore, dates given by Arnold et al.  are largely based on Carranza et al. , which in contrast to other molecular clock studies relies on only a single calibration point for their estimates (the age of the Canary island El Hierro to calibrate the node between Gallotia caesaris caesaris and G. c. gomerensis). In addition, both of the above studies rely on the method of nonparametric rate smoothing  which may not properly account for rate variation as it has a tendency to overfit data, particularly for regions of the tree with short branches . Therefore, thoroughly performed divergence estimates for Lacertidae, particularly for the African radiation, are still needed.
In this article, we estimate evolutionary relationships and divergence dates for the major lineages of Lacertidae with the goal of forming biogeographical hypotheses for their origin and subsequent spread throughout Africa. We construct a molecular phylogeny for the family using published nuclear and mitochondrial gene sequences, which for the first time are combined in a total evidence approach. We use multiple, well-spaced dates from fossil taxa within and outside of the family as independent calibrations. To account for uncertainty in paleontological dates we use flexible priors, meaning that calibrations are expressed as probabilistic distributions with minimum and maximum bounds [15, 16]. The use of "soft" bounds is advantageous over simple point calibrations, as potential errors in fossil dating and identification, as well as the lag time between speciation and appearance of a fossil descendent, are statistically incorporated into the prior distribution [15–17]. Additionally, we test the reliability of our proposed calibrations by excluding individual priors and evaluating posterior estimates.
All molecular clock analyses are performed in a newly available software application for Bayesian analysis called TreeTime . Like MrBayes  and BEAST , TreeTime uses a Metropolis-Coupled Markov Chain Monte Carlo (MCMCMC) method for Bayesian phylogenetic sampling. TreeTime simultaneously estimates tree topology and diversification dates and therefore does not require a starting tree topology, making it particularly appropriate for groups with uncertain phylogenies. Prior information on tree topology can be input by specifying two taxa, A and B, so that only trees in which at least one branch separates A from B are permitted. The user can also specify differently distributed priors for the time of the split between A and B. Within TreeTime, we implement four different models of rate variation and compare their performance using Bayes factor analysis. Phylogenetic relationships among the different genera are compared to previous studies, and age estimates from the different models are evaluated against available data from geology, climatology and the fossil record. Finally, we use our results as a platform for evaluating alternative hypotheses for the origins of Lacertidae and interpret our findings in light of paleogeographic trends in the Cenozoic.
Taxon sampling and alignment
GenBank accession numbers for mitochondrial and nuclear gene sequences used in the phylogenetic analysis of Lacertidae.
Alignments were performed separately for each gene using ClustalW  and manually corrected in SEAVIEW . A total of 15–20 base pairs (bp) of 16S that could not be aligned unambiguously were excluded from the analysis. Final gene lengths are 254 bp 16S, 327 bp 12S, 281 bp Cytb, 1012 bp Rag-1 and 375 bp C-mos. To test for incongruence among genes, a partition homogeneity test  was conducted in PAUP* 4.0b10 . The test (100 replicates of random addition heuristic search option with tree-bisection-reconnection branch swapping) indicated significant heterogeneity among genes (p = 0.01). However, since a growing number of studies indicate that incongruence tests are not reliable indicators of data set combinability  and no strongly supported nodes were in conflict with previous studies, genes were concatenated into a multigene data set of 2249 bp. Following a total evidence approach , the following analyses were conducted on the combined data to maximize the amount of characters and explanatory power of the available data. As a test of our combined approach, we also analyzed partitioned mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) sequences for one of the relaxed clock models (Uncorrelated lognormal with 10% prior probability distributions, described below). These values were then compared to results from the concatenated data set to explore possible biases associated with the different genomes.
Phylogenetic and molecular clock analyses
Divergence dates for Lacertidae were estimated under four different Bayesian molecular clock models. Minimum constraints for five nodes were chosen based on evidence from the fossil record. In a conservative approach, the oldest age of the stratigraphic layer in which a fossil was found was used to represent the earliest occurrence of that lineage, and potential calibrations were limited to fossils that are reliably assigned to extant clades. Calibrated nodes are: (i) Sphenodon punctatus – Cnemidophorus tigris, 228.0 Mya, based on the earliest identified rhynchocephalian from the late Triassic [Carnian; ], and corroborated by the oldest-known fossil squamate, Tikiguania, from the Carnian of India , (ii) Cnemidophorus tigris – Rhineura floridana, 113.0 Mya, corresponding to the oldest known teiid, Ptilotodon, from the lower Cretaceous [Aptian-Albian; ], (iii) Rhineura floridana – Gallotia galloti, 64.2 Mya, based on the fossil rhineurid Plesiorhineura from the Paleocene [Torrejonian; ], and (iv) Timon lepidus – Dalmatolacerta oxycephala, 5.3 Mya based on the Pliocene "Lacerta ruscinensis" from Roussillon, France, whose fossil remains are indistinguishable from the modern T. lepidus presently living in the same area .
In addition to estimating divergence dates, we evaluate the reliability of our proposed fossil calibrations by systematically removing individual priors and comparing posterior estimates. Specifically, we test the accuracy of the dates proposed for amphisbaenians (Plesiorhineura, 64.2 Mya) and teiids (Ptilotodon, 113.0 Mya) using three different treatments. In the first treatment, both the amphisbaenian and teiid are excluded so that only the oldest date (Rhynchocephalia, 228 Mya) and youngest date (Lacerta ruscinensis, 5.3 Mya) remain. In the second and third treatments, only the amphisbaenian or teiid is removed, respectively. If a calibration is accurate, provided the remaining calibrations are reliable and the data and model are appropriate, the posterior estimate should remain within the prior range even in the absence of the fossil constraint. If the calibration is poor, the posterior should move away from the prior . This approach also allows us to compare our results with other studies using similar combinations of fossil calibrations to date the origins of squamate groups [see [10, 11]].
For nucleotide data, all models are nested in the General Time Reversible model of sequence evolution with a proportion of invariant sites and gamma distributed rate heterogeneity (GTR+I+Γ), as determined by jModelTest 0.1.1 [33, 34]. For each analysis, the MCMC was run for 50,500,000 steps each chain and sampled every 500,000 steps. The first 1,000,000 steps of each run were discarded as burnin. To couple the four parallel chains we used a heating coefficient of 0.3. This resulted in a sample of size of 100 from the posterior distribution, taken from the cold chain.
MCMC calculations were performed in the program TreeTime, freely available at http://www.zi.biologie.uni-muenchen.de/evol/statgen/software/treetime. Within that program, the following models were implemented:
MC: Strict molecular clock model , assumes a fixed rate of evolution along all branches of the tree.
CPP: Compound Poisson Process , in which points of rate change are interspersed along branches. Following each substitution event, the current rate is modified according to a Poisson process with an adaptive intensity, which determines the a priori distribution of the number of changes. Rate modulations are gamma distributed, such that the expectancy value of the product of multiple rate changes is equal to 1.
ULN: Uncorrelated lognormal distributed model of Drummond et al. , in which the evolutionary rate of each branch is independently drawn from a lognormal distribution. There is no autocorrelation of rates between neighbouring branches. Parameters within the model determine the expectancy value and variance of rates. A smaller variance indicates a smaller deviation from the strict molecular clock, since rates of change are similar across branches.
DM: Dirichlet model . The a priori distribution of evolutionary rates at the branches follows a dirichlet distribution. Parameters within the model determine the variance of rates. The smaller the variance, the smaller the deviation from a strict molecular clock. The average evolutionary rate across branches is kept constant, so that only relative differences between rates are considered.
As an independent evaluation of our results, we also calculate divergence dates for Lacertidae under the ULN model in BEAST , an alternative program for Bayesian analysis. Identical model parameters were used in the two programs with the following exceptions: 1) In addition to priors for calibrated nodes, BEAST requires a prior for the distribution of divergence dates, for which we chose the Yule process. 2) BEAST estimates the equilibrium distribution of nucleotides only once at the beginning of the analysis, TreeTime samples these estimates continuously. 3) In BEAST the molecular clock is relaxed by varying molecular rates of the substitution model among branches, for which reason the rates are dependent on the time scale of the tree. TreeTime compresses or stretches the lengths of branches in the tree, given in molecular time units, by rate multipliers with a mean of one.
Finally, we test the performance of alternative clock relaxations on our data by computing Bayes factors, a Bayesian alternative to likelihood ratio tests. Bayes factors calculate the ratio of marginal likelihoods between two given models by integrating over all possible parameter values (as opposed to estimating the maximum likelihood for each parameter). In a comparison between models M1 and M2, a Bayes factor >10 on a logarithmic scale indicates that M1 is more strongly supported by the data under consideration than M2 . A significant advantage of Bayes factors over likelihood ratio tests is that they automatically penalize models with increasing complexity, and thus guard against overfitting. Furthermore, by using the strict molecular clock as a reference, they allow for a general comparison among any number of independent models .
In all phylogenetic analyses, relationships within Gallotinae were identical, however some differences in tree topology exist regarding the subfamily Lacertinae. The Eurasian subclade Lacertini is represented by a comb-like topology, where generic relationships are largely unresolved. In contrast, the African radiation is split into two monophyletic groups corresponding to Saharo-Eurasian and Ethiopian (Africa south of the Sahara Desert) distributions [e.g. [7, 39]], with the latter divided into southern and east African subclades. Phylogenetic resolution is generally higher among African genera, with the only differences between trees limited to the placement of Acanthodactylus boskianus as sister taxon to Mesalina in all cases except for the CPP 20% model, and Poromera as outgroup to the inclusive Nucras/Latastia/Heliobolus in all cases except for the strict molecular clock (MC), where it instead branches from the base of the Ethiopian clade. These slight variations in topology have no affect on relationships among the major lineages, so they are not discussed further here. Overall, all trees are essentially in agreement and in the following discussion we refer to the single tree shown in Figure 2.
Divergence dates for Lacertidae estimated from each of the molecular clock analyses are listed in Additional file 1: Lacertid clade ages. For nearly all of the relaxed clock models, the origin of modern lacertids, as indicated by the split between Gallotinae and Lacertinae, is estimated to be in the Paleocene (56–58 Mya), with the initial radiation of the African clade occurring in the mid-Eocene (44–46 Mya). Within the Eremiadini, the separation of the Saharo-Eurasian and Ethiopian clades occurred after their split from the Lacertini, 40–43 Mya. The subfamily Gallotinae diverged into its component genera, Gallotia and Psammodromus, during the Oligocene, 29–32 Mya.
Natural logarithm of Bayes factors for the molecular clock models Compound Poisson Process (CPP), Dirichlet Model (DM), Uncorrelated lognormal (ULN), and the strict Molecular Clock (MC), based on the concatenated data set.
In addition to our original divergences calculated with a 10% maximum soft bound, we expanded probability ranges to within 20% of the minimum date. Doubling prior bounds increased divergence estimates for all nodes, as well as widening confidence intervals (Additional file 1). For example, the original bounds for the Amphisbaenia-Lacertidae split were (64.2, 70.6) and the posterior estimate from to the ULN model was 68.5–83.3 Mya. When prior bounds were increased to 20%, the prior range became (64.2, 77) and the posterior estimate increased to 77.2–100.2 Mya (Additional file 1, node 3). The smallest changes resulting from this increase occurred at the Sphenodon-Squamata and Teiidae-Amphisbaenia nodes, which increased by an average of 3.2% and 7.5%, respectively. The largest change occurred at the Timon/Dalmatolacerta node, where the divergence date increased by an average of 43.6% across all models, more than double that of any other posterior expansion. Effects were most dramatic in the CPP model, which without exception produced the largest increase in divergence estimates and standard deviations when prior distributions were expanded to 20%. However, because the CPP model is unreliable for our data (see Bayes factors, Table 2), we ignore these dates in the final discussion.
Partitioned data sets
Selective deletion of calibration points
In general, our results confirm recent molecular based studies of lacertid phylogeny, including their sister relationship with amphisbaenians. Low taxon sampling within the subfamily Gallotinae hinders any phylogenetic interpretations for the group, apart from being a monophyletic clade that forms the outgroup to the remaining lacertid taxa. The Palearctic clade Lacertini forms a large polytomy in the strict consensus tree that includes the Asian grass lizard Takydromus sexlineatus (Figure 2). Relationships among the remaining Palearctic genera are poorly resolved and are therefore not discussed further here, though some biogeographical implications for the clade in general are considered below.
Several well supported sister groups occur within the African subclade Eremiadini (Figure 2). The Saharo-Eurasian group consists of the mainly North African/Asian genera Ophisops, Acanthodactylus, Eremias and Mesalina, as well as the Middle Eastern Omanosaura and equatorial African genus Adolfus. Within the Ethiopian clade, the taxon pairs Pedioplanis/Tropidosaura and Meroles/Ichnotropis form a distinct southern African cluster, while the remaining genera are distributed primarily in east Africa (with the exception of the singular species of Poromera in western central Africa). The restricted distribution of Poromera in equatorial forests may be the result of high levels of extinction some time in the past. Indeed, Africa's rainforests have diminished greatly since the early Cenozoic, and during the last 30 My a trend toward increasing aridity, coupled with repeated glacial phases, has left only small remnants of the once extensive blocks of rainforest .
Divergence estimates for Lacertidae and their evolutionary implications
Mean dates for the origin of Squamata based on the DM and ULN model (236.9, 238.2 Mya) fall well within estimates given by Vidal and Hedges  (221–251 Mya) based on nine nuclear genes, two of which are included in the present study. The split between amphisbaenians and lacertids, on the other hand, is not as well supported. Although their sister relationship is corroborated under all models, node ages vary by over 10 My in the 10% and 20% analyses, the latter case placing the split almost 30 My earlier than the earliest known rhineurid. Previous studies by Vidal and Hedges  and Wiens et al.  give much older dates for amphisbaenians, pushing their origin back to the late Jurassic-early Cretaceous. It should be noted, however, that Wiens et al.  use a different date to calibrate the Amphisbaenia-Lacertidae split based on an older fossil from the early Cretaceous (98 Mya), Hodzhakulia magna [42, 43]. This specimen consists only of incomplete maxillaries and dentaries and its purported amphisbaenian affinities have long been in doubt [44, 45], making it problematic as a calibration point.
Overall, our dates for the origin of modern lacertids are much earlier than previous estimates, placing them in the late Paleocene, 58–56 Mya. Within the Lacertidae, the majority of divergences occur in the mid- to late Eocene after the Eremiadini split from their palearctic sister clade. The separation of the African clade into its Saharo-Eurasian and Ethiopian genera occurs shortly after, and they continue to diversify until well into the mid-Miocene, some 10 Mya. The relatively young ages of the African lineages are somewhat surprising given the high levels of species richness found in desert clades. Increased rates of speciation in desert lineages may be due to selection pressures experienced in extreme environments. Adaptations to xeric habitat favoring 'r-selected' strategies (e.g. reproducing and dying quickly) could promote a shift towards shortened generation times, thus accelerating diversification [7, 46, 47]. Unfortunately, very little is known about the ecology of desert lacertids, making it difficult to determine factors underlying their biogeographic patterns. However, recent studies indicate that physiological and life history variables, such as generation time, metabolic rate, body size and clutch size, influence mutation rates in terrestrial vertebrates [48, 49], and may affect rates of molecular evolution in reptiles as well .
Historical biogeography of Lacertidae
Most authors agree that lacertids originated in Europe, as indicated by the mainly European distribution of the basal Gallotinae . According to our most reliable model (ULN), the majority of the lacertid radiation occurred in the mid-Eocene, 43–46 Mya. During that time, Europe was an archipelago of larger and smaller islands separated by shallow bodies of water . The appearance of land bridges in the Eocene as well as increasing aridity are thought to have played an important role in terrestrial vertebrate migration, and evidence for faunal exchange between Europe and Africa can be seen in the fossil records of mammals and alligators [52, 53]. A notable transition in fossil assemblages of squamate reptiles also occurs around the early Eocene in Europe, with large increases in diversity occurring at both the family and species level [ and references therein]. Unfortunately, the fossil record for African squamates during that time, particularly for small-bodied lizards, is poor  so that comparable estimates of lacertid diversity are unavailable. However, both the warming trend during the late Paleocene-early Eocene and low sea levels presumably made intercontinental dispersal feasible for a wide range of terrestrial vertebrates , possibly via land bridges or rafting.
Mayer and Benyr  and Arnold et al.  proposed the colonization of Africa by Lacertidae in the Miocene over the land bridge connecting Arabia and Africa, which remained up until the early Pliocene . Although our dates for the initial radiation of African lacertids conflict with this hypothesis, this geological event could still have played an important role in the dispersal of certain members of the Saharo-Eurasian clade. Within that group, only Adolfus and Holaspis (the latter of which is absent in our analysis) are truly African in distribution, while the remaining genera are palearctic. Distributions of Acanthodactylus, Mesalina, and Ophisops in Africa are mainly restricted to the northern Atlantic coast, and the majority of their species, along with Eremias, are found in the Middle East and Asia. This pattern suggests that the Saharo-Eurasian lineage may have originated in Eurasia and only partially left the Palearctic, as opposed to secondarily recolonizing the Middle East and Asia from Africa. In fact, the land bridge could have been crossed in the other direction, with the ancestors of modern Acanthodactylus, Mesalina, and Ophisops entering Africa from Arabia once the continents established secondary contact.
An alternative colonization scenario is that the African lineage split from the Lacertini in Europe prior to migrating to Africa, and then only later radiated into its component lineages after reaching the African continent. Discovery of a fossil lacertid in Europe with African-like qualities would support this hypothesis. Interestingly, the Baltic amber lizard Succinilacerta  from mid-Eocene Poland was for some time assigned to the south African genus Nucras [59–61], suggesting that it resembles an African lacertid, at least superficially. Unfortunately, most of the diagnostic features separating the European and African clades, including features of the clavicle, tail, ulnar nerve and hemipenis, are not externally visible in preserved specimens . Detailed investigation of this fossil, for example using X-ray Computed Tomography, could reveal internal structures assigning it to one of the modern clades. Other alternatives to fossil evidence may be found in additional taxon sampling for molecular studies. For example, inclusion of the basal African species Atlantolacerta andreanskyi in future molecular clock analyses could place a clearer temporal framework around the early evolutionary history of Eremiadini.
Reliability of fossil calibrations
To test the reliability of the oldest squamate calibrations, the Teiidae-Amphisbaenia and the Amphisbaenia-Lacertidae splits, we alternately excluded each of them and compared their respective outcomes. Our results indicate that when all calibrations are combined, the teiid calibration does not have a significant effect on the estimated divergences for lacertids, since its exclusion only marginally alters the ages of the other splits (Figure 5). At the same time, when the teiid calibration is removed, its own divergence from lacertids/amphisbaenians becomes 17.4% older than the oldest-known teiid fossil, whereas removal of the amphisbaenian calibration causes the same split to be 64% older than the prior. Although this result may support the use of fewer calibrations, it should be noted that the use of only two fossil constraints led to unrealistic estimates, pushing the origin of Squamata well into the Permian.
With respect to future studies, we suggest that it may be preferable to constrain calibration points individually depending on the quality of the fossils themselves. For some clades, the quality of different fossils in terms of stratigraphic age or reliable phylogenetic position may be highly variable, with some being easier to constrain confidently based on prior knowledge than others. In such cases, the application of qualitative phylogenetic and stratigraphic criteria as suggested by Reisz & Müller [5, 62] and Müller & Reisz  may be combined with exponential probability distributions, such that in case of a "good" fossil calibration, the soft bound spans the estimated temporal range in which the split must have occurred. Conversely, in cases of more questionable fossil dates, a 10% or 20% (or any other) upper bound may be applied.
Estimation of evolutionary ages for crown clades such as the lizard family Lacertidae may be hampered by multiple sources of uncertainty, including unknown phylogenetic relationships, lack of an adequate fossil record, and variable evolutionary rates. These are not uncommon obstacles in molecular dating, however they must still be addressed within a statistical framework. Our results highlight the advantages of a Bayesian approach. The methods we describe allow incorporation of prior information in the form of multiple fossil calibrations, while allowing for statistical flexibility and the evaluation of alternative clock models using Bayes factors. We also support the use of a total evidence approach, in which all available molecular data is combined. Particularly when implemented with multiple calibrations in a Bayesian framework, the simultaneous analysis of multiple loci provides independent constraints on the evolutionary model, thereby avoiding potential biases associated with a single gene or genome . Finally, we stress the importance of communication between paleontologists and molecular biologists in establishing suitable calibrations for more than just the major clades of Metazoa or Tetrapoda. Access to accurate information on divergence dates and paleontological material will allow biologists with diverse study systems to investigate topics such as evolutionary diversification, rates and patterns of morphological change, and historical biogeography at finer phylogenetic scales. In this regard, identifying groups needing additional study and developing plans to enable that study should be a top priority for paleontologists to position themselves as important contributors to the field of molecular dating.
Thanks to LT Lancaster (University of California, Santa Barbara) and DB Miles (Ohio University) for comments on the manuscript, and to T Liebrecht (Museum für Naturkunde, Humboldt Universität, Berlin) for help constructing figures. We are also grateful to two anonymous referees for their helpful suggestions. This work was supported by a Max Planck Institute fellowship to CAH.
- Zuckerkandl E, Pauling L: Evolutionary Divergence and Convergence in Proteins. Evolving Genes and Proteins. Edited by: Bryson V, Vogel H. 1965, New York: Academic Press, 97-166.Google Scholar
- Bromham L, Penny D: The modern molecular clock. Nature Rev Gen. 2003, 4: 216-224.View ArticleGoogle Scholar
- Smith AB, Peterson JJ: Dating the time of origin of major clades: molecular clocks and the fossil record. Annu Rev Earth Planet Sci. 2002, 30: 65-88.View ArticleGoogle Scholar
- Graur D, Martin W: Reading the entrails of chickens: molecular timescales of evolution and the illusion of precision. Trends Genet. 2004, 20: 80-86.View ArticlePubMedGoogle Scholar
- Reisz RR, Müller J: Molecular timescales and the fossil record: a paleontological perspective. Trends Genet. 2004, 20: 237-241.View ArticlePubMedGoogle Scholar
- Benton MJ, Donoghue PCJ: Paleontological evidence to date the tree of life. Mol Biol Evol. 2007, 24 (1): 26-53.View ArticlePubMedGoogle Scholar
- Arnold EN, Arribas O, Carranza S: Systematics of the Palaearctic and Oriental lizard tribe Lacertini (Squamata: Lacertidae: Lacertinae), with descriptions of eight new genera. Zootaxa. 2007, 1430: 1-86.Google Scholar
- Evans SE: At the feet of dinosaurs: the early evolutionary history and radiation of lizards. Biol Rev. 2003, 78: 513-551.View ArticlePubMedGoogle Scholar
- Arnold EN: Competition, evolutionary change and montane distributions. The Evolving Biosphere. Edited by: Forey PL. 1981, British Museum (Natural History), London, 217-228.Google Scholar
- Vidal N, Hedges SB: The phylogeny of squamate reptiles (lizards, snakes, and amphisbaenians) inferred from nine nuclear protein-coding genes. C R Biologies. 2005, 328: 1000-1008.View ArticlePubMedGoogle Scholar
- Wiens JJ, Brandley MC, Reeder TW: Why does a trait evolve multiple times within a clade? Repeated evolution of snakelike body form in squamate reptiles. Evolution. 2006, 60: 123-141.PubMedGoogle Scholar
- Carranza S, Arnold EN, Amat F: DNA phylogeny of Lacerta (Iberolacerta) and other lacertine lizards (Reptilia: Lacertidae): did competition cause long-term mountain restriction?. System Biodivers. 2004, 2 (1): 57-77.View ArticleGoogle Scholar
- Sanderson MJ: A nonparameteric approach to estimating divergence times in the absence of rate constancy. Mol Biol Evol. 1997, 14: 1218-1231.View ArticleGoogle Scholar
- Sanderson MJ: R8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003, 19: 301-302.View ArticlePubMedGoogle Scholar
- Yang Z, Rannala B: Bayesian estimation of species divergence times under a molecular clock using multiple fossil calibrations with soft bounds. Mol Biol Evol. 2006, 23: 212-226.View ArticlePubMedGoogle Scholar
- Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLOS Biol. 2006, 4 (5): e88-PubMed CentralView ArticlePubMedGoogle Scholar
- Magallón SA: Dating lineages: molecular and paleontological approaches to the temporal framework of clades. Int J Plant Sci. 2004, 165 (Suppl 4): S7-S21.View ArticleGoogle Scholar
- Himmelmann L, Metzler D: TreeTime: An extensible C++ software package for Bayesian phylogeny reconstruction with time-calibration. Bioinformatics. 2009, [http://www.zi.biologie.uni-muenchen.de/evol/statgen/software/treetime]Google Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-755.View ArticlePubMedGoogle Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-PubMed CentralView ArticlePubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignments through sequence weighting, position specific gap penalties and weight matrix choice. Nucl Acids Res. 1994, 22: 4673-4680.PubMed CentralView ArticlePubMedGoogle Scholar
- Galtier N, Gouy M, Gautier C: SeaView and Phylo_win, two graphic tools for sequence alignment and molecular phylogeny. Comput Applic Biosci. 1996, 12: 543-548.Google Scholar
- Farris JS, Källersjö M, Kluge A, Buit C: Testing significance of congruence. Cladistics. 1994, 10: 315-320.View ArticleGoogle Scholar
- Swofford DL: PAUP* Phylogeny Analysis Using Parsimony (*and other methods), version 40b10. Sinauer Association, Inc, Sunderland, Massachusetts. 2002Google Scholar
- Yoder AD, Irwin JA, Payseur BA: Failure of the ILD to determine data combinability for slow loris phylogeny. Syst Biol. 2001, 50: 408-424.View ArticlePubMedGoogle Scholar
- Kluge AG: Total evidence or taxonomic congruence: cladistics or consensus classification. Cladistics. 1998, 14 (2): 151-158.View ArticleGoogle Scholar
- Sues HD, Olsen PE: Triassic vertebrates of Gondwanan aspect from the Richmond Basin of Virginia. Science. 1990, 249: 1020-1023.View ArticlePubMedGoogle Scholar
- Datta PM, Ray S: Earliest lizard from the Late Triassic (Carnian) of India. J Vert Paleontol. 2006, 26: 795-800.View ArticleGoogle Scholar
- Nydam RL, Cifelli RL: Lizards from the Lower Cretaceous (Aptian-Albian) Antlers and Cloverly formations. J Vertebr Paleontol. 2002, 22: 286-298.View ArticleGoogle Scholar
- Sullivan RM: A new middle Paleocene (Torrejonian) rhineurid amphisbaenian, Plesiorhineura tsentasi new genus, new species, from the San Juan Basin, New Mexico. J Paleontol. 1985, 59: 1481-1485.Google Scholar
- Estes R: Handbuch der Paläoherpetologie, Part 10: Sauria terrestria, Amphisbaenia. 1983, Stuttgart: Gustav Fischer VerlagGoogle Scholar
- Sanders KL, Lee MSY: Evaluating molecular clock calibrations using Bayesian analyses with soft and hard bounds. Biol Lett. 2007, 3: 275-279.PubMed CentralView ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704.View ArticlePubMedGoogle Scholar
- Posada D: jModelTest: Phylogenetic model averaging. Mol Biol Evol. 2008, 25: 1253-1256.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Larget B, Swofford D: A compound Poisson process for relaxing the molecular clock. Genetics. 2000, 154: 1879-1892.PubMed CentralPubMedGoogle Scholar
- Himmelmann L: Bayessche Methoden zur Schätzung von Stammbäumen mit Verzweigungszeitpunkten aus molekularen Daten. PhD thesis. 2008, Goethe-University Frankfurt am Main, Germany, Department of Computer Science and MathematicsGoogle Scholar
- Kass RE, Raftery A: Bayes factors. J Am Stat Assoc. 1995, 90: 773-795.View ArticleGoogle Scholar
- Lepage T, Bryant D, Philippe H, Lartillot N: A general comparison of relaxed molecular clock models. Mol Biol Evol. 2007, 24: 2669-2680.View ArticlePubMedGoogle Scholar
- Mayer W, Pavlicev M: The phylogeny of the family Lacertidae (Reptilia) based on nuclear DNA sequences: convergent adaptations to arid habitats within the subfamily Eremiainae. Mol Phylogenet Evol. 2007, 44: 1155-1163.View ArticlePubMedGoogle Scholar
- Townsend TM, Larson A, Louis E, Macey JR: Molecular phylogenetics of Squamata: the position of snakes, amphisbaenians, and dibamids, and the root of the squamate tree. Syst Biol. 2004, 53: 735-757.View ArticlePubMedGoogle Scholar
- Wolfe JA: Distribution of major vegetation types during the Tertiary. Geophys Monogr. 1985, 32: 357-375.Google Scholar
- Nessov LA: Rare bony fishes, terrestrial lizards and mammals from the lagoonal zone of the littoral lowlands of the Cretaceous of the Kyzylkumy. Yearbook of the All-Union Palaeontological Society, Leningrad. 1985, 28: 199-219.Google Scholar
- Gao K: Sineoamphisbaena phylogenetic relationships discussed: reply. Can J Earth Sci. 1997, 34: 886-889.View ArticleGoogle Scholar
- Wu X, Russell AP, Brinkman DB: Discussion: phylogenetic relationships of Sineoamphisbaena hexatabularis: further considerations. Can J Earth Sci. 1997, 4 (6): 883-885.View ArticleGoogle Scholar
- Kearney M: Systematics of the Amphisbaenia (Lepidosauria: Squamata) based on morphological evidence from recent and fossil forms. Herpetological Monographs. 2003, 17: 1-74.View ArticleGoogle Scholar
- McKinney ML, McNamara KJ: Heterochrony the Evolution of Ontogeny. 1991, New York: Plenum PressGoogle Scholar
- Smith SA, Donoghue MJ: Rates of molecular evolution are linked to life history in flowering plants. Science. 2008, 322: 86-89.View ArticlePubMedGoogle Scholar
- Martin AP, Palumbi SR: Body size, metabolic rate, generation time, and the molecular clock. Proc Natl Acad Sci USA. 1993, 90: 4087-4091.PubMed CentralView ArticlePubMedGoogle Scholar
- Rand DM: Thermal habit, metabolic rate and the evolution of mitochondrial DNA. Trends Ecol Evol. 1994, 9: 125-131.View ArticlePubMedGoogle Scholar
- Bromham L: Molecular clocks in reptiles: life history influences rates of molecular evolution. Mol Biol Evol. 2002, 19 (3): 302-309.View ArticlePubMedGoogle Scholar
- Franzen JL: Mammalian faunal turnover in the Eocene of central Europe. Causes and Consequences of Globally Warm Climates in the Early Paleogene. Edited by: Wing SL, Gingerich D, Schmitz B, Thomas E. 2003, Colorado: Geological Society of America, 455-461.View ArticleGoogle Scholar
- Rossmann T, Müller J, Forst M: First evidence of an alligatorid from Africa and its implications for crocodilian palaeobiogeography. Neues Jahrbuch für Geologie und Paläontologie Monatshefte. 2000, 12: 651-668.Google Scholar
- Augé M: Lacertilian faunal change across the Paleocene-Eocene boundary in Europe. Causes and Consequences of Globally Warm Climates in the Early Paleogene. Edited by: Wing SL, Gingerich D, Schmitz B, Thomas E. 2003, Colorado: Geological Society of America, 441-454.View ArticleGoogle Scholar
- Popov SV, Rögl F, Rozanov AY, Steininger FF, Shcherba IG, Kovac M: Lithological-Paleogeographic maps of the Paretethys. – 10 maps Late Eocene to Pliocene. Courier Forschungsinstitut Senckenberg. 2004, 250: 1-46.Google Scholar
- Gheerbrant E, Rage J-C: Paleobiogeography of Africa: how distinct from Gondwana and Laurasia?. Palaeogeogr, Palaeoclimatol, Palaeoecol. 2006, 241: 224-246.View ArticleGoogle Scholar
- Mayer W, Benyr G: Albumin-Evolution und Phylogenese in der Familie Lacertidae. Annual Naturhistorisches Museum Wien. 1994, 96B: 621-648.Google Scholar
- Whybrow P: Geological and faunal evidence from Arabia from mammal "migrations" between Asia and Africa during the Miocene. Cour Forsch Inst Senckenberg. 1984, 69: 189-198.Google Scholar
- Böhme W, Weitschat W: Redescription of the Eocene lacertid lizard Nucras succinea Boulenger, 1917 from Baltic amber and its allocation to Succinilacerta n gen. Mitteilungen des Geologisch-Paläontologischen Instituts der Universität Hamburg. 1998, 81: 203-222.Google Scholar
- Klebs R: Über Bernsteineinschlüsse im allgemeinen und die Coleopteren meiner Bernsteinsammlung. Schr physik ökonom Ges Königsberg. 1910, 51: 217-242.Google Scholar
- Boulenger GA: A revision of the lizards of the genus Nucras, Gray. Ann S Afr Mus. 1917, 5: 196-216.Google Scholar
- Boulenger GA: Monograph of the Lacertidae 1. 1920, London: Trustees of the British Museum (Natural History), 1-352.Google Scholar
- Reisz RR, Müller J: The comparative method for evaluating fossil calibration dates: a reply to Hedges and Kumar. Trends in Genetics. 2004, 20: 596-597.View ArticlePubMedGoogle Scholar
- Müller J, Reisz RR: Four well-constrained calibration points from the vertebrate fossil record for molecular clock estimates. BioEssays. 2005, 27: 1069-1075.View ArticlePubMedGoogle Scholar
- Yang Z, Yoder AD: Comparison of likelihood and Bayesian methods for estimating divergence times using multiple gene loci and calibration points, with application to a radiation of cute-looking mouse lemur species. Syst Biol. 2003, 52: 705-716.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.