- Open Access
CircPVT1 promotes proliferation of lung squamous cell carcinoma by binding to miR-30d/e
Journal of Experimental & Clinical Cancer Research volume 40, Article number: 193 (2021)
Circular RNAs (circRNAs) are a new type of extensive non-coding RNAs that regulate the activation and progression of different human diseases, including cancer. However, information on the underlying mechanisms and clinical significance of circRNAs in lung squamous cell carcinoma (LUSC) remains scant.
The expression profile of RNAs in 8 LUSC tissues, and 9 healthy lung tissues were assayed using RNA sequencing (RNA-seq) techniques. Real-time quantitative polymerase chain reaction (qRT-PCR) was used to profile the expression of circPVT1 and its relationship with the prognosis of LUSC, i.e., survival analysis. Moreover, in vitro and in vivo experiments were performed to evaluate the impacts of circPVT1 on the growth of tumors. RNA pull-down tests, mass spectrometry, dual-luciferase reporter assessment, and RNA immune-precipitation tests were further conducted to interrogate the cross-talk between circPVT1, HuR, or miR-30d/e in LUSC.
Our data showed that circPVT1 was upregulated in LUSC tissues, serum, and cell lines. LUSC patients with higher circPVT1 expression exhibited shorter survival rates. The in vivo and in vitro data revealed that circPVT1 promotes the proliferation of LUSC cells. Additionally, mechanistic analysis showed that HuR regulated circPVT1. On the other hand, circPVT1 acted as a competing endogenous RNA (ceRNA) of miR-30d and miR-30e in alleviating the suppressive influences of miR-30d and miR-30e on its target cyclin F (CCNF).
CircPVT1 promotes LUSC progression via HuR/circPVT1/miR-30d and miR-30e/CCNF cascade. Also, it acts as a novel diagnostic biomarker or treatment target of individuals diagnosed with LUSC.
Lung cancer (LC) accounts for the highest rate of cancer incidence and mortalities across the globe. Notably, 2 types of LCs have been reported, including non-small cell lung cancer (NSCLC) and small cell lung cancer (SCLC), accounting for up to 85, and 15% of all cases, respectively . NSCLC patients present with dismal clinical outcomes, where, only 15% survive for 5 years or more . The most prevalent types of NSCLC include lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC), which accounts for approximately 90% of all cases. Moreover, LUSC has been associated with annual fatalities of more than 400, 000 . Unlike LUAD, the spectrum of targeted therapies achievable for LUSC is limited. To enhance the prognosis of LUSC, there is a desperate and urgent need for a profound understanding of the molecular mechanisms of LUSC and the development of individualized treatment strategies.
Circular RNAs (circRNAs) are a new type of extensive non-coding RNAs, regulating pivotal functions in distinct kinds of cancers. They generate non-polar or non- polyadenylated tailed covalently closed continuous looped structures, a phenomenon that enhances their stability and abundance compared to their well-established linear coding RNAs of similar cellular genes [4, 5]. In addition to tissues, circRNAs have been documented in various body fluids, including urine, saliva, and blood, indicating that they can be applied as non-invasive circulating diagnostic biomarkers for cancer. Previous reports showed that circRNAs modulate the expression of genes during transcription, post-transcription, and translation. The circRNAs in cytoplasm possess microRNA (miRNA) response elements (MREs), which competitively combined with miRNAs, downregulating their expression, reducing their function, hence promoting expressions of miRNAs targeted genes. This implies that circRNAs are linked to diverse biological processes promoting the progression of cancer, i.e., cell proliferation, migration, infiltration, and apoptosis.
Herein, we constructed RNA expression data from LUSC tissues, then established numerous candidates of circRNAs (at least 2 unique back spliced reads). Subsequently, we characterized one circRNA derived from the PVT1 locus (circPVT1), which is often up-regulated in LUSC individuals. Functional assessments revealed that HuR regulated the expression level of circPVT1 and circPVT1 improved cell proliferation via sponging miR-30d/e members. Moreover, circPVT1 was documented as an independent prognostic biomarker for survival in individuals diagnosed with LUSC.
For RNA-seq analysis, total RNAs from 8 LUSC tissues and 9 healthy lung tissues, including 8 matched vicinal healthy tissues and 1 other healthy lung tissue, were acquired from the subjects. Cancer samples and healthy tissues were collected from individuals diagnosed with LUSC admitted to the second affiliated hospital of Xi’an Jiaotong University (Xi’an China). All the participants recruited had not undergone previous surgery, radiotherapy, chemotherapy, targeted or immunotherapy and the tissues were examined histologically. Clinically pathological features comprising age, differentiation grade, gender, site of the tumor, TNM stage and tumor size are shown in Table 1. The patients were followed-up between the time of surgery and disease progression date, demise, or the last clinical assessment. All the patients provided informed consent. Besides, this study was approved by the Clinical Research Ethics Committee of the second affiliated hospital of Xi’an Jiaotong University.
For the generation of the RNA library, processing of the total RNA (tRNA) was conducted using the TruSeq Stranded Total RNA with Ribo-Zero Gold (Illumina, USA). Following the removal of rRNA, preparation of the strand-distinct RNA-seq libraries was conducted based on the manufacturer’s instructions. Briefly, fragmentation of the ribosome free RNA was conducted. Thereafter, 1st strand cDNA was generated from the fragmented RNA, followed by the preparation of the 2nd strand cDNA using the random hexamers, and dUTP, respectively. Subsequently, end repairing of the generated cDNA fragments was conducted, followed by poly-A tailing, and eventually ligation with indexed adapters. Purification of these ligated cDNA libraries was performed, then uracil DNA glycosylase treatment was applied to degrade the 2nd strand cDNA. Library enrichment of the 1st strand purified cDNA products was performed via 13–16 amplification cycles in PCR. Library quality was evaluated using a Bioanalyzer 2200 (Agilent, Santa Clara, CA), then sequencing was conducted on a HiSeq X machine (Illumina, San Diego, CA) via a run of 150 bp pair-ends.
Cell culture and treatments
The A549, H520, H226, SKMES-1, and H1270 cells were obtained from American Type Culture Collection (ATCC). Then, the cells were incubated in RPMI1640, enriched with fetal bovine serum (FBS) (10%) at 37 °C and 5% CO2. Blocking of transcription was accomplished by the introduction of DMSO (2 mg/ml) (Sigma-Aldrich, St. Louis, MO, USA), which was the control of the cell growth medium.
RNA purification and qRT–PCR
The TRIzol reagent (Life Technologies, Carlsbad, CA) was used in the purification of tRNA from the complete-cell lysates. Afterward, quantification of the mature miRNA was conducted via the TaqMan MicroRNA assays (Life Technologies) as an internal standard. On the other hand, cDNA generation from the 500 ng tRNA was performed using the PrimeScript RT Master Mix kit (Takara, Dalian, China) to quantify the mRNA and circRNA. qRT-PCR tests were conducted using the SYBR Premix Ex Taq II (Takara).
Divergent primers complementary to the circRNA distal ends were adopted in the quantification of the circRNA. The absolute RNA quantity was determined, the pure PCR products generated from amplification of cDNA complementary to the circPVT1 and purified, were serially diluted to form a standard curve.
The full-length sequence of circPVT1 was synthesized by chemical synthesis, then the sequence was constructed into the GV486 vector (Genechem Co.,Ltd., Shanghai) using the two restriction sites of KpnI / BamHI. There are two intron sequences on both sides of the KpnI / BamHI digestion site, which can promote the formation of circRNA between the two intron sequences. The full-length sequence of HuR and CCNF was synthesized, then the sequence was constructed into the PCDNA3.1vector (Jinkairui, Wuhan, China).
RiboBio (Guangzhou, China) synthesized the siRNA (Additional file 1), miRNA mimics, and miRNA inhibitors. Transfection of the cells was conducted using the Lipofectamine RNAiMax generated from Life Technologies.
The H520, H2170 cell proliferation was assessed using the CCK-8 kit (Doindo, Japan). Triplicates of 3.5 × 103 transfected cells were incubated in 100 ml in 96-well dishes. At 0 h, 24 h, 48 h and 72 h, 10 ml of the CCK-8 reagent was added to every well, followed by 2 h incubation at 37 °C. The automatic microplate reader (Synergy4; BioTek, Winooski, VT, USA) was set at 450 nm to acquire the optical density of each well.
The cell-light EdU DNA Cell Proliferation Kit (obtained from RiboBio, Shanghai, China) was used to inspect the 5-ethynyl-20-deoxyuridine (EdU). H520 or H2170 (1 × 104) cells were each planted in 96-well dishes for insertion with si-circPVT1, vector-circPVT1, mock or negative control (NC) oligonucleotide via transfection. After incubation for 48 h at 5% CO2 and 37 °C, 50 mM EdU was added to the cells followed by another 2 h of incubation. The cells were fixed using paraformaldehyde (4%), stained using the Apollo Dye Solution for cell proliferation, whereas DAPI staining was performed for nucleic acids. The rate of cell proliferation was computed as outlined in manufacturer provided instructions. The images were visualized under a fluorescence microscope (Olympus FSX100).
Luciferase reporter assay
CircPVT1, miR-30d, and miR-30e sequences and their correspondent mutant forms without docking sites were constructed and sub-cloned into luciferase enzyme reporter vector, termed circPVT1-WT, circPVT1-Mut, miR-30d-WT, miR-30d –Mut, miR-30e-WT, miR-30e –Mut. Sequencing was applied to verify these plasmids. Moreover, the relative luciferase enzyme activity was inspected using the Dual-Luciferase Assay Kit (obtained from Promega, Madison, WI, USA) following the protocol provided by the manufacturer.
Further, 1 × 103 cells were plated onto 60 mm culture plates. Then, the cells were incubated at 37 °C and 5% CO2 for 10 days. Subsequently, colony staining was performed using crystal violet for 30 min (Sigma-Aldrich) and the cells were enumerated.
Western blot analysis
Cell homogenization was performed on ice using the lysis constituting Hepes (50 mM, pH 7.5), protease inhibitor cocktail (Roche), EDTA (5 mM, pH 8.0), NaCl (150 mM), β-glycerophosphate (20 mM), NP40 (0.5%), sodium orthovanadate (0.1 mM), dithiothreitol (DTT) (1 mM), NaF (50 mM), PMSF (1 mM) and MLUADl2 (10 mM). Clarification of the cell lysates was performed via 10 min centrifugation at 4 °C. Thereafter, 30μg/lane of the proteins were loaded into wells and fractionated using SDS-PAGE gels (10%) and transfer-embedded onto nitrocellulose membranes. Conjugation of the immunoblots with primary rabbit polyclonal anti-CCNF (Abcam, USA), and mouse monoclonal anti-GAPDH (Calbiochem) was conducted. The immunostained bands were developed, and the inspected via a chemiluminescent strategy.
RNA immunoprecipitation (RIP) evaluation
The Magna RIP kit (Millipore, Billerica, MA, USA) was used in the RIP assay following the manufacturer’s instructions. H520 cells were lysed using the complete RNA lysis buffer. Then, the cell lysates with magnetic beads probed anti-HuR (Abcam, USA) or negative control IgG antibody (Abcam, USA) were incubated at 4 °C for 4 h. The beads were rinsed using the washing buffer. Eventually, immunoprecipitated RNA was purified, followed by enrichment to evaluate the targeted RNAs via qRT-PCR.
CircPVT1 probe (5′-GCCAAAAGATCAGGCCTCAAGCCCAGCTGAGCGCCG.
GATGGAACG-3′-Biotin), and the control probe were synthesized (Jinkairui, Wuhan, China). H520 cells were lysed using the lysis buffer. Subsequently, the lysed cells were incubated with the respective circPVT1 probes. Then, the cell lysates were incubated with the streptavidin-coated magnetic beads to pull down the biotin-conjugated RNA complex. The beads were rinsed, and the RNA complex was purified using TRIzol (Takara, Dalian, China). Lastly, the immunoprecipitated RNA and protein were analyzed by miRNA sequencing and proteomic identification (LC-MS/MS).
The Xi’an Jiaotong University Animal Care and Use Committee approved all the mice experiments. Besides, the mice experiment protocols adhered to the guidelines of the National Institutes of Health. Stably down expressed cell lines were established through transfection of H520 cells and selection using puromycin. H520cells were infected using the lentiviruses (Jikai, Shanghai, China) comprising NC and si-cPVT1. Regarding the xenograft assays, diverse types of H520 cells (2 × 106) were subcutaneously administered into nude mice. The tumor volume and body weight were inspected once every week. The mice were sacrificed after 7 weeks, then the tumors were harvested for subsequent experimentation.
For the IHC assay, the paraffin-embedded segments were probed overnight at 4 °C with the primary antibodies against CCNF (1:100) (Abcam, CA). Then, the secondary antibodies were conjugated via incubation at 37 °C for 1 h. The HRP-conjugated streptavidin solution was incubated with the segments for 10 min, followed by staining using DAB.
Statistical analyses were performed using GraphPad Prism 8. Data were indicated as mean ± standard deviation (SD). The student’s t-test, or χ2 test was used to evaluate the difference between study groups. Moreover, the Kaplan-Meier approach was applied to establish the survival rates. Furthermore, the relationship linking study groups was determined through the Pearson correlation. Notably, P < 0.05 denoted statistical significance.
Differential circRNA expression in LUSC
To elucidate the expression patterns of circRNA and mRNA in LUSC, RNA-seq tools were used to sequence RNA from LUSC and normal lung tissues. Then, the DE-Seq algorithm was applied in filtering the differentially expressed genes. Besides, significant differences in expression were determined by P-value and FDR evaluation as follows: Fold Change > 2 or < 0.5, P-value < 0.05, FDR < 0.05.
Relative to the normal tissues, the expression analysis revealed that the circRNA transcripts were differentially expressed in cancerous tissues (Fig. 1a, b and c). Among the 12 differentially expressed circRNAs, 7 were upmodulated while 5 were downmodulated in LUSC tissues in contrast with the normal tissues. Among the 7 up-modulated circRNA, circPVT1 demonstrated the most significant effect of cell proliferation, hence, circPVT1 was selected as a research target.
Upregulation of CircPVT1 in LUSC samples correlates with the prognosis of LUSC patients
To establish whether the expression of circPVT1 was changed in LUSC, the expression of circPVT1 was assessed in 104 pairs of LUSC tissues and 110 pairs of LUSC serum using qRT-PCR. In contrast with the vicinal normal tissues, our data showed upregulated circPVT1 expressions in LUSC tissues (Fig. 1d) and serum (Fig. 1f). Similarly, the levels of circPVT1 expressions in lung cancer cell lines (H2170, H520, SKMES1, H226, 95c, 95d, A549, and PC9) were upregulated in contrast with the healthy human bronchial epithelial cell line (BEAS-2B) (Fig. 1i).
Additionally, the diagnostic significance of the expression of circPVT1 in tissues and serum of LUSC patients were evaluated. Consequently, for the circPVT1 expression in tissues, the ROC analysis exhibited an area under curve (AUC) of 0.774, and a diagnostic sensitivity and specificity of 97.1 and 51% (95% CI = 0.711–0.837), respectively (Fig. 1e). For serum circPVT1 expressions, the AUC was 0.789, and the diagnostic sensitivity and specificity were 91.3 and 60.6% (95% CI = 0.726–0.852), respectively (Fig. 1g). This indicated that serum circPVT1 potentially acts as a biomarker for the diagnosis of LUSC. To investigate the relationship between circPVT1 and patient prognosis, a cohort of 104 LUSC subjects with clinicopathological features and survival data was analyzed. As a result, the content of circPVT1 expression was directly linked to tumor node metastasis stage (TNM) (P = .0001), lymph node metastasis (P = .0014) and tumor size ((P = .0468), but not with other clinicopathological characteristics including age, or gender (Table 1). Multivariate analysis identified circPVT1 expression as an independent predictor of OS (Table 2). Besides, LUSC subjects with lower circPVT1 expression exhibited a better overall survival relative to those with elevated circPVT1 expression (Fig. 1h). Collectively, these data indicate that circPVT1 was overexpressed in LUSC tissues and serum, also, higher circPVT1 expression caused a poor prognosis of LUSC patients.
CircPVT1 promotes cell proliferation in vitro
Further, the functional role of circPVT1 in LUSC progression was investigated. First, the expression of circPVT1 in H520 and H2170 cell lines were silenced by transfecting with circPVT1 siRNA. To evaluate the role of circPVT1 in LUSC cells, 2 siRNAs were designed against circPVT1 and used to silence circPVT1 without influencing PVT1 level. The circular transcript expression vector circPVT1 was effectively constructed, as it could upregulate circPVT1 expression level rather than mRNA (Fig. 2a). CCK-8 (Fig. 2b) and EdU (Fig. 2c, d) assays data illustrated that si-circPVT1 could inhibited cell proliferation in H520 or H2170 cell lines, whereas over-expression of circPVT1 might promoted cell proliferation in H520 or H2170 cells lines. Further, colony formation evaluations revealed significantly improved cell cloning potentials of H520 or H2170 cells. This improvement was related to the upregulation of circPVT1 but was significantly impaired by the down-modulation of circPVT1 (Fig. 2e, f). These data suggested that circPVT1 promotes LUSC cell proliferation.
CircPVT1 is actively triggered by HuR in LUSC cells
Mounting research evidence confirmed that multiple key proteins bind to circRNAs and mediate RNA dysregulation in human cancer cells. To search for proteins interacting with the circPVT1, RNA pull-down experiments were conducted. Mass spectrometry discovered 112 proteins associated with circPVT1, including HuR (Fig. 3a, b). HuR binds and stabilize cancer-associated mRNAs. However, whether HuR affects circPVT1 remains unclarified. To confirm the potential interaction between HUR and circPVT1, HuR specific antibody was used to perform RNA precipitation. Through RIP and qRT–PCR, a direct interaction was found between HuR and circPVT1 (Fig. 3c). Further, 2 siRNAs were designed against HuR (Fig. 3d) and downregulation of HuR suppressed the expression of circPVT1 (Fig. 3e). The vector HuR was effectively constructed, as it could upregulate HuR expression level rather than mRNA (Fig. 3f) and HuR overexpression upregulated the expression of circPVT1 (Fig. 3g). Whereas, circPVT1 did not change the level of HuR (Fig. 3h). Therefore, it was inferred that the up-modulation of circPVT1 in LUSC might have been triggered by HuR.
CircPVT1 sponges the miR-30d/e
Given that circRNA has been documented as a miRNA sponge and circPVT1 is stably expressed in the cytoplasm, the plausibility of circPVT1 binding to the miRNAs was investigated. To identify the functional miRNAs potentially interacting with circPVT1 in LUSC cells, circPVT1 specific probe was used to perform RIP as described previously . Through RIP and RNA-seq (Fig. 4a), 244 downregulated RNA and 66 upregulated RNA were identified. Then, the potential miRNAs were screened which were predicted by starBase 3.0.
Among the candidate miRNAs predicted by starBase 3.0 and RNA-seq, 9 miRNAs were analyzed, comprising miR-140, miR-30a, miR542, miR-30b, miR-365a, miR-30c, miR-30e, miR-365b, and miR-30d (Fig. 4b). Then, a luciferase enzyme assay was conducted. A mimic of each miRNA was co-inserted with the luciferase reporters into H520 cells via co-transfection. As expected, relative to the control, miR-30d and miR-30e reduced the luciferase reporter activity in both H520 cells (Fig. 4c) and T293 cells (FigS1). MiR-30d and miR-30e belong to the miR-30 family and share similar miRNA seed sequences. Besides, both of these miRNAs have been documented as tumor repressors [7, 8]. Since only one potential docking site of the miR-30d/e was identified within the circPVT1 sequence, the possible circPVT1 target site was mutated with the inclusion of each miRNA sequence in the 3’UTR (Fig. 4d). We established that the transfection of H520 cells with miR-30d and miR-30e had an apparent distinct influence on luciferase enzyme activity when the correspondent target regions were mutated from the luciferase reporter (Fig. 4e). qRT-PCR analysis revealed that miR-30d or miR-30e did not show significant changes after circPVT1 was overexpressed or silenced (Fig. 4f) and miR-30d or miR-30e did not decrease the level of circPVT1 (Additional file 3, FigS2).
Functionally, CCK8 assay was used to verify whether miR30d/e could reverse the function of circPVT1 on cell proliferation. As showing in Fig. 4g, downregulation of circPVT1 inhibited cell proliferation, treatment with miR-30d and 30e inhibitors promoted cell proliferation, a combination of si-cPVT1 with miR-30d or miR-30e inhibitors partially reduced the effect of miR-30d/e inhibitors to promote cell proliferation and si-cPVT1 to inhibit cell proliferation. Si-cPVT1 combined with miR-30d inhibitors and miR-30e inhibitors significantly reduced the ability of miR30d or miR-30e inhibitors to promote cell proliferation and si-cPVT1to inhibit cell proliferation. The result revealed that downregulation of circPVT1 significantly suppressed the proliferation of H520cells, and the suppression could be blocked by miR-30d or miR-30e. Our results confirmed the hypothesis that circPVT1 regulates cell proliferation by combining with miR-30d and miR-30e.
CircPVT1 suppresses the proliferation of LUSC by sponging the miR-30d/e to regulate CCNF
The mRNA expression analysis of RNA-seq for elucidating the expression patterns of mRNA in LUSC revealed that 2441 mRNAs were differentially expressed. Out of the 2441 mRNAs, 835 were upregulated while 1606 mRNAs were down-regulated (Fig. 5a). To investigate whether circPVT1 exerts its effect by sponging the activity of the miR-30d/e, both Starbase 3.0 interaction prediction maps and upregulation shifts from RNA-seq were used to draw 105 RNAs based on the RNA binding with miR-30d and miR-30e predicted by starBase 3.0, upregulated RNA in LUSC. Protein-protein interaction (PPI) network was constructed using the STRING repository via the Cytoscape software. Consequently, a remarkable module with 16 nodes was established via the MCODE plugin. In total, 22 key genes were ranked by combining 6 methods as highlighted by cytoHubba (Fig. 5b). Out of the 20 genes (Additional file 4, FigS3), cyclin F (CCNF) and DEP domain-containing protein 1B (DEPDC1B) were consistently down-regulated after inhibition circPVT1 expression (Fig. 5c). Using qRT-PCR, the expression of these 2 genes after conditional knockdown of miR-30d or miR-30e was evaluated. Our results revealed that silencing of miR-30d or miR-30e improves the expression of CCNF in H520 cells (Fig. 5d).
Further, the luciferase reporter assays were used to evaluate whether miR-30d or miR-30e directly binds to CCNF 3′ UTR in H520 cells (Fig. 5e). The findings illustrated that overexpression of miR-30d or miR-30e reduced the luciferase activity of a wild-type 3’UTR of CCNF reporter, and not the mutant one (Fig. 5f, g). Using the cancer genome atlas (TCGA) repository, CCNF was inversely linked to the expression of miR-30d (Fig. 5h) or miR-30e (Fig. 5i) in LUSC. On the other hand, western blot illustrated that down-expression of circ-PVT1 significantly repressed the expression of CCNF, and the suppression could be blocked by miR-30d or (and) miR-30e (Fig. 5j). Collectively, these findings suggested that circPVT1 regulates CCNF via the sponge activity of miR-30d and miR-30e in LUSC.
CCNF promotes cell proliferation in vitro
So far, the role of CNNF in NSCLS remains unreported, as such, this work further investigated the functional role of CCNF in LUSC progression. Three siRNAs against CCNF and 1 vector CCNF were effectively constructed and their knockdown and overexpression efficiencies were verified by qRT–PCR and western blot (Fig. 6a). Data on CCK-8(Fig. 6b and FigS4) and EdU assays (Fig. 6c) illustrated that siCCNF inhibited cell proliferation, whereas over-expression of CCNF might promoted cell proliferation in H520 or H2170 cells lines. Further, colony formation evaluations illustrated that the cell cloning potentials of H520 or H2170 cells were significantly impaired by down-modulation of CCNF and was upregulated by the upregulation of CCNF (Fig. 6d). IHC findings revealed that CCNF was remarkably up modulated in LUSC tissues (Fig. 6e). This was also confirmed TCGA repository (Fig. 6f). Functionally, CCK8 assay revealed that downregulation of miR-30d and miR-30e remarkably promoted the proliferation of H520 cells, a promotion which could be blocked by CCNF (Fig. 6g). Our findings verified that CCNF promotes LUSC cell proliferation and partly reverses the effect of promoting cell proliferation of miR30d and miR30e.
CircPVT1 suppresses growth of xenograft tumor in vivo
To verify the above findings, H520 cells were subcutaneously injected into nude mice. As shown in Fig. 7a, the development of tumors incubated from si-circPVT1 stably transfected cells was significantly slower than those tumors incubated from si-NC-transfected cells. Three weeks after cell inoculation, the average tumor volume, not the body weight, in the si-circPVT1 group was considerably lower compared to the NC group (Fig. 7b, c). Besides, through qRT-qPCR and western blot analysis, the expression contents of circPVT1 (Fig. 7e) and CCNF (Fig. 7d.f) were reduced in the tumors from the si-circPVT1 group.
In recent years, numerous circRNAs have been documented in tissues and cell lines. Based on cell/tissue-distinct and development stage-distinct expression and molecular structure of the circRNAs, they regulate various biological processes. Besides, they are potential diagnostic biomarkers or treatment targets for cancer, relative to linear transcripts . Reports indicate that a few circRNAs act as oncogenes or tumor repressors in colon cancer , gastric cancer , breast cancer , liver cancer  and other types of cancer. Nonetheless, so far, only a few circRNAs have been completely characterized and the biological role of most of the circRNAs remains largely unclear.
Drugs targeting the ALK receptor tyrosine kinase (ALK), EGFR mutations, or ROS1 rearrangements have improved the outcomes of LUAD patients. This highlights the significance of molecular, genomic profiling and the subtyping of lung tumors [14,15,16]. However, in contrast with LUAD, only a few drug-gable genomic changes have been documented for LUSC, hence dismal LUSC patient outcomes. Therefore, understanding the molecular mechanisms underlying the development of LUSC might reveal the identification of novel therapeutic approaches. We utilized RNA-seq to acquire the expression patterns of circRNA and mRNA in LUSC tissues. Then, we uncovered circPVT1, which was upregulated in LUSC tissues and serum, and whose expression correlated with clinical-pathological data.
Previous studies have documented the role of circPVT1 in the progression of cancer and other diseases. For instance, Chen et al. reported that circPVT1 improves cell proliferation and acts as a prognostic index in gastric cancer . Verduci et al. documented an oncogenic function of circPVT1 in head and neck squamous cell carcinoma through a mutant p53/YAP/TEAD transcription-competent complex . Besides, Panda et al. opined that circPVT1 is senescence linked circular RNA, which is elevated in multiplying cells but diminished in senescent cells and sequesters let-7 to facilitate proliferation . Elsewhere, Li et al discovered that transcription factor c-Fos regulate the expression of circPVT1, which mediates the progression of NSCLC. The NSCLC progression is achieved through the modulation of the E2F2 signaling cascade by binding to a miR-125b . On the other hand, Qin et al revealed that circPVT1 acts as a competing endogenous RNA of miR-497 and indirectly modulates the expression of Bcl-2 . While LUAD and LUSC are distinct in their histological, molecular, and clinical presentation, all the published experiments were performed in adenocarcinoma cells. We hypothesized that circPVT1 might be acting via different mechanisms in squamous cell carcinoma. In line with previous reports, our study reveals that circPVT1 serves as an oncogene in LUSC by different mechanisms unlike in LUAD.
The oncogenic influence of HuR was primarily attributed to its docking and stabilization of cancer-correlated mRNAs in the cytoplasm. Yang et al. discovered that circ-HuR downregulates the expression of HuR in gastric cancer . CircPABPN1 blocks HuR binding to ATG16L1 mRNA and lowers ATG16L1 production . Lan et al. established that the downregulation of lncOCC-1 triggers an escalation of HuR in Caco-2 cells . Abdelmohsen et al. reported that masse circRNAs docks HuR in HeLa cells, CircPABPN1 being the most predominant. Further investigations illustrated that HuR did not impact CircPABPN1 richness, however, high contents of CircPABPN1 repressed HuR docking to PABPN1 mRNA . Previous studies revealed that ncRNA (circRNA or lncRNA) regulates the expression of HuR. However, no data was demonstrating the effect of HuR on ncRNA abundance. For the first time, this study showed that HuR regulates the level of circPVT1, but the specific mechanism needs further investigations.
To get a profound understanding of the function of circPVT1 in LUSC, circRIP, RNA-seq, and bioinformatics were conducted to establish the downstream effector pathway after circPVT1 induction. Dual-luciferase reporter proved that circPVT1 could directly interact with miR-30d and miR-30e. Previous reports have revealed that miR-30d was a tumor suppressor in prostate cancer, ovarian cancer among other types of cancer [26, 27]. Besides, it is a potential non-invasive diagnostic biomarker in cervical cancer, rectal cancer, and NSCLC [28,29,30]. In NSCLC, miR-30d represses tumor cell proliferation and motility by directly targeting CCNE2 and NFIB [31, 32]. On the other hand, miR-30e acted as an oncogene in hepatocellular carcinoma, breast cancer among other types of cancer [33, 34]. In NSCLC, it inhibited cell proliferation and infiltration by directly targeting SOX9, PTPN1, and USP22-mediated Sirt1/JAK/STAT3 signaling [8, 35, 36]. Additionally, miR-30e was presented to be a potential non-invasive diagnostic biomarker for NSCLC . Our data illustrated that circPVT1 acts as an oncogene by sponging miR-30d and miR-30e in LUSC. Further, we highlighted the value of the cross-talk between circPVT1 and miR-30d or miR-30e in tumorigenesis and development of LUSC.
Reports hypothesize that circRNA acts as a ceRNA and modulates the expression of a miRNA target gene. Through bioinformatics, we found that miR-30d and miR-30e might play a vital role by sponging key genes. Based on qRT-PCR and dual-luciferase reporter assays, we discovered that miR-30d or miR-30e directly target the 3′-untranslated regions of CCNF. Besides, down-expression of miR-30d or miR-30e caused the up-regulation of CCNF at mRNA and protein levels. CCNF is the founding member of the F-box family of proteins, the substrate recognition subunits of Skp1-Cul1-F-box protein (SCF) ubiquitin ligase complexes. CCNF degradation permits rapid aggregation of RRM2 in response to genotoxic stress . A secondary response that involves p53R2, a paralog of RRM2, whose transcription was induced by p53 . P53R2 is likely defective in approximately 50% of human cancers bearing p53 mutations, making CCNF closely related to the development of cancer. In malignant germ cell tumors, CCNF participated in the LIN28/let-7 pathway as an oncogene . Acute myeloid leukemia patients with long-term proliferation demonstrated altered expression of CCNF . Using the TCGA dataset, CCNF was significantly up modulated in LUSC tissues. So far, the role of CNNF in NSCLS remain unreported. As such, we conducted in vitro experiments to confirm the role of CCNF in the proliferation of LUSC cells, which was following a previous study on the role of CCNF in cancers and the CCNF expression in LUSC in TCGA. To confirm the interaction between circPVT1 and CCNF, down-regulation of circPVT1 downregulated CCNF at both the mRNA and protein expression levels. Moreover, these impacts could be partly alleviated by miR-30d or miR-30e inhibitors. The data support our hypothesis that circPVT1 acts as a ceRNA to improve CCNF-mediated proliferation by sponging miR-30d and miR-30e in LUSC.
In conclusion, our findings revealed that elevated expression of circPVT1 is a potential independent prognostic marker for LUSC. Moreover, we demonstrated that HuR regulates the level of circPVT1 and circPVT1 might sponge miR-30d and miR-30e to modulate CCNF expression, thereby causing tumorigenesis and development of LUSC. We further inferred that circPVT1 is a potential diagnostic marker and a potential therapy target for LUSC. Finally, the regulatory network involving HuR/circPVT1/miR-30d and miR-30e/CCNF cascade potentially provides a profound understanding of the mechanism of pathogenesis and progression of LUSC.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its additional information files.
Competing endogenous RNA
Lung squamous cell carcinoma
Real-time quantitative polymerase chain reaction
Small cell lung cancer
Non-small cell lung cancer
MicroRNA response elements
The Cancer Genome Atlas
DEP domain containing protein 1B
ALK receptor tyrosine kinase
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018;68(1):7–30. https://doi.org/10.3322/caac.21442.
Chen Z, Fillmore CM, Hammerman PS, Kim CF, Wong KK. Non-small-cell lung cancers: a heterogeneous set of diseases. Nat Rev Cancer. 2014;14(8):535–46. https://doi.org/10.1038/nrc3775.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424. https://doi.org/10.3322/caac.21492.
Jeck WR, Sorrentino JA, Wang K, Slevin MK, Burd CE, Liu J, et al. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA. 2013;19(2):141–57. https://doi.org/10.1261/rna.035667.112.
Suzuki H, Tsukahara T. A view of pre-mRNA splicing from RNase R resistant RNAs. Int J Mol Sci. 2014;15(6):9331–42. https://doi.org/10.3390/ijms15069331.
Su X, Wang H, Ge W, Yang M, Hou J, Chen T, et al. An in vivo method to identify microRNA targets not predicted by computation algorithms: p21 targeting by miR-92a in Cancer. Cancer Res. 2015;75(14):2875–85. https://doi.org/10.1158/0008-5472.CAN-14-2218.
Zhuang H, Wu F, Wei W, Dang Y, Yang B, Ma X, et al. Glycine decarboxylase induces autophagy and is downregulated by miRNA-30d-5p in hepatocellular carcinoma. Cell Death Dis. 2019;10(3):192. https://doi.org/10.1038/s41419-019-1446-z.
Zhuang L, Shou T, Li K, Gao CL, Duan LC, Fang LZ, et al. MicroRNA-30e-5p promotes cell growth by targeting PTPN13 and indicates poor survival and recurrence in lung adenocarcinoma. J Cell Mol Med. 2017;21(11):2852–62. https://doi.org/10.1111/jcmm.13198.
Li X, Yang L, Chen LL. The biogenesis, functions, and challenges of circular RNAs. Mol Cell. 2018;71(3):428–42. https://doi.org/10.1016/j.molcel.2018.06.034.
Ju HQ, Zhao Q, Wang F, et al. A circRNA signature predicts postoperative recurrence in stage II/III colon cancer. Embo Mol Med. 2019;11:e10168.
Tang W, Fu K, Sun H, Rong D, Wang H, Cao H. CircRNA microarray profiling identifies a novel circulating biomarker for detection of gastric cancer. Mol Cancer. 2018;17(1):137. https://doi.org/10.1186/s12943-018-0888-8.
Yang R, Xing L, Zheng X, Sun Y, Wang X, Chen J. The circRNA circAGFG1 acts as a sponge of miR-195-5p to promote triple-negative breast cancer progression through regulating CCNE1 expression. Mol Cancer. 2019;18(1):4. https://doi.org/10.1186/s12943-018-0933-7.
Liang WC, Wong CW, Liang PP, Shi M, Cao Y, Rao ST, et al. Translation of the circular RNA circbeta-catenin promotes liver cancer cell growth through activation of the Wnt pathway. Genome Biol. 2019;20(1):84. https://doi.org/10.1186/s13059-019-1685-4.
Hanna N, Johnson D, Temin S, Masters G. Systemic therapy for stage IV non-small-cell lung Cancer: American Society of Clinical Oncology clinical practice guideline update summary. J Oncol Pract. 2017;13(12):832–7. https://doi.org/10.1200/JOP.2017.026716.
The Cancer Genome Atlas Research Network. Comprehensive genomic characterization of squamous cell lung cancers. Nature. 2012;489(7417):519–25. https://doi.org/10.1038/nature11404.
Chen F, Zhang Y, Parra E, Rodriguez J, Behrens C, Akbani R, et al. Multiplatform-based molecular subtypes of non-small-cell lung cancer. Oncogene. 2017;36(10):1384–93. https://doi.org/10.1038/onc.2016.303.
Chen J, Li Y, Zheng Q, Bao C, He J, Chen B, et al. Circular RNA profile identifies circPVT1 as a proliferative factor and prognostic marker in gastric cancer. Cancer Lett. 2017;388:208–19. https://doi.org/10.1016/j.canlet.2016.12.006.
Verduci L, Ferraiuolo M, Sacconi A, Ganci F, Vitale J, Colombo T, et al. The oncogenic role of circPVT1 in head and neck squamous cell carcinoma is mediated through the mutant p53/YAP/TEAD transcription-competent complex. Genome Biol. 2017;18(1):237. https://doi.org/10.1186/s13059-017-1368-y.
Panda AC, Grammatikakis I, Kim KM, de S, Martindale JL, Munk R, et al. Identification of senescence-associated circular RNAs (SAC-RNAs) reveals senescence suppressor CircPVT1. Nucleic Acids Res. 2017;45(7):4021–35. https://doi.org/10.1093/nar/gkw1201.
Li X, Zhang Z, Jiang H, Li Q, Wang R, Pan H, et al. Circular RNA circPVT1 promotes proliferation and invasion through sponging miR-125b and activating E2F2 signaling in non-small cell lung Cancer. Cell Physiol Biochem. 2018;51(5):2324–40. https://doi.org/10.1159/000495876.
Qin S, Zhao Y, Lim G, Lin H, Zhang X, Zhang X. Circular RNA PVT1 acts as a competing endogenous RNA for miR-497 in promoting non-small cell lung cancer progression. Biomed Pharmacother. 2019;111:244–50. https://doi.org/10.1016/j.biopha.2018.12.007.
Yang F, Hu A, Li D, Wang J, Guo Y, Liu Y, et al. Circ-HuR suppresses HuR expression and gastric cancer progression by inhibiting CNBP transactivation. Mol Cancer. 2019;18(1):158. https://doi.org/10.1186/s12943-019-1094-z.
Li XX, Xiao L, Chung HK, Ma XX, Liu X, Song JL, et al. circPABPN1 Interaction between HuR and Modulates Autophagy in the Intestinal Epithelium by Altering ATG16L1 Translation. Mol Cell Biol. 2020;40(6):e00492–19. https://doi.org/10.1128/MCB.00492-19.
Lan Y, Xiao X, He Z, Luo Y, Wu C, Li L, et al. Long noncoding RNA OCC-1 suppresses cell growth through destabilizing HuR protein in colorectal cancer. Nucleic Acids Res. 2018;46(11):5809–21. https://doi.org/10.1093/nar/gky214.
Abdelmohsen K, Panda AC, Munk R, Grammatikakis I, Dudekula DB, de S, et al. Identification of HuR target circular RNAs uncovers suppression of PABPN1 translation by CircPABPN1. RNA Biol. 2017;14(3):361–9. https://doi.org/10.1080/15476286.2017.1279788.
Ye Z, Li J, Han X, Hou H, Chen H, Zheng X, et al. TET3 inhibits TGF-beta1-induced epithelial-mesenchymal transition by demethylating miR-30d precursor gene in ovarian cancer cells. J Exp Clin Cancer Res. 2016;35(1):72. https://doi.org/10.1186/s13046-016-0350-y.
Lin ZY, Chen G, Zhang YQ, He HC, Liang YX, Ye JH, et al. MicroRNA-30d promotes angiogenesis and tumor growth via MYPT1/c-JUN/VEGFA pathway and predicts aggressive outcome in prostate cancer. Mol Cancer. 2017;16(1):48. https://doi.org/10.1186/s12943-017-0615-x.
Zheng M, Hou L, Ma Y, Zhou L, Wang F, Cheng B, et al. Exosomal let-7d-3p and miR-30d-5p as diagnostic biomarkers for non-invasive screening of cervical cancer and its precursors. Mol Cancer. 2019;18(1):76. https://doi.org/10.1186/s12943-019-0999-x.
Bjornetro T, Redalen KR, Meltzer S, et al. An experimental strategy unveiling exosomal microRNAs 486-5p, 181a-5p and 30d-5p from hypoxic tumour cells as circulating indicators of high-risk rectal cancer. J Extracell Vesicles. 2019;8(1):1567219. https://doi.org/10.1080/20013078.2019.1567219.
Markou A, Sourvinou I, Vorkas PA, Yousef GM, Lianidou E. Clinical evaluation of microRNA expression profiling in non small cell lung cancer. Lung Cancer. 2013;81(3):388–96. https://doi.org/10.1016/j.lungcan.2013.05.007.
Chen D, Guo W, Qiu Z, Wang Q, Li Y, Liang L, et al. MicroRNA-30d-5p inhibits tumour cell proliferation and motility by directly targeting CCNE2 in non-small cell lung cancer. Cancer Lett. 2015;362(2):208–17. https://doi.org/10.1016/j.canlet.2015.03.041.
Wu Y, Zhang J, Hou S, Cheng Z, Yuan M. Non-small cell lung cancer: miR-30d suppresses tumor invasion and migration by directly targeting NFIB. Biotechnol Lett. 2017;39(12):1827–34. https://doi.org/10.1007/s10529-017-2428-9.
Deng L, Tang J, Yang H, Cheng C, Lu S, Jiang R, et al. MTA1 modulated by miR-30e contributes to epithelial-to-mesenchymal transition in hepatocellular carcinoma through an ErbB2-dependent pathway. Oncogene. 2017;36(28):3976–85. https://doi.org/10.1038/onc.2016.491.
Croset M, Pantano F, Kan C, et al. miRNA-30 family members inhibit breast Cancer invasion, Osteomimicry, and bone destruction by directly targeting multiple bone metastasis-associated genes. Cancer Res. 2018;78(18):5259–73. https://doi.org/10.1158/0008-5472.CAN-17-3058.
Cui Y, Zhao L, Zhao S, Guo T, Li F, Li Z, et al. MicroRNA-30e inhibits proliferation and invasion of non-small cell lung cancer via targeting SOX9. Hum Cell. 2019;32(3):326–33. https://doi.org/10.1007/s13577-018-0223-0.
Xu G, Cai J, Wang L, Jiang L, Huang J, Hu R, et al. MicroRNA-30e-5p suppresses non-small cell lung cancer tumorigenesis by regulating USP22-mediated Sirt1/JAK/STAT3 signaling. Exp Cell Res. 2018;362(2):268–78. https://doi.org/10.1016/j.yexcr.2017.11.027.
D'Angiolella V, Donato V, Forrester FM, Jeong YT, Pellacani C, Kudo Y, et al. Cyclin F-mediated degradation of ribonucleotide reductase M2 controls genome integrity and DNA repair. Cell. 2012;149(5):1023–34. https://doi.org/10.1016/j.cell.2012.03.043.
Tanaka H, Arakawa H, Yamaguchi T, Shiraishi K, Fukuda S, Matsui K, et al. A ribonucleotide reductase gene involved in a p53-dependent cell-cycle checkpoint for DNA damage. Nature. 2000;404(6773):42–9. https://doi.org/10.1038/35003506.
Murray MJ, Saini HK, Siegler CA, Hanning JE, Barker EM, van Dongen S, et al. LIN28 expression in malignant germ cell tumors downregulates let-7 and increases oncogene levels. Cancer Res. 2013;73(15):4872–84. https://doi.org/10.1158/0008-5472.CAN-12-2085.
Hatfield KJ, Reikvam H, Bruserud O. Identification of a subset of patients with acute myeloid leukemia characterized by long-term in vitro proliferation and altered cell cycle regulation of the leukemic cells. Expert Opin Ther Targets. 2014;18(11):1237–51. https://doi.org/10.1517/14728222.2014.957671.
We thank Dr. Chen Huang for her comments that improve the final version of the manuscript.
This work was supported by National Natural Science Foundation of China (No. 81672300), Shaanxi Provincial Natural Science Foundation (NO.2019SF-224) and the project of the free exploration and innovation project of Xi’an Jiaotong University (NO.xjj2018144).
Ethics approval and consent to participate
All tissue samples were in accordance with the ethical guidelines of the second affiliated hospital of Xi’an Jiaotong University. Written informed consents were obtained from all patients. The ethics approval statements for animal work were provided by Xi’an Jiaotong University Animal Care and Use Committee.
Consent for publication
All authors give consent for the publication of the manuscript in Journal of Experimental & Clinical Cancer Research.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Sequences of primers used in this study. Table S2. Sequences of siRNAs used in this study.
Additional file 2: Figure S1.
A luciferase reporter assay was used to detect the luciferase activity of LUC-cPVT1 in T293 cells transfected with miRNA mimics to identify miRNAs that bind to the circPVT1 sequence. Data are shown as mean ± SD, ***P < 0.001.
Additional file 3: Figure S2.
qRT-PCR analysis of circPVT1 expression in H520 cells after treatment with miR-30d minics or miR-30e minics or miR-30d inhibitors or miR-30e inhibitors.
Additional file 4: Figure S3.
qRT-PCR analysis of mRNA expression in H520 cells after treatment with si-cPVT1–1 and si-cPVT1–2.
Additional file 5: Figure S4.
Assessment of proliferation in H520 cells transfected with NC, si-CCNF-1, siCCNF-2, si-CCNF-3 by CCK-8 assay. *P < 0.05;**P < 0.01; ***P < 0.001.
Rights and permissions
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Shi, J., Lv, X., Zeng, L. et al. CircPVT1 promotes proliferation of lung squamous cell carcinoma by binding to miR-30d/e. J Exp Clin Cancer Res 40, 193 (2021). https://doi.org/10.1186/s13046-021-01976-w
- Lung squamous cell carcinoma