Mitochondrial DNA point mutations and large-scale deletions in maternally inherited metabolic disorders: Genotype-phenotype correlations

Introduction

The human mitochondrial genome is a 16,569 base-pair circular molecule encoding 13 polypeptides essential for oxidative phosphorylation, alongside 22 transfer RNAs and 2 ribosomal RNAs [1] . Since its complete sequencing in 1981, over 300 pathogenic mtDNA variants have been catalogued, producing an extraordinary range of clinical manifestations that can affect virtually any organ system [2] . The maternal inheritance pattern and polyplasmic nature of mtDNA create a unique genetic landscape where mutant and wild-type molecules coexist within individual cells, a state known as heteroplasmy [3] . Clinical expression of mtDNA mutations typically follows a threshold effect: symptoms manifest when the proportion of mutant molecules exceeds a tissue-specific critical level, generally between 60% and 90% [4] . But this threshold isn't fixed. It varies with mutation type, affected gene, and the metabolic demands of the tissue in question. Highly oxidative tissues such as brain, heart, and skeletal muscle tend to be affected at lower heteroplasmy levels than metabolically less demanding organs [5] . Despite extensive case reports and small series, systematic analyses correlating precise heteroplasmy measurements across multiple tissues with standardized clinical severity assessments remain limited [6] . This research examined 128 patients with confirmed pathogenic mtDNA variants to establish quantitative genotype-phenotype relationships across the most common mutation categories and clinical phenotypes encountered in mitochondrial medicine.

Material and Methods

Material

Patient recruitment occurred at the Mitochondrial Disease Clinic of Wellington Institute of Biomedical Sciences from February 2022 to September 2023. The cohort comprised 128 individuals (74 female, 54 male; age range 4-71 years, mean 38.2±16.7 years) with Genetically confirmed pathogenic mtDNA variants. Inclusion criteria required molecular confirmation of a recognized pathogenic mtDNA mutation or deletion and availability of at Least two tissue samples for heteroplasmy quantification. Peripheral blood, skeletal muscle biopsy, and urinary epithelial cell samples were collected from each participant following institutional ethics committee approval (Protocol WIB-2022-047, approved January 2022).

Extraction protocol: Total DNA was extracted from blood using the QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) and from muscle tissue using the DNeasy Blood and Tissue Kit following manufacturer protocols. Urinary epithelial cells were pelleted by centrifugation at 3,000 × g for ten minutes and processed using the same tissue extraction kit. DNA concentration and purity were assessed by NanoDrop spectrophotometry (Thermo Fisher Scientific, Waltham, MA, USA), requiring A260/A280 ratios between 1.8 and 2.0 for inclusion [7] .

Methods: Point mutations were detected and quantified using fluorescent restriction fragment length polymorphism analysis and confirmed by Sanger sequencing [8] . Largescale deletions were identified through long-range PCR and characterized by breakpoint sequencing. Heteroplasmy quantification employed last-cycle fluorescent PCR with capillary electrophoresis, providing resolution of ±2% [9] .

Clinical severity was assessed using a standardized Newcastle Mitochondrial Disease Scale, which assigns composite scores from 1 (mild) to 10 (severe) based on neurological, cardiac, endocrine, ophthalmological, and muscular involvement [10] . Statistical analysis used Spearman's rank correlation for heteroplasmy-severity associations and Kruskal-Wallis test with Dunn's post-hoc correction for between-group comparisons.

Quality control

All heteroplasmy measurements were performed in duplicate, and samples with inter-assay variation exceeding 5% were re-analyzed. Positive controls with known heteroplasmy levels (constructed by mixing cloned wildtype and mutant mtDNA at defined ratios) were included in each analytical run. External quality assessment was performed through participation in the European Molecular Genetics Quality Network mitochondrial proficiency testing scheme [11] .

Results

The 128 patients carried five major mutation categories: m.3243A>G (n=38), m.8344A>G (n=24), m.8993T>G (n=19), single large-scale deletions (n=31), and multiple deletions (n=16). Table 1 presents the distribution of mutations across clinical phenotypes with associated heteroplasmy data.

Comprehensive interpretation

The overall correlation between blood heteroplasmy and clinical severity (r² = 0.73) confirms the expected doseresponse relationship, though considerable individual variation exists within each mutation group. The m.8993T>G mutation showed the narrowest heteroplasmy range but highest median severity, reflecting the sharp biochemical threshold associated with ATPase 6 dysfunction. Patients with the m.3243A>G mutation presented the broadest phenotypic spectrum, ranging from isolated diabetes with deafness at low heteroplasmy to full MELAS syndrome at high levels. Figure 3 maps the locations of the principal disease-associated mutations on the mitochondrial genome.

Discussion

These findings confirm that heteroplasmy level serves as a useful but imperfect predictor of clinical severity in mtDNA disorders. The 73% variance explained by heteroplasmy alone leaves a substantial portion attributable to nuclear genetic background, environmental modifiers, and stochastic factors influencing mitochondrial segregation during development [12] . The tissue-specific differences in heteroplasmy observed here, with muscle consistently exceeding blood by 20-30 percentage points for tRNA mutations, have practical implications for diagnostic testing. Blood-based screening may underestimate disease burden, and muscle biopsy remains the gold standard for accurate heteroplasmy assessment.

The clear phenotypic boundaries associated with the m.8993T>G mutation contrast sharply with the continuous spectrum seen for m.3243A>G, suggesting fundamental differences in how these mutations affect cellular bioenergetics. The m.8993T>G variant directly impairs ATP synthase function, creating a binary-like effect above its threshold. In contrast, the m.3243A>G tRNA-Leu mutation affects translation of all 13 mtDNA-encoded proteins to varying degrees, producing a more graded cellular phenotype that may intersect with tissue-specific compensatory mechanisms [13] .

Conclusions

This research established quantitative genotype-phenotype correlations across 128 patients with five categories of pathogenic mtDNA variants. Heteroplasmy level in blood correlated positively with clinical severity but varied substantially in its predictive value across mutation types. The m.8993T>G mutation showed the tightest genotypephenotype coupling, while m.3243A>G and single deletion patients displayed the greatest clinical variability. Tissuespecific heteroplasmy differences underscore the importance of analyzing multiple sample types for accurate diagnostic assessment. These correlations can inform genetic counseling, prognostic estimation, and clinical trial stratification in mitochondrial medicine. Future efforts should integrate nuclear genome variation and metabolomic profiling to improve individual outcome prediction.

Acknowledgement

Acknowledgements Funding SourcesThis research was funded by the Health Research Council of New Zealand through the Rare Disease Research Program. Institutional supportThe authors thank Wellington Institute of Biomedical Sciences for providing access to the molecular diagnostics laboratory and clinical facilities. Contributions not qualifying for authorshipWe acknowledge the genetic counselors at the Mitochondrial Disease Clinic for patient coordination and Dr. Rachel Ngata for histological analysis of muscle biopsies.

References

1. Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson AR, Drouin J, et al. Sequence and organization of the human mitochondrial genome. Nature. 1981;290(5806):457-465.
2. Lott MT, Leipzig JN, Derbeneva O, Xie HM, Chalkia D, Sarmady M, et al. mtDNA variation and analysis using MITOMAP and MITOMASTER. Curr Protoc Bioinformatics. 2013;44:1.23.1-1.23.26.
3. Stewart JB, Chinnery PF. The dynamics of mitochondrial DNA heteroplasmy: implications for human health and disease. Nat Rev Genet. 2015;16(9):530-542.
4. Rossignol R, Faustin B, Rocher C, Malgat M, Mazat JP, Letellier T. Mitochondrial threshold effects. Biochem J. 2003;370(Pt 3):751-762.
5. Schon EA, DiMauro S, Hirano M. Human mitochondrial DNA: roles of inherited and somatic mutations. Nat Rev Genet. 2012;13(12):878-890.
6. Gorman GS, Schaefer AM, Ng Y, Gomez N, Blakely EL, Alston CL, et al. Prevalence of nuclear and mitochondrial DNA mutations related to adult mitochondrial disease. Ann Neurol. 2015;77(5):753759.
7. Dimmock D, Tang LY, Schmitt ES, Wong LJ. Quantitative evaluation of the mitochondrial DNA depletion syndrome. Clin Chem. 2010;56(7):11191127.
8. White HE, Durston VJ, Seller A, Fratter C, Harvey JF, Cross NC. Accurate detection and quantitation of heteroplasmic mitochondrial point mutations by pyrosequencing. Genet Test. 2005;9(3):190-199.
9. Bai RK, Wong LJ. Detection and quantification of heteroplasmic mutant mitochondrial DNA by real-time amplification refractory mutation system quantitative PCR analysis. Methods Mol Biol. 2004;197:157-169.
10. Schaefer AM, Phoenix C, Elson JL, McFarland R, Chinnery PF, Turnbull DM. Mitochondrial disease in adults: a scale to monitor progression and treatment. Neurology. 2006;66(12):1932-1934.
11. Blackwood JK, Whittaker RG, Blakely EL, Alston CL, Turnbull DM, Taylor RW. Investigation and diagnosis of pathogenic mitochondrial DNA mutations in human urothelial cells. Biochem Soc Trans. 2010;38(6):16331637.
12. Chinnery PF, Howell N, Lightowlers RN, Turnbull DM. Molecular pathology of MELAS and MERRF: relationship between mutation load and clinical phenotypes. Brain. 1997;120(Pt 10):1713-1721.
13. Tatuch Y, Christodoulou J, Feigenbaum A, Clarke JT, Wherret J, Smith C, et al. Heteroplasmic mtDNA mutation (T→G) at position 8993 causes Leigh disease when the proportion of mutant mtDNA is high. Am J Hum Genet. 1992;50(4):852-858.
14. McFarland R, Taylor RW, Turnbull DM. A neurological perspective on mitochondrial disease. Lancet Neurol. 2010;9(8):829-840.