Super-Enhancer Associated Five-Gene Risk Score Model Predicts Overall Survival in Multiple Myeloma Patients

Construction and experimental validation of a signature for predicting prognosis and immune infiltration analysis of glioma based on disulfidptosis-related lncRNAs

Backgrounds

Disulfidptosis, a newly discovered mechanism of programmed cell death, is believed to have a unique role in elucidating cancer progression and guiding cancer therapy strategies. However, no studies have yet explored this mechanism in glioma.

Methods

We downloaded data on glioma patients from online databases to address this gap. Subsequently, we identified disulfidptosis-related genes from published literature and verified the associated lncRNAs.

Results

Through univariate, multivariate, and least absolute shrinkage and selection operator (LASSO) regression algorithms analyses, we identified 10 lncRNAs. These were then utilized to construct prognostic prediction models, culminating in a risk-scoring signature. Reliability and validity tests demonstrated that the model effectively discerns glioma patients' prognosis outcomes. We also analyzed the relationship between the risk score and immune characteristics, and identified several drugs that may be effective for high-risk patients. In vitro experiments revealed that LINC02525 could enhances glioma cells' migration and invasion capacities. Additionally, knocking down LINC02525 was observed to promote glioma cell disulfidptosis.

Conclusion

This study delves into disulfidptosis-related lncRNAs in glioma, offering novel insights into glioma therapeutic strategies.

Complete miRNA-15/16 loss in mice promotes hematopoietic progenitor expansion and a myeloid-biased hyperproliferative state

Dysregulated apoptosis and proliferation are fundamental properties of cancer, and microRNAs (miRNA) are critical regulators of these processes. Loss of miR-15a/16-1 at chromosome 13q14 is the most common genomic aberration in chronic lymphocytic leukemia (CLL). Correspondingly, the deletion of either murine miR-15a/16-1 or miR-15b/16-2 locus in mice is linked to B cell lymphoproliferative malignancies. However, unexpectedly, when both miR-15/16 clusters are eliminated, most double knockout (DKO) mice develop acute myeloid leukemia (AML). Moreover, in patients with CLL, significantly reduced expression of miR-15a, miR-15b, and miR-16 associates with progression of myelodysplastic syndrome to AML, as well as blast crisis in chronic myeloid leukemia. Thus, the miR-15/16 clusters have a biological relevance for myeloid neoplasms. Here, we demonstrate that the myeloproliferative phenotype in DKO mice correlates with an increase of hematopoietic stem and progenitor cells (HSPC) early in life. Using single-cell transcriptomic analyses, we presented the molecular underpinning of increased myeloid output in the HSPC of DKO mice with gene signatures suggestive of dysregulated hematopoiesis, metabolic activities, and cell cycle stages. Functionally, we found that multipotent progenitors (MPP) of DKO mice have increased self-renewing capacities and give rise to significantly more progeny in the granulocytic compartment. Moreover, a unique transcriptomic signature of DKO MPP correlates with poor outcome in patients with AML. Together, these data point to a unique regulatory role for miR-15/16 during the early stages of hematopoiesis and to a potentially useful biomarker for the pathogenesis of myeloid neoplasms.

A novel signature constructed by super-enhancer-related genes for the prediction of prognosis in hepatocellular carcinoma and associated with immune infiltration

Background

Super-enhancer (SE) refers to a regulatory element with super transcriptional activity, which can enrich transcription factors and drive gene expression. SE-related genes play an important role in the pathogenesis of malignant tumors, including hepatocellular carcinoma (HCC).

Methods

The SE-related genes were obtained from the human super-enhancer database (SEdb). Data from the transcriptome analysis and related clinical information with HCC were obtained from The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) database. The upregulated SE-related genes from TCGA-LIHC were identified by the DESeq2R package. Multivariate Cox regression analysis was used to construct a four-gene prognostic signature. According to the median risk score, HCC patients were divided into high-risk and low-risk group patients.

Results

The Kaplan-Meier (KM) curve showed that a significantly worse prognosis was found for the high-risk group (P<0.001). In the TCGA-LIHC dataset, the area under the curve (AUC) values were 0.737, 0.662, and 0.667 for the model predicting overall survival (OS) over 1-, 3-, and 5- years, respectively, indicating the good prediction ability of our prediction model. This model's prognostic value was further validated in the LIRI-JP dataset and HCC samples (n=65). Furthermore, we found that higher infiltration level of M0 macrophages and upregulated of CTLA4 and PD1 in the high-risk group, implying that immunotherapy could be effective for those patients.

Conclusion

These results provide further evidence that the unique SE-related gene model could accurately predict the prognosis of HCC.

Superenhancers as master gene regulators and novel therapeutic targets in brain tumors

Transcriptional deregulation, a cancer cell hallmark, is driven by epigenetic abnormalities in the majority of brain tumors, including adult glioblastoma and pediatric brain tumors. Epigenetic abnormalities can activate epigenetic regulatory elements to regulate the expression of oncogenes. Superenhancers (SEs), identified as novel epigenetic regulatory elements, are clusters of enhancers with cell-type specificity that can drive the aberrant transcription of oncogenes and promote tumor initiation and progression. As gene regulators, SEs are involved in tumorigenesis in a variety of tumors, including brain tumors. SEs are susceptible to inhibition by their key components, such as bromodomain protein 4 and cyclin-dependent kinase 7, providing new opportunities for antitumor therapy. In this review, we summarized the characteristics and identification, unique organizational structures, and activation mechanisms of SEs in tumors, as well as the clinical applications related to SEs in tumor therapy and prognostication. Based on a review of the literature, we discussed the relationship between SEs and different brain tumors and potential therapeutic targets, focusing on glioblastoma.

Hypoxia–Immune-Related Gene SLC19A1 Serves as a Potential Biomarker for Prognosis in Multiple Myeloma

Background

Multiple myeloma (MM) remains an incurable malignant tumor of plasma cells. Increasing evidence has reported that hypoxia and immune status contribute to the progression of MM. In this research, the prognostic value of the hypoxia–immune-related gene SLC19A1 in MM was evaluated by bioinformatics analysis.

Method

RNA-sequencing (RNA-seq) data along with clinical information on MM were downloaded from the Gene Expression Omnibus (GEO) database. Consistent clustering analysis and ESTIMATE algorithms were performed to establish the MM sample subgroups related to hypoxia and immune status, respectively, based on the GSE24080 dataset. The differentially expressed analysis was performed to identify the hypoxia–immune-related genes. Subsequently, a hypoxia–immune-gene risk signature for MM patients was constructed by univariate and multivariate Cox regression analyses, which was also verified in the GSE4581 dataset. Furthermore, the mRNA expression of SLC19A1 was determined using qRT-PCR in 19 MM patients, and the correlations between the genetic expression of SLC19A1 and clinical features were further analyzed.

Result

A total of 47 genes were identified as hypoxia–immune-related genes for MM. Among these genes, SLC19A1 was screened to construct a risk score model that had better predictive power for MM. The constructed prognostic signature based on SLC19A1 was verified in the GSE4581 dataset. All independent prognostic factors (age, β₂-microglobulin, LDH, albumin, MRI, and gene risk score) were used to develop a nomogram that showed a better performance for predicting the survival probability of MM patients for 1–5 years. Furthermore, SLC19A1 was highly expressed in newly diagnosed and relapsed MM patients, and high expression of SLC19A1 is correlated with higher bone marrow aspiration plasma cells and β₂-microglobulin levels in MM patients.

Conclusion

In conclusion, our results suggest that SLC19A1 is aberrantly expressed in MM and highly expressed SLC19A1 might be a biomarker correlated with inferior prognosis. More importantly, we identified SLC19A1 as a hypoxia–immune-related gene in MM. Future functional and mechanistic studies will further clarify the roles of SLC19A1 in MM.

Clinical M2 macrophages-related genes to aid therapy in pancreatic ductal adenocarcinoma

Background

Increasing evidence supports that infiltration M2 Macrophages act as pivotal player in tumor progression of pancreatic ductal adenocarcinoma (PDAC). Nonetheless, comprehensive analysis of M2 Macrophage infiltration and biological roles of hub genes (FAM53B) in clinical outcome and immunotherapy was lack.

Method

The multiomic data of PDAC samples were downloaded from distinct datasets. CIBERSORT algorithm was performed to uncover the landscape of TIME. Weighted gene co-expression network analysis (WGCNA) was performed to identify candidate module and significant genes associated with M2 Macrophages. Kaplan-Meier curve and receiver operating characteristic (ROC) curves were applied for prognosis value validation. Mutation data was analyzed by using “maftools” R package. Gene set variation analysis (GSVA) was employed to assign pathway activity estimates to individual sample. Immunophenoscore (IPS) was implemented to estimate immunotherapeutic significance of risk score. The half-maximal inhibitory concentration (IC50) of chemotherapeutic drugs was predicted by using the pRRophetic algorithm. Finally, quantitative real-time polymerase chain reaction was used to determine FAM53B mRNA expression and TIMER database was utilized to uncover its possible role in immune infiltration of PDAC.

Results

Herein, 17,932 genes in 234 samples (214 tumor and 20 normal) were extracted from three platforms. Taking advantage of WGCNA, significant module (royalblue) and 135 candidate genes were considered as M2 Macrophages-related genes. Subsequently, risk signature including 5 hub genes was developed by multiple analysis, which exhibited excellent prognostic performance. Besides, comprehensive prognostic nomogram was constructed to quantitatively estimate risk. Then, intrinsic link between risk score with tumor mutation burden (TMB) was explored. Additionally, risk score significantly correlated with diversity of tumor immune microenvironment (TIME). PDAC samples within different risk presented diverse signaling pathways activity and experienced significantly distinct sensitivity to administering chemotherapeutic or immunotherapeutic agents. Finally, the biological roles of FAM53B were revealed in PDAC.

Conclusions

Taken together, comprehensive analyses of M2 Macrophages profiling will facilitate prognostic prediction, delineating complexity of TIME, and contribute insight into precision therapy for PDAC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02289-w.

An interactive nomogram based on clinical and molecular signatures to predict prognosis in multiple myeloma patients

Although novel drugs and treatments have been developed and improved, multiple myeloma (MM) is still recurrent and difficult to cure. In the present study, the magenta module containing 400 hub genes was determined from the training dataset of GSE24080 through weighted gene co-expression network analysis (WGCNA). Then, using the least absolute shrinkage and selection operator (Lasso) analysis, a fifteen-gene signature was firstly selected and the predictive performance for overall survival (OS) was favorable, which was identified by Receiver Operating Characteristic (ROC) curves. The risk score model was constructed based on survival-associated fifteen genes from the Lasso model, which classified MM patients into high-risk and low-risk groups. Areas under the curve (AUC) of ROC curve and log-rank test showed that the high-risk group was correlated to the dismal survival outcome of MM patients, which was also identified in testing dataset of GSE9782. The calibration plot, the AUC value of the ROC curve and Concordance-index showed that the interactive nomogram with risk score could favorably predict the probability of multi-year OS of MM patients. Therefore, it may help clinicians make a precise therapeutic decision based on the easy-to-use tool of the nomogram.

Identification of a three-gene-based prognostic model in multiple myeloma using bioinformatics analysis

Background

Multiple myeloma (MM), the second most hematological malignancy, has high incidence and remains incurable till now. The pathogenesis of MM is poorly understood. This study aimed to identify novel prognostic model for MM on gene expression profiles.

Methods

Gene expression datas of MM (GSE6477, GSE136337) were downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) in GSE6477 between case samples and normal control samples were screened by the limma package. Meanwhile, enrichment analysis was conducted, and a protein-protein interaction (PPI) network of these DEGs was established by STRING and cytoscape software. Co-expression modules of genes were built by Weighted Correlation Network Analysis (WGCNA). Key genes were identified both from hub genes and the DEGs. Univariate and multivariate Cox congression were performed to screen independent prognostic genes to construct a predictive model. The predictive power of the model was evaluated by Kaplan–Meier curve and time-dependent receiver operating characteristic (ROC) curves. Finally, univariate and multivariate Cox regression analyse were used to investigate whether the prognostic model could be independent of other clinical parameters.

Results

GSE6477, including 101 case and 15 normal control, were screened as the datasets. A total of 178 DEGs were identified, including 59 up-regulated and 119 down-regulated genes. In WGCNA analysis, module black and module purple were the most relevant modules with cancer traits, and 92 hub genes in these two modules were selected for further analysis. Next, 47 genes were chosen both from the DEGs and hub genes as key genes. Three genes (LYVE1, RNASE1, and RNASE2) were finally screened by univariate and multivariate Cox regression analyses and used to construct a risk model. In addition, the three-gene prognostic model revealed independent and accurate prognostic capacity in relation to other clinical parameters for MM patients.

Conclusion

In summary, we identified and constructed a three-gene-based prognostic model that could be used to predict overall survival of MM patients.

Super-enhancers: a new frontier for epigenetic modifiers in cancer chemoresistance

Although new developments of surgery, chemotherapy, radiotherapy, and immunotherapy treatments for cancer have improved patient survival, the emergence of chemoresistance in cancer has significant impacts on treatment effects. The development of chemoresistance involves several polygenic, progressive mechanisms at the molecular and cellular levels, as well as both genetic and epigenetic heterogeneities. Chemotherapeutics induce epigenetic reprogramming in cancer cells, converting a transient transcriptional state into a stably resistant one. Super-enhancers (SEs) are central to the maintenance of identity of cancer cells and promote SE-driven-oncogenic transcriptions to which cancer cells become highly addicted. This dependence on SE-driven transcription to maintain chemoresistance offers an Achilles' heel for chemoresistance. Indeed, the inhibition of SE components dampens oncogenic transcription and inhibits tumor growth to ultimately achieve combined sensitization and reverse the effects of drug resistance. No reviews have been published on SE-related mechanisms in the cancer chemoresistance. In this review, we investigated the structure, function, and regulation of chemoresistance-related SEs and their contributions to the chemotherapy via regulation of the formation of cancer stem cells, cellular plasticity, the microenvironment, genes associated with chemoresistance, noncoding RNAs, and tumor immunity. The discovery of these mechanisms may aid in the development of new drugs to improve the sensitivity and specificity of cancer cells to chemotherapy drugs.