Mitochondrial DNA as a tool for reconstructing past life-history traits in mammals

Revised phylogeny of mouflon based on expanded sampling of mitogenomes

Mouflons are flagship species of the Mediterranean islands where they persist. Once thought to be the remnants of a European wild sheep population, archaeology suggests they were introduced by humans to the islands of Cyprus in the Early Neolithic (~10,000 years ago) and later to Corsica and Sardinia. Their status as truly wild animals remains a subject of debate. To investigate the phylogenetic relationship between these island populations and other domestic and wild sheep from the Mediterranean region, we sequenced 50 mitogenomes of mouflons from Sardinia and Corsica, and modern and ancient Sardinian domestic sheep. A total of 68 additional publicly available mitogenomes were included in the comparative analysis and used to reconstruct the phylogeny of sheep and its closest wild relative, the mouflon (Ovis gmelini). Our study analyzed the evolutionary relationships within the C-E-X and haplogroup B clusters, showing that: a) Cyprus mouflons are more related to Anatolian and Iranian mouflons belonging to the wild haplogroup X, which seems to be basal to the domestic C and E haplogroups; b) Corsican and Sardinian mouflon arise from basal lineages associated with the early European expansion of domestic sheep. These results highlight the phylogenetic distinctiveness of the mouflon populations from the Mediterranean islands, suggesting a revision of their systematic classification and an update of the nomenclature for Sardinian and Corsican mouflons from the current status of subspecies of domestic sheep (Ovis aries musimon) to subspecies of their wild relatives (Ovis gmelini musimon) which would facilitate conservation efforts.

When paleontology meets genomics: complete mitochondrial genomes of two saber-toothed cats' species (Felidae: Machairodontinae)

Within the Machairodontinae subfamily, commonly referred to as saber-toothed cats, it is worth noting that only two species, namely Homotherium latidens, recognized as the scimitar-toothed cat, and Smilodon populator, renowned as the saber-toothed tiger, possess partial mitochondrial genomes accessible in the NCBI database. These sequences, however, do not include the mitogenome control region (mtDNA control region) and have several gaps in their genes, including protein-coding genes (PCGs) that are widely used in phylogenetic analysis. In this study, we aimed to obtain a complete assembly of the mitogenomes of these two species from next-generation sequencing data available at NCBI's SRA. The de novo assemblies showed complete mitogenomes with 17,323bp (H. latidens) and 16,769 bp (S. populator), both with 13 PCGs, 22tRNAs, two rRNAs and the mtDNA control region, with all genes following the standard order and position of most vertebrate mitogenomes. Despite being generally very similar to previous studies, our phylogeny and molecular dating reveals an earliest divergence between North American and North Sea H. latidens specimens which may be related to an Early Pleistocene migration across Beringia.

Reconstructing the History of Variation in Effective Population Size along Phylogenies

The nearly neutral theory predicts specific relations between effective population size (Ne) and patterns of divergence and polymorphism, which depend on the shape of the distribution of fitness effects (DFE) of new mutations. However, testing these relations is not straightforward, owing to the difficulty in estimating Ne. Here, we introduce an integrative framework allowing for an explicit reconstruction of the phylogenetic history of Ne, thus leading to a quantitative test of the nearly neutral theory and an estimation of the allometric scaling of the ratios of nonsynonymous over synonymous polymorphism (πN/πS) and divergence (dN/dS) with respect to Ne. As an illustration, we applied our method to primates, for which the nearly neutral predictions were mostly verified. Under a purely nearly neutral model with a constant DFE across species, we find that the variation in πN/πS and dN/dS as a function of Ne is too large to be compatible with current estimates of the DFE based on site frequency spectra. The reconstructed history of Ne shows a 10-fold variation across primates. The mutation rate per generation u, also reconstructed over the tree by the method, varies over a 3-fold range and is negatively correlated with Ne. As a result of these opposing trends for Ne and u, variation in πS is intermediate, primarily driven by Ne but substantially influenced by u. Altogether, our integrative framework provides a quantitative assessment of the role of Ne and u in modulating patterns of genetic variation, while giving a synthetic picture of their history over the clade.

Long-Lived Species of Bivalves Exhibit Low MT-DNA Substitution Rates

Bivalves represent valuable taxonomic group for aging studies given their wide variation in longevity (from 1–2 to >500 years). It is well known that aging is associated to the maintenance of Reactive Oxygen Species homeostasis and that mitochondria phenotype and genotype dysfunctions accumulation is a hallmark of these processes. Previous studies have shown that mitochondrial DNA mutation rates are linked to lifespan in vertebrate species, but no study has explored this in invertebrates. To this end, we performed a Bayesian Phylogenetic Covariance model of evolution analysis using 12 mitochondrial protein-coding genes of 76 bivalve species. Three life history traits (maximum longevity, generation time and mean temperature tolerance) were tested against 1) synonymous substitution rates (dS), 2) conservative amino acid replacement rates (Kc) and 3) ratios of radical over conservative amino acid replacement rates (Kr/Kc). Our results confirm the already known correlation between longevity and generation time and show, for the first time in an invertebrate class, a significant negative correlation between dS and longevity. This correlation was not as strong when generation time and mean temperature tolerance variations were also considered in our model (marginal correlation), suggesting a confounding effect of these traits on the relationship between longevity and mtDNA substitution rate. By confirming the negative correlation between dS and longevity previously documented in birds and mammals, our results provide support for a general pattern in substitution rates.

Population demographic history and adaptability of the vulnerable Lolokou Sucker Frog

Amphibians are experiencing worldwide declines due to increasing anthropogenetic disturbances. However, the genetic variability and hence adaptability are still unknown for most frogs. We integrated the mitochondrial (ND2 gene), nuclear (TYR gene) and major histocompatibility complex (MHC) loci, to clarify the demographic patterns and immune-gene diversity of the Lolokou Sucker Frog (Amolops loloensis). Demographic analysis of the ND2 and TYR genes suggested that the Lolokou Sucker Frog experienced a population expansion within the last 10,000 years. High MHC diversity was detected, which has likely resulted from positive selection, indicating the current diversity bodes well for the species' adaptive potential to pathogenic challenges. These findings broaden our knowledge on the population history and evolution adaptation of the reclusive torrent frog, and conservation implications are provided.

Mitochondrial genomics reveals the evolutionary history of the porpoises (Phocoenidae) across the speciation continuum

Historical variation in food resources is expected to be a major driver of cetacean evolution, especially for the smallest species like porpoises. Despite major conservation issues among porpoise species (e.g., vaquita and finless), their evolutionary history remains understudied. Here, we reconstructed their evolutionary history across the speciation continuum. Phylogenetic analyses of 63 mitochondrial genomes suggest that porpoises radiated during the deep environmental changes of the Pliocene. However, all intra-specific subdivisions were shaped during the Quaternary glaciations. We observed analogous evolutionary patterns in both hemispheres associated with convergent evolution to coastal versus oceanic environments. This suggests that similar mechanisms are driving species diversification in northern (harbor and Dall's) and southern species (spectacled and Burmeister's). In contrast to previous studies, spectacled and Burmeister's porpoises shared a more recent common ancestor than with the vaquita that diverged from southern species during the Pliocene. The low genetic diversity observed in the vaquita carried signatures of a very low population size since the last 5,000 years. Cryptic lineages within Dall's, spectacled and Pacific harbor porpoises suggest a richer evolutionary history than previously suspected. These results provide a new perspective on the mechanisms driving diversification in porpoises and an evolutionary framework for their conservation.

A mixed relaxed clock model

Over recent years, several alternative relaxed clock models have been proposed in the context of Bayesian dating. These models fall in two distinct categories: uncorrelated and autocorrelated across branches. The choice between these two classes of relaxed clocks is still an open question. More fundamentally, the true process of rate variation may have both long-term trends and short-term fluctuations, suggesting that more sophisticated clock models unfolding over multiple time scales should ultimately be developed. Here, a mixed relaxed clock model is introduced, which can be mechanistically interpreted as a rate variation process undergoing short-term fluctuations on the top of Brownian long-term trends. Statistically, this mixed clock represents an alternative solution to the problem of choosing between autocorrelated and uncorrelated relaxed clocks, by proposing instead to combine their respective merits. Fitting this model on a dataset of 105 placental mammals, using both node-dating and tip-dating approaches, suggests that the two pure clocks, Brownian and white noise, are rejected in favour of a mixed model with approximately equal contributions for its uncorrelated and autocorrelated components. The tip-dating analysis is particularly sensitive to the choice of the relaxed clock model. In this context, the classical pure Brownian relaxed clock appears to be overly rigid, leading to biases in divergence time estimation. By contrast, the use of a mixed clock leads to more recent and more reasonable estimates for the crown ages of placental orders and superorders. Altogether, the mixed clock introduced here represents a first step towards empirically more adequate models of the patterns of rate variation across phylogenetic trees.This article is part of the themed issue 'Dating species divergences using rocks and clocks'.

Darwinism for the Genomic Age: Connecting Mutation to Diversification

A growing body of evidence suggests that rates of diversification of biological lineages are correlated with differences in genome-wide mutation rate. Given that most research into differential patterns of diversification rate have focused on species traits or ecological parameters, a connection to the biochemical processes of genome change is an unexpected observation. While the empirical evidence for a significant association between mutation rate and diversification rate is mounting, there has been less effort in explaining the factors that mediate this connection between genetic change and species richness. Here we draw together empirical studies and theoretical concepts that may help to build links in the explanatory chain that connects mutation to diversification. First we consider the way that mutation rates vary between species. We then explore how differences in mutation rates have flow-through effects to the rate at which populations acquire substitutions, which in turn influences the speed at which populations become reproductively isolated from each other due to the acquisition of genomic incompatibilities. Since diversification rate is commonly measured from phylogenetic analyses, we propose a conceptual approach for relating events of reproductive isolation to bifurcations on molecular phylogenies. As we examine each of these relationships, we consider theoretical models that might shine a light on the observed association between rate of molecular evolution and diversification rate, and critically evaluate the empirical evidence for these links, focusing on phylogenetic comparative studies. Finally, we ask whether we are getting closer to a real understanding of the way that the processes of molecular evolution connect to the observable patterns of diversification.

Monte Carlo algorithms for Brownian phylogenetic models

MOTIVATION

Brownian models have been introduced in phylogenetics for describing variation in substitution rates through time, with applications to molecular dating or to the comparative analysis of variation in substitution patterns among lineages. Thus far, however, the Monte Carlo implementations of these models have relied on crude approximations, in which the Brownian process is sampled only at the internal nodes of the phylogeny or at the midpoints along each branch, and the unknown trajectory between these sampled points is summarized by simple branchwise average substitution rates.

RESULTS

A more accurate Monte Carlo approach is introduced, explicitly sampling a fine-grained discretization of the trajectory of the (potentially multivariate) Brownian process along the phylogeny. Generic Monte Carlo resampling algorithms are proposed for updating the Brownian paths along and across branches. Specific computational strategies are developed for efficient integration of the finite-time substitution probabilities across branches induced by the Brownian trajectory. The mixing properties and the computational complexity of the resulting Markov chain Monte Carlo sampler scale reasonably with the discretization level, allowing practical applications with up to a few hundred discretization points along the entire depth of the tree. The method can be generalized to other Markovian stochastic processes, making it possible to implement a wide range of time-dependent substitution models with well-controlled computational precision.

AVAILABILITY

The program is freely available at www.phylobayes.org.