LncRNA UCA1, Upregulated in CRC Biopsies and Downregulated in Serum Exosomes, Controls mRNA Expression by RNA-RNA Interactions.

Long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) contribute to the onset of many neoplasias through RNA-RNA competitive interactions; in addition, they could be secreted by cancer cells into biological fluids, suggesting their potential diagnostic application. By analyzing the expression of 17 lncRNAs and 31 circRNAs in biopsies and serum exosomes from colorectal cancer (CRC) patients through qRT-PCR, we detected CCAT1, CCAT2, HOTAIR, and UCA1 upregulation and CDR1AS, MALAT1, and TUG1 downregulation in biopsies. In serum exosomes, UCA1 was downregulated, while circHIPK3 and TUG1 were upregulated. Combined receiver operating characteristic (ROC) curves of TUG1:UCA1 and circHIPK3:UCA1 showed high values of sensitivity and specificity. Through in vitro (i.e., RNA silencing and mitogen-activated protein kinase [MAPK] inhibition) and in silico analyses (i.e., expression correlation and RNA-RNA-binding prediction), we found that UCA1 could (1) be controlled by MAPKs through CEBPB; (2) sequester miR-135a, miR-143, miR-214, and miR-1271, protecting ANLN, BIRC5, IPO7, KIF2A, and KIF23 from microRNA (miRNA)-induced degradation; and (3) interact with mRNA 3'-UTRs, preventing miRNA binding. UCA1 and its co-regulated antisense LINC01764 could interact and reciprocally mask their own miRNA-binding sites. Functional enrichment analysis of the RNA-RNA network controlled by UCA1 suggested its potential involvement in cellular migration. The UCA1 regulatory axis would represent a promising target to develop innovative RNA-based therapeutics against CRC.

Introduction

Non-coding RNAs (ncRNAs) are a highly heterogeneous group of untranslated RNA molecules having regulatory functions within eukaryotic cells. Based on their size, ncRNAs can be classified as (1) long non-coding RNAs (lncRNAs), longer than 200 nt; and (2) small non-coding RNAs, including microRNAs (miRNAs), whose length may be 200 nt or less. LncRNAs regulate a plethora of fundamental biological processes (i.e., cell cycle, chromatin remodeling, histone modifications, splicing, and apoptosis) at different levels (transcriptional, post-transcriptional, translational, post-translational, and epigenetic). However, the exact molecular mechanisms of action have been elucidated just for a few of them.

A particular class of lncRNAs is represented by circular RNAs (circRNAs): these are covalently closed RNA molecules produced from protein-coding genes through a splicing-like process, called backsplicing, during which the free ends of the transcript are bound in the typical head-to-tail junction.3, 4 The molecular mechanisms of action of lncRNAs and circRNAs are under intense investigation, but evidence suggests their involvement in a complex competitive endogenous RNA (ceRNA) network based on RNA-RNA interactions, including those between ncRNAs and mRNAs.5, 6 Indeed, lncRNAs and circRNAs may act as miRNA sponges: that is, they can bind and sequester miRNAs through miRNA-binding sites included within their sequences.4, 7 Moreover, lncRNAs could directly bind the 3′-UTR of mRNAs, protecting them against miRNA-induced degradation.

In recent years, the discovery of ncRNAs has led to the investigation of their involvement in physiological and pathological processes. Moreover, ncRNAs have also been detected in biological fluids, often packaged within membrane-bound vesicles, such as exosomes.9, 10, 11, 12 Exosomes are small membrane vesicles characteristically cup shaped, with a diameter of 30–100 nm, containing proteins and RNAs (mRNAs, miRNAs, and lncRNAs);13, 14 recently, the presence of double-stranded genomic DNA within circulating exosomes has been reported. Exosomes are secreted by several, if not all, cell types, and they have been detected in all biological fluids analyzed to date. Their main function is to mediate cell-to-cell communication through their molecular cargo, both in physiological and pathological conditions. Several studies have reported a differential expression of ncRNAs within serum exosomes of patients affected by different diseases, including tumors; this has suggested the application of circulating ncRNAs for diagnostic and prognostic approaches.

Colorectal cancer (CRC) is one of the most common diseases in developed countries, and it represents one of the tumors showing the most relevant dysregulation of ncRNA functions.1, 19, 20, 21, 22 In the present study, we analyzed the expression of lncRNAs and circRNAs of oncological interest in tissues and serum exosomes from CRC patients, aiming to (1) increase our knowledge of the molecular bases of CRC, and (2) identify potential new biomarkers for its diagnosis and prognosis. This analysis led us to focus our attention on the lncRNA UCA1: through a combination of in vitro analyses (i.e., RNA silencing and inhibition of the mitogen-activated protein kinase [MAPK] pathway) and in silico approaches (i.e., expression correlation analysis and prediction of RNA-RNA binding), we reconstructed the ceRNA network controlled by UCA1 to get further insights into its molecular signaling.

Results

NcRNA Expression in CRC Cell Lines and Exosomes

To verify their potential synthesis in CRC biopsies, the expression of selected lncRNAs and circRNAs was investigated through Real-Time PCR in two CRC cell lines: HCT-116 and Caco-2. The results of this preliminary screening showed that almost all analyzed molecules are expressed in both cell lines (Figure S1A). We performed the same analysis on exosomes secreted by the same cells, with the aim of evaluating the presence of lncRNAs and circRNAs within these vesicles. Except for circZRAMB1, all tested molecules were detected within the exosomes secreted by both CRC cell lines, even if at very low expression levels. In Figure S1B, we show the expression of exosomal ncRNAs characterized by a Ct value <35 in at least one cell line.

NcRNA Expression in FFPE CRC Biopsies and Serum Exosomes

To pinpoint the ncRNA molecules to be assayed in the entire cohort of patients, we performed a trial experiment on 3 formalin-fixed, paraffin-embedded (FFPE) biopsy pairs and serum exosomes compared to 3 healthy individuals (data not shown). Based on their expression (Ct < 35), we selected for FFPE biopsy analysis 15 of 17 lncRNAs (CCAT1, CCAT2, CRNDE, H19, HOTAIR, MALAT1, MEG3, MIR17HG, PCAT1, PCAT6, PTENP1, TUG1, UCA1, WRAP53, and ZEB2AS1) and 6 of 31 circRNAs (CDR1AS, circCAMSAP1, circGON4L, circSPECC1, circZKSCAN1, and circZRAMB1). Real-Time PCR analysis on 20 FFPE CRC biopsies and their normal adjacent tissues (NATs) showed significant differential expression of 6 lncRNAs and 1 circRNA (Table 1). More specifically, 4 lncRNAs (CCAT1, CCAT2, HOTAIR, and UCA1) were upregulated, while 2 lncRNAs (MALAT1 and TUG1) and 1 circRNA (CDR1AS) were downregulated in CRC tumor tissues compared to NATs.

Table 1

lncRNA and circRNA Expression in CRC Biopsies Compared to Normal Adjacent Tissues and in Serum Exosomes of CRC Patients Compared to Healthy Individuals

	Biopsies	Exosomes
CCAT1	8.37 (0.0007)a	NA
CCAT2	4.67 (0.001)a	−1.13 (0.38)
CDR1AS	−3.3 (0.002)a	NA
circCAMSAP1	−1.09 (0.46)	−1.39 (0.42)
circGON4L	−1.13 (0.28)	NA
circHIPK3	NA	3.82 (0.036)a
circLRP6	NA	1.45 (0.21)
circSPECC1	−1.12 (0.31)	NA
circZKSCAN1	1.19 (0.23)	NA
circZRANB1	1.07 (0.29)	NA
CRNDE	1.23 (0.37)	−1.37 (0.45)
H19	−1.41 (0.75)	NA
HOTAIR	9.17 (0.0008)a	1.42 (0.16)
LIT1	NA	−1.15 (0.69)
MALAT1	−1.8 (0.004)a	1.41 (0.35)
MEG3	−1.73 (0.93)	NA
MR17HG	1.32 (0.63)	NA
PCAT1	−1.61 (0.76)	−1.18 (0.31)
PCAT6	1.44 (0.81)	−1.19 (0.57)
PTENP1	1.31 (0.83)	NA
TUG1	−2.04 (0.036)a	1.69 (0.029)a
UCA1	7.14 (0.001)a	−4.22 (0.03)a
WRAP53	−1.26 (0.83)	NA
ZEB2AS1	1.14 (0.27)	−1.36 (0.95)

The t test p-value (in parentheses) and median fold change are shown for each ncRNA. NA, not analyzed.

Statistically significant.

Similarly, we selected for exosome analysis 10 lncRNAs (CCAT2, CRNDE, HOTAIR, LIT1, MALAT1, PCAT1, PCAT6, TUG1, UCA1, and ZEB2AS1) and 3 circRNAs (circCAMSAP1, circHIPK3, and circLRP6). Real-Time PCR analysis showed that 2 lncRNAs and 1 circRNA were differentially expressed (DE) in serum exosomes of CRC patients compared to healthy individuals: in particular, UCA1 was downregulated, while circHIPK3 and TUG1 were upregulated (Table 1) (Figure 1A). We computed Pearson’s correlation coefficient (r) to evaluate the correlation of expression between DE ncRNAs in serum exosomes of all CRC patients and healthy individuals and in serum exosomes of CRC patients only. The Pearson coefficients showed a significant negative correlation between TUG1 and UCA1 in the entire cohort (Pearson = −0.45, p-value = 0.005), stronger in the CRC patient group (Pearson = −0.65, p-value = 0.004) (Figures 1B and 1C). Such a statistical correlation was not found in CRC tissues. No statistical correlation was observed between UCA1 and circHIPK3 and between TUG1 and circHIPK3.

Figure 1

Boxplots Showing Differential Expression of ncRNAs in CRC Biopsies and Serum Exosomes

(A) DE ncRNAs in serum exosomes of CRC patients compared to healthy individuals. (B) Pearson’s correlation coefficient (r) showed a negative correlation between TUG1 and UCA1 expression in serum exosomes of CRC patients and healthy individuals (C) and in serum exosomes of CRC patients only (D). *p ≤ 0.05 and **p ≤ 0.005.

Evaluation of Diagnostic Accuracy of DE ncRNAs in Serum Exosomes

We computed receiver operating characteristic (ROC) curves of DE ncRNAs using ΔCt values with respect to the endogenous control GAPDH. The circHIPK3 ROC curve showed an area under the curve (AUC) of 0.771 (95% confidence intervals [CIs], 0.508–0.936; p = 0.02) with 71% sensitivity and 80% specificity (ΔCt cutoff value ≤ 8.881197) (Figure 2A); the UCA1 curve showed an AUC of 0.719 (95% CIs, 0.533–0.863; p = 0.01) with 100% sensitivity and 43% specificity (ΔCt cutoff value ≤ 3.148002) (Figure 2B); the TUG1 curve was not significant (p = 0.2). The TUG1:UCA1 combined ROC curve, designed using ΔCt values calculated as CtTUG1 – CtUCA1, showed an AUC of 0.814 (95% CIs, 0.627–0.933; p = 0.0001) with 93% sensitivity and 64% specificity (ΔCt cutoff value ≤ 5.855948) (Figure 2C). The circHIPK3:UCA1 combined ROC curve, computed using ΔCt values calculated as CtcircHIPK3 – CtUCA1, showed an AUC of 0.9 (95% CIs, 0.648–0.992; p < 0.0001) with 100% sensitivity and 70% specificity (ΔCt cutoff value ≤ 12.651507) (Figure 2D).

Figure 2

ROC Curves of DE ncRNAs in Serum Exosomes of CRC Patients Compared to Healthy Individuals

ROC curves showing the sensitivity and specificity as CRC biomarkers of (A) circHIPK3 and (B) UCA1. (C) Combined ROC curve of TUG:UCA1 (CtTUG1 – CtUCA1). (D) Combined ROC curve of circHIPK3:UCA1 (CtcircHIPK3 – CtUCA1).

Analysis of GEO Datasets Confirms the Deregulation of DE ncRNAs

To validate the differential expression of ncRNAs observed in CRC FFPE biopsies and serum exosomes, we reanalyzed the expression of these ncRNAs in 15 CRC GEO datasets. The ncRNA expression data extracted from GEO also confirmed our data in larger cohorts of CRC tissues and normal colon mucosa. This analysis revealed a positive correlation between the up- or downregulation of each DE ncRNA and CRC progression and clinical and genetic features (such as tumor stage and size, relapse, metastasis, and CRC-associated mutations) (Figures 3A–3D). In addition, a dataset on peripheral blood of CRC patients compared to healthy individuals was analyzed, confirming the upregulation of TUG1 and the downregulation of UCA1 observed in our cohort (Figure 3E). Specifically, the same deregulation trend of both lncRNAs was observed in two different groups of CRC patients, divided according to Dukes classification: (1) Dukes stages A and B, and (2) Dukes stages C and D. In the same dataset, the negative correlation of expression between TUG1 and UCA1 (here described in serum exosomes) was confirmed (Pearson = −0.43, p-value = 0.01) (Figure 3F). Dataset screening provided no results for CCAT1, CCAT2, CDR1AS, and circHIPK3, because no specific probes for these ncRNAs were present in the analyzed datasets.

Figure 3

Validation of ncRNA Dysregulation in Tissues and Peripheral Blood of CRC Patients by GEO Dataset Screening

(A–D) GEO dataset analysis confirmed the upregulation of (A) HOTAIR and (B) UCA1 and the downregulation of (C) MALAT1 and (D) TUG1 in CRC tissues compared to normal colon mucosa. This analysis revealed a statistically significant positive association between ncRNA up- or downregulation and CRC progression and clinical and genetic features. Statistical significance is shown as –log10 (p value). Each dataset is represented by a different bar. MSI, microsatellite instability; CIN, chromosomal instability. (E) The dataset GEO: GSE10715 confirmed the upregulation of TUG1 and the downregulation of UCA1. The same trend is observed in two different groups of CRC patients, composed of (1) Dukes stages A and B, and (2) Dukes stages C and D. (F) Significant negative correlation of expression between TUG1 and UCA1.

Inhibition of MAPKs Affects ncRNA Expression

As the MAPK pathway plays a master control role in various CRC-related biological processes, we evaluated the effects of MAPK inhibition on the expression of ncRNAs showing altered expression in CRC in our study (Figures S2A and S2B). The expression of ncRNAs was investigated in HCT-116 cells at 6, 12, 24, and 48 hr after treatment with U0126. The results showed that the expression of DE ncRNAs was affected by the inhibition of the MAPK pathway (Table 2). More specifically, all ncRNAs, excluding HOTAIR and TUG1, significantly decreased their expression at 6 hr, while their expression levels increased at 24 hr, except for CCAT1, which was consistently downregulated at 6, 12, and 24 hr but upregulated at 48 hr.

Table 2

ncRNA Expression after U0126 Treatment

ncRNA	6 hr	12 hr	24 hr	48 hr
CCAT1	−3.6 (0.013)a	−2.34 (0.026)a	−2.48 (0.022)a	2.32 (0.015)a
CCAT2	−3.56 (0.05)a	−1.15 (0.39)	2.33 (0.001)a	2.5 (0.015)a
CDR1AS	−9.88 (0.002)a	−1.29 (0.24)	4.88 (0.004)a	−1.02 (0.36)
HOTAIR	−1.89 (0.18)	1.23 (0.53)	3.27 (0.01)a	2.8 (0.002)a
circHIPK3	−3.13 (0.023)a	1.18 (0.41)	2.26 (0.005)a	1.95 (0.033)a
MALAT1	−3.69 (0.017)a	1.43 (0.2)	2.61 (0.001)a	1.36 (0.13)
TUG1	−1.13 (0.068)	1.31 (0.45)	2.38 (0.005)a	2.77 (0.0004)a
UCA1	−4.29 (0.016)a	1.58 (0.12)	2.36 (0.006)a	2.14 (0.003)a

Median fold change and t test p-value (in parentheses) are shown.

Statistically significant.

To corroborate our findings on other experimental CRC models, we investigated within GEO datasets the relationship between the expression of post-treatment deregulated ncRNAs and the activation of the MAPK pathway. More precisely, we computed the expression correlation (Pearson coefficient) between each DE ncRNA and HSPA5, a transcriptional target of phospho-Erk (p-Erk). According to our data, we expected a positive correlation between HSPA5 and ncRNAs showing decreased expression at 6 hr post-treatment. Results showed both positive and negative correlations in different datasets, but UCA1 data were more consistent with a positive linear correlation with HSPA5, while MALAT1 showed high heterogeneity (Table S1). These data would hint to a hypothetical control of UCA1 expression induced by the activation of MAPK signaling.

Identification of Regulation Axes of UCA1

Because of its asymmetric deregulation between tissues and serum exosomes and its potential association with the MAPK pathway, we focused our subsequent in silico and in vitro analyses on UCA1 to get a better insight into its biomolecular roles and related regulatory circuits. We overlapped different tracks from the University of California, Santa Cruz (UCSC) Genome Browser, which had been derived from Encyclopedia of DNA Elements (ENCODE) experiments. This allowed us to identify the transcription factor-binding sites (TFBSs) for CEBPB, TEAD4, TFAP2A, and TFAP2C within 1 kb upstream of the UCA1 transcription start site. The expression of these TFs in CRC tumor tissues and their correlation of expression with UCA1 were investigated in GEO and The Cancer Genome Atlas (TCGA) datasets. CEBPB and TEAD4 showed the most significant upregulation and positive correlation with UCA1 (Table 3; Table S2). The expression of CEBPB and TEAD4 was also investigated after MAPK inhibition, showing a decrease for both TFs at 6 and 12 hr after treatment and an increase at 24 and 48 hr (Table 3). Importantly, the expression trends of CEBPB and UCA1 showed a significant linear correlation (Spearman = 0.66, p value = 0.0004), suggesting that CEBPB may be involved in UCA1 transcription regulation.

Table 3

Expression of TFs, miRNAs, mRNAs, and LINC01764 after U0126 Treatment

miRNAs Targeting mRNAs Co-expressed with UCA1	mRNAs/miRNAs	6 hr	12 hr	24 hr	48 hr
–	miR-135a	1.38 (0.24)	−1.03 (0.82)	1.12 (0.73)	−1.73 (0.14)
–	miR-143	1.29 (0.15)	1.5 (0.028)a	1.56 (0.027)a	1.16 (0.43)
–	miR-214	−1.94 (0.88)	−2.34 (0.23)	1.13 (0.61)	3.14 (0.02)a
–	miR-1271	1.23 (0.21)	1.03 (0.71)	1.31 (0.91)	−1.36 (0.14)
miR-135a, miR-214, miR-1271	ANLN	−1.38 (0.41)	−1.02 (0.6)	−1.21 (0.26)	−4.19 (0.039)a
miR-135a, miR-143	BIRC5	−4.19 (0.032)a	1.23 (0.61)	1.74 (0.23)	4.57 (0.008)a
miR-135a, miR-143, miR-214	CD46	1.55 (0.23)	−1.76 (0.05)a	2.13 (0.014)a	1.27 (0.07)
miR-135a, miR-214	DEK	−4.3 (0.014)a	−1.15 (0.5)	−1.31 (0.18)	−1.26 (0.2)
miR-143	DNMT3A	−1.87 (0.11)	1.41 (0.75)	1.92 (0.009)a	1.15 (0.72)
miR-143, miR-1271	HMMR	−2.13 (0.17)	−1.01 (0.59)	1.34 (0.005)a	−2.62 (0.07)
miR-214	IPO7	−1.13 (0.41)	−1.35 (0.25)	1.77 (0.09)	1.08 (0.84)
miR-135a, miR-143, miR-1271	KIF23	−1.87 (0.19)	−1.25 (0.28)	1.2 (0.17)	−2.53 (0.05)a
–	LINC01764	−2.65 (0.08)	−1.22 (0.23)	2.83 (0.004)a	4.8 (0.008)a
miR-143	MACC1	−3.11 (0.05)a	1.4 (0.034)a	2.55 (0.0006)a	1.54 (0.06)
miR-143	MMP7	5.01 (0.25)	2.37 (0.06)	2.2 (0.09)	1.19 (0.83)
miR-135a	MYC	−5.29 (0.013)a	−4.98 (0.023)a	1.25 (0.08)	2.21 (0.09)
–	CEBPB	−3.39 (0.01)a	−1.08 (0.55)	2.08 (0.001)a	1.72 (0.046)a
–	TEAD4	−1.99 (0.049)a	−1.28 (0.21)	3.43 (0.002)a	1.64 (0.07)

Median fold change and t test p-value (in parentheses) are shown.

Statistically significant.

To investigate the RNA-based signaling of UCA1, we hypothesized that UCA1 could regulate mRNA decay by sequestering miRNAs, thus preventing them from repressing their mRNA targets. Accordingly, we expected that in the module UCA1:miRNAs:mRNAs, UCA1 and mRNAs would show a significant co-expression. Based on sequence complementarity between miRNA seed regions and the lncRNA sequence, the miRcode tool allowed us to identify miRNAs that could be sponged by UCA1. Among them, we selected as hypothetical targets of UCA1 4 miRNAs (miR-135a, miR-143, miR-214, and miR-1271), showing from 1 to 4 binding sites within the UCA1 sequence and reduced expression in CRC according to GEO dataset screening (Table S3). Experimentally validated targets of the 4 miRNAs were retrieved from StarBase, miRTarbase, and the literature (Table S4), while their expression in CRC was assessed by Oncomine and analyzed in GEO and TCGA expression datasets. In Table S5, we list the 14 miRNA targets selected, which are characterized by upregulation in CRC and have positive correlation with UCA1.

The expression of miR-135a, miR-143, miR-214, and miR-1271 was investigated after inhibition of the MAPK pathway: a statistically significant upregulation of miR-143 was observed at 12 and 24 hr after treatment, while miR-214 showed the same deregulation trend at 48 hr after treatment (Table 3). Similarly, we assayed the expression of mRNA targets of miR-143 and miR-214 after the same treatment, showing significant alterations of ANLN, BIRC5, CD46, DEK, DNMT3A, KIF23, HMMR, and MACC1 at different time points (Table 3). A significant positive correlation was shown between expression levels of UCA1 and BIRC5 (Spearman = 0.39, p-value = 0.05), CD46 (Spearman = 0.79, p-value < 0.0001), DNMT3A (Spearman = 0.75, p-value < 0.0001), HMMR (Spearman = 0.4, p-value = 0.05), and MACC1 (Spearman = 0.43, p-value = 0.03). Although IPO7 expression was not significantly altered by U0126 treatment, it was strongly positively correlated to UCA1 expression (Spearman = 0.76, p-value < 0.0001). These data would suggest that expression correlation between UCA1 and some of its potential targets persisted also after the in vitro perturbation induced by U0126 treatment. Despite the fact that miR-135a expression was not altered by U0126 treatment, we investigated the expression of MYC, which was selected as a target of miR-135a. It has been reported that MYC expression is promoted by CCAT1 and CCAT2,24, 25 whose expression was altered by MAPK inhibition. MYC showed a significant downregulation at 6 and 12 hr post-treatment, similar to the expression trend of CCAT1 and CCAT2 (Table 3). Moreover, MYC and CCAT2 showed a significant positive correlation of expression (Spearman = 0.89, p-value < 0.0001).

We found that both long and short variants of LINC01764, a non-coding antisense on the UCA1 locus, include a region longer than 100 nt showing sequence complementarity on the 3′-UTR of UCA1. More specifically, the LINC01764 short variant would bind the UCA1 region harboring the miR-214-binding site, while the complementary region between UCA1 and the LINC01764 long variant overlapped the binding sites of miR-135a and miR-143 on the LINC01764 sequence (Figure 4A). The expression of LINC01764 after MAPK inhibition was decreased at 6 and 12 hr after treatment and significantly increased at 24 and 48 hr (Table 3). Moreover, LINC01764 levels were parallel to the expression of UCA1 during the same time course (Spearman = 0.42, p-value = 0.03). These data would suggest a potential physical interaction between UCA1 and LINC01764, which could reciprocally mask their own miRNA-binding sites and, accordingly, explain the expression correlation observed between UCA1 and its antisense.

Figure 4

Prediction of UCA1 RNA Interactions

(A) Complementarity between UCA1 and the two variants of its antisense RNA (LINC01764). The white areas show the complementary regions between UCA1 3′-UTR and the long and short variants of LINC01764; the black boxes represent the miRNA-binding sites overlapping complementary regions. (B) Graphical representation of free energy (kJ) computed from the binding of UCA1 and the 3′-UTRs of its co-expressed mRNAs. Gray-coded scale, based on free energy values, is used to depict the 3′-UTRs of mRNAs. White regions depict the binding regions with minimum free energy, and, accordingly, they represent the most reliable binding regions. The boxes represent miRNAs binding to minimum free energy regions: white boxes, tumor suppressor miRNAs in CRC; black boxes, onco-miRNAs in CRC. Free energy scale limits: (1) ANLN: minimum = −22.8, maximum = −0.28; (2) BIRC5: minimum = −27.32, maximum = −2.2; (3) IPO7: minimum = −30.17, maximum = 0.82; (4) KIF2A: minimum = −15.11, maximum = 3.01; and (5) KIF23: minimum = −15.01, maximum = −0.38.

UCA1 Target Expression after In Vitro Silencing

The miRNAs potentially sponged by UCA1 were moderately upregulated at 24 hr post-transfection (PT) with ASOs (antisense oligonucleotides). In particular, miR-214 and miR-1271 showed the highest fold change values at 24 hr PT (Table 4). Moreover, UCA1 expression levels showed a significant negative correlation with miR-135a (Spearman = −0.8, p-value = 0.01) and miR-143 (Spearman = −0.52, p-value = 0.02). In the same way, expression analysis was also performed on experimentally validated miRNA targets: ANLN, BIRC5, IPO7, KIF2A, and KIF23 showed a statistically significant downregulation at 48 hr PT (Table 4). These results would show a potential positive control of UCA1 on ANLN, BIRC5, IPO7, KIF2A, and KIF23 expression. Although LINC01764 was downregulated at 24 and 48 hr after UCA1 silencing, its variation was not statistically significant, but it did show a significant positive correlation with UCA1 expression (Spearman = 0.47, p-value = 0.04).

Table 4

Expression of miRNAs, mRNAs, and LINC01764 after UCA1 Silencing

miRNAs Targeting mRNAs Co-expressed with UCA1	miRNAs/mRNAs	24 hr	48 hr
–	miR-135a	1.26 (0.25)	1.24 (0.26)
–	miR-143	1.11 (0.56)	−1.15 (0.39)
–	miR-214	1.61 (0.3)	−1.12 (0.81)
–	miR-1271	1.76 (0.34)	−3 (0.27)
miR-135a, miR-214, miR-1271	ANLN	−1.23 (0.34)	−1.56 (0.003)a
miR-135a, miR-143	BIRC5	−1.21 (0.05)a	−1.43 (0.009)a
miR-135a	BZW2	1.01 (0.59)	−1.29 (0.072)
miR-135a, miR-143, miR-214	CD46	−1.19 (0.22)	−1.02 (0.75)
miR-135a, miR-214	DEK	−1.29 (0.3)	−1.09 (0.58)
miR-143	DNMT3A	−1.07 (0.19)	−1.34 (0.13)
miR-143, miR-1271	HMMR	−1.32 (0.21)	−1.12 (0.37)
miR-214	IPO7	−1.35 (0.19)	−1.72 (0.004)a
miR-1271	KIF2A	−1.32 (0.38)	−1.9 (0.008)a
miR-135a, miR-143, miR-1271	KIF23	−1.26 (0.076)	−1.74 (0.007)a
–	LINC01764	−1.49 (0.16)	−1.18 (0.58)
miR-143	MACC1	−1.21 (0.23)	−1.32 (0.18)
miR-143	MMP7	−1.22 (0.56)	1.34 (0.072)
miR-135a	MYC	1.08 (0.42)	1.24 (0.52)
miR-135a	TRIP13	−1.05 (0.27)	−1.2 (0.61)

Median fold change and t test p-value (in parentheses) are shown.

Statistically significant.

Prediction of Direct Binding between UCA1 and the 3′-UTR of ANLN, BIRC5, IPO7, KIF2A, and KIF23

Additionally, we tested the hypothesis that mRNA downregulation induced by UCA1 silencing could be the effect of its direct binding on mRNA 3′-UTRs, which would protect them from miRNA binding and, thus, increase mRNA stability. Using the IntaRNA tool, we computed the potential RNA-RNA interactions among UCA1 and the 3′-UTRs of mRNAs downregulated PT. We found that binding of UCA1 on the 3′-UTRs of mRNAs showed several regions with strong negative free energy (NFE) (i.e., index of very effective binding), which overlapped several miRNA-binding sites identified by the microrna.org database (Figure 4B; Table S6). More specifically, we detected 17, 5, 16, 10, and 10 miRNAs whose binding sites lay in such NFE regions for ANLN, BIRC5, IPO7, KIF2A, and KIF23, respectively. About 80% of these 44 miRNAs exert a tumor suppressor role in CRC, according to GEO datasets (Table S7). These findings could suggest that UCA1, by binding the 3′-UTRs of mRNAs, would mask binding sites of miRNAs with a potential tumor-suppressive role in CRC, thus decreasing miRNA-induced degradation of mRNAs.

Functional Enrichment Analysis of the ceRNA Network Controlled by UCA1

To evaluate the biological functions potentially regulated by the ceRNA network controlled by UCA1, we computationally created a network of molecular interactions among Anln, Birc5, Ipo7, Kif2a, and Kif23 proteins and their first interactors. We performed functional enrichment analysis on all network nodes, showing a statistical over-representation of pathways associated with CRC carcinogenesis (e.g., MAPK, TGF-β, and Wnt pathways) or regulating cancer-related processes, such as cell cycle, apoptosis, adhesion, and cytoskeleton organization and migration; interestingly, processes related to mRNA stability and processing were also over-represented (Figure S3; Table S8).

Discussion

NcRNA Dysregulation Is Associated with the CRC Phenotype

Our data show that 6 lncRNAs (CCAT1, CCAT2, HOTAIR, MALAT1, TUG1, and UCA1) and 1 circRNA (CDR1AS) are differentially expressed in tumor tissues compared to normal mucosa; in particular, we observed the upregulation of CCAT1, CCAT2, HOTAIR, and UCA1 and the downregulation of CDR1AS, MALAT1, and TUG1. We confirmed these data by statistical re-analysis of expression microarray datasets from large cohorts of CRC patients with probes specific for ncRNAs. Interestingly, upregulation of both HOTAIR and UCA1 was statistically associated with several worse prognostic factors (e.g., high-grade tumors, tumor size, microsatellite instability, and BRAF and P53 mutations). Similarly, downregulation of MALAT1 and TUG1 was related to tumor size, microsatellite instability, and young age at diagnosis and recurrence interval. These data strongly suggest that the alteration of these lncRNAs could contribute to the progression of CRC.

High HOTAIR and UCA1 levels had been previously observed in several tumors, including CRC. HOTAIR upregulation has also shown a positive correlation with poor prognosis, epithelial-mesenchymal transition (EMT), and metastasis in CRC patients.26, 27, 28 UCA1 expression is increased in tumor tissues of patients with lymph node metastasis compared to non-metastatic ones, and its upregulation correlates with tumor stage, poor prognosis, and metastasis. Although the downregulation of MALAT1 and TUG1 observed here was in agreement with their decreased expression in several other tumors,31, 32, 33, 34, 35 it differed from other reports about CRC, in which the upregulation of MALAT1 and TUG1 was related to poor prognosis and metastasis.36, 37, 38, 39 It is noteworthy that the published data report a very complex scenario about MALAT1 and TUG1 expression in different tumors, and, in some cases, opposite functions have been reported in the same cancer model.33, 35, 40, 41 It is possible that the molecular functions involving lncRNAs are strongly dependent on the specific molecular context in which lncRNAs act, which is highly heterogeneous also within the same cancer type. The expression dysregulation reported here for CCAT1, CCAT2, and CDR1AS has not been supported by statistical re-analysis of CRC microarray datasets because no specific probe for these ncRNAs was available. Previous studies reported the upregulation of CCAT1 and CCAT2 in CRC tissues and their involvement in proliferation, invasion, and migration;25, 42, 43, 44 a recent paper reported the upregulation of CDR1AS in CRC tissues compared to normal colon mucosa.

Interestingly, the expression of all DE ncRNAs analyzed was modified by drug-induced inhibition of the MAPK pathway, whose hyperactivation is considered a key event for CRC onset and progression.46, 47 It is well known that the molecular effects of U0126 appear a few hours after treatment; accordingly, we observed a global downregulation all ncRNAs at 6 hr after treatment, except for HOTAIR and TUG1. On the other hand, ncRNAs showed a global upregulation at 24 and 48 hr after treatment (when the drug effect had likely terminated), which could represent a compensation effect due to pathway reactivation. The early down-modulation of ncRNAs after MAPK inhibition partially matched the positive expression correlation between HSPA5 (transcriptional target of p-Erk) and MALAT1 and UCA1. These data would suggest a potential transcriptional control of these ncRNAs from MAPK signaling, which seems more consistent for UCA1. Indeed, very reliable TFBSs for CEBPB and TEAD4 were found on the upstream regulatory region of UCA1. Dysregulation of CEBPB and TEAD4 during U0126 treatment was analogous to the UCA1 deregulation trend; moreover, expression of CEBPB and UCA1 showed a strong significant linear correlation. This would suggest that CEBPB may act as a transcription activator for UCA1 through the MAPK pathway. This hypothesis is in agreement with several papers reporting that CEBPB harbors a highly conserved MAPK consensus site, whose phosphorylation operated by MAPKs regulates its activity.48, 49, 50 Accordingly, CEBPB has been associated with several tumor models, including CRC, where it has shown increased expression compared to normal mucosa.51, 52

NcRNAs in Serum Exosomes of CRC Patients May Be Applied in Non-invasive Diagnosis

Expression analysis of ncRNAs in serum exosomes of CRC patients showed the upregulation of circHIPK3 and TUG1 and the downregulation of UCA1. circHIPK3 expression has been investigated in several cancer cell lines, including CRC cells, showing higher levels compared to the linear HIPK3 transcripts. Its presence in serum exosomes was previously reported, but no association with CRC was observed. We describe for the first time the presence of TUG1 within exosomes from human serum and its deregulation in CRC patient serum. The possible application of circulating TUG1 as a diagnostic and prognostic biomarker has already been reported in other cancer models.55, 56 To date, UCA1 deregulation in serum exosomes of CRC patients has never been reported; on the contrary, its upregulation was observed in CRC patient plasma. However, the analysis of the dataset GEO: GSE10715 allowed us to validate our results in an independent cohort of 19 CRC patients and 11 healthy individuals, confirming the upregulation of TUG1 and the downregulation of UCA1, as well as their negative expression correlation. Our data showed an opposite deregulation trend in CRC tissues and in serum exosomes for TUG1, which is downregulated within tumor cells and upregulated in serum exosomes. This observation may suggest that tumor cells protect themselves from TUG1 tumor-suppressive activity by secreting the lncRNA via exosomes. Similar to TUG1, UCA1 showed an asymmetric distribution between tumor tissues, where the lncRNA is upregulated, and serum exosomes, where it is downregulated. This opposite deregulation would suggest that the oncogenic function of UCA1 is crucial for tumor progression, inducing tumor cells to retain the lncRNA by limiting its secretion through exosomes.

ROC curves showed good accuracy for 2 of 3 DE ncRNAs (circHIPK3: 71% sensitivity, 80% specificity; UCA1: 100% sensitivity, 43% specificity), with significant p-values. On the contrary, the TUG1 ROC curve showed no statistical significance. The goal of using diagnostic and prognostic biomarkers is to detect a disease with the highest possible accuracy, and a single biomarker is often not sufficient. For this reason, the use of a signature of multiple biomarkers is preferable, as it would increase diagnostic and prognostic accuracy. We thus computed combined ROC curves considering two ncRNAs with opposite trends of expression at a time (ΔCt values obtained as CtTUG1 – CtUCA1 and as CtcircHIPK3 – CtUCA1 were used to compute two different curves).60, 61 Both combined ROC curves were statistically significant, showing lower p-values, higher AUC values, and a better combination of sensitivity and specificity values than single ROC curves, thus demonstrating an increased diagnostic accuracy. These data support the hypothesis that the combination of two ncRNAs may diagnose CRC with higher sensitivity and specificity compared to a single biomarker. Further analyses on larger cohorts will be necessary to evaluate the diagnostic accuracy of these molecules as CRC biomarkers.

UCA1 Potentially Contributes to CRC Pathogenesis by Regulating an RNA Network

Several studies reported that UCA1 may act as an miRNA sponge in different in vitro models. It has been reported that enforced expression of UCA1 resulted in decreased expression of its miRNA targets, or, vice versa, reduced expression of UCA1 caused increased levels of miRNAs.62, 63, 64, 65 Based on this evidence, we investigated the effect of UCA1 silencing on its hypothetical miRNA targets, miR-135a, miR-214, miR-143, and miR-1271, observing a slight trend of upregulation for all miRNAs at 24 hr after transfection. We also analyzed the expression of experimentally validated mRNA targets of these miRNAs, selected because of their positive correlation of expression with UCA1.

Among miRNA targets, ANLN, IPO7, KIF2A, and KIF23 showed a significant reduced expression at 48 hr after transfection, while BIRC5 levels were decreased both at 24 and 48 hr. It is conceivable that miRNAs sponged by lncRNAs are no longer available for binding the 3′-UTR of mRNA targets, without alterations of their expression levels. This may explain why no statistically significant alterations were observed in miRNA expression after UCA1 silencing. According to the hypothesis of miRNA sponging, when UCA1 levels decrease, miRNAs are no longer sequestered and are free to bind their mRNA targets, which, in turn, show decreased expression. Indeed, the downregulation of ANLN, IPO7, KIF2A, and KIF23 was observed at 48 hr after UCA1 silencing, while BIRC5 levels were already decreased at 24 hr. Intriguingly, such downregulation could be due to a direct protection on mRNAs performed by UCA1, which could bind their 3′-UTR and prevent miRNA binding. Following this second hypothesis, we computationally found reliable interactions between UCA1 and the 3′-UTRs of ANLN, BIRC5, IPO7, KIF2A, and KIF23 that overlapped several miRNA-binding sites of tumor-suppressor miRNAs.

Finally, we hypothesized that UCA1 and its antisense LINC01764 could be involved in a mechanism of reciprocal regulation of expression, due to a mutual protection from binding of miR-214 on UCA1 and miR-143 and miR-135a on the LINC01764 long isoform. Accordingly, LINC01764 and UCA1 showed a significant positive correlation of expression after U0126 treatment and UCA1 silencing. Taken together, this evidence would suggest the existence of a co-regulation mechanism of UCA1 and LINC01764 based on RNA interactions. The biological role of LINC01764 still remains unknown.

The UCA1 network showed strong over-representation of processes related to actin and cytoskeleton organization, which underlie the migration phenomenon. Accordingly, UCA1 was previously reported to induce migration in cancer.66, 67, 68 Moreover, functional enrichment analysis showed the over-representation of processes regulating mRNA stability and processing. These results support our hypothesis, suggesting that UCA1 controls key biological processes promoting CRC progression by regulating the expression of ANLN, BIRC5, IPO7, KIF2A, and KIF23.

Conclusions

Based on our hypotheses, UCA1 may be an RNA regulator of CRC progression by controlling a ceRNA network, thus fostering the upregulation of ANLN, BIRC5, IPO7, KIF2A, and KIF23 in two ways: (1) by sponging miRNAs negatively regulating their expression, and (2) by directly binding the 3′-UTRs of mRNAs and protecting them from miRNA-mediated degradation. The existence of this complex RNA-based regulatory signaling, which controls cancer-related pathways, paves the way to innovative therapeutic applications of RNA-based anticancer strategies targeting UCA1.

Materials and Methods

Cell Lines

HCT-116 and Caco-2 cell lines, both derived from primary tumors, were obtained from the Interlab Cell Line Collection (ICLC), an International Certified Repository Authority within the IRCCS Azienda Ospedaliera Universitaria San Martino-IST Istituto Nazionale per la Ricerca sul Cancro (Genova, Italy). HCT-116 cells were cultured in RPMI-1640 medium (Gibco, Thermo Fisher Scientific, Waltham, MA), supplemented with 10% fetal bovine serum (FBS) (Gibco) and 2 mM L-glutamine (Lonza, Basel, Switzerland); Caco-2 cells were cultured in DMEM (Sigma-Aldrich, St. Louis, MO), supplemented with 20% FBS (Gibco), 2 mM L-glutamine (Lonza), and 1% non-essential amino acids (NEAA) (Lonza); both media were supplemented with 1% penicillin/streptomycin (10,000 U/mL) (Gibco). Cells were cultured at 37°C and 5% CO2.

Isolation of Exosomes Secreted by CRC Cell Lines

Exosomes secreted by both cell lines were isolated from culture media 3 days after seeding in 175-cm2 flasks (6 × 106 cells/flask), as previously reported. Briefly, medium was harvested and centrifuged at 300 × g (4°C, 10 min) on a Beckman J-6M/E centrifuge (rotor TY.JS 5.2) (Beckman Coulter, Brea, CA) to pellet debris. The supernatant was then centrifuged at 16,500 × g (4°C, 30 min) on a Beckman L8 70M ultracentrifuge (rotor SW28), and then filtered through a 0.2-μm filter. The final supernatant was centrifuged at 120,000 × g (4°C, 70 min) with the same ultracentrifuge. Exosome pellet was directly lysed for RNA isolation. FBS added to culture medium had been previously deprived of exosomes through the same protocol, removing the exosome pellet and harvesting the FBS supernatant for subsequent use.

RNA Isolation from CRC FFPE Samples and Serum Exosomes

20 CRC patients were recruited at Azienda Ospedaliero-Universitaria Policlinico Vittorio Emanuele (Catania, Italy). FFPE tumor tissues and distal normal colon mucosa were isolated at the Section of Anatomic Pathology, Department G.F. Ingrassia, University of Catania (Catania, Italy). RNA was extracted from FFPE samples through a PureLink FFPE RNA Isolation Kit (Invitrogen), according to the manufacturer’s instructions. Serum of these 20 CRC patients and 20 unaffected individuals was collected through venous blood sampling after having signed informed consent for the use of their blood, in accordance with the Declaration of Helsinki. This study was approved by the Ethical Committee of “Vittorio Emanuele” Hospital (Azienda Ospedaliero-Universitaria Policlinico Vittorio Emanuele, Catania). To separate serum from cellular components, blood samples were centrifuged at 15,000 × g (4°C, 15 min), and then the supernatant was harvested and centrifuged again under the same conditions to completely remove cells or debris, as these could alter subsequent analysis. Supernatant serum was harvested and mixed to ExoQuick (System Biosciences, Palo Alto, CA) to precipitate nanovesicles, according to the manufacturer’s instructions. Finally, the exosome pellet was directly lysed with TriZol for total RNA isolation. RNA was quantified by GenQuant pro spectrophotometer (Biochrom, Cambridge, UK) and Qubit (Invitrogen, Thermo Fisher Scientific, Waltham, MA), and treated with DNase (Invitrogen), according to the manufacturer’s instructions. Notably, both procedures used for exosome isolation (ultracentrifugation for CRC cell exosomes and Exoquick for serum exosomes) do not allow the discrimination of exosomes from high-density and low-density lipoproteins (HDL and LDL, respectively), because all of these particles have a similar density.

RNA Isolation from Cell Lines and Exosomes

Total RNA was extracted from cell lines and exosomes with TriZol (Invitrogen), according to the manufacturer’s instructions, and quantified by GenQuant pro spectrophotometer (Biochrom) and Qubit fluorescence quantification system (Invitrogen).

PCR Primer Design

By literature mining, we selected for this study 17 lncRNAs involved in CRC or other tumors and regulating cancer-related processes, such as cell cycle, apoptosis, gene expression, and splicing. We designed specific PCR primers for the selected lncRNAs and two housekeeping genes used for normalization: PPIA (peptidylprolyl isomerase A) and GAPDH (glyceraldehyde 3-phosphate dehydrogenase). Primer design was performed by using the online tool Primer-Blast (https://www.ncbi.nlm.nih.gov/tools/primer-blast/). Primer pairs are shown in Table S9. PCR primers for 31 circRNAs were retrieved from the paper published by Memczak and colleagues.

NcRNA Expression Analysis by Real-Time PCR

We investigated the expression of selected lncRNAs and circRNAs in CRC patient tissues and serum exosomes and in CRC cell lines and exosomes through Real-Time PCR by using Power SYBR Green RNA-to-CT 1-Step kit (Applied Biosystems, Thermo Fisher Scientific, Waltham, MA), according to the manufacturer’s instructions. All Real-Time PCR reactions were performed on a 7900HT Fast Real-Time PCR System (Applied Biosystems). DE lncRNAs and circRNAs were identified through SDS RQ Manager 1.2 software (Applied Biosystems); normalization was performed using PPIA and GAPDH for tumor tissues and serum exosomes, respectively. miRNA expression was assessed through TaqMan microRNA assays and TaqMan Universal Master Mix II, no UNG (Applied Biosystems); miR-24 was used as endogenous control. Expression fold changes of DE ncRNAs were calculated by applying the 2−ΔΔCt method. Statistical analysis was performed by paired two-tailed t test to compare tissue sample ΔCts, while unpaired two-tailed t test was used for ΔCts from the following: (1) serum exosomes, (2) U0126 treatment, and (3) UCA1 silencing. Statistical significance was established at a p-value ≤ 0.05.

ROC Curve Analysis

To evaluate the diagnostic accuracy of ncRNAs deregulated in patient serum exosomes, ROC curves were computed through the software MedCalc v15.11.4. The AUC and 95% CIs were calculated to assess the accuracy of sensitivity and specificity parameters and find the appropriate cutoff point. Statistical significance was established at a p-value ≤ 0.05. Combined ROC curves were also computed by using ΔCts calculated as CtncRNA1 – CtncRNA2, combining two ncRNAs showing an opposite trend of expression.

CRC Cell Line Treatment with MAPK Inhibitor U0126

MAPK inhibition was achieved by treating the HCT-116 cell line with U0126, a highly selective ATP-non-competitive inhibitor of both Mek1 and Mek2 kinases. Cells (3.2 × 104 per well) were seeded in 24-well plates and cultured in serum starvation conditions (no FBS) for 24 hr; successively, cells were treated with 25 μM U0126 (MEK1/2 inhibitor, Merck, Darmstadt, Germany). Cells were exposed to drug for 6, 12, 24, and 48 hr. Control samples were treated with an equivalent volume of DMSO (solvent of the drug used for treatment). All experiments were performed in biological triplicates. Cell viability was assessed through the MTT [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide] assay at 24 hr post-treatment.

Protein Extraction and Western Blot Analysis

MAPK pathway inhibition was verified by analyzing p-ERK levels through western analysis. We performed protein extraction through radioimmunoprecipitation assay (RIPA) buffer and quantification through Qubit. Blotting was performed by iBlot Dry Blotting System (Invitrogen). Proteins were detected using the ECL Plus Western Blotting detection (Amersham Biosciences, Little Chalfont, Buckinghamshire, UK). Membranes were incubated overnight with monoclonal antibodies to p44/42 MAPK (Erk1/2) and phospho-p44/p42 MAPK (Erk1/2) (Cell Signaling, Leiden, the Netherlands); beta-actin (Abcam, Cambridge, MA) was used as loading control. Bound antibodies were detected with peroxidase-labeled goat antirabbit IgG-horseradish peroxidase (HRP) secondary antibody (Santa Cruz Biotechnology, Dallas, TX).

In Vitro Silencing of UCA1

In vitro silencing of lncRNA UCA1 was performed in the HCT-116 cell line through antisense oligonucleotide (ASO) technology. Cells (5 × 104 per well) were seeded in 24-well plates and simultaneously reverse transfected with 100 nM UCA1 specific Antisense LNA GapmeR (Exiqon, Vedbaek, Denmark) and Lipofectamine RNAiMAX Reagent (Invitrogen), according to the manufacturers’ instructions. Cells treated with a scramble molecule were used as controls. All experiments were performed in biological triplicates. Cells were harvested at 24 and 48 hr after transfection and lysed with Trizol. Specific PCR primers were designed and used to evaluate the expression of UCA1 target mRNAs by Real-Time PCR (Table S9).

Computational Analysis

The results obtained from CRC biopsies were validated in independent cohorts represented by datasets of microarray experiments from GEO DataSets (https://www.ncbi.nlm.nih.gov/gds/). We reanalyzed the expression of DE lncRNAs in 14 CRC datasets; also, results obtained in serum exosomes were confirmed by analyzing a dataset on peripheral blood samples from CRC patients. All analyzed datasets are listed in Table S10. Potential TFBSs harbored on the UCA1 promoter and upstream regulatory region (1 kb) were retrieved from ENCODE tracks (i.e., DNase hypersensitive and chromatin immunoprecipitation sequencing [ChIP-seq] on TFs and H3K27Ac, H3K4me1, and H3K4me3 histone marks) mapped on the UCSC Genome Browser (https://genome.ucsc.edu/). The existence of a positive or negative correlation of expression between the identified TFs and UCA1 was evaluated in expression datasets showing increased expression of UCA1 in CRC tumor tissues compared to normal mucosa.

Prediction of miRNAs sponged by UCA1 was performed through miRcode (http://www.mircode.org/). MiRNA targets experimentally validated in CRC models were retrieved by StarBase (http://starbase.sysu.edu.cn/), miRTarbase (http://mirtarbase.mbc.nctu.edu.tw/php/index.php), and literature data mining; their expression in CRC was investigated in CRC datasets by Oncomine (https://www.oncomine.org/). MiRNAs showing binding sites within the UCA1 sequence and acting as tumor suppressors in CRC were selected for subsequent analysis. Among the experimentally validated targets of selected miRNAs, those characterized by a positive correlation of expression with UCA1 were selected (see also the Results). Correlation analysis (Pearson or Spearman correlation, p-value ≤ 0.05) was performed on GEO and TCGA datasets. RNA-RNA interactions were predicted by using the IntaRNA tool (http://rna.informatik.uni-freiburg.de/IntaRNA/Input.jsp). MiRNAs binding the 3′-UTRs of selected targets were predicted by microrna.org (http://34.236.212.39/microrna/home.do). A molecular network was built and visualized by using the BisoGenet plug-in of Cytoscape 3.6 (http://www.cytoscape.org).78, 79 The FatiGO tool (http://v4.babelomics.org) was applied to determine the statistical over-representation (Fisher’s exact test, p-value < 0.05) of Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Reactome, and BioCarta terms within the molecular network.

Author Contributions

M.P., M.R., C.D.P., F.B., A. Cappellani, and A.B. designed and conceived the experiments. R.C. and A.M. obtained and characterized biological samples from patients. C.B. performed the experiments. C.B., D. Brex, A. Caponnetto, M.C., M.S., and D. Barbagallo contributed to the acquisition, analysis, and interpretation of data. M.P., C.B., and M.R. wrote the paper. All authors read and approved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.