Differential circular RNAs expression in ovary during oviposition in honey bees

Circular RNAs (circRNAs) are non-coding RNAs newly identified and play important roles in RNA regulation. The mechanism and function of circRNAs have been reported in some species. However, little is known regarding circRNAs in honey bees. In this study, we analyzed circRNAs through bioinformatics, and predicted 12,211 circRNAs in the ovary of honey bee queens. 1340, 175 and 100 circRNAs were differentially expressed in comparisons of egg-laying queens vs virgin queens, egg-laying inhibited queens vs egg-laying queens and egg-laying recovery queens vs egg-laying inhibited queens. Further, functional annotation of differentially expressed circRNAs revealed several pathways that are closely related to ovary activation and oviposition, including insulin secretion and calcium signaling pathways. Moreover, the potential interactions among circRNAs, miRNAs, lncRNAs and mRNAs were investigated. Ame_circ_0005197 and ame_circ_0016640 were observed to sponge several reproductive related miRNAs. These findings demonstrate that circRNAs have potential effects in ovary activation and oviposition of honey bees.     

Systematic identification of lncRNAs and circRNAs-associated ceRNA networks in human lumbar disc degeneration

ABSTRACT

Lumbar disc degeneration (LDD) is a common cause of low back and neck pain. The molecular mechanisms underlying LDD, however, are unclear. Noncoding RNAs have been reported to participate in human diseases. We investigated a series of public datasets (GSE67566, GSE56081 and GSE63492) and identified 568 mRNAs, 55 microRNAs (miRNAs), 765 long noncoding RNAs (lncRNAs), and 586 circular RNAs (circRNAs) that were expressed differently in LDD than in normal discs. We constructed lncRNAs and circRNAs regulated competing endogenous RNAs (ceRNA) networks in LDD. Four lncRNAs, DANCR, CASK-AS1, SCARNA2, and LINC00638), and three circRNAs, hsa_circ_0005139, hsa_circ_0037858, and hsa_circ_0087890, were identified as key regulators of LDD progression. We found that hsa-miR-486-5p regulated the crosstalk among circRNA hsa_circ_0000189, lncRNA DANCR and 6 mRNAs, PYCR2, TOB1, ARHGAP5, RBPJ, CD247, SLC34A1. Gene ontology (GO) analysis demonstrated that these differently expressed lncRNAs and circRNAs were involved in cellular component organization or biogenesis, gene expression and negative regulation of metabolic processes. Our findings provide useful information for exploring new mechanisms for LDD and candidates for therapeutic targets.     

The novel circular RNA CircMef2c is positively associated with muscle growth in Nile tilapia

Muscle growth in teleosts is a complex biological process orchestrated by numerous protein-coding genes and non-coding RNAs. A few recent studies suggest that circRNAs are involved in teleost myogenesis, but the molecular networks involved remain poorly understood. In this study, an integrative omics approach was used to determine myogenic circRNAs in Nile tilapia by quantifying and comparing the expression profile of mRNAs, miRNAs, and circRNAs in fast muscle from full-sib fish with distinct growth rates. There were 1947 mRNAs, 9 miRNAs, and 4 circRNAs differentially expressed between fast- and slow-growing individuals. These miRNAs can regulate myogenic genes and have binding sites for the novel circRNA circMef2c. Our data indicate that circMef2c may interact with three miRNAs and 65 differentially expressed mRNAs to form multiple competing endogenous RNA networks that regulate growth, thus providing novel insights into the role of circRNAs in the regulation of muscle growth in teleosts.     

Construction of lncRNA-associated ceRNA networks to identify prognostic lncRNA biomarkers for glioblastoma

Long noncoding RNAs (lncRNAs) serve as competitive endogenous RNAs (ceRNAs) that play significant regulatory roles in the pathogenesis of tumors. However, the role of lncRNAs, especially the lncRNA-related ceRNA regulatory network, in glioblastoma (GBM) has not been fully elucidated. The goal of the current study was to construct lncRNA-microRNA-mRNA-related ceRNA networks for further investigation of their mechanism of action in GBM. We downloaded data from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases and identified differential lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) associated with GBM. A ceRNA network was constructed and analyzed to examine the relationship between lncRNAs and patients’ overall survival. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGGs) were used to analyze the related mRNAs to indirectly explain the mechanism of action of lncRNAs. The potential effective drugs for the treatment of GBM were identified using the connectivity map (CMap). After integrated analysis, we obtained a total of 210 differentially expressed lncRNAs, 90 differentially expressed miRNAs, and 2508 differentially expressed mRNAs (DEmRNAs) from the TCGA and GEO databases. Using these differential genes, we constructed a lncRNA-associated ceRNA network. Six lncRNAs in the ceRNA network were associated with the overall survival of patients with GBM. Through KEGG analysis, it was found that the DEmRNAs involved in the network are related to cancer-associated pathways, for instance, mitogen-activated protein kinase and Ras signaling pathways. CMap analysis revealed four small-molecule compounds that could be used as drugs for the treatment of GBM. In this study, a multi-database joint analysis was used to construct a lncRNA-related ceRNA network to help identify the regulatory functions of lncRNAs in the pathogenesis of GBM.     

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.     

circRNA-associated ceRNA network construction reveals the circRNAs involved in the progression and prognosis of breast cancer

Recently, increasing evidences show that circular RNAs (circRNAs) are important regulators of various diseases, especially cancer. However, the regulatory role and the potential mechanism of action of circRNAs in breast cancer remain largely unknown. In this study, weighted gene co-expression network analysis was conducted with the differentially expressed miRNAs and mRNAs in breast cancer from The Cancer Genome Atlas database to identify the key modules associated with the carcinogenesis of breast cancer. In the significant turquoise and brown modules, 22 miRNAs and 1877 mRNAs were identified, respectively. Then, We compared and predicted the target genes and performed survival analysis to identify the miRNAs and mRNAs related to the prognosis of breast cancer. A circRNA-related competitive endogenous RNA network was identified by database co-screening, and deleted in liver cancer 1 (DLC1) was identified as a key gene. Finally, to assess how genes in key modules and key genes contribute to the development of breast cancer, relevant pathway information was obtained through DAVID and Gene Set Enrichment Analysis. These data demonstrated that three circRNAs (hsa-circ-0083373, hsa-circ-0083374, and hsa-circ-0083375) that regulate DLC1 expression via hsa-mir-511 and are involved in the pathogenesis and development of breast cancer.     

Comprehensive analysis of lncRNA–miRNA–mRNA during proliferative phase of rat liver regeneration

This study aims to reveal the regulatory mechanism of lncRNAs–miRNAs–mRNAs network during the proliferative phase of liver regeneration (LR). High-throughput sequencing technology was performed, and a total of 1,738 differentially expressed lncRNAs (DE lncRNAs), 167 known differentially expressed miRNAs (DE miRNAs), and 2,727 differentially expressed mRNAs were identified. Then, the target DE lncRNAs and DE mRNAs regulated by the same miRNAs were screened and a ceRNA regulatory network containing 32 miRNAs, 107 lncRNAs, and 270 mRNAs was constructed. Insulin signaling pathway, pyrimidine metabolism, axon guidance, carbohydrate digestion and absorption, and pyruvate metabolism were significantly enriched in the network. Through literature review and the regulatory relationship between lncRNAs and miRNAs, nine core lncRNAs were identified, which might play important roles during the proliferative phase of rat LR. This study analyzed lncRNA–miRNA–mRNA regulatory network for the first time during the proliferative phase of rat LR, providing clues for exploring the mechanism of LR and the treatment of liver diseases.     

Screening and bioinformatics analysis of mRNA,long non-coding RNA and circular RNA expression profiles in mucoepidermoid carcinoma of salivary gland

Mucoepidermoid carcinoma (MEC) of salivary gland is a disease characterized by high rate of diatant metastasis, and associated with poor outcomes. However, the molecular mechanisms underlying the MEC remain poorly understand. Here, we simultaneously detected, for the first time, the expression profiles of mRNAs, lncRNAs, and circRNAs in four pairs of MEC and matched non-carcinoma tissues by microarrays. A total of 3612?mRNA, 3091 lncRNAs, and 284 circRNAs were altered during the pathogenesis of MEC. The functions of these differentially expressed RNAs were predicted by Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Co-expression networks of lncRNA-mRNA and circRNA-miRNA were conducted to uncovered the hidden ceRNA mechanisms. Moreover, NONHSAT154433.1 that associated with ADAM12 and hsa_circ_0012342 were further screened and confirmed using qRT-PCR analysis. In conclusion, this study provides a systematic perspective on the potential function of non-coding RNAs (ncRNAs) in the molecular mechanisms of MEC. Among these, NONHSAT154433.1 and hsa_circ_0012342 might be served as potential prognostic biomarkers and therapeutic target of MEC.