Next generation sequencing to characterize mitochondrial genomic DNA heteroplasmy

Mitochondrial DNA mutations in extremely preterm infants with bronchopulmonary dysplasia

Bronchopulmonary dysplasia (BPD) is a serious chronic lung disease affecting extremely preterm infants. While mitochondrial dysfunction has been investigated in various medical conditions, limited research has explored mitochondrial DNA (mtDNA) gene mutations, specifically in BPD. This study aimed to evaluate mitochondrial mtDNA gene mutations in extremely preterm infants with BPD. In this prospective observational study, we enrolled a cohort of extremely preterm infants diagnosed with BPD. Clinical data were collected to provide comprehensive patient profiles. Peripheral blood mononuclear cells were isolated from whole-blood samples obtained within a defined timeframe. Subsequently, mtDNA extraction and sequencing using next-generation sequencing technology were performed to identify mtDNA gene mutations. Among the cohort of ten extremely preterm infants with BPD, mtDNA sequencing revealed the presence of mutations in seven patients, resulting in a total of twenty-one point mutations. Notably, many of these mutations were identified in loci associated with critical components of the respiratory chain complexes, vital for proper mitochondrial function and cellular energy production. This pilot study provides evidence of mtDNA point mutations in a subset of extremely preterm infants with BPD. These findings suggest a potential association between mitochondrial dysfunction and the pathogenesis of BPD. Further extensive investigations are warranted to unravel the mechanisms underlying mtDNA mutations in BPD.

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

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

Use of Next-Generation Sequencing for Identifying Mitochondrial Disorders

Mitochondria are major contributors to ATP synthesis, generating more than 90% of the total cellular energy production through oxidative phosphorylation (OXPHOS): metabolite oxidation, such as the β-oxidation of fatty acids, and the Krebs’s cycle. OXPHOS inadequacy due to large genetic lesions in mitochondrial as well as nuclear genes and homo- or heteroplasmic point mutations in mitochondrially encoded genes is a characteristic of heterogeneous, maternally inherited genetic disorders known as mitochondrial disorders that affect multisystemic tissues and organs with high energy requirements, resulting in various signs and symptoms. Several traditional diagnostic approaches, including magnetic resonance imaging of the brain, cardiac testing, biochemical screening, variable heteroplasmy genetic testing, identifying clinical features, and skeletal muscle biopsies, are associated with increased risks, high costs, a high degree of false-positive or false-negative results, or a lack of precision, which limits their diagnostic abilities for mitochondrial disorders. Variable heteroplasmy levels, mtDNA depletion, and the identification of pathogenic variants can be detected through genetic sequencing, including the gold standard Sanger sequencing. However, sequencing can be time consuming, and Sanger sequencing can result in the missed recognition of larger structural variations such as CNVs or copy-number variations. Although each sequencing method has its own limitations, genetic sequencing can be an alternative to traditional diagnostic methods. The ever-growing roster of possible mutations has led to the development of next-generation sequencing (NGS). The enhancement of NGS methods can offer a precise diagnosis of the mitochondrial disorder within a short period at a reasonable expense for both research and clinical applications.

Comparison of whole genome sequencing and targeted sequencing for mitochondrial DNA

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

Identification of unique and shared mitochondrial DNA mutations in neurodegeneration and cancer by single-cell mitochondrial DNA structural variation sequencing (MitoSV-seq)

Background

Point mutations and structural variations (SVs) in mitochondrial DNA (mtDNA) contribute to many neurodegenerative diseases. Technical limitations and heteroplasmy, however, have impeded their identification, preventing these changes from being examined in neurons in healthy and disease states.

Methods

We have developed a high-resolution technique—Mitochondrial DNA Structural Variation Sequencing (MitoSV-seq)—that identifies all types of mtDNA SVs and single-nucleotide variations (SNVs) in single neurons and novel variations that have been undetectable with conventional techniques.

Findings

Using MitoSV-seq, we discovered SVs/SNVs in dopaminergic neurons in the Ifnar1^−/− murine model of Parkinson disease. Further, MitoSV-seq was found to have broad applicability, delivering high-quality, full-length mtDNA sequences in a species-independent manner from human PBMCs, haematological cancers, and tumour cell lines, regardless of heteroplasmy. We characterised several common SVs in haematological cancers (AML and MDS) that were linked to the same mtDNA region, MT-ND5, using only 10 cells, indicating the power of MitoSV-seq in determining single-cancer-cell ontologies. Notably, the MT-ND5 hotspot, shared between all examined cancers and Ifnar1^−/− dopaminergic neurons, suggests that its mutations have clinical value as disease biomarkers.

Interpretation

MitoSV-seq identifies disease-relevant mtDNA mutations in single cells with high resolution, rendering it a potential drug screening platform in neurodegenerative diseases and cancers.

Funding

The Lundbeck Foundation, Danish Council for Independent Research-Medicine, and European Union Horizon 2020 Research and Innovation Programme.

Investigating Human Mitochondrial Genomes in Single Cells

Mitochondria host multiple copies of their own small circular genome that has been extensively studied to trace the evolution of the modern eukaryotic cell and discover important mutations linked to inherited diseases. Whole genome and exome sequencing have enabled the study of mtDNA in a large number of samples and experimental conditions at single nucleotide resolution, allowing the deciphering of the relationship between inherited mutations and phenotypes and the identification of acquired mtDNA mutations in classical mitochondrial diseases as well as in chronic disorders, ageing and cancer. By applying an ad hoc computational pipeline based on our MToolBox software, we reconstructed mtDNA genomes in single cells using whole genome and exome sequencing data obtained by different amplification methodologies (eWGA, DOP-PCR, MALBAC, MDA) as well as data from single cell Assay for Transposase Accessible Chromatin with high-throughput sequencing (scATAC-seq) in which mtDNA sequences are expected as a byproduct of the technology. We show that assembled mtDNAs, with the exception of those reconstructed by MALBAC and DOP-PCR methods, are quite uniform and suitable for genomic investigations, enabling the study of various biological processes related to cellular heterogeneity such as tumor evolution, neural somatic mosaicism and embryonic development.

The special considerations of gene therapy for mitochondrial diseases

The recent success of gene therapy across multiple clinical trials has inspired a great deal of hope regarding the treatment of previously intractable genetic diseases. This optimism has been extended to the prospect of gene therapy for mitochondrial disorders, which are not only particularly severe but also difficult to treat. However, this hope must be tempered by the reality of the mitochondrial organelle, which possesses specific biological properties that complicate genetic manipulation. In this perspective, we will discuss some of these complicating factors, including the unique pathways used to express and import mitochondrial proteins. We will also present some ways in which these challenges can be overcome by genetic manipulation strategies tailored specifically for mitochondrial diseases.

Detection of mitochondrial DNA (mtDNA) mutations

The maternally inherited mitochondrial DNA (mtDNA) is a circular 16,569bp double stranded DNA that encodes 37 genes, 24 of which (2 rRNAs and 22 tRNAs) are necessary for transcription and translation of 13 polypeptides that are all subunits of respiratory chain. Pathogenic mutations in mtDNA cause respiratory chain dysfunction, and are the underlying defect in an ever-increasing number of mtDNA-related encephalomyopathies with distinct phenotypes. In this chapter, we present an overview of mtDNA mutations and describe the molecular techniques currently employed in our laboratory to detect two types of mtDNA mutations: single-large-scale rearrangements and point mutations.

Leishmania Mitochondrial Genomes: Maxicircle Structure and Heterogeneity of Minicircles

The mitochondrial DNA (mtDNA), which is present in almost all eukaryotic organisms, is a useful marker for phylogenetic studies due to its relative high conservation and its inheritance manner. In Leishmania and other trypanosomatids, the mtDNA (also referred to as kinetoplast DNA or kDNA) is composed of thousands of minicircles and a few maxicircles, catenated together into a complex network. Maxicircles are functionally similar to other eukaryotic mtDNAs, whereas minicircles are involved in RNA editing of some maxicircle-encoded transcripts. Next-generation sequencing (NGS) is increasingly used for assembling nuclear genomes and, currently, a large number of genomic sequences are available. However, most of the time, the mitochondrial genome is ignored in the genome assembly processes. The aim of this study was to develop a pipeline to assemble Leishmania minicircles and maxicircle DNA molecules, exploiting the raw data generated in the NGS projects. As a result, the maxicircle molecules and the plethora of minicircle classes for Leishmania major, Leishmania infantum and Leishmania braziliensis have been characterized. We have observed that whereas the heterogeneity of minicircle sequences existing in a single cell hampers their use for Leishmania typing and classification, maxicircles emerge as an extremely robust genetic marker for taxonomic studies within the clade of kinetoplastids.

Genotypes of chronic progressive external ophthalmoplegia in a large adult-onset cohort

Chronic progressive external ophthalmoplegia (CPEO) is a common presentation of mitochondrial disease. We performed a retrospective evaluation of the molecular genetic testing and genotype-phenotype correlations in a large cohort of adult-onset CPEO patients (N = 111). One hundred percent of patients tested had at least one mitochondrial DNA (mtDNA) deletion. Genetic testing of nuclear genes encoding mitochondrial proteins identified pathogenic/likely pathogenic variants likely to be associated with CPEO in 7.6% of patients. As expected, the nuclear gene most associated with DNA variation was POLG. A single likely pathogenic mitochondrial DNA variant (m.12278T>C) was identified in two unrelated patients. No significant differences were noted in the clinical phenotypes of patients with pathogenic or likely pathogenic nuclear variants in comparison to those with negative nuclear gene testing. Analysis of deletion size and heteroplasmy in muscle-derived mtDNA showed significant correlations with age of symptom onset but not disease severity (number of canonical CPEO features). Results suggest that smaller mtDNA deletions (p = 0.0127, r² = 0.1201) and higher heteroplasmy of single mtDNA deletions (p = 0.0112, r² = 0.2483) are associated with an earlier age of onset in CPEO patients.

Next-generation sequencing identifies novel mitochondrial variants in pituitary adenomas

PURPOSE

Disrupted mitochondrial functions and genetic variants of mitochondrial DNA (mtDNA) have been observed in different human neoplasms. Next-generation sequencing (NGS) can be used to detect even low heteroplasmy-level mtDNA variants. We aimed to investigate the mitochondrial genome in pituitary adenomas by NGS.

METHODS

We analysed 11 growth hormone producing and 33 non-functioning [22 gonadotroph and 11 hormone immunonegative] pituitary adenomas using VariantPro™ Mitochondrion Panel on Illumina MiSeq instrument. Revised Cambridge Reference Sequence (rCRS) of the mtDNA was used as reference. Heteroplasmy was determined using a 3% cutoff.

RESULTS

496 variants were identified in pituitary adenomas with overall low level of heteroplasmy (7.22%). On average, 35 variants were detected per sample. Samples harbouring the highest number of variants had the highest Ki-67 indices independently of histological subtypes. We identified eight variants (A11251G, T4216C, T16126C, C15452A, T14798C, A188G, G185A, and T16093C) with different prevalences among different histological groups. T16189C was found in 40% of non-recurrent adenomas, while it was not present in the recurrent ones. T14798C and T4216C were confirmed by Sanger sequencing in all 44 samples. 100% concordance was found between NGS and Sanger method.

CONCLUSIONS

NGS is a reliable method for investigating mitochondrial genome and heteroplasmy in pituitary adenomas. Out of the 496 detected variants, 414 have not been previously reported in pituitary adenoma. The high number of mtDNA variants may contribute to adenoma genesis, and some variants (i.e., T16189C) might associate with benign behaviour.

A novel variant m.641A>T in the mitochondrial MT-TF gene is associated with epileptic encephalopathy in adolescent

We present a 14-year-old girl with loss of motor functions, tetraplegia, epilepsy and nystagmus, caused by a novel heteroplasmic m.641A>T transition in an evolutionary conserved region of mitochondrial genome, affecting the aminoacyl stem of mitochondrial tRNA-Phe. In silico prediction, respirometry, Western blot and enzymatic analyses in skin fibroblasts support the pathogenicity of the m.641A>T substitution. This is the 18th MT-TF point mutation associated with a mitochondrial disorder. The onset and the severity of the disease, however, is unique in this case and broadens the clinical picture related to mutations of mitochondrial tRNA-Phe.

Cause or casualty: The role of mitochondrial DNA in aging and age-associated disease

The mitochondrial genome (mtDNA) represents a tiny fraction of the whole genome, comprising just 16.6 kilobases encoding 37 genes involved in oxidative phosphorylation and the mitochondrial translation machinery. Despite its small size, much interest has developed in recent years regarding the role of mtDNA as a determinant of both aging and age-associated diseases. A number of studies have presented compelling evidence for key roles of mtDNA in age-related pathology, although many are correlative rather than demonstrating cause. In this review we will evaluate the evidence supporting and opposing a role for mtDNA in age-associated functional declines and diseases. We provide an overview of mtDNA biology, damage and repair as well as the influence of mitochondrial haplogroups, epigenetics and maternal inheritance in aging and longevity.

Bioinformatics Tools and Databases to Assess the Pathogenicity of Mitochondrial DNA Variants in the Field of Next Generation Sequencing

The development of next generation sequencing (NGS) has greatly enhanced the diagnosis of mitochondrial disorders, with a systematic analysis of the whole mitochondrial DNA (mtDNA) sequence and better detection sensitivity. However, the exponential growth of sequencing data renders complex the interpretation of the identified variants, thereby posing new challenges for the molecular diagnosis of mitochondrial diseases. Indeed, mtDNA sequencing by NGS requires specific bioinformatics tools and the adaptation of those developed for nuclear DNA, for the detection and quantification of mtDNA variants from sequence alignment to the calling steps, in order to manage the specific features of the mitochondrial genome including heteroplasmy, i.e., coexistence of mutant and wildtype mtDNA copies. The prioritization of mtDNA variants remains difficult, relying on a limited number of specific resources: population and clinical databases, and in silico tools providing a prediction of the variant pathogenicity. An evaluation of the most prominent bioinformatics tools showed that their ability to predict the pathogenicity was highly variable indicating that special efforts should be directed at developing new bioinformatics tools dedicated to the mitochondrial genome. In addition, massive parallel sequencing raised several issues related to the interpretation of very low mtDNA mutational loads, discovery of variants of unknown significance, and mutations unrelated to patient phenotype or the co-occurrence of mtDNA variants. This review provides an overview of the current strategies and bioinformatics tools for accurate annotation, prioritization and reporting of mtDNA variations from NGS data, in order to carry out accurate genetic counseling in individuals with primary mitochondrial diseases.

Biparental Inheritance of Mitochondrial DNA in Humans

Although there has been considerable debate about whether paternal mitochondrial DNA (mtDNA) transmission may coexist with maternal transmission of mtDNA, it is generally believed that mitochondria and mtDNA are exclusively maternally inherited in humans. Here, we identified three unrelated multigeneration families with a high level of mtDNA heteroplasmy (ranging from 24 to 76%) in a total of 17 individuals. Heteroplasmy of mtDNA was independently examined by high-depth whole mtDNA sequencing analysis in our research laboratory and in two Clinical Laboratory Improvement Amendments and College of American Pathologists-accredited laboratories using multiple approaches. A comprehensive exploration of mtDNA segregation in these families shows biparental mtDNA transmission with an autosomal dominantlike inheritance mode. Our results suggest that, although the central dogma of maternal inheritance of mtDNA remains valid, there are some exceptional cases where paternal mtDNA could be passed to the offspring. Elucidating the molecular mechanism for this unusual mode of inheritance will provide new insights into how mtDNA is passed on from parent to offspring and may even lead to the development of new avenues for the therapeutic treatment for pathogenic mtDNA transmission.

Whole Exome Sequencing Is the Preferred Strategy to Identify the Genetic Defect in Patients With a Probable or Possible Mitochondrial Cause

Mitochondrial disorders, characterized by clinical symptoms and/or OXPHOS deficiencies, are caused by pathogenic variants in mitochondrial genes. However, pathogenic variants in some of these genes can lead to clinical manifestations which overlap with other neuromuscular diseases, which can be caused by pathogenic variants in non-mitochondrial genes as well. Mitochondrial pathogenic variants can be found in the mitochondrial DNA (mtDNA) or in any of the 1,500 nuclear genes with a mitochondrial function. We have performed a two-step next-generation sequencing approach in a cohort of 117 patients, mostly children, in whom a mitochondrial disease-cause could likely or possibly explain the phenotype. A total of 86 patients had a mitochondrial disorder, according to established clinical and biochemical criteria. The other 31 patients had neuromuscular symptoms, where in a minority a mitochondrial genetic cause is present, but a non-mitochondrial genetic cause is more likely. All patients were screened for pathogenic variants in the mtDNA and, if excluded, analyzed by whole exome sequencing (WES). Variants were filtered for being pathogenic and compatible with an autosomal or X-linked recessive mode of inheritance in families with multiple affected siblings and/or consanguineous parents. Non-consanguineous families with a single patient were additionally screened for autosomal and X-linked dominant mutations in a predefined gene-set. We identified causative pathogenic variants in the mtDNA in 20% of the patient-cohort, and in nuclear genes in 49%, implying an overall yield of 68%. We identified pathogenic variants in mitochondrial and non-mitochondrial genes in both groups with, obviously, a higher number of mitochondrial genes affected in mitochondrial disease patients. Furthermore, we show that 31% of the disease-causing genes in the mitochondrial patient group were not included in the MitoCarta database, and therefore would have been missed with MitoCarta based gene-panels. We conclude that WES is preferable to panel-based approaches for both groups of patients, as the mitochondrial gene-list is not complete and mitochondrial symptoms can be secondary. Also, clinically and genetically heterogeneous disorders would require sequential use of multiple different gene panels. We conclude that WES is a comprehensive and unbiased approach to establish a genetic diagnosis in these patients, able to resolve multi-genic disease-causes.

Reducing the information gap on Loricarioidei (Siluriformes) mitochondrial genomics

Background

The genetic diversity of Neotropical fish fauna is underrepresented in public databases. This distortion is evident for the order Siluriformes, in which the suborders Siluroidei and Loricarioidei share equivalent proportion of species, although far less is known about the genetics of the latter clade, endemic to the Neotropical Region. Recently, this information gap was evident in a study about the structural diversity of fish mitochondrial genomes, and hampered a precise chronological resolution of Siluriformes. It has also prevented molecular ecology investigations about these catfishes, their interactions with the environment, responses to anthropogenic changes and potential uses.

Results

Using high-throughput sequencing, we provide the nearly complete mitochondrial genomes for 26 Loricariidae and one Callichthyidae species. Structural features were highly conserved. A notable exception was identified in the monophyletic clade comprising species of the Hemiancistrus, Hypostomini and Peckoltia-clades, a ~60 nucleotide-long deletion encompassing the seven nucleotides at the 3′ end of the Conserved Sequence Block (CSB) D of the control region. The expression of mitochondrial genes followed the usual punctuation pattern. Heteroplasmic sites were identified in most species. The retrieved phylogeny strongly corroborates the currently accepted tree, although bringing to debate the relationship between Schizolecis guntheri and Pareiorhaphis garbei, and highlighting the low genetic variability within the Peckoltia-clade, an eco-morphologically diverse and taxonomically problematic group.

Conclusions

Herein we have launched the use of high-throughput mitochondrial genomics in the studies of the Loricarioidei species. The new genomic resources reduce the information gap on the molecular diversity of Neotropical fish fauna, impacting the capacity to investigate a variety of aspects of the molecular ecology and evolution of these fishes. Additionally, the species showing the partial CSB-D are candidate models to study the replication and transcription of vertebrate mitochondrial genome.

Electronic supplementary material

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

Current strategies towards therapeutic manipulation of mtDNA heteroplasmy

Mitochondrial disease is a multifactorial disorder involving both nuclear and mitochondrial genomes. Over the past 20 years, great progress was achieved in the field of gene editing which raised the possibility of partial or complete elimination of mutant mtDNA that causes disease phenotypes. Each cell contains thousands of copies of mtDNA which can be either wild-type (WT) or mutant, a condition called heteroplasmy. As there are multiple copies of mtDNA inside a cell, the percentage of mutant mtDNA can vary and a directional shift in the heteroplasmy ratio towards an increase of WT mtDNA copies would have therapeutic value. Gene editing tools have been adapted to translocate to mitochondria and were able to change heteroplasmy in a predictable manner. These include mitochondrial targeted restriction endonucleases, Zinc-finger nucleases, and TAL-effector nucleases. These procedures could also be adapted to reduce the levels of mutant mtDNA in embryos, offering an option to the controversial mitochondrial replacement techniques during in vitro fertilization. The current strategies to induce heteroplasmy shift of mtDNA and its implications will be comprehensively discussed.

The ability of human nuclear DNA to cause false positive low-abundance heteroplasmy calls varies across the mitochondrial genome

Background

Low-abundance mutations in mitochondrial populations (mutations with minor allele frequency ≤ 1%), are associated with cancer, aging, and neurodegenerative disorders. While recent progress in high-throughput sequencing technology has significantly improved the heteroplasmy identification process, the ability of this technology to detect low-abundance mutations can be affected by the presence of similar sequences originating from nuclear DNA (nDNA). To determine to what extent nDNA can cause false positive low-abundance heteroplasmy calls, we have identified mitochondrial locations of all subsequences that are common or similar (one mismatch allowed) between nDNA and mitochondrial DNA (mtDNA).

Results

Performed analysis revealed up to a 25-fold variation in the lengths of longest common and longest similar (one mismatch allowed) subsequences across the mitochondrial genome. The size of the longest subsequences shared between nDNA and mtDNA in several regions of the mitochondrial genome were found to be as low as 11 bases, which not only allows using these regions to design new, very specific PCR primers, but also supports the hypothesis of the non-random introduction of mtDNA into the human nuclear DNA.

Conclusion

Analysis of the mitochondrial locations of the subsequences shared between nDNA and mtDNA suggested that even very short (36 bases) single-end sequencing reads can be used to identify low-abundance variation in 20.4% of the mitochondrial genome. For longer (76 and 150 bases) reads, the proportion of the mitochondrial genome where nDNA presence will not interfere found to be 44.5 and 67.9%, when low-abundance mutations at 100% of locations can be identified using 417 bases long single reads. This observation suggests that the analysis of low-abundance variations in mitochondria population can be extended to a variety of large data collections such as NCBI Sequence Read Archive, European Nucleotide Archive, The Cancer Genome Atlas, and International Cancer Genome Consortium.

Electronic supplementary material

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

mit-o-matic: a comprehensive computational pipeline for clinical evaluation of mitochondrial variations from next-generation sequencing datasets

The human mitochondrial genome has been reported to have a very high mutation rate as compared with the nuclear genome. A large number of mitochondrial mutations show significant phenotypic association and are involved in a broad spectrum of diseases. In recent years, there has been a remarkable progress in the understanding of mitochondrial genetics. The availability of next-generation sequencing (NGS) technologies have not only reduced sequencing cost by orders of magnitude but has also provided us good quality mitochondrial genome sequences with high coverage, thereby enabling decoding of a number of human mitochondrial diseases. In this study, we report a computational and experimental pipeline to decipher the human mitochondrial DNA variations and examine them for their clinical correlation. As a proof of principle, we also present a clinical study of a patient with Leigh disease and confirmed maternal inheritance of the causative allele. The pipeline is made available as a user-friendly online tool to annotate variants and find haplogroup, disease association, and heteroplasmic sites. The "mit-o-matic" computational pipeline represents a comprehensive cloud-based tool for clinical evaluation of mitochondrial genomic variations from NGS datasets. The tool is freely available at http://genome.igib.res.in/mitomatic/.

Mitochondrial mutations in subjects with psychiatric disorders

A considerable body of evidence supports the role of mitochondrial dysfunction in psychiatric disorders and mitochondrial DNA (mtDNA) mutations are known to alter brain energy metabolism, neurotransmission, and cause neurodegenerative disorders. Genetic studies focusing on common nuclear genome variants associated with these disorders have produced genome wide significant results but those studies have not directly studied mtDNA variants. The purpose of this study is to investigate, using next generation sequencing, the involvement of mtDNA variation in bipolar disorder, schizophrenia, major depressive disorder, and methamphetamine use. MtDNA extracted from multiple brain regions and blood were sequenced (121 mtDNA samples with an average of 8,800x coverage) and compared to an electronic database containing 26,850 mtDNA genomes. We confirmed novel and rare variants, and confirmed next generation sequencing error hotspots by traditional sequencing and genotyping methods. We observed a significant increase of non-synonymous mutations found in individuals with schizophrenia. Novel and rare non-synonymous mutations were found in psychiatric cases in mtDNA genes: ND6, ATP6, CYTB, and ND2. We also observed mtDNA heteroplasmy in brain at a locus previously associated with schizophrenia (T16519C). Large differences in heteroplasmy levels across brain regions within subjects suggest that somatic mutations accumulate differentially in brain regions. Finally, multiplasmy, a heteroplasmic measure of repeat length, was observed in brain from selective cases at a higher frequency than controls. These results offer support for increased rates of mtDNA substitutions in schizophrenia shown in our prior results. The variable levels of heteroplasmic/multiplasmic somatic mutations that occur in brain may be indicators of genetic instability in mtDNA.

Accurate mitochondrial DNA sequencing using off-target reads provides a single test to identify pathogenic point mutations

PURPOSE

Mitochondrial disorders are a common cause of inherited metabolic disease and can be due to mutations affecting mitochondrial DNA or nuclear DNA. The current diagnostic approach involves the targeted resequencing of mitochondrial DNA and candidate nuclear genes, usually proceeds step by step, and is time consuming and costly. Recent evidence suggests that variations in mitochondrial DNA sequence can be obtained from whole-exome sequence data, raising the possibility of a comprehensive single diagnostic test to detect pathogenic point mutations.

METHODS

We compared the mitochondrial DNA sequence derived from off-target exome reads with conventional mitochondrial DNA Sanger sequencing in 46 subjects.

RESULTS

Mitochondrial DNA sequences can be reliably obtained using three different whole-exome sequence capture kits. Coverage correlates with the relative amount of mitochondrial DNA in the original genomic DNA sample, heteroplasmy levels can be determined using variant and total read depths, and-providing there is a minimum read depth of 20-fold-rare sequencing errors occur at a rate similar to that observed with conventional Sanger sequencing.

CONCLUSION

This offers the prospect of using whole-exome sequence in a diagnostic setting to screen not only all protein coding nuclear genes but also all mitochondrial DNA genes for pathogenic mutations. Off-target mitochondrial DNA reads can also be used to assess quality control and maternal ancestry, inform on ethnic origin, and allow genetic disease association studies not previously anticipated with existing whole-exome data sets.

A systematic approach to assessing the clinical significance of genetic variants

Molecular genetic testing informs diagnosis, prognosis, and risk assessment for patients and their family members. Recent advances in low-cost, high-throughput DNA sequencing and computing technologies have enabled the rapid expansion of genetic test content, resulting in dramatically increased numbers of DNA variants identified per test. To address this challenge, our laboratory has developed a systematic approach to thorough and efficient assessments of variants for pathogenicity determination. We first search for existing data in publications and databases including internal, collaborative and public resources. We then perform full evidence-based assessments through statistical analyses of observations in the general population and disease cohorts, evaluation of experimental data from in vivo or in vitro studies, and computational predictions of potential impacts of each variant. Finally, we weigh all evidence to reach an overall conclusion on the potential for each variant to be disease causing. In this report, we highlight the principles of variant assessment, address the caveats and pitfalls, and provide examples to illustrate the process. By sharing our experience and providing a framework for variant assessment, including access to a freely available customizable tool, we hope to help move towards standardized and consistent approaches to variant assessment.

Molecular diagnosis of autosomal dominant polycystic kidney disease using next-generation sequencing

Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in PKD1 and PKD2. However, genetic analysis is complicated by six PKD1 pseudogenes, large gene sizes, and allelic heterogeneity. We developed a new clinical assay for PKD gene analysis using paired-end next-generation sequencing (NGS) by multiplexing individually bar-coded long-range PCR libraries and analyzing them in one Illumina MiSeq flow cell. The data analysis pipeline has been optimized and automated with Unix shell scripts to accommodate variant calls. This approach was validated using a cohort of 25 patients with ADPKD previously analyzed by Sanger sequencing. A total of 250 genetic variants were identified by NGS, spanning the entire exonic and adjacent intronic regions of PKD1 and PKD2, including all 16 pathogenic mutations. In addition, we identified three novel mutations in a mutation-negative cohort of 24 patients with ADPKD previously analyzed by Sanger sequencing. This NGS method achieved sensitivity of 99.2% (95% CI, 96.8%-99.9%) and specificity of 99.9% (95% CI, 99.7%-100.0%), with cost and turnaround time reduced by as much as 70%. Prospective NGS analysis of 25 patients with ADPKD demonstrated a detection rate comparable with Sanger standards. In conclusion, the NGS method was superior to Sanger sequencing for detecting PKD gene mutations, achieving high sensitivity and improved gene coverage. These characteristics suggest that NGS would be an appropriate new standard for clinical genetic testing of ADPKD.

Assessment of clinical analytical sensitivity and specificity of next-generation sequencing for detection of simple and complex mutations

Background

Detecting mutations in disease genes by full gene sequence analysis is common in clinical diagnostic laboratories. Sanger dideoxy terminator sequencing allows for rapid development and implementation of sequencing assays in the clinical laboratory, but it has limited throughput, and due to cost constraints, only allows analysis of one or at most a few genes in a patient. Next-generation sequencing (NGS), on the other hand, has evolved rapidly, although to date it has mainly been used for large-scale genome sequencing projects and is beginning to be used in the clinical diagnostic testing. One advantage of NGS is that many genes can be analyzed easily at the same time, allowing for mutation detection when there are many possible causative genes for a specific phenotype. In addition, regions of a gene typically not tested for mutations, like deep intronic and promoter mutations, can also be detected.

Results

Here we use 20 previously characterized Sanger-sequenced positive controls in disease-causing genes to demonstrate the utility of NGS in a clinical setting using standard PCR based amplification to assess the analytical sensitivity and specificity of the technology for detecting all previously characterized changes (mutations and benign SNPs). The positive controls chosen for validation range from simple substitution mutations to complex deletion and insertion mutations occurring in autosomal dominant and recessive disorders. The NGS data was 100% concordant with the Sanger sequencing data identifying all 119 previously identified changes in the 20 samples.

Conclusions

We have demonstrated that NGS technology is ready to be deployed in clinical laboratories. However, NGS and associated technologies are evolving, and clinical laboratories will need to invest significantly in staff and infrastructure to build the necessary foundation for success.