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

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

Mycobacterium leprae diversity and population dynamics in medieval Europe from novel ancient genomes

Background: Hansen’s disease (leprosy), widespread in medieval Europe, is today mainly prevalent in tropical and subtropical regions with around 200,000 new cases reported annually. Despite its long history and appearance in historical records, its origins and past dissemination patterns are still widely unknown. Applying ancient DNA approaches to its major causative agent, Mycobacterium leprae, can significantly improve our understanding of the disease’s complex history. Previous studies have identified a high genetic continuity of the pathogen over the last 1500 years and the existence of at least four M. leprae lineages in some parts of Europe since the Early Medieval period. Results: Here, we reconstructed 19 ancient M. leprae genomes to further investigate M. leprae’s genetic variation in Europe, with a dedicated focus on bacterial genomes from previously unstudied regions (Belarus, Iberia, Russia, Scotland), from multiple sites in a single region (Cambridgeshire, England), and from two Iberian leprosaria. Overall, our data confirm the existence of similar phylogeographic patterns across Europe, including high diversity in leprosaria. Further, we identified a new genotype in Belarus. By doubling the number of complete ancient M. leprae genomes, our results improve our knowledge of the past phylogeography of M. leprae and reveal a particularly high M. leprae diversity in European medieval leprosaria. Conclusions: Our findings allow us to detect similar patterns of strain diversity across Europe with branch 3 as the most common branch and the leprosaria as centers for high diversity. The higher resolution of our phylogeny tree also refined our understanding of the interspecies transfer between red squirrels and humans pointing to a late antique/early medieval transmission. Furthermore, with our new estimates on the past population diversity of M. leprae, we gained first insights into the disease’s global history in relation to major historic events such as the Roman expansion or the beginning of the regular transatlantic long distance trade. In summary, our findings highlight how studying ancient M. leprae genomes worldwide improves our understanding of leprosy’s global history and can contribute to current models of M. leprae’s worldwide dissemination, including interspecies transmissions

Mycobacterium leprae diversity and population dynamics in medieval Europe from novel ancient genomes.

Funder: Max-Planck SocietyFunder: St John’s College, CambridgeFunder: Fondation Raoul FollereauFunder: University of Zurich’s University Research Priority Program “Evolution in Action: From Genomes to Ecosystems”Funder: the Senckenberg Centre for Human Evolution and Palaeoenvironment (S-HEP) at the University of TübingenBackgroundHansen's disease (leprosy), widespread in medieval Europe, is today mainly prevalent in tropical and subtropical regions with around 200,000 new cases reported annually. Despite its long history and appearance in historical records, its origins and past dissemination patterns are still widely unknown. Applying ancient DNA approaches to its major causative agent, Mycobacterium leprae, can significantly improve our understanding of the disease's complex history. Previous studies have identified a high genetic continuity of the pathogen over the last 1500 years and the existence of at least four M. leprae lineages in some parts of Europe since the Early Medieval period.ResultsHere, we reconstructed 19 ancient M. leprae genomes to further investigate M. leprae's genetic variation in Europe, with a dedicated focus on bacterial genomes from previously unstudied regions (Belarus, Iberia, Russia, Scotland), from multiple sites in a single region (Cambridgeshire, England), and from two Iberian leprosaria. Overall, our data confirm the existence of similar phylogeographic patterns across Europe, including high diversity in leprosaria. Further, we identified a new genotype in Belarus. By doubling the number of complete ancient M. leprae genomes, our results improve our knowledge of the past phylogeography of M. leprae and reveal a particularly high M. leprae diversity in European medieval leprosaria.ConclusionsOur findings allow us to detect similar patterns of strain diversity across Europe with branch 3 as the most common branch and the leprosaria as centers for high diversity. The higher resolution of our phylogeny tree also refined our understanding of the interspecies transfer between red squirrels and humans pointing to a late antique/early medieval transmission. Furthermore, with our new estimates on the past population diversity of M. leprae, we gained first insights into the disease's global history in relation to major historic events such as the Roman expansion or the beginning of the regular transatlantic long distance trade. In summary, our findings highlight how studying ancient M. leprae genomes worldwide improves our understanding of leprosy's global history and can contribute to current models of M. leprae's worldwide dissemination, including interspecies transmissions

Population genomics of the Viking world.

The maritime expansion of Scandinavian populations during the Viking Age (about AD 750-1050) was a far-flung transformation in world history1,2. Here we sequenced the genomes of 442 humans from archaeological sites across Europe and Greenland (to a median depth of about 1×) to understand the global influence of this expansion. We find the Viking period involved gene flow into Scandinavia from the south and east. We observe genetic structure within Scandinavia, with diversity hotspots in the south and restricted gene flow within Scandinavia. We find evidence for a major influx of Danish ancestry into England; a Swedish influx into the Baltic; and Norwegian influx into Ireland, Iceland and Greenland. Additionally, we see substantial ancestry from elsewhere in Europe entering Scandinavia during the Viking Age. Our ancient DNA analysis also revealed that a Viking expedition included close family members. By comparing with modern populations, we find that pigmentation-associated loci have undergone strong population differentiation during the past millennium, and trace positively selected loci-including the lactase-persistence allele of LCT and alleles of ANKA that are associated with the immune response-in detail. We conclude that the Viking diaspora was characterized by substantial transregional engagement: distinct populations influenced the genomic makeup of different regions of Europe, and Scandinavia experienced increased contact with the rest of the continent

Evolutionary Analysis of Ancient DNA Sequences

Studies of genetic data from old specimens – ancient DNA – offer a unique possibility to explore evolutionary history. Ancient DNA allows examination of past populations, extinct prehistoric species, and even large-scale analyses of past ecosystems. Since the beginnings of ancient DNA research in the mid 1980s, the field has undergone major technological and methodological transitions. However, investigations of the theory and methods used to analyse ancient DNA sequences have not kept up with the rapid pace of data generation. This thesis investigates a number of issues associated with the evolutionary analysis of ancient DNA. First, I built a model of post-mortem DNA degradation and examined its potential effects on the resulting sequence data (Chapter 2). Second, I conducted a large-scale study of time-dependent rates in order to look into their causes, characteristics, and ubiquity (Chapter 3). Third, I explored the impact of incorporating uncertainties associated with sample age into phylogenetic analysis (Chapters 4 and 5). Fourth, I evaluated the performance of Bayesian skyline plots in detecting recent population bottlenecks (Chapter 6). Finally, I combined several of these developments in a biogeographical analysis of the extinct cave bear (Chapter 7). My thesis shows that many of the common practices in ancient DNA studies, such as employing sequence-authentication criteria or ignoring age uncertainty in phylogenetic analyses, appear to be effective and to produce reliable results. However, there are areas in which more caution is needed when interpreting the results of DNA analyses. These include estimations of evolutionary rates, which are highly sensitive to the calibration points that are used; and the poor performance of Bayesian skyline reconstructions of recent population-size changes. The results of my studies improve our understanding of ancient DNA research and will serve as a useful guide for future evolutionary analyses of ancient DNA sequences

Evolutionary Analysis of Ancient DNA Sequences

Studies of genetic data from old specimens – ancient DNA – offer a unique possibility to explore evolutionary history. Ancient DNA allows examination of past populations, extinct prehistoric species, and even large-scale analyses of past ecosystems. Since the beginnings of ancient DNA research in the mid 1980s, the field has undergone major technological and methodological transitions. However, investigations of the theory and methods used to analyse ancient DNA sequences have not kept up with the rapid pace of data generation. This thesis investigates a number of issues associated with the evolutionary analysis of ancient DNA. First, I built a model of post-mortem DNA degradation and examined its potential effects on the resulting sequence data (Chapter 2). Second, I conducted a large-scale study of time-dependent rates in order to look into their causes, characteristics, and ubiquity (Chapter 3). Third, I explored the impact of incorporating uncertainties associated with sample age into phylogenetic analysis (Chapters 4 and 5). Fourth, I evaluated the performance of Bayesian skyline plots in detecting recent population bottlenecks (Chapter 6). Finally, I combined several of these developments in a biogeographical analysis of the extinct cave bear (Chapter 7). My thesis shows that many of the common practices in ancient DNA studies, such as employing sequence-authentication criteria or ignoring age uncertainty in phylogenetic analyses, appear to be effective and to produce reliable results. However, there are areas in which more caution is needed when interpreting the results of DNA analyses. These include estimations of evolutionary rates, which are highly sensitive to the calibration points that are used; and the poor performance of Bayesian skyline reconstructions of recent population-size changes. The results of my studies improve our understanding of ancient DNA research and will serve as a useful guide for future evolutionary analyses of ancient DNA sequences

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

Phylogenetic estimation of timescales using ancient DNA:the effects of temporal sampling scheme and uncertainty in sample ages

In recent years, ancient DNA has increasingly been used for estimating molecular timescales, particularly in studies of substitution rates and demographic histories. Molecular clocks can be calibrated using temporal information from ancient DNA sequences. This information comes from the ages of the ancient samples, which can be estimated by radiocarbon dating the source material or by dating the layers in which the material was deposited. Both methods involve sources of uncertainty. The performance of Bayesian phylogenetic inference depends on the information content of the data set, which includes variation in the DNA sequences and the structure of the sample ages. Various sources of estimation error can reduce our ability to estimate rates and timescales accurately and precisely. We investigated the impact of sample-dating uncertainties on the estimation of evolutionary timescale parameters using the software BEAST. Our analyses involved 11 published data sets and focused on estimates of substitution rate and root age. We show that, provided that samples have been accurately dated and have a broad temporal span, it might be unnecessary to account for sample-dating uncertainty in Bayesian phylogenetic analyses of ancient DNA. We also investigated the sample size and temporal span of the ancient DNA sequences needed to estimate phylogenetic timescales reliably. Our results show that the range of sample ages plays a crucial role in determining the quality of the results but that accurate and precise phylogenetic estimates of timescales can be made even with only a few ancient sequences. These finding

Empirical calibrated radiocarbon sampler: a tool for incorporating radiocarbon‐date and calibration error into Bayesian phylogenetic analyses of ancient DNA

Studies of DNA from ancient samples provide a valuable opportunity to gain insight into past evolutionary and demographic processes. Bayesian phylogenetic methods can estimate evolutionary rates and timescales from ancient DNA sequences, with the ages of the samples acting as calibrations for the molecular clock. Sample ages are often estimated using radiocarbon dating, but the associated measurement error is rarely taken into account. In addition, the total uncertainty quantified by converting radiocarbon dates to calendar dates is typically ignored. Here, we present a tool for incorporating both of these sources of uncertainty into Bayesian phylogenetic analyses of ancient DNA. This empirical calibrated radiocarbon sampler (ECRS) integrates the age uncertainty for each ancient sequence over the calibrated probability density function estimated for its radiocarbon date and associated error. We use the ECRS to analyse three ancient DNA data sets. Accounting for radiocarbon-dating and calibration error appeared to have little impact on estimates of evolutionary rates and related parameters for these data sets. However, analyses of other data sets, particularly those with few or only very old radiocarbon dates, might be more sensitive to using artificially precise sample ages and should benefit from use of the ECRS