Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 06 February 2020
Sec. Molecular and Cellular Oncology

A 6 lncRNA-Based Risk Score System for Predicting the Recurrence of Colon Adenocarcinoma Patients

\nHaojie Yang&#x;Haojie Yang1Hong-Cheng Lin&#x;Hong-Cheng Lin2Hua Liu&#x;Hua Liu1Dan GanDan Gan1Wei JinWei Jin1Can CuiCan Cui1Yixin YanYixin Yan3Yiming QianYiming Qian3Changpeng Han
Changpeng Han1*Zhenyi Wang
Zhenyi Wang1*
  • 1Department of Coloproctology, Yueyang Hospital of Integrated Traditional Chinese and Western Medicine, Shanghai University of Traditional Chinese Medicine, Shanghai, China
  • 2Department of Coloproctology, The Sixth Affiliated Hospital of Sun Yat-sen University (Gastrointestinal & Anal Hospital of Sun Yat-sen University), Guangzhou, China
  • 3Department of Emergency Medicine, Yueyang Hospital of Integrated Traditional Chinese and Western Medicine, Shanghai University of Traditional Chinese Medicine, Shanghai, China

Colon adenocarcinoma (COAD) is a common type of colon cancer, and post-operative recurrence and metastasis may occur in COAD patients. This study is designed to build a risk score system for COAD patients. The Cancer Genome Atlas (TCGA) dataset of COAD (the training set) was downloaded, and GSE17538 and GSE39582 (the validation sets) from Gene Expression Omnibus database were obtained. The differentially expressed RNAs (DERs) were analyzed by limma package. Using survival package, the independent prognosis-associated long non-coding RNAs (lncRNAs) were selected for constructing risk score system. After the independent clinical prognostic factors were screened out using survival package, a nomogram survival model was constructed using rms package. Furthermore, competitive endogenous RNA (ceRNA) regulatory network and enrichment analyses separately were performed using Cytoscape software and DAVID tool. Totally 404 DERs between recurrence and non-recurrence groups were identified. Based on the six independent prognosis-associated lncRNAs (including H19, KCNJ2-AS1, LINC00899, LINC01503, PRKAG2-AS1, and SRRM2-AS1), the risk score system was constructed. After the independent clinical prognostic factors (Pathologic M, pathologic T, and RS model status) were identified, the nomogram survival model was built. In the ceRNA regulatory network, there were three lncRNAs, four miRNAs, and 77 mRNAs. Additionally, PPAR signaling pathway and hedgehog signaling pathway were enriched for the mRNAs in the ceRNA regulatory network. The risk score system and the nomogram survival model might be used for predicting COAD recurrence. Besides, PPAR signaling pathway and hedgehog signaling pathway might affect the recurrence of COAD patients.

Introduction

As a common malignancy occurring in the colon, colon cancer is divided into adenocarcinoma, undifferentiated carcinoma, and mucinous adenocarcinoma (1). Colon cancer is induced mainly by a high-fat low-fiber diet, and its symptoms include hematochezia, purulent stools, the change of bowel habits, abdominal mass, ileus, bellyache, and anemia (2). Colon cancer ranks third among the gastrointestinal cancers, and its morbidity and mortality rates in men are higher than those in women (3). The 5-years survival rate of colon cancer after radical operation is about 60–70%, and the main reasons leading to the failure of surgical treatment are post-operative recurrence and metastasis (4, 5). Therefore, the discovery of valuable molecular markers is important for the early prediction and treatment of colon cancer recurrence.

Long non-coding RNA (lncRNA) exerts a key regulatory effect in the development and progression of some diseases (6). In recent years, a variety of lncRNAs have been reported to be dysregulated in colon cancer. For example, the lncRNA small ubiquitin-like modifier 1 pseudogene 3 (SUMO1P3) is overexpressed in colon cancer, and is positively related to angiogenesis, metastases, advanced histological stages, and unfavorable outcome of colon cancer patients (7). Through activating the Wnt/β-catenin signaling, the lncRNA breast cancer anti-estrogen resistance 4 (BCAR4) accelerates the progression of colon cancer by promoting cell proliferation and suppressing cell apoptosis (8). The lncRNA cytoskeleton regulator RNA (CYTOR) facilitates epithelial-to-mesenchymal transition (EMT) and metastasis in colon cancer via interacting with β-catenin, which predicts poor prognosis and may be promising therapeutic target of colon cancer (9, 10). These results suggested that lncRNAs might be used as biomarkers for predicting the prognosis of colon cancer.

Along with the development of bioinformatics, several researchers have developed lncRNA-related signatures for predicting prognosis of colon cancer. Xue et al. proposed a two-lncRNA expression signature to predict survival of patients with colon adenocarcinoma (11). Xing et al. developed a 14-lncRNA prognostic signature for patients with colon adenocarcinoma (12). Lv et al. identified a five-lncRNA prognostic signature to predict the survival of colon cancer patients (13). These studies proved the strong power of prognostic prediction value of the lncRNA signatures for colon cancer patients. However, seldom studies were conducted to develop lncRNA signatures for predicting recurrence of colon cancer.

As a competitive endogenous RNA (ceRNA), lncRNA can mediate gene expression through acting as microRNA (miRNA) sponge (14, 15). The up-regulation of the lncRNA metastasis associated lung adenocarcinoma transcript 1 (MALAT1) mediates high mobility group box 1 (HMGB1) expression in colon cancer via competing with miR-129-5p, and the MALAT1/miR-129-5p/HMGB1 axis serves as a key prognostic marker in the development of the tumor (16). In this study, the lncRNAs, miRNAs, and mRNAs with differential expression between recurrence and non-recurrence colon adenocarcinoma (COAD) samples were analyzed. The recurrence prognosis-associated lncRNAs were screened, and then the independent prognosis-associated lncRNAs were further selected for constructing risk score system. Moreover, nomogram survival model construction, ceRNA regulatory network construction and enrichment analysis were conducted. Our findings might be conducive to predicting the recurrence of COAD patients.

Results

Differential Expression Analysis

A total of 13,834 mRNAs, 827 lncRNAs, and 1,037 miRNAs were annotated from the TCGA transcriptomic RNA and miRNA datasets. The 310 COAD samples in the TCGA dataset were classified into recurrence (66 samples) and non-recurrence (244 samples) groups. Under the defined thresholds, 404 DERs (including 357 DE-mRNAs (122 down-regulated and 235 up-regulated), 26 DE-lncRNAs (eight down-regulated and 18 up-regulated), and 21 DE-miRNAs (eight down-regulated and 13 up-regulated) between recurrence and non-recurrence groups were screened out (Figure 1A). Based on the expression of the DERs, the clustering heatmap is drew and presented in Figure 1B.

FIGURE 1
www.frontiersin.org

Figure 1. The screening results of differentially expressed RNAs (DERs). (A) The volcano plot (left; the horizontal dashed line represents false discovery rate (FDR) < 0.05, and the vertical dashed lines represent |log2 fold change (FC)| > 0.263; the blue dots indicate DERs; FC: fold change) and the histogram showing the proportional distribution of different kinds of DERs (right; pink and blue separately represent up-regulation and down-regulation; lncRNA: long non-coding RNA, miRNA: microRNA); (B) The hierarchical clustering heatmap (in sample strip, blue and pink separately represent recurrence samples and non-recurrence samples).

Construction of Risk Score System

Based on the univariable Cox regression analysis, a total of 21 DE-lncRNAs were found to be significantly related to recurrence prognosis (Table 1). From the 21 recurrence prognosis-associated lncRNAs, six independent prognosis-associated lncRNAs (including H19 imprinted maternally expressed transcript, H19; KCNJ2 antisense RNA 1, KCNJ2-AS1; long intergenic non-protein coding RNA 899, LINC00899; long intergenic non-protein coding RNA 1503, LINC01503; PRKAG2 antisense RNA 1, PRKAG2-AS1; and SRRM2 antisense RNA 1, SRRM2-AS1) were further identified (Table 2).

TABLE 1
www.frontiersin.org

Table 1. Univariable Cox regression analysis identified 21 DE-lncRNAs related to recurrence prognosis.

TABLE 2
www.frontiersin.org

Table 2. The information of the six independent prognosis-associated long non-coding RNAs (lncRNAs).

Subsequently, the risk score system based on the independent prognosis-associated lncRNAs was constructed, and the relevant formula was as follow:

RS=(0.1647)×ExpH19 + (-1.1438)×ExpKCNJ2-AS1     + (1.4922)×ExpLINC00899 + (0.8525)×ExpLINC01503     + (-0.8100)×ExpPRKAG2-AS1 + (3.8523)×ExpSRRM2-AS1

The median of the RSs of samples were calculated, and the samples in the TCGA dataset and the validation sets separately were divided into high and low risk group. Then, the correlation between the grouping and the actual recurrence prognosis information was evaluated using KM curves. The results showed that the grouping based on the risk score system had significant correlations with the actual recurrence prognosis in the TCGA dataset and the validation sets (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. The Kaplan-Meier (KM) curves and receiver operating characteristic (ROC) curves showing the correlations between the grouping based on the risk score system and the actual recurrence prognosis. (A) The KM curves for The Cancer Genome Atlas (TCGA) dataset; (B) The KM curves for the validation set GSE17538; (C) The KM curves for the validation set GSE39582; (D) The ROC curves for the TCGA dataset and the validation sets; In KM curves, blue and red curves separately represent low risk group and high risk group. In ROC curves, red, black, and blue curves represent the TCGA dataset, the validation set GSE17538, and the validation set GSE39582, respectively. HR, hazard ratio; AUC, area under the receiver operating characteristic curve.

Establishment and Validation of Nomogram Survival Model

By using univariable and multivariable Cox regression analyses, pathologic M, pathologic T, and RS model status were selected as the independent clinical prognostic factors in the COAD samples (Table 3). The COAD samples in lower pathologic M and pathologic T stages had better recurrence prognosis, which was consistent with clinical facts (Figure 3).

TABLE 3
www.frontiersin.org

Table 3. The selection of independent clinical prognostic factors.

FIGURE 3
www.frontiersin.org

Figure 3. The Kaplan-Meier (KM) curves showed the associations between pathologic M/pathologic T and recurrence prognosis. (A) The KM curves for pathologic M (blue and red curves separately represent the samples in pathologic M0 stage and pathologic M1 stage); (B) The KM curves for pathologic T (black, red, blue, and purples curves represent the samples in pathologic T1 stage, pathologic T2 stage, pathologic T3 stage, and pathologic T4 stage, respectively). HR, hazard ratio.

Moreover, a nomogram survival model was built based on the three independent prognostic factors (Figure 4A). The 3-years survival probability and 5-years survival probability of patients could be easily calculated based on their pathologic M, pathologic T and RS. Additionally, the nomogram-predicted 3-years survival probability/5-years survival probability was further compared with the actual 3-years survival probability/5-years survival probability recorded in TCGA and results suggested that there is high agreement between nomogram-predicted probability of recurrence free survival and the actual recurrence free survival (Figure 4B), indicating the good performance of nomogram survival model.

FIGURE 4
www.frontiersin.org

Figure 4. Nomogram survival model. (A) The nomogram survival model. (B) The graph showing the consistency of the nomogram-predicted survival probability and the actual survival probability (the horizontal and vertical axes represent the predicted survival probability and the actual survival probability, respectively; black and red segments separately represent the 5-years survival rate and the 3-years survival rate in the group with the highest consistency).

CeRNA Regulatory Network Construction and Enrichment Analysis

The lncRNA-miRNA regulatory network was constructed, which had seven nodes (three lncRNAs and four miRNAs) (Figure 5). Meanwhile, the miRNA-mRNA regulatory network was built, involving 81 nodes (four miRNAs and 77 mRNAs) (Figure 6). Based on the lncRNA-miRNA-mRNA relationships, the ceRNA regulatory network (84 nodes, including three lncRNAs (two up-regulated and one down-regulated), four miRNAs (two up-regulated and two down-regulated), and 77 mRNAs (21 up-regulated and 56 down-regulated) was constructed (Figure 7).

FIGURE 5
www.frontiersin.org

Figure 5. The long non-coding RNA (lncRNA)-microRNA (miRNA) regulatory network. Squares and triangles represent lncRNAs and miRNAs, respectively. Red and green separately represent up-regulation and down-regulation. FC, fold change.

FIGURE 6
www.frontiersin.org

Figure 6. The microRNA (miRNA)-mRNA regulatory network. Circles and triangles represent mRNAs and miRNAs, respectively. Red and green separately represent up-regulation and down-regulation. FC, fold change.

FIGURE 7
www.frontiersin.org

Figure 7. The competing endogenous RNA (ceRNA) regulatory network. Squares, circles, and triangles represent long non-coding RNAs (lncRNAs), mRNAs, and microRNAs (miRNAs), respectively. Red and green separately represent up-regulation and down-regulation. Red and black edges represent lncRNA-miRNA and miRNA-mRNA pairs, respectively. FC, fold change.

Based on DAVID online tool, 13 GO functional terms (such as regulation of neuron differentiation, p-value = 1.68E-03; regulation of neurogenesis, p-value = 3.74E-03; and inorganic anion transport, p-value = 5.45E-03) and five KEGG pathways (such as PPAR signaling pathway, p-value = 3.82E-03; neuroactive ligand-receptor interaction, p-value = 1.06E-02; intestinal immune network for IgA production, p-value = 2.00E-02; and hedgehog signaling pathway, p-value = 2.25E-02) were enriched for the mRNAs in the ceRNA regulatory network (Table 4).

TABLE 4
www.frontiersin.org

Table 4. The enrichment results for the mRNAs in the competitive endogenous RNA (ceRNA) regulatory network.

Discussion

In the present study, a total of 404 DERs (including 357 DE-mRNAs, 26 DE-lncRNAs, and 21 DE-miRNAs) between recurrence and non-recurrence groups were identified. After the 21 recurrence prognosis-associated lncRNAs were screened out, six independent prognosis-associated lncRNAs (including H19, KCNJ2-AS1, LINC00899, LINC01503, PRKAG2-AS1, and SRRM2-AS1) were further selected. Based on the independent prognosis-associated lncRNAs, the RS system was constructed. Pathologic M, pathologic T, and RS model status were identified as the independent clinical prognostic factors, and then the nomogram survival model was built. This nomogram survival model might be used to predict 3/5 years recurrence free survival in future clinical practice, by combining pathologic M, pathologic T and the RS calculated from the expression level of the 6 lncRNAs detected by qRT-PCR from surgical specimens.

High mobility group AT-hook 1 (HMGA1) is inhibited by H19 short hairpin RNA (shRNA) and is promoted by miR-138 inhibitor, and the H19-miR138-HMGA1 pathway plays important roles in mediating the invasion and migration of colon cancer (17). H19 expression is significantly up-regulated in immunodeficient mice induced by colon cancer cells, and H19 may be taken as a novel therapeutic target in colon cancer (18). H19 can suppress vitamin D receptor (VDR) expression via miR-675-5p, and increased H19 leads to the resistance to 1,25(OH)2D3 treatment in the advanced colon cancer (19). LINC00899 is elevated in the serum and bone marrow of acute myeloid leukemia (AML) patients, therefore, serum LINC00899 may be a promising marker for the early diagnosis and prognosis of AML (20). These suggested that H19 and LINC00899 might be involved in the recurrence prognosis of COAD. Thirteen lncRNAs (including KCNJ2-AS1) are selected as prognostic biomarkers, based on which a prognostic signature is constructed for predicting the disease free survival in gastric cancer patients (21). Potassium voltage-gated channel subfamily J member 2 (KCNJ2)/Kir2.1 belongs to the inwardly rectifying potassium channel family, which regulates drug resistance and cell growth through mediating mitochondrial 37S ribosomal protein MRP1 (MRP1)/ATP binding cassette subfamily C member 1 (ABCC1) expression in small-cell lung cancer (22). Up-regulated LINC01503 contributes to cell proliferation, growth, invasion, and migration in esophageal squamous cell carcinoma (ESCC), which may be used as a marker of aggressive ESCC (23). Overexpressed LINC01503 regulates forkhead box K1 (FOXK1) expression by competing with miR-4492, which promotes cell proliferation and invasion in colorectal cancer (24). Therefore, KCNJ2-AS1 and LINC01503 might also be correlated with the recurrence prognosis of COAD.

The genetic variation in a candidate pathway contributes to the risk of both colon and rectal cancers, and protein kinase AMP-activated non-catalytic subunit gamma 2 (PRKAG2), mechanistic target of rapamycin kinase (FRAP1), TSC complex subunit 2 (TSC2), and protein kinase AMP-activated catalytic subunit alpha 1 (PRKAA1) genes involved in the pathway are significantly related to colon cancer (25, 26). Through influencing the alternative splicing of downstream genes, the S346F mutation in serine/arginine repetitive matrix 2 (SRRM2) promotes the susceptibility of papillary thyroid carcinoma (27). Erb-b2 receptor tyrosine kinase 2 (ERBB2) overexpression is correlated with increased metastasis of human cancers, and the depletion of SRRM2, splicing factor, arginine/serine-rich 1 (SFRS1), SFRS9, and SFRS10 proteins reduced the migration rate of ovarian cancer cells overexpressing ERBB2 (28). These indicated that PRKAG2-AS1 and SRRM2-AS1 might play roles in the recurrence prognosis of COAD.

Through interacting with T cell transcription factor-4 (Tcf-4) and beta-catenin, peroxisome proliferator activated receptor gamma (PPAR gamma) may determine colon cell fate and serve as a target of the Wnt pathway in colon cancer cells (29). Hedgehog signaling pathway functions in gastrointestinal development and affects the formation of multiple tumors, which predicts unfavorable outcomes in colon cancer patients (30). The mRNAs in the ceRNA regulatory network were involved in PPAR signaling pathway and hedgehog signaling pathway, suggesting that PPAR signaling pathway and hedgehog signaling pathway might be related to the recurrence prognosis of COAD.

However, there are some limitations in this study. First, the differential expression of the 6 lncRNAs was identified from RNA sequence data. Experiment validation of these lncRNAs in colon cancer patients should be conducted in further research. Second, the ceRNAs related with these lncRNAs should also be validated in in vitro and in vivo studies.

In conclusion, 357 DE-mRNAs, 26 DE-lncRNAs, and 21 DE-miRNAs were identified between recurrence and non-recurrence groups. Besides, the risk score system (involving H19, KCNJ2-AS1, LINC00899, LINC01503, PRKAG2-AS1, and SRRM2-AS1) and the nomogram survival model might be useful for the prediction of COAD recurrence. Moreover, PPAR signaling pathway and hedgehog signaling pathway might be correlated with the recurrence prognosis of COAD patients.

Materials and Methods

Data Source

From The Cancer Genome Atlas (TCGA) database (https://cancergenome.nih.gov/), the transcriptomic RNA (including 512 samples) and miRNA (including 461 samples) expression data of COAD (downloaded in May 15, 2019; platform: Illumina HiSeq 2000 RNA Sequencing) was extracted. The two groups of samples were paired according to the sample numbers, and then were corresponded to the downloaded clinical information. Finally, a total of 310 paired COAD samples with recurrence prognosis information were screened and used as the training set.

Taking “colon adenocarcinoma” and “Homo sapiens” as searching words, the eligible datasets were selected from Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/) according to the following criteria: (1) the samples were solid tissue samples from COAD patients; (2) the total sample size was no <200, and the number of COAD samples was not under 150; (3) the COAD samples had recurrence prognosis information. After searching, two datasets (GSE17538, 244 samples (including 145 COAD samples with recurrence prognosis information), platform: GPL570 Affymetrix Human Genome U133 Plus 2.0 Array, the validation set 1; and GSE39582, 585 samples (including 519 COAD samples with recurrence prognosis information), platform: GPL570 Affymetrix Human Genome U133 Plus 2.0 Array, the validation set 2) that met the requirements were downloaded.

Differential Expression Analysis

Using HUGO Gene Nomenclature Committee (HGNC) database (http://www.genenames.org/) (31), the lncRNAs and mRNAs in the datasets were recognized. Based on the source information, the samples in the TCGA dataset were classified into recurrence and non-recurrence groups. The R package limma (version 3.34.7, https://bioconductor.org/packages/release/bioc/html/limma.html) (32) was applied for exploring the differentially expressed RNAs [DERs, including differentially expressed lncRNAs (DE-lncRNAs), differentially expressed mRNAs (DE-mRNAs), and differentially expressed miRNAs (DE-miRNAs)] between recurrence and non-recurrence groups. The false discovery rate (FDR) < 0.05 and |log2 fold change (FC)| > 0.263 were taken as the screening criteria. Using the R package pheatmap (version 1.0.8, https://cran.r-project.org/web/packages/pheatmap/index.html) (33), the expression values of the DERs in TCGA dataset were conducted with bidirectional hierarchical clustering.

Construction of Risk Score System

Based on the recurrence information of the 310 COAD samples in the TCGA dataset, the DE-lncRNAs having significant associations with recurrence prognosis were selected combined with the univariable Cox regression analysis in the survival package (version 2.41-1, http://bioconductor.org/packages/survivalr/) in R (34). The threshold was set as log-rank p-value < 0.05.

From the recurrence prognosis-associated lncRNAs, the independent prognosis-associated lncRNAs were further selected using the multivariable Cox regression analysis in the survival package in R (34). Afterwards, the risk score system was built based on the expression levels and independent prognostic coefficients of the independent prognosis-associated lncRNAs. The risk scores (RSs) of the COAD samples were calculated using the formula below:

Risk score (RS) = βlncRNA×ExplncRNA

β lncRNA indicates the independent prognostic coefficient of independent prognosis-associated lncRNA, and ExplncRNA stands for the expression level of independent prognosis-associated lncRNA.

The median of the RSs of the COAD samples in the TCGA dataset were calculated and used to divide the samples into high risk group and low risk group. Using the Kaplan-Meier (KM) curve in the R package survival (34), the association between the grouping and the actual recurrence prognosis information was analyzed. Similarly, the prediction efficiencies of the risk score system in the validation sets were evaluated.

Construction of Nomogram Survival Model

Using the univariable and multivariable Cox regression analyses in the R package survival (34), the independent clinical prognostic factors in the COAD samples in the TCGA dataset were identified. The log-rank p-value < 0.05 was selected as the significant threshold. Using the R package rms (version 5.1-2, https://cran.r-project.org/web/packages/rms/index.html) (35), nomogram survival model was constructed based on the independent clinical prognostic factors and the predicted risk information.

CeRNA Regulatory Network Analysis and Enrichment Analysis

The DIANA-LncBase v2 database (http://carolina.imis.athena-innovation.gr/diana_tools/web/index.php?r=lncbasev2%2Findex-experimental) (36) was used to search the relationships between the DE-miRNAs and the independent prognosis-associated lncRNAs. The pairs involving miRNAs and lncRNAs with opposite expression directions were selected, based on which the lncRNA-miRNA regulatory network was visualized using Cytoscape software [version 3.6.1, https://cytoscape.org/ (37)].

The starBase database (version 2.0, http://starbase.sysu.edu.cn/, including the information of RNA22, PITA, targetScan, picTar, and miRanda databases) (38) was applied for searching the targets of the miRNAs implicated in the lncRNA-miRNA regulatory network. The regulatory relationships included by at least two of RNA22, PITA, targetScan, picTar, and miRanda databases were taken as the miRNA-target pairs. Subsequently, the pairs involving negatively correlated miRNAs and DE-mRNAs were further screened. Moreover, Cytoscape software (37) was utilized to visualize the miRNA-mRNA regulatory network.

The lncRNA-miRNA-mRNA regulatory relationships were obtained through integrating the lncRNA-miRNA and miRNA-mRNA pairs. Based on the lncRNA-miRNA-mRNA regulatory relationships, the ceRNA regulatory network was built by Cytoscape software (37). Using DAVID online tool (version 6.8, https://david.ncifcrf.gov/) (39), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) functional enrichment analyses for the mRNAs in the ceRNA regulatory network were carried out. The p-value < 0.05 was selected as the threshold of enrichment significance.

Data Availability Statement

The datasets generated for this study can be found in The Cancer Genome Atlas (https://portal.gdc.cancer.gov/); the NCBI Gene Expression Omnibus (GSE17538, GSE39582).

Author Contributions

CH and ZW: funding acquisition. HY, CH, and ZW: investigation. HY, H-CL, HL, DG, YQ, CH, and ZW: methodology. H-CL, CH, and ZW: project administration. HY, DG, WJ, CC, and YY: software. HL: supervision. H-CL: validation. HY, HL, CH, and ZW: writing–original draft.

Funding

This study was supported by National Natural Science Foundation of China Grant 81874468, Peak Research Team Project in Shanghai University of Traditional Chinese Medicine, Three Years Plan of Action to further accelerate the Development of Traditional Chinese Medicine in Shanghai under Grant ZY (2018-2020)-CCCX-2004-09.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Huaping Liu for critical reading of this manuscript and providing helpful suggestions.

References

1. Obaro AE, Burling DN, Plumb AA. Colon cancer screening with CT colonography: logistics, cost effectiveness, efficiency and progress. Br J Radiol. (2018) 91:20180307. doi: 10.1259/bjr.20180307

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Sunaabbcda V. Quality of life and barriers to symptom management in colon cancer. Eur J Oncol Nurs. (2012) 16:276–80. doi: 10.1016/j.ejon.2011.06.011

CrossRef Full Text | Google Scholar

3. Mcguire S. World cancer report 2014. Geneva, Switzerland: World Health Organization, International Agency for Research on Cancer, WHO Press, 2015. Adv Nutr. (2016) 7:418. doi: 10.3945/an.116.012211

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Mejri N, Dridi M, El BH, Labidi S, Daoud N, Boussen H. Tumor location impact in stage II and III colon cancer: epidemiological and outcome evaluation. J Gastrointest Oncol. (2018) 9:263–8. doi: 10.21037/jgo.2017.12.02

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Mody K, Bekaii-Saab T. Clinical trials and progress in metastatic colon cancer. Surg Oncol Clin N Am. (2018) 27:349–65. doi: 10.1016/j.soc.2017.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Bhan A, Soleimani M, Mandal SS. Long noncoding RNA and cancer: a new paradigm. Cancer Res. (2017) 77:3965. doi: 10.1158/0008-5472.CAN-16-2634

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Zhang LM, Wang P, Liu XM, Zhang YJ. LncRNA SUMO1P3 drives colon cancer growth, metastasis and angiogenesis. Am J Transl Res. (2017) 9:5461–72.

PubMed Abstract | Google Scholar

8. Ouyang S, Zheng X, Zhou X, Chen Z, Yang X, Xie M. LncRNA BCAR4 promotes colon cancer progression via activating Wnt/β-catenin signaling. Oncotarget. (2017) 8:92815–26. doi: 10.18632/oncotarget.21590

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Liang J, Wei X, Liu Z, Cao D, Tang Y, Zou Z, et al. Long noncoding RNA CYTOR in cancer: a TCGA data review. Clin Chim Acta. (2018) 483:227. doi: 10.1016/j.cca.2018.05.010

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Yue B, Liu C, Sun H, Liu M, Song C, Cui R, et al. A positive feed-forward loop between LncRNA-CYTOR and Wnt/β-catenin signaling promotes metastasis of colon cancer. Mol Ther. (2018) 26:1287–98. doi: 10.1016/j.ymthe.2018.02.024

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Xue W, Li J, Wang F, Han P, Liu Y, Cui B. A long non-coding RNA expression signature to predict survival of patients with colon adenocarcinoma. Oncotarget. (2017) 8:101298–308. doi: 10.18632/oncotarget.21064

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Xing Y, Zhao Z, Zhu Y, Zhao L, Zhu A, Piao D. Comprehensive analysis of differential expression profiles of mRNAs and lncRNAs and identification of a 14-lncRNA prognostic signature for patients with colon adenocarcinoma. Oncol Rep. (2018) 39:2365–75. doi: 10.3892/or.2018.6324

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Lv J, Guo Y, Yan L, Lu Y, Liu D, Niu J. Development and validation of a five-lncRNA signature with prognostic value in colon cancer. J Cell Biochem. (2019). doi: 10.1002/jcb.29518. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

14. An Y, Furber KL, Ji S. Pseudogenes regulate parental gene expression via ceRNA network. J Cell Mol Med. (2017) 21:185–92. doi: 10.1111/jcmm.12952

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Chiu HS, Martínez MR, Bansal M, Subramanian A, Golub TR, Yang X, et al. High-throughput validation of ceRNA regulatory networks. BMC Genomics. (2017) 18:418. doi: 10.1186/s12864-017-3790-7

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Jie Y, Zhao H. LncRNA MALAT1 induces colon cancer development by regulating miR-129-5p/HMGB1 axis. J Cell Physiol. (2017) 233:6750–7. doi: 10.1002/jcp.26383

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Yang Q, Wang X, Tang C, Chen X, He J. H19 promotes the migration and invasion of colon cancer by sponging miR-138 to upregulate the expression of HMGA1. Int J Oncol. (2017) 50:1801. doi: 10.3892/ijo.2017.3941

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Hu Q, Wang YB, Zeng P, Yan GQ, Xin L, Hu XY. Expression of long non-coding RNA (lncRNA) H19 in immunodeficient mice induced with human colon cancer cells. Eur Rev Med Pharmacol Sci. (2016) 20:4880–4.

PubMed Abstract | Google Scholar

19. Chen S, Bu D, Ma Y, Zhu J, Chen G, Sun L, et al. H19 overexpression induces resistance to 1,25(OH)2D3 by targeting VDR through miR-675-5p in colon cancer cells. Neoplasia. (2017) 19:226–36. doi: 10.1016/j.neo.2016.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Wang Y, Li Y, Song HQ, Sun GW. Long non-coding RNA LINC00899 as a novel serum biomarker for diagnosis and prognosis prediction of acute myeloid leukemia. Eur Rev Med Pharmacol Sci. (2018) 22:7364–70. doi: 10.26355/eurrev_201811_16274

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Cheng C, Wang Q, Zhu M, Liu K, Zhang Z. Integrated analysis reveals potential long non-coding RNA biomarkers and their potential biological functions for disease free survival in gastric cancer patients. Cancer Cell Int. (2019) 19:123. doi: 10.1186/s12935-019-0846-6

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Liu H, Huang J, Peng J, Wu X, Zhang Y, Zhu W, et al. Upregulation of the inwardly rectifying potassium channel Kir2.1 (KCNJ2) modulates multidrug resistance of small-cell lung cancer under the regulation of miR-7 and the Ras/MAPK pathway. Mol Cancer. (2015) 14:59. doi: 10.1186/s12943-015-0298-0

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Xie JJ, Jiang YY, Jiang Y, Li CQ, Lim MC, An O, et al. Super-enhancer-driven long non-coding RNA LINC01503, regulated by TP63, is over-expressed and oncogenic in squamous cell carcinoma. Gastroenterology. (2018) 154:2137–51 e31. doi: 10.1053/j.gastro.2018.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Lu SR, Li Q, Lu JL, Liu C, Xu X, Li JZ. Long non-coding RNA LINC01503 promotes colorectal cancer cell proliferation and invasion by regulating miR-4492/FOXK1 signaling. Exp Ther Med. (2018) 16:4879–85. doi: 10.3892/etm.2018.6775

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Slattery ML, Herrick JS, Lundgreen A, Fitzpatrick FA, Curtin K, Wolff RK. Genetic variation in a metabolic signaling pathway and colon and rectal cancer risk: mTOR, PTEN, STK11, RPKAA1, PRKAG2, TSC1, TSC2, PI3K and Akt1. Carcinogenesis. (2010) 31:1604–11. doi: 10.1093/carcin/bgq142

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Slattery ML, Abbie L, Herrick JS, Caan BJ, Potter JD, Wolff RK. Diet and colorectal cancer: analysis of a candidate pathway using SNPS, haplotypes, and multi-gene assessment. Nutr Cancer. (2011) 63:1226–34. doi: 10.1080/01635581.2011.607545

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Tomsic J, He H, Akagi K, Liyanarachchi S, Pan Q, Bertani B, et al. A germline mutation in SRRM2, a splicing factor gene, is implicated in papillary thyroid carcinoma predisposition. Sci Rep. (2015) 5:10566. doi: 10.1038/srep10566

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Mukherji M, Brill LM, Ficarro SB, Hampton GM, Schultz PG. A phosphoproteomic analysis of the ErbB2 receptor tyrosine kinase signaling pathways. Biochemistry. (2006) 45:15529–40. doi: 10.1021/bi060971c

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Jansson EA, Alexandra A, Gediminas G, Chun IK, Denise K, Velmurugesan A, et al. The Wnt/beta-catenin signaling pathway targets PPARgamma activity in colon cancer cells. Proc Natl Acad Sci USA. (2005) 102:1460. doi: 10.1073/pnas.0405928102

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Xu M, Li X, Liu T, Leng A, Zhang G. Prognostic value of hedgehog signaling pathway in patients with colon cancer. Med Oncol. (2012) 29:1010–6. doi: 10.1007/s12032-011-9899-7

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Povey S, Lovering R, Bruford E, Wright M, Lush M, Wain H, et al. The HUGO Gene Nomenclature Committee (HGNC). Hum Genet. (2001) 109:678. doi: 10.1007/s00439-001-0615-0

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Smyth GK. limma: Linear models for microarray data. In: Gentleman R, Carey VJ, Huber W, Irizarry RA, Dudoit S, editors. Bioinformatics & Computational Biology Solutions Using R & Bioconductor. New York, NY: Springer (2011). p. 397–420. doi: 10.1007/0-387-29362-0_23

CrossRef Full Text | Google Scholar

33. Wang L, Cao C, Ma Q, Zeng Q, Wang H, Cheng Z, et al. RNA-seq analyses of multiple meristems of soybean: novel and alternative transcripts, evolutionary and functional implications. BMC Plant Biol. (2014) 14:169. doi: 10.1186/1471-2229-14-169

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Diao R, Mu X, Wang T, Li S. Risk score based on ten lncRNA-mRNA expression predicts the survival of stage II–III colorectal carcinoma. PLoS ONE. (2017) 12:e0182908. doi: 10.1371/journal.pone.0182908

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Eng KH, Emily S, Kayla M. On representing the prognostic value of continuous gene expression biomarkers with the restricted mean survival curve. Oncotarget. (2015) 6:36308–18. doi: 10.18632/oncotarget.6121

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Paraskevopoulou MD, Vlachos IS, Karagkouni D, Georgakilas G, Kanellos I, Vergoulis T, et al. DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts. Nucleic Acids Res. (2016) 44:D231–8. doi: 10.1093/nar/gkv1270

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Kohl M, Wiese S, Warscheid B. Cytoscape: software for visualization and analysis of biological networks. Methods Mol Biol. (2011) 696:291–303. doi: 10.1007/978-1-60761-987-1_18

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Jun-Hao L, Shun L, Hui Z, Liang-Hu Q, Jian-Hua Y. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. (2014) 42:D92. doi: 10.1093/nar/gkt1248

CrossRef Full Text | Google Scholar

39. Wei HD, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. (2009) 37:1–13. doi: 10.1093/nar/gkn923

CrossRef Full Text | Google Scholar

Keywords: colon adenocarcinoma, differential expression analysis, risk score system, nomogram survival model, competitive endogenous RNA

Citation: Yang H, Lin H-C, Liu H, Gan D, Jin W, Cui C, Yan Y, Qian Y, Han C and Wang Z (2020) A 6 lncRNA-Based Risk Score System for Predicting the Recurrence of Colon Adenocarcinoma Patients. Front. Oncol. 10:81. doi: 10.3389/fonc.2020.00081

Received: 10 October 2019; Accepted: 16 January 2020;
Published: 06 February 2020.

Edited by:

Peng Wang, Fudan University Shanghai Cancer Center, China

Reviewed by:

Tongzheng Liu, Jinan University, China
Chunli Zheng, Northwest University, China

Copyright © 2020 Yang, Lin, Liu, Gan, Jin, Cui, Yan, Qian, Han and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Changpeng Han, 13918314533@163.com; Zhenyi Wang, yyyygangchangke@163.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.