Skip to main content

CRISPR-Cas13d mediates robust RNA virus interference in plants

Abstract

Background

CRISPR-Cas systems endow bacterial and archaeal species with adaptive immunity mechanisms to fend off invading phages and foreign genetic elements. CRISPR-Cas9 has been harnessed to confer virus interference against DNA viruses in eukaryotes, including plants. In addition, CRISPR-Cas13 systems have been used to target RNA viruses and the transcriptome in mammalian and plant cells. Recently, CRISPR-Cas13a has been shown to confer modest interference against RNA viruses. Here, we characterized a set of different Cas13 variants to identify those with the most efficient, robust, and specific interference activities against RNA viruses in planta using Nicotiana benthamiana.

Results

Our data show that LwaCas13a, PspCas13b, and CasRx variants mediate high interference activities against RNA viruses in transient assays. Moreover, CasRx mediated robust interference in both transient and stable overexpression assays when compared to the other variants tested. CasRx targets either one virus alone or two RNA viruses simultaneously, with robust interference efficiencies. In addition, CasRx exhibits strong specificity against the target virus and does not exhibit collateral activity in planta.

Conclusions

Our data establish CasRx as the most robust Cas13 variant for RNA virus interference applications in planta and demonstrate its suitability for studying key questions relating to virus biology.

Introduction

Plant viruses are obligate parasites that rely mostly on host cells to complete their life (infection) cycle. They infect many economically important crops, leading to a significant decrease in crop yields worldwide. It is estimated that plant diseases cause 10 to 15% reduction in global crop yields each year, and 47% of this loss is caused by viruses [1,2,3]. Therefore, plant viruses threaten world agriculture and the food security of the rapidly growing world population. Most viruses infecting plants are RNA viruses, which comprise diverse groups and subgroups that are classified based on phylogenetic relationships determined by sequence homologies among the conserved virus genes, such as RNA-dependent RNA polymerase (RdRp), coat protein (CP), and movement protein (MP) [4]. Viruses with positive-sense, single-stranded RNA genomes represent most of these viruses. Therefore, we urgently need to develop effective antiviral agents to preserve crop yields from such viruses [5]. Various strategies, based on genome and genetic engineering approaches, have been developed to combat viruses, including directly inhibiting virus replication or targeting host factors required by viruses to replicate and persist [6]. RNA interference (RNAi), for example, is an innate antiviral immunity mechanism that has been successfully used to combat various plant viruses [7, 8]. Nevertheless, the availability of such antiviral strategies is still limited to specific virus groups, and many viruses, through evolution, mutate readily and have developed various counter-defense mechanisms, leading to rapid emergence of new viruses and nullifying the available antiviral approaches.

Bacterial and archaeal species employ CRISPR-Cas systems to fend off foreign genetic elements from invading phages and nucleic acids [9]. CRISPR-Cas adaptive immunity involves three main stages: (1) adaptation and spacer acquisition, where a piece of the invader genome is incorporated into the CRISPR array proximal to the 5′ leader sequence; (2) biogenesis of the CRISPR array into pre-CRISPR RNA (crRNA) to provide targeting specificity and processing into mature crRNA; and (3) interference, where the crRNA binds and guides the effector protein to the invader’s genome for cleavage or degradation [10]. These systems are divided into two classes: class I systems, which rely on multiple Cas proteins that mediate the target interference, and class II systems, which are represented by single, multidomain effector proteins. The class II CRISPR/Cas systems include types II, V, and VI. Types II and V comprise endonucleases that operate at the DNA level. On the other hand, type VI CRISPR/Cas systems are distinct in that they exclusively target RNA molecules [10, 11].

Because class II CRISPR-Cas systems are composed of a single effector to mediate the interference and immunity against the invader’s genomes, their engineering is simple; therefore, they are widely adopted in genome and transcriptome engineering. In recent years, class II CRISPR/Cas systems have been extensively developed to provide powerful and versatile tools for various genome engineering purposes [12,13,14]. Several studies have interrogated genomic and metagenomic sequencing data through the use of signature proteins to identify novel RNA-guided DNA and RNA endonucleases [11, 15]. Class II type II CRISPR-Cas9 targets DNA and has been used for genome engineering across eukaryotic species [16]. The impressive success of the DNA-targeting CRISPR/Cas systems in engineering and specifically editing the genome of different eukaryotic cells has led to the exploitation of this technology as a programmable antiviral defense strategy to confer resistance to many eukaryotic viruses, including human viruses [17,18,19]. In plants, CRISPR/Cas9 efficiently confers viral resistance to host plants. We and others have shown that CRISPR/Cas9 mediates strong interference against various DNA viruses. We have also shown that a single guide RNA targeting a conserved intergenic region in multiple Gemini viruses guides Cas9 to target different viruses simultaneously, providing interference against these DNA viruses [20,21,22,23]. These studies demonstrate the enormous potential of CRISPR/Cas9 as a promising strategy against plant DNA viruses [17]. Nevertheless, such DNA-targeting CRISPR modalities cannot target RNA molecules and therefore limit their applicability for targeted interference against RNA viruses.

Class II type VI CRISPR-Cas systems are RNA-guided and RNA-targeting machineries that provide prokaryotes with immunity against RNA [24] and DNA phages [25]. All type VI CRISPR-Cas systems have a single effector protein known as the Cas13 effector, formerly the C2c2 CRISPR-Cas system [24]. Cas13 effectors possess two distinct catalytic activities; one is RNase activity provided by the two higher eukaryotic and prokaryotic nucleotide binding domains (HEPN), which is required for target RNA degradation, and the other is RNase activity that catalyzes the processing and maturation of the pre-crRNA into mature crRNA [26, 27]. Intriguingly, these systems, once activated through the specific binding of the target RNA, confer collateral “promiscuous” activity and non-specifically degrade the RNA transcripts in the cell [24]. This collateral activity was harnessed to develop highly sensitive pathogen detection methods [28,29,30,31].

CRISPR/LshCas13a from Leptotrichia shahii was the first Cas13 orthologue to be harnessed for programmable RNA-targeting activities [24]. We previously developed LshCas13a to engineer immunity against the tobacco turnip mosaic RNA virus (TuMV) in different plant species, including Nicotiana benthamiana and Arabidopsis thaliana [32, 33]. We demonstrated that LshCas13a mediates specific RNA virus targeting in plants, albeit with moderate efficiency. In addition, we demonstrated that LshCas13a can process pre-crRNA to generate mature and functional crRNAs that guide LshCas13a to degrade the targeted virus. Subsequently, LshCas13a has become the prominent Cas13 orthologue in engineering immunity against different RNA viruses in different plant species, including monocot and dicot plants [34, 35]. Since the first report of LshCas13a, several studies have identified more variants of Cas13 proteins belonging to different Cas13 families, which have been classified into four type VI subtypes (A–D) [36,37,38]. These other variants have shown more robust catalytic activity and specificity compared to LshCas13a and have been harnessed for different RNA targeting and manipulations [28, 36, 38,39,40]. In addition, these studies have shown that these Cas13 variants exhibit variable efficiencies for transcriptome targeting and engineering. LwaCas13a, for example, has been reported to mediate more robust RNA-targeting activity than the LshCas13 system but requires a stabilizer fusion, for example, msfGFP for efficient interference activity [38]. Subsequent studies have identified PspCas13b to be more efficient for RNA targeting in mammalian cells compared to the LwaCas13a protein, and PspCas13b does not require a stabilizer protein for its activity [39]. Recently, a new and additional Cas13 subtype has been identified, type VI-D, which shows minimal sequence identity with previous Cas13 effectors. Interestingly, Cas13d effector proteins are much smaller than all previously reported subtypes (median size ~ 300 amino acids, 26% smaller than other Cas13 proteins) [36, 40] and similar to type VI-B; the majority of the type VI-D orthologues contain associated accessory proteins harboring WYL domains [36]. Notably, experimental characterization of the Cas13d activity of the Ruminococcus flavefaciens (CasRx) orthologue has shown consistently robust activity and specificity in mammalian cells compared to the previously characterized Cas13 proteins, namely LwaCas13a and PspCas13b [40]. These enzymes will enable diverse RNA manipulations for different purposes.

Because most viruses infecting plants possess RNA genomes, and reasoning that other uncharacterized Cas13 proteins might exhibit more robust activity in plants against RNA viruses, we sought to identify a more robust RNA-targeting CRISPR/Cas13 system by characterizing a set of Cas13 family members to assess their virus RNA-targeting activity in planta. Here, we investigated nine different Cas13 variants for robust virus interference to establish the use of this adaptive immunity mechanism for engineering plant immunity against single and multiple RNA viruses. We identified CasRx as the most effective Cas13 orthologue for RNA virus targeting in N. benthamiana. CasRx confers robust RNA virus interference compared to other variants in transient assays targeting a highly replicating RNA virus. In addition, transgenic plants overexpressing CasRx exhibited efficient interference against TuMV RNA compared to plants overexpressing other Cas13 variants, including LwaCas13a and PspCas13b, using cognate, position-matched crRNAs. Importantly, the catalytically deactivated (dCasRx) version could not mediate any virus targeting and, thus, inhibition, indicating that the HEPN-dependent catalytic activity of CasRx is required for virus RNA degradation. Finally, using synchronous coinfection of two non-competing RNA viruses, we show that CasRx RNA virus targeting is highly specific to the targeted virus and lacks collateral activities in planta. In addition, CasRx mediated efficient multiplexed virus interference by targeting two different RNA viruses simultaneously. This work establishes the use of CasRx for efficient RNA virus interference applications in plants.

Results

Construction of different CRISPR-Cas13 variants for expression in planta

Different CRISPR-Cas13 variants have been identified through computational pipelines searching for effector proteins possessing two HEPN domains. Cas13a, Cas13b, and Cas13d have been identified and tested for RNA manipulation in vitro and in vivo [36,37,38, 40,41,42]. Here, we investigated different Cas13 orthologues to assess their catalytic activities for RNA virus interference in plants. To this end, we generated binary vectors for in planta expression of multiple Cas13 proteins of three different Cas13 subtypes, including subtypes A, B, and D (Fig. 1a). In addition to the LshCas13a that we previously employed [32, 33], we selected one orthologue of Cas13a from Leptotrichia wadei (LwaCas13a), two orthologues of Cas13b from Bergeyella zoohelcum (BzCas13b) and from Prevotella sp. P5-125 (PspCas13b), and one orthologue of Cas13d from Ruminococcus flavefaciens XPD3002 (CasRx). The selection of these different Cas13 orthologues was based on several reports showing improved RNA-targeting activity of these variants in comparison to LshCas13a protein in vitro and in mammalian cells [37, 38, 40, 42]. Different fusions were made to these Cas13 variants to improve their stability or control their subcellular localization for optimum activity against the target transcripts. The different Cas13 protein open reading frames (ORFs) were fused to either a nuclear localization signal (NLS) or a nuclear export signal (NES), or not fused to any localization signal. In addition, all proteins were fused with an HA tag to facilitate protein detection (Fig. 1b). The different Cas13 orthologues were cloned into the pK2GW7 binary vector, and their expression was driven by the 35S cauliflower mosaic virus promoter. Because the PspCas13b orthologue has been previously shown to be most active for RNA interference with the NES signal [39], we proceeded with the PspCas13b-NES module without considering an NLS variant. For expression of Cas13’s cognate crRNAs conferring specificity against the green fluorescent protein (GFP) reporter gene in the target virus or against a virus protein-coding sequence, we employed our tobacco rattle virus (TRV)-based system to transiently and systemically express the crRNAs in plants [43], where the corresponding crRNAs were cloned under the pea early browning virus promoter, pPEBV, in the TRV-RNA2 genome.

Fig. 1
figure 1

CRISPR/Cas13 orthologues for RNA virus interference. a Schematic representation of CRISPR/Cas13 variants and their respective crRNA structures. Structural representations of different Cas13 subtypes with their corresponding crRNA structures are shown with the estimated average size of the Cas13 protein under each subtype. HEPN, high eukaryotic and prokaryotic nucleotide binding domains; aa, amino acid; nt, nucleotides. b Schematic of different Cas13 protein variants used in this study. The different Cas13 variants are shown with their different fusions. NLS, nuclear localization signal; NES, nuclear export signal; GFP, green fluorescent protein. c Schematic of TRBO-based Cas13 mediated RNA virus interference. The highly replicating plant TRBO RNA virus expressing GFP protein was used as a reporter system to screen for efficient Cas13 activity in transient assays via crRNA delivered through TRV system, targeting two different regions within the GFP sequence of the virus RNA genome. 35S, CaMV 35S promoter; pPEBV, pea early browning virus promoter; RNA1/RNA2, genomes of the tobacco rattle virus (TRV)

Screening of Cas13 variants for efficient RNA virus interference in transient assays

To assess the CRISPR/Cas13 activity in planta, generation of transgenic plants overexpressing the CRISPR/Cas13 machinery is normally needed. However, such an approach is time-consuming and challenging when multiple and different CRISPR/Cas13 variants need to be evaluated, especially when considering the evaluation of multiple and different crRNAs for each variant. Therefore, the development of a transient assay that allows the assessment of different Cas13 variants and crRNAs simultaneously and systemically prior to the generation of transgenic plants would be valuable. We adopted a tobacco mosaic virus (TMV)-RNA-based overexpression (TRBO-G) system expressing the GFP gene to serve as a reporter construct to measure Cas13 interference activity in transient assays in plants. TRBO-G is a positive-sense, single-stranded RNA plant-infecting virus that has been engineered to replace the CP-encoding sequence with a GFP-encoding sequence, rendering the virus unable to move systemically in the infected plants, but retaining its ability to highly and efficiently replicate and express GFP protein in the infected (infiltrated) leaves (Additional file 1: Figure S1) [44]. This reporter system was employed to provide a unique and easy-to-monitor system for virus infection and replication, thereby facilitating the assessment and comparison of the interference activity of the different CRISPR/Cas13 systems against RNA viruses in plants (Fig. 1c).

To assay for the interference activity of the different Cas13 variants, we co-delivered each of the Cas13 variant with the position-matched Cas13’s cognate crRNAs targeting GFP coding sequence within the TRBO-G genome, and the TRBO-G-expressing construct into leaves of wild type N. benthamiana plants via Agro-infiltration. We used a non-specific crRNA with no complementarity to the TRBO-G genome as a negative control. We tested two position-matched crRNAs (GFP-T1 and GFP-T2) with each of the Cas13 variants and measured the virus-expressed GFP signal 3 days post-infiltration. Notably, we observed a broad range in the interference levels between the different Cas13 variants. Our data show that both CasRx variants (CasRx-NLS and CasRx-NES) exhibited the highest level of interference against TRBO-GFP using GFP-T1 and GFP-T2 crRNAs compared to the other Cas13 orthologues. In addition, other variants including LwaCas13a-NLS, LwaCas13a-NES, and PspCas13p-NES exhibited an efficient interference against TRBO-GFP using also the two GFP-targeting crRNAs (Fig. 2a). In addition, as another control, we conducted one set of experiments in which we delivered the crRNAs and the target virus (TRBO-GFP) in the inoculated leaves, but without the Cas13 protein variants. We did not observe any virus interference, indicating that the virus interference observed in Fig. 2a is exclusively Cas13-crRNA dependent (Additional file 1: Figure S2). We repeated these transient assays in at least three independent experiments using the GFP signal as an indicator of RNA virus interference, and our data showed that CasRx is highly efficient compared to other Cas13 variants and NS-crRNA controls (Fig. 2b).

Fig. 2
figure 2

TRBO-GFP-based screening of different Cas13 orthologues for efficient virus interference. a GFP monitoring to assess the Cas13-mediated virus interference activities in Agro-infiltrated wild type N. benthamiana leaves in transient assays. Images were taken 3 days post-infiltration. NS, non-specific crRNA. b Relative fluorescence intensity quantification of leaf images. For each Cas13 variant, GFP signal intensity of each targeting crRNA (GFP-T1 and GFP-T2) is shown relative to the non-targeting (NS) crRNA. Values shown as mean ± SEM (n = 3). c Western blot analysis of the abundance of the virus expressed GFP protein to confirm the Cas13-mediated TRBO-GFP virus interference. Protein blots were developed with anti-GFP antibody. α-GFP, anti-GFP antibody. Ponceau staining served as loading control. d Quantification of the western blot data. For each Cas13 variant, the abundance of the GFP protein with the targeting crRNAs (GFP-T1 and GFP-T2) is shown relative to the non-targeting (NS) crRNA. Error bars indicate SEM (n = 3). e RT-qPCR analysis of TRBO-GFP knockdown with different Cas13 variants using the two position-matched crRNAs. For each Cas13 variant, knockdown efficiency of each targeting crRNA (GFP-T1 and GFP-T2) is shown relative to the non-targeting (NS) crRNA. Values shown as mean ± SEM (n = 3)

Next, we conducted molecular analyses to assess the RNA virus interference via immunoblotting of the GFP protein of TRBO-GFP virus to validate the observed reduction in the GFP signal. Our immunoblotting data corroborated our phenotypic data (GFP signal in inoculated leaves) and showed that LwaCas13a-NLS, LwaCas13a-NES, and PspCas13b-NES mediated strong interference against the RNA virus. Intriguingly, CasRx mediated the highest and most efficient interference when compared to these or other Cas13 variants, as indicated with the low GFP protein level of the targeted virus (Fig. 2c, d).

Because the CRISPR-Cas13 machinery targets the RNA transcripts of the virus for degradation, we wanted to assess the interference at the RNA level. Therefore, we conducted reverse transcription-quantitative PCR (RT-qPCR) to assess the interference activities of different Cas13 variants. Consistent with the level of GFP signal, the RT-qPCR data showed that CasRx mediated robust and highly efficient RNA interference when compared to other Cas13 variants (Fig. 2e). These data corroborated our phenotypic and immunoblotting data and identified CasRx as the most efficient Cas13 variant for RNA interference applications.

CasRx outperforms LwaCas13a and PspCas13b variants in targeting the virus replicase for RNA interference

In our previous interference experiments using crRNA against the GFP sequence of the virus genome, we identified three Cas13 orthologues with the highest level of interference activity against the TRBO-GFP RNA genome, namely LwaCas13a, PspCas13b, and CasRx. However, CasRx mediated the highest interference efficiency among these variants. Next, we wanted to test whether this interference activity and efficiency are applicable to other virus sequences and whether targeting other (conserved) virus sequences would result in similar interference efficiencies. Therefore, we designed three crRNAs targeting the replicase gene of the virus (Rep-T1, Rep-T2, and Rep-T3) and tested the activity of LwaCas13a, PspCas13b, and CasRx variants for RNA virus interference. Our data showed that CasRx mediated the highest level of interference activity, followed by LwaCas13a proteins (Fig. 3a). Interestingly, and despite its strong interference activity against the GFP sequence of the virus, PspCas13b did not mediate an efficient interference against the replicase transcripts using the three crRNAs (Fig. 3a). We also found that different crRNAs targeting different regions of the replicase transcripts exhibited variable interference efficiencies, where Rep-T3 crRNA mediated better interference relative to the other crRNAs (Fig. 3a), probably due to different factors affecting target accessibility [24]. Strikingly, CasRx mediated strong interference with all crRNAs used to target the replicase transcripts, with CasRx-NES mediating better interference compared to the CasRx-NLS variant (Fig. 3a). Independent repetition of these experiments against the virus replicase transcripts showed similar and consistent results, and confirmed that CasRx mediated robust RNA virus interference efficiencies when compared to other Cas13 variants (Fig. 3b).

Fig. 3
figure 3

Characterization of LwaCas13a, PspCas13b, and CasRx activity against essential and conserved genomic region of RNA viruses. a GFP monitoring to assess the Cas13-mediated virus interference activities in Agro-infiltrated wild type N. benthamiana leaves in transient assays. Images were taken 3 days post-infiltration. NS, non-specific crRNA; Rep, replicase. b Relative fluorescence intensity quantification of leaf images. For each Cas13 variant, GFP signal intensity of each targeting crRNA (Rep-T1, Rep-T2, and Rep-T3) is shown relative to the non-targeting (NS) crRNA. Values shown as mean ± SEM (n = 3). c Western blot analysis of the abundance of the virus expressed GFP protein to confirm the Cas13-mediated TRBO-GFP virus interference. Protein blots were developed with anti-GFP and anti-HA antibodies. α-GFP, anti-GFP antibody, α-HA, anti-HA antibody. Ponceau staining served as loading control. d Quantification of the western blot data. For each Cas13 variant, the abundance of the GFP protein with the targeting crRNAs (Rep-T1, Rep-T2, and Rep-T3) is shown relative to the non-targeting (NS) crRNA. Error bars indicate SEM (n = 3). e RT-qPCR analysis of TRBO-GFP knockdown with different Cas13 variants using the three position-matched crRNAs. Transcript levels are shown relative to leaves inoculated with only TRBO-GFP vector. Values shown as mean ± SEM (n = 3)

Subsequently, we conducted molecular analyses using immunoblotting to assess the virus interference activities of these Cas13 variants. The immunoblotting data showed that targeting the Rep-T3 region within the virus replicase sequence provided strong interference with both LwaCas13a variants relative to the other two crRNAs, with improved efficiency with the cytoplasmic variant of LwaCas13a (LwaCas13a-NES) as evident with the Rep-T2 crRNA (Fig. 3c). Consistent with the phenotypic results of the GFP signal, no reduction in the GFP protein level was detected with the PspCas13b-NES protein in comparison to the NS control, with exception of the Rep-T1 crRNA that was better than other crRNAs (Fig. 3c). Analyzing the interference activities of CasRx showed that CasRx mediated highly efficient RNA interference using the three crRNAs targeting the replicase sequences. CasRx-NLS was efficient with the three crRNAs, but the highest efficiency was observed when Rep-T3 crRNA was used (Fig. 3c). Notably, CasRx-NES was efficient in conferring RNA virus interference against the RNA virus using all three crRNAs (Fig. 3c). Of note, although the protein levels of CasRx-NES were much lower compared to CasRx-NLS, the interference efficiencies of CasRx-NES were much stronger and more robust, indicating the importance of the localization of the Cas13 variant for virus interference activities (Fig. 3c). We repeated these experiments three times with similar results, and our combined data showed that CasRx mediated robust RNA virus interference using three crRNAs against the replicase transcripts (Fig. 3d).

To further validate these observations, RT-qPCR was performed to detect the accumulation of the TRBO-GFP genome. Consistent with the phenotypic and immunoblotting data, RT-qPCR results showed a significant reduction in the accumulation of the TRBO-GFP genomic RNA with Rep-T3 crRNA of both LwaCas13a and CasRx variants. Notably, CasRx was able to mediate efficient virus targeting with all three crRNAs in comparison with LwaCas13a or PspCas13b variants, with increased efficiency observed with the CasRx-NES variant (Fig. 3e). These data indicate that CasRx is capable of mediating robust and consistent efficiency in targeting different genomic sequences of the virus genome.

CasRx confers highly efficient interference against TuMV

Next, we tested whether the interference activities of the Cas13 variants conducted in transient assays in inoculated leaves were similar when using a different virus capable of systemic spread throughout the plant. To this end, we conducted a set of transient assays in wild-type N. benthamiana plants using the TuMV-GFP virus as our target. We co-delivered the Cas13 variants (LwaCas13a-NLS, LwasCas13a-NES, PspCas13b-NES, CasRx-NLS, and CasRx-NES) identified from our previous experiments to confer strong RNA virus interference activities, with their cognate crRNAs targeting different genomic sequences of virus genome, and TuMV-GFP as a target. Cas13-mediated interference against the TuMV-GFP genome would result in attenuated replication and, thus, spread of the virus, which can be measured by monitoring the spread of the virus-expressed GFP to systemic leaves [33]. Our data showed that LwaCas13a-NES mediated strong interference against TuMV-GFP by limiting its systemic spread when HC-Pro and GFP-T2 crRNAs were used (Additional file 1: Figure S3). In addition, LwaCas13a-NLS mediated interference activities with the same crRNAs, but to a lower efficiency when compared to the LwaCas13a-NES variant. Moreover, PspCas13b-NES mediated modest interference when compared to LwaCas13a-NES with these position-matched crRNAs (Additional file 1: Figure S3). The interference activities were variable based on the crRNA that was used, with HC-Pro and GFP-T2 crRNAs providing the highest activity in conferring virus interference. Relative to the other tested Cas13 variants, CasRx exhibited the strongest and most efficient interference activities with CasRx-NLS and CasRx-NES using all crRNAs (Additional file 1: Figure S3). However, CP crRNA was the least efficient when compared to other crRNAs, and CasRx-NES was better than CasRx-NLS with all crRNAs tested, corroborating our previous data and indicating that CasRx was the most efficient Cas13 variant for RNA virus interference (Additional file 1: Figure S3).

Next, we tested whether these interference activities would be observed in plants stably expressing Cas13 variants, compared to the transient expression experiments performed above. To this end, we generated stable expression lines of N. benthamiana plants expressing Cas13 variants, including LwaCas13a-NLS, LwaCas13a-NES, PspCas13b-NES, and CasRx-NLS, using Agrobacterium tumefaciens-mediated transformation to assess the interference efficiencies of these Cas13 variants. We confirmed the expression of Cas13 proteins in the generated permanent lines, and our western blotting analysis show that the different Cas13 proteins were expressed efficiently in N. benthamiana plants, and the correct proteins sizes were detected (Additional file 1: Figure S4). Using Cas13 overexpressing lines with relatively similar level of protein expression, we conducted our RNA interference experiments using TuMV-GFP as the target virus and delivered different crRNAs (GFP-T1, GFP-T2, HC-Pro, and CP, as well as NS-crRNA) via the TRV system. All these Cas13 variants exhibited RNA virus interference, but with different efficiencies among the different crRNAs (Fig. 4a). Consistent with our previous data, LwaCas13a-NES was better with all crRNA tested than LwaCas13a-NLS, and both of these variants exhibited better virus interference than PspCas13b-NES (Fig. 4a). Most importantly, we found that the CasRx variant resulted in the highest level of virus interference using different crRNAs compared to other Cas13 variants (Fig. 4a).

Fig. 4
figure 4

CasRx mediates efficient interference against TuMV-GFP virus by preventing its systemic spread. a Cas13-mediated interference against the GFP-expressing TuMV virus in plants. N. benthamiana plants expressing LwaCas13a, PspCas13b, and CasRx proteins were co-infiltrated with TRV (expressing crRNAs targeting the TuMV-GFP virus and NS crRNA) and TuMV-GFP. At 7 dpi, plants were imaged under UV light for GFP signal monitoring. b RT-qPCR analysis of TuMV-GFP knockdown with different Cas13 variants using the different position-matched crRNAs. For each Cas13 variant, knockdown efficiency of each targeting crRNA is shown relative to the non-targeting (NS) crRNA. Values shown as mean ± SEM (n = 3)

We next investigated the virus interference activities of these Cas13 variants at the molecular level. Therefore, we conducted RT-qPCR assays to quantify the RNA message of the virus in the systemic leaves of the infected plants of different Cas13 variants and crRNAs. Consistent with the TuMV-GFP signal observed in the phenotypic data, our RT-qPCR results showed that CasRx mediated the highest interference efficiency in conferring RNA virus interference across all crRNAs used compared to other variants, with comparable knockdown efficiency of GFP-T1 crRNA between LwaCas13a-NES and CasRx-NLS variants (Fig. 4b). Therefore, we concluded that CasRx conferred much higher efficiency interference compared to other Cas13 variants.

CasRx catalytic activity is required for RNA virus interference

Our RNA virus interference transient and stable overexpression assays indicated that CasRx mediated high efficient virus interference with different crRNAs targeting TRBO-GFP or TuMV-GFP viruses. We wanted to investigate whether the catalytic activity of CasRx is required to confer interference against RNA viruses. It is possible that binding of CasRx to the virus genome may attenuate its replication and translation, leading to interference activities mediated by binding, but not the catalytic degradation of the RNA transcripts [45]. Therefore, we used a catalytically deactivated CasRx (dCasRx) variant harboring mutations in the catalytic HEPN domains of the CasRx-NLS protein (R239A/H244A/R858A/H863A). These mutations have been reported to abrogate the catalytic cleavage activity but not the binding activity of the CasRx protein [40]. Subsequently, we co-delivered the catalytically active variants of CasRx (CasRx-NLS and CasRx-NES) and the catalytically deactivated variant of CasRx (dCasRx) into the leaves of N. benthamiana plants with TRBO-GFP virus as a target as well as GFP-T1, GFP-T2, and NS crRNAs. Our data revealed that only the catalytically active variants of CasRx mediated virus interference and that the dCasRx variant did not show any interference activities (Fig. 5a, b).

Fig. 5
figure 5

CasRx catalytic activity is required for RNA virus interference. a Illustration of the targeting activity of the catalytically active and catalytically deactivated (dCasRx) CasRx variants against the targeted TRBO-GFP virus. b GFP monitoring to assess the Cas13-mediated virus interference activities in Agro-infiltrated wild type N. benthamiana leaves in transient assays. Images were taken 3 days post-infiltration. NS, non-specific crRNA. c Western blot analysis of the abundance of the virus expressed GFP protein to confirm the Cas13-mediated TRBO-GFP virus interference. Protein blots were developed with anti-GFP and anti-HA antibodies. α-GFP, anti-GFP antibody, α-HA, anti-HA antibody. Ponceau staining served as loading control. d RT-qPCR analysis of TRBO-GFP knockdown with the catalytically active and catalytically deactivated CasRx variants. For each CasRx variant, knockdown efficiency of each targeting crRNA (GFP-T1 and GFP-T2) is shown relative to the non-targeting (NS) crRNA. Values shown as mean ± SEM (n = 3)

To validate this observation, we conducted immunoblotting and RT-qPCR assays to detect and quantify the virus-expressed GFP protein as well as the accumulation of the virus RNA genome with either the catalytically active CasRx-NLS and CasRx-NES variants or the catalytically deactivated (dCasRx) variant. Consistent with the phenotypic data, we observed a significant reduction in the abundance of the virus GFP protein with the catalytically active CasRx variants compared to the dCasRx (Fig. 5c). In addition, a significant knockdown of the targeted RNA virus was observed with the catalytically active CasRx variants with both crRNAs, while no decrease in the RNA level of the targeted TRBO-GFP virus genome was observed with dCasRx (Fig. 5d). These data indicated that CasRx-mediated virus interference is dependent on the catalytic activity of CasRx and suggesting that binding of dCasRx to the virus genome does not inhibit its replication.

CasRx exhibits multiplexed and specific virus-targeting activity with no observed collateral activity

We have shown that CasRx is a robust Cas13 variant for RNA virus interference applications in planta. Specificity of CasRx is an important issue in its interference applications against viruses. Therefore, we sought to test the specificity of CasRx against different RNA viruses. We developed a virus-based reporter system that relies on synchronous coinfection and co-replication of two different viral vectors. The two vectors need to be derived from different RNA viruses so they can non-competitively replicate and exist in the same cell [46]. Therefore, besides our TRBO RNA virus, we adopted chimeric potato virus X expressing GFP (PVX-GFP), which has been found to be non-competing with the TMV-based vectors [46]. To build an easy-to-monitor system to determine CasRx specificity, we replaced the GFP sequence in our TRBO-GFP chimeric virus with an enhanced blue fluorescent protein (EBFP)-encoding sequence, resulting in the generation of a TRBO-BFP construct. We tested the engineered TRBO-BFP construct and detected a blue fluorescent signal instead of the GFP signal, indicating that our modification was successful. Subsequently, we co-delivered CasRx-NES and CasRx-NLS variants, PVX-GFP and TRBO-BFP, into the leaves of N. benthamiana with crRNAs providing specificity to either virus. Theoretically, by delivering virus-specific crRNAs, only the titer of the targeted virus will be mitigated, indicating specific interference against the targeted virus. However, if the interference activity is not specific, interference of both viruses will occur (Fig. 6a). Our data showed that using CasRx variants to target the TRBO-BFP virus, the interference activity occurred against the TRBO-BFP virus alone, but the co-delivered PVX-GFP signal was high, indicating the absence of interference activity against the non-targeted virus (Fig. 6b, c). Similarly, when we targeted the PVX-GFP virus, only the GFP signal was mitigated, and not the BFP signal generated by TRBO-BFP (Fig. 6d, e). Quantification of the viruses’ genomic RNA by RT-qPCR showed robust knockdown of the targeted viruses, whereas no change in the abundance of the non-targeted virus or viruses with non-specific (NS) crRNA was observed, indicating the specificity of the targeted virus interference (Fig. 6). Next, we assessed the capability of CasRx to mediate multiplexed knockdown of different RNA viruses. Therefore, we co-delivered each of the catalytically active CasRx variants with both PVX-GFP and TRBO-BFP viruses into the leaves of N. benthamiana with simultaneous delivery of two different crRNAs, each providing specificity to either virus (GFP-T2 against PVX-GFP, and Rep-T3 against TRBO-BFP virus). Our data show that both CasRx variants were able to exhibit efficient targeting of both viruses simultaneously using the two (GFP-T2 and Rep-T3) crRNAs, while no virus interference was observed when NS crRNA was used (Additional file 1: Figure S5a, S5b). These results indicate the high specificity of CasRx in targeting RNA viruses, and that CasRx is amenable to multiplex targeting of different RNA viruses in plants.

Fig. 6
figure 6

Specific virus targeting via CasRx with no observed collateral activity. a Schematic of CasRx targeting of TRBO-BFP and PVX-GFP viruses in synchronous co-infected leaves with the two different viral vectors, and the possible outcomes. b BFP monitoring to assess the Cas13-mediated virus interference activities in Agro-infiltrated wild type N. benthamiana leaves in transient assays (left). Images were taken 3 days post-infiltration. NS, non-specific crRNA; Rep, replicase. RT-qPCR analysis of TRBO-BFP knockdown with Rep-T3 crRNA used against TRBO-BFP virus (right). The transcript levels of TRBO-BFP are shown relative to the transcript level of TRBO-BFP virus targeted with (NS) crRNA. Values shown as mean ± SEM (n = 3). c BFP and GFP monitoring to assess the Cas13-mediated virus interference activities in Agro-infiltrated wild type N. benthamiana leaves in transient assays (left). Images were taken 3 days post-infiltration. NS, non-specific crRNA; Rep, replicase. RT-qPCR analysis of TRBO-BFP and PVX-GFP knockdown with Rep-T3 crRNA (right). The transcript levels of TRBO-BFP are shown relative to the transcript level of the non-targeted virus (PVX-GFP), and the transcript level of both viruses is shown relative to the transcript level of both viruses targeted with (NS) crRNA. Values shown as mean ± SEM (n = 3). d GFP monitoring to assess the Cas13-mediated virus interference activities in Agro-infiltrated wild type N. benthamiana leaves in transient assays (left). Images were taken 3 days post-infiltration. NS, non-specific crRNA. RT-qPCR analysis of PVX-GFP knockdown with GFP-T2 crRNA used against PVX-GFP virus (right). The transcript levels of PVX-GFP are shown relative to the transcript level of PVX-GFP virus targeted with (NS) crRNA. Values shown as mean ± SEM (n = 3). e BFP and GFP monitoring to assess the Cas13-mediated virus interference activities in Agro-infiltrated wild type N. benthamiana leaves in transient assays (left). Images were taken 3 days post-infiltration. NS, non-specific crRNA. RT-qPCR analysis of TRBO-BFP and PVX-GFP knockdown with GFP-T2 crRNA (right). The transcript levels of PVX-GFP are shown relative to the transcript level of the non-targeted virus (TRBO-BFP), and the transcript level of both viruses is shown relative to the transcript level of both viruses targeted with (NS) crRNA. Values shown as mean ± SEM (n = 3)

Discussion

Prokaryotic adaptive immunity systems have been harnessed for engineering genomes and transcriptomes across eukaryotic species [47]. Different class II type VI CRISPR-Cas13 systems have been employed for transcriptome engineering and nucleic acid detection in mammalian cells. Despite the fact that several Cas13 variants have been reported with variable crRNA scaffolds, sizes, and importantly catalytic activities, CRISPR-LshCas13a has been the only Cas13 system used for RNA virus interference activities in plants. Here, we characterized different Cas13 variants for robust interference against RNA viruses in transient assays and stable N. benthamiana Cas13-overexpression plants. We employed RNA virus interference activity reporter systems and conducted a comprehensive screen to identify the most catalytically active Cas13 variant in conferring RNA interference. Our data revealed that CasRx exhibits the most robust catalytic activity among the variants tested.

In this study, we constructed different Cas13 variants, including Cas13a (LshCas13a and LwaCas13a), Cas13b (PspCas13b and BzCas13b), and Cas13d (CasRx), and their corresponding crRNA scaffolds for in planta expression. Transient assays using TRBO-GFP as the interference target showed that CasRx was the most efficient in conferring interference using GFP-T1 and GFP-T2 crRNAs. In addition, LwaCas13a and PspCas13b exhibited strong efficiencies compared to LshCas13a and BzCas13b variants. Interestingly, although we used position-matched crRNAs with all the tested Cas13 proteins, these crRNAs exhibited variable efficiencies among the Cas13 variants. This might result from the diverse molecular behaviors, molecular architectures, and protospacer flanking sequence (PFS) requirements for these different Cas13 variants [48]. We selected the top three variants for further analysis and tested them in targeting a conserved virus sequence of the replicase transcript of TMV. Our data showed that CasRx exhibited robust virus interference using the three crRNAs that were designed against the replicase sequence. However, other variants exhibited interference, and LwaCas13a was more efficient than PspCas13b. The fact that CasRx targeted the different regions of the replicase sequence with high efficiency indicates the competency and suitability of this system for efficient targeting of other essential sequences in other RNA viruses, and offers the potential development of this system for further RNA virus targeting and manipulation. Molecular analysis corroborated our phenotypic data and indicated that in general, cytoplasmic-localized Cas13 (Cas13-NES) provided better interference than the NLS-localized variants (Cas13-NLS), which could be explained by the simultaneous localization of the cytoplasmic-localized Cas13 proteins and the targeted virus in the plant cytoplasm, where the targeted viruses replicate.

To test whether these variants limit systemic spread of the virus with similar efficiencies, we used TuMV-GFP as a target and co-delivered the Cas13 variants with their corresponding crRNAs. CasRx provided robust interference against the TuMV virus, thus significantly limited the spread of the TuMV-GFP virus into systemic leaves. In addition, the cytoplasmic-localized LwaCas13a provided strong interference compared to PspCas13b. To corroborate these data, we generated N. benthamiana overexpression lines of these Cas13 variants and conducted interference assays. Our data showed that CasRx is indeed the most robust Cas13 variant to confer virus interference followed by the LwaCas13a-NES variant. These data indicate that CasRx is an excellent candidate for RNA virus interference applications in plants.

Previous studies harnessing CRISPR-Cas proteins, such as FnCas9 from Francisella novicia, to target RNA viruses have shown that the binding, but not the cleavage, ability of these proteins is sufficient to inhibit virus replication in mammalian and plant cells [45, 49]. Our data show that the HEPN-dependent catalytic activity of CasRx is essential for efficient virus interference, which in turn indicate that the observed CasRx-mediated virus interference results from targeted viral genomic degradation. However, the ability of targeting the dCasRx to specifically bind to mRNA transcripts may have various uses in RNA transcript manipulations. For example, fusion of dCasRx to a splicing effector has been used to manipulate alternative splicing in mammalian cells [40]. Similarly, development of dCasRx as a specific and programmable RNA-binding platform in plants would offer various applications, including virus transcript tracking and imaging to study the biology of virus replication. Furthermore, the generation of chimeric fusions of dCasRx and other functional domains may be beneficial for developing CasRx-based regulatory reagents for imaging of endogenous mRNA and non-coding RNA, fine-tuning and regulation of alternative splicing, tracking and trafficking, translational control, and regulation under diverse growth, environmental, and developmental conditions and cues [50, 51].

The specificity of CasRx for in vivo application is an important point to consider. Therefore, we evaluated the specificity of CasRx in the context of virus targeting in planta, and our data show that the CasRx targeting activity is highly specific to the targeted virus, where the transcript abundance of the co-delivered, non-targeted, virus was not affected. Cas13 variants possess robust collateral activities in vitro, enabling the development of pathogen detection systems that have revolutionized molecular diagnostics [52, 53]. Although the collateral activities of Cas13 variants are an essential part of the broader immunity conferred by CRISPR-Cas13 systems in prokaryotes [25], such activity is not observed in human and plant cells [33, 38]. Here, we could infer from these results that CasRx, and probably other Cas13 variants, lack such collateral activity in plants. Nevertheless, further experiments and analysis are needed to corroborate this conclusion. It remains to be tested whether accessory proteins, including CXS28 and WYL, can stimulate the activity of CasRx against the targeted transcripts and induce collateral activities in planta. This interference system can further be improved by determining the molecular basis of how the type VI system mediates antiviral defenses within their own hosts (spacer acquisition and adaptation) and how these effectors might function in eukaryotic hosts. These questions warrant detailed studies to understand the molecular functions, and the basis of the catalytic activities, of these proteins for broader applications in endowing immunity to eukaryotes.

In addition, we have shown that CasRx can be used to target either one virus alone or two RNA viruses simultaneously, thereby extending the applicability of this interference system. Such system might have the potential use in crop species for virus interference of mixed infections of RNA viruses that occur in the field and under natural conditions. It is worth noting that recent studies have shown the ability of DNA viruses to evolve as a result of Cas9 targeting [21, 54]. Therefore, further studies are needed to test the ability of RNA viruses to evolve or escape the CRISPR Cas13 machinery in plants. Interestingly, while preparing this manuscript, a recent study harnessed different CRISPR-Cas13 systems, including LwaCas13a and PspCas13b, for interference against several human RNA viruses. Importantly, the study showed that Cas13 targeting of the virus genome does not result in mutations at the crRNA target sites, suggesting that viruses targeted by CRISPR-Cas13 systems may not be able to evolve resistance to CRISPR-Cas13 machinery [55].

Conclusion

In conclusion, CRISPR-CasRx enables programmable RNA virus interference in plants and has the potential to be optimized and developed in the future for transcriptome engineering applications. CRISPR-CasRx expands the power of CRISPR systems for virus interference applications and may help to address basic plant virology questions.

Methods

Vector construction

Construction of Cas13 orthologues for in planta expression

Nine different Cas13 variants were constructed for activity screening in plants and were cloned in binary vectors for in planta expression. The DNA sequences of the LshCas13a variants were plant codon optimized and subsequently cloned into binary vectors, as described previously [33]. Similarly, plant codon-optimized versions of BzCas13b variants with fusions of the 3xHA tag at the N terminus were designed and synthesized (Blue Heron Biotech), either with or without the NLS-encoding sequence at the C terminus, and were flanked by attL1 and attL2 sites to facilitate Gateway cloning into binary vectors (Additional file 1: Sequence S1, S2). Subsequently, BzCas13b-NLS and BzCas13b (no NLS) were subcloned into pK2GW7 binary vectors using LR Gateway recombination cloning to generate pK2GW7-BzCas13b, in which expression is controlled by the constitutive 35S cauliflower mosaic virus promoter (Additional file 1: Map S1). For the other Cas13 orthologues, including LwaCas13a-NLS/NES, PspCas13b-NES, CasRx-NLS, and dCasRx-NLS, the DNA sequences were amplified from plasmids provided by Addgene (Addgene plasmid numbers 91902/105815, 103862, 109049, and 109050, respectively) using primers that add a CACC sequence to the 5′ end of the amplified sequence to facilitate subsequent directional cloning into pENTR/D-TOPO vectors (Invitrogen) (Additional file 1: Table S2). Subsequently, Cas13 sequences were transferred from the pENTR/D-TOPO entry vectors into pK2GW7 binary vectors using Gateway recombination reactions. To generate the CasRx-NES version, the CasRx ORF sequence was amplified from the CasRx-NLS vector (plasmid number 109049) using primers that remove NLS sequences from the N and C termini of CasRx-NLS ORF, leading to the removal of the HA tag and GFP-encoding sequences, and BsaI sites were added at both ends of the isolated ORF (Additional file 1: Table S2). To add nuclear export sequences (NES), DNA sequences encoding NES and NES-1xHA tags were purchased as single-stranded DNA (ssDNA) oligonucleotides from Sigma with overhangs complementary to the overhangs generated from digestion of the BsaI sites at the two ends of the isolated CasRx ORF sequence (Additional file 1: Table S2). Phosphorylated and annealed double-stranded DNA (dsDNA) oligonucleotides were ligated into the BsaI-digested DNA sequence of the CasRx ORF to produce the 5′-NES-CasRx-NES-HA-3′ sequence. Subsequently, the ligation product of 5′-NES-CasRx-NES-HA-3′ was PCR amplified and subcloned into the pENTR/D-TOPO vector, followed by transfer into the pK2GW7 binary vector using LR Gateway recombination cloning. All the Cas13 variant sequences were confirmed in pK2GW7 binary vectors with Sanger sequencing.

Generation of the TRBO-BFP construct

To generate the TRBO infectious clone expressing BFP, the EBFP-encoding sequence was amplified from the LeGO-EBFP2 vector (Addgene plasmid number 85213) with primers that add AseI and NotI restriction sites to the 5′ and 3′ ends of the amplified BFP sequence, respectively (Additional file 1: Table S2). The GFP-encoding sequence was removed from the pJL-TRBO-G clone (Addgene plasmid number 80083) by digestion with AseI and NotI and replaced by the ligation of the AseI- and NotI-digested PCR product of BFP into the digested vector.

Construction of crRNA expression constructs

To construct the crRNA expression clones, the orthologue’s corresponding direct repeat with the targeting or non-targeting (NS) sequences was purchased from Sigma as ssDNA oligonucleotides (all crRNA sequences used in this study are listed in Additional file 1: Table S1). The phosphorylated and annealed dsDNA fragments of crRNAs were ligated into the TRV RNA2 genome under the pPEBV promoter using XbaI and either XhoI or BamHI restriction sites.

Plant material

Two- to 3-week-old wild-type Nicotiana benthamiana plants grown under long-day conditions (16 h light, 8 h dark at 25 °C) were used for all transient experiments. To generate transgenic Nicotiana benthamiana plants overexpressing different Cas13 orthologues, Agrobacterium tumefaciens GV310 strains harboring the selected constructs were transformed into N. benthamiana leaf discs following previous published protocols [56]. Generated plants (T1) were screened for the expression of the Cas13 protein by immunoblotting using an anti-HA antibody, and seeds harvested from the Cas13-positive plants were screened on kanamycin-containing MS (Murashige and Skoog) media. Surviving seedlings were used for virus-targeting experiments.

Agroinfiltration of N. benthamiana leaves

The different constructs (including all pK2GW7-Cas13 constructs, the TRV RNA1 genome, the engineered TRV RNA2 genomes harboring crRNA sequences under the pPEBV promoter, and the different infectious clones of the RNA viruses, including TRBO-GFP, TRBO-BFP, PVX-GFP, or TuMV-GFP) were individually electroporated into A. tumefaciens strain GV3101. Single colonies grown overnight in selective medium were centrifuged and suspended in infiltration medium (10 mM MES [pH 5.7], 10 mM MgCl2, and 200 μM acetosyringone) and incubated at ambient temperature for 2 h. For infiltration into leaves of N. benthamiana plants, cultures were mixed at a final OD600 of 0.5 for pK2GW7-Cas13, 0.1 for RNA1 and RNA2-crRNAs, 0.05 for TuMV-GFP and PVX-GFP, and 0.005 for TRBO-GFP and TRBO-BFP, relative to the suspended culture grown overnight. Healthy and fully developed leaves were selected for experiments, and agroinfiltration was done with a 1-ml needleless syringe.

Fluorescent protein intensity assays and quantitative fluorescence intensity analysis of images

For GFP and/or BFP expression in transient studies, infiltrated leaves were cut and then photographed after 3 days under a handheld UV light in the dark. For the systemic TuMV-expressed GFP in wild-type and transgenic plants, the whole plants were photographed after 6 to 7 days post-infiltration (dpi). For the quantitative measurements of the GFP signal intensity in images, the mean pixel values of the TRBO-GFP gene expression images were analyzed with ImageJ software to estimate GFP expression across samples. For each sample (each target of Cas13 variants), three biological replicates, represented as three different infiltrated leaves from three different plants, were analyzed. For each leaf, the GFP signal intensities were obtained as the corrected total cell fluorescence (CTCF) following this formula: CTCF = integrated density − (area of selected cell × mean fluorescence of background readings). The CTCF values of all three biological replicates for each sample were then averaged, and the normalized average of each target was compared to the normalized average CTCF value of the NS target. All values were limited to max = 100.

Immunoblot analysis

Total proteins were extracted from 100 mg of sample using extraction buffer (100 mM Tris-Cl, pH 8, 150 mM NaCl, 0.6% IGEPAL, 1 mM EDTA, and 3 mM DTT) with protease inhibitors (PMSF, leupeptin, aprotinin, pepstatin, antipain, chymostatin, Na2VO3, NaF, MG132, and MG115). Proteins were separated on a 10% polyacrylamide gel. Immunoblot analysis was conducted using mouse α-GFP antibody (1:3000, Invitrogen) for virus-expressed GFP and rat α-HA (1:1000) antibody for detection of Cas13 proteins. The antigens were detected by chemiluminescence using an ECL-detecting reagent (Thermo Scientific). The quantitative graphs summarizing the average of three independent biological replicates were produced by calculating the densiometric data generated by the relative quantification of protein bands from immunoblot membranes using ImageJ software. The average data of three independent immunoblot replicates were calculated as previously described [57]. Results were presented as fold change relative to the NS and/or virus-only control samples.

RNA extraction and RT-qPCR for analysis of viral RNA genomes

Total RNA was extracted from infiltrated leaves (in transient assays) or systemic leaves (in Cas13 overexpression lines) using Direct-zol RNA Miniprep Kits (Zymo Research) following the manufacturer’s instructions. Viral RNA was quantified with one-step RT-qPCR using the iTaq Universal SYBR Green One-Step Kit (Bio-Rad). All RT-qPCR reactions were performed in 10 μl reactions with three technical replicates in either 384-well or 96-well format and read out using a StepOnePlus™ Real-Time PCR System (Applied Biosystem) for the 96-well plates or the CFX384 Touch™ Real-Time PCR detection system (Bio Rad) for the 384-well plates. Expression levels were calculated by subtracting the housekeeping reference gene (tobacco PP2A [58]) cycle threshold (Ct) values from target Ct values to normalize for total input, resulting in ΔCt levels. Relative transcript abundance was computed as 2−ΔΔCt [59]. All samples were performed as three biological replicates.

Availability of data and materials

Key clones used in this study are available from the Laboratory of Genome Engineering and Synthetic Biology at KAUST.

References

  1. Mahy, M.v.R.B., Desk Encyclopedia of Plant and Fungal Virology. Academic Press, virology 632. San Diego: Elsevier; 2009.

  2. Boualem A, Dogimont C, Bendahmane A. The battle for survival between viruses and their host plants. Curr Opin Virol. 2016;17:32–8.

    Article  CAS  PubMed  Google Scholar 

  3. Anderson PK, et al. Emerging infectious diseases of plants: pathogen pollution, climate change and agrotechnology drivers. Trends Ecol Evol. 2004;19(10):535–44.

    Article  PubMed  Google Scholar 

  4. King AMQ, et al. Changes to taxonomy and the International Code of Virus Classification and Nomenclature ratified by the International Committee on Taxonomy of Viruses (2018). Arch Virol. 2018;163(9):2601–31.

    Article  CAS  PubMed  Google Scholar 

  5. Zaidi SS, et al. Engineering plant immunity: using CRISPR/Cas9 to generate virus resistance. Front Plant Sci. 2016;7:1673.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Kaufmann SHE, et al. Host-directed therapies for bacterial and viral infections. Nat Rev Drug Discov. 2018;17(1):35–56.

    Article  CAS  PubMed  Google Scholar 

  7. Agius C, et al. RNA silencing and antiviral defense in plants. Methods Mol Biol. 2012;894:17–38.

    Article  CAS  PubMed  Google Scholar 

  8. Obbard DJ, et al. The evolution of RNAi as a defence against viruses and transposable elements. Philos Trans R Soc Lond Ser B Biol Sci. 2009;364(1513):99–115.

    Article  CAS  Google Scholar 

  9. Barrangou R, Marraffini LA. CRISPR-Cas systems: prokaryotes upgrade to adaptive immunity. Mol Cell. 2014;54(2):234–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Makarova KS, et al. Evolution and classification of the CRISPR-Cas systems. Nat Rev Microbiol. 2011;9(6):467–77.

    Article  CAS  PubMed  Google Scholar 

  11. Shmakov S, et al. Discovery and functional characterization of diverse class 2 CRISPR-Cas systems. Mol Cell. 2015;60(3):385–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Adli M. The CRISPR tool kit for genome editing and beyond. Nat Commun. 2018;9(1):1911.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Cong L, et al. Multiplex genome engineering using CRISPR/Cas systems. Science. 2013;339:819–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Wang H, La Russa M, Qi LS. CRISPR/Cas9 in genome editing and beyond. Annu Rev Biochem. 2016;85:227–64.

    Article  CAS  PubMed  Google Scholar 

  15. Shmakov S, et al. Diversity and evolution of class 2 CRISPR-Cas systems. Nat Rev Microbiol. 2017;15(3):169–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Hsu PD, Lander ES, Zhang F. Development and applications of CRISPR-Cas9 for genome engineering. Cell. 2014;157:1262–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Mahas A, Mahfouz M. Engineering virus resistance via CRISPR-Cas systems. Curr Opin Virol. 2018;32:1–8.

    Article  CAS  PubMed  Google Scholar 

  18. Soppe JA, Lebbink RJ. Antiviral goes viral: harnessing CRISPR/Cas9 to combat viruses in humans. Trends Microbiol. 2017;25(10):833–50.

    Article  CAS  PubMed  Google Scholar 

  19. Price AA, Grakoui A, Weiss DS. Harnessing the prokaryotic adaptive immune system as a eukaryotic antiviral defense. Trends Microbiol. 2016;24(4):294–306.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Ali Z, et al. CRISPR/Cas9-mediated viral interference in plants. Genome Biol. 2015;16:238.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Ali Z, et al. CRISPR/Cas9-mediated immunity to Geminiviruses: differential interference and evasion. Sci Rep. 2016;6:26912.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Ali Z, et al. Pea early-browning virus-mediated genome editing via the CRISPR/Cas9 system in Nicotiana benthamiana and Arabidopsis. Virus Res. 2017.

  23. Tashkandi M, et al. Engineering resistance against tomato yellow leaf curl virus via the CRISPR/Cas9 system in tomato. Plant Signal Behav. 2018;13(10):e1525996.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Abudayyeh OO, et al. C2c2 is a single-component programmable RNA-guided RNA-targeting CRISPR effector. Science. 2016;353(6299):aaf5573.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Meeske AJ, Nakandakari-Higa S, Marraffini LA. Cas13-induced cellular dormancy prevents the rise of CRISPR-resistant bacteriophage. Nature. 2019;570:241–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Liu L, et al. Two distant catalytic sites are responsible for C2c2 RNase activities. Cell. 2017;168(1–2):121–34 e12.

    Article  CAS  PubMed  Google Scholar 

  27. East-Seletsky A, et al. Two distinct RNase activities of CRISPR-C2c2 enable guide-RNA processing and RNA detection. Nature. 2016;538(7624):270–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Gootenberg JS, et al. Nucleic acid detection with CRISPR-Cas13a/C2c2. Science. 2017;356(6336):438–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Gootenberg JS, et al. Multiplexed and portable nucleic acid detection platform with Cas13, Cas12a, and Csm6. Science. 2018;360(6387):439–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Sashital DG. Pathogen detection in the CRISPR-Cas era. Genome Med. 2018;10(1):32.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Abudayyeh OO, et al. Nucleic acid detection of plant genes using CRISPR-Cas13. CRISPR J. 2019;2:165–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Aman R, et al. Engineering RNA virus interference via the CRISPR/Cas13 machinery in Arabidopsis. Viruses. 2018;10(12).

    Article  CAS  PubMed Central  Google Scholar 

  33. Aman R, et al. RNA virus interference via CRISPR/Cas13a system in plants. Genome Biol. 2018;19(1):1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Zhan X, et al. Generation of virus-resistant potato plants by RNA genome targeting. Plant Biotechnol J. 2019;17:1814–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Zhang T, et al. Establishing CRISPR/Cas13a immune system conferring RNA virus resistance in both dicot and monocot plants. Plant Biotechnol J. 2019;17:1185–87.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Yan WX, et al. Cas13d is a compact RNA-targeting type VI CRISPR effector positively modulated by a WYL-domain-containing accessory protein. Mol Cell. 2018;70(2):327–39 e5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Cox DBT, et al. RNA editing with CRISPR-Cas13. Science. 2017;358:1019–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Abudayyeh OO, et al. RNA targeting with CRISPR-Cas13. Nature. 2017;550(7675):280–4.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Cox DBT, et al. RNA editing with CRISPR-Cas13: Science; 2017;358:1019–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Konermann S, et al. Transcriptome engineering with RNA-targeting type VI-D CRISPR effectors. Cell. 2018;173(3):665–76 e14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Makarova KS, et al. An updated evolutionary classification of CRISPR-Cas systems. Nat Rev Microbiol. 2015;13(11):722–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Smargon AA, et al. Cas13b is a type VI-B CRISPR-associated RNA-guided RNase differentially regulated by accessory proteins Csx27 and Csx28. Mol Cell. 2017;65(4):618–30 e7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Ali Z, et al. Efficient virus-mediated genome editing in plants using the CRISPR/Cas9 system. Mol Plant. 2015;8(8):1288–91.

    Article  CAS  PubMed  Google Scholar 

  44. Lindbo JA. TRBO: a high-efficiency tobacco mosaic virus RNA-based overexpression vector. Plant Physiol. 2007;145(4):1232–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Price AA, et al. Cas9-mediated targeting of viral RNA in eukaryotic cells. Proc Natl Acad Sci U S A. 2015;112(19):6164–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Giritch A, et al. Rapid high-yield expression of full-size IgG antibodies in plants coinfected with noncompeting viral vectors. Proc Natl Acad Sci U S A. 2006;103(40):14701–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Komor AC, Badran AH, Liu DR. CRISPR-based technologies for the manipulation of eukaryotic genomes. Cell. 2017;168(1–2):20–36.

    Article  CAS  PubMed  Google Scholar 

  48. O'Connell MR. Molecular mechanisms of RNA targeting by Cas13-containing type VI CRISPR-Cas systems. J Mol Biol. 2019;431(1):66–87.

    Article  CAS  PubMed  Google Scholar 

  49. Zhang T, et al. Establishing RNA virus resistance in plants by harnessing CRISPR immune system. Plant Biotechnol J. 2018;16:1415–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Mahas A, Neal Stewart C Jr, Mahfouz MM. Harnessing CRISPR/Cas systems for programmable transcriptional and post-transcriptional regulation. Biotechnol Adv. 2017;36:295–310.

    Article  CAS  PubMed  Google Scholar 

  51. Ali Z, Mahas A, Mahfouz M. CRISPR/Cas13 as a tool for RNA interference. Trends Plant Sci. 2018;23:374–8.

    Article  CAS  PubMed  Google Scholar 

  52. Chertow DS. Next-generation diagnostics with CRISPR. Science. 2018;360(6387):381–2.

    Article  PubMed  Google Scholar 

  53. Maxmen A. Faster, better, cheaper: the rise of CRISPR in disease detection. Nature. 2019;566(7745):437.

    Article  CAS  PubMed  Google Scholar 

  54. Mehta D, et al. Linking CRISPR-Cas9 interference in cassava to the evolution of editing-resistant geminiviruses. Genome Biol. 2019;20(1):80.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Freije CA, et al. Programmable inhibition and detection of RNA viruses using Cas13. Mol Cell. 2019. https://doi.org/10.1016/j.molcel.2019.09.013.

  56. A simple and general method for transferring genes into plants. Science, 1985. 227(4691): p. 1229–31.

  57. Taylor SC, Posch A. The design of a quantitative western blot experiment. Biomed Res Int. 2014;2014:361590.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. Liu D, et al. Validation of reference genes for gene expression studies in virus-infected Nicotiana benthamiana using quantitative real-time PCR. PLoS One. 2012;7(9):e46451.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We wish to thank members of the Laboratory for Genome Engineering and Synthetic Biology at King Abdullah University of Science and Technology for helpful discussions and comments on the manuscript.

Peer review information: Yixin Yao was the primary editor(s) on this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.

Review history

The review history is available as Additional file 2.

Funding

This publication is based upon work supported by the King Abdullah University of Science and Technology (KAUST) Office of Sponsored Research (OSR) under Award No. OSR-2015-CRG4-2647.

Author information

Authors and Affiliations

Authors

Contributions

MM conceived the research. MM, RA, and AM designed the research. AM and RA performed the research. MM, AM, and RA analyzed the data. AM, RA, and MM wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Magdy Mahfouz.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

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: Figure S1.

TRBO-GFP virus genome. Figure S2. Cas13 proteins are required for RNA virus interference. Figure S3. CasRx mediates efficient interference against TuMV-GFP virus by preventing its systemic spread in wild type N. benthamiana plants. Figure S4. Confirmation of Cas13 protein expression in permanent N. benthamiana lines. Figure S5. Multiplexed targeting of TRBO-BFP and PVX-GFP viruses through delivery of two different crRNAs with CasRx. Sequence S1. BzCas13b amino acid sequence (3xHA-BzCas13b-NLS). Sequence S2. BzCas13b full-length plant codon optimized DNA sequence (3x-HA-BzCas13b-NLS). Map S1. BzCas13b in pK2GW7 (HA-BzCas13b-NLS). Table S1. crRNA sequences used in this study. Table S2. primers used in this study.

Additional file 2.

Review history.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mahas, A., Aman, R. & Mahfouz, M. CRISPR-Cas13d mediates robust RNA virus interference in plants. Genome Biol 20, 263 (2019). https://doi.org/10.1186/s13059-019-1881-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13059-019-1881-2

Keywords