Congruence test of molecular clock calibration hypotheses based on Bayes factor comparisons

HyRAD-X Exome Capture Museomics Unravels Giant Ground Beetle Evolution

Advances in phylogenomics contribute toward resolving long-standing evolutionary questions. Notwithstanding, genetic diversity contained within more than a billion biological specimens deposited in natural history museums remains recalcitrant to analysis owing to challenges posed by its intrinsically degraded nature. Yet that tantalizing resource could be critical in overcoming taxon sampling constraints hindering our ability to address major evolutionary questions. We addressed this impediment by developing phyloHyRAD, a new bioinformatic pipeline enabling locus recovery at a broad evolutionary scale from HyRAD-X exome capture of museum specimens of low DNA integrity using a benchtop RAD-derived exome-complexity-reduction probe set developed from high DNA integrity specimens. Our new pipeline can also successfully align raw RNAseq transcriptomic and ultraconserved element reads with the RAD-derived probe catalog. Using this method, we generated a robust timetree for Carabinae beetles, the lack of which had precluded study of macroevolutionary trends pertaining to their biogeography and wing-morphology evolution. We successfully recovered up to 2,945 loci with a mean of 1,788 loci across the exome of specimens of varying age. Coverage was not significantly linked to specimen age, demonstrating the wide exploitability of museum specimens. We also recovered fragmentary mitogenomes compatible with Sanger-sequenced mtDNA. Our phylogenomic timetree revealed a Lower Cretaceous origin for crown group Carabinae, with the extinct Aplothorax Waterhouse, 1841 nested within the genus Calosoma Weber, 1801 demonstrating the junior synonymy of Aplothorax syn. nov., resulting in the new combination Calosoma  burchellii (Waterhouse, 1841) comb. nov. This study compellingly illustrates that HyRAD-X and phyloHyRAD efficiently provide genomic-level data sets informative at deep evolutionary scales.

Dispersal versus vicariance in the Aegean: combining molecular and morphological phylogenies of eastern Mediterranean Dendarus (Coleoptera: Tenebrionidae) sheds new light on the phylogeography of the Aegean area

The Aegean archipelago, as an ‘evolutionary laboratory of nature’, is an ideal model for research in phylogeography. In this area, the darkling beetles of the genus Dendarus (distributed from Morocco to the Caucasus) exhibit a high level of diversity with 36 species, 27 of which are island endemics. However, their taxonomy is complex and unstable, having undergone continuous revision to address extensive morphological and ecological plasticity. Here, we examine the phylogenetic relationships of 23 species from Greece and Turkey, using mitochondrial and nuclear DNA sequences and 61 morphological characters, to unveil their phylogeny in the Aegean. This helps to clarify phylogeographic scenarios and historical processes that shaped the observed patterns. The analyses reveal 13 distinct lineages with several para- and polyphyletic cases that correspond to three major phylogroups [south/south-east Aegean (D. foraminosus complex, D. rhodius, D. sporadicus, D. wettsteini); central to north Aegean, Turkey and mainland Greece (D. crenulatus, D. moesiacus group, D. sinuatus complex, D. stygius) and mainland Greece (D. messenius, D. paganettii)], indicating the need for further taxonomic re-evaluation. Lineage topology and phylogeography suggest a spatial and temporal sequence of geographic isolation, following either a vicariant or a dispersal model coincident with major palaeogeographic separations in the Aegean.

Reliable Confidence Intervals for RelTime Estimates of Evolutionary Divergence Times

Confidence intervals (CIs) depict the statistical uncertainty surrounding evolutionary divergence time estimates. They capture variance contributed by the finite number of sequences and sites used in the alignment, deviations of evolutionary rates from a strict molecular clock in a phylogeny, and uncertainty associated with clock calibrations. Reliable tests of biological hypotheses demand reliable CIs. However, current non-Bayesian methods may produce unreliable CIs because they do not incorporate rate variation among lineages and interactions among clock calibrations properly. Here, we present a new analytical method to calculate CIs of divergence times estimated using the RelTime method, along with an approach to utilize multiple calibration uncertainty densities in dating analyses. Empirical data analyses showed that the new methods produce CIs that overlap with Bayesian highest posterior density intervals. In the analysis of computer-simulated data, we found that RelTime CIs show excellent average coverage probabilities, that is, the actual time is contained within the CIs with a 94% probability. These developments will encourage broader use of computationally efficient RelTime approaches in molecular dating analyses and biological hypothesis testing.

Dates and rates in grape's plastomes: evolution in slow motion

The family Vitaceae includes the domesticated grapevine (Vitis vinifera), one of the most economically important crops in the world. Despite the importance of Vitaceae, there is still considerable controversy surrounding their phylogenetic relationships and evolutionary timescales. Moreover, variation in rates of molecular evolution among Vitaceae remains mostly unexplored. The present research aims to fill these knowledge gaps through the analysis of plastome sequences. Thirteen newly sequenced grape plastomes are presented and their phylogenetic relationships examined. Divergence times and absolute substitution rates are inferred under different molecular clocks by the analysis of 95 non-coding plastid regions and 43 representative accessions of the major lineages of Vitaceae. Furthermore, the phylogenetic informativeness of non-coding plastid regions is investigated. We find strong evidence in favor of the random local clock model and rate heterogeneity within Vitaceae. Substitution rates decelerate in Ampelocissus, Ampelopsis, Nekemias, Parthenocissus, Rhoicissus, and Vitis, with genus Vitis showing the lowest values up to a minimum of ~ 4.65 × 10^-11 s/s/y. We suggest that liana-like species of Vitaceae evolve slower than erect growth habit plants and we invoke the "rate of mitosis hypothesis" to explain the observed pattern of the substitution rates. We identify a reduced set of 20 non-coding regions able to accurately reconstruct the phylogeny of Vitaceae and we provide a detailed description of all 152 non-coding regions identified in the plastomes of subg. Vitis. These polymorphic regions will find their applications in phylogenetics, phylogeography, and population genetics as well in grapes identification through DNA barcoding techniques.

Consequences of Secondary Calibrations on Divergence Time Estimates

Secondary calibrations (calibrations based on the results of previous molecular dating studies) are commonly applied in divergence time analyses in groups that lack fossil data; however, the consequences of applying secondary calibrations in a relaxed-clock approach are not fully understood. I tested whether applying the posterior estimate from a primary study as a prior distribution in a secondary study results in consistent age and uncertainty estimates. I compared age estimates from simulations with 100 randomly replicated secondary trees. On average, the 95% credible intervals of node ages for secondary estimates were significantly younger and narrower than primary estimates. The primary and secondary age estimates were significantly different in 97% of the replicates after Bonferroni corrections. Greater error in magnitude was associated with deeper than shallower nodes, but the opposite was found when standardized by median node age, and a significant positive relationship was determined between the number of tips/age of secondary trees and the total amount of error. When two secondary calibrated nodes were analyzed, estimates remained significantly different, and although the minimum and median estimates were associated with less error, maximum age estimates and credible interval widths had greater error. The shape of the prior also influenced error, in which applying a normal, rather than uniform, prior distribution resulted in greater error. Secondary calibrations, in summary, lead to a false impression of precision and the distribution of age estimates shift away from those that would be inferred by the primary analysis. These results suggest that secondary calibrations should not be applied as the only source of calibration in divergence time analyses that test time-dependent hypotheses until the additional error associated with secondary calibrations is more properly modeled to take into account increased uncertainty in age estimates.

Constraining the Deep Origin of Parasitic Flatworms and Host-Interactions with Fossil Evidence

Novel fossil discoveries have contributed to our understanding of the evolutionary appearance of parasitism in flatworms. Furthermore, genetic analyses with greater coverage have shifted our views on the coevolution of parasitic flatworms and their hosts. The putative record of parasitic flatworms is consistent with extant host associations and so can be used to put constraints on the evolutionary origin of the parasites themselves. The future lies in new molecular clock analyses combined with additional discoveries of exceptionally preserved flatworms associated with hosts and coprolites. Besides direct evidence, the host fossil record and biogeography have the potential to constrain their evolutionary history, albeit with caution needed to avoid circularity, and a need for calibrations to be implemented in the most conservative way. This might result in imprecise, but accurate divergence estimates for the evolution of parasitic flatworms.

Calibration uncertainty in molecular dating analyses: there is no substitute for the prior evaluation of time priors

Calibration is the rate-determining step in every molecular clock analysis and, hence, considerable effort has been expended in the development of approaches to distinguish good from bad calibrations. These can be categorized into a priori evaluation of the intrinsic fossil evidence, and a posteriori evaluation of congruence through cross-validation. We contrasted these competing approaches and explored the impact of different interpretations of the fossil evidence upon Bayesian divergence time estimation. The results demonstrate that a posteriori approaches can lead to the selection of erroneous calibrations. Bayesian posterior estimates are also shown to be extremely sensitive to the probabilistic interpretation of temporal constraints. Furthermore, the effective time priors implemented within an analysis differ for individual calibrations when employed alone and in differing combination with others. This compromises the implicit assumption of all calibration consistency methods, that the impact of an individual calibration is the same when used alone or in unison with others. Thus, the most effective means of establishing the quality of fossil-based calibrations is through a priori evaluation of the intrinsic palaeontological, stratigraphic, geochronological and phylogenetic data. However, effort expended in establishing calibrations will not be rewarded unless they are implemented faithfully in divergence time analyses.

The impact of calibration and clock-model choice on molecular estimates of divergence times

Phylogenetic estimates of evolutionary timescales can be obtained from nucleotide sequence data using the molecular clock. These estimates are important for our understanding of evolutionary processes across all taxonomic levels. The molecular clock needs to be calibrated with an independent source of information, such as fossil evidence, to allow absolute ages to be inferred. Calibration typically involves fixing or constraining the age of at least one node in the phylogeny, enabling the ages of the remaining nodes to be estimated. We conducted an extensive simulation study to investigate the effects of the position and number of calibrations on the resulting estimate of the timescale. Our analyses focused on Bayesian estimates obtained using relaxed molecular clocks. Our findings suggest that an effective strategy is to include multiple calibrations and to prefer those that are close to the root of the phylogeny. Under these conditions, we found that evolutionary timescales could be estimated accurately even when the relaxed-clock model was misspecified and when the sequence data were relatively uninformative. We tested these findings in a case study of simian foamy virus, where we found that shallow calibrations caused the overall timescale to be underestimated by up to three orders of magnitude. Finally, we provide some recommendations for improving the practice of molecular-clock calibration.