Dataset Information

MapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters.

ABSTRACT:

Motivation

Ancient DNA (aDNA) molecules in fossilized bones and teeth, coprolites, sediments, mummified specimens and museum collections represent fantastic sources of information for evolutionary biologists, revealing the agents of past epidemics and the dynamics of past populations. However, the analysis of aDNA generally faces two major issues. Firstly, sequences consist of a mixture of endogenous and various exogenous backgrounds, mostly microbial. Secondly, high nucleotide misincorporation rates can be observed as a result of severe post-mortem DNA damage. Such misincorporation patterns are instrumental to authenticate ancient sequences versus modern contaminants. We recently developed the user-friendly mapDamage package that identifies such patterns from next-generation sequencing (NGS) sequence datasets. The absence of formal statistical modeling of the DNA damage process, however, precluded rigorous quantitative comparisons across samples.

Results

Here, we describe mapDamage 2.0 that extends the original features of mapDamage by incorporating a statistical model of DNA damage. Assuming that damage events depend only on sequencing position and post-mortem deamination, our Bayesian statistical framework provides estimates of four key features of aDNA molecules: the average length of overhangs (?), nick frequency (?) and cytosine deamination rates in both double-stranded regions ( ) and overhangs ( ). Our model enables rescaling base quality scores according to their probability of being damaged. mapDamage 2.0 handles NGS datasets with ease and is compatible with a wide range of DNA library protocols.

Availability

mapDamage 2.0 is available at ginolhac.github.io/mapDamage/ as a Python package and documentation is maintained at the Centre for GeoGenetics Web site (geogenetics.ku.dk/publications/mapdamage2.0/).

Supplementary information

Supplementary data are available at Bioinformatics online.

SUBMITTER: Jonsson H

PROVIDER: S-EPMC3694634 | biostudies-literature | 2013 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters.

Jónsson Hákon H Ginolhac Aurélien A Schubert Mikkel M Johnson Philip L F PL Orlando Ludovic L

Bioinformatics (Oxford, England) 20130423 13

<h4>Motivation</h4>Ancient DNA (aDNA) molecules in fossilized bones and teeth, coprolites, sediments, mummified specimens and museum collections represent fantastic sources of information for evolutionary biologists, revealing the agents of past epidemics and the dynamics of past populations. However, the analysis of aDNA generally faces two major issues. Firstly, sequences consist of a mixture of endogenous and various exogenous backgrounds, mostly microbial. Secondly, high nucleotide misincorp ...[more]

PMID: 23613487

Similar Datasets

Project description:BackgroundThe Approximate Bayesian Computation (ABC) approach has been used to infer demographic parameters for numerous species, including humans. However, most applications of ABC still use limited amounts of data, from a small number of loci, compared to the large amount of genome-wide population-genetic data which have become available in the last few years.ResultsWe evaluated the performance of the ABC approach for three 'population divergence' models - similar to the 'isolation with migration' model - when the data consists of several hundred thousand SNPs typed for multiple individuals by simulating data from known demographic models. The ABC approach was used to infer demographic parameters of interest and we compared the inferred values to the true parameter values that was used to generate hypothetical "observed" data. For all three case models, the ABC approach inferred most demographic parameters quite well with narrow credible intervals, for example, population divergence times and past population sizes, but some parameters were more difficult to infer, such as population sizes at present and migration rates. We compared the ability of different summary statistics to infer demographic parameters, including haplotype and LD based statistics, and found that the accuracy of the parameter estimates can be improved by combining summary statistics that capture different parts of information in the data. Furthermore, our results suggest that poor choices of prior distributions can in some circumstances be detected using ABC. Finally, increasing the amount of data beyond some hundred loci will substantially improve the accuracy of many parameter estimates using ABC.ConclusionsWe conclude that the ABC approach can accommodate realistic genome-wide population genetic data, which may be difficult to analyze with full likelihood approaches, and that the ABC can provide accurate and precise inference of demographic parameters from these data, suggesting that the ABC approach will be a useful tool for analyzing large genome-wide datasets.

Project description:Efforts to mitigate global warming are often justified through calculations of the economic damages that may occur absent mitigation. The earliest such damage estimates were speculative mathematical representations, but some more recent studies provide empirical estimates of damages on economic growth that accumulate over time and result in larger damages than those estimated previously. These heightened damage estimates have been used to suggest that limiting global warming this century to 1.5 °C avoids tens of trillions of 2010 US$ in damage to gross world product relative to limiting global warming to 2.0 °C. However, in order to estimate the net effect on gross world product, mitigation costs associated with decarbonizing the world's energy systems must be subtracted from the benefits of avoided damages. Here, we follow previous work to parameterize the aforementioned heightened damage estimates into a schematic global climate-economy model (DICE) so that they can be weighed against mainstream estimates of mitigation costs in a unified framework. We investigate the net effect of mitigation on gross world product through finite time horizons under a spectrum of exogenously defined levels of mitigation stringency. We find that even under heightened damage estimates, the additional mitigation costs of limiting global warming to 1.5 °C (relative to 2.0 °C) are higher than the additional avoided damages this century under most parameter combinations considered. Specifically, using our central parameter values, limiting global warming to 1.5 °C results in a net loss of gross world product of roughly forty trillion US$ relative to 2 °C and achieving either 1.5 °C or 2.0 °C require a net sacrifice of gross world product, relative to a no-mitigation case, though 2100 with a 3%/year discount rate. However, the benefits of more stringent mitigation accumulate over time and our calculations indicate that stabilizing warming at 1.5 °C or 2.0 °C by 2100 would eventually confer net benefits of thousands of trillions of US$ in gross world product by 2300. The results emphasize the temporal asymmetry between the costs of mitigation and benefits of avoided damages from climate change and thus the long timeframe for which climate change mitigation investment pays off.