ORIGINAL RESEARCH article

Front. Genet., 05 March 2019

Sec. Behavioral and Psychiatric Genetics

Volume 10 - 2019 | https://doi.org/10.3389/fgene.2019.00175

SNP Variation of RELN Gene and Schizophrenia in a Chinese Population: A Hospital-Based Case–Control Study

  • 1. Department of Psychiatry, Shenzhen Kangning Hospital, Shenzhen Mental Health Center, Shenzhen, China

  • 2. Department of Psychiatry, Shantou University Medical College, Shantou, China

  • 3. Faculty of Forensic Medicine, Zhongshan School of Medicine, Sun Yat-sen University, Guangzhou, China

  • 4. Institute of Biomedical and Pharmaceutical Sciences, Guangdong University of Technology, Guangzhou, China

  • 5. Laboratory of Genomic and Precision Medicine, Wuxi School of Medicine, Jiangnan University, Wuxi, China

  • 6. Guangdong Province Translational Forensic Medicine Engineering Technology Research Center, Sun Yat-sen University, Guangzhou, China

Article metrics

View details

7

Citations

3,7k

Views

1,2k

Downloads

Abstract

Aims: We aimed to explore whether RELN contributes to the vulnerability and severity of clinical symptoms of schizophrenia (SZ) in a Chinese population.

Methods: The following were conducted in an adult Han Chinese population from southern China: case–control association analyses of 30 representative single nucleotide polymorphisms (SNPs) that were screened according to specific programs based on bioinformatics tools and former research and quantitative trait locus analyses with SNPs and psychiatric symptoms evaluated with the positive and negative symptoms scale.

Results: A 4-SNP haplotype consisting of rs362814, rs39339, rs540058, and rs661575 was found to be significantly associated with SZ even after Bonferroni correction (χ2 = 29.024, p = 6.42E-04, pBonf = 0.017), and the T-C-T-C haplotype was a protective factor for SZ (OR = 0.050, 95% CI = 0.004–0.705). Moreover, the 4-SNP haplotype showed a significant association with G16 (active social avoidance) after false discovery rate correction (χ2 = 28.620, p = 1.697E-04, pFDR = 0.025). In addition, P7 (hostility) was related to the haplotype comprising rs2229864, rs2535764, and rs262355 (χ2 = 31.424, p = 2.103E-05, padjustment = 0.019) in quantitative trait loci analyses.

Conclusion: Overall, this study showed several positive associations between RELN and SZ, as well as psychiatric symptoms, which not only supports the proposition that RELN is a susceptibility gene for SZ but also provides information on a genotype-phenotype correlation for SZ in a Chinese population.

Introduction

Schizophrenia (SZ, Online Mendelian Inheritance in Man [OMIM] 181500) is a common and serious psychiatric disorder with a lifetime prevalence estimate of 4.0 (1.6–12.1) per 1,000 individuals (Esan et al., 2012). SZ has been reported to be a predominantly genetic disorder, in which heritability is estimated to be 80% (Lichtenstein et al., 2009).

The human RELN gene (OMIM 600514) maps to chromosome 7q22 and encodes reelin, a large secreted glycoprotein, that is thought to be critical for cell positioning and neuronal migration by controlling cell-cell interactions during brain development (Hartfuss et al., 2003; Huang and D’Arcangelo, 2008; Niu et al., 2008). Recently, growing evidence has shown that the reelin protein might also be associated with neurotransmission, memory formation and synaptic plasticity (Herz and Chen, 2006), which have been demonstrated to also be damaged in SZ patients (Falkai et al., 2015). Decrease trends in reelin expression in patients with SZ have been found in brain and blood tissues (Impagnatiello et al., 1998; Guidotti et al., 2000; Nabil Fikri et al., 2017). Hence, low level or dysfunctional of reelin protein may cause deficits in neuronal development and cognitive function in adults and may play pathogenic roles in neuropsychiatric illnesses, such as SZ. This proposition was further supported by anatomical and in vivo studies that included the reeler mice and the heterozygous reeler mouse (HRM) model (Hill et al., 2006; Qiu et al., 2006; Tueting et al., 2006). The role of RELN as a potential risk for SZ has also been suggested by genetic association analysis, especially in GWASs (Shifman et al., 2008; Li et al., 2011, 2013; Ovadia and Shifman, 2011; Zhou et al., 2016). Over the past decade, dozens of SNPs in RELN gene loci have been reported to be associated with the onset and/or severity of clinical symptoms of SZ (Kahler et al., 2008; Shifman et al., 2008; Wedenoja et al., 2008, 2010; Need et al., 2009; Ben-David et al., 2010; Liu et al., 2010; Kuang et al., 2011; Li et al., 2013). However, the results remain controversial (Tost and Weinberger, 2011; Li et al., 2013; Bocharova et al., 2017).

In this study, we assumed that RELN may be related to the onset of SZ and the severity of some clinical symptoms. A case–control study had been performed in the Han Chinese population from southern China. To thoroughly understand the genetic basis of mental symptoms in SZ rather than only focus on verifying previous positive results and searching for new susceptible SNPs in RELN, we specifically performed quantitative trait locus (QTL) analyses in addition to qualitative association studies. Due to the high clinical and genetic heterogeneity of SZ, different clinical subtypes of SZ may be related to different genetic bases, we specifically collected patients with paranoid or undifferentiated SZ who experienced first onset or recurrence after drug withdrawal at least 1 month to improve the consistency of research subjects and to minimize the influence of drug on the scores of PANSS.

Materials and Methods

Subjects

The patient sample consisted of 102 unrelated SZ patients (46 females and 56 males; mean age, 32.42 ± 10.08 years) recruited from the Third People’s Hospital of Zhongshan City during 2012.11–2013.06. Patients were diagnosed by at least two psychiatrists according to the criteria of the Diagnostic and Statistical Manual of Mental Disorders (Fourth Edition, DSM-IV). Detailed information on clinical features, such as the onset time, symptoms and family history of mental illness, were obtained according to DIGS and FIGS. The PANSS was used to evaluate the severity of psychosis symptoms of SZ patients. It is one of the most common instruments applied to evaluate the severity of clinical symptoms in the world (Kay et al., 1987). PANSS comprises 33 items, including 30 psychopathological items which are usually divided into the positive subscale (7 items), the negative subscale (7 items) and the general psychopathology subscale (16 items) and 3 complementary attack risk items. The item level of the 33 items ranges from 1 to 7, with 1 equaling “no symptoms” (Kay et al., 1987). The Chinese version of PANSS, with acceptable validity and reliability, was used to evaluate the severity of psychosis symptoms of SZ patients in this study (Si et al., 2004). Paranoid or undifferentiated SZ patients who were at first onset, having never been treated, or at recurrence, having not taken any antipsychotics for at least 1 month, were enrolled in. Among them, patients with the both-negative type (i.e., the number of items with a score greater than or equal to 4 was less than three on both positive and negative subscales in PANSS) were excluded. All patients with ambiguous diagnoses or accompanied by neurological diseases, organic mental disorders or other symptomatic psychoses, or who had concomitant severe somatic disease, cancer, pregnancy or lactation, an immune or endocrine system disease, or a history of alcoholism or substance abuse were excluded.

A total of 169 healthy controls (75 females and 94 males; mean age 33.27 ± 8.63 years) matched with the patients in sex, age, and birthplace at the same time were randomly recruited from the volunteers who came to perform physical examination in the Affiliated Hospital of Sun Yat-sen University during 2013.06–2013.11. Volunteers were asked to provide detailed information on medical and family psychiatric histories. Subjects were excluded due to any of the following: positive family histories (first-degree relatives) of psychiatric illness; substance abuse; abnormal birth; febrile convulsions; juvenile adoption or residence with a single-parent family; pregnancy or lactation.

All participants were unrelated Han Chinese born and living in the southern China (mostly come from Guangdong province and nearby cities), and all their biological grandparents were Han Chinese ancestry. They provided written informed consent for participation, and the research protocol was approved by the Ethical Committee for Genetic Studies of Shantou University.

SNP Selection

The SNPs were based on previous studies (Shen et al., 2006; Doss et al., 2010; Wu and Jiang, 2013) and/or screened by SNP functional prediction softwares, including SIFT (Kumar et al., 2009), PolyPhen-2 (Adzhubei et al., 2010), and FASTSNP (Yuan et al., 2006). Information on the SNPs was obtained from dbSNP, HapMap and other human genome databases. A SNP would be approved only when its MAF was greater than or equal to 0.01 in the Chinese Beijing population in the HapMap database. The detailed strategies for SNP selection were as follows. (1) “Positive” SNPs that have been reported to be associated with SZ previously were selected once their MAFs were greater than or equal to 0.01 in the Chinese Beijing population in the HapMap database. (2) A two-step method was executed for the selection of exonic SNPs (Figure 1A). First, PolyPhen-2, the SIFT online service and the FASTSNP online service were employed for functional analysis of exonic SNPs. Next, by comparing the outputs of the three software programs, repeated SNPs with MAFs equal to or greater than 0.01 were chosen. (3) For intronic SNPs, a four-step method was conducted using the FASTSNP, SNP functional prediction (Xu et al., 2007), F-SNP online services and HapMap database successively, and SNPs with MAFs greater than or equal to 0.05 were preserved (Figure 1B). (4) By using FASTSNP, UTRscan (Shen et al., 2006), and UCSC Variant Annotation Integrator online software (Meyer et al., 2013), similar functional analyses of exonic SNPs were conducted among UTR and upstream SNPs of RELN, but a higher MAF threshold of 0.05 was implemented. Rs7341475, a “positive” SNP thought to be associated with SZ (Shifman et al., 2008; Ben-David et al., 2010) was not included in this study, as there was no significant evidence implying that rs7341475 may associate with SZ in the Han Chinese population (Li et al., 2011, 2013).

FIGURE 1

FIGURE 1

The general screening flows of exonic (A) and intronic (B) SNPs in the SNP selection phase.

SNP Genotyping

For each participant, 5 ml peripheral venous blood was collected after signing informed consent and before drug used. Genomic DNA was extracted using the phenol-chloroform method within 1–2 weeks and saved in −20°C refrigerator. SNP genotyping was performed with the Sequenom MassARRAY iPLEX Gold platform (Sequenom, San Diego, CA, United States) (Ikeda et al., 2008). Primers were designed using the MassARRAY Assay Design 3.1 software. Two amplification primer systems were built with the balance of capturing as many SNP loci as possible while limiting costs. Genomic fragments containing SNPs were amplified by polymerase chain reaction (PCR) in a total reaction volume of 5 μl, which included 20–50 ng of genomic DNA, in two 384-well plates using the ABI GeneAmp® 9700 384 Dual. Purified and specific genotyping primers were used to amplify target sites. Genotypes were automatically called by MassARRAY Typer 4.0 and verified manually. To ensure the accuracy of genotyping, 32 samples were duplicated for quality control, and no genotyping errors were found. Additionally, each 384-well plate had four blank controls. Any individual whose missing genotypes were greater than 50% was excluded from further statistical analysis.

Statistical Analysis

Statistical analyses including HWE, the description of MAF, single marker association of allele, the standardized measure of linkage disequilibrium (LD) coefficients (D’), haplotype block and haplotype frequency and association within the block were all assessed using Haploview v4.2 software (version 4.2, Broad Institute of MIT and Harvard, Cambridge, MA, United States) (Gabriel et al., 2002; Barrett et al., 2005). The haplotype frequency was estimated using the expectation maximization (EM) algorithm. The criterion for significances in Haploview v4.2 was set at p smaller than 0.05 for tests of expected HWE and corrected p-value after 10000 permutations.

The SPSS 19.0 were used to analyzing the demographic and genotypic frequency distribution of samples using t-tests or chi-square tests. Logistic regression analyses of population, sex, and age were used to generate population stratification assignments for all individuals. The criterion for significance in SPSS 19.0 was set at p smaller than 0.05. Genotype and multiple SNPs qualitative traits association tests between patients and controls were performed with Unphased 3.1.7 (Dudbridge, 2008). The overall test of association in Unphased 3.1.7 is a likelihood ratio test. And the odds ratio (OR) and 95% confidence interval (95% CI) were calculated automatically to evaluate the effects of alleles and haplotypes. To further investigate the possibility of complexity in the genotype-phenotype relationship, QTL analyses were also performed with Unphased v3.1.7 using 33 factor scores and 4 scale scores (i.e., the total scale score, the positive scale score, the negative scale score, and the general psychopathology scale score) from PANSS as quantitative traits. In QTL analyses, the addVal (i.e., additive genetic value) instead of OR was calculated automatically to show the estimated additive genetic value for special haplotype. For all multiple tests performed in Unphased 3.1.7, Bonferroni correction and FDR correction were employed to correct p-values and minimize the influence of type II errors. The corrected p-value threshold of 0.05 was used for significant after Bonferroni correction (pBonf = α × m, α is the desired overall alpha level of 0.05, and m is the number of hypotheses) and FDR correction (the corrected p-value was shown as pFDR) which was performed using R projection Version 3.1.21. A power analysis was performed using the GPower software for this study (Faul et al., 2009).

Results

SNP Selection and Genotyping

Initially, thirty-nine SNPs were screened, including twenty-nine “positive” SNPs reported previously and ten “functional” SNPs, from more than ten thousand RELN loci. Screening flows and the number of targeted SNPs in exonic and intronic zones in the RELN gene are shown in Figure 1. Six SNPs (i.e., rs362746, rs12705169, rs885995, rs11761011, rs16872603, and rs2237628) were ruled out for poor specificity or missing data, and three SNPs (i.e., rs3025962, rs73714410, and rs7811571) were monomorphic in the samples. Therefore, 30 SNPs in the RELN gene were genotyped successfully in 100 SZ patients (45 females and 55 males; mean age, 32.170 ± 10.000 years) and 163 healthy controls (74 females and 89 males; mean age, 33.350 ± 8.585 years). No differences were found between patients and controls in age (t = 1.015, df = 261, p = 0.311) and gender (χ2 = 0.004, df = 1, p = 1.000). The average call rate of all SNPs was 99.5% in total sample. No SNP deviation from HWE was found in the controls (Table 1). Among 30 SNPs, there are 2 exonic SNPs and 28 intronic SNPs. The rough distribution of the 30 SNPs is shown in Figure 2.

Table 1

dbSNP IDAllele (D/d)aCase
Control
Genotype (p)eAllele (P)e
nbHWE (p)cGenotype
MAFdnbHWE (p)cGenotype
MAFd
DDDdddDDDddd
rs11764507G/A1000.068722260.1701610.3811065230.1800.0620.768
rs17157643T/A1001.000782110.1151620.7341372500.0770.2130.145
rs6465938C/T1000.0464038220.4101630.3886083200.3770.0480.455
rs39339A/C1000.416772030.1301620.7141263510.1140.2830.589
rs155333G/A1000.1235434120.2901620.5467872120.2960.2010.878
rs262355T/A990.0715532120.2831620.7008070120.2900.1730.858
rs1510846T/C950.8342949170.4371620.2105072400.4690.4080.478
rs11496125C/T981.0003647150.3931620.4476571260.3800.8500.764
rs10435342C/T1000.277851320.0851630.9951431910.0640.5990.376
rs661575C/T990.1315134140.3081630.0728160220.3190.9580.794
rs607755A/G991.0003747150.3891600.1824788250.4310.4290.342
rs563264A/G1001.00095500.0251620.4431501110.0400.3860.356
rs540058T/C981.00093500.0261610.4461491110.0400.4040.370
rs727708A/G1000.9433051190.4451601.0003680440.4750.1840.076
rs12705141A/T1000.664742330.1451610.7471104560.1770.5280.338
rs2299356A/G1000.4063453130.3951610.8804782320.4530.3440.190
rs362691fG/C1001.000841510.0851631.0001213930.1380.1530.067
rs123714C/T1000.7635240100.3001630.614916480.2450.2470.169
rs123713C/T1000.7635040100.3001610.899906290.2480.3390.196
rs144525A/G990.536692910.1571610.3801015640.1990.4480.227
rs362626C/A1001.0004046140.3701610.3076182180.3660.6280.935
rs362814A/T990.6134143150.3691610.4976280190.3660.5210.959
rs362813T/C1000.9612649250.4951630.6944778380.4720.8570.614
rs362731T/C1001.0002650240.4901620.7594778370.4690.8560.642
rs362726T/C980.8792451230.4951631.0005679280.4140.1890.072
rs2535764C/T990.066682470.1921630.3521064890.2020.6470.769
rs362719C/A1000.9263647170.4051590.1317462230.3400.2580.132
rs2299334G/A1000.680663220.1801631.0001045270.2020.6120.527
rs11976900A/G990.486573480.2531610.1399254150.2610.9510.833
rs2229864fC/T1001.000682930.1751610.21410447100.2080.5450.354

Genotype frequencies, HWE tests, and single-SNP association analyses of 30 RELN SNPs in the SZ and control cases.

aMajor/minor allele, major, and minor alleles are denoted by D and d, respectively. bNumbers of samples that are successfully genotyped in patients and controls, respectively. cp-Values of the Hardy-Weinberg equilibrium (HWE) tests from Haploview v4.2. dThe minor allele frequency of 30 selected SNPs in this study. eUncorrected p-values of allelic and genotypic association analysis in this study. fTwo exonic SNPs selected in this report.

FIGURE 2

FIGURE 2

The genomic structure of RELN and the physical locations of the 30 tested SNPs. Among them, rs362691 and rs2229864 are located in the exonic areas (vertical bars), and the others are located in the intronic areas (spaces).

Single-SNP Association Analyses

Detailed information regarding the allelic and genotypic frequencies of the 30 SNPs from patients and controls was obtained (Table 1). No single polymorphism was found to be significantly associated with SZ neither in the total sample (Table 1) nor when grouped by sex (Table 2). The weak genotypic association of rs6465938 (χ2 = 6.087, df = 2, p = 0.048) did not withstand the FDR correction (pFDR = 0.713).

Table 2

SNPs IDsAllele
Genotype
Man
Woman
Man
Woman
χ2dfpχ2dfpχ2dfpχ2dfp
rs104353420.18210.670.34710.5560.56420.7540.29720.862
rs114961250.10910.7410.02310.880.15720.9250.35720.836
rs117645070.02910.8640.01710.8962.35120.3093.80820.149
rs119769000.08310.7730.00610.9360.07420.9640.23620.889
rs1237132.08210.1490.12810.7212.95620.2284.12920.127
rs1237142.51410.1130.12810.7212.86420.2394.12920.127
rs127051410.12810.722.40810.1210.36820.8322.15820.34
rs1445250.3810.5380.87310.352.94120.231.58520.453
rs15108461.1210.290.05610.8131.24720.5362.27820.32
rs1553330.03910.8440.00210.9673.36720.1860.48520.785
rs171576432.82910.0930.11110.7393.60320.1650.12710.721
rs22298640.00410.9511.33810.2470.56120.7552.54320.28
rs22993340.64510.4220.01310.911.07320.5850.07720.962
rs22993560.12310.7252.31210.1280.62620.7315.61720.06
rs25357640.810.3710.07910.7793.57320.1680.29520.863
rs2623550.110.7520.02510.8744.16120.1250.36720.832
rs3626260.46510.4950.14210.7060.66720.7164.03120.133
rs3626911.33910.2472.30410.1291.90320.3862.81220.245
rs3627190.72810.3930.98610.3211.57120.4561.02520.599
rs3627261.64110.21.0110.3153.04220.2191.53420.464
rs3627310.83410.3610.14110.7072.8720.2381.3520.509
rs3628130.95510.3290.14110.7072.54120.2811.3520.509
rs3628140.09410.759010.9970.19720.9064.08320.13
rs393390.91510.3390.1910.6631.34220.5112.8920.236
rs5400582.50410.1140.25310.6152.72220.2560.25810.612
rs5632642.68210.1020.25610.6132.87120.2380.2610.61
rs6077550.31410.5750.54210.4621.1220.5710.75620.685
rs64659380.36710.5450.19310.663.5620.1693.3820.184
rs6615752.41310.123.33910.0682.16620.3383.05820.217
rs7277080.63610.4253.32110.0681.45820.4824.56520.102

The sex differences in single-SNP association analyses of 30 RELN SNPs in SZ and control cases.

Haplotype Association Analyses

Five blocks were captured using Haploview v4.2 software in the total sample (Figure 3). The frequency and haplotype association within the block were shown in Supplementary Table S1. No haplotype analysis within the block survived before 10000 permutations. Haplotype associations of 2-SNP, 3-SNP, and 4-SNP were analyzed using Unphased v3.1.7 software. When performing logistic regression analyses, only birthplace shown the population stratification [p = 0.019, Exp (B) = 0.526, 95% CI 0.308–0.898]. To reduce the stratification, we employed birthplace, sex, and age as covariates in the association analyses. Two 3-SNP haplotype blocks (one consisting of rs362814, rs39339, and rs540058, and the other consisting of rs362691, rs362719, and rs362726) were identified before multiple adjustment (χ2 = 19.406, df = 6, p = 0.004, pBonf = 0.099 and χ2 = 20.621, df = 7, p = 0.004, pBonf = 0.122, respectively). The haplotype with four SNPs (rs362814, rs39339, rs540058, and rs661575) was significantly associated with SZ (χ2 = 29.024, df = 9, p = 6.42E-04), even after Bonferroni and FDR correction (both padjustment = 0.017). And the T-C-T-C haplotype of the four SNPs was more common in the controls (0.764% in case versus 6.173% in controls, OR = 0.050, 95% CI = 0.004–0.705). No sex effect was found in multiple SNP analyses (date not shown).

FIGURE 3

FIGURE 3

Linkage disequilibrium (LD) among the 30 RELN SNPs in the case and control groups. The pairwise LD R2 values of the sample set are illustrated in the matrix. The dark color indicates relatively strong LD. The R2 values of five blocks covering the RELN gene were larger than 0.75, indicating reasonable haplotype blocks.

QTL Association Analyses

For the QTL analyses, we used smoking, birthplace, age, and sex as covariates, as smoking may be associated with the severity of the negative symptoms (Zhang et al., 2012). Several factors rather than scale score were significantly associated with 3-SNP or 4-SNP haplotypes after FDR correction (Table 3, more information shown in Supplementary Table S2). The T-C-T haplotype and C-T-T haplotype of the 3-SNP consisting of rs2229864, rs2535764 and rs262355 were related to the severity of P7 (hostility) (AddVal = 0.748, 95% CI = 0.291–1.204 and AddVal = 0.718, 95% CI = 0.230–1.206). No significant result was found in SNP or 2-SNP QTL association analyses neither in total sample nor grouped by sex (data not shown). We did not performed 3-SNP and 4-SNP QTL association analyses in male or female patients due to the small sample.

Table 3

PANSS itemSNP1 IDSNP2 IDSNP3 IDSNP4 IDχ2dfppBonfpFDR
P7rs2229864rs2535764rs26235531.42462.10E-050.0190.019
G11rs12705141rs144525rs1510846rs229933436.900105.89E-050.0520.025
G16rs362814rs39339rs540058rs66157528.62071.70E-040.1510.025
P7rs2229864rs2535764rs262355rs36262637.772121.67E-040.1490.025
G14rs2229864rs2535764rs262355rs36262636.621122.57E-040.2290.026
S3rs2229864rs2535764rs262355rs36262633.360114.60E-040.410.041
P7rs17157643rs2229864rs2535764rs26235532.06089.09E-050.0810.025
S3rs17157643rs2229864rs2535764rs26235529.61182.48E-040.2210.026
G14rs10435342rs11496125rs11764507rs1197690033.28091.19E-040.1060.025
P4rs10435342rs11496125rs11764507rs1197690025.77475.53E-040.4930.045
P7rs155333rs17157643rs2229864rs253576430.95981.43E-040.1270.025
S3rs155333rs17157643rs2229864rs253576429.44582.65E-040.2360.026

Positive results of haplotype QTL association analyses.

P4: excitement; P7: hostility; G11: poor attention; G14: Poor impulse control; G16: active social avoidance; S3: emotional instability; pBonf: corrected p-values by Bonferroni correction; pFDR: corrected p-values by FDR correction.

Power

The GPower program was used to perform the power calculation. The size of this sample revealed a power of 98.873% to detect a significant association (α < 0.05) when given an effect size index of 0.5 (corresponding to a “medium” gene effect). However, the power narrowed down to 47.033% when the effect size index was 0.2 (corresponding to a “small” gene effect).

Discussion

Schizophrenia is a complex genetic disease with diverse clinical symptoms. Here, we performed qualitative and quantitative trait association analyses in a Han Chinese population from southern China to verify that RELN was a susceptibility gene for SZ. The results of multiple SNP analyses confirmed the genotype-phenotype relationship between RELN and SZ.

A 4-SNP haplotype consisting of rs362814, rs39339, rs540058, and rs661575 was observed to be significantly associated with SZ even after the Bonferroni correction, and the T-C-T-C haplotype may be a protective factor for SZ. Moreover, when performing QTL analysis, this 4-SNP haplotype had a significant association with G16 after FDR correction. These results suggested the importance of the four SNPs, which was supported by previous studies. A study of Han population origin from southwestern China demonstrated that the allele of rs362814 was more common in SZ cases and the A and T alleles took part in the building of a risk haplotype and a protective haplotype, respectively (Li et al., 2013). In our study, the T allele of rs362814 also took part in building a protective haplotype in the participants, but with different SNPs. Rs39339 was slightly significant in allelic association tests before the Bonferroni correction and nominally significant in a combined analysis in a Scandinavian population (Kahler et al., 2008). Rs540058 was reported to be associated with the severity of positive symptoms of SZ (Wedenoja et al., 2010), and rs661575 was associated with visual learning and memory in Finnish families (Wedenoja et al., 2008). The four SNPs are intronic and not in high linkage disequilibrium, so they are unlikely to have direct functional impacts on RELN and are unlikely to be co-inherited. However, SNPs in different functional areas may affect gene expression in different manners and have different effects on gene function (Manolio, 2010). Moreover, up to 52% of all SNPs that were associated with disease were intronic SNPs (Chen et al., 2010; Manolio, 2010). Hence, further studies on large samples are needed to verify the associations between the four SNPs and SZ, and rs362814 might be associated with SZ in an allele-dependent manner.

Rs2229864 is a synonymous mutation located in the 50th exon of the RELN gene. Rs2229864 was recognized to be detrimental in the screening process and was imbalanced in the allelic expression of RELN in SZ (Ovadia and Shifman, 2011). Rs2535764 and rs262355 are both intronic SNPs and were slightly significant in the allelic association test only before the Bonferroni correction (Kahler et al., 2008). Both rs262355 and rs2229864 failed to obtain definite positive results in subsequent studies (Li et al., 2011). However, a weak association with P7 and haplotype consisting of rs2229864, rs2535764, and rs262355 withstood both the FDR and Bonferroni correction in the QTL analysis, which demonstrates that the relationship may be genuine rather than a false positive in this study. In addition, 4-SNP haplotypes consisting of the three SNPs and rs362626 or rs17157643 were associated with P7, G14 and S3, or P7 and S3, which all reflect the impulsive and dangerous behavior of SZ patients in PANSS. Both single item scores and scale scores were used as quantitative traits instead of synthesized scale scores in this study to deeply understand the relationships between the severity of clinical symptoms and genetic bases as well as the category of symptoms and genetics. Overall, we presumed that rs2229864, rs2535764, and rs262355 may be closely related to the risk for SZ patients, and subsequent studies with large samples will be especially necessary.

Another interesting phenomenon was that the AA/GG ratio of rs727708 seemingly had opposite trends in the SZ patients and controls, although a significant difference was not found in the allele or genotype association analyses. Both rs727708 and rs540058 were thought to be related to the severity of positive symptoms (Wedenoja et al., 2010), so rs727708 may be related to both the susceptibility to and severity of SZ.

Both rs362691 and rs6465938 failed in single and multiple analyses. Rs362691 is a missense mutation located in the 22th exon of the RELN gene and can cause the Leu-Val amino acid change. Rs362691 had also been recognized as a detrimental mutation in the screening process. Researchers have found that rs362691 was not only associated with autism spectrum disorders (ASDs) (Wang et al., 2014) but might also take part in the influence of the RELN gene on the cognitive functions of healthy people (Baune et al., 2010). Haplotypes consisting of rs362691 rather than a SNP played roles in susceptibility to SZ in the Chinese Va population (Yang et al., 2013). A subtle but significant difference in the genotype frequency distribution of rs6465938 was found before multiple corrections. Notably, rs6465938 was not in HWE only in the patient samples. Moreover, rs6465938 showed a nominal association in a Scandinavian population (Kahler et al., 2008). Therefore, negative results in this study should not be sufficient evidence to deny the possible associations between rs362691 or rs6465938 and SZ.

When subjects were divided by sex, none of these SNPs was found to be significantly associated with SZ in this study. Due to the small sample size, QTL analyses of different sexes were not performed. Many studies have demonstrated that males and females show different clinical factors, such as age of onset, treatment, cognitive function and clinical symptomatology (Tang et al., 2007; Ochoa et al., 2012; Zhang et al., 2012; Thorup et al., 2014), as well as genetic information (Shifman et al., 2008; Liu et al., 2010). We agree with the notion that male SZ patients are different from female patients in terms of both clinical and genetic conditions. Our negative findings are mainly due to the small sample size and low statistical power, so studies with larger sample sizes would contribute to identifying sex differences in patients with SZ.

There are several possible explanations for the discrepancies between studies, such as genetic heterogeneity in distinct ethnic populations (Caucasian versus East Asian race, Han versus Va populations), the heterogeneity of SZ (subtypes and disease categories), population stratification, environmental exposure, cultures and diets (Li et al., 2011, 2013). Although we tried to reduce the effects of region and SZ subtype, the major limitation of the small sample size due to the rigorous standards still exists. The enrollment of patients principally from an outpatient service instead of an inpatient department may be another reason, as the overall PANSS scores may tend to be lower for outpatients than for hospitalized patients.

Conclusion

In summary, the current investigation showed qualitative and quantitative trait associations between genetic variants in RELN and SZ in a southern Han Chinese population. Several SNPs showed significant associations with the pathogenesis of SZ and/or the severity of psychiatric symptoms principally in multi-SNP analyses. Considering the limitation of our work, further investigations of genetic susceptibility among larger samples and inpatients are required to elucidate the role of the RELN polymorphisms in the risk and sex difference of SZ as well as the severity of clinic symptoms in the future.

Statements

Data availability statement

The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

Author contributions

HZ, Y-WS, and J-HC contributed to the conception and design of the study, and provided the approval for publication of the content. XL wrote the protocol and managed the literature searches and SNP screening. SC and LX organized the database and performed the statistical analysis. XL and SC wrote sections of the manuscript. All authors contributed to the manuscript revision and read and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Numbers 81871535 and 81273350).

Acknowledgments

We are grateful to all the voluntary donors of DNA samples in this study. We thank all the field workers who participated in the sample collection and diagnostic assessment, especially Mr. Ting-Yun Jiang, Ms. Wen-Wei Zhang, and Mr. Lei Liu. We acknowledge the contributions of Ms. Qiu-Lin Liu, Miss Lin Zhou, and Miss Lu Zong for the laboratory work to help with the statistical analysis.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00175/full#supplementary-material

Abbreviations

  • DIGS

    the Diagnostic Interview for Genetic Studies

  • DSM-IV

    the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition

  • FASTSNP

    function analysis and selection tool for single nucleotide polymorphisms

  • FDR

    false discovery rate

  • FIGS

    the Familial Interview for Genetic Studies

  • GWAS

    genome wide association study

  • HWE

    Hardy-Weinberg equilibrium

  • MAF

    minor allele frequency

  • OMIM

    Online Mendelian Inheritance in Man

  • OR

    odds ratios

  • PANSS

    the positive and negative symptoms scale

  • QTL

    quantitative trait loci

  • SIFT

    Sorting Tolerant From Intolerant

  • SNPs

    single nucleotide polymorphisms

  • SZ

    schizophrenia

  • UTR

    untranslated region.

References

  • 1

    AdzhubeiI. A.SchmidtS.PeshkinL.RamenskyV. E.GerasimovaA.BorkP.et al (2010). A method and server for predicting damaging missense mutations.Nat. Methods7248249. 10.1038/nmeth0410-248

  • 2

    BarrettJ. C.FryB.MallerJ.DalyM. J. (2005). Haploview: analysis and visualization of LD and haplotype maps.Bioinformatics21263265. 10.1093/bioinformatics/bth457

  • 3

    BauneB. T.KonradC.SuslowT.DomschkeK.BirosovaE.SehlmeyerC.et al (2010). The reelin (RELN) gene is associated with executive function in healthy individuals.Neurobiol. Learn. Mem.94446451. 10.1016/j.nlm.2010.08.002

  • 4

    Ben-DavidE.ShifmanS.International Schizophrenia Consortium (2010). Further investigation of the association between rs7341475 and rs17746501 and schizophrenia.Am. J. Med. Genet. B Neuropsychiatr. Genet.153B, 12441247. 10.1002/ajmg.b.31093

  • 5

    BocharovaA. V.StepanovV. A.MarusinA. V.KharkovV. N.VagaitsevaK. V.FedorenkoO. Y.et al (2017). Association study of genetic markers of schizophrenia and its cognitive endophenotypes.Russ. J. Genet.53139146. 10.1134/S1022795417010033

  • 6

    ChenR.DavydovE. V.SirotaM.ButteA. J. (2010). Non-synonymous and synonymous coding SNPs show similar likelihood and effect size of human disease association.PLoS One5:e13574. 10.1371/journal.pone.0013574

  • 7

    DossG. P.RajasekaranR.ArjunP.SethumadhavanR. (2010). Prioritization of candidate SNPs in colon cancer using bioinformatics tools: an alternative approach for a cancer biologist.Interdiscip. Sci.2320346. 10.1007/s12539-010-0003-3

  • 8

    DudbridgeF. (2008). Likelihood-based association analysis for nuclear families and unrelated subjects with missing genotype data.Hum. Hered.668798. 10.1159/000119108

  • 9

    EsanO. B.OjagbemiA.GurejeO. (2012). Epidemiology of schizophrenia–an update with a focus on developing countries.Int. Rev. Psychiatry24387392. 10.3109/09540261.2012.725219

  • 10

    FalkaiP.RossnerM. J.SchulzeT. G.HasanA.BrzozkaM. M.MalchowB.et al (2015). Kraepelin revisited: schizophrenia from degeneration to failed regeneration.Mol. Psychiatry20671676. 10.1038/mp.2015.35

  • 11

    FaulF.ErdfelderE.BuchnerA.LangA. G. (2009). Statistical power analyses using GPower 3.1: tests for correlation and regression analyses.Behav. Res. Methods4111491160. 10.3758/BRM.41.4.1149

  • 12

    GabrielS. B.SchaffnerS. F.NguyenH.MooreJ. M.RoyJ.BlumenstielB.et al (2002). The structure of haplotype blocks in the human genome.Science29622252229. 10.1126/science.1069424

  • 13

    GuidottiA.AutaJ.DavisJ. M.Di-Giorgi-GereviniV.DwivediY.GraysonD. R.et al (2000). Decrease in reelin and glutamic acid decarboxylase67 (GAD67) expression in schizophrenia and bipolar disorder: a postmortem brain study.Arch. Gen. Psychiatry5710611069. 10.1001/archpsyc.57.11.1061

  • 14

    HartfussE.ForsterE.BockH. H.HackM. A.LeprinceP.LuqueJ. M.et al (2003). Reelin signaling directly affects radial glia morphology and biochemical maturation.Development13045974609. 10.1242/dev.00654

  • 15

    HerzJ.ChenY. (2006). Reelin, lipoprotein receptors and synaptic plasticity.Nat. Rev. Neurosci.7850859. 10.1038/nrn2009

  • 16

    HillJ. J.HashimotoT.LewisD. A. (2006). Molecular mechanisms contributing to dendritic spine alterations in the prefrontal cortex of subjects with schizophrenia.Mol. Psychiatry11557566. 10.1038/sj.mp.4001792

  • 17

    HuangC. C.D’ArcangeloG. (2008). “The Reelin Gene and Its Functions in Brain Development,” in Reelin Glycoprotein: Structure, Biology and Roles in Health and Disease, ed.FatemiS. H. (New York, NY: Springer Science), 114.

  • 18

    IkedaS.SasazukiS.NatsukawaS.ShauraK.KoizumiY.KasugaY.et al (2008). Screening of 214 single nucleotide polymorphisms in 44 candidate cancer susceptibility genes: a case-control study on gastric and colorectal cancers in the Japanese population.Am. J. Gastroenterol.10314761487. 10.1111/j.1572-0241.2008.01810.x

  • 19

    ImpagnatielloF.GuidottiA. R.PesoldC.DwivediY.CarunchoH.PisuM. G.et al (1998). A decrease of reelin expression as a putative vulnerability factor in schizophrenia.Proc. Natl. Acad. Sci. U.S.A.951571815723. 10.1073/pnas.95.26.15718

  • 20

    KahlerA. K.DjurovicS.KulleB.JonssonE. G.AgartzI.HallH.et al (2008). Association analysis of schizophrenia on 18 genes involved in neuronal migration: MDGA1 as a new susceptibility gene.Am. J. Med. Genet. B Neuropsychiatr. Genet.147B, 10891100. 10.1002/ajmg.b.30726

  • 21

    KayS. R.FiszbeinA.OplerL. A. (1987). The positive and negative syndrome scale (PANSS) for schizophrenia.Schizophr. Bull.13261276. 10.1093/schbul/13.2.261

  • 22

    KuangW. J.SunR. F.ZhuY. S.LiS. B. (2011). A new single-nucleotide mutation (rs362719) of the reelin (RELN) gene associated with schizophrenia in female Chinese Han.Genet. Mol. Res.1016501658. 10.4238/vol10-3gmr1343

  • 23

    KumarP.HenikoffS.NgP. C. (2009). Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm.Nat. Protoc.410731081. 10.1038/nprot.2009.86

  • 24

    LiM.LuoX. J.XiaoX.ShiL.LiuX. Y.YinL. D.et al (2013). Analysis of common genetic variants identifies RELN as a risk gene for schizophrenia in Chinese population.World J. Biol. Psychiatry149199. 10.3109/15622975.2011.587891

  • 25

    LiW.SongX.ZhangH.YangY.JiangC.XiaoB.et al (2011). Association study of RELN polymorphisms with schizophrenia in Han Chinese population.Prog. Neuropsychopharmacol. Biol. Psychiatry3515051511. 10.1016/j.pnpbp.2011.04.007

  • 26

    LichtensteinP.YipB. H.BjorkC.PawitanY.CannonT. D.SullivanP. F.et al (2009). Common genetic determinants of schizophrenia and bipolar disorder in swedish families: a population-based study.Lancet373234239. 10.1016/S0140-6736(09)60072-6

  • 27

    LiuY.ChenP. L.McGrathJ.WolyniecP.FallinD.NestadtG.et al (2010). Replication of an association of a common variant in the Reelin gene (RELN) with schizophrenia in Ashkenazi Jewish women.Psychiatr. Genet.20184186. 10.1097/YPG.0b013e32833a220b

  • 28

    ManolioT. A. (2010). Genomewide association studies and assessment of the risk of disease.N. Engl. J. Med.363166176. 10.1056/NEJMra0905980

  • 29

    MeyerL. R.ZweigA. S.HinrichsA. S.KarolchikD.KuhnR. M.WongM.et al (2013). The UCSC genome browser database: extensions and updates 2013.Nucleic Acids Res.41D64D69. 10.1093/nar/gks1048

  • 30

    Nabil FikriR. M.NorlelawatiA. T.Nour El-HudaA. R.HanisahM. N.KartiniA.NorsidahK.et al (2017). Reelin (RELN) DNA methylation in the peripheral blood of schizophrenia.J. Psychiatr. Res.882837. 10.1016/j.jpsychires.2016.12.020

  • 31

    NeedA. C.GeD.WealeM. E.MaiaJ.FengS.HeinzenE. L.et al (2009). A genome-wide investigation of SNPs and CNVs in schizophrenia.PLoS Genet.5:e1000373. 10.1371/journal.pgen.1000373

  • 32

    NiuS.YabutO.D’ArcangeloG. (2008). The Reelin signaling pathway promotes dendritic spine development in hippocampal neurons.J. Neurosci.281033910348. 10.1523/JNEUROSCI.1917-08.2008

  • 33

    OchoaS.UsallJ.CoboJ.LabadX.KulkarniJ. (2012). Gender differences in schizophrenia and first-episode psychosis: a comprehensive literature review.Schizophr. Res. Treatment.2012:916198. 10.1155/2012/916198

  • 34

    OvadiaG.ShifmanS. (2011). The genetic variation of RELN expression in schizophrenia and bipolar disorder.PLoS One6:e19955. 10.1371/journal.pone.0019955

  • 35

    QiuS.KorwekK. M.Pratt-DavisA. R.PetersM.BergmanM. Y.WeeberE. J. (2006). Cognitive disruption and altered hippocampus synaptic function in reelin haploinsufficient mice.Neurobiol. Learn. Mem.85228242. 10.1016/j.nlm.2005.11.001

  • 36

    ShenJ.DeiningerP. L.ZhaoH. (2006). Applications of computational algorithm tools to identify functional SNPs in cytokine genes.Cytokine356266. 10.1016/j.cyto.2006.07.008

  • 37

    ShifmanS.JohannessonM.BronsteinM.ChenS. X.CollierD. A.CraddockN. J.et al (2008). Genome-wide association identifies a common variant in the reelin gene that increases the risk of schizophrenia only in women.PLoS Genet.4:e28. 10.1371/journal.pgen.0040028

  • 38

    SiT.YangJ.ShuL.WangX.KongQ.ZhouM.et al (2004). The reliability,validity of PANSS and its Implication.Chin. Ment. Health J.14547.

  • 39

    TangY. L.GillespieC. F.EpsteinM. P.MaoP. X.JiangF.ChenQ.et al (2007). Gender differences in 542 Chinese inpatients with schizophrenia.Schizophr. Res.978896. 10.1016/j.schres.2007.05.025

  • 40

    ThorupA.AlbertN.BertelsenM.PetersenL.JeppesenP.Le QuackP.et al (2014). Gender differences in first-episode psychosis at 5-year follow-up–two different courses of disease? Results from the OPUS study at 5-year follow-up.Eur. Psychiatry294451. 10.1016/j.eurpsy.2012.11.005

  • 41

    TostH.WeinbergerD. R. (2011). RELN rs7341475 and schizophrenia risk: confusing, yet somehow intriguing.Biol. Psychiatry69:e19. 10.1016/j.biopsych.2010.10.022

  • 42

    TuetingP.DoueiriM. S.GuidottiA.DavisJ. M.CostaE. (2006). Reelin down-regulation in mice and psychosis endophenotypes.Neurosci. Biobehav. Rev.3010651077. 10.1016/j.neubiorev.2006.04.001

  • 43

    WangZ.HongY.ZouL.ZhongR.ZhuB.ShenN.et al (2014). Reelin gene variants and risk of autism spectrum disorders: an integrated meta-analysis.Am. J. Med. Genet. B Neuropsychiatr. Genet.165B, 192200. 10.1002/ajmg.b.32222

  • 44

    WedenojaJ.LoukolaA.Tuulio-HenrikssonA.PaunioT.EkelundJ.SilanderK.et al (2008). Replication of linkage on chromosome 7q22 and association of the regional Reelin gene with working memory in schizophrenia families.Mol. Psychiatry13673684. 10.1038/sj.mp.4002047

  • 45

    WedenojaJ.Tuulio-HenrikssonA.SuvisaariJ.LoukolaA.PaunioT.PartonenT.et al (2010). Replication of association between working memory and Reelin, a potential modifier gene in schizophrenia.Biol. Psychiatry67983991. 10.1016/j.biopsych.2009.09.026

  • 46

    WuJ.JiangR. (2013). Prediction of deleterious nonsynonymous single-nucleotide polymorphism for human diseases.Sci. World J.2013:675851. 10.1155/2013/675851

  • 47

    XuZ.KaplanN. L.TaylorJ. A. (2007). Tag SNP selection for candidate gene association studies using HapMap and gene resequencing data.Eur. J. Hum. Genet.1510631070. 10.1038/sj.ejhg.5201875

  • 48

    YangX. B.KangC.LiuH.YangJ. (2013). Association study of the reelin (RELN) gene with Chinese Va schizophrenia.Psychiatr. Gene.23:138. 10.1097/YPG.0b013e32835d705c

  • 49

    YuanH. Y.ChiouJ. J.TsengW. H.LiuC. H.LiuC. K.LinY. J.et al (2006). FASTSNP: an always up-to-date and extendable service for SNP function analysis and prioritization.Nucleic Acids Res.34W635W641. 10.1093/nar/gkl236

  • 50

    ZhangX. Y.ChenD. C.XiuM. H.HaileC. N.SunH.LuL.et al (2012). Cigarette smoking and cognitive function in Chinese male schizophrenia: a case-control study.PLoS One7:e36563. 10.1371/journal.pone.0036563

  • 51

    ZhouZ.HuZ.ZhangL.LiuH.LiuZ.DuJ.et al (2016). Identification of RELN variation p.Thr3192Ser in a Chinese family with schizophrenia.Sci. Rep.6:24327. 10.1038/srep24327

Summary

Keywords

RELN, schizophrenia, SNP, association analysis, quantitative trait loci analyses, psychiatric symptoms, bioinformatics tools

Citation

Luo X, Chen S, Xue L, Chen J-H, Shi Y-W and Zhao H (2019) SNP Variation of RELN Gene and Schizophrenia in a Chinese Population: A Hospital-Based Case–Control Study. Front. Genet. 10:175. doi: 10.3389/fgene.2019.00175

Received

21 August 2018

Accepted

18 February 2019

Published

05 March 2019

Volume

10 - 2019

Edited by

Weihua Yue, Peking University Sixth Hospital, China

Reviewed by

Xiaoping Wang, Second Xiangya Hospital of Central South University, China; Qiang Wang, West China Hospital of Sichuan University, China

Updates

Copyright

*Correspondence: Jian-Huan Chen, Yan-Wei Shi, Hu Zhao,

Co-first authors

This article was submitted to Behavioral and Psychiatric Genetics, a section of the journal Frontiers in Genetics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics