- Case report
- Open Access
- Open Peer Review
Novel likely pathogenic variants in TMEM126A identified in non-syndromic autosomal recessive optic atrophy: two case reports
BMC Medical Geneticsvolume 20, Article number: 62 (2019)
Reports on autosomal recessive optic atrophy (arOA) are sparse and so far, only one gene has been specifically associated with non-syndromic arOA, namely TMEM126A. To date, all reports of pathogenic TMEM126A variants are from affected individuals of Maghrebian origin, who all carry an identical nonsense variant. Here we report two novel variants in the TMEM126A gene from non-Maghreb individuals, both found in affected individuals with an arOA phenotype.
We report three affected individuals from two families. The proband of family A, a 24-year-old Turkish woman, was diagnosed with visual loss in early childhood but a diagnosis of optic atrophy was only made at 14 years. A diagnostic gene panel revealed a splice donor variant (c.86 + 2 T > C) in homozygous state in the TMEM126A gene. Analysis of this variant based on RNA from whole blood revealed a single aberrant transcript lacking exon 2, presumably representing a functional null allele. Two siblings from family B, a 16-year old Iraqi girl and her 14-year old brother, were diagnosed with optic atrophy in early childhood. A missense variant p.(S36 L) in the TMEM126A gene was identified in homozygous state in a gene panel-based diagnostic setting in both siblings. This missense variant is ultra rare in the general population, affects a highly evolutionarily conserved amino acid and segregates with the disease within the family. The three probands reported in this study had a relatively mild clinical course without any evidence of a syndromic (e.g. neurological) comorbidity, which is in line with previous studies.
We provide additional evidence for the implication of biallelic pathogenic TMEM126A variants in arOA. Our findings extend both the mutational spectrum and geographic presence of TMEM126A in arOA. Screening of the entire gene should be considered in affected individuals presenting with features resembling arOA and also from non-Maghrebian descent.
Inherited optic neuropathies are rare ocular diseases that are characterized by early-onset, slowly progressive and bilateral vision loss, central scotoma and color vision disturbances due to dysfunction and degeneration of retinal ganglion cells and the optic nerve . Most of the published genes and pathogenic variants have been associated with autosomal-dominant optic atrophy (adOA, MIM 165500) or maternally inherited Leber hereditary optic neuropathy (LHON, MIM 535000). To this day, only one gene has been specifically associated with non-syndromic autosomal-recessive atrophy (arOA), namely TMEM126A . This gene has been mapped to chromosome 11q14.1, covers a genomic region of 8.5 kB (GenBank NM_032273) and is composed of five exons. So far, all affected individuals reported in the literature originate from the same region, the Maghreb, and are homozygous for a recurrent nonsense variant (p.R55*) in the TMEM126A gene [2,3,4].
There is an intimate relationship between optic atrophy disease genes and mitochondrial function, since virtually all publicly known disease genes encode proteins that are localized in mitochondria and/or play a role in mitochondrial homeostasis. The exact function of the TMEM126A protein remains unknown. However, extensive studies performed by Hanein and coworkers have shown that it is a protein of the mitochondrial inner membrane with its N- and C-termini facing the mitochondrial matrix . While affected individuals with pathogenic variants in mitochondrial DNA often clinically present with a multi-systemic disorder, pathogenic variants in nuclear-encoded mitochondrial proteins may be associated with single-system dysfunctions such as adOA caused by pathogenic variants in OPA1, a nuclear encoded mitochondrial large GTPase of the dynamin family that plays a central role in mitochondrial dynamics and cristae junction maintenance [6,7,8]. Progressive visual loss is the most common symptom in these affected individuals resulting from a dysfunction and eventually loss of retinal ganglion cells and their axons. Recent findings, however, showed that the disease spectrum due to pathogenic OPA1 variants is much broader and includes a variety of syndromic diseases such as adOAplus featuring hearing loss, cerebellar ataxia, ptosis, and peripheral neuropathy or a form of Behr syndrome with additional neuromuscular deficits, partly even with biallelic inheritance [9, 10]. While individuals with pathogenic variants in TMEM126A have been initially reported to be affected with non-syndromic optic atrophy , subsequent studies have shown that some affected individuals also presented with mild extra-ocular symptoms such as auditory neuropathy  or sensory-motor axonal neuropathy .
In this study, we report three probands from two families harboring two novel putative pathogenic variants in TMEM126A. In-depth neurological examination revealed no syndromic involvement in these affected individuals. Our results emphasize the need for sequencing the entire TMEM126A gene in affected individuals with arOA instead of focusing on the recurrent (p.R55*) variant, and testing also affected individuals outside from Maghreb. Moreover, the identification of a first missense variant may provide a hint on protein domains crucial for the mitochondrial function of TMEM126A.
Family A was recruited and clinically examined at the Centre for Ophthalmology and at the Department of Neurodegenerative Diseases, University of Tübingen, Germany. Ethical approval (dated October 23, 2017; Project number 637/2017BO1) and informed consent for all individuals was obtained. DNA from whole blood was extracted from all family members using standard protocols. Total RNA was isolated from whole blood drawn in PAXgene tubes using the PAXgene blood RNA Kit (Qiagen, Hilden, Germany).
Family B was recruited and clinically examined at the Center for Ophthalmology and at the Institute of Human Genetics, University Medical Center Hamburg-Eppendorf, Hamburg, Germany. Ethical approval (dated March 1, 2016; Project number PV3802) and informed consent for all individuals including research purposes was obtained. DNA from whole blood was extracted from all family members using standard protocols.
Ophthalmic examination of the affected probands of family A (II:1) and family B (II:1 and II:2) included fundus photography, best corrected visual acuity (BCVA, Snellen), measurement of intraocular pressure, perimetry and spectral domain optical coherence tomography (SD-OCT; Heidelberg Engineering GmbH, Heidelberg, Germany). Details on methodology have already been described [10, 11]. Neurological examination of the affected probands of family A (II:1) and family B (II:1, II:2) included assessment of cognitive function, primary sensory modalities, reflexes, coordination and gait. Nerve conduction studies were performed in subject A (II:1). Magnetic resonance imaging (MRI) of the head was also performed in proband II:1 of family A and proband II:1 of family B.
Panel sequencing and data filtering
Sequencing libraries were prepared starting from DNA using the SureSelectXT workflow (Agilent, Santa Clara, CA) and a custom-design enrichment kit (CeGaT EYE, version 8) for the index patients. Library preparation and capture was performed according to the manufacturer’s instructions and paired-end sequencing was performed on a HiSeq4000 instrument (Illumina, San Diego, CA) (for family A) or on a NovaSeq6000 instrument (Illumina, San Diego, CA) (for family B) with 2 × 100 base pairs (bp) read length, yielding 5 Gbp and 5.3 Gbp of data, respectively. After demultiplexing (Illumina bcl2fastq 2.19), adapters were trimmed with Skewer version 0.1.116 . Trimmed raw reads were aligned to the human genome (hg19) with the Burrows-Wheeler Aligner (BWA-mem version 0.7.2) . Reads likely arising from PCR duplication as well as reads that aligned at more than one locus were discarded (CeGaT internal software). Average coverage on target was 522x (99.2% >10x, 99.53% >30x) for family A and 1145x (99.7% >10x, 99.7% >30x) for family B. Sequence variants were called (VarScan 2.4.2, CeGaT extended version) with a minimum variant allele frequency of 5%. Calls resulting from technical artefacts were removed (CeGaT internal software). Resulting variants were annotated with population frequencies from dbSNP (release 149 and 150) , ExAC version 0.3.1  (for family A), gnomAD version 2.0.1  (for family B) and an internal database (CeGaT), with functional predictions from dbNSFP (3.0c and 3.4c) , with publications from HGMD (16.4 and 17.2) , and with transcript information from Ensembl , RefSeq , and CCDS . Variants were filtered to remove frequent variants (minor allele frequency 1.5%). All variants were manually assessed before inclusion in the final report. Analysis of family members (segregation analysis) was performed by standard Sanger sequencing.
Total RNA isolated from whole blood was reverse transcribed using the Transcriptor First Strand cDNA Synthesis Kit (Roche, Mannheim, Germany) according to the manufacturer’s recommendations. The cDNA was PCR amplified using a forward primer located in exon 1 of TMEM126A (GTGGCTGAGGAAGGAGGAG) and a reverse primer located in exon 4 (CAGTCCACTCCGTGTTATGG). Sequence analysis of RT-PCR products was performed as described previously .
The affected proband from family A is a 24-year-old Turkish woman, born to unaffected parents (Fig. 1, II:1 in the pedigree). Consanguinity was not reported. Family history revealed three healthy sisters. At age 6, the proband noticed a decrease in visual acuity but a diagnosis of optic atrophy was only made at 14 years of age. Fundus examination at 24 years of age revealed pale discs (see Fig. 2a + b). Visual acuity was reduced to 20/100 (OD) and 20/200 (OS). Optical coherence tomography showed a markedly reduced peripapillary retinal nerve fiber layer in the temporal sector on both eyes (see Fig. 3a + b). Perimetry showed centrocecal relative scotomas consistent with papillary atrophy (see Fig. 4a).
Neurological in-depth examination by a neurological expert with long-standing experience in examining neurogenetic patients (M.S.) revealed no clinical signs of damage to any neurological system, in particular not of the pyramidal, basal ganglia, cerebellar, spinal or peripheral nerve system. Nerve conduction studies showed normal results for sural and tibial nerves. Also her cognitive abilities were within the full normal range, with the proband finishing regular school and working full-time on the regular job market as a cashier in a grocery store.
The affected probands from family B are a 16-year old Iraqi girl and her 14-year old brother, born to unaffected consanguineous parents (Fig. 1, Family B, I:1 and I:2). Family history reported four healthy sisters (II:3–5; II:7) and one healthy brother (II:6). A grandfather on the mother’s side was reported to suffer from blindness at the age of 20 years. The girl was initially diagnosed with poor vision and achromatopsia in early childhood. The clinical diagnosis was revised to optic atrophy when visual evoked potentials showed reduced potentials and delayed latency periods at the age of 14. Fundus examination at age 15 revealed temporal pallor of the optic discs and the macular reflex was diminished (see Fig. 2c + d). Visual acuity at this point was reduced to 20/100 (proximity) and 50/100 (distance). Peripapillary retinal nerve fiber layer showed temporal thinning similar to proband A (see Fig. 3c + d). Perimetry showed central and paracentral defects. Additionally, arcuated defects were recorded within the temporal visual fields (see Fig. 4b + c) that could not be attributed to the physiognomy of the proband. The boy was diagnosed during infancy with reduced visual acuity in one eye and was reported to suffer from nyctalopia. Ophthalmological examination revealed slight temporal pallor of the optic discs and the macular reflex was diminished (see Fig. 2e + f). Visual acuity at this point was 25/100 (proximity) and 10/100 (distance). Peripapillary retinal nerve fiber layer showed temporal thinning (see Fig. 3e + f). Perimetry showed paracentral defects (see Fig. 4d + e).
Neurological in-depth examination of both children by a trained neuropediatrician and a skilled clinical geneticist revealed no signs of neurological disease or dysfunction, in particular not of the pyramidal, basal ganglia, cerebellar, spinal or peripheral nerve system. Their cognitive abilities were within the full normal range, with both affected individuals effectively attending regular schooling.
DNA of both index patients was analyzed using a panel-based sequencing approach in a diagnostic setting. The gene panel targets 16 genes that are associated with optic atrophy, namely OPA1, OPA3, TMEM126A, WFS1, MFN2, TIMM8A, SPG7, NR2F1, ACO2, RTN4IP1, AFG3L2, C12ORF65, SLC25A46, CISD2, DNM1L, and YME1L1. No rare and potentially disease-causing exonic variants explaining the disease phenotype were found except for a canonical splice donor site variant in TMEM126A in family A and a missense variant in the same gene in family B. Both variants segregate with the disease in the families as can be seen in the pedigrees shown in Fig. 1.
The splice donor variant c.86 + 2 T > C found in family A, which is not present in large population databases (dbSNP, gnomAD browser), directly changes the invariant 5′ splice core GT sequence of intron 2, ablating the site, with anticipated effects on mRNA splicing. Subsequent cDNA analysis based on RNA isolated from whole blood revealed a single transcript with complete skipping of exon 2 (see Fig. 5). As exon 2 contains the translation start site, the variant c.86 + 2 T > C presumably results in the failure of translation from its natural start site. There is a shorter isoform of TMEM126A, using an ATG in exon 3, which codes for a protein of 125 aa (Ensembl Transcript ENST00000528105.5). This isoform is observed in normal controls, though the larger isoform (Ensembl Transcript ENST00000304511.6) with the ATG in exon 2 is much more abundant (see Fig. 5a). This larger isoform is completely lacking in the proband homozygous for the c.86 + 2 T > C variant.
The p.(S36 L) variant identified in family B is extremely rare in the general population with a minor allele frequency of 8.122e-6 (sourced from gnomAD). Two different online tools were used to assess the potential pathogenicity of the p.(S36 L) variant, namely the UMD predictor (version March 2018)  and the Combined Annotation Dependent Depletion tool (version 1.2) . Within the UMD predictor, scores range from 0 to 100, with scores > 74 indicating pathogenicity of a variant. The Combined Annotation Dependent Depletion tool uses a different scoring system with scores > 20 predicting pathogenicity of a variant. The analysis of the p.(S36 L) variant using both tools resulted in scores of 93 and 25, respectively, thereby emphasizing the putative pathogenicity of this novel missense variant.
Discussion and conclusions
The homozygous nonsense variant (p.R55*) in TMEM126A was proposed as the genetic defect underlying arOA in a large multiplex inbred Algerian family . Clinical features included early-onset severe bilateral deficiency in visual acuity, optic disc pallor, and central scotoma. Subsequently, additional cases, all of Maghrebian descent, carrying the same founder mutation were reported [2,3,4].
In the current manuscript, we provide further evidence for the implication of biallelic pathogenic TMEM126A variants in arOA, as we diagnosed two additional cases with presumably pathogenic variants in this gene. The index patient of family A, of Turkish descent, carried a novel canonical splice site variant c.86 + 2 T > C, while two affected siblings from family B, of Iraqi descent, carried a novel missense variant, p.(S36 L). Affected members from both families presented with typical symptoms of optic atrophy but revealed no syndromic involvement.
TMEM126A (NP_115649) is predicted to contain four transmembrane (TM) domains with the first being formed by amino acid (aa) residues 39–57 . The function of the N-terminus (aa 1–38) is completely unknown. One possibility could be that it acts as a mitochondrial targeting leader sequence. However, mitochondrial presequences are usually rich in positively charged aa residues, they have the potential to form amphiphilic alpha-helices and contain consensus motifs for cleavage sites , neither of which holds true for the N-terminus of TMEM126A. In fact, it has been shown that it is the second transmembrane domain (aa 66–85) that promotes localization to the mitochondria . While it seems obvious that the p.(S36 L) variant does not influence the localization to the mitochondria, it may exert a pathogenic effect through another mechanism. Missense variants may not only act at the protein level but also at the nucleotide level by interfering with the correct assembly of the pre-mRNA splicing machinery. However, in silico analysis of the p.(S36 L) variant using the Human Splicing Finder version 3.0  revealed no difference between mutant and reference sequence. Due to the lack of functional data, and according to the guidelines for the classification of sequence variants by the American College of Medical Genetics and Genomics and the Association for Molecular Pathology , the p.(S36 L) variant has to be categorized as variant of unclear significance, despite the fact that it represents an ultra-rare variant in the general population, segregates with the disease in the family and affects an evolutionary conserved amino acid residue. On the other hand, we could demonstrate that the canonical splice site variant c.86 + 2 T > C leads to the expression of only the short transcript isoform which, if translated, would lack the first seventy amino acids of the canonical protein. This would impair the formation of the second transmembrane domain and we hypothesize that this shortened protein would not be located to the mitochondria.
In conclusion, more information on the effect of TMEM126A alterations on mitochondrial function is needed to confirm the exact pathogenic mechanisms that cause arOA in the probands of this study. Taking our data and the previous publications into account, TMEM126A remains the most promising candidate for non-syndromic arOA and should be considered in affected individuals with a corresponding phenotype (non-syndromic optical atrophy with or without minor extra-ocular findings such as cardiomyopathy, hearing loss, MRI alterations), also in affected individuals outside Maghreb.
Autosomal-dominant optic atrophy
Autosomal recessive optic atrophy
Leber Hereditary Optic Neuropathy
Magnetic resonance imaging
Optical coherence tomography
Optic atrophy 1
Retinal nerve fiber layer
Reverse transcription polymerase chain reaction
Transmembrane protein 126A
Carelli V, La Morgia C, Ross-Cisneros FN, Sadun AA. Optic neuropathies: the tip of the neurodegeneration iceberg. Hum Mol Genet. 2017;26(R2):R139–50.
Hanein S, Perrault I, Roche O, Gerber S, Khadom N, Rio M, Boddaert N, Jean-Pierre M, Brahimi N, Serre V, et al. TMEM126A, encoding a mitochondrial protein, is mutated in autosomal-recessive nonsyndromic optic atrophy. Am J Hum Genet. 2009;84:493–8.
Meyer E, Michaelides M, Tee LJ, Robson AG, Rahman F, Pasha S, Luxon LM, Moore AT, Maher ER. Nonsense mutation in TMEM126A causing autosomal recessive optic atrophy and auditory neuropathy. Mol Vis. 2010;16:650–64.
Désir J, Coppieters F, Van Regemorter N, De Baere E, Abramowicz M, Cordonnier M. TMEM126A mutation in a Moroccan family with autosomal recessive optic atrophy. Mol Vis. 2012;18:1849–57.
Hanein S, Garcia M, Fares-Taie L, Serre V, De Keyzer Y, Delaveau T, Perrault I, Delphin N, Gerber S, Schmitt A, et al. TMEM126A is a mitochondrial located mRNA (MLR) protein of the mitochondrial inner membrane. Biochim Biophys Acta. 2013;1830:3719–33.
Delettre C, Lenaers G, Griffoin JM, Gigarel N, Lorenzo C, Belenguer P, Pelloquin L, Grosgeorge J, Turc-Carel C, Perret E, et al. Nuclear gene OPA1, encoding a mitochondrial dynamin-related protein, is mutated in dominant optic atrophy. Nat Genet. 2000;26:207–10.
Cipolat S, Martins de Brito O, Dal Zilio B, Scorrano L. OPA1 requires mitofusin 1 to promote mitochondrial fusion. Proc Natl Acad Sci U S A. 2004;101:15927–32.
Olichon A, Baricault L, Gas N, Guillou E, Valette A, Belenguer P, Lenaers G. Loss of OPA1 perturbates the mitochondrial inner membrane structure and integrity, leading to cytochrome c release and apoptosis. J Biol Chem. 2003;278:7743–6.
Amati-Bonneau P, Valentino ML, Reynier P, Gallardo ME, Bornstein B, Boissière A, Campos Y, Rivera H, de la Aleja JG, Carroccia R, et al. OPA1 mutations induce mitochondrial DNA instability and optic atrophy ‘plus’ phenotypes. Brain. 2008;131:338–51.
Bonifert T, Karle KN, Tonagel F, Batra M, Wilhelm C, Theurer Y, Schoenfeld C, Kluba T, Kamenisch Y, Carelli V, et al. Pure and syndromic optic atrophy explained by deep intronic OPA1 mutations and an intralocus modifier. Brain. 2014;137:2164–77.
Kernstock C, Beisse F, Wiethoff S, Mast A, Krapp E, Grund R, Dietzsch J, Lagrèze W, Fischer D, Schiefer U. Assessment of functional and morphometric endpoints in patients with non-arteritic anterior ischemic optic neuropathy (NAION). Graefes Arch Clin Exp Ophthalmol. 2014;252:515–21.
Jiang H, Lei R, Ding SW, Zhu S. Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinformatics. 2014;15:182.
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25:1754–60.
Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29:308–11.
Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, O'Donnell-Luria AH, Ware JS, Hill AJ, Cummings BB, et al. Exome aggregation consortium. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016;536:285–91.
Liu X, Jian X, Boerwinkle E. dbNSFP: a lightweight database of human non-synonymous SNPs and their functional predictions. Hum Mutat. 2011;32:894–9.
Stenson PD, Mort M, Ball EV, Evans K, Hayden M, Heywood S, Hussain M, Phillips AD, Cooper DN. The human gene mutation database: towards a comprehensive repository of inherited mutation data for medical research, genetic diagnosis and next-generation sequencing studies. Hum Genet. 2017;136:665–77.
Aken BL, Achuthan P, Akanni W, Amode MR, Bernsdorff F, Bhai J, Billis K, Carvalho-Silva D, Cummins C, Clapham P, et al. Ensembl 2017. Nucleic Acids Res. 2017;45(D1):D635–42.
O‘Leary NA, Wright MW, Brister JR, Ciufo S, Haddad D, McVeigh R, Rajput B, Robbertse B, Smith-White B, Ako-Adjei D, et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016;44(D1):D733–45.
Pruitt KD, Harrow J, Harte RA, Wallin C, Diekhans M, Maglott DR, Searle S, Farrell CM, Loveland JE, Ruef BJ, et al. The consensus coding sequence (CCDS) project: identifying a common protein-coding gene set for the human and mouse genomes. Genome Res. 2009;19:1316–23.
Weisschuh N, Wissinger B, Gramer E. A splice site mutation in the PAX6 gene which induces exon skipping causes autosomal dominant inherited aniridia. Mol Vis. 2012;18:751–7.
Salgado D, Desvignes JP, Rai G, Blanchard A, Miltgen M, Pinard A, Lévy N, Collod-Béroud G, Béroud C. UMD-predictor: a high-throughput sequencing compliant system for pathogenicity prediction of any human cDNA substitution. Hum Mutat. 2016;37:439–46.
Kircher M, Witten DM, Jain P, O'Roak BJ, Cooper GM, Shendure J. A general framework for estimating the relative pathogenicity of human genetic variants. Nat Genet. 2014;46:310–5.
Habib SJ, Neupert W, Rapaport D. Analysis and prediction of mitochondrial targeting signals. Methods Cell Biol. 2007;80:761–81.
Desmet FO, Hamroun D, Lalande M, Collod-Béroud G, Claustres M, Béroud C. Human splicing finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res. 2009;37:e67.
Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, Grody WW, Hegde M, Lyon E, Spector E, et al. ACMG laboratory quality assurance committee. 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:405–24.
We thank the families for their participation in the study.
We acknowledge support by the Deutsche Forschungsgemeinschaft and the Open Access Publishing Fund of the University of Tübingen. Matthis Synofzik is supported by the Else-Kröner Fresenius Stiftung and the grant 779257 “Solve-RD” from the EU Horizon 2020 research and innovation program. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The data that support the findings of this study are available from CeGaT GmbH and Praxis für Humangenetik Tübingen but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of CeGaT GmbH and Praxis für Humangenetik Tübingen.
Ethics approval and consent to participate
Samples from all affected individuals and family members were recruited in accordance with the principles of the Declaration of Helsinki and were obtained with written informed consent accompanying the patients’ samples. For the probands of family B who are underage, informed consent for participation in this study was obtained from the probands’ parents. The study was approved by the institutional review board of the Ethics Committee of the University Hospital of Tübingen (study number 637/2017BO1, dated October 23, 2017) and the Ethics Committee of the Hamburg Chamber of Physicians (study number PV3802, dated March 1, 2016).
Consent for publication
Both families have given their written consents for the case report to be published. For the probands of family B who are underage, informed consent for publication of this study was obtained from the probands’ parents.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.