Analysis of the whole mitochondrial genome: translation of the Ion Torrent Personal Genome Machine system to the diagnostic bench?

High frequency of mitochondrial DNA rearrangements in the peripheral blood of adults with intellectual disability

BACKGROUND

Mitochondrial DNA (mtDNA) rearrangements are recognised factors in mitochondrial disorders and ageing, but their involvement in neurodevelopmental disorders, particularly intellectual disability (ID) and autism spectrum disorder (ASD), remains poorly understood. Previous studies have reported mitochondrial dysfunction in individuals with both ID and ASD. The aim of this study was to investigate the prevalence of large-scale mtDNA rearrangements in ID and ID with comorbid ASD (ID-ASD).

METHOD

We used mtDNA-targeted next-generation sequencing and the MitoSAlt high-throughput computational pipeline in peripheral blood samples from 76 patients with ID (mean age 52.5 years, 37% female), 59 patients with ID-ASD (mean age 41.3 years, 46% female) and 32 healthy controls (mean age 42.4 years, 47% female) from Catalonia.

RESULTS

The study revealed a high frequency of mtDNA rearrangements in patients with ID, with 10/76 (13.2%) affected individuals. However, the prevalence was significantly lower in patients with ID-ASD 1/59 (1.7%) and in HC 1/32 (3.1%). Among the mtDNA rearrangements, six were identified as deletions (median size 6937 bp and median heteroplasmy level 2.3%) and six as duplications (median size 10 455 bp and median heteroplasmy level 1.9%). One of the duplications, MT-ATP6 m.8765-8793dup (29 bp), was present in four individuals with ID with a median heteroplasmy level of 3.9%.

CONCLUSIONS

Our results show that mtDNA rearrangements are frequent in patients with ID, but not in those with ID-ASD, when compared to HC. Additionally, MitoSAlt has demonstrated high sensitivity and accuracy in detecting mtDNA rearrangements, even at very low heteroplasmy levels in blood samples. While the high frequency of mtDNA rearrangements in ID is noteworthy, the role of these rearrangements is currently unclear and needs to be confirmed with further data, particularly in post-mitotic tissues and through age-matched control studies.

Mitochondrial dysfunction and oxidative stress in Alzheimer's disease, and Parkinson's disease, Huntington's disease and Amyotrophic Lateral Sclerosis -An updated review

Misfolded proteins in the central nervous system can induce oxidative damage, which can contribute to neurodegenerative diseases in the mitochondria. Neurodegenerative patients face early mitochondrial dysfunction, impacting energy utilization. Amyloid-ß and tau problems both have an effect on mitochondria, which leads to mitochondrial malfunction and, ultimately, the onset of Alzheimer's disease. Cellular oxygen interaction yields reactive oxygen species within mitochondria, initiating oxidative damage to mitochondrial constituents. Parkinson's disease, linked to oxidative stress, α-synuclein aggregation, and inflammation, results from reduced brain mitochondria activity. Mitochondrial dynamics profoundly influence cellular apoptosis via distinct causative mechanisms. The condition known as Huntington's disease is characterized by an expansion of polyglutamine, primarily impactingthe cerebral cortex and striatum. Research has identified mitochondrial failure as an early pathogenic mechanism contributing to HD's selective neurodegeneration. The mitochondria are organelles that exhibit dynamism by undergoing fragmentation and fusion processes to attain optimal bioenergetic efficiency. They can also be transported along microtubules and regulateintracellular calcium homeostasis through their interaction with the endoplasmic reticulum. Additionally, the mitochondria produce free radicals. The functions of eukaryotic cells, particularly in neurons, have significantly deviated from the traditionally assigned role of cellular energy production. Most of them areimpaired in HD, which may lead to neuronal dysfunction before symptoms manifest. This article summarises the most important changes in mitochondrial dynamics that come from neurodegenerative diseases including Alzheimer's, Parkinson's, Huntington's and Amyotrophic Lateral Sclerosis. Finally, we discussed about novel techniques that can potentially treat mitochondrial malfunction and oxidative stress in four most dominating neuro disorders.

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.

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.

SARS-CoV-2 genetic variations associated with COVID-19 pathogenicity

In this study, we performed genome-wide association analyses on SARS-CoV-2 genomes to identify genetic mutations associated with pre-symptomatic/asymptomatic COVID-19 cases. Various potential covariates and confounding factors of COVID-19 severity, including patient age, gender and country, as well as virus phylogenetic relatedness were adjusted for. In total, 3021 full-length genomes of SARS-CoV-2 generated from original clinical samples and whose patient status could be determined conclusively as either ‘pre-symptomatic/asymptomatic’ or ‘symptomatic’ were retrieved from the GISAID database. We found that the mutation 11 083G>T, located in the coding region of non-structural protein 6, is significantly associated with asymptomatic COVID-19. Patient age is positively correlated with symptomatic infection, while gender is not significantly correlated with the development of the disease. We also found that the effects of the mutation, patient age and gender do not vary significantly among countries, although each country appears to have varying baseline chances of COVID-19 symptom development.

Identification of full-length circular nucleic acids using long-read sequencing technologies

Unlike the traditional perception in genomic DNA or linear RNA, circular nucleic acids are a class of functional biomolecules with a circular configuration and are often observed in nature. These circular molecules encompass the full spectrum of size and play an important role in organisms, making circular nucleic acids research worthy. Due to the low abundance of most types of circular nucleic acids and the disadvantages of short-read sequencing platforms, accurate and full-length circular nucleic acid sequencing and identification is difficult. In this review, we have provided insights into full-length circular nucleic acid detection methods using long-read sequencing technologies, with a focus on the experimental and bioinformatics strategies to obtain accurate sequences.

Applying genomic and transcriptomic advances to mitochondrial medicine

Next-generation sequencing (NGS) has increased our understanding of the molecular basis of many primary mitochondrial diseases (PMDs). Despite this progress, many patients with suspected PMD remain without a genetic diagnosis, which restricts their access to in-depth genetic counselling, reproductive options and clinical trials, in addition to hampering efforts to understand the underlying disease mechanisms. Although they represent a considerable improvement over their predecessors, current methods for sequencing the mitochondrial and nuclear genomes have important limitations, and molecular diagnostic techniques are often manual and time consuming. However, recent advances in genomics and transcriptomics offer realistic solutions to these challenges. In this Review, we discuss the current genetic testing approach for PMDs and the opportunities that exist for increased use of whole-genome NGS of nuclear and mitochondrial DNA (mtDNA) in the clinical environment. We consider the possible role for long-read approaches in sequencing of mtDNA and in the identification of novel nuclear genomic causes of PMDs. We examine the expanding applications of RNA sequencing, including the detection of cryptic variants that affect splicing and gene expression and the interpretation of rare and novel mitochondrial transfer RNA variants.

A novel nonsense variant in MT-CO3 causes MELAS syndrome

Both mitochondrial and nuclear gene mutations can cause cytochrome c oxidase (COX, complex Ⅳ) dysfunction, leading to mitochondrial diseases. Although numerous diseases caused by defects of the COX subunits or COX assembly factors have been documented, clinical cases directly related to mitochondrial cytochrome c oxidase subunit 3 gene (MT-CO3) mutations are relatively rare. Here, we report a 47-year-old female patient presented with mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) syndrome. Muscle pathology revealed ragged-red fibres and remarkable COX-deficient muscle fibres. Muscle mitochondrial DNA sequencing analysis identified a novel MT-CO3 variant (m.9553G>A) that changed a highly conserved amino acid to a stop codon (p.Trp116*). This variant was heteroplasmic in multiple tissues, where the mutation load was 13% in oral epithelial cells, 89% in muscle samples, and not detectable in the peripheral blood lymphocytes. Single muscle fiber PCR analysis showed clear segregation of the mutation load with COX deficient fibres. Western blot analysis of the muscle samples revealed a significant decrease in the levels of COX1, COX2, COX3, COX4 and UQCRC2. COX respiration activity was remarkably reduced (58.84%) relative to the controls according to spectrophotometric assays. Taken together, our results indicated that this m.9553G>A variant may be responsible for the MELAS symdrome in the proband by affecting the stability and function of COX. The study expands the clinical and molecular spectrum of COX3-specific mitochondrial diseases.

Splice-Break: exploiting an RNA-seq splice junction algorithm to discover mitochondrial DNA deletion breakpoints and analyses of psychiatric disorders

Deletions in the 16.6 kb mitochondrial genome have been implicated in numerous disorders that often display muscular and/or neurological symptoms due to the high-energy demands of these tissues. We describe a catalogue of 4489 putative mitochondrial DNA (mtDNA) deletions, including their frequency and relative read rate, using a combinatorial approach of mitochondria-targeted PCR, next-generation sequencing, bioinformatics, post-hoc filtering, annotation, and validation steps. Our bioinformatics pipeline uses MapSplice, an RNA-seq splice junction detection algorithm, to detect and quantify mtDNA deletion breakpoints rather than mRNA splices. Analyses of 93 samples from postmortem brain and blood found (i) the 4977 bp ‘common deletion’ was neither the most frequent deletion nor the most abundant; (ii) brain contained significantly more deletions than blood; (iii) many high frequency deletions were previously reported in MitoBreak, suggesting they are present at low levels in metabolically active tissues and are not exclusive to individuals with diagnosed mitochondrial pathologies; (iv) many individual deletions (and cumulative metrics) had significant and positive correlations with age and (v) the highest deletion burdens were observed in major depressive disorder brain, at levels greater than Kearns–Sayre Syndrome muscle. Collectively, these data suggest the Splice-Break pipeline can detect and quantify mtDNA deletions at a high level of resolution.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

eKLIPse: a sensitive tool for the detection and quantification of mitochondrial DNA deletions from next-generation sequencing data

PURPOSE

Accurate detection of mitochondrial DNA (mtDNA) alterations is essential for the diagnosis of mitochondrial diseases. The development of high-throughput sequencing technologies has enhanced the detection sensitivity of mtDNA pathogenic variants, but the detection of mtDNA rearrangements, especially multiple deletions, is still poorly processed. Here, we present eKLIPse, a sensitive and specific tool allowing the detection and quantification of large mtDNA rearrangements from single and paired-end sequencing data.

METHODS

The methodology was first validated using a set of simulated data to assess the detection sensitivity and specificity, and second with a series of sequencing data from mitochondrial disease patients carrying either single or multiple deletions, related to pathogenic variants in nuclear genes involved in mtDNA maintenance.

RESULTS

eKLIPse provides the precise breakpoint positions and the cumulated percentage of mtDNA rearrangements at a given gene location with a detection sensitivity lower than 0.5% mutant. eKLIPse software is available either as a script to be integrated in a bioinformatics pipeline, or as user-friendly graphical interface to visualize the results through a Circos representation ( https://github.com/dooguypapua/eKLIPse ).

CONCLUSION

Thus, eKLIPse represents a useful resource to study the causes and consequences of mtDNA rearrangements, for further genotype/phenotype correlations in mitochondrial disorders.

Analytical parameters and validation of homopolymer detection in a pyrosequencing-based next generation sequencing system

Background

Current technologies in next-generation sequencing are offering high throughput reads at low costs, but still suffer from various sequencing errors. Although pyro- and ion semiconductor sequencing both have the advantage of delivering long and high quality reads, problems might occur when sequencing homopolymer-containing regions, since the repeating identical bases are going to incorporate during the same synthesis cycle, which leads to uncertainty in base calling. The aim of this study was to evaluate the analytical performance of a pyrosequencing-based next-generation sequencing system in detecting homopolymer sequences using homopolymer-preintegrated plasmid constructs and human DNA samples originating from patients with cystic fibrosis.

Results

In the plasmid system average correct genotyping was 95.8% in 4-mers, 87.4% in 5-mers and 72.1% in 6-mers. Despite the experienced low genotyping accuracy in 5- and 6-mers, it was possible to generate amplicons with more than a 90% adequate detection rate in every homopolymer tract. When homopolymers in the CFTR gene were sequenced average accuracy was 89.3%, but varied in a wide range (52.2 – 99.1%). In all but one case, an optimal amplicon-sequencing primer combination could be identified. In that single case (7A tract in exon 14 (c.2046_2052)), none of the tested primer sets produced the required analytical performance.

Conclusions

Our results show that pyrosequencing is the most reliable in case of 4-mers and as homopolymer length gradually increases, accuracy deteriorates. With careful primer selection, the NGS system was able to correctly genotype all but one of the homopolymers in the CFTR gene. In conclusion, we configured a plasmid test system that can be used to assess genotyping accuracy of NGS devices and developed an accurate NGS assay for the molecular diagnosis of CF using self-designed primers for amplification and sequencing.

Electronic supplementary material

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

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

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

A 2 bp deletion in the mitochondrial ATP 6 gene responsible for the NARP (neuropathy, ataxia, and retinitis pigmentosa) syndrome

Mitochondrial (mt) DNA-associated NARP (neurogenic muscle weakness, ataxia, and retinitis pigmentosa) syndrome is due to mutation in the MT-ATP6 gene. We report the case of a 18-year-old man who presented with deafness, a myoclonic epilepsy, muscle weakness since the age of 10 and further developed a retinitis pigmentosa and ataxia. The whole mtDNA analysis by next-generation sequencing revealed the presence of the 2 bp microdeletion m.9127-9128 del AT in the ATP6 gene at 82% heteroplasmy in muscle and to a lower load in blood (10-20%) and fibroblasts (50%). Using the patient's fibroblasts, we demonstrated a 60% reduction of the oligomycin-sensitive ATPase hydrolytic activity, a 40% decrease in the ATP synthesis and determination of the mitochondrial membrane potential using the fluorescent probe tetramethylrhodamine, ethyl ester indicated a significant reduction in oligomycin sensitivity. In conclusion, we demonstrated that this novel AT deletion in the ATP6 gene is pathogenic and responsible for the NARP syndrome.

Accurate and comprehensive analysis of single nucleotide variants and large deletions of the human mitochondrial genome in DNA and single cells

Massive parallel sequencing (MPS) can accurately quantify mitochondrial DNA (mtDNA) single nucleotide variants (SNVs), but no MPS methods are currently validated to simultaneously and accurately establish the breakpoints and frequency of large deletions at low heteroplasmic loads. Here we present the thorough validation of an MPS protocol to quantify the load of very low frequency, large mtDNA deletions in bulk DNA and single cells, along with SNV calling by standard methods. We used a set of well-characterized DNA samples, DNA mixes and single cells to thoroughly control the study. We developed a custom script for the detection of mtDNA rearrangements that proved to be more accurate in detecting and quantifying deletions than pre-existing tools. We also show that PCR conditions and primersets must be carefully chosen to avoid biases in the retrieved variants and an increase in background noise, and established a lower detection limit of 0.5% heteroplasmic load for large deletions, and 1.5 and 2% for SNVs, for bulk DNA and single cells, respectively. Finally, the analysis of different single cells provided novel insights into mtDNA cellular mosaicism.

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.

A Comparison and Integration of MiSeq and MinION Platforms for Sequencing Single Source and Mixed Mitochondrial Genomes

Single source and multiple donor (mixed) samples of human mitochondrial DNA were analyzed and compared using the MinION and the MiSeq platforms. A generalized variant detection strategy was employed to provide a cursory framework for evaluating the reliability and accuracy of mitochondrial sequences produced by the MinION. The feasibility of long-read phasing was investigated to establish its efficacy in quantitatively distinguishing and deconvolving individuals in a mixture. Finally, a proof-of-concept was demonstrated by integrating both platforms in a hybrid assembly that leverages solely mixture data to accurately reconstruct full mitochondrial genomes.

Massive parallel sequencing of human whole mitochondrial genomes with Ion Torrent technology: an optimized workflow for Anthropological and Population Genetics studies

Investigation of human mitochondrial DNA variation patterns and phylogeny has been extensively used in Anthropological and Population Genetics studies and sequencing the whole mitochondrial genome is progressively becoming the gold standard. Among the currently available massive parallel sequencing technologies, Ion Torrent™ semiconductor sequencing represents a promising approach for such studies. Nevertheless, an experimental protocol conceived to enable the achievement of both as high as possible yield and of the most homogeneous sequence coverage through the whole mitochondrial genome is still not available. The present work was thus aimed at improving the overall performance of whole mitochondrial genomes Ion Torrent™ sequencing, with special focus on the capability to obtain robust coverage and highly reliable variants calling. For this purpose, a series of cost-effective modifications in standard laboratory workflows was fine-tuned to optimize them for medium- and large-scale population studies. A total of 54 human samples were thus subjected to sequencing of the whole mitochondrial genome with the Ion Personal Genome Machine™ System in four distinct experiments and using Ion 314 chips. Seven of the selected samples were also characterized by means of conventional Sanger sequencing for the sake of comparison. Obtained results demonstrated that the implemented optimizations had definitely improved sequencing outputs in terms of both variants calling efficiency and coverage uniformity, enabling to setup an effective and accurate protocol for whole mitochondrial genome sequencing and a considerable reduction in experimental time consumption and sequencing costs.

Replication Errors Made During Oogenesis Lead to Detectable De Novo mtDNA Mutations in Zebrafish Oocytes with a Low mtDNA Copy Number

Of all pathogenic mitochondrial DNA (mtDNA) mutations in humans, ∼25% is de novo, although the occurrence in oocytes has never been directly assessed. We used next-generation sequencing to detect point mutations directly in the mtDNA of 3-15 individual mature oocytes and three somatic tissues from eight zebrafish females. Various statistical and biological filters allowed reliable detection of de novo variants with heteroplasmy ≥1.5%. In total, we detected 38 de novo base substitutions, but no insertions or deletions. These 38 de novo mutations were present in 19 of 103 mature oocytes, indicating that ∼20% of the mature oocytes carry at least one de novo mutation with heteroplasmy ≥1.5%. This frequency of de novo mutations is close to that deducted from the reported error rate of polymerase gamma, the mitochondrial replication enzyme, implying that mtDNA replication errors made during oogenesis are a likely explanation. Substantial variation in the mutation prevalence among mature oocytes can be explained by the highly variable mtDNA copy number, since we previously reported that ∼20% of the primordial germ cells have a mtDNA copy number of ≤73 and would lead to detectable mutation loads. In conclusion, replication errors made during oogenesis are an important source of de novo mtDNA base substitutions and their location and heteroplasmy level determine their significance.

Mitochondrial genome architecture in non-alcoholic fatty liver disease

Non-alcoholic fatty liver disease (NAFLD) is associated with mitochondrial dysfunction, a decreased liver mitochondrial DNA (mtDNA) content, and impaired energy metabolism. To understand the clinical implications of mtDNA diversity in the biology of NAFLD, we applied deep-coverage whole sequencing of the liver mitochondrial genomes. We used a multistage study design, including a discovery phase, a phenotype-oriented study to assess the mutational burden in patients with steatohepatitis at different stages of liver fibrosis, and a replication study to validate findings in loci of interest. We also assessed the potential protein-level impact of the observed mutations. To determine whether the observed changes are tissue-specific, we compared the liver and the corresponding peripheral blood entire mitochondrial genomes. The nuclear genes POLG and POLG2 (mitochondrial DNA polymerase-γ) were also sequenced. We observed that the liver mtDNA of patients with NAFLD harbours complex genomes with a significantly higher mutational (1.28-fold) rate and degree of heteroplasmy than in controls. The analysis of liver mitochondrial genomes of patients with different degrees of fibrosis revealed that the disease severity is associated with an overall 1.4-fold increase in mutation rate, including mutations in genes of the oxidative phosphorylation (OXPHOS) chain. Significant differences in gene and protein expression patterns were observed in association with the cumulative number of OXPHOS polymorphic sites. We observed a high degree of homology (∼98%) between the blood and liver mitochondrial genomes. A missense POLG p.Gln1236His variant was associated with liver mtDNA copy number. In conclusion, we have demonstrated that OXPHOS genes contain the highest number of hotspot positions associated with a more severe phenotype. The variability of the mitochondrial genomes probably originates from a common germline source; hence, it may explain a fraction of the 'missing heritability' of NAFLD. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Mitochondrial DNA heteroplasmy in the emerging field of massively parallel sequencing

Long an important and useful tool in forensic genetic investigations, mitochondrial DNA (mtDNA) typing continues to mature. Research in the last few years has demonstrated both that data from the entire molecule will have practical benefits in forensic DNA casework, and that massively parallel sequencing (MPS) methods will make full mitochondrial genome (mtGenome) sequencing of forensic specimens feasible and cost-effective. A spate of recent studies has employed these new technologies to assess intraindividual mtDNA variation. However, in several instances, contamination and other sources of mixed mtDNA data have been erroneously identified as heteroplasmy. Well vetted mtGenome datasets based on both Sanger and MPS sequences have found authentic point heteroplasmy in approximately 25% of individuals when minor component detection thresholds are in the range of 10-20%, along with positional distribution patterns in the coding region that differ from patterns of point heteroplasmy in the well-studied control region. A few recent studies that examined very low-level heteroplasmy are concordant with these observations when the data are examined at a common level of resolution. In this review we provide an overview of considerations related to the use of MPS technologies to detect mtDNA heteroplasmy. In addition, we examine published reports on point heteroplasmy to characterize features of the data that will assist in the evaluation of future mtGenome data developed by any typing method.

A bumpy ride on the diagnostic bench of massive parallel sequencing, the case of the mitochondrial genome

The advent of massive parallel sequencing (MPS) has revolutionized the field of human molecular genetics, including the diagnostic study of mitochondrial (mt) DNA dysfunction. The analysis of the complete mitochondrial genome using MPS platforms is now common and will soon outrun conventional sequencing. However, the development of a robust and reliable protocol is rather challenging. A previous pilot study for the re-sequencing of human mtDNA revealed an uneven coverage, affecting predominantly part of the plus strand. In an attempt to address this problem, we undertook a comparative study of standard and modified protocols for the Ion Torrent PGM system. We could not improve strand representation by altering the recommended shearing methodology of the standard workflow or omitting the DNA polymerase amplification step from the library construction process. However, we were able to associate coverage bias of the plus strand with a specific sequence motif. Additionally, we compared coverage and variant calling across technologies. The same samples were also sequenced on a MiSeq device which showed that coverage and heteroplasmic variant calling were much improved.

Whole-mitochondrial genome sequencing in primary open-angle glaucoma using massively parallel sequencing identifies novel and known pathogenic variants

PURPOSE

The aim of this study was to determine whether mutations in mitochondrial DNA play a role in high-pressure primary open-angle glaucoma (OMIM 137760) by analyzing new data from massively parallel sequencing of mitochondrial DNA.

METHODS

Glaucoma patients with high-tension primary open-angle glaucoma and ethnically matched and age-matched control subjects without glaucoma were recruited. The entire human mitochondrial genome was amplified in two overlapping fragments by long-range polymerase chain reaction and used as a template for massively parallel sequencing on an Ion Torrent Personal Genome Machine. All variants were confirmed by conventional Sanger sequencing.

RESULTS

Whole-mitochondrial genome sequencing was performed in 32 patients with primary open-angle glaucoma from India (n = 16) and Ireland (n = 16). In 16 of the 32 patients with primary open-angle glaucoma (50% of cases), there were 22 mitochondrial DNA mutations consisting of 7 novel mutations and 8 previously reported disease-associated sequence variants. Eight of 22 (36.4%) of the mitochondrial DNA mutations were in complex I mitochondrial genes.

CONCLUSION

Massively parallel sequencing using the Ion Torrent Personal Genome Machine with confirmation by Sanger sequencing detected a pathogenic mitochondrial DNA mutation in 50% of the primary open-angle glaucoma cohort. Our findings support the emerging concept that mitochondrial dysfunction results in the development of glaucoma and, more specifically, that complex I defects play a significant role in primary open-angle glaucoma pathogenesis.

Clonal expansion of early to mid-life mitochondrial DNA point mutations drives mitochondrial dysfunction during human ageing

Age-related decline in the integrity of mitochondria is an important contributor to the human ageing process. In a number of ageing stem cell populations, this decline in mitochondrial function is due to clonal expansion of individual mitochondrial DNA (mtDNA) point mutations within single cells. However the dynamics of this process and when these mtDNA mutations occur initially are poorly understood. Using human colorectal epithelium as an exemplar tissue with a well-defined stem cell population, we analysed samples from 207 healthy participants aged 17-78 years using a combination of techniques (Random Mutation Capture, Next Generation Sequencing and mitochondrial enzyme histochemistry), and show that: 1) non-pathogenic mtDNA mutations are present from early embryogenesis or may be transmitted through the germline, whereas pathogenic mtDNA mutations are detected in the somatic cells, providing evidence for purifying selection in humans, 2) pathogenic mtDNA mutations are present from early adulthood (<20 years of age), at both low levels and as clonal expansions, 3) low level mtDNA mutation frequency does not change significantly with age, suggesting that mtDNA mutation rate does not increase significantly with age, and 4) clonally expanded mtDNA mutations increase dramatically with age. These data confirm that clonal expansion of mtDNA mutations, some of which are generated very early in life, is the major driving force behind the mitochondrial dysfunction associated with ageing of the human colorectal epithelium.