Next Article in Journal
Searching for New Beneficial Bacterial Isolates of Wild Raspberries for Biocontrol of Phytopathogens-Antagonistic Properties and Functional Characterization
Next Article in Special Issue
Roles and Mechanisms of the Long Noncoding RNAs in Cervical Cancer
Previous Article in Journal
The Topographical Optimization of 3D Microgroove Pattern Intervals for Ligamentous Cell Orientations: In Vitro
Previous Article in Special Issue
Long Non-Coding RNA Associated with Cholesterol Homeostasis and Its Involvement in Metabolic Diseases
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 24
10.3390/ijms21249359
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of the 3-lncRNA Signature as a Prognostic Biomarker for Colorectal Cancer

by
Shuzhen Liu
,
Qing Cao
,
Guoyan An
,
Bianbian Yan
and
Lei Lei
*
Key Laboratory of Resource Biology and Biotechnology in Western China, Ministry of Education, School of Medicine, Northwest University, Xi’an 710069, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(24), 9359; https://doi.org/10.3390/ijms21249359
Submission received: 11 November 2020 / Revised: 1 December 2020 / Accepted: 3 December 2020 / Published: 8 December 2020
(This article belongs to the Special Issue LncRNAs in Health and Disease)
Browse Figures
Versions Notes

Abstract

:
Colorectal cancer (CRC) is one of the most common malignant carcinomas in the world, and metastasis is the main cause of CRC-related death. However, the molecular network involved in CRC metastasis remains poorly understood. Long noncoding RNA (lncRNA) plays a vital role in tumorigenesis and may act as a competing endogenous RNA (ceRNA) to affect the expression of mRNA by suppressing miRNA function. In this study, we identified 628 mRNAs, 144 lncRNAs, and 25 miRNAs that are differentially expressed (DE) in metastatic CRC patients compared with nonmetastatic CRC patients from the Cancer Genome Atlas (TCGA) database. Functional enrichment analyses confirmed that the identified DE mRNAs are extensively involved in CRC tumorigenesis and migration. By bioinformatics analysis, we constructed a metastasis-associated ceRNA network for CRC that includes 28 mRNAs, 12 lncRNAs, and 15 miRNAs. We then performed multivariate Cox regression analysis on the ceRNA-related DE lncRNAs and identified a 3-lncRNA signature (LINC00114, LINC00261, and HOTAIR) with the greatest prognostic value for CRC. Clinical feature analysis and functional enrichment analysis further proved that these three lncRNAs are involved in CRC tumorigenesis. Finally, we used Transwell, Cell Counting Kit (CCK)-8, and colony formation assays to clarify that the inhibition of LINC00114 promotes the migratory, invasive, and proliferative abilities of CRC cells. The results of the luciferase assay suggest that LINC00114 is the direct target of miR-135a, which also verified the ceRNA network. In summary, this study provides a metastasis-associated ceRNA network for CRC and suggests that the 3-lncRNA signature may be a useful candidate for the diagnosis and prognosis of CRC.

1. Introduction

Worldwide, the incidence of colorectal cancer (CRC) ranks third and the mortality rate ranks second. Despite a decreasing trend of CRC mortality in many developed countries, the highest CRC incidence rates occur in Europe, Australia, and Northern America. Moreover, some developing countries have shown increases in both incidence and mortality in recent decades [1]. Metastasis is the main cause of CRC-related death, but the molecular mechanisms remain poorly understood [2]. Therefore, investigating the molecular networks that regulate metastasis is imperative for improving CRC early diagnosis and therapy.
Except for a few transcripts with protein-coding capacity, the vast majority of transcripts in humans are noncoding RNAs, which comprise two classes: small noncoding RNAs (sncRNAs) and long noncoding RNAs (lncRNAs) [3,4]. MicroRNAs (miRNAs) are sncRNAs with a length of approximately 21–24 nucleotides (nt) that target the 3′ untranslated region of mRNAs to suppress gene expression [5,6]. In contrast, lncRNAs are greater than 200 nt in length, act as translational regulators to regulate the degradation and stability of mRNAs, and play critical roles in the progression of many cancers, including CRC [7,8,9]. For instance, high expression levels of MALAT1, HOTAIR, and H19 correlate tightly with metastasis and a poor prognosis in CRC patients [10,11,12]. Salmena et al. reported that lncRNAs can act as endogenous “miRNA sponges” to inhibit miRNA function by sharing miRNA response elements (MREs) [13]. Increasing evidence indicates that molecular networks play important roles in a wide variety of human diseases [14]. Dysregulation of the competing endogenous RNA (ceRNA) network could promote tumorigenesis, including CRC [15,16].
The biological function of lncRNAs and the ceRNA network in CRC have been characterized. A recent study demonstrated that the lncRNAs KCNQ1OT1 and WT1-AS interact with most of the miRNAs in the ceRNA network and, thus, may play important roles in CRC progression [17]. Moreover, LINC00400 and LINC00355 contribute to the regulation of the ceRNA network in CRC [18]. In this study, we aim to identify a lncRNA signature as a biomarker for the diagnosis and prognosis of CRC. We first applied bioinformatics methods to construct a metastasis-associated ceRNA network. Then, multivariate Cox regression analysis was used to identify a 3-lncRNA signature according to CRC patient survival. Furthermore, we verified the function of LINC00114 using Transwell, CCK-8, and colony formation assays with CRC cell lines. This study reveals a metastasis-associated ceRNA network in CRC that might provide useful candidates for CRC pathogenesis.

2. Results

2.1. Metastasis-Specific mRNAs, miRNAs, and lncRNAs in CRC Patients

To investigate metastasis-specific mRNAs, miRNAs, and lncRNAs in CRC tissues, we used RNA expression profiles from the TCGA dataset, including 459 nonmetastatic CRC samples and 87 metastatic CRC samples (Figure 1). The results revealed 628 DE mRNAs between nonmetastatic and metastatic colorectal tumor tissues, of which 354 were upregulated and 274 were downregulated. We also obtained 144 dysregulated lncRNAs and 25 aberrantly expressed miRNAs using the “edgeR” package in R.

2.2. Functional Enrichment Analysis of DE mRNA

To identify the functions of the DE mRNAs, the 628 dysregulated mRNAs were divided into two groups (354 upregulated and 274 downregulated) for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses using DAVID and KOBAS bioinformatics resources, respectively (Figure 2). The top 15 GO and KEGG pathways indicated that the upregulated genes participate in “GO:0045815~positive regulation of gene expression, epigenetic”, “GO:0006342~chromatin silencing”, “Transcriptional misregulation in cancer”, and “Wnt signaling pathway”, which are considered to be associated with cancer metastasis [19,20,21,22]. In addition, “GO:0032200~telomere organization” and “Viral carcinogenesis” are related to oncogenesis [23]. Among the top 15 GO and KEGG pathways of downregulated mRNAs, “GO:0008544~epidermis development”, “GO:0030216~keratinocyte differentiation”, “Natural killer cell-mediated cytotoxicity”, “Chemical carcinogenesis”, and “Metabolism of xenobiotics by cytochrome P450” are related to tumorigenesis and cancer cell migration [24,25,26,27]. These findings indicate that these DE mRNAs might be associated with CRC metastasis.

2.3. Construction of the ceRNA Network

MiRNAs play a major role in ceRNA networks by suppressing the expression of target mRNAs. Salmena et al. reported that lncRNAs could act as miRNA sponges by occupying their binding sites and inhibiting their functions [13]. In this study, we used miRDB and Targetscan to predict target mRNAs of 25 DE miRNAs. We then compared the predicted mRNAs with 628 DE mRNAs and selected 28 mRNAs to construct the ceRNA network (Table S1). We also studied whether these 25 miRNAs target 144 DE lncRNAs. According to the miRcode results, 9 miRNAs are predicted to interact with 12 lncRNAs (Table S2).
Based on the above evidence (Tables S1 and S2), we next constructed a ceRNA network using Cytoscape v3.7. Ultimately, 28 mRNAs, 15 miRNAs (Table S3), and 12 lncRNAs (Table S4) are involved in the ceRNA network. Figure 3 shows that 4 downregulated miRNAs, 10 upregulated mRNAs, and 5 upregulated lncRNAs are involved in one ceRNA network. Meanwhile, 11 upregulated miRNAs, 18 downregulated mRNAs, and 7 downregulated lncRNAs are involved in another ceRNA network.

2.4. Construction of the 3-lncRNA Prognostic Model

The 12 lncRNAs in the ceRNA network were further analyzed according to CRC patient survival. We randomly divided CRC patients into a training group (n = 272) and a testing group (n = 272), and no significant differences were found between the two groups (Table S5). By multivariate Cox regression analysis, three lncRNAs (LINC00114, LINC00261, and HOTAIR) displayed a remarkable prognostic value in the training set (Table S6). HOTAIR was associated with high risk, whereas LINC00114 and LINC00261 were protective (Figure 4A). CRC patients were classified into low-risk (n = 136) and high-risk (n = 136) sets based on the risk score formula (risk score = 0.1184 × HOTAIR − 0.2257 × LINC00114 − 0.1477 × LINC00261), and Kaplan–Meier analysis indicated a positive correlation with overall survival (OS) for the low-risk group compared with the high-risk group (p = 0.002; Figure 4B). Additionally, the accuracy of the prognostic model was evaluated by receiver operating characteristic (ROC) curve analysis, and the area under the ROC curve (AUC) at 5 years for OS was 0.716 (Figure 4C).
In order to verify the predictive ability of the 3-lncRNA model, Kaplan–Meier analysis and ROC curve analyses were further used on the testing group. According to the previous risk score formula, patients in the testing group were divided into low-risk and high-risk groups (Figure 4D). There were statistically significant differences between the low-risk group and the high-risk group in both Kaplan–Meier analysis (p = 0.02; Figure 4E) and ROC curve analysis (AUC = 0.649; Figure 4F). These findings demonstrate that the 3-lncRNA signature may act as an independent biomarker that predicts CRC survival prognosis.

2.5. Clinical Feature Analysis of the Key lncRNAs

To further research the three lncRNAs, we analyzed the association between the expression of lncRNAs and clinical features of CRC. Kaplan–Meier survival curves indicated that HOTAIR was negatively correlated with OS (p = 0.012), while LINC00261 and LINC00114 were positively correlated with OS (p = 0.044, p = 0.0006, respectively; Figure 5A). Although there was no significant difference in the expression of HOTAIR in normal tissues and tumors (Figure 5B), HOTAIR was upregulated in metastatic CRC tissues and stage 4 tissues compared with nonmetastatic samples and stage 1 tissues, respectively (Figure 5C,D). In contrast, the expression of LINC00261 showed no significance between normal tissues and tumors and was downregulated in metastatic CRC tissues and stage 4 tissues (Figure 5B–D). The expression of LINC00114 was downregulated in tumor tissues, metastatic CRC tissues, and stage 4 tissues (Figure 5B–D).

2.6. Functional Enrichment Analysis of the Key lncRNAs

We further investigated the function of the three key lncRNAs according to the principle that coexpressed genes might share roles in related biological processes [28,29]. First, we searched for mRNAs that are coexpressed with the key lncRNAs based on the web resource Multi-Experiment Matrix (MEM) [30], which searches a large number of microarray datasets and uses rank aggregation to detect globally coexpressed genes. The top 200 mRNAs that are potentially coexpressed with HOTAIR and LINC00261 are shown in Table S7, while LINC00114 could not be searched for on MEM.
We then used GO and KEGG pathway enrichment analyses to assess the enrichment of these mRNAs. The top 15 biological processes in the GO and KEGG pathways showed that HOTAIR participates in “GO:0016477~cell migration”, “Hippo signaling pathway”, “Ubiquitin-mediated proteolysis”, and “TGF-beta signaling pathway” (Figure 6A,B), which are correlated with tumorigenesis and cancer metastasis [31,32,33]. The mRNAs positively coexpressed with LINC00261 are mainly involved in “GO:0042158~lipoprotein biosynthetic process”, “GO:0008152~metabolic process”, “Drug metabolism—cytochrome P450”, and “Chemical carcinogenesis” (Figure 6C,D), suggesting that LINC00261 may regulate tumorigenesis [24,34] and also coincide with functional enrichment of metastasis-associated downregulated mRNAs (Figure 2D).
We performed GSEA to further determine the pathways in which LINC00114 affects CRC. The gene sets “KRAS.300_UP.V1_UP” (p < 0.0001), “TGF_BETA_SIGNALING” (p = 0.002), and “PTEN_DN.V1_UP” (p < 0.0001) were remarkably enriched with low levels of LINC00114, while “GNF2_MSH6” (p = 0.018) was enriched with high levels of LINC00114 (Figure 6E–H). Many studies have shown that both activation of KRAS [35] and TGF-β pathways [36] and loss of PTEN [37] and MSH6 [38] are the driving force of tumorigenesis [39]. Therefore, these results suggest that LINC00114 might be involved in CRC progression as a tumor suppressor.

2.7. LINC00114 Represses CRC Cell Migration and Proliferation

To further determine whether LINC00114 impacts CRC tumorigenesis, we validated the effect of LINC00114 knockdown on the migration and proliferation of CRC cells. After we confirmed the suppression efficiency of the three LINC00114 siRNAs in HCT-116 and LOVO cells by quantitative RT-PCR, siLINC00114-3 and siLINC00114-1 were selected for subsequent functional experiments, respectively (Figure 7A,B). Transwell assays showed that compared to the negative control, suppression of LINC00114 by siRNA promoted both migration and invasion abilities (Figure 7C,D). CCK-8 and colony formation assays were used to assess the role of LINC00114 in CRC proliferation. Compared to the negative control, inhibition of LINC00114 significantly promoted CRC cell proliferation (Figure 8A–D). These results indicate that LINC00114 represses the migratory, invasive, and proliferative abilities of CRC cells.
The ceRNA networks revealed that LINC00114 may target hsa-mir-34b and hsa-mir-135a. By Pearson correlation analysis, we found that LINC00114 was negatively associated with hsa-mir-135a (Figure 9A). This result was also confirmed by transfecting miR-135a mimics in LOVO cells, which indicated that overexpression of miR-135a suppressed the expression of LINC00114 (Figure 9B). To figure out whether LINC00114 is a direct target of miR-135a, we constructed plasmids with a wild-type fragment of LINC00114 (wt-LINC00114) and a mutant fragment of LINC00114 (mt-LINC00114). Luciferase assays showed that coexpression of miR-135a mimics luciferase activity and wt-LINC00114 markedly reduced luciferase activity compared with the control (Figure 9C). These data suggest that LINC00114 is the direct target of miR-135a, which also verifies the ceRNA network.

3. Discussion

CRC is one of the most common tumors worldwide, and it was the second leading cause of mortality in 2018 [1]. Metastasis is the main reason for the poor overall 5-year survival rate of advanced cases [2]. Recently, an increasing number of articles have shown that dysregulated lncRNA expression is related to CRC tumorigenesis and metastasis [10,11,12]. The ceRNA network hypothesis suggests that lncRNAs are able to regulate mRNA expression as endogenous “miRNA sponges” by binding to miRNA MREs [13]. However, the role of ceRNAs in CRC metastasis has remained unclear. The present study reveals the metastasis-associated ceRNA network in CRC, which provides new biomarkers for predicting CRC survival.
In this study, we identified DE mRNAs, lncRNAs, and miRNAs between metastatic CRC patients and nonmetastatic CRC patients. We used bioinformatics methods to construct a metastasis-associated ceRNA network that includes 28 mRNAs, 15 miRNAs, and 12 lncRNAs. By using the MIENTURNET web tool, we found that most of the DE miRNAs are remarkably enriched in pathways associated with tumorigenesis [40]. The results of GO and KEGG pathway analyses indicated that upregulated mRNAs mainly participate in “Transcriptional misregulation in cancer” and “Wnt signaling pathway”. Many studies have shown that invasive edge-localized CRC cells exhibit strong nuclear β-catenin staining [41] and preferentially express Wnt target genes, which confer invasive metastatic capacity [42,43]. In contrast, downregulated mRNAs were enriched in “GO:0008544~epidermis development” and “GO:0030216~keratinocyte differentiation”, in which small proline-rich protein 2F (SPRR2F), small proline-rich protein 1B (SPRR1B), small proline-rich protein 2G (SPRR2G), late cornified envelope 3D (LCE3D), and late cornified envelope 1E (LCE1E) are involved. A previous study showed that overexpression of SPRR proteins increases local aggressiveness but prevents metastasis in cholangiocarcinoma [44]. Therefore, the downregulation of these genes promotes the metastatic capacity of CRC. These data indicate that the metastasis-related DE mRNAs are involved in CRC tumorigenesis and migration.
LncRNAs play critical roles in the progression of cancers and may act as candidate indicators to predict survival prognosis. A previous study showed that a 4-lncRNA signature (LINC01555, RP11-610P16.1, RP11-108K3.1, and LINC01207) had the most significant correlation with OS in colon adenocarcinoma patients [45]. In this study, we determined dysregulated lncRNA profiles in metastatic CRC and constructed metastasis-associated ceRNA networks in large CRC cohorts. We further found that 3 lncRNAs (LINC00114, LINC00261, and HOTAIR) showed a significant prognostic value for CRC based on multivariate Cox regression analysis. HOTAIR has been reported to increase in primary tumors and metastases, and its expression level is a powerful predictor of metastasis and mortality [10,20,46,47,48]. Additionally, a recent study showed that HOTAIR controls Rab22a expression in ovarian cancer by sponging miR-373 [49] and regulates HER2 expression in gastric cancer by sponging miR-331-3p [50]. In our study, HOTAIR was found to be negatively associated with OS and was upregulated in metastatic CRC and stage 4 disease. Moreover, GO and KEGG pathway analyses showed that the mRNAs that are positively coexpressed with HOTAIR participate in “GO:0016477~cell migration”, “Transcriptional misregulation in cancer”, “Hippo signalling pathway”, and “TGF-beta signaling pathway”, which are considered to be correlated with CRC tumorigenesis and metastasis [31,32,33].
LINC00261 also plays a notable role in CRC pathogenesis and metastasis. Recent research has demonstrated that overexpression of LINC00261 relieves cisplatin resistance in CRC cells by promoting apoptosis and inhibiting migration [51]. Several studies have shown that LINC00261 exerts tumor-suppressive effects in gastric cancer and pancreatic cancer [52,53]. Moreover, LINC00261 competes with miR-132-3p to regulate the expression of BCL2L11 in endometriosis [54]. In the present study, LINC00261 was positively correlated with OS and downregulated in metastatic CRC tissues and stage 4 disease. Interestingly, LINC00261 coexpressed genes were mainly clustered in “GO:0008152~metabolic process”, “Drug metabolism—cytochrome P450”, and “Chemical carcinogenesis”, which also coincided with the functional enrichment of metastasis-associated downregulated mRNAs. These results indicate that overexpression of LINC00261 may inhibit CRC tumorigenesis and metastasis.
We used Transwell, CCK-8, and colony formation assays to verify the function of LINC00114 in CRC cell lines. The results indicated that suppression of LINC00114 promoted migratory, invasive, and proliferative abilities compared to the negative control, which was consistent with in-silico analysis. However, a recent study showed that LINC00114 depletion inhibited CRC cell proliferation [55]. As they did not investigate the role of LINC00114 in the migration of CRC cells, we will further study the functional mechanism of LINC00114 in the future.
Based on the ceRNA network, LINC00114 may compete with two key DE miRNAs (miR-135a and miR-34b) to regulate the expression of target genes. Several studies have confirmed that miR-135a and miR-34b are highly expressed in CRC tissues and promote CRC cell proliferation, migration, and invasion [56,57,58]. In this article, we showed that LINC00114 is the direct target of miR-135a by luciferase assay.
In conclusion, we constructed a metastasis-associated ceRNA network in CRC and proved that a 3-lncRNA signature that includes LINC00114, LINC00261, and HOTAIR is an independent factor for predicting CRC prognosis. The current study provides insight into the identification of potential biomarkers for the diagnosis and prognosis of CRC.

4. Materials and Methods

4.1. Study Population and Identification of Differential Expression

CRC transcriptome profiles (level 3) and corresponding clinical information were downloaded. These data included 51 normal tissues and 649 CRC tissues. According to the TNM staging systems, the number of M1 (metastatic cases) and M0 (nonmetastatic cases) tumor tissues were 87 and 459, respectively (unknown information for 103 tumor tissues).
We divided the CRC samples into 2 groups (M0 and M1) and used the “edgeR” package in R software to analyze DE RNAs (mRNA, miRNA, and lncRNA) [59]. Compared to the M0 group, we obtained 628 DE mRNAs, 144 DE lncRNAs, and 25 DE miRNAs, with thresholds of |logFC| > 1 and FDR (false discovery rate) < 0.05.

4.2. Functional Enrichment Analysis

To explore the functions of the identified DE mRNAs, 628 DE mRNAs were divided into an upregulated group and a downregulated group for Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. DAVID was used to analyze GO biological processes, and KOBAS was used to analyze KEGG pathways. The top 15 pathways were visualized using the R package ggplot2.
To clarify the function of the DE lncRNAs, we first searched for coexpressed mRNAs based on the web resource Multi-Experiment Matrix (MEM). The top 200 potentially coexpressed mRNAs for HOTAIR and LINC00261 were selected for GO and KEGG analyses.
Gene set enrichment analysis (GSEA) was performed to investigate the function of LINC00114. CRC gene sets were obtained from the TCGA database and uploaded to GSEA software (v3.0). The phenotype label was high-level of LINC00114 versus low-level of LINC00114. The permutation number was 1000, and a ranked-list metric was created by counting the signal-to-noise ratio.

4.3. Construction of the ceRNA Network

We used metastasis-associated DE mRNAs, lncRNAs, and miRNAs to construct the ceRNA network. miRDB and Targetscan were used to predict mRNA-miRNA interactions. In addition, miRcode was performed to construct lncRNA–miRNA interactions. We retained intersections with the DE mRNAs, miRNAs, and lncRNAs to build the ceRNA network, as the flow chart in Figure 1 illustrates. Cytoscape v3.7.1 was used to visualize the ceRNA network.

4.4. Definition of the Prognostic Model and the Survival Curve

The 544 CRC patients were randomly divided into a training group (n = 272) and a testing group (n = 272). Multivariate Cox regression analysis was performed to identify the prognostic model for the lncRNAs in the ceRNA network. Based on the Akaike information criterion (AIC), the best mathematical model was established using the “survival” R package [60]. A risk score for predicting overall survival (OS) was generated by including the expression level of each lncRNA (explncRNA) weighted by the regression coefficient (βlncRNA), as follows:
Risk score = βlncRNA1 × explncRNA1 + βlncRNA2 × explncRNA2 + βlncRNAn × explncRNAn.
According to the risk score of the lncRNA prognostic model, patients were divided into a low-risk group and a high-risk group. Kaplan–Meier analysis was utilized to generate the OS curve for the two groups. The accuracy of the prognostic model within 5 years was evaluated by time-dependent receiver operating characteristic (ROC) curve analysis using the “survivalROC” R package [61].

4.5. Cell Lines and Transfection

HCT-116 and LOVO human CRC cells were purchased from American Type Culture Collection (ATCC). A short interfering RNA (siRNA) for LINC00114 and a negative control were transfected at 100 nM according to the instructions for the X-tremeGENE siRNA transfection reagent (Roche, Mannheim, Germany). The following siRNA sequences were used in this study: siLINC00114-1, sense, 5′-GGAGCCGACUUGCUUAGAATT-3′; antisense: 5′-UUCUAAGCAAGUCGGCUCCTT-3′; siLINC00114-2, sense, 5′-GCGCAUCUCCUAAGAACUUTT-3′, antisense, 5′-AAGUUCUUAGGAGAUGCGCTT-3′; siLINC00114-3, sense, 5′-GCAUGAAGACUAAAUCCAUTT-3′, antisense: 5′-AUGGAUUUAGUCUUCAUGCTT-3′.

4.6. RNA Extraction and Quantitative RT-PCR

Total RNA was extracted from transfected HCT-116 cells using TRIzol (Invitrogen, Carlsbad, CA, USA). After cDNA synthesis (TAKARA, Dalian, China), quantitative real-time PCR assays were performed using the SYBR Premix Ex Taq™ II kit (TAKARA, Dalian, China) with a CFX96 Touch Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). The following primer sequences were used in this study: LINC00114, forward, 5′-ACTTGTGGCATCATAGCAG-3′, reverse, 5′-AAGTGAAACAGTTCCCTTCT-3′; GAPDH, forward, 5′-TGGAAATCCCATCACCATCT-3′, reverse, 5′-TGGACTCCACGACGTACTCA-3′.

4.7. Transwell Assays

Transwell assays (Corning, MA, USA), without or with Matrigel (BD Biosciences, Bedford, MA, USA), were used to determine cell migration or invasion ability, respectively. Briefly, 5 × 104 HCT-116 or LOVO cells with serum-free medium were placed into the upper chamber, with serum-containing medium in the bottom chamber. After 48 h, the migrated cells were fixed in 4% paraformaldehyde and then stained with crystal violet.

4.8. Cell Proliferation

After 48 h of transfection, HCT-116 or LOVO cells were seeded at a density of 1000 cells per well into 96-well plates, with 6 replicates. Briefly, 10% of Cell Counting Kit (CCK)-8 solution (Dojindo, Kyushu, Japan) was incubated at 37 °C for 3 h and measured at a wavelength of 450 nm.
For colony formation assays, we seeded 1000 transfected HCT-116 or LOVO cells per well in 6-well plates. Fourteen days later, the number of colonies was counted. Each assay was performed in triplicate.

4.9. Dual-Luciferase Assay

The LINC00114 wild-type fragment (wt-LINC00114) and mutant fragment (mt-LINC00114) were constructed in the pmiRGLO plasmid (Promega, Madison, WI, USA). For dual-luciferase assay, the method has been described previously [62].

4.10. Statistical Analysis

GraphPad Prism 6 (GraphPad Software Inc., La Jolla, CA, USA) was used for nonpaired t-tests and Pearson’s correlation analyses. R Studio (R version 3.3.2) was used for bioinformatics analyses.

Supplementary Materials

Supplementary Materials can be found at https://www.mdpi.com/1422-0067/21/24/9359/s1. Table S1: 13 DE miRNAs target with 28 DE mRNAs in metastasis CRC; Table S2: 9 DE miRNAs interact with 12 DE lncRNAs in metastasis CRC; Table S3: Metastasis CRC specific DE miRNAs in ceRNA network; Table S4: Metastasis CRC specific DE lncRNAs in ceRNA network; Table S5: Clinical covariates in the training and testing sets; Table S6: 3-lncRNA risk score model; Table S7: The top 200 mRNAs co-expressed with LINC00261 and HOTAIR.

Author Contributions

S.L. and L.L. contributed to the conception of the project and the design of all experiments. S.L., Q.C., G.A., and B.Y. carried out all bioinformatic analyses. Experiments were conducted by S.L. and Q.C. L.L. wrote the main body of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the National Natural Science Foundation of China (grant no. 31401081), the Shaanxi Province Natural Science Basic Research Program (grant no. 2019JQ-538), and the Opening Foundation of Key Laboratory of Resource Biology and Biotechnology in Western China (Northwest University), Ministry of Education (No.ZSK2018002).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [Green Version]
  2. Tauriello, D.V.F.; Calon, A.; Lonardo, E.; Batlle, E. Determinants of metastatic competency in colorectal cancer. Mol. Oncol. 2017, 11, 97–119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Consortium, E.P. Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature 2007, 447, 799–816. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Fang, Y.; Fullwood, M.J. Roles, Functions, and Mechanisms of Long Non-coding RNAs in Cancer. Genom. Proteom. Bioinform. 2016, 14, 42–54. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Lee, R.C.; Feinbaum, R.L.; Ambros, V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell 1993, 75, 843–854. [Google Scholar] [CrossRef]
  6. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  7. Huarte, M. The emerging role of lncRNAs in cancer. Nat. Med. 2015, 21, 1253. [Google Scholar] [CrossRef] [PubMed]
  8. Han, D.; Wang, M.; Ma, N.; Xu, Y.; Jiang, Y.; Gao, X. Long noncoding RNAs: Novel players in colorectal cancer. Cancer Lett. 2015, 361, 13–21. [Google Scholar] [CrossRef]
  9. Poliseno, L.; Salmena, L.; Zhang, J.; Carver, B.; Haveman, W.J.; Pandolfi, P.P. A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature 2010, 465, 1033–1038. [Google Scholar] [CrossRef] [Green Version]
  10. Kogo, R.; Shimamura, T.; Mimori, K.; Kawahara, K.; Imoto, S.; Sudo, T.; Tanaka, F.; Shibata, K.; Suzuki, A.; Komune, S.; et al. Long Noncoding RNA HOTAIR Regulates Polycomb-Dependent Chromatin Modification and Is Associated with Poor Prognosis in Colorectal Cancers. Cancer Res. 2011, 71, 6320–6326. [Google Scholar] [CrossRef] [Green Version]
  11. Liang, W.-C.; Fu, W.-M.; Wong, C.-W.; Wang, Y.; Wang, W.-M.; Hu, G.-X.; Zhang, L.; Xiao, L.-J.; Wan, D.C.-C.; Zhang, J.-F.; et al. The lncRNA H19 promotes epithelial to mesenchymal transition by functioning as miRNA sponges in colorectal cancer. Oncotarget 2015, 6, 22513–22525. [Google Scholar] [CrossRef] [Green Version]
  12. Zheng, H.-T.; Shi, D.-B.; Wang, Y.-W.; Li, X.-X.; Xu, Y.; Tripathi, P.; Gu, W.-L.; Cai, G.; Cai, S.-J. High expression of lncRNA MALAT1 suggests a biomarker of poor prognosis in colorectal cancer. Int. J. Clin. Exp. Pathol. 2014, 7, 3174–3181. [Google Scholar] [PubMed]
  13. Salmena, L.; Poliseno, L.; Tay, Y.; Kats, L.; Pandolfi, P.P. A ceRNA hypothesis: The Rosetta Stone of a hidden RNA language? Cell 2011, 146, 353–358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Silverman, E.K.; Schmidt, H.H.H.W.; Anastasiadou, E.; Altucci, L.; Angelini, M.; Badimon, L.; Balligand, J.-L.; Benincasa, G.; Capasso, G.; Conte, F.; et al. Molecular networks in Network Medicine: Development and applications. Wiley Interdiscip. Rev. Syst. Biol. Med. 2020, 12, e1489. [Google Scholar] [CrossRef] [PubMed]
  15. Conte, F.; Fiscon, G.; Chiara, M.; Colombo, T.; Farina, L.; Paci, P. Role of the long non-coding RNA PVT1 in the dysregulation of the ceRNA-ceRNA network in human breast cancer. PLoS ONE 2017, 12, e0171661. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Russo, F.; Fiscon, G.; Conte, F.; Rizzo, M.; Paci, P.; Pellegrini, M. Interplay between long noncoding RNAs and MicroRNAs in cancer. In Computational Cell Biology; Humana Press: New York, NY, USA, 2018; pp. 75–92. [Google Scholar]
  17. Li, F.; Li, Q.; Wu, X. Construction and analysis for differentially expressed long non-coding RNAs and MicroRNAs mediated competing endogenous RNA network in colon cancer. PLoS ONE 2018, 13, e0192494. [Google Scholar] [CrossRef] [PubMed]
  18. Yuan, W.; Li, X.; Liu, L.; Wei, C.; Sun, D.; Peng, S.; Jiang, L. Comprehensive analysis of lncRNA-associated ceRNA network in colorectal cancer. Biochem. Biophys. Res. Commun. 2019, 508, 374–379. [Google Scholar] [CrossRef]
  19. Wang, Y.; Shang, Y. Epigenetic control of epithelial-to-mesenchymal transition and cancer metastasis. Exp. Cell Res. 2013, 319, 160–169. [Google Scholar] [CrossRef]
  20. Gupta, R.A.; Shah, N.; Wang, K.C.; Kim, J.; Horlings, H.M.; Wong, D.J.; Tsai, M.-C.; Hung, T.; Argani, P.; Rinn, J.L. Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature 2010, 464, 1071. [Google Scholar] [CrossRef]
  21. Ell, B.; Kang, Y. Transcriptional control of cancer metastasis. Trends Cell Biol. 2013, 23, 603–611. [Google Scholar] [CrossRef] [Green Version]
  22. Yu, M.; Ting, D.T.; Stott, S.L.; Wittner, B.S.; Ozsolak, F.; Paul, S.; Ciciliano, J.C.; Smas, M.E.; Winokur, D.; Gilman, A.J.; et al. RNA sequencing of pancreatic circulating tumour cells implicates WNT signalling in metastasis. Nature 2012, 487, 510–513. [Google Scholar] [CrossRef]
  23. Jafri, M.A.; Ansari, S.A.; Alqahtani, M.H.; Shay, J.W. Roles of telomeres and telomerase in cancer, and advances in telomerase-targeted therapies. Genome Med. 2016, 8, 69. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kumarakulasingham, M.; Rooney, P.H.; Dundas, S.R.; Telfer, C.; Melvin, W.T.; Curran, S.; Murray, G.I. Cytochrome p450 profile of colorectal cancer: Identification of markers of prognosis. Clin. Cancer Res. 2005, 11, 3758–3765. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Koizumi, K.; Hojo, S.; Akashi, T.; Yasumoto, K.; Saiki, I. Chemokine receptors in cancer metastasis and cancer cell-derived chemokines in host immune response. Cancer Sci. 2007, 98, 1652–1658. [Google Scholar] [CrossRef] [PubMed]
  26. Hanna, N. Expression of metastatic potential of tumor cells in young nude mice is correlated with low levels of natural killer cell-mediated cytotoxicity. Int. J. Cancer 1980, 26, 675–680. [Google Scholar] [CrossRef]
  27. Golan, T.; Messer, A.R.; Amitai-Lange, A.; Melamed, Z.; Ohana, R.; Bell, R.E.; Kapitansky, O.; Lerman, G.; Greenberger, S.; Khaled, M.; et al. Interactions of Melanoma Cells with Distal Keratinocytes Trigger Metastasis via Notch Signaling Inhibition of MITF. Mol. Cell 2015, 59, 664–676. [Google Scholar] [CrossRef] [Green Version]
  28. Wolfe, C.J.; Kohane, I.S.; Butte, A.J. Systematic survey reveals general applicability of “guilt-by-association” within gene coexpression networks. BMC Bioinform. 2005, 6, 227. [Google Scholar] [CrossRef] [Green Version]
  29. Stuart, J.M.; Segal, E.; Koller, D.; Kim, S.K. A gene-coexpression network for global discovery of conserved genetic modules. Science 2003, 302, 249–255. [Google Scholar] [CrossRef] [Green Version]
  30. Adler, P.; Kolde, R.; Kull, M.; Tkachenko, A.; Peterson, H.; Reimand, J.; Vilo, J. Mining for coexpression across hundreds of datasets using novel rank aggregation and visualization methods. Genome Biol. 2009, 10, R139. [Google Scholar] [CrossRef] [Green Version]
  31. Pan, D. The hippo signaling pathway in development and cancer. Dev. Cell 2010, 19, 491–505. [Google Scholar] [CrossRef] [Green Version]
  32. Kang, T.; Wei, Y.; Honaker, Y.; Yamaguchi, H.; Appella, E.; Hung, M.-C.; Piwnica-Worms, H. GSK-3β targets Cdc25A for ubiquitin-mediated proteolysis, and GSK-3β inactivation correlates with Cdc25A overproduction in human cancers. Cancer Cell 2008, 13, 36–47. [Google Scholar] [CrossRef] [Green Version]
  33. Wakefield, L.M.; Roberts, A.B. TGF-β signaling: Positive and negative effects on tumorigenesis. Curr. Opin. Genet. Dev. 2002, 12, 22–29. [Google Scholar] [CrossRef]
  34. Santos, C.R.; Schulze, A. Lipid metabolism in cancer. FEBS J. 2012, 279, 2610–2623. [Google Scholar] [CrossRef] [PubMed]
  35. Pylayeva-Gupta, Y.; Grabocka, E.; Bar-Sagi, D. RAS oncogenes: Weaving a tumorigenic web. Nat. Rev. Cancer 2011, 11, 761–774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Bierie, B.; Moses, H.L. TGF-β and cancer. Cytokine Growth Factor Rev. 2006, 17, 29–40. [Google Scholar] [CrossRef] [PubMed]
  37. Sartore-Bianchi, A.; Di Nicolantonio, F.; Nichelatti, M.; Molinari, F.; De Dosso, S.; Saletti, P.; Martini, M.; Cipani, T.; Marrapese, G.; Mazzucchelli, L.; et al. Multi-Determinants Analysis of Molecular Alterations for Predicting Clinical Benefit to EGFR-Targeted Monoclonal Antibodies in Colorectal Cancer. PLoS ONE 2009, 4, e7287. [Google Scholar] [CrossRef] [PubMed]
  38. Baglietto, L.; Lindor, N.M.; Dowty, J.G.; White, D.M.; Wagner, A.; Garcia, E.B.G.; Vriends, A.H.J.T.; Cartwright, N.R.; Barnetson, R.A.; Farrington, S.M.; et al. Risks of Lynch Syndrome Cancers for MSH6 Mutation Carriers. J. Natl. Cancer Inst. 2010, 102, 193–201. [Google Scholar] [CrossRef] [PubMed]
  39. Armaghany, T.; Wilson, J.D.; Chu, Q.; Mills, G. Genetic alterations in colorectal cancer. Gastrointes. Cancer Res. 2012, 5, 19–27. [Google Scholar]
  40. Licursi, V.; Conte, F.; Fiscon, G.; Paci, P. MIENTURNET: An interactive web tool for microRNA-target enrichment and network-based analysis. BMC Bioinform. 2019, 20, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Brabletz, T.; Jung, A.; Hermann, K.; Günther, K.; Hohenberger, W.; Kirchner, T. Nuclear Overexpression of the Oncoprotein β-Catenin in Colorectal Cancer is Localized Predominantly at the Invasion Front. Pathol. Res. Pract. 1998, 194, 701–704. [Google Scholar] [CrossRef]
  42. Gavert, N.; Conacci-Sorrell, M.; Gast, D.; Schneider, A.; Altevogt, P.; Brabletz, T.; Ben-Ze’Ev, A. L1, a novel target of β-catenin signaling, transforms cells and is expressed at the invasive front of colon cancers. J. Cell Biol. 2005, 168, 633–642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Vignjevic, D.M.; Schoumacher, M.; Gavert, N.; Janssen, K.-P.; Jih, G.; Laé, M.; Louvard, D.; Ben-Ze’Ev, A.; Robine, S. Fascin, a Novel Target of β-Catenin-TCF Signaling, Is Expressed at the Invasive Front of Human Colon Cancer. Cancer Res. 2007, 67, 6844–6853. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Specht, S.; Isse, K.; Nozaki, I.; Lunz, J.G.; Demetris, A.J. SPRR2A expression in cholangiocarcinoma increases local tumor invasiveness but prevents metastasis. Clin. Exp. Metastasis 2013, 30, 877–890. [Google Scholar] [CrossRef] [PubMed]
  45. Zeng, J.-H.; Liang, L.; He, R.-Q.; Tang, R.-X.; Cai, X.-Y.; Chen, J.-Q.; Luo, D.-Z.; Chen, G. Comprehensive investigation of a novel differentially expressed lncRNA expression profile signature to assess the survival of patients with colorectal adenocarcinoma. Oncotarget 2017, 8, 16811–16828. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Yang, Z.; Zhou, L.; Wu, L.-M.; Lai, M.-C.; Xie, H.-Y.; Zhang, F.; Zheng, S. Overexpression of Long Non-coding RNA HOTAIR Predicts Tumor Recurrence in Hepatocellular Carcinoma Patients Following Liver Transplantation. Ann. Surg. Oncol. 2011, 18, 1243–1250. [Google Scholar] [CrossRef]
  47. Hajjari, M.; Salavaty, A. HOTAIR: An oncogenic long non-coding RNA in different cancers. Cancer Biol. Med. 2015, 12, 1–9. [Google Scholar]
  48. Teschendorff, A.E.; Lee, S.-H.; Jones, A.; Fiegl, H.; Kalwa, M.; Wagner, W.; Chindera, K.; Evans, I.; Dubeau, L.; Orjalo, A.; et al. HOTAIR and its surrogate DNA methylation signature indicate carboplatin resistance in ovarian cancer. Genome Med. 2015, 7, 108. [Google Scholar] [CrossRef] [Green Version]
  49. Zhang, Z.; Cheng, J.; Wu, Y.; Qiu, J.; Sun, Y.; Tong, X. LncRNA HOTAIR controls the expression of Rab22a by sponging miR-373 in ovarian cancer. Mol. Med. Rep. 2016, 14, 2465–2472. [Google Scholar] [CrossRef] [Green Version]
  50. Liu, X.-H.; Sun, M.; Nie, F.; Ge, Y.-B.; Zhang, E.-B.; Yin, D.-D.; Kong, R.; Chen, J.; Lu, K.-H.; Li, J.-H.; et al. Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-331-3p in gastric cancer. Mol. Cancer 2014, 13, 92. [Google Scholar] [CrossRef] [Green Version]
  51. Wang, Z.; Yang, L.; Wu, L.; Mao, H.; Zhou, Y.; Zhang, P.; Dai, G. Long non-coding RNA LINC00261 sensitizes human colon cancer cells to cisplatin therapy. Braz. J. Med Biol. Res. 2018, 51, e6793. [Google Scholar] [CrossRef] [Green Version]
  52. Yu, Y.; Li, L.; Zheng, Z.; Chen, S.; Chen, E.; Hu, Y. Long non-coding RNA linc00261 suppresses gastric cancer progressionviapromoting Slug degradation. J. Cell. Mol. Med. 2017, 21, 955–967. [Google Scholar] [CrossRef]
  53. Zhang, B.; Li, C.; Sun, Z. Long non-coding RNA LINC00346, LINC00578, LINC00673, LINC00671, LINC00261, and SNHG9 are novel prognostic markers for pancreatic cancer. Am. J. Transl. Res. 2018, 10, 2648–2658. [Google Scholar] [PubMed]
  54. Wang, H.; Sha, L.; Huang, L.; Yang, S.; Zhou, Q.; Luo, X.; Shi, B. LINC00261 functions as a competing endogenous RNA to regulate BCL2L11 expression by sponging miR-132-3p in endometriosis. Am. J. Transl. Res. 2019, 11, 2269–2279. [Google Scholar] [PubMed]
  55. Lv, L.; He, L.; Chen, S.; Yu, Y.; Che, G.; Tao, X.; Wang, S.; Jian, Z.; Zhang, X. Long Non-coding RNA LINC00114 Facilitates Colorectal Cancer Development Through EZH2/DNMT1-Induced miR-133b Suppression. Front. Oncol. 2019, 9, 1383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Nagel, R.; Le Sage, C.; Diosdado, B.; Van Der Waal, M.; Vrielink, J.A.O.; Bolijn, A.; Meijer, G.A.; Agami, R. Regulation of the Adenomatous Polyposis Coli Gene by the miR-135 Family in Colorectal Cancer. Cancer Res. 2008, 68, 5795–5802. [Google Scholar] [CrossRef] [Green Version]
  57. Zhou, W.; Li, X.; Liu, F.; Xiao, Z.; He, M.; Shen, S.-R.; Liu, S. MiR-135a promotes growth and invasion of colorectal cancer via metastasis suppressor 1 in vitro. Acta Biochim. Biophys. Sin. 2012, 44, 838–846. [Google Scholar] [CrossRef] [Green Version]
  58. Hiyoshi, Y.; Schetter, A.J.; Okayama, H.; Inamura, K.; Anami, K.; Nguyen, G.H.; Horikawa, I.; Hawkes, J.E.; Bowman, E.D.; Leung, S.-Y.; et al. Increased MicroRNA-34b and -34c Predominantly Expressed in Stromal Tissues Is Associated with Poor Prognosis in Human Colon Cancer. PLoS ONE 2015, 10, e0124899. [Google Scholar] [CrossRef] [Green Version]
  59. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  60. Therneau, T.; Lumley, T. Survival: Survival Analysis, Including Penalised Likelihood; R package Version 2.35-7; R Foundation for Statistical Computing: Vienna, Austria, 2009. [Google Scholar]
  61. Heagerty, P.J.; Lumley, T.; Pepe, M.S. Time-dependent ROC curves for censored survival data and a diagnostic marker. Biometrics 2000, 56, 337–344. [Google Scholar] [CrossRef]
  62. Lei, L.; Zhao, X.; Liu, S.; Cao, Q.; Yan, B.; Yang, J. MicroRNA-3607 inhibits the tumorigenesis of colorectal cancer by targeting DDI2 and regulating the DNA damage repair pathway. Apoptosis 2019, 24, 662–672. [Google Scholar] [CrossRef]
Ijms 21 09359 g001 550
Figure 1. Flow chart of construction of competing endogenous RNA (ceRNA) network.
Figure 1. Flow chart of construction of competing endogenous RNA (ceRNA) network.
Ijms 21 09359 g001
Ijms 21 09359 g002 550
Figure 2. Functions of differentially expressed (DE) mRNAs were analyzed with DAVID and KOBAS bioinformatics resources, respectively. (A,B) Top 15 Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) results of upregulated mRNAs, respectively. (C,D) Top 15 GO and KEGG results of downregulated mRNAs, respectively. Color from blue to red indicates −log10 (p) from low to high. X-axes represent the number of genes involved in each pathway.
Figure 2. Functions of differentially expressed (DE) mRNAs were analyzed with DAVID and KOBAS bioinformatics resources, respectively. (A,B) Top 15 Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) results of upregulated mRNAs, respectively. (C,D) Top 15 GO and KEGG results of downregulated mRNAs, respectively. Color from blue to red indicates −log10 (p) from low to high. X-axes represent the number of genes involved in each pathway.
Ijms 21 09359 g002
Ijms 21 09359 g003 550
Figure 3. Metastasis-associated ceRNA network for colorectal cancer (CRC). (A) Green squares represent downregulated miRNAs, red ellipses represent upregulated mRNAs, and purple diamonds represent upregulated lncRNAs. (B) Red squares represent upregulated miRNAs, green ellipses represent downregulated mRNAs, and blue diamonds represent downregulated lncRNAs.
Figure 3. Metastasis-associated ceRNA network for colorectal cancer (CRC). (A) Green squares represent downregulated miRNAs, red ellipses represent upregulated mRNAs, and purple diamonds represent upregulated lncRNAs. (B) Red squares represent upregulated miRNAs, green ellipses represent downregulated mRNAs, and blue diamonds represent downregulated lncRNAs.
Ijms 21 09359 g003
Ijms 21 09359 g004 550
Figure 4. Prognostic evaluation of the 3-lncRNA (long noncoding RNA) signature. (A) CRC patients in the training set were divided into high- and low-risk groups using the risk score generated from the 3-lncRNA signature. The heatmap shows the expression profiles of the 3 lncRNAs in CRC patients. Color from green to red indicates an expression level from low to high, respectively. (B) The Kaplan–Meier curves of the 3-lncRNA signature for overall survival among CRC patients in the training set (p = 0.002). (C) Receiver operating characteristic (ROC) curve for predicting 5-year survival in CRC patients in the training set based on risk scores. (D) CRC patients in the testing set were divided into high- and low-risk groups using the risk score generated from the 3-lncRNA signature. (E) The Kaplan–Meier curves of the 3-lncRNA signature for overall survival among CRC patients in the testing set (p = 0.02). (F) ROC curve for predicting 5-year survival in CRC patients in the testing set based on risk scores.
Figure 4. Prognostic evaluation of the 3-lncRNA (long noncoding RNA) signature. (A) CRC patients in the training set were divided into high- and low-risk groups using the risk score generated from the 3-lncRNA signature. The heatmap shows the expression profiles of the 3 lncRNAs in CRC patients. Color from green to red indicates an expression level from low to high, respectively. (B) The Kaplan–Meier curves of the 3-lncRNA signature for overall survival among CRC patients in the training set (p = 0.002). (C) Receiver operating characteristic (ROC) curve for predicting 5-year survival in CRC patients in the training set based on risk scores. (D) CRC patients in the testing set were divided into high- and low-risk groups using the risk score generated from the 3-lncRNA signature. (E) The Kaplan–Meier curves of the 3-lncRNA signature for overall survival among CRC patients in the testing set (p = 0.02). (F) ROC curve for predicting 5-year survival in CRC patients in the testing set based on risk scores.
Ijms 21 09359 g004
Ijms 21 09359 g005 550
Figure 5. Clinical feature analysis of HOTAIR, LINC00261, and LINC00114. (A) Kaplan–Meier analysis of CRC overall survival according to the three lncRNAs. Percent survival is shown on the y-axis, and time is shown on the x-axis. Red lines represent patients with higher expression levels of lncRNAs, and blue lines represent patients with lower expression levels of lncRNAs (p < 0.05 log-rank test). (B) LncRNA expression analysis in CRC tissues and normal tissues. (C) LncRNA expression analysis in metastatic CRC tissues (M1) and nonmetastatic samples (M0). (D) LncRNA expression analysis in stages 1–4 of CRC. Significant differences were assessed by nonpaired t-tests. n.s., not significant; * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 5. Clinical feature analysis of HOTAIR, LINC00261, and LINC00114. (A) Kaplan–Meier analysis of CRC overall survival according to the three lncRNAs. Percent survival is shown on the y-axis, and time is shown on the x-axis. Red lines represent patients with higher expression levels of lncRNAs, and blue lines represent patients with lower expression levels of lncRNAs (p < 0.05 log-rank test). (B) LncRNA expression analysis in CRC tissues and normal tissues. (C) LncRNA expression analysis in metastatic CRC tissues (M1) and nonmetastatic samples (M0). (D) LncRNA expression analysis in stages 1–4 of CRC. Significant differences were assessed by nonpaired t-tests. n.s., not significant; * p < 0.05, ** p < 0.01, *** p < 0.001.
Ijms 21 09359 g005
Ijms 21 09359 g006 550
Figure 6. Enrichment analysis of HOTAIR, LINC00261, and LINC00114. (A,B) Top 15 GO and KEGG results for coexpressed genes of HOTAIR. (C,D) Top 15 GO and KEGG results of the coexpressed genes of LINC00261. Color from blue to red indicates −log10 (p) from low to high. X-axes represent the number of genes involved in each pathway. (E,F) The gene sets “KRAS.300_UP.V1_UP” (p < 0.0001) and “TGF_BETA_SIGNALING” (p = 0.002) were significantly enriched with low levels of LINC00114. (G) The gene sets “PTEN_DN.V1_UP” (p < 0.0001), which included genes upregulated upon knockdown of PTEN, were significantly enriched with low levels of LINC00114. (H) The gene set “GNF2_MSH6” (p = 0.018) was enriched with high levels of LINC00114.
Figure 6. Enrichment analysis of HOTAIR, LINC00261, and LINC00114. (A,B) Top 15 GO and KEGG results for coexpressed genes of HOTAIR. (C,D) Top 15 GO and KEGG results of the coexpressed genes of LINC00261. Color from blue to red indicates −log10 (p) from low to high. X-axes represent the number of genes involved in each pathway. (E,F) The gene sets “KRAS.300_UP.V1_UP” (p < 0.0001) and “TGF_BETA_SIGNALING” (p = 0.002) were significantly enriched with low levels of LINC00114. (G) The gene sets “PTEN_DN.V1_UP” (p < 0.0001), which included genes upregulated upon knockdown of PTEN, were significantly enriched with low levels of LINC00114. (H) The gene set “GNF2_MSH6” (p = 0.018) was enriched with high levels of LINC00114.
Ijms 21 09359 g006
Ijms 21 09359 g007 550
Figure 7. LINC00114 represses CRC cell migration. (A,B) LINC00114 expression levels were determined in negative control or siRNA-transfected HCT-116 or LOVO cells by qRT-PCR. Expression levels were normalized to GAPDH, and the results were expressed as the fold increase relative to negative control. (C,D) Transwell assays were used to determine the role of LINC00114 siRNA in HCT-116 or LOVO cell migration and invasion. All data are expressed as mean ± standard deviation (SD) of three independent experiments; * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 7. LINC00114 represses CRC cell migration. (A,B) LINC00114 expression levels were determined in negative control or siRNA-transfected HCT-116 or LOVO cells by qRT-PCR. Expression levels were normalized to GAPDH, and the results were expressed as the fold increase relative to negative control. (C,D) Transwell assays were used to determine the role of LINC00114 siRNA in HCT-116 or LOVO cell migration and invasion. All data are expressed as mean ± standard deviation (SD) of three independent experiments; * p < 0.05, ** p < 0.01, *** p < 0.001.
Ijms 21 09359 g007
Ijms 21 09359 g008 550
Figure 8. LINC00114 represses CRC cell proliferation. (A–D). Cell Counting Kit (CCK)-8 and colony formation assays were used to assess the role of LINC00114 siRNA in HCT-116 or LOVO cell proliferation. All data are expressed as mean ± SD of three independent experiments; ** p < 0.01, *** p < 0.001.
Figure 8. LINC00114 represses CRC cell proliferation. (A–D). Cell Counting Kit (CCK)-8 and colony formation assays were used to assess the role of LINC00114 siRNA in HCT-116 or LOVO cell proliferation. All data are expressed as mean ± SD of three independent experiments; ** p < 0.01, *** p < 0.001.
Ijms 21 09359 g008
Ijms 21 09359 g009 550
Figure 9. LINC00114 is the direct target of miR-135a. (A) Pearson’s correlation analysis showed that miR-135a was negatively associated with LINC00114. (B) The levels of LINC00114 were determined in negative control or miR-135a-mimic-transfected LOVO cells by qPCR. (C) Schematic illustration shows the miR-135a binding site and the reverse complementary sequence mutations in LINC00114. Luciferase assay was used to identify the miR-135a binding site. All data are expressed as mean ± SD of three independent experiments; * p < 0.05, *** p < 0.001.
Figure 9. LINC00114 is the direct target of miR-135a. (A) Pearson’s correlation analysis showed that miR-135a was negatively associated with LINC00114. (B) The levels of LINC00114 were determined in negative control or miR-135a-mimic-transfected LOVO cells by qPCR. (C) Schematic illustration shows the miR-135a binding site and the reverse complementary sequence mutations in LINC00114. Luciferase assay was used to identify the miR-135a binding site. All data are expressed as mean ± SD of three independent experiments; * p < 0.05, *** p < 0.001.
Ijms 21 09359 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, S.; Cao, Q.; An, G.; Yan, B.; Lei, L. Identification of the 3-lncRNA Signature as a Prognostic Biomarker for Colorectal Cancer. Int. J. Mol. Sci. 2020, 21, 9359. https://doi.org/10.3390/ijms21249359

AMA Style

Liu S, Cao Q, An G, Yan B, Lei L. Identification of the 3-lncRNA Signature as a Prognostic Biomarker for Colorectal Cancer. International Journal of Molecular Sciences. 2020; 21(24):9359. https://doi.org/10.3390/ijms21249359

Chicago/Turabian Style

Liu, Shuzhen, Qing Cao, Guoyan An, Bianbian Yan, and Lei Lei. 2020. "Identification of the 3-lncRNA Signature as a Prognostic Biomarker for Colorectal Cancer" International Journal of Molecular Sciences 21, no. 24: 9359. https://doi.org/10.3390/ijms21249359

APA Style

Liu, S., Cao, Q., An, G., Yan, B., & Lei, L. (2020). Identification of the 3-lncRNA Signature as a Prognostic Biomarker for Colorectal Cancer. International Journal of Molecular Sciences, 21(24), 9359. https://doi.org/10.3390/ijms21249359

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop