Characterization of NIST human mitochondrial DNA SRM-2392 and SRM-2392-I standard reference materials by next generation sequencing

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.

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.

Assessment of ForenSeq mtDNA Whole Genome Kit for forensic application

Mitochondrial DNA (mtDNA) is an indispensable genetic marker in forensic genetics. The emergence and development of massively parallel sequencing (MPS) makes it possible to obtain complete mitochondrial genome sequences more quickly and accurately. The study evaluated the advantages and limitations of the ForenSeq mtDNA Whole Genome Kit in the practical application of forensic genetics by detecting human genomic DNA standards and thirty-three case samples. We used control DNA with different amount to determine sensitivity of the assay. Even when the input DNA is as low as 2.5 pg, most of the mitochondrial genome sequences could still be covered. For the detection of buccal swabs and aged case samples (bloodstains, bones, teeth), most samples could achieve complete coverage of mitochondrial genome. However, when ancient samples and hair samples without hair follicles were sequenced by the kit, it failed to obtain sequence information. In general, the ForenSeq mtDNA Whole Genome Kit has certain applicability to forensic low template and degradation samples, and these results provide the data basis for subsequent forensic applications of the assay. The overall detection process and subsequent analysis are easy to standardize, and it has certain application potential in forensic cases.

Identification of new variants in MTRNR1 and MTRNR2 genes using whole mitochondrial genome sequencing in a Taiwanese family with MERRF (myoclonic epilepsy with ragged-red fibers) syndrome

Mitochondrial encephalomyopathy is a multi-system disorder mostly caused by inborn errors of the oxidative phosphorylation (OXPHOS) system and usually manifested as complex neurological disorder and muscle weakness. Myoclonic epilepsy with ragged-red fibers (MERRF) syndrome is one of the major subtypes of mitochondrial disease associated with the m.8344A>G mutation in mitochondrial tRNALys gene. In addition to the symptoms in central nervous and muscle systems, a portion of the patients may develop hearing loss, which has been linked to the genetic mutations of mitochondrial DNA (mtDNA) especially in the mitochondrial ribosome RNA (rRNA) gene. Despite a great number of studies focusing on the consequences of mtDNA mutations, the mechanism of pathogenesis of these overt diseases has remained unclear, and there is no specific and effective treatment for MERRF syndromes. In this study, we developed a high-quality mtDNA sequencing method by next generation sequencing technology to search for the additional pathogenic variations of mtDNA from skin fibroblasts of four members in a Taiwanese family with MERRF syndrome. Through uncovering the signatures of all mtDNA variants in the MERRF family, we identified novel mtDNA variants in the genes encoding mitochondrial 12S and 16S rRNAs. The finding from this study will give us further insight into the molecular mechanisms driving the phenotypic variability and timing of onset of the MERRF syndrome.

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.

Application of a custom haplotype caller to analyze sequence-based data of 56 microhaplotypes

Microhaplotypes (microhaps) are recently introduced markers that aim to complement the limitations of conventional forensic markers such as short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs). With the potential of microhaps in forensics becoming clearer through massively parallel sequencing (MPS), MPS-based studies on microhaps are being actively reported. However, simpler workflow schemes for the generation and analysis of MPS data are still required to facilitate the practical application of MPS in forensics. In this study, we developed an in-house MPS panel that simultaneously amplifies 56 microhaps and a custom haplotype caller, Visual Microhap. The developed tool works on a web browser and provides four analysis options to extract SNP-based haplotypes from sequence-based data obtained by STRait Razor 3.0. To demonstrate the utility of the MPS panel and data analysis workflow scheme, we also analyzed 56 microhaps of 286 samples from four populations (African-American, Caucasian, Hispanic, and Korean). The average effective number of alleles (A_e) for the four groups was 3.45, ranging from 1.74 to 6.98. Forensic statistical parameters showed that this microhap panel is more powerful than conventional autosomal STRs for human identification. Meanwhile, the 56-plex panel mostly comprised microhaps with high Ae; however, the four populations were grossly distinguishable from each other by cluster analysis. Consequently, the developed in-house MPS panel for 56 microhaps and the adopted workflow using open-source tools can increase the utility of microhap MPS in forensic research and practice.

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.

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.

Human Mitochondrial Control Region and mtGenome: Design and Forensic Validation of NGS Multiplexes, Sequencing and Analytical Software

Forensic mitochondrial DNA (mtDNA) analysis conducted using next-generation sequencing (NGS), also known as massively parallel sequencing (MPS), as compared to Sanger-type sequencing brings modern advantages, such as deep coverage per base (herein referred to as read depth per base pair (bp)), simultaneous sequencing of multiple samples (libraries) and increased operational efficiencies. This report describes the design and developmental validation, according to forensic quality assurance standards, of end-to-end workflows for two multiplexes, comprised of ForenSeq mtDNA control region and mtDNA whole-genome kits the MiSeq FGx^TM instrument and ForenSeq universal analysis software (UAS) 2.0/2.1. Polymerase chain reaction (PCR) enrichment and a tiled amplicon approach target small, overlapping amplicons (60–150 bp and 60–209 bp for the control region and mtGenome, respectively). The system provides convenient access to data files that can be used outside of the UAS if desired. Studies assessed a range of environmental and situational variables, including but not limited to buccal samples, rootless hairs, dental and skeletal remains, concordance of control region typing between the two multiplexes and as compared to orthogonal data, assorted sensitivity studies, two-person DNA mixtures and PCR-based performance testing. Limitations of the system and implementation considerations are discussed. Data indicated that the two mtDNA multiplexes, MiSeq FGx and ForenSeq software, meet or exceed forensic DNA quality assurance (QA) guidelines with robust, reproducible performance on samples of various quantities and qualities.

Assessment of Illumina® Human mtDNA Genome assay: workflow evaluation with development of analysis and interpretation guidelines

Mitochondrial DNA (mtDNA) is a small but significant part of the human genome, whose applicability potential has gradually increased with the advent of massively parallel sequencing (MPS) technology. Knowledge of the particular workflow, equipment, and reagents used, along with extensive usage of negative controls to monitor all preparation steps constitute the prerequisites for confident reporting of results. In this study, we performed an assessment of Illumina® Human mtDNA Genome assay on MiSeq FGx™ instrument. Through analysis of several types of negative controls, as well as mtDNA positive controls, we established thresholds for data analysis and interpretation, consisting of several components: minimum read depth (220 reads), minimum quality score (41), percentage of minor allele sufficient for analysis (3.0%), percentage of minor allele sufficient for interpretation (6.0%), and percentage of major allele sufficient for homoplasmic variant call (97.0%). Based on these criteria, we defined internal guidelines for analysis and interpretation of mtDNA results obtained by MPS. Our study shows that the whole mtDNA assay on MiSeq FGx™ produces repeatable and reproducible results, independent of the analyst, which are also concordant with Sanger-type sequencing results for mtDNA control region, as well as with MPS results produced by NextSeq®. Overall, established thresholds and interpretation guidelines were successfully applied for the sequencing of complete mitochondrial genomes from high-quality samples. The underlying principles and proposed methodology on the definition of internal laboratory guidelines for analysis and interpretation of MPS results may be applicable to similar MPS workflows, e.g. targeting good-quality samples in forensic genetics and molecular diagnostics.

Association Study Between Genetic Variation in Whole Mitochondrial Genome and Ischemic Stroke

Mitochondrial DNA (mtDNA) affects the mitochondrial function, which is potentially related to susceptibility to ischemic stroke (IS). However, study on IS genetics by whole mitochondrial genome sequencing has not been extensively explored. Therefore, a two-stage study was designed to explore the relationship between the whole mitochondrial genome variants and IS. In the first stage, whole mitochondrial genomes of 52 IS patients and 55 controls were sequenced by next-generation sequencing. Fifty-three mtDNA mutation sites which may be related to the pathogenesis of IS were discovered. Nine unreported mtDNA variation sites were found for the first time. In the second larger Chinese cohort, we confirmed that m.T195C and m.T12338C in the mitochondrial D-loop region were the protective factors of IS, especially m.T195C and m.C311T in the LAA subtype. In conclusion, our study provided population genetic information and a reference for IS-relevant research, with wide applications in diagnosis, therapeutic treatments and prediction of IS.

STRait Razor Online: An enhanced user interface to facilitate interpretation of MPS data

Since 2013, STRait Razor has enabled analysis of massively parallel sequencing (MPS) data from various marker systems such as short tandem repeats, single nucleotide polymorphisms, insertion/deletions, and mitochondrial DNA. In this paper, STRait Razor Online (SRO), available at https://www.unthsc.edu/straitrazor, is introduced as an interactive, Shiny-based user interface for primary analysis of MPS data and secondary analysis of STRait Razor haplotype pileups. This software can be accessed from any common browser via desktop, tablet, or smartphone device. SRO is available also as a standalone application and open-source R script available at https://github.com/ExpectationsManaged/STRaitRazorOnline. The local application is capable of batch processing of both fastq files and primary analysis output. Processed batches generate individual report folders and summary reports at the locus- and haplotype-level in a matter of minutes. For example, the processing of data from ∼700 samples generated with the ForenSeq Signature Preparation Kit from allsequences.txt to a final table can be performed in ∼40 min whereas the Excel-based workbooks can take 35-60 h to compile a subset of the tables generated by SRO. To facilitate analysis of single-source, reference samples, a preliminary triaging system was implemented that calls potential alleles and flags loci suspected of severe heterozygote imbalance. When compared to published, manually curated data sets, 98.72 % of software-assigned allele calls without manual interpretation were consistent with curated data sets, 0.99 % loci were presented to the user for interpretation due to heterozygote imbalance, and the remaining 0.29 % of loci were inconsistent due to the analytical thresholds used across the studies.

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.

The lot-to-lot variability in the mitochondrial genome of controls

Current research in the biomedical field has illustrated how cell lines used as reference standards can change over time and, more importantly, can affect research and diagnostic results obtained from these cell lines. With the use of increasingly sensitive and highly resolving technologies (e.g., massively parallel sequencing), forensic scientists must be aware of and account for potential variability in the cell lines used as controls in their validation studies and day-to-day casework. In this study, multiple lot numbers from four commonly-used control cell line DNAs were sequenced with massively parallel sequencing on the Ion S5. The variability among these different lots was evaluated, and the effect on forensic laboratory work discussed.

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.

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.