Citation: Ramulifho, E.; Rey, C. A Coiled-Coil Nucleotide-Binding Domain Leucine-Rich Repeat Receptor Gene MeRPPL1 Plays a Role in the Replication of a Geminivirus in Cassava. Viruses 2024, 16, 941. https://doi.org/10.3390/v16060941 Academic Editor: Francisco Murilo Zerbini Received: 19 January 2024 Revised: 14 May 2024 Accepted: 30 May 2024 Published: 11 June 2024 Copyright: © 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/). viruses Article A Coiled-Coil Nucleotide-Binding Domain Leucine-Rich Repeat Receptor Gene MeRPPL1 Plays a Role in the Replication of a Geminivirus in Cassava Elelwani Ramulifho 1,2 and Chrissie Rey 1,* 1 Plant Biotechnology Laboratory, School of Molecular and Cell Biology, University of the Witwatersrand, Johannesburg 2001, South Africa; ramulifhoelelwani@gmail.com 2 Germplasm Development, Agricultural Research Council, Small Grain Institute, Bethlehem 9700, South Africa * Correspondence: chrissie.rey@wits.ac.za Abstract: Disease resistance gene (R gene)-encoded nucleotide-binding leucine-rich repeat proteins (NLRs) are critical players in plant host defence mechanisms because of their role as receptors that recognise pathogen effectors and trigger plant effector-triggered immunity (ETI). This study aimed to determine the putative role of a cassava coiled-coil (CC)-NLR (CNL) gene MeRPPL1 (Manes.12G091600) (single allele) located on chromosome 12 in the tolerance or susceptibility to South African cassava mosaic virus (SACMV), one of the causal agents of cassava mosaic disease (CMD). A transient protoplast system was used to knock down the expression of MeRPPL1 by clustered regularly interspaced short palindromic repeats-CRISPR-associated protein 9 (CRISPR-Cas9). The MeRPPL1-targeting CRISPR vectors and/or SACMV DNA A and DNA B infectious clones were used to transfect protoplasts isolated from leaf mesophyll cells from the SACMV-tolerant cassava (Manihot esculenta) cultivar TME3. The CRISPR/Cas9 silencing vector significantly reduced MeRPPL1 expression in protoplasts whether with or without SACMV co-infection. Notably, SACMV DNA A replication was higher in protoplasts with lower MeRPPL1 expression levels than in non-silenced protoplasts. Mutagenesis studies revealed that protoplast co-transfection with CRISPR-MeRPPL1 silencing vector + SACMV and transfection with only SACMV induced nucleotide substitution mutations that led to altered amino acids in the highly conserved MHD motif of the MeRPPL1- translated polypeptide. This may abolish or alter the regulatory role of the MHD motif in controlling R protein activity and could contribute to the increase in SACMV-DNA A accumulation observed in MeRPPL1-silenced protoplasts. The results herein demonstrate for the first time a role for a CNL gene in tolerance to a geminivirus in TME3. Keywords: NLR; geminivirus; SACMV; cassava mosaic disease; MeRPPL1; CRISPR/Cas9; resistance gene; TME3; tolerance 1. Introduction Cassava (Manihot esculenta Crantz) serves as both an industrial crop and a staple food source for more than 800 million people worldwide [1]. However, its production is threatened by viruses, belonging to the genus Begomovirus within the Geminiviridae family [2], that cause cassava mosaic disease (CMD) [3]. Cassava mosaic begomoviruses (CMBs) rank among the most damaging viruses, causing losses of at least 25 million tons annually in Africa, India, and Sri Lanka [4,5]. Amongst several CMB species is the South African cassava mosaic virus (SACMV), which is a bipartite (DNA A and B) circular ssDNA begomovirus [6]. CMD has effectively resulted in yield losses of up to 82% [4] by significantly reducing the number and size of tubers, depending on factors such as virus virulence, the susceptibility or tolerance/resistance of the host, and other environmental factors [7]. Infection presents as deformed/curly leaves with a yellow/pale-green mosaic and general reduced plant growth. Conventional strategies, such as phytosanitation, Viruses 2024, 16, 941. https://doi.org/10.3390/v16060941 https://www.mdpi.com/journal/viruses https://doi.org/10.3390/v16060941 https://creativecommons.org/ https://creativecommons.org/licenses/by/4.0/ https://creativecommons.org/licenses/by/4.0/ https://www.mdpi.com/journal/viruses https://www.mdpi.com https://orcid.org/0000-0002-7409-9790 https://doi.org/10.3390/v16060941 https://www.mdpi.com/journal/viruses https://www.mdpi.com/article/10.3390/v16060941?type=check_update&version=1 Viruses 2024, 16, 941 2 of 15 cultural practices, and whitefly vector control for controlling CMD have failed over the past years [5]. However, breeding programs have successfully generated cassava genotypes with some tolerance/resistance to CMD, exhibiting higher yields compared to susceptible varieties [8]. However, breeding programs have successfully generated cassava genotypes with some tolerance/resistance to CMD, exhibiting higher yields compared to susceptible genotypes [8]. Three cassava QTL markers (CMD1, CMD2, and CMD3) are associated with conferring resistance/partial resistance to CMD. Tropical M. esculenta (TME) [9,10] exhibits CMD2-associated monogenic dominant resistance. Landraces/cultivars with this marker display near-immunity phenotypes against most cassava geminiviruses [11]. Rabbi et al. [10] mapped the dominant monogenic resistance CMD2 locus and reported that all markers linked to qualitative resistance occur in the same chromosome region. Further genome-wide association studies (GWAS) predicted the location of genes associated with CMD2 on chromosome 8 and further identified two possible epistatic loci and multiple resistance alleles at the major CMD2 QTL [12]. The update of the cassava genome assembly from version Manihot esculenta v4.1 to v6.1/v8.1 relocated the CMD2 QTL on chromosome 12 (https://phytozome.jgi.doe.gov, accessed on 8 January 2024). While recent mapping studies identified the CMD2 QTL on chromosome 12, the actual gene(s) involved in triggering CMD2-associated immunity, and their associated molecular mechanisms, remain unknown. Plant immunity to viral pathogens comprises three major mechanisms, namely pathogen-triggered immunity (PTI), effector-triggered immunity (ETI) and innate immunity (RNA silencing) [13]. Monogenic dominant resistance is often associated with ETI triggered by resistance genes. These resistance (R) genes do not always but often encode nucleotide-binding leucine-rich repeat (NLR) protein receptors that are critical players as pathogen effector recognition receptors. NLRs perceive virus effectors and trigger the activation of signal transduction cascades resulting in ETI, suggesting their crucial role in defence response mechanisms [14,15]. Resistance genes associated specifically with plant viruses were reported in an overall review by Sett et al. [16]. Specifically, those induced by geminiviruses are CYR1, associated with mung bean yellow mosaic India virus (MYMIV) [17]; Ty-2, associated with tomato yellow leaf curl virus (TYLCV) [18,19]; Sw-5a, associated with tomato leaf curl New Delhi virus [20]; and Ty-6, associated with tomato yellow leaf curl virus and tomato mottle virus [21]. While a dominant monogenic resistance locus CMD2/CMD3 is reported for cassava mosaic geminiviruses [9], no specific genes specifically associated with this QTL have been confirmed in vivo. Intracellular R gene-translated proteins are represented mainly by the NLR receptors [16]. The natural Ty-2 gene that confers high resistance to the geminivirus tomato yellow leaf curl virus (TYLCV) encodes an NLR protein named TYNBS1, which activates ETI via the recognition of a Rep/C1 gene-encoded geminiviral effector [18,19]. More recently, a nonsynonymous single nucleotide polymorphism in DNA polymerase δ subunit 1 (MePOLD1) located within this CMD2 region has also been linked to CMD2-type resistance. Virus-induced gene silencing of MePOLD1 in a CMD-susceptible cassava variety produced a recovery phenotype typical of CMD2-type resistance [22]. Despite progress in cassava breeding for CMD resistance, specific cassava NLR genes involved in the molecular mechanisms of immunity/tolerance have not been identified or experimentally demonstrated. While several NLR resistance gene analogues (RGAs) were identified in seven African cassava landraces [23], no in planta functional studies on these were performed. The cassava transcriptome [24], proteome [25], and sequenced cassava genome data [26,27] revealed putative R genes/proteins that could be involved in the recognition of CMD-associated geminiviral pathogens. Furthermore, an NLR resistance gene with a single allele, MeRPPL1, recently relocated to cassava chromosome 12, was suggested to be involved in CMD tolerance in cassava cultivar TME3 [28]. From all the identified NLR genes with the potential to confer resistance to viruses, only a handful have been cloned and functionally characterised in different plants, including the orphan crop cassava. This research aimed to ascertain whether MeRPPL1 is associated with tolerance to SACMV in CMD-tolerant TME3 protoplasts. Clustered regularly interspaced short https://phytozome.jgi.doe.gov Viruses 2024, 16, 941 3 of 15 palindromic repeats (CRISPR) editing of MeRPPL1 was performed in an established cassava protoplast system. The results herein demonstrate, for the first time, a role for a resistance gene MeRPPL1 in a CMD-tolerant cassava host in response to a geminivirus. 2. Materials and Methods 2.1. MeRPPL1 Expression Knockdown gRNA Cassette Design and Synthesis Two gRNA targets were identified from the genomic sequence of Manes.12G091600 (herein referred to as MeRPPL1)genomic sequence using CRISPR version 4.7 [29]. The gR- NAs were used to design a duplex gRNA cassette in the following arrangement: Arabidopsis thaliana U6-26 promoter + pre-tRNA + crRNA 1 + SpCas9 tracrRNA + pre-tRNA + crRNA 2 + SpCas9 tracrRNA + terminator (Figure S1). The designed duplex silencing cassette was synthesised in vitro and cloned into a pBluescript II SK(+) plasmid at Inqaba Biotech (Pretoria, South Africa). 2.2. Cloning the CRISPR Cassettes into pCAMBIA 1380 Delivery Vector The gRNA cassette was digested from the pBluescript II SK(+) plasmid using re- striction enzymes HindIII and BglII and cloned into the delivery vector pCAMBIA 1380 to produce pC1380-MeRPPL1. A Cas9 cassette, consisting of Tobacco mosaic virus pro- moter + Cas9 + eGFP + Hsp terminator from the pL1m-f2-p35s-cas9-egfp-nucleo-thsp Golden Gate vector, was cloned to produce pC1380-Cas9 (Figure S2). Cloning was per- formed using a 5:1 (insert:vector) molar ratio. The cassettes were ligated into delivery vectors using the T4 DNA ligase enzyme (ThemoFisher Scientific, Waltham, MA, USA) according to the manufacturer’s manual, then transformed into chemically competent E. coli DH5-alpha cells using the snap freeze method. Restriction enzyme (HindIII and BglII) digestions and sequencing confirmed the positive ligation of the cassettes into the delivery vector. 2.3. CRISPR-Mediated MeRPPL1 Expression Knockdown: Protoplast Isolation and Transfection Methods with minor modifications were based on a previous study where the TME3 protoplast-SACMV system was optimised [30]. Protoplasts were isolated from 4-week-old leaf mesophyll cells of cassava TME3 plantlets propagated in vitro on Murashige and Skoog (MS) medium [31] by the shoot multiplication technique in sterile tissue culture bottles in environmentally controlled growth conditions (16:8 h (light: dark); 26–28 ◦C). Three independent biological replicate experiments were performed. For the first step, approximately 500 mg of fully expanded leaves were transversely sliced into thin strips and immediately plasmolysed in 10 mL CPW9M medium in order to facilitate subsequent cell wall degradation to release the protoplasts [30]. This CPW9M medium is composed of 0.5 M mannitol, 27.2 mg KH2PO4, 100 mg KNO3, 150 mg CaCl2, 250 mg MgSO4, 2.5 mg Fe2(SO4)3.6H20, 0.6 mg KI, 0.00025 mg CuSO4 per litre; pH 5.8. Incubation was performed for one hour at room temperature. CPW9M medium was replaced with 10 mL enzyme solution (5 mM MES, 2-(N-morpholino) ethanesulfonic acid, MES, 0.4% cellulase, 0.2% macerozyme) before vacuum infiltrating the leaf strips at 80 kPa for 30 min (step 2). The mixture was left to gently shake at 40 rpm for 24 h (overnight) in the dark at room temperature. Protoplasts were then released (step 3) by shaking at 80 rpm for 5 min, followed by filtering through a 75 µm sieve into a 50 mL round-bottom tube. Protoplasts were collected by centrifuging at 100 g for 10 min, washed twice with CPW9M medium, and resuspended in MMg solution (0.4 M mannitol, 15 mM MgCl2, 4 mM MES, pH 5.8) (step 4). The integrity and viability of protoplasts at 24 h after steps 3 and 4 were assessed using the Olympus BX 63 OM/FM microscope (Olympus Scientific Solutions, Waltham, MA, USA) and Evans Blue dye staining, respectively. Protoplasts were quantified using the BD AccuriTM C6 flow cytometer (BD Biosciences, Franklin Lakes, NJ, USA). One hundred and fifty protoplasts were counted for viability prior to treatments (at 24 h post-treatment; step 4) from the three independent biological experiments. On average, 135 out of 150 (90%) Viruses 2024, 16, 941 4 of 15 counted protoplasts were viable and intact. The ratio of viable to non-viable protoplasts was 90:10. Three independent biological experiments (n = 3) were performed, and for each biological experiment, three technical replicates were performed. For each treatment, the protoplasts were diluted to 104 protoplasts/mL. Protoplasts were transfected with (i) 15 µg of the pCAMBIA 1380-derived CRISPR vectors only (pC1380-MeRPPL1 + pC1380-Cas9); (ii) CRISPR vectors and 4 µg of each of the pBIN19-SACMV-DNA A and pBIN19-SACMV- DNA B infectious clones [32]; and (iii) SACMV-A and B infectious clones only. Control protoplasts were non-transfected or transfected with the pCAMBIA1380 vector only (mock). Previous well-established infectivity studies with SACMV infectious clones in planta [6] and a previous study in TME3 protoplasts [30] demonstrated that pBIN9 has no effect on cassava and was not included as an additional control. Polyethylene glycol 4000 (PEG4000) to a final concentration of 25% was immediately added to the transfected protoplasts and incubated at 28 ◦C for 20 min. The transfection mixture was then diluted with three volumes of W5 solution (154 mM NaCl, 125 mM CaCl2, 5 mM KCl, 2 mM MES, pH 5.8) and washed twice by centrifuging at 100× g for 2 min. Protoplasts were resuspended in 300 µL of WI solution (4 mM MES, 0.5 M mannitol, 20 mM KCl, pH 5.8) and incubated in the dark overnight at 28 ◦C to induce gene expression. Protoplasts were collected by centrifuging at 100× g for 5 min before being used in subsequent experiments. The viability of the protoplasts at 24 hpt was in a similar range for all transfections [(i) the CRISPR silencing vector only control (RC), (ii) SACMV (S), or both virus and CRISPR vector (RCS); pCAMBIA1350 only (mock)]. 2.4. RNA Extraction and RT-qPCR Gene Expression Analysis RNA was extracted from protoplasts at 24 hpt using QIAzol Lysis Reagent (QIAGEN) and cleaned using GeneJET RNA Cleanup and Concentration Micro Kit according to the manufacturer’s protocol. RevertAid First Strand cDNA Synthesis Kit (ThermoFisher Scientific) was used to synthesise first-strand cDNA using an Oligo(dT)18 primer, according to the manufacturer’s protocol. For relative MeRPPL1 expression, RT-qPCR was performed in triplicate with cDNA as a template using Maxima SYBR Green/ROX qPCR Master Mix (2X) with optimised MeRPPL1 primers (Fwd: 5′ agctcttgtgtatgggtgcc 3′, Rev: 5′ atggtgcccataatggaggc 3′). Log 2 fold change in the expression of MeRPPL1 was normalised to 18S rRNA gene [29]. Basal MeRPPL1 expression level was measured in non-transfected cassava TME3 protoplasts and protoplasts following transfection with SACMV, the CRISPR silencing vector, and the CRISPR + SACMV infectious clones. Expression was relative to mock, which was the pCAMBIA1380 vector only. Gene expression fold change was calculated using the log 2(∆∆CT) method [33]. Student’s t-test was used to assess statistical differences in MeRPPL1 expression between treatments. 2.5. DNA Extraction and Viral Replication Quantification DNA was extracted from transfected protoplasts using a User-Developed Protocol to Isolate Genomic DNA from the Interphase and Organic (phenol) Phase of Tissue Samples Treated with QIAzol Lysis Reagent (QIAGEN, Germantown, MD, USA). The DNA was cleaned up using Thermo Scientific GeneJET Gel Extraction and DNA Cleanup Micro Kit according to the manufacturer’s protocol before it was fast digested with the restriction enzyme DpnI (ThermoFisher Scientific, Massachusetts, USA). The quality of the Dpn-I digested DNA was visually checked by gel electrophoresis (1.5% agarose gel) to confirm that only newly replicated SACMV DNA A was measured and not input DNA from the transfection. Twenty-five nanograms of the digested DNA was used as a template for SACMV DNA A relative viral load quantification using qPCR. Quantification was performed in triplicates using Maxima SYBR Green/ROX qPCR Master Mix (2X) (ThermoFisher Scien- tific, Massachusetts, USA) with SACMV coat protein (CP) (on DNA A) primers (Fwd: 5′ acgtccgtcgcaagtac 3′, Rev: 5′ attgtcatgtcgaatagtacg 3′). As optimised in previous cassava Viruses 2024, 16, 941 5 of 15 studies [23,29], 18S rRNA was used as a reference gene for qPCR quantifications. The SACMV viral load was calculated using the log 2(∆∆CT) method [33]. Student’s t-test was used to assess statistical differences between samples. 2.6. Sequencing, Mutagenesis Analysis, and Protein Structure Prediction DpnI-digested DNA samples from the treatments (Section 2.3) were sent for the sequencing of MeRPPL1 at Inqaba Biotech (Pretoria, South Africa) using the ABI 3500XL Genetic Analyzer. SnapGene software was used to assess sequence variations (https://www.snapgene.com/, accessed on 8 January 2024), while MEGAX software (version 2.0) was used for multiple sequence alignment using the ClustalW algorithm [34]. Target gene homologs from wild- type AM560-2 (Accession NC_035172.1) and Manihot esculenta cultivar TME3 (RefSeq ID RSFT01000007, GenBank ID GCA_003957995.1, unpublished data) were used as references for the alignments. The frequency of clones with the altered sequence was obtained by expressing the number of amplicons from the polyclonal mix with sequence alteration as a ratio of 10 amplicons sequenced. AlphaFold.ipynb (https://colab.research.google.com/, accessed on 8 January 2024) [35] was used to predict the secondary structure of the proteins of translated amino acids from the sequenced DNA results. The monomer model was used for individual peptides. The run_relax option was selected to run AlphaFold and generate the predicted protein structures. 3. Results 3.1. MeRPPL1 gRNA Cassette Confirmation The MeRPPL1 gene is a single allele which was shown to reside in a locus on a chromosome 12 scaffold which is also associated with CMD2-dominant resistance in cas- sava (Figure S3). It was also selected from the transcriptomic data of the CMD-tolerant cassava cultivar TME3 infected with SACMV. Upregulation was observed at 32 days post- infection (dpi) [23], suggesting a putative role in response to SACMV. Phytozome (https:// phytozome-next.jgi.doe.gov/report/gene/Mesculenta_v8_1/Manes.12G091600, accessed on 8 January 2024) was accessed in January 2021, to retrieve the MeRPPL1 sequence [phy- tozome genome ID: 305; NCBI taxonomy ID: 3983 A Pathogen Receptor Genes (PRGs) database (PRGdb) version 3.0 (http://www.prgdb.org/prgdb/plants/?id=375, accessed on 8 January 2024)] to confirm that MeRPPL1 is a nucleotide-binding leucine-rich repeat (NLR) gene. The MeRPPL1 gRNA cassette (Figure 1A) was cloned into a pCAMBIA 1380 vector (Figure 1B) and confirmed by restriction enzyme digestion and Sanger sequencing. The vector was digested with HindIII and BglII, resulting in the expected band of approximately 760 bp by gel electrophoresis. 1 2910371 942 crRNA 1 crRNA 2 MeRPPL1 genomic sequence TTTTTTtracrRNA tracrRNA Arabidopsis U6-26 Promoter Pre-tRNA crRNA 1 gRNA 2 TerminatorcrRNA 2 gRNA cassette Pre-tRNA gRNA 1 394 bp 78 bp 20 bp 6 bp77 bp 78 bp 77 bp20 bp pCAMBIA 1380 HindIII BglII A B Figure 1. (A) MeRPPL1 gene (2910 nt), illustrating the positions (371 nt and 942 nt) for the design of the CRISPR silencing constructs crRNA1 and crRNA2. Below: Representation of the MeRPPL1 gRNA https://www.snapgene.com/ https://colab.research.google.com/ https://phytozome-next.jgi.doe.gov/report/gene/Mesculenta_v8_1/Manes.12G091600 https://phytozome-next.jgi.doe.gov/report/gene/Mesculenta_v8_1/Manes.12G091600 http://www.prgdb.org/prgdb/plants/?id=375 Viruses 2024, 16, 941 6 of 15 silencing cassette (760 bp): (L to R) the Arabidopsis U6-26 promoter; gRNA 1 and gRNA2 silencing constructs and terminator. Each gRNA is composed of a pre-tRNA, crRNA (20 bp), and tracrRNA (77 bp) region. The 5′ and 3 have a HindIII and BglI restriction site, respectively. (B) pCAMBIA 1380 vector (8861 bp) showing restriction sites HindIII and BglI where the MeRPPL1 gRNA cassette was cloned into. 3.2. Protoplast Isolation and Viability Protoplasts were isolated from 4-week-old leaf mesophyll cells of cassava TME3 plantlets propagated in vitro by the shoot multiplication technique in sterile tissue culture bottles in environmentally controlled conditions. When analysed under a light microscope, spherical intact and irregularly shaped protoplasts were observed (Figure 2A). The diam- eters between 16.1 and 38.9 µm were measured from at least ten spherical protoplasts. Chloroplasts were observed around the edge of protoplasts’ cell membranes. Before trans- fection, protoplasts were assessed for viability using Evans Blue dye, and at least 90% were viable (Figure 2B). On average, 135 out of 150 counted untreated protoplasts were viable and intact (ratio 90 viable:10 non-viable). The average viability of the protoplasts at 24 hpt was in a similar range (~90%) after all transfections [(i) pCAMBIA1350; SACMV only (pBIN19-SACMV-DNA A and pBIN19-SACMV-DNA B infectious clones); CRISPR vector pC1380-MeRPPL1 + pC1380-Cas9 only; and CRISPR vector + SACMV)]. The transfection rate (90%) was similar to a previous study in TME3 protoplasts [90%]. Viruses 2023, 15, x FOR PEER REVIEW 7 of 17 Figure 2. Analyses of the quality and quantity of protoplasts isolated from cassava TME3 landrace. (A) Cassava protoplasts isolated from leaf mesophyll cells. Red arrows show intact spherical pro- toplasts with chloroplasts around the edge of the cell membrane. (B) Evans Blue dye protoplast viability assessment on isolated protoplasts visualised under a light microscope. Viable cells do not take up the blue dye, while non-viable cells take up the blue dye due to compromised cell mem- branes (leaky). Hence, since the protoplasts do not appear blue (only the background appeared pale blue from the dye), this demonstrates that the protoplasts are viable. The scale bar on proto- plast images is 50 μm. (C) Protoplasts’ cell count and quality were measured using flow cytometry density measurement based on side scatter (SSC-A) and forward scatter (FSC-A). 3.3. Expression of MeRPPL1 in CRISPR/Cas9-Treated Protoplasts In this study, RT-qPCR was performed at 24 h post-transfection (hpt). To measure the relative expression levels of CRISPR-targeted MeRPPL1 in TME3 protoplasts, re- al-time qPCR was performed on cDNA using MeRPPL1 primers, with the 18S rRNA gene for normalisation. MeRPPL1 expression in non-transfected leaf protoplasts was low (log2fold −9.8). RT-qPCR comparisons in expression between transfection with CRISPR vectors only (RC), CRISPR vectors and SACMV (RCS), and SACMV only (S) were un- dertaken (Figure 3 top). The highest log2 fold change (3.1 log2 fold increase) of MeRPPL1 expression normalised to 18S rRNA was recorded in protoplasts that were only trans- fected with SACMV infectious clones (S), suggesting that SACMV triggers a MeRPPL1 response in SACMV-tolerant TME3-derived protoplasts. A decrease (log2 fold 1.5) in MeRPPL1 expression was observed in protoplasts co-transfected with both the CRISPR/Cas9 vectors and SACMV infectious clones (RCS) compared to protoplasts transfected only with SACMV (S) (Figure 3 top). While the ex- pression of MeRPPL1 was not abolished, its expression was significantly down-regulated by CRISPR. These results showed that SACMV may play a role in reducing/mitigating Figure 2. Analyses of the quality and quantity of protoplasts isolated from cassava TME3 landrace. (A) Cassava protoplasts isolated from leaf mesophyll cells. Red arrows show intact spherical proto- plasts with chloroplasts around the edge of the cell membrane. (B) Evans Blue dye protoplast viability assessment on isolated protoplasts visualised under a light microscope. Viable cells do not take up the blue dye, while non-viable cells take up the blue dye due to compromised cell membranes (leaky). Hence, since the protoplasts do not appear blue (only the background appeared pale blue from the Viruses 2024, 16, 941 7 of 15 dye), this demonstrates that the protoplasts are viable. The scale bar on protoplast images is 50 µm. (C) Protoplasts’ cell count and quality were measured using flow cytometry density measurement based on side scatter (SSC-A) and forward scatter (FSC-A). In addition, flow cytometry was used to measure the integrity of protoplasts depend- ing on their forward scatter (FSC) and side scatter (SSC) gating [36], and a 3.8 × 105/µL cell count was recorded (Figure 2C). 3.3. Expression of MeRPPL1 in CRISPR/Cas9-Treated Protoplasts In this study, RT-qPCR was performed at 24 h post-transfection (hpt). To measure the relative expression levels of CRISPR-targeted MeRPPL1 in TME3 protoplasts, real-time qPCR was performed on cDNA using MeRPPL1 primers, with the 18S rRNA gene for normalisation. MeRPPL1 expression in non-transfected leaf protoplasts was low (log2fold −9.8). RT-qPCR comparisons in expression between transfection with CRISPR vectors only (RC), CRISPR vectors and SACMV (RCS), and SACMV only (S) were undertaken (Figure 3 top). The highest log2 fold change (3.1 log2 fold increase) of MeRPPL1 expression normalised to 18S rRNA was recorded in protoplasts that were only transfected with SACMV infectious clones (S), suggesting that SACMV triggers a MeRPPL1 response in SACMV-tolerant TME3-derived protoplasts. A decrease (log2 fold 1.5) in MeRPPL1 expression was observed in protoplasts co- transfected with both the CRISPR/Cas9 vectors and SACMV infectious clones (RCS) com- pared to protoplasts transfected only with SACMV (S) (Figure 3 top). While the expression of MeRPPL1 was not abolished, its expression was significantly down-regulated by CRISPR. These results showed that SACMV may play a role in reducing/mitigating the impact of the CRISPR vectors on MeRPPL1 expression. This is further supported by the fact that the transfection of protoplasts with CRISPR vectors (RC) only, without SACMV, experienced the lowest (log2-fold = 0.9) MeRPPL1 gene expression knockdown compared to SACMV- transfected (S) protoplasts and CRISPR vectors + SACMV (RCS). Alternatively, it may well be that SACMV’s ability to upregulate the expression of MeRPPL1 remains even when CRISPR is acting to reduce it. 3.4. South African Cassava Mosaic Virus Replication Quantification Relative log2 fold change in SACMV replication (viral load) was measured by qPCR of the coat protein gene (AC1) on DNA A at 24 h post-transfection using the log 2(∆∆CT) method (Figure 3 below). Protoplasts that were transfected with only the pBIN19-SACMV- DNA A and pBIN19-SACMV-DNA B infectious clones (S) had significantly higher (2.5-fold) SACMV DNA A accumulation in comparison to co-inoculation with the pC1380-pC1380- Cas9 vector (RC) + SACMV infectious clones (S) (RCS), where transfected protoplasts had a significantly (p ≤ 0.05) lower (1.6-fold) SACMV viral load. Viruses 2023, 15, x FOR PEER REVIEW 8 of 17 the impact of the CRISPR vectors on MeRPPL1 expression. This is further supported by the fact that the transfection of protoplasts with CRISPR vectors (RC) only, without SACMV, experienced the lowest (log2-fold = 0.9) MeRPPL1 gene expression knockdown compared to SACMV-transfected (S) protoplasts and CRISPR vectors + SACMV (RCS). Alternatively, it may well be that SACMV’s ability to upregulate the expression of MeRPPL1 remains even when CRISPR is acting to reduce it. 3.4. South African Cassava Mosaic Virus Replication Quantification Relative log2 fold change in SACMV replication (viral load) was measured by qPCR of the coat protein gene (AC1) on DNA A at 24 h post-transfection using the log 2(ΔΔCT) method (Figure 3 below). Protoplasts that were transfected with only the pBIN19-SACMV-DNA A and pBIN19-SACMV-DNA B infectious clones (S) had signifi- cantly higher (2.5-fold) SACMV DNA A accumulation in comparison to co-inoculation with the pC1380-pC1380-Cas9 vector (RC) + SACMV infectious clones (S) (RCS), where transfected protoplasts had a significantly (p ≤ 0.05) lower (1.6-fold) SACMV viral load. Figure 3. Above (top): log 2 fold change in expression of MeRPPL1 normalised to 18S rRNA gene in cassava TME3 protoplasts following transfection with the CRISPR vectors. Expression was rela- tive to mock which was the pCAMBIA1380 vector only. Below (bottom): Quantification of SACMV measured as the log2 fold change in the CP gene (on DNA A) following CRISPR knockdown of MeRPPL1 expression. S = TME3 protoplasts transfected with SACMV only (pBIN19-SACMV-DNA A and pBIN19-SACMV-DNA B infectious clones); RC = CRISPR vector pC1380-MeRPPL1 + pC1380-Cas9 only; RCS = CRISPR vector + SACMV. Error bars represent means ± S.D. (n = 3). One A 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 RCS RC SR M eR PP L1 Lo g2 Fo ld ch an ge *** * * 0.0 1.0 2.0 3.0 RCS S CP Lo g2 Fo ld C ha ng e * * Figure 3. Cont. Viruses 2024, 16, 941 8 of 15 Viruses 2023, 15, x FOR PEER REVIEW 8 of 17 the impact of the CRISPR vectors on MeRPPL1 expression. This is further supported by the fact that the transfection of protoplasts with CRISPR vectors (RC) only, without SACMV, experienced the lowest (log2-fold = 0.9) MeRPPL1 gene expression knockdown compared to SACMV-transfected (S) protoplasts and CRISPR vectors + SACMV (RCS). Alternatively, it may well be that SACMV’s ability to upregulate the expression of MeRPPL1 remains even when CRISPR is acting to reduce it. 3.4. South African Cassava Mosaic Virus Replication Quantification Relative log2 fold change in SACMV replication (viral load) was measured by qPCR of the coat protein gene (AC1) on DNA A at 24 h post-transfection using the log 2(ΔΔCT) method (Figure 3 below). Protoplasts that were transfected with only the pBIN19-SACMV-DNA A and pBIN19-SACMV-DNA B infectious clones (S) had signifi- cantly higher (2.5-fold) SACMV DNA A accumulation in comparison to co-inoculation with the pC1380-pC1380-Cas9 vector (RC) + SACMV infectious clones (S) (RCS), where transfected protoplasts had a significantly (p ≤ 0.05) lower (1.6-fold) SACMV viral load. Figure 3. Above (top): log 2 fold change in expression of MeRPPL1 normalised to 18S rRNA gene in cassava TME3 protoplasts following transfection with the CRISPR vectors. Expression was rela- tive to mock which was the pCAMBIA1380 vector only. Below (bottom): Quantification of SACMV measured as the log2 fold change in the CP gene (on DNA A) following CRISPR knockdown of MeRPPL1 expression. S = TME3 protoplasts transfected with SACMV only (pBIN19-SACMV-DNA A and pBIN19-SACMV-DNA B infectious clones); RC = CRISPR vector pC1380-MeRPPL1 + pC1380-Cas9 only; RCS = CRISPR vector + SACMV. Error bars represent means ± S.D. (n = 3). One A 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 RCS RC SR M eR PP L1 Lo g2 F ol d ch an ge *** * * 0.0 1.0 2.0 3.0 RCS S CP Lo g2 F ol d Ch an ge * * Figure 3. Above (top): log 2 fold change in expression of MeRPPL1 normalised to 18S rRNA gene in cassava TME3 protoplasts following transfection with the CRISPR vectors. Expression was relative to mock which was the pCAMBIA1380 vector only. Below (bottom): Quantification of SACMV measured as the log2 fold change in the CP gene (on DNA A) following CRISPR knockdown of MeRPPL1 expression. S = TME3 protoplasts transfected with SACMV only (pBIN19-SACMV-DNA A and pBIN19-SACMV-DNA B infectious clones); RC = CRISPR vector pC1380-MeRPPL1 + pC1380- Cas9 only; RCS = CRISPR vector + SACMV. Error bars represent means ± S.D. (n = 3). One and three asterisks indicate statistically significant differences between treated samples at p ≤ 0.05 and 0.001, respectively, using Student’s t-test. 3.5. CRISPR/Cas9-Mediated Mutations in MeRPPL1 DpnI-treated DNA extracted from transfected protoplasts was sent to Inqaba Biotech (Pretoria, South Africa) for sequencing of the MeRPPL1 gene in order to determine nu- cleotide mutations induced by CRISPR, and consequently, the derived translated amino acid (aa) sequence alterations. This is critical, as altered amino acids would lead to changes in the folding/secondary structure of the MeRPPL1 protein, affecting its interaction with SACMV effectors. The sequencing results of ten amplicons from a polyclonal mix revealed different nucleic acid mutations caused by the CRISPR/Cas9 vectors. Substitution muta- tions occurred randomly at a 90–100% frequency, while addition mutations upstream of MeRPPL1 crRNA2 in pC1380-MeRPPL1 + pC1380-Cas9 + SACMV (TME3 RCS) samples occurred at a 10% frequency (Figure 4A). These nucleotide mutations resulted in 40–100% frequencies where amino acid sequences were altered. Interestingly, mutations in both TME3 transfected with the CRISPR silencing vector and SACMV (TME3 RCS) and SACMV only (TME3 S) led to amino acid changes in the MHD domain (from position 177) of the MeRPPL1 peptide, where methionine (M); histidine (H); and aspartic acid (D) amino acids were substituted with isoleucine (I); tyrosine (Y) or methionine (M); and serine (S), cysteine (C), or isoleucine (I), respectively. MeRPPL1 translated peptide from only the CRISPR vector-transfected protoplasts (TME3 RC) retained the conserved MHD domain (Figure 4B). From the examination of the entire MeRPPL1 sequenced region, it was evident that substi- tution mutations occurred randomly throughout the gene region leading to several changes in amino acid sequences of the peptide when compared to the peptide sequence of the wild-type reference genome AM560-2 and untreated TME3 MeRPPL1 homologs (Figure 4C). Interestingly, the infection of TME3 protoplasts with only SACMV also resulted in random substitution mutations that altered genes and translated peptide sequences. Viruses 2024, 16, 941 9 of 15Viruses 2023, 15, x FOR PEER REVIEW 10 of 17 Figure 4. (A) A summary of the types of mutations observed in MeRPPL1 in the SACMV-tolerant TME3 cultivar, and their frequency (derived from sequencing 10 amplicons per genotype), induced by gRNA1 or gRNA2 in RC (silencing vector) or RSC (silencing vector + SACMV)-treated proto- plasts. Random substitution mutations resulting in amino acid changes occurred the most fre- quently (80–100% following gRNA1 treatment and 40–65% with gRNA2), while deletion mutations (yellow bar below graph) were not observed. Addition mutations (10% frequency) were only ob- served in TME3 RCS protoplasts treated with gRNA2. (B) Nucleotide mutations resulting in amino acid mutations in the conserved MHD domain (position 165–180) induced by CRISPR vector (RC), RC + SACMV infectious clones-edited MeRPPL1 (TME3 RCS), and SACMV-infected only (TME3 S). Mutations include substitutions and deletions. (C) Amino acid (aa) sequence alignment of MeRPPL1 in cassava TME3 cultivar protoplasts subjected to the different treatments. Amino acid deletions and substitutions were most frequently observed. Mutations in MeRPPL1 were deter- mined by aligning the sequences of treated samples to MeRPPL1 homologs from wild-type refer- ence genome AM560-2 and TME3 (GenBank ID GCA_003957995.1) using the ClustalW algorithm in MEGA X software (version 2.0). Since alterations in amino acids affect the secondary structure of proteins, and hence their function, amino acid sequences (from Figure 4C) were used to predict secondary protein structures of MeRPPL1 following the different treatments, with reference to wild-type AM560-2 and TME3 MeRPPL1 homologs, using AlphaFold.ipynb (Figure 5; Figures S4–S7). CRISPR-induced mutations in amino acids affected the secondary struc- ture of the MeRPPL1 predicted protein to varying levels, resulting in model prediction scores between very low (pLDDT < 50) and very high (pLDDT > 90) [37]. MeRPPL1 from TME3 RC displayed a secondary structure more similar to that of the AM560-2/TME3 reference protein compared to TME3 RCS and TME3 S. The protein secondary structure was mostly made up of residues that were modelled to high (pLDDT between 70 and 90) and very high (pLDDT > 90) confidence levels, making it a good backbone prediction that 165 177 207 90 149 149 207 266 266 325 325 385 A C B 0 20 40 60 80 100 TME3 RC gRNA1 TME3 RC gRNA2 TME3 RCS gRNA 1 TME3 RCS gRNA 2 Fr eq ue nc y (% ) Substitution mutation Deletion mutation Addition mutation Altered amino acid sequence Figure 4. (A) A summary of the types of mutations observed in MeRPPL1 in the SACMV-tolerant TME3 cultivar, and their frequency (derived from sequencing 10 amplicons per genotype), induced by gRNA1 or gRNA2 in RC (silencing vector) or RSC (silencing vector + SACMV)-treated proto- plasts. Random substitution mutations resulting in amino acid changes occurred the most frequently (80–100% following gRNA1 treatment and 40–65% with gRNA2), while deletion mutations (yellow bar below graph) were not observed. Addition mutations (10% frequency) were only observed in TME3 RCS protoplasts treated with gRNA2. (B) Nucleotide mutations resulting in amino acid muta- tions in the conserved MHD domain (position 165–180) induced by CRISPR vector (RC), RC + SACMV infectious clones-edited MeRPPL1 (TME3 RCS), and SACMV-infected only (TME3 S). Mutations include substitutions and deletions. (C) Amino acid (aa) sequence alignment of MeRPPL1 in cassava TME3 cultivar protoplasts subjected to the different treatments. Amino acid deletions and substi- tutions were most frequently observed. Mutations in MeRPPL1 were determined by aligning the sequences of treated samples to MeRPPL1 homologs from wild-type reference genome AM560-2 and TME3 (GenBank ID GCA_003957995.1) using the ClustalW algorithm in MEGA X software (version 2.0). Since alterations in amino acids affect the secondary structure of proteins, and hence their function, amino acid sequences (from Figure 4C) were used to predict secondary pro- tein structures of MeRPPL1 following the different treatments, with reference to wild-type AM560-2 and TME3 MeRPPL1 homologs, using AlphaFold.ipynb (Figure 5; Figures S4–S7). CRISPR-induced mutations in amino acids affected the secondary structure of the MeRPPL1 predicted protein to varying levels, resulting in model prediction scores between very low (pLDDT < 50) and very high (pLDDT > 90) [37]. MeRPPL1 from TME3 RC displayed a sec- ondary structure more similar to that of the AM560-2/TME3 reference protein compared to TME3 RCS and TME3 S. The protein secondary structure was mostly made up of residues that were modelled to high (pLDDT between 70 and 90) and very high (pLDDT > 90) confidence levels, making it a good backbone prediction that also characterised binding Viruses 2024, 16, 941 10 of 15 sites. TME3 RCS had a few residues with good backbone prediction (pLDDT between 70 and 90). The remainder of the protein structure, and that of the TME3 S model, were dominated by regions with low and very low confidence, with pLDDT scores between 50 and 70 and less than 50, respectively. Low-confidence regions still appeared to have a backbone prediction; however, this should be interpreted cautiously. Very low confidence scores also often predict a ribbon-like structure, which can be a strong indication of protein structure disorder. This suggests that silencing CRISPR vectors + SACMV and SACMV caused mutations leading to highly disordered proteins. These mutations could poten- tially affect the functioning of the protein and its binding capacity to SACMV virulence factors/proteins and/or other cellular interacting proteins. Viruses 2023, 15, x FOR PEER REVIEW 11 of 17 also characterised binding sites. TME3 RCS had a few residues with good backbone pre- diction (pLDDT between 70 and 90). The remainder of the protein structure, and that of the TME3 S model, were dominated by regions with low and very low confidence, with pLDDT scores between 50 and 70 and less than 50, respectively. Low-confidence regions still appeared to have a backbone prediction; however, this should be interpreted cau- tiously. Very low confidence scores also often predict a ribbon-like structure, which can be a strong indication of protein structure disorder. This suggests that silencing CRISPR vectors + SACMV and SACMV caused mutations leading to highly disordered proteins. These mutations could potentially affect the functioning of the protein and its binding capacity to SACMV virulence factors/proteins and/or other cellular interacting proteins. Figure 5. MeRPPL1 protein secondary structure predictions derived from deduced amino acids. The structures were predicted using translated amino acid sequences from the Sanger sequencing results and alignment in Figure 4C. The reference sequence is the MeRPPL1 homolog from wild-type AM560-2 and TME3 (GenBank ID GCA_003957995.1). TME3 RC = TME3 protoplasts transfected with CRISPR vectors only. TME3 RCS = TME3 protoplasts co-transfected with CRISPR vectors + SACMV infectious clones. TME3 S = TME3 protoplasts transfected with SACMV infec- tious clones. pLDDT corresponds to the model’s predicted score and is used to colour-code the residues of the predicted model in the 3D structure viewer [37]. pLDDT > 90 are expected to be modelled to high accuracy; between 70 and 90 are expected to be modelled well; between 50 and 70 are low confidence; less than 50 is a reasonably strong predictor of disorder. 4. Discussion While many studies on resistance proteins (NLRs) have been undertaken in exper- imental plant hosts such as Arabidopsis thaliana and Nicotiana benthamiana, and in a range of annual plants, there are few studies on NLRs in perennial crops, including cassava. Additionally, most studies have been performed on pathogen-resistant plants, not vi- Figure 5. MeRPPL1 protein secondary structure predictions derived from deduced amino acids. The structures were predicted using translated amino acid sequences from the Sanger sequencing results and alignment in Figure 4C. The reference sequence is the MeRPPL1 homolog from wild-type AM560-2 and TME3 (GenBank ID GCA_003957995.1). TME3 RC = TME3 protoplasts transfected with CRISPR vectors only. TME3 RCS = TME3 protoplasts co-transfected with CRISPR vectors + SACMV infectious clones. TME3 S = TME3 protoplasts transfected with SACMV infectious clones. pLDDT corresponds to the model’s predicted score and is used to colour-code the residues of the predicted model in the 3D structure viewer [37]. pLDDT > 90 are expected to be modelled to high accuracy; between 70 and 90 are expected to be modelled well; between 50 and 70 are low confidence; less than 50 is a reasonably strong predictor of disorder. 4. Discussion While many studies on resistance proteins (NLRs) have been undertaken in experi- mental plant hosts such as Arabidopsis thaliana and Nicotiana benthamiana, and in a range of annual plants, there are few studies on NLRs in perennial crops, including cassava. Additionally, most studies have been performed on pathogen-resistant plants, not virus- tolerant hosts. Despite some recent progress [38–41], the molecular mechanisms of cassava Viruses 2024, 16, 941 11 of 15 mosaic disease (CMD) resistance or tolerance are not fully comprehended. It is reported from transcriptome/proteome data that the upregulation of R gene expression [24,28] and the overexpression of some R proteins [25] are induced upon the invasion of cassava by SACMV; however, no laboratory data have been provided to date to support these findings. This study demonstrates experimentally, for the first time in the orphan crop cassava, the involvement of a CNL-type R gene, MeRPPL1, in response to a geminivirus associated with CMD tolerance. The CRISPR/Cas9 vector successfully induced mutations in MeRPPL1 in the SACMV-tolerant TME3 cultivar, resulting in a significant decrease in the expression of the MeRPPL1 gene compared to TME3 protoplasts transfected with SACMV only. This re- sulted in an increase in SACMV-DNA A replication, providing evidence that this CNL gene is involved in the response to SACMV. Notably, the highest log2 fold increase in MeRPPL1 expression relative to 18S rRNA was recorded in protoplasts that were transfected with only SACMV infectious clones (S), suggesting that SACMV triggers a significant MeRPPL1 response in SACMV-tolerant TME3-derived protoplasts. The fact that MeRPPL1 maps to chromosome 12 associated with the CMD2 resistance-associated QTL (Figure S3) made it an informed candidate for this study. A cassava protoplast system was selected, as virus DNA accumulation in cassava leaf mesophyll and tobacco BY-2 protoplasts has previously been observed from 6 to 24 h [30] and 36 to 48 h [42] post-transfection, respectively. Additionally, protoplasts lack a cell wall, making it easy for the viral DNA to enter the cell by transfusion across the cell membrane and replicate inside the host cells [43]. Furthermore, in a previous study, the stability of transient expression in TME3 protoplasts of a MeE3 ligase was verified by the detection of eGFP expression at 24 hpt and also showed that at least 90% of protoplasts had been successfully transformed [30]. The expression of MeRPPL1 in this study was low in untreated cassava TME3 protoplasts, and in tomato and potato, approximately 10% of all R genes were shown to be differentially expressed (either up or down) in non-infected plants [44]. Differences in R gene expression can also be tissue-specific. In this study, young leaves just below the apical meristem of cassava TME3 were selected for protoplast generation, as they represent actively dividing cells/growth that geminiviruses target due to their replication cycle. Since NLRs are well known for their functional role as receptors for viral pathogen effectors [45,46], which induce downstream effector-triggered immunity (ETI), we propose that the knockdown of the expression of the CNL, MeRPPL1, negatively affected ETI, resulting in an increase in SACMV replication in cassava protoplasts. A more recent study [47] also identified a novel LRR-RLK protein that contributes to tolerance to tobacco mosaic virus (TMV) in N. tabacum. Interestingly, the level of this protein was significantly elevated in TMV-infected tolerant N. tabacum alongside an E3 ubiquitin ligase (E3L). It was also shown in SACMV-infected cassava TME3 protoplasts that an E3L was upregulated and is also involved in tolerance to SACMV [30]. The lack of studies using the CRISPR/Cas system to target plant defence response genes in protoplasts limits further comparisons. Future in planta studies will be performed to confirm the observations in this study. The NLR gene, Ty-2 in tomato, is among the few identified genes shown to confer resistance to begomoviruses in plants [18]. Other reported geminivirus-associated NLRs are CRY1 (CC-NBS-LRR) and Sw-5a (SD-CC-NBS-LRR) in Vigna radiate and Solanum lycoper- sicum, respectively [16]. While Ty-2 in tomatoes plays a role in resistance to tomato yellow leaf curl virus (TYLCV), the type of gene has not been demonstrated [16]. MeRRPL1 is the first NLR gene reported in cassava to respond to a CMD-associated geminivirus. However, the Avr protein/effector associated with CMD1 or CMD2-associated resistance in cassava has not been identified. It is suggested that SACMV Rep or CP encoded by the AC1 or AV1 gene, respectively, on DNA A could be candidate effectors for NLR-mediated resistance; however, further studies are warranted. Further research on the effect of MeRPPL1 expres- sion knockdown on the quantification of protein translation also needs to be determined from in planta studies. Furthermore, molecular mechanisms/responses in plants, such as resistance, tolerance and biotic stress, often display similar molecular signatures/features, Viruses 2024, 16, 941 12 of 15 and in perennial crops, such as cassava, it has been proposed that CMD tolerance may also be associated with biotic stress responses [48]. Nucleotide mutations induced by CRISPR/Cas9 vector and SACMV, singly or to- gether, led to mutations in the MeRPPL1 protein in treated TME3 protoplasts, resulting in about 40–90% deduced amino acid alterations in MeRPPL1 peptides (Figure 4). Such muta- tions may have affected the detection of SACMV DNA-derived effectors by the MeRPPL1 protein or reduced its functionality, possibly leading to an attenuated or delayed ETI defence response. This would have led to increased SACMV-DNA A replication, which was observed in this study. Similar results in different hosts have been reported, where geminiviruses co-infected with CRISPR vectors induced random mutations in targeted genes/proteins, resulting in a loss of gene function. Some examples include studies on TYLCV in tomato [49,50], SACMV in cassava [30], and cotton leaf curl geminiviruses in cotton [51]. A previous study in TME3 protoplasts also revealed that SACMV DNA accumulation induced multiple mutations in a MeE3L homolog [30]. Resistance genes have been shown to be highly susceptible to mutations caused by viral and bacterial pathogens [49–52]. For example, Potato virus X induced random mutations in twenty-five CC-NLR R genes in potato [52]. It was also demonstrated in this study that the disruption of the central nucleotide-binding (NBS) domain’s highly conserved and strictly ordered MHD motif in the Rx R gene occurred. Notably, mutations caused by SACMV also resulted in several amino acid alterations in the MHD motif of the MeRPPL1 protein in CRISPR vectors + SACMV and SACMV-only treated protoplasts (Figure 4B). The MHD motif is found in the ARC2 region of the NB-ARC domain. It functions in nucleotide-dependent conformational changes, which lead to the exchange of ADP for ATP and the subsequent activation of the defence mechanisms [53]. Mutations in this motif have been linked to the auto-activation of many NLR proteins in the absence of a viral effector [54–57]. The histidine (H) residue was noted to be crucial to the functioning of this motif in transducing pathogen perception by LRR [58,59]. In this study, histidine was substituted with tyrosine (Y) in RCS treatments, and cysteine (C) in S-treated TME3 protoplasts (Figure 4B), sug- gesting that these substitutions could potentially affect the functioning of this MeRPPL1 protein as a receptor, hence disrupting further downstream defence responses. This is not the first report of the disruption of the MHD motif in R proteins in cassava TME3 landrace. It was previously reported that SACMV-induced amino acid substitutions in planta resulted in the disruption or total removal of the MHD motif in some R gene ac- cessions [28]. Further, mutations in R proteins could lead to alterations in the predicted secondary structure, which would affect the conformation of the protein and how the protein interacts with other cellular proteins and molecules. Herein, mutations were shown to affect the secondary structure of the MeRPPL1 predicted proteins in TME3, as illustrated in Figure 5. From these results, it is proposed that in TME3 protoplasts, the mutations induced by CRISPR and SACMV in the MHD motif and other regions of the MeRPPL1 gene cause secondary structure alterations of MeRPPL1 proteins which would have disrupted binding of SACMV effector proteins or binding and interactions with other cellular proteins involved in the activation of signal pathways for defence. Further studies on the effects of protein alterations would prove informative. 5. Conclusions The CRISPR/Cas9 system successfully induced the gene expression knockdown of MeRPPL1 in TME3 protoplasts. CRISPR-mediated introduction of random point mutations or deletions resulted in polypeptide amino acid alterations in the MeRPPL1 protein. We demonstrate that these two events contribute to an increase in SACMV DNA replication. Therefore, we conclude that this CC-NLR (CNL) gene, MeRPPL1, plays a role in cassava TME3 tolerance to SACMV. Amino acid alterations in the highly conserved MHD motif in MeRPPL1 proteins may have hindered the recognition and binding of the SACMV effector due to alterations in the predicted secondary structure of the MeRPPL1 protein. Thus, the activation of downstream signalling and defence pathways would be disrupted, leading Viruses 2024, 16, 941 13 of 15 to high SACMV-DNA A replication in TME3 protoplasts. Further identification of the SACMV effector, and in planta studies, would prove informative. Investigation into the intersection between innate immunity (PTI), R protein triggering of ETI, and biotic stress is warranted and would reveal some clues for the future management and engineering of resistance to geminiviruses in crops. Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/v16060941/s1, Figure S1: CRISPR-Cas9 gRNA cassette for MeRPPL1 knockdown; Figure S2: pC1380-Cas9 vector; Figure S3: Map of MeRPPL1 (Manes.12G091600) lo- cus (red text) on chromosome 12 of the Manihot esculenta v8.1 genome; Figure S4: Predicted per- residue confidence metric (pLDDT) and Predicted Aligned Error of MeRPPL1 protein from references TME3/AM560-2; Figure S5: Predicted per-residue confidence metric (pLDDT) and Predicted Aligned Error of MeRPPL1 protein from TME3 RC sample; Figure S6: Predicted per-residue confidence metric (pLDDT) and Predicted Aligned Error of MeRPPL1 protein from TME3 RCS sample; Figure S7: Pre- dicted per-residue confidence metric (pLDDT) and Predicted Aligned Error of MeRPPL1 protein from TME3 S sample. Author Contributions: Conceptualisation, E.R. and C.R.; methodology, E.R.; formal analysis, E.R.; investigation, E.R.; resources, C.R.; data curation, E.R.; writing—original draft preparation, E.R.; writing—review and editing, E.R. and C.R.; supervision, C.R.; project administration, E.R. and C.R.; funding acquisition, C.R. All authors have read and agreed to the published version of the manuscript. Funding: This research and a bursary was supported by C.R. through the National Research Founda- tion South Africa (nrf.ac.za), grant number 104829 awarded in 2017. E.R. was also supported as a PhD student through the Agricultural Research Council (https://www.arc.za). Data Availability Statement: All relevant data are available in the manuscript. Conflicts of Interest: The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. References 1. Lebot, V. Tropical Root and Tuber Crops: Cassava, Sweet Potato, Yams, Aroids; CABI Publishing: Wallingford, UK, 2008; ISBN 9781845934248. 2. Zerbini, F.M.; Briddon, R.W.; Idris, A.; Martin, D.P.; Moriones, E.; Navas-Castillo, J.; Rivera-Bustamante, R.; Roumagnac, P.; Varsani, A.; ICTV Report Consortium. ICTV Virus Taxonomy Profile: Geminiviridae. J. Gen. Virol. 2017, 98, 131–133. [CrossRef] 3. Rey, C.; Vanderschuren, H. Cassava Mosaic and Brown Streak Diseases: Current Perspectives and Beyond. Annu. Rev. Virol. 2017, 4, 429–452. [CrossRef] 4. Owor, B.; Legg, J.P.; Okao-Okuja, G.; Obonyo, R.; Ogenga-Latigo, M.W. The Effect of Cassava Mosaic Geminiviruses on Symptom Severity, Growth and Root Yield of a Cassava Mosaic Virus Disease-Susceptible Cultivar in Uganda. Ann. Appl. Biol. 2004, 145, 331–337. [CrossRef] 5. Legg, J.P.; Lava Kumar, P.; Makeshkumar, T.; Tripathi, L.; Ferguson, M.; Kanju, E.; Ntawuruhunga, P.; Cuellar, W. Cassava Virus Diseases: Biology, Epidemiology, and Management. Adv. Virus Res. 2015, 91, 85–142. 6. Berrie, L.C.; Palmer, K.; Rybicki, E.P.; Hidayat, S.H.; Maxwell, D.P.; Rey, M.E.C. A New Isolate of African Cassava Mosaic Virus in South Africa. Afr. J. Root Tuber Crop. 1997, 2, 49–52. 7. Alabi, O.J.; Kumar, P.L.; Naidu, R.A. Cassava Mosaic Disease: A curse to food security in sub-Saharan Africa. APSnet Features 2011, 7, 1–17. 8. Fondong, V.N.; Rey, C. Recent Biotechnological Advances in the Improvement of Cassava. In Cassava; Waisundara, V., Ed.; InTech: London, UK, 2018. [CrossRef] 9. Akano, A.O.; Dixon, A.G.O.; Mba, C.; Barrera, E.; Fregene, M. Genetic Mapping of a Dominant Gene Conferring Resistance to Cassava Mosaic Disease. Theor. Appl. Genet. 2002, 105, 521–525. [CrossRef] 10. Rabbi, I.Y.; Hamblin, M.T.; Kumar, P.L.; Gedil, M.A.; Ikpan, A.S.; Jannink, J.L.; Kulakow, P.A. High-Resolution Mapping of Resistance to Cassava Mosaic Geminiviruses in Cassava Using Genotyping-by-Sequencing and Its Implications for Breeding. Virus Res. 2014, 186, 87–96. [CrossRef] 11. Okogbenin, E.; Egesi, C.N.; Olasanmi, B.; Ogundapo, O.; Kahya, S.; Hurtado, P.; Marin, J.; Akinbo, O.; Mba, C.; Gomez, H.; et al. Molecular Marker Analysis and Validation of Resistance to Cassava Mosaic Disease in Elite Cassava Genotypes in Nigeria. Crop Sci. 2012, 52, 2576–2586. [CrossRef] https://www.mdpi.com/article/10.3390/v16060941/s1 https://www.mdpi.com/article/10.3390/v16060941/s1 https://www.arc.za https://doi.org/10.1099/JGV.0.000738 https://doi.org/10.1146/annurev-virology-101416-041913 https://doi.org/10.1111/j.1744-7348.2004.tb00390.x https://doi.org/10.5772/intechopen.70758 https://doi.org/10.1007/s00122-002-0891-7 https://doi.org/10.1016/j.virusres.2013.12.028 https://doi.org/10.2135/cropsci2011.11.0586 Viruses 2024, 16, 941 14 of 15 12. Wolfe, M.D.; Rabbi, I.Y.; Egesi, C.; Hamblin, M.; Kawuki, R.; Kulakow, P.; Lozano, R.; Del Carpio, D.P.; Ramu, P.; Jannink, J. Genome-Wide Association and Prediction Reveals Genetic Architecture of Cassava Mosaic Disease Resistance and Prospects for Rapid Genetic Improvement. Plant Genome 2016, 9, plantgenome2015.11.0118. [CrossRef] 13. Wu, X.; Valli, A.; García, J.A.; Zhou, X.; Cheng, X. The Tug-of-War between Plants and Viruses: Great Progress and Many Remaining Questions. Viruses 2019, 11, 203. [CrossRef] 14. Ding, L.N.; Li, Y.T.; Wu, Y.Z.; Li, T.; Geng, R.; Cao, J.; Zhang, W.; Tan, X.L. Plant Disease Resistance-Related Signaling Pathways: Recent Progress and Future Prospects. Int. J. Mol. Sci. 2022, 23, 16200. [CrossRef] [PubMed] 15. Ngou, B.P.M.; Ahn, H.-K.; Ding, P.; Jones, J.D.G. Mutual Potentiation of Plant Immunity by Cell-Surface and Intracellular Receptors. Nature 2021, 592, 110–115. [CrossRef] 16. Sett, S.; Prasad, A.; Prasad, M. Resistance Genes on the Verge of Plant–Virus Interaction. Trends Plant Sci. 2022, 27, 1242–1252. [CrossRef] 17. Maiti, S.; Paul, S.; Pal, A. Isolation, Characterization, and Structure Analysis of a Non-TIR-NBS-LRR Encoding Candidate Gene from MYMIV-Resistant Vigna mungo. Mol. Biotechnol. 2012, 52, 217–233. [CrossRef] [PubMed] 18. Shen, X.; Yan, Z.; Wang, X.; Wang, Y.; Arens, M.; Du, Y.; Visser, R.G.F.; Kormelink, R.; Bai, Y.; Wolters, A.M.A. The NLR Protein Encoded by the Resistance Gene Ty-2 Is Triggered by the Replication-Associated Protein Rep/C1 of Tomato Yellow Leaf Curl Virus. Front. Plant Sci. 2020, 11, 545306. [CrossRef] [PubMed] 19. Yamaguchi, H.; Ohnishi, J.; Saito, A.; Ohyama, A.; Nunome, T.; Miyatake, K.; Fukuoka, H. An NB-LRR Gene, TYNBS1, Is Responsible for Resistance Mediated by the Ty-2 Begomovirus Resistance Locus of Tomato. Theor. Appl. Genet. 2018, 131, 1345–1362. [CrossRef] 20. Sharma, N.; Sahu, P.; Prasad, A.; Muthamilarasan, M.; Waseem, M.; Khan, Y.; Thakur, J.K.; Chakraborty, S.; Prasad, S. The Sw5a gene confers resistance to ToLCNDV and triggers an HR response after direct AC4 effector recognition. Proc. Natl. Acad. Sci. USA 2021, 118, e2101833118. [CrossRef] 21. Gill, U.; Scott, J.W.; Shekasteband, R.; Ogundiwin, E.; Schuit, C.; Francis, D.M.; Sim, S.-C.; Smith, H.; Hutton, S.F. Ty-6, a Major Begomovirus Resistance Gene on Chromosome 10, Is Effective against Tomato Yellow Leaf Curl Virus and Tomato Mottle Virus. Theor. Appl. Genet. 2019, 132, 1543–1554. [CrossRef] 22. Lim, Y.W.; Mansfeld, B.N.; Schläpfer, P.; Gilbert, K.B.; Narayanan, N.N.; Qi, W.; Wang, Q.; Zhong, Z.; Boyher, A.; Gehan, J.; et al. Mutations in DNA polymerase δ subunit 1 co-segregate with CMD2-type resistance to Cassava Mosaic Geminiviruses. Nat. Commun. 2022, 13, 3933. [CrossRef] 23. Gedil, M.; Kumar, M.; Igwe, D. Isolation and Characterization of Resistant Gene Analogs in Cassava, Wild Manihot Species, and Castor Bean (Ricinus communis). Afr. J. Biotechnol. 2012, 11, 15111–15123. [CrossRef] 24. Allie, F.; Pierce, E.J.; Okoniewski, M.J.; Rey, C. Transcriptional Analysis of South African Cassava Mosaic Virus-Infected Susceptible and Tolerant Landraces of Cassava Highlights Differences in Resistance, Basal Defense and Cell Wall Associated Genes during Infection. BMC Genom. 2014, 15, 1006. [CrossRef] [PubMed] 25. Ramulifho, E.; Rey, M.E.C. Proteome Mapping of South African Cassava Mosaic Virus-Infected Susceptible and Tolerant Landraces of Cassava. Proteomes 2021, 9, 41. [CrossRef] [PubMed] 26. Prochnik, S.; Marri, P.R.; Desany, B.; Rabinowicz, P.D.; Kodira, C.; Mohiuddin, M.; Rodriguez, F.; Fauquet, C.; Tohme, J.; Harkins, T.; et al. The Cassava Genome: Current Progress, Future Directions. Trop. Plant Biol. 2012, 5, 88–94. [CrossRef] [PubMed] 27. Kuon, J.E.; Qi, W.; Schläpfer, P.; Hirsch-Hoffmann, M.; Von Bieberstein, P.R.; Patrignani, A.; Poveda, L.; Grob, S.; Keller, M.; Shimizu-Inatsugi, R.; et al. Haplotype-Resolved Genomes of Geminivirus-Resistant and Geminivirus-Susceptible African Cassava Cultivars. BMC Biol. 2019, 17, 75. [CrossRef] [PubMed] 28. Louis, B.; Rey, C. Resistance Gene Analogs Involved in Tolerant Cassava–Geminivirus Interaction That Shows a Recovery Phenotype. Virus Genes 2015, 51, 393–407. [CrossRef] [PubMed] 29. Concordet, J.P.; Haeussler, M. CRISPOR: Intuitive Guide Selection for CRISPR/Cas9 Genome Editing Experiments and Screens. Nucleic Acids Res. 2018, 46, W242–W245. [CrossRef] [PubMed] 30. Chatukuta, P.; Rey, M.E.C. A Cassava Protoplast System for Screening Genes Associated with the Response to South African Cassava Mosaic Virus. Virol. J. 2020, 17, 184. [CrossRef] [PubMed] 31. Murashige, T.; Skoog, F. A Revised Medium for Rapid Growth and Bio Assays with Tobacco Tissue Cultures. Physiol. Plant. 1962, 15, 473–497. [CrossRef] 32. Berrie, L.C.; Rybicki, E.P.; Rey, M.E.C. Complete Nucleotide Sequence and Host Range of South African Cassava Mosaic Virus: Further Evidence for Recombination amongst Begomoviruses. J. Gen. Virol. 2001, 82, 53–58. [CrossRef] 33. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2-∆∆CT Method. Methods 2001, 25, 402–408. [CrossRef] [PubMed] 34. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X Version 2.0. Bioinformatics 2007, 23, 2947–2948. [CrossRef] [PubMed] 35. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Žídek, A.; Potapenko, A.; et al. Highly Accurate Protein Structure Prediction with AlphaFold. Nature 2021, 596, 583–589. [CrossRef] [PubMed] 36. Haynes, J.L. Principles of Flow Cytometry. Cytometry 1988, 9, 7–17. [CrossRef] [PubMed] 37. Mariani, V.; Biasini, M.; Barbato, A.; Schwede, T. LDDT: A Local Superposition-Free Score for Comparing Protein Structures and Models Using Distance Difference Tests. Bioinformatics 2013, 29, 2722–2728. [CrossRef] https://doi.org/10.3835/plantgenome2015.11.0118 https://doi.org/10.3390/v11030203 https://doi.org/10.3390/IJMS232416200 https://www.ncbi.nlm.nih.gov/pubmed/36555841 https://doi.org/10.1038/s41586-021-03315-7 https://doi.org/10.1016/j.tplants.2022.07.003 https://doi.org/10.1007/s12033-011-9488-1 https://www.ncbi.nlm.nih.gov/pubmed/22205242 https://doi.org/10.3389/fpls.2020.545306 https://www.ncbi.nlm.nih.gov/pubmed/33013967 https://doi.org/10.1007/s00122-018-3082-x https://doi.org/10.1073/pnas.2101833118 https://doi.org/10.1007/s00122-019-03298-0 https://doi.org/10.1038/s41467-022-31414-0 https://doi.org/10.5897/AJB12.2373 https://doi.org/10.1186/1471-2164-15-1006 https://www.ncbi.nlm.nih.gov/pubmed/25412561 https://doi.org/10.3390/proteomes9040041 https://www.ncbi.nlm.nih.gov/pubmed/34842800 https://doi.org/10.1007/s12042-011-9088-z https://www.ncbi.nlm.nih.gov/pubmed/22523606 https://doi.org/10.1186/s12915-019-0697-6 https://www.ncbi.nlm.nih.gov/pubmed/31533702 https://doi.org/10.1007/s11262-015-1246-1 https://www.ncbi.nlm.nih.gov/pubmed/26370397 https://doi.org/10.1093/nar/gky354 https://www.ncbi.nlm.nih.gov/pubmed/29762716 https://doi.org/10.1186/s12985-020-01453-4 https://www.ncbi.nlm.nih.gov/pubmed/33228712 https://doi.org/10.1111/j.1399-3054.1962.tb08052.x https://doi.org/10.1099/0022-1317-82-1-53 https://doi.org/10.1006/meth.2001.1262 https://www.ncbi.nlm.nih.gov/pubmed/11846609 https://doi.org/10.1093/BIOINFORMATICS/BTM404 https://www.ncbi.nlm.nih.gov/pubmed/17846036 https://doi.org/10.1038/s41586-021-03819-2 https://www.ncbi.nlm.nih.gov/pubmed/34265844 https://doi.org/10.1002/CYTO.990090804 https://www.ncbi.nlm.nih.gov/pubmed/3076370 https://doi.org/10.1093/BIOINFORMATICS/BTT473 Viruses 2024, 16, 941 15 of 15 38. Houngue, J.A.; Zandjanakou-Tachin, M.; Ngalle, H.B.; Pita, J.S.; Cacaï, G.H.T.; Ngatat, S.E.; Bell, J.M.; Ahanhanzo, C. Evaluation of Resistance to Cassava Mosaic Disease in Selected African Cassava Cultivars Using Combined Molecular and Greenhouse Grafting Tools. Physiol. Mol. Plant Pathol. 2019, 105, 47–53. [CrossRef] [PubMed] 39. Thuy, C.T.; Lopez-Lavalle, L.A.; Vu, N.A.; Hy, N.H.; Nhan, P.T.; Ceballos, H.; Newby, J.; Tung, N.B.; Hien, N.T.; Tuan, L.N.; et al. Identifying New Resistance to Cassava Mosaic Disease and Validating Markers for the CMD2 Locus. Agriculture 2021, 11, 829. [CrossRef] 40. Sheat, S.; Winter, S. Developing Broad-Spectrum Resistance in Cassava against Viruses Causing the Cassava Mosaic and the Cassava Brown Streak Diseases. Front. Plant Sci. 2023, 14, 1042701. [CrossRef] [PubMed] 41. Codjia, E.D.; Olasanmi, B.; Agre, P.A.; Uwugiaren, R.; Ige, A.D.; Rabbi, I.Y. Selection for Resistance to Cassava Mosaic Disease in African Cassava Germplasm Using Single Nucleotide Polymorphism Markers. S. Afr. J. Sci. 2022, 118, 1–7. [CrossRef] 42. Vanitharani, R.; Chellappan, P.; Fauquet, C.M. Short Interfering RNA-Mediated Interference of Gene Expression and Viral DNA Accumulation in Cultured Plant Cells. Proc. Natl. Acad. Sci. USA 2003, 100, 9632–9636. [CrossRef] 43. Houston, K.; Tucker, M.R.; Chowdhury, J.; Shirley, N.; Little, A. The Plant Cell Wall: A Complex and Dynamic Structure As Revealed by the Responses of Genes under Stress Conditions. Front. Plant Sci. 2016, 7, 984. [CrossRef] 44. von Dahlen, J.K.; Schulz, K.; Nicolai, J.; Rose, L.E. Global expression patterns of R-genes in tomato and potato. Front. Plant Sci. 2023, 14, 1216795. [CrossRef] 45. Gouveia, B.C.; Calil, I.P.; Machado, J.P.B.; Santos, A.A.; Fontes, E.P.B. Immune Receptors and Co-Receptors in Antiviral Innate Immunity in Plants. Front. Microbiol. 2017, 7, 2139. [CrossRef] 46. Yuan, M.; Ngou, B.P.M.; Ding, P.; Xin, X.F. PTI-ETI Crosstalk: An Integrative View of Plant Immunity. Curr. Opin. Plant Biol. 2021, 62, 102030. [CrossRef] 47. Wang, J.; Hao, F.; Song, K.; Jin, W.; Fu, B.; Wei, Y.; Shi, Y.; Guo, H.; Liu, W. Identification of a Novel NtLRR-RLK and Biological Pathways That Contribute to Tolerance of TMV in Nicotiana tabacum. Mol. Plant-Microbe Interact. 2020, 33, 996–1006. [CrossRef] 48. Freeborough, W.; Gentle, N.; Rey, M.E.C. WRKY Transcription Factors in Cassava Contribute to Regulation of Tolerance and Susceptibility to Cassava Mosaic Disease through Stress Responses. Viruses 2021, 13, 1820. [CrossRef] 49. Tashkandi, M.; Ali, Z.; Aljedaani, F.; Shami, A.; Mahfouz, M.M. Engineering Resistance against Tomato Yellow Leaf Curl Virus via the CRISPR/Cas9 System in Tomato. Plant Signal. Behav. 2018, 13, e1525996. [CrossRef] 50. Pramanik, D.; Shelake, R.M.; Park, J.; Kim, M.J.; Hwang, I.; Park, Y.; Kim, J.Y. CRISPR/Cas9-Mediated Generation of Pathogen- Resistant Tomato against Tomato Yellow Leaf Curl Virus and Powdery Mildew. Int. J. Mol. Sci. 2021, 22, 1878. [CrossRef] 51. Mubarik, M.S.; Wang, X.; Khan, S.H.; Ahmad, A.; Khan, Z.; Amjid, M.W.; Razzaq, M.K.; Ali, Z.; Azhar, M.T. Engineering Broad-Spectrum Resistance to Cotton Leaf Curl Disease by CRISPR-Cas9 Based Multiplex Editing in Plants. GM Crop. Food 2021, 12, 647–658. [CrossRef] 52. Bendahmane, A.; Farnham, G.; Moffett, P.; Baulcombe, D. Constitutive Gain-of-Function Mutants in a Nucleotide Binding Site-Leucine Rich Repeat Protein Encoded at the Rx Locus of Potato. Plant J. 2002, 32, 195–204. [CrossRef] 53. DeYoung, B.; Innes, R. Plant NBS-LRR Proteins in Pathogen Sensing and Host Defense. Nat. Immunol. 2006, 7, 1243–1249. [CrossRef] 54. Howles, P.; Lawrence, G.; Finnegan, J.; McFadden, H.; Ayliffe, M.; Dodds, P.; Ellis, J. Autoactive Alleles of the Flax L6 Rust Resistance Gene Induce Non-Race-Specific Rust Resistance Associated with the Hypersensitive Response. Mol. Plant-Microbe Interact. 2005, 18, 570–582. [CrossRef] 55. de la Fuente van Bentem, S.; Vossen, J.; de Vries, K.; van Wees, S.; Tameling, W.; Dekker, H.; de Koster, C.; Haring, M.; Takken, F.; Cornelissen, B. Heat Shock Protein 90 and Its Co-Chaperone Protein Phosphatase 5 Interact with Distinct Regions of the Tomato I-2 Disease Resistance Protein. Plant J. 2005, 43, 284–298. [CrossRef] 56. Dodds, P.; Lawrence, G.; Catanzariti, A.; Teh, T.; Wang, C.; Ayliffe, M.; Kobe, B.; Ellis, J. Direct Protein Interaction Underlies Gene-for-Gene Specificity and Coevolution of the Flax Resistance Genes and Flax Rust Avirulence Genes. Proc. Natl. Acad. Sci. USA 2006, 103, 8888–8893. [CrossRef] 57. van Ooijen, G.; Mayr, G.; Kasiem, M.; Albrecht, M.; Cornelissen, B.; Takken, F. Structure-Function Analysis of the NB-ARC Domain of Plant Disease Resistance Proteins. J. Exp. Bot. 2008, 59, 1383–1397. [CrossRef] 58. Tameling, W.I.L.; Elzinga, S.D.J.; Darmin, P.S.; Vossen, J.H.; Takken, F.L.W.; Haring, M.A.; Cornelissen, B.J.C. The Tomato R Gene Products I-2 and Mi-1 Are Functional ATP Binding Proteins with ATPase Activity. Plant Cell 2002, 14, 2929–2939. [CrossRef] 59. McHale, L.; Tan, X.; Koehl, P.; Michelmore, R.W. Plant NBS-LRR Proteins: Adaptable Guards. Genome Biol. 2006, 7, 212. [CrossRef] 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. https://doi.org/10.1016/j.pmpp.2018.07.003 https://www.ncbi.nlm.nih.gov/pubmed/31007373 https://doi.org/10.3390/agriculture11090829 https://doi.org/10.3389/fpls.2023.1042701 https://www.ncbi.nlm.nih.gov/pubmed/36778712 https://doi.org/10.17159/sajs.2022/11607 https://doi.org/10.1073/PNAS.1733874100 https://doi.org/10.3389/FPLS.2016.00984 https://doi.org/10.3389/fpls.2023.1216795 https://doi.org/10.3389/fmicb.2016.02139 https://doi.org/10.1016/J.PBI.2021.102030 https://doi.org/10.1094/MPMI-12-19-0343-R https://doi.org/10.3390/v13091820 https://doi.org/10.1080/15592324.2018.1525996 https://doi.org/10.3390/ijms22041878 https://doi.org/10.1080/21645698.2021.1938488 https://doi.org/10.1046/J.1365-313X.2002.01413.X https://doi.org/10.1038/NI1410 https://doi.org/10.1094/MPMI-18-0570 https://doi.org/10.1111/J.1365-313X.2005.02450.X https://doi.org/10.1073/PNAS.0602577103 https://doi.org/10.1093/JXB/ERN045 https://doi.org/10.1105/tpc.005793 https://doi.org/10.1186/GB-2006-7-4-212 Introduction Materials and Methods MeRPPL1 Expression Knockdown gRNA Cassette Design and Synthesis Cloning the CRISPR Cassettes into pCAMBIA 1380 Delivery Vector CRISPR-Mediated MeRPPL1 Expression Knockdown: Protoplast Isolation and Transfection RNA Extraction and RT-qPCR Gene Expression Analysis DNA Extraction and Viral Replication Quantification Sequencing, Mutagenesis Analysis, and Protein Structure Prediction Results MeRPPL1 gRNA Cassette Confirmation Protoplast Isolation and Viability Expression of MeRPPL1 in CRISPR/Cas9-Treated Protoplasts South African Cassava Mosaic Virus Replication Quantification CRISPR/Cas9-Mediated Mutations in MeRPPL1 Discussion Conclusions References