The long noncoding RNA Vax2os1 controls the cell cycle progression of photoreceptor progenitors in the mouse retina

Long noncoding RNAs (lncRNAs) are emerging as regulators of many basic cellular pathways. Several lncRNAs are selectively expressed in the developing retina, although little is known about their functional role in this tissue. Vax2os1 is a retina-specific lncRNA whose expression is restricted to the mouse ventral retina. Here we demonstrate that spatiotemporal misexpression of Vax2os1 determines cell cycle alterations in photoreceptor progenitor cells. In particular, the overexpression of Vax2os1 in the developing early postnatal mouse retina causes an impaired cell cycle progression of photoreceptor progenitors toward their final committed fate and a consequent delay of their differentiation processes. At later developmental stages, this perturbation is accompanied by an increase of apoptotic events in the photoreceptor cell layer, in comparison with control retinas, without affecting the proper cell layering in the adult retina. Similar results are observed in mouse photoreceptor-derived 661W cells in which Vax2os1 overexpression results in an impairment of the cell cycle progression rate and cell differentiation. Based on these results, we conclude that Vax2os1 is involved in the control of cell cycle progression of photoreceptor progenitor cells in the ventral retina. Therefore, we propose Vax2os1 as the first example of lncRNA that acts as a cell cycle regulator in the mammalian retina during development.     

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.     

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.     

Long Non-coding RNA Derived from lncRNA–mRNA Co-expression Networks Modulates the Locust Phase Change

Long non-coding RNAs (lncRNAs) regulate various biological processes ranging from gene expression to animal behavior. Although protein-coding genes, microRNAs, and neuropeptides play important roles in the regulation of phenotypic plasticity in migratory locust, empirical studies on the function of lncRNAs in this process remain limited. Here, we applied high-throughput RNA-seq to compare the expression patterns of lncRNAs and mRNAs in the time course of locust phase change. We found that lncRNAs responded more rapidly at the early stages of phase transition. Functional annotations demonstrated that early changed lncRNAs employed different pathways in isolation and crowding phases to cope with changes in the population density. Two overlapping hub lncRNA loci in the crowding and isolation networks were screened for functional verification. One of them, LNC1010057, was validated as a potential regulator of locust phase change. This work offers insights into the molecular mechanism underlying locust phase change and expands the scope of lncRNA functions in animal behavior.     

Identification of messenger and long noncoding RNAs associated with gallbladder cancer via gene expression profile analysis

The present study aimed to investigate the long noncoding RNAs (lncRNAs) and messenger RNAs (mRNAs) involved in the progression of gallbladder cancer and explore the potential physiopathologic mechanisms of gallbladder cancer in terms of competing endogenous RNAs (ceRNAs). The original lncRNA and mRNA expression profile data (nine gallbladder cancer tissues samples and nine normal gallbladder samples) in GSE76633 was downloaded from the Gene Expression Omnibus database. Differentially expressed mRNAs and lncRNAs between gallbladder cancer tissue and normal control were selected and the pathways in which they are involved were analyzed using bioinformatics analyses. MicroRNAs (miRNAs) were also predicted based on the differentially expressed mRNAs. Finally, the co-expression relation between lncRNA and mRNA was analyzed and the ceRNA network was constructed by combining the lncRNA-miRNA, miRNA-mRNA, and lncRNA-mRNA pairs. Overall, 373 significantly differentially expressed mRNAs and 47 lncRNAs were identified between cancer and normal tissue samples. The upregulated genes were significantly enriched in the extracellular matrix (ECM)-receptor interaction pathway, while the downregulated genes were involved in the complement and coagulation cascades. Altogether, 128 co-expression relations between lncRNA and mRNA were obtained. In addition, 196 miRNA-mRNA regulatory relations and 145 miRNA-lncRNA relation pairs were predicted. Finally, the lncRNA-miRNA-gene ceRNA network was constructed by combining the three types of relation pairs, such as XLOC_011309-miR-548c-3p-SPOCK1 and XLOC_012588-miR-765-CEACAM6. mRNAs and lncRNAs may be involved in gallbladder cancer progression via ECM-receptor interaction pathways and the complement and coagulation cascades. Moreover, ceRNAs such as XLOC_011309-miR-548c-3p-SPOCK1 and XLOC_012588-miR-765-CEACAM6 can also be implicated in the pathogenesis of gallbladder cancer.     

Identification of lncRNA–miRNA–mRNA regulatory network associated with epithelial ovarian cancer cisplatin-resistant

To construct a long noncoding RNA (lncRNA)–microRNA (miRNA)–messenger RNA (mRNA) regulatory network related to epithelial ovarian cancer (EOC) cisplatin-resistant, differentially expressed genes (DEGs), differentially expressed lncRNAs (DELs), and differentially expressed miRNAs (DEMs) between MDAH and TOV-112D cells lines were identified. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to analyze the biological functions of DEGs. Downstream mRNAs or upstream lncRNAs for miRNAs were analyzed at miRTarBase 7.0 or DIANA-LncBase V2, respectively. A total of 485 significant DEGs, 85 DELs, and 5 DEMs were identified. Protein–protein interaction (PPI) network of DEGs contrains 81 nodes and 141 edges was constructed, and 25 hub genes related to EOC cisplatin-resistant were identified. Subsequently, a lncRNA–miRNA–mRNA regulatory network contains 4 lncRNAs, 4 miRNAs, and 35 mRNAs was established. Taken together, our study provided evidence concerning the alteration genes involved in EOC cisplatin-resistant, which will help to unravel the mechanisms underlying drug resistant.     

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.     

Microarray analysis of long-noncoding RNAs and mRNA expression profiles in human steroid-induced avascular necrosis of the femoral head

This study performed the first microarray analysis of long-noncoding RNA (lncRNA) and mRNA expression profiles in human steroid-induced avascular necrosis of the femoral head (SAVNFH). Expression levels of lncRNAs and mRNAs in three human SAVNFH samples and three human femoral head fracture samples (controls) were detected using third-generation lncRNA microarrays (KangChen Biotech, Shanghai, China). The fold change, false discovery rate, and P value were utilized to filter genes with significant differential expression in the SAVNFH samples compared with the control samples. In total, there were 1179 upregulated and 3214 downregulated lncRNAs (P2. zerofold, P < 0.05). Meanwhile, 1092 upregulated and 565 downregulated mRNAs were found in the SAVNFH samples compared with the control samples. Then, quantitative real-time polymerase chain reaction was used to confirm the previous microarray results using 8 and 20 selected dysregulated lncRNAs and mRNAs, respectively, and the results generally confirmed the microarray findings. Finally, we used Gene Ontology (GO) and pathway analysis to investigate the functions of the altered mRNAs and their associated GO terms and biological pathways. The Immune system process term (GO:0002376) was the most significantly upregulated GO term, and the Regulation of blood coagulation term (GO:0030193) was the most significantly downregulated GO term in the biological process category for the SAVNFH samples. “Hematopoietic cell lineage - Homo sapiens (human) (Pathway ID: hsa04640)” and “Complement and coagulation cascades - Homo sapiens (human) (Pathway ID: hsa04610)” were the most significantly up- and downregulated pathways in the SAVNFH samples compared with the controls. In conclusion, the differential expression of lncRNAs and mRNAs may be correlated with the pathogenesis of SAVNFH, and these significantly dysregulated lncRNAs and mRNAs may function through networks or participate in several specific biological processes. Further research is needed to understand their exact functions and mechanisms in SAVNFH.     

Integrative analysis of competing endogenous RNA networks reveals the functional lncRNAs in heart failure

Heart failure has become one of the top causes of death worldwide. It is increasing evidence that lncRNAs play important roles in the pathology processes of multiple cardiovascular diseases. Additionally, lncRNAs can function as ceRNAs by sponging miRNAs to affect the expression level of mRNAs, implicating in numerous biological processes. However, the functional roles and regulatory mechanisms of lncRNAs in heart failure are still unclear. In our study, we constructed a heart failure‐related lncRNA‐mRNA network by integrating probe re‐annotation pipeline and miRNA‐target interactions. Firstly, some lncRNAs that had the central topological features were found in the heart failure‐related lncRNA‐mRNA network. Then, the lncRNA‐associated functional modules were identified from the network, using bidirectional hierarchical clustering. Some lncRNAs that involved in modules were demonstrated to be enriched in many heart failure‐related pathways. To investigate the role of lncRNA‐associated ceRNA crosstalks in certain disease or physiological status, we further identified the lncRNA‐associated dysregulated ceRNA interactions. And we also performed a random walk algorithm to identify more heart failure‐related lncRNAs. All these lncRNAs were verified to show a strong diagnosis power for heart failure. These results will help us to understand the mechanism of lncRNAs in heart failure and provide novel lncRNAs as candidate diagnostic biomarkers or potential therapeutic targets.     

Analysis of long noncoding RNA expression profiles in the whole blood of neonates with hypoxic-ischemic encephalopathy

Long noncoding RNAs (lncRNAs) have been shown to participate in many biological processes. To investigate the expression profiles of lncRNAs and their potential functions in neonatal hypoxic-ischemic encephalopathy (HIE), we detected the lncRNA and messenger RNA (mRNA) expression in the peripheral blood samples from HIE patients and controls using a microarray. A total of 376 lncRNAs and 126 mRNAs were differentially expressed between the HIE and the non-HIE samples (fold change > 2). Quantitative real-time polymerase chain reaction was used to validate the microarray data. Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analysis were performed to determine the gene function. Furthermore, the lncRNA-mRNA coexpression network was generated to predict the potential targets of lncRNAs. In conclusion, our study first demonstrated the differential expression profiles of lncRNAs in the whole blood of infants with HIE and may provide a new view of the distinct lncRNA functions in HIE.     

ZEA exerts toxic effects on reproduction and development by mediating Dio3os in mouse endometrial stromal cells

Zearalenone (ZEA) and imprinted long noncoding RNAs (lncRNAs) are both closely related to reproduction and development. However, whether they have connections in regulating reproduction and development is not clear yet. The aim of this research is to investigate their relationship. lncRNA microarray was performed to analyze differentially expressed genes, and real‐time quantitative polymerase chain reaction (PCR) was used to verify the accuracy of microarray analysis. Meanwhile, the technologies of rapid amplification of cDNA ends, RNA fluorescence in situ hybridization and bioinformatics were adopted to characterize the selected lncRNA. Analysis of lncRNA microarray showed lncRNAs and messenger RNAs related to reproduction and development were significantly differently expressed, and Dio3os was probably the target lncRNA. Then, the experiment of real‐time quantitative PCR verified the accuracy of microarray data. Characterization of Dio3os showed Dio3os, an antisense lncRNA with 2312 bp and 15 open reading frames, was enriched in the cytoplasm. Our findings suggest ZEA probably exerts toxic effects on reproduction and development by mediating Dio3os.     

Time-Series Clustering of lncRNA-mRNA Expression during the Adipogenic Transdifferentiation of Porcine Skeletal Muscle Satellite Cells

Skeletal muscle satellite cells (SMSCs), which are multifunctional muscle-derived stem cells, can differentiate into adipocytes. Long-chain non-coding RNA (lncRNA) has diverse biological functions, including the regulation of gene expression, chromosome silencing, and nuclear transport. However, the regulatory roles and mechanism of lncRNA during adipogenic transdifferentiation in muscle cells have not been thoroughly investigated. Here, porcine SMSCs were isolated, cultured, and induced for adipogenic differentiation. The expressions of lncRNA and mRNA at different time points during transdifferentiation were analysed using RNA-seq analysis. In total, 1005 lncRNAs and 7671 mRNAs showed significant changes in expression at differential differentiation stages. Time-series expression analysis showed that the differentially expressed (DE) lncRNAs and mRNAs were clustered into 5 and 11 different profiles with different changes, respectively. GO, KEGG, and REACTOME enrichment analyses revealed that DE mRNAs with increased expressions during the trans-differentiation were mainly enriched in the pathways for lipid metabolism and fat cell differentiation. The genes with decreased expressions were mainly enriched in the regulation of cell cycle and genetic information processing. In addition, 1883 DE mRNAs were regulated by 193 DE lncRNAs, and these genes were related to the controlling in cell cycle mainly. Notably, three genes in the fatty acid binding protein (FABP) family significantly and continuously increased during trans-differentiation, and 15, 13, and 11 lncRNAs may target FABP3, FABP4, and FABP5 genes by cis- or trans-regulation, respectively. In conclusion, these studies identify a set of new potential regulator for adipogenesis and cell fate and help us in better understanding the molecular mechanisms of trans-differentiation.     

Comprehensive analysis of a ceRNA network reveals potential prognostic cytoplasmic lncRNAs involved in HCC progression

The aberrant expression of long noncoding RNAs (lncRNAs) has drawn increasing attention in the field of hepatocellular carcinoma (HCC) biology. In the present study, we obtained the expression profiles of lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) in 371 HCC tissues and 50 normal tissues from The Cancer Genome Atlas (TCGA) and identified hepatocarcinogenesis-specific differentially expressed genes (DEGs, log fold change ≥ 2, FDR < 0.01), including 753 lncRNAs, 97 miRNAs, and 1,535 mRNAs. Because the specific functions of lncRNAs are closely related to their intracellular localizations and because the cytoplasm is the main location for competitive endogenous RNA (ceRNA) action, we analyzed not only the interactions among these DEGs but also the distributions of lncRNAs (cytoplasmic, nuclear or both). Then, an HCC-associated deregulated ceRNA network consisting of 37 lncRNAs, 10 miRNAs, and 26 mRNAs was constructed after excluding those lncRNAs located only in the nucleus. Survival analysis of this network demonstrated that 15 lncRNAs, 3 miRNAs, and 16 mRNAs were significantly correlated with the overall survival of HCC patients (p < 0.01). Through multivariate Cox regression and lasso analysis, a risk score system based on 13 lncRNAs was constructed, which showed good discrimination and predictive ability for HCC patient survival time. This ceRNA network-construction approach, based on lncRNA distribution, not only narrowed the scope of target lncRNAs but also provided specific candidate molecular biomarkers for evaluating the prognosis of HCC, which will help expand our understanding of the ceRNA mechanisms involved in the early development of HCC.     

Alterations of long noncoding RNAs and mRNAs in extracellular vesicles derived from the murine heart post‐ischemia–reperfusion injury

Extracellular vesicles (EVs) play important roles in cardiovascular diseases by delivering their RNA cargos. However, the features and possible role of the lncRNAs and mRNAs in cardiac EVs during ischemia–reperfusion (IR) remain unclear. Therefore, we performed RNA sequencing analysis to profile the features of lncRNAs and mRNAs and predicted their potential functions. Here, we demonstrated that the severity of IR injury was significantly correlated with cardiac EV production. RNA sequencing identified 73 significantly differentially expressed (DE) lncRNAs (39 upregulated and 34 downregulated) and 720 DE‐mRNAs (317 upregulated and 403 downregulated). Gene Ontology (GO) and pathway analysis were performed to predict the potential functions of the DE‐lncRNAs and mRNAs. The lncRNA‐miRNA‐mRNA ceRNA network showed the possible functions of DE‐lncRNAs with DE‐mRNAs which are enriched in the pathways of T cell receptor signalling pathway and cell adhesion molecules. Moreover, the expressions of ENSMUST00000146010 and ENSMUST00000180630 were negatively correlated with the severity of IR injury. A significant positive correlation was revealed between TCONS_00010866 expression and the severity of the cardiac injury. These findings revealed the lncRNA and mRNA profiles in the heart derived EVs and provided potential targets and pathways involved in cardiac IR injury.     

Integrated analysis of a competing endogenous RNA network reveals key long noncoding RNAs as potential prognostic biomarkers for hepatocellular carcinoma

Growing evidence has revealed that long noncoding RNAs (lncRNAs) have an important impact on tumorigenesis and tumor progression via a mechanism involving competing endogenous RNAs (ceRNAs). However, their use in predicting the survival of a patient with hepatocellular carcinoma (HCC) remains unclear. The aim of this study was to develop a novel lncRNA expression–based risk score system to accurately predict the survival of patients with HCC. In our study, using expression profiles downloaded from The Cancer Genome Atlas database, the differentially expressed messenger RNAs (mRNAs), lncRNAs, and microRNAs (miRNAs) were explored in patients with HCC and normal liver tissues, and then a ceRNA network constructed. A risk score system was established between lncRNA expression of the ceRNA network and overall survival (OS) or recurrence-free survival (RFS); it was further analyzed for associations with the clinical features of patients with HCC. In HCC, 473 differentially expressed lncRNAs, 63 differentially expressed miRNAs, and 1417 differentially expressed mRNAs were detected. The ceRNA network comprised 41 lncRNA nodes, 12 miRNA nodes, 24 mRNA nodes, and 172 edges. The lncRNA expression–based risk score system for OS was constructed based on six lncRNAs (MYLK-AS1, AL359878.1, PART1, TSPEAR-AS1, C10orf91, and LINC00501), while the risk score system for RFS was based on four lncRNAs (WARS2-IT1, AL359878.1, AL357060.1, and PART1). Univariate and multivariate Cox analyses showed the risk score systems for OS or RFS were significant independent factors adjusted for clinical factors. Receiver operating characteristic curve analysis showed the area under the curve for the risk score system was 0.704 for OS, and 0.71 for RFS. Our result revealed a lncRNA expression–based risk score system for OS or RFS can effectively predict the survival of patients with HCC and aid in good clinical decision-making.