Full mtGenome reference data: Development and characterization of 588 forensic-quality haplotypes representing three U.S. populations

Though investigations into the use of massively parallel sequencing technologies for the generation of complete mitochondrial genome (mtGenome) profiles from difficult forensic specimens are well underway in multiple laboratories, the high quality population reference data necessary to support full mtGenome typing in the forensic context are lacking. To address this deficiency, we have developed 588 complete mtGenome haplotypes, spanning three U.S. population groups (African American, Caucasian and Hispanic) from anonymized, randomly-sampled specimens. Data production utilized an 8-amplicon, 135 sequencing reaction Sanger-based protocol, performed in semi-automated fashion on robotic instrumentation. Data review followed an intensive multi-step strategy that included a minimum of three independent reviews of the raw data at two laboratories; repeat screenings of all insertions, deletions, heteroplasmies, transversions and any additional private mutations; and a check for phylogenetic feasibility. For all three populations, nearly complete resolution of the haplotypes was achieved with full mtGenome sequences: 90.3–98.8% of haplotypes were unique per population, an improvement of 7.7–29.2% over control region sequencing alone, and zero haplotypes overlapped between populations. Inferred maternal biogeographic ancestry frequencies for each population and heteroplasmy rates in the control region were generally consistent with published datasets. In the coding region, nearly 90% of individuals exhibited length heteroplasmy in the 12418-12425 adenine homopolymer; and despite a relatively high rate of point heteroplasmy (23.8% of individuals across the entire molecule), coding region point heteroplasmies shared by more than one individual were notably absent, and transversion-type heteroplasmies were extremely rare. The ratio of nonsynonymous to synonymous changes among point heteroplasmies in the protein-coding genes (1:1.3) and average pathogenicity scores in comparison to data reported for complete substitutions in previous studies seem to provide some additional support for the role of purifying selection in the evolution of the human mtGenome. Overall, these thoroughly vetted full mtGenome population reference data can serve as a standard against which the quality and features of future mtGenome datasets (especially those developed via massively parallel sequencing) may be evaluated, and will provide a solid foundation for the generation of complete mtGenome haplotype frequency estimates for forensic applications.     

Development of forensic-quality full mtGenome haplotypes: Success rates with low template specimens

Forensic mitochondrial DNA (mtDNA) testing requires appropriate, high quality reference population data for estimating the rarity of questioned haplotypes and, in turn, the strength of the mtDNA evidence. Available reference databases (SWGDAM, EMPOP) currently include information from the mtDNA control region; however, novel methods that quickly and easily recover mtDNA coding region data are becoming increasingly available. Though these assays promise to both facilitate the acquisition of mitochondrial genome (mtGenome) data and maximize the general utility of mtDNA testing in forensics, the appropriate reference data and database tools required for their routine application in forensic casework are lacking. To address this deficiency, we have undertaken an effort to: (1) increase the large-scale availability of high-quality entire mtGenome reference population data, and (2) improve the information technology infrastructure required to access/search mtGenome data and employ them in forensic casework.Here, we describe the application of a data generation and analysis workflow to the development of more than 400 complete, forensic-quality mtGenomes from low DNA quantity blood serum specimens as part of a U.S. National Institute of Justice funded reference population databasing initiative. We discuss the minor modifications made to a published mtGenome Sanger sequencing protocol to maintain a high rate of throughput while minimizing manual reprocessing with these low template samples. The successful use of this semi-automated strategy on forensic-like samples provides practical insight into the feasibility of producing complete mtGenome data in a routine casework environment, and demonstrates that large (>2 kb) mtDNA fragments can regularly be recovered from high quality but very low DNA quantity specimens. Further, the detailed empirical data we provide on the amplification success rates across a range of DNA input quantities will be useful moving forward as PCR-based strategies for mtDNA enrichment are considered for targeted next-generation sequencing workflows.     

Massively parallel sequencing of complete mitochondrial genomes from hair shaft samples

Though shed hairs are one of the most commonly encountered evidence types, they are among the most limited in terms of DNA quantity and quality. As a result, DNA testing has historically focused on the recovery of just about 600 base pairs of the mitochondrial DNA control region. Here, we describe our success in recovering complete mitochondrial genome (mtGenome) data (∼16,569 bp) from single shed hairs. By employing massively parallel sequencing (MPS), we demonstrate that particular hair samples yield DNA sufficient in quantity and quality to produce 2–3 kb mtGenome amplicons and that entire mtGenome data can be recovered from hair extracts even without PCR enrichment. Most importantly, we describe a small amplicon multiplex assay comprised of sixty-two primer sets that can be routinely applied to the compromised hair samples typically encountered in forensic casework. In all samples tested here, the MPS data recovered using any one of the three methods were consistent with the control Sanger sequence data developed from high quality known specimens. Given the recently demonstrated value of complete mtGenome data in terms of discrimination power among randomly sampled individuals, the possibility of recovering mtGenome data from the most compromised and limited evidentiary material is likely to vastly increase the utility of mtDNA testing for hair evidence.     

Reprint of: Evaluation of next generation mtGenome sequencing using the Ion Torrent Personal Genome Machine (PGM)

Insights into the human mitochondrial phylogeny have been primarily achieved by sequencing full mitochondrial genomes (mtGenomes). In forensic genetics (partial) mtGenome information can be used to assign haplotypes to their phylogenetic backgrounds, which may, in turn, have characteristic geographic distributions that would offer useful information in a forensic case. In addition and perhaps even more relevant in the forensic context, haplogroup-specific patterns of mutations form the basis for quality control of mtDNA sequences. The current method for establishing (partial) mtDNA haplotypes is Sanger-type sequencing (STS), which is laborious, time-consuming, and expensive. With the emergence of Next Generation Sequencing (NGS) technologies, the body of available mtDNA data can potentially be extended much more quickly and cost-efficiently. Customized chemistries, laboratory workflows and data analysis packages could support the community and increase the utility of mtDNA analysis in forensics. We have evaluated the performance of mtGenome sequencing using the Personal Genome Machine (PGM) and compared the resulting haplotypes directly with conventional Sanger-type sequencing. A total of 64 mtGenomes (>1 million bases) were established that yielded high concordance with the corresponding STS haplotypes (<0.02% differences). About two-thirds of the differences were observed in or around homopolymeric sequence stretches. In addition, the sequence alignment algorithm employed to align NGS reads played a significant role in the analysis of the data and the resulting mtDNA haplotypes. Further development of alignment software would be desirable to facilitate the application of NGS in mtDNA forensic genetics.     

Reprint of: Evaluation of next generation mtGenome sequencing using the Ion Torrent Personal Genome Machine (PGM)

Insights into the human mitochondrial phylogeny have been primarily achieved by sequencing full mitochondrial genomes (mtGenomes). In forensic genetics (partial) mtGenome information can be used to assign haplotypes to their phylogenetic backgrounds, which may, in turn, have characteristic geographic distributions that would offer useful information in a forensic case. In addition and perhaps even more relevant in the forensic context, haplogroup-specific patterns of mutations form the basis for quality control of mtDNA sequences. The current method for establishing (partial) mtDNA haplotypes is Sanger-type sequencing (STS), which is laborious, time-consuming, and expensive. With the emergence of Next Generation Sequencing (NGS) technologies, the body of available mtDNA data can potentially be extended much more quickly and cost-efficiently. Customized chemistries, laboratory workflows and data analysis packages could support the community and increase the utility of mtDNA analysis in forensics. We have evaluated the performance of mtGenome sequencing using the Personal Genome Machine (PGM) and compared the resulting haplotypes directly with conventional Sanger-type sequencing. A total of 64 mtGenomes (>1 million bases) were established that yielded high concordance with the corresponding STS haplotypes (<0.02% differences). About two-thirds of the differences were observed in or around homopolymeric sequence stretches. In addition, the sequence alignment algorithm employed to align NGS reads played a significant role in the analysis of the data and the resulting mtDNA haplotypes. Further development of alignment software would be desirable to facilitate the application of NGS in mtDNA forensic genetics.     

Evaluation of next generation mtGenome sequencing using the Ion Torrent Personal Genome Machine (PGM)

Insights into the human mitochondrial phylogeny have been primarily achieved by sequencing full mitochondrial genomes (mtGenomes). In forensic genetics (partial) mtGenome information can be used to assign haplotypes to their phylogenetic backgrounds, which may, in turn, have characteristic geographic distributions that would offer useful information in a forensic case. In addition and perhaps even more relevant in the forensic context, haplogroup-specific patterns of mutations form the basis for quality control of mtDNA sequences. The current method for establishing (partial) mtDNA haplotypes is Sanger-type sequencing (STS), which is laborious, time-consuming, and expensive. With the emergence of Next Generation Sequencing (NGS) technologies, the body of available mtDNA data can potentially be extended much more quickly and cost-efficiently. Customized chemistries, laboratory workflows and data analysis packages could support the community and increase the utility of mtDNA analysis in forensics. We have evaluated the performance of mtGenome sequencing using the Personal Genome Machine (PGM) and compared the resulting haplotypes directly with conventional Sanger-type sequencing. A total of 64 mtGenomes (>1 million bases) were established that yielded high concordance with the corresponding STS haplotypes (<0.02% differences). About two-thirds of the differences were observed in or around homopolymeric sequence stretches. In addition, the sequence alignment algorithm employed to align NGS reads played a significant role in the analysis of the data and the resulting mtDNA haplotypes. Further development of alignment software would be desirable to facilitate the application of NGS in mtDNA forensic genetics.     

Mitochondrial DNA heteroplasmy in the emerging field of massively parallel sequencing

Long an important and useful tool in forensic genetic investigations, mitochondrial DNA (mtDNA) typing continues to mature. Research in the last few years has demonstrated both that data from the entire molecule will have practical benefits in forensic DNA casework, and that massively parallel sequencing (MPS) methods will make full mitochondrial genome (mtGenome) sequencing of forensic specimens feasible and cost-effective. A spate of recent studies has employed these new technologies to assess intraindividual mtDNA variation. However, in several instances, contamination and other sources of mixed mtDNA data have been erroneously identified as heteroplasmy. Well vetted mtGenome datasets based on both Sanger and MPS sequences have found authentic point heteroplasmy in approximately 25% of individuals when minor component detection thresholds are in the range of 10–20%, along with positional distribution patterns in the coding region that differ from patterns of point heteroplasmy in the well-studied control region. A few recent studies that examined very low-level heteroplasmy are concordant with these observations when the data are examined at a common level of resolution. In this review we provide an overview of considerations related to the use of MPS technologies to detect mtDNA heteroplasmy. In addition, we examine published reports on point heteroplasmy to characterize features of the data that will assist in the evaluation of future mtGenome data developed by any typing method.     

mtDNA sequence diversity of Hazara ethnic group from Pakistan

The present study was undertaken to investigate mitochondrial DNA (mtDNA) control region sequences of Hazaras from Pakistan, so as to generate mtDNA reference database for forensic casework in Pakistan and to analyze phylogenetic relationship of this particular ethnic group with geographically proximal populations. Complete mtDNA control region (nt 16024-576) sequences were generated through Sanger Sequencing for 319 Hazara individuals from Quetta, Baluchistan. The population sample set showed a total of 189 distinct haplotypes, belonging mainly to West Eurasian (51.72%), East & Southeast Asian (29.78%) and South Asian (18.50%) haplogroups. Compared with other populations from Pakistan, the Hazara population had a relatively high haplotype diversity (0.9945) and a lower random match probability (0.0085). The dataset has been incorporated into EMPOP database under accession number EMP00680. The data herein comprises the largest, and likely most thoroughly examined, control region mtDNA dataset from Hazaras of Pakistan.     

Concordance and reproducibility of a next generation mtGenome sequencing method for high-quality samples using the Illumina MiSeq

Sanger-type sequencing (STS) of mitochondrial DNA (mtDNA), specifically the control region (CR), is routinely employed in forensics in human identification and missing persons scenarios. Yet next-generation sequencing (NGS) has the potential to overcome some of the major limitations of STS processing, permitting reasonable paths forward for full mitochondrial genome (mtGenome) sequencing, while also offering higher-throughput and higher sensitivity capabilities. To establish the accuracy and reproducibility of NGS for the development of mtDNA data, 90 DNA extracts that were previously used to generate forensic quality full mtGenomes using STS were sequenced using Nextera XT library preparation and the Illumina MiSeq. Using the same amplicon product, replicate library sets were generated and sequenced at different laboratories, and analysis was performed in replicate using the CLC Genomics Workbench. Both sequencing sets resulted in 99.998% of positions with greater than 10X coverage when 96 samples (including controls) were multiplexed. Overall, 99.9996% concordance was observed between the NGS data and the STS data for the full mtGenome. The only “discordant” calls involved low level point heteroplasmies, with the differences resulting from stochastic variation and/or the increased sensitivity of NGS. Higher sensitivity also allowed for the detection of a mixed sample previously not detected with STS. Additionally, variant calls were reproducible between sequencing sets and between software analysis versions with the variant frequency only differing by 0.23% and 0.01%, respectively. Further validation studies and specialized software functionality tailored to forensic practice should facilitate the incorporation of NGS processing into standard casework applications. The data herein comprise the largest, and likely most thoroughly examined, complete mtGenome STS-NGS concordance dataset available.     

A rapid mtDNA assay of 22 SNPs in one multiplex reaction increases the power of forensic testing in European Caucasians

We have developed a multiplex mitochondrial (mtDNA) assay of 21 coding region single nucleotide polymorphisms (SNPs) and one control region SNP outside hypervariable region 1 (HVR1) and hypervariable region 2 (HVR2) that can be amplified in a single reverse touchdown polymerase chain reaction. Single base extension using the SNaPshot technique is also carried out as one multiplex. Besides the nine major European haplogroups (i.e. H, I, J, K, T, U, V, W, and X), 16 additional subclades (i.e. N1, X2, X2b, U2′-4/7′-9′, J/T, J1, J1c, HV, H1, H1a1, H1c, H3, H4, H6a, H7a H10) can be detected and classified into a phylogenetic mtDNA tree. By analyzing 130 Caucasoid samples from Germany, 36 different haplotypes were found resulting in a power of discrimination of 93.2%. Although 49% of all samples belonged to superhaplogroup H, the most common haplotype, i.e., haplogroup-specific SNPs plus haplogroup unspecific SNPs, had a frequency of only 18%. This assay is applicable for high-throughput mtDNA analysis and forensic mass screening. It will give additional information to the common control region sequencing of HVR1 and HVR2. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Phylogenetic classification of Japanese mtDNA assisted by complete mitochondrial DNA sequences

We investigated control and coding region polymorphisms in mitochondrial DNA (mtDNA) in 100 unrelated individuals from a Japanese population and determined the basal phylogenetic haplogroup lineages in all samples under updated information. Many of the basal phylogenetic haplogroup lineages assigned on East Asian mtDNA haplogroups corresponded to those previously established. However, new haplogroup lineages such as M7a2a, M7a2b, M7a2*, M7c1b, M11b2*, G2b*, D4c1b1a, D4g2b, A4*, A9, N9b*, B4d1, B4d2, and F1e were identified and established by complete sequencing. Although sequence comparison of the 1.15-kb control region identified 84 mitochondrial haplotypes, examination of coding region polymorphisms increased the total number of haplotypes to 91. Determination of the basal haplogroup lineages increased the discrimination power of mtDNA polymorphisms for personal identification and their usefulness in determining geographic origin in forensic casework in Japanese and other East Asian populations.     

Mitochondrial control region diversity in Sindhi ethnic group of Pakistan

The entire mitochondrial DNA control region (nt 16024–576) of 88 unrelated individuals of Sindhi ethnic group residing in different parts of Sindh province of Pakistan was sequenced. Out of 66 different observed haplotypes 50 were unique and 16 were shared by more than one individual. Results showed admixture of mtDNA pool constituting the haplogroups derived mainly from South Asia (47.6%) and West Eurasian (35.7%) whereas the contribution of the African haplogroup was very small (2.4%). High values of genetic diversity (0.992), power of discrimination (0.981) and low value of random match probability (0.018) indicates that mtDNA analysis for this population can effectively be used for forensic casework. The results are valuable contribution towards building mtDNA population variation database for this particular ethnic group from Pakistan.     

Mitochondrial DNA variation in the Sindh population of Pakistan

The present study was undertaken to investigate the control region of mitochondrial DNA for forensic discrimination and to explore the ethno-linguistic affiliations among ethnic groups of Sindh province, Pakistan. A total of 115 individuals, from six major ethnic/isonym groups, namely, Bijarani, Chandio, Ghallu, Khoso, Nasrani and Solangi, have been studied. We identified 88 haplotypes, defined by the particular set of nucleotides; of these, 70 haplotypes were unique in the investigated population. In addition, 82% sequences were observed once, 12% twice and 5.2% thrice. The most common South Asian haplogroup in six ethnic groups of Sindh, are; M (42%) and R (6.9%), whereas West Eurasian haplogroups were N (6.9%), W (6.9%), J (1.7%), U (23.4%), H (9.5%) and T (0.86%). A random match probability between two unrelated individuals was found to be 0.06%, while genetic diversity varied from 0.991 to 0.998. The high nucleotide diversity and the low random match probability of the mtDNA control region make it a beneficial tool for forensic discrimination as well as useful to the evolutionary biologist. This work is an important contribution towards establishing a National Mitochondrial DNA Database in Pakistan. Currently, the Human Genetic department at the University of Health Sciences Lahore has the depository.     

Resolving mitochondrial haplogroups B2 and B4 with next-generation mitogenome sequencing to distinguish Native American from Asian haplotypes

Mitochondrial haplogroup information can be useful in forensic contexts that rely primarily on mitochondrial DNA (mtDNA) testing, which often involve limited or degraded DNA. Due to the phylogeographic patterning of mtDNA in human populations, mitochondrial haplogroups are indicative of maternal ancestry (as mtDNA is a maternally inherited marker). In certain circumstances, maternal ancestry inferred from mitochondrial haplogrouping could be beneficial to forensic investigations. For example, ancestry information could assist in the identification of unknown service members from past conflicts, such as the World War II Battle of Tarawa involving American and Japanese forces. In this context, it could be useful to distinguish Native American mtDNA from Asian mtDNA to bolster the anthropological and circumstantial evidence leading to an identification or foreign national determination. Although most of the founding Native American haplogroups contain diagnostic variants in the mitochondrial control region (CR), haplogroup B2 does not, and this makes it more difficult to distinguish B2 from the parental B4 and closely related B4b haplogroups found in Asia. In this paper, the amount of mtDNA information required to distinguish Native American haplotypes from Asian haplotypes within haplogroup B was examined. Fifty-six samples belonging to subtypes of B2 and B4 were sequenced for the entire mitogenome. Haplogroups were estimated from three ranges of mitochondrial DNA (HV1 and 2, CR, and full mitogenome). Half of the samples could not be precisely haplogrouped without full mitogenome data, although enough variants were often provided to make an accurate B2 versus B4 distinction. Native American B2 haplotypes were distinguishable using CR data alone in 82% of samples, though the remaining samples required full mitogenome data for haplogroup B2 designation. The use of full mitogenome data consistently enables accurate haplogroup determination, and opens the possibility for gaining information on maternal ancestry.     

Homogeneity in mitochondrial DNA control region sequences in Swedish subpopulations

In order to promote mitochondrial DNA (mtDNA) testing in Sweden we have typed 296 Swedish males, which will serve as a Swedish mtDNA frequency database. The tested males were taken from seven geographically different regions representing the contemporary Swedish population. The complete mtDNA control region was typed and the Swedish population was shown to have high haplotype diversity with a random match probability of 0.5%. Almost 47% of the tested samples belonged to haplogroup H and further haplogroup comparison with worldwide populations clustered the Swedish mtDNA data together with other European populations. AMOVA analysis of the seven Swedish subregions displayed no significant maternal substructure in Sweden (F _ST = 0.002). Our conclusion from this study is that the typed Swedish individuals serve as good representatives for a Swedish forensic mtDNA database. Some caution should, however, be taken for individuals from the northernmost part of Sweden (provinces of Norrbotten and Lapland) due to specific demographic conditions. Furthermore, our analysis of a small sample set of a Swedish Saami population confirmed earlier findings that the Swedish Saami population is an outlier among European populations.     

A collaborative EDNAP exercise on SNaPshot™-based mtDNA control region typing

A collaborative European DNA Profiling (EDNAP) Group exercise was undertaken to assess the performance of an earlier described SNaPshot™-based screening assay (denoted mini-mtSNaPshot) (Weiler et al., 2016) [1] that targets 18 single nucleotide polymorphism (SNP) positions in the mitochondrial (mt) DNA control region and allows for discrimination of major European mtDNA haplogroups. Besides the organising laboratory, 14 forensic genetics laboratories were involved in the analysis of 13 samples, which were centrally prepared and thoroughly tested prior to shipment. The samples had a variable complexity and comprised straightforward single-source samples, samples with dropout or altered peak sizing, a point heteroplasmy and two-component mixtures resulting in one to five bi-allelic calls. The overall success rate in obtaining useful results was high (97.6%) given that some of the participating laboratories had no previous experience with the typing technology and/or mtDNA analysis. The majority of the participants proceeded to haplotype inference to assess the feasibility of assigning a haplogroup and checking phylogenetic consistency when only 18 SNPs are typed. To mimic casework procedures, the participants compared the SNP typing data of all 13 samples to a set of eight mtDNA reference profiles that were described according to standard nomenclature (Parson et al., 2014) [2], and indicated whether these references matched each sample or not. Incorrect scorings were obtained for 2% of the comparisons and derived from a subset of the participants, indicating a need for training and guidelines regarding mini-mtSNaPshot data interpretation.     

Hungarian mtDNA population databases from Budapest and the Baranya county Roma

To facilitate forensic mtDNA testing in Hungary, we have generated control region databases for two Hungarian populations: 211 individuals were sampled from the urban Budapest population and 208 individuals were sampled from a Romani (“gypsy”) population in Baranya county. Sequences were generated using a highly redundant approach to minimize potential database errors. The Budapest population had high sequence diversity with 180 lineages, 183 polymorphic positions, and a random match probability of 1%. In contrast, the Romani population exhibited low sequence diversity, with only 56 lineages, 109 segregating sites, and a random match probability of 8.8%. The mtDNA haplogroup compositions of the two populations were also distinct, with the large proportion of haplogroup M samples (35%) in the Roma the most obvious difference between the two populations. These factors highlight the importance of considering population structure when generating reference databases for forensic testing purposes. Comparisons between our Romani population sample and other published data indicate the need for heightened caution when sampling and using mtDNA databases of small endogamous populations. The Romani populations that we compared showed significant departures from genetic uniformity. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Mitochondrial DNA analysis of a Viking age mass grave in Sweden

In 1998, a Viking Age mass grave was discovered and excavated at St. Laurence´s churchyard in Sigtuna, Sweden. The excavated bones underwent osteoarchaeological analysis and were assigned to at least 19 individuals. Eleven skeletons showed sharp force trauma from bladed weapons. Mass graves are an unusual finding from this time period, making the burial context extraordinary. To investigate a possible maternal kinship among the individuals, bones and teeth from the skeletal remains were selected for mitochondrial DNA (mtDNA) analysis. Sanger sequencing of short stretches of the hypervariable segments I and II (HVS-I and HVS-II) was performed. A subset of the samples was also analysed by massively parallel sequencing analysis (MPS) of the entire mtDNA genome using the Precision ID mtDNA Whole Genome Panel. A total of 15 unique and three shared mtDNA profiles were obtained. Based on a combination of genetic and archaeological data, we conclude that a minimum of 20 individuals was buried in the mass grave. The majority of the individuals were not maternally related. However, two possible pairs of siblings or mother-child relationships were identified. All individuals were assigned to West Eurasian haplogroups, with a predominance of haplogroup H. Although the remains showed an advanced level of DNA degradation, the combined use of Sanger sequencing and MPS with the Precision ID mtDNA Whole Genome Panel revealed at least partial mtDNA data for all samples.     

Haplotype diversity in mitochondrial DNA hypervariable region in a population of southeastern Brazil

Brazilian population derives from Native Amerindians, Europeans, and Africans. Southeastern Brazil is the most populous region of the country. The present study intended to characterize the maternal genetic ancestry of 290 individuals from southeastern (Brazil) population. Thus, we made the sequencing of the three hypervariable regions (HV1, HV2, and HV3) of the mitochondrial DNA (mtDNA). The statistical analyses were made using Arlequin software, and the median-joining haplotype networks were generated using Network software. The analysis of three hypervariable regios showed 230 (79.3 %) unique haplotypes and the most common haplotype was “263G” carried by 12 (4.1 %) individuals. The strikingly high variability generated by intense gene flow is mirrored in a high sequence diversity (0.9966?±?0.0010), and the probability of two random individuals showing identical mtDNA haplotypes were 0.0068. The analysis of haplogroup distribution revealed that 36.9 % (n?=?107) presented Amerindian haplogroups, 35.2 % (n?=?102) presented African haplogroups, 27.6 % (n?=?80) presented European haplogroups, and one (0.3 %) individual presented East Asian haplogroup, evidencing that the southeastern population is extremely heterogeneous and the coexistence of matrilineal lineages with three different phylogeographic origins. The genetic diversity found in the mtDNA control region in the southeastern Brazilian population reinforces the importance of increased national database in order to be important and informative in forensic cases.