A high-coverage genome sequence from an archaic Denisovan individual - PubMed (original) (raw)

. 2012 Oct 12;338(6104):222-6.

doi: 10.1126/science.1224344. Epub 2012 Aug 30.

Martin Kircher, Marie-Theres Gansauge, Heng Li, Fernando Racimo, Swapan Mallick, Joshua G Schraiber, Flora Jay, Kay Prüfer, Cesare de Filippo, Peter H Sudmant, Can Alkan, Qiaomei Fu, Ron Do, Nadin Rohland, Arti Tandon, Michael Siebauer, Richard E Green, Katarzyna Bryc, Adrian W Briggs, Udo Stenzel, Jesse Dabney, Jay Shendure, Jacob Kitzman, Michael F Hammer, Michael V Shunkov, Anatoli P Derevianko, Nick Patterson, Aida M Andrés, Evan E Eichler, Montgomery Slatkin, David Reich, Janet Kelso, Svante Pääbo

Affiliations

A high-coverage genome sequence from an archaic Denisovan individual

Matthias Meyer et al. Science. 2012.

Abstract

We present a DNA library preparation method that has allowed us to reconstruct a high-coverage (30×) genome sequence of a Denisovan, an extinct relative of Neandertals. The quality of this genome allows a direct estimation of Denisovan heterozygosity indicating that genetic diversity in these archaic hominins was extremely low. It also allows tentative dating of the specimen on the basis of "missing evolution" in its genome, detailed measurements of Denisovan and Neandertal admixture into present-day human populations, and the generation of a near-complete catalog of genetic changes that swept to high frequency in modern humans since their divergence from Denisovans.

PubMed Disclaimer

Figures

Figure 1

Figure 1

For single-stranded library preparation, ancient DNA molecules are dephosphorylated and heat-denatured. A biotinylated adaptor oligonucleotide is ligated to 3′-ends of the molecules, which are immobilized on streptavidin-coated beads and copied by extension of a primer hybridized to the adaptor. One strand of a double-stranded adaptor is then ligated to the newly synthesized strand. Finally, the beads are destroyed by heat to release the library molecules (not shown).

Figure 2

Figure 2

Average sequence divergence and branch length differences between the Denisovan genome and 11 present-day humans represented as a tree. Divergence is reported as fraction of the branch leading from human to the common ancestor with chimpanzee, and in years, assuming a human-chimpanzee divergence of 6.5 million years ago.

Figure 3

Figure 3

Maximum likelihood tree relating the Denisovan genome and the genomes of eleven present-day humans, allowing one migration event (shown in yellow).

Figure 4

Figure 4

(A) Sharing of derived alleles among present-day humans, Denisovans and Neandertals. We compare all possible pairs of 11 present-day humans {H1, H2} in their “_D-_statistics”, which measure the rate at which they share derived alleles with Denisovans (x-axis) and Neandertals (y-axis). Each point reports ±1 standard error bars from a Block Jackknife. _D_-statistics are color-coded by geographic region. The _D-_statistic is not the same as the mixture proportion; it is also affected by quantities like the amount of shared genetic drift between the samples being compared. (B) Sharing of derived alleles that are absent in Africans among present-day humans, Denisovans and Neandertals. We enhance the power of the _D-_statistics by restricting to sites where 35 sub-Saharan African samples have the ancestral allele, and pooling modern humans by region to increase resolution (bars again give one standard error). Eastern non-African populations have significantly more archaic ancestry than European populations (Z=5.3 and Z=4.8 for the tests based on the Denisovan and Neandertal _D-_statistics, respectively).

Figure 5

Figure 5

(A) Heterozygosity shown by the distribution of the number of bases matching the human reference genome at sites sampled to 20-fold coverage. The Y-axis is scaled to show the peak representing heterozygous sites in the center. (B) Inference of population size change over time using variation in the time since the most recent common ancestors across the genome shows that Denisovans have had a small population size over the last few hundred thousand years compared with modern humans, but a similar demographic history earlier. The y-axis specifies a number proportional to the population size Ne. The x-axis specifies time in units of divergence per base pair (along the top in years, assuming rates of 0.5×10−9 to 1.0×10−9 per year). Thin red lines around the Denisovan curve represent 100 bootstraps, thus showing the uncertainty of the inference. (C) The small population size in Denisovans is reflected in a greater accumulation of non-synonymous sites (normalized by the number of synonymous sites), whether measured in terms of heterozygous sites in Denisovans vs. modern humans (ratio 2.0 –2.5), or the accumulation of divergent sites on the Denisovan lineage divided by modern human lineages (ratio 1.5–2.0). The analysis is restricted to non-synonymous sites predicted to have a possibly or probably damaging effect on protein structure or function.

Similar articles

Cited by

References

    1. Green RE, et al. A draft sequence of the Neandertal genome. Science. 2010;328:710. - PMC - PubMed
    1. Reich D, et al. Genetic history of an archaic hominin group from Denisova Cave in Siberia. Nature. 2010;468:1053. - PMC - PubMed
    1. Hublin JJ. Out of Africa: modern human origins special feature: the origin of Neandertals. Proc Natl Acad Sci USA. 2009;106:16022. - PMC - PubMed
    1. Krause J, et al. The complete mitochondrial DNA genome of an unknown hominin from southern Siberia. Nature. 2010;464:894. - PMC - PubMed
    1. Gibbons A. Paleoanthropology. Who were the Denisovans? Science. 2011;333:1084. - PubMed

Publication types

MeSH terms

Grants and funding

LinkOut - more resources