Identification and validation of potential key long noncoding RNAs in sorafenib-resistant hepatocellular carcinoma cells

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.

HMOX1 Attenuates the Sensitivity of Hepatocellular Carcinoma Cells to Sorafenib via Modulating the Expression of ABC Transporters

Hepatocellular carcinoma (HCC) represents a common malignancy, and mechanisms of acquired sorafenib resistance during the treatment of HCC patients remain elusive. The present study performed integrated bioinformatics analysis and explored the potential action of heme oxygenase 1 (HMOX1) in sorafenib-resistant HCC cells. Differentially expressed genes (DEGs) of the sorafenib-resistant group as compared to the sorafenib-sensitive group from GSE140202 and GSE143233 were extracted. Fifty common DEGs between GSE140202 and GSE143233 were extracted. Ten hub genes were identified from the protein-protein interaction network based on common DEGs. Experimental results revealed the upregulation of HMOX1 in sorafenib-resistant HCC cells. HMOX1 silence promoted the sensitivity to sorafenib in sorafenib-resistant HCC cells; overexpression of HMOX1 attenuated the sensitivity. In addition, HMOX1 silence downregulated the mRNA expression of ABC transporters in sorafenib-resistant HCC cells, while HMOX1 overexpression upregulated mRNA expression of ABC transporter expression in HCC cells. Further analysis also revealed that high expression of HMOX1 was associated with shorter OS and DSS in HCC patients. In conclusion, our analysis identified ten hub genes associated with sorafenib resistance in HCC. Further validation studies demonstrated that HMOX1 promoted sorafenib resistance of HCC cells via modulating ABC transporter expression.

Integrated Transcriptomic Analysis Revealed Hub Genes and Pathways Involved in Sorafenib Resistance in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC), a high mortality malignancy, has become a worldwide public health concern. Acquired resistance to the multikinase inhibitor sorafenib challenges its clinical efficacy and the survival benefits it provides to patients with advanced HCC. This study aimed to identify critical genes and pathways associated with sorafenib resistance in HCC using integrated bioinformatics analysis. Differentially expressed genes (DEGs) were identified using four HCC gene expression profiles (including 34 sorafenib-resistant and 29 sorafenib-sensitive samples) based on the robust rank aggregation method and R software. Gene ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) online tool. A protein–protein interaction (PPI) network was constructed using the Search Tool for the Retrieval of Interacting Genes (STRING), and small molecules reversing sorafenib resistance were searched for using the connectivity map (CMAP) database. Pearson correlation and survival analyses of hub genes were performed using cBioPortal and Gene Expression Profiling and Interactive Analysis (GEPIA). Finally, the expression levels of hub genes in sorafenib-resistant HCC cells were verified using quantitative polymerase chain reaction (q-PCR). A total of 165 integrated DEGs (66 upregulated and 99 downregulated in sorafenib resistant samples compared sorafenib sensitive ones) primarily enriched in negative regulation of endopeptidase activity, extracellular exosome, and protease binding were identified. Some pathways were commonly shared between the integrated DEGs. Seven promising therapeutic agents and 13 hub genes were identified. These findings provide a strategy and theoretical basis for overcoming sorafenib resistance in HCC patients.