Introduction
Treatment-resistant schizophrenia (TRS) refers to approximately one-third of schizophrenia patients who do not adequately alleviate their psychotic symptoms despite standard antipsychotic treatment [Reference Conley and Kelly1]. Patients with TRS have higher rates of unemployment, poorer quality of life, and poorer social and occupational functioning than individuals who respond to treatment [Reference Iasevoli, Giordano, Balletta, Latte, Formato and Prinzivalli2]. Clozapine is the only antipsychotic recommended for TRS, which is effective in about 60% of cases [Reference Kumra, Kranzler, Gerbino-Rosen, Kester, DeThomas and Kafantaris3] and improves indicators of morbidity and mortality [Reference Conley and Kelly1]. Nevertheless, studies have shown that identification difficulties with TRS led to delays in clozapine prescription, which in turn was associated with reduced responsiveness of patients to clozapine [Reference Shah, Iwata, Brown, Kim, Sanches and Takeuchi4]. This makes early identification of TRS critical and the ascertainment of biomarkers of TRS a priority for the field of schizophrenia research.
Evidence suggests that earlier age of schizophrenia onset is a robust predictor of TRS, and that male gender, autumn/winter birth, poor premorbid functioning, and rural upbringing may also contribute [Reference Smart, Kepinska, Murray and MacCabe5]. However, a gap in the literature exists in the genetics of TRS, particularly biomarkers. To date, there is considerable heterogeneity in the genetic findings associated with TRS. A family history of schizophrenia is likely linked to developing TRS [Reference Kowalec, Lu, Sariaslan, Song, Ploner and Dalman6]. Candidate gene research investigating the involvement of specific targets in TRS mostly clustered around the serotoninergic and dopaminergic systems, as well as on systems involved in oxidative stress and inflammation [Reference Pisanu and Squassina7]. Conversely, a small sample study of TRS, defined by American Psychiatric Association criteria, did not find any significant association among the 384 candidate loci [Reference Teo, Zai, Borlido, Tomasetti, Strauss and Shinkai8]. The collective interpretation of these results is made difficult by the slightly different recruitment TRS criteria [Reference Hodgson, McGuffin and Lewis9]. The underrepresentation of individuals with treatment-resistant psychiatric symptoms in studies similarly reduced statistical efficiency [Reference Hofer, Radner, Edlinger, Kemmler, Rettenbacher and Fleischhacker10]. Another critical question in the TRS field is whether TRS represents a more severe form of schizophrenia, or if it represents a distinct subtype of schizophrenia with a different symptom profile and pathophysiology. Some clinical, imaging, biological, and genetic evidence supports the latter [Reference Nucifora, Woznica, Lee, Cascella and Sawa11, Reference Pardiñas, Smart, Willcocks, Holmans, Dennison and Lynham12]. Given the potentially complex genetic architecture of this trait and the discrepancies in the clinical definition of TRS, there remains disagreement on the best approach to maximize the informativity and power of genetic studies on TRS. Our work was designed as a proteomics analysis to identify functional brain proteins that distinguish TRS from schizophrenia, in an attempt to provide a new perspective on this issue by exploring the expression products.
Proteomic techniques are increasingly being used to screen for biomarkers of schizophrenia, while no studies have yet applied proteomics to TRS. Proteome-wide association study (PWAS) captures any variant affecting the coding regions of genes, and then assigns each protein-coding gene functional affecting scores, and is widely utilized to robustly prioritize candidate genes [Reference Wingo, Gerasimov, Liu, Duong, Vattathil and Lori13]. PWAS is a high-throughput approach where proteomic studies detect fewer expressed proteins than expressed genes detected by the transcriptome, but protein expression provides an accurate functional profile and reflects the complex interactions between genes and the environment, presenting an unbiased picture of the current physiological state. The importance of those interactions has been increasing in the research of neurological diseases [Reference Hofer, Radner, Edlinger, Kemmler, Rettenbacher and Fleischhacker10]. We performed two independent PWAS to validate each other, and functional enrichment and annotation analysis was conducted based on the results. Colocalization analysis was performed to identify variant loci that were causally related to the expressed proteins. These methods may identify biomarkers and discern the specific mechanisms underlying TRS, thereby offering proof of its classification as a subtype of schizophrenia and facilitating the early detection of TRS.
Materials and Methods
GWAS summary data
The GWAS summary data of TRS and non-treatment-resistant schizophrenia (non-TRS) were derived from a recently published study [Reference Pardiñas, Smart, Willcocks, Holmans, Dennison and Lynham12]. Patients with TRS were derived from the CLOZUK1 and CLOZUK2 cohorts, with a total sample size of 10,501 individuals [Reference Pardinas, Holmans, Pocklington, Escott-Price, Ripke and Carrera14]. All patients with TRS were prescribed clozapine following at least two failed trials of antipsychotics and in accordance with National Institute for Health and Care Excellence guidelines for TRS. Patients with non-TRS were derived from 34 studies that were included in the meta-analysis by the Schizophrenia Working Group of the Psychiatric Genomics Consortium (PGC), with a total sample size of 20,325 patients with schizophrenia [15]. Fourteen studies with clinical records identified and excluded individuals with TRS, and cases from the remaining 20 studies without clinical records were conservatively included in our analysis as individuals with non-TRS. Owing to the number of different datasets and genotyping arrays involved in this analysis, the processing of TRS and non-TRS GWAS samples was performed separately on the data generated by the original study [Reference Choi and O’Reilly16, Reference Tiihonen, Mittendorfer-Rutz, Majak, Mehtälä, Hoti and Jedenius17]. Both of these imputations used the SHAPEIT/IMPUTE2 pipeline. Single Nucleotide Polymorphisms (SNPs) were restricted to minor allele frequencies (MAF) of 5% or higher and called in at least 20,000 combined samples, and any strand-ambiguous markers (A/T and G/C) with MAF ≥ 40% were discarded in both datasets.
Human brain proteome reference weights for PWAS
Two human brain proteome reference weight datasets were obtained from recent publicly available studies. The discovery dataset was derived from Religious Orders Study and Rush Memory and Aging Project (ROS/MAP) cohorts, recruiting 391 individuals from two longitudinal clinical-pathologic cohort studies of aging and Alzheimer’s disease (AD) [Reference Wingo, Fan, Duong, Gerasimov, Dammer and Liu18]. After quality control, 8,356 proteins from 376 subjects were included in our analysis. Among these, 262 were female and the average age at death was 89 years. The final clinical diagnosis included no cognitive impairment, mild cognitive impairment (MCI), AD dementia, or other causes of dementia. The confirmation dataset was derived from the Banner Sun Health Research Institute (Banner) containing 198 individuals [Reference Wingo, Liu, Gerasimov, Gockley, Logsdon and Duong19]. After quality control, 11,518 proteins from 152 subjects were quantified. Of these, 87 were female and the average age at death was 85 years. Only individuals with a final diagnosis of AD or normal cognition were included. Both of the proteomic reference datasets utilized the same quality control procedures to identify and control the effects of clinical covariates (i.e., age, gender, and final clinical diagnosis of cognitive status) before estimating protein weights. Details can be found in the Supplementary aterial.
Proteome-wide association study analysis
PWAS analysis was performed by integrating the TRS GWAS data with two brain proteomes using the FUSION pipeline (http://gusevlab.org/projects/fusion/). Briefly, we utilized FUSION to calculate protein weights in both the discovery and confirmation datasets individually. Subsequently, we combined the genetic effect of TRS (TRS GWAS z-score) with precomputed protein weights by calculating the linear sum of z-score × weight of independent SNPs to perform the PWAS of TRS. Only proteins identified in the discovery dataset and replicated in the confirmation dataset were considered TRS-associated proteins. The linkage disequilibrium reference panel routinely utilized 1,190,321 HapMap SNPs from 489 individuals of European descent from the 1000 Genomes Project in FUSION. We implemented 2,000 permutations to control the potential impact of multiple testing on PWAS results. The proteins with permutated P < 0.05 were considered significant. The design is presented in the Supplementary Figure S1.
Colocalization analysis
Our research performed a colocalization analysis of all genes identified by the two-stage PWAS using the coloc R package. We evaluated the colocalization status of a gene by calculating the PP that the genetic and functional associations derived from a shared causal SNP (PP4). Genes with PP4 > 0.75 are considered to be colocalized.
Functional enrichment and annotation analysis
The GO annotation and KEGG pathway enrichment analyses of the genes identified by PWAS were performed by the DAVID tool (https://david.ncifcrf.gov/). GO enrichment analysis includes biological process (BP), cellular component (CC), and molecular function (MF) analysis. KEGG database is a bioinformatics resource for mining significantly metabolic pathways enriched in the gene list. The purpose was to extract important GO terms and KEGG pathways, which can depict the characteristics of TRS. The results were considered statistically significant if P < 0.05.
Results
PWAS results of TRS
In PWAS, two significant proteins were identified in the discovery reference dataset and replicated in the confirmation reference dataset, including CPT2 (P PWAS-ROS/MAP = 4.15 × 10−2 and P PWAS-Banner = 3.38 × 10−3) and APOL2 (P PWAS-ROS/MAP = 4.49 × 10−3 and P PWAS-Banner = 8.26 × 10−3). Specifically, in the discovery stage, a total of 24 TRS-related candidate proteins were identified using the ROS/MAP reference dataset, such as PRDX1 (P PWAS-ROS/MAP = 1.00 × 10−3), APOL2 (P PWAS-ROS/MAP = 4.49 × 10−3), and WARS2 (P PWAS-ROS/MAP = 6.41 × 10−3). In the confirmation stage, a total of 19 TRS-related candidate proteins were identified using the Banner reference dataset, such as CPT2 (P PWAS-Banner = 3.38 × 10−3) and ABCC1 (P PWAS-Banner = 4.25 × 10−3). Statistically significant genes identified in PWAS analysis are shown in Table 1 and Figure 1.
Significant proteins identified by PWAS.
Note: Only significant permutation P values are presented.
Abbreviation: PWAS, proteome-wide association study.
Manhattan plots of significant human brain proteins identified in proteome-wide association study (PWAS). Each point corresponds to a single test of association between a gene and phenotype, plotted according to genomic position on the x-axis and the strength of association (−log10 P-value) on the y-axis. Two common statistically significant proteins resulting from the analysis were mapped out.
Colocalization analysis
Of all 41 proteins identified by PWAS, colocalization analysis identified three genes that are causal for TRS and encoded functional proteins, including CCDC91 (PP4 = 0.981), PRDX1 (PP4 = 0.894), and WARS2 (PP4 = 0.757). The results of the colocalization analysis of the genes identified by PWAS are presented in Table 1.
Functional enrichment and annotation analysis
GO enrichment analysis of all 41 genes identified by PWAS results are shown in Table 2. DAVID detected 14 GO terms, such as mitochondrion (GO:0005739, P < 0.0001), Golgi apparatus (GO:0005794, P = 0.0145), oxidoreductase activity (GO:0016620, P = 0.0338), and protein domain specific binding (GO:0019904, P = 0.0117). For pathway enrichment analysis of the genes identified by PWAS, DAVID detected only one candidate pathway for TRS, metabolic pathways (hsa01100, P = 0.0312).
GO enrichment analysis results of TRS-associated genes identified by PWAS.
Abbreviation: PWAS, proteome-wide association study.
Discussion
In this work, we performed TRS-associated PWAS and found a total of 41 proteins that were differentially expressed, 2 of which were duplicated in the PWAS analysis of the discovery and confirmation dataset, suggesting a prominent role in the BPs of TRS, including CPT2 and APOL2. Owing to polygenic inheritance, a complex trait is often influenced by multiple genes with similar functions as annotated in gene pathways. We extended PWAS results from gene-based analysis to pathway-based analysis, identifying 14 GO terms and 1 candidate pathway for TRS. Aiming to understand the mechanisms driving GWAS risk loci, we performed a colocalization analysis of PWAS-identified genes and identified three variants that were causally related to expression in the human dorsolateral prefrontal cortex.
TRS may be associated with abnormalities in the β-oxidative metabolic pathway of long-chain fatty acids in mitochondria. We discovered a novel TRS-related protein CPT2. Carnitine palmitoyl transferases 2 (CPT2) is the core protein of a catalytically active multiprotein complex localized in the inner mitochondrial membrane, assisting long-chain fatty acids to enter the mitochondrial matrix for oxidation and energy. The role of the lipid regulatory system on TRS has been supported by studies. One work found that clozapine altered the activity of the AMPK-ACC-CPT1 pathway, a central pathway of lipid metabolism, to affect the lipid compositions of the neuronal membranes in the rat frontal cortex [Reference Kim, Kim, Yu, Park, Kang and Ahn20]. Abnormalities in membrane lipid composition have also been reported in the frontal cortex of patients with schizophrenia [Reference Narayan, Head, Gilmartin, Dean and Thomas21]. Beta-oxidation may be potentially linked to the pathogenesis of TRS. The metabolite of clozapine is capable of interacting with a wide range of neurotransmitter receptors, suggesting that TRS may have a neurobiological etiology [Reference Pardinas, Owen and Walters22]. One study found that deletion of CPT2 in the nervous system leads to elevated expression of β-oxidation enzymes [Reference White, Lee, Choi, Chu, Scafidi and Wolfgang23]. Individuals with genetic disorders in mitochondrial fatty acid β-oxidation may suffer from neurological disorders, including seizures, encephalopathies, and cortical atrophy [Reference Tein, Vajsar, MacMillan and Sherwood24–Reference Merritt, Norris and Kanungo26]. CPT activity is present in almost all brain regions especially the brainstem [Reference Virmani, Pinto, Bauermann, Zerelli, Diedenhofen and Binienda27], and carnitine shuttle and β-oxidation genes are expressed primarily in astrocytes and neural stem cells [Reference Xie, Jones, Deeney, Hur and Bankaitis28], suggesting that CTP2 deficiency may involve the central nervous system. Our GO and pathway enrichment analysis also support the CTP2-centred bio-metabolic processes.
TRS may be linked to cholesterol transport and homeostasis. The other brain protein we identified was apolipoprotein L2 (APOL2). The apolipoprotein family of proteins facilitates the tightly regulated delivery of lipids and lipophilic substrates to specific cells in the brain, as well as regulating signal transduction pathways [Reference Li, Cam and Bu29]. APOL2, mainly localized at the endoplasmic reticulum, is implicated in cholesterol biosynthesis and trafficking and is thought to mediate cell death induced by interferon-gamma or viral infection, indicating a role in inflammatory processes [Reference Galindo-Moreno, Iurlaro, El Mjiyad, Diez-Perez, Gabaldon and Munoz-Pinedo30]. Dysregulation of the inflammatory response system has been associated with the pathophysiology of schizophrenia [Reference Muller, Riedel, Ackenheil and Schwarz31]. Prior works have found the gene APOL2 was upregulated in the brains of schizophrenic patients, and polymorphisms in this gene were linked to schizophrenia risk [Reference Takahashi, Cui, Han, Fagerness, Galloway and Shen32]. Differential expression of APOL2 has also been observed in individuals with substance use disorders, including cocaine, cannabis, and phencyclidine [Reference Lehrmann33]. APOL2 is highly expressed in some brain regions, including the hippocampus, intralobular white matter, and the medulla [Reference Luo, Jung, Longley, Rosoff, Charlet and Muench34]. However, the biological function of APOL2 in the brain remains unclear.
Three gene variants identified by colocalization analysis implied a potential association with TRS pathogenesis, including CCDC91, PRDX1, and WARS2. CCDC91 is highly expressed in the central nervous system, located in the nucleoplasm and trans-Golgi networks and is involved in subcellular transport and protein localization in the Golgi complex. One study showed an increased incidence of copy number variation in CCDC91 in bipolar disorder patients [Reference Chen, Li, Wang, Liu, Zhang and Zhao35]. Its interacting partners, GGA1 and GGA2, have been implicated in the pathophysiology of AD through interactions with β-amyloid precursors [Reference Herskowitz, Offe, Deshpande, Kahn, Levey and Lah36]. Also, our work identified one GO annotation associated with CCDC91: identical protein binding. PRDX is a protein family with antioxidant enzyme activity that reduces hydrogen peroxide and alkyl hydroperoxides in cells. Numerous studies have demonstrated the anti-inflammatory and anti-apoptotic effects of PRDX1 [Reference Yang, Huang, Yuan, Zhang, Gong and Chen37]. It appears to have neuroprotective activities in neuronal cells, which reduce reactive oxygen species-mediated cell death in schizophrenia [Reference Lee38]. Antipsychotic drugs affect PRDX1 expression. PRDX1 expression was increased in haloperidol-treated C6 cells but decreased in C6 cells treated with risperidone and clozapine [Reference Chen39]. WARS2 encodes mitochondrial tryptophan-tRNA synthetase, a homologous class Ic enzyme. The clinical spectrum associated with WARS2 defects appears to be quite broad, including clinical features (cardiomyopathy, movement disorders, retinitis pigmentosa, optic atrophy, hypoglycemia, etc.) as well as the age of onset and clinical course [Reference Wortmann, Timal, Venselaar, Wintjes, Kopajtich and Feichtinger40]. Also, WARS2 mutations cause dopa-responsive early-onset parkinsonism and progressive myoclonus ataxia [Reference Skorvanek, Rektorova, Mandemakers, Wagner, Steinfeld and Orec41]. These works support the contention that the biparental loss-of-function WARS2 variants cause mitochondrial dysfunction and disease. The effect of variants in these three genes on TRS needs further study.
Functional enrichment and annotation analysis revealed several mitochondrial-related results implicating mitochondrial function in the pathogenesis of TRS. Some evidence has confirmed that mitochondrial dysfunction is an important pathological factor in schizophrenia, including decreased mitochondrial respiration due to altered complex I activity [Reference Ben-Shachar42], motor deficits, altered mitochondrial dynamics [Reference Flippo and Strack43], increased levels of mitochondrial DNA mutations [Reference Gonçalves, Cappi, Hagen, Sequeira, Vawter and Derkach44], and decreased cognitive abilities in mitochondrial diseases [Reference Geng, Liu, Xiong, Ding, Jiang and Yang45]. One study identified a strong correlation between a TRS susceptibility gene and mitochondrial dysfunction, which correlates with the dysregulation of NRG-1/mTOR/miR143-3p signaling [Reference Sun, Zhang, Cong, Chen, Yi and Huang46]. More attention should be paid to the role of mitochondria in TRS.
Our results revolve more around lipid oxidation and inflammation, which tends to be, but is not fully explained by the “inflammation and oxidative stress” hypothesis of the neurobiological mechanisms of TRS. More evidence is being mined to support this hypothesis. A study showed elevated lipid peroxidation in patients with TRS compared to treatment-responsive schizophrenia patients and healthy controls. This exacerbated peroxidation process in TRS may reflect a deeper abnormality in the fatty acid content of synaptic membranes, leading to the dysfunction of neurons as well as its microenvironment [Reference Medina-Hernandez, Ramos-Loyo, Luquin, Sanchez, Garcia-Estrada and Navarro-Ruiz47]. Early neuroinflammation and chronic hyperactivation are thought to contribute to schizophrenia; high levels of inflammation may also play a role in treatment resistance [Reference Howes and McCutcheon48]. There is no universally accepted or defined mechanism for TRS, limited by different criteria for TRS or small study sample sizes. Other hypotheses have been proposed regarding the neurobiological mechanisms of TRS, including differences in dopaminergic function, glutamate dysregulation, and serotonin dysregulation. Previous studies have provided genetic evidence for different hypotheses, indicating that TRS may develop through multiple pathways or change in potential mechanisms at specific times [Reference Ruderfer, Charney, Readhead, Kidd, Kähler and Kenny49]. These theories are not mutually exclusive, and combining several pathways may contribute to the neurobiology of TRS [Reference Howes and Kapur50]. These results are insufficient to highlight the distinctiveness of TRS from schizophrenia and more characteristic markers are needed to understand the BPs and prove the heterogeneity of TRS.
Our work is the first TRS-related study with proteomic analysis based on a large sample of GWAS summary data, providing an accurate functional profile that presents an unbiased picture of current physiological status. Despite this, the study has certain limitations. TRS is an underreported diagnosis, and although our definition of the phenotype is in line with international criteria, we acknowledge that there may still be individuals with treatment-resistant symptoms in non-TRS datasets. This could result in imperfect phenotypes and misclassifications that weaken our findings and reduce the exploratory power of the brain proteins analyzed. The study sample is mainly of European descent and the conclusions may not pertain to non-European countries. Furthermore, we performed a cross-sectional proteomics study, so no exploration of longitudinal changes in biomarkers was available. More research is needed to test our results in different patient populations and different phases of illness. Additionally, it is worth noting that the age and brain tissue preparation technique of the individuals included in the proteome reference datasets may have a slight impact on the final protein expression results. Finally, some proteins identified in the discovery of PWAS failed to be replicated in the confirmation PWAS, which is attributed to the limited sample size and the stochastic nature of high-throughput proteomic sequencing.
In conclusion, we performed PWAS analysis and identified two TRS-associated brain proteins, CPT2 and APOL2. Colocalization analysis based on PWAS results identified three variants that were causally related to protein expression, including CCDC91, PRDX1, and WARS2. Our results cautiously support that the pathological mechanism of TRS is linked to lipid oxidation and inflammation, where mitochondria-related functions may play a role.
Supplementary material
The supplementary material for this article can be found at http://doi.org/10.1192/j.eurpsy.2023.20.
Data availability statement
The datasets can be downloaded from the Psychiatric Genomics Consortium website (http://pgc.unc.edu).
Author contribution
Material preparation, data collection, and analysis were performed by Wenming Wei and Huijie Zhang. The first draft of the manuscript was written by Wenming Wei. The figures and tables were made by Bolun Cheng, Xiaoyue Qin, and Dan He. The literature searches were performed by Na Zhang, Yijing Zhao, Qingqing Cai, Xiaoge Chu, Sirong Shi, Liu Huan, Yan Wen, and Yumeng Jia. The study design was performed by Feng Zhang.
Financial support
This work was supported by the National Natural Scientific Foundation of China [81922059] and the Natural Science Basic Research Plan in Shaanxi Province of China [2021JCW-08].
Consent to participate
Written informed consent was obtained from each participant of previously published GWASs before data collection.
Consent for publication
All the authors have read and approved the final version of the manuscript.
Competing interest
The authors report no financial interests or potential conflicts of interest.
Ethical standard
This study is based on publicly available summarized data. The protocol and data collection were approved by the ethics committee of each genome-wide association study.
Open access
Comments
No Comments have been published for this article.