Typing of mitochondrial DNA coding region SNPs of forensic and anthropological interest using SNaPshot minisequencing

Multiplexed SNP typing of ancient DNA clarifies the origin of Andaman mtDNA haplogroups amongst South Asian tribal populations

The issue of errors in genetic data sets is of growing concern, particularly in population genetics where whole genome mtDNA sequence data is coming under increased scrutiny. Multiplexed PCR reactions, combined with SNP typing, are currently under-exploited in this context, but have the potential to genotype whole populations rapidly and accurately, significantly reducing the amount of errors appearing in published data sets. To show the sensitivity of this technique for screening mtDNA genomic sequence data, 20 historic samples of the enigmatic Andaman Islanders and 12 modern samples from three Indian tribal populations (Chenchu, Lambadi and Lodha) were genotyped for 20 coding region sites after provisional haplogroup assignment with control region sequences. The genotype data from the historic samples significantly revise the topologies for the Andaman M31 and M32 mtDNA lineages by rectifying conflicts in published data sets. The new Indian data extend the distribution of the M31a lineage to South Asia, challenging previous interpretations of mtDNA phylogeography. This genetic connection between the ancestors of the Andamanese and South Asian tribal groups approximately 30 kya has important implications for the debate concerning migration routes and settlement patterns of humans leaving Africa during the late Pleistocene, and indicates the need for more detailed genotyping strategies. The methodology serves as a low-cost, high-throughput model for the production and authentication of data from modern or ancient DNA, and demonstrates the value of museum collections as important records of human genetic diversity.

Extended guidelines for mtDNA typing of population data in forensic science

Mitochondrial DNA analysis has become a vital niche in forensic science as it constitutes a powerful technique for low quality and low quantity DNA samples. For the forensic field it is important to employ standardized procedures based on scientific grounds, in order to have mtDNA evidence be accepted in court. Here, we modify and extend recommendations that were spelled out previously in the absence of solid knowledge about the worldwide phylogeny. Refinement of those earlier guidelines became necessary in regard to sample selection, amplification and sequencing strategies, as well as a posteriori quality control of mtDNA profiles. The notation of sequence data should thus reflect this growing knowledge.

Coding region mitochondrial DNA SNPs: targeting East Asian and Native American haplogroups

We have developed a single PCR multiplex SNaPshot reaction that consists of 32 coding region SNPs that allows (i) increasing the discrimination power of the mitochondrial DNA (mtDNA) typing in forensic casework, and (ii) haplogroup assignments of mtDNA profiles in both human population studies (e.g. anthropological) and medical research. The selected SNPs target the East Asian phylogeny, including its Native American derived branches. We have validated this multiplex assay by genotyping a sample of East Asians (Taiwanese) and Native Americans (Argentineans). In addition to the coding SNP typing, we have sequenced the complete control region for the same samples. The genotyping results (control region plus SNaPshot profiles) are in good agreement with previous human population genetic studies (based on e.g. complete sequencing) and the known mtDNA phylogeny. We observe that the SNaPshot method is reliable, rapid, and cost effective in comparison with other techniques of multiplex SNP genotyping. We discuss the advantages of our SNP genotyping selection with respect to previous attempts, and we highlight the importance of using the known mtDNA phylogeny as a framework for SNP profile interpretation and as a tool to minimize genotyping errors.

Forensic mitochondrial coding region analysis for increased discrimination using pyrosequencing technology

Analysis of mitochondrial DNA (mtDNA) is very useful when nuclear DNA analysis fails due to degradation or insufficient amounts of DNA in forensic analysis. However, mtDNA analysis has a lower discrimination power compared to what can be obtained by nuclear DNA (nDNA) analysis, potentially resulting in multiple individuals showing identical mtDNA types in the HVI/HVII region. In this study, the increase in discrimination by analysis of mitochondrial coding regions has been evaluated for identical or similar HVI/HVII sequences. A pyrosequencing-based system for coding region analysis, comprising 17 pyrosequencing reactions performed on 15 PCR fragments, was utilised. This assay was evaluated in 135 samples, resulting in an average read length of 81 nucleotides in the pyrosequencing analysis. In the sample set, a total of 52 coding region SNPs were identified, of which 18 were singletons. In a group of 60 samples with 0 or 1 control region difference from the revised Cambridge reference sequence (rCRS), only 12 samples could not be resolved by at least two differences using the pyrosequencing assay. Thus, the use of this pyrosequencing-based coding region assay has the potential to substantially increase the discriminatory power of mtDNA analysis.

Atlantic sturgeons (Acipenser sturio, Acipenser oxyrinchus): American females successful in Europe

Recent molecular data on the maternally inherited mitochondrial (mt) DNA have challenged the traditional view that the now extinct Baltic sturgeon population belonged to the European sturgeon Acipenser sturio. Instead, there is evidence that American sea sturgeon Acipenser oxyrinchus historically immigrated into the Baltic Sea. In this study, we test the hypothesis that A. oxyrinchus introgressed into, rather than replaced, the A. sturio population in the Baltic. We established four single nucleotide polymorphisms (SNPs) in the nuclear MHC II antigen gene with a species-specific SNP pattern. Using an ancient DNA approach and two independent lines of molecular evidence (sequencing of allele-specific clones, SNaPshot), we detected both A. sturio and A. oxyrinchus alleles in the available museum material of the now extinct Baltic sturgeon population. The hybrid nature of the Baltic population was further confirmed by very high levels of heterozygosity. It had been previously postulated that the immigration of the cold-adapted A. oxyrinchus into the Baltic occurred during the Medieval Little Ice Age, when temperature likely dropped below the degree inducing spawning in A. sturio. Under this scenario, our new findings suggest that the genetic mosaic pattern in the Baltic sturgeon population (oxyrinchus mtDNA, sturio and oxyrinchus MHC alleles) is possibly caused by sex-biased introgression where spawning was largely restricted to immigrating American females, while fertilization was predominantly achieved by abundant local European males. The hybrid nature of the former Baltic sturgeon population should be taken into account in the current reintroduction measures.

Origin and expansion of haplogroup H, the dominant human mitochondrial DNA lineage in West Eurasia: the Near Eastern and Caucasian perspective

More than a third of the European pool of human mitochondrial DNA (mtDNA) is fragmented into a number of subclades of haplogroup (hg) H, the most frequent hg throughout western Eurasia. Although there has been considerable recent progress in studying mitochondrial genome variation in Europe at the complete sequence resolution, little data of comparable resolution is so far available for regions like the Caucasus and the Near and Middle East-areas where most of European genetic lineages, including hg H, have likely emerged. This gap in our knowledge causes a serious hindrance for progress in understanding the demographic prehistory of Europe and western Eurasia in general. Here we describe the phylogeography of hg H in the populations of the Near East and the Caucasus. We have analyzed 545 samples of hg H at high resolution, including 15 novel complete mtDNA sequences. As in Europe, most of the present-day Near Eastern-Caucasus area variants of hg H started to expand after the last glacial maximum (LGM) and presumably before the Holocene. Yet importantly, several hg H subclades in Near East and Southern Caucasus region coalesce to the pre-LGM period. Furthermore, irrespective of their common origin, significant differences between the distribution of hg H sub-hgs in Europe and in the Near East and South Caucasus imply limited post-LGM maternal gene flow between these regions. In a contrast, the North Caucasus mitochondrial gene pool has received an influx of hg H variants, arriving from the Ponto-Caspian/East European area.

Developing multiplexed SNP assays with special reference to degraded DNA templates

This protocol describes a single nucleotide polymorphism (SNP) genotyping strategy for highly degraded DNA, using a two-stage multiplex whereby multiple fragments are first amplified in a single exponential reaction and the products of this PCR are added to a linear single-base-extension reaction. It utilizes the analytical power of a capillary electrophoresis system to simultaneously type all the target sites. The protocol is specifically written for use with severely fragmented templates, typical of ancient DNA, and can be adapted to widely used detection platforms. The addition of the single-phase genotyping step avoids the need for the re-amplification and cloning of PCR products, while providing its own controls for the detection of contamination and allelic drop-out. This protocol can facilitate the routine analysis of up to 52 SNP markers (haploid or diploid) in 96 samples in a single day, and is recommended for the authentication of data in all areas of DNA research (population and medical genetics, forensics, ancient DNA).

East Asian mtDNA haplogroup determination in Koreans: Haplogroup-level coding region SNP analysis and subhaplogroup-level control region sequence analysis

The present study analyzed 21 coding region SNP markers and one deletion motif for the determination of East Asian mitochondrial DNA (mtDNA) haplogroups by designing three multiplex systems which apply single base extension methods. Using two multiplex systems, all 593 Korean mtDNAs were allocated into 15 haplogroups: M, D, D4, D5, G, M7, M8, M9, M10, M11, R, R9, B, A, and N9. As the D4 haplotypes occurred most frequently in Koreans, the third multiplex system was used to further define D4 subhaplogroups: D4a, D4b, D4e, D4g, D4h, and D4j. This method allowed the complementation of coding region information with control region mutation motifs and the resultant findings also suggest reliable control region mutation motifs for the assignment of East Asian mtDNA haplogroups. These three multiplex systems produce good results in degraded samples as they contain small PCR products (101-154 bp) for single base extension reactions. SNP scoring was performed in 101 old skeletal remains using these three systems to prove their utility in degraded samples. The sequence analysis of mtDNA control region with high incidence of haplogroup-specific mutations and the selective scoring of highly informative coding region SNPs using the three multiplex systems are useful tools for most applications involving East Asian mtDNA haplogroup determination and haplogroup-directed stringent quality control.

Admixture in Mexico City: implications for admixture mapping of type 2 diabetes genetic risk factors

Admixture mapping is a recently developed method for identifying genetic risk factors involved in complex traits or diseases showing prevalence differences between major continental groups. Type 2 diabetes (T2D) is at least twice as prevalent in Native American populations as in populations of European ancestry, so admixture mapping is well suited to study the genetic basis of this complex disease. We have characterized the admixture proportions in a sample of 286 unrelated T2D patients and 275 controls from Mexico City and we discuss the implications of the results for admixture mapping studies. Admixture proportions were estimated using 69 autosomal ancestry-informative markers (AIMs). Maternal and paternal contributions were estimated from geographically informative mtDNA and Y-specific polymorphisms. The average proportions of Native American, European and, West African admixture were estimated as 65, 30, and 5%, respectively. The contributions of Native American ancestors to maternal and paternal lineages were estimated as 90 and 40%, respectively. In a logistic model with higher educational status as dependent variable, the odds ratio for higher educational status associated with an increase from 0 to 1 in European admixture proportions was 9.4 (95%, credible interval 3.8-22.6). This association of socioeconomic status with individual admixture proportion shows that genetic stratification in this population is paralleled, and possibly maintained, by socioeconomic stratification. The effective number of generations back to unadmixed ancestors was 6.7 (95% CI 5.7-8.0), from which we can estimate that genome-wide admixture mapping will require typing about 1,400 evenly distributed AIMs to localize genes underlying disease risk between populations of European and Native American ancestry. Sample sizes of about 2,000 cases will be required to detect any locus that contributes an ancestry risk ratio of at least 1.5.

Dissection of mitochondrial superhaplogroup H using coding region SNPs

Haplogroup H (hg H) includes about 40-50% of the West Eurasian mitochondrial DNA (mtDNA) samples investigated so far. In order to enhance discrimination within this haplogroup we selected 45 coding region SNPs that allow to ascribe samples to the main phylogenetic branches of super hg HV (that embraces hg H) and, in particular, to H sublineages with a much finer resolution than previous studies. SNP selection was carried out using the most up-to-date available literature on population and forensic genetics and extended by means of phylogenetic analysis of complete or coding region genomes (<430) and control region sequences. A meticulous inspection of the H phylogeny led us to the observation of various but uncharacterized subclades of hg H. The selected SNPs were amplified in two PCR-multiplex reactions and subsequently targeted in three single-base extension multiplex reactions. A total of 2214 West Eurasian samples were screened for hg H specific loci 2706 and 7028, of which 859 fell in hg H and were further subjected to subhaplogroup typing. We observed 35 different subhaplogroups in total, 33 of which were found at frequencies below 5%. This assay can be used as a prescreening tool in forensic casework for rapid discrimination between divergent lineages (very effective for high-volume crime cases) or as discriminatory assay, when identical hg H haplotypes were obtained by control region sequencing.

Phylogeographic investigations: the role of trees in forensic genetics

The human mitochondrial DNA (mtDNA) genome is commonly analyzed in various disciplines, such as population, medical, and forensic genetics, but conceptual and scientific exchange between them is still limited. Here we review several aspects of the mtDNA phylogeny that are particularly--but not exclusively--of interest to the forensic community. Among the issues that arise, we emphasize the importance of integrating evolutionary concepts into the forensic routine. We also discuss topics such as mtDNA mutation-rate heterogeneity and the weight of evidence, ethnic affiliations of mtDNA profiles, and the abuse of reference databases. Finally, we show the usefulness of coding-region variation in a forensic context.

Forensic mass screening using mtDNA

At the forensic autopsy of a sexual murder victim, some trace hairs, possibly belonging to the perpetrator, were saved. Initially, the analysis of a pubic hair shaft only revealed the presence of the mitochondrial (mt) DNA haplotype profile consisting of the (CA)(6) allele and the complete hypervariable region 1 (HV1) and 2 (HV2) sequence. Later, typing of some further telogene trace hairs, which had been stored for several years, yielded a nuclear short tandem repeat (STR) profile. We used both the mtDNA haplotype and the STR profile to start a DNA mass screening project involving 2,335 male citizens of the relevant communities. MtDNA screening was carried out by using the CA repeat amplification in combination with an SNP typing procedure based on the restriction site analysis of amplified d-loop sequences. The aim of our paper is to put mass screening with mtDNA up for discussion.

Introduction of an single nucleodite polymorphism-based "Major Y-chromosome haplogroup typing kit" suitable for predicting the geographical origin of male lineages

The European Consortium "High-throughput analysis of single nucleotide polymorphisms for the forensic identification of persons--SNPforID", has performed a selection of candidate Y-chromosome single nucleotide polymorphisms (SNPs) for making inferences on the geographic origin of an unknown sample. From more than 200 SNPs compiled in the phylogenetic tree published by the Y-Chromosome Consortium, and looking at the population studies previously published, a package of 29 SNPs has been selected for the identification of major population haplogroups. A "Major Y-chromosome haplogroup typing kit" has been developed, which allows the multiplex amplification of all 29 SNPs in a single reaction. Allele genotyping was performed with a single base extension reaction (minisequencing) detected by CE. The validation of the multiplex was performed in a total of 1126 unrelated males distributed among 12 worldwide populations. The approach takes advantage of the specific geographic distribution of the Y-chromosome haplogroups and demonstrates the utility of binary polymorphisms to infer the origin of a male lineage.

A microarray system for genotyping 150 single nucleotide polymorphisms in the coding region of human mitochondrial DNA

We established a genotyping system for a panel of 150 SNPs in the coding regions of mitochondrial DNA based on multiplex tag-array minisequencing. We show the feasibility of this system for simultaneous identification of individuals and prediction of the geographical origin of the mitochondrial DNA population lineage of the sample donors by genotyping the panel of SNPs in 265 samples representing nine different populations from Africa, Europe, and Asia. Nearly 40,000 genotypes were produced in the study, with an overall genotyping success rate of 95% and accuracy close to 100%. The gene diversity value of the panel of 150 SNPs was 0.991, compared to 0.995 for sequencing 500 nucleotides of the hypervariable regions I and II of mtDNA. For 17 individuals with identical sequences in the hypervariable regions of mtDNA, our panel of SNPs increased the power of discrimination. We observed 144 haplotypes that correspond to previously determined mitochondrial "haplogroups," and they allowed prediction of the origin of the maternal population lineage of 97% of the analyzed samples.

Subtyping mtDNA haplogroup H by SNaPshot minisequencing and its application in forensic individual identification

Sequence variation of the hypervariable segments (HVS) I/II of mitochondrial DNA (mtDNA) and the haplogroup affiliation were determined in a sample of 271 Italian subjects. This analysis showed that 42% of the individuals could be ascribed to H, the most frequent haplogroup in European Caucasian populations. This fraction was then screened for specific single nucleotide polymorphisms located in the coding region to identify H subclades H1-H15. We set up two multiplex polymerase chain reactions and specific SNaPshot assays to investigate the frequency distribution of these subgroups in our population sample and to examine their usefulness in discriminating among commonly shared HVS I/II sequences. This allowed the assignment of a large portion of the mtDNAs ( approximately 70%) to specific subhaplogroups, with H1 and H5 being the most represented. About two-thirds of the individuals sharing common HVS I/II sequences were subdivided and ascribed to specific H subhaplogroups with a significant reduction of the frequencies of the most common mtDNA haplotypes. Haplogroup H subtyping could thus be extremely useful in forensic identification when many samples have to be analysed and compared, avoiding excessive time-consuming and labor-intensive sequencing analysis.

Evaluating the forensic informativeness of mtDNA haplogroup H sub-typing on a Eurasian scale

The impact of phylogeographic information on mtDNA forensics has been limited to the quality control of published sequences and databases. In this work we use the information already available on Eurasian mtDNA phylogeography to guide the choice of coding-region SNPs for haplogroup H. This sub-typing is particularly important in forensics since, even when sequencing both HVRI and HVRII, the discriminating power is low in some Eurasian populations. We show that a small set (eight) of coding-region SNPs resolves a substantial proportion of the identical haplotypes, as defined by control-region variation alone. Moreover, this SNP set, while substantially increasing the discriminating efficiency in most Eurasian populations by roughly equal amounts, discloses population-specific profiles.

High-resolution mtDNA evidence for the late-glacial resettlement of Europe from an Iberian refugium

The advent of complete mitochondrial DNA (mtDNA) sequence data has ushered in a new phase of human evolutionary studies. Even quite limited volumes of complete mtDNA sequence data can now be used to identify the critical polymorphisms that define sub-clades within an mtDNA haplogroup, providing a springboard for large-scale high-resolution screening of human mtDNAs. This strategy has in the past been applied to mtDNA haplogroup V, which represents <5% of European mtDNAs. Here we adopted a similar approach to haplogroup H, by far the most common European haplogroup, which at lower resolution displayed a rather uninformative frequency distribution within Europe. Using polymorphism information derived from the growing complete mtDNA sequence database, we sequenced 1580 base pairs of targeted coding-region segments of the mtDNA genome in 649 individuals harboring mtDNA haplogroup H from populations throughout Europe, the Caucasus, and the Near East. The enhanced genealogical resolution clearly shows that sub-clades of haplogroup H have highly distinctive geographical distributions. The patterns of frequency and diversity suggest that haplogroup H entered Europe from the Near East approximately 20,000-25,000 years ago, around the time of the Last Glacial Maximum (LGM), and some sub-clades re-expanded from an Iberian refugium when the glaciers retreated approximately 15,000 years ago. This shows that a large fraction of the maternal ancestry of modern Europeans traces back to the expansion of hunter-gatherer populations at the end of the last Ice Age.

SNPs in forensic genetics: a review on SNP typing methodologies

There is an increasing interest in single nucleotide polymorphism (SNP) typing in the forensic field, not only for the usefulness of SNPs for defining Y chromosome or mtDNA haplogroups or for analyzing the geographical origin of samples, but also for the potential applications of autosomal SNPs. The interest of forensic researchers in autosomal SNPs has been attracted due to the potential advantages in paternity testing because of the low mutation rates and specially in the analysis of degraded samples by use of short amplicons. New SNP genotyping methods, chemistries and platforms are continuously being developed and it is often difficult to be keeping up to date and to decide on the best technology options available. This review offers to the reader a state of the art of SNP genotyping technologies with the advantages and disadvantages of the different chemistries and platforms for different forensic requirements.

A DNA microarray system for forensic SNP analysis

Forensic DNA analysis is routinely performed using polymorphic short tandem repeat (STR) markers. However, for degraded or minute DNA samples, analysis of autosomal single nucleotide polymorphisms (SNPs) in short fragments might be more successful. Furthermore, sequencing of mitochondrial DNA (mtDNA) is often performed on highly degraded or scarce samples due to the high copy number of mtDNA in each cell. Due to the increasing number of complete mtDNA genome sequences available, the limited discrimination power of an mtDNA analysis, may be increased by analysis of coding region polymorphisms in addition to the non-coding variation. Since sequence analysis of the coding region would require more material than generally present in forensic samples, an alternative SNP analysis approach is required. We have developed a one-colour microarray-based SNP detection system for limited forensic materials. The method is based on minisequencing in solution prior to hybridisation to universal tag-arrays. In a first outline of a forensic chip, a combination of 12 nuclear and 21 mitochondrial SNP markers are analysed simultaneously. The mitochondrial markers on the chip are polymorphisms within the hypervariable region as well as in the coding region. Even though the number of markers in the current system is limited, it can easily be extended to yield a greater power of discrimination. When fully developed, microarray analysis provides a promising system for efficient sensitive SNP analysis of forensic samples in the future.

Phylogeny of mitochondrial DNA macrohaplogroup N in India, based on complete sequencing: implications for the peopling of South Asia

To resolve the phylogeny of the autochthonous mitochondrial DNA (mtDNA) haplogroups of India and determine the relationship between the Indian and western Eurasian mtDNA pools more precisely, a diverse subset of 75 macrohaplogroup N lineages was chosen for complete sequencing from a collection of >800 control-region sequences sampled across India. We identified five new autochthonous haplogroups (R7, R8, R30, R31, and N5) and fully characterized the autochthonous haplogroups (R5, R6, N1d, U2a, U2b, and U2c) that were previously described only by first hypervariable segment (HVS-I) sequencing and coding-region restriction-fragment-length polymorphism analysis. Our findings demonstrate that the Indian mtDNA pool, even when restricted to macrohaplogroup N, harbors at least as many deepest-branching lineages as the western Eurasian mtDNA pool. Moreover, the distribution of the earliest branches within haplogroups M, N, and R across Eurasia and Oceania provides additional evidence for a three-founder-mtDNA scenario and a single migration route out of Africa.

The molecular dissection of mtDNA haplogroup H confirms that the Franco-Cantabrian glacial refuge was a major source for the European gene pool

Complete sequencing of 62 mitochondrial DNAs (mtDNAs) belonging (or very closely related) to haplogroup H revealed that this mtDNA haplogroup--by far the most common in Europe--is subdivided into numerous subhaplogroups, with at least 15 of them (H1-H15) identifiable by characteristic mutations. All the haplogroup H mtDNAs found in 5,743 subjects from 43 populations were then screened for diagnostic markers of subhaplogroups H1 and H3. This survey showed that both subhaplogroups display frequency peaks, centered in Iberia and surrounding areas, with distributions declining toward the northeast and southeast--a pattern extremely similar to that previously reported for mtDNA haplogroup V. Furthermore, the coalescence ages of H1 and H3 (~11,000 years) are close to that previously reported for V. These findings have major implications for the origin of Europeans, since they attest that the Franco-Cantabrian refuge area was indeed the source of late-glacial expansions of hunter-gatherers that repopulated much of Central and Northern Europe from ~15,000 years ago. This has also some implications for disease studies. For instance, the high occurrence of H1 and H3 in Iberia led us to re-evaluate the haplogroup distribution in 50 Spanish families affected by nonsyndromic sensorineural deafness due to the A1555G mutation. The survey revealed that the previously reported excess of H among these families is caused entirely by H3 and is due to a major, probably nonrecent, founder event.

Disuniting uniformity: a pied cladistic canvas of mtDNA haplogroup H in Eurasia

It has been often stated that the overall pattern of human maternal lineages in Europe is largely uniform. Yet this uniformity may also result from an insufficient depth and width of the phylogenetic analysis, in particular of the predominant western Eurasian haplogroup (Hg) H that comprises nearly a half of the European mitochondrial DNA (mtDNA) pool. Making use of the coding sequence information from 267 mtDNA Hg H sequences, we have analyzed 830 mtDNA genomes, from 11 European, Near and Middle Eastern, Central Asian, and Altaian populations. In addition to the seven previously specified subhaplogroups, we define fifteen novel subclades of Hg H present in the extant human populations of western Eurasia. The refinement of the phylogenetic resolution has allowed us to resolve a large number of homoplasies in phylogenetic trees of Hg H based on the first hypervariable segment (HVS-I) of mtDNA. As many as 50 out of 125 polymorphic positions in HVS-I were found to be mutated in more than one subcluster of Hg H. The phylogeographic analysis revealed that sub-Hgs H1*, H1b, H1f, H2a, H3, H6a, H6b, and H8 demonstrate distinct phylogeographic patterns. The monophyletic subhaplogroups of Hg H provide means for further progress in the understanding of the (pre)historic movements of women in Eurasia and for the understanding of the present-day genetic diversity of western Eurasians in general.