Primary mitochondrial disorders and mimics: Insights from a large French cohort.

OBJECTIVE: The objective of this study was to evaluate the implementation of NGS within the French mitochondrial network, MitoDiag, from targeted gene panels to whole exome sequencing (WES) or whole genome sequencing (WGS) focusing on mitochondrial nuclear-encoded genes. METHODS: Over 2000 patients suspected of Primary Mitochondrial Diseases (PMD) were sequenced by either targeted gene panels, WES or WGS within MitoDiag. We described the clinical, biochemical, and molecular data of 397 genetically confirmed patients, comprising 294 children and 103 adults, carrying pathogenic or likely pathogenic variants in nuclear-encoded genes. RESULTS: The cohort exhibited a large genetic heterogeneity, with the identification of 172 distinct genes and 253 novel variants. Among children, a notable prevalence of pathogenic variants in genes associated with oxidative phosphorylation (OXPHOS) functions and mitochondrial translation was observed. In adults, pathogenic variants were primarily identified in genes linked to mtDNA maintenance. Additionally, a substantial proportion of patients (54% (42/78) and 48% (13/27) in children and adults, respectively), undergoing WES or WGS testing displayed PMD mimics, representing pathologies that clinically resemble mitochondrial diseases. INTERPRETATION: We reported the largest French cohort of patients suspected of PMD with pathogenic variants in nuclear genes. We have emphasized the clinical complexity of PMD and the challenges associated with recognizing and distinguishing them from other pathologies, particularly neuromuscular disorders. We confirmed that WES/WGS, instead of panel approach, was more valuable to identify the genetic basis in patients with "possible" PMD and we provided a genetic testing flowchart to guide physicians in their diagnostic strategy.

CHCHD10S59L/+ mouse model: Behavioral and neuropathological features of frontotemporal dementia.

CHCHD10-related disease causes a spectrum of clinical presentations including mitochondrial myopathy, cardiomyopathy, amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). We generated a knock-in mouse model bearing the p.Ser59Leu (S59L) CHCHD10 variant. Chchd10S59L/+ mice have been shown to phenotypically replicate the disorders observed in patients: myopathy with mtDNA instability, cardiomyopathy and typical ALS features (protein aggregation, neuromuscular junction degeneration and spinal motor neuron loss). Here, we conducted a comprehensive behavioral, electrophysiological and neuropathological assessment of Chchd10S59L/+ mice. These animals show impaired learning and memory capacities with reduced long-term potentiation (LTP) measured at the Perforant Pathway-Dentate Gyrus (PP-DG) synapses. In the hippocampus of Chchd10S59L/+ mice, neuropathological studies show the involvement of protein aggregates, activation of the integrated stress response (ISR) and neuroinflammation in the degenerative process. These findings contribute to decipher mechanisms associated with CHCHD10 variants linking mitochondrial dysfunction and neuronal death. They also validate the Chchd10S59L/+ mice as a relevant model for FTD, which can be used for preclinical studies to test new therapeutic strategies for this devastating disease.

A Case Report of SYNE1 Deficiency-Mimicking Mitochondrial Disease and the Value of Pangenomic Investigations.

Mitochondrial disorders are characterized by a huge clinical, biochemical, and genetic heterogeneity, which poses significant diagnostic challenges. Several studies report that more than 50% of patients with suspected mitochondrial disease could have a non-mitochondrial disorder. Thus, only the identification of the causative pathogenic variant can confirm the diagnosis. Herein, we describe the diagnostic journey of a family suspected of having a mitochondrial disorder who were referred to our Genetics Department. The proband presented with the association of cerebellar ataxia, COX-negative fibers on muscle histology, and mtDNA deletions. Whole exome sequencing (WES), supplemented by a high-resolution array, comparative genomic hybridization (array-CGH), allowed us to identify two pathogenic variants in the non-mitochondrial SYNE1 gene. The proband and her affected sister were found to be compound heterozygous for a known nonsense variant (c.13258C>T, p.(Arg4420Ter)), and a large intragenic deletion that was predicted to result in a loss of function. To our knowledge, this is the first report of a large intragenic deletion of SYNE1 in patients with cerebellar ataxia (ARCA1). This report highlights the interest in a pangenomic approach to identify the genetic basis in heterogeneous neuromuscular patients with the possible cause of mitochondrial disease. Moreover, even rare copy number variations should be considered in patients with a phenotype suggestive of SYNE1 deficiency.

Mitochondria, a Key Target in Amyotrophic Lateral Sclerosis Pathogenesis.

Mitochondrial dysfunction occurs in numerous neurodegenerative diseases, particularly amyotrophic lateral sclerosis (ALS), where it contributes to motor neuron (MN) death. Of all the factors involved in ALS, mitochondria have been considered as a major player, as secondary mitochondrial dysfunction has been found in various models and patients. Abnormal mitochondrial morphology, defects in mitochondrial dynamics, altered activities of respiratory chain enzymes and increased production of reactive oxygen species have been described. Moreover, the identification of CHCHD10 variants in ALS patients was the first genetic evidence that a mitochondrial defect may be a primary cause of MN damage and directly links mitochondrial dysfunction to the pathogenesis of ALS. In this review, we focus on the role of mitochondria in ALS and highlight the pathogenic variants of ALS genes associated with impaired mitochondrial functions. The multiple pathways demonstrated in ALS pathogenesis suggest that all converge to a common endpoint leading to MN loss. This may explain the disappointing results obtained with treatments targeting a single pathological process. Fighting against mitochondrial dysfunction appears to be a promising avenue for developing combined therapies in the future.

Clinical and genetic analysis further delineates the phenotypic spectrum of ALDH1A3-related anophthalmia and microphthalmia.

Biallelic pathogenic variants in ALDH1A3 are responsible for approximately 11% of recessively inherited cases of severe developmental eye anomalies. Some individuals can display variable neurodevelopmental features, but the relationship to the ALDH1A3 variants remains unclear. Here, we describe seven unrelated families with biallelic pathogenic ALDH1A3 variants: four compound heterozygous and three homozygous. All affected individuals had bilateral anophthalmia/microphthalmia (A/M), three with additional intellectual or developmental delay, one with autism and seizures and three with facial dysmorphic features. This study confirms that individuals with biallelic pathogenic ALDH1A3 variants consistently manifest A/M, but additionally display neurodevelopmental features with significant intra- and interfamilial variability. Furthermore, we describe the first case with cataract and highlight the importance of screening ALDH1A3 variants in nonconsanguineous families with A/M.

A novel WFS1 variant associated with severe diabetic retinopathy in Wolfram syndrome type 1.

BACKGROUND: Wolfram syndrome type 1 is a rare neurodegenerative disorder including diabetes insipidus, diabetes mellitus, optic atrophy, and deafness, with variable additional findings. The phenotypic spectrum is very heterogeneous, with non-autoimmune juvenile-onset diabetes and optic atrophy as minimal criteria for the diagnosis. Biallelic mutations in the WFS1 gene are the causative genetic anomaly for the syndrome, with, however, no evident genotype-phenotype correlation. Among the clinical features of the disease, diabetic retinopathy depicts a rarely reported microvascular complication. In this report, we describe the clinical and genetic findings in a 26-year-old patient presenting with Wolfram syndrome and severe diabetic retinopathy. METHODS: The mutation screening was performed by polymerase chain reaction followed by Sanger sequencing of the entire coding sequence of the WFS1 gene. RESULTS: A novel homozygous missense variant c.1901A>T (p.Lys634Met) was found in the proband and classified as probably pathogenic according to the American College of Medical Genetics and Genomics. CONCLUSIONS: The molecular study of the WFS1 gene is essential for the diagnostic confirmation, to provide appropriate genetic counseling and a mutational screening in the at-risk relatives. The c.1901A>T (p.Lys634 Met) is a novel variant that could be responsible for a severe form of Wolfram syndrome with early and proliferative diabetic retinopathy.

Splicing variants in NARS2 are associated with milder phenotypes and intra-familial variability.

Biallelic rare variants in NARS2 that encode the mitochondrial asparaginyl-tRNA synthetase are associated with a wide spectrum of clinical phenotypes ranging from severe neurodegenerative disorders to isolated mitochondrial myopathy or deafness. To date, only a small number of patients with NARS2 variants have been reported, and possible genotype-phenotype correlations are still lacking. Here, we present three siblings who had an early-onset hearing loss, while one developed severe symptoms in adulthood associated with early intellectual impairment, refractory seizures, moderate axonal sensorimotor neuropathy, and atypical psychiatric symptoms. Biochemical analysis revealed impairment of the activity and assembly of the respiratory chain complexes in this patient's muscle and fibroblasts. Whole Exome Sequencing allowed identification of a heterozygous variant NM_024678.5(NARS2):c.822G > C (p.Gln274His) that is known to be pathogenic and to affect splicing of the NARS2 gene, but was unable to detect a second variant in this gene. Coverage analysis and Sanger sequencing led to identification of a novel intronic deletion NM_024678.5(NARS2):c.922-21_922-19del in the three siblings in trans with the c.822G > C. Functional analysis by RT-PCR showed that this deletion was causing aberrant splicing and led to exon 9 skipping in NARS2 mRNA in patient fibroblasts. Our work expands the phenotype and genotype spectrum of NARS2-related disorders. We provide evidence of the pathogenic effect of a novel intronic deletion in the NARS2 gene and report on additional adult patients with a large intrafamilial variability associated with splice variants in this gene. More specifically, we detail the phenotype of the oldest living patient to date with NARS2 variants and, for the first time, we report the psychiatric symptoms associated with this gene. Our work confirms the complexity of genotype-phenotype correlation in patients with pathogenic NARS2 variants.

Molecular and clinical descriptions of patients with GABAA receptor gene variants (GABRA1, GABRB2, GABRB3, GABRG2): A cohort study, review of literature, and genotype-phenotype correlation.

OBJECTIVE: γ-Aminobutyric acid (GABA)A -receptor subunit variants have recently been associated with neurodevelopmental disorders and/or epilepsy. The phenotype linked with each gene is becoming better known. Because of the common molecular structure and physiological role of these phenotypes, it seemed interesting to describe a putative phenotype associated with GABAA -receptor-related disorders as a whole and seek possible genotype-phenotype correlations. METHODS: We collected clinical, electrophysiological, therapeutic, and molecular data from patients with GABAA -receptor subunit variants (GABRA1, GABRB2, GABRB3, and GABRG2) through a national French collaboration using the EPIGENE network and compared these data to the one already described in the literature. RESULTS: We gathered the reported patients in three epileptic phenotypes: 15 patients with fever-related epilepsy (40%), 11 with early developmental epileptic encephalopathy (30%), 10 with generalized epilepsy spectrum (27%), and 1 patient without seizures (3%). We did not find a specific phenotype for any gene, but we showed that the location of variants on the transmembrane (TM) segment was associated with a more severe phenotype, irrespective of the GABAA -receptor subunit gene, whereas N-terminal variants seemed to be related to milder phenotypes. SIGNIFICANCE: GABAA -receptor subunit variants are associated with highly variable phenotypes despite their molecular and physiological proximity. None of the genes described here was associated with a specific phenotype. On the other hand, it appears that the location of the variant on the protein may be a marker of severity. Variant location may have important weight in the development of targeted therapeutics.

CHCHD10 and SLP2 control the stability of the PHB complex: a key factor for motor neuron viability.

CHCHD10 is an amyotrophic lateral sclerosis/frontotemporal dementia gene that encodes a mitochondrial protein whose precise function is unclear. Here we show that Coiled-Coil-Helix-Coiled-Coil-Helix Domain Containing protein 10 interacts with the Stomatin-Like Protein 2 and participates in the stability of the prohibitin complex in the inner mitochondrial membrane. By using patient fibroblasts and mouse models expressing the same CHCHD10 variant (p.Ser59Leu), we show that Stomatin-Like Protein 2 forms aggregates with prohibitins, found in vivo in the hippocampus and as aggresome-like inclusions in spinal motor neurons of Chchd10S59L/+ mice. Affected cells and tissues display instability of the prohibitin complex, which participates at least in part in the activation of the OMA1 cascade with OPA1 processing leading to mitochondrial fragmentation, abnormal mitochondrial cristae morphogenesis and neuronal death found in spinal cord and the hippocampus of Chchd10S59L/+ animals. Destabilization of the prohibitin complex leads to the instability of the mitochondrial contact site and cristae organizing the system complex, probably by the disruption of OPA1-mitofilin interaction. Thus, Stomatin-Like Protein 2/prohibitin aggregates and destabilization of the prohibitin complex are critical in the sequence of events leading to motor neuron death in CHCHD10S59L-related disease.

GenomeMixer and TRUST: Novel bioinformatics tools to improve reliability of Non-Invasive Prenatal Testing (NIPT) for fetal aneuploidies.

Non-invasive prenatal testing (NIPT) screens for common fetal chromosomal abnormalities through analysis of circulating cell-free DNA in maternal blood by massive parallel sequencing. NIPT reliability relies on both the estimation of the fetal fraction (ff) and on the sequencing depth (sd) but how these parameters are linked is unknown. Several bioinformatics tools have been developed to determine the ff but there is no universal ff threshold applicable across diagnostics laboratories. Thus, we developed two tools allowing the implementation of a strategy for NIPT results validation in clinical practice: GenomeMixer, a semi-supervised approach to create synthetic sequences and to estimate confidence intervals for NIPT validation and TRUST to estimate the reliability of NIPT results based on confidence intervals found in this study. We retrospectively validated these new tools on 2 cohorts for a total of 1439 samples with 31 confirmed aneuploidies. Through the analysis of the interrelationship between ff, sd and chromosomal aberration detection, we demonstrate that these parameters are profoundly connected and cannot be considered independently. Our tools take in account this critical relationship to improve NIPT reliability and facilitate cross laboratory standardization of this screening test.

Cardiac dyspnea risk zones in the South of France identified by geo-pollution trends study.

The incidence of cardiac dyspnea (CD) and the distribution of pollution in the south of France suggests that environmental pollution may have a role in disease triggering. CD is a hallmark symptom of heart failure leading to reduced ability to function and engage in activities of daily living. To show the impact of short-term pollution exposure on the increment of CD emergency room visits, we collected pollutants and climate measurements on a daily basis and 43,400 events of CD in the Région Sud from 2013 to 2018. We used a distributed lag non-linear model (DLNM) to assess the association between air pollution and CD events. We divided the region in 357 zones to reconciliate environmental and emergency room visits data. We applied the DLNM on the entire region, on zones grouped by pollution trends and on singular zones. Each pollutant has a significant effect on triggering CD. Depending on the pollutant, we identified four shapes of exposure curves to describe the impact of pollution on CD events: early and late effect for NO2; U-shape and rainbow-shape (or inverted U) for O3; all the four shapes for PM10. In the biggest cities, O3 has the most significant association along with the PM10. In the west side, a delayed effect triggered by PM10 was found. Zones along the main highway are mostly affected by NO2 pollution with an increase of the association for a period up to 9 days after the pollution peak. Our results can be used by local authorities to set up specific prevention policies, public alerts that adapt to the different zones and support public health prediction-making. We developed a user-friendly web application called Health, Environment in PACA Region Tool (HEART) to collect our results. HEART will allow citizens, researchers and local authorities to monitor the impact of pollution trends on local public health.

NiPTUNE: an automated pipeline for noninvasive prenatal testing in an accurate, integrative and flexible framework.

Noninvasive prenatal testing (NIPT) consists of determining fetal aneuploidies by quantifying copy number alteration from the sequencing of cell-free DNA (cfDNA) from maternal blood. Due to the presence of cfDNA of fetal origin in maternal blood, in silico approaches have been developed to accurately predict fetal aneuploidies. Although NIPT is becoming a new standard in prenatal screening of chromosomal abnormalities, there are no integrated pipelines available to allow rapid, accurate and standardized data analysis in any clinical setting. Several tools have been developed, however often optimized only for research purposes or requiring enormous amount of retrospective data, making hard their implementation in a clinical context. Furthermore, no guidelines have been provided on how to accomplish each step of the data analysis to achieve reliable results. Finally, there is no integrated pipeline to perform all steps of NIPT analysis. To address these needs, we tested several tools for performing NIPT data analysis. We provide extensive benchmark of tools performances but also guidelines for running them. We selected the best performing tools that we benchmarked and gathered them in a computational pipeline. NiPTUNE is an open source python package that includes methods for fetal fraction estimation, a novel method for accurate gender prediction, a principal component analysis based strategy for quality control and fetal aneuploidies prediction. NiPTUNE is constituted by seven modules allowing the user to run the entire pipeline or each module independently. Using two cohorts composed by 1439 samples with 31 confirmed aneuploidies, we demonstrated that NiPTUNE is a valuable resource for NIPT analysis.

Improved detection of mitochondrial DNA instability in mitochondrial genome maintenance disorders.

PURPOSE: Diseases caused by defects in mitochondrial DNA (mtDNA) maintenance machinery, leading to mtDNA deletions, form a specific group of disorders. However, mtDNA deletions also appear during aging, interfering with those resulting from mitochondrial disorders. METHODS: Here, using next-generation sequencing (NGS) data processed by eKLIPse and data mining, we established criteria distinguishing age-related mtDNA rearrangements from those due to mtDNA maintenance defects. MtDNA deletion profiles from muscle and urine patient samples carrying pathogenic variants in nuclear genes involved in mtDNA maintenance (n = 40) were compared with age-matched controls (n = 90). Seventeen additional patient samples were used to validate the data mining model. RESULTS: Overall, deletion number, heteroplasmy level, deletion locations, and the presence of repeats at deletion breakpoints were significantly different between patients and controls, especially in muscle samples. The deletion number was significantly relevant in adults, while breakpoint repeat lengths surrounding deletions were discriminant in young subjects. CONCLUSION: Altogether, eKLIPse analysis is a powerful tool for measuring the accumulation of mtDNA deletions between patients of different ages, as well as in prioritizing novel variants in genes involved in mtDNA stability.

TDP-43 and PINK1 mediate CHCHD10S59L mutation-induced defects in Drosophila and in vitro.

Mutations in coiled-coil-helix-coiled-coil-helix domain containing 10 (CHCHD10) can cause amyotrophic lateral sclerosis and frontotemporal dementia (ALS-FTD). However, the underlying mechanisms are unclear. Here, we generate CHCH10S59L-mutant Drosophila melanogaster and HeLa cell lines to model CHCHD10-associated ALS-FTD. The CHCHD10S59L mutation results in cell toxicity in several tissues and mitochondrial defects. CHCHD10S59L independently affects the TDP-43 and PINK1 pathways. CHCHD10S59L expression increases TDP-43 insolubility and mitochondrial translocation. Blocking TDP-43 mitochondrial translocation with a peptide inhibitor reduced CHCHD10S59L-mediated toxicity. While genetic and pharmacological modulation of PINK1 expression and activity of its substrates rescues and mitigates the CHCHD10S59L-induced phenotypes and mitochondrial defects, respectively, in both Drosophila and HeLa cells. Our findings suggest that CHCHD10S59L-induced TDP-43 mitochondrial translocation and chronic activation of PINK1-mediated pathways result in dominant toxicity, providing a mechanistic insight into the CHCHD10 mutations associated with ALS-FTD.

Accumulation of amyloid precursor protein C-terminal fragments triggers mitochondrial structure, function, and mitophagy defects in Alzheimer's disease models and human brains.

Several lines of recent evidence indicate that the amyloid precursor protein-derived C-terminal fragments (APP-CTFs) could correspond to an etiological trigger of Alzheimer's disease (AD) pathology. Altered mitochondrial homeostasis is considered an early event in AD development. However, the specific contribution of APP-CTFs to mitochondrial structure, function, and mitophagy defects remains to be established. Here, we demonstrate in neuroblastoma SH-SY5Y cells expressing either APP Swedish mutations, or the ß-secretase-derived APP-CTF fragment (C99) combined with ß- and γ-secretase inhibition, that APP-CTFs accumulation independently of Aß triggers excessive mitochondrial morphology alteration (i.e., size alteration and cristae disorganization) associated with enhanced mitochondrial reactive oxygen species production. APP-CTFs accumulation also elicit basal mitophagy failure illustrated by enhanced conversion of LC3, accumulation of LC3-I and/or LC3-II, non-degradation of SQSTM1/p62, inconsistent Parkin and PINK1 recruitment to mitochondria, enhanced levels of membrane and matrix mitochondrial proteins, and deficient fusion of mitochondria with lysosomes. We confirm the contribution of APP-CTFs accumulation to morphological mitochondria alteration and impaired basal mitophagy in vivo in young 3xTgAD transgenic mice treated with γ-secretase inhibitor as well as in adeno-associated-virus-C99 injected mice. Comparison of aged 2xTgAD and 3xTgAD mice indicates that, besides APP-CTFs, an additional contribution of Aß to late-stage mitophagy activation occurs. Importantly, we report on mitochondrial accumulation of APP-CTFs in human post-mortem sporadic AD brains correlating with mitophagy failure molecular signature. Since defective mitochondria homeostasis plays a pivotal role in AD pathogenesis, targeting mitochondrial dysfunctions and/or mitophagy by counteracting early APP-CTFs accumulation may represent relevant therapeutic interventions in AD.

Multi-Omics Approaches to Improve Mitochondrial Disease Diagnosis: Challenges, Advances, and Perspectives.

Mitochondrial diseases (MD) are rare disorders caused by deficiency of the mitochondrial respiratory chain, which provides energy in each cell. They are characterized by a high clinical and genetic heterogeneity and in most patients, the responsible gene is unknown. Diagnosis is based on the identification of the causative gene that allows genetic counseling, prenatal diagnosis, understanding of pathological mechanisms, and personalized therapeutic approaches. Despite the emergence of Next Generation Sequencing (NGS), to date, more than one out of two patients has no diagnosis in the absence of identification of the responsible gene. Technologies currently used for detecting causal variants (genetic alterations) is far from complete, leading many variants of unknown significance (VUS) and mainly based on the use of whole exome sequencing thus neglecting the identification of non-coding variants. The complexity of human genome and its regulation at multiple levels has led biologists to develop several assays to interrogate the different aspects of biological processes. While one-dimension single omics investigation offers a peek of this complex system, the combination of different omics data allows the discovery of coherent signatures. The community of computational biologists and bioinformaticians, in order to integrate data from different omics, has developed several approaches and tools. However, it is difficult to understand which suits the best to predict diverse phenotypic outcome. First attempts to use multi-omics approaches showed an improvement of the diagnostic power. However, we are far from a complete understanding of MD and their diagnosis. After reviewing multi-omics algorithms developed in the latest years, we are proposing here a novel data-driven classification and we will discuss how multi-omics will change and improve the diagnosis of MD. Due to the growing use of multi-omics approaches in MD, we foresee that this work will contribute to set up good practices to perform multi-omics data integration to improve the prediction of phenotypic outcomes and the diagnostic power of MD.

Accuracy of prenatal screening for congenital heart disease in population: A retrospective study in Southern France.

Congenital heart diseases (CHDs) are the most common congenital malformations. The objective of our study was to evaluate the prenatal screening accuracy of congenital heart disease (CHD) in Southern France and to evaluate the impact of a prenatal diagnosis on pregnancies outcomes and neonatal outcomes. We performed a bicentric, retrospective observational study in the southern region over 4 years was conducted between 1 January 2014 and 31 December 2017. All foetuses and children under one year of age with CHD monitored in the UTHs (University Teaching Hospitals) in Marseille and Nice were included. CHD cases were divided into 3 groups: group 1, those with no possible options for anatomical repair; group 2, those with anatomical repair possibilities but that may require neonatal cardiologic management; and group 3, those with anatomical repair possibilities that do not require an emergency neonatal procedure. Among the 249070 deliveries during the study period, 677 CHD cases were included in the study. The overall prenatal screening rate was 71.5%. The screening rates were 97.8%, 63.6%, and 65.9% for groups 1, 2 and 3, respectively. Among group 2 CHD cases, 80% of the transpositions of the great arteries, 56% of the aortic coarctations, and 20% of the total anomalous pulmonary venous returns were detected during the prenatal period. A genetic anomaly was found in 16% of CHD cases. The overall mortality rate was 11.3% with a higher death rate in cases of prenatal screening (17.2% versus 2.1%; p < 0.001). However, when focusing only on children who died of CHD, prenatal screening did not create an impact (56.6% versus 100%, p = 0,140). Our data showed that the prenatal screening rate of CHD appears satisfactory in Southern France. Nevertheless, it could be improved for some CHD. This study did not find any benefit in terms of mortality from prenatal screening for CHD.

Single-fiber studies for assigning pathogenicity of eight mitochondrial DNA variants associated with mitochondrial diseases.

Whole mitochondrial DNA (mtDNA) sequencing is now systematically used in clinical laboratories to screen patients with a phenotype suggestive of mitochondrial disease. Next Generation Sequencing (NGS) has significantly increased the number of identified pathogenic mtDNA variants. Simultaneously, the number of variants of unknown significance (VUS) has increased even more, thus challenging their interpretation. Correct classification of the variants' pathogenicity is essential for optimal patient management, including treatment and genetic counseling. Here, we used single muscle fiber studies to characterize eight heteroplasmic mtDNA variants, among which were three novel variants. By applying the pathogenicity scoring system, we classified four variants as "definitely pathogenic" (m.590A>G, m.9166T>C, m.12293G>A, and m.15958A>T). Two variants remain "possibly pathogenic" (m.4327T>C and m.5672T>C) but should these be reported in a different family, they would be reclassified as "definitely pathogenic." We also illustrate the contribution of single-fiber studies to the diagnostic approach in patients harboring pathogenic variants with low level heteroplasmy.