Hypoxia-induced PTTG3P contributes to colorectal cancer glycolysis and M2 phenotype of macrophage

Abstract Long noncoding RNAs (lncRNAs) play critical factors in tumor progression and are ectopically expressed in malignant tumors. Until now, lncRNA pituitary tumor-transforming 3, pseudogene (PTTG3P) biological function in colorectal cancer (CRC) further needs to be clarified. qRT-PCR was used to measure the PTTG3P level and CCK-8, glucose uptake, lactate assay, adenosine triphosphate (ATP) assay, extracellular acidification rate (ECAR) assay, and xenograft mice model were adopted to evaluate the glycolysis and proliferation, and macrophage polarization were determined in CRC cells. Xenograft experiments were utilized to analyze tumor growth. Ectopic expression of PTTG3P was involved in CRC and related to dismal prognosis. Through gain- and loss-of-function approaches, PTTG3P enhanced cell proliferation and glycolysis through YAP1. Further, LDHA knockdown or glycolysis inhibitor (2-deoxyglucose (2-DG), 3-BG) recovered from PTTG3P-induced proliferation. And PTTG3P overexpression could facilitate M2 polarization of macrophages. Silenced PTTG3P decreased the level of inflammatory cytokines TNF-α, IL-1β and IL-6, and low PTTG3P expression related with CD8+ T, NK, and TFH cell infiltration. Besides, hypoxia-inducible factor-1α (HIF1A) could increase PTTG3P expression by binding to the PTTG3P promoter region. Hypoxia-induced PTTG3P contributes to glycolysis and M2 phenotype of macrophage, which proposes a novel approach for clinical treatment.


Introduction
Colorectal cancer (CRC) remains one of the major triggers of deaths from malignant tumors. Globally more than 1 million people suffer from CRC every year [1]. As of 2012, CRC is the fourth cause of cancer deaths after lung, stomach, and liver cancer [2]. Despite improvements in diagnosis and combined treatment, patients with CRC have an even worse prognosis, especially in advanced patients. Therefore, it is quite urgent to clarify the mechanism, even potential approaches for the therapeutic intervention of CRC.
Accumulating evidence has shown that pseudogene, a type of long noncoding RNA (lncRNA), exhibits pivotal functions. It is estimated that the human genome has more than 18000 pseudogenes [3]. And pseudogene has emerged as key regulators of important biological processes involved in the development of human cancers. For instance, increased CYP4Z1 expression promotes tumor angiogenesis and growth in breast cancer partly via PI3K/Akt and ERK1/2 activation [4]. Nanog regulates primitive hematopoiesis by directly repressing critical erythroid lineage specifiers [5]. PTENP1 can exert a growth-suppressive role by regulating cellular levels of PTEN [6]. NAIL identified in colitis patients that regulates inflammation through NFkB which is important for the inflammation [7], and cross-talk of NFkB and HIF signaling pathways in mediating the occurrence of inflammation. LncRNA PTTG3P (pituitary tumor-transforming 3, pseudogene, NR 002734), located at chromosome 8q13.1, was first reported in the study of the human pituitary tumor transforming gene (hPTTG) family in 2000 [8]. PTTG3P plays a vital role in CRC, Liu et al. [9] revealed that lncRNA PTTG3P was up-regulated in CRC cancer tissues by lncRNA databases, and could function as an indicator for the prognosis in CRC. However, its biological function in M2 phenotype of macrophage is yet to be illustrated in CRC.
Our data discovered that PTTG3P predicts poor prognosis in patients with CRC. Further study revealed that PTTG3P facilitates cell growth by regulating hypoxia-inducible factor-1α (HIF1A)/PTTG3P/YAP1 axis, and induces M2 polarization of macrophages.

Clinical samples
A total of 120 patients with CRC were enrolled from the Affiliated Hospital of Youjiang Medical University for nationalities, the Central Hospital of Shenyang Medical Hospital, and Cancer Hospital of China Medical University between March 2010 and November 2015. The including criteria were as follows: patients with a definite pathological diagnosis; no patients were received with chemotherapy or radiotherapy before surgery, and the characteristics of cases were thoroughly noted. The tumor and paired non-tumor tissues were also collected after lesion excision with 30 min and stored in liquid nitrogen, then transferred to a −80 • C refrigerator. And the characteristics of cases were thoroughly noted. All the CRC patients have signed informed consent before utilizing the clinical resources for investigation aims. The study was approved by the Ethics Committee of Youjiang Medical University for nationalities and Cancer Hospital of China Medical University.

Cell lines culture
Five human CRC cell lines (HT-29, SW620, HCT-8, SW480, and HCT116) and normal human intestinal epithelial cell lines (FHC, NCM460) were obtained from ATCC (Manassas, VA, U.S.A.). And cultured according to their instructions. All cells were cultured in an incubator according to their instructions at 37 • C and in a humidified atmosphere with 5% CO 2 .

Total RNA isolation, qRT-PCR, and transfection
The expression levels of RNA were calculated by the qRT-PCR system. Total RNA was extracted by TRIzol Reagent (Invitrogen), and 1 μg of total RNA was reverse transcribed using the PrimeScript RT Reagent Kit (Perfect Real-Time; Takara). pcDNA3.1-PTTG3P, PTTG3P-containing lentiviral sequence vector (sh-PTTG3P) were purchased from GeneChem Corporation (Shanghai, China). CRC cells were transfected with plasmids in the presence of Lipofectamine 3000 (Invitrogen). After 48 h of transfection, cells were gathered for further use in the following experiments. The gene expression quantity was calculated using the 2 − C t method. The oligonucleotides were transfected into CRC cells using Lipofectamine 2000 (Invitrogen, U.S.A.) and the transfection efficiency was confirmed by qRT-PCR. The details are listed in Supplementary Tables S1 and S2.

Cell proliferation assay
Cell viability assay was carried out to analyze cell proliferation. Cell viability was estimated using CCK8 (CK04, DO-JINDO, Beijing, China), based on the manufacturer's instructions. Cells were seeded in 96-well culture plates. After incubation for the indicated time, a CCK-8 reagent (10 μl) was added to each well. Cell viability was measured with a microplate reader for the absorbance at a wavelength of 450 nm.

Flow cytometry of apoptosis
CRC cells in six-well plates were rinsed in phosphate buffer saline (PBS), and then were trypsinized and resuspended in 100 μl binding buffer added with 2.5 μl of fluorescein isothiocyanate (FITC) conjugated Annexin V and 1 μl of PI (Invitrogen). Fifteen minutes later, flow cytometry (BD Biosciences) was utilized for apoptotic cells.

Chromatin immunoprecipitation assay
CRC cells (2 × 10 6 ) were used for chromatin immunoprecipitation (ChIP) assay according to the manufacturer's protocol of hIP)assay kit (Millipore). The resulting precipitated DNA specimens were analyzed using PCR to amplify fractions of the PTTG3P promoter. The PCR products were resolved electrophoretically on a 2% agarose gel and visualized using Ethidium Bromide staining.

Animal study
HCT116 cells were transfected with sh-PTTG3P. A total of 1 × 10 7 indicated cells were subcutaneously injected into 4-week-old male nude mice. Tumor volume measured every 5 days. Isoflurane gas anesthesia was used in mice. 2-Deoxyglucose (2-DG) using in the study was injected into the abdominal cavity (1000 mg/kg, injected into the abdominal cavity). Oxaliplatin using in the study was injected into the abdominal cavity (5 mg/kg twice per week). After 35 days, the 20 mice were sacrificed by excessive carbon dioxide asphyxiation euthanasia, and the tumor weight was measured. The animal study was carried out in Youjiang Medical University following the Guide for the Care and Use of Laboratory Animals of the NIH. The present study had been approved by the Committee on the Ethics of Youjiang Medical University and China Medical University.

Statistical analysis
All the data were shown as the mean + − standard deviation, at least three independent experiments. The difference between two independent groups was analyzed by a two-tailed Student's t test, while multigroup comparison was made by ANOVA. Expression correlation between genes was analyzed by Pearson correlation analysis. Survival analysis was conducted using the Kaplan-Meier method and analyzed by the log-rank test. SPSS 22.0 (SPSS Inc., Chicago, IL, U.S.A.) was used to conduct statistical analyses, and differences were ensured when P-value was <0.05.

PTTG3P is highly expressed in CRC
To evaluate potential lncRNAs involved in mediating CRC progression, we examined the lncRNA expression profile (GSE 84983) (Supplementary Figure S1A). Comparison between CRC tumor tissues and adjacent normal tissues, we focused on the up-regulated lncRNAs (fold change > 5, P<0.01), for these lncRNAs might be oncogenes and therapeutic targets. LncRNA PTTG3P was one of the most up-regulated and chosen for consideration (Supplementary Figure S1B). Then, we found that PTTG3P had rarely the ability to code proteins, using the open-reading frames (ORFs) Finder and conserved domain database. Moreover, five other different online metrics got the same conclusion (Supplementary Table S3). Additionally, we identified no valid Kozak consensus sequence in PTTG3P [10], indicating that PTTG3P was an lncRNA with no protein-coding potential. Then we explored the subcellular location of PTTG3P by using lncRNA subcellular localization predictor software (lncLocator, http://www.csbio.sjtu.edu.cn/bioinf/lncLocator/) (Supplementary Figure S1C), suggesting PTTG3P was mainly localized to the cytoplasm, and subcellular fractionation confirmed the prediction (Supplementary Figure S1D).
Further, we explored PTTG3P expression in a cohort of 120 paired and non-tumor tissues of CRC, the clinicopathologic characteristics are demonstrated in Table 1. Significantly, the PTTG3P level was overexpressed in CRC tissues compared with their counterparts ( Figure 1D,E), which was in accordance with the results of The Cancer Genome Atlas (TCGA) database ( Figure 1F,G). Besides, high PTTG3P expression was observed in other malignant tumors ( Figure 1H). Also, our specimens confirmed PTTG3P overexpression in stomach adenocarcinoma (STAD) and esophageal squamous cell carcinoma (ESCA) ( Figure 1I,J). Intriguingly, there was no actionable EGFR, VEGFR or RAS mutations, indicating that higher expressed PTTG3P may be driven by oncogenic event (Supplementary Figure S1E-G). Altogether, these data revealed that PTTG3P was elevated in CRC and might be an oncogene.

High PTTG3P level correlates with poor prognosis
To identify the connection between the level of PTTG3P and clinicopathologic features, we divided the cases into PTTG3P low-and high-expression groups based on the median expression. Up-regulated PTTG3P was positively linked with tumor size (P=0.02) and differentiation (P=0.01), but not with age (P=0.86), gender (P=0.74), tumor invasion depth (P=0.28), lymph node metastasis (P=0.09) or vessel invasion (P=0.06) ( Table 1). Moreover, the PTTG3P expression in stage III-IV was higher than that in stage I-II tissues, and PTTG3P was expressed much more highly in advanced CRC samples than early CRC tissue ( Figure 2A). Additionally, Kaplan-Meier survival curves illustrated that patients with highly expressed PTTG3P had poorer survival time ( Figure 2B). Further, we determined the prognostic ability of PTTG3P in CRC. As shown in Table 2, univariate analyses suggested highly expressed PTTG3P was associated with a dramatic risk of death (P<0.01). Multivariate analysis demonstrated that PTTG3P expression was an independent prognostic factor (P<0.01). A model that incorporated the independent predictor was proposed as the nomogram (Supplementary Figure S2). Subsequently, the ROC curve was carried out to evaluate the diagnostic capacity of PTTG3P in CRC tissues compared with normal counterparts, the area under the ROC curve (AUC) was 0.776 (95% CI: 0.733-0.819) ( Figure 2C). Thus, these data suggested that high expression of PTTG3P predicted a worse prognosis and may serve as a clinical biomarker for CRC patients.

PTTG3P is caused by metabolic stress and promotes glycolysis and proliferation in CRC
To investigate the biological function of PTTG3P, we transfected the PTTG3P overexpressed plasmids and shRNA targeting PTTG3P into HT-29 and HCT116 cells, respectively ( Figure 2D). By determining PTTG3P expression via   gene set enrichment analysis (GSEA), TCGA profiles, we found that PTTG3P level was positively correlated with the glycolysis by affecting genes in glycolysis regulation ( Figure 2E). PTTG3P knockdown restrained the mRNA level of GLUT-1, ALDOA, PKM2 and LDHA, and the effect of sh-PTTG3P on glycolytic gene transcription could be rescued by PTTG3P re-expression ( Figure 2F). And glucose deprivation is a well-known feature of solid tumors. Subsequently, we wonder whether PTTG3P participated in cell survival under metabolic stress, then we carried out several experiments with different glucose concentration and glycolysis inhibitor 2-DG to make a condition of glucose deprivation. Obviously, PTTG3P expression was increased by glucose deprivation or 2-DG treatment in either dose-or time-dependent manner ( Figure 2G,H). Thus, we elucidated that PTTG3P could play a crucial role in the progression of metabolic stress.
Next, we performed the glucose uptake analysis, ATP analysis, lactate production analysis, and discovered that sh-PTTG3P repressed these phenomena. In contrast, PTTG3P overexpression boosted glucose uptake ( Figure 3A), lactate production ( Figure 3B), and ATP accumulation ( Figure 3C). Additionally, we calculated the level of ECAR, sh-PTTG3P notably repressed glycolytic capacity and vice versa ( Figure 3D). Also, we found that silenced PTTG3P suppressed the proliferation, facilitated apoptosis of HCT116 cells, whereas up-regulated PTTG3P increased the proliferation, inhibited apoptosis of HT-29 cells according to the CCK-8 assay and flow cytometry analysis ( Figure 3E,F). In vivo, highly expressed PTTG3P efficiently increased the tumor growth ( Figure 3G,H). We then explored whether glycolysis played a vital role in cell proliferation and tumor growth. Notably, the glycolic inhibitors 2-DG and 3-BP or depletion of LDHA, which catalyzed the final step of glycolysis, could partly abrogate cancer cell proliferation and tumor growth ( Figure 3I-K). Clinically, Oxaliplatin is used for the treatment of CRC. Previously, it is reported that suppression of glycolysis is an effective strategy to block cell proliferation and conquer drug resistance. As shown in  Figure 3L,M, PTTG3P depletion and Oxaliplatin played a synergistic role in emancipating tumor growth. As a taken, PTTG3P knockdown plus Oxaliplatin is a promising therapy for CRC.

PTTG3P regulates Hippo signaling pathway in CRC
In order to elucidate which pathway is involved in PTTG3P-mediated CRC progression, GSEA in the published TCGA CRC database was explored. And we suggested that PTTG3P expression was associated with the YAP1-activated gene signatures, indicating that Hippo signaling pathway might be involved ( Figure 4A). Then the hub genes in the Hippo pathway, including LATS1/2, MST1/2 and YAP1, and Hippo pathway target genes, such as CDX2, FOXM1, CTGF and CYR61, were checked in sh-PTTG3P HCT116 cells. Subsequently, PTTG3P knockdown impaired the level of YAP1, FOXM1 and CTGF ( Figure 4B). Moreover, silenced PTTG3P obviously decreased the enrichment of H3K27Ac at the YAP1 promoter, while that of H3K27me3 was increased ( Figure 4C). Next, we explored the TCGA database and drew the PTTG3P co-expression heat map. Then, we discovered that PTTG3P expression was correlated with genes in the Hippo pathway (YAP1, TEAD1-3), genes in the phenotype of proliferation (PCNA, MKI67, MCM2, MCM3, MCM5), gene in the phenotype of apoptosis (BAX, CASP1, CASP3, CASP10), gene in the phenotype of autophagy (ATG5), and genes in the phenotype of cell cycle (CDK1, CCDN1, CCND2, CCNB1), but not with MST1, an upstream factor of YAP1 ( Figure 4D). It is commonly acknowledged that YAP1, a crucial factor in the Hippo pathway, involves in cell proliferation and suppresses apoptotic genes, and YAP1 was highly expressed in CRC (Supplementary Figure S1H,I), and associated with advanced characteristics of CRC (Supplementary Table S4). Further, YAP1 had a higher diagnostic value (AUC = 0.793, 95% CI: 0.729-0.858) from the TCGA database (Supplementary Figure S1J). Additionally, we performed rescue assays in HT-29 cells. PTTG3P OE plus YAP1 KD could reverse the bioeffect of PTTG3P. Besides, we applied YAP1 inhibitors, CA3 (a novel specific YAP1 inhibitor) and verteporfin (an inhibitor of YAP1/TEAD interaction), got the same conclusion ( Figure 4E-H). Intriguingly, the treatment of Hippo pathway inhibitor, XMU-MP-1 (inhibiting MST1/2), could not recover the effect of PTTG3P on proliferation, apoptosis, and tumor growth ( Figure 4I,J). In brief, all the data uncovered that PTTG3P hedges the key factor MST1/2, while modulates YAP1 in the Hippo pathway to exhibit pivotal functions in CRC progression.

PTTG3P promotes M2 phenotype polarization of macrophage
Tumor microenvironment (TME) is composed of tumor cells and surrounding non-tumor stromal cells, mainly including tumor-associated macrophages (TAMs), endothelial cells, and carcinoma-associated fibroblasts (CAFs) [11,12]. With the occurrence of tumors, tumor cells secrete various chemokines to recruit monocytes to infiltrate tumor tissues and further promote their M2-type polarization. Importantly, M2-like TAMs can in turn accelerate tumor growth, promote tumor cell invasion and metastasis, and inhibit immune killing to promote tumor progression [13].

PTTG3P is regulated by HIF1A under hypoxic conditions
We next investigated the upstream factors which induced the elevated level of PTTG3P in CRC cells.
Then, we treated CRC cells with hypoxia or CoCl 2 (hypoxia chemical inducer) and found the PTTG3P expression in HT-29 cells was obviously increased as well as the elevation of HIF1A ( Figure 5F). While, depletion of HIF1A strikingly ameliorated PTTG3P expression in both normoxia and hypoxia conditions ( Figure 5G,H). We also indicated that HIF1A and PTTG3P had a positive correlation ( Figure 5E,I). HIF1A was associated with advanced characteristics of CRC from the TCGA-COAD database (Supplementary Table S5). Besides, HIF1A could partly rescue the effect of PTTG3P KD ( Figure 5J-L). In addition, ChIP analysis indicated that HIF1A enriched in the PTTG3P promoter region ( Figure 5M). Thus, the HIF1A/PTTG3P/YAP1 axis played a crucial role in CRC progression.

PTTG3P plays vital functions in CRC immunology
Cancer cells have high glucose uptake and glycolysis, resulting in a low level of glucose in the tumor, thus inhibiting the production of IFN-r by CD8 + T cells in the tumor. The tumor immune environment broadly participates in different malignant tumors, including CRC. Recently, the treatment of CRC with immune checkpoint inhibitors (ICIs) has provided a potential clinical treatment. Interestingly, our findings proposed that low PTTG3P expression relates with CD8 + T, NK and TFH cells infiltration in the microenvironment of CRC, not with Treg or macrophages infiltration, based on the TCGA database ( Figure 6A-F). And the results of ELISA showed that the level of inflammatory cytokines TNF-α, IL-1β and IL-6 were decreased with PTTG3P depletion ( Figure 6G). In addition, PTTG3P was correlated with TNF-α, IL-1β and IL-6 from the TCGA-COAD database ( Figure 5E).

Discussion
Pseudogenes are nonfunctional segments of DNA that resemble functional genes. Most arise as superfluous copies of functional genes, either directly by DNA duplication or indirectly by reverse transcription of an mRNA transcript. For instance, lncRNA HK2P1 increased the lactate production and glucose uptake in endometrial stromal cells [14]. lncRNA PTENP1 suppressed the PI3K/AKT signaling pathway and hindered the progression of HCC [15]. However, the role of PTTG3P in regulating CRC glycolysis has not been fully elucidated. Our study uncovered that PTTG3P boosted cell proliferation and glycolysis via the miR-1271-5p/PTTG3P/YAP1 axis. Our study verified that PTTG3P is highly expressed and has a potential diagnostic value, with an AUC of 0.776 (95% CI: 0.733-0.819) in CRC. Clinically, high PTTG3P expression considerably associates with tumor size and TNM stage as well as shorter survival time. These results confirmed that PTTG3P serves as a valuable prognostic biomarker and aids innovatively efficient therapies for CRC patients. Additionally, our findings stand in line with other research, Liu et al. [9] reported that PTTG3P was remarkably up-regulated in CRC tumor samples than that in normal samples. Zhou et al. [16] revealed that PTTG3P is a valuable resource for identification in HCC progression and is useful for biomarker development. Weng et al. [17] certified that PTTG3P facilitates cell proliferation, migration and invasion, and might serve as a new promising strategy for gastric cancer. Recently, PTTG3P expression has a relationship with breast cancer [18] and pancreatic cancer [19]. Thus, the oncogenic role of PTTG3P in malignant tumors is strongly suggested.
Malignant tumors could undergo glycolysis at a higher speed than non-tumor tissue controls [20][21][22]. This phenomenon is known as the Warburg effect [23]. The Warburg hypothesis demonstrates that malignant tumor is fundamentally caused by mitochondrial metabolism disorder. Doherty et al. [22] found that tumor lactate levels correlate with increased metastasis, tumor recurrence, and poor outcome. And targeting lactate metabolism is a prospective method for cancer therapeutics. Furthermore, cancer cells with a high level of glycolysis and acid resistance have an energetic growth advantage, which facilitates unrestrained proliferation and invasion. In our study, we explored gainand loss-of-function approaches in HT-29 and HCT116 cells and found PTTG3P ablation resulted in the inhibition of CRC cell glycolysis by regulating numerous genes linked with metabolic pathways, whereas the opposite outcome was observed after enforced expression of PTTG3P. Nowadays, a ketogenic diet was used to constrain glycolysis to starve cancer cells, adjusting mitochondrial metabolism [24]. Here, we also proposed that the biological mechanism of PTTG3P on boosting cell proliferation might resist apoptosis.
Hippo signaling pathway has become increasingly important in human cancer [25], the key regulator YAP1 has been certified to be up-regulated in breast cancer, CRC and liver cancer [26], and YAP1 could promote cell growth [27][28][29] and inhibit apoptosis [30]. Clinically, YAP1 could be a target for the development of cancer drugs [31]. Yi et al. [32] suggested that inhibiting TEAD-YAP1 interactions or block the binding function of WW domains is a pharmacologically viable strategy against the YAP1 oncoprotein. In our presented study, we discovered that PTTG3P activates Hippo signaling pathway by promoting YPA1, FOXM1 and CTGF, not MST1/2, and rescue assay consolidates this by using the Hippo pathway inhibitor, XMU-MP-1 (inhibiting MST1/2).
These years, cancer has been considered to be a complex system including the TME. TAMs are the most common immune-related stromal cells in the TME, and communication between cancer cells and TAMs is crucial for the progression of epithelial ovarian cancer (EOC) [33]. Recently, a great amount of studies elucidated that TAMs has the property of M2 macrophages that are related to cancer progression. Our finding suggested that PTTG3P expression was elevated in M2 macrophages, not in M1 macrophages or unpolarized macrophages. Next, enforced PTTG3P markedly increased the M2 macrophages markers (CD206 and MRC-2), while PTTG3P depletion attenuated these markers. Hence, PTTG3P might play a vital role in M2 polarization progression.
During cancer progression, tumor cells acquire comprehensive metabolic reprogramming, and tissue hypoxia is a prominent feature of solid tumors leading to cell metabolism adaptive changes. HIF1A is a key oxygen-regulated transcriptional activator, playing a fundamental role in the adaptation of tumor cells to hypoxia by up-regulating the transcription of target genes related to multiple biological processes, including cell survival, proliferation, angiogenesis and anti-apoptosis [34,35]. In the microenvironment of oxygen-glucose shortage, HIF1A is a major factor in cancer survival. Zhou et al. [36] proposed that HIF1A activated lncRNA RAET1K modulated glycolysis in hepatocellular carcinoma cells via miR-100-5p. Tiwari et al. [37] discovered that HIF1A might act as a tumor suppressor by preventing the expression of PPP1R1B and subsequent degradation of the p53 protein in pancreatic cancer cells. Our study indicated that HIF1A increased PTTG3P expression at the transcription level and enriched in the promoter region of PTTG3P, and could partly rescue the bioeffect of PTTG3P. Last, our findings proposed that low PTTG3P expression relates with CD8 + T, NK and TFH cells infiltration in the microenvironment of CRC, not with Treg or macrophages infiltration from the TCGA database. And we found that a decreased TNF-α, IL-1β and IL-6 level was accompanied by ablation of PTTG3P.
In summary, we found hypoxia-induced PTTG3P was an oncogene in CRC, and higher PTTG3P expression was related to a dismal prognosis. PTTG3P might facilitate cell proliferation, glycolysis and M2 phenotype of macrophage. Thus, PTTG3P plus HIF1A should be together adopted as critical targets for the prevention and therapy of CRC, illuminating some light on the understanding of the lncRNAs in CRC progression.

Data Availability
The datasets used and analyzed in the current study are available from the corresponding authors on reasonable request.