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

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.

Genetic Analysis of Mingrelians Reveals Long-Term Continuity of Populations in Western Georgia (Caucasus)

To elucidate the population history of the Caucasus, we conducted a survey of genetic diversity in Samegrelo (Mingrelia), western Georgia. We collected DNA samples and genealogical information from 485 individuals residing in 30 different locations, the vast majority of whom being Mingrelian speaking. From these DNA samples, we generated mitochondrial DNA (mtDNA) control region sequences for all 485 participants (female and male), Y-short tandem repeat haplotypes for the 372 male participants, and analyzed all samples at nearly 590,000 autosomal single nucleotide polymorphisms (SNPs) plus around 33,000 on the sex chromosomes, with 27,000 SNP removed for missingness, using the GenoChip 2.0+ microarray. The resulting data were compared with those from populations from Anatolia, the Caucasus, the Near East, and Europe. Overall, Mingrelians exhibited considerable mtDNA haplogroup diversity, having high frequencies of common West Eurasian haplogroups (H, HV, I, J, K, N1, R1, R2, T, U, and W. X2) and low frequencies of East Eurasian haplogroups (A, C, D, F, and G). From a Y-chromosome standpoint, Mingrelians possessed a variety of haplogroups, including E1b1b, G2a, I2, J1, J2, L, Q, R1a, and R1b. Analysis of autosomal SNP data further revealed that Mingrelians are genetically homogeneous and cluster with other modern-day South Caucasus populations. When compared with ancient DNA samples from Bronze Age archaeological contexts in the broader region, these data indicate that the Mingrelian gene pool began taking its current form at least by this period, probably in conjunction with the formation of a distinct linguistic community.

High Mitochondrial Haplotype Diversity Found in Three Pre-Hispanic Groups from Colombia

The analysis of mitochondrial DNA (mtDNA) hypervariable region (HVR) sequence data from ancient human remains provides valuable insights into the genetic structure and population dynamics of ancient populations. mtDNA is particularly useful in studying ancient populations, because it is maternally inherited and has a higher mutation rate compared to nuclear DNA. To determine the genetic structure of three Colombian pre-Hispanic populations and compare them with current populations, we determined the haplotypes from human bone remains by sequencing several mitochondrial DNA segments. A wide variety of mitochondrial polymorphisms were obtained from 33 samples. Our results support a high population heterogeneity among pre-Hispanic populations in Colombia.

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.

Genetic analysis of the mitochondrial DNA control region in Tai-Kadai-speaking Dong population in southwest China

BACKGROUND

Dong people in Southwest China are officially recognised as an ethnic group, but there has been a lack of population genetic research on this group, especially based on mitochondrial DNA data.

AIM

To study the sequences and haplogroups of the mitochondrial DNA control region in a typical Dong population, and to provide help for the construction of a forensic mitochondrial DNA analysis reference database in East Asia.

SUBJECTS AND METHODS

The sequences of the mitochondrial DNA control region were analysed in 200 individuals of Dong in Guizhou. The haplotype frequencies, haplogroup distribution and paired Fst values of Guizhou Dong and 51 other populations in the world were calculated and explained to explore the genetic polymorphism and population relationships.

RESULTS

A total of 180 haplotypes were detected, with frequencies of 0.005-0.02. All haplotypes were assigned to 97 different haplogroups. The haplotype diversity and random matching probability were 0.998643 and 0.00635, respectively. The paired Fst values and correlation p values of 52 populations showed that the Guizhou Dong had the closest genetic relationship with the Henan Han and the Guizhou Miao in China, and were closest to the Punjab population in Pakistan and the Kashmiri population when compared with the world populations.

CONCLUSIONS

Our study was based on the matrilineal genetic structure of Guizhou Dong to study mitochondrial DNA, which was helpful to promote the establishment of the forensic DNA reference database in East Asia and provide reference for anthropological research.

Mitochondrial Control Region Variants Related to Breast Cancer

Breast cancer has an important incidence in the worldwide female population. Although alterations in the mitochondrial genome probably play an important role in carcinogenesis, the actual evidence is ambiguous and inconclusive. Our purpose was to explore differences in mitochondrial sequences of cases with breast cancer compared with control samples from different origins. We identified 124 mtDNA sequences associated with breast cancer cases, of which 86 were complete and 38 were partial sequences. Of these 86 complete sequences, 52 belonged to patients with a confirmed diagnosis of breast cancer, and 34 sequences were obtained from healthy mammary tissue of the same patients used as controls. From the mtDNA analysis, two polymorphisms with significant statistical differences were found: m.310del (rs869289246) in 34.6% (27/78) of breast cancer cases and 61.7% (21/34) in the controls; and m.315dup (rs369786048) in 60.2% (47/78) of breast cancer cases and 38.2% (13/34) in the controls. In addition, the variant m.16519T>C (rs3937033) was found in 59% of the control sequences and 52% of the breast cancer sequences with a significant statistical difference. Polymorphic changes are evolutionarily related to the haplogroup H of Indo-European and Euro-Asiatic origins; however, they were found in all non-European breast cancers.

Post hoc deconvolution of human mitochondrial DNA mixtures by EMMA 2 using fine-tuned Phylotree nomenclature

•

MtDNA mixtures are observed frequently and difficult to deconvolute.

•

Most previous methods require raw data or quantitative information.

•

EMMA 2 produces valid splittings from consensus sequences of any sequencing technology.

•

EMMA 2 can deconvolute 2 and 3 person mixtures in a fast and traceable way.

In this paper we present a new algorithm for splitting (partial) human mitogenomes into components with high similarity to haplogroup motifs of Phylotree. The algorithm reads a (partial) mitogenome coded by the differences to the reference (rCRS) and outputs the estimated haplogroups of the putative components. The algorithm requires no special information on the raw data of the sequencing process and is therefore suited for the post hoc analysis of mixtures of any sequencing technology. The software EMMA 2 implementing the algorithm will be made available via the EMPOP (https://empop.online) database and extends the nine years old software EMMA for haplogrouping single mitogenomes to mixtures with at most three components.

Evolution and dispersal of mitochondrial DNA haplogroup U5 in Northern Europe: insights from an unsupervised learning approach to phylogeography

Background

We combined an unsupervised learning methodology for analyzing mitogenome sequences with maximum likelihood (ML) phylogenetics to make detailed inferences about the evolution and diversification of mitochondrial DNA (mtDNA) haplogroup U5, which appears at high frequencies in northern Europe.

Methods

Haplogroup U5 mitogenome sequences were gathered from GenBank. The hierarchal Bayesian Analysis of Population Structure (hierBAPS) method was used to generate groups of sequences that were then projected onto a rooted maximum likelihood (ML) phylogenetic tree to visualize the pattern of clustering. The haplogroup statuses of the individual sequences were assessed using Haplogrep2.

Results

A total of 23 hierBAPS groups were identified, all of which corresponded to subclades defined in Phylotree, v.17. The hierBAPS groups projected onto the ML phylogeny accurately clustered all haplotypes belonging to a specific haplogroup in accordance with Haplogrep2. By incorporating the geographic source of each sequence and subclade age estimates into this framework, inferences about the diversification of U5 mtDNAs were made. Haplogroup U5 has been present in northern Europe since the Mesolithic, and spread in both eastern and western directions, undergoing significant diversification within Scandinavia. A review of historical and archeological evidence attests to some of the population interactions contributing to this pattern.

Conclusions

The hierBAPS algorithm accurately grouped mitogenome sequences into subclades in a phylogenetically robust manner. This analysis provided new insights into the phylogeographic structure of haplogroup U5 diversity in northern Europe, revealing a detailed perspective on the diversity of subclades in this region and their distribution in Scandinavian populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08572-y.

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.

Mitochondrial genome sequencing with short overlapping amplicons on MiSeq FGx system

With the development and maturation of massively parallel sequencing (MPS) technology, the mitochondrial genome (mitogenome) sequencing is increasingly applied in the forensic field. In this study, we employed the strategy of short overlapping amplicons for the whole mitogenome, library preparation with tagmentation using the Nextera® XT DNA Library Preparation Kit, sequencing on the MiSeq FGxTM Forensic Genomics System and analyzing data using the mitochondrial(mtDNA) MSR Plug-in and the mtDNA Variant Analyzer. A total of 27 libraries and 56 libraries were sequenced in a run using MiSeq Reagent Kit v2 and v3, respectively. Results showed more than 1800 × of averaged depth of coverage (DoC) at each position. Concordant haplotypes of 9947 A and 2800 M were obtained at 32 variants. Cross-reactivity was observed with 1 ng primate DNA and 10 ng non-primate DNA but could be easily distinguished. Full and accurate variants were obtained from at least 50 pg input DNA and from minor contributors between 19:1 and 1:19 mixed ratios with known reference profiles. More than 86% variants were detected from ≥200-bp degraded samples but its haplotype was assigned to more ancestral haplogroup. Further, a total of 3 962 variants were observed at 613 nucleotide positions from 103 Xibe mitogenomes with 25:1 ratio of transitions to transversions. Two new transversions (C13735A and A14755C) and two tri-alleles at nps 9824 and 16092 were identified. There were 103 unique mitogenome haplotypes from 103 Chinese Xibe that were assigned to 79 haplogroups. Haplogroup D was the preponderant top-level haplogroup in Xibe followed by F, B, M, A, N, G, C, Z, Y, HV and J. Random match probability (RMP) and haplotype diversity (HD) of the whole mitogenome was calculated as 0.0097 and 1.0000, respectively. Compared with HVS-I only, RMP decreased 33.56%, while the number of haplotypes and HD increased 15.73% and 0.49%, respectively. Principal component analysis (PCA) showed that Xibe was clustered to East and Southeast Asian. As a whole, this MPS strategy is suitable for the whole mitogenome sequencing especially for degraded samples and can facilitate generating mitogenome data to support the routine application in forensic sciences. EMP00726 is the first whole mitogenome dataset from Xibe contributed to the EMPOP. Supplemental data for this article are available online at.

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.

A benchmarking of human mitochondrial DNA haplogroup classifiers from whole-genome and whole-exome sequence data

The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. We also assessed the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short-read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

Genetic characterization of a collection of Tsantsas from Ecuadorian museums

Tsantsas are shrunken human heads originally made for ceremonial purposes by Amazonian indigenous groups of the Shuar and Achuar family, previously called Jivaroan tribes. A significant demand of these objects during the first half of the 20th century led to the manufacture of counterfeit shrunken heads for commercial purposes. For museums where these collections are held, as well as for the indigenous groups who claim their ownership, it is important to identify the origin and authenticity of these tsantsas. We hypothesized that a collection of 14 tsantsas from 3 different museum collections in Ecuador are human and aimed to characterize their sex and potential origin. We amplified the amelogenin gene and performed a high resolution melting analysis to determine their human origin and characterize their sex. We also analyzed a fragment (16209-16402) from the HVR-1 region to identify the mtDNA haplogroups present in the tsantsa collection. Our exploratory results show that all the tsantsas are human and that the collection is comprised of 13 males and 1 female. A total of seven mtDNA haplogroups were found among the tsantsa collection using the mtDNA EMPOP database. These results show a predominance of the Amerindian mtDNA haplogroups B, C and D. Additional principal component analysis, genetic distance tree and haplotype network analyses suggest a relationship between the tsantsa specimens and Native American groups.

Sweet tooth: DNA profiling of a cranium from an isolated retained root fragment

Sampling of healthy multi‐rooted teeth is recommended for the genetic identification of human skeletal remains. However, this may not always be possible, as in the reported case consisting of an isolated human cranium found in an aggregate crushing and processing plant in Piedmont, Northwest Italy. The cranium displayed significant weathering, suggesting a post‐mortem interval of several years, and was edentulous with the exception of the apical root fragment of the upper left canine, consequence of an antemortem horizontal fracture. Prolonged decalcification of the root fragment followed by powder‐free DNA extraction from ~10 mg of root tip tissue led to the recovery of >10 ng of high molecular weight human DNA, in comparison with ~0.01 ng of DNA per mg of bone powder obtained from the petrous portion of the temporal bone. Quantity and quality of DNA isolated from apical tooth tissue enabled multiple genotyping, including a reportable female STR profile, mitochondrial DNA analysis, and ancestry‐informative insertion/deletion polymorphisms. Although the cranium remained unidentified after DNA comparisons, our findings confirm that apical tooth tissue is a promising source of DNA, easily obtained through a powder‐free extraction protocol. Results also indicate that root tips should not be overlooked in challenging identification cases, even in the presence of compromised tooth specimens.

HaploGrouper: a generalized approach to haplogroup classification

MOTIVATION

We introduce HaploGrouper, a versatile software to classify haplotypes into haplogroups on the basis of a known phylogenetic tree. A typical use case for this software is the assignment of haplogroups to human mitochondrial DNA (mtDNA) or Y-chromosome haplotypes. Existing state-of-the-art haplogroup-calling software is typically hard-wired to work only with either mtDNA or Y-chromosome haplotypes from humans.

RESULTS

HaploGrouper exhibits comparable accuracy in these instances and has the advantage of being able to assign haplogroups to any kind of haplotypes from any species-given an extant annotated phylogenetic tree defined by sequence variants.

AVAILABILITY AND IMPLEMENTATION

The software is available at the following URL https://gitlab.com/bio_anth_decode/haploGrouper.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

The analysis of ancestry with small-scale forensic panels of genetic markers

In the last 10 years, forensic genetic analysis has been extended beyond identification tests that link a suspect to crime scene evidence using standard DNA profiling, to new supplementary tests that can provide information to investigators about a suspect in the absence of a database hit or eyewitness testimony. These tests now encompass the prediction of physical appearance, ancestry and age. In this review, we give a comprehensive overview of the full range of DNA-based ancestry inference tests designed to work with forensic contact traces, when the level of DNA is often very low or highly degraded. We outline recent developments in the design of ancestry-informative marker sets, forensic assays that use capillary electrophoresis or massively parallel sequencing, and the statistical analysis frameworks that examine the test profile and compares it to reference population variation. Three casework ancestry analysis examples are described which were successfully accomplished in the authors' laboratory, where the ancestry information obtained was critical to the outcome of the DNA analyses made.

Searching the undetected mtDNA variants in forensic MPS data

The efficiency of MPS in forensic mtDNA analysis has been thoroughly proven, although a reliable and well established data evaluation still remains a critical point. Numerous bioinformatics tools have been developed, but most of them require specific operating systems and high costs, while free open-source programs with user-friendly interfaces are few. In this study, 43 full mtGenomes were sequenced using the Ion Personal Genome Machine™ (PGM™) System and analyzed utilizing the plug-in Variant Caller (TVC) of the Ion Torrent Software Suite and the mtDNA-Server (mDS), a free web-based mitochondrial analysis tool for MPS data. The outcomes of these two different analysis tools were compared to variants noted after manual inspection of the aligned reads performed using Integrative Genomics Viewer (IGV). The comparison highlighted the presence of thirty-nine discordant variant calls, which were resolved by Sanger sequencing that confirmed the presence of all variants, except for 7 deletions. The combined adoption of IGV and Sanger type sequencing confirmatory steps, in addition of TVC and mDS analysis, resulted in a more accurate variants assignment with the detection of 32 additional true polymorphisms, which were noted in the final dataset. Regarding the heteroplasmy issue, out of a total of thirty heteroplasmic variants, twenty-eight were detected by the TVC, while the mDS detected twenty-two. Overall, none of the used bioinformatics tools were the perfect choice and a secondary analysis with an expert's opinion in complete mtGenome MPS data evaluation is still required in forensic genetic analysis.

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.

Base specific variation rates at mtDNA positions 16093 and 16183 in human hairs

Small variations between haplotypes detected in different tissues from the same individual have been previously described. These differences complicate the interpretation of mtDNA results in real forensic casework. mtDNA haplotypes detected in hair strands collected at the crime scene have to be frequently compared with haplotypes of reference samples (buccal swabs) from victims or suspects. Nucleotide position 16093 is a well-known hot spot where differences can accumulate between different tissues of the same individual. Intra individual variation was also detected at positions 16182 and 16183 in haplotypes showing an uninterrupted HV1 poly-C stretch (with 16189C). In order to better characterize the type of variation in these positions between buccal cells and hair strands from the same individual, we have performed Sanger sequencing in 25-28 hair strands (411 in total) from 15 individuals showing either an uninterrupted HV1 polyC-stretch (16189C) or 16093C/Y in their buccal cells. The results have been evaluated by also taking into account our previous results published in [19]. We have found that no variation among hair strands was detected in individuals showing T16093 in buccal cells, while variation in hair strands (T16093, 16093C and 16093Y) were detected in individuals showing 16093C or 16093Y in buccal cells. Regarding nucleotide positions 16182 and 16183 in combination with an uninterrupted polyC-stretch, no variation was detected in hairs from individuals showing A16182 16183C in their buccal cells. In contrast, individuals A16182 A16183 showed hair strands with A16182 16183 M and A16182 16183C. And finally, individuals with 16182C 16183C showed some variation in a small amount of their hair strands (some hairs with 16182 M 16183C). These results can be relevant for forensic practitioners when comparing reference samples with hair strands, which is the type of sample most tested by using mtDNA analysis in forensic casework.

Mitochondrial DNA in human identification: a review

Mitochondrial DNA (mtDNA) presents several characteristics useful for forensic studies, especially related to the lack of recombination, to a high copy number, and to matrilineal inheritance. mtDNA typing based on sequences of the control region or full genomic sequences analysis is used to analyze a variety of forensic samples such as old bones, teeth and hair, as well as other biological samples where the DNA content is low. Evaluation and reporting of the results requires careful consideration of biological issues as well as other issues such as nomenclature and reference population databases. In this work we review mitochondrial DNA profiling methods used for human identification and present their use in the main cases of humanidentification focusing on the most relevant issues for forensics.

Big data analysis of human mitochondrial DNA substitution models: a regression approach

BACKGROUND

We study Phylotree, a comprehensive representation of the phylogeny of global human mitochondrial DNA (mtDNA) variations, to better understand the mtDNA substitution mechanism and its most influential factors. We consider a substitution model, where a set of genetic features may predict the rate at which mtDNA substitutions occur. To find an appropriate model, an exhaustive analysis on the effect of multiple factors on the substitution rate is performed through Negative Binomial and Poisson regressions. We examine three different inclusion options for each categorical factor: omission, inclusion as an explanatory variable, and by-value partitioning. The examined factors include genes, codon position, a CpG indicator, directionality, nucleotide, amino acid, codon, and context (neighboring nucleotides), in addition to other site based factors. Partitioning a model by a factor's value results in several sub-models (one for each value), where the likelihoods of the sub-models can be combined to form a score for the entire model. Eventually, the leading models are considered as viable candidates for explaining mtDNA substitution rates.

RESULTS

Initially, we introduce a novel clustering technique on genes, based on three similarity tests between pairs of genes, supporting previous results regarding gene functionalities in the mtDNA. These clusters are then used as a factor in our models. We present leading models for the protein coding genes, rRNA and tRNA genes and the control region, showing it is disadvantageous to separate the models of transitions/transversions, or synonymous/non-synonymous substitutions. We identify a context effect that cannot be attributed solely to protein level constraints or CpG pairs. For protein-coding genes, we show that the substitution model should be partitioned into sub-models according to the codon position and input codon; additionally we confirm that gene identity and cluster have no significant effect once the above factors are accounted for.

CONCLUSIONS

We leverage the large, high-confidence Phylotree mtDNA phylogeny to develop a new statistical approach. We model the substitution rates using regressions, allowing consideration of many factors simultaneously. This admits the use of model selection tools helping to identify the set of factors best explaining the mutational dynamics when considered in tandem.

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.

Mitochondrial genomes reveal an east to west cline of steppe ancestry in Corded Ware populations

From around 4,000 to 2,000 BC the forest-steppe north-western Pontic region was occupied by people who shared a nomadic lifestyle, pastoral economy and barrow burial rituals. It has been shown that these groups, especially those associated with the Yamnaya culture, played an important role in shaping the gene pool of Bronze Age Europeans, which extends into present-day patterns of genetic variation in Europe. Although the genetic impact of these migrations from the forest-steppe Pontic region into central Europe have previously been addressed in several studies, the contribution of mitochondrial lineages to the people associated with the Corded Ware culture in the eastern part of the North European Plain remains contentious. In this study, we present mitochondrial genomes from 23 Late Eneolithic and Bronze Age individuals, including representatives of the north-western Pontic region and the Corded Ware culture from the eastern part of the North European Plain. We identified, for the first time in ancient populations, the rare mitochondrial haplogroup X4 in two Bronze Age Catacomb culture-associated individuals. Genetic similarity analyses show close maternal genetic affinities between populations associated with both eastern and Baltic Corded Ware culture, and the Yamnaya horizon, in contrast to larger genetic differentiation between populations associated with western Corded Ware culture and the Yamnaya horizon. This indicates that females with steppe ancestry contributed to the formation of populations associated with the eastern Corded Ware culture while more local people, likely of Neolithic farmer ancestry, contributed to the formation of populations associated with western Corded Ware culture.

Mitochondrial DNA control region diversity in a population from Parana state-increasing the Brazilian forensic database

The entire mitochondrial DNA (mtDNA) control region (nucleotide position 16024-576) sequences were obtained through Sanger sequencing method for 122 individuals from Parana state, South of Brazil. We observed a total of 108 different haplotypes of which 97 were unique and 11 were shared by more than one individual. The haplogroups were classified according to the updated mtDNA phylogeny, by EMMA (estimating mitochondrial haplogroups using a maximum likelihood approach). Our results revealed the predominance of Amerindian haplogroups with a frequency of 49.2% of the population sample, followed by European lineages with 38.5% and 12.3% of African lineages. Parana population sample set presented a high haplotype diversity (0.9976) and the random match probability was 0.0106. The phylogenetical findings and the diversity indices confirm the high genetic heterogeneity of this population and suggest a high informativeness of mtDNA analyses in forensic cases. The population data will contribute to increase the Brazilian mtDNA database for forensic purposes and it is available through EMPOP (European DNA Profiling Group mitochondrial DNA population database) under the accession number EMP00714.

Validation of haplotype-specific extraction for separating a mitochondrial DNA model mixture and application to simulated casework

Haplotype-specific extraction (HSE) is a new field of application for the separation of mitochondrial DNA (mtDNA) mixtures and is developed to identify the mtDNA haplotypes of the contributors subsequently by sequencing. Here we show the validation of HSE with an exemplary mitochondrial DNA mixture into its individual haplotypes according to our laboratory standards. These specify several critical areas of assay performance to be tested, such as sensitivity, robustness and mixture studies comprising varying proportions of their components,degraded samples and samples of different qualities and material. Wereport the successful and unambiguous analysis of the exemplary separated mitochondrial DNA mixture under various conditions as well as simulated casework samples, which manifest as mixed nucleotide calls at single base positions previously. Here we demonstrate that the HSE assay is high sensitive, stable against degradation and applicable in a wide range of sample qualities. Based on our findings from the validation study, we believe that this assay has great potential power and may be useful for distinguishing among the mtDNA of individuals and their geographical origin in mixed DNA samples.

Concordance of mitochondrial DNA sequencing methods on bloodstains using Ion PGM™

In this study, the complete mitochondrial genome (mtGenome) of six samples from three forensic cases was sequenced using the Ion Torrent Personal Genome Machine (PGM). The analyzed samples from forensic cases included bloodstains from several materials, such as gauze, Flinder's Technology Associates (FTA) cards and swabs. The age of the samples ranged from two months to twelve years. The complete mtGenomes were amplified using the tiling sequencing strategy which divided the whole mtGenome into 162 amplicons. All amplicons were successfully recovered. A phylogenetic analysis was performed to determine the accuracy of the PGM data, and which were compared to partial Sanger-based sequencing data. The average coverage of the PGM data were above 4000× in all case samples, and 99.86% concordance was observed using both sequencing methods. In conclusion, we demonstrate the ability to recover the complete mtGenome from bloodstains with relatively poor DNA quality by PGM. Moreover, the results are concordant with Sanger sequencing data. This new method has potential use in forensic practice.

Developmental validation of a Nextera XT mitogenome Illumina MiSeq sequencing method for high-quality samples

Generating mitochondrial genome (mitogenome) data from reference samples in a rapid and efficient manner is critical to harnessing the greater power of discrimination of the entire mitochondrial DNA (mtDNA) marker. The method of long-range target enrichment, Nextera XT library preparation, and Illumina sequencing on the MiSeq is a well-established technique for generating mitogenome data from high-quality samples. To this end, a validation was conducted for this mitogenome method processing up to 24 samples simultaneously along with analysis in the CLC Genomics Workbench and utilizing the AQME (AFDIL-QIAGEN mtDNA Expert) tool to generate forensic profiles. This validation followed the Federal Bureau of Investigation's Quality Assurance Standards (QAS) for forensic DNA testing laboratories and the Scientific Working Group on DNA Analysis Methods (SWGDAM) validation guidelines. The evaluation of control DNA, non-probative samples, blank controls, mixtures, and nonhuman samples demonstrated the validity of this method. Specifically, the sensitivity was established at ≥25 pg of nuclear DNA input for accurate mitogenome profile generation. Unreproducible low-level variants were observed in samples with low amplicon yields. Further, variant quality was shown to be a useful metric for identifying sequencing error and crosstalk. Success of this method was demonstrated with a variety of reference sample substrates and extract types. These studies further demonstrate the advantages of using NGS techniques by highlighting the quantitative nature of heteroplasmy detection. The results presented herein from more than 175 samples processed in ten sequencing runs, show this mitogenome sequencing method and analysis strategy to be valid for the generation of reference data.

A review of bioinformatic methods for forensic DNA analyses

Short tandem repeats, single nucleotide polymorphisms, and whole mitochondrial analyses are three classes of markers which will play an important role in the future of forensic DNA typing. The arrival of massively parallel sequencing platforms in forensic science reveals new information such as insights into the complexity and variability of the markers that were previously unseen, along with amounts of data too immense for analyses by manual means. Along with the sequencing chemistries employed, bioinformatic methods are required to process and interpret this new and extensive data. As more is learnt about the use of these new technologies for forensic applications, development and standardization of efficient, favourable tools for each stage of data processing is being carried out, and faster, more accurate methods that improve on the original approaches have been developed. As forensic laboratories search for the optimal pipeline of tools, sequencer manufacturers have incorporated pipelines into sequencer software to make analyses convenient. This review explores the current state of bioinformatic methods and tools used for the analyses of forensic markers sequenced on the massively parallel sequencing (MPS) platforms currently most widely used.

Optimized mtDNA Control Region Primer Extension Capture Analysis for Forensically Relevant Samples and Highly Compromised mtDNA of Different Age and Origin

The analysis of mitochondrial DNA (mtDNA) has proven useful in forensic genetics and ancient DNA (aDNA) studies, where specimens are often highly compromised and DNA quality and quantity are low. In forensic genetics, the mtDNA control region (CR) is commonly sequenced using established Sanger-type Sequencing (STS) protocols involving fragment sizes down to approximately 150 base pairs (bp). Recent developments include Massively Parallel Sequencing (MPS) of (multiplex) PCR-generated libraries using the same amplicon sizes. Molecular genetic studies on archaeological remains that harbor more degraded aDNA have pioneered alternative approaches to target mtDNA, such as capture hybridization and primer extension capture (PEC) methods followed by MPS. These assays target smaller mtDNA fragment sizes (down to 50 bp or less), and have proven to be substantially more successful in obtaining useful mtDNA sequences from these samples compared to electrophoretic methods. Here, we present the modification and optimization of a PEC method, earlier developed for sequencing the Neanderthal mitochondrial genome, with forensic applications in mind. Our approach was designed for a more sensitive enrichment of the mtDNA CR in a single tube assay and short laboratory turnaround times, thus complying with forensic practices. We characterized the method using sheared, high quantity mtDNA (six samples), and tested challenging forensic samples (n = 2) as well as compromised solid tissue samples (n = 15) up to 8 kyrs of age. The PEC MPS method produced reliable and plausible mtDNA haplotypes that were useful in the forensic context. It yielded plausible data in samples that did not provide results with STS and other MPS techniques. We addressed the issue of contamination by including four generations of negative controls, and discuss the results in the forensic context. We finally offer perspectives for future research to enable the validation and accreditation of the PEC MPS method for final implementation in forensic genetic laboratories.

Mitogenomic diversity in Russians and Poles

Complete mtDNA genome sequencing improves molecular resolution for distinguishing variation between individuals and populations, but there is still deficiency of mitogenomic population data. To overcome this limitation, we used Sanger-based protocol to generate complete mtDNA sequences of 376 Russian individuals from six populations of European part of Russia and 100 Polish individuals from northern Poland. Nearly complete resolution of mtDNA haplotypes was achieved - about 97% of haplotypes were unique both in Russians and Poles, and no haplotypes overlapped between them when indels were considered. While European populations showed a low, but statistically significant level of between-population differentiation (Fst=0.66%, p=0), Russians demonstrate lack of between-population differences (Fst=0.22%, p=0.15). Results of the Bayesian skyline analysis of Russian mitogenomes demonstrate not only post-Last Glacial Maximum expansion, but also rapid population growth starting from about 4.3kya (95% CI: 2.9-5.8kya), i.e. in the Bronze Age. This expansion strongly correlates with the Kurgan model established by archaeologists and confirmed by paleogeneticists.

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.

Analysis of human mitochondrial DNA sequences from fecally polluted environmental waters as a tool to study population diversity

Mitochondrial signature sequences have frequently been used to study human population diversity around the world. Traditionally, this requires obtaining samples directly from individuals which is cumbersome, time consuming and limited to the number of individuals that participated in these types of surveys. Here, we used environmental DNA extracts to determine the presence and sequence variability of human mitochondrial sequences as a means to study the diversity of populations inhabiting in areas nearby a tropical watershed impacted with human fecal pollution. We used high-throughput sequencing (Illumina) and barcoding to obtain thousands of sequences from the mitochondrial hypervariable region 2 (HVR2) and determined the different haplotypes present in 10 different water samples. Sequence analyses indicated a total of 19 distinct variants with frequency greater than 5%. The HVR2 sequences were associated with haplogroups of West Eurasian (57.6%), Sub-Saharan African (23.9%), and American Indian (11%) ancestry. This was in relative accordance with population census data from the watershed sites. The results from this study demonstrates the potential value of mitochondrial sequence data retrieved from fecally impacted environmental waters to study the population diversity of local municipalities. This environmental DNA approach may also have other public health implications such as tracking background levels of human mitochondrial genes associated with diseases. It may be possible to expand this approach to other animal species inhabiting or using natural water systems.

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.

Considering DNA damage when interpreting mtDNA heteroplasmy in deep sequencing data

Resolution of mitochondrial (mt) DNA heteroplasmy is now possible when applying a massively parallel sequencing (MPS) approach, including minor components down to 1%. However, reporting thresholds and interpretation criteria will need to be established for calling heteroplasmic variants that address a number of important topics, one of which is DNA damage. We assessed the impact of increasing amounts of DNA damage on the interpretation of minor component sequence variants in the mtDNA control region, including low-level mixed sites. A passive approach was used to evaluate the impact of storage conditions, and an active approach was employed to accelerate the process of hydrolytic damage (for example, replication errors associated with depurination events). The patterns of damage were compared and assessed in relation to damage typically encountered in poor quality samples. As expected, the number of miscoding lesions increased as conditions worsened. Single nucleotide polymorphisms (SNPs) associated with miscoding lesions were indistinguishable from innate heteroplasmy and were most often observed as 1-2% of the total sequencing reads. Numerous examples of miscoding lesions above 2% were identified, including two complete changes in the nucleotide sequence, presenting a challenge when assessing the placement of reporting thresholds for heteroplasmy. To mitigate the impact, replication of miscoding lesions was not observed in stored samples, and was rarely seen in data associated with accelerated hydrolysis. In addition, a significant decrease in the expected transition:transversion ratio was observed, providing a useful tool for predicting the presence of damage-induced lesions. The results of this study directly impact MPS analysis of minor sequence variants from poorly preserved DNA extracts, and when biological samples have been exposed to agents that induce DNA damage. These findings are particularly relevant to clinical and forensic investigations.

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.

Hairy matters: MtDNA quantity and sequence variation along and among human head hairs

Hairs from the same donor have been found to differ in mtDNA sequence within and among themselves and from other tissues, which impacts interpretation of results obtained in a forensic setting. However, little is known on the magnitude of this phenomenon and published data on systematic studies are scarce. We addressed this issue by generating mtDNA control region (CR) profiles of >450 hair fragments from 21 donors by Sanger-type sequencing (STS). To mirror forensic scenarios, we compared hair haplotypes from the same donors to each other, to the corresponding buccal swab reference haplotypes and analyzed several fragments of individual hairs. We also investigated the effects of hair color, donor sex and age, mtDNA haplogroup and chemical treatment on mtDNA quantity, amplification success and variation. We observed a wide range of individual CR sequence variation. The reference haplotype was the only or most common (≥75%) hair haplotype for most donors. However, in two individuals, the reference haplotype was only found in about a third of the investigated hairs, mainly due to differences at highly variable positions. Similarly, most hairs revealed the reference haplotype along their entire length, however, about a fifth of the hairs contained up to 71% of segments with deviant haplotypes, independent of the longitudinal position. Variation affected numerous positions, typically restricted to the individual hair and in most cases heteroplasmic, but also fixed (i.e. homoplasmic) substitutions were observed. While existing forensic mtDNA interpretation guidelines were found still sufficient for all comparisons to reference haplotypes, some comparisons between hairs from the same donor could yield false exclusions when those guidelines are strictly followed. This study pinpoints the special care required when interpreting mtDNA results from hair in forensic casework.

Haplotype diversity in mitochondrial genome in a Chinese Han population

Investigations into the use of mitochondrial genome (mtGenome) typing by massively parallel sequencing technologies are well underway in many areas, including forensic genetics. Previous studies have demonstrated that mtGenome sequencing data generated from Ion torrent personal genome machine (PGM) system were highly viable and reliable in forensic research. In this study, 145 whole mtGenomes from unrelated Chinese Han population were sequenced using the Ion PGM system. Results showed that 145 distinct haplotypes were obtained at a relatively high coverage with limited strand bias. The distribution of variants across the entire mtGenomes was illustrated and 70.74% of the variants were observed outside of the control region. An overall increase in the number of unique haplotypes as well as haplotype diversity were observed by detection of mtGenome compared with hypervariable region I/II (HV I/II) and control region (CR). This study demonstrates the substantially higher degree of haplotype resolution with whole-mtGenome sequences in comparison to HV I/II or CR that historically targeted for forensic testing, which shows the potential value of mtGenome typing in forensic testing in the future.

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.

Maternal Genetic Composition of a Medieval Population from a Hungarian-Slavic Contact Zone in Central Europe

The genetic composition of the medieval populations of Central Europe has been poorly investigated to date. In particular, the region of modern-day Slovakia is a blank spot in archaeogenetic research. This paper reports the study of mitochondrial DNA (mtDNA) in ancient samples from the 9^th–12^th centuries originating from the cemeteries discovered in Nitra-Šindolka and Čakajovce, located in western Slovakia (Central Europe). This geographical region is interesting to study because its medieval multi-ethnic population lived in the so-called contact zone of the territory of the Great Moravian and later Hungarian state formations. We described 16 different mtDNA haplotypes in 19 individuals, which belong to the most widespread European mtDNA haplogroups: H, J, T, U and R0. Using comparative statistical and population genetic analyses, we showed the differentiation of the European gene pool in the medieval period. We also demonstrated the heterogeneous genetic characteristics of the investigated population and its affinity to the populations of modern Europe.