Hsa_circRNA_102002 facilitates metastasis of papillary thyroid cancer through regulating miR-488-3p/HAS2 axis

As important modulators in various physiological processes, circular RNAs (circRNAs) have been increasingly demonstrated in tumors, including papillary thyroid cancer (PTC). Hsa_circRNA_102002 (circ_102002) is a circRNA derived from alternative splicing of ubiquitin-specific peptidase 22 (USP22) transcript, the role of which needs further investigation. Our results suggested the upregulation of circ_102002 in PTC tissues and cells, and its promoting effects on epithelial–mesenchymal transition (EMT) and cell migration. Mechanism studies showed that circ_102002 could sponge microRNA-488-3p (miR-488-3p) and downregulate its expression. The target relationship between miR-488-3p and hyaluronic acid synthetase 2 (HAS2) in PTC was systematically studied. In addition, our results showed that HAS2 overexpression could restore the inhibited cell EMT and migration. Moreover, the inhibitory effect of downregulation of circ_102002 on PTC growth was evaluated in a mouse xenograft model, which involved miR-488-3p and HAS2 regulation. These findings about the signal axis of circ_102002/miR-488-3p/HAS2 may further elucidate the PTC pathogenesis and improve clinical treatment.


Introduction
Papillary thyroid cancer (PTC) is the most common endocrine malignant tumor, representing~90% of all thyroid cancers [1]. Although the detection and treatment methods have made some progress and improved the prognosis, metastases commonly involve cervical lymph nodes and lung [2]. Lung metastasis occurs in 7-30% in children and adolescents (≤20 years), only~4% in adult population [3,4]. These PTC patients still suffered from poor outcomes. Therefore, exploration of the molecular mechanism related to the pathogenesis of PTC metastasis may offer alternative intervention targets and improve the therapeutic strategies.
CircRNAs have been identified as endogenous noncoding RNAs, which are widely found in eukaryotes and possess tissue specificity [5][6][7]. Characterized as a covalently closed continuous loop, circRNA is resistant to cleavage of RNA exoribonuclease, making it a regulator of microRNA (miRNA) binding, protein interaction and regulatory splicing [8,9]. The characteristics and functions of circRNAs are increasingly demonstrated in PTC. For example, circ_0067934 was highly expressed in PTC tissues, while circ_0067934 knockdown inhibited epithelial-mesenchymal transition (EMT) and PI3K/AKT signaling pathways [10]. In addition, circRNA_102171 promotes PTC progression by modulating the activation of CTNNBIP1-dependent β-catenin pathway [11]. Hsa_circRNA_102002, derived from alternative splicing of ubiquitin-specific peptidase 22 (USP22) transcript, is profiled to be upregulated in PTC tissues by microarray [12]. Although USP22 has been reported to be related to glioma progression [13], colorectal cancer drug resistance [14], and PTC metastasis [15], the role of hsa_circRNA_102002 in PTC has not been fully recognized. Serving as epigenetic regulators, circRNAs usually functions by sponging miRNAs, which further regulate the expression of complementary messenger RNAs and impart cancer development. Recent studies indicated that miR-488-3p functioned as a tumor suppressor in esophageal squamous cell carcinoma by targeting zinc finger and BTB domain containing 2 (ZBTB2) and activating p53 pathway [16]. Moreover, miR-488-3p could regulate protein kinase, DNA-activated, catalytic subunit (PRKDC) and sensitize malignant melanoma cells to cisplatin [17]. However, the function of miR-488-3p and potential target genes in PTC are still unknown.
As a member of hyaluronic acid synthetase (HAS), HAS2 is a transmembrane glycosyltransferase and a major ECM proteoglycan, which plays an important role in the high expression of macromolecular hyaluronic acid during the malignant transformation [18]. In the tumor microenvironment, macromolecular hyaluronic acid provides a hydrated matrix to create gaps in the extracellular matrix that facilitate tumor cell migration [19]. Recent studies have found that hyaluronic acid is involved in a number of cancer processes, including multidrug resistance and angiogenesis. HAS2 was found to promote ZEB1-mediated EMT and metastasis in breast cancer [19]. In addition, HAS2 could activate the PI3K/Akt pathway and increase the invasion capacity of breast cancer cells by inhibiting timp-1 and promoting FAK phosphorylation [20]. Moreover, the function of HAS2 in promoting tumor progression have been demonstrated in pancreatic cancer [21], colorectal cancer [22] etc. However, the regulatory effects and mechanisms of HAS2 in PTC remains to be elucidated.
In this study, we determined the circ_102002 expression in PTC tissues and explored the effect of circ_102002 on EMT and cell migration. Potential target miRNAs by circ_102002 were predicted by starBase and Circular RNA Interactome. Based on target prediction analysis, we found that miR-488-3p might target the 3ʹ-UTR of HAS2. Therefore, it can be speculated that circ_102002 might promote PTC EMT and metastasis by sponging miR-488-3p and upregulating HAS2 expression. The regulatory manner of circ_102002/miR-488-3p/HAS2 axis in PTC was systematically studied, which may help to develop new candidates for PTC diagnosis and therapy.  RNA extraction, microarray analysis, and quantitative real-time PCR (qRT-PCR) Total RNAs were extracted from PTC cells and tissues using TRIzol reagent (Thermo Fisher Scientific, Rockford, IL, USA). The RNA quality was evaluated by electrophoresis in agarose and observed by Agilent 2100 Bioanalyzer (Agilent Technologies, Chandler, USA). Microarray of circRNAs was performed in BGI.Tech (Beijing, China) and analyzed according to the Arraystar. Nucleus and cytoplasm were separated using an extraction Kit (Inventbiotech, Beijing, China). Complementary DNA was synthesized from total RNA by using PrimeScript RT Reagent Kit (Takara, Shiga, Japan). qRT-PCR assays for genes were conducted by SYBR Premix Ex Taq (Takara). SYBR PrimeScript miRNA RT-PCR Kit (Takara) was applied to quantity miRNA expressions. GAPDH (glyceraldehyde-3-phosphate  20 15 dehydrogenase) is an enzyme in the classic glycolysis reaction, accounting for 10-20% of total proteins. The GAPDH gene has a highly conserved sequence and its expression in the same cell or tissue is generally constant. To normalize circ_102002 and HAS2 mRNA levels, GAPDH was selected as endogenous control. U6 snRNA was used to normalize miRNA expressions. After the Ct values of GAPDH or U6 in each sample were obtained by qRT-PCR, the stability of the reference genes was analyzed by internal reference screening software (geNorm, NormFinder and Bestkeeper). The amount of RNA to relative to the reference gene was calculated using the equation 2 -ΔΔCT . The corresponding primer sequences were listed in Table S1.

Overall survival analysis
After qRT-PCR analysis, the average circ_102002 level was calculated and used as a criterion to evaluate the expression of circ_102002. If circ_102002 level exceeds the average level, high expression level was identified. When the circ_102002 level was under the average value, low expression level was classified. According to patients' follow-up data, overall survival curves were analyzed by the Kaplan-Meier method using log-rank test.

In-situ hybridization (ISH)
An ISH kit (Thermo Fisher Scientific) was used to perform fluorescence detection of circ_102002 according to the manufacturer's instructions. Tissues were first fixed in paraformaldehyde and then embedded using paraffin. After that, sections (5 μm) were cut. For tissue ISH, circ_102002 probe was labeled with peroxidase (Thermo Fisher Scientific). For cell ISH, cells were fixed in paraformaldehyde, blocked in pre-hybridization buffer, incubated with a circ_102002 probe conjugated by FITC (Exiqon, Vedbaek, Denmark) in hybridization buffer, and finally stained by DAPI (Sigma-Aldrich). Fluorescence images were photographed by a confocal microscope (Nikon, Tokyo, Japan).
During the experiments of ISH, tissue sections hybridizated with mismatch control probe instead of circ_102002 probe were used as negative controls.

Immunofluorescence (IF) and immunohistochemistry (IHC) assays
Tissue sections were first cut. Since methylene bridge was generated in the process of formaldehyde fixation, the proteins were cross-linked to screen the antigen sites. The sections were immersed in citrate repair solution, and then heated in the microwave until boiling. The subboiling temperature (95-98°C) was maintained for 10 min, and the antigens were reexposed by heat to restore the original

Western blotting
Protein was extracted from PTC cells by RIPA lysis buffer (Sigma-Aldrich). Protein quantification was performed by a BCA protein assay kit (Thermo Fisher Scientific). After that, 50 µg total proteins from each sample) were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), and transferred onto polyvinylidene fluoride membranes (Millipore, Bedford, MA, USA), Then, immunoblot analyses was conducted by incubation with primary antibodies (overnight at 4°C), including Vimentin, Slug, Twist, MMP2, MMP9, HAS2, p-FAK, FAK, p-AKT, and AKT. After that, membranes were incubated with a horseradish peroxidase-conjugated secondary antibodyand and visualized using the enhanced chemiluminescence kit (Santa Cruz, Dallas, TX, USA). GAPDH was used as internal control and protein bands were quantified using Quantity One software (Bio-Rad, Berkeley, CA, USA). For antibody information, see Supplementary Table S2.

Wound scratch assay
Linear scratch wounds were conducted on the cell monolayer. After 24 h, an inverted microscope (Nikon) was subjected to take photos of cell migration. Then, the distances of the scrape were determined using an Image J software (http://rsb.info.nih.gov/ij/).

Transwell assay
For cell migration assay, SNU-790 or TPC-1 cells (1 × 10 5 ) were cultured in the top chamber of an insert (Corning Costar Co., Cambridge, MA, USA). For cell invasion assay, the insert was pre-coated by Matrigel (working concentration: 2 mg/ml, BD Biosciences, San Jose, CA). In the top chamber, cell culture was performed using serum-free medium. Cell culture medium with FBS (10%) was applied to induce cells migrating to the lower chamber of the insert. Twelve hours later, the migrated and invaded cells were fixed. The matrigel and cells in the upper chamber were swabbed away by a cotton swab. Moreover, cells underside of the membrane were then stained with 0.1% crystal violet, and imaged by a microscope (Nikon). The number of cells in ten different fields was counted and compared with the control group.

RNA binding protein immunoprecipitation (RIP) assay
RIP assay was performed as previously reported [24]. Ago2 antibodies or IgG were used to incubate the cell lysates of SNU-790 or TPC-1 cells. RNA was extracted from the pull down complexes, and relative levels of miR-488-3p, miR-625-5p, and miR-1197 were then determined using qRT-PCR.

Plasmid construction
The sequence of circ_102002 was amplified using a human genomic DNA and inserted into the p-MIR-reporter plasmid (Promega, Madison, WI, USA). The mutant luciferase plasmid (binding site: UGGCCUUUCA substituted by AUCGG AAAGU) was also used as a control. According to the methods above, the p-MIR-reporter plasmid containing the 3ʹ-UTR of HAS2 was obtained. In addition, the corresponding mutant luciferase plasmid (binding site: UCAAGG AAAAGUUCUUUCA substituted by UCAAGAAAA GUUGAAAGU) was constructed as a negative control. Primers sequences were listed in Table S3.

Statistical analysis
Experiments were performed in triplicate and repeated for at least three times independently. Samples and animals were randomly allocated to experimental groups and processed, and was blind with the investigator. Data were checked and they obeyed a normal distribution. Results were presented as the mean ± standard error of the mean (SEM). Student's t-test was used to analyze the statistical differences between two groups, and one-way analysis of variance was used for multiple groups. The correlations between circ_102002 and miR-488-3p, circ_102002 and HAS2 mRNA, miR-488-3p and HAS2 mRNA were, respectively, determined by Pearson correlation analysis. p < 0.01 was classified as significant. */# p < 0.05; **/# # p < 0.01.

Results
Hsa_circ_102002 was upregulated in papillary thyroid carcinoma tissues and cells To explore the role of circRNAs in papillary thyroid cancer, circRNA microarray was used to determine differently expressed circRNAs in PTC tissues (tumor, n = 3) and normal adjacent tissues (non-tumor, n = 3). Among them, hsa_circ_102002 was found to be the most significant differentially expressed circRNA in PTC samples compared to normal samples (Fig. 1a). After that, 50 pairs of PTC tissues and normal adjacent tissues were used to analyze hsa_-circ_102002 expression by qRT-PCR. The results showed that hsa_circ_102002 was remarkably upregulated in PTC tissues (Fig. 1b). As shown in Fig. 1c, PTC patients with high hsa_circ_102002 expression had a relatively lower survival rate. Patients in high hsa_circ_102002 expression group had a higher probability of lymph node metastasis and distant metastases, and their 5-year overall survival rate was only 50%. While, patients in low hsa_circ_102002 expression group were not prone to metastases, and their 5year survival rate was close to 90%. In addition, the results of ISH also showed that hsa_circ_102002 was highly expressed in PTC tissues (Fig. 1d). Similarly, the levels of hsa_circ_102002 in PTC cells, such as TPC-1, KTC-1, CAL62 and SNU-790, were higher than that in normal thyroid cells (HTori-3). Since SNU-790 or TPC-1 cells expressed a relatively lower or higher hsa_circ_102002 expression, they were used in the following study (Fig. 1e). RNase R is a ribonuclease and can degrade linear RNA, but is inefficient for circular RNA. As shown in Fig. 1f, RNase R treatment could degrade the linear transcript of USP22, but did not work for hsa_circ_102002. Subcellular localization of hsa_circ_102002 in SNU-790 or TPC-1 cells indicated that hsa_circ_102002 was mainly expressed in cytoplasm (Fig. 1g). Altogether, these results suggested that hsa_circ_102002 is upregulated in PTC tissues and cells, which may exert an important function in PTC progression.

Hsa_circ_102002 promoted epithelial-mesenchymal transformation of PTC cells
To study hsa_circ_102002's function, SNU-790 cells stably overexpressing hsa_circ_102002 and TPC-1 cells stably downregulating hsa_circ_102002 were constructed and confirmed by q-PCR (Fig. 2a). As shown in Fig. 2b, the cellular morphology of Lv-circ_102002 SNU-790 cells changed to a more obvious spindle shape, and the intercellular mass increased compared to Lv-NC SNU-790 cells.
The cellular morphology showed the opposite trend when downregulating hsa_circ_102002 in TPC-1 cells. IF analysis indicated that the expression of E-cadherin was downregulated and the level of N-cadherin was upregulated in Lv-circ_102002 SNU-790 cells. Accordingly, opposite phenomenon was observed in sh-circ_102002 TPC-1 cells (Fig. 2c). In addition, western blot analysis showed that the expression of E-cadherin was reduced and the level of Ncadherin was increased in Lv-circ_102002 SNU-790 cells, while the change trends were reversed in sh-circ_102002 TPC-1 cells. Furthermore, the levels of mesenchymal phenotype biomarkers (Vimentin, Slug and Twist) were increased in Lv-circ_102002 SNU-790 cells, and were reduced in sh-circ_102002 TPC-1 cells (Fig. 2d). All these results suggested that hsa_circ_102002 could promote epithelial-mesenchymal transformation of PTC cells.

Hsa_circ_102002 could promote cell migration in PTC
The effects of circ_102002 on cell migration were studied by wound scratch healing, transwell and western blot assays. The results indicated that the migration in Lv-circ_102002 SNU-790 cells was promoted, while inhibited in sh-circ_102002 TPC-1 cells (Fig. 3a). Moreover, the number of migrated and invaded Lv-circ_102002 SNU-790 cells was significantly increased, and that of sh-circ_102002 TPC-1 cells was markedly reduced (Fig. 3b). The protein levels of MMP2 and MMP9 were significantly upregulated in Lv-circ_102002 SNU-790 cells, whereas decreased in sh-circ_102002 TPC-1 cells (Fig. 3c). These data above suggested that circ_102002 could promote cell migration in PTC, and its downregulation exerted an inverse effect.

Circ_102002's effects on EMT and migration of PTC cells were mediated by HAS2
As mentioned above, our results revealed that circ_102002 promote EMT and cell migration in PTC by sponging miR-488-3p. Moreover, miR-488-3p could target HAS2 in PTC. Thus, whether the effects of circ_102002 were mediated by HAS2 was studied. As shown in Fig. 6a, Ecadherin expression was increased and N-cadherin level was downregulated in sh-circ_102002 TPC-1 cells, which could be restored by upregulating HAS2. Results from wound scratch healing assay indicated that the mobility of sh-circ_102002 TPC-1 cells was inhibited and HAS2 overexpression rescued the effect caused by circ_102002 inhibition (Fig. 6b). Transwell assays results showed that the number of migrated and invaded sh-circ_102002 TPC-1 cells was markedly reduced, which was restored by HAS2 overexpression (Fig. 6c). In addition, western blot analysis Fig. 4 Hsa_circ_102002 could sponge and negatively regulated miR-488-3p. a The intersection of potential miRNAs that could be sponged by hsa_circ_102002, as predicted by Circular RNA Interactome and starbase. b Relative enrichment levels of miR-488-3p, miR-625-5p and miR-1197 in Lv-NC or Lv-circ_102002 SNU-790 cells, as determined using qRT-PCR. c Relative expression levels of miR-488-3p, miR-625-5p, and miR-1197 in Lv-NC or Lv-circ_102002 SNU-790 cells and sh-NC or sh-circ_102002 TPC-1 cells, as determined using qRT-PCR. d Relative miR-488-3p levels in SNU-790 cells transfected with NC inh or miR-488-3p inh and in TPC-1 cells transfected with NC mimics or miR-488-3p mimics, as determined using qRT-PCR. e A schematic about the design of luciferase assay. The predicted seed-recognition site in the corresponding circ_102002 sequence are marked out. The circ_102002 reporter with mutant binding site was used as a negative control. f Relative luciferase activity of the circ_102002 reporter plasmid in SNU-790 cells upon miR-488-3p inh or NC inh transfection and TPC-1 cells upon miR-488-3p mimics or NC mimics transfection. The mutant circ_102002 reporter was used as a negative control. g Relative miR-488-3p levels in PTC tissues (tumor, n = 50) and normal adjacent tissues (non-tumor, n = 50), as determined using qRT-PCR. h Pearson's correlation analysis of the relative expressions between miR-488-3p and circ_102002 (mean ± SEM, *p < 0.05, **p < 0.01).
showed that protein levels of HAS2, p-FAK, and p-Akt were significantly downregulated in sh-circ_102002 TPC-1 cells, while, HAS2 overexpression restore expression levels of p-FAK and p-Akt (Fig. 6d). These findings indicated that circ_102002 affected EMT and migration of PTC cells by regulating HAS2 expression.

Circ_102002 downregulation inhibited tumor metastasis in vivo
On account of that circ_102002 inhibition suppress the EMT and cell migration in vitro, the role of circ_102002 was then studied in a mouse model. To construct metastatic mouse model, sh-NC and sh-circ_102002 TPC-1 cells with a luciferase label were intravenously injected. The luciferase activities in sh-NC and sh-circ_102002 TPC-1 cells exhibited similar prior to the cells being injected into mice (Fig. 7a). Twenty-days days post-implantation, a significant lower degree of lung metastasis were found in mice from sh-circ_102002 group, as determined by an in vivo imaging system (Fig. 7b). The corresponding lungs were excised and histopathological analysis showed that the number of nodules in the lung from sh-circ_102002 group mice was significantly less than that from sh-NC group mice (Fig. 7c). IHC results in Fig. 7d demonstrated that circ_102002 downregulation resulted in decreased expressions of HAS2, p-FAK, and p-Akt in lung tissues of sh-circ_102002 group of mice. Then, upregulated E-cadherin and decreased N-cadherin expression were determined. Altogether, these findings revealed that circ_102002 inhibition could downregulate HAS2 level, suppress the phosphorylation of FAK and Akt, regulate expressions of HAS2, E-cadherin and Ncadherin, and finally inhibit PTC metastasis in vivo (Fig. 8).

Discussion
We not only demonstrated the simulative effects of circ_102002 on EMT and cell migration, but also identified a regulatory relationship between circ_102002 and miR-488-3p/HAS2. Recently, increasing number of circRNAs are identified and becoming a research hotspot, some of which have been investigated in PTC. For example, circ-ITCH overexpression inhibits PTC cell proliferation and invasion and promotes apoptosis through the miR-22-3p/ CBL/β-catenin pathway [27]. In addition, CircZFR promotes the expression of C8orf4 through sponging miR-1261 and contributes to PTC cell proliferation and invasion [28]. These findings suggested that circRNAs usually sponge miRNAs and regulate their target genes. It was predicted that circ_102002 might sponge miR-488-3p, miR-625-5p and miR-1197. Since miR-488-3p was most enriched by circ_102002 probe using RIP assay and downregulated by circ_102002 overexpression. The present study subsequently found that circ_102002 exerts its regulatory functions by downregulating miR-8075. Based on the fact that one circRNA might interact with several molecules at the same time, it is worth further investigation whether there are other miRNAs regulated by circ_102002. As a newly identified miRNA, miR-488-3p was found as a tumor suppressor in esophageal squamous cell carcinoma and melanoma by respectively targeting ZBTB2 and PRKDC [16,17]. The features of ZBTB2 have not been reported in PTC. PRKDC encodes DNA-dependent protein kinase catalytic subunit (DNA-PKcs), and its expression was identified positively correlated with the radiosensitivity of PTC [29]. This study showed that miR-488-3p expression was downregulated in PTC and inhibited EMT and cell migration. Among the potential target genes of miR-488-3p, HAS2 stood out because of the significant inhibition efficiency on luciferase reporters caused by miR-488-3p. Therefore, HAS2 is identified as the main target of miR-488-3p in this study. Whether circ_102002/miR-488-3p affects the radiosensitivity of PTC by regulating PRKDC requires further investigations. Furthermore, some noncoding RNAs were found to act as ceRNAs for miR-488-3p. Recently, miR-488-3p has been reported to be sponged by circ_CCDC66 in Hirschsprung's disease [30], circ_0020123 in non-small cell lung cancer [31], and lncRNA PVT1 in glioma [32] etc. Therefore, it can be speculated that miR-488-3p function as regulatory hub in the progression of these diseases.
As a synthetic enzyme for hyaluronan, HAS2 affects various processes of epithelial-mesenchymal transition and cell migration [20]. HAS2 has been demonstrated to promote tumor progression in some cancers. For instance, HAS2 synthesized hyaluronic acid, enhanced ZEB1 expression and accelerated EMT in breast cancer [19]. HAS2-mediated cancer progression was also indicated in lung cancer [33] and squamous cell carcinoma [34]. Moreover, combined use of HAS2 and HYAL-1 hyaluronidase was able to predict bladder cancer diagnosis and prognosis [35]. So far, the role of HAS2 in PTC has not been reported. In the presented study, the results revealed that overexpression of HAS2 could restore the inhibition of EMT and cell migration caused by circ_102002 downregulation. This process was related to the FAK/PI3K/Akt pathway and was mediated by miRNA-488-3p. Some Fig. 5 Hsa_circ_102002 negatively regulates miR-488-3p through upregulating HAS2. a The intersection of potential genes that could be targeted by miR-488-3p, as predicted by Targetscan, mirTarBase, miRDB and starBase. b Relative luciferase activities of the 3ʹ-UTR reporter plasmids (HAS2, LRRC41, BRPF1, EIF3A, VPS4B, ZNF281, ARPP19, PKNOX1, ANAPC16, and STK4) were measured in TPC-1 cells upon transfection of miR-488-3p mimics or NC mimics. c Seeds match for miR-488-3p in the 3ʹ-UTR of HAS2. The predicted seedrecognition sites in the HAS2 mRNA sequence and the corresponding miR-488-3p sequence are marked in blue. d Relative luciferase activity of the HAS2 3ʹ-UTR reporter plasmid was measured in SNU-790 cells after expressing NC inh or miR-488-3p inh. The mutant HAS2 3ʹ-UTR reporter was used as a negative control. e Representative IHC staining of GAPDH and HAS2 in the PTC tissues (tumor, n = 50) and normal adjacent tissues. Original magnification, ×40 and ×400. f Relative HAS2 mRNA expression levels in PTC tissues (tumor, n = 50) and normal adjacent tissues (non-tumor, n = 50), as determined using qRT-PCR. g Pearson's correlation analysis of the relative expressions between circ_102002 and HAS2 mRNA, miR-488-3p and HAS2 mRNA in PTC patients. h Protein expression levels of HAS2 in Lv-NC + NC mimics, Lv-circ_102002, Lv-circ_102002+miR-488-3p mimics SNU-790 cells, and in sh-NC + NC inh, sh-circ_102002, sh-circ_102002+miR-488-3p inh TPC-1 cells, as determined using western blotting. Bands were quantified and shown in histogram (mean ± SEM, *p < 0.05, **p < 0.01). miRNAs have been reported to target HAS2 in physiopathologic process, which including miR-26b in ovarian granulosa cell apoptosis [36], miR-574 in oocyte maturation [37], and miR-101-3p in synovial sibroblasts [38] etc. However, there was no report about miRNA-regulating the expression of HAS2 during cancer development. Our results for the first time demonstrated miR-488-3p could negatively regulate HAS2 expression and form a molecular signal axis of circ_102002/miR-488-3p/HAS2 in PTC.
The emerging roles for circRNA and miRNAs as tumor suppressors or oncogenes in tumor-related pathways highlight new perspectives for the molecular pathological mechanism, and also offer novel therapeutic tools for cancer treatment [39]. Some miRNAs-based anticancer therapies have entered into clinical testing. For example, MRX34, a miR-34 mimic targeting forkhead box P1 (FOXP1) and BCL2, entered a multicenter phase I trial in 2013 for patients with primary liver cancer and small cell lung cancer. However, this trial was terminated due to immunerelated adverse events, which demands new design and efficient delivery [40]. Recently, two phase I clinical trials of the miR-29 mimics targeting histone deacetylase 4 (HDAC4) and MMPs (matrix metalloproteinases) in scleroderma patients and miR-155 inhibitors targeting SHIP (SH2-containing-inositol phosphatase) and RHOA (Ras homolog gene family,member A) in cutaneous T-cell lymphoma patients were initiated and placed great expectations. RNA therapies, when used in combination with other cancer treatments, might have synergistic effects, thereby reducing the dosage and minimizing off-target effects [41,42]. In this Fig. 6 Hsa_circ_102002/HAS2 promoted the epithelial-mesenchymal transformation and migration of PTC cells through regulating FAK/ Akt passway. a Representative images of E-cadherin and N-cadherin in sh-NC + Vector, sh-circ_102002+Vector, and sh-circ_102002 + HAS2 TPC-1 cells, as determined using IF. (red, E-cadherin; green, N-cadherin; blue, DAPI). b Wound scratch healing assay of sh-NC + Vector, sh-circ_102002+Vector, and sh-circ_102002 + HAS2 TPC-1 cells. Quantification of the wound-healing assay was shown as histograms.
c Representative migration and invasion assay images of sh-NC + Vector, sh-circ_102002+Vector, and sh-circ_102002 + HAS2 TPC-1 cells. The migrated cells were quantified and shown as histograms. d Protein expression levels of HAS2, p-FAK, FAK, p-Akt, and Akt in sh-NC + Vector, sh-circ_102002 + Vector, and sh-circ_102002 + HAS2 TPC-1 cells, as determined using western blotting. Bands were quantified and shown in histogram (mean ± SEM, *p < 0.05, **/##p < 0.01). study, the circRNA-miRNA code of circ_102002/miR-488-3p and its effects on epithelial-mesenchymal transition and cell migration were systematically elucidated. Tumor sitespecific delivery of modulators of circ_102002/miR-488-3p/HAS2 axis or combination with other approaches will offer new opportunities to design innovative therapeutic strategies for PTC.
In summary, circ_102002 was overexpressed in PTC clinical tissues and cells. EMT and cell migration were promoted by circ_102002 overexpression and inhibited by circ_102002 downregulation. Circ_102002 mainly sponged miR-488-3p and negatively regulated its expression. In addition, HAS2 was identified as a target of miR-488-3p in PTC, which was related to the FAK/PI3K/Akt pathway. Furthermore, the inhibitory effect of downregulation of circ_102002 on PTC growth was evaluated in a mouse xenograft model, which involved miR-488-3p and HAS2 regulation. These findings about the signal axis of Fig. 7 Knockdown of hsa_circ_102002 reduced lung metastasis of PTC in vivo. a Relative luciferase activity in sh-NC or sh-circ_102002 TPC-1 cells prior to the cells being injected into mice. b Tumor metastasis progression was measured by in vivo luciferase imaging at days 21 after intravenous injection of sh-NC or sh-circ_102002 TPC-1 cells. c Representative images and H&E staining analysis of lungs dissected from mice with intravenous injection of sh-NC or sh-circ_102002 TPC-1 cells. Quantification of the lung tumor nodules was shown as histograms. d IHC staining of HAS2, p-FAK, p-Akt, E-cadherin, and N-cadherin in the lungs dissected from mice with intravenous injection of sh-NC or sh-circ_102002 TPC-1 cells.

Data availability
All data generated or analyzed during this study are included in this published article. Informed consent Written informed consents were obtained from patients.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons. org/licenses/by/4.0/.