ITPR1 gene p.Val1553Met mutation in Russian family with mild Spinocerebellar ataxia
© Shadrina et al. 2016
Received: 5 November 2015
Accepted: 6 December 2015
Published: 13 January 2016
Spinocerebellar ataxias (SСAs) are a highly heterogeneous group of inherited neurological disorders. The symptoms of ataxia vary in individual patients and even within the same SCA subtype. A study of a four-generation family with autosomal dominant (AD) non-progressive SCA with mild symptoms was conducted. The genotyping of this family revealed no frequent pathogenic mutations. So the objective of this study was to identify the genetic causes of the disease in this family with the technology of whole-exome sequencing (WES).
Methods and results
WES, candidate variant analysis with further Sanger sequencing, mRNA secondary structure prediction, and RSCU analysis were performed; a heterozygous missense mutation in ITPR1 was identified.
Our study confirms the fact that ITPR1 gene plays a certain role in the pathogenesis of SCAs, and, therefore, we suggest that c.4657G>A p.Val1553Met) is a disease-causing mutation in the family studied.
KeywordsSpinocerebellar ataxia ITPR1 Whole-exome sequencing Sanger sequencing
Spinocerebellar ataxias (SСAs) are a highly heterogeneous group of inherited neurological disorders. They are characterized by significant clinical polymorphism resulting from combination of cerebellar ataxia with additional non-cerebellar symptoms, which seriously complicates clinical differentiation . The symptoms of ataxia vary in individual patients and even within the same SCA subtype. In general, affected individuals exhibit progressive cerebellar ataxia variably associated with dysarthria, oculomotor abnormalities, pyramidal and extrapyramidal signs, epilepsy, and mental impairment. Neuronal loss is observed predominantly in the cerebellum and the brainstem, and neuroimaging (CT or MRI) demonstrates the atrophy of those regions. These pathologies can be caused by autosomal dominant (AD), autosomal recessive, and X-linked mutations.
There are many different SCA classifications based on both clinical and anatomical principles. Harding suggested a classification according to the mode of inheritance and clinical signs . The actual numbering corresponds to the order of gene description. This molecular classification is currently the most accepted by the scientific community. To date, 38 types of AD SCAs have been identified , which are difficult to distinguish clinically. There are several common forms of the disease in each geographical region (such as SCA1, SCA2, SCA3 and some others in Russia) with dynamic mutations resulting in most severe clinical manifestations, and these forms are usually diagnosed in clinical practice . Other forms of AD SCAs are rare, occur predominantly in a small number of families, exhibit relatively mild manifestations, and are difficult to define because of some phenotypic overlap between them.
Here, a study of a four-generation family with AD non-progressive SCA with mild symptoms was conducted. The genotyping of this family revealed an absence of pathogenic mutations that are common for the Russian population. In such cases for which the screening of known mutations is negative, whole-exome sequencing (WES) can be used in order to find new or rare mutations in the coding regions of the genome. We applied this technology in our family with AD SCA.
WES was performed in two affected individuals (II-2 and IV-1) and in unaffected family member (II-3) at the NGS Core Facility of the Estonian Genome Center, University of Tartu, Estonia. Exome capture was performed with the TruSeq™ Exome Enrichment Kit (Illumina, USA), according to the manufacturer’s protocol. The captured libraries were sequenced with 100-bp paired-end reads on Illumina platform HiSeq2000. The in-house pipeline of the Estonian Genome Center, University of Tartu, was used for the alignment to the reference genome and variant calling, as described previously . Because the SCA was an autosomal dominant hereditary disorder and only one parent was affected, we focused on heterozygous nonsynonymous, frameshift, and canonical splice-site variants that were absent from public datasets, including dbSNP137 , the 1000 Genomes Project , and the NHLBI Exome Sequencing Project (ESP) Exome Variant Server database . Variants that were shared by two affected patients and were absent in a healthy family member were considered to be potentially related to the disease. Nonsynonymous amino acid variants were analyzed using SIFT , PolyPhen-2 , fathmm , MutationAssessor , and MutationTaster , to assess any potentially damaging effects. Protein-altering SNPs that were predicted to be damaging by at least two methods were considered as candidate causative variants. Relevant mutations in all SCA genes were then prioritized manually.
Validation of selected variants
Candidate variant analysis was performed by Sanger sequencing of all available family members. PCR primers for ITPR1 [NG_016144.1], COL6A3 [NG_008676.1], SCN9A [NG_012798.1], SYNE1 [NG_012855.1], and COL9A1 [NG_011654.1] variants were designed using NCBI Primer-BLAST  (sequences and PCR conditions are available upon request). The purification of PCR products was performed using the Axygen Gel Extraction Kit. Sanger sequencing was performed in both forward and reverse directions on an ABI 3730 DNA Analyzer. Sequences were analyzed using the Sequence Scanner Software v.1.0 and AlignX® module for Vector NTI.
Substitution characterization: mRNA secondary structure prediction and RSCU analysis
The analysis of the secondary structure of the ITPR1 mRNA and of associated Gibbs free energy (ΔG) values was performed essentially as described previously [15–17]. Gibbs free energy calculation was carried out using the mfold software  based on the nearest neighbor free energy model. The RNA structures with the lowest free energy for a given mRNA fragment (51 nt in length) were selected and the Gibbs free energy difference (ΔΔG = mutantΔG – WTΔG) was calculated (with the mutation of interest positioned in the middle of the mRNA fragment analyzed). It was shown previously that analysis of short mRNA fragments (~25–75 nt in length) was the most predictive [15–17].
Calculations of the differences in relative synonymous codon usage (RSCU) values  were carried out as described previously [15–17]. Briefly, the RSCU value for a mutant variant codon and that for the corresponding WT codon were calculated. ΔRSCU = RSCUmutant – RSCUwildtype represents a change in the RSCU values as a consequence of the mutation in the ITPR1 gene and may be indicative of a change in translation rate around a particular codon.
Clinical features of the family
According to the proband’s evidence, a similar disease affected his father (II-2) between 15 and 20 years of age. He had an unsteady gait and scanning speech throughout his life, but was able to live independently and died at the age of 81. The proband’s son (IV-1, 27 years of age) started walking independently only at 3 years of age and often fell during his childhood. Having graduated from high school and university with good results, he pursued studies in hi-tech; now he works successfully as an IT administrator. He also succeeded in learning a second language (English). The patient’s clinical features are very similar to those of his father, but his gait difficulties are more severe. The proband’s brother, who is 52 years old, also suffered from the same disorder and exhibited similar clinical features starting in adolescence.
Next-generation sequencing (NGS) data
A total of 46507918, 50734956, and 40536854 reads were obtained for II:2, IV:1 and II:3, respectively. The average target coverage was 34.0× for II:2, 37.3× IV:1, and 29.0× for II:3.
Identification of disease-causing mutations
WES identified 947 variants that were heterozygous in two affected patients and reference homozygous in one unaffected family member. These heterozygous variants were filtered through the databases SIFT, PolyPhen-2, fathmm, MutationAssessor, and MutationTaster; 12 matching variants were found that had a minor allele frequency (MAF) not exceeding 0.01. Among them, five variants were nonsynonymous and were predicted as damaging by at least two prediction methods.
Validation of identified substitution
Summary of variants revealed by whole-exome sequencing in the patients with SCA
Predicted non-functional (low)
Predicted non-functional (low)
Predicted functional (medium)
Predicted functional (medium)
Predicted non-functional (low)
mRNA secondary structure prediction and RSCU analysis
As mentioned previously, SCAs represent a heterogeneous group of diseases from both the clinical and genetic points of view, as they include forms that are rare and difficult to diagnose. Therefore, a NGS analysis was needed to detect the underlying genetic variation. We performed an analysis of a family with a mild SCA clinical phenotype to identify the genetic variants that may contribute to the development of disease. One possibly pathogenic variant, p.Val1553Met, was identified in exon 37 of the ITPR1 gene. We analyzed this mutation in all family members for whom DNA was available, and in the Russian population controls. The segregation of the identified substitution with the disease and its absence in the control Russian population as well as in the 1000 Genomes Project and the NHLBI ESP Exome Variant Server database were shown.
Prediction tests for this substitution indicated its pathogenic significance, as it was predicted as being probably damaging by PolyPhen-2 and as disease causing by MutationTaster.
Moreover, a recent genetic data analysis of patients with hereditary diseases revealed that not only the information contained at the amino acid sequence level of the target protein (mutant vs wild type) but also the information contained at the mRNA/nucleotide level may serve as a predictor of disease development and severity [16, 17]. In the case of blood disorders such as hemophilia A and B, which are caused by mutations in the f8 and f9 genes, respectively, it has been shown that statistically significant differences between severe and nonsevere disease-causing mutations can be observed at the mRNA level, i.e., severe disease-causing mutations tend to lead to substantial changes in the free energy of small mRNA fragments (harboring the site of mutation), and tend to lead to substantial changes in RSCU values. We evaluated the stability of the mRNA and the RCSU values for the possibly pathogenic variant identified, and found that it leads to a change in the structure and stability of the mRNA (ΔΔG = 0.4), as well as changes in RSCU (ΔRSCU = –0.85). Reported data indicate that such changes may affect the local translation rates, which can in turn perturb protein folding .
Thus, all listed data indicate that the G4684A substitution identified here, which results in a p.Val1553Met amino acid change, is disease causing and leads to the development of a form of SCA with a mild clinical phenotype in the family studied. This was supported by the fact that the involvement of ITPR1 gene in the etiology of SCAs, SCA15 (former SCA16) and SCA29, was already shown [20, 22–24]. The same substitution was described previously in an Australian family, and the authors found another p.N602D substitution in the ITPR1 gene in a Canadian family . The neurological syndrome in the Canadian and Australian families was similar, with obvious anticipation and very slow clinical deterioration over many years. The lack of cognitive retardation in childhood and adolescence and the lack of clinically significant cognitive impairment in adulthood were important differences between our case and those occurring in the Australian and Canadian families.
To date, little is known about the role of the ITPR1 gene in the pathogenesis of SCAs; however, we know that it encodes the type 1 inositol 1,4,5-trisphosphate (IP3) receptor, which is an intracellular IP3-gated calcium channel that modulates intracellular calcium signaling . It releases calcium from the endoplasmic reticulum by binding to specific receptors that are coupled to calcium channels. These receptors are abundant in neuronal and non-neuronal tissues. Val1553 is located within the binding domain for carbonic anhydrase VIII (CA8). Its precise function is unknown, but it is assumed that CA8 competes with IP3 for binding to IP3R and inhibits the release of Ca from cells. It was supposed that the p.Val1553Met affects the cellular homeostasis of calcium, which is important for neuronal cells . In addition, there is evidence of a direct physical interaction between ITPR1 and VCP proteins . This protein is involved in the processes of neurodegeneration: mutations in VCP cause a pleiotropic degenerative disorder called multisystem proteinopathy, which can manifest clinically as classical amyotrophic lateral sclerosis, frontotemporal dementia, inclusion body myopathy, or a combination of these disorders [28, 29]. Furthermore, VCP interacts directly with ATXN3, which plays an important role in the pathogenesis of SCAs [17, 18]. Therefore, we can assume that changes in the amino acid structure of ITPR1 may indirectly affect the activity of VCP or ATXN3, which may contribute to the development of a neurodegenerative process in this type of ataxia.
Thus, all existing data on the ITPR1 genetic variation allows to presume that this gene plays an important role in the pathogenesis of SCAs, and, therefore, we suggest that the p.Val1553Met variant is a disease-causing mutation in the family studied.
This study was supported by the program of the Federal Agency of Scientific Organization “Molecular and Cellular Biology”; the Leading Scientific Schools supporting programs SS_4294.2012.4 and SS_2628.2012.4; the targeted Financing from the Estonian Ministry of Education and Research (grant SF0180142s08), the EU FP7 (grant #313010), the Development Fund of the University of Tartu (grant SP1GVARENG), and the European Regional Development Fund to the Centre of Excellence in Genomics (EXCEGEN; grant 3.2.0304.11-0312). Data analyzes were carried out in part in the High Performance Computing Center of University of Tartu. The authors declare no potential conflicts of interest with respect to the authorship and/or publication of this article.
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.
- Matilla-Duenas A, Corral-Juan M, Volpini V, Sanchez I. The spinocerebellar ataxias: clinical aspects and molecular genetics. Adv Exp Med Biol. 2012;724:351–74. doi:https://doi.org/10.1007/978-1-4614-0653-2_27.PubMedView ArticleGoogle Scholar
- Harding AE. The clinical features and classification of the late onset autosomal dominant cerebellar ataxias. A study of 11 families, including descendants of the ‘the Drew family of Walworth’. Brain. 1982;105(Pt 1):1–28.PubMedView ArticleGoogle Scholar
- Hereditary Ataxias: Dominant [database on the Internet]. Neuromuscular Disease Center 2015. Available from: http://neuromuscular.wustl.edu/ataxia/domatax.html. Accessed: 28.01.2015.
- Rossi M, Perez-Lloret S, Doldan L, Cerquetti D, Balej J, Millar Vernetti P, et al. Autosomal dominant cerebellar ataxias: a systematic review of clinical features. Eur J Neurol. 2014;21(4):607–15.PubMedView ArticleGoogle Scholar
- Nikopensius T, Annilo T, Jagomagi T, Gilissen C, Kals M, Krjutskov K, et al. Non-syndromic tooth agenesis associated with a nonsense mutation in ectodysplasin-A (EDA). J Dent Res. 2013;92(6):507–11. doi:https://doi.org/10.1177/0022034513487210.PubMedView ArticleGoogle Scholar
- dbSNP Short Genetic Variations [database on the Internet]2014. Available from: http://www.ncbi.nlm.nih.gov/projects/SNP/index.html. Accessed: 19.09.2014.
- 1000 Genomes [database on the Internet]2014. Available from: http://www.1000genomes.org/home. Accessed: 6.11.2014.
- NHLBI Exome Sequencing Project (ESP) Exome Variant Server [database on the Internet]2014. Available from: http://evs.gs.washington.edu/EVS/. Accessed: 03.11.2014.
- SIFT [database on the Internet]. Available from: http://sift.jcvi.org/www/SIFT_chr_coords_submit.html. Accessed: Aug 2011.
- PolyPhen-2 prediction of functional effects of human nsSNPs [database on the Internet]. Available from: http://genetics.bwh.harvard.edu/pph2/. Accessed: 15.02.2012.
- fathmm, Functional Analysis through Hidden Markov Models [database on the Internet]. Available from: http://fathmm.biocompute.org.uk/. Accessed: May 2014.
- MutationAssessor [database on the Internet]. Available from: http://mutationassessor.org/. Accessed: 12.09.2009.
- Schwarz JM, Cooper DN, Schuelke M, Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age. Nat Methods. 2014;11(4):361–2. doi:https://doi.org/10.1038/nmeth.2890.PubMedView ArticleGoogle Scholar
- Primer-BLAST. http://www.ncbi.nlm.nih.gov/tools/primer-blast/. Accessed 2015.
- Edwards NC, Hing ZA, Perry A, Blaisdell A, Kopelman DB, Fathke R, et al. Characterization of coding synonymous and non-synonymous variants in ADAMTS13 using ex vivo and in silico approaches. PLoS One. 2012;7(6):e38864. doi:https://doi.org/10.1371/journal.pone.0038864.PubMedPubMed CentralView ArticleGoogle Scholar
- Hamasaki-Katagiri N, Salari R, Simhadri VL, Tseng SC, Needlman E, Edwards NC, et al. Analysis of F9 point mutations and their correlation to severity of haemophilia B disease. Haemophilia. 2012;18(6):933–40. doi:https://doi.org/10.1111/j.1365-2516.2012.02848.x.PubMedView ArticleGoogle Scholar
- Hamasaki-Katagiri N, Salari R, Wu A, Qi Y, Schiller T, Filiberto AC, et al. A gene-specific method for predicting hemophilia-causing point mutations. J Mol Biol. 2013;425(21):4023–33. doi:https://doi.org/10.1016/j.jmb.2013.07.037.PubMedPubMed CentralView ArticleGoogle Scholar
- Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003;31(13):3406–15.PubMedPubMed CentralView ArticleGoogle Scholar
- Sharp PM, Tuohy TM, Mosurski KR. Codon usage in yeast: cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Res. 1986;14(13):5125–43.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang L, Chardon JW, Carter MT, Friend KL, Dudding TE, Schwartzentruber J, et al. Missense mutations in ITPR1 cause autosomal dominant congenital nonprogressive spinocerebellar ataxia. Orphanet J Rare Dis. 2012;7:67. doi:https://doi.org/10.1186/1750-1172-7-67.PubMedPubMed CentralView ArticleGoogle Scholar
- Komar AA. A pause for thought along the co-translational folding pathway. Trends Biochem Sci. 2009;34(1):16–24. doi:https://doi.org/10.1016/j.tibs.2008.10.002.PubMedView ArticleGoogle Scholar
- Hara K, Fukushima T, Suzuki T, Shimohata T, Oyake M, Ishiguro H, et al. Japanese SCA families with an unusual phenotype linked to a locus overlapping with SCA15 locus. Neurology. 2004;62(4):648–51.PubMedView ArticleGoogle Scholar
- Hara K, Shiga A, Nozaki H, Mitsui J, Takahashi Y, Ishiguro H, et al. Total deletion and a missense mutation of ITPR1 in Japanese SCA15 families. Neurology. 2008;71(8):547–51. doi:https://doi.org/10.1212/01.wnl.0000311277.71046.a0.PubMedView ArticleGoogle Scholar
- Iwaki A, Kawano Y, Miura S, Shibata H, Matsuse D, Li W, et al. Heterozygous deletion of ITPR1, but not SUMF1, in spinocerebellar ataxia type 16. J Med Genet. 2008;45(1):32–5.PubMedView ArticleGoogle Scholar
- Berridge MJ. Inositol trisphosphate and calcium signalling. Nature. 1993;361(6410):315–25. doi:https://doi.org/10.1038/361315a0.PubMedView ArticleGoogle Scholar
- Hirota J, Ando H, Hamada K, Mikoshiba K. Carbonic anhydrase-related protein is a novel binding protein for inositol 1,4,5-trisphosphate receptor type 1. Biochem J. 2003;372(Pt 2):435–41. doi:https://doi.org/10.1042/BJ20030110.PubMedPubMed CentralView ArticleGoogle Scholar
- Pearce MM, Wormer DB, Wilkens S, Wojcikiewicz RJ. An endoplasmic reticulum (ER) membrane complex composed of SPFH1 and SPFH2 mediates the ER-associated degradation of inositol 1,4,5-trisphosphate receptors. J Biol Chem. 2009;284(16):10433–45. doi:https://doi.org/10.1074/jbc.M809801200.PubMedPubMed CentralView ArticleGoogle Scholar
- Johnson JO, Mandrioli J, Benatar M, Abramzon Y, Van Deerlin VM, Trojanowski JQ, et al. Exome sequencing reveals VCP mutations as a cause of familial ALS. Neuron. 2010;68(5):857–64. doi:https://doi.org/10.1016/j.neuron.2010.11.036.PubMedPubMed CentralView ArticleGoogle Scholar
- Watts GD, Wymer J, Kovach MJ, Mehta SG, Mumm S, Darvish D, et al. Inclusion body myopathy associated with Paget disease of bone and frontotemporal dementia is caused by mutant valosin-containing protein. Nat Genet. 2004;36(4):377–81. doi:https://doi.org/10.1038/ng1332.PubMedView ArticleGoogle Scholar