Circular RNA circIPO11 drives self-renewal of liver cancer initiating cells via Hedgehog signaling
Molecular Cancer volume 20, Article number: 132 (2021)
Hepatocellular carcinoma (HCC) is one of the most intractable tumors in the world due to its high rate of recurrence and heterogeneity. Liver cancer initiating cells also called cancer stem cells (CSCs) play a critical role in resistance against typical therapy and high tumor-initiating potential. However, the role of the novel circular RNA (circRNA) circIPO11 in the maintenance of liver cancer initiating cells remains elusive.
CircRNAs highly conserved in humans and mice were identified from 3 primary HCC samples by circRNA array. The expression and function of circIPO11 were further evaluated by Northern blot, limiting dilution xenograft analysis, chromatin isolation by RNA purification-PCR assay (ChIRP) and HCC patient-derived tumor cells (PDC) models. CircIpo11 knockout (KO) mice were generated by a CRISPR/Cas9 technology.
CircIPO11 is highly expressed in HCC tumor tissues and liver CSCs. CircIPO11 is required for the self-renewal maintenance of liver CSCs to initiate HCC development. Mechanistically, circIPO11 recruits TOP1 to GLI1 promoter to trigger its transcription, leading to the activation of Hedgehog signaling. Moreover, GLI1 is also highly expressed in HCC tumor tissues and liver CSCs, and TOP1 expression levels positively correlate with the metastasis, recurrence and survival of HCC patients. Additionally, circIPO11 knockout in mice suppresses the progression of chemically induced liver cancer development.
Our findings reveal that circIPO11 drives the self-renewal of liver CSCs and promotes the propagation of HCC via activating Hedgehog signaling pathway. Antisense oligonucleotides (ASOs) against circIPO11 combined with TOP1 inhibitor camptothecin (CPT) exert synergistic antitumor effect. Therefore, circIPO11 and the Hedgehog signaling pathway may provide new potential targets for the treatment of HCC patients.
Hepatocellular carcinoma (HCC), a major type of primary liver cancer, is the second leading cause of cancer-related deaths globally . Because of its increasing morbidity and mortality, liver cancer has become a growing global health-care problem . But the effective treatments for liver cancer are limited. Traditional surgical resection, radiotherapy and chemotherapy are prone to disease recurrence, resulting in poor prognosis . New treatments come from an improved understanding of HCC tumorigenesis. Two main biological characteristics of HCC are high recurrence and heterogeneity . Studies have shown that a group of stem/progenitor cells known as cancer initiating cells or cancer stem cells (CSCs) can contribute to organization hierarchy and lead to heterogeneity . Unlike bulk tumor cells, these CSCs display the ability to self-renew, differentiate, and form new tumors, accounting for resistance to traditional treatments and the high recurrence rate of HCC . Recently, plenty of cell surface markers have been identified to separate liver CSCs, such as CD13, CD133, CD90 and EpCAM and so on . However, the mechanism of maintaining liver CSC self-renewal remains elusive.
CSCs are similar to embryonic or adult tissue stem cells in that they rely on activation of highly conserved stemness signaling pathways to maintain their self-renewal and regeneration . Accumulating evidence suggests that alteration of stemness signaling pathways, including Hedgehog (Hh), Wnt and Notch signaling, can promote tumor initiation and progression . Of note, aberrant activation of these pathways in stem cells may evoke uncontrolled cell proliferation and abnormal differentiation, leading to tumor initiation. Moreover, their reactivation can induce tumor reprogramming, resulting in the onset of CSC phenotypes. These stemness pathways also contribute to the hepatocarcinogenesis and recurrence . We previously isolated a rare subset from HCC cells using two surface markers, CD13 and CD133, and defined CD13+CD133+ cells as liver CSCs, CD13−CD133− cells as non-CSCs [11, 12]. Moreover, we also confirmed that several long noncoding RNAs (lncRNAs) can regulate the self-renewal of liver CSCs via Wnt and BMP signaling [13, 14].
Circular RNAs (circRNAs), a class of single-stranded covalent closed RNAs, are generated by backsplicing derived from precursor mRNA (pre-mRNA). As a new class of non-coding RNAs, some circRNAs are highly conserved across species with cell-specific and tissue-specific characteristics. Recent studies have shown that some specific circRNAs that accumulate in the brain are conserved from human to fruit fly . Due to difficultly degraded by RNaseR and lack of polyadenylated (poly A) tails compared with linear RNA, circRNAs are rarely detected in next-generation RNA sequencing (RNA-seq) profiling that is usually enriched in poly(A) + RNA. Until now, most of the well-characterized circRNAs have been identified and established database such as circBase (http://circrna.org/). CircRNAs functions in extensive biological processes have been elucidated by diverse mechanisms. Emerging studies have illustrated that some circRNAs are involved in neuropsychiatric disorder, innate immune responses, hematopoiesis and cell proliferation [16,17,18]. Interestingly, circ-CDYL, specifically up-regulated in the early stage of HCC, can be used as a potential biomarker for early surveillance of liver cancer . However, it is largely unknown whether circRNAs participate in the regulation of liver CSC self-renewal. The ability of circRNAs to regulate gene expression makes them potential drug targets. Targeting circLONP2 via anti-sense oligonucleotides (ASOs) showed excellent therapeutic effects on colorectal carcinoma . However, whether circRNAs can be used as drug candidates in liver cancer is still worth exploring.
DNA topoisomerases have an essential role in chromatin dynamics, transcription, replication, DNA damage repair and genomic stability via introducing transient DNA breaks for relaxing supercoiled DNA . For example, topoisomerase 1 (TOP1) binds DNA substrates to form a TOP1 cleavage complex (TOP1cc), which causes transient DNA nicks to relax supercoiled DNAs. Camptothecin (CPT) and its derivatives inhibit TOP1 by binding at the enzyme-DNA interface leading to blockade of religation of TOP1cc . Previous studies showed that inhibition of topoisomerases can increase chemosensitivity of HCC cells . However, how topoisomerases function in liver CSCs remains unclear. Here we identified a conserved circRNA circIPO11 (originated from IPO11 gene transcript, circbase symbol hsa_circ_0007915) that is highly expressed in tumor tissues and liver CSCs. CircIPO11 associates with TOP1 to initiate GLI1 transcription and activates Hh signaling for self-renewal maintenance of liver CSCs. Moreover, combination administration of ASOs against circIPO11 and TOP1 inhibitor CPT has effective therapeutic effect on HCC tumors.
Materials and methods
Antibodies and reagents
Anti-histone H3 (catalog 4499) antibody was purchased from Cell Signaling Technology. Anti-EEA1 (catalog sc-6415) antibody was from Santa Cruz Biotechnology. Anti-TOP1 (catalog ab109374) antibody was from Abcam. Anti-β-actin (catalog A1978), anti-Flag (catalog F1804) antibodies were from Sigma-Aldrich. Anti-GLI1 (catalog 66,905–1-Ig), anti-IPO11 (catalog 14,403–1-AP), anti-SMO (catalog 20,787–1-AP), anti-CD13 (catalog 14,553–1-AP) and anti-Nanog (catalog 14,295–1-AP) antibodies were purchased from Proteintech Group Inc. Anti-PTCH1 (catalog LS-C114391) antibody was from LifeSpan BioSciences. HRP-conjugated secondary antibodies were from Sungene Biotech. FITC-conjugated CD13 antibody (catalog 11–0138) was obtained from eBioscience. PE-conjugated CD133 (catalog 130–098-826) was from Miltenyi Biotec. Alexa405, Alexa488 and Alexa594-conjugated secondary antibodies were obtained from Invitrogen. bFGF (catalog GF446-50UG) was obtained from Millipore. EGF (catalog E5036-200UG) and DAPI (catalog 28,718–90-3) were from Sigma-Aldrich. N2 supplement (catalog 17,502–048) and B27 (catalog 17,504–044) were obtained from Invitrogen. The Dual Luciferase Reporter Gene Assay Kit (catalog RG027) was purchased from Beyotime. D-Luciferin (catalog 40902ES01) was obtained from Yeasen BioTechnologies. Biotin RNA Labeling Mix (catalog 11,685,597,910) and T7 RNA polymerase (catalog 10,881,767,001) were from Roche. The LightShift™ Chemiluminescent RNA EMSA kit (catalog 20,158) and Chemiluminescent Nucleic Acid Detection Module (catalog 89,880) were from Thermo Scientific. Ultra-low attachment plates (catalog 3471) were purchased from Corning 70μm company.
Cell lines and HCC samples
Human HCC cell lines Huh7 and PLC/PRF/5 (PLC) were provided by Dr. Zeguang Han (Shanghai Jiao Tong University School of Medicine, Shanghai, China). All these cell lines were maintained in DMEM medium supplemented with 10% FBS, 100 μg/ml penicillin G, and 100 U/ml streptomycin. Human HCC samples were obtained from consenting patients at the Department of Hepatobiliary Surgery, PLA General Hospital (Beijing, China). HCC samples were treated within 2 h after resection. Necrotic tissue was removed, and tumor was cut into 1 mm3 with scissors. Then digested with collagenase IV for 45 min (0.05% collagenase IV, 0.05% proteinase, 0.01% DNase), shaken every 15 min. A 70 μm cell filter was used and centrifuged at 50 g for 1 min. Supernatants were collected and centrifuged at 150 g for 8 min. And HCC cells were enriched at the bottom. After red cell elimination, HCC primary cells were obtained and used for sphere cultivation. Oncosphere cells were cultured in DMEM/F12 serum-free medium supplemented with 20 ng/ml bFGF, 20 ng/ml EGF, N2 and B27.
CircIpo11−/− mice were generated using CRISPR/Cas9 approaches . The intronic sequence between exon 3 and exon 4 of Ipo11 gene was targeted by sgRNAs (Table S2). After confirming cleavage efficiency of these sgRNAs in 293 T cells, in vitro transcribed mRNAs containing the sgRNA/Cas9 sequences were injected into zygotes through microinjection. Approximate 200 ~ 250 zygotes in C57BL/6 N background were implanted into the uterus of surrogate mice, from which F0 mice were generated. Genotypes of F0 mice were subsequently verified by PCR and DNA sequencing. Identification of primers were listed in Table S3. WT allele had a PCR length of about 222 bp and deficient allele had a PCR length of about 336 bp. F0 mice were crossed to generate circIpo11−/− mice. BALB/c nude mice were from Beijing Vital River Laboratory Animal Technology. Immunodeficient B-NDG mice were from Biocytogen. Mice were sacrificed when they developed tumors larger than 15 mm in diameter or skin ulceration. Mice were fed under specific pathogen-free (SPF) grade with approval by the Institutional Animal Care and Use Committees at the Institute of Biophysics, Chinese Academy of Sciences. For DEN-induced liver tumors model, Mice were injected intraperitoneally with 25 mg/kg of diethylnitrosamine (DEN, Sigma, 73,861) at 2 weeks old. Livers were collected and analyzed after indicated times of DEN injection.
CRISPR/Cas9 knockout system
GLI1, TOP1, and GLI1 promoter-deletion cells were established using CRISPR/Cas9 approaches provided by Zhang’s lab . All sgRNAs were designed by online CRISPR Design Tool (http://crispr.mit.edu/). In brief, for GLI1 or TOP1 targeting, sgRNAs were designed to target the fourth exon of GLI1 (exon 4) or TOP1 (exon 3), and cloned into LentiCRISPRv2 (puro, catalog 52,961, Addgene). We used LentiCRISPRv2, pVSVg (catalog 8454, Addgene), and psPAX2 (catalog 12,260, Addgene) plasmids to produce CRISPR–Cas9 lentivirus. For GLI1-deleted and TOP1-deleted cells, single sgRNA was used. For GLI1 promoter-deleted cells, a pair of sgRNAs targeting sequences flanking the binding sequence (− 2050 bp ~ − 1850 bp) was used. sgRNA sequences targeting GLI1, TOP1, and GLI1 promoter were listed in Table S2. We generated lentivirus in 293 T with lipo3000 for 2 days. After the supernatant was filtered with 70μm strainer and mixed with equal volume fresh DMEM (10%FBS), HCC cells were infected for 12 h, followed by puromycin selection and monoclonalization. The monoclines were propagated and confirmed by DNA electrophoresis and sequencing. The correctly targeted clones were obtained and used for subsequent experiments. For rescue assay, circIPO11 or GLI1–overexpressing lentivirus was co-infected, followed by GFP sorting. Finally, GLI1 KO CSCs and non-CSCs were isolated by FACS for further research.
Overexpression of circIPO11, ΔcircIPO11 and GLI1
For circIPO11 overexpression, the exonic sequences of circIPO11 (Exon4 and 5 of IPO11 gene) and complementary sequences were needed. The intronic sequences including complementary sequences were obtained from genomic DNA. Then the full sequences were cloned into the pBPLV-GFP vector. For ΔcircIPO11 overexpression, the TOP1-binding region (31-94 nt) was deleted. For GLI1 overexpression, pBPLV-GFP-GLI1 plasmids were constructed. Then these plasmids and their package plasmids (pBPLV: VSVG: pLp1: pLp2 = 5: 2.8: 4.2: 2μg) were transfected into 293 T cells for virus generation. After supernatants was filtered with 70 μm strainer and mixed with equal volume fresh DMEM (10%FBS), HCC cells were infected for 12 h, followed by GFP sorting.
Immunohistochemical staining was performed as previously described . In brief, paraffin-embedded tumor tissue sections were deparaffinized in xylene and rehydrated in graded alcohol (100, 100, 95, 85, and 70% alcohols). The slides were processed for antigen retrieval in Tris-EDTA buffer (10 mM, pH 8.0), 100 °C for 15 min, and then cooled down naturally. Treatment with 3% hydrogen peroxide (H2O2) for 10 min blocked the activity of endogenous peroxidase. After blocking with 10% donkey serum for 30 min, the sections were incubated overnight at 4 °C with the primary antibodies. After washing with PBS, HRP-conjugated secondary antibodies were incubated at RT for 1 h. Subsequently detection was performed using DAB, counterstained with hematoxylin, dehydrated and mounted.
IVIS in vivo imaging
IVIS in vivo imaging was performed as described . In brief, anaesthetized B-NDG mice and open the abdominal cavity to expose the liver. 1 × 106 luciferase transfected Huh7 cells were in situ injected into the liver with insulin needle. The wound was sutured and coated with antibiotics, and fed with acidified water (PH3.0). After 2 weeks, D-fluorescein potassium salt was injected intraperitoneally of the mice. After 5 min, the mice were anesthetized with isoflurane. Then the luciferase activity was detected in the IVIS-100 imaging system. The collected images are analyzed using Living Image 4.3 software (PerkinElmer).
Chromatin immunoprecipitation (ChIP) assay
ChIP assay was performed according to the protocol (Uptate Biotechnology, Inc.). In brief, oncosphere cells were treated with 1% formaldehyde for 10 min at 37 °C, and then lysed by SDS buffer for 10 min on ice, followed by ultrasonic (30% of maximum power value) to get 200 ~ 1000 bp DNA fragments. The samples were precleared with salmon sperm DNA/protein agarose beads for 0.5 h in rotor at 4 °C, and then incubated with the anti-TOP1 antibodies overnight. The enrichments after elution were analyzed by qPCR. Primer sequences are shown in the Table S3.
Chromatin isolation by RNA purification-PCR assay (ChIRP)
ChIRP assay were performed as described . In brief, cross-linked with 1% glutaraldehyde, oncosphere cells were cracked by SDS buffer for 10 min on ice. Subsequently chromatins were sonicated into fragments 200 ~ 1000 bp in length. Antisense probes for circIPO11 were designed by Biosearch Probe Designer (targeting junction sequence of circIPO11). Two sets of biotin-labeled probes (Probeset-circ and Probeset-Ctrl) were added into samples for incubation at 4 °C overnight. Combinative chromatin fragments were enriched by Streptavidin C1 magnetic beads before being purified for qPCR examination.
Fixed by 4% paraformaldehyde (PFA) for 20 min, HCC cells were permeated by 1% triton X-100 for 30 min. After blocking with 10% donkey serum for 30 min, primary antibodies were added and incubated 2 h at 37 °C. After washing with PBS, fluorescence-conjugated secondary antibodies were added and incubated at 37 °C for 1 h. After sealing, confocal microscopy (Nikon A1R+) was performed for observation.
Measurement of circIPO11 copy numbers
The linearized plasmid pcDNA3-circIPO11 was continuously diluted to generate a standard curve of circIPO11 by qRT-PCR. The copy number of the diluted pcDNA3-circIPO11 was analyzed via DNA/RNA Copy Number Calculator from website (http://endmemo.com/bio/dnacopynum.php). To measure the circIPO11 copy number in spheres and CSCs, total RNA extracted from 3 × 105 cells of each HCC line was reverse transcribed into cDNAs, followed by qPCR analysis. The copy number could be measured from the standard curve.
Cell proliferation assay
2000 ~ 5000 indicated cells were transplanted into a well of 6-well plates and cultured in a cell incubator until the cells grew to 75%. After washing with PBS, the cells were fixed with 4% PFA for 10 min, and incubated with 0.1% crystal violet for 30 min at 37 °C. Plates were washed gently with distilled water, and placed in the oven to air-dry. After photographed, 33% acetic acid was added to each well for decolorization, and then the absorbance value was measured at 570 nm after sufficient shaking.
Sphere formation assay
Oncospheres were cultured as previous described . Indicated 2000 HCC cells or 5000 HCC primary cells were grown in low attachment 6-well plates and cultured in sphere formation medium (DMEM/F12 supplemented with 20 ng/ml EGF, 20 ng/ml bFGF, B27 supplement, N2 supplement). Numbers of spheres were counted under an optical microscope after 1 ~ 2 weeks and photographs were taken. For non-sphere cell separation, we collected cell culture containing non-sphere cells and sphere cells, and let stand for 10 min. Pellets were sphere cells, and supernatants were non-sphere cells. The sphere cells were collected for subsequent experiments. Supernatants were then collected by centrifugation at 2000 g for 10 min. Pellets were non-sphere cells and used directly for subsequent assays.
Generation of PDC models
Mouse PDC models were generated as described . For tumor-initiating assays, gradient dilutions (10, 102, 103, and 104) of indicated Huh7 cells were injected with 100 μl Matrigel subcutaneously into the BALB/c nude mice. The percentages of tumor-free mice were calculated 3 months later. Tumorigenic cell frequencies in circIPO11-depleted and overexpressed tumor cells were analyzed with a limiting dilution assay (http://bioinf.wehi.edu.au/software/elda/). Eight parallel groups were used per sample. For tumor-formation assays, xenografts were produced by injecting 1 × 106 oncosphere cells subcutaneously into BALB/c nude mice. Tumor sizes were measured every 4 days. Five mice were used per sample. Twenty-five milligram per kilogram ASOs were injected around the tumor every 2 days (16, 18, 20, 22, and 24). Mice were killed at 35th day after HCC primary cells injection.
RNA pulldown and mass spectrometry assay
Biotin-labeled RNAs (circIPO11 junction sequences, antisense and IPO11 intron sequences) were transcribed with the T7 RNA polymerase and biotin RNA labeling mix. For RNA pulldown assay , probes were incubated with HCC oncosphere lysates overnight at 4 °C. Then the streptavidin-conjugated agarose beads were added for enrichment. After washed with RIPA buffer, biotin-enriched proteins were separated by SDS–PAGE, and visualized by silver staining. Differential bands in SDS-PAGE gels were collected for mass spectrometry (LTQ Orbitrap XL, Thermo).
Isolation of liver CSCs was performed as previously described . Briefly, HCC sample CSCs were labeled with FITC-conjugated CD13 and PE-conjugated CD133 antibodies and isolated by flow cytometry. To analyze survival and proliferation rates of CSCs, Ki-67(17–5699-41, eBioscience), AnnexinV-APC and 7-AAD apoptosis detection kit (AO2001-11A-H, Sungene Biotech) were used according to the manufacturer’s instructions. FlowJo 7.6 software was used for subsequent analysis.
CircRNAs sequencing and microarray assay
For circRNAs sequencing, RNA was isolated with Trizol method from HCC primary tumor and peri-tumor tissues, and purified with RNease Mini kit. Then circRNA were sequenced using human circRNA Array (8x15K, Arraystar) was performed. For identification of circRNA downstream target genes, total RNA was isolated with Trizol from circIPO11-depleted or control oncosphere cells. Then NimbleGen microarray analysis was performed.
RNA EMSA assay
EMSA experiments were performed using a LightShift Chemiluminescent RNA EMSA Kit (Thermo Scientific). Recombinant TOP1 protein was obtained by purification (Constructed into plasmids pGEX-6P-1 and 3xFlag). RNA fragments of circIPO11 were transcribed with T7 and labeled with biotin according to standard protocols. Probes and recombinant proteins were incubated in binding buffer and mobility shift assay was performed using native gel electrophoresis.
DNase I sensibility assay
Cell nuclei were extracted according to the protocol from a nuclei-isolation kit (Sigma-Aldrich). Then cell nuclei were digested with 1 U DNase I for 5 min at 37 °C. After stopping the digestion, total DNA was extracted and followed by qRT-PCR.
Total RNA was extracted from indicated samples using Trizol, then subjected to electrophoresis on 2% denaturing agarose gel with 1% formaldehyde for 1.5 h. Samples were transferred to positively charged NC membranes with 20 × SSC buffer. After UV cross-linking (265 nm ultraviolet with energy of 200,000 μJ/cm2) and prehybridization, membranes were incubated with biotin-labeled probes at 65 °C for 16 ~ 20 h. After washed with washing buffer, biotin signals were detected with Chemiluminescent Nucleic Acid Detection Module according to the manufacturer’s instructions. For detecting circRNAs only, junction sequences were used for probes.
Psoralen photobinding assay
Psoralen photobinding assay was performed as described. Briefly, the CPT was added into oncosphere cells, and then treated with psoralen at 37 °C for 8 min. Cross-linked with UV, the oncosphere cells were collected and lysed at 55 °C for 5 h. The DNA was extracted, and followed by ultrasonic to get DNA fragments. The cross-linked and uncross-linked segments were separated by PCR, and then analyzed by qPCR.
Data were analyzed with an unpaired Student’s t-test using the GraphPad Prism 8 software and Excel 2016. p < 0.05 was considered significant (*p < 0.05; **p < 0.01; ***p < 0.001); NS, non-significant. For survival analysis, the Kaplan-Meier survival analysis was used. All flow cytometry data were analyzed with FlowJo 7.6.
CircIPO11 is highly expressed in HCC tumors and liver CSCs
To identify the function circRNAs in the tumorigenesis of HCC, we conducted circRNA transcriptome analysis of human HCC tumor and peri-tumor tissues with human circRNA array (Arraystar). Of note, only a small fraction (10–20%) was highly conserved between humans and mice . We further screened out 46 differentially expressed circRNAs as well as highly conserved in humans and mice (Fig. 1A), including 34 upregulated and 12 downregulated circRNAs. Top six upregulated circRNAs in HCC tumors were listed on the right (Fig. 1A). Biological characteristics of these circRNAs were further validated by sequencing (Fig. S1A), PCR (Fig. S1B), RNase R digestion (Fig. S1C) and actinomycin D treatment (Fig. S1D). In order to explore the roles of these circRNAs in liver CSCs, we depleted these circRNAs in HCC samples and analyzed their sphere formation capacity. Among these circRNAs, circIPO11 depletion significantly decreased the ability of oncosphere formation (Fig. 1B). CircIPO11 was highly expressed in liver cancer tissues by a transcriptomic sequencing , which is consistent with our results. CircIPO11, located on human chromosome 5, consists of the fourth and fifth exons of IPO11 gene (Fig. 1C). Combining circRNA Array and circBase, IPO11 gene transcribed nine circRNA transcripts due to variable cyclization (Fig. S1E). Compared with other cycled transcripts, circIPO11 transcript was abundant in tumors (Fig. S1F). In addition, circIPO11 complementary upstream and downstream sequences were highly conserved in various species (Fig. S1G-J). Of note, circIPO11 was more substantially expressed in tumor than peri-tumor tissues (Fig. 1D), oncospheres than non-sphere cells (Fig. 1E), and CSCs than non-CSCs (Fig. 1F). These results were further validated by Northern blot. In addition, we found that copy numbers of circIPO11 in liver oncospheres and CSCs were about 200 per cell (Fig. 1G). Moreover, circIPO11 was mainly localized in the nucleus via in situ hybridization of HCC samples (Fig. 1H), nuclear-cytoplasmic separation assay (Fig. 1I) and Immunofluorescence staining of human tumor tissues (Fig. 1J). Altogether, circIPO11, a conserved circRNA, is highly expressed in HCC tumors and liver CSCs that participates in the maintenance of liver CSC self-renewal.
CircIPO11 depletion impairs liver CSC self-renewal maintenance and inhibits HCC propagation
In order to further explore the role of circIPO11 in liver CSCs, we targeted the head-to-tail junction by lentivirus-mediated short hairpin RNA (shRNA) to obtain circIPO11 depleted HCC cells and confirmed by qRT-PCR (Fig. S2A). Of note, circIPO11 depletion did not affect intracellular levels of its parental gene IPO11 (Fig. S2B, C). We observed that circIPO11 depletion significantly reduced primary, secondary, and tertiary sphere formation capacity. In addition, overexpression of circIPO11 could rescue the number of oncospheres (Fig. 2A, B). We found that the proportion of CD13+CD133+ cells decreased significantly when circIPO11 was depleted (Fig. S2D, E).
We next wanted to determine the function of circIPO11 in vivo. Extremely limiting dilutions of different cancer cell populations, and then transplants serial tumor into immunodeficient mice to measure their capacity to reform secondary tumors is the gold standard to assess CSC potential . Through limiting dilution xenograft analysis, circIPO11 depletion remarkably impaired CSC self-renewal and reduced tumor initiation (Fig. 2C). By limiting dilution analysis, circIPO11 knockdown significantly reduced self-renewal capacity of liver CSCs (Fig. S2F). In addition, we subcutaneously injected 1 × 106 circIPO11 depletion or control cells into BALB/c nude mice. CircIPO11 depletion cells remarkably reduced tumor propagation (Fig. 2D). To assess the role of circIPO11 on orthotopic liver tumor development, circIPO11 depletion and control cells were transfected with plasmids containing luciferase (Luc). We observed that circIPO11 depletion dramatically reduced tumor growth in situ (Fig. 2E). To further determine the role of circIPO11 in liver CSCs, we isolated liver CSCs (CD13+CD133+) and non-CSCs (CD13−CD133−) from Huh7-Luc cells, followed by depletion of circIPO11 and transplantation into livers of immunodeficient mice. We observed that circIPO11 knockdown dramatically suppressed xenograft tumor growth in CSCs but not in non-CSCs (Fig. S2G). Inhibition of circIPO11 displayed poor proliferation capacity via clone formation assay (Fig. 2F). Then, we used two effective antisense oligonucleotides (ASOs) against circIPO11 junction to investigate its therapeutic roles. We found ASOs treatment of HCC cell line could significantly inhibit circIPO11 expression (Fig. 2G). We observed that the two ASOs significantly inhibited tumor growth in BALB/c nude mice (Fig. 2H). These data indicate that circIPO11 promotes the self-renewal maintenance of liver CSCs.
Given that circIPO11 was highly homologous in humans and mice, mouse ortholog circIpo11 also has two exons originated from its parental gene Ipo11 gene on chr 13 (Fig. S3A). Mouse circIpo11 was further validated by sequencing (Fig. S3B). Since intronic complementary sequences flanking exons are necessary for formation of exonic circRNAs, the exonic circRNAs do not produce if either of intronic complementary sequences is missing . Moreover, intronic complementary sequences verified by mini-gene assay played an important role in formation of circIpo11 (Fig. S3C, D). Based on this, we next generated circIpo11 knockout (KO) mice by deleting the upstream reverse complementary sequence via a CRISPR/Cas9 technology (Fig. S3E). CircIpo11 was completely deleted in mice by northern bolt (Fig. 2I) and the upstream complementary pairing region was knocked out and confirmed by PCR (Fig. S3F). Meanwhile, the inability of the divergent primer to amplify the circIpo11 fragment further indicated that circIpo11 was successfully KO in mice by PCR (Fig. S3G). Similarly, circIpo11 KO did not affect the expression of Ipo11 gene (Fig. S3H). Of note, circIpo11 KO mice had no apparent phenotypes, and displayed normal liver development. Diethylnitrosamine (DEN) has been commonly used to chemically induced HCC. We injected mice intraperitoneally injected with DEN. At 8th and 12th month after DEN treatment, sacrificed mice for macroscopic examination showed that circIpo11−/− liver tumor foci were dramatically reduced in size compared with circIpo11+/+ littermate control mice (Fig. 2J). In addition, livers from DEN-treated circIpo11−/− mice displayed much smaller and fewer nodular hyperplastic areas in histological sections than those of DEN-treated WT mice (Fig. 2K-L). The proportion of liver weight vs. body weight is another important factor that reflects tumor malignance . Ratios of liver weights vs. body weights of circIpo11−/− mice were much lower than WT mice (Fig. 2M). Collectively, circIpo11 KO impairs the self-renewal maintenance of liver CSCs and suppresses tumor propagation as well.
Overexpression of circIPO11 promotes liver CSC self-renewal and tumorigenic capacity
We next overexpressed circIPO11 in HCC cells (Fig. 3A, Fig. S2H). Of note, circIPO11 overexpression did not affect the expression level of its parental gene IPO11 (Fig. S2I, J). CircIPO11 overexpression dramatically increased the capacity of oncosphere formation (Fig. 3B). Serial passage assays illustrated that circIPO11 overexpression improved the capacity of liver CSC self-renewal (Fig. 3C). Though limiting dilution, subcutaneous tumor formation, and in vivo imaging assay, we found that circIPO11 overexpression could significantly increase the proportion of tumor-bearing mice (Fig. 3D) and promote HCC progression of in situ (Fig. 3E, F). Meanwhile, circIPO11 overexpression significantly increased self-renewal capacity of CSCs by limiting dilution analysis (Fig. S2K). In parallel, circIPO11 overexpression enhanced proliferation of HCC cells (Fig. 3G). Taken together, circIPO11 overexpression promotes self-renewal and tumorigenic capacity of liver CSCs.
CircIPO11 associates with TOP1 in liver CSCs
We then investigated the molecular mechanism by which circIPO11 regulated liver CSCs. Through circBank and NCBI website, none of the open reading frame (ORF) was potentially existed, suggesting that circIPO11 is a noncoding RNA. Moreover, there was also little miRNA binding sites in circIPO11, indicating that circIPO11 does not serve as microRNA sponges. In addition, there was no change of its parental gene IPO11 by knockdown or overexpressing circIPO11 (Fig. S2B, C, Fig. S2G, H). Interestingly, based on circBank, circIPO11 was predicted to have protein-binding potential. On the other hand, circIPO11 harbored some stem-loop structures through prediction analysis by the RNAfold Web Services (Fig. S4A), suggesting that circIPO11 may associate with proteins to perform its function. RNA pull down assay using biotin-labeled probes was performed for searching potential associated proteins of circIPO11 in liver oncospheres. After mass spectrometry analysis, we identified topoisomerases 1 (TOP1), relaxing supercoiled DNA during gene replication and transcription, as a candidate binding protein (Fig. 4A and Fig. S4B). The interaction of circIPO11 with TOP1 was further validated by Western blot (Fig. 4B). In addition, we evaluated whether circIPO11 interacted with TOP1 by structure prediction. CircIPO11 secondary structure was submitted to RNA Composer to generate the 3D structure (PDB: 1a36). NPDock was then used to calculate the molecular docking between circIPO11 and TOP1, which indicated that circIPO11 perfectly docked TOP1 (Fig. S4A). Through RNA immunoprecipitation (RIP) assay, TOP1 was able to precipitate circIPO11 in oncosphere lysates, rather than other variable cyclizations derived from IPO11 (Fig. 4C). Co-localization of circIPO11 with TOP1 was further confirmed by immunofluorescence staining in oncosphere cells (Fig. 4D). And we found they were mainly located in the nucleus. The NTD domain of TOP1 protein was necessary to bind circIPO11 via domain mapping assays (Fig. 4E), suggesting circIPO11 directly binds to TOP1. In addition, the segment 31–94 nt of circIPO11 was sufficient to bind TOP1 through fragment mapping (Fig. 4F). Finally, the connection of circIPO11 segment (31–94 nt) with TOP1 was further verified by RNA electrical mobility shift assay (EMSA) (Fig. 4G).
We next overexpressed circIPO11 lacking the TOP1-binding region (ΔcircIPO11) in HCC cells. The capacity of oncosphere formation (Fig. 4H) and cell proliferation (Fig. S4C) was almost abrogated in ΔcircIPO11 overexpressed HCC cells, suggesting that the interaction of circIPO11 with TOP1 is necessary for the self-renewal of liver CSCs. Collectively, we conclude that circIPO11 interacts with TOP1 in liver CSCs.
CircIPO11 recruits TOP1 onto GLI1 promoter to activate its expression
To further identify target genes of circIPO11, we performed transcriptome microarray analysis using circIPO11 depleted and control oncosphere cells (Fig. 5A, Fig. S5A). Top 13 downregulated transcription factors (TFs) were selected to further verify their expression in circIPO11 depleted oncosphere cells by qRT-PCR (Fig. 5B). We then obtained these 13 TFs silenced HCC cell lines via shRNAs and performed sphere formation assays. Among them, knockdown of GLI1 dramatically affected sphere formation ability in Huh7 cell lines (Fig. 5C). GLI1 (GLI family zinc finger 1), acts exclusively as an activator of Hh signaling, since it lacks a repressor domain compared to GLI2 and GLI3 . The strength and downstream gene levels of Hh pathway are regulated by the ratio of GLI1 activators and inhibitors forms. Hh/GLI signaling mediates various cellular activities by activating distinct target gene sets, especially in the field of stem cells , including hair follicle stem cells, cancer stem cells, and follicle stem cell. We then performed chromatin immunoprecipitation (ChIP) assays in HCC oncospheres with anti-TOP1 antibody. We analyzed a 3 kb region upstream of the GLI1 gene transcription start site (TSS). We found that TOP1 was enriched at the − 2050 to − 1850 bp region of the GLI1 promoter (Fig. 5D). In addition, dual luciferase reporter assay also revealed TOP1 bound to the same locus of GLI1 promoter. (Figs. S5B, C).
To further explore the role of circIPO11 in regulating GLI1 expression, we performed ChIRP-PCR assay. We found that circIPO11 transcripts were mainly distributed on the − 2050 to − 1850 bp region of the GLI1 promoter as same as TOP1 did (Fig. 5E). Moreover, dual luciferase reporter assays further verified circIPO11 depletion dramatically decreased luciferase activity of the − 2050 to − 1850 bp region, illustrating circIPO11 could recruit TOP1 to GLI1 promoter (Fig. 5H). Furthermore, circIPO11 depletion significantly reduced TOP1 enrichment on the GLI1 promoter (Fig. 5I). In addition, circIPO11 overexpression remarkedly increased TOP1 enrichment on the GLI1 promoter (Fig. S4D). Camptothecin (CPT), a classic inhibitor of TOP1 , was used to treat cells with different doses and detected the expression of GLI1 by qRT-PCR and Western blot. We found that when TOP1 was inhibited by 10 μM CPT, and the intracellular level of GLI1 was completely repressed (Fig. 5F). Similarly, through psoralen photobinding assay , we observed that there was a negative superhelix structure on the region of GLI1 promoter -2 kb (Fig. 5G). While circIPO11 and TOP1 bound to this region, the negative superhelix structure would be released to initiate the transcription of GLI1. We next tested chromatin accessibility of GLI1 promoter through DNase I sensibility assay. We found that circIPO11 depletion dramatically wrinkled chromatin accessibility of the GLI1 promoter (Fig. 5J). In contrast, circIPO11 overexpression substantially opened the chromatin accessibility of the GLI1 promoter (Fig. S5E). Next, we constructed TOP1 deleted HCC cell lines by CRISPR/Cas9 technology and identified them by Western blot (Fig. S5F). And then dual luciferase reporter assays further verified TOP1 depletion dramatically decreased the luciferase activity of the − 2050 to − 1850 bp region (Fig. S5G). In general, these data suggest that circIPO11 recruits TOP1 to the GLI1 promoter to initiate its expression.
GLI1 initiates Hh signaling activation for liver CSC self-renewal
To further test the downstream signaling pathway of circIPO11 in liver CSCs, integrative gene set enrichment analysis (GSEA) of transcriptome microarray data were conducted. We noticed that the target genes of Hh signaling were remarkably suppressed with circIPO11 depletion, suggesting Hh signaling was involved in the regulation of liver CSCs (Fig. 6A). We recently showed that Hh signaling pathway plays an important role in the regulation of liver CSCs . We then constructed GLI1 knockout HCC cells using CRISPR/Cas9 technology and identified them by sequencing (Fig. S6A, B). GLI1 protein was completely abrogated by Western blot detection (Fig. S6C). Of note, GLI1 deletion dramatically changed the expression levels of Hh signaling target genes and downstream stemness gene (Fig. 6B, C). As expected, GLI1 KO cells could reduce oncosphere formation (Fig. 6D) and cell proliferation capacity (Fig. S6D). Since CPT is a classical inhibitor against TOP1, we then used CPT to treat GLI1 KO cells for oncosphere formation assays. We noticed that treatment of CPT could not impair GLI1 deletion-induced oncosphere formation in vitro (Fig. 6D). Unexpectedly, we found an interesting phenomenon that tumors were significantly smaller through subcutaneous tumor formation experiment after treatment with CPT (Fig. 6E). We speculated that the inhibition of TOP1 by CPT might also have damage on normal cells, thus causing the storm of immune factors and increasing the killing effect on tumors, finally leading to the reduction of tumors in vivo. To further examine the role of GLI1 in liver CSCs, we isolated liver CSCs and non-CSCs from Huh7 cells, followed by depletion of GLI1 and subcutaneous transplantation into BALB/c nude mice. We found that GLI1 KO dramatically inhibited xenograft tumor growth in CSCs (Fig. S6E). In addition, we noticed that proliferation-related transcription factors and stemness factors were significantly upregulated in circIPO11 overexpressed tumor cells (Fig. S6F). Finally, we depleted GLI1 in HCC cells with versus without circIPO11 overexpression, followed oncosphere formation assays. We found that GLI1 depletion in circIPO11 overexpressed HCC cells could not promote oncosphere formation capacity compared with control HCC cells.
Furthermore, the TOP1 binding region of GLI1 promoter (GLI1-P-KO) was deleted in HCC cells using CRISPER/Cas9 system (Fig. S6G). We observed that GLI1-P-KO decreased sphere formation and cell colony formation, whereas GLI1 overexpression rescued this phenomenon (Fig. 6F, G). However, circIPO11 overexpression did not improve sphere formation in GLI1-P-KO cells (Fig. S6H), indicating the TOP1 binding region of GLI1 promoter was required for circIPO11 function. Taken together, circIPO11-mediated GLI1 upregulation enhances Hh signaling activation that initiates the self-renewal maintenance of liver CSCs and tumor propagation.
GLI1 and TOP1 expression are positively related to HCC severity
To further examine the function of GLI1 in liver CSCs, we detected the expression levels of GLI1 in HCC tumors and peri-tumor tissues, oncespheres and non-sphere cells, as well as liver CSCs and non-CSCs. We found that GLI1 was highly expressed in tumor tissues, liver CSCs and oncosphere cells (Fig. 7A-C). High expression of GLI1 in HCC samples was further validated by TCGA database (Fig. 7D). The expression of GLI1 was positively correlated with the poor prognosis of HCC patients (Fig. 7E). Of note, GLI1, mainly localized in the nucleus, was highly expressed in HCC samples via in situ hybridization (Fig. 7F). These data suggest that GLI1 expression levels are positively correlated with HCC severity. CircIPO11 depletion decreased the expression of GLI1 (Fig. 7G). Moreover, GLI1 overexpression in circIPO11 depleted cells remarkably enhanced GLI1 expression and sphere formation (Fig. 7G, H). To further determine the correlation of circIPO11, TOP1 and GLI1, we examined their expression levels among 29 HCC samples. We observed that circIPO11 expression levels were positively related to expression levels of GLI1 and TOP1 (Fig. S7A, B). Moreover, GLI1 expression levels were positively correlated with TOP1 expression levels (Fig. S7C). In addition, circIPO11 knockdown and GLI1 KO dramatically suppressed proliferation of CSCs, but not affected cell death (Fig. S7D-G).
We further analyzed the expression of TOP1 with Wang’s cohort (GSE14520). We noticed that TOP1 was also highly expressed in HCC tumor tissues (Fig. 7I). In addition, the expression of TOP1 was significantly correlated with HCC patients’ survival, metastasis and recurrence (Fig. 7J-M), indicating that TOP1 expression levels are related to the severity of HCC. Taken together, GLI1 and TOP1 expression are positively correlated with HCC severity.
Cancer stem cells (CSCs) have been considered to be the source of tumor initiation, drug resistance, metastasis, and recurrence . Some circRNAs play critical roles in many cellular processes, particularly in cancer. circGprc5a, is highly expressed in bladder CSCs and regulates the self-renewal of CSCs through a coding peptide. However, how circRNAs regulate liver CSC stemness remains enclusive. Here we identified a conserved circRNA named circIPO11 that is significantly increased in liver CSCs and necessary for the self-renewal maintenance of liver CSCs. Mechanistically, circIPO11 recruits TOP1 onto GLI1 promoter to trigger its expression, leading to activation of Hh signaling. Moreover, GLI1 is also highly expressed in liver CSCs. TOP1 expression is positively correlated with the survival, metastasis and recurrence of HCC patients. Targeting circIPO11 and TOP1 has synergistic antitumor effects on HCC patient-derived tumor cells (PDC) models.
Compared with traditional linear RNAs, circRNAs are more stable and hard to degrade owing to their covalently closed circular structure, which is not affected by RNA exonuclease . Because circRNAs have a very long half-life and tissue-specific features, their functions are more complicated than expected . Due to the moderate conservation of most circRNAs, it is difficult to use animal models to research their physiological roles and test their therapeutic effects under pathological effects. In this study, we therefore focus on identifying highly conserved circRNAs and exploring their roles in tumorigenesis. To identify evolutionarily conserved circRNAs, we innovatively performed homology screening from differential circRNAs compared with tumor and peri-tumor tissues. We found that circIPO11 was highly conserved across various species. Its mouse ortholog is circIpo11 from murine Ipo11 gene including 2 exons. To explore the function of circRNA in vivo, we generated circIpo11 knockout (KO) mice by deleting the upstream reverse complementary sequence using CRISPR/Cas9 technology. Intriguingly, circIpo11−/− liver tumor foci were dramatically reduced in size than circIpo11+/+ littermates, suggesting that circIPO11 play a critical role in tumorigenesis and progression of liver cancer. Up to date, the relationship between circRNA and its linear mRNA is still controversial. Some studies showed that circRNA competes with its parental RNA  and others have suggested that circRNAs can promote the transcription of their parental genes in cis . In addition, some circRNAs, highly specific in brain, are independent of their linear isoforms generated from their parental genes [15, 41]. In this study, we showed that circIPO11 expression is not significant correlation with its parental gene IPO11 in liver CSCs, suggesting a critical role of circIPO11 in the regulation of liver CSCs. CircIPO11 is an exon circular RNA consisting of the fourth and fifth exons of IPO11 gene. The function of IPO11 is controversial, which has been reported as a tumor suppressor . However, this observation was most challenged by IPO11 as a required factor for β-catenin-mediated transcription in colorectal cancer . We found that IPO11 expression is comparable between HCC tumor and peri-tumor tissues and has no significant correlation with HCC patients’ prognosis. Additionally, circIPO11 depletion does not affect the expression of its parental gene IPO11 and other IPO11-derived circRNAs. We thus conclude that circIPO11 regulates the self-renewal of liver CSCs in trans.
Reorganization of chromatin and DNA topological pressure are intrinsic to transcription. Accumulated data on the dynamics of DNA and DNA-binding elements during transcription provide new insights into the mechanism of gene expression . DNA topoisomerases induce either single or double stranded DNA breaks to regulate DNA topological status during replication, transcription, recombination, and chromatin remodeling . On the one hand, TOP1 acts as an elongation factor by relaxing topological pressure in the process of transcription, such as transcription of the Hsp70 genes in Drosophila . On the other hand, TOP1 initiates transcription by binding to the promoter of target gene that impacts its expression . In addition, there exists a dynamic superhelix structure spread ~ 1.5 kilobases upstream of the transcription starting site of the target genes, as the binding site of TOP1 . Recent studies have shown that inhibition of TOP1 in HCC can improve the sensitivity to chemical drugs [23, 48], but the relationship between TOP1 and liver CSCs is still unknown. Here we showed that circIPO11 can interact with TOP1 and GLI1 promoter to produce relaxing supercoiled DNA, leading to the transcription of GLI1 gene.
The Hh signaling pathway is involved in numerous biological processes including embryonic development, tissue homeostasis, regeneration and healing. Aberrant Hh signaling may cause various human malignancies and CSCs [9, 49]. Hh pathway is inhibited by the unliganded form of receptor patched1 (PTCH1) in the absence of Hh ligands, including Sonic (Shh), Desert (Dhh), and Indian (Ihh). Following ligands engagement, PTCH1 promotes the activation and ciliary entry of the transmembrane protein Smoothened (Smo), which activates GLI transcription factors to initiate tumorigenesis. GLI1, as a downstream effector of Hh signaling, can promote the transcription of some crucial stemness genes such as NANOG. In addition, GLI1 expression is epigenetically regulated by chromatin remodeling complexes. Some studies have shown that ZNF521 can interact with the NuRD complex onto GLI promoter to activate its expression . However, how Hh signaling is activated and regulates the progression of HCC is still rarely explored. Our findings reveal that Hh signaling mediated by activation of GLI1 plays a critical role in the self-renewal of liver CSCs and liver tumorigenesis.
Some studies have reported circRNAs as novel potential biomarkers for cancer. For example, exosomal circNRIP1 serves as a biomarker for gastric cancer diagnosis . Therapies targeting cancer-associated circRNAs are also under investigation. ASO therapy has been applied in neural diseases in non-human primates and human clinical trials through intrathecal administration, without serious adverse reactions . ASO therapy has been shown to be safe in human clinical trials by intrathecal administration . In this study, we revealed that ASOs against circIPO11 can dramatically reduce tumor size, which indicates that targeting circIPO11 by ASOs could be potential therapeutic strategy for HCC patients. However, the efficient delivery of these oligonucleotides to tumors and the eliciting of an immune response still require further investigation.
In summary, our findings reveal that circIPO11 initiates the self-renewal of liver CSCs and promotes the propagation of HCC. Mechanistically, we first discovered that circIPO11 recruits TOP1 to GLI1 promoter to trigger its transcription, leading to the activation of Hedgehog signaling. Our findings provide a novel potential target for the treatment of HCC patients.
Availability of data and materials
All the data obtained and/or analyzed during the current study were available from the corresponding authors on reasonable request.
Cancer stem cells
Long noncoding RNAs
Open reading frame
Electrical mobility shift assay
GLI family zinc finger 1
Gene set enrichment analysis
Transcription start site
RNA immunoprecipitation assay
Chromatin immunoprecipitation assay
Chromatin isolation by RNA purification-PCR assay
Patient-derived tumor cells models
- Poly A:
Fluorescent in situ hybridization
Hematoxylin and Eosin
Topoisomerase cleavage complex
Marquardt JU, Andersen JB, Thorgeirsson SS. Functional and genetic deconstruction of the cellular origin in liver cancer. Nat Rev Cancer. 2015;15:653–67. https://doi.org/10.1038/nrc4017.
Torre LA, et al. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65:87–108. https://doi.org/10.3322/caac.21262.
Ji J, et al. MicroRNA expression, survival, and response to interferon in liver cancer. N Engl J Med. 2009;361:1437–47. https://doi.org/10.1056/NEJMoa0901282.
Visvader JE. Cells of origin in cancer. Nature. 2011;469:314–22. https://doi.org/10.1038/nature09781.
Meacham CE, Morrison SJ. Tumour heterogeneity and cancer cell plasticity. Nature. 2013;501:328–37. https://doi.org/10.1038/nature12624.
Kreso A, Dick JE. Evolution of the cancer stem cell model. Cell Stem Cell. 2014;14:275–91. https://doi.org/10.1016/j.stem.2014.02.006.
Sun JH, Luo Q, Liu LL, Song GB. Liver cancer stem cell markers: progression and therapeutic implications. World J Gastroenterol. 2016;22:3547–57. https://doi.org/10.3748/wjg.v22.i13.3547.
Kaur G, Sharma P, Dogra N, Singh S. Eradicating Cancer stem cells: concepts, issues, and challenges. Curr Treat Options in Oncol. 2018;19:20. https://doi.org/10.1007/s11864-018-0533-1.
Takebe N, et al. Targeting notch, Hedgehog, and Wnt pathways in cancer stem cells: clinical update. Nat Rev Clin Oncol. 2015;12:445–64. https://doi.org/10.1038/nrclinonc.2015.61.
Chen J, et al. The microtubule-associated protein PRC1 promotes early recurrence of hepatocellular carcinoma in association with the Wnt/beta-catenin signalling pathway. Gut. 2016;65:1522–34. https://doi.org/10.1136/gutjnl-2015-310625.
Wang Y, et al. The long noncoding RNA lncTCF7 promotes self-renewal of human liver cancer stem cells through activation of Wnt signaling. Cell Stem Cell. 2015;16:413–25. https://doi.org/10.1016/j.stem.2015.03.003.
Wu J, et al. The long non-coding RNA LncHDAC2 drives the self-renewal of liver cancer stem cells via activation of Hedgehog signaling. J Hepatol. 2019;70:918–29. https://doi.org/10.1016/j.jhep.2018.12.015.
Zhu P, et al. Lnc-beta-Catm elicits EZH2-dependent beta-catenin stabilization and sustains liver CSC self-renewal. Nat Struct Mol Biol. 2016;23:631–9. https://doi.org/10.1038/nsmb.3235.
Wang Y, et al. LncRNA HAND2-AS1 promotes liver cancer stem cell self-renewal via BMP signaling. The EMBO journal. 2019;38:e101110. https://doi.org/10.15252/embj.2018101110.
Rybak-Wolf A, et al. Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed. Mol Cell. 2015;58:870–85. https://doi.org/10.1016/j.molcel.2015.03.027.
Piwecka M, et al. Loss of a mammalian circular RNA locus causes miRNA deregulation and affects brain function. Science. 2017;357. https://doi.org/10.1126/science.aam8526.
Liu CX, et al. Structure and Degradation of Circular RNAs Regulate PKR Activation in Innate Immunity. Cell. 2019;177:865–80 e821. https://doi.org/10.1016/j.cell.2019.03.046.
Xia P, et al. A Circular RNA Protects Dormant Hematopoietic Stem Cells from DNA Sensor cGAS-Mediated Exhaustion. Immunity. 2018;48:688–701 e687. https://doi.org/10.1016/j.immuni.2018.03.016.
Wei Y, et al. A noncoding regulatory RNAs network driven by Circ-CDYL acts specifically in the early stages hepatocellular carcinoma. Hepatology. 2020;71:130–47. https://doi.org/10.1002/hep.30795.
Han K, et al. CircLONP2 enhances colorectal carcinoma invasion and metastasis through modulating the maturation and exosomal dissemination of microRNA-17. Mol Cancer. 2020;19:60. https://doi.org/10.1186/s12943-020-01184-8.
Pommier Y, Sun Y, Huang SN, Nitiss JL. Roles of eukaryotic topoisomerases in transcription, replication and genomic stability. Nat Rev Mol Cell Biol. 2016;17:703–21. https://doi.org/10.1038/nrm.2016.111.
Pommier Y, et al. Interaction of an Alkylating Camptothecin Derivative with a DNA-Base at Topoisomerase I-DNA Cleavage Sites. P Natl Acad Sci USA. 1995;92:8861–5. https://doi.org/10.1073/pnas.92.19.8861.
Wang N, et al. MiR-23a-mediated inhibition of topoisomerase 1 expression potentiates cell response to etoposide in human hepatocellular carcinoma. Mol Cancer. 2013;12:119. https://doi.org/10.1186/1476-4598-12-119.
Wang Y, et al. Long noncoding RNA lncHand2 promotes liver repopulation via c-met signaling. J Hepatol. 2018;69:861–72. https://doi.org/10.1016/j.jhep.2018.03.029.
Zhu P, et al. LncGata6 maintains stemness of intestinal stem cells and promotes intestinal tumorigenesis. Nat Cell Biol. 2018;20:1134–44. https://doi.org/10.1038/s41556-018-0194-0.
Zhu P, et al. LncBRM initiates YAP1 signalling activation to drive self-renewal of liver cancer stem cells. Nat Commun. 2016;7:13608. https://doi.org/10.1038/ncomms13608.
Dong R, Ma XK, Chen LL, Yang L. Increased complexity of circRNA expression during species evolution. RNA Biol. 2017;14:1064–74. https://doi.org/10.1080/15476286.2016.1269999.
Ren S, Xin Z, Xu Y, Xu J, Wang G. Construction and analysis of circular RNA molecular regulatory networks in liver cancer. Cell Cycle. 2017;16:2204–11. https://doi.org/10.1080/15384101.2017.1346754.
Boumahdi S, et al. SOX2 controls tumour initiation and cancer stem-cell functions in squamous-cell carcinoma. Nature. 2014;511:246–50. https://doi.org/10.1038/nature13305.
Patop IL, Wust S, Kadener S. Past, present, and future of circRNAs. EMBO J. 2019;38:e100836. https://doi.org/10.15252/embj.2018100836.
Kondylis V, et al. NEMO prevents Steatohepatitis and hepatocellular carcinoma by inhibiting RIPK1 kinase activity-mediated hepatocyte apoptosis. Cancer Cell. 2015;28:830. https://doi.org/10.1016/j.ccell.2015.11.007.
Hui CC, Angers S. Gli proteins in development and disease. Annu Rev Cell Dev Biol. 2011;27:513–37. https://doi.org/10.1146/annurev-cellbio-092910-154048.
Hsu YC, Li L, Fuchs E. Transit-amplifying cells orchestrate stem cell activity and tissue regeneration. Cell. 2014;157:935–49. https://doi.org/10.1016/j.cell.2014.02.057.
Rialdi A, et al. Topoisomerase 1 inhibition suppresses inflammatory genes and protects from death by inflammation. Science. 2016;352:aad7993. https://doi.org/10.1126/science.aad7993.
Kouzine F, et al. Transcription-dependent dynamic supercoiling is a short-range genomic force. Nat Struct Mol Biol. 2013;20:396–403. https://doi.org/10.1038/nsmb.2517.
Gupta PB, Chaffer CL, Weinberg RA. Cancer stem cells: mirage or reality? Nat Med. 2009;15:1010–2. https://doi.org/10.1038/nm0909-1010.
Zhang XO, et al. Complementary sequence-mediated exon circularization. Cell. 2014;159:134–47. https://doi.org/10.1016/j.cell.2014.09.001.
Szabo L, Salzman J. Detecting circular RNAs: bioinformatic and experimental challenges. Nat Rev Genet. 2016;17:679–92. https://doi.org/10.1038/nrg.2016.114.
Xu XL, et al. CircRNA inhibits DNA damage repair by interacting with host gene. Mol Cancer. 2020;19:128. https://doi.org/10.1186/s12943-020-01246-x.
Li Z, et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. 2015;22:256–64. https://doi.org/10.1038/nsmb.2959.
Chen LL. The expanding regulatory mechanisms and cellular functions of circular RNAs. Nat Rev Mol Cell Bio. 2020;21:475–90. https://doi.org/10.1038/s41580-020-0243-y.
Chen M, et al. The nuclear transport receptor Importin-11 is a tumor suppressor that maintains PTEN protein. J Cell Biol. 2017;216:641–56. https://doi.org/10.1083/jcb.201604025.
Mis M, et al. IPO11 mediates betacatenin nuclear import in a subset of colorectal cancers. J Cell Biol. 2020;219. https://doi.org/10.1083/jcb.201903017.
Sainsbury S, Niesser J, Cramer P. Structure and function of the initially transcribing RNA polymerase II-TFIIB complex. Nature. 2013;493:437–40. https://doi.org/10.1038/nature11715.
Wang JC. Cellular roles of DNA topoisomerases: a molecular perspective. Nat Rev Mol Cell Bio. 2002;3:430–40. https://doi.org/10.1038/nrm831.
Zobeck KL, Buckley MS, Zipfel WR, Lis JT. Recruitment timing and dynamics of transcription factors at the Hsp70 loci in living cells. Mol Cell. 2010;40:965–75. https://doi.org/10.1016/j.molcel.2010.11.022.
Onishi Y, Kawano Y. Rhythmic binding of topoisomerase I impacts on the transcription of Bmal1 and circadian period. Nucleic Acids Res. 2012;40:9482–92. https://doi.org/10.1093/nar/gks779.
Zhang P, et al. MicroRNA-139 suppresses hepatocellular carcinoma cell proliferation and migration by directly targeting topoisomerase I. Oncol Lett. 2019;17:1903–13. https://doi.org/10.3892/ol.2018.9746.
Cazet AS, et al. Targeting stromal remodeling and cancer stem cell plasticity overcomes chemoresistance in triple negative breast cancer. Nat Commun. 2018;9:2897. https://doi.org/10.1038/s41467-018-05220-6.
Scicchitano S, et al. The stem cell-associated transcription co-factor, ZNF521, interacts with GLI1 and GLI2 and enhances the activity of the sonic hedgehog pathway. Cell Death Dis. 2019;10:715. https://doi.org/10.1038/s41419-019-1946-x.
Zhang X, et al. Circular RNA circNRIP1 acts as a microRNA-149-5p sponge to promote gastric cancer progression via the AKT1/mTOR pathway. Mol Cancer. 2019;18:20. https://doi.org/10.1186/s12943-018-0935-5.
Kordasiewicz HB, et al. Sustained therapeutic reversal of Huntington's disease by transient repression of huntingtin synthesis. Neuron. 2012;74:1031–44. https://doi.org/10.1016/j.neuron.2012.05.009.
Miller TM, et al. An antisense oligonucleotide against SOD1 delivered intrathecally for patients with SOD1 familial amyotrophic lateral sclerosis: a phase 1, randomised, first-in-man study. Lancet Neurol. 2013;12:435–42. https://doi.org/10.1016/S1474-4422(13)70061-9.
We thank Drs. Xiang Ding, Junying Jia, Yihui Xu, Xiang Shi and Shu Meng for technical support. We also thank Jing Li (Cnkingbio Company Ltd., Beijing, China) for technical support.
This work was supported by the National Natural Science Foundation of China (82130088, 81921003, 82022058, 92042302, 31930036, 91940305, 81872413, 81772646, 32070533); National Key R&D Program of China (2020YFA0803501, 2019YFA0508501), Strategic Priority Research Programs of the Chinese Academy of Sciences (XDB19030203), and Biological Resources Program of Chinese Academy of Sciences.
Ethics approval and consent to participate
All experiment involving human tumor tissues and mice were approved by the institutional committee of Institute of Biophysics, Chinese Academy of Sciences (SYXK2018–28).
Consent for publication
All authors give consent for the publication of manuscript in Molecular Cancer.
The authors have declared that no conflict of interest exists.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Characterization of upregulated circIPO11 in HCC tumor tissues. (A) Validation of circRNAs by DNA sequencing. PCR products with divergent primers were sequenced. Divergent primers were listed in Table S3. (B) Using complementary DNA (cDNA) and genomic DNA (gDNA) as templates, circIPO11 was amplified with divergent (grey arrowheads) and convergent primers (black arrowheads). Gapdh was used as a positive control. (C) Total RNAs extracted from HCC oncospheres were digested with or without 2 U/μg RNase R for 1 h at 37 °C, followed by qPCR analysis, IPO11 and ACTB were used as positive controls. Data are shown as means ± SD. (D) HCC cell lines were treated with 2 μg/ml actinomycin D for 14 h, and whole RNAs were extracted for qPCR analysis. IPO11 and ACTB were used as positive controls. Data are shown as means ± SD. (E) Schematic representation of nine circRNAs derived from IPO11 transcript due to variable cyclizations. (F) Expression levels of the variable cycled circRNAs were examined by qRT-PCR. CircIPO11 was the only circular RNA that highly expressed in tumor, oncospheres and liver CSCs. Data are shown as means ± SD. (G) Sequence conservation analysis of circIPO11 in vertebrates from zebrafish to human. The scores listed in the table were identified by the WATER algorithm. (H, I) Conservation analysis of upstream (H) and downstream (I) complementary intron sequence of circIPO11 in vertebrates. (J) Match information of complementary introns of circIPO11 in vertebrates. U, upstream complementary intron sequences; D, downstream complementary intron sequences. **P < 0.01; *** P < 0.001 by two-tailed Student’s t test. Data are representative of at least three independent experiments. Figure S2. CircIPO11 is required for the self-renewal of liver CSCs. (A) CircIPO11 expression was detected by qRT-PCR in HCC cell lines and HCC primary tumor cells. Data are shown as means ± SD. (B, C) Expression levels of IPO11 in circIPO11-depleted HCC cells were detected by Western blot (B) and qRT-PCR (C). β-actin was used as a loading control. Data are shown as means ± SD. (D) CD13 and CD133 subpopulations in HCC primary tumor cells were detected by flow cytometry. (E) CD13+CD133+ subpopulations were sorted from circIPO11-depleted and control HCC cell lines, followed by flow cytometry analysis. (F) Tumorigenic cell frequencies in circIPO11-depleted and empty vector control (shVec) cells were analyzed with a limiting dilution assay (http://bioinf.wehi.edu.au/software/elda/). Data are shown as means and 95% confidence intervals (n = 8). (G) Orthotopic liver tumor imaging of CSCs and non-CSCs with and without circIPO11 transduced with Luc vectors. Representative images are shown (left panel), and statistical results are shown as means ± SD (right panel). n = 3 for each group. (H) CircIPO11 was overexpressed in HCC cell lines and HCC primary tumor cells. Data are shown as means ± SD. (I, J) Expression of circIPO11 did not affect the expression of its parental gene IPO11 by qRT-PCR (I) and Western blot (J). β-actin was used as a loading control. Data are shown as means ± SD. (K) Tumorigenic cell frequencies in circIPO11-overexpressing and control (oeVec) cells were determined with limiting dilution assay (http://bioinf.wehi.edu.au/software/elda/). Data are shown as means and 95% confidence intervals (n = 8). *** P < 0.001 by two-tailed Student’s t test. Data are representative of at least three independent experiments. Figure S3. Generation of circIPO11 knockout mice. (A) Schematic representation of mouse circIPO11. Convergent primers for linear Ipo11 were denoted in black arrowheads, and divergent primers agaisnt circIPO11 were denoted in red arrowheads. E1, exon 1. (B) Validation of circIPO11 by DNA sequencing in mice. PCR products with divergent primers were sequenced. (C, D) Schematic representation of mini-gene assay was shown (C). circ△up: missing upstream intronic complementary sequences flanking exon 4. circ△down: missing downstream intronic complementary sequences flanking exon 5. Identification of complementary intron sequences were required for circIPO11 formation, followed by qRT-PCR (D). (E) Schematic representation for circIpo11 KO mice. The upstream complementary sequences (#1) were deleted using CRISPR/Cas9 system. Corresponding sgRNAs were shown in Table S2. (F) CircIpo11−/− mice were verified by agarose gel electrophoresis. Genome DNA of circIpo11+/+ and circIpo11−/− mice was used for PCR template. WT allele had a PCR length of about 222 bp and deficient allele had a PCR length of about 336 bp. (G) Using complementary DNA (cDNA) and genomic DNA (gDNA) as templates, circIpo11 was detected in WT or circIpo11−/− mice with divergent (grey arrowheads) and convergent primers (black arrowheads). (H) Ipo11 expression levels were detected by Western blot in circIpo11+/+ and circIpo11−/− livers. β-actin was used as a loading control. *P < 0.05; **P < 0.01; *** P < 0.001 by two-tailed Student’s t test. Data are representative of at least three independent experiments. Figure S4. TOP1 interacts with circIPO11. (A) CircIPO11 was predicted to interact with TOP1 via NPdock website. CircIPO11 secondary structure was predicted according to minimum free energy (MFE). Red color indicates strong confidence of the prediction. 3D structure of TOP1 is derived from PDB database. (B) MS profiles of TOP1, corresponding peptide sequences were listed on the top of diagram. (C) CircIPO11 and ΔcircIPO11 were overexpressed in HCC primary cells and HCC cell lines, followed by clone formation assay. Representative images (left panel), and statistical results (right panel) are shown as means ± SD. *P < 0.05 by two-tailed Student’s t test. Data are representative of at least three independent experiments. Figure S5. circIPO11 overexpression recruits TOP1 on GLI1 promoter. (A) Heat map of downstream transcription factors in circIPO11 depleted and Ctrl HCC cells. Top 13 downregulated transcription factors were shown. (B) Different segments of GLI1 promoter were constructed into pGL3 vector for dual luciferase reporter assays. (C) Overexpression of circIPO11 enhanced luciferase activity. Luciferase activities were normalized to control luciferase activity. Results are shown as means ± SD. (D) ChIP-qPCR analysis of TOP1 enrichment on the − 2050 to − 1850 bp region of GLI1 promoter in circIPO11-overexpressing and control cells. Indicated HCC oncosphere cells were used for ChIP assay. Results are shown as means ± SD. (E) DNase I digestion assays were performed in circIPO11-overexpressing and control cells. Results are shown as means ± SD. (F) TOP1 was completely deleted in TOP1 KO cells through CRISPR/Cas9 approach. β-actin was used as a loading control. (G) The binding region of TOP1 to GLI1 promoter was confirmed by luciferase reporter assay. Results are shown as means ± SD. *P < 0.05; **P < 0.01; by two-tailed Student’s t test. Data are representative of at least three independent experiments. Representation for GLI1 knockout mice (upper panel). GLI1 deficiency was confirmed by DNA sequencing. (C) GLI1 KO efficiency in HCC cell lines were detected by Western blot. (D) GLI1-deleted or control HCC cells were transfered to 6-well plates for clone formation assay. (E) 5 × 105 CSC and non-CSC cells with GLI1-deleted or control were subcutaneously injected into BALB/c nude mice. Results are shown as means ± SD. n = 4 for each group. (F) Quantitative measurement of GLI1 activation in circIPO11-overexpression. n = 3 for each group. Results are shown as means ± SD. (G) Schematic diagram of deleting the binding region of TOP1 and GLI1 promoter by CRISPR/Cas9. Two sgRNAs were designed to target the flanking sequence of binding region. (H) CircIPO11 overexpression had no effect on sphere formation ability by GLI1 promoter deletion. Representative images and statistical results are shown. Scale bar, 500 μm. *P < 0.05; **P < 0.01, *** P < 0.001 by two-tailed Student’s t test. Data are representative of at least three independent experiments. Figure S6. GLI1 knockout impairs the stemness of liver CSCs. (A, B) Schematic representation for GLI1 knockout mice (upper panel). GLI1 deficiency was confirmed by DNA sequencing. (C) GLI1 KO efficiency in HCC cell lines were detected by Western blot. (D) GLI1-deleted or control HCC cells were transfered to 6-well plates for clone formation assay. (E) 5x105 CSC and non-CSC cells with GLI1-deleted or control were subcutaneously injected into BALB/c nude mice. Results are shown as means ± SD. n = 4 for each group. (F) Quantitative measurement of GLI1 activation in circIPO11-overexpression. n = 3 for each group. Results are shown as means ± SD. (G) Schematic diagram of deleting the binding region of TOP1 and GLI1 promoter by CRISPR/Cas9. Two sgRNAs were designed to target the flanking sequence of binding region. (H) CircIPO11 overexpression had no effect on sphere formation ability by GLI1 promoter deletion. Representative images and statistical results are shown. Scale bar, 500 μm. *P < 0.05; **P < 0.01, *** P < 0.001 by two-tailed Student’s t test. Data are representative of at least three independent experiments.Figure S7. CircIPO11 depletion and GLI1 knockout inhibit CSC proliferation. (A-C) CircIPO11, TOP1 and GLI1 mRNA levels were detected by qRT-PCR, followed by correlation analysis. Black dots represent HCC samples (n = 29). (D, E) CSC (CD13+CD133+) subpopulations non-CSC (CD13−CD133−) were sorted from HCC cell lines and transfected with respective vectors, followed by flow cytometry analysis. Ki-67 was stained for CSC proliferation (D), and statistical results are shown as means ± SD (E). n = 3 for each group. (F, G) CSC subpopulations and non-CSC were sorted from Huh7 cells and transfected with respective vectors, followed by flow cytometry analysis. 7-AAD and Annexin V was used for detection of cell death (F), and statistical results are shown as means ± SD (G). n = 3 for each group. Table S1. shRNA sequences used in this study. Table S2. sgRNA sequences used in this study. Table S3. PCR primers used in this study.
About this article
Cite this article
Gu, Y., Wang, Y., He, L. et al. Circular RNA circIPO11 drives self-renewal of liver cancer initiating cells via Hedgehog signaling. Mol Cancer 20, 132 (2021). https://doi.org/10.1186/s12943-021-01435-2