Strategies for complete mitochondrial genome sequencing on Ion Torrent PGM™ platform in forensic sciences

Mitochondrial genome sequencing with ForenSeq™ mtDNA Whole Genome Kit

Mitochondrial DNA (mtDNA) possesses unique genetic characteristics and plays a crucial role in forensic DNA analysis. Based on the massively parallel sequencing (MPS) technology alongside the short overlapping amplicon method, the ForenSeq™ mtDNA Whole Genome Kit is specifically designed for mtDNA analysis. In this study, we employ the ForenSeq™ mtDNA Whole Genome Kit on the MiSeq FGx® Sequencing System for mitochondrial genome (mtGenome) sequencing across nine consecutive runs and assess its MPS performance, such as read depth (RD), forward/reverse strand bias (SB), and mtGenome coverage. Furthermore, we conduct internal validations to evaluate its routine application in forensic sciences, including sensitivity, repeatability, concordance, degraded samples, inhibitor samples, case-type samples, and contamination. As a result, the Real-Time Analysis (RTA) and Universal Analysis Software (UAS) demonstrate proficient run metrics and MPS performance when 12-14 libraries are sequenced within a standard flow cell, achieving > 80 % of reads passing filter, > 80 % bases with ≥Q30, > 5000 × of the average RD, ∼1.0 of the average SB, > 70 % of the inter-amplicon balance, and > 99 % of the mtGenome coverage. The five most vulnerable amplicons, exhibiting low RD and high SB, are identified as nucleotide positions (nps) 1094-1177, 5858-5975, 6109-6149, 6718-6810, and 7021-7090. For tertiary data analysis, the substitutions are accurately reported by UAS, while insertions and deletions (indels), point heteroplasmies (PHPs), and/or length heteroplasmies (LHPs) still necessitate manual inspection. On average, 40 variants were found in 60 samples, ranging from 27 to 54. A total of 2426 variants are observed at 491 nps. Moreover, the workflow can yield repeatable and reproducible results, generate complete mtGenome profiles from ≥ 2 pg input gDNA for high quality samples/control DNA or ≥ 0.5 cm hair shafts, and recover more/complete mtGenome information from severely degraded samples (degradation index >10) and various types of case samples. If two rounds of purification are conducted, it can more effectively remove additional reaction components and enhance data recovery from the mtGenome, especially for low-input samples. The negative controls in three runs cover some reads, but these contaminations cannot compromise the mitochondrial analyses. In conclusion, the ForenSeq™ mtDNA Whole Genome Kit, including 234 short overlapping amplicons with an average size of 131 bp, can meet forensic needs on the whole mtGenome sequencing in real scenarios. In addition, the ten insights gained from this study may serve as a valuable reference for forensic scientists who are utilizing this kit.

Quantitation of human mitochondrial DNA and whole mtGenomes sequencing of fingernail/hair shaft samples

The analysis of mitochondrial DNA (mtDNA) is suitable for fingernail/hair shaft samples in forensic science. A successful mtDNA forensic analysis depends on the quantity and quality of the mtDNA. The application of massively parallel sequence techniques to the analysis of mtDNA has the potential to improve the recovery of genetic information from difficult forensic specimens and to increase the discrimination potential of mtDNA by capturing and comparing full mtGenomes. In this study, we constructed mtDNA-specific standard curves using real-time PCR. The 105-base pair target sequence facilitates the amplification of degraded DNA and is minimally homologous to non-human mtDNA. The results indicated that the assay enables the absolute quantification of down to 10 copies of mtDNA and provides a dynamic range of eight orders of magnitude. Whole mtGenome sequencing experiments demonstrated that as few as 2 000 mtDNA copies resulted in a successful full region amplification and sequencing. Moreover, the frequency of point heteroplasmy from one donor showed that hairs from the same donor have been found to differ within and among themselves and from other tissues, which could impact the interpretation of the results obtained in a forensic investigation.

Concurrent genotyping of mitochondrial DNA and nuclear DNA in rootless hair shafts and blood samples for enhanced analysis

Hair is an important type of biological evidence at crime scenes. However, the highly degraded nature of DNA fragments in hair shafts poses challenges for the detection of nuclear DNA (nuDNA) through capillary electrophoresis-based short tandem repeat (STR) genotyping. In this study, an all-in-one multiplex system named MGIEasy Signature Identification Library Prep Kit (MGI Tech, China) was employed to the simultaneous genotyping of both mitochondrial DNA (mtDNA) and nuDNA in hair shafts. This system is based on massively parallel sequencing (MPS) technology and encompasses Amelogenin, STRs, single nucleotide polymorphisms (SNPs) and mtDNA hypervariable regions (HVRs) in a single reaction. A total of 370 hair shafts, together with 180 blood samples as the references, were examined. The mtDNA analysis of 110 unrelated blood samples unveiled a total of 150 homoplasmic variants and 105 distinct haplotypes, revealing population polymorphisms in the Guangdong Han Chinese. The study also delved into the detection of mtDNA heteroplasmy, revealing 8.18 % and 16.36 % of individuals with point heteroplasmies (PHPs) in blood and hair shaft samples, respectively. Additionally, hair shafts with DNA extracted using the Investigator method yielded higher average depth of coverage (DoC), lower drop-out rate for SNP genotyping, higher nuDNA genotyped rates and success rates, than those using the MinElute method. In the longitudinal research, a gradual decrease in the total DoC of mtDNA fragments was observed along the length of the hair shaft, from the proximal root to the distal end. In contrast, the DoC of nuDNA exhibited a relatively stable pattern along the length of the hair shafts. The study contributes valuable insights into the simultaneous detection of nuDNA and mtDNA in hair shafts, emphasizing the need for optimized DNA extraction and detection methods for these highly degraded samples.

Whole mitochondrial genome analyses of Han population from Shandong of China using massively parallel sequencing

Introduction

Mitochondrial DNA (mtDNA) has been extensively utilized in archeology, human evolutionary genetics, and forensic genetic for over three decades, primarily due to its maternal inheritance and relatively high mutation rate. The Chinese Han, the largest and most widely distributed ethnic group in China, have been the focus of numerous genetic studies. However, the forensic parameters and genetic structure of the Shandong Han, specifically in relation to the whole mitochondrial genome, remain undocumented.

Methods

We performed whole mitochondrial genome sequencing on 141 unrelated Han individuals from Shandong province using massively parallel sequencing.

Results

A total of 135 unique mtDNA haplotypes were identified, classified into 105 haplogroups, resulting in a haplotype diversity value of 0.9993. The discriminatory capacity of whole mitochondrial genome was calculated at 0.9574, compared to 0.8936 when only the control region was analyzed. The majority of the haplogroups observed were specific to East Asian lineages, including D4, D5 and F1. Population comparisons revealed that the modern Shandong Han share genetic connections with ancient populations from the Yellow River and West Liao River basins. Additionally, the Shandong Han may have integrated a significant number of maternal lineages from other regions during their development. The demographic expansion of the Shandong Han is estimated to have occurred approximately 9,000 years ago, corresponding to the Neolithic period, a time of significant cultural and technological development.

Discussion

The dataset generated in this study is available in the EMPOP database under the accession number EMP00886 and will serve as an important mtDNA reference for forensic casework in China. The study of whole mitochondrial genome based on the analysis of matrilineal genetic structure of the Shandong Han population can help to enrich the forensic mtDNA reference database in East Asia and provide reference for future archeology and forensic genetics research.

Forensic efficiency evaluation of a mtDNA whole genome sequencing system constructed with long fragment amplification strategy on DNA nanoball sequencing platform

Mitochondrial DNA (mtDNA) is an important genetic marker for degraded biological sample identification, maternal pedigree tracing, and population genetic structure study owing to its characteristics of high copy number, anti-degradable ring structure, and maternal inheritance. Whole mtDNA genome sequencing is an optimal method for the analysis of mtDNA polymorphism and heterogeneity because it allows for the comprehensive use of maternal genetic information. However, because of lacking quantitative evaluations for sequencing data, the scientific interpretation standards for mtDNA sequencing results of the previously used sequencing systems are often different, and false positive or false negative results are prone to occur when faced with the interference of nuclear genomic DNA, or the heterogeneities of mtDNA sequence and structure. In this study, we evaluated a novel mtDNA whole genome sequencing system using long fragment amplification strategy on the DNA nanoball (DNB) sequencing platform. This system demonstrated high sequencing quality and specific mtDNA sequencing efficiencies on positive control DNA and FTA bloodstain samples, as the average Q20 and Q30 values of the corresponding samples were 97.17 % and 91.93 %; 97.37 % and 92.48 %, respectively. The mean mapping percentages for the reference sequences of whole genome DNA (wgDNA), mtDNA, and nuclear genomic DNA (ngDNA) in the corresponding samples were 99.98 %, 99.97 %, 0.03 %, and 99.91 %, 99.40 %, 0.60 %; respectively. The average error calling rates for the bases A, C, G, and T of the whole mtDNA genome were 0.2519 %, 0.2550 %, 0.2906 %; and 0.2392 %, respectively. The efficacy of heteroplasmy identification was assessed using a set of theoretical sites with predetermined rates. These sites were created by combining the samples with known mtDNA haplotypes in certain proportions. The absolute errors between observed and theoretical heteroplasmy values were 89.59 %, 74.68 %, 50.20 %, 12.65 %, 8.31 %, and 4.85 %, while the theoretical heteroplasmy values were 5 %, 10 %, 20 %, 80 %, 90 %, and 95 %, respectively. The absolute error exhibited relative stability when the mtDNA sequencing depth exceeded 500×. Furthermore, the system sequencing efficiency was also confirmed among different kinds of samples, and these samples included natural samples (e.g., peripheral blood samples preserved on FTA cards for 2 and 11 years, and on filter paper for 6 and 9 years), degraded samples, sensitivity samples, samples derived from various bodily fluids, and maternal pedigree samples. In summary, the whole mtDNA genome sequencing system used for forensic identification demonstrated high performance in analyzing mtDNA sequence information, and showed significant prospects for forensic application and maternal genetic research.

Forensic analysis and sequence variation of 133 STRs in the Hakka population

Introduction: Short Tandem Repeats (STRs) are highly valuable genetic markers in forensic science. However, the conventional PCR-CE technique has limitations, and the emergence of massively parallel sequencing (MPS) technology presents new opportunities for STR analysis. Yet, there is limited research on Chinese population diversity using MPS. Methods: In this study, we obtained genotype data for 52 A-STRs and 81 Y-STRs from the Hakka population in Meizhou, Guangdong, China, using the Forensic Analysis System Multiplecues SetB Kit on the MGISEQ-2000 platform. Results: Our findings demonstrate that these 133 STRs are highly efficient for forensic applications within the Meizhou Hakka population. Statistical analysis revealed Hobs values ranging from 0.61306 to 0.91083 and Hexp values ranging from 0.59156 to 0.91497 for A-STRs based on length polymorphism. For sequence polymorphism, Hobs values ranged from 0.61306 to 0.94586, and Hexp values fluctuated between 0.59156 and 0.94487. The CPE values were 1-5.0779620E-21 and 1-3.257436E-24 for length and sequence polymorphism, respectively, while the CPD values were 1-1.727007E-59 and 1-5.517015E-66, respectively. Among the 80 Y-STR loci, the HD values for length and sequence polymorphism were 0.99764282 and 0.99894195, respectively. The HMP values stood at 0.00418102 and 0.00288427, respectively, and the DC values were 0.75502742 and 0.83363803, respectively. For the 52 A-STR loci, we identified 554 and 989 distinct alleles based on length and sequence polymorphisms, respectively. For the 81 Y-STR loci, 464 and 652 unique alleles were detected at the length and sequence level, respectively. Population genetic analysis revealed that the Meizhou Hakka population has a close kinship relationship with the Asian populations THI and KOR based on length polymorphism data of A-STRs. Conversely, based on length polymorphism data of Y-STRs, the Meizhou Hakka population has the closest kinship relationship with the Henan Han population. Discussion: Overall, the variation information of repeat region sequences significantly enhances the forensic identification efficacy of STR genetic markers, providing an essential database for forensic individual and paternity testing in this region. Additionally, the data generated by our study will serve as a vital resource for research into the genetic structure and historical origins of the Meizhou Hakka population.

Mitochondrial DNA in Human Diversity and Health: From the Golden Age to the Omics Era

Mitochondrial DNA (mtDNA) is a small fraction of our hereditary material. However, this molecule has had an overwhelming presence in scientific research for decades until the arrival of high-throughput studies. Several appealing properties justify the application of mtDNA to understand how human populations are-from a genetic perspective-and how individuals exhibit phenotypes of biomedical importance. Here, I review the basics of mitochondrial studies with a focus on the dawn of the field, analysis methods and the connection between two sides of mitochondrial genetics: anthropological and biomedical. The particularities of mtDNA, with respect to inheritance pattern, evolutionary rate and dependence on the nuclear genome, explain the challenges of associating mtDNA composition and diseases. Finally, I consider the relevance of this single locus in the context of omics research. The present work may serve as a tribute to a tool that has provided important insights into the past and present of humankind.

High-resolution genotyping of 58 STRs in 635 Northern Han Chinese with MiSeq FGx ® Forensic Genomics System

Sequence polymorphisms were characterized at 27 autosomal STRs (A-STRs), 7 X chromosomal STRs (X-STRs), and 24 Y chromosomal STRs (Y-STRs) in 635 Northern Han Chinese with the ForenSeq DNA Signature Prep Kit on the MiSeq FGx Forensic Genomics System. Since repeat region (RR) and flanking region (FR) variation can be detected by massively parallel sequencing (MPS), the increase in the number of unique alleles and the average of gene diversity was 78.18% and 3.51% between sequence and length, respectively. A total of 74 novel RR variants were identified at 33 STRs compared with STRSeq and previous studies, and 13 FR variants (rs1770275883, rs2053373277, rs2082557941, rs1925525766, rs1926380862, rs1569322793, rs2051848492, rs2051848696, rs2016239814, rs2053269960, rs2044518192, rs2044536444, and rs2089968964) were first submitted to dbSNP. Also, 99.94% of alleles were concordant between the ForenSeq DNA Signature Prep Kit and commercial CE kits. Discordance resulted from the low performance at D22S1045 and occasionally at DYS392, flanking region deletions at D7S820 and DXS10074, and the strict alignment algorithm at DXS7132. Null alleles at DYS505 and DYS448 and multialleles at DYS387S1a/b, DYS385a/b, DYS448, DYS505, DXS7132, and HPRTB were validated with other MPS and CE kits. Thus, a high-resolution sequence-based (SB) and length-based (LB) allele frequencies dataset from Northern Han Chinese has been established already. As expected, forensic parameters increased significantly on combined power of discrimination (PD) and combined power of exclusion (PE) at A-STRs, mildly on combined PD and combined mean exclusion chance (MEC) at X-STRs, and barely on discrimination capacity (DC) at Y-STRs. Additionally, MiSeq FGx quality metrics and MPS performance were evaluated in this study, which presented the high-quality of the dataset at 20 consecutive runs, such as ≥ 60% bases with a quality score of 20 or higher (%≥ Q20), > 60% of effective reads, > 2000 × of depth of coverage (DoC), ≥ 60% of allele coverage ratio (ACR) or heterozygote balance, ≥ 70% of inter-locus balance, and ≤ 0.4 of the absolute value of observed minus expected heterozygosity (|H_exp - H_obs|). In conclusion, MiSeq FGx can help us generate a high-resolution and high-quality dataset for human identification and population genetic studies.

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.

A global snapshot of current opinions of next-generation sequencing technologies usage in forensics

The future of forensic DNA testing is being shaped by the research and usage of next-generation systems, which have increased the multiplexing capabilities of the field and the type and amount of genetic data that can be utilized for investigations. The NGS adoption for casework has been slow, albeit the plethora of data that has been published. This study evaluated the current opinions on sequencing in forensics. A 20-question online-survey focusing on NGS knowledge, training, and usage was distributed to 6001 forensic DNA researchers and practitioners worldwide. A total of 367 responses were obtained from all continents (North/South America (69.8%), Europe (21.2%), Asia (5.5%), Oceania (2.5%), and Africa (1%)). The respondents consisted of 50% practitioners, 31% researchers, and 19% both. Of these, 38% already own a next-gen sequencing instrument, and 13% are planning to purchase one. Overall, there exists an extensive knowledge on next-gen sequencing within the forensic community, including among laboratories that have not yet implemented this high-throughput technology in their workflows. Current usage focuses primarily on SNP analysis for investigative leads and mitochondrial DNA analysis while future applications included both STR and SNP testing applied to general casework. The major overall concerns respondents have for implementing a sequencing instrument include limited funding, staffing, lack of time, and the cost-effectiveness of providing this service. Specific technical concerns that the respondents had are the lack of training, statistical applications, bioinformatics support, and of rigorous guidelines and recommendations. Most of the respondents do believe there will be a technology shift from using CE only to the use of NGS on casework in 5-10 years. In addition, around 66% of respondents believe that it is moderately to very likely that the court will accept sequencing analysis. Sixteen percent fell in the middle, and the remaining 15% believe it is more unlikely, with 3% of respondents believing it is very unlikely. In conclusion, this work outlines current analytical challenges experienced by the global forensic DNA community and addresses different strategies for the implementation of next-gen sequencing technologies in casework.

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.

The frequency of the known mitochondrial variants associated with drug-induced toxicity in a Korean population

BACKGROUND

Few studies have annotated the whole mitochondrial DNA (mtDNA) genome associated with drug responses in Asian populations. This study aimed to characterize mtDNA genetic profiles, especially the distribution and frequency of well-known genetic biomarkers associated with diseases and drug-induced toxicity in a Korean population.

METHOD

Whole mitochondrial genome was sequenced for 118 Korean subjects by using a next-generation sequencing approach. The bioinformatic pipeline was constructed for variant calling, haplogroup classification and annotation of mitochondrial mutation.

RESULTS

A total of 681 variants was identified among all subjects. The MT-TRNP gene and displacement loop showed the highest numbers of variants (113 and 74 variants, respectively). The m.16189T > C allele, which is known to reduce the mtDNA copy number in human cells was detected in 25.4% of subjects. The variants (m.2706A > G, m.3010A > G, and m.1095T > C), which are associated with drug-induced toxicity, were observed with the frequency of 99.15%, 30.51%, and 0.08%, respectively. The m.2150T > A, a genotype associated with highly disruptive effects on mitochondrial ribosomes, was identified in five subjects. The D and M groups were the most dominant groups with the frequency of 34.74% and 16.1%, respectively.

CONCLUSIONS

Our finding was consistent with Korean Genome Project and well reflected the unique profile of mitochondrial haplogroup distribution. It was the first study to annotate the whole mitochondrial genome with drug-induced toxicity to predict the ADRs event in clinical implementation for Korean subjects. This approach could be extended for further study for validation of the potential ethnic-specific mitochondrial genetic biomarkers in the Korean population.

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.

Developmental Validation of a MPS Workflow with a PCR-Based Short Amplicon Whole Mitochondrial Genome Panel

For the adoption of massively parallel sequencing (MPS) systems by forensic laboratories, validation studies on specific workflows are needed to support the feasibility of implementation and the reliability of the data they produce. As such, the whole mitochondrial genome sequencing methodology-Precision ID mtDNA Whole Genome Panel, Ion Chef, Ion S5, and Converge-has been subjected to a variety of developmental validation studies. These validation studies were completed in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) validation guidelines and assessed reproducibility, repeatability, accuracy, sensitivity, specificity to human DNA, and ability to analyze challenging (e.g., mixed, degraded, or low quantity) samples. Intra- and inter-run replicates produced an average maximum pairwise difference in variant frequency of 1.2%. Concordance with data generated with traditional Sanger sequencing and an orthogonal MPS platform methodology was used to assess accuracy, and generation of complete and concordant haplotypes at DNA input levels as low as 37.5 pg of nuclear DNA or 187.5 mitochondrial genome copies illustrated the sensitivity of the system. Overall, data presented herein demonstrate that highly accurate and reproducible results were generated for a variety of sample qualities and quantities, supporting the reliability of this specific whole genome mitochondrial DNA MPS system for analysis of forensic biological evidence.

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.

Pathogenic Variant Filtering for Mitochondrial Genome Haplotype Reporting

Given the enhanced discriminatory power of the mitochondrial DNA (mtDNA) genome (mitogenome) over the commonly sequenced control region (CR) portion, the scientific merit of mitogenome sequencing is generally accepted. However, many laboratories remain beholden to CR sequencing due to privacy policies and legal requirements restricting the use of disease information or coding region (codR) information. In this report, we present an approach to obviate the reporting of sensitive codR data in forensic haplotypes. We consulted the MitoMap database to identify 92 mtDNA codR variants with confirmed pathogenicity. We determined the frequencies of these pathogenic variants in literature-quality and forensic-quality databases to be very low, at 1.2% and 0.36%, respectively. The observed effect of pathogenic variant filtering on random match statistics in 2488 forensic-quality mitogenome haplotypes from four populations was nil. We propose that pathogenic variant filtering should be incorporated into variant calling algorithms for mitogenome haplotype reporting to maximize the discriminatory power of the locus while minimizing the reveal of sensitive genetic information.

Comparison of CE- and MPS-based analyses of forensic markers in a single cell after whole genome amplification

Whole genome amplification (WGA) allows for multiple genetic analyses with low template DNA, such as DNA derived from a single cell. WGA could increase the amount of input DNA from the pg to the μg level. However, there are no studies comparing the performance of forensic markers with DNA from a single cell after WGA evaluated on both capillary electrophoresis (CE) and massively parallel sequencing (MPS) platforms. In this study, cell lines consisting of female cultured B-lymphoblastoid cells and karyocytes from male venous blood were segregated into one, two, three and five cells. Including the references with the bulk cells, all samples were generated by WGA with the multiple displacement amplification (MDA) strategy in triplicate and genotyped on CE and MPS platforms. Allele balance, stutter ratio, accuracy, repeatability and concordance of short tandem repeat (STR) markers were used to evaluate the genotyping performance on both platforms. Additionally, the sequence coverage ratio (SCR) and SNP genotypes were evaluated for sequence information generated from the MPS. Heterozygous loci showed high allele balance, with an overall average allele balance ratio larger than 0.79 on the CE and 0.75 on the MPS platforms for the venous blood cell samples; the cultured B-lymphoblastoid cell samples had ratios of 0.62 and 0.70, respectively. The stutter ratio of every source and cell number from both cell line samples were very close, ranging from 5.3%-7.2% for autosomal STRs and approximately 10 % of Y chromosomal STRs on the CE platform. The average stutter, allele, and sequence-based and length-based noise ratios were 6.6 %, 88 %, 4.7 % and 0.7 %, respectively, in the single male cell sample. SNPs also showed high consistency and intralocus balance. Our study indicated that WGA with MDA strategy works relatively well of STR and SNP genotyping with low copy number samples on CE and MPS, even with one-cell sample.

Validation of NGS for mitochondrial DNA casework at the FBI Laboratory

As a first step towards integrating next generation sequencing (NGS) technology into the FBI Laboratory's operational casework, the PowerSeq™ CRM Nested System, an NGS-based mitochondrial DNA (mtDNA) control region assay, was developmentally and internally validated. The validation studies were conducted in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) Validation Guidelines for Forensic DNA Analysis Methods, and the FBI's Quality Assurance Standards (QAS) for Forensic DNA Testing Laboratories. The assay was shown to be highly reproducible, with variant frequencies across intra and inter-run replicates of the same sample differing, on average, by just 0.3% for substitutions and point heteroplasmies and 1.5% for insertions and deletions. The assay was also shown to be extremely sensitive, yielding complete control region sequence data from as few as 2000 copies of mtDNA. This is a more than 20-fold increase in sensitivity when compared to the FBI Laboratory's current Sanger sequencing-based protocols and, based on mtDNA quantitation values of samples routinely encountered in mtDNA casework, suggests that the percentage of questioned samples from which full control region data can be recovered will increase from our current 20% to approximately 90% success with NGS technology. In addition, the assay requires on average only 30% of the extract volume typically required to develop control region profiles from degraded samples via Sanger sequencing. Overall, these studies establish the reliability of the PowerSeq™ CRM Nested System for accurate mtDNA control region typing and can serve as a model for laboratories seeking to validate NGS protocols for forensic mtDNA analysis.

Whole Mitochondrial DNA Sequencing Analysis in 47 Han Populations in Southwest China

Background

Mitochondrial DNA (mtDNA) sequencing has been used in many areas, including forensic genetics. Due to the rapid development of sequencing technology, whole mtDNA sequencing is now possible and may be used in epidemiological and forensic studies. This study aimed to use whole mtDNA sequencing to investigate 47 Chongqing Han populations in southwest China and the diversity in the mtGenome reference data.

Material/Methods

The mtDNA of 47 Chongqing Han populations was generated using the Ion Torrent Personal Genome Machine (PGM) system. The extent of the effects of the mtDNA on the subpopulations was investigated and compared with six other populations from published studies. Pairwise fixation index (FST), a measure of population differentiation due to genetic structure, were calculated. Analysis of molecular variance (AMOVA) was performed, and 1257 hypervariable region data sets were added to the principal component analysis (PCA).

Results

The whole mtDNA sequencing data of 47 southwest Chinese Han populations were successfully recovered. Expanding the sequencing rage increased the discrimination power of mtDNA from three-times to 25-times based on different populations. The subpopulation effects showed 20 times the differences in match probability when compared with south China regions.

Conclusions

Whole mtDNA sequencing distinguished between individuals from 47 Chongqing Han populations in southwest China and has potential applications that include high-quality forensic identification.

Full mtDNA genome sequencing of Brazilian admixed populations: A forensic-focused evaluation of a MPS application as an alternative to Sanger sequencing methods

The use of Massive Parallel Sequencing (MPS) techniques have been proposed by the forensic community as an alternative to Sanger sequencing methods in routine forensic casework analysis regarding mitochondrial DNA (mtDNA). Interesting features of MPS include high throughput, ability to simultaneously genotype a significant number of samples by barcoding techniques, processing automation, reduced time and costs, among others. Advantages include the capability of generating full mtDNA genome sequences versus usual techniques, usually limited to hypervariable or control regions exclusively. In this work, 96 reference single-source samples from three different Brazilian cities were subjected to full mtDNA genome sequencing by MPS techniques using an early-access version of Precision ID mtDNA Whole Genome Panel on an Ion Torrent PGM platform (Thermo Fisher Scientific, Waltham, MA, USA). Complete, high-quality sequences were obtained and sequencing performance was evaluated via four different metrics. As a subset of evaluated samples have been previously submitted for Sanger sequencing of the control region, a comparative analysis of both methods' results was conducted in order to compare technique adequacy within a forensic context. Even though this study is one of the first to report full mtDNA genome sequences for Brazilian admixed populations, the observed haplotypes exhibit a predominance of Native American and African maternal lineages in the studied sample set, reproducing results described in the literature for control regions only. Interpopulation analysis among Brazilian and 26 worldwide populations was also carried out. The results indicate that MPS-generated full mtDNA genome sequences may have great utility in forensic real casework applications, with a pronounced gain of genetic information and discrimination power provided by coding region evaluation and the enhanced capacity of heteroplasmies determination. Database construction and other relevant factors concerning implementation of such techniques in Brazilian forensic laboratories are also discussed.

Analysis and interpretation of mixture DNA using AS-PCR of mtDNA

Semi-nested PCR with allele-specific (AS) primers and sequencing of mitochondrial DNA (mtDNA) were performed to analyze and interpret DNA mixtures, especially when biological materials were degraded or contained a limited amount of DNA. SNP-STR markers were available to identify the minor DNA component using AS-PCR; moreover, SNPs in mtDNA could be used when the degraded or limited amounts of DNA mixtures were not successful with SNP-STR markers. Five pairs of allele-specific primers were designed based on three SNPs (G15043A, T16362C, and T16519C). The sequence of mtDNA control region of minor components was obtained using AS-PCR and sequencing. Sequences of the amplification fragments were aligned and compared with the sequences of known suspects or databases. When this assay was used with the T16362C and T16519C SNPs, we found it to be highly sensitive for detecting small amounts of DNA (∼30 pg) and analyzing DNA mixtures of two contributors, even at an approximately 1‰ ratio of minor and major components. An exception was tests based on the SNP G15043A, which required approximately 300 pg of a 1% DNA mixture. In simulated three contributor DNA mixtures (at rate of 1:1:1), control region fragments from each contributor were detected and interpreted. AS-PCR combined with semi-nested PCR was successfully used to identify the mtDNA control region of each contributor, providing biological evidence for excluding suspects in forensic cases, especially when biological materials were degraded or had a limited amount of DNA.

Detection of Innate and Artificial Mitochondrial DNA Heteroplasmy by Massively Parallel Sequencing: Considerations for Analysis

BACKGROUND

Mitochondrial heteroplasmy, the co-existence of different mitochondrial polymorphisms within an individual, has various forensic and clinical implications. But there is still no guideline on the application of massively parallel sequencing (MPS) in heteroplasmy detection. We present here some critical issues that should be considered in heteroplasmy studies using MPS.

METHODS

Among five samples with known innate heteroplasmies, two pairs of mixture were generated for artificial heteroplasmies with target minor allele frequencies (MAFs) ranging from 50% to 1%. Each sample was amplified by two-amplicon method and sequenced by Ion Torrent system. The outcomes of two different analysis tools, Torrent Suite Variant Caller (TVC) and mtDNA-Server (mDS), were compared.

RESULTS

All the innate heteroplasmies were detected correctly by both analysis tools. Average MAFs of artificial heteroplasmies correlated well to the target values. The detection rates were almost 90% for high-level heteroplasmies, but decreased for low-level heteroplasmies. TVC generally showed lower detection rates than mDS, which seems to be due to their own computation algorithms which drop out some reference-dominant heteroplasmies. Meanwhile, mDS reported several unintended low-level heteroplasmies which were suggested as nuclear mitochondrial DNA sequences. The average coverage depth of each sample placed on the same chip showed considerable variation. The increase of coverage depth had no effect on the detection rates.

CONCLUSION

In addition to the general accuracy of the MPS application on detecting heteroplasmy, our study indicates that the understanding of the nature of mitochondrial DNA and analysis algorithm would be crucial for appropriate interpretation of MPS results.

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.

Massively parallel sequencing-enabled mixture analysis of mitochondrial DNA samples

The mitochondrial genome has a number of characteristics that provide useful information to forensic investigations. Massively parallel sequencing (MPS) technologies offer improvements to the quantitative analysis of the mitochondrial genome, specifically the interpretation of mixed mitochondrial samples. Two-person mixtures with nuclear DNA ratios of 1:1, 5:1, 10:1, and 20:1 of individuals from different and similar phylogenetic backgrounds and three-person mixtures with nuclear DNA ratios of 1:1:1 and 5:1:1 were prepared using the Precision ID mtDNA Whole Genome Panel and Ion Chef, and sequenced on the Ion PGM or Ion S5 sequencer (Thermo Fisher Scientific, Waltham, MA, USA). These data were used to evaluate whether and to what degree MPS mixtures could be deconvolved. Analysis was effective in identifying the major contributor in each instance, while SNPs from the minor contributor's haplotype only were identified in the 1:1, 5:1, and 10:1 two-person mixtures. While the major contributor was identified from the 5:1:1 mixture, analysis of the three-person mixtures was more complex, and the mixed haplotypes could not be completely parsed. These results indicate that mixed mitochondrial DNA samples may be interpreted with the use of MPS technologies.

Performance of the Early Access AmpliSeq™ Mitochondrial Panel with degraded DNA samples using the Ion Torrent™ platform

The Early Access AmpliSeq™ Mitochondrial Panel amplifies whole mitochondrial genomes for phylogenetic and kinship identifications, using Ion Torrent™ technology. There is currently limited information on its performance with degraded DNA, a common occurrence in forensic samples. This study evaluated the performance of the Panel with DNA samples degraded in vitro, to mimic conditions commonly found in forensic investigations. Purified DNA from five individuals was heat-treated at five time points each (125°C for 0, 30, 60, 120, and 240 min; total n = 25). The quality of DNA was assessed via a real-time DNA assay of genomic DNA and prepared for massively parallel sequencing on the Ion Torrent™ platform. Mitochondrial sequences were obtained for all samples and had an amplicon coverage averaging between 66X to 2803X. Most amplicons (157/162) displayed high coverages (452 ± 333X), while reads with less than 100X coverage were recorded in five amplicons only (90 ± 5X). Amplicon coverage was decreased with prolonged heating. At 72% strand balance, reads were well balanced between forward and reverse strands. Using a coverage threshold of ten reads per SNP, complete sequences were recovered in all samples and resolved kinship and, haplogroup relations. Additionally, the HV1 and HV2 regions of the reference and 240-min heat-treated samples (n = 10) were Sanger-sequenced for concordance. Overall, this study demonstrates the efficacy of a novel forensic Panel that recovers high quality mitochondrial sequences from degraded DNA samples.

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.

Assessment of mitochondrial DNA heteroplasmy detected on commercial panel using MPS system with artificial mixture samples

Mitochondrial DNA (mtDNA) heteroplasmy is a potential genetic marker for forensic mtDNA analysis as well as phylogenic studies. Frequency of mtDNA heteroplasmy has been investigated in different populations through massively parallel sequencing (MPS) analysis, revealing various levels of frequency based on different MPS systems. For accurate heteroplasmy identification, it is essential to explore reliable detection threshold on various MPS systems. In addition, software solutions and pipelines need to be evaluated to analyze data effectively. In this study, heteroplasmy analysis was conducted on a commercially available mtDNA analysis system developed for forensic caseworks with artificially mixed DNA samples known for ratios and variant positions for assessment. mtDNA heteroplasmy > 10% was detectable with Torrent Variant Caller (TVC) while lower levels were identified using GeneMarker® HTS specialized software for minor variant detection. This study implies that analytical parameters and tools need to be optimized and evaluated for low-level heteroplasmy identification. Automated system with simple and efficient workflow is needed for forensic caseworks.

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.

Entire Mitochondrial DNA Sequencing on Massively Parallel Sequencing for the Korean Population

Mitochondrial DNA (mtDNA) genome analysis has been a potent tool in forensic practice as well as in the understanding of human phylogeny in the maternal lineage. The traditional mtDNA analysis is focused on the control region, but the introduction of massive parallel sequencing (MPS) has made the typing of the entire mtDNA genome (mtGenome) more accessible for routine analysis. The complete mtDNA information can provide large amounts of novel genetic data for diverse populations as well as improved discrimination power for identification. The genetic diversity of the mtDNA sequence in different ethnic populations has been revealed through MPS analysis, but the Korean population not only has limited MPS data for the entire mtGenome, the existing data is mainly focused on the control region. In this study, the complete mtGenome data for 186 Koreans, obtained using Ion Torrent Personal Genome Machine (PGM) technology and retrieved from rather common mtDNA haplogroups based on the control region sequence, are described. The results showed that 24 haplogroups, determined with hypervariable regions only, branched into 47 subhaplogroups, and point heteroplasmy was more frequent in the coding regions. In addition, sequence variations in the coding regions observed in this study were compared with those presented in other reports on different populations, and there were similar features observed in the sequence variants for the predominant haplogroups among East Asian populations, such as Haplogroup D and macrohaplogroups M9, G, and D. This study is expected to be the trigger for the development of Korean specific mtGenome data followed by numerous future studies.

Comparison of Two Massively Parallel Sequencing Platforms using 83 Single Nucleotide Polymorphisms for Human Identification

The potential of Massively Parallel Sequencing (MPS) technology to vastly expand the capabilities of human identification led to the emergence of different MPS platforms that use forensically relevant genetic markers. Two of the MPS platforms that are currently available are the MiSeq^® FGx™ Forensic Genomics System (Illumina) and the HID-Ion Personal Genome Machine (PGM)™ (Thermo Fisher Scientific). These are coupled with the ForenSeq™ DNA Signature Prep kit (Illumina) and the HID-Ion AmpliSeq™ Identity Panel (Thermo Fisher Scientific), respectively. In this study, we compared the genotyping performance of the two MPS systems based on 83 SNP markers that are present in both MPS marker panels. Results show that MiSeq^® FGx™ has greater sample-to-sample variation than the HID-Ion PGM™ in terms of read counts for all the 83 SNP markers. Allele coverage ratio (ACR) values show generally balanced heterozygous reads for both platforms. Two and four SNP markers from the MiSeq^® FGx™ and HID-Ion PGM™, respectively, have average ACR values lower than the recommended value of 0.67. Comparison of genotype calls showed 99.7% concordance between the two platforms.