Skip to main content

Genotyping human ancient mtDNA control and coding region polymorphisms with a multiplexed Single-Base-Extension assay: the singular maternal history of the Tyrolean Iceman



Progress in the field of human ancient DNA studies has been severely restricted due to the myriad sources of potential contamination, and because of the pronounced difficulty in identifying authentic results. Improving the robustness of human aDNA results is a necessary pre-requisite to vigorously testing hypotheses about human evolution in Europe, including possible admixture with Neanderthals. This study approaches the problem of distinguishing between authentic and contaminating sequences from common European mtDNA haplogroups by applying a multiplexed Single-Base-Extension assay, containing both control and coding region sites, to DNA extracted from the Tyrolean Iceman.


The multiplex assay developed for this study was able to confirm that the Iceman's mtDNA belongs to a new European mtDNA clade with a very limited distribution amongst modern data sets. Controlled contamination experiments show that the correct results are returned by the multiplex assay even in the presence of substantial amounts of exogenous DNA. The overall level of discrimination achieved by targeting both control and coding region polymorphisms in a single reaction provides a methodology capable of dealing with most cases of homoplasy prevalent in European haplogroups.


The new genotyping results for the Iceman confirm the extreme fallibility of human aDNA studies in general, even when authenticated by independent replication. The sensitivity and accuracy of the multiplex Single-Base-Extension methodology forms part of an emerging suite of alternative techniques for the accurate retrieval of ancient DNA sequences from both anatomically modern humans and Neanderthals. The contamination of laboratories remains a pressing concern in aDNA studies, both in the pre and post-PCR environments, and the adoption of a forensic style assessment of a priori risks would significantly improve the credibility of results.


It is often hoped that the fields of ancient DNA (aDNA) and population genetics will converge to provide novel insights into human prehistory, but there have been very few convincing aDNA studies conducted on human remains, due to the difficulty of differentiating between authentic results and contaminating DNA [14]. This problem is particularly acute when the putative endogenous sequences are similar to possible contaminants, and for this reason the most credible research has been restricted to the investigation of remains with highly distinctive, individual or population-specific, genetic signatures, e.g. [58]. Consequently, a number of high-profile human aDNA studies have been questioned on the grounds of insufficient proof of authenticity [1, 912].

Traditionally, the molecular tool for the identification of both endogenous human DNA and potential contaminant sequences has been the comparatively high levels of polymorphism found in the first and second hypervariable sections (HVS1, HVS2) of the mitochondrial (mtDNA) control region. The rapid rate of evolution in this short stretch of the molecule increases the likelihood of encountering multiple distinctive single nucleotide polymorphisms (SNPs) in single overlapping polymerase chain reaction (PCR) products. However, in practice, the utility of the control region for this purpose is often compromised by elevated levels of homoplasy, where identical sequences can be shared between different haplotypes. This has lead to an increasing use of hierarchical coding region SNPs to better discriminate between separate haplotypes with identical HVS1 sequences [13]. Whilst this is relatively unproblematic for modern DNA samples, in the case of ancient samples, such as the Iceman, the production of control and coding region sequences in separate PCRs makes it difficult to unequivocally rule out contamination as the source of individual results [12, 14, 15]. So, paradoxically, the most credible examples of this approach are those where the control region motif is already distinctive in its own right, and the provision of a single additional SNP often provides no significant additional resolution, as with the five Amerindian mtDNA haplogroups, e.g. [16, 17].

The most critical factor affecting the accurate genotyping of mtDNA data comprised of both SNPs and control region sequences is the reliability of interpreting unlinked PCR products. These may be subject to sporadic contamination and differential amplification of endogenous and contaminant DNA. One way to address this issue for SNPs is to simultaneously characterise multiple polymorphic sites in a two-stage multiplex reaction using a linear Single-Base-Extension (SBE) mini-sequencing reaction [1822]. These can overcome many of the problems inherent in providing phylogenetic resolution for human aDNA samples and, in principle, may be able to differentiate between endogenous DNA and a possible contaminant from the same mtDNA haplogroup. Here, we extend this principle to genotype both control region polymorphisms and hierarchical coding regions SNPs within a single multiplexed SBE reaction (K-plex). This is used to obtain a detailed genotype of the Tyrolean Iceman's mtDNA using three separate extracts made from samples of bone and faecal material. The robustness of the results are explored through controlled contamination experiments using SBE reactions and quantitative PCR, and by providing an evaluation of the a priori possibility of obtaining an identical genotype from contaminants originating within the laboratory. The discussion includes a detailed reanalysis of the original study of the Iceman's mtDNA [14] and the conclusions serve to highlight the problems of achieving authentic results from European aDNA samples in general, and the fallibility of independent replication, in particular.


Cloned Sequences

The results from initial PCR amplification and direct sequencing of the three new DNA extracts are reported in Figure 1 and Table 1 and characterise the Iceman's diagnostic control region polymorphisms. The first hyper-variable section (HVS1) between nucleotide positions (nps) 16040 and 16401, revealed three consistent polymorphisms (16224C-16311C-16362C), relative to the revised Cambridge Reference Sequence (rCRS) [23, 24], in all clones from all samples (see additional file 1 showing clones covering the key sites nps 16311 and 16362 from all three extracts). There were no positive results obtained from the corresponding extraction blanks and negative controls used in each amplification. The results suggest that the mtDNA belongs to K1* by virtue of matching the revised Cambridge Reference Sequence (rCRS) [24] at nps 497 (K1a), 498 (K1c) and 5913 (K1b), whilst being derived at nps 1189 and 10398 (both K1). Importantly, the presence of the transition at np 8137 is unique within contemporary hg K [13], and the mutation at np 16362 was not reported in the original study of the Iceman's mtDNA [14].

Figure 1

PCR results for the three extracts of the Iceman's DNA (HB49, HB50, BPE1). The results place the Iceman's mtDNA into a new (undefined) phylogenetic position (K1*) whilst also reporting the transitions at nps 8137 and 16362. Human mtDNA mutations are scored relative to the revised Cambridge Reference Sequence (rCRS) [23, 24], so the Tyrolean Iceman's mtDNA is 'ancestral' at the three positions (C497T, G5913A, 498delC) that define K1a, K1b, and K1c. The sites tested in, 1 [14], 2 [15] and 3this study, are shown to compare the phylogenetic discrimination achieved. Critically, the 16362C polymorphism was not reported in the original study [14] and not tested for in [15]. A recent study [25] explored the entire mtDNA genome of the Iceman by singleplex PCR and has reported all of the polymorphisms listed in Table 1.

Table 1 Singleplex PCR results for the bone (BPE1, HB50) and faecal (HB49) extracts of the Iceman's mtDNA, showing all polymorphisms obtained in separate sets of cloned sequences.

K1 SBE multiplex

Due to the difficulty inherent in proving that results matching the reference sequence at nps 497, 498 and 5913 were not due to contamination [12], a multiplex SBE assay (K1-plex) was designed to test for these three HVS1 polymorphisms (16224C-16311C-16362C), plus nine hierarchical SNPs from both the control and coding regions. The multiplex results revealed transitions at nps 1189, 1811, 9698, 10398 and 14167 (Table 2), indicating membership of the human mtDNA sub-haplogroup K1 [13], as previously suggested by [15, 25]. In addition, because the Iceman's sequence matches the rCRS at nps 497, 498, and 5913 (diagnostic of K1a, K1c, and K1b, respectively), we can extend previous suggestions that the Iceman's genotype does not belong to any of the three known clades of K1, and should be referred to a new (undefined) paraphyletic clade, of K1 [25] (Figure 1). A transition at np 8137, previously discovered by singleplex assay of the Iceman's mtDNA, was confirmed by the multiplex SBE results, and is currently unique within hg K [13]. All multiplex results were cross-checked by cloned single-plex PCR results to investigate the potential for false positives in the SBE reaction. The allelic state of each of the 11 SNPs matched the multiplex results (details of all singleplexes are given in Table 1).

Table 2 Multiplex Single Base Extension results for the HB50 bone extracts of the Iceman's DNA using the K1 multiplex.

Crucially, the actual SBE traces (Figure 2) show that the results originated from a single effective population of DNA templates because there were no additional alleles at any of the sites tested. The use of multiplexed SBE also ensured that those hg K1 diagnostic sites matching the rCRS (nps 497, 498, 5913) were not a result of contamination, or sample mix-up, because all were typed concurrently in a single assay, together with other SNPs (16362C, 8137T) diagnostic of the Iceman's mtDNA. Real-time PCR quantification (see additional file 2) indicated that the multiplex assay was initiated with the equivalent of ~36,000 molecules of mtDNA capable of amplifying 60 bp and that inhibition was not detected. This effective copy number is consistent with the absence of repetitive damage motifs amongst the clones and the lack of allelic drop-out in the SBE results [26].

Figure 2

Representative electropherogram of the K1 multiplex assay conducted on the HB50 aDNA extract made from the Iceman's bone, reporting all three HVS1 mutations revealed in the singleplex assays (16224–16311–16362) and eight additional phylogenetically informative SNPs. The SBE results are differentiated by a combination of the alleles and their electrophoretic mobility, determined by the length of the primers (see Table S4). All alleles are called relative to the strand of origin targeted in the SBE reaction; those from the reverse strand should be complemented to be consistent with the rCRS (see Table S2). Positions 5913 and 498 (the smaller peak of the two) have coalesced slightly due to departures from the predicted electrophoretic mobility of the SBE primers. Minor peaks to the left and immediately adjacent of some SNP positions represent background noise, due to variation in length of some SBE primers and can be ignored. There is some background noise evident on the A system but only the peak at the 16362 position is in the correct position to represent contamination (see Fig. 3).

Contamination controls

To further test the potential effect of contamination on the results, a second K1 multiplex assay was performed with the addition of various quantities of a second DNA extract (hg J) to imitate the role of modern contaminant. This produced clearly observable double peaks (representing different alleles at the same site) across a large range of dilutions (Figure 3), confirming the presence of two populations of templates. The lowest quantity of contaminant detected by the multiplex assay was measured by real-time quantitative PCR to be equivalent to ~4% of the total starting copy number. As the proportion of the contaminant increases, so do the secondary peaks, but even at levels of ~20% it is still possible to unambiguously distinguish the two separate populations of templates. This result strongly supports the conclusion from the cloned HVS1 sequences (Figure S1) that the same single endogenous population of templates is represented in the three extracts of the Iceman's DNA.

Figure 3

Results of the contamination assay showing the G/A system from three different multiplex SBE reactions (MP1, MP2, MP3), initiated with the Iceman's DNA mixed with three different amounts of an hg J contaminant. Because the contaminant has ancestral alleles at nps 16224, 16311 and 16362, this causes additional peaks to appear in the A system (highlighted). The small peak to the left of 16362A in each case is background noise (i.e. not due to primer length variation which would produce multiple peaks of varying size, see the G system in Figure 2. This does not interfere with the base calling because it is outside the range of either the 00497 or 16362 SBE primers. The sensitivity of this method is demonstrated by its ability to detect the second population of templates (hg J) at the level of ~4%; even at levels of ~20% contamination, the results for the Iceman's mtDNA are unambiguous and unaffected by the contaminant. The ability to test multiple sites in a single reaction, in the presence of contaminants, without the need to clone individual PCR products, represents a significant advantage of the SBE methodology.

We also typed the mtDNA of all workers involved in the study, from the removal of the samples from the mummy in Italy to the extraction of the DNA in Oxford and no individual had mtDNA belonging to hg K. Although positive results from extraction blanks and PCR negative controls did not occur during the present study, a human contaminant was identified amongst earlier work conducted in the Oxford isolation facility. This was from a modern source of hg K mtDNA with the control region motif 16224C-16311C (i.e. without 16362C) and was eliminated by modification to laboratory protocols regarding common consumables (see Methods).

This hg K contaminant, together with all aDNA extracts with putative European mtDNA prepared in the isolation facility, is listed in additional file 3 but these sequences include no examples of hg K with 16362C. Therefore, at the time that the research into the Iceman's mtDNA was conducted, the probability of obtaining a matching hg K contaminant in the Henry Wellcome Ancient Biomolecules Centre appears to have been extremely low.


The significance and authenticity of the transition T16362C

Revisiting the mtDNA of the Tyrolean Iceman has raised the issue of how to differentiate between conflicting but highly similar sets of aDNA results, in this case from independent studies of the same individual [[14, 15, 25], this study]. Although each study obtained hg K mtDNA in their PCR products, only two reported the 16362C polymorphism [[25], this study]. Therefore, there would appear to be multiple sources of hg K mtDNA templates involved because of this difference in the control region motifs. It is important to note that this single base change (T>C) cannot derive from post-mortem damage to the mtDNA molecule because this results in C>T changes only, from the deamination of cytosine to uracil, which is read as a thymine, and appears as a G>A on the opposite strand [27]. Therefore, an erroneous T>C transition could only derive from an artefact of PCR, but the chances of this occurring in multiple amplifications from our three samples, and for the same artefact becoming fixed in every PCR product, appear to be remote. The lack of variation at all the sites we tested, both in cloned sequences and SBE assays, is only consistent with the presence of a single population of mtDNA templates from all three extracts with the control region motif of 16224C-16311C-16362C, for which the a priori probability of these deriving from a common contaminant, either of the laboratory or the sampling procedure, has been demonstrated to be extremely low.

Phylogenetic evidence of authenticity

The likelihood of the current results representing the endogenous DNA of the Iceman is substantially increased by the haplotype falling outside of K1a/K1b/K1c, combined with the unique hg K transition C8137T, matching the results of [25]. A Genbank search reveals two contemporary coding region genomes (SO44 from [28] and Oc05 from [29]) that can be securely assigned to the same paraphyletic K1* position as the Iceman and which share three coding region polymorphisms (see Figure 4). Only one, a Dutch sample (SO44 from [28]) has a published control region sequence, but this also displays the T16362C transition. We propose that these two genomes define a new clade K1d defined by the coding region motif of 4856C-5592G-8901G. This finding provides some support for the proposal of a single European paraphyletic clade of mtDNA hg K1 defined by the 16362C polymorphism [25]. However, the use of a single mutational hotspot (i.e. np 16362, which occurs at least eight times on the hg K background alone) for this purpose is usually discouraged [30]. If all these genomes belong to a single clade, then convention suggests that it should be named K1d, following the cladistic rules for the hierarchical ordering of haplogroups and sub-haplogroups formalised by [31].

Figure 4

Phylogenetic tree of the main branches of human mtDNA hg K1 showing the western European paraphyletic lineages falling outside of K1a, K1b and K1c. All mutations are transitions unless indicated by a letter. Dotted lines indicate samples for which replicated whole genome data is not available; SNPs denoted in bold are those tested in the K1 multiplex; the tree is rooted in African hg L3. All lineages with 16362C are grouped by parsimony into a single potential European paraphyletic clade. The clade of K1d is defined by three coding region mutations shared between SO44 and Oc5 from [29]. Together with a single example of K1* from Iran [13] (not shown), these represent the sum total of hg K1 mtDNA lineages currently known to fall outside of K1a\K1b\K1c, amongst a worldwide survey of over 1000 haplogroup K samples (1Howell et al., 2003 [28]; 2this study; 3 Genbank direct submission).

The perils of contamination

The difference between the conclusions of the original study and the current consensus of the Iceman's mtDNA control region might appear to be trivial. Yet the implications are actually of the utmost importance and warrant a detailed reanalysis of the original data. This is because the absence of the 16362C allele amongst the results of Handt et al [14], suggests that they genotyped a different hg K individual. Reanalysis of the original data reveals the presence of at least two populations of templates (Figure 5a). One of these is a variant of human mtDNA hg J (defined by the motif 16069T-16126C-16390A) [21, 32], whilst the other observable haplotype (16093C-16224C-16311C) is a common occurrence within hg K [13].

Figure 5

a/b Graphical representation of the sequences obtained for the Iceman from Handt et al. 1994. (a) Cloned HVS1 control region sequences from [14] (the nucleotide positions are listed in the top row, whilst the first column refers to the number of each clone). There are two clearly identifiable European mtDNA haplotypes belonging to hg K (16093C-16224C-16311C, clones 1–11 and potentially 19–43), and hg J (16069T-16126C-16390A, clones 12–16). The absence of both 16093C and 16126C from the shortest (100 bp) amplification product (clones 44–55) indicates the presence of at least three populations of templates. Positions of template switching ('jumping PCR') are indicated by brown highlighting, and evident by the replacement of the expected mutation either by the rCRS (no letter), or a polymorphism from a different template. Critically, there are no examples of the 16362C polymorphism in any templates covering this position. (b) Results from the independent replication [37] which were only indirectly sequenced (not cloned), identifying just the two control region transitions basal to hg K (16224C-16311C). This suggests that these results derived from a second hg K contaminant, differing from those of Figure 5a through the lack of the 16093C transition, and from the Tyrolean Iceman by the absence of the key T>C transition at np 16362.

Appropriate Molecular Behaviour

The published consensus sequence excluded the transition detected at np 16093 on the basis that it did not feature in the 100 bp product [14]. This short sequence was argued to have been more likely to represent the endogenous DNA on the basis of appropriate molecular behaviour [14]. This concept predicts that the shortest amplification product, in the presence of modern contaminants, is more likely to represent the endogenous DNA [33, 34] (but see also [35]). As we have demonstrated that the mtDNA of our Iceman extracts does not contain the 16093C allele, it is conceivable that the 100 bp product of Handt et al. [14] did originate from DNA endogenous to the sample. Support for, or against, this conclusion can be sought from the clones for sequences including np 16362. It was recently suggested that the lack of the 16362C allele amongst the Handt et al. sequences may partly be attributable to the production of chimeras amongst these longer amplicons [25]. The obvious mixture of alleles at nps 16093, 16311 and 16390 (clones 1, 4, 6, 7, 10, 12, 16, 17, 19–21, and 24–26; Figure 5a) demonstrate that either these sites were subject to post-mortem damage or hg K mtDNA templates were involved in 'jumping PCR' events with contaminants [36]. Yet, of all the clones covering the upper portion of the control region sequence (1–17, and 28–34), none contain the key transition at np 16362. This suggests that the extraction actually contained a different source of hg K mtDNA to the single mtDNA sequence recovered in [25] and the present study.

The fallibility of independent replication

Given the absence of 16093C and 16362C amongst the replication results of Handt et al., the only plausible conclusion is that these derived from a second, entirely different, hg K source. These results, combined with a unique haplotype in [25] and the present study, strongly suggests that both laboratories in the original study suffered from hg K contamination. The mtDNA profiles of the staff in the primary laboratory of the first study were not provided [14]; although a subsequent publication of the replicated results [37] disclosed that a member of the secondary laboratory staff was hg K (16224C-16311C), substantially increasing the likelihood for contamination occurring. It is important to note that the independent replication step, in this case, did not technically fail. Instead, the analysis did not manage to distinguish between these two non-identical sources of hg K mtDNA, neither of which appear to have derived from the Iceman. The fallibility of the original, independently replicated, study emphasizes why unlinked singleplex PCR results may not be sufficient proof of authenticity. This is particularly problematic where novel aDNA genomes are defined by long stretches of sequence matching the rCRS, and motifs including C>T changes that can arise as a product of post-mortem modifications to DNA templates [25].

Recommendations for future human aDNA research

Criteria of authenticity

We suggest that any contamination risk assessment in human aDNA research would be improved through the incorporation of an analysis of the a priori probability that the extract will be contaminated. Human aDNA studies are now routinely providing genetic characterisation of individuals who have been involved in the manipulation of a sample, e.g. [17]. For an accurate evaluation of the true risk, in a way that would satisfy common international standards in forensic science, it is also necessary to supply data on false positives and other previous amplification of similar genotypes from the laboratories concerned. The presence of hg K mtDNA amongst the common reagents of the laboratory in the present study confirms how consumables can be a source of contamination [38], with obvious implications for the efficacy of independent replication. The findings of the present study provide clear evidence why, as with any scientific investigation, all reasonable efforts to falsify human aDNA results should be made, even when controls such as independent replication are known to be fallible under certain circumstances [39]. This is clearly the case for human aDNA research carried out on European remains in European laboratories.

Identifying DNA haplotypes using multiple SNPs

The use of the K1-plex demonstrates the potential power of discrimination achievable by linking control and coding region sites in a single SBE multiplex. Although the fortuitous circumstances of a unique maternal lineage cannot be expected in the majority of human aDNA studies, this kind of multiplexed approach should in many situations be able to identify where DNA templates from two or more discrete individuals are present. So, the standards applied here should not be taken as a reference for other research using human DNA extracts that fall squarely into more mainstream situations. Nevertheless, the linking of all loci in one multiplex SBE reaction can still improve the security of interpretation, whilst assisting in the identification of possible contamination events and helping to conserve irreplaceable aDNA extracts.


This study has highlighted the potential fallibility of the interpretation of human aDNA results, and the need for a forensic style assessment of risks associated with the production of data. Independent replication remains a core component of authentication strategies for human aDNA studies, and it is paramount that this test is applied critically to function properly. More stringent controls are required to move the field forward, both a priori and a posteriori, for the production and interpretation of data; in particular, the disclosure of all possible sources of contamination in aDNA studies using generic hominin primers. In PCR-based studies, this data should include genetic profiles of the workers involved for whichever target loci are being investigated.

The use of multiplex SBE assays offer one solution for authenticating human aDNA results and screening novel genome data. In the present study the K1 multiplex accurately genotyped the Iceman's mtDNA and has independently confirmed it belongs to a novel lineage of mtDNA hg K. A novel aspect of this approach is the ability to link SNPs from both the control and coding regions of the mtDNA genome in a single reaction, reducing the possibility of contamination affecting the results. In a controlled contamination experiment, the SBE methodology correctly identified the endogenous genotype of the Iceman and the relative quantities of the modern mtDNA contaminant.


Two samples of bone removed from the ileum, and one of faecal matter from the rectum were digested in temporally separated extractions (4 ml total volume) by a previously published protocol [22]. One extraction blank was produced for each sample, and PCRs also included negative controls.

The specialist Ancient Biomolecules Centre (ABC) at Oxford was used to extract DNA, and set up PCR reactions. The ABC is physically isolated, subject to stringent anti-contamination procedures, equipped with positive air pressure and UV lighting, and has a DNA laboratory and equipment (glove box, instruments, full body suits, protective masks, etc.) dedicated solely to ancient human specimens. The thermal cycling reactions and subsequent downstream work took place in a separate laboratory located in the Department of Zoology.

Initial singleplex PCRs on the control region were conducted with primers listed in additional file 4 and cloned by previously published protocols [22]. Several hierarchical SNPs from the human mtDNA phylogeny together with haplogroup K specific SNPs were identified [13, 40] using additional singleplex PCR primers designed to be compatible under multiplex PCR conditions [26].

The three control region SNPs, located in the singleplex PCRs, were combined with key SNPs from major clades of the hg K phylogeny, to form the K1 multiplex. A selection of higher order hierarchical SNPs from the human phylogeny were also included to check the consistency of results across additional loci [13, 40]. This integration of all diagnostic sites into a single assay reduces the range of explanations required when results do not conform to expectations because the possibility of sample mix up is eliminated.

The initial multiplex PCR was performed using 2 μ l of DNA extract in a 25 μ l reaction volume using primers listed in additional file 5 by a previously published protocol [22]. SBE reactions were performed in a total volume of 5 μ l with 0.5 μ l purified PCR product, 2.5 μ l SNaPshot™ reaction mix (AB) and 0.5 μ l SBE primer mix (SBE primers listed in additional file 6). Negative controls were included in both stages of the assay to assess the potential for contamination.

Automated allele calls were made using macros constructed in Genotyper 3.7 (AB). Primer concentrations were adjusted to give balanced amounts of products and peak heights, concentrating on the SBE reaction [19]. All primers were synthesized and reverse-phase HPLC purified by GmbH (Ulm, Germany). A detailed protocol for the optimization of SBE multiplexes is available elsewhere [26].

Quantitative PCR (qPCR) was conducted on the HB50 aDNA bone extract and a modern DNA sample (belonging to the mitochondrial hg J) in 25 μ l reaction volumes, using SYBR Green PCR master-mix (AB), 300 nM human-specific primers, synthesized template standards (see additional file 7) and an ABI Prism 7000 Sequence Detection System (95°C for 10 min, followed by 40 cycles of 95°C for 15 sec, and 60°C for 60 sec). All samples were run in triplicate with two 1:4 serial dilutions. The standards were also run in triplicate with 1:9 dilutions ranging from 106-102 copies per μ l. The adjusted R2 value of all assays performed was equal to a minimum of 99%, and PCR efficiency, in each case, was calculated with the formula [41] to be ~2 (see additional file 2).

Once the effective copy number of the modern DNA was known, a range of dilutions were run against the HB50 extract in a relative quantification, using the same primers and conditions, until the reported copy number of both the ancient and modern extracts were equivalent. The multiplex contamination experiments were then performed using successive serial dilutions of the hg J contaminant mixed with constant quantities of the aDNA extract to seed the multiplex PCR. This procedure provides a detailed measure of the sensitivity of the SBE methodology for identifying and estimating the amount of contamination present in an aDNA extract through the incorporation of additional peaks at segregating sites.

All contamination experiments, thermal cycling reactions and downstream work took place in a separate laboratory located in the Department of Zoology. All relevant personnel in this study had their mtDNA typed. The two researchers who took the samples from the Iceman have mtDNA matching hg H and hg J1b. All aDNA work in Oxford (extractions, singleplex and multiplex PCRs) were carried out by PE, who is hg H.

Figures were prepared using Inkscape Vector Illustrator compiled for use in Ubuntu Linux. The real-time PCR plot was created using the R statistical package The manuscript was written using the LyX document preparation system


  1. 1.

    Bandelt HJ: Mosaics of ancient mitochondrial DNA: positive indicators of nonauthenticity. Eur J Hum Genet. 2005, 13 (10): 1106-1112. 10.1038/sj.ejhg.5201476.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Caramelli D, Lalueza-Fox C, Vernesi C, Lari M, Casoli A, Mallegni M, Chiarelli B, Dupanloup I, Bertranpetit J, Barbujani G, Bertorelle G: Evidence for a genetic discontinuity between Neandertals and 24,000-year-old anatomically modern Europeans. Proc Nautl Acad Sci USA. 2003, 100 (11): 6593-6597. 10.1073/pnas.1130343100.

    CAS  Article  Google Scholar 

  3. 3.

    Pääbo S, Poinar H, Serre D, Jaenicke-Despres V, Hebler J, Rohland N, Kuch M, Krause J, Vigilant L, Hofreiter M: Genetic analyses from ancient DNA. Annu Rev Genet. 2004, 38: 645-679. 10.1146/annurev.genet.37.110801.143214.

    Article  PubMed  Google Scholar 

  4. 4.

    Serre D, Langaney A, Chech M, Teschler-Nicola M, Paunovic M, Mennecier P, Hofreiter M, Possnert G, Pääbo S: No evidence of Neandertal mtDNA contribution to early modern humans. PLoS Biol. 2004, 2 (3): E57-10.1371/journal.pbio.0020057.

    PubMed Central  Article  PubMed  Google Scholar 

  5. 5.

    Ricaut FX, Fedoseeva A, Keyser-Tracqui C, Crubezy E, Ludes B: Ancient DNA analysis of human neolithic remains found in northeastern Siberia. Am J Phys Anthropol. 2005, 126 (4): 458-462. 10.1002/ajpa.20257.

    Article  PubMed  Google Scholar 

  6. 6.

    Kemp B, Malhi R, McDonough J, Bolnick D, Eshleman J, Rickards O, Martinez-Labarga C, Johnson J, Lorenz J, Dixon E, Fifield T, Heaton T, Worl R, Smith D: Genetic analysis of early holocene skeletal remains from Alaska and its implications for the settlement of the Americas. Am J Phys Anthropol. 2007, 132: 605-621. 10.1002/ajpa.20543.

    Article  PubMed  Google Scholar 

  7. 7.

    Barik SS, Sahani R, Prasad BVR, Endicott P, Metspalu M, Sarkar BN, Bhatacharya S, Annapooma PCH, Sreenath J, Sun D, Sanchez JJ, Ho SYW, Chandrasekar A, Rao VR: Detailed mtDNA Genotypes Permit a Reassessment of the Settlement and Population Structure of the Andaman Islands. Am J Phys Anthropol. 2008, 136: 19-27. 10.1002/ajpa.20773.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Gilbert MT, Kivisild T, Grønnow B, Andersen PK, Metspalu E, Reidla M, Tamm E, Axelsson E, Götherström A, Campos PF, Rasmussen M, Metspalu M, Higham TF, Schwenninger JL, Nathan R, De Hoog CJ, Koch A, Møller LN, Andreasen C, Meldgaard M, Villems R, Bendixen C, Willerslev E: Paleo-Eskimo mtDNA genome reveals matrilineal discontinuity in Greenland. Science. 2008, 320: 1787-1789. 10.1126/science.1159750.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Cooper A: Reply to Stoneking: ancient DNA-how do you really know when you have it?. Am J Hum Genet. 1997, 60 (4): 1001-1003.

    PubMed Central  CAS  PubMed  Google Scholar 

  10. 10.

    Cooper A, Rambaut A, Macaulay V, Willerslev E, Hansen AJ, Stringer C: Human origins and ancient human DNA. Science. 2001, 292 (5522): 1655-1656. 10.1126/science.292.5522.1655.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Gilbert MTP, Barnes I, Collins MJ, Smith C, Eklund J, Goudsmit J, Poinar H, Cooper A: Long-term survival of ancient DNA in Egypt: response to Zink and Nerlich (2003). Am J Phys Anthropol. 2005, 128: 110-114. 10.1002/ajpa.20045.

    Article  PubMed  Google Scholar 

  12. 12.

    Endicott P, Sanchez JJ, Metspalu E, Behar DM, Kivisild T: The Unresolved Location of Ötzi's mtDNA Within Haplogroup K. Am J Phys Anthropol. 2007, 132 (4): 590-591. 10.1002/ajpa.20501.

    Article  PubMed  Google Scholar 

  13. 13.

    Behar DM, Metspalu E, Kivisild T, Achilli A, Hadid Y, Tzur S, Pereira L, Amorim A, Quintana-Murci L, Majamaa K, Herrnstadt C, Howell N, Balanovsky O, Kutuev I, Pshenichnov A, Gurwitz D, Bonne-Tamir B, Torroni A, Villems R, Skorecki K: The matrilineal ancestry of Ashkenazi Jewry: portrait of a recent founder event. Am J Hum Genet. 2006, 78 (3): 487-497. 10.1086/500307.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  14. 14.

    Handt O, Richards M, Trommsdorff M, Kilger C, Simanainen J, Georgiev O, Bauer K, Stone A, Hedges R, Schaffner W, Utermann G, Sykes B, Pääbo S: Molecular genetic analyses of the Tyrolean Ice Man. Science. 1994, 264 (5166): 1775-1778. 10.1126/science.8209259.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Rollo F, Ermini L, Luciani S, Marota I, Olivieri C, Luiselli D: Fine characterization of the Iceman's mtDNA haplogroup. Am J Phys Anthropol. 2006, 130 (4): 557-564. 10.1002/ajpa.20384.

    Article  PubMed  Google Scholar 

  16. 16.

    Stone A, Stoneking M: mtDNA analysis of a prehistoric Oneota population: implications for the peopling of the New World. Am J Hum Genet. 1998, 62: 1153-1170. 10.1086/301838.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  17. 17.

    Gilbert M, Jenkins D, Götherstrom A, Naveran N, Sanchez J, Hofreiter M, Thomsen P, Binladen J, Higham T, Yohe R, Parr R, Cummings L, Willerslev E: DNA from Pre-Clovis Human Coprolites in Oregon, North America. Science. 2008, 320:

    Google Scholar 

  18. 18.

    Syävnen AC: From gels to chips: "minisequencing" primer extension for analysis of point mutations and single nucleotide polymorphisms. Hum Mutat. 1999, 13: 1-10. 10.1002/(SICI)1098-1004(1999)13:1<1::AID-HUMU1>3.0.CO;2-I.

    Article  Google Scholar 

  19. 19.

    Norton N, Williams NM, Williams HJ, Spurlock G, Kirov G, Morris DW, Hoogendoorn B, Owen MJ, O'Donovan MC: Universal, robust, highly quantitative SNP allele frequency measurement in DNA pools. Hum Genet. 2002, 110: 110-478. 10.1007/s00439-002-0706-6.

    Article  Google Scholar 

  20. 20.

    Coble MD, Just RS, O'Callaghan JE, Letmanyi IH, Peterson CT, Irwin JA, Parsons TJ: Single nucleotide polymorphisms over the entire mtDNA genome that increase the power of forensic testing in Caucasians. Int J Legal Med. 2004, 118: 137-146. 10.1007/s00414-004-0427-6.

    Article  PubMed  Google Scholar 

  21. 21.

    Brandstätter A, Salas A, Niederstätter H, Gassner C, Carracedo A, Parson W: Dissection of mitochondrial superhaplogroup H using coding region SNPs. Electrophoresis. 2006, 27 (13): 2541-2550. 10.1002/elps.200500772.

    Article  PubMed  Google Scholar 

  22. 22.

    Endicott P, Metspalu M, Stringer C, Macaulay V, Cooper A, Sanchez JJ: Multiplexed SNP typing of ancient DNA clarifies the origin of Andaman mtDNA haplogroups amongst South Asian tribal populations. PLoS One. 2006, 1: e81-10.1371/journal.pone.0000081.

    PubMed Central  Article  PubMed  Google Scholar 

  23. 23.

    Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson AR, Drouin J, Eperon IC, Nierlich DP, Roe BA, Sanger F, Schreier PH, Smith AJ, Staden R, Young IG: Sequence and organization of the human mitochondrial genome. Nature. 1981, 290 (5806): 457-465. 10.1038/290457a0.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N: Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet. 1999, 23 (2): 147-10.1038/13779.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Ermini L, Olivieri C, Rizzi E, Corti G, Bonnal R, Soares P, Luciani S, Marota I, De Bellis G, Richards MB, Rollo F: Complete mitochondrial genome sequence of the Tyrolean Iceman. Curr Biol. 2008, 18: 1687-1693. 10.1016/j.cub.2008.09.028.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Sanchez JJ, Endicott P: Developing Multiplex SNP assays with special reference to degraded DNA templates. Nature Protocols. 2006, 1 (3): 1370-1378. 10.1038/nprot.2006.247.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Brotherton P, Endicott P, Sanchez JJ, Beaumont M, Barnett R, Austin J, Cooper A: Novel high-resolution characterization of ancient DNA reveals C U-type base modification events as the sole cause of post mortem miscoding lesions. Nucleic Acids Res. 2007, 35 (17): 5717-5728. 10.1093/nar/gkm588.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  28. 28.

    Howell N, Oostra RJ, Bolhuis PA, Spruijt L, Clarke LA, Mackey DA, Preston G, Herrnstadt C: Sequence analysis of the mitochondrial genomes from Dutch pedigrees with Leber hereditary optic neuropathy. Am J Hum Genet. 2003, 72 (6): 1460-1469. 10.1086/375537.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  29. 29.

    Kivisild T, Shen P, Wall DP, Do B, Sung R, Davis K, Passarino G, Underhill PA, Scharfe C, Torroni A, Scozzari R, Modiano D, Coppa A, de Knijff P, Feldman M, Cavalli-Sforza LL, Oefner PJ: The role of selection in the evolution of human mitochondrial genomes. Genetics. 2006, 172: 373-387. 10.1534/genetics.105.043901.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  30. 30.

    Kivisild T, Metspalu M, Bandelt HJ, Richards M, Villems R: The World mtDNA Phylogeny. 2006, Springer-Verlag; Berlin, Heidelberg, Volume 18 of Nucleic Acids and Molecular Biology, 149-179.

    Google Scholar 

  31. 31.

    Richards M, Macaulay V, Torroni A, Bandelt HJ, Sykes B: Phylogeography of mitochondrial DNA in western Europe. Ann Hum Genet. 1998, 62: 241-260. 10.1046/j.1469-1809.1998.6230241.x.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Parson W, Parsons TJ, Scheithauer R, Holland MM: Population Data for 101 Austrian Caucasian mitchondrial DNA d-loop sequences: Application of mtDNA sequence analysis to a forensic case. Int J Legal Med. 1998, 111: 124-132. 10.1007/s004140050132.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Handt O, Krings M, Ward RH, Pääbo S: The retrieval of ancient human DNA sequences. Am J Hum Genet. 1996, 59 (2): 368-376.

    PubMed Central  CAS  PubMed  Google Scholar 

  34. 34.

    Handt O, Höss M, Krings M, Pääbo S: Ancient DNA: methodological challenges. Experientia. 1994, 50 (6): 524-529. 10.1007/BF01921720.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Sampietro ML, Gilbert MTB, Lao O, Caramelli D, Lari M, Bertranpetit J, Lalueza-Fox C: Tracking down Human Contamination in Ancient Human Teeth. Mol Biol Evol. 2006, 23 (9): 1801-1807. 10.1093/molbev/msl047.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Pääbo S, Irwin DM, Wilson AC: DNA damage promotes jumping between templates during enzymatic amplification. J Biol Chem. 1990, 265: 4718-4721.

    PubMed  Google Scholar 

  37. 37.

    Richards M, Sykes B, Hedges R: Authenticating DNA Extracted From Ancient Skeletal Remains. J Arch Sci. 1995, 22: 291-299. 10.1006/jasc.1995.0031.

    Article  Google Scholar 

  38. 38.

    Hummel S: Ancient DNA Typing Methods, Strategies and Applications. 2003, Springer; Berlin

    Google Scholar 

  39. 39.

    Gilbert MTP, Bandelt HJ, Hofreiter M, Barnes I: Assessing ancient DNA studies. Trends Ecol Evol. 2005, 20 (10): 541-544. 10.1016/j.tree.2005.07.005.

    Article  PubMed  Google Scholar 

  40. 40.

    Palanichamy MG, Sun C, Agrawal S, Bandelt HJ, Kong QP, Khan F, Wang CY, Chaudhuri TK, Palla V, Zhang YP: Phylogeny of mitochondrial DNA macrohaplogroup N in India, based on complete sequencing: implications for the peopling of South Asia. Am J Hum Genet. 2004, 75 (6): 966-978. 10.1086/425871.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  41. 41.

    Pfaffl MW: A new mathematical model for the relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

Download references


AC provided the laboratory facilities for the work. PE was supported by the National Environment Council of Great Britain, and Magdalen College, Oxford; JJS by a grant from Plan Nacional I+D+I (MEC. 18.08.43B.780); AC by the Wellcome and Leverhulme trusts; IP and PP were supported by the Health Care Assessorate of the Autonomous Province of Bolzano and the South Tyrolean Sparkasse Foundation.

Author information



Corresponding author

Correspondence to Phillip Endicott.

Additional information

Authors' contributions

PP and AC were responsible for the original idea of the study. PE conceived, designed and performed the experiments; analysed and interpreted the data; drafted the manuscript. JJS participated in the design of the study; designed the multiplex, and helped to draft the manuscript. IP performed analysis of modern DNA samples and edited the manuscript. PB helped with the design of the study and helped to draft the manuscript. JB helped to design and run the real-time PCR quantitation and helped to draft the manuscript. EV supplied the samples. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Cloned HVS1 sequences from the Iceman. HVR1 sequences covering diagnostic SNPs at nps 16311 and 16362, obtained from three extracts of the Iceman's DNA (BPE1, HB49, HB50). (EPS 507 KB)

Additional file 2: qPCR results for the Iceman (extract HB50). The blue diamonds are the various dilutions of the reference standard, and the red ones represent the Iceman extract. This gives the estimated starting copy number per microlitre of extract (~18 k). (EPS 196 KB)

Additional file 3: Table of mtDNA haplotypes. List of mtDNA haplotypes of all results from aDNA extracts an negative controls (NC) with putative European human mtDNA prepared in the ABC isolation facility in Oxford. rCRS means that there were no mutations observed relative to the revised Cambridge Reference Sequence. (PDF 510 KB)

Additional file 4: Table of Singleplex PCR primers used in this study. The first and fourth columns give the 5' ends of the forward and reverse primers, respectively. (PDF 676 KB)

Additional file 5: K1 multiplex PCR primers. The fourth and the sixth columns give the 5' ends of the forward and reverse primers, respectively, whilst the second column indicates the position of the site being targeted in the SBE reaction. (PDF 509 KB)

Additional file 6: K1 multiplex SBE primers. The name of the Single-Base-Extension primer indicates the site being tested and the strand orientation, whilst the last column gives the possible ancestral and derived alleles. (PDF 538 KB)

Additional file 7: Quantitative PCR primers and synthetic standards. The name of the primer indicates the position immediately 3' of the primer sequence. (PDF 495 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Endicott, P., Sanchez, J.J., Pichler, I. et al. Genotyping human ancient mtDNA control and coding region polymorphisms with a multiplexed Single-Base-Extension assay: the singular maternal history of the Tyrolean Iceman. BMC Genet 10, 29 (2009).

Download citation


  • Multiplex Polymerase Chain Reaction
  • Control Region Sequence
  • Revise Cambridge Reference Sequence
  • Code Region Polymorphism
  • aDNA Study