Mitochondrial DNA heteroplasmy in the emerging field of massively parallel sequencing

Novel blood and tissue-based mitochondrial D-loop mutations detected in an Iranian NAFLD patient cohort

Non-alcoholic fatty liver disease (NAFLD) is an increasingly prevalent chronic liver disease characterized by an elusive etiology. In its advanced stages, this condition can pose life-threatening implications. Mitochondrial dysfunction due to its impact on hepatic lipid homeostasis, cytokine release, ROS production, and cell death, contributes to the pathogenesis of NAFLD. Previous research reveals a direct link between NAFLD genetic predictors and mitochondrial dysfunction. The emphasis on the D-loop stems from its association with impaired mtDNA replication, underscoring its crucial role in NAFLD progression. We included 38 Iranian NAFLD patients (comprising 16 patients with non-alcoholic fatty liver [NAFL] and 22 patients with non-alcoholic steatohepatitis [NASH]), with matched blood and liver tissue samples collected from each to compare variations in the mitochondrial D-loop sequence within samples. The mitochondrial DNA (mtDNA) D-loop region was amplified using PCR, and variations were identified through sequencing. The resultant sequences were compared with the reference sequence of human mtDNA available in the MITOMAP Database for comparative analysis. In this study, 97 somatic mutations in the mtDNA D-loop region were identified in NAFLD patients. Our study revealed significant difference between the NAFLD patients and control group in 13 detected mutations (P ≤ 0.05). Novel mutations were discovered in hepatic tissues, while mutation 16220-16221ins C was found in both tissues and blood. A significant difference was found in the distribution of D310 and mt514-mt523 (CA)n repeat variations between NAFLD patients and the control group (P < 0.001). C to T and T to C transitions were the prevalent substitution among patients. Identification of the 16220-16221ins C mutation in both blood and tissue samples from NAFLD patients holds substantial promise as a potential diagnostic marker. However, further research is imperative to corroborate these findings.

A PCR-independent approach for mtDNA enrichment and next-generation sequencing: comprehensive evaluation and clinical application

BACKGROUND

Sequencing the mitochondrial genome has been increasingly important for the investigation of primary mitochondrial diseases (PMD) and mitochondrial genetics. To overcome the limitations originating from PCR-based mtDNA enrichment, we set out to develop and evaluate a PCR-independent approach in this study, named Pime-Seq (PCR-independent mtDNA enrichment and next generation Sequencing).

RESULTS

By using the optimized mtDNA enrichment procedure, the mtDNA reads ratio reached 88.0 ± 7.9% in the sequencing library when applied on human PBMC samples. We found the variants called by Pime-Seq were highly consistent among technical repeats. To evaluate the accuracy and reliability of this method, we compared Pime-Seq with lrPCR based NGS by performing both methods simultaneously on 45 samples, yielding 1677 concordant variants, as well as 146 discordant variants with low-level heteroplasmic fraction, in which Pime-Seq showed higher reliability. Furthermore, we applied Pime-Seq on 4 samples of PMD patients retrospectively, and successfully detected all the pathogenic mtDNA variants. In addition, we performed a prospective study on 192 apparently healthy pregnant women during prenatal screening, in which Pime-Seq identified pathogenic mtDNA variants in 4 samples, providing extra information for better health monitoring in these cases.

CONCLUSIONS

Pime-Seq can obtain highly enriched mtDNA in a PCR-independent manner for high quality and reliable mtDNA deep-sequencing, which provides us an effective and promising tool for detecting mtDNA variants for both clinical and research purposes.

MitoH3: Mitochondrial Haplogroup and Homoplasmic/Heteroplasmic Variant Calling Pipeline for Alzheimer's Disease Sequencing Project

Background

Mitochondrial DNA (mtDNA) is a double-stranded circular DNA and has multiple copies in each cell. Excess heteroplasmy, the coexistence of distinct variants in copies of mtDNA within a cell, may lead to mitochondrial impairments. Accurate determination of heteroplasmy in whole-genome sequencing (WGS) data has posed a significant challenge because mitochondria carrying heteroplasmic variants cannot be distinguished during library preparation. Moreover, sequencing errors, contamination, and nuclear mtDNA segments can reduce the accuracy of heteroplasmic variant calling.

Objective

To efficiently and accurately call mtDNA homoplasmic and heteroplasmic variants from the large-scale WGS data generated from the Alzheimer's Disease Sequencing Project (ADSP), and test their association with Alzheimer's disease (AD).

Methods

In this study, we present MitoH3-a comprehensive computational pipeline for calling mtDNA homoplasmic and heteroplasmic variants and inferring haplogroups in the ADSP WGS data. We first applied MitoH3 to 45 technical replicates from 6 subjects to define a threshold for detecting heteroplasmic variants. Then using the threshold of 5% ≤variant allele fraction≤95%, we further applied MitoH3 to call heteroplasmic variants from a total of 16,113 DNA samples with 6,742 samples from cognitively normal controls and 6,183 from AD cases.

Results

This pipeline is available through the Singularity container engine. For 4,311 heteroplasmic variants identified from 16,113 samples, no significant variant count difference was observed between AD cases and controls.

Conclusions

Our streamlined pipeline, MitoH3, enables computationally efficient and accurate analysis of a large number of samples.

Systematic evaluation of the Precision ID GlobalFiler™ NGS STR panel v2 using single-source samples of various quantity and quality and mixed DNA samples

Massively parallel sequencing (MPS) techniques were developed approximately 15 years ago. Meanwhile, several MPS kits for forensic identification, phenotypic information, ancestry, and mitochondrial DNA analysis have been developed and their use has been established. Sequencing short tandem repeats (STRs) has certain advantages over the currently used length-based genotyping methods, which are based on PCR amplification followed by capillary electrophoresis (CE). MPS is more discriminative and includes the possibility of testing high numbers of targets (> 100), different types of markers [STRs and single nucleotide polymorphisms (SNPs)], as well as the use of smaller amplicons (< 300 bp). This study evaluated in 24 experimental runs the Precision ID GlobalFiler™ NGS STR panel v2 from ThermoFisher, which targets 31 autosomal STRs, amelogenin, and three Y-markers (one STR, SRY, and Yindel). Single-source samples were used in 18 experimental runs, for systematic evaluation. These included assessing library preparation benchmark conditions, limited DNA input, as well as testing repeatability, number of samples per run, and degraded DNA samples. Full profiles were consistently obtained from as little as 50 pg DNA input. Using the optional recovery PCR method improved outcomes for samples with low DNA input. Full profiles were also obtained from severely degraded DNA samples with degradation indices (DI) of > 60. In addition, six experimental runs were performed testing various two-person mixtures with mixture ratios ranging from 1:20 to 20:1. Major and minor contributors were distinguishable by their read counts (coverage), because less DNA input yielded lower read counts, analogous to the traditional CE technology, where less DNA produces lower peak heights. Mixture ratios of approximately 1:1 were indistinguishable, while a greater imbalance, i.e., higher mixture ratios, made the mixture more distinguishable between major and minor contributors. Based on this information, the highest success rate of correctly deconvoluted four-allelic loci was from mixtures with 1:3 ratios. At higher mixture ratios, the drop-out rate of the minor contributor increased, reducing the number of four-allelic loci.

Targeted enrichment of whole-genome SNPs from highly burned skeletal remains

Genetic assessment of highly incinerated and/or degraded human skeletal material is a persistent challenge in forensic DNA analysis, including identifying victims of mass disasters. Few studies have investigated the impact of thermal degradation on whole-genome single-nucleotide polymorphism (SNP) quality and quantity using next-generation sequencing (NGS). We present whole-genome SNP data obtained from the bones and teeth of 27 fire victims using two DNA extraction techniques. Extracts were converted to double-stranded DNA libraries then enriched for whole-genome SNPs using unpublished biotinylated RNA baits and sequenced on an Illumina NextSeq 550 platform. Raw reads were processed using the EAGER (Efficient Ancient Genome Reconstruction) pipeline, and the SNPs filtered and called using FreeBayes and GATK (v. 3.8). Mixed-effects modeling of the data suggest that SNP variability and preservation is predominantly determined by skeletal element and burn category, and not by extraction type. Whole-genome SNP data suggest that selecting long bones, hand and foot bones, and teeth subjected to temperatures <350°C are the most likely sources for higher genomic DNA yields. Furthermore, we observed an inverse correlation between the number of captured SNPs and the extent to which samples were burned, as well as a significant decrease in the total number of SNPs measured for samples subjected to temperatures >350°C. Our data complement previous analyses of burned human remains that compare extraction methods for downstream forensic applications and support the idea of adopting a modified Dabney extraction technique when traditional forensic methods fail to produce DNA yields sufficient for genetic identification.

Mitochondrial point heteroplasmy: insights from deep-sequencing of human replicate samples

BACKGROUND

Human mitochondrial heteroplasmy is an extensively investigated phenomenon in the context of medical diagnostics, forensic identification and molecular evolution. However, technical limitations of high-throughput sequencing hinder reliable determination of point heteroplasmies (PHPs) with minor allele frequencies (MAFs) within the noise threshold.

RESULTS

To investigate the PHP landscape at an MAF threshold down to 0.1%, we sequenced whole mitochondrial genomes at approximately 7.700x coverage, in multiple technical and biological replicates of longitudinal blood and buccal swab samples from 11 human donors (159 libraries in total). The results obtained by two independent sequencing platforms and bioinformatics pipelines indicate distinctive PHP patterns below and above the 1% MAF cut-off. We found a high inter-individual prevalence of low-level PHPs (MAF < 1%) at polymorphic positions of the mitochondrial DNA control region (CR), their tissue preference, and a tissue-specific minor allele linkage. We also established the position-dependent potential of minor allele expansion in PHPs, and short-term PHP instability in a mitotically active tissue. We demonstrate that the increase in sensitivity of PHP detection to minor allele frequencies below 1% within a robust experimental and analytical pipeline, provides new information with potential applicative value.

CONCLUSIONS

Our findings reliably show different mutational loads between tissues at sub-1% allele frequencies, which may serve as an informative medical biomarker of time-dependent, tissue-specific mutational burden, or help discriminate forensically relevant tissues in a single person, close maternal relatives or unrelated individuals of similar phylogenetic background.

Heteroplasmy and Individual Mitogene Pools: Characteristics and Potential Roles in Ecological Studies

The mitochondrial genome (mitogenome or mtDNA), the extrachromosomal genome, is a multicopy circular DNA with high mutation rates due to replication and repair errors. A mitochondrion, cell, tissue, organ, or an individual body may hold multiple variants, both inherited and developed over a lifetime, which make up individual mitogene pools. This phenomenon is also called mtDNA heteroplasmy. MtDNA variants influence cellular and tissular functions and are consequently subjected to selection. Although it has long been recognized that only inheritable germline heteroplasmies have evolutionary significance, non-inheritable somatic heteroplasmies have been overlooked since they directly affect individual fitness and thus indirectly affect the fate of heritable germline variants. This review focuses on the characteristics, dynamics, and functions of mtDNA heteroplasmy and proposes the concept of individual mitogene pools to discuss individual genetic diversity from multiple angles. We provide a unique perspective on the relationship between individual genetic diversity and heritable genetic diversity and guide how the individual mitogene pool with novel genetic markers can be applied to ecological research.

Review: Computational analysis of human skeletal remains in ancient DNA and forensic genetics

Degraded DNA is used to answer questions in the fields of ancient DNA (aDNA) and forensic genetics. While aDNA studies typically center around human evolution and past history, and forensic genetics is often more concerned with identifying a specific individual, scientists in both fields face similar challenges. The overlap in source material has prompted periodic discussions and studies on the advantages of collaboration between fields toward mutually beneficial methodological advancements. However, most have been centered around wet laboratory methods (sampling, DNA extraction, library preparation, etc.). In this review, we focus on the computational side of the analytical workflow. We discuss limitations and considerations to consider when working with degraded DNA. We hope this review provides a framework to researchers new to computational workflows for how to think about analyzing highly degraded DNA and prompts an increase of collaboration between the forensic genetics and aDNA fields.

Whole mitochondrial genome analysis in highland Tibetans: further matrilineal genetic structure exploration

Introduction: The Qinghai-Tibet Plateau is one of the last terrestrial environments conquered by modern humans. Tibetans are among the few high-altitude settlers in the world, and understanding the genetic profile of Tibetans plays a pivotal role in studies of anthropology, genetics, and archaeology. Methods: In this study, we investigated the maternal genetic landscape of Tibetans based on the whole mitochondrial genome collected from 145 unrelated native Lhasa Tibetans. Molecular diversity indices, haplotype diversity (HD), Tajima's D and Fu's Fs were calculated and the Bayesian Skyline Plot was obtained to determining the genetic profile and population fluctuation of Lhasa Tibetans. To further explore the genetic structure of Lhasa Tibetans, we collected 107 East Asian reference populations to perform principal component analysis (PCA), multidimensional scaling (MDS), calculated Fst values and constructed phylogenetic tree. Results: The maternal genetic landscape of Tibetans showed obvious East Asian characteristics, M9a (28.28%), R (11.03%), F1 (12.41%), D4 (9.66%), N (6.21%), and M62 (4.14%) were the dominant haplogroups. The results of PCA, MDS, Fst and phylogenetic tree were consistent: Lhasa Tibetans clustered with other highland Tibeto-Burman speakers, there was obvious genetic homogeneity of Tibetans in Xizang, and genetic similarity between Tibetans and northern Han people and geographically adjacent populations was found. In addition, specific maternal lineages of Tibetans also be determined in this study. Discussion: In general, this study further shed light on long-time matrilineal continuity on the Tibetan Plateau and the genetic connection between Tibetans and millet famers in the Yellow River Basin, and further revealed that multiple waves of population interaction and admixture during different historical periods between lowland and highland populations shaped the maternal genetic profile of Tibetans.

Complete Mitochondrial DNA Genome Variation in the Swedish Population

The development of complete mitochondrial genome (mitogenome) reference data for inclusion in publicly available population databases is currently underway, and the generation of more high-quality mitogenomes will only enhance the statistical power of this forensically useful locus. To characterize mitogenome variation in Sweden, the mitochondrial DNA (mtDNA) reads from the SweGen whole genome sequencing (WGS) dataset were analyzed. To overcome the interference from low-frequency nuclear mtDNA segments (NUMTs), a 10% variant frequency threshold was applied for the analysis. In total, 934 forensic-quality mitogenome haplotypes were characterized. Almost 45% of the SweGen haplotypes belonged to haplogroup H. Nearly all mitogenome haplotypes (99.1%) were assigned to European haplogroups, which was expected based on previous mtDNA studies of the Swedish population. There were signature northern Swedish and Finnish haplogroups observed in the dataset (e.g., U5b1, W1a), consistent with the nuclear DNA analyses of the SweGen data. The complete mitogenome analysis resulted in high haplotype diversity (0.9996) with a random match probability of 0.15%. Overall, the SweGen mitogenomes provide a large mtDNA reference dataset for the Swedish population and also contribute to the effort to estimate global mitogenome haplotype frequencies.

Quantifying human genome parameters in aging

Healthy human longevity is a global goal of the world health system. Determining the causes and processes influencing human longevity is the primary fundamental goal facing the scientific community. Currently, the main efforts of the scientific community are aimed at identifying the qualitative characteristics of the genome that determine the trait. At the same time, when evaluating qualitative characteristics, there are many challenges that make it difficult to establish associations. Quantitative traits are burdened with such problems to a lesser extent, but they are largely overlooked in current genomic studies of aging and longevity. Although there is a wide repertoire of quantitative trait analyses based on genomic data, most opportunities are ignored by authors, which, along with the inaccessibility of published data, leads to the loss of this important information. This review focuses on describing quantitative traits important for understanding aging and necessary for analysis in further genomic studies, and recommends the inclusion of the described traits in the analysis. The review considers the relationship between quantitative characteristics of the mitochondrial genome and aging, longevity, and age-related neurodegenerative diseases, such as the frequency of extensive mitochondrial DNA (mtDNA) deletions, mtDNA half-life, the frequency of A>G replacements in the mtDNA heavy chain, the number of mtDNA copies; special attention is paid to the mtDNA methylation sign. A separate section of this review is devoted to the correlation of telomere length parameters with age, as well as the association of telomere length with the amount of mitochondrial DNA. In addition, we consider such a quantitative feature as the rate of accumulation of somatic mutations with aging in relation to the lifespan of living organisms. In general, it may be noted that there are quite serious reasons to suppose that various quantitative characteristics of the genome may be directly or indirectly associated with certain aspects of aging and longevity. At the same time, the available data are clearly insufficient for definitive conclusions and the determination of causal relationships.

Population analysis of complete mitogenomes for 334 samples from El Salvador

The use of mitochondrial DNA (mtDNA) in the field of forensic genetics is widely spread mainly due to its advantages when identifying highly degraded samples. In this sense, massive parallel sequencing has made the analysis of the whole mitogenome more accessible, noticeably increasing the informativeness of mtDNA haplotypes. The civil war (1980-1992) in El Salvador caused many deaths and disappearances (including children) all across the country and the economic and social instability after the war forced many people to emigration. For this reason, different organizations have collected DNA samples from relatives with the aim of identifying missing people. Thus, we present a dataset containing 334 complete mitogenomes from the Salvadoran general population. To the best of our knowledge, this is the first publication of a nationwide forensic-quality complete mitogenome database of any Latin American country. We found 293 different haplotypes, with a random match probability of 0.0041 and 26.6 mean pairwise differences, which is similar to other Latin American populations, and which represent a marked improvement from the values obtained with just control region sequences. These haplotypes belong to 54 different haplogroups, being 91% of them of Native American origin. Over a third (35.9%) of the individuals carried at least a heteroplasmic site (excluding length heteroplasmies). Ultimately, the present database aims to represent mtDNA haplotype diversity in the general Salvadoran populations as a basis for the identification of people that disappeared during or after the civil war.

Exploring rare differences in mitochondrial genome between MZ twins using Ion Torrent semiconductor sequencing

Monozygotic (MZ) twins are considered to be genetically identical in that they have the same genomic DNA sequences in theory, and thus cannot be differentiated using forensic standard STR-based DNA profiling. However, a recent study employed deep sequencing to explore extremely rare mutations in the nuclear genome and reported that the mutation analysis could be applied to differentiate between MZ twins. Compared with the nuclear genome, the mitochondrial DNA (mtDNA) exhibits higher mutation rates due to fewer DNA repair mechanisms in the mitochondrial genome (mtGenome) and the lack of proofreading capability of the mtDNA polymerase. In a previous study, we used Illumina ultra-deep sequencing to describe point heteroplasmy (PHP) and nucleotide variant of the mtGenomes in venous blood samples of MZ twins. In the present study, we characterized minor differences of the mtGenomes in three tissue samples from seven sets of MZ twins using Ion Torrent semiconductor sequencing (Thermo Fisher Ion S5 XL system) and commercialized mtGenome sequencing kit (Precision ID mtDNA Whole Genome Panel). PHP was observed in blood samples from one set of MZ twins and in saliva samples from two sets of twins, but it presented in hair shaft samples from all seven sets of MZ twins. Overall, the coding region of the mtGenome exhibits more PHPs than the control region. The results of this study have further attested the competence of mtGenome sequencing in differentiating between MZ twins, and that among the three kinds of samples tested, hair shaft is more likely to accumulate minor differences in the mtGenomes of MZ twins.

Whole Mitochondrial Genome Detection and Analysis of Two- to Four-Generation Maternal Pedigrees Using a New Massively Parallel Sequencing Panel

Mitochondrial DNA (mtDNA) is an effective genetic marker in forensic practice, especially for aged bones and hair shafts. Detection of the whole mitochondrial genome (mtGenome) using traditional Sanger-type sequencing is laborious and time-consuming. Additionally, its ability to distinguish point heteroplasmy (PHP) and length heteroplasmy (LHP) is limited. The application of massively parallel sequencing in mtDNA detection helps researchers to study the mtGenome in-depth. The ForenSeq mtDNA Whole Genome Kit, which contains a total of 245 short amplicons, is one of the multiplex library preparation kits for the mtGenome. We used this system to detect the mtGenome in the blood samples and hair shafts of thirty-three individuals from eight two-generation pedigrees, one three-generation pedigree, and one four-generation pedigree. High-quality sequencing results were obtained. Ten unique mtGenome haplotypes were observed in the mothers from the ten pedigrees. A total of 26 PHPs were observed using the interpretation threshold of 6%. Eleven types of LHPs in six regions were evaluated in detail. When considering homoplasmic variants only, consistent mtGenome haplotypes were observed between the twice-sequenced libraries and between the blood and hair shafts from the same individual and among maternal relatives in the pedigrees. Four inherited PHPs were observed, and the remainder were de novo/disappearing PHPs in the pedigrees. Our results demonstrate the effective capability of the ForenSeq mtDNA Whole Genome Kit to generate the complete mtGenome in blood and hair shafts, as well as the complexity of mtDNA haplotype comparisons between different types of maternal relatives when heteroplasmy is considered.

DNA and protein analyses of hair in forensic genetics

Hair is one of the most common pieces of biological evidence found at a crime scene and plays an essential role in forensic investigation. Hairs, especially non-follicular hairs, are usually found at various crime scenes, either by natural shedding or by forcible shedding. However, the genetic material in hairs is usually highly degraded, which makes forensic analysis difficult. As a result, the value of hair has not been fully exploited in forensic investigations and trials. In recent years, with advances in molecular biology, forensic analysis of hair has achieved remarkable strides and provided crucial clues in numerous cases. This article reviews recent developments in DNA and protein analysis of hair and attempts to provide a comprehensive solution to improve forensic hair analysis.

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.

Evaluating the suitability of current mitochondrial DNA interpretation guidelines for multigenerational whole mitochondrial genome comparisons

Sanger sequencing of the mitochondrial DNA (mtDNA) control region was previously the only method available for forensic casework involving degraded samples from skeletal remains. The introduction of Next Generation Sequencing (NGS) has transformed genetic data generation and human identification using mtDNA. Whole mitochondrial genome (mtGenome) analysis is now being introduced into forensic laboratories around the world to analyze historical remains. Research into large pedigrees using the mtGenome is critical to evaluate currently available interpretation guidelines for mtDNA analysis, which were developed for comparisons using the control region. This study included mtGenomes from 225 individuals from the last four generations of the Norfolk Island (NI) genetic isolate pedigree consisting of 49 distinct maternal lineages. The data from these individuals were arranged into 2339 maternally related pairs separated by up to 18 meioses. Our results show that 97.3% of maternally related pairs were concordant at all nucleotide positions, resulting in the correct interpretation of “Cannot Exclude”; 2.7% of pairs produced an “Inconclusive” result, and there were no instances of false exclusion. While these results indicate that existing guidelines are suitable for multigenerational whole mtGenome analysis, we recommend caution be taken when classifying heteroplasmic changes as differences for human identification. Our data showed the classification of heteroplasmic changes as differences increases the prevalence of inconclusive identification by 6%, with false exclusions observed in 0.34% of pairs examined. Further studies of multigenerational pedigrees, however, are needed to validate mtGenome interpretation guidelines for historical case work to more fully utilize emerging advancements.

Whole mitochondrial genome analysis of the Daur ethnic minority from Hulunbuir in the Inner Mongolia Autonomous Region of China

BACKGROUND

Mitochondrial DNA (mtDNA) variations are often associated with bioenergetics, disease, and speciation and can be used to track the history of women. Although advances in massively parallel sequencing (MPS) technology have greatly promoted our understanding of the population's history (especially genome-wide data and whole Y chromosome sequencing), the whole mtDNA sequence of many important groups has not been fully studied. In this study, we employed whole mitogenomes of 209 healthy and unrelated individuals from the Daur group, a Mongolic-speaking representative population of the indigenous groups in the Heilongjiang River basin (also known as the Amur River basin).

RESULTS

The dataset presented 127 distinct mtDNA haplotypes, resulting in a haplotype diversity of 0.9933. Most of haplotypes were assigned to eastern Eurasian-specific lineages, such as D4 (19.62%), B4 (9.09%), D5 (7.66%) and M7 (4.78%). Population comparisons showed that the Daurians do have certain connections with the ancient populations in the Heilongjiang River basin but the matrilineal genetic composition of the Daur group was also greatly influenced by other non-Mongolic groups from neighboring areas.

CONCLUSIONS

Collectively, the whole mtDNA data generated in the present study will augment the existing mtDNA database. Our study provides genetic links between the Daur population and the aborigine peoples from Siberia and the Amur-Ussuri Region. But on the whole, compared with other Mongolic-speaking groups, the modern Daur population is closer to the East Asian ancestry group.

Benchmarking the Effectiveness and Accuracy of Multiple Mitochondrial DNA Variant Callers: Practical Implications for Clinical Application

Mitochondrial DNA (mtDNA) mutations contribute to human disease across a range of severity, from rare, highly penetrant mutations causal for monogenic disorders to mutations with milder contributions to phenotypes. mtDNA variation can exist in all copies of mtDNA or in a percentage of mtDNA copies and can be detected with levels as low as 1%. The large number of copies of mtDNA and the possibility of multiple alternative alleles at the same DNA nucleotide position make the task of identifying allelic variation in mtDNA very challenging. In recent years, specialized variant calling algorithms have been developed that are tailored to identify mtDNA variation from whole-genome sequencing (WGS) data. However, very few studies have systematically evaluated and compared these methods for the detection of both homoplasmy and heteroplasmy. A publicly available synthetic gold standard dataset was used to assess four mtDNA variant callers (Mutserve, mitoCaller, MitoSeek, and MToolBox), and the commonly used Genome Analysis Toolkit “best practices” pipeline, which is included in most current WGS pipelines. We also used WGS data from 126 trios and calculated the percentage of maternally inherited variants as a metric of calling accuracy, especially for homoplasmic variants. We additionally compared multiple pathogenicity prediction resources for mtDNA variants. Although the accuracy of homoplasmic variant detection was high for the majority of the callers with high concordance across callers, we found a very low concordance rate between mtDNA variant callers for heteroplasmic variants ranging from 2.8% to 3.6%, for heteroplasmy thresholds of 5% and 1%. Overall, Mutserve showed the best performance using the synthetic benchmark dataset. The analysis of mtDNA pathogenicity resources also showed low concordance in prediction results. We have shown that while homoplasmic variant calling is consistent between callers, there remains a significant discrepancy in heteroplasmic variant calling. We found that resources like population frequency databases and pathogenicity predictors are now available for variant annotation but still need refinement and improvement. With its peculiarities, the mitochondria require special considerations, and we advocate that caution needs to be taken when analyzing mtDNA data from WGS data.

Exploring statistical weight estimates for mitochondrial DNA matches involving heteroplasmy

Massively parallel sequencing (MPS) of mitochondrial (mt) DNA allows forensic laboratories to report heteroplasmy on a routine basis. Statistical approaches will be needed to determine the relative frequency of observing an mtDNA haplotype when including the presence of a heteroplasmic site. Here, we examined 1301 control region (CR) sequences, collected from individuals in four major population groups (European, African, Asian, and Latino), and covering 24 geographically distributed haplogroups, to assess the rates of point heteroplasmy (PHP) on an individual and nucleotide position (np) basis. With a minor allele frequency (MAF) threshold of 2%, the data was similar across population groups, with an overall PHP rate of 37.7%, and the majority of heteroplasmic individuals (77.3%) having only one site of heteroplasmy. The majority (75.2%) of identified PHPs had an MAF of 2-10%, and were observed at 12.6% of the nps across the CR. Both the broad and phylogenetic testing suggested that in many cases the low number of observations of heteroplasmy at any one np results in a lack of statistical association. The posterior frequency estimates, which skew conservative to a degree depending on the sample size in a given haplogroup, had a mean of 0.152 (SD 0.134) and ranged from 0.031 to 0.83. As expected, posterior frequency estimates decreased in accordance with 1/n as the sample size (n) increased. This provides a proposed conservative statistical framework for assessing haplotype/heteroplasmy matches when applying an MPS technique in forensic cases and will allow for continual refinement as more data is generated, both within the CR and across the mitochondrial genome.

The Value of Whole-Genome Sequencing for Mitochondrial DNA Population Studies: Strategies and Criteria for Extracting High-Quality Mitogenome Haplotypes

Whole-genome sequencing (WGS) data present a readily available resource for mitochondrial genome (mitogenome) haplotypes that can be utilized for genetics research including population studies. However, the reconstruction of the mitogenome is complicated by nuclear mitochondrial DNA (mtDNA) segments (NUMTs) that co-align with the mtDNA sequences and mimic authentic heteroplasmy. Two minimum variant detection thresholds, 5% and 10%, were assessed for the ability to produce authentic mitogenome haplotypes from a previously generated WGS dataset. Variants associated with NUMTs were detected in the mtDNA alignments for 91 of 917 (~8%) Swedish samples when the 5% frequency threshold was applied. The 413 observed NUMT variants were predominantly detected in two regions (nps 12,612–13,105 and 16,390–16,527), which were consistent with previously documented NUMTs. The number of NUMT variants was reduced by ~97% (400) using a 10% frequency threshold. Furthermore, the 5% frequency data were inconsistent with a platinum-quality mitogenome dataset with respect to observed heteroplasmy. These analyses illustrate that a 10% variant detection threshold may be necessary to ensure the generation of reliable mitogenome haplotypes from WGS data resources.

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.

Weaving Mitochondrial DNA and Y-Chromosome Variation in the Panamanian Genetic Canvas

The Isthmus of Panama was a crossroads between North and South America during the continent’s first peopling (and subsequent movements) also playing a pivotal role during European colonization and the African slave trade. Previous analyses of uniparental systems revealed significant sex biases in the genetic history of Panamanians, as testified by the high proportions of Indigenous and sub-Saharan mitochondrial DNAs (mtDNAs) and by the prevalence of Western European/northern African Y chromosomes. Those studies were conducted on the general population without considering any self-reported ethnic affiliations. Here, we compared the mtDNA and Y-chromosome lineages of a new sample collection from 431 individuals (301 males and 130 females) belonging to either the general population, mixed groups, or one of five Indigenous groups currently living in Panama. We found different proportions of paternal and maternal lineages in the Indigenous groups testifying to pre-contact demographic events and genetic inputs (some dated to Pleistocene times) that created genetic structure. Then, while the local mitochondrial gene pool was marginally involved in post-contact admixtures, the Indigenous Y chromosomes were differentially replaced, mostly by lineages of western Eurasian origin. Finally, our new estimates of the sub-Saharan contribution, on a more accurately defined general population, reduce an apparent divergence between genetic and historical data.

From Forensics to Clinical Research: Expanding the Variant Calling Pipeline for the Precision ID mtDNA Whole Genome Panel

Despite a multitude of methods for the sample preparation, sequencing, and data analysis of mitochondrial DNA (mtDNA), the demand for innovation remains, particularly in comparison with nuclear DNA (nDNA) research. The Applied Biosystems™ Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific, USA) is an innovative library preparation kit suitable for degraded samples and low DNA input. However, its bioinformatic processing occurs in the enterprise Ion Torrent Suite™ Software (TSS), yielding BAM files aligned to an unorthodox version of the revised Cambridge Reference Sequence (rCRS), with a heteroplasmy threshold level of 10%. Here, we present an alternative customizable pipeline, the PrecisionCallerPipeline (PCP), for processing samples with the correct rCRS output after Ion Torrent sequencing with the Precision ID library kit. Using 18 samples (3 original samples and 15 mixtures) derived from the 1000 Genomes Project, we achieved overall improved performance metrics in comparison with the proprietary TSS, with optimal performance at a 2.5% heteroplasmy threshold. We further validated our findings with 50 samples from an ongoing independent cohort of stroke patients, with PCP finding 98.31% of TSS's variants (TSS found 57.92% of PCP's variants), with a significant correlation between the variant levels of variants found with both pipelines.

Internal validation and improvement of mitochondrial genome sequencing using the Precision ID mtDNA Whole Genome Panel

With the recent advances in next-generation sequencing (NGS), mitochondrial whole-genome sequencing has begun to be applied to the field of the forensic biology as an alternative to the traditional Sanger-type sequencing (STS). However, experimental workflows, commercial solutions, and output data analysis must be strictly validated before being implemented into the forensic laboratory. In this study, we performed an internal validation for an NGS-based typing of the entire mitochondrial genome using the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific) on the Ion S5 sequencer (Thermo Fisher Scientific). Concordance, repeatability, reproducibility, sensitivity, and heteroplasmy detection analyses were assessed using the 2800 M and 9947A standard control DNA as well as typical casework specimens, and results were compared with conventional Sanger sequencing and another NGS sequencer in a different laboratory. We discuss the strengths and limitations of this approach, highlighting some issues regarding noise thresholds and heteroplasmy detection, and suggesting solutions to mitigate these effects and improve overall data interpretation. Results confirmed that the Precision ID Whole mtDNA Genome Panel is highly reproducible and sensitive, yielding useful full mitochondrial DNA sequences also from challenging DNA specimens, thus providing further support for its use in forensic practice.

Mitochondrial DNA in forensic use

Genetic analysis of mitochondrial DNA (mtDNA) has always been a useful tool for forensic geneticists, mainly because of its ubiquitous presence in biological material, even in the absence of nuclear DNA. Sequencing, however, is not a skill that is part of the routine forensic analysis because of the relative rarity of requests, and the need for retention of necessary skill sets and associated accreditation issues. While standard Sanger sequencing may be relatively simple, many requests are made in the face of compromised biological samples. Newer technologies, provided through massively parallel sequencing (MPS), will increase the opportunity for scientists to include this tool in their routine, particularly for missing person investigations. MPS has also enabled a different approach to sequencing that can increase sensitivity in a more targeted approach. In these circumstances it is likely that only a laboratory that specialises in undertaking forensic mtDNA analysis will be able to take these difficult cases forward, more so because reviews of the literature have revealed significantly high levels of typing errors in publications reporting mtDNA sequences. The forensic community has set out important guidelines, not only in the practical aspects of analysis, but also in the interpretation of that sequence to ensure that accurate comparisons can be made. Analysis of low-level, compromised and ancient DNA is not easy, however, as contamination is extremely difficult to eliminate and circumstances leading to sequencing errors are all too easily introduced. These problems, and solutions, are discussed in the article in relation to several historic cases.

Massively parallel sequencing of human skeletal remains in Vietnam using the precision ID mtDNA control region panel on the Ion S5™ system

Mitochondrial DNA (mtDNA) analysis using Sanger sequencing has been a routine practice for the identification of human skeletal remains. However, this process is usually challenging since DNA from the remains is highly degraded and at low concentration. Recently, the advent and implementation of massively parallel sequencing (MPS) have been offered the ability to improve mtDNA sequence data for forensic analysis. To assess the utility of the Ion S5™ system - an MPS platform for mtDNA analysis in challenging samples, we sequenced the mitochondrial control region of 52 age-old skeletal remains. Using the Precision ID mtDNA Control Region Panel, 50 full and two partial control region haplotypes at relatively high mean coverage of 2494 × were achieved for variant calling. Further variant analysis at 10% threshold for point heteroplasmy showed high degradation degree in terms of DNA damage in our bone samples. A higher point heteroplasmy threshold of 20% was required to diminish most of background noise caused by the damage. The results from this study indicated the potential application of the Ion S5™ system in sequencing degraded samples in Vietnam and provided valuable data sources for forensic analyses in the future.

Comparison of whole genome sequencing and targeted sequencing for mitochondrial DNA

Mitochondrial dysfunction has emerged to be associated with a broad spectrum of diseases, and there is an increasing demand for accurate detection of mitochondrial DNA (mtDNA) variants. Whole genome sequencing (WGS) has been the dominant sequencing approach to identify genetic variants in recent decades, but most studies focus on variants on the nuclear genome. Whole genome sequencing is also costly and time consuming. Sequencing specifically targeted for mtDNA is commonly used in the diagnostic settings and has lower costs. However, there is a lack of pairwise comparisons between these two sequencing approaches for calling mtDNA variants on a population basis. In this study, we compared WGS and mtDNA-targeted sequencing (targeted-seq) in analyzing mitochondrial DNA from 1499 participants recruited into the Severe Asthma Research Program (SARP). Our study reveals that targeted-sequencing and WGS have comparable capacity to determine genotypes and to call haplogroups and homoplasmies on mtDNA. However, there exists a large variability in calling heteroplasmies, especially for low-frequency heteroplasmies, which indicates that investigators should be cautious about heteroplasmies acquired from different sequencing methods. Further research is highly desired to improve variant detection methods for mitochondrial DNA.

Interpreting NUMTs in forensic genetics: Seeing the forest for the trees

Nuclear mitochondrial DNA (mtDNA) segments (NUMTs) were discovered shortly after sequencing the first human mitochondrial genome. They have earlier been considered to represent archaic elements of ancient insertion events, but modern sequencing technologies and growing databases of mtDNA and NUMT sequences confirm that they are abundant and some of them phylogenetically young. Here, we build upon mtDNA/NUMT review articles published in the mid 2010 s and focus on the distinction of NUMTs and other artefacts that can be observed in aligned sequence reads, such as mixtures (contamination), point heteroplasmy, sequencing error and cytosine deamination. We show practical examples of the effect of the mtDNA enrichment method on the representation of NUMTs in the mapped sequence data and discuss methods to bioinformatically filter NUMTs from mtDNA reads.

Contamination detection in sequencing studies using the mitochondrial phylogeny

Within-species contamination is a major issue in sequencing studies, especially for mitochondrial studies. Contamination can be detected by analyzing the nuclear genome or by inspecting polymorphic sites in the mitochondrial genome (mtDNA). Existing methods using the nuclear genome are computationally expensive, and no appropriate tool for detecting sample contamination in large-scale mtDNA data sets is available. Here we present haplocheck, a tool that requires only the mtDNA to detect contamination in both targeted mitochondrial and whole-genome sequencing studies. Our in silico simulations and amplicon mixture experiments indicate that haplocheck detects mtDNA contamination accurately and is independent of the phylogenetic distance within a sample mixture. By applying haplocheck to The 1000 Genomes Project Consortium data, we further evaluate the application of haplocheck as a fast proxy tool for nDNA-based contamination detection using the mtDNA and identify the mitochondrial copy number within a mixture as a critical component for the overall accuracy. The haplocheck tool is available both as a command-line tool and as a cloud web service producing interactive reports that facilitates the navigation through the phylogeny of contaminated samples.

Analyzing Low-Level mtDNA Heteroplasmy-Pitfalls and Challenges from Bench to Benchmarking

Massive parallel sequencing technologies are promising a highly sensitive detection of low-level mutations, especially in mitochondrial DNA (mtDNA) studies. However, processes from DNA extraction and library construction to bioinformatic analysis include several varying tasks. Further, there is no validated recommendation for the comprehensive procedure. In this study, we examined potential pitfalls on the sequencing results based on two-person mtDNA mixtures. Therefore, we compared three DNA polymerases, six different variant callers in five mixtures between 50% and 0.5% variant allele frequencies generated with two different amplification protocols. In total, 48 samples were sequenced on Illumina MiSeq. Low-level variant calling at the 1% variant level and below was performed by comparing trimming and PCR duplicate removal as well as six different variant callers. The results indicate that sensitivity, specificity, and precision highly depend on the investigated polymerase but also vary based on the analysis tools. Our data highlight the advantage of prior standardization and validation of the individual laboratory setup with a DNA mixture model. Finally, we provide an artificial heteroplasmy benchmark dataset that can help improve somatic variant callers or pipelines, which may be of great interest for research related to cancer and aging.

Reducing the Number of Mismatches between Hairs and Buccal References When Analysing mtDNA Heteroplasmic Variation by Massively Parallel Sequencing

In forensics, mitochondrial DNA (mtDNA) analysis is foremost applied to rootless hairs often lacking detectable nuclear DNA. Sanger sequencing is the routine mtDNA method in most forensic laboratories, even though interpretation of mixed samples and heteroplasmic sites can be challenging. Individuals may hold cells with low-level heteroplasmy variants below the detection threshold and other cells where this minor variant is the major one. This difference may be interpreted as a mismatch between reference and evidentiary trace samples, such as buccal specimens and rootless hairs. Such mismatches may be solved by Massively Parallel Sequencing (MPS), allowing more sensitive quantitative analysis for mixed positions than Sanger. The mtDNA control region was analysed in buccal reference samples from 26 individuals and 475 corresponding hairs by MPS and compared to Sanger sequencing data generated on the same samples. With MPS, mixed contributions down to 3% were regarded, leading to a substantial increase in the frequency of heteroplasmy. Our results demonstrate that previously reported mismatches between buccal reference and hair shaft samples by Sanger are detected as low-level heteroplasmy by MPS. A detailed overview of buccal and hair heteroplasmy is provided and implications for MPS-based mtDNA analysis in the context of forensic cases are discussed.

The Genetic Structure of Chinese Hui Ethnic Group Revealed by Complete Mitochondrial Genome Analyses Using Massively Parallel Sequencing

Mitochondrial DNA (mtDNA), coupled with maternal inheritance and relatively high mutation rates, provides a pivotal way for us to investigate the formation histories of populations. The Hui minority with Islamic faith is one of the most widely distributed ethnic groups in China. However, the exploration of Hui's genetic architecture from the complete mitochondrial genome perspective has not been detected yet. Therefore, in this study, we employed the complete mitochondrial genomes of 98 healthy and unrelated individuals from Northwest China, as well as 99 previously published populations containing 7274 individuals from all over the world as reference data, to comprehensively dissect the matrilineal landscape of Hui group. Our results demonstrated that Hui group exhibited closer genetic relationships with Chinese Han populations from different regions, which was largely attributable to the widespread of haplogroups D4, D5, M7, B4, and F1 in these populations. The demographic expansion of Hui group might occur during the Late Pleistocene. Finally, we also found that Hui group might have gene exchanges with Uygur, Tibetan, and Tajik groups in different degrees and retained minor genetic imprint of European-specific lineages, therefore, hinting the existence of multi-ethnic integration events in shaping the genetic landscape of Chinese Hui group.

Platinum-Quality Mitogenome Haplotypes from United States Populations

A total of 1327 platinum-quality mitochondrial DNA haplotypes from United States (U.S.) populations were generated using a robust, semi-automated next-generation sequencing (NGS) workflow with rigorous quality control (QC). The laboratory workflow involved long-range PCR to minimize the co-amplification of nuclear mitochondrial DNA segments (NUMTs), PCR-free library preparation to reduce amplification bias, and high-coverage Illumina MiSeq sequencing to produce an average per-sample read depth of 1000 × for low-frequency (5%) variant detection. Point heteroplasmies below 10% frequency were confirmed through replicate amplification, and length heteroplasmy was quantitatively assessed using a custom read count analysis tool. Data analysis involved a redundant, dual-analyst review to minimize errors in haplotype reporting with additional QC checks performed by EMPOP. Applying these methods, eight sample sets were processed from five U.S. metapopulations (African American, Caucasian, Hispanic, Asian American, and Native American) corresponding to self-reported identity at the time of sample collection. Population analyses (e.g., haplotype frequencies, random match probabilities, and genetic distance estimates) were performed to evaluate the eight datasets, with over 95% of haplotypes unique per dataset. The platinum-quality mitogenome haplotypes presented in this study will enable forensic statistical calculations and thereby support the usage of mitogenome sequencing in forensic laboratories.

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.

Damage patterns observed in mtDNA control region MPS data for a range of template concentrations and when using different amplification approaches

Massively parallel sequencing (MPS) of mitochondrial (mt) DNA allows practitioners the ability to fully resolve heteroplasmic sites. In forensic DNA analysis, identifying heteroplasmy (a naturally occurring mixture of two mtDNA profiles) can provide additional mtDNA profile information which can lead to an increase in the discrimination potential of an mtDNA match between an evidentiary sample and reference source. Forensic samples such as hair and skeletal remains, especially older, more compromised samples, can often exhibit DNA damage. Because both damage and heteroplasmy can manifest as a mixture of two nucleotides, it is important to differentiate between the two conditions when interpreting mtDNA MPS data. In this study, DNA damage was applied under controlled conditions to samples containing a range of template concentrations, including some with identified heteroplasmy. Damage was applied via storage in water at room temperature on samples diluted before or after storage to mimic low template scenarios. Damage was assessed with respect to the following areas: mtDNA quantification and degradation ratios, MPS read depth, MPS profile results, overall damage rates, and the interpretation of heteroplasmy. Datasets were generated to assess and compare two different amplification and library preparation strategies: the Promega PowerSeq™ CRM Nested System kit and a 1.16 kb target amplicon of the entire mtDNA control region followed by a Nextera® XT library preparation. The results of this study provide an evaluation of the Promega 10-plex MPS procedure as an improved process to mitigate the impact of mtDNA damage on low template samples. Some of the negative effects of damage observed in this study were a decrease in mtDNA yield by 20-30% and lower quality MPS sequencing results. These effects were observed more frequently when samples were diluted prior to inducing damage, illustrating that low template samples are more susceptible to damage. The findings of this study will assist forensic laboratories in differentiating between damage and heteroplasmy, which is essential when developing robust mtDNA MPS interpretation guidelines such as setting appropriate reporting thresholds.

Marked mitochondrial genetic variation in individuals and populations of the carcinogenic liver fluke Clonorchis sinensis

Clonorchiasis is a neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, and is often associated with a malignant form of bile duct cancer (cholangiocarcinoma). Although some aspects of the epidemiology of clonorchiasis are understood, little is known about the genetics of C. sinensis populations. Here, we conducted a comprehensive genetic exploration of C. sinensis from endemic geographic regions using complete mitochondrial protein gene sets. Genomic DNA samples from C. sinensis individuals (n = 183) collected from cats and dogs in China (provinces of Guangdong, Guangxi, Hunan, Heilongjiang and Jilin) as well as from rats infected with metacercariae from cyprinid fish from the Russian Far East (Primorsky Krai region) were deep sequenced using the BGISEQ-500 platform. Informatic analyses of mitochondrial protein gene data sets revealed marked genetic variation within C. sinensis; significant variation was identified within and among individual worms from distinct geographical locations. No clear affiliation with a particular location or host species was evident, suggesting a high rate of dispersal of the parasite across endemic regions. The present work provides a foundation for future biological, epidemiological and ecological studies using mitochondrial protein gene data sets, which could aid in elucidating associations between particular C. sinensis genotypes/haplotypes and the pathogenesis or severity of clonorchiasis and its complications (including cholangiocarcinoma) in humans.

Clonorchiasis is an important neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, which can induce malignant bile duct cancer (cholangiocarcinoma). Little precise information is available on the biology, epidemiology and population genetics of C. sinensis. For this reason, we explored here the genetic composition of C. sinensis populations in distinct endemic areas in China and Russia. Using a deep sequencing-informatic approach, we revealed marked mitochondrial genetic variation within and between individuals and populations of C. sinensis, with no particular affiliation with geographic or host origin. These molecular findings and the methodology established should underpin future genetic studies of C. sinensis causing human clonorchiasis and associated complications (cancer) as well as transmission patterns in endemic regions.

Ancestry and phenotype predictions from touch DNA using massively parallel sequencing

Direct PCR can be used to successfully generate full STR profiles from DNA present on the surface of objects. STR profiles are only of use in cases where a potential donor profile is available for comparison, and DNA is of sufficient DNA quality and quantity to generate a reliable profile. Often, no donor information is available and only trace DNA is present on items. As a result, alternative techniques are required to generate genetic data that can provide investigative leads. Massively parallel sequencing (MPS) offers the ability to detect trace levels of DNA and improve DNA analysis success from touched items. Here, we present the first application of direct PCR coupled with MPS to generate forensic intelligence SNP data from latent DNA. The panels assessed are (1) the HIrisplex System that targets 24 SNPs to simultaneously predict hair and eye, and (2) the Precision ID Ancestry Panel that targets 165 autosomal SNPs indicative of biogeographic ancestry. For each panel, we analysed 60 touched samples across five individuals and four substrates (glass slide, fuse, zip-lock bag and wire) using Ion AmpliSeq Library Preparation Kit on the automated Ion Chef System and Ion Torrent PGM. We examine the SNP recovery, concordance with reference samples and the genotype reproducibility from different substrates and donors. The results demonstrate the application of this approach for obtaining informative genetic from trace amounts of DNA.

Genomic Landscape of the Mitochondrial Genome in the United Arab Emirates Native Population

In order to assess the genomic landscape of the United Arab Emirates (UAE) mitogenome, we sequenced and analyzed the complete genomes of 232 Emirate females mitochondrial DNA (mtDNA) within and compared those to Africa. We investigated the prevalence of haplogroups, genetic variation, heteroplasmy, and demography among the UAE native population with diverse ethnicity and relatively high degree of consanguinity. We identified 968 mtDNA variants and high-resolution 15 haplogroups. Our results show that the UAE population received enough gene flow from Africa represented by the haplogroups L, U6, and M1, and that 16.8% of the population has an eastern provenance, depicted by the U haplogroup and the M Indian haplogroup (12%), whereas western Eurasian and Asian haplogroups (R, J, and K) represent 11 to 15%. Interestingly, we found an ancient migration present through the descendant of L (N1 and X) and other sub-haplogroups (L2a1d and L4) and (L3x1b), which is one of the oldest evolutionary histories outside of Africa. Our demographic analysis shows no population structure among populations, with low diversity and no population differentiation. In addition, we show that the transmission of mtDNA in the UAE population is under purifying selection with hints of diversifying selection on ATP8 gene. Last, our results show a population bottleneck, which coincides with the Western European contact (1400 ybp). Our study of the UAE mitogenomes suggest that several maternal lineage migratory episodes liking African–Asian corridors occurred since the first modern human emerges out of Africa.

Massively parallel sequencing of mitogenome sequences reveals the forensic features and maternal diversity of tai-kadai-speaking hlai islanders

As a single maternally inherited locus, human mitochondrial DNA (mtDNA) is geographically arranged and plays a key role in forensic applications. Hlai population has been evidenced as the most typical and unmixed representative of the Tai-Kadai-speaking populations via genome-wide analyses. However, forensic features and maternal diversity of the complete mitogenomes in this Tai-Kadai ancestrally related population are scarce. Thus, we sequenced the complete mitogenomes in 127 Hainan Hlais and found 109 distinct haplotypes belonging to 43 terminal haplogroups resulting in the haplotype diversity of 0.9970. Our results of comprehensive population comparisons showed that Hlai islanders had a close genetic affinity with Tai-Kadai-speaking populations from Southeast Asia, which is consistent with the back-migration of Chinese Neolithic farmers into this region via the inland route. Besides, maternally genetic evidence further revealed a close genetic relationship between Tai-Kadai-speaking and Austronesian-speaking populations when only East Asian dataset was considered, which is consistent with the common origin from Yangtze rice farmers and then spread southward along the inland and coastal routes, respectively. In the reconstructed phylogenetic tree and median-joining networks, the vast majority of Hlais were clustered in exclusive clades, which demonstrated that Hlai people probably had undergone founder effect or genetic bottleneck in their history, and remained genetically isolated for a long time. Collectively, Hainan Hlai did not exhibit detectable maternal gene flow from surrounding or incoming populations. Mitogenome information generated in this study is a contribution in mitigating the underrepresentation of Chinese data in forensic mitogenetics and will assist geography-, metapopulation-, as well as phylogeny-based queries.

Characterization of background noise in MiSeq MPS data when sequencing human mitochondrial DNA from various sample sources and library preparation methods

Improved resolution of massively parallel sequencing (MPS) allows for the characterization of mitochondrial (mt) DNA heteroplasmy to levels previously unattainable with traditional sequencing approaches. An essential criterion for the reporting of heteroplasmy is the ability of the MPS method to distinguish minor sequence variants (MSVs) from system noise, or error. Therefore, an assessment of the background noise in the MPS method is desirable to identify the point at which reliable data can be reported. Substitution and sequence specific error (SSE) was evaluated for a variety of sample types and two library preparations. Substitution error rates ranged from 0.18 to 0.49 per 100 nucleotides with C positions generally having the highest rate of misincorporation. Comparison of error rates across sample types indicated a significant increase for samples with damaged DNA. The positions of error were varied across datasets (pairwise concordance 0-68%), but had greater consistency within the damaged samples (80-96%). The most commonly observed motif preceding error in forward reads was CCG, while GGT was most common in reverse reads, both consistent with previous findings. The findings illustrate that for datasets containing samples with damaged DNA, reporting thresholds for heteroplasmy may have to be modified and individual sites with error levels exceeding thresholds should be scrutinized. Collectively, the shifting error profiles observed across the various sample types and library preparation methods demonstrates the need for an assessment of error under these varying circumstances. Characterization of the applicable background noise will help to ensure that thresholds are reliably set for detection of true MSVs.

Forensically relevant phylogeographic evaluation of mitogenome variation in the Basque Country

The Basque Country has been the focus of population (genetic) and evolutionary studies for decades, as it represents an interesting evolutionary feature: it is the only European country where a non-Indo-European language is still spoken today and, for which there are no known living or extinct relatives. Early studies that were based on anatomical and serological methods, along with subsequent molecular genetic investigations, contain controversial interpretations of their data. Additionally, the analysis of mitochondrial DNA, which is maternally inherited and thus suitable for the examination of the maternal phylogeny of the population, was the focus of some studies. Early mtDNA studies were however restricted to the information provided by the control region or its hypervariable segments only. These are known to harbour little phylogenetic information, particularly for haplogroup H that is dominant in Westeurasian populations including the Basques. Later studies analysed complete mitogenome sequences. Their information content is however limited, either because the number of samples was low, or because these studies only considered particular haplogroups. In this study we present the full mitogenome sequences of 178 autochthonous Basque individuals that were carefully selected based on their familial descent and discuss the observed phylogenetic signals in the light of earlier published findings. We confirm the presence of Basque-specific mtDNA lineages and extend the knowledge of these lineages by providing data on their distribution in comparison to other Basque and non-Basque populations. This dataset improves our understanding of the Basque mtDNA phylogeny and serves as a high-quality dataset that is provided via EMPOP for forensic genetic purposes.

Pronounced somatic bottleneck in mitochondrial DNA of human hair

Heteroplasmy is the presence of variable mitochondrial DNA (mtDNA) within the same individual. The dynamics of heteroplasmy allele frequency among tissues of the human body is not well understood. Here, we measured allele frequency at heteroplasmic sites in two to eight hairs from each of 11 humans using next-generation sequencing. We observed a high variance in heteroplasmic allele frequency among separate hairs from the same individual—much higher than that for blood and cheek tissues. Our population genetic modelling estimated the somatic bottleneck during embryonic follicle development of separate hairs to be only 11.06 (95% confidence interval 0.6–34.0) mtDNA segregating units. This bottleneck is much more drastic than somatic bottlenecks for blood and cheek tissues (136 and 458 units, respectively), as well as more drastic than, or comparable to, the germline bottleneck (equal to 25–32 or 7–10 units, depending on the study). We demonstrated that hair undergoes additional genetic drift before and after the divergence of mtDNA lineages of individual hair follicles. Additionally, we showed a positive correlation between donor's age and variance in heteroplasmy allele frequency in hair. These findings have important implications for forensics and for our understanding of mtDNA dynamics in the human body.

This article is part of the theme issue ‘Linking the mitochondrial genotype to phenotype: a complex endeavour’.

A simple method for sequencing the whole human mitochondrial genome directly from samples and its application to genetic testing

Next-generation sequencing (NGS) is a revolutionary sequencing technology for analyzing genomes. However, preprocessing methods for mitochondrial DNA (mtDNA) sequencing remain complex, and it is required to develop an authenticated preprocessing method. Here, we developed a simple and easy preprocessing method based on isothermal rolling circle mtDNA amplification using commercially available reagents. Isothermal amplification of mtDNA was successfully performed using both nanoliter quantities of plasma directly and 25 ng of total DNA extracted from blood or tissue samples. Prior to mtDNA amplification, it was necessary to treat the extracted total DNA with Exonuclease V, but it was not required to treat plasma. The NGS libraries generated from the amplified mtDNA provided sequencing coverage of the entire human mitochondrial genome. Furthermore, the sequencing results successfully detected heteroplasmy in patient samples, with called mutations and variants matching those from previous, independent, Sanger sequencing analysis. Additionally, a novel single nucleotide variant was detected in a healthy volunteer. The successful analysis of mtDNA using very small samples from patients is likely to be valuable in clinical medicine, as it could reduce patient discomfort by reducing sampling-associated damage to tissues. Overall, the simple and convenient preprocessing method described herein may facilitate the future development of NGS-based clinical and forensic mtDNA tests.

Comparison of whole mitochondrial genome variants between hair shafts and reference samples using massively parallel sequencing

Hair shafts are one of the most common types of evidence at crime scenes, and mitochondrial DNA (mtDNA) has been analyzed as a valuable genetic marker for hair shafts in forensic casework. However, the mtDNA analysis strategy may vary according to the quantity and quality of DNA extracted from a forensic sample and the available massively parallel sequencing (MPS) platform in laboratories. Forensic practitioners often have to interpret mtDNA sequences exhibiting point heteroplasmy (PHP) that are analyzed using different analytical methods. In the present study, the whole mitochondrial genome (mtGenome) variants of hair shaft samples obtained from 20 donors, which were sampled in duplicate and stored at room temperature for > 1 year, were analyzed using the Precision ID mtDNA Whole Genome Panel and Ion S5 system. The whole mtGenome variants of 20 blood and 20 buccal swab samples (reference samples) from the hair shaft donors were analyzed using the Nextera XT DNA Library Prep Kit and MiSeq System. A total of 20 unique mtGenome haplotypes were observed, and 56 PHP variants were identified across the 4 sets of tissue. When the major nucleotide of PHP was considered, 16 of 20 haplotypes of the hair shaft samples matched those of the corresponding blood and buccal swab samples. In four donors, the major nucleotide of PHP was inverted at one nucleotide position between the hair shaft and reference samples. However, the data obtained on MPS, showing high PHP resolution, provided substantial information to avoid false exclusion when comparing two haplotypes containing PHP with inverted major nucleotides. In conclusion, the present study demonstrates the utility of MPS in forensic casework in the comparative analysis of mtGenome variants containing PHP.

Whole mitochondrial genome analysis of highland Tibetan ethnicity using massively parallel sequencing

Mitochondrial DNA (mtDNA) is a key player in numerous multifaceted and intricate biological processes and plays a pivotal role in dissecting the peopling of different populations, due to its maternally inherited property and comparatively high mutation rate. In this study, 119 Tibetan individuals from the Muli Tibetan Autonomous County of China (average altitude above 3,000 m) were employed in mitochondrial genome (mitogenome) sequencing by massively parallel sequencing (MPS) techniques using the Precision ID mtDNA Whole Genome Panel on an Ion S5XL system. The dataset presented 88 distinct haplotypes, resulting in the haplotype diversity of 0.9909. The majority of haplotypes were assigned to East Asian lineages and the distribution of haplogroups of Muli Tibetan significantly differed from reference Tibetan populations. The maximum parsimony phylogeny reconstructed by 119 newly generated mitogenomes revealed 12 major Muli Tibetan lineages. Intriguingly, a Sherpa-specific sub-haplogroup A15c1 with the lack of mutations at 4216 and 15,924 was discerned in our dataset, which suggested that the maternal gene pool of Sherpas may derive from Tibetan populations. The shared haplogroups between Muli Tibetan and lowland Han Chinese hinted that these lineages may derive from non-Tibetans and have already differentiated before their arrival on the Tibetan Plateau. Furthermore, extensive pairwise population comparisons displayed that Muli Tibetan had a closer genetic relationship with ethnically or linguistically close Nyingtri Tibetan, Nyingtri Lhoba and Chamdo Tibetan populations. Genetic affinity was also observed between the Muli Tibetan and North Han Chinese. Collectively, the results generated in this study enriched the existing forensic mtDNA database and raised additional interest in the application of whole mitogenome sequencing in forensic investigations.