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

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, p Bonf = 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, p FDR = 0.025). In addition, P7 (hostility) was related to the haplotype comprising rs2229864, rs2535764, and rs262355 (χ2 = 31.424, p = 2.103E-05, p adjustment = 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

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.2[1]. A power analysis was performed using the G∗Power 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

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

dbSNP ID	Allele (D/d)a	Case						Control						Genotype (p)e	Allele (P)e
		nb	HWE (p)c	Genotype			MAFd	nb	HWE (p)c	Genotype			MAFd
				DD	Dd	dd				DD	Dd	dd
rs11764507	G/A	100	0.068	72	22	6	0.170	161	0.381	106	52	3	0.180	0.062	0.768
rs17157643	T/A	100	1.000	78	21	1	0.115	162	0.734	137	25	0	0.077	0.213	0.145
rs6465938	C/T	100	0.046	40	38	22	0.410	163	0.388	60	83	20	0.377	0.048	0.455
rs39339	A/C	100	0.416	77	20	3	0.130	162	0.714	126	35	1	0.114	0.283	0.589
rs155333	G/A	100	0.123	54	34	12	0.290	162	0.546	78	72	12	0.296	0.201	0.878
rs262355	T/A	99	0.071	55	32	12	0.283	162	0.700	80	70	12	0.290	0.173	0.858
rs1510846	T/C	95	0.834	29	49	17	0.437	162	0.210	50	72	40	0.469	0.408	0.478
rs11496125	C/T	98	1.000	36	47	15	0.393	162	0.447	65	71	26	0.380	0.850	0.764
rs10435342	C/T	100	0.277	85	13	2	0.085	163	0.995	143	19	1	0.064	0.599	0.376
rs661575	C/T	99	0.131	51	34	14	0.308	163	0.072	81	60	22	0.319	0.958	0.794
rs607755	A/G	99	1.000	37	47	15	0.389	160	0.182	47	88	25	0.431	0.429	0.342
rs563264	A/G	100	1.000	95	5	0	0.025	162	0.443	150	11	1	0.040	0.386	0.356
rs540058	T/C	98	1.000	93	5	0	0.026	161	0.446	149	11	1	0.040	0.404	0.370
rs727708	A/G	100	0.943	30	51	19	0.445	160	1.000	36	80	44	0.475	0.184	0.076
rs12705141	A/T	100	0.664	74	23	3	0.145	161	0.747	110	45	6	0.177	0.528	0.338
rs2299356	A/G	100	0.406	34	53	13	0.395	161	0.880	47	82	32	0.453	0.344	0.190
rs362691f	G/C	100	1.000	84	15	1	0.085	163	1.000	121	39	3	0.138	0.153	0.067
rs123714	C/T	100	0.763	52	40	10	0.300	163	0.614	91	64	8	0.245	0.247	0.169
rs123713	C/T	100	0.763	50	40	10	0.300	161	0.899	90	62	9	0.248	0.339	0.196
rs144525	A/G	99	0.536	69	29	1	0.157	161	0.380	101	56	4	0.199	0.448	0.227
rs362626	C/A	100	1.000	40	46	14	0.370	161	0.307	61	82	18	0.366	0.628	0.935
rs362814	A/T	99	0.613	41	43	15	0.369	161	0.497	62	80	19	0.366	0.521	0.959
rs362813	T/C	100	0.961	26	49	25	0.495	163	0.694	47	78	38	0.472	0.857	0.614
rs362731	T/C	100	1.000	26	50	24	0.490	162	0.759	47	78	37	0.469	0.856	0.642
rs362726	T/C	98	0.879	24	51	23	0.495	163	1.000	56	79	28	0.414	0.189	0.072
rs2535764	C/T	99	0.066	68	24	7	0.192	163	0.352	106	48	9	0.202	0.647	0.769
rs362719	C/A	100	0.926	36	47	17	0.405	159	0.131	74	62	23	0.340	0.258	0.132
rs2299334	G/A	100	0.680	66	32	2	0.180	163	1.000	104	52	7	0.202	0.612	0.527
rs11976900	A/G	99	0.486	57	34	8	0.253	161	0.139	92	54	15	0.261	0.951	0.833
rs2229864f	C/T	100	1.000	68	29	3	0.175	161	0.214	104	47	10	0.208	0.545	0.354

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

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

SNPs IDs	Allele						Genotype
	Man			Woman			Man			Woman
	χ2	df	p	χ2	df	p	χ2	df	p	χ2	df	p
rs10435342	0.182	1	0.67	0.347	1	0.556	0.564	2	0.754	0.297	2	0.862
rs11496125	0.109	1	0.741	0.023	1	0.88	0.157	2	0.925	0.357	2	0.836
rs11764507	0.029	1	0.864	0.017	1	0.896	2.351	2	0.309	3.808	2	0.149
rs11976900	0.083	1	0.773	0.006	1	0.936	0.074	2	0.964	0.236	2	0.889
rs123713	2.082	1	0.149	0.128	1	0.721	2.956	2	0.228	4.129	2	0.127
rs123714	2.514	1	0.113	0.128	1	0.721	2.864	2	0.239	4.129	2	0.127
rs12705141	0.128	1	0.72	2.408	1	0.121	0.368	2	0.832	2.158	2	0.34
rs144525	0.38	1	0.538	0.873	1	0.35	2.941	2	0.23	1.585	2	0.453
rs1510846	1.12	1	0.29	0.056	1	0.813	1.247	2	0.536	2.278	2	0.32
rs155333	0.039	1	0.844	0.002	1	0.967	3.367	2	0.186	0.485	2	0.785
rs17157643	2.829	1	0.093	0.111	1	0.739	3.603	2	0.165	0.127	1	0.721
rs2229864	0.004	1	0.951	1.338	1	0.247	0.561	2	0.755	2.543	2	0.28
rs2299334	0.645	1	0.422	0.013	1	0.91	1.073	2	0.585	0.077	2	0.962
rs2299356	0.123	1	0.725	2.312	1	0.128	0.626	2	0.731	5.617	2	0.06
rs2535764	0.8	1	0.371	0.079	1	0.779	3.573	2	0.168	0.295	2	0.863
rs262355	0.1	1	0.752	0.025	1	0.874	4.161	2	0.125	0.367	2	0.832
rs362626	0.465	1	0.495	0.142	1	0.706	0.667	2	0.716	4.031	2	0.133
rs362691	1.339	1	0.247	2.304	1	0.129	1.903	2	0.386	2.812	2	0.245
rs362719	0.728	1	0.393	0.986	1	0.321	1.571	2	0.456	1.025	2	0.599
rs362726	1.641	1	0.2	1.01	1	0.315	3.042	2	0.219	1.534	2	0.464
rs362731	0.834	1	0.361	0.141	1	0.707	2.87	2	0.238	1.35	2	0.509
rs362813	0.955	1	0.329	0.141	1	0.707	2.541	2	0.281	1.35	2	0.509
rs362814	0.094	1	0.759	0	1	0.997	0.197	2	0.906	4.083	2	0.13
rs39339	0.915	1	0.339	0.19	1	0.663	1.342	2	0.511	2.89	2	0.236
rs540058	2.504	1	0.114	0.253	1	0.615	2.722	2	0.256	0.258	1	0.612
rs563264	2.682	1	0.102	0.256	1	0.613	2.871	2	0.238	0.26	1	0.61
rs607755	0.314	1	0.575	0.542	1	0.462	1.12	2	0.571	0.756	2	0.685
rs6465938	0.367	1	0.545	0.193	1	0.66	3.56	2	0.169	3.38	2	0.184
rs661575	2.413	1	0.12	3.339	1	0.068	2.166	2	0.338	3.058	2	0.217
rs727708	0.636	1	0.425	3.321	1	0.068	1.458	2	0.482	4.565	2	0.102

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

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

Positive results of haplotype QTL association analyses.

PANSS item	SNP1 ID	SNP2 ID	SNP3 ID	SNP4 ID	χ2	df	p	pBonf	pFDR
P7	rs2229864	rs2535764	rs262355		31.424	6	2.10E-05	0.019	0.019
G11	rs12705141	rs144525	rs1510846	rs2299334	36.900	10	5.89E-05	0.052	0.025
G16	rs362814	rs39339	rs540058	rs661575	28.620	7	1.70E-04	0.151	0.025
P7	rs2229864	rs2535764	rs262355	rs362626	37.772	12	1.67E-04	0.149	0.025
G14	rs2229864	rs2535764	rs262355	rs362626	36.621	12	2.57E-04	0.229	0.026
S3	rs2229864	rs2535764	rs262355	rs362626	33.360	11	4.60E-04	0.41	0.041
P7	rs17157643	rs2229864	rs2535764	rs262355	32.060	8	9.09E-05	0.081	0.025
S3	rs17157643	rs2229864	rs2535764	rs262355	29.611	8	2.48E-04	0.221	0.026
G14	rs10435342	rs11496125	rs11764507	rs11976900	33.280	9	1.19E-04	0.106	0.025
P4	rs10435342	rs11496125	rs11764507	rs11976900	25.774	7	5.53E-04	0.493	0.045
P7	rs155333	rs17157643	rs2229864	rs2535764	30.959	8	1.43E-04	0.127	0.025
S3	rs155333	rs17157643	rs2229864	rs2535764	29.445	8	2.65E-04	0.236	0.026

Power

The G∗Power 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.

Data Availability

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.

Conflict of Interest Statement

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.