- Research article
- Open Access
Association of CASC18/miR-20a-3p/TGFB2 ceRNA axis with occult lymph node metastasis in tongue squamous cell carcinoma
Molecular Medicine volume 27, Article number: 85 (2021)
Tongue squamous cell carcinoma (TSCC) ranks as the most prevalent malignancy in the oral cavity. TSCC patients with occult lymph node metastasis (OLNM) are thought to be at risk of worse outcome. However, regulatory mechanisms underlying OLNM remain less investigated.
In the present study, CASC18/miR-20a-3p/TGFB2 axis was identified and evaluated by bioinformatic and qRT-PCR analyses. Effects of CASC18 knockdown on cell migration and invasion were determined by wound healing and transwell assays. Western blot, ELISA, RNA pulldown and luciferase reporter assays were performed for mechanism verification.
CASC18 was identified up-regulating in TSCC tumours, and especially in those from patients with OLNM. Importantly, we found higher CASC18 expression was positively correlated with the presence of OLNM and worse outcome of TSCC patients. Furthermore, we demonstrated that CASC18 knockdown repressed cell migration and invasion through inhibiting epithelial-mesenchymal transition, which could be partly rescued by miR-20a-3p inhibitor. Regarding the molecular mechanism, we further confirmed that CASC18 functioned as a ceRNA to sponge miR-20a-3p to enhanceTGFB2 expression and secretion.
In conclusion, we have reported a novel CASC18/miR-20a-3p/TGFB2 ceRNA axis in OLNM of TSCC. Our findings will contribute to a deeper understanding of the molecular mechanism of OLNM in TSCC, and facilitate the development of diagnostic methods for assisting treatment decision-making.
Tongue squamous cell carcinoma (TSCC) is the most prevalent malignancy in the oral cavity, with an increasing incidence in the past decades (Bray et al. 2018; Miller et al. 2019; Gulland and Gulland 2016). The presence of lymph node metastasis is one of the most accurate predictors of prognosis for TSCC patients. With advances in detection and treatment, TSCC outcomes have improved in recent years (Jansen et al. 2018). However, approximately 20–40 % of TSCC patients with a clinical N0 (cN0) neck are subsequently diagnosed with micro-metastasis in the lymph nodes after careful pathological examination (Ganly et al. 2012; Bilde et al. 2008). These patients harboring occult lymph node metastasis (OLNM) would miss the opportunities of therapeutic neck dissection unless their clinical N stage is accurately assessed prior to treatment. Thus, dissecting the molecular mechanisms involved in OLNM is of great clinical importance.
Recent high-throughput studies are contributing to an improved understanding of the molecular mechanism concerning various phenotypes, such as carcinogenesis (Sun et al. 2018), metastasis (Tian et al. 2017), and chemoresistance (Dai et al. 2019). These cancer-associated phenotypes are at least partly controlled by differentially expressed genes (DEGs), including long non-coding RNA (lncRNA), microRNA (miRNA), and protein-coding gene (PCG). Recently, these three types of molecules can be linked together according to the competing endogenous RNA (ceRNA) regulatory hypothesis (Salmena et al. 2011). Some lncRNAs can act as miRNA sponges to regulate the expression of PCG through competitively binding to shared miRNAs. Although increasing evidence has been well documented that ceRNA network plays pivotal roles in many cancers (Qi et al. 2015), even in TSCC, current knowledge for roles of ceRNA regulation in OLNM of TSCC is not enough.
In the present study, we first identified the ceRNA axes mediated by aberrantly expressed lncRNAs in TSCC with OLNM via mining the sub-dataset of TCGA-HNSC project. Functional enrichment analysis was performed to reveal the potential roles of this OLNM-related lncRNA. Through bioinformatic prediction and in vitro experiments, we further confirmed that CASC18 functions as a ceRNA to regulate cell migration and invasion via miR-20a-3p/TGFB2 axis. Our results provide novel insights into the ceRNA regulation mechanism and reveal potential roles of CASC18 in OLNM of TSCC.
Materials and methods
Tissues and cell line
34 FFPE samples of TSCC were collected with the consent of patients, and the experiments were approved by the ethics committee of Hunan Cancer Hospital (KYJJ-2020-222), Changsha, China. TSCC cell line CAL27 was purchased from Shanghai Genechem Co., Ltd, and cultured in DMEM with 10 % FBS (Gibco, Gaithersburg, MD,) as routine. For RNAi, siRNAs (si-CASC18#1: ACCTAGTCATACATCCTAA; si-CASC18#2: CCTTATCCCTGGATGGAAA) were transfected into CAL27 cells with Lipofectamine 2000 (Invitrogen, Carlsbad, CA) for 48 h.
Acquisition and DEG screening of RNA-seq data
The gene expression data and clinical information of TCGA-HNSC subset were downloaded from The Cancer Genome Atlas (TCGA; cancergenome.nih.gov), and data of TSCC patients were filtered. Subsequently, the DESeq2 (version 1.22) package (https://bioconductor.org/packages/DESeq2/) was used to screen for TSCC- and OLNM-related DEGs with the cut-off values of |log2FC|≥1 (where FC is fold change)and adjusted p < 0.05 (Love et al. 2014). The ggplot2 (version 3.3, cran.r-project.org/package = ggplot2) and ggpubr (version 0.2.3, https://rpkgs.datanovia.com/ggpubr) packages were used to display the results as volcano plots (Maag and Maag 2018; Ito et al. 2013). Then, the intersection of the DEGs was calculated, and the result was visualized as Venn diagrams using the VennDiagram package (https://CRAN.R-project.org/package=VennDiagram) in R (version 3.4.3; www.r-project.org) (Chen et al. 2011) . To identify the correlation between differentially expressed lncRNAs and PCGs, the Pearson correlation coefficient (PCC) value was subsequently calculated for each pair of lncRNA-PCG. Pairs with |PCC| > 0.45 (p < 0.05) and corrected p-value < 0.05 were considered as statistically significant correlation.
Functional enrichment analysis
The WEB-based GEne SeT AnaLysis Toolkit (WebGestalt, version 2019, http://www.webgestalt.org/) (Liao et al. 2019) was applied to map candidate PCGs, correlated with differentially expressed lncRNAs, onto their associated biological annotation, with Gene Ontology (GO; www.geneontology.org) (Gene Ontology Consortium 2017; Ashburner et al. 2000). All significantly enriched terms, with adjusted p < 0.05, were visualized by using the ggcorrplot (version 0.1. 3, CRAN.R-project.org/package = ggcorrplot) package.
Fluorescence in situ hybridization (FISH)
FISH was executed to determine the subcellular location of CASC18 in CAL27 cells. Fluorogenic probes targeting CASC18, 18 S and U6 were designed and synthesized by RiboBio (Guangzhou, China). Hybridization was carried out overnight according to the manufacturer’s instructions. All fluorescent images were captured with the same exposure setting and illumination.
For Formalin-Fixed and Paraffin-Embedded (FFPE) samples, the total RNA was extracted from at least four 10-µm FFPE scrolls using the RNeasy FFPE Kits (QIAGEN, Germantown, MD), and then reverse transcribed to cDNA using PrimeScriptTM RT-PCR Kit (Takara, Dalian, China) in accordance with the manuals. The expression level of CASC18 was detected, and samples were divided into low and high expression groups by the median value of this lncRNA. Subsequently, comparison of VIM (vimentin) and CDH1 (E-cadherin) expression between two groups were performed. For cells, nuclear and cytoplasmic RNA were isolated using the PARIS™ Kit (Ambion), and the total RNA was extracted by using Trizol reagent (Invitrogen, Carlsbad, CA). For miRNA quantitation, reverse transcriptions were performed using the PrimeScript RT Reagent Kit (Takara, Dalian, China) with specific stem-loop primers, designed and synthesized by RiboBio (Guangzhou, China). qRT-PCR was performed using SYBR® Premix DimerEraser™ (Takara, Dalian, China) in Roche LightCycler 480 II Real-Time PCR system (Roche, Basel, Switzerland). The threshold cycle value (Ct) of each product was determined and normalized against that of the internal control of β-actin (for lncRNA/PCG) or U6 (for miRNA), and the differences were compared by t-test using SPSS 23.0, with p < 0.05 considered as statistically significant. Primer pairs for gene expression quantitation were as follows: CASC18-F: 5′-TGTATAGTCTAGCCAAGTCC-3′, CASC18-R: 5′-ATTTCAGCCATCTTCAGTCCC-3′, β-actin-F: 5′-ACCCTGAAGTACCCCATCGAG-3′, β-actin-R: 5′-AGCACAGCCTGGATAGCAAC-3′, VIM-F: 5′-CGAGGAGAGCAGGATTTCTC-3′, VIM-R: 5′-GGTATCAACCAGAGGGAGTGA-3′, CDH1-F: 5′-TGCCCAGAAAATGAAAAAGG-3′, CDH1-R: 5′-TGCCCAGAAAATGAAAAAGG-3′.
The concentration of TGF-β2 in CAL27 culture supernatants was measured by using ELISA kit (Cloud-Clone corp., Wuhan, China), according to the manufacturer’s instructions.
Luciferase reporter assay
The fragments of CASC18-wt (wildtype) and CASC18-mut (miR-30a-3p binding site mutation) were ligated into pcDNA3.1 vector to form two plasmids. Vectors for Luciferase reporter assay were generated based on the pmirGLO Dual-Luciferase miRNA Target Expression Vector (Promega, Madison, WI). Then, pmirGLO-TGFB2 3′-UTR-wt (wildtype) and pmirGLO-TGFB2 3′-UTR-mut (miR-20a-3p binding site mutation) were constructed, and the plasmids was co-transfected with miR-20a-3p mimics, negative control, pcDNA3.1-CASC18-wt, and pcDNA3.1-CASC18-mut into CAL27 cells with Lipofectamine 2000 according to the guidelines, respectively. The relative luciferase activity was measured with the Dual-Luciferase Reporter Assay System (Promega, Madison, WI).
The biotinylated DNA probe complementary to CASC18 and was synthesized (GenePharm, Shanghai, China), with a random probe as the control. M-280 streptavidin magnetic beads (Sigma-Aldrich, St. Louis, MO) were coated by incubating with the probes, respectively. The cell lysates of CAL27 transfected with si-CASC18#2 or si-NC were prepared after 48 h. Lysates were then incubated with the probe-coated beads at 4 °C overnight, and molecules interacted with CASC18 were captured after washing. The bound RNAs were subsequently purified using TRIzol and the miR-20-3p abundance was finally measured by qRT-PCR.
Wound healing assay
As cells at 80 % confluence in a six-well culture dish, the monolayer was scraped with a sterile 200 µl pipette tip. The cells were cultured in a culture medium with 3 % FBS and 1 % penicillin/streptomycin (TransGen, Beijing, China). Images were captured at 0, 24, and 48 h after the scratch was made. Image-Pro Plus v6.0 image analysis software was used to analyse the cell wound healing rate as described in our previous study (Dong et al. 2019).
Transwell invasion assay
The transwell invasion assay was performed using the transwell (Corning, NY, USA) and matrigel (BD Biosciences, San Jose, CA) as described previously (Justus et al. 2014), with slight modifications. Cells were added into the upper chamber. After 36 h of incubation at 37 °C in 5 % CO2, cells invading through the matrigel and membrane to the lower surface were fixed and stained with Crystal violet for 10 min. The number of invaded cells were counted in 6 randomly selected visual fields under the microscope.
The harvested cells were lysed in RIPE buffer at 4 ℃ for 30 min and centrifuged at 15,000×g for 15 min for the protein sediment. Equal amount of proteins was loaded to SDS-polyacrylamide gel electrophoresis (SDS–PAGE) for separation, and then transferred to PVDF membranes (Millipore, Billerica, MA) according to the standard protocol. Subsequently, the membrane was incubated with primary antibodies (ProteinTech Group Inc., Chicago, IL) against TGF-β2 (1:2000 dilution), E-cadherin (1:1,000 dilution), N-cadherin (1:1000 dilution), Vimentin (1:1500 dilution) and β-actin (FN; 1:2000 dilution) overnight at 4 °C. After incubating with the corresponding secondary antibody for 1 h at room temperature, the immunoblotting signals were measure with enhanced chemiluminescence (ECL) kit (Pierce, Rockford, IL).
Data were presented as mean ± SD from at least three separate experiments, with a p-value < 0.05 considered statistically significant. The Student t test, Wilcox test, log-rank test, or Fisher exact test was used for comparisons between groups. The Kaplan–Meier method was used to estimate overall survival.
Identification of DEGs associated with OLNM in TSCC
A total of 53 cN0M0 TSCC patients were selected from the TCGA-HNSC cohort through filtration of TNM stage and anatomic sites. To understand expression differences between TSCC patients with and without OLNM, we carried out differential gene expression analysis and identified 28 differentially expressed lncRNAs (Fig. 1A). Among them, 6 genes were up-regulated and 22 genes were down-regulated (Additional file 1: Table S1). Then, expression data of 13 pairs of TSCC tissues and corresponding adjacent non-cancerous tongue tissues were obtained from TCGA-HNSC cohort. Comparing the TSCC tissues with the normal tongue tissues, a total of 1224 DEGs were identified, comprising 746 up-regulated lncRNAs and 478 down-regulated lncRNAs (Fig. 1B, Additional file 1: Table S2). Subsequently, when the DEGs were investigated for overlap, a total of 9 consistently aberrant genes were identified (Fig. 1C).
Exploration of potential ceRNA axes by bioinformatics analysis
The correlation between each lncRNA and PCG expression was evaluated by the Pearson correlation coefficient (PCC) analysis in the TCGA-TSCC subset. A total of 1758 PCGs were found correlated with at least one of 9 common lncRNAs (PCC > 0.45, p < 0.05). Subsequently, these 1758 PCGs were subjected to GO enrichment analysis. As shown in Fig. 2A, they were enriched significantly into 15 diverse GO terms, being categorized into three functional groups: Biological process (BP), cellular component (CC) and molecular function (MF). Among these terms, “cell-cell adhesion via plasma-membrane adhesion molecules” was the term related to cancer metastasis, and abnormal expression of genes in “epidermis development” may partly contribute to tumorigenesis of TSCC. Then, overlapping genes, TGFB2, DSG1, DSC1 and DSC2, were identified between mentioned two terms by the intersecting analysis (Fig. 2B). The correlation of 9 lncRNAs to 4 PCGs was visualized in Fig. 2C. With a multi-step filtration strategy (Fig. 3), we then investigated the potential regulatory axes between each lncRNAs and PCGs. A total of 539 miRNA-lncRNA pairs were predicted by lncBase computational methods. After removal of miRNAs (1) expressed in less than 50 % TCGA-TSCC samples, (2) with interaction score < 0.85, (3) with < 5 binding sites in the target lncRNA, 12 candidate miRNAs were selected. By using miRWalk computational methods, 1525 miRNA-mRNA pairs were obtained. A list of 220 miRNAs was collected after deleting miRNAs expressed in less than 50 % TCGA-TSCC samples. Overlapping the lncRNA-related miRNA and mRNA-related miRNA, we selected 4 candidate miRNAs and formed 9 triple gene axes. Interestingly, only the CASC18/miR-20a-3p/TGFB2 axis met lncRNA-mRNA positive correlation, and thus was selected for further validation.
Expression and characterization of CASC18 in TSCC
In TCGA-TSCC cohort, the significantly lower level of CASC18 expression was detected in normal tongue tissues when compared with TSCC (p = 0.007) but also corresponding tumor tissues (p = 0.040), respectively (Fig. 4A, B). Moreover, Fig. 4C showed that CASC18 obviously enhanced in tumor samples of TSCC patients with OLNM (p = 0.001). Results of PCC analysis (Fig. 4D) showed markedly positive correlation between CASC18 and TGFB2 (PCC = 0.487, p < 0.01). What is more important, increased CASC18 expression is associated with worse prognosis of patients in TCGA-TSCC cohort (Fig. 4E). To further investigate the subcellular localization of this lncRNA, fluorescence in situ hybridization (FISH) and quantitative real-time PCR (qRT-PCR) assays were performed. The results demonstrated that CASC18 was predominantly localized in the cytoplasm (Fig. 4F, G), implying that CASC18 can act as a sponge of miRNA and function as a ceRNA in cytoplasm. Then, we further examined CASC18 expression in an independent cohort of 34 TSCC FFPE samples, and determined the correlation between their CASC18 expression and clinicopathologic characteristics. As shown in Table 1, high CASC18 expression was positively associated with TSCC tumors characterized by the Pathologic N stage (p = 0.006) as well as OLNM (p = 0.006).
Knockdown of CASC18 inhibited cell migration and invasion in vitro
To investigate the biological function of CASC18, we interfered its expression by transient transfecting of siRNA in CAL27 cells. Obviously, si-CASC18#2 displayed the most significant suppression on silencing CASC18 (Fig. 5A). Wound healing assay showed that CASC18 knockdown repressed cell migration in CAL27 cells, and si-CASC18#2 displayed a greater inhibitory effect (Fig. 5B). Likewise, the transwell assay also indicated the suppressive effect of CASC18 downregulation on cell invasion (Fig. 5C). However, MTT assay did not show a significant effect of CASC18 expression on cell proliferation (data not shown). Considering the association of epithelial-mesenchymal transition (EMT) with cell mobility and invasion, as well as the key role of TGF-β pathway in EMT, we further detected the expression level of TGF-β2 and EMT-related markers by western blot. Our result revealed that expression of TGF-β2, as well as N-cadherin and Vimentin, which are mesenchymal markers, was decreased; whereas expression of E-cadherin, an epithelial marker, was alternatively increased after CASC18 RNAi (Fig. 5D). Additionally, results of qRT-PCR indicated that VIM was highly expressed in FFPE samples carrying high CASC18 expression, while CDH1 was highly expressed in those carrying low CASC18 expression (Fig. 5E). Besides, we found that CASC18 knockdown significantly decreased TGF-β2 secretion in CAL27 culture supernatants (Fig. 5F). These results suggested the correlation among CASC18, TGFB2 and EMT in promotion of CAL27 cell migration and invasion.
CASC18 acts as a miR-20a-3p sponge to prompt CAL27 cell migration and invasion via regulating TGFB2 expression
To compare levels of miR-20a-3p in CAL27 cells after CASC18 knockdown, we transfected CAL27 cell with siRNAs against CASC18 and si-NC. The expression of miR-20a-3p unchanged with the decrease of CASC18 (Fig. 6 A), indicating that CASC18 did not regulate miR-20a-3p expression. We then constructed luciferase plasmids containing a TGFB2 wildtype 3′UTR (TGFB2 3′UTR-wt) or a mutated TGFB2 3′UTR (TGFB2 3′UTR-mut), whose miR-20a-3p binding site mutated. Luciferase assay showed that miR-20a-3p inhibitor significantly enhanced luciferase activity of TGFB2 3′UTR-wt, but it presented less effect on the TGFB2 3′UTR-mut (Fig. 6B). Moreover, CASC18 knockdown could also decrease the luciferase activity of TGFB2 3′UTR-wt, and this change can be partially eliminated by treating with the miR-20a-3p inhibitor. What is more important, miR-20a-3p inhibitor and CASC18 knockdown had no influence on the luciferase activity of TGFB2 3′UTR-mut, indicating regulation in a miR-20a-3p binding-site dependent manner. To next determine the binding property between miR-20a-3p and CASC18, RNA-pull down assay was performed. The biotinylated probe was firstly used to pulldown CASC18 in CAL27. (Fig. 6C, left panel). Subsequently, the enrichment of miR-20a-3p, which was captured by CASC18, was measured in CAL27. However, silencing of CASC18 significantly decreased the abundance of miR-20a-3p (Fig. 6C, right panel), indicating the binding interaction between miR-20a-3p and CASC18. The above results suggest that CASC18 elevates TGFB2 expression and secretion via sponging miR-20a-3p in TSCC cells. To further validate the role of CASC18/mir-20a-3p/TGFB2 axis in CAL27 cell migration and invasion, we performed the rescue experiments. As expected, exogenous miR-20a-3p inhibition significantly deprived the suppression role of cell mobility by CASC18 downregulation in CAL27, when compared with CASC18 RNAi alone (Fig. 7A, B). Moreover, results of western blot and Enzyme-linked immunosorbent assay (ELISA) revealed that expression and secretion of TGF-β2 were alternatively increased after miR-20a-3p inhibition in CASC18 knockdown CAL27 (Fig. 7C, D). Taken together, our results suggest that CASC18 stimulates EMT via relieving miR-20a-3p-mediated TGFB2 repression, and thus plays an ultimate role in regulating TSCC cell migration and invasion (Fig. 8).
For TSCC patients harbouring OLNM, prophylactic neck dissection was not routinely performed, attributing to the misclassification of a cN0 neck. Consequently, most of these patients suffered from the disease recurrence and metastases with a worse outcome. However, the molecular mechanisms involved are not clearly defined. In the present study, we found CASC18 expression was associated with the status of OLNM in TSCC cohorts, and acted as a ceRNA to modulate cell migration and invasion via miR-20a-3p/TGFB2 axis.
In the last decade, considerable studies have identified that lncRNAs played critical roles in carcinogenesis and metastasis (Lv et al. 2021; Yao et al. 2019; Poursheikhani et al. 2021). To data, the ceRNA hypothesis sparked new areas of lncRNA research on gene expression regulation, but some controversy still surrounded it (Tay et al. 2011; Smillie et al. 2018). The prediction of ceRNA axes, which is wholly dependent on the web tool for identifying shared MRE identification, is an overwhelming limitation of studies assessing ceRNA regulation(Thomson et al. 2016). This strategy is common in bioinformatics studies without experimental approach (Liu et al. 2020; Tang et al. 2020), and high false positive rate of their results is inevitable. Additionally, the cytoplasm abundance of miRNA is critical for its efficacy, and some lncRNA is expressed much higher for miRNAs with lower expression (Di Leva et al. 2013). Accordingly, we screened out the OLNM-related ceRNA axes using a multi-step filtration that takes into account the existence and number of MREs in lncRNA, positive correlation between lncRNA and associated-PCG, expression level of miRNA in TSCC samples, and computational interaction score, to discriminate signal from noise. Together with the intersection of PCGs enriched in both the terms of “cell-cell adhesion via plasma-membrane adhesion molecules” and “epidermis development”, only the CASC18/miR-20a-3p/TGFB2 axis, who fulfilled all the condition, was selected for further validation.
Although CASC18 named as Cancer Susceptibility 18, this gene was previously reported involved in neural cell differentiation and oocyte growth (Bouckenheimer et al. 2018; Mehravar et al. 2017), instead of cancer. To verify the function of the lncRNA in OLNM of TSCC, we first investigated the expression pattern and prognosis relationship of CASC18 in TCGA-TSCC subset. We found CASC18 is a prognosis factor up-regulated in tumour samples and especially in those of TSCC patients with OLNM, partly supporting the roles of CASC18 in lymph node metastases. Several models, including orthotopic tumour model (Eun et al. 2020), subcutaneous tumour xenograft model (Knudsen et al. 2021), tumour cell tail vein injection for metastasis model (Zhao et al. 2021), were extensively applied for cancer research. Regrettably, appropriate animal model for OLNM is still unavailable, and thus the association of CASC18 expression with OLNM was not verified in animal model in the present study. However, we found high CASC18 expression was positively correlated with the presence of OLNM in an independent TSCC cohort, which provides a clinical evidence for the association between CASC18 expression and OLNM in TSCC.
According to the predicted ceRNA axis, it has been easy to come away with the idea that CASC18 was involved in the regulation of EMT, due to the representative role of TGF-β signals in this process (Hao et al. 2019). EMT defines dynamic changes in cellular organization from epithelial to mesenchymal phenotypes, which is commonly considered to facilitate cell migration and invasion (Yang et al. 2020), particularly in cancers (Brabletz et al. 2018). Our in vitro experiments indicated that CASC18 influenced TGFB2 expression and secretion, and affected CAL27 cell migration and invasion. Consistent with previous studies (Thieme 2020; Jiang et al. 2021), EMT was inhibited with the attenuation of TGFB2 expression in CASC8-knockdown CAL27 cells. These data suggested that CASC18 regulated cell mobility via TGFB2-mediated EMT in TSCC cells. In order to further detect the underlying mechanism, we then investigated the subcellular localization of CASC18. As expected, this lncRNA was mainly distributed in the cytoplasm of CAL27 cells, in accordance with the basic requirement of ceRNA regulation (Salmena et al. 2011). To note, target-directed miRNA decay (TDMD) mechanism is another regulation mode analogous to ceRNA in some way, emphasizing the degradation of miRNA mediated by its binding target and ubiquitin ligase ZSWIM8 (Han et al. 2020; Shi et al. 2020; Sheu-Gruttadauria et al. 2019). To exclude the possibility of TDMD, we investigated miR-20a-3p expression after CASC18 knockdown, and our results indicated the stable expression of the miRNA in CAL27. Accordingly, we further verified the promising ceRNA regulation of CASC18/miR-20a-3p/TGFB2 axis. We confirmed that TGFB2 was a direct target of miR-20a-3p, and could be decreased by CASC18 knockdown in a miR-20a-3p-binding dependent manner by using luciferase reporter assays. Results of CASC18 pulldown assays further showed the direct binding of CASC18 and miR-20a-3p. Furthermore, we found that CASC18 knockdown suppressed the TGFB2 protein levels and secretion, but these effects were partly reversed by miR-20a-3p inhibitor treatment in CAL27 cells. Similarly, miR-20a-3p inhibitor treatment partly rescued the suppression effect of CASC18 knockdown on CAL27 cell migration and invasion. Taken together, our findings indicate that CASC18 functions as miR-20a-3p sponge to regulate cell migration and invasion via stimulating TGFB2-mediated EMT in CAL27.
In conclusion, we have uncovered a novel CASC18/miR-20a-3p/TGFB2 axis and elaborated its involvements in OLNM of TSCC patients. Our study provided pioneering insights and fundamental clues for future thorough understanding of the molecular mechanism of TSCC OLNM, which might facilitate the development of improved diagnostic methods for treatment decision-making.
Availability of data and materials
The datasets analysed during the current study are available in TCGA (cancergenome.nih.gov) repository. TCGA belongs to the public database and the patients involved in the database have obtained ethical approval. TCGA allows users to download relevant data for free for research and publish relevant articles.
Ashburner M, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.
Bilde A, et al. Need for intensive histopathologic analysis to determine lymph node metastases when using sentinel node biopsy in oral cancer. Laryngoscope. 2008;118:408–14.
Bouckenheimer J, et al. Differential long non-coding RNA expression profiles in human oocytes and cumulus cells. Sci Rep. 2018;8:2202.
Brabletz T, Kalluri R, Angela Nieto M, Weinberg RA. EMT in cancer. Nat Rev Cancer. 2018;18:128-+.
Bray F, et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424.
Chen H, Boutros PC. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinform. 2011;12:35.
Dai P, et al. Screening candidate microRNA-mRNA network for predicting the response to chemoresistance in osteosarcoma by bioinformatics analysis. J Cell Biochem. 2019;120:16798–810.
Di Leva G, et al. Estrogen mediated-activation of miR-191/425 cluster modulates tumorigenicity of breast cancer cells depending on estrogen receptor status. PLoS Genet. 2013;9:e1003311.
Dong S, et al. HOXD-AS1 promotes the epithelial to mesenchymal transition of ovarian cancer cells by regulating miR-186-5p and PIK3R3. J Exp Clin Cancer Res. 2019;38:110.
Eun Y-G, Yoon YJ, Won KY, Lee YC. Circulating tumor DNA in saliva in an orthotopic head and neck cancer mouse model. Anticancer Res. 2020;40:191–9.
Ganly I, Patel S, Shah J. Early stage squamous cell cancer of the oral tongue–clinicopathologic features affecting outcome. Cancer. 2012;118:101–11.
Gene Ontology Consortium. Expansion of the gene ontology knowledgebase and resources. Nucleic Acids Res 2017;45: D331-D338.
Gulland A. Oral cancer rates rise by two thirds. BMJ. 2016;355:i6369.
Han J, et al. A ubiquitin ligase mediates target-directed microRNA decay independently of tailing and trimming. Science 2020;370.
Hao Y, Baker D, Ten Dijke P. TGF-beta-mediated epithelial-mesenchymal transition and cancer metastasis. Int J Mol Sci. 2019;20:2767.
Ito K, Murphy D. Application of ggplot2 to pharmacometric graphics. CPT: Pharmacomet Syst Pharmacol. 2013;2:e79.
Jansen L, et al. Differences in incidence and survival of oral cavity and pharyngeal cancers between Germany and the United States depend on the HPV-association of the cancer site. Oral Oncol. 2018;76:8–15.
Jiang XM, et al. TGFbeta2-mediated epithelial-mesenchymal transition and NF-kappaB pathway activation contribute to osimertinib resistance. Acta Pharmacol Sin. 2021;42:451–9.
Justus CR, Leffler N, Ruiz-Echevarria M, Yang LV. In vitro cell migration and invasion assays. J Visualized Exp: JoVE 2014;51046.
Knudsen ES, et al. Targeting dual signalling pathways in concert with immune checkpoints for the treatment of pancreatic cancer. Gut. 2021;70:127–38.
Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B. WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res. 2019;47:W199-205.
Liu H, et al. The ceRNA network has potential prognostic value in clear cell renal cell carcinoma: a study based on TCGA database. Biomed Res Int. 2020;4830847.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Lv E, Sheng J, Yu C, Rao D, Huang W. LncRNA influence sequential steps of hepatocellular carcinoma metastasis. Biomed Pharmacother. 2021;136:111224.
Maag JLV. gganatogram: an R package for modular visualisation of anatograms and tissues based on ggplot2. F1000Research . 2018;7:1576.
Mehravar M, et al. Introduction of novel splice variants for CASC18 gene and its relation to the neural differentiation. Gene. 2017;603:27–33.
Miller KD, et al. Cancer treatment and survivorship statistics, 2019. CA Cancer J Clin. 2019;69:363–85.
Poursheikhani A, Abbaszadegan MR, Kerachian MA. Mechanisms of long non-coding RNA function in colorectal cancer tumorigenesis. Asia Pac J Clin Oncol. 2021;17:7–23.
Qi X, et al. ceRNA in cancer: possible functions and clinical implications. J Med Genet. 2015;52:710–8.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta stone of a hidden RNA language? Cell. 2011;146:353–8.
Sheu-Gruttadauria J, et al. Structural basis for target-directed MicroRNA degradation. Mol Cell. 2019;75:1243-55 e1247.
Shi CY, et al. The ZSWIM8 ubiquitin ligase mediates target-directed microRNA degradation. Science 2020;370.
Smillie CL, Sirey T, Ponting CP. Complexities of post-transcriptional regulation and the modeling of ceRNA crosstalk. Crit Rev Biochem Mol Biol. 2018;53:231–45.
Sun W, Qiu Z, Huang W, Cao M. Gene expression profiles and proteinprotein interaction networks during tongue carcinogenesis in the tumor microenvironment. Mol Med Rep. 2018;17:165–71.
Tang H, Wang Z, Shao Q, Wang Y, Yang Q. Comprehensive analysis of competing endogenous RNA (ceRNA) network based on RNAs differentially expressed in lung adenocarcinoma using The Cancer Genome Atlas (TCGA) database. Med Sci Monit. 2020;26:e922676.
Tay Y, et al. Coding-independent regulation of the tumor suppressor PTEN by competing endogenous mRNAs. Cell. 2011;147:344–57.
Thieme R, et al. TGFB1 and TGFB2 Mediated Epithelial-Mesenchymal Transition in Esophageal Adenocarcinoma Cells. Oncology Research Treatment. 2020;43:68–8.
Thomson DW, Dinger ME. Endogenous microRNA sponges: evidence and controversy. Nat Rev Genet. 2016;17:272–83.
Tian F, Zhao J, Fan X, Kang Z. Weighted gene co-expression network analysis in identification of metastasis-related genes of lung squamous cell carcinoma based on the Cancer Genome Atlas database. J Thorac Dis. 2017;9:42–53.
Yang J, et al. Guidelines and definitions for research on epithelial-mesenchymal transition. Nat Rev Mol Cell Biol. 2020;21:341–52.
Yao R-W, Wang Y, Chen L-L. Cellular functions of long noncoding RNAs. Nat Cell Biol. 2019;21:542–51.
Zhao B, et al. (2021) A novel oncotherapy strategy, direct thrombin inhibitors suppress progression, dissemination and spontaneous metastasis in non-small cell lung cancer. Br J Pharmacol.
We acknowledge TCGA database for providing their platforms and contributors for uploading their meaningful datasets.
The present study was supported in part by the National Natural Science Foundation of China (Grant no. 81703005), the Natural Science Foundation of Hunan Province (Grant nos. 2017JJ3195, 2018JJ3317 and 2021JJ30427), the Research Project of Health and Family Planning Commission of Hunan Province (Grant no. C2019080 and Grant no.202102082037) and the Key Clinical Specialty Construction Project (Head&Neck Surgery) of Provincial Health Commission of Hunan province.
Ethics approval and consent to participate
The present study was approved by the ethics committee of Hunan Cancer Hospital (KYJJ-2020-222), Changsha, China.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhou, B., Zhou, Y., Liu, Y. et al. Association of CASC18/miR-20a-3p/TGFB2 ceRNA axis with occult lymph node metastasis in tongue squamous cell carcinoma. Mol Med 27, 85 (2021). https://doi.org/10.1186/s10020-021-00345-9
- Occult lymph node metastasis
- Tongue squamous cell carcinoma