Next Article in Journal / Special Issue
Intranasal Boosting with Spike Fc-RBD of Wild-Type SARS-CoV-2 Induces Neutralizing Antibodies against Omicron Subvariants and Reduces Viral Load in the Nasal Turbinate of Mice
Previous Article in Journal
Rooting and Dating Large SARS-CoV-2 Trees by Modeling Evolutionary Rate as a Function of Time
Previous Article in Special Issue
Nonthermal Biocompatible Plasma Inactivation of Coronavirus SARS-CoV-2: Prospects for Future Antiviral Applications
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 15
Issue 3
10.3390/v15030686
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Efficient CRISPR-Cas13d-Based Antiviral Strategy to Combat SARS-CoV-2

by
Mouraya Hussein
1,
Zaria Andrade dos Ramos
1,
Monique A. Vink
1,
Pascal Kroon
1,
Zhenghao Yu
1,
Luis Enjuanes
2,
Sonia Zuñiga
2,
Ben Berkhout
1 and
Elena Herrera-Carrillo
1,*
1
Laboratory of Experimental Virology, Department of Medical Microbiology, Amsterdam UMC, Academic Medical Center, University of Amsterdam, 1105 AZ Amsterdam, The Netherlands
2
Department of Molecular and Cell Biology, National Center of Biotechnology (CNB-CSIC), Campus Universidad Autónoma de Madrid, 28049 Madrid, Spain
*
Author to whom correspondence should be addressed.
Viruses 2023, 15(3), 686; https://doi.org/10.3390/v15030686
Submission received: 17 January 2023 / Revised: 27 February 2023 / Accepted: 28 February 2023 / Published: 6 March 2023
(This article belongs to the Special Issue Advances in Antiviral Agents against SARS-CoV-2 and Its Variants)
Browse Figures
Versions Notes

Abstract

:
The current SARS-CoV-2 pandemic forms a major global health burden. Although protective vaccines are available, concerns remain as new virus variants continue to appear. CRISPR-based gene-editing approaches offer an attractive therapeutic strategy as the CRISPR-RNA (crRNA) can be adjusted rapidly to accommodate a new viral genome sequence. This study aimed at using the RNA-targeting CRISPR-Cas13d system to attack highly conserved sequences in the viral RNA genome, thereby preparing for future zoonotic outbreaks of other coronaviruses. We designed 29 crRNAs targeting highly conserved sequences along the complete SARS-CoV-2 genome. Several crRNAs demonstrated efficient silencing of a reporter with the matching viral target sequence and efficient inhibition of a SARS-CoV-2 replicon. The crRNAs that suppress SARS-CoV-2 were also able to suppress SARS-CoV, thus demonstrating the breadth of this antiviral strategy. Strikingly, we observed that only crRNAs directed against the plus-genomic RNA demonstrated antiviral activity in the replicon assay, in contrast to those that bind the minus-genomic RNA, the replication intermediate. These results point to a major difference in the vulnerability and biology of the +RNA versus −RNA strands of the SARS-CoV-2 genome and provide important insights for the design of RNA-targeting antivirals.

1. Introduction

The outbreak of Coronavirus Disease 2019 (COVID-19) at the end of 2019 has created a profound global health burden with pandemic-scale numbers of infected and fatal cases [1,2,3]. COVID-19 is caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) [4,5]. As of August 2022, the World Health Organization has confirmed some 572 million human cases worldwide and 6.4 million (1.1%) deaths. The pandemic spurred research towards the development of protective vaccines at an unprecedented pace, resulting in the approval of four different vaccines by the U.S. Food and Drug Administration (FDA) and the European Medicines Agency (EMA) [6,7,8,9,10,11]. However, new SARS-CoV-2 variants continue to appear and some are designated variants of concern (VOC) that may escape from these vaccines, e.g., the recent Omicron lineage [12,13]. This means that it remains prudent to develop alternative antiviral strategies in parallel. This therapeutic track should ideally consider the genetic variation among SARS-CoV-2 strains and thus prepare for a future zoonotic outbreak of an unrelated or yet unknown coronavirus [14,15,16]. The FDA has approved different antiviral drugs for COVID-19 treatment (Remdesivir and Baricitinib) and authorized others for emergency use (Molnupiravir and Paxlovid). However, there is currently no antiviral treatment that directly attacks the viral RNA genome [17,18]. To this end, the innate RNA interference (RNAi) mechanism was used to develop a novel therapeutic option against the virus [19,20,21]. We and others proposed an alternative antiviral approach based on the CRISPR genome editing tool for the specific recognition and degradation of the SARS-CoV-2 RNA genome [22,23,24].
The current COVID-19 pandemic, which started in 2019, represents the third zoonotic outbreak of a pathogenic coronavirus (SARS-CoV-2) in a relatively short timeframe after MERS-CoV in 2012 [25] and SARS-CoV in 2003 [26,27]. Moreover, the continuous emergence of new virus variants in the current SARS-CoV-2 pandemic is urging us to design antiviral strategies that are not only robust and specific, but ideally also capable of covering the genetic variation present in the virus variants. The clustered regularly interspaced short palindromic repeat (CRISPR)-Cas toolbox provides a relatively straightforward platform to design antivirals that target viral RNA or DNA in a sequence-specific manner [28,29,30]. We opted to target the SARS-CoV-2 RNA genome with the CRISPR-Cas13 system because it selectively cleaves RNAs that are complementary to the designed CRISPR-associated RNA (crRNA) and because it theoretically allows for an RNA attack anywhere in the cell (nucleus, cytosol and other subcellular compartments such as double-membrane vesicles (DMVs) that support viral RNA synthesis).
Four different RNA-targeting Cas13 orthologues have been described (Cas13a, b, c and d) [31]. We selected Cas13d as it combines superior activity with a small size, which is beneficial for gene delivery by means of viral vectors [32,33,34]. Moreover, Cas13d lacks sequence constraints besides the crRNA target sequence, meaning that in principle any viral RNA sequence can be attacked. CRISPR-Cas13d uses a customizable ~23-nucleotide (nt) spacer to direct the endonuclease to specific RNA molecules for targeted RNA degradation [31,35,36,37]. We propose to attack either the SARS-CoV-2 +RNA or –RNA strand. Once a host cell is infected by SARS-CoV-2, the plus-genomic RNA (+gRNA) genome is released into the cytoplasm to act as mRNA for translation of the polyproteins pp1a and pp1ab (Figure 1A). The polyproteins are proteolytically processed into the different non-structural proteins (nsps) with different functions. Several viral proteins constitute the replication/transcription complex (RTC) that associates with DMVs in the perinuclear regions. That is where a full-length complementary minus-genomic RNA (−gRNA) is synthetized, which is subsequently used as template for the production of new +gRNA molecules. Additionally, a set of 3′ co-terminal sub-genomic mRNAs (sgmRNAs) with an identical 5′ leader sequence are produced from intermediate minus sub-genomic RNAs (sgRNAs) via an orchestrated process of discontinuous transcription (Figure 1A). Thus, the +gRNA and all sgmRNAs share an identical leader sequence in the 5′UTR (untranslated region) and an identical 3′UTR, making these terminal domains ideal targets for an antiviral attack on all sgmRNAs. In theory, both the viral + and −RNA strands can be attacked, but the latter may be an easier target because of the much lower copy number in the infected cell.
We set out to find optimal target sequences in the SARS-CoV-2 RNA genome for Cas13d (Figure 1B). We selected the most conserved targets [33,34] that cover all VOCs known to date and SARS-CoV [26,38]. An additional rationale for targeting conserved sequences is that viral escape will be more difficult at sites that encode important viral sequences, and escape by mutation will likely cause a drop in replication capacity [39]. For instance, the highly conserved ribosomal frameshift site has been an attractive target for anti-SARS-CoV-2 strategies, which resulted in the inhibition of frameshifting during RNA translation and virus replication [40,41,42]. Based on the alignment of the complete viral genomes of SARS-CoV and the current SARS-CoV-2 strains, we composed a CRISPR platform comprised of two sets of 29 crRNAs targeting highly conserved regions on the + or −RNA strand. We aimed at obtaining crRNAs with a broad antiviral potential against multiple human coronaviruses and at gaining mechanistic insights into factors that could affect the target RNA vulnerability based on the specific biology of human coronaviruses.

2. Materials and Methods

2.1. Design of crRNAs for a Cas13d-Mediated Attack on SARS-CoV-2

Cas13d employs a customizable crRNA that directs the RNA endonuclease to a specific sequence on the target RNA for cleavage and subsequent degradation [36]. We first ran a computational model (R script) developed by Wessels et al. along the SARS-CoV-2 genome of the Wuhan Hu-1 reference strain (MN908947) to identify 23 nt crRNAs with a high predicted efficacy [43] (Figure 2A). Then, we selected conserved target sequences by aligning the reference SARS-CoV-2 genome with 31,576 viral genome sequences, including the VOCs Alpha, Beta, Gamma, Delta and Omicron as well as the de-escalated variants Lambda and Mu (https://www.gisaid.org/ (accessed on 8 December 2021, Table S1) [22]. In addition, SARS-CoV-2 sequences were aligned with the complete genomes of SARS-CoV (260 strains). The Shannon entropy of the aligned sequences was calculated as a measure of the genetic variability per nucleotide position and plotted with Bioedit in Figure 2A [44,45]. The Shannon entropy ranges from 0 to 1, where values close to zero represent low genetic diversity (Table S2). Based on the generated list of crRNAs from the R script and the conserved target regions selected from the multiple whole-genome alignments, a set of 29 crRNAs was selected with no predicted off-target effect based on alignment with the complete transcriptome of human cells (Table 1). The selected target sites on the SARS-CoV-2 RNA genome are indicated in Figure 2A. We selected predicted weak and strong crRNAs to validate the power of the on-target prediction tool developed by Wessels et al. (Table 1).
Two sets of crRNAs were designed against either the + or the -RNA strand of SARS-CoV-2. We selected targets in the 5′UTR and 3′UTR as the genomic and sub-genomic RNAs share common sequences in these termini. We included crRNAs targeting the transcription-regulating sequence (TRS) present the 5′-leader sequence (TRS-Leader) and at the gene junctions (TRS-B), which we named 5′UTR-1 and 5′UTR-2, as these are critical for the key process of discontinuous transcription and present in the +gRNA as well as sgmRNAs [43,44]. Another highly conserved element at the 3′ end of the coronavirus genome is the s2m motif, which was recently proposed as a promising target for antivirals [45,46]. We also selected target sequences that encode critical enzymes involved in virus replication and transcription, e.g., the helicase and RNA-dependent RNA polymerase (RdRp) (Table 1).

2.2. Plasmid Construction

The expression plasmid pLentiRNACRISPR_005-hU6-DR_BsmBI-EFS-RfxCas13d-NLS-2A-Puro-WPRE (Addgene, #138147) expresses the RfxCas13d endonuclease in mammalian cells. This construct was kindly donated by Neville Sanjana [47]. The EF-1a core promoter was used for expression of the CaS13d nuclease and the U6 polymerase-3 promoter was used for crRNA expression. Oligonucleotides encoding SARS-CoV-2-targeting crRNAs and non-targeting control crRNA (NC) were ligated into the Esp3i site of the pLentiRNACRISPR_005-hU6-DR_BsmBI-EFS-RfxCas13d-NLS-2A-Puro-WPRE vector. All crRNAs are listed in Table S3, including a non-targeting control crRNA that does not target any sequence in the human and coronavirus genomes. All constructs were sequence-verified using the BigDye Terminator Cycle Sequencing kit (ABI). For sequencing of crRNA constructs a denaturation temperature of 98 °C was used in the presence of 1M Betaine to disrupt structures in the DNA template. The luciferase reporter plasmid was constructed by the insertion of a 250 bp fragment into the Xbal site of the pGL3-control plasmid (Promega; GenBank: U47296.2) to create EcoRI and PstI sites (Figure 2B) [48]. DNA fragments encoding the luciferase reporter with multiple CRISPR targets were synthesized by Integrated DNA Technologies (IDT, Coralville, Iowa) and cloned, one for +gRNA as well as sgmRNAs targets and one for −gRNA as well as sgRNAs targets. In the context of our experiments, we will refer to +gRNA and sgmRNAs as +RNA targets and –gRNA and sgRNAs as –RNA targets.
The SARS-CoV-2 infectious cDNA clone based on the SARS-CoV-2 Wuhan Hu-1 strain (MN908947) was assembled in a bacterial artificial chromosome (BAC) using an approach described previously [49]. To generate a SARS-CoV-2 replicon that can be used in a BSL-2 level laboratory, the structural and genus-specific genes were removed and replaced by a reporter gene (Figure 2C). The nucleocapsid (N) gene was maintained as it is required for efficient RNA synthesis [50]. A 4028 bp fragment flanked by SanDI and RsrII restriction sites was chemically synthesized (IDT, Coralville, Iowa). This DNA fragment encodes part of the SARS-CoV-2 genome (nt 20085 to 21555 corresponding to ORF1b), the mNeonGreen (mNG) reporter gene under control of the transcription regulating sequence (TRS) of the membrane (M) gene (TRS-M, nt 26431 to 26522), and another SARS-CoV-2 fragment (nt 28177 to 29533 encoding the N protein and its TRS). The synthetic fragment was inserted into the SanDI and RsrII restriction sites of pBAC-SARS-CoV2 to generate the replicon pBAC-SARSCoV-2-mNG (Figure 2C). The SARS-CoV replicon pBAC-SARS-CoV-REP was previously generated for the SARS-CoV Urbani strain (AY278741) [51]. pBAC-SARS-CoV-REP encodes a modified version of Montastraea cavernosa Green Fluorescent Protein (hmGFP) (Figure 2C).

2.3. Cell Culture and DNA Transfection

Human embryonic kidney (HEK) 293T cells were cultured in DMEM (Life Technologies, Invitrogen, Carlsbad, CA, USA) supplemented with 10% fetal calf serum (FCS), penicillin (100 U/mL), streptomycin (100 mg/mL) and 1% L-glutamine. The cells were cultured in a humidified chamber at 37 °C and 5% CO2.
For luciferase assays, HEK293T cells were seeded one day before transfection in 24-well plates at a density of 1.4 × 105 cells per well in 0.5 mL media. Cells were transfected with 100 ng of the firefly luciferase expression plasmid, 1 ng of Renilla luciferase expression plasmid (pRL) and 300 ng of CRISPR-Cas13d/crRNA vector using Lipofectamine 2000 reagent (Invitrogen) according to the manufacturer’s instructions [52]. For the titration transfection, cells were transfected with 100 ng of firefly luciferase expression plasmid, 1 ng of pRL and 75, 150 or 300 ng of CRISPR-Cas13d/crRNA. pBS plasmid was supplemented to reach an equal amount of DNA (300 ng) per sample. The knockdown efficiency was quantified by measuring the firefly luciferase (Luc) fluorescence signal at two days post-transfection according to the manufacturer’s instructions [53]. A non-targeting crRNA served as negative control (NC), and the luciferase level scored for this construct was set at 100%. We performed three independent transfections, each in duplicate. The ratio between firefly and Renilla luciferase activity was used for normalization of the experimental variation such as differences in the transfection efficiency. The luciferase data were subsequently corrected for between-session variation using Factor Correction v10.5 [54]. The resulting six values were used to calculate the standard deviation shown as error bars. Data were analyzed with the Prism software (GraphPad Prism 9.1.0). One-way ANOVA was used for all statistical analyses: * p ≤ 0.05; ** p ≤ 0.01; *** p ≤ 0.001; **** p ≤ 0.0001.

2.4. SARS-CoV-2 Replicon Assay

HEK293T cells were plated a day before transfection in 24-well plates at a density of 1.4 × 105 cells per well in 0.5 mL media. For transfection, we used 300 ng vector encoding Cas13d and crRNA, 1.43 μg pBAC-SARS-CoV-REP or pBAC-SARS-CoV-2-mNG, and Lipofectamine 2000 according to the manufacturer’s instructions (Invitrogen) [52]. At two days post-transfection, cells were washed twice with FACS (Fluorescence Activated Cell Sorting) buffer (PBS supplemented with 5% FCS and 2mM EDTA). The number of GFP- or mNeonGreen-positive cells was measured with flow cytometry and analyzed using FlowJo™ v10.7 [50,55]. Three independent experiments were performed in duplicate. The resulting six values were factor-corrected for between-session variation and used to calculate the standard deviation shown as error bars. Data were analyzed with the Prism software (GraphPad Prism 9.1.0). One-way ANOVA was used for all statistical analyses: * p ≤ 0.05; ** p ≤ 0.01; *** p ≤ 0.001; **** p ≤ 0.0001.

3. Results

3.1. Design of crRNAs against Highly Conserved SARS-CoV-2 RNA Sequences

We ran the R script algorithm developed by Wessels et al. [47] along the MN908947 reference genome of SARS-CoV-2 to identify 23 nt crRNAs with a high predicted efficacy. This yielded an initial collection of 28,749 candidate targets for Cas13d. Next, Bioedit was used to identify the crRNA targets that are most conserved among different SARS-CoV-2 variants and SARS-CoV. The R script places the candidate crRNAs in four quartiles based on the predicted efficacy (Q1–Q4), with Q1 for crRNAs of modest activity and Q4 for the most active crRNAs [47]. We included crRNA candidates of all four quartiles to validate the predictive value of the R script. The selected crRNAs target important open reading frames: six non-structural proteins (nsps), five RNA-dependent RNA polymerase (RdRp; nsp12), five N protein and four viral helicase (nsp13). In addition, important replication signals were also targeted: in the leader (two), 3′-pseudoknot (three), ribosomal frameshift signal (two) and 3′-UTR (two), including the highly conserved stem-loop II motif (s2m) (Table 1, Figure 2A). In total, we selected 29 viral targets, for which crRNAs were designed to target either the +RNA or -RNA strand of the SARS-CoV-2 genome, making a total of 58 crRNAs.

3.2. Targeting of SARS-CoV-2 RNA by Cas13d

To test the silencing efficiency of the designed crRNAs, HEK293T cells were co-transfected with a single Cas13d-crRNA plasmid (300 ng) and one of two luciferase reporter constructs that represent either the viral + or − strand (Figure 2B). A fixed amount of Renilla plasmid was included to control for variation in the transfection efficiency. A non-targeting crRNA that targets neither SARS-CoV-2 nor human sequences was used as negative control (NC), for which the luciferase activity was set at 100%. Most crRNAs against the viral +RNA strand induced a significant reduction in the Luc signal of at least 50% (Figure 3A), but not all were active, e.g., nsp4-1, nsp4-2, RdRp-3, helicase-2 and N-3. We obtained comparable results for the set of –RNA-strand-targeting crRNAs: most crRNAs demonstrated good knockdown activity, expect for nsp3-1, RdRp-2 and RdRp-3 (Figure 3B). Overall, we obtained efficient Cas13d-mediated inhibition with the majority of the designed crRNAs, for both the + and the − strand of the SARS-CoV-2 RNA genome.
With this data set, we validated the accuracy of the in silico crRNA prediction tool of Wessels for the 29 crRNA targets in the + RNA strand. We included eight crRNAs with poor predicted activity (nsp4-1, nsp4-2, nsp6-2, N-1, RdRp-3, RdRP-4, helicase-2 and helicase-4) and the remaining 21 crRNAs with optimal predicted activity (Table 1). The experimentally determined knockdown activity did not accurately correspond to the predicted crRNA activity. For instance, N-3 was predicted to have good silencing activity but showed poor activity and nsp6-2 is a good crRNA, but was predicted to be a weak one. The correlation between the predicted guide score generated by the Wessels algorithm and the measured crRNA activity for all 29 selected target sequences was calculated, showing a significant negative correlation coefficient of −0.48 (p = 0.0079). Therefore, we suggest that it remains of critical importance to experimentally validate the predicted crRNA candidates.

3.3. A Combinatorial Antiviral crRNA Approach

We selected the 14 best (+) strand-targeting crRNAs for a dose–response test in HEK293T cells co-transfected with a fixed amount of the Luc reporter (100 ng) and increasing amounts of the Cas13d-crRNA construct: 75, 150 and 300 ng (Figure 4). These results confirmed the inhibitory potency of these crRNAs, which act in a dose-dependent manner.
Like any other virus, SARS-CoV-2 will acquire spontaneous mutations in its genome during virus replication that could trigger viral escape. Viral escape will more likely occur when the crRNA inhibitor is applied as “mono-therapy”, but a combination of two or multiple crRNAs can result in additive or possibly even synergistic inhibition. A combinatorial attack will also increase the genetic threshold for the acquisition of crRNA resistance as more mutations will be required, which explains the clinical success of a combination drug therapy for HIV infection [56,57,58,59]. We selected the best crRNAs (5′UTR-1, nsp3-1, slippery, RdRp-2, helicase-3 and s2m) for a combinatorial approach and tested all possible combinations of two crRNAs in HEK293T cells co-transfected with the Luc reporter (100 ng), the Renilla plasmid (1 ng) and 300 ng of the CRISPR-Cas13d/crRNA constructs. Note that for a dual attack 150 ng of each crRNA construct was transfected, thus keeping a fixed amount of crRNA plasmid. The luciferase expression was measured at two days post-transfection and compared to the non-targeting control crRNA (NC), which was set at 100% (Figure 5). Overall, the dual crRNA combinations profoundly reduced the luciferase level to values significantly lower than those measured for the individual crRNAs. No additive effect was measured for 5′ UTR-1, possibly because it already showed strong inhibition as monotherapy.

3.4. The Impact of Cas13d on Intracellular SARS-CoV-2 Replication

We next tested the inhibitory capacity of these crRNAs on SARS-CoV-2 replication using the replicon system. Replicons are sub-genomic elements that replicate the viral RNA genome inside cells without the production of infectious virus and cell-to-cell viral spread. To generate a single-cycle SARS-CoV-2 replicon, an infectious cDNA clone was assembled in a bacterial artificial chromosome (BAC) as described [49], with the structural genes replaced by the mNeonGreen reporter gene. We selected the best + and − crRNA inhibitors and tested their inhibitory capacity in HEK293T cells co-transfected with the SARS-CoV-2 replicon (Figure 6A). Two days post-transfection the mNeonGreen signal was quantitated with flow cytometry, which reflects the level of SARS-CoV-2 RNA replication. The fluorescence measured for the control crRNA (NC) was set at 100%. All five crRNAs targeting the + strand SARS-CoV-2 RNA caused a robust reduction in the mNeonGreen signal (Figure 6A). Interestingly, such suppression was not observed for four out of five crRNAs that target the -RNA strand, with the exception of the crRNA against the slippery sequence that showed modest inhibitory activity. Thus, it seems that -RNA targeting crRNAs that were fully active in the Luc assay lose their knockdown activity in the replicon assay.
Next, we evaluated the breadth of these antiviral crRNAs against the replicon of the related SARS-CoV. In this case, we actually expected optimal breadth as we selected targets that are in fact identical between SARS-CoV-2 and SARS-CoV. We measured the fluorescence signal at two days post-transfection and observed the same trend as observed for SARS-CoV-2. The crRNAs that target the +RNA strand caused a robust reduction in the fluorescent signal, demonstrating the breadth of this antiviral approach (Figure 6B). No inhibition was apparent for crRNAs that target the -RNA strand, with the exception of crRNA-slippery. Thus, we could reproduce the apparent loss of activity for most -RNA targeting crRNAs against the SARS-CoV replicon. Importantly, crRNAs that suppress the SARS-CoV-2 replicon were also able to suppress SARS-CoV, thus demonstrating the breadth of this antiviral strategy.

3.5. The Impact of Cas13 on the Replication Intermediate minus Sense Transcripts

The lack of inhibitory activity of most crRNAs that target the -RNA strand in the replicon system is striking. To explain this, we considered the potential differences in the biology of the + and -RNA strands during coronavirus replication. A first difference is their abundance in infected cells. The -RNA strands are up to 100-fold less abundant than +RNA strands (38,69). Second, as a consequence of this difference in abundance, it is likely that most − strands will be annealed to the surplus of + strands. Such complexation is likely as both RNA strands accumulate in the same DMV compartment (69–74). DMVs are a hallmark of coronaviruses, which employ these membrane-associated structures for the replication of their RNA genome. It is possible that RNA duplex formation restricts the accessibility of all − strands for Cas13d attack, but not the surplus + strands that remain single-stranded (Figure 7A). To test this scenario, we designed a disrupted luciferase reporter with all SARS-CoV-2 + strand targets to produce a mimic of the viral + transcripts. We also made a luciferase reporter with all − targets. Both constructs were transfected in HEK293T cells with the plasmid (300 ng) expressing Cas13d and a crRNA that targets the -RNA strand. We mimicked the surplus of + strands by mixing an increasing amount of the + construct with a fixed amount of the − construct: 1:5, 1:25, 1:50, 1:75 and 1:100 (ng). We included a condition without the + construct (1:0) and compared the results to those obtained with the non-targeting crRNA (NC). We observed a gradual loss of activity of the -RNA targeting crRNA with titration of more +RNA construct (Figure 7B). Overall, this result supports the hypothesis that viral −RNA strands become inaccessible to the CRISPR-Cas13d endonuclease when forced into a stable RNA duplex by the expression of excess +RNA strands.

4. Discussion

The SARS-CoV-2 pandemic has triggered a search for antiviral drugs and other antiviral approaches. We show that the novel CRISPR-Cas13d endonuclease tool can be instructed for a sequence-specific attack on the viral RNA genome. We analyzed the SARS-CoV-2 genome to find the best crRNA targets of 23 nt based on our extensive expertise with sequence-specific antiviral strategies [60,61,62]. Viruses have a remarkably high mutational capacity, which creates genetic diversity that facilitates rapid adaptation to new conditions (e.g., antiviral drug pressure or a new host upon zoonotic transfer). We therefore considered viral targets of which the sequence is well-conserved among the current SARS-CoV-2 isolates. As we should prepare for a potential future outbreak of related viruses, we screened for targets that are conserved between SARS-CoV-2 and SARS-CoV. We designed a set of crRNAs to target the viral +RNA strand, but also the matching crRNAs that target the complementary -RNA strand. First, we experimentally validated the activity of these crRNAs in simple cell culture model systems with a luciferase reporter construct. The knockdown efficiency in transiently transfected cells was determined and the best crRNA antivirals were selected. This experimental validation has shown that the R script used to predict the on-target activity of crRNAs did not accurately predict the antiviral efficacy of the Cas13d-crRNAs. Therefore, experimental tests remain important in order to select the best crRNA candidates. We tested different combinations of the most potent crRNAs. The knockdown efficiency of dual crRNA combinations was significantly higher than those measured for the individual crRNAs.
Next, we used a non-infectious replicon system that contains all viral genes necessary for intracellular RNA replication and an exogenous reporter gene to quantitate the inhibition of viral RNA replication by the co-transfected crRNA constructs. CRISPR-Cas13d efficiently inhibited RNA replication when the crRNA targets the viral +RNA strand (Figure 6), but four out of five crRNAs that target the viral −RNA strand lost their inhibitory activity in the replicon assay. Consistent with these results, a recent study on antiviral siRNAs reported inactivity of those that target the viral −RNA strand [63]. We hypothesize that the differential vulnerability of viral + and −RNA strands can be explained by a profound difference in their abundance, which means that the less abundant −RNA transcripts will be annealed to the excess of viral +RNA strands inside the cells (Figure 7A). This scenario requires that both RNA transcripts be in the same location inside the infected cell, which is indeed the case as viral RNA replication occurs in DMVs in the cytoplasm. As a result of +/− duplex formation, Cas13d-mediated targeting of the viral −RNA strand is likely blocked. We mimicked this scenario with two luciferase reporters that express transcripts with a complementary sequence and demonstrated that crRNA inhibitory activity against a reporter is lost upon expression of the complementary RNA strand (Figure 7B). This result also reveals that viral RNA replication intermediates that are present inside DMVs are accessible to the Cas13d endonuclease.
Despite this dominant suppressive effect of RNA duplex formation on crRNA attack, we scored modest inhibition for the crRNA against the slippery sequence in the -RNA. This target site is unusual as it adopts a relatively stable RNA hairpin structure in the +gRNA [40,64,65] and possibly also a “mirror-image structure” in the -RNA. Such intramolecular base pairing may explain the special behavior of this target as it may interfere with intermolecular base pairing and formation of the +/−RNA duplex that facilitates crRNA annealing and Cas endonuclease cleavage. On the other hand, we observed good inhibition of several 5′/3′-UTR targets that are part of important RNA secondary structures (Figure 7A). Recent efforts to assess the influence of RNA secondary structure on Cas13 efficiency have shown a strong cleavage preference for single-stranded RNAs with 18 or more nucleotides available for base pairing to the crRNA [66]. We argue that local intramolecular base-pairing can interfere with the formation of the perfect intermolecular +/-RNA duplex and thus expose parts of the + or −RNA strand for Cas13d attack. Another factor that may influence RNA accessibility for Cas13d binding is the binding of proteins of viral or cellular origin to the viral RNA, thus forming a steric block for Cas13d annealing.
The breadth of the crRNAs was explored by testing the same crRNAs against both SARS-CoV-2 and SARS-CoV replicons. The crRNAs caused a robust reduction in the fluorescent signal, demonstrating the breadth of this antiviral approach. New SARS-CoV-2 variants will likely continue to emerge as this virus seems to be in the process of becoming endemic and there is the sustained, but low risk of yet another pandemic in the future of a “new” coronavirus. Thus, it remains of paramount importance to develop broadly active antivirals. We thus checked whether the final selection of 29 target viral sequences carry mutations in the new VOCs and confirmed that most crRNA targets are well conserved among all seven SARS-CoV-2 variants. Overall, the selected targets are conserved in at least 98% of SARS-CoV-2 variants, with the exception of 3’Pseudoknot-2 and s2m, which acquired mismatches in the Alpha and Delta variants, respectively. We have observed two insertions in the target region of 3’Pseudoknot-2 in the seed region (nucleotides 15–21) and a single point mutation for s2m. Thus, a reduced targeting efficiency of crRNA 3’Pseudoknot-2 against the Alpha variant is expected, which may not be the case for crRNA s2m with the Delta variant. The next step will therefore be to experimentally confirm the efficacy of our panel of crRNAs against these VOCs, including the highly mutated Omicron SARS-CoV-2 variant for which we found no mismatches with the designed crRNAs [67].
In conclusion, we designed and selected a panel of very active antiviral crRNAs that target highly conserved genome sequences of SARS-CoV-2 and confirmed that these antivirals also exhibit activity against the related SARS-CoV. Our crRNA design has proven to be successful as we maintain complementarity with the new SARS-CoV-2 VOCs. These results emphasize the importance of the rational design of gene-editing tools to develop antiviral strategies in order to prepare for future zoonotic outbreaks of novel pathogenic coronaviruses.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15030686/s1. Table S1: SARS-CoV-2 genome sequences available in the GISAID database; Table S2: Shannon entropy per nucleotide position; Table S3: crRNAs oligonucleotides.

Author Contributions

E.H.-C. designed the project; M.H., Z.A.d.R., M.A.V. and P.K. conducted the experiments; E.H.-C., B.B., M.H. and Z.Y. analyzed and studied the data, S.Z. and L.E. provided materials, E.H.-C., B.B. and M.H. wrote the original draft, B.B., S.Z. and L.E. reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Aspasia-NWO 015.015.040 grant.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data underlying this article are available in the article and in its online Supplementary Material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Spinelli, A.; Pellino, G. COVID-19 pandemic: Perspectives on an unfolding crisis. Br. J. Surg. 2020, 107, 785–787. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ahn, D.-G.; Shin, H.-J.; Kim, M.-H.; Lee, S.; Kim, H.-S.; Myoung, J.; Kim, B.-T.; Kim, S.-J. Current Status of Epidemiology, Diagnosis, Therapeutics, and Vaccines for Novel Coronavirus Disease 2019 (COVID-19). J. Microbiol. Biotechnol. 2020, 30, 313–324. [Google Scholar] [CrossRef] [PubMed]
  3. Chen, Y.; Liu, Q.; Guo, D. Emerging coronaviruses: Genome structure, replication, and pathogenesis. J. Med. Virol. 2020, 92, 418–423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Lai, C.C.; Shih, T.P.; Ko, W.C.; Tang, H.J.; Hsueh, P.R. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): The epidemic and the challenges. Int. J. Antimicrob. Agents 2020, 55, 105924. [Google Scholar] [CrossRef]
  5. Phan, T. Novel coronavirus: From discovery to clinical diagnostics. Infect. Genet. Evol. 2020, 79, 104211. [Google Scholar] [CrossRef]
  6. Torjesen, I. Covid-19: AstraZeneca vaccine is approved in EU with no upper age limit. BMJ 2021, 372, n295. [Google Scholar] [CrossRef]
  7. Mahase, E. Covid-19: Pfizer vaccine efficacy was 52% after first dose and 95% after second dose, paper shows. BMJ 2020, 371, m4826. [Google Scholar] [CrossRef]
  8. Meo, S.A.; Bukhari, I.A.; Akram, J.; Meo, A.S.; Klonoff, D.C. COVID-19 vaccines: Comparison of biological, pharmacological characteristics and adverse effects of Pfizer/BioNTech and Moderna Vaccines. Eur. Rev. Med. Pharmacol. Sci. 2021, 25, 1663–1669. [Google Scholar]
  9. Mahase, E. Covid-19: Moderna applies for US and EU approval as vaccine trial reports 94.1% efficacy. BMJ 2020, 371, m4709. [Google Scholar] [CrossRef]
  10. Burgos, R.M.; Badowski, M.E.; Drwiega, E.; Ghassemi, S.; Griffith, N.; Herald, F.; Johnson, M.; Smith, R.O.; Michienzi, S.M. The race to a COVID-19 vaccine: Opportunities and challenges in development and distribution. Drugs Context 2021, 10, 1–10. [Google Scholar] [CrossRef]
  11. Livingston, E.H.; Malani, P.N.; Creech, C.B. The Johnson & Johnson Vaccine for COVID-19. JAMA 2021, 325, 1575. [Google Scholar] [CrossRef]
  12. Choi, J.Y.; Smith, D. SARS-CoV-2 Variants of Concern. Yonsei Med. J. 2021, 62, 961–968. [Google Scholar] [CrossRef] [PubMed]
  13. Araf, Y.; Akter, F.; Tang, Y.D.; Fatemi, R.; Parvez MS, A.; Zheng, C.; Hossain, M.G. Omicron variant of SARS-CoV-2: Genomics, transmissibility, and responses to current COVID-19 vaccines. J. Med. Virol. 2022, 94, 1825–1832. [Google Scholar] [CrossRef] [PubMed]
  14. Wang, P.F.; Casner, R.G.; Nair, M.S.; Wang, M.; Yu, J.; Cerutti, G.; Liu, L.; Kwong, P.D.; Huang, Y.; Shapiro, L.; et al. Increased resistance of SARS-CoV-2 variant P.1 to antibody neutralization. Cell Host Microbe 2021, 29, 747–751.e4. [Google Scholar] [CrossRef]
  15. Abdool Karim, S.S.; de Oliveira, T. New SARS-CoV-2 Variants—Clinical, Public Health, and Vaccine Implications. N. Engl. J. Med. 2021, 384, 1866–1868. [Google Scholar] [CrossRef]
  16. Jogalekar, M.P.; Veerabathini, A.; Gangadaran, P. SARS-CoV-2 variants: A double-edged sword? Exp. Biol. Med. 2021, 246, 15353702211014146. [Google Scholar] [CrossRef]
  17. Ng, T.I.; Correia, I.; Seagal, J.; DeGoey, D.A.; Schrimpf, M.R.; Hardee, D.J.; Noey, E.L.; Kati, W.M. Antiviral Drug Discovery for the Treatment of COVID-19 Infections. Viruses 2022, 14, 961. [Google Scholar] [CrossRef] [PubMed]
  18. FDA. Coronavirus (COVID-19) | Drugs. 2022. Available online: https://www.fda.gov/drugs/emergency-preparedness-drugs/coronavirus-covid-19-drugs (accessed on 22 December 2022).
  19. Traube, F.R.; Stern, M.; Tolke, A.J.; Rudelius, M.; Mejias-Perez, E.; Raddaoui, N.; Kummerer, B.M.; Douat, C.; Streshnev, F.; Albanese, M.; et al. Suppression of SARS-CoV-2 Replication with Stabilized and Click-Chemistry Modified siRNAs. Angew. Chem. Int. Ed. Engl. 2022, 61, e202204556. [Google Scholar] [CrossRef] [PubMed]
  20. Baldassi, D.; Ambike, S.; Feuerherd, M.; Cheng, C.C.; Peeler, D.J.; Feldmann, D.P.; Porras-Gonzalez, D.L.; Wei, X.; Keller, L.A.; Kneidinger, N.; et al. Inhibition of SARS-CoV-2 replication in the lung with siRNA/VIPER polyplexes. J. Control. Release 2022, 345, 661–674. [Google Scholar] [CrossRef]
  21. Friedrich, M.; Pfeifer, G.; Binder, S.; Aigner, A.; Vollmer Barbosa, P.; Makert, G.R.; Fertey, J.; Ulbert, S.; Bodem, J.; Konig, E.M.; et al. Selection and Validation of siRNAs Preventing Uptake and Replication of SARS-CoV-2. Front. Bioeng. Biotechnol. 2022, 10, 801870. [Google Scholar] [CrossRef]
  22. Hussein, M.; Andrade dos Ramos, Z.; Berkhout, B.; Herrera-Carrillo, E. In Silico Prediction and Selection of Target Sequences in the SARS-CoV-2 RNA Genome for an Antiviral Attack. Viruses 2022, 14, 385. [Google Scholar] [CrossRef]
  23. Abbott, T.R.; Dhamdhere, G.; Liu, Y.; Lin, X.; Goudy, L.; Zeng, L.; Chemparathy, A.; Chmura, S.; Heaton, N.S.; Debs, R.; et al. Development of CRISPR as an Antiviral Strategy to Combat SARS-CoV-2 and Influenza. Cell 2020, 181, 865–876.e12. [Google Scholar] [CrossRef]
  24. Wang, L.; Zhou, J.; Wang, Q.; Wang, Y.; Kang, C. Rapid design and development of CRISPR-Cas13a targeting SARS-CoV-2 spike protein. Theranostics 2021, 11, 649–664. [Google Scholar] [CrossRef]
  25. Zaki, A.M.; Van Boheemen, S.; Bestebroer, T.M.; Osterhaus, A.D.M.E.; Fouchier, R.A.M. Isolation of a Novel Coronavirus from a Man with Pneumonia in Saudi Arabia. N. Engl. J. Med. 2012, 367, 1814–1820. [Google Scholar] [CrossRef]
  26. Li, W.; Shi, Z.; Yu, M.; Ren, W.; Smith, C.; Epstein, J.H.; Wang, H.; Crameri, G.; Hu, Z.; Zhang, H.; et al. Bats Are Natural Reservoirs of SARS-Like Coronaviruses. Science 2005, 310, 676–679. [Google Scholar] [CrossRef] [PubMed]
  27. Guan, Y.J.; Zheng, B.J.; He, Y.Q.; Liu, X.L.; Zhuang, Z.X.; Cheung, C.L.; Luo, S.W.; Li, P.H.; Zhang, L.J.; Butt, K.M.; et al. Isolation and Characterization of Viruses Related to the SARS Coronavirus from Animals in Southern China. Science 2003, 302, 276–278. [Google Scholar] [CrossRef] [Green Version]
  28. De Buhr, H.; Lebbink, R.J. Harnessing CRISPR to combat human viral infections. Curr. Opin. Immunol. 2018, 54, 123–129. [Google Scholar] [CrossRef] [PubMed]
  29. Adli, M. The CRISPR tool kit for genome editing and beyond. Nat. Commun. 2018, 9, 1911. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Freije, C.A.; Sabeti, P.C. Detect and destroy: CRISPR-based technologies for the response against viruses. Cell Host Microbe 2021, 29, 689–703. [Google Scholar] [CrossRef]
  31. O’Connell, M.R. Molecular Mechanisms of RNA Targeting by Cas13-containing Type VI CRISPR–Cas Systems. J. Mol. Biol. 2018, 431, 66–87. [Google Scholar] [CrossRef]
  32. Konermann, S.; Lotfy, P.; Brideau, N.J.; Oki, J.; Shokhirev, M.N.; Hsu, P.D. Transcriptome Engineering with RNA-Targeting Type VI-D CRISPR Effectors. Cell 2018, 173, 665–676.e14. [Google Scholar] [CrossRef] [Green Version]
  33. Yan, W.X.; Chong, S.; Zhang, H.; Makarova, K.S.; Koonin, E.V.; Cheng, D.R.; Scott, D.A. Cas13d Is a Compact RNA-Targeting Type VI CRISPR Effector Positively Modulated by a WYL-Domain-Containing Accessory Protein. Mol. Cell 2018, 70, 327–339.e5. [Google Scholar] [CrossRef] [Green Version]
  34. Zhang, C.; Konermann, S.; Brideau, N.J.; Lotfy, P.; Wu, X.B.; Novick, S.J.; Strutzenberg, T.; Griffin, P.R.; Hsu, P.D.; Lyumkis, D. Structural Basis for the RNA-Guided Ribonuclease Activity of CRISPR-Cas13d. Cell 2018, 175, 212–223.e17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Lotfi, M.; Rezaei, N. CRISPR/Cas13: A potential therapeutic option of COVID-19. Biomed. Pharmacother. 2020, 131, 110738. [Google Scholar] [CrossRef] [PubMed]
  36. Abudayyeh, O.O.; Gootenberg, J.S.; Essletzbichler, P.; Han, S.; Joung, J.; Belanto, J.J.; Verdine, V.; Cox, D.B.T.; Kellner, M.J.; Regev, A.; et al. RNA targeting with CRISPR-Cas13. Nature 2017, 550, 280–284. [Google Scholar] [CrossRef] [Green Version]
  37. Cox, D.B.T.; Gootenberg, J.S.; Abudayyeh, O.O.; Franklin, B.; Kellner, M.J.; Joung, J.; Zhang, F. RNA editing with CRISPR-Cas13. Science 2017, 358, 1019–1027. [Google Scholar] [CrossRef] [Green Version]
  38. Marra, M.A.; Jones, S.J.M.; Astell, C.R.; Holt, R.A.; Brooks-Wilson, A.; Butterfield, Y.S.N.; Khattra, J.; Asano, J.K.; Barber, S.A.; Chan, S.Y.; et al. The Genome Sequence of the SARS-Associated Coronavirus. Science 2003, 300, 1399–1404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. ter Brake, O.; von Eije, K.J.; Berkhout, B. Probing the sequence space available for HIV-1 evolution. Aids 2008, 22, 1875–1877. [Google Scholar] [CrossRef] [PubMed]
  40. Plant, E.P.; Sims, A.C.; Baric, R.S.; Dinman, J.D.; Taylor, D.R. Altering SARS Coronavirus Frameshift Efficiency Affects Genomic and Subgenomic RNA Production. Viruses 2013, 5, 279–294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Ahn, D.-G.; Lee, W.; Choi, J.-K.; Kim, S.-J.; Plant, E.P.; Almazán, F.; Taylor, D.R.; Enjuanes, L.; Oh, J.-W. Interference of ribosomal frameshifting by antisense peptide nucleic acids suppresses SARS coronavirus replication. Antivir. Res. 2011, 91, 1–10. [Google Scholar] [CrossRef]
  42. Ren, P.X.; Shang, W.-j.; Yin, W.-c.; Ge, H.; Wang, L.; Zhang, X.-l.; Li, B.-q.; Li, H.-l.; Xu, Y.-c.; Xu, E.H.; et al. A multi-targeting drug design strategy for identifying potent anti-SARS-CoV-2 inhibitors. Acta Pharmacol. Sin. 2021, 43, 483–493. [Google Scholar] [CrossRef] [PubMed]
  43. Wessels, H.-H.; Méndez-Mancilla, A.; Guo, X.; Legut, M.; Daniloski, Z.; Sanjana, N.E. Massively parallel Cas13 screens reveal principles for guide RNA design. Nat. Biotechnol. 2020, 38, 722–727. [Google Scholar] [CrossRef] [PubMed]
  44. Schmitt, A.; Herzel, H. Estimating the Entropy of DNA Sequences. J. Theor. Biol. 1997, 188, 369–377. [Google Scholar] [CrossRef]
  45. Shannon, C.E. The mathematical theory of communication. 1963. MD Comput. 1997, 14, 306–317. [Google Scholar]
  46. Malone, B.; Urakova, N.; Snijder, E.J.; Campbell, E.A. Structures and functions of coronavirus replication-transcription complexes and their relevance for SARS-CoV-2 drug design. Nat. Rev. Mol. Cell Biol. 2022, 23, 21–39. [Google Scholar] [CrossRef] [PubMed]
  47. Wang, D.; Wandel, M.P.; Bandyra, K.J.; Ulferts, R.; Wu, M.; Dendooven, T.; Yang, X.; Doyle, N.; Oerum, S.; Beale, R.; et al. The SARS-CoV-2 subgenome landscape and its novel regulatory features. Mol. Cell. 2021, 81, 2135–2147.e5. [Google Scholar] [CrossRef]
  48. Robertson, M.P.; Igel, H.; Baertsch, R.; Haussler, D.; Ares, M.; Scott, W.G. The Structure of a Rigorously Conserved RNA Element within the SARS Virus Genome. PLoS Biol. 2004, 3, e5. [Google Scholar] [CrossRef] [PubMed]
  49. Lulla, V.; Wandel, M.P.; Bandyra, K.J.; Ulferts, R.; Wu, M.; Dendooven, T.; Yang, X.; Doyle, N.; Oerum, S.; Beale, R.; et al. Targeting the Conserved Stem Loop 2 Motif in the SARS-CoV-2 Genome. J. Virol. 2021, 95, e0066321. [Google Scholar] [CrossRef]
  50. Westerhout, E.M.; Ooms, M.; Vink, M.; Das, A.T.; Berkhout, B. HIV-1 can escape from RNA interference by evolving an alternative structure in its RNA genome. Nucleic Acids Res. 2005, 33, 796–804. [Google Scholar] [CrossRef] [Green Version]
  51. Almazan, F.; DeDiego, M.L.; Sola, I.; Zuñiga, S.; Nieto-Torres, J.L.; Márquez-Jurado, S.; Andrés, G.; Enjuanes, L. Engineering a Replication-Competent, Propagation-Defective Middle East Respiratory Syndrome Coronavirus as a Vaccine Candidate. mBio 2013, 4, e00650-13. [Google Scholar] [CrossRef] [Green Version]
  52. Almazán, F.; Galán, C.; Enjuanes, L. The Nucleoprotein Is Required for Efficient Coronavirus Genome Replication. J. Virol. 2004, 78, 12683–12688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Almazán, F.; DeDiego, M.L.; Galán, C.; Escors, D.; Álvarez, E.; Ortego, J.; Sola, I.; Zuñiga, S.; Alonso, S.; Moreno, J.L.; et al. Construction of a Severe Acute Respiratory Syndrome Coronavirus Infectious cDNA Clone and a Replicon To Study Coronavirus RNA Synthesis. J. Virol. 2006, 80, 10900–10906. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Scientific, I.T.F. Lipofectamine 2000 Transfection Reagent Protocol. Available online: https://www.thermofisher.com/document-connect/document-connect.html?url=https://assets.thermofisher.com/TFS-Assets%2FLSG%2Fmanuals%2FLipofectamine_2000_Reag_protocol.pdf (accessed on 3 March 2023).
  55. Promega. Dual-Luciferase® Reporter Assay System Technical Manual. Available online: https://nld.promega.com/products/luciferase-assays/reporter-assays/dual_luciferase-reporter-assay-system/?catNum=E1910 (accessed on 3 March 2023).
  56. Ruijter, J.M.; Thygesen, H.H.; Schoneveld, O.J.; Das, A.T.; Berkhout, B.; Lamers, W.H. Factor correction as a tool to eliminate between-session variation in replicate experiments: Application to molecular biology and retrovirology. Retrovirology 2006, 3, 2. [Google Scholar] [CrossRef] [Green Version]
  57. Feng, X.; Zhang, X.; Jiang, S.; Tang, Y.; Cheng, C.; Krishna, P.A.; Wang, X.; Dai, J.; Zeng, J.; Xia, T.; et al. A DNA-based non-infectious replicon system to study SARS-CoV-2 RNA synthesis. Comput. Struct. Biotechnol. J. 2022, 20, 5193–5202. [Google Scholar] [CrossRef]
  58. Wang, G.; Zhao, N.; Berkhout, B.; Das, A.T. A Combinatorial CRISPR-Cas9 Attack on HIV-1 DNA Extinguishes All Infectious Provirus in Infected T Cell Cultures. Cell Rep. 2016, 17, 2819–2826. [Google Scholar] [CrossRef] [Green Version]
  59. Darcis, G.; Binda, C.S.; Klaver, B.; Herrera-Carrillo, E.; Berkhout, B.; Das, A.T. The Impact of HIV-1 Genetic Diversity on CRISPR-Cas9 Antiviral Activity and Viral Escape. Viruses 2019, 11, 255. [Google Scholar] [CrossRef] [Green Version]
  60. Herrera-Carrillo, E.; Berkhout, B. The impact of HIV-1 genetic diversity on the efficacy of a combinatorial RNAi-based gene therapy. Gene Ther. 2015, 22, 485–495. [Google Scholar] [CrossRef]
  61. ter Brake, O.; Hooft, K.; Liu, Y.P.; Centlivre, M.; von Eije, K.J.; Berkhout, B. Lentiviral Vector Design for Multiple shRNA Expression and Durable HIV-1 Inhibition. Mol. Ther. 2008, 16, 557–564. [Google Scholar] [CrossRef]
  62. Herrera-Carrillo, E.; Berkhout, B. Gene Therapy Strategies to Block HIV-1 Replication by RNA Interference. Gene Ther. HIV Chronic Infect. 2015, 848, 71–95. [Google Scholar] [CrossRef]
  63. Herrera-Carrillo, E.; Berkhout, B. Attacking HIV-1 RNA versus DNA by sequence-specific approaches: RNAi versus CRISPR-Cas. Biochem. Soc. Trans. 2016, 44, 1355–1365. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Herrera-Carrillo, E.; Harwig, A.; Berkhout, B. Silencing of HIV-1 by AgoshRNA molecules. Gene Ther. 2017, 24, 453–461. [Google Scholar] [CrossRef] [PubMed]
  65. Ambike, S.; Cheng, C.-C.; Feuerherd, M.; Velkov, S.; Baldassi, D.; Afridi, S.Q.; Porras-Gonzalez, D.; Wei, X.; Hagen, P.; Kneidinger, N.; et al. Targeting genomic SARS-CoV-2 RNA with siRNAs allows efficient inhibition of viral replication and spread. Nucleic Acids Res. 2022, 50, 333–349. [Google Scholar] [CrossRef] [PubMed]
  66. Kelly, J.A.; Olson, A.N.; Neupane, K.; Munshi, S.; Emeterio, J.S.; Pollack, L.; Woodside, M.T.; Dinman, J.D. Structural and functional conservation of the programmed -1 ribosomal frameshift signal of SARS coronavirus 2 (SARS-CoV-2). J. Biol. Chem. 2020, 295, 10741–10748. [Google Scholar] [CrossRef] [PubMed]
  67. Plant, E.P.; Pérez-Alvarado, G.C.; Jacobs, J.L.; Mukhopadhyay, B.; Hennig, M.; Dinman, J.D. A Three-Stemmed mRNA Pseudoknot in the SARS Coronavirus Frameshift Signal. PLoS Biol. 2005, 3, e172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Viruses 15 00686 g001 550
Figure 1. Replication cycle of SARS-CoV-2 and genome organization. (A) The SARS-CoV-2 replication cycle. Binding of the viral S protein to the ACE-2 receptor triggers SARS-CoV-2 infection and the plus genomic RNA (+gRNA) is translated into polyprotein (pp)1a and pp1ab. These proteins are proteolytically cleaved to generate 16 non-structural proteins (nsps), including some that form the replication/transcription complex that drives the synthesis of minus genomic RNA (−gRNA). The −gRNA is converted into genomic +gRNA that is packaged into new virion particles. Discontinuous transcription generates a set of 3′ co-terminal sub-genomic mRNAs (sgmRNAs) with identical 5′-leader and 3′-trailer ends. The sgmRNAs are translated into structural and accessory proteins that are required for the assembly of infectious virions, which takes place at the ER-Golgi intermediate compartment (ERGIC). Nascent virions are released from the cell via exocytosis. (B) The SARS-CoV-2 RNA genome. Human coronaviruses contain the largest viral genome (27–31 kb) among the RNA viruses and they share a similar genome organization. At the 5′-terminus two large overlapping open reading frames (ORFs), ORF 1a and ORF 1b, encode non-structural proteins and the 3′-terminal ORFs encode four structural proteins: spike (S), envelope (E), membrane (M) and nucleocapsid (N). The SARS-CoV-2 genome encodes eight accessory proteins: 3, 6, 7a, 7b, 8, 9, 14 and 10. Accessory proteins are indicated in grey.
Figure 1. Replication cycle of SARS-CoV-2 and genome organization. (A) The SARS-CoV-2 replication cycle. Binding of the viral S protein to the ACE-2 receptor triggers SARS-CoV-2 infection and the plus genomic RNA (+gRNA) is translated into polyprotein (pp)1a and pp1ab. These proteins are proteolytically cleaved to generate 16 non-structural proteins (nsps), including some that form the replication/transcription complex that drives the synthesis of minus genomic RNA (−gRNA). The −gRNA is converted into genomic +gRNA that is packaged into new virion particles. Discontinuous transcription generates a set of 3′ co-terminal sub-genomic mRNAs (sgmRNAs) with identical 5′-leader and 3′-trailer ends. The sgmRNAs are translated into structural and accessory proteins that are required for the assembly of infectious virions, which takes place at the ER-Golgi intermediate compartment (ERGIC). Nascent virions are released from the cell via exocytosis. (B) The SARS-CoV-2 RNA genome. Human coronaviruses contain the largest viral genome (27–31 kb) among the RNA viruses and they share a similar genome organization. At the 5′-terminus two large overlapping open reading frames (ORFs), ORF 1a and ORF 1b, encode non-structural proteins and the 3′-terminal ORFs encode four structural proteins: spike (S), envelope (E), membrane (M) and nucleocapsid (N). The SARS-CoV-2 genome encodes eight accessory proteins: 3, 6, 7a, 7b, 8, 9, 14 and 10. Accessory proteins are indicated in grey.
Viruses 15 00686 g001
Viruses 15 00686 g002 550
Figure 2. Targeting of SARS-CoV-2 RNA with CRISPR/Cas13d. (A) SARS-CoV-2 RNA genome diversity and the crRNA targets. Distribution of crRNAs along the SARS-CoV-2 RNA genome (MN908947). The Shannon entropy along the RNA genome varies from 0 to 1, where lower genetic diversity gives values closer to zero.The position of the crRNAs targeting conserved RNA sequences of SARS-CoV-2 is indicated with black triangles. The SARS-CoV-2 RNA genome diversity is shown by plotting Shannon entropy along the whole genome. The values vary between 0 and 1, where more conserved sequences give lower values (closer to 0); (B) Luciferase multitarget reporter construct. To measure the efficiency of the designed crRNAs, a pGL3 control plasmid-based multi-target luciferase reporter construct was designed with SARS-CoV-2 conserved target sequences cloned downstream of the firefly luciferase gene. Two constructs were designed: one with the + and one with the −RNA target sequences. SV40: simian virus 40 promoter, pA: polyA sequence; (C) Engineering a SARS-CoV-2 replicon. Scheme of the SARS-CoV-2 cDNA cloned in a BAC (upper panel) and the SARS-CoV-2 replicon (lower panel). The letters indicate the viral genes: S, spike; E, envelope; M, membrane; N, nucleocapsid; or the reporter gene mNeonGreen (mNG). UTR, untranslated regions; CMV, cytomegalovirus promoter; pA, polyA sequence; Rz, hepatitis delta virus ribozyme; BGH, bovine growth hormone polyadenylation and termination signals.
Figure 2. Targeting of SARS-CoV-2 RNA with CRISPR/Cas13d. (A) SARS-CoV-2 RNA genome diversity and the crRNA targets. Distribution of crRNAs along the SARS-CoV-2 RNA genome (MN908947). The Shannon entropy along the RNA genome varies from 0 to 1, where lower genetic diversity gives values closer to zero.The position of the crRNAs targeting conserved RNA sequences of SARS-CoV-2 is indicated with black triangles. The SARS-CoV-2 RNA genome diversity is shown by plotting Shannon entropy along the whole genome. The values vary between 0 and 1, where more conserved sequences give lower values (closer to 0); (B) Luciferase multitarget reporter construct. To measure the efficiency of the designed crRNAs, a pGL3 control plasmid-based multi-target luciferase reporter construct was designed with SARS-CoV-2 conserved target sequences cloned downstream of the firefly luciferase gene. Two constructs were designed: one with the + and one with the −RNA target sequences. SV40: simian virus 40 promoter, pA: polyA sequence; (C) Engineering a SARS-CoV-2 replicon. Scheme of the SARS-CoV-2 cDNA cloned in a BAC (upper panel) and the SARS-CoV-2 replicon (lower panel). The letters indicate the viral genes: S, spike; E, envelope; M, membrane; N, nucleocapsid; or the reporter gene mNeonGreen (mNG). UTR, untranslated regions; CMV, cytomegalovirus promoter; pA, polyA sequence; Rz, hepatitis delta virus ribozyme; BGH, bovine growth hormone polyadenylation and termination signals.
Viruses 15 00686 g002
Viruses 15 00686 g003 550
Figure 3. The efficiency of the designed anti-SARS-CoV-2 crRNAs. (A,B) The efficiency of 29 crRNAs targeting either the + or the -RNA strand with Cas13d was tested in HEK293T cells by co-transfecting them with the luciferase construct encoding SARS-CoV-2 +RNA (A) and −RNA (B) target sequences, respectively. To quantify viral gene expression, luciferase level was measured at two days after transfection. Mean values (±SD) of three experiments in duplicate are shown. The average luciferase levels are expressed as relative percentage level (%), setting the negative control (NC), a crRNA not targeting any SARS-CoV-2 encoding sequence, at 100%. Statistical significance was determined using one-way ANOVA, * p ≤ 0.05.
Figure 3. The efficiency of the designed anti-SARS-CoV-2 crRNAs. (A,B) The efficiency of 29 crRNAs targeting either the + or the -RNA strand with Cas13d was tested in HEK293T cells by co-transfecting them with the luciferase construct encoding SARS-CoV-2 +RNA (A) and −RNA (B) target sequences, respectively. To quantify viral gene expression, luciferase level was measured at two days after transfection. Mean values (±SD) of three experiments in duplicate are shown. The average luciferase levels are expressed as relative percentage level (%), setting the negative control (NC), a crRNA not targeting any SARS-CoV-2 encoding sequence, at 100%. Statistical significance was determined using one-way ANOVA, * p ≤ 0.05.
Viruses 15 00686 g003
Viruses 15 00686 g004 550
Figure 4. Titration of SARS-CoV-2 targeting crRNAs. The efficiency of 14 crRNAs was tested in HEK293T cells by co-transfecting them with the luciferase construct encoding SARS-CoV-2 +RNA target sequences. A titration of the crRNA constructs was performed (75, 150 and 300 ng). To quantify viral gene expression, luciferase level was measured at two days after transfection. Mean values (±SD) of three experiments in duplicate are shown. The average luciferase levels are expressed as percentage (%) of luciferase expression, setting the negative control (NC) at 100%. Statistical significance was determined using one-way ANOVA, * p ≤ 0.05.
Figure 4. Titration of SARS-CoV-2 targeting crRNAs. The efficiency of 14 crRNAs was tested in HEK293T cells by co-transfecting them with the luciferase construct encoding SARS-CoV-2 +RNA target sequences. A titration of the crRNA constructs was performed (75, 150 and 300 ng). To quantify viral gene expression, luciferase level was measured at two days after transfection. Mean values (±SD) of three experiments in duplicate are shown. The average luciferase levels are expressed as percentage (%) of luciferase expression, setting the negative control (NC) at 100%. Statistical significance was determined using one-way ANOVA, * p ≤ 0.05.
Viruses 15 00686 g004
Viruses 15 00686 g005 550
Figure 5. The anti-SARS-CoV-2 efficiency of single versus dual crRNA therapy. The efficiency of six crRNAs was tested in HEK293T cells by co-transfecting them with the luciferase construct encoding SARS-CoV-2 +RNA target sequences. The crRNAs were tested as single inhibitor versus a dual-combinatorial approach. To quantify viral gene expression, luciferase level was measured at two days after transfection. Mean values (±SD) of three experiments in duplicates are shown. The average luciferase levels are expressed as percentage (%) of luciferase expression, setting the negative control (NC) at 100%. Statistical significance was determined using one-way ANOVA, **** p ≤ 0.0001.
Figure 5. The anti-SARS-CoV-2 efficiency of single versus dual crRNA therapy. The efficiency of six crRNAs was tested in HEK293T cells by co-transfecting them with the luciferase construct encoding SARS-CoV-2 +RNA target sequences. The crRNAs were tested as single inhibitor versus a dual-combinatorial approach. To quantify viral gene expression, luciferase level was measured at two days after transfection. Mean values (±SD) of three experiments in duplicates are shown. The average luciferase levels are expressed as percentage (%) of luciferase expression, setting the negative control (NC) at 100%. Statistical significance was determined using one-way ANOVA, **** p ≤ 0.0001.
Viruses 15 00686 g005
Viruses 15 00686 g006 550
Figure 6. The efficiency of crRNAs targeting SARS-CoV-2 and SARS-CoV replicon. (A) The efficiency of five selected single crRNAs targeting SARS-CoV-2 sequences was tested in HEK293T cells by co-transfecting them with a SARS-CoV-2 replicon. crRNAs targeting both SARS-CoV-2 + and −RNA sequences were included. Viral RNA replication was quantified by measuring levels of mNeonGreen. Mean values (±SD) of three experiments in duplicates are shown. The average mNeonGreen levels are expressed as percentage (%) of fluorescent cells, setting the negative control (NC) at 100%; (B) The same experimental set-up was used with a SARS-CoV replicon. Here, viral RNA replication was quantified by measuring levels of GFP. Statistical significance was determined using one-way ANOVA, **** p ≤ 0.0001.
Figure 6. The efficiency of crRNAs targeting SARS-CoV-2 and SARS-CoV replicon. (A) The efficiency of five selected single crRNAs targeting SARS-CoV-2 sequences was tested in HEK293T cells by co-transfecting them with a SARS-CoV-2 replicon. crRNAs targeting both SARS-CoV-2 + and −RNA sequences were included. Viral RNA replication was quantified by measuring levels of mNeonGreen. Mean values (±SD) of three experiments in duplicates are shown. The average mNeonGreen levels are expressed as percentage (%) of fluorescent cells, setting the negative control (NC) at 100%; (B) The same experimental set-up was used with a SARS-CoV replicon. Here, viral RNA replication was quantified by measuring levels of GFP. Statistical significance was determined using one-way ANOVA, **** p ≤ 0.0001.
Viruses 15 00686 g006
Viruses 15 00686 g007 550
Figure 7. The impact of Cas13 on the intermediate minus sense transcripts. (A) Schematic representation of differences between +RNA and –RNA strands during SARS-CoV-2 RNA replication. +RNA strands (black arrows) are far more abudant than the –RNA strands (black dashed arrows). The less abundant –RNA strands are present as double-stranded replicative forms during RNA replication. (B) The efficiency of crRNAs targeting -RNA targets when annealed to complementary +RNA targets. The efficiency of six crRNAs was tested in HEK293T cells by co-transfecting them with a mixture of different ratios of − and + luciferase constructs (1:5, 1:25, 1:50, 1:75 and 1:100 of − luciferase construct: + luciferase construct) or with the − luciferase construct alone (1:0). To quantify viral gene expression, luciferase level was measured at two days after transfection. Mean values (±SD) of three experiments in duplicates are shown. The average luciferase levels are expressed as percentage (%) of luciferase expression, setting the negative control (NC) at 100%. Statistical significance was determined using one-way ANOVA, * p ≤ 0.05.
Figure 7. The impact of Cas13 on the intermediate minus sense transcripts. (A) Schematic representation of differences between +RNA and –RNA strands during SARS-CoV-2 RNA replication. +RNA strands (black arrows) are far more abudant than the –RNA strands (black dashed arrows). The less abundant –RNA strands are present as double-stranded replicative forms during RNA replication. (B) The efficiency of crRNAs targeting -RNA targets when annealed to complementary +RNA targets. The efficiency of six crRNAs was tested in HEK293T cells by co-transfecting them with a mixture of different ratios of − and + luciferase constructs (1:5, 1:25, 1:50, 1:75 and 1:100 of − luciferase construct: + luciferase construct) or with the − luciferase construct alone (1:0). To quantify viral gene expression, luciferase level was measured at two days after transfection. Mean values (±SD) of three experiments in duplicates are shown. The average luciferase levels are expressed as percentage (%) of luciferase expression, setting the negative control (NC) at 100%. Statistical significance was determined using one-way ANOVA, * p ≤ 0.05.
Viruses 15 00686 g007
Table
Table 1. crRNAs targeting SARS-CoV-2 RNA.
Table 1. crRNAs targeting SARS-CoV-2 RNA.
NameNucleotide PositioncrRNA SequenceTarget SequenceGuide Score
5’UTR-144–66GAGAACAGATCTACAAGAGATCGCGATCTCTTGTAGATCTGTTCTC0.46
5’UTR-253–75GTTCGTTTAGAGAACAGATCTACGTAGATCTGTTCTCTAAACGAAC0.45
nsp3-15666–5688GACATCATAACAAAAGGTGACTCGAGTCACCTTTTGTTATGATGTC0.89
nsp3-27529–7551CCATTAACAATAGTTGTACATTCGAATGTACAACTATTGTTAATGG0.42
nsp4-19584–9606AATGTCAAGTACAAGTAAATAACGTTATTTACTTGTACTTGACATT0.13
nsp4-29901–9923CTCCACTAAAATACTTGTACTTATAAGTACAAGTATTTTAGTGGAG0.15
nsp6-111603–11625CCTAAGAAACAATAAACTAGCATATGCTAGTTTATTGTTTCTTAGG0.38
nsp6-211675–11697TAAACACCAAGAGTCAGTCTAAATTTAGACTGACTCTTGGTGTTTA0.24
slippery13463–13485ACTTACACCGCAAACCCGTTTAATTAAACGGGTTTGCGGTGTAAGT0.37
frameshift 13468–13490GCTGCACTTACACCGCAAACCCGCGGGTTTGCGGTGTAAGTGCAGC0.68
RdRP-113769–13791TTGACGTGATATATGTGGTACCATGGTACCACATATATCACGTCAA0.77
RdRP-214507–14529GCTATGTAAGTTTACATCCTGATATCAGGATGTAAACTTACATAGC0.64
RdRP-315052–15074TTAAGATTCATTTGAGTTATAGTACTATAACTCAAATGAATCTTAA0.21
RdRP-415452–15474ACCTGGTTTAACATATAGTGAACGTTCACTATATGTTAAACCAGGT0.29
RdRP-515622–15644CTATTTCTATAGAGACACTCATATATGAGTGTCTCTATAGAAATAG0.35
Helicase-117011–17033ACATTGCTAGAAAACTCATCTGATCAGATGAGTTTTCTAGCAATGT0.51
Helicase-217221–17243GCAGGTATAATTCTACTACATTTAAATGTAGTAGAATTATACCTGC0.22
Helicase-317479–17501AAATATTCTGGTTCTAGTGTGCCGGCACACTAGAACCAGAATATTT0.44
Helicase-416348–16370ACAGACAAGACTAATTTATGTGATCACATAAATTAGTCTTGTCTGT0.17
N-128409–28431AACCAAGACGCAGTATTATTGGGCCCAATAATACTGCGTCTTGGTT0.25
N-228434–28456CTTGCCATGTTGAGTGAGAGCGGCCGCTCTCACTCAACATGGCAAG0.61
N-328513–28535GGTAGTAGCCAATTTGGTCATCTAGATGACCAAATTGGCTACTACC0.58
N-428544–28566TCACCACCACGAATTCGTCTGGTACCAGACGAATTCGTGGTGGTGA0.59
N-529096–29118GTTTGTTCTGGACCACGTCTGCCGGCAGACGTGGTCCAGAACAAAC0.51
3’Pseudoknot-129543–29565TAGCCCATCTGCCTTGTGTGGTCGACCACACAAGGCAGATGGGCTA0.43
3’Pseudoknot-229586–29608GAGTAGACTATATATCGTAAACGCGTTTACGATATATAGTCTACTC0.41
3’Pseudoknot-329638–29660AGTTAACTACATCTACTTGTGCTAGCACAAGTAGATGTAGTTAACT0.38
s2m29742–29764CTGTACACTCGATCGTACTCCGCGCGGAGTACGATCGAGTGTACAG0.62
3’UTR29787–29809ACATTAGGGCTCTTCCATATAGGCCTATATGGAAGAGCCCTAATGT0.38
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hussein, M.; Andrade dos Ramos, Z.; Vink, M.A.; Kroon, P.; Yu, Z.; Enjuanes, L.; Zuñiga, S.; Berkhout, B.; Herrera-Carrillo, E. Efficient CRISPR-Cas13d-Based Antiviral Strategy to Combat SARS-CoV-2. Viruses 2023, 15, 686. https://doi.org/10.3390/v15030686

AMA Style

Hussein M, Andrade dos Ramos Z, Vink MA, Kroon P, Yu Z, Enjuanes L, Zuñiga S, Berkhout B, Herrera-Carrillo E. Efficient CRISPR-Cas13d-Based Antiviral Strategy to Combat SARS-CoV-2. Viruses. 2023; 15(3):686. https://doi.org/10.3390/v15030686

Chicago/Turabian Style

Hussein, Mouraya, Zaria Andrade dos Ramos, Monique A. Vink, Pascal Kroon, Zhenghao Yu, Luis Enjuanes, Sonia Zuñiga, Ben Berkhout, and Elena Herrera-Carrillo. 2023. "Efficient CRISPR-Cas13d-Based Antiviral Strategy to Combat SARS-CoV-2" Viruses 15, no. 3: 686. https://doi.org/10.3390/v15030686

APA Style

Hussein, M., Andrade dos Ramos, Z., Vink, M. A., Kroon, P., Yu, Z., Enjuanes, L., Zuñiga, S., Berkhout, B., & Herrera-Carrillo, E. (2023). Efficient CRISPR-Cas13d-Based Antiviral Strategy to Combat SARS-CoV-2. Viruses, 15(3), 686. https://doi.org/10.3390/v15030686

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop