Background sequence characteristics influence the occurrence and severity of disease-causing mtDNA mutations

Mitochondrial DNA disease discovery through evaluation of genotype and phenotype data: The Solve-RD experience

The diagnosis of mitochondrial DNA (mtDNA) diseases remains challenging with next-generation sequencing, where bioinformatic analysis is usually more focused on the nuclear genome. We developed a workflow for the evaluation of mtDNA diseases and applied it in a large European rare disease cohort (Solve-RD). A semi-automated bioinformatic pipeline with MToolBox was used to filter the unsolved Solve-RD cohort for rare mtDNA variants after validating this pipeline on exome datasets of 42 individuals previously diagnosed with mtDNA variants. Variants were filtered based on blood heteroplasmy levels (≥1%) and reported association with disease. Overall, 10,157 exome and genome datasets from 9,923 affected individuals from 9,483 families within Solve-RD met the quality inclusion criteria. 136 mtDNA variants in 135 undiagnosed individuals were prioritized using the filtering approach. A focused MitoPhen-based phenotype similarity scoring method was tested in a separate genetically diagnosed "phenotype test cohort" consisting of nuclear gene and mtDNA diseases using a receiving operator characteristic evaluation. We applied the MitoPhen-based phenotype similarity score of >0.3, which was highly sensitive for detecting mtDNA diseases in the phenotype test cohort, to the filtered cohort of 135 undiagnosed individuals. This aided the prioritization of 34 out of 37 (92%) individuals who received confirmed and likely causative mtDNA disease diagnoses. The phenotypic evaluation was limited by the quality of input data in some individuals. The overall pipeline led to an additional diagnostic yield of 0.4% in a cohort where mitochondrial disease was not initially suspected. This highlights the value of our mtDNA analysis pipeline in diverse datasets.

Mitochondrial DNA variation and intervertebral disc degeneration: a genotypic analysis in a South African cohort

BACKGROUND

Non-communicable diseases are multifactorial in that they can be caused by genetic factors, age, sex and poor lifestyle choices. They are estimated to account for 71% of deaths globally with 80% of these deaths occurring in low- and middle-income countries. This is particularly true for Intervertebral Disc Degeneration associated with mitochondrial dysfunction. Interestingly, mitochondrial dysfunction can arise from mutations in both the nuclear and the mitochondrial genomes. The present study, therefore, aimed to determine if there is an association between mitochondrial DNA mutations associated with mitochondrial dysfunction and disc degeneration in a South African cohort, and in addition, generate genetic data for understudied mutations in African populations.

METHODS AND RESULTS

Mutations were selected using a systematic literature review. DNA was collected using buccal swabs and extracted using a standard salt-lysis protocol. Mass-array genotyping was done for previously reported as well as novel mutations. GenAlEx (version 6.5), RStudio and SHEsis were used for statistical analyses. Although no significant associations were found, the identified polymorphic mutations C16223T, A10398G and A8536G were found to have higher mutant allele frequencies in case individuals indicating that had a larger cohort been used, significance may have been observed.

CONCLUSIONS

This study was able to generate genotypic information for a South African cohort for both reported and understudied mutations. Furthermore, the identification of higher mutant allele frequencies for C16223T, A10398G and A8536G highlights the importance of considering these mutations in future studies using a larger cohort.

Novel pathogenic mtDNA variants in Chinese children with neurological mitochondrial disorders

OBJECTIVE

Pathogenic variations in the mitochondrial genome are tightly linked to neurological mitochondrial disorders in children. However, the mutation spectrum of mitochondrial DNA (mtDNA) in the Chinese population remains incomplete. Therefore, the primary objective of our study was to comprehensively characterize pathogenic mtDNA variants in Chinese children with mitochondrial disorders at clinical, molecular, and functional levels.

METHODS

Between February 2019 and September 2023, we analyzed pathogenic mtDNA variants in a cohort of over 600 Chinese children suspected of having mitochondrial disorders. Whole-exome sequencing (WES) and whole-mtDNA sequencing were performed on the cohort.

RESULTS

We identified 54 pathogenic or likely pathogenic mtDNA variants in 227 Chinese children with neurological mitochondrial disorders. Among the eight novel heteroplasmic variants detected in seven patients, in silico analyses suggested likely pathogenic features. Functional analyses using either primary fibroblasts or cybrid cells carrying different mutant loads of mtDNA variants showed impaired mitochondrial respiration, ATP generation, and mitochondrial membrane potential in five of the eight novel variants, including m.4275G>A, m.10407G>A, m.5828G>A, m.3457G>A, and m.13112T>C. The m.8427T>C variant was identified as a rare polymorphism because, despite being located at MT-ATP8, it does not affect both the assembly and activity of mitochondrial complex V in cells carrying homoplasmic m.8427T>C variation. Transcriptome profiling further confirmed the pathogenic contributions of these five variants by altering mitochondrial pathways.

CONCLUSION

In summary, we revisited the mtDNA mutation spectrum in Chinese children with mitochondrial disorders, and identified five novel pathogenic mtDNA variants with functional verification that are related to neurological mitochondrial disorders in children.

Mitochondrial genome variants associated with amyotrophic lateral sclerosis and their haplogroup distribution

INTRODUCTION/AIMS

Amyotrophic lateral sclerosis (ALS) may be familial or sporadic, and twin studies have revealed that even sporadic forms have a significant genetic component. Variants in 55 nuclear genes have been associated with ALS and although mitochondrial dysfunction is observed in ALS, variants in mitochondrial genomes (mitogenomes) have not yet been tested for association with ALS. The aim of this study was to determine whether mitogenome variants are associated with ALS.

METHODS

We conducted a genome-wide association study (GWAS) in mitogenomes of 1965 ALS patients and 2547 controls.

RESULTS

We identified 51 mitogenome variants with p values <10-7, of which 13 had odds ratios (ORs) >1, in genes RNR1, ND1, CO1, CO3, ND5, ND6, and CYB, while 38 variants had OR <1 in genes RNR1, RNA2, ND1, ND2, CO2, ATP8, ATP6, CO3, ND3, ND4, ND5, ND6, and CYB. The frequencies of haplogroups H, U, and L, the most frequent in our ALS data set, were the same in different onset sites (bulbar, limb, spinal, and axial). Also, intra-haplogroup GWAS revealed unique ALS-associated variants in haplogroups L and U.

DISCUSSION

Our study shows that mitogenome single nucleotide variants (SNVs) are associated with ALS and suggests that these SNVs could be included in routine genetic testing for ALS and that mitochondrial replacement therapy has the potential to serve as a basis for ALS treatment.

Spectrum of rare and common mitochondrial DNA variations from 1029 whole genomes of self-declared healthy individuals from India

Mitochondrial disorders are a class of heterogeneous disorders caused by genetic variations in the mitochondrial genome (mtDNA) as well as the nuclear genome. The spectrum of mtDNA variants remains unexplored in the Indian population. In the present study, we have cataloged 2689 high confidence single nucleotide variants, small insertions and deletions in mtDNA in 1029 healthy Indian individuals. We found a major proportion (76.5 %) of the variants being rare (AF<=0.005) in the studied population. Intriguingly, we found two 'confirmed' pathogenic variants (m.1555 A>G and m.14484 T>C) with a frequency of ∼1 in 250 individuals in our dataset. The high carrier frequency underscores the need for screening of the mtDNA pathogenic mutations in newborns in India. Interestingly, our analysis also revealed 202 variants in our dataset which have been 'reported' in disease cases as per the MITOMAP database. Additionally, we found the frequency of haplogroup M (52.2 %) to be the highest among all the 18 top-level haplogroups found in our dataset. In comparison to the global population datasets, 20 unique mtDNA variants are found in the Indian population. We hope the whole genome sequencing based compendium of mtDNA variants along with their allele frequencies and heteroplasmy levels in the Indian population will drive additional genome scale studies for mtDNA. Furthermore, the identification of clinically relevant variants in our dataset will aid in better clinical interpretation of the variants in mitochondrial disorders.

mtDNA "nomenclutter" and its consequences on the interpretation of genetic data

Population-based studies of human mitochondrial genetic diversity often require the classification of mitochondrial DNA (mtDNA) haplotypes into more than 5400 described haplogroups, and further grouping those into hierarchically higher haplogroups. Such secondary haplogroup groupings (e.g., "macro-haplogroups") vary across studies, as they depend on the sample quality, technical factors of haplogroup calling, the aims of the study, and the researchers' understanding of the mtDNA haplogroup nomenclature. Retention of historical nomenclature coupled with a growing number of newly described mtDNA lineages results in increasingly complex and inconsistent nomenclature that does not reflect phylogeny well. This "clutter" leaves room for grouping errors and inconsistencies across scientific publications, especially when the haplogroup names are used as a proxy for secondary groupings, and represents a source for scientific misinterpretation. Here we explore the effects of phylogenetically insensitive secondary mtDNA haplogroup groupings, and the lack of standardized secondary haplogroup groupings on downstream analyses and interpretation of genetic data. We demonstrate that frequency-based analyses produce inconsistent results when different secondary mtDNA groupings are applied, and thus allow for vastly different interpretations of the same genetic data. The lack of guidelines and recommendations on how to choose appropriate secondary haplogroup groupings presents an issue for the interpretation of results, as well as their comparison and reproducibility across studies. To reduce biases originating from arbitrarily defined secondary nomenclature-based groupings, we suggest that future updates of mtDNA phylogenies aimed for the use in mtDNA haplogroup nomenclature should also provide well-defined and standardized sets of phylogenetically meaningful algorithm-based secondary haplogroup groupings such as "macro-haplogroups", "meso-haplogroups", and "micro-haplogroups". Ideally, each of the secondary haplogroup grouping levels should be informative about different human population history events. Those phylogenetically informative levels of haplogroup groupings can be easily defined using TreeCluster, and then implemented into haplogroup callers such as HaploGrep3. This would foster reproducibility across studies, provide a grouping standard for population-based studies, and reduce errors associated with haplogroup nomenclatures in future studies.

Mitochondrial DNA disorders in neuromuscular diseases in diverse populations

Neuromuscular features are common in mitochondrial DNA (mtDNA) disorders. The genetic architecture of mtDNA disorders in diverse populations is poorly understood. We analysed mtDNA variants from whole-exome sequencing data in neuromuscular patients from South Africa, Brazil, India, Turkey and Zambia. In 998 individuals, there were two definite diagnoses, two possible diagnoses and eight secondary findings. Surprisingly, common pathogenic mtDNA variants found in people of European ancestry were very rare. Whole-exome or -genome sequencing from undiagnosed patients with neuromuscular symptoms should be re-analysed for mtDNA variants, but the landscape of pathogenic mtDNA variants differs around the world.

Feasibility and impact of haplogroup matching for mitochondrial replacement treatment

Mitochondrial replacement technology (MRT) aims to reduce the risk of serious disease in children born to women who carry pathogenic mitochondrial DNA (mtDNA) variants. By transplanting nuclear genomes from eggs of an affected woman to enucleated eggs from an unaffected donor, MRT creates new combinations of nuclear and mtDNA. Based on sets of shared sequence variants, mtDNA is classified into ~30 haplogroups. Haplogroup matching between egg donors and women undergoing MRT has been proposed as a means of reducing mtDNA sequence divergence between them. Here we investigate the potential effect of mtDNA haplogroup matching on clinical delivery of MRT and on mtDNA sequence divergence between donor/recipient pairs. Our findings indicate that haplogroup matching would limit the availability of egg donors such that women belonging to rare haplogroups may have to wait > 4 years for treatment. Moreover, we find that intra-haplogroup sequence variation is frequently within the range observed between randomly matched mtDNA pairs. We conclude that haplogroup matching would restrict the availability of MRT, without necessarily reducing mtDNA sequence divergence between donor/recipient pairs.

Mitochondrial heteroplasmic shifts reveal a positive selection of breast cancer

BACKGROUND

Breast cancer is, despite screening, not always detected early enough and is together with other tumor types known to shed genetic information in circulation. Unlike single-copy nuclear DNA, mitochondrial DNA (mtDNA) copies range from 100s to 10,000s per cell, thus providing a potentially alternative to identify potential missing cancer information in circulation at an early stage.

METHODS

To characterize mitochondrial mutation landscapes in breast cancer, whole mtDNA sequencing and bioinformatics analyses were performed on 86 breast cancer biopsies and 50 available matched baseline cancer-free whole blood samples from the same individuals, selected from a cohort of middle-aged women in Sweden. To determine whether the mutations can be detected in blood plasma prior to cancer diagnosis, we further designed a nested case-control study (n = 663) and validated the shortlisted mutations using droplet digital PCR.

RESULTS

We detected different mutation landscapes between biopsies and matched whole blood samples. Compared to whole blood samples, mtDNA from biopsies had higher heteroplasmic mutations in the D-loop region (P = 0.02), RNR2 (P = 0.005), COX1 (P = 0.037) and CYTB (P = 0.006). Furthermore, the germline mtDNA mutations had higher heteroplasmy level than the lost (P = 0.002) and de novo mutations (P = 0.04). The nonsynonymous to synonymous substitution ratio (dN/dS) was higher for the heteroplasmic mutations (P = 7.25 × 10-12) than that for the homoplasmic mutations, but the de novo (P = 0.06) and lost mutations (P = 0.03) had lower dN/dS than the germline mutations. Interestingly, we found that the critical regions for mitochondrial transcription: MT-HSP1 (odds ratio [OR]: 21.41), MT-TFH (OR: 7.70) and MT-TAS2 (OR: 3.62), had significantly higher heteroplasmic mutations than the rest of the D-loop sub-regions. Finally, we found that the presence of mt.16093T > C mutation increases 67% risk of developing breast cancer.

CONCLUSIONS

Our findings show that mitochondrial genetic landscape changes during cancer pathogenesis and positive selection of mtDNA heteroplasmic mutations in breast cancer. Most importantly, the mitochondrial mutations identified in biopsies can be traced back in matched plasma samples and could potentially be used as early breast cancer diagnostic biomarkers.

Profiling mitochondrial DNA mutations in tumors and circulating extracellular vesicles of triple-negative breast cancer patients for potential biomarker development

Early detection and recurrence prediction are challenging in triple-negative breast cancer (TNBC) patients. We aimed to develop mitochondrial DNA (mtDNA)-based liquid biomarkers to improve TNBC management. Mitochondrial genome (MG) enrichment and next-generation sequencing mapped the entire MG in 73 samples (64 tissues and 9 extracellular vesicles [EV] samples) from 32 metastatic TNBCs. We measured mtDNA and cardiolipin (CL) contents, NDUFB8, and SDHB protein expression in tumors and in corresponding circulating EVs. We identified 168 nonsynonymous mtDNA mutations, with 73% (123/186) coding and 27% (45/168) noncoding in nature. Twenty percent of mutations were nucleotide transversions. Respiratory complex I (RCI) was the key target, which harbored 44% (74/168) of the overall mtDNA mutations. A panel of 11 hotspot mtDNA mutations was identified among 19%-38% TNBCs, which were detectable in the serum-derived EVs with 82% specificity. Overall, 38% of the metastatic tumor-signature mtDNA mutations were traceable in the EVs. An appreciable number of mtDNA mutations were homoplasmic (18%, 31/168), novel (14%, 23/168), and potentially pathogenic (9%, 15/168). The overall and RCI-specific mtDNA mutational load was higher in women with African compared to European ancestry accompanied by an exclusive abundance of respiratory complex (RC) protein NDUFB8 (RCI) and SDHB (RCII) therein. Increased mtDNA (p < 0.0001) content was recorded in both tumors and EVs along with an abundance of CL (p = 0.0001) content in the EVs. Aggressive tumor-signature mtDNA mutation detection and measurement of mtDNA and CL contents in the EVs bear the potential to formulate noninvasive early detection and recurrence prediction strategies.

Association of mitochondrial DNA variation with high myopia in a Han Chinese population

High myopia (HM), which is characterized by oxidative stress, is one of the leading causes of visual impairment and blindness across the world. Family and population genetic studies have uncovered nuclear-genome variants in proteins functioned in the mitochondria. However, whether mitochondrial DNA mutations are involved in HM remains unexplored. Here, we performed the first large-scale whole-mitochondrial genome study in 9613 HM cases and 9606 control subjects of Han Chinese ancestry for identifying HM-associated mitochondrial variants. The single-variant association analysis identified nine novel genetic variants associated with HM reaching the entire mitochondrial wide significance level, including rs370378529 in ND2 with an odds ratio (OR) of 5.25. Interestingly, eight out of nine variants were predominantly located in related sub-haplogroups, i.e. m.5261G > A in B4b1c, m.12280A > G in G2a4, m.7912G > A in D4a3b, m.94G > A in D4e1, m.14857 T > C in D4e3, m.14280A > G in D5a2, m.16272A > G in G2a4, m.8718A > G in M71 and F1a3, indicating that the sub-haplogroup background can increase the susceptible risk for high myopia. The polygenic risk score analysis of the target and validation cohorts indicated a high accuracy for predicting HM with mtDNA variants (AUC = 0.641). Cumulatively, our findings highlight the critical roles of mitochondrial variants in untangling the genetic etiology of HM.

Spectrum of germline and somatic mitochondrial DNA variants in Tuberous Sclerosis Complex

Tuberous Sclerosis Complex (TSC) is caused by loss of function variants in either TSC1 or TSC2 and is characterized by broad phenotypic heterogeneity. Currently, there is limited knowledge regarding the role of the mitochondrial genome (mtDNA) in TSC pathogenesis. In this study, we aimed to determine the prevalence and spectrum of germline and somatic mtDNA variants in TSC and identify potential disease modifiers. Analysis of mtDNA amplicon massively parallel sequencing (aMPS) data, off-target mtDNA from whole-exome sequencing (WES), and/or qPCR, revealed mtDNA alterations in 270 diverse tissues (139 TSC-associated tumors and 131 normal tissue samples) from 199 patients and six healthy individuals. Correlation of clinical features to mtDNA variants and haplogroup analysis was done in 102 buccal swabs (age: 20-71 years). No correlation was found between clinical features and either mtDNA variants or haplogroups. No pathogenic variants were identified in the buccal swab samples. Using in silico analysis, we identified three predicted pathogenic variants in tumor samples: MT-ND4 (m.11742G>A, p. Cys328Tyr, VAF: 43%, kidney angiomyolipoma), MT-CYB (m.14775T>C, p. Leu10Pro, VAF: 43%, LAM abdominal tumor) and MT-CYB (m.15555C>T, p. Pro270Leu, VAF: 7%, renal cell carcinoma). Large deletions of the mitochondrial genome were not detected. Analysis of tumors from 23 patients with corresponding normal tissue did not reveal any recurrent tumor-associated somatic variants. The mtDNA/gDNA ratio between tumors and corresponding normal tissue was also unchanged. Overall, our findings demonstrate that the mitochondrial genome is highly stable across tissues and within TSC-associated tumors.

Analysis of Human Clinical Mutations of Mitochondrial ND1 in a Bacterial Model System for Complex I

The most common causes of mitochondrial dysfunction and disease include mutations in subunits and assembly factors of Complex I. Numerous mutations in the mitochondrial gene ND1 have been identified in humans. Currently, a bacterial model system provides the only method for rapid construction and analysis of mutations in homologs of human ND1. In this report, we have identified nine mutations in human ND1 that are reported to be pathogenic and are located at subunit interfaces. Our hypothesis was that these mutations would disrupt Complex I assembly. Seventeen mutations were constructed in the homologous nuoH gene in an E. coli model system. In addition to the clinical mutations, alanine substitutions were constructed in order to distinguish between a deleterious effect from the introduction of the mutant residue and the loss of the original residue. The mutations were moved to an expression vector containing all thirteen genes of the E. coli nuo operon coding for Complex I. Membrane vesicles were prepared and rates of deamino-NADH oxidase activity and proton translocation were measured. Samples were also tested for assembly by native gel electrophoresis and for expression of NuoH by immunoblotting. A range of outcomes was observed: Mutations at four of the sites allow normal assembly with moderate activity (50−76% of wild type). Mutations at the other sites disrupt assembly and/or activity, and in some cases the outcomes depend upon the amino acid introduced. In general, the outcomes are consistent with the proposed pathogenicity in humans.

Precision mitochondrial medicine

Mitochondria play a key role in cell homeostasis as a major source of intracellular energy (adenosine triphosphate), and as metabolic hubs regulating many canonical cell processes. Mitochondrial dysfunction has been widely documented in many common diseases, and genetic studies point towards a causal role in the pathogenesis of specific late-onset disorder. Together this makes targeting mitochondrial genes an attractive strategy for precision medicine. However, the genetics of mitochondrial biogenesis is complex, with over 1,100 candidate genes found in two different genomes: the nuclear DNA and mitochondrial DNA (mtDNA). Here, we review the current evidence associating mitochondrial genetic variants with distinct clinical phenotypes, with some having clear therapeutic implications. The strongest evidence has emerged through the investigation of rare inherited mitochondrial disorders, but genome-wide association studies also implicate mtDNA variants in the risk of developing common diseases, opening to door for the incorporation of mitochondrial genetic variant analysis in population disease risk stratification.

Modulating mitochondrial DNA mutations: factors shaping heteroplasmy in the germ line and somatic cells

Until recently it was thought that most humans only harbor one type of mitochondrial DNA (mtDNA), however, deep sequencing and single-cell analysis has shown the converse - that mixed populations of mtDNA (heteroplasmy) are the norm. This is important because heteroplasmy levels can change dramatically during transmission in the female germ line, leading to high levels causing severe mitochondrial diseases. There is also emerging evidence that low level mtDNA mutations contribute to common late onset diseases such as neurodegenerative disorders and cardiometabolic diseases because the inherited mutation levels can change within developing organs and non-dividing cells over time. Initial predictions suggested that the segregation of mtDNA heteroplasmy was largely stochastic, with an equal tendency for levels to increase or decrease. However, transgenic animal work and single-cell analysis have shown this not to be the case during germ-line transmission and in somatic tissues during life. Mutation levels in specific mtDNA regions can increase or decrease in different contexts and the underlying molecular mechanisms are starting to be unraveled. In this review we provide a synthesis of recent literature on the mechanisms of selection for and against mtDNA variants. We identify the most pertinent gaps in our understanding and suggest ways these could be addressed using state of the art techniques.

The human mitochondrial genome contains a second light strand promoter

The human mitochondrial genome must be replicated and expressed in a timely manner to maintain energy metabolism and supply cells with adequate levels of adenosine triphosphate. Central to this process is the idea that replication primers and gene products both arise via transcription from a single light strand promoter (LSP) such that primer formation can influence gene expression, with no consensus as to how this is regulated. Here, we report the discovery of a second light strand promoter (LSP2) in humans, with features characteristic of a bona fide mitochondrial promoter. We propose that the position of LSP2 on the mitochondrial genome allows replication and gene expression to be orchestrated from two distinct sites, which expands our long-held understanding of mitochondrial gene expression in humans.

Mitochondrial DNA variation in Parkinson’s disease: Analysis of “out-of-place” population variants as a risk factor

Mitochondrial DNA (mtDNA), a potential source of mitochondrial dysfunction, has been implicated in Parkinson’s disease (PD). However, many previous studies investigating associations between mtDNA population variation and PD reported inconsistent or contradictory findings. Here, we investigated an alternative hypothesis to determine whether mtDNA variation could play a significant role in PD risk. Emerging evidence suggests that haplogroup-defining mtDNA variants may have pathogenic potential if they occur “out-of-place” on a different maternal lineage. We hypothesized that the mtDNA of PD cases would be enriched for out-of-place variation in genes encoding components of the oxidative phosphorylation complexes. We tested this hypothesis with a unique dataset comprising whole mitochondrial genomes of 70 African ancestry PD cases, two African ancestry control groups (n = 78 and n = 53) and a replication group of 281 European ancestry PD cases and 140 controls from the Parkinson’s Progression Markers Initiative cohort. Significantly more African ancestry PD cases had out-of-place variants than controls from the second control group (P < 0.0125), although this association was not observed in the first control group nor the replication group. As the first mtDNA study to include African ancestry PD cases and to explore out-of-place variation in a PD context, we found evidence that such variation might be significant in this context, thereby warranting further replication in larger cohorts.

Human clinical mutations in mitochondrially encoded subunits of Complex I can be successfully modeled in E. coli

Respiratory Complex I is the site of a large fraction of the mutations that appear to cause mitochondrial disease. Seven of its subunits are mitochondrially encoded, and therefore, such mutants are particularly difficult to construct in cell-culture model systems. We have selected 13 human clinical mutations found in ND2, ND3, ND4, ND4L, ND5 and ND6 that are generally found at subunit interfaces, and not in critical residues. These mutations have been modeled in E. coli subunits of Complex I, nuoN, nuoA, nuoM, nuoK, nuoL, and nuoJ, respectively. All mutants were expressed from a plasmid encoding the entire nuo operon, and membrane vesicles were analyzed for deamino-NADH oxidase activity, and proton translocation activity. ND5 mutants were also analyzed using a time-delayed expression system, recently described by this lab. Other mutants were analyzed for the ability to associate in subcomplexes, after expression of subsets of the genes. For most mutants there was a positive correlation between those that were previously determined to be pathogenic, or likely to be pathogenic, and those that we found with compromised Complex I activity or subunit interactions in E. coli. In conclusion, this approach provides another way to explore the deleterious effects of human mitochondrial mutations, and it can contribute to molecular understanding of such mutations.

Genetic Diversity, Heteroplasmy, and Recombination in Mitochondrial Genomes of Daphnia pulex, Daphnia pulicaria, and Daphnia obtusa

Genetic variants of mitochondrial DNA at the individual (heteroplasmy) and population (polymorphism) levels provide insight into their roles in multiple cellular and evolutionary processes. However, owing to the paucity of genome-wide data at the within-individual and population levels, the broad patterns of these two forms of variation remain poorly understood. Here, we analyze 1,804 complete mitochondrial genome sequences from Daphnia pulex, Daphnia pulicaria, and Daphnia obtusa. Extensive heteroplasmy is observed in D. obtusa, where the high level of intraclonal divergence must have resulted from a biparental-inheritance event, and recombination in the mitochondrial genome is apparent, although perhaps not widespread. Global samples of D. pulex reveal remarkably low mitochondrial effective population sizes, <3% of those for the nuclear genome. In addition, levels of population diversity in mitochondrial and nuclear genomes are uncorrelated across populations, suggesting an idiosyncratic evolutionary history of mitochondria in D. pulex. These population-genetic features appear to be a consequence of background selection associated with highly deleterious mutations arising in the strongly linked mitochondrial genome, which is consistent with polymorphism and divergence data suggesting a predominance of strong purifying selection. Nonetheless, the fixation of mildly deleterious mutations in the mitochondrial genome also appears to be driving positive selection on genes encoded in the nuclear genome whose products are deployed in the mitochondrion.

Heteroplasmic mitochondrial DNA variants in cardiovascular diseases

Mitochondria are implicated in the pathogenesis of cardiovascular diseases (CVDs) but the reasons for this are not well understood. Maternally-inherited population variants of mitochondrial DNA (mtDNA) which affect all mtDNA molecules (homoplasmic) are associated with cardiometabolic traits and the risk of developing cardiovascular disease. However, it is not known whether mtDNA mutations only affecting a proportion of mtDNA molecules (heteroplasmic) also play a role. To address this question, we performed a high-depth (~1000-fold) mtDNA sequencing of blood DNA in 1,399 individuals with hypertension (HTN), 1,946 with ischemic heart disease (IHD), 2,146 with ischemic stroke (IS), and 723 healthy controls. We show that the per individual burden of heteroplasmic single nucleotide variants (mtSNVs) increases with age. The age-effect was stronger for low-level heteroplasmies (heteroplasmic fraction, HF, 5-10%), likely reflecting acquired somatic events based on trinucleotide mutational signatures. After correcting for age and other confounders, intermediate heteroplasmies (HF 10-95%) were more common in hypertension, particularly involving non-synonymous variants altering the amino acid sequence of essential respiratory chain proteins. These findings raise the possibility that heteroplasmic mtSNVs play a role in the pathophysiology of hypertension.

Scinderin promotes fusion of electron transport chain dysfunctional muscle stem cells with myofibers

Muscle stem cells (MuSCs) experience age-associated declines in number and function, accompanied by mitochondrial electron transport chain (ETC) dysfunction and increased reactive oxygen species (ROS). The source of these changes, and how MuSCs respond to mitochondrial dysfunction, is unknown. We report here that in response to mitochondrial ROS, murine MuSCs directly fuse with neighboring myofibers; this phenomenon removes ETC-dysfunctional MuSCs from the stem cell compartment. MuSC-myofiber fusion is dependent on the induction of Scinderin, which promotes formation of actin-dependent protrusions required for membrane fusion. During aging, we find that the declining MuSC population accumulates mutations in the mitochondrial genome, but selects against dysfunctional variants. In the absence of clearance by Scinderin, the decline in MuSC numbers during aging is repressed; however, ETC-dysfunctional MuSCs are retained and can regenerate dysfunctional myofibers. We propose a model in which ETC-dysfunctional MuSCs are removed from the stem cell compartment by fusing with differentiated tissue.

Is osteoarthritis a mitochondrial disease? What is the evidence

PROPOSE OF REVIEW

To summarize the evidence that suggests that osteoarthritis (OA) is a mitochondrial disease.

RECENT FINDINGS

Mitochondrial dysfunction together with mtDNA damage could contribute to cartilage degradation via several processes such as: (1) increased apoptosis; (2) decreased autophagy; (3) enhanced inflammatory response; (4) telomere shortening and increased senescence chondrocytes; (5) decreased mitochondrial biogenesis and mitophagy; (6) increased cartilage catabolism; (7) increased mitochondrial fusion leading to further reactive oxygen species production; and (8) impaired metabolic flexibility.

SUMMARY

Mitochondria play an important role in some events involved in the pathogenesis of OA, such as energy production, the generation of reactive oxygen and nitrogen species, apoptosis, authophagy, senescence and inflammation. The regulation of these processes in the cartilage is at least partially controlled by retrograde regulation from mitochondria and mitochondrial genetic variation. Retrograde regulation through mitochondrial haplogroups exerts a signaling control over the nuclear epigenome, which leads to the modulation of nuclear genes, cellular functions and development of OA. All these data suggest that OA could be considered a mitochondrial disease as well as other complex chronic disease as cancer, cardiovascular and neurologic diseases.

Genetic landscape of human mitochondrial genome using whole genome sequencing

Increasing evidences suggest that mitochondrial dysfunction is implicated in diseases and aging, and whole-genome sequencing (WGS) is the most unbiased method in analyzing the mitochondrial genome (mtDNA). However, the genetic landscape of mtDNA in the Chinese population has not been fully examined. Here, we described the genetic landscape of mtDNA using WGS data from Chinese individuals (n = 3241). We identified 3892 mtDNA variants, of which 3349 (86%) were rare variants. Interestingly, we observed a trend toward extreme heterogeneity of mtDNA variants. Our study observed a distinct purifying selection on mtDNA, which inhibits the accumulation of harmful heteroplasmies at the individual level: (1) mitochondrial dN/dS ratios were much less than 1; (2) the dN/dS ratio of heteroplasmies was higher than homoplasmies; (3) heteroplasmies had more indels and predicted deleterious variants than homoplasmies. Furthermore, we found that haplogroup M (20.27%) and D (20.15%) had the highest frequencies in the Chinese population, followed by B (18.51%) and F (16.45%). The number of variants per individual differed across haplogroup groups, with a higher number of homoplasmies for the M lineage. Meanwhile, mtDNA copy number was negatively correlated with age but positively correlated with the female sex. Finally, we developed an mtDNA variation database of Chinese populations called MTCards (http://genemed.tech/mtcards/) to facilitate the query of mtDNA variants in this study. In summary, these findings contribute to different aspects of understanding mtDNA, providing a better understanding of the genetic basis of mitochondrial-related diseases.

MitoPhen database: a human phenotype ontology-based approach to identify mitochondrial DNA diseases

Diagnosing mitochondrial disorders remains challenging. This is partly because the clinical phenotypes of patients overlap with those of other sporadic and inherited disorders. Although the widespread availability of genetic testing has increased the rate of diagnosis, the combination of phenotypic and genetic heterogeneity still makes it difficult to reach a timely molecular diagnosis with confidence. An objective, systematic method for describing the phenotypic spectra for each variant provides a potential solution to this problem. We curated the clinical phenotypes of 6688 published individuals with 89 pathogenic mitochondrial DNA (mtDNA) mutations, collating 26 348 human phenotype ontology (HPO) terms to establish the MitoPhen database. This enabled a hypothesis-free definition of mtDNA clinical syndromes, an overview of heteroplasmy-phenotype relationships, the identification of under-recognized phenotypes, and provides a publicly available reference dataset for objective clinical comparison with new patients using the HPO. Studying 77 patients with independently confirmed positive mtDNA diagnoses and 1083 confirmed rare disease cases with a non-mitochondrial nuclear genetic diagnosis, we show that HPO-based phenotype similarity scores can distinguish these two classes of rare disease patients with a false discovery rate <10% at a sensitivity of 80%. Enriching the MitoPhen database with more patients will improve predictions for increasingly rare variants.

Cell reprogramming shapes the mitochondrial DNA landscape

Individual induced pluripotent stem cells (iPSCs) show considerable phenotypic heterogeneity, but the reasons for this are not fully understood. Comprehensively analysing the mitochondrial genome (mtDNA) in 146 iPSC and fibroblast lines from 151 donors, we show that most age-related fibroblast mtDNA mutations are lost during reprogramming. However, iPSC-specific mutations are seen in 76.6% (108/141) of iPSC lines at a mutation rate of 8.62 × 10-5/base pair. The mutations observed in iPSC lines affect a higher proportion of mtDNA molecules, favouring non-synonymous protein-coding and tRNA variants, including known disease-causing mutations. Analysing 11,538 single cells shows stable heteroplasmy in sub-clones derived from the original donor during differentiation, with mtDNA variants influencing the expression of key genes involved in mitochondrial metabolism and epidermal cell differentiation. Thus, the dynamic mtDNA landscape contributes to the heterogeneity of human iPSCs and should be considered when using reprogrammed cells experimentally or as a therapy.

Clinical Bioinformatics in Precise Diagnosis of Mitochondrial Disease

Clinical bioinformatics system is well-established for diagnosing genetic disease based on next-generation sequencing, but requires special considerations when being adapted for the next-generation sequencing-based genetic diagnosis of mitochondrial diseases. Challenges are caused by the involvement of mitochondrial DNA genome in disease etiology. Heteroplasmy and haplogroup are key factors in interpreting mitochondrial DNA variant effects. Data resources and tools for analyzing variant and sequencing data are available at MSeqDR, MitoMap, and HmtDB. Revised specifications of the American College of Medical Genetics/Association of Molecular Pathology standards and guidelines for mitochondrial DNA variant interpretation are proposed by the MSeqDr Consortium and community experts.

Extreme heterogeneity of human mitochondrial DNA from organelles to populations

Contrary to the long-held view that most humans harbour only identical mitochondrial genomes, deep resequencing has uncovered unanticipated extreme genetic variation within mitochondrial DNA (mtDNA). Most, if not all, humans contain multiple mtDNA genotypes (heteroplasmy); specific patterns of variants accumulate in different tissues, including cancers, over time; and some variants are preferentially passed down or suppressed in the maternal germ line. These findings cast light on the origin and spread of mtDNA mutations at multiple scales, from the organelle to the human population, and challenge the conventional view that high percentages of a mutation are required before a new variant has functional consequences.

Mitochondrial genome variation in male LHON patients with the m.11778G > A mutation

Leber hereditary optic neuropathy (LHON) is a mitochondrial disorder with symptoms limited to a single tissue, optic nerve, resulting in vision loss. In the majority of cases it is caused by one of three point mutations in mitochondrial DNA (mtDNA) but their presence is not sufficient for disease development, since ~50% of men and ~10% women who carry them are affected. Thus additional modifying factors must exist. In this study, we use next generation sequencing to investigate the role of whole mtDNA variation in male Polish patients with LHON and m.11778G > A, the most frequent LHON mutation. We present a possible association between mtDNA haplogroup K and variants in its background, a combination of m.3480A > G, m.9055G > A, m.11299 T > C and m.14167C > T, and LHON mutation. These variants may have a negative effect on m.11778G > A increasing its penetrance and the risk of LHON in the Polish population. Surprisingly, we did not observe associations previously reported for m.11778G > A and LHON in European populations, particularly for haplogroup J as a risk factor, implying that mtDNA variation is much more complex. Our results indicate possible contribution of novel combination of mtDNA genetic factors to the LHON phenotype.

Specifications of the ACMG/AMP standards and guidelines for mitochondrial DNA variant interpretation

Mitochondrial DNA (mtDNA) variant pathogenicity interpretation has special considerations given unique features of the mtDNA genome, including maternal inheritance, variant heteroplasmy, threshold effect, absence of splicing, and contextual effects of haplogroups. Currently, there are insufficient standardized criteria for mtDNA variant assessment, which leads to inconsistencies in clinical variant pathogenicity reporting. An international working group of mtDNA experts was assembled within the Mitochondrial Disease Sequence Data Resource Consortium and obtained Expert Panel status from ClinGen. This group reviewed the 2015 American College of Medical Genetics and Association of Molecular Pathology standards and guidelines that are widely used for clinical interpretation of DNA sequence variants and provided further specifications for additional and specific guidance related to mtDNA variant classification. These Expert Panel consensus specifications allow for consistent consideration of the unique aspects of the mtDNA genome that directly influence variant assessment, including addressing mtDNA genome composition and structure, haplogroups and phylogeny, maternal inheritance, heteroplasmy, and functional analyses unique to mtDNA, as well as specifications for utilization of mtDNA genomic databases and computational algorithms.

Ancient mitochondrial DNA pathogenic variants putatively associated with mitochondrial disease

Mitochondrial DNA variants associated with diseases are widely studied in contemporary populations, but their prevalence has not yet been investigated in ancient populations. The publicly available AmtDB database contains 1443 ancient mtDNA Eurasian genomes from different periods. The objective of this study was to use this data to establish the presence of pathogenic mtDNA variants putatively associated with mitochondrial diseases in ancient populations. The clinical significance, pathogenicity prediction and contemporary frequency of mtDNA variants were determined using online platforms. The analyzed ancient mtDNAs contain six variants designated as being "confirmed pathogenic" in modern patients. The oldest of these, m.7510T>C in the MT-TS1 gene, was found in a sample from the Neolithic period, dated 5800-5400 BCE. All six have well established clinical association, and their pathogenic effect is corroborated by very low population frequencies in contemporary populations. Analysis of the geographic location of the ancient samples, contemporary epidemiological trends and probable haplogroup association indicate diverse spatiotemporal dynamics of these variants. The dynamics in the prevalence and distribution is conceivably result of de novo mutations or human migrations and subsequent evolutionary processes. In addition, ten variants designated as possibly or likely pathogenic were found, but the clinical effect of these is not yet well established and further research is warranted. All detected mutations putatively associated with mitochondrial disease in ancient mtDNA samples are in tRNA coding genes. Most of these mutations are in a mt-tRNA type (Model 2) that is characterized by loss of D-loop/T-loop interaction. Exposing pathogenic variants in ancient human populations expands our understanding of their origin and prevalence dynamics.

Whole Mitochondrial Genome Analysis in Serbian Cases of Leber's Hereditary Optic Neuropathy

Leber’s hereditary optic neuropathy (LHON) is a maternally inherited disorder that affects central vision in young adults and is typically associated with mitochondrial DNA (mtDNA) mutations. This study is based on a mutational screening of entire mtDNA in eight Serbian probands clinically and genetically diagnosed with LHON and four of their family members, who are asymptomatic mutation carriers. All obtained sequence variants were compared to human mtDNA databases, and their potential pathogenic characteristics were assessed by bioinformatics tools. Mitochondrial haplogroup analysis was performed by MITOMASTER. Our study revealed two well-known primary LHON mutations, m.11778G>A and m.3460G>A, and one rare LHON mutation, m.8836A>G. Various secondary mutations were detected in association with the primary mutations. MITOMASTER analysis showed that the two well-known primary mutations belong to the R haplogroup, while the rare LHON m.8836A>G was detected within the N1b haplogroup. Our results support the need for further studies of genetic background and its role in the penetrance and severity of LHON.

Mitochondrial Diseases: A Diagnostic Revolution

Mitochondrial disorders have emerged as a common cause of inherited disease, but are traditionally viewed as being difficult to diagnose clinically, and even more difficult to comprehensively characterize at the molecular level. However, new sequencing approaches, particularly whole-genome sequencing (WGS), have dramatically changed the landscape. The combined analysis of nuclear and mitochondrial DNA (mtDNA) allows rapid diagnosis for the vast majority of patients, but new challenges have emerged. We review recent discoveries that will benefit patients and families, and highlight emerging questions that remain to be resolved.

What can a comparative genomics approach tell us about the pathogenicity of mtDNA mutations in human populations?

Mitochondrial disorders are heterogeneous, showing variable presentation and penetrance. Over the last three decades, our ability to recognize mitochondrial patients and diagnose these mutations, linking genotype to phenotype, has greatly improved. However, it has become increasingly clear that these strides in diagnostics have not benefited all population groups. Recent studies have demonstrated that patients from genetically understudied populations, in particular those of black African heritage, are less likely to receive a diagnosis of mtDNA disease. It has been suggested that haplogroup context might influence the presentation and penetrance of mtDNA disease; thus, the spectrum of mutations that are associated with disease in different populations. However, to date there is only one well-established example of such an effect: the increased penetrance of two Leber's hereditary optic neuropathy mutations on a haplogroup J background. This paper conducted the most extensive investigation to date into the importance of haplogroup context on the pathogenicity of mtDNA mutations. We searched for proven human point mutations across 726 multiple sequence alignments derived from 33 non-human species absent of disease. A total of 58 pathogenic point mutations arise in the sequences of these species. We assessed the sequence context and found evidence of population variants that could modulate the phenotypic expression of these point mutations masking the pathogenic effects seen in humans. This supports the theory that sequence context is influential in the presentation of mtDNA disease and has implications for diagnostic practices. We have shown that our current understanding of the pathogenicity of mtDNA point mutations, primarily built on studies of individuals with haplogroups HVUKTJ, will not present a complete picture. This will have the effect of creating a diagnostic inequality, whereby individuals who do not belong to these lineages are less likely to receive a genetic diagnosis.

Rectifying long-standing misconceptions about the ρ statistic for molecular dating

When divided by a given mutation rate, the ρ (rho) statistic provides a simple estimator of the age of a clade within a phylogenetic tree by averaging the number of mutations from each sample in the clade to its root. However, a long-standing critique of the use of ρ in genetic dating has been quite often cited. Here we show that the critique is unfounded. We demonstrate by a formal mathematical argument and illustrate with a simulation study that ρ estimates are unbiased and also that ρ and maximum likelihood estimates do not differ in any systematic fashion. We also demonstrate that the claim that the associated confidence intervals commonly estimate the uncertainty inappropriately is flawed since it relies on a means of calculating standard errors that is not used by any other researchers, whereas an established expression for the standard error is largely unproblematic. We conclude that ρ dating, alongside approaches such as maximum likelihood (ML) and Bayesian inference, remains a useful tool for genetic dating.

New mtDNA Association Model, MutPred Variant Load, Suggests Individuals With Multiple Mildly Deleterious mtDNA Variants Are More Likely to Suffer From Atherosclerosis

The etiology of common complex diseases is multifactorial, involving both genetic, and environmental factors. A role for mitochondrial dysfunction and mitochondrial DNA (mtDNA) variation has been suggested in the pathogenesis of common complex traits. The aim of this study was to investigate a potential role of mtDNA variants in the development of obesity, diabetes, and atherosclerosis in the Polish population. Whole mtDNA sequences from 415 Polish individuals representing three disease cohorts and a control group were obtained using high-throughput sequencing. Two approaches for the assessment of mtDNA variation were applied, traditional mitochondrial haplogroup association analysis and the mutational or variant load model using the MutPred pathogenicity prediction algorithm for amino acid substitutions in humans. We present a possible association between mildly deleterious mtDNA variant load and atherosclerosis that might be due to having more than one likely mildly deleterious non-synonymous substitution. Moreover, it seems largely dependent upon a few common haplogroup associated variants with MutPred score above 0.5.

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.

Whole sequence of the mitochondrial DNA genome of Kearns Sayre Syndrome patients: Identification of deletions and variants

Mitochondria both produce the energy of the cell as ATP via respiration and regulate cellular metabolism. Accordingly, any deletion or mutation in the mitochondrial DNA (mtDNA) may result in a disease. One of these diseases is Kearns Sayre syndrome (KSS), described for the first time in 1958, where different large-scale deletions of different sizes and at different positions have been reported in the mitochondrial genome of patients with similar clinical symptoms. In this study, sequences of the mitochondrial genome of three patients with clinic features of KSS were analyzed. Our results revealed the position, heteroplasmy percentage, size of deletions, and their haplogroups. Two patients contained deletions reported previously and one patient showed a new deletion not reported previously. These results display for the first time a systematic analysis of mtDNA variants in the whole mtDNA genome of patients with KSS to help to understand their association with the disease.

Oldies but Goldies mtDNA Population Variants and Neurodegenerative Diseases

mtDNA is transmitted through the maternal line and its sequence variability, which is population specific, is assumed to be phenotypically neutral. However, several studies have shown associations between the variants defining some genetic backgrounds and the susceptibility to several pathogenic phenotypes, including neurodegenerative diseases. Many of these studies have found that some of these variants impact many of these phenotypes, including the ones defining the Caucasian haplogroups H, J, and Uk, while others, such as the ones defining the T haplogroup, have phenotype specific associations. In this review, we will focus on those that have shown a pleiotropic effect in population studies in neurological diseases. We will also explore their bioenergetic and genomic characteristics in order to provide an insight into the role of these variants in disease. Given the importance of mitochondrial population variants in neurodegenerative diseases a deeper analysis of their effects might unravel new mechanisms of disease and help design new strategies for successful treatments.

Mitochondrial Genomics: A complex field now coming of age

Purpose of review

The groundwork for mitochondrial medicine was laid 30 years ago with identification of the first disease-causing mitochondrial DNA (mtDNA) mutations in 1988. Three decades later, mutations in nearly 300 genes involving every possible mode of inheritance within both nuclear and mitochondrial genomes are now recognized to collectively comprise the largest class of inherited metabolic disease affecting at least 1 in 4,300 individuals across all ages. Significant progress has been made in recent years to improve understanding of mitochondrial biology and disease pathophysiology.

Recent findings

Markedly improved understanding of the highly diverse molecular etiologies of multi-systemic phenotypes in primary mitochondrial disease has resulted from massively parallel genomic sequencing technologies and improved bioinformatic resources that enable identification in individual patients of their disease's precise genetic etiology. Key informatics resources of particular utility to the mitochondrial disease genomics community have been developed, including: (1) Mitocarta 2.0 repository of 1200+ verified mitochondria-localized proteins, (2) MITOMAP Web resource of curated mtDNA genome variants, and (3) Mitochondrial Disease Sequence Data Resource (MSeqDR) that centralizes Web curation and annotation of mitochondrial disease genes and variants in both genomes, ontology-defined phenotypes, and access to many analytic tools to support genomic data mining and interpretation. Gene and mutation-based disease categorization has proven particularly useful to identify the full clinical spectrum of disease that may affect a given individual.

Summary

Extensive genomic advances, both in technologic platforms and bioinformatics resources, have facilitated dramatic improvement in the accurate recognition and understanding of primary mitochondrial disease.