1. African cassava landraces susceptible and tolerant to SACMV are amenable to enzymatic protoplast isolation and PEG-mediated transformation
Protoplasts were chosen for this transient gene expression study because they can conveniently and efficiently be transformed with several DNA constructs simultaneously, and they allow higher resolution imaging compared to cells in intact tissue (Faracoet al., 2011). Additionally, they can be used for high-throughput efficient screening of candidate genes, and those genes that show an effect can then be silenced by virus induced gene silencing (VIGS) in planta, which takes considerably longer (3-4 months) and requires more complicated procedures for the non-model host cassava. Leaf mesophyll was used as the source of protoplasts because SACMV exerts its effects mainly in the leaves, where symptoms arise.Leaf mesophyll protoplasts therefore provide functional information (Faracoet al., 2011) relating to the effect of SACMV at leaf tissue level. Round and irregularly-shaped protoplasts of different sizes were observed (Fig.1B-D), although spherical leaf mesophyll protoplasts dominated, as generally reported (Wu et al., 2017).
A previously determined enzyme concentration (1.6% cellulase, 0.8% macerozyme) that is suitable for obtaining the optimum number of viable protoplasts was used for leaf cell wall digestion (Wu et al., 2017). The viability of protoplasts in this study was at least 85% although cassava protoplast viability of up to 95% has been reported (International Plant Research Institute, 1984). The long digestion period (16 h) was ruled out as the cause of death for ~15% of protoplasts as this reportedly does not induce serious damage in protoplasts (Tang, 1982). It has been reported that micro-propagated plants grown in vitro lack epicuticular wax and thus allow rapid enzyme penetration (Kumar and Rao, 2018). Cassava, however, has a thick epicuticular layer which necessitates the long digestion period of 16h compared to 0.3-1h for Arabidopsis (Wu et al., 2009). Macerating enzymes such as macerase are known to cause wound reactions in protoplasts because of their degradation of the cell wall, which may lead to necrosis (Ishii, 1988), while cellulase is known to exert inadequate enzymatic activity at low concentrations and higher concentrations have no benefit or detriment (Uchimiya and Murashige, 1974). Therefore, a balance between digestion enzyme concentration and viability is essential in order to obtain the optimum yield of viable good quality protoplasts. Cassava protoplast viability presented herein may differ from previously reported percentages possibly due to differences inthe cassava genotypes used, and in this particular case could be due to the particular physiological characteristics of the African cassava landraces from which protoplasts were derived.
Although lower than the previously reported yields of 4.4x107 protoplasts/g FW leaves from M. esculentacv. South China 8 (Wu et al., 2017) and 1.9x107 protoplasts/g FW leaves from M. esculenta cv. M. Thai 8 (Anthony et al., 1995), the protoplast yields in this present study (4.90-6.36x106/g FW) were sufficient to provide the recommended number of protoplasts (104-107) required for each transfection (Yooet al., 2007). Both pre-treatment of leaves in the dark (24-72 h) (Shahin and Shepard, 1980) and vacuum infiltration before enzyme digestion (Nanjareddy, 2016) have been shown to intensify enzyme penetration in bean leaves. We found that pre-treatment resulted in release of much undesirable plant debris alongside protoplasts, and that vacuum infiltration did indeed help increase protoplast yield. Flow cytometry indicated a high concentration of protoplasts compared to irregularly-shaped debris. Based on microscopy images of purified protoplasts (Fig.1B-D), the irregular debris outside the flow cytometry gated area was deemed to be free chloroplasts, plasmolysed cells, undigested cell wall fragments and other aggregates arising from the long digestion period of leaf material.
Stability of protoplast transformation was confirmed by expression of eGFP from the CRISPR construct at 24 hpt following PEG-mediated transfection of 104 protoplasts with 15 µg CRISPR construct and 4 µg SACMV infectious clones. Protoplasts were deemed unsuitable for further analysis from 36 hpt as they rapidly lost viability. PEG-mediated plant protoplast transfection with plasmid DNA is a well-established procedure (Hayashimotoet al., 1990; Locatelli et al., 2003) and a popular protocol uses 10 μg DNA to transfect 2x104 Arabidopsis protoplasts (Yooet al., 2007). At least 5 µg of plasmid DNA have previously been used to transform 106 tobacco protoplasts with African cassava mosaic virus (ACMV) (Ermaket al., 1993) and Cowpea mosaic comovirus (Wellinket al., 1993). Highly efficient co-expression of multiple constructs in plant protoplasts has been reported (Chen et al., 2006; Walter et al., 2004) and virus infectious clones have also been used in conjunction with other plasmid constructs for co-inoculation of plant protoplasts (Cheng et al., 2009). Nicotianatabacum protoplasts have been co-transformed with 5 µg eGFP construct, 3 μg siRNA and 4 μg each of ACMV or East African cassava mosaic virus (EACMV) with DNA/RNA extraction at 36 and 48 hpt (Vanitharaniet al., 2003). To our knowledge, the present study is the first report of co-transformation of cassava protoplasts with a CRISPR construct and geminivirus infectious clones.The number of constructs used for transformation can be increased for high-throughput studies to target multiple genes simultaneously since it has been shown that transformation efficiency is independent of plasmid amount (Locatelli et al., 2003).
The percentage mutation frequency of 50-80% at targeted sites in the cassava genome is comparable to the 70% and 60% attained in cassava (Wu et al., 2017) andN. tabacum(Lin et al., 2018) protoplasts respectively, and considerably less than the 100% previously obtained in cassava plants (Odipioet al., 2017).CRISPR targets with GC content greater than 50% are known to achieve higher efficiency than those with less than 50% (Ma et al., 2015), and this may explain why gRNA1 target sequence (55% GC content) has considerably higher mutation frequency than gRNA2 target sequence (50% GC content).
2. The M. esculentaT200 MeE3L encodes a truncated RING-less protein due to a nonsense mutation
The frameshift resulting from the insertion mutation in susceptible T200 MeE3L introduces a stop codon upstream of the RING domain, thus encoding a truncated protein in which the C3HC4-type RING finger motif is absent. Essentially, the E3 ligase domain in the T200 MeE3L would not be translated because of this mutation, making its potential protein product non-functional with respect to this E3 ligase activity. We were able to amplify the T200 MeE3L exon from cDNA, showing that this gene is transcribed, but we have no evidence of its translation or lack thereof. It is proposed that the loss of this substantial portion of the T200 MeE3L C-terminal region would not only avert E3 ligase activity but may also alter the spatial chemical conformation necessary for any other interactions (such as binding the ubiquitin-conjugated E2 and substrate) to occur. The protocol in this present study therefore enablesgene sequence comparison among wild type cassava genotypes, to determine the presence of single nucleotide polymorphisms (SNPs) or other forms of mutations, and to form a basis for in planta exploration of phenotypic differences among cassava genotypes.
The AM560-2, TME3, and cv.60444 cassava genotypes have previously been shown to cluster together under their nearest ancestor, Heveabrasiliensis(Bredesonet al., 2016). The evolutionary history of T200, a southern African landrace, is unknown and phylogenetic analysis suggests that the T200 MeE3L evolved after the cv.60444 and TME3 variants. It is known that wild plants in natural ecosystems co-evolve with their virus partners. While it is recognized that there is a relationship between virus virulence/pathogenicity and co-adaptation to plant hosts (Sacristan and Garcia-Arenal, 2008), information regarding how viruses apply selective pressure to alter plant susceptibility is not known. A study of Drosophila and its host-specific viruses found that coevolution may cause sustained genetic variation in susceptibility (Duxbury et al., 2019). This may explain why a southern African cassava landrace is highly susceptible to SACMV that appears to have migrated south from its origin, suspected to be in east Africa or the south-west Indian Ocean islands such as Madagascar (De Bruynet al., 2016) that geographically separated from the African continent (Lefeuvreet al. 2007). South African cassava mosaic virus is a recombinant between East African cassava mosaic virus and two other unknown geminiviruses which contributed the AC4 and IR regions (Berrieet al., 2001), and moved southwards into Mozambique, Zimbabwe and South Africa where it may have encountered the T200 landrace. Subsequent to its first discovery in South Africa, SACMV has been reported in Zimbabwe (Briddonet al. 2004) and Madagascar (De Bruynet al., 2016). It is known that infection with anew recombinant begomovirus requires the host to adjust to minor or major differences in virus-host interactions (Montes et al., 2019).It is suggested that the T200 landrace and SACMV may still be in the process of co-adaption, which would explain why T200 exhibits extreme susceptibility to SACMV.We speculate that the MeE3L is either a paralog in T200 or it was introgressed from a wild relative in southern Africa.
3. SACMV DNA accumulation in cassava protoplasts is genotype-dependent
Quantitative PCR is a well-established method for precise quantitation of viral DNA amount in infected tissue and it requires a host reference gene with stable expression patterns under experimental conditions as the internal control for correct data normalisation (Moreno et al., 2011). Data from the qPCR measurement of SACMV DNA accumulation (relative to the 18S rRNA gene) show that SACMV DNA accumulates in cassava protoplasts, correlating well with previous reports of geminivirus DNA accumulation in planta and in vitro. Quantitative detection of African cassava mosaic virus and East African cassava mosaic virus using qPCR has been reported (Ottiet al., 2013) and SACMV titre, in particular, has been assayed in planta in Arabidopsis (Pierce and Rey, 2013) and cassava (Allieet al., 2014). Replication of the geminivirus, Cassava brown streak virus, in cassava leaf mesophyll protoplasts has been assayed at 6 hpt (Anjanappaet al., 2016) and it has been reported that there was significant viral DNA accumulation in tobacco BY-2 protoplasts 36 and 48 hpt by co-inoculating with infectious ACMV and EACMV clones and siRNA (Vanitharaniet al., 2003). The present study is the first to report accumulation of SACMV DNA in cassava protoplasts.
Based on previously reported in planta evidence (Allie et al., 2014) and the known presence of a CMD2 locus in tolerant TME3 (Akanoet al., 2002),it was expected that SACMV DNA accumulation would be genotype-dependent and significantly lower in TME3 than in the model cv.60444 and susceptible T200 protoplasts. Interestingly, there was differential SACMV accumulation in CRISPR-transformed cassava protoplasts expressing the gene-edited MeE3L. The upregulation of SACMV DNA accumulation in susceptible T200 and tolerant TME3 in the presence of mutant MeE3L suggests a role for MeE3L as one of the host genes involved in the response to SACMV infection. CRISPR-associated modification of MeE3L may enhance SACMV DNA accumulation in susceptible T200 and tolerant TME3 by interfering with the ubiquitin proteasome system-dependent tolerance/resistance response mechanisms of cassava.
4. Viral activity and gene editing of MeE3L affect the expression of MeE3L
Geminiviruses elude plant defense mechanisms by hijacking and redirecting ubiquitination, and interfering with responses regulated by ubiquitin E3 ligases (including responses to jasmonates, auxins, gibberellins, ethylene,abscisic acid) (Lozano-Duran et al., 2011). It follows then that alterations to E3 ligase genomic sequences may alter E3 ligase expression patterns during viral infection, as viruses are known to modulate RNA levels to enhance infection (Verchot, 2016). Both plant viruses and CRISPR systems are known to induce mutations in the genome (Machida et al., 2004; Cougotet al., 2005; Sander and Joung, 2014), and the employment of both against the MeE3L would provide an indication whether MeE3L may be involved in the plant’s response to SACMV.
Previously, plant E3 ligases have been shown to be induced by viral infection (Lai et al., 2009; Czosneket al., 2013; Chen et al., 2018) and plant defence elicitors (Libaultet al., 2007; Sadanandomet al., 2012).It is known that geminiviruses interact with plant E3 ligases and induce their up- or down-regulation to promote infection or undergo degradation (Lai et al., 2009; Lozano-Duran et al., 2011; Lozano-Duran and Bejarano, 2011; Shen et al., 2016). Results presented herein indicate that in TME3 protoplasts, MeE3L expression is upregulated during SACMV infection. The concurrent CRISPR-mediated gene editing of MeE3L and infection with SACMV appears to induce increased expression of the MeE3L, suggesting that MeE3L’s specific base sequence is important for the interaction between the virus and the plant host. The muted response of the T200 MeE3L to all treatments was expected given its nonsense mutation which silences the RING domain responsible for E3 ligase activity. However, the muted response of model cv.60444 MeE3L was unexpected and suggests that this MeE3L sequence variant is not responsive to SACMV infection.
The MeE3L homolog sequences in SACMV-infected protoplasts reflect a silenced RING domain, suggesting that SACMV may possibly induce silencing of the RING domain in order to achieve full infection of the host. The concomitant increase in SACMV DNA accumulation and gene-editedMeE3L in TME3 points to the response of MeE3L to SACMV being more directed at advancing susceptibility. There is evidence for geminiviral (Tomato yellow leaf curl sardinia virus) silencing of a plant E3 ligase, RHF2A, to promote infection (Lozano-Duran et al., 2011) and impairment of plant defenceduring Cabbage leaf curl virus (CaLCuV) infection due toinhibition of a RING E3 ligase (Sahuet al., 2013). The present study provides further evidence that geminiviruses may interfere with activity of a plant E3 ligase in plants.
5. SACMV’s interaction with a tolerant cassava genotype induces numerous mutations in MeE3L
Functions of E3 ligases in regulating immunity systems are orchestrated at the interface of host-virus interactions (Zhang et al., 2018) and some of these interactions occur in the nucleus (Kushwahaet al., 2017). Sequencing of genomic MeE3L from SACMV-infected TME3 protoplasts revealed multiple random single base mutations along the length of MeE3L,which translate to amino acid substitution (Fig. 4E). While these mutations do not alter the reading frame, they are predicted to silence the whole protein and not just the RING domain. The resulting disordered protein would presumably not only lack RING E3 ligase activity, but also the E2 and substrate binding activity. These mutations were present and similar in all 10 genomic DNA amplicons derived from the polyclonal mixes of each of 3 biological replicates. Similar mutations encoding multiple stop codons have been observed in an Argonaute 4-encoding gene (Manes.18g121900)from SACMV-transformed tolerant TME3 protoplasts (unpublished data; Chatukuta and Rey), indicating that other host genes may be similarly affected by SACMV infection.
Interestingly, the discovery of mutations in genomic DNA presented herein possibly point to a yet unknown geminivirus-induced host mechanism for genome editing. Geminiviruses are known to induce the expression of genes related to repair of double-stranded breaks (DSBs) and DNA synthesis (Lozano et al., 2015), and to promote somatic homologous recombination (Richter et al., 2014).Some E3 ligases and viral proteins can localise to the nucleus, such as the tobacco E3 ligase, NtHUB1 which has a nuclear localisation sequence, is recruited by geminiviral Rep protein, and co-localises and interacts with the Rep protein to monoubiquitinate cellular chromatin and thus enable infection (Kushwahaet al., 2017). The viral coat protein, CP, also has a nuclear localisation signal, can localise in the nucleolus and nucleoplasm, and facilitates entry of ssDNA into the nucleus (Wang et al., 2017; Kumar, 2019). However, the mechanisms for SACMV-mediated gene mutation induction in cassava protoplasts are yet to be investigated.
6. The response of MeE3L to SACMV is virus- and host-specific
Ubiquitin ligases are abundant in plants and provide substrate specificity to target particular proteins. In Arabidopsis alone, RING E3 ligases make up 499 out of over 1,500 E3 ligases (Mazzucotelliet al., 2006). A comparison of E3 ligase and E3 ligase complex-associated gene expression during other plant geminivirus infection studies (Additional File 3) was conducted to determine whether MeE3L’s response to SACMV is geminivirus-specific or host-dependent.
In susceptible cassava, E3 ligase expression is downregulated during SACMV infection at early, middle and late time points (12, 32 and 64 days post infection (dpi)) but there is no differential expression of E3 ligases in tolerant cassava at any time point (Allie et al., 2014). However, no differential expression of E3 ligases is recorded during SACMV infection of Arabidopsis which is susceptible (Pierce and Rey, 2013). A study of transcriptomic responses to geminivirusTomato leaf curl New Delhi virus (ToLCNDV) infection in potato found that five E3 ligases in the susceptible cultivar and two in the tolerant cultivar are upregulated at 30 dpi (Jeevalathaet al., 2017) while the geminivirusTomato yellow leaf curl virus (TYLCSV) has been shown to induce upregulation of E3 ligases in susceptible tomato at 42 dpi, except in the case of a CUL1 which is downregulated (Miozziet al., 2014). A transcriptome study of Arabidopsis during geminivirusCaLCuV infection found that, out of 1570 E3 ligases, 149 were up-regulated and 23 were downregulated (Ascencio-Ibanez et al., 2008; Lozano-Duran et al., 2013). The CaLCuV AC2 protein, in particular, induces downregulation of two E3 ligases in Arabidopsis (Liu et al., 2014). These findings, together with the current study, prove that plant E3 ligase responses to geminivirus infection are neither uniform nor similar, but they vary according to the specific geminivirus and host involved in the interaction.
Responses of cassava to the ssRNApotyvirusesCassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) with respect to E3 ligase expression variably show both downregulation and upregulation in the susceptible varieties. Interestingly, there is no differential expression of E3 ligases in resistant cassava varieties except in Kaleso where a CUL1 is upregulated and a RZPF34 is downregulated (Maruthiet al, 2014; Anjanappaet al., 2016; Amugeet al., 2017). This variable expression of E3 ligases with respect to the virus in the same host suggests that while responses to viral infection are host-dependent, they are also modulated according to the particular virus infecting the plant.
SACMV infection in planta is associated with occurrence of severe symptoms leading to persistent severe infection in susceptible T200 and mild symptoms with recovery at 67 days post infection (dpi) in tolerant TME3 (Allie et al., 2014). The SACMV-induced genetic mutations and differential expression of MeE3L post-infection in TME3 and T200 indicate that it is one of the significant genes involved in the plant’s response to the virus.In planta proteome data from our laboratory shows that during SACMV infection, an E3 ligase (Manes.08G075100) is upregulated in susceptible T200 and downregulated in tolerant TME3 cassava plants at 32 and 67 dpi (unpublished data; Rey), supporting indications from the protoplast system that E3 ligases are responsive to SACMV infection.
7. Limitations
This protocol presented herein suffers some limitations due to the independent cell nature of protoplasts and the short-lived viability of cassava protoplasts in particular. It cannot be used to study cell wall-related genes, cell-to-cell signalling, intercellular movement, long-term responses, or long-term stability of CRISPR-induced gene edits. Further, the use of mesophyll protoplasts may not be suitable to correlate responses in other organs such as roots or flowers. The widely used T7 endonuclease I (T7EI) assay for detecting gene editing activity produced inconclusive results for this present study, and therefore gene editing was indicated by restriction digestion and confirmed by sequencing. It has been reported that CRISPR-Cas9 activity is more accurately reflected by Next Generation Sequencing (NGS) (Sentmanatet al., 2018).In planta validation of results from this protoplast protocol, particularly overexpression and virus-induced gene silencing (VIGS)of targeted genes, as well as functional and interaction studies, must be conducted to confirm the specific roles played by candidate genes in the host-virus interaction.