Construction of competitive endogenous RNA network reveals regulatory role of long non‐coding RNAs in type 2 diabetes mellitus

Increasing epidemic of type 2 diabetes mellitus (T2DM) and its comorbidities makes it urgent to understand the pathogenesis and regulatory mechanism. However, little is known about the regulatory role of lncRNAs in diabetes. Here, we constructed a T2DM‐related competitive endogenous RNA (ceRNA) network (DMCN) to explore biological function of lncRNAs during the development of diabetes mellitus. This network contained 351 nodes including 98 mRNAs, 86 microRNAs and 167 lncRNAs. Functional analysis showed that the mRNAs in DMCN were annotated into some diabetes‐related pathways. Furthermore, mTOR‐centred subnetwork was extracted and ncRNA‐involved mTOR pathway was established. Finally, we validated that NEAT1 was potentially communicated with mTOR signalling target protein mLST8 via the association with miR‐181b. These findings provide significant insight into lncRNA regulatory network in T2DM.     

Genome-wide identification of a competing endogenous RNA network in cholangiocarcinoma

Cholangiocarcinoma (CCA) is the second widespread liver tumor with relatively poor survival. Increasing evidence in recent studies showed long noncoding RNAs (lncRNAs) exert a crucial impact on the development and progression of CCA based on the mechanism of competing endogenous RNAs (ceRNAs). However, functional roles and regulatory mechanisms of lncRNA-regulated ceRNA in CCA, are only partially understood. The expression profile of messenger RNAs (mRNAs), lncRNAs, and microRNAs (miRNAs) downloaded from The Cancer Genome Atlas were comprehensively investigated. Differential expression of these three types of RNA between CCA and corresponding precancerous tissues were screened out for further analysis. On the basis of interactive information generated from miRDB, miRTarBase, TargetScan, and miRcode public databases, we then constructed an mRNA-miRNA-lncRNA regulatory network. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology analyses were conducted to identify the biological function of the ceRNA network involved in CCA. As a result, 2883 mRNAs, 136 miRNAs, and 993 lncRNAs were screened out as differentially expressed RNAs in CCA. In addition, a ceRNA network in CCA was constructed, composing of 50 up and 27 downregulated lncRNAs, 14 up and 7 downregulated miRNAs, 29 up and 25 downregulated mRNAs. Finally, gene set enrichment and pathway analysis indicated our CCA-specific ceRNA network was related with cancer-related pathway and molecular function. In conclusion, our research identified a novel lncRNA-related ceRNA network in CCA, which might act as a potential therapeutic target for patients with CCA.     

Identification of a potential prognostic lncRNA-miRNA-mRNA signature in endometrial cancer based on the competing endogenous RNA network

Endometrial cancer is one of the most common gynecological malignant tumors. The roles of competing endogenous RNAs (ceRNAs) in this disease, however, remain unclear. In this study, we constructed a ceRNA network to reveal the core ceRNAs in endometrial cancer. Differentially expressed genes were summarized from The Cancer Genome Atlas database, whereupon 140 genes were identified for building the network. Further correlation, survival, and enrichment analyses suggested that these genes may help towards elucidating the molecular mechanisms of endometrial cancer. After validation of the findings with the GSE17025 data set, LINC00958, microRNA-761, and DOLPP1 were highlighted as the critical genes in the ceRNA network. Our work suggests that LINC00958 may regulate DOLPP1 by “sponging” miR-761 in endometrial cancer.     

The clinical significance of circulating DSCAM-AS1 in patients with ER-positive breast cancer and construction of its competitive endogenous RNA network

Molecular Biology Reports - Long Non-Coding RNAs (lncRNAs), with diagnostic and therapeutic applications in malignancies, are newly described tumour-related molecules. Here, we reported the...     

Construction of a three-component regulatory network of transcribed ultraconserved regions for the identification of prognostic biomarkers in gastric cancer

Altered expression and functional roles of the transcribed ultraconserved regions (T-UCRs), as genomic sequences with 100% conservation between the genomes of human, mouse, and rat, in the pathophysiology of neoplasms has already been investigated. Nevertheless, the relevance of the functions for T-UCRs in gastric cancer (GC) is still the subject of inquiry. In the current study, we first used a genome-wide profiling approach to analyze the expression of T-UCRs in GC patients. Then, we constructed a three-component regulatory network and investigated potential diagnostic and prognostic values of the T-UCRs. The Cancer Genome Atlas Stomach Adenocarcinoma (TCGA-STAD) dataset was used as a resource for the RNA-sequencing data. FeatureCounts was utilized to quantify the number of reads mapped to each T-UCR. Differential expression analysis was then conducted using DESeq2. In the following, interactions between T-UCRs, microRNAs (miRNAs), and messenger RNAs (mRNAs) were combined into a three-component network. Enrichment analyses were performed and a protein–protein interaction (PPI) network was constructed. The R Survival package was utilized to identify survival-related significantly differentially expressed T-UCRs (DET-UCRs). Using an in-house cohort of GC tissues, expression of two DET-UCRs was furthermore experimentally verified. Our results showed that several T-UCRs were dysregulated in TCGA-STAD tumoral samples compared to nontumoral counterparts. The three-component network was constructed which composed of DET-UCRs, miRNAs, and mRNAs nodes. Functional enrichment and PPI network analyses revealed important enriched signaling pathways and gene ontologies such as “pathway in cancer” and regulation of cell proliferation and apoptosis. Five T-UCRs were significantly correlated with the overall survival of GC patients. While no expression of uc.232 was observed in our in-house cohort of GC tissues, uc.343 showed an increased expression, although not statistically significant, in gastric tumoral tissues. The constructed three-component regulatory network of T-UCRs in GC presents a comprehensive understanding of the underlying gene expression regulation processes involved in tumor development and can serve as a basis to investigate potential prognostic biomarkers and therapeutic targets.     

Construction of a lncRNA/pseudogene-hsa-miR-30d-5p-GJA1 regulatory network related to metastasis of pancreatic cancer

Pancreatic cancer, the most lethal malignant tumor, is notorious for its poor prognosis and metastatic potential. Non-coding RNAs (ncRNAs) are reported to play key roles in cancer metastasis. In this study, miRNA and gene expression profiles between metastatic pancreatic cancer cell M8 and its parental cell BxPC.3 were determined. Using differential expression analysis, survival analysis, target gene prediction, pathway enrichment analysis, intersection analysis and correlation analysis, hsa-miR-30d-5p/GJA1 axis was identified as the most potential pathway involved in metastasis of pancreatic cancer. Subsequently, two upstream lncRNAs (HELLPAR and OIP-AS1) and four upstream pseudogenes (AC093616.1, AC009951.1, TMEM183B and PABPC1P4) of hsa-miR-30d-5p/GJA1 axis were predicted and were then identified via assessment of RNA-RNA expression relationship. Furthermore, CTNNA1, CTNNB1 and CTNND1 were regarded as three crucial molecules to be participated in hsa-miR-30d-5p/GJA1-mediated metastatic potential in pancreatic cancer. In conclusion, we established a novel lncRNA/pseudogene-hsa-miR-30d-5p-GJA1 regulatory network linked to metastasis of pancreatic cancer.     

Comprehensive analysis of microRNA-messenger RNA regulatory network in gemcitabine-resistant bladder cancer cells

Chemotherapy is still a standard treatment of unresectable bladder cancer or distant metastases. The chemotherapy resistance always occurs after a period of treatment indicating poor prognosis. The current study aimed to explore the molecular mechanism of chemoresistance in bladder cancer cells. The gene expression profiles of GSE77883, including three untreated T24 cells samples and three gemcitabine-resistant T24 cells samples, was downloaded from Gene Expression Omnibus database. The screening of differentially expressed genes (DEGs), gene function analysis, and interaction prediction between microRNAs (miRNAs) and DEGs were performed by R software. The protein-protein interaction (PPI) and miRNA-DEGs networks were constructed and visualized by Cytoscape software. Then, the small molecules, with potential synergistic or antagonistic effects to gemcitabine resistance, were identified using the Connectivity Map database. Finally, gemcitabine-resistant T24 cell line was established and key genes were validated by quantitative real-time polymerase chain reaction (qRT-PCR). In total, 536 upregulated and 513 downregulated genes were screened and mainly enriched in oxidative stress response and signaling pathways related to extracellular matrix–receptor interaction and focal adhesion. PPI network showed interleukin 6, tumor necrosis factor, kinesin family member 11, and BUB1 mitotic checkpoint serine/threonine kinase B were key genes. The miRNA-DEGs regulatory networks included 18 miRNAs and 185 DEGs, including miR-182-5p, miR-590-3p, miR-320a and serum- and glucocorticoid-regulated kinase 1 (SGK1). Then, the related key genes and miRNAs were confirmed by qRT-PCR. Furthermore, 81 small molecules with antagonistic or synergistic effect to GEM were screened. We have investigated the molecular mechanisms driving GEM-resistance in bladder cancer cells that would contribute to the development of chemotherapy for advanced bladder cancer.     

Construction and analysis of circular RNA molecular regulatory networks in liver cancer

Liver cancer is the sixth most prevalent cancer, and the third most frequent cause of cancer-related deaths. Circular RNAs (circRNAs), a kind of special endogenous ncRNAs, have been coming back to the forefront of cancer genomics research. In this study, we used a systems biology approach to construct and analyze the circRNA molecular regulatory networks in the context of liver cancer. We detected a total of 127 differentially expressed circRNAs and 3,235 differentially expressed mRNAs. We selected the top-5 upregulated circRNAs to construct a circRNA-miRNA-mRNA network. We enriched the pathways and gene ontology items and determined their participation in cancer-related pathways such as p53 signaling pathway and pathways involved in angiogenesis and cell cycle. Quantitative real-time PCR was performed to verify the top-five circRNAs. ROC analysis showed circZFR, circFUT8, circIPO11 could significantly distinguish the cancer samples, with an AUC of 0.7069, 0.7575, and 0.7103, respectively. Our results suggest the circRNA-miRNA-mRNA network may help us further understand the molecular mechanisms of tumor progression in liver cancer, and reveal novel biomarkers and therapeutic targets.     

Construction of a competing endogenous RNA network to analyse glucose-6-phosphate dehydrogenase dysregulation in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a common malignant tumour with high rates of morbidity and mortality worldwide. Therefore, it is of great significance to find new molecular markers for HCC diagnosis and treatment. G6PD is known to be dysregulated in a variety of tumours. In addition, the ceRNA network plays a crucial role in the occurrence and development of HCC. However, the mechanism by which the ceRNA network regulates G6PD in HCC remains unclear. We used TCGA-LIHC data to analyse the possibility of using G6PD as an independent prognostic marker. Univariate Cox proportional hazards regression, multivariate Cox proportional hazards regression, and receiver operating characteristic curve analysis were used to analyse the influence of G6PD overexpression on the prognosis of HCC patients. We also analysed the biological function of G6PD, its effect on the immune microenvironment, and drug sensitivity. Finally, we constructed a ceRNA network of lncRNAs/miR-122-5p/G6PD to explore the regulatory mechanism of G6PD. G6PD was highly expressed in HCC, was related to pathological stage and poor prognosis, and could be used as an independent prognostic indicator of HCC. The expression of G6PD was closely related to the immune microenvironment of HCC. In addition, the expression of G6PD in HCC could be regulated by the ceRNA network. Therefore, G6PD can be used as an immunotherapy target to improve the survival and prognosis of HCC patients, and the ceRNA regulatory network of G6PD has potential diagnostic and therapeutic value for HCC.     

microRNA调控网络与肿瘤

microRNA(miRNA)是一类长度约为22nt的小分子非编码RNA,对基因表达进行转录后调控。基因调控网络由各种回路,如前馈和反馈回路组成。它通过将外界刺激整合成合适的输出信号,控制细胞从增殖到死亡几乎所有的生命活动。已知基因调控网络的失调与肿瘤的发生发展密切相关。新近的实验证据表明,miRNA广泛参与基因调控网络中的反馈和前馈回路。该文主要介绍miRNA在基因表达调控网络中的作用及其与肿瘤发生发展的关系。     

Reconstruction and analysis of the aberrant lncRNA‐miRNA‐mRNA network based on competitive endogenous RNA in CESC

A growing body of studies has demonstrated that long non‐coding RNA (lncRNA) are regarded as the primary section of the ceRNA network. This is thought to be the case owing to its regulation of protein‐coding gene expression by functioning as miRNA sponges. However, functional roles and regulatory mechanisms of lncRNA‐mediated ceRNA in cervical squamous cell carcinoma (CESC), as well as their use for potential prediction of CESC prognosis, remains unknown. The aberrant expression profiles of mRNA, lncRNA, and miRNA of 306 cervical squamous cancer tissues and three adjacent cervical tissues were obtained from the TCGA database. A lncRNA‐mRNA‐miRNA ceRNA network in CESC was constructed. Meanwhile, Gene Ontology (GO) and KEGG pathway analysis were performed using Cytoscape plug‐in BinGo and DAVID database. We identified a total of 493 lncRNA, 70 miRNA, and 1921 mRNA as differentially expressed profiles. An aberrant lncRNA‐mRNA‐miRNA ceRNA network was constructed in CESC, it was composed of 50 DElncRNA, 18 DEmiRNA, and 81 DEmRNA. According to the overall survival analysis, 3 out of 50 lncRNA, 10 out of 81 mRNA, and 1 out of 18 miRNA functioned as prognostic biomarkers for patients with CESC (P value < 0.05). We extracted the sub‐network in the ceRNA network and found that two novel lncRNA were recognized as key genes. These included lncRNA MEG3 and lncRNA ADAMTS9‐AS2. The present study provides a new insight into a better understanding of the lncRNA‐related ceRNA network in CESC, and the novel recognized ceRNA network will help us to improve our understanding of lncRNA‐mediated ceRNA regulatory mechanisms in the pathogenesis of CESC.     

Large-scale profiling and identification of potential regulatory mechanisms for allelic gene expression in colorectal cancer cells

Allelic variation in gene expression is common in humans and this variation is associated with phenotypic variation. In this study, we employed high-density single nucleotide polymorphism (SNP) chips containing 13,900 exonic SNPs to identify genes with allelic gene expression in cells from colorectal cancer cell lines. We found 2 monoallelically expressed genes (ERAP2 and MYLK4), 32 genes with an allelic imbalance in their expression, and 13 genes showing allele substitution by RNA editing. Among a total of 34 allelically expressed genes in colorectal cancer cells, 15 genes (44.1%) were associated with cis-acting eQTL, indicating that large portions of allelically expressed genes are regulated by cis-acting mechanisms of gene expression. In addition, potential regulatory variants present in the proximal promoter regions of genes showing either monoallelic expression or allelic imbalance were not tightly linked with coding SNPs, which were detected with allelic gene expression. These results suggest that multiple rare variants could be involved in the cis-acting regulatory mechanism of allelic gene expression. In the comparison with allelic gene expression data from Centre d'Etude du Polymorphisme Humain (CEPH) family B cells, 12 genes showed B-cell specific allelic imbalance and 1 noncoding SNP showed colorectal cancer cell-specific allelic imbalance. In addition, different patterns of allele substitution were observed between B cells and colorectal cancer cells. Overall, our study not only indicates that allelic gene expression is common in colorectal cancer cells, but our study also provides a better understanding of allele-specific gene expression in colorectal cancer cells.