Long noncoding RNA MEG3 play an important role in osteosarcoma development through sponging microRNAs

More and more evidence indicate long noncoding RNAs (lncRNAs) as competing endogenous RNAs (ceRNAs) to indirectly regulate messenger RNAs (mRNAs) by acting as microRNA (miRNA) sponges, which represents a novel layer of gene regulation that plays a critical role in the development of cancers. However, functional roles and regulatory mechanisms of lncRNA-mediated ceRNAs network in osteosarcoma are still largely unknown. Here, we comprehensively compared the expression profiles of mRNAs, lncRNAs, and miRNAs between osteosarcoma and normal samples from the Gene Expression Omnibus (GEO) to elaborate related latent mechanisms. Two lncRNAs, ie, LINC01560 and MEG3, were identified to be aberrantly expressed. Importantly, MEG3 was considered as a promising diagnostic biomarker and therapeutic target for patients with osteosarcoma according to the Kaplan-Meier analysis of another independent osteosarcoma data set from the Cancer Genome Atlas (P = 0.05). Eventually, we successfully established a dysregulated lncRNA-related ceRNA network, including one osteosarcoma-specific lncRNA, three miRNAs and four mRNAs. In conclusion, this study should be beneficial for improving our understanding of the lncRNA-mediated ceRNA regulatory mechanisms in the pathogenesis of osteosarcoma and providing it with novel candidate diagnostic and therapeutic biomarkers.     

Expression profile of long non-coding RNAs in cervical spondylotic myelopathy of rats by microarray and bioinformatics analysis

IntroductionIt has been reported that a wide range of long non-coding RNAs (lncRNAs) are implicated in numerous diseases such as tumor, cardiopathy and neurological disorders. Identifying the differentially expressed (DE) profile of lncRNAs in cervical spondylotic myelopathy (CSM) is essential to understand the mechanisms of CSM.MethodsMicroarray assay, quantitative real-time PCR (qRT-PCR) and bioinformatics analysis were employed to reveal the DE profile and potential functions of lncRNAs in CSM.ResultsMicroarray analysis displayed the DE profiles of lncRNAs and mRNAs in rats between the CSM group and the control (CON) group. Thereinto, 1266 DE lncRNAs (738 up-regulation and 528 down-regulation) and 847 mRNAs (487 up-regulation and 360 down-regulation) with >1.1 fold change (FC) were finally identified. Moreover, 17 lncRNAs (13 up-regulation and 4 down-regulation) and 18 mRNAs (13 up-regulation and 5 down-regulation) were found deregulated by >2 FC. Further bioinformatics analysis showed the most remarkable biological processes among up-regulated RNAs contain cellular response to interferon-beta, inflammatory response and innate immune response, which may involve in CSM. Besides, related DE mRNAs of 17 DE lncRNAs in the genome were related to signaling pathway about NOD-like receptor, TNF, and apoptosis. In addition, a co-expression network of lncRNA-mRNA was established for analyzing the biological roles of lncRNAs. Among these, we found a ceRNA network related to CSM. Finally, the expressions of the DE lncRNAs and ceRNA network confirmed by qRT-PCR were in agreement with microarray data.ConclusionsOur study revealed the DE profiles of lncRNAs and mRNAs for CSM. Those dysregulated RNAs may represent potential therapeutic targets of CSM for further study.     

Genome-wide differential expression of long noncoding RNAs and mRNAs in ovarian follicles of two different chicken breeds

Bashang long-tail chickens are an indigenous breed with dual purpose in China (meat and eggs) but have low egg laying performance. To improve the low egg laying performance, a genome-wide analysis of mRNAs and long noncoding RNAs (lncRNAs) from Bashang long-tail chickens and Hy-Line brown layers was performed. A total of 16,354 mRNAs and 8691 lncRNAs were obtained from ovarian follicles. Between the breeds, 160 mRNAs and 550 lncRNAs were found to be significantly differentially expressed. Integrated network analysis suggested some differentially expressed genes were involved in ovarian follicular development through oocyte meiosis, progesterone-mediated oocyte maturation, and cell cycle. The impact of lncRNAs on cis and trans target genes, indicating some lncRNAs may play important roles in ovarian follicular development. The current results provided a catalog of chicken ovarian follicular lncRNAs and genes for further study to understand their roles in regulation of egg laying performance.     

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.     

Orchestrating a biomarker panel with lncRNAs and mRNAs for predicting survival in pancreatic ductal adenocarcinoma

The low survival of patients with pancreatic ductal adenocarcinoma (PDAC) makes the treatment of this disease one of the most challenging task in modern medicine. Here, by mining a large‐scale cancer genome atlas data set of pancreatic cancer tissues, we identified 21 long noncoding RNAs (lncRNAs) that significantly associated with overall survival in patients with PDAC (P < .01). Further analysis revealed that 8 lncRNAs turned out to be independently correlated with patients’ overall survival, and the risk score could be calculated based on their expression. To obtain a better predicting power, we integrated lncRNA data with a total of 410 differently expressed messenger RNAs (mRNAs) screened from PDAC and normal tissues in gene expression omnibus (GEO) database. The integration resulted in a much better panel including 8 lncRNAs (RP3.470B24.5, CTA.941F9.9, RP11.557H15.3, LINC00960, AP000479.1, LINC00635, LINC00636, and AC073133.1) and 8 mRNAs (DHRS9, ONECUT1, OR8D4, MT1M, TCN1, MMP9, DPYSL3, and TTN) to predict prognosis. A functional evaluation showed that these lncRNAs might play roles in pancreatic secretion, cell adhesion, and proteolysis. Using normal and pancreatic cancer cell lines, we confirmed that a majority of identified lncRNAs and mRNAs showed altered expressions in pancreatic cancer cells. Especially, LINC01589, LINC00960, TCN1, and MT1M showed a profoundly increased expression in pancreatic cancer cells, which suggests their potentially important role in pancreatic cancer. The results of our work indicate that lncRNAs have vital roles in PADC and provide new insights to integrate multiple kinds of markers in clinical practices.     

Genome-wide identification and characterization of putative lncRNAs in the diamondback moth, Plutella xylostella (L.)

Long non-coding RNAs (lncRNAs) are of particular interest because of their contributions to many biological processes. Here, we present the genome-wide identification and characterization of putative lncRNAs in a global insect pest, Plutella xylostella. A total of 8096 lncRNAs were identified and classified into three groups. The average length of exons in lncRNAs was longer than that in coding genes and the GC content was lower than that in mRNAs. Most lncRNAs were flanked by canonical splice sites, similar to mRNAs. Expression profiling identified 114 differentially expressed lncRNAs during the DBM development and found that majority were temporally specific. While the biological functions of lncRNAs remain uncharacterized, many are microRNA precursors or competing endogenous RNAs involved in micro-RNA regulatory pathways. This work provides a valuable resource for further studies on molecular bases for development of DBM and lay the foundation for discovery of lncRNA functions in P. xylostella.     

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.     

Identification of functional lncRNAs in gastric cancer by integrative analysis of GEO and TCGA data

Gastric cancer (GC) is a prevalent malignant cancer of digestive system, identification of novel diagnostic and prognostic biomarkers for GC is urgently demanded. The aim of this study was to determine potential long noncoding RNAs (lncRNAs) associated with the pathogenesis and prognosis of GC. Raw noncoding RNA microarray data (GSE53137, GSE70880, and GSE99417) was downloaded from Gene Expression Omnibus (GEO) database. Differentially expressed genes between GC and adjacent normal gastric tissue samples were screened by an integrated analysis of multiple gene expression profile after gene reannotation and batch normalization. Differentially expressed genes were further confirmed by The Cancer Genome Atlas (TCGA) database. Competing endogenous RNA (ceRNA) network, Gene Ontology term and Kyoto Encyclopedia of Genes and Genomes pathway, survival analysis were extensively applied to identify hub lncRNAs and discover potential biomarkers related to diagnosis and prognosis of GC. In total of 246 integrated differential genes including 15 lncRNAs and 241 messenger RNAs (mRNAs) were obtained after intersections of differential genes between GEO and TCGA database. ceRNA network comprised of three lncRNAs (UCA1, HOTTIP, and HMGA1P4), 26 microRNAs (miRNAs) and 72 mRNAs. Functional analysis revealed that three lncRNAs were mainly dominated in cell cycle and cellular senescence. Survival analysis showed that HMGA1P4 was statistically related to the overall survival rate. For the first time, we identified that HMGA1P4, a target of miR-301b/miR-508, is involved in cell cycle and senescence process by regulating CCNA2 in GC. Finally, the expression levels of three lncRNAs were validated to be upregulated in GC tissues. Thus, three lncRNAs including UCA1, HOTTIP, and HMGA1P4 may contribute to GC development and their potential functions might be associated with the prognosis of GC.     

Expression profiling of long noncoding RNAs associated with vasculogenic mimicry in osteosarcoma

Osteosarcoma (OS) is the most common highly malignant bone tumor in teens. Vasculogenic mimicry (VM) is defined as de novo extracellular matrix-rich vascular-like networks formed by highly aggressive tumor cells. We previously reported the presence of VM and it is an unfavorable prognostic factor in OS patients. Long noncoding RNAs (lncRNAs) are aberrantly expressed in OS and involved in cancer cell VM. However, lncRNAs in VM formation of OS have not been investigated. We, therefore, profiled the expression of lncRNAs in highly aggressive OS cell line 143B compared with its parental poorly aggressive cell line HOS. The differentially expressed (DE) lncRNAs and messenger RNA (mRNAs) were subjected to constructed lncRNA-mRNA coexpressed network. The top-ranked hub gene lncRNA n340532 knockdown 143B cells were used for in vitro and in vivo VM assays. The annotation of DE lncRNAs was performed according to the coexpressed mRNAs by Gene Ontology and pathway analysis. A total of 1360 DE lncRNAs and 1353 DE mRNAs were screened out. lncRNA MALAT1 and FTX, which have known functions related to VM formation and tumorigenesis were identified in our data. The coexpression network composed of 226 lncRNAs and 118 mRNAs in which lncRNA n340532 had the highest degree number. lncRNA n340532 knockdown reduced VM formation in vitro. The suppression of n340532 also exhibited potent anti-VM and antimetastasis effect in vivo, suggesting its potential role in OS VM and metastasis. Furthermore, n340532 coexpressed with 10 upregulation mRNAs and 3 downregulation mRNAs. The enriched transforming growth factor-β signaling pathway, angiogenesis and so forth were targeted by those coexpressed mRNAs, implying n340532 may facilitate VM formation in OS through these pathways and gene functions. Our findings provide evidence for the potential role of lncRNAs in VM formation of OS that could be used in the clinic for anti-VM therapy in OS.     

The exploration of new therapeutic targets for HPV-negative head and neck squamous cell cancer through the construction of a ceRNA network and immune microenvironment analysis

Previous studies have shown that human papillomavirus (HPV)-negative patients with head and neck squamous cell cancer (HNSCC) suffer from an unsatisfactory prognosis. Long noncoding RNAs (lncRNAs) have been verified to participate in many biological processes, including regulating gene expression as competing endogenous RNAs (ceRNAs), while few studies focused the ceRNA network regulation mechanism in patients with HPV-negative HNSCC tumor. Meanwhile, the immune microenvironment may be critical in the development and prognosis of HPV-negative tumors. Our study aimed to further investigate the pathogenesis and potential biomarkers for the diagnosis, therapy and prognosis of HPV-negative HNSCC through a ceRNA network. Comprehensively analyzing the sequencing data of lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) in The Cancer Genome Atlas HNSCC dataset, we constructed a differentially expressed ceRNA network containing 131 lncRNAs, 35 miRNAs and 162 mRNAs. Then, survival analysis in the network was cited to explore the prognostic biomarkers. Eight mRNAs, nine lncRNAs, and one miRNA were identified to be associated with prognosis. Neuropilin (NRP) binding function, retinoid X receptor (RXR) binding, and the vascular endothelial growth factor (VEGF) signaling pathway were associated with the enrichment analysis, and they also related to the immune microenvironment. Combined with the analysis of the immune microenvironment differences, we obtained new targeted therapies using an RXR agonist, or a combination of the VEGF monoclonal antibody and an NRP antagonist, which may provide a promising future for HPV-negative HNSCC patients.     

Aberrantly expressed long non‐coding RNAs in air pollution‐induced congenital defects

Air pollution has been a serious public health issue over the past few decades particularly in developing countries. Air pollution exposure during pregnancy poses potential threat to offspring as the deleterious substances might pass through placenta to alter foetal development. A growing number of studies have demonstrated that long non‐coding RNAs (lncRNAs) participate in the development of many diseases, including congenital defects. Here, we used RNA sequencing to identify differentially expressed lncRNAs in air pollution‐exposed rat embryos compared with control group. Our data suggested that 554 lncRNAs (216 up‐regulated and 338 down‐regulated) were significantly differentially expressed in the air pollution‐exposed embryos. Moreover, potential cellular functions of these deregulated lncRNAs were predicted via KEGG signal pathway/GO enrichment analyses, which suggested the possible involvements of neurological process, sensory perception of smell and the G‐protein signalling pathway. Furthermore, potential functional network of deregulated lncRNAs and their correlated mRNAs in the development of congenital spinal abnormality was established. Our data suggested that lncRNAs may play a vital role in the pathophysiology of air pollution‐exposed congenital spinal malformation.     

An 11-lncRNA expression could be potential prognostic biomarkers in head and neck squamous cell carcinoma

The aim of our study is to construct the competing endogenous RNA (ceRNA) network of head and neck squamous cell carcinoma (HNSCC) and identify key long noncoding RNAs (lncRNAs) to predict prognosis. The genes whose expression were differentially in HNSCC and normal tissues were explored by the Cancer Genome Atlas database. The ceRNA network was constructed by the Cytoscape software. The lncRNAs which could estimate the overall survival were explored from Cox proportional hazards regression. There are 1997, 589, and 82 mRNAs, lncRNAs, and miRNAs whose expression were statistically significant different, respectively. Then, the network between miRNA and mRNA or miRNA and lncRNA was constructed by miRcode, miRDB, TargetScan, and miRanda. Five mRNAs, 10 lncRNAs, and 3 miRNAs were associated with overall survival. Then, 11-lncRNAs were found to be prognostic factors. Therefore, our research analyzed the potential signature of novel 11-lncRNA as candidate prognostic biomarker from the ceRNA network for patients with HNSCC.     

Identification of functional lncRNAs in atrial fibrillation by integrative analysis of the lncRNA-mRNA network based on competing endogenous RNAs hypothesis

A mounting body of evidence has suggested that long noncoding RNAs (lncRNAs) play critical roles in human diseases by acting as competing endogenous RNAs (ceRNAs). However, the functions and ceRNA mechanisms of lncRNAs in atrial fibrillation (AF) remain to date unclear. In this study, we constructed an AF-related lncRNA-mRNA network (AFLMN) based on ceRNA theory, by integrating probe reannotation pipeline and microRNA (miRNA)-target regulatory interactions. Two lncRNAs with central topological properties in the AFLMN were first obtained. By using bidirectional hierarchical clustering, we identified two modules containing four lncRNAs, which were significantly enriched in many known pathways of AF. To elucidate the ceRNA interactions in certain disease or normal condition, the dysregulated lncRNA-mRNA crosstalks in AF were further analyzed, and six hub lncRNAs were obtained from the network. Furthermore, random walk analysis of the AFLMN suggested that lncRNA RP11-296O14.3 may function importantly in the pathological process of AF. All these eight lncRNAs that were identified from previous steps (RP11-363E7.4, GAS5, RP11-410L14.2, HAGLR, RP11-421L21.3, RP11-111K18.2, HOTAIRM1, and RP11-296O14.3) exhibited a strong diagnostic power for AF. The results of our study provide new insights into the functional roles and regulatory mechanisms of lncRNAs in AF, and facilitate the discovery of novel diagnostic biomarkers or therapeutic targets.