PDIA5 is Correlated With Immune Infiltration and Predicts Poor Prognosis in Gliomas

Comprehensive analysis of lncRNA-miRNA-mRNA ceRNA network and key genes in granulosa cells of patients with biochemical primary ovarian insufficiency

Primary ovarian insufficiency (POI) is a common condition leading to the pathological decline of ovarian function in women of reproductive age, resulting in amenorrhea, hypogonadism, and infertility. Biochemical premature ovarian insufficiency (bPOI) is an intermediate stage in the pathogenesis of POI in which the fertility of patients has been reduced. Previous studies suggest that granulosa cells (GCs) play an essential role in the pathogenesis of POI, but their pathogenetic mechanisms remain unclear. To further explore the potential pathophysiological mechanisms of GCs in POI, we constructed a molecular long non-coding RNA (lncRNA)-microRNA (miRNA)-messenger RNA (mRNA) network using GC expression data collected from biochemical premature ovarian failure (bPOI) patients in the GEO database. We discovered that the GCs of bPOI patients had differential expression of 131 mRNAs, 191 lncRNAs, and 28 miRNAs. By systematic network analysis, we identified six key genes, including SRSF1, PDIA5, NEURL1B, UNK, CELF2, and CFL2, and five hub miRNAs, namely hsa-miR-27a-3p, hsa-miR-24-3p, hsa-miR-22-3p, hsa-miR-129-5p, and hsa-miR-17-5p, and the results suggest that the expression of these key genes may be regulated by two hub miRNAs, hsa-miR-27a-3p and hsa-miR-17-5p. Additionally, a POI model in vitro was created to confirm the expression of a few important genes. In this study, we discovered a unique lncRNA-miRNA-mRNA network based on the ceRNA mechanism in bPOI for the first time, and we screened important associated molecules, providing a partial theoretical foundation to better understand the pathogenesis of POI.

BCL2A1 is associated with tumor-associated macrophages and unfavorable prognosis in human gliomas

B-cell lymphoma 2-related protein A1 (BCL2A1) is a member of the BCL-2 family. Previous studies have shown that BCL2A1 is closely related to the tumorigenesis and resistance to chemotherapy of multiple solid tumors, such as breast cancer. However, the expression pattern and potential biological function of BCL2A1 in glioma remain unknown. For the first time, we found that the expression of BCL2A1 was higher in human glioma tissues than in normal brain tissues (NBTs) in both public datasets and an in-house cohort. High BCL2A1 expression was associated with advanced WHO grade, IDH 1/2 wild type and the mesenchymal (ME) subtype, and its overexpression in glioma predicted resistance to temozolomide (TMZ) chemotherapy and unfavorable prognosis. In addition, Gene set enrichment analysis (GSEA), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that BCL2A1 was significantly correlated with the immune response and immune-related pathways, and BCL2A1 expression was positively correlated with microenvironmental parameters (immune, stromal, and ESTIMATE scores) and macrophage infiltration. Interestingly, bioinformatic prediction and immunohistochemical/immunofluorescence staining analysis revealed that BCL2A1 expression was obviously associated with the tumor-associated macrophages (TAMs) markers CD68 and CCL2. Notably, knockdown of BCL2A1 significantly inhibited cell proliferation of U87 and U251 in vitro, induced smaller tumor size and prolonged survival time of mice in vivo. Co-culture experiments of macrophages and GBM cells showed that BCL2A1 knockdown inhibited macrophage migration. Meanwhile, knockdown of BCL2A1 was associated with low expression of CD68 and CCL2 in intracranial xenograft model. This may suggest that BCL2A1 promotes the progression of glioma and influences the prognosis of patients by participating in TAMs infiltration. In conclusion, these findings suggest that BCL2A1 could serve as a promising prognostic indicator and immunotherapy target in gliomas.

The evaluation of six genes combined value in glioma diagnosis and prognosis

PURPOSE

Glioma is the most common and fatal type of brain tumour. Owing to its aggressiveness and lethality, early diagnosis and prediction of patient survival are very important. This study aimed to identify key genes and biomarkers for glioma that can guide clinicians in making rapid diagnosis and prognostication.

METHODS

Data mining of The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), Repository of Molecular Brain Neoplasia Data, and Genotype-Tissue Expression Project brain expression data revealed significantly differentially expressed genes (DEGs), and the risk scores of individual patients were calculated. WGCNA was utilized to screen for genes most related to clinical diagnosis. Prognostic genes associated with glioma were selected via combining the LASSO regression with univariate and multivariate Cox regression and protein-protein interaction network analyses. Then, a nomogram was constructed. And CGGA dataset was utilized to validated. The protein expression levels of the signature were detected using the human protein atlas. Drug response prediction was carried out using the package "pRRophetic".

RESULTS

A six-gene signature (KLF6, CHI3L1, SERPINE1, ANGPT2, TGFBR1, and PTX3) was identified and used to stratify patients into low- and high-risk groups. Survival, ROC curve, and Cox analyses clarified that the six hub genes were a favourable independent prognostic factor for patients with glioma. A nomogram was set up by integrating clinical parameters with risk signatures, showing high precision for predicting 2-, 3-, 4-, 5-years survival. In addition, the expression of most genes was consistent with protein expression. Furthermore, the sensitivity to the top ten drugs in the GDSC database of the high-risk group was significantly higher than the low-risk group.

CONCLUSION

Based on genetic profiles and clinicopathological features, including age, grade, isocitrate dehydrogenase mutation status, we constructed a comprehensive prognostic model for patients with glioma. These signatures can be regarded as biomarkers to predict the prognosis of gliomas, possibly providing more therapeutic strategies for future clinical research.

HELLPAR/RRM2 axis related to HMMR as novel prognostic biomarker in gliomas

BACKGROUND

Gliomas are the most frequent type of central nervous system tumor, accounting for more than 70% of all malignant CNS tumors. Recent research suggests that the hyaluronan-mediated motility receptor (HMMR) could be a novel potential tumor prognostic marker. Furthermore, mounting data has highlighted the important role of ceRNA regulatory networks in a variety of human malignancies. The complexity and behavioural characteristics of HMMR and the ceRNA network in gliomas, on the other hand, remained unknown.

METHODS

Transcriptomic expression data were collected from TCGA, GTEx, GEO, and CGGA database.The relationship between clinical variables and HMMR was analyzed with the univariate and multivariate Cox regression. Kaplan-Meier method was used to assess OS. TCGA data are analyzed and processed, and the correlation results obtained were used to perform GO, GSEA, and ssGSEA. Potentially interacting miRNAs and lncRNAs were predicted by miRWalk and StarBase.

RESULTS

HMMR was substantially expressed in gliomas tissues compared to normal tissues. Multivariate analysis revealed that high HMMR expression was an independent predictive predictor of OS in TCGA and CGGA. Functional enrichment analysis found that HMMR expression was associated with nuclear division and cell cycle. Base on ssGSEA analysis, The levels of HMMR expression in various types of immune cells differed significantly. Bioinformatics investigation revealed the HEELPAR-hsa-let-7i-5p-RRM2 ceRNA network, which was linked to gliomas prognosis. And through multiple analysis, the good predictive performance of HELLPAR/RRM2 axis for gliomas patients was confirmed.

CONCLUSION

This study provides multi-layered and multifaceted evidence for the importance of HMMR and establishes a HMMR-related ceRNA (HEELPAR-hsa-let-7i-5p-RRM2) overexpressed network related to the prognosis of gliomas.

IP-score correlated to endogenous tumour antigen peptide processing: A candidate clinical response score algorithm of immune checkpoint inhibitors therapy in multiple cohorts

The processing of endogenous tumour antigen peptides was essential for anti-tumour immunity in the tumour microenvironment. A high degree of Endogenous tumour antigen peptide processing has been demonstrated to improve the prognosis of carcinoma patients. However, there is insufficient evidence to prove its effect on the clinical response to immune checkpoint inhibitor therapy. To undertake a more in-depth analysis of the effects of the aforementioned genes on immunotherapy, we constructed a gene set evaluation score system relevant to tumour endogenous antigen peptide therapy using the GSVA approach. This rating mechanism is known as IP score (IPs). Immediately afterwards, we used the TCGA pan-cancer cohorts to conduct a comprehensive analysis of 6 genes in the IPs, and the analysis results showed that these six genes were related to the proportion of CD8⁺ T lymphocytes in a variety of solid tumours. As a prognostic protective factor for solid tumours, patients had better prognosis outcomes in the group with high expression levels of the above genes. We analysed the differential expression of six genes between immune checkpoint inhibitor treatment response and disease progression groups using several treatment cohorts. The results revealed that after treatment with PD-1 or CTLA4 inhibitors, the expression levels of the above six genes were comparatively high in the effective group, but the expression of the signature genes was dramatically downregulated in the ICI-insensitive groups. This indicates that the 6 genes are related to the clinical response to ICI treatment. Finally, we used the GSVA method to evaluate the above signatures, and the results showed that PDCD1, CTAL4, CD274 and LAG3 were significantly higher expressed in the IPs high-expression group; therefore, based on the processing of endogenous antigenic peptides in tumours, a predictive score of clinical response to immune checkpoint inhibitor therapy composed of 6 genes(PSMB8/PSMB9/PSMB10/PSME1/PSME2/IRF1) was constructed, and the role of each independent variable in the signature in the solid tumour microenvironment and the impact on ICI treatment were comprehensively analysed. This study provides a candidate evaluation score for predicting clinical response to immune checkpoint inhibitor therapy.

RUNX1/CD44 axis regulates the proliferation, migration, and immunotherapy of gliomas: A single-cell sequencing analysis

Background

Glioma is one of the most common, primary, and lethal adult brain tumors because of its extreme aggressiveness and poor prognosis. Several recent studies relevant to the immune function of CD44, a transmembrane glycoprotein as a significant hyaluronic acid receptor, have achieved great success, revealing the critical role of CD44 in immune infiltration in gliomas. The overexpression of CD44 has been verified to correlate with cancer aggressiveness and migration, while the clinical and immune features of CD44 expression have not yet been thoroughly characterized in gliomas.

Methods

Molecular and clinical data of glioma collected from publicly available genomic databases were analyzed.

Results

CD44 was up-expressed in malignant gliomas, notably in the 1p/19q non-codeletion cases, isocitrate dehydrogenase (IDH) wild-type, and mesenchymal subtypes in GBM samples. CD44 expression level strongly correlates with stromal and immune cells, mainly infiltrating the glioma microenvironment by single-cell sequencing analysis. Meanwhile, CD44 can be a promising biomarker in predicting immunotherapy responses and mediating the expression of PD-L1. Finally, RUNX1/CD44 axis could promote the proliferation and migration of gliomas.

Conclusions

Therefore, CD44 was responsible for glioma growth and progression. It could potentially lead to a novel target for glioma immunotherapy or a prognostic biomarker.

Single-cell and microarray chip analysis revealed the underlying pathogenesis of ulcerative colitis and validated model genes in diagnosis and drug response

The morbidity rate of ulcerative colitis (UC) in the world is increasing year by year, recurrent episodes of diarrhea, mucopurulent and bloody stools, and abdominal pain are the main symptoms, reducing the quality of life of the patient and affecting the productivity of the society. In this study, we sought to develop robust diagnostic biomarkers for UC, to uncover potential targets for anti-TNF-ɑ drugs, and to investigate their associated pathway mechanisms. We collected single-cell expression profile data from 9 UC or healthy samples and performed cell annotation and cell communication analysis. Revealing the possible pathogenesis of ulcerative colitis by Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA) analysis. Based on the disease-related modules obtained from weighted correlation network analysis (WGCNA) analysis, we used Lasso regression analysis and random forest algorithm to identify the genes with the greatest impact on disease (EPB41L3, HSD17B3, NDRG1, PDIA5, TRPV3) and further validated the diagnostic value of the model genes by various means. To further explore the relationship and mechanism between model genes and drug sensitivity, we collected gene expression profiles of 185 UC patients before receiving anti-tumor necrosis factor drugs, and we performed functional analysis based on the results of differential analysis between NR tissues and R tissues, and used single-sample GSEA (ssGSEA) and CIBERSORT algorithms to explore the important role of immune microenvironment on drug sensitivity. The results suggest that our model is not only helpful in aiding diagnosis, but also has implications for predicting drug efficacy; in addition, model genes may influence drug sensitivity by affecting immune cells. We suggest that this study has developed a diagnostic model with higher specificity and sensitivity, and also provides suggestions for clinical administration and drug efficacy prediction.

Machine learning-based identification of SOX10 as an immune regulator of macrophage in gliomas

Gliomas, originating from the glial cells, are the most lethal type of primary tumors in the central nervous system. Standard treatments like surgery have not significantly improved the prognosis of glioblastoma patients. Recently, immune therapy has become a novel and effective option. As a conserved group of transcriptional regulators, the Sry-type HMG box (SOX) family has been proved to have a correlation with numerous diseases. Based on the large-scale machine learning, we found that the SOX family, with significant immune characteristics and genomic profiles, can be divided into two distinct clusters in gliomas, among which SOX10 was identified as an excellent immune regulator of macrophage in gliomas. The high expression of SOX10 is related to a shorter OS in LGG, HGG, and pan-cancer groups but benefited from the immunotherapy. It turned out in single-cell sequencing that SOX10 is high in neurons, M1 macrophages, and neural stem cells. Also, macrophages are found to be elevated in the SOX10 high-expression group. SOX10 has a positive correlation with macrophage cytokine production and negative regulation of macrophages' chemotaxis and migration. In conclusion, our study demonstrates the outstanding cluster ability of the SOX family, indicating that SOX10 is an immune regulator of macrophage in gliomas, which can be an effective target for glioma immunotherapy.

Identification of CD73 as a Novel Biomarker Encompassing the Tumor Microenvironment, Prognosis, and Therapeutic Responses in Various Cancers

CD73 is essential in promoting tumor growth by prohibiting anti-tumor immunity in many cancer types. While the mechanism remains largely unknown, our paper comprehensively confirmed the onco-immunological characteristics of CD73 in the tumor microenvironment (TME) of pan-cancer. This paper explored the expression pattern, mutational profile, prognostic value, tumor immune infiltration, and response to immunotherapy of CD73 in a continuous cohort of cancers through various computational tools. The co-expression of CD73 on cancer cells, immune cells, and stromal cells in the TME was also detected. Especially, we examined the correlation between CD73 and CD8⁺ (a marker of T cell), CD68⁺ (a marker of macrophage), and CD163⁺ (a marker of M2 macrophage) cells using multiplex immunofluorescence staining of tissue microarrays. CD73 expression is significantly associated with a patient's prognosis and could be a promising predictor of these cancers. High CD73 levels are strongly linked to immune infiltrations, neoantigens, and immune checkpoint expression in the TME. In particular, enrichment signaling pathway analysis demonstrated that CD73 was obviously related to activation pathways of immune cells, including T cells, macrophages, and cancer-associated fibroblasts (CAFs). Meanwhile, single-cell sequencing algorithms found that CD73 is predominantly co-expressed on cancer cells, CAFs, M2 macrophages, and T cells in several cancers. In addition, we explored the cellular communication among 14 cell types in glioblastoma (GBM) based on CD73 expression. Based on the expression of CD73 as well as macrophage and T cell markers, we predicted the methylation and enrichment pathways of these markers in pan-cancer. Furthermore, a lot of therapeutic molecules sensitive to these markers were predicted. Finally, potential anticancer inhibitors, immunotherapies, and gene therapy responses targeting CD73 were identified from a series of immunotherapy cohorts. CD73 is closely linked to clinical prognosis and immune infiltration in many cancers. Targeting CD73-dependent signaling pathways may be a promising therapeutic strategy for future tumor immunotherapy.

Repurposing live attenuated trivalent MMR vaccine as cost-effective cancer immunotherapy

It has long been known that oncolytic viruses wield their therapeutic capability by priming an inflammatory state within the tumor and activating the tumor immune microenvironment, resulting in a multifaceted antitumor immune response. Vaccine-derived viruses, such as measles and mumps, have demonstrated promising potential for treating human cancer in animal models and clinical trials. However, the extensive cost of manufacturing current oncolytic viral products makes them far out of reach for most patients. Here by analyzing the impact of intratumoral (IT) administrations of the trivalent live attenuated measles, mumps, and rubella viruses (MMR) vaccine, we unveil the cellular and molecular basis of MMR-induced anti-cancer activity. Strikingly, we found that IT delivery of low doses of MMR correlates with tumor control and improved survival in murine hepatocellular cancer and colorectal cancer models via increased tumor infiltration of CD8+ granzyme B+ T-cells and decreased macrophages. Moreover, our data indicate that MMR activates key cellular effectors of the host’s innate and adaptive antitumor immunity, culminating in an immunologically coordinated cancer cell death. These findings warrant further work on the potential for MMR to be repurposed as safe and cost-effective cancer immunotherapy to impact cancer patients globally.

Identify the immune characteristics and immunotherapy value of CD93 in the pan-cancer based on the public data sets

CD93 is a transmembrane receptor that is mainly expressed on endothelial cells. A recent study found that upregulated CD93 in tumor vessels is essential for tumor angiogenesis in several cancers. However, the underlying mechanisms are largely unexplored. Our present research systematically analyzed the characteristics of CD93 in tumor immunotherapy among 33 cancers. CD93 levels and co-expression of CD93 on cancer and stromal cells were detected using public databases and multiple immunofluorescence staining. The Kaplan-Meier (KM) analysis identified the predictive role of CD93 in these cancer types. The survival differences between CD93 mutants and WT, CNV groups, and methylation were also investigated. The immune landscape of CD93 in the tumor microenvironment was analyzed using the SangerBox, TIMER 2.0, and single-cell sequencing. The immunotherapy value of CD93 was predicted through public databases. CD93 mRNA and protein levels differed significantly between cancer samples and adjacent control tissues in multiply cancer types. CD93 mRNA expression associated with patient prognosis in many cancers. The correlation of CD93 levels with mutational status of other gene in these cancers was also analyzed. CD93 levels significantly positively related to three scores (immune, stromal, and extimate), immune infiltrates, immune checkpoints, and neoantigen expression.. Additionally, single-cell sequencing revealed that CD93 is predominantly co-expressed on tumor and stromal cells, such as endothelial cells, cancer-associated fibroblasts (CAFs), neutrophils, T cells, macrophages, M1 and M2 macrophages. Several immune-related signaling pathways were enriched based on CD93 expression, including immune cells activation and migration, focal adhesion, leukocyte transendothelial migration, oxidative phosphorylation, and complement. Multiple immunofluorescence staining displayed the relationship between CD93 expression and CD8, CD68, and CD163 in these cancers. Finally, the treatment response of CD93 in many immunotherapy cohorts and sensitive small molecules was predicted from the public datasets. CD93 expression is closely associated with clinical prognosis and immune infiltrates in a variety of tumors. Targeting CD93-related signaling pathways in the tumor microenvironment may be a novel therapeutic strategy for tumor immunotherapy.

Recent developments in application of single-cell RNA sequencing in the tumour immune microenvironment and cancer therapy

The advent of single-cell RNA sequencing (scRNA-seq) has provided insight into the tumour immune microenvironment (TIME). This review focuses on the application of scRNA-seq in investigation of the TIME. Over time, scRNA-seq methods have evolved, and components of the TIME have been deciphered with high resolution. In this review, we first introduced the principle of scRNA-seq and compared different sequencing approaches. Novel cell types in the TIME, a continuous transitional state, and mutual intercommunication among TIME components present potential targets for prognosis prediction and treatment in cancer. Thus, we concluded novel cell clusters of cancer-associated fibroblasts (CAFs), T cells, tumour-associated macrophages (TAMs) and dendritic cells (DCs) discovered after the application of scRNA-seq in TIME. We also proposed the development of TAMs and exhausted T cells, as well as the possible targets to interrupt the process. In addition, the therapeutic interventions based on cellular interactions in TIME were also summarized. For decades, quantification of the TIME components has been adopted in clinical practice to predict patient survival and response to therapy and is expected to play an important role in the precise treatment of cancer. Summarizing the current findings, we believe that advances in technology and wide application of single-cell analysis can lead to the discovery of novel perspectives on cancer therapy, which can subsequently be implemented in the clinic. Finally, we propose some future directions in the field of TIME studies that can be aided by scRNA-seq technology.

MS4A6A is a new prognostic biomarker produced by macrophages in glioma patients

MS4A6A has been recognized as being associated with aging and the onset of neurodegenerative disease. However, the mechanisms of MS4A6A in glioma biology and prognosis are ill-defined. Here, we show that MS4A6A is upregulated in glioma tissues, resulting in unfavorable clinical outcomes and poor responses to adjuvant chemotherapy. Multivariate Cox regression analysis suggested that MS4A6A expression can act as a strong and independent predictor for glioma outcomes (CGGA1: HR: 1.765, p < 0.001; CGGA2: HR: 2.626, p < 0.001; TCGA: HR: 1.415, p < 0.001; Rembrandt: HR: 1.809, p < 0.001; Gravendeel: HR: 1.613, p < 0.001). A protein–protein interaction (PPI) network revealed that MS4A6A might be coexpressed with CD68, CD163, and macrophage-specific signatures. Enrichment analysis showed the innate immune response and inflammatory response to be markedly enriched in the high MS4A6A expression group. Additionally, single-cell RNA sequencing (scRNA-seq) analysis revealed distinctive expression features for MS4A6A in macrophages in the glioma immune microenvironment (GIME). Immunofluorescence staining confirmed colocalization of CD68/MS4A6A and CD163/MS4A6A in macrophages. Correlation analysis revealed that MS4A6A expression is positively related to the tumor mutation burden (TMB) of glioma, displaying the high potential of applying MS4A6A to evaluate responsiveness to immunotherapy. Altogether, our research indicates that MS4A6A upregulation may be used as a promising and effective indicator for adjuvant therapy and prognosis assessment.

Pan-Cancer Analysis of the Immunological Role of PDIA5: A Potential Target for Immunotherapy

The aberrant protein disulfide isomerase A5 (PDIA5) expression was relevant to the poor prognosis of patients with human cancers. However, its relationship with the epigenetic and genetic alterations and its effect on tumor immunity is still lacking. In the present study, we comprehensively analyzed the immune infiltration role of PDIA5 in human cancers based on large-scale bioinformatics analyses and in vitro experiments. Obvious DNA methylation and moderate alteration frequency of PDIA5 were observed in human cancers. The expression level of PDIA5 was significantly correlated with infiltrated immune cells, immune pathways, and other immune signatures. We found that cancer cells and macrophages exhibited high PDIA5 expression in human cancers using the single-cell RNA sequencing analysis. We also demonstrated the interaction between PDIA5 and immune cells in glioblastoma multiforme (GBM). Multiplex immunofluorescence staining showed the upregulated expression level of PDIA5 and the increased number of M2 macrophage markers-CD163 positive cells in pan-cancer samples. Notably, PDIA5 silencing resulted in upregulated expression of PD-L1 and SPP1 in U251 cells. Silencing of PDIA5 in hepG2 cells, U251 cells, and PC3 cells contributed to a decline in their ability of proliferation, clone formation, and invasion and inhibited the migration of cocultured M2 macrophages. Additionally, PDIA5 also displayed predictive value in the immunotherapy response of both murine and human cancer cohorts. Overall, our findings indicated that PDIA5 might be a potential target for immunotherapies in cancers.

Proteomic Analysis Reveals Molecular Differences in the Development of Gastric Cancer

Gastric cancer (GC) is the 3rd leading cause of death from cancer and the 5th most common cancer worldwide. The detection rate of GC among Tibetans is significantly higher than that in Han Chinese, probably due to differences in their living habits, dietary structure, and environment. Despite such a high disease burden, the epidemiology of gastric cancer has not been studied in this population. Molecular markers are required to aid the diagnosis and treatment of GC. In this study, we collected gastric tissue samples from patients in Tibet with chronic nonatrophic gastritis (CNAG) (n = 6), chronic atrophic gastritis (CAG) (n = 7), gastric intraepithelial neoplasia (GIN) (n = 4), and GC (n = 5). The proteins in each group were analyzed using coupled label-free mass spectrometry. In addition, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, and protein interaction networks were used to analyze the differentially expressed proteins (DEPs) among groups. DEPs were quantified in comparisons of GC versus CNAG (223), GC versus GIN (100), and GIN versus CNAG (341). GO and KEGG analyses showed that the DEPs were mainly associated with immunity (GC versus CNAG) and cancer proliferation and metastasis (GC versus GIN, and GIN versus CNAG). Furthermore, the expression levels of cell proliferation and cytoskeleton-related proteins increased consistently during cancer development, such as ITGA4, DDC, and CPT1A; thus, they are potential diagnostic markers. These results obtained by proteomics analysis could improve our understanding of cancer biology in GC and provide a rich resource for data mining and discovering potential immunotherapy targets.

A Novel Thrombosis-Related Signature for Predicting Survival and Drug Compounds in Glioblastoma

Glioblastoma is the most common primary tumor in the central nervous system, and thrombosis-associated genes are related to its occurrence and progression. Univariate Cox and LASSO regression analysis were utilized to develop a new prognostic signature based on thrombosis-associated genes. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and HALLMARK were used for functional annotation of risk signature. ESTIMATE, MCP-counter, xCell, and TIMER algorithms were used to quantify immune infiltration in the tumor microenvironment. Genomics of Drug Sensitivity in Cancer (GDSC) was used for selecting potential drug compounds. Risk signature based on thrombosis-associated genes shows moderate performance in prognosis prediction. The functional annotation of the risk signature indicates that the signaling pathways related to the cell cycle, apoptosis, tumorigenesis, and immune suppression are rich in the high-risk group. Somatic mutation analysis shows that tumor-suppressive gene TP53 and oncogene PTEN have higher expression in low-risk and high-risk groups, respectively. Potential drug compounds are explored in risk score groups and show higher AUC values in the low-risk score group. A nomogram with valuable prognostic factors exhibits high sensitivity in predicting the survival outcome of GBM patients. Our research screens out multiple thromboses-associated genes with remarkable clinical significance in GBM and further develops a meaningful prognostic risk signature predicting drug sensitivity and survival outcome.

Role of CD68 in tumor immunity and prognosis prediction in pan-cancer

CD68 plays a critical role in promoting phagocytosis; however, the function of CD68 in tumor immunity and prognosis remains unknown. We analyzed CD68 expression among 33 tumor and normal tissues from The Cancer Genome Atlas and Genotype-Tissue Expression datasets. The relationship between CD68 expression and cancer prognosis, immune infiltration, checkpoint markers, and drug response was explored. Upregulated CD68 levels were observed in various cancer types, which were verified through tumor tissue chips using immunohistochemistry. High levels of CD68 in tumor samples correlated with an adverse prognosis in glioblastoma, kidney renal clear cell carcinoma, lower-grade glioma, liver hepatocellular carcinoma, lung squamous cell carcinoma, thyroid carcinoma, and thymoma and a favorable prognosis in kidney chromophobe. The top three negatively enriched Kyoto Encyclopedia of Genes and Genomes terms in the high CD68 subgroup were chemokine signaling pathway, cytokine-cytokine receptor interaction, and cell adhesion molecule cams. The top negatively enriched HALLMARK terms included complement, allograft rejection, and inflammatory response. A series of targeted drugs and small-molecule drugs with promising therapeutic effects were predicted. The clinical prognosis and immune infiltration of high expression levels of CD68 differ across tumor types. Inhibiting CD68-dependent signaling could be a promising therapeutic strategy for immunotherapy in many tumor types.

Prognostic and Diagnostic Values of Semaphorin 5B and Its Correlation With Tumor-Infiltrating Immune Cells in Kidney Renal Clear-Cell Carcinoma

Background: Semaphorin 5B (SEMA5B) has been described to be involved in the development and progression of cancer. However, the potential diagnostic and prognosis roles and its correlation with tumor-infiltrating immune cells in KIRC have not been clearly reported yet. Methods: The mRNA level of SEMA5B was analyzed via the TCGA and GTEx database as well as the CCLE dataset and verified by GSE53757 and GSE40435 datasets. Meanwhile, the protein level of SEMA5B was analyzed by CPTAC and validated by HPA. The diagnostic value of SEMA5B was analyzed according to the TCGA database and validated by GSE53757, GSE46699, and GSE11024 + GSE46699 datasets. Then, the survival analysis was conducted using GEPIA2. R software (v3.6.3) was applied to investigate the relevance between SEMA5B and immune checkpoints and m6A RNA methylation regulator expression. The correlation between SEMA5B and MMRs and DNMT expression and tumor-infiltrating immune cells was explored via TIMER2. Co-expressed genes of SEMA5B were assessed by cBioPortal, and enrichment analysis was conducted by Metascape. The methylation analysis was conducted with MEXPRESS and MethSurv online tools. Gene set enrichment analysis (GSEA) was applied to annotate the biological function of SEMA5B. Results: SEMA5B was significantly upregulated at both the mRNA and protein levels in KIRC. Further analysis demonstrated that the mRNA expression of SEMA5B was significantly correlated with gender, age, T stage, pathologic stage, and histologic grade. High levels of SEMA5B were found to be a favorable prognostic factor and novel diagnostic biomarker for KIRC. SEMA5B expression was shown to be significantly associated with the abundance of immune cells in KIRC. Also, SEMA5B expression was significantly correlated with the abundance of MMR genes, DNMTs, and m6A regulators in KIRC. Enrichment analysis indicated that the co-expressed genes may involve in crosslinking in the extracellular matrix (ECM). GSEA disclosed that SYSTEMIC_LUPUS_ERYTHEMATOSUS and NABA_ECM_REGULATORS were prominently enriched in the SEMA5B low-expression phenotype. Finally, the methylation analysis demonstrated a correlation between hypermethylation of the SEMA5B gene and a poor prognosis in KIRC. Conclusion: Increased SEMA5B expression correlated with immune cell infiltration, which can be served as a favorable prognostic factor and a novel diagnostic biomarker for KIRC.

Antigen Presentation Machinery Signature-Derived CALR Mediates Migration, Polarization of Macrophages in Glioma and Predicts Immunotherapy Response

Immunogenicity, influenced by tumor antigenicity and antigen presenting efficiency, critically determines the effectiveness of immune checkpoint inhibitors. The role of immunogenicity has not been fully elucidated in gliomas. In this study, a large-scale bioinformatics analysis was performed to analyze the prognostic value and predictive value of antigen presentation machinery (APM) signature in gliomas. ssGSEA algorithm was used for development of APM signature and LASSO regression analysis was used for construction of APM signature-based risk score. APM signature and risk score showed favorable performance in stratifying survival and predicting tumorigenic factors of glioma patients. APM signature and risk score were also associated with different genomic features in both training cohort TCGA and validating cohort CGGA. Furthermore, APM signature-based risk score was independently validated in three external cohorts and managed to predict immunotherapy response. A prognostic nomogram was constructed based on risk score. Risk score-derived CALR was found to mediate the invasion and polarization of macrophages based on the coculture of HMC3 and U251 cells. CALR could significantly predict immunotherapy response. In conclusion, APM signature and APM signature-based risk score could help promote the clinical management of gliomas.

CMTM Family Genes Affect Prognosis and Modulate Immunocytes Infiltration in Grade II/III Glioma Patients by Influencing the Tumor Immune Landscape and Activating Associated Immunosuppressing Pathways

Lower-grade glioma (LGG) is one of the most common primary tumor types in adults. The chemokine-like factor (CKLF)-like Marvel transmembrane domain-containing (CMTM) family is widely expressed in the immune system and can modulate tumor progression. However, the role of the CMTM family in LGG remains unknown. A total of 508 LGG patients from The Cancer Genome Atlas (TCGA) database were used as a training cohort, and 155 LGG patients from the Chinese Glioma Genome Atlas (CGGA) array database, 142 LGG patients from the CGGA RNA-sequencing database, and 168 LGG patients from the GSE108474 database were used as the validation cohorts. Patients were subdivided into two groups using consensus clustering. The ENET algorithm was applied to build a scoring model based on the cluster model. Finally, ESTIMATE, CIBERSORT, and xCell algorithms were performed to define the tumor immune landscape. The expression levels of the CMTM family genes were associated with glioma grades and isocitrate dehydrogenase (IDH) status. Patients in cluster 2 and the high-risk score group exhibited a poor prognosis and were enriched with higher grade, wild-type IDH (IDH-WT), 1p19q non-codeletion, MGMT promoter unmethylation, and IDH-WT subtype. Patients in cluster 1 and low-risk score group were associated with high tumor purity and reduced immune cell infiltration. Enrichment pathways analysis indicated that several essential pathways involved in tumor progression were associated with the expression of CMTM family genes. Importantly, PD-1, PD-L1, and PD-L2 expression levels were increased in cluster 2 and high-risk groups. Therefore, the CMTM family contributes to LGG progression through modulating tumor immune landscape.

Functions of RNF Family in the Tumor Microenvironment and Drugs Prediction in Grade II/III Gliomas

Increasing evidence has demonstrated that RING finger (RNF) proteins played a vital role in cellular and physiological processes and various diseases. However, the function of RNF proteins in low-grade glioma (LGG) remains unknown. In this study, 138 RNF family members revealed their role in LGG. The TCGA database was used as the training cohort; two CGGA databases and GSE108474 were selected as external validation cohorts. Patients were grouped into cluster 1 and cluster 2, both in the training and validation cohorts, using consensus clustering analysis. The prognosis of patients in cluster 1 is significantly better than that in cluster 2. Meanwhile, biofunction prediction was further introduced to explore the potential mechanisms that led to differences in survival outcomes. Patients in Cluster 2 showed more complicated immunocytes infiltration and highly immunosuppressive features than cluster 1. Enrichment pathways such as negative regulation of mast cell activation, DNA replication, mismatch repair, Th17 cell differentiation, antigen processing and presentation, dendritic cell antigen processing and presentation, dendritic cell differentiation were also enriched in cluster 2 patients. For the last, the main contributors were distinguished by employing a machine learning algorithm. A lot of targeted and small molecule drugs that are sensitive to patients in cluster 2 were predicted. Importantly, we discovered TRIM8, DTX2, and TRAF5 as the most vital contributors from the RNF family, which were related to immune infiltration in LGG tumor immune landscape. In this study, we demonstrated the predicted role of RNF proteins in LGG. In addition, we found out three markers among RNF proteins that are closely related to the immune aspects of LGG, which might serve as novel therapeutic targets for immunotherapy in the future.

Integrin Alpha-2 as a Potential Prognostic and Predictive Biomarker for Patients With Lower-Grade Glioma

Diffuse gliomas are the most common malignant brain tumors with the highest mortality and recurrence rate in adults. Integrin alpha-2 (ITGA2) is involved in a series of biological processes, including cell adhesion, stemness regulation, angiogenesis, and immune/blood cell functions. The role of ITGA2 in lower-grade gliomas (LGGs) is not well defined. Firstly, we downloaded RNA sequencing and relevant clinical information from The Cancer Genome Atlas cohort, the Chinese Glioma Genome Atlas cohort, and related immune cohorts. Next, prognosis analysis, difference analysis, clinical model construction, enrichment analysis, and immune infiltration analysis are performed for this study. These analyses indicated that ITGA2 may have clinical application value and research value in LGG immunotherapy. We also detected the mRNA and protein expression of ITGA2 in three LGG cell lines and normal glial cells using quantitative real-time polymerase chain reaction assay and western blot assay. Our study not only offers a novel target for LGG immunotherapy but also can better comprehend the mechanism of the development and progression of patients with LGG. This study revealed that ITGA2 may be a potential prognostic and predictive biomarker for LGG, which can bring new insights into targeted immunotherapy.

TNFSF13 Is a Novel Onco-Inflammatory Marker and Correlates With Immune Infiltration in Gliomas

Existing therapeutic strategies for gliomas are restricted; hence, exploration for novel diagnostic indicator and treatment is essential. Here, we performed bioinformatic analyses for TNFSF13 (also known as APRIL), a proliferation-inducing ligand of the tumor necrosis factor (TNF) superfamily, aiming to assess its potential for predicting glioma patient’s prognosis and targeted therapy. TNFSF13 expression was upregulated in the increase of tumor grades based on Xiangya cohort. In high TNFSF13 gliomas, somatic mutation was proved to correlate with amplification of EGFR and deletion of CDKN2A; while mutation of IDH1 was more frequently observed in low TNFSF13 group. We also confirmed the positive correlation between TNFSF13 and infiltrating immune and stromal cells in glioma microenvironment. Further, TNFSF13 was found to be involved in immunosuppression via diverse immunoregulation pathways and was associated with other immune checkpoints and inflammation. Single-cell sequencing revealed an abundant expression of TNFSF13 in neoplastic cells and M2 macrophages, which TNFSF13 might potentially regulate the cell communication via IL-8, C3, and CD44. Lastly, TNFSF13 mediated the activities of transcription factors including FOXO3, MEIS2, and IRF8. Our analyses demonstrated the relevance between TNFSF13 and glioma progress and indicated the potential of TNFSF13 as a novel diagnostic onco-inflammatory biomarker and immunotherapy target of gliomas.

Novel Immune Infiltrating Cell Signature Based on Cell Pair Algorithm Is a Prognostic Marker in Cancer

Tumor-infiltrating immune cells (TIICs) have become an important source of markers for predicting the clinical outcomes of cancer patients. However, measurements of cellular heterogeneity vary due to the frequently updated reference genomes and gene annotations. In this study, we systematically collected and evaluated the infiltration pattern of 65 immune cells. We constructed the Immune Cell Pair (ICP) score based on the cell pair algorithm in 3,715 samples and across 12 independent cancer types, among which, the ICP score from six cancer types was further validated in 2,228 GEO samples. An extensive tumorigenic and immunogenomic analysis was subsequently conducted. As a result, the ICP score showed a robust reliability and efficacy in predicting the survival of patients with gliomas, in pan-cancer samples, and six independent cancer types. Notably, the ICP score was correlated with the genomic alteration features in gliomas. Moreover, the ICP score exhibited a remarkable association with multiple immunomodulators that could potentially mediate immune escape. Finally, the ICP score predicted immunotherapeutic responses with a high sensitivity, allowing a useful tool for predicting the overall survival and guiding immunotherapy for cancer patients.

Immune Infiltrating Cells-Derived Risk Signature Based on Large-scale Analysis Defines Immune Landscape and Predicts Immunotherapy Responses in Glioma Tumor Microenvironment

The glioma tumor microenvironment (TME), composed of several noncancerous cells and biomolecules is known for its complexity of cancer-immune system interaction. Given that, novel risk signature is required for predicting glioma patient responses to immunotherapy. In this study, we systematically evaluated the TME infiltration pattern of 2877 glioma samples. TME phenotypes were determined using the Partitioning Around Medoid method. Machine learning including SVM-RFE and Principal component analysis (PCA) were used to construct a TME scoring system. A total of 857 glioma samples from four datasets were used for external validation of the TME-score. The correlation of TME phenotypes and TME-scores with diverse clinicopathologic characteristics, genomic features, and immunotherapeutic efficacy in glioma patients was determined. Immunohistochemistry staining for the M2 macrophage marker CD68 and CD163, mast cell marker CD117, neutrophil marker CD66b, and RNA sequencing of glioma samples from the XYNS cohort were performed. Two distinct TME phenotypes were identified. High TME-score correlated with a high number of immune infiltrating cells, elevated expression of immune checkpoints, increased mutation rates of oncogenes, and poor survival of glioma patients. Moreover, high TME-score exhibited remarkable association with multiple immunomodulators that could potentially mediate immune escape of cancer. Thus, the TME-score showed the potential to predict the efficacy of anti-PD-1 immunotherapy. Univariate and multivariate analyses demonstrated the TME-score to be a valuable prognostic biomarker for gliomas. Our study demonstrated that TME could potentially influence immunotherapy efficacy in melanoma patients whereas its role in immunotherapy of glioma patients remains unknown. Therefore, a better understanding of the TME landscape in gliomas would promote the development of novel immunotherapy strategies against glioma.

HOXA5 Is Recognized as a Prognostic-Related Biomarker and Promotes Glioma Progression Through Affecting Cell Cycle

Glioma is malignant tumor derives from glial cells in the central nervous system. High-grade glioma shows aggressive growth pattern, and conventional treatments, such as surgical removal and chemo-radiotherapy, archive limitation in the interference of this process. In this work, HOXA5, from the HOX family, was identified as a glioma cell proliferation-associated factor by investigating its feature in the TCGA and CGGA data set. High HOXA5 expression samples contain unfavorable clinical features of glioma, including IDH wild type, un-methylated MGMT status, non-codeletion 1p19q status, malignant molecular subtype. Survival analysis indicates that high HOXA5 expression samples are associated with worse clinical outcome. The CNVs and SNPs profile difference further confirmed the enrichment of glioma aggressive related biomarkers. In the meantime, the activation of DNA damage repair-related pathways and TP53-related pathways is also related to HOXA5 expression. In cell lines, U87MG and U251, by interfering HOXA5 expression significantly inhibit glioma progression and apoptosis, and cell cycle is arrested at the G2/M phase. Collectively, increased HOXA5 expression can promote glioma progression via affecting glioma cell proliferation.

The molecular feature of macrophages in tumor immune microenvironment of glioma patients

Background

Gliomas are one of the most common types of primary tumors in central nervous system. Previous studies have found that macrophages actively participate in tumor growth.

Methods

Weighted gene co-expression network analysis was used to identify meaningful macrophage-related gene genes for clustering. Pamr, SVM, and neural network were applied for validating clustering results. Somatic mutation and methylation were used for defining the features of identified clusters. Differentially expressed genes (DEGs) between the stratified groups after performing elastic regression and principal component analyses were used for the construction of MScores. The expression of macrophage-specific genes were evaluated in tumor microenvironment based on single cell sequencing analysis. A total of 2365 samples from 15 glioma datasets and 5842 pan-cancer samples were used for external validation of MScore.

Results

Macrophages were identified to be negatively associated with the survival of glioma patients. Twenty-six macrophage-specific DEGs obtained by elastic regression and PCA were highly expressed in macrophages at single-cell level. The prognostic value of MScores in glioma was validated by the active proinflammatory and metabolic profile of infiltrating microenvironment and response to immunotherapies of samples with this signature. MScores managed to stratify patient survival probabilities in 15 external glioma datasets and pan-cancer datasets, which predicted worse survival outcome. Sequencing data and immunohistochemistry of Xiangya glioma cohort confirmed the prognostic value of MScores. A prognostic model based on MScores demonstrated high accuracy rate.

Conclusion

Our findings strongly support a modulatory role of macrophages, especially M2 macrophages in glioma progression and warrants further experimental studies.

Gene Model Related to m6A Predicts the Prognostic Effect of Immune Infiltration on Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a highly aggressive solid tumor. Because most studies have focused on the intrinsic carcinogenic pathways of tumors, we focused on the relationship between N6-methyladenosine (m6A) and the prognosis of HNSCC in the tumor immune microenvironment. We downloaded RNA-seq data from the TCGA dataset and used univariate Cox regression to screen m6A-related lncRNAs. The expression value of LASSO-screened genes was the sum of LASSO regression coefficients. We then evaluated relationships between the risk score and cellular components or cellular immune response. Differences in immune response under various algorithms were visualized with heat maps. The GSVA package in R was used to analyze GO, BP, KEGG, and hallmark gene sets of immune checkpoint clusters and immune checkpoint scores. The GSEA analysis was performed with the cluster profile package, yielding 21 m6A genes. Related lncRNAs were screened with Pearson's correlations, and the resulting 442 lncRNAs were screened using single-factor analysis. Eight lncRNAs closely related to prognosis were identified through survival random forest. Survival analysis showed that patients with a high risk score had a poor prognosis. Low- and high-risk-score groups differed significantly in m6A gene expression. Prognostic scores from different algorithms were significantly correlated with B cells, T cells, and memory cells in the immune microenvironment. Expression of immune checkpoints and signal pathways differed significantly across risk-score groups, suggesting that m6A could mediate lncRNA-induced immune system dysfunction and affect HNSCC development. A comprehensive study of tumor-cell immune characteristics should provide more insight into the complex immune microenvironment, thus contributing to the development of new immunotherapeutic agents.

Molecular Characterization of AEBP1 at Transcriptional Level in Glioma

Background

Glioma is the most common malignant tumor of the brain in adult patients. The standardized treatment protocol is based on surgical therapy, supplemented with radiotherapy and chemotherapy. However, the prognosis is still unsatisfied. Chemoresistance is one of the most important reason for the poor prognosis of glioma patients. It has confirmed that glioma stem cell (GSC) is one of the reasons for chemoresistance.

Methods

In this study, three datasets (GSE23806, COSMIC, and TCGA) were used to perform the analysis to search for the key genes related to GSC, temozolomide (TMZ) resistance, and prognosis. The key gene for further research was selected by reviewing the previous studies. The selected gene investigated the relation between expression levels and clinical characteristics in both TCGA and CGGA dataset. The bioinformatics analysis was performed by Gene Ontology (GO) analysis. The survival analysis was performed by Kaplan–Meier survival analysis.

Results

AE binding protein 1 (AEBP1) was selected for further analysis. AEBP1 was overexpressed in GSCs and TMZ resistance cells. In both TCGA and CGGA dataset, the results showed that the expression level of AEBP1 was increased in glioblastoma (GBM) samples, IDH wild-type samples, and MGMT promoter unmethylated samples. Meanwhile, AEBP1 expression was positively related to several GSC markers. GO analysis showed that AEBP1 was related to immune response, cell adhesion, apoptotic process, inflammatory response, positive regulation of cell proliferation, angiogenesis, response to drug, and response to hypoxia. The survival analysis showed that the overexpressed level of AEBP1 was correlated with short survival time in both glioma and GBM patients.

Conclusion

In summary, AEBP1 was related with GSC-induced TMZ resistance. Our study showed that AEBP1 might be an oncogene and a new effective therapeutic target for the treatment of glioma.

A comprehensive prognostic signature for glioblastoma patients based on transcriptomics and single cell sequencing

PURPOSE

Glioblastoma (GBM) is the most common and deadly brain tumor. We aimed to reveal potential prognostic GBM marker genes, elaborate their functions, and build an effective a prognostic model for GBM patients.

METHODS

Through data mining of The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA), we screened for significantly differentially expressed genes (DEGs) to calculate risk scores for individual patients. Published data of somatic mutation and copy number variation profiles were analyzed for distinct genomic alterations associated with risk scores. In addition, single-cell sequencing was used to explore the biological functions of the identified prognostic marker genes. By combining risk scores and other clinical features, we built a comprehensive prognostic GBM model.

RESULTS

Seven DEGs (CLEC5A, HOXC6, HOXA5, CCL2, GPRASP1, BSCL2 and PTX3) were identified as being prognostic for GBM. Expression of these genes was confirmed in different GBM cell lines using real-time PCR. Risk scores calculated from the seven DEGs revealed prognostic value irrespective of other clinical factors, including IDH mutation status, and were negatively correlated with TP53 expression. The prognostic genes were found to be associated with tumor proliferation and progression based on pseudo-time analysis in neoplastic cells. A final prognostic model was developed and validated with a good performance, especially in geriatric GBM patients.

CONCLUSIONS

Using genetic profiles, age, IDH mutation status, and chemotherapy and radiotherapy, we constructed a comprehensive prognostic model for GBM patients. The model has a good performance, especially in geriatric GBM patients.

The Predictive Value of Monocytes in Immune Microenvironment and Prognosis of Glioma Patients Based on Machine Learning

Gliomas are primary malignant brain tumors. Monocytes have been proved to actively participate in tumor growth. Weighted gene co-expression network analysis was used to identify meaningful monocyte-related genes for clustering. Neural network and SVM were applied for validating clustering results. Somatic mutation and copy number variation were used for defining the features of identified clusters. Differentially expressed genes (DEGs) between the stratified groups after performing elastic regression and principal component analyses were used for the construction of risk scores. Monocytes were associated with glioma patients’ survival and exhibited high predictive value. The prognostic value of risk score in glioma was validated by the abundant expression of immune checkpoint and metabolic profile. Additionally, high risk score was positively associated with the expression of immunogenic and antigen presenting factors, which indicated high immune infiltration. A prognostic model based on risk score demonstrated high accuracy rate of receiver operating characteristic curves. Compared with previous studies, our research dissected functional roles of monocytes from large-scale analysis. Findings of our analyses strongly support an immune modulatory and prognostic role of monocytes in glioma progression. Notably, monocyte could be an effective predictor for therapy responses of glioma patients.