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

Cladistic estimates of evolutionary rates focused on palaeontological datasets using TNT

We describe a protocol for estimating evolutionary rates from phylogenetic trees based on parsimony character optimization. The rate estimation is conducted through a TNT script and the results are analysed in a script for the software environment R. The TNT script allows analysing multiple optimal topologies, considering optimization ambiguity, and alternative time-calibrations or pre-calibrated trees. The R script summarizes estimated rates on a consensus tree and plots the variation of evolutionary rates through time, jointly with the phylogenetic diversity and a new metric (clade completeness index) that measures the distribution of missing data along the tree. We present results for simulated and empirical analyses, and evaluate the impact of missing data and alternative calibration methods in rate estimates. We found that while missing data can lower the nominal values of evolutionary rates, the overall pattern of rate variation through time remained robust. Empirical cases highlight different scenarios, such as datasets in which peaks of evolutionary rates can be coupled or decoupled from diversification dynamics (phylogenetic diversity) and cases in which missing data may influence the variation of estimated evolutionary rates. We conclude with recommendations for using this protocol and interpreting the results of parsimony-based rate estimates.

From cyanobacteria and cyanophages to chloroplasts: the fate of the genomes of oxyphototrophs and the genes encoding photosystem II proteins

Chloroplasts are the result of endosymbiosis of cyanobacterial organisms with proto-eukaryotes. The psbA, psbD and psbO genes are present in all oxyphototrophs and encode the D1/D2 proteins of photosystem II (PSII) and PsbO, respectively. PsbO is a peripheral protein that stabilizes the O2-evolving complex in PSII. Of these genes, psbA and psbD remained in the chloroplastic genome, while psbO was transferred to the nucleus. The genomes of selected cyanobacteria, chloroplasts and cyanophages carrying psbA and psbD, respectively, were analysed. The highest density of genes and coding sequences (CDSs) was estimated for the genomes of cyanophages, cyanobacteria and chloroplasts. The synonymous mutation rate (rS) of psbA and psbD in chloroplasts remained almost unchanged and is lower than that of psbO. The results indicate that the decreasing genome size in chloroplasts is more similar to the genome reduction observed in contemporary endosymbiotic organisms than in streamlined genomes of free-living cyanobacteria. The rS of atpA, which encodes the α-subunit of ATP synthase in chloroplasts, suggests that psbA and psbD, and to a lesser extent psbO, are ancient and conservative and arose early in the evolution of oxygenic photosynthesis. The role of cyanophages in the evolution of oxyphototrophs and chloroplastic genomes is discussed.

Scalable Bayesian Divergence Time Estimation With Ratio Transformations

Divergence time estimation is crucial to provide temporal signals for dating biologically important events from species divergence to viral transmissions in space and time. With the advent of high-throughput sequencing, recent Bayesian phylogenetic studies have analyzed hundreds to thousands of sequences. Such large-scale analyses challenge divergence time reconstruction by requiring inference on highly correlated internal node heights that often become computationally infeasible. To overcome this limitation, we explore a ratio transformation that maps the original $N-1$ internal node heights into a space of one height parameter and $N-2$ ratio parameters. To make the analyses scalable, we develop a collection of linear-time algorithms to compute the gradient and Jacobian-associated terms of the log-likelihood with respect to these ratios. We then apply Hamiltonian Monte Carlo sampling with the ratio transform in a Bayesian framework to learn the divergence times in 4 pathogenic viruses (West Nile virus, rabies virus, Lassa virus, and Ebola virus) and the coralline red algae. Our method both resolves a mixing issue in the West Nile virus example and improves inference efficiency by at least 5-fold for the Lassa and rabies virus examples as well as for the algae example. Our method now also makes it computationally feasible to incorporate mixed-effects molecular clock models for the Ebola virus example, confirms the findings from the original study, and reveals clearer multimodal distributions of the divergence times of some clades of interest.

Conceptual and empirical bridges between micro- and macroevolution

Explaining broad molecular, phenotypic and species biodiversity patterns necessitates a unifying framework spanning multiple evolutionary scales. Here we argue that although substantial effort has been made to reconcile microevolution and macroevolution, much work remains to identify the links between biological processes at play. We highlight four major questions of evolutionary biology whose solutions require conceptual bridges between micro and macroevolution. We review potential avenues for future research to establish how mechanisms at one scale (drift, mutation, migration, selection) translate to processes at the other scale (speciation, extinction, biogeographic dispersal) and vice versa. We propose ways in which current comparative methods to infer molecular evolution, phenotypic evolution and species diversification could be improved to specifically address these questions. We conclude that researchers are in a better position than ever before to build a synthesis to understand how microevolutionary dynamics unfold over millions of years.

Ancient dolphin genomes reveal rapid repeated adaptation to coastal waters

Parallel evolution provides strong evidence of adaptation by natural selection due to local environmental variation. Yet, the chronology, and mode of the process of parallel evolution remains debated. Here, we harness the temporal resolution of paleogenomics to address these long-standing questions, by comparing genomes originating from the mid-Holocene (8610-5626 years before present, BP) to contemporary pairs of coastal-pelagic ecotypes of bottlenose dolphin. We find that the affinity of ancient samples to coastal populations increases as the age of the samples decreases. We assess the youngest genome (5626 years BP) at sites previously inferred to be under parallel selection to coastal habitats and find it contained coastal-associated genotypes. Thus, coastal-associated variants rose to detectable frequencies close to the emergence of coastal habitat. Admixture graph analyses reveal a reticulate evolutionary history between pelagic and coastal populations, sharing standing genetic variation that facilitated rapid adaptation to newly emerged coastal habitats.

Survivors and colonizers: Contrasting biogeographic histories reconciled in the Antarctic freshwater copepod Boeckella poppei

Two main hypotheses have been proposed to explain the contemporary distribution of Antarctic terrestrial biota. We assess whether the current distribution of maritime Antarctic populations of the freshwater copepod Boeckella poppei is the result of (1) a post-Last Glacial Maximum (LGM) colonization, or whether (2) the species survived in regional glacial refugia throughout the LGM and earlier glaciations. Using 438 specimens from 34 different sampling sites across Southern South America, South Georgia, South Orkney Islands, South Shetland Islands, and the Antarctic Peninsula, we analyzed mitochondrial and nuclear sequences to uncover patterns of genetic diversity and population structure. We also performed median-joining haplotype network, phylogenetic reconstruction, and divergence time analyses. Finally, we evaluated past demographic changes and historical scenarios using the Approximate Bayesian Computation (ABC) method. Our data support the existence of two clades with different and contrasting biogeographic histories. The first clade has been present in maritime Antarctica since at least the mid-Pleistocene, with the South Orkney Islands the most likely refugial area. The second clade has a broader distribution including southern South America, South Georgia, South Shetland Islands, and the Antarctic Peninsula. The ABC method identified long-distance dispersal (LDD) colonization event(s) from southern South America to South Georgia and the maritime Antarctic after the LGM deglaciation, supporting more recent colonization of Antarctic locations. The current Antarctic and sub-Antarctic distribution of B. poppei is likely derived from two independent biogeographic events. The combination of both (1) post-LGM colonization from southern South America and (2) longer-term persistence in in situ regional refugia throughout glacial periods challenges current understanding of the biogeographic history of Antarctic freshwater biota. Re-colonization of ice-impacted Antarctic areas would have occurred following a LDD and Establishment model, pointing to the existence of possible post-dispersal barriers, despite widely assumed high passive dispersal capacity in freshwater invertebrates.

Detection of Ancient Viruses and Long-Term Viral Evolution

The COVID-19 outbreak has reminded us of the importance of viral evolutionary studies as regards comprehending complex viral evolution and preventing future pandemics. A unique approach to understanding viral evolution is the use of ancient viral genomes. Ancient viruses are detectable in various archaeological remains, including ancient people's skeletons and mummified tissues. Those specimens have preserved ancient viral DNA and RNA, which have been vigorously analyzed in the last few decades thanks to the development of sequencing technologies. Reconstructed ancient pathogenic viral genomes have been utilized to estimate the past pandemics of pathogenic viruses within the ancient human population and long-term evolutionary events. Recent studies revealed the existence of non-pathogenic viral genomes in ancient people's bodies. These ancient non-pathogenic viruses might be informative for inferring their relationships with ancient people's diets and lifestyles. Here, we reviewed the past and ongoing studies on ancient pathogenic and non-pathogenic viruses and the usage of ancient viral genomes to understand their long-term viral evolution.

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.

Causes and Consequences of Apparent Timescaling Across All Estimated Evolutionary Rates

Evolutionary rates play a central role in connecting micro- and macroevolution. All evolutionary rate estimates, including rates of molecular evolution, trait evolution, and lineage diversification, share a similar scaling pattern with time: The highest rates are those measured over the shortest time interval. This creates a disconnect between micro- and macroevolution, although the pattern is the opposite of what some might expect: Patterns of change over short timescales predict that evolution has tremendous potential to create variation and that potential is barely tapped by macroevolution. In this review, we discuss this shared scaling pattern across evolutionary rates. We break down possible explanations for scaling into two categories, estimation error and model misspecification, and discuss how both apply to each type of rate. We also discuss the consequences of this ubiquitous pattern, which can lead to unexpected results when comparing ratesover different timescales. Finally, after addressing purely statistical concerns, we explore a few possibilities for a shared unifying explanation across the three types of rates that results from a failure to fully understand and account for how biological processes scale over time.

Molecular archeology of human viruses

The evolution of human-virus associations is usually reconstructed from contemporary patterns of genomic diversity. An intriguing, though still rarely implemented, alternative is to search for the genetic material of viruses in archeological and medical archive specimens to document evolution as it happened. In this chapter, we present lessons from ancient DNA research and incorporate insights from virology to explore the potential range of applications and likely limitations of archeovirological approaches. We also highlight the numerous questions archeovirology will hopefully allow us to tackle in the near future, and the main expected roadblocks to these avenues of research.

Ape Origins of Human Malaria

African apes harbor at least twelve Plasmodium species, some of which have been a source of human infection. It is now well established that Plasmodium falciparum emerged following the transmission of a gorilla parasite, perhaps within the last 10,000 years, while Plasmodium vivax emerged earlier from a parasite lineage that infected humans and apes in Africa before the Duffy-negative mutation eliminated the parasite from humans there. Compared to their ape relatives, both human parasites have greatly reduced genetic diversity and an excess of nonsynonymous mutations, consistent with severe genetic bottlenecks followed by rapid population expansion. A putative new Plasmodium species widespread in chimpanzees, gorillas, and bonobos places the origin of Plasmodium malariae in Africa. Here, we review what is known about the origins and evolutionary history of all human-infective Plasmodium species, the time and circumstances of their emergence, and the diversity, host specificity, and zoonotic potential of their ape counterparts.

Genomes of Anguillid Herpesvirus 1 Strains Reveal Evolutionary Disparities and Low Genetic Diversity in the Genus Cyprinivirus

Anguillid herpesvirus 1 (AngHV-1) is a pathogen of eels and a member of the genus Cyprinivirus in the family Alloherpesviridae. We have compared the biological and genomic features of different AngHV-1 strains, focusing on their growth kinetics in vitro and genetic content, diversity, and recombination. Comparisons based on three core genes conserved among alloherpesviruses revealed that AngHV-1 exhibits a slower rate of change and less positive selection than other cypriniviruses. We propose that this may be linked to major differences in host species and corresponding epidemiological circumstances. Efforts to derive evolutionary rate estimates for cypriniviruses under various theoretical models were ultimately unrewarding. We highlight the potential value of future collaborative efforts towards generating short-term evolutionary rate estimates based on known sequence sampling dates. Finally, we revealed that there is significantly less genetic diversity in core gene sequences within cyprinivirus species clades compared to species in the family Herpesviridae. This suggests that cyprinivirus species may have undergone much more vigorous purifying selection post species clade divergence. We discuss whether this may be linked to biological and anthropogenic factors or to sampling bias, and we propose that the comparison of short-term evolutionary rates between species may provide further insights into these differences.

Paso doble: A two-step Late Pleistocene range expansion in the Tyrrhenian tree frog Hyla sarda

The Tyrrhenian tree frog, Hyla sarda, is an amphibian endemic to the Tyrrhenian islands (Western Mediterranean). Previous investigations of its Pleistocene evolutionary history suggested that it colonised the northern portion of its current range, through a spatial diffusion process from the Sardinia island, during the last glaciation. However, southern and northern portions of the species' range experienced markedly different climatic conditions during the Late Pleistocene, suggesting the possibility of an unusual two-step process of demographic expansion. Here, we use Bayesian phylogeographic approaches to locate the ancestral area in Sardinia and to characterise better the demographic component of this expansion event. These analyses located the ancestral area for H. sarda populations along the central-eastern coast of the Sardinia island, within an area previously shown to host suitable bioclimatic conditions for H. sarda populations throughout the Late Pleistocene. Historical demographic reconstructions clearly showed that a two-step process of demographic growth fits well the data, with northern populations expanding later than Sardinia populations. The harsher climatic conditions occurred in northern islands during the glacial epoch, as compared to Sardinia, likely delayed tree frog colonisation of northern territories, and the associated demographic growth.

Ancient Bacterial Genomes Reveal a High Diversity of Treponema pallidum Strains in Early Modern Europe

Syphilis is a globally re-emerging disease, which has marked European history with a devastating epidemic at the end of the 15^th century. Together with non-venereal treponemal diseases, like bejel and yaws, which are found today in subtropical and tropical regions, it currently poses a substantial health threat worldwide. The origins and spread of treponemal diseases remain unresolved, including syphilis' potential introduction into Europe from the Americas. Here, we present the first genetic data from archaeological human remains reflecting a high diversity of Treponema pallidum in early modern Europe. Our study demonstrates that a variety of strains related to both venereal syphilis and yaws-causing T. pallidum subspecies were already present in Northern Europe in the early modern period. We also discovered a previously unknown T. pallidum lineage recovered as a sister group to yaws- and bejel-causing lineages. These findings imply a more complex pattern of geographical distribution and etiology of early treponemal epidemics than previously understood.

Measles virus and rinderpest virus divergence dated to the sixth century BCE

Many infectious diseases are thought to have emerged in humans after the Neolithic revolution. Although it is broadly accepted that this also applies to measles, the exact date of emergence for this disease is controversial. We sequenced the genome of a 1912 measles virus and used selection-aware molecular clock modeling to determine the divergence date of measles virus and rinderpest virus. This divergence date represents the earliest possible date for the establishment of measles in human populations. Our analyses show that the measles virus potentially arose as early as the sixth century BCE, possibly coinciding with the rise of large cities.

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.

A practical introduction to sequentially Markovian coalescent methods for estimating demographic history from genomic data

A common goal of population genomics and molecular ecology is to reconstruct the demographic history of a species of interest. A pair of powerful tools based on the sequentially Markovian coalescent have been developed to infer past population sizes using genome sequences. These methods are most useful when sequences are available for only a limited number of genomes and when the aim is to study ancient demographic events. The results of these analyses can be difficult to interpret accurately, because doing so requires some understanding of their theoretical basis and of their sensitivity to confounding factors. In this practical review, we explain some of the key concepts underpinning the pairwise and multiple sequentially Markovian coalescent methods (PSMC and MSMC, respectively). We relate these concepts to the use and interpretation of these methods, and we explain how the choice of different parameter values by the user can affect the accuracy and precision of the inferences. Based on our survey of 100 PSMC studies and 30 MSMC studies, we describe how the two methods are used in practice. Readers of this article will become familiar with the principles, practice, and interpretation of the sequentially Markovian coalescent for inferring demographic history.

Mode and tempo of human hepatitis virus evolution

Human viral hepatitis, a major cause of morbidity and mortality worldwide, is caused by highly diverse viruses with different genetic, ecological, and pathogenetic features. Technological advances that allow throughput sequencing of viral genomes, as well as the development of computational tools to analyze such genome data, have largely expanded our knowledge on the host range and evolutionary history of human hepatitis viruses. Thus, with the exclusion of hepatitis D virus, close or distant relatives of these human pathogens were identified in a number of domestic and wild mammals. Also, sequences of human viral strains isolated from different geographic locations and over different time-spans have allowed the application of phylogeographic and molecular dating approaches to large viral phylogenies. In this review, we summarize the most recent insights into our understanding of the evolutionary events and ecological contexts that determined the origin and spread of human hepatitis viruses.

Viral quasispecies

Viral quasispecies refers to a population structure that consists of extremely large numbers of variant genomes, termed mutant spectra, mutant swarms or mutant clouds. Fueled by high mutation rates, mutants arise continually, and they change in relative frequency as viral replication proceeds. The term quasispecies was adopted from a theory of the origin of life in which primitive replicons) consisted of mutant distributions, as found experimentally with present day RNA viruses. The theory provided a new definition of wild type, and a conceptual framework for the interpretation of the adaptive potential of RNA viruses that contrasted with classical studies based on consensus sequences. Standard clonal analyses and deep sequencing methodologies have confirmed the presence of myriads of mutant genomes in viral populations, and their participation in adaptive processes. The quasispecies concept applies to any biological entity, but its impact is more evident when the genome size is limited and the mutation rate is high. This is the case of the RNA viruses, ubiquitous in our biosphere, and that comprise many important pathogens. In virology, quasispecies are defined as complex distributions of closely related variant genomes subjected to genetic variation, competition and selection, and that may act as a unit of selection. Despite being an integral part of their replication, high mutation rates have an upper limit compatible with inheritable information. Crossing such a limit leads to RNA virus extinction, a transition that is the basis of an antiviral design termed lethal mutagenesis.

Recalibrating the molecular clock for Arctic marine invertebrates based on DNA barcodes

Divergence times for species assemblages of Arctic marine invertebrates have often been estimated using a standard rate (1.4%/MY) of molecular evolution calibrated using a single sister pair of tropical crustaceans. Because rates of molecular evolution vary among taxa and environments, it is essential to obtain clock calibrations from northern lineages. The recurrent opening and closure of the Bering Strait provide an exceptional opportunity for clock calibration. Here, we apply the iterative calibration approach to investigate patterns of molecular divergence among lineages of northern marine molluscs and arthropods using publicly available sequences of the cytochrome c oxidase subunit I (COI) gene and compare these results with previous estimates of trans-Bering divergences for echinoderms and polychaetes. The wide range of Kimura two-parameter (K2P) divergences among 73 trans-Bering sister pairs (0.12%–16.89%) supports multiple pulses of migration through the Strait. Overall, the results indicate a rate of K2P divergence of 3.2%/MY in molluscs, 5%–5.2%/MY in arthropods, and 3.5%–4.7%/MY in polychaetes. While these rates are considerably higher than the often-adopted 1.4%/MY rate, they are similar to calibrations (3%–5%/MY) in several other studies of marine invertebrates. This upward revision in rates means there is a need both to reevaluate the evolutionary history of marine lineages and to reexamine the impact of prior climatic changes upon the diversification of marine life.

Modeling site-specific amino-acid preferences deepens phylogenetic estimates of viral sequence divergence

Molecular phylogenetics is often used to estimate the time since the divergence of modern gene sequences. For highly diverged sequences, such phylogenetic techniques sometimes estimate surprisingly recent divergence times. In the case of viruses, independent evidence indicates that the estimates of deep divergence times from molecular phylogenetics are sometimes too recent. This discrepancy is caused in part by inadequate models of purifying selection leading to branch-length underestimation. Here we examine the effect on branch-length estimation of using models that incorporate experimental measurements of purifying selection. We find that models informed by experimentally measured site-specific amino-acid preferences estimate longer deep branches on phylogenies of influenza virus hemagglutinin. This lengthening of branches is due to more realistic stationary states of the models, and is mostly independent of the branch-length extension from modeling site-to-site variation in amino-acid substitution rate. The branch-length extension from experimentally informed site-specific models is similar to that achieved by other approaches that allow the stationary state to vary across sites. However, the improvements from all of these site-specific but time homogeneous and site independent models are limited by the fact that a protein’s amino-acid preferences gradually shift as it evolves. Overall, our work underscores the importance of modeling site-specific amino-acid preferences when estimating deep divergence times—but also shows the inherent limitations of approaches that fail to account for how these preferences shift over time.

A Post-pleistocene Calibrated Mutation Rate from Insect Museum Specimens

Quantifying the age of recent species divergence events can be challenging in the absence of calibration points within many groups. The katydid species Neoconocephalus lyristes provides the opportunity to calibrate a post-Pleistocene, taxa specific mutation rate using a known biogeographic event, the Mohawk-Hudson Divide. DNA was extracted from pinned museum specimens of N. lyristes from both Midwest and Atlantic populations and the mitochondrial gene COI sequenced using primers designed from extant specimens. Coalescent analyses using both strict and relaxed molecular clock models were performed in BEAST v1.8.2. The assumption of a strict molecular clock could not be rejected in favor of the relaxed clock model as the distribution of the standard deviation of the clock rate strongly abutted zero. The strict molecular clock model resulted in an intraspecific calculated mutation rate of 14.4-17.3 %/myr, a rate substantially higher than the common rates of sequence evolution observed for insect mitochondrial DNA sequences. The rate, however, aligns closely with mutation rates estimated from other taxa with similarly recent lineage divergence times.

The Role of aDNA in Understanding the Coevolutionary Patterns of Human Sexually Transmitted Infections

Analysis of pathogen genome data sequenced from clinical and historical samples has made it possible to perform phylogenetic analyses of sexually transmitted infections on a global scale, and to estimate the diversity, distribution, and coevolutionary host relationships of these pathogens, providing insights into pathogen emergence and disease prevention. Deep-sequenced pathogen genomes from clinical studies and ancient samples yield estimates of within-host and between-host evolutionary rates and provide data on changes in pathogen genomic stability and evolutionary responses. Here we examine three groups of pathogens transmitted mainly through sexual contact between modern humans to provide insight into ancient human behavior and history with their pathogens. Exploring ancient pathogen genomic divergence and the ancient viral-host parallel evolutionary histories will help us to reconstruct the origin of present-day geographical distribution and diversity of clinical pathogen infections, and will hopefully allow us to foresee possible environmentally induced pathogen evolutionary responses. Lastly, we emphasize that ancient pathogen DNA research should be combined with modern clinical pathogen data, and be equitable and provide advantages for all researchers worldwide, e.g., through shared data.

Genetic structure and demographic inference of the regular sea urchin Sterechinus neumayeri (Meissner, 1900) in the Southern Ocean: The role of the last glaciation

One of the most relevant characteristics of the extant Southern Ocean fauna is its resiliency to survive glacial processes of the Quaternary. These climatic events produced catastrophic habitat reductions and forced some marine benthic species to move, adapt or go extinct. The marine benthic species inhabiting the Antarctic upper continental shelf faced the Quaternary glaciations with different strategies that drastically modified population sizes and thus affected the amount and distribution of intraspecific genetic variation. Here we present new genetic information for the most conspicuous regular sea urchin of the Antarctic continental shelf, Sterechinus neumayeri. We studied the patterns of genetic diversity and structure in this broadcast-spawner across three Antarctic regions: Antarctic Peninsula, the Weddell Sea and Adélie Land in East Antarctica. Genetic analyses based on mitochondrial and nuclear markers suggested that S. neumayeri is a single genetic unit around the Antarctic continent. The species is characterized by low levels of genetic diversity and exhibits a typical star-like haplotype genealogy that supports the hypothesis of a single in situ refugium. Based on two mutation rates standardized for this genus, the Bayesian Skyline plot analyses detected a rapid demographic expansion after the Last Glacial Maximum. We propose a scenario of rapid postglacial expansion and recolonization of Antarctic shallow areas from a less ice-impacted refugium where the species survived the LGM. Considering the patterns of genetic diversity and structure recorded in the species, this refugium was probably located in East Antarctica.

Evolutionary Analysis Provides Insight Into the Origin and Adaptation of HCV

Hepatitis C virus (HCV) belongs to the Hepacivirus genus and is genetically heterogeneous, with seven major genotypes further divided into several recognized subtypes. HCV origin was previously dated in a range between ∼200 and 1000 years ago. Hepaciviruses have been identified in several domestic and wild mammals, the largest viral diversity being observed in bats and rodents. The closest relatives of HCV were found in horses/donkeys (equine hepaciviruses, EHV). However, the origin of HCV as a human pathogen is still an unsolved puzzle. Using a selection-informed evolutionary model, we show that the common ancestor of extant HCV genotypes existed at least 3000 years ago (CI: 3192–5221 years ago), with the oldest genotypes being endemic to Asia. EHV originated around 1100 CE (CI: 291–1640 CE). These time estimates exclude that EHV transmission was mainly sustained by widespread veterinary practices and suggest that HCV originated from a single zoonotic event with subsequent diversification in human populations. We also describe a number of biologically important sites in the major HCV genotypes that have been positively selected and indicate that drug resistance-associated variants are significantly enriched at positively selected sites. HCV exploits several cell-surface molecules for cell entry, but only two of these (CD81 and OCLN) determine the species-specificity of infection. Herein evolutionary analyses do not support a long-standing association between primates and hepaciviruses, and signals of positive selection at CD81 were only observed in Chiroptera. No evidence of selection was detected for OCLN in any mammalian order. These results shed light on the origin of HCV and provide a catalog of candidate genetic modulators of HCV phenotypic diversity.

Impact of the tree prior on estimating clock rates during epidemic outbreaks

Genetic sequencing data of pathogens allow one to quantify the evolutionary rate together with epidemiological dynamics using Bayesian phylodynamic methods. Such tools are particularly useful for obtaining a timely understanding of newly emerging epidemic outbreaks. During the West African Ebola virus disease epidemic, an unusually high evolutionary rate was initially estimated, promoting discussions regarding the potential danger of the strain quickly evolving into an even more dangerous virus. We show here that such high evolutionary rates are not necessarily real but can stem from methodological biases in the analyses. While most analyses of epidemic outbreak data are performed such that these biases may be present, we suggest a solution to overcome these biases in the future.

Bayesian phylogenetics aims at estimating phylogenetic trees together with evolutionary and population dynamic parameters based on genetic sequences. It has been noted that the clock rate, one of the evolutionary parameters, decreases with an increase in the sampling period of sequences. In particular, clock rates of epidemic outbreaks are often estimated to be higher compared with the long-term clock rate. Purifying selection has been suggested as a biological factor that contributes to this phenomenon, since it purges slightly deleterious mutations from a population over time. However, other factors such as methodological biases may also play a role and make a biological interpretation of results difficult. In this paper, we identify methodological biases originating from the choice of tree prior, that is, the model specifying epidemiological dynamics. With a simulation study we demonstrate that a misspecification of the tree prior can upwardly bias the inferred clock rate and that the interplay of the different models involved in the inference can be complex and nonintuitive. We also show that the choice of tree prior can influence the inference of clock rate on real-world Ebola virus (EBOV) datasets. While commonly used tree priors result in very high clock-rate estimates for sequences from the initial phase of the epidemic in Sierra Leone, tree priors allowing for population structure lead to estimates agreeing with the long-term rate for EBOV.

Estimating evolutionary rates in giant viruses using ancient genomes

Pithovirus sibericum is a giant (610 Kpb) double-stranded DNA virus discovered in a purportedly 30,000-year-old permafrost sample. A closely related virus, Pithovirus massiliensis, was recently isolated from a sewer in southern France. An initial comparison of these two virus genomes assumed that P. sibericum was directly ancestral to P. massiliensis and gave a maximum evolutionary rate of 2.60 × 10^-5 nucleotide substitutions per site per year (subs/site/year). If correct, this would make pithoviruses among the fastest-evolving DNA viruses, with rates close to those seen in some RNA viruses. To help determine whether this unusually high rate is accurate we utilized the well-known negative association between evolutionary rate and genome size in DNA microbes. This revealed that a more plausible rate estimate for Pithovirus evolution is ∼2.23 × 10^-6 subs/site/year, with even lower estimates obtained if evolutionary rates are assumed to be time-dependent. Hence, we estimate that Pithovirus has evolved at least an order of magnitude more slowly than previously suggested. We then used our new rate estimates to infer a time-scale for Pithovirus evolution. Strikingly, this suggests that these viruses could have diverged at least hundreds of thousands of years ago, and hence have evolved over longer time-scales than previously suggested. We propose that the evolutionary rate and time-scale of pithovirus evolution should be reconsidered in the light of these observations and that future estimates of the rate of giant virus evolution should be carefully examined in the context of their biological plausibility.

Origin and dispersal of Hepatitis E virus

Hepatitis E virus (HEV, genus Orthohepevirus) is a common cause of hepatitis worldwide. Human-infecting HEV strains (Orthohepevirus A) include human-restricted and enzootic genotypes. Viruses in the Orthohepevirus A species also infect rabbits (HEV-3ra), camels, and swine. Using a selection-informed method, we dated the origin of the Orthohepevirus genus at least 21 million years ago, whereas the Orthohepevirus A species originated in Asia, most likely from a human-infecting ancestor that existed ~4500 to 6800 years ago. In this period, the appearance of large human settlements probably facilitated HEV emergence and spread. The earliest events in Orthohepevirus A evolutionary history involved the separation of the enzootic and human-restricted genotypes, as well as the split of the camel-infecting genotypes, which occurred during the time-frame of camel domestication. The place and timing of HEV-3ra divergence also correspond to the circumstances of rabbit domestication. This study clarifies the origin and historical events underlying HEV dispersal.

Pre- versus post-mass extinction divergence of Mesozoic marine reptiles dictated by time-scale dependence of evolutionary rates

The fossil record of a major clade often starts after a mass extinction even though evolutionary rates, molecular or morphological, suggest its pre-extinction emergence (e.g. squamates, placentals and teleosts). The discrepancy is larger for older clades, and the presence of a time-scale-dependent methodological bias has been suggested, yet it has been difficult to avoid the bias using Bayesian phylogenetic methods. This paradox raises the question of whether ecological vacancies, such as those after mass extinctions, prompt the radiations. We addressed this problem by using a unique temporal characteristic of the morphological data and a high-resolution stratigraphic record, for the oldest clade of Mesozoic marine reptiles, Ichthyosauromorpha. The evolutionary rate was fastest during the first few million years of ichthyosauromorph evolution and became progressively slower over time, eventually becoming six times slower. Using the later slower rates, estimates of divergence time become excessively older. The fast, initial rate suggests the emergence of ichthyosauromorphs after the end-Permian mass extinction, matching an independent result from high-resolution stratigraphic confidence intervals. These reptiles probably invaded the sea as a new ecosystem was formed after the end-Permian mass extinction. Lack of information on early evolution biased Bayesian clock rates.

Symbiosis in eukaryotic evolution

Fifty years ago, Lynn Margulis, inspiring in early twentieth-century ideas that put forward a symbiotic origin for some eukaryotic organelles, proposed a unified theory for the origin of the eukaryotic cell based on symbiosis as evolutionary mechanism. Margulis was profoundly aware of the importance of symbiosis in the natural microbial world and anticipated the evolutionary significance that integrated cooperative interactions might have as mechanism to increase cellular complexity. Today, we have started fully appreciating the vast extent of microbial diversity and the importance of syntrophic metabolic cooperation in natural ecosystems, especially in sediments and microbial mats. Also, not only the symbiogenetic origin of mitochondria and chloroplasts has been clearly demonstrated, but improvement in phylogenomic methods combined with recent discoveries of archaeal lineages more closely related to eukaryotes further support the symbiogenetic origin of the eukaryotic cell. Margulis left us in legacy the idea of 'eukaryogenesis by symbiogenesis'. Although this has been largely verified, when, where, and specifically how eukaryotic cells evolved are yet unclear. Here, we shortly review current knowledge about symbiotic interactions in the microbial world and their evolutionary impact, the status of eukaryogenetic models and the current challenges and perspectives ahead to reconstruct the evolutionary path to eukaryotes.

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.

Model misspecification confounds the estimation of rates and exaggerates their time dependency

While welcoming the comment of Ho et al. (2015), we find little that undermines the strength of our criticism, and it would appear they have misunderstood our central argument. Here we respond with the purpose of reiterating that we are (i) generally critical of much of the evidence presented in support of the time-dependent molecular rate (TDMR) hypothesis and (ii) specifically critical of estimates of μ derived from tip-dated sequences that exaggerate the importance of purifying selection as an explanation for TDMR over extended timescales. In response to assertions put forward by Ho et al. (2015), we use panmictic coalescent simulations of temporal data to explore a fundamental assumption for tip-dated tree shape and associated mutation rate estimates, and the appropriateness and utility of the date randomization test. The results reveal problems for the joint estimation of tree topology, effective population size and μ with tip-dated sequences using BEAST. Given the simulations, BEAST consistently obtains incorrect topological tree structures that are consistent with the substantial overestimation of μ and underestimation of effective population size. Data generated from lower effective population sizes were less likely to fail the date randomization test yet still resulted in substantially upwardly biased estimates of rates, bringing previous estimates of μ from temporally sampled DNA sequences into question. We find that our general criticisms of both the hypothesis of time-dependent molecular evolution and Bayesian methods to estimate μ from temporally sampled DNA sequences are further reinforced.