Skip to main content

Circular RNA ZNF800 (hsa_circ_0082096) regulates cancer stem cell properties and tumor growth in colorectal cancer

Abstract

Background

Cancer stem cells form a rare cell population in tumors that contributes to metastasis, recurrence and chemoresistance in cancer patients. Circular RNAs (circRNAs) are post-transcriptional regulators of gene expression that sponge targeted microRNA (miRNAs) to affect a multitude of downstream cellular processes. We previously showed in an expression profiling study that circZNF800 (hsa_circ_0082096) was up-regulated in cancer stem cell-enriched spheroids derived from colorectal cancer (CRC) cell lines.

Methods

Spheroids were generated in suspension spheroidal culture. The ZNF800 mRNA, pluripotency stem cell markers and circZNF800 levels were determined by quantitative RT-PCR. CircZNF800-miRNA interactions were shown in RNA pulldown assays and the miRNA levels determined by stem-loop qRT-PCR. The effects of circZNF800 on cell proliferation were tested by EdU staining followed by flowcytometry. Expression of stem cell markers CD44/CD133, Lgr5 and SOX9 was demonstrated in immunofluorescence microscopy. To manipulate the cellular levels of circZNF800, circZNF800 over-expression was achieved via transfection of in vitro synthesized and circularized circZNF800, and knockdown attained using a CRISPR-Cas13d-circZNF800 vector system. Xenografted nude mice were used to demonstrate effects of circZNF800 over-expression and knockdown on tumor growth in vivo.

Results

CircZNF800 was shown to be over-expressed in late-stage tumor tissues of CRC patients. Data showed that circZNF800 impeded expression of miR-140-3p, miR-382-5p and miR-579-3p while promoted the mRNA levels of ALK/ACVR1C, FZD3 and WNT5A targeted by the miRNAs, as supported by alignments of seed sequences between the circZNF800-miRNA, and miRNA-mRNA paired interactions. Analysis in CRC cells and biopsied tissues showed that circZNF800 positively regulated the expression of intestinal stem cell, pluripotency and cancer stem cell markers, and promoted CRC cell proliferation, spheroid and colony formation in vitro, all of which are cancer stem cell properties. In xenografted mice, circZNF800 over-expression promoted tumor growth, while circZNF800 knockdown via administration of CRISPR Cas13d-circZNF800 viral particles at the CRC tumor sites impeded tumor growth.

Conclusions

CircZNF800 is an oncogenic factor that regulate cancer stem cell properties to lead colorectal tumorigenesis, and may be used as a predictive marker for tumor progression and the CRISPR Cas13d-circZNF800 knockdown strategy for therapeutic intervention of colorectal cancer.

Peer Review reports

Background

Cancer stem cells constitute a rare population in the heterogeneous cell population in tumors [1, 2]. It has been proposed that cancer stem cells are derived from normal stem cells residing in tissues or organs on assaults by mutations, or stem cell-like properties are acquired in terminally differentiated cells induced by de-differentiation mutations [2, 3]. Cancer stem cells are characterized by enhanced self-renewal abilities and tumorigenic properties, including unhindered cell division capacities [2]. Clinically, cancer stem cells contribute to metastasis, increased tumor aggressiveness, recurrence and chemo- and irradiation-resistance in therapeutic treatment of cancer patients [3, 4].

In colorectal cancer (CRC), the study model of this work, oncogenic mutations in the intestinal stem cell (ISC) population give rise to CRC stem cells to participate in CRC tumorigenesis [2, 5]. Introducing adenomatous polyposis coli mutations reportedly lead to adenoma only in LGR5+ (leucine-rich repeat containing G protein-coupled receptor 5) stem cells, which resembles ISC, but not in LGR5 cell population, marking LGR5 a cancer stem cell marker in CRC [6]. The unhindered CRC stem cell proliferation as a result of gain-of-function mutations in the intestinal stem cells in the crypts of intestines has become a major challenge for therapeutic interventions [7, 8].

The rarity of cancer stem cells is often addressed using cancer cell-derived spheroids, which are enriched in cancer stem cells and are a convenient system for studying cancer stem cell properties and for elucidating the role of cancer stem cells in the tumorigenesis process [9]. Spheroidal culture is an anchorage-independent culture that positively selects cells that resist anoikis, which is programmed cell death induced on loss of cell-to-matrix interactions [10, 11]. Cells resistant to anoikis have been shown to exhibit cancer stem cell-like properties in various cancers [12]. The resulting spheroids can be propagated for many passages to increasingly enrich cells with enhanced cancer stem cell-like properties [13].

Of the many molecular mechanisms that regulate tumorigenesis, post-transcriptional regulation by circular RNA (circRNA) is important it is at the apex that controls other downstream regulatory mechanisms [14]. CircRNAs are back-splicing products of the pre-mRNA of a host transcript; multiple isoforms of different combinations of exon and intron sequences may be derived from a single host transcript, targeting and sponging off another group of post-transcriptional regulators, the microRNA (miRNA) [14,15,16]. CircRNA also functions via interacting with RNA-binding proteins, or in being translated into functional proteins to impact on downstream biochemical and cellular pathways [14, 16, 17]. Dysregulation of circRNAs in cancer has been associated with the enhancement of various cancer hallmarks [18].

In our previous expression profiling analysis of circRNAs in spheroids derived from two CRC cell lines, HCT-15 and WiDr, hsa_circ_0082096 was identified as one of top four up-regulated circRNAs in the spheroids compared with the parental CRC cells; bioinformatics analysis further predicts association of hsa_circ_0082096 with stem-cell properties [19]. The present work aimed to experimentally validate and to elucidate the role of hsa_circ_0082096 in modulating cancer stem cell properties and in the tumorigenesis process of colorectal cancer.

Materials and methods

Cell and spheroidal cultures and CRC patient tissues

All colorectal cancer (CRC) cell lines were obtained from American Type Culture Collection and cultured as suggested. The cell lines were regularly monitored for mycoplasma contamination. CRC spheroids were generated as previously described [19]. Primary CRC and adjacent normal paired tissues were obtained with consents from patients who underwent colectomy in Taipei Veterans General Hospital, Taiwan, between 2019–2022. in accordance with the International Ethical Guidelines for Biomedical Research Involving Human Subjects and was approved by the Medical Ethical Committee of Taipei Veterans General Hospital, Taiwan (IRB-2019–08-009B). The clinical CRC samples (n = 30) included 22 cases of patients with age ≥ 60 years old and 8 cases of age < 60, and 17 male and 13 female patients.

RNA extraction, quantitative RT-PCR and sequencing analysis

RNA extraction using TRIzol reagent (Invitrogen, Thermo Fisher Scientific, USA), cDNA synthesis using SuperScript® III Reverse Transcriptase Kit (Invitrogen), quantitative reverse transcription-PCR (qRT-PCR) using SYBR® Select Master Mix (Invitrogen) using glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as an internal control, and Sanger sequencing were as described [19]. The primers used are shown in Suppl. Table S1A.

RNA pulldown and stem-loop qRT-PCR

The RNA pulldown protocol was modified from a previous work [20]. CRC cell lysates were collected using radio-immunoprecipitation assay buffer (RIPA) buffer and quantified. An aliquot of 2,000 µg cell lysate was incubated with 500 ng biotinylated circZNF800, or a scrambled control RNA. The mixture was gently agitated at 4 ºC for 4 h. Streptavidin Mag Sepharose beads (Cytiva, Massachusetts, USA) were incubated with the lysate overnight. After three washes using the bead washing buffer supplemented with RNase inhibitor provided by the vendor, the beads were resuspended in TRIzol and RNA was extracted and quantified by stem-loop qRT-PCR [21]. RNA (1 µg) was reverse transcribed with the stem-loop primers (Suppl. Table S1B) and qRT-PCR was performed as described in “RNA extraction, quantitative RT-PCR and sequencing analysis” section above, using U6 as the PCR control. Relative miRNA expression levels were calculated using the comparative Cт (ΔΔCт).

Actinomycin D and RNase R treatments

CRC cells seeded at 5 × 105 cells/well in a 6-well plate were incubated overnight. After treatment with 2.5 µg/mL of actinomycin D (Sigma-Aldrich, Merck, Germany), the cells were harvested for RNA extraction. For RNase R treatment, RNA extracted were treated with RNase R (Lucigen, Thermo Fisher Scientific, USA) based on the manufacturer’s instructions. The RNA was pelleted using lithium chloride precipitation (Thermo Fisher Scientific) and was resuspended in RNase-free water for further analysis.

CircZNF800 synthesis in vitro

CircZNF800 synthesis in vitro was done using a permutated group I self-catalytic of td (thymidylate synthase) gene intron plasmid, driven by T7 RNA polymerase, a gift from Professor Manuel Ares [22]. The circZNF800 cDNA, or the green fluorescence protein (GFP), sequence was cloned at the AflII and MfeI restriction sites between the 3’ and 5’ splice sites of the plasmid. In vitro transcription was performed using the HiScribe™ T7 Quick High Yield RNA Synthesis Kit, (New England Biolabs, Ipswich, MA, USA) per manufacturer’s instructions. The resulting RNA was purified by lithium chloride precipitation. For modified RNA, cyanine-3-uridine-5'-triphosphate (Cy3-UTP) was added at a 1:4 ratio. RNA circularization was induced by first heating the RNA to 65 ºC for 3 min and placed on the ice for 10 min. The RNA was mixed with guanosine-5'-triphosphate (GTP) and MgCl2 at a final concentration of 10 mM in 1 × polynucleotide kinase reaction buffer (PNK) buffer (New England Biolabs). The mixture was heated at 55 ºC for 35 min. The RNA was purified using lithium chloride precipitation and resuspended in nuclease-free water for analysis in an agarose-formaldehyde gel (2%) in 3-(N-morpholino) propanesulfonic acid (MOPS) buffer and deionized. The agarose gel was visualized under UV light.

Preparation of CRISPR-Cas13d circZNF800 constructs

The CRISPR Cas13d plasmid was purchased from Addgene (#138148, Massachusetts, USA). Two CRISPR RNA sequences, crRNA1 and crRNA2, were designed and synthesized based on the backsplice junction of circZNF800 and were cloned into the BsmBI restriction site of the plasmid. The constructs were sequence verified. A control construct, crSC, which harbored a scrambled sequence, was also prepared.

DNA and RNA transfection

CRC cells were seeded at 5 × 105 cells/well in a 6-well plate and incubated overnight. The DNA or RNA transfection was performed using Lipofectamine 2000 (Invitrogen) per manufacturer’s instructions.

RNA FISH and dual immunofluorescent assay

A cyanine-5 (Cy5)-tagged circZNF800 antisense RNA was commercially synthesized (MDBio Inc, Taiwan), and was mixed with the hybridization buffer at a final concentration of 200 nM in an overnight incubation. RNA FISH protocol was modified from previous works [23]. After hybridization and washing (see Suppl. Table S1C for probe information), the slides were stained with 4’,6-diamidino-2-phenylindole (DAPI) (Invitrogen, Thermo Fisher) at 1:10,000 dilutions in phosphate buffered saline (PBS) for 10 min. Fluorescent images were captured using a FV3000 Confocal Laser Scanning Microscope (Olympus, Japan). For dual RNA-FISH and protein immunofluorescent staining, upon overnight incubation with the circRNA probe as described above, the tissue slides were rinsed with saline sodium citrate (SCC) buffer before incubation with antibodies for Lgr5 (ab219107, Abcam, UK), Sox9 (14–9765-82, eBioscience, Thermo Fisher Scientific) or Ki-67 (#9129, Cell Signalling Technology, MA, USA). Goat anti-mouse IgG secondary antibody conjugated with Alexa Fluor 546 and goat anti-rabbit IgG secondary antibody conjugated with Alexa Fluor 488 (Invitrogen) were used. The fluorescence images were captured as described above.

Flowcytometric analysis of protein markers and cell proliferation rate

CRC cells were harvested using Accutase (Gibco, Thermo Fisher) in 5 × 105 cells aliquots. The antibodies used were CD44 (#3570, Cell Signalling Technology), CD133 (#5860, Cell Signalling Technology), Lgr5 and SOX9. Cell proliferation rate was assessed using a commercial EdU Staining Proliferation Kit (fluor 488) (Abcam, UK) per manufacturer’s instructions. The treated cells were collected through centrifugation and washed with once with Permeabilization Buffer. The cell pellet was resuspended in 300 µl PBS and data acquisition was achieved by flowcytometry as described above.

Colony forming assay

CRISPR Cas13d-transfected CRC cells were seeded at 2.5 × 103 cells in a 6-well plate. The cells were cultured in the respective media and replenished every other day. The cells were allowed to form colonies for 10 – 14 days. Crystal violet staining for 20 min at room temperature was used to visualize the colonies formed.

Analysis of the effects of circZNF800 over-expression and knockdown in mice

BALB/c nude male mice (strain CAnN.Cg-Foxn1nu/CrlBltW) were purchased from Bio-LASCO Taiwan, Taiwan, and were housed in cages with free access to food source and water. The animal work was conducted according to the Guidelines for Laboratory Animals in the Taipei Veterans General Hospital, with the institutional ethical review board approval (IACUC 2012–079).

To evaluate the effects of circZNF800 over-expression on tumor growth, HCT-15 and WiDr cells in PBS were first injected subcutaneously at 2.5 × 106 cells per dorsal flank of nude mice until tumors were formed to an average size of ~ 500 mm3. The mice were divided into the vehicle (saline), circGFP RNA- and circZNF800-treated groups (n = 3 per group). RNA (10 µg) in TransIT®-QR Delivery Solution (Mirus Bio, Wisconsin, USA), or 10 µL saline, was diluted in 90 µL of TransIT®-QR Delivery Solution to a final volume of 100 µL. The mixture was injected nto the tumor sites within 15 min of preparation (n = 3). The treatment was performed 4 times over 10 days, after which the tumors were left to grow for 20 days, and tumor sizes were measured every 4 days using a caliper. The mice were sacrificed and the tumors were harvested for imaging.

To analyze the effects of circZNF800 knockdown in mice, circZNF800 was stably knocked down in HCT-15 and WiDr cells by lentivirus transduction of CRISPR-Cas13d circZNF800 knockdown constructs crRNA1 and crRNA2 and the control crSC. The cells were suspended in PBS for injection at 2.5 × 106 cells per flank of nude mice. The tumors were left to grow for 24–30 days while tumor sizes were monitored every 3–5 days.

Analysis of the effects of intratumoral circZNF800 knockdown in mice

For intratumoral circZNF800 knockdown, purified and quantified lentiviruses wereobtained from the RNAi Core, Academia Sinica, Taiwan. Each flank was administered 1 × 108 lentivirus particles carrying the control crSC or the targeting crRNA1 or crRNA2 (n = 2). The treatment was performed 4 times over 10 days. The tumors were left to grow for 16–22 days in WiDr and HCT-15 cells, respectively, due to different growth rates. Tumor sizes were monitored before the tumors were harvested for imaging.

Bioinformatics analysis

Messenger RNA sequences were derived from GenBank, NCBI; https://www.ncbi.nlm.nih.gov/genbank/); circRNA sequences and isoforms were derived from circBase (http://www.circbase.org) and circBank (http://www.circbank.cn); miRNA sequences were derived from miRBase (https://www.mirbase.org). CircRNA-miRNA interactions were analyzed using the CircInteractome database (https://circinteractome.nia.nih.gov). MiRNA-mRNA interactions were analyzed by using the TargetScanHuman 7.0 database (https://www.targetscan.org/vert_70/). All bioinformatics analyses using the above databases were done in January–February 2023.

Statistical analysis

All experiments were performed in triplicates. The triplicated data were evaluated for statistical significance using the standard paired Student’s t-test. Comparisons were made between the control group and the treated group, unless specified otherwise. The Prism GraphPad software was used to evaluate the significance levels and for graph construction.

Results

CircZNF800 expression is up-regulated in in CRC-derived spheroids and CRC tumor tissues

In our previous study, we have identified hsa_circ_0082096 as one of top four up-regulated circRNAs in spheroids derived from the HCT-15 and WiDr CRC cell lines [19]. Based on circBase and circBank databases, the host transcript of hsa_circ_0082096 is ZNF800 (zinc-finger protein 800). The ZNF800 gene (NM_176814 and XM_01151585.4) is composed of seven exons, with an alternative exon 6a or 6b, and is predicted to generate eight circRNA isoforms, including hsa_circ_0082096, based on expression profiling studies (Suppl. Table S2) [19, 24,25,26]. The hsa_circ_0082096 isoform that we reported earlier is derived from the back-splicing of exons 4 & 5 of the ZNF800 transcript 1,837 nucleotides in length and is designated as hsa_circZNF800, or circZNF800 for short (Fig. 1A). It is noted in the various expression profiling reports that circZNF800 is also expressed in the chronic myeloid leukemia cells K562 and in SK-N-SH-RA, a neuroblastoma SK-N-SH cells differentiated with retinoic acid, and in normal fetal lung fibroblast AG04450 cells, EBV-transformed lymphoblastoid cells, chondrocytes and in multiple brain tissues and neuronal cells (Suppl. Table S2). Taken together, expression profiling reports, including our previous work, have predicted participation of circZNF800 in tumorigenesis, cancer stem cells and the functioning of brain tissues. CircZNF800 in CRC cells was confirmed by RNase R resistance (Fig. 1B) and significantly longer half-life than the host transcript ZNF800 on actinomycin D challenge (Fig. 1C) in the two CRC cell lines tested.

Fig. 1
figure 1

CircZNF800 expression is up-regulated in in CRC-derived spheroids and CRC tumor tissues. A Schematic illustration of circZNF800 biogenesis via backsplicing of exons 4 and 5 of the ZNF800 pre-mRNA (NC_000007.14). Two ZNF800 alternative transcripts carry either exon 6a (NM_176814.5) or 6b (XM_011515859.4) are shown. F and R indicate the forward and reverse primers used in sequencing of the backsplicing junction. B, C Stability analysis of circZNF800 in WiDr cells on RNase R (B) and actinomycin D (C) treatments. The relative RNA levels were determined by qRT-PCR. *P < 0.05 and **P < 0.01 relative to those of the untreated sample (B) or the ZNF800 transcript (C). D Morphology of the CRC-derived spheroids. E Up-regulated circZNF800 expression in spheroids derived from CRC cell lines. F Passage (P)-dependent up-regulated circZNF800 expression in CRC cell-derived spheroids and on serum-induced differentiation of the spheroids at passage 5 (Diff. P5). G, H Co-localization of circZNF800 and intestinal stem cell markers LGR5 and SOX9 in the adjacent normal tissue (G) or Ki-67 high proliferative marker (H) in a CRC tumor tissue in dual RNA FISH-immunofluorescent assays. Bars represent 100 µm. I Significant circZNF800 over-expression in paired tumor and normal tissues of CRC patients, and in late-stage CRC tumors. Paired tissue samples (n = 30) were subjected to qRT-PCR analysis to determine the circZNF800 and ZNF800 expression levels (left two panels); the circZNF800 expression data from the same patient cohort were further grouped according to pathological stages 1 & 2 (n = 17) and 3 & 4 (n = 13) (right two panels) for circZNF800 expression analysis. *p < 0.05, **p < 0.01 and ***P < 0.001 relative to the control samples as indicated or the adjacent normal tissues

Since our previous report on circZNF800 in two CRC cell lines [19], we confirmed here that circZNF800 expression was consistently up-regulated in spheroids generated from a panel of four other CRC cell lines (Fig. 1D and E). Furthermore, circZNF800 up-regulation was passage-dependent, and was down-regulated on serum-induced differentiation in passage 5 cells (Fig. 1F), supporting involvement of circZNF800 in modulating CRC stem cell properties in maintaining the spheroidal phenotype. The cancer stem cell population in CRC is thought to have arisen from a mutated ISC population [5], and that LGR5 and SOX9 are ISC markers that reside in the crypt base columnar cells of the tissue [27]. Dual immunofluorescence assays showed that circZNF800 co-resided in the LGR5high/SOX9high cell population in non-pathological colorectal tissues (Fig. 1G), supporting association of circZNF800 with ISC. Dual immunofluorescence assays further showed co-location of circZNF800 and the Ki-67 high proliferation marker in the same tumor cell population (Fig. 1H) [28]. In clinical biopsies, circZNF800 expression, determined by qRT-PCR, was significantly higher in CRC tumor tissues than in the paired adjacent normal tissues analyzed (n = 30), whereas the ZNF800 transcript levels were not significantly altered between the tumor and normal tissues (Fig. 1I, left two panels). However, when the circZNF800 expression data were further grouped according to pathological stages 1 & 2 (n = 17) and 3 & 4 (n = 13), significant up-regulated circZNF800 expression was found in the late stages 3/4 tumor tissues, but not in the early stages 1/2 samples (Fig. 1I; right two panels). Taken together, the data show that up-regulated circZNF800 expression is associated with cancer stem cell properties and CRC tumorigenesis.

Strategies for circZNF800 over-expression (OE) and knockdown (KD) in cells

To further dissect molecular and cellular events affected by circZNF800, experimental approaches to manipulate cellular levels of circZNF800 were developed. For over-expression (OE), an in vitro scheme of generating self-circularizing circular RNA using permutated T4 phage Td gene was used [22, 29]. The RNA products obtained were ~ 1,800 nucleotides as expected and were RNase R resistant supporting circularity (Fig. 2A; see also Suppl. Fig. S1). The in vitro transcribed and circularized circZNF800 was confirmed to carry exons 4 and 5 sequences by sequencing. To assess cellular uptake, a modified UTP-conjugated circZNF800 RNA tagged with fluorescence Cy3 was transfected in WiDr cells and was visualized by fluorescence microscopy 12 h post-transfection. The fluorescence images showed that the tagged circZNF800 overlapped with the WiDr cells (Fig. 2B), suggesting successful circZNF800 uptake by the cells.

Fig. 2
figure 2

CircZNF800 sponges miR-140-3p, miR-382-5p & miR-579-3p to up-regulate expression of ALK7, FZD3 & WNT5A. A, B Properties of in vitro transcribed and circularized circZNF800: RNase R resistance shown in a 2% formaldehyde-agarose gel (A) and co-localization of circZNF800-Cy3 RNA transfected in WiDr cells observed under a fluorescence microscope (B). C, D For CRISPR-Cas13d knockdown of circZNF800, the two CRISPR RNA (crRNA) constructs, crRNA1 and crRNA2 (C), and the knockdown efficiency (D) are shown. In (D), CrSC was a crRNA construct with scrambled sequence used as a reference control. Relative expression levels of circZNF800 and ZNF800 post-crRNA transfection were determined by qRT-PCR, normalizing to the levels of the control crSC-transfected cells. **p < 0.01 and ***p < 0.001 were values relative to the CrSC cells. E Mapping of the miRNA seed sequence-interacting sites on circZNF800 (top panel) and the 3’-UTR sequences of the targeted transcripts (bottom panel) (see Suppl. Figure S1 for sequence alignment details). The relative positions of the miRNAs, the GenBank accession numbers of the mRNA sequences and the 3’-UTR lengths are shown. CDS, coding sequence; An, polyA tail. F RNA pulldown assays using biotinylated circZNF800 in HCT-15 and WiDr cells. The pulled-down RNA products were quantified in stem loop qRT-PCR (see Suppl. Table S1 for primer details) using a scramble RNA as a control and normalized to 1.0. G, H CircZNF800 over-expression (OE) sponges while knockdown (KD) up-regulates expression of the targeted miRNAs (G) and mRNAs (H) in CRC cells relative to the controls without OE and KD transfection, normalized to 1.0 as shown by the horizontal control (Ctrl) dashed lines. OE and KD were achieved by transfection of the in vitro synthesized circZNF800 and the CRIPSR constructs, respectively, in CRC cells for 48 h before RNA quantification by stem loop qRT-PCR relative to the control values (horizontal Ctrl dashed lines). *P < 0.05, **P < 0.01 and ***P < 0.001 were values relative to the specified controls

For circZF800 knockdown (KD), the CRISPR-Cas13d system, which induces specific circRNA degradation, was used [30, 31]. Two CRISPR RNA (crRNA) constructs, crRNA1 and crRNA2, that span different sequences of the circZNF800 back-splice junction, were designed (Fig. 2C). A similar construct with a scrambled sequence, designated crSC, was used as a control. Both crRNAs were able to induce ~ 40% significant knockdown of circZNF800 in two CRC cell lines with little effects on the ZNF800 transcript compared with the crSC control-treated cells (Fig. 2D).

CircZNF800 sponges miR-140-3p, miR-382-5p and miR-579-3p and up-regulates expression of ALK7/ACR1C, FDZ3 and WNT5A

In our previous expression profiling analysis, miR-140-3p, miR-382-5p and miR-579-3p were predicted to be sponged by circZNF800 to collectively affect expression of the ALK7/ACR1C (activin A receptor type 1C), FDZ3 (frizzled class receptor 3) and WNT5A (WNT family member 5A) genes in spheroids derived from CRC cell lines [19]. The predicted circZNF800-miRNA-mRNA interactions are supported by sequence alignments of the seed sequences of the circZNF800-miRNA and miRNA-mRNA pairs using the CircInteractome data (Fig. 2E). While miR-382-5p is mapped in the short exon 4, the other two miRNAs are mapped in exon 6 of the ZNF800 transcript (Fig. 2E, top panel; Suppl. Fig. S2A for sequence alignment details). On miRNA-mRNA alignments using the TargetScan database, the 3’-untranslated region (3’-UTR) of the 7.15-kb ALK7 transcript is found to be targeted specifically by the seed sequence of miR-382-5p, while the 4.38-kb 3’-UTR sequence of WNT5A is targeted by the seed sequences of miR-140-3p once and miR-579-3p at four sites (Fig. 2E). Interestingly, the 11.25-kb 3’-UTR of the FZD3 transcript is targeted three times by miR-382-5p and ten times by miR-579-3p, which are clustered in the 5’-half of the UTR sequence (Fig. 2E). The seed sequences of the miRNA-mRNA interactions are 7–8 nucleotides in length and in perfect sequence homology, supporting with high confidence the bioinformatics-predicted miRNA-mRNA interactions (Suppl. Fig. S2B). To validate the predicted circZNF800-miRNA interactions, biotinylated-circZNF800 and a scrambled RNA control generated in transcription in vitro were used in RNA pulldown assays in HCT-15 and WiDr cells. The pulled-down RNA products were purified and subjected to stem-loop qRT-PCR analysis. The results indicated that the circZNF800-pulldown products were enriched in the predicted miRNAs in the HCT-15 and WiDr cells compared to the scrambled RNA control (Fig. 2F).

Using the over-expression and knockdown methodologies described in Fig. 2A-D above, ectopic circZNF800 OE in CRC cells was shown to result in significant down-regulation of the expression levels of all three miRNAs while circZNF800 KD up-regulated the expression of the miRNAs compared to the control (Fig. 2G), collaborating with circZNF800 sponging of the miRNAs. Furthermore, circZNF800 OE significantly up-regulated the mRNA levels of ALK7, FZD3 and WNT5A while the reverse was observed on circZNF800 KD (Fig. 2H), supporting circZNF800-miRNA and miRNA-mRNA interactions (Fig. 2E). The data collectively show that by sponging miR-140-3p, miR-382-5p and miR-579-3p, circZNF800 promotes ALK7, FZD3 and WNT5A expression, possibly contributing to cancer stem cell properties of CRC cells.

CircZNF800 up-regulates CSC markers in CRC cells

The pluripotency markers NANOG and the Yamanaka factors OCT4, SOX2, KLF4 and c-MYC have previously been shown to be linked with enhanced cancer stem cell properties in cancer cells [32,33,34]. On circZNF800 over-expression, OCT4, SOX2 and NANGO were significantly up-regulated in the two CRC cell lines tested in qRT-PCR analysis; on circZNF800 knockdown, the same three pluripotency factors were also down-regulated (Fig. 3A). The data support circZNF800-mediated modulation of expression of the core regulators, OCT4, SOX2 and NANOG, in maintaining stemness properties in the CRC population. The impacts of circZNF800 OE and KD on the expression levels of the CSC markers, CD133 and CD44, and the intestinal stem cell (ISC) markers, LGR5 and SOX9 (see above) were also investigated by flowcytometry (Fig. 3B and C). CD133+ cells isolated from CRC cell populations have been shown to mimic properties of CRC cancer stem cells and CD44+ cells exhibit enhanced tumorigenicity and cell proliferation in CRC [35, 36]. The results showed that the expression levels of CD133 and CD44 were up-regulated on circZNF800 OE, and down-regulated on circZNF800 KD (Fig. 3B). The expression levels of LGR5 and SOX9 were also up- or down-regulated on circZNF800 OE and KD, respectively (Fig. 3C), echoing the previous observation that LRG5 and SOX9 are enriched in the normal colorectal tissue (Fig. 1G). These results support circZNF800 regulation of the expression of these cancer stem cell markers in CRC cells.

Fig. 3
figure 3

CircZNF800 up-regulates expression of CSC markers in CRC cells. Effects of circZNF800 over-expression (OE) (left panels) and knockdown (KD) (right panels) on the expression levels of the pluripotency markers determined by qRT-PCR (A), the CSC markers CD44 and CD133 (B) and the intestinal stem cell markers LGR5 and SOX9 determined by flowcytometry (C). Representative flowcytometry histograms are shown. Half-offset plots are used to illustrate the shifts of the markers. Scales of 1–0 units used in all cases were algorithmically determined by FlowJo software. Bar charts represent triplicated quantitative analysis. *p < 0.05 and **p < 0.01 were values relative to the scrambled CrSC or vehicle control group

CircZNF800 promotes cell proliferation and sphere and colony formation in CRC cells in vitro

The cellular effects of circZNF800 over-expression and knockdown on other cancer stem cell properties were next investigated. The effects on cell proliferation were evaluated by EdU assay and flowcytometric analysis. The results showed increased EdU-positive cell population on circZNF800 over-expression compared to the vehicle and the mock circularized green fluorescence protein (GFP) RNA, circGFP; the data were supported by a reverse effect on knocking down circZNF800 (Fig. 4A). The EdU assay data indicated circZNF800-mediated promotion of cell proliferation, in line with the clinical findings that circZNF800 was up-regulated in CRC tumor samples (Fig. 1I). Furthermore, circZNF800 knockdown in CRC cells led to impeded spheroid and colony forming abilities (Fig. 4B and C), supporting that silencing circZNF800 suppresses cancer stem cell phenotype and cellular transformation in vitro. It was also noted that the spheroids formed on circZNF800 knockdown were smaller and irregular in morphology, indicating compromised cellular processes. These data collectively demonstrate that circZNF800 modulates cellular processes pertaining to cancer stem cell properties, leading to enhanced tumor growth in vitro. Due to the high toxicity observed in circZNF800 over-expression cells on extended period of culture before data collection, the colony and spheroid formation assays were not performed on circZNF800 ectopic over-expression cells.

Fig. 4
figure 4

circZNF800 promotes cell proliferation and sphere and colony formation in CRC cells in vitro. A Effects of circZNF800 on cell proliferation. On circZNF800 over-expression (OE) (left panels) or knockdown (KD) (right panels), the cells were treated with EdU and analyzed by flowcytometry. Representative flowcytometry histographs are shown. Bar charts represent triplicated quantitative analysis. *P < 0.05, **P < 0.01, ***P < 0.001 relative to the respective controls. B Spheroid formation on circZNF800 knockdown. Images were taken after 10 days in the spheroid culture conditions; spheroid quantification is shown (bar represents 100 µm). Colony quantification is also shown. C Colony formation of CRC cells on circZNF800 knockdown. Images show the resulting colonies on staining with a crystal violet solution

CircZNF800 over-expression promotes while knockdown suppresses tumor growth in vivo

The effects of circZNF800-mediated effects observed in vitro were further assessed in vivo in xenograft mice (Fig. 5). For over-expression analysis, WiDr cells were first transplanted subcutaneously to nude mice, and the tumors were allowed to grow to ~ 500 mm3 before injecting the in vitro synthesized circZNF800 RNA at the tumor sites four times over 7 days. Growth of the tumor were monitored up to 20 days starting from the first RNA injection (Fig. 5A, top panel). CircZNF800-treated mice showed significant enhanced tumor growth rate and developed larger tumors compared to the two control groups injected with saline (vehicle group) or circularized EGPF (Fig. 5A). The data echo the in vitro data that circZNF800 promotes cellular proliferation (Fig. 4A).

Fig. 5
figure 5

CircZNF800 over-expression promotes while knockdown suppresses tumor growth in vivo. A CircZNF800 over-expression promotes tumor growth in vivo. WiDr cells were first transplanted subcutaneously into nude mice. The mice were divided into three groups: the vehicle group (saline) and the two groups injected either with circZNF800 RNA or circEGFP RNA (n = 3 for each group). When the tumors reached ~ 500 mm3, each group was injected four times over 10 days with 10 µg RNA for each group. Tumors were harvested on day 20 from the start of the RNA injection (schematic panel). Tumor growth rates and images of the harvested subcutaneous CRC tumors on day 20 of treatment are shown. B CircZNF800 knockdown suppresses tumor growth in vivo. Stably knockdown circZNF800 CRC cells were subcutaneously transplanted in nude mice and the tumors were harvested on day 24 (HCT-15 cells) or on day 30 (WiDr cells) post-cell injections (n = 3). Tumor growth rates and images of the harvested subcutaneous CRC tumors are shown. C Intratumoral administration of the circZNF800-knockdown lentiviral CRISPR-Cas13d constructs retards tumor growth in mice. HCT-15 or WiDr cells were subcutaneously transplanted in nude mice and the tumors were allowed to grow to an average size of ~ 300 mm3. Purified lentivirus constructs of either the control crSC, or crRNA1 or crRNA2 were injected at the tumor sites twice a week for 10 days. The tumors were allowed to continue to grow up to 16 or 22 days for WiDr or HCT-15, respectively, before harvesting for analysis (n = 2). *p < 0.05 and **p < 0.01 were relative to the crSC or vehicle-treated group

The effects of circZNF800 knockdown on tumor growth were assessed in two different experimental approaches. In the first approach, stable circZNF800-knockdown HCT-15 and WiDr cells were first generated using the CRISPR Cas 13d-based crRNA1 or crRNA2, and crSC negative control carrying a scrambled sequence; subsequently, the KD CRC cells were subcutaneously transplanted in nude mice, and the tumor growth was monitored for up to 24 days for HCT-15 or 30 days for WiDr cells due to different growth rates of the two cell lines (Fig. 5B, top panel). In both cell lines, tumor growth on the crRNA1- and crRNA2-treated cells was significantly slower than that in the scrambled control crSC cells (Fig. 5B).

In the second approach, the effects of circZNF800 knockdown in tumors in nude mice developed from subcutaneously grafted CRC cells were investigated. When the tumor grew to ~ 300 mm3, each tumor site was injected four times each with 1 × 108 purified lentiviral particles harboring crRNA1 or crRNA2, or the EGFP or crSC control constructs (Fig. 5C, top panel). Due to the different growth rates, the HCT-15 and WiDr-injected mice were monitored for 22 or 16 days, respectively, by measuring the tumor size regularly before the mice were sacrificed. Using the crSC as a control reference, tumor growth in the crRNA1- and crRNA2-treated mice was significantly retarded (Fig. 5C). Despite the small number of mice treated with each crRNA (n = 2 each), two KD constructs were used and each in two different cell lines (total n = 8), the data collectively support that intratumoral circZNF800 knockdown leads to suppressed tumor growth in the mice. The in vivo data obtained are in line with the in vitro data that circZFN800-knockdown suppresses cell proliferation, colony and spheroid formation (Fig. 4B-C). The data may be a first hint to indicate the possibility of delivering crRNA1- and crRNA2-like viral constructs for therapeutic treatment of colorectal cancer in humans.

Discussion

Expression of circZNF800 is previously shown to be up-regulated in CRC-derived spheroids [19]. We show here that circZNF800 is up-regulated in expression in late-stage CRC tumors (Fig. 1I), which may be enriched in cancer stem cells, consistent with circZNF800 over-expression in the spheroids. By manipulating the cellular levels of circZNF800 via over-expressing or knockdown of the circRNA (Fig. 2A-D), circZNF800 was shown to suppress expression of miR-140-3p, miR-382-5p and miR-579-3p and to up-regulate expression of the ALK7/ACVR1C, FZD3 and WNT5A, as predicted by bioinformatics analysis (Fig. 2E-H). CircZNF800 was further shown in in vitro experiments to positively modulate the expression of pluripotency and cancer stem cell markers (Fig. 3). Functionally, the circRNA promoted cell proliferation, spheroid and colony formation (Fig. 4). The in vitro data are supported by in vivo experiments in which ectopic circZNF800 over-expression led to accelerated while circZNF800 knockdown retarded tumor growth in xenograft mice (Fig. 5). The data collectively support the hypothesis that circZNF800 plays a major role in the CRC tumorigenesis process via modulating cancer stem cell properties, as summarized in Fig. 6.

Fig. 6
figure 6

Schematic summary of the predicted role of circZNF800 in regulating cancer stem cell properties via multiple miRNA-mRNA axes in CRC tumorigenesis. The up- (green) and downward (red)-pointing thick arrows indicate up- and down-regulated expression, respectively, of the indicated parameters demonstrated in this work. Briefly, circZNF800 in late-stage CRC or CRC-derived spheroidal cells suppresses miR-140-3p, miR-382-5p and miR-579-3p to up-regulate expression of ALK7, FZD3 and WNT5A (in grey boxes). On the other hand, circZNF800 up-regulates the expression of cancer stem cell markers (in white boxes) and various cellular processes (green boxes) leading to enhanced stemness properties of CRC cells in the tumorigenesis process. The dashed-line arrow indicates the predicted connection of the circZNF800-miRNA-mRNA axes with the affected downstream cellular events

Dysregulated expression of circRNA has been implicated in the tumorigenesis of various types of cancer [37]. Data presented in this study indicate that over-expression circZNF800 contributes to CRC tumorigenesis via exerting its detrimental regulatory effects in the intestinal stem cells and the high proliferative cell population (Fig. 1G and H), consistent with circZNF800-mediated enhancement of cancer stem cell-like properties in CRC. Involvement of other circRNAs in mediating stemness has been reported in other cell types, including induced pluripotent and cancer stem cells [38]. Over-expression of circBIRC6 promotes and maintains the pluripotency state of human embryonic stem cells by regulating the expression levels of OCT4 and NANGO [39]. Likewise, circNOTCH1 over-expression in gastric cancer stem cells is linked with maintenance of stemness [40]. In bladder cancer, however, circGPRC5A is shown to be translated into a GPRC5A protein to regulate the bladder cancer stem cell phenotype in promoting tumorigenicity [41]. On the other hand, the host transcript ZNF800 of circZNF800 has been linked to regulation of apoptosis, migration and invasion of lung cancer cells, and acts as an adipogenesis master regulator associated with cardio-metabolic traits [42, 43].

The three miRNAs sponged by circZNF800 and the three miRNA-targeted mRNAs have been reported to modulate cancer stem cell-related properties, including cellular proliferation and hence, tumor progression, migration and invasion in metastasis, relapses and chemosensitivity in other cancer types, as summarized in Table 1. Besides circZNF800, only miR-382-5p is also reported to be sponged by another circular RNA circUBAP2 to regulate GPRC5A (G protein-coupled receptor class C group 5 member A) expression in glioma [44]. Moreover, miR-382-5p is also most versatile in also targeting four other mRNAs, YB-1, ANXA3, KLF12 and HIPK3, besides GPR5A (Table 1) [45,46,47]. In the targeted mRNA group, the most prevalent transcript associated with effects on cancer stemness is WNT5A that acts via the WNT/β-catenin signaling pathway (Table 1) [48,49,50,51,52]. Collectively, the circZNF800-regulated miRNAs and mRNAs reported in various cancer types affect cancer cell stem properties via the WNT/β-catenin, PI3K/AKT and TGFβ-SMAD signaling pathways, amongst others (Table 1) [44,45,46,47,48,49,50,51,52,53,54,55,56,57]. In short, circZNF800 acts as a switch that negatively controls the function of a specific group of downstream miRNAs, which, in turn, regulate expression of targeted downstream genes to promote cancer stemness, contributing to the many features of human cancers.

Table 1 Other circRNA-miRNA-mRNA axes and affected signaling pathway and cancer stemness reported for the circZNF800-regulated circRNAs and mRNAs

We also demonstrated here that CRISPR Cas13d-mediated circZNF800 knockdown was highly specific and induces suppressed cell proliferation and colony formation in vitro and tumor growth in vivo in CRC cells, as has been reported for the silencing other circRNAs in other cancer types [30, 31, 58]. The high sequence specificity of the CRISPR Cas13d system is assured by the observation that a single mismatch in the CRISPR RNA-targeting sequence of a circRNA could significantly reduce the knockdown efficiency, thus, avoiding non-specific knockdown [30]. Integrating the CRISPR Cas13d-mediated circZNF800 knockdown system reported here with other more human compatible and effective delivery vectors, such as the adeno-associated virus, could be clinically feasible in intervention treatment of colorectal cancer [59].

Conclusions

We show in this work that circZNF800 is over-expressed in late-stage colorectal tumors and in cancer stem cell-enriched spheroids derived from CRC cells. CircZNF800 is shown to sponge a small set of three miRNAs to up-regulate three targeted downstream mRNAs. The multiple circZNF800-miRNA-mRNA axes form chains of regulatory events likely to be related to enhanced expression of cancer stem cell markers and associated cellular functions. Our findings highlight the critical roles of circZNF800 in promoting cancer stem cell properties and contributing to the CRC tumorigenesis process. Silencing the circZNF800 suppresses the cancer stem cell phenotype in vitro and impedes tumor growth in mice and may form the basis for developing strategies for therapeutic treatment of colorectal cancer.

Availability of data and materials

All analyzed data are included in this published article and its supplementary information files. The original data are available upon reasonable request to the corresponding authors.

Abbreviations

ACVR1C:

Activin A Receptor Type 1C

AKT:

AKT serine/threonine kinase

ALK7:

Activin A receptor type 1C

circRNA:

Circular RNA

circZNF800:

Circular RNA ZNF800

CRC:

Colorectal cancer

CRISPR:

Clustered Regularly Interspaced Short Palindromic Repeats

crRNA:

CRISPR RNA

crSC :

CRISPR RNA with scrambled sequence

Cy-5:

Cyanine-5

EGFP:

Enhanced green fluorescence protein

FACS:

Fluorescence Activated Cell Sorter

FISH:

in situ Hybridization

FZD3:

Frizzled class receptor 3

GPRC5A:

G protein-coupled receptor class C group 5 member A

ISC:

Intestinal stem cell

KD:

Knockdown

LGR5:

Leucine rich repeat containing G protein-coupled receptor 5

miRNA:

MicroRNA

NOTCH1:

Notch receptor 1

OE:

Over-expression

SOX9:

SRY-box transcription factor 9

UTR:

Untranslated region

VAPA:

VAMP associated protein A

WNT5A:

WNT family member 5A

ZNF800:

Zinc-finger protein 800

References

  1. Batlle E, Clevers H. Cancer stem cells revisited. Nat Med. 2017;23:1124–34.

    Article  CAS  PubMed  Google Scholar 

  2. Munro MJ, Wickremesekera SK, Peng L, Tan ST, Itinteang T. Cancer stem cells in colorectal cancer: a review. J Clin Pathol. 2018;71:110–6.

    Article  CAS  PubMed  Google Scholar 

  3. Kreso A, Dick JE. Evolution of the cancer stem cell model. Cell Stem Cell. 2014;14:275–91.

    Article  CAS  PubMed  Google Scholar 

  4. Zhou Y, Xia L, Wang H, Oyang L, Su M, Liu Q, Lin J, Tan S, Tian Y, Liao Q, et al. Cancer stem cells in progression of colorectal cancer. Oncotarget. 2018;9:33403–15.

    Article  PubMed  Google Scholar 

  5. van der Heijden M, Vermeulen L. Stem cells in homeostasis and cancer of the gut. Mol Cancer. 2019;18:66.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Barker N, Ridgway RA, van Es JH, van de Wetering M, Begthel H, van den Born M, Danenberg E, Clarke AR, Sansom OJ, Clevers H. Crypt stem cells as the cells-of-origin of intestinal cancer. Nature. 2009;457:608–11.

    Article  CAS  PubMed  Google Scholar 

  7. Huels DJ, Sansom OJ. Stem vs non-stem cell origin of colorectal cancer. Br J Cancer. 2015;113:1–5.

    Article  CAS  PubMed Central  Google Scholar 

  8. Agliano A, Calvo A, Box C. The challenge of targeting cancer stem cells to halt metastasis. Semin Cancer Biol. 2017;44:25–42.

    Article  CAS  PubMed  Google Scholar 

  9. Ishiguro T, Ohata H, Sato A, Yamawaki K, Enomoto T, Okamoto K. Tumor-derived spheroids: Relevance to cancer stem cells and clinical applications. Cancer Sci. 2017;108:283–9.

    Article  CAS  PubMed Central  Google Scholar 

  10. Fang D, Nguyen TK, Leishear K, Finko R, Kulp AN, Hotz S, Van Belle PA, Xu X, Elder DE, Herlyn M. A tumorigenic subpopulation with stem cell properties in melanomas. Cancer Res. 2005;65:9328–37.

    Article  CAS  PubMed  Google Scholar 

  11. Fekir K, Dubois-Pot-Schneider H, Desert R, Daniel Y, Glaise D, Rauch C, Morel F, Fromenty B, Musso O, Cabillic F, et al. Retrodifferentiation of human tumor hepatocytes to stem cells leads to metabolic reprogramming and chemoresistance. Cancer Res. 2019;79:1869–83.

    Article  CAS  PubMed  Google Scholar 

  12. Adeshakin FO, Adeshakin AO, Afolabi LO, Yan D, Zhang G, Wan X. Mechanisms for modulating anoikis resistance in cancer and the relevance of metabolic reprogramming. Front Oncol. 2021;11: 626577.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. De Angelis ML, Francescangeli F, Zeuner A, Baiocchi M. Colorectal cancer stem cells: An overview of evolving methods and concepts. Cancers (Basel). 2021;13:5910.

    Article  PubMed  Google Scholar 

  14. Huang C-J, Choo KB. Circular RNA- and microRNA-mediated post-transcriptional regulation of preadipocyte differentiation in adipogenesis: from expression profiling to signaling pathway. Intl J Mol Sci. 2023;24:4549.

    Article  CAS  Google Scholar 

  15. Zhou W-Y, Cai Z-R, Liu J, Wang D-S, Ju H-Q, Xu R-H. Circular RNA: metabolism, functions and interactions with proteins. Mol Cancer. 2020;19:172.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Panda AC. Circular RNAs act as miRNA sponges. Adv Exp Med Biol. 2018;1087:67–79.

    Article  CAS  PubMed  Google Scholar 

  17. Zang J, Lu D, Xu A. The interaction of circRNAs and RNA binding proteins: An important part of circRNA maintenance and function. J Neurosci Res. 2020;98:87–97.

    Article  CAS  PubMed  Google Scholar 

  18. Su M, Xiao Y, Ma J, Tang Y, Tian B, Zhang Y, Li X, Wu Z, Yang D, Zhou Y, et al. Circular RNAs in cancer: emerging functions in hallmarks, stemness, resistance and roles as potential biomarkers. Mol Cancer. 2019;18:90.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Rengganaten V, Huang CJ, Tsai PH, Wang ML, Yang YP, Lan YT, Fang WL, Soo S, Ong HT, Cheong SK, et al. Mapping a circular RNA-microRNA-mRNA-signaling regulatory axis that modulates stemness properties of cancer stem cell populations in colorectal cancer spheroid cells. Intl J Mol Sci. 2020;21:7864.

    Article  CAS  Google Scholar 

  20. Panda AC, Martindale JL, Gorospe M. Affinity pulldown of biotinylated RNA for detection of protein-RNA complexes. Bio-Protoc. 2016;6:e2062.

    PubMed  Google Scholar 

  21. Kramer MF. Stem-loop RT-qPCR for miRNAs. Curr Protoc Mol Biol Chapter. 2011;15:10.

    Google Scholar 

  22. Ford E, Ares M Jr. Synthesis of circular RNA in bacteria and yeast using RNA cyclase ribozymes derived from a group I intron of phage T4. Proc Natl Acad Sci USA. 1994;91:3117–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Meng C, Zhao X, Lao J. A modified immunofluorescence in situ hybridization method to detect long non-coding RNAs and proteins in frozen spinal cord sections. Exp Ther Med. 2018;15:4623–8.

    PubMed  PubMed Central  Google Scholar 

  24. Salzman J, Chen RE, Olsen MN, Wang PL, Brown PO. Cell-type specific features of circular RNA expression. PLoS Genet. 2013;9:e1003777.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Rybak-Wolf A, Stottmeister C, Glazar P, Jens M, Pino N, Giusti S, Hanan M, Behm M, Bartok O, Ashwal-Fluss R, et al. Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed. Mol Cell. 2015;58:870–85.

    Article  CAS  PubMed  Google Scholar 

  26. Maass PG, Glazar P, Memczak S, Dittmar G, Hollfinger I, Schreyer L, Sauer AV, Toka O, Aiuti A, Luft FC, et al. A map of human circular RNAs in clinically relevant tissues. J Mol Med. 2017;95:1179–89.

    Article  CAS  PubMed  Google Scholar 

  27. Roche KC, Gracz AD, Liu XF, Newton V, Akiyama H, Magness ST. SOX9 maintains reserve stem cells and preserves radioresistance in mouse small intestine. Gastroenterology. 2015;149:1553–63.

    Article  CAS  PubMed  Google Scholar 

  28. Nishimura R, Osako T, Okumura Y, Hayashi M, Toyozumi Y, Arima N. Ki-67 as a prognostic marker according to breast cancer subtype and a predictor of recurrence time in primary breast cancer. Exp Ther Med. 2010;1:747–54.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Wesselhoeft RA, Kowalski PS, Anderson DG. Engineering circular RNA for potent and stable translation in eukaryotic cells. Nat Commun. 2018;9:2629.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Zhang Y, Nguyen TM, Zhang XO, Wang L, Phan T, Clohessy JG, Pandolfi PP. Optimized RNA-targeting CRISPR/Cas13d technology outperforms shRNA in identifying functional circRNAs. Genome Biol. 2021;22:41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Li S, Li X, Xue W, Zhang L, Yang LZ, Cao SM, Lei YN, Liu CX, Guo SK, Shan L, et al. Screening for functional circular RNAs using the CRISPR-Cas13 system. Nat Methods. 2021;18:51–9.

    Article  PubMed  Google Scholar 

  32. Zhang J, Espinoza LA, Kinders RJ, Lawrence SM, Pfister TD, Zhou M, Veenstra TD, Thorgeirsson SS, Jessup JM. NANOG modulates stemness in human colorectal cancer. Oncogene. 2013;32:4397–405.

    Article  CAS  PubMed  Google Scholar 

  33. Amini S, Fathi F, Mobalegi J, Sofimajidpour H, Ghadimi T. The expressions of stem cell markers: Oct4, Nanog, Sox2, nucleostemin, Bmi, Zfx, Tcl1, Tbx3, Dppa4, and Esrrb in bladder, colon, and prostate cancer, and certain cancer cell lines. Anat Cell Biol. 2014;47:1–11.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Munro MJ, Wickremesekera SK, Peng L, Marsh RW, Itinteang T, Tan ST. Cancer stem cell subpopulations in primary colon adenocarcinoma. PLoS One. 2019;14: e0221963.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Gisina AM, Kim YS, Potashnikova DM, Tvorogova AV, Yarygin KN, Lupatov AY. Proliferative activity of colorectal cancer cells with different levels of CD133 expression. Bull Exp Biol Med. 2019;167:541–5.

    Article  CAS  Google Scholar 

  36. Marhaba R, Klingbeil P, Nuebel T, Nazarenko I, Buechler MW, Zoeller M. CD44 and EpCAM: cancer-initiating cell markers. Curr Mol Med. 2008;8:784–804.

    Article  CAS  PubMed  Google Scholar 

  37. Vo JN, Cieslik M, Zhang Y, Shukla S, Xiao L, Zhang Y, Wu YM, Dhanasekaran SM, Engelke CG, Cao X, et al. The landscape of circular RNA in cancer. Cell. 2019;176:869–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Lu HJ, Li J, Yang G, Yi CJ, Zhang D, Yu F, Ma Z. Circular RNAs in stem cells: from basic research to clinical implications. Biosci Rep. 2022;42:BSR20212510.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Yu CY, Li TC, Wu YY, Yeh CH, Chiang W, Chuang CY, Kuo HC. The circular RNA circBIRC6 participates in the molecular circuitry controlling human pluripotency. Nat Commun. 2017;8:1149.

    Article  PubMed  Google Scholar 

  40. Zhao X, Zhong Q, Cheng X, Wang S, Wu R, Leng X, Shao L. miR-449c-5p availability is antagonized by circ-NOTCH1 for MYC-induced NOTCH1 upregulation as well as tumor metastasis and stemness in gastric cancer. J Cell Biochem. 2020;121:4052–63.

    Article  CAS  PubMed  Google Scholar 

  41. Gu C, Zhou N, Wang Z, Li G, Kou Y, Yu S, Feng Y, Chen L, Yang J, Tian F. circGprc5a promoted bladder oncogenesis and metastasis through Gprc5a-targeting peptide. Mol Ther Nucleic Acids. 2018;13:633–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Civelek M, Wu Y, Pan C, Raulerson CK, Ko A, He A, Tilford C, Saleem NK, Stancakova A, Scott LJ, et al. genetic regulation of adipose gene expression and cardio-metabolic traits. Am J Hum Genet. 2017;100:428–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Zhuo E, Cai C, Liu W, Li K, Zhao W. Downregulated microRNA-140-5p expression regulates apoptosis, migration and invasion of lung cancer cells by targeting zinc finger protein 800. Oncol Lett. 2020;20:390.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Wang J, Li T, Wang B. Circ-UBAP2 functions as sponges of miR-1205 and miR-382 to promote glioma progression by modulating STC1 expression. Cancer Med. 2021;10:1815–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Wan X, Guo D, Zhu Q, Qu R. microRNA-382 suppresses the progression of pancreatic cancer through the PI3K/Akt signaling pathway by inhibition of Anxa3. Am J Physiol Gastrointest Liver Physiol. 2020;319:G309–22.

    Article  CAS  PubMed  Google Scholar 

  46. Yao H, Xia D, Li ZL, Ren L, Wang MM, Chen WS, Hu ZC, Yi GP, Xu L. MiR-382 functions as tumor suppressor and chemosensitizer in colorectal cancer. Biosci Rep. 2019;39:BSR20180441.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Xu M, Jin H, Xu CX, Sun B, Song ZG, Bi WZ, Wang Y. miR-382 inhibits osteosarcoma metastasis and relapse by targeting Y box-binding protein 1. Mol Ther. 2015;23:89–98.

    Article  CAS  PubMed  Google Scholar 

  48. Borcherding N, Kusner D, Kolb R, Xie Q, Li W, Yuan F, Velez G, Askeland R, Weigel RJ, Zhang W. Paracrine WNT5A signaling inhibits expansion of tumor-initiating cells. Cancer Res. 2015;75:1972–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Kurayoshi M, Oue N, Yamamoto H, Kishida M, Inoue A, Asahara T, Yasui W, Kikuchi A. Expression of Wnt-5a is correlated with aggressiveness of gastric cancer by stimulating cell migration and invasion. Cancer Res. 2006;66:10439–48.

    Article  CAS  PubMed  Google Scholar 

  50. Zhao Y, Dai Q, Fu X, Chen Q, Tang Y, Gao X, Zhou Q. CircVAPA exerts oncogenic property in non-small cell lung cancer by the miR-876-5p/WNT5A axis. J Gene Med. 2021;23:e3325.

    Article  CAS  PubMed  Google Scholar 

  51. Qin L, Yin YT, Zheng FJ, Peng LX, Yang CF, Bao YN, Liang YY, Li XJ, Xiang YQ, Sun R, et al. WNT5A promotes stemness characteristics in nasopharyngeal carcinoma cells leading to metastasis and tumorigenesis. Oncotarget. 2015;6:10239–52.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Chen X, Jia C, Jia C, Jin X, Gu X. MicroRNA-374a inhibits aggressive tumor biological behavior in bladder carcinoma by suppressing Wnt/beta-catenin signaling. Cell Physiol Biochem. 2018;48:815–26.

    Article  CAS  PubMed  Google Scholar 

  53. Wu S, Wang H, Pan Y, Yang X, Wu D. miR-140-3p enhances cisplatin sensitivity and attenuates stem cell-like properties through repressing Wnt/beta-catenin signaling in lung adenocarcinoma cells. Exp Ther Med. 2020;20:1664–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Kalhori MR, Irani S, Soleimani M, Arefian E, Kouhkan F. The effect of miR-579 on the PI3K/AKT pathway in human glioblastoma PTEN mutant cell lines. J Cell Biochem. 2019;120:16760–74.

    Article  CAS  PubMed  Google Scholar 

  55. Principe M, Chanal M, Karam V, Wierinckx A, Mikaélian I, Gadet R, Auger C, Raverot V, Jouanneau E, Vasiljevic A, et al. ALK7 expression in prolactinoma is associated with reduced prolactin and increased proliferation. Endocr Relat Cancer. 2018;25:795–806.

    Article  CAS  PubMed  Google Scholar 

  56. Asnaghi L, White DT, Key N, Choi J, Mahale A, Alkatan H, Edward DP, Elkhamary SM, Al-Mesfer S, Maktabi A, et al. ACVR1C/SMAD2 signaling promotes invasion and growth in retinoblastoma. Oncogene. 2019;38:2056–75.

    Article  CAS  PubMed  Google Scholar 

  57. Xia L, Wu L, Bao J, Li Q, Chen X, Xia H, Xia R. Circular RNA circ-CBFB promotes proliferation and inhibits apoptosis in chronic lymphocytic leukemia through regulating miR-607/FZD3/Wnt/beta-catenin pathway. Biochem Biophys Res Commun. 2018;503:385–90.

    Article  CAS  PubMed  Google Scholar 

  58. Pamudurti NR, Patop IL, Krishnamoorthy A, Ashwal-Fluss R, Bartok O, Kadener S. An in vivo strategy for knockdown of circular RNAs. Cell Discov. 2020;6:52.

    Article  CAS  PubMed Central  Google Scholar 

  59. Wells DJ. Systemic AAV gene therapy close to clinical trials for several neuromuscular diseases. Mol Ther. 2017;25:834–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Professor Manuel Ares, University of California, Santa Cruz, USA, for the permutated group I self-catalytic of td gene intron plasmid.

Funding

This study was funded by the Ministry of Science and Technology, Taiwan (NSTC111-2320-B-A49-028-MY3 and NSTC112-2321-B-A49-007), Taipei Veterans General Hospital (V112E-002–4, V112C-135), the IBMS CRC Research Program of the Institute of Biomedical Sciences, Academia Sinica, Taiwan (IBMS-CRC109-P04), the Veterans Affairs Council (110VACS-003, 110VACS-007), VGHUST Joint Research Program (VGHUST112-G1-5-1), Center for Intelligent Drug System and Smart Bio-device (IDS2B) from the Featured Areas Research Center Program within the framework of the Higher Education Sprout Project by the Ministry of Education, Ministry of Economy Affairs (112-EC-17-A-20-S6-026), and Ministry of Education, Higher Education SPROUT project for Cancer Progression Research Center (111W31101) and Cancer and Immunology Research Center (112W31101), awarded to SHC.

Author information

Authors and Affiliations

Authors

Contributions

SHC, VR, KBC, CJH & HTO designed the experiments and discussed the data. VR conducted the experiments. MLW, YC, PHT consulted on experimental design and data analysis. YTL provided the CRC clinical samples. KBC, VR & CJH wrote the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Shih-Hwa Chiou or Kong Bung Choo.

Ethics declarations

Ethics approval and consent to participate

Primary CRC and adjacent normal paired tissues were obtained from patients who underwent colectomy in Taipei Veterans General Hospital, Taiwan, and informed consent was obtained from all participants, study was in accordance with the International Ethical Guidelines for Biomedical Research Involving Human Subjects and was approved by the Medical Ethical Committee of Taipei Veterans General Hospital, Taiwan (IRB-2019–08-009B). The mouse work was conducted according to the Guidelines for Laboratory Animals in the Taipei Veterans General Hospital, and the institutional ethical review board of Taipei Veterans General Hospital approved the study (IACUC 2012–079), and was in accordance with ARRIVE guidelines."

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Suppl. file 1: Table S1.

List of PCR primers and probes used.

Additional file 2: Suppl. file 2: Table S2.

Expression profiling of circZNF800 isoforms in cancer and non-cancer cells and tissues.

Additional file 3: Suppl. file 3: Fig. S1.

Uncropped and unedited RNA gel electrophoresis for in vitro transcribed and circularized circZNF800.

Additional file 4: Suppl. file 4: Fig. S2.

CircZNF800-miRNA & miRNA-mRNA seed sequence alignments.

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rengganaten, V., Huang, CJ., Wang, ML. et al. Circular RNA ZNF800 (hsa_circ_0082096) regulates cancer stem cell properties and tumor growth in colorectal cancer. BMC Cancer 23, 1088 (2023). https://doi.org/10.1186/s12885-023-11571-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-11571-1

Keywords