Hostname: page-component-7bb8b95d7b-l4ctd Total loading time: 0 Render date: 2024-10-04T15:15:45.900Z Has data issue: false hasContentIssue false

Polygenic effects of schizophrenia on hippocampal grey matter volume and hippocampus–medial prefrontal cortex functional connectivity

Published online by Cambridge University Press:  06 June 2019

Shu Liu
Affiliation:
MSc Student, Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences School of Artificial Intelligence, University of Chinese Academy of Sciences, China
Ang Li
Affiliation:
School of Artificial Intelligence, University of Chinese Academy of Sciences, China PhD Student, Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences
Yong Liu
Affiliation:
School of Artificial Intelligence, University of Chinese Academy of Sciences, China Professor, Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences
Hao Yan
Affiliation:
Associate Professor, Peking University Sixth Hospital, Institute of Mental Health Key Laboratory of Mental Health, Ministry of Health (Peking University), China
Meng Wang
Affiliation:
School of Artificial Intelligence, University of Chinese Academy of Sciences, China PhD Student, Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences
Yuqing Sun
Affiliation:
School of Artificial Intelligence, University of Chinese Academy of Sciences, China PhD Student, Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences
Lingzhong Fan
Affiliation:
School of Artificial Intelligence, University of Chinese Academy of Sciences, China Professor, Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences
Ming Song
Affiliation:
School of Artificial Intelligence, University of Chinese Academy of Sciences, China Associate Professor, Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences
Kaibin Xu
Affiliation:
School of Artificial Intelligence, University of Chinese Academy of Sciences, China PhD Student, Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences
Jun Chen
Affiliation:
Associate Professor, Department of Radiology, Renmin Hospital of Wuhan University, China
Yunchun Chen
Affiliation:
Associate Professor, Department of Psychiatry, Xijing Hospital, The Fourth Military Medical University, China
Huaning Wang
Affiliation:
Associate Professor, Department of Psychiatry, Xijing Hospital, The Fourth Military Medical University, China
Hua Guo
Affiliation:
Professor, Zhumadian Psychiatric Hospital, China
Ping Wan
Affiliation:
Professor, Zhumadian Psychiatric Hospital, China
Luxian Lv
Affiliation:
Professor, Department of Psychiatry, Henan Mental Hospital, The Second Affiliated Hospital of Xinxiang Medical University Henan Key Lab of Biological Psychiatry, Xinxiang Medical University, China
Yongfeng Yang
Affiliation:
Henan Key Lab of Biological Psychiatry, Xinxiang Medical University, China Attending Doctor, Department of Psychiatry, Henan Mental Hospital, The Second Affiliated Hospital of Xinxiang Medical University
Peng Li
Affiliation:
Key Laboratory of Mental Health, Ministry of Health (Peking University), China Associate Professor, Peking University Sixth Hospital, Institute of Mental Health
Lin Lu
Affiliation:
Key Laboratory of Mental Health, Ministry of Health (Peking University), China Professor, Peking University Sixth Hospital, Institute of Mental Health
Jun Yan
Affiliation:
Key Laboratory of Mental Health, Ministry of Health (Peking University), China Professor, Peking University Sixth Hospital, Institute of Mental Health
Huiling Wang
Affiliation:
Professor, Department of Radiology, Renmin Hospital of Wuhan University, China
Hongxing Zhang
Affiliation:
Henan Key Lab of Biological Psychiatry, Xinxiang Medical University, China Professor, Department of Psychiatry, Henan Mental Hospital, The Second Affiliated Hospital of Xinxiang Medical University
Huawang Wu
Affiliation:
Attending Doctor, Guangzhou Brain Hospital, The Affiliated Brain Hospital of Guangzhou Medical University, China
Yuping Ning
Affiliation:
Professor, Guangzhou Brain Hospital, The Affiliated Brain Hospital of Guangzhou Medical University, China
Dai Zhang
Affiliation:
Key Laboratory of Mental Health, Ministry of Health (Peking University), China Professor, Peking University Sixth Hospital, Institute of Mental Health
Tianzi Jiang
Affiliation:
School of Artificial Intelligence, University of Chinese Academy of Sciences, China Professor, Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences
Bing Liu*
Affiliation:
School of Artificial Intelligence, University of Chinese Academy of Sciences, China Professor, Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences
*
Correspondence: Bing Liu, Brainnetome Center, Institute of Automation, Chinese Academy of Sciences, 95 Zhong Guan Cun East Road, Hai Dian District, Beijing100190, China. Email: bliu@nlpr.ia.ac.cn
Rights & Permissions [Opens in a new window]

Abstract

Background

Schizophrenia is a complex mental disorder with high heritability and polygenic inheritance. Multimodal neuroimaging studies have also indicated that abnormalities of brain structure and function are a plausible neurobiological characterisation of schizophrenia. However, the polygenic effects of schizophrenia on these imaging endophenotypes have not yet been fully elucidated.

Aims

To investigate the effects of polygenic risk for schizophrenia on the brain grey matter volume and functional connectivity, which are disrupted in schizophrenia.

Method

Genomic and neuroimaging data from a large sample of Han Chinese patients with schizophrenia (N = 509) and healthy controls (N = 502) were included in this study. We examined grey matter volume and functional connectivity via structural and functional magnetic resonance imaging, respectively. Using the data from a recent meta-analysis of a genome-wide association study that comprised a large number of Chinese people, we calculated a polygenic risk score (PGRS) for each participant.

Results

The imaging genetic analysis revealed that the individual PGRS showed a significantly negative correlation with the hippocampal grey matter volume and hippocampus–medial prefrontal cortex functional connectivity, both of which were lower in the people with schizophrenia than in the controls. We also found that the observed neuroimaging measures showed weak but similar changes in unaffected first-degree relatives of patients with schizophrenia.

Conclusions

These findings suggested that genetically influenced brain grey matter volume and functional connectivity may provide important clues for understanding the pathological mechanisms of schizophrenia and for the early diagnosis of schizophrenia.

Type
Papers
Copyright
Copyright © The Royal College of Psychiatrists 2019

Studies have shown that schizophrenia involves a highly heritable form of psychosis1,Reference Ripke, O'Dushlaine, Chambert, Moran, Kähler and Akterin2 and that the effect of each identified common genetic variant is very small.Reference Ripke, O'Dushlaine, Chambert, Moran, Kähler and Akterin2 Furthermore, the effects of specific genetic variants on the clinical diagnosis are sometimes inconsistent across studies.Reference Kishi, Fukuo, Okochi, Kawashima, Kitajima and Inada3 This may be because the phenotypic characteristics of schizophrenia that are modulated by polygenic factors are heterogeneous and complex.Reference Purcell, Wray, Stone, Visscher, O'Donovan and Sullivan4 To capture the polygenic nature of complex disorders, a polygenic risk score (PGRS) – which is based on the additive effects of a great number of weak disease-related genetic variantsReference Purcell, Wray, Stone, Visscher, O'Donovan and Sullivan4 – can be used. In this way, the PGRS can effectively assess the effects of these disease-related genetic variants on various phenotypes. Neuroimaging studies have provided sufficient evidence that there are substantial structural and functional brain abnormalities in people with schizophrenia. A systematic meta-analysis found reduced grey matter volume (GMV) in people with schizophrenia in many cortical and subcortical regions such as the frontal lobe, temporal lobe, anterior cingulate and hippocampus.Reference Shepherd, Laurens, Matheson, Carr and Green5 Impaired functional activities and connectivity in some brain regions have also been widely reported. For example, a study found that the functional connectivity between the hippocampus and prefrontal cortex was impaired in people with schizophrenia.Reference Benetti, Mechelli, Picchioni, Broome, Williams and McGuire6 However, whether and how the changes in these brain measures are genetically influenced remains unclear.

Previous studies have reported that the schizophrenia PGRS is associated with many types of neuroimaging endophenotypes, such as cortical gyrificationReference Liu, Zhang, Cui, Qin, Tao and Li7 and working memory-related brain activity.Reference Walton, Geisler, Lee, Hass, Turner and Liu8Reference Walton, Turner, Gollub, Manoach, Yendiki and Ho10 However, to the best of our knowledge few studies have reported the abnormal GMV and functional connectivity that are associated with the PGRS for schizophrenia,Reference Wang, Liu, Zhang, Cui, Yu and Jiang11 so our aim was to investigate whether and how schizophrenia polygenic risk variants influence these brain measures. All participants were of Chinese Han ancestry, so PGRSs were derived from a meta-analysis of a genome-wide association study (GWAS) comprising a large number of Chinese people.Reference Li, Chen, Yu, He, Xu and Zhang12 We investigated the impaired GMV and functional connectivity in people with schizophrenia compared with healthy controls and further examined the effects of schizophrenia polygenic risk on these brain measures. Finally, we collected data from unaffected first-degree relatives (who are at high genetic risk for schizophrenia) to find further evidence for the identified biomarkers related to the PGRS.

Methods

Participants

A total of 1011 unrelated and extensively evaluated people of Chinese Han ancestry (509 people with schizophrenia and 502 healthy controls) with both high-quality genetic and structural magnetic resonance imaging (MRI) neuroimaging data were included in this study (Table 1). They were recruited from six hospitals (Peking University Sixth Hospital, Beijing Huilongguan Hospital, Xijing Hospital, Henan Mental Hospital, Renmin Hospital of Wuhan University and Zhumadian Psychiatric Hospital). In accordance with the DSM-IV (1994) criteria for schizophrenia, all the people with schizophrenia were diagnosed consistently by two veteran senior psychiatrists using the Structured Clinical Interview for DSM-IV-TR Axis I Disorders (SCID-I, patient edition). In the Positive and Negative Syndrome Scale (PANSS) test, the total scores of all the people with schizophrenia were higher than 60, and the scores in at least 3 positive items were higher than 4. The exclusion criteria for the patients included: other psychiatric and cognitive disorders at the time of the study or by history; severe physical diseases, such as diabetes, thyroid diseases, hypertension or cardiac diseases; a history of epilepsy; a DSM-IV diagnosis of alcohol or drug dependence; electroconvulsive therapy within the past 6 months; a suicide attempt or symptoms of severe excitement and agitation; and being pregnant or breastfeeding. Healthy controls were also carefully screened by the SCID-I, non-patient edition. Furthermore, the healthy participants were also excluded if they or their first- or second-degree relatives had ever been diagnosed with any psychiatric disorders. After receiving a complete description of the study, all the participants and/or their legal guardians provided written informed consent. The project was approved by the Medical Research Ethics Committee of Peking University Sixth Hospital (number 2010052). The Chinese revision of the Wechsler Adult Intelligence Scale test, including the Symbol Digital Modalities Test (SDMT), Forward Digital Span (forward-DS) and Backward Digital Span (backward-DS) tests (Table 1), was also used to assess the cognitive ability of both the people with schizophrenia and the healthy controls.

Table 1 Demographic and clinical characteristics

Values are means (s.d.) if nothing else is specified. All group comparisons were made with t-test except for gender (χ2). SDMT, Symbol Digital Modalities Test; forward-DS, Forward Digital Span; backward-DS, Backward Digital Span; PANSS, Positive and Negative Syndrome Scale; CPZ-eq, chlorpromazine equivalents.

a. The medication dosage was converted to CPZ-eq according to previous references.

To further explore the clinical significance of the endophenotypes we found, we also recruited 95 unaffected first-degree relatives (39 males, mean age 34.72 ± 9.57, age range 18.17–59.75) of the participants with schizophrenia and collected genomic and imaging data. Apart from being a first-degree relative of a person diagnosed with schizophrenia, the other inclusion and exclusion criteria were the same as those of the healthy controls.

Genotype processing

We used the EZgene Blood gDNA Miniprep Kit to extract genomic DNA from the whole blood of each participant and carried out whole-genome genotyping on Illumina Human OmniZhongHua-8 BeadChips. We then performed the genotype quality control on the Linux system using PLINK version 1.0.7.Reference Purcell, Neale, Todd-Brown, Thomas, Ferreira and Bender13 First, we removed the individuals whose missing genotype rates were greater than 0.05. We then used the pairwise identity by descent estimate to identify individuals who could possibly be related. For pairs of individuals who had more similar genotypes than we would have expected by chance in a random sample, we removed the one with the greater missing rate. Next, we used single nucleotide polymorphism (SNP)-level filtering to remove SNPs with missing genotype rates greater than 0.05, a minor allele frequency less than 0.01 and a significant departure from Hardy–Weinberg equilibrium (P < 0.001). To avoid the confounding effects of population stratification,Reference Walton, Geisler, Lee, Hass, Turner and Liu8 we performed a principal component analysis on the Linux system by using EIGENSTRAT 5.0.2Reference Patterson, Price and Reich14,Reference Price, Patterson, Plenge, Weinblatt, Shadick and Reich15 on a linkage disequilibrium pruned set of autosomal SNPs. This set was obtained by carrying out linkage disequilibrium pruning with PLINK (r 2 < 0.05) and removing five long-range linkage disequilibrium regions with the HapMap phase 3 reference data-set.Reference Thorisson, Smith, Krishnan and Stein16 After getting 10 principal components we excluded the outliers of the samples with >6 s.d. Ungenotyped SNPs were imputed on the Linux system using SHAPEIT version 2 (r790)Reference Delaneau, Marchini and Zagury17 and IMPUTE2Reference Howie, Donnelly and Marchini18 with the 1000 Genomes Phase 1 reference data-set. Further analyses focused on autosomal SNPs with imputation quality scores greater than 0.8.

Computation of PGRS

The PGRS for schizophrenia was created using the ‘score’ utility in PLINK. The PGRS computation method was developed by Purcell and colleagues,Reference Purcell, Wray, Stone, Visscher, O'Donovan and Sullivan4 as described in full detail in Walton and colleagues.Reference Walton, Turner, Gollub, Manoach, Yendiki and Ho10 The PGRS for each participant was calculated by summing the number of risk alleles weighted by the strength of the association of each SNP with schizophrenia. The strength of the association for each SNP was measured by the risk allele effect size (natural log of the odds ratio) reported by the meta-analysis of the GWAS comprising a large number of Chinese individuals.Reference Li, Chen, Yu, He, Xu and Zhang12 The participants included in our study were obviously independent of the participants from the meta-analysis. Overall, 17 PGRSs for each participant were obtained with the following different SNP inclusion thresholds: P T < 0.5, P T < 0.4, P T < 0.3, P T < 0.2, P T < 0.1, P T < 0.05, P T < 0.01, P T < 0.005, P T < 0.001, P T < 5 × 10−4, P T < 1 × 10−4, P T < 1 × 10−5, P T < 1 × 10−6, P T < 1 × 10−7, P T < 1 × 10−8, P T < 1 × 10−9 and P T < 1 × 10−10.

Image acquisition

Participants were scanned on eight 3.0 Tesla Magnetic Resonance Scanners. Specifically, four sites used Siemens Trim 3T scanners and the other four sites used GE SIGNA 3T scanners. No hardware or system upgrade was carried out on any of the scanners during the scanning period. A high-resolution, three-dimensional, T1-weighted brain volume (BRAVO) MRI sequence was performed using a protocol with a matrix size of 256 × 256, resolution of 1 × 1 mm2, inversion time of 1100 ms and slice thickness of 1 mm. A total of 192 sagittal slices were acquired on the Siemens scanners and 188 sagittal slices on the GE scanners. A total of 820 of the original participants, including 409 people with schizophrenia (209 males, mean age 27.78 ± 6.54, age range 17.50–45.41) and 411 healthy controls (197 males, mean age 27.88 ± 6.00, age range 17.75–44.58) had sufficient resting-state functional imaging data and were thus included for further imaging analyses. Resting-state functional imaging was performed using a single-shot, gradient-echo, echo-planar imaging (EPI) sequence that was sensitive to blood oxygen level-dependent contrast with the following parameters: repetition time 2000 ms, echo time 30 ms, flip angle 90°, matrix size 64 × 64, resolution of axial slice 3.4375 × 3.4375 mm2, slice thickness 4 mm, gap between slices 0.6 mm. Only participants from one site had a 6-minute scan (180 time points), whereas the participants from the other sites had an 8-minute scan (240 time points).

Image processing

During the calculation of the voxel-based morphometry (VBM), T1 images were processed using MATLAB for Windows with SPM8 (https://www.fil.ion.ucl.ac.uk/spm/software/spm8/) implemented in the VBM Toolbox 8 (http://dbm.neuro.uni-jena.de/vbm/). Except for applying the affine regularisation according to the International Consortium for Brain Mapping template for East Asian brains, we mainly applied standard routines and default parameters. The T1 images were finally segmented into grey matter, white matter and cerebrospinal fluid and then the grey matter images were smoothed with a 6 mm full width at half maximum isotropic Gaussian kernel.

The data-set was screened for artefacts, head motion (translation >3 mm or degree >3°), registration and normalisation quality. We discarded the first 10 time points of the EPI images. The preprocessing steps of the EPI images included: (a) slice timing correction; (b) head motion correction; (c) rigid-body registration of the T1 image to the EPI mean image; (d) spatial normalization of EPI images to Montreal Neurological Institute (MNI) standard space by using the T1 image and subsequent resampling to 3 × 3 × 3 mm3; (f) removing the noise from whole brain signals, head motions and linear trends; and (g) smoothing with a 6 mm Gaussian kernel. All the processes were performed using a MATLAB-based pipeline toolbox BRANT for Windows (https://github.com/kbxu/brant) in which slice timing, realign, co-register, normalise and smooth are functions called from SPM.

For the preprocessed EPI images, the voxel-wise functional connectivity to the whole brain from a specific region of interest was also computed with the BRANT toolbox. The functional connectivity map for each individual was obtained by computing the Pearson's correlation coefficient between the average time series in the region of interest and all voxels in the brain.

Statistical analysis

We performed a voxel-wise, two-sample t-test to measure the difference in the GMV between the people with schizophrenia and healthy controls across the whole brain, including age, gender and sites as covariates. Because the participants were from the same population, we did not include the principal components from population stratification as covariates. We focused on whether the impaired GMV was influenced by polygenic risk for schizophrenia, so we defined the regions that were significantly different between the two groups as a mask and then used a voxel-wise multiple regression model to investigate the association between the schizophrenia PGRS and the GMV. Age, gender, sites and diagnosis were again added as covariates. Next, the observed region was extracted and defined as the region of interest and the functional connectivity from the region of interest to the whole brain was computed. Similar statistical methods were applied to the functional connectivity analysis, including a two-sample t-test and a multiple regression model. We used SPM to complete the above analyses, and the multiple comparisons were corrected by the AlphaSim method using the Resting-State fMRI Data Analysis Toolkit for Windows (REST; http://restfmri.net/forum). Finally, a Pearson's correlation analysis was performed to test the association of the PANSS scores with the PGRS-related GMV and functional connectivity in the people with schizophrenia with age, gender and sites as covariates. Similarly, as with all the samples (schizophrenia group and healthy control group combined), we obtained the Pearson's correlation of the PGRS-related GMV and functional connectivity with the cognitive abilities, using age, gender, sites and diagnosis as covariates.

Results

PGRS differences between people with schizophrenia and healthy controls

The participants included 509 people with schizophrenia (275 males, mean age 27.78 ± 6.48, age range 17.50–45.41) and 502 healthy controls (255 males, mean age 27.69 ± 5.91, age range 17.75–44.58) (Table 1). We compared 17 PGRSs between the schizophrenia group and the healthy control group and found that the PGRS in the people with schizophrenia was significantly higher than that in the healthy controls under each threshold. In addition, the PGRS that was computed with a threshold of 0.1 showed the most significant difference between these two groups (Supplementary Fig. 1 available at https://doi.org/10.1192/bjp.2019.127 and Table 1). Therefore, this PGRS was then used in the rest of the analyses because it best explained the difference in genetic architecture between the patients and the controls.

Association between the schizophrenia PGRS and disrupted GMV

We first investigated the difference in GMV between the people with schizophrenia and the healthy controls. We found decreased GMV in many brain regions, including the thalamus, frontal lobe, cingulate cortex, basal ganglia and hippocampus (Fig. 1a), but observed no increased GMV in the patients compared with the healthy controls. We then analysed the association between the PGRS and the GMV in the brain regions that showed significant case–control differences in GMV and found that a smaller right hippocampal GMV (peak voxel MNI coordinate: x = 21, y = −13, z = −24) was significantly correlated with a higher PGRS (Fig. 1b).

Fig. 1 Association between the schizophrenia polygenic risk score (PGRS) and impaired grey matter volume (GMV). (a) A two-sample t-test of the GMV between schizophrenia and healthy control groups with AlphaSim correction (single voxel P < 0.005, the corrected threshold P < 0.05 and cluster size threshold >161 voxels), the negative t-value representing the GMV of this region was significantly decreased in people with schizophrenia. (b) Multiple regression analysis testing the association between PGRS and disrupted GMV with AlphaSim correction (single voxel P < 0.005, the corrected threshold P < 0.05 and cluster size threshold >132 voxels).

Association between the schizophrenia PGRS and disrupted functional connectivity

We selected the right hippocampus as the region of interest to calculate each individual's hippocampal functional connectivity to the whole brain and found similar functional connectivity patterns in the schizophrenia group and the healthy controls (Fig. 2a and Fig. 1b). Using a two-sample t-test with AlphaSim correction, the functional connectivities from the hippocampus to the thalamus, middle/superior frontal cortex, precuneus and middle temporal cortex were significantly stronger in the patients than in the controls. Conversely, the functional connectivities to the prefrontal cortex, superior temporal cortex, insula lobe and postcentral/precentral gyrus were significantly weaker in the patients (Fig. 2c). We then analysed the association between the PGRS and the disrupted functional connectivities and found a relatively weak correlation between the PGRS and the hippocampus–left/right medial prefrontal cortex (mPFC, peak voxel MNI coordinate: x = −12, y = 51, z = 9) functional connectivity with AlphaSim correction (Fig. 2d).

Fig. 2 Association between the schizophrenia polygenic risk score (PGRS) and disrupted functional connectivity. (a) Mean graph of the functional connectivity between the hippocampus and the whole brain in people with schizophrenia. (b) Mean graph of the functional connectivity between the hippocampus and the whole brain in the healthy control group. (c) Two-sample t-test results for the functional connectivity in the schizophrenia and healthy control groups with AlphaSim correction (single voxel P < 0.005, corrected threshold P < 0.05 and cluster size threshold >21 voxels), the negative t-value representing the functional connectivity was decreased in the people with schizophrenia and increased in the healthy controls. (d) Multiple regression analysis testing the association between the PGRS and the disrupted functional connectivity with AlphaSim correction (single voxel P < 0.01, the corrected threshold P < 0.05 and cluster size threshold >56 voxels). L, left; R, right.

PGRS, hippocampal GMV and hippocampus–mPFC functional connectivity in unaffected first-degree relatives

To further investigate whether the PGRS-related brain measures reflect fundamental neural circuits that are influenced by genetic factors and could be regarded as plausible biomarkers for the early identification of schizophrenia, we compared the values of the PGRS, hippocampal GMV and hippocampus–left/right mPFC functional connectivity of unaffected first-degree relatives with the healthy control group and the people with schizophrenia, using analysis of variance (Fig. 3). We found that the PGRS (P T < 0.1) for the unaffected first-degree relatives was significantly greater than the PGRS for the healthy control group and considerably smaller than the PGRS for the people with schizophrenia (Fig. 3a). In addition, we found that the hippocampal GMV and hippocampus–mPFC functional connectivity in the unaffected first-degree relatives were slightly greater than those in the people with schizophrenia, and both of them were significantly smaller than those in the healthy control group (Fig. 3, b and c). However, we did not find the significant association of the PGRS for schizophrenia with the hippocampal GMV and hippocampus–mPFC functional connectivity in the unaffected first-degree relatives.

Fig. 3 Polygenic risk score (PGRS), right hippocampal grey matter volume (GMV) and hippocampus–medial prefrontal cortex (mPFC) functional connectivity in the healthy control (HC) group, unaffected first-degree relatives (high-risk individuals, HR) and people with schizophrenia (SZ). (a) Distribution of the PGRS. (b) Distribution of the right hippocampal GMV. (c) Distribution of the hippocampus–mPFC functional connectivity (FC). *P < 0.05, **P < 0.01, ***P < 0.001.

Discussion

In this study we found that the PGRS in the people with schizophrenia was significantly greater than that in the controls. Our imaging genetic analysis then revealed a genetic effect of the PGRS on the impaired hippocampal GMV and hippocampus–mPFC functional connectivity in a large sample of Han Chinese participants that included both the schizophrenia and healthy control groups. Specifically, the PGRS showed a significantly negative correlation with the hippocampal GMV and hippocampus–mPFC functional connectivity, both of which were smaller in the people with schizophrenia than in the controls. We further found that these brain measures (including the hippocampal GMV and hippocampus–mPFC functional connectivity) showed similar changes in unaffected first-degree relatives as in patients. These findings suggested that the genetically influenced hippocampal GMV and hippocampus–mPFC functional connectivity might be important biomarkers for the early diagnosis of schizophrenia.

We found that, regardless of which threshold was chosen, the PGRS in the schizophrenia group was significantly greater than that in the healthy control group, a finding which is consistent with a previous study.Reference Kauppi, Westlye, Tesli, Bettella, Brandt and Mattingsdal19 As a measure of cumulative genetic risk, the finding confirmed that the PGRSs successfully captured the difference in genetic structure between the patients and the controls. Moreover, the threshold in this study was chosen to be 0.1 because it best explained the case–control difference in the genetic structure according to our results from the two-sample t-test (Supplementary Fig. 1). This was also in line with previous studies which reported that the thresholds 0.5, 0.1 and 0.05 could maximally capture the heritability of schizophrenia.1,Reference Tesli, Espeseth, Bettella, Mattingsdal, Aas and Melle20,Reference Ripke, Neale, Corvin, Walters, Farh and Holmans21 In addition, the computation of our PGRSs was based on a meta-analysis of Chinese GWAS samples and Psychiatry Genomics Consortium (PGC2) GWAS samples from Li and colleagues.Reference Li, Chen, Yu, He, Xu and Zhang12 In fact, they have conducted polygenic scoring analyses to predict the case–control status for Chinese participants and found the PGRS derived from the combined data-set (Chinese plus PGC2) explained a larger proportion of the variance than the data-set with only PGC2 or Chinese samples. The limited sample size may reduce the statistical power of genome-wide meta-analysis using only Chinese samples.

By imaging genetic analyses, we found that the higher PGRS was associated with decreased right hippocampal GMV, which was found to be significantly smaller in people with schizophrenia than in healthy controls. In fact, we found that the bilateral hippocampus and other regions included in the default and salience networks showed impaired GMV in the people with schizophrenia, in line with many previous studies.Reference Kalmady, Venkatasubramanian, Shivakumar, Gautham, Subramaniam and Jose22,Reference Manoliu, Riedl, Zherdin, Mühlau, Schwerthöffer and Scherr23 Among these regions, we only found that the right hippocampal GMV had a significantly negative correlation with the PGRS. Therefore, we speculated that the impaired right hippocampal GMV is influenced by genetic factors. Among the limited studies about the associations between schizophrenia-related genetic variants and volume measures, to the best of our knowledge only hippocampal volume was found to be genetically influenced.Reference Smeland, Wang, Frei, Li, Hibar and Franke24 One recent study reported they did not find overlapping evidence of genetic variation between schizophrenia and eight subcortical volume measures, including the hippocampus,Reference Franke, van Hulzen, Arias-Vasquez, Bralten, Hoogman and Klein25 but they still found a negative tendency for a genetic correlation between hippocampal volume and schizophrenia (P = 0.068).Reference Franke, van Hulzen, Arias-Vasquez, Bralten, Hoogman and Klein25 Furthermore, the PGRS in our study was computed using a meta-analysis of previously unused Chinese GWAS samples and PGC2 GWAS samples,Reference Li, Chen, Yu, He, Xu and Zhang12 and the further analysis was also based on a sample of Chinese Han participants. In addition, unlike the research of Franke and colleagues, which considered the whole hippocampal volume,Reference Franke, van Hulzen, Arias-Vasquez, Bralten, Hoogman and Klein25 we performed a voxel-wise regression analysis and found a subfield of the right hippocampus that correlated with the schizophrenia PGRS. This was also in line with a previous study that found that not the total hippocampal volume but a subfield volume was influenced by the schizophrenia PGRS.Reference Papiol, Popovic, Keeser, Hasan, Schneider-Axmann and Degenhardt26 We think that the slightly inconsistent results may be due to differences in ethnicity and methodologies. We further tested the association between impaired hippocampal GMV and cognitive performance and the PANSS scores (Supplementary Table 2). We found that a smaller right hippocampal GMV was significantly associated with a worse SDMT performance (r = 0.067, one-tailed P = 0.019) and a greater PANSS negative score (r = −0.083, one-tailed P = 0.032).

We further observed that the hippocampus–mPFC functional connectivity was reduced in people with schizophrenia and that the lower connectivity was correlated with a higher PGRS. Previous studies have found impaired activity and structure of the mPFC in schizophrenia.Reference Bleich-Cohen, Kupchik, Gruberger, Kotler and Hendler32 A study using an animal model reported that alteration of the functional connectivity between the hippocampus and the mPFC possibly contributed to schizophrenia-like symptoms.Reference Kjaerby, Hovelsø, Dalby and Sotty33 Moreover, deficits in overnight memory consolidation in schizophrenia were found to be associated with hippocampus–mPFC connectivity.Reference Genzel, Dresler, Cornu, Jäger, Konrad and Adamczyk34 Most importantly, the hippocampus and mPFC are two major hubs in the default mode network, a brain system that is reported to modulate internally directed thought, watchfulness and maintenance of self-perception.Reference Pomarol-Clotet, Canales-Rodríguez, Salvador, Sarró, Gomar and Vila35 This network has attracted considerable interest in neuropsychiatric disorders, including schizophrenia.Reference Swanson, Eichele, Pearlson, Kiehl, Yu and Calhoun36 However, until now, no study has revealed that the hippocampus–mPFC functional connectivity was related to the schizophrenia polygenic information. Some genetic variants associated with schizophrenia, such as the glutamate–cysteine ligase catalytic gene,Reference Xin, Mekle, Fournier, Baumann, Ferrari and Alameda37 neuregulen-1/-2Reference Yan, Shamir, Skirzewski, Leiva-Salcedo, Kwon and Karavanova38 and catechol O-methyltransferase messenger RNAReference Yan, Shamir, Skirzewski, Leiva-Salcedo, Kwon and Karavanova38 have been found to influence mPFC phenotypes such as glutathione metabolism dysregulation, extracellular dopamine dysregulation and reduced grey matter concentration. To some extent, these studies support our conclusion that the connectivity from the mPFC is influenced by a schizophrenia genetic mechanism. In addition, we also tested the association between the hippocampus–mPFC functional connectivity and the cognitive abilities and PANSS scores, but found no significant association (Supplementary Table 2).

Our investigation of unaffected individuals at high familial risk of developing schizophrenia found that their identified brain measures were similar to those in the patients with schizophrenia and that their PGRSs fell between the schizophrenia and healthy control groups. These results further supported our finding that the hippocampal GMV and hippocampus–mPFC functional connectivity are important imaging endophenotypes, findings which could help in the early identification and diagnosis of schizophrenia. Currently, schizophrenia is primarily diagnosed using clinical symptoms, and no reliable biological marker has been identified to improve early detection.Reference Neilson, Bois, Clarke, Hall, Johnstone and Owens39 Our study in unaffected first-degree relatives seems to provide complementary evidence for PGRS-related biomarkers and may help to link vulnerability genes to clinical syndromes. In line with our findings, many previous studies found that unaffected first-degree relatives show similar but milder neuroimaging abnormities compared with schizophrenia patientsReference Zhao, Tian, Yan, Yue, Yan and Zhang40 and a higher PGRS than healthy controls.Reference Neilson, Bois, Clarke, Hall, Johnstone and Owens39 But we did not find the significant association between PGRS and the observed neuroimaging abnormalities in unaffected first-degree relatives. The main reason for this may be that these first-degree relatives are of a relatively older age and consist of a smaller sample size. Moreover, some of them are siblings of the patients and others are parents of the patients. Overall, our current findings, accompanied by the supporting evidence from unaffected first-degree relatives, indicate that hippocampal GMV and hippocampus–mPFC functional connectivity could be traits or vulnerability markers for schizophrenia.

Although our study reported plausible schizophrenia polygenic effects on hippocampal GMV and hippocampus–mPFC functional connectivity, several issues still need to be considered and addressed in future studies. First, although the sample size of our study was quite large, the findings still need to be replicated in independent samples. Second, the aim of our study was to identify plausible endophenotypes, so a relatively weak multiple comparison correction method, AlphaSim, was used. However, the results in a cohort of unaffected first-degree relatives further confirmed the clinical significance of our findings. Third, the present study was based on participants of Han Chinese ancestry, so the findings may need to be tested in other ethnic populations. Finally, because the majority of patients included in our study were taking psychotropic drugs during MRI scanning and we had no lifestyle records, the current results did not control for the potential effects of these drugs on people with schizophrenia. To examine whether our observed neuroimaging measures were caused by drug exposure unique to the patient group, we carefully consulted every patient's medication record and found that only 22 of all patients were first-episode drug naive. We then compared the 22 first-episode, drug-naive patients with those of 22 randomly selected healthy controls for 1000 permutation tests and found hippocampal GMV and hippocampus–mPFC functional connectivity in first-episode, drug-naïve patients were significantly smaller than those in healthy controls (P = 0.0186 and 0.0363, respectively). This result indicated that our current finding may be potential biomarkers for schizophrenia. However, because of the small sample size of first-episode, drug-naive patients and no exclusion of the influence of lifestyle factors, future study is needed to further validate our findings.

Supplementary material

Supplementary material is available online at https://doi.org/10.1192/bjp.2019.127.

Funding

This work was supported by the National Key Research and Development Plan (grant number 2016YFC0904300), the Strategic Priority Research Program of Chinese Academy of Science (grant number XDBS01020200) and the Natural Science Foundation of China (grant numbers 31300934 and 81771451).

Acknowledgements

We thank all the patients and controls who took part in this research, as well as the clinical staff who facilitated their involvement.

Footnotes

Declaration of interest None.

References

1Schizophrenia Psychiatric Genome-Wide Association Study (GWAS) Consortium. Genome-wide association study identifies five new schizophrenia loci. Nat Genet 2011; 43: 969–76.10.1038/ng.940CrossRefGoogle Scholar
2Ripke, S, O'Dushlaine, C, Chambert, K, Moran, JL, Kähler, AK, Akterin, S, et al. Genome-wide association analysis identifies 13 new risk loci for schizophrenia. Nat Genet 2013; 45: 1150–9.10.1038/ng.2742CrossRefGoogle Scholar
3Kishi, T, Fukuo, Y, Okochi, T, Kawashima, K, Kitajima, T, Inada, T.Serotonin 6 receptor gene and schizophrenia: case-control study and meta-analysis. Hum Psychopharmacol 2012; 27: 63–9.CrossRefGoogle ScholarPubMed
4Purcell, SM, Wray, NR, Stone, JL, Visscher, PM, O'Donovan, MC, Sullivan, PF, et al. Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature 2009; 460: 748–52.Google ScholarPubMed
5Shepherd, AM, Laurens, KR, Matheson, SL, Carr, VJ, Green, MJ.Systematic meta-review and quality assessment of the structural brain alterations in schizophrenia. Neurosci Biobehav Rev 2012; 36: 1342–56.CrossRefGoogle Scholar
6Benetti, S, Mechelli, A, Picchioni, M, Broome, M, Williams, S, McGuire, P.Functional integration between the posterior hippocampus and prefrontal cortex is impaired in both first episode schizophrenia and the at risk mental state. Brain 2009; 132: 2426–36.CrossRefGoogle ScholarPubMed
7Liu, B, Zhang, X, Cui, Y, Qin, W, Tao, Y, Li, J, et al. Polygenic risk for schizophrenia influences cortical gyrification in 2 independent general populations. Schizophr Bull 2017; 43: 673–80.Google ScholarPubMed
8Walton, E, Geisler, D, Lee, PH, Hass, J, Turner, JA, Liu, J, et al. Prefrontal inefficiency is associated with polygenic risk for schizophrenia. Schizophr Bull 2014; 40: 1263–71.10.1093/schbul/sbt174CrossRefGoogle ScholarPubMed
9Miller, JA, Scult, MA, Conley, ED, Chen, Q, Weinberger, DR, Hariri, AR.Effects of schizophrenia polygenic risk scores on brain activity and performance during working memory subprocesses in healthy young adults. Schizophr Bull 2018; 44: 844–53.10.1093/schbul/sbx140CrossRefGoogle ScholarPubMed
10Walton, E, Turner, J, Gollub, RL, Manoach, DS, Yendiki, A, Ho, BC, et al. Cumulative genetic risk and prefrontal activity in patients with schizophrenia. Schizophr Bull 2013; 39: 703–11.CrossRefGoogle ScholarPubMed
11Wang, C, Liu, B, Zhang, X, Cui, Y, Yu, C, Jiang, T.Multilocus genetic profile in dopaminergic pathway modulates the striatum and working memory. Sci Rep 2018; 8: 613.Google ScholarPubMed
12Li, Z, Chen, J, Yu, H, He, L, Xu, Y, Zhang, D, et al. Genome-wide association analysis identifies 30 new susceptibility loci for schizophrenia. Nat Genet 2017; 49: 1576–83.CrossRefGoogle Scholar
13Purcell, S, Neale, B, Todd-Brown, K, Thomas, L, Ferreira, MAR, Bender, D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007; 81: 559–75.10.1086/519795CrossRefGoogle ScholarPubMed
14Patterson, N, Price, AL, Reich, D.Population structure and eigenanalysis. PLoS Genet 2006; 2: 2074–93.CrossRefGoogle ScholarPubMed
15Price, AL, Patterson, NJ, Plenge, RM, Weinblatt, ME, Shadick, NA, Reich, D.Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet 2006; 38: 904–9.10.1038/ng1847CrossRefGoogle ScholarPubMed
16Thorisson, GA, Smith, AV, Krishnan, L, Stein, LD.The international hapmap project web site. Genome Res 2005; 15: 1592–3.10.1101/gr.4413105CrossRefGoogle ScholarPubMed
17Delaneau, O, Marchini, J, Zagury, JF.A linear complexity phasing method for thousands of genomes. Nat Methods 2012; 9: 179–81.10.1038/nmeth.1785CrossRefGoogle Scholar
18Howie, BN, Donnelly, P, Marchini, J.A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet 2009; 5: e1000529.10.1371/journal.pgen.1000529CrossRefGoogle ScholarPubMed
19Kauppi, K, Westlye, LT, Tesli, M, Bettella, F, Brandt, CL, Mattingsdal, M, et al. Polygenic risk for schizophrenia associated with working memory-related prefrontal brain activation in patients with schizophrenia and healthy controls. Schizophr Bull 2015; 41: 736–43.10.1093/schbul/sbu152CrossRefGoogle ScholarPubMed
20Tesli, M, Espeseth, T, Bettella, F, Mattingsdal, M, Aas, M, Melle, I, et al. Polygenic risk score and the psychosis continuum model. Acta Psychiatr Scand 2014; 130: 311–7.10.1111/acps.12307CrossRefGoogle ScholarPubMed
21Ripke, S, Neale, BM, Corvin, A, Walters, JTR, Farh, KH, Holmans, PA, et al. Biological insights from 108 schizophrenia-associated genetic loci. Nature 2014; 511: 421–7.Google Scholar
22Kalmady, SV, Venkatasubramanian, G, Shivakumar, V, Gautham, S, Subramaniam, A, Jose, DA, et al. Relationship between interleukin-6 gene polymorphism and hippocampal volume in antipsychotic-naïve schizophrenia: evidence for differential susceptibility? PLoS ONE 2014; 9: e96021.CrossRefGoogle ScholarPubMed
23Manoliu, A, Riedl, V, Zherdin, A, Mühlau, M, Schwerthöffer, D, Scherr, M, et al. Aberrant dependence of default mode/central executive network interactions on anterior insular salience network activity in schizophrenia. Schizophr Bull 2014; 40: 428–37.10.1093/schbul/sbt037CrossRefGoogle Scholar
24Smeland, OB, Wang, Y, Frei, O, Li, W, Hibar, DP, Franke, B, et al. Genetic overlap between schizophrenia and volumes of hippocampus, putamen, and intracranial volume indicates shared molecular genetic mechanisms. Schizophr Bull 2018; 44: 854–64.CrossRefGoogle ScholarPubMed
25Franke, B, van Hulzen, KJE, Arias-Vasquez, A, Bralten, J, Hoogman, M, Klein, M, et al. Genetic influences on schizophrenia and subcortical brain volumes: large-scale proof of concept. Nat Neurosci 2016; 19: 420–31.CrossRefGoogle ScholarPubMed
26Papiol, S, Popovic, D, Keeser, D, Hasan, A, Schneider-Axmann, T, Degenhardt, F, et al. Polygenic risk has an impact on the structural plasticity of hippocampal subfields during aerobic exercise combined with cognitive remediation in multi-episode schizophrenia. Transl Psychiatry 2017; 7: e1159.10.1038/tp.2017.131CrossRefGoogle ScholarPubMed
27Riccitelli, GC, Pagani, E, Rodegher, M, Colombo, B, Preziosa, P, Falini, A, et al. Imaging patterns of gray and white matter abnormalities associated with PASAT and SDMT performance in relapsing-remitting multiple sclerosis. Mult Scler J 2017; 25: 204–16.10.1177/1352458517743091CrossRefGoogle ScholarPubMed
28Fukuzako, H, Fukuzako, T, Hashiguchi, T, Hokazono, Y, Takeuchi, K, Hirakawa, K, et al. Reduction in hippocampal formation volume is caused mainly by its shortening in chronic schizophrenia: assessment by MRI. Biol Psychiatry 1996; 39: 938–45.CrossRefGoogle ScholarPubMed
29Guo, X, Li, J, Wang, J, Fan, X, Hu, M, Shen, Y, et al. Hippocampal and orbital inferior frontal gray matter volume abnormalities and cognitive deficit in treatment-naive, first-episode patients with schizophrenia. Schizophr Res 2014; 152: 339–43.10.1016/j.schres.2013.12.015CrossRefGoogle ScholarPubMed
30Haukvik, UK, Westlye, LT, Mørch-Johnsen, L, Jørgensen, KN, Lange, EH, Dale, AM, et al. In vivo hippocampal subfield volumes in schizophrenia and bipolar disorder. Biol Psychiatry 2015; 77: 581–8.10.1016/j.biopsych.2014.06.020CrossRefGoogle ScholarPubMed
31Bachevalier, J, Alvarado, MC, Malkova, L.Memory and socioemotional behavior in monkeys after hippocampal damage incurred in infancy or in adulthood. Biol Psychiatry 1999; 46: 329–39.CrossRefGoogle ScholarPubMed
32Bleich-Cohen, M, Kupchik, M, Gruberger, M, Kotler, M, Hendler, T.Never resting region–mPFC in schizophrenia. Schizophr Res 2012; 140: 155–8.10.1016/j.schres.2012.06.030CrossRefGoogle Scholar
33Kjaerby, C, Hovelsø, N, Dalby, NO, Sotty, F.Phencyclidine administration during neurodevelopment alters network activity in prefrontal cortex and hippocampus in adult rats. J Neurophysiol 2017; 118: 1002–11.10.1152/jn.00081.2017CrossRefGoogle ScholarPubMed
34Genzel, L, Dresler, M, Cornu, M, Jäger, E, Konrad, B, Adamczyk, M, et al. Medial prefrontal-hippocampal connectivity and motor memory consolidation in depression and schizophrenia. Biol Psychiatry 2015; 77: 177–86.10.1016/j.biopsych.2014.06.004CrossRefGoogle Scholar
35Pomarol-Clotet, E, Canales-Rodríguez, EJ, Salvador, R, Sarró, S, Gomar, JJ, Vila, F, et al. Medial prefrontal cortex pathology in schizophrenia as revealed by convergent findings from multimodal imaging. Mol Psychiatry 2010; 15: 823–30.10.1038/mp.2009.146CrossRefGoogle ScholarPubMed
36Swanson, N, Eichele, T, Pearlson, G, Kiehl, K, Yu, Q, Calhoun, VD.Lateral differences in the default mode network in healthy controls and patients with schizophrenia. Hum Brain Mapp 2011; 32: 654–64.CrossRefGoogle ScholarPubMed
37Xin, L, Mekle, R, Fournier, M, Baumann, PS, Ferrari, C, Alameda, L, et al. Genetic polymorphism associated prefrontal glutathione and its coupling with brain glutamate and peripheral redox status in early psychosis. Schizophr Bull 2016; 42: 1185–96.CrossRefGoogle ScholarPubMed
38Yan, L, Shamir, A, Skirzewski, M, Leiva-Salcedo, E, Kwon, OB, Karavanova, I, et al. Neuregulin-2 ablation results in dopamine dysregulation and severe behavioral phenotypes relevant to psychiatric disorders. Mol Psychiatry 2018; 23: 1233–43.CrossRefGoogle ScholarPubMed
39Neilson, E, Bois, C, Clarke, TK, Hall, L, Johnstone, EC, Owens, DGC, et al. Polygenic risk for schizophrenia, transition and cortical gyrification: a high-risk study. Psychol Med 2018; 48: 1532–39.CrossRefGoogle ScholarPubMed
40Zhao, X, Tian, L, Yan, J, Yue, W, Yan, H, Zhang, D.Abnormal rich-club organization associated with compromised cognitive function in patients with schizophrenia and their unaffected parents. Neurosci Bull 2017; 33: 445–54.CrossRefGoogle ScholarPubMed
Figure 0

Table 1 Demographic and clinical characteristics

Figure 1

Fig. 1 Association between the schizophrenia polygenic risk score (PGRS) and impaired grey matter volume (GMV). (a) A two-sample t-test of the GMV between schizophrenia and healthy control groups with AlphaSim correction (single voxel P < 0.005, the corrected threshold P < 0.05 and cluster size threshold >161 voxels), the negative t-value representing the GMV of this region was significantly decreased in people with schizophrenia. (b) Multiple regression analysis testing the association between PGRS and disrupted GMV with AlphaSim correction (single voxel P < 0.005, the corrected threshold P < 0.05 and cluster size threshold >132 voxels).

Figure 2

Fig. 2 Association between the schizophrenia polygenic risk score (PGRS) and disrupted functional connectivity. (a) Mean graph of the functional connectivity between the hippocampus and the whole brain in people with schizophrenia. (b) Mean graph of the functional connectivity between the hippocampus and the whole brain in the healthy control group. (c) Two-sample t-test results for the functional connectivity in the schizophrenia and healthy control groups with AlphaSim correction (single voxel P < 0.005, corrected threshold P < 0.05 and cluster size threshold >21 voxels), the negative t-value representing the functional connectivity was decreased in the people with schizophrenia and increased in the healthy controls. (d) Multiple regression analysis testing the association between the PGRS and the disrupted functional connectivity with AlphaSim correction (single voxel P < 0.01, the corrected threshold P < 0.05 and cluster size threshold >56 voxels). L, left; R, right.

Figure 3

Fig. 3 Polygenic risk score (PGRS), right hippocampal grey matter volume (GMV) and hippocampus–medial prefrontal cortex (mPFC) functional connectivity in the healthy control (HC) group, unaffected first-degree relatives (high-risk individuals, HR) and people with schizophrenia (SZ). (a) Distribution of the PGRS. (b) Distribution of the right hippocampal GMV. (c) Distribution of the hippocampus–mPFC functional connectivity (FC). *P < 0.05, **P < 0.01, ***P < 0.001.

Supplementary material: File

Liu et al. supplementary material

Figure S1 and Tables S1-S2

Download Liu et al. supplementary material(File)
File 308.2 KB
Submit a response

eLetters

No eLetters have been published for this article.