
==== Front
BMC Med Genomics
BMC Med Genomics
BMC Medical Genomics
1755-8794
BioMed Central London

2003
10.1186/s12920-024-02003-5
Research
Identification of autophagy-related genes as potential biomarkers correlated with immune infiltration in bipolar disorder: a bioinformatics analysis
Cao Dong 1234
Liu Yafang 14
Mei Jinghong 1
Yu Shuailong 1
Zeng Cong 1
Zhang Jing zhjing63@mail.sysu.edu.cn

14
Li Yujuan liyuj@mail.sysu.edu.cn

124
1 grid.12981.33 0000 0001 2360 039X Department of Anesthesiology, Brain Research Center, Guangdong Province Key Laboratory of Malignant Tumor Epigenetics and Gene Regulation, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, 510120 China
2 https://ror.org/0064kty71 grid.12981.33 0000 0001 2360 039X Guangdong Province Key Laboratory of Brain Function and Disease, Zhongshan School of Medicine, Sun Yat-sen University, Guangzhou, 510080 China
3 grid.12981.33 0000 0001 2360 039X Department of Anesthesiology, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, guangzhou, China
4 grid.12981.33 0000 0001 2360 039X Department of Anesthesiology, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, No.107 Yanjiang West Road, Guangzhou, 510120 China
13 9 2024
13 9 2024
2024
17 23130 10 2023
2 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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-nc-nd/4.0/.
Background

Bipolar disorder (BPD) is a kind of manic and depressive phase alternate episodes of serious mental illness, and it is correlated with well-documented cortical brain abnormalities. Emerging evidence supports that autophagy dysfunction in neuronal system contributes to pathophysiological changes in neurological disease. However, the role of autophagy in bipolar disorder has rarely been elucidated. This study aimed to identify the autophagy-related gene as a potential biomarker Correlated to immune infiltration in BPD.

Methods

The microarray dataset GSE23848 and autophagy-related genes (ARGs) were downloaded. Differentially expressed genes (DEGs) between normal and BPD samples were screened using the R software. Machine learning algorithms were performed to screen the significant candidate biomarker from autophagy-related differentially expressed genes (ARDEGs). The correlation between the screened ARDEGs and infiltrating immune cells was explored through correlation analysis.

Results

In this study, the autophagy pathway was abundantly enriched and activated in BPD, as indicated by Pathway enrichment analysis. We identified 16 ARDEGs in BPD compared to the normal group. A signature of 4 ARDEGs (ERN1, ATG3, CTSB, and EIF2AK3) was screened. ROC analysis showed that the above genes have good diagnostic performance. In addition, immune correlation analysis considered that the above four genes significantly correlated with immune cells in BPD.

Conclusions

Autophagy - immune cell axis mediates pathophysiological changes in BPD. Four important ARDEGs are prospective to be potential biomarkers associated with immune infiltration in BPD and helpful for the prediction or diagnosis of BPD.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-024-02003-5.

Keywords

Autophagy
Bipolar disorder
Immune
Machine learning algorithms
Macrophage
Monocyte
http://dx.doi.org/10.13039/501100014857 National Natural Science Foundation of China-Guangdong Joint Fund 2018A0303130272 Science and Technology Planning Project of Guangzhou, China20210201020 issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcIntroduction

Bipolar disorder (BPD) is a severe mental illness in which a person’s mood swings back and forth between mania and depression. Existing research shows that BPD affects more than 1% of the global population and is one of the major causes of psychological change in young people [1]. The latest data from the World Mental Health Surveys involving 156,331 respondents across 29 countries from 2001 to 2022 reveal that the probability of the first episode peaks around the age of 15, with a median onset age of approximately 20 for bipolar disorder [2]. In addition, patients with bipolar disorder show variable courses and often leading to cognitive impairment and significantly reduced quality of life [3, 4], and tend to get worse as the disease progresses [5, 6], and in more severe cases, individuals with bipolar disorder may face life-threatening outcomes, with a heightened association with incidents involving harm, such as suicide, interpersonal violence, and road injuries, as well as cardiovascular diseases, notably ischemic heart disease [7].

The Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition, Text Revision outlines stringent criteria for the diagnosis of bipolar disorder, necessitating the presence of at least one episode of mania or hypomania [8]. This nuanced clinical framework underscores the intricacies associated with the characterization and identification of bipolar pathology. The challenges faced by clinicians encompass discerning subtle variations in mood episodes, disentangling bipolar manifestations from related mood disorders, and negotiating the inherent subjectivity in symptom reporting. Despite the meticulous delineation of these clinical parameters, the molecular diagnostic landscape for individuals with BPD remains characterized by uncertainty.

In some studies of recurrent affective bipolar disorder, progressive brain structure and function changes have been observed in patients [9]. Recurrent episodes of long-term illness have been linked to reduced thickness in parts of the cerebral cortex, such as the prefrontal cortex, which may play an essential role in stress regulation [10]. Epigenetic mechanisms [11], mitochondrial dysfunction, pathways that maintain neuroplasticity, increased inflammation, and oxidative stress are all thought to contribute to neural progression in bipolar disorder [9].

An increasing body of evidence suggests that BPD is correlated with the physiological processes mediated by autophagy [12–14]. Autophagy is primarily involved in protein degradation pathways and plays a vital role in helping the body remove protein aggregates and misfolded proteins from healthy cells [15]. According to the present researches, changes in autophagy-related factors in the prefrontal cortex are involved in regulating depression-like behaviors [16].The Akt-mTOR pathway, which plays a crucial role in the regulation of autophagy, is reduced in the prefrontal cortex in BPD patients, and decreased activity leads to cognitive impairment associated with changes in synaptic connectivity and function [12]. Moreover, imaging, biochemical and genetic research has suggested that mitochondrial dysfunction is the core of BPD features [17]. Previous studies have demonstrated a complex interplay between autophagy and mitochondria, particularly in stress responses such as oxidative stress [18] and endoplasmic reticulum stress [19, 20]. In the context of bipolar disorder, dysregulation of these stress responses is implicated in the pathophysiology of the condition [21]. The intricate signaling pathways activated by stressors contribute to the dysregulation of autophagy and mitochondrial dynamics, potentially playing a crucial role in the progression of bipolar disorder.

Numerous studies have noted that immune disorders accompany patients with BPD (e.g., increased production of pro-inflammatory cytokines, activation of monocytes and macrophages, reduced T cell numbers or activity) [22–26]. Moreover, autophagy interacts with immune system. Autophagy can control the immune response by regulating the activity of immune cells and the production of cytokines [27]. In turn, autophagy is also significantly affected by a wide variety of immune cells and cytokines [28]. Consequently, it becomes crucially critical to examine the relationship between immune cells and autophagy-related genes (ARGs) in BPD.

Nevertheless, understanding the pathophysiology and molecular mechanisms leading to BPD is still limited, and more unambiguous evidence for the correlation between BPD and autophagy should be provided further. Early diagnosis is quite difficult for professionals who did not previously comprehend the longitudinal evolution of the disease because there are no exact biomarkers available.

This research sought to link the ARG signature to likely diagnosis of BPD. Furthermore, the immune status of BPD patients was further explored to reveal the potential immune infiltration pattern of ARGs in BPD, which may be helpful for diagnosing and treating BPD.

Materials and methods

Gathering and processing of gene expression omnibus (GEO) datasets

The GEO database was used in this investigation to explore the gene expression profile of BPD [29]. Data from 20 BPD group and 15 control group samples’ peripheral blood gene expression matrices were included in the GSE23848 dataset, and the GPL6106 platform was used to evaluate the aforementioned information. The composition of BPD samples includes both patients undergoing medication treatment and those without medication treatment. Considering the diversity within the actual patient population, we did not conduct separate analyses, aiming to provide a more comprehensive reflection of the characteristics of bipolar disorder. All data employed in the study originated from GEO, such that ethical approval and informed consent were not required. The workflow chart is illustrated in Fig. 1.

Fig. 1 Workflow of the research. Abbreviations are defined as follows, Gene Expression Omnibus database (GEO), differentially expressed genes (DEGs), Autophagy-related differentially expressed genes (ARDEGs), protein-protein interaction (PPI), Principal Component Analysis (PCA), Receiver Operating Characteristic (ROC)

Selection and functional enrichment of DEGs in BPD patients

The GEO2R online analytic tool was used to find the DEGs between the BPD and normal groups [29]. In the identification of DEGs associated with BPD, our approach focused on an adjusted P-value threshold of less than 0.05. Unlike conventional methodologies, we opted not to incorporate a specific Log2 fold change threshold in our criteria. This decision was grounded in the consideration that BPD, as a complex and multifaceted disorder, may involve subtle yet biologically significant changes in gene expression. By exclusively utilizing the adjusted P-value criterion, we aimed to capture a comprehensive view of gene expression alterations without imposing restrictions on the magnitude of fold changes. This approach aligns with the intricate nature of psychiatric disorders, where nuanced expression changes can hold biological relevance. Based on the “clusterProfiler” package in the R programming language, enrichment studies for the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were carried out.

Acquisition of autophagy-related differential genes (ARDEGs)

We were able to extract 221 genes associated with autophagy from the human autophagy database (HADb, http://www.autophagy.lu/). The “venn” package in R was used to find the genes that are differentially expressed in autophagy. Principal Component Analysis (PCA) was employed to provide a comprehensive exploration and analysis of the identified genes.

Screening of potential biomarkers for diagnosis of bipolar disorder based on machine learning algorithm and construction of PPI network

A support vector machine (SVM), random forest (RF), and generalized linear model (GLM) were established depending on GSE23848. The response variable was BPD diagnosis, while the explanatory variable was ARDEGs. The three models were examined using the interpretative capabilities of the “DALEX” package in R, and the residual distribution was improved using the test set. This research additionally examines the significance of these factors, screens the core genes employing the PPI network, and derives four significant explanatory variables. Finally, the “glmnet” package in R was used to build the Least Absolute Shrinkage and Selection Operator (LASSO) analysis with penalty settings for 10-fold cross-validation, which was utilized to validate the screened genes.

Evaluation of expression levels and diagnostic significance of crucial genes identified

The expression levels of the crucial genes were examined using the Wilcoxon rank-sum test. Additionally, a receiver operating characteristic curve (ROC) analysis was conducted to assess the diagnostic performance of these genes in differentiating between samples with BPD and healthy controls.

Immune status assessment

The ESTIMATE algorithm, employed in this study, was utilized to compute immune scores and stromal scores in the GSE23848 dataset for DEGs between the control and BPD groups. ESTIMATE (Estimation of Stromal and Immune cells in samples using Expression data) is a computational approach designed to infer the presence of stromal and immune cells within samples based on gene expression profiles.

Immune infiltration

A total of 27 immune-related cells and related genes were collected (for details, see Supplementary material_1) [30]. The R package “GSVA” was used to analyze the differential immune cell infiltration in the GSE23848 dataset; “ggpubr” was used for visualization; and “corrplot” was used to analyze the relationship between four crucially significant genes and immune cells. The Pearson test was performed to calculate the correlation coefficient, and the significance test.

Prediction of transcription factors

Using ChIP-X Enrichment Analysis (ChEA3), the transcription factors associated with ARDEGs were predicted [31]. A tool for TFs enrichment analysis and evaluating TFs of user-submitted gene sets, the ChEA3 database contains several gene banks from various sources.

Statistical analysis

Statistical analysis was conducted using the R programming language, specifically version 4.3.1. A significance threshold of p < 0.05 was employed, indicating statistical significance.

Result

Differentially expressed genes and their functional enrichment in peripheral blood of BPD patients

In order to investigate if autophagy plays a role in the development of BPD, we analyzed the DEGs and their functional enrichment in peripheral blood in BPD patients compared to healthy individuals. 400 up-regulated genes and 599 down-regulated genes, totaling 999 DEGs, were found in GSE23848 (Fig. 2A). The expression profile of all DEGs was shown in the heatmap (Fig. 2B). Furthermore, the 999 DEGs were subjected to GO and KEGG functional enrichment analysis. The genes mentioned above, located inside of cytosol, which mostly display catalytic activity and nucleic acid binding, are positive regulators of the cellular protein metabolism process (Fig. 2C). According to the KEGG enrichment analysis (Fig. 2D), the highest-level classical pathways correlated with differentially expressed genes (DEGs) included autophagy, apoptosis, and Alzheimer’s disease modulation. These findings suggest a close association between bipolar disorder (BPD) and autophagy, with a potential pivotal role of the immune system in BPD.

Fig. 2 Differentially expression analysis and Gene Ontology and KEGG pathway enrichment analysis. (A) Volcano plot of significant DEGs between BPD and Ctrl samples (BPD: Bipolar disorder; Ctrl: Control). (B) A heatmap of the DEGs. (C) The results of GO are presented in cnetplot. (D) The results of KEGG are presented in Lollipop illustration

Acquisition of autophagy-related differential genes (ARDEGs)

The intersection of 999 DEGs and 221 Autophagy relative genes were further taken to screen out autophagy relative genes from DEGs. On that basis, 16 differential ARGs were obtained (Fig. 3A). Subsequently, ARDEGs were classified through PCA analysis. The results of the principal component analysis revealed that 16 genes exhibit a clear distinction between BPD and control group samples. This observation suggests the potential significance of these genes in the diagnosis of BPD; however, it should be noted that experimental confirmation is required to validate their diagnostic relevance (Fig. 3B). Figure 3C and D present the expression profiles of these 16 ARDEGs in different samples, displayed in the form of a heatmap and violin plot, respectively.

Fig. 3 Identification of ARDEGs. (A) Venn diagram was generated to obtain the intersection of DEGs in GSE23848 dataset and ARGs originating from the Human Autophagy Database covered 16 ARDEGs. (B) PCA analysis was performed to classify infiltrating ARDEGs between BPD and normal peripheral blood samples. The expression levels of ARDEGs are presented in the heatmap (C) and vioplot (D). *p < 0.05, **p < 0.01, ***p < 0.001

Selection of signature genes via support vector machine model algorithms, LASSO analysis and PPI network

The 16 ARDEGs above were employed as crucial genes to build three models to further narrow the scope of ARDEGs. As illustrated in Fig. 4(A-B), the SVM model is considered the optimal match for minimizing sample residues. Five significant explanatory variables (i.e., KLHL24, ERN1, ATG3, CTSB, and EIF2AK3) were selected from the SVM model (Fig. 4C). Next, the functional interaction of 16 ARDEGs was studied using STRING 11.5 (http://string-db.org). The KLHL24 gene without protein interaction was deleted, and four essential genes were obtained (Fig. 4D). Finally, on the basis of 16 ARDEGs expression matrices, the regression analysis model of LASSO was established, and the LASSO model was optimized and selected with the optimal λ value by using the minimum standard of 10-fold cross-validation (Fig. 4E). It was decided to employ a value of = 0.0115532 with log (λ)= -4.460793. The log (λ) sequence was plotted against a coefficient profile (Fig. 4F). ERN1, ATG3, CTSB, and EIF2AK3 coefficients were 0.37669785, 0.33988395, 0.12711732, and − 0.22049001, respectively.

Fig. 4 Construction and assessment of RF, GLM and SVM model. (A) Cumulative residual distribution map of the sample. (B) Boxplots of the residuals of the sample. Red dot stands for root mean square of residuals. (C) The significance of the variables in the SVM model. (D) Functional protein association network of ARDEGs constructed using STRING dataset. (E) LASSO coefficient profiles. (F) Ten time cross-validation for tuning parameter selection in the LASSO model

Evaluation of the selected genes’ diagnostic efficacy and expression levels in BPD

Subsequently, we performed further analysis on the four significant variables that the SVM model had filtered. The chromosomal locations of ERN1, ATG3, CTSB, and EIF2AK3 were shown in Fig. 5B. Patients with BPD reported higher peripheral blood levels of ERN1, ATG3, and CTSB than patients without the condition. However, compared to control group patients, EIF2AK3 expression in peripheral blood was lower in BPD patients (Fig. 5A and C). The following four markers can effectively discriminate between BPD and control group patients, as demonstrated in the PCA analyses in Fig. 5D, demonstrating that they are essential for diagnosing BPD. In this regard, ROC analysis verified the diagnostic validity of the above four genes for BPD. AUC greater than 0.800 could diagnose BPD with excellent specificity and sensitivity. The AUC values of ERN1, ATG3, CTSB, and EIF2AK3 were 0.827 (95%CI 0.933-0.600), 0.867 (95%CI 0.933 − 0.750), 0.862 (95%CI 0.867 − 0.850) and 0.882 (95%CI 0.733–0.950), respectively, As shown in Fig. 5E-H. Furthermore, a multi-gene combined diagnosis model was built, and the result indicated that the AUC of those models are all greater than 0.9 and the 4-gene model was 1 through logistic regression analysis (Fig. 5I).

Fig. 5 Relative expression level of the 4 ARDEGs. (A) Heat map of the expression pattern of the 4 ARDEGs. (B) The chromosomal locations of the 4 ARDEGs (C) The relative expression level of the 4 ARDEGs between BPD and control group from GSE23848 dataset. (D) Principal component analysis shows that the four genes aforementioned can clearly distinguished BPD and control group. (E-H) The GSE23848 dataset was used to validate the diagnostic effectiveness of the four important ARDEGs by ROC analysis. (I) The diagnostic performance of multi-gene combined model was calculated based on the 4 genes’ expression. *p < 0.05, **p < 0.01, ***p < 0.001

Distribution of Immune Infiltration Pattern in BPD patients

Existing research has suggested that autophagy is significantly correlated with the immune system [27, 32–35]. To gain insights into the interaction between immunity and autophagy in BPD, we evaluated the immune status of DEGs in patients with BPD compared with healthy people based on the ESTIMATE algorithm. As shown in Fig. 6, BPD patients had a lower immune score (Fig. 6A. Wilcoxon test, p = 9.8E-07) and higher stromal scores (Fig. 6B. Wilcoxon test, p = 3.1E-07), which indicated immune dysfunction might participate in the occurrence and development of BDP patients.

Fig. 6 Assessment of immune status and immune cell infiltration patterns in BPD samples and control group samples. The immune score (A) and the stromal score (B) of DEGs in patients with bipolar disorder compared to healthy people based on the ESTIMATE algorithm. (C) Heatmap of 21 immune cell subpopulations in each BPD and control group sample. (D) Correlation heatmap of 21 immune cells. (E) Boxplot showing the differentially infiltrated immune cells between the two groups. (F) Correlation among four important ARDEGs and 21 differential immune cells. *p < 0.05, **p < 0.01, ***p < 0.001

Subsequently, we utilized the R package ‘GSVA’ to assess the relative proportions of the 27 immune cell biomarkers in each sample, aiming to gain a more comprehensive understanding of the immune microenvironment’s functionality in BPD. The analysis revealed statistically significant associations for 21 distinct types of immune cells. Using a heatmap and a boxplot, we could compare 21 immune cell infiltrations across samples BPD and normal. We found a variety of immune cell infiltration abundances that were statistically significant between BPD and control group samples and a correlation between the above differentially associated immune cells. In this regard, we found that the four essential ARDEGs above were correlated with differentially related immune cells, such as ERN1 showed a positive correlation with Macrophage (R = 0.38, P = 7.83E-07) and Monocyte (R = 0.39, P = 3.68E-07), ATG3 showed a positive correlation with Macrophage (R = 0.51, P = 1.43E-07) and Monocyte (R = 0.56, P = 1.16E-09), and CTSB showed a positive correlation with Macrophage (R = 0.70, P = 5.82E-13) and Monocyte (R = 0.69, P = 4.62E-12). In addition, EIF2AK3 showed a negative correlation with Macrophage (R=-0.70, P = 3.36E-15) and Monocyte (R=-0.65, P = 1.02E-13). (Fig. 6F)

Biological function analysis and transcription factor regulation mechanism of the above four genes

We conducted GO and KEGG functional enrichment analyses for these four putative biomarkers in order to learn more about the biological processes and regulatory mechanisms of ERN1, ATG3, CTSB, and EIF2AK3. The result indicated that the above markers were not only correlated with autophagy, but they were likely to regulate pathophysiological processes (e.g., apoptosis (adjust. P < 0.001) and Alzheimer’s disease (adjust. P < 0.01)). (Fig. 7A-D)

Fig. 7 Functional correlation analysis, and Prediction of transcription factors for ARDEGs. (A, B) The results of KEGG were presented by bubble and circle charts. (C, D) The results of GO were shown by bubble and circle graphs. (E) The TF enrichment analysis of four important ARDEGs. Blue diamond indicate ARDEGs, and red circles indicate predicted TFs (Supplementary material_2). (F-I) Top10 TFs network of four ARDEGs

Furthermore, the transcription factors regulating the above four genes were predicted through the CHA3E website to study the regulatory mechanisms of ERN1, ATG3, CTSB, and EIF2AK3, and 1632 transcription factors were obtained in total. A total of 377 transcription factors were obtained by removing blank values (Fig. 7E, Please refer to Supplementary material_2 for detailed content). The above four genes were compared with the TF target gene set ChEA3 library assembled from multiple orthogonal omics datasets. The Fisher precision test (background size 20000) was used to compare the 4 screened genes with the TF target gene set to determine which TF might be most closely related to the 4 genes and to screen out the top 10 core transcription factors most related to each gene. (Fig. 7F-I. Table 1)

Table 1 The top 10 transcription factors of prediction results for the 4 ARDEGs

ERN1/Rank	TFs	Mean Rank	CTSB/Rank	TFs	Mean Rank	
1	BHLHA15	19	1	NFE2L2	71.4	
2	IRF4	62	2	ZNF581	71.67	
3	MAFB	113.3	3	GLMP	105	
4	CEBPA	129.4	4	MAFB	113.3	
5	CEBPB	141	5	MITF	114.8	
6	PPARD	144.3	6	CEBPA	129.4	
7	ZNF710	159	7	PAX8	132.3	
8	ZNF217	167.8	8	ATF5	137	
9	XBP1	175.3	9	CEBPB	141	
10	ELF1	178.2	10	PPARD	144.3	
ATG3/Rank			EIF2AK3/Rank			
1	C11ORF95	31	1	BHLHA15	19	
2	IRF4	62	2	IRF4	62	
3	GZF1	68.33	3	GZF1	68.33	
4	DUX4	109.5	4	MITF	114.8	
5	MITF	114.8	5	ZBTB32	140.3	
6	CEBPA	129.4	6	CEBPB	141	
7	CEBPB	141	7	IKZF2	147	
8	KDM5B	192.4	8	ZNF710	159	
9	USF1	202	9	SPI1	165	
10	NFE2	213.4	10	ZBTB42	165.5	
This provides a ranking of transcription factors whose putative transcriptional targets are most closely similar to the query set. Integrated results, which take into account results from all libraries, are sorted in ascending order by score. Lower scores indicate more relevancy to the transcription factor

Discussion

In the clinical setting, the precise and timely diagnosis of BPD continues to pose a formidable challenge, regardless of the overt manifestation of symptoms. Therefore, the diagnosis and treatment of BPD still need to be vigorously promoted. This work identifies useful diagnostic biomarkers for BPD and sheds light on the landscape of autophagy associated with the disorder. Additionally, this study examines the relationship between immune cells and ARDEGs, which are involved in the onset and progression of BPD. In a nutshell, this study used bioinformatics analysis to first pinpoint four ARDEGs as possible biomarkers associated with immune infiltration in BPD patients.

We integrated the differentially expressed genes of patients with bipolar disorder and the gene set of the Human Autophagy Database. Based on machine learning methods, we screened four important ARDEGs (ERN1, ATG3, CTSB, and EIF2AK3). An essential sensor of unfolded proteins in the ER is the signaling enzyme endoplasmic reticulum nuclear signal 1 (ERN1) [36]. According to a recent study, exosomal microRNA-124-3p from bone marrow mesenchymal stem cells may reduce ERN1 levels and thereby lessen nerve damage in spinal cord ischemia reperfusion injury [37]. Unfolded protein response (UPR) results in neuronal apoptosis by regulating ERN1/IRE1 or EIF2AK3/PERK pathways [38]. In mental diseases such schizophrenia, depression and post-traumatic stress disorder, aberrant neuronal death is mediated by endoplasmic reticulum stress [39]. The endoplasmic reticulum stress pathway involves ERN1, which has a strong correlation with autophagy [40–42]. Our results indicated that ERN1 was up-regulated in patients with BPD compared with normal, suggesting a complex interplay between endoplasmic reticulum function, autophagy, and BPD. Autophagy-related gene 3 (ATG3) refers to an enzyme primarily known for its role in the lipidation of LC3B, and it is of critical significance to autophagy. More recently, dysfunctional mitochondria in early BPD have been identified by co-localizing mitochondria (Hsp60) with autophagosomes (LC3B) in BPD patient cells [43]. This finding suggested that ATG3 was implicated in the development of BPD, and this research verified that the expression level of ATG3 in blood samples of BPD patients was greater than that of the control group. A member of the cysteine cathepsin family, CTSB (cathepsin B), contributes to neuropathological alterations in Alzheimer’s disease and traumatic brain injury in addition to being strongly linked to autophagy [44–46]. The hippocampus and cortical areas of healthy individuals’ brains were discovered to have significant levels of CTSB expression [47]. As revealed by the transcriptome analysis of different brain regions in the mouse limbic system, CTSB is closely correlated with mood [48]. Our investigation discovered that BPD patients had much greater levels of CTSB expression in their peripheral blood than the control group, indicating that CTSB may have a special function in BPD. EIF2AK3(Eukaryotic translation initiation factor 2α kinase 3) may be involved in the pathogenesis of BPD by regulating ER stress. The down-regulation of PERK/EIF2AK3 leads to increase of endoplasmic reticulum stress and impaired apoptosis induction, antioxidant response and autophagy impaired flux [49]. Our study found that the expression of EIF2AK3 in the blood samples of BPD patients was lower than that of the control group. In brief, it is speculated that the four ARDEGs may become the peripheral blood monitoring markers for BPD. In order to confirm our suspicions, ROC curves were created to assess the four aforementioned genes’ diagnostic effectiveness. In this study, the AUCs of the above ARDEGs were all greater than 0.8, and the combination of multiple genes achieved a higher diagnostic efficiency.

There is growing proof that immune dysfunction and inflammation play a significant role in the development of bipolar disorder [9, 50]. Actually, a number of autoimmune processes can have an impact on the brain, and this may explain various psychiatric problems. As evidenced by pathogenic microglial over-activation and elevated levels of pro-inflammatory cytokines, manic and depressive episodes are associated with the activation of neuroinflammation mechanisms [51–53]. 27 immune-associated cells from the BPD and healthy groups were examined for immune cell infiltration in this investigation. The expression of activated CD8 T cell, effector memory CD4 T cell, gamma delta T cell, immature B cell, macrophage, monocyte, macrophage, and Type 1 T helper cells in BPD was significantly up-regulated. However, the expression of myeloid-derived suppressor cell (MDSC), activated CD4 T cell (activated dendritic cell), memory B cell, mast cell, plasmacytoid dendritic cell, regulatory T cell, T follicular helper cell was significantly down-regulated. In existing research, significant changes in the circulating frequency of T cells and their subsets, playing a key role in cellular immune response, have been identified in patients with BPD [22–24]. The result indicated that patients with bipolar disorder showed a significant reduction in regulatory T cells and a significant bias towards high levels of Type 1 T helper cells rather than Type 2 T helper cells [25]. Drexhage, R.C. et al. verified that BPD patients had considerably higher levels of proinflammatory monocyte expression [54]. High levels of heterogeneity and plasticity allow macrophages to develop into M1 and M2 macrophages under various circumstances. Macrophages mostly differentiate into the pro-inflammatory M1 type in the peripheral blood microenvironment of BPD [26]. Additionally, dendritic cells have been linked to the genesis of BPD. They exhibit abnormalities in BPD and totally recover or even activate after receiving lithium in vivo [55, 56]. In light of the fact that our study’s findings agree with other studies, bioinformatics analysis is used to highlight the importance of the aforementioned immune cells in the pathogenesis of BPD.

Based on the significance of autophagy and immunological infiltration in BPD, the link between four efficient autophagy-related biomarkers (i.e., ERN1, ATG3, CTSB, and EIF2AK3) and various immune cells in BPD was further investigated. The results for Macrophage and Monocyte, which have a strong correlation with the four ARDEGs, are the most interesting. Macrophage and monocyte demonstrated positive correlations with ERN1, ATG3, and CTSB. Conversely, Macrophage and Monocyte had a negative connection with EIF2AK3. Existing research has suggested that ERN1 and EIF2AK3 mediate oxidative stress and induce macrophage autophagy [57], and they have a specific function in the control of monocyte alterations brought on by inflammation [58]. CTSB is mainly expressed in macrophage and promotes collagen synthesis in the infiltrated region of macrophage [59]. Thus, the notion that the autophagy-related gene-immune cell axis plays a crucial role in the early pathogenesis of BPD was put out. This theory still requires experimental confirmation.

Discussion section: limitations

In the pursuit of scientific inquiry, it is imperative to acknowledge and critically evaluate the limitations inherent in our study. By transparently addressing these limitations, we aim to provide a comprehensive understanding of the scope and potential impact of our findings.

One notable limitation of our study lies in the reliance on a single database, namely the Gene Expression Omnibus (GEO). While GEO offers valuable and relevant gene expression data related to bipolar disorder, the exclusivity of this database introduces a potential constraint. The lack of multiple database inclusion may limit the generalizability of our findings and could impact the robustness of the conclusions.

This study conducted a thorough exploration of the GEO database and considered datasets such as GSE46449, GSE62191, and GSE124326 for multifaceted validation of our experiments. However, discrepancies in sample composition and experimental design were identified, rendering these datasets unsuitable for meeting the specific objectives of our research. To elaborate, GSE46449 employed leukocyte samples, whereas our training dataset GSE23848 utilized peripheral blood tissue samples. This divergence may introduce distinct gene expression patterns, making GSE46449 inadequate for robustly validating our hypotheses. Similarly, GSE62191 utilized cells from the cerebral cortex, representing a significantly different tissue source compared to our study population, and was consequently excluded. Furthermore, the datasets from GSE124326 involved subjects under intense lithium treatment, introducing a potential confounding factor related to lithium therapy. To ensure the precision of our study in reflecting the distinct characteristics of bipolar disorder, we chose to exclude datasets with these specific limitations.

Furthermore, in this study, we included a total of 20 BPD patient samples and 15 control group patient samples. However, it is crucial to explicitly acknowledge that this relatively small sample size may introduce potential bias, impacting the generalizability of our research findings. Due to the limitation in sample size, the statistical power of our study may be constrained, preventing a comprehensive representation of the entire patient population. It is noteworthy that, particularly for diagnostic purposes, the small sample size could affect the accuracy of our understanding of the patient group. Future research endeavors should focus on enlarging the sample size, employing larger-scale study designs to ensure the reliability and generalizability of the results.

Conclusion

The autophagy-immune cell axis takes on a critical significance in the early pathogenesis of BPD. Four important ARDEGs (i.e., ERN1, ATG3, CTSB, and EIF2AK3) serve as potential biomarkers correlated with immune infiltration in BPD and be helpful for the prediction or diagnosis of BPD.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Acknowledgements

We acknowledge the GEO database for providing their platforms as well as their contributors for uploading meaningful datasets.

Author contributions

YJL and JZ contributed to the study concepts and study design and helped in revising the manuscript. YJL provided funding acquisition. DC and YFL supervised the research, analyzed the data, and wrote the manuscript. JHM, SLY and CZ performed data management, and bioinformatics analysis. All authors read and approved the final manuscript.

Funding

This work was supported by the Natural Science Foundation of Guangdong Province, China (No.2018A0303130272); the Science and Technology Planning Project of Guangzhou, China (20210201020).

Data availability

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

List of Abbreviations

BPD Bipolar disorder

GEO Gene Expression Omnibus

DEGs Differentially expressed genes

ARDEGs Autophagy-related differentially expressed genes

PPI protein-protein interaction

ROC receiver operating characteristic curve

GO Gene Ontology

KEGG Kyoto Encyclopedia of Genes and Genomes

LASSO Least absolute shrinkage and selection operator

ERN1 Endoplasmic reticulum nuclear signal 1

ATG3 Autophagy-related gene 3

CTSB cathepsin B

EIF2AK3 eukaryotic translation initiation factor 2 alpha kinase 3

PCA Principal Component Analysis

RF random forest

SVM support vector machine

GLM generalized linear model

TF Transcription fact

ChEA3 ChIP-X Enrichment Analysis 3

KLHL24 Kelch like family member 24

AUC Area Under Curve

UPR Unfolded protein response

ER Endoplasmic reticulum

LC3B Microtubule associated protein 1 light chain 3 beta

Hsp60 Heat shock protein family D (Hsp60) member 1

Publisher’s note

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

Dong Cao, Yafang Liu and Jinghong Mei contributed equally to this work.
==== Refs
References

1. Alonso J Petukhova M Vilagut G Chatterji S Heeringa S Üstün TB Days out of role due to common physical and mental conditions: results from the WHO World Mental Health surveys Mol Psychiatry 2011 16 12 1234 46 10.1038/mp.2010.101 20938433
Alonso J, Petukhova M, Vilagut G, Chatterji S, Heeringa S, Üstün TB, et al. Days out of role due to common physical and mental conditions: results from the WHO World Mental Health surveys. Mol Psychiatry. 2011;16(12):1234–46.20938433 10.1038/mp.2010.101
2. McGrath JJ Al-Hamzawi A Alonso J Altwaijri Y Andrade LH Bromet EJ Age of onset and cumulative risk of mental disorders: a cross-national analysis of population surveys from 29 countries Lancet Psychiatry 2023 10 9 668 81 10.1016/S2215-0366(23)00193-1 37531964
McGrath JJ, Al-Hamzawi A, Alonso J, Altwaijri Y, Andrade LH, Bromet EJ, et al. Age of onset and cumulative risk of mental disorders: a cross-national analysis of population surveys from 29 countries. Lancet Psychiatry. 2023;10(9):668–81.37531964 10.1016/S2215-0366(23)00193-1
3. Martinez-Aran A Vieta E Torrent C Sanchez-Moreno J Goikolea JM Salamero M Functional outcome in bipolar disorder: the role of clinical and cognitive factors Bipolar Disord 2007 9 1–2 103 13 10.1111/j.1399-5618.2007.00327.x 17391354
Martinez-Aran A, Vieta E, Torrent C, Sanchez-Moreno J, Goikolea JM, Salamero M, et al. Functional outcome in bipolar disorder: the role of clinical and cognitive factors. Bipolar Disord. 2007;9(1–2):103–13.17391354 10.1111/j.1399-5618.2007.00327.x
4. Grande I Goikolea JM de Dios C González-Pinto A Montes JM Saiz-Ruiz J Occupational disability in bipolar disorder: analysis of predictors of being on severe disablement benefit (PREBIS study data) Acta Psychiatr Scand 2013 127 5 403 11 10.1111/acps.12003 22924855
Grande I, Goikolea JM, de Dios C, González-Pinto A, Montes JM, Saiz-Ruiz J, et al. Occupational disability in bipolar disorder: analysis of predictors of being on severe disablement benefit (PREBIS study data). Acta Psychiatr Scand. 2013;127(5):403–11.22924855 10.1111/acps.12003
5. Rosa AR González-Ortega I González-Pinto A Echeburúa E Comes M Martínez-Àran A One-year psychosocial functioning in patients in the early vs. late stage of bipolar disorder Acta Psychiatr Scand 2012 125 4 335 41 10.1111/j.1600-0447.2011.01830.x 22283440
Rosa AR, González-Ortega I, González-Pinto A, Echeburúa E, Comes M, Martínez-Àran A, et al. One-year psychosocial functioning in patients in the early vs. late stage of bipolar disorder. Acta Psychiatr Scand. 2012;125(4):335–41.22283440 10.1111/j.1600-0447.2011.01830.x
6. Elshahawi HH Essawi H Rabie MA Mansour M Beshry ZA Mansour AN Cognitive functions among euthymic bipolar I patients after a single manic episode versus recurrent episodes J Affect Disord 2011 130 1–2 180 91 10.1016/j.jad.2010.10.027 21074274
Elshahawi HH, Essawi H, Rabie MA, Mansour M, Beshry ZA, Mansour AN. Cognitive functions among euthymic bipolar I patients after a single manic episode versus recurrent episodes. J Affect Disord. 2011;130(1–2):180–91.21074274 10.1016/j.jad.2010.10.027
7. Melo A, Dippenaar I, Johnson S, Weaver N, de Assis Acurcio F, Malta D, et al. All-cause and cause-specific mortality among people with severe mental illness in Brazil’s public health system, 2000-15: a retrospective study. The lancet Psychiatry; 2022.
8. Kaltenboeck A Winkler D Kasper S Bipolar and related disorders in DSM-5 and ICD-10 CNS Spectr 2016 21 4 318 23 10.1017/S1092852916000079 27378177
Kaltenboeck A, Winkler D, Kasper S. Bipolar and related disorders in DSM-5 and ICD-10. CNS Spectr. 2016;21(4):318–23.27378177 10.1017/S1092852916000079
9. Berk M Kapczinski F Andreazza AC Dean OM Giorlando F Maes M Pathways underlying neuroprogression in bipolar disorder: focus on inflammation, oxidative stress and neurotrophic factors Neurosci Biobehav Rev 2011 35 3 804 17 10.1016/j.neubiorev.2010.10.001 20934453
Berk M, Kapczinski F, Andreazza AC, Dean OM, Giorlando F, Maes M, et al. Pathways underlying neuroprogression in bipolar disorder: focus on inflammation, oxidative stress and neurotrophic factors. Neurosci Biobehav Rev. 2011;35(3):804–17.20934453 10.1016/j.neubiorev.2010.10.001
10. Hibar DP Westlye LT Doan NT Jahanshad N Cheung JW Ching CRK Cortical abnormalities in bipolar disorder: an MRI analysis of 6503 individuals from the ENIGMA bipolar disorder Working Group Mol Psychiatry 2018 23 4 932 42 10.1038/mp.2017.73 28461699
Hibar DP, Westlye LT, Doan NT, Jahanshad N, Cheung JW, Ching CRK, et al. Cortical abnormalities in bipolar disorder: an MRI analysis of 6503 individuals from the ENIGMA bipolar disorder Working Group. Mol Psychiatry. 2018;23(4):932–42.28461699 10.1038/mp.2017.73
11. Post RM Epigenetic basis of sensitization to stress, affective episodes, and stimulants: implications for illness progression and prevention Bipolar Disord 2016 18 4 315 24 10.1111/bdi.12401 27346321
Post RM. Epigenetic basis of sensitization to stress, affective episodes, and stimulants: implications for illness progression and prevention. Bipolar Disord. 2016;18(4):315–24.27346321 10.1111/bdi.12401
12. Vanderplow AM Eagle AL Kermath BA Bjornson KJ Robison AJ Cahill ME Akt-mTOR hypoactivity in bipolar disorder gives rise to cognitive impairments associated with altered neuronal structure and function Neuron 2021 109 9 1479 e966 10.1016/j.neuron.2021.03.008 33765445
Vanderplow AM, Eagle AL, Kermath BA, Bjornson KJ, Robison AJ, Cahill ME. Akt-mTOR hypoactivity in bipolar disorder gives rise to cognitive impairments associated with altered neuronal structure and function. Neuron. 2021;109(9):1479–e966.33765445 10.1016/j.neuron.2021.03.008
13. Scaini G Barichello T Fries GR Kennon EA Andrews T Nix BR TSPO upregulation in bipolar disorder and concomitant downregulation of mitophagic proteins and NLRP3 inflammasome activation Neuropsychopharmacology 2019 44 7 1291 9 10.1038/s41386-018-0293-4 30575805
Scaini G, Barichello T, Fries GR, Kennon EA, Andrews T, Nix BR, et al. TSPO upregulation in bipolar disorder and concomitant downregulation of mitophagic proteins and NLRP3 inflammasome activation. Neuropsychopharmacology. 2019;44(7):1291–9.30575805 10.1038/s41386-018-0293-4
14. Kim Y Santos R Gage FH Marchetto MC Molecular mechanisms of bipolar disorder: Progress made and Future challenges Front Cell Neurosci 2017 11 30 10.3389/fncel.2017.00030 28261061
Kim Y, Santos R, Gage FH, Marchetto MC. Molecular mechanisms of bipolar disorder: Progress made and Future challenges. Front Cell Neurosci. 2017;11:30.28261061 10.3389/fncel.2017.00030
15. Glick D Barth S Macleod KF Autophagy: cellular and molecular mechanisms J Pathol 2010 221 1 3 12 10.1002/path.2697 20225336
Glick D, Barth S, Macleod KF. Autophagy: cellular and molecular mechanisms. J Pathol. 2010;221(1):3–12.20225336 10.1002/path.2697
16. Jin Y Pang H Zhao L Zhao F Cheng Z Liu Q Ginseng total saponins and Fuzi total alkaloids exert antidepressant-like effects in ovariectomized mice through BDNF-mTORC1, autophagy and peripheral metabolic pathways Phytomedicine 2022 107 154425 10.1016/j.phymed.2022.154425 36137328
Jin Y, Pang H, Zhao L, Zhao F, Cheng Z, Liu Q, et al. Ginseng total saponins and Fuzi total alkaloids exert antidepressant-like effects in ovariectomized mice through BDNF-mTORC1, autophagy and peripheral metabolic pathways. Phytomedicine. 2022;107:154425.36137328 10.1016/j.phymed.2022.154425
17. Kato T Kato N Mitochondrial dysfunction in bipolar disorder Bipolar Disord 2000 2 3 Pt 1 180 90 10.1034/j.1399-5618.2000.020305.x 11256685
Kato T, Kato N. Mitochondrial dysfunction in bipolar disorder. Bipolar Disord. 2000;2(3 Pt 1):180–90.11256685 10.1034/j.1399-5618.2000.020305.x
18. Lee J Giordano S Zhang J Autophagy, mitochondria and oxidative stress: cross-talk and redox signalling Biochem J 2012 441 2 523 40 10.1042/BJ20111451 22187934
Lee J, Giordano S, Zhang J. Autophagy, mitochondria and oxidative stress: cross-talk and redox signalling. Biochem J. 2012;441(2):523–40.22187934 10.1042/BJ20111451
19. Bouman L Schlierf A Lutz AK Shan J Deinlein A Kast J Parkin is transcriptionally regulated by ATF4: evidence for an interconnection between mitochondrial stress and ER stress Cell Death Differ 2011 18 5 769 82 10.1038/cdd.2010.142 21113145
Bouman L, Schlierf A, Lutz AK, Shan J, Deinlein A, Kast J, et al. Parkin is transcriptionally regulated by ATF4: evidence for an interconnection between mitochondrial stress and ER stress. Cell Death Differ. 2011;18(5):769–82.21113145 10.1038/cdd.2010.142
20. Senft D Ronai ZA UPR, autophagy, and mitochondria crosstalk underlies the ER stress response Trends Biochem Sci 2015 40 3 141 8 10.1016/j.tibs.2015.01.002 25656104
Senft D, Ronai ZA. UPR, autophagy, and mitochondria crosstalk underlies the ER stress response. Trends Biochem Sci. 2015;40(3):141–8.25656104 10.1016/j.tibs.2015.01.002
21. Menéndez-Valle I, Cachán-Vega C, Boga JA, González-Blanco L, Antuña E, Potes Y et al. Differential Cellular Interactome in Schizophrenia and bipolar disorder-discriminatory biomarker role. Antioxid (Basel). 2023;12(11).
22. Barbosa IG Machado-Vieira R Soares JC Teixeira AL The immunology of bipolar disorder Neuroimmunomodulation 2014 21 2–3 117 22 10.1159/000356539 24557044
Barbosa IG, Machado-Vieira R, Soares JC, Teixeira AL. The immunology of bipolar disorder. Neuroimmunomodulation. 2014;21(2–3):117–22.24557044 10.1159/000356539
23. Breunis MN Kupka RW Nolen WA Suppes T Denicoff KD Leverich GS High numbers of circulating activated T cells and raised levels of serum IL-2 receptor in bipolar disorder Biol Psychiatry 2003 53 2 157 65 10.1016/S0006-3223(02)01452-X 12547472
Breunis MN, Kupka RW, Nolen WA, Suppes T, Denicoff KD, Leverich GS, et al. High numbers of circulating activated T cells and raised levels of serum IL-2 receptor in bipolar disorder. Biol Psychiatry. 2003;53(2):157–65.12547472 10.1016/S0006-3223(02)01452-X
24. Wu W Zheng YL Tian LP Lai JB Hu CC Zhang P Circulating T lymphocyte subsets, cytokines, and immune checkpoint inhibitors in patients with bipolar II or major depression: a preliminary study Sci Rep 2017 7 40530 10.1038/srep40530 28074937
Wu W, Zheng YL, Tian LP, Lai JB, Hu CC, Zhang P, et al. Circulating T lymphocyte subsets, cytokines, and immune checkpoint inhibitors in patients with bipolar II or major depression: a preliminary study. Sci Rep. 2017;7:40530.28074937 10.1038/srep40530
25. do Prado CH Rizzo LB Wieck A Lopes RP Teixeira AL Grassi-Oliveira R Reduced regulatory T cells are associated with higher levels of Th1/TH17 cytokines and activated MAPK in type 1 bipolar disorder Psychoneuroendocrinology 2013 38 5 667 76 10.1016/j.psyneuen.2012.08.005 22989476
do Prado CH, Rizzo LB, Wieck A, Lopes RP, Teixeira AL, Grassi-Oliveira R, et al. Reduced regulatory T cells are associated with higher levels of Th1/TH17 cytokines and activated MAPK in type 1 bipolar disorder. Psychoneuroendocrinology. 2013;38(5):667–76.22989476 10.1016/j.psyneuen.2012.08.005
26. Brambilla P Bellani M Isola M Bergami A Marinelli V Dusi N Increased M1/decreased M2 signature and signs of Th1/Th2 shift in chronic patients with bipolar disorder, but not in those with schizophrenia Transl Psychiatry 2014 4 7 e406 10.1038/tp.2014.46 24984193
Brambilla P, Bellani M, Isola M, Bergami A, Marinelli V, Dusi N, et al. Increased M1/decreased M2 signature and signs of Th1/Th2 shift in chronic patients with bipolar disorder, but not in those with schizophrenia. Transl Psychiatry. 2014;4(7):e406.24984193 10.1038/tp.2014.46
27. Ma Y Galluzzi L Zitvogel L Kroemer G Autophagy and cellular immune responses Immunity 2013 39 2 211 27 10.1016/j.immuni.2013.07.017 23973220
Ma Y, Galluzzi L, Zitvogel L, Kroemer G. Autophagy and cellular immune responses. Immunity. 2013;39(2):211–27.23973220 10.1016/j.immuni.2013.07.017
28. Jiang GM Tan Y Wang H Peng L Chen HT Meng XJ The relationship between autophagy and the immune system and its applications for tumor immunotherapy Mol Cancer 2019 18 1 17 10.1186/s12943-019-0944-z 30678689
Jiang GM, Tan Y, Wang H, Peng L, Chen HT, Meng XJ, et al. The relationship between autophagy and the immune system and its applications for tumor immunotherapy. Mol Cancer. 2019;18(1):17.30678689 10.1186/s12943-019-0944-z
29. Barrett T Wilhite SE Ledoux P Evangelista C Kim IF Tomashevsky M NCBI GEO: archive for functional genomics data sets–update Nucleic Acids Res 2013 41 Database issue D991 5 23193258
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res. 2013;41(Database issue):D991–5.23193258
30. Charoentong P Finotello F Angelova M Mayer C Efremova M Rieder D Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade Cell Rep 2017 18 1 248 62 10.1016/j.celrep.2016.12.019 28052254
Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 2017;18(1):248–62.28052254 10.1016/j.celrep.2016.12.019
31. Keenan AB Torre D Lachmann A Leong AK Wojciechowicz ML Utti V ChEA3: transcription factor enrichment analysis by orthogonal omics integration Nucleic Acids Res 2019 47 W1 W212 24 10.1093/nar/gkz446 31114921
Keenan AB, Torre D, Lachmann A, Leong AK, Wojciechowicz ML, Utti V, et al. ChEA3: transcription factor enrichment analysis by orthogonal omics integration. Nucleic Acids Res. 2019;47(W1):W212–24.31114921 10.1093/nar/gkz446
32. Clarke AJ Simon AK Autophagy in the renewal, differentiation and homeostasis of immune cells Nat Rev Immunol 2019 19 3 170 83 10.1038/s41577-018-0095-2 30531943
Clarke AJ, Simon AK. Autophagy in the renewal, differentiation and homeostasis of immune cells. Nat Rev Immunol. 2019;19(3):170–83.30531943 10.1038/s41577-018-0095-2
33. Deretic V Saitoh T Akira S Autophagy in infection, inflammation and immunity Nat Rev Immunol 2013 13 10 722 37 10.1038/nri3532 24064518
Deretic V, Saitoh T, Akira S. Autophagy in infection, inflammation and immunity. Nat Rev Immunol. 2013;13(10):722–37.24064518 10.1038/nri3532
34. Levine B Mizushima N Virgin HW Autophagy in immunity and inflammation Nature 2011 469 7330 323 35 10.1038/nature09782 21248839
Levine B, Mizushima N, Virgin HW. Autophagy in immunity and inflammation. Nature. 2011;469(7330):323–35.21248839 10.1038/nature09782
35. Deretic V Autophagy in inflammation, infection, and immunometabolism Immunity 2021 54 3 437 53 10.1016/j.immuni.2021.01.018 33691134
Deretic V. Autophagy in inflammation, infection, and immunometabolism. Immunity. 2021;54(3):437–53.33691134 10.1016/j.immuni.2021.01.018
36. Minchenko DO Kharkova AP Hubenia OV Minchenko OH Insulin receptor, IRS1, IRS2, INSIG1, INSIG2, RRAD, and BAIAP2 gene expressions in glioma U87 cells with ERN1 loss of function: effect of hypoxia and glutamine or glucose deprivation Endocr Regul 2013 47 1 15 26 10.4149/endo_2013_01_15 23363253
Minchenko DO, Kharkova AP, Hubenia OV, Minchenko OH. Insulin receptor, IRS1, IRS2, INSIG1, INSIG2, RRAD, and BAIAP2 gene expressions in glioma U87 cells with ERN1 loss of function: effect of hypoxia and glutamine or glucose deprivation. Endocr Regul. 2013;47(1):15–26.23363253 10.4149/endo_2013_01_15
37. Li R Zhao K Ruan Q Meng C Yin F Bone marrow mesenchymal stem cell-derived exosomal microRNA-124-3p attenuates neurological damage in spinal cord ischemia-reperfusion injury by downregulating Ern1 and promoting M2 macrophage polarization Arthritis Res Ther 2020 22 1 75 10.1186/s13075-020-2146-x 32272965
Li R, Zhao K, Ruan Q, Meng C, Yin F. Bone marrow mesenchymal stem cell-derived exosomal microRNA-124-3p attenuates neurological damage in spinal cord ischemia-reperfusion injury by downregulating Ern1 and promoting M2 macrophage polarization. Arthritis Res Ther. 2020;22(1):75.32272965 10.1186/s13075-020-2146-x
38. Wan H Wang Q Chen X Zeng Q Shao Y Fang H WDR45 contributes to neurodegeneration through regulation of ER homeostasis and neuronal death Autophagy 2020 16 3 531 47 10.1080/15548627.2019.1630224 31204559
Wan H, Wang Q, Chen X, Zeng Q, Shao Y, Fang H, et al. WDR45 contributes to neurodegeneration through regulation of ER homeostasis and neuronal death. Autophagy. 2020;16(3):531–47.31204559 10.1080/15548627.2019.1630224
39. Hetz C Saxena S ER stress and the unfolded protein response in neurodegeneration Nat Rev Neurol 2017 13 8 477 91 10.1038/nrneurol.2017.99 28731040
Hetz C, Saxena S. ER stress and the unfolded protein response in neurodegeneration. Nat Rev Neurol. 2017;13(8):477–91.28731040 10.1038/nrneurol.2017.99
40. Ogata M Hino S Saito A Morikawa K Kondo S Kanemoto S Autophagy is activated for cell survival after endoplasmic reticulum stress Mol Cell Biol 2006 26 24 9220 31 10.1128/MCB.01453-06 17030611
Ogata M, Hino S, Saito A, Morikawa K, Kondo S, Kanemoto S, et al. Autophagy is activated for cell survival after endoplasmic reticulum stress. Mol Cell Biol. 2006;26(24):9220–31.17030611 10.1128/MCB.01453-06
41. Yorimitsu T Nair U Yang Z Klionsky DJ Endoplasmic reticulum stress triggers autophagy J Biol Chem 2006 281 40 30299 304 10.1074/jbc.M607007200 16901900
Yorimitsu T, Nair U, Yang Z, Klionsky DJ. Endoplasmic reticulum stress triggers autophagy. J Biol Chem. 2006;281(40):30299–304.16901900 10.1074/jbc.M607007200
42. Li T Su L Zhong N Hao X Zhong D Singhal S Salinomycin induces cell death with autophagy through activation of endoplasmic reticulum stress in human cancer cells Autophagy 2013 9 7 1057 68 10.4161/auto.24632 23670030
Li T, Su L, Zhong N, Hao X, Zhong D, Singhal S, et al. Salinomycin induces cell death with autophagy through activation of endoplasmic reticulum stress in human cancer cells. Autophagy. 2013;9(7):1057–68.23670030 10.4161/auto.24632
43. Marques AP, Resende R, Silva DF, Batista M, Pereira D, Wildenberg B et al. Mitochondrial alterations in fibroblasts of early stage bipolar disorder patients. Biomedicines. 2021;9(5).
44. Sun Y Rong X Lu W Peng Y Li J Xu S Translational study of Alzheimer’s disease (AD) biomarkers from brain tissues in AβPP/PS1 mice and serum of AD patients J Alzheimers Dis 2015 45 1 269 82 10.3233/JAD-142805 25502766
Sun Y, Rong X, Lu W, Peng Y, Li J, Xu S, et al. Translational study of Alzheimer’s disease (AD) biomarkers from brain tissues in AβPP/PS1 mice and serum of AD patients. J Alzheimers Dis. 2015;45(1):269–82.25502766 10.3233/JAD-142805
45. Boutté AM Hook V Thangavelu B Sarkis GA Abbatiello BN Hook G Penetrating traumatic Brain Injury triggers dysregulation of cathepsin B protein levels Independent of Cysteine protease activity in brain and cerebral spinal fluid J Neurotrauma 2020 37 13 1574 86 10.1089/neu.2019.6537 31973644
Boutté AM, Hook V, Thangavelu B, Sarkis GA, Abbatiello BN, Hook G, et al. Penetrating traumatic Brain Injury triggers dysregulation of cathepsin B protein levels Independent of Cysteine protease activity in brain and cerebral spinal fluid. J Neurotrauma. 2020;37(13):1574–86.31973644 10.1089/neu.2019.6537
46. Man SM Kanneganti TD Regulation of lysosomal dynamics and autophagy by CTSB/cathepsin B Autophagy 2016 12 12 2504 5 10.1080/15548627.2016.1239679 27786577
Man SM, Kanneganti TD. Regulation of lysosomal dynamics and autophagy by CTSB/cathepsin B. Autophagy. 2016;12(12):2504–5.27786577 10.1080/15548627.2016.1239679
47. Hsu A Podvin S Hook V Lysosomal cathepsin protease gene expression profiles in the human brain during Normal Development J Mol Neurosci 2018 65 4 420 31 10.1007/s12031-018-1110-6 30008074
Hsu A, Podvin S, Hook V. Lysosomal cathepsin protease gene expression profiles in the human brain during Normal Development. J Mol Neurosci. 2018;65(4):420–31.30008074 10.1007/s12031-018-1110-6
48. Czibere L Baur LA Wittmann A Gemmeke K Steiner A Weber P Profiling trait anxiety: transcriptome analysis reveals cathepsin B (Ctsb) as a novel candidate gene for emotionality in mice PLoS ONE 2011 6 8 e23604 10.1371/journal.pone.0023604 21897848
Czibere L, Baur LA, Wittmann A, Gemmeke K, Steiner A, Weber P, et al. Profiling trait anxiety: transcriptome analysis reveals cathepsin B (Ctsb) as a novel candidate gene for emotionality in mice. PLoS ONE. 2011;6(8):e23604.21897848 10.1371/journal.pone.0023604
49. Saptarshi N Porter LF Paraoan L PERK/EIF2AK3 integrates endoplasmic reticulum stress-induced apoptosis, oxidative stress and autophagy responses in immortalised retinal pigment epithelial cells Sci Rep 2022 12 1 13324 10.1038/s41598-022-16909-6 35922637
Saptarshi N, Porter LF, Paraoan L. PERK/EIF2AK3 integrates endoplasmic reticulum stress-induced apoptosis, oxidative stress and autophagy responses in immortalised retinal pigment epithelial cells. Sci Rep. 2022;12(1):13324.35922637 10.1038/s41598-022-16909-6
50. Rege S Hodgkinson SJ Immune dysregulation and autoimmunity in bipolar disorder: synthesis of the evidence and its clinical application Aust N Z J Psychiatry 2013 47 12 1136 51 10.1177/0004867413499077 23908311
Rege S, Hodgkinson SJ. Immune dysregulation and autoimmunity in bipolar disorder: synthesis of the evidence and its clinical application. Aust N Z J Psychiatry. 2013;47(12):1136–51.23908311 10.1177/0004867413499077
51. Chen MH Chang WC Hsu JW Huang KL Tu PC Su TP Correlation of proinflammatory cytokines levels and reduced gray matter volumes between patients with bipolar disorder and unipolar depression J Affect Disord 2019 245 8 15 10.1016/j.jad.2018.10.106 30359810
Chen MH, Chang WC, Hsu JW, Huang KL, Tu PC, Su TP, et al. Correlation of proinflammatory cytokines levels and reduced gray matter volumes between patients with bipolar disorder and unipolar depression. J Affect Disord. 2019;245:8–15.30359810 10.1016/j.jad.2018.10.106
52. Naaldijk YM Bittencourt MC Sack U Ulrich H Kinins and microglial responses in bipolar disorder: a neuroinflammation hypothesis Biol Chem 2016 397 4 283 96 10.1515/hsz-2015-0257 26859499
Naaldijk YM, Bittencourt MC, Sack U, Ulrich H. Kinins and microglial responses in bipolar disorder: a neuroinflammation hypothesis. Biol Chem. 2016;397(4):283–96.26859499 10.1515/hsz-2015-0257
53. Rosenblat JD, McIntyre RS. Bipolar disorder and Immune Dysfunction: epidemiological findings, proposed pathophysiology and clinical implications. Brain Sci. 2017;7(11).
54. Drexhage RC Hoogenboezem TH Versnel MA Berghout A Nolen WA Drexhage HA The activation of monocyte and T cell networks in patients with bipolar disorder Brain Behav Immun 2011 25 6 1206 13 10.1016/j.bbi.2011.03.013 21443944
Drexhage RC, Hoogenboezem TH, Versnel MA, Berghout A, Nolen WA, Drexhage HA. The activation of monocyte and T cell networks in patients with bipolar disorder. Brain Behav Immun. 2011;25(6):1206–13.21443944 10.1016/j.bbi.2011.03.013
55. Knijff EM Ruwhof C de Wit HJ Kupka RW Vonk R Akkerhuis GW Monocyte-derived dendritic cells in bipolar disorder Biol Psychiatry 2006 59 4 317 26 10.1016/j.biopsych.2005.06.041 16165108
Knijff EM, Ruwhof C, de Wit HJ, Kupka RW, Vonk R, Akkerhuis GW, et al. Monocyte-derived dendritic cells in bipolar disorder. Biol Psychiatry. 2006;59(4):317–26.16165108 10.1016/j.biopsych.2005.06.041
56. Leite Dantas R, Freff J, Ambrée O, Beins EC, Forstner AJ, Dannlowski U et al. Dendritic cells: neglected modulators of Peripheral Immune responses and neuroinflammation in Mood disorders? Cells. 2021;10(4).
57. Chaurasia M Gupta S Das A Dwarakanath BS Simonsen A Sharma K Radiation induces EIF2AK3/PERK and ERN1/IRE1 mediated pro-survival autophagy Autophagy 2019 15 8 1391 406 10.1080/15548627.2019.1582973 30773986
Chaurasia M, Gupta S, Das A, Dwarakanath BS, Simonsen A, Sharma K. Radiation induces EIF2AK3/PERK and ERN1/IRE1 mediated pro-survival autophagy. Autophagy. 2019;15(8):1391–406.30773986 10.1080/15548627.2019.1582973
58. Brenner C Galluzzi L Kepp O Kroemer G Decoding cell death signals in liver inflammation J Hepatol 2013 59 3 583 94 10.1016/j.jhep.2013.03.033 23567086
Brenner C, Galluzzi L, Kepp O, Kroemer G. Decoding cell death signals in liver inflammation. J Hepatol. 2013;59(3):583–94.23567086 10.1016/j.jhep.2013.03.033
59. Morrone C, Smirnova NF, Jeridi A, Kneidinger N, Hollauer C, Schupp JC et al. Cathepsin B promotes collagen biosynthesis, which drives bronchiolitis obliterans syndrome. Eur Respir J. 2021;57(5).
