Haplogrouping mitochondrial DNA sequences in Legal Medicine/Forensic Genetics

Deciphering the maternal ancestral lineage of Greek Cypriots, Armenian Cypriots and Maronite Cypriots

Cyprus was conquered from several populations because of its special geographical location. In this study, 406 unrelated Cypriot samples were tested based on their mitochondrial DNA. In more detail, 185 were Greek Cypriots, 114 Armenian Cypriots and 107 Maronite Cypriots. This is the first time where the mitochondrial DNA of Greek Cypriots, Armenian Cypriots and Maronite Cypriots is compared with the aim of characterizing the maternal ancestry of Cypriots. The control region of the mtDNA is the most informative in terms of studying maternal ancestry and consists of three hypervariable regions (HVS-I, HVS-II, HVS-III). The hypervariable regions can provide important information regarding the maternal ancestor of the tested samples. The entire control region of the mtDNA was used to determine the mitotypes and subsequently the haplogroups of all the Cypriot DNA samples. Based on the aforementioned analyses, Greek Cypriots were found to be genetically closer to Armenian Cypriots, while Greek Cypriots and Armenian Cypriots showed moderate genetic differentiation with Maronite Cypriots. The most prevalent haplogroups among Cypriots were haplogroups H and U, while R0 is common but in different frequencies for Greek Cypriots, Armenian Cypriots and Maronite Cypriots. It is proposed that the maternal ancestor may have originated during the Neolithic period and/or the Bronze age.

Analysis of maternal genetic structure of mitochondrial DNA control region from Tai-Kadai-speaking Buyei population in southwestern China

BACKGROUND

Even though the Buyei are a recognised ethnic group in southwestern China, there hasn't been much work done on forensic population genetics, notably using mitochondrial DNA. The sequences and haplogroups of mitochondrial DNA control regions of the Buyei peoples were studied to provide support for the establishment of a reference database for forensic DNA analysis in East Asia.

METHODS AND RESULTS

The mitochondrial DNA control region sequences of 200 Buyei individuals in Guizhou were investigated. The haplotype frequencies and haplogroup distribution of the Buyei nationality in Guizhou were calculated. At the same time, the paired Fst values of the study population and other populations around the world were computed, to explore their genetic polymorphism and population relationship. A total of 179 haplotypes were detected in the Buyei population, with frequencies of 0.005-0.015. All haplotypes were assigned to 89 different haplogroups. The haplotype diversity and random matching probability were 0.999283 and 0.0063, respectively. The paired Fst genetic distances and correlation p-values among the 54 populations revealed that the Guizhou Buyei was most closely related to the Henan Han and the Guizhou Miao, and closer to the Hazara population in Pakistan and the Chiang Mai population.

CONCLUSIONS

The study of mitochondrial DNA based on the maternal genetic structure of the Buyei nationality in Guizhou will benefit the establishment of an East Asian forensic DNA reference database and provide a reference for anthropological research in the future.

HaploCart: Human mtDNA haplogroup classification using a pangenomic reference graph human mtDNA haplogroup inference

Current mitochondrial DNA (mtDNA) haplogroup classification tools map reads to a single reference genome and perform inference based on the detected mutations to this reference. This approach biases haplogroup assignments towards the reference and prohibits accurate calculations of the uncertainty in assignment. We present HaploCart, a probabilistic mtDNA haplogroup classifier which uses a pangenomic reference graph framework together with principles of Bayesian inference. We demonstrate that our approach significantly outperforms available tools by being more robust to lower coverage or incomplete consensus sequences and producing phylogenetically-aware confidence scores that are unbiased towards any haplogroup. HaploCart is available both as a command-line tool and through a user-friendly web interface. The C++ program accepts as input consensus FASTA, FASTQ, or GAM files, and outputs a text file with the haplogroup assignments of the samples along with the level of confidence in the assignments. Our work considerably reduces the amount of data required to obtain a confident mitochondrial haplogroup assignment.

Whole mitochondrial genome analysis of Tai-Kadai-speaking populations in Southwest China

As a single matrilineal gene, human mitochondrial DNA plays a very important role in the study of population genetics. The whole mitogenome sequences of 287 individuals of the Tai-Kadai-speaking population in Guizhou were obtained. It was discovered that there were 82, 104, and 94 haplotypes in 83 Bouyei individuals, 107 Dong individuals, and 97 Sui individuals, respectively; and the haplotype diversity in Bouyei, Dong, and Sui groups was 1.000 ± 0.02, 0.9993 ± 0.0015, and 0.999 ± 0.002, respectively. The result of neutrality tests of the Tai-Kadai-speaking population in Guizhou showed significant negative values, and the analysis of mismatch distribution showed an obvious unimodal distribution. The results implied that Guizhou Tai-Kadai-speaking populations had high genetic diversities and may have experienced recent population expansion. In addition, the primary haplogroups of studied populations were M*, F, B, D, and R*, implying that they may origin from Southern China. The matrilineal genetic structure of the Tai-Kadai-speaking populations in Guizhou was analyzed by merging the mitogenome data of 79 worldwide populations as reference data. The results showed that there were close relationships between studied populations and other Tai-Kadai as well as some Austronesian populations in East and Southeast Asia. Overall, the mitogenome data generated in this study will provide important data for the study of genetic structure of Tai-Kadai speaking populations.

Mitochondrial Haplogroup Classification of Ancient DNA Samples Using Haplotracker

Mitochondrial DNA haplogroup classification is used to study maternal lineage of ancient human populations. The haplogrouping of ancient DNA is not easy because the DNA is usually found in small pieces in limited quantities. We have developed Haplotracker, a straightforward and efficient high-resolution haplogroup classification tool optimized specifically for ancient DNA samples. Haplotracker offers a user-friendly input interface for multiple mitochondrial DNA sequence fragments in a sample. It provides accurate haplogroup classification with full-length mitochondrial genome sequences and provides high-resolution haplogroup predictions for some fragmented control region sequences using a novel algorithm built on Phylotree mtDNA Build 17 (Phylotree) and our haplotype database (n = 118,869). Its performance for accuracy was demonstrated to be high through haplogroup classification using 8,216 Phylotree full-length and control region mitochondrial DNA sequences compared with HaploGrep 2, one of the most accurate current haplogroup classifiers. Haplotracker provides a novel haplogroup tracking solution for fragmented sequences to track subhaplogroups or verify the haplogroups efficiently. Using Haplotracker, we classified mitochondrial haplogroups to the final subhaplogroup level in nine ancient DNA samples extracted from human skeletal remains found in 2,000-year-old elite Xiongnu cemetery in Northeast Mongolia. Haplotracker can be freely accessed at https://haplotracker.cau.ac.kr.

The mitochondrial DNA control region sequences from the Chinese Sui population of southwestern China

BACKGROUND

Sui people are officially recognised people living in southwest China, but there has been a lack of genetic research, especially based on mitochondrial DNA data.

AIM

To study the sequences and haplogroups of the mitochondrial DNA control region in a typical Sui population, with the aim of helping to promote the establishment of a forensic DNA analysis reference database in East Asia.

SUBJECTS AND METHODS

We analysed 201 Sui individuals and observed the sequences of the mitochondrial DNA control region. We calculated and explained the haplotype frequencies, haplogroup distribution and pairwise Fst values between the Sui and 47 other populations in the world, in order to explore genetic polymorphisms and population relationships.

RESULTS

161 haplotypes were found in the Sui population, with frequencies of 0.0049-0.0199. All samples were assigned to 80 different haplogroups. The haplotype diversity and random matching probability were 0.999938 and 0.024729, respectively. The pairwise Fst values and correlation p-values of 48 populations showed that the Sui population was most closely related to the Miao population in Guizhou and the Han population in Henan, and closer to the Punjab population and Pukhtunkhwa population in Pakistan, and was significantly different from the other 43 groups. Compared with the other 43 groups, it is relatively isolated.

CONCLUSION

Our results show that the study of mitochondrial DNA based on the analysis of matrilineal genetic structure of the Sui population can help to promote the establishment of a forensic DNA reference database in East Asia and provide reference for future anthropological research.

Fine-Tuning Phylogenetic Alignment and Haplogrouping of mtDNA Sequences

In this paper, we present a new algorithm for alignment and haplogroup estimation of mitochondrial DNA (mtDNA) sequences. Based on 26,011 vetted full mitogenome sequences, we refined the 5435 original haplogroup motifs of Phylotree Build 17 without changing the haplogroup nomenclature. We adapted 430 motifs (about 8%) and added 966 motifs for yet undetermined subclades. In summary, this led to an 18% increase of haplogroup defining motifs for full mitogenomes and a 30% increase for the mtDNA control region that is of interest for a variety of scientific disciplines, such as medical, population and forensic genetics. The new algorithm is implemented in the EMPOP mtDNA database and is freely accessible.

Developmental validation of the MGIEasy Signature Identification Library Prep Kit, an all-in-one multiplex system for forensic applications

Analyzing genetic markers in nuclear and mitochondrial genomes is helpful in various forensic applications, such as individual identifications and kinship analyses. However, most commercial kits detect these markers separately, which is time-consuming, laborious, and more error-prone (mislabelling, contamination, ...). The MGIEasy Signature Identification Library Prep Kit (hereinafter "MGIEasy identification system"; MGI Tech, Shenzhen, China) has been designed to provide a simple, fast, and robust way to detect appropriate markers in one multiplex PCR reaction: 52 autosomal STRs, 27 X-chromosomal STRs, 48 Y-chromosomal STRs, 145 identity-informative SNPs, 53 ancestry-informative SNPs, 29 phenotype-informative SNPs, and the hypervariable regions of mitochondrial DNA (mtDNA). Here, we validated the performance of MGIEasy identification system following the guidelines of the Scientific Working Group on DNA Analysis Methods (SWGDAM), assessing species specificity, sensitivity, mixture identification, stability under non-optimal conditions (degraded samples, inhibitor contamination, and various substrates), repeatability, and concordance. Libraries prepared using MGIEasy identification system were sequenced on a MGISEQ-2000 instrument (MGI Tech). MGIEasy-derived STR, SNP, and mtDNA genotypes were highly concordant with CE-based STR genotypes (99.79%), MiSeq FGx-based SNP genotypes (99.78%), and Sanger-based mtDNA genotypes (100%), respectively. This system was strongly human-specific, resistant to four common PCR inhibitors, and reliably amplified both low quantities of DNA (as low as 0.125 ng) and degraded DNA (~ 150 nt). Most of the unique alleles from the minor contributor were detected in 1:10 male-female and male-male mixtures; some minor Y-STR alleles were even detected in 1:1000 male-female mixtures. MGIEasy also successfully directly amplified markers from blood stains on FTA cards, filter papers, and swabs. Thus, our results demonstrated that MGIEasy identification system was suitable for use in forensic analyses due to its robust and reliable performance on samples of varying quality and quantity.

Match statistics for sequence-based alleles in profiles from forensic PCR-mps kits

The first autosomal sequence-based allele (aka SNP-STR haplotype) frequency database for forensic massively parallel sequencing (MPS) has been published, thereby removing one of the remaining barriers to implementing MPS in casework. The database was developed using a specific set of flank trim sites. If different trim sites or different kits with different primers are used for casework, then SNP-STR haplotypes may be detected that do not have frequencies in the database. We describe a procedure to address calculation of match probabilities when casework samples are generated using an MPS kit with different trim sites than those present in the relevant population frequency database. The procedure provides a framework for comparison of any MPS kit or database combination while also accommodating comparison of MPS and CE profiles.

The mitochondrial DNA control region sequences from the Chinese Miao population of southeastern China

Background: Miao people are an officially recognised ethnic group living in southwest China, but have seldom been studied genetically, especially with respect to mtDNA data.Aim: To investigate the sequences and haplogroups of the mtDNA control region in a typical Miao population, with the aim of providing a good start for the expansion of the East Asian mtDNA reference database for forensic DNA analysis.Subjects and methods: We analysed 203 Miao individuals, looking at mtDNA control region sequences. We calculated and illustrated the haplotype frequencies, haplogroup distribution and pairwise Fst values between the Miao and six other worldwide populations to explore genetic polymorphisms and population relationships.Results: We observed 121 haplotypes with corresponding frequencies ranging from 0.0049 to 0.0690 in the Miao population. All the samples were assigned to 71 different haplogroups. The haplotype diversity and the random match probability were estimated to be 0.9844 and 0.0204, respectively. The pairwise Fst values and associated p values among seven populations suggest that the Miao population has significant differences to the other six populations, and is relatively isolated compared with them.Conclusions: Our results suggest that frequency estimates for mtDNA haplotypes in Miao ethnic groups should be determined independently rather than being pooled with other populations.

A single nucleotide polymorphism panel for individual identification and ancestry assignment in Caucasians and four East and Southeast Asian populations using a machine learning classifier

Single nucleotide polymorphism (SNP) profiling is an effective means of individual identification and ancestry inferences in forensic genetics. This study established a SNP panel for the simultaneous individual identification and ancestry assignment of Caucasian and four East and Southeast Asian populations. We analyzed 220 SNPs (125 autosomal, 17 X-chromosomal, 30 Y-chromosomal, and 48 mitochondrial SNPs) of the DNA samples from 563 unrelated individuals of five populations (89 Caucasian, 234 Taiwanese Han, 90 Filipino, 79 Indonesian and 71 Vietnamese) and 18 degraded DNA samples. Informativeness for assignment (In) was used to select ancestry informative SNPs (AISNPs). A machine learning classifier, support vector machine (SVM), was used for ancestry assignment. Of the 220 SNPs, 62 were individual identification SNPs (IISNPs) (51 autosomal and 11 X-chromosomal SNPs) and 191 were AISNPs (100 autosomal, 13 X-chromosomal, 30 Y-chromosomal, and 48 mitochondrial SNPs). The 51 autosomal IISNPs offered cumulative random match probabilities (cRMPs) ranging from 1.56 × 10^-21 to 3.16 × 10^-22 among these five populations. Using AISNPs with the SVM, the overall accuracy rate of ancestry inference achieved in the testing dataset between Caucasian, Taiwanese Han, and Filipino populations was 88.9%, whereas it was 70.0% between Caucasians and each of the four East and Southeast Asian populations. For the 18 degraded DNA samples with incomplete profiling, the accuracy rate of ancestry assignment was 94.4%. We have developed a 220-SNP panel for simultaneous individual identification and ethnic origin differentiation between Caucasian and the four East and Southeast Asian populations. This SNP panel may assist with DNA analysis of forensic casework.

Next generation database search algorithm for forensic mitogenome analyses

Mitochondrial DNA (mtDNA) variation is being reported relative to the corrected version of the first sequenced human mitochondrial genome. A review of the existing literature across disciplines that employ mtDNA demonstrates that insertions and deletions are not reported in a standardized way. This may lead to false exclusions of identical sequences, unidentified matches in missing persons mtDNA databases, biased mtDNA database frequency estimates and overestimation of the genetic evidence. Seven years ago we introduced alignment-free database search software (SAM) and implemented it into the mtDNA database EMPOP (https://empop.online) to produce reliable and conservative frequency estimates that are required in the forensic context. However, ambiguity remained in how laboratories have been reporting mitotypes, as often more than one single alignment of a given mtDNA sequence was feasible. In order to overcome this limitation we here describe a concept and provide software for producing stable, harmonized phylogenetic alignment of mtDNA sequences for database searches. The new software SAM 2 will be made available via EMPOP and provide the user with the already established conservative frequency estimates. In addition, SAM 2 offers the rCRS-coded haplotype of a given mtDNA sequence following the established and widely accepted phylogenetic alignment. This provides the user with feedback on how mitotypes are stored in EMPOP and how they should be reported in order to harmonize nomenclature. Finally, this approach does not only permit reliable mtDNA nomenclature in forensics but invites related disciplines to take advantage of a standardized way of reporting mtDNA variation, thus closing the ranks between different genetic fields and supporting dialogue and collaboration between mtDNA scholars from various disciplines.

Development of a new screening method to determine the main 52 mitochondrial haplogroups through a single minisequencing reaction

This work presents the design, development and optimization of a screening method based on single-base extension sequencing to simultaneously analyze a panel of 52 mitochondrial SNPs. This enables to recognize the main mitochondrial haplogroups and to discriminate even between lineages from the same phylogenetic branch that diverged in different continents. The unavailability of individuals harboring infrequent variants was a limitation to optimize the panel. To overcome this, we have modified DNA by site-directed mutagenesis to create the unavailable allelic variants. This allowed us to verify the reliability of this panel and its usefulness to be applied in biomedicine, forensic and population genetic studies.

A phylogenetic approach for haplotype analysis of sequence data from complex mitochondrial mixtures

Massively parallel (next-generation) sequencing provides a powerful method to analyze DNA from many different sources, including degraded and trace samples. A common challenge, however, is that many forensic samples are often known or suspected mixtures of DNA from multiple individuals. Haploid lineage markers, such as mitochondrial (mt) DNA, are useful for analysis of mixtures because, unlike nuclear genetic markers, each individual contributes a single sequence to the mixture. Deconvolution of these mixtures into the constituent mitochondrial haplotypes is challenging as typical sequence read lengths are too short to reconstruct the distinct haplotypes completely. We present a powerful computational approach for determining the constituent haplotypes in massively parallel sequencing data from potentially mixed samples. At the heart of our approach is an expectation maximization based algorithm that co-estimates the overall mixture proportions and the source haplogroup for each read individually. This approach, implemented in the software package mixemt, correctly identifies haplogroups from mixed samples across a range of mixture proportions. Furthermore, our method can separate fragments in a mixed sample by the most likely originating contributor and generate reconstructions of the constituent haplotypes based on known patterns of mtDNA diversity.

MtDNA As a Cancer Marker: A Finally Closed Chapter?

Sequence alterations of the mitochondrial DNA (mtDNA) have been identified in many tu-mor types. Their nature is not entirely clear. Somatic mutation or shifts of heteroplasmic mtDNA vari-ants may play a role. These sequence alterations exhibit a sufficient frequency in all tumor types investi-gated thus far to justify their use as a tumor marker. This statement is supported by the high copy num-ber of mtDNA, which facilitates the detection of aberrant tumor-derived DNA in bodily fluids. This will be of special interest in tumors, which release a relatively high number of cells into bodily fluids, which are easily accessible, most strikingly in urinary bladder carcinoma. Due to the wide distribution of the observed base substitutions, deletions or insertions within the mitochondrial genome, high efforts for whole mtDNA sequencing (16.5 kb) from bodily fluids would be required, if the method would be in-tended for initial tumor screening. However, the usage of mtDNA for sensitive surveillance of known tumor diseases is a meaningful option, which may allow an improved non-invasive follow-up for the urinary bladder carcinoma, as compared to the currently existing cytological or molecular methods. Fol-lowing a short general introduction into mtDNA, this review demonstrates that the scenario of a sensi-tive cancer follow-up by mtDNA-analysis deserves more attention. It would be most important to inves-tigate precisely in the most relevant tumor types, if sequencing approaches in combination with simple PCR-assays for deletions/insertions in homopolymeric tracts has sufficient sensitivity to find most tu-mor-derived mtDNAs in bodily fluids.

Updating the African human mitochondrial DNA tree: Relevance to forensic and population genetics

Analysis of human mitochondrial DNA (mtDNA) variation plays an important role in forensic genetic investigations, especially in degraded biological samples and hair shafts. There are many issues of the mtDNA phylogeny that are of special interest to the forensic community, such as haplogroup classification or the post hoc investigation of potential errors in mtDNA datasets. We have analyzed >2200 mitogenomes of African ancestry with the aim of improving the known worldwide phylogeny. More than 300 new minor subclades were identified, and the Time to the Most Recent Common Ancestor (TMRCA) was estimated for each node of the phylogeny. Phylogeographic details are provided which might also be relevant to forensic genetics. The present study has special interest for forensic investigations because current analysis and interpretation of mtDNA casework rest on a solid worldwide phylogeny, as is evident from the role that phylogeny plays in popular resources in the field (e.g. PhyloTree), software (e.g. Haplogrep 2), and databases (e.g. EMPOP). Apart from this forensic genetic interest, we also highlight the impact of this research in anthropological studies, such as those related to the reconstruction of the transatlantic slave trade.

Massive parallel sequencing of human whole mitochondrial genomes with Ion Torrent technology: an optimized workflow for Anthropological and Population Genetics studies

Investigation of human mitochondrial DNA variation patterns and phylogeny has been extensively used in Anthropological and Population Genetics studies and sequencing the whole mitochondrial genome is progressively becoming the gold standard. Among the currently available massive parallel sequencing technologies, Ion Torrent™ semiconductor sequencing represents a promising approach for such studies. Nevertheless, an experimental protocol conceived to enable the achievement of both as high as possible yield and of the most homogeneous sequence coverage through the whole mitochondrial genome is still not available. The present work was thus aimed at improving the overall performance of whole mitochondrial genomes Ion Torrent™ sequencing, with special focus on the capability to obtain robust coverage and highly reliable variants calling. For this purpose, a series of cost-effective modifications in standard laboratory workflows was fine-tuned to optimize them for medium- and large-scale population studies. A total of 54 human samples were thus subjected to sequencing of the whole mitochondrial genome with the Ion Personal Genome Machine™ System in four distinct experiments and using Ion 314 chips. Seven of the selected samples were also characterized by means of conventional Sanger sequencing for the sake of comparison. Obtained results demonstrated that the implemented optimizations had definitely improved sequencing outputs in terms of both variants calling efficiency and coverage uniformity, enabling to setup an effective and accurate protocol for whole mitochondrial genome sequencing and a considerable reduction in experimental time consumption and sequencing costs.

EMPOP-quality mtDNA control region sequences from Kashmiri of Azad Jammu & Kashmir, Pakistan

The mitochondrial DNA (mtDNA) control region (nucleotide position 16024-576) sequences were generated through Sanger sequencing method for 317 self-identified Kashmiris from all districts of Azad Jammu & Kashmir Pakistan. The population sample set showed a total of 251 haplotypes, with a relatively high haplotype diversity (0.9977) and a low random match probability (0.54%). The containing matrilineal lineages belonging to three different phylogeographic origins of Western Eurasian (48.9%), South Asian (47.0%) and East Asian (4.1%). The present study was compared to previous data from Pakistan and other worldwide populations (Central Asia, Western Asia, and East & Southeast Asia). The dataset is made available through EMPOP under accession number EMP00679 and will serve as an mtDNA reference database in forensic casework in Pakistan.

Meta-Analysis of Mitochondrial DNA Variation in the Iberian Peninsula

The Iberian Peninsula has been the focus of attention of numerous studies dealing with mitochondrial DNA (mtDNA) variation, most of them targeting the control region segment. In the present study we sequenced the control region of 3,024 Spanish individuals from areas where available data were still limited. We also compiled mtDNA haplotypes from the literature involving 4,588 sequences and 28 population groups or small regions. We meta-analyzed all these data in order to shed further light on patterns of geographic variation, taking advantage of the large sample size and geographic coverage, in contrast with the atomized sampling strategy of previous work. The results indicate that the main mtDNA haplogroups show primarily clinal geographic patterns across the Iberian geography, roughly along a North-South axis. Haplogroup HV0 (where haplogroup U is nested) is more prevalent in the Franco Cantabrian region, in good agreement with previous findings that identified this area as a climate refuge during the Last Glacial Maximum (LGM), prior to a subsequent demographic re-expansion towards Central Europe and the Mediterranean. Typical sub-Saharan and North African lineages are slightly more prevalent in South Iberia, although at low frequencies; this pattern has been shaped mainly by the transatlantic slave trade and the Arab invasion of the Iberian Peninsula. The results also indicate that summary statistics that aim to measure molecular variation, or AMOVA, have limited sensitivity to detect population substructure, in contrast to patterns revealed by phylogeographic analysis. Overall, the results suggest that mtDNA variation in Iberia is substantially stratified. These patterns might be relevant in biomedical studies given that stratification is a common cause of false positives in case-control mtDNA association studies, and should be also considered when weighting the DNA evidence in forensic casework, which is strongly dependent on haplotype frequencies.

Analysis of uni and bi-parental markers in mixture samples: Lessons from the 22nd GHEP-ISFG Intercomparison Exercise

Since 1992, the Spanish and Portuguese-Speaking Working Group of the ISFG (GHEP-ISFG) has been organizing annual Intercomparison Exercises (IEs) coordinated by the Quality Service at the National Institute of Toxicology and Forensic Sciences (INTCF) from Madrid, aiming to provide proficiency tests for forensic DNA laboratories. Each annual exercise comprises a Basic (recently accredited under ISO/IEC 17043: 2010) and an Advanced Level, both including a kinship and a forensic module. Here, we show the results for both autosomal and sex-chromosomal STRs, and for mitochondrial DNA (mtDNA) in two samples included in the forensic modules, namely a mixture 2:1 (v/v) saliva/blood (M4) and a mixture 4:1 (v/v) saliva/semen (M8) out of the five items provided in the 2014 GHEP-ISFG IE. Discrepancies, other than typos or nomenclature errors (over the total allele calls), represented 6.5% (M4) and 4.7% (M8) for autosomal STRs, 15.4% (M4) and 7.8% (M8) for X-STRs, and 1.2% (M4) and 0.0% (M8) for Y-STRs. Drop-out and drop-in alleles were the main cause of errors, with laboratories using different criteria regarding inclusion of minor peaks and stutter bands. Commonly used commercial kits yielded different results for a micro-variant detected at locus D12S391. In addition, the analysis of electropherograms revealed that the proportions of the contributors detected in the mixtures varied among the participants. In regards to mtDNA analysis, besides important discrepancies in reporting heteroplasmies, there was no agreement for the results of sample M4. Thus, while some laboratories documented a single control region haplotype, a few reported unexpected profiles (suggesting contamination problems). For M8, most laboratories detected only the haplotype corresponding to the saliva. Although the GHEP-ISFG has already a large experience in IEs, the present multi-centric study revealed challenges that still exist related to DNA mixtures interpretation. Overall, the results emphasize the need for further research and training actions in order to improve the analysis of mixtures among the forensic practitioners.

HaploGrep 2: mitochondrial haplogroup classification in the era of high-throughput sequencing

Mitochondrial DNA (mtDNA) profiles can be classified into phylogenetic clusters (haplogroups), which is of great relevance for evolutionary, forensic and medical genetics. With the extensive growth of the underlying phylogenetic tree summarizing the published mtDNA sequences, the manual process of haplogroup classification would be too time-consuming. The previously published classification tool HaploGrep provided an automatic way to address this issue. Here, we present the completely updated version HaploGrep 2 offering several advanced features, including a generic rule-based system for immediate quality control (QC). This allows detecting artificial recombinants and missing variants as well as annotating rare and phantom mutations. Furthermore, the handling of high-throughput data in form of VCF files is now directly supported. For data output, several graphical reports are generated in real time, such as a multiple sequence alignment format, a VCF format and extended haplogroup QC reports, all viewable directly within the application. In addition, HaploGrep 2 generates a publication-ready phylogenetic tree of all input samples encoded relative to the revised Cambridge Reference Sequence. Finally, new distance measures and optimizations of the algorithm increase accuracy and speed-up the application. HaploGrep 2 can be accessed freely and without any registration at http://haplogrep.uibk.ac.at.

The Genomic Legacy of the Transatlantic Slave Trade in the Yungas Valley of Bolivia

During the period of the Transatlantic Slave Trade (TAST) some enslaved Africans were forced to move to Upper Peru (nowadays Bolivia). At first they were sent to Potosí, but later to the tropical Yungas valley where the Spanish colonizers established a so-called “hacienda system” that was based on slave labor, including African-descendants. Due to their isolation, very little attention has been paid so far to ‘Afro-Bolivian’ communities either within the research field of TAST or in genetic population studies. In this study, a total of 105 individuals from the Yungas were sequenced for their mitochondrial DNA (mtDNA) control region, and mitogenomes were obtained for a selected subset of these samples. We also genotyped 46 Ancestry Informative Markers (AIM) in order to investigate continental ancestry at the autosomal level. In addition, Y-chromosome STR and SNP data for a subset of the same individuals was also available from the literature. The data indicate that the partitioning of mtDNA ancestry in the Yungas differs significantly from that in the rest of the country: 81% Native American, 18% African, and 1% European. Interestingly, the great majority of ‘Afro-descendant’ mtDNA haplotypes in the Yungas (84%) concentrates in the locality of Tocaña. This high proportion of African ancestry in the Tocaña is also manifested in the Y-chromosome (44%) and in the autosomes (56%). In sharp contrast with previous studies on the TAST, the ancestry of about 1/3 of the ‘Afro-Bolivian’ mtDNA haplotypes can be traced back to East and South East Africa, which may be at least partially explained by the Arab slave trade connected to the TAST.

Mitogenomes from The 1000 Genome Project reveal new Near Eastern features in present-day Tuscans

BACKGROUND

Genetic analyses have recently been carried out on present-day Tuscans (Central Italy) in order to investigate their presumable recent Near East ancestry in connection with the long-standing debate on the origins of the Etruscan civilization. We retrieved mitogenomes and genome-wide SNP data from 110 Tuscans analyzed within the context of The 1000 Genome Project. For phylogeographic and evolutionary analysis we made use of a large worldwide database of entire mitogenomes (>26,000) and partial control region sequences (>180,000).

RESULTS

Different analyses reveal the presence of typical Near East haplotypes in Tuscans representing isolated members of various mtDNA phylogenetic branches. As a whole, the Near East component in Tuscan mitogenomes can be estimated at about 8%; a proportion that is comparable to previous estimates but significantly lower than admixture estimates obtained from autosomal SNP data (21%). Phylogeographic and evolutionary inter-population comparisons indicate that the main signal of Near Eastern Tuscan mitogenomes comes from Iran.

CONCLUSIONS

Mitogenomes of recent Near East origin in present-day Tuscans do not show local or regional variation. This points to a demographic scenario that is compatible with a recent arrival of Near Easterners to this region in Italy with no founder events or bottlenecks.

Genomic insights on the ethno-history of the Maya and the 'Ladinos' from Guatemala

Background

Guatemala is a multiethnic and multilingual country located in Central America. The main population groups separate ‘Ladinos’ (mixed Native American-African-Spanish), and Native indigenous people of Maya descent. Among the present-day Guatemalan Maya, there are more than 20 different ethnic groups separated by different languages and cultures. Genetic variation of these communities still remains largely unexplored. The principal aim of this study is to explore the genetic variability of the Maya and ‘Ladinos’ from Guatemala by means of uniparental and ancestry informative markers (AIMs).

Results

Analyses of uniparental genetic markers indicate that Maya have a dominant Native American ancestry (mitochondrial DNA [mtDNA]: 100%; Y-chromosome: 94%). ‘Ladino’, however, show a clear gender-bias as indicated by the large European ancestry observed in the Y-chromosome (75%) compared to the mtDNA (0%). Autosomal polymorphisms (AIMs) also mirror this marked gender-bias: (i) Native American ancestry: 92% for the Maya vs. 55% for the ‘Ladino’, and (ii) European ancestry: 8% for the Maya vs. 41% for the ‘Ladino’. In addition, the impact of the Trans-Atlantic slave trade on the present-day Guatemalan population is very low (and only occurs in the ‘Ladino’; mtDNA: 9%; AIMs: 4%), in part mirroring the fact that Guatemala has a predominant orientation to the Pacific Ocean instead of a Caribbean one. Sequencing of entire Guatemalan mitogenomes has led to improved Native American phylogeny via the addition of new haplogroups that are mainly observed in Mesoamerica and/or the North of South America.

Conclusions

The data reveal the existence of a fluid gene flow in the Mesoamerican area and a predominant unidirectional flow towards South America, most likely occurring during the Pre-Classic (1800 BC-200 AD) and the Classic (200–1000 AD) Eras of the Mesoamerican chronology, coinciding with development of the most distinctive and advanced Mesoamerican civilization, the Maya. Phylogenetic features of mtDNA data also suggest a demographic scenario that is compatible with moderate local endogamy and isolation in the Maya combined with episodes of gene exchange between ethnic groups, suggesting an ethno-genesis in the Guatemalan Maya that is recent and supported on a cultural rather than a biological basis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1339-1) contains supplementary material, which is available to authorized users.

DomeTree: a canonical toolkit for mitochondrial DNA analyses in domesticated animals

Mitochondrial DNA (mtDNA) is widely used in various genetic studies of domesticated animals. Many applications require comprehensive knowledge about the phylogeny of mtDNA variants. Herein, we provide the most up-to-date mtDNA phylogeny (i.e. haplogroup tree or matrilineal genealogy) and a standardized hierarchical haplogroup nomenclature system for domesticated cattle, dogs, goats, horses, pigs, sheep, yaks and chickens. These high-resolution mtDNA haplogroup trees based on 1240 complete or near-complete mtDNA genome sequences are available in open resource DomeTree (http://www.dometree.org). In addition, we offer the software MitoToolPy (http://www.mitotool.org/mp.html) to facilitate the mtDNA data analyses. We will continuously and regularly update DomeTree and MitoToolPy.

Phy-Mer: a novel alignment-free and reference-independent mitochondrial haplogroup classifier

MOTIVATION

All current mitochondrial haplogroup classification tools require variants to be detected from an alignment with the reference sequence and to be properly named according to the canonical nomenclature standards for describing mitochondrial variants, before they can be compared with the haplogroup determining polymorphisms. With the emergence of high-throughput sequencing technologies and hence greater availability of mitochondrial genome sequences, there is a strong need for an automated haplogroup classification tool that is alignment-free and agnostic to reference sequence.

RESULTS

We have developed a novel mitochondrial genome haplogroup-defining algorithm using a k-mer approach namely Phy-Mer. Phy-Mer performs equally well as the leading haplogroup classifier, HaploGrep, while avoiding the errors that may occur when preparing variants to required formats and notations. We have further expanded Phy-Mer functionality such that next-generation sequencing data can be used directly as input.

AVAILABILITY AND IMPLEMENTATION

Phy-Mer is publicly available under the GNU Affero General Public License v3.0 on GitHub (https://github.com/danielnavarrogomez/phy-mer).

CONTACT

Xiaowu_Gai@meei.harvard.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genetic portrait of Lisboa immigrant population from Angola with mitochondrial DNA

Portugal has been considered a country of emigrants, nevertheless in the past decades the number of immigrants has grown throughout all the country. This migratory flux has contributed to a raise of heterogeneity at multiple levels. According to statistical data, at the end of 2012 the total number of Angolan immigrants in Portugal equalled about 20,000 individuals. A territorial predominance has been found for the metropolitan region of Lisboa. Angola is a country located in the Atlantic coast of Africa. The presence of Bantu people and the colonisation by Portuguese people on Angolan territory are considered to be the major modulators of the genetic patterns in Angola. Mitochondrial DNA is known for its features that enable an approach to the study of human origin and evolution, as well to the different migration pathways of populations. This genetic marker can also contribute to ascertaining the identity of individuals in forensic cases. The main aim of this study was to determine the genetic structure of the Angolan immigrant population living in Lisboa. Therefore, a total of 173 individuals, inhabitants in Lisboa, nonrelated and with Angolan ancestry were studied. Total control region of mitochondrial DNA was amplified from position 16,024 to position 576 using two pairs of primers - L15997/H016 and L16555/H639. The majority of the identified haplotypes belong to mtDNA lineages known to be specific of the sub-Saharan region. Our results show that this immigrant population inhabitant in Lisboa presents a genetic profile that is characteristic of African populations. This study also demonstrates the genetic diversity that this immigrant population introduces in Lisboa. This does not contradict the historical data concerning colonization of Angola, since this was made mainly by male European individuals, who did not contribute with their maternal information of mtDNA. Lisboa immigrant population from Angola can be accessed via EMPOP dataset with accession number EMPOP662.

Cuba: exploring the history of admixture and the genetic basis of pigmentation using autosomal and uniparental markers

We carried out an admixture analysis of a sample comprising 1,019 individuals from all the provinces of Cuba. We used a panel of 128 autosomal Ancestry Informative Markers (AIMs) to estimate the admixture proportions. We also characterized a number of haplogroup diagnostic markers in the mtDNA and Y-chromosome in order to evaluate admixture using uniparental markers. Finally, we analyzed the association of 16 single nucleotide polymorphisms (SNPs) with quantitative estimates of skin pigmentation. In the total sample, the average European, African and Native American contributions as estimated from autosomal AIMs were 72%, 20% and 8%, respectively. The Eastern provinces of Cuba showed relatively higher African and Native American contributions than the Western provinces. In particular, the highest proportion of African ancestry was observed in the provinces of Guantánamo (40%) and Santiago de Cuba (39%), and the highest proportion of Native American ancestry in Granma (15%), Holguín (12%) and Las Tunas (12%). We found evidence of substantial population stratification in the current Cuban population, emphasizing the need to control for the effects of population stratification in association studies including individuals from Cuba. The results of the analyses of uniparental markers were concordant with those observed in the autosomes. These geographic patterns in admixture proportions are fully consistent with historical and archaeological information. Additionally, we identified a sex-biased pattern in the process of gene flow, with a substantially higher European contribution from the paternal side, and higher Native American and African contributions from the maternal side. This sex-biased contribution was particularly evident for Native American ancestry. Finally, we observed that SNPs located in the genes SLC24A5 and SLC45A2 are strongly associated with melanin levels in the sample.

mitoSAVE: mitochondrial sequence analysis of variants in Excel

The mitochondrial genome (mtGenome) contains genetic information amenable to numerous applications such as medical research, population and evolutionary studies, and human identity testing. However, inconsistent nomenclature assignment makes haplotype comparison difficult and can lead to false exclusion of potentially useful profiles. Massively Parallel Sequencing (MPS) is a platform for sequencing large datasets and potentially whole populations with relative ease. However, the data generated are not easily parsed and interpreted. With this in mind, mitoSAVE has been developed to enable fast conversion of Variant Call Format (VCF) files. mitoSAVE is an Excel-based workbook that converts data within the VCF into mtDNA haplotypes using phylogenetically-established nomenclature as well as rule-based alignments consistent with current forensic standards. mitoSAVE is formatted for human mitochondrial genome; however, it can easily be adapted to support other reasonably small genomes.

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.

Mitochondrial DNA variants in obesity

Heritability estimates for body mass index (BMI) variation are high. For mothers and their offspring higher BMI correlations have been described than for fathers. Variation(s) in the exclusively maternally inherited mitochondrial DNA (mtDNA) might contribute to this parental effect. Thirty-two to 40 mtDNA single nucleotide polymorphisms (SNPs) were available from genome-wide association study SNP arrays (Affymetrix 6.0). For discovery, we analyzed association in a case-control (CC) sample of 1,158 extremely obese children and adolescents and 435 lean adult controls. For independent confirmation, 7,014 population-based adults were analyzed as CC sample of n = 1,697 obese cases (BMI ≥ 30 kg/m2) and n = 2,373 normal weight and lean controls (BMI<25 kg/m2). SNPs were analyzed as single SNPs and haplogroups determined by HaploGrep. Fisher's two-sided exact test was used for association testing. Moreover, the D-loop was re-sequenced (Sanger) in 192 extremely obese children and adolescents and 192 lean adult controls. Association testing of detected variants was performed using Fisher's two-sided exact test. For discovery, nominal association with obesity was found for the frequent allele G of m.8994G/A (rs28358887, p = 0.002) located in ATP6. Haplogroup W was nominally overrepresented in the controls (p = 0.039). These findings could not be confirmed independently. For two of the 252 identified D-loop variants nominal association was detected (m.16292C/T, p = 0.007, m.16189T/C, p = 0.048). Only eight controls carried the m.16292T allele, five of whom belonged to haplogroup W that was initially enriched among these controls. m.16189T/C might create an uninterrupted poly-C tract located near a regulatory element involved in replication of mtDNA. Though follow-up of some D-loop variants still is conceivable, our hypothesis of a contribution of variation in the exclusively maternally inherited mtDNA to the observed larger correlations for BMI between mothers and their offspring could not be substantiated by the findings of the present study.

Reviewing population studies for forensic purposes: Dog mitochondrial DNA

The identification of dog hair through mtDNA analysis has become increasingly important in the last 15 years, as it can provide associative evidence connecting victims and suspects. The evidential value of an mtDNA match between dog hair and its potential donor is determined by the random match probability of the haplotype. This probability is based on the haplotype’s population frequency estimate. Consequently, implementing a population study representative of the population relevant to the forensic case is vital to the correct evaluation of the evidence. This paper reviews numerous published dog mtDNA studies and shows that many of these studies vary widely in sampling strategies and data quality. Therefore, several features influencing the representativeness of a population sample are discussed. Moreover, recommendations are provided on how to set up a dog mtDNA population study and how to decide whether or not to include published data. This review emphasizes the need for improved dog mtDNA population data for forensic purposes, including targeting the entire mitochondrial genome. In particular, the creation of a publicly available database of qualitative dog mtDNA population studies would improve the genetic analysis of dog traces in forensic casework.

The case for the continuing use of the revised Cambridge Reference Sequence (rCRS) and the standardization of notation in human mitochondrial DNA studies

Since the determination in 1981 of the sequence of the human mitochondrial DNA (mtDNA) genome, the Cambridge Reference Sequence (CRS), has been used as the reference sequence to annotate mtDNA in molecular anthropology, forensic science and medical genetics. The CRS was eventually upgraded to the revised version (rCRS) in 1999. This reference sequence is a convenient device for recording mtDNA variation, although it has often been misunderstood as a wild-type (WT) or consensus sequence by medical geneticists. Recently, there has been a proposal to replace the rCRS with the so-called Reconstructed Sapiens Reference Sequence (RSRS). Even if it had been estimated accurately, the RSRS would be a cumbersome substitute for the rCRS, as the new proposal fuses--and thus confuses--the two distinct concepts of ancestral lineage and reference point for human mtDNA. Instead, we prefer to maintain the rCRS and to report mtDNA profiles by employing the hitherto predominant circumfix style. Tree diagrams could display mutations by using either the profile notation (in conventional short forms where appropriate) or in a root-upwards way with two suffixes indicating ancestral and derived nucleotides. This would guard against misunderstandings about reporting mtDNA variation. It is therefore neither necessary nor sensible to change the present reference sequence, the rCRS, in any way. The proposed switch to RSRS would inevitably lead to notational chaos, mistakes and misinterpretations.

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 64mtGenomes (>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.

Indian signatures in the westernmost edge of the European Romani diaspora: new insight from mitogenomes

In agreement with historical documentation, several genetic studies have revealed ancestral links between the European Romani and India. The entire mitochondrial DNA (mtDNA) of 27 Spanish Romani was sequenced in order to shed further light on the origins of this population. The data were analyzed together with a large published dataset (mainly hypervariable region I [HVS-I] haplotypes) of Romani (N=1,353) and non-Romani worldwide populations (N>150,000). Analysis of mitogenomes allowed the characterization of various Romani-specific clades. M5a1b1a1 is the most distinctive European Romani haplogroup; it is present in all Romani groups at variable frequencies (with only sporadic findings in non-Romani) and represents 18% of their mtDNA pool. Its phylogeographic features indicate that M5a1b1a1 originated 1.5 thousand years ago (kya; 95% CI: 1.3-1.8) in a proto-Romani population living in Northwest India. U3 represents the most characteristic Romani haplogroup of European/Near Eastern origin (12.4%); it appears at dissimilar frequencies across the continent (Iberia: ≈ 31%; Eastern/Central Europe: ≈ 13%). All U3 mitogenomes of our Iberian Romani sample fall within a new sub-clade, U3b1c, which can be dated to 0.5 kya (95% CI: 0.3-0.7); therefore, signaling a lower bound for the founder event that followed admixture in Europe/Near East. Other minor European/Near Eastern haplogroups (e.g. H24, H88a) were also assimilated into the Romani by introgression with neighboring populations during their diaspora into Europe; yet some show a differentiation from the phylogenetically closest non-Romani counterpart. The phylogeny of Romani mitogenomes shows clear signatures of low effective population sizes and founder effects. Overall, these results are in good agreement with historical documentation, suggesting that cultural identity and relative isolation have allowed the Romani to preserve a distinctive mtDNA heritage, with some features linking them unequivocally to their ancestral Indian homeland.

A generalized model to estimate the statistical power in mitochondrial disease studies involving 2×k tables

Background

Mitochondrial DNA (mtDNA) variation (i.e. haplogroups) has been analyzed in regards to a number of multifactorial diseases. The statistical power of a case-control study determines the a priori probability to reject the null hypothesis of homogeneity between cases and controls.

Methods/Principal Findings

We critically review previous approaches to the estimation of the statistical power based on the restricted scenario where the number of cases equals the number of controls, and propose a methodology that broadens procedures to more general situations. We developed statistical procedures that consider different disease scenarios, variable sample sizes in cases and controls, and variable number of haplogroups and effect sizes. The results indicate that the statistical power of a particular study can improve substantially by increasing the number of controls with respect to cases. In the opposite direction, the power decreases substantially when testing a growing number of haplogroups. We developed mitPower (http://bioinformatics.cesga.es/mitpower/), a web-based interface that implements the new statistical procedures and allows for the computation of the a priori statistical power in variable scenarios of case-control study designs, or e.g. the number of controls needed to reach fixed effect sizes.

Conclusions/Significance

The present study provides with statistical procedures for the computation of statistical power in common as well as complex case-control study designs involving 2×k tables, with special application (but not exclusive) to mtDNA studies. In order to reach a wide range of researchers, we also provide a friendly web-based tool – mitPower – that can be used in both retrospective and prospective case-control disease studies.

Concept for estimating mitochondrial DNA haplogroups using a maximum likelihood approach (EMMA)

The assignment of haplogroups to mitochondrial DNA haplotypes contributes substantial value for quality control, not only in forensic genetics but also in population and medical genetics. The availability of Phylotree, a widely accepted phylogenetic tree of human mitochondrial DNA lineages, led to the development of several (semi-)automated software solutions for haplogrouping. However, currently existing haplogrouping tools only make use of haplogroup-defining mutations, whereas private mutations (beyond the haplogroup level) can be additionally informative allowing for enhanced haplogroup assignment. This is especially relevant in the case of (partial) control region sequences, which are mainly used in forensics. The present study makes three major contributions toward a more reliable, semi-automated estimation of mitochondrial haplogroups. First, a quality-controlled database consisting of 14,990 full mtGenomes downloaded from GenBank was compiled. Together with Phylotree, these mtGenomes serve as a reference database for haplogroup estimates. Second, the concept of fluctuation rates, i.e. a maximum likelihood estimation of the stability of mutations based on 19,171 full control region haplotypes for which raw lane data is available, is presented. Finally, an algorithm for estimating the haplogroup of an mtDNA sequence based on the combined database of full mtGenomes and Phylotree, which also incorporates the empirically determined fluctuation rates, is brought forward. On the basis of examples from the literature and EMPOP, the algorithm is not only validated, but both the strength of this approach and its utility for quality control of mitochondrial haplotypes is also demonstrated.

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.