PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 11750.
Published online 2017 September 18. doi:  10.1038/s41598-017-10507-7
PMCID: PMC5603592

The common variants implicated in microstructural abnormality of first episode and drug-naïve patients with schizophrenia

Abstract

Both post-mortem and neuroimaging studies have identified abnormal white matter (WM) microstructure in patients with schizophrenia. However, its genetic underpinnings and relevant biological pathways remain unclear. In order to unravel the genes and the pathways associated with abnormal WM microstructure in schizophrenia, we recruited 100 first-episode, drug-naïve patients with schizophrenia and 140 matched healthy controls to conduct genome-wide association analysis of fractional anisotropy (FA) value measured using diffusing tensor imaging (DTI), followed by multivariate association study and pathway enrichment analysis. The results showed that one intergenic SNP (rs11901793), which is 20 kb upstream of CXCR7 gene on chromosome 2, was associated with the total mean FA values with genome-wide significance (p = 4.37 × 10−8), and multivariate association analysis identified a strong association between one region-specific SNP (rs10509852), 400 kb upstream of SORCS1 gene on chromosome 10, and the global trait of abnormal WM microstructure (p = 1.89 × 10−7). Furthermore, one pathway that is involved in cell cycle regulation, REACTOME_CHROMOSOME _MAINTENANCE, was significantly enriched by the genes that were identified in our study (p = 1.54 × 10−17). In summary, our study provides suggestive evidence that abnormal WM microstructure in schizophrenia is associated with genes that are likely involved in diverse biological signals and cell-cycle regulation although further replication in a larger independent sample is needed.

Introduction

Schizophrenia is a chronic mental disorder that affects about 1% of the population all over the world, causing a heavy burden to the family and society1. Although family and twin studies showed that schizophrenia has a high heritability rate of 80%, its prevalence cannot be explained by the model of monogenic disorder2,3. Schizophrenia is a complex illness that is caused by multiple genetic variants of small to modest effect size4,5. Over the last decade, there have been a number of case-control genome-wide association studies (GWAS) of schizophrenia with different sample size in different populations. Whereas some genes have been repeatedly identified, they explained only a small fraction of schizophrenia prevalence6,7. With the missing heritability yet to be explored, some studies tried to expand their sample size to increase the power of GWAS8, which, although is promising, can be limited by funding and sample source for replication. Meanwhile, it has been recognized that schizophrenia is a complex disease without reliable biomarkers and its diagnosis has been largely based on the clinical descriptions of psychiatric practitioners, which can be biased by many confounding factors, such as the personal experience of clinicians, lack of longitudinal observation etc.

To address such challenges, one of the possible strategies is to use quantitative traits (QTs) of schizophrenia, which lay between the genes and the clinical diagnosis, to increase the power of study9,10. In fact, it has been shown that complex disorder per se displays the feature of QT involving multiple genes, implying that disorder’s genetic liability is distributed quantitatively rather than qualitatively11. Indeed, studies using measurable and objective QTs as proxy phenotypes to map genetic variants of schizophrenia have come up with promising and telling findings. For instance, Potkin et al. used task-related blood oxygen level–dependent (BOLD) signal from functional magnetic resonance imaging (fMRI) as the QTs to map susceptible variants in GWAS of schizophrenia. They identified that a number of variants in 3 genes and chromosomal regions (ROBO1-ROBO2, TNIK, and CTXN3-SLC12A2) were significantly associated with BOLD signal differences in schizophrenia12. Dickinson et al. used g score, which is a composite score that combines 6 neuro-cognitive dimensions, as the QTs in a genome association study to detect potential genetic variants associated with schizophrenia. This approach has led to unveil the effects of SCN2A on the deficits in general cognitive ability in schizophrenia and the genetic directional trend was further validated by MRI scanning and mRNA expression in schizophrenia13. In Chinese population, Wang et al. chose grey matter (GM) volume as QTs to map common variants related to schizophrenia and the results showed that the variants in the HS3ST5, TBXAS1 and PIK3C2G genes were significantly associated with GM volume changes in hOC3vL, vermisL10 and vermisR114. Wang et al. later used functional connectivity as the QTs in GWAS of schizophrenia and found that the CHRM3 gene was significantly associated with functional connectivity deficits in schizophrenia15. These studies highlighted the efficacy of QTs as the endophenotype in exploring the pathogenesis of schizophrenia.

Diffusion tensor imaging (DTI) is a noninvasive neuroimaging technique that yields information about the microstructure of brain white matter (WM) and the phenotype has been found to be heritable16,17. Moreover, abnormal WM microstructure has been repeatedly identified in both first-episode and chronic schizophrenic patients18,19. According to the study by Nesvag et al., WM abnormality was also related to the response to antipsychotic treatment20. Although there were several genetic studies of microstructure WM abnormality in schizophrenia21,22, no GWAS of schizophrenia using WM measures as QTs has been conducted. Thus, our study was carried out with the following aims: firstly, to investigate WM microstructure deficits in first-episode, drug-naïve patients with schizophrenia; secondly, to map common variants/genes underlying the microstructure deficits in patients with schizophrenia; finally, to annotate biological pathway enriched by mapped common variants/genes which are implicated in WM microstructural abnormality in patients with schizophrenia.

Materials and Methods

Participants

One hundred first-episode patients with schizophrenia were recruited from the West China Hospital of Sichuan University, and 140 demographically matched healthy controls without family history of mental disorders (recruited from the local neighborhood) were included in the study. Patients were interviewed with Structured Clinical Interview for Diagnostic and Statistical Manual of Mental Disorders (SCID, Patient version) by trained psychiatrists, and the diagnosis of schizophrenia was made based on criteria in Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV). Healthy controls were screened with SCID (Non-patient version) to ensure the absence of psychiatric illness. Patients were followed up for at least six months in order to confirm the diagnosis of schizophrenia. All patients did not use any antipsychotic medication at the time of clinical and MRI scanning. Both patients and controls were excluded if they had one of following conditions: (1) organic cerebral diseases; (2) neurological diseases; (3) severe physical diseases; and (4) axis I and II diagnosis other than schizophrenia according to DSM-IV criteria. Informed written consent was obtained from all participants after explaining the purpose and procedure of the study. This study was approved by the Institutional Review Broad (IRB) of West China Hospital, Sichuan University and all the procedures of sample collection and data analysis described here are in strict accordance with relevant guidelines and regulations.

MRI data acquisition and preprocessing

Whole brain diffusion-weighted images were recorded along 15 gradient directions (b = 1000 s/mm2, number of excitations = 2) together with one unweighted (b = 0) image (42 images in total). Each image was acquired using a single-shot spin echo planar imaging (EPI) sequence [repetition time (TR) = 10000 ms, echo time (TE) = 70.8 ms, slice thickness = 3.0 mm with no gap, field of view = 240 mm2, matrix size = 256 × 256, voxel resolution = 0.94 × 0.94 × 3 mm3]. Two experienced neuroradiologists qualitatively inspected the raw MRI data. No gross abnormalities were observed for any subject. DTI data were preprocessed using the FMRIB’s Diffusion Toolbox (FDT) within FSL (http://www.fmrib.ox.ac.uk/fsl). The preprocessing included the following steps: (1) The effects of head motion and image distortion caused by eddy currents were corrected by applying an affine alignment to register all other diffusion images to the b0 images in the original DTI data. (2) The non-brain tissue and background noise were removed from b0 image using BET23. The diffusion tensor for each voxel was estimated by the multivariate linear fitting algorithm, and the tensor matrix was diagonalized to obtain its three pairs of eigenvalues (L1, L2, and L3) and eigenvectors. (3) The voxel-wise value of fractional anisotropy (FA) was calculated for each subject. Finally, (4) FA maps of all subjects were normalized to a 2 × 2 × 2 mm3 Montreal Neurological Institute standard space (MNI152) and smoothed with a 6 mm full-width at the half- maximum Gaussian kernel.

Genotyping, quality control and imputation

The genotype data has been used by our previous studies and the details of genotyping and quality control have been published elsewhere14,15. In brief, DNA extracted from whole blood sample was genotyped using HumanOmniZhongHua-8 BeadChip (Illumina, San Diego, CA) and the raw genotype data underwent systematic quality control steps (QCs) including missing rate per genotype (< 5%) and per individuals (< 3%), minor allele frequency (MAF, < 0.01), Hardy-Weinberg equilibrium tests (p < 0.001), cryptic relatedness and heterozygosity rate. In total, 109,923 SNPs and four individuals (two patients and two healthy controls) were excluded for not meeting the criteria of quality controls. The steps described above were implemented in PLINK1.0724. Furthermore, the standard pipeline was used to detect outliers of population stratification by using EIGENSTART25 with the first two components of PCA included as covariates in the subsequent analyses.

Using IMPUTE226, and with the 1000 Genomes Project phase I dataset as a reference, we implemented imputation of the missing genotypes based on pre-phased haplotypes from SHAPEIT227. In addition to same quality control steps (missing rate per genotype and per individuals, MAF and Hardy-Weinberg equilibrium tests), we excluded the imputed genotypes with INFO less than 0.6. After quality control, 7 patients and 6 controls were excluded due to low genotyping quality and cryptic relatedness.

Statistical analysis

Comparison of FA values between patients and healthy controls

The voxel-wise comparisons of FA maps were performed using ANCOVA in SPM (Statistical Parametric Mapping) with age and sex being included as covariates. We reported the results at q < 0.01 after multiple corrections using the false discovery rate (FDR) method.

Univariate QT analysis and functional annotation

In order to detect genetic loci associated with FA values in brain regions identified from above, a mixed linear model in MixABEL of ABEL packages28 was used to test association of FA values in the brain regions showing significant inter-group difference with the interaction effect of genotype × group under the assumption of additive dosage effect; controlling for age, sex and the first two components from population PCA. Following the QT analysis, the top associated SNPs were annotated in terms of histone modification29, regulatory motif alterations30 and Regulome DB scores using both Haploreg31 and Regulome DB32.

In order to detect potential effect of gene-gene interaction (epistasis) on the QTs in our study, an analysis using model-based multifactor dimension reduction (MB-MDR) was conducted with two order combination of the top SNPs identified from univariate QT analysis as epistasis variables and covariates such as age, gender and main effect of single SNP being controlled for. In brief, MB-MDR consists of 3 steps: (1) an association test for each multilocus cell with the phenotype: a logistic regression for binary traits, a linear regression for was performed. Each genotype cell is then assigned to one of three categories, high risk (H), low risk (L) or no evidence (0), as a result of the association test with a liberal threshold of 0.1 for assigning the genotype to a risk category, H or L. (2) Two new association analysis, one for each risk category, high and low, on the outcome variable was implemented, with the corresponding unadjusted p value for each test (PH and PL) being generated. The minimum between PH and PL is chosen as minimum P value. (3) The significance of the specified models was decided through 1000 permutations on the maximum Wald statistic from step 233.

Multivariate association study

To identify the shared genetic effects within multiple regions of abnormal white matter microstructure in schizophrenia, we conducted a multivariate genotype-phenotype analysis using efficient multivariate genotype-phenotype analysis for genome-wide association studies (Trait-based Association Test that uses Extended Simes procedure, TATES)34. This analysis presents a global trait-based p value (PT) by combining the p values derived from standard univariate GWAS and correcting for the observed correctional structure between the phenotypes (i.e. multiple abnormal white matter microstructure in the present study). The global trait-based results from TATES were then used to conduct gene-based testing using Gene-Based Association Test Using Extended Simes Procedure (GATES)35 in Knowledge-based mining system for Genome-wide Genetic studies (KGG) platform. After mapping SNPs onto genes according to the gene coordinate information from NCBI (SNPs within 5 kilobase pairs of each gene were also assigned into the gene) and weighted by LD coefficient between SNPs, GATES rapidly combine the p values of SNPs within a gene and generate a gene-based p value by dividing independent p values of effective number among the SNPs within a gene by independent p values of the effective number among the top SNPs.

Pathway enrichment analysis using KGG

In order to further understand and interpret the association results from our study, we explored the biological pathway enriched by the genes associated with WM microstructure abnormality using hybrid set-based test (HYST) combining the extended Simes’ test and scaled chi-square test in KGG platform36. To reduce type I error, a Bonferroni adjustment for multiple testing was conducted using a significant threshold of 0.05/N, where N = number of gene sets in GSEA (gene set enrichment analysis) pathway database (0.05/17779 ~ = 2.81×10–6).

Results

Demographic characteristics

All participants were Han Chinese. As shown in Table 1, there was no significant difference observed in age, years of education and gender between the patients and healthy controls. The median duration of illness for patients was 5.5 months (SD: 26.7). All patients were drug-naïve at the time of clinical evaluation and MRI scanning.

Table 1
Summary of demographic and clinic characteristics of subjects.

Comparison of white matter microstructure in patients with schizophrenia and controls

Two hundred and twenty-seven subjects (93 cases and 134 controls) had been scanned by MRI successfully. As shown in Figure 1, patients had reduced FA values in 5 brain regions, i.e. left anterior cingulate cortex (ACC-L), right anterior cingulate cortex (ACC-R), left inferior parietal cortex (IPC-L), left posterior cingulate cortex (PCC-L) and right posterior cingulate cortex (PCC-R), in comparison with healthy controls. The differences survived after FDR multiple corrections. The total mean FA value, which was calculated by combining and averaging FA values of all the brain regions, was significantly lower in patients compared to healthy controls. Therefore, we included total mean FA values along with FA values of 5 brain regions as QTs in subsequent association analysis.

Figure 1
DTI analysis revealing reductions of fractional anisotropy (FA) in schizophrenia patients compared to healthy controls in left anterior cingulate cortex, right anterior cingulate cortex, left inferior parietal cortex, left posterior cingulate cortex and ...

GWAS analysis of WM microstructural abnormalities in schizophrenia

In total, 227 subjects with high-quality genotypes (6,055,918 SNPs) remained and there was no notable population stratification within our sample. Following univariate association analysis, we found an association of a genotyped SNP (rs11901793) within 6.2 kb upstream of CXCR7 gene on chromosome 2 with the total mean FA value that passed the genome-wide threshold significance (p = 4.37 × 10–8). The other genotyped SNP, rs10509852, which is located in the intergenic region within 400 kb downstream of SORCS1 on chromosome 10 showed the second strongest signal after association analysis with FA values of PCC-L (p = 5.21 × 10−8). The regional plottings of SNP with the strongest association for each of the six traits were displayed in Figure Figure22 and the summary of results are shown in Supplementary Table 1. Accordingly, our power was 76% for an effect size of 0.41, at an allele frequency of 0.4237. Besides, gene-gene interaction analysis using MB-MDR observed a significant high-risk effect of interaction between rs10509852 (SORCS1) and rs11901793 (CXCR7) on the FA values in PCC-L in both cases (Wald statistic = 14.45, permutation p = 0.008) and controls (Wald statistic = 14.38, permutation p = 0.008) after 1000 permutations, highlighting the involvement of epistasis between these two genes in white matter microstructure in PCC-L independent of diagnosis.

Figure 2
The regional association plot using web-interface of LOCUSZOOM, (a) ACC-L ; (b) ACC-R ; (c) IPC-L ; (d) PCC-L ; (e) PCC-R ; (f) total mean FA values.

Moreover, functional annotation of 6 top associated SNPs indicated that these SNPs played an important role in transcription factor binding and transcription-enhancing in the tissues related to neurodevelopment and immunology (Table (Table22).

Table 2
Functional annotation of the risk SNPs associated with WM microstructural abnormalities in schizophrenia.

Multivariate analysis and pathway enrichment

Besides, the same SNP, rs10509852, yielded the most significant global trait-based p value (p = 1.89 × 10−7) after combining p values from univariate association analysis of FA values from 5 brain regions (ACC-L, ACC-R, IPC-L, PCC-L, PCC-R) using TATES (Figure 3). Furthermore, the lambda (λ) score following association study indicated that all confounding factors regarding genotyping and population structure have been well-adjusted which is illustrated through Q-Q plot in Supplementary Figure 1. The gene-based test using GATES showed that 8 genes (Table 3) consisting of four non-coding micro-RNA genes (miRNA-27a, miRNA24-2, LOC284454, and miRNA-23a) and four protein-coding genes (TEP1, PDZD9, MPP4, and UQCRC2) passed the significance threshold (p < 0.05), details of these genes are shown in Supplementary Table 2.

Figure 3
Manhattan plot of TATES test, green highlighted points are SNPs located within 400 kb of SORCS1 gene.
Table 3
Gene-based association using GATE.

Furthermore, hybrid set-based test (HYST) in KGG identified a cell cycle pathway (REACTOME_CHROMOSOME_MAINTENANCE, p = 8.68 × 10−22) to be significantly enriched and remained significant after the Bonferroni correction for multiple comparisons (Corrected p = 1.54 × 10−17). Figure Figure44 displayed the Reactome F1 network and showed that out of 104 genes constituting the pathway, 38 genes have nominally significant gene-wise p values in our study (Details of 38 genes are displayed in Supplementary Table 2).

Figure 4
Reactome F1 network visualization of enriched pathway (REACTOME_CHROMOSOME_MAINTENANCE), nodes with thick purple circle in diagram were hit genes.

Discussion

To our knowledge, this is the first study to use FA values as QTs in GWAS of schizophrenia. In our study, after inter-group comparison, FA values in 5 brain regions (ACC-L, ACC-R, IPC-L, PCC-L, PCC-R) and the total mean FA values showed a significant difference between first-episode, drug-naïve schizophrenia patients and controls after multiple corrections. Moreover, using microstructural abnormality as QTs, we also found that two variants near/within CXCR7 and SORCS1 genes, eight genes (four non-coding micro-RNA genes including miRNA-27a, miRNA24-2, LOC284454 and miRNA-23a and four protein-coding genes including TEP1, PDZD9, MPP4 and UQCRC2), and one cell cycle pathway (REACTOME_CHROMOSOME_MAINTENANCE) were implicated in WM microstructural abnormality of schizophrenia in the present study.

Of the five brain regions, ACC-L, ACC-R, PCC-R, and PCC-L belong to the cingulum bundle (CB) which connects the limbic structure with frontal and temporal cortex. Previous studies have also identified the alterations of CB in schizophrenia. Kubicki et al. found that the mean FA values in bilateral CB in schizophrenic patients was decreased while compared to the controls38. Voineskos et al. examined the microstructural integrity of frontotemporal and interhemispheric WM tracts in schizophrenia across the adult lifespan, while younger patients had lower FA values in the right CB in relation to the younger controls, this difference disappeared in the older group of patients39. In addition, our study also found that FA values in the left inferior parietal lobe were significantly lower in patients than controls and the finding is consistent with those of Ardekani et al. and Rowland et al.40,41. It is worth noting that, patients with schizophrenia in our study showed significantly lower total mean FA values. In fact, schizophrenia has long been presumed to be the illness with comprehensive deficits of WM42,43 and our findings derived from the total mean FA values implies that focal WM deficit might appear based on the contextual whole-brain WM alterations.

In the present study, we found that an intergenic SNP (rs11901793) flanking CXCR7 gene (6.2 kb 3′) was associated with the total mean FA values with genome-wide significance. CXCR7 is a chemokine receptor which modulates cell migration in multiple biological contexts including brain development, especially in the signaling modulation of GABaergic interneurons4446. It has been found in various studies that the deficits in the development of GABaergic neurons were involved in the pathogenesis of schizophrenia, especially the cognitive impairment in schizophrenia4749. Moreover, an intergenic SNP (rs10509852) flanking SORCS1 gene on chromosome 10 was associated with both (PCC-L) and the global trait. SORCS1 is one of the transmembrane receptor family with a homologous domain to the yeast vacuolar sorting protein VPS10p (Sortilin) and it has been demonstrated to bind with the neuropeptides and neurotensins5052. Intriguingly, one study showed that Sortilin played an important role in BDNF sorting that regulated the secretory pathway by interacting specifically with BDNF in a region encompassing the methionine substitution, and colocalized with BDNF in secretory granules in neurons53. Indeed, both BDNF and neuronal cell death have been found to be involved in the pathophysiology of schizophrenia. Our findings provided further insights into understanding the mechanism of schizophrenia, especially the one associated with these seemingly diverse biological signals. Intriguingly, the majority of the top associated SNPs in our study was found to confer their impact through enhancing the different transcription factors in tissues critical of neurodevelopment and immunology, such as derived neuronal progenitor cultured cells (ES-deriv), cortex derived primary cultured nuerosphere (ES-deriv, Neurosph) and T-cell. Such a pattern shows another piece of evidence that neurodevelopmental and immunological pathways are implicated in the pathogenesis of schizophrenia. A future study could focus on the biological function of the protein-coding sequence to which these transcription factors bind and also on the effect of environmental factors on the function of these transcription factors, especially the chromatin states.

In our gene-based analysis, we found eight genes that were significantly associated with the global trait of FA values from five correlated brain regions (ACC-L, ACC-R, IPC-L, PCC-L, PCC-R) and survived the multiple corrections after pooling effect size of all variants of each gene together. Of these eight genes, TEP1, telomerase-associated protein 1, is responsible for catalyzing the addition of new telomeres to chromosomes and prevents developing neurons from DNA damage-induced cell death54,55. MPP4, membrane protein, palmitoylated 4, functions as scaffold proteins that contribute to the cell polarity and organize signal transducers at the neuronal synapse membrane56. Although it remains unclear on the pathophysiological role of these two genes, our study provided evidence to support the involvement of these two genes in schizophrenia. Future research on the biological functions of genes should be justified. Another four genes which were significantly associated with the abnormal white matter microstructure in schizophrenia at the gene-based level encode the non-coding micro-RNA (miRNA-27a, miRNA24-2, LOC284454 and miRNA23a). In recent years, micro-RNAs have gained increasing attention in the etiology of schizophrenia. In fact, over half of the miRNAs identified so far have been shown to be highly or exclusively expressed in the brain, and some of them have been found to be abnormally expressed in patients with schizophrenia57,58. Of interest, one of the important regulatory targets of miRNA23a gene is CXCL12, of which CXCR7 is the main receptor59.

Furthermore, our pathway analysis identified an enriched cell cycle pathway, REACTOME_CHROMOSOME_MAINTENANCE, which was associated with WM microstructure abnormality in schizophrenia. The previous study has shown that accurate neuronal cell cycle is one of the most critical components for neuronal development, deviation of which can cause premature cell cycle exit60,61. Using unbiased molecular profiling analysis, Wang et al. obtained skin fibroblasts from the living patients with schizophrenia and identified perturbations in cell cycle proteins. They further validated their results at expression and functional levels in the same study62. Moreover, for 38 genes enriched in this pathway, some of them have already been found to function in neurodevelopment and our results provided additional evidence for supporting the role of cell cycle pathway in schizophrenia, but further studies were still required to validate our present findings.

Although the results of our study are promising and reliable, we fully acknowledge that there are a few limitations. Firstly, we conducted studies in schizophrenic patients as a whole and did not further stratify the patients based on their clinical features. Although the traits related to the disease as a whole can avoid the confounding factors arisen from the subjective categorization of clinical symptomology, future studies should focus on the genetic basis of WM microstructure abnormality in different subgroups of schizophrenia. Secondly, due to the complexity of schizophrenia and modest sample size of our study, any definitive conclusions should be precluded, and future studies in a larger sample are needed to validate our findings.

In conclusion, our results suggest that WM abnormalities of schizophrenia patients mainly present in CB region and parietal cortex. These abnormalities are associated with the genes that are likely to be involved in diverse biological signals and highly enriched in the biological pathway of cell cycle. In addition, our study suggests that the strategy of using neuroimaging measures as QTs of schizophrenia could significantly improve the power of GWAS to identify the susceptible genes/pathways of schizophrenia.

Electronic supplementary material

Acknowledgements

This work was partly funded by National Nature Science Foundation of China Key Project 81630030 and 81130024 (to T.Li.); National Natural Science Foundation of China/Research Grants Council of Hong Kong Joint Research Scheme 81461168029 (to T.Li.), National Key Research and Development Program of the Ministry of Science and Technology of China 2016YFC0904300 (to T.Li.).

Author Contributions

Author Contributions

H.Y.R. and Q.W. conducted the genome-wide association analysis, the gene-set analysis, the multivariate analysis and wrote the main manuscript; W.L. conducted FA comparisons between the patients and the controls; C.C.Z., Y.F.L., X.J.L., C.H.H. and M.L.L. collected data including DNA samples and neuroimaging data; W.D., L.S.Z., and Y.C.W. extracted DNA from the blood samples; X.H.M., X.H., and F.D. supervised the study; T.L. designed and supervised the study.

Notes

Competing Interests

The authors declare that they have no competing interests.

Footnotes

H. Y. Ren and Q. Wang contributed equally to this work.

Electronic supplementary material

Supplementary information accompanies this paper at doi:10.1038/s41598-017-10507-7

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

X. Hu, moc.361@uhxxh.

T. Li, nc.ude.ucs@xhoatil.

References

1. Torrey EF. Prevalence studies in schizophrenia. The British Journal of Psychiatry. 1987;150:598–608. doi: 10.1192/bjp.150.5.598. [PubMed] [Cross Ref]
2. McGuffin P, Farmer AE, Gottesman II, Murray RM, Reveley AM. Twin concordance for operationally defined schizophrenia: confirmation of familiality and heritability. Archives of General Psychiatry. 1984;41:541–545. doi: 10.1001/archpsyc.1984.01790170015002. [PubMed] [Cross Ref]
3. Sawa A, Snyder SH. Schizophrenia: diverse approaches to a complex disease. Science. 2002;296:692–695. doi: 10.1126/science.1070532. [PubMed] [Cross Ref]
4. Consortium SPG-WAS. Genome-wide association study identifies five new schizophrenia loci. Nature genetics. 2011;43:969–976. doi: 10.1038/ng.940. [PMC free article] [PubMed] [Cross Ref]
5. O’donovan MC, et al. Identification of loci associated with schizophrenia by genome-wide association and follow-up. Nature genetics. 2008;40:1053–1055. doi: 10.1038/ng.201. [PubMed] [Cross Ref]
6. Yue W-H, et al. Genome-wide association study identifies a susceptibility locus for schizophrenia in Han Chinese at 11p11. 2. Nature genetics. 2011;43:1228–1231. doi: 10.1038/ng.979. [PubMed] [Cross Ref]
7. Stefansson H, et al. Common variants conferring risk of schizophrenia. Nature. 2009;460:744–747. [PMC free article] [PubMed]
8. Ripke S, et al. Biological insights from 108 schizophrenia-associated genetic loci. Nature. 2014;511:421. doi: 10.1038/nature13595. [PMC free article] [PubMed] [Cross Ref]
9. Manolio TA, et al. Finding the missing heritability of complex diseases. Nature. 2009;461:747–753. doi: 10.1038/nature08494. [PMC free article] [PubMed] [Cross Ref]
10. Almasy L, Blangero J. Endophenotypes as quantitative risk factors for psychiatric disease: rationale and study design. American journal of medical genetics. 2001;105:42–44. doi: 10.1002/1096-8628(20010108)105:1<42::AID-AJMG1055>3.0.CO;2-9. [PubMed] [Cross Ref]
11. Plomin R, Haworth CM, Davis OS. Common disorders are quantitative traits. Nature Reviews Genetics. 2009;10:872–878. doi: 10.1038/nrg2670. [PubMed] [Cross Ref]
12. Potkin SG, et al. A genome-wide association study of schizophrenia using brain activation as a quantitative phenotype. Schizophrenia bulletin. 2009;35:96–108. doi: 10.1093/schbul/sbn155. [PMC free article] [PubMed] [Cross Ref]
13. Dickinson D, et al. Differential effects of common variants in SCN2A on general cognitive ability, brain physiology, and messenger RNA expression in schizophrenia cases and control individuals. JAMA psychiatry. 2014;71:647–656. doi: 10.1001/jamapsychiatry.2014.157. [PMC free article] [PubMed] [Cross Ref]
14. Wang Q, et al. Genome-wide association analysis with gray matter volume as a quantitative phenotype in first-episode treatment-naive patients with schizophrenia. PloS one. 2013;8:e75083. doi: 10.1371/journal.pone.0075083. [PMC free article] [PubMed] [Cross Ref]
15. Wang Q, et al. The CHRM3 gene is implicated in abnormal thalamo-orbital frontal cortex functional connectivity in first-episode treatment-naive patients with schizophrenia. Psychological medicine. 2016;46:1523. doi: 10.1017/S0033291716000167. [PubMed] [Cross Ref]
16. Camchong J, Lim KO, Sponheim SR, MacDonald AW., III Frontal white matter integrity as an endophenotype for schizophrenia: diffusion tensor imaging in monozygotic twins and patients’ nonpsychotic relatives. Frontiers in human neuroscience. 2009;3:35. doi: 10.3389/neuro.09.035.2009. [PMC free article] [PubMed] [Cross Ref]
17. Bertisch H, Li D, Hoptman MJ, DeLisi LE. Heritability estimates for cognitive factors and brain white matter integrity as markers of schizophrenia. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics. 2010;153:885–894. doi: 10.1002/ajmg.b.31054. [PMC free article] [PubMed] [Cross Ref]
18. Zeng B, et al. Abnormal white matter microstructure in drug-naive first episode schizophrenia patients before and after eight weeks of antipsychotic treatment. Schizophrenia research. 2016;172:1–8. doi: 10.1016/j.schres.2016.01.051. [PubMed] [Cross Ref]
19. Ellison-Wright I, Bullmore E. Meta-analysis of diffusion tensor imaging studies in schizophrenia. Schizophrenia research. 2009;108:3–10. doi: 10.1016/j.schres.2008.11.021. [PubMed] [Cross Ref]
20. Nesvåg R, et al. Regional thinning of the cerebral cortex in schizophrenia: effects of diagnosis, age and antipsychotic medication. Schizophrenia research. 2008;98:16–28. doi: 10.1016/j.schres.2007.09.015. [PubMed] [Cross Ref]
21. Mallas, E. et al. The impact of CACNA1C gene, and its epistasis with ZNF804A, on white matter microstructure in health, schizophrenia and bipolar disorder(1). Genes, brain, and behavior, doi:10.1111/gbb.12355 (2016). [PubMed]
22. Mounce J, et al. Association of GRM3 polymorphism with white matter integrity in schizophrenia. Schizophrenia research. 2014;155:8–14. doi: 10.1016/j.schres.2014.03.003. [PMC free article] [PubMed] [Cross Ref]
23. Chen Y, Norton D, Stromeyer C. Prolonged temporal interaction for peripheral visual processing in schizophrenia: evidence from a three-flash illusion. Schizophrenia research. 2014;156:190–196. doi: 10.1016/j.schres.2014.04.018. [PMC free article] [PubMed] [Cross Ref]
24. Purcell S, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. American journal of human genetics. 2007;81:559–575. doi: 10.1086/519795. [PubMed] [Cross Ref]
25. Price AL, et al. Principal components analysis corrects for stratification in genome-wide association studies. Nature genetics. 2006;38:904–909. doi: 10.1038/ng1847. [PubMed] [Cross Ref]
26. Howie BN, Donnelly P, Marchini J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS genetics. 2009;5:e1000529. doi: 10.1371/journal.pgen.1000529. [PMC free article] [PubMed] [Cross Ref]
27. Delaneau O, Marchini J, Zagury JF. A linear complexity phasing method for thousands of genomes. Nat Methods. 2012;9:179–181. doi: 10.1038/nmeth.1785. [PubMed] [Cross Ref]
28. Aulchenko YS, Ripke S, Isaacs A, van Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007;23:1294–1296. doi: 10.1093/bioinformatics/btm108. [PubMed] [Cross Ref]
29. Kundaje A, et al. Integrative analysis of 111 reference human epigenomes. Nature. 2015;518:317–330. doi: 10.1038/nature14248. [PMC free article] [PubMed] [Cross Ref]
30. Ernst J, et al. Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011;473:43–49. doi: 10.1038/nature09906. [PMC free article] [PubMed] [Cross Ref]
31. Ward LD, Kellis M. HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic acids research. 2011;40:D930–D934. doi: 10.1093/nar/gkr917. [PMC free article] [PubMed] [Cross Ref]
32. Boyle AP, et al. Annotation of functional variation in personal genomes using RegulomeDB. Genome research. 2012;22:1790–1797. doi: 10.1101/gr.137323.112. [PubMed] [Cross Ref]
33. Calle ML, Urrea V, Malats N, Van Steen K. mbmdr: an R package for exploring gene–gene interactions associated with binary or quantitative traits. Bioinformatics. 2010;26:2198–2199. doi: 10.1093/bioinformatics/btq352. [PubMed] [Cross Ref]
34. van der Sluis, S., Posthuma, D. & Dolan, C. V. TATES: Efficient Multivariate Genotype-Phenotype Analysis for Genome-Wide Association Studies. PLoS genetics9, doi:ARTN e100323510.1371/journal.pgen.1003235 (2013). [PMC free article] [PubMed]
35. Li MX, Gui HS, Kwan JSH, Sham PC. GATES: A Rapid and Powerful Gene-Based Association Test Using Extended Simes Procedure. American journal of human genetics. 2011;88:283–293. doi: 10.1016/j.ajhg.2011.01.019. [PubMed] [Cross Ref]
36. Li M-X, Kwan JS, Sham PC. HYST: a hybrid set-based test for genome-wide association studies, with application to protein-protein interaction-based association analysis. The American Journal of Human Genetics. 2012;91:478–488. doi: 10.1016/j.ajhg.2012.08.004. [PubMed] [Cross Ref]
37. Faul F, Erdfelder E, Lang A-G, Buchner A. G* Power 3: A flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behavior research methods. 2007;39:175–191. doi: 10.3758/BF03193146. [PubMed] [Cross Ref]
38. Kubicki M, et al. Cingulate fasciculus integrity disruption in schizophrenia: a magnetic resonance diffusion tensor imaging study. Biological psychiatry. 2003;54:1171–1180. doi: 10.1016/S0006-3223(03)00419-0. [PMC free article] [PubMed] [Cross Ref]
39. Voineskos AN, et al. Diffusion tensor tractography findings in schizophrenia across the adult lifespan. Brain: a journal of neurology. 2010;133:1494–1504. doi: 10.1093/brain/awq040. [PMC free article] [PubMed] [Cross Ref]
40. Ardekani BA, Nierenberg J, Hoptman MJ, Javitt DC, Lim KO. MRI study of white matter diffusion anisotropy in schizophrenia. Neuroreport. 2003;14:2025–2029. doi: 10.1097/00001756-200311140-00004. [PubMed] [Cross Ref]
41. Rowland LM, et al. White matter alterations in deficit schizophrenia. Neuropsychopharmacology: official publication of the American College of Neuropsychopharmacology. 2009;34:1514–1522. doi: 10.1038/npp.2008.207. [PMC free article] [PubMed] [Cross Ref]
42. Hao Y, et al. White matter integrity of the whole brain is disrupted in first-episode schizophrenia. Neuroreport. 2006;17:23–26. doi: 10.1097/01.wnr.0000195664.15090.46. [PubMed] [Cross Ref]
43. Liu Y, et al. Disrupted small-world networks in schizophrenia. Brain: a journal of neurology. 2008;131:945–961. doi: 10.1093/brain/awn018. [PubMed] [Cross Ref]
44. Sánchez-Alcañiz JA, et al. Cxcr7 controls neuronal migration by regulating chemokine responsiveness. Neuron. 2011;69:77–90. doi: 10.1016/j.neuron.2010.12.006. [PubMed] [Cross Ref]
45. Wang Y, et al. CXCR4 and CXCR7 have distinct functions in regulating interneuron migration. Neuron. 2011;69:61–76. doi: 10.1016/j.neuron.2010.12.005. [PMC free article] [PubMed] [Cross Ref]
46. Zhu B, et al. CXCL12 Enhances Human Neural Progenitor Cell Survival Through a CXCR7‐and CXCR4‐Mediated Endocytotic Signaling Pathway. Stem Cells. 2012;30:2571–2583. doi: 10.1002/stem.1239. [PMC free article] [PubMed] [Cross Ref]
47. Hashimoto T, et al. Gene expression deficits in a subclass of GABA neurons in the prefrontal cortex of subjects with schizophrenia. The Journal of Neuroscience. 2003;23:6315–6326. [PubMed]
48. Fung, S. J. et al. Expression of interneuron markers in the dorsolateral prefrontal cortex of the developing human and in schizophrenia. American Journal of Psychiatry (2010). [PubMed]
49. Volk, D. W. & Lewis, D. A. Early developmental disturbances of cortical inhibitory neurons: contribution to cognitive deficits in schizophrenia. Schizophrenia bulletin, sbu111 (2014). [PMC free article] [PubMed]
50. Marcusson EG, Horazdovsky BF, Cereghino JL, Gharakhanian E, Emr SD. The sorting receptor for yeast vacuolar carboxypeptidase Y is encoded by the VPS10 gene. Cell. 1994;77:579–586. doi: 10.1016/0092-8674(94)90219-4. [PubMed] [Cross Ref]
51. Jacobsen L, et al. Activation and functional characterization of the mosaic receptor SorLA/LR11. Journal of Biological Chemistry. 2001;276:22788–22796. doi: 10.1074/jbc.M100857200. [PubMed] [Cross Ref]
52. Mazella J, et al. The 100-kDa neurotensin receptor is gp95/sortilin, a non-G-protein-coupled receptor. Journal of Biological Chemistry. 1998;273:26273–26276. doi: 10.1074/jbc.273.41.26273. [PubMed] [Cross Ref]
53. Chen Z-Y, et al. Sortilin controls intracellular sorting of brain-derived neurotrophic factor to the regulated secretory pathway. The Journal of neuroscience. 2005;25:6156–6166. doi: 10.1523/JNEUROSCI.1017-05.2005. [PMC free article] [PubMed] [Cross Ref]
54. Lu C, Fu W, Mattson MP. Telomerase protects developing neurons against DNA damage-induced cell death. Developmental Brain Research. 2001;131:167–171. doi: 10.1016/S0165-3806(01)00237-1. [PubMed] [Cross Ref]
55. Yan L, Wu S, Zhang S, Ji G, Gu A. Genetic variants in telomerase reverse transcriptase (TERT) and telomerase-associated protein 1 (TEP1) and the risk of male infertility. Gene. 2014;534:139–143. doi: 10.1016/j.gene.2013.11.008. [PubMed] [Cross Ref]
56. Aartsen WM, et al. Mpp4 recruits Psd95 and Veli3 towards the photoreceptor synapse. Human molecular genetics. 2006;15:1291–1302. doi: 10.1093/hmg/ddl047. [PubMed] [Cross Ref]
57. Perkins DO, et al. microRNA expression in the prefrontal cortex of individuals with schizophrenia and schizoaffective disorder. Genome biology. 2007;8:R27. doi: 10.1186/gb-2007-8-2-r27. [PMC free article] [PubMed] [Cross Ref]
58. Miller BH, et al. MicroRNA-132 dysregulation in schizophrenia has implications for both neurodevelopment and adult brain function. Proceedings of the National Academy of Sciences. 2012;109:3125–3130. doi: 10.1073/pnas.1113793109. [PubMed] [Cross Ref]
59. Arabanian LS, et al. miRNA-23a mediates post-transcriptional regulation of CXCL12 in bone marrow stromal cells. haematologica, haematol. 2014;2013:097675. [PubMed]
60. Ohnuma S-i, Philpott A, Harris WA. Cell cycle and cell fate in the nervous system. Current opinion in neurobiology. 2001;11:66–73. doi: 10.1016/S0959-4388(00)00175-6. [PubMed] [Cross Ref]
61. Nguyen L, Besson A, Roberts JM, Guillemot F. Coupling cell cycle exit, neuronal differentiation and migration in cortical neurogenesis. Cell Cycle. 2006;5:2314–2318. doi: 10.4161/cc.5.20.3381. [PubMed] [Cross Ref]
62. Wang L, et al. Expression profiling of fibroblasts identifies cell cycle abnormalities in schizophrenia. Journal of proteome research. 2009;9:521–527. doi: 10.1021/pr900867x. [PubMed] [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group