Association of the CACNA2D2 gene with schizophrenia in Chinese Han population

Yingli Fu; Na Zhou; Wei Bai; Yaoyao Sun; Xin Chen; Yueying Wang; Mingyuan Zhang; Changgui Kou; Yaqin Yu; Qiong Yu

doi:10.7717/peerj.8521

Association of the CACNA2D2 gene with schizophrenia in Chinese Han population

Yingli Fu^1,2, Na Zhou³, Wei Bai¹, Yaoyao Sun¹, Xin Chen¹, Yueying Wang¹, Mingyuan Zhang¹, Changgui Kou¹, Yaqin Yu¹, Qiong Yu ¹

1Department of Epidemiology and Statistics, School of Public Health, Jilin University, Changchun, Jilin, China

2Division of Clinical Research, First Hospital of Jilin University, Changchun, Jilin, China

3Department of Pharmacy, School and Hospital of Stomatology, Jilin University, Changchun, Jilin, China

DOI: 10.7717/peerj.8521

Published: 2020-02-05
Accepted: 2020-01-07
Received: 2019-07-24

Academic Editor: Bao-Liang Zhong

Subject Areas: Genetics, Epidemiology
Keywords: Schizophrenia, CACNA2D2, SNPS, Haplotype

Copyright: © 2020 Fu et al.
Licence: This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ) and either DOI or URL of the article must be cited.

Cite this article: Fu Y, Zhou N, Bai W, Sun Y, Chen X, Wang Y, Zhang M, Kou C, Yu Y, Yu Q. 2020. Association of the CACNA2D2 gene with schizophrenia in Chinese Han population. PeerJ 8:e8521 https://doi.org/10.7717/peerj.8521

The authors have chosen to make the review history of this article public.

Abstract

Background

Schizophrenia (SCZ) is a severely complex psychiatric disorder in which ~80% can be explained by genetic factors. Single nucleotide polymorphisms (SNPs) in calcium channel genes are potential genetic risk factors for a spectrum of psychiatric disorders including SCZ. This study evaluated the association between SNPs in the voltage-gated calcium channel auxiliary subunit alpha2delta 2 gene (CACNA2D2) and SCZ in the Han Chinese population of Northeast China.

Methods

A total of 761 SCZ patients and 775 healthy controls were involved in this case-control study. Three SNPs (rs3806706, rs45536634 and rs12496815) of CACNA2D2 were genotyped by the MALDI-TOF-MS technology. Genotype distribution and allele frequency differences between cases and controls were tested by Chi-square (χ²) in males and females respectively using SPSS 24.0 software. Linkage disequilibrium and haplotype analyses were conducted using Haploview4.2. The false discovery rate correction was utilized to control for Type I error by R3.2.3.

Results

There was a significant difference in allele frequencies (χ² = 9.545, P_adj = 0.006) and genotype distributions (χ² = 9.275, P_adj = 0.006) of rs45536634 between female SCZ patients and female healthy controls after adjusting for multiple comparisons. Minor allele A (OR = 1.871, 95% CI [1.251–2.798]) and genotype GA + AA (OR = 1.931, 95% CI [1.259–2.963]) were associated with an increased risk of SCZ. Subjects with haplotype AG consisting of rs45536634 and rs12496815 alleles had a higher risk of SCZ (OR = 1.91, 95% CI [1.26–2.90]) compared those with other haplotypes.

Conclusions

This study provides evidence that CACNA2D2 polymorphisms may influence the susceptibility to SCZ in Han Chinese women.

Introduction

Schizophrenia (SCZ) is a severely debilitating psychiatric disorder characterized by positive and negative symptoms as well as cognitive dysfunction (Allen et al., 2008; Koike et al., 2014). The lifetime risk of SCZ is approximately 1% across the world (Mayilyan, Weinberger & Sim, 2008), and the lifetime prevalence of adults in China was 0.6% (Huang et al., 2019). SCZ has a marked impact on the life quality and accounts for approximately 2.8% of the global burden of diseases reported by the World Health Organization in 2001, and the prevalence of SCZ disability was 0.41% in China (Liu et al., 2015).

Genetic and environmental factors may combine to increase SCZ risk (Plomin, Owen & McGuffin, 1994). Genetic factors explain ~80% of the risk for SCZ, and the risk of the SCZ decrease as the parental relationship recedes (Zhu et al., 2009). Despite vigorous genome-wide association studies have been conducted to elucidate the common genetic variations associated with the susceptibility to SCZ (Børglum et al., 2014; Riley et al., 2010), the etiology of SCZ remains obscure, suggesting that additional studies are required to discover the “missing heritability.”

The voltage-gated calcium channel auxiliary subunit alpha2delta 2 gene (CACNA2D2), located in 3p21.31 and highly expressed in the brain, encodes a calcium channel protein. Voltage-gated calcium channels are widely distributed throughout the brain and mediate the intracellular Ca2+ influx of synaptic action potentials (Guan et al., 2016). In recent years, the gene encoding voltage-gated calcium channel subunit attracts wide attention in the field of SCZ pathogenesis. SNPs in calcium channel genes have been identified as genetic risk factors for a spectrum of psychiatric disorders (Cross-Disorder Group of the Psychiatric Genomics Consortium, 2013; Cross-Disorder Group of the Psychiatric Genomics Consortium, 2014). For example, several studies showed that SNPs of CACNA1C were significantly associated with SCZ, and the finding has been confirmed in different populations (Gassó et al., 2016; Zhang et al., 2017; Zhu & Li, 2019). Although there are no reported studies on the relationship between CACNA2D2 and SCZ, CACNA2D2 may cause other psychiatric disorders (Berridge, 2013) and severe neurological diseases (Strupp et al., 2005). In this study, we conducted a genetic association study to examine the association between SNPs of CACNA2D2 and SCZ by a case-control study.

Materials and Methods

Study population

A total of 761 patients with SCZ from the Mental Hospital of Changchun and 775 healthy controls from the physical examination center of the First Hospital of Jilin University were recruited. All participants were Han Chinese. The patients were diagnosed according to the criteria of the International Statistical Classification of Disease and Related Health Problems, Tenth Revision (ICD-10) independently by at least two experienced psychiatrists. Healthy controls matched the patients by gender and age. All controls had no personal or family history of mental illness. This study was performed under protocols approved by the Ethics Committee of Jilin University, China (2014-05-01). All subjects signed written informed consent before participating in this study, and all experiments were performed in accordance with relevant guidelines and regulations.

DNA extraction and SNP selection

Peripheral blood of five mL was collected from each subject, and the genomic DNA was extracted from blood samples using a commercial DNA extraction kit (Kangwei Biotech Company, Beijing, China). SNPs located in the promoter and 3′ untranslated region (UTR) of CACNA2D2 were searched in NCBI-SNP (https://www.ncbi.nlm.nih.gov/snp/) and Ensembl (https://asia.ensembl.org/) databases. We predicted the function of these SNPs in SNPinfo (https://snpinfo.niehs.nih.gov/) and searched for minor allele frequency (MAF) each SNP in 1,000 Genomes. Finally, we chose three SNPs (MAF > 0.05) located in promoter or 3′UTR regions (rs3806706 in the promoter region and rs45536634 and rs12496815 in the 3′ UTR) of CACNA2D2 and predicted to be located in transcription factor binding sites. SNP genotyping was performed using matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF-MS). SNP genotyping reactions were performed in a 384-well Spectro-CHIP using a Mass Array nanodispenser (Sequenom Inc., San Diego, CA, USA). The primers for genotyping were designed by AssayDesigner3.1 and were listed in Table 1.

Table 1:

Primers for polymerase chain reaction.

SNP	Primer sequence (5′–3′)
rs12496815	F: ACGTTGGATGTGGTTTTGGCACCAGTGCGT
	R: ACGTTGGATGTGGCACCCAAATCACATCTC
rs3806706	F: ACGTTGGATGTGAGCTCAACAGCTGCCTTC
	R: ACGTTGGATGGTCCAGCAAACAGGTAAGAG
rs45536634	F: ACGTTGGATGCAATGTATGTCAAGGGCCTG
	R: ACGTTGGATGGAGTCCCACTTAGTGCTCTG

DOI: 10.7717/peerj.8521/table-1

Statistical analysis

Pearson’s Chi-square (χ²) test and Student’s t-test were used to test the distribution of sex and age between case and control groups, respectively. The distributions of allele and genotype were analyzed using χ² tests. The odds ratio (OR) was used to estimate the relative risk of SCZ associated with genotypes with minor alleles. The Type I error due to multiple testing was corrected by the FDR method. All the above analyses were performed using Software SPSS 24.0 (IBM SPSS, IBM Corp, Armonk, NY, USA), and R version 3.2.3 was used for FDR corrections. The Hardy–Weinberg equilibrium (HWE) test was conducted in the case and control group separately by Goodness of fit χ² test using online software SNPStats (https://www.snpstats.net/snpstats/start.htm). The linkage disequilibrium (LD) between SNPs was estimated in both females and males separately using Haploview 4.2 (Barrett et al., 2005), and the haplotype analysis was further performed using Haploview. The statistical power for each SNP was calculated according to the MAF of each SNP (rs45536634: 0.073, rs3806706: 0.378 and rs12496815: 0.388). The prevalence of SCZ (1%) was estimated by Quanto 1.2.4 (Gauderman, 2002). The OR was set from 1.4 to 2.0. All tests were two-sided and P_adj-value less than 0.05 was considered to be statistically significant.

Results

Demographic characteristics

A total of 1,536 subjects were included in this study, comprised of 761 SCZ patients (58.2% males, Mean age = 34.61 ± 12.02 years) and 775 healthy controls (56.2% males, Mean age = 34.74 ± 11.41 years). There was no significant difference either in sex (χ² = 0.681, P = 0.409) or age (t = 0.221, P = 0.825) between patients and healthy controls. The results of the HWE test were shown in Table 2. All SNPs were in accordance with the HWE in both cases and controls (P > 0.05).

Table 2:

Test of HWE for case and control groups.

Gene	SNP	Case				Control
Gene	SNP	Ho	He	χ²	P	Ho	He	χ²	P
CACNA2D2	rs45536634	0.172	0.171	0.024	0.9	0.127	0.137	2.42	0.12
	rs12496815	0.469	0.497	2.126	0.1	0.498	0.499	0.002	0.968
	rs3806706	0.458	0.437	1.66	0.2	0.403	0.422	1.583	0.208

DOI: 10.7717/peerj.8521/table-2

Note:

Ho, observed heterozygosity; He, expected heterozygosity.

The distribution of alleles and genotypes in males and females

The detection rate of rs45536634, rs12496815 and rs3806706 were 97%, 92% and 98%, respectively. Genotype and allele frequencies of females were depicted in Table 3. A significant difference (P_ajd = 0.012) was observed in allele frequencies of rs45536634 between female SCZ patients and female healthy controls. Subjects who carried minor allele A had a 1.9 times higher risk of SCZ than those homozygous for the major G allele. Similarly, a significant difference (P_ajd = 0.006) was observed in the genotype distribution in females, and subjects with the minor allele (GA + AA) has an increased risk of SCZ when compared those with genotype GG (OR = 1.931, 95% CI [1.259–2.963]). These associations were found only in females, but not in males. In the male group, there was a difference between cases and controls in the genotype distribution of rs3806706 (P = 0.02); however, after adjusting for multiple testing, the difference was not significant (P_ajd = 0.12) (Table 4).

Table 3:

Genotype and allele distributions of CACNA2D2 SNPs in female.

SNPs	Genotype	Case	Control	χ²	P	P_adj	OR [95% CI]
rs3806706	GG	140	158	0.032	0.858	0.858	1
	GC + CC	165	181				1.184 [0.689–2.035]
	Allele
	G	417	460	0.039	0.843	0.858	1
	C	193	218				0.977 [0.772–1.235]
rs45536634	GG	227	297	9.275	0.002	0.006*	1
	GA + AA	62	42				1.931 [1.259–2.963]*
	Allele			9.545	0.002	0.006*
	G	513	635				1
	A	65	43				1.871 [1.251–2.798]*
rs12496815	GG	62	69	1.421	0.233	0.395	1
	GA + AA	190	268				0.789 [0.534–1.165]
	Allele
	G	247	308	1.268	0.263	0.395	1
	A	257	366				0.876 [0.695–1.103]

DOI: 10.7717/peerj.8521/table-3

Notes:

P_adj represent P corrected by FDR.

* P_adj < 0.05.

Table 4:

Genotype and allele distributions of CACNA2D2 SNPs in male.

SNPs	Genotype	Case	Control	χ²	P	P_adj	OR [95% CI]
rs3806706	GG	190	225	5.373	0.02	0.12	1
	GC + CC	241	208				1.372 [1.050–1.793]
	Allele
	G	580	617	3.185	0.074	0.172	1
	C	282	249				1.205 [0.982–1.478]
rs45536634	GG	361	369	0.127	0.721	0.865	1
	GA + AA	67	64				1.070 [0.738–1.552]
	Allele			0.00025	0.987	0.987
	G	786	795				1
	A	70	71				1.997 [0.707–1.407]
rs12496815	GG	85	105	0.66	0.416	0.624	1
	GA + AA	303	327				1.145 [0.826–1.586]
	Allele
	G	347	423	2.953	0.086	0.172	1
	A	429	441				1.186 [0.976–1.440]

DOI: 10.7717/peerj.8521/table-4

Note:

P_adj represent P corrected by FDR.

LD and haplotype analysis

As shown in Fig. 1, the LD analysis of rs45536634 and rs12496815 in CACNA2D2 showed that the D’ values were equal to one in both female and male groups. According to the results of LD analysis, haplotype association analyses of rs45536634–rs12496815 were conducted in females and males respectively, and the results are shown in Table 5. Three common haplotypes were estimated to have a frequency > 1%, and haplotype AG was significantly associated with SCZ (OR = 1.91, P_adj = 0.0096) in females.

Linkage disequilibrium (LD) of SNPs within CACNA2D2 in female (A) and male (B) D’ values were used toestimate the LD between pairwise SNPs. — Figure 1: Linkage disequilibrium (LD) of SNPs within *CACNA2D2* in female (A) and male (B) D’ values were used toestimate the LD between pairwise SNPs.

Download full-size image

DOI: 10.7717/peerj.8521/fig-1

Table 5:

Association between haplotypes and schizophrenia by sex.

Haplotype	Male (frequency)					Female (frequency)
Haplotype	Control	Case	OR [95% CI]	P	P_adj	Control	Case	OR [95% CI]	P	P_adj
rs45536634–rs12496815
GA	0.5105	0.5518	1	–	–	0.5424	0.5069	1	–	–
GG	0.4075	0.3666	0.84 [0.68–1.03]	0.086	0.28	0.3942	0.3807	1.02 [0.80–1.31]	0.87	0.87
AG	0.082	0.0816	0.93 [0.66–1.31]	0.66	0.7	0.0634	0.1125	1.91 [1.26–2.90]	0.0026	0.026*
AA	0	0		–	–	0	0		–	–
rs124996815–rs3806706
GG	0.4462	0.3933	1			0.4098	0.4196	1	–	–
AG	0.2659	0.2796	1.18 [0.92–1.51]	0.18	0.3	0.2686	0.264	0.96 [0.71–1.29]	0.77	0.87
AC	0.2439	0.2752	1.27 [1.00–1.61]	0.046	0.23	0.2742	0.2473	0.88 [0.66–1.16]	0.37	0.74
GC	0.044	0.0519	1.33 [0.77–2.29]	0.31	0.388	0.0473	0.0691	1.43 [0.77–2.62]	0.26	0.65
rs124996815–rs3806706–rs45536634
GGG	0.3712	0.3148	1	–	–	0.3492	0.3297	1	–	–
AGG	0.266	0.2798	1.23 [0.95–1.59]	0.12	0.28	0.2683	0.2598	1.04 [0.76–1.43]	0.79	0.87
ACG	0.2439	0.2735	1.31 [1.03–1.69]	0.031	0.23	0.2741	0.2496	0.97 [0.72–1.31]	0.85	0.87
GGA	0.0749	0.0782	1.23 [0.83–1.81]	0.3	0.288	0.061	0.0941	1.67 [1.02–2.72]	0.042	0.167
GCG	0.0369	0.0502	1.58 [0.87–2.87]	0.14	0.28	0.045	0.0482	1.18 [0.60–2.34]	0.63	0.87
rare	0.0071	0.0035	0.72 [0.13–3.85]	0.7	0.7	0.0024	0.0186	8.05 [1.01–64.35]	0.05	0.167

DOI: 10.7717/peerj.8521/table-5

Note:

* P_adj < 0.05.

Statistical power

The statistical power for rs45536634, rs12496815 and rs3806706 were 0.675–0.999, 0.872–0.999 and 0.878–0.999, respectively, if the OR varied from 1.4 to 2.0.

Discussion

The association between variants of a number of genes and SCZ has been reported in previous studies. To the best of our knowledge, this is the first report of a significant association between rs45536634 of CACNA2D2 and SCZ in females of the Northeast Han Chinese population.

It is known that Ca2+ ion represents one of the most important second messengers in the brain and plays an essential role in neuronal development, synaptic transmission and plasticity, besides regulating various metabolic pathways (Striessnig et al., 2006). Notably, as demonstrated by several studies, the Ca2+ homeostasis disorder is associated with many pathological mechanisms, especially those related to neurodegenerative disorders, such as SCZ, Alzheimer’s disease, and bipolar disorder (Berridge, 2013; Sulzer & Surmeier, 2013). CACNA2D2 encodes the Alpha 2 delta 2 subunit of voltage-gated calcium channel (Tedeschi et al., 2016) which is a key signaling element, allowing changes in membrane potential to control a large number of Ca2+ dependent neurotransmitter release and neuronal plasticity in electrically excitable cells (Striessnig et al., 2006). A study conducted by Villela et al. (2016) showed that the copy number change in CACNA2D2 was a risk factor for Alzheimer’s disease.

Moreover, CACNA1C, a gene in the same family as CACNA2D2, has been repeatedly confirmed as one of the susceptibility genes for SCZ in various populations (Guan et al., 2014; He et al., 2014). In addition, Zhang et al. (2018) conducted a review research on calcium channel genes associated with SCZ in the Han Chinese population and found that CACNA1C, CACNB2, CACNA2D1 and CACNA2D3 were related to SCZ.

The sex-specific molecular phenotype of SCZ was observed in previous studies. A study conducted by Oumaima et al. indicated that minor alleles of SNPs in genes LTA and TNFA were over-represented in male SCZ patients but not in female schizophrenia patients (Inoubli et al., 2018). Jemli et al. (2017) researched the association between the functional polymorphism of IFNGR2 with SCZ and found the IFNGR2 Q64R polymorphism was associated with SCZ in males. The study conducted by Yang et al. (2016) showed that the genotypes and allele distributions of rs3087494 in PLA2G12A were significantly associated with SCZ in males, but not in females. Another study focused on sex-specific molecular phenotypes found that eight genes showed a differential expression in female and male schizophrenia patients (Ramsey et al., 2013).

For a better understanding of the association between SCZ and CACNA2D2, a more in-depth investigation—haplotype analysis was carried out to determine whether the combination of specific alleles was associated with the SCZ risk. The AG haplotype, consisted of rs44536634 and rs12496815 alleles, was correlated with an increased risk of SCZ in Han Chinese women. The haplotype analysis not only confirmed the association between rs44536634 and SCZ, but also supported that rs44536634 allele A was associated with an increased risk of SCZ in females. Furthermore, our research provided an evidence to support the distinct molecular phenotypes of SCZ patients with different gender, as reported in previous studies (Ben Nejma et al., 2013; Jemli et al., 2017; Ramsey et al., 2013).

In this study, several limitations should be considered. Firstly, our study was performed at a single center and only three SNPs were analyzed, it ignored SNPs in other genes that may be associated with SCZ. Secondly, the study is limited to interpreting the causal relationship between genetic risk factors and SCZ as this is a cross-sectional study. Furthermore, the representative of this study was limited to the adults of Northeast China because the samples were collected from Jilin Province. Finally, we lost some demographic covariates of the controls when analyzing the association between CACNA2D2 SNP polymorphism and SCZ patients due to the difficulty of demographic characteristics collection. Large-scale examination with more demographic characteristics is warranted to further examine the association between CACNA2D2 and SCZ.

Conclusion

The sample size of the present study was sufficient for detecting the effect of CACNA2D2 variants on SCZ. This study demonstrated that CACNA2D2 polymorphisms might influence the susceptibility to SCZ in Han Chinese women from Northeast China. The findings support the hypothesis that CACNA2D2 may represent a novel susceptibility gene for SCZ in females. Functional genomics studies should be performed in future to validate the function of SCZ-associated CACNA2D2 variants.

Supplemental Information

Raw data.

DOI: 10.7717/peerj.8521/supp-1

Download

[1] Allen NC, Bagade S, McQueen MB, Ioannidis JP, Kavvoura FK, Khoury MJ, Tanzi RE, Bertram L. 2008. Systematic meta-analyses and field synopsis of genetic association studies in schizophrenia: the SzGene database. Nature Genetics 40(7):827-834

[2] Barrett JC, Fry B, Maller J, Daly MJ. 2005. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21(2):263-265

[3] Ben Nejma M, Zaabar I, Zaafrane F, Thabet S, Mechri A, Gaha L, Ben Salem K, Bel Hadj Jrad B. 2013. A gender-specific association of interleukin 1 receptor antagonist polymorphism with schizophrenia susceptibility. Acta Neuropsychiatrica 25(6):349-355

[4] Berridge MJ. 2013. Dysregulation of neural calcium signaling in Alzheimer disease, bipolar disorder and schizophrenia. Prion 7(1):2-13

[5] Børglum AD, Demontis D, Grove J, Pallesen J, Hollegaard MV, Pedersen CB, Hedemand A, Mattheisen M, Investigators G, Uitterlinden A, Nyegaard M, Ørntoft T, Wiuf C, Didriksen M, Nordentoft M, Nöthen MM, Rietschel M, Ophoff RA, Cichon S, Yolken RH, Hougaard DM, Mortensen PB, Mors O. 2014. Genome-wide study of association and interaction with maternal cytomegalovirus infection suggests new schizophrenia loci. Molecular Psychiatry 19(3):325-333

[6] Cross-Disorder Group of the Psychiatric Genomics Consortium. 2014. Biological insights from 108 schizophrenia-associated genetic loci. Nature 511(7510):421-427

[7] Cross-Disorder Group of the Psychiatric Genomics Consortium. 2013. Identification of risk loci with shared effects on five major psychiatric disorders: a genome-wide analysis. Lancet 381(9875):1371-1379

[8] Gassó P, Sánchez-Gistau V, Mas S, Sugranyes G, Rodríguez N, Boloc D, De la Serna E, Romero S, Moreno D, Moreno C, Díaz-Caneja CM, Lafuente A, Castro-Fornieles J. 2016. Association of CACNA1C and SYNE1 in offspring of patients with psychiatric disorders. Psychiatry Research 245:427-435

[9] Gauderman WJ. 2002. Sample size requirements for association studies of gene–gene interaction. American Journal of Epidemiology 155(5):478-484

[10] Guan F, Zhang T, Liu X, Han W, Lin H, Li L, Chen G, Li T. 2016. Evaluation of voltage-dependent calcium channel γ gene families identified several novel potential susceptible genes to schizophrenia. Scientific Reports 6(1):24914

[11] Guan F, Zhang B, Yan T, Li L, Liu F, Li T, Feng Z, Zhang B, Liu X, Li S. 2014. MIR137 gene and target gene CACNA1C of miR-137 contribute to schizophrenia susceptibility in Han Chinese. Schizophrenia Research 152(1):97-104

[12] He K, An Z, Wang Q, Li T, Li Z, Chen J, Li W, Wang T, Ji J, Feng G, Lin H, Yi Q, Shi Y. 2014. CACNA1C, schizophrenia and major depressive disorder in the Han Chinese population. British Journal of Psychiatry 204(1):36-39

[13] Huang Y, Wang Y, Wang H, Liu Z, Yu X, Yan J, Yu Y, Kou C, Xu X, Lu J, Wang Z, He S, Xu Y, He Y, Li T, Guo W, Tian H, Xu G, Xu X, Ma Y, Wang L, Wang L, Yan Y, Wang B, Xiao S, Zhou L, Li L, Tan L, Zhang T, Ma C, Li Q, Ding H, Geng H, Jia F, Shi J, Wang S, Zhang N, Du X, Du X, Wu Y. 2019. Prevalence of mental disorders in China: a cross-sectional epidemiological study. Lancet Psychiatry 6(3):211-224

[14] Inoubli O, Jemli A, Ben Fredj S, Mechri A, Gaha L, Bel Hadj Jrad B. 2018. Haplotypes of TNF α/β genes associated with sex-specific paranoid schizophrenic risk in tunisian population. Disease Markers 2018(1):3502564

[15] Jemli A, Inoubli O, Trifa F, Mechri A, Zaafrane F, Gaha L, Jrad BBH. 2017. IFNGR2 genetic polymorphism associated with sex-specific paranoid schizophrenia risk. Nordic Journal of Psychiatry 71(1):42-47

[16] Koike S, Bundo M, Iwamoto K, Suga M, Kuwabara H, Ohashi Y, Shinoda K, Takano Y, Iwashiro N, Satomura Y, Nagai T, Natsubori T, Tada M, Yamasue H, Kasai K. 2014. A snapshot of plasma metabolites in first-episode schizophrenia: a capillary electrophoresis time-of-flight mass spectrometry study. Translational Psychiatry 4(4):e379

[17] Liu T, Zhang L, Pang L, Li N, Chen G, Zheng X. 2015. Schizophrenia-related disability in China: prevalence, gender, and geographic location. Psychiatric Services 66(3):249-257

[18] Mayilyan KR, Weinberger DR, Sim RB. 2008. The complement system in schizophrenia. Drug News & Perspectives 21(4):200-210

[19] Plomin R, Owen MJ, McGuffin P. 1994. The genetic basis of complex human behaviors. Science 264(5166):1733-1739

[20] Ramsey JM, Schwarz E, Guest PC, Van Beveren NJ, Leweke FM, Rothermundt M, Bogerts B, Steiner J, Bahn S. 2013. Distinct molecular phenotypes in male and female schizophrenia patients. PLOS ONE 8(11):e78729

[21] Riley B, Thiselton D, Maher BS, Bigdeli T, Wormley B, McMichael GO, Fanous AH, Vladimirov V, O’Neill FA, Walsh D, Kendler KS. 2010. Replication of association between schizophrenia and ZNF804A in the Irish case-control study of Schizophrenia sample. Molecular Psychiatry 15(1):29-37

[22] Striessnig J, Koschak A, Sinnegger-Brauns MJ, Hetzenauer A, Nguyen NK, Busquet P, Pelster G, Singewald N. 2006. Role of voltage-gated L-type Ca²⁺ channel isoforms for brain function. Biochemical Society Transactions 34(5):903-909

[23] Strupp M, Kalla R, Freilinger T, Dichgans M, Brandt T. 2005. Dysfunction of the brain calcium channel CaV2.1 in absence epilepsy and episodic ataxia—a comment. Brain 128(6):E32

[24] Sulzer D, Surmeier DJ. 2013. Neuronal vulnerability, pathogenesis, and Parkinson’s disease. Movement Disorders 28(6):715-724

[25] Tedeschi A, Dupraz S, Laskowski CJ, Xue J, Ulas T, Beyer M, Schultze JL, Bradke F. 2016. The calcium channel subunit Alpha2delta2 suppresses axon regeneration in the adult CNS. Neuron 92(2):419-434

[26] Villela D, Suemoto CK, Pasqualucci CA, Grinberg LT, Rosenberg C. 2016. Do copy number changes in CACNA2D2, CACNA2D3, and CACNA1D constitute a predisposing risk factor for Alzheimer’s disease? Frontiers in Genetics 7(e74203):107

[27] Yang G, Xu H, Zhang H, Yu Q, Wu Y, Shi J, Rao W, You Y, Kou C. 2016. Association between PLA2G12A Polymorphisms and Schizophrenia in a Han Chinese population from Northeast China. PLOS ONE 11(7):e0159584