Mitochondrial genome sequence analysis: a custom bioinformatics pipeline substantially improves Affymetrix MitoChip v2.0 call rate and accuracy

Mitochondrial genome study in blood of maternally inherited ALS cases

BACKGROUND

ALS is a heterogeneous disease in which different factors such as mitochondrial phenotypes act in combination with a genetic predisposition. This study addresses the question of whether homoplasmic (total mitochondrial genome of a sample is affected) and/or heteroplasmic mutations (wildtype and mutant mitochondrial DNA molecules coexist) might play a role in familial ALS. Blood was drawn from familial ALS patients with a possible maternal pattern of inheritance according to their pedigrees, which was compared to blood of ALS patients without maternal association as well as age-matched controls. In two cohorts, we analyzed the mitochondrial genome from whole blood or isolated white blood cells and platelets using a resequencing microarray (Affymetrix MitoChip v2.0) that is able to detect homoplasmic and heteroplasmic mitochondrial DNA mutations and allows the assessment of low-level heteroplasmy.

RESULTS

We identified an increase in homoplasmic ND5 mutations, a subunit of respiratory chain complex I, in whole blood of ALS patients that allowed maternal inheritance. This effect was more pronounced in patients with bulbar onset. Heteroplasmic mutations were significantly increased in different mitochondrial genes in platelets of patients with possible maternal inheritance. No increase of low-level heteroplasmy was found in maternal ALS patients.

CONCLUSION

Our results indicate a contribution of homoplasmic ND5 mutations to maternally associated ALS with bulbar onset. Therefore, it might be conceivable that specific maternally transmitted rather than randomly acquired mitochondrial DNA mutations might contribute to the disease process. This stands in contrast with observations from Alzheimer's and Parkinson's diseases showing an age-dependent accumulation of unspecific mutations in mitochondrial DNA.

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.

Baseline Characteristics of Mitochondrial DNA and Mutations Associated With Short-Term Posttreatment CD4+T-Cell Recovery in Chinese People With HIV

Background

Mitochondrial DNA (mtDNA) profiles and contributions of mtDNA variants to CD4+T-cell recovery in Euramerican people living with HIV (PLWH) may not be transferred to East-Asian PLWH, highlighting the need to consider more regional studies. We aimed to identify mtDNA characteristics and mutations that explain the variability of short-term CD4+T-cell recovery in East-Asian PLWH.

Method

Eight hundred fifty-six newly reported antiretroviral therapy (ART)-naïve Chinese PLWH from the Comparative HIV and Aging Research in Taizhou (CHART) cohort (Zhejiang Province, Eastern China) were enrolled. MtDNA was extracted from peripheral whole blood of those PLWH at HIV diagnosis, amplified, and sequenced using polymerase chain reaction and gene array. Characterization metrics such as mutational diversity and momentum were developed to delineate baseline mtDNA mutational patterns in ART-naïve PLWH. The associations between mtDNA genome-wide single nucleotide variants and CD4+T-cell recovery after short-term (within ~48 weeks) ART in 724 PLWH were examined using bootstrapping median regressions.

Results

Of 856 participants, 74.18% and 25.82% were male and female, respectively. The median age was 37 years; 94.51% were of the major Han ethnicity, and 69.04% and 28.62% were of the heterosexual and homosexual transmission, respectively. We identified 2,352 types of mtDNA mutations and mtDNA regions D-loop, ND5, CYB, or RNR1 with highest mutational diversity or volume. Female PLWH rather than male PLWH at the baseline showed remarkable age-related uptrends of momentum and mutational diversity as well as correlations between CD4+T <200 (cells/μl) and age-related uptrends of mutational diversity in many mtDNA regions. After adjustments of important sociodemographic and clinical variables, m.1005T>C, m.1824T>C, m.3394T>C, m.4491G>A, m.7828A>G, m.9814T>C, m.10586G>A, m.12338T>C, m.13708G>A, and m.14308T>C (at the Bonferroni-corrected significance) were negatively associated with short-term CD4+T-cell recovery whereas m.93A>G, m.15218A>G, and m.16399A>G were positively associated with short-term CD4+T-cell recovery.

Conclusion

Our baseline mtDNA characterization stresses the attention to East-Asian female PLWH at risk of CD4+T-cell loss-related aging and noncommunicable chronic diseases. Furthermore, mtDNA variants identified in regression analyses account for heterogeneity in short-term CD4+T-cell recovery of East-Asian PLWH. These results may help individualize the East-Asian immune recovery strategies under complicated HIV management caused by CD4+T-cell loss.

Mitochondrial DNA haplogroups and domain-specific neurocognitive performance in adults with HIV

Neurocognitive (NC) impairment (NCI) is an important cause of morbidity in persons with HIV (PWH). In the high-energy environment of the central nervous system, mitochondria contribute to neuroinflammation and aging, which may ultimately drive the pathogenesis of neurodegenerative diseases. Mitochondrial DNA (mtDNA) haplogroups are associated with health outcomes in PWH. For example, we previously observed less global NCI in Hispanic ancestry PWH having mtDNA haplogroup B. Another study reported increased NCI among PWH having African subhaplogroup L2a. We therefore analyzed NC domains in relation to these haplogroups in CNS HIV Antiretroviral Therapy Effects Research (CHARTER), a multi-site observational neuro-HIV study. Haplogroups were assigned using mtDNA sequence in 1016 PWH. Outcomes were NCI, defined by domain deficit score and mean T-scores (TS) for seven NC domains. Ancestry-stratified analyses of NC performance included Wilcoxon rank sum, χ2, and Fisher's exact tests. Multivariable regression adjusted for NC comorbidity, antiretroviral therapy use, and nadir CD4+ T cells. Among 98 Hispanic ancestry PWH, executive function, learning, and recall performance were better with haplogroup B (N = 17) than other haplogroups. With adjustment for covariates, haplogroup B remained associated with better executive function (p = 0.04) and recall TS (p = 0.03). PWH with haplogroup B had fewer impaired domains than other haplogroups (p < 0.01). Subhaplogroup L2a (N = 89) was associated with greater NCI in learning, recall, and working memory among 478 PWH of African ancestry, and had more impaired domains than other subhaplogroups (p < 0.01). These findings may inform risk stratification for NCI and studies to define mechanisms by which mtDNA variation may influence NCI in PWH.

Biocomplexity and Fractality in the Search of Biomarkers of Aging and Pathology: Mitochondrial DNA Profiling of Parkinson's Disease

Increasing evidence implicates mitochondrial dysfunction in the etiology of Parkinson's disease (PD). Mitochondrial DNA (mtDNA) mutations are considered a possible cause and this mechanism might be shared with the aging process and with other age-related neurodegenerative disorders such as Alzheimer's disease (AD). We have recently proposed a computerized method for mutated mtDNA characterization able to discriminate between AD and aging. The present study deals with mtDNA mutation-based profiling of PD. Peripheral blood mtDNA sequences from late-onset PD patients and age-matched controls were analyzed and compared to the revised Cambridge Reference Sequence (rCRS). The chaos game representation (CGR) method, modified to visualize heteroplasmic mutations, was used to display fractal properties of mtDNA sequences and fractal lacunarity analysis was applied to quantitatively characterize PD based on mtDNA mutations. Parameter β, from the hyperbola model function of our lacunarity method, was statistically different between PD and control groups when comparing mtDNA sequence frames corresponding to GenBank np 5713-9713. Our original method, based on CGR and lacunarity analysis, represents a useful tool to analyze mtDNA mutations. Lacunarity parameter β is able to characterize individual mutation profile of mitochondrial genome and could represent a promising index to discriminate between PD and aging.

Peripheral Blood Mitochondrial DNA Copy Number Obtained From Genome-Wide Genotype Data Is Associated With Neurocognitive Impairment in Persons With Chronic HIV Infection

BACKGROUND

Mitochondrial DNA (mtDNA) copy number varies by cell type and energy demands. Blood mtDNA copy number has been associated with neurocognitive function in persons without HIV. Low mtDNA copy number may indicate disordered mtDNA replication; high copy number may reflect a response to mitochondrial dysfunction. We hypothesized that blood mtDNA copy number estimated from genome-wide genotyping data is related to neurocognitive impairment (NCI) in persons with HIV.

METHODS

In the CNS HIV Antiretroviral Therapy Effects Research (CHARTER) study, peripheral blood mtDNA copy number was obtained from genome-wide genotyping data as a ratio of mtDNA single-nucleotide polymorphism probe intensities relative to nuclear DNA single-nucleotide polymorphisms. In a multivariable regression model, associations between mtDNA copy number and demographics, blood cell counts, and HIV disease and treatment characteristics were tested. Associations of mtDNA copy number with the global deficit score (GDS), GDS-defined NCI (GDS ≥ 0.5), and HIV-associated neurocognitive disorder (HAND) diagnosis were tested by logistic regression, adjusting for potential confounders.

RESULTS

Among 1010 CHARTER participants, lower mtDNA copy number was associated with longer antiretroviral therapy duration (P < 0.001), but not with d-drug exposure (P = 0.85). mtDNA copy number was also associated with GDS (P = 0.007), GDS-defined NCI (P < 0.001), and HAND (P = 0.002). In all analyses, higher mtDNA copy number was associated with poorer cognitive performance.

CONCLUSIONS

Higher mtDNA copy number estimated from peripheral blood genotyping was associated with worse neurocognitive performance in adults with HIV. These results suggest a connection between peripheral blood mtDNA and NCI, and may represent increased mtDNA replication in response to mitochondrial dysfunction.

Nuclear-Mitochondrial interactions influence susceptibility to HIV-associated neurocognitive impairment

HIV-associated neurocognitive impairment (NCI) is a term established to capture a wide spectrum of HIV related neurocognitive deficits ranging in severity from asymptomatic to dementia. The genetic underpinnings of this complex phenotype are incompletely understood. Mitochondrial function has long been thought to play a role in neurodegeneration, along with iron metabolism and transport. In this work, we aimed to characterize the interplay of mitochondrial DNA (mtDNA) haplogroup and nuclear genetic associations to NCI phenotypes in the CHARTER cohort, encompassing 1025 individuals of European-descent, African-descent, or admixed Hispanic. We first employed a polygenic modeling approach to investigate the global effect of previous marginally associated nuclear SNPs, and to examine how the polygenic effect of these SNPs is influenced by mtDNA haplogroups. We see evidence of a significant interaction between nuclear SNPs en masse and mtDNA haplogroups within European-descent and African-descent individuals. Subsequently, we performed an analysis of each SNP by mtDNA haplogroup, and detected significant interactions between two nuclear SNPs (rs17160128 and rs12460243) and European haplogroups. These findings, which require validation in larger cohorts, indicate a potential new role for nuclear-mitochondrial DNA interactions in susceptibility to NCI and shed light onto the pathophysiology of this neurocognitive phenotype.

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.

Recent Advances in Detecting Mitochondrial DNA Heteroplasmic Variations

The co-existence of wild-type and mutated mitochondrial DNA (mtDNA) molecules termed heteroplasmy becomes a research hot point of mitochondria. In this review, we listed several methods of mtDNA heteroplasmy research, including the enrichment of mtDNA and the way of calling heteroplasmic variations. At the present, while calling the novel ultra-low level heteroplasmy, high-throughput sequencing method is dominant while the detection limit of recorded mutations is accurate to 0.01% using the other quantitative approaches. In the future, the studies of mtDNA heteroplasmy may pay more attention to the single-cell level and focus on the linkage of mutations.

Biocomplexity and Fractality in the Search of Biomarkers of Aging and Pathology: Focus on Mitochondrial DNA and Alzheimer's Disease

Alzheimer’s disease (AD) represents one major health concern for our growing elderly population. It accounts for increasing impairment of cognitive capacity followed by loss of executive function in late stage. AD pathogenesis is multifaceted and difficult to pinpoint, and understanding AD etiology will be critical to effectively diagnose and treat the disease. An interesting hypothesis concerning AD development postulates a cause-effect relationship between accumulation of mitochondrial DNA (mtDNA) mutations and neurodegenerative changes associated with this pathology. Here we propose a computerized method for an easy and fast mtDNA mutations-based characterization of AD. The method has been built taking into account the complexity of living being and fractal properties of many anatomic and physiologic structures, including mtDNA. Dealing with mtDNA mutations as gaps in the nucleotide sequence, fractal lacunarity appears a suitable tool to differentiate between aging and AD. Therefore, Chaos Game Representation method has been used to display DNA fractal properties after adapting the algorithm to visualize also heteroplasmic mutations. Parameter β from our fractal lacunarity method, based on hyperbola model function, has been measured to quantitatively characterize AD on the basis of mtDNA mutations. Results from this pilot study to develop the method show that fractal lacunarity parameter β of mtDNA is statistically different in AD patients when compared to age-matched controls. Fractal lacunarity analysis represents a useful tool to analyze mtDNA mutations. Lacunarity parameter β is able to characterize individual mutation profile of mitochondrial genome and appears a promising index to discriminate between AD and aging.

European Mitochondrial DNA Haplogroups are Associated with Cerebrospinal Fluid Biomarkers of Inflammation in HIV Infection

Background:

Mitochondrial DNA (mtDNA) haplogroups are ancestry-related patterns of single-nucleotide polymorphisms that are associated with differential mitochondrial function in model systems, neurodegenerative diseases in HIV-negative populations, and chronic complications of HIV infection, including neurocognitive impairment. We hypothesized that mtDNA haplogroups are associated with neuroinflammation in HIV-infected adults.

Methods:

CNS HIV Antiretroviral Therapy Effects Research (CHARTER) is a US-based observational study of HIV-infected adults who underwent standardized neurocognitive assessments. Participants who consented to DNA collection underwent whole blood mtDNA sequencing, and a subset also underwent lumbar puncture. IL-6, IL-8, TNF-α (high-sensitivity), and IP-10 were measured in cerebrospinal fluid (CSF) by immunoassay. Multivariable regression of mtDNA haplogroups and log-transformed CSF biomarkers were stratified by genetic ancestry using whole-genome nuclear DNA genotyping (European [EA], African [AA], or Hispanic ancestry [HA]), and adjusted for age, sex, antiretroviral therapy (ART), detectable CSF HIV RNA, and CD4 nadir. A total of 384 participants had both CSF cytokine measures and genetic data (45% EA, 44% AA, 11% HA, 22% female, median age 43 years, 74% on ART).

Results:

In analyses stratified by the 3 continental ancestry groups, no haplogroups were significantly associated with the 4 biomarkers. In the subgroup of participants with undetectable plasma HIV RNA on ART, European haplogroup H participants had significantly lower CSF TNF-α (P = 0.001).

Conclusions:

Lower CSF TNF-α may indicate lower neuroinflammation in the haplogroup H participants with well-controlled HIV on ART.

A Clinical, Neuropathological and Genetic Study of Homozygous A467T POLG-Related Mitochondrial Disease

Mutations in the nuclear gene POLG (encoding the catalytic subunit of DNA polymerase gamma) are an important cause of mitochondrial disease. The most common POLG mutation, A467T, appears to exhibit considerable phenotypic heterogeneity. The mechanism by which this single genetic defect results in such clinical diversity remains unclear. In this study we evaluate the clinical, neuropathological and mitochondrial genetic features of four unrelated patients with homozygous A467T mutations. One patient presented with the severe and lethal Alpers-Huttenlocher syndrome, which was confirmed on neuropathology, and was found to have a depletion of mitochondrial DNA (mtDNA). Of the remaining three patients, one presented with mitochondrial encephalomyopathy, lactic acidosis and stroke-like episodes (MELAS), one with a phenotype in the Myoclonic Epilepsy, Myopathy and Sensory Ataxia (MEMSA) spectrum and one with Sensory Ataxic Neuropathy, Dysarthria and Ophthalmoplegia (SANDO). All three had secondary accumulation of multiple mtDNA deletions. Complete sequence analysis of muscle mtDNA using the MitoChip resequencing chip in all four cases demonstrated significant variation in mtDNA, including a pathogenic MT-ND5 mutation in one patient. These data highlight the variable and overlapping clinical and neuropathological phenotypes and downstream molecular defects caused by the A467T mutation, which may result from factors such as the mtDNA genetic background, nuclear genetic modifiers and environmental stressors.

Mitochondrial DNA Haplogroups and Neurocognitive Impairment During HIV Infection

BACKGROUND

Neurocognitive impairment (NCI) remains an important complication in persons infected with human immunodeficiency virus (HIV). Ancestry-related mitochondrial DNA (mtDNA) haplogroups have been associated with outcomes of HIV infection and combination antiretroviral therapy (CART), and with neurodegenerative diseases. We hypothesize that mtDNA haplogroups are associated with NCI in HIV-infected adults and performed a genetic association study in the CNS HIV Antiretroviral Therapy Effects Research (CHARTER) cohort.

METHODS

CHARTER is an observational study of ambulatory HIV-infected adults. Haplogroups were assigned using mtDNA sequence, and principal components were derived from ancestry-informative nuclear DNA variants. Outcomes were cross-sectional global deficit score (GDS) as a continuous measure, GDS impairment (GDS ≥ 0.50), and HIV-associated neurocognitive disorder (HAND) using international criteria. Multivariable models were adjusted for comorbidity status (incidental vs contributing), current CART, plasma HIV RNA, reading ability, and CD4 cell nadir.

RESULTS

Haplogroups were available from 1027 persons; median age 43 years, median CD4 nadir 178 cells/mm(3), 72% on CART, and 46% with HAND. The 102 (9.9%) persons of genetically determined admixed Hispanic ancestry had more impairment by GDS or HAND than persons of European or African ancestry (P < .001 for all). In multivariate models including persons of admixed Hispanic ancestry, those with haplogroup B had lower GDS (β = -0.34; P = .008) and less GDS impairment (odds ratio = 0.16; 95% confidence interval, .04, .63; P = .009) than other haplogroups. There were no significant haplogroup associations among persons of European or African ancestry.

CONCLUSIONS

In these mostly CART-treated persons, mtDNA haplogroup B was associated with less NCI among persons of genetically determined Hispanic ancestry. mtDNA variation may represent an ancestry-specific factor influencing NCI in HIV-infected persons.

Contribution of non-reference alleles in mtDNA of Alzheimer's disease patients

Many observations suggest that mutations of mitochondrial DNA (mtDNA) could be responsible for the neurodegenerative changes of Alzheimer's disease (AD). Here we examined the signal intensity of the four alleles of each mtDNA nucleotide position (np) in whole blood of AD patients and age-matched controls using MitoChip v2.0 array. Our analysis identified 270 significantly different nps which, with one exception, showed an increased contribution of non-reference alleles in AD patients. Principal component analysis (PCA) and cluster analysis showed that five of these nps could discriminate AD from control subjects with 80% of cases correctly classified. Our data support the hypothesis of mtDNA alterations as an important factor in the etiology of AD.

The case for the continuing use of the revised Cambridge Reference Sequence (rCRS) and the standardization of notation in human mitochondrial DNA studies

Since the determination in 1981 of the sequence of the human mitochondrial DNA (mtDNA) genome, the Cambridge Reference Sequence (CRS), has been used as the reference sequence to annotate mtDNA in molecular anthropology, forensic science and medical genetics. The CRS was eventually upgraded to the revised version (rCRS) in 1999. This reference sequence is a convenient device for recording mtDNA variation, although it has often been misunderstood as a wild-type (WT) or consensus sequence by medical geneticists. Recently, there has been a proposal to replace the rCRS with the so-called Reconstructed Sapiens Reference Sequence (RSRS). Even if it had been estimated accurately, the RSRS would be a cumbersome substitute for the rCRS, as the new proposal fuses--and thus confuses--the two distinct concepts of ancestral lineage and reference point for human mtDNA. Instead, we prefer to maintain the rCRS and to report mtDNA profiles by employing the hitherto predominant circumfix style. Tree diagrams could display mutations by using either the profile notation (in conventional short forms where appropriate) or in a root-upwards way with two suffixes indicating ancestral and derived nucleotides. This would guard against misunderstandings about reporting mtDNA variation. It is therefore neither necessary nor sensible to change the present reference sequence, the rCRS, in any way. The proposed switch to RSRS would inevitably lead to notational chaos, mistakes and misinterpretations.

Prevalence of rare mitochondrial DNA mutations in mitochondrial disorders

Background

Mitochondrial DNA (mtDNA) diseases are rare disorders whose prevalence is estimated around 1 in 5000. Patients are usually tested only for deletions and for common mutations of mtDNA which account for 5–40% of cases, depending on the study. However, the prevalence of rare mtDNA mutations is not known.

Methods

We analysed the whole mtDNA in a cohort of 743 patients suspected of manifesting a mitochondrial disease, after excluding deletions and common mutations. Both heteroplasmic and homoplasmic variants were identified using two complementary strategies (Surveyor and MitoChip). Multiple correspondence analyses followed by hierarchical ascendant cluster process were used to explore relationships between clinical spectrum, age at onset and localisation of mutations.

Results

7.4% of deleterious mutations and 22.4% of novel putative mutations were identified. Pathogenic heteroplasmic mutations were more frequent than homoplasmic mutations (4.6% vs 2.8%). Patients carrying deleterious mutations showed symptoms before 16 years of age in 67% of cases. Early onset disease (<1 year) was significantly associated with mutations in protein coding genes (mainly in complex I) while late onset disorders (>16 years) were associated with mutations in tRNA genes. MTND5 and MTND6 genes were identified as ‘hotspots’ of mutations, with Leigh syndrome accounting for the large majority of associated phenotypes.

Conclusions

Rare mitochondrial DNA mutations probably account for more than 7.4% of patients with respiratory chain deficiency. This study shows that a comprehensive analysis of mtDNA is essential, and should include young children, for an accurate diagnosis that is now accessible with the development of next generation sequencing technology.

Mitochondrial disease genetic diagnostics: optimized whole-exome analysis for all MitoCarta nuclear genes and the mitochondrial genome

Discovering causative genetic variants in individual cases of suspected mitochondrial disease requires interrogation of both the mitochondrial (mtDNA) and nuclear genomes. Whole-exome sequencing can support simultaneous dual-genome analysis, although currently available capture kits do not target the mtDNA genome and provide insufficient capture for some nuclear-encoded mitochondrial genes. To optimize interrogation of nuclear and mtDNA genes relevant to mitochondrial biology and disease, a custom SureSelect "Mito-Plus" whole-exome library was formulated by blending RNA "baits" from three separate designs: (A) Agilent Technologies SureSelectXT 50 Mb All Exon PLUS Targeted Enrichment Kit, (B) 16-gene nuclear panel targeting sequences for known MitoCarta proteins not included in the 50 Mb All Exon design, and (C) sequences targeting the entire mtDNA genome. The final custom formulations consisted of a 1:1 ratio of nuclear baits to which a 1 to 1,000-fold diluted ratio of mtDNA genome baits were blended. Patient sample capture libraries were paired-end sequenced on an Illumina HiSeq 2000 system using v3.0 SBS chemistry. mtDNA genome coverage varied depending on the mtDNA:nuclear blend ratio, where a 1:100 ratio provided optimal dual-genome coverage with 10X coverage for over 97.5% of all targeted nuclear regions and 1,000X coverage for 99.8% of the mtDNA genome. mtDNA mutations were reliably detected to at least an 8% heteroplasmy level, as discriminated both from sequencing errors and potential contamination from nuclear mtDNA transcripts (Numts). The "1:100 Mito-Plus Whole-Exome" Agilent capture kit offers an optimized tool for whole-exome analysis of nuclear and mtDNA genes relevant to the diagnostic evaluation of mitochondrial disease.