Introduction
The current nosology classifies psychotic disorders in categories of non-affective (i.e. schizophrenia and schizoaffective disorder) and affective psychosis (i.e. bipolar and major depressive disorder with psychotic features) based on the type and course of symptomatology (American Psychiatric Association, 2013). However, the diagnosis of psychotic disorders at first is unstable longitudinally (Heslin et al., Reference Heslin, Lomas, Lappin, Donoghue, Reininghaus, Onyejiaka and Doody2015), and there is significant comorbidity between non-affective psychotic disorders with either bipolar (BD) or major depressive disorders (MDD) (Laursen, Agerbo, & Pedersen, Reference Laursen, Agerbo and Pedersen2009; Upthegrove, Marwaha, & Birchwood, Reference Upthegrove, Marwaha and Birchwood2016). Genome-wide association studies (GWAS) have also revealed shared genetic risks among these disorders (Cross-Disorder Group of the Psychiatric Genomics Consortium, 2019). Recent evidence suggests that the existing categorical classification system may be too simplistic to unravel the complex heterogeneity of the psychotic phenomena. Instead, non-categorical and transdiagnostic dimensional approaches are recommended (Quattrone et al., Reference Quattrone, Di Forti, Gayer-Anderson, Ferraro, Jongsma, Tripoli and Reininghaus2019; Quattrone et al., Reference Quattrone, Ferraro, Tripoli, Cascia, Quigley, Quattrone and Forti2021a, Reference Quattrone, Reininghaus, Richards, Tripoli, Ferraro, Quattrone and D'Andrea2021b). Item response modeling and factor analysis have been utilized to explore these non-categorical approaches. Indeed, psychotic disorders can be statistically conceptualized using a bifactorial model, with one general factor reflecting the mood-psychosis spectrum and a set of specific factors representing dimensions that can encompass positive, negative, disorganization, manic and depressive symptoms (Quattrone et al., Reference Quattrone, Di Forti, Gayer-Anderson, Ferraro, Jongsma, Tripoli and Reininghaus2019; Quattrone et al., Reference Quattrone, Ferraro, Tripoli, Cascia, Quigley, Quattrone and Forti2021a, Reference Quattrone, Reininghaus, Richards, Tripoli, Ferraro, Quattrone and D'Andrea2021b).
The transdiagnostic multidimensional model of psychopathology can be integrated into the continuity conceptualization of psychosis (van Os, Linscott, Myin-Germeys, Delespaul, & Krabbendam, Reference van Os, Linscott, Myin-Germeys, Delespaul and Krabbendam2009). The psychosis continuum or extended psychosis phenotype model (Linscott & van Os, Reference Linscott and van Os2010; Linscott & van Os, Reference Linscott and van Os2013) postulates that all psychotic phenomena exist on a continuum. At one end are individuals with psychotic disorders who experience severe and frequent psychotic symptoms. The healthier end of this continuum, however, constitutes individuals of the general population who experience attenuated psychotic symptoms expressed in the form of psychotic experiences (PEs) that do not fulfill the diagnosis criteria (Ragazzi et al., Reference Ragazzi, Shuhama, Sinval, Marôco, Corsi-Zuelli, Roza and Del-Ben2020). PEs are usually transient, although the prevalence is much higher (~7%) than the clinical phenotype (~3%), and they can become persistent to pose longitudinal risks for the development of psychotic disorders and common mental disorders (van Os & Reininghaus, Reference van Os and Reininghaus2016).
Several lines of evidence suggest that immune dysfunction is implicated in the pathophysiology of psychotic disorders (Corsi-Zuelli & Deakin, Reference Corsi-Zuelli and Deakin2021; Goldsmith, Rapaport, & Miller, Reference Goldsmith, Rapaport and Miller2016), MDD and BD (Çakici et al., Reference Çakici, Sutterland, Penninx, Dalm, de Haan and van Beveren2020; Goldsmith et al., Reference Goldsmith, Rapaport and Miller2016). A low level of non-resolving inflammation is already present in first-episode patients, even before psychopharmacological treatment initiation (Çakici et al., Reference Çakici, Sutterland, Penninx, Dalm, de Haan and van Beveren2020; Pillinger et al., Reference Pillinger, Osimo, Brugger, Mondelli, McCutcheon and Howes2019). In addition, schizophrenia GWAS show the strongest association with the major histocompatibility complex, which encodes proteins essential for the functioning of the immune system (Trubetskoy et al., Reference Trubetskoy, Pardiñas, Qi, Panagiotaropoulou, Awasthi, Bigdeli and van Os2022).
Only a few studies examined the contribution of blood immune proteins to clinical and subclinical symptom domains across the psychosis continuum, and most studies focused only on C-reactive protein (CRP) and interleukin (IL)-6, with mixed findings. For example, in early-stage psychosis, increased CRP levels were positively related to the severity of positive but not negative symptoms (Fernandes et al., Reference Fernandes, Steiner, Bernstein, Dodd, Pasco, Dean and Berk2016), and higher IL-6 levels were positively associated with negative but not positive symptoms (Stojanovic et al., Reference Stojanovic, Martorell, Montalvo, Ortega, Monseny, Vilella and Labad2014). However, these CRP and IL-6 findings were not replicated in a later study of early-stage schizophrenia (Krynicki et al., Reference Krynicki, Dazzan, Pariante, Barnes, Vincent and Roberts2021). In addition, depressive symptoms were positively associated (Noto et al., Reference Noto, Ota, Santoro, Ortiz, Rizzo, Higuchi and Brietzke2015), weakly associated (Krynicki et al., Reference Krynicki, Dazzan, Pariante, Barnes, Vincent and Roberts2021) or unrelated (Borovcanin et al., Reference Borovcanin, Jovanovic, Radosavljevic, Djukic Dejanovic, Bankovic, Arsenijevic and Lukic2012; Noto et al., Reference Noto, Maes, Nunes, Ota, Rossaneis, Verri and Noto2019) to increased IL-6 in recent-onset schizophrenia. The inconsistencies across studies could be related to differences in study design and methods, including assessment tools, sample sizes, age of assessment, and illness stage.
Given that PEs are common in the general population and share many etiological and pathophysiological backgrounds with psychotic disorders (Monshouwer et al., Reference Monshouwer, ten Have, Tuithof, van Dorsselaer, Bak, Gunter and de Graaf2022; Staines et al., Reference Staines, Healy, Coughlan, Clarke, Kelleher, Cotter and Cannon2022), examining their biological correlates could aid in comprehending the pathophysiology of the psychotic phenomena. This exploration also offers the potential to identify promising targets for early intervention, with the added benefit of reducing the influence of confounding factors related to illness duration and prolonged exposure to antipsychotics. However, the investigation of immune proteins and dimensions of PEs in non-clinical samples is scarce. A population-based study of adolescents found that raised CRP levels measured at age 16 were longitudinally associated with both positive and negative subclinical psychotic symptoms at age 17 (Khandaker et al., Reference Khandaker, Stochl, Zammit, Lewis, Dantzer and Jones2021). Another study from the same cohort found that raised IL-6 but not CRP levels measured during childhood at age nine correlated with negative symptom dimension at age 24 (Perry, Zammit, Jones, & Khandaker, Reference Perry, Zammit, Jones and Khandaker2021). We are not aware of other studies.
We aimed to investigate how immune proteins map onto a transdiagnostic multidimensional framework across the psychosis continuum in first-episode psychosis patients (FEP), their unaffected siblings, and community controls. Previous studies primarily focused on total scores of positive and negative symptoms, often neglecting immune proteins other than CRP or IL-6. In our study, we employed item response modeling and factor analysis to robustly estimate the transdiagnostic expression of psychosis. We then examined associations between the estimated transdiagnostic dimensions with a range of immune proteins, including trait (tumor necrosis factor (TNF)-α, interferon(IFN)-γ), state (IL-6, IL-1β), and regulatory (transforming growth factor (TGF)-β, IL-10, IL-4) (Capuzzi, Bartoli, Crocamo, Clerici, & Carrà, Reference Capuzzi, Bartoli, Crocamo, Clerici and Carrà2017; Miller, Buckley, Seabolt, Mellor, & Kirkpatrick, Reference Miller, Buckley, Seabolt, Mellor and Kirkpatrick2011), along with high-sensitive (hs)-CRP. To evaluate the role of disease endophenotypes without the interference of illness confounding factors, we considered first-degree unaffected siblings, who have an elevated risk of developing the disorder (ten-fold to 10% increased morbidity risk) (Lichtenstein et al., Reference Lichtenstein, Yip, Björk, Pawitan, Cannon, Sullivan and Hultman2009). Finally, given the high genetic liability to develop schizophrenia (~80%) (Sullivan, Kendler, & Neale, Reference Sullivan, Kendler and Neale2003), we explored if any significant associations would vary according to the degree of polygenic liability, as summarized into a schizophrenia polygenic risk score (SZ-PRS).
We predicted that higher levels of trait, state, and hs-CRP but lower regulatory markers would be associated with transdiagnostic dimensions along the psychosis continuum, with the unaffected siblings constituting an intermediate group between patients and controls. We also expected that significant associations would depend on a higher degree of SZ-PRS across the three groups.
Methods
We used data from the Schizophrenia and Other Psychosis Translational Research: Environment and Molecular Biology (STREAM) study, an incidence and case–sibling–control investigation conducted in the Ribeirão Preto catchment area (São Paulo, Brazil) (Del-Ben et al., Reference Del-Ben, Shuhama, Loureiro, Ragazzi, Zanatta, Tenan and Menezes2019), which integrates the international multicenter consortium EU-GEI (Gayer-Anderson et al., Reference Gayer-Anderson, Jongsma, Di Forti, Quattrone, Velthorst, de Haan and Morgan2020). The present study is cross-sectional, and we only included the Brazilian participants because blood cytokine measurement was a protocol specific from Brazil.
The study was approved by the local Research Ethics Committee (12606/2012), and written informed consent was obtained from all the participants. All data were collected by clinical psychologists and nurses with experience in mental health, with weekly supervision from senior psychiatrists.
Participants
All the STREAM participants were recruited from Ribeirão Preto city (São Paulo, Brazil) and the surrounding 25 municipalities during a three-year recruitment period (April 2012 until March 2015). The sample comprised FEP patients, their unaffected siblings, and community controls (Del-Ben et al., Reference Del-Ben, Shuhama, Loureiro, Ragazzi, Zanatta, Tenan and Menezes2019).
As part of the incidence study, we screened all individuals referred to mental health services with a suspected FEP within the Ribeirão Preto catchment area during the study period. Inclusion criteria were the following: (a) aged between 16–64 years; (b) presentation with a clinical diagnosis of FEP (International Classification of Diseases (ICD)-10 diagnoses F20-33), even if longstanding; (c) residing within the Ribeirão Preto catchment area during the study period. Exclusion criteria were: (a) previous contact with psychiatric services for psychosis; (b) FEP related to other medical conditions (ICD-10: F09); (c) transient psychotic symptoms resulting from substance intoxication/withdrawal (ICD-10: F1X.5); (d) severe learning disabilities (intelligent quotient <50 or diagnosis of intellectual disability, ICD-10: F70-79).
First-degree biological siblings of patients aged 16–64 years and residing within the same catchment area of cases were invited to participate in the study following the patients’ agreement. Community controls aged 16–64 years were recruited using a quota-sampling method according to the Brazilian Official Census Bureau 2010 (www.ibge.gov.br) with stratification by age and sex to ensure representativeness of the Ribeirão Preto catchment area's population at risk. Siblings of patients and community controls were only excluded if they had a lifetime history of psychosis or had ever been treated for psychosis. Per the EU-GEI protocol, those with a current or previous history of other mental disorders were not excluded to ensure the representativeness of the catchment area under study (Gayer-Anderson et al., Reference Gayer-Anderson, Jongsma, Di Forti, Quattrone, Velthorst, de Haan and Morgan2020).
For the present investigation, we discussed with a rheumatologist to identify medical conditions that could potentially influence the biological variables of interest. Consequently, we excluded participants who, either at the time of blood collection or during the past month, were receiving treatment with anti-inflammatory or immunosuppressive medications, were pregnant or lactating, and had chronic or acute immunological or inflammatory medical conditions. These excluded conditions encompassed but were not limited to nephropathy, urinary tract infection, dengue fever, rheumatic fever, human immunodeficiency virus, syphilis, Crohn's disease, throat infection, multiple sclerosis, pneumonia, hidradenitis suppurative; for details see (Corsi-Zuelli et al., Reference Corsi-Zuelli, Loureiro, Shuhama, Fachim, Menezes, Louzada-Junior and Del-Ben2020; Corsi-Zuelli et al., Reference Corsi-Zuelli, Marques, da Roza, Loureiro, Shuhama, Di Forti and Del-Ben2022a, Reference Corsi-Zuelli, Schneider, Santos-Silva, Loureiro, Shuhama, Menezes and Del-Ben2022b). To further minimize the effect of confounding by infection or chronic inflammatory illness, we excluded participants with hs-CRP levels >10 mg/L (FEP n = 14; unaffected sibling n = 6; community control n = 21), as suggested in the literature (Perry et al., Reference Perry, Zammit, Jones and Khandaker2021). A total of 134 FEP patients, 66 unaffected siblings, and 235 community controls had complete data on immune proteins and dimensions of psychopathology. However, genotype data were missing for 24 FEP patients, 14 unaffected siblings, and 35 community controls. These participants were excluded from the analyses investigating associations between symptom dimensions and immune proteins, leaving 362 participants.
Sociodemographic and clinical assessment
We used the Medical Research Council Sociodemographic Schedule (Mallet, Reference Mallet1997) for the acquisition of sociodemographic and clinical data, the Structured Clinical Interview for DSM-IV clinical version (Del-Ben et al., Reference Del-Ben, Vilela, Crippa, Hallak, Labate and Zuardi2001; First, Reference First1997) for diagnosis assessment, the Brief Psychiatry Rating Scale (BPRS) (Crippa, Sanches, Hallak, Loureiro, & Zuardi, Reference Crippa, Sanches, Hallak, Loureiro and Zuardi2001; Overall & Gorham, Reference Overall and Gorham1962) to measure illness severity, and the Nottingham Onset Schedule (Singh et al., Reference Singh, Cooper, Fisher, Tarrant, Lloyd, Banjo and Jones2005) for registration of psychosis onset date and pharmacological treatment starting date. Given the influence of smoking on inflammation (Yuan, Chen, Xia, Dai, & Liu, Reference Yuan, Chen, Xia, Dai and Liu2019), we used the Composite International Diagnostic Interview (Robins et al., Reference Robins, Wing, Wittchen, Helzer, Babor, Burke and Regier1988) to register the history of tobacco smoking during the last 12 months relative to the interview, and added this variable as a covariate in our analyses.
Subclinical and clinical dimensions of psychopathology
The Community Assessment of Psychic Experiences
For unaffected siblings of patients and community controls, three dimensions of PEs (positive, negative, and depressive) and a general dimension were evaluated using the Community Assessment of Psychic Experiences (CAPE) (Ragazzi et al., Reference Ragazzi, Shuhama, Sinval, Marôco, Corsi-Zuelli, Roza and Del-Ben2020; Stefanis et al., Reference Stefanis, Hanssen, Smirnis, Avramopoulos, Evdokimidis, Stefanis and Van Os2002). CAPE is a self-report instrument consisting of 42 items rated on a four-point Likert scale (1–4; never to nearly always) used to evaluate the frequency of psychotic experiences (PEs). Although it is a self-report instrument, the evaluation of PEs in our studied site occurred with the close supervision of a mental health specialist to ensure the quality of assessment.
The Operational Criteria System
For FEP patients, we used the Operational CRITeria (OPCRIT) system (McGuffin, Farmer, & Harvey, Reference McGuffin, Farmer and Harvey1991) to assess psychopathology and generate symptom dimensions (general, delusion, hallucination, negative, disorganization, manic, depressive). The OPCRIT system consists of a checklist with 59 items related to examining the mental state that can be filled using different sources, including clinical interviews or case records, allowing the establishment of diagnoses based on algorithms.
Blood immune proteins
Cytokines and hs-CRP were measured in plasma. For cytokines, we used panels assessing the inflammatory (IL-1β, IL-6, TNF-α, IFN-γ) and compensatory systems (TGF-β, IL-10, and IL-4), quantified using the Milliplex MAP human cytokine/chemokine magnetic bead panels. For more information, see Supplementary Material and our previous publications (Corsi-Zuelli et al., Reference Corsi-Zuelli, Loureiro, Shuhama, Fachim, Menezes, Louzada-Junior and Del-Ben2020; Corsi-Zuelli et al., Reference Corsi-Zuelli, Marques, da Roza, Loureiro, Shuhama, Di Forti and Del-Ben2022a, Reference Corsi-Zuelli, Schneider, Santos-Silva, Loureiro, Shuhama, Menezes and Del-Ben2022b). hs-CRP levels were measured using the latex immunoturbidimetric method according to the manufacturer's instructions with a detection limit of 0.1 mg/L (Wiener Laboratorios, Rosario, Argentina).
Schizophrenia Polygenic Risk Score
Genotyping of the present sample was carried out at the MRC Centre for Neuropsychiatric Genetics and Genomics in Cardiff (UK) using a custom Illumina HumanCoreExome-24 BeadChip array which included 570 038 genetic variants (Quattrone et al., Reference Quattrone, Ferraro, Tripoli, Cascia, Quigley, Quattrone and Forti2021a, Reference Quattrone, Reininghaus, Richards, Tripoli, Ferraro, Quattrone and D'Andrea2021b), and detailed in the Supplementary Material.
Statistical analyses
Demographic and clinical data
We used Pearson's Chi-square test to compare differences in categorical data (self-reported ethnicity, years of study, and tobacco smoking) across FEP patients, their unaffected siblings and community controls. For continuous data, quantile–quantile (Q-Q) plots and histograms were examined for inspection of the Gaussian Distribution. Parametric (CAPE) and non-parametric data (age, body mass index (BMI), and BPRS) were analyzed using a One-way Analysis of Variance or Kruskal–Wallis test, respectively, with Bonferroni post-hoc when appropriate.
Subclinical and clinical dimensions of psychopathology
The estimation of subclinical and clinical dimensions of psychopathology is fully described in previous publications of the EU-GEI consortium (Quattrone et al., Reference Quattrone, Di Forti, Gayer-Anderson, Ferraro, Jongsma, Tripoli and Reininghaus2019; Quattrone et al., Reference Quattrone, Ferraro, Tripoli, Cascia, Quigley, Quattrone and Forti2021a, Reference Quattrone, Reininghaus, Richards, Tripoli, Ferraro, Quattrone and D'Andrea2021b). Briefly, two bifactor models were estimated in Mplus, version 7.4 (Muthén & Muthén, Reference Muthén and Muthén2012) based on the associations among self-ratings of psychotic experiences in non-clinical samples (CAPE data) and observer ratings of psychotic symptoms in patients (OPCRIT data). Model fit statistics Log-Likelihood (LL), Akaike information criterion (AIC), Bayesian information criterion (BIC), and Sample-size Adjusted BIC (SABIC) of the bifactor solutions were compared with other solutions (i.e. unidimensional, multidimensional and hierarchical models). Then, the reliability and strength indices of the bifactor models were estimated using McDonald's omega (ω) (Rodriguez, Reise, & Haviland, Reference Rodriguez, Reise and Haviland2016), omega hierarchical (ω H)) (Rodriguez et al., Reference Rodriguez, Reise and Haviland2016), and index H (Hancock & Mueller, Reference Hancock and Mueller2001). Using the ‘FSCORE’ function in Mplus, factor scores were computed for (1) one general psychotic experience dimension and three specific dimensions of positive, negative, and depressive psychotic experiences for non-clinical samples; and (2) one general psychotic symptom dimension and six specific dimensions of delusions, hallucinations, negative, disorganization, manic, and depressive symptoms for patients. Please refer to (Quattrone et al., Reference Quattrone, Di Forti, Gayer-Anderson, Ferraro, Jongsma, Tripoli and Reininghaus2019; Quattrone et al., Reference Quattrone, Ferraro, Tripoli, Cascia, Quigley, Quattrone and Forti2021a, Reference Quattrone, Reininghaus, Richards, Tripoli, Ferraro, Quattrone and D'Andrea2021b) for further details.
Subclinical and clinical dimensions of psychopathology and immune proteins
We evaluated associations between subclinical and clinical dimensions of psychopathology with immune proteins separately in FEP patients, unaffected siblings, and community controls. For that, we used generalized linear models with robust estimators. In the FEP group, outcome variables were one general and six specific dimensions of delusion, hallucination, negative, disorganization, manic, and depressive symptoms derived from the OPCRIT. For unaffected siblings and community controls, outcomes variables were one general and three specific dimensions of positive, negative, and depressive experiences derived from the CAPE. Independent variables included data on eight immune proteins, both inflammatory (IL-1β, IL-6, TNF-α, IFN-γ, hsCRP) and regulatory (TGF-β, IL-10, and IL-4). Immune protein data were natural log-transformed and z-scored. Thus, regression coefficients and 95% confidence intervals (95% CI) represent changes per standard deviation (s.d.) of the exposure. All analyses were adjusted for sex, age, tobacco smoking, body mass index, SZ-PRS and five ancestry principal components. In the patient group, we additionally adjusted for the duration of pharmacological treatment.
To correct for multiple comparisons, we used the False Discovery Rate (FDR), Benjamin–Hochberg (5%) method, computed using R p.adjust function.
Schizophrenia Polygenic Risk Score (SZ-PRS)
To evaluate mean differences in SZ-PRS across the groups, values were standardized (z-scored), and we used a One-Way Analysis of Covariance (ANCOVA), including five ancestry principal components as covariates. To account for potential confounding, we additionally adjusted for age, sex, tobacco smoking, and BMI. All analyses concerning SZ-PRS were considered exploratory, given the limited sample size.
Exploratory analysis: Schizophrenia Polygenic Risk Score contribution to associations between dimensions of psychopathology and immune proteins
In the event of significant associations between clinical and subclinical dimensions of psychopathology with immune proteins, we explored how the associations would vary according to polygenic liability to schizophrenia. Study participants were stratified into terciles of SZ-PRS (top = high-risk group, middle = no genetic risk group, bottom = low genetic risk group), and the generalized linear models described before were re-run in each group.
Main analyses were performed using SPSS v.28 and R v1.3.1093, and all results were considered significant at p ⩽ 0.05 (two-sided).
Results
Sample characteristics
Differences among the groups on socio-demographic data are detailed in Table 1.
1Pearson's Chi-square test; 2Kruskal–Wallis (with Bonferroni post-hoc when appropriate).
Pairwise comparison: aFEP v. Siblings; bFEP v. Controls; cSiblings v Controls.
SD, Standard deviation; FEP, First-episode psychosis; BPRS, Brief Psychiatric Rating Scale.
Significant results are highlighted in bold.
The median of pharmacological treatment in the FEP group was thirteen weeks, and the patients were relatively stable at the time of assessment (BPRS mean [s.d.]: 8.8 [6.8]).
The Community Assessment of Psychic Experiences: group comparison
Unaffected siblings and community controls did not differ in any of the four CAPE dimensions (p > 0.05 for all). These results are detailed in online Supplementary Table 1.
Exploratory analysis: Schizophrenia Polygenic Risk Score group comparison
We found that FEP participants and their unaffected siblings did not differ in their SZ-PRS mean values. However, both FEP participants and their unaffected siblings had higher SZ-PRS than community controls. The results remained similar after including sex, age, tobacco smoking, and BMI in the model (Fig. 1).
Subclinical and clinical dimensions of psychopathology by immune proteins
Unaffected siblings
High IFN-γ (B [s.e.] = 0.84 [0.25]; FDR = 0.008) and low IL-4 levels (B [s.e.] = −0.28 [0.09] FDR = 0.008) were associated with the CAPE general dimension. Low IL-6 (B [s.e.] = −0.20 [0.08]; FDR = 0.044) and high TNF-α (B [s.e.] = 0.29 [0.11]; FDR = 0.044) were associated with the CAPE positive dimension. In addition, both low IL-6 (B [s.e.] = −0.22 [0.08]; FDR = 0.021) and IL-1β (B [s.e.] = −0.30 [0.09]; FDR = 0.008), but high IFN-γ levels (B [s.e.] = 0.54 [0.15]; FDR = 0.008) were related to the CAPE depressive dimension. There was a marginal association between low IL-10 levels and the CAPE negative dimension, but this did not survive FDR correction (B [s.e.] = −0.07 [0.11]; p = 0.052; FDR = 0.416). These results are presented in Figs 2A-D and online Supplementary Table S2.
Community controls
Similar to unaffected siblings, the CAPE general dimension was associated with high IFN-γ levels (B [s.e.] = 0.19 [0.07]; FDR = 0.024) in community controls. We also observed that IL-1β was negatively related to the general dimension, although this association did not survive FDR correction (B [s.e.] = −0.18 [0.08]; FDR = 0.124). In addition, low TGF-β related to the positive dimension (B [s.e.] = −0.14 [0.05]; FDR = 0.048). No other significant results were found in community controls. These results are presented in Figs 2A-D and online Supplementary Table S3.
First episode psychosis
We found that higher IFN-γ levels were associated with the manic symptom dimension (B [s.e.] = 0.29 [0.09]; FDR = 0.008) and lower TNF-α levels were associated with the depressive symptom dimension (B [s.e.] = −0.31 [0.10]; FDR = 0.008). Other associations did not survive the FDR correction, including a positive trend between IL-1β and the negative symptom dimension (B [s.e.] = 0.25 [0.10]; p = 0.013; FDR = 0.117), a positive and negative trend between TNF-α (B [s.e.] = 0.23 [0.11]; p = 0.031; FDR = 0.093) and IL-10 (B [s.e.] = −0.23 [0.10]; p = 0.017; FDR = 0.077) with the manic symptom dimension, respectively, and a positive trend between IL-6 and the depressive dimension (B [s.e.] = 0.17 [0.08]; p = 0.032; FDR = 0.144). These results are presented in Figs 3A-G and online Supplementary Table S4.
Exploratory analyses: Variability of associations by Schizophrenia Polygenic Risk Score
Previous associations that survived FDR correction were evaluated according to participants' SZ-PRS terciles and are detailed in online Supplementary Table 5.
Unaffected siblings
More than a third of siblings were classified with high polygenic liability (n = 18; 34.6%) and around the same proportion with low polygenic liability (n = 17; 32.7%).
Regarding the general dimension of PEs, we observed that the negative association with IL-4 only occurred in individuals at the top tercile of SZ-PRS (B [s.e.] = −0.21 [0.08]; p = 0.005) but not in those at the middle (p > 0.05) and bottom (p > 0.05) terciles. The IFN-γ association was significant and positively related to the general dimension in both those at the top (B [s.e.] = 0.53 [0.25]; p = 0.032) and bottom (B [s.e.] = 1.42 [0.57]; p = 0.012) terciles.
Regarding the positive dimension of PEs, associations with IL-6 differed according to the degree of SZ-PRS. Individuals at the top tercile presented a positive association (B [s.e.] = 0.18 [0.08]; p = 0.030), while a negative association was observed for those at the bottom tercile (B [s.e.] = −0.28 [0.14]; p = 0.048). Concerning TNF-α, significant positive associations were seen in those at the top tercile of SZ-PRS only (B [s.e.] = 0.59 [0.13]; p < 0.01).
Regarding the depressive dimension of PEs, associations with IFN-γ differed according to the degree of SZ-PRS. Those at the top tercile presented a positive association (B [s.e.] = 0.46 [0.21]; p = 0.029), while a negative association was observed for those at the bottom tercile (B [s.e.] = −1.73 [0.50]; p < 0.01). The IL-6 and IL-β associations with the depressive dimension were irrespective of the degree of SZ-PRS.
Community controls
More than half of community controls were classified as having a low polygenic liability to schizophrenia (n = 110; 55%), and the minority with high polygenic liability (n = 32; 16%).
In community controls, the previously observed association between IFN-γ and CAPE general dimension was only significant in those at the top tercile of SZ-PRS (B [s.e.] = 0.66 [0.14]; p < 0.01) but not in those in the middle or bottom terciles (p > 0.05 for both).
However, the previously negative association observed for TGF-β and the positive dimension was irrespective of the degree of SZ-PRS.
First-episode psychosis
More than a third of patients were classified with high polygenic liability (n = 38; 34.5%), and around the same proportion with low polygenic liability (n = 33; 30%).
The observed association between IFN-γ and manic symptoms dimension was significant and positive for those at the top (B [s.e.] = 0.29 [0.10]; p = 0.003) and bottom (B [s.e.] = 0.26 [0.13]; p = 0.040) terciles. TNF-α association with the OPCRIT depressive dimension was irrespective of the degree of SZ-PRS.
Discussion
Our study sought to investigate the involvement of immune proteins in the multidimensional expression of psychosis and the contribution of schizophrenia polygenic liability in the associations. We provide evidence that both inflammatory and regulatory immune proteins map onto dimensions of psychopathology in both clinical and community samples with some variability by polygenic liability to schizophrenia. We also noticed that among non-clinical samples, most associations were seen in first-degree relatives of patients, supporting a role for familial liability. The results align with the view that psychosis exists on a continuum and that PEs experienced by the general population have similar biological substrates to clinical disorders. Our findings also reinforce that immune dysregulation goes beyond the traditional nosology by impacting the pathogenesis of symptom domains that cut across diagnostic categories, as reflected in associations between immune proteins with the depressive dimension in FEP patients and unaffected siblings and manic dimension in patients.
Immune dysfunction in the continuity of psychosis
Meta-analyses of cross-sectional studies have demonstrated elevated levels of cytokines and other inflammatory proteins in the blood of patients with psychosis and depression, compared with controls (Çakici et al., Reference Çakici, Sutterland, Penninx, Dalm, de Haan and van Beveren2020; Goldsmith et al., Reference Goldsmith, Rapaport and Miller2016). Large prospective cohort studies indicate that heightened concentrations of IL-6 during childhood or adolescence precede the onset of symptoms and are associated with an increased risk for both disorders in adulthood (Giollabhui, Ng, Ellman, & Alloy, Reference Giollabhui, Ng, Ellman and Alloy2021; Khandaker, Pearson, Zammit, Lewis, and Jones, Reference Khandaker, Pearson, Zammit, Lewis and Jones2014). Nevertheless, there is scarce evidence regarding the association between a immune protein profile and specific patterns of symptomatology across the psychosis spectrum. Our bifactorial model allowed us to capture for the first time that alterations in IL-6 and immune proteins upstream the IL-6 pathway (e.g. IL-1β and TNF-α) extend to subclinical positive and depressive symptoms, particularly in vulnerable individuals such as first-degree relatives, confirming a transdiagnostic role for IL-6 and its pathways along the continuum.
Siblings of patients may represent endophenotypes of the disease because they can display biological dysfunction similar to the probands and subclinical psychiatric symptoms without fulfilling the diagnostic criteria (Lichtenstein et al., Reference Lichtenstein, Yip, Björk, Pawitan, Cannon, Sullivan and Hultman2009; Noyan et al., Reference Noyan, Erdağ, Tüzün, Yaylım, Küçükhüseyin, Hakan and Üçok2021). One important consideration is that, in our study, the unaffected siblings did not differ from community controls in the mean scores for general or specific dimensions of PEs. A hypothetical explanation is that the recruited siblings were the caregivers and, therefore, healthy relatives. However, our findings are consistent with a study showing that rates of PEs do not significantly differ between relatives and non-relatives of probands (Landin-Romero et al., Reference Landin-Romero, McKenna, Romaguera, Álvarez-Moya, Sarró, Aguirre and Pomarol-Clotet2016). On the other hand, our SCZ-PRS findings confirm a significant polygenic liability to schizophrenia among siblings, evidenced by higher means in patients and their siblings compared to controls, consistent with prior research (Ahangari et al., Reference Ahangari, Gentry, Nguyen, Kirkpatrick, Verrelli and Riley2022a; Ahangari et al., Reference Ahangari, Kirkpatrick, Nguyen, Gillespie, Kendler, Bacanu and Riley2022b; Ohi et al., Reference Ohi, Nishizawa, Shimada, Kataoka, Hasegawa, Shioiri and Ikeda2020, Reference Ohi, Nishizawa, Sugiyama, Takai, Fujikane, Kuramitsu and Shioiri2022).
The only previous study that simultaneously assessed clinical and subclinical symptom dimensions and their association with immune proteins found elevated IL-1β levels associated with the severity of psychotic symptoms in established schizophrenia, with a similar trend in their siblings. Different from our study, these measurements were obtained through total scores of the Comprehensive Assessment of Symptoms, for evaluation of clinical symptoms in patients, and the History and the Structured Interview for Schizotypy-Revised, for the evaluation of subclinical symptoms in siblings (Noyan et al., Reference Noyan, Erdağ, Tüzün, Yaylım, Küçükhüseyin, Hakan and Üçok2021). In our study, the negative association between IL-1β with the subclinical depressive dimension, and between IL-6 with the subclinical positive and depressive dimensions in vulnerable siblings may offer distinct yet complementary interpretations. One is that IL-1β and IL-6 are designated state immune markers (Capuzzi et al., Reference Capuzzi, Bartoli, Crocamo, Clerici and Carrà2017; Miller et al., Reference Miller, Buckley, Seabolt, Mellor and Kirkpatrick2011). Hence, their levels are not stable and may fluctuate along the continuum, with associations turning positive only when transiting to full-blown psychosis or during acute states. In fact, we identified positive associations between IL-6 levels and depressive symptoms and between IL-1β and negative symptoms in FEP (although these did not survive FDR correction). One important limitation in our patient sample is that, although these were early-stage psychosis, they were not free of antipsychotic treatment. Meta-analyses show that levels of state markers such as IL-6 and IL-1β tend to decrease with antipsychotic treatment (Capuzzi et al., Reference Capuzzi, Bartoli, Crocamo, Clerici and Carrà2017; Miller et al., Reference Miller, Buckley, Seabolt, Mellor and Kirkpatrick2011). The other interpretation may relate to different actions and pathways of these immune proteins. For example, IL-6 ‘trans-signalling’ involves IL-6 binding to its soluble receptor and is responsible for most of its pathogenic actions (Rose-John, Winthrop, & Calabrese, Reference Rose-John, Winthrop and Calabrese2017). However, it is widely known that IL-6 has physiological roles for neuronal development and function mediated mainly through its protective membrane-bound ‘classic signalling’ (Kummer, Zeidler, Kalpachidou, & Kress, Reference Kummer, Zeidler, Kalpachidou and Kress2021). The negative associations in siblings could, therefore, reflect a failure in the protective IL-6 pathway, which may function as a gateway for future dysregulation of the pathogenic ‘trans-signalling’ pathway and subsequent manifestation of full-blown psychosis (Hartwig, Borges, Horta, Bowden, & Davey Smith, Reference Hartwig, Borges, Horta, Bowden and Davey Smith2017).
Different from previous studies, we were also able to evaluate the contribution of other immune proteins along the continuum. An interesting discovery was the positive association between the general dimension with IFN-γ in both unaffected siblings and community controls. This association was stronger for first-degree relatives, expanding to the depressive dimension in this group. The general dimension is the covariance among all items and can be a useful representation of the mood-psychosis spectrum (Quattrone et al., Reference Quattrone, Di Forti, Gayer-Anderson, Ferraro, Jongsma, Tripoli and Reininghaus2019). Our findings seem to align with a hypothesis that IFN-γ may be a potential trait marker for psychosis, as longitudinally evaluated in drug-naïve FEP and acute exacerbations (Capuzzi et al., Reference Capuzzi, Bartoli, Crocamo, Clerici and Carrà2017; Miller et al., Reference Miller, Buckley, Seabolt, Mellor and Kirkpatrick2011). TNF-α, which we found to be positively associated with the positive dimension in siblings, has also been designated as a trait marker for psychosis when evaluated in drug-naïve FEP and acute exacerbations (Capuzzi et al., Reference Capuzzi, Bartoli, Crocamo, Clerici and Carrà2017; Miller et al., Reference Miller, Buckley, Seabolt, Mellor and Kirkpatrick2011). Interestingly, in our minimally treated and relatively stable FEP patients, only the two previously described trait markers remained significant after FDR correction, IFN-γ with the manic dimension and TNF-α with the depressive dimension. These significant findings between cytokines and affective dimensions in FEP align with evidence that psychotic disorders, MDD, and BD share commonalities in genetic, inflammatory and clinical features that converge in a transdiagnostic conceptualization (Wei et al., Reference Wei, Womer, Sun, Zhu, Sun, Duan and Wang2022).
Finally, we found a negative association between IL-4 and the general dimension in siblings, in addition to a negative association between TGF-β and the positive dimension in controls. These findings would align with a T-cell imbalance hypothesis of psychosis, which claims that T helper lymphocytes type 2 (Th-2)-producing anti-inflammatory IL-4 (Kim et al., Reference Kim, Myint, Lee, Han, Lee, Kim and Leonard2004) and protective regulatory T cells producing TGF-β (Corsi-Zuelli & Deakin, Reference Corsi-Zuelli and Deakin2021; Corsi-Zuelli et al., Reference Corsi-Zuelli, Deakin, de Lima, Qureshi, Barnes, Upthegrove and Del-Ben2021) fail to counteract the inflammatory actions of pathogenic Th1 cells and monocytes, resulting in disinhibited production of trait cytokines such as IFN-γ and TNF-α, and state cytokines such as IL-6 and IL-1β.
Altogether, our results indicate that innate immune proteins related to the IL-6 pathway (IL-6, TNF-α, IL-1β) and adaptive T cell-related proteins (IFN-γ, IL-4, TGF-β) may be biological contributors to the multidimensional mood-psychosis spectrum. Our hypothetical model is summarized in Fig. 4.
Polygenic threshold and the immune-psychosis spectrum
Our exploratory analyses suggested that associations between immune proteins and symptom dimensions significantly varied according to the degree of polygenic liability to schizophrenia across the three studied subgroups. This variability occurred for both innate and adaptive cytokines but was mainly observed for the psychotic dimension, in which most of the significant associations occurred primarily in those at the top tercile of SZ-PRS. Intriguing, in the sibling group, those at the highest genetic risk displayed a positive association between the positive dimension and both IL-6 and TNF-α, whereas those at the lowest genetic risk had a negative association with IL-6 in this dimension. These different directions of associations across the genetic groups – with positive associations only seen in the highest genetic group – support our previously discussed hypothesis of positive associations reflecting dysregulation in cytokines’ pathogenic pathways only in those at higher risk for the disease. Longitudinal studies are needed to clarify if siblings with high SZ-PRS and raised immune proteins are those who later transition to psychosis.
Associations with the affective dimensions tended to be irrespective of SZ-PRS, particularly between the depressive dimension and innate cytokines (IL-6, IL-1β, and TNF-α). These findings suggest that positive and general dimensions of psychosis may have a greater genetic and familial basis in their association with both innate and adaptive immune proteins than affective dimensions. The independence of the depressive dimension from SZ-PRS in its association with innate cytokines may indicate the influence of other but still unknown factors that may be shared across disorders and impact the functioning of the innate immune system (e.g. environmental factors such as stress) (Corsi-Zuelli et al., Reference Corsi-Zuelli, Marques, da Roza, Loureiro, Shuhama, Di Forti and Del-Ben2022a, Reference Corsi-Zuelli, Schneider, Santos-Silva, Loureiro, Shuhama, Menezes and Del-Ben2022b).
Finally, we noted an interdependency of the adaptive cytokine IFN-γ to SZ-PRS in its positive association with the general dimension (community controls and siblings) and affective dimension (particularly manic in FEP and depression in siblings). Noteworthy, in the control group, for which the sample size (n = 200) was slightly greater compared to FEP patients (n = 110) and siblings (n = 52), significant associations were only found for those at the top tercile of SZ-PRS, but not in those in the middle or bottom terciles. In addition, in the FEP group, associations for those in the top tercile were slightly greater than those at the bottom tercile. While these findings are exploratory, they appear to align with the concept of IFN-γ being a trait marker in psychosis, irrespective of the stage of illness. Our study now offers novel evidence that extends to the transdiagnostic continuity spectrum.
Altogether, our exploratory results suggest that genetic liability to schizophrenia influences the immune phenotype of transdiagnostic dimensions in psychosis.
Strengths and limitations
Our study is the first to investigate associations between immune proteins along the extended psychosis phenotype using multidimensional models to help unravel the heterogeneity of the psychotic phenomena. Instead of relying on total scores, our item response modeling and factor analysis enabled us to generate one general factor representing the mood-psychosis spectrum, and a set of specific factors to estimate the transdiagnostic expression of psychosis. Many of the generated factors had not been evaluated in their associations with immune proteins before. In addition, previous studies were mainly focused on CRP and IL-6. We expanded that by including both regulatory and inflammatory cytokines. Moreover, we examined for the first time the contribution of polygenic liability to schizophrenia in the immune-psychosis spectrum. Finally, our community controls were recruited based on the representativeness of the catchment area's population at risk, and we additionally recruited unaffected siblings of patients to investigate disease endophenotypes.
Important limitations should be acknowledged. First, although we considered first-stage individuals with psychosis to minimize the masking effects of long-term treatment and illness duration, the patients were not free of pharmacological treatment. The inclusion of their first-degree siblings was an advantage to that limitation, although we recognize that our sibling sample was smaller. Furthermore, even though we used transdiagnostic dimensional approaches, we recognize limitations. For example, we could not capture the full extent of the negative symptom dimension with the OPCRIT and CAPE. The structure of negative symptoms comprises at least two subdomains of diminished (a) expressivity and (b) motivation and pleasure, each encompassing distinct mechanisms (Strauss et al., Reference Strauss, Horan, Kirkpatrick, Fischer, Keller, Miski and Carpenter2013), and that would require more specific assessment scales. Finally, we treated the SZ-PRS analyses as exploratory because our total sample size was not extensive. There are several important considerations for future research. Firstly, the findings from this study should be replicated and validated in larger samples. Item response modeling and factor analysis have proven valuable in our investigation, and their application in future studies could help establish the robustness and generalizability of the observed associations. Moreover, longitudinal studies are essential to shed light on the progression of psychosis and the potential transition from non-clinical to clinical states.
Conclusions and future perspectives
Our study delves into the intricate relationship between the immune system and the multidimensional expression of psychosis. By examining both clinical and non-clinical samples, we have identified significant associations between innate and adaptive inflammatory and regulatory immune proteins with the transdiagnostic multidimensional expression of psychosis, observing some variability by polygenic liability to schizophrenia. Our research contibutes to the understanding of the extended psychosis phenotype, oppening new avenues for further exploration.
Comprehending the trajectory of individuals with high polygenic risk scores and specific immune profiles could be an important step in predicting and preventing the onset of psychosis. Understanding the immune system, particularly at early stages or before illness onset, may shed light on mechanisms of disease pathophysiology to improve diagnosis, prognosis and development of novel personalized treatments.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0033291724000199.
Acknowledgements
We acknowledge Giuliana Bertozi for the technical support during cytokine measurement.
Authors’ contributions
F.C.Z. interpreted critically all the results of the study for important intellectual content to write a first draft of the manuscript; revised the manuscript critically; performed statistical analyses; measured cytokines in human plasma with technical supervision; created the figures and tables and responded to reviewers’ questions. D.Q revised the manuscript critically, generated clinical and subclinical dimensions and computed the SZ-PRS for statistical analyses. C.M.L. and T.C.C.R. critically revised the manuscript and participated in data collection. R.S. contributed to the management of data collected and ethical aspects. P.R.M., P.L.J., and C.M.D.B. conceptualized the clinical work and critically revised the manuscript.
Funding statement
This work received financial support from the São Paulo Research Foundation – FAPESP, Brazil (grant number 2012/05178-0), the National Council for Scientific and Technological Development – CNPq, Brazil (grant number 476945/2012-7), and the Center for Research in Inflammatory Diseases (CRID grant number 2013/08216-2). P.L.-J., C.M.D-B., and P.R.M. are recipients of fellowships from CNPq; F.C-Z, C.M.L, and R.S. acknowledge fellowships from FAPESP (21/07448-3 and 19/13229-2, 20/15752-1, and 2013/11167-3 respectively). DQ is a recipient of a grant from MRC/UKRI (MRC CARP grant MR/W030608/1).
Competing interests
None.