Comprehensive Analysis of Differentially Expressed Long Noncoding RNA-mRNA in the Adenoma-Carcinoma Sequence of DNA Mismatch Repair Proficient Colon Cancer

Li, Wenkun; Li, Qian; Ge, Jiang; Wang, Yun; Li, Nanshan; Lao, Yueqiong; Wang, Yadan; Liu, Kuiliang; Guo, Chunmei; Lin, Wu; Jiang, Guojun; Wei, Nan; Wang, Canghai; Liu, Hong; Wu, Jing

doi:https://doi.org/10.1155/2021/9977695

Journal of Oncology

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Disclosure Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 9977695 | https://doi.org/10.1155/2021/9977695

Comprehensive Analysis of Differentially Expressed Long Noncoding RNA-mRNA in the Adenoma-Carcinoma Sequence of DNA Mismatch Repair Proficient Colon Cancer

Wenkun Li,¹Qian Li,¹Jiang Ge,¹Yun Wang,²Nanshan Li,¹Yueqiong Lao,²Yadan Wang,¹Kuiliang Liu,¹Chunmei Guo,¹Wu Lin,¹Guojun Jiang,¹and Nan Wei¹ et al.

Academic Editor: Sumanta Chatterjee

Received14 Mar 2021

Accepted08 May 2021

Published09 Jun 2021

Abstract

DNA proficient mismatch repair colon cancer (pMMR CC) is the most common subtype of sporadic CC. We aimed to investigate the role of long noncoding RNAs (lncRNAs) in pMMR CC carcinogenesis. In the present study, we conducted transcriptomic analysis of lncRNAs-mRNAs in five low-grade intraepithelial neoplasia (LGIN), five high-grade intraepithelial neoplasia (HGIN), four pMMR CC, and five normal control (NC) tissues. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment pathway, and coexpression network analyses were performed to elucidate the functions of lncRNAs and mRNAs as well as their interactions. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to validate five dysregulated lncRNAs in a large set of colon tissues. Receiver-operating characteristic (ROC) curves were employed to evaluate the performance of the candidate lncRNAs. A set of 5783 differentially expressed lncRNAs and 4483 differentially expressed mRNAs were detected among the LGIN, HGIN, pMMR CC, and NC samples. These differentially expressed lncRNAs and mRNAs were assigned to 275 significant GO terms and 179 significant KEGG enriched pathways. qRT-PCR confirmed that the expression of five selected lncRNAs (ENST00000521815, ENST00000603052, ENST00000609220, NR_026543, and ENST00000545920) were consistent with the microarray data. ROC analysis showed that four lncRNAs (ENST00000521815, ENST00000603052, ENST00000609220, and NR_026543) had larger area under the ROC curve (AUC) values compared to serum carcinoembryonic antigens, thereby distinguishing NC from pMMR CC. In conclusion, several lncRNAs play various roles in the adenoma-carcinoma sequence and may serve as potential biomarkers for the early diagnosis of pMMR CC.

1. Introduction

Colorectal cancer (CRC) is the third most commonly diagnosed cancer worldwide, comprising over 1.8 million new cases in 2018 [1]. More than 70% of CRC cases are located in the colon, which is also called colon cancer (CC) [2]. According to the latest data, CRC ranks fifth in both incidence and mortality among cancers in China [3].

CRC is a heterogeneous disease with a progressive accumulation of genetic and epigenetic alterations [4–6]. About two-thirds of CRC cases arise with unknown contributions from germline factors or a significant family history of cancer or inflammatory bowel disease, defined as sporadic CRC. Three major genetic mechanisms underlie the pathogenesis of sporadic CRC [7], namely, chromosomal instability (CIN), microsatellite instability (MSI), and the epigenomic CpG island methylator phenotype (CIMP). CIN, the first described and the most common pathway, accounts for 80%–85% of sporadic CRC [8]. CIN CRCs are nonhypermutated with DNA proficient mismatch repair (pMMR) functions [7]. The mechanism responsible for CIN is different from that responsible for MSI and CIMP. Understanding the molecular mechanisms underlying CRC is critically important for the clinical prognosis and therapeutic response. A meta-analysis of 63 eligible studies (10126 patients with CRC) showed a worse prognosis for CRC patients with CIN/pMMR [9]. Additionally, stage II CRC patients with defective mismatch repair (dMMR) or high-level MSI (MSI-H) have better prognosis but no benefits from fluorouracil chemotherapy [10, 11].

More than 90% of colorectal tumorigenesis cases follow the adenoma-carcinoma sequence (ACS), with multiple gene mutations and abnormal activation of signaling pathways [12, 13]. Mutations in the adenomatous polyposis coli (APC) gene occur early during colorectal tumorigenesis, followed by the activation of the KRAS gene and the inactivation of the TP53 gene, which are often associated with CIN [14].

Long noncoding RNAs (lncRNAs) are a class of noncoding RNAs with more than 200 nucleotides that are not translated into proteins. LncRNAs can regulate gene expression at different levels, including epigenetically, transcriptionally, and posttranscriptionally [15]. Decreased expression of the lncRNA SATB2-AS1 promotes metastasis and influences the tumor microenvironment in CRC by regulating SATB2, thereby resulting in poor survival [16]. Overexpression of the lncRNA UCA1 contributes to the immune escape of cancer cells and protects PDL1 expression from miRNA repression in gastric carcinoma, serving as a potentially novel target for immunotherapy [17]. However, the differential expression and the role of lncRNAs in ACS for CRC remain to be elucidated.

Previous studies on the lncRNA expression profiles of CC were generally performed on all cancers in the colon together including rectal cancer, thereby neglecting the heterogeneity of molecular subtypes of CC. Although often linked together, CC and rectal cancer are different when it comes to treatment. Thus far, the differential expression profile for lncRNAs in the normal colonic mucosa-pMMR adenoma-pMMR sporadic CC sequence has not been reported. In the present study, we aimed to characterize the expression profile of lncRNAs in the malignant evolution process for pMMR sporadic CC using transcriptome microarray technology. The findings provide useful candidates for the diagnosis and treatment of CC.

2. Materials and Methods

2.1. Clinical Samples

A total of 244 sporadic colonic tissues, including 63 adenomas of low-grade intraepithelial neoplasia (LGIN), 32 adenomas of high-grade intraepithelial neoplasia (HGIN), 66 pMMR CC, eight dMMR CC, and 75 normal control (NC) samples from adjacent (≥10 cm) hyperplastic or inflammatory polyps were sampled at Beijing Shijitan Hospital (Beijing, China) from 2018 to 2019. The samples were immediately placed in RNAlater solution (Cat. No. 76104, Qiagen Co, GmBH, Germany) for 24 hours at 4°C. Freshly frozen samples were stored at −80°C for RNA extraction. None of the patients received targeted therapy, chemotherapy, radiotherapy, or intervention therapy. This study was approved by the ethics committee of Beijing Shijitan Hospital (No.: 2018-59) and written informed consent was obtained from all individual participants included in the study.

2.2. Screening of pMMR Samples and Total RNA Isolation

As in our previous study [18], multiplex PCR and immunohistochemistry were used to determine the MMR status. pMMR colonic tissues were defined as tissue samples with a low frequency of MSI (MSI-L) or microsatellite stability (MSS, no marker tested). Total RNA was extracted from the tissue samples using RNAiso Plus (Code No. 9109, TAKARA, Japan). RNA integrity was confirmed with a NanoDrop ND-2000 spectrophotometer (Thermo Fisher Scientific, Rochester, NY, USA) and Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, US) and further purified with an RNeasy mini kit (Cat. # 74106, Qiagen, GmBH, Germany) and RNase-Free DNase Set (Cat. #79254, Qiagen Co, GmBH, Germany). RNA samples with integrity ≥7.0 and a 28S : 18S ratio ≥0.7 were used for further experiments.

2.3. Microarrays Analysis

The lncRNA and mRNA expression profiles in multistage colonic mucosa tissues (five NC, five LGIN, five HGIN, and four pMMR CC samples) were determined with the Agilent custom SBC Human (4180 K) lncRNA Microarray V6.0 (Product No. G4862A-074348; Agilent Technologies, Santa Clara, CA). Total RNA was amplified and labeled with the Low Input Quick Amp WT Labeling Kit (Cat. No. 5190-2943, Agilent Technologies, Santa Clara, CA, USA) following the manufacturer’s instructions. Labeled cRNAs were purified with the RNeasy mini kit (Cat. No. 74106, Qiagen, GmBH, Germany). Each slide was hybridized with 1.65 μg Cy3-labeled cRNA using a Gene Expression Hybridization Kit (Cat. No. #5188-5242, Agilent Technologies, Santa Clara, CA, US) in a hybridization oven (Cat. No. G2545 A, Agilent Technologies, Santa Clara, CA, USA) for 17 h. The slides were washed in staining dishes (Cat. No. 121, Thermo Shandon, Waltham, MA, US) with the Gene Expression Wash Buffer Kit (Cat. No. 5188-5327, Agilent Technologies, Santa Clara, CA, US) following the manufacturer’s instructions and scanned with an Agilent Microarray Scanner (Cat. No. G2565CA, Agilent Technologies, Santa Clara, CA, USA) with default settings (Dye channel: Green, Scan resolution = 3 μm, PMT 100%, 20 bit). The data were extracted with Feature Extraction Software 10.7 (Agilent Technologies, Santa Clara, CA, USA).

2.4. Data Processing and Identification of DEGs

2.4.1. Series Test of Cluster (STC) Analysis

Raw data were normalized with the Quantile algorithm in GeneSpring Software 11.0 (Agilent Technologies, Santa Clara, CA, USA). Cluster analysis was performed to systematically and intuitively display the relatedness between samples. Differential genes expressions (DEGs) among the four groups were filtered with the F test in the random variance model (RVM) using values < 0.05 and Benjamini-Hochberg false discovery rate (FDR) < 0.05 as the significant cutoff criteria [19]. STC analysis was employed to identify the most probable set of clusters responsible for the observed ACS series by profiling the dynamic nature of the temporal gene expression. NC, LGIN, HGIN, and pMMR CC were set as different points to identify significant trending models related to multistage colonic mucosa tissue and their associated DEGs by STC analysis. Clusters with values < 0.05/26 were considered statistically significant.

2.4.2. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Pathway Analysis

Gene Ontology (GO) analysis was performed to annotate DEGs in three categories: molecular functions, biological processes, and cellular components. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was performed to identify the functional pathways and pathways for molecular interactions, reactions, and relation networks. The statistically significant threshold was set as value < 0.01 and FDR <0.05.

2.5. Construction of Coexpression and Transcription Factor (TF) Regulatory Network

To identify coexpressed lncRNA-mRNA pairs and hub DEGS, the lncRNA-mRNA coexpression network was constructed based on the gene expression value and Pearson's correlation coefficient (PCC) between the differentially expressed lncRNAs and mRNAs. |PCC| ≥0.8 and value < 0.01 were considered as statistically significantly relevant. In order to distinguish the TF regulatory network in the stepwise process encompassing CC progression, the interactions were predicted starting from a TF by searching the conserved TF binding sites within a putative promoter area 2000 bp upstream and 500 bp downstream of the transcriptional initiation site of coexpressed genes.

2.6. Validation by Quantitative Real-Time PCR (qRT-PCR)

Based on types of lncRNAs, aberrant expression (fold change), the changes in coexpression degrees, and the results of CPC, CNCI, PFAM, and PhyloCSF scores, 5 lncRNAs were selected to validate the results of microarray. The expression levels of the lncRNAs were examined by qRT-PCR in 225 colorectal tissue samples. qRT-PCR was performed with a QuantStudio^TM 6 Flex Real-Time PCR System (Thermo Fisher, Rochester, NY, USA). Primers are listed in Additional file 1 and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as the endogenous control. Each qRT-PCR reaction (in 10 μL) contained 5 μL 2X SYBR Green PCR Mix (ABI, USA), 0.4 μM of each of the forward and reverse primers, and 5 ng of template cDNA. The PCR program was as follows: initial denaturing at 95°C for 5 min followed by 40 cycles at 94°C for 30 s and 60°C for 1 min and fluorescence acquisition at 60°C for 1 min. PCR amplifications were performed in triplicate for each sample.

2.7. Statistical Analysis

Computer-based calculations were conducted using SPSS version 20.0 (SPSS Inc., Chicago, IL, USA) and MedCalc version 19.1 (MedCalc Software Bvba, Ostend, Belgium). The differences in the expression of selected lncRNAs among multiple groups were compared using ANOVA. Further pairwise comparisons were performed with the least significance difference (LSD) test. The differences in the expression of lncRNAs between the pMMR CC group and the dMMR CC group were compared with the nonparametric Mann–Whitney U test. The relationship between lncRNAs and their target coding genes were determined by Spearman correlation analysis. Receiver-operating characteristic (ROC) curves and the area under the ROC curve (AUC) were employed to evaluate the diagnostic accuracy of candidate lncRNAs serum carcinoembryonic antigen (CEA), obtained from Hospital Information System. The data were presented as the mean with standard deviation (SD) for normally distributed data and the median with the interquartile range for skewed data. |Fold change| ≥5 was used as the significant threshold to screen differentially expressed lncRNA and mRNA. All values were two-sided, and FDR was calculated for multiple comparisons. Differences with were considered as statistically significant.

3. Results

3.1. Basic Characteristics of the Microarray Data

The total RNA extracted from pMMR colonic tissue samples, consisting of five NC, five LGIN, five HGIN, and four pMMR CC samples, were hybridized to a transcriptome microarray. The quality control results for RNA/microarray are summarized in Additional file 1 and Additional file 2, respectively. Cluster analysis of normalized data indicated that the samples from each group were well separated when all DEGs were considered (Figure 1), suggesting that the microarray data were reliable for further bioinformatics analysis.

3.2. Dysregulated lncRNAs and mRNAs in the Colonic Adenoma-Carcinoma Sequence

A total of 5783 lncRNAs and 4483 mRNAs were differentially expressed among the NC, LGIN, HGIN, and pMMR CC samples. In order to generate dynamic and significant trending models of DEGs related to colon mucosa malignant transformation, STC analysis was used to further identify the clusters of DEGs with similar expression patterns. As shown in Figure 2(a), the differential expression profile consisted of 26 differentially expressed lncRNA clusters, among which 15 clusters (profile no. 0, 1, 2, 3, 4, 5, 9, 12, 17, 20, 21, 22, 23, 24, and 25) showed significant expression trends (). Moreover, the lncRNA expression levels in clusters 21, 22, 24, and 25 were stable or gradually elevated, while the lncRNA expression in clusters 0, 1, 3, 4, 9, and 12 showed a decreasing trend. Cluster 23 contained 207 lncRNAs that gradually increased from NC to adenoma and were suppressed in CC (Figure 2(b)). However, lncRNA expression in cluster 3 was opposite that in cluster 23. As shown in Figure 2(a), a total of nine mRNA expression clusters, specifically clusters 2, 4, 5, 13, 20, 21, 22, 24, and 25, showed significant expression trends (), among which clusters 21, 22, 24, and 25 showed a stable or gradual elevation. The results suggested that adenoma is an intermediate step from normal tissue to CC, and certain lncRNAs may play important roles during the dynamic process in colonic mucosal protruding lesions.

(a)

(b)

Figure 2

STC analysis of the dysregulated DEGs in multistage colonic mucosa tissues. (a) Twenty-six clusters for the expression pattern of DEGs, 15 lncRNA expression patterns, and nine mRNA expression patterns showed significant values (colored boxes, /26). Each box represents a model expression profile, the lower number in the profile box is the value, and the upper number is the model profile number. (b) Cluster 25 (173 lncRNAs) gradually increased, cluster 0 (87 lncRNAs) showed a decreasing trend, cluster 23 contained 207 lncRNAs that gradually increased from NC to adenoma and then decreased in CC, and cluster 25 (146 mRNAs) gradually increased. Horizontal axis represents different stages of colonic mucosal protruding lesions and the vertical axis shows the time series for the lncRNA/mRNA expression levels. r normalized with log transformation.

3.3. Functional and Pathway Enrichment Analysis

DEGs from significant expression trends were subjected to GO and KEGG pathway analyses to identify their potential functions and mechanisms in the dynamic process for colonic mucosal protruding lesions (Figure 3).

(a)

(b)

The top 20 GO terms included multicellular organismal development, cell division, cell adhesion, cell cycle, cell differentiation, mitosis, proteolysis, and angiogenesis. KEGG enrichment analysis revealed that the DEGs were mainly involved in several signaling pathways, including metabolism, cancer, cell cycle, PI3K-Akt signaling, and transcriptional misregulation in cancer (Figure 3(b)).

3.4. Construction of the lncRNA-mRNA Coexpression networks

To discover the significant molecular mechanisms for the lncRNAs associated with tumorigenesis in colonic ACS, lncRNAs from elevated clusters [20–23] or reduced clusters (0, 1, 3, 4, 9, and 12) and mRNAs with significant GO terms and KEGG enrichment pathways were selected to construct elevated-lncRNA-mRNA coexpression networks (coexpression network E) and decreased-lncRNAs-mRNA coexpression networks (coexpression network D), respectively. The top 10 lncRNA/mRNA in terms by degree are listed in Tables 1 and 2. Coexpression network E contained 714 nodes and 1711 edges (Figure 4(a)), among which the hub nodes with the highest degrees were NOTCH4 (mRNA, degree = 29), INHBA (mRNA, degree = 26), TSTA3 (mRNA, degree = 29), lnc-NPRL3-1:1 (lncRNA, degree = 22), NR_024431 (lncRNA, degree = 21), and ENST00000564984 (lncRNA, degree = 17). Coexpression network D contained 598 nodes and 1908 edges (Figure 4(b)) and the DEGs with the highest degrees were AQP8 (mRNA, degree = 47), STRADB (mRNA, degree = 47), ENST00000435912 (lncRNA, degree = 15), NR_024431 (lncRNA, degree = 21), and ENST00000564984 (lncRNA, degree = 17) (Tables 1 and 3).

(a)

(b)

(c)

(d)

Figure 4

lncRNA-mRNA coexpression networks. (a-b) LncRNAs with elevated or decreased profiles and mRNAs for significant terms from GO and KEGG enrichment analysis were, respectively, selected to construct coexpression networks E and D in colonic ACS. (c-d) LncRNAs with significantly elevated/decreased model profiles and angiogenesis-related GO and KEGG items were, respectively, selected to construct A-coexpression networks E and D in colonic ACS. Nodes represent DEGs (green for mRNAs and yellow for lncRNAs). Lines represent interactions and node size represents degree value.

Angiogenesis is a crucial step in tumor growth and progression. The angiogenic switch has been observed at the adenoma stage in ACS [24]. We constructed angiogenesis-related lncRNA-mRNA coexpression networks (A-coexpression networks E and D) with DEGs from significantly elevated/decreased clusters as well as GO terms and KEGG pathways related to angiogenesis. A-coexpression network E consisted of 649 nodes and 2299 edges (Figure 4(c)) and the top hub genes were NOTCH4 (mRNA, degree = 35), MAD2L1 (mRNA, degree = 34), E2F7 (mRNA, degree = 34), lnc-ZBTB20-2:1 (lncRNA, degree = 22), lnc-CENPH-2:1 (lncRNA, degree = 20), and lnc-PAQR4-2:1 (lncRNA, degree = 19) (Tables 2 and 3). The top hub genes in A-coexpression network D (Figure 4(d)) were CAMK2B (mRNA, degree = 61), STRADB (mRNA, degree = 60), LTK (mRNA, degree = 57), ENST00000513255 (lncRNA, degree = 12), NR_024605 (lncRNA, degree = 11), and NONHSAT057082 (lncRNA, degree = 11) (Tables 2 and 3). All the DEGs primarily participated in angiogenesis as the key genes in each network.

3.5. Construction of TF Regulatory Networks

To investigate the mechanism of gene regulation at the transcriptional level, we constructed two interaction networks (TF-DEG network (Figure 5(a)) and TF-angiogenesis-DEG network (Figure 5(b))) between TFs and DEGs that originated from the coexpression network or A-coexpression network. The TF-DEG network contained 99 regulation models and 402 nodes involving 85 TFs (Figure 5(a)), among which three TFs, namely, LUN-1 (TF, degree = 218), Tel-2 (TF, degree = 22), and 1-Oct (TF, degree = 20), were regulated over 20 nodes (Table 4). LUN-1 had the highest degree and might play more important roles in the TF regulatory networks. The lncRNAs regulated by the most TFs were NR_103548, lnc-ARRDC3-1:16, ENST00000513626, and ENST00000516496 (regulated by six TFs). In the TF-angiogenesis-DEG network, due to the important role of LUN-1 (TF, degree = 73) and Tel-2 (TF, degree = 18) in the progression of CC, a subnetwork was constructed showing their effects on angiogenesis functions (Figure 5(b) and Table 4).

(a)

(b)

3.6. qRT-PCR Verification for the Candidate Genes

According to the probe signal value, type, DEG fold change, degree value, and protein coding ability, five LncRNAs were selected from the networks above for validation in the expanded clinical samples (70 NC, 58 LIGN, 27 HIGN, 62 pMMR CC, and 8 dMMRCC samples). LncRNAs ENST00000545920, ENST00000521815, ENST00000609220, and ENST00000603052 presented a sequentially increasing trend in expression with tumor progression (Figures 6(a)–6(d)), whereas NR_026543 showed a sequentially decreasing trend (Figure 6(e)). No statistically significant differences were observed in the expressions of NR_026543 and ENST00000521815 among LGIN, HGIN, and pMMR CC. However, the expression levels for each group were significantly elevated or decreased compared with those of the NC samples, suggesting that these genes could be used to distinguish normal from lesion tissue. Furthermore, the expressions of ENST00000603052, ENST00000521815, and ENST00000609220 were significantly different in pMMR CC and dMMR CC tissue samples (Figures 6(f)–6(j)), indicating that these three genes could distinguish pMMR CC from dMMR CC.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

Figure 6

qRT-PCR verification of the candidate genes. ENST00000545920 (a), ENST00000521815 (b), ENST00000609220 (c), and ENST00000603052 (d) presented a sequentially ascending trend in expression with tumor progression, whereas NR_026543 (e) showed a sequentially decreasing trend. The expression of ENST00000521815 (f), ENST00000603052 (h), and ENST00000609220 (i) showed statistically significant differences between pMMR CC and dMMR CC. However, the expression of ENST00000609220 (g) and NR_026543 (j) did not show statistically significant differences between two groups.

3.7. Diagnostic Efficacy of Selected lncRNAs

The diagnostic efficacy of the above selected five lncRNAs was evaluated as potential biomarkers for pMMR CC diagnosis compared to carcinoembryonic antigen (CEA), a traditional serum tumor biomarker. The specificities of ENST00000521815, ENST00000603052, ENST00000609220, NR_026543, ENST00000545920, and CEA were 0.981, 1, 0.962, 0.932, 0.981, and 0.904, and the sensitivities were 0.931, 0.864, 0.746, 0.846, 0.576, and 0.593, respectively. ROC analysis showed significantly higher AUC for lncRNAs ENST00000521815, ENST00000603052, ENST00000609220, and NR_026543 compared with CEA (0.785) (, 0.0002, 0.0028, and 0.0011, respectively) (Figures 7(a)–7(d)). No statistically significant difference in the AUC was observed between ENST00000545920 (0.794) and CEA () (Figure 7(e)).

(a)

(b)

(c)

(d)

(e)

Figure 7

Diagnostic efficacy of selected lncRNAs. Compared with CEA, ENST00000521815 (a), ENST00000603052 (b), ENST00000609220 (c), and (d) NR_026543 had significantly larger AUCs (AUCs 0.957, 0.958, 0.933, and 0.949 versus 0.785, , 0.0002, 0.0028, and 0.0011, respectively). However, there were no statistically significant differences between the AUC of ENST00000545920 (0.794) and that of CEA () (e). The ROC curve for CEA is shown in blue, while the ROC curves for the validation of the lncRNAs are shown in green.

4. Discussion

Among the human transcriptome, more than 90% are noncoding RNAs (ncRNAs), including microRNAs (miRNAs), lncRNAs, and circular RNAs (circRNAs) [20]. As a relatively new field, accumulated studies have uncovered the crucial roles of lncRNAs in oncogenicity [21, 25]. With the improvement of microarray and next-generation sequencing technology, large-scale aberrantly expressed lncRNAs have been discovered in multiple cancers, including breast cancer [22], gastric cancer [23], CRC [26], lung cancer [27], and pancreatic cancer [28].

A large proportion of colon adenomas are able to transform into adenocarcinoma through a dynamic process with the accumulation of multiple mutations and preternatural activation of signaling pathways [13]. Altered lncRNA expression patterns for microsatellite-stable colorectal cancer indicated that several lncRNAs were continuously upregulated/downregulated during the adenoma-adenocarcinoma sequence [29]. Unlike previous studies that focused on one or more lncRNAs in CRC [30–32], we first conducted a systematic and comprehensive analysis using microarrays to reveal the differentially expressed profiles of lncRNAs and mRNAs in the malignant evolutionary process for pMMR CC. Thousands of dysregulated transcripts were identified in the ACS, including 5783 lncRNAs and 4483 mRNAs. Subsequently, STC analysis revealed that adenoma is an intermediate step from normal tissue to CC and certain lncRNAs participate in this continuous process. As the dysregulation can be detected already in adenomas, these lncRNAs can be used as biomarkers for the screening and early detection of CC.

GO annotation and KEGG pathway analysis were performed to determine the functions underlying the differentially expressed lncRNAs in ACS. GO annotation suggested that aberrant transcripts were predominantly enriched in multicellular organismal development, cell division, cell adhesion, cell cycle, cell differentiation, mitosis, proteolysis, and angiogenesis. KEGG enrichment analysis revealed that the dysregulated transcripts mainly participated in metabolic pathways, pathways in cancer, cell cycle, PI3K-Akt signaling pathway, and purine metabolism. Cell cycle involvement was identified in both analyses. The uncontrolled tumor cell cycle is an essential feature of cancer and is caused by aberrant cell cycle proteins including cyclin-dependent kinases, Aurora kinases, and Polo-like kinases. Exploiting key molecules in the cell cycle may provide a new perspective for cancer therapy [33].

To further explore the underlying functions of lncRNAs in ACS, coexpression of lncRNAs with their associated coding genes was performed based on the degree of correlation. Among the top 10 elevated and decreased lncRNAs according to the degree, NR_029374 with 12 degrees was upregulated, consistent with the results obtained by Sun et al. [34] who reported that NR_029374 was upregulated in CC and correlated with poor overall survival. NR_029374 promoted the migration, invasion, and metastasis of CC through the miR-30-5p/SOX9 axis and this finding was supported by previous results [35]. NR_029374 was upregulated in hepatocellular carcinoma (HCC), markedly promoting HCC proliferation, migration, and angiogenesis [35].

Angiogenesis plays a key role in tumorigenesis and development by delivering nutrients and evacuating metabolic wastes [36]. Recently, several lncRNAs and mRNAs were demonstrated to participate in tumor angiogenesis. Upregulation of the lncRNA FLANC promoted angiogenesis by upregulating and prolonging the half-life of phosphorylated STAT3/VEGFA in CRC [37]. In the present study, we discovered that NOTCH4 (mRNA, degree = 35), a member of the Notch family of transmembrane receptors, had the highest degree. Wu et al. reported the association of overexpressed NOTCH4 with CRC survival [38]. The NOTCH4 expression level was especially higher in non-small-cell lung cancer tissues than in the NC tissues, related to the tumor size and TNM stage [39]. Active NOTCH4 might inhibit endothelial sprouting in vitro and vivo [40]. Further study is required to determine the expression level of NOTCH4 in a large population of pMMR CC patients and the underlying interactions with lncRNAs in angiogenesis. Several lncRNAs have demonstrated interactions with TFs in tumorigenesis and progression. For instance, lncRNA SNHG15 obstructed the ubiquitination and degradation of Slug, a fast-turnover TF critical for controlling cancer cell invasion and metastasis, thereby promoting CC progression [41]. The transcription of lncRNA LINC01503 was activated by TF TP63, resulting in shorter survival times for patients with esophageal squamous cell carcinoma [42]. In the current study, the lncRNA-TF network analysis revealed that LUN-1 had the highest degree in both the TF-DEG network and the TF-angiogenesis-DEG network.

qRT-PCR analysis showed the upregulation of ENST00000545920, ENST00000521815, ENST00000609220, and ENST00000603052 and the downregulation of NR_026543 with tumor progression. Furthermore, ENST00000603052, ENST00000521815, and ENST00000609220 can be used to identify pMMR CC and dMMR CC. Previous studies revealed the upregulation of ENST00000545920, also known as lnc-SNHG1, in CRC [41], non-small-cell lung cancer [43], gastric cancer [44], and glioma [45] to promote cancer cell growth by interacting with EZH2 in the nucleus and miR-154-5p in the cytoplasm [46]. Similar results have been observed in other tumors, including ENST00000521815, also known as CASC19, which possesses an oncogenic function through targeting miR-140-5p/CEMIP in CRC progression [47]. The downregulation of NR_026543 was associated with liver metastasis and poor prognosis for CC by binding with the miR-203 promoter [48]. ENST00000609220 was found to be highly expressed in CRC tissue, promoting CRC growth and metastasis through the miR-206/YAP1 axis [49]. To date, there is little knowledge on the function of ENST00000603052 in cancer.

Ultimately, the diagnostic efficacy of the identified lncRNAs was evaluated against serum CEA. Most of the selected lncRNAs achieved a slightly higher AUC in distinguishing CC from normal tissue. Further study is needed to validate these lncRNAs as biomarkers for the early diagnosis of CC by measuring their expression in serum samples to assess the consistency in tissue samples. Monitoring expression changes in these genes may also provide a new strategy to assess disease progression.

There are certain limitations in the present study. First, the sample size was relatively small, thereby limiting result reliability. We increased the number of clinical samples to confirm the function of the lncRNAs. Second, we only constructed expression profiles for colon tissue. Additional studies are required to evaluate tissues from other tumors. Last, compared with high-throughput sequencing, microarray is incapable of identifying novel lncRNAs.

5. Conclusions

In conclusion, we verified a series of differentially expressed lncRNAs and mRNAs in a normal mucosa-adenoma-pMMR CC sequence in the current study. The potential roles of these RNAs were predicted through bioinformatics analyses. The findings may provide novel insights for the diagnosis and therapeutic strategy of pMMR CC. Further studies are required to provide robust validation evidence for the functions of lncRNAs in CC.

Data Availability

All the data generated or analyzed during this study are included in this published article and its supplementary information files.

Disclosure

This manuscript was sent to Research Square as a preprint (https://www.researchsquare.com/article/rs-47627/v1.pdf); however, it has not been published elsewhere in whole or in part. The funding bodies had no role in the study design, data collection, data interpretation, or writing of this manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Wenkun Li and Qian Li contributed equally to this work.

Acknowledgments

This work was supported by the Digestive Medical Coordinated Development Center of Beijing Hospitals Authority under Grant no. XXZ015, Medical and Health Public Foundation of Beijing under Grant no. YWJKJJHKYJJ-B17262-067, and the Science and Technology Development Project of China State Railway Group under Grant no. N2019Z004.

Supplementary Materials

Additional file 1. Quality control results of sample RNA. Additional file 2. Quality control results of microarrays. (Supplementary Materials)

References

F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
C. E. Bailey, C.-Y. Hu, Y. N. You et al., “Increasing disparities in the age-related incidences of colon and rectal cancers in the United States, 1975-2010,” JAMA Surgery, vol. 150, no. 1, pp. 17–22, 2015.
View at: Publisher Site | Google Scholar
W. Chen, R. Zheng, P. D. Baade et al., “Cancer statistics in China, 2015,” CA: A Cancer Journal for Clinicians, vol. 66, no. 2, pp. 115–132, 2016.
View at: Publisher Site | Google Scholar
R. Dienstmann, L. Vermeulen, J. Guinney, S. Kopetz, S. Tejpar, and J. Tabernero, “Consensus molecular subtypes and the evolution of precision medicine in colorectal cancer,” Nature Reviews Cancer, vol. 17, no. 2, pp. 79–92, 2017.
View at: Publisher Site | Google Scholar
E. R. Fearon, “Molecular genetics of colorectal cancer,” Annual Review of Pathology: Mechanisms of Disease, vol. 6, no. 1, pp. 479–507, 2011.
View at: Publisher Site | Google Scholar
E. R. Fearon and B. Vogelstein, “A genetic model for colorectal tumorigenesis,” Cell, vol. 61, no. 5, pp. 759–767, 1990.
View at: Publisher Site | Google Scholar
J. M. Carethers and B. H. Jung, “Genetics and genetic biomarkers in sporadic colorectal cancer,” Gastroenterology, vol. 149, no. 5, Article ID e1173, pp. 1177–1190, 2015.
View at: Publisher Site | Google Scholar
W. M. Grady and J. M. Carethers, “Genomic and epigenetic instability in colorectal cancer pathogenesis,” Gastroenterology, vol. 135, no. 4, pp. 1079–1099, 2008.
View at: Publisher Site | Google Scholar
A. Walther, R. Houlston, and I. Tomlinson, “Association between chromosomal instability and prognosis in colorectal cancer: a meta-analysis,” Gut, vol. 57, no. 7, pp. 941–950, 2008.
View at: Publisher Site | Google Scholar
A. B. Benson 3rd, A. P. Venook, L Cederquist et al., “Colon cancer, version 1.2017, NCCN clinical practice guidelines in oncology,” Journal of the National Comprehensive Cancer Network: JNCCN, vol. 15, no. 3, pp. 370–398, 2017.
View at: Publisher Site | Google Scholar
D. J. Sargent, S. Marsoni, G. Monges et al., “Defective mismatch repair as a predictive marker for lack of efficacy of fluorouracil-based adjuvant therapy in colon cancer,” Journal of Clinical Oncology, vol. 28, no. 20, pp. 3219–3226, 2010.
View at: Publisher Site | Google Scholar
A. Leslie, F. A. Carey, N. R. Pratt, and R. J. C. Steele, “The colorectal adenoma-carcinoma sequence,” British Journal of Surgery, vol. 89, no. 7, pp. 845–860, 2002.
View at: Publisher Site | Google Scholar
Z. Wu, Z. Liu, W. Ge et al., “Analysis of potential genes and pathways associated with the colorectal normal mucosa-adenoma-carcinoma sequence,” Cancer Medicine, vol. 7, no. 6, pp. 2555–2566, 2018.
View at: Publisher Site | Google Scholar
H. Brenner, M. Kloor, and C. P. Pox, “Colorectal cancer,” The Lancet, vol. 383, no. 9927, pp. 1490–1502, 2014.
View at: Publisher Site | Google Scholar
Z. Liu, L. Ren, J. Tian, N. Liu, Y. Hu, and P. Zhang, “Comprehensive analysis of long noncoding RNAs and messenger RNAs expression profiles in patients with marjolin ulcer,” Medical Science Monitor: International Medical Journal of Experimental and Clinical Research, vol. 24, pp. 7828–7840, 2018.
View at: Publisher Site | Google Scholar
M. Xu, X. Xu, B. Pan et al., “LncRNA SATB2-AS1 inhibits tumor metastasis and affects the tumor immune cell microenvironment in colorectal cancer by regulating SATB2,” Molecular Cancer, vol. 18, no. 1, p. 135, 2019.
View at: Publisher Site | Google Scholar
C. J. Wang, C. C. Zhu, J. Xu et al., “The lncRNA UCA1 promotes proliferation, migration, immune escape and inhibits apoptosis in gastric cancer by sponging anti-tumor miRNAs,” Molecular Cancer, vol. 18, no. 1, p. 115, 2019.
View at: Publisher Site | Google Scholar
Q. Li, N. Li, Y. Lao et al., “Variable levels of long noncoding RNA expression in DNA mismatch repair-proficient early-stage colon cancer,” Digestive Diseases and Sciences, vol. 62, no. 5, pp. 1235–1245, 2017.
View at: Publisher Site | Google Scholar
B. Yoav and H. Yosef, “Controlling the false discovery rate: a practical and powerful approach to multiple testing,” Journal of the Royal Statistical Society: Series B (Methodological), vol. 57, pp. 289–300, 1995.
View at: Publisher Site | Google Scholar
X. Yu, H. Zheng, G. Tse, M. T. Chan, and W. K. Wu, “Long non-coding RNAs in melanoma,” Cell Proliferation, vol. 51, no. 4, Article ID e12457, 2018.
View at: Publisher Site | Google Scholar
Y. Fang and M. J. Fullwood, “Roles, functions, and mechanisms of long non-coding RNAs in cancer,” Genomics, Proteomics and Bioinformatics, vol. 14, no. 1, pp. 42–54, 2016.
View at: Publisher Site | Google Scholar
Y.-Z. Jiang, Y.-R. Liu, X.-E. Xu et al., “Transcriptome analysis of triple-negative breast cancer reveals an integrated mRNA-lncRNA signature with predictive and prognostic value,” Cancer Research, vol. 76, no. 8, pp. 2105–2114, 2016.
View at: Publisher Site | Google Scholar
K. Zhang, H. Shi, H. Xi et al., “Genome-Wide lncRNA microarray profiling identifies novel circulating lncRNAs for detection of gastric cancer,” Theranostics, vol. 7, no. 1, pp. 213–227, 2017.
View at: Publisher Site | Google Scholar
C. A. Staton, A. S. A. Chetwood, I. C. Cameron, S. S. Cross, N. J. Brown, and M. W. R. Reed, “The angiogenic switch occurs at the adenoma stage of the adenoma carcinoma sequence in colorectal cancer,” Gut, vol. 56, no. 10, pp. 1426–1432, 2007.
View at: Publisher Site | Google Scholar
A. Sanchez Calle, Y. Kawamura, Y. Yamamoto, F. Takeshita, and T. Ochiya, “Emerging roles of long non-coding RNA in cancer,” Cancer Science, vol. 109, no. 7, pp. 2093–2100, 2018.
View at: Publisher Site | Google Scholar
A. Yamada, P. Yu, W. Lin, Y. Okugawa, C. R. Boland, and A. Goel, “A RNA-Sequencing approach for the identification of novel long non-coding RNA biomarkers in colorectal cancer,” Scientific Reports, vol. 8, no. 1, p. 575, 2018.
View at: Publisher Site | Google Scholar
Z. Peng, J. Wang, B. Shan et al., “Genome-wide analyses of long noncoding RNA expression profiles in lung adenocarcinoma,” Scientific Reports, vol. 7, no. 1, Article ID 15331, 2017.
View at: Publisher Site | Google Scholar
S. Müller, S. Raulefs, P. Bruns et al., “Next-generation sequencing reveals novel differentially regulated mRNAs, lncRNAs, miRNAs, sdRNAs and a piRNA in pancreatic cancer,” Molecular Cancer, vol. 14, p. 94, 2015.
View at: Publisher Site | Google Scholar
A. N. Bartley, H. Yao, B. A. Barkoh et al., “Complex patterns of altered MicroRNA expression during the adenoma-adenocarcinoma sequence for microsatellite-stable colorectal cancer,” Clinical Cancer Research, vol. 17, no. 23, pp. 7283–7293, 2011.
View at: Publisher Site | Google Scholar
Z. Lan, X. Yao, K. Sun, A. Li, S. Liu, and X. Wang, “The interaction between lncRNA SNHG6 and hnRNPA1 contributes to the growth of colorectal cancer by enhancing aerobic glycolysis through the regulation of alternative splicing of PKM,” Frontiers in Oncology, vol. 10, p. 363, 2020.
View at: Publisher Site | Google Scholar
H. Wu, M. Wei, X. Jiang et al., “lncRNA PVT1 promotes tumorigenesis of colorectal cancer by stabilizing miR-16-5p and interacting with the VEGFA/VEGFR1/AKT Axis,” Molecular Therapy—Nucleic Acids, vol. 20, pp. 438–450, 2020.
View at: Publisher Site | Google Scholar
X. Yao, Z. Lan, Q. Lai, A. Li, S. Liu, and X. Wang, “LncRNA SNHG6 plays an oncogenic role in colorectal cancer and can be used as a prognostic biomarker for solid tumors,” Journal of Cellular Physiology, vol. 235, pp. 7620–7634, 2020.
View at: Google Scholar
T. Otto and P. Sicinski, “Cell cycle proteins as promising targets in cancer therapy,” Nature Reviews Cancer, vol. 17, no. 2, pp. 93–115, 2017.
View at: Publisher Site | Google Scholar
T. Sun, Z. Liu, R. Zhang et al., “Long non-coding RNA LEF1-AS1 promotes migration, invasion and metastasis of colon cancer cells through miR-30-5p/SOX9 Axis,” OncoTargets and Therapy, vol. 13, pp. 2957–2972, 2020.
View at: Publisher Site | Google Scholar
H. Dong, P. Jian, M. Yu, and L. Wang, “Silencing of long noncoding RNA LEF1-AS1 prevents the progression of hepatocellular carcinoma via the crosstalk with microRNA-136-5p/WNK1,” Journal of Cellular Physiology, vol. 235, 2020.
View at: Publisher Site | Google Scholar
J. Folkman, “Role of angiogenesis in tumor growth and metastasis,” Seminars in Oncology, vol. 29, no. 6, pp. 15–18, 2002.
View at: Publisher Site | Google Scholar
M. Pichler, C. Rodriguez-Aguayo, S. Y. Nam et al., “Therapeutic potential of FLANC, a novel primate-specific long non-coding RNA in colorectal cancer,” Gut, vol. 69, pp. 1818–1831, 2020.
View at: Google Scholar
G. Wu, Z. Chen, J. Li et al., “NOTCH4 is a novel prognostic marker that correlates with colorectal cancer progression and prognosis,” Journal of Cancer, vol. 9, no. 13, pp. 2374–2379, 2018.
View at: Publisher Site | Google Scholar
Y. Wang, R. Yang, X. Wang et al., “Evaluation of the correlation of vasculogenic mimicry, Notch4, DLL4, and KAI1/CD82 in the prediction of metastasis and prognosis in non-small cell lung cancer,” Medicine, vol. 97, no. 52, Article ID e13817, 2018.
View at: Publisher Site | Google Scholar
K. G. Leong, X. Hu, L. Li et al., “Activated Notch4 inhibits angiogenesis: role of β1-integrin activation,” Molecular and Cellular Biology, vol. 22, no. 8, pp. 2830–2841, 2002.
View at: Publisher Site | Google Scholar
H. Jiang, T. Li, Y. Qu et al., “Long non-coding RNA SNHG15 interacts with and stabilizes transcription factor Slug and promotes colon cancer progression,” Cancer Letters, vol. 425, pp. 78–87, 2018.
View at: Publisher Site | Google Scholar
J.-J. Xie, Y.-Y. Jiang, Y. Jiang et al., “Super-enhancer-driven long non-coding RNA LINC01503, regulated by TP63, is over-expressed and oncogenic in squamous cell carcinoma,” Gastroenterology, vol. 154, no. 8, Article ID e2131, pp. 2137–2151, 2018.
View at: Publisher Site | Google Scholar
Q. Lu, S. Shan, Y. Li, D. Zhu, W. Jin, and T. Ren, “Long noncoding RNASNHG1promotes non‐small cell lung cancer progression by up‐regulatingMTDH viasponging miR‐145-5p,” The FASEB Journal, vol. 32, no. 7, pp. 3957–3967, 2018.
View at: Publisher Site | Google Scholar
W. Guo, J. Huang, P. Lei, L. Guo, and X. Li, “LncRNA SNHG1 promoted HGC-27 cell growth and migration via the miR-140/ADAM10 axis,” International Journal of Biological Macromolecules, vol. 122, pp. 817–823, 2019.
View at: Publisher Site | Google Scholar
H. Li, Y. Xue, J. Ma et al., “SNHG1 promotes malignant biological behaviors of glioma cells via microRNA-154-5p/miR-376b-3p- FOXP2- KDM5B participating positive feedback loop,” Journal of Experimental and Clinical Cancer Research : CR, vol. 38, no. 1, p. 59, 2019.
View at: Publisher Site | Google Scholar
M. Xu, X. Chen, K. Lin et al., “The long noncoding RNA SNHG1 regulates colorectal cancer cell growth through interactions with EZH2 and miR-154-5p,” Molecular Cancer, vol. 17, no. 1, p. 141, 2018.
View at: Publisher Site | Google Scholar
X.-D. Wang, J. Lu, Y.-S. Lin, C. Gao, and F. Qi, “Functional role of long non-coding RNA CASC19/miR-140-5p/CEMIP axis in colorectal cancer progressionin vitro,” World Journal of Gastroenterology, vol. 25, no. 14, pp. 1697–1714, 2019.
View at: Publisher Site | Google Scholar
L. Wang, Z. Wei, K. Wu et al., “Long noncoding RNA B3GALT5-AS1 suppresses colon cancer liver metastasis via repressing microRNA-203,” Aging, vol. 10, no. 12, pp. 3662–3682, 2018.
View at: Publisher Site | Google Scholar
X. Dong, Z. Yang, H. Yang, D. Li, and X. Qiu, “Long non-coding RNA mir4435-2HG promotes colorectal cancer proliferation and metastasis through miR-206/YAP1 Axis,” Frontiers in Oncology, vol. 10, p. 160, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Wenkun Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

263

Downloads

898

Citations