Decoding competing endogenous RNA networks for cancer biomarker discovery

LINC01089 in cancer: multifunctional roles and therapeutic implications

LINC01089 is a prime example of a long non-coding RNA that plays a pivotal role in the progression of human cancers. The gene encoding this lncRNA is located on 12q24.31. LINC01089 has been demonstrated to exert tumor-suppressive effects in various cancers, including colorectal cancer, gastric cancer, lung cancer, ovarian cancer, cervical cancer, papillary thyroid carcinoma, breast cancer, and osteosarcoma. However, its role in hepatocellular carcinoma shows significant discrepancies across different studies. In this review, we systematically explore the functions of LINC01089 in human cancers through bioinformatics analysis, clinical studies, animal models, and fundamental experimental research. Furthermore, we delve into the biological mechanisms and functions of LINC01089, and discuss its potential as a future biomarker and therapeutic target in detail.

The Traces of Dysregulated lncRNAs-Associated ceRNA Axes in Retinoblastoma: A Systematic Scope Review

PURPOSE

Long non-coding RNAs are an essential component of competing endogenous RNA regulatory axes and play their role by sponging microRNAs and interfering with the regulation of gene expression. Because of the broadness of competing endogenous RNA interaction networks, they may help investigate treatment targets in complicated disorders.

METHODS

This study performed a systematic scoping review to assess verified loops of competing endogenous RNAs in retinoblastoma, emphasizing the competing endogenous RNAs axis related to long non-coding RNAs. We used a six-stage approach framework and the PRISMA guidelines. A systematic search of seven databases was done to locate suitable papers published before February 2022. Two reviewers worked independently to screen articles and collect data.

RESULTS

Out of 363 records, fifty-one articles met the inclusion criteria, and sixty-three axes were identified in desired articles. The majority of the research reported several long non-coding RNAs that were experimentally verified to act as competing endogenous RNAs in retinoblastoma: XIST/NEAT1/MALAT1/SNHG16/KCNQ1OT1, respectively. At the same time, around half of the studies investigated unique long non-coding RNAs.

CONCLUSIONS

Understanding the many features of this regulatory system may aid in elucidating the unknown etiology of Retinoblastoma and providing novel molecular targets for therapeutic and clinical applications.

Tweety homolog 3 promotes colorectal cancer progression through mutual regulation of histone deacetylase 7

Colorectal cancer (CRC) is one of the leading cancers worldwide, with metastasis being a major cause of high mortality rates among patients. In this study, dysregulated gene Tweety homolog 3 (TTYH3) was identified by Gene Expression Omnibus database. Public databases were used to predict potential competing endogenous RNAs (ceRNAs) for TTYH3. Quantitative real-time polymerase chain reaction, western blot, and immunohistochemistry were utilized to analyze TTYH3 and histone deacetylase 7 (HDAC7) levels. Luciferase assays confirmed miR-1271-5p directly targeting the 3' untranslated regions of TTYH3 and HDAC7. In vitro experiments such as transwell and human umbilical vein endothelial cell tube formation, as well as in vivo mouse models, were conducted to assess the biological functions of TTYH3 and HDAC7. We discovered that upregulation of TTYH3 in CRC promotes cell migration by affecting the Epithelial-mesenchymal transition pathway, which was independent of its ion channel activity. Mechanistically, TTYH3 and HDAC7 functioned as ceRNAs, reciprocally regulating each other's expression. TTYH3 competes for binding miR-1271-5p, increasing HDAC7 expression, facilitating CRC metastasis and angiogenesis. This study reveals the critical role of TTYH3 in promoting CRC metastasis through ceRNA crosstalk, offering new insights into potential therapeutic targets for clinical intervention.

Long Non-Coding RNAs as Determinants of Thyroid Cancer Phenotypes: Investigating Differential Gene Expression Patterns and Novel Biomarker Discovery

Thyroid Cancer (TC) is the most common endocrine malignancy, with increasing incidence globally. Papillary thyroid cancer (PTC), a differentiated form of TC, accounts for approximately 90% of TC and occurs predominantly in women of childbearing age. Although responsive to current treatments, recurrence of PTC by middle age is common and is much more refractive to treatment. Undifferentiated TC, particularly anaplastic thyroid cancer (ATC), is the most aggressive TC subtype, characterized by it being resistant and unresponsive to all therapeutic and surgical interventions. Further, ATC is one of the most aggressive and lethal malignancies across all cancer types. Despite the differences in therapeutic needs in differentiated vs. undifferentiated TC subtypes, there is a critical unmet need for the identification of molecular biomarkers that can aid in early diagnosis, prognosis, and actionable therapeutic targets for intervention. Advances in the field of cancer genomics have enabled for the elucidation of differential gene expression patterns between tumors and healthy tissue. A novel category of molecules, known as non-coding RNAs, can themselves be differentially expressed, and extensively contribute to the up- and downregulation of protein coding genes, serving as master orchestrators of regulated and dysregulated gene expression patterns. These non-coding RNAs have been identified for their roles in driving carcinogenic patterns at various stages of tumor development and have become attractive targets for study. The identification of specific genes that are differentially expressed can give insight into mechanisms that drive carcinogenic patterns, filling the gaps of deciphering molecular and cellular processes that modulate TC subtypes, outside of well-known driver mutations.

Machine Learning Gene Signature to Metastatic ccRCC Based on ceRNA Network

Clear-cell renal-cell carcinoma (ccRCC) is a silent-development pathology with a high rate of metastasis in patients. The activity of coding genes in metastatic progression is well known. New studies evaluate the association with non-coding genes, such as competitive endogenous RNA (ceRNA). This study aims to build a ceRNA network and a gene signature for ccRCC associated with metastatic development and analyze their biological functions. Using data from The Cancer Genome Atlas (TCGA), we constructed the ceRNA network with differentially expressed genes, assembled nine preliminary gene signatures from eight feature selection techniques, and evaluated the classification metrics to choose a final signature. After that, we performed a genomic analysis, a risk analysis, and a functional annotation analysis. We present an 11-gene signature: SNHG15, AF117829.1, hsa-miR-130a-3p, hsa-mir-381-3p, BTBD11, INSR, HECW2, RFLNB, PTTG1, HMMR, and RASD1. It was possible to assess the generalization of the signature using an external dataset from the International Cancer Genome Consortium (ICGC-RECA), which showed an Area Under the Curve of 81.5%. The genomic analysis identified the signature participants on chromosomes with highly mutated regions. The hsa-miR-130a-3p, AF117829.1, hsa-miR-381-3p, and PTTG1 were significantly related to the patient's survival and metastatic development. Additionally, functional annotation resulted in relevant pathways for tumor development and cell cycle control, such as RNA polymerase II transcription regulation and cell control. The gene signature analysis within the ceRNA network, with literature evidence, suggests that the lncRNAs act as "sponges" upon the microRNAs (miRNAs). Therefore, this gene signature presents coding and non-coding genes and could act as potential biomarkers for a better understanding of ccRCC.

CircZCCHC2 decreases pirarubicin sensitivity and promotes triple-negative breast cancer development via the miR-1200/TPR axis

Triple-negative breast cancer (TNBC) has attracted attention due to its poor prognosis and limited treatment options. The mechanisms underlying the association between circular RNAs (circRNAs) and the occurrence and development of TNBC remain unclear. CircZCCHC2 is observed to be upregulated in TNBC cells, tissues, and plasma exosomes. Knockdown of circZCCHC2 inhibited the proliferation, migration, invasion, and epithelial-mesenchymal transition of TNBC cells in vitro and in vivo. Pirarubicin (THP) treatment downregulated circZCCHC2, and circZCCHC2 affected the sensitivity to THP. CircZCCHC2/miR-1200/translocated promoter region, the nuclear basket protein (TPR) pathway was cascaded and verified. It is demonstrated that circZCCHC2 plays a crucial role in the malignant progression of TNBC via the miR-1200/TPR axis, thereby activating the RAS-RAF-MEK-ERK pathway. The present results indicate that circZCCHC2 has the potential to serve as a novel prognostic biomarker for TNBC.

Hsa_circ_0007990 promotes breast cancer growth via inhibiting YBX1 protein degradation to activate E2F1 transcription

Breast cancer (BC) is the most commonly diagnosed malignant tumour in females worldwide. Although remarkable advances in early detection and treatment strategies have led to decreased mortality, recurrence and metastasis remain the major causes of cancer death in BC patients. Increasing evidence has demonstrated that circular RNAs (circRNAs) play critical roles in cancer progression. However, the detailed biological functions and molecular mechanisms of circRNAs in BC are unclear. The aim of this study was to investigate the possible role of circRNAs in the progression of BC. Differentially expressed circRNAs in BC were identified by integrating breast tumour-associated somatic CNV data and circRNA high-throughput sequencing. Aberrant hsa_circ_0007990 expression and host gene copy number were detected in BC cell lines via quantitative polymerase chain reaction (qPCR). The expression level of hsa_circ_0007990 in BC tissues was validated by in situ hybridization (ISH). Loss- and gain-of-function experiments were performed in vitro and in vivo, respectively, to explore the potential biological function of hsa_circ_0007990 in BC. The underlying mechanisms of hsa_circ_0007990 were investigated through MS2 RNA pull-down, RNA immunoprecipitation, RNA fluorescence in situ hybridization, immunofluorescence, chromatin immunoprecipitation and luciferase reporter assays. The levels of hsa_circ_0007990 were elevated in BC tissues and cell lines, an effect that was partly due to host gene copy number gains. Functional assays showed that hsa_circ_0007990 promoted BC cell growth. Mechanistically, hsa_circ_0007990 could bind to YBX1 and inhibit its degradation by preventing ubiquitin/proteasome-dependent degradation, thus enhancing the expression of the cell cycle-associated gene E2F1. Rescue experiments suggested that hsa_circ_0007990 promoted BC progression through YBX1. In general, our study demonstrated that hsa_circ_0007990 modulates the ubiquitination and degradation of YBX1 protein and further regulates E2F1 expression to promote BC progression. We explored the possible function and molecular mechanism of hsa_circ_0007990 in BC and identified a novel candidate target for the treatment of BC.

The anti-hepatocellular carcinoma effect of Aidi injection was related to the synergistic action of cantharidin, formononetin, and isofraxidin through BIRC5, FEN1, and EGFR

ETHNOPHARMACOLOGICAL RELEVANCE

Aidi injection (ADI) is a popular anti-tumor Chinese patent medicine, widely used in clinics for the treatment of hepatocellular carcinoma (HCC) with remarkable therapeutic effects through multiple targets and pathways. However, the scientific evidence of the synergistic role of the complex chemical component system and the potential mechanism for treating diseases are ignored and remain to be elucidated.

AIM OF THE STUDY

This study aimed to elucidate and verify the cooperative association between the potential active ingredient of ADI, which is of significance to enlarge our understanding of its anti-HCC molecular mechanisms.

MATERIALS AND METHODS

Firstly, the anti-HCC effect of ADI was evaluated in various HCC cells and the zebrafish xenograft model. Subsequently, a variety of bioinformatic technologies, including network pharmacology, weighted gene co-expression network analysis (WGCNA), meta-analysis of gene expression profiles, and pathway enrichment analysis were performed to construct the competitive endogenous RNA (ceRNA) network of ADI intervention in HCC and to establish the relationship between the critical targets/pathways and the key corresponding components, which were involved in ADI against HCC in a synergistic way and were validated by molecular biology experiments.

RESULTS

ADI exerted remarkable anti-HCC in vitro cells and in vivo zebrafish model, especially that the Hep 3B2.1-7 cell showed substantial sensibility to ADI. The ceRNA network revealed that the EGFR/PI3K/AKT signaling pathway was identified as the promising pathway. Furthermore, the meta-analysis also demonstrated the critical role of BIRC5 and FEN1 as key targets. Finally, the synergistic effect of ADI was revealed by discovering the inhibitory effect of cantharidin on BIRC5, formononetin on FEN1 and EGFR, as well as isofraxidin on EGFR.

CONCLUSION

Our study unveiled that the incredible protective effect of ADI on HCC resulted from the synergistic inhibition effect of cantharidin, formononetin, and isofraxidin on multiple targets/pathways, including BIRC5, FEN1, and EGFR/PI3K/AKT, respectively, providing a scientific interpretation of ADI against HCC and a typical example of pharmacodynamic evaluation of other proprietary Chinese patent medicine.

Comprehensive analysis of lncRNA-miRNA-mRNA ceRNA network and key genes in granulosa cells of patients with biochemical primary ovarian insufficiency

Primary ovarian insufficiency (POI) is a common condition leading to the pathological decline of ovarian function in women of reproductive age, resulting in amenorrhea, hypogonadism, and infertility. Biochemical premature ovarian insufficiency (bPOI) is an intermediate stage in the pathogenesis of POI in which the fertility of patients has been reduced. Previous studies suggest that granulosa cells (GCs) play an essential role in the pathogenesis of POI, but their pathogenetic mechanisms remain unclear. To further explore the potential pathophysiological mechanisms of GCs in POI, we constructed a molecular long non-coding RNA (lncRNA)-microRNA (miRNA)-messenger RNA (mRNA) network using GC expression data collected from biochemical premature ovarian failure (bPOI) patients in the GEO database. We discovered that the GCs of bPOI patients had differential expression of 131 mRNAs, 191 lncRNAs, and 28 miRNAs. By systematic network analysis, we identified six key genes, including SRSF1, PDIA5, NEURL1B, UNK, CELF2, and CFL2, and five hub miRNAs, namely hsa-miR-27a-3p, hsa-miR-24-3p, hsa-miR-22-3p, hsa-miR-129-5p, and hsa-miR-17-5p, and the results suggest that the expression of these key genes may be regulated by two hub miRNAs, hsa-miR-27a-3p and hsa-miR-17-5p. Additionally, a POI model in vitro was created to confirm the expression of a few important genes. In this study, we discovered a unique lncRNA-miRNA-mRNA network based on the ceRNA mechanism in bPOI for the first time, and we screened important associated molecules, providing a partial theoretical foundation to better understand the pathogenesis of POI.

A Five-gene Signature based on MicroRNA for Predicting Prognosis and Immunotherapy in Stomach Adenocarcinoma

AIMS

We aimed to classify molecular subtypes and establish a prognostic gene signature based on miRNAs for the prognostic prediction and therapeutic response in Stomach adenocarcinoma (STAD).

BACKGROUND

STAD is a common diagnosed gastrointestinal malignancy and its heterogeneity is a big challenge that influences prognosis and precision therapies. Present study was designed to classify molecular subtypes and construct a prognostic gene signature based on miRNAs for the prognostic prediction and therapeutic response in STAD.

OBJECTIVE

The objective of this study is to investigate the molecular subtypes and prognostic model for STAD.

METHODS

A STAD specific miRNA-messenger RNA (mRNA) competing endogenous RNA (ceRNA) network was generated using the RNA-Seq and miRNA expression profiles from The Cancer Genome Atlas (TCGA) database, in which miRNA-related mRNAs were screened. Molecular subtypes were then determined using miRNA-related genes. Through univariate Cox analysis and multivariate regression analysis, a prognostic model was established in GSE84437 Train dataset and validated in GSE84437 Test, TCGA, GSE84437 and GSE66229 datasets. Immunotherapy datasets were employed for assessing the performance of the risk model. Finally, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was applied to validate the expression of hub genes used for the risk score signature.

RESULTS

We constructed a ceRNA network containing 84 miRNAs and 907 mRNAs and determined two molecular subtypes based on 26 genes from the intersection of TCGASTAD and GSE84437 datasets. Subtype S2 had poor prognosis, lower tumor mutational burden, higher immune score and lower response to immunotherapy. Subtype S1 was more sensitive to Sorafenib, Pyrimethamine, Salubrinal, Gemcitabine, Vinorelbine and AKT inhibitor VIII. Next, a five-gene signature was generated and its robustness was validated in Test and external datasets. This risk model also had a good prediction performance in immunotherapy datasets.

CONCLUSION

This study promotes the underlying mechanisms of miRNA-based genes in STAD and offers directions for classification. A five-gene signature accurately predicts the prognosis and helps therapeutic options.

Identification of AP002498.1 and LINC01871 as prognostic biomarkers and therapeutic targets for distant metastasis of colorectal adenocarcinoma

BACKGROUND

Increasing evidence suggests that lncRNA (Long non-coding RNA, lncRNA)-mediated ceRNA (competing endogenous RNA, ceRNA) networks are involved in the occurrence and progression of colorectal cancer (CRC). However, the roles of the lncRNA-miRNA-mRNA ceRNA network in distant metastasis of CRC are still unclear.

METHODS

In this study, we constructed a specific ceRNA network to identify potential biomarkers and therapeutic targets for distant metastasis of CRC. Specifically, RNA-Seq data from The Cancer Genome Atlas (TCGA) were used to screen for differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) related to metastasis. After validation and selection by qRT-PCR and univariate and multivariate analysis of the metastasis- and prognosis-related lncRNAs, the regulated microRNAs (miRNAs) and coexpressed mRNAs were used to construct a ceRNA network for distant metastasis of CRC.

RESULTS

Two key distant metastasis-related DElncRNAs, AP002498.1 and LINC01871, were identified by univariate and multivariate analysis in combination with analyses of clinical data and expression levels. Furthermore, lncRNA-associated ceRNA subnetworks were constructed from the predicted miRNAs and 13 coexpressed DEmRNAs (SERPINA1, ITLN1, REG4, L1TD1, IGFALS, MUC5B, CIITA, CXCL9, CXCL10, GBP4, GNLY, IDO1, and NOS2). The AP002498.1- and LINC01871-associated ceRNA subnetworks regulated the expression of the target genes SERPINA1 and MUC5B and GNLY, respectively, through the associated miRNAs.

CONCLUSION

The DElncRNA AP002498.1 and the LINC01871/miR-4644 and miR-185-5p/GNLY axes were identified as being closely associated with distant metastasis and could represent independent prognostic biomarkers or therapeutic targets in colorectal adenocarcinoma.

New understandings of the genetic regulatory relationship between non-coding RNAs and m6A modification

One of the most frequent epigenetic modifications of RNA in eukaryotes is N6 methyladenosine (m6A), which is mostly present in messenger RNAs. Through the influence of several RNA processing stages, m6A modification is a crucial approach for controlling gene expression, especially in cancer progression. It is universally acknowledged that numerous non-coding RNAs (ncRNAs), such as microRNAs, circular RNAs, long non-coding RNAs, and piRNAs, are also significantly affected by m6A modification, and the complex genetic regulatory relationship between m6A and ncRNAs plays a pivotal role in the development of cancer. The connection between m6A modifications and ncRNAs offers an opportunity to explore the oncogene potential regulatory mechanisms and suggests that m6A modifications and ncRNAs could be vital biomarkers for multiple cancers. In this review, we discuss the mechanisms of interaction between m6A methylation and ncRNAs in cancer, and we also summarize diagnostic and prognostic biomarkers for clinical cancer detection. Furthermore, our article includes some methodologies for identifying m6A sites when assessing biomarker potential.

Construction of lncRNA-miRNA-mRNA regulatory network in severe asthmatic bronchial epithelial cells: A bioinformatics study

Asthma is a chronic respiratory disease caused by environment-host interactions. Bronchial epithelial cells (BECs) are the first line of defense against environmental toxins. However, the mechanisms underlying the role of BECs in severe asthma (SA) are not yet fully understood. Long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) have been shown to play important roles in the regulation of gene expression in the pathogenesis of SA. In this study, bioinformatics was used for the first time to reveal the lncRNA-miRNA-mRNA regulatory network of BECs in SA. Five mRNA datasets of bronchial brushing samples from patients with SA and healthy controls (HC) were downloaded from the Gene Expression Omnibus (GEO) database. A combination of the Venn diagram and robust rank aggregation (RRA) method was used to identify core differentially expressed genes (DEGs). Protein-protein interaction (PPI) analysis of core DEGs was performed to screen hub genes. The miRDB, miRWalk, and ENCORI databases were used to predict the miRNA-mRNA relationships, and the ENCORI and starBase v2.0 databases were used to predict the upstream lncRNAs of the miRNA-mRNA relationships. Four core DEGs were identified: carcinoembryonic antigen-related cell adhesion molecule 5 (CEACAM5), interleukin-1 receptor type 2 (IL1R2), trefoil factor 3 (TFF3), and vascular endothelial growth factor A (VEGFA). These 4 core DEGs indicated that SA was not significantly associated with sex. Enrichment analysis showed that the MAPK, Rap1, Ras, PI3K-Akt and Calcium signaling pathways may serve as the principal pathways of BECs in SA. A lncRNA-miRNA-mRNA regulatory network of the severe asthmatic bronchial epithelium was constructed. The top 10 competing endogenous RNAs (ceRNAs) were FGD5 antisense RNA 1 (FGD5-AS1), metastasis associated lung adenocarcinoma transcript 1 (MALAT1), X inactive specific transcript (XIST), HLA complex group 18 (HCG18), small nucleolar RNA host gene 16 (SNHG16), has-miR-20b-5p, has-miR-106a-5p, hsa-miR-106b-5p, has-miR-519d-3p and Fms related receptor tyrosine kinase 1 (FLT1). Our study revealed a potential mechanism for the lncRNA-miRNA-mRNA regulatory network in BECs in SA.

RNA-RNA competitive interactions: a molecular civil war ruling cell physiology and diseases

The idea that proteins are the main determining factors in the functioning of cells and organisms, and their dysfunctions are the first cause of pathologies, has been predominant in biology and biomedicine until recently. This protein-centered view was too simplistic and failed to explain the physiological and pathological complexity of the cell. About 80% of the human genome is dynamically and pervasively transcribed, mostly as non-protein-coding RNAs (ncRNAs), which competitively interact with each other and with coding RNAs generating a complex RNA network regulating RNA processing, stability, and translation and, accordingly, fine-tuning the gene expression of the cells. Qualitative and quantitative dysregulations of RNA-RNA interaction networks are strongly involved in the onset and progression of many pathologies, including cancers and degenerative diseases. This review will summarize the RNA species involved in the competitive endogenous RNA network, their mechanisms of action, and involvement in pathological phenotypes. Moreover, it will give an overview of the most advanced experimental and computational methods to dissect and rebuild RNA networks.

Novel insights on perils and promises of miRNA in understanding colon cancer metastasis and progression

Colorectal cancer (CRC) is the third highest frequent malignancy and ultimate critical source of cancer-associated mortality around the world. Regardless of latest advances in molecular and surgical targeted medicines that have increased remedial effects in CRC patients, the 5-year mortality rate for CRC patients remains dismally low. Evidence suggests that microRNAs (miRNAs) execute an essential part in the development and spread of CRC. The miRNAs are a type of short non-coding RNA that exhibited to control the appearance of tumor suppressor genes and oncogenes. miRNA expression profiling is already being utilized in clinical practice as analytical and prognostic biomarkers to evaluate cancer patients' tumor genesis, advancement, and counteraction to drugs. By modulating their target genes, dysregulated miRNAs are linked to malignant characteristics (e.g., improved proliferative and invasive capabilities, cell cycle aberration, evasion of apoptosis, and promotion of angiogenesis). This review presents an updated summary of circulatory miRNAs, tumor-suppressive and oncogenic miRNAs, and the potential reasons for dysregulated miRNAs in CRC. Further we will explore the critical role of miRNAs in CRC drug resistance.

Crosstalk between long noncoding RNA and microRNA in Cancer

miRNAs and lncRNAs play a central role in cancer-associated gene regulations. The dysregulated expression of lncRNAs has been reported as a hallmark of cancer progression, acting as an independent prediction marker for an individual cancer patient. The interplay of miRNA and lncRNA decides the variation of tumorigenesis that could be mediated by acting as sponges for endogenous RNAs, regulating miRNA decay, mediating intra-chromosomal interactions, and modulating epigenetic components. This paper focuses on the influence of crosstalk between lncRNA and miRNA on cancer hallmarks such as epithelial-mesenchymal transition, hijacking cell death, metastasis, and invasion. Other cellular roles of crosstalks, such as neovascularization, vascular mimicry, and angiogenesis were also discussed. Additionally, we reviewed crosstalk mechanism with specific host immune responses and targeting interplay (between lncRNA and miRNA) in cancer diagnosis and management.

CircEGFR reduces the sensitivity of pirarubicin and regulates the malignant progression of triple-negative breast cancer via the miR-1299/EGFR axis

Circular RNAs (circRNAs) have been found to be involved in cancer progression and chemotherapy sensitivity. However, the biological function of circRNAs in triple-negative breast cancer (TNBC) and its effect on the sensitivity to pirarubicin (THP) chemotherapy are still unclear. CircEGFR (hsa_circ_0080220) was screened and verified by bioinformatics analysis, proving it was highly expressed in TNBC cell lines, patient tissues, and plasma exosomes, and was associated with poor prognosis of patients. The expression level of circEGFR in patient tissue has potential diagnostic value to distinguish TNBC tissue from normal breast tissue. In vitro studies confirmed that overexpression of circEGFR promoted the proliferation, migration, invasion, and EMT of TNBC cells and decreased the sensitivity of THP treatment while silencing circEGFR showed the opposite effect. The circEGFR/miR-1299/EGFR pathway was cascaded and verified. CircEGFR regulated malignant progression of TNBC by regulating EGFR via sponging miR-1299. THP can inhibit the malignant phenotype of MDA-MB-231 cells by downregulating the expression of circEGFR. In vivo studies confirmed that overexpression of circEGFR can promote tumor growth and EMT and reduce tumor sensitivity to THP treatment. Silencing circEGFR inhibited the malignant progression of the tumor. These results revealed circEGFR is a promising biomarker for TNBC diagnosis, therapeutic and prognosis.

MiR-200/183 family-mediated module biomarker for gastric cancer progression: an AI-assisted bioinformatics method with experimental functional survey

BACKGROUND

Gastric cancer (GC) is a major cancer burden throughout the world with a high mortality rate. The performance of current predictive and prognostic factors is still limited. Integrated analysis is required for accurate cancer progression predictive biomarker and prognostic biomarkers that help to guide therapy.

METHODS

An AI-assisted bioinformatics method that combines transcriptomic data and microRNA regulations were used to identify a key miRNA-mediated network module in GC progression. To reveal the module's function, we performed the gene expression analysis in 20 clinical samples by qRT-PCR, prognosis analysis by multi-variable Cox regression model, progression prediction by support vector machine, and in vitro studies to elaborate the roles in GC cells migration and invasion.

RESULTS

A robust microRNA regulated network module was identified to characterize GC progression, which consisted of seven miR-200/183 family members, five mRNAs and two long non-coding RNAs H19 and CLLU1. Their expression patterns and expression correlation patterns were consistent in public dataset and our cohort. Our findings suggest a two-fold biological potential of the module: GC patients with high-risk score exhibited a poor prognosis (p-value < 0.05) and the model achieved AUCs of 0.90 to predict GC progression in our cohort. In vitro cellular analyses shown that the module could influence the invasion and migration of GC cells.

CONCLUSIONS

Our strategy which combines AI-assisted bioinformatics method with experimental and clinical validation suggested that the miR-200/183 family-mediated network module as a "pluripotent module", which could be potential marker for GC progression.

LINC00511, a future star for the diagnosis and therapy of digestive system malignant tumors

The digestive system malignant tumors (DSMTs), mainly consist of digestive tract and digestive gland tumors, become an inescapable culprit to hazard human health worldwide. Due to the huge hysteresis in the cognitive theories of DSMTs occurrence and progression, advances in medical technology have not improved the prognosis. Therefore, more studies on a variety of tumor-associated molecular biomarkers and more detailed disclosure on potential regulatory networks are urgently needed to facilitate the diagnostic and therapeutic strategies of DSMTs. With the development of cancer bioinformatics, a special type of endogenous RNA involved in multi-level cellular function regulation rather than encoding protein, is categorized as non-coding RNAs (ncRNAs) and becomes a hotspot issue in oncology. Among them, long non-coding RNAs (lncRNAs), transcription length > 200 nt, show obvious superiority in both research quantity and dimension compared to microRNAs (miRNAs) and circular RNAs (circRNAs). As a recently discovered lncRNA, LINC00511 has been confirmed to be closely associated with DSMTs and might be exploited as a novel biomarker. In the present review, the comprehensive studies of LINC00511 in DSMTs are summarized, as well as the underlying molecular regulatory networks. In addition, deficiencies in researches are point out and discussed. The Cumulative oncology studies provide a fully credible theoretical basis for identifying the regulatory role of LINC00511 in human DSMTs. LINC00511, proved to be an oncogene in DSMTs, might be defined as a potential biomarker for diagnosis and prognosis evaluation, as well as a rare therapeutic target.

Hsa_circ_0000129 drives tumor growth via sequestering miR-485-3p and upregulating SPIN1 in breast cancer

BACKGROUND

Breast cancer (BC) is second cancer frequently occurring worldwide. Circular RNA hsa_circ_0000129 (circ_0000129) exerts a tumor-promoting effect in BC. Nevertheless, the molecular mechanisms mediated by the upregulation of circ_0000129 during BC progression are not well understood.

METHODS

Forty-five BC patients were recruited for the research. Changes in circ_0000129 levels were detected with quantitative reverse transcription-polymerase chain reaction. Cell proliferation, apoptosis, migration, invasion, and angiopoiesis were determined by cell counting, 5-ethynyl-2'-deoxyuridine (EdU), flow cytometry, transwell, and tube formation assays. Protein levels were detected by western blot analysis. The regulatory mechanism of circ_0000129 was predicted by bioinformatics analysis and validated by dual-luciferase reporter and RNA immunoprecipitation assays. In vivo experiments were carried out to verify the function of circ_0000129.

RESULTS

Circ_0000129 was overexpressed in BC samples and cell lines. Functionally, circ_0000129 silencing reduced cell proliferation, migration, invasion, and promoted cell apoptosis, as well as induced HUVEC angiopoiesis in vitro. Furthermore, circ_0000129 knockdown decreased BC cell growth in mouse xenograft models. Mechanically, circ_0000129 interacted with miR-485-3p to mediate the inhibiting effect of miR-485-3p on SPIN1. Silenced miR-485-3p expression weakened the inhibiting effect of circ_0000129 knockdown on BC cell malignant behaviors. Also, forced SPIN1 expression weakened miR-485-3p upregulation mediated effects on BC cell malignant behaviors.

CONCLUSION

Circ_0000129 acted as a miR-485-3p sponge molecular to mediate expression, thus promoting BC progression.

Identification of a key ceRNA network associated with ferroptosis in gastric cancer

Ferroptosis, a newly discovered irondependent form of regulated cell death caused by excessive accumulation of lipid peroxides, is linked to the development and treatment response of various types of cancer, including gastric cancer (GC). Noncoding RNAs (ncRNAs), as key regulators in cancer, have both oncogenic and tumor suppressive roles. However, studies on ferroptosis-related ncRNA networks in GC are still lacking. Here, we first identified 61 differentially expressed genes associated with ferroptosis in GC by computing and analyzing gene expression profile of tumor and normal tissues for GC. Then, upstream lncRNAs and miRNAs interacting with them were found through miRNet and miRBase databases, and hub lncRNAs and miRNAs were obtained through topological analysis. Finally, the ceRNA regulatory network linked to ferroptosis in GC was established, which includes two ferroptosis marker genes (TXNIP and TSC22D3), one driver gene (GABARAPL1), and one suppressor gene (CAV1). Kaplan-Meier survival analysis showed that changes in the expression of these genes were associated with the survival of GC patients. Furthermore, our study revealed that this ceRNA network may influence the progression of GC by regulating ferroptosis process. These results will help experimental researchers to design an experiment study to further explore the roles of this regulatory network in GC ferroptosis.

The roles of lncRNA functions and regulatory mechanisms in the diagnosis and treatment of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most frequent and deadly type of liver cancer. While the underlying molecular mechanisms are poorly understood, it is documented that lncRNAs may play key roles. Many HCC-associated lncRNAs have been linked to HBV and HCV infection, mediating gene expression, cell growth, development, and death. Studying the regulatory mechanisms and biological functions of HCC-related lncRNAs will assist our understanding of HCC pathogenesis as well as its diagnosis and management. Here, we address the potential of dysregulated lncRNAs in HCC as diagnostic and therapeutic biomarkers, and we evaluate the oncogenic or tumor-suppressive properties of these lncRNAs.

Identification of Hypoxia-Related Prognostic Signature and Competing Endogenous RNA Regulatory Axes in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a common type of liver cancer and one of the highly lethal diseases worldwide. Hypoxia plays an important role in the development and prognosis of HCC. This study aimed to construct a new hypoxia-related prognosis signature and investigate its potential ceRNA axes in HCC. RNA profiles and hypoxia genes were downloaded, respectively, from the Cancer Genome Atlas hepatocellular carcinoma database and Gene Set Enrichment Analysis website. Cox regression analyses were performed to select the prognostic genes and construct the risk model. The ENCORI database was applied to build the lncRNA-miRNA-mRNA prognosis-related network. The TIMER and CellMiner databases were employed to analyze the association of gene expression in ceRNA with immune infiltration and drug sensitivity, respectively. Finally, the co-expression analysis was carried out to construct the potential lncRNA/miRNA/mRNA regulatory axes. We obtained a prognostic signature including eight hypoxia genes (ENO2, KDELR3, PFKP, SLC2A1, PGF, PPFIA4, SAP30, and TKTL1) and further established a hypoxia-related prognostic ceRNA network including 17 lncRNAs, six miRNAs, and seven mRNAs for hepatocellular carcinoma. Then, the analysis of immune infiltration and drug sensitivity showed that gene expression in the ceRNA network was significantly correlated with the infiltration abundance of multiple immune cells, the expression level of immune checkpoints, and drug sensitivity. Finally, we identified three ceRNA regulatory axes (SNHG1/miR-101-3p/PPFIA4, SNHG1/miR-101-3p/SAP30, and SNHG1/miR-101-3p/TKTL1) associated with the progression of HCC under hypoxia. Here, we constructed a prognosis gene signature and a ceRNA network related to hypoxia for hepatocellular carcinoma. Among the ceRNA network, six highly expressed lncRNAs (AC005540.1, AC012146.1, AC073529.1, AC090772.3, AC138150.2, AL390728.6) and one highly expressed mRNA (PPFIA4) were the potential biomarkers of hepatocellular carcinoma which we firstly reported. The three predicted hypoxia-related regulatory axes may play a vital role in the progression of hepatocellular carcinoma.

Comprehensive analysis of ESR1-related ceRNA axis as a novel prognostic biomarker in hepatocellular carcinoma

Aims: To further understand, detect and treat hepatocellular carcinoma (HCC), it is urgent to conduct more in-depth research on the mechanism of sex-associated differences. Materials & methods: We established a ceRNA triple regulatory axis associated with ESR1 in HCC and performed expression, survival and nuclear-cytoplasmic localization analyses. In addition to this, we performed methylation analysis and immune infiltration analysis of the ceRNA axis. Results: We constructed the LINC01018/hsa-miR-197-3p/GNA14 (lncRNA/miRNA/mRNA) ceRNA axis to further explain the mechanism of sex-related prognosis in the development of HCC and to provide new insights into candidate biomarkers for targeted therapies. Conclusion: Our study is an innovative attempt at demonstrating the mechanism underlying the prognosis associated with sex differences in HCC by constructing a ceRNA axis (LINC01018/hsa-miR-197-3p/GNA14).

Long noncoding RNA SNHG5 promotes podocyte injury via the microRNA-26a-5p/TRPC6 pathway in diabetic nephropathy

Podocyte injury is a characteristic pathological hallmark of diabetic nephropathy (DN). However, the exact mechanism of podocyte injury in DN is incompletely understood. This study was conducted using db/db mice and immortalized mouse podocytes. High-throughput sequencing was used to identify the differentially expressed long noncoding RNAs in kidney of db/db mice. The lentiviral shRNA directed against long noncoding RNA small nucleolar RNA host gene 5 (SNHG5) or microRNA-26a-5p (miR-26a-5p) agomir was used to treat db/db mice to regulate the SNHG5/miR-26a-5p pathway. Here, we found that the expression of transient receptor potential canonical type 6 (TRPC6) was significantly increased in injured podocytes under the condition of DN, which was associated with markedly decreased miR-26a-5p. We determined that miR-26a-5p overexpression ameliorated podocyte injury in DN via binding to 3'-UTR of Trpc6, as evidenced by the markedly reduced activity of luciferase reporters by miR-26a-5p mimic. Then, the upregulated SNHG5 in podocytes and kidney in DN was identified, and it was proved to sponge to miR-26a-5p directly using luciferase activity, RNA immunoprecipitation, and RNA pull-down assay. Knockdown of SNHG5 attenuated podocyte injury in vitro, accompanied by an increased expression of miR-26a-5p and decreased expression of TRPC6, demonstrating that SNHG5 promoted podocyte injury by controlling the miR-26a-5p/TRPC6 pathway. Moreover, knockdown of SNHG5 protects against podocyte injury and progression of DN in vivo. In conclusion, SNHG5 promotes podocyte injury via the miR-26a-5p/TRPC6 pathway in DN. Our findings provide novel insights into the pathophysiology of podocyte injury and a potential new therapeutic strategy for DN.

A Novel Strategy for Identifying NSCLC MicroRNA Biomarkers and Their Mechanism Analysis Based on a Brand-New CeRNA-Hub-FFL Network

Finding reliable miRNA markers and revealing their potential mechanisms will play an important role in the diagnosis and treatment of NSCLC. Most existing computational methods for identifying miRNA biomarkers only consider the expression variation of miRNAs or rely heavily on training sets. These deficiencies lead to high false-positive rates. The independent regulatory model is an important complement to traditional models of co-regulation and is more impervious to the dataset. In addition, previous studies of miRNA mechanisms in the development of non-small cell lung cancer (NSCLC) have mostly focused on the post-transcriptional level and did not distinguish between NSCLC subtypes. For the above problems, we improved mainly in two areas: miRNA identification based on both the NOG network and biological functions of miRNA target genes; and the construction of a 4-node directed competitive regulatory network to illustrate the mechanisms. NSCLC was classified as lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) in this work. One miRNA biomarker of LUAD (miR-708-5p) and four of LUSC (miR-183-5p, miR-140-5p, miR-766-5p, and miR-766-3p) were obtained. They were validated using literature and external datasets. The ceRNA-hub-FFL involving transcription factors (TFs), microRNAs (miRNAs), mRNAs, and long non-coding RNAs (lncRNAs) was constructed. There were multiple interactions among these components within the net at the transcriptional, post-transcriptional, and protein levels. New regulations were revealed by the network. Meanwhile, the network revealed the reasons for the previous conflicting conclusions on the roles of CD44, ACTB, and ITGB1 in NSCLC, and demonstrated the necessity of typing studies on NSCLC. The novel miRNA markers screening method and the 4-node directed competitive ceRNA-hub-FFL network constructed in this work can provide new ideas for screening tumor markers and understanding tumor development mechanisms in depth.

Dysregulation of Pseudogenes/lncRNA-Hsa-miR-1-3p-PAICS Pathway Promotes the Development of NSCLC

Objective

Non-small cell lung cancer (NSCLC) explains about 80 percent of whole lung cancers, and its 5-year survival rate is impoverished, as when people are first diagnosed, 68% of whom are identified at a dangerous stage. The molecular mechanisms of NSCLC are still being explored.

Methods

GSE18842 and GSE19804 were exerted to scan for diversely expressed genes (DEGs) in NSCLC, and then we used GEPIA for the validation of DEGs expression. The prognostic values were determined through Kaplan–Meier analysis. Three target prediction databases indicated potential microRNAs (miRNAs), while miRNet predicted hsa-miR-1-3p′s upstream long non-coding RNAs (lncRNAs) and pseudogenes. UALCAN was utilized to identify the co-expressed genes of PAICS, while enrichment analysis on them was managed with Enrichr.

Results

We initially found that the gene expression level of cyclin B1 (CCNB1), cyclin-dependent kinases1 (CDK1), and phosphoribosylaminoimidazole succinocarboxamide synthetase (PAICS) had a notable increase in NSCLC. We predicted 6, 10, and 7 microRNAs to target CCNB1, CDK1, and PAICS, respectively. Among miRNA-mRNA (microRNA-messenger RNA) pairs, we deduced that the hsa-miR-1-PAICS axis was the most potential one to inhibit the occurrence of NSCLC. We also noted that the hsa-miR-1-3p-PAICS axis participated in regulating the process of mitosis with mechanical functions. Moreover, we identified 5 pseudogenes and 33 long non-coding RNAs (lncRNAs) that might inhibit the hsa-miR-1-3p-PAICS axis in NSCLC.

Conclusions

The pseudogene/lncRNA-hsa-miR-1-3p-PAICS is very important in NSCLC on the basis of this study, thus providing us with effective treatments and promising biomarkers for the diagnosis of NSCLC.

Comprehensive analysis of GINS subunits prognostic value and ceRNA network in sarcoma

Background: The GINS complex, composed of GINS1/2/3/4 subunits, is an essential structure of Cdc45-MCM-GINS (CMG) helicase and plays a vital role in establishing the DNA replication fork and chromosome replication. Meanwhile, GINS genes have been associated with the poor prognosis of various malignancies. However, the abnormal expression of GINS genes and their diagnostic and prognostic value in sarcomas (SARC) remain unclear.

Methods: Oncomine, Gene Expression Profiling Interactive Analysis (GEPIA), Kaplan-Meier Plotter, Cancer cell line encyclopedia (CCLE), The University of Alabama at Birmingham Cancer Data Analysis Portal (UALCAN), R studio, and Tumor Immune Estimation Resource (TIMER) were used to analyze the expression profiles, prognostic value, biological function, ceRNA, and immune infiltration associated with GINS genes in sarcomas.

Results: We found that GINS1/2/3/4 genes exhibited significantly upregulated transcription levels in SARC samples compared to non-tumor tissues and exhibited high expression levels in sarcoma cell lines. In addition, SARC patients with increased expression levels of GINS1/2/3/4 showed poorer survival rates. Immune infiltration analysis showed that GINS subunits were closely associated with the infiltration of immune cells in sarcomas.

Conclusion: Our research identified GINS subunits as potential diagnostic and prognostic biological targets in SARC and elucidated their underlying effects in the genesis and progression of SARC. These results may provide new opportunities and research directions for targeted sarcoma therapy.

Pan-cancer identification of the relationship of metabolism-related differentially expressed transcription regulation with non-differentially expressed target genes via a gated recurrent unit network

The transcriptome describes the expression of all genes in a sample. Most studies have investigated the differential patterns or discrimination powers of transcript expression levels. In this study, we hypothesized that the quantitative correlations between the expression levels of transcription factors (TFs) and their regulated target genes (mRNAs) serve as a novel view of healthy status, and a disease sample exhibits a differential landscape (mqTrans) of transcription regulations compared with healthy status. We formulated quantitative transcription regulation relationships of metabolism-related genes as a multi-input multi-output regression model via a gated recurrent unit (GRU) network. The GRU model was trained using healthy blood transcriptomes and the expression levels of mRNAs were predicted by those of the TFs. The mqTrans feature of a gene was defined as the difference between its predicted and actual expression levels. A pan-cancer investigation of the differentially expressed mqTrans features was conducted between the early- and late-stage cancers in 26 cancer types of The Cancer Genome Atlas database. This study focused on the differentially expressed mqTrans features, that did not show differential expression in the actual expression levels. These genes could not be detected by conventional differential analysis. Such dark biomarkers are worthy of further wet-lab investigation. The experimental data also showed that the proposed mqTrans investigation improved the classification between early- and late-stage samples for some cancer types. Thus, the mqTrans features serve as a complementary view to transcriptomes, an OMIC type with mature high-throughput production technologies, and abundant public resources.

Prognostic Roles of ceRNA Network-Based Signatures in Gastrointestinal Cancers

Gastrointestinal cancers (GICs) are high-incidence malignant tumors that seriously threaten human health around the world. Their complexity and heterogeneity make the classic staging system insufficient to guide patient management. Recently, competing endogenous RNA (ceRNA) interactions that closely link the function of protein-coding RNAs with that of non-coding RNAs, such as long non-coding RNA (lncRNA) and circular RNA (circRNA), has emerged as a novel molecular mechanism influencing miRNA-mediated gene regulation. Especially, ceRNA networks have proven to be powerful tools for deciphering cancer mechanisms and predicting therapeutic responses at the system level. Moreover, abnormal gene expression is one of the critical breaking events that disturb the stability of ceRNA network, highlighting the role of molecular biomarkers in optimizing cancer management and treatment. Therefore, developing prognostic signatures based on cancer-specific ceRNA network is of great significance for predicting clinical outcome or chemotherapy benefits of GIC patients. We herein introduce the current frontiers of ceRNA crosstalk in relation to their pathological implications and translational potentials in GICs, review the current researches on the prognostic signatures based on lncRNA or circRNA-mediated ceRNA networks in GICs, and highlight the translational implications of ceRNA signatures for GICs management. Furthermore, we summarize the computational approaches for establishing ceRNA network-based prognostic signatures, providing important clues for deciphering GIC biomarkers.

The role of N6-methyladenosine-modified non-coding RNAs in the pathological process of human cancer

Non-coding RNAs (ncRNAs) account for the majority of the widespread transcripts of mammalian genomes. They rarely encode proteins and peptides, but their regulatory role is crucial in numerous physiological and pathological processes. The m6A (N6-methyladenosine) modification is one of the most common internal RNA modifications in eukaryotes and is associated with all aspects of RNA metabolism. Accumulating researches have indicated a close association between m6A modification and ncRNAs, and suggested m6A-modified ncRNAs played a crucial role in tumor progression. The correlation between m6A modification and ncRNAs offers a novel perspective for investigating the potential mechanisms of cancer pathological processes, which suggests that both m6A modification and ncRNAs are critical prognostic markers and therapeutic targets in numerous malignancies. In the present report, we summarized the interaction between m6A modification and ncRNA, emphasizing how their interaction regulates pathological processes in cancer.

CircRNA-Associated CeRNAs Regulatory Axes in Retinoblastoma: A Systematic Scoping Review

Retinoblastoma (RB) is one of the most common childhood cancers caused by RB gene mutations (tumor suppressor gene in various patients). A better understanding of molecular pathways and the development of new diagnostic approaches may lead to better treatment for RB patients. The number of studies on ceRNA axes is increasing, emphasizing the significance of these axes in RB. Circular RNAs (circRNAs) play a vital role in competing endogenous RNA (ceRNA) regulatory axes by sponging microRNAs and regulating gene expression. Because of the broadness of ceRNA interaction networks, they may assist in investigating treatment targets in RB. This study conducted a systematic scoping review to evaluate verified loops of ceRNA in RB, focusing on the ceRNA axis and its relationship to circRNAs. This scoping review was carried out using a six-step strategy and the Prisma guideline, and it involved systematically searching the publications of seven databases. Out of 363 records, sixteen articles were entirely consistent with the defined inclusion criteria and were summarized in the relevant table. The majority of the studies focused on the circRNAs circ_0000527, circ_0000034, and circTET1, with approximately two-fifths of the studies focusing on a single circRNA. Understanding the many features of this regulatory structure may help elucidate RB’s unknown causative factors and provide novel molecular potential therapeutic targets and medical fields.

Silencing of Circ_0135889 Restrains Proliferation and Tumorigenicity of Human Neuroblastoma Cells

INTRODUCTION

Neuroblastoma (NB) is the most common extracranial solid tumor in infants and young children. Circular ribonucleic acid (RNA) hsa_circ_0135889 (circ_0135889; hsa_circ argonaute 2 _001) is highly expressed in multiple cancer tissues, including NB. However, its role in tumor progression of NB was unclear.

METHODS

Real-time quantitative polymerace chain reaction was used to detect RNA expression, and western blotting, or immunohistochemistry was used to measure protein expression. Functional experiments were performed using 3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium, 5-Ethynyl-2'- deoxyuridine, Annexin V-fluorescein isothiocyanate/propidium iodide, and transwell assays, as well as xenograft tumor model. The intermolecular interaction was predicted by online databases and confirmed by dual-luciferase reporter assay and RNA pull-down assay.

RESULTS

Circ_0135889 and neuronal differentiation 1 (NEUROD1) were upregulated whilst microRNA (miR)-127-5p, was downregulated in NB tumors and immortalized NB cells. Silencing of circ_0135889 could suppress cell proliferation, migration and invasion, but enhance apoptosis rate of NB cells in vitro. More importantly, circ_0135889 depletion inhibited xenograft tumor growth of NB cells. Circ_0135889 was a sponge for miR-127-5p, and inhibition of miR-127-5p counteracted the inhibitory impact of circ_0135889 knockdown on the malignant behaviors of NB cells. Moreover, NEUROD1 was a direct target of miR-127-5p, and miR-127-5p exerted the anti-tumor role in NB cells by targeting NEUROD1. Furthermore, circ_0135889 regulated NEUROD1expression by sponging miR-127-5p.

CONCLUSIONS

Circ_0135889 promoted the tumorigenicity of NB by regulating miR-127-5p/NEUROD1 axis, which might provide a promising therapeutic target for NB.

A Proposed TUSC7/miR-211/Nurr1 ceRNET Might Potentially be Disturbed by a cer-SNP rs2615499 in Breast Cancer

Evidence and in silico analyses showed that TUSC7, miR-211, and Nurr1 may be involved in BC pathogenesis by ceRNET signaling axis. This study aimed to investigate the potential role of TUSC7/miR-211/Nurr1 ceRNET and rs2615499 variant as a novel cer-SNP in BC subjects. The expression assays were conducted by qPCR on tumor tissues (n = 50), tumor-adjacent normal tissues (TANTs) (n = 50), and clinically healthy control tissues (n = 50). The expression of TUSC7 and Nurr1 significantly decreased, but the level of miR-211 significantly increased in tumor tissues compared to TANTs and healthy normal tissues. Altered expression of TUSC7 and miR-211 was associated with poor prognosis of patients. The Nurr1 exhibited a double-edged sword-like activity in BC. In addition, TUSC7, Nurr1, and miR-211 expressions were significantly related to a novel BC-associated rs2615499 (A > C) located in the miR-211 binding site on Nurr1 3'-UTR. In the second part of the study, a case-control study was performed on BC patients (n = 100) and matched healthy controls (n = 100). The genomic DNA was isolated and genotyping was performed using Tetra-Primer ARMS PCR. The CC and AC genotypes were associated with higher expression levels of Nurr1 and worse outcomes of the disease. Our findings revealed that TUSC7 functions as a tumor suppressor in BC potentially via miR-211/Nurr1, which might be disturbed by the cer-SNP rs2615499. However, functional studies are needed to validate these results.

Multivariate competing endogenous RNA network characterization for cancer MicroRNA biomarker discovery: a novel bioinformatics model with application to prostate cancer metastasis

Background

MicroRNAs (miRNAs) are post-transcriptional regulators with potential as biomarkers for cancer management. Data-driven competing endogenous RNA (ceRNA) network modeling is an effective way to decipher the complex interplay between miRNAs and spongers. However, there are currently no general rules for ceRNA network-based biomarker prioritization.

Methods and results

In this study, a novel bioinformatics model was developed by integrating gene expression with multivariate miRNA-target data for ceRNA network-based biomarker discovery. Compared with traditional methods, the structural vulnerability in the human long non-coding RNA (lncRNA)–miRNA–messenger RNAs (mRNA) network was comprehensively analyzed, and the single-line regulatory or competing mode among miRNAs, lncRNAs, and mRNAs was characterized and quantified as statistical evidence for miRNA biomarker identification. The application of this model to prostate cancer (PCa) metastasis identified a total of 12 miRNAs as putative biomarkers from the metastatic PCa-specific lncRNA–miRNA–mRNA network and nine of them have been previously reported as biomarkers for PCa metastasis. The receiver operating characteristic curve and cell line qRT-PCR experiments demonstrated the power of miR-26b-5p, miR-130a-3p, and miR-363-3p as novel candidates for predicting PCa metastasis. Moreover, PCa-associated pathways such as prostate cancer signaling, ERK/MAPK signaling, and TGF-β signaling were significantly enriched by targets of identified miRNAs, indicating the underlying mechanisms of miRNAs in PCa carcinogenesis.

Conclusions

A novel ceRNA-based bioinformatics model was proposed and applied to screen candidate miRNA biomarkers for PCa metastasis. Functional validations using human samples and clinical data will be performed for future translational studies on the identified miRNAs.

Identification of Differentially Expressed and Prognostic lncRNAs for the Construction of ceRNA Networks in Lung Adenocarcinoma

Background

Long noncoding RNAs (lncRNAs) could function as competitive endogenous RNAs (ceRNAs) to competitively adsorb microRNAs (miRNAs), thereby regulating the expression of their target protein-coding mRNAs. In this study, we aim to identify more effective diagnostic and prognostic markers for lung adenocarcinoma (LUAD).

Methods

We obtained differentially expressed lncRNAs (DElncRNAs), miRNAs (DEmiRNAs), and mRNAs (DEmRNAs) for LUAD by using The Cancer Genomes Atlas (TCGA) portal. Weighted gene coexpression network analysis (WGCNA) was performed to unveil core gene modules associated with LUAD. The Cox proportional hazards model was performed to determine the prognostic significance of DElncRNAs. The diagnostic and prognostic significance of DElncRNAs was further verified based on the receiver operating characteristic curve (ROC). Cytoscape was used to construct the ceRNA networks comprising the lncRNAs-miRNAs-mRNAs axis based on the correlation obtained from the miRcode, miRDB, and TargetScan.

Results

Compared with normal lung tissues, 2355 DElncRNAs, 820 DEmiRNAs, and 17289 DEmRNAs were identified in LUAD tissues. We generated 8 WGCNA core modules in the lncRNAs coexpression network, 5 modules in the miRNAs, and 12 modules in the mRNAs coexpression network, respectively. One lncRNA module (blue) consisting of 441 lncRNAs, two miRNA modules (blue and turquoise) containing 563 miRNAs, and one mRNA module (turquoise), which consisted of 15162 mRNAs, were mostly significantly related to LUAD status. Furthermore, 67 DEmRNAs were found to be tumor-associated as well as the target genes of the DElncRNAs-DEmiRNAs axis. Survival analyses showed that 6 lncRNAs (LINC01447, WWC2-AS2, OGFRP1, LINC00942, LINC01168, and AC005863.1) were significantly correlated with the prognosis of LUAD patients. Ultimately, the potential ceRNA networks including 6 DElncRNAs, 4 DEmiRNAs, and 22 DEmRNAs were constructed.

Conclusion

Our study indicated that 6 DElncRNAs had the possibilities as diagnostic and prognostic biomarkers for LUAD. The lncRNA-mediated ceRNA networks might provide novel insights into the molecular mechanisms of LUAD progression.

The role of non-coding RNA in the diagnosis and treatment of Helicobacter pylori-related gastric cancer, with a focus on inflammation and immune response

Helicobacter pylori (H. pylori) is one of the globally recognized causative factors of gastric cancer (GC). Currently, no definite therapy and drugs for H. pylori-related GC have been widely acknowledged although H. pylori infection could be eradicated in early stage. Inflammation and immune response are spontaneous essential stages during H. pylori infection. H pylori may mediate immune escape by affecting inflammation and immune response, leading to gastric carcinogenesis. As an important component of transcriptome, non-coding RNAs (ncRNAs) have been proven to play crucial roles in the genesis and development of H. pylori-induced GC. This review briefly described the effects of ncRNAs on H. pylori-related GC from the perspective of inflammation and immune response, as well as their association with inflammatory reaction and immune microenvironment. We aim to explore the potential of ncRNAs as markers for the early diagnosis, prognosis, and treatment of H. pylori-related GC. The ncRNAs involved in H. pylori-related GC may all hold promise as novel therapeutic targets for immunotherapy.

HFBD: a biomarker knowledge database for heart failure heterogeneity and personalized applications

MOTIVATION

Heart failure (HF) is a cardiovascular disease with a high incidence around the world. Accumulating studies have focused on the identification of biomarkers for HF precision medicine. To understand the HF heterogeneity and provide biomarker information for the personalized diagnosis and treatment of HF, a knowledge database collecting the distributed and multiple-level biomarker information is necessary.

RESULTS

In this study, the HF biomarker knowledge database (HFBD) was established by manually collecting the data and knowledge from literature in PubMed. HFBD contains 2618 records and 868 HF biomarkers (731 single and 137 combined) extracted from 1237 original articles. The biomarkers were classified into proteins, RNAs, DNAs, and the others at molecular, image, cellular and physiological levels. The biomarkers were annotated with biological, clinical and article information as well as the experimental methods used for the biomarker discovery. With its user-friendly interface, this knowledge database provides a unique resource for the systematic understanding of HF heterogeneity and personalized diagnosis and treatment of HF in the era of precision medicine.

AVAILABILITY

The platform is openly available at http://sysbio.org.cn/HFBD/.

Competing Endogenous RNAs in Cervical Carcinogenesis: A New Layer of Complexity

MicroRNAs (miRNAs) regulate gene expression by binding to complementary sequences within target mRNAs. Apart from working ‘solo’, miRNAs may interact in important molecular networks such as competing endogenous RNA (ceRNA) axes. By competing for a limited pool of miRNAs, transcripts such as long noncoding RNAs (lncRNAs) and mRNAs can regulate each other, fine-tuning gene expression. Several ceRNA networks led by different lncRNAs—described here as lncRNA-mediated ceRNAs—seem to play essential roles in cervical cancer (CC). By conducting an extensive search, we summarized networks involved in CC, highlighting the major impacts of such dynamic molecular changes over multiple cellular processes. Through the sponging of distinct miRNAs, some lncRNAs as HOTAIR, MALAT1, NEAT1, OIP5-AS1, and XIST trigger crucial molecular changes, ultimately increasing cell proliferation, migration, invasion, and inhibiting apoptosis. Likewise, several lncRNAs seem to be a sponge for important tumor-suppressive miRNAs (as miR-140-5p, miR-143-3p, miR-148a-3p, and miR-206), impairing such molecules from exerting a negative post-transcriptional regulation over target mRNAs. Curiously, some of the involved mRNAs code for important proteins such as PTEN, ROCK1, and MAPK1, known to modulate cell growth, proliferation, apoptosis, and adhesion in CC. Overall, we highlight important lncRNA-mediated functional interactions occurring in cervical cells and their closely related impact on cervical carcinogenesis.

miRSM: an R package to infer and analyse miRNA sponge modules in heterogeneous data

In molecular biology, microRNA (miRNA) sponges are RNA transcripts which compete with other RNA transcripts for binding with miRNAs. Research has shown that miRNA sponges have a fundamental impact on tissue development and disease progression. Generally, to achieve a specific biological function, miRNA sponges tend to form modules or communities in a biological system. Until now, however, there is still a lack of tools to aid researchers to infer and analyse miRNA sponge modules from heterogeneous data. To fill this gap, we develop an R/Bioconductor package, miRSM, for facilitating the procedure of inferring and analysing miRNA sponge modules. miRSM provides a collection of 50 co-expression analysis methods to identify gene co-expression modules (which are candidate miRNA sponge modules), four module discovery methods to infer miRNA sponge modules and seven modular analysis methods for investigating miRNA sponge modules. miRSM will enable researchers to quickly apply new datasets to infer and analyse miRNA sponge modules, and will consequently accelerate the research on miRNA sponges.

Circular RNA circPIKfyve acts as a sponge of miR-21-3p to enhance antiviral immunity through regulating MAVS in teleost fish

Circular RNAs (circRNAs) are a class of widespread and diverse covalently closed circular endogenous RNAs that exert crucial functions in regulating gene expression in mammals. However, the function and regulation mechanism of circRNAs in lower vertebrates are still unknown. Here, we discovered a novel circRNA derived from PIKfyve, named circPIKfyve, that is related to the antiviral responses in teleost fish. The results showed that circPIKfyve plays essential roles in host antiviral immunity and inhibition of SCRV replication. Moreover, we also found that the antiviral effect inhibited by miR-21-3p could be reversed with the addition of circPIKfyve. In mechanism, our data revealed that circPIKfyve is a competitive endogenous RNA (ceRNA) of MAVS by sponging miR-21-3p, leading to activation of NF-κB/IRF3 pathway, which then enhance the innate antiviral responses. In addition, we firstly found that RNA binding protein QKI is involved in the formation and regulation of circPIKfyve. Our results provided a strong basis that circRNAs to play a regulatory role in antiviral immune responses in teleost fish.Importance: Here, we identified a novel circRNA, namely, circPIKfyve, that can act as a key regulator of the innate immune response in teleost fish. circPIKfyve acts as a molecular sponge by competitive adsorbing of miR-21-3p, thereby increasing the abundance of MAVS and activating the downstream NF-κB/IRF3 pathway to enhance the antiviral response. In addition, this study was the first to find that QKI protein is involved in regulating the formation of circPIKfyve in fish. The overall results of this study suggest that circPIKfyve plays an active regulatory role in the antiviral immune response of teleost fish.

CeNet Omnibus: an R/Shiny application to the construction and analysis of competing endogenous RNA network

Background

The competing endogenous RNA (ceRNA) regulation is a newly discovered post-transcriptional regulation mechanism and plays significant roles in physiological and pathological progress. CeRNA networks provide global views to help understand the regulation of ceRNAs. CeRNA networks have been widely used to detect survival biomarkers, select candidate regulators of disease genes, and predict long noncoding RNA functions. However, there is no software platform to provide overall functions from the construction to analysis of ceRNA networks.

Results

To fill this gap, we introduce CeNet Omnibus, an R/Shiny application, which provides a unified framework for the construction and analysis of ceRNA network. CeNet Omnibus enables users to select multiple measurements, such as Pearson correlation coefficient (PCC), mutual information (MI), and liquid association (LA), to identify ceRNA pairs and construct ceRNA networks. Furthermore, CeNet Omnibus provides a one-stop solution to analyze the topological properties of ceRNA networks, detect modules, and perform gene enrichment analysis and survival analysis. CeNet Omnibus intends to cover comprehensiveness, high efficiency, high expandability, and user customizability, and it also offers a web-based user-friendly interface to users to obtain the output intuitionally.

Conclusion

CeNet Omnibus is a comprehensive platform for the construction and analysis of ceRNA networks. It is highly customizable and outputs the results in intuitive and interactive. We expect that CeNet Omnibus will assist researchers to understand the property of ceRNA networks and associated biological phenomena. CeNet Omnibus is an R/Shiny application based on the Shiny framework developed by RStudio. The R package and detailed tutorial are available on our GitHub page with the URL https://github.com/GaoLabXDU/CeNetOmnibus.

Construction and Analysis of Survival-Associated Competing Endogenous RNA Network in Lung Adenocarcinoma

Increasing evidence has shown that noncoding RNAs play significant roles in the initiation, progression, and metastasis of tumours via participating in competing endogenous RNA (ceRNA) networks. However, the survival-associated ceRNA in lung adenocarcinoma (LUAD) remains poorly understood. In this study, we aimed to investigate the regulatory mechanisms underlying ceRNA in LUAD to identify novel prognostic factors. mRNA, lncRNA, and miRNA sequencing data obtained from the GDC data portal were utilized to identify differentially expressed (DE) RNAs. Survival-related RNAs were recognized using univariate Kaplan-Meier survival analysis. We performed functional enrichment analysis of survival-related mRNAs using the clusterProfiler package of R and STRING. lncRNA-miRNA and miRNA-mRNA interactions were predicted based on miRcode, Starbase, and miRanda. Subsequently, the survival-associated ceRNA network was constructed for LUAD. Multivariate Cox regression analysis was used to identify prognostic factors. Finally, we acquired 15 DE miRNAs, 49 DE lncRNAs, and 843 DE mRNAs associated with significant overall survival. Functional enrichment analysis indicated that survival-related DE mRNAs were enriched in cell cycle. The survival-associated lncRNA-miRNA-mRNA ceRNA network was constructed using five miRNAs, 49 mRNAs, and 21 lncRNAs. Furthermore, seven hub RNAs (LINC01936, miR-20a-5p, miR-31-5p, TNS1, TGFBR2, SMAD7, and NEDD4L) were identified based on the ceRNA network. LINC01936 and miR-31-5p were found to be significant using the multifactorial Cox regression model. In conclusion, we successfully constructed a survival-related lncRNA-miRNA-mRNA ceRNA regulatory network in LUAD and identified seven hub RNAs, which provide novel insights into the regulatory molecular mechanisms associated with survival of LUAD, and identified two independent prognostic predictors for LUAD.

Long Non-Coding RNA ZNF667-AS1 Knockdown Curbs Liver Metastasis in Acute Myeloid Leukemia by Regulating the microRNA-206/AKAP13 Axis

BACKGROUND

Zinc finger protein 667-antisense RNA 1 (ZNF667-AS1), a long non-coding RNA (lncRNA), plays important parts in tumorigenesis and development of esophageal squamous cell carcinoma, but its function in acute myeloid leukemia (AML) is unknown. Our goal here was to probe the functional mechanism of ZNF667-AS1 in AML by mediating microRNA-206 (miR-206)/A-kinase anchoring protein 13 (AKAP13) axis.

MATERIALS AND METHODS

The bone marrow samples from AML patients and controls were selected for microarray analysis to select significantly upregulated lncRNAs. Next, effects of ZNF667-AS1 on cell aggressiveness of AML were assessed after delivery of cells with siRNA against ZNF667-AS1. Subcellular fractionation location assay and FISH experiments were used to determine ZNF667-AS1 localization in cells. Dual-luciferase experiments detect the targeting relationships among ZNF667-AS1, miR-206 and AKAP13. Finally, tumor growth and metastasis were evaluated in vivo to determine the relevance of ZNF667-AS1/miR-206/AKAP13 axis.

RESULTS

The expression of ZNF667-AS1 was upregulated in AML patients, which predicted poor prognosis. Downregulation of ZNF667-AS1 reduced cell proliferation, invasion, tumorigenesis and metastasis. miR-206 inhibitor reversed the repressive role of ZNF667-AS1 knockdown in cell proliferation, invasion and tumorigenesis, while AKAP13 silencing flattened the stimulative role of miR-206 inhibitor in AML malignant aggressiveness. Mechanistically, we demonstrated that ZNF667-AS1 functioned as a molecular sponge for miR-206. In addition, we observed that Wnt/β-catenin pathway was suppressed by ZNF667-AS1 knockdown.

CONCLUSION

ZNF667-AS1 potentiated AML progression by targeting the miR-206/AKAP13 axis. This indicates ZNF667-AS 1 inhibition may act as a prospective therapeutic option for the treatment of AML.

Biomarker Discovery for the Carcinogenic Heterogeneity Between Colon and Rectal Cancers Based on lncRNA-Associated ceRNA Network Analysis

Background

Colorectal cancer (CRC) is one of the leading causes of cancer death worldwide. Emerging evidence has revealed that risk factors and metastatic patterns differ greatly between colon and rectal cancers. However, the molecular mechanism underlying their pathogenic differences remains unclear. Therefore, we here aimed to identify non-coding RNA biomarkers based on lncRNA-associated ceRNA network (LceNET) to elucidate the carcinogenic heterogeneity between colon and rectal cancers.

Methods

A global LceNET in human was constructed by employing experimental evidence-based miRNA-mRNA and miRNA-lncRNA interactions. Then, four context-specific ceRNA networks related to cancer initiation and metastasis were extracted by mapping differentially expressed lncRNAs, miRNAs and mRNAs to the global LceNET. Notably, a novel network-based bioinformatics model was proposed and applied to identify lncRNA/miRNA biomarkers and critical ceRNA triplets for understanding the carcinogenic heterogeneity between colon and rectal cancers. Moreover, the identified biomarkers were further validated by their diagnostic/prognostic performance, expression pattern and correlation analysis.

Results

Based on network modeling, lncRNA KCNQ1OT1 (AUC>0.85) and SNHG1 (AUC>0.94) were unveiled as common diagnostic biomarkers for the initiation and metastasis of colon and rectal cancers. qRT-PCR analysis uncovered that these lncRNAs had significantly higher expression level in CRC cell lines with high metastatic potential. In particular, KCNQ1OT1 and SNHG1 function in colon and rectal cancers via different ceRNA mechanisms. For example, KCNQ1OT1/miR-484/ANKRD36 axis was involved in the initiation of colon cancer, while KCNQ1OT1/miR-181a-5p/PCGF2 axis was implicated in the metastasis of rectal cancer; the SNHG1/miR-484/ORC6 axis played a role in colon cancer, while SNHG1/miR-423-5p/EZH2 and SNHG1/let-7b-5p/ATP6V1F axes participated in the initiation and metastasis of rectal cancer, respectively. In these ceRNA triplets, miR-484, miR-181a-5p, miR-423-5p and let-7b-5p were identified as miRNA biomarkers with excellent distinguishing ability between normal and tumor tissues, and ANKRD36, PCGF2, EZH2 and ATP6V1F were closely related to the prognosis of corresponding cancer.

Conclusion

The landscape of lncRNA-associated ceRNA network not only facilitates the exploration of non-coding RNA biomarkers, but also provides deep insights into the oncogenetic heterogeneity between colon and rectal cancers, thereby contributing to the optimization of diagnostic and therapeutic strategies of CRC.

Inferences of individual differences in response to tripterysium glycosides across patients with Rheumatoid arthritis using a novel ceRNA regulatory axis

BACKGROUND

To identify biomarkers for guiding therapy and predicting clinical response of Tripterysium Glycosides Tablets (TGT) treatment is an urgent task due to individual differences in TGT response across rheumatoid arthritis (RA) patients. Competing endogenous RNA (ceRNA) regulatory system may influence drug response with involvement in diverse biological processes. Herein, we aimed to identify a TGT response-related ceRNA axis.

METHODS

A TGT response-related ceRNA axis was screened according to clinical cohort-based RNA expression profiling, lncRNA-mRNA coexpression, and ceRNA network analyses. Its clinical relevance was evaluated by computational modeling. Regulatory mechanisms of ceRNA axis were also experimentally investigated.

RESULTS

The ceRNA regulatory axis combined with lncRNA ENST00000494760, miR-654-5p, and C1QC was identified as a candidate biomarker for RA patients' response to TGT. Both ENST00000494760 and C1QC mRNA expression were significantly lower, while miR-654-5p expression was dramatically higher in TGT responders than nonresponders. Its clinical relevance was verified by computational modeling based on both independent clinical validation cohort and collagen-induced arthritis (CIA) mice. Mechanistically, miR-654-5p directly bound to the 3'-untranslated region of both ENST00000494760 and C1QC mRNA to inhibit their expression. Moreover, miR-654-5p suppressed C1QC mRNA expression, but ENST00000494760 bound to miR-654-5p and relieved its repression on C1QC mRNA, leading to RA aggressive progression and weak TGT response.

CONCLUSIONS

LncRNA ENST00000494760 overexpression may sponge miR-654-5p to promote C1QC expression in RA patients. This novel ceRNA axis may serve as a biomarker for screening the responsive RA patients to TGT treatment, which will allow improved personalized healthcare.

Data-driven translational prostate cancer research: from biomarker discovery to clinical decision

Prostate cancer (PCa) is a common malignant tumor with increasing incidence and high heterogeneity among males worldwide. In the era of big data and artificial intelligence, the paradigm of biomarker discovery is shifting from traditional experimental and small data-based identification toward big data-driven and systems-level screening. Complex interactions between genetic factors and environmental effects provide opportunities for systems modeling of PCa genesis and evolution. We hereby review the current research frontiers in informatics for PCa clinical translation. First, the heterogeneity and complexity in PCa development and clinical theranostics are introduced to raise the concern for PCa systems biology studies. Then biomarkers and risk factors ranging from molecular alternations to clinical phenotype and lifestyle changes are explicated for PCa personalized management. Methodologies and applications for multi-dimensional data integration and computational modeling are discussed. The future perspectives and challenges for PCa systems medicine and holistic healthcare are finally provided.