PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2010; 5(5): e10218.
Published online 2010 May 3. doi:  10.1371/journal.pone.0010218
PMCID: PMC2862705

A Reduced Number of mtSNPs Saturates Mitochondrial DNA Haplotype Diversity of Worldwide Population Groups

Vincent Macaulay, Editor

Abstract

Background

The high levels of variation characterising the mitochondrial DNA (mtDNA) molecule are due ultimately to its high average mutation rate; moreover, mtDNA variation is deeply structured in different populations and ethnic groups. There is growing interest in selecting a reduced number of mtDNA single nucleotide polymorphisms (mtSNPs) that account for the maximum level of discrimination power in a given population. Applications of the selected mtSNP panel range from anthropologic and medical studies to forensic genetic casework.

Methodology/Principal Findings

This study proposes a new simulation-based method that explores the ability of different mtSNP panels to yield the maximum levels of discrimination power. The method explores subsets of mtSNPs of different sizes randomly chosen from a preselected panel of mtSNPs based on frequency. More than 2,000 complete genomes representing three main continental human population groups (Africa, Europe, and Asia) and two admixed populations (“African-Americans” and “Hispanics”) were collected from GenBank and the literature, and were used as training sets. Haplotype diversity was measured for each combination of mtSNP and compared with existing mtSNP panels available in the literature. The data indicates that only a reduced number of mtSNPs ranging from six to 22 are needed to account for 95% of the maximum haplotype diversity of a given population sample. However, only a small proportion of the best mtSNPs are shared between populations, indicating that there is not a perfect set of “universal” mtSNPs suitable for all population contexts. The discrimination power provided by these mtSNPs is much higher than the power of the mtSNP panels proposed in the literature to date. Some mtSNP combinations also yield high diversity values in admixed populations.

Conclusions/Significance

The proposed computational approach for exploring combinations of mtSNPs that optimise the discrimination power of a given set of mtSNPs is more efficient than previous empirical approaches. In contrast to precedent findings, the results seem to indicate that only few mtSNPs are needed to reach high levels of discrimination power in a population, independently of its ancestral background.

Introduction

Variations in human mtDNA molecules have been deeply investigated in several fields of biomedical research such as forensic genetics, molecular anthropology, and disease studies [1], [2], [3], [4], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19]. There are 100 to 10,000 copies of the mtDNA genome per cell, and each of them consists of circular molecules of about 16,569 base pairs (bps). The mtDNA is inherited exclusively from the mother. Each mtDNA genome can be divided into two main parts: the control and the coding region. The control region occupies about 1,200 bp of the molecule and contains, among other, regulatory elements related to the replication of mtDNA or gene expression. It is usually divided into two segments characterised by their high (average) mutation rate, namely the first and second hypervariable segments (HVS-I/II). The coding region encodes 37 densely packed genes (13 for proteins, 22 for transfer RNA [tRNA], and two for subunits of ribosomal RNA [rRNA]) that are needed to maintain the correct function of the mitochondrion.

The study of mtDNA variability has been approached using different methodologies. Analysis of restriction fragment length polymorphism (RFLP) sites and screening procedures such as heteroduplex analysis (HD) and single strand conformation polymorphisms (SSCP) have been extensively used in the past and are still used by many laboratories [2], [20], [21]. Since the mid-1990s, sequencing the HVS-I (and sporadically the HVS-II), coupled with the analysis of selected RFLP sites, has been the most common strategy for analysing mtDNA variation. Nowadays, there are more than 130,000 HVS-I mtDNAs from different human population groups reported in the literature.

The interest in analysing complete mtDNA genomes is growing as indicated by the more than 6,700 complete genomes available in the literature and in GenBank, most of them reported in the past five years. Sequencing complete genomes is now benefiting from improvements to the sequencing chemistry and the higher sophistication of automatic sequencers. However, this analysis is technically complex and costly and, therefore, unfeasible not only for most biomedical applications at high-throughput scales, but also for those applications depending on low quality or small amounts of DNA (e.g., forensic samples). This is the main reason why the majority of laboratories just target the control region (usually the HVS-I), complemented by analysing selected coding region SNPs. The traditional screening approaches (RFLP, SSCP) for genotyping mtSNPs are now being replaced by lower cost mini-sequencing techniques that allow multiplexing several polymorphisms in single reactions [22], [23], [24], [25]. In particular, forensic geneticists are especially interested in developing mini-sequencing assays that interrogate a reduce numbers of mtSNPs per reaction (e.g., multiplexes of five to 20 mtSNPs), because the technique can perform well with degraded or low copy number samples [26]; many evidentiary samples only allow a single PCR reaction, and high-throughput mtSNP techniques are unsuitable for sub-optimal samples [27]. Ideally, targeted mtSNPs should retain the maximum level of discrimination power in a single mtDNA test and, therefore, a careful selection of mtSNPs for mini-sequencing assays is a key step in the process. Generally, this selection is based on phylogenetic criteria, where mtSNPs are chosen from the phylogeny to represent the main branches of the phylogenetic tree (which define haplogroups); however, this strategy does not necessarily optimise the discrimination power of a particular set of mtSNPs. Alternatively, mtSNPs can be selected according to their mutation rate; mtSNPs with the highest mutation rates tend to yield the maximum diversity values. However, current positional mutation rates are under suspicion because the methods employed to compute them could be flawed [28], [29], and only recently a new proposal of site-specific mutation rates has been published aimed at overcoming the problems of past approaches [29]. It is impossible to decide a priori which one of these two approaches is more efficient for maximizing discrimination power. It is possible that a combination of both rationales could perform better. On the other hand, a particular combination of mtSNPs could yield good results in a specific population context, but might be unsuitable in a different population group (say Europeans versus Native Americans). Ideally, we could also envisage selecting a universal mtSNP panel that could generate reasonable discrimination power independently of the population group considered. Moreover, this panel would be particularly useful when dealing with highly admixed populations (e.g., the US, South American admixed populations, highly cosmopolitan cities).

With the latter applications in mind, we aimed to explore the optimum combinations of mtSNPs needed to maximize the diversity values of a given population group, taking into account the premise that multiplexing techniques (excluding high-throughput platforms) only allow genotyping a moderate number of mtSNPs (usually a maximum of 20 to 30). The method employed here is based on an algorithm that allows the exploration of the full set of combinations arising from a given set of known mtSNPs. We then evaluate the combinations yielding the highest values of diversity and the best candidate mtSNPs within these combinations. Various biomedical applications will also be discussed.

Material and Methods

Complete genome database

The database used for the simulation experiments was built based on the following criteria:

  1. Selecting complete genomes from the literature and GenBank capable of being used as proxies of human population samples representing main continental regions. This criterion filters out those available complete genomes that have been analysed based on phylogenetic criteria or, in particular, patients in (mtDNA) disease studies [1], [30], [31], [32], [33], [34], [35], [36];
  2. The main population groups should be represented by at least 300 complete genomes such that most of the considered ‘speedy’ mutations (sensu [37]) can be polymorphic with a minimum allele frequency (MAF) >5% in the whole database; and
  3. The compiled database should represent at least three main continental groups, namely, Africans, Asians, and Europeans.

According to these criteria, we collected the following datasets: (i) N = 309 from [5], representing the African subset; (ii) N = 672 Japanese to represent the Asian subset [38]; and (iii) N = 241 from [24] and N = 192 [39] representing the European subset. In addition, N = 326 individuals belonging to the American haplogroups A2, B2, C1, D1, and X2a (see [3], [10] and references therein) were also collected for representing a Native American subset. Apart from the mentioned relatively homogenous population groups, two admixed population samples from the US, ‘African-American’ and ‘Hispanic’ datasets (N = 140 and N = 125 respectively) from [40], were also used in the simulation experiments.

We are aware of potential sequencing errors affecting some of these datasets and have tried to disregard suspicious complete genomes or use corrected versions of the reported datasets (not necessarily those available in GenBank) [29], [41]. For instance, some of the Tanaka's complete genomes were corrected in Kong et al. [42] but the flawed versions of these genomes are still in GenBank [41].

As usual, mtSNPs are referred using the revised Cambridge Reference Sequence or rCRS [43].

Panel of SNPs

The approach followed in this study is based on exploring all (or a subset of all) possible combinations that arise from combining a given set of n mtSNPs candidates taken m at a time (in what follows, m value). Given the fact that the number of mtSNP variants considered in this study is above 3,200, it is computationally impossible to explore the entire universe of combinations arising from such a large number of variants (for instance, the number of possible combinations C(100,20) is >1022). To overcome this problem, we selected the 394 mtSNPs from the whole dataset fulfilling a MAF >0.05. By definition this criterion eliminates rare mtSNPs (MAF≤0.05) under the premise that these SNPs cannot substantially contribute to increase diversity levels of populations. These variants are mainly ‘private’ to single genomes (singleton mutations usually located at the tips of the mtDNA phylogenies); therefore, these variants cannot be extrapolated to independent or larger samples due to ascertainment bias. On the contrary (see above), the selected mtSNPs (and in particular those that better contribute to increase the discrimination power) are polymorphic in different human populations; either (i) because they mutated before the divergence of major population groups, or (ii) because they have a high mutation rate (mutational hotspots). Variants that are known to be problematic from a genotyping point of view were disregarded (16182C, 16183C, 16193+C, variants around 310, length variation around positions 523–524, etc.).

Programming

For the panel of mtSNPs considered and all possible m values, we computed two diversity indices in every population dataset: (a) the ‘normalised’ number of haplotypes, defined as H = h/N, where h represents the number of different haplotypes in the dataset defined by a given m and N is the dataset sample size; and (b) the haplotype diversity, defined as HD = 1 − An external file that holds a picture, illustration, etc.
Object name is pone.0010218.e001.jpg, where pi is the haplotype frequency of the i haplotype resulting from each m combination of mtSNPs. We are aware to the fact that H does not scale out the dependence on sample size [44], [45], [46], and therefore the results on H cannot be extrapolated across samples of different sizes.

The main drawback of this study is the extremely high computational cost needed to explore the entire universe of all possible mtSNP combinations. We run a parallelisable script on a shared memory node system with the SMP NUMA architecture and a cluster containing over 2500 processors named Finis Terrae, which is available at the Supercomputing Centre of Galicia (CESGA; http://www.cesga.es/). Considering that studying each of the m values is completely independent from the rest, we dispatched each individual analysis run among the available cores of the cluster, reducing the overall running time to the largest m value analysis as if performed alone.

We ran up to 10,000 iterations for each m. We recorded the values of H and HD and derived their mean and standard deviation values. We also recorded the maximum values of diversity obtained for each m among the full number of iterations.

For every dataset, we also computed the maximum diversity values (HMAX and HDMAX) that could be obtained considering the complete genome information.

Results

Considerations about computational limitations

The large amount of mtSNPs considered in the different mtSNP panels tested made it computationally unfeasible to explore all possible combinations of n mtSNPs taken m at a time. We, therefore, developed an algorithm that creates subsets of m mtSNPs by sampling without replacement and equal probability each studied panel, while recording the best mtSNP combinations for every m value (the ones giving the maximum diversity values of H and HD). Other sampling strategies were used, such as those based on a stepwise mtSNP selection (find the best single mtSNP and then the best one to add it, and so on) but the one presented here was the most efficient in terms of maximizing the values of H and HD (data not shown). However, we first assessed the impact of analysing only a given number of combinations (iterations) of mtSNPs on the estimation of the parameters of interest. We observed that the trend in H or HD values with different m values was reasonably stable when exploring at least 10,000 mtSNP combinations (Figures S1, S2, S3, S4, S5). Based on this simulation, we decided to run 10,000 iterations for every population scenario considered in this study.

We are aware that this approach will probably prevent capture of the best candidate combination of mtSNPs (the one providing the highest real value of H and HD), but this is not the main goal of the simulations. The simulation approach could be considered a success if the best combination of mtSNPs obtained from it performed much better than traditional procedures. On the other hand, note that the standard deviation values of H and HD for the different m values (Figure 1) were very low, indirectly indicating not only that these values do not fluctuate significantly on different iterations, but also that the best values should be close to the ones obtained using 10,000 iterations.

Figure 1
Values of H (left panel) and HD (right panel) computed in different population groups.

Performance of different mtSNP combinations on diversity values

The shapes of the H and HD curves in Figure 1 respond to the same expected phenomenon: increased values of m are accompanied by progressive increments in the values of H and HD. However, the m values needed to reach a plateau are different for H and HD. While H grows slowly by increasing the number of mtSNPs analysed, HD reaches its maximum quickly at lower mtSNP m values. Thus, somehow unexpectedly, the HD curves (Figures 1) indicate that only a small number of mtSNPs is required to saturate the diversity values to 95% of the HDMAX independently of the population database employed (Figures 1 and Table 1). In other words, the addition of new mtSNPs beyond a certain point does not apparently contribute to an increase in discrimination power.

Table 1
Diversity values (H and HD) for the SNPs considered in different articles published in the literature and comparison with those obtained in the present study.

The mathematical difference between HMAX and HDMAX and the mean values of H and HD are more pronounced for intermediate m values (Figures 1) because these m values admit more different combinations of mtSNPs, thereby yielding a wider range of H and HD values. Concordant with this observation is the shape of the standard deviations, indicating that the maximum fluctuations of H and HD values also occur around these intermediate m values.

Evaluation of the best candidate mtSNPs

All the mtSNPs were ranked in a list according to the number of times each of them appears in a combination that maximizes the HD values; the first position being assigned to the one that appears most often (Table S1).

To some extent, the best mtSNPs in the different population groups considered in this study overlap (Table S1 and Table S2). For instance, transition T16519C appears to be the best mtSNP in virtually all the different population groups. Other well-known mutational hotspots, T152C, T195C, T16189C, and T16311C, occupy different positions in the ranking (which in part mirrors the stochastic nature of the simulation approach used in the present study), but all are listed among the top 15 top mtSNPs in the three main continental groups when looking at the HD values (Table 2 and S1).

Table 2
Excerpt of the data in Table S1 showing the top five mtSNP that are shared between the top 15 mtSNPs in the three main continental groups.

The scores in Table S1 are useful for those analysts (e.g. forensic geneticists) wishing to empirically design a test panel of mtSNPs in a given population context. Thus, for instance, from Figures 1 (and Table 1), we know that nine mtSNPs sufficiently account for 95% of the haplotype diversity in a given African sample. From Table S1, it can be observed that the top nine polymorphisms in the African sample are for the values on HD (sorted by their position in the ranking): T16519C, T152C, T16189C, A189G, T16093C, G16129A, T195C, T16311C, and G143A.

Diversity accounted by existing mtSNP panels and the simulation-based approach

A number of different mtSNP panels have been proposed in the literature to date, and these provide a suitable framework for evaluating the efficiency of the simulation-based approach used in this study.

Table 1 summarises the values of H and HD for the different mtSNP panels proposed in the literature and those used in this study, all of them evaluated using the same collections of complete genomes (see M&M). Some of the panels reported in the literature were designed for specific population groups, mainly Europeans and Asians, so these panels behave worse in population samples with different genetic backgrounds. For instance, the panel of 45 mtSNPs proposed by Brandstätter et al. [25] was designed to increase the discrimination power within the typically European haplogroup H, but the discrimination power of these mtSNPs in Asians or Native Americans is much lower (Table 1). The opposite example is the mtSNP panel provided in Álvarez-Iglesias et al. [47].

The results summarised in Table 1 clearly indicate that almost all the panels proposed in the literature yield lower values of H and HD than the mtSNPs combinations inferred from the present study (with the caveat that H is not appropriate for inter-population comparisons). In addition, the number of mtSNPs needed to reach 95% of the HMAX and HDMAX are very low for most of the population groups (ranging from 10 to 22); while all the panels available in the literature yield values of diversity well below 95% (Table 1).

Evolutionary nature of the most discriminating mtSNPs

Not surprisingly, the top five mtSNPs that better contribute to increase the diversity values in the different population groups are non-synonymous or are located in non-coding or un-translated regions of the mtDNA molecule (Table 2). In other words, these mtSNPs are almost polymorphic in different human population groups (universal) because they do not seem to be as subjected to the effect of stabilizing or purifying natural selection as coding region variants. These top mtSNPs are in fact located in the control region (non-coding), are transitions, and roughly match the top mutational hotspots reported by [29] (Table 2).

Discussion

Several fields of research (population, medical, and forensic genetic) are interested in the analysis of mtSNPs for several genetic applications. Some of these applications rest on the ability of a selected (low) number of mtSNPs providing a high discrimination power when analysing human population or casework evidentiary samples. We have used a simulation-based approach to evaluate the discrimination power of mtSNPs in different population contexts, including samples representing the main continental regions and admixed population groups.

Our simulated approach indicates that no more than a dozen mtSNPs is sufficient to account for ~95% of the maximum level of HD diversity for almost all population groups. However, admixed populations, such as ‘African-Americans’, need to double this amount (~22 mtSNPs) to reach similar values of diversity. The top mtSNP variants are mutational hotspots mainly located in the control region. In addition, only a small proportion of the best mtSNPs are shared between population groups, indicating thatthere is not a perfect set of ‘universal’ mtSNPs suitable for all population contexts.

Today, there are two well-known strategies commonly used for selecting mtSNPs aimed to account for the highest levels of diversity in population groups, namely, the phylogenetic-based approach and mtSNP selection based on mutational hotspots. The results of the present study seem to indicate however that these traditional strategies perform worse than the simulation-based approach developed in the present study, and that in reality, it is a combination of highly mutable mtSNPs and haplogroup diagnostic sites that optimized the ability of a given mtSNP panels to account for the highest levels of diversity.

The strategy proposed in this study is particularly relevant for forensic studies, where small panels of mtSNPs are frequently demanded in routine casework, but we can also foresee other biomedical applications. For instance, mtSNP panels could be used to evaluate mtDNA instability in studies on tumours [48], [49], [50], where patients belong to different population groups.

Some caveats should be added about the potential bias arising from the limited number of complete genomes considered in this study. Some potential ascertainment bias could arise when computing HD using a limited number of complete genomes. The reproducibility of these results will only be possible with the availability of independent complete genome datasets. A bootstrap strategy could be used instead, but the high computational demands of this procedure makes it unfeasible. However, it can tentatively be said that there are few considerations that allow us to predict reproducible results when applied to independent complete genome datasets. These considerations are that (a) the best mtSNPs have by definition an MAF above 5%, (b) various of the best mtSNPs overlap in different population groups, and (c) most of the best mtSNPs have a high mutation rate. Therefore, these mtSNPs appear in different parts of the worldwide phylogeny and are not restricted to any particular population or ethnic group.

This study has demonstrated that future proposals for mtSNP panels aimed at obtaining a high discrimination power could be considered in the light of the simulation approach proposed here. The phylogenetic approach, although essential for most mtDNA studies (e.g. [19], [51], [52], [53]), is probably not the best tool for predicting the discrimination power of a particular set of a mtSNP panel, but can still be useful for understanding the biological nature of a selected panel of mtSNPs and assist in its design. The mtSNP panels proposed in the literature do not perform as well as those suggested by this study.

Supporting Information

Figure S1

Effect of size iteration (number of mtDNA combinations explored from the full universe of possible combinations) for the estimation of H and HD in the African dataset. Only the mtSNPs overlapping in all the population datasets were used.

(0.31 MB TIF)

Figure S2

Effect of size iteration for the estimation of H and HD in the ‘African-American’ dataset. Only the mtSNPs overlapping in all the population datasets were used.

(0.25 MB TIF)

Figure S3

Effect of size iteration for the estimation of H and HD in the Asian dataset. Only the mtSNPs overlapping in all the population datasets were used.

(0.29 MB TIF)

Figure S4

Effect of size iteration for the estimation of H and HD in the European dataset. Only the mtSNPs overlapping in all the population datasets were used.

(0.27 MB TIF)

Figure S5

Effect of size iteration for the estimation of H and HD in the ‘Hispanic’ dataset. Only the mtSNPs overlapping in all the population datasets were used.

(0.24 MB TIF)

Table S1

Scores of all the mtSNPs in the different population samples (AFR  = Africans, AFR-AM  = “African-Americans”, ASI  = Asians, EUR  = Europeans, and HIS  = “Hispanics”; ALL  = all the complete genomes considered as a single group) based on HD. First, the number of times a particular mtSNP shows up when computing the maximum values of HD in every iteration and m values is recorded. The mtSNPs are sorted from those that appear more times in the different iterations to those that never appear. The final score is given according to their relative position in this ranking; from 1 (received by the best mtSNP) to n (number of mtSNP in each panel). The columns indicate, in this order, the position of the mtSNP according to the rCRS [48], the mutational change (all are transitions unless a suffix indicates a transversion or an indel specified), the population sample set (as indicated above), the rCRS variant, the nature of the mutational change (transition, transvertion, or indel), the MapLocus, a shorthand of the locus, the locus description, the coding or non-coding condition of the mutational change, the changed-codon, the amino acid change, the relative position of the mtSNP within genes, the relative position of the mtSNP within codons, and their synonymous or non-synonymous condition.

(0.12 MB XLS)

Table S2

Scores of all the mtSNPs in the different population samples based on H. More details in legend of Table S1.

(0.12 MB XLS)

Acknowledgments

We would like to acknowledge the CESGA for their supercomputing availability and support and Yong-Gang Yao and an anonymous reviewer for critically reading and commenting on the manuscript.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: Two grants from Fundacion de Investigacion Medica Mutua Madrile a (2006/CL370 and 2008/CL444), a grant Grupos Emerxentes from Xunta de Galicia (2008/XA122) and “Fondos Feder”, and a grant from the Ministerio de Ciencia e Innovacion (SAF2008-02971) awarded to A. Salas supported this project. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Macaulay V, Hill C, Achilli A, Rengo C, Clarke D, et al. Single, rapid coastal settlement of Asia revealed by analysis of complete mitochondrial genomes. Science. 2005;308:1034–1036. [PubMed]
2. Macaulay V, Richards M, Hickey E, Vega E, Cruciani F, et al. The emerging tree of West Eurasian mtDNAs: a synthesis of control-region sequences and RFLPs. Am J Hum Genet. 1999;64:232–249. [PubMed]
3. Perego UA, Achilli A, Angerhofer N, Accetturo M, Pala M, et al. Distinctive Paleo-Indian migration routes from Beringia marked by two rare mtDNA haplogroups. Curr Biol. 2009;19:1–8. [PubMed]
4. Herrnstadt C, Elson JL, Fahy E, Preston G, Turnbull DM, et al. Reduced-median-network analysis of complete mitochondrial DNA coding-region sequences from the major African, Asian, and European haplogroups. Am J Hum Genet. 2002;70:1152–1171. [PubMed]
5. Behar DM, Villems R, Soodyall H, Blue-Smith J, Pereira L, et al. The dawn of human matrilineal diversity. Am J Hum Genet. 2008;82:1130–1140. [PubMed]
6. Torroni A, Achilli A, Macaulay V, Richards M, Bandelt H-J. Harvesting the fruit of the human mtDNA tree. Trends Genet. 2006;22:339–345. [PubMed]
7. Salas A, Richards M, De la Fé T, Lareu MV, Sobrino B, et al. The making of the African mtDNA landscape. Am J Hum Genet. 2002;71:1082–1111. [PubMed]
8. Richards M, Macaulay V, Hickey E, Vega E, Sykes B, et al. Tracing European founder lineages in the Near Eastern mtDNA pool. Am J Hum Genet. 2000;67:1251–1276. [PubMed]
9. Baudouin SV, Saunders D, Tiangyou W, Elson JL, Poynter J, et al. Mitochondrial DNA and survival after sepsis: a prospective study. Lancet. 2005;366:2118–2121. [PubMed]
10. Achilli A, Perego UA, Bravi CM, Coble MD, Kong Q-P, et al. The phylogeny of the four pan-American MtDNA haplogroups: implications for evolutionary and disease studies. PLoS ONE. 2008;3:e1764. [PMC free article] [PubMed]
11. Kivisild T, Reidla M, Metspalu E, Rosa A, Brehm A, et al. Ethiopian mitochondrial DNA heritage: tracking gene flow across and around the gate of tears. Am J Hum Genet. 2004;75:752–770. [PubMed]
12. Kivisild T, Rootsi S, Metspalu M, Mastana S, Kaldma K, et al. The genetic heritage of the earliest settlers persists both in Indian tribal and caste populations. Am J Hum Genet. 2003;72:313–332. [PubMed]
13. Torroni A, Bandelt H-J, D'Urbano L, Lahermo P, Moral P, et al. mtDNA analysis reveals a major late Paleolithic population expansion from southwestern to northeastern Europe. Am J Hum Genet. 1998;62:1137–1152. [PubMed]
14. Torroni A, Bandelt H-J, Macaulay V, Richards M, Cruciani F, et al. A signal, from human mtDNA, of postglacial recolonization in Europe. Am J Hum Genet. 2001;69:844–852. [PubMed]
15. Chinnery PF, Howell N, Andrews RM, Turnbull DM. Clinical mitochondrial genetics. J Med Genet. 1999;36:425–436. [PMC free article] [PubMed]
16. Chinnery PF, Howell N, Andrews RM, Turnbull DM. Mitochondrial DNA analysis: polymorphisms and pathogenicity. J Med Genet. 1999;36:505–510. [PMC free article] [PubMed]
17. Chinnery PF, Johnson MA, Wardell TM, Singh-Kler R, Hayes C, et al. The epidemiology of pathogenic mitochondrial DNA mutations. Ann Neurol. 2000;48:188–193. [PubMed]
18. Carelli V, Giordano C, d'Amati G. Pathogenic expression of homoplasmic mtDNA mutations needs a complex nuclear-mitochondrial interaction. Trends Genet. 2003;19:257–262. [PubMed]
19. Salas A, Bandelt H-J, Macaulay V, Richards MB. Phylogeographic investigations: the role of trees in forensic genetics. Forensic Sci Int. 2007;168:1–13. [PubMed]
20. Barros F, Lareu MV, Salas A, Carracedo A. Rapid and enhanced detection of mitochondrial DNA variation using single-strand conformation analysis of superposed restriction enzyme fragments from polymerase chain reaction-amplified products. Electrophoresis. 1997;18:52–54. [PubMed]
21. Salas A, Rasmussen EM, Lareu MV, Morling N, Carracedo Á. Fluorescent SSCP of overlapping fragments (FSSCP-OF): a highly sensitive method for the screening of mitochondrial DNA variation. Forensic Sci Int. 2001;124:97–103. [PubMed]
22. Álvarez-Iglesias V, Barros F, Carracedo Á, Salas A. Minisequencing mitochondrial DNA pathogenic mutations. BMC Med Genet. 2008;9:26. [PMC free article] [PubMed]
23. Quintáns B, Álvarez-Iglesias V, Salas A, Phillips C, Lareu MV, et al. Typing of mitochondrial DNA coding region SNPs of forensic and anthropological interest using SNaPshot minisequencing. Forensic Sci Int. 2004;140:251–257. [PubMed]
24. Coble MD, Just RS, O'Callaghan JE, Letmanyi IH, Peterson CT, et al. 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. [PubMed]
25. Brandstätter A, Salas A, Niederstätter H, Gassner C, Carracedo Á, et al. Dissection of mitochondrial superhaplogroup H using coding region SNPs. Electrophoresis. 2006;27:2541–2550. [PubMed]
26. Mosquera-Miguel A, Álvarez-Iglesias V, Lareu MV, Carracedo Á, Salas A. Testing the performance of mtSNP minisequencing in forensic samples. Forensic Sci Int Genet. 2009;3:261–264. [PubMed]
27. Cerezo M, Černý V, Carracedo Á, Salas A. Applications of MALDI-TOF MS to large-scale human mtDNA population-based studies. Electrophoresis. 2009;30:3665–3673. [PubMed]
28. Bandelt H-J, Kong Q-P, Richards M, Macaulay V. Estimation of mutation rates and coalescence times: some caveats. In: Bandelt H-J, Richards M, Macaulay V, editors. Human mitochondrial DNA and the evolution of Homo sapiens. Berlin: Springer-Verlag; 2006. pp. 47–90.
29. Soares P, Ermini L, Thomson N, Mormina M, Rito T, et al. Correcting for purifying selection: an improved human mitochondrial molecular clock. Am J Hum Genet. 2009;84:740–759. [PubMed]
30. Soares P, Trejaut JA, Loo JH, Hill C, Mormina M, et al. Climate change and postglacial human dispersals in southeast Asia. Mol Biol Evol. 2008;25:1209–1218. [PubMed]
31. Achilli A, Rengo C, Battaglia V, Pala M, Olivieri A, et al. Saami and Berbers–an unexpected mitochondrial DNA link. Am J Hum Genet. 2005;76:883–886. [PubMed]
32. Olivieri A, Achilli A, Pala M, Battaglia V, Fornarino S, et al. The mtDNA legacy of the Levantine early Upper Palaeolithic in Africa. Science. 2006;314:1767–1770. [PubMed]
33. Kivisild T, Tolk H-V, Parik J, Wang Y, Papiha SS, et al. The emerging limbs and twigs of the East Asian mtDNA tree. Mol Biol Evol. 2002;19:1737–1751. [PubMed]
34. Brisighelli F, Capelli C, Álvarez-Iglesias V, Onofri V, Paoli G, et al. The Etruscan timeline: a recent Anatolian connection. Eur J Hum Genet. 2009;17:693–696. [PMC free article] [PubMed]
35. Álvarez-Iglesias V, Mosquera-Miguel A, Cerezo M, Quintáns B, Zarrabeitia MT, et al. New population and phylogenetic features of the internal variation within mitochondrial DNA macro-haplogroup R0. PLoS ONE. 2009;4:e5112. [PMC free article] [PubMed]
36. Quintana-Murci L, Quach H, Harmant C, Luca F, Massonnet B, et al. Maternal traces of deep common ancestry and asymmetric gene flow between Pygmy hunter-gatherers and Bantu-speaking farmers. Proc Natl Acad Sci U S A. 2008;105:1596–1601. [PubMed]
37. Bandelt H-J, Quintana-Murci L, Salas A, Macaulay V. The fingerprint of phantom mutations in mitochondrial DNA data. Am J Hum Genet. 2002;71:1150–1160. [PubMed]
38. Tanaka M, Cabrera VM, González AM, Larruga JM, Takeyasu T, et al. Mitochondrial genome variation in eastern Asia and the peopling of Japan. Genome Res. 2004;14:1832–1850. [PubMed]
39. Finnilä S, Lehtonen MS, Majamaa K. Phylogenetic network for European mtDNA. Am J Hum Genet. 2001;68:1475–1484. [PubMed]
40. Just RS, Diegoli TM, Saunier JL, Irwin JA, Parsons TJ. Complete mitochondrial genome sequences for 265 African American and U.S. “Hispanic” individuals. Forensic Sci Int Genet. 2008;2:e45–48. [PubMed]
41. Yao Y-G, Salas A, Logan I, Bandelt H-J. mtDNA data-mining in GenBank needs surveying. Am J Hum Genet in press. 2009. [PubMed]
42. Kong Q-P, Salas A, Sun C, Fuku N, Tanaka M, et al. Distilling artificial recombinants from large sets of complete mtDNA genomes. PLoS ONE. 2008;3:e3016. [PMC free article] [PubMed]
43. Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, et al. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet. 1999;23:147. [PubMed]
44. Ewens WJ. The sampling theory of selectively neutral alleles. Theoretical Population Biology. 1972;3:87–112. [PubMed]
45. Egeland T, Salas A. Statistical evaluation of haploid genetic evidence. TOForensicSJ. 2008;1:4–11.
46. Egeland T, Salas A. Estimating haplotype frequency and coverage of databases. PLoS ONE. 2008;3:e3988. [PMC free article] [PubMed]
47. Álvarez-Iglesias V, Jaime JC, Carracedo Á, Salas A. Coding region mitochondrial DNA SNPs: targeting East Asian and Native American haplogroups. Forensic Sci Int Genet. 2007;1:44–55. [PubMed]
48. Salas A, Yao Y-G, Macaulay V, Vega A, Carracedo Á, et al. A critical reassessment of the role of mitochondria in tumorigenesis. PLoS Med. 2005;2:e296. [PMC free article] [PubMed]
49. Vega A, Salas A, Gamborino E, Sobrido MJ, Macaulay V, et al. mtDNA mutations in tumors of the central nervous system reflect the neutral evolution of mtDNA in populations. Oncogene. 2004;23:1314–1320. [PubMed]
50. Cerezo M, Bandelt H-J, Martin-Guerrero I, Ardanaz M, Vega A, et al. High mitochondrial DNA stability in B-cell chronic lymphocytic leukemia. PLoS One. 2009;4:e7902. [PMC free article] [PubMed]
51. Salas A, Carracedo Á, Macaulay V, Richards M, Bandelt H-J. A practical guide to mitochondrial DNA error prevention in clinical, forensic, and population genetics. Biochem Biophys Res Commun. 2005;335:891–899. [PubMed]
52. Yao Y-G, Bravi CM, Bandelt H-J. A call for mtDNA data quality control in forensic science. Forensic Sci Int. 2004;141:1–6. [PubMed]
53. Bandelt H-J, Olivieri A, Bravi C, Yao Y-G, Torroni A, et al. ‘Distorted’ mitochondrial DNA sequences in schizophrenic patients. Eur J Hum Genet. 2007;15:400–402; author reply 402-404. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science