Skip to main content

Wide diagnostic and genotypic spectrum in patients with suspected mitochondrial disease



Mitochondrial Diseases (MDs) are a diverse group of neurometabolic disorders characterized by impaired mitochondrial oxidative phosphorylation and caused by pathogenic variants in more than 400 genes. The implementation of next-generation sequencing (NGS) technologies helps to increase the understanding of molecular basis and diagnostic yield of these conditions. The purpose of the study was to investigate diagnostic and genotypic spectrum in patients with suspected MD. The comprehensive analysis of mtDNA variants using Sanger sequencing was performed in the group of 83 unrelated individuals with clinically suspected mitochondrial disease. Additionally, targeted next generation sequencing or whole exome sequencing (WES) was performed for 30 patients of the study group.


The overall diagnostic rate was 21.7% for the patients with suspected MD, increasing to 36.7% in the group of patients where NGS methods were applied. Mitochondrial disease was confirmed in 11 patients (13.3%), including few classical mitochondrial syndromes (MELAS, MERRF, Leigh and Kearns-Sayre syndrome) caused by pathogenic mtDNA variants (8.4%) and MDs caused by pathogenic variants in five nDNA genes. Other neuromuscular diseases caused by pathogenic variants in seven nDNA genes, were confirmed in seven patients (23.3%).


The wide spectrum of identified rare mitochondrial or neurodevelopmental diseases proves that MD suspected patients would mostly benefit from an extensive genetic profiling allowing rapid diagnostics and improving the care of these patients.


Mitochondrial diseases (MDs) are a clinically heterogeneous group of genetic disorders caused by defects affecting the function of the mitochondrial respiratory chain and mitochondrial oxidative phosphorylation pathway. The estimated collective prevalence of MDs in adults is approximately 1 in 4,300 (1 in 11,500 in children) [1].

MD may occur at any age and may involve any organ and tissue, most commonly affecting organs that require high energy levels, such as the nervous system, heart, eyes, skeletal muscles, liver, and kidney. This results in a diversity of clinical presentations and a wide spectrum of clinical symptoms including psychomotor delay or regression, epilepsy, cerebellar ataxia, encephalopathy and stroke-like episodes, dystonia, myopathy, exercise intolerance, sensorineural deafness, ophthalmoplegia, ptosis, optic atrophy, pigmentary retinopathy, cardiomyopathy, diabetes mellitus, etc. [1,2,3] Most MDs are progressive and have poor prognosis with high morbidity and mortality. Although management and treatment is limited to supportive care for the vast majority of MDs, a growing portfolio of pharmacological and genetic therapies are in preclinical development or reaching clinical trials [4].

In addition to the clinical heterogeneity, the genetics of MDs is extremely complex, as it is estimated that more than 1,500 mitochondrial proteins are involved in mitochondrial structure, biogenesis and proper respiratory chain functioning [5]. The majority of these proteins are encoded by nuclear DNA (nDNA) and are essential for mtDNA maintenance, mitochondrial dynamics, coenzyme Q10 biosynthesis, assembling of the respiratory chain complexes, activity and turnover, but 13 essential subunits of respiratory chain complexes are encoded by mitochondrial genome. Therefore, mitochondria are under the control of two genomes, and MDs result from pathogenic variants either in nDNA or in mtDNA [6]. 1136 human genes are currently linked to mitochondrial proteome and more than 400 genes are recognized as playing role in MDs following every inheritance pattern, but the list of candidate genes is still growing up [7].

In mtDNA, the variant may not be present in all of the multiple mtDNA copies per cell, the coexistence of wild-type and mutated mtDNA molecules is known as heteroplasmy. Mutant load of mtDNA may vary from tissue to tissue and over time, contributing to the highly variable disease expression. mtDNA variants are maternally inherited with unpredictable levels of heteroplasmy or can also be sporadic or secondary to nuclear genes’ mutations [8].

Given the considerable clinical and genetic heterogeneity associated with MDs, the diagnosis of MDs is really challenging. Some well-defined clinical syndromes are associated with specific mtDNA mutations, however, in the majority of cases, clinical manifestations are diverse, with overlap between conditions, making specific MDs difficult to identify and differentiate clinically. MD patients sharing the same genotype may have limited phenotypic congruence, and conversely, the similar phenotypes can be caused by different genetic mechanisms [9]. Moreover, other genetic and aquired diseases can lead to a secondary respiratory chain deficiency. Currently, there are no reliable biomarkers to diagnose MDs. Diagnosis is typically defined by a complex combination of clinical assessment, blood and cerebrospinal fluid metabolite profiles, brain imaging, tissue histology and enzymology, and specific nuclear gene and/or mtDNA sequencing [8, 10]. Some of those clinical criteria have been formalized into scoring systems [11,12,13]. Nevertheless, those diagnostic methods have been costly, invasive, time-consuming, and in many cases have failed to provide a molecular diagnosis. The identification of the responsible gene is a prerequisite for proper genetic counseling, prenatal or preimplantation genetic diagnosis and sometimes indicates treatment options [3, 6].

Genome-wide approaches including whole exome sequencing (WES) and whole genome sequencing (WGS) in complement with the transcriptome analysis by RNA sequencing (RNA-seq) have overcome the need to prioritize candidate genes for targeted analysis in MDs. Despite the great advantage of these approaches, the diagnostic success rate of massively parallel sequencing is still reliant on careful clinical and biochemical characterization of patients, and to date, more than one out of two patients has no diagnosis, as genetic defect remains unknown [3, 6, 8].

The aim of the present study was to investigate diagnostic and genotypic spectrum in a group of Lithuanian patients with suspected MDs. The study was performed during the process of reforming our diagnostic approach from traditional diagnostic methods to NGS-based diagnostics. Therefore, traditional diagnostic approaches, as well as NGS-based metods were used in the present study.


The study included 83 unrelated patients (50 females and 33 males) with suspected MD. All patients were evaluated and scored based on a modified Nijmegen disease severity score [11, 12]. The structure of the studied group according to the age at presentation of the clinical symptoms and mitochondrial disease criteria (MDC) scores is shown in Table 1.

Table 1 The structure of the study group

mtDNA screening for large-scale rearrangements and complete mtDNA sequencing was performed for all patient samples. Nuclear mitochondrial genes were investigated using targeted gene panel sequencing or WES technology in 9 and 21 cases, respectively (30 cases in total accounting 36% of all patients).

In total, likely causal variants were identified in 18/83 cases (21.7%), conclusive molecular MD or other diagnosis was obtained in 17 of these cases (20.5%) and probable diagnosis in one case (1.2%). Five novel, sixteen known pathogenic or likely pathogenic variants and one known variant of uncertain significance were identified. The overall diagnostic rate according to used investigation methods, the age of disease onset, MDC scores, established disorder (mitochondrial or other disease) and genome changes (mtDNA or nDNA) is summarized in Table 2.

Table 2 Overall diagnostic rate of the study

Mitochondrial disease genes

Mitochondrial genome analysis confirmed MD diagnosis in eight patients. Six known pathogenic variants in MT-TL1, MT-ATP6 and MT-TK causing Mitochondrial encephalomyopathy, lactic acidosis and stroke-like episodes (MELAS), Leigh and Myoclonus epilepsy associated with ragged-red fibers (MERRF) syndromes, respectively, were identified (Table 3). In addition, one patient was confirmed as carrying common pathogenic mtDNA variant causing Leber hereditary optic neuropathy (LHON) although the finding did not fit the clinical phenotype of the patient and NGS analysis was further performed to specify the diagnosis. Furthermore, mtDNA deletion / duplication analysis identified a novel single 5888 bp deletion involving 13 mtDNA genes and confirmed the clinical diagnosis of Kearns-Sayre syndrome (KSS) in one patient [14].

MD diagnosis was also confirmed in three (10.0%) of the 30 patients tested for pathogenic variants in nuclear genes. Pathogenic variants were detected in TRMU, SURF1 and PNPLA8 genes. One additional patient (3.3%) received a probable diagnosis with only one heterozygous pathogenic variant found in RRM2B and pathogenic variants in BTD gene. Although the second RRM2B variant was not identified for this patient, additional mtDNA quantification showed low mtDNA copy number in muscle biopsy (Table 3). In all cases, an autosomal recessive mode of inheritance was observed.

Table 3 Variants identified in mitochondrial disease genes in the patient group

Genes not obviously related to mitochondrial function

Seven patients of 83 received a confirmed genetic diagnosis of other diseases as pathogenic or likely pathogenic variants were identified in genes with no previous evidence of causing a MD. Three cases had dominant inheritance (two autosomal and X-linked), and 3 cases were autosomal recessive (one homozygous and two compound heterozygous). Limb-girdle muscular dystrophy was diagnosed in one patient after identifying two likely causal variants in ANO5 gene, but common mitochondrial variant m.11778G > A associated with LHON was also identified (Table 4).

Table 4 Variants identified in genes not obviously related to mitochondrial function in the patient group


In past decade, rapid development of advanced technologies in genetics and genomics allowed to dramatically increase the identification of the genetic causes of mitochondrial disorders.

This study was performed during a period when traditional research methods were being replaced by NGS technologies, therefore different molecular genetic methods were applied to study a group of patients with clinically suspected MD. The genetic diagnosis was conclusive in 17 of 83 cases, inconclusive in 1 of 83 cases, giving the overall diagnostic rate of 21.7%. Unfortunately, for the rest of the patients the diagnosis still remains undetermined.

The present study allowed to identify pathogenic variants in either mtDNA or nDNA genes. Some of the diagnoses were clinically suspected and consistent with the phenotype. For example, a patient with clinically diagnosed Kearns-Sayre syndrome was found to have a de novo heteroplasmic single 5888 bp mtDNA deletion m.6069_11956del [14]. Common pathogenic mtDNA variants were identified as causing MELAS and Leigh syndromes. However, some of the findings expanded the knowledge of possible phenotypic presentation associated with some genes. For instance, PNPLA8 gene encodes mitochondrial phospholipase iPLA2γ which is involved in mitochondrial membrane lipid metabolism and is required for efficient bioenergetic function [15]. Biallelic pathogenic variants in PNPLA8 gene have been recently associated with a severe mitochondrial neurodegenerative disease in children, manifesting with microcephaly, hypotonia, weakness, epilepsy, global developmental delay, poor weight gain, and lactic acidosis [16, 17]. The patient reported in this study presented with milder phenotype consisting of adolescence onset ataxia and progressive sensorimotor polyneuropathy [18]. Thus, variability in the manifestation of PNPLA8-related disease seemingly exists.

In this study, one patient received a probable diagnosis with only one heterozygous pathogenic variant identified for autosomal recessive MD. We were able to identify only one previously reported [19] heterozygous pathogenic variant c.414_415del p.(Tyr138Ter) in RRM2B gene. The 2-month-old boy suffered from severe weakness, feeding difficulties, progressive deterioration of respiratory function, hearing loss, and generalized hypotonia with minimal spontaneous movements in legs and significantly reduced deep tendon reflexes. Lactate was elevated. A muscle biopsy showed 20% of ragged-red fibers, > 95% fibers negative for cytochrome-c oxidase, and abnormal deposits of lipids. The disease progressed rapidly and the boy died at the age of 3 months. Additionally, partial biotinidase deficiency caused by pathogenic variants in BTD gene was confirmed for the patient [20]. However, the patient’s phenotype was consistent with autosomal recessive infantile-onset RRM2B-related mitochondrial DNA depletion syndrome [21]. The second pathogenic variant in RRM2B was not identified, as pathogenic variants in deep introns or other regulatory sequences were not covered by targeted gene NGS. RNA sequencing approach could be the choice to further test the patient to identify variants causing aberrant splicing, aberrant or even monoallelic expression. For the patients with WES-inconclusive results, transcriptome study by RNA-seq has shown the increase of diagnostic yield by 10–16% for rare mitochondrial and other Mendelian diseases [22, 23].

Other disease was confirmed for seven patients. One patient showed multisystem involvement suggesting the presence of MD, yet, biallelic likely pathogenic variants were identified in YARS1 gene. YARS1 gene encodes tyrosyl-tRNA synthetase that catalyzes the aminoacylation of the amino acid tyrosine to its corresponding tRNA(Tyr) [24]. This activity is essential for translation and protein synthesis. Pathogenic variants in YARS1 gene have been recently associated with infantile-onset multisystem neurologic, endocrine, and pancreatic disease. In other case, the male patient was referred due to ocular symptoms and diagnosed as a carrier of common familial mitochondrial DNA variant m.11778G>A resulting Leber hereditary optic neuropathy. Later he was referred with isolated hyperCKemia, further genetic analysis identified biallelic variants in ANO5 gene. Several other neurodevelopmental disorders, manifesting with developmental delay, hypotonia, movement disorders, mild dysmorphic features were identified. Most of pathogenic variants identified in CACNA1A, DDX3X, KIF1A, and SETD5 genes were de novo with autosomal or X-linked dominant inheritance pattern. The unspecific clinical and biochemical findings of those pediatric patients allowed the clinicians to consider these patients as having possible or even probable MD, but in individuals who do not present with a recognisable clinical syndrome, the suspicion of MD becomes part of the differential diagnosis of complex neurogenetic disorders. WES may be an indispensable tool to differentiate and reach the final diagnosis of those rare cases.

The limited diagnostic rate of the study can be explained by several reasons. In the study, mtDNA was tested in all patients, the identified diagnostic rate for point mtDNA variants using Sanger sequencing method was 7.5%. Such diagnostic rate could be explained by the fact that pathogenic variants in mtDNA account for only part of the causes of MDs. Moreover, the sensitivity of Sanger sequencing method is limited to detect pathogenic variants at low levels of heteroplasmy (e.g. if not appropriate tissue is tested), therefore some of the variants may remain undetected. mtDNA deletion / duplication analysis identified only one single large-scale deletion and confirmed diagnosis of Kearns-Sayre syndrome. Few patients with chronic progressive external ophthalmoplegia (CPEO) were included in the study, therefore, more positive results were expected. However, mtDNA deletions were not identified in those patients, and possible reason for that could be the tested tissue. It has been shown that deletions and duplications are often undetectable in blood samples in patients over 20 years of age [25], so other tissues tested may have been a better option in these patients.

Pathogenic variants in nuclear genes were tested in 30 patients of the study group, and the detection rate of 36.7% was estimated for NGS methods. The present data are comparable with other data reported in the literature, where a diagnostic rate of targeted gene panels for suspected MD is indicated as 7–31%, while WES achieved diagnostic yields of 25–70% [26, 27]. The identification of pathogenic variants may face limitations due to untested DNA regions, including limited gene panel lists or uncovered regulatory regions in WES. Technical limitations of NGS, such as insufficient coverage of certain regions, may also lead to missed variants. Additionaly, interpreting clinical significance of the identified variants and evaluating gene functions is challenging. Some variants, particularly those associated with a specific disease phenotype, are pathogenic, while others require segregation analysis. WES of family trios or functional studies may be more appropriate for analyzing variants of uncertain significance [28].

Studies show that higher diagnostic rates were obtained in groups with complete phenotyping and more accurate selection of patients with suspected MD [8, 19, 29]. In the study, the patients were categorized as having definite, probable, or possible MD score. Genetic diagnosis was established in only 8.7% of patients with possible MD, while the majority of diagnoses were conclusive in patients with probable (40.0%) and definite (33.3%) MD. A significant limitation of the current study was the small proportion of patients tested using NGS methods. The lower diagnostic rate in patients with definite MD compared to probable MD may be explained by the fact that only 3 of 12 (25%) patients with definite MD were tested using NGS methods compared to 12 of 25 (48%) patients with probable MD and 15 of 46 (33%) patients with possible MD. MD was confirmed mostly in patients classified as having definite or probable MD, while other diseases were diagnosed mostly in patients classified as having probable or possible MD. Patients classified as having probable MD had a high prevalence of clinical symptoms suggesting mitochondrial disease, even if later confirmed to have other disease. Our results showed that patients with higher MDC were more likely to receive not only MD diagnosis, but also an overall genetic diagnosis. It is likely that the scoring system classifies not only primary MD patients but also patients with some other neurodevelopmental diseases, presenting with overlapping clinical symptoms of MD or possible secondary mitochondrial dysfunction.


In the study, the overall diagnostic rate was 21.7% for the patients with suspected MD, increasing to 36.7% in the group of patients where NGS techniques were applied. The data suggest that WES should be used as a first-choice genetic diagnostic method for the patients with multisystem, neurodevelopmental or neuromuscular manifestation, as many patients harbor disease causing variants in both MD genes and genes not obviously related to mitochondrial function, and overlapping clinical phenotypes might not always indicate MD.



Eighty three unrelated patients were included in the study: 58 pediatric (under 16 years of age) and 25 adult (over 16 years of age). All patients were evaluated and scored based on a modified mitochondrial disease criteria (MDC) scale [11, 12]. Although the Nijmegen mitochondrial disease criteria were developed in children with mitochondrial disease, they were used for our cohort as the clinical presentations included in these diagnostic criteria reflect the phenotypic variability. Most of the adults enrolled in the study had multi-systemic disorders (≥ 3 affected organs), characteristically affecting tissues with high-energy demands, such as central nervous system, skeletal muscle, eye, and heart. Some patients had a recognizable mitochondrial syndrome.

Ethical considerations

DNA samples and data collection were performed in accordance to the Declaration of Helsinki, the study protocol was approved by the Vilnius Regional Biomedical Research Ethics Committee of Lithuania. Written informed consent was obtained from all patients or their parents before involvement in the study.

Sample preparation

Molecular genetic testing for the patients was performed using total DNA extracted from peripheral blood samples via standard procedures using the phenol-chloroform-isoamyl alcohol extraction method. In few cases DNA samples extracted from urine epithelial cells or muscle biopsies were also available for the study.

mtDNA deletion / duplication analysis

Long-range polymerase chain reactions (LR-PCR) using two different pairs of primers and multiplex ligation-dependent probe amplification (MLPA) method using Mitochondria Salsa MLPA Kit P125 (MRC-Holland, Amsterdam, The Netherlands) were applied for mtDNA deletion / duplication analysis in all patients. Detailed analysis protocols were previously published [14].

mtDNA sequencing

For complete mitochondrial genome sequencing nine overlapping fragments of approximately 1700–3000 bp in length were amplified using Phusion Hot Start II DNA Polymerase according to the manufacturer’s protocol (Thermo Fisher Scientific, USA). Oligonucleotide sequences for PCR amplification are provided in Table S1 (Additional file 1). Each fragment was sequenced with 3–5 forward and 3–5 reverse primers. The sequences of the 62 oligonucleotide sequencing primers are provided in Table S2 (Additional file 1). Sanger sequencing was performed using the BigDye™ Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, USA), and automatic genetic analyser ABI PRISM 3130xl (Applied Biosystems, USA) according to the manufacturer’s protocol.

nDNA panel sequencing

Targeted amplicon NGS for mitochondrial and other neuromuscular disorders was applied for 9 patients with suspected MD. The library preparation was performed using an Ion AmpliSeq Library Kit 2.0 and Ion AmpliSeq™ Neurological Research Panel consisting of 752 genes associated with neurological disorders and including 99 genes related to mitochondrial diseases (Thermo Fisher Scientific, USA; Table S3 (Additional file 1)). Enrichment of exonic sequences was performed with an Ion PGM™ Hi-Q™ View OT2 Kit (Thermo Fisher Scientific, USA) on the Ion OneTouch™ 2 System (Life Technologies, Thermo Fisher Scientific, USA) and sequenced on an Ion PGM™ System (Life Technologies, Thermo Fisher Scientific, USA) using Ion PGM™ Hi-Q™ View Sequencing Kit (Thermo Fisher Scientific, USA) according to the manufacturer’s protocol. Mapping and variants calling were performed using the Ion Torrent Suite™ Server v. 5.0.2 (Thermo Fisher Scientific, USA).

Whole-exome sequencing

Whole exome sequencing (WES) was performed for 21 patients. Primary steps of raw sequencing data i.e. demultiplexing and trimming of the adaptors were performed by the subcontracting NGS provider (CeGaT GmbH, Tübingen, Germany) using high-throughput next-generation Illumina (Illumina, Inc., San Diego, CA, USA) platform. Obtained sequencing data (FASTQ, BAM, VCF files) analysis was further performed in our laboratory using validated in-house bioinformatic pipeline.

Bioinformatic analysis of sequence variants

mtDNA variants

mtDNA sequencing data obtained using Sanger sequencing were analysed using BioEdit v.7.2 The obtained sequences were compared with the revised Cambridge Reference Sequence (rCRS, GenBank number NC_012920.1).

All mtDNA variants obtained using WES were filtered out. mtDNA sequencing read depth in WES data was on average 400–600. mtDNA variants present at lower than 20% heteroplasmy level were not detected using standard analysis algorithm. The GenBank NC_012920.1 reference sequence was used for analysis and reporting of mtDNA variants.

Homoplasmic and heteroplasmic variants were evaluated on the basis of their allele frequency according to Mitomap Frequency, gnomAD 3.1 Frequency and Helix Frequency provided in MITOMAP database [30]. Variants with a population allele frequency < 0.1% were considered significant (with some exceptional known pathogenic variants present in populations at higher frequencies). Variant pathogenicity was evaluated according to criteria proposed for mtDNA variants [31,32,33]. The pathogenicity of mtDNA variants was evaluated using MITOMAP, HmtVar manually curated database [34], NCBI ClinVar [35], Varsome [36] databases as well as scientific literature. The pathogenicity of mitochondrial tRNA variants was evaluated using MitoTIP in silico prediction algorithm [37]. APOGEE [38] meta-predictor and other common in silico tools (SIFT [39], PolyPhen-2 [40], etc.) were used to predict pathogenicity of the variants in protein‐coding genes.

Nuclear DNA variants

Annotation and post-annotation of obtained sequencing data (FASTQ, BAM, VCF files) from both, panel and WES sequencing, were further post-processed on site using validated in-house bioinformatic pipeline (including ANNOVAR [41]). Reads were mapped/aligned to the Human NCBI Build GRCh37 (hg19/2009) reference genome, duplicated reads marked before variant calling and annotation. In WES data, CNVs were annotated separately from SNP/indels using AnnotSV program. Prioritization and interpretation of genome variants were performed using genomic tools and databases provided by ANNOVAR program in up to ± 6 bp of flanking regions and up to 2% minor allele frequency (MAF). Criteria provided by the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) [42], in silico tools and databases (e.g., SIFT, PolyPhen-2, GERP++ [43], CADD [44], ExAC [45], GnomAD [46], 1000 Genome Project data [47], NCBI dbSNP [48], NCBI ClinVar, Varsome, etc.), and the relevant literature were used to assess the pathogenicity of detected variants. Clinically significant variants were reported following HGVS recommendations [49].

Data Availability

The technical data supporting the findings of this study are available within the paper and its Supplementary Information. Oligonucleotide sequences for PCR amplification are provided in Table S1 (Additional file 1). The sequences of the oligonucleotide sequencing primers are provided in Table S2 (Additional file 1). The list of genes tested using targeted amplicon NGS (Ion AmpliSeq™ Neurological Research Panel) is provided in Table S3 (Additional file 1). Clinical and genetic data that support the findings of this study are not openly available due to reasons of sensitivity and are available from the corresponding author upon reasonable request.



Autosomal dominant


Autosomal recessive


Chronic progressive external ophthalmoplegia


Deoxyribonucleic acid




Kearns-Sayre syndrome


Leber hereditary optic neuropathy




Mitochondrial disease


Mitochondrial disease criteria


Mitochondrial encephalomyopathy, lactic acidosis and stroke-like episodes


Myoclonus epilepsy associated with ragged-red fibers


Mitochondrial DNA


Nuclear DNA


Next-generation sequencing


RNA sequencing


Ribonucleic acid


Whole exome sequencing


Whole genome sequencing




  1. Macken WL, Vandrovcova J, Hanna MG, Pitceathly RDS. Applying genomic and transcriptomic advances to mitochondrial medicine. Nat Rev Neurol. 2021;17(4):215–30.

    Article  PubMed  Google Scholar 

  2. Gorman GS, Chinnery PF, DiMauro S, Hirano M, Koga Y, McFarland R, et al. Mitochondrial diseases. Nat Rev Dis Primers. 2016;2(1).

  3. Barcia G, Assouline Z, Magen M, Pennisi A, Rötig A, Munnich A, et al. Improving post-natal detection of mitochondrial DNA mutations. Expert Rev Mol Diagn. 2020;20(10):1003–8.

    Article  CAS  PubMed  Google Scholar 

  4. Pitceathly RDS, Keshavan N, Rahman J, Rahman S. Moving towards clinical trials for mitochondrial diseases. J Inherit Metab Dis. 2021;44(1):22–41.

    Article  PubMed  Google Scholar 

  5. Plutino M, Chaussenot A, Rouzier C, Ait-El-Mkadem S, Fragaki K, Paquis-Flucklinger V, et al. Targeted next generation sequencing with an extended gene panel does not impact variant detection in mitochondrial diseases. BMC Med Genet. 2018;19(1).

  6. Labory J, Fierville M, Ait-El-Mkadem S, Bannwarth S, Paquis-Flucklinger V, Bottini S. Multi-omics approaches to improve mitochondrial disease diagnosis: Challenges, advances, and perspectives. Front Mol Biosci. 2020;7:590842.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Rath S, Sharma R, Gupta R, Ast T, Chan C, Durham TJ, et al. MitoCarta3.0: an updated mitochondrial proteome now with sub-organelle localization and pathway annotations. Nucleic Acids Res. 2021;49(D1):D1541–7.

    Article  CAS  PubMed  Google Scholar 

  8. Riley LG, Cowley MJ, Gayevskiy V, Minoche AE, Puttick C, Thorburn DR, et al. The diagnostic utility of genome sequencing in a pediatric cohort with suspected mitochondrial disease. Genet Med. 2020;22(7):1254–61.

    Article  CAS  PubMed  Google Scholar 

  9. Álvarez-Mora MI, Sánchez A, Rodríguez-Revenga L, Corominas J, Rabionet R, Puig S, et al. Diagnostic yield of next-generation sequencing in 87 families with neurodevelopmental disorders. Orphanet J Rare Dis. 2022;17(1):60.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Kerr M, Hume S, Omar F, Koo D, Barnes H, Khan M, et al. MITO-FIND: a study in 390 patients to determine a diagnostic strategy for mitochondrial disease. Mol Genet Metab. 2020;131(1–2):66–82.

    Article  CAS  PubMed  Google Scholar 

  11. Wolf NI, Smeitink JAM. Mitochondrial disorders: a proposal for consensus diagnostic criteria in infants and children. Neurology. 2002;59(9):1402–5.

    Article  PubMed  Google Scholar 

  12. Morava E, van den Heuvel L, Hol F, de Vries MC, Hogeveen M, Rodenburg RJ, et al. Mitochondrial disease criteria: diagnostic applications in children. Neurology. 2006;67(10):1823–6.

    Article  CAS  PubMed  Google Scholar 

  13. Witters P, Saada A, Honzik T, Tesarova M, Kleinle S, Horvath R, et al. Revisiting mitochondrial diagnostic criteria in the new era of genomics. Genet Med. 2018;20(4):444–51.

    Article  CAS  PubMed  Google Scholar 

  14. Grigalionienė K, Burnytė B, Balkelienė D, Ambrozaitytė L, Utkus A. Kearns-Sayre syndrome case. Novel 5,9 kb mtDNA deletion. Mol Genet Genomic Med. 2023;11(1):e2059.

    Article  CAS  PubMed  Google Scholar 

  15. Stelzer G, Rosen N, Plaschkes I, ZimMerman S, Twik M, Fishilevich S, et al. The GeneCards suite: from gene data mining to disease genome sequence analyses. Curr Protoc Bioinformatics. 2016;54(1):1.

    Article  Google Scholar 

  16. Saunders CJ, Moon SH, Liu X, Thiffault I, Coffman K, LePichon J-B, et al. Loss of function variants in human PNPLA8 encoding calcium-independent phospholipase A2 γ recapitulate the mitochondriopathy of the homologous null mouse. Hum Mutat. 2015;36(3):301–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Shukla A, Saneto RP, Hebbar M, Mirzaa G, Girisha KM. A neurodegenerative mitochondrial disease phenotype due to biallelic loss-of‐function variants in PNPLA8 encoding calcium‐independent phospholipase A2γ. Am J Med Genet A. 2018;176(5):1232–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Burnyte B, Vilimiene R, Grigalioniene K, Adomaitiene I, Utkus A. Cerebellar ataxia and peripheral neuropathy in a family with PNPLA8-associated disease. Neurol Genet. 2023;9(3):e200068.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Pronicka E, Piekutowska-Abramczuk D, Ciara E, Trubicka J, Rokicki D, Karkucińska-Więckowska A, et al. New perspective in diagnostics of mitochondrial disorders: two years’ experience with whole-exome sequencing at a national paediatric centre. J Transl Med. 2016;14(1).

  20. Wolf B, Biotinidase Deficiency. Adam MP, EVerman DB, Mirzaa GM, Pagon RA, Wallace SE, Bean LJ, et al. editors. PubMed. Seattle (WA): University of Washington, Seattle; 1993. Accessed 19 Apr 2023.

  21. Lim AZ, McFarland R, Taylor RW, Gorman GS. RRM2B Mitochondrial DNA Maintenance Defects. University of Washington, Seattle; 2021. Accessed 19 Apr 2023.

  22. Kremer LS, Bader DM, Mertes C, Kopajtich R, Pichler G, Iuso A, et al. Genetic diagnosis of mendelian disorders via RNA sequencing. Nat Commun. 2017;8(1):15824.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Yépez VA, Gusic M, Kopajtich R, Mertes C, Smith NH, Alston CL, et al. Clinical implementation of RNA sequencing for mendelian disease diagnostics. Genome Med. 2022;14(1):38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Entry. - *603623 - Tyrosyl-tRNA synthetase 1; YARS1. OMIM.; Accessed 17 Apr 2023.

  25. Mavraki E, Labrum R, Sergeant K, Alston CL, Woodward C, Smith C et al. Best practice guidelines for the molecular diagnosis of mitochondrial disease. Acgs uk com

  26. Stenton SL, Prokisch H. Genetics of mitochondrial diseases: Identifying mutations to help diagnosis. EBioMedicine. 2020 [cited 2023 Apr 18];56(102784):102784;

  27. Schuler BA, Nelson ET, Koziura M, Cogan JD, Hamid R, Phillips JA 3rd. Lessons learned: next-generation sequencing applied to undiagnosed genetic diseases. J Clin Invest. 2022;132(7).

  28. Forny P, Footitt E, Davison JE, Lam A, Woodward CE, Batzios S, et al. Diagnosing mitochondrial disorders remains challenging in the omics era. Neurol Genet. 2021;7(3):e597.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Puusepp S, Reinson K, Pajusalu S, Murumets Ü, Õiglane-Shlik E, Rein R, et al. Effectiveness of whole exome sequencing in unsolved patients with a clinical suspicion of a mitochondrial disorder in Estonia. Mol Genet Metab Rep. 2018;15:80–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. MITOMAP: A Human Mitochondrial Genome Database. Accessed 14 Apr 2023.

  31. McCormick EM, Lott MT, Dulik MC, Shen L, Attimonelli M, Vitale O, et al. Specifications of the ACMG/AMP standards and guidelines for mitochondrial DNA variant interpretation. Hum Mutat. 2020;41(12):2028–57.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Bris C, Goudenege D, Desquiret-Dumas V, Charif M, Colin E, Bonneau D, et al. Bioinformatics tools and databases to assess the pathogenicity of mitochondrial DNA variants in the field of next generation sequencing. Front Genet. 2018;9:632.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Castellana S, Rónai J, Mazza T. MitImpact: an exhaustive collection of pre-computed pathogenicity predictions of human mitochondrial non-synonymous variants. Hum Mutat. 2015;36(2):E2413–22.

    Article  CAS  PubMed  Google Scholar 

  34. Preste R, Vitale O, Clima R, Gasparre G, Attimonelli M. HmtVar: a new resource for human mitochondrial variations and pathogenicity data. Nucleic Acids Res. 2019;47(D1):D1202–10.

    Article  PubMed  Google Scholar 

  35. Landrum MJ, Lee JM, Benson M, Brown GR, Chao C, Chitipiralla S, et al. ClinVar: improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 2018;46(D1):D1062–7.

    Article  CAS  PubMed  Google Scholar 

  36. Kopanos C, Tsiolkas V, Kouris A, Chapple CE, Albarca Aguilera M, Meyer R, et al. VarSome: the human genomic variant search engine. Bioinformatics. 2019;35(11):1978–80.

    Article  CAS  PubMed  Google Scholar 

  37. Sonney S, Leipzig J, Lott MT, Zhang S, Procaccio V, Wallace DC, et al. Predicting the pathogenicity of novel variants in mitochondrial tRNA with MitoTIP. PLoS Comput Biol. 2017;13(12):e1005867.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Castellana S, Fusilli C, Mazzoccoli G, Biagini T, Capocefalo D, Carella M, et al. High-confidence assessment of functional impact of human mitochondrial non-synonymous genome variations by APOGEE. PLoS Comput Biol. 2017;13(6):e1005628.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009;4(7):1073–81.

    Article  CAS  PubMed  Google Scholar 

  40. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7(4):248–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405–24.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Davydov EV, Goode DL, Sirota M, Cooper GM, Sidow A, Batzoglou S. Identifying a high fraction of the human genome to be under selective constraint using GERP++. PLoS Comput Biol. 2010;6(12):e1001025.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Rentzsch P, Witten D, Cooper GM, Shendure J, Kircher M. CADD: predicting the deleteriousness of variants throughout the human genome. Nucleic Acids Res. 2019;47(D1):D886–94.

    Article  CAS  PubMed  Google Scholar 

  45. Karczewski KJ, Weisburd B, Thomas B, Solomonson M, Ruderfer DM, Kavanagh D, et al. The ExAC browser: displaying reference data information from over 60 000 exomes. Nucleic Acids Res. 2017;45(D1):D840–5.

    Article  CAS  PubMed  Google Scholar 

  46. Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020;581(7809):434–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. The 1000 Genomes Project Consortium, Auton A, Abecasis GR, Altshuler DM et al. (Co-Chair), Durbin RM (Co-Chair), Abecasis GR,. A global reference for human genetic variation. Nature. 2015;526(7571):68–74;

  48. Sherry ST, Ward M, Sirotkin K. dbSNP-database for single nucleotide polymorphisms and other classes of minor genetic variation. Genome Res. 1999;9(8):677–9.

    Article  CAS  PubMed  Google Scholar 

  49. Sequence Variant Nomenclature. Accessed 14 Apr 2023.

Download references


We are grateful to all of the patients and their families for their kind cooperation.


The reasearch was funded by the Taiwan-Latvian-Lithuanian Collaboration Project ‘Functional model for the mitochondrial disease evaluation and biomarker development’, grant number TAP LLT-02/2015, Research Council of Lithuania.

Author information

Authors and Affiliations



Study concepts and design: B.B., K.G. Clinical information collection: B.B. Data acquisition: K.G. Data analysis/interpretation: K.G., L.A. Manuscript preparation: K.G. Manuscript editing/revision/review: L.A., B.B., A.U. Principal investigator: A.U. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kristina Grigalionienė.

Ethics declarations

Ethics approval and consent to participate

This study was performed in accordance to the Declaration of Helsinki, the study protocol was approved by the Vilnius Regional Biomedical Research Ethics Committee of Lithuania. Written informed consent was obtained from all patients or their parents before involvement in the study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Additional file 1: Table S1

. Oligonucleotide sequences for PCR amplification of mtDNA. Table S2. Oligonucleotide sequencing primers for mtDNA Sanger sequencing. Table S3. A list of 751 targeted genes associated with mitochondrial and other neuromuscular disorders tested using targeted gene next generation sequencing technology (Ion AmpliSeq™ Neurological Research Panel). MD-associated genes (99) are bolded

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Grigalionienė, K., Burnytė, B., Ambrozaitytė, L. et al. Wide diagnostic and genotypic spectrum in patients with suspected mitochondrial disease. Orphanet J Rare Dis 18, 307 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: