Prolonged decay of molecular rate estimates for metazoan mitochondrial DNA

Diversity of Land Snail Tribe Helicini (Gastropoda: Stylommatophora: Helicidae): Where Do We Stand after 20 Years of Sequencing Mitochondrial Markers?

Sequences of mitochondrial genes revolutionized the understanding of animal diversity and continue to be an important tool in biodiversity research. In the tribe Helicini, a prominent group of the western Palaearctic land snail fauna, mitochondrial data accumulating since the 2000s helped to newly delimit genera, inform species-level taxonomy and reconstruct past range dynamics. We combined the published data with own unpublished sequences and provide a detailed overview of what they revealed about the diversity of the group. The delimitation of Helix is revised by placing Helix godetiana back in the genus and new synonymies are suggested within the genera Codringtonia and Helix. The spatial distribution of intraspecific mitochondrial lineages of several species is shown for the first time. Comparisons between species reveal considerable variation in distribution patterns of intraspecific lineages, from broad postglacial distributions to regions with a fine-scale pattern of allopatric lineage replacement. To provide a baseline for further research and information for anyone re-using the data, we thoroughly discuss the gaps in the current dataset, focusing on both taxonomic and geographic coverage. Thanks to the wealth of data already amassed and the relative ease with which they can be obtained, mitochondrial sequences remain an important source of information on intraspecific diversity over large areas and taxa.

The Tien Shan vole (Microtus ilaeus; Rodentia: Cricetidae) as a new species in the Late Pleistocene of Europe

Grey voles (subgenus Microtus) represent a complex of at least seven closely related and partly cryptic species. The range of these species extends from the Atlantic to the Altai Mountains, but most of them occur east of the Black Sea. Using ancient DNA analyses of the Late Pleistocene specimens, we identified a new mtDNA lineage of grey voles in Europe. Phylogenetic analysis of mitochondrial DNA cytochrome b sequences from 23 voles from three caves, namely, Emine-Bair-Khosar (Crimea, Ukraine), Cave 16 (Bulgaria), and Bacho Kiro (Bulgaria), showed that 14 specimens form a previously unrecognized lineage, sister to the Tien Shan vole. The average sequence divergence of this lineage and the extant Tien Shan vole was 4.8%, which is similar to the divergence of grey vole forms, which are considered distinct species or being on the verge of speciation; M. arvalis and M. obscurus or M. mystacinus and M. rossiaemeridionalis. We estimated the time to the most recent common ancestor of the grey voles to be 0.66 Ma, which is over twice the recent estimates, while the divergence of the extant Tien Shan vole and the new lineage to be 0.29 Ma. Our discovery suggests that grey voles may have been more diversified in the past and that their ranges may have differed substantially from current ones. It also underlines the utility of ancient DNA to decipher the evolutionary history of voles.

Genetic Diversity of Ligidium Isopods in Hokkaido and Niigata, Northern Japan, Based on Mitochondrial DNA Analysis

The genetic diversity of the genus Ligidium in Hokkaido and Niigata, northern Japan, was investigated by analyzing the cytochrome c oxidase subunit 1 (CO1) region in the mitochondrial DNA (mtDNA). The genetic diversity in Hokkaido was much lower than that in Niigata. Nine different operational taxonomic units (OTUs) were identified. Only a single OTU, most likely Ligidium japonicum, was found in Hokkaido, whereas all nine OTUs were found in Niigata. Using the mtDNA evolutionary rate determined for the marine invertebrate Haptosquilla pulchella (Miers, 1880), population expansion for OTU1 in Hokkaido was estimated to have occurred at 12,600 years BP, suggesting that Ligidium underwent a bottleneck due to glacial cooling, and the population then expanded after postglacial warming. Assuming that the expansion of the OTU1 population occurred at 9600 years BP, when the sea surface temperature rose offshore of Tokachi in the Northwestern Pacific, the evolutionary rate (µ) of the mtDNA CO1 region in Ligidium is calculated as: 0.087 (95% confidence intervals: min: 0.042-max: 0.12) (substitutions/site/million years). The presence of a haplotype common to Hokkaido and Niigata implies that the haplotype migrated across the Tsugaru Strait. Considering that geological evidence indicates that the Tsugaru Strait was continuously present even during the last glacial maximum when the sea level was at its lowest, accidental transport by human beings or animals might have been critical to the migration of Ligidium.

When parasites persist: tapeworms survive host extinction and reveal waves of dispersal across Beringia

Investigations of intercontinental dispersal between Asia and North America reveal complex patterns of geographic expansion, retraction and isolation, yet historical reconstructions are largely limited by the depth of the record that is retained in patterns of extant diversity. Parasites offer a tool for recovering deep historical insights about the biosphere, improving the resolution of past community-level interactions. We explored biogeographic hypotheses regarding the history of dispersal across Beringia, the region intermittently linking Asia and North America, through large-scale multi-locus phylogenetic analyses of the genus Schizorchis, an assemblage of host-specific cestodes in pikas (Lagomorpha: Ochotonidae). Our genetic data support palaeontological evidence for two separate geographic expansions into North America by Ochotona in the late Tertiary, a history that genomic evidence from extant pikas does not record. Pikas descending from the first colonization of Miocene age persisted into the Pliocene, subsequently coming into contact with a second wave of Nearctic colonists from Eurasia before going extinct. Spatial and temporal overlap of historically independent pika populations provided a window for host colonization, allowing persistence of an early parasite lineage in the contemporary fauna following the extinction of its ancestral hosts. Empirical evidence for ancient 'ghost assemblages' of hosts and parasites demonstrates how complex mosaic faunas are assembled in the biosphere through episodes of faunal mixing encompassing parasite lineages across deep and shallow time.

Mitochondrial DNAs provide insight into trypanosome phylogeny and molecular evolution

Background

Trypanosomes are single-celled eukaryotic parasites characterised by the unique biology of their mitochondrial DNA. African livestock trypanosomes impose a major burden on agriculture across sub-Saharan Africa, but are poorly understood compared to those that cause sleeping sickness and Chagas disease in humans. Here we explore the potential of the maxicircle, a component of trypanosome mitochondrial DNA to study the evolutionary history of trypanosomes.

Results

We used long-read sequencing to completely assemble maxicircle mitochondrial DNA from four previously uncharacterized African trypanosomes, and leveraged these assemblies to scaffold and assemble a further 103 trypanosome maxicircle gene coding regions from published short-read data. While synteny was largely conserved, there were repeated, independent losses of Complex I genes. Comparison of pre-edited and non-edited genes revealed the impact of RNA editing on nucleotide composition, with non-edited genes approaching the limits of GC loss. African tsetse-transmitted trypanosomes showed high levels of RNA editing compared to other trypanosomes. The gene coding regions of maxicircle mitochondrial DNAs were used to construct time-resolved phylogenetic trees, revealing deep divergence events among isolates of the pathogens Trypanosoma brucei and T. congolense.

Conclusions

Our data represents a new resource for experimental and evolutionary analyses of trypanosome phylogeny, molecular evolution and function. Molecular clock analyses yielded a timescale for trypanosome evolution congruent with major biogeographical events in Africa and revealed the recent emergence of Trypanosoma brucei gambiense and T. equiperdum, major human and animal pathogens.

Counterbalancing the time-dependent effect on the human mitochondrial DNA molecular clock

BACKGROUND

The molecular clock is an important genetic tool for estimating evolutionary timescales. However, the detection of a time-dependent effect on substitution rate estimates complicates its application. It has been suggested that demographic processes could be the main cause of this confounding effect. In the present study, I propose a new algorithm for estimating the coalescent age of phylogenetically related sequences, taking into account the observed time-dependent effect on the molecular rate detected by others.

RESULTS

By applying this method to real human mitochondrial DNA trees with shallow and deep topologies, I obtained significantly older molecular ages for the main events of human evolution than were previously estimated. These ages are in close agreement with the most recent archaeological and paleontological records favoring the emergence of early anatomically modern humans in Africa 315 ± 34 thousand years ago (kya) and the presence of recent modern humans outside of Africa as early as 174 ± 48 thousand years ago. Furthermore, during the implementation process, I demonstrated that in a population with fluctuating sizes, the probability of fixation of a new neutral mutant depends on the effective population size, which is in better accordance with the fact that under the neutral theory of molecular evolution, the fate of a molecular mutation is mainly determined by random drift.

CONCLUSIONS

I suggest that the demographic history of populations has a more decisive effect than purifying selection and/or mutational saturation on the time-dependent effect observed for the substitution rate, and I propose a new method that corrects for this effect.

Bayesian Inference of Evolutionary Histories under Time-Dependent Substitution Rates

Many factors complicate the estimation of time scales for phylogenetic histories, requiring increasingly complex evolutionary models and inference procedures. The widespread application of molecular clock dating has led to the insight that evolutionary rate estimates may vary with the time frame of measurement. This is particularly well established for rapidly evolving viruses that can accumulate sequence divergence over years or even months. However, this rapid evolution stands at odds with a relatively high degree of conservation of viruses or endogenous virus elements over much longer time scales. Building on recent insights into time-dependent evolutionary rates, we develop a formal and flexible Bayesian statistical inference approach that accommodates rate variation through time. We evaluate the novel molecular clock model on a foamy virus cospeciation history and a lentivirus evolutionary history and compare the performance to other molecular clock models. For both virus examples, we estimate a similarly strong time-dependent effect that implies rates varying over four orders of magnitude. The application of an analogous codon substitution model does not implicate long-term purifying selection as the cause of this effect. However, selection does appear to affect divergence time estimates for the less deep evolutionary history of the Ebolavirus genus. Finally, we explore the application of our approach on woolly mammoth ancient DNA data, which shows a much weaker, but still important, time-dependent rate effect that has a noticeable impact on node age estimates. Future developments aimed at incorporating more complex evolutionary processes will further add to the broad applicability of our approach.

Unexpected low genetic variation in the South American hystricognath rodent Lagostomus maximus (Rodentia: Chinchillidae)

The South American plains vizcacha, Lagostomus maximus inhabits primarily the Pampean and adjoining Espinal, Monte and Chaquenean regions of Argentina. In order to study the population genetic structure of L. maximus, a fragment of 560 bp of the mitochondrial DNA hypervariable region 1from 90 individuals collected from the 3 subspecies and 8 groups along Argentina was amplified and analyzed. We found 9 haplotypes. The haplotype network did not show an apparent phylogeographical signal. Although low levels of genetic variation were found in all the subspecies and groups analyzed, a radiation of L. maximus would have occurred from the North and Center of the Pampean region toward the rest of its geographic range in Argentina. Low levels of genetic diversity, the existence of a single genetically distinct population in Argentina and changes of its effective size indicate that metapopulation processes and changes in human population dynamics during the late-Holocene were important factors shaping the population genetic structure of L. maximus in Argentina.

Loss and gain of sexual reproduction in the same stick insect

The outcome of competition between different reproductive strategies within a single species can be used to infer selective advantage of the winning strategy. Where multiple populations have independently lost or gained sexual reproduction it is possible to investigate whether the advantage is contingent on local conditions. In the New Zealand stick insect Clitarchus hookeri, three populations are distinguished by recent change in reproductive strategy and we determine their likely origins. One parthenogenetic population has established in the United Kingdom and we provide evidence that sexual reproduction has been lost in this population. We identify the sexual population from which the parthenogenetic population was derived, but show that the UK females have a post‐mating barrier to fertilisation. We also demonstrate that two sexual populations have recently arisen in New Zealand within the natural range of the mtDNA lineage that otherwise characterizes parthenogenesis in this species. We infer independent origins of males at these two locations using microsatellite genotypes. In one population, a mixture of local and nonlocal alleles suggested males were the result of invasion. Males in another population were most probably the result of loss of an X chromosome that produced a male phenotype in situ. Two successful switches in reproductive strategy suggest local competitive advantage for outcrossing over parthenogenetic reproduction. Clitarchus hookeri provides remarkable evidence of repeated and rapid changes in reproductive strategy, with competitive outcomes dependent on local conditions.

Estimation of Evolutionary Rates of Mitochondrial DNA in Two Japanese Wood Mouse Species Based on Calibrations with Quaternary Environmental Changes

Reliable estimates of evolutionary rates of mitochondrial DNA might allow us to build realistic evolutionary scenarios covering broad time scales based on phylogenetic inferences. In the present study, we sought to obtain estimates of evolutionary rates in murine rodents using calibrations against historical biogeographic events. We first assumed that land-bridge-like structures that appeared intermittently at glacial maxima with 100,000-year intervals shaped the divergence patterns of cytochrome b (Cytb) sequences (1140 bp) of the larger Japanese wood mouse Apodemus speciosus. The comparison of sequences from peripheral remote islands that are separated from one another by deep straits allowed us to estimate mitochondrial DNA evolutionary rates (substitutions/site/million years) to be 0.027 to 0.036, with presumed calibrations from 140,000, 250,000, 350,000, and 440,000 years ago. Second, we addressed rapid expansion events inferred from analyses of the Cytb sequences of the lesser Japanese wood mouse A. argenteus. We detected five expansion signals in the dataset and established three categories based on the expansion parameter tau values: 3.9, 5.6-5.7, and 7.8-8.1. Considering that the climate became warmer 15,000, 53,000, and 115,000 years ago after preceding periods of rapid cooling, we calculated evolutionary rates to be 0.114, 0.047, and 0.031, respectively. This preliminary concept of the evolutionary rates on a time scale from 15,000 to 440,000 years ago for the wood mouse should be refined and tested in other species of murine rodents, including mice and rats.

A comparison of methods for estimating substitution rates from ancient DNA sequence data

Background

Phylogenetic analysis of DNA from modern and ancient samples allows the reconstruction of important demographic and evolutionary processes. A critical component of these analyses is the estimation of evolutionary rates, which can be calibrated using information about the ages of the samples. However, the reliability of these rate estimates can be negatively affected by among-lineage rate variation and non-random sampling. Using a simulation study, we compared the performance of three phylogenetic methods for inferring evolutionary rates from time-structured data sets: regression of root-to-tip distances, least-squares dating, and Bayesian inference. We also applied these three methods to time-structured mitogenomic data sets from six vertebrate species.

Results

Our results from 12 simulation scenarios show that the three methods produce reliable estimates when the substitution rate is high, rate variation is low, and samples of similar ages are not all grouped together in the tree (i.e., low phylo-temporal clustering). The interaction of these factors is particularly important for least-squares dating and Bayesian estimation of evolutionary rates. The three estimation methods produced consistent estimates of rates across most of the six mitogenomic data sets, with sequence data from horses being an exception.

Conclusions

We recommend that phylogenetic studies of ancient DNA sequences should use multiple methods of inference and test for the presence of temporal signal, among-lineage rate variation, and phylo-temporal clustering in the data.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1192-3) contains supplementary material, which is available to authorized users.

Ice age unfrozen: severe effect of the last interglacial, not glacial, climate change on East Asian avifauna

BACKGROUND

The glacial-interglacial cycles in the Pleistocene caused repeated range expansion and contraction of species in several regions in the world. However, it remains uncertain whether such climate oscillations had similar impact on East Asian biota, despite its widely recognized importance in global biodiversity. Here we use both molecular and ecological niche profiles on 11 East Asian avian species with various elevational ranges to reveal their response to the late Pleistocene climate changes.

RESULTS

The ecological niche models (ENM) consistently showed that these avian species might substantially contract their ranges to the south during the Last Interglacial period (LIG) and expanded their northern range margins through the Last Glacial Maximum (LGM), leading to the LGM ranges observed for all 11 species. Consistently, coalescent simulations based on 25-30 nuclear genes retrieved signatures of significant population growth through the last glacial period across all species studied. Climate statistics suggested that high climatic variability during the LIG and a relatively mild climate at the LGM potentially explained the historical population dynamics of these birds.

CONCLUSIONS

This is the first study based on multiple species and both lines of ecological niche profiles and genetic data to characterize the unique response of East Asian biota to late Pleistocene climate. The present study highlights regional differences in the evolutionary consequence of climate change during the last glacial cycle and implies that global warming might pose a great risk to species in this region given potentially higher climatic variation in the future analogous to that during the LIG.

Phylogeography of the Chydorus sphaericus Group (Cladocera: Chydoridae) in the Northern Palearctic

The biodiversity and the biogeography are still poorly understood for freshwater invertebrates. The crustacean Chydorus sphaericus-brevilabris complex (Cladocera: Chydoridae) is composed of species that are important components of Holarctic freshwater food webs. Recent morphological and genetic study of the complex has indicated a substantial species diversity in the northern hemisphere. However, we know little of the geographic boundaries of these novel lineages. Moreover, a large section of the Palearctic remains unexamined at the genetic level. Here we attempt to address the biodiversity knowledge gap for the Chydorus sphaericus group in the central Palearctic and assess its diversity and biogeographic boundaries. We sequenced nuclear (ITS-2) and mitochondrial (COI) gene regions of Chydorus specimens across the Palearctic and compared them with already available Holarctic sequences. We detected six main clades in the C. sphaericus group in the Palearctic, of which two of the groups are novel. Three of the more divergent clades are geographically widespread. The central portion of Eurasia (the Yenisey River basin) appears to be a narrow zone of secondary contact for phylogroups that expanded from European and Beringian refugia. As such, the previously unsampled central Palearctic represents an important region for understanding the evolutionary consequences of Pleistocene climatic oscillations on the Chydorus sphaericus group.

Genome-scale rates of evolutionary change in bacteria

Estimating the rates at which bacterial genomes evolve is critical to understanding major evolutionary and ecological processes such as disease emergence, long-term host–pathogen associations and short-term transmission patterns. The surge in bacterial genomic data sets provides a new opportunity to estimate these rates and reveal the factors that shape bacterial evolutionary dynamics. For many organisms estimates of evolutionary rate display an inverse association with the time-scale over which the data are sampled. However, this relationship remains unexplored in bacteria due to the difficulty in estimating genome-wide evolutionary rates, which are impacted by the extent of temporal structure in the data and the prevalence of recombination. We collected 36 whole genome sequence data sets from 16 species of bacterial pathogens to systematically estimate and compare their evolutionary rates and assess the extent of temporal structure in the absence of recombination. The majority (28/36) of data sets possessed sufficient clock-like structure to robustly estimate evolutionary rates. However, in some species reliable estimates were not possible even with ‘ancient DNA’ data sampled over many centuries, suggesting that they evolve very slowly or that they display extensive rate variation among lineages. The robustly estimated evolutionary rates spanned several orders of magnitude, from approximately 10⁻⁵ to 10⁻⁸ nucleotide substitutions per site year⁻¹. This variation was negatively associated with sampling time, with this relationship best described by an exponential decay curve. To avoid potential estimation biases, such time-dependency should be considered when inferring evolutionary time-scales in bacteria.

Mitogenomic analysis of a 50-generation chicken pedigree reveals a rapid rate of mitochondrial evolution and evidence for paternal mtDNA inheritance

Mitochondrial genomes represent a valuable source of data for evolutionary research, but studies of their short-term evolution have typically been limited to invertebrates, humans and laboratory organisms. Here we present a detailed study of 12 mitochondrial genomes that span a total of 385 transmissions in a well-documented 50-generation pedigree in which two lineages of chickens were selected for low and high juvenile body weight. These data allowed us to test the hypothesis of time-dependent evolutionary rates and the assumption of strict maternal mitochondrial transmission, and to investigate the role of mitochondrial mutations in determining phenotype. The identification of a non-synonymous mutation in ND4L and a synonymous mutation in CYTB, both novel mutations in Gallus, allowed us to estimate a molecular rate of 3.13 × 10(-7) mutations/site/year (95% confidence interval 3.75 × 10(-8)-1.12 × 10(-6)). This is substantially higher than avian rate estimates based upon fossil calibrations. Ascertaining which of the two novel mutations was present in an additional 49 individuals also revealed an instance of paternal inheritance of mtDNA. Lastly, an association analysis demonstrated that neither of the point mutations was strongly associated with the phenotypic differences between the two selection lines. Together, these observations reveal the highly dynamic nature of mitochondrial evolution over short time periods.

Time-dependent estimates of molecular evolutionary rates: evidence and causes

We are writing in response to a recent critique by Emerson & Hickerson (2015), who challenge the evidence of a time-dependent bias in molecular rate estimates. This bias takes the form of a negative relationship between inferred evolutionary rates and the ages of the calibrations on which these estimates are based. Here, we present a summary of the evidence obtained from a broad range of taxa that supports a time-dependent bias in rate estimates, with a consideration of the potential causes of these observed trends. We also describe recent progress in improving the reliability of evolutionary rate estimation and respond to the concerns raised by Emerson & Hickerson (2015) about the validity of rates estimated from time-structured sequence data. In doing so, we hope to dispel some misconceptions and to highlight several research directions that will improve our understanding of time-dependent biases in rate estimates.

Large-scale mitogenomics enables insights into Schizophora (Diptera) radiation and population diversity

True flies are insects of the order Diptera and encompass one of the most diverse groups of animals on Earth. Within dipterans, Schizophora represents a recent radiation of insects that was used as a model to develop a pipeline for generating complete mitogenomes using various sequencing platforms and strategies. 91 mitogenomes from 32 different species were sequenced and assembled with high fidelity, using amplicon, whole genome shotgun or single molecule sequencing approaches. Based on the novel mitogenomes, we estimate the origin of Schizophora within the Cretaceous-Paleogene (K-Pg) boundary, about 68.3 Ma. Detailed analyses of the blowfly family (Calliphoridae) place its origin at 22 Ma, concomitant with the radiation of grazing mammals. The emergence of ectoparasitism within calliphorids was dated 6.95 Ma for the screwworm fly and 2.3 Ma for the Australian sheep blowfly. Varying population histories were observed for the blowfly Chrysomya megacephala and the housefly Musca domestica samples in our dataset. Whereas blowflies (n = 50) appear to have undergone selective sweeps and/or severe bottlenecks in the New World, houseflies (n = 14) display variation among populations from different zoogeographical zones and low levels of gene flow. The reported high-throughput mitogenomics approach for insects enables new insights into schizophoran diversity and population history of flies.