Robust analysis of novel mRNA-lncRNA cross talk based on ceRNA hypothesis uncovers carcinogenic mechanism and promotes diagnostic accuracy in esophageal cancer

Identification and diagnostic potential of hsa_circ_101303 in colorectal cancer: unraveling a regulatory network

BACKGROUND

The role of novel circular RNAs (circRNAs) in colorectal cancer (CRC) remains to be determined. This study aimed to identify a novel circRNA involved in CRC pathogenesis, assess its diagnostic value, and construct a regulatory network.

METHODS

Differential expression analysis was conducted using circRNA datasets to screen for differentially expressed circRNAs. The expression of selected circRNAs was validated in external datasets and clinical samples. Diagnostic value of plasma circRNA levels in CRC was assessed. A competing endogenous RNA (ceRNA) network was constructed for the circRNA using TCGA dataset.

RESULTS

Analysis of datasets revealed that hsa_circ_101303 was significantly overexpressed in CRC tissues compared to normal tissues. The upregulation of hsa_circ_101303 in CRC tissues was further confirmed through the GSE138589 dataset and clinical samples. High expression of hsa_circ_101303 was associated with advanced N stage, M stage, and tumor stage in CRC. Plasma levels of hsa_circ_101303 were markedly elevated in CRC patients and exhibited moderate diagnostic ability for CRC (AUC = 0.738). The host gene of hsa_circ_101303 was also found to be related to the TNM stage of CRC. Nine miRNAs were identified as target miRNAs for hsa_circ_101303, and 27 genes were identified as targets of these miRNAs. Subsequently, a ceRNA network for hsa_circ_101303 was constructed to illustrate the interactions between the nine miRNAs and 27 genes.

CONCLUSIONS

The study identifies hsa_circ_101303 as a highly expressed circRNA in CRC, which is associated with the progression of the disease. Plasma levels of hsa_circ_101303 show promising diagnostic potential for CRC. The ceRNA network for hsa_circ_101303 provides valuable insights into the regulatory mechanisms underlying CRC.

A review on the role of ADAMTS9-AS2 in different disorders

Recent decade has seen a tremendous progress in identification of the role of different long non-coding RNAs (lncRNAs) in human pathologies. ADAMTS9-AS2 is an example of lncRNAs with different roles in human disorders. It is mostly acknowledged as a tumor suppressor lncRNA in different types of cancers. However, it has been reported to be up-regulated in tongue squamous cell carcinoma, salivary adenoid cystic carcinoma and glioblastoma. Moreover, ADAMTS9-AS2 is possibly involved in the pathoetiology of pulpitis, acute ischemic stroke, type 2 diabetes and its complications. This lncRNA sponges miR-196b-5p, miR-223-3p, miR-130a-5p, miR-600, miR-223-3p, miR-27a-3p, miR-32, miR-143-3p, miR-143-3p and miR-182-5p in order to regulate downstream mRNAs. This review aims at summarization of the role of ADAMTS9-AS2 in different disorders with a particular focus on its diagnostic and prognostic values.

Inhibitory role of LINC00332 in gastric cancer progression through regulating cell EMT and stemness

OBJECTIVE

Gastric cancer (GC) is one of the most common lethal malignancies worldwide. The molecular mechanisms underlying GC early detection are poorly understood. Identifying potential coding and non-coding markers and related pathways in the GC progression is essential. Some Long non-coding RNAs (lncRNAs) reportedly play vital roles during gastric GC development. However, the clinical significance and biological function of LINC00332 in GC remain largely unclear.

METHODS

The gene expression patterns of GC from an RNAseq dataset (GSE122401) were retrieved from the Gene Expression Omnibus (GEO) database to recognize differentially expressed genes (DEGs) and lncRNAs (DELs) between normal and GC samples through several bioinformatic analysis. The expression of LINC00332 and MMP-13 as a target gene was quantified in fresh frozen tissues obtained from GC patients. In addition, we investigated the potential function of LINC00332 in silico and in vitro.

RESULTS

The expressions of LINC00332 and MMP-13 were significantly downregulated and upregulated in GC tissues, respectively. A significant inverse correlation between LINC00332 and MMP-13 mRNA expression was observed in tumor samples. The mRNA expression level of mesenchymal markers, stem cell factors, and MMP genes were significantly decreased after the LINC00332 ectopic expression, while epithelial markers expression was significantly increased. The LINC00332 overexpression markedly repressed proliferation, migration, and invasion and did not induce apoptosis in AGS cells. In addition, LINC00332 overexpression notably promoted the E-cadherin protein expression. Moreover, LINC00332 significantly decreased the cisplatin resistance.

CONCLUSION

Our findings indicated that LINC00332 may be a critical anti-EMT factor and provided a new efficient therapeutic strategy for GC treatment.

Construction of a Prognostic Model for Hypoxia-Related LncRNAs and Prediction of the Immune Landscape in the Digestive System Pan-Cancer

Digestive system pan-cancer is a general term for digestive system tumors including colorectal carcinoma (CRC), esophageal carcinoma (ESCA), stomach adenocarcinoma (STAD), and liver hepatocellular carcinoma (LIHC). Since the anatomical location, function and metabolism are closely related, there may be similarities in development and progression of these tumors. Hypoxia is the consequence of an imbalance between oxygen demand and supply, and intracellular hypoxia is associated with malignant progression, treatment resistance, and poor prognosis in tumors. Therefore, an urgent and challenging task is to investigate the molecular mechanisms associated with hypoxia in digestive system pan-cancer for the prognosis and treatment of digestive tract tumors. In this study, we identified 18 hypoxia-related lncRNAs (HRlncRNAs) by co-expression analysis between hypoxia genes and lncRNAs from digestive system pan-cancer. Six HRlncRNAs were then obtained using lasso regression and multivariate cox analysis to construct a prognostic model. Next, the Akaike information criterion (AIC) values for 3-year receiver operating curve (ROC) were counted to determine the cut-off point and establish an optimal model to distinguish between high- or low-risk groups among patients with digestive system pan-cancer. To evaluate the stability of the prognosis model, we validated it in terms of survival outcomes, clinicopathological stage, tumor-infiltrating immune cells, immune checkpoint inhibitors (ICIs) and anticancer drugs sensitivity. The results suggested that high- risk group had a worse prognosis and a more positive association with tumor-infiltrating immune cells such as B cells, cancer-associated fibroblasts, endothelial cells, monocytes, macrophages and bone marrow dendritic cells in digestive system pan-cancer. Immune checkpoint inhibitors (ICIs) related biomarkers discovered that high-risk group was positively correlated with high expression of HAVCR2 in digestive system pan-cancer. The anticancer drugs sensitivity analysis showed that the high-risk group was associated with the lower half-inhibitory centration (IC50) of Imatinib in digestive system pan-cancer. In conclusion, the prognostic model of HRlncRNAs showed a promising clinical prediction value and may provide a useful reference for the diagnosis and treatment of the digestive system tumors.

ADAMTS9-AS2 regulates PPP1R12B by adsorbing miR-196b-5p and affects cell cycle-related signaling pathways inhibiting the malignant process of esophageal cancer

This study explored the mechanism that ADAMTS9-AS2/miR-196b-5p/PPP1R12B/cell cycle pathway axis in inhibiting the malignant progression of esophageal cancer (EC), providing a new idea for targeted molecular therapy of EC. The expression data of EC tissue were acquired from TCGA database. The target lncRNA, downstream miRNA and its target gene were determined by bioinformatics analysis. ADAMTS9-AS2, miR-196b-5p and PPP1R12B levels in EC tissue and cells were assayed through qRT-PCR. Western blot was applied to assess protein level of PPP1R12B in cells and tissues, as well as protein expression of CDK1, cyclin A2, cyclin B1 and Plk1 in EC cells. Cell proliferation was assayed via CCK-8 assay. Cell cycle distribution was analyzed by flow cytometry. Cell migratory and invasive abilities were measured through scratch healing and transwell assays. Pearson correlation analysis was utilized to analyze relationship among ADAMTS9-AS2, miR-196b-5p and PPP1R12B. RIP was introduced to assess binding among the three. Dual-luciferase assay was utilized to verify targeted binding sites. The tumor formation in nude mice assay was utilized to detect tumorigenesis of EC cells in vivo. ADAMTS9-AS2 was significantly lowly expressed while miR-196b-5p was increased in EC tissue and cells. ADAMTS9-AS2 bound to miR-196b-5p and constrained its expression. Overexpressed ADAMTS9-AS2 inhibited EC cell malignant progression via downregulating miR-196b-5p, while overexpressed miR-196b-5p reversed this inhibitory effect. ADAMTS9-AS2 modulated PPP1R12B level by competitively inhibiting miR-196b-5p. PPP1R12B played a modulatory role in EC by inhibiting cell cycle pathway. Overexpressed ADAMTS9-AS2 regulated the tumor-forming ability of EC cells in vivo through miR-196b-5p/PPP1R12B/cell cycle signaling pathway axis. ADAMTS9-AS2 downregulated PPP1R12B by adsorbing miR-196b-5p, so as to regulate the cell cycle signaling pathway to inhibit EC malignant progression.

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.

Differential expression of long non-coding RNAs as diagnostic markers for lung cancer and other malignant tumors

Due to advances in chip and sequencing technology, several types and numbers of long non-coding RNAs (lncRNAs) have been identified. LncRNAs are defined as non-protein-coding RNA molecules longer than 200 nucleotides, and are now thought as a new frontier in the study of human malignant diseases including NSCLC. Diagnosis of numerous malignant tumors has been closely linked to the differential expression of certain lncRNAs. LncRNAs are involved in gene expression regulation at multiple levels of epigenetics, transcriptional regulation, and post-transcriptional regulation. Mutations, deletions, or abnormal expression levels lead to physiological abnormalities, disease occurrence and are closely associated with human tumor diseases. LncRNAs play a crucial role in cancerous processes as either oncogenes or tumor suppressor genes. The expression of lncRNAs can regulate tumor cell in the proliferation, migration, apoptosis, cycle, invasion, and metastasis. As such, lncRNAs are potential diagnostic and treatment targets for cancer. And that, tumor biomarkers need to be detectable in easily accessible body samples, should be characterized by high specificity and sufficient sensitivity. Herein, it is significant clinical importance to screen and supplement new biomarkers for early diagnosis of lung cancer. This study aimed at systematically describing lncRNAs from five aspects based on recent studies: concepts, classification, structure, molecular mechanism, signal pathway, as well as review lncRNA implications in malignant tumor.

Overexpression of cyclin-dependent kinase 1 in esophageal squamous cell carcinoma and its clinical significance

Cyclin‐dependent kinase 1 (CDK1) plays a significant role in certain malignancies. However, it remains unclear whether CDK1 plays a role in esophageal squamous cell carcinoma (ESCC). The aim of this study was to analyze the expression and clinical value of CDK1 in ESCC. CDK1 protein in 151 ESCC tissues and 138 normal esophageal tissues was detected by immunohistochemistry. RNA‐seq of eight pairs of ESCC and adjacent esophageal specimens was performed to evaluate the levels of CDK1 mRNA. Microarray and external RNA‐seq data from 664 cases of ESCC and 1733 cases of control tissues were used to verify the difference in CDK1 expression between the two groups. A comprehensive analysis of all data was performed to evaluate the difference in CDK1 between ESCC tissues and control tissues. Further, functional enrichment analyses were performed based on differentially expressed genes (DEGs) of ESCC and co‐expressed genes (CEGs) of CDK1. In addition, a lncRNA‐miRNA‐CDK1 network was constructed. The expression of CDK1 protein was obviously increased in ESCC tissues (3.540 ± 2.923 vs. 1.040 ± 1.632, P < 0.001). RNA‐seq indicated that the mRNA level of CDK1 was also highly expressed in ESCC tissues (5.261 ± 0.703 vs. 2.229 ± 1.161, P < 0.0001). Comprehensive analysis revealed consistent up‐regulation of CDK1 (SMD = 1.41; 95% CI 1.00–1.83). Further, functional enrichment analyses revealed that the functions of these genes were mainly concentrated in the cell cycle. A triple regulatory network of PVT1‐hsa‐miR‐145‐5p/hsa‐miR‐30c‐5p‐CDK1 was constructed using in silico analysis. In summary, overexpression of CDK1 is closely related to ESCC tumorigenesis.

Upregulated lncRNA HCG18 in Patients with Non-Alcoholic Fatty Liver Disease and Its Regulatory Effect on Insulin Resistance

PURPOSE

Non-alcoholic fatty liver disease (NAFLD) is a disease associated with genetic-environmental-metabolic stress, which severely damages the liver function of patients. This study aimed to explore the significance and probable functions of HCG18 in NAFLD.

PATIENTS AND METHODS

The expression of HCG18 and miR-197-3p was tested by qRT-PCR. The clinical signification of HCG18 was provided by the ROC curve and Pearson correlation. The corresponding mechanism was punctuated by the luciferase reporter assay and HFD-managed mice.

RESULTS

HCG18 expression was higher in the patients with NAFLD than in controls and in individuals with HOMA-IR score ≥2.5 than those with HOMA-IR score <2.5. HCG18 expression in NAFLD patients was related to BMI, HOMA-IR, ALT, FBG, TC, and TG. HCG18 showed satisfactory predictive accuracy in differentiating NAFLD patients and patients with HOMA-IR ≥2.5. Besides, HCG18 had protective impacts on blood glucose and fat deposition but not on body weight. MiR-197-3p is a direct gene of HCG18, and a reverse correlation was found between miR-197-3p and HCG18. Furthermore, miR-197-3p regulated the influence of HCG18 on insulin resistance and lipid accumulation.

CONCLUSION

Increased levels of HCG18 might be an alternate indicator for NAFLD patients. The HCG18-miR-197-3p axis exerted effects on the progression of fat sedimentation and glucose disorder in NAFLD.

ADAMTS9-AS2 and CADM2 expression and association with the prognosis in esophageal squamous cell carcinoma

Background: We investigated whether ADAMTS9-AS2 and CADM2 were related to esophageal squamous cell carcinoma (ESCC). Methodology: ESCC microarray datasets and reverse transcriptase qualitative PCR were used to analyze ADAMTS9-AS2 and CADM2 expression. Results: The GSE120356 and GSE33810 datasets identified ADAMTS9-AS2 and CADM2 as the candidates and ADAMTS9-AS2 and CADM2 expression was downregulated in ESCC. ADAMTS9-AS2 and CADM2 were positively correlated with ESCC. ADAMTS9-AS2 and CADM2 expression could discriminate ESCC from normal tissue. Five-year overall survival was shorter in underexpressed ADAMTS9-AS2 patients, and CADM2 expression level was related to 5-year overall survival. ADAMTS9-AS2 and CADM2 expression were independent prognosis indicators in ESCC patients. Conclusion: Our findings shed new light on the clinical significance of ADAMTS9-AS2 and CADM2 in ESCC carcinogenesis.

Construction of a 14-lncRNA risk score system predicting survival of children with acute myelocytic leukemia

Acute myelocytic leukemia (AML) is a frequent type of acute leukemia. The present study was performed to build a risk score system for the prognostic prediction of AML. AML RNA-sequencing data from samples from 111 children were downloaded from The Cancer Genome Atlas database. Using the DEseq and edgeR packages, the differentially expressed long non-coding RNAs (DE-lncRNAs) between bad and good prognosis groups were identified. A survival package was used to screen prognosis-associated lncRNAs and clinical factors. The optimal lncRNA combination was selected using the penalized package, and the risk-score system was built and evaluated. After the lncRNA-mRNA expression correlation network was constructed, the potential pathways involving the key lncRNAs were enriched using Gene Set Enrichment Analysis. Among the 61 DE-lncRNAs, 48 lncRNAs were significantly associated with prognosis. Relapse was an independent prognostic factor. The optimal 14-lncRNA risk score system was constructed. After 730 differentially expressed mRNAs were identified between the good and bad prognosis groups divided using a prognostic index, the lncRNA-mRNA expression correlation network was constructed. Enrichment analysis showed that semaphorin-3C [SEMA3C; regulated by probable leucine-tRNA ligase, mitochondrial (LARS2-AS1)] and secreted frizzled-related protein 5 [SFRP5; mediated by WASH complex subunit 5 (WASHC5)-antisense RNA 1 (AS1)] were involved in axon guidance and the Wnt signaling pathway, respectively. A 14-lncRNA (including paired box protein Pax8-AS1 and MYB AS1) risk-score system might be effective in predicting the prognosis of AML. Axon guidance (involving SEMA3C and LARS2-AS1) and the Wnt signaling pathway (involving SFRP5 and WASHC5-AS1) might be two important pathways affecting the prognosis of AML.

Synthesis and Biological Evaluation of HDAC Inhibitors With a Novel Zinc Binding Group

Vorinostat (SAHA) with great therapeutic potential has been approved by the FDA for the treatment of cutaneous T-cell lymphoma as the first HDACs inhibitor, but the drawbacks associated with hydroxamic acid group (poor stability, easy metabolism, weak binding ability to class IIa isozymes, and poor selectivity) have been exposed during the continuous clinical application. Based on the pharmacophore of HDAC inhibitors, two series of compounds with novel zinc binding group (ZBG) were designed and synthesized, and the antitumor bioactivities were evaluated in four human cancer cell lines (A549, Hela, HepG2, and MCF-7). Among the synthesized compounds, compounds a6, a9, a10, b8, and b9 exhibited promising inhibitory activities against the selected tumor cell lines, especially compounds a9 and b8 on Hela's cytostatic activity (a9: IC₅₀ = 11.15 ± 3.24 μM; b8: IC₅₀ = 13.68 ± 1.31 μM). The enzyme inhibition assay against Hela extracts and HDAC1&6 subtypes showed that compound a9 had a certain broad-spectrum inhibitory activity, while compound b8 had selective inhibitory activity against HDAC6, which was consistent with Western blot results. In addition, the inhibitory mechanism of compounds a9 and b8 in HDAC1&6 were both compared through computational approaches, and the binding interactions between the compounds and the enzymes target were analyzed from the perspective of energy profile and conformation. In summary, the compounds with novel ZBG exhibited certain antitumor activities, providing valuable hints for the discovery of novel HDAC inhibitors.

The cox-filter method identifies respective subtype-specific lncRNA prognostic signatures for two human cancers

BACKGROUND

The most common histological subtypes of esophageal cancer are squamous cell carcinoma (ESCC) and adenocarcinoma (EAC). It has been demonstrated that non-marginal differences in gene expression and somatic alternation exist between these two subtypes; consequently, biomarkers that have prognostic values for them are expected to be distinct. In contrast, laryngeal squamous cell cancer (LSCC) has a better prognosis than hypopharyngeal squamous cell carcinoma (HSCC). Likewise, subtype-specific prognostic signatures may exist for LSCC and HSCC. Long non-coding RNAs (lncRNAs) hold promise for identifying prognostic signatures for a variety of cancers including esophageal cancer and head and neck squamous cell carcinoma (HNSCC).

METHODS

In this study, we applied a novel feature selection method capable of identifying specific prognostic signatures uniquely for each subtype - the Cox-filter method - to The Cancer Genome Atlas esophageal cancer and HSNCC RNA-Seq data, with the objectives of constructing subtype-specific prognostic lncRNA expression signatures for esophageal cancer and HNSCC.

RESULTS

By incorporating biological relevancy information, the lncRNA lists identified by the Cox-filter method were further refined. The resulting signatures include genes that are highly related to cancer, such as H19 and NEAT1, which possess perfect prognostic values for esophageal cancer and HNSCC, respectively.

CONCLUSIONS

The Cox-filter method is indeed a handy tool to identify subtype-specific prognostic lncRNA signatures. We anticipate the method will gain wider applications.

Integrated analysis highlights multiple long non‑coding RNAs and their potential roles in the progression of human esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is a prevalent aggressive malignant tumor with poor prognosis. Investigations into the molecular changes that occur as a result of the disease, as well as identification of novel biomarkers for its diagnosis and prognosis, are urgently required. Long non‑coding RNAs (lncRNAs) have been reported to play a critical role in tumor progression. The present study performed data mining analyses for ESCC via an integrated study of accumulated datasets and identification of the differentially expressed lncRNAs from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. The identified intersection of differentially expressed genes (lncRNAs, miRNAs and mRNAs) in ESCC tissues between the GEO and TCGA datasets was investigated. Based on these intersected lncRNAs, the present study constructed a competitive endogenous RNA (ceRNA) network of lncRNAs in ESCC. A total of 81 intersection lncRNAs were identified; 67 of these were included in the ceRNA network. Functional analyses revealed that these 67 key lncRNAs primarily dominated cellular biological processes. The present study then analyzed the associations between the expression levels of these 67 key lncRNAs and the clinicopathological characteristics of the ESCC patients, as well as their survival time using TCGA. The results revealed that 31 of these lncRNAs were associated with tumor grade, tumor‑node‑metastasis (TNM) stage and lymphatic metastasis status (P<0.05). In addition, 15 key lncRNAs were demonstrated to be associated with survival time (P<0.05). Finally, 5 key lncRNAs were selected for validation of their expression levels in 30 patients newly diagnosed with ESCC via reverse transcription‑quantitative PCR (RT‑qPCR). The results suggested that the fold changes in the trends of up‑ and downregulation between GEO, TCGA and RT‑qPCR were consistent. In addition, it was also demonstrated that a select few of these 5 key lncRNAs were significantly associated with TNM stage and lymph node metastasis (P<0.05). The results of the clinically relevant analysis and the aforementioned bioinformatics were similar, hence proving that the bioinformatics analysis used in the present study is credible. Overall, the results from the present study may provide further insight into the functional characteristics of lncRNAs in ESCC through bioinformatics integrative analysis of the GEO and TCGA datasets, and reveal potential diagnostic and prognostic biomarkers for ESCC.