Abstract

DNA sequence determination

The entire internal portion of the phalanx sample was used for DNA extraction in our clean-room facility, where procedures to minimize contamination from present-day human DNA are rigorously implemented^24,25 (Supplementary Information section 1). The DNA was treated with two enzymes: uracil-DNA-glycosylase, which removes uracil residues from DNA to leave abasic sites²⁶, and endonuclease VIII, which cuts DNA at the 5′ and 3′ sides of abasic sites. Subsequent incubation with T4 polynucleotide kinase and T4 DNA polymerase was used to generate 5′-phosphorylated blunt ends that are amenable to adaptor ligation. Because the great majority of uracil residues occur close to the ends of ancient DNA molecules, this procedure leads to only a moderate reduction in average length of the molecules in the library, but a several-fold reduction in uracil-derived nucleotide misincorporation²⁷.

Two independent sequencing libraries (SL3003 and SL3004) were created from the DNA, using a modified Illumina protocol²⁸ where a polymerase chain reaction (PCR) is used to add a 7-nucleotide index (in this case 5′-GTCGACT-3′) to the library molecules. This index ensures that the libraries are not contaminated by other sequencing libraries when they are taken out of the clean room to be sequenced²⁹. The libraries were sequenced on the Illumina Genome Analyser IIx platform for 101 cycles from each end of the molecules and an additional 7 cycles for determination of the index until almost every unique sequence in the libraries had been seen multiple times, that is, almost every clone present in the libraries has been sequenced (Supplementary Information section 1). Bases were called using the machine-learning algorithm Ibis³⁰ and an overlap of at least 11 bases was required for paired-end reads to be fused to full-molecule-size DNA sequences that were further analysed. This results in a greatly reduced error rate²⁷, although it removes the few molecules that are above 191 nucleotides in length from analysis (~0.1% in SL3003 and ~0.2% in SL3004). Sequences were mapped using the program BWA³¹ to the human (hg18/NCBI 36) and the chimpanzee (panTro2/CGSC 2.1) genomes as well as to the inferred ancestral genome of these species (from the six-way Enredo-Pecan-Ortheus alignment)³². PCR duplicates were identified and used to further increase sequence accuracy by calling consensus sequences.

A total of 82,227,320 sequences mapped uniquely (mapping quality $30) to the human genome, yielding about 5.2 gigabases of DNA sequences (1.9-fold genomic coverage), and 72,304,848 sequences mapped uniquely to the chimpanzee genome. When the substitutions inferred to have occurred on the Denisova and the present-day human lineages were compared, the relative numbers of different classes of nucleotide substitutions are remarkably similar, and the excess number of candidate substitutions on the Denisova lineage relative to the present-day human lineage is only 1.7-fold (Supplementary Fig. 2.2 and Supplementary Table 2.4). This reflects an improvement in error rate over the Neanderthal genome by over an order of magnitude⁸ and is mainly due to the enzymatic removal of uracil residues from the Denisova DNA²⁷. We estimate that most errors in the Denisova DNA sequences are due to low genomic coverage rather than to any features typical of ancient DNA.

Human DNA contamination estimates

Although rigorous measures to prevent contamination of the experiments by DNA from present-day humans were implemented at all laboratory steps, it is impossible to completely prevent contamination because bone samples as well as reagents may be contaminated before they enter the clean-room facility. To estimate the levels of contamination in the sequences produced we used three approaches (Supplementary Information section 3).

First, we estimated the level of mtDNA contamination using 276 sequence positions where the Denisova mtDNA differs from >99% of present-day human mtDNAs. For library SL3003, we observed 7,433 unique sequences that covered such positions and 7,421 were of the Denisova type. For library SL3004 the corresponding numbers were 5,042 and 5,036, indicating that the mtDNA contamination in the libraries is on the order of 0.2% (95% confidence interval (CI): 0.1–0.3%) and 0.1% (CI: 0.1–0.3%), respectively.

Second, we identified sequences that are unique to the Y chromosome⁸. If the individual from whom the phalanx derives is female, the number of such sequences represents the extent of male DNA contamination. We found zero and three such Y chromosomal sequences in the two libraries, respectively, whereas 1,449 and 696 are expected if the individual is male. Thus, the bone derives from a female and male DNA contamination in the two libraries is on the order of 0.00% (CI: 0.00–0.25%) and 0.43% (CI: 0.09–1.26%), respectively.

Third, to estimate the extent of nuclear DNA contamination we used one library to identify positions where the Denisova individual carries an ancestral, that is, chimpanzee-like, sequence variant that among present-day humans is derived and not known to vary. We then examined sequences that map at these positions in the other library and determined if they carry the ancestral sequence or the derived sequence. Observation of a derived sequence in the second library could be due to one of three possibilities: that the DNA fragment in question comes from present-day human contamination; that the Denisova individual is heterozygous at the position in question; or that there has been a sequencing error. We implemented a maximum likelihood method that uses the number of independent observations of ancestral and derived states across positions to co-estimate contamination along with heterozygosity and sequencing error as nuisance parameters (Supplementary Information section 3). From this analysis, both libraries are inferred to have contamination rates of less than 1%.

Ancestral features and duplications

The Denisova draft genome sequence allows features that are ancestral in the Denisova genome and derived in present-day humans to be identified. We previously described a set of 10.5 million single nucleotide differences and about half a million insertion/deletions (indels) inferred to be due to changes that occurred on the human lineage since the split from the common ancestor with the chimpanzee⁸. Of these, 4,267,431 (40.5%) single nucleotide differences and 105,372 (22.0%) indels are covered by the Denisova sequences. We identified 129 inferred amino substitutions and 14 indels in the coding sequences of genes where the Denisova individual carries the ancestral alleles at positions where present-day humans carry derived alleles and are not known to vary (Supplementary Information section 4). We also identified 90 such sites in 5′ untranslated regions (UTRs), 392 in 39 UTRs, two in microRNA genes and 104 in human accelerated regions. When we compared the Denisova and Neanderthal genomes we found that they carry the same assigned state at single nucleotide differences in 87.9% of the ancestral positions and 97.7% of the derived positions. The results for indels are similar: 87.6% for ancestral states and 98.6% for the derived states (Supplementary Table 4.3).

We analysed the segmental duplication content of the Denisova genome by detecting regions with an excess read depth (Supplementary Information section 5). In a three-way comparison of Denisova, Neanderthal and present-day human genomes, we found an excess of private Denisova duplications (2.27 megabases (Mb)) compared with duplications that were private in Neanderthals (0.60 Mb) or present-day humans (1.32 Mb). These regions were identified based on signatures of both excess read depth and increased sequence divergence, making them unlikely to be artefacts. We also identified two regions where the duplication architecture of Denisova is more similar to that of chimpanzee than to that of either Neanderthals or present-day humans, including two chromosomal regions associated with neurological disease in humans: spinal muscular atrophy on 5q13 (including SMN2, one of the most recent gene duplications in the human lineage) and neuropsychiatric disease on 16p12.1.

Relationship to Neanderthals and modern humans

A fundamental question is whether the Denisova individual is an outgroup to Neanderthals and modern humans, as the mtDNA suggests¹⁹, whether it is a sister group to Neanderthals or to modern humans, or whether it falls within the range of variation of either of these two groups. We addressed this by estimating the divergence between the Denisova and the human genome reference sequence as a fraction of the divergence between present-day humans and the common ancestor shared with the chimpanzee. To do this, we scored the frequency with which the Denisova genome carries the human versus the chimpanzee state at positions where the human and chimpanzee reference genomes differ; assuming constant evolutionary rates (Supplementary Information section 2). We restricted this analysis to the parts of the human reference genome that are of African ancestry³³ as gene flow from Neanderthals to non-Africans⁸ could otherwise complicate these analyses. The Denisova genome diverged from the reference human genome 11.7% (CI: 11.4–12.0%) of the way back along the lineage to the human–chimpanzee ancestor. For the Vindija Neanderthal, the divergence is 12.2% (CI: 11.9–12.5%). Thus, whereas the divergence of the Denisova mtDNA to present-day human mtDNAs is about twice as deep as that of Neanderthal mtDNA¹⁹, the average divergence of the Denisova nuclear genome from present-day humans is similar to that of Neanderthals.

A possible explanation for the similar divergence of the Denisova individual and Neanderthals from present-day Africans is that they both descend from a common ancestral population that separated earlier from ancestors of present-day humans. Such a scenario would predict a closer relationship between the Denisova individual and Neanderthals than between either of them and present-day humans. To test this prediction, we estimated the divergence between pairs of seven ancient and modern genomes (Denisova, Neanderthals, French, Han, Papuan, Yoruba and San), using an approach where we correct for error rates in each genome based on the assumption that each has the same number of true differences from chimpanzee (Supplementary Information section 6). The average divergence between Denisova and Vindija Neanderthals is estimated to be 9.84% of the way to the chimpanzee–human ancestor; that is, less than the average 12.38% divergence of both from present-day Africans. Assuming 6.5 million years for human–chimpanzee divergence, this implies that DNA sequences of Neanderthals and the Denisova individual diverged on average 640,000 years ago, and from present-day Africans 804,000 years ago.

To analyse further the relationship of the Denisova individual and Neanderthals, we aligned Denisova, Neanderthal and Yoruba sequences to the chimpanzee genome, picked a single sequence at random to represent each group, and examined sites where two copies of a derived and one copy of an ancestral allele were observed. Sequencing errors are expected to make a negligible contribution at such sites. The number of sites where the Denisova individual and Neanderthal cluster to the exclusion of the Yoruba and chimpanzee is 46,362, compared with an average of 22,012 sites for the other two possible patterns (Yoruba and Denisova, or Yoruba and Neanderthal). This excess of sites where Denisova and Neanderthal cluster supports the view that the Denisova individual and Neanderthals share a common history since separating from the ancestors of modern humans (Supplementary Information section 6).

A Neanderthal-specific bottleneck

The fact that the Denisova nuclear genome on average shares a more recent common ancestor with Neanderthal than with present-day humans raises the question of whether the overall DNA sequence divergence of the Denisova individual falls inside the group morphologically and geographically defined as Neanderthals, or if it represents a sister group to Neanderthals.

To investigate this question, we took advantage of the fact that in addition to the three individuals from Vindija Cave, Croatia, from which most of the Neanderthal genome sequences were produced, we have determined nuclear DNA sequences from three further Neanderthal individuals from Russia, Spain and Germany⁸. Of these, the 60,000–70,000-year-old skeleton of a Neanderthal child found in Mezmaiskaya Cave, Russia, is both oldest and geographically closest to the Denisova individual. Using the 56 Mb of autosomal DNA sequences determined from this specimen⁸, we estimate that the DNA sequence divergence between the Vindija and Mezmaiskaya Neanderthals corresponds to a date of 140,000 ± 33,000 years ago (Supplementary Information section 6) (Fig. 1). This remarkably low divergence—which is about one-third of the closest pair of present-day humans that we analysed—is in agreement with the observation that diversity among Neanderthal mtDNAs is low relative to present-day humans¹⁰ and indicates that the Vindija and Mezmaiskaya Neanderthals descend from a common ancestral population that experienced a drastic bottleneck since separating from the ancestors of the Denisova individual.

Open in a separate window

Figure 1

A neighbour-joining tree based on pairwise autosomal DNA sequence divergences for five ancient and five present-day hominins

Vindija 33.16, Vindija 33.25 and Vindija 33.26 refer to the catalogue numbers of the Neanderthal bones.

To understand further the bottleneck in the history of Vindija and Mezmaiskaya Neanderthals, we examined four-way alignments of the Vindija Neanderthal genome sequence, the Mezmaiskaya Neanderthal, the Denisova individual and the chimpanzee genome. At transversion substitutions where two copies of the derived alleles are observed, we detect 924 substitutions that cluster the Vindija and Mezmaiskaya Neanderthals, 80 that cluster Vindija and Denisova, and 81 that cluster Mezmaiskaya and Denisova. This corresponds to at least a 65% probability that the DNA sequences in the Neanderthals share a common ancestor more recently than their split from the ancestor of the Denisova individual (Supplementary Information section 7). It is much higher than the 15–20% probability associated with the ‘Out of Africa’ bottleneck common to present-day non-Africans³⁴. If we replace the Mezmaiskaya Neanderthal in this analysis with a Neanderthal from El Sidron, Spain, or from Feldhofer, Germany, results are qualitatively similar although numbers are smaller (Supplementary Information section 7). Thus, we conclude that late Neanderthals across a broad geographical range have a population history distinct from that of the Denisova individual in that they share a strong population bottleneck not experienced by the ancestors of the Denisova individual. We call the group to which this individual belonged Denisovans in analogy to Neanderthals, as Denisovans are described for the first time based on molecular data from Denisova Cave just as Neanderthals were first described based on skeletal remains retrieved in the Neander Valley in Germany.

No Denisovan gene flow into all Eurasians

We have previously shown that Vindija Neanderthals share more derived alleles with non-Africans than with Africans, consistent with Neanderthals contributing 1–4% of the genomes of present-day humans across Eurasia⁸. To investigate the extent to which the Denisova individual shares this pattern, we examined alignments of sets of four genomes, each consisting of an African (Yoruba or San), a Eurasian (French or Han), an archaic hominin (Neanderthal or Denisovan) and the chimpanzee. We randomly sampled one allele from each of the three hominins, and counted all transversion differences between the African and the Eurasian where the archaic individual carries the derived allele (the ‘D statistics’ of ref. 8). Neanderthals match the French genome on average 4.6 ± 0.7% more often than they match the Yoruba genome (Table 1). Although the Denisova individual also matches the French more than the Yoruba genome, this skew is significantly less strong at 1.8 ± 0.5%. The estimates of D statistics were quantitatively consistent (within two standard deviations) for all other choices of Eurasian and African populations (Table 1). These findings indicate that the archaic component of the Eurasian gene pool is less closely related to the Denisova individual than to Neanderthals.

Table 1

Sharing of derived alleles between present-day and archaic hominins

Sample H1	Sample H2	Source of data for H1 and H2	D(H1, H2, Neanderthal, chimpanzee)					D(H1, H2, Denisova, chimpanzee)
			n BABA	n ABBA	D (%)	s.e. (%)	Z-score	n BABA	n ABBA	D (%)	s.e. (%)	Z-score
Eurasian/Eurasian*
French	Han	Ref. 8	17,214	17,602	−1.1	0.8	−1.4	27,250	27,265	0.0	0.6	0.0
Karitiana	Sardinian	This study	1,116	1,085	1.4	2.1	0.7	1,559	1,627	−2.1	1.8	−1.2
Karitiana	Cambodian	This study	1,683	1,707	−0.7	1.8	−0.4	2,371	2,460	−1.8	1.5	−1.2
Karitiana	Mongolian	This study	1,128	1,195	−2.9	2.2	−1.3	1,765	1,742	0.7	1.8	0.4
Sardinian	Cambodian	This study	2,592	2,670	−1.5	1.5	−1.0	3,935	3,925	0.1	1.2	0.1
Sardinian	Mongolian	This study	1,966	2,027	−1.5	1.6	−0.9	3,036	3,057	−0.3	1.3	−0.3
Cambodian	Mongolian	This study	2,811	2,804	0.1	1.4	0.1	4,442	4,342	1.1	1.2	1.0
African/African*
San	Yoruba	Ref. 8	23,690	23,855	−0.3	0.6	−0.6	39,042	39,019	0.0	0.5	0.1
Melanesian/Melanesian*
Papuan2	Bougainville	This study	3,351	3,284	1.0	1.3	0.8	5,319	5,140	1.7	1.1	1.5
Eurasian/African*
French	San	Ref. 8	25,242	22,982	4.7	0.6	7.6†	39,838	38,495	1.7	0.5	3.4†
French	Yoruba	Ref. 8	21,794	19,890	4.6	0.7	6.9†	34,262	33,078	1.8	0.5	3.6†
Han	San	Ref. 8	25,081	22,470	5.5	0.6	8.5†	38,815	37,439	1.8	0.5	3.4†
Han	Yoruba	Ref. 8	21,741	19,412	5.7	0.7	7.9†	33,182	32,184	1.5	0.5	2.8
Karitiana	Mbuti	This study	1,577	1,473	3.4	1.9	1.8	2,368	2,360	0.2	1.5	0.1
Sardinian	Mbuti	This study	2,562	2,400	3.3	1.5	2.2	4,028	3,784	3.1	1.2	2.6
Cambodian	Mbuti	This study	4,235	3,641	7.5	1.2	6.5†	6,329	5,850	3.9	1.0	4.0†
Mongolian	Mbuti	This study	3,077	2,765	5.3	1.4	3.9†	4,514	4,505	0.1	1.1	0.1
Eurasian/Melanesian*
French	Papuan1	Ref. 8	15,523	15,548	−0.1	0.8	−0.1	23,509	25,470	−4.0	0.7	−5.7†
Han	Papuan1	Ref. 8	15,059	14,677	1.3	0.9	1.5	22,262	24,198	−4.2	0.7	−5.8†
Karitiana	Papuan2	This study	1,522	1,658	−4.3	1.9	−2.2	2,201	2,641	−9.1	1.6	−5.8†
Karitiana	Bougainville	This study	1,577	1,717	−4.3	1.8	−2.4	2,229	2,671	−9.0	1.5	−5.9†
Sardinian	Papuan2	This study	2,447	2,647	−3.9	1.5	−2.6	3,714	4,150	−5.5	1.2	−4.5†
Sardinian	Bougainville	This study	2,531	2,762	−4.4	1.5	−3.0	3,877	4,336	−5.6	1.1	−4.9†
Cambodian	Papuan2	This study	3,713	3,891	−2.3	1.3	−1.8	5,457	6,272	−6.9	1.1	−6.5†
Cambodian	Bougainville	This study	3,847	3,994	−1.9	1.2	−1.6	5,751	6,333	−4.8	1.0	−4.7†
Mongolian	Papuan2	This study	2,783	2,852	−1.2	1.5	−0.8	4,192	4,758	−6.3	1.2	−5.3†
Mongolian	Bougainville	This study	2,813	3,066	−4.3	1.5	−2.9	4,234	4,847	−6.8	1.1	−6.0†
Melanesian/African*
Papuan1	San	Ref. 8	21,985	20,366	3.8	0.7	5.1†	35,923	32,841	4.5	0.6	7.2†
Papuan1	Yoruba	Ref. 8	19,107	17,646	4.0	0.8	4.9†	30,995	28,186	4.7	0.6	7.4†
Papuan2	Mbuti	This study	3,832	3,324	7.1	1.3	5.4†	6,124	5,233	7.8	1.1	7.2†
Bougainville	Mbuti	This study	4,216	3,596	7.9	1.2	6.8†	6,498	5,633	7.1	1.1	6.7†

Open in a separate window

We present the D statistic D(H₁, H₂, X, chimpanzee), the normalized difference between the number of sites at which the derived allele in an archaic read from X matches human sample H₁ (n_BABA) and human sample H₂ (n_BABA); thus, its value is D = (n_BABA − n_ABBA)/(n_BABA + n_ABBA). We restrict to autosomal transversion substitutions, compute standard errors (s.e.) from a block jackknife, and highlight (dagger symbol) the D statistics that are more than Z> 3 s.d. from zero. Both Neanderthals and Denisovans match Eurasians more than the Africans, but the signals are consistently and significantly stronger when X = Neanderthal than when X = Denisova. The slight numerical differences with Table 4 of ref. 8 are due to differences in the data filtering. Here we restrict to comparisons of present-day human samples that were sequenced by the same protocol (the five individuals sequenced in ref. 8, or the seven in this study); Supplementary Table 8.2 presents the complete set of pairwise comparisons.

^*Comparison.

^†D statistics that are more than Z> 3 s.d. from zero.

We also examined 13 genomic regions that we previously identified as candidates for a contribution of archaic genetic material into non-Africans, based on their deeper genetic divergences in non-Africans than in Africans⁸. Using ‘tag SNPs’ that are informative about whether a haplotype is from the lineage unique to non-Africans, we find that the Denisova individual matches the deeply diverged non-African haplotype in 6 cases, whereas Neanderthals do so in 11 cases (Supplementary Information section 7). Thus, both Neanderthals and Denisovans are more related than would be expected by chance to these genomic segments, but the signal in Denisovans is weaker.

These analyses indicate that Neanderthals are more closely related than Denisovans to the population that contributed to the gene pool of the ancestors of present-day Eurasians. The fact that Eurasians share some additional affinity with the Denisova individual relative to Africans is compatible with a scenario in which Denisovans shared some of their history with Neanderthals before the gene flow from Neanderthals into modern humans occurred.

Denisovan gene flow into the ancestors of Melanesians

Although the Denisova individual derives from a population that was not directly involved in the gene flow from Neanderthals to Eurasians, it is possible that Denisovans admixed with the ancestors of present-day people in some parts of the Old World. To investigate this, we analysed the relationship of the Denisova genome to the genomes of 938 present-day humans from 53 populations that have been geno-typed at 642,690 single nucleotide polymorphisms (SNPs)³⁵. We scored each of these present-day humans based on their relative proximity to Neanderthals and the Denisova individual at positions where we have high-quality data for both the Neanderthal and Denisova genomes (Supplementary Information section 8). Using the means of the 53 populations, the first two principal components separate the populations into three groups (Fig. 2): first, the 7 sub-Saharan African populations; second, a group of 44 non-African populations as well as one north African group; and third, Papuan and Bougainville populations from Melanesia. When individuals from selected populations are analysed separately, the Papuan and Bougainville islanders remain distinct from almost all individuals outside Africa (Supplementary Fig. 8.1b). Thus, with respect to their relationship to Neanderthals and Denisovans, the Melanesian populations stand out relative to other non-African populations.

Open in a separate window

Figure 2

Relationship of present-day populations to the Denisova individual and Neanderthals based on 255,077 SNPs

Principal component analysis of the means of 53 present-day human populations projected onto the top two principal components defined by Denisova, Neanderthal and chimpanzee. The seven ‘African’ populations are San, Mbuti, Biaka, Bantu Kenya, Bantu South Africa, Yoruba and Mandenka; the ‘Non-African’ populations are 44 diverse groups from outside Africa except for Papuan and Bougainville islanders.

To explore this further, we analysed the relationship of the Denisova genome to the genomes of five present-day humans that we previously sequenced to about fivefold coverage⁸ (a Yoruba and a San genome from Africa, a French genome from Europe, a Han genome from China and a Papuan genome from Melanesia), as well as seven present-day humans that we sequenced to 1–2-fold coverage for this study (a Mbuti genome from Africa, a Sardinian genome from Europe, a Mongolian genome from Central Asia, a Cambodian genome from South-East Asia, an additional Papuan genome from Melanesia, a Bougainville islander genome from Melanesia, and a Karitiana genome from South America) (Supplementary Information section 9). We used the D statistic⁸ to test if various pairs of present-day humans share equal numbers of derived alleles with the Denisova individual. To do this, we restricted comparisons to pairs of present-day humans sequenced at the same time to minimize the chance that differences in sample processing could affect the results. We find that the fivefold coverage Papuan individual shares 4.0 ± 0.7% more alleles with the Denisova individual than does the French individual, and we observed a similar skew in all 10 comparisons of Melanesian and other non-African populations (Table 1). When we stratified the data by base substitution class and chromosome, the D statistics are qualitatively unchanged (Supplementary Information section 10). Similarly, the D statistics are consistent for all depths of read coverage, indicating that mapping errors, for example due to segmental duplications, are not likely to explain these results. Finally, differences in sequencing error rate across samples cannot explain the observed D statistics (Supplementary Information section 10).

Under the assumption that gene flow explains these observations, we determined the direction of this gene flow by asking whether Melanesians and other Eurasians share derived alleles with Africans equally often. If the gene flow was entirely into the ancestors of the Denisovan individual, we would not expect this to affect the relationship of Africans to Melanesians and other Eurasians and thus we would expect them to share derived alleles equally often with Africans. However, we find that derived alleles in Africans match Melanesians 3.4 ± 0.4% less often than other non-Africans (Z = 10.8). Because this skew is seen without using Denisovan data it cannot be explained by gene flow into Denisovans or, for example, by contamination of the Denisova sample by present-day Melanesian DNA. Thus, at least some of the putative gene flow must have been into Melanesians (Supplementary Information section 8).

When we compare the skew in the fraction of derived alleles shared with the two archaic hominins to what would be expected for individuals of 100% Neanderthal or Denisova ancestry, respectively (Supplementary Information section 8 and ref. 8), we estimate that 2.5 ± 0.6% of the genomes of non-African populations derive from Neanderthals, in agreement with our previous estimate of 1–4%⁸. In addition, we estimate that 4.8 ± 0.5% of the genomes of Melanesians derive from Denisovans. Altogether, as much as 7.4 ± 0.8% of the genomes of Melanesians may thus derive from recent admixture with archaic hominins.

A model of population history

To understand the implications of the relationships observed among the Denisova individual, the Neanderthals and present-day humans, we fit the D statistics described in the previous sections to a parameterized model of population history. The D statistics for the Denisova individual differ in two important ways from those for the Neanderthal. First, the Denisova individual shares fewer derived alleles with either the French or Han Chinese populations than do the Neanderthals. Second, the Denisova individual shares more derived alleles with the Papuans than do the Neanderthals. We are able to fit the data with a model that assumes the Denisovans are a sister group of Neanderthals with a population divergence time of one-half to two-thirds of the time to the common ancestor of Neanderthals and humans. After the divergence of the Denisovans from Neanderthals, there was gene flow from Neanderthals into the ancestors of all present-day non-Africans. Later there was admixture between the Denisovans and the ancestors of Melanesians that did not affect other non-African populations. This model is illustrated in Fig. 3 and is described in detail in Supplementary Information section 11.

Open in a separate window

Figure 3

A model of population history compatible with the data

N denotes effective population size, t denotes time of population separation, f denotes amount of gene flow and t_GF denotes time of gene flow.

Other, more complex models could also explain the data. For example, a model that invokes only gene flow from Denisovans to Melanesian ancestors outside Africa and assumes four subpopulations in Africa that existed between the times of the origin of Denisovan and Neanderthal ancestors and the ancestors of present-day Eurasians could also fit the data (Supplementary Fig. 11.4). However, because barriers to gene flow between such subpopulations would have to persist for hundreds of thousands of years to create the observed patterns, such a model is less plausible on biological grounds than a model that invokes two instances of gene flow outside Africa.

Discordance of mtDNA and nuclear histories

The population history indicated by the nuclear genome is different from that indicated by the mtDNA phylogeny. There are two possible explanations for this. One is that the mtDNA lineage was introduced into Denisovan ancestors by admixture from another hominin lineage for which we have no data. The other is that the discordance is the result of ‘incomplete lineage sorting’, that is, the random assortment of genetic lineages due to genetic drift which may have allowed a divergent mtDNA lineage to survive in Denisovans by chance while becoming lost in Neanderthals and modern humans. A large ancestral population size makes incomplete lineage sorting more likely to occur. In Supplementary Information section 11, we show that given reasonable assumptions about the size of the ancestral populations, the discordance of the mtDNA phylogeny with that indicated by the nuclear DNA can be explained either by a small amount of admixture from another archaic hominin or by incomplete lineage sorting. Thus, the data do not allow us to favour one hypothesis over the other.

A tooth from Denisova Cave

In 2000, a hominin tooth was discovered in layer 11.1 of the south gallery of Denisova Cave (Fig. 4a, b). The tooth is from a young adult and therefore from another individual than the phalanx which stems from a juvenile (Supplementary Information section 12). To elucidate the relationship of the tooth to the individual from which the phalanx is derived, we extracted DNA from 50 mg of dentin from the root of the tooth and prepared a sequencing library (Supplementary Information section 13). About 0.17% of random DNA sequences determined from this library aligned to the human genome, whereas the rest is likely to represent microbial contamination common in ancient bones. We therefore used a novel DNA capture approach³⁶ to isolate mtDNA sequences from the sequencing library. A total of 15,094 sequences were identified which allowed the complete mtDNA genome to be assembled at an average coverage of 58-fold. This sequence differs at two positions from the mtDNA of the phalanx whereas it differs at about 380 positions from both Neanderthal and present-day humans. The time since the most recent common ancestor of the two mtDNAs from Denisova Cave is estimated to be 7,500 years, with a 95% upper bound of 16,000 years (Supplementary Information section 13). We conclude that the tooth and the phalanx derive from two different individuals that are probably from the same hominin population.

Open in a separate window

Figure 4

Morphology of the Denisova molar

a, b, Occlusal (a) and mesial (b) views. c, Comparison of the Denisova molar to diverse third molars, in a biplot of the mesiodistal and buccolingual lengths (in mm). AMH, anatomically modern humans; SH, Sima de los Huesos. Supplementary Fig. 12.1 presents a similar comparison to second molars.

Morphology of the Denisova molar

The tooth is an almost complete left, probably third, but possibly second, upper molar (Fig. 4b). The crown is trapezoidal and tapers strongly distally, with bulging lingual and buccal walls giving the tooth an inflated appearance (Supplementary Information section 12). The roots are short but robust and strongly flaring.

Overall, the tooth is very large (mesiodistal diameter, 13.1 mm; buccolingual, 14.7 mm). As a third molar, it is outside the range of normal size variation of all fossil taxa of the genus Homo, with the exception of H. habilis and H. rudolfensis, and comparable to Australopithecines (Fig. 4c). Compared to second molars, it is larger than Neanderthals or early modern humans, but similar to H. erectus and H. habilis (Supplementary Fig. 12.1).

Besides size, it is also distinguished from most Neanderthal third molars by the absence of hypocone reduction, and from both second and third Neanderthal molars by the presence of a large talon basin and the strong flare of the crown. Furthermore, it lacks the lingual hypocone projection seen in all Neanderthal first and many second molars, and has strongly diverging roots, unlike the closely spaced and frequently fused roots of Neanderthals.

It is of particular interest to compare the Denisova molar to Middle Pleistocene hominins from China, where H. erectus and other archaic forms, sometimes interpreted as H. heidelbergensis, may have survived until recently. Unfortunately, very few of these fossils preserve third upper molars. Of the few examples that are available, most differ from the Denisova molar by their strongly reduced size. Second molars are more frequent than third molars, and most have a trapezoidal shape like Denisova, but they do not have the lingually skewed position of the hypocone and metacone and the strong basal flare of the crown.

The Denisova molar supports the DNA evidence that the Denisovan population is distinct from late Neanderthals as well as from modern humans. In fact, the primitive traits of the Denisova tooth suggest that Denisovans may have been separated from the Neanderthal lineage before Neanderthal dental features are documented in Western Eurasia (>300,000 years BP) (Supplementary Information section 12), although we cannot exclude the possibility that the Denisovan dental morphology results from a reversion.

Stratigraphy and dating

The small size of both the phalanx and the tooth precludes direct radiocarbon dating. We instead dated seven bone fragments found close to the hominin remains in layer 11 in the east and south galleries. To ensure that they were associated with human occupation of the cave we chose bones that have evidence of human modification, including a rib with regular incisions and a bone projectile point blank generally associated with Upper Palaeolithic cultural assemblages. In the south gallery, where modified bones were not available, we used herbivore bones (Supplementary Information section 12).

Four of the seven dates are infinite dates older than 50,000 years BP (uncalibrated), whereas three are finite dates between 16,000 and 30,000 years BP (Supplementary Table 12.1). The rib with incisions and the projectile point blank are about 30,000 and 23,000 years BP, respectively. Together with three previous dates²³ this shows that layer 11 contains cultural remains from at least two different time periods, one period older than 50,000 years BP and one more recent period. However, the stratigraphy is complicated by the discovery of a wedge-shaped area close to the area where the phalanx was found that is likely to be disturbed (Supplementary Information section 12). Hominin remains large enough to allow direct radiocarbon dates may eventually be discovered in the cave, but a reasonable hypothesis is that the phalanx and molar belong to the older occupation.

Discussion

The molecular preservation of the Denisova phalanx is exceptional in that the fraction of endogenous relative to microbial DNA is about 70%. By contrast, in all Neanderthal remains studied so far the relative abundance of endogenous DNA is below 5%, and typically below 1%. Furthermore, the average length of hominin DNA fragments in the Denisova phalanx is 58 base pairs (bp) (SL3003) and 74 bp (SL3004) in spite of the enzymatic treatment that removes uracil residues and decreases the average fragment size, whereas in most well-preserved Neanderthal samples it is 50 bp or smaller without this treatment. Thus, although many Neanderthals are preserved under conditions apparently similar to those in Denisova Cave, the Denisova phalanx is one of few bones found in temperate conditions that are as well preserved as many permafrost remains^37,38. It is not clear why this is. It is not due to some condition that affects all hominin remains in Denisova Cave because the fraction of endogenous DNA in the tooth is 0.17%; that is, typical of other Late Pleistocene hominin remains. It is possible that a rapid desiccation of the tissue after death, which would limit degradation of the DNA by endogenous enzymes as well as microbial growth, has allowed this exceptional preservation.

The Denisova individual and the population to which it belonged carry some exceptionally archaic molecular (mtDNA) as well as morphological (dental) features. Nevertheless, the picture that emerges from analysis of the nuclear genome is one where the Denisova population is a sister group to Neanderthals. Three possibilities could account for how such archaic features have come to be present in Denisovans. One possibility is that these features were retained in Denisovans but became lost in modern humans and Neanderthals. A second, not mutually exclusive, possibility is that they entered the Denisova population through gene flow from some even more diverged hominin. Although such gene flow cannot be detected with the current mtDNA and nuclear DNA data, further sequencing of other hominin remains may in the future allow testing for it. A third possibility that could account for the apparently archaic dental morphology, but not the mtDNA, is a reversal to ancestral traits.

After they diverged from one another, Denisovans and Neanderthals had largely separate population histories as shown by a number of observations. First, patterns of allele sharing indicate that Denisovan ancestors did not contribute genes at a detectable level to present-day people all over Eurasia whereas Neanderthals did⁸. Thus, Neanderthals at some point interacted with ancestors of present-day Eurasians independently of Denisovans. Second, the genetic diversity of Neanderthals across their geographical range in the last thirty or forty thousand years of their history was extremely low, indicating that they experienced one or more strong genetic bottlenecks independently of the Denisovans. Third, our results indicate that Denisovans but not Neanderthals contributed genes to ancestors of present-day Melanesians. Fourth, the dental morphology shows no evidence of any derived features seen in Neanderthals. In fact, dental remains from the Sima de los Huesos of Atapuerca, for which ages between 350,000 and 600,000 years have been proposed^39,40, already carry Neanderthal-like morphological features that are not seen in the Denisova molar.

An interesting question is how widespread Denisovans were. A possibility is that they lived in large parts of East Asia at the time when Neanderthals were present in Europe and western Asia. One observation compatible with this possibility is that Denisovan relatives seem to have contributed genes to present-day Melanesians but not to present-day populations which currently live much closer to the Altai region such as Han Chinese or Mongolians (Table 1). Thus, they have at least at some point been present in an area where they interacted with the ancestors of Melanesians and this was presumably not in southern Siberia. Further studies of both molecular and morphological features of hominin remains across Asia should clarify how widespread Denisovans were and how they were related to archaic hominins other than Neanderthals.

The Denisova individual belongs to a hominin group that shares a common ancestor with Neanderthals but has a distinct population history. We define this group based on genomic evidence and call it Denisovans, but refrain from any formal Linnaean taxonomic designations that would indicate species or subspecies status for either Neanderthals or Denisovans. In our view, these results show that on the Eurasian mainland there existed at least two forms of archaic hominins in the Upper Pleistocene: a western Eurasian form with morphological features that are commonly used to define them as Neanderthals, and an eastern form to which the Denisova individuals belong. In the future, when more complete genomes from these and other archaic hominins will be sequenced from remains that allow more morphological features to assessed, their relationships will become even better understood. This will be an important endeavour as the emerging picture of Upper Pleistocene hominin evolution is one in which gene flow among different hominin groups was common.

Supplementary Material

Acknowledgements

Footnotes

References