Two Novel SNPs in ATXN3 3’ UTR May Decrease Age at Onset of SCA3/MJD in Chinese Patients

Abstract

Spinocerebellar ataxia type 3 (SCA3), or Machado—Joseph disease (MJD), is an autosomal dominantly-inherited disease that produces progressive problems with movement. It is caused by the expansion of an area of CAG repeats in a coding region of ATXN3. The number of repeats is inversely associated with age at disease onset (AO) and is significantly associated with disease severity; however, the degree of CAG expansion only explains 50 to 70% of variance in AO. We tested two SNPs, rs709930 and rs910369, in the 3’ UTR of ATXN3 gene for association with SCA3/MJD risk and with SCA3/MJD AO in an independent cohort of 170 patients with SCA3/MJD and 200 healthy controls from mainland China. rs709930 genotype frequencies were statistically significantly different between patients and controls (p = 0.001, α = 0.05). SCA3/MJD patients carrying the rs709930 A allele and rs910369 T allele experienced an earlier onset, with a decrease in AO of approximately 2 to 4 years. The two novel SNPs found in this study might be genetic modifiers for AO in SCA3/MJD.

About the Authors:

Zhe Long

¶‡ These authors contributed equally to this work.

Affiliation: Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan 410008, P. R. China

Zhao Chen

¶‡ These authors contributed equally to this work.

Affiliation: Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan 410008, P. R. China

Chunrong Wang

Affiliations Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan 410008, P. R. China, Department of Neurology, Xiangtan Central Hospital, Xiangtan, Hunan 411100, P. R. China

Fengzhen Huang

Affiliations Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan 410008, P. R. China, Department of Neurology & Institute of Translational Medicine at University of South China, the First People’s Hospital of Chenzhou, Chenzhou, Hunan 423000, P. R. China

Huirong Peng

Affiliation: Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan 410008, P. R. China

Xuan Hou

Affiliation: Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan 410008, P. R. China

Dongxue Ding

Affiliation: Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan 410008, P. R. China

Wei Ye

Affiliation: Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan 410008, P. R. China

Junling Wang

Affiliations Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan 410008, P. R. China, Key Laboratory of Hunan Province in Neurodegenerative Disorders, Central South University, Changsha, Hunan 410008, P. R. China, State Key Lab of Medical Genetics, Central South University, Changsha, Hunan 410078, P. R. China

Qian Pan

Affiliation: State Key Lab of Medical Genetics, Central South University, Changsha, Hunan 410078, P. R. China

Jiada Li

Affiliation: State Key Lab of Medical Genetics, Central South University, Changsha, Hunan 410078, P. R. China

Kun Xia

Affiliation: State Key Lab of Medical Genetics, Central South University, Changsha, Hunan 410078, P. R. China

Beisha Tang

Affiliations Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan 410008, P. R. China, Key Laboratory of Hunan Province in Neurodegenerative Disorders, Central South University, Changsha, Hunan 410008, P. R. China, State Key Lab of Medical Genetics, Central South University, Changsha, Hunan 410078, P. R. China

Tetsuo Ashizawa

Affiliation: Department of Neurology, University of Florida, Gainesville, Florida, United States of America

Hong Jiang

* E-mail: [email protected]

Affiliations Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan 410008, P. R. China, Key Laboratory of Hunan Province in Neurodegenerative Disorders, Central South University, Changsha, Hunan 410008, P. R. China, State Key Lab of Medical Genetics, Central South University, Changsha, Hunan 410078, P. R. China

Introduction

Spinocerebellar ataxia type 3 (SCA3), or Machado—Joseph disease (MJD), is an autosomal dominantly-inherited polyglutamine (polyQ) disease. It presents heterogeneous clinical features, such as progressive cerebellar ataxia, external ophthalmoplegia, dysarthria, dysphagia, pyramidal signs, dystonia, rigidity, and distal muscle atrophy [1]. In general, SCA3/MJD is the most common autosomal dominant spinocerebellar ataxia in China, accounting for 62.09% of cases [2]. It is caused by the expansion of CAG repeats within the ATXN3 coding region: normal individuals have 12 to 40 repeats, while patients have 51 to 86. The age at onset (AO) of SCA3/MJD ranges from 4 years to 75 years; mean AO is around 40 years [3, 4]. In patients, the number of CAG repeats is inversely associated with AO of SCA3/MJD and significantly affect disease severity [5, 6]. However, only 50 to 70% of variance in AO can be ascribed to the degree of CAG expansion, suggesting other factors may modulate AO [7, 8].

In a previous study, we used the miRCURY LNA Array (Exiqon, Denmark) to assay miRNA expression levels in serums from SCA3/MJD patients, finding that the miR-25 expression levels were significantly different than those of healthy controls [9]. Subsequently, in vitro, we confirmed that miR-25 directly targets sequences at position 259–266 of the ATXN3 3’UTR and down-regulates the aggregation of polyQ-expanded ataxin-3 protein (officially accepted by FEBS Lett, 2014). Accordingly, we used 170 SCA3/MJD patients and 200 healthy controls to look for mirSNPs in the miR-25 gene or in miR-25 target sequences in the 3’ UTR of the ATXN3 gene related to the results above. We did not find any miR-25 related mirSNPs. However, we did find two SNPs, rs910369 and rs709930, located 107 base pairs upstream and 108 base pairs downstream from the miR-25 target sequences in the 3’ UTR of the ATXN3 gene respectively: the two SNPs were in complete linkage disequilibrium (D’ = 1.000), as identified by SHEsis (http://analysis.bio-x.cn/myAnalysis.php). We tested the two SNPs for a statistical association with SCA3/MJD, specifically with AO of SCA3/MJD patients.

Methods

Subjects and materials

In our study, we recruited 170 unrelated Chinese SCA3/MJD patients (97 males and 73 females) from the outpatient neurology clinic of Xiangya Hospital, Central South University, in Hunan, China. Mean age of patients was 42.91±11.18 years old (range from 18 to 75 years old). This case-control study included 200 healthy controls matched for age, gender, ethnicity and area of residence. Written informed consent was obtained from all individuals. The study was approved by the Expert Committee of Xiangya Hospital of Central South University (equivalent to an Institutional Review Board).

Clinical data were collected, including AO, clinical presentations, and duration of SCA3/MJD (Table 1). Genomic DNA was extracted from all samples according to a standard protocol [10]. The (CAG)n tract size was determined for all SCA3/MJD patients by T-vector cloning and direct DNA sequencing. The miR-25 target site in the 3’ UTR of ATXN3 were amplified using a pair of primers designed using Primer3 (http://bioinfo.ut.ee/primer3-0.4.0/): 5’-GGCAGCTGTGACCATGTC-3’ (forward) and 5’-AGCATCTGGGAAAGCACATG-3’ (reverse). The amplification reactions contained 1 μL genomic DNA (50 ng/μL), 0.1 μL LA Taq DNA Polymerase (TaKaRa, Japan), 0.4 μL dNTPs, 0.2 μL of each primer, 3.1 μL sterile water and 5 μL 2×GC Buffer I (TaKaRa, Japan), for a total of 10 μL. Genomic DNA was amplified by polymerase chain reaction (PCR) performed in Mastercyclers (Eppendorf AG, 22331 Hamburg, Germany) under the following conditions: initial denaturation at 95.0 C for 5 minutes, followed by 35 cycles of 94.0 C for 30 seconds, 63.5 C for 30 seconds, and 72.0 C for 40 seconds. PCR products were sequenced on an ABI 3730XL DNA Analyzer (Applied Biosystems, Foster City, CA, USA). The reference gDNA sequence was obtained from the UCSC Genome Browser (http://genome.ucsc.edu/). All sequences were compared to the referenced gDNA sequence, using Chromas software (Technelysium Pty Ltd, Brisbane, Australia).

[Figure omitted. See PDF.]

Table 1. Clinical presentations of SCA3/MJD patients.

https://doi.org/10.1371/journal.pone.0117488.t001

Statistics

Given their complete correlation, we tested only one of the two SNPs for association, rs709930. Frequencies of the rs709930 genotypes and alleles in SCA3/MJD patients and healthy controls were calculated, and analyzed for statistically significant differences using the Chi-square test. The difference of the distribution of CAG repeat lengths between the presence and the absence of the rs709930 A allele was analyzed using the Mann-Whitney U test. The risk of developing SCA3/MJD before age 36 years (mean AO for the present series: 36.29±9.81) among patients with the rs709930 A allele was estimated as an odds ratio by logistic regression analysis, with AO≤36 years versus AO＞36 years as the dependent variable. We used multivariate linear regression analysis to estimate the effect of several potential variables on AO, including the (CAG)n tract sizes, the genotypes of the rs709930, and gender. AO for the rs709930 genotypes was adjusted for the mean CAG repeat lengths in the expanded ATXN3 allele after fitting a linear regression model. All analyses were performed using the SPSS Statistics software (version 17.0). A P-value ≤0.05 was considered statistically significant.

Results

Three rs709930 genotypes (G/G, G/A, A/A) appeared both in SCA3/MJD patients and healthy controls. rs709930 genotypes and alleles frequencies were statistically significantly different between patients and controls (both p = 0.001,α = 0.05) (Table 2).

[Figure omitted. See PDF.]

Table 2. Genotype and allele frequencies of rs709930 in SCA3/MJD patients and healthy controls.

https://doi.org/10.1371/journal.pone.0117488.t002

AO of SCA3/MJD was negatively correlated with the CAG repeat lengths in the expanded ATXN3 gene, confirming previous findings [1] (present series, r = －0.694, p = 0.000) (Fig. 1). However, no significant difference was detected in the distribution of the CAG repeat lengths between the presence and the absence of the rs709930 A allele (p = 0.628). A statistically significant difference was found between the AO≤36 years group and the AO＞36 years group, demonstrating the association between the presence of the A allele and an earlier AO (p = 0.001, OR = 2.660, 95% CI = [1.506, 4.700]) that is not mediated by a change in the number of CAG repeats.

[Figure omitted. See PDF.]

Fig 1. The association of expanded CAG repeats in the expanded ATXN3 gene with age at onset (AO) in SCA3/MJD patients.

The X-axis indicates the expanded CAG repeat lengths and the Y-axis denotes AO in years. AO of SCA3/MJD is inversely correlated with the length of CAG repeat (r = －0.694, p = 0.000).

https://doi.org/10.1371/journal.pone.0117488.g001

Given the relatively earlier onset observed for carriers of the rs709930 A allele, the role of presence or absence of the rs709930 A allele was analyzed. When rs709930 A allele status was taken into consideration, the percentage of onset variance explained increased from 44.6% to 46.5% (F = 74.482, p = 0.000). In this series of patients, after adjusting for the size of expanded CAG repeats, the presence of rs709930 A allele decreased the AO by approximately 2–4 years (Table 3). When gender was considered as a variable, the model was not significantly improved.

[Figure omitted. See PDF.]

Table 3. Age at onset in SCA3/MJD patients with rs709930.

https://doi.org/10.1371/journal.pone.0117488.t003

Discussion

In this study, we found statistically significant differences between the frequencies of rs709930 alleles and genotypes SCA3/MJD patients and healthy controls. We also identified an association between AO of SCA3/MJD and rs709930 (and rs910369) genotype. When the status of rs709930 A allele (and rs910369 T allele) was taken into account, the percentage of total variance explained increased by nearly 2%. However, in our study, the expanded CAG repeats that still came to less than 50% of AO variance, which is lower than previously reported [7, 8]; differences in ancestry and sample size may be responsible for the discrepancy. In addition to the CAG repeat lengths in expanded alleles, AO of SCA3/MJD was significantly earlier in patients with the rs709930 A/A and G/A genotypes (rs910369 T/T and G/T genotype) than those with the rs709930 G/G genotype (rs910369 G/G genotype). Moreover, the presence of rs709930 A allele (and rs910369 T allele) decreased the AO by nearly 2–4 years, which showed that variance of AO might be affected by SNPs of disease-causing genes or disease-associated genes as modifiers.

In addition to variation in the number of CAG repeats in causative genes, the variability of AO in PolyQ diseases can also be explained by the effects of modifier sister genes: ATXN2, ATN1 and HTT in SCA3/MJD; ATXN1 and ATXN3 in SCA6; and ATXN3 and TBP in SCA7 [7], as well as other modifiers like APOE in SCA3/MJD [11, 12]. Additionally, investigation on variation in the 5’ UTR region of ATXN3 was performed to evaluate the potential to influence SCA3/MJD phenotype, however, no improvement on the explanation of AO variance was observed [13]. Here, we firstly found that two novel SNPs in 3’UTR of ATXN3, so called “modifying SNPs”, might have effect on AO of SCA3/MJD, which may provide new insights to explore other modifying factors in PolyQ diseases. Due to a relatively small sample size, these findings should be validated based on a larger dataset, together with functional assays to elucidate the underlying mechanisms. The upcoming investigation on other non-coding region of ATXN3 and other causative genes in PolyQ diseases to explore more “modifying SNPs” may be helpful to expand the spectrum of modifiers in PolyQ diseases, which would be beneficial to better understanding of the relationship among causative genes, genetic modifiers, and phenotypes.

Acknowledgments

The authors are grateful to all subjects for their participation in our study.

Author Contributions

Conceived and designed the experiments: ZL ZC HJ. Performed the experiments: ZL. Analyzed the data: ZL ZC CRW HRP. Contributed reagents/materials/analysis tools: CRW XH FZH HRP DXD WY JLW QP JDL KX BST TA. Wrote the paper: ZL ZC.

Citation: Long Z, Chen Z, Wang C, Huang F, Peng H, Hou X, et al. (2015) Two Novel SNPs in ATXN3 3’ UTR May Decrease Age at Onset of SCA3/MJD in Chinese Patients. PLoS ONE 10(2): e0117488. https://doi.org/10.1371/journal.pone.0117488

References

1. Costa MC, Paulson HL (2012) Toward understanding Machado-Joseph disease. Prog Neurobiol 97: 239–257. pmid:22133674

2. Wang J, Shen L, Lei L, Xu Q, Zhou J, et al. (2011) Spinocerebellar ataxias in mainland China: an updated genetic analysis among a large cohort of familial and sporadic cases. Zhong Nan Da Xue Xue Bao Yi Xue Ban 36: 482–489. pmid:21743138

3. Carvalho DR, La Rocque-Ferreira A, Rizzo IM, Imamura EU, Speck-Martins CE (2008) Homozygosity enhances severity in spinocerebellar ataxia type 3. Pediatr Neurol 38: 296–299. pmid:18358414

4. Sequeiros J, Coutinho P (1993) Epidemiology and clinical aspects of Machado-Joseph disease. Adv Neurol 61: 139–153. pmid:8421964

5. Maciel P, Gaspar C, DeStefano AL, Silveira I, Coutinho P, et al. (1995) Correlation between CAG repeat length and clinical features in Machado-Joseph disease. Am J Hum Genet 57: 54–61. pmid:7611296

6. Durr A, Stevanin G, Cancel G, Duyckaerts C, Abbas N, et al. (1996) Spinocerebellar ataxia 3 and Machado-Joseph disease: clinical, molecular, and neuropathological features. Ann Neurol 39: 490–499. pmid:8619527

7. Tezenas DMS, Durr A, Bauer P, Figueroa KP, Ichikawa Y, et al. (2014) Modulation of the age at onset in spinocerebellar ataxia by CAG tracts in various genes. Brain 137(Pt 9):2444–2455. pmid:24972706

8. van de Warrenburg BP, Hendriks H, Durr A, van Zuijlen MC, Stevanin G, et al. (2005) Age at onset variance analysis in spinocerebellar ataxias: a study in a Dutch-French cohort. Ann Neurol 57: 505–512. pmid:15747371

9. Shi Y, Huang F, Tang B, Li J, Wang J, et al. (2014) MicroRNA profiling in the serums of SCA3/MJD patients. Int J Neurosci 124: 97–101. pmid:23879331

10. Miller SA, Dykes DD, Polesky HF (1988) A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res 16: 1215. pmid:3344216

11. Peng H, Wang C, Chen Z, Sun Z, Jiao B, et al. (2014) The APOE epsilon2 allele may decrease the age at onset in patients with spinocerebellar ataxia type 3 or Machado-Joseph disease from the Chinese Han population. Neurobiol Aging 35: 2115–2179.

12. Bettencourt C, Raposo M, Kazachkova N, Cymbron T, Santos C, et al. (2011) The APOE epsilon2 allele increases the risk of earlier age at onset in Machado-Joseph disease. Arch Neurol 68: 1580–1583. pmid:22159055

13. Bettencourt C, Raposo M, Kazachkova N, Santos C, Kay T, et al. (2012) Sequence analysis of 5′ regulatory regions of the Machado-Joseph disease gene (ATXN3). Cerebellum 11: 1045–1050. pmid:22422287