- Open Access
CircRNA 0009043 suppresses non-small-cell lung cancer development via targeting the miR-148a-3p/DNAJB4 axis
Biomarker Research volume 10, Article number: 61 (2022)
Circular RNAs (circRNAs) are important regulators of the development and progression of non-small-cell lung cancer (NSCLC) and many other malignancies. The functional importance of circ_0009043 in NSCLC, however, has yet to be established.
The expression of circ_0009043, miR-148a-3p, and DnaJ heat shock protein family (Hsp40) member B4 (DNAJB4) in NSCLC cells was assessed via qPCR. The proliferative activity of these cells was examined through EdU uptake and CCK-8 assays, while flow cytometry approaches were used to examine apoptotic cell death rates. Protein expression was measured through Western immunoblotting. Interactions between miR-148a-3p and circ_0009043 or DNAJB4 were detected through RNA immunoprecipitation (RIP) and dual-luciferase reporter assays. The in vivo importance of circ_0009043 as a regulator of oncogenic activity was assessed using murine xenograft models.
Both NSCLC cells and tissue samples were found to exhibit circ_0009043 upregulation, and lower circ_0009043 expression levels were found to be related to poorer NSCLC patient overall survival. Knocking down circ_0009043 resulted in the enhancement of NSCLC cell proliferative activity and the suppression of apoptotic tumor cell death in vitro, while also driving more rapid in vivo tumorigenesis. Mechanistically, circ_0009043 was found to function as a molecular sponge that sequestered miR-148a-3p, which was in turn able to directly suppress DNAJB4 expression. When miR-148a-3p was overexpressed, this reversed the impact of knocking down circ_0009043 on the apoptotic death and proliferation of NSCLC cells. Conversely, miR-148a-3p inhibition resulted in the suppression of NSCLC cell apoptosis and the enhancement of tumor cell growth, while the downregulation of DNAJB4 reversed these changes.
Circ_0009043 acts as a tumor suppressor in NSCLC cells, promoting DNAJB4 upregulation via the sequestration of miR-148a-3p.
Lung cancer remains the deadliest form of cancer in the world, accounting for approximately 1.6 million deaths annually . In 2020 alone, 228,820 new cases of lung cancer were reported in the USA together with 135,720 deaths . Non-small cell lung cancer (NSCLC) is the most prevalent lung cancer histological subtype, making up ~ 85% of these cases. The 5-year survival rate for NSCLC patients is poor at just 19%, in part as many patients are diagnosed when the disease is relatively advanced and late-stage treatments for this cancer are lacking [3, 4]. The diagnosis and treatment of NSCLC are challenging owing to the high degree of neoplastic diversity within tumors and the challenges inherent in collecting the volumes of tissue necessary for reliable pathological or molecular diagnostic analyses and characterization [5, 6]. Genetic changes in tumors can also shift the associated molecular marker profiles over the course of disease. Biopsy or cytology samples are the primary materials used to diagnose most patients with advanced NSCLC at present, as resected tumor tissue availability is limited . As such, there is a clear need for the identification of more reliable biomarkers of NSCLC and the further characterization of the mechanisms governing the pathogenesis of this disease so as to better guide the diagnosis and clinical treatment of patients such that the most effective surgical and therapeutic interventions can be administered to improve survival outcomes.
Circular RNAs (circRNAs) are a recently identified group of covalently closed non-coding RNAs formed through reverse splicing of their 3′ and 5′ ends . Owing to their structural characteristics, circRNAs are very stable, in addition to being abundantly expressed and highly conserved. Indeed, circRNAs, which can be derived from intergenic, exonic, or intronic regions of the genome, often exhibit tissue- and developmental stage-specific patterns of expression [8, 9]. Functionally, these circRNAs can serve as scaffolds for protein complexes such that they can regulate the expression of genes including the parental genes from which they are derived. Moreover, they can influence alternative splicing events, modulate interactions between RNAs and proteins, and function as molecular sponges to sequester microRNAs (miRNAs) in a sequence-specific manner . A growing body of evidence suggests that circRNAs are critical regulators of many forms of cancer including thyroid cancer, osteosarcoma, breast cancer, pancreatic cancer, hepatocellular carcinoma, and gastric cancer wherein they can act in an oncogenic manner, offering value as promising clinical biomarkers and/or viable targets for therapeutic intervention [11,12,13,14,15]. The recently identified circRNA_0009043 species of circRNA has been reported to be significantly downregulated in grade 1–3 endometrial cancer tissues relative to normal paracancerous endometrial tissue . The functional importance of this circRNA in NSCLC, however, has yet to be defined.
The present study was developed to explore the relationship between circRNA_0009043 and NSCLC pathogenesis in vitro and in vivo. Functionally, this circRNA was found to inhibit NSCLC cell proliferation while promoting the apoptotic death of these cells. Mechanistically, circRNA_0009043 was found to regulate the miR-148a-3p / DNAJB4 axis so as to suppress NSCLC development.
All specimens were obtained from the Department of Thoracic Surgery, Hunan Provincial People’s Hospital. A total of 132 pairs NSCLC tissues and compared normal tissues were confirmed pathologically. More specifically, tumor and adjacent normal tissues were obtained from the patients who were diagnosed with NSCLC and underwent initial surgery between September 2014 and June 2016. These patients had not received any systemic treatment, including radio-chemical therapy, before sampling. The protocol was approved by the ethical review committee of Hunan Provincial People’s Hospital, and patient consent was obtained before the samples were taken.
Cell line culture and transfection
The human NSCLC cell lines (A549 cells and HCC827 cells) were provided by the American Type Culture Collection (ATCC, Rockville, MD). Cells received incubation treatment in RPMI-1640 (Solarbio, Beijing, China) and dulbecco’s modified eagle medium (DMEM; Solarbio) composed of 1% penicillin/streptomycin (Solarbio) and 10% fetal bovine serum (FBS; Solarbio) under 5% CO2 at 37 °C. A pcircRNA2.2-circ_0009043 overexpression plasmid and empty pcircRNA2.2 plasmid, the pcDNA 3.1 plasmid specific to DNAJB4 (termed OE- DNAJB4) and empty pcDNA3.1 plasmid were purchased from Guangzhou RiBoBio Co., Ltd. (Guangzhou, China). Circ_0009034 short hairpin (sh) RNA (shcirc) (AATTCAAAAAACCATGAAGCAAAATCAAGTGATCTCTTGAAGTAGCACAACATTCTCCACCCG), scrambled shRNA (shCTRL), DNAJB4 short hairpin (sh) RNA (shDNAJB4) (AATTCAAAAAACAGAAGCTTATGAAGTATTGAGTTCTCTTGAAGTAGCACAACATTCTCCACCCG), miR-148a-3p mimics, NC miRNA mimic (mimic NC), miR-148a-3p inhibitor, and inhibitor NC were designed and synthesized by GenePharma Co., Ltd. (Shanghai, China). Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) was employed to transfect abovementioned plasmids into A549 cells and HCC827 cells. 48 h after the transfection, we harvested cells for later study.
Cell proliferation experiment
The proliferation of HCC cells was tested through Cell Counting Kit-8 (CCK-8) assay, the transfected A549 and HCC827 cells (1 × 104 cells/well) were added into 96-well plates and maintained for 48 h. Next, 10 μL CCK-8 reagent (Sigma-Aldrich) was added into each well and cultured for further 2 h. After that, a microplate reader was adopted to examine the optical density value at 450 nm.
Cell proliferation was also assessed using an EDU assay kit (Ribobio, Guangzhou, China). For this purpose, after seeding 1 × 106 cells into each well of confocal plates, a 2-h incubation was performed at 37 °C with 50 μM of EDU buffer. Cell fixing was then carried out for half an hour using 4% of formaldehyde before cell permeabilization for 20 min using 0.1% Triton X-100. Eventually, after adding the EDU solution, cell nuclei were stained with Hoechst prior to visualization under a fluorescence microscope.
In vivo xenograft experiments
Nude female mice from Vital River Laboratories (Beijing, China) were housed under specific pathogen-free conditions in the Department of Laboratory Animal Science of Hunan Provincial People’s Hospital. The Ethics Review Committee of the Department of Laboratory Animal Science of Hunan Provincial People’s Hospital approved this study. All animals were randomly assigned to five groups (n = 6/group) that were then subcutaneously implanted in the left flank with A549 cells (1 × 107) that were either WT cells or cells that had been stably transfected with the pcircRNA2.2, pcircRNA2.2-circ_0009043 vector, NC‑shRNA, or circ_0009043‑shRNA constructs. Tumor volume was measured once per week using calipers (V = length × width2 × 0.5) for four weeks, after which mice were euthanized. Tumors were then excised, fixed for 24 h using 4% paraformaldehyde, stained via TUNEL assay, and subjected to downstream immunohistochemical staining for DNAJB4, followed by imaging with a brightfield microscope (Olympus, Tokyo, Japan).
Quantitative real-time reverse transcription PCR
Tissue samples and cells were collected, and total RNA was extracted using the TRIzol reagent (Invitrogen, Waltham, MA, USA) for quantitative real-time reverse transcription PCR (qRT-PCR) detection. The cDNA obtained by reverse transcription using the Promega GoScript reverse transcription system (Promega Madison, WI, USA) was subjected to quantitative real-time PCR in a 20 µL PCR amplification system using the ABI7500 real-time quantitative PCR instrument. For miRNAs, miRNA first-strand cDNA was constructed using the stem-loop method  using the cDNA Synthesis Kit (R601, Novabio, Shanghai, China). Stem-loop primer 5’-CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGACAAAGTT-3’ was purchased from General Biol (Anhui, China). The differences in the expression levels of different genes were analysed using 2−ΔΔCt relative quantification. The primers were circ_0009043: (F: 5’-TCCGCAAACATTCAGACAAA-3’ and R: 5’-GCTTCAGCTCTTCCATTGCT-3’); miR-148a-3p: (F: 5’-ACACTCCAGCTGGGTCAGTGCACTACAGAACT-3’ and R: 5’-CTCAACTGGTGTCGTGGA-3’); DNAJB4: (F: 5’-TGGGGACGCTGTTTTCTTTTAC-3’ and R: 5’-CTCCTGCTCCTCCTTTCAACC-3’); U6: (F: 5’-CTCGCTTCGGCAGCACA-3’ and R: 5’-AACGCTTCACGAATTTGCGT-3’); and GAPDH: (F: 5’-CTCCTCCTGTTCGACAGTCAGC-3’ and R: 5’-CCCAATACGACCAAATCCGTT-3’).
After isolation from tissue or cell samples, total proteins were separated via SDS-PAGE and transferred onto PVDF membranes that were blocked with an appropriate solution for 1 h at room temperature. Blots were then probed for 2 h at room temperature with primary antibodies specific for Bax (1:1000, BA0315-2, Boster), Bcl-2 (1:1200, A00040-1, Boster), Cytochrome C (1:5000, ab133504, Abcam), DNAJB4 (0.1 µg/ml, ab254641, Abcam), and GAPDH (1:10,000, BA2913, Boster). Blots were then probed for 1 h at room temperature using Goat Anti-Rabbit IgG H&L (HRP) (1: 2000, ab6721, Abcam), after which they were washed, incubated in developer solution, and exposed to X-ray film to detect protein bands.
Nuclear and cytoplasmic fractions from A549 and HCC827 cells were isolated using a PARIS Kit (Invitrogen) based on provided directions. Circ_0009043 abundance in these two fractions was assessed via qPCR, with U6 and GAPDH being used as respective normalization controls for nuclear and cytoplasmic transcripts.
Dual-luciferase reporter gene assay
Potential miRNA binding partners for circ_0009043 were predicted using an online tool (https://circinteractome.nia.nih.gov/). Luciferase reporter constructs were then prepared for circ_0009043 and miRNAs of interest, with HEK293T cells being transfected using WT or mutant circ_0009043 reporter constructs and transfected with appropriate miRNA mimic or control constructs. At 48 h post-transfection, firefly and Renilla luciferase activity was assessed, with results being reported as the ratio of these two activity levels.
RNA immunoprecipitation (RIP)
A Magna RIP RNA‐Binding Protein Immunoprecipitation Kit (Millipore, MA, USA) was used based on provided directions. Total RNA and species-appropriate IgG (n = 3) were utilized as control samples to confirm the specificity of the detected RNA signal. For the anti‐AGO2 RIP assay, NSCLC cells were transfected using miRNA mimics after which RIP was performed using anti-AGO2 (2 µg/mL, clone number: ab32381; Abcam) at 48 h post-transfection.
After incubation for 48 h, cells that have been transfected and were seeded for the counting in the 96-well plates were then removed from the plates by using trypsin reagent in the absence of ethylenediaminetetraacetic acid, and cells were then rinsed with cooled PBS solution. After that, centrifugation was done, and the PBS solution was discarded. Then annexin V–fluorescein isothiocyanate (FITC) apoptosis discovery kit (BioLegend, Inc.) was applied to detect the status apoptotic of the obtained cells. Resuspension of cells was then carried out in 100 μL of buffer (1 × binding buffer). Then resuspended cells were further treated with 5 μL propidium iodide and annexin V–FITC in the darkness for the 15-min duration, and at room temperature. Finally, for the analysis of the double-stained cells, the FACSCalibur Flow Cytometry instrument (BD Biosciences, San Jose, CA, USA) was used. Version 2.9 of the CellQuest (BD Biosciences) was used for the evaluation of the proportion of apoptotic cells.
Immunohistochemical (IHC) staining
Tissue sections were deparaffinized, subjected to antigen retrieval, treated to quench endogenous peroxidase activity, blocked using goat serum (Gibco), and incubated with primary anti-DNAJB4 (1 µg/ml, ab254641, Abcam) and secondary antibodies (abcam; #ab205718; 1:1000). Sections were then imaged using a phase-contrast microscope (Leica, Cat. #DMI 1).
Terminal deoxynucleotidyl transferase-medimed dUTP nick-endlabeling (TUNEL) assay
A one-step TUNEL kit (C1089, Beyotime Institute) was used as in prior reports to assess apoptotic cell death. Briefly, NSCLC tissue Sect. (3 μm) were deparaffinized, treated for 30 min with DNase- protease K (20 μg/mL) at 37℃, rinsed using PBS, and stained in the dark for 1 h with 50 μL of TUNEL reaction mixture at 37℃. Sections were then rinsed using PBS and imaged with a fluorescence confocal microscope (Zeiss LSM710, Germany). TUNEL-positive cells in 10 random fields of view in size sections were then counted using Image J (Bio-Rad Laboratories, CA, USA).
SPSS 18.0 (SPSS, Inc.) was used to analyze all data. Data were compared between groups using t-tests and one-way ANOVAs with Dunnett’s post hoc test as appropriate. Pearson’s correlation coefficient was used to assess the relationships between circ_0009043, miR-148a-3p, and DNAJB4 expression. Data are given as means ± SD from triplicate experiments, with P < 0.05 as the significance threshold.
Circ_0009043 is downregulated in NSCLC and predicted a poor prognosis
As shown in Fig. 1A, circular properties of the circ_0009043 were recognized through divergent primers as well as convergent primers. Then, we assessed its expression in tissue sections from NSCLC patients by FISH. Intriguingly, circ_0009043 was downregulated in the NSCLC as compared to the tissues of the normal lung (Fig. 1B). Moreover, the expression of circ_0009043 in the NSCLC tissues and in paired normal tissues was examined applying qRT-PCR approach. The results displayed that circ_0009043 level was downregulated significantly in the lung cancer tissues than in the adjacent nontumor tissues (Fig. 1C). Then, the 132 NSCLC tissues were separated into two sub-groups based upon circ_0009043 expression (the median value as cut-off). Kaplan–Meier curve was plotted for the analysis of the survival rate. The results indicated that circ_0009043 lower group of ovarian patients have a shorter survival rate (Fig. 1D). Moreover, the clinical relevance of this was then evaluated, showing that LINC00326 levels were linked with TNM stage, tumor differentiation, and lymphatic metastasis in NSCLC patients (Table 1). These results specified that circ_0009043 is clearly downregulated in NSCLC and predicted a poor prognosis. Afterwards, the cytoplasm and nucleus segmentation and RNA-FISH analyses verified that the circ_0009043 was predominantly localized in cell cytoplasm instead of their nucleus. RT-qPCR was performed to analyze fractions of the cytoplasmic as well as nuclear RNA in order to ascertain the circ_0009043 localization in the NSCLC cells. The obtained findings displayed that the circ_0009043 was predominantly detected within the cytoplasm of NSCLC cells (Fig. 1E), indicating that circ_0009043 had a primary effect on the cytoplasm of NSCLC cells.
Overexpression of circ_0009043 inhibits the proliferation, while accelerating apoptosis of the NSCLC cells
First, we found circ_0009043 was expressed highest in HCC827 cell line and lowest A549 cell line in different NSCLC cell lines, including A549, HCC827, HCI-H1299, NCI-H1568, NCI-H1650, and HCI-H358 (Data not shown). Thus, we selected HCC827 cell line and A549 cell line for further studies. The transfection efficacies of the circ_0009043 overexpression vector and shRNA targeting circ_0009043 were validated through the qRT-PCR in the A549 and HCC827 cells (Fig. 2A). After this, the CCK-8 assay was performed for determining the effect of circ_0009043 on cell proliferation. Figure 2B represents that the cell proliferation was remarkably inhibited in the circ_0009043 overexpression group than that of the negative control group (OE-CTRL) in the A549 as well as HCC827 cells. In the same way, circ_0009043 knockdown cells get proliferated faster as compared to the shCTRL group. In parallel, as presented in Fig. 2C, the Edu assay showed significantly less positive cells after circ_0009043 overexpression in A549 and in HCC827 cells. These obtained results indicate that circ_0009043 suppressed the NSCLC cells proliferation. Apoptosis is a key contributing factor that regulates the malignancy in the NSCLC. For investigating the cancer suppressor activity of circ_0009043 through apoptosis regulation, we scrutinized the apoptosis through FCM assay. The flow cytometry (FC) results displayed that the rate of apoptosis was decreased in the shCirc_0009043 transfected cells, and increased in the circ_0009043 overexpression cells as observed in Fig. 2D-E. The data obtained revealed that circ_0009043 accelerates the in vitro apoptosis of NSCLC cells. Additionally, the expression levels of Bax and Cytochrome C were gets decreased in the A549 as well as HCC827 cells transfected with sh-circ_0009043 as compared to sh-CTRL groups. On the contrary, circ_0009043 upregulation led to an increased level of the Bcl-2 as compared to the control groups (Fig. 2F-G). Collectively, circ_0009043 exerted a tumour-suppressing role in NSCLC.
Circ_0009043 functioned as a sponge for miR-148a-3p in NSCLC cells
We then identified the miR-148a-3p expression in the NSCLC as compared with the tissues of normal lungs and observed that the miR-148a-3p was upregulated considerably in NSCLC tissues (Fig. 3A). Similarly, an inverse expression correlation was observed between circ_0009043 and miR-148a-3p in NSCLC tissues (Fig. 3B). Additionally, we observed that miR-148a-3p might target circ_0009043 by using the online prediction software Starbase (https://starbase.sysu.edu.cn/). The active binding site for the miR-148a-3p in the 3’UTR of circ_0009043 was recognized (Fig. 3C). Furthermore, the transfection of circ_0009043 overexpression vector downregulated miR-148a-3p level in the A549 as well as HCC827 cells. On the contrary, the transfection of shRNA targeting circ_0009043 upregulated miR-148a-3p level (Fig. 3D). RIP experiments displayed that the miR-148a-3p coprecipitated with circ_0009043, confirming the strong physical interaction between the circ_0009043 and miR-148a-3p (Fig. 3E-F). Moreover, Luciferase reporter assay represented that the miR-148a-3p overexpression suppressed the WT- circ_0009043 reporter activity (Fig. 3G). Therefore, these findings demonstrated that circ_0009043 establishes a direct connection with the miR-148a-3p and consequently, it mirrors the role of a molecular sponge for that of the miR-148a-3p in the NSCLC cells.
MiR-148a-3p promotes proliferation, while inhibits apoptosis in the A549 and HCC827 cells
We further evaluate the miR-148a-3p role in the A549 and HCC827 cells. The transfection efficacies of the miR-148a-3p inhibitor and miR-148a-3p mimic were examined in HCC827 and A549 cells as represented in Fig. 4A. Then, a CCK-8 assay was conducted for determining the influence of miR-148a-3p on the proliferation of the cells. The proliferation of cells was enhanced significantly in the miR-148a-3p overexpression group by comparing to the proliferation in the negative control group in the A549 and HCC827 cells, as observed in Fig. 4B. Likewise, cells transfected with the miR-148a-3p inhibitor proliferated slower comparatively than those transfected with inhibitor NC cells. Consistently, the EdU assay revealed that upregulation of the miR-148a-3p enhanced the NSCLC cells proliferation as observed in Fig. 4C. The obtained results clearly demonstrated that the miR-148a-3p promoted the NSCLC cells proliferation. Afterwards, we examined apoptosis via FCM assay. By comparing with the NC group, the results of FC displayed that the rate of apoptosis decreased in miR-148a-3p mimics transfected cells and increased in the miR-148a-3p downregulated cells (Fig. 4D-E). These data showed that miR-148a-3p prevented apoptosis of the NSCLC cells in vitro. At last, we also carried out the western blotting for investigating the changes in the expression of Bax, Bcl-2, and also Cytochrome C. The upregulated levels of protein of BAX and Cytochrome C induced by miR-148a-3p knockdown demonstrated substantial downregulation upon transfection with miR-148a-3p mimics (Fig. 4F-G). These results are further confirmed that miR-148a-3p acts as an oncogene in NSCLC.
MiR-148a-3p Reversed the Regulatory Effect of circ_0009043 on A549 and HCC827 cells
To clarify further the interaction of the circ_0009043 and miR-148a-3p, we performed carefully a series of rescue experiments. The upregulated level of the miR-148a-3p in A549 and HCC827 cells transfected shRNA targeting circ_0009043 was reversed through miR-148a-3p downregulation (Fig. 5A). Downregulation of circ_0009043 promoted the viability in A549 and HCC827 cells that was partially reversed via co-transfection of the miR-148a-3p inhibitor (Fig. 5B). Meanwhile, EdU assay showed that circ_0009043 downregulation mediated acceleration effect on the proliferation capability of the A549 and HCC827 cells get weakened after miR-148a-3p downregulated (Fig. 5C). Identically, the reduced apoptosis cell numbers in the A549, as well as HCC827 cells due to the downregulation of circ_0009043, were elevated after the co-transfection of the miR-148a-3p inhibitor (Fig. 5D). In addition, miR-148a-3p downregulation attenuated the shCirc facilitated upregulation of the E-cadherin, while downregulation of the N-cadherin, vimentin, and slug expression (Fig. 5E). Thus, it is believed that the downregulation of miR-148a-3p leads to reverse the regulatory effects of the circ_0009043 knockdown on proliferative and apoptosis abilities of NSCLC cells.
DNAJB4 is a direct target gene of miR-148a-3p
For elucidating the fundamental mechanism of circ_0009043/miR-148a-3p axis in regulating the progression of NSCLC cells, we predicted the potential target genes for the miR-148a-3p through starBase. Consequently, we investigated the DNAJB4 as a potential target gene for the miR-148a-3p (Fig. 6A). Then, 66 tissues pairs of the clinical NSCLC and the matched adjacent were carefully chosen for detecting the expression levels of the circ_0009043 mRNA applying RT-PCR. Circ_0009043 mRNA expression in the NSCLC tissues was considerably downregulated than that observed in adjacent tissues (P < 0.0001; Fig. 6B), also recognized through IHC staining (Fig. 6C). Additionally, DNAJB4 protein in the NSCLC tissues was inspected using WB and it was detected that DNAJB4 expression was usually lower in the NSCLC tumor tissue than observed in adjacent tissue (Fig. 6D). The RT-qPCR assay was also carried out for detecting the expression level of DNAJB4 in NSCLC tissues as well as adjacent nontumor tissues (Normal). The obtained results demonstrated that the DNAJB4 level was noticeably decreased in NSCLC tissues as compared to Normal (Fig. 6E).
Similarly, an inverse expression correlation between DNAJB4 and the miR-148a-3p was notably detected, while DNAJB4 expression is a positive correlation with circ_0009043 (Fig. 6F-G). Then, we transfected circ_0009043 into A549 and HCC827 cells and clearly observed that circ_0009043 overexpression particularly enhanced the mRNA level of DNAJB4 in the A549 and HCC827 cells (Fig. 6H). Furthermore, miR-148a-3p mimic transfection significantly suppressed the level of DNAJB4 in the A549 as well as HCC827 cells (Fig. 6I), supporting miR-148a-3p targets DNAJB4. Also, we observed that the miR-148a-3p possibly targets DNAJB4 through bioinformatics investigation. Similarly, we recognized the active binding sites for the miR-148a-3p in 3’UTR of DNAJB4 (Fig. 6J). The Luciferase reporter assay specified that the miR-148a-3p overexpression suppressed the activity of the WT-DNAJB4 reporter (Fig. 6K). On the basis of the above results, it was speculated that DNAJB4 acts as a targeted gene of the miR-148a-3p.
Overexpression of DNAJB4 inhibits the proliferation but promotes apoptosis of NSCLC cells
In the next step, we evaluate the DNAJB4 role in the NSCLC cells. RT-qPCR, as well as western blotting, revealed upregulation in the expression of DNAJB4 in NSCLC cells transfected with DNAJB4 overexpression vectors. In addition to that, sh DNAJB4 led to a decreased level of DNAJB4 as also confirmed through qRT-PCR as well as western blotting (Fig. 7A, B). Then, CCK8 and EdU assays were carried out for determining the effect of DNAJB4 on the proliferation of the cell. Figure 7C represents that the proliferation of cells was considerably suppressed in the DNAJB4 overexpression group as compared to observed in the negative control group (OE-CTRL) in A549 as well as HCC827 cells. Likewise, DNAJB4 -Knockdown proliferation of cells is faster as compared with the shCTRL group. Meanwhile, Edu positive numbers of A549, as well as HCC827 cells transfected with DNAJB4 overexpression vector (OE- DNAJB4), were considerably lower than those that are transfected with OE-CTRL. An opposite process was noticed in the sh DNAJB4 group (shDNAJB4) from the process observed in the shCTRL group in A549 as well as HCC827 cells (Fig. 7D). Moreover, by comparing with the OE-CTRL group, the FC results demonstrate that the rate of apoptosis gets increased in the cells transfected with DNAJB4 overexpression vectors while decreased in the DNAJB4 downregulated cells (Fig. 7E). In the last, we carried out the western blotting for determining the changes that occur in the expression of Bcl-2, Bax, and Cytochrome C. These findings showed that DNAJB4 downregulation led to the upregulation of Bcl-2, while downregulation of Bax, and Cytochrome C expression. However, opposite results can be observed by transfection of overexpression DNAJB4 vector in A549 and HCC827 cells (Fig. 7F). These obtained results collectively provided evidence that DNAJB4 inhibits the proliferation but promotes apoptosis of NSCLC cells.
DNAJB4 Reversed the Regulatory Effect of miR-148a-3p on A549 and HCC827 cells
For exposing further the contact between circ_0009043, miR-148a-3p and DNAJB4, a series of rescue experiments were conducted. The transfection of the miR-148a-3p inhibitor reduced the viability in the A549 and NCI-HCC827 cells which were partially reversed by the co-transfection of shDNAJB4 (Fig. 8A). Meanwhile, the EdU assay displayed that the miR-148a-3p inhibitor facilitated the inhibition effect on the cell proliferation of the A549 and HCC827 cells were weakened after DNAJB4 downregulated (Fig. 8B). Identically, the increased apoptosis cell numbers in the A549, as well as HCC827 cells due to the transfection of the miR-148a-3p inhibitor, were suppressed after the knockdown of DNAJB4 (Fig. 8C). Moreover, DNAJB4 downregulation attenuated the miR-148a-3p inhibitor mediated upregulation of Bax, while downregulation of Bcl-2 expression (Fig. 8D). Thus, it is believed that the downregulation of DNAJB4 reversed the regulatory effects of the miR-148a-3p inhibitor on the proliferative and apoptosis abilities of NSCLC cells.
Circ_0009043 inhibits tumor growth via targeting the miR-148a-3p / DNAJB4 pathway in vivo
We applied finally the xenograft tumor model for investigating the influence of circ_0009043 on tumor growth in vivo. pcDNA4.0 vector, pcDNA4.0- circ_0009043 vector, shCTRL, or shCirc transfected A549 or the HCC827 cells were injected subcutaneously in the left flank of the athymic nude mice, respectively. Like the previous cell experiments results, the obtained result suggested that the tumors of the circ_0009043 knockdown group (shCirc) were considerably bigger than that of the shCTRL group (P < 0.0001), while the tumors of circ_0009043 overexpression group (OE-Circ) were considerably smaller from that of the negative control group (OE-CTRL, P < 0.01) as depicted in Fig. 9A-D. Meanwhile, circ_0009043 mRNA expression in the 5 tumor groups established through differently treated cells was investigated through RT-PCR. These results demonstrated that the expression of the circ_0009043 was significantly enhanced by OE-Circ transfection while downregulated with shCirc transfection (Fig. 9E). However, the transfection of circ_0009043 overexpression vectors downregulated miR-148a-3p level in tumor the tissues. On the contrary, the transfection of shCirc_0009043 upregulated miR-148a-3p level (Fig. 9F). Furthermore, WB assay and IHC staining of the xenograft tissues obtained from the nude mice confirmed that the knockdown of circ_0009043 can effectively prevent the expression of DNAJB4, while circ_0009043 upregulation enhanced the expression of DNAJB4 expression (Fig. 9G-H). Moreover, TUNEL assay represented that circ_0009043 overexpression (OE-Circ) resulted in much more apoptosis of cells than that in the control (OE-CTRL). However, the NSCLC cell’s apoptosis considerably decreased in the circ_0009043 knockdown group (shCirc), as compared to the shCTRL group (Fig. 9I-J). Western blotting was also carried out for evaluating the changes that occur in Bax, Bcl-2, and Cytochrome C expression. These results are almost consistent with the previously mentioned experiment on the cells (Fig. 9K). The overall results demonstrated that overexpression of circ_0009043 could suppress the growth as well as differentiation of NSCLC. We further confirmed that circRNA 0,009,043 targets the miR-148a-3p / DNAJB4 pathway to suppress the NSCLC in vivo.
Although initially disregarded as unimportant byproducts of aberrant splicing events, the advent of high-throughput sequencing and bioinformatics approaches has revealed that circRNAs are widely expressed in a diverse array of cell types and tissues in which they play important functional roles . Differential circRNA expression has been reported in many cancers, and particular circRNA species have been shown to influence cellular invasivity, metastatic progression, proliferation, and differentiation, in addition to being linked to disease status and patient prognostic outcomes [8, 19, 20]. In cervical cancer, for example, circRNA_0000285 and hsa_circ_0000745 have been shown to offer value as prognostic biomarkers [21, 22], while circ_0081001 may offer diagnostic value in osteosarcoma . Here, circ_0009043 was found to be downregulated in NSCLC tissues and cells, thus suggesting it to be an important regulator of NSCLC development.
In prior reports, circ_0009043 was found to be upregulated in grade 3 noncancerous endometrial tissues relative to grade 1–2 noncancerous endometrial tissues . The present study was the first to report on the functional importance of circ_0009043 in NSCLC. Specifically, circ_0009043 downregulation was evident in NSCLC and such downregulation was predictive of poorer NSCLC patient overall survival. Functionally, circ_0009043 knockdown augmented the proliferation of NSCLC cells, and in vivo xenograft assays confirmed that the upregulation of circ_0009043 was sufficient to suppress NSCLC tumor growth. Together, these results thus support a role for circ_0009043 as a tumor suppressor circRNA that interferes with NSCLC development.
Despite their lack of coding potential, miRNAs can prevent translation or induce mRNA degradation through sequence-specific binding to the 3’-untranslated regions of target mRNAs [24, 25]. Specific miRNAs have been shown to regulate the pathogenesis of many diseases and to influence normal tissue homeostasis [26, 27]. Functionally, miRNAs play a central role in tumor cell proliferation, metastasis, development, and invasion [28, 29]. Many circRNAs contain multiple miRNA binding sites such that they can serve as de facto miRNA sponges, sequestering these non-coding RNAs and thereby influencing cancer development [8, 10, 30]. In gastric cancer, for example, circ-PVT1 is upregulated and sequesters miR-124-3p, thereby upregulating ZEB1 expression and promoting tumor cell resistance to paclitaxel treatment . Circ-VANGL1 overexpression in bladder cancer can also drive oncogenesis via the miR-605-3p-VANGL1 axis . Notable roles for circRNAs as competing endogenous RNAs have also been reported in breast, ovarian, and colon cancers [33,34,35]. Potential miRNA binding partners of circ_0009043 were thus identified, with bioinformatics analysis revealing miR-148a-3p binding sites within this circRNA. Luciferase reporter assays then confirmed the ability of circ_0009043 to bind miR-148a-3p. This was noteworthy given that miR-148a-3p has previously been reported to exert tumor suppressor activity in several cancers [36, 37], and differential miR-148a-3p expression has been observed in various malignancies such that it is upregulated in certain cancer types while downregulated in others [38,39,40]. Functionally, miR-148a-3p has been reported to play a role in shaping cellular development and differentiation, shaping tumor development through the regulation of particular signaling pathways . Previously, miR-148a-3p upregulation has been reported in glioma and osteosarcoma [42, 43], consistent with its importance in oncogenic settings. Here, miR-148a-3p was found to promote NSCLC cell proliferative activity while preventing apoptotic death, functioning at least in part via targeting the DNAJB4 mRNA and thus preventing it from negatively regulating NSCLC progression.
HSPs are frequently identified as promising therapeutic targets and relevant biomarkers in various forms of cancer [44, 45]. DNAJB4 (HLJ1) is an HSP40 family member that has previously been identified as a tumor suppressor gene in lung, colon, and breast cancer patients [46,47,48], and that has also been specifically identified as a promising therapeutic target and putative biomarker in NSCLC . Consistently, the results of this study revealed that circ_0009043 overexpression resulted in DNAJB4 upregulation, whereas circ_0009043 downregulation had the opposite effect. Notably, the overexpression of miR-148a-3p was sufficient to reverse the impact of circ_0009043 on DNAJB4 upregulation.
Together, these results revealed that circ_000904 is downregulated in NSCLC cells and that, when expressed, this circRNA can serve as a molecular sponge to sequester miR-148a-3p, thereby leading to DNAJB4 upregulation and the suppression of NSCLC tumor growth (Fig. 10). These results demonstrate the tumor suppressor activity exerted by circ_0009043 in NSCLC cells and suggest that it may offer value as a diagnostic/prognostic biomarker in patients with NSCLC as well as a viable target for therapeutic intervention.
Availability of data and materials
The data and study materials that support the findings of this study will be available to other researchers from the corresponding authors on reasonable request.
Non-small-cell lung cancer
Heat shock protein
DnaJ heat shock protein family (Hsp40) member B4
Fetal bovine serum
Cell Counting Kit-8
Terminal deoxynucleotidyl transferase-medimed dUTP nick-endlabeling
Quantitative real-time reverse transcription PCR
Mao Y, Yang D, He J, Krasna MJ. Epidemiology of lung cancer. Surg Oncol Clin N Am. 2016;25(3):439–45.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30.
Herbst RS, Morgensztern D, Boshoff C. The biology and management of non-small cell lung cancer. Nature. 2018;553(7689):446–54.
Duma N, Santana-Davila R, Molina JR. Non-small cell lung cancer: epidemiology, screening, diagnosis, and treatment. Mayo Clin Proc. 2019;94(8):1623–40.
Jonna S, Subramaniam DS. Molecular diagnostics and targeted therapies in non-small cell lung cancer (NSCLC): an update. Discov Med. 2019;27(148):167–70.
Alexander M, Kim SY, Cheng H. Update 2020: management of non-small cell lung cancer. Lung. 2020;198(6):897–907.
Bauer TL, Berkheim DB. Bronchoscopy: diagnostic and therapeutic for non-small cell lung cancer. Surg Oncol Clin N Am. 2016;25(3):481–91.
Kristensen LS, Andersen MS, Stagsted LVW, Ebbesen KK, Hansen TB, Kjems J. The biogenesis, biology and characterization of circular RNAs. Nat Rev Genet. 2019;20(11):675–91.
Chen LL. The expanding regulatory mechanisms and cellular functions of circular RNAs. Nat Rev Mol Cell Biol. 2020;21(8):475–90.
Panda AC. Circular RNAs Act as miRNA sponges. Adv Exp Med Biol. 2018;1087:67–79.
Qiu L, Xu H, Ji M, Shang D, Lu Z, Wu Y, et al. Circular RNAs in hepatocellular carcinoma: biomarkers, functions and mechanisms. Life Sci. 2019;231:116660.
Misir S, Hepokur C, Aliyazicioglu Y, Enguita FJ. Circular RNAs serve as miRNA sponges in breast cancer. Breast Cancer. 2020;27(6):1048–57.
Li R, Jiang J, Shi H, Qian H, Zhang X, Xu W. CircRNA: a rising star in gastric cancer. Cell Mol Life Sci. 2020;77(9):1661–80.
Rong Z, Xu J, Shi S, Tan Z, Meng Q, Hua J, et al. Circular RNA in pancreatic cancer: a novel avenue for the roles of diagnosis and treatment. Theranostics. 2021;11(6):2755–69.
Li Z, Li X, Xu D, Chen X, Li S, Zhang L, et al. An update on the roles of circular RNAs in osteosarcoma. Cell Prolif. 2021;54(1):e12936.
Ye F, Tang QL, Ma F, Cai L, Chen M, Ran XX, et al. Analysis of the circular RNA transcriptome in the grade 3 endometrial cancer. Cancer Manag Res. 2019;11:6215–27.
Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, et al. Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005;33(20):e179.
Li S, Teng S, Xu J, Su G, Zhang Y, Zhao J, et al. Microarray is an efficient tool for circRNA profiling. Brief Bioinform. 2019;20(4):1420–33.
Patop IL, Kadener S. circRNAs in Cancer. Curr Opin Genet Dev. 2018;48:121–7.
Yu T, Wang Y, Fan Y, Fang N, Wang T, Xu T, et al. CircRNAs in cancer metabolism: a review. J Hematol Oncol. 2019;12(1):90.
Chen RX, Liu HL, Yang LL, Kang FH, Xin LP, Huang LR, et al. Circular RNA circRNA_0000285 promotes cervical cancer development by regulating FUS. Eur Rev Med Pharmacol Sci. 2019;23(20):8771–8.
Jiao J, Zhang T, Jiao X, Huang T, Zhao L, Ma D, et al. hsa_circ_0000745 promotes cervical cancer by increasing cell proliferation, migration, and invasion. J Cell Physiol. 2020;235(2):1287–95.
Wei W, Ji L, Duan W, Zhu J. Circular RNA circ_0081001 knockdown enhances methotrexate sensitivity in osteosarcoma cells by regulating miR-494-3p/TGM2 axis. J Orthop Surg Res. 2021;16(1):50.
Lu TX, Rothenberg ME. MicroRNA. J Allergy Clin Immunol. 2018;141(4):1202–7.
Mohr AM, Mott JL. Overview of microRNA biology. Semin Liver Dis. 2015;35(1):3–11.
Kalayinia S, Arjmand F, Maleki M, Malakootian M, Singh CP. MicroRNAs: roles in cardiovascular development and disease. Cardiovasc Pathol. 2021;50:107296.
Wu Y, Li Q, Zhang R, Dai X, Chen W, Xing D. Circulating microRNAs: biomarkers of disease. Clin Chim Acta. 2021;516:46–54.
Lee YS, Dutta A. MicroRNAs in cancer. Annu Rev Pathol. 2009;4:199–227.
Hayes J, Peruzzi PP, Lawler S. MicroRNAs in cancer: biomarkers, functions and therapy. Trends Mol Med. 2014;20(8):460–9.
Zhang M, Han Y, Zhai Y, Ma X, An X, Zhang S, et al. Integrative analysis of circRNAs, miRNAs, and mRNAs profiles to reveal ceRNAs networks in chicken intramuscular and abdominal adipogenesis. BMC Genomics. 2020;21(1):594.
Liu YY, Zhang LY, Du WZ. Circular RNA circ-PVT1 contributes to paclitaxel resistance of gastric cancer cells through the regulation of ZEB1 expression by sponging miR-124–3p. Biosci Rep. 2019;39(12).
Zeng Z, Zhou W, Duan L, Zhang J, Lu X, Jin L, et al. Circular RNA circ-VANGL1 as a competing endogenous RNA contributes to bladder cancer progression by regulating miR-605-3p/VANGL1 pathway. J Cell Physiol. 2019;234(4):3887–96.
Chen Y, Yang F, Fang E, Xiao W, Mei H, Li H, et al. Circular RNA circAGO2 drives cancer progression through facilitating HuR-repressed functions of AGO2-miRNA complexes. Cell Death Differ. 2019;26(7):1346–64.
Gan X, Zhu H, Jiang X, Obiegbusi SC, Yong M, Long X, et al. CircMUC16 promotes autophagy of epithelial ovarian cancer via interaction with ATG13 and miR-199a. Mol Cancer. 2020;19(1):45.
Sang Y, Chen B, Song X, Li Y, Liang Y, Han D, et al. circRNA_0025202 regulates tamoxifen sensitivity and tumor progression via regulating the miR-182-5p/FOXO3a axis in breast cancer. Mol Ther. 2019;27(9):1638–52.
Friedrich M, Pracht K, Mashreghi MF, Jack HM, Radbruch A, Seliger B. The role of the miR-148/-152 family in physiology and disease. Eur J Immunol. 2017;47(12):2026–38.
Li Y, Deng X, Zeng X, Peng X. The role of mir-148a in cancer. J Cancer. 2016;7(10):1233–41.
Wang Y, Hu Y, Guo J, Wang L. miR-148a-3p suppresses the proliferation and invasion of esophageal cancer by targeting DNMT1. Genet Test Mol Biomarkers. 2019;23(2):98–104.
Wang W, Dong J, Wang M, Yao S, Tian X, Cui X, et al. miR-148a-3p suppresses epithelial ovarian cancer progression primarily by targeting c-Met. Oncol Lett. 2018;15(5):6131–6.
Dybos SA, Flatberg A, Halgunset J, Viset T, Rolfseng T, Kvam S, et al. Increased levels of serum miR-148a-3p are associated with prostate cancer. APMIS. 2018;126(9):722–31.
Murata T, Takayama K, Katayama S, Urano T, Horie-Inoue K, Ikeda K, et al. miR-148a is an androgen-responsive microRNA that promotes LNCaP prostate cell growth by repressing its target CAND1 expression. Prostate Cancer Prostatic Dis. 2010;13(4):356–61.
Hua D, Mo F, Ding D, Li L, Han X, Zhao N, et al. A catalogue of glioblastoma and brain MicroRNAs identified by deep sequencing. OMICS. 2012;16(12):690–9.
Ma W, Zhang X, Chai J, Chen P, Ren P, Gong M. Circulating miR-148a is a significant diagnostic and prognostic biomarker for patients with osteosarcoma. Tumour Biol. 2014;35(12):12467–72.
Wu J, Liu T, Rios Z, Mei Q, Lin X, Cao S. Heat shock proteins and cancer. Trends Pharmacol Sci. 2017;38(3):226–56.
Saini J, Sharma PK. Clinical, prognostic and therapeutic significance of heat shock proteins in cancer. Curr Drug Targets. 2018;19(13):1478–90.
Acun T, Doberstein N, Habermann JK, Gemoll T, Thorns C, Oztas E, et al. HLJ1 (DNAJB4) gene is a novel biomarker candidate in breast cancer. OMICS. 2017;21(5):257–65.
Liu Y, Zhou J, Zhang C, Fu W, Xiao X, Ruan S, et al. HLJ1 is a novel biomarker for colorectal carcinoma progression and overall patient survival. Int J Clin Exp Pathol. 2014;7(3):969–77.
Chen HW, Lee JY, Huang JY, Wang CC, Chen WJ, Su SF, et al. Curcumin inhibits lung cancer cell invasion and metastasis through the tumor suppressor HLJ1. Cancer Res. 2008;68(18):7428–38.
Wang CC, Lin SY, Lai YH, Liu YJ, Hsu YL, Chen JJ. Dimethyl sulfoxide promotes the multiple functions of the tumor suppressor HLJ1 through activator protein-1 activation in NSCLC cells. PLoS ONE. 2012;7(4):e33772.
The research is supported by grants from Hunan Provincial Health Commission general funding project (202102080487).
Ethics approval and consent to participate
The protocol was approved by the ethical review committee of Hunan Provincial People’s Hospital, and patient consent was obtained before the samples were taken. And The Ethics Review Committee of the Department of Laboratory Animal Science of Hunan Provincial People’s Hospital approved this study.
Consent to publication
The authors have no relevant financial or non-financial interests to disclose.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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
She, K., Yu, S., He, S. et al. CircRNA 0009043 suppresses non-small-cell lung cancer development via targeting the miR-148a-3p/DNAJB4 axis. Biomark Res 10, 61 (2022). https://doi.org/10.1186/s40364-022-00407-y
- Circular RNA