Timing the origin of New World monkeys

The origin of New World monkeys (Infraorder Platyrrhini) has been an extensively debated issue. In this study, we analyzed mitochondrial genomes from Cebus (Platyrrhini), Homo, Hylobates, Pan, Pongo (Hominoids), Macaca, Papio (Cercopithecoids), and Tarsius (outgroup) to investigate this matter. Two distinct methodologies were employed on mitochondrial genes to estimate divergence times: the traditional likelihood ratio test performed in ML analyses of individual and concatenated gene sequences and the recent multigene Bayesian approach. Using the Cercopithecoid-Hominoid split as calibration point (25 MYA), our results show consistently that Platyrrhines split from Catarrhines at around 35 MYA. Although the main focus of the study is New World monkey origins, we have also estimated other primate divergence times: Homo-Pan at 5-7 MYA; Pongo-(Homo/Pan) at 13-16 MYA; Hylobates-(Pongo/Homo/Pan) at 15-19 MYA; and Macaca-Papio at 10-12 MYA. Our estimate for the origin of New World monkeys is in agreement with the hypothesis of a transatlantic journey from Africa to South America, as suggested by the fossil record.     

Estimation of primate speciation dates using local molecular clocks

Protein-coding genes of the mitochondrial genomes from 31 mammalian species were analyzed to estimate the speciation dates within primates and also between rats and mice. Three calibration points were used based on paleontological data: one at 20-25 MYA for the hominoid/cercopithecoid divergence, one at 53-57 MYA for the cetacean/artiodactyl divergence, and the third at 110-130 MYA for the metatherian/eutherian divergence. Both the nucleotide and the amino acid sequences were analyzed, producing conflicting results. The global molecular clock was clearly violated for both the nucleotide and the amino acid data. Models of local clocks were implemented using maximum likelihood, allowing different evolutionary rates for some lineages while assuming rate constancy in others. Surprisingly, the highly divergent third codon positions appeared to contain phylogenetic information and produced more sensible estimates of primate divergence dates than did the amino acid sequences. Estimated dates varied considerably depending on the data type, the calibration point, and the substitution model but differed little among the four tree topologies used. We conclude that the calibration derived from the primate fossil record is too recent to be reliable; we also point out a number of problems in date estimation when the molecular clock does not hold. Despite these obstacles, we derived estimates of primate divergence dates that were well supported by the data and were generally consistent with the paleontological record. Estimation of the mouse-rat divergence date, however, was problematic.     

The phylogeny of the hominoid primates,as indicated by DNA-DNA hybridization

Summary The living hominoid primates are Man, the chimpanzees, the Gorilla, the Orangutan, and the gibbons. The cercopithecoids (Old World monkeys) are the sister group of the hominoids. The composition of the Hominoidea is not in dispute, but a consensus has not yet been reached concerning the phylogenetic branching pattern and the dating of divergence nodes. We have compared the single-copy nuclear DNA sequences of the hominoid genera using DNA-DNA hybridization to produce a complete matrix of delta T₅₀H values. The data show that the branching sequence of the lineages, from oldest to most recent, was: Old World monkeys, gibbons, Orangutan, Gorilla, chimpanzees, and Man. The calibration of the delta T₅₀H scale in absolute time needs further refinement, but the ranges of our estimates of the datings of the divergence nodes are: Cercopithecoidea, 27–33 million years ago (MYA); gibbons, 18–22 MYA; Orangutan, 13–16 MYA; Gorilla, 8–10 MYA; and chimpanzees-Man, 6.3–7.7 MYA.     

Fossil calibration of molecular divergence infers a moderate mutation rate and recent radiations for pinus

Silent mutation rate estimates for Pinus vary 50-fold, ranging from angiosperm-like to among the slowest reported for plants. These differences either reflect extraordinary genomic processes or inconsistent fossil calibration, and they have important consequences for population and biogeographical inferences. Here we estimate mutation rates from 4 Pinus species that represent the major lineages using 11 nuclear and 4 chloroplast loci. Calibration was tested at the divergence of Pinus subgenera with the oldest leaf fossil from subg. Strobus (Eocene; 45 MYA) or a recently published subg. Strobus wood fossil (Cretaceous; 85 MYA). These calibrations place the origin of Pinus 190-102 MYA and give absolute silent rate estimates of 0.70-1.31x10(-9) and 0.22-0.42x10(-9).site-1.year-1 for the nuclear and chloroplast genomes, respectively. These rates are approximately 4- to 20-fold slower than angiosperms, but unlike many previous estimates, they are more consistent with the high per-generation deleterious mutation rates observed in pines. Chronograms from nuclear and chloroplast genomes show that the divergence of subgenera accounts for about half of the time since Pinus diverged from Picea, with subsequent radiations occurring more recently. By extending the sampling to encompass the phylogenetic diversity of Pinus, we predict that most extant subsections diverged during the Miocene. Moreover, subsect. Australes, Ponderosae, and Contortae, containing over 50 extant species, radiated within a 5 Myr time span starting as recently as 18 MYA. An Eocene divergence of pine subgenera (using leaf fossils) does not conflict with fossil-based estimates of the Pinus-Picea split, but a Cretaceous divergence using wood fossils accommodates Oligocene fossils that may represent modern subsections. Because homoplasy and polarity of character states have not been tested for fossil pine assignments, the choice of fossil and calibration node represents a significant source of uncertainty. Based on several lines of evidence (including agreement with ages inferred using calibrations outside of Pinus), we conclude that the 85 MYA calibration at the divergence of pine subgenera provides a reasonable lower bound and that further refinements in age and mutation rate estimates will require a synthetic examination of pine fossil history.     

The mitochondrial genome of spotted green pufferfish Tetraodon nigroviridis (Teleostei: Tetraodontiformes) and divergence time estimation among model organisms in fishes

We determined the whole mitochondrial genome sequence for spotted green pufferfish, Tetraodon nigroviridis (Teleostei: Tetraodontiformes). The genome (16,488 bp) contained 37 genes (two ribosomal RNA genes, 22 transfer RNA genes, and 13 protein-coding genes) plus control region as found in other vertebrates, with the gene order identical to that of typical vertebrates. The sequence was used to estimate phylogenetic relationships and divergence times among major lineages of fishes, including representative model organisms in fishes. We employed partitioned Bayesian approaches for these two analyses using two datasets that comprised concatenated amino acid sequences from 12 protein-coding genes (excluding the ND6 gene) and concatenated nucleotide sequences from the 12 protein-coding genes (without 3rd codon positions), 22 transfer RNA genes, and two ribosomal RNA genes. The resultant trees from the two datasets were well resolved and largely congruent with those from previous studies, with spotted green pufferfish being placed in a reasonable phylogenetic position. The approximate divergence times between spotted green pufferfish and model organisms in fishes were 85 million years ago (MYA) vs. torafugu, 183 MYA vs. three-spined stickleback, 191 MYA vs. medaka, and 324 MYA vs. zebrafish, all of which were about twice as old as the divergence times estimated by their earliest occurrences in fossil records.     

Pleistocene diversification of living squirrel monkeys (Saimiri spp.) inferred from complete mitochondrial genome sequences

In order to enhance our understanding of the evolutionary history of squirrel monkeys (Saimiri spp.), we newly sequenced and analyzed data from seven complete mitochondrial genomes representing six squirrel monkey taxa. While previous studies have lent insights into the taxonomy and phylogeny of the genus, phylogenetic relationships and divergence date estimates among major squirrel monkey clades remain unclear. Using maximum likelihood and Bayesian procedures, we inferred a highly resolved phylogenetic tree with strong support for a sister relationship between Saimiri boliviensis and all other Saimiri, for monophyly of Saimiri oerstedii and Saimiri sciureus sciureus, and for Saimiri sciureus macrodon as the sister lineage to the S. oerstedii/S. s. sciureus clade. We inferred that crown lineages for extant squirrel monkeys diverged around 1.5 million years ago (MYA) in the Pleistocene Epoch, with other major clades diverging between 0.9 and 1.1 MYA. Our results suggest a relatively recent timeline of squirrel monkey evolution and challenge previous conceptions about the diversification of the genus and its expansion into Central America.     

Slow molecular clocks in Old World monkeys,apes, and humans

Two longstanding issues on the molecular clock hypothesis are studied in this article. First, is there a global molecular clock in mammals? Although many authors have observed unequal rates of nucleotide substitution among mammalian lineages, some authors have proposed a global clock for all eutherians, i.e., a single global rate of 2.2 x 10(-9) substitutions per nucleotide site per year. We reexamine this issue using noncoding, nonrepetitive DNA from Old World monkeys (OWMs), chimpanzee, and human. First, using the minimal date of 6 MYA for the human-chimpanzee divergence and more than 2.5 million base pairs of genomic sequences from human and chimpanzee, we estimate a maximal rate of 0.99 x 10(-9) for noncoding, nonrepetitive genomic regions for these two species. This estimate is less than half of the proposed global rate and much smaller than the commonly used rate (3.5 x 10(-9)) for eutherians. Further, using a minimal date of 23 MYA for the human-OWM divergence, we estimate a maximal rate of 1.5 x 10(-9) for both introns and fourfold degenerate sites in humans and OWMs. In addition, with the New World monkey (NWM) lineage as an outgroup, we estimate that the rate of substitution in introns is 30% higher in the OWM lineage than in the human lineage. Clearly, there is no global molecular clock in eutherians. Second, although many studies have indicated considerable variation in the mutation rate among regions of the mammalian genome, a recent study proposed a uniform rate. Using new and existing intron sequence data from higher primates, we find significant rate variation among genomic regions and a positive correlation between the rate of substitution and the GC content, refuting the claim of a uniform rate.     

Catarrhine primate divergence dates estimated from complete mitochondrial genomes: concordance with fossil and nuclear DNA evidence

Accurate divergence date estimates improve scenarios of primate evolutionary history and aid in interpretation of the natural history of disease-causing agents. While molecule-based estimates of divergence dates of taxa within the superfamily Hominoidea (apes and humans) are common in the literature, few such estimates are available for the Cercopithecoidea (Old World monkeys), the sister taxon of the hominoids in the primate infraorder Catarrhini. To help fill this gap, we have sequenced the entire mitochondrial DNA (mtDNA) genomes from a representative of three cercopithecoid tribes, Cercopithecini (Chlorocebus aethiops), Colobini (Colobus guereza), and Presbytini (Trachypithecus obscurus), and analyzed these new data together with other catarrhine mtDNA genomes available in public databases. Molecular divergence date estimates are dependent on calibration points gleaned from the paleontological record. We defined criteria for the selection of good calibration points and identified three points meeting these criteria: Homo-Pan, 6.0 Ma; Pongo-hominines, 14.0 Ma; hominoid/cercopithecoid, 23.0 Ma. Because a uniform molecular clock does not fit the catarrhine mtDNA data, we estimated divergence dates using a penalized likelihood and a Bayesian method, both of which take into account the effects of rate differences on lineages, phylogenetic tree structure, and multiple calibration points. The penalized likelihood method applied to the coding regions of the mtDNA genome yielded the following divergence date estimates, with approximate 95% confidence intervals: cercopithecine-colobine, 16.2 (14.4-17.9) Ma; colobin-presbytin, 10.9 (9.6-12.3) Ma; cercopithecin-papionin, 11.6 (10.3-12.9) Ma; and Macaca-Papio, 9.8 (8.6-10.9) Ma. Within the hominoids, the following dates were inferred: hylobatid-hominid, 16.8 (15.0-18.5) Ma; Gorilla-Homo+Pan, 8.1 (7.1-9.0) Ma; Pongo pygmaeus pygmaeus-P. p. abelii, 4.1 (3.5-4.7) Ma; and Pan troglodytes-P. paniscus, 2.4 (2.0-2.7) Ma. These dates were similar to those found using penalized likelihood on other subsets of the data, but slightly younger than several of the Bayesian estimates.     

Parasites reveal movement of bats between the New and Old Worlds

The global distribution of bat taxa indicates that the Atlantic and Pacific Oceans are effective barriers to movement between the Old and New Worlds. For instance, one of the major suborders, Yinpterochiroptera, has an exclusively Old World distribution, and within the other, Yangochiroptera, no species and only five genera are common to both. However, as bats are sometimes blown out to sea, and have colonised isolated islands, occasional natural movement between the New and Old Worlds does appear to be possible. Here we identify new genotypes of a blood parasite, Trypanosoma dionisii, in Old World bats that are closely related to South American strains. Using highly conservative calibration points, divergence of Old and New World strains is estimated to have occurred 3.2-5.0 million years ago (MYA), depending on the method used (upper 95% CL for maximum time 11.4MYA). The true date of divergence is likely to be considerably more recent. These results demonstrate that taxon-specific parasites can indicate historical movements of their hosts, even where their hosts may have left no lasting phylogenetic footprint.     

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.     

The origin and diversification of eukaryotes: problems with molecular phylogenetics and molecular clock estimation

Determining the relationships among and divergence times for the major eukaryotic lineages remains one of the most important and controversial outstanding problems in evolutionary biology. The sequencing and phylogenetic analyses of ribosomal RNA (rRNA) genes led to the first nearly comprehensive phylogenies of eukaryotes in the late 1980s, and supported a view where cellular complexity was acquired during the divergence of extant unicellular eukaryote lineages. More recently, however, refinements in analytical methods coupled with the availability of many additional genes for phylogenetic analysis showed that much of the deep structure of early rRNA trees was artefactual. Recent phylogenetic analyses of a multiple genes and the discovery of important molecular and ultrastructural phylogenetic characters have resolved eukaryotic diversity into six major hypothetical groups. Yet relationships among these groups remain poorly understood because of saturation of sequence changes on the billion-year time-scale, possible rapid radiations of major lineages, phylogenetic artefacts and endosymbiotic or lateral gene transfer among eukaryotes. Estimating the divergence dates between the major eukaryote lineages using molecular analyses is even more difficult than phylogenetic estimation. Error in such analyses comes from a myriad of sources including: (i) calibration fossil dates, (ii) the assumed phylogenetic tree, (iii) the nucleotide or amino acid substitution model, (iv) substitution number (branch length) estimates, (v) the model of how rates of evolution change over the tree, (vi) error inherent in the time estimates for a given model and (vii) how multiple gene data are treated. By reanalysing datasets from recently published molecular clock studies, we show that when errors from these various sources are properly accounted for, the confidence intervals on inferred dates can be very large. Furthermore, estimated dates of divergence vary hugely depending on the methods used and their assumptions. Accurate dating of divergence times among the major eukaryote lineages will require a robust tree of eukaryotes, a much richer Proterozoic fossil record of microbial eukaryotes assignable to extant groups for calibration, more sophisticated relaxed molecular clock methods and many more genes sampled from the full diversity of microbial eukaryotes.     

Maximum-likelihood divergence date estimates based on rRNA gene sequences suggest two scenarios of Trypanosoma cruzi intraspecific evolution.

The phylogenetic relationships of Trypanosoma cruzi strains were inferred using maximum-likelihood from complete 18S rDNA sequences and D7-24Salpha rDNA regions from 20 representative strains of T. cruzi. For this we sequenced the 18S rDNA of 14 strains and the D7-24Salpha rDNA of four strains and aligned them to previously published sequences. Phylogenies inferred from these data sets identified four groups, named Riboclades 1, 2, 3, and 4, and a basal dichotomy that separated Riboclade 1 from Riboclades 2, 3, and 4. Substitution models and other parameters were optimized by hierarchical likelihood tests, and our analysis of the 18S rDNA molecular clock by the likelihood ratio test suggests that a taxa subset encompassing all 2,150 positions in the alignment supports rate constancy among lineages. The present analysis supports the notion that divergence dates of T. cruzi Riboclades can be estimated from 18S rDNA sequences and therefore, we present alternative evolutionary scenarios based on two different views of T. cruzi intraspecific divergence. The first assumes a faster evolutionary rate, which suggests that the divergence between T. cruzi I and II and the extant strains occurred in the Tertiary period (37-18 MYA). The other, which supports the hypothesis that the divergence between T. cruzi I and II occurred in the Cretaceous period (144-65 MYA) and the divergence of the extant strains occurred in the Tertiary period of the Cenozoic era (65-1.8 MYA), is consistent with our previously proposed hypothesis of divergence by geographical isolation and mammalian host coevolution.     

Phylogenetics of Perissodactyla and Tests of the Molecular Clock

Two mitochondrial genes, the protein-coding cytochrome c oxidase subunit II (COII) gene and a portion of the 12S rRNA gene, were used for phylogenetic investigation of the mammalian order Perissodactyla. The primary objective of the study was to utilize the extensive fossil record of perissodactyls for calibrating molecular clocks and comparing estimates of divergence times using both genes and two fossil calibration points. Secondary objectives included clarification of previously unresolved relationships within Tapiridae and comparison of the results of separate and combined analyses of two genes. Analyses included several perissodactyl lineages representing all three families (Tapiridae, Equidae, and Rhinocerotidae), most extant genera, all four species of tapirs, two to four species of rhinoceros, and two species of Equus. The application of a relatively recent fossil calibration point and a relatively ancient calibration point produced greatly different estimates of evolutionary rates and divergence times for both genes, even though a relative rates test did not find significant rate differences among taxa. A likelihood-ratio test, however, rejected a molecular clock for both genes. Neither calibration point produced estimates of divergence times consistent with paleontological evidence over a range of perissodactyl radiations. The combined analysis of both genes produces a well-resolved phylogeny with Perissodactyla that conforms to traditional views of interfamilial relationships and supports monophyly of neotropical tapirs. Combining the data sets increases support for most nodes but decreases the support for a neotropical tapir clade because the COII and 12S rRNA data sets are in conflict for tapir relationships. Received: 6 January 1999 / Accepted: 2 August 1999     

Integrating fossil preservation biases in the selection of calibrations for molecular divergence time estimation

The selection of fossil data to use as calibration age priors in molecular divergence time estimates inherently links neontological methods with paleontological theory. However, few neontological studies have taken into account the possibility of a taphonomic bias in the fossil record when developing approaches to fossil calibration selection. The Sppil-Rongis effect may bias the first appearance of a lineage toward the recent causing most objective calibration selection approaches to erroneously exclude appropriate calibrations or to incorporate multiple calibrations that are too young to accurately represent the divergence times of target lineages. Using turtles as a case study, we develop a Bayesian extension to the fossil selection approach developed by Marshall (2008. A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibrations points. Am. Nat. 171:726-742) that takes into account this taphonomic bias. Our method has the advantage of identifying calibrations that may bias age estimates to be too recent while incorporating uncertainty in phylogenetic parameter estimates such as tree topology and branch lengths. Additionally, this method is easily adapted to assess the consistency of potential calibrations to any one calibration in the candidate pool.     

Calibration of a molecular clock in tits (Paridae)--do nucleotide substitution rates of mitochondrial genes deviate from the 2% rule?

The ongoing debate on the reliability of avian molecular clocks is actually based on only a small number of calibrations carried out under different assumptions with respect to the choice and constraints of calibration points or to the use of substitution models. In this study, we provide substitution rate estimates for two mitochondrial genes, cytochrome b and the control region, and age estimates for lineage splits within four subgenera of tits (Paridae: Parus, Cyanistes, Poecile and Periparus). Overall sequence divergence between cytochrome b lineages covers a range of 0.4-1.8% per million years and is thus consistent with the frequently adopted approximation for a sequence divergence between avian lineages of 1.6-2% per my. Overall rate variation is high and encompasses the 2% value in a 95% CI for model corrected data. Mean rate estimates for cytochrome b range between 1.9 and 8.9 x 10(-3) substitutions per site per lineage. Local rates differ significantly between taxonomic levels with lowest estimates for haplotype lineages. At the population/subspecies level mean sequence divergence between lineages matches the 2% rule best for most cytochrome b datasets (1.5-1.9% per my) with maximum estimates for small isolated populations like those of the Canarian P. teneriffae complex (up to 3.9% per my). Overall rate estimates for the control region range at similar values like those for cytochrome b (2.7-8.8 x 10(-3), 0.5-1.8% per my), however, within some subgenera mean rates are higher than those for cytochrome b for uncorrected sequence data. The lowest rates for both genes were calculated for coal tits of subgenus Periparus (0.04-0.6% per my). Model-corrected sequence data tend to result in higher rate estimates than uncorrected data. Increase of the gamma shape parameter goes along with a significant decrease of rate and partly age estimates, too. Divergence times for earliest deep splits within tit subgenera Periparus and Parus were dated to the mid Miocene at 10-14my bp. Most recent splits between east and west Palearctic taxa of blue, willow and great tits were dated to the Pliocene/Pleistocene boundary with the earliest estimates based on model-corrected trees. Relaxation of the Messinian calibration point leads to more recent divergence times for North African coal and blue tit populations during the mid Pliocene. Despite a relatively broad age constraint for the split between Nearctic and Palearctic Poecile due to the Pliocene re-opening of the Bering Strait, the split between chickadees and willow tits is dated considerably earlier than in former studies to the upper bound of the age constraint at 7.4 my BP.     

Correlating Bayesian date estimates with climatic events and domestication using a bovine case study

The tribe Bovini contains a number of commercially and culturally important species, such as cattle. Understanding their evolutionary time scale is important for distinguishing between post-glacial and domestication-associated population expansions, but estimates of bovine divergence times have been hindered by a lack of reliable calibration points. We present a Bayesian phylogenetic analysis of 481 mitochondrial D-loop sequences, including 228 radiocarbon-dated ancient DNA sequences, using a multi-demographic coalescent model. By employing the radiocarbon dates as internal calibrations, we co-estimate the bovine phylogeny and divergence times in a relaxed-clock framework. The analysis yields evidence for significant population expansions in both taurine and zebu cattle, European aurochs and yak clades. The divergence age estimates support domestication-associated expansion times (less than 12 kyr) for the major haplogroups of cattle. We compare the molecular and palaeontological estimates for the Bison-Bos divergence.     

Assessing the sensitivity of divergence time estimates to locus sampling,calibration points,and model priors in a RAD-seq phylogeny of Carex section Schoenoxiphium

Restriction site-associated DNA sequencing (RAD-seq) and related methods have become relatively common approaches to resolve species-level phylogeny. It is not clear, however, whether RAD-seq data matrices are well suited to relaxed clock inference of divergence times, given the size of the matrices and the abundance of missing data. We investigated the sensitivity of Bayesian relaxed clock estimates of divergence times to alternative analytical decisions on an empirical RAD-seq phylogenetic matrix. We explored the relative contribution of secondary calibration strategies, amount of missing data, and the data partition analyzed to overall variance in divergence times inferred using BEAST MCMC analyses of Carex section Schoenoxiphium (Cyperaceae)—a recent radiation for which we have nearly complete species sampling of RAD-seq data. The crown node for Schoenoxiphium was estimated to be 15.22 (9.56–21.18) Ma using a single calibration point and low missing data, 11.93 (8.07–16.03) Ma using multiple calibration points and low missing data, and 8.34 (5.41–11.22) using multiple calibrations but high missing data. We found that using matrices with more than half of the individuals with missing data inferred younger mean ages for all nodes. Moreover, we have found that our molecular clock estimates are sensitive to the positions of the calibration(s) in our phylogenetic tree (using matrices with low missing data), especially when only a single calibration was applied to estimate divergence times. These results argue for sensitivity analyses and caution in interpreting divergence time estimates from RAD-seq data.     

Integration of Bayesian molecular clock methods and fossil-based soft bounds reveals early Cenozoic origin of African lacertid lizards

Background  

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.     

Birth-and-death evolution in primate MHC class I genes: divergence time estimates

The major histocompatibility complex (MHC) is a multigene family that mediates the host immune response by helping T lymphocytes to recognize and respond to foreign antigens. The high degree of polymorphism and a quick turnover of the genetic loci make the evolution of MHC genes an intriguing subject of study. To understand the evolutionary pattern of this multigene family, we studied the phylogeny and divergence times of six functional MHC class I loci from primate species. On the phylogenetic trees, locus F occupies the most basal position among these loci. Our results suggest that the F locus diverged from the other MHC class I loci about 46-66 MYA. The major diversification of the other class I loci was estimated to have occurred at about 35-49 MYA, which is before the time of separation of Old World-New World monkeys. The gene duplication leading to the classical C locus in great apes appears to have occurred about 21-28 MYA. At approximately the same time the duplication of the B locus occurred in macaques. The oldest allelic lineages of A, B, and C loci in humans seem to have appeared at least 14-19, 10-15, and 13-17 MYA, respectively. Our phylogenetic analysis supports the hypothesis that the nonclassical locus F has diverged from the rest of class I loci very early in primate evolution. The overall phylogenetic pattern observed among class I genes is consistent with the model of birth-and-death evolution.     

X-chromosomal window into the evolutionary history of the guenons (Primates: Cercopithecini)

Molecular studies of the guenons suggest that the arboreal Cercopithecus species form a monophyletic group within the tribe Cercopithecini. However, the evolutionary relationships among these arboreal congeners remain poorly resolved. The present work marks the first attempt to reconstruct the history of this group through the phylogenetic analysis of long nuclear sequences. We surveyed 19 guenons and seven outgroup taxa for a approximately 9.3kb fragment of X-chromosomal DNA homologous to a portion of human Xq13.3. Parsimony and maximum likelihood analyses of these sequences consistently recover two strongly-supported patterns within the arboreal Cercopithecus clade: (1) a clustering of members of the cephus and mitis species groups, and (2) a monophyletic aggregate including the mona, neglectus, and diana species groups. Although guenons occasionally hybridize in the wild, interbreeding forms of different species groups do not cluster together as sister-taxa in the X-chromosomal tree, suggesting that the two clades inferred here are not reticulate patterns due to recent gene flow. These clades are most likely the result of either ancestral hybridization or true phylogenetic history. We advocate the latter explanation because the same two aggregates (cephus/mitis and mona/neglectus/diana) are recovered, albeit with weak support, by a number of earlier analyses. Finally, X-chromosomal divergence dates are estimated for a number of nodes in the guenon radiation. The divergence of guenon and papionin lineages at 11.5 (+/-1.3) MYA appears to be a particularly robust estimate since it is inferred from both mitochondrial and X-chromosomal studies, each using different fossil calibration points.