Identification of mitochondrial function‐associated lncRNAs in septic mice myocardium

The present study aimed to analyze long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression profiles in septic mice heart and to identify potential lncRNAs and mRNAs that be responsible for cardiac mitochondrial dysfunction during sepsis. Mice were treated with 10 mg/kg of lipopolysaccharides to induce sepsis. LncRNAs and mRNAs expression were evaluated by using lncRNA and mRNA microarray or real‐time polymerase chain reaction technique. LncRNA‐mRNA coexpression network assay, Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed. The results showed that 1275 lncRNAs were differentially expressed in septic myocardium compared with those in the control group. A total of 2769 mRNAs were dysregulated in septic mice heart, most of which are mainly related to the process of inflammation, mitochondrial metabolism, oxidative stress, and apoptosis. Coexpression network analysis showed that 14 lncRNAs were highly correlated with 11 mitochondria‐related differentially expressed mRNA. Among all lncRNAs and their cis‐acting mRNAs, 41 lncRNAs‐mRNA pairs (such as NONMMUG004378 and Apaf1 gene) were enriched in GO terms and KEGG pathways. In summary, we gained some specific lncRNAs and their potential target mRNAs that might be involved in mitochondrial dysfunction in septic myocardium. These findings provide a panoramic view of lncRNA and might allow developing new treatment strategies for sepsis.     

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.     

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.     

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.     

Long noncoding RNA HOXA-AS2 regulates the expression of SCN3A by sponging miR-106a in breast cancer

Breast cancer is the most commonly diagnosed cancer that affects women worldwide. This study aimed to investigate the competing endogenous RNAs (ceRNAs) mechanism in breast cancer. Microarray data were downloaded from the University of California Santa Cruz (UCSC) Xena database. The limma package was used to screen the differentially expressed messenger RNAs (DEMs) and differentially expressed long noncoding RNAs (DELs). Subsequently, functional analysis was performed using DAVID tool. After constructing the protein-protein interaction (PPI) network, we identified the major gene modules using the Cytoscape software. Univariate survival analysis in the survival package was performed. Finally, the ceRNA regulatory network was constructed to identify the critical genes. A total of 1380 DEMs and 345 DELs were identified in breast cancer samples compared with normal samples. Functional enrichment analysis showed that DEMs were mainly involved in cell division, and cell cycle. We screened four major gene modules and identified the hub nodes in these functional modules. Several DEMs (including FABP7, C4BPA, and LAMB3) and three long noncoding RNAs (lncRNAs) (LINC00092, SLC26A4.AS1, and COLCA1) exhibited significant correlation with patients' survival outcomes. In the ceRNA network, the lncRNA HOXA-AS2 regulated the expression level of SCN3A by interacting with hsa-miR-106a-5p. Thus, our study investigated the ceRNA mechanism in breast cancer. The results showed that lncRNA HOXA-AS2 might modulate the expression of SCN3A by sponging miR-106a in breast cancer.     

Understanding pathogen–host interplay by expression profiles of lncRNA and mRNA in the liver of Echinococcus multilocularis-infected mice

Almost all Echinococcus multilocularis (Em) infections occur in the liver of the intermediate host, causing a lethal zoonotic helminthic disease, alveolar echinococcosis (AE). However, the long non-coding RNAs (lncRNAs) expression profiles of the host and the potential regulatory function of lncRNA during Em infection are poorly understood. In this study, the profiles of lncRNAs and mRNAs in the liver of mice at different time points after Em infection were explored by microarray. Thirty-one differentially expressed mRNAs (DEMs) and 68 differentially expressed lncRNAs (DELs) were found continuously dysregulated. These DEMs were notably enriched in “antigen processing and presentation”, “Th1 and Th2 cell differentiation” and “Th17 cell differentiation” pathways. The potential predicted function of DELs revealed that most DELs might influence Th17 cell differentiation and TGF-β/Smad pathway of host by trans-regulating SMAD3, STAT1, and early growth response (EGR) genes. At 30 days post-infection (dpi), up-regulated DEMs were enriched in Toll-like and RIG-I-like receptor signaling pathways, which were validated by qRT-PCR, Western blotting and downstream cytokines detection. Furthermore, flow cytometric analysis and serum levels of the corresponding cytokines confirmed the changes in cell-mediated immunity in host during Em infection that showed Th1 and Th17-type CD4⁺ T-cells were predominant at the early infection stage whereas Th2-type CD4⁺ T-cells were significantly higher at the middle/late stage. Collectively, our study revealed the potential regulatory functions of lncRNAs in modulating host Th cell subsets and provide novel clues in understanding the influence of Em infection on host innate and adaptive immune response.     

Analysis of competing endogenous RNA network identifies a poorly differentiated cancer-specific RNA signature for hepatocellular carcinoma

Plenty of evidence has suggested that long noncoding RNAs (lncRNAs) play a vital role in competing endogenous RNA (ceRNA) networks. Poorly differentiated hepatocellular carcinoma (PDHCC) is a malignant phenotype. This paper aimed to explore the effect and the underlying regulatory mechanism of lncRNAs on PDHCC as a kind of ceRNA. Additionally, prognosis prediction was assessed. A total of 943 messenger RNAs (mRNAs), 86 miRNAs, and 468 lncRNAs that were differentially expressed between 137 PDHCCs and 235 well-differentiated HCCs were identified. Thereafter, a ceRNA network related to the dysregulated lncRNAs was established according to bioinformatic analysis and included 29 lncRNAs, 9 miRNAs, and 96 mRNAs. RNA-related overall survival (OS) curves were determined using the Kaplan-Meier method. The lncRNA ARHGEF7-AS2 was markedly correlated with OS in HCC (P = .041). Moreover, Cox regression analysis revealed that patients with low ARHGEF7-AS2 expression were associated with notably shorter survival time (P = .038). In addition, the area under the curve values of the lncRNA signature for 1-, 3-, and 5-year survival were 0.806, 0.741, and 0.701, respectively. Furthermore, a lncRNA nomogram was established, and the C-index of the internal validation was 0.717. In vitro experiments were performed to demonstrate that silencing ARHGEF7-AS2 expression significantly promoted HCC cell proliferation and migration. Taken together, our findings shed more light on the ceRNA network related to lncRNAs in PDHCC, and ARHGEF7-AS2 may be used as an independent biomarker to predict the prognosis of HCC.     

Screening of significant biomarkers related with prognosis of liver cancer by lncRNA-associated ceRNAs analysis

This study aimed to identify significant biomarkers related to the prognosis of liver cancer using long noncoding RNA (lncRNA)-associated competing endogenous RNAs (ceRNAs) analysis. Differentially expressed mRNA and lncRNAs between liver cancer and paracancerous tissues were screened, and the functions of these mRNAs were predicted by gene ontology and pathway enrichment analyses. A ceRNA network consisting of differentially expressed mRNAs and lncRNAs was constructed. LncRNA FENDRR and lncRNA HAND2-AS1 were hub nodes in the ceRNA network. A risk score assessment model consisting of eight genes (PDE2A, ESR1, FBLN5, ALDH8A1, AKR1D1, EHHADH, ADRA1A, and GNE) associated with prognosis were developed. Multivariate Cox regression suggested that both pathologic_T and risk group could be regarded as independent prognostic factors. Furthermore, a nomogram model consisting of pathologic_T and risk group showed a good prediction ability for predicting the survival rate of liver cancer patients. The nomogram model consisting of pathologic_T and a risk score assessment model could be regarded as an independent factor for predicting prognosis of liver cancer.     

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.     

Competing endogenous RNA network and prognostic nomograms for hepatocellular carcinoma patients who underwent R0 resection

The prognosis of hepatocellular carcinoma (HCC) after R0 resection is unsatisfactory due to the high rate of recurrence. In this study, we investigated the recurrence-related RNAs and the underlying mechanism. The long noncoding RNA (lncRNA), microRNA (miRNA), and messenger RNA (mRNA) expression data and clinical information of 247 patients who underwent R0 resection patients with HCC were obtained from The Cancer Genome Atlas. Comparing the 1-year recurrence group (n = 56) with the nonrecurrence group (n = 60), we detected 34 differentially expressed lncRNAs (DElncRNAs), five DEmiRNAs, and 216 DEmRNAs. Of these, three DElncRNAs, hsa-mir-150-5p, and 11 DEmRNAs were selected for constructing the competing endogenous RNA (ceRNA) network. Next, two nomogram models were constructed based separately on the lncRNAs and mRNAs that were further selected by Cox and least absolute shrinkage and selection operator regression analysis. The two nomogram models that showed a high prediction accuracy for disease-free survival with the concordance indexes at 0.725 and 0.639. Further functional enrichment analysis of DEmRNAs showed that the mRNAs in the ceRNA network and nomogram models were associated with immune pathways. Hence, we constructed a hsa-mir-150-5p-centric ceRNA network and two effective nomogram prognostic models, and the related RNAs may be useful as potential biomarkers for predicting recurrence in patients with HCC.     

Comprehensive bioinformatics analysis identifies lncRNA HCG22 as a migration inhibitor in esophageal squamous cell carcinoma

Esophageal cancer is one of the most lethal malignancies worldwide, and esophageal squamous cell carcinoma (ESCC) is the dominant histological type. However, the long noncoding RNA (lncRNA) alterations in ESCC have not been elucidated to date. In this study, reliable databases from Gene Expression Omnibus (GEO), which analyzed lncRNA expression in ESCC tumor tissues and adjacent normal tissues were searched, and common differentially expressed lncRNAs and genes were analyzed. Next, cis- trans analysis was performed to predict the underlying relationships between altered lncRNAs and mRNAs, and the lncRNA-mRNA regulatory network was established. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of altered lncRNA-related genes were performed. The promising lncRNA HCG22 was validated by quantitative polymerase chain reaction (qPCR), and clinicopathological data were collected to identify the relationship between lncRNA HCG22 expression level and clinical features. Finally, Transwell assays were performed to explore the biological functions of lncRNA HCG22 in ESCC cells. Two hundred forty-one lncRNAs and 835 mRNAs were observed to be remarkably altered between ESCC tumor tissues and adjacent normal tissues. The lncRNA-mRNA regulatory network showed the coexpression association between lncRNA HCG22 and SPINK7 and ADAMTS12. GO and KEGG analyses showed that HCG22 and ADAMTS12 had potential biological functions in the cell migration of ESCC. The downregulation of lncRNA HCG22 in ESCC tumor tissues was validated by qPCR, and the clinicopathological data showed a noticeable correlation between lncRNA HCG22 expression level and the ESCC differentiational degree and clinical TNM stage. Kaplan-Meier analysis showed that patients with ESCC having low lncRNA HCG22 expression in ESCC tissues had considerably shorter overall survival compared with patients with ESCC having high lncRNA HCG22 expression. Following Transwell assays confirmed the migratory role of lncRNA HCG22 in ESCC cells. In conclusion, lncRNA HCG22 was downregulated in ESCC tissues and can be a migration inhibitor of ESCC cells, and SPINK7 and ADAMTS12 are promising to be the regulatory targets of lncRNA HCG22.     

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.