Skip to main content

Transcriptional level of inflammation markers associates with short-term brain structural changes in first-episode schizophrenia

Abstract

Background

Inflammation has been implicated in the pathology of schizophrenia and may cause neuronal cell death and dendrite loss. Neuroimaging studies have highlighted longitudinal brain structural changes in patients with schizophrenia, yet it is unclear whether this is related to inflammation. We aim to address this question, by relating brain structural changes with the transcriptional profile of inflammation markers in the early stage of schizophrenia.

Methods

Thirty-eight patients with first-episode schizophrenia and 51 healthy controls were included. High-resolution T1-weighted magnetic resonance imaging (MRI) and clinical assessments were performed at baseline and 2 ~ 6 months follow-up for all subjects. Changes in the brain structure were analyzed using surface-based morphological analysis and correlated with the expression of immune cells-related gene sets of interest reported by previous reviews. Transcriptional data were retrieved from the Allen Human Brain Atlas. Furthermore, we examined the brain structural changes and peripheral inflammation markers in association with behavioral symptoms and cognitive functioning in patients.

Results

Patients exhibited accelerated cortical thickness decrease in the left frontal cortices, less decrease or an increase in the superior parietal lobule and right lateral occipital lobe, and increased volume in the bilateral pallidum, compared with controls. Changes in cortical thickness correlated with the transcriptional level of monocyte across cortical regions in patients (r = 0.54, p < 0.01), but not in controls (r =  − 0.05, p = 0.76). In addition, cortical thickness change in the left superior parietal lobule positively correlated with changes in digital span-backward test scores in patients.

Conclusions

Patients with schizophrenia exhibit regional-specific cortical thickness changes in the prefrontal and parietooccipital cortices, which is related to their cognitive impairment. Inflammation may be an important factor contributing to cortical thinning in first-episode schizophrenia. Our findings suggest that the immunity-brain-behavior association may play a crucial role in the pathogenesis of schizophrenia.

Peer Review reports

Background

Schizophrenia is a severe mental disorder ranked among the top 20 causes of disability worldwide [1]. The prolonged illness course and impaired cognitive functioning are the main reasons contributing to the disease burden [2, 3]. Yet, the pathological mechanism of schizophrenia remains to be unclear. Imaging genetics study provides an opportunity to quantify disease-related neuroanatomical deviations and to elucidate the potential biological mechanism underlying these changes. In recent years, imaging genetics analysis has been increasingly applied to schizophrenia to reveal its potential pathological pathway [4, 5].

The suggestion that immune dysfunction may contribute to the pathophysiology of schizophrenia has a long history [6]. The vulnerability-stress-inflammation model [7] proposes that genetic risks and early life exposures, such as stress or infection, may promote a chronic pro-inflammatory state and subsequently lead to neurotransmission disturbances and brain structural changes that are relevant for schizophrenia [7, 8]. Supporting this hypothesis, therapeutic studies indicated a beneficial effect of using anti-inflammatory medications as an add-on treatment in the early stage of schizophrenia [7]. Studies of inflammation indicated that genetically predicted IL-6 was associated with brain volume in the middle temporal gyrus and may potentially be involved in the pathology of schizophrenia [9]. Likewise, studies also found overlapped genetic loci of schizophrenia and cortical morphometry to be enriched in immunologic gene sets [4]. However, few longitudinal studies have examined whether treatment-associated brain structural changes are related to neuroinflammation or gene expression of inflammatory markers.

Studies have indicated that progressive brain structural changes have been shown in the general population and patients with schizophrenia [10, 11]. In schizophrenia, previous studies found subcortical enlargement and cortical thinning in widespread brain regions [11, 12]. The frontal and temporal cortices were the most affected regions with excessive cortical thinning in patients than controls [11, 13]. Several large-scale, longitudinal studies have been conducted on schizophrenia [14], and a few studies focused on the early stage of illness when the majority of patients were medication-naïve and experienced substantial symptom alleviation [15, 16]. However, most studies did not recruit controls for follow-up, and it is difficult to distinguish the effects of disease progression and natural aging.

Considering that surface area and cortical thickness reflect different aspects of neural development processes [17] and may be differently affected in schizophrenia, we separately analyzed cortical thickness, surface area, and subcortical volume using surface-based morphometry (SBM). The aim of this study is threefold: (1) to determine the early-stage brain structural changes after a short period of treatment in patients with schizophrenia, (2) to investigate factors (inflammation and antipsychotic treatment) that may contribute to brain structural changes by correlating structural changes with transcription level of gene sets of interest, and (3) to evaluate the association of brain structural changes with changes in symptoms and cognitive performances in patients.

Methods

Participants

In this study, 38 first-episode patients with schizophrenia were recruited during their hospitalization in the Department of Psychiatry, Xijing Hospital, between May 2015 and October 2017 [18] (Table 1). The diagnosis was based on the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition [19]. Patients were included if this was the first time they were in the hospital or sought outpatient clinical help for their mental disorders. Additional inclusion criteria for patients require no more than 2 weeks of cumulative exposure to antipsychotics [20]. The exclusion criteria were listed in the Additional file 1: Supplementary Methods [20, 21]. All patients received second-generation antipsychotic medications determined by the clinicians during the study period, and 36 out of 38 patients received combined rTMS treatment during hospitalization (Additional file 1: Supplementary Methods) [22]. After a period of treatment (mean treatment duration: 148.29 ± 56.20 days), patients received a follow-up MRI scanning and assessment of clinical symptoms and cognitive functions.

Table 1 Demographic and clinical profile of 38 first-episode patients with schizophrenia and 51 healthy controls

Fifty-one healthy controls were recruited from the local community through advertisement and social media, with gender- and age-matched to the patient group (Table 1). The healthy control group received MRI scanning and cognitive assessment at recruitment and follow-up periods (131.51 ± 4.46 days). Written informed consent was signed by all subjects or parents of subjects under 18 years of age, after understanding this study in detail.

Clinical assessment, cognitive tests, and laboratory tests

Symptom severity of the enrolled patients was evaluated using the Positive and Negative Syndrome Scale (PANSS) at baseline and follow-up time points [23]. PANSS subscale scores were calculated including positive symptoms, negative symptoms, and general psychopathology. We used ΔPANSS to indicate the degree of symptom relief in patients after treatment (T1) and before treatment (T0): ∆PANSS = (PANSS_T1 − PANSS_T0)/(PANSS_T0 − 30). The Wechsler Adult Intelligence Scale (WAIS), revised in China [24], was applied to patients and controls to evaluate their cognitive functioning. We performed the WAIS digital span (forward and backward) test and digital symbol coding test on all participants. In addition to the MRI scans, we also collected routine laboratory records of blood samples from the clinical database. Specifically, 2 ml of intravenous blood was collected from the upper forearm in an EDTA vacuum tube under fasting condition after admission, and the samples were measured using a Sysmex XS Automatic Hematology Analyzer. The data on leukocyte count and percentage were extracted from baseline laboratory tests to analyze the inflammatory indicators.

MRI acquisition and preprocessing

High-resolution T1-weighted MRI was acquired using a GE Discovery MR750 3.0 T scanner located in the Department of Radiology, and all subjects underwent T2WI scans to rule out organic diseases. Participants showing gross artifacts and/or excessive head motion were excluded, following protocols established previously [21]. The MRI acquisition was accomplished at baseline and follow-up, at the same time as the clinical assessment. The scanning parameters were repetition time 8.2 ms, echo time 3.2 ms, flip angle 12°, field of view 256 × 256 mm, matrix 256 × 256, slice thickness 1 mm, and sagittal slices 196. The image preprocessing followed the analysis process of FreeSurfer (version 6.0.0, https://surfer.nmr.mgh.harvard.edu) and was performed on the Linux operating system. The brief preprocessing steps included correcting the uniformity of the bias field; removing the skull; segmenting the brain tissue into gray matter, white matter, and cerebrospinal fluid; cortical surface reconstruction; registration of individual space to the standard space by nonlinear transformation; the registered image was segmented into 68 cortices based on Desikan-Killiany atlas and 14 subcortical regions. For each subject, 150 regional structural measurements were obtained, including 68 regional cortical thicknesses, 68 surface areas, and 14 subcortical volumes.

Linear mixed-effect model on brain structural changes

The statistical analysis was performed on MATLAB (R2020b) platform. A linear mixed-effect model was applied to simultaneously test the group differences (group), time point (time), and their interaction effect on global and regional morphological measurements (Y). Age and sex were included as covariates. When comparing regional measurements, we additionally included global measurements to adjust for global differences between individuals. Each subject was given a subject ID (SubID) which was included as the random factor. Multiple testing effect was corrected using false discovery rate (FDR) correction with q < 0.05 for regional measurement comparisons.

Brain transcriptional data

Brain transcriptional data were obtained from the Allen Human Brain Atlas (AHBA) [25] (http://human.brain-map.org/), including 20,734 genes across the brain from six healthy human donors (five males, age range of 24 ~ 57). For each tissue sample, gene expression data represented by 58,692 probes were quantified, normalized, and averaged across probes following a previous study [26]. We used gene expression data from the left hemisphere as they are available from all donors. Tissue samples were mapped to the Desikan-Killiany atlas [27] based on Euclidean distance to the nearest voxel in the MNI 152 template [26]. For each gene, the expression level was averaged across the six donors and normalized to z score across all cortical regions, resulting in a normalized gene expression matrix of size 34 × 20,734 (region × genes) for further analysis. Detailed processing steps are described in Additional file 1: Supplementary Methods [25,26,27,28].

To identify inflammation-related genes, we selected genes from a previous systematic review on leukocytes and subtypes (monocyte, lymphocyte, neutrophil, eosinophil, and basophil) [29]. To reduce potential bias in gene set selection, we further verified our results using inflammation-related genes reported in another review [30]. Genes related to antipsychotic treatment response were selected from a genome-wide association study (GWAS) on antipsychotic treatment in schizophrenia [31]. Included genes are listed in Additional file 1: Table S1.

Association between brain structural changes and peripheral inflammation markers

We next examined the associations between brain structural changes and baseline peripheral inflammation markers in patients. Inflammation markers were quantified including leukocyte count and percentage of leukocyte subtypes (monocyte, lymphocyte, neutrophil, eosinophil, and basophil) (Table 2). As some inflammation markers do not follow a normal distribution (Additional file 1: Table S2), Spearman’s rank correlation was performed between inflammation markers and structural changes in the brain regions that showed a significant group main effect or group by time interaction effect. We further examined the associations between inflammation markers and subcortical volume changes, as some previous studies reported an association between subcortical volume and peripheral inflammatory markers [32], but not consistently [33].

Table 2 Leukocyte count and leukocyte subtypes (percentage) in patients

Association between brain structural changes and transcriptional profiles

We next investigated whether brain structural changes were associated with gene expression patterns related to inflammation and antipsychotic response (Additional file 1: Table S1). For each set of genes of interest, normalized gene expression was averaged across samples. The normalized gene expression maps were cross-correlated with group-averaged cortical thickness and surface area changes in patients. The same correlation analysis was also performed in controls to examine whether the effects were unique to patients.

As recent studies pointed out that spatial autocorrelation of brain maps results in false-positive effects [34], we verified our correlation results between brain and transcriptional maps using generative null models proposed by Burt et al. [35]. Specifically, 5000 surrogate maps preserving spatial autocorrelation of brain structural changes in patients were generated using the BrainSMASH toolkit [35]. The null distribution of correlations between gene expression and surrogate brain maps was calculated and compared with our empirical results. The two-sided p value was determined as the proportion of random permutations that exceeded the empirical correlation coefficient between gene expression and brain maps.

Brain structural changes in relation to clinical profile and cognitive functioning

To determine the behavioral relevance of brain structural changes in patients, we performed Pearson correlation between brain structural changes in regions showing a significant group or group by time interaction effect and clinical factors: PANSS subscale score changes and olanzapine equivalent antipsychotic dosage. We also examined the correlations between brain structural changes and WAIS digital symbol coding and digital span (forward and backward) score changes in both patients and controls. Multiple testing (11 regions × 7 behavioral items) was corrected using the FDR at q < 0.05.

Peripheral inflammation markers in relation to clinical profile and cognitive functioning

Spearman’s rank correlation was performed between inflammation markers (leukocyte and leukocyte subtypes) and patients’ behavioral scores (PANSS and WAIS subscales) at baseline and follow-up.

Results

Demographic characteristics and clinical outcomes

Demographic and clinical information of all participants were shown in Table 1. Controls were matched with patients on age (p = 0.09) and sex (p = 0.37), and controls had a higher education level (p = 0.01) and shorter follow-up interval (p = 0.04) than patients. The patient’s symptom severity at baseline and follow-up period was illustrated in Table 3. The majority of patients (33 out of 38) had an improvement in their symptoms, with PANSS total score reduced by more than 30% (Additional file 1: Fig. S1). On cognitive assessments, patients showed significantly lower scores on WAIS digital symbol coding and digital span (forward and backward) tests compared to controls (all p < 0.001). Patients did not show significant improvements in WAIS scores at follow-up. Information on patients who were medicated before admission is listed in Additional file 1: Table S3.

Table 3 PANSS score and cognitive functions before and after treatments

Group comparisons of brain structural changes

Global effects

In both patients and controls, brain structural measurements decreased at the follow-up time point compared to baseline (Table 4). Global volume, cortical thickness, and surface area on average decreased by 1.08%, 0.87%, and 0.65% in patients and 0.97%, 0.76%, and 0.24% in controls, respectively (Fig. 1). Linear mixed-effect model showed a significant main effect of time on global volume (t =  − 2.34, p = 0.02) and cortical thickness (t =  − 2.36, p = 0.02), but not for surface area (t =  − 0.52, p = 0.61). No significant group main effect or group by time interaction effect was found in global measurements.

Table 4 Global structural measurements and change ratio in patients and controls
Fig. 1
figure 1

Percentage of cortical thickness (upper panel) and surface area (lower panel) changes in patients with schizophrenia (left) and healthy controls (HC). Blue color indicates a decreased value at follow-up compared to baseline measurements, whereas red color indicates an increased value

Regional effects

For regional structural changes, we analyzed the group by time interaction effects and group main effects on 150 regional measurements (14 subcortical volumes, 68 cortical thicknesses, and 68 surface areas). For group by time interaction effects, several regional cortical thicknesses showed a significant result, including the left superior and middle frontal gyri, left caudal anterior cingulate gyrus, bilateral superior parietal lobules, right inferior parietal lobule, right lateral occipital lobe, and bilateral pallidum. For significant group main effects, the right superior frontal gyrus showed significantly lower cortical thickness whereas the right entorhinal cortex showed higher cortical thickness in patients than in controls (Fig. 2 and Table 5). No significant group or group by time interaction effects were found for surface area comparison. We conducted a sub-group analysis of patients who were medicated at baseline, which showed similar results (Additional file 1: Table S4).

Fig. 2
figure 2

Significant group by time interaction effect (A) and group effect (B) on brain structural measurements. T-statistics from linear mix-effects model was plotted on the left, with blue color indicating lower values or accelerated decrease in patients, while red indicates the opposite. For significant regions, cortical thickness and volume in patients (SZ) and controls (HC) at baseline (BL) and follow-up (FU) were plotted on the right.  Abbreviations: L: left, R: right, SFG: superior frontal gyrus, cauACC: caudal anterior cingulate cortex, rosMFG, rostral middle frontal gyrus, SPL: superior parietal lobule, IPL: inferior parietal lobule, latOL: lateral occipital lobe, EC: entorhinal cortex

Table 5 Significant group by time interaction and group main effects on regional measurements (FDR significant results in bold font)

Association between brain structural changes and peripheral inflammation markers

We first tested whether the brain regions that showed a significant group main effect or group by time interaction effect would correlate with patients’ peripheral inflammation levels. Next, we performed an exploratory analysis in the subcortical regions. We did not observe a significant association between structural changes and inflammation markers.

Association between brain structural changes and transcriptional profiles

Next, we investigated whether brain structural changes correlated with the gene expression profile of inflammation and antipsychotic treatment (Additional file 1: Fig. S2). Our analysis revealed a positive correlation between cortical thickness changes and gene expression of monocyte [29] in patients (r = 0.54, p = 8.8 × 10−4, FDR corrected) but not in healthy controls (r =  − 0.05, p = 0.76) (Fig. 3 and Table 6). This result was replicated using inflammation genes from another study, which showed a positive correlation between cortical thickness changes and gene expression of monocyte in patients (r = 0.51, p = 0.002) [30]. We further tested whether this micro–macro association was driven by spatial autocorrelation of neuroimaging measurements. Our empirical results exceeded the null distribution of correlations between 5000 surrogate maps and monocyte gene expression (two-tail p = 0.001 for [29] and p = 0.006 for [30]). The correlation between cortical thickness changes and antipsychotic genes in patients is r = 0.37, p = 0.03, but did not survive FDR correction [31]. We also conducted a sub-group analysis of patients who were medicated at baseline, which showed a significant correlation with monocyte gene expression level (Additional file 1: Table S5).

Fig. 3
figure 3

Association between brain structural changes and transcriptional profiles of inflammation and antipsychotic response. A Transcriptional data were retrieved from the Allen Human Brain Atlas (AHBA) [25] (http://human.brain-map.org/) and calculated for gene sets of interests [29,30,31] (Additional file 1: Table S1). Brain structural changes (cortical thickness and surface area) were calculated for patients with schizophrenia (SZ) and healthy controls (HC) as illustrated in Fig. 1. B A significant correlation was found between monocyte-related gene expression and cortical thickness changes in patients with schizophrenia (r = 0.54, p < 0.01) but not in healthy controls (r =  − 0.05, p = 0.76). C The normalized monocyte gene expression profile from a systematic review on leukocyte and subtypes [29]. D Association between monocyte gene expression and cortical thickness changes in patients is significantly higher than the null model generated by correlations between 5000 surrogate maps and monocyte gene expression profile (pspin = 0.001)

Table 6 Correlation between structural changes and gene expression of antipsychotic treatment response and inflammation in patients (SZ) and controls (HC)

Brain structural changes in relation to clinical factors and cognitive functioning

For WAIS tests, we found that the left superior parietal lobule thickness change correlated with WAIS digital span (backward) change in patients (r = 0.57, p = 3.4 × 10−4) but not in controls (r = 0.09, p = 0.54) (Fig. 4). For PANSS score changes and medication dosage, we did not find any significant correlations with brain structural changes.

Fig. 4
figure 4

Correlation analysis between brain structure changes and cognitive score changes in patients (A) and controls (B). A significant association was found between the left superior parietal lobule (SPL) cortical thickness change and WAIS digital span (backward) change in patients (r = 0.57, p = 3.4 × 10−4) but not in controls (r = 0.09, p = 0.54)

Correlation between patients’ peripheral inflammation markers and clinical profile

We found an uncorrected association between monocyte percentage and PANSS positive score (Spearman’s r = 0.44, p = 0.006), as well as between monocyte percentage and digital span (backward) score (Spearman’s r = 0.38, p = 0.02) at baseline (Table 7). Yet, these correlations did not survive multiple corrections.

Table 7 Spearman’s correlation between patients’ inflammation markers and clinical scores

Discussion

This longitudinal study examined the patterns of brain structural changes in patients with first-episode schizophrenia and investigated the contributing factors (inflammation and antipsychotic treatment) to brain structural changes. We found that patients had accelerated pallidum enlargement and frontal cortical thinning but preserved parietal and occipital cortical thickness as compared to controls. Among these structural changes, cortical thickness change in the left superior parietal lobule positively correlated with cognitive performance changes in patients. Notably, we observed a positive correlation between the gene expression level of monocyte and cortical thinning in patients. Taken together, our results identified early-stage brain structural changes and their correlation with the transcriptional profile of inflammatory markers, providing preliminary evidence of the underlying biological process that may contribute to longitudinal brain structural changes in schizophrenia.

Immunological abnormalities have long been reported to be involved in schizophrenia. In previous schizophrenia studies, higher peripheral monocyte count could be relevant to the pathophysiology [36] and has been considered a possible marker of microglia activation [37]. Inflammation in the brain involved both monocyte and microglia. Microglial cells and chemokines might recruit monocytes circulating in the peripheral blood into the brain [38], where neuroinflammation was exacerbated by these monocytes. The investigation of monocyte transcriptome demonstrated a shift in monocyte phenotype in different stages of schizophrenia [39]. In the current study, the transcriptional profile of monocyte was directly related to the longitudinal change of cortical thickness in patients, suggesting an association of monocyte gene expression with the effect of treatment and this disease itself on the brain.

However, we did not observe a significant association between structural changes and inflammation markers. Previous studies revealed elevated peripheral inflammatory levels in patients with schizophrenia and their correlation with increased basal ganglia volume [40]. The effect of peripheral inflammatory markers on brain structural changes was further supported by genetic studies implicating multiple biological processes such as neural development and synaptic transmission [9]. Leukocyte count, even within the normal range, was correlated with polygenic scores of schizophrenia [41]. Since the antipsychotic medication was related to the decrease of cytokines [42], and basal ganglia volume change was quite common after medication [43], one possible mechanism of antipsychotics in the treatment of psychosis could be mitigating neuroinflammation.

Although antipsychotics can alleviate symptoms of schizophrenia, some studies raised concern that they may contribute to additional cortical thinning [44]. Subsequent studies argued that cortical thinning after medication was not linked to clinical or cognitive deterioration [45] and might even improve patients’ prefrontal functional activity and cognitive control ability [13]. In a macaque monkey study, antipsychotic-induced volume reduction was related to reduced glial cell number, whereas the number of neurons remained unchanged [46]. In our study, we found the transcriptional profile of monocyte positively correlated with cortical thickness changes after treatment, indicating the important role of inflammation, especially monocyte, in cortical thickness changes in the early stage of schizophrenia.

Several other factors may contribute to structural changes in patients with schizophrenia, such as medication and neuromodulation treatment. Studies found that medication duration and dosage correlated with cortical thinning in patients [11, 45], with effects varied across different types of antipsychotics [14]. As a non-invasive adjuvant therapy, repetitive transcranial magnetic stimulation (rTMS) may prevent volume reduction or cortical thinning [47] and increase brain metabolism at the stimulation site [48]. The temporoparietal junction is a commonly applied target site.

We also analyzed the factors that might contribute to these structural changes, including medication, times of rTMS, and peripheral inflammatory markers. Although given the fact that treatment factors were able to cause structural changes, we did not find a significant correlation between brain changes and the medication dose, duration, or times of rTMS as former studies reported [11, 45]. Three reasons might contribute to this. First, it might be due to our lack of comparison to unmediated patients, since the disease course would also lead to cortical thinning even in prodromal high-risk individuals [49]. Thus, when calculating the effect of treatment, the effect was not linearly related to structural changes due to the interference of the disease course. Second, it could also be on account of the combined treatment of multiple antipsychotics for each patient, because some studies argued that different types of antipsychotics had varied influences on structural changes as well [14]. Lastly, different patients had different responses to antipsychotics, and our limited sample size might not suffice to prove the correlation with treatment factors.

Among the brain regions showing significant group or group by time interaction effects, we found a positive correlation between superior parietal lobule thickness change and WAIS digital span score change in patients, suggesting the cognitive relevance of this region. The superior parietal lobule belongs to the association cortex and is involved in high-order processes like cognition, information integration, and self-awareness [50]. Structural and functional abnormalities of the frontal-parietal network were frequently reported and have been related to cognition impairments in schizophrenia. Previous studies have reported that low-frequency rTMS could regulate functional connectivity within the frontal-parietal network and improve cognitive functioning [47]. Our study supported the role of the superior parietal lobule in working memory and suggested this potential region as a stimulation target site for cognitive remediation.

In the study, we found a regional-specific effect of accelerated frontal cortical thinning and pallidum enlargement as well as preserved parietal and occipital cortices that may be due to medication and rTMS, respectively. Our results were consistent with former cross-sectional studies [11,12,13, 43] and provided a longitudinal view in patients with first-episode schizophrenia who were unmedicated and responsive to treatments. Remarkably, consistent with previous meta-analysis [11] and longitudinal study [45], cortical thickness represented more prominent changes in regional analysis rather than surface area, which might reflect different mechanisms underlying these two measurements. It is assumed that the developmental trajectories of surface area were predominantly influenced by genetic factors [11, 17]. Meanwhile, cortical thickness can be affected by additional environmental or neurodegenerative factors (disease, inflammation, treatment, aging, etc.) even in adulthood [17]. Our study supported that cortical thickness might be more sensitive to treatment and inflammation effects.

Several limitations need to be considered when interpreting our results. First, transcriptomic data were not collected in this sample but retrieved from the external AHBA dataset, leading to false-positive effects in micro–macro association due to spatial autocorrelation. Using generative null models, the correlation results between brain and transcriptional maps were verified. Second, our findings of brain structural changes reflect mixed effects of treatment and disease progression on the brain, which cannot be easily separated apart. However, we tried to control for the time effect by including follow-up MRI measurements of healthy controls in linear mixed-effect models. Third, the sample size of this study was moderate. Future studies with a larger sample size or using a meta-analytic approach could further confirm our findings.

Conclusions

In summary, this study found brain structural changes in patients with first-episode schizophrenia after short-term treatment. The cortical thickness change in patients was related to the gene expression level of inflammation, and cortical thickness in the superior parietal lobule correlated with ameliorated cognitive impairments in patients. These results provided empirical evidence of potential biological processes underlying brain structural changes and suggested possible treatment targets to improve cognitive function in patients. More importantly, identified immunity-brain-behavior associations would contribute to our understanding of the pathogenesis of schizophrenia.

Availability of data and materials

Data used in the present study can be accessed upon request from the corresponding authors. The code of the project is available at https://github.com/xchang007/SCZ_treatment/tree/main.

Abbreviations

AHBA:

Allen Human Brain Atlas

FDR:

False discovery rate

GMV:

Gray matter volume

GWAS:

Genome-wide association study

rTMS:

Repetitive transcranial magnetic stimulation

SBM:

Surface-based morphometry

References

  1. James SL, Abate D, Abate KH, Abay SM, Abbafati C, Abbasi N, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018;392(10159):1789–858.

  2. Wu X-S, Yan T-C, Wang X-Y, Cao Y, Liu X-F, Fu Y-F, et al. Magnetic resonance imaging-guided and navigated individualized repetitive transcranial magnetic stimulation for cognitive impairment in schizophrenia. Neurosci Bull. 2021;37(9):1365–9.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Zhao Q, Li J, Xiao Y, Cao H, Wang X, Zhang W, et al. Distinct neuroanatomic subtypes in antipsychotic-treated patients with schizophrenia classified by the predefined classification in a never-treated sample. Psychoradiology. 2021;1(4):212–24.

    Article  Google Scholar 

  4. Cheng W, Frei O, van der Meer D, Wang Y, O’Connell KS, Chu Y, et al. Genetic association between schizophrenia and cortical brain surface area and thickness. JAMA Psychiat. 2021;78(9):1020–30.

    Article  Google Scholar 

  5. Patel Y, Parker N, Shin J, Howard D, French L, Thomopoulos SI, et al. Virtual histology of cortical thickness and shared neurobiology in 6 psychiatric disorders. JAMA Psychiat. 2021;78(1):47–63.

    Article  Google Scholar 

  6. Buckley PF. Neuroinflammation and schizophrenia. Curr Psychiatry Rep. 2019;21(8):72.

    Article  PubMed  Google Scholar 

  7. Müller N. Inflammation in schizophrenia: pathogenetic aspects and therapeutic considerations. Schizophr Bull. 2018;44(5):973–82.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Khandaker GM, Cousins L, Deakin J, Lennox BR, Yolken R, Jones PB. Inflammation and immunity in schizophrenia: implications for pathophysiology and treatment. Lancet Psychiatry. 2015;2(3):258–70.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Williams JA, Burgess S, Suckling J, Lalousis PA, Batool F, Griffiths SL, et al. Inflammation and brain structure in schizophrenia and other neuropsychiatric disorders: a mendelian randomization study. JAMA Psychiat. 2022;79(5):498–507.

    Article  Google Scholar 

  10. Bethlehem RAI, Seidlitz J, White SR, Vogel JW, Anderson KM, Adamson C, et al. Brain charts for the human lifespan. Nature. 2022;604(7906):525–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. van Erp TGM, Walton E, Hibar DP, Schmaal L, Jiang W, Glahn DC, et al. Cortical brain abnormalities in 4474 individuals with schizophrenia and 5098 control subjects via the Enhancing Neuro Imaging Genetics Through Meta Analysis (ENIGMA) Consortium. Biol Psychiatry. 2018;84(9):644–54.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Di Sero A, Jørgensen KN, Nerland S, Melle I, Andreassen OA, Jovicich J, et al. Antipsychotic treatment and basal ganglia volumes: exploring the role of receptor occupancy, dosage and remission status. Schizophr Res. 2019;208:114–23.

    Article  PubMed  Google Scholar 

  13. Lesh TA, Tanase C, Geib BR, Niendam TA, Yoon JH, Minzenberg MJ, et al. A multimodal analysis of antipsychotic effects on brain structure and function in first-episode schizophrenia. JAMA Psychiat. 2015;72(3):226–34.

    Article  Google Scholar 

  14. van Haren NE, Schnack HG, Cahn W, van den Heuvel MP, Lepage C, Collins L, et al. Changes in cortical thickness during the course of illness in schizophrenia. Arch Gen Psychiatry. 2011;68(9):871–80.

    Article  PubMed  Google Scholar 

  15. Andreasen NC, Nopoulos P, Magnotta V, Pierson R, Ziebell S, Ho B-C. Progressive brain change in schizophrenia: a prospective longitudinal study of first-episode schizophrenia. Biol Psychiatry. 2011;70(7):672–9.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Yang C, Tang J, Liu N, Yao L, Xu M, Sun H, et al. The effects of antipsychotic treatment on the brain of patients with first-episode schizophrenia: a selective review of longitudinal MRI studies. Front Psychiatry. 2021;12:593703.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Grasby KL, Jahanshad N, Painter JN, Colodro-Conde L, Bralten J, Hibar DP, et al. The genetic architecture of the human cerebral cortex. Science. 2020;367(6484):eaay6690.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Cui L-B, Yin H. The Xi’an Schizophrenia Imaging Lab (SIL) data and ten years of MRI study on schizophrenia. Psychoradiology. 2022;2(2):54–5.

    Article  Google Scholar 

  19. Stein DJ, Phillips KA, Bolton D, Fulford KW, Sadler JZ, Kendler KS. What is a mental/psychiatric disorder? From DSM-IV to DSM-V. Psychol Med. 2010;40(11):1759–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Cui L-B, Liu L, Wang H-N, Wang L-X, Guo F, Xi Y-B, et al. Disease definition for schizophrenia by functional connectivity using radiomics strategy. Schizophr Bull. 2018;44(5):1053–9.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Cui L-B, Wei Y, Xi Y-B, Griffa A, De Lange SC, Kahn RS, et al. Connectome-based patterns of first-episode medication-naïve patients with schizophrenia. Schizophr Bull. 2019;45(6):1291–9.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Hoffman RE, Hawkins KA, Gueorguieva R, Boutros NN, Rachid F, Carroll K, et al. Transcranial magnetic stimulation of left temporoparietal cortex and medication-resistant auditory hallucinations. Arch Gen Psychiatry. 2003;60(1):49–56.

    Article  PubMed  Google Scholar 

  23. Kay SR, Fiszbein A, Opler LA. The Positive and Negative Syndrome Scale (PANSS) for schizophrenia. Schizophr Bull. 1987;13(2):261–76.

    Article  CAS  PubMed  Google Scholar 

  24. Gong YX. Wechsler Adult Intelligence Scale-Revised in China. Changsha: Hunan Medical College; 1982.

    Google Scholar 

  25. Hawrylycz MJ, Lein ES, Guillozet-Bongaarts AL, Shen EH, Ng L, Miller JA, et al. An anatomically comprehensive atlas of the adult human brain transcriptome. Nature. 2012;489(7416):391–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Wei Y, de Lange SC, Scholtens LH, Watanabe K, Ardesch DJ, Jansen PR, et al. Genetic mapping and evolutionary analysis of human-expanded cognitive networks. Nat Commun. 2019;10(1):4839.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Desikan RS, Ségonne F, Fischl B, Quinn BT, Dickerson BC, Blacker D, et al. An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage. 2006;31(3):968–80.

    Article  PubMed  Google Scholar 

  28. Arnatkeviciute A, Fulcher BD, Fornito A. A practical guide to linking brain-wide gene expression and neuroimaging data. Neuroimage. 2019;189:353–67.

    Article  PubMed  Google Scholar 

  29. Nalls MA, Couper DJ, Tanaka T, van Rooij FJA, Chen M-H, Smith AV, et al. Multiple loci are associated with white blood cell phenotypes. PLoS Genet. 2011;7(6):e1002113.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Keller MF, Reiner AP, Okada Y, van Rooij FJA, Johnson AD, Chen M-H, et al. Trans-ethnic meta-analysis of white blood cell phenotypes. Hum Mol Genet. 2014;23(25):6944–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Yu H, Yan H, Wang L, Li J, Tan L, Deng W, et al. Five novel loci associated with antipsychotic treatment response in patients with schizophrenia: a genome-wide association study. Lancet Psychiatry. 2018;5(4):327–38.

    Article  PubMed  Google Scholar 

  32. Lizano P, Lutz O, Xu Y, Rubin LH, Paskowitz L, Lee AM, et al. Multivariate relationships between peripheral inflammatory marker subtypes and cognitive and brain structural measures in psychosis. Mol Psychiatry. 2021;26(7):3430–43.

    Article  CAS  PubMed  Google Scholar 

  33. Laskaris L, Mancuso S, Shannon Weickert C, Zalesky A, Chana G, Wannan C, et al. Brain morphology is differentially impacted by peripheral cytokines in schizophrenia-spectrum disorder. Brain Behav Immun. 2021;95:299–309.

    Article  CAS  PubMed  Google Scholar 

  34. Alexander-Bloch AF, Shou H, Liu S, Satterthwaite TD, Glahn DC, Shinohara RT, et al. On testing for spatial correspondence between maps of human brain structure and function. Neuroimage. 2018;178:540–51.

    Article  PubMed  Google Scholar 

  35. Burt JB, Helmer M, Shinn M, Anticevic A, Murray JD. Generative modeling of brain maps with spatial autocorrelation. Neuroimage. 2020;220:117038.

    Article  PubMed  Google Scholar 

  36. Jackson AJ, Miller BJ. Meta-analysis of total and differential white blood cell counts in schizophrenia. Acta Psychiatr Scand. 2020;142(1):18–26.

    Article  CAS  PubMed  Google Scholar 

  37. Mazza MG, Capellazzi M, Lucchi S, Tagliabue I, Rossetti A, Clerici M. Monocyte count in schizophrenia and related disorders: a systematic review and meta-analysis. Acta Neuropsychiatr. 2020;32(5):229–36.

    Article  PubMed  Google Scholar 

  38. Sacks D, Baxter B, Campbell BCV, Carpenter JS, Cognard C, Dippel D, et al. Multisociety consensus quality improvement revised consensus statement for endovascular therapy of acute ischemic stroke. Int J Stroke. 2018;13(6):612–32.

    PubMed  Google Scholar 

  39. Melbourne JK, Rosen C, Chase KA, Feiner B, Sharma RP. Monocyte transcriptional profiling highlights a shift in immune signatures over the course of illness in schizophrenia. Front Psychiatry. 2021;12:649494.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Bishop JR, Zhang L, Lizano P. Inflammation subtypes and translating inflammation-related genetic findings in schizophrenia and related psychoses: a perspective on pathways for treatment stratification and novel therapies. Harv Rev Psychiatry. 2022;30(1):59–70.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Sealock JM, Lee YH, Moscati A, Venkatesh S, Voloudakis G, Straub P, et al. Use of the PsycheMERGE Network to investigate the association between depression polygenic scores and white blood cell count. JAMA Psychiat. 2021;78(12):1365–74.

    Article  Google Scholar 

  42. Tourjman V, Kouassi É, Koué M-È, Rocchetti M, Fortin-Fournier S, Fusar-Poli P, et al. Antipsychotics’ effects on blood levels of cytokines in schizophrenia: a meta-analysis. Schizophr Res. 2013;151(1–3):43–7.

    Article  PubMed  Google Scholar 

  43. Hashimoto N, Ito YM, Okada N, Yamamori H, Yasuda Y, Fujimoto M, et al. The effect of duration of illness and antipsychotics on subcortical volumes in schizophrenia: analysis of 778 subjects. Neuroimage Clin. 2018;17:563–9.

    Article  PubMed  Google Scholar 

  44. Ho B-C, Andreasen NC, Ziebell S, Pierson R, Magnotta V. Long-term antipsychotic treatment and brain volumes: a longitudinal study of first-episode schizophrenia. Arch Gen Psychiatry. 2011;68(2):128–37.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Gutiérrez-Galve L, Chu EM, Leeson VC, Price G, Barnes TR, Joyce EM, et al. A longitudinal study of cortical changes and their cognitive correlates in patients followed up after first-episode psychosis. Psychol Med. 2015;45(1):205–16.

    Article  PubMed  Google Scholar 

  46. Konopaske GT, Dorph-Petersen K-A, Pierri JN, Wu Q, Sampson AR, Lewis DA. Effect of chronic exposure to antipsychotic medication on cell numbers in the parietal cortex of macaque monkeys. Neuropsychopharmacology. 2007;32(6):1216–23.

    Article  CAS  PubMed  Google Scholar 

  47. Xie Y, Guan M, Wang Z, Ma Z, Wang H, Fang P, et al. rTMS induces brain functional and structural alternations in schizophrenia patient with auditory verbal hallucination. Front Neurosci. 2021;15:722894.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Horacek J, Brunovsky M, Novak T, Skrdlantova L, Klirova M, Bubenikova-Valesova V, et al. Effect of low-frequency rTMS on electromagnetic tomography (LORETA) and regional brain metabolism (PET) in schizophrenia patients with auditory hallucinations. Neuropsychobiology. 2007;55(3–4):132–42.

    Article  PubMed  Google Scholar 

  49. Jung WH, Kim JS, Jang JH, Choi J-S, Jung MH, Park J-Y, et al. Cortical thickness reduction in individuals at ultra-high-risk for psychosis. Schizophr Bull. 2011;37(4):839–49.

    Article  PubMed  Google Scholar 

  50. Sydnor VJ, Larsen B, Bassett DS, Alexander-Bloch A, Fair DA, Liston C, et al. Neurodevelopment of the association cortices: patterns, mechanisms, and implications for psychopathology. Neuron. 2021;109(18):2820–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Prof. Xin-Jiang Wang and Prof. Feng Cao for their continued scientific support and Dr. Guusje Collin and Prof. Yinchuan Jin for the additional information about neuroinflammation.

Funding

This work was supported by grants 82102138, 82271949, 82222031, 82220614, and 82171873 from the National Natural Science Foundation of China; grants 2020M683739 and 2019TQ0130 from the Project funded by China Postdoctoral Science Foundation; and grant 21YF1402400 from Shanghai Sailing Program.

Author information

Authors and Affiliations

Authors

Contributions

Guarantor of integrity of the entire study: L.-B.C. Study concepts/study design: L.-B.C., F.G., and X.C. Data acquisition: L.-B.C., Y.-F.F., X.-F.L., S.-W.Z., Y.-W.G., J.-W.F., and W.-JW. Data analysis/interpretation: all authors. Manuscript drafting or manuscript revision for important intellectual content: all authors. Approval of the final version of the submitted manuscript: all authors. Agrees to ensure any questions related to the work are appropriately resolved: all authors. Administrative or material support: L.-B.C. and H.Y. Literature research: L.-B.C., X.-Y.W., Y.-F.F., B.D.L., H.G., and X.C. Clinical studies: L.-B.C. and W.-J.W. Experimental studies: L.-B.C., X.-Y.W., Y.-F.F., X.-F.L., Y.W., S.-W.Z., Y.-W.G., J.-W.F., and X.C. Statistical analysis: Y.W. and X.C. Manuscript editing: L.-B.C., X.-Y.W., Y.-F.F., and X.C. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Long-Biao Cui, Fanglin Guan or Xiao Chang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of the First Affiliated Hospital of Fourth Military Medical University (Xijing Hospital) and complied with the Declaration of Helsinki.

Consent for publication

All authors have read the manuscript and provided consent for publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Methods.

Exclusion criteria. rTMS treatment. Brain transcriptional data. Table S1. Genes of interest related to inflammation and antipsychotic treatment response. Table S2. Shapiro-Wilk and Shapiro-Francia normality tests for inflammation markers. Table S3. Antipsychotics received at baseline for patients with schizophrenia. Table S4. Significant group by time interaction, group main effects, and time main effects on regional measurements among patients who were medicated at baseline (FDR significant results in bold font). Table S5. Correlation between structural changes and gene expression of antipsychotic treatment response and inflammation among patients who were medicated at baseline. Fig. S1. Change on PANSS score at baseline (T0) and follow-up period (T1). Fig. S2. Transcriptional level of gene sets of interest.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cui, LB., Wang, XY., Fu, YF. et al. Transcriptional level of inflammation markers associates with short-term brain structural changes in first-episode schizophrenia. BMC Med 21, 250 (2023). https://doi.org/10.1186/s12916-023-02963-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-023-02963-y

Keywords