Pyroptosis in glioblastoma: A crucial regulator of the tumour immune microenvironment and a predictor of prognosis

Recent studies have shown that pyroptosis, an inflammatory form of cell death, has a dual role in tumorigenesis and tumour progression and affects the prognosis of patients; however, the role of pyroptosis in glioblastoma (GBM) is still unclear. In this study, based on GBM patients'' data from two independent cohorts, we performed a comprehensive analysis of the expression and prognostic value of 33 pyroptosis‐associated genes (PAGs) in GBM, as well as their role in the tumour immune microenvironment (TIME) of GBM. We identified 29 PAGs that were differentially expressed between GBM and normal brain tissue, 18 of which were upregulated in GBM tissue. Most of the 33 PAGs were strongly correlated with the levels of immune cell infiltration. Based on the 33 PAGs, the GBM samples can be divided into two clusters (C1‐C2), with C1 having a ‘hot’ but immunosuppressive TIME and C2 having a ‘cold’ TIME, suggesting different immunotherapeutic responses in the two clusters. In addition, we identified four PAGs that were strongly associated with GBM prognosis and constructed a risk model based on these four PAGs. This risk model is an independent prognostic factor for GBM patients, and there is a different immune status between high‐ and low‐risk groups. In conclusion, this study demonstrates that pyroptosis is closely associated with the prognosis and TIME of GBM and provides an important basis for further studies on the relationship between pyroptosis and GBM.     

Predicting the prognosis of glioma by pyroptosis‐related signature

Glioma is the most common malignant primary brain tumour. It is of great significance for the prognosis and personalized treatment of glioma patients to accurate identification of glioma based on biomarkers. Pyroptosis, a kind of programmed cell death, is closely related to tumour progression and tumour immune microenvironment. However, the role of pyroptosis in glioma remained unclear. Herein, we used glioma clinical and expression data from TCGA and CGGA to explore the relationship between pyroptosis and glioma. We first summarized the incidence of copy number variations and somatic mutations of 33 pyroptosis‐related genes and explored prognostic correlation of these genes. Based on pyroptosis‐related genes, three molecular subgroups of glioma related to prognosis were identified. We also found that each subgroup has unique immune and biological behaviours characteristics. Finally, based on 7 pyroptosis‐related genes (CASP3, CASP4, CASP6, CASP8, CASP9, PRKACA and ELANE), we constructed a prognosis model by Lasso and Cox regression, which had a strong predictive power for the overall survival in CGGA test cohort (p < 0.05). In summary, we explored the role of pyroptosis‐related genes in gliomas and the association of these genes with tumour immunity. We found the biomarkers valuable to diagnosis and prognosis, hence, provide reference to the development and treatment of tumorigenesis in glioma.     

Pyroptosis‐related molecular classification and immune microenvironment infiltration in breast cancer: A novel therapeutic target

The underlying role of pyroptosis in breast cancer (BC) remains unknown. Herein, we investigated the correlations of 33 pyroptosis‐related genes (PRGs) with immune checkpoints and immune cell infiltrations in BC patients based on The Cancer Genome Atlas cohort (n = 996) and Gene Expression Omnibus cohort (n = 3,262). Enrichment analysis revealed that these PRGs mainly functioned in pyroptosis, inflammasomes and regulation of autophagy pathway. Four prognostic independent PRGs (CASP9, TIRAP, GSDMC and IL18) were identified. Then, cluster 1/2 was recognized using consensus clustering for these four PRGs. Patients from cluster 1 had a favourable prognosis and diverse immune cell infiltrations. A nomogram was developed based on age, TNM stage, tumour subtype and pyroptosis score. Patients with the high‐risk group exhibited worse 5‐year OS, and the result was consistent in the external cohort. Additionally, high‐risk group patients were associated with downregulated immune checkpoint expression. Further analysis suggested that the high‐risk group patients were associated with a higher IC50 of paclitaxel, doxorubicin, cisplatin, methotrexate and vinorelbine. In summarizing, the pyroptosis score‐based nomogram might serve as an independent prognostic predictor and could guide medication for chemotherapy. Additionally, it may bring novel insight into the regulation of tumour immune microenvironment in BC and help to achieve precision immunotherapy.     

Circadian clock genes promote glioma progression by affecting tumour immune infiltration and tumour cell proliferation

ObjectivesCircadian rhythm controls complicated physiological activities in organisms. Circadian clock genes have been related to tumour progression, but its role in glioma is unknown. Therefore, we explored the relationship between dysregulated circadian clock genes and glioma progression.Materials and MethodsSamples were divided into different groups based on circadian clock gene expression in training dataset (n = 672) and we verified the results in other four validating datasets (n = 1570). The GO and GSEA enrichment analysis were conducted to explore potential mechanism of how circadian clock genes affected glioma progression. The single‐cell RNA‐Seq analysis was conducted to verified previous results. The immune landscape was evaluated by the ssGSEA and CIBERSORT algorithm. Cell proliferation and viability were confirmed by the CCK8 assay, colony‐forming assay and flow cytometry.ResultsThe cluster and risk model based on circadian clock gene expression can predict survival outcome. Samples were scoring by the least absolute shrinkage and selection operator regression analysis, and high scoring tumour was associated with worse survival outcome. Samples in high‐risk group manifested higher activation of immune pathway and cell cycle. Tumour immune landscape suggested high‐risk tumour infiltrated more immunocytes and more sensitivity to immunotherapy. Interfering TIMELESS expression affected circadian clock gene expression, inhibited tumour cell proliferation and arrested cell cycle at the G0/G1 phase.ConclusionsDysregulated circadian clock gene expression can affect glioma progression by affecting tumour immune landscape and cell cycle. The risk model can predict glioma survival outcome, and this model can also be applied to pan‐cancer.     

Prediction of a competing endogenous RNA co‐expression network as a prognostic marker in glioblastoma

Due to its high proliferation capacity and rapid intracranial spread, glioblastoma (GBM) has become one of the least curable malignant cancers. Recently, the competing endogenous RNAs (ceRNAs) hypothesis has become a focus in the researches of molecular biological mechanisms of cancer occurrence and progression. However, there is a lack of correlation studies on GBM, as well as a lack of comprehensive analyses of GBM molecular mechanisms based on high‐throughput sequencing and large‐scale sample sizes. We obtained RNA‐seq data from The Cancer Genome Atlas (TCGA) and Genotype‐Tissue Expression (GTEx) databases. Further, differentially expressed mRNAs were identified from normal brain tissue and GBM tissue. The similarities between the mRNA modules with clinical traits were subjected to weighted correlation network analysis (WGCNA). With the mRNAs from clinical‐related modules, a survival model was constructed by univariate and multivariate Cox proportional hazard regression analyses. Thereafter, we carried out Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Finally, we predicted interactions between lncRNAs, miRNAs and mRNAs by TargetScan, miRDB, miRTarBase and starBase. We identified 2 lncRNAs (NORAD, XIST), 5 miRNAs (hsa‐miR‐3613, hsa‐miR‐371, hsa‐miR‐373, hsa‐miR‐32, hsa‐miR‐92) and 2 mRNAs (LYZ, PIK3AP1) for the construction of a ceRNA network, which might act as a prognostic biomarker of GBM. Combined with previous studies and our enrichment analysis results, we hypothesized that this ceRNA network affects immune activities and tumour microenvironment variations. Our research provides novel aspects to study GBM development and treatment.     

RNA‐binding protein RBM47 stabilizes IFNAR1 mRNA to potentiate host antiviral activity

The type I interferon (IFN‐I, IFN‐α/β)‐mediated immune response is the first line of host defense against invading viruses. IFN‐α/β binds to IFN‐α/β receptors (IFNARs) and triggers the expression of IFN‐stimulated genes (ISGs). Thus, stabilization of IFNARs is important for prolonging antiviral activity. Here, we report the induction of an RNA‐binding motif‐containing protein, RBM47, upon viral infection or interferon stimulation. Using multiple virus infection models, we demonstrate that RBM47 has broad‐spectrum antiviral activity in vitro and in vivo. RBM47 has no noticeable impact on IFN production, but significantly activates the IFN‐stimulated response element (ISRE) and enhances the expression of interferon‐stimulated genes (ISGs). Mechanistically, RBM47 binds to the 3''UTR of IFNAR1 mRNA, increases mRNA stability, and retards the degradation of IFNAR1. In summary, this study suggests that RBM47 is an interferon‐inducible RNA‐binding protein that plays an essential role in enhancing host IFN downstream signaling.     

HNF1B inhibits cell proliferation via repression of SMAD6 expression in prostate cancer

Prostate cancer is the most common malignancy in men in developed countries. In previous study, we identified HNF1B (Hepatocyte Nuclear Factor 1β) as a downstream effector of Enhancer of zeste homolog 2 (EZH2). HNF1B suppresses EZH2‐mediated migration of two prostate cancer cell lines via represses the EMT process by inhibiting SLUG expression. Besides, HNF1B expression inhibits cell proliferation through unknown mechanisms. Here, we demonstrated that HNF1B inhibited the proliferation rate of prostate cancer cells. Overexpression of HNF1B in prostate cancer cells led to the arrest of G1 cell cycle and decreased Cyclin D1 expression. In addition, we re‐explored data from ChIP‐sequencing (ChIP‐seq) and RNA‐sequencing (RNA‐seq), and demonstrated that HNF1B repressed Cyclin D1 via direct suppression of SMAD6 expression. We also identified CDKN2A as a HNF1B‐interacting protein that would contribute to HNF1B‐mediated repression of SMAD6 expression. In summary, we provide the novel mechanisms and evidence in support HNF1B as a tumour suppressor gene for prostate cancer.     

Modulation of lncRNA H19 enhances resveratrol‐inhibited cancer cell proliferation and migration by regulating endoplasmic reticulum stress

The phytoalexin resveratrol exhibits anti‐tumour activity in many types of cancer. In this study, we showed that resveratrol suppressed the survival of gastric tumour cells both in vivo and in vitro. Resveratrol promoted apoptosis, autophagy and endoplasmic reticulum (ER) stress in a dose‐dependent manner. RNA‐seq analysis showed that multiple cell death signalling pathways were activated after resveratrol treatment, while the use of ER stress activators (tunicamycin and thapsigargin) in combinatorial with resveratrol led to further inhibition of cancer cell survival. Results also showed that resveratrol altered the expression of several long non‐coding RNAs (lncRNAs), including MEG3, PTTG3P, GAS5, BISPR, MALAT1 and H19. Knockdown of H19 in resveratrol‐treated cells further enhanced the effects of resveratrol on apoptosis, ER stress and cell cycle S‐phase arrest. Furthermore, the migratory ability of resveratrol‐treated cells was dramatically decreased after H19 knockdown. In conclusion, resveratrol inhibited cancer cell survival, while knockdown of lncRNA H19 resulted in increased sensitivity to resveratrol therapy.     

Characterizing heterogeneity of non‐small cell lung tumour microenvironment to identify signature prognostic genes

Growing evidence has highlighted the immune response as an important feature of carcinogenesis and therapeutic efficacy in non‐small cell lung cancer (NSCLC). This study focused on the characterization of immune infiltration profiling in patients with NSCLC and its correlation with survival outcome. All TCGA samples were divided into three heterogeneous clusters based on immune cell profiles: cluster 1 (''low infiltration'' cluster), cluster 2 (''heterogeneous infiltration'' cluster) and cluster 3 (''high infiltration'' cluster). The immune cells were responsible for a significantly favourable prognosis for the ''high infiltration'' community. Cluster 1 had the lowest cytotoxic activity, tumour‐infiltrating lymphocytes and interferon‐gamma (IFN‐γ), as well as immune checkpoint molecules expressions. In addition, MHC‐I and immune co‐stimulator were also found to have lower cluster 1 expressions, indicating a possible immune escape mechanism. A total of 43 differentially expressed genes (DEGs) that overlapped among the groups were determined based on three clusters. Finally, based on a univariate Cox regression model, prognostic immune‐related genes were identified and combined to construct a risk score model able to predict overall survival (OS) rates in the validation datasets.     

Novel long non‐coding RNA CYB561‐5 promotes aerobic glycolysis and tumorigenesis by interacting with basigin in non‐small cell lung cancer

Abnormally expressed long non‐coding RNAs (lncRNAs) have been recognized as potential diagnostic biomarkers or therapeutic targets in non‐small cell lung cancer (NSCLC). The role of the novel lnc‐CYB561‐5 in NSCLC and its specific biological activity remain unknown. In this study, lncRNAs highly expressed in NSCLC tissue samples compared with paired adjacent normal tissue samples and atypical adenomatous hyperplasia were identified by RNA‐seq analysis. Lnc‐CYB561‐5 is highly expressed in human NSCLC and is associated with a poor prognosis in lung adenocarcinoma. In vivo, downregulation of lnc‐CYB561‐5 significantly decreases tumour growth and metastasis. In vitro, lnc‐CYB561‐5 knockdown treatment inhibits cell migration, invasion and proliferation ability, as well as glycolysis rates. In addition, RNA pulldown and RNA immunoprecipitation (RIP) assays show that basigin (Bsg) protein interacts with lnc‐CYB561‐5. Overall, this study demonstrates that lnc‐CYB561‐5 is an oncogene in NSCLC, which is involved in the regulation of cell proliferation and metastasis. Lnc‐CYB561‐5 interacts with Bsg to promote the expression of Hk2 and Pfk1 and further lead to metabolic reprogramming of NSCLC cells.     

Prognostic value and immune cell infiltration of hypoxic phenotype‐related gene signatures in glioblastoma microenvironment

Glioblastoma (GBM) is a malignant intracranial tumour with the highest proportion and lethality. It is characterized by invasiveness and heterogeneity. However, the currently available therapies are not curative. As an essential environmental cue that maintains glioma stem cells, hypoxia is considered the cause of tumour resistance to chemotherapy and radiation. Growing evidence shows that immunotherapy focusing on the tumour microenvironment is an effective treatment for GBM; however, the current clinicopathological features cannot predict the response to immunotherapy and provide accurate guidance for immunotherapy. Based on the ESTIMATE algorithm, GBM cases of The Cancer Genome Atlas (TCGA) data set were classified into high‐ and low‐immune/stromal score groups, and a four‐gene tumour environment‐related model was constructed. This model exhibited good efficiency at forecasting short‐ and long‐term prognosis and could also act as an independent prognostic biomarker. Additionally, this model and four of its genes (CLECL5A, SERPING1, CHI3L1 and C1R) were found to be associated with immune cell infiltration, and further study demonstrated that these four genes might drive the hypoxic phenotype of perinecrotic GBM, which affects hypoxia‐induced glioma stemness. Therefore, these might be important candidates for immunotherapy of GBM and deserve further exploration.     

Ginsenoside Rh2 mitigates doxorubicin‐induced cardiotoxicity by inhibiting apoptotic and inflammatory damage and weakening pathological remodelling in breast cancer‐bearing mice

ObjectivesThere are presently a few viable ways to reduce cardiotoxicity of doxorubicin (Dox). The combination of chemotherapy agents with natural compounds delivers greater efficacy and reduces adverse effects in recent researches for cancer treatment. Here, we examined the potential effect of ginsenoside Rh2 on a Dox‐based regimen in chemotherapy treatment.Materials and MethodsHuman breast tumour (MDA‐MB‐231) xenograft nude mice, human cardiac ventricle fibroblasts, and human umbilical vein endothelial cells (HUVEC) were employed in the present study. Histology, immunohistochemistry, immunofluorescence, western blot, antibody array, and RNA‐sequencing analyses were utilized to assess the protective effect of Rh2 on cardiotoxicity induced by Dox and the underlying mechanisms.ResultsRh2‐reduced cardiotoxicity by inhibiting the cardiac histopathological changes, apoptosis and necrosis, and consequent inflammation. Pathological remodelling was attenuated by reducing fibroblast to myofibroblast transition (FMT) and endothelial–mesenchymal transition (EndMT) in hearts. RNA‐sequencing analysis showed that Dox treatment predominantly targets cell cycle and attachment of microtubules and boosted tumour necrosis, chemokine and interferon‐gamma production, response to cytokine and chemokine, and T cell activation, whereas Rh2 regulated these effects. Intriguingly, Rh2 also attenuated fibrosis via promoting senescence in myofibroblasts and reversing established myofibroblast differentiation in EndMT.ConclusionsRh2 regulates multiple pathways in the Dox‐provoked heart, proposing a potential candidate for cancer supplement and therapy‐associated cardiotoxicity.

Doxorubicin is extensively reported to induce severe cardiotoxicity in clinical applications. Our work proposed a natural herbal compound, ginsenoside Rh2, as a potential candidate for attenuating this side effect. Rh2 significantly inhibited cardiac apoptosis and necrosis, inflammation, and pathological remodelling in Dox‐challenged hearts.     

TATA and paused promoters active in differentiated tissues have distinct expression characteristics

Core promoter types differ in the extent to which RNA polymerase II (Pol II) pauses after initiation, but how this affects their tissue‐specific gene expression characteristics is not well understood. While promoters with Pol II pausing elements are active throughout development, TATA promoters are highly active in differentiated tissues. We therefore used a genomics approach on late‐stage Drosophila embryos to analyze the properties of promoter types. Using tissue‐specific Pol II ChIP‐seq, we found that paused promoters have high levels of paused Pol II throughout the embryo, even in tissues where the gene is not expressed, while TATA promoters only show Pol II occupancy when the gene is active. The promoter types are associated with different chromatin accessibility in ATAC‐seq data and have different expression characteristics in single‐cell RNA‐seq data. The two promoter types may therefore be optimized for different properties: paused promoters show more consistent expression when active, while TATA promoters have lower background expression when inactive. We propose that tissue‐specific genes have evolved to use two different strategies for their differential expression across tissues.