ABSTRACT
Objective:
As the most prevalent type of thyroid malignancy, papillary thyroid carcinoma (PTC) accounts for over 80% of all thyroid cancers. Circular RNAs (circRNAs) have been found to regulate multiple cancers, including PTC.
Materials and methods:
Quantitative real-time polymerase chain reaction (qRT-PCR) and western blotting were used to analyse RNA and protein levels. Fluorescence in situ hybridization (FISH) was used to detect the distribution of the target genes. Functional experiments and animal experiments were implemented to analyse the biological functions of target genes in vitro and in vivo. Luciferase reporter, RNA pulldown, RNA binding protein immunoprecipitation (RIP) and mRNA stability assays were used to probe the underlying mechanisms.
Results:
CircSEMA6Awas found to be upregulated in PTC tissues and cells, and its circular structure was verified. CircSEMA6A promotes PTC cell migration and invasion. Moreover, circSEMA6A functions as a competing endogenous RNA (ceRNA) to upregulate proline rich and Gla domain 4 (PRRG4) expression by sponging microRNA-520h (miR-520h). CircSEMA6A recruits ELAV1 to stabilize PRRG4 mRNA and drives PTC progression via PRRG4.
Conclusion:
CircSEMA6A upregulates PRRG4 by targeting miR-520h and recruiting ELAVL1 to affect the invasion and migration of PTC cells, offering insight into the molecular mechanisms of PTC.
Keywords
Papillary thyroid carcinoma; circSEMA6A; miR-520h; ELAVL1; PRRG4
INTRODUCTION
Over the last four decades, the incidence of thyroid cancer (TC) has increased across the globe (11 Prete A, Borges de Souza P, Censi S, Muzza M, Nucci N, Sponziello M. Update on Fundamental Mechanisms of Thyroid Cancer. Front Endocrinol (Lausanne). 2020;11:102.). Papillary thyroid carcinoma (PTC), as the most frequent histotype, accounts for approximately 80%-90% of all TCs (22 Deng XY, Gan XX, Feng JH, Cai WS, Wang XQ, Shen L, et al. ALDH5A1 acts as a tumour promoter and has a prognostic impact in papillary thyroid carcinoma. Cell Biochem Funct. 2021;39(2):317-25.). Surgical resection, radioiodine therapy, and levothyroxine treatment have been used to treat PTC patients (33 Wu G, Zhou W, Lin X, Sun Y, Li J, Xu H, et al. circRASSF2 Acts as ceRNA and Promotes Papillary Thyroid Carcinoma Progression through miR-1178/TLR4 Signaling Pathway. Mol Ther Nucleic Acids. 2020;19:1153-63.). PTC features an overall good prognosis, and its overall 5-year survival rate is as high as 97%. However, PTC can become life-threatening once it metastasizes or becomes invasive (44 Yao Y, Chen X, Yang H, Chen W, Qian Y, Yan Z, et al. Hsa_circ_0058124 promotes papillary thyroid cancer tumorigenesis and invasiveness through the NOTCH3/GATAD2A axis. J Exp Clin Cancer Res. 2019;38(1):318.). Thus, clarifying the molecular mechanisms underlying the migration and invasion of PTC is of great significance to the improvement of PTC treatment.
CircRNA is a type of RNA with a covalently closed loop structure. It features a stable structure with tissue specificity. Many circRNAs play important roles in the occurrence and development of multiple human cancers (55 Kristensen LS, Hansen TB, Venø MT, Kjems J. Circular RNAs in cancer: opportunities and challenges in the field. Oncogene. 2018;37(5):555-65.,66 Lei M, Zheng G, Ning Q, Zheng J, Dong D. Translation and functional roles of circular RNAs in human cancer. Mol Cancer. 2020;19(1):30.). A large body of evidence has shown that many circRNAs participate in the regulation of PTC. For instance, circ_0005273 drives the progression of PTC by modulating SOX2 expression (77 Zhang W, Zhang H, Zhao X. circ_0005273 promotes thyroid carcinoma progression by SOX2 expression. Endocr Relat Cancer. 2020;27(1):11-21.); circFNDC3B regulates PTC progression by targeting the miR-1178/TLR4 pathway (88 Wu G, Zhou W, Pan X, Sun Z, Sun Y, Xu H, et al. Circular RNA Profiling Reveals Exosomal circ_0006156 as a Novel Biomarker in Papillary Thyroid Cancer. Mol Ther Nucleic Acids. 2020;19:1134-44.); circ-PSD3 targets the miR-7-5p/METTL7B axis to aggravate PTC cell proliferation and invasion (99 Zhu J, Wang Y, Yang C, Feng Z, Huang Y, Liu P, et al. circ-PSD3 promoted proliferation and invasion of papillary thyroid cancer cells via regulating the miR-7-5p/METTL7B axis. J Recept Signal Transduct Res. 2022;42(3):251-60.); and circBACH2 serves as a regulator in PTC by sequestering miR-139-5p to regulate LMO4 expression levels (1010 Cai X, Zhao Z, Dong J, Lv Q, Yun B, Liu J, et al. Circular RNA circBACH2 plays a role in papillary thyroid carcinoma by sponging miR-139-5p and regulating LMO4 expression. Cell Death Dis. 2019;10(3):184.).
The role of circSEMA6A, the focus of this study, has rarely been reported. Its host gene, SEMA6A, can modulate multiple cancers according to previous studies. For instance, SEMA6A can act as a prognostic biomarker for glioblastoma (1111 Zhao J, Tang H, Zhao H, Che W, Zhang L, Liang P. SEMA6A is a prognostic biomarker in glioblastoma. Tumour Biol. 2015;36(11):8333-40.). In addition, SEMA6A modulates the growth and survival of human melanoma cells (1212 Loria R, Bon G, Perotti V, Gallo E, Bersani I, Baldassari P, et al. Sema6A and Mical1 control cell growth and survival of BRAFV600E human melanoma cells. Oncotarget. 2015;6(5):2779-93.). However, the roles of circSEMA6A and SEMA6Ahave never been reported in PTC.
PRRG4 was also a research object in our study. It has been found to facilitate the metastasis of breast cancer by recruiting NEDD4 and downregulating Robo1 (1313 Zhang L, Qin Y, Wu G, Wang J, Cao J, Wang Y, et al. PRRG4 promotes breast cancer metastasis through the recruitment of NEDD4 and downregulation of Robo1. Oncogene. 2020;39(49):7196-208.). Moreover, Liu and cols. noted that PRRG4 can act as a potential prognostic marker in cholangiocarcinoma (1414 Liu J, Liu W, Li H, Deng Q, Yang M, Li X, et al. Identification of key genes and pathways associated with cholangiocarcinoma development based on weighted gene correlation network analysis. PeerJ. 2019;7:e7968.). However, the role of PRRG4 has never been explored in PTC.
In this study, we decided to investigate the roles of circSEMA6A and PRRG4 in PTC and their molecular mechanisms in cell migration and invasion. It is urgent to deepen our understanding of the mechanisms underlying PTC to develop treatments for PTC patients.
MATERIALS AND METHODS
Cell line culture and vector construction
Nthy-ori 3-1, IHH-4, TPC-1 and GLAG-66 cells were purchased from WHELAB (China). These cell lines were placed in RPMI-1640 medium supplemented with 10% FBS. As verified by Cellosaurus (https://web.expasy.org/cellosaurus/), all the cell lines were free of contamination. Full-length sequences of PRRG4 were inserted into the pcDNA3.1 vector for the construction of overexpression vectors, with the vector itself as a negative control (NC). Small interfering RNAs (siRNAs) against circSEMA6A and ELAVL1 as well as si-NC were constructed for interference.
qRT-PCR
All RNAs were extracted from Nthy-ori 3-1, IHH-4, TPC-1 and GLAG-66 cells using TRIzol reagent. Afterwards, the total RNA was converted into cDNA by reverse transcription. All samples were assessed through qPCR using a qRT-PCR kit. The reaction conditions were as follows: 95°C for 10 min, 95 °C for 30s, and 60 °C for 30s. The results were quantified according to the 2-ΔΔCt method. There were three replicates per condition.
Western blotting
Total protein was isolated from PTC cells using a total protein extraction kit. After that, the samples were subjected to separation by SDS/PAGE and then transferred to polyvinylidene fluoride (PVDF) membranes. Afterwards, these membranes were blocked with 5% skim milk for 2 h at room temperature. Subsequently, they were subjected to overnight incubation with the primary antibodies at 4 °C, followed by 1 h of incubation with the secondary antibodies at room temperature. Finally, images we retaken and assessed. We used β-actin as our internal reference. The primary antibody used included anti-ELAVL1. There were three replicates per condition.
Transwell assay
Transwell chambers coated with or without Matrigel were used to test the invasion and migration of transfected PTC cells. The transfected cells were plated into the upper chamber with the addition of 0.2 mL serum-free medium, while the lower chamber was supplemented with 0.6 mL medium containing 10% FBS. After 24 h, the cells in the upper chamber were removed. The PTC cells left in the lower chambers were counted using a microscope with the help of 0.2% crystal violet solution for staining. There were three replicates per condition.
Wound healing assay
The wound healing assay was carried out to detect PTC cell migration. A pipette was used to scratch the cell layers. The average width of the wound was quantified 24 h later. There were three replicates per condition.
FISH
PTC cells were subjected to fixation in 4% PFA for 15 min, and then permeabilization of the cells was conducted using 0.5% Triton X-100 for 15 min at 4 °C. Next, the digoxigenin (DIG)-labelledcircSEMA6A probe or control probe was incubated with PTC cells for 4 h at 55 °C. The signal was detected via Hoechst-conjugated anti-digoxin antibodies. A laser confocal microscope was employed to capture images. Cell nuclei were counterstained with DAPI. There were three replicates per condition.
RIP assay
Magnetic beads were used for conjugation with anti-IgG, anti-AGO2 or anti-ELAVL1. PTC cells were lysed for the assay. The cell lysate was incubated with the conjugated magnetic beads. The precipitated RNAs were subjected to extraction and assessed by qPCR. There were three replicates per condition.
RNA pulldown assay
The structure buffer was added to 1 μg of biotin-labelled circSMEA6A/miR-520h. Then, biotinylated RNAs were heated at 95 °C for 2 min, ice-bathed for 3 min, and left for 30 min at room temperature. Afterwards, 15 μL streptavidin beads were incubated with biotin-labelled and denatured RNA for 2h at 4 °C. Three groups (input, Bio-NC, Bio-circSMEA6A/miR-520h) were set up, and PTC cells in each group were lysed. The magnetic beads were mixed with the cell lysate, followed by overnight incubation at 4 °C. After incubation, RNA was extracted, followed by qRT-PCR and western blot analyses. There were three replicates per condition.
Luciferase reporter assay
To perform the assay, pmirGLO-circSMEA6A/PRRG4 3’UTR-Wt was constructed by subcloning the sequences of circSMEA6A/PRRG4 3’UTR into the pmirGLO reporter vector. Similarly, the circSMEA6A/PRRG4 3’UTR containing mutated binding sites in the miR-520h seed region was inserted into the pmirGLO vector to create pmirGLO-circSMEA6A/PRRG4 3’UTR-Mut. PTC cells were placed into 96-well plates. These reporter vectors were cotransfected with miR-520hmimics or NC mimics into PTC cells. The vectors themselves were used as the negative controls. The luciferase activity was tested with Renilla luciferase as the internal reference. There were three replicates per condition.
Animal experiments
BALB/c nude mice (4 weeks old; female) were divided into three groups, with 3 mice in each group. These animal experiments were permitted by the Ethics Committee of the First Affiliated Hospital of Hebei North University. GLAG-66 cells (1 × 106) stably transfected with si-NC, si-circSEMA6A-1 or si-circSEMA6A-2 were subcutaneously injected into the right flanks of mice. From day 7 to day 28, the tumour size was measured every 3 days. On day 28, the mice were sacrificed. The tumours were removed and weighed.
mRNA stability assay
For this assay, 50 mMα-amanitin was used to treat PTC cells to inhibit RNA transcription. qRT-PCR was used to examinePRRG4and β-actin expression levels at 0, 6, 12, 18 and 24 h. 18S rRNA, which is not affected by α-amanitin, served as the internal reference. There were three replicates per condition.
Statistical analysis
Student’s t test and one-way/two-way analysis of variance (ANOVA) were applied to assess differences between groups. A P value lower than 0.05 was indicative of statistical significance.
RESULTS
CircSEMA6A, featuring a circular structure, is upregulated in PTC
We retrieved the Gene Expression Omnibus (GEO) dataset GSE173299 to screen circRNAs that were differentially overexpressed in PTC tissues relative to adjacent tissues. As indicated in the volcano plot, we identified three circRNAs, circRNA_103923, circRNA_405716 and circRNA_014234 (Figure 1A). The host genes of circRNA_103923, circRNA_405716 andcircRNA_014234 are SEMA6A, CHAF1A and S100A2, respectively. Hence, they are termed circSEMA6A, circCHAF1A and circS100A2. According to a literature review, circCHAF1A is overexpressed in glioma (1515 Jiang Y, Wang Z, Ying C, Hu J, Zeng T, Gao L. FMR1/circCHAF1A/miR-211-5p/HOXC8 feedback loop regulates proliferation and tumorigenesis via MDM2-dependent p53 signaling in GSCs. Oncogene. 2021;40(24):4094-110.), while circSEMA6A and circS100A2 have never been reported to be involved in diseases. Next, we conducted qRT-PCR to examine the expression levels of three candidates in a human thyroid follicular epithelial cell line (Nthy-ori 3-1) and PTC cell lines (IHH-4, TPC-1 and GLAG-66). The results revealed that circSEMA6A was most overexpressed in PTC cell lines (Figure 1B). Therefore, circSEMA6A was selected as the focus of our study. Moreover, due to the relatively high expression of circSEMA6A, TPC-1 and GLAG-66 cells were chosen for follow-up assays. Next, PCR-agarose gel electrophoresis was performed to verify the circular structure of circSEMA6A. As indicated in Figure 1C, circSEMA6A was amplified by convergent primers and divergent primers in complementary DNA (cDNA), but in genomic DNA (gDNA), it was only amplified by convergent primers. This verified the circular structure of circSEMA6A. For further verification, qRT-PCR was implemented to determine circSEMA6A and SEMA6A mRNA levels after treatment with ribonuclease R (RNase R). SEMA6A mRNA was digested by RNase R, while circSEMA6A expression remained almost unchanged (Figure 1D). The results further validated the circular structure of circSEMA6A, as the structure of circRNA is more stable than that of linear RNA. Taken together, these findings indicate that circSEMA6A, featuring a circular structure, is upregulated in PTC.
CircSEMA6A facilitates PTC progression in vitro and in vivo
Next, we performed functional assays to investigate the influence of circSEMA6A on PTC cell migration and invasion. First, we performed qRT-PCR to evaluate the expression of circSEMA6Aafter si-circSEMA6A-1/2/3 transfection. These plasmids suppressed the expression of circSEMA6A (Figure 2A). Owing to their higher efficiency, si-circSEMA6A-1 and si-circSEMA6A-2 were chosen for the following experiments. Afterwards, we performed a wound healing assay in PTC cells and found that the wound width was increased after the transfection of si-circSEMA6A-1/2 (Figure 2B). Transwell assays were implemented in PTC cells to detect the impacts of circSEMA6A on cell migration and invasion. The number of migrated cells and invaded cells was reduced after circSEMA6A depletion (Figure 2C). In addition, we performed qRT-PCR to evaluate the mRNA levels of the invasion-related proteins MMP2 and MMP9 after knockdown of circSEMA6A. From the results, we found that MMP2 and MMP9 were downregulated by circSEMA6Aknockdown (Figure 2D). Additionally, we conducted animal experiments to determine the functions of circSEMA6A in vivo. As displayed by the growth curve, tumour growth was suppressed after circSEMA6A knockdown. After the resection of tumours from mice, we measured the volume and weight of tumours. CircSEMA6A interference led to a decrease in tumour volume and weight. In summary, circSEMA6A aggravates PTC cell progression in vitro and in vivo.
CircSEMA6A sponges miR-520h via ceRNA mode
We next probed the regulatory mechanisms of circSEMA6A. We first conducted a FISH assay to detect the subcellular distribution of circSEMA6A in PTC cells. As indicated in Figure 3A, circSEMA6A was located in both the cytoplasm and nucleus of PTC cells but mainly in the cytoplasm. Moreover, a RIP assay was performed to show the enrichment of circSEMA6A in AGO2 complexes (Figure 3B). This indicated its existence in the RNA-induced silencing complex (RISC), which is the major workplace for miRNAs. Based on the results of Figure 3A and Figure 3B, we hypothesized that circSEMA6A functions as a ceRNA. Subsequently, we utilized the starBase database (http://starbase.sysu.edu.cn/) and CircInteractome (https://circinteractome.nia.nih.gov/) to forecast miRNA binding with circSEMA6A. As indicated in the Venn diagram, miR-520h was predicted as a binding partner by both starBase and CircInteractome. Moreover, sequences of the binding site between miR-520h and circSEMA6A were obtained (Figure 3C). To prove the interaction between circSEMA6A and miR-520h, we carried out an RNA pulldown assay to determine the enrichment of circSEMA6A in samples pulled down by Bio-miR-520h-Wt and Bio-miR-520h-Mut (sequences shown in Figure 3C). The results revealed that circSEMA6A was most abundant in the Bio-miR-520h-Wt group, verifying the interaction (Figure 3D). We conducted qRT-PCR to probe the miR-520h expression level after miR-520h mimic transfection and found that miR-520h expression was upregulated (Figure 3E). Afterwards, a luciferase reporter assay was conducted to further validate the interaction between circSEMA6A and miR-520h. The results showed that the cotransfection of miR-520h mimics reduced pmirGLO-circSEMA6A-Wt luciferase activity, while that of pmirGLO-circSEMA6A-Mut (designed by the sequences shown in Figure 3C) remained unchanged (Figure 3F). We then implemented qRT-PCR to investigate the expression of miR-520h after the knockdown of circSEMA6A. The results revealed that miR-520h expression remained unchanged after circSEMA6Aknockdown (Figure 3G). Next, we carried out assays to examine the influence of miR-520h on PTC cell progression. Wound healing, Transwell and qRT-PCR assays revealed that PTC cell migration and invasion were hampered after miR-520h overexpression. Taken together, these findings indicate that circSEMA6A sponges miR-520h via the ceRNA mode.
CircSEMA6A sponges miR-520h to upregulate PRRG4 expression
We next probed the downstream mRNA targeted by miR-520h. We used starBase to predict the mRNA binding to miR-520h under certain conditions (CLIP-Data>=5, pan-Cancer>=10). PRRG4, ATP1B1 and GCNT2 were selected as candidates (Figure 4A). According to a literature review, PRRG4 is highly expressed in breast cancer and promotes its metastasis (1313 Zhang L, Qin Y, Wu G, Wang J, Cao J, Wang Y, et al. PRRG4 promotes breast cancer metastasis through the recruitment of NEDD4 and downregulation of Robo1. Oncogene. 2020;39(49):7196-208.); ATP1B1 is expressed at low levels in PTC tissues (1616 Rosignolo F, Sponziello M, Durante C, Puppin C, Mio C, Baldan F, et al. Expression of PAX8 Target Genes in Papillary Thyroid Carcinoma. PloS One. 2016;11(6):e0156658.), which is not in line with our research and thus was excluded; and GCNT2 can induce the epithelial-mesenchymal transition of oesophageal squamous cell carcinoma cells and enhance cell migration and invasion (1717 Peng F, He Q, Cheng C, Pan J. GCNT2 induces epithelial-mesenchymal transition and promotes migration and invasion in esophageal squamous cell carcinoma cells. Cell Biochem Funct. 2019;37(1):42-51.). There is no research report on PRRG4 and GCNT2 in PTC. Therefore, we first analysed the expression of PRRG4 and GCNT2 in PTC tissues by using bioinformatics. UALCAN predictions (http://ualcan.path.uab.edu/) showed that only PRRG4 was significantly overexpressed in PTC tissues (Figure 4B). Furthermore, we implemented qRT-PCR to examine PRRG4 and GCNT2 expression levels in TPC-1 cells after the knockdown of circSEMA6A and found that PRRG4 expression was much lower than GCNT2 expression (Figure 4C). Based on the results of Figure 4B-C, we selected PRRG4 as the focus of our study. We obtained the binding sites between PRRG4 and miR-520h from starBase and then performed an RNA pulldown assay in PTC cells to evaluate the binding between PRRG4 and miR-520h. The results showed that the PRRG4 3’UTR was more enriched in the Bio-miR-520h-Wt group than in the Bio-miR-520h-Mut group, verifying the interaction (Figure 4D). Afterwards, a luciferase reporter assay was conducted in PTC cells to further validate the interaction between PRRG4 and miR-520h. The results revealed that the cotransfection of miR-520h mimics reduced pmirGLO-PRRG4 3’UTR-Wt luciferase activity, while that of pmirGLO-PRRG4 3’UTR-Mut (designed by the sequences shown in Figure 4D) remained unchanged (Figure 4E). Next, we conducted qRT-PCR in PTC cells to examine the expression of PRRG4 after transfection with the miR-520h inhibitor. The results showed that PRRG4 was overexpressed with miR-520h inhibitor treatment, indicating the successful construction of the miR-520h inhibitor (Figure 4F). Subsequently, a RIP assay in PTC cells was used to detect the enrichment of PRRG4 in the anti-AGO2 group after the knockdown of miR-520h. PRRG4 enrichment was decreased after miR-520h knockdown (Figure 4G). Next, the regulatory axis of circSEMA6A, PRRG4 and miR-520h was explored by rescue experiments. We carried out qRT-PCR to evaluate PRRG4 expression in PTC cells after transfection with si-NC, si-circSEMA6A-1, si-circSEMA6A-1+inhibitor-NC or si-circSEMA6A-1+miR-520h inhibitor. The results showed that PRRG4 expression was decreased withcircSEMA6A knockdown and was then partially rescued by miR-520h interference (Figure 4H). Taken together, circSEMA6A sponges miR-520h to upregulate PRRG4 expression.
CircSEMA6A recruits ELAVL1 to stabilize PRRG4 mRNA
As indicated in the rescue experiments, PRRG4 expression was only partially rescued by miR-520h interference. Therefore, we speculated that circSEMA6A can also regulate the expression of PRRG4 through other mechanisms. It has been reported that circRNAs regulate downstream targets by recruiting RNA binding proteins (RBPs) (1818 Liu Z, Wang Q, Wang X, Xu Z, Wei X, Li J. Circular RNA cIARS regulates ferroptosis in HCC cells through interacting with RNA binding protein ALKBH5. Cell Death Discov. 2020;6:72.,1919 Chen RX, Chen X, Xia LP, Zhang JX, Pan ZZ, Ma XD, et al. N(6)-methyladenosine modification of circNSUN2 facilitates cytoplasmic export and stabilizes HMGA2 to promote colorectal liver metastasis. Nat Commun. 2019;10(1):4695.). Hence, we conducted experiments to study whether circSEMA6A regulates the expression of PRRG4 by recruiting RBPs. First, we performed bioinformatics analysis to identify the RBPs interacting with circSEMA6A. ELAVL1 was predicted by both starBase and circAtlas 2.0 (http://circatlas.biols.ac.cn/), as indicated in the Venn diagram (Figure 5A). ELAVL1 can stabilize mRNA by binding to its 3’UTR (2020 Chen X, Li A, Sun BF, Yang Y, Han YN, Yuan X, et al. 5-methylcytosine promotes pathogenesis of bladder cancer through stabilizing mRNAs. Nat Cell Biol. 2019;21(8):978-90.). We performed a RIP assay to investigate the interaction between ELAVL1 and circSEMA6A. The results revealed that circSEMA6A was abundant in the samples precipitated by anti-ELAVL1 (Figure 5B). For further validation, we implemented an RNA pulldown assay and found that ELAVL1 protein was enriched in the Bio-circSEMA6A group (Figure 5C). Afterwards, we performed RIP and RNA pulldown to probe the binding of the PRRG4 3’UTR to ELAVL1. The results revealed that the PRRG4 3’UTR was abundant in the samples precipitated by anti-ELAVL1 and that the ELAVL1 protein was enriched in PRRG4 3’UTR-sense (Figure 5D-E). We performed a RIP assay to detect the influence of circSEMA6A on the interaction between the PRRG4 3’UTR and ELAVL1. CircSEMA6A decreased the ELAVL1-induced enrichment of the PRRG4 3’UTR (Figure 5F). Next, qRT-PCR was implemented to examine ELAVL1 expression levels after interference with circSEMA6A. CircSEMA6A overexpression did not affect the expression of ELAVL1 (Figure 5G). We performed qRT-PCR to investigate the expression of ELAVL1 after si-ELAVL1-1/2/3 transfection. These plasmids suppressed the expression of ELAVL1 (Figure 5H). Owing to the higher efficiency, we chose si-ELAVL1-1 and si-ELAVL1-2 for the ensuing assays. Next, we conducted qRT-PCR to examine the expression of PRRG4 after the knockdown of ELAVL1. ELAVL1 knockdown led to a decrease inPRRG4 expression (Figure 5I). Finally, we conducted qRT-PCR to assess the stability of PRRG4 in PTC cells treated withα-amanitin. The results showed that the mRNA stability of PRRG4 was decreased after ELAVL1knockdown (Figure 5J). Taken together, these findings indicate that circSEMA6A recruits ELAVL1 to stabilize PRRG4 mRNA.
CircSEMA6A regulates PTC cell migration and invasion via PRRG4
Our results clearly indicated that circSEMA6A upregulates PRRG4 by targeting miR-520h and recruiting ELAVL1, so we next designed rescue experiments for further validation. We carried out qRT-PCR to assess PRRG4 expression levelsafterpcDNA3.1-PRRG4 transfection. pcDNA3.1-PRRG4enhanced the expression of PRRG4 (Figure 6A). We next performed a wound healing assay in TPC-1 cells after the transfection of si-NC, si-circSEMA6A-1, si-circSEMA6A-1+pcDNA3.1 or si-circSEMA6A-1+pcDNA3.1-PRRG4 and found that the increase in wound width induced by circSEMA6A was counteracted after the transfection of pcDNA3.1-PRRG4 (Figure 6B). Transwell assays were implemented in TPC-1 cells to determine cell migration and invasion after transfection. The number of migrated cells and invaded cells was reduced by circSEMA6A knockdown, and this change was reversed by PRRG4 overexpression (Figure 6C). We then performed qRT-PCR to evaluate the mRNA levels of MMP2 and MMP9 after transfection. From the results, we found that MMP2 and MMP9 were downregulated by circSEMA6A knockdown, this effect was reversed byPRRG4 overexpression (Figure 6D). In summary, circSEMA6A regulates PTC cell migration and invasion via PRRG4.
DISCUSSION
Previous studies indicated that circRNAs can modulate PTC cell migration and invasion. CircFAT1 (e2) is overexpressed in PTC cells and facilitates PTC cell migration and invasion (2121 Liu J, Li H, Wei C, Ding J, Lu J, Pan G, et al. circFAT1(e2) Promotes Papillary Thyroid Cancer Proliferation, Migration, and Invasion via the miRNA-873/ZEB1 Axis. Comput Math Methods Med. 2020;2020:1459368.). Circ_007293 is abundant in PTC cells and drives cell migration and invasion (2222 Lin Q, Qi Q, Hou S, Chen Z, Jiang N, Zhang L, et al. Exosomal circular RNA hsa_circ_007293 promotes proliferation, migration, invasion, and epithelial-mesenchymal transition of papillary thyroid carcinoma cells through regulation of the microRNA-653-5p/paired box 6 axis. Bioengineered. 2021;12(2):10136-49.). In our study, we used the GEO database to screen out circRNAs overexpressed in PTC tissues. Afterwards, we carried out qRT-PCR to evaluate the expression of circRNAs in PTC cell lines and identified circSEMA6A, which was overexpressed in PTC cell lines. Furthermore, we used wound healing, Transwell and qRT-PCR assays to validate that circSEMA6A promotes PTC cell migration and invasion. Animal experiments verified the functions of circSEMA6A in vivo.
Previously, ceRNA networks have been widely explored in PTA. Circ-PRKCI has been reported to propel PTC cell progression via themiR-335/E2F3 ceRNA axis (2323 Liu Y, Chen G, Wang B, Wu H, Zhang Y, Ye H. Silencing circRNA protein kinase C iota (circ-PRKCI) suppresses cell progression and glycolysis of human papillary thyroid cancer through circ-PRKCI/miR-335/E2F3 ceRNA axis. Endocr J. 2021;68(6):713-27.); circFOXM1 plays a regulatory role in PTC by sequestering miR-1179 to regulate the expression of HMGB1 (2424 Ye M, Hou H, Shen M, Dong S, Zhang T. Circular RNA circFOXM1 Plays a Role in Papillary Thyroid Carcinoma by Sponging miR-1179 and Regulating HMGB1 Expression. Mol Ther Nucleic Acids. 2020;19:741-50.); and circRPS28 regulates cell proliferation and migration by functioning as a ceRNA to target the miR-345-5p/FZD8 regulatory axis (2525 Mao Y, Huo Y, Li J, Zhao Y, Wang Y, Sun L, et al. circRPS28 (hsa_circ_0049055) is a novel contributor for papillary thyroid carcinoma by regulating cell growth and motility via functioning as ceRNA for miR-345-5p to regulate frizzled family receptor 8 (FZD8). Endocr J. 2021;68(11):1267-81.). In line with previous studies, circSEMA6A was proven to function as a ceRNA based on the results of FISH and RIP assays. We used bioinformatics to predict the downstream target miRNAs and mRNAs of circSEMA6A, miR-520h and PRRG4. MiR-520h can be sponged by PART1 to overexpress CTNNB1, thereby regulating the development of colorectal cancer (2626 Zhou T, Wu L, Ma N, Tang F, Zong Z, Chen S. LncRNA PART1 regulates colorectal cancer via targeting miR-150-5p/miR-520h/CTNNB1 and activating Wnt/β-catenin pathway. Int J Biochem Cell Biol. 2020;118:105637.). CircCSPP1 enhances prostate cancer cell progression by sequestering miR-520h to regulate EGR1 (2727 Lu J, Zhong C, Luo J, Shu F, Lv D, Liu Z, et al. HnRNP-L-regulated circCSPP1/miR-520h/EGR1 axis modulates autophagy and promotes progression in prostate cancer. Mol Ther Nucleic Acids. 2021;26:927-44.). MiR-520h overexpression hampers pancreatic cancer cell migration and invasion and is targeted by LINC00657 to modulate CKS1B in pancreatic cancer cells (2828 Li P, Wang H, Tang Y, Sun S, Ma Y, Xu Y, et al. Knockdown of LINC00657 inhibits the viability, migration and invasion of pancreatic cancer cells by regulating the miR-520h/CKS1B axis. Exp Ther Med. 2021;22(4):1142.). PRRG4 is the downstream mRNA of miR-1299 in oesophageal squamous cell carcinoma (2929 Fang J, Ji WH, Wang FZ, Xie TM, Wang L, Fu ZF, et al. Circular RNA hsa_circ_0000700 promotes cell proliferation and migration in Esophageal Squamous Cell Carcinoma by sponging miR-1229. J Cancer. 2021;12(9):2610-23.). In this study, we carried out mechanistic assays, proving that circSEMA6A acts as a sponge for miR-520h and that miR-520h can bind to PRRG4. Functional experiments proved that miR-520h overexpression hinders PTC cell migration and invasion. Rescue experiments revealed that miR-520h knockdown partially rescued the inhibitory effects of circSEMA6A knockdown on PRRG4 expression and that circSEMA6A promoted PTC cell migration and invasion via PRRG4. In summary, circSEMA6A functions as a ceRNA to upregulatePRRG4 expression by sequestering miR-520h, thereby facilitating PTC cell migration and invasion.
The potential RBPs of circSEMA6A and ELAVL1 were predicted by bioinformatics. It has been reported that ELAVL1 interacts with TIMM44 mRNA to regulate its stability in ovarian cancer (3030 Yu X, Li Y, Ding Y, Zhang H, Ding N, Lu M. HuR Promotes Ovarian Cancer Cell Proliferation by Regulating TIMM44 mRNA Stability. Cell Biochem Biophys. 2020;78(4):447-53.). In addition, ELAVL1 can interact with thePTBP1 3’UTR to enhance the stabilization of PTBP1incolorectal cancer (3131 Chen J, Wu Y, Luo X, Jin D, Zhou W, Ju Z, et al. Circular RNA circRHOBTB3 represses metastasis by regulating the HuR-mediated mRNA stability of PTBP1 in colorectal cancer. Theranostics. 2021;11(15):7507-26.). We investigated whether ELAVL1 can modulate the stability of PRRG4 mRNA in PTC cells. RIP and RNA pulldown assays were used to detect the interaction between ELAVL1 and circSEMA6A and ELAVL1 and the PRRG4 3’UTR. We then performed a RIP assay and proved that circSEMA6A inhibits the interaction between ELAVL1 and the PRRG4 3’UTR. As indicated by the mRNA stability assay, ELAVL1 knockdown led to a decrease in PRRG4 mRNA stability.
In conclusion, circSEMA6A upregulates PRRG4 by targeting miR-520h and recruiting ELAVL1 to promote PTC cell invasion and migration. The current study identified a novel circRNA whose host gene is SEMA6A, which promotes PTC cell progression. Furthermore, our study first investigated the mechanisms by whichcircSEMA6A exerts its functions. The findings in our study may aid in the development of PTC treatments. Investigating the underlying molecular mechanisms of circSEMA6Ain PTC might contribute to the treatment of PTC patients. However, clinicopathological analysis is needed to explore its clinical applicability. We will study the association of circSEMA6A with PTC patients in our future study.
-
Ethic approval: these animal experiments were permitted by the Ethics Committee of the First Affiliated Hospital of Hebei North University.
Acknowledgement:
we appreciate all the participants who provide supports for the study.
REFERENCES
-
1Prete A, Borges de Souza P, Censi S, Muzza M, Nucci N, Sponziello M. Update on Fundamental Mechanisms of Thyroid Cancer. Front Endocrinol (Lausanne). 2020;11:102.
-
2Deng XY, Gan XX, Feng JH, Cai WS, Wang XQ, Shen L, et al. ALDH5A1 acts as a tumour promoter and has a prognostic impact in papillary thyroid carcinoma. Cell Biochem Funct. 2021;39(2):317-25.
-
3Wu G, Zhou W, Lin X, Sun Y, Li J, Xu H, et al. circRASSF2 Acts as ceRNA and Promotes Papillary Thyroid Carcinoma Progression through miR-1178/TLR4 Signaling Pathway. Mol Ther Nucleic Acids. 2020;19:1153-63.
-
4Yao Y, Chen X, Yang H, Chen W, Qian Y, Yan Z, et al. Hsa_circ_0058124 promotes papillary thyroid cancer tumorigenesis and invasiveness through the NOTCH3/GATAD2A axis. J Exp Clin Cancer Res. 2019;38(1):318.
-
5Kristensen LS, Hansen TB, Venø MT, Kjems J. Circular RNAs in cancer: opportunities and challenges in the field. Oncogene. 2018;37(5):555-65.
-
6Lei M, Zheng G, Ning Q, Zheng J, Dong D. Translation and functional roles of circular RNAs in human cancer. Mol Cancer. 2020;19(1):30.
-
7Zhang W, Zhang H, Zhao X. circ_0005273 promotes thyroid carcinoma progression by SOX2 expression. Endocr Relat Cancer. 2020;27(1):11-21.
-
8Wu G, Zhou W, Pan X, Sun Z, Sun Y, Xu H, et al. Circular RNA Profiling Reveals Exosomal circ_0006156 as a Novel Biomarker in Papillary Thyroid Cancer. Mol Ther Nucleic Acids. 2020;19:1134-44.
-
9Zhu J, Wang Y, Yang C, Feng Z, Huang Y, Liu P, et al. circ-PSD3 promoted proliferation and invasion of papillary thyroid cancer cells via regulating the miR-7-5p/METTL7B axis. J Recept Signal Transduct Res. 2022;42(3):251-60.
-
10Cai X, Zhao Z, Dong J, Lv Q, Yun B, Liu J, et al. Circular RNA circBACH2 plays a role in papillary thyroid carcinoma by sponging miR-139-5p and regulating LMO4 expression. Cell Death Dis. 2019;10(3):184.
-
11Zhao J, Tang H, Zhao H, Che W, Zhang L, Liang P. SEMA6A is a prognostic biomarker in glioblastoma. Tumour Biol. 2015;36(11):8333-40.
-
12Loria R, Bon G, Perotti V, Gallo E, Bersani I, Baldassari P, et al. Sema6A and Mical1 control cell growth and survival of BRAFV600E human melanoma cells. Oncotarget. 2015;6(5):2779-93.
-
13Zhang L, Qin Y, Wu G, Wang J, Cao J, Wang Y, et al. PRRG4 promotes breast cancer metastasis through the recruitment of NEDD4 and downregulation of Robo1. Oncogene. 2020;39(49):7196-208.
-
14Liu J, Liu W, Li H, Deng Q, Yang M, Li X, et al. Identification of key genes and pathways associated with cholangiocarcinoma development based on weighted gene correlation network analysis. PeerJ. 2019;7:e7968.
-
15Jiang Y, Wang Z, Ying C, Hu J, Zeng T, Gao L. FMR1/circCHAF1A/miR-211-5p/HOXC8 feedback loop regulates proliferation and tumorigenesis via MDM2-dependent p53 signaling in GSCs. Oncogene. 2021;40(24):4094-110.
-
16Rosignolo F, Sponziello M, Durante C, Puppin C, Mio C, Baldan F, et al. Expression of PAX8 Target Genes in Papillary Thyroid Carcinoma. PloS One. 2016;11(6):e0156658.
-
17Peng F, He Q, Cheng C, Pan J. GCNT2 induces epithelial-mesenchymal transition and promotes migration and invasion in esophageal squamous cell carcinoma cells. Cell Biochem Funct. 2019;37(1):42-51.
-
18Liu Z, Wang Q, Wang X, Xu Z, Wei X, Li J. Circular RNA cIARS regulates ferroptosis in HCC cells through interacting with RNA binding protein ALKBH5. Cell Death Discov. 2020;6:72.
-
19Chen RX, Chen X, Xia LP, Zhang JX, Pan ZZ, Ma XD, et al. N(6)-methyladenosine modification of circNSUN2 facilitates cytoplasmic export and stabilizes HMGA2 to promote colorectal liver metastasis. Nat Commun. 2019;10(1):4695.
-
20Chen X, Li A, Sun BF, Yang Y, Han YN, Yuan X, et al. 5-methylcytosine promotes pathogenesis of bladder cancer through stabilizing mRNAs. Nat Cell Biol. 2019;21(8):978-90.
-
21Liu J, Li H, Wei C, Ding J, Lu J, Pan G, et al. circFAT1(e2) Promotes Papillary Thyroid Cancer Proliferation, Migration, and Invasion via the miRNA-873/ZEB1 Axis. Comput Math Methods Med. 2020;2020:1459368.
-
22Lin Q, Qi Q, Hou S, Chen Z, Jiang N, Zhang L, et al. Exosomal circular RNA hsa_circ_007293 promotes proliferation, migration, invasion, and epithelial-mesenchymal transition of papillary thyroid carcinoma cells through regulation of the microRNA-653-5p/paired box 6 axis. Bioengineered. 2021;12(2):10136-49.
-
23Liu Y, Chen G, Wang B, Wu H, Zhang Y, Ye H. Silencing circRNA protein kinase C iota (circ-PRKCI) suppresses cell progression and glycolysis of human papillary thyroid cancer through circ-PRKCI/miR-335/E2F3 ceRNA axis. Endocr J. 2021;68(6):713-27.
-
24Ye M, Hou H, Shen M, Dong S, Zhang T. Circular RNA circFOXM1 Plays a Role in Papillary Thyroid Carcinoma by Sponging miR-1179 and Regulating HMGB1 Expression. Mol Ther Nucleic Acids. 2020;19:741-50.
-
25Mao Y, Huo Y, Li J, Zhao Y, Wang Y, Sun L, et al. circRPS28 (hsa_circ_0049055) is a novel contributor for papillary thyroid carcinoma by regulating cell growth and motility via functioning as ceRNA for miR-345-5p to regulate frizzled family receptor 8 (FZD8). Endocr J. 2021;68(11):1267-81.
-
26Zhou T, Wu L, Ma N, Tang F, Zong Z, Chen S. LncRNA PART1 regulates colorectal cancer via targeting miR-150-5p/miR-520h/CTNNB1 and activating Wnt/β-catenin pathway. Int J Biochem Cell Biol. 2020;118:105637.
-
27Lu J, Zhong C, Luo J, Shu F, Lv D, Liu Z, et al. HnRNP-L-regulated circCSPP1/miR-520h/EGR1 axis modulates autophagy and promotes progression in prostate cancer. Mol Ther Nucleic Acids. 2021;26:927-44.
-
28Li P, Wang H, Tang Y, Sun S, Ma Y, Xu Y, et al. Knockdown of LINC00657 inhibits the viability, migration and invasion of pancreatic cancer cells by regulating the miR-520h/CKS1B axis. Exp Ther Med. 2021;22(4):1142.
-
29Fang J, Ji WH, Wang FZ, Xie TM, Wang L, Fu ZF, et al. Circular RNA hsa_circ_0000700 promotes cell proliferation and migration in Esophageal Squamous Cell Carcinoma by sponging miR-1229. J Cancer. 2021;12(9):2610-23.
-
30Yu X, Li Y, Ding Y, Zhang H, Ding N, Lu M. HuR Promotes Ovarian Cancer Cell Proliferation by Regulating TIMM44 mRNA Stability. Cell Biochem Biophys. 2020;78(4):447-53.
-
31Chen J, Wu Y, Luo X, Jin D, Zhou W, Ju Z, et al. Circular RNA circRHOBTB3 represses metastasis by regulating the HuR-mediated mRNA stability of PTBP1 in colorectal cancer. Theranostics. 2021;11(15):7507-26.
Publication Dates
-
Publication in this collection
07 June 2024 -
Date of issue
2024
History
-
Received
01 Dec 2021 -
Accepted
07 Sept 2022