The preservation of endogenous DNA in fossils is correlated with the amount, composition, and chemical preservation of amino acids (
Poinar et al. 1996). We find that endogenous DNA can be amplified from Pleistocene remains when the amino acid content is more than 30,000 parts per million (ppm), the ratio of glycine to aspartic acid between two and ten, and the aspartic acid racemization (i.e., the stereoisomeric D/L ratio) less than 0.10 (
Poinar et al. 1996;
Krings et al. 1997, 2000;
Schmitz et al. 2002; data not shown). We analyzed the amino acid preservation of 24 Neandertal and 40 early modern human fossils (
Table S1). Several important Neandertal fossils, such as La Ferrassie and Krapina, as well as important modern human fossils, such as Veternica, proved to be too poorly preserved to be likely to allow DNA retrieval. Thus, further destructive sampling of these specimens was not considered justified. However, four Neandertal and five early modern human fossils fulfilled the above criteria for amino acid preservation and were thus expected to contain endogenous DNA (; ). These samples were geographically well distributed across Europe () and included remains whose morphology is typical of Neandertals (e.g., La Chapelle-aux-Saints) and of modern humans (La Madeleine, Cro-Magnon). They also included samples that have sometimes been considered “transitional” between Neandertals and modern humans, based on their morphological features: Vindija (
Smith 1984) and Mladeč (
Frayer 1986,
Frayer 1992;
Wolpoff 1999).
| Table 1DNA Retrieved from Late Pleistocene Fossils in This Study |
If low amounts of DNA are preserved in a specimen, some extracts will fail to contain DNA molecules by chance (
Hofreiter et al. 2001a). Therefore, except in the case of Mladeč 2, in which the amount of material available permitted only two extractions, we extracted each of the four Neandertal and the five early modern human samples three times. For each extraction, amplifications were performed using two primer pairs: (i) “hominoid primers” that amplify homologous mtDNA sequences from the previously determined Neandertals and contemporary modern humans, as well as African great apes; (ii) “Neandertal primers” that, under the conditions used, amplify only Neandertal mtDNAs even in the presence of a large excess of modern human DNA (
Krings et al. 2000;
Schmitz et al. 2002). Since authentic ancient DNA is typically highly degraded, both primer pairs were designed to amplify short mtDNA fragments (72 and 31 bp, respectively, excluding primers). In each of these fragments, two substitutions allow the discrimination of previously determined Neandertal mtDNA sequences from contemporary modern human sequences. The sensitivity of both primer pairs is similar, as shown by the fact that they are both able to amplify single template molecules as judged from nucleotide misincorporation patterns (
Hofreiter et al. 2001a). In order to determine the nature of the DNA sequences amplified, each amplification product was cloned and approximately 30 clones were sequenced for each “hominoid product” and ten clones for each “Neandertal product.”
When amplified with the hominoid primers, all Neandertal and all early modern human remains yielded modern human DNA sequences (see ). In addition, five cave bear teeth from Vindija, Croatia, and one from Gamssulzen, Austria, extracted in parallel with the hominid samples, all yielded human sequences. This confirms previous results in showing that most, if not all, ancient remains yield human DNA sequences when amplification conditions that allow single DNA molecules to be detected are used (
Hofreiter et al. 2001b). For three Neandertal and all five modern human remains, several different mtDNA sequences were retrieved from individual extractions, and in the case of one Neandertal and one modern human, at least two of the sequences were also found in an independent extraction from the same specimen. Additionally, one of the cave bear teeth yielded a human sequence found in two independent extracts. Thus, the fact that a DNA sequence is found in two independent extracts is a necessary, but not sufficient, criterion of authenticity when human remains are analyzed. This implies that in the absence of further technical improvements, it is impossible to produce undisputable human mtDNA sequences from ancient human remains. In addition to DNA sequences identical to those previously amplified from present-day humans, the Neandertal bones Vi-77 and Vi-80 from Vindija yielded four out of 89 and 73 out of 85 mtDNA sequences, respectively, that were identical to previously determined Neandertal sequences. Thus, these two specimens contain a proportion of Neandertal-like mtDNA sequences (i.e., sequences that carry two substitutions that differentiate Neandertal mtDNA sequences from modern human mtDNA sequences as described above) that is high enough to detect using primers that amplify also modern human DNA.
When amplified with Neandertal-specific primers, Neandertal-like mtDNA sequences were amplified from two independent extractions from all Neandertal fossils (see ; ). For one of these, Vi-80 from Vindija, DNA preservation was sufficient to allow the retrieval of longer fragments and thus the reconstruction of 357 bp of the hypervariable region I (see Supporting Information section;
Figure S1). This mtDNA sequence was identical to that retrieved from another bone from the same locality (Vi-75;
Krings et al. 2000). In contrast to the Neandertal remains, none of the early modern human extracts yielded any amplification products with the Neandertal primers, although these remains are similar in chemical preservation to the Neandertal remains (see ).
Thus, all Neandertal remains analyzed yielded mtDNA sequences that are not found in the human mtDNA gene pool today but are similar to those found in four previously published Neandertals (
Krings et al. 1997,
Krings et al. 2000;
Ovchinnikov et al. 2000;
Schmitz et al. 2002) (see ). This is compatible with results suggesting that the extent of Neandertal mtDNA diversity was similar to that of current humans and lower than that of the great apes (
Krings et al. 2000;
Schmitz et al. 2002). It is noteworthy that this result is not an artifact created by discarding “modern-like” mtDNA sequences amplified from Neandertals (
Trinkaus 2001), since all Neandertal remains with good biomolecular preservation yield “Neandertal-like” mtDNA sequence. Furthermore, none of the five early modern humans yields “Neandertal-like” mtDNA sequences in spite of the fact that these remains are as well preserved in terms of amino acids as the Neandertal remains. Thus, we fail to detect any evidence of mtDNA gene flow from Neandertals to early modern humans or from early modern humans to Neandertals.
However, a relevant question is what extent of gene flow between Neandertals and early modern humans the current data allow us to exclude. In this regard, it is of relevance that the five early modern humans analyzed lived much closer in time to the Neandertals than do contemporary individuals. The probability that mtDNA sequences potentially contributed to modern humans by Neandertals were lost by drift (
Nordborg 1998) or swamped by continuous influx of modern human mtDNAs (
Enflo et al. 2001) in the Neandertal gene pool is therefore much smaller than when contemporary humans are analyzed (e.g.,
Relethford 1999). In fact, the five early modern humans analyzed almost double the amount of information about the Upper Pleistocene mtDNA gene pool since, under a model of constant effective population size, all contemporary humans trace their mtDNA ancestors back to only four to seven mtDNA lineages 20,000 to 30,000 years ago (A;
Figure S2), while all other mtDNA sequences present in the gene pool at that time have been lost by random genetic drift. Since the probability is very low (
p < 0.007) that one or more of the five early modern humans analyzed here are among these few ancestors of current humans, the five Upper Pleistocene individuals can be added to the ancestors of the current mtDNA gene pool to allow us to ask what extent of Neandertal mtDNA contribution to early modern humans can be statistically excluded using the coalescent. Under the model of a constant human effective population size (
Tavare 1984;
Nordborg 1998) of 10,000 over time (A), any contribution of Neandertal mtDNA to modern humans 30,000 years ago larger than 25% can be excluded at the 5% level (
Figure S3). A more realistic scenario may be that the spread of modern humans was accompanied by an increase in population size before and during their migration out of Africa and subsequent colonization of western Eurasia (see B). In that case, the Neandertal contribution that can be excluded is smaller (i.e., less gene flow could have taken place), but that depends critically on when and how the expansion occurred. Finally, under the unlikely scenario that population size was constant during the migration out of Africa and colonization of Europe and expanded only after a putative merging with Neandertals, the Neandertal contribution could have been larger, but this also depends on the nature of the growth (see C).