Background: Colorectal cancer is one of the most common malignant tumors worldwide. A various of neurotransmitter receptors have been found to be expressed in tumor cells, and the activation of these receptors may promote tumor growth and metastasis. This study aimed to construct a novel neurotransmitter receptor-related genes signature to predict the survival, immune microenvironment, and treatment response of colorectal cancer patients.Methods: RNA-seq and clinical data of colorectal cancer from The Cancer Genome Atlas database and Gene Expression Omnibus were downloaded. Neurotransmitter receptor-related gene were collected from publicly available data sources. The Weighted Gene Coexpression Network Analysis (WGCNA), Least Absolute Shrinkage and Selection Operator (LASSO) logistic regression, Support Vector Machine-Recursive Feature Elimination (SVM-RFE), and Random Forest (RF) algorithms were employed to construct the Neurotransmitter receptor-related gene prognostic signature. Further analyses, functional enrichment, CIBERSORTx, The Tumor Immune Single Cell Center (TISCH), survival analysis, and CellMiner, were performed to analyze immune status and treatment responses. Quantitative real-time polymerase chain reaction (qRT-PCR) assays were carried out to confirm the expression levels of prognostic genes.Results: By combining machine learning algorithm and WGCNA, we identified CHRNA3, GABRD, GRIK3, and GRIK5 as Neurotransmitter receptor-related prognostic genes signature. Functional enrichment analyses showed that these genes were enriched with cellular metabolic-related pathways, such as organic acid, inorganic acid, and lipid metabolism. CIBERSORTx and Single cell analysis showed that the high expression of genes were positively correlated with immunosuppressive cells infiltration, and the genes were mainly expressed in cancer-associated fibroblasts and endothelial cells. A nomogram was further built to predict overall survival (OS). The expression of CHRNA3, GABRD, GRIK3, and GRIK5 in cancer cells significantly impacted their response to chemotherapy.Conclusion: A neurotransmitter receptor-related prognostic gene signature was developed and validated in the current study, giving novel sights of neurotransmitter in predicting the prognostic and improving the treatment of CRC.

A merged human kidney scRNA-seq dataset made from the following scRNA-seq studies:

Lake, B.B. et al. A single-nucleus RNA-sequencing pipeline to decipher the molecular anatomy and pathophysiology of human kidneys. Nat Commun 10, 2832 (2019).

Liao, J., Yu, Z., Chen, Y. et al. Single-cell RNA sequencing of human kidney. Sci Data 7, 4 (2020).

Menon, R. et al. Single cell transcriptomics identifies focal segmental glomerulosclerosis remission endothelial biomarker. JCI Insight 5, e133267 (2020).

Wu, H. et al. Single-cell transcriptomics of a human kidney allograft biopsy specimen defines a diverse inflammatory response. J Am Soc Nephrol 29: 2069–2080 (2018).

Young, M. D. et al. Single-cell transcriptomes from human kidneys reveal the cellular identity of renal tumors. Science 361, 594–599 (2018).

Background: Despite the progress in early diagnosis and treatment, prognosis of pancreatic adenocarcinoma (PAAD) is still poor. Basic leucine zipper and W2 domain-containing protein 1 (BZW1) and protein 2 (BZW2) are attached to the basic leucine zipper (bZIP) superfamily. Recently, BZW1 was identified as an important role in glycolysis of PAAD. However, the comprehensive reports about BZW1/2 in PAAD are not sufficient.Methods: RNA-seq data in the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases were retrospectively analyzed. We explored the expression of BZW1/2 in PAAD tissues and the associations between BZW1/2 and prognosis. In addition, the potential roles of BZW1/2 in tumor microenvironment (TME) of PAAD were analyzed. Finally, clinicopathological data of 49 patients with PAAD in our institution were collected. Immunohistochemistry was used to determine the expression of BZW1/2 in PAAD samples.Results: BZW1 and BZW2 were upregulated in PAAD tissues compared to normal tissues (p < 0.05). The expression of BZW1/2 were not significantly correlated with gender, grade and stage of PAAD (p > 0.05). High expression of BZW2 was an independent predictor for poor prognosis of PAAD (HR 1.834, 95%CI 1.303–2.581, p = 0.001). And a nomogram to predict overall survival (OS) of PAAD was established with a C-index of 0.685. BZW1 and BZW2 expression were positively associated with T cell mediated immune response to tumor cell and Th2 cells in xCell database. Tumor Immune Single-Cell Hub (TISCH) analyses indicated that BZW1 and BZW2 were mainly expressed in B cells and malignant cells. External cohort furtherly validated that high expression of BZW1 and BZW2 were predictors for poor prognosis of PAAD.Conclusion: We found that BZW1 and BZW2 are highly expressed in malignant cells and B cells in the TME of PAAD. BZW2 is an independent predictor for OS of PAAD. BZW1 and BZW2 expression are positively associated with T cell mediated immune response to tumor cell and Th2 cells in PAAD.

BackgroundAlthough the functional damage of the diabetic pancreas can affect the postoperative recovery of pancreatic cancer patients, there is no significant difference in the prognosis of pancreatic cancer patients with a history of diabetes and ordinary pancreatic cancer patients. There is still no practical theory to explain this phenomenon.Materials and MethodThe mRNA expression profile data of 141 cases and 51 cases with clinical data of diabetes status were obtained from the TCGA database and the GEO database, respectively. The CRA001160 data set was obtained in the TISCH database. The Seurat was used to process single-cell expression profile sequencing data. The Cibersortx was used to construct a feature matrix of single-cell sequencing data and to deconvolve Bulk-RNAseq data to obtain each pancreatic cancer patients’ tumour invasion score. TIDE was used to assess the immune escape potential of the tumour. MiRNet was used to construct the miRNA-mRNA regulatory network.ResultCompared with regular pancreatic cancer patients, the immune-related signal transduction pathways in diabetic pancreatic cancer patients are in an activated state. In patients with diabetic pancreatic cancer, the infiltration score of CD8+ T cells is high, and the infiltration score of corresponding malignant tumour cells is low. The Bayesian classifier can distinguish diabetic pancreatic cancer patients from non-diabetic pancreatic cancer patients based on 10 signature genes. The miRNA-mRNA regulatory network suggests that regulation by miRNA can influence mRNA expression and thus prognostic survival of pancreatic cancer patients.ConclusionThe activation of inflammatory-related signalling pathways in diabetic pancreatic cancer patients increases the immune infiltration of CD8+ T cells in cancer patients and reduces the development of malignant tumour tissues. The expression of 10 signature genes allowed the diagnosis of diabetic and non-diabetic pancreatic cancer patients. The miRNA-mRNA regulatory network may be the main cause of the differences in the tumour inflammatory microenvironment between the two groups of patients. These findings help us further understand the immune microenvironment of patients with diabetic pancreatic cancer.

Background: Despite the progress in early diagnosis and treatment, prognosis of pancreatic adenocarcinoma (PAAD) is still poor. Basic leucine zipper and W2 domain-containing protein 1 (BZW1) and protein 2 (BZW2) are attached to the basic leucine zipper (bZIP) superfamily. Recently, BZW1 was identified as an important role in glycolysis of PAAD. However, the comprehensive reports about BZW1/2 in PAAD are not sufficient.Methods: RNA-seq data in the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases were retrospectively analyzed. We explored the expression of BZW1/2 in PAAD tissues and the associations between BZW1/2 and prognosis. In addition, the potential roles of BZW1/2 in tumor microenvironment (TME) of PAAD were analyzed. Finally, clinicopathological data of 49 patients with PAAD in our institution were collected. Immunohistochemistry was used to determine the expression of BZW1/2 in PAAD samples.Results: BZW1 and BZW2 were upregulated in PAAD tissues compared to normal tissues (p < 0.05). The expression of BZW1/2 were not significantly correlated with gender, grade and stage of PAAD (p > 0.05). High expression of BZW2 was an independent predictor for poor prognosis of PAAD (HR 1.834, 95%CI 1.303–2.581, p = 0.001). And a nomogram to predict overall survival (OS) of PAAD was established with a C-index of 0.685. BZW1 and BZW2 expression were positively associated with T cell mediated immune response to tumor cell and Th2 cells in xCell database. Tumor Immune Single-Cell Hub (TISCH) analyses indicated that BZW1 and BZW2 were mainly expressed in B cells and malignant cells. External cohort furtherly validated that high expression of BZW1 and BZW2 were predictors for poor prognosis of PAAD.Conclusion: We found that BZW1 and BZW2 are highly expressed in malignant cells and B cells in the TME of PAAD. BZW2 is an independent predictor for OS of PAAD. BZW1 and BZW2 expression are positively associated with T cell mediated immune response to tumor cell and Th2 cells in PAAD.

Background: Hexokinase 3 (HK3) is one of the key enzymes involved in glucose phosphorylation (the first step in most glucose metabolic pathways). Many studies have demonstrated the vital role of dysregulation of HK3 in several tumors. However, there is a need for in-depth characterization of the role of HK3 in glioblastoma multiforme (GBM).Methods: All data were sourced from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA). Kaplan-Meier analysis and univariate regression were applied for survival analysis. Gene set enrichment analysis (GSEA) was used for enrichment analysis. Tumor Immune Single Cell Hub (TISCH) database was applied for single-cell analysis. Tumor Immune Dysfunction and Exclusion (TIDE) analysis was applied to evaluate the immune response.Results: HK3 expression was upregulated in GBM and correlated with poor prognosis. The high HK3 expression group was primarily enriched in adaptive immune response, chemokine signaling pathway, and cytokine-cytokine receptor interaction. The high HK3 expression group showed significantly greater enrichment of the majority of immune cells and immune-related pathways. HK3 showed significant correlation with most immune cells, especially macrophages (p < .001, R = .81). TISCH analysis showed that HK3 was predominantly expressed in macrophages in most cancers. HK3 showed significant correlation with most immune-related genes, such as PD-1 (p < .001, R = .41), PDL-1 (p < .001, R = .27), and CTLA-4 (p < .001, R = .29). TIDE analysis revealed that the low HK3 expression group has a lower TIDE score and may benefit from immunotherapy. Drug sensitivity analysis showed that patients with high HK3 expression frequently showed drug resistance.Conclusion: HK3 was associated with poor prognosis and may serve as a biomarker of macrophages in GBM. HK3 was also associated with immune response and drug resistance. Our findings may provide novel insights for GBM immunotherapy.

Background: Hexokinase 3 (HK3) is one of the key enzymes involved in glucose phosphorylation (the first step in most glucose metabolic pathways). Many studies have demonstrated the vital role of dysregulation of HK3 in several tumors. However, there is a need for in-depth characterization of the role of HK3 in glioblastoma multiforme (GBM).Methods: All data were sourced from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA). Kaplan-Meier analysis and univariate regression were applied for survival analysis. Gene set enrichment analysis (GSEA) was used for enrichment analysis. Tumor Immune Single Cell Hub (TISCH) database was applied for single-cell analysis. Tumor Immune Dysfunction and Exclusion (TIDE) analysis was applied to evaluate the immune response.Results: HK3 expression was upregulated in GBM and correlated with poor prognosis. The high HK3 expression group was primarily enriched in adaptive immune response, chemokine signaling pathway, and cytokine-cytokine receptor interaction. The high HK3 expression group showed significantly greater enrichment of the majority of immune cells and immune-related pathways. HK3 showed significant correlation with most immune cells, especially macrophages (p < .001, R = .81). TISCH analysis showed that HK3 was predominantly expressed in macrophages in most cancers. HK3 showed significant correlation with most immune-related genes, such as PD-1 (p < .001, R = .41), PDL-1 (p < .001, R = .27), and CTLA-4 (p < .001, R = .29). TIDE analysis revealed that the low HK3 expression group has a lower TIDE score and may benefit from immunotherapy. Drug sensitivity analysis showed that patients with high HK3 expression frequently showed drug resistance.Conclusion: HK3 was associated with poor prognosis and may serve as a biomarker of macrophages in GBM. HK3 was also associated with immune response and drug resistance. Our findings may provide novel insights for GBM immunotherapy.

BackgroundSuccinylation, a key post-translational modification, plays a crucial role in metabolic regulation and tumor progression. However, its influence on the tumor immune microenvironment and its prognostic implications remain unclear. A systematic pan-cancer analysis of succinylation-related molecular activities is needed.MethodsBulk transcriptomic, single-cell RNA sequencing, and spatial transcriptomic data across pan-cancer from TCGA, GEO, TISCH, and multiple other databases were analyzed. Succinylation scores were calculated using Gene Set Variation Analysis (GSVA). The interactions between succinylation scores, immune infiltration, tumor microenvironment, tumor mutational burden, and immunotherapy response were assessed. A succinylation-based prognostic model was constructed and validated in colorectal cancer (CRC) cohorts. PCED1A protein expression was evaluated by immunohistochemistry and Western blotting. The function of PCED1A in CRC was investigated through in vitro experiments.ResultsSuccinylation scores were significantly altered in multiple tumor types. Higher succinylation scores correlated with mitochondrial oxidative phosphorylation, while lower succinylation scores were linked to immune cell differentiation. Spatial transcriptomic analysis showed a negative correlation between succinylation scores and immune cell activity in tumor-adjacent regions. A prognostic model consisting of 11 succinylation-related genes (ATP6V1C2, CAPS, DAPK1, P4HA1, PCED1A, RASL10B, AGT, EREG, HYAL1, SARAF, and SLC4A4) was developed. High-risk patients exhibited significantly shorter overall survival. PCED1A was upregulated in CRC and positively associated with SIRT5. Overexpression of PCED1A promoted intracellular protein desuccinylation, along with enhanced CRC cell proliferation, migration, and invasion.ConclusionOur analysis demonstrates that succinylation-related molecular activities display distinct expression patterns across cancers, which are associated with metabolic regulation, immune modulation, and disease prognosis. The succinylation-based prognostic model provides a novel risk stratification tool for CRC, while PCED1A-dependent succinylation regulation may serve as a potential therapeutic target.

Background: ARLs, which are a class of small GTP-binding proteins, play a crucial role in facilitating tumor tumorigenesis and development. ARL4C, a vital member of the ARLs family, has been implicated in the progression of tumors, metastatic dissemination, and development of resistance to therapeutic drugs. Nevertheless, the precise functional mechanisms of ARL4C concerning tumor prognosis and immunotherapy drug susceptibility remain elusive.Methods: By combining the GTEx and TCGA databases, the presence of ARL4C was examined in 33 various types of cancer. Immunohistochemistry and immunofluorescence staining techniques were utilized to confirm the expression of ARL4C in particular tumor tissues. Furthermore, the ESTIMATE algorithm and TIMER2.0 database were utilized to analyze the tumor microenvironment and immune infiltration associated with ARL4C. The TISCH platform facilitated the utilization of single-cell RNA-seq datasets for further analysis. ARL4C-related immune escape was investigated using the TISMO tool. Lastly, drug sensitivity analysis was conducted to assess the sensitivity of different types of tumors to compounds based on the varying levels of ARL4C expression.Results: The study found that ARL4C was highly expressed in 23 different types of cancer. Moreover, the presence of high ARL4C expression was found to be associated with a poor prognosis in BLCA, COAD, KIRP, LGG, and UCEC. Notably, ARL4C was also expressed in immune cells, and its high expression was found to be correlated with cancer immune activation. Most importantly, the drug sensitivity analysis revealed a positive correlation between ARL4C expression and the heightened sensitivity of tumors to Staurosporine, Midostaurin, and Nelarabine.Conclusion: The findings from our study indicate that the expression level of ARL4C may exert an influence on cancer development, prognosis, and susceptibility to immunotherapy drugs. In addition, the involvement of ARL4C in the tumor immune microenvironment has expanded the concept of ARL4C-targeted immunotherapy.

Background: Hexokinase 3 (HK3) is one of the key enzymes involved in glucose phosphorylation (the first step in most glucose metabolic pathways). Many studies have demonstrated the vital role of dysregulation of HK3 in several tumors. However, there is a need for in-depth characterization of the role of HK3 in glioblastoma multiforme (GBM).Methods: All data were sourced from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA). Kaplan-Meier analysis and univariate regression were applied for survival analysis. Gene set enrichment analysis (GSEA) was used for enrichment analysis. Tumor Immune Single Cell Hub (TISCH) database was applied for single-cell analysis. Tumor Immune Dysfunction and Exclusion (TIDE) analysis was applied to evaluate the immune response.Results: HK3 expression was upregulated in GBM and correlated with poor prognosis. The high HK3 expression group was primarily enriched in adaptive immune response, chemokine signaling pathway, and cytokine-cytokine receptor interaction. The high HK3 expression group showed significantly greater enrichment of the majority of immune cells and immune-related pathways. HK3 showed significant correlation with most immune cells, especially macrophages (p < .001, R = .81). TISCH analysis showed that HK3 was predominantly expressed in macrophages in most cancers. HK3 showed significant correlation with most immune-related genes, such as PD-1 (p < .001, R = .41), PDL-1 (p < .001, R = .27), and CTLA-4 (p < .001, R = .29). TIDE analysis revealed that the low HK3 expression group has a lower TIDE score and may benefit from immunotherapy. Drug sensitivity analysis showed that patients with high HK3 expression frequently showed drug resistance.Conclusion: HK3 was associated with poor prognosis and may serve as a biomarker of macrophages in GBM. HK3 was also associated with immune response and drug resistance. Our findings may provide novel insights for GBM immunotherapy.

Breast cancer (BC) is a malignant tumor that occurs in breast tissue. This project aims to predict the prognosis of BC patients using genes related to hypoxia and endoplasmic reticulum stress (ERS). RNA-seq and clinical data for BC were downloaded from TCGA and GEO databases. Hypoxia and ERS-related genes were collected from the Genecards database. Univariate/multivariate Cox regression and Lasso regression analyses were used to screen genes and construct prognostic models. Patients were divided into high-risk (HR) and low-risk (LR) groups based on risk scores. The CIBERSORT algorithm was used to analyze differences in immune infiltration between the two groups. The mutations of the two groups were analyzed statistically. The CellMiner database was used for drug prediction and the TISCH database for single-cell sequencing analysis. We screened 8 feature genes to construct a prognostic model. Patients in the HR group had a remarkably worse prognosis. TP53 exhibited a higher mutation frequency in the HR group. CIBERSORT analysis uncovered a remarkable increase in the infiltration levels of Macrophages M0 and Tregs in cancer patients and HR patients. Drug sensitivity prediction demonstrated that the expression of IVL was greatly negatively linked with the sensitivity of COLCHICINE. PTGS2 had a remarkably negative correlation with the Vincristine sensitivity. The prognostic model based on 8 hypoxia and ERS-related genes can predict the survival, immune status, and potential drugs of BC patients, bringing a new perspective on individualized treatment.

Glioma is a highly common pathological brain tumor. Misfolded protein response, which is strongly associated with the growth of cancerous tumors, is mediated by the gene, endoplasmic reticulum degradation-enhancing alpha-mannosidase-like protein 2. However, this gene has not been linked to glioma. To assess the same, we used The Cancer Genome Atlas, Chinese Glioma Genome Atlas, and Genotype-Tissue Expression datasets. The gene was overexpressed in gliomas. This overexpression was linked to unfavorable clinical characteristics, such as the World Health Organization grade, isocitrate dehydrogenase mutation, and the combined loss of the short arm chromosome 1 and the long arm of chromosome 19. Quantitative polymerase chain reaction experiments and immunohistochemistry on clinical samples from our institution verified the gene’s expression and clinical importance. The Human Protein Atlas website verified the messenger ribonucleic acid expression of the gene in glioma cell lines, and immunohistochemistry verified the presence of its protein. A previous survival study indicated that its high expression is substantially related to a bad prognosis. It was identified as an independent predictor of primary glioma prognosis using multivariate Cox regression analysis. To forecast individual survival, we created a nomogram based on this (concordance-index = 0.847). Additionally, functional annotation demonstrated its major role in the control of the extracellular matrix and immune system. The scratch assay and transwell migration assay confirmed the decreased invasive ability of U251 glioma cells with the gene knockdown. Its increased expression was found to be related to the extent of macrophage infiltration using the CIBERSORT, ESTIMATE, Single-sample Gene Set Enrichment Analysis, and Tumor Immune Single-Cell Hub (TISCH) algorithms. The Tumor Immune Dysfunction and Exclusion algorithm revealed that the gene can accurately predict the response of immunotherapy (area under the receiver operating characteristic curve = 0.857). Further, isocitrate dehydrogenase 1 mutation is typically more frequent when the gene expression is high. Finally, five medicines targeting this gene were discovered utilizing the molecular docking program and drug sensitivity analysis of the RNAactDrug website. Low expression of the gene inhibited glioma cell invasion. Therefore, the gene is helpful for the diagnosis, prognosis, and case-specific immunotherapy of glioma.

BackgroundGlia maturation factor-γ (GMFG) regulates actin cytoskeletal organization and promotes the invasion of cancer cells. However, its expression pattern and molecular function in gliomas have not been clearly defined.MethodsIn this study, public datasets comprising 2,518 gliomas samples were used to explore GMFG expression and its correlation with malignancy in gliomas. Immunohistochemistry (IHC) staining was performed to determine the expression of GMFG in gliomas using an in-house cohort that contained 120 gliomas samples. Gene ontology enrichment analysis was conducted using the DAVID tool. The correlation between GMFG expression and immune cell infiltration was evaluated using TIMER, Tumor Immune Single-Cell Hub (TISCH) database, and IHC staining assays. The Kaplan–Meier analysis was performed to determine the prognostic role of GMFG and its association with temozolomide (TMZ) response in gliomas.ResultsThe GMFG expression was higher in gliomas compared with non-tumor brain tissues both in public datasets and in-house cohort. High expression of GMFG was significantly associated with WHO grade IV, IDH 1/2 wild-type, and mesenchymal (ME) subtypes. Bioinformatic prediction and IHC analysis revealed that GMFG expression obviously correlated with the macrophage marker CD163 in gliomas. Moreover, both lower grade glioma (LGG) and glioblastoma multiforme (GBM) patients with high GMFG expression had shorter overall survival than those with low GMFG expression. These results indicate that GMFG may be a therapeutic target for the treatment of such patients. Patients with low GMFG expression who received chemotherapy had a longer survival time than those with high GMFG expression. For patients who received ion radiotherapy (IR) only, the GMFG expression level had no effect on the overall survival neither in CGGA and TCGA datasets.ConclusionThe GMFG is a novel prognostic biomarker for patients with both LGG and GBM. Increased GMFG expression is associated with tumor-associated macrophages (TAMs) infiltration and with a bad response to TMZ treatment.

BackgroundR3HDM1, an RNA binding protein with one R3H domain, remains uncharacterized in terms of its association with tumor progression, malignant cell regulation, and the tumor immune microenvironment. This paper aims to fill this gap by analyzing the potential of R3HDM1 in diagnosis, prognosis, chemotherapy, and immune function across various cancers.MethodsData was collected from the Firehost database (http://gdac.broadinstitute.org) to obtain the TCGA pan-cancer queue containing tumor and normal samples. Additional data on miRNA, TCPA, mutations, and clinical information were gathered from the UCSC Xena database (https://xenabrowser.net/datapages/). The mutation frequency and locus of R3HDM1 in the TCGA database were examined using the cBioPortal. External validation through GEO data was conducted to assess the differential expression of R3HDM1 in different cancers. Protein expression levels were evaluated using the Clinical Proteomics Tumor Analysis Alliance (CPTAC). The differential expression of R3HDM1 was verified in lung adenocarcinoma cell lines and normal lung glandular epithelial cells via RT-qPCR. Cell migration and proliferation experiments were conducted by knocking down the expression of R3HDM1 in two lung adenocarcinoma cell lines using small interfering RNA. The biological role of R3HDM1 in pan-cancer was explored using the GSEA method. Multiple immune infiltration algorithms from the TIMER2.0 database was employed to investigate the correlation between R3HDM1 expression and the tumor immune microenvironment. Validation of transcriptome immune infiltration was based on 140 single-cell datasets from the TISCH database. The study also characterized a pan-cancer survival profile and analyzed the differential expression of R3HDM1 in different molecular subtypes. The relationship between R3HDM1 and drug resistance was investigated using four chemotherapy data sources: CellMiner, GDSC, CTRP and PRISM. The impact of chemicals on the expression of R3HDM1 was explored through the CTD database.ResultThe study revealed differential expression of R3HDM1 in various tumors, indicating its potential as an early diagnostic marker. Changes in somatic copy number (SCNA) and DNA methylation were identified as factors contributing to abnormal expression levels. Additionally, the study found that R3HDM1 expression is associated with clinical features, metabolic pathways, and important pathways related to metastasis and the immune system. High expression of R3HDM1 was linked to poor prognosis across different tumors and altered drug sensitivity. Furthermore, the expression of R3HDM1 showed significant correlations with immune modulatory molecules and biomarkers of lymphocyte subpopulation infiltration. Finally, the study highlighted four chemicals that could influence the expression of R3HDM1.ConclusionOverall, this study proposes that R3HDM1 expression is a promising biomarker for predicting the prognosis of cancer, especially lung adenocarcinoma, and the efficacy of immunotherapy, demonstrating the rationale for further exploration in the development of anti-tumor therapies.

BackgroundDERL3 has been implicated as an essential element in the degradation of misfolded lumenal glycoproteins induced by endoplasmic reticulum (ER) stress. However, the correlation of DERL3 expression with the malignant phenotype of lung adenocarcinoma (LUAD) cells is unclear and remains to be elucidated. Herein, we investigated the interaction between the DERL3 and LUAD pathological process.MethodsThe Cancer Genome Atlas (TCGA) database was utilized to determine the genetic alteration of DERL3 in stage I LUAD. Clinical LUAD samples including carcinoma and adjacent tissues were obtained and were further extracted to detect DERL3 mRNA expression via RT-qPCR. Immunohistochemistry was performed to evaluate the protein expression of DERL3 in LUAD tissues. The GEPIA and TIMER website were used to evaluate the correlation between DERL3 and immune cell infiltration. We further used the t-SNE map to visualize the distribution of DERL3 in various clusters at the single-cell level via TISCH database. The potential mechanisms of the biological process mediated by DERL3 in LUAD were conducted via KEGG and GSEA.ResultsIt was indicated that DERL3 was predominantly elevated in carcinoma compared with adjacent tissues in multiple kinds of tumors from the TCGA database, especially in LUAD. Immunohistochemistry validated that DERL3 was also upregulated in LUAD tissues compared with adjacent tissues from individuals. DERL3 was preliminarily found to be associated with immune infiltration via the TIMER database. Further, the t-SNE map revealed that DERL3 was predominantly enriched in plasma cells of the B cell population. It was demonstrated that DERL3 high-expressed patients presented significantly worse response to chemotherapy and immunotherapy. GSEA and KEGG results indicated that DERL3 was positively correlated with B cell activation and unfolded protein response (UPR).ConclusionOur findings indicated that DERL3 might play an essential role in the endoplasmic reticulum-associated degradation (ERAD) process in LUAD. Moreover, DERL3 may act as a promising immune biomarker, which could predict the efficacy of immunotherapy in LUAD.