Gene Network Analysis of Hepatocellular Carcinoma Identifies Modules Associated with Disease Progression, Survival, and Chemo Drug Resistance

Precision models in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) represents a global health challenge, and ranks among one of the most prevalent and deadliest cancers worldwide. Therapeutic advances have expanded the treatment armamentarium for patients with advanced HCC, but obstacles remain. Precision oncology, which aims to match specific therapies to patients who have tumours with particular features, holds great promise. However, its implementation has been hindered by the existence of numerous 'HCC influencers' that contribute to the high inter-patient heterogeneity. HCC influencers include tumour-related characteristics, such as genetic alterations, immune infiltration, stromal composition and aetiology, and patient-specific factors, such as sex, age, germline variants and the microbiome. This Review delves into the intricate world of HCC, describing the most innovative model systems that can be harnessed to identify precision and/or personalized therapies. We provide examples of how different models have been used to nominate candidate biomarkers, their limitations and strategies to optimize such models. We also highlight the importance of reproducing distinct HCC influencers in a flexible and modular way, with the aim of dissecting their relative contribution to therapy response. Next-generation HCC models will pave the way for faster discovery of precision therapies for patients with advanced HCC.

Cancer drug resistance as learning of signaling networks

Drug resistance is a major cause of tumor mortality. Signaling networks became useful tools for driving pharmacological interventions against cancer drug resistance. Signaling datasets now cover the entire human cell. Recently, network adaptation became understood as a learning process. We review rapidly increasing evidence showing that the development of cancer drug resistance can be described as learning of signaling networks. During drug adaptation, the network forgets drug-affected pathways by desensitization and relearns by strengthening alternative pathways. Thus, resistant cancer cells develop a drug resistance memory. We show that all key players of cellular learning (i.e., IDPs, protein translocation, microRNAs/lncRNAs, scaffolding proteins and epigenetic/chromatin memory) have important roles in the development of cancer drug resistance. Moreover, all of them are central components of the epithelial-mesenchymal transition leading to metastases and resistance. Phenotypic plasticity was recently listed as a hallmark of cancer. We review how network plasticity induces rare, pre-existent drug-resistant cells in the absence of drug treatment. Key network methods assessing the development of drug resistance and network pharmacological interventions against drug resistance are summarized. Finally, we highlight the class of cellular memory drugs affecting cellular learning and forgetting, and we summarize current challenges to prevent or break drug resistance using network models.

Identification of differentially expressed mRNA/lncRNA modules in acutely regorafenib-treated sorafenib-resistant Huh7 hepatocellular carcinoma cells

The multikinase inhibitor sorafenib is the standard first-line treatment for advanced hepatocellular carcinoma (HCC), but many patients become sorafenib-resistant (SR). This study investigated the efficacy of another kinase inhibitor, regorafenib (Rego), as a second-line treatment. We produced SR HCC cells, wherein the PI3K-Akt, TNF, cAMP, and TGF-beta signaling pathways were affected. Acute Rego treatment of these cells reversed the expression of genes involved in TGF-beta signaling but further increased the expression of genes involved in PI3K-Akt signaling. Additionally, Rego reversed the expression of genes involved in nucleosome assembly and epigenetic gene expression. Weighted gene co-expression network analysis (WGCNA) revealed four differentially expressed long non-coding RNA (DElncRNA) modules that were associated with the effectiveness of Rego on SR cells. Eleven putative DElncRNAs with distinct expression patterns were identified. We associated each module with DEmRNAs of the same pattern, thus obtaining DElncRNA/DEmRNA co-expression modules. We discuss the potential significance of each module. These findings provide insights and resources for further investigation into the potential mechanisms underlying the response of SR HCC cells to Rego.

A gene network database for the identification of key genes for diagnosis, prognosis, and treatment in sepsis

Sepsis and sepsis-related diseases cause a high rate of mortality worldwide. The molecular and cellular mechanisms of sepsis are still unclear. We aim to identify key genes in sepsis and reveal potential disease mechanisms. Six sepsis-related blood transcriptome datasets were collected and analyzed by weighted gene co-expression network analysis (WGCNA). Functional annotation was performed in the gProfiler tool. DSigDB was used for drug signature enrichment analysis. The proportion of immune cells was estimated by the CIBERSORT tool. The relationships between modules, immune cells, and survival were identified by correlation analysis and survival analysis. A total of 37 stable co-expressed gene modules were identified. These modules were associated with the critical biology process in sepsis. Four modules can independently separate patients with long and short survival. Three modules can recurrently separate sepsis and normal patients with high accuracy. Some modules can separate bacterial pneumonia, influenza pneumonia, mixed bacterial and influenza A pneumonia, and non-infective systemic inflammatory response syndrome (SIRS). Drug signature analysis identified drugs associated with sepsis, such as testosterone, phytoestrogens, ibuprofen, urea, dichlorvos, potassium persulfate, and vitamin B12. Finally, a gene co-expression network database was constructed ( https://liqs.shinyapps.io/sepsis/ ). The recurrent modules in sepsis may facilitate disease diagnosis, prognosis, and treatment.

Using Weighted Gene Co-Expression Network Analysis to Identify Increased MND1 Expression as a Predictor of Poor Breast Cancer Survival

Objective

We used bioinformatics analysis to identify potential biomarker genes and their relationship with breast cancer (BC).

Materials and Methods

We used a weighted gene co-expression network analysis (WGCNA) to create a co-expression network based on the top 25% genes in the GSE24124, GSE33926, and GSE86166 datasets obtained from the Gene Expression Omnibus. We used the DAVID online platform to perform GO and KEGG pathway enrichment analyses and the Cytoscape CytoHubba plug-in to screen the potential genes. Then, we related the genes to prognostic values in BC using the Oncomine, GEPIA, and Kaplan–Meier Plotter databases. Findings were validated by immunohistochemical (IHC) staining in the Human Protein Atlas and the TCGA-BRCA cohort. LinkedOmics identified the interactive expressions of hub genes. We used UALCAN to evaluate the methylation levels of these hub genes. MethSurv and SurvivalMeth were used to assess the multilevel prognostic value. Finally, we assessed hub gene association with immune cell infiltration using TIMER.

Results

The mRNA levels of MKI67, UBE2C, GTSE1, CCNA2, and MND1 were significantly upregulated in BC, whereas ESR1, THSD4, TFF1, AGR2, and FOXA1 were significantly downregulated. The DNA methylation signature analysis showed a better prognosis in the low-risk group. Further subgroup analyses revealed that MND1 might serve as an independent risk factor for unfavorable BC prognosis. Additionally, MND1 expression levels positively correlate with the immune infiltration statuses of CD4+ T cells, CD8+ T cells, B cells, neutrophils, dendritic cells, and macrophages.

Conclusion

Our results indicate that the ten hub genes may be involved in BC’s carcinogenesis, development, or metastasis, and MND1 may be a potential biomarker and therapeutic target for BC.