Recessive VARS2 mutation underlies a novel syndrome with epilepsy, mental retardation, short stature, growth hormone deficiency, and hypogonadism
© The Author(s). 2017
Received: 30 August 2017
Accepted: 6 November 2017
Published: 14 November 2017
The Correction to this article has been published in Human Genomics 2017 11:33
Most mitochondrial and cytoplasmic aminoacyl-tRNA synthetases (aaRSs) are encoded by nuclear genes. Syndromic disorders resulting from mutation of aaRSs genes display significant phenotypic heterogeneity. We expand aaRSs-related phenotypes through characterization of the clinical and molecular basis of a novel autosomal-recessive syndrome manifesting severe mental retardation, ataxia, speech impairment, epilepsy, short stature, microcephaly, hypogonadism, and growth hormone deficiency.
A G>A variant in exon 29 of VARS2 (c.3650G>A) (NM_006295) was identified in the index case. This homozygous variant was confirmed by Sanger sequencing and segregated with disease in the family studied. The c.3650G>A change results in alteration of arginine to histidine at residue 1217 (R1217H) of the mature protein and is predicted to be pathogenic.
These findings contribute to a growing list of aaRSs disorders, broadens the spectrum of phenotypes attributable to VARS2 mutations, and provides new insight into genotype-phenotype correlations among the mitochondrial synthetase genes.
Mental retardation (MR) affects approximately 1–3% of the human population. Both genetic and environmental factors are responsible for causing MR, but in almost 60% of the cases, the actual cause cannot be identified precisely [1, 2]. In developed countries, most severe forms of intellectual disability are thought to have a genetic cause . Mutations in more than 400 genes have been linked to intellectual disability, but most of these mutations have a very low prevalence and their phenotypes are often indistinguishable .
A number of single-gene disorders result in mental retardation. Many of these are associated with atypical or dysmorphic physical characteristics. Such conditions include fragile X syndrome, neurofibromatosis, tuberous sclerosis, Noonan’s syndrome, Cornelia de Lange’s syndrome, and Angelman syndrome(AS) . AS is usually considered as a diagnosis in individuals with severe developmental delay, in combination with seizures, ataxia, hypermotoric behaviors, and absence of speech. Ninety percent of these individuals have an identifiable molecular defect that results in the loss of expression of maternally inherited UBE3A (encoding E6AP) in the brain . Over the last few years, many novel syndromes that mimic AS have been delineated and it is now recognized that some individuals who were clinically diagnosed with “test-negative AS” have one of these AS-like syndromes instead. Their phenotypes include severe intellectual disability, motor delay, variable degrees of speech delay ranging from being completely non-verbal to being able to speak in short sentences, autistic and maladaptive behaviors, short attention span, and seizures .
We report a novel autosomal recessive syndromic disorder with clinical resemblance to AS. The clinical spectrum includes severe mental retardation, ataxia, speech impairment, epilepsy, short stature, microcephaly, flap occiput, protruding tongue, prognathia, tremulous movement of limbs, apparent happy demeanor, easily excitable personality, and excessive chewing mouth behaviors. In addition, there were severe growth hormone deficiency, hypogonadism, and severe osteomalacia. Whole exome sequencing of the affected patients identified a homozygous mutation in VARS2.
Subjects and nucleic acid extraction
Four patients from a consanguineous extended family from the same region of Saudi Arabia were examined by a neurologist (AAS), in the department of Neurosciences, King Faisal Specialist Hospital and Research Center (KFSHRC). They were diagnosed with an Angelman-like syndrome. All subjects were enrolled under a KFSHRC IRB-approved (RAC# 13-MED2056-20) protocol with full-written informed consent. Genomic DNA was extracted from peripheral blood samples using standard procedures (Flexi Gene DNA Handbook, Qiagen). Samples were quantitated spectrophotometrically and stored at − 20 °C.
Linkage and homozygosity mapping
All participating individuals (affected and unaffected) were genotyped using an Affymetrix Axiom array (Affymetrix, Santa Clara, CA, USA) following the manufacturer’s protocol (http://www.affymetrix.com/support/technical/manuals.affx). Resulting genotypes were analyzed for shared runs of homozygosity (ROH) using autoSNPa (http://dna.leeds.ac.uk/autosnpa/). Linkage analysis was performed using the Allegro module of Easy Linkage assuming autosomal recessive inheritance and 100% penetrance (http://genetik.charite.de/hoffmann/easyLINKAGE/index.html).
Whole exome sequencing and analysis
DNA was amplified to obtain a whole exome Ion Proton AmpliSeq library which was further used for emulsion PCR on an Ion OneTouch System followed by an enrichment process using an Ion OneTouch ES, both procedures following the manufacturer’s instructions (Life Technologies, Carlsbad, CA, USA). The template-positive Ion PI Ion Sphere particles were processed for sequencing on the Ion Proton sequencer (Life Technologies, Carlsbad, CA, USA). Approximately 15–17 Gb of sequence was generated per sequencing run. Reads were mapped to UCSC hg19 (http://genome.ucsc.edu/) and variants identified using the Ion Torrent pipeline (Life Technologies, Carlsbad, CA, USA). The resultant variant caller file (vcf) was filtered to include only homozygous variants within the pre-determined ROH (shared by affecteds only) that were absent from dbSNP, 1000 genomes, approximately 2000 Saudi exomes. Included variants were further selected based upon pathogenicity predicted by Polyphen2 (http://genetics.bwh.harvard.edu/pph2/) and SIFT (http://sift.jcvi.org/). Potential causative variants were confirmed by Sanger sequencing and screened further for familial segregation based upon autosomal recessive inheritance. For confirmation, an amplicon incorporating the candidate variant was amplified and sequenced using a BigDye Terminator kit (Applied Biosystems, Foster City, CA) and run on an ABI 3730xl automated sequencer (Applied Biosystems, Foster City, CA). SeqScape v.2.6 software (Applied Biosystems, Foster City, CA) was used to align sequence data and identify variants.
Computational structural analysis of mutants
Sequences were retrieved from the Uniprot database. BLAST and SwissModel  were used to search for suitable structural templates in the Protein Data Bank (PDB). SwissModel, RaptorX , and I-Tasser  were used to produce homology models. Models were manually inspected, and mutations evaluated, using the Pymol program (pymol.org). Disorder and secondary structure elements were predicted using RaptorX. Transmembrane helices were predicted using Phobius . Functional information was compiled from various resources, including Uniprot , InterPro , and publications associated with the model templates used.
Laboratory investigations of the index case showed that he had a normal blood count, normal liver function test, and normal electrolytes and renal function test. His peripheral smear morphology was normal. His calcium was low at 2.3 mmol/l, while his phosphate amylase was normal. He had high alkaline phosphatase at 288 U/L (normal range is 40–135). The patient had normal synacthen test and a low IGF level (26 μg/l) (normal range between 116 and 358). IGF binding protein-3 was 1.8 mg/l, which is extremely low (normal range is 3.4–7.8 mg/l). His FSH was 1.8 IU/L and LH was 0.5 IU/L, both in the lower range. The prolactin level and TSH was normal. Testosterone was extremely low < 0.1 nmol/l (normal range is 9.9–27.8). Free androgen index is low < 0.1% (normal range is 14.8–94.8). The 25-hydroxy vitamin D level is undetectable in blood.
Twenty-four-hour urine collection for oligosaccharides was normal, and mucopolysaccharide was negative in the urine. Chromosomal analysis was 46,XY and attributed as no apparent chromosomal abnormalities. FISH analysis to test SNRPN deletion within the 15q11.2 locus for AS was negative. The patient also had a growth hormone stimulation test, and the result showed a flat response with a fasting growth hormone of 0.6 m U/L. After 10 min, it was 0.6 m U/L and after 1 h was 1.3 m U/L. This indicated severe growth hormone failure.
A skeletal survey was also done for the patient, and the result was in keeping with metabolic bone disease with evidence of severe osteomalacia and multiple looser zones and fractures, microcephaly, scoliosis, and kyphosis. Serum homocysteine was normal. MRI of the brain was performed and concluded that there were relative atrophic changes seen involving the cerebellum; otherwise, there was no evidence of white matter disease, cortical matter disease, or cortex abnormality.
A muscle biopsy was studied for the index patient. It was interpreted as having type II fiber atrophy with glycogen excess associated with a mild increase in mitochondria. Similar findings were observed in his brother and sister. There were no deficiencies in mitochondrial electron transport chain enzymes. Citrate synthase activity was increased, suggesting increased mitochondrial proliferation.
Linkage and homozygosity mapping
Linkage analysis (Fig. 1b) localized the disease to a region on chromosome 6 defined proximally by rs743646 and distally by rs34811366 (chr6:4,536,241-6,670,705) with a maximum LOD score of 4.20. Based upon runs of homozygosity shared by all affecteds and absent in all unaffected individuals (Fig. 1c; red rectangle), we refined the critical interval to a 2.1 Mb region (chr6:4,231,133-6,018,243) containing approximately 60 RefSeq genes (hg19).
Whole exome sequencing and analysis
Structure and function
Aminoacyl-tRNA synthetase enzymes (aaRSs) are involved in translation at the second step in creating a protein from a gene’s instructions. The mRNA interacts with the ribosomal complex which reads the sequence of mRNA bases. Each sequence codes for one particular amino acid, and tRNA assembles the protein, one amino acid at a time. Subsequently, aaRSs help attach the amino acid to the tRNA .
A typical animal cell contains 37 cytoplasmic and mitochondrial synthetase genes, all of which are encoded on nuclear chromosomes. Mitochondrial synthetases are given the same name as cytosolic synthetases, appended with the numeral 2 (e.g., HARS and HARS2) . Mitochondrial and cytoplasmic aaRSs are encoded by distinct nuclear genes, with the exception of GARS and KARS, which are present in both cellular compartments . Mt-tRNA genes are located in three transcription units: the short H-strand unit (with the rRNA region and two tRNA genes) is transcribed twice as fast as the L-strand that contains eight tRNA genes and more frequently the entire H-strand unit that produces 14 tRNAs . Currently, pathological mutations that cause human diseases have been identified in many genes that encode cytoplasmic and mitochondrial tRNA-synthetases. The first cytoplasmic aaRSs mutation associated with a human disease, Charcot-Marie-Tooth, was discovered in glycyl-tRNA synthetase (GARS) . Typically, defects in each mtARS have been identified in one tissue-specific disease (most commonly affecting the brain) or in one syndrome .
VARS2 encodes mitochondrial valyl-tRNA (Val-tRNA) synthetase, which is involved in mitochondrial translation . The deduced 993-amino acid protein has an N-terminal mitochondrial targeting signal with a predicted cleavage site after residue 29. VARS2 has characteristics of a class I mitochondrial aminoacyl-tRNA synthetase, including a classical Rossmann fold .
A homozygous mutation of VARS2 was reported in a sporadic case with microcephaly and epilepsy . It was associated with isolated deficiency of the mitochondrial respiratory chain (MRC) complex I. The disorder was secondary to a homozygous c.1100C>T transition in VARS2, resulting in a threonine 367 to isoleucine (T367I) substitution at a highly conserved residue. The affected patient was an 8-year-old with psychomotor delay, facial dysmorphisms, and microcephaly. At 4 years of age, there was onset of partial seizures. His brain MRI displayed hyperintense lesions in the periventricular regions, the insula, and the right fronto-temporal cortex. Brain proton magnetic resonance spectroscopy revealed a lactate peak in the frontal white matter. Muscle computerized tomography showed bilateral quadriceps hypotrophy, but no evident histological alteration was observed in a muscle biopsy .
In a British boy with combined oxidative phosphorylation deficiency-20, a compound heterozygous mutation of VARS2: a c.1135G>A transition, resulting in an alanine 379 to threonine (A379T) substitution, and a c.1877C-A transversion, resulting in an alanine 626 to aspartic acid substitution (A626D). He presented with progressive external ophthalmoplegia, ataxia, and combined oxidative phosphorylation deficiencies .
Genetic and phenotypic heterogeneity
Diodato et al.
Taylor et al.
NM_001167734.1: c.1100C>T, p.Thr367Ile (Hom)
c.1877C>A, p.Ala626Asp (Het)
Combined oxidative phosphorylation deficiencies
Progressive external ophthalmoplegia
Hyperintense lesions in the periventricular regions, the insula, and the right frontotemporal cortex
Evidence of severe growth hormone deficiency
In our VARS2 family, there were no apparent deficiencies in mitochondrial electron transport chain enzymes. However, citrate synthase activity was increased, suggesting increased mitochondrial proliferation. We identified a homozygous change VARS2:NM_006295:c.3650G>A:R1217H, segregating with disease in our family. This mutation is in the anticodon-binding domain of tRNA synthetase, while the reported isolated VARS2 homozygous mutation  is located in the catalytic core domain of valyl-tRNA synthetase.
The lack of clear biochemical phenotypes of OXPHOS or mitochondrial protein synthesis defects in skin fibroblasts and myoblasts from most of mutant aaRSs patients prevented the use of cultured cells for functional studies. In a few RARS2 patients, fibroblasts showed OXPHOS defects, with RC enzymes variably affected; YARS2 patients presented deficient mitochondrial protein synthesis in myotubes but not in cultured fibroblasts .
In summary, we describe a novel autosomal recessive syndromic disorder associated with a mutation in VARS2 (mitochondrial valyl tRNA synthetase). Clinical features included severe mental retardation, ataxia, speech impairment, epilepsy, short stature, microcephaly, dysmorphia, an easily excitable personality, and excessive chewing mouth behaviors. In addition, there was clinical and chemical evidence of severe growth hormone deficiency, hypogonadism, and severe osteomalacia. Our findings further expand the phenotypic spectrum of Aminoacyl-tRNA synthetase enzymes.
The authors would like to thank the family members for their participation in this study. We thank the Genotyping and Sequencing Core Facilities of King Faisal Specialist Hospital and Research Centre for their contributions. We would also like to acknowledge the support of King Abdulaziz City for Science and Technology (KACST) through the Saudi Human Genome Project (SHGP).
This work was funded through grants from King Faisal Specialist Hospital and Research Centre (RAC 2140029) and King Abdulaziz City for Science and Technology (KACST#13-MED2056-20).
Availability of data and materials
The datasets analyzed during the current study are available from the corresponding author.
AA contributed to the conception and design, acquisition and analysis of data, and drafting of the manuscript. BA contributed to the acquisition and analysis of data and preparation of figures. EG contributed to the acquisition and analysis of data and preparation of figures. DJ contributed to the acquisition and analysis of data and preparation of figures. FB contributed to the acquisition and analysis of data. BFM contributed to the conception and design, acquisition and analysis of data, and drafting of the manuscript and figures. STA contributed to the acquisition and analysis of data. DM contributed to the conception and design, acquisition and analysis of data, and drafting of the manuscript and figures. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Written informed consent was obtained from all enrolled participants in an IRB-approved project (RAC 2140029). The study was carried out in accordance with the declaration of Helsinki under IRB guidelines of the King FaisalSpecialist Hospital and Research Centre (KFSHRC), Riyadh, Saudi Arabia.
Consent for publication
Written informed consent was obtained to publish the video included in the supplement.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Al-Mujadi H, AR AR, Katzarov MG, Dehrab NA, Batra YK, Al-Qattan AR. Preemptive gabapentin reduces postoperative pain and opioid demand following thyroid surgery. Can J Anaesth. 2006;53:268–73.View ArticlePubMedGoogle Scholar
- Fernell E. Aetiological factors and prevalence of severe mental retardation in children in a Swedish municipality: the possible role of consanguinity. Dev Med Child Neurol. 1998;40:608–11.View ArticlePubMedGoogle Scholar
- Ropers HH. Genetics of early onset cognitive impairment. Annu Rev Genomics Hum Genet. 2010;11:161–87.View ArticlePubMedGoogle Scholar
- de Ligt J, Willemsen MH, van Bon BW, Kleefstra T, Yntema HG, Kroes T, Vulto-van Silfhout AT, Koolen DA, de Vries P, Gilissen C, et al. Diagnostic exome sequencing in persons with severe intellectual disability. N Engl J Med. 2012;367:1921–9.View ArticlePubMedGoogle Scholar
- Fryns JP, de Ravel TJ. London Dysmorphology Database, London Neurogenetics Database and Dysmorphology Photo Library on CD-ROM [Version 3] 2001R. M. Winter, M. Baraitser, Oxford University Press, ISBN 019851-780, pound sterling 1595. Hum Genet. 2002;111:113.View ArticlePubMedGoogle Scholar
- Ramsden SC, Clayton-Smith J, Birch R, Buiting K. Practice guidelines for the molecular analysis of Prader-Willi and Angelman syndromes. BMC Med Genet. 2010;11:70.View ArticlePubMedPubMed CentralGoogle Scholar
- Tan WH, Bird LM, Thibert RL, Williams CA. If not Angelman, what is it? A review of Angelman-like syndromes. Am J Med Genet A. 2014;164A:975–92.View ArticlePubMedGoogle Scholar
- Arnold K, Bordoli L, Kopp J, Schwede T. The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006;22:195–201.View ArticlePubMedGoogle Scholar
- Kallberg M, Margaryan G, Wang S, Ma J, Xu J. RaptorX server: a resource for template-based protein structure modeling. Methods Mol Biol. 2014;1137:17–27.View ArticlePubMedGoogle Scholar
- Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER Suite: protein structure and function prediction. Nat Methods. 2015;12:7–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Kall L, Krogh A, Sonnhammer EL. A combined transmembrane topology and signal peptide prediction method. J Mol Biol. 2004;338:1027–36.View ArticlePubMedGoogle Scholar
- UniProt C. UniProt: a hub for protein information. Nucleic Acids Res. 2015;43:D204–12.View ArticleGoogle Scholar
- Hunter S, Jones P, Mitchell A, Apweiler R, Attwood TK, Bateman A, Bernard T, Binns D, Bork P, Burge S, et al. InterPro in 2011: new developments in the family and domain prediction database. Nucleic Acids Res. 2012;40:D306–12.View ArticlePubMedGoogle Scholar
- Francklyn C, Musier-Forsyth K, Martinis SA. Aminoacyl-tRNA synthetases in biology and disease: new evidence for structural and functional diversity in an ancient family of enzymes. RNA. 1997;3:954–60.PubMedPubMed CentralGoogle Scholar
- Abbott JA, Francklyn CS, Robey-Bond SM. Transfer RNA and human disease. Front Genet. 2014;5:158.View ArticlePubMedPubMed CentralGoogle Scholar
- Diodato D, Melchionda L, Haack TB, Dallabona C, Baruffini E, Donnini C, Granata T, Ragona F, Balestri P, Margollicci M, et al. VARS2 and TARS2 mutations in patients with mitochondrial encephalomyopathies. Hum Mutat. 2014;35:983–9.View ArticlePubMedPubMed CentralGoogle Scholar
- King MP, Attardi G. Post-transcriptional regulation of the steady-state levels of mitochondrial tRNAs in HeLa cells. J Biol Chem. 1993;268:10228–37.PubMedGoogle Scholar
- Antonellis A, Ellsworth RE, Sambuughin N, Puls I, Abel A, Lee-Lin SQ, Jordanova A, Kremensky I, Christodoulou K, Middleton LT, et al. 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–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Euro L, Konovalova S, Asin-Cayuela J, Tulinius M, Griffin H, Horvath R, Taylor RW, Chinnery PF, Schara U, Thorburn DR, et al. Structural modeling of tissue-specific mitochondrial alanyl-tRNA synthetase (AARS2) defects predicts differential effects on aminoacylation. Front Genet. 2015;6:21.View ArticlePubMedPubMed CentralGoogle Scholar
- Bonnefond L, Fender A, Rudinger-Thirion J, Giege R, Florentz C, Sissler M. Toward the full set of human mitochondrial aminoacyl-tRNA synthetases: characterization of AspRS and TyrRS. Biochemistry. 2005;44:4805–16.View ArticlePubMedGoogle Scholar
- Taylor RW, Pyle A, Griffin H, Blakely EL, Duff J, He L, Smertenko T, Alston CL, Neeve VC, Best A, et al. Use of whole-exome sequencing to determine the genetic basis of multiple mitochondrial respiratory chain complex deficiencies. JAMA. 2014;312:68–77.View ArticlePubMedGoogle Scholar
- Diodato D, Ghezzi D, Tiranti V. The mitochondrial aminoacyl tRNA synthetases: genes and syndromes. Int J Cell Biol. 2014;2014:787956.View ArticlePubMedPubMed CentralGoogle Scholar