Further evidence for the genetic association between CACNA1I and schizophrenia
© The Author(s) 2017
Received: 12 November 2017
Accepted: 17 December 2017
Published: 2 January 2018
Recent large-scale genome-wide association studies (GWAS) have showed that the neuronal calcium signaling has pivotal roles in schizophrenia (SCZ) in populations of European of ancestry. However, it is not known if calcium signaling pathway genes are also associated with SCZ in Han Chinese population.
Here we investigated the association between genetic variants in three calcium signaling pathway genes (CACNB2, CACNA1C and CACNA1I) and SCZ in 1615 SCZ cases and 1597 controls.
A single nucleotide polymorphism (SNP) (rs4522708) in CACNA1I is significantly associated with SCZ in our Chinese sample (ORA allele = 1.19, corrected P = 0.042), suggesting that CACNA1I may also be a risk gene for SCZ in Chinese population. Of note, the risk allele (A allele) of SNP rs4522708 is same in European and Chinese populations. Meta-analysis of Chinese and European samples further strengthened the association of rs4522708 with SCZ (ORA allele = 1.074, P = 6.26 × 10−11). Expression analysis showed that CACNA1I was significantly up-regulated in hippocampus of SCZ cases compared with controls, implying that dysregulation of CACNA1I may have a role in schizophrenia pathogenesis.
Our study suggests that CACNA1I is a risk gene for SCZ in Chinese population and provides further evidence that supports the potential role of neuronal calcium signaling in schizophrenia.
Schizophrenia (SCZ) is a severe mental disorder that affects about 1% of the world’s population . The core symptoms of SCZ include delusions, hallucinations and cognitive dysfunction [2, 3]. Adoption, family and twin studies indicated that SCZ has a strong genetic component . The heritability of SCZ was estimated about 80% [4, 5], suggesting that genetic factor plays a crucial role in SCZ. Though numerous genetic linkage and association studies have been carried out in different ethnic populations and multiple susceptibility genes have been identified [6–8], only limited risk genes have been consistently replicated in diverse populations. The advent of genome-wide association studies (GWAS) provides an opportunity to dissect the genetic basis of SCZ. In the past decade, multiple GWAS have been performed in world’s populations and multiple promising risk variants and genes have been reported [9–13]. In 2014, the schizophrenia working group of the psychiatric genomics consortium (PGC2) conducted a large-scale meta-analysis of GWAS through using 36,989 SCZ cases and 113,075 controls . Over 100 SCZ-associated loci have been identified by PGC. A detailed characterization of the risk loci showed that calcium signaling pathway genes may play pivotal roles in SCZ . Two calcium signaling pathway genes (including CACNA1C and CACNA1I) showed robust association with SCZ in PGC2, strongly implying these genes may be involved in SCZ pathogenesis. In addition, another calcium signaling pathway gene (CACNB2) was also reported to be associated with SCZ . A recent trans-ancestry meta-analyses of Chinese and PGC2 samples (a total of 43,175 cases and 65,166 controls) further supports the association between CACNB2, CACNA1C, CACNA1I genes and schizophrenia . These findings strongly suggest that neuronal calcium signaling pathway is involved in schizophrenia. To further investigate if these calcium signaling pathway genes are also associated with SCZ in Han Chinese population, we performed a replication study through recruiting 1615 SCZ cases and 1597 controls. We first genotyped three SNPs (one SNP in CACNB2, one SNP in CACNA1C and one SNP in CACNA1I) and tested their associations with SCZ in a Chinese sample. We then explored the expression level of these three genes in schizophrenia cases and controls. Our study indicates that a SNP in CACNA1I is significantly associated with SCZ in Chinese population, suggesting that CACNA1I is also a risk gene for SCZ in Han Chinese.
Materials and methods
We included 1615 SCZ cases and 1597 controls in this study. All of SCZ cases were from local mental health hospitals (or centers) and some of the patients have been reported in our previously studies [16, 17]. Diagnosis was finished by at least two experienced psychiatrists with DSM IV criteria. 61.3% of the SCZ cases were females and the ages of the SCZ cases range from 13 to 81 years (35.1 ± 13.26 years). Patients with history of head injuries and drug abuse were excluded in this study. Healthy controls were local residents and screened for lifetime history of mental illness. 63.9% of the controls were females and the ages of the controls range from 17 to 74 years (38.08 ± 12.31 years). Written informed consents were obtained from all the studied subjects and this the research was approved by the internal review board of Kunming Institute of Zoology, CAS. More detailed information about the sample description can be found in our previous studies [16, 17].
As the goal of this study is to explore if the calcium signaling pathway genes reported in European populations are also associated with SCZ in Han Chinese, we selected three previously reported risk genes for SCZ, including CACNB2 , CACNA1C and CACNA1I . All of these three genes encode voltage-dependent calcium channel subunits. CACNB2 encodes voltage-dependent L-type calcium channel subunit beta-2. CACNA1C encodes Cav1.2, an alpha-1 subunit of a voltage-dependent calcium channel. And CACNA1I encodes Cav3.3, an alpha-1I subunit of a voltage-dependent calcium channel. For CACNB2, we selected SNP rs4748478, which showed significant association with SCZ in PGC2 (P = 8.0 × 10−3). For CACNA1C, we selected SNP rs1006737, which frequently reported to be associated with SCZ in previously studies [18–22]. For CACNA1I, we selected SNP rs4522708, which is also significantly associated with SCZ in PGC2. Of note, rs1006737 and rs4522708 reached genome-wide significance level in PGC2 (P = 1.09 × 10−16 and P = 2.41 × 10−16, respectively) .
SNP genotyping was performed using SNaPShot method as described previously . Single nucleotide extension strategy was used for genotyping. Briefly, DNA fragments containing the SNP site were amplified first using PCR. The PCR products were then treated with Shrimp Alkaline Phosphatase (SAP) and Exonuclease I (ExoI), which remove primers and unincorporated dNTPs. The cleaned PCR products were then used as templates for single nucleotide extension. The genotyping primers were designed to stop just one base upstream of the interest SNP. As ddNTPs were used, so the polymerase extends the primer by one nucleotide, adding a single ddNTP to its 3′ end. ABI 3730 was used to determine which base was added. More detailed information about SNaPShot genotyping method can be found in our previous study .
We used PLINK (v1.9)  to perform the association test (chi-square allelic test with 1 df). Bonferroni correction was used to correct the P values. Meta-analysis was performed using metafor package (http://www.metafor-project.org/)  implemented in R as described previously  and Fixed-effect model was used.
Expression analysis of CACNA1I in SCZ cases and controls
To explore if expression of these CACNA1I was changed in SCZ cases compared with controls, we examined the expression level of CACNA1I in SCZ subjects and controls using GES53987 . In brief, 19 SCZ cases and 19 controls were included in GSE53987 and gene expression in hippocampus was measured. The raw expression values were processed as described in the original study  and student t test was used to determine if expression of CACNA1I is significantly changed in SCZ cases compared with controls.
Linkage disequilibrium analysis
Linkage disequilibrium (LD) analysis were performed using genotype information from the 1000 genomes project . As all of the subjects included in this study were Han Chinese, we only used the genotype information of 103 Chinese subjects (CHB) from the 1000 genomes project to construct and calculate LD values (r ) among the SNPs. Haploview program was used to calculated LD .
Expression analysis of CACNA1I in human tissues
To explore the expression of CACNA1I in human tissues, we downloaded the RNA-seq-based expression data from the study of Fagerberg et al.  Briefly, expression data from 27 different tissues of 95 human individuals were used.
A SNP in CACNA1I is significantly associated with SCZ in Chinese population
Association significance between the studied SNPs and schizophrenia
SNP rs4522708 has two alleles, A and G. The risk allele of rs4522708 in our sample is A allele (ORA allele = 1.19, P = 0.0145). Of note, we found the risk allele of rs4522708 is also the A allele in PGC2 (ORA allele = 1.071, P = 2.41 × 10−10) (35,476 schizophrenia cases and 46,839 controls). We thus performed a meta-analysis through combing our samples (1615 cases and 1597 controls) with samples from PGC2 (35,476 schizophrenia cases and 46,839 controls). A total of 37,091 SCZ cases and 48,436 controls were included in the meta-analysis. Heterogeneity test showed there was no heterogeneity in the combined samples (P > 0.05). The meta-analysis (fixed-effect model) further strengthened the association between rs4522708 and SCZ (ORA allele = 1.074, two-tailed P = 6.26 × 10−11). Taken together, these results suggest that rs4522708 may represent an authentic risk variant for SCZ.
CACNA1I was significantly up-regulated in SCZ cases compared with controls
CACNA1I was highly expressed in human brains compared with other tissues
Recent large-scale GWAS of schizophrenia (PGC2) have frequently reported that calcium signaling pathway genes were associated with schizophrenia. Of note, genetic variants in CACNA1C and CACNA1I reached genome-wide significance level in PGC2. In addition, previous study also showed that CACNB2 gene was a risk gene for SCZ. Though the associations between these three calcium signaling pathway genes and SCZ have been frequently reported, we noticed that most of the studies were from populations of European ancestry. It is not known if these three calcium signaling pathway genes were also associated with SCZ in Han Chinese. In this study, we genotyped three SNPs in CACNB2, CACNA1C and CACNA1I in 1615 SCZ cases and 1597 controls (All of the subjects were Han Chinese). We found that a SNP in CACNA1I is significantly associated with SCZ in our Chinese sample, suggesting CACNA1I may also represent a risk gene for SCZ in Chinese population. To further explore the potential role of CACNA1I in SCZ, we examined expression analysis and found that expression of CACNA1I was significantly up-regulated in SCZ cases compared with controls, implying dysregulation of CACNA1I may have a role in SCZ.
There were several limitations of this study. First, only one SNP was selected for each of the three calcium signaling pathway genes in this study. Thus, the SNP coverage is limited. We could not exclude if other SNPs in CACNB2 and CACNA1C were associated with SCZ. Second, the sample size is relatively small in this study. Further replication study with larger sample size is needed to validate our results. Taken together, our study showed that CACNA1I is also associated with SCZ in Chinese population. Our results provide further evidence that support the association of CACNA1I with SCZ.
We thank all of the participants in the study. This work was supported by the Key Research Project of Yunnan Province (2017FA008 to X.J.L).
XJL conceived and designed the study. YJX, DH and LW performed the genotyping and analyzed the data. XJL wrote the paper. All authors reviewed the manuscript. All authors read and approved the final manuscript.
Competing of interest
The authors declare no competing financial 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.
- Saha S, Chant D, Welham J, McGrath J. A systematic review of the prevalence of schizophrenia. PLoS Med. 2005;2:e141.View ArticlePubMedPubMed CentralGoogle Scholar
- Andreasen NC. Symptoms, signs, and diagnosis of schizophrenia. Lancet. 1995;346:477–81.View ArticlePubMedGoogle Scholar
- Owen MJ, Sawa A, Mortensen PB. Schizophrenia. Lancet. 2016;388:86–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Sullivan PF, Kendler KS, Neale MC. Schizophrenia as a complex trait: evidence from a meta-analysis of twin studies. Arch Gen Psychiatry. 2003;60:1187–92.View ArticlePubMedGoogle Scholar
- Sullivan PF, Daly MJ, O'Donovan M. Genetic architectures of psychiatric disorders: the emerging picture and its implications. Nat Rev Genet. 2012;13:537–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Lewis CM, Levinson DF, Wise LH, DeLisi LE, Straub RE, Hovatta I, et al. Genome scan meta-analysis of schizophrenia and bipolar disorder, part II: schizophrenia. Am J Hum Genet. 2003;73:34–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Ng MY, Levinson DF, Faraone SV, Suarez BK, DeLisi LE, Arinami T, et al. Meta-analysis of 32 genome-wide linkage studies of schizophrenia. Mol Psychiatry. 2009;14:774–85.View ArticlePubMedGoogle Scholar
- Allen NC, Bagade S, McQueen MB, Ioannidis JP, Kavvoura FK, Khoury MJ, et al. Systematic meta-analyses and field synopsis of genetic association studies in schizophrenia: the SzGene database. Nat Genet. 2008;40:827–34.View ArticlePubMedGoogle Scholar
- O'Donovan MC, Craddock N, Norton N, Williams H, Peirce T, Moskvina V, et al. Identification of loci associated with schizophrenia by genome-wide association and follow-up. Nat Genet. 2008;40:1053–5.View ArticlePubMedGoogle Scholar
- Yue WH, Wang HF, Sun LD, Tang FL, Liu ZH, Zhang HX, et al. Genome-wide association study identifies a susceptibility locus for schizophrenia in Han Chinese at 11p11.2. Nat Genet. 2011;43:1228–31.View ArticlePubMedGoogle Scholar
- Shi Y, Li Z, Xu Q, Wang T, Li T, Shen J, et al. Common variants on 8p12 and 1q24.2 confer risk of schizophrenia. Nat Genet. 2011;43:1224–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Ripke S, O'Dushlaine C, Chambert K, Moran JL, Kahler AK, Akterin S, et al. Genome-wide association analysis identifies 13 new risk loci for schizophrenia. Nat Genet. 2013;45:1150–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Schizophrenia Working Group of the Psychiatric Genomics Consortium. Biological insights from 108 schizophrenia-associated genetic loci. Nature. 2014;511:421–7.View ArticlePubMed CentralGoogle Scholar
- Juraeva D, Haenisch B, Zapatka M, Frank J, Witt SH, Muhleisen TW, et al. Integrated pathway-based approach identifies association between genomic regions at CTCF and CACNB2 and schizophrenia. PLoS Genet. 2014;10:e1004345.View ArticlePubMedPubMed CentralGoogle Scholar
- Li Z, Chen J, Yu H, He L, Xu Y, Zhang D, et al. Genome-wide association analysis identifies 30 new susceptibility loci for schizophrenia. Nat Genet. 2017;49:1576–83.View ArticlePubMedGoogle Scholar
- Luo XJ, Diao HB, Wang JK, Zhang H, Zhao ZM, Su B. Association. Of haplotypes spanning PDZ-GEF2, LOC728637 and ACSL6 with schizophrenia in Han Chinese. J Med Genet. 2008;45:818–26.View ArticlePubMedGoogle Scholar
- Li M, Luo XJ, Xiao X, Shi L, Liu XY, Yin LD, et al. Allelic differences between Han Chinese and Europeans for functional variants in ZNF804A and their association with schizophrenia. Am J Psychiatry. 2011;168:1318–25.View ArticlePubMedGoogle Scholar
- Green EK, Grozeva D, Jones I, Jones L, Kirov G, Caesar S, et al. The bipolar disorder risk allele at CACNA1C also confers risk of recurrent major depression and of schizophrenia. Mol Psychiatry. 2010;15:1016–22.View ArticlePubMedGoogle Scholar
- Nyegaard M, Demontis D, Foldager L, Hedemand A, Flint TJ, Sorensen KM, et al. CACNA1C (rs1006737) is associated with schizophrenia. Mol Psychiatry. 2010;15:119–21.View ArticlePubMedGoogle Scholar
- He K, An Z, Wang Q, Li T, Li Z, Chen J, et al. CACNA1C, schizophrenia and major depressive disorder in the Han Chinese population. Br J Psychiatry. 2014;204:36–9.View ArticlePubMedGoogle Scholar
- Zheng F, Zhang Y, Xie W, Li W, Jin C, Mi W, et al. Further evidence for genetic association of CACNA1C and schizophrenia: new risk loci in a Han Chinese population and a meta-analysis. Schizophr Res. 2014;152:105–10.View ArticlePubMedGoogle Scholar
- Ivorra JL, Rivero O, Costas J, Iniesta R, Arrojo M, Ramos-Rios R, et al. Replication of previous genome-wide association studies of psychiatric diseases in a large schizophrenia case-control sample from Spain. Schizophr Res. 2014;159:107–13.View ArticlePubMedGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Viechtbauer W. Conducting Meta-Analyses in R with the metafor Package. J Statistical Software. 2010;36 https://doi.org/10.18637/jss.v18036.i18603.
- Luo XJ, Li M, Huang L, Steinberg S, Mattheisen M, Liang G, et al. Convergent lines of evidence support CAMKK2 as a schizophrenia susceptibility gene. Mol Psychiatry. 2014;19:774–83.View ArticlePubMedGoogle Scholar
- Lanz TA, Joshi JJ, Reinhart V, Johnson K, Grantham LE, 2nd, Volfson D. STEP Levels are unchanged in pre-frontal cortex and associative striatum in post-mortem human brain samples from subjects with schizophrenia, bipolar disorder and major depressive disorder. PLoS One 2015; 10: e0121744.Google Scholar
- Abecasis GR, Altshuler D, Auton A, Brooks LD, Durbin RM, Gibbs RA, et al. A map of human genome variation from population-scale sequencing. Nature. 2010;467:1061–73.View ArticlePubMedGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–5.View ArticlePubMedGoogle Scholar
- Fagerberg L, Hallstrom BM, Oksvold P, Kampf C, Djureinovic D, Odeberg J, et al. Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. Mol Cell Proteomics. 2014;13:397–406.View ArticlePubMedGoogle Scholar
- Xu W, Liu Y, Chen J, Guo Q, Liu K, Wen Z, et al. Genetic risk between the CACNA1I gene and schizophrenia in Chinese Uygur population. Hereditas. 2017;155:5.View ArticlePubMedPubMed CentralGoogle Scholar