PCR-Free Enrichment of Mitochondrial DNA from Human Blood and Cell Lines for High Quality Next-Generation DNA Sequencing

CRISPR Screening in Tandem with Targeted mtDNA Damage Reveals WRNIP1 Essentiality

A major impediment to the characterization of mtDNA repair mechanisms, in comparison to nuclear DNA repair mechanisms, is the difficulty of specifically addressing mitochondrial damage. Using a mitochondria-penetrating peptide, we can deliver DNA-damaging agents directly to mitochondria, bypassing the nuclear compartment. Here, we describe the use of a mtDNA-damaging agent in tandem with CRISPR/Cas9 screening for the genome-wide discovery of factors essential for mtDNA damage response. Using mitochondria-targeted doxorubicin (mtDox) we generate mtDNA double-strand breaks (mtDSBs) specifically in this organelle. Combined with an untargeted Dox screen, we identify genes with significantly greater essentiality during mitochondrial versus nuclear DNA damage. We characterize the essentially of our top hit - WRNIP1 - observed here for the first time to respond to mtDNA damage. We further investigate the mitochondrial role of WRNIP1 in innate immune signaling and nuclear genome maintenance, outlining a model that experimentally supports mitochondrial turnover in response to mtDSBs.

Enhanced mitochondrial genome analysis: bioinformatic and long-read sequencing advances and their diagnostic implications

INTRODUCTION

Primary mitochondrial diseases (PMDs) comprise a large and heterogeneous group of genetic diseases that result from pathogenic variants in either nuclear DNA (nDNA) or mitochondrial DNA (mtDNA). Widespread adoption of next-generation sequencing (NGS) has improved the efficiency and accuracy of mtDNA diagnoses; however, several challenges remain.

AREAS COVERED

In this review, we briefly summarize the current state of the art in molecular diagnostics for mtDNA and consider the implications of improved whole genome sequencing (WGS), bioinformatic techniques, and the adoption of long-read sequencing, for PMD diagnostics.

EXPERT OPINION

We anticipate that the application of PCR-free WGS from blood DNA will increase in diagnostic laboratories, while for adults with myopathic presentations, WGS from muscle DNA may become more widespread. Improved bioinformatic strategies will enhance WGS data interrogation, with more accurate delineation of mtDNA and NUMTs (nuclear mitochondrial DNA segments) in WGS data, superior coverage uniformity, indirect measurement of mtDNA copy number, and more accurate interpretation of heteroplasmic large-scale rearrangements (LSRs). Separately, the adoption of diagnostic long-read sequencing could offer greater resolution of complex LSRs and the opportunity to phase heteroplasmic variants.

The Mighty NUMT: Mitochondrial DNA Flexing Its Code in the Nuclear Genome

Nuclear-mitochondrial DNA segments (NUMTs) are mitochondrial DNA (mtDNA) fragments that have been inserted into the nuclear genome. Some NUMTs are common within the human population but most NUMTs are rare and specific to individuals. NUMTs range in size from 24 base pairs to encompassing nearly the entire mtDNA and are found throughout the nuclear genome. Emerging evidence suggests that the formation of NUMTs is an ongoing process in humans. NUMTs contaminate sequencing results of the mtDNA by introducing false positive variants, particularly heteroplasmic variants present at a low variant allele frequency (VAF). In our review, we discuss the prevalence of NUMTs in the human population, the potential mechanisms of de novo NUMT insertion via DNA repair mechanisms, and provide an overview of the existing approaches for minimizing NUMT contamination. Apart from filtering known NUMTs, both wet lab-based and computational methods can be used to minimize the contamination of NUMTs in analyses of human mtDNA. Current approaches include: (1) isolating mitochondria to enrich for mtDNA; (2) applying basic local alignment to identify NUMTs for subsequent filtering; (3) bioinformatic pipelines for NUMT detection; (4) k-mer-based NUMT detection; and (5) filtering candidate false positive variants by mtDNA copy number, VAF, or sequence quality score. Multiple approaches must be applied in order to effectively identify NUMTs in samples. Although next-generation sequencing is revolutionizing our understanding of heteroplasmic mtDNA, it also raises new challenges with the high prevalence and individual-specific NUMTs that need to be handled with care in studies of mitochondrial genetics.

Heteroplasmic Mutant Load Differences in Mitochondrial DNA-Associated Leigh Syndrome

BACKGROUND

Mitochondrial DNA (mtDNA)-associated Leigh syndrome is influenced by mutant pathogenicity and corresponding heteroplasmic loads; however, the manner in which heteroplasmic mutant load affects patient phenotypes and the relationship between mutant types and heteroplasmic mutant loads remain unknown. We aimed to investigate the distribution of the mutant load of different mtDNA mutations in a single-center cohort.

METHODS

We used next-generation sequencing to confirm mtDNA mutations in 31 patients with Leigh syndrome. Subsequently, we counted the number of mtDNA reads to quantitatively analyze the heteroplasmic mutant load and categorize the patients according to the mtDNA mutations they harbored. Confirmed cases of mtDNA-associated Leigh syndrome were classified according to the mutations observed in six genes and 10 nucleotides.

RESULTS

Of the 31 patients with Leigh syndrome, 27 harbored known pathogenic mutations. We discovered that MT-ATP6 was the most commonly mutated gene (n = 13 patients), followed by MT-ND3 (n = 7) and MT-ND5 (n = 4). MT-ATP6 had a significantly higher mutant load than MT-ND3 and MT-ND5 (P < 0.001, each). By contrast, MT-ND5 had a significantly lower mutant load than MT-ND3 (P = 0.007). Notably, the mutation loads varied significantly among patients carrying the MT-ATP6, MT-ND3, and MT-ND5 mutations.

CONCLUSIONS

Our study illustrated the heteroplasmic diversity and phenotypic expression threshold of mutated mitochondrial genes in mtDNA-associated Leigh syndrome. The results provide promising insights into the genotype-phenotype correlation in mtDNA-associated Leigh syndrome that are expected to guide the development of tailored treatments for Leigh syndrome.

Mito-SiPE is a sequence-independent and PCR-free mtDNA enrichment method for accurate ultra-deep mitochondrial sequencing

The analysis of somatic variation in the mitochondrial genome requires deep sequencing of mitochondrial DNA. This is ordinarily achieved by selective enrichment methods, such as PCR amplification or probe hybridization. These methods can introduce bias and are prone to contamination by nuclear-mitochondrial sequences (NUMTs), elements that can introduce artefacts into heteroplasmy analysis. We isolated intact mitochondria using differential centrifugation and alkaline lysis and subjected purified mitochondrial DNA to a sequence-independent and PCR-free method to obtain ultra-deep (>80,000X) sequencing coverage of the mitochondrial genome. This methodology avoids false-heteroplasmy calls that occur when long-range PCR amplification is used for mitochondrial DNA enrichment. Previously published methods employing mitochondrial DNA purification did not measure mitochondrial DNA enrichment or utilise high coverage short-read sequencing. Here, we describe a protocol that yields mitochondrial DNA and have quantified the increased level of mitochondrial DNA post-enrichment in 7 different mouse tissues. This method will enable researchers to identify changes in low frequency heteroplasmy without introducing PCR biases or NUMT contamination that are incorrectly identified as heteroplasmy when long-range PCR is used.

Long read mitochondrial genome sequencing using Cas9-guided adaptor ligation

The mitochondrial genome (mtDNA) is an important source of disease-causing genetic variability, but existing sequencing methods limit understanding, precluding phased measurement of mutations and clear detection of large sporadic deletions. We adapted a method for amplification-free sequence enrichment using Cas9 cleavage to obtain full length nanopore reads of mtDNA. We then utilized the long reads to phase mutations in a patient with an mtDNA-linked syndrome and demonstrated that this method can map age-induced mtDNA deletions. We believe this method will offer deeper insight into our understanding of mtDNA variation.

Genotype-phenotype analysis of MT-ATP6-associated Leigh syndrome

OBJECTIVES

Mitochondrial DNA (mtDNA)-associated Leigh syndrome (LS) is characterized by maternal inheritance, and the heteroplasmic mutant load of mtDNA pathogenic variants is known to affect clinical phenotypes. Among mtDNA pathogenic variants, variants of the MT-ATP6 gene account for most of reported cases. In this report, we aimed to describe the clinical and genetic findings of MT-ATP6-associated LS patients diagnosed at a single tertiary institution in Korea.

METHODS

Thirteen patients with genetically confirmed MT-ATP6-associated LS were selected. We reviewed each patient's clinical findings, including general characteristics, biochemical parameters, brain MR images, muscle biopsy results, and heteroplasmic mutant load over a long-term follow-up period.

RESULTS

MT-ATP6-associated LS was of predominantly early onset (age <2 years), although we identified 2 late-onset (>60 months) LS patients. The heteroplasmic mutant load estimated by next-generation sequencing was 96%-100% in all nucleotide change groups. Compared with other forms of MT-ATP6-associated LS, the m.8993T>G point mutation elicited a significantly higher rate of symptom onset before 2 years of age. Brain MRI showed bilateral basal ganglia involvement in all patients, followed by cerebral atrophy, brainstem and thalamus involvement, and cerebellar atrophy. After follow-up (median 7.2 years, range 1.4 to 11.5 years), LS with m.8993T>G point mutations had a slightly more severe clinical progression compared with other forms of MT-ATP6-associated LS.

CONCLUSIONS

MT-ATP6-associated LS patients presented with a broad spectrum of clinical diagnoses and had a very high heteroplasmic mutant load. This study provides valuable data on MT-ATP6-associated LS that will inform subsequent studies on LS.

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.

Association Between Epilepsy and Leigh Syndrome With MT-ND3 Mutation, Particularly the m.10191T>C Point Mutation

Background and Purpose: Recent advances in molecular genetic testing have led to a rapid increase in the understanding of the genetics of Leigh syndrome. Several studies have suggested that Leigh syndrome with MT-ND3 mutation is strongly associated with epilepsy. This study focused on the epilepsy-related characteristics of Leigh syndrome with MT-ND3 mutation identified in a single tertiary hospital in South Korea. Methods: We selected 31 patients with mitochondrial DNA (mtDNA) mutations who were genetically diagnosed with mtDNA-associated Leigh syndrome. Among them, seven patients with MT-ND3 mutations were detected. We reviewed various clinical findings such as laboratory findings, brain images, electroencephalography data, seizure types, seizure frequency, antiepileptic drug use history, and current seizure status. Results: The nucleotide changes in the seven patients with the Leigh syndrome with MT-ND3 mutation were divided into two groups: m.10191T>C and m.10158T>C. Six of the seven patients were found to have the m.10191T>C mutations. The median value of the mutant load was 82.5%, ranging from 57.9 to 93.6%. No particular tendency was observed for the first symptom or seizure onset or mutant load. The six patients with the m.10191T>C mutation were diagnosed with epilepsy. Three of these patients were diagnosed with Lennox-Gastaut syndrome (LGS). Conclusion: We reported a very strong association between epilepsy and MT-ND3 mutation in Leigh syndrome, particularly the m.10191T>C mutation. The possibility of an association between the epilepsy phenotype of the m.10191T>C mutation and LGS was noted.

Extension of Mitogenome Enrichment Based on Single Long-Range PCR: mtDNAs and Putative Mitochondrial-Derived Peptides of Five Rodent Hibernators

Enriching mitochondrial DNA (mtDNA) for sequencing entire mitochondrial genomes (mitogenomes) can be achieved by single long-range PCR. This avoids interference from the omnipresent nuclear mtDNA sequences (NUMTs). The approach is currently restricted to the use of samples collected from humans and ray-finned fishes. Here, we extended the use of single long-range PCR by introducing back-to-back oligonucleotides that target a sequence of extraordinary homology across vertebrates. The assay was applied to five hibernating rodents, namely alpine marmot, Arctic and European ground squirrels, and common and garden dormice, four of which have not been fully sequenced before. Analysis of the novel mitogenomes focussed on the prediction of mitochondrial-derived peptides (MDPs) providing another level of information encoded by mtDNA. The comparison of MOTS-c, SHLP4 and SHLP6 sequences across vertebrate species identified segments of high homology that argue for future experimentation. In addition, we evaluated four candidate polymorphisms replacing an amino acid in mitochondrially encoded subunits of the oxidative phosphorylation (OXPHOS) system that were reported in relation to cold-adaptation. No obvious pattern was found for the diverse sets of mammalian species that either apply daily or multiday torpor or otherwise cope with cold. In summary, our single long-range PCR assay applying a pair of back-to-back primers that target a consensus sequence motif of Vertebrata has potential to amplify (intact) mitochondrial rings present in templates from a taxonomically diverse range of vertebrates. It could be promising for studying novel mitogenomes, mitotypes of a population and mitochondrial heteroplasmy in a sensitive, straightforward and flexible manner.

Dystrophin Dp71 Subisoforms Localize to the Mitochondria of Human Cells

Duchenne muscular dystrophy (DMD) is a fatal muscle wasting disease caused by deficiency in dystrophin, a protein product encoded by the DMD gene. Mitochondrial dysfunction is now attracting much attention as a central player in DMD pathology. However, dystrophin has never been explored in human mitochondria. Here, we analyzed dystrophin in cDNAs and mitochondrial fractions of human cells. Mitochondrial fraction was obtained using a magnetic-associated cell sorting (MACS) technology. Dystrophin was analyzed by reverse transcription (RT)-PCR and western blotting using an antibody against the dystrophin C-terminal. In isolated mitochondrial fraction from HEK293 cells, dystrophin was revealed as a band corresponding to Dp71b and Dp71ab subisoforms. Additionally, in mitochondria from HeLa, SH-SY5Y, CCL-136 and HepG2 cells, signals for Dp71b and Dp71ab were revealed as well. Concomitantly, dystrophin mRNAs encoding Dp71b and Dp71ab were disclosed by RT-PCR in these cells. Primary cultured myocytes from three dystrophinopathy patients showed various levels of mitochondrial Dp71 expression. Coherently, levels of mRNA were different in all cells reflecting the protein content, which indicated predominant accumulation of Dp71. Dystrophin was demonstrated to be localized to human mitochondrial fraction, specifically as Dp71 subisoforms. Myocytes derived from dystrophinopathy patients manifested different levels of mitochondrial Dp71, with higher expression revealed in myocytes from Becker muscular dystrophy (BMD) patient-derived myocytes.

Mitochondrial DNA variation and cancer

Variation in the mitochondrial DNA (mtDNA) sequence is common in certain tumours. Two classes of cancer mtDNA variants can be identified: de novo mutations that act as 'inducers' of carcinogenesis and functional variants that act as 'adaptors', permitting cancer cells to thrive in different environments. These mtDNA variants have three origins: inherited variants, which run in families, somatic mutations arising within each cell or individual, and variants that are also associated with ancient mtDNA lineages (haplogroups) and are thought to permit adaptation to changing tissue or geographic environments. In addition to mtDNA sequence variation, mtDNA copy number and perhaps transfer of mtDNA sequences into the nucleus can contribute to certain cancers. Strong functional relevance of mtDNA variation has been demonstrated in oncocytoma and prostate cancer, while mtDNA variation has been reported in multiple other cancer types. Alterations in nuclear DNA-encoded mitochondrial genes have confirmed the importance of mitochondrial metabolism in cancer, affecting mitochondrial reactive oxygen species production, redox state and mitochondrial intermediates that act as substrates for chromatin-modifying enzymes. Hence, subtle changes in the mitochondrial genotype can have profound effects on the nucleus, as well as carcinogenesis and cancer progression.

Prediction of mitochondrial genome-wide variation through sequencing of mitochondrion-enriched extracts

Although mitochondrial DNA (mtDNA) haplotype variation is often applied for estimating population dynamics and phylogenetic relationships, economical and generalized methods for entire mtDNA genome enrichment prior to high-throughput sequencing are not readily available. This study demonstrates the utility of differential centrifugation to enrich for mitochondrion within cell extracts prior to DNA extraction, short-read sequencing, and assembly using exemplars from eight maternal lineages of the insect species, Ostrinia nubilalis. Compared to controls, enriched extracts showed a significant mean increase of 48.2- and 86.1-fold in mtDNA based on quantitative PCR, and proportion of subsequent short sequence reads that aligned to the O. nubilalis reference mitochondrial genome, respectively. Compared to the reference genome, our de novo assembled O. nubilalis mitochondrial genomes contained 82 intraspecific substitution and insertion/deletion mutations, and provided evidence for correction of mis-annotated 28 C-terminal residues within the NADH dehydrogenase subunit 4. Comparison to a more recent O. nubilalis mtDNA assembly from unenriched short-read data analogously showed 77 variant sites. Twenty-eight variant positions, and a triplet ATT codon (Ile) insertion within ATP synthase subunit 8, were unique within our assemblies. This study provides a generalizable pipeline for whole mitochondrial genome sequence acquisition adaptable to applications across a range of taxa.

STAMP: a multiplex sequencing method for simultaneous evaluation of mitochondrial DNA heteroplasmies and content

Human mitochondrial genome (mtDNA) variations, such as mtDNA heteroplasmies (the co-existence of mutated and wild-type mtDNA), have received increasing attention in recent years for their clinical relevance to numerous diseases. But large-scale population studies of mtDNA heteroplasmies have been lagging due to the lack of a labor- and cost-effective method. Here, we present a novel human mtDNA sequencing method called STAMP (sequencing by targeted amplification of multiplex probes) for measuring mtDNA heteroplasmies and content in a streamlined workflow. We show that STAMP has high-mapping rates to mtDNA, deep coverage of unique reads and high tolerance to sequencing and polymerase chain reaction errors when applied to human samples. STAMP also has high sensitivity and low false positive rates in identifying artificial mtDNA variants at fractions as low as 0.5% in genomic DNA samples. We further extend STAMP, by including nuclear DNA-targeting probes, to enable assessment of relative mtDNA content in the same assay. The high cost-effectiveness of STAMP, along with the flexibility of using it for measuring various aspects of mtDNA variations, will accelerate the research of mtDNA heteroplasmies and content in large population cohorts, and in the context of human diseases and aging.

Complexities and pitfalls in analyzing and interpreting mitochondrial DNA content in human cancer

Mutations in the human mitochondrial genome have been observed in all types of human cancer, indicating that mutations might contribute to tumorigenesis, metastasis, recurrence, or drug response. This possibility is appealing because of the known shift from oxidative metabolism to glycolysis, known as the Warburg effect, that occurs in malignancy. Mitochondrial DNA (mtDNA) mutations could either be maternally inherited and predispose to cancer (germ line mutations) or occur sporadically in the mtDNA of specific tissues (tissue- or tumor-specific somatic mutations) and contribute to the tumor initiation and progression process. High-throughput sequencing technologies now enable comprehensive detection of mtDNA variation in tissues and bodily fluids, with the potential to be used as an early detection tool that may impact the treatment of cancer. Here, we discuss insights into the roles of mtDNA mutations in carcinogenesis, highlighting the complexities involved in the analysis and interpretation of mitochondrial genomic content, technical challenges in studying their contribution to pathogenesis, and the value of mtDNA mutations in developing early detection, diagnosis, prognosis, and therapeutic strategies for cancer.

N6-Deoxyadenosine Methylation in Mammalian Mitochondrial DNA

N⁶-Methyldeoxyadenosine (6mA) has recently been shown to exist and play regulatory roles in eukaryotic genomic DNA (gDNA). However, the biological functions of 6mA in mammals have yet to be adequately explored, largely due to its low abundance in most mammalian genomes. Here, we report that mammalian mitochondrial DNA (mtDNA) is enriched for 6mA. The level of 6mA in HepG2 mtDNA is at least 1,300-fold higher than that in gDNA under normal growth conditions, corresponding to approximately four 6mA modifications on each mtDNA molecule. METTL4, a putative mammalian methyltransferase, can mediate mtDNA 6mA methylation, which contributes to attenuated mtDNA transcription and a reduced mtDNA copy number. Mechanistically, the presence of 6mA could repress DNA binding and bending by mitochondrial transcription factor (TFAM). Under hypoxia, the 6mA level in mtDNA could be further elevated, suggesting regulatory roles for 6mA in mitochondrial stress response. Our study reveals DNA 6mA as a regulatory mark in mammalian mtDNA.

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

Background

The 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.

Method

We 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).

Results

The 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.

Conclusion

We 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.

Sensitive detection of mitochondrial DNA variants for analysis of mitochondrial DNA-enriched extracts from frozen tumor tissue

Large variation exists in mitochondrial DNA (mtDNA) not only between but also within individuals. Also in human cancer, tumor-specific mtDNA variation exists. In this work, we describe the comparison of four methods to extract mtDNA as pure as possible from frozen tumor tissue. Also, three state-of-the-art methods for sensitive detection of mtDNA variants were evaluated. The main aim was to develop a procedure to detect low-frequent single-nucleotide mtDNA-specific variants in frozen tumor tissue. We show that of the methods evaluated, DNA extracted from cytosol fractions following exonuclease treatment results in highest mtDNA yield and purity from frozen tumor tissue (270-fold mtDNA enrichment). Next, we demonstrate the sensitivity of detection of low-frequent single-nucleotide mtDNA variants (≤1% allele frequency) in breast cancer cell lines MDA-MB-231 and MCF-7 by single-molecule real-time (SMRT) sequencing, UltraSEEK chemistry based mass spectrometry, and digital PCR. We also show de novo detection and allelic phasing of variants by SMRT sequencing. We conclude that our sensitive procedure to detect low-frequent single-nucleotide mtDNA variants from frozen tumor tissue is based on extraction of DNA from cytosol fractions followed by exonuclease treatment to obtain high mtDNA purity, and subsequent SMRT sequencing for (de novo) detection and allelic phasing of variants.

Detection and quantification of mitochondrial DNA deletions from next-generation sequence data

BACKGROUND

Chromosomal deletions represent an important class of human genetic variation. Various methods have been developed to mine "next-generation" sequencing (NGS) data to detect deletions and quantify their clonal abundances. These methods have focused almost exclusively on the nuclear genome, ignoring the mitochondrial chromosome (mtDNA). Detecting mtDNA deletions requires special care. First, the chromosome's relatively small size (16,569 bp) necessitates the ability to detect extremely focal events. Second, the chromosome can be present at thousands of copies in a single cell (in contrast to two copies of nuclear chromosomes), and mtDNA deletions may be present on only a very small percentage of chromosomes. Here we present a method, termed MitoDel, to detect mtDNA deletions from NGS data.

RESULTS

We validate the method on simulated and real data, and show that MitoDel can detect novel and previously-reported mtDNA deletions. We establish that MitoDel can find deletions such as the "common deletion" at heteroplasmy levels well below 1%.

CONCLUSIONS

MitoDel is a tool for detecting large mitochondrial deletions at low heteroplasmy levels. The tool can be downloaded at http://mendel.gene.cwru.edu/laframboiselab/ .

MitoRS, a method for high throughput, sensitive, and accurate detection of mitochondrial DNA heteroplasmy

Background

Mitochondrial dysfunction is linked to numerous pathological states, in particular related to metabolism, brain health and ageing. Nuclear encoded gene polymorphisms implicated in mitochondrial functions can be analyzed in the context of classical genome wide association studies. By contrast, mitochondrial DNA (mtDNA) variants are more challenging to identify and analyze for several reasons. First, contrary to the diploid nuclear genome, each cell carries several hundred copies of the circular mitochondrial genome. Mutations can therefore be present in only a subset of the mtDNA molecules, resulting in a heterogeneous pool of mtDNA, a situation referred to as heteroplasmy. Consequently, detection and quantification of variants requires extremely accurate tools, especially when this proportion is small. Additionally, the mitochondrial genome has pseudogenized into numerous copies within the nuclear genome over the course of evolution. These nuclear pseudogenes, named NUMTs, must be distinguished from genuine mtDNA sequences and excluded from the analysis.

Results

Here we describe a novel method, named MitoRS, in which the entire mitochondrial genome is amplified in a single reaction using rolling circle amplification. This approach is easier to setup and of higher throughput when compared to classical PCR amplification. Sequencing libraries are generated at high throughput exploiting a tagmentation-based method. Fine-tuned parameters are finally applied in the analysis to allow detection of variants even of low frequency heteroplasmy. The method was thoroughly benchmarked in a set of experiments designed to demonstrate its robustness, accuracy and sensitivity. The MitoRS method requires 5 ng total DNA as starting material. More than 96 samples can be processed in less than a day of laboratory work and sequenced in a single lane of an Illumina HiSeq flow cell. The lower limit for accurate quantification of single nucleotide variants has been measured at 1% frequency.

Conclusions

The MitoRS method enables the robust, accurate, and sensitive analysis of a large number of samples. Because it is cost effective and simple to setup, we anticipate this method will promote the analysis of mtDNA variants in large cohorts, and may help assessing the impact of mtDNA heteroplasmy on metabolic health, brain function, cancer progression, or ageing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3695-5) contains supplementary material, which is available to authorized users.