Detection Of Heteroplasmy Research Articles

The quality and detection limits of mitochondrial heteroplasmy by long read nanopore sequencing

This study evaluates long-read and short-read sequencing for mitochondrial DNA (mtDNA) heteroplasmy detection. 592,315 bootstrapped datasets generated from two single-nucleotide mismatched ultra-deep sequenced mtDNA samples were used to assess basecalling error and accuracy, limit of heteroplasmy detection, and heteroplasmy detection across various coverage depths. Results showed high Phred scores of data with GC-rich sequence bias for long reads. Limit of detection of 12% heteroplasmy was identified, showing strong correlation (R2 ≥ 0.955) with expected heteroplasmy but underreporting tendency of high-level variants. Nanopore sequencing shows potential for direct applicability in mitochondrial diseases diagnostics, but stringent validation processes to ensure diagnostic result quality are required.

Low-level point heteroplasmy detection in human mitogenomes amplified with different polymerases and sequenced on MiSeq FGx platform.

Massively parallel sequencing of mitogenomes usually requires prior amplification. The PCR step may influence the quality of the data obtained, especially when low-level heteroplasmy detection is applied. The aim of this study was to compare the reliability of two different DNA polymerases in detecting homoplasmic and heteroplasmic substitutions in human mitogenomes. Mitogenomes of five samples were amplified with Long PCR Enzyme Mix from Fermentas or TaKaRa LA Taq DNA Polymerase from TaKaRa. Then, NexteraTM XT DNA libraries were sequenced on MiSeq FGx platform (Illumina). mtDNA substitutions were called for alternative variants above the 1% level. All homoplasmic substitutions detected in amplicons generated with polymerases studied here and sequenced on MiSeq FGx system were consistently identified as homoplasmies with alternative sequencing methods. TaKaRa LA Taq DNA Polymerase was found to be less accurate in low-level heteroplasmy detection than Long PCR Enzyme Mix enzyme as more false negative and false positive results were observed for minority variants called above the 1% level. Nevertheless, both PCR systems studied can be successfully used to detect authentic mtDNA substitutions, for which minority variants exceed the 3.61% level assuming at least 10,000x coverage and sequencing Nextera XT DNA libraries on MiSeq FGx machine. The accuracy and sensitivity of point heteroplasmy detection with the MiSeq FGx instrument varies on polymerase used for mtDNA amplification. Therefore, it is recommended to validate the laboratory protocols used for mtDNA substitution detection prior to their implementation for the forensic or medical genetics purposes.

PHFinder: assisted detection of point heteroplasmy in Sanger sequencing chromatograms.

Heteroplasmy is the presence of two or more organellar genomes (mitochondrial or plastid DNA) in an organism, tissue, cell or organelle. Heteroplasmy can be detected by visual inspection of Sanger sequencing chromatograms, where it appears as multiple peaks of fluorescence at a single nucleotide position. Visual inspection of chromatograms is both consuming and highly subjective, as heteroplasmy is difficult to differentiate from background noise. Few software solutions are available to automate the detection of point heteroplasmies, and those that are available are typically proprietary, lack customization or are unsuitable for automated heteroplasmy assessment in large datasets. Here, we present PHFinder, a Python-based, open-source tool to assist in the detection of point heteroplasmies in large numbers of Sanger chromatograms. PHFinder automatically identifies point heteroplasmies directly from the chromatogram trace data. The program was tested with Sanger sequencing data from 100 humpback whales (Megaptera novaeangliae) tissue samples with known heteroplasmies. PHFinder detected most (90%) of the known heteroplasmies thereby greatly reducing the amount of visual inspection required. PHFinder is flexible and enables explicit specification of key parameters to infer double peaks (i.e., heteroplasmies).

CmVCall: An automated and adjustable nanopore analysis pipeline for heteroplasmy detection of the control region in human mitochondrial genome

Genetic associations between human mitochondrial DNA (mtDNA) heteroplasmy and mitochondrial diseases, aging, and cancer have been elaborated, contributing a lot to the further understanding of mtDNA polymorphic spectrum in anthropology, population, and forensic genetics. In the past decade, heteroplasmy detection using Sanger sequencing and next generation sequencing (NGS) was hampered by the former’s inefficiency and the latter’s inherent bias due to amplification and mapping of short reads, respectively. Nanopore sequencing stands out for its ability to yield long contiguous segments of DNA, providing a new insight into heterogeneity authentication. In addition to MinION from Oxford Nanopore Technologies, an alternative nanopore sequencer QNome (Qitan Technology) has also been applied to various biological research and the forensic applicability of this platform has been proved recently. In this study, we evaluated the performance of four commonly used variant callers in the heterogeneity authentication of the control region of human mtDNA based on simulations of different ratios generated by mixing QNome nanopore sequencing reads of two synthetic sequences. Then, an open-source and python-based nanopore analytics pipeline, CmVCall was developed and incorporated multiple programs including reads filtering, removal of nuclear mitochondrial sequences (NUMTs), alignment, optional ‘Correction’ mode, and heterogeneity identification. CmVCall can achieve high precision, accuracy, and recall of 100%, 99.9%, and 92.3% with a 5% heteroplasmy level in ‘Correction’ mode. Moreover, blood, saliva, and hair shaft samples from monozygotic (MZ) twins were used for heterogeneity evaluation and comparison with the NGS data. Results of MZ twin samples showed that CmVCall could identify more point heteroplasmy sites, revealing significant levels of inter- and intra-individual mtDNA polymorphism. In conclusion, we believe that this analysis pipeline will lay a solid foundation for the development of a comprehensive nanopore analysis pipeline targeting the whole mitochondrial genome.

Nuclear DNA segments homologous to mitochondrial DNA are obstacles for detecting heteroplasmy in sugar beet (Beta vulgaris L.).

Heteroplasmy, the coexistence of multiple mitochondrial DNA (mtDNA) sequences in a cell, is well documented in plants. Next-generation sequencing technology (NGS) has made it feasible to sequence entire genomes. Thus, NGS has the potential to detect heteroplasmy; however, the methods and pitfalls in heteroplasmy detection have not been fully investigated and identified. One obstacle for heteroplasmy detection is the sequence homology between mitochondrial-, plastid-, and nuclear DNA, of which the influence of nuclear DNA segments homologous to mtDNA (numt) need to be minimized. To detect heteroplasmy, we first excluded nuclear DNA sequences of sugar beet (Beta vulgaris) line EL10 from the sugar beet mtDNA sequence. NGS reads were obtained from single plants of sugar beet lines NK-195BRmm-O and NK-291BRmm-O and mapped to the unexcluded mtDNA regions. More than 1000 sites exhibited intra-individual polymorphism as detected by genome browsing analysis. We focused on a 309-bp region where 12 intra-individual polymorphic sites were closely linked to each other. Although the existence of DNA molecules having variant alleles at the 12 sites was confirmed by PCR amplification from NK-195BRmm-O and NK-291BRmm-O, these variants were not always called by six variant-calling programs, suggesting that these programs are inappropriate for intra-individual polymorphism detection. When we changed the nuclear DNA reference, a numt absent from EL10 was found to include the 309-bp region. Genetic segregation of an F2 population from NK-195BRmm-O x NK-291BRmm-O supported the numt origin of the variant alleles. Using four references, we found that numt detection exhibited reference dependency, and extreme polymorphism of numts exists among sugar beet lines. One of the identified numts absent from EL10 is also associated with another intra-individual polymorphic site in NK-195mm-O. Our data suggest that polymorphism among numts is unexpectedly high within sugar beets, leading to confusion about the true degree of heteroplasmy.

Capability of the iSeq 100 sequencing system from Illumina to detect low-level substitutions in the human mitochondrial genome

The significance of mtDNA heteroplasmy in forensic and medical genetics has increased recently because massively parallel sequencing (MPS) technologies enable more accurate and precise detection of minority nucleotide variants. Recent reports have shown that detection of low-level substitutions may depend on library preparation or sequencing protocol, and can vary for different MPS platforms. The MiSeq (Illumina) and Ion S5 (Thermo Fisher Scientific) are mainly used for heteroplasmy detection, but no data are available regarding the iSeq 100, an Illumina platform of the smallest throughput. Notably, unlike the other systems, the machine utilizes sequencing by synthesis one-channel chemistry to determine DNA sequences. Thus, it is important to verify the capability of the iSeq 100 system to determine mitochondrial haplotypes and detect heteroplasmic substitutions. In this study, previously determined entire mitochondrial genomes were sequenced with the iSeq 100 system. Each mitogenome was sequenced twice, giving approximately 2000x and 10,000x coverage. All homoplasmic mutations and minority variants above the 19 % level detected with the iSeq 100 system were also observed after dideoxy sequencing. Moreover, all heteroplasmic substitutions above the 2 % level were consistently detected with SBS one-channel chemistry. However, detection of low-level mtDNA variants may require additional, confirmatory experiments. In summary, the iSeq 100 system enables reproducible and accurate sequencing of human mitochondrial genomes. Detection of mtDNA minority variants depends on the laboratory protocol and sequencing platform used, but homoplasmic mutations and heteroplasmy above the 2 % level can be correctly detected with the iSeq 100 system.

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.

First mitochondrial genome-wide association study with metabolomics.

In the era of personalized medicine with more and more patient-specific targeted therapies being used, we need reliable, dynamic, faster and sensitive biomarkers both to track the causes of disease and to develop and evolve therapies during the course of treatment. Metabolomics recently has shown substantial evidence to support its emerging role in disease diagnosis and prognosis. Aside from biomarkers and development of therapies, it is also an important goal to understand the involvement of mitochondrial DNA (mtDNA) in metabolic regulation, aging and disease development. Somatic mutations of the mitochondrial genome are also heavily implicated in age-related disease and aging. The general hypothesis is that an alteration in the concentration of metabolite profiles (possibly conveyed by lifestyle and environmental factors) influences the increase of mutation rate in the mtDNA and thereby contributes to a range of pathophysiological alterations observed in complex diseases. We performed an inverted mitochondrial genome-wide association analysis between mitochondrial nucleotide variants (mtSNVs) and concentration of metabolites. We used 151 metabolites and the whole sequenced mitochondrial genome from 2718 individuals to identify the genetic variants associated with metabolite profiles. Because of the high coverage, next-generation sequencing-based analysis of the mitochondrial genome allows for an accurate detection of mitochondrial heteroplasmy and for the identification of variants associated with the metabolome. The strongest association was found for mt715G > A located in the MT-12SrRNA with the metabolite ratio of C2/C10:1 (P-value = 6.82*10−09, β = 0.909). The second most significant mtSNV was found for mt3714A > G located in the MT-ND1 with the metabolite ratio of phosphatidylcholine (PC) ae C42:5/PC ae C44:5 (P-value = 1.02*10−08, β = 3.631). A large number of significant metabolite ratios were observed involving PC aa C36:6 and the variant mt10689G > A, located in the MT-ND4L gene. These results show an important interconnection between mitochondria and metabolite concentrations. Considering that some of the significant metabolites found in this study have been previously related to complex diseases, such as neurological disorders and metabolic conditions, these associations found here might play a crucial role for further investigations of such complex diseases. Understanding the mechanisms that control human health and disease, in particular, the role of genetic predispositions and their interaction with environmental factors is a prerequisite for the development of safe and efficient therapies for complex disorders.

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.

Validation of whole genome sequencing from dried blood spots

BackgroundDried blood spots (DBS) are a relatively inexpensive source of nucleic acids and are easy to collect, transport, and store in large-scale field surveys, especially in resource-limited settings. However, their performance in whole-genome sequencing (WGS) relative to that of venous blood DNA has not been analyzed for various downstream applications.MethodsThis study compares the WGS performance of DBS paired with venous blood samples collected from 12 subjects.ResultsResults of standard quality checks of coverage, base quality, and mapping quality were found to be near identical between DBS and venous blood. Concordance for single-nucleotide variants, insertions and deletions, and copy number variants was high between these two sample types. Additionally, downstream analyses typical of population-based studies were performed, such as mitochondrial heteroplasmy detection, haplotype analysis, mitochondrial copy number changes, and determination of telomere lengths. The absolute mitochondrial copy number values were higher for DBS than for venous blood, though the trend in sample-to-sample variation was similar between DBS and blood. Telomere length estimates in most DBS samples were on par with those from venous blood.ConclusionDBS samples can serve as a robust and feasible alternative to venous blood for studies requiring WGS analysis.

Validation and clinical performance of a combined nuclear-mitochondrial next-generation sequencing and copy number variant analysis panel in a Canadian population.

Diagnosing mitochondrial disorders is a challenge due to the heterogeneous clinical presentation and large number of associated genes. A custom next generation sequencing (NGS) panel was developed incorporating the full mitochondrial genome (mtDNA) plus 19 nuclear genes involved in structural mitochondrial defects and mtDNA maintenance. This assay is capable of simultaneously detecting small gene sequence variations and larger copy number variants (CNVs) in both the nuclear and mitochondrial components along with heteroplasmy detection down to 5%. We describe technical validations of this panel and its implementation for clinical testing in a Canadian reference laboratory, and report its clinical performance in the initial 950 patients tested. Using this assay, we demonstrate a diagnostic yield of 18.1% of patients with known pathogenic variants. In addition to the common 5 kb mtDNA deletion, we describe significant contribution of pathogenic CNVs in both the mitochondrial genome and nuclear genes in this patient population.

Approaches to Whole Mitochondrial Genome Sequencing on the Oxford Nanopore MinION.

Traditional approaches for interrogating the mitochondrial genome often involve laborious extraction and enrichment protocols followed by Sanger sequencing. Although preparation techniques are still demanding, the advent of next-generation or massively parallel sequencing has made it possible to routinely obtain nucleotide-level data with relative ease. These short-read sequencing platforms offer deep coverage with unparalleled read accuracy in high-complexity genomic regions but encounter numerous difficulties in the low-complexity homopolymeric sequences characteristic of the mitochondrial genome. The inability to discern identical units within monomeric repeats and resolve copy-number variations for heteroplasmy detection results in suboptimal genome assemblies that ultimately complicate downstream data analysis and interpretation of biological significance. Oxford Nanopore Technologies offers the ability to generate long-read sequencing data on a pocket-sized device known as the MinION. Nanopore-based sequencing is scalable, portable, and theoretically capable of sequencing the entire mitochondrial genome in a single contig. Furthermore, the recent development of a nanopore protein with dual reader heads allows for clear identification of nucleotides within homopolymeric stretches, significantly increasing resolution throughout these regions. The unrestricted read lengths, superior homopolymeric resolution, and affordability of the MinION device make it an attractive alternative to the labor-intensive, time-consuming, and costly mainstay deep-sequencing platforms. This article describes three approaches to extract, prepare, and sequence mitochondrial DNA on the Oxford Nanopore MinION device. Two of the workflows include enrichment of mitochondrial DNA prior to sequencing, whereas the other relies on direct sequencing of native genomic DNA to allow for simultaneous assessment of the nuclear and mitochondrial genomes. © 2019 by John Wiley & Sons, Inc. Basic Protocol: Enrichment-free mitochondrial DNA sequencing Alternate Protocol 1: Mitochondrial DNA sequencing following enrichment with polymerase chain reaction (PCR) Alternate Protocol 2: Mitochondrial DNA sequencing following enrichment with PCR-free hybridization capture Support Protocol 1: DNA quantification and quality assessment using the Agilent 4200 TapeStation System Support Protocol 2: AMPure XP bead clean-up Support Protocol 3: Suggested data analysis pipeline.

Impact of the sequencing method on the detection and interpretation of mitochondrial DNA length heteroplasmy

Advancements in sequencing technologies allow for rapid and efficient analysis of mitochondrial DNA (mtDNA) in forensic laboratories, which is particularly beneficial for specimens with limited nuclear DNA. Next generation sequencing (NGS) offers higher throughput and sensitivity over traditional Sanger-type sequencing (STS) as well as the ability to quantitatively analyze the data. Changes in sample preparation, sequencing method and analysis required for NGS may alter the mtDNA haplotypes compared to previously generated STS data. Thus, the present study aimed to characterize the impact of different sequencing workflows on the detection and interpretation of length heteroplasmy (LHP), a particularly complicated aspect of mtDNA analysis. Whole mtDNA genome (mitogenome) data were generated for 16 high-quality samples using well-established Illumina and Ion methods, and the NGS data were compared to previously-generated STS mtDNA control region data. Although the mitogenome haplotypes were concordant with the exception of length and low-level variants (<30 % variant frequency), LHP in the hypervariable segment (HVS) polycytosine regions (C-tracts) differed across sequencing methods. Consistent with previous studies, LHP in HVS1 was observed in samples with nine or more consecutive cytosines (Cs) and eight Cs in the HVS2 region in the STS data. The Illumina data produced a similar pattern of LHP as the STS data, whereas the Ion data were noticeably different. More complex LHP (i.e. more length molecules) was observed in the Ion data, as length variation occurred in multiple homopolymer stretches within the targeted HVS regions. Further, the STS dominant or major molecule (MM) differed from the Ion MM in 11 (37 %) of the 30 regions evaluated and six instances (20 %) in Illumina data. This is of particular interest, as the MM is used by many forensic laboratories to report the HVS C-tract in the mtDNA haplotype. In general, the STS MMs were longer than the Illumina MMs, while the Ion MMs were the shortest. The higher rate of homopolymer indels in Ion data likely contributed to these differences. Supplemental analysis with alternative approaches demonstrated that the LHP pattern may also be altered by the bioinformatic tool and workflow used for data interpretation. The broader application of NGS in forensic laboratories will undoubtedly result in the use of varying sample preparation and sequencing methods. Based on these findings, minor LHP differences are expected across sequencing workflows, and it will be important that C-tract indels continue to be ignored for forensic queries and comparisons.

Sensitivity of mitochondrial DNA heteroplasmy detection using Next Generation Sequencing

Although the use of Next Generation Sequencing (NGS) in mitochondrial DNA (mtDNA) studies related to forensic and human genetics has contributed to the report of heteroplasmy at very low levels (lower than 1% and even 0.5%), their detection is not a straightforward process. Our purpose is to establish mitochondrial heteroplasmy detection limits, generating mixed bases at low frequencies by the PCR co-amplification of mtDNA and a nuclear insertion of mitochondrial origin (NUMT). A primer set that presents a perfect annealing with both mitochondrial and nuclear DNA was used to amplify the mitochondrial DNA region located between positions 6739 and 8910 and the corresponding region located inside a NUMT located in chromosome 1 (GRCh38.p12 Primary Assembly 631910-634079) that presents an identity of 98% with the corresponding region of mtDNA in two samples. Independent Nextera XT® (Illumina) NGS libraries were generated for each sample and sequenced in different MiSeq (Illumina) runs. Non-identical and identical positions between individual mtDNA and NUMT sequences were detected, and heteroplasmy detections limits were established: a) with a minor allele frequency <1.5%, false positive and negative can occur, even with a depth higher than 3000X; b) with a minor allele frequency >3%, no presence of false positive and negative were detected with a depth of ~1000X; and c) with a minor allele frequency between 1.5% and 3%, a minimal depth of 3000X was necessary to avoid false positives or negatives detection. Our results demonstrate an alternative strategy to establish a reliable limit of mitochondrial heteroplasmy detection.

Mitochondrial DNA data analysis strategies that inform MPS-based forensic casework implementation

Massively parallel sequencing (MPS) of human mitochondrial control region or whole genome is now faster, easier and more informative than was possible with Sanger sequencing and capillary electrophoresis. This study focuses on the aspects of mtDNA data that are important to consider when implementing MPS-based mtDNA typing into casework operations. This includes evaluation of data from the ForenSeq mtDNA typing system that indicate how interpretation concepts integral to Sanger sequencing of mtDNA apply directly to MPS-based mtDNA typing. In addition, there are new interpretation considerations that may be helpful when determining standard operating procedures and casework interpretation guidelines. Sequencing data are presented that may assist with SOP determination relative to sensitivity assessments that may help inform interpretation criteria, including thresholds, relative to input DNA template concentration ranges, number of pooled samples for MiSeq FGx sequencing runs and downstream bioinformatic demultiplexing. Given the quantitative nature of digital MPS data, these limit of detection studies may also assist when laboratories wish to set criteria for heteroplasmy detection (within an individual) or for mixed DNA sample interpretation. Assessment of variables will be provided across sample types such as hairs and PCR/sequencing inhibition when typing skeletal and dental remains to assist with fundamental understanding of ForenSeq mtDNA typing system limitations as well as metrics that may assist when making data quality determinations.

Evaluation of mitogenome sequence concordance, heteroplasmy detection, and haplogrouping in a worldwide lineage study using the Precision ID mtDNA Whole Genome Panel

The emergence of Massively Parallel Sequencing technologies enabled the analysis of full mitochondrial (mt)DNA sequences from forensically relevant samples that have, so far, only been typed in the control region or its hypervariable segments. In this study, we evaluated the performance of a commercially available multiplex-PCR-based assay, the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific), for the amplification and sequencing of the entire mitochondrial genome (mitogenome) from even degraded forensic specimens. For this purpose, more than 500 samples from 24 different populations were selected to cover the vast majority of established superhaplogroups. These are known to harbor different signature sequence motifs corresponding to their phylogenetic background that could have an effect on primer binding and, thus, could limit a broad application of this molecular genetic tool. The selected samples derived from various forensically relevant tissue sources and were DNA extracted using different methods. We evaluated sequence concordance and heteroplasmy detection and compared the findings to conventional Sanger sequencing as well as an orthogonal MPS platform. We discuss advantages and limitations of this approach with respect to forensic genetic workflow and analytical requirements.

Mitochondrial DNA enrichment reduced NUMT contamination in porcine NGS analyses.

Genetic associations between mitochondrial DNA (mtDNA) and economic traits have been widely reported for pigs, which indicate the importance of mtDNA. However, studies on mtDNA heteroplasmy in pigs are rare. Next generation sequencing (NGS) methodologies have emerged as a promising genomic approach for detection of mitochondrial heteroplasmy. Due to the short reads, flexible bioinformatic analyses and the contamination of nuclear mitochondrial sequences (NUMTs), NGS was expected to increase false-positive detection of heteroplasmy. In this study, Sanger sequencing was performed as a gold standard to detect heteroplasmy with a detection sensitivity of 5% in pigs and then one whole-genome sequencing method (WGS) and two mtDNA enrichment sequencing methods (Capture and LongPCR) were carried out. The aim of this study was to determine whether mitochondrial heteroplasmy identification from NGS data was affected by NUMTs. We find that WGS generated more false intra-individual polymorphisms and less mapping specificity than the two enrichment sequencing methods, suggesting NUMTs indeed led to false-positive mitochondrial heteroplasmies from NGS data. In addition, to accurately detect mitochondrial diversity, three commonly used tools-SAMtools, VarScan and GATK-with different parameter values were compared. VarScan achieved the best specificity and sensitivity when considering the base alignment quality re-computation and the minimum variant frequency of 0.25. It also suggested bioinformatic workflow interfere in the identification of mtDNA SNPs. In conclusion, intra-individual polymorphism in pig mitochondria from NGS data was confused with NUMTs, and mtDNA-specific enrichment is essential before high-throughput sequencing in the detection of mitochondrial genome sequences.

Mitochondrial and nuclear disease panel (Mito-aND-Panel): Combined sequencing of mitochondrial and nuclear DNA by a cost-effective and sensitive NGS-based method.

BackgroundThe diagnosis of mitochondrial disorders is challenging because of the clinical variability and genetic heterogeneity of these conditions. Next‐Generation Sequencing (NGS) technology offers a robust high‐throughput platform for nuclear and mitochondrial DNA (mtDNA) analyses.MethodWe developed a custom Agilent SureSelect Mitochondrial and Nuclear Disease Panel (Mito‐aND‐Panel) capture kit that allows parallel enrichment for subsequent NGS‐based sequence analysis of nuclear mitochondrial disease‐related genes and the complete mtDNA genome. Sequencing of enriched mtDNA simultaneously with nuclear genes was compared with the separated sequencing of the mitochondrial genome and whole exome sequencing (WES).ResultsThe Mito‐aND‐Panel permits accurate detection of low‐level mtDNA heteroplasmy due to a very high sequencing depth compared to standard diagnostic procedures using Sanger sequencing/SNaPshot and WES which is crucial to identify maternally inherited mitochondrial disorders.ConclusionWe established a NGS‐based method with combined sequencing of the complete mtDNA and nuclear genes which enables a more sensitive heteroplasmy detection of mtDNA mutations compared to traditional methods. Because the method promotes the analysis of mtDNA variants in large cohorts, it is cost‐effective and simple to setup, we anticipate this is a highly relevant method for sequence‐based genetic diagnosis in clinical diagnostic applications.

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.

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

BackgroundMitochondrial 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.MethodsAmong 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.ResultsAll 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.ConclusionIn 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.