Luwen Zhang1,2,3†
Zhen Li1,2,3†Qing Liu1,2,3Minglong Shao1,2,3Fuping Sun1,2,3Xi Su1,2,3Meng Song1,2,3Yan Zhang1,2,3Minli Ding1
Yanli Lu1Jiewei Liu4Yongfeng Yang1,2,3
Ming Li1,2,4*
Wenqiang Li1,2,3*
Luxian Lv1,2,3,5*- 1Henan Mental Hospital, The Second Affiliated Hospital of Xinxiang Medical University, Xinxiang, China
- 2Henan Key Lab of Biological Psychiatry, Xinxiang Medical University, Xinxiang, China
- 3International Joint Research Laboratory for Psychiatry and Neuroscience of Henan, Xinxiang, China
- 4Key Laboratory of Animal Models and Human Disease Mechanisms of the Chinese Academy of Sciences and Yunnan Province, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China
- 5Henan Province People’s Hospital, Zhengzhou, China
Objectives: Schizophrenia (SZ) is a complex psychiatric disorder with high heritability, and genetic components are thought to be pivotal risk factors for this illness. The glutamate decarboxylase 1 gene (GAD1) was hypothesized to be a candidate risk locus for SZ given its crucial role in the GABAergic neurotransmission system, and previous studies have examined the associations of single nucleotide polymorphisms (SNPs) spanning the GAD1 gene with SZ. However, inconsistent results were obtained. We hence examined the associations between GAD1 SNPs and SZ in two independent case-control samples of Han Chinese ancestry.
Materials and Methods: Two Han Chinese SZ case-control samples, referred as the discovery sample and the replication sample, respectively, were recruited for the current study. The discovery sample comprised of 528 paranoid SZ cases (with age of first onset ≥ 18) and 528 healthy controls; the independent replication sample contained 1,256 early onset SZ cases (with age of first onset < 18) and 2,661 healthy controls. Logistic regression analysis was performed to examine the associations between GAD1 SNPs and SZ.
Results: Ten SNPs covering GAD1 gene were analyzed in the discovery sample, and two SNPs showed nominal associations with SZ (rs2241165, P = 0.0181, OR = 1.261; rs2241164, P = 0.0225, OR = 1.219). SNP rs2241164 was also nominally significant in the independent replication sample (P = 0.0462, OR = 1.110), and the significance became stronger in a subsequent meta-analysis combining both discovery and replication samples (P = 0.00398, OR = 1.138). Nevertheless, such association could not survive multiple corrections, although the effect size of rs2241164 was comparable with other SZ risk loci identified in genome-wide association studies (GWAS) in Han Chinese population. We also examined the associations between GAD1 SNPs and SZ in published datasets of SZ GWAS in East Asians and Europeans, and no significant associations were observed.
Conclusion: We observed weak associations between GAD1 SNPs and risk of SZ in Han Chinese populations. Further analyses in larger Han Chinese samples with more detailed phenotyping are necessary to elucidate the genetic correlation between GAD1 SNPs and SZ.
Introduction
Schizophrenia (SZ) is a severe mental illness that affects approximately 1% of human populations and imposes a high burden on family and society (Huang et al., 2019; McCutcheon et al., 2020). This illness is characterized by positive symptoms (auditory hallucinations, delusions, and disorganized speech), negative symptoms (blunted affect and social withdrawal), and cognitive dysfunctions (Marder and Cannon, 2019). Although the etiology of SZ remains largely unclear, it is believed to be resulted from a combinatorial impact of complex genetic and environmental factors (Gottesman and Shields, 1967). Studies have also found that genetic factor play important roles in the onset and development of SZ (Sullivan et al., 2003; Lichtenstein et al., 2009).
Recently, the Schizophrenia Working Group of the Psychiatric Genomics Consortium (PGC) has performed the SZ genome-wide association studies (GWAS) with continuously accumulating samples of European ancestry (Schizophrenia Working Group of the Psychiatric Genomics Consortium,, 2014; Pardinas et al., 2018), and identified substantial single nucleotide polymorphisms (SNPs) significantly associated with risk of SZ. In the genetically divergent Han Chinese or East Asian populations, GWAS have also reported multiple genome-wide significant loci associated with SZ (Li et al., 2017; Yu et al., 2017; Lam et al., 2019). Overall, these GWAS have provided essential information regarding SZ pathogenesis, and have confirmed some classical pathological hypotheses of the illness, such as the dopamine hypothesis and neurodevelopmental hypothesis.
Glutamic acid decarboxylase (GAD) catalyzes synthesis of the inhibitory neurotransmitter gamma-amino butyric acid (GABA) in the brain (Akbarian and Huang, 2006). GABA is essentially implicated in the pathophysiology of SZ, particularly in three core symptoms: learning, memory, and executive functions (Yoon et al., 2010; Schmidt and Mirnics, 2015). Molecules related to the GABA signaling transduction pathway are widely interacting with other neurotransmitters and signaling pathways implicated in the development of SZ (Represa and Ben-Ari, 2005; Fazzari et al., 2010). To the mammalian, GABA is synthesized by two isoforms of glutamic acid decarboxylase, and GAD1 is the rate-limiting enzyme (Erlander et al., 1991; Szabo et al., 1994). Some autopsy studies have found that both mRNA and protein levels of GAD1 are lower in multiple cortical areas and hippocampus of SZ patients (Akbarian et al., 1995; Fatemi et al., 2005; Hashimoto et al., 2008; Curley et al., 2011; Hyde et al., 2011; Joshi et al., 2012; Ruzicka et al., 2015; Atagun et al., 2017; Tao et al., 2018; Kilonzo et al., 2020). For example, dysfunction of GABAergic neurotransmission has been observed in SZ through investigations of human postmortem brain tissues and in vivo levels of GABA (de Jonge et al., 2017); alternative splicing of the GAD1 gene also likely contributes to GABA dysfunction in the brains of schizophrenics (Zhang et al., 2021). Therefore, GAD1 is likely involved in the pathogenesis of SZ. To date, several previous studies have explored the genetic associations between GAD1 and SZ in humans (Addington et al., 2005; Straub et al., 2007). A study in European populations revealed that the GAD1 gene was associated with childhood-onset SZ and reduction in the volume of gray matter cortex of the cerebrum (Addington et al., 2005). However, few studies have examined the link between GAD1 and SZ in Han Chinese populations.
In this study, using two independent case-control samples recruited from Han Chinese populations, we aimed to systematically investigate the associations of GAD1 SNPs with SZ, and hope to provide evidence for the involvement of this gene in the heritable risk of SZ in Han Chinese subjects.
Materials and Methods
Ethnic Statement and Study Approval
The protocol was approved by the Ethical Committee of the Second Affiliated Hospital of Xinxiang Medical University. Written informed consents were obtained from all participants after the objectives and procedures of the study were fully explained.
Description About Discovery Sample
In the discovery sample, 900 SZ patients were recruited from inpatients of the Second Affiliated Hospital of Xinxiang Medical University, whereas and 900 healthy controls were recruited from employees of the hospital, visitors of the physical examination center and college students of the Xinxiang Medical University. All participants were of Han Chinese origin from northern Henan Province. Diagnoses were made according to the consensus of at least two experienced psychiatrists following the diagnostic criteria for SZ in the Diagnostic and Statistical Manual of Mental Disorders-Fourth Edition IV (DSM-IV) (1994). To ensure inter-rater consistency, all participating psychiatrists were trained every 6 months in which diagnoses and test results were compared using videotaped demonstration interviews. Individuals with complex medical complications history, organic brain diseases, other major psychiatric disorders, or substance dependence were excluded. To attain the aim of achieving a gender-equalized number of patients with paranoid SZ, all female patients with paranoid SZ who did not fulfill the exclusion criteria was included as well as an equal number of male patients with the same diagnosis. The age of first manifestation of positive symptoms was defined as the age at onset of SZ (White et al., 2006). Age at onset was derived from the Comprehensive Assessment of Symptoms and History (CASH) (Andreasen et al., 1992), and all cases had onset of SZ in adulthood (defined as age of first onset ≥ 18). The controls were evaluated by psychiatrists using a simple unstructured interview to exclude those with a personal or family history of mental or neurological diseases. Cases (528 cases) and controls (528 cases) included for the current study were well matched in age, gender (264 males and 264 females pre group, mean age 27.73 ± 8.01 years old, respectively), and ethnicity.
GAD1 SNPs were selected for the present analyses based on the following criteria: (1) all SNPs covering the genomic region of GAD1 gene were subjected to functional analysis using the FASTSNP online service (Yuan et al., 2006), and only SNPs with highly ranked risk and a minor allele frequency (MAF) ≥ 0.05 in the Chinese Beijing population according to the HapMap database were selected (International HapMap Consortium, 2003); (2) tag SNPs were chosen based on the aggressive tagging algorithm (r2 ≥ 0.80, MAF ≥ 0.05) using genotype data from the HapMap dataset as implemented in Haploview v4.1 (Carlson et al., 2004; Barrett et al., 2005); (3) SNPs selected with the above two criteria and had Illumina design scores over 0.6 were proceeded for further analyses.
Peripheral venous blood was drawn from each participant and collected with vacutainer tubes containing anticoagulant ethylene diamine tetracetic acid (EDTA). Genomic DNA was extracted using Blood Genome DNA Extraction Kit (QIAGEN, Hiden, Germany). SNP genotyping was performed using the Illumina Golden Gate assay on a Bead Station 500G Genotyping System. Genomic DNA samples from cases and controls (250 ng/subject) were randomly sorted and genotyped following the Illumina protocol. All assays were performed blind to diagnosis or genotype. Genotype calls were made using the Genotyping module of the Bead Studio 2.0 software (Illumina, Inc., San Diego, United States). Data was examined for cluster separation using quality scores generated by the Illumina software. SNPs with a Gen Train score b0.4 or a cluster separation score b0.6 were considered poorly genotyped and excluded. For quality control analysis, 96 individuals were also genotyped by Sanger DNA sequencing to ensure the quality of Illumina Golden Gate assay. SNPs significantly deviated from Hardy-Weinberg equilibrium (P < 0.05) in the control samples were excluded from further analysis.
Description About Replication Sample
Summary statistics of the SNPs spanning GAD1 from a recent GWAS of early onset SZ (defined as age of first onset < 18) in 1,256 SZ cases and 2,661 healthy controls of Han Chinese origin (Guo et al., 2021) were retrieved for the replication analysis. SZ cases were assessed with Structured Clinical Interview for DSM-IV Axis I Disorders (SCID), and diagnosis was made based on the Diagnostic and Statistical Manual of Mental Disorders IV (DSM-IV) criteria by at least two experienced psychiatrists. The age at first manifestation of positive symptoms was defined as the age at onset of SZ (White et al., 2006). Age at onset was derived from the Comprehensive Assessment of Symptoms and History (CASH) (Andreasen et al., 1992), and all cases were early onset (defined as age of first onset < 18). SZ cases with neurological diseases (including multiple sclerosis and epilepsy), pervasive developmental disorder (including mental retardation), and other psychiatric disorders (including autism, bipolar disorder, depression, and attention deficit hyperactivity disorder), or with a history of drug abuse and head injury were excluded. The mean onset age of the SZ cases was 14.57 ± 2.27 years old, and 41% cases were males. All controls (55% were males) were healthy volunteers at a mean age of 28.60 ± 7.01 years without family history of any psychiatric disorders. There were no overlapped individuals between the discovery and replication samples, and detailed descriptions about the replication samples can be found in the original study (Guo et al., 2021).
The replication sample subjects were genotyped using Illumina Genome-Wide Asian Screening Array (ASA) Chip. Quality control (QC) analyses were conducted following a previously published pipeline (Anderson et al., 2010). After QC, autosomal biallelic SNPs were proceeded for genotype imputation using minimac3 software (Das et al., 2016) with the 1000 Genomes Project Phase 3 data as the reference panel (Genomes Project Consortium Auton et al., 2015). Eventually, autosomal biallelic SNPs with imputation quality score INFO > 0.8, MAF > 0.01, call rate > 95%, and Hardy-Weinberg equilibrium (HWE) P > 1.00 × 10–5 in controls and p > 1.00 × 10–10 in cases were remained. A total of 182 SNPs spanning the coding region of GAD1 gene as well as the ± 20 kb flanking regions were chosen for the present study.
Statistical Analysis
Hardy-Weinberg equilibrium (HWE), linkage disequilibrium (LD) coefficient (r2) and haplotype block were assessed using Haploview v4.1 (Barrett et al., 2005). Logistic regression analysis by adjusting for the significant PCs and sex in the discovery and replication samples was conducted using PLINK (version 1.9) (Purcell et al., 2007) to examine associations of 182 GAD1 SNPs with SZ. SNPs with a P-value of 2.75 × 10–4 (≈0.05/182) was considered statistically significant. For the meta-analysis combining discovery and replication samples, we used odds ratio (OR) and standard error (SE) of each sample to calculate the pooled OR and the overall 95% confidence intervals (CIs) using an inverse-variance-weighted method under either the fixed-effects model or random-effects model. The Higgins and Thompson I2 index was used for assessing heterogeneity across samples (Higgins et al., 2003). SNPs showing pronounced heterogeneity (I2> 75%) underwent meta-analysis in a random-effects model considering the possible test-statistic inflation; SNPs with negligible heterogeneity (I2≤ 75%) underwent analysis in a fixed-effects model (Winkler et al., 2014). Power analysis was performed using the Power and Sample Size Program software, and the observed OR (1.30) of the most significant SNP in a previous Han Chinese SZ GWAS (i.e., according to SNP rs115070292, P = 4.96 × 10–10, OR = 1.30) (Yu et al., 2017) was applied in the power analysis (Dupont and Plummer, 1990). Regional plots were made using LocusZoom1 (Pruim et al., 2010).
Results
We herein performed a two-stage analysis, firstly in the discovery sample and then in an independent replication sample, to explore the association of GAD1 gene with SZ in Han Chinese subjects. The discovery sample included 528 SZ patients and 528 controls of Han Chinese ancestry, and ten common SNPs covering GAD1 were genotyped with a genotyping success rate of 99.82%. We have previously examined the population structure of our discovery sample (Li et al., 2012, 2013) and observed no evidence of population stratification. The replication sample comprised of 1,256 SZ cases and 2,661 healthy controls from Han Chinese populations with negligible population stratification (Guo et al., 2021). A total of 182 SNPs spanning the GAD1 gene locus were included.
We conducted a power analysis of the total sample using the following assumptions:1,784 SZ cases and 3,189 controls, as well as two-tailed P = 2.75 × 10–4. Observed OR of the most significant SNP (i.e., rs115070292, OR = 1.30) from a recent Han Chinese SZ GWAS (Yu et al., 2017) was applied as the input to calculate the statistical power. The present sample size revealed a 78.3% power of detecting a significant allelic association when the allele frequency in general populations reached 0.50.
In the discovery sample (528 SZ cases and 528 controls), the genotypic distributions of ten SNPs were examined. Nine SNPs were in Hardy-Weinberg equilibrium (HWE P > 0.05) in controls, and one SNP rs3791860 with HWE P < 1.00 × 10–5 in both cases and controls was dropped out from further analyses (Supplementary Table 1). LD structures of these SNPs in the SZ cases and controls of our discovery sample were shown in Figure 1. Logistic regression analysis found that two SNPs (rs2241165 and rs2241164) showed nominal associations with SZ (rs2241165, P = 0.0181, OR = 1.261; rs2241164, P = 0.0225, OR = 1.219; Table 1), but neither of them was significant after multiple correction (which is P = 2.75 × 10–4).
Figure 1. The linkage disequilibrium (LD) map of the nine GAD1 SNPs in the SZ cases and controls of our discovery sample. The LD of pairwise SNPs was calculated using r2 algorithm implemented in the Haploview program.
Table 1. Association results of the GAD1 SNPs in our discovery and replication samples.
We then examined the same nine SNPs within GAD1 in the independent replication sample of 1,256 early onset SZ cases and 2,661 controls. We found that rs2241164 was also nominally associated with SZ in the replication sample (P = 0.0462, OR = 1.110), but it did not survive multiple corrections (Table 1). We also retrieved 173 additional SNPs spanning GAD1 to examine their associations with SZ in the replication sample. Although some SNPs exhibited nominal associations with SZ (P < 0.05), none of the SNP-associations could survive multiple correction (Figure 2).
Figure 2. Regional association plots for GAD1 locus in the replication sample. LD information was derived from East Asian individuals in 1000 Genomes Project Phase 3. The LD is defined based on the SNP rs2241164.
To maximize the statistical power of our present sample, we also performed a meta-analysis of the nine SNPs in the total sample combining both discovery and replication samples. While the association between rs2241164 and SZ became stronger after the meta-analysis (P = 0.00398, OR = 1.138, Table 1), it could not survive multiple correction. The other nine SNPs showed no evidence of associations with SZ in the pooled samples (all P > 0.05, Table 1), suggesting that they were unlikely risk factors for SZ at least in Han Chinese populations.
Discussion
A classical etiological hypothesis of SZ is the glutamatergic neuronal dysfunction (Moghaddam and Javitt, 2012), and correlation between abnormal neural transmission in GABAergic circuits and SZ has been confirmed by accumulating studies (Devor et al., 2017). In vertebrates, the GABA synthesis enzyme GAD exists in two isoforms (GAD67 and GAD65) that are encoded by two distinct genes GAD1 and GAD2, respectively (Erlander et al., 1991). GAD67 was found to be the main rate-limiting enzyme of GABA, and down-regulation of the GAD1 gene was hence hypothesized to cause decrease of this inhibitory neurotransmitter and alterations of GABAergic circuits in the brain (Akbarian and Huang, 2006; Brown et al., 2015). Indeed, altered expression of GAD1 has been frequently reported in the SZ brains (Bharadwaj et al., 2013; Mitchell et al., 2015).
We thus sought to analyze the correlation between GAD1 SNPs and SZ in two independent Han Chinese SZ case-control samples. We have observed that certain SNPs spanning GAD1 showed nominal associations with SZ, but they were not statistically significant after multiple corrections either in each individual sample or in the pooled samples. Therefore, GAD1 is unlikely a robust risk gene for SZ in Han Chinese populations. Similarly, non-significant associations between GAD1 SNPs and SZ were reported in a previous Japanese case-control sample (Ikeda et al., 2007). To further verify these results in a larger sample, we also examined SNPs spanning GAD1 in a large GWAS of SZ in East Asians (22,778 cases and 35,362 controls primarily recruited from Han Chinese populations), but no significant SNPs surpassing the P-value threshold of 2.75 × 10–4 (Supplementary Figure 1) were identified (Lam et al., 2019).
We then examined the associations of GAD1 SNPs and SZ susceptibility in European populations. Although several candidate gene studies have reported positive associations between GAD1 and SZ (Addington et al., 2005; Straub et al., 2007), this genetic locus and its SNPs were not highlighted either in Europeans (40,675 cases and 64,643 controls) or in diverse ancestry populations (69,369 cases and 236,642 controls) in two recent large SZ GWAS (Supplementary Figures 2, 3; Pardinas et al., 2018; Ripke et al., 2020). Since these GWAS have undoubtedly sufficient statistical power in the identification of significant risk loci for SZ, we may conclude that the contribution of GAD1 in the genetic susceptibility of SZ in general populations is minimal.
While our analysis did not reveal robust statistical associations between GAD1 SNPs and SZ, it raised a possibility that a gene with putative roles in pathophysiological processes linked with SZ may not always be a strong genetic risk factor (Sullivan, 2007, 2017). In fact, GAD1 is not the only gene supporting this contention. Majority of the genes, which encode proteins with pivotal roles in physiological pathways related to SZ, were extensively analyzed during the earlier era of candidate gene studies in SZ genetics, yet these genes did not show significant associations in later SZ GWAS (Farrell et al., 2015; Johnson et al., 2017). Nevertheless, this inconsistency does not deny potential impact of these genes on SZ. Modulation of their mRNA levels or epigenetic signatures, by undetermined genetic or environmental factors, may still confer risk of SZ. For example, Tao et al. have reported significant differences of GAD1 mRNA expression and epigenetic modifications between SZ cases and controls (Tao et al., 2018).
Notably, inconsistencies between genetic studies in smaller samples and in large GWAS are often observed, i.e., SNPs of certain genes (such GAD1) may be significant in smaller samples, but not significant in large GWAS cohorts. This may be a result of the etiological and phenotypic heterogeneity of the illness. Specifically, such SNPs/genes may participate in the pathogenesis of certain subtype(s) of SZ, whereas these subjects may take up only a limited proportion in large GWAS studies that usually do not stratify samples based on phenotypic characteristics. In the present study, the discovery sample comprised of paranoid SZ cases, and the replication sample included early onset SZ patients. Therefore, the non-significance of the results may actually reflect involvement of GAD1 in certain subtype(s) of SZ remain to be characterized. Future analyses exploring the correlation between GAD1 and various subtypes of SZ in larger Han Chinese samples with detailed phenotyping are thus warranted.
Conclusion
In summary, our study revealed a weak association between GAD1 SNPs and risk of SZ in Han Chinese populations. Further analyses of SNPs in this gene in different certain subtypes of SZ with detailed phenotyping may provide useful insights about the role of this gene in SZ pathogenesis.
Data Availability Statement
The datasets generated for this study are available on https://www.oebiotech.com/index.php?c=show&id=352.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of the Second Affiliated Hospital of Xinxiang Medical University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
WL and LL designed the study protocol. LZ, ZL, QL, MiS, and YY managed the literature searches and analyses. YY, FS, MeS, XS, and YZ conducted sample selection and data management. YZ, MD, YL, and WL undertook the statistical analysis. ML and JL provided advice and discussion during the design and revision of the present study. LZ, WL, ML, and YY wrote the first draft of the manuscript. All authors contributed to and have approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (81971252 to LL; 81903688 to MiS; U1904130 to WL; 82001407 to XS), the National Key Research and Development Program of China (2016YFC1307001 to LL), the Henan Clinical Research Center for Mental Disorders (2019-zxkfkt-004 to YY), the Natural Science Foundation of Henan Province (202300410318 to YZ), the Medical Science and Technology Research Project of Henan Province (LHGJ20200528 to LZ; LHGJ20190482 to QL), the Science and Technology Research Project of Henan Province (212102310589 to YY), the High Scientific and Technological Research Fund of Xinxiang Medical University (2017ZDCG-04 to LL), and the National Undergraduate Training Program for Innovation and Entrepreneurship (202010472026 to LL).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank the patients, their families and the healthy volunteers for their participation, and the physicians who helped us to take clinical data and blood samples in the Second Affiliated Hospital of Xinxiang Medical University.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2021.677153/full#supplementary-material
Footnotes
References
Addington, A. M., Gornick, M., Duckworth, J., Sporn, A., Gogtay, N., Bobb, A., et al. (2005). GAD1 (2q31.1), which encodes glutamic acid decarboxylase (GAD67), is associated with childhood-onset schizophrenia and cortical gray matter volume loss. Mol. Psychiatry 10, 581–588. doi: 10.1038/sj.mp.4001599
Akbarian, S., and Huang, H. S. (2006). Molecular and cellular mechanisms of altered GAD1/GAD67 expression in schizophrenia and related disorders. Brain Res. Rev. 52, 293–304. doi: 10.1016/j.brainresrev.2006.04.001
Akbarian, S., Kim, J. J., Potkin, S. G., Hagman, J. O., Tafazzoli, A., and Bunney, W. E. Jr., et al. (1995). Gene expression for glutamic acid decarboxylase is reduced without loss of neurons in prefrontal cortex of schizophrenics. Arch. Gen. Psychiatry 52, 258–266. doi: 10.1001/archpsyc.1995.03950160008002
Anderson, C. A., Pettersson, F. H., Clarke, G. M., Cardon, L. R., Morris, A. P., and Zondervan, K. T. (2010). Data quality control in genetic case-control association studies. Nat. Protoc. 5, 1564–1573. doi: 10.1038/nprot.2010.116
Andreasen, N. C., Flaum, M., and Arndt, S. (1992). The comprehensive assessment of symptoms and history (CASH). an instrument for assessing diagnosis and psychopathology. Arch. Gen. Psychiatry 49, 615–623. doi: 10.1001/archpsyc.1992.01820080023004
Atagun, M. I., Sikoglu, E. M., Soykan, C., Serdar Suleyman, C., Ulusoy-Kaymak, S., Caykoylu, A., et al. (2017). Perisylvian GABA levels in schizophrenia and bipolar disorder. Neurosci. Lett. 637, 70–74. doi: 10.1016/j.neulet.2016.11.051
Barrett, J. C., Fry, B., Maller, J., and Daly, M. J. (2005). Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21, 263–265. doi: 10.1093/bioinformatics/bth457
Bharadwaj, R., Jiang, Y., Mao, W., Jakovcevski, M., Dincer, A., Krueger, W., et al. (2013). Conserved chromosome 2q31 conformations are associated with transcriptional regulation of GAD1 GABA synthesis enzyme and altered in prefrontal cortex of subjects with schizophrenia. J. Neurosci. 33, 11839–11851. doi: 10.1523/jneurosci.1252-13.2013
Brown, J. A., Ramikie, T. S., Schmidt, M. J., Baldi, R., Garbett, K., Everheart, M. G., et al. (2015). Inhibition of parvalbumin-expressing interneurons results in complex behavioral changes. Mol. Psychiatry 20, 1499–1507. doi: 10.1038/mp.2014.192
Carlson, C. S., Eberle, M. A., Rieder, M. J., Yi, Q., Kruglyak, L., and Nickerson, D. A. (2004). Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am. J. Hum. Genet 74, 106–120. doi: 10.1086/381000
Curley, A. A., Arion, D., Volk, D. W., Asafu-Adjei, J. K., Sampson, A. R., Fish, K. N., et al. (2011). Cortical deficits of glutamic acid decarboxylase 67 expression in schizophrenia: clinical, protein, and cell type-specific features. Am. J. Psychiatry 168, 921–929. doi: 10.1176/appi.ajp.2011.11010052
Das, S., Forer, L., Schonherr, S., Sidore, C., Locke, A. E., Kwong, A., et al. (2016). Next-generation genotype imputation service and methods. Nat. Genet 48, 1284–1287. doi: 10.1038/ng.3656
de Jonge, J. C., Vinkers, C. H., Hulshoff Pol, H. E., and Marsman, A. (2017). GAB aergic mechanisms in schizophrenia: linking postmortem and in vivo studies. Front. Psychiatry 8:118.
Devor, A., Andreassen, O. A., Wang, Y., Maki-Marttunen, T., Smeland, O. B., Fan, C. C., et al. (2017). Genetic evidence for role of integration of fast and slow neurotransmission in schizophrenia. Mol. Psychiatry 22, 792–801. doi: 10.1038/mp.2017.33
Dupont, W. D., and Plummer, W. D. Jr. (1990). Power and sample size calculations. a review and computer program. Control. Clin. Trials 11, 116–128.
Erlander, M. G., Tillakaratne, N. J., Feldblum, S., Patel, N., and Tobin, A. J. (1991). Two genes encode distinct glutamate decarboxylases. Neuron 7, 91–100. doi: 10.1016/0896-6273(91)90077-d
Farrell, M. S., Werge, T., Sklar, P., Owen, M. J., Ophoff, R. A., O’donovan, M. C., et al. (2015). Evaluating historical candidate genes for schizophrenia. Mol. Psychiatry 20, 555–562. doi: 10.1038/mp.2015.16
Fatemi, S. H., Stary, J. M., Earle, J. A., Araghi-Niknam, M., and Eagan, E. (2005). GABAergic dysfunction in schizophrenia and mood disorders as reflected by decreased levels of glutamic acid decarboxylase 65 and 67 kDa and Reelin proteins in cerebellum. Schizophr. Res. 72, 109–122. doi: 10.1016/j.schres.2004.02.017
Fazzari, P., Paternain, A. V., Valiente, M., Pla, R., Lujan, R., Lloyd, K., et al. (2010). Control of cortical GABA circuitry development by Nrg1 and ErbB4 signalling. Nature 464, 1376–1380. doi: 10.1038/nature08928
Genomes Project Consortium, Auton, A., Brooks, L. D., Durbin, R. M., Garrison, E. P., Kang, H. M., et al. (2015). A global reference for human genetic variation. Nature 526, 68–74. doi: 10.1038/nature15393
Gottesman, I., and Shields, J. (1967). A polygenic theory of schizophrenia. Proc. Natl. Acad. Sci. U.S.A. 58, 199–205. doi: 10.1073/pnas.58.1.199
Guo, S., Liu, J., Li, W., Yang, Y., Lv, L., Xiao, X., et al. (2021). Genome wide association study identifies four loci for early onset schizophrenia. Transl. Psychiatry 11:248.
Hashimoto, T., Arion, D., Unger, T., Maldonado-Aviles, J. G., Morris, H. M., Volk, D. W., et al. (2008). Alterations in GABA-related transcriptome in the dorsolateral prefrontal cortex of subjects with schizophrenia. Mol. Psychiatry 13, 147–161. doi: 10.1038/sj.mp.4002011
Higgins, J. P., Thompson, S. G., Deeks, J. J., and Altman, D. G. (2003). Measuring inconsistency in meta-analyses. BMJ 327, 557–560. doi: 10.1136/bmj.327.7414.557
Huang, Y., Wang, Y., Wang, H., Liu, Z., Yu, X., Yan, J., et al. (2019). Prevalence of mental disorders in China: a cross-sectional epidemiological study. Lancet Psychiatry 6, 211–224.
Hyde, T. M., Lipska, B. K., Ali, T., Mathew, S. V., Law, A. J., Metitiri, O. E., et al. (2011). Expression of GABA signaling molecules KCC2, NKCC1, and GAD1 in cortical development and schizophrenia. J. Neurosci. 31, 11088–11095. doi: 10.1523/jneurosci.1234-11.2011
Ikeda, M., Ozaki, N., Yamanouchi, Y., Suzuki, T., Kitajima, T., Kinoshita, Y., et al. (2007). No association between the glutamate decarboxylase 67 gene (GAD1) and schizophrenia in the Japanese population. Schizophr. Res. 91, 22–26. doi: 10.1016/j.schres.2006.12.020
International HapMap Consortium (2003). The international hapmap project. Nature 426, 789–796. doi: 10.1038/nature02168
Johnson, E. C., Border, R., Melroy-Greif, W. E., De Leeuw, C. A., Ehringer, M. A., and Keller, M. C. (2017). No evidence that schizophrenia candidate genes are more associated with schizophrenia than noncandidate genes. Biol. Psychiatry 82, 702–708. doi: 10.1016/j.biopsych.2017.06.033
Joshi, D., Fung, S. J., Rothwell, A., and Weickert, C. S. (2012). Higher gamma-aminobutyric acid neuron density in the white matter of orbital frontal cortex in schizophrenia. Biol. Psychiatry 72, 725–733. doi: 10.1016/j.biopsych.2012.06.021
Kilonzo, V. W., Sweet, R. A., Glausier, J. R., and Pitts, M. W. (2020). Deficits in glutamic acid decarboxylase 67 immunoreactivity, parvalbumin interneurons, and perineuronal nets in the inferior colliculus of subjects with schizophrenia. Schizophr. Bull. 46, 1053–1059. doi: 10.1093/schbul/sbaa082
Lam, M., Chen, C. Y., Li, Z., Martin, A. R., Bryois, J., Ma, X., et al. (2019). Comparative genetic architectures of schizophrenia in East Asian and European populations. Nat. Genet. 51, 1670–1678.
Li, W., Wang, X., Zhao, J., Lin, J., Song, X. Q., Yang, Y., et al. (2012). Association study of myelin transcription factor 1-like polymorphisms with schizophrenia in Han Chinese population. Genes Brain Behav. 11, 87–93. doi: 10.1111/j.1601-183x.2011.00734.x
Li, W., Yang, Y., Lin, J., Wang, S., Zhao, J., Yang, G., et al. (2013). Association of serotonin transporter gene (SLC6A4) polymorphisms with schizophrenia susceptibility and symptoms in a Chinese-Han population. Prog. Neuropsychopharmacol. Biol. Psychiatry 44, 290–295. doi: 10.1016/j.pnpbp.2013.04.003
Li, Z., Chen, J., Yu, H., He, L., Xu, Y., Zhang, D., et al. (2017). Genome-wide association analysis identifies 30 new susceptibility loci for schizophrenia. Nat. Genet 49, 1576–1583.
Lichtenstein, P., Yip, B. H., Bjork, C., Pawitan, Y., Cannon, T. D., Sullivan, P. F., et al. (2009). Common genetic determinants of schizophrenia and bipolar disorder in Swedish families: a population-based study. Lancet 373, 234–239. doi: 10.1016/s0140-6736(09)60072-6
McCutcheon, R. A., Reis Marques, T., and Howes, O. D. (2020). Schizophrenia-an overview. JAMA Psychiatry 77, 201–210.
Mitchell, A. C., Jiang, Y., Peter, C., and Akbarian, S. (2015). Transcriptional regulation of GAD1 GABA synthesis gene in the prefrontal cortex of subjects with schizophrenia. Schizophr. Res. 167, 28–34. doi: 10.1016/j.schres.2014.10.020
Moghaddam, B., and Javitt, D. (2012). From revolution to evolution: the glutamate hypothesis of schizophrenia and its implication for treatment. Neuropsychopharmacology 37, 4–15. doi: 10.1038/npp.2011.181
Pardinas, A. F., Holmans, P., Pocklington, A. J., Escott-Price, V., Ripke, S., Carrera, N., et al. (2018). Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection. Nat. Genet 50, 381–389.
Pruim, R. J., Welch, R. P., Sanna, S., Teslovich, T. M., Chines, P. S., Gliedt, T. P., et al. (2010). LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26, 2336–2337. doi: 10.1093/bioinformatics/btq419
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet 81, 559–575. doi: 10.1086/519795
Represa, A., and Ben-Ari, Y. (2005). Trophic actions of GABA on neuronal development. Trends Neurosci. 28, 278–283. doi: 10.1016/j.tins.2005.03.010
Ripke, S., Walters, J. T., and O’donovan, M. C. Schizophrenia Working Group of the Psychiatric Genomics Consortium (2020). Mapping genomic loci prioritises genes and implicates synaptic biology in schizophrenia. MedRxiv 2020.09.12.20192922. doi: 10.1101/2020.09.12.20192922
Ruzicka, W. B., Subburaju, S., and Benes, F. M. (2015). Circuit- and diagnosis-specific DNA methylation changes at gamma-Aminobutyric acid-related genes in postmortem human hippocampus in schizophrenia and bipolar disorder. JAMA Psychiatry 72, 541–551. doi: 10.1001/jamapsychiatry.2015.49
Schizophrenia Working Group of the Psychiatric Genomics Consortium (2014). Biological insights from 108 schizophrenia-associated genetic loci. Nature 511, 421–427. doi: 10.1038/nature13595
Schmidt, M. J., and Mirnics, K. (2015). Neurodevelopment, GABA system dysfunction, and schizophrenia. Neuropsychopharmacology 40, 190–206. doi: 10.1038/npp.2014.95
Straub, R. E., Lipska, B. K., Egan, M. F., Goldberg, T. E., Callicott, J. H., Mayhew, M. B., et al. (2007). Allelic variation in GAD1 (GAD67) is associated with schizophrenia and influences cortical function and gene expression. Mol. Psychiatry 12, 854–869. doi: 10.1038/sj.mp.4001988
Sullivan, P. F. (2007). Spurious genetic associations. Biol. Psychiatry 61, 1121–1126. doi: 10.1016/j.biopsych.2006.11.010
Sullivan, P. F. (2017). How good were candidate gene guesses in schizophrenia genetics? Biol. Psychiatry 82, 696–697. doi: 10.1016/j.biopsych.2017.09.004
Sullivan, P. F., Kendler, K. S., and Neale, M. C. (2003). Schizophrenia as a complex trait: evidence from a meta-analysis of twin studies. Arch. Gen. Psychiatry 60, 1187–1192.
Szabo, G., Katarova, Z., and Greenspan, R. (1994). Distinct protein forms are produced from alternatively spliced bicistronic glutamic acid decarboxylase mRNAs during development. Mol. Cell Biol. 14, 7535–7545.
Tao, R., Davis, K. N., Li, C., Shin, J. H., Gao, Y., Jaffe, A. E., et al. (2018). GAD1 alternative transcripts and DNA methylation in human prefrontal cortex and hippocampus in brain development, schizophrenia. Mol. Psychiatry 23, 1496–1505.
White, T., Ho, B. C., Ward, J., O’leary, D., and Andreasen, N. C. (2006). Neuropsychological performance in first-episode adolescents with schizophrenia: a comparison with first-episode adults and adolescent control subjects. Biol. Psychiatry 60, 463–471.
Winkler, T. W., Day, F. R., Croteau-Chonka, D. C., Wood, A. R., Locke, A. E., Magi, R., et al. (2014). Quality control and conduct of genome-wide association meta-analyses. Nat. Protoc. 9, 1192–1212.
Yoon, J. H., Maddock, R. J., Rokem, A., Silver, M. A., Minzenberg, M. J., Ragland, J. D., et al. (2010). GABA concentration is reduced in visual cortex in schizophrenia and correlates with orientation-specific surround suppression. J. Neurosci. 30, 3777–3781.
Yu, H., Yan, H., Li, J., Li, Z., Zhang, X., Ma, Y., et al. (2017). Common variants on 2p16.1, 6p22.1 and 10q24.32 are associated with schizophrenia in Han Chinese population. Mol. Psychiatry 22, 954–960.
Yuan, H. Y., Chiou, J. J., Tseng, W. H., Liu, C. H., Liu, C. K., Lin, Y. J., et al. (2006). FASTSNP: an always up-to-date and extendable service for SNP function analysis and prioritization. Nucleic Acids Res. 34, W635–W641.
Keywords: GAD1, schizophrenia, single-nucleotide polymorphisms, Han Chinese, genetic association
Citation: Zhang L, Li Z, Liu Q, Shao M, Sun F, Su X, Song M, Zhang Y, Ding M, Lu Y, Liu J, Yang Y, Li M, Li W and Lv L (2021) Weak Association Between the Glutamate Decarboxylase 1 Gene (GAD1) and Schizophrenia in Han Chinese Population. Front. Neurosci. 15:677153. doi: 10.3389/fnins.2021.677153
Received: 07 March 2021; Accepted: 26 May 2021;
Published: 21 June 2021.
Edited by:
Yunpeng Wang, University of Oslo, NorwayReviewed by:
Weihua Yue, Peking University Sixth Hospital, ChinaMark Logue, Boston University, United States
Copyright © 2021 Zhang, Li, Liu, Shao, Sun, Su, Song, Zhang, Ding, Lu, Liu, Yang, Li, Li and Lv. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ming Li, limingkiz@mail.kiz.ac.cn; Wenqiang Li, lwq781603@163.com; Luxian Lv, lvx928@126.com
†These authors have contributed equally to this work