Skip to main content

Engineered circular guide RNAs boost CRISPR/Cas12a- and CRISPR/Cas13d-based DNA and RNA editing

Abstract

Background

The CRISPR/Cas12a and CRISPR/Cas13d systems are widely used for fundamental research and hold great potential for future clinical applications. However, the short half-life of guide RNAs (gRNAs), particularly free gRNAs without Cas nuclease binding, limits their editing efficiency and durability.

Results

Here, we engineer circular free gRNAs (cgRNAs) to increase their stability, and thus availability for Cas12a and Cas13d processing and loading, to boost editing. cgRNAs increases the efficiency of Cas12a-based transcription activators and genomic DNA cleavage by approximately 2.1- to 40.2-fold for single gene editing and 1.7- to 2.1-fold for multiplexed gene editing than their linear counterparts, without compromising specificity, across multiple sites and cell lines. Similarly, the RNA interference efficiency of Cas13d is increased by around 1.8-fold. In in vivo mouse liver, cgRNAs are more potent in activating gene expression and cleaving genomic DNA.

Conclusions

CgRNAs enable more efficient programmable DNA and RNA editing for Cas12a and Cas13d with broad applicability for fundamental research and gene therapy.

Background

The CRISPR-Cas nucleases are widely used for DNA and RNA editing in human and other cells and organisms and have broad applications in fundamental biological research and translational medicine [1, 2]. Three types of Cas nucleases are most commonly used: Cas9, Cas12a (formerly Cpf1), and Cas13, exemplified by Streptococcus pyogenes Cas9 (SpCas9) [3, 4], Acidaminococcus sp. Cas12a (AsCas12a) and Lachnospiraceae bacterium ND2006 Cas12a (LbCas12a) [5, 6], and Ruminococcus flavefaciens XPD3002 Cas13d (RfxCas13d, also known as CasRx) [7, 8], respectively.

Both Cas9 and Cas12a edit DNA. Different from Cas9, Cas12a exhibits several unique features. First, Cas12a exhibits higher specificity than Cas9, which enables more precise gene editing for therapeutic applications [9,10,11]. Second, Cas12a possesses crRNA self-processing capability [5, 12], enabling multiplexed gene editing with a single gRNA transcript [13,14,15], which is also a feature of Cas13 [7]. Third, rather than a G-rich protospacer adjacent motif (PAM) for Cas9, Cas12a recognizes a T-rich PAM, which makes it a good complement to Cas9 and thus broadens the genomic targeting scope [5]. Forth, Cas12a cuts target DNA with a single RuvC domain to generate sticky ends instead of with both RuvC and HNH domains for Cas9 to generate blunt ends [5, 16]. Finally, Cas12a is able to trans-cleave single-strand DNA, making it a powerful tool for nucleic acid detection [17, 18].

From the point of specificity, Cas12a-based DNA editing induces less off-target effect with safety risk than Cas9 when applied for genomic DNA sequence change, gene activation, or gene repression. As well, Cas13-based RNA editing is generally reversible and tunable without causing permanent genomic DNA changes. Therefore, Cas12a and Cas13 may hold certain advantages over Cas9 in therapeutic applications. However, relatively low editing efficiency hinders fulfilling their promise of therapeutic editing.

Both Cas nuclease and gRNA may affect editing efficiency. For example, Cas12a engineering is an effective strategy to enhance activity, like the engineered variants, Cas12a-Plus, enAsCas12a, AsCas12a-HF, AsCas12a-Ultra, and hyperCas12a [19,20,21]. However, it is believed that gRNA may play a more important role to limit editing efficiency because RNA is well known for its high sensitivity to be rapidly degraded by endonucleases and exonucleases. As proof, it has been reported that gRNAs especially free gRNAs (without Cas protein binding) are extremely unstable [22]. Chemical modification at the gRNA ends to reduce degradation by exonucleases is able to significantly enhance gene editing both in vitro and in vivo [23,24,25]. Therefore, it is essential to increase the stability and thus the abundance of gRNAs for efficient processing and loading by Cas12a and Cas13 to enhance editing efficiency.

Circular RNA is a highly stable RNA species because its covalently closed ring structure is resistant to degradation by exonucleases [26,27,28,29,30]. And circular ADAR-recruiting RNAs have been recently reported to increase the efficiency of A-to-I RNA editing [31, 32]. Therefore, we utilized circular free gRNAs and the gRNA self-processing ability of Cas12a and CasRx to increase the efficiency of Cas12a-based transcription activators as well as the DNA and RNA cleavage efficiency of Cas12a and CasRx.

Results

Stabilization of gRNA by circularization in human cells

To circularize gRNAs in mammalian cells, we adopted the elegant Tornado expression system [29]. As shown in Fig. 1a, we flanked the gRNAs by Twister ribozymes, which undergo autocatalytic cleavage, leaving termini to be ligated by the ubiquitous endogenous RNA ligase RtcB to yield circular gRNAs (cgRNAs). With gRNA self-processing ability, Cas12a and CasRx then cleaved and loaded the cgRNAs. Via mFold prediction [33], the Sp1 linker was designed to maintain the gRNA structure to be correctly recognized and processed by Cas nucleases. To visualize the expression of gRNAs, the Broccoli RNA aptamer was integrated into the circular RNA via an F30 3-way junction, which could bind the fluorophore DFHBI-1T and activate green fluorescence [29, 34]. Live staining of transfected HEK293T cells showed that circular RNAs were more abundant than their linear counterparts (Fig. 1b). Reverse transcription PCR using outward-facing primers which selectively amplified only the circular gRNAs demonstrated that gRNAs were circularized in cells (Fig. 1c). In addition, using actinomycin D treatment to block RNA transcription, we checked the stability of the circular gRNA and the linear counterparts and found that the former is much more stable than the latter (Fig. 1d). Therefore, all the above data demonstrated that circularization increased the stability of gRNAs in human cells.

Fig. 1
figure 1

Circularization increases the stability of gRNAs in human cells. a Schematic of circular guide RNAs (cgRNAs). b Broccoli fluorescence revealed the abundance of circular RNAs in cells. HEK293T cells were transfected with indicated plasmids encoding linear or circular RNAs and live stained with DFHBI-1T 48 hours (hrs) after transfection. c Reverse transcription PCR (RT-PCR) revealed circularization of RNAs in cells. HEK293T cells were transfected with indicated plasmids encoding linear or circular RNAs, and RNA was harvested 72 hrs after transfection, followed by RT-PCR with indicated outward-facing primers. d Stability of circular gRNA in cells. HEK293T cells were treated with actinomycin D for 1, 3, 6, 9, 18 hrs starting at 24 hrs post-transfection with plasmids encoding linear or circular RNAs, and RNA was harvested for quantitative RT-PCR analysis. n = 3 independent experiments

Circular gRNAs promote Cas12a-based transcriptional activation

To test whether cgRNA could enhance Cas nuclease-based editing, we first analyzed the transcription efficiency of dCas12a-based gene activators. In a doxycycline-inducible dLbCas12a-p300 knock-in (KI) HEK293T cell line (Additional file 1: Fig. S1), transiently transfected plasmids encoding cgRNAs (C-Sp1-F30 and C-Sp1) activated IL1RN and HBG gene expression more potently than the linear counterparts (L-Sp1-F30 and L-Sp1), as well as commonly used mature gRNAs (U6 + 27) and unprocessed gRNAs (Pre, spacer flanked by two scaffold sequences) (Fig. 2a). All the gRNAs were driven by the polymerase III promoter U6 + 27 cassette, which has been reported to improve the stability of small interfering RNA [35]. Time-course and dose-dependent analyses showed that cgRNAs exhibited better durability and performance than U6 + 27 and Pre gRNAs (Fig. 2b, c). Because linker sequences were essential to maintain the structure and function of cgRNAs, we optimized the length and component of linker sequences. Through screening several digital libraries generated by random sequences with RNAfold and mFold prediction [33, 36], 8 linkers were selected out for the wet-experiment test, and C-linker7 (C-L7) was found to be the best one, which was improved 5.5–12.2-fold changes than U6+27 (Fig. 2d). As expected, about 3.3–5.7-fold changes increased efficiency was also observed in MCF7 cells (Additional file 1: Fig. S2a). We further demonstrated that cgRNAs were applicable for other dLbCas12a-based activators, including dLbCas12a-VPR, -SunTag-VP64, and -SunTag-VPR (Fig. 2e–g; Additional file 1: Fig. S2b). Finally, RNA-seq analyses showed that only the IL1RN target gene was significantly activated, indicating a high specificity for cgRNA-directed gene activation (Fig. 2h; Additional file 1: Fig. S2c). In summary, cgRNAs could enhance the activity and maintain the high specificity of Cas12a-based gene activators.

Fig. 2
figure 2

Circular guide RNAs increase the transcription efficiency of Cas12a-based activators. a cgRNA-directed gene activation in a dLbCas12a-p300 knock-in (KI) HEK293T cell line. b Time-course analyses of cgRNA-directed gene activation in the KI cells. c Dose-dependent analyses of cgRNA-directed gene activation in the dLbCas12a-p300 knock-in HEK293T cells. d Gene activation guided by cgRNAs with different linkers in the KI cells. e–g cgRNA-directed gene activation with a variant of dLbCas12a-based gene activators in HEK293T cells transiently co-transfected with indicated activator-encoded and gRNA-encoded plasmids. h The specificity of cgRNA-directed gene activation. Gene expression plot generated from RNA-seq data from the KI HEK293T cells transfected with U6 + 27 linear gRNAs or C-L7 cgRNAs targeting mNeonGreen (control) or IL1RN. R indicates Pearson’s correlation coefficient. The average of three biological replicates was shown. For a–g, quantitative RT-PCR revealed relative mRNA expression of IL1RN and HBG. Mean values are presented with S.D., n = 3 independent experiments. For each experiment, fold changes of mRNA expression in tested samples versus that in the U6 + 27 linear mNeonGreen gRNA were shown. *p <0.05, **p <0.01, ***p <0.001, one-way ANOVA test

Circular gRNAs increase the DNA cleavage efficiency of Cas12a

Next, we tested whether cgRNA could enhance Cas12a-mediated DNA cleavage. We designed different forms of gRNAs targeting the same site within the mNeonGreen gene and co-transfected these plasmids with LbCas12a-mCherry into a mNeonGreen reporter HEK293T cell line. Fluorescence-activated cell sorting (FACS) analyses showed that C-L7 for LbCas12a significantly enhanced DNA editing activity, with about 1.4- and 1.5-fold than U6+27 and Pre gRNA (Fig. 3a, b; Additional file 1: Fig. S3a). Next, using Tag-seq [37] and a doxycycline-induced LbCas12a knock-in HEK293T cell line (Additional file 1: Fig. S3b-d), transiently transfected plasmids encoding cgRNAs showed higher efficiency and slight lower specificity compared with linear U6 + 27 gRNAs across 14 genomic sites (Fig. 3c–f; Additional file 1: Fig. S3e), with about 1.75- and 2.10-fold for C-Sp1 and C-L7 in efficiency, respectively (Fig. 3c, d), and with 3, 5, and 4 off-target sites as well as specificity index (the ratio of total on-target reads to the on-target reads plus the off-target reads) of 0.98, 0.95, and 0.94 for U6 + 27, C-Sp1, and C-L7, respectively (Fig. 3e, f; Additional file 1: Fig. S3e). Similar results were observed in MCF7 cells (Additional file 1: Fig. S4).

Fig. 3
figure 3

Circular guide RNAs increase the DNA cleavage efficiency of Cas12a. a Fluorescence-activated cell sorting (FACS) analyses of the mNeonGreen reporter cells 4 days after co-transfection with LbCas12a-P2A-mCherry and mNeonGreen-targeting-gRNA plasmids. b The cleavage efficiency was quantified by the cell ratio of mNeonGreen mCherry+ / mCherry+ in the FACS assays. n = 3 independent experiments. NT gRNA, non-targeting gRNA, which recognized no site in the human genome and transcriptome. ***p <0.001, one-way ANOVA test. c–f The efficiency and specificity of different gRNAs-directed DNA cleavage at 14 sites in a LbCas12a knock-in HEK293T cell line revealed by Tag-seq. The gRNA reference as well as the on-target and off-target sites was shown on the left, and sequencing read counts were shown to the right of each site (c). Efficiency comparison between different gRNAs (d). The total number of off-target sites detected for the 14 sites (e). Specificity index (value was calculated by the ratio of total on-target reads to the on-target reads plus the off-target reads within the 14 sites) (f)

Next, we performed comparison between cgRNA of other extended structuralized gRNAs. According to reports, extending the 5′ end of the crRNA by 9 nucleotides (gRNA + 9) and 59 nucleotides (gRNA + 59) could enhance the gene editing efficiency of the Cas12a ribonucleoproteins complexes (RNP) [38]. And the crRNA with a U4AU4 3′-overhang was more favorable binding to Cas12a to improve the activity [39]. In the doxycycline-inducible dLbCas12a-p300 knock-in HEK293T cell line, transiently transfected plasmids encoding gRNAs, cgRNAs (C-L7-d27, without 27 nucleotides at the 5′ end, and C-L7) showed more potently activation than the 5′ end extended gRNA (gRNA + 9 and gRNA + 59) with about 8.8–11.4-fold change in gene IL1RN and 4.9–6.3-fold change in gene HBG, as well as 3′ end extended gRNA (gRNA + T4AT4) with about 4.0–7.6-fold change in gene IL1RN and 6.5–9.1-fold change in gene HBG (Additional file 1: Fig. S5a). Moreover, we tested DNA cleavage and found that cgRNA (C-L7-d27 and C-L7) performed 1.6–1.9-fold change activity than 5′ end extended gRNA (+ 9 and + 59) and 1.5–1.7-fold change activity than 3′ end extended gRNA (gRNA + T4AT4) to disrupt mNeonGreen in the HEK293T reporter cell line (Additional file 1: Fig. S5b, c). Likewise, cgRNA (C-L7-d27 and C-L7) performed about 1.6–2.0-fold change activity than 5′ end extended gRNA (+ 9 and + 59) and about 1.3–1.5-fold change activity than 3′ end extended gRNA (gRNA + T4AT4) to disrupt endogenous gene VEGFA in the doxycycline-inducible LbCas12a knock-in HEK293T cell line (Additional file 1: Fig. S5d, e).

AsCas12a is another popular Cas12a nuclease and thus we tested whether cgRNAs were applicable to it. The transcriptional activator dAsCas12a-VPR activated IL1RN and HBG more efficiently when coupled with cgRNAs, including C-Sp1-F30, C-Sp1, and 4 cgRNAs with different linkers which were screened out of digital libraries via RNAfold and mFold prediction (Additional file 1: Fig. S6a, b). Since the cgRNA C-L1 performed best among all the tested 6 cgRNAs for gene activation, we tested its performance in DNA cleavage and found that C-L1 outperformed the counterpart linear gRNAs to disrupt mNeonGreen in the HEK293T reporter cell line (Additional file 1: Fig. S6c-e).

Circular gRNAs improve Cas12a-mediated multiplexed gene activation and DNA cleavage

Multiplexed editing with a single gRNA transcript is a unique feature of Cas12a over Cas9; therefore, we tested whether cgRNAs were compatible or even improved multiplexed editing. The gRNAs targeting the promoter regions of NTT, IL1RN, and HBG were co-expressed within a single linear (Pre-NIH) or circular (C-L7-NIH) transcript and C-L7-NIH significantly improved dLbCas12a-VPR-mediated transcriptional activation of each of the three genes (Additional file 1: Fig. S7a). Similarly, Deep-seq revealed that the circular transcript (C-L7-CVDER) containing 5 gRNAs targeting CD47-S3, VEGFA-S1, DNMT1-S3, EMX1, and RUNX1 was more potent to direct LbCas12a-mediated DNA cleavage than the linear transcript (Pre-CVDER, Additional file 1: Fig. S7b). Tag-seq also showed that C-L7-CVDER exhibited a 1.67-fold efficiency and slight lower specificity when compared to Pre-CVDER (Additional file 1: Fig. S7c-g).

In summary, all the above results demonstrated that circular gRNAs were able to increase the efficiency and maintain almost equal specificity of LbCas12a- and AsCas12a-based gene activation and DNA cleavage.

Circular gRNAs increase the RNA cleavage efficiency of CasRx

Further, we tested whether cgRNAs could enhance CasRx-mediated RNA cleavage. To this end, we selected 4 circular RNA backbones with different linkers out of digital libraries via RNAfold and mFold prediction and did similar experiments in the mNeonGreen reporter cell line. FACS analyses and reverse transcription PCR (RT-PCR) showed that cgRNAs significantly increased RNA cleavage efficiency with about 1.2–1.4-fold change (Fig. 4a–c; Additional file 1: Fig. S8a). For endogenous genes (STAT3, NF2, B4GALNT1, KRAS, and RPL4) in HEK293T cells, C-L1 cgRNAs were also observed to repress their expression more potently than the linear counterpart Pre gRNAs with about 1.2–2.8-fold change (Fig. 4d). Similar results were observed in MCF7 cells (Additional file 1: Fig. S8b). Moreover, RNA-seq analyses showed that NF2 (the target gene) was significantly decreased in the C-L1 cgRNA sample compared to the Pre and NT gRNA samples with 2.4-fold and 3.6-fold change, respectively (Fig. 4e; Additional file 1: Fig. S8c). Of note, 64 and 126 differentially expressed genes (DEGs) were observed when comparing the C-L1 cgRNA group to the NT gRNA group and to the Pre gRNA group, respectively, indicating nonspecific and collateral cutting as well as the associated secondary effects [40, 41].

Fig. 4
figure 4

Circular guide RNAs increase the RNA cleavage efficiency of CasRx. a, b The RNA cleavage efficiency of RfxCas13d (CasRx) on the mNeonGreen reporter gene. FACS analyses of the mNeonGreen knock-in HEK293T cell line 48 hrs after co-transfection with CasRx-P2A-mCherry and mNeonGreen-targeting-gRNA plasmids (a). The cleavage efficiency was quantified by the cell ratio of mNeonGreen mCherry+ / mCherry+ and mean fluorescence intensity (MFI) of mNeonGreen in transfected cells (mCherry positive) in the FACS assays (b). c Relative degradation of mNeonGreen transcripts induced by CasRx with circular or linear gRNAs. The mRNA expression levels were determined by RT-PCR. d The RNA cleavage efficiency of CasRx on endogenous genes. HEK293T cells were co-transfected with CasRx-P2A-mCherry and gRNA plasmids, and mCherry+ cells were sorted out by FACS for RNA extraction and quantitative RT-PCR analyses. n = 3 independent experiments. NT, Non-targeting gRNA, which recognized no site in the human genome and transcriptome. e The specificity of cgRNA-directed gene degradation. Gene expression plot generated from RNA-seq data from HEK293T cells transfected with non-targeting cgRNA or Pre linear gRNAs or C-L1 cgRNAs targeting NF2. R indicates Pearson’s correlation coefficient. The average of three biological replicates was shown. For b,c, **p <0.01, ***p <0.001, one-way ANOVA test. For d, **p <0.01, ***p <0.001, Student’s t test

Collateral effect is a feature of Cas13 protein which may hindered their application in vivo [41]. To evaluate the effect, we first tested the trans-cleavage activity (cleavage of mCherry) of CasRx when targeting the exogenously expressed mNeonGreen gene. As expected, cgRNAs showed significantly decreased mNeonGreen fluorescence intensity, and with a similar trans-cleavage activity compared with linear U6 +27 gRNA (Additional file 1: Fig. S9). Next, we tested the trans-cleavage activity of CasRx when targeting endogenous genes. Consistent with the previous report [40], RPL4 gRNA induced dramatic collateral effects. However, only slight collateral degradation was observed when targeting STAT3 and NF2, and no apparent collateral degradation was observed for B4GALNT1 and KRAS. More importantly, we observed no significant difference between cgRNAs and linear Pre gRNA (Additional file 1: Fig. S10).

Circular gRNAs enhance Cas12a-based transcriptional activation in vivo

Next, we explored whether cgRNAs could improve gene activation in vivo. To this end, we constructed the TRE-Luciferase-pA plasmid, which contained 6 gRNA binding sites within the TRE promoter region. The TRE-Luciferase-pA plasmid and dLbCas12a-VPR, as well as a control cgRNA targeting mNeonGreen or a linear gRNA or a C-L7 cgRNA targeting the TRE promoter, were co-delivered to mouse liver via hydrodynamic tail vein injection (HTVI) [42], and from the next day, luciferase activity was live examined each day for a consecutive 8 days (Fig. 5a). Similar to the in vitro results, cgRNA outperformed linear gRNA at each time point with about 9.9–32.6-fold change (Fig. 5b, c).

Fig. 5
figure 5

Circular gRNAs enhance the activation efficiency of dLbCas12a-VPR in vivo. a Experiment design for dLbCas12a-VPR to activate Luciferase expression in mouse liver. gCtrl, control gRNA, which targeted mNeonGreen. HTVI, hydrodynamic tail vein injection. b Representative bioluminescence imaging results at day 4 for all the 8 groups of mice. c Quantification of bioluminescence imaging detected for a consecutive 8 days. *p <0.05, **p <0.01, ***p <0.001, Student’s t test

Circular gRNAs enhance the DNA cleavage efficiency of Cas12a in vivo

Finally, we explored whether cgRNAs could improve DNA cleavage in vivo. To this end, we adopted the Cas-N57 system developed by our group to induce tumorigenesis in mouse liver [42]. Because KRAS, TP53, and PTEN mutants are the major drivers of intrahepatic cholangiocarcinoma (ICC) [43], we used LbCas12a-N57 to simultaneously insert KrasG12D into the Rosa26 site and disrupt Trp53 and Pten in the mouse liver. As shown in Fig. 6a, the KrasG12D donor and LbCas12a-N57 plasmids, as well as the plasmid encoding a C-L7 circular multiplexed transcript or a Pre linear multiplexed transcript containing 2 gRNAs for Pten, 2 gRNAs for Trp53, and 1 gRNA for Rosa26 or encoding a control cgRNA targeting mNeonGreen, were co-delivered to mouse liver via hydrodynamic tail vein injection and tumorigenesis was examined 7 weeks after injection. The C-L7 cgRNA group showed more tumor nodules and the liver weighed more (Fig. 6b, c). H&E and immunochemical staining of the tumor nodules showed pathological features of bile duct differentiation and expression of the ICC marker cytokeratin 19 (Fig. 6d). Sanger sequencing and Tracking of indels by decomposition (TIDE) analyses of genomic PCR amplicons from tumor nodules showed that the C-L7 cgRNA induced more indels than the linear Pre gRNA (Fig. 6e). In addition, Sanger sequencing also demonstrated targeted insertion of KrasG12D into the Rosa26 locus (Fig. 6f, g). Collectively, these results suggested that circular gRNA was more potent than linear gRNA to cleave genomic DNA in vivo.

Fig. 6
figure 6

Circular gRNAs enhance the DNA cleavage efficiency of LbCas12a in vivo. a Experiment design for LbCas12a-N57 to induce liver tumor in adult mice. gCtrl, control gRNA, which targeted mNeonGreen. HTVI, hydrodynamic tail vein injection. b Image analysis of mouse liver harvested 7 weeks after injection. c Quantification of liver weight and surface liver tumor nodules per mouse. n = 7. d Representative images of H&E and IHC staining of Ck19 in mouse liver tumors. Scale bar, 100 μm. e Sanger sequencing results and tracking of indels by decomposition (TIDE) analyses of tumor DNA for Pten and Trp53 targeted sites. Blue shadow denoted the gRNA recognizing sites. f Genomic PCR of targeted integration of KRASG12D donor in tumors. Primers for the 5′-junction and 3′-junction were R26–5-F/R26–5-R and R26–3-F/R26–3-R, respectively. g Verification of the presence of the KrasG12D mutation in the tumors by Sanger sequencing. *p <0.05, **p <0.01, Student’s t test

Discussion

In this study, using the Tornado expression system [29], we generated cgRNAs, which were much more stable than linear gRNAs. cgRNAs increased the transcription efficiency of dLbCas12a-based gene activators by up to ~40-fold, including dLbCas12a-p300, dLbCas12a-VPR, dLbCas12a-SunTag-VP64, and dLbCas12a-SunTag-VPR systems. Further, cgRNA increased the genomic DNA cleavage efficiency of LbCas12a by ~2-fold. Apart from single gene editing, multiplexed gene editing with a single gRNA transcript was also enhanced by cgRNAs. More importantly, enhanced activity did not compromise specificity. And the enhancement phenomenon was also observed with AsCas12a, another popular Cas12a nuclease. Similarly, the RNA interference efficiency of CasRx was increased by ~2-fold when directed by cgRNAs and similar collateral effect was maintained. And finally, in mouse liver, cgRNAs were more potent to activate gene expression and were able to enhance Cas12a-mediated KrasG12D insertion and Pten and Trp53 disruption to promote tumorigenesis.

Owing to the delivery challenge associated with in vivo therapeutic gene editing, target cells generally uptake much less nucleic acid or ribonucleoprotein (RNP) encoding the editing tools than the cultured cells for in vitro or ex vivo gene editing [44]. Therefore, developing editing tools with high efficiency is an essential approach to fulfill the therapeutic efficacy of in vivo gene editing, apart from improving the delivery method itself. In this study, we engineered cgRNAs to enhance the efficiency of Cas12a-based gene activation and genomic DNA cleavage. The cgRNA system performed well for cultured cells, achieving durable editing for at least 8 days and efficient editing with extremely low amount of cgRNAs (Fig. 2b, c). Encouraged by these results in vitro, we performed gene activation and gene disruption experiments in adult mouse liver and found that cgRNAs significantly outperformed the linear gRNA counterparts (Figs. 5 and 6). These observations indicated a great potential for cgRNAs to be used for therapeutic gene editing in vivo in the future.

Because the secondary structure of RNA is essential for both ribozyme self-cleavage for circularization and Cas12a and CasRx self-processing and loading, the linker 1 and linker 2 between the ribozymes and the gRNA need to be optimized to maintain the correct RNA structure (Fig. 1a). We adopted adenine- and cytosine-rich (AC-rich) RNA linkers as they have been widely used as flexible RNA linkers [30, 32]. To further increase the performance, we generated several digital cgRNA libraries with variant linkers and screened out cgRNAs with the predicted most stable correct structure via RNAfold and mFold prediction and verified by wet experiments (Figs. 2d and 4a, b, Additional file 1: Fig. S6b, Fig. S11-13, and Linker design in the “ Methods” section). Restricted by the server computing power, we only screened limited digital cgRNA libraries. It is likely to find out better cgRNA linkers to further enhance editing efficiency after more digital libraries were screened. Moreover, we would suggest confirming the cgRNA structure by RNAfold and mFold prediction when designing a particular gRNA since the spacer sequence might interfere with the RNA secondary structure.

With gRNA self-processing capability, Cas12a and Cas13 are theoretically unable to form protein/gRNA complex with circular gRNA, and circularization of free gRNA is thus a proper approach to boost their editing activity. Without gRNA self-processing capability, it would be hard for Cas9 to load circular gRNA or to execute efficient editing when bound to circular gRNA. Our preliminary data indicated that Cas9 and circular gRNA might form a complex to edit genes with less efficiency. Extensive optimization of the gRNA scaffold and circular RNA linker facilitated by AI-assisted RNA and protein structure prediction might solve this problem in the near future [45, 46].

Conclusion

In summary, we engineered free circular gRNA to boost programmable DNA and RNA editing for Cas12a and CasRx, which might have broad applications for fundamental biological research and translational medicine.

Methods

Plasmid construction

For U6+27 gRNA-expressing constructs, the DNA sequences including the hU6+27 promoter, a gRNA scaffold, and a gRNA insert site were synthesized and cloned into the pBluSKM vector. For Pre gRNA-expressing constructs, a second gRNA scaffold was inserted downstream of the gRNA insert site in the U6+27 gRNA-expressing constructs. For linear Sp1-F30 gRNA-expressing constructs, the DNA sequences including a 5' ligation sequence, a 5' 42-nt spacer, an F30 3-way junction, a broccoli sequence, two gRNA scaffolds with a gRNA insert site, a 3' 41-nt spacer, and a 3' ligation sequence were synthesized and cloned into the pBluSKM vector with the hU6+27 promoter. For circular sp1-F30 gRNA-expressing constructs, a Twister P3 U2A ribozyme sequence and a Twister P1 ribozyme sequence were inserted to flank the linear sp1-F30 gRNA-expressing constructs. For linear and circular sp1 gRNA-expressing constructs, the F30 3-way junction and the broccoli sequences were deleted in the corresponding F30 plasmids. For U6 gRNA-expressing constructs, 27-nt in the 3′ end of U6+27 was deleted in the corresponding U6+27 plasmids. For C-L7-d27 gRNA-expressing constructs, 27-nt in the 3′ end of U6+27 was deleted in the corresponding C-L7-d27 plasmids. For +9 or +59 gRNA-expressing constructs, a 9-nt or 59-nt 5′ extended sequence was added by Gibson Assembly in the corresponding U6 plasmids. For + T4AT4 gRNA-expressing constructs, the 3′ extended sequence (TTTTATTTT) was added by Gibson Assembly in the corresponding U6 plasmids.

For dLbCas12a-p300 and LbCas12a knock-in plasmids, the sequences of dLbCas12a-p300 and LbCas12a were amplified from the pCAG-dLbCas12a-p300-mCherry and the pCAG-LbCas12a-mCherry plasmids, respectively, then cloned into the pBlue-AAVS1-Puro-Cas9:p300-M2rtTA-AAVS1 plasmid [47] to replace Cas9:p300, next, two AAVS1 gRNA targeting sequences were inserted to flank the whole knock-in fragment.

For the plasmid pEF1α-CasRx-mCherry, the CasRx sequence was amplified from plasmid EF1a-CasRx-2A-EGFP (addgene, #109049), and used to replace LbCas12a sequence in plasmid pCAG-LbCas12a-mCherry.

For the luciferase reporter, the Tet operator containing six gRNA binding sites was amplified from the plasmid pBlue-AAVS1-Puro-Cas9:p300-M2rtTA-AAVS1 and the sequence of luciferase-polyA was synthesized, and then the two sequences were sub-cloned into the pBluSKM vector.

For the KrasG12D donor, the LbCas12a gRNA targeting sequence within the Rosa26 locus was synthesized and replaced the original sequence in the plasmid pBlue-Rosa26-IR-T2A-Puro- pEF1α-KrasG12D-IR-Rosa26 [42].

The plasmid pCAG-dLbCas12a-p300-mCherry, pCAG-dLbCas12a-VPR-mCherry, pCAG-dLbCas12a-10×Suntag-mCherry, pHRdSV40-scFv-GCN4-sfGFP-VP94-GB1-NLS, pHRdSV40-scFv-GCN4-sfGFP-VPR-GB1-NLS, pCAG-LbCas12a-mCherry, pCAG-LbCas12a-N57-mCherry, pCAG-dAsCas12a-VPR-mCherry, pCAG-AsCas12a-mCherry and pBlue-Rosa26-IR-T2A-Puro- pEF1α-KrasG12D-IR-Rosa26 were described previously [37, 42, 48, 49].

All sgRNAs were designed through https://benchling.com/ and ligated to the corresponding sgRNA expression plasmid. The sequences of all sgRNAs, linkers, and extended structures are listed in Additional file 2: Table S1. All constructs were verified through Sanger sequencing.

Linker design

To increase editing efficiency, different linkers were designed by screening several digital libraries. Briefly, linker libraries were constructed with 5' and 3' linkers containing a 10-nt polyAC sequence and 5–7nt random bases, and the structure with each linker pair was predicted by RNAfold. The ones predicted to contain correct ribozyme structure and gRNA scaffold structure were selected out and ranked by the Gibbs free energy change (ΔG). The top 10–20 circular RNAs in the rank list were divided into groups according to structural similarity. The structure of group members was predicted by mFold, another RNA structure prediction tool. In each group, the circular RNA with the lowest ΔG (most stable) and correct structure predicted by both RNAfold and mFold was selected for the wet-experiment test.

Cell culture and transfection

Cell lines were obtained from ATCC and regularly checked for mycoplasma. HEK293T and B16 cells were cultured in DMEM medium (Life Technologies), and MCF7 cells were maintained in RPMI 1640 medium (Life Technologies) at 37 °C under 5% CO2. All growth media were supplemented with 100 U/mL penicillin, 100 μg/mL streptomycin (Life Technologies), and 10% FBS.

For transient transfection experiments, plasmids were transfected into cells by polyethylenimine (PEI). Cells were plated in 24-well plates and transfected when reached about 70% confluence. For the mNeonGreen reporter HEK293T cells, 15.6.ng gRNA plasmid, 46.9 ng Cas12a / CasRx plasmid, and 62.5ng pBluSKM plasmid were mixed and co-transfected into cells to test the DNA or RNA cleavage of mNeonGreen. For the Tet-dLbCas12a-p300 knock-in HEK293T cells, 15.6ng gRNA plasmid and 484.4ng pBluSKM plasmid were mixed and co-transfected into cells to test target gene activation. For non-KI HEK293T or MCF7 cells, 15.6ng gRNA plasmid, 46.9ng LbCas12a based activator plasmid, and 437.5ng pBluSKM plasmid were mixed and co-transfected. For the Tet-LbCas12a knock-in HEK293T cells, 62.5ng gRNA plasmid and 437.5ng pBluSKM plasmid were mixed and co-transfected to test endogenous gene cleavage. For CasRx-mediated endogenous gene degradation, 62.5ng gRNA plasmid, 187.5ng CasRx plasmid, and 250ng pBluSKM plasmid were mixed and co-transfected into HEK293T cells or MCF7 cells. Two days (for CasRx) or 4 days (for Cas12a) after transfection, cells were harvested for gene expression analyses, flow cytometry, Deep-seq, Tag-seq, RNA-seq, or TIDE analyses etc.

The Tet-dLbCas12a-p300 and Tet-LbCas12a knock-in HEK293T cell lines were obtained by transfecting corresponding plasmids and selecting positive clones with 1μg/ ml puromycin and were added 2μg/ ml doxycycline to induce Cas protein expression.

For testing the stability of gRNAs, 5 μg/ ml actinomycin D was added to cells 24 h after transfections, then cells were harvested at indicated time points to analyze linear or circular gRNA content.

Quantitative real-time PCR

Total RNA was extracted from cells using TRIzol reagent (Invitrogen) according to the manufacturer’s protocol. Briefly, the total RNA was extracted from cells by adding 500 μl TRIzol and 100μl chloroform, after centrifugation at 13,000 rpm for 10min at 4°C, the supernatant was transferred to a 1.5-mL RNase-free tube. RNA was purified by precipitation with isopropanol and 75% ethanol. Five hundred nanogram RNA was reverse transcribed using Prime Script™ RT Reagent Kit (TAKARA). cDNA was diluted 10-fold, and 2.0µl diluted cDNA was used for each RT-PCR reaction with TB Green Premix Ex Taq II Kit (TAKARA) and Lightcycler 96 (Roche). The primers are listed in Additional file 2: Table S1.

RNA-seq

Gene activation of dLbCas12a-p300 activator and RNA cleavage of CasRx were determined by RNA-seq as described previously [48]. Briefly, total RNA was isolated from cells using TRIzol reagent (Invitrogen), purified by magnetic beads with Oligo(dT), and random fragmented by fragmentation buffer. The first strand was synthesized by six-base random hexamers, and then followed by the second strand synthesis. After purification, terminal repair, and dA-tailing addition, and adaptor addition, double-strand cDNA was amplified by PCR to complete library construction. After quality verified using Qubit 3.0, Agilent 2100 Bioanalyzer and agarose gel electrophoresis, libraries were sequenced by Illumina HiSeq instrument with a 150-bp paired-end module. Significant differentially expressed genes were defined with a false discovery rate (FDR) < 0.05 and a fold change > 2.

Tag-seq

The detailed procedures for Tag-seq library construction and analyses were described in our previous work [37]. Briefly, genomic DNA (gDNA) was isolated from cells using TIANamp Genomic DNA Kit (TIANGEN). Then the genomic DNA was fragmented, end-repaired, and ligated by dA-tail. The library was constructed with Nested PCR and sequenced with Illumina HiSeq instrument with a 150-bp paired-end module. The sequencing data were analyzed through Tag filtering, quality control, read alignment, PCR duplicate consolidation, and identification of RGN-mediated off-target cleavage sites. The Tag-seq data analysis pipeline is available at GitHub (https://github.com/zhoujj2013/Tag-seq).

Deep-seq

Multiplexed DNA cleavage efficiency of LbCas12a was determined by Deep-seq as described previously [48]. Briefly, the primers with forward and reverse indexes were used to amplify on-target sites in the first-round PCR. Then, equal amounts of the first PCR products were mixed for a second round PCR with the P5- and P7-containing primers to generate the libraries. The library was sequenced by Illumina HiSeq instrument with 150-bp paired-end reads. Pooled samples were demultiplexed by the indexes within the primers for the first-round PCR. Sequencing reads were trimmed, mapped to, and aligned with the genome reference. Indels were called using the R package Genomic Alignments [50]. The primers are listed in Additional file 2: Table S1.

FACS analysis

In the mNeonGreen reporter assay, cells were harvested 2 or 4 days after transfection and resuspended in 400 μl FACS buffer (1× DPBS, 0.2% BSA, 2mM EDTA), then loaded onto a flow cytometer (BD Fortessa, CA, USA) to detect mNeonGreen- and mCherry-positive cells. The cleavage efficiency of Cas12a or CasRx was calculated as the proportion of mNeonGreen-negative cells and the mean fluorescence intensity (MFI) of mNeonGreen within the transfected cells (mCherry positive). For the collateral effects assay of CasRx, MFI of mCherry and mNeonGreen of total cells were analyzed. For examining the cleavage efficiency of CasRx when targeting endogenous genes, transfected cells (mCherry positive) were sorted out for RNA extraction by MoFlo XDP flow cytometry sorter 48h after transfection.

DFHBI-1T staining and microscopy

Cells were observed 2 days after transfection. Thirty minutes before imaging, the culture medium was changed to FluoroBrite medium (ThermoFisher) with 40 μM DFHBI-1T and 0.1 μg/ml Hoechst. Live cell fluorescence images were acquired on a Nikon microscope.

Western blotting

Cells were lysed in 2X SDS loading buffer (200 mM β-mercaptoethanol, 100 mM Tris pH 6.8, 20% glycerol, 4% SDS, 0.05% bromophenol blue). The lysates were separated by SDS-PAGE and transferred onto the NC membrane, followed by blocking with 5% milk in TBST solution and incubation with primary antibody overnight at 4 °C and secondary antibodies for 1h at room temperature. Finally, the NC membrane was incubated with Immobilon Western Chemiluminescent HRP Substrate (Millipore) and imaged by Gel Imager System. Antibodies included Anti-HA-tag antibody (MBL, M180-3) for dLbCas12a-p300, dLbCas12a-VPR, dLbCas12a-SunTag, scFv-sfGFP-VP64, scFv-sfGFP-VPR, LbCas12a, dAsCas12a-VPR, AsCas12a, and CasRx. Anti-Tubulin antibody (Proteintech, 66240), GAPDH antibody (Proteintech, 60004), and HRP-conjugated horse anti-mouse IgG secondary antibody (CST 7076S).

H&E staining and immunohistochemistry

Samples were fixed overnight in 4% paraformaldehyde at 4°C, embedded with paraffin, and then sliced into 5-μm sections. For H&E staining, the sections were rehydrated with gradient ethanol and stained with hematoxylin and eosin. For immunohistochemistry, the rehydrated sections were boiled for 15 min to retrieve antigen. After the endogenous peroxidase was blocked for 15 min, the sections were sequentially incubated with a primary antibody (Anti-CK19, ab133496, Abcam) overnight at 4 °C, a secondary antibody (HRP-conjugated anti-rabbit IgG secondary antibody, 7074S, CST) for 1h at room temperature, and the chromogenic substrate for 20 min. Finally, the sections were counterstained with hematoxylin, dehydrated, and sealed with neutral resins.

Mice

C57BL/6 female mice were purchased from Guangdong Animal Center. In the luciferase reporter assay, each mouse aged 6 weeks was injected with a total of 26 μg plasmids (10 μg luciferase reporter plasmid, 12 μg pCAG-dLbCas12a-VPR-mCherry plasmid, and 4 μg sgRNA plasmid were mixed in 1.8-2.0 ml 0.9% sterile saline) via hydrodynamic tail vein injection (HTVI). Fifteen minutes before imaging, each group which included three mice injected with different gRNAs was treated with 200 μl 15mg/ml D-luciferin potassium salt (Beyotime) and then imaged with an exposure time of 10 s. In the cancer modeling assay, each mouse aged 6 weeks was injected with a total of 55 μg plasmids (30 μg KrasG12D donor plasmid, 20 μg pCAG-LbCas12a-N57-mCherry plasmid, and 5 μg sgRNA plasmid were mixed in 1.8–2.0 ml 0.9% sterile saline) via hydrodynamic tail vein injection. After 7 weeks, the mice were sacrificed for assessment. The livers were harvested for weighing, genomic DNA extraction, H&E staining, and immunohistochemistry.

Availability of data and materials

All sgRNAs, linkers, and primer sequences in this study are available in the Additional file 2: Table S1. RNA-seq, Tag-seq, and Deep-seq data have been deposited on the National Center for Biotechnology Information database (accession nos. PRJNA830337 [51]). All code for RNA-seq analysis pipeline in this study is available at https://github.com/YuchenLiu1621/circular_RNA and https://doi.org/10.5281/zenodo.7991279 [52]. All code for Deep-seq analysis pipeline in this study is available at https://github.com/TZH0511/deepseq and https://doi.org/10.5281/zenodo.7992742 [53]. All code for Tag-seq analysis pipeline in this study is available at https://github.com/zhoujj2013/Tag-seq and https://doi.org/10.5281/zenodo.4679460 [54]. Any updates will also be published on Zenodo and GitHub. All uncropped versions of the gel and microscopy images are available in FigShare (https://doi.org/10.6084/m9.figshare.23268527 [55]).

References

  1. Pickar-Oliver A, Gersbach CA. The next generation of CRISPR-Cas technologies and applications. Nat Rev Mol Cell Biol. 2019;20:490–507.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Anzalone AV, Koblan LW, Liu DR. Genome editing with CRISPR-Cas nucleases, base editors, transposases and prime editors. Nat Biotechnol. 2020;38:824–44.

    Article  CAS  PubMed  Google Scholar 

  3. Cong L, Ran FA, Cox D, Lin S, Barretto R, Habib N, et al. Multiplex genome engineering using CRISPR/Cas systems. Science. 2013;339:819–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Mali P, Yang L, Esvelt KM, Aach J, Guell M, DiCarlo JE, Norville JE, Church GM. RNA-guided human genome engineering via Cas9. Science. 2013;339:823–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Zetsche B, Gootenberg JS, Abudayyeh OO, Slaymaker IM, Makarova KS, Essletzbichler P, et al. Cpf1 is a single RNA-guided endonuclease of a class 2 CRISPR-Cas system. Cell. 2015;163:759–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Safari F, Zare K, Negahdaripour M, Barekati-Mowahed M, Ghasemi Y. CRISPR Cpf1 proteins: structure, function and implications for genome editing. Cell Biosci. 2019;9:36.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Konermann S, Lotfy P, Brideau NJ, Oki J, Shokhirev MN, Hsu PD. Transcriptome engineering with RNA-targeting type VI-D CRISPR effectors. Cell. 2018;173:665-676 e614.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Buchman AB, Brogan DJ, Sun R, Yang T, Hsu PD, Akbari OS. Programmable RNA targeting using CasRx in flies. CRISPR J. 2020;3:164–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Yang M, Wei H, Wang Y, Deng J, Tang Y, Zhou L, Guo G, Tong A. Targeted disruption of V600E-Mutant BRAF gene by CRISPR-Cpf1. Mol Ther Nucleic Acids. 2017;8:450–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Kleinstiver BP, Tsai SQ, Prew MS, Nguyen NT, Welch MM, Lopez JM, McCaw ZR, Aryee MJ, Joung JK. Genome-wide specificities of CRISPR-Cas Cpf1 nucleases in human cells. Nat Biotechnol. 2016;34:869–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Kim D, Kim J, Hur JK, Been KW, Yoon SH, Kim JS. Genome-wide analysis reveals specificities of Cpf1 endonucleases in human cells. Nat Biotechnol. 2016;34:863–8.

    Article  CAS  PubMed  Google Scholar 

  12. Fonfara I, Richter H, Bratovic M, Le Rhun A, Charpentier E. The CRISPR-associated DNA-cleaving enzyme Cpf1 also processes precursor CRISPR RNA. Nature. 2016;532:517–21.

    Article  CAS  PubMed  Google Scholar 

  13. Campa CC, Weisbach NR, Santinha AJ, Incarnato D, Platt RJ. Multiplexed genome engineering by Cas12a and CRISPR arrays encoded on single transcripts. Nat Methods. 2019;16:887–93.

    Article  CAS  PubMed  Google Scholar 

  14. Zetsche B, Heidenreich M, Mohanraju P, Fedorova I, Kneppers J, DeGennaro EM, et al. Multiplex gene editing by CRISPR-Cpf1 using a single crRNA array. Nat Biotechnol. 2017;35:31–4.

    Article  CAS  PubMed  Google Scholar 

  15. Huang H, Huang G, Tan Z, Hu Y, Shan L, Zhou J, et al. Engineered Cas12a-Plus nuclease enables gene editing with enhanced activity and specificity. BMC Biol. 2022;20:91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Nishimasu H, Ran FA, Hsu PD, Konermann S, Shehata SI, Dohmae N, Ishitani R, Zhang F, Nureki O. Crystal structure of Cas9 in complex with guide RNA and target DNA. Cell. 2014;156:935–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Chen JS, Ma E, Harrington LB, Da Costa M, Tian X, Palefsky JM, Doudna JA. CRISPR-Cas12a target binding unleashes indiscriminate single-stranded DNase activity. Science. 2018;360:436–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Ooi KH, Liu MM, Tay JWD, Teo SY, Kaewsapsak P, Jin S, et al. An engineered CRISPR-Cas12a variant and DNA-RNA hybrid guides enable robust and rapid COVID-19 testing. Nat Commun. 2021;12:1739.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Zhang L, Zuris JA, Viswanathan R, Edelstein JN, Turk R, Thommandru B, et al. AsCas12a ultra nuclease facilitates the rapid generation of therapeutic cell medicines. Nat Commun. 2021;12:3908.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Kleinstiver BP, Sousa AA, Walton RT, Tak YE, Hsu JY, Clement K, et al. Engineered CRISPR-Cas12a variants with increased activities and improved targeting ranges for gene, epigenetic and base editing. Nat Biotechnol. 2019;37:276–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Guo LY, Bian J, Davis AE, Liu P, Kempton HR, Zhang X, et al. Multiplexed genome regulation in vivo with hyper-efficient Cas12a. Nat Cell Biol. 2022;24:590–600.

  22. Ma H, Tu LC, Naseri A, Huisman M, Zhang S, Grunwald D, Pederson T. CRISPR-Cas9 nuclear dynamics and target recognition in living cells. J Cell Biol. 2016;214:529–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. McMahon MA, Prakash TP, Cleveland DW, Bennett CF, Rahdar M. Chemically modified Cpf1-CRISPR RNAs mediate efficient genome editing in mammalian cells. Mol Ther. 2018;26:1228–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Hendel A, Bak RO, Clark JT, Kennedy AB, Ryan DE, Roy S, et al. Chemically modified guide RNAs enhance CRISPR-Cas genome editing in human primary cells. Nat Biotechnol. 2015;33:985–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Yin H, Song CQ, Suresh S, Wu Q, Walsh S, Rhym LH, et al. Structure-guided chemical modification of guide RNA enables potent non-viral in vivo genome editing. Nat Biotechnol. 2017;35:1179–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Kristensen LS, Andersen MS, Stagsted LVW, Ebbesen KK, Hansen TB, Kjems J. The biogenesis, biology and characterization of circular RNAs. Nat Rev Genet. 2019;20:675–91.

    Article  CAS  PubMed  Google Scholar 

  27. Enuka Y, Lauriola M, Feldman ME, Sas-Chen A, Ulitsky I, Yarden Y. Circular RNAs are long-lived and display only minimal early alterations in response to a growth factor. Nucleic Acids Res. 2016;44:1370–83.

    Article  CAS  PubMed  Google Scholar 

  28. Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495:333–8.

    Article  CAS  PubMed  Google Scholar 

  29. Litke JL, Jaffrey SR. Highly efficient expression of circular RNA aptamers in cells using autocatalytic transcripts. Nat Biotechnol. 2019;37:667–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. 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 

  31. Katrekar D, Yen J, Xiang Y, Saha A, Meluzzi D, Savva Y, Mali P. Efficient in vitro and in vivo RNA editing via recruitment of endogenous ADARs using circular guide RNAs. Nat Biotechnol. 2022;40:938–45.

  32. Yi Z, Qu L, Tang H, Liu Z, Liu Y, Tian F, et al. Engineered circular ADAR-recruiting RNAs increase the efficiency and fidelity of RNA editing in vitro and in vivo. Nat Biotechnol. 2022;40:946–55.

  33. Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003;31:3406–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Filonov GS, Moon JD, Svensen N, Jaffrey SR. Broccoli: rapid selection of an RNA mimic of green fluorescent protein by fluorescence-based selection and directed evolution. J Am Chem Soc. 2014;136:16299–308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Paul CP, Good PD, Winer I, Engelke DR. Effective expression of small interfering RNA in human cells. Nat Biotechnol. 2002;20:505–8.

    Article  CAS  PubMed  Google Scholar 

  36. Lorenz R, Bernhart SH, Honer Zu Siederdissen C, Tafer H, Flamm C, Stadler PF, Hofacker IL. ViennaRNA Package 2.0. Algorithms Mol Biol. 2011;6:26.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Huang H, Hu Y, Huang G, Ma S, Feng J, Wang D, Lin Y, Zhou J, Rong Z. Tag-seq: a convenient and scalable method for genome-wide specificity assessment of CRISPR/Cas nucleases. Commun Biol. 2021;4:830.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Park HM, Liu H, Wu J, Chong A, Mackley V, Fellmann C. Extension of the crRNA enhances Cpf1 gene editing in vitro and in vivo. Nat Commun. 2018;9:3313.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Bin Moon S, Lee JM, Kang JG, Lee NE, Ha DI, Kim DY, et al. Highly efficient genome editing by CRISPR-Cpf1 using CRISPR RNA with a uridinylate-rich 3’-overhang. Nat Commun. 2018;9:3651.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Tong H, Huang J. High-fidelity Cas13 variants for targeted RNA degradation with minimal collateral effects. Nat Biotechnol. 2023;41:108–19.

    Article  CAS  PubMed  Google Scholar 

  41. Shi P, Murphy MR, Aparicio AO, Kesner JS, Fang Z, Chen Z, Trehan A, Guo Y, Wu X. Collateral activity of the CRISPR/RfxCas13d system in human cells. Commun Biol. 2023;6:334.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Ma S, Wang X, Hu Y, Lv J, Liu C, Liao K, et al. Enhancing site-specific DNA integration by a Cas9 nuclease fused with a DNA donor-binding domain. Nucleic Acids Res. 2020;48:10590–601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Ong CK, Subimerb C, Pairojkul C, Wongkham S, Cutcutache I, Yu W, et al. Exome sequencing of liver fluke-associated cholangiocarcinoma. Nat Genet. 2012;44:690–3.

    Article  CAS  PubMed  Google Scholar 

  44. van Haasteren J, Li J, Scheideler OJ, Murthy N, Schaffer DV. The delivery challenge: fulfilling the promise of therapeutic genome editing. Nat Biotechnol. 2020;38:845–55.

    Article  PubMed  Google Scholar 

  45. Fu L, Cao Y, Wu J, Peng Q, Nie Q, Xie X. UFold: fast and accurate RNA secondary structure prediction with deep learning. Nucleic Acids Res. 2022;50:e14–e14.

    Article  CAS  PubMed  Google Scholar 

  46. Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, et al. Highly accurate protein structure prediction with AlphaFold. Nature. 2021;596:583–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Ma S, Lv J, Sun J, Tang P, Li H, Zhou H, Zhang Z, Lin Y, Rong Z. iKA-CRISPR hESCs for inducible and multiplex orthogonal gene knockout and activation. FEBS Lett. 2018;592:2238–47.

    Article  CAS  PubMed  Google Scholar 

  48. Zhang X, Lv S, Luo Z, Hu Y, Peng X, Lv J, et al. MiniCAFE, a CRISPR/Cas9-based compact and potent transcriptional activator, elicits gene expression in vivo. Nucleic Acids Res. 2021;49:4171–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Zhang X, Wang W, Shan L, Han L, Ma SF, Zhang Y, Hao BT, Lin Y, Rong ZL. Gene activation in human cells using CRISPR/Cpf1-p300 and CRISPR/Cpf1-SunTag systems. Protein Cell. 2018;9:380–3.

    Article  CAS  PubMed  Google Scholar 

  50. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Zhang X, Wang X, Lv J, Huang H, Wang J, Zhuo M, et al. Effect of depletion of MEN1 or MLL on gene expression in the lung tissues of mice. PRJNA830337. BioProject. https://www.ncbi.nlm.nih.gov/bioproject/830337. (2022).

  52. Zhang X, Wang X, Lv J, Huang H, Wang J, Zhuo M, et al. Engineered circular guide RNAs boost CRISPR/Cas12a- and CRISPR/Cas13d-based DNA and RNA editing. https://doi.org/10.5281/zenodo.7991279. (2023).

  53. Zhang X, Wang X, Lv J, Huang H, Wang J, Zhuo M, et al. Engineered circular guide RNAs boost CRISPR/Cas12a- and CRISPR/Cas13d-based DNA and RNA editing. https://doi.org/10.5281/zenodo.7992742. (2023).

  54. Zhang X, Wang X, Lv J, Huang H, Wang J, Zhuo M, et al. Engineered circular guide RNAs boost CRISPR/Cas12a- and CRISPR/Cas13d-based DNA and RNA editing. https://doi.org/10.5281/zenodo.4679460. (2023).

  55. Zhang X, Wang X, Lv J, Huang H, Wang J, Zhuo M, et al. Engineered circular guide RNAs boost CRISPR/Cas12a and CRISPR/Cas13d based DNA and RNA editing. Figshare. https://doi.org/10.6084/m9.figshare.23268527. (2023).

Download references

Acknowledgements

We thank every member of Dr Zhili Rong’s lab, Dr Xiaoyang Zhao’s lab, and Dr Ying Lin’s lab for their helpful discussion and suggestions.

Review history

The review history is available as Additional file 3.

Peer review information

Tim Sands was the primary editor of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.

Funding

This work was funded by the National Key R&D Program of China (2022YFA0806300 to X.Z.2 (Xiaoyang Zhao) and Z.R.), the National Natural Science Foundation of China (82072329 to Y.L., 82070002 to Z.R., 82200072 to S.M., and 82102443 to T.H.), the GuangDong Basic and Applied Basic Research Foundation (2023A1515012269 and 2022A1515011091 to Y.L., 2021A1515110878 to S.M., 2022A1515111046 to X.Z.1 (Xin Zhang), and 2021B1515140031 to Z.R.), and the Fellowship of China Postdoctoral Science Foundation (2021M701634 to X.Z.1, 2022T150300 and 2021M691473 to S.M.)

Author information

Authors and Affiliations

Authors

Contributions

Z.R., X.Z.1 (Xin Zhang), X.Z.2 (Xiaoyang Zhao), and Y.L.1 (Ying Lin) conceived the idea, designed the experiments, analyzed the data, and wrote the manuscript. X.Z.1 performed most experiments. X.W. constructed plasmids and performed CasRx experiments. H.H. performed the Tag-seq assay. J.W. performed computer-assistant linker design and screening. Z.T. and Y.L.2 (Yuchen Liu) performed bioinformatics analysis for RNA-seq. J.L.1 (Jie Lv), M.Z., G.H., J.L.2 (Jiawei Liu), M.L., Q.L., and L.L. constructed plasmids. S.M. and T.H. performed some gene activation and genomic DNA cleavage experiments.

Ethics declarations

Ethics approval and consent to participate

All experiments involving mice were approved by the Institutional Animal Care and Use Committee of Southern Medical University (IACUC approval number: L2019018).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have 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:

 Fig. S1 Establishment of the dLbCas12a-p300knock-in HEK293T cell line. Fig. S2 Circular gRNAs increase thetranscription efficiency of LbCas12a-based activators. Fig. S3 Establishment of the LbCas12a knock-in HEK293Tcell line and off-target analysis of LbCas12a with different gRNAs. Fig. S4 Circular gRNAs improve the DNA cleavage efficiency of LbCas12ain MCF7 cells. Fig. S5 comparison between cgRNA with other extendedstructuralized gRNAs. Fig. S6 Circular gRNAs improve the gene expressionor DNA cleavage of AsCas12a-based effectors. Fig. S7 Multiplexedgene activation and cleavage guided by cgRNAs. Fig. S8 Efficient andspecific RNA cleavage activity of CasRX with cgRNA. Fig. S9 Thetrans-cleavage activity of CasRx–mediated exogenous transcripts degradation. Fig.S10 The trans-cleavage activity of CasRx–mediated endogenous transcriptsdegradation. Fig. S11 The structures of circular gRNAs for LbCas12awith different linkers targeting IL1RN predicted by mFold. Fig. S12The structures of circular gRNAs for LbCas12a with different linkerstargeting IL1RN predicted by mFold. Fig. S13 The structures ofcircular gRNAs for CasRx with different linkers targeting STAT3predicted by mFold.

Additional file 2:

 Table S1 The sequences of sgRNAs, primers and linkers.

Additional file 3.

 Peer review history.

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

Zhang, X., Wang, X., Lv, J. et al. Engineered circular guide RNAs boost CRISPR/Cas12a- and CRISPR/Cas13d-based DNA and RNA editing. Genome Biol 24, 145 (2023). https://doi.org/10.1186/s13059-023-02992-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13059-023-02992-z

Keywords