Oncogenic Landscape of Somatic Mutations Perturbing Pan-Cancer lncRNA-ceRNA Regulation

LncRNAs: the good, the bad, and the unknown

Long non-coding RNAs (lncRNAs) are significant contributors in maintaining genomic integrity through epigenetic regulation. LncRNAs can interact with chromatin-modifying complexes in both cis and trans pathways, drawing them to specific genomic loci and influencing gene expression via DNA methylation, histone modifications, and chromatin remodeling. They can also operate as building blocks to assemble different chromatin-modifying components, facilitating their interactions and gene regulatory functions. Deregulation of these molecules has been associated with various human diseases, including cancer, cardiovascular disease, and neurological disorders. Thus, lncRNAs are implicated as potential diagnostic indicators and therapeutic targets. This review discusses the current understanding of how lncRNAs mediate epigenetic control, genomic integrity, and their putative functions in disease pathogenesis.

Multi-Omics Mining of lncRNAs with Biological and Clinical Relevance in Cancer

In this review, we provide a general overview of the current panorama of mining strategies for multi-omics data to investigate lncRNAs with an actual or potential role as biological markers in cancer. Several multi-omics studies focusing on lncRNAs have been performed in the past with varying scopes. Nevertheless, many questions remain regarding the pragmatic application of different molecular technologies and bioinformatics algorithms for mining multi-omics data. Here, we attempt to address some of the less discussed aspects of the practical applications using different study designs for incorporating bioinformatics and statistical analyses of multi-omics data. Finally, we discuss the potential improvements and new paradigms aimed at unraveling the role and utility of lncRNAs in cancer and their potential use as molecular markers for cancer diagnosis and outcome prediction.

A new pyroptosis-related signature for predicting the immune status and injury of liver ischemia-reperfusion

OBJECTIVE

Pyroptosis is a type of programmed cell death. This study aimed to explore the roles of key pyroptosis-related genes in liver ischemia-reperfusion injury.

METHODS

After collection and standardization of the transcriptome data from GSE12720 database, differentially expressed pyroptosis-related genes were identified. The risk genes screened by a random forest model were used to establish the line graph model. Consensus clustering was used to classify all samples according to the differentially expressed pyroptosis-related genes. Single-sample Gene Set Enrichment Analysis (ssGSEA) was performed to investigate the immune cell infiltration after hepatic ischemia-reperfusion. Cytoscape was used to visualize the regulatory network of transcription factor (TF)-microRNA (miRNA)-target genes.

RESULTS

We identified 18 significantly and differentially expressed pyroptosis-related genes between the disease and normal samples. Among these 18 genes, IL1β was positively correlated with CXCL8 (r = 0.791) and BIRC3 (r = 0.78), while ADORA3 was negatively correlated with GZMB (r = -0.567) and CXCL8 (r = -0.566). Furthermore, the random forest model constructed using the top 10 pyroptosis-related genes could predict the risk of hepatic ischemia-reperfusion. Importantly, the decision curve analysis showed that patients could benefit from the risk prediction model. Moreover, we found that the expression of TXNIP, IRF1, and GJA1 was the mostly regulated by miRNAs, while the expression of BIRC3, NFκB1, and TXNIP was regulated by the TF RELA. RELA had the most hub genes involved in the regulation.

CONCLUSION

Our study provides an overview of the expression landscape and the functional significance of pyroptosis-related genes in liver ischemia-reperfusion. Our findings also shed light on the clinical application of pyroptosis-related genes in the treatment of hepatic ischemia-reperfusion injury.

The Oncogenic and Tumor Suppressive Long Non-Coding RNA–microRNA–Messenger RNA Regulatory Axes Identified by Analyzing Multiple Platform Omics Data from Cr(VI)-Transformed Cells and Their Implications in Lung Cancer

Chronic exposure to hexavalent chromium (Cr(VI)) causes lung cancer in humans, however, the underlying mechanism has not been well understood. Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) are commonly studied non-coding RNAs. miRNAs function mainly through interaction with the 3′-untranslated regions of messenger RNAs (mRNAs) to down-regulate gene expression. LncRNAs have been shown to function as competing endogenous RNAs (ceRNAs) to sponge miRNAs and regulate gene expression. It is now well accepted that lncRNAs and miRNAs could function as oncogenes or tumor suppressors. Dysregulations of lncRNAs and miRNAs have been shown to play important roles in cancer initiation, progression, and prognosis. To explore the mechanism of Cr(VI) lung carcinogenesis, we performed lncRNA, mRNA, and miRNA microarray analysis using total RNAs from our previously established chronic Cr(VI) exposure malignantly transformed and passage-matched control human bronchial epithelial BEAS-2B cells. Based on the differentially expressed lncRNAs, miRNAs, and mRNAs between the control (BEAS-2B-Control) and Cr(VI)-transformed (BEAS-Cr(VI)) cells and by using the lncRNA–miRNA interaction and miRNA target prediction algorithms, we identified three oncogenic (HOTAIRM1/miR-182-5p/ERO1A, GOLGA8B/miR-30d-5p/RUNX2, and PDCD6IPP2/miR-23a-3p/HOXA1) and three tumor suppressive (ANXA2P1/miR-20b-5p/FAM241A (C4orf32), MIR99AHG/miR-218-5p/GPM6A, and SH3RF3-AS1/miR-34a-5p/HECW2) lncRNA–miRNA–mRNA regulatory axes. Moreover, the relevance of these three oncogenic and three tumor suppressive lncRNA–miRNA–mRNA regulatory axes in lung cancer was explored by analyzing publicly available human lung cancer omics datasets. It was found that the identified three oncogenic lncRNA–miRNA–mRNA regulatory axes (HOTAIRM1/miR-182-5p/ERO1A, GOLGA8B/miR-30d-5p/RUNX2, and PDCD6IPP2/miR-23a-3p/HOXA1) and the three tumor suppressive lncRNA–miRNA–mRNA regulatory axes (ANXA2P1/miR-20b-5p/FAM241A (C4orf32), MIR99AHG/miR-218-5p/GPM6A, and SH3RF3-AS1/miR-34a-5p/HECW2) have significant diagnostic and prognosis prediction values in human lung cancer. In addition, our recent studies showed that Cr(VI)-transformed cells display cancer stem cell (CSC)-like properties. Further bioinformatics analysis identified the oncogenic lncRNA–miRNA–mRNA regulatory axes as the potential regulators of cancer stemness. In summary, our comprehensive analysis of multiple platform omics datasets obtained from Cr(VI)-transformed human bronchial epithelial cells identified several oncogenic and tumor suppressive lncRNA–miRNA–mRNA regulatory axes, which may play important roles in Cr(VI) carcinogenesis and lung cancer in general.

Long noncoding RNA LIPH-4 promotes esophageal squamous cell carcinoma progression by regulating the miR-216b/IGF2BP2 axis

Introduction

Esophageal squamous cell carcinoma (ESCC) represents a major malignancy with poor clinical outcomes. Long noncoding RNAs (lncRNAs) are known to regulate the development and progression of multiple cancers. However, how lncRNAs are involved in ESCC is currently undefined.

Methods

LIPH-4 levels in ESCC tissue specimens and cells were assessed by qRT-PCR. The biological function of LIPH-4 was examined in cell and animal studies, applying CCK-8, EdU, colony formation and flow cytometry assays as well as xenograft model experiments. The underlying mechanisms of action of LIPH-4 were explored through bioinformatics, luciferase reporter assay, RNA-immunoprecipitation assay and immunoblot.

Results

We identified a novel lncRNA, LIPH-4, which showed elevated amounts in ESCC tissues and positive correlations with increased tumor size and poor prognosis in ESCC patients. Functional studies showed that LIPH-4 promoted the growth, mediated cell cycle progression and inhibited apoptosis in ESCC cells in vitro, and promoted tumor growth in mice. In terms of mechanism, LIPH-4 could bind to miR-216b and act as a competing endogenous RNA (ceRNA) to induce the expression of miR-216’s target gene IGF2BP2. LIPH-4 played an oncogenic role in ESCC through the miR-216b/IGF2BP2 axis.

Conclusions

This study suggested that LIPH-4 functions as a novel oncogenic lncRNA by acting as a ceRNA for miR-216b to regulate IGF2BP2, indicating LIPH-4 likely constitutes a prognostic biomarker and therapeutic target in ESCC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40364-022-00408-x.

LncRNA-IMAT1 Promotes Invasion of Meningiomas by Suppressing KLF4/hsa-miR22-3p/Snai1 Pathway

Malignant meningiomas often show invasive growth that makes complete tumor resection challenging, and they are more prone to recur after radical resection. Invasive meningioma associated transcript 1 (IMAT1) is a long noncoding RNA located on Homo sapiens chromosome 17 that was identified by our team based on absolute expression differences in invasive and non-invasive meningiomas. Our studies indicated that IMAT1 was highly expressed in invasive meningiomas compared with non-invasive meningiomas. In vitro studies showed that IMAT1 promoted meningioma cell invasion through the inactivation of the Krüppel-like factor 4 (KLF4)/hsa-miR22-3p/Snai1 pathway by acting as a sponge for hsa-miR22-3p, and IMAT1 knockdown effectively restored the tumor suppressive properties of KLF4 by preserving its tumor suppressor pathway. In vivo experiments confirmed that IMAT1 silencing could significantly inhibit the growth of subcutaneous tumors and prolong the survival period of tumor-bearing mice. Our findings demonstrated that the high expression of IMAT1 is the inherent reason for the loss of the tumor suppressive properties of KLF4 during meningioma progression. Therefore, we believe that IMAT1 may be a potential biological marker and treatment target for meningiomas.

Long non-coding RNA lncHUPC1 induced by FOXA1 promotes tumor progression by inhibiting apoptosis via miR-133b/SDCCAG3 in prostate cancer

Long non-coding RNAs (lncRNAs) were confirmed to be involved in regulating various malignant behaviors of tumor cells in prostate cancer (PCa). Using The Cancer Genome Atlas (TCGA) prostate adenocarcinoma datasets, several endogenous competing RNA (ceRNA) networks of lncRNA/miRNA/mRNA associated with the progression-free survival (PFS) and Gleason score (GS) were identified using bioinformatics analysis. lncRNA AC004447.2 (lncHUPC1, ENSG00000269131)/miR-133b/serologically defined colon cancer antigen-3 (SDCCAG3) was a newly identified ceRNA network that affected cell growth and apoptosis in PCa. Using q-PCR, lncHUPC1 and SDCCAG3 were found to be up-regulated in PCa cells, while miR-133b was down-regulated. The same results were found in tissue samples from 70 PCa cases. It was confirmed that the knockdown of lncHUPC1 increased the expression of miR-133b and decreased that of SDCCAG3, which further increased apoptosis and inhibited cell growth, while the miR-133b inhibitor partially reversed these effects. After transfection with miR-133b mimic after lncHUPC1-knockdown, the expression of miR-133b increased while that of SDCCAG3 reduced, and the apoptosis of the cells was more obvious and the growth of the cells was slower. Therefore, lncHUPC1 was confirmed to regulate SDCCAG3 by binding to miR-133b. Additionally, we found that the transcription factor Forkhead Box A1 (FOXA1) directly bound to the promoter of lncHUPC1 to activate it. In conclusion, the ceRNA network of lncHUPC1/miR-133b/SDCCAG3 affected the growth and apoptosis of PCa cells, and FOXA1 may be involved in the process as a transcription factor of lncHUPC1.

Long Non-Coding RNAs Might Regulate Phenotypic Switch of Vascular Smooth Muscle Cells Acting as ceRNA: Implications for In-Stent Restenosis

Coronary in-stent restenosis is a late complication of angioplasty. It is a multifactorial process that involves vascular smooth muscle cells (VSMCs), endothelial cells, and inflammatory and genetic factors. In this study, the transcriptomic landscape of VSMCs’ phenotypic switch process was assessed under stimuli resembling stent injury. Co-cultured contractile VSMCs and endothelial cells were exposed to a bare metal stent and platelet-derived growth factor (PDGF-BB) 20 ng/mL. Migratory capacity (wound healing assay), proliferative capacity, and cell cycle analysis of the VSMCs were performed. RNAseq analysis of contractile vs. proliferative VSMCs was performed. Gene differential expression (DE), identification of new long non-coding RNA candidates (lncRNAs), gene ontology (GO), and pathway enrichment (KEGG) were analyzed. A competing endogenous RNA network was constructed, and significant lncRNA–miRNA–mRNA axes were selected. VSMCs exposed to “stent injury” conditions showed morphologic changes, with proliferative and migratory capacities progressing from G0-G1 cell cycle phase to S and G2-M. RNAseq analysis showed DE of 1099, 509 and 64 differentially expressed mRNAs, lncRNAs, and miRNAs, respectively. GO analysis of DE genes showed significant enrichment in collagen and extracellular matrix organization, regulation of smooth muscle cell proliferation, and collagen biosynthetic process. The main upregulated nodes in the lncRNA-mediated ceRNA network were PVT1 and HIF1-AS2, with downregulation of ACTA2-AS1 and MIR663AHG. The PVT1 ceRNA axis appears to be an attractive target for in-stent restenosis diagnosis and treatment.

Network Integration Analysis and Immune Infiltration Analysis Reveal Potential Biomarkers for Primary Open-Angle Glaucoma

Primary open-angle glaucoma (POAG) is a progressive optic neuropathy and its damage to vision is irreversible. Therefore, early diagnosis assisted by biomarkers is essential. Although there were multiple researches on the identification of POAG biomarkers, few studies systematically revealed the transcriptome dysregulation mechanism of POAG from the perspective of pre- and post-transcription of genes. Here, we have collected multiple sets of POAG's aqueous humor (AH) tissue transcription profiles covering long non-coding RNA (lncRNA), mRNA and mircoRNA (miRNA). Through differential expression analysis, we identified thousands of significant differentially expressed genes (DEGs) between the AH tissue of POAG and non-glaucoma. Further, the DEGs were used to construct a competing endogenous RNA (ceRNA) regulatory network and 1,653 qualified lncRNA-miRNA-mRNA regulatory units were identified. Two ceRNA regulatory subnets were identified based on the random walk algorithm and revealed to be involved in the regulation of multiple complex diseases. At the pre-transcriptional regulation level, a transcriptional regulatory network was constructed and three transcription factors (FOS, ATF4, and RELB) were identified to regulate the expression of multiple genes and participate in the regulation of T cells. Moreover, we revealed the immune desert status of AH tissue for POAG patients based on immune infiltration analysis and identified a specific AL590666.2-hsa-miR-339-5p-UROD axis can be used as a biomarker of POAG. Taken together, the identification of regulatory mechanisms and biomarkers will contribute to the individualized diagnosis and treatment for POAG.

Characterizing the Copy Number Variation of Non-Coding RNAs Reveals Potential Therapeutic Targets and Prognostic Markers of LUSC

Lung squamous cell carcinoma (LUSC) has a poor clinical prognosis and a lack of available targeted therapies. Therefore, there is an urgent need to identify novel prognostic markers and therapeutic targets to assist in the diagnosis and treatment of LUSC. With the development of high-throughput sequencing technology, integrated analysis of multi-omics data will provide annotation of pathogenic non-coding variants and the role of non-coding sequence variants in cancers. Here, we integrated RNA-seq profiles and copy number variation (CNV) data to study the effects of non-coding variations on gene regulatory network. Furthermore, the 372 long non-coding RNAs (lncRNA) regulated by CNV were used as candidate genes, which could be used as biomarkers for clinical application. Nine lncRNAs including LINC00896, MCM8-AS1, LINC01251, LNX1-AS1, GPRC5D-AS1, CTD-2350J17.1, LINC01133, LINC01121, and AC073130.1 were recognized as prognostic markers for LUSC. By exploring the association of the prognosis-related lncRNAs (pr-lncRNAs) with immune cell infiltration, GPRC5D-AS1 and LINC01133 were highlighted as markers of the immunosuppressive microenvironment. Additionally, the cascade response of pr-lncRNA-CNV-mRNA-physiological functions was revealed. Taken together, the identification of prognostic markers and carcinogenic regulatory mechanisms will contribute to the individualized treatment for LUSC and promote the development of precision medicine.

Regulation of cancer progression by circRNA and functional proteins

Circular RNAs (circRNAs) are closed back-splicing products of precursor mRNA in eukaryotes. Compared with linear mRNAs, circRNAs have a special structure and stable expression. A large number of studies have provided different regulatory mechanisms of circRNAs in tumors. Challenges exist in understanding the control of circRNAs because of their sequence overlap with linear mRNA. Here, we survey the most recent progress regarding the regulation of circRNA biogenesis by RNA-binding proteins, one of the vital functional proteins. Furthermore, substantial circRNAs exert compelling biological roles by acting as protein sponges, by being translated themselves or regulating posttranslational modifications of proteins. This review will help further explore more types of functional proteins that interact with circRNA in cancer and reveal other unknown mechanisms of circRNA regulation.

Identification of Mutation Landscape and Immune Cell Component for Liver Hepatocellular Carcinoma Highlights Potential Therapeutic Targets and Prognostic Markers

Liver hepatocellular carcinoma (LIHC) is a primary malignancy, and there is a lack of effective treatment for advanced patients. Although numerous studies exist to reveal the carcinogenic mechanism of LIHC, few studies have integrated multi-omics data to systematically analyze pathogenesis and reveal potential therapeutic targets. Here, we integrated genomic variation data and RNA-seq profiles obtained by high-throughput sequencing to define high- and low-genomic instability samples. The mutational landscape was reported, and the advanced patients of LIHC were characterized by high-genomic instability. We found that the tumor microenvironment underwent metabolic reprograming driven by mutations accumulate to satisfy tumor proliferation and invasion. Further, the co-expression network identifies three mutant long non-coding RNAs as potential therapeutic targets, which can promote tumor progression by participating in specific carcinogenic mechanisms. Then, five potential prognostic markers (RP11-502I4.3, SPINK5, CHRM3, SLC5A12, and RP11-467L13.7) were identified by examining the association of genes and patient survival. By characterizing the immune landscape of LIHC, loss of immunogenicity was revealed as a key factor of immune checkpoint suppression. Macrophages were found to be significantly associated with patient risk scores, and high levels of macrophages accelerated patient mortality. In summary, the mutation-driven mechanism and immune landscape of LIHC revealed by this study will serve precision medicine.

Integrating Genetic and Transcriptomic Data to Reveal Pathogenesis and Prognostic Markers of Pancreatic Adenocarcinoma

Pancreatic adenocarcinoma (PAAD) is one of the deadliest malignancies and mortality for PAAD have remained increasing under the conditions of substantial improvements in mortality for other major cancers. Although multiple of studies exists on PAAD, few studies have dissected the oncogenic mechanisms of PAAD based on genomic variation. In this study, we integrated somatic mutation data and gene expression profiles obtained by high-throughput sequencing to characterize the pathogenesis of PAAD. The mutation profile containing 182 samples with 25,470 somatic mutations was obtained from The Cancer Genome Atlas (TCGA). The mutation landscape was generated and somatic mutations in PAAD were found to have preference for mutation location. The combination of mutation matrix and gene expression profiles identified 31 driver genes that were closely associated with tumor cell invasion and apoptosis. Co-expression networks were constructed based on 461 genes significantly associated with driver genes and the hub gene FAM133A in the network was identified to be associated with tumor metastasis. Further, the cascade relationship of somatic mutation-Long non-coding RNA (lncRNA)-microRNA (miRNA) was constructed to reveal a new mechanism for the involvement of mutations in post-transcriptional regulation. We have also identified prognostic markers that are significantly associated with overall survival (OS) of PAAD patients and constructed a risk score model to identify patients’ survival risk. In summary, our study revealed the pathogenic mechanisms and prognostic markers of PAAD providing theoretical support for the development of precision medicine.

Integrative network analysis of N6 methylation-related genes reveal potential therapeutic targets for spinal cord injury

The diagnosis of the severity of spinal cord injury (SCI) and the revelation of potential therapeutic targets are crucial for urgent clinical care and improved patient outcomes. Here, we analyzed the overall gene expression data in peripheral blood leukocytes during the acute injury phase collected from Gene Expression Omnibus (GEO) and identified six m6A regulators specifically expressed in SCI compared to normal samples. LncRNA-mRNA network analysis identified AKT2/3 and PIK3R1 related to m6A methylation as potential therapeutic targets for SCI and constructed a classifier to identify patients of SCI to assist clinical diagnosis. Moreover, FTO (eraser) and RBMX (reader) were found to be significantly down-regulated in SCI and the functional gene co-expressed with them was found to be involved in the signal transduction of multiple pathways related to nerve injury. Through the construction of the drug-target gene network, eight key genes were identified as drug targets and it was emphasized that fostamatinib can be used as a potential drug for the treatment of SCI. Taken together, our study characterized the pathogenesis and identified a potential therapeutic target of SCI providing theoretical support for the development of precision medicine.