Long non-coding RNA CCL14-AS suppresses invasiveness and lymph node metastasis of colorectal cancer cells by regulating MEP1A
Cancer Cell International volume 23, Article number: 27 (2023)
Long non-coding RNAs (lncRNAs) play important roles in the biology of colorectal cancer (CRC). There are several lncRNAs associated with invasion and metastasis have been characterized in CRC. However, studies focusing on the precise molecular mechanisms by which lncRNAs function in lymph node (LN) metastasis in CRC are still limited.
In this study, by analyzing TCGA dataset, we identified that AC244100.2 (termed CCL14-AS), a novel lncRNA enriched in the cytoplasm, was negatively correlated with LN metastasis and unfavorable prognosis of CRC. In situ hybridization was used to examine CCL14-AS expression in clinical CRC tissues. Various functional experiments including migration assay and wound-healing assay were used to investigate the effects of CCL14-AS on CRC cells migration. The nude mice popliteal lymph node metastasis model assay further confirmed the effects of CCL14-AS in vivo.
CCL14-AS expression was significantly downregulated in CRC tissues compared to adjacent normal tissues. In addition, low CCL14-AS expression was correlated with advanced T classification, LN metastasis, distant metastasis, and shorter disease-free survival of CRC patients. Functionally, CCL14-AS overexpression inhibited the invasiveness of CRC cells in vitro and LN metastasis in nude mice. On the contrary, knockdown of CCL14-AS promoted the invasiveness and LN metastasis abilities of CRC cells. Mechanistically, CCL14-AS downregulated the expression of MEP1A via interacting with MEP1A mRNA and reduced its stability. Overexpression of MEP1A rescued the invasiveness and LN metastasis abilities in CCL14-AS-overexpressing CRC cells. Moreover, the expression levels of CCL14-AS was negatively correlated with that of MEP1A in CRC tissues.
We identified a novel lncRNA, CCL14-AS, as a potential tumor suppressor in CRC. Our findings supported a model in which the CCL14-AS/MEP1A axis serves as critical regulator in CRC progression, suggesting a novel biomarker and therapeutic target in advanced CRC.
Colorectal cancer (CRC), the third most diagnosed cancer and the second leading cause of cancer-associated mortality worldwide, accounts for 10% of all new cancer cases diagnosis and 9.4% of all cancer deaths in 2020 . Approximately 20% CRC patients are with stage IV disease which is present with metastasis when diagnosed  and over half of patients with advanced CRC die from recurrence and metastasis. Early diagnosis and treatment of patients with invasive CRC tumor in early stages is vital to reduce the mortality of CRC patients . However, the potential risk of lymph node (LN) metastasis is still an important clinical consideration, especially in patients with early-stage CRC . LN metastasis occurs in 6.8–19.6% early-stage CRC patients [5, 6]. In CRC patients without distant metastasis, LN metastasis is a significant factor for predicting survival. In addition, regional LN metastasis is an essential step in tumor cell dissemination . Therefore, it is of great value to identify novel biomarkers and potential targets for LN metastasis in CRC.
Long non-coding RNAs (lncRNAs) are defined as transcripts longer than 200 nucleotides that have no or limited protein-coding capacity [8, 9]. With the great progress in understanding the role of lncRNAs, increasing evidence have shown that lncRNAs are key regulators of cancer pathways and as biomarkers of diseases . LncRNAs may contribute to carcinogenesis and metastasis in cancer by regulating protein-coding genes expression at both the transcriptional and post-transcriptional levels [10,11,12]. For examples, lncRNA H19, HOTAIR, and MALAT1 were widely reported to be involved in carcinogenesis and metastasis of human cancers[13,14,15,16,17,18,19,20,21]. LncRNA H19 has been reported to play an oncogenic role by promoting cell proliferation, motility, invasion, and metastasis, in many types of cancers [13,14,15,16,17]. LncRNA HOTAIR is upregulated in multiple cancers and is a powerful predictor of tumor aggressiveness and metastasis [3, 18,19,20]. HOTAIR promotes invasion and metastasis in a manner dependent on Polycomb repressive complex 2 (PRC2) . MALAT1 suppresses breast cancer metastasis by binding and inactivating the pro-metastatic transcription factor TEAD . LncRNAs, such as HOXD-AS1, LINC00460, ZEB1-AS1, MIR17HG, have been demonstrated to have essential roles in regulating the proliferation, migration, invasion, and metastasis of CRC cells [23,24,25,26]. However, studies focusing on the precise molecular mechanisms by which lncRNAs function in LN metastatic CRC are still limited and would be critical for exploring novel biomarkers and therapeutic targets.
In this context, by analyzing TCGA dataset, we identified a novel lncRNA AC244100.2 that was negatively correlated with LN metastasis and unfavorable prognosis of CRC. AC244100.2 is an antisense lncRNA of C–C motif chemokine ligand 14 (CCL14) and overlaps with the 5ʹ UTR and first coding exon of the CCL14 mRNA. Thus, we termed it as CCL14-AS. CCL14-AS is located at human chromosome 17q12 and is transcribed into a 3878 bp transcript. The function and underlying mechanism of CCL14-AS in cancer progression are unclear. In the current study, we sought to figure out the expression pattern and biological role of CCL14-AS in the progression of CRC. We found that CCL14-AS was enriched in the cytoplasm of CRC cells. Overexpression of CCL14-AS significantly weakened migration in vitro and LN metastasis ability of CRC cells in nude mice. Moreover, CCL14-AS downregulated the expression of MEP1A via interacting with MEP1A mRNA and reduced its stability. Overexpression of MEP1A rescued the migration and LN metastasis abilities in CCL14-AS-overexpressing CRC cells. Our data demonstrated that CCL14-AS acted as a potential tumor suppressor in CRC. CCL14-AS/MEP1A axis may serve as promising regulator in CRC progression, suggesting a novel biomarker and therapeutic target in advanced CRC.
Materials and methods
Patients and specimens
The tissue microarrays containing a total of 96 CRC samples were collected and diagnosed in the Department of Pathology, Nanfang Hospital Southern Medical University. None of these patients received radiotherapy and chemotherapy before surgery. The clinical features and pathological information were obtained retrospectively from the medical records of these patients. In addition, eight pairs of fresh tumor tissues and non-cancerous adjacent tissues from patients with CRC were obtained from the operating room of Nanfang Hospital. The fresh biopsies were stored in liquid nitrogen before use. Written-informed consent was obtained from all participants prior to sample collection and the use of clinical materials for research purposes were approved by the Southern Medical University Institutional Board (Guangzhou, China).
Cell lines and cell culture
The human CRC cell lines and human immortalized colon epithelial cell line NCM460 cells were purchased from the American Type Culture Collection and stored in the Department of Pathology, Nanfang Hospital Southern Medical University. RKO, SW837, DLD1, COLO205 and LoVo cells were cultured in RPMI 1640 medium (Gibco) with 10% fetal bovine serum (FBS) (Gibco), 100 U/mL penicillin and 100 μg/ml streptomycin. SW480 and SW620 cells were cultured in Leibovitz's L-15 medium (Gibco) with 10% FBS (Gibco), 100 U/ml penicillin and 100 μg/ml streptomycin. NCM460 and HCT116 cells were cultured in McCoy’s 5A medium (Gibco) supplemented with 10% FBS (Gibco),100 U/mL penicillin and 100 μg/mL streptomycin. Cells were cultured at a 37 °C humidified incubator with 5% CO2 and periodically tested mycoplasma negative.
Lentivirus production and construction of stable cell lines
To evaluate the expression of CCL14-AS, CCL14-AS constructs were generated by cloning full-length human CCL14-AS cDNA into Ubi-MCS-SV40-puromycin vector. For deletion of CCL14-AS, two short hairpin RNA (shRNA) (sequences: CCL14-AS shRNA#1: 5ʹ-GCTCTGGAATTAGCCATTT-3ʹ; CCL14-AS shRNA#2: 5ʹ-GCTCCAACATAATGACAAA-3ʹ) were cloned into hU6-MCS-Ubiquitin-EGFP-IRES-puromycin. The cDNA encoding the ORF of MEP1A was amplified by PCR and inserted into pLVX-IRES-Neo lentivirus vector. The following primers were used for MEP1A cloning: Forward primer, 5ʹ-CCGGTGAATTCCTCGAGCCACCATGGCTTGGATTAGATCCACTTG-3ʹ; Reverse primer: 5ʹ-GGCGGGATCCTACTTGTCATCGTCGTCCTTGTAATCCTTCCTTGGCCTTTGGGAAAGG-3ʹ. Lentivirus were generated by co-transfecting 293 T cells with CCL14-AS or MEP1A-overexpressing plasmids and two packing vectors (psPAX and pMD2G) using Lipofectamine 2000 (Invitrogen). The supernatants of 293 T cultured medium after 48 h, filtering through 0.45-mm filters (Millipore, Temecula, CA, USA) and concentrated for further using. The expression of CCL14-AS in stable cell lines was detected by qRT-PCR. The expression of MEP1A was detected by qRT-PCR and western blot.
RNA isolation, reverse transcription (RT) and real-time PCR
RNA isolation was carried out as previously described . qRT-PCR was performed at least three times in triplicate using SYBR GREEN MIXTURE (Vazyme, Q311-03) and CFX96 Touch Real-Time PCR detection system. The data were normalized to the geometric mean of the housekeeping gene GAPDH and calculated using the 2−ΔΔCT method. The following primer sequences were used for amplification: CCL14-AS Forward primer, 5ʹ- AGCATGAGTGGTCTTTAATTCAAA -3ʹ; CCL14-AS reverse, 5ʹ-GGCTGCCATTCCCTTCTT-3ʹ; MEP1A Forwardprimer, 5ʹ-TGACAGCACAGGCAATGTTCGC-3ʹ; MEP1A reverse, 5ʹ-GTCGCCTTTTGTGCCCTGGAAA-3ʹ; GAPDH Forward primer, 5ʹ-AACGGGAAGCTTGTCATCAA-3ʹ; GAPDH reverse, 5ʹ-TGGACTCCACGACGTACTCA-3ʹ.
Western blot was performed as described previously  using anti-MEP1A (Bioss, bs-6056R) and anti-GFP (Abcam, ab290). Anti-GAPDH and anti-α-tubulin antibody (Sigma, #9026) were used as the internal control.
Immunohistochemical (IHC) staining was performed as previously described using a specific MEP1A antibody (Bioss, bs-6056R) . After reviewed and scored independently by two observers, MEP1A protein expression levels were quantified based on the intensity of staining (0–3). The staining intensity was graded according to the following criteria: 0 (no staining); 1 (weak staining = light yellow), 2 (moderate staining = yellow brown), and 3 (strong staining = brown). Determine the best cut-off value: intensity of staining ≥ 2 is used to define tumors with high MEP1A expression, and index ≤ 1 is used to define tumors with low MEP1A expression.
RNA fluorescence in situ hybridization (FISH)
The probe of CCL14-AS was designed and synthesized by GenePharma (Shanghai, China). The sequence is: 5′-CGTTCTTTGGATCGTGAGGATGATTTGTCAGCTCAGAGGAACAACTAGGGAAAGGCTTCTATCAAATACCAACTGTCAAAAGCTTAGGGGAAAAGTACAGGAACTTCC-3′. The hybridization and fluorescence detection were performed using FISH Kits from GenePharma (F22201/100, Shanghai, China) according to the manufacturer’s protocol. Finally, the signal was observed and photographed under a confocal microscope or fluorescent microscopy (Olympus, Japan). ImageJ software was used to measure the fluorescence intensity and distribution. The average of intensity and standard deviation were calculated and divided into low, medium, and high groups.
Cellular proliferation, colony formation assays
For cell proliferation assays, a total of 3000 cells/well were plated into 96-well plates and cultured for 24 h, then added 20 µl of MTT in each well and incubated at 37 °C for 4 h. After removing the supernatant, in each well, we added 150 µl of DMSO and incubated at 37 °C for 10 min. Finally, we examined the absorbance at 570 nm with a microplate reader (800TS, Biotek, USA). For the clone-formation experiment, 1000 cells/well were seeded into 6-well plates and cultured for 2 weeks. Clones were fixed in 4% paraformaldehyde (PFA) and stained with 0.5% crystal violet. colonies were counted manually and each single colony consists of at least 50 cells.
Cell migration and wound healing assays
Cell migration assays were performed using transwell. Cells (1 × 105) in culture medium containing 1% FBS were added to the top chamber of 24-well transwell chambers plates (8.0 μm, Corning, NY, USA.) and culture medium with 10% FBS was added in the lower chamber as a chemoattractant.
After incubation for 24 h, the underside cells of the membrane were fixed, stained with 0.1% Crystal Violet, and then the percentages of migrated cells were counted (10 random 200 × fields per well). Three independent experiments were performed, and the data were presented as the mean ± SEM. For wound healing assays, a total of 1 × 106 cells/well were plated into 6-well plates and cultured to about 80% confluence. After serum starvation for 24 h, scratches were performed in the middle slides using 10 μl pipette tip for each well. The closure of the gap distance was observed and photographed at 0 h, 12 h, 24 h, 36 h and 48 h under an inverted light microscope and quantitatively evaluated using ImageJ. Each experiment was repeated at least three times.
Isolation of cytoplasmic and nuclear
The nuclear and cytosolic fractions of LoVo cells were isolated with Nuclear/Cytosol Extraction Kit (Bestbio, China) following the manufacturer’s protocol.
RNA sequencing analysis
Total RNA was isolated from LoVo stable cells expressing vector or CCL14-AS. The sequencing libraries were generated using NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, USA) following manufacturer’s recommendations and index codes were added to attribute sequences to each sample. The samples were paired end sequenced with a read length of 150 bp on an Illumina Novaseq platform. The process of sequencing was controlled by Illumina Data collection software. The library construction and sequencing were performed at GenePharma (Shanghai, China). After pre-treatment of the raw reads (filtering and QC) and alignment to the genome, the expression level of genes was determined based on the value of FPKM (Fragments Per Kilobase of transcript per Million mapped reads). Differential expression analysis of two groups was performed using the DESeq2 R package. Genes with an adjusted P-value < 0.05 found by DESeq2 were assigned as differentially expressed. Gene Ontology (GO) enrichment analysis of differentially expressed genes was implemented by using the DAVID GO database to search for enriched pathways.
RNA immunoprecipitation (RIP)
pMS2-GFP (Addgene) was transfected into stably overexpress CCL14-AS cells. After 48 h, cells were used to perform RIP experiments using a GFP antibody (5 µg per reaction, Abcam, ab290). LoVo cells were lysed in lysis Buffer with RNase inhibitor (Beyotime, 1:1000), DMSF (1:100) and Protease inhibitor (Beyotime,1:100). The protein A/G sepharose beads were washed by lysis Buffer, and the cell extracts were incubated with these beads bound to GFP antibody or control IgG at 4 °C for 2 h. The beads were washed and divided two parts. One part was incubated with proteinase K to remove protein and performed RNA isolation. Other part was used for western blot. Finally, qRT-PCR analysis was performed on the purified RNA, and western blot analysis for Anti-IgG and anti-GFP are performed to detect the quality of the experiment.
Cells were treated with the transcriptional inhibitor actinomycin D (Sigma, USA) for 0 min, 20 min, 40 min, 1 h, 2 h, 4 h, 6 h, 8 h, 10 h and 12 h, respectively. Cells were harvested at different time points, and RNA was isolated. Then, qRT-PCR analysis was used to detect the level of MEAP1 mRNA as mentioned above.
Mouse popliteal lymphatic metastasis model
The female BALB/c nude mice (4 to 6-week-old) were obtained from the Animal Center of Southern Medical University, Guangzhou, China. All animal experiments were approved by the Animal Care and Use Committee at Southern Medical University. Groups of LoVo-Vector, LoVo-CCL14-AS, LoVo-CCL14-AS-MEP1A, HCT116-Scramble and HCT116-shCCL14-AS cells (5 × 106 per mouse) were injected into the footpads of mice correspondingly. After 8 weeks, the mice were sacrificed, and their footpad tumors and inguinal lymph nodes were detached. These tissues were fixed with formaldehyde and paraffin-embedded, and then 5-mm sections were cut and stained with haematoxylin and eosin (H&E).
All statistical analyses were conducted by using SPSS version 20.0 (SPSS Inc, IL, USA) and GraphPad Prism Version 8.0 (GraphPadSoftware, Inc.USA). The differential expression between the two groups was calculated using Student’s t test. One-way ANOVA test was used to compare the differences of data from multiple groups, and LSD-T was used for multiple comparisons. The relationship between pathological parameters was compared using two independent samples nonparametric test and a chi-square test. The Kaplan–Meier method was used for analyzing the survival curves of CRC patients in low- and high-CCL14-AS or MEP1A expression groups and the log-rank test log-rank test was used to compare differences. A p value of < 0.05 was considered significant. (* p < 0.05, ** p < 0.01, *** p < 0.001).
CCL14-AS is downregulated in CRC and associated with favorable prognosis
To identify novel biomarkers which can effectively predict CRC metastasis, The Cancer Genome Atlas (TCGA) datasets with expression profiles of lncRNAs in exoRBase CRC data were analyzed. The relative expressions of lncRNAs were analyzed by comparing cancer samples with normal samples, and metastatic (M1) samples with non-metastatic (M0) samples. Fold change ≥ 0.5 and p< 0.01 were considered differently expressed genes. LncRNAs differentially expressed in CRC tissues compared with normal tissues (Additional file 4: Table S1, cancer-related lncRNAs) and differently expressed in LN negative samples compared with LN positive samples (Additional file 5: Table S2, LN-related lncRNAs) were identified. Univariate Cox proportional-hazards regression model were performed to analyze the lncRNA that were associated with patient survival (Additional file 6: Table S3, survival-related lncRNAs). Integrating these three groups of lncRNAs triangulated on a single lncRNA, AC244100.2 (Fig. 1A). As AC244100.2 is an antisense lncRNA of C–C motif chemokine ligand 14 (CCL14) and overlaps with the 5ʹ UTR and first coding exon of the CCL14 mRNA, we termed it as CCL14-AS (Additional file 1: Fig S1A). The coding potential of CCL14-AS was calculated using the National Center for Biotechnology Information (NCBI) ORF finder, showing that CCL14-AS fails to translate protein of more than 85 amino acids (Additional file 7: Table S4). Analyses using the Coding Potential Assessment Tool (CPAT) showed that protein coding probability of CCL14-AS is 2.69% (Additional file 1: Figure S1B). This result suggested that CCL14-AS is a noncoding RNA. The secondary structure of CCL14-AS was predicted using online software RNAfold (Additional file 1: Fig S1C).
To evaluate the potential role CCL14-AS in CRC, the expression pattern of CCL14-AS and its clinical relevance were analyzed in TCGA COADREAD dataset and clinical CRC samples. In the TCGA COADREAD dataset, CCL14-AS was significantly downregulated in CRC tissues as compared with normal tissues (Fig. 1B, C). In addition, low CCL14-AS expression level is significantly correlated with lymph node metastasis (Fig. 1D) and poor Overall survival rate (OS) in TCGA dataset (Fig. 1E) in TCGA CRC patients. Next, the expression of CCL14-AS was validated in 8 paired clinical CRC samples by qRT-PCR, showing that CCL14-AS was dramatically downregulated in CRC tissues relative to adjacent non-cancerous tissues (Fig. 1F). CCL14-AS expression was further examined via in situ hybridization (ISH) in an expanded cohort of 96 paraffin-embedded CRC samples. The results showed that CCL14-AS is highly expressed in normal intestine tissues. However, it was dramatically downregulated in CRC tissues (Fig. 1G). Moreover, low CCL14-AS expression was significantly correlated with advanced T classification, lymph node metastasis, and distant metastasis (Fig. 1H). Furthermore, Kaplan–Meier analysis showed that low CCL14-AS expression was significantly associated with poor overall survival rate (OS) in this cohort of CRC samples (Fig. 1I). These results showed that CCL14-AS is a novel lncRNA associated with CRC progression.
CCL14-AS suppresses invasiveness and LN metastasis of CRC cells
The endogenous expression of CCL14-AS in CRC cell lines were examined by qRT-PCR. The results showed that the CCL14-AS expression level was much higher in the human immortalized colon epithelial cell NCM460 than in CRC cells (Fig. 2A). To investigate the biological functions of CCL14-AS on the progression of CRC, CCL14-AS was stably overexpressed in LoVo and SW620 cell lines, which showed low endogenous expression of CCL14-AS (Fig. 2B). MTT assays and colony formation assays revealed that overexpression of CCL14-AS does not affect the proliferation of CRC cells (Additional file 2: Figure S2A–S2B). However, wound healing assays and migration assays showed that overexpression of CCL14-AS dramatically decreased the migration and invasion of CRC cells (Fig. 2C, D). As the CCL14-AS expression levels are closely correlated with LN metastasis both in TCGA COAD datasets and in clinical CRC samples, we examined the role of CCL14-AS in LN metastasis potential using the nude mice popliteal lymph node metastasis model. Expectedly, overexpression of CCL14-AS significantly inhibited LN metastasis, as determined the volumes and numbers of the metastatic lymph nodes (Fig. 2E, F).
Next, we validated the role of CCL14-AS in invasiveness and LN metastasis by knocking down endogenous CCL14-AS in HCT116 and RKO CRC cell lines, which showed relatively high expression of CCL14-AS (Fig. 3A). MTT assays and colony formation assays showed knockdown of CCL14-AS has no effects on the proliferation of both RKO and HCT116 cells (Additional file 2: Figure S2C-S2D). However, wound healing assays and migration assays showed that knockdown of CCL14-AS dramatically increased the migration and invasion of CRC cells (Fig. 3B–D). Furthermore, knockdown of CCL14-AS significantly promoted the LN metastasis in nude mice (Fig. 3E, F). Taken together, these data showed that lncRNA CCL14-AS suppresses invasiveness and LN metastasis of CRC.
CCL14-AS interacts with MEP1A mRNA and reduces its stability
Next, we sought to unveil the mechanism by which CCL14-AS suppressed invasiveness and LN metastasis. The subcellular location of CCL14-AS in CRC cells was examined by FISH, showing that CCL14-AS was mainly localized in the cytoplasm of SW620 cells (Fig. 4A). Overexpression of CCL14-AS in SW620 cells further increased the amount of cytoplasm enriched CCL14-AS (Fig. 4A). Nuclear and cytoplasmic RNA fractions were prepared from LoVo cells to confirm the subcellular localization of CCL14-AS. The results also demonstrated that CCL14-AS was mainly enriched in the cytoplasm (Fig. 4B, C). These findings were further confirmed by online software LncLocator analyses (Fig. 4D), indicating that CCL14-AS is a cytoplasic-enriched lncRNA in CRC cells.
LncRNAs located in the cytoplasm can mediate signal transduction pathways and translational programs by posttranscriptional control of gene expression . For examples, cytoplasmic lncRNAs can modulate protein posttranslational modifications , or regulate mRNA translation and stability [31, 32]. To explore the potential downstream target genes of CCL14-AS, RNA-seq was performed. Differentially expressed genes (DEGs) analysis revealed a series of genes with dramatic difference in their expression. GO analysis showed that the DEGs were significantly enriched in RNA catabolic processes, cell-substrate adherent junction, extracellular matrix binding, regulation of translation initiation, ameboidal-type cell migration, and angiogenesis (Fig. 4E). Among these DEGs genes, LYVE1 (Lymphatic Vessel Endothelial Hyaluronan Receptor 1) and MEP1A (Meprin A Subunit Alpha) were the top 2 downregulated genes, while KRT40 (Keratin 40) and POMC (Pro-opiomelanocortin) were the top 2 upregulated genes upon CCL14-AS overexpression (Fig. 4F, G). Software IntaRNA  was applied to predict whether these 4 genes could potentially interact with CCL14-AS. This in silico analysis revealed that CCL14-AS was predicted to interact with MEP1A mRNA and LYVE1 mRNA (Additional file 3: Fig S3). Since the LYVE1 is mainly expressed in endothelial cells and immune cells (THE HUMAN PROTEIN ATLAS), we focus on MEP1A for further investigation in the current study. Consistent with the RNA-seq results, qRT-PCR displayed that overexpression of CCL14-AS significantly decreased, while knockdown of CCL14-AS increased the expression of MEP1A, at both mRNA and protein levels (Fig. 4H–J).
To further elucidate the mechanisms regulating MEP1A expression by CCL14-AS in CRC, we performed RNA immunoprecipitation (RIP) experiments. MS2-RIP experiments revealed that CCL14-AS interacted with MEP1A mRNA (Fig. 4K). RNA-stability analysis showed that CCL14-AS significantly promoted the degradation of MEP1A mRNA (Fig. 4L). Altogether, these data suggested MEP1A is a target mRNA of CCL14-AS in the cytoplasm of CRC cells.
CCL14-AS suppresses CRC migration and LN metastasis by regulating MEP1A
It has been reported that MEP1A is upregulated in CRC. Overexpression of MEP1A promotes migration and metastasis ability of CRC cells . We hypothesized that MEP1A might be a key CCL14-AS downstream regulator that promotes progression of CRC. To confirm this hypothesis, MEP1A were overexpressed in CCL14-AS overexpressing cells (Fig. 5A-B). Wound-healing assays and migration assays showed that the ectopic expression of MEP1A significantly reduced the repressive effects on migration and invasion of CRC cells mediated by CCL14-AS in vitro (Fig. 5C–E). Popliteal lymph node metastasis model demonstrated that ectopic expression of MEP1A significantly reduces the repressive effects of lymph node metastasis mediated by CCL14-AS (Fig. 5F, G).
MEP1A is up-regulated in CRC and negative correlates with CCL14-AS
To assess a potential link between CCL14-AS and MEP1A expression in human CRC, we analyzed the expression levels of CCL14-AS and MEP1A in 8 pairs of matched CRC samples and adjacent noncancerous tissues. The results showed that both the mRNA level and protein level of MEP1A are significantly higher in CRC tissues than in the adjacent tissues (Fig. 6A). Statistical analyses showed that high MEP1A mRNA expression was negatively correlated with CCL14-AS expression in these samples (Fig. 6B). By analyzing the MEP1A protein expression levels using IHC, we found that MEP1A was significantly elevated in CRC tissues as compared in adjacent non-cancerous tissues, which is in consistent with a previous study  (Fig. 6C). Moreover, the analysis of CRC tissue microarray containing 96 cases of CRC tissues revealed the significant negative correlation between CCL4-AS expression and MEP1A expression (Fig. 6D).
Numerous studies have revealed that upregulation or downregulation of lncRNA expression has been found in CRC patients. These dysregulated lncRNAs may serve as oncogenes or tumor suppressors. For example, lnc-GNAT1-1 was downregulated in CRC. Low lnc-GNAT1-1 expression was associated with unfavorable clinicopathological features and a poor survival of CRC patients. Knockdown of lnc-GNAT1-1 inhibited the aggressive phenotypes of CRC cell lines and their ability of metastasis to the liver . LINC01082 was downregulated in CRC tissues. Ectopic expression of LINC01082 suppressed the proliferation, migration, and invasion of CRC cells . On the contrary, there are other lncRNAs function as oncogenes in CRC. For example, CCAT1 and CCAT2 expression were upregulated in CRC. CCAT1 alone or in combination with CCAT2 associates with poor survival in CRC patients, and serves as an important prognostic biomarker in CRC . NEAT1 expression is upregulated in CRC and correlates with advanced disease stages and poor overall survival. Overexpression of NEAT1 promotes cell viability and metastasis ability of CRC cells . Studies focusing on the precise molecular mechanisms by which lncRNAs function in LN metastatic CRC are still limited.
In the current study, we identified CCL14-AS,a novel lncRNA, that is associated with disease progression and metastasis in CRC.CCL14-AS expression is downregulated in CRC tissues as compared with adjacent normal intestine tissues. Low CCL14-AS expression was associated with advanced T classification, LN metastasis, distant metastasis, and shorter disease-free survival of CRC patients. Such novel data suggest that the levels of CCL14-AS may be valuable for prognosis of CRC if validated. In addition, CCL14-AS overexpression significantly weakened the invasiveness and LN metastasis ability of CRC cells. These data supported that CCL14-AS may serve as a potential tumor suppressor and a biomarker for LN metastasis in CRC. Moreover, CCL14-AS downregulated MEP1A expression through decreasing its mRNA stability. Furthermore, overexpression MEP1A rescued the CCL14-AS mediated suppression of invasiveness and LN metastasis of CRC cells. Accordingly, CCL14-AS may be a new therapeutic target and interventional strategy to descrease MEP1A may be valuable for control of CRC.
As a novel lncRNA, the functional role and underlying mechanisms of CCL14-AS are unclear. To address this question, RNA-seq was performed to identify the pathways that may affected by CCL14-AS in CRC cells. Interestingly, the DEGs by CCL14-AS overexpression were significantly enriched in RNA catabolic processes, cell-substrate adherent junction, extracellular matrix binding, regulation of translation initiation, ameboidal-type cell migration, and angiogenesis. These results strongly suggested a role of CCL14-AS in invasion and metastasis of tumor cells. Further analysis indicated that LYVE1 and MEP1A might be downstream effectors of CCL14-AS. As a specific marker for lymphatic endothelial cells [42, 43], LYVE1 could be an interesting candidate of CCL14-AS downstream target that mediate invasion and LN metastasis. In addition, it has been reported that high LYVE-1-positive lymphatic vessel numbers were associated with poor prognosis in breast cancer . Since LYVE1 was mainly localized in the endothelial cells while not in epithelial cells, we considered that it might not be a major functional effector downstream of CCL14-AS in CRC cells. However, the regulation of LYVE1 by CCL14-AS could exist in other types of cells. MEP1A is a membrane-bound and secreted astacin metalloproteinases belongs to the metzincin superfamily . Abnormal MEP1A expression has been implicated in several cancers, including CRC [34, 46,47,48], breast cancer , and hepatocellular carcinoma . As a protease that is known to cleave extracellular matrix components in vitro, MEP1A may contribute to tumor progression by facilitating migration, invasion, and metastasis of tumor cells. In CRC, MEP1A activated the EGFR signaling by promoting the shedding of the ligand, thereby activating the ERK/ZEB1 pathway to promote proliferation and migration of CRC cells . These studies suggested that MEP1A is an important regulator that involves in CRC invasion and metastasis. However, the mechanisms that mediated the high expression level of MEP1A in CRC is unclear. In the current study, we showed that CCL14-AS interacted with MEP1A mRNA and promoted its degradation, leading to downregulation of MEP1A. The silence or downregulation of CCL14-AS increased the stability of MEP1A mRNA, resulting in upregulation of MEP1A in CRC. Fascinatedly, CCL14-AS have no effect on CRC cells proliferation. This result may show that CCL14-AS/MEP1A is not the only way to regulate proliferation of CRC cells. There may be the regulatory network to regulate proliferation of CRC cells.
The subcellular localization of lncRNA is closely related to its biological function . Cytoplasm-enriched lncRNAs are widely involved in mediating signal transduction pathways and translational programs . They have essential roles in regulating mRNA translation [31, 32, 52] and in modulating protein posttranslational modifications [30, 53]. They may also serve as competing endogenous RNAs  or precursors of microRNAs . There are a few cytoplasm-enriched lncRNAs that target mRNA transcripts and modulate their stability. For example, lncRNA BACE1-AS increases BACE1 mRNA stability in Alzheimer's disease . Overexpression of 1/2-sbsRNAs promotes the degradation of the target mRNAs, while depletion of 1/2-sbsRNAs led to the stabilization and upregulation of the target mRNAs [32, 57]. In the present study, we demonstrated that CCL14-AS was enriched in the cytoplasm of CRC cells. MS2-RIP assays revealed a direct interaction between CCL14-AS and MEP1A mRNA. Overexpression of CCL14-AS lead to increased degradation and downregulation of MEP1A. Interestingly, found that double‑stranded RNA‑binding protein Staufen homolog 1 (STAU1) is one of CCL14-AS binding proteins were predicted by catRAPID. STAU1-mediated mRNA decay (SMD) is a common mRNA degradation process in mammalian cells. 1/2-sbsRNAs mediate mRNA decay by recruiting STAU1 to target mRNAs . Thus, we speculated that CCL14-AS may promote MEP1A mRNA decay through recruiting STAU1 to MEP1A mRNA.
In summary, we identify a novel lncRNA, CCL14-AS, involved in the progression of CRC. CCL14-AS is downregulated in CRC and inhibits the invasiveness and LN metastasis in CRC. We also provide evidence that the metalloprotease MEP1A may represent a downstream effector of CCL14-AS in the cancer cells. CCL14-AS is mainly localized in the cytoplasm of CRC cells and interacts with MEP1A mRNA to promote the degradation of MEP1A mRNA. Our results suggested that CCL14-AS may be a potential diagnostic biomarker and therapeutic target for CRC (Fig. 6E).
Availability of data and materials
The datasets generated or analyzed for this study will be made available by the corresponding author on reasonable request.
C–C motif chemokine ligand 14
Coding potential assessment tool
Haematoxylin and eosin
Long non-coding RNAs
In situ hybridization
Lymphatic vessel endothelial hyaluronan receptor 1
Meprin A subunit alpha
The national center for biotechnology information
Overall survival rate
Polycomb repressive complex 2
Short hairpin RNA
Staufen homolog 1
STAU1-mediated mRNA decay
The cancer genome atlas
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 Countries. CA Cancer J Clin. 2021;71:209–49.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2015. CA Cancer J Clin. 2015;65:5–29.
Qu X, Alsager S, Zhuo Y, Shan B. HOX transcript antisense RNA (HOTAIR) in cancer. Cancer Lett. 2019;454:90–7.
Yasue C, Chino A, Takamatsu M, Namikawa K, Ide D, Saito S, et al. Pathological risk factors and predictive endoscopic factors for lymph node metastasis of T1 colorectal cancer: a single-center study of 846 lesions. J Gastroenterol. 2019;54:708–17.
Lee YJ, Huh JW, Shin JK, Park YA, Cho YB, Kim HC, et al. Risk factors for lymph node metastasis in early colon cancer. Int J Colorectal Dis. 2020;35:1607–13.
Ha RK, Han KS, Sohn DK, Kim BC, Hong CW, Chang HJ, et al. Histopathologic risk factors for lymph node metastasis in patients with T1 colorectal cancer. Ann Surg Treat Res. 2017;93:266–71.
Pereira ER, Jones D, Jung K, Padera TP. The lymph node microenvironment and its role in the progression of metastatic cancer. Semin Cell Dev Biol. 2015;38:98–105.
Ulitsky I, Bartel DP. lincRNAs: genomics, evolution, and mechanisms. Cell. 2013;154:26–46.
Kopp F, Mendell JT. Functional classification and experimental dissection of long noncoding RNAs. Cell. 2018;172:393–407.
Huarte M. The emerging role of lncRNAs in cancer. Nat Med. 2015;21:1253–61.
Liu SJ, Dang HX, Lim DA, Feng FY, Maher CA. Long noncoding RNAs in cancer metastasis. Nat Rev Cancer. 2021;21:446–60.
Slack FJ, Chinnaiyan AM. The role of non-coding RNAs in oncology. Cell. 2019;179:1033–55.
Hu XT, Xing W, Zhao RS, Tan Y, Wu XF, Ao LQ, et al. HDAC2 inhibits EMT-mediated cancer metastasis by downregulating the long noncoding RNA H19 in colorectal cancer. J Exp Clin Cancer Res. 2020;39:270.
Li A, Mallik S, Luo H, Jia P, Lee DF, Zhao Z. H19, a long non-coding RNA, mediates transcription factors and target genes through interference of microRNAs in pan-cancer. Mol Ther Nucleic Acids. 2020;21:180–91.
Zhang Y, Huang W, Yuan Y, Li J, Wu J, Yu J, et al. Long non-coding RNA H19 promotes colorectal cancer metastasis via binding to hnRNPA2B1. J Exp Clin Cancer Res. 2020;39:141.
Raveh E, Matouk IJ, Gilon M, Hochberg A. The H19 Long non-coding RNA in cancer initiation, progression and metastasis—a proposed unifying theory. Mol Cancer. 2015;14:184.
Vennin C, Spruyt N, Robin YM, Chassat T, Le Bourhis X, Adriaenssens E. The long non-coding RNA 91H increases aggressive phenotype of breast cancer cells and up-regulates H19/IGF2 expression through epigenetic modifications. Cancer Lett. 2017;385:198–206.
Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, Wong DJ, et al. Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010;464:1071–6.
Liu XH, Liu ZL, Sun M, Liu J, Wang ZX, De W. The long non-coding RNA HOTAIR indicates a poor prognosis and promotes metastasis in non-small cell lung cancer. BMC Cancer. 2013;13:464.
Huang KB, Zhang SP, Zhu YJ, Guo CH, Yang M, Liu J, et al. Hotair mediates tumorigenesis through recruiting EZH2 in colorectal cancer. J Cell Biochem. 2019;120:6071–7.
Amodio N, Raimondi L, Juli G, Stamato MA, Caracciolo D, Tagliaferri P, et al. MALAT1: a druggable long non-coding RNA for targeted anti-cancer approaches. J Hematol Oncol. 2018;11:63.
Kim J, Piao HL, Kim BJ, Yao F, Han Z, Wang Y, et al. Long noncoding RNA MALAT1 suppresses breast cancer metastasis. Nat Genet. 2018;50:1705–15.
Yang MH, Zhao L, Wang L, Ou-Yang W, Hu SS, Li WL, et al. Nuclear lncRNA HOXD-AS1 suppresses colorectal carcinoma growth and metastasis via inhibiting HOXD3-induced integrin beta3 transcriptional activating and MAPK/AKT signalling. Mol Cancer. 2019;18:31.
Hou P, Meng S, Li M, Lin T, Chu S, Li Z, et al. LINC00460/DHX9/IGF2BP2 complex promotes colorectal cancer proliferation and metastasis by mediating HMGA1 mRNA stability depending on m6A modification. J Exp Clin Cancer Res. 2021;40:52.
Xu J, Meng Q, Li X, Yang H, Xu J, Gao N, et al. Long noncoding RNA MIR17HG promotes colorectal cancer progression via miR-17-5p. Cancer Res. 2019;79:4882–95.
Zhao S, Guan B, Mi Y, Shi D, Wei P, Gu Y, et al. LncRNA MIR17HG promotes colorectal cancer liver metastasis by mediating a glycolysis-associated positive feedback circuit. Oncogene. 2021;40:4709–24.
Wang S, Qiu J, Liu L, Su C, Qi L, Huang C, et al. CREB5 promotes invasiveness and metastasis in colorectal cancer by directly activating MET. J Exp Clin Cancer Res. 2020;39:168.
Liao WT, Jiang D, Yuan J, Cui YM, Shi XW, Chen CM, et al. HOXB7 as a prognostic factor and mediator of colorectal cancer progression. Clin Cancer Res. 2011;17:3569–78.
Bridges MC, Daulagala AC, Kourtidis A. LNCcation: lncRNA localization and function. J Cell Biol. 2021. https://doi.org/10.1083/jcb.202009045.
Lin A, Li C, Xing Z, Hu Q, Liang K, Han L, et al. The LINK-A lncRNA activates normoxic HIF1alpha signalling in triple-negative breast cancer. Nat Cell Biol. 2016;18:213–24.
Carrieri C, Cimatti L, Biagioli M, Beugnet A, Zucchelli S, Fedele S, et al. Long non-coding antisense RNA controls Uchl1 translation through an embedded SINEB2 repeat. Nature. 2012;491:454–7.
Gong C, Maquat LE. lncRNAs transactivate STAU1-mediated mRNA decay by duplexing with 3ʹ UTRs via Alu elements. Nature. 2011;470:284–8.
Mann M, Wright PR, Backofen R. IntaRNA 2.0: enhanced and customizable prediction of RNA-RNA interactions. Nucleic Acids Res. 2017;45:435–9.
Wang X, Chen J, Wang J, Yu F, Zhao S, Zhang Y, et al. Metalloproteases meprin-a (MEP1A) is a prognostic biomarker and promotes proliferation and invasion of colorectal cancer. BMC Cancer. 2016;16:383.
Ye C, Shen Z, Wang B, Li Y, Li T, Yang Y, et al. A novel long non-coding RNA lnc-GNAT1-1 is low expressed in colorectal cancer and acts as a tumor suppressor through regulating RKIP-NF-kappaB-Snail circuit. J Exp Clin Cancer Res. 2016;35:187.
Xiong W, Qin J, Cai X, Xiong W, Liu Q, Li C, et al. Overexpression LINC01082 suppresses the proliferation, migration and invasion of colon cancer. Mol Cell Biochem. 2019;462:33–40.
Ozawa T, Matsuyama T, Toiyama Y, Takahashi N, Ishikawa T, Uetake H, et al. CCAT1 and CCAT2 long noncoding RNAs, located within the 8q.24.21 ‘gene desert’, serve as important prognostic biomarkers in colorectal cancer. Ann Oncol. 2017;28:1882–8.
Luo Y, Chen JJ, Lv Q, Qin J, Huang YZ, Yu MH, et al. Long non-coding RNA NEAT1 promotes colorectal cancer progression by competitively binding miR-34a with SIRT1 and enhancing the Wnt/beta-catenin signaling pathway. Cancer Lett. 2019;440–441:11–22.
Lin R, Maeda S, Liu C, Karin M, Edgington TS. A large noncoding RNA is a marker for murine hepatocellular carcinomas and a spectrum of human carcinomas. Oncogene. 2007;26:851–8.
Gutschner T, Hammerle M, Eissmann M, Hsu J, Kim Y, Hung G, et al. The noncoding RNA MALAT1 is a critical regulator of the metastasis phenotype of lung cancer cells. Cancer Res. 2013;73:1180–9.
Ji Q, Zhang L, Liu X, Zhou L, Wang W, Han Z, et al. Long non-coding RNA MALAT1 promotes tumour growth and metastasis in colorectal cancer through binding to SFPQ and releasing oncogene PTBP2 from SFPQ/PTBP2 complex. Br J Cancer. 2014;111:736–48.
Banerji S, Ni J, Wang SX, Clasper S, Su J, Tammi R, et al. LYVE-1, a new homologue of the CD44 glycoprotein, is a lymph-specific receptor for hyaluronan. J Cell Biol. 1999;144:789–801.
Jackson DG, Prevo R, Clasper S, Banerji S. LYVE-1, the lymphatic system and tumor lymphangiogenesis. Trends Immunol. 2001;22:317–21.
Bono P, Wasenius VM, Heikkila P, Lundin J, Jackson DG, Joensuu H. High LYVE-1-positive lymphatic vessel numbers are associated with poor outcome in breast cancer. Clin Cancer Res. 2004;10:7144–9.
Sterchi EE, Stocker W, Bond JS. Meprins, membrane-bound and secreted astacin metalloproteinases. Mol Aspects Med. 2008;29:309–28.
Lottaz D, Maurer CA, Noel A, Blacher S, Huguenin M, Nievergelt A, et al. Enhanced activity of meprin-alpha, a pro-migratory and pro-angiogenic protease, in colorectal cancer. PLoS ONE. 2011;6: e26450.
Lottaz D, Maurer CA, Hahn D, Buchler MW, Sterchi EE. Nonpolarized secretion of human meprin alpha in colorectal cancer generates an increased proteolytic potential in the stroma. Cancer Res. 1999;59:1127–33.
Minder P, Bayha E, Becker-Pauly C, Sterchi EE. Meprinalpha transactivates the epidermal growth factor receptor (EGFR) via ligand shedding, thereby enhancing colorectal cancer cell proliferation and migration. J Biol Chem. 2012;287:35201–11.
Matters GL, Manni A, Bond JS. Inhibitors of polyamine biosynthesis decrease the expression of the metalloproteases meprin alpha and MMP-7 in hormone-independent human breast cancer cells. Clin Exp Metastasis. 2005;22:331–9.
OuYang HY, Xu J, Luo J, Zou RH, Chen K, Le Y, et al. MEP1A contributes to tumor progression and predicts poor clinical outcome in human hepatocellular carcinoma. Hepatology. 2016;63:1227–39.
Rashid F, Shah A, Shan G. Long Non-coding RNAs in the Cytoplasm. Genomics Proteomics Bioinform. 2016;14:73–80.
Huarte M, Guttman M, Feldser D, Garber M, Koziol MJ, Kenzelmann-Broz D, et al. A large intergenic noncoding RNA induced by p53 mediates global gene repression in the p53 response. Cell. 2010;142:409–19.
Liu B, Sun L, Liu Q, Gong C, Yao Y, Lv X, et al. A cytoplasmic NF-kappaB interacting long noncoding RNA blocks IkappaB phosphorylation and suppresses breast cancer metastasis. Cancer Cell. 2015;27:370–81.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi P. A ceRNA hypothesis: the rosetta stone of a hidden RNA language? Cell. 2011;146:353–8.
Cai X, Cullen BR. The imprinted H19 noncoding RNA is a primary microRNA precursor. RNA. 2007;13:313–6.
Faghihi MA, Modarresi F, Khalil AM, Wood DE, Sahagan BG, Morgan TE, et al. Expression of a noncoding RNA is elevated in Alzheimer’s disease and drives rapid feed-forward regulation of beta-secretase. Nat Med. 2008;14:723–30.
Park E, Maquat LE. Staufen-mediated mRNA decay. Wiley Interdiscip Rev RNA. 2013;4:423–35.
This work was supported by Grants from the National Natural Science Foundation of China (No. 81872401, 81672886) and the Basic and Applied Basic Research of Guangzhou Basic Research Program (SL2022A04J02033).
Ethics approval and consent to participate
The studies involving the collection of clinical materials were reviewed and approved by The Institutional Committee of Southern Medical University and Nanfang Hospital (Guangzhou, China). The ethics approval statements for animal work were provided by the Committee on the Ethics of Animal Experiments of Southern Medical University.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Figure S1. Bioinformatics analysis of the fundamental properties of CCL14-AS. A Model diagram of CCL14-AS in the genome. B Online software the Coding Potential Assessment Tool (CPAT) predict the Coding ability of CCL14-AS. Blue is noncoding RNA, while red is the coding gene. C The predicted secondary structure (left) and mountain plot (right) representing the MFE (Red), thermodynamic ensemble of RNA (Green), and centroid (Blue) structures of CCL14-AS.
: Figure S2. CCL14-AS had no effect on proliferation of CRC cells. A–B The LoVo and SW620 stable cells overexpressing CCL14-AS were subjected to MTT assays A and Colony-formation assays B. in. C–D The HCT116 and RKO stable cells expressing shRNA against CCL14-AS were subjected to MTT assays C and Colony-formation assays D. in. ns: no significance, * p < 0.05.
: Figure S3. LncRRIsearch database predicts the binding capacity of CCL14-AS to MEP1A or LYVE1 mRNAs. A-B Prediction of binding capacity of CCL14-AS to MEPA1 A or LYVE1 B by LncRRIsearch database.
: Table S1. List of cancer-related lncRNAs (lncRNAs differently expressed in CRC tissues compared with).
: Table S2. List of metastasis-related lncRNAs (differently expressed lncRNA in metastatic samples compared with non-metastatic samples（log2FC<-0.5, adj.P.val<0.001)).
: Table S3. List of survival-related lncRNAs.
: Table S4. The coding potential of CCL14-AS predicted by the National Center for Biotechnology.
About this article
Cite this article
Li, M., Huang, C., Wu, Y. et al. Long non-coding RNA CCL14-AS suppresses invasiveness and lymph node metastasis of colorectal cancer cells by regulating MEP1A. Cancer Cell Int 23, 27 (2023). https://doi.org/10.1186/s12935-023-02866-1
- Long noncoding RNA
- Colorectal cancer