Rates of evolution in ancient DNA from Adélie penguins

Cooling perspectives on the risk of pathogenic viruses from thawing permafrost

Climate change is inducing wide-scale permafrost thaw in the Arctic and subarctic, triggering concerns that long-dormant pathogens could reemerge from the thawing ground and initiate epidemics or pandemics. Viruses, as opposed to bacterial pathogens, garner particular interest because outbreaks cannot be controlled with antibiotics, though the effects can be mitigated by vaccines and newer antiviral drugs. To evaluate the potential hazards posed by viral pathogens emerging from thawing permafrost, we review information from a diverse range of disciplines. This includes efforts to recover infectious virus from human remains, studies on disease occurrence in polar animal populations, investigations into viral persistence and infectivity in permafrost, and assessments of human exposure to the enormous viral diversity present in the environment. Based on currently available knowledge, we conclude that the risk posed by viruses from thawing permafrost is no greater than viruses in other environments such as temperate soils and aquatic systems.

Ancient and Modern Genomes Reveal Microsatellites Maintain a Dynamic Equilibrium Through Deep Time

Microsatellites are widely used in population genetics, but their evolutionary dynamics remain poorly understood. It is unclear whether microsatellite loci drift in length over time. This is important because the mutation processes that underlie these important genetic markers are central to the evolutionary models that employ microsatellites. We identify more than 27 million microsatellites using a novel and unique dataset of modern and ancient Adélie penguin genomes along with data from 63 published chordate genomes. We investigate microsatellite evolutionary dynamics over 2 timescales: one based on Adélie penguin samples dating to ∼46.5 ka and the other dating to the diversification of chordates aged more than 500 Ma. We show that the process of microsatellite allele length evolution is at dynamic equilibrium; while there is length polymorphism among individuals, the length distribution for a given locus remains stable. Many microsatellites persist over very long timescales, particularly in exons and regulatory sequences. These often retain length variability, suggesting that they may play a role in maintaining phenotypic variation within populations.

Wild pedigrees inform mutation rates and historic abundance in baleen whales

Phylogeny-based estimates suggesting a low germline mutation rate (μ) in baleen whales have influenced research ranging from assessments of whaling impacts to evolutionary cancer biology. We estimated μ directly from pedigrees in four baleen whale species for both the mitochondrial control region and nuclear genome. The results suggest values higher than those obtained through phylogeny-based estimates and similar to pedigree-based values for primates and toothed whales. Applying our estimate of μ reduces previous genetic-based estimates of preexploitation whale abundance by 86% and suggests that μ cannot explain low cancer rates in gigantic mammals. Our study shows that it is feasible to estimate μ directly from pedigrees in natural populations, with wide-ranging implications for ecological and evolutionary research.

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.

Temporal dynamics of mildly deleterious nonsynonymous substitutions in mitochondrial gene sequences in rodents and moles

We have previously estimated the evolutionary rate (number of substitutions/site/million years) of mitochondrial cytochrome b gene (Cytb) sequences in rodents and moles to be about 0.11 at more recent divergence times of tens of thousands of years, and to decrease rapidly to about 0.03 at more distant divergence times. Because this time dependency is thought to be caused by the removal of mildly deleterious substitutions in later generations, we focused in this study on the abundance of nonsynonymous substitutions. We collected 23 haplogroups of Cytb with signals of late Quaternary population expansion events from rodents and moles and categorized them into three groups for comparison based on predicted expansion start time: 5,000-15,000 years ago (Group I), ca. 53,000 years ago (Group II) and 130,000-230,000 years ago (Group III). We counted the numbers of nonsynonymous and synonymous substitutions in all haplogroups. The rates of nonsynonymous substitutions were lowest in Groups II and III (0.08-0.22), whereas those in Group I varied markedly. We further classified Group I into two subgroups based on high (0.29-0.43) and low (0.09-0.20) nonsynonymous substitution rates, which were likely to be associated with the start of the expansion within 10,000 years and at around 15,000 years ago, respectively. The Group II and III networks had two- or three-step star-shaped structures and tended to exhibit frequent and less frequent nonsynonymous substitutions on exterior and interior branches, respectively. Based on temporal dynamics, nonsynonymous mitochondrial DNA (mtDNA) substitutions in small mammals accounted for at most 40% of all substitutions during the early evolutionary stage and then rapidly declined, dropping to approximately 15%. The results of this study provide a good explanation of the time-dependent trend in the mtDNA evolution rate predicted in previous work.

Biomolecular analyses reveal the age, sex and species identity of a near-intact Pleistocene bird carcass

Ancient remains found in permafrost represent a rare opportunity to study past ecosystems. Here, we present an exceptionally well-preserved ancient bird carcass found in the Siberian permafrost, along with a radiocarbon date and a reconstruction of its complete mitochondrial genome. The carcass was radiocarbon dated to approximately 44–49 ka BP, and was genetically identified as a female horned lark. This is a species that usually inhabits open habitat, such as the steppe environment that existed in Siberia at the time. This near-intact carcass highlights the potential of permafrost remains for evolutionary studies that combine both morphology and ancient nucleic acids.

Nicolas Dussex et al. identify a 44,000–49,000 year old bird found in Siberian permafrost as a female horned lark using ancient DNA. This exceptionally well-preserved specimen illustrates the potential contribution to science of permafrost deposits, such as the study of ecology and evolution of ancient ecosystems, calibration of molecular clocks, and furthering our understanding of processes such as biological regulation and gene expression in relation to climate change.

Unlocking the potential of ancient fish DNA in the genomic era

Fish are the most diverse group of vertebrates, fulfil important ecological functions and are of significant economic interest for aquaculture and wild fisheries. Advances in DNA extraction methods, sequencing technologies and bioinformatic applications have advanced genomic research for nonmodel organisms, allowing the field of fish ancient DNA (aDNA) to move into the genomics era. This move is enabling researchers to investigate a multitude of new questions in evolutionary ecology that could not, until now, be addressed. In many cases, these new fields of research have relevance to evolutionary applications, such as the sustainable management of fisheries resources and the conservation of aquatic animals. Here, we focus on the application of fish aDNA to (a) highlight new research questions, (b) outline methodological advances and current challenges, (c) discuss how our understanding of fish ecology and evolution can benefit from aDNA applications and (d) provide a future perspective on how the field will help answer key questions in conservation and management. We conclude that the power of fish aDNA will be unlocked through the application of continually improving genomic resources and methods to well-chosen taxonomic groups represented by well-dated archaeological samples that can provide temporally and/or spatially extensive data sets.

Insights into the molecular phylogeny of Rhaphidophoridae, an ancient, worldwide lineage of Orthoptera

We investigated the molecular phylogenetic divergence and historical biogeography of cave crickets belonging to the family Rhaphidophoridae (Orthoptera, Ensifera). We used taxa representative of most of the regions embraced by the family, considering samples of Macropathinae from Gondwana land (i.e., Tasmania, Australia, New Zealand, South Africa, and South America); Aemodogryllinae and Rhaphidophorinae from Southern-eastern Asia (i.e., India, Bhutan, China, Philippines and the Sulawesi islands); Dolichopodainae and Troglophilinae from the Mediterranean region and Ceuthophilinae from North America. Based on previous papers, we carried out an analysis of both mitochondrial and nuclear DNA sequences considering the ribosomal RNA units 12S, 16S, 18S, and 28S. To reconstruct phylogeny, we use cladistics, Maximum Likelihood (ML), and Bayesian analyses. All phylogenetic analyses showed the same highly supported topology generally congruent with the classical systematic arrangement at the level of each sub-family but strongly disagree with previous affinity hypotheses between sub-families based on morphological characters. Our results reveal a close affinity between Asiatic and Gondwanian taxa from one hand and between North American and Mediterranean ones from the other hand. Dating estimates indicated that Rhaphidophoridae originated in the Cretaceous period during the Mesozoic era with the ancestral area located both in the northern and southern hemisphere. A possible biogeographic scenario, reconstructed using S-DEC with RASP software, suggested that the current distribution of Rhaphidophoridae might be explained by a combination of both dispersal and vicariance events occurred especially in the ancestral populations. The radiation of Rhaphidophoridae started within the Pangaea, where the ancestor of Rhaphidophoridae occurred throughout an ancestral area including Australia, North America, and the Mediterranean region. The opening of the Atlantic Ocean promoted the divergence of North American and Mediterranean lineages while the differentiation of the southern lineages, spread from Australia, appears to be related to the fragmentation of Gondwana land.

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.

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.

Comparative genome-wide polymorphic microsatellite markers in Antarctic penguins through next generation sequencing

Microsatellites are valuable molecular markers for evolutionary and ecological studies. Next generation sequencing is responsible for the increasing number of microsatellites for non-model species. Penguins of the Pygoscelis genus are comprised of three species: Adélie (P. adeliae), Chinstrap (P. antarcticus) and Gentoo penguin (P. papua), all distributed around Antarctica and the sub-Antarctic. The species have been affected differently by climate change, and the use of microsatellite markers will be crucial to monitor population dynamics. We characterized a large set of genome-wide microsatellites and evaluated polymorphisms in all three species. SOLiD reads were generated from the libraries of each species, identifying a large amount of microsatellite loci: 33,677, 35,265 and 42,057 for P. adeliae, P. antarcticus and P. papua, respectively. A large number of dinucleotide (66,139), trinucleotide (29,490) and tetranucleotide (11,849) microsatellites are described. Microsatellite abundance, diversity and orthology were characterized in penguin genomes. We evaluated polymorphisms in 170 tetranucleotide loci, obtaining 34 polymorphic loci in at least one species and 15 polymorphic loci in all three species, which allow to perform comparative studies. Polymorphic markers presented here enable a number of ecological, population, individual identification, parentage and evolutionary studies of Pygoscelis, with potential use in other penguin species.

A latitudinal phylogeographic diversity gradient in birds

High tropical species diversity is often attributed to evolutionary dynamics over long timescales. It is possible, however, that latitudinal variation in diversification begins when divergence occurs within species. Phylogeographic data capture this initial stage of diversification in which populations become geographically isolated and begin to differentiate genetically. There is limited understanding of the broader implications of intraspecific diversification because comparative analyses have focused on species inhabiting and evolving in restricted regions and environments. Here, we scale comparative phylogeography up to the hemisphere level and examine whether the processes driving latitudinal differences in species diversity are also evident within species. We collected genetic data for 210 New World bird species distributed across a broad latitudinal gradient and estimated a suite of metrics characterizing phylogeographic history. We found that lower latitude species had, on average, greater phylogeographic diversity than higher latitude species and that intraspecific diversity showed evidence of greater persistence in the tropics. Factors associated with species ecologies, life histories, and habitats explained little of the variation in phylogeographic structure across the latitudinal gradient. Our results suggest that the latitudinal gradient in species richness originates, at least partly, from population-level processes within species and are consistent with hypotheses implicating age and environmental stability in the formation of diversity gradients. Comparative phylogeographic analyses scaled up to large geographic regions and hundreds of species can show connections between population-level processes and broad-scale species-richness patterns.

Invader or resident? Ancient-DNA reveals rapid species turnover in New Zealand little penguins

The expansion of humans into previously unoccupied parts of the globe is thought to have driven the decline and extinction of numerous vertebrate species. In New Zealand, human settlement in the late thirteenth century AD led to the rapid demise of a distinctive vertebrate fauna, and also a number of 'turnover' events where extinct lineages were subsequently replaced by closely related taxa. The recent genetic detection of an Australian little penguin (Eudyptula novaehollandiae) in southeastern New Zealand may potentially represent an additional 'cryptic' invasion. Here we use ancient-DNA (aDNA) analysis and radiocarbon dating of pre-human, archaeological and historical Eudyptula remains to reveal that the arrival of E. novaehollandiae in New Zealand probably occurred between AD 1500 and 1900, following the anthropogenic decline of its sister taxon, the endemic Eudyptula minor. This rapid turnover event, revealed by aDNA, suggests that native species decline can be masked by invasive taxa, and highlights the potential for human-mediated biodiversity shifts.

Post-glacial habitat release and incipient speciation in the genus Delphinus

The role of ecological and changing environmental factors in the radiation of species diversity is a fundamental question in evolutionary biology. Of particular interest is the potential for these factors to determine the boundary between what we would consider differentiation among populations and incipient speciation. Dolphins in the genus Delphinus provide a useful test case, exhibiting morphological variation in beak length, coloration and body size across their wide geographic distribution, and in particular among coastal and more pelagic habitats. Two species have been proposed, D. delphis and D. capensis, but morphologically similar allopatric populations are not monophyletic, indicating that the mostly coastal 'long-beaked' D. capensis form is not a single globally distributed species. However, the sympatric populations in the Eastern North Pacific currently designated as these two species are both morphologically and genetically differentiated. Here we use microsatellite DNA and mitochondrial DNA markers to investigate the evolutionary mechanisms that led to this incipient speciation event. We used coalescent and assignment methods to investigate the timing and extent of reproductive isolation. Our data indicate that although there is some level of on-going gene flow, the putative species found in the Eastern North Pacific are reciprocally monophyletic. The timing of isolation appears to be associated with regional changes in paleoceanographic conditions within the Holocene timeframe.

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.

Sequence capture by hybridization to explore modern and ancient genomic diversity in model and nonmodel organisms

The recent expansion of next-generation sequencing has significantly improved biological research. Nevertheless, deep exploration of genomes or metagenomic samples remains difficult because of the sequencing depth and the associated costs required. Therefore, different partitioning strategies have been developed to sequence informative subsets of studied genomes. Among these strategies, hybridization capture has proven to be an innovative and efficient tool for targeting and enriching specific biomarkers in complex DNA mixtures. It has been successfully applied in numerous areas of biology, such as exome resequencing for the identification of mutations underlying Mendelian or complex diseases and cancers, and its usefulness has been demonstrated in the agronomic field through the linking of genetic variants to agricultural phenotypic traits of interest. Moreover, hybridization capture has provided access to underexplored, but relevant fractions of genomes through its ability to enrich defined targets and their flanking regions. Finally, on the basis of restricted genomic information, this method has also allowed the expansion of knowledge of nonreference species and ancient genomes and provided a better understanding of metagenomic samples. In this review, we present the major advances and discoveries permitted by hybridization capture and highlight the potency of this approach in all areas of biology.

The influence of historical climate changes on Southern Ocean marine predator populations: a comparative analysis

The Southern Ocean ecosystem is undergoing rapid physical and biological changes that are likely to have profound implications for higher-order predators. Here, we compare the long-term, historical responses of Southern Ocean predators to climate change. We examine palaeoecological evidence for changes in the abundance and distribution of seabirds and marine mammals, and place these into context with palaeoclimate records in order to identify key environmental drivers associated with population changes. Our synthesis revealed two key factors underlying Southern Ocean predator population changes; (i) the availability of ice-free ground for breeding and (ii) access to productive foraging grounds. The processes of glaciation and sea ice fluctuation were key; the distributions and abundances of elephant seals, snow petrels, gentoo, chinstrap and Adélie penguins all responded strongly to the emergence of new breeding habitat coincident with deglaciation and reductions in sea ice. Access to productive foraging grounds was another limiting factor, with snow petrels, king and emperor penguins all affected by reduced prey availability in the past. Several species were isolated in glacial refugia and there is evidence that refuge populations were supported by polynyas. While the underlying drivers of population change were similar across most Southern Ocean predators, the individual responses of species to environmental change varied because of species specific factors such as dispersal ability and environmental sensitivity. Such interspecific differences are likely to affect the future climate change responses of Southern Ocean marine predators and should be considered in conservation plans. Comparative palaeoecological studies are a valuable source of long-term data on species' responses to environmental change that can provide important insights into future climate change responses. This synthesis highlights the importance of protecting productive foraging grounds proximate to breeding locations, as well as the potential role of polynyas as future Southern Ocean refugia.

Proliferation of East Antarctic Adélie penguins in response to historical deglaciation

Background

Major, long-term environmental changes are projected in the Southern Ocean and these are likely to have impacts for marine predators such as the Adélie penguin (Pygoscelis adeliae). Decadal monitoring studies have provided insight into the short-term environmental sensitivities of Adélie penguin populations, particularly to sea ice changes. However, given the long-term nature of projected climate change, it is also prudent to consider the responses of populations to environmental change over longer time scales. We investigated the population trajectory of Adélie penguins during the last glacial-interglacial transition to determine how the species was affected by climate warming over millennia. We focussed our study on East Antarctica, which is home to 30 % of the global population of Adélie penguins.

Methods

Using mitochondrial DNA from extant colonies, we reconstructed the population trend of Adélie penguins in East Antarctica over the past 22,000 years using an extended Bayesian skyline plot method. To determine the relationship of East Antarctic Adélie penguins with populations elsewhere in Antarctica we constructed a phylogeny using mitochondrial DNA sequences.

Results

We found that the Adélie penguin population expanded 135-fold from approximately 14,000 years ago. The population growth was coincident with deglaciation in East Antarctica and, therefore, an increase in ice-free ground suitable for Adélie penguin nesting. Our phylogenetic analysis indicated that East Antarctic Adélie penguins share a common ancestor with Adélie penguins from the Antarctic Peninsula and Scotia Arc, with an estimated age of 29,000 years ago, in the midst of the last glacial period. This finding suggests that extant colonies in East Antarctica, the Scotia Arc and the Antarctic Peninsula were founded from a single glacial refuge.

Conclusions

While changes in sea ice conditions are a critical driver of Adélie penguin population success over decadal and yearly timescales, deglaciation appears to have been the key driver of population change over millennia. This suggests that environmental drivers of population trends over thousands of years may differ to drivers over years or decades, highlighting the need to consider millennial-scale trends alongside contemporary data for the forecasting of species’ abundance and distribution changes under future climate change scenarios.

Fast Dating Using Least-Squares Criteria and Algorithms

Phylogenies provide a useful way to understand the evolutionary history of genetic samples, and data sets with more than a thousand taxa are becoming increasingly common, notably with viruses (e.g., human immunodeficiency virus (HIV)). Dating ancestral events is one of the first, essential goals with such data. However, current sophisticated probabilistic approaches struggle to handle data sets of this size. Here, we present very fast dating algorithms, based on a Gaussian model closely related to the Langley–Fitch molecular-clock model. We show that this model is robust to uncorrelated violations of the molecular clock. Our algorithms apply to serial data, where the tips of the tree have been sampled through times. They estimate the substitution rate and the dates of all ancestral nodes. When the input tree is unrooted, they can provide an estimate for the root position, thus representing a new, practical alternative to the standard rooting methods (e.g., midpoint). Our algorithms exploit the tree (recursive) structure of the problem at hand, and the close relationships between least-squares and linear algebra. We distinguish between an unconstrained setting and the case where the temporal precedence constraint (i.e., an ancestral node must be older that its daughter nodes) is accounted for. With rooted trees, the former is solved using linear algebra in linear computing time (i.e., proportional to the number of taxa), while the resolution of the latter, constrained setting, is based on an active-set method that runs in nearly linear time. With unrooted trees the computing time becomes (nearly) quadratic (i.e., proportional to the square of the number of taxa). In all cases, very large input trees (>10,000 taxa) can easily be processed and transformed into time-scaled trees. We compare these algorithms to standard methods (root-to-tip, r8s version of Langley–Fitch method, and BEAST). Using simulated data, we show that their estimation accuracy is similar to that of the most sophisticated methods, while their computing time is much faster. We apply these algorithms on a large data set comprising 1194 strains of Influenza virus from the pdm09 H1N1 Human pandemic. Again the results show that these algorithms provide a very fast alternative with results similar to those of other computer programs. These algorithms are implemented in the LSD software (least-squares dating), which can be downloaded from http://www.atgc-montpellier.fr/LSD/, along with all our data sets and detailed results. An Online Appendix, providing additional algorithm descriptions, tables, and figures can be found in the Supplementary Material available on Dryad at http://dx.doi.org/10.5061/dryad.968t3.

Ancient population genomics and the study of evolution

Recently, the study of ancient DNA (aDNA) has been greatly enhanced by the development of second-generation DNA sequencing technologies and targeted enrichment strategies. These developments have allowed the recovery of several complete ancient genomes, a result that would have been considered virtually impossible only a decade ago. Prior to these developments, aDNA research was largely focused on the recovery of short DNA sequences and their use in the study of phylogenetic relationships, molecular rates, species identification and population structure. However, it is now possible to sequence a large number of modern and ancient complete genomes from a single species and thereby study the genomic patterns of evolutionary change over time. Such a study would herald the beginnings of ancient population genomics and its use in the study of evolution. Species that are amenable to such large-scale studies warrant increased research effort. We report here progress on a population genomic study of the Adélie penguin (Pygoscelis adeliae). This species is ideally suited to ancient population genomic research because both modern and ancient samples are abundant in the permafrost conditions of Antarctica. This species will enable us to directly address many of the fundamental questions in ecology and evolution.

Mercury Stable Isotopes in Ornithogenic Deposits As Tracers of Historical Cycling of Mercury in Ross Sea, Antarctica

Production of methylmercury (MeHg) in ocean waters and its bioaccumulation in marine organisms are critical processes controlling the fate and toxicity of mercury (Hg). However, these processes are not well understood in the Antarctic, where high levels of MeHg are observed in the subsurface ocean (100-1000 m). We explored the use of Hg stable isotope compositions in historical and modern biological deposits as a new approach for discerning Hg sources and tracing MeHg cycling in the ocean and bioaccumulation in marine biota. We found similar mass independent isotope fractionation (MIF) of Hg between a sediment profile containing historical penguin and seal feces deposits from coastal Antarctica and modern penguin and seal feces, suggesting that penguin and seal feces were the dominant sources of Hg to the sediments at different time periods. Furthermore, sediments dominated by seal feces displayed a significantly lower MIF slope (Δ(199)Hg/Δ(201)Hg) than those dominated by penguin feces despite similar extents of MIF. Since seals forage at greater depths (>400 m) than penguins (<100 m), the high MIF values and lower Δ(199)Hg/Δ(201)Hg in seal feces suggest that a significant fraction of MeHg accumulated by seals was produced in situ in the subsurface ocean from residual inorganic Hg(II) that sank from the euphotic zone after partial photoreduction. Our results suggest that in situ Hg methylation can be an important source of MeHg for marine biota, and Hg isotope compositions in biological archives can be valuable tracers of MeHg cycling.

Prolonged decay of molecular rate estimates for metazoan mitochondrial DNA

Evolutionary timescales can be estimated from genetic data using the molecular clock, often calibrated by fossil or geological evidence. However, estimates of molecular rates in mitochondrial DNA appear to scale negatively with the age of the clock calibration. Although such a pattern has been observed in a limited range of data sets, it has not been studied on a large scale in metazoans. In addition, there is uncertainty over the temporal extent of the time-dependent pattern in rate estimates. Here we present a meta-analysis of 239 rate estimates from metazoans, representing a range of timescales and taxonomic groups. We found evidence of time-dependent rates in both coding and non-coding mitochondrial markers, in every group of animals that we studied. The negative relationship between the estimated rate and time persisted across a much wider range of calibration times than previously suggested. This indicates that, over long time frames, purifying selection gives way to mutational saturation as the main driver of time-dependent biases in rate estimates. The results of our study stress the importance of accounting for time-dependent biases in estimating mitochondrial rates regardless of the timescale over which they are inferred.

Major transitions in human evolution revisited: a tribute to ancient DNA

The origin and diversification of modern humans have been characterized by major evolutionary transitions and demographic changes. Patterns of genetic variation within modern populations can help with reconstructing this ∼200 thousand year-long population history. However, by combining this information with genomic data from ancient remains, one can now directly access our evolutionary past and reveal our population history in much greater detail. This review outlines the main recent achievements in ancient DNA research and illustrates how the field recently moved from the polymerase chain reaction (PCR) amplification of short mitochondrial fragments to whole-genome sequencing and thereby revisited our own history. Ancient DNA research has revealed the routes that our ancestors took when colonizing the planet, whom they admixed with, how they domesticated plant and animal species, how they genetically responded to changes in lifestyle, and also, which pathogens decimated their populations. These approaches promise to soon solve many pending controversies about our own origins that are indecipherable from modern patterns of genetic variation alone, and therefore provide an extremely powerful toolkit for a new generation of molecular anthropologists.

The future of ancient DNA: Technical advances and conceptual shifts

Technological innovations such as next generation sequencing and DNA hybridisation enrichment have resulted in multi-fold increases in both the quantity of ancient DNA sequence data and the time depth for DNA retrieval. To date, over 30 ancient genomes have been sequenced, moving from 0.7× coverage (mammoth) in 2008 to more than 50× coverage (Neanderthal) in 2014. Studies of rapid evolutionary changes, such as the evolution and spread of pathogens and the genetic responses of hosts, or the genetics of domestication and climatic adaptation, are developing swiftly and the importance of palaeogenomics for investigating evolutionary processes during the last million years is likely to increase considerably. However, these new datasets require new methods of data processing and analysis, as well as conceptual changes in interpreting the results. In this review we highlight important areas of future technical and conceptual progress and discuss research topics in the rapidly growing field of palaeogenomics.

Selective constraints determine the time dependency of molecular rates for human nuclear genomes

In contrast to molecular rates for neutral mitochondrial sequences, rates for constrained sites (including nonsynonymous sites, D-loop, and RNA) in the mitochondrial genome are known to vary with the time frame used for their estimation. Here, we examined this issue for the nuclear genomes using single-nucleotide polymorphisms (SNPs) from six complete human genomes of individuals belonging to different populations. We observed a strong time-dependent distribution of nonsynonymous SNPs (nSNPs) in highly constrained genes. Typically, the proportion of young nSNPs specific to a single population was found to be up to three times higher than that of the ancient nSNPs shared between diverse human populations. In contrast, this trend disappeared, and a uniform distribution of young and old nSNPs was observed in genes under relaxed selective constraints. This suggests that because mutations in constrained genes are highly deleterious, they are removed over time, resulting in a relative overabundance of young nSNPs. In contrast, mutations in genes under relaxed constraints are nearly neutral, which leads to similar proportions of young and old SNPs. These results could be useful to researchers aiming to select appropriate genes or genomic regions for estimating evolutionary rates and species or population divergence times.

The de novo assembly of mitochondrial genomes of the extinct passenger pigeon (Ectopistes migratorius) with next generation sequencing

The information from ancient DNA (aDNA) provides an unparalleled opportunity to infer phylogenetic relationships and population history of extinct species and to investigate genetic evolution directly. However, the degraded and fragmented nature of aDNA has posed technical challenges for studies based on conventional PCR amplification. In this study, we present an approach based on next generation sequencing to efficiently sequence the complete mitochondrial genome (mitogenome) of two extinct passenger pigeons (Ectopistes migratorius) using de novo assembly of massive short (90 bp), paired-end or single-end reads. Although varying levels of human contamination and low levels of postmortem nucleotide lesion were observed, they did not impact sequencing accuracy. Our results demonstrated that the de novo assembly of shotgun sequence reads could be a potent approach to sequence mitogenomes, and offered an efficient way to infer evolutionary history of extinct species.

Ancient DNA studies: new perspectives on old samples

In spite of past controversies, the field of ancient DNA is now a reliable research area due to recent methodological improvements. A series of recent large-scale studies have revealed the true potential of ancient DNA samples to study the processes of evolution and to test models and assumptions commonly used to reconstruct patterns of evolution and to analyze population genetics and palaeoecological changes. Recent advances in DNA technologies, such as next-generation sequencing make it possible to recover DNA information from archaeological and paleontological remains allowing us to go back in time and study the genetic relationships between extinct organisms and their contemporary relatives. With the next-generation sequencing methodologies, DNA sequences can be retrieved even from samples (for example human remains) for which the technical pitfalls of classical methodologies required stringent criteria to guaranty the reliability of the results. In this paper, we review the methodologies applied to ancient DNA analysis and the perspectives that next-generation sequencing applications provide in this field.

Phylogeny and evolutionary history of the silkworm

The silkworm, Bombyx mori, played an important role in the old Silk Road that connected ancient Asia and Europe. However, to date, there have been few studies of the origins and domestication of this species using molecular methods. In this study, DNA sequences of mitochondrial and nuclear loci were used to infer the phylogeny and evolutionary history of the domesticated silkworm and its relatives. All of the phylogenetic analyses indicated a close relationship between the domesticated silkworm and the Chinese wild silkworm. Domestication was estimated to have occurred about 4100 years ago (ya), and the radiation of the different geographic strains of B. mori about 2000 ya. The Chinese wild silkworm and the Japanese wild silkworm split about 23600 ya. These estimates are in good agreement with the fossil evidence and historical records. In addition, we show that the domesticated silkworm experienced a population expansion around 1000 ya. The divergence times and the population dynamics of silkworms presented in this study will be useful for studies of lepidopteran phylogenetics, in the genetic analysis of domestic animals, and for understanding the spread of human civilizations.

Adélie penguins and temperature changes in Antarctica: a long-term view

During the summer months, Adélie penguins represent the dominant biomass of terrestrial Antarctica. Literally millions of individuals nest in ice-free areas around the coast of the continent. Hence, these modern populations of Adélie penguins have often been championed as an ideal biological indicator of ecological and environmental changes that we currently face. In addition, Adélie penguins show an extraordinary record of sub-fossil remains, dating back to the late Pleistocene. At this time, temperatures were much lower than now. Hence, this species offers unique long-term information, at both the genomic and ecological levels, about how a species has responded to climate change over more than 40 000 years.

Resurrecting ancient animal genomes: the extinct moa and more

Recently two developments have had a major impact on the field of ancient DNA (aDNA). First, new advances in DNA sequencing, in combination with improved capture/enrichment methods, have resulted in the recovery of orders of magnitude more DNA sequence data from ancient animals. Second, there has been an increase in the range of tissue types employed in aDNA. Hair in particular has proven to be very successful as a source of DNA because of its low levels of contamination and high level of ancient endogenous DNA. These developments have resulted in significant advances in our understanding of recently extinct animals: namely their evolutionary relationships, physiology, and even behaviour. Hair has been used to recover the first complete ancient nuclear genome, that of the extinct woolly mammoth, which then facilitated the expression and functional analysis of haemoglobins. Finally, we speculate on the consequences of these developments for the possibility of recreating extinct animals.

Phylogenetic analysis of four nuclear protein-encoding genes largely corroborates the traditional classification of Bivalvia (Mollusca)

Revived interest in molluscan phylogeny has resulted in a torrent of molecular sequence data from phylogenetic, mitogenomic, and phylogenomic studies. Despite recent progress, basal relationships of the class Bivalvia remain contentious, owing to conflicting morphological and molecular hypotheses. Marked incongruity of phylogenetic signal in datasets heavily represented by nuclear ribosomal genes versus mitochondrial genes has also impeded consensus on the type of molecular data best suited for investigating bivalve relationships. To arbitrate conflicting phylogenetic hypotheses, we evaluated the utility of four nuclear protein-encoding genes-ATP synthase β, elongation factor-1α, myosin heavy chain type II, and RNA polymerase II-for resolving the basal relationships of Bivalvia. We sampled all five major lineages of bivalves (Archiheterodonta, Euheterodonta [including Anomalodesmata], Palaeoheterodonta, Protobranchia, and Pteriomorphia) and inferred relationships using maximum likelihood and Bayesian approaches. To investigate the robustness of the phylogenetic signal embedded in the data, we implemented additional datasets wherein length variability and/or third codon positions were eliminated. Results obtained include (a) the clade (Nuculanida+Opponobranchia), i.e., the traditionally defined Protobranchia; (b) the monophyly of Pteriomorphia; (c) the clade (Archiheterodonta+Palaeoheterodonta); (d) the monophyly of the traditionally defined Euheterodonta (including Anomalodesmata); and (e) the monophyly of Heteroconchia, i.e., (Palaeoheterodonta+Archiheterodonta+Euheterodonta). The stability of the basal tree topology to dataset manipulation is indicative of signal robustness in these four genes. The inferred tree topology corresponds closely to those obtained by datasets dominated by nuclear ribosomal genes (18S rRNA and 28S rRNA), controverting recent taxonomic actions based solely upon mitochondrial gene phylogenies.

The influence of rate heterogeneity among sites on the time dependence of molecular rates

Molecular evolutionary rate estimates have been shown to depend on the time period over which they are estimated. Factors such as demographic processes, calibration errors, purifying selection, and the heterogeneity of substitution rates among sites (RHAS) are known to affect the accuracy with which rates of evolution are estimated. We use mathematical modeling and Bayesian analyses of simulated sequence alignments to explore how mutational hotspots can lead to time-dependent rate estimates. Mathematical modeling shows that underestimation of molecular rates over increasing time scales is inevitable when RHAS is ignored. Although a gamma distribution is commonly used to model RHAS, we show that when the actual RHAS deviates from a gamma-like distribution, rates can either be under- or overestimated in a time-dependent manner. Simulations performed under different scenarios of RHAS confirm the mathematical modeling and demonstrate the impacts of time-dependent rates on estimates of divergence times. Most notably, erroneous rate estimates can have narrow credibility intervals, leading to false confidence in biased estimates of rates, and node ages. Surprisingly, large errors in estimates of overall molecular rate do not necessarily generate large errors in divergence time estimates. Finally, we illustrate the correlation between time-dependent rate patterns and differential saturation between quickly and slowly evolving sites. Our results suggest that data partitioning or simple nonparametric mixture models of RHAS significantly improve the accuracy with which node ages and substitution rates can be estimated.

Postglacial recolonization of eastern Blacknose Dace,Rhinichthys atratulus(Teleostei: Cyprinidae), through the gateway of New England

During the last ice age, much of North America far south as 40°N was covered by glaciers (Hewitt 2000). About 20,000 years ago, as the glaciers retreated, the hydrologic landscape changed dramatically creating waterways for fish dispersal. The number of populations responsible for recolonization and the regions from which they recolonized are unknown for many freshwater fishes living in New England and southeastern Canada. The Blacknose Dace,Rhinichthys atratulus, is one of the freshwater fish species that recolonized this region. We hypothesize that the earliest deglaciated region, modern-day Connecticut, was recolonized byR. atratulusvia a single founding event by a single population. In this paper, we test this hypothesis phylogenetically with regard to the major drainage basins within Connecticut. The mitochondrial DNA exhibits low nucleotide diversity, high haplotype diversity, and a dominant haplotype found across the state. A small percentage of individuals in the Housatonic drainage basin, however, share a haplotype with populations in New York drainage basins, a haplotype not found elsewhere in Connecticut's drainage basins. We calculated a range for the rate of divergence for NADH dehydrogenase subunit 2 (nd2) and control region (ctr) of 4.43-6.76% and 3.84-8.48% per million years (my), respectively. While this range is higher than the commonly accepted rate of 2% for mitochondrial DNA, these results join a growing list of publications finding high rates of divergence for various taxa (Peterson and Masel 2009). The data support the conclusion that Connecticut as a whole was recolonized initially by a single founding event that came from a single refugium. Subsequently, the Housatonic basin alone experienced a secondary recolonization event.

Evaluating the impact of post-mortem damage in ancient DNA: a theoretical approach

The growth of ancient DNA research has offered exceptional opportunities and raised great expectations, but has also presented some considerable challenges. One of the ongoing issues is the impact of post-mortem damage in DNA molecules. Nucleotide alterations and DNA strand breakages lead to a significant decrease in the quantity of DNA molecules of useful length in a sample and to errors in the final DNA sequences obtained. We present a model of age-dependent DNA damage and quantify the influence of that damage on subsequent steps in the sequencing process, including the polymerase chain reaction and cloning. Calculations using our model show that deposition conditions, rather than the age of a sample, have the greatest influence on the level of DNA damage. In turn, this affects the probability of interpreting an erroneous (possessing damage-derived mutations) sequence as being authentic. We also evaluated the effect of post-mortem damage on real data sets using a Bayesian phylogenetic approach. According to our study, damage-derived sequence alterations appear to have little impact on the final DNA sequences. This indicates the effectiveness of current methods for sequence authentication and validation.