Studying the effect of environmental change on biotic evolution: past genetic contributions, current work and future directions

Phylo-Epigenetics in Phylogeny Analyses and Evolution

Long-standing, continuous blurring and controversies in the field of phylogenetic interspecies relations, associated with insufficient explanations for dynamics and variability of speeds of evolution in mammals, hint at a crucial missing link. It has been suggested that transgenerational epigenetic inheritance and the concealed mechanisms behind play a distinct role in mammalian evolution. Here, a comprehensive sequence alignment approach in hominid species, i.e., Homo sapiens, Homo neanderthalensis, Denisovan human, Pan troglodytes, Pan paniscus, Gorilla gorilla, and Pongo pygmaeus, comprising conserved CpG islands of housekeeping genes, uncover evidence for a distinct variability of CpG dinucleotides. Applying solely these evolutionary consistent and inconsistent CpG sites in a classic phylogenetic analysis, calibrated by the divergence time point of the common chimpanzee (P. troglodytes) and the bonobo or pygmy chimpanzee (P. paniscus), a "phylo-epigenetic" tree has been generated, which precisely recapitulates branch points and branch lengths, i.e., divergence events and relations, as they have been broadly suggested in the current literature, based on comprehensive molecular phylogenomics and fossil records of many decades. It is suggested here that CpG dinucleotide changes at CpG islands are of superior importance for evolutionary developments. These changes are successfully inherited through the germ line, determining emerging methylation profiles, and they are a central component of transgenerational epigenetic inheritance. It is hidden in the DNA, what will happen on it later.

Potential for bias and low precision in molecular divergence time estimation of the Canopy of Life: an example from aquatic bird families

Uncertainty in divergence time estimation is frequently studied from many angles but rarely from the perspective of phylogenetic node age. If appropriate molecular models and fossil priors are used, a multi-locus, partitioned analysis is expected to equally minimize error in accuracy and precision across all nodes of a given phylogeny. In contrast, if available models fail to completely account for rate heterogeneity, substitution saturation and incompleteness of the fossil record, uncertainty in divergence time estimation may increase with node age. While many studies have stressed this concern with regard to deep nodes in the Tree of Life, the inference that molecular divergence time estimation of shallow nodes is less sensitive to erroneous model choice has not been tested explicitly in a Bayesian framework. Because of available divergence time estimation methods that permit fossil priors across any phylogenetic node and the present increase in efficient, cheap collection of species-level genomic data, insight is needed into the performance of divergence time estimation of shallow (<10 MY) nodes. Here, we performed multiple sensitivity analyses in a multi-locus data set of aquatic birds with six fossil constraints. Comparison across divergence time analyses that varied taxon and locus sampling, number and position of fossil constraint and shape of prior distribution showed various insights. Deviation from node ages obtained from a reference analysis was generally highest for the shallowest nodes but determined more by temporal placement than number of fossil constraints. Calibration with only the shallowest nodes significantly underestimated the aquatic bird fossil record, indicating the presence of saturation. Although joint calibration with all six priors yielded ages most consistent with the fossil record, ages of shallow nodes were overestimated. This bias was found in both mtDNA and nDNA regions. Thus, divergence time estimation of shallow nodes may suffer from bias and low precision, even when appropriate fossil priors and best available substitution models are chosen. Much care must be taken to address the possible ramifications of substitution saturation across the entire Tree of Life.

Back to BaySICS: a user-friendly program for Bayesian Statistical Inference from Coalescent Simulations

Inference of population demographic history has vastly improved in recent years due to a number of technological and theoretical advances including the use of ancient DNA. Approximate Bayesian computation (ABC) stands among the most promising methods due to its simple theoretical fundament and exceptional flexibility. However, limited availability of user-friendly programs that perform ABC analysis renders it difficult to implement, and hence programming skills are frequently required. In addition, there is limited availability of programs able to deal with heterochronous data. Here we present the software BaySICS: Bayesian Statistical Inference of Coalescent Simulations. BaySICS provides an integrated and user-friendly platform that performs ABC analyses by means of coalescent simulations from DNA sequence data. It estimates historical demographic population parameters and performs hypothesis testing by means of Bayes factors obtained from model comparisons. Although providing specific features that improve inference from datasets with heterochronous data, BaySICS also has several capabilities making it a suitable tool for analysing contemporary genetic datasets. Those capabilities include joint analysis of independent tables, a graphical interface and the implementation of Markov-chain Monte Carlo without likelihoods.

Using classical population genetics tools with heterochroneous data: time matters!

BACKGROUND

New polymorphism datasets from heterochroneous data have arisen thanks to recent advances in experimental and microbial molecular evolution, and the sequencing of ancient DNA (aDNA). However, classical tools for population genetics analyses do not take into account heterochrony between subsets, despite potential bias on neutrality and population structure tests. Here, we characterize the extent of such possible biases using serial coalescent simulations.

METHODOLOGY/PRINCIPAL FINDINGS

We first use a coalescent framework to generate datasets assuming no or different levels of heterochrony and contrast most classical population genetic statistics. We show that even weak levels of heterochrony ( approximately 10% of the average depth of a standard population tree) affect the distribution of polymorphism substantially, leading to overestimate the level of polymorphism theta, to star like trees, with an excess of rare mutations and a deficit of linkage disequilibrium, which are the hallmark of e.g. population expansion (possibly after a drastic bottleneck). Substantial departures of the tests are detected in the opposite direction for more heterochroneous and equilibrated datasets, with balanced trees mimicking in particular population contraction, balancing selection, and population differentiation. We therefore introduce simple corrections to classical estimators of polymorphism and of the genetic distance between populations, in order to remove heterochrony-driven bias. Finally, we show that these effects do occur on real aDNA datasets, taking advantage of the currently available sequence data for Cave Bears (Ursus spelaeus), for which large mtDNA haplotypes have been reported over a substantial time period (22-130 thousand years ago (KYA)).

CONCLUSIONS/SIGNIFICANCE

Considering serial sampling changed the conclusion of several tests, indicating that neglecting heterochrony could provide significant support for false past history of populations and inappropriate conservation decisions. We therefore argue for systematically considering heterochroneous models when analyzing heterochroneous samples covering a large time scale.

Fire and ice: genetic structure of the Uinta ground squirrel (Spermophilus armatus) across the Yellowstone hotspot

The range of the Uinta ground squirrel, Spermophilus armatus, is centred over one of the most tectonically active regions today, the Yellowstone hotspot. We document the role of Quaternary tectonic and climatic history on the genetic structure of this species by screening museum and extant individuals throughout its range. Phylogeographic, divergence time, and demographic analyses of partial mitochondrial cytochrome b and control region DNA sequences yield insight into the cadence of evolution across three spatiotemporal scales: (i) a relatively deep intraspecific divergence of S. armatus into three lineages coincident with the last major volcanic eruption in the region and maintained by the Snake River Plain; (ii) demographic expansion in two lineages corresponding to the time of last deglaciation of the region; and (iii) a recent (< 50 years) local extinction of the third lineage coincident with climatic change and conversion of habitat for agricultural purposes in eastern Idaho. Beyond these inferences, our study highlights the unique value of museum material to phylogeography, and shows that small mammal recolonization of previously glaciated montane 'islands' differs from northward postglacial expansion observed in areas previously covered by continental ice sheets. Montane 'islands' may harbour high genetic diversity because of admixture and recurrent expansion/extinction.