Compound heterozygous mutations in glycyl-tRNA synthetase are a proposed cause of systemic mitochondrial disease
© McMillan et al.; licensee BioMed Central Ltd. 2014
Received: 20 October 2013
Accepted: 12 March 2014
Published: 26 March 2014
Glycyl-tRNA synthetase (GARS) is an aminoacyl-tRNA synthetase (ARS) that links the amino acid glycine to its corresponding tRNA prior to protein translation and is one of three bifunctional ARS that are active within both the cytoplasm and mitochondria. Dominant mutations in GARS cause rare forms of Charcot-Marie-Tooth disease and distal spinal muscular atrophy.
We report a 12-year old girl who presented with clinical and biochemical features of a systemic mitochondrial disease including exercise-induced myalgia, non-compaction cardiomyopathy, persistent elevation of blood lactate and alanine and MRI evidence of mild periventricular leukomalacia. Using exome sequencing she was found to harbor compound heterozygous mutations within the glycyl-tRNA synthetase (GARS) gene; c.1904C > T; p.Ser635Leu and c.1787G > A; p.Arg596Gln. Each mutation occurred at a highly conserved site within the anticodon binding domain.
Our findings suggest that recessive mutations in GARS may cause systemic mitochondrial disease. This phenotype is distinct from patients with previously reported dominant mutations in this gene, thereby expanding the spectrum of disease associated with GARS dysregulation.
KeywordsGlycyl-tRNA synthetase Amino acyl-tRNA synthetase Cardiomyopathy Charcot-Marie-tooth disease
Aminoacyl-tRNA synthetases (ARS) are nuclear-encoded, ubiquitously expressed enzymes that are essential for cytosolic and mitochondrial protein translation . Each ARS links a specific amino acid to its corresponding tRNA, a process known as aminoacylation or “charging” of tRNA. Aminoacylation must take place before protein translation can occur. ARS are being increasingly recognized as having important secondary functions that include regulation of transcription, translation, splicing and apoptosis . ARS mutations have been implicated in a wide range of human diseases. There are 37 different ARS that are divided into three groups depending upon the site of tRNA aminoacylation. Cytoplasm-specific ARS are active within the cell nucleus and cytoplasm. The tRNA-amino acid complexes are transported within the cytoplasm to ribosomes where they facilitate protein translation. Mitochondria-specific ARS are a distinct group of nuclear-encoded proteins that are imported into the mitochondria to carry out tRNA aminoacylation at that site. Three bifunctional ARS have been described that are active in both the cytoplasm and mitochondria.
Glycyl-tRNA synthetase (GARS) is a bifunctional ARS. Structurally, the functional enzyme exists as a homodimer with three functional domains . GARS mutations have been linked to autosomal-dominant Charcot-Marie-Tooth disease type 2D (CMT2D) and distal spinal muscular atrophy type 5 (dSMA5) [3, 4]. Numerous mechanisms have been postulated as to how dysregulation of GARS may give rise to an axonal neuropathy including loss of function, protein aggregation, loss of a secondary ‘housekeeping’ function, and/or mitochondrial dysfunction . The latter is particularly appealing in light of the activity of GARS within mitochondria  as well as the knowledge that other inherited neuropathies have been linked to mitochondrial dysfunction: MFN2 mutations causing CMT2A  and GDAP1 mutations causing CMT4A or CMT2K .
We report a girl who presented with clinical and biochemical features of a systemic mitochondrial disease and using exome sequencing we identified compound heterozygous mutations within the GARS gene. We propose that her phenotype is secondary to recessive mutations in GARS, thereby expanding the phenotypic spectrum associated with mutation of this gene.
A 12 year old girl initially presented at 6 years of age with exercise intolerance. She was born to non-consanguineous healthy parents. She has two younger, healthy brothers although one has autistic spectrum disorder. She presented with shortness of breath with low-intensity aerobic exercise such as jogging or biking. She had no chest pain or palpitations. Cardiac examination at 7 years old noted a normal clinical examination however her electrocardiogram revealed biventricular hypertrophy. Subsequent echocardiogram and cardiac MRI identified thickening of the posterior left ventricle (LV) wall, apex and septum consistent with a non-compaction cardiomyopathy. Bilateral ventricular systolic and diastolic function was normal with a LV ejection fraction of 63%. Stress testing confirmed a normal baseline heart rate, blood pressure and a normal response to exercise. Pulmonary function tests were normal. Biochemical testing revealed normal serum creatine kinase (CK) with a slight elevation of serum troponin T. Genetic testing included normal SCN5A sequencing, chromosomal microarray and hypertrophic and dilated cardiomyopathy panel (GeneDx, Gaithersburg, MD).
Over the next year, she reported exercise-induced myalgia. She had no muscle weakness, cramping or pigmenturia. She could perform short bursts of anaerobic activity without difficulty however, sustained activity would elicit muscle pain. Biochemical testing was abnormal on multiple occasions including: plasma lactate (2.3 – 4.6 mmol/L; normal 0.5 - 2.2 mmol/L) and plasma alanine (603 - 841 μmol /L; normal 152 - 547 μmol/L). Her acylcarnitine profile and carnitine levels (free and total) were normal as was urine organic acid analysis. Repeated serum CK, liver and renal function were normal. Neuromuscular assessment at 9 years of age showed her cranial nerves, muscle power, reflexes, sensory testing and coordination to be within normal limits. Gower manoeuvre was negative and gait was normal. Electrodiagnostic testing confirmed normal right median and sural nerve sensory responses and normal right median and tibial nerve motor responses. Concentric needle electromyography of her right quadriceps was normal.
Muscle biopsy of the quadriceps performed at 9 years of age revealed a preponderance of type 1 fibers although no other microscopic, histochemical or ultrastructural abnormalities were apparent. Muscle respiratory chain enzyme testing and muscle mitochondrial DNA sequencing was normal. MRI of the brain at 10 years of age revealed abnormal T2 and T2FLAIR hyperintensity in the periventricular and trigonal white matter bilaterally. MR spectroscopy of the basal ganglia and subcortical white matter was normal. MRI of the proximal leg muscles was unremarkable. Treatment with ubiquinone, vitamin B50 complex, levocarnitine and ubiquinol were started at 10 years old without apparent clinical effect. Creatine monohydrate was added several months later after which she reported a sustained, subjective clinical improvement in her exercise tolerance.
Neuromuscular evaluation at 12 years old was significant for bilateral extensor hallucis longus and extensor digitorum brevis weakness (4/5). Strength testing of all other muscles was within normal limits. Cranial nerve testing, deep tendon reflexes, sensory testing and coordination were within normal limits. Repeat neurophysiological testing revealed a reduction in her right common peroneal nerve motor response due to a slight CMAP amplitude reduction of 2.1 mV (normal >2.4 mV). Motor responses at the left peroneal, right tibial, right median and ulnar nerves were normal. Sensory responses at the right median, ulnar, superficial peroneal and sural nerves were normal. Concentric needle EMG of the right tibialis anterior and medial gastrocnemius was normal. Her most recent cardiac evaluation at 12 years old identified two new findings: electrocardiogram identified a new subclinical Wolf-Parkinson-White pre-excitation that was not noted on prior studies as well as evidence for LV diastolic dysfunction on echocardiogram.
Her medical history was otherwise unremarkable. She was born at term with no complications. Early milestones were appropriate; she sat at 6 months old, pulled to stand by 12 months and walked independently by 18 months of age. Her growth parameters were stable; height (just <50th percentile) and weight (just < 25th percentile) following along her percentile curves from infancy. She has never had any seizures, headaches or endocrine dysfunction. She has no oculobulbar symptoms and no sensory or autonomic dysfunction. Her visual acuity and hearing were normal. She excelled academically, achieving high grades in a gifted program.
Exome sequencing and analysis
We followed standard manufacturer protocols to perform target capture with the Agilent SureSelect All Exon 50 MB (V3) exome enrichment kit and sequencing of 100 bp paired end reads on Illumina Hiseq 2000, which generated over 12.4 Gb of data for the proband. We removed adaptor sequences and quality trimmed reads using the Fastx toolkit (http://hannonlab.cshl.edu/fastx_toolkit/) and then used a custom script to ensure that only read pairs with both mates present were subsequently used. Reads were aligned to hg19 with BWA 0.5.9  and indel realignment was done using the GATK . Duplicate reads were then marked using Picard (http://picard.sourceforge.net/) and excluded from downstream analyses. We assessed coverage of consensus coding sequence (CCDS) bases using the GATK, which showed that all samples had >91% of CCDS bases covered by at least 10 reads, and >85% of CCDS bases covered by at least 20 reads. Single nucleotide variants (SNVs) and short insertions and deletions (indels) were called using samtools mpileup  with the extended base alignment quality (BAQ) adjustment (-E), and were then quality filtered to require at least 20% of reads supporting the variant call. Variants were annotated using both Annovar  and custom scripts to identify whether they affected protein coding sequence, and whether they had previously been seen in dbSNP132, the 1000 genomes dataset (Nov 2011), the NHLBI GO exomes, or in approximately 500 exomes previously sequenced at our center.
We identified our patient to have a sequence variant in the MIB1 (mindbomb E3 ubiquitin protein ligase 1) gene which has been linked to non-compaction cardiomyopathy. Our patient’s MIB1 sequence variant has not been previously reported but did occur in a highly conserved region of this gene. We confirmed that the proband’s father carries the same MIB1 sequence variant although he shows no evidence of cardiomyopathy on clinical examination or echocardiogram.
Sanger sequencing was used to validate the variants in GARS and to evaluate segregation in the family. Blood samples were obtained and DNA was extracted from the patient as well as her parents and two unaffected siblings. PCR amplification and sequencing was performed with primers 5′CAGATGATCCACCTACCTCAG3′ and 5′ATAACACAGGAAACTGGTTTGTC-3′ to test the c.1787G > A variant. PCR amplification was performed using 5′AGTGAAGATTTGGATTCCCG-3′ and 5′GGACTTGAGAATCTGGGCTC3′ primers and Sanger sequencing was done using 5′AAGAAGCAGTACACATTTCTAAG-3 and 5′GTAAGACAGTAGTTAGATAAC-3 primers, to test the c.1904C > T variant.
Sanger sequencing confirmed the presence of these variants in the proband. The parents were each heterozygous for one of the mutations; c.1787G > A was inherited from her mother and c.1904C > T was inherited from her father (Figure 1). The c.1904C > T mutation has previously been reported to be disease-causing [4, 14], the reported patient exhibited a clinical phenotype characterized by adolescent-onset foot deformity necessitating an orthopedic referral at the age 27 years old.
Examination of family members
Given the findings of recessive mutations in GARS and the knowledge that heterozygous mutations in this gene can cause disease, the family was evaluated in detail. The proband’s father reported no weakness or sensory deficits when examined at 55 years old. His clinical examination was entirely within normal limits. His electrodiagnostic testing, however, revealed evidence of a mild sensorimotor polyneuropathy with axonal features. Bilateral sural and superficial peroneal nerve sensory amplitudes were low. His common peroneal and tibial motor amplitudes were within normal limits. Concentric needle EMG of his right extensor digitorum brevis and abductor hallucis revealed fibrillation potentials and positive sharp waves in addition to chronic neurogenic changes. Needle EMG of his right tibialis anterior also revealed chronic neurogenic changes. The proband’s mother reported no functional difficulty. Her clinical examination, echocardiogram and nerve conduction studies at age 47 years old were within normal limits. Needle EMG was not performed. The proband’s two younger brothers had no significant findings on clinical examination or nerve conduction testing.
Aminoacyl tRNA synthetase mutations are emerging as an important cause of rare childhood and adult diseases. Autosomal dominant mutations within several ARS are known to cause distal motor neuropathy or polyneuropathies, including cytoplasmic ARS such as alanyl-tRNA synthetase (AARS) , tyrosyl-tRNA (YARS)  and lysyl-tRNA synthetase (KARS)  as well as the bifunctional GARS . Additionally, an increasing number of autosomal recessive ARS mutations have been linked to severe clinical phenotypes affecting the central nervous system and other highly metabolically active tissues. Mutations within the mitochondrial aspartyl-tRNAsynthetase (DARS2) have been linked to a syndrome of progressive spastic ataxia with MRI evidence of diffuse subcortical leukoencephalopathy with brainstem and spinal cord involvement [18, 19]. Mutations within another mitochondrial ARS, tyrosyl-tRNA synthetase (YARS2) have been identified in a patient presenting with lactic acidosis, sideroblastic anemia and myopathy associated with a severe deficiency in respiratory chain enzyme function . Our patient demonstrated similar, albeit milder, features of periventricular leukoencephalopathy, lactic acidosis and myalgia suggesting evidence of mitochondrial dysfunction. Each of the two GARS mutations occurred at a highly conserved site within the anticodon binding domain; mutations within the anticodon domain have been previously linked to more severe, earlier onset forms of disease .
Dominant mutations in GARS are a recognized cause of both CMT2D and dSMA5 and our findings suggest that recessive mutations in this gene are associated with mitochondrial disease. Mouse modeling supports our clinical observation that there are fundamental differences in disease phenotype and mechanism in autosomal dominant versus autosomal recessive GARS-associated disease . While mice that were heterozygous for a missense mutation in GARS show a dominant-negative or toxic gain-of-function effect, mice that are homozygous for GARS mutations appear to show some additional loss-of-function effect that is not observed in the heterozygous state . Finally, we cannot exclude the possibility that elements of our proband’s phenotype could secondary to more than GARS dysregulation and that a second concomitant disease is also present. The sequence variant in the MIB1 gene for example raises suspicion that the non-compaction cardiomyopathy could be related to dysfunction in this gene and may or may not be related to the GARS gene dysfunction. However, even if the cardiomyopathy is due to this second gene, there is still substantial evidence for mitochondrial dysfunction in this girl given her persistent lactic acidosis, elevated serum alanine, exercise-induced myalgia and white matter changes on MRI brain.
Our findings suggest that recessive mutations in GARS cause systemic mitochondrial disease. Given our patient’s clinical phenotype, the location of each of the two mutations in GARS and data from the mouse models, we postulate that the clinical phenotype of the patient reported may be explained by a loss-of-function mechanism due to impaired glycyl-tRNA aminoacylation. This phenotype reported here and associated with recessive mutations is distinct from patients with previously reported dominant mutations in this gene, thereby expanding the spectrum of disease associated with GARS dysregulation.
Members of the study family consisted of the proband, mother, father, and siblings. Parents provided informed consent for themselves and their children to be enrolled in the Finding of Rare Disease Genes (FORGE) Canada study. The Research Ethics Board of the Children’s Hospital of Eastern Ontario approved this study in accordance with the Declaration of Helsinki. A copy of the written consent is available for review by the Editor of this journal.
HJM and MTG wrote the manuscript. MTG, PC and KB designed and coordinated the study. HJM, SL, and MTG provided subspecialist consultation services, serial clinical examinations and diagnostic testing. AS and DEB completed Sanger sequencing, modeling and conservation of GARS mutations. JS and JM carried out analysis of the next-generation sequencing data. All authors read and approved the final manuscript.
Sorting Intolerant From Tolerant
Compound muscle action potential.
The authors would like to thank the family for their cooperation and permission to publish these findings. Funding was provided by the Government of Canada through Genome Canada, the Canadian Institutes of Health Research (CIHR) and the Ontario Genomics Institute (OGI-049). Funding was also provided by Genome Québec and Genome British Columbia. KMB is supported by a Clinical Investigatorship Award from the CIHR Institute of Genetics. This work was selected for study by the FORGE Canada Steering Committee which consists of K. Boycott (University of Ottawa), J. Friedman (University of British Columbia), J. Michaud (University of Montreal), F. Bernier (University of Calgary), M. Brudno (University of Toronto), B. Fernandez (Memorial University), B. Knoppers (McGill University), M. Samuels (Université de Montréal), and S. Scherer (University of Toronto).
- Antonellis A, Green ED: The role of aminoacyl-tRNA synthetases in genetic diseases. Annu Rev Genomics Hum Genet. 2008, 9: 87-107. 10.1146/annurev.genom.9.081307.164204.View ArticlePubMedGoogle Scholar
- Park SG, Ewalt KL, Kim S: Functional expansion of aminoacyl-tRNA synthetases and their interacting factors: new perspectives on housekeepers. Trends Biochem Sci. 2005, 30: 569-574. 10.1016/j.tibs.2005.08.004.View ArticlePubMedGoogle Scholar
- Antonellis A, Ellsworth RE, Sambuughin N, Puls I, Abel A, Lee-Lin SQ, Jordanova A, Kremensky I, Christodoulou K, Middleton LT, Sivakumar K, Ionasescu V, Funalot B, Vance JM, Goldfarb LG, Fischbeck KH, Green ED: Glycyl tRNA synthetase mutations in charcot-marie-tooth disease type 2D and distal spinal muscular atrophy type V. Am J Hum Genet. 2003, 72: 1293-1299. 10.1086/375039.View ArticlePubMedPubMed CentralGoogle Scholar
- James PA, Cader MZ, Muntoni F, Childs AM, Crow YJ, Talbot K: Severe childhood SMA and axonal CMT due to anticodon binding domain mutations in the GARS gene. Neurology. 2006, 67: 1710-1712. 10.1212/01.wnl.0000242619.52335.bc.View ArticlePubMedGoogle Scholar
- Motley WW, Talbot K, Fischbeck KH: GARS axonopathy: not every neuron’s cup of tRNA. Trends Neurosci. 2010, 33: 59-66. 10.1016/j.tins.2009.11.001.View ArticlePubMedPubMed CentralGoogle Scholar
- Züchner S, Mersiyanova IV, Muglia M, Bissar-Tadmouri N, Rochelle J, Dadali EL, Zappia M, Nelis E, Patitucci A, Senderek J, Parman Y, Evgrafov O, Jonghe PD, Takahashi Y, Tsuji S, Pericak-Vance MA, Quattrone A, Battaloglu E, Polyakov AV, Timmerman V, Schröder JM, Vance JM: Mutations in the mitochondrial GTPase mitofusin 2 cause charcot-marie-tooth neuropathy type 2A. Nat Genet. 2004, 36: 449-451. 10.1038/ng1341.View ArticlePubMedGoogle Scholar
- Noack R, Frede S, Albrecht P, Henke N, Pfeiffer A, Knoll K, Dehmel T, Meyer Zu Hörste G, Stettner M, Kieseier BC, Summer H, Golz S, Kochanski A, Wiedau-Pazos M, Arnold S, Lewerenz J, Methner A: Charcot-Marie-Tooth disease CMT4A: GDAP1 increases cellular glutathione and the mitochondrial membrane potential. Hum Mol Genet. 2012, 21: 150-162. 10.1093/hmg/ddr450.View ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.View ArticlePubMedPubMed CentralGoogle Scholar
- McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA: The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20: 1297-1303. 10.1101/gr.107524.110.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, 1000 Genome Project Data Processing Subgroup: The sequence alignment/Map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079. 10.1093/bioinformatics/btp352.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang K, Li M, Hakonarson H: ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38: e164-10.1093/nar/gkq603.View ArticlePubMedPubMed CentralGoogle Scholar
- 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: 1073-1081. 10.1038/nprot.2009.86.View ArticlePubMedGoogle Scholar
- Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR: A method and server for predicting damaging missense mutations. Nat Methods. 2010, 7: 248-249. 10.1038/nmeth0410-248.View ArticlePubMedPubMed CentralGoogle Scholar
- Leiden Open Variation Database. http://www.dmd.nl/nmdb2/variants.php?select_db=GARS&action=view&view=0027532%2C0000006%2C1 [Cited 17-June-2013
- Zhao Z, Hashiguchi A, Hu J, Sakiyama Y, Okamoto Y, Tokunaga S, Zhu L, Shen H, Takashima H: Alanyl-tRNA synthetase mutation in a family with dominant distal hereditary motor neuropathy. Neurology. 2012, 78: 1644-1649. 10.1212/WNL.0b013e3182574f8f.View ArticlePubMedPubMed CentralGoogle Scholar
- Jordanova A, Irobi J, Thomas FP, Van Dijck P, Meerschaert K, Dewil M, Dierick I, Jacobs A, De Vriendt E, Guergueltcheva V, Rao CV, Tourney I, Gondim FA, D’Hooghe M, Van Gerwen V, Callaerts P, Van Den Bosch L, Timmermans JP, Robberecht W, Gettemans J, Thevelein JM, De Jonghe P, Kremensky I, Timmerman V: Disrupted function and axonal distribution of mutant tyrosyl-tRNA synthetase in dominant intermediate charcot-marie-tooth neuropathy. Nat Genet. 2006, 38: 197-202. 10.1038/ng1727.View ArticlePubMedGoogle Scholar
- McLaughlin HM, Sakaguchi R, Liu C, Igarashi T, Pehlivan D, Chu K, Iyer R, Cruz P, Cherukuri PF, Hansen NF, Mullikin JC, Biesecker LG, Wilson TE, Ionasescu V, Nicholson G, Searby C, Talbot K, Vance JM, Züchner S, Szigeti K, Lupski JR, Hou YM, Green ED, Antonellis A, NISC Comparative Sequencing Program: Compound heterozygosity for loss-of-function lysyl-tRNA synthetase mutations in a patient with peripheral neuropathy. Am J Hum Genet. 2010, 87: 560-566. 10.1016/j.ajhg.2010.09.008.View ArticlePubMedPubMed CentralGoogle Scholar
- Tzoulis C, Tran GT, Gjerde IO, Aasly J, Neckelmann G, Rydland J, Varga V, Wadel-Andersen P, Bindoff LA: Leukoencephalopathy with brainstem and spinal cord involvement caused by a novel mutation in the DARS2 gene. J Neurol. 2012, 259: 292-296. 10.1007/s00415-011-6176-9.View ArticlePubMedGoogle Scholar
- Synofzik M, Schicks J, Lindig T, Biskup S, Schmidt T, Hansel J, Lehmann-Horn F, Schöls L: Acetazolamide-responsive exercise-induced episodic ataxia associated with a novel homozygous DARS2 mutation. J Med Genet. 2011, 48: 713-715. 10.1136/jmg.2011.090282.View ArticlePubMedGoogle Scholar
- Sasarman F, Nishimura T, Thiffault I, Shoubridge EA: A novel mutation in YARS2 causes myopathy with lactic acidosis and sideroblastic anemia. Hum Mutat. 2012, 33: 1201-1206. 10.1002/humu.22098.View ArticlePubMedGoogle Scholar
- Motley WW, Seburn KL, Nawaz MH, Miers KE, Cheng J, Antonellis A, Green ED, Talbot K, Yang XL, Fischbeck KH, Burgess RW: Charcot-Marie-Tooth-linked mutant GARS is toxic to peripheral neurons independent of wild-type GARS levels. PLoS Genet. 2011, 7: e1002399-10.1371/journal.pgen.1002399.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/15/36/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.