Construction and analysis of dysregulated lncRNA-associated ceRNA network in colorectal cancer

Constructing lncRNA-miRNA-mRNA networks specific to individual cancer patients and finding prognostic biomarkers

BACKGROUND

The competitive endogenous RNA (ceRNA) hypothesis suggests that microRNAs (miRNAs) mediate a regulatory relation between long noncoding RNAs (lncRNAs) and messenger RNAs (mRNAs) which share similar miRNA response elements (MREs) to bind to the same miRNA. Since the ceRNA hypothesis was proposed, several studies have been conducted to construct a network of lncRNAs, miRNAs and mRNAs in cancer. However, most cancer-related ceRNA networks are intended for representing a general relation of RNAs in cancer rather than for a patient-specific relation. Due to the heterogeneous nature of cancer, lncRNA-miRNA-mRNA interactions can vary in different patients.

RESULTS

We have developed a new method for constructing a ceRNA network of lncRNAs, miRNAs and mRNAs, which is specific to an individual cancer patient and for finding prognostic biomarkers consisting of lncRNA-miRNA-mRNA triplets. We tested our method on extensive data sets of three types of cancer (breast cancer, liver cancer, and lung cancer) and obtained potential prognostic lncRNA-miRNA-mRNA triplets for each type of cancer.

CONCLUSIONS

Analysis of expression patterns of the RNAs involved in the triplets and survival rates of cancer patients revealed several interesting findings. First, even for the same cancer type, prognostic lncRNA-miRNA-mRNA triplets can be different depending on whether lncRNA and mRNA show opposite or similar expression patterns. Second, prognostic lncRNA-miRNA-mRNA triplets are often more predictive of survival rates than RNA pairs or individual RNAs. Our approach will be useful for constructing patient-specific lncRNA-miRNA-mRNA networks and for finding prognostic biomarkers from the networks.

Computational identification of potential inhibitors targeting cdk1 in colorectal cancer

Introduction: Despite improved treatment options, colorectal cancer (CRC) remains a huge public health concern with a significant impact on affected individuals. Cell cycle dysregulation and overexpression of certain regulators and checkpoint activators are important recurring events in the progression of cancer. Cyclin-dependent kinase 1 (CDK1), a key regulator of the cell cycle component central to the uncontrolled proliferation of malignant cells, has been reportedly implicated in CRC. This study aimed to identify CDK1 inhibitors with potential for clinical drug research in CRC. Methods: Ten thousand (10,000) naturally occurring compounds were evaluated for their inhibitory efficacies against CDK1 through molecular docking studies. The stability of the lead compounds in complex with CDK1 was evaluated using molecular dynamics simulation for one thousand (1,000) nanoseconds. The top-scoring candidates' ADME characteristics and drug-likeness were profiled using SwissADME. Results: Four hit compounds, namely, spiraeoside, robinetin, 6-hydroxyluteolin, and quercetagetin were identified from molecular docking analysis to possess the least binding scores. Molecular dynamics simulation revealed that robinetin and 6-hydroxyluteolin complexes were stable within the binding pocket of the CDK1 protein. Discussion: The findings from this study provide insight into novel candidates with specific inhibitory CDK1 activities that can be further investigated through animal testing, clinical trials, and drug development research for CRC treatment.

Constructing Integrative ceRNA Networks and Finding Prognostic Biomarkers in Renal Cell Carcinoma

Inspired by a newly discovered gene regulation mechanism known as competing endogenous RNA (ceRNA) interactions, several computational methods have been proposed to generate ceRNA networks. However, most of these methods have focused on deriving restricted types of ceRNA interactions such as lncRNA-miRNA-mRNA interactions. Competition for miRNA-binding occurs not only between lncRNAs and mRNAs but also between lncRNAs or between mRNAs. Furthermore, a large number of pseudogenes also act as ceRNAs, thereby regulate other genes. In this study, we developed a general method for constructing integrative networks of all possible interactions of ceRNAs in renal cell carcinoma (RCC). From the ceRNA networks we derived potential prognostic biomarkers, each of which is a triplet of two ceRNAs and miRNA (i.e., ceRNA-miRNA-ceRNA). Interestingly, some prognostic ceRNA triplets do not include mRNA at all, and consist of two non-coding RNAs and miRNA, which have been rarely known so far. Comparison of the prognostic ceRNA triplets to known prognostic genes in RCC showed that the triplets have a better predictive power of survival rates than the known prognostic genes. Our approach will help us construct integrative networks of ceRNAs of all types and find new potential prognostic biomarkers in cancer.

Long Non-Coding RNAs in Colorectal Cancer: Navigating the Intersections of Immunity, Intercellular Communication, and Therapeutic Potential

This comprehensive review elucidates the intricate roles of long non-coding RNAs (lncRNAs) within the colorectal cancer (CRC) microenvironment, intersecting the domains of immunity, intercellular communication, and therapeutic potential. lncRNAs, which are significantly involved in the pathogenesis of CRC, immune evasion, and the treatment response to CRC, have crucial implications in inflammation and serve as promising candidates for novel therapeutic strategies and biomarkers. This review scrutinizes the interaction of lncRNAs with the Consensus Molecular Subtypes (CMSs) of CRC, their complex interplay with the tumor stroma affecting immunity and inflammation, and their conveyance via extracellular vesicles, particularly exosomes. Furthermore, we delve into the intricate relationship between lncRNAs and other non-coding RNAs, including microRNAs and circular RNAs, in mediating cell-to-cell communication within the CRC microenvironment. Lastly, we propose potential strategies to manipulate lncRNAs to enhance anti-tumor immunity, thereby underlining the significance of lncRNAs in devising innovative therapeutic interventions in CRC.

A Blood-Based Immune Gene Signature with Prognostic Significance in Localized Prostate Cancer

Prostate cancer (PCa) is one of the most common male cancers worldwide and one of the deadliest if unsuccessfully treated. Τhe need for reliable, easily accessible immune-related molecular biomarkers that could be combined with clinically defined criteria, including PSA and Gleason score, to accurately predict PCa patients' clinical outcomes is emerging. Herein, we describe for the first time a blood-identified immune-related gene signature comprising eight upregulated multi-functional genes associated with poor prognosis. Next-generation sequencing (NGS) analysis of PCa patients' peripheral blood samples revealed a more than three-fold upregulation of each of the eight genes as compared to samples originating from healthy donors. The construction of gene and protein interaction networks revealed different extents of the functional implications of these genes in the regulation of cell proliferation and immune responses. Analysis of the available data from The Cancer Genome Atlas (TCGA) regarding gene expression and survival of prostate adenocarcinoma (PRAD) and pan-cancer (PANCAN) patients revealed that intra-tumoral upregulation of this eight-gene signature (8-GS) was associated with poor 5-year progression-free intervals in PCa patients, even in those with high Gleason scores, and also with an unfavorable prognosis for cancer patients irrespective of the cancer type and even in the early stages. These observations suggest that further investigation of the 8-GS prospectively in randomized clinical trials, in which clinical benefit in terms of evaluating time to disease progression can be assessed, is warranted.

Identification of aberrantly expressed lncRNAs and ceRNA networks in multiple myeloma: a combined high-throughput sequencing and microarray analysis

Background

This study aimed to explore the potential effects of long non-coding RNAs (lncRNAs) in multiple myeloma (MM) patients using two detection methods: high-throughput sequencing and microarray.

Methods

In this study, lncRNAs were detected in 20 newly diagnosed MM patients, with 10 patients analyzed by whole transcriptome-specific RNA sequencing and 10 patients analyzed by microarray (Affymetrix Human Clariom D). The expression levels of lncRNAs, microRNAs, and messenger RNAs (mRNAs) were analyzed, and the differentially expressed lncRNAs identified by both methods were selected. The significant differentially expressed lncRNAs were further validated using PCR.

Results

This study established the aberrant expression of certain lncRNAs involved in the occurrence of MM, with AC007278.2 and FAM157C showing the most significant differences. The top 5 common pathways identified by the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were the chemokine signaling pathway, inflammatory mediator regulation, Th17 cell differentiation, apoptosis, and NF-kappa B signaling pathway. Furthermore, three microRNAs (miRNAs) (miR-4772-3p, miR-617, and miR-618) were found to constitute competing endogenous RNA (ceRNA) networks in both sequencing and microarray analyses.

Conclusions

By the combination analysis, our understanding of lncRNAs in MM will be increased significantly. More overlapping differentially expressed lncRNAs were found to predict therapeutic targets precisely.

Finding miRNA-RNA Network Biomarkers for Predicting Metastasis and Prognosis in Cancer

Despite remarkable progress in cancer research and treatment over the past decades, cancer ranks as a leading cause of death worldwide. In particular, metastasis is the major cause of cancer deaths. After an extensive analysis of miRNAs and RNAs in tumor tissue samples, we derived miRNA-RNA pairs with substantially different correlations from those in normal tissue samples. Using the differential miRNA-RNA correlations, we constructed models for predicting metastasis. A comparison of our model to other models with the same data sets of solid cancer showed that our model is much better than the others in both lymph node metastasis and distant metastasis. The miRNA-RNA correlations were also used in finding prognostic network biomarkers in cancer patients. The results of our study showed that miRNA-RNA correlations and networks consisting of miRNA-RNA pairs were more powerful in predicting prognosis as well as metastasis. Our method and the biomarkers obtained using the method will be useful for predicting metastasis and prognosis, which in turn will help select treatment options for cancer patients and targets of anti-cancer drug discovery.

A Prognostic Cuproptosis-Related LncRNA Signature for Colon Adenocarcinoma

Background

Cuproptosis, a recently discovered form of cell death, is caused by copper levels exceeding homeostasis thresholds. Although Cu has a potential role in colon adenocarcinoma (COAD), its role in the development of COAD remains unclear.

Methods

In this study, 426 patients with COAD were extracted from the Cancer Genome Atlas (TCGA) database. The Pearson correlation algorithm was used to identify cuproptosis-related lncRNAs. Using the univariate Cox regression analysis, the least absolute shrinkage and selection operator (LASSO) was used to select cuproptosis-related lncRNAs associated with COAD overall survival (OS). A risk model was established based on the multivariate Cox regression analysis. A nomogram model was used to evaluate the prognostic signature based on the risk model. Finally, mutational burden and sensitivity analyses of chemotherapy drugs were performed for COAD patients in the low- and high-risk groups.

Result

Ten cuproptosis-related lncRNAs were identified and a novel risk model was constructed. A signature based on ten cuproptosis-related lncRNAs was an independent prognostic predictor for COAD. Mutational burden analysis suggested that patients with high-risk scores had higher mutation frequency and shorter survival.

Conclusion

Constructing a risk model based on the ten cuproptosis-related lncRNAs could accurately predict the prognosis of COAD patients, providing a fresh perspective for future research on COAD.

The expression of long non-coding RNA LINC01389, LINC00365, RP11-138J23.1, and RP11-354K4.2 in gastric cancer and their impacts on EMT

BACKGROUND

Epithelial cancers acquire the epithelial to mesenchymal transition (EMT), which leads tumor cells to invade and metastasize to adjacent and distant tissues. The mechanisms involved in EMT phenotype are controlled by numerous markers as well as signalling pathways. Recently, long non-coding RNAs (lncRNAs) were introduced that play the regulatory role in EMT via crosstalk with EMT-related transcription factors and signalling pathways. The present study aimed to investigate the expression of four lncRNAs in human GC and elucidate their probable role in EMT procedure and the pathogenesis of gastric cancer (GC).

METHODS

The expression profile of lncRNAs (LINC01389, LINC00365, RP11-138J23.1, and RP11-354K4.2) and mRNAs (TWIST1, MMP13, MAML1, CD44s, and SALL4) between eighty-three GC and adjacent non-cancerous tissues were assessed by quantitative real-time PCR.

RESULTS

The significant downregulation of LINC00365 (66.3%) and RP11-354K4.2 (62.7%) were observed in GC samples; while the upregulation of LINC01389, RP11-138J23.1, TWIST1, MMP13, MAML1, CD44s, and SALL4 were found in 67.5%, 45.8%, 56.6%, 44.6%, 59%, 55.4%, and 62.7% tumors samples at the mRNA level, respectively. Dysregulation of these lncRNAs and EMT-related markers was significantly related to each other in a variety of clinicopathological features of patients (P < 0.05), indicating positive correlations between LINC01389, LINC00365, RP11-138J23.1, and RP11-354K4.2 with EMT status in GC.

CONCLUSION

These EMT-regulating lncRNAs may play a key role in transforming gastric epithelial to mesenchymal phenotype and can be novel therapeutic targets for GC. Our results highlight the importance of discovering new lncRNAs involved in gastric carcinogenesis. Detailed molecular mechanisms of these noncoding-coding markers in GC are urgently required.

Systematic Pan-Cancer Analysis Identifies CDK1 as an Immunological and Prognostic Biomarker

Cyclin-dependent kinase 1 (CDK1) plays an important role in cancer development, progression, and the overall process of tumorigenesis. However, no pan-cancer analysis has been reported for CDK1, and the predictive role of CDK1 in immune checkpoint inhibitors (ICIs) therapy response remains unexplored. Thus, in this study, we first investigated the potential oncogenic role of CDK1 in 33 tumors by multidimensional bioinformatics analysis based on The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets. Bioinformatic analysis and immunohistochemical experiments confirmed that CDK1 is significantly upregulated in most common cancers and is strongly associated with prognosis. Further analysis indicated that CDK1 may influence tumor immunity mainly by mediating the degree of tumor infiltration of immune-associated cells, and the effect of CDK1 on immunity is diverse across tumor types in tumor microenvironment. CDK1 was also positively correlated with tumor mutational burden (TMB) and microsatellite instability (MSI) in certain cancer types, linking its expression to the assessment of possible treatment response. The results of the pan-cancer analysis study showed that the CDK1 gene was positively associated with the expression of three classes of RNA methylation regulatory proteins, and affects RNA function through multiple mechanisms of action and plays an important role in the posttranscriptional regulation of the tumor microenvironment. These findings shed light on the role of the CDK1 gene in cancer progression and provide information to further study the CDK1 gene as a potential target for pan-cancer.

Comprehensive analysis of ceRNA network of ERCC4 in colorectal cancer

Objective

ERCC4 is one of the most significant molecules of Nucleotide Excision Repair (NER), which has been researched due to its high expression in colorectal cancer (CRC). This study aimed to find out the ceRNA (competitive endogenous RNA) network of ERCC4 in CRC.

Methods and Materials

Pan cancer mRNA expression of ERCC4 was evaluated using TCGA database. The protein expression of ERCC4 was evaluated based on the Human Protein Atlas (HPA). We screened DElncRNAs and DEmiRNAs in two groups of ERCC4^high and ERCC4^low expression in CRC. Then a lncRNA-miRNA-ERCC4 regulatory network was constructed based on DElncRNAs and DEmiRNAs using Starbase database and visualized by Cytoscape software. Kaplan-Meier analysis was performed to evaluate the prognostic value of the ceRNA network. Further, RT-PCR was performed to validate the expression of the representative molecules in the ceRNA network in CRC and normal tissues. The relationship between drug sensitivity and these molecules were also evaluated using RNAactDrug database.

Results

ERCC4 was overexpressed in a variety of tumors at mRNA levels, including CRC. High expression of ERCC4 was also observed on protein level in CRC. A total of 1,885 DElncRNAs and 68 DEmiRNAs were identified from CRC samples in ERCC4^high and ERCC4^low expression groups. Predicted by the Starbase database, we got interacting miRNAs and lncRNAs of ERCC4 from the DEmiRNAs and DElncRNAs, and a lncRNA-miRNA-ERCC4 regulatory network was constructed. Kaplan-Meier survival curves results showed that miR-200c-3p (hazard ratio [HR] = 0.62, P = 0.032), MALAT1 (HR = 1.54, P = 0.016), and AC005520.2 (hazard ratio [HR] = 1.75, P = 0.002) were significantly associated with the prognosis of CRC. After validation by RT-PCR, we found that ERCC4 and MALAT1 were up-regulated in CRC compared with normal tissues, while miR-200c-3p was down-regulated. A strong negative correlation was observed between MALAT1 and miR-200c-3p. Drug sensitivity analysis showed that ERCC4, miR-200c and MALAT1 were all associated with Cisplatin.

Conclusion

We constructed a ceRNA network of ERCC4 in CRC, of which the MALAT1-miR-200c-3p-ERCC4 axis may be involved in the development, prognosis and chemotherapy sensitivity of CRC. These findings might provide novel clues and insights on the molecular mechanisms of ERCC4 and NER pathway in CRC.

LncRNA HCG18 upregulates TRAF4/TRAF5 to facilitate proliferation, migration and EMT of epithelial ovarian cancer by targeting miR-29a/b

Background

Although long noncoding RNA HLA complex group 18 (lncRNA HCG18) has been suggested to regulate cell growth in several tumours, the function of HCG18 in epithelial ovarian cancer (EOC) and its mechanism are still unclear.

Methods

shRNAs were applied to reduce HCG18 and related genes. For overexpression of miRNA, a miRNA mimic was transfected into cells. Quantitative real-time PCR (qRT–PCR) was used to detect levels of HCG18, miR-29a/b, and mRNAs. MTT, colony formation, wound healing and Transwell assays were used to evaluate cell proliferation, migration and invasion, respectively. A luciferase reporter assay was utilized to evaluate NF-κB activity and the binding of miRNAs with HCG18 or TRAF4/5. BALB nude mice injected with cells stably expressing shHCG18 or shNC were used for in vivo modelling. Subcutaneous tumour growth was monitored in nude mice, and immunohistochemistry (IHC) was used to determine expression of the proliferation marker Ki67.

Results

Abnormal expression of HCG18 and miR-29a/b was observed in EOC tissues. Knockdown of HCG18 using shRNA inhibited proliferation, migration, EMT and the proinflammatory pathway in EOC cells. miR-29a/b mimics and TRAF4/5 knockdown exhibited effects similar to HCG18 knockdown. Further experiments suggested that HCG18 directly targets miR-29a/b and upregulates TRAF4/5 expression, which are inhibited by targeting miR-29a/b. Moreover, overexpression of TRAF4/5 antagonized the inhibitory effect of HCG18 knockdown, suggesting that they are involved in HCG18-mediated oncogenic effects. Silencing HCG18 reduced tumour size and levels of Ki67 and TRAF4/5 while increasing miR-29a/b levels in vivo.

Conclusions

Taken together, our data revealed an oncogenic signalling pathway mediated by HCG18 in ovarian cell lines, which functions as a ceRNA of miR-29a/b and thus derepresses expression levels of TRAF4/5, facilitating NF-κB pathway-mediated promotion of EOC cell proliferation and migration.

Supplementary Information

The online version contains supplementary material available at 10.1186/s10020-021-00415-y.

Identification of Novel Biomarkers for Predicting Prognosis and Immunotherapy Response in Head and Neck Squamous Cell Carcinoma Based on ceRNA Network and Immune Infiltration Analysis

Objectives

Patients with head and neck squamous cell carcinoma (HNSCC) have poor prognosis and show poor responses to immune checkpoint (IC) inhibitor (ICI) therapy. Competing endogenous RNA (ceRNA) networks, tumor-infiltrating immune cells (TIICs), and ICIs may influence tumor prognosis and response rates to ICI therapy. This study is aimed at identifying prognostic and IC-related biomarkers and key TIIC signatures to improve prognosis and ICI therapy response in HNSCC patients.

Methods and Results

Ninety-five long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and 1746 mRNAs were identified using three independent methods. We constructed a ceRNA network and estimated the proportions of 22 immune cell subtypes. Ten ceRNAs were related to prognosis according to Kaplan–Meier analysis. Two risk signatures based, respectively, on nine ceRNAs (ANLN, CFL2, ITGA5, KDELC1, KIF23, NFIA, PTX3, RELT, and TMC7) and three immune cell types (naïve B cells, neutrophils, and regulatory T cells) via univariate Cox regression, least absolute shrinkage and selection operator, and multivariate Cox regression analyses could accurately and independently predict the prognosis of HNSCC patients. Key mRNAs in the ceRNA network were significantly correlated with naïve B cells and regulatory T cells and with stage, grade, and immune and molecular subtype. Eight IC genes exhibited higher expression in tumor tissues and were correlated with eight key mRNAs in the ceRNA network in HNSCC patients with different HPV statuses according to coexpression and TIMER 2.0 analyses. Most drugs were effective in association with expression of these key signatures (ANLN, CFL2, ITGA5, KIF23, NFIA, PTX3, RELT, and TMC7) based on GSCALite analysis. The prognostic value of key biomarkers and associations between key ceRNAs and IC genes were validated using online databases. Eight key ceRNAs were confirmed to predict response to ICI in other cancers based on TIDE analysis.

Conclusions

We constructed two risk signatures to accurately predict prognosis in HNSCC. Key IC-related signatures may be associated with response to ICI therapy. Combinations of ICIs with inhibitors of eight key mRNAs may improve survival outcomes of HNSCC patients.

Construction of a long noncoding RNA-based competing endogenous RNA network and prognostic signatures of left- and right-side colon cancer

Background

Cancers located on the right and left sides of the colon have distinct clinical and molecular characteristics. This study aimed to explore the regulatory mechanisms of location-specific long noncoding RNAs (lncRNAs) as competing endogenous RNAs (ceRNAs) in colon cancer and identify potential prognostic biomarkers.

Method

Differentially expressed lncRNAs (DELs), miRNAs (DEMs), and genes (DEGs) between right- and left-side colon cancers were identified by comparing RNA sequencing profiles. Functional enrichment analysis was performed for the DEGs, and a ceRNA network was constructed. Associations between DELs and patient survival were examined, and a DEL-based signature was constructed to examine the prognostic value of these differences. Clinical colon cancer tissues and Gene Expression Omnibus (GEO) datasets were used to validate the results.

Results

We identified 376 DELs, 35 DEMs, and 805 DEGs between right- and left-side colon cancers. The functional enrichment analysis revealed the functions and pathway involvement of DEGs. A ceRNA network was constructed based on 95 DEL–DEM–DEG interactions. Three DELs (LINC01555, AC015712, and FZD10-AS1) were associated with the overall survival of patients with colon cancer, and a prognostic signature was established based on these three DELs. High risk scores for this signature indicated poor survival, suggesting that the signature has prognostic value for colon cancer. Examination of clinical colon cancer tissues and GEO dataset analysis confirmed the results.

Conclusion

The ceRNA regulatory network suggests roles for location-specific lncRNAs in colon cancer and allowed the development of an lncRNA-based prognostic signature, which could be used to assess prognosis and determine treatment strategies in patients with colon cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-01901-3.

Inhibition of CDK1 Reverses the Resistance of 5-Fu in Colorectal Cancer

Introduction

Although the survival rate of colorectal cancer (CRC) patients can be improved by surgery, radiotherapy, and chemotherapy, the resistance to 5-fluorouracil (5-Fu) affects the effect of chemotherapy and the prognosis of patients. An increasing number of studies showed that 5-Fu resistance was the main reason for the failure of colorectal cancer treatment. The poor prognosis of colorectal cancer greatly harms people’s health. This study aimed to clarify the correlation between cyclin-dependent kinase 1 (CDK1) and 5-Fu-induced tumor resistance.

Materials and Methods

Cell proliferation and invasion experiments showed that down-regulation of CDK1 inhibited fluorouracil-resistant CRC cell proliferation. The expression level of CDK1 was detected in 5-Fu-resistant CRC cells in vitro. Tumor growth was inhibited by down-regulation of CDK1 in tumor xenograft mouse models.

Results

We found that CDK1 was highly expressed in tumor tissues, especially in fluorouracil-resistant tissues. We also confirmed that the differential expression of 5-Fu in tumor tissues was related to tumor site, lymph node metastasis and stage. CDK1 promoted migration, invasion and inhibited apoptosis in 5-Fu-resistant CRC cells. Down-regulation of CDK1 inhibited fluorouracil-resistant CRC cell proliferation and tumorigenesis in vivo.

Conclusion

High expression of CDK1 may lead to poor clinical prognosis, and inhibition of CDK1 enhances 5-Fu sensitivity in CRC. Our research suggested that CDK1 may be used to predict 5-Fu efficacy and as a therapeutic target for CRC.

Genome-Wide Analysis of the FOXA1 Transcriptional Network Identifies Novel Protein-Coding and Long Noncoding RNA Targets in Colorectal Cancer Cells

Differentiation status of tumors is correlated with metastatic potential and malignancy. FOXA1 (forkhead box A1) is a transcription factor known to regulate differentiation in certain tissues. Here, we investigate FOXA1 function in human colorectal cancer (CRC). We found that FOXA1 is robustly expressed in the normal human colon but significantly downregulated in colon adenocarcinoma. Applying FOXA1 chromatin immunoprecipitation coupled with deep sequencing and transcriptome analysis upon FOXA1 knockdown in well-differentiated CRC cells and FOXA1 overexpression in poorly differentiated CRC cells, we identified novel protein-coding and lncRNA genes regulated by FOXA1. Among the numerous novel FOXA1 targets we identified, we focused on CEACAM5, a tumor marker and facilitator of cell adhesion. We show that FOXA1 binds to a distal enhancer downstream of CEACAM5 and strongly activates its expression. Consistent with these data, we show that FOXA1 inhibits anoikis in CRC cells. Collectively, our results uncover novel protein-coding and noncoding targets of FOXA1 and suggest a vital role of FOXA1 in enhancing CEACAM5 expression and anoikis resistance in CRC cells.

Integration analysis of long non-coding RNA (lncRNA) role in tumorigenesis of colon adenocarcinoma

Background

Colon adenocarcinoma (COAD) is one of the most common gastrointestinal cancers globally. Molecular aberrations of tumor suppressors and/or oncogenes are the main contributors to tumorigenesis. However, the exact underlying mechanisms of COAD pathogenesis are clearly not known yet. In this regard, there is an urgent need to indicate promising potential diagnostic and prognostic biomarkers in COAD patients.

Methods

In the current study, level 3 RNA-Seq and miR-Seq data and corresponding clinical data of colon adenocarcinoma (COAD) were retrieved from the TCGA database. The “limma” package in R software was utilized to indicate the differentially expressed genes. For in silico functional analysis, GO and KEGG signaling pathways were conducted. PPI network was constructed based on the STRING online database by Cytoscape 3.7.2. A ceRNA network was also constructed by “GDCRNATools” package in R software. Kaplan-Meier survival analysis (log-rank test) and ROC curve analysis were used to indicate the diagnostic and prognostic values of the biomarkers.

Results

The differential expression data demonstrated that 2995 mRNAs, 205 lncRNAs, and 345 miRNAs were differentially expressed in COAD. The GO and KEGG pathway analysis indicated that the differentially expressed mRNAs were primarily enriched in canonical processes in cancer. The PPI network showed that the CDKN2A, CCND1, MYC, E2F, CDK4, BRCA2, CDC25B, and CDKN1A proteins were the critical hubs. In addition, the Kaplan-Meier analysis revealed that 215 mRNAs, 14 lncRNAs, and 39 miRNAs were associated with overall survival time in the patients. Also, the ceRNA network data demonstrated that three lncRNAs including MIR17HG, H19, SNHG1, KCNQ1OT1, MALAT1, GAS5, SNHG20, OR2A1-AS1, and MAGI2-AS3 genes were involved in the development of COAD.

Conclusions

Our data suggested several promising lncRNAs in the diagnosis and prognosis of patients with COAD.

A Novel Circular RNA, circ_0005394, Predicts Unfavorable Prognosis and Contributes to Hepatocellular Carcinoma Progression by Regulating miR-507/E2F3 and miR-515-5p/CXCL6 Signaling Pathways

Objective

Circular RNAs (circRNAs) play a key role in cancer development and progression. Previously, circ_0005394 was found to be highly expressed in hepatocellular carcinoma (HCC) screened by circRNA microarray. However, the research with regard to the functions and mechanisms of circ_0005394 in HCC remains unknown.

Materials and Methods

The expression of circ_0005394 in HCC was measured by qRT-PCR. The clinical relevance was evaluated by Fisher's exact test, Kaplan-Meier curves, and Cox regression model. Gain/loss-of function assays were performed to elucidate the functions of circ_0005394 in Huh-7 and HepG2 cells. Dual-luciferase reporter assay was applied to reveal the mechanism of circ_0005394.

Results

circ_0005394 expression was higher in HCC tissues and cells than noncancerous samples and normal cell line, respectively. High expression of circ_0005394 was associated with larger tumor size, more advanced TNM stages, and poorer overall survival for the patients with HCC. Gain/loss-of function assays demonstrated its oncogenic role in cell growth, apoptosis, migration, and invasion. Mechanistically, miR-507 and miR-515-5p could be sponged by circ_0005394. Furthermore, E2F Transcription Factor 3 (E2F3) and C-X-C motif chemokine ligand 6 (CXCL6) were confirmed as the target of miR-507 and miR-515-5p, respectively. Rescue assay indicated that circ_0005394 facilitated HCC growth and invasion by regulating miR-507/E2F3 and miR-515-5p/CXCL6 signaling pathways.

Conclusion

This study uncovered an important role of circ_0005394 in regulating HCC progression, providing a novel perspective for clarifying its pathogenesis.

Identification and Validation of an Immune-Related RNA Signature to Predict Survival of Patients With Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a heterogeneous disease characterized by different molecular subgroups and clinical features. Therefore, it is important to uncover reliable molecular biomarkers for distinguishing different risk patient subgroup. Here, we conducted a multi-omics analysis to examine the joint predictive power of a multi-type RNA signature in the prognosis of HNSCC patients through integration analysis of mRNA, miRNA, and lncRNA expression profiles and clinical data in a large number of HNSCC patients. A multi-type RNA signature (15SigRS) was constructed which can classify patients into the high-risk group and low-risk group with the significantly different outcome [hazard ratio (HR) = 2.718, 95% confidence interval (CI), 2.258–3.272, p < 0.001] in the discovery data set, and subsequently validated in the Cancer Genome Atlas (TCGA) testing data set (HR = 1.299, 95% CI, 1.170–1.442, p < 0.001) and another independent GSE65858 data set (HR = 1.077, 95% CI, 1.016–1.143, p = 0.013). Further multivariate Cox regression analysis and stratification analysis demonstrated the independence of predictive performance of the 15SigRS relative to conventional clinicopathological factors. Furthermore, the 15SigRS has a prior performance in prognostic prediction than other single RNA type-based signatures. Functional analysis suggested that the 15SigRS are involved in immune- or metabolism-related KEGG pathways. In summary, our study demonstrated the potential application of mixed RNA types as molecular markers for predicting the outcome of cancer patients.

LINC00365 promotes colorectal cancer cell progression through the Wnt/β-catenin signaling pathway

In the past decade, substantial evidence established that long noncoding RNAs are serious about mediating the evolution of malignancies. In previous studies, LINC00365, which has not been reported in colorectal cancer (CRC), was selected using the bioinformatics analysis in GSE109454 and GSE41655 data sets. However, the function and mechanism of LINC00365 are still obscure. In our study, LINC00365 was found upregulated in CRC specimens and intimately connected with the prognosis of patients with CRC. In addition, LINC00365 overexpression enhances the cell abilities of proliferation, migration, and invasion in vitro. Meanwhile, mechanistic studies showed that LINC00365 might involve in CRC cell progression by mediating the Wnt/β-catenin pathway. Furthermore, LINC00365 upregulation increased CDK1 protein expression. In conclusion, this study suggests that LINC00365 acts as a vital part in facilitating CRC progression and might play as a therapeutic target for patients with CRC.

miR-548c-5p inhibits colorectal cancer cell proliferation by targeting PGK1

Accumulating studies have implicated that microRNAs (miRNAs) are involved in the pathogenesis of colorectal cancer (CRC). However, the role of miR-548c-5p, a novel identified miRNA in malignancies, in colorectal carcinogenesis remains largely unknown. The present study is aimed to investigate the effect and molecular mechanism of miR-548c-5p in CRC by a sequence of cellular experiments. miR-548c-5p was significantly downregulated, whereas phosphoglycerate kinase 1 (PGK1), a key enzyme for glycolysis, was obviously upregulated in peripheral blood mononuclear cells and cancer tissues from patients with CRC. Besides, miR-548c-5p and PGK1 were negatively associated with each other. The luciferase reporter assay revealed that PGK1 was a targeted gene of miR-548c-5p. Moreover, the proliferation and generation of inflammatory cytokines (TNF-α and IL-6) were significantly inhibited in miR-548c-5p-overexpressed SW480 CRC cells stimulated by lipopolysaccharide (LPS). Accordingly, miR-548c-5p may serve as a cancer suppressor in CRC by targeting PGK1.