Presence of the integrated endogenous banana streak virus (eBSV) in the B genome of plantain (AAB) is a major challenge for breeding and dissemination of hybrids. As the eBSV activates into infectious viral particles under stress, the progenitor Musa balbisiana and its derivants, having at least one B genome, cannot be used as parents for crop improvement. Here, we report a strategy to inactivate the eBSV by editing the virus sequences. The regenerated genome-edited events of Gonja Manjaya showed mutations in the targeted sites with the potential to prevent proper transcription or/and translational into functional viral proteins. Seventy-five percent of the edited events remained asymptomatic in comparison to the non-edited control plants under water stress conditions, confirming inactivation of eBSV into infectious viral particles. This study paves the way for the improvement of B genome germplasm and its use in breeding programs to produce hybrids that can be globally disseminated.
Banana streak virus (BSV) is a plant pathogenic badnavirus of the family Caulimoviridae, affecting production of banana (Musa spp.). It is an unenveloped, non-covalently closed, bacilliform double-stranded DNA (dsDNA) virus with a monopartite genome of ~7.2–7.8 kb encoding three open-reading frames (ORFs). ORF1 encodes a small protein of unknown function that is associated with virions1. OFR2 encodes a protein of ~14 kDa, which is involved in virion assembly due to the presence of predicted N-terminal coiled-coil domain1. ORF3 contains domains associated with movement, the virus capsid, aspartic protease, reverse transcriptase, and ribonuclease H functions2.
BSV was first identified in West Africa in 1958 and is now reported in most banana and plantain growing countries3. Like other badnaviruses, BSV sequences integrate into the genome of Musa spp.4. BSV is a complex of different viruses belonging to the pararetroviruses and classified as endogenous pararetroviruses when they are integrated into a host genome. The BSV integrated in the banana host genome are known as endogenous BSV (eBSV).
Banana cultivars are polyploid clones derived from Musa accuminata (A genome) or/and Musa balbisiana (B genome). The eBSV sequences are mainly integrated in the B genome derived from M. balbisiana5. Many economically important sub-groups of banana, such as plantain (AAB), an important staple food in Africa, contain at least one B genome. Multiple copies of eBSV viral sequences have become integrated as direct and inverted tandem repeats at a single locus in the B genome of the host during viral infection5. When the banana plants are stressed, the eBSV recombines to produce a functional episomal viral genome and infectious viral particles and as a result the plant develops disease symptoms. In vitro culture for propagation and hybridization through conventional breeding may also trigger its activation. Consequently, BSV is considered as a major constraint in banana breeding programs, restricting the use of the diploid progenitor M. balbisiana or its derivants carrying a B genome as parents for introgression of desirable agronomic traits6. It also restricts germplasm movement of genotypes with the B genome worldwide due to this potential activation of eBSV into the episomal infectious form of virus. Therefore, it is crucial to design a strategy to irreversibly silence the latent eBSV in the B genome. In this study, CRISPR/Cas9-based genome-editing technology was applied to inactivate eBSV sequences in the host plantain genome. The genome-edited events of plantain Gonja Manjaya were generated with mutations in the targeted sites of integrated eBSV sequences in the host genome. Seventy-five percent (6 out of 8) of the edited events tested remained asymptomatic in comparison to the control non-edited plants under water stress conditions, confirming inhibition of eBSV and reversal of its ability to be converted into infectious viral particles in the edited lines.
Confirmation of presence of eBSV in genome of Gonja Manjaya
Gonja Manjaya, a false horn plantain (AAB genome) commonly grown in East and Central Africa, was initiated from an asymptomatic field plant and used in this study to knockout the integrated sequences of eBSV. The BSV strain Obino l’Ewai (BSOLV), is widely distributed among plantain plantations in Africa. BSOLV exhibits in two forms, episomal (Fig. 1a, GenBank accession KJ013506) and eBSOLV5 (Fig. 1b, GenBank accession AP009334) integrated in a recombined form at a single locus in the B genome of plantain. The presence of BSOLV or eBSOLV sequences was confirmed in the mother stock of in vitro plantlets by PCR analysis using the three different primer pairs (P1–P3) (Fig. 1c), which amplify both the endogenous and episomal forms of the virus. The presence of the episomal form of BSOLV was ruled out by PCR analysis using primer pair P4 that specifically amplifies BSOLV but not the integrated eBSOLV7. The P4 primer pair (Fig. 1a) binds in the ORF3 region of BSOLV and generates a 1336 bp fragment in the presence of the replicating viral genome. This primer pair also binds to eBSOLV, but the position and orientation of the two primers’ binding precludes PCR amplification (Fig. 1b). The amplified fragment diagnostic for the episomal BSOLV was also observed in the symptomatic plants of plantain Agbagba, collected from the field as positive controls (PC1–2) and a symptomatic glasshouse plant of Gonja Manjaya (PC3). However, no amplification was observed in any of the in vitro Gonja Manjaya plantlets, similar to the asymptomatic control plant of Cavendish Williams (AAA) (Fig. 1d).
The presence of eBSOLV in the host genome was further confirmed by Southern blot analysis using a BSOLV-specific probe (termed SBP, Fig. 1b–e). The genomic DNA was restricted with HindIII, which has only one site in the episomal BSOLV and three sites in the eBSOLV (Fig. 1a, b). As expected based on the genomic sequence of eBSOLV integrated in the reference genome of M. balbisiana (Accession AP009334) (Fig. 1b), two hybridization signals (about 13.5 and 5.4 kb) were observed in all four samples tested (Fig. 1e). As predicted, a band of 13.5 kb was obtained by restriction of two HindIII sites within eBSOLV DNA sequence and a second fragment of about 5.4 kb was obtained by cleavage at one restriction site in eBSOLV DNA sequence and a second site in the flanking genome of Gonja Manjaya bordering the integrated site.
gThe in vitro plantlet 3 with the expected two bands on Southern hybridization (Fig. 1e) and PCR product characteristics of integrated eBSOLV, but no episomal BSOLV (Fig. 1d) was used for developing embryogenic cell suspensions. These cell suspensions were used to generate genome-edited plants after confirming the absence of episomal BSOLV and the presence of integrated eBSOLV in their genome by PCR analysis.
Preparation of CRISPR/Cas9 plasmid construct
The gRNAs were designed for each CRISPR target from the most conserved sequences of the BSOLV (GenBank accession KJ013506) and eBSOLV (GenBank accession AP009334 and HE983609). The gRNAs [gRNA1 (ORF1)-TTGAGCAAAGAAGACGTTG, gRNA2 (ORF2)-GCTGGAACAACTGGTGACT and gRNA3 (ORF3)-ATGGAGTTCATATGATCAT] were designed for the three CRISPR targets (S1–S3) in the three ORFs of BSOLV and eBSOLV (Fig. 1a, b). All the three gRNAs were integrated into the T-DNA region of a Cas9 expression vector pDe-Cas9-Kan8 to generate the plasmid, pMR215-BSV1 (Fig. 2a).
Generation and validation of transgenic events
The CRISPR/Cas9 construct was delivered into the embryogenic cell suspensions of Gonja Manjaya via Agrobacterium tumefaciens. A total of 100 kanamycin-resistant, well-rooted plantlets were generated (Fig. 2b–d). Fifty random regenerants were selected and screened for the presence of the Cas9 gene by PCR using primers P5 (Fig. 2a, Supplementary Table 1). The transgenic events displayed a 906 bp fragment, as expected, confirming the presence of the Cas9 gene (Supplementary Figure 1a).
The well-rooted plantlets of Cas9 PCR-positive events and non-transgenic control were acclimatized and potted in the greenhouse (Fig. 2e–f).
Molecular analysis to confirm editing in ORF1 and ORF2
Three gRNAs (targeting sequences S1 in ORF1, S2 in ORF2, and S3 in ORF3) were introduced into triploid Gonja Manjaya (AAB). The eBSOLV is known to integrate at a single locus on chromosome 1 of the B genome based on the M. balbisiana sequence (GenBank accession AP009334 and HE983609). There is only one infectious allele of eBSOLV6. BSOLV sequences are duplicated in the same or opposite direction resulting in several copies of ORFs (either complete or partial) within the integrated eBSOLV (Fig. 1b). As a result of these duplications, each guide (S1–S3) has three targets in the BSOLV sequences integrated in the plant genome (Fig. 1b).
Twenty independent transgenic events (1, 4, 7, 17, 29, 30, 33, 36, 48, 53, 66, 72, 76, 78, 81, 85, 87, 94, 97, and 100), confirmed for the presence of the Cas9 gene, were analyzed for a shift in the size of PCR amplicons as an indication for mutations in target sites S1 and S2. PCR analysis was performed using primer pair P1 flanking the gRNA1 (in S1 targets) and gRNA2 (in S2 targets) (Fig. 1, Supplementary Figure 1b). The amplicons of several of these independently derived events showed a predicted band shift of 197 bp, suggesting deletions in the region between S1 and S2 as expected when Cas9 cuts both target sites simultaneously (Supplementary Figure 1b). Seven (36, 66, 81, 85, 87, 94, and 97) out of 20 events showed shorter amplicons (band shift by 197 bp) in addition to amplicons similar in size to those of wild type control plants suggesting that they might be heterozygous for mutations in all the copies of target sites S1 and S2. Two events (81 and 85) showed a third amplicon with a smaller deletion indicating these plants contained additional genome rearrangements. Two other events (48 and 76) had only the smaller amplicons suggesting all target sites might have the same mutation, while the remaining nine events (1, 4, 7, 17, 29, 53, 72, 78, and 100) showed amplicons similar to the wild type control plant (Supplementary Figure 1b), indicating that they might have no indels or small indels.
A T7 endonuclease I (T7E1) assay was performed to check for the presence of small indels in the nine transgenic events (1, 4, 7, 17, 29, 53, 72, 78, and 100), which did not show any shift in the size of amplicon compared with the wild type control. T7E1-digested fragments were detected in all the events except for events 7 and 78, confirming the mutated sites in these events (Supplementary Figure 1c).
All the 20 independent events (1, 4, 7, 17, 29, 30, 33, 36, 48, 53, 66, 72, 76, 78, 81, 85, 87, 94, 97, and 100) tested by PCR analysis were further analyzed for targeted mutations by sequencing of target sites S1 and S2. The PCR products amplified using P1 primers were purified and cloned, and 10 individual clones per event sequenced. Sequencing of the cloned PCR products further confirmed the targeted mutations in all the edited events except for 3 events (7, 33, and 78) confirming a mutation efficiency of 85% in the targeted site S1 and S2 of transgenic plants (Fig. 3). All 17 of these independently derived mutant plants showed either only deletions (∆1 bp to ∆198 bp) or deletions (∆1 bp to ∆17 bp) accompanied by insertion (+1) at the predicted Cas9 cleavage site, 3 bp upstream of PAM motif (Fig. 3). Ten (36, 48, 66, 72, 76, 81, 85, 87, 97, and 100) out of 20 events tested showed deletions of 198 bp, suggesting that Cas9 cuts 4 bp upstream of the PAM for gRNA1 and 3 bp upstream of the PAM for gRNA2. Simultaneous cutting at S1 and S2 site was apparently repaired to generate the deletion of 198 nucleotides in a coding region, eliminating 66 codons. In other events (1, 4, 17, 29, 30, and 53), the deletions were smaller. Events 17 and 94 showed small deletions at both target sites, indicating Cas9 cleavage and independent mutagenic repair events at two nearby target sites in the same eBSOLV. Events 29, 30, and 53 also showed single nucleotide insertion in addition to deletions at one target site. Only one nucleotide replacement was observed in event 1 (Fig. 3). The majority of edited events (15/17) showed mutations at both S1 and S2 targeted sites. All the edited plants contained indels resulting in mutations with the ability to prevent proper transcription of eBSOLV genes or/and translational to non-functional viral proteins.
Only two edited events (48 and 76) showed the same indels in all the 10 clones sequenced, which was in agreement with the result of band shift PCR analysis showing one amplicon (Fig. 3). The other events showed 2‒3 different versions of sequences of the target sites indicating that these events might be chimeric, with mutations happening at different stages of plant regeneration, or indicating differences in mutations at the various copies of the target sites in a single eBSOLV. The sequencing result of events 48 and 76 also confirmed knockout of all the three copies of target sequences of eBSOLV from the plantain genome.
Event 76 was further tested for uniformity and stability of mutations. Sequencing of three different leaves from the same plantlet showed the same deletion (∆198 bp) indicating uniform mutations in event 76. The sequencing of daughter plants of event 76 after many sub-culturing cycles also showed the same deletions (∆198 bp), thus, confirming the stability of mutations during multiple rounds of vegetative propagation.
All four different wild type non-transgenic control plants, generated from the same embryogenic cell line used in the experiments along with the mutant plants, showed wild type sequences with no indels confirming that all the mutations among the different transgenic events were due to the action of CRISPR/Cas9.
Molecular analysis to confirm edits in ORF3 of eBSOLV
The PCR product for all 20 events (1, 4, 7, 17, 29, 30, 33, 36, 48, 53, 66, 72, 76, 78, 81, 85, 87, 94, 97, and 100) in which the target S3 site was amplified using the P3 primer pair were purified and cloned. Ten individual clones per event were analyzed for mutations at the S3 target site by DNA sequencing (Fig. 4). The sequence analysis showed a variety of mutations (deletion or insertion or deletion accompanied by insertion) in 14 events (4, 7, 17, 29, 30, 33, 36, 48, 76, 81, 85, 87, 94, and 100) while the remaining six events (1, 53, 66, 72, 78, and 97) did not show any indels in the S3 site (Fig. 4) indicating editing efficiency of 70% at site S3. For this target site, only smaller deletions ranging between 1 and 7 bp and a single bp insertion were observed. Most of the mutations were observed at 3 bp upstream of PAM, as expected. All the 14 edited events showed indels leading to reading frame shift and/or other mutations that have potential to inactivate the eBSOLV genes (Fig. 4).
Out of the 20 events analyzed for editing in three different target sequences, 11 events (17, 29, 30, 36, 48, 76, 81, 85, 87, 94, and 100) showed editing in all the targeted sites (S1–S3), two events (7 and 33) showed editing only in S3 site with no indels in S1 and S2 site, and four events (53, 66, 72, and 97) showed mutations in S1/S2 sites but no indels in S3 site. Event 1 showed indels at only the S1 site and no indels in S2 and S3 sites, whereas event 4 showed mutations in S1 and S3 sites but no indels in S2 site. There was only one event (78) that did not show indels in any of the targeted sites indicating a very high mutation efficiency of 95% using three gRNAs.
In searches for potential off-target mutations for each gRNA using the CRISPR web-tool Breaking-Cas (http://bioinfogp.cnb.csic.es/tools/breakingcas), a total of 66 potential off-targets were discovered. Because plantain contains two genomes, A and B, the gRNA + PAM sequences were blasted against the respective banana genomes (http://banana-genome-hub.southgreen.fr). Putative 58 off-target sites were identified in the A-genome from M. accuminata and 40 in the B-genome from M. balbisiana (Table 1). Based on the banana genome, all these potential off-target sites had matches of 63‒95% to the gRNA+PAM sequences (Supplementary Table 2). Seven putative off-target loci with high sequence similarity to the gRNAs or/and containing a Cas9 PAM sequence (NGG), were selected for sequencing and detailed examination of potential off-target effects in nine edited events (1, 17, 48, 66, 76, 81, 85, 87, and 97) (Table 2). Of the seven selected genomic loci, only one of them contained PAM motif (NGG), which could potentially be targeted by gRNA–Cas9.
Sequencing of seven potential off-target sites in nine events showed point mutation at one putative off-target locus (Ma01_t10610.1) in the genome of just one event (72). This off-target site was in the A genome sharing 86% nucleotide identity with the target site. The site did not contain the same PAM sequence (GGG) as the intended target, but the guide might potentially be able to recognize an NAG PAM sequence in close proximity. These results indicate a very low probability for potential off-target mutations.
Glasshouse evaluation of genome-edited plants
Eight genome-edited events (17, 48, 66, 76, 81, 85, 87, and 97) along with three wild type control plants of Gonja Manjaya were potted, acclimatized, and grown in the greenhouse. The genome-edited plants and wild type control plants, developed from the same embryogenic cell line, were found to be phenotypically similar with no growth abnormalities (Fig. 2e, f).
After 2 months, the plants were subjected to water stress for 2 weeks and observed daily for development of BSV symptoms for 60 days. Disease symptoms such as broken or continuous streaks of yellow, chlorotic, black, or brown color on the leaf appeared in all three of the wild type control plants tested, as would be expected if there was activation of eBSOLV into episomal infectious BSOLV (Fig. 5a). However, six (17, 48, 76, 81, 85, and 87) out of eight genome-edited events tested remained asymptomic and the remaining two events (66 and 97) showed only moderate symptoms (Fig. 5). The presence of the infectious episomal form of BSOLV was confirmed by PCR using primers7 (employed in Fig. 1d) that detect BSOLV but not eBSOLV. PCR analysis showed amplification of BSOLV in symptomatic plants of genome-edited plants 66 and 97 and in the stressed wild type Gonja Manjaya plant, as well as the symptomatic field plant sample of plantain Agbagba (Fig. 5b). However, BSOLV was not detected in asymptomatic genome-edited (17, 48, 76, 81, 85, and 87) events (Fig. 5b). This result was further validated by qPCR using the same primers7. In these analyses, the cycle threshold (CT) values for all the samples tested ranged from 23 to 35, whereas the stressed wild type Gonja Manjaya plant showed a lower CT of 23.6 (Fig. 5c). The symptomatic field plant of Agbagba (PC1) and symptomatic-edited plants (66 and 97) showed a CT of about 28 confirming the presence of BSOLV. The asymptomatic-edited plants (17, 48, 76, 81, 85, and 87) showed CT values similar or higher to the asymptomatic plants of Cavendish Williams (CW) confirming the absence of BSOLV (Fig. 5c).
The events 66 and 97 showed mutations in S1 and S2 target sites but no mutations in the S3 target site, however the asymptomatic edited events (17, 48, 76, 81, 85, and 87) showed indels in all the targeted sites suggesting either that knockouts in all the ORFs are important for complete inactivation of the virus, or that mutations in S3 are especially effective. This result indicates that targeted mutagenesis of eBSOLV can lead to its stable inactivation, preventing the emergence of infectious BSV.
Banana and plantains are important staple food crops cultivated across tropical and subtropical countries. Genetic improvement of banana by conventional breeding is key for the development of improved varieties with disease and pest resistance and higher yields. BSV is one of the major challenges in banana breeding and dissemination of plantain hybrids due to the presence of integrated eBSV in the B genome4. These proviruses can be activated into the infectious episomal BSV by various stress conditions. Thus, major epidemics caused by BSV are not due to natural transmission through mealy bugs or through use of infected planting materials, but, rather due to activation of eBSV9 under stress conditions, such as in vitro propagation, hybridization or/and unfavorable conditions like water and temperature stress. As a result, BSV has become a hurdle to overcome in plantain breeding. Breeders cannot use diploid M. balbisiana (BB) in breeding programs, even though it has good attributes like hardiness, strong root system, and tolerance to biotic and abiotic stresses. Therefore, it is necessary to make the B genome of banana free of the activatable sequences of eBSV. With this in mind, in this study we knocked out the integrated eBSV by targeted mutagenesis using CRISPR/Cas9.
CRISPR-mediated genome editing for plant virus resistance is mainly reported for ssDNA viruses like geminiviruses10,11,12 and ssRNA viruses13,14. CRISPR/Cas9 technology can also be applied to either dsDNA viruses or retroviruses that have dsDNA as part of their life cycle15. Currently, this has not been demonstrated for any plant viruses, however, application of CRISPR/Cas9 to interrupt both episomal and integrated dsDNA viruses has been demonstrated in many human viruses16,17,18. HIV is a retrovirus that integrates its viral DNA into the host DNA and has potential to reactivate causing the disease, HIV-AIDS. Recently, inactivation of HIV-1, was achieved by knocking out proviral DNA integrated in latently infected cells by targeting long terminal repeat flanking sequences or multiple regulatory genes within the HIV-1 provirus19,20. Similarly, in this study, we used CRISPR/Cas9 to knockout the integrated dsDNA of BSV from the banana genome using gRNAs targeting multiple repeats of all the three ORFs of BSV (Fig. 1).
Our results showed high mutation efficiency of 95% in the targeted sites (S1–S3) using multiple gRNAs, which is in agreement with a recent study in which 100% editing frequency was reported for banana using a polycistronic gRNA plasmid construct targeting the PDS gene21. In our study, the frequency of mutations differed between the target sites. About 85% of events (17/20) showed indels in the S1 and S2 target sites and 59% of these mutants showed big deletions of 198 bp as a result of simultaneous cleavage of both target sites by Cas9. Simultaneous action of two gRNAs targeting sequences in close proximity has been shown to generate large deletions in other plant species as well22,23,24,25. A comparatively lower editing efficiency of 70% with small indels was observed with a single gRNA-targeting site S3. Similarly, a lower editing efficiency of 59% with small indels was also reported previously in banana using a single gRNA targeting the PDS gene26.
The sequencing of target sites S1 and S2 of edited events 48 and 76 showed the same indels in all the clones sequenced, confirming knockout of all the three copies of ORF1 and ORF2 in the integrated eBSOLV sequence from the plantain genome. The eBSOLV has multiple copies of ORFs inserted at the same locus on chromosome 15 (Fig. 1). All the edited events showed indels with mutations in the targeted sites in ORFs resulting in the potential defects in gene transcription or mRNA translation.
Our results confirmed the uniform and stable inheritance of targeted mutations by CRISPR/Cas9 to clonal progeny. A similar result of stable inheritance of a disease resistance trait was observed in the field-testing of transgenic banana27.
A common limitation of Cas9-mediated genome editing is the possibility of off-target mutation of regions similar to the target site28. In order to minimize the off-target effects in this study, the gRNAs were selected based on their specificity to the target site and minimal potential off targets in the Musa genome. Our results showed very low potential off-target site mutations indicating that CRISPR/Cas9-induced mutagenesis in banana is precise and efficient.
The evaluation of eight edited events under water stress conditions in the glasshouse showed BSV symptoms in only two events, whereas the other six events remained asymptomatic. PCR and qPCR analysis of these plants confirmed the activation of eBSOLV into infectious episomal BSOLV in the two symptomatic-edited events and stressed control wild type plants, whereas no infectious virus was detected in the six asymptomatic-edited events. The qPCR assay is reported to be specific for detection of episomal BSOLV7.
The asymptomatic events have indels in all the three ORFs, whereas symptomatic events had deletions only in ORF1 and ORF2. This result suggests that ORF3, encoding essential proteins, such as movement protein, coat protein, aspartic protease, reverse transcriptase, and ribonuclease H, may play a crucial role in symptom development. ORF3 encodes a polyprotein, which is cleaved post-translationally by aspartic protease29 into functional versions of each of the ORF3-encoded proteins. In this study, the mutations created in the aspartic protease gene should have disrupted the functions of the protease and thus, the functions of other ORF3 genes as well. However, there is need for multiple knockouts in all three ORFs to completely prevent activation of eBSV as indels only in one ORF may result in partial inactivation and might allow viral escape. Although the edited plants with only mutations in the ORF3 region were not tested for activation of eBSOLV in this study, plants with mutations only in ORF3 may also exhibit partial inactivation of the virus as ORF2 encodes for a protein that is involved in virion assembly1. Therefore, it appears that simultaneous generation of mutations in all three ORFs is an important strategy for completely inactivating the provirus and preventing the occurrence of symptoms. Indeed, regarding HIV-1, a retrovirus whose provirus is similar to eBSV, studies have demonstrated the strong suppression of HIV-1 using a multiplex approach19,20,30,31.
BSV is a complicated virus and mechanisms for its integration, activation, and subsequent episomal infection are still largely unknown5. There is only one infectious allele6 of eBSOLV integrated in the chromosome 1 of B genome of Gonja Manjaya. However, three or more copies (complete or partial) of each ORF were present either in the same or opposite orientation within the eBSOLV sequences (Fig. 1b).
It is interesting that a few events showing indels in some clones and wild type sequence in other clones showed a knockout phenotype with no BSV symptom development under stress conditions. Although this triploid genome carries only a single copy of the integrated sequence, that integrated sequence carries three copies of each target, in various orientations. The symptom-free plants may carry mutations copies of the target that are essential for production of a functional episomal virus. Further analysis of these plants may enable us to distinguish between essential and disposable copies of the target site.
In summary, we have successfully applied the CRISPR/Cas9 system to edit the integrated sequence of eBSOLV in plantain cultivar Gonja Manjaya. Our results also demonstrated that multiplexing of CRISPR/Cas9 technology is very efficient in creating precise deletions in banana. Although preliminary results showed stable mutations and silencing of the reactivation of eBSOLV, more robust evaluation of promising edited plants in the field for several generations and with new BSV infections, due to natural transmission through mealy bugs in the field conditions, is required.
This is the first report to our knowledge to demonstrate the knockout of integrated endogenous DNA sequences of the pararetrovirus in a host plant genome. We have demonstrated CRISPR/Cas9-based targeted mutagenesis can permanently inactivate endogenous eBSV, and presents a promising model for the inactivation of other endogenous viral genomes. This study paves the way for editing of banana germplasm with B genome(s) that can be used as one or both parents in the breeding programs. This is a particularly important strategy for both improving plantains and enabling global dissemination of the resulting hybrids with improved B genome.
Plantain cultivar Gonja Manjaya (genome AAB), showing no symptoms of BSV under field conditions, was collected from the field in Uganda and used in this study. The in vitro culture was established from the sucker collected from banana field as described32. The outer layer of leaves and corm tissue were removed from the suckers to obtain an explant of 5 × 10 cm long. The explants were washed with water and surface sterilized by immersing in 70% ethanol for 5 min followed by soaking in 15% sodium hypochlorite solution with few drops of Tween 20 for 10 min. After surface sterilization, the explants were rinsed 3–4 times with sterile distilled water. The explants were further trimmed by removing the outer sheaths to obtain the size of 1 × 1 cm. The explant was divided into two pieces through the meristem and cultured in the proliferation medium [4.4 g L−1 Murashige and Skoog medium including vitamins (Duchefa, M0222), 10 mg L−1 ascorbic acid, 5 mg L−1 6-benzylaminopurine, 30 g L−1 sucrose, 3 g L−1 gelrite, pH 5.8] at 26 ± 2 °C under 16 h photoperiod furnished with fluorescent tube providing light of 94 μmol m−2 s−1. The explants were transferred to fresh medium every 2 weeks. Once the in vitro culture was initiated, it was maintained by sub-culturing the shoots to fresh medium every 4 weeks. The shoots regenerated on proliferation medium were transferred to rooting medium [4.4 g L−1 Murashige and Skoog medium including vitamins (Duchefa, M0222), 10 mg L−1 ascorbic acid, 1 mg L−1 indole-3-butyric acid, 30 g L−1 sucrose, 3 g L−1 gelrite, pH 5.8] for development of well-rooted plantlets.
The in vitro mother stock plantlets were checked for the presence of BSOLV or eBSOLV in the host by PCR using three pairs of primers P1–P3 (Supplementary Table 1, Fig. 1). Genomic DNA was extracted from leaves of four in vitro plantlets of Gonja Manjaya plants using CTAB33. About 500 mg of fresh leaf samples were ground to fine powder using mortar and pestle, and resuspended in 500 μL of extraction buffer (100 mM Tris–HCl, pH 8.0; 20 mM EDTA pH 8.0; 1.4 M NaCl; 2% PVP-10; 0.8% cetyltrimethylammonium bromide, and 0.001% mercaptoethanol). The samples were vortex vigorously, incubated at 65 °C for 30 min and then centrifuged at 15,000 rpm for 10 min. The supernatant was collected and mixed with an equal volume of chloroform:isoamyl alcohol (24:1) and centrifuged at 15,000 rpm for 10 min. This step was repeated. After a second chloroform:isoamyl alcohol extraction, the supernatant was collected and mixed with an equal volume of pre-chilled isopropanol and incubated at room temperature for 10 min for precipitation of DNA. Nucleic acids were pelleted by centrifugation at 15,000 rpm for 10 min. The DNA pellets were washed with 70% ethanol, air-dried and resuspended in 100 μL of TE buffer containing RNase. DNA was reprecipitated with 1 mL of absolute ethanol for 30 min and centrifuged at 12,000 rpm for 10 min. The supernatant was discarded and DNA pellet was air-dried and dissolve in 100 μL nuclease-free water. The concentration and quality of the DNA was checked with nanodrop.
PCR was performed in a 20 µL reaction volume containing 1 µL genomic DNA (20 ng µL−1), 10 µL of HotStarTaq master mix, 1 µL of 10 µM of each primer, and 7 µL nuclease-free water. PCR reaction conditions were: initial denaturation step at 95 °C for 15 min, followed by 34 cycles of denaturation at 94 °C for 30 s, annealing at 57 °C for 30 s, extension at 72 °C for 1 min, and final extension at 72 °C for 10 min. PCR amplicons were resolved on 1% agarose gel.
The presence of episomal BSV was ruled out by PCR amplification using primer set P4 (Supplementary Table 1), which amplifies only episomal BSOLV and not integrated eBSOLV7 (Fig. 1). The symptomatic field grown plants of plantain Agbagba (AAB) showing BSV symptoms and symptomatic Gonja Manjaya plants from the glasshouse were used as positive controls and asymptomatic in vitro plantlet of Cavendish Williams (AAA) was used as a negative control. PCR was performed with 20 ng of genomic DNA with an initial denaturation step at 95 °C for 15 min, followed by 34 cycles of denaturation at 94 °C for 30 s, annealing at 55 °C for 30 s, extension at 72 °C for 1 min, and final extension at 72 °C for 10 min. After amplification, 5 µL of PCR product was resolved on 1% agarose gel stained with gel red.
The integration of eBSOLV was further confirmed by Southern blot hybridization. The genomic DNA of four in vitro plantlets of Gonja Manjaya was restricted with HindIII, resolved on 0.8% agarose gel, and transferred to a nylon membrane. HindIII cuts at a single site within the BSOLV genome and three sites in eBSOLV (Fig. 1a, b). The blot was hybridized with DIG-labeled probe specific to the ORF3 region of BSOLV (Fig. 1) and detected by CDP star chemiluminiscent substrate.
After validation for the presence of eBSOLV and absence of BSOLV, the in vitro plantlet 3 was multiplied and used for developing embryogenic cell suspensions using multiple buds following the method as described32. The embryogenic callus was induced by culturing slices (1–2 mm thick) of multiple buds on callus induction medium [4.4 g L−1 Murashige and Skoog medium including vitamins (Duchefa, M0222), 10 mg L−1 ascorbic acid, 1 mg L−1 2,4-dichlorophenoxyacetic acid, 0.2 mg L−1 zeatin, 30 g L−1 sucrose, 3 g L−1 gelrite, pH 5.8) in the dark at 26 ± 2 °C for 3–6 months. The embryogenic calli were transferred to liquid callus induction medium in conical flasks and incubated in the dark at 26 ± 2 °C, with agitation at 90 rpm. Once the homogenous embryogenic cell suspensions were established, they were maintained and multiplied by regular sub-culturing every 2 weeks and used for generating genome-edited plants.
Guide RNA design and CRISPR/Cas9 plasmid construction
CRISPR guide RNA design: Complete genome sequence of nine badnaviruses, four of which were different BSV strains [banana streak virus strain Obino l’Ewai (BSOLV, GenBank accession KJ013506 and JQ409539), banana streak VN virus (GenBank accession KJ013510), banana streak virus Acuminata Vietnam (NCBI reference sequence NC_007003), banana streak Mysore virus (NCBI reference sequence NC_006955)] affecting banana and five other badnaviruses [rubus yellow net virus (GenBank accession KF241951), gooseberry vein banding-associated virus (GenBank accession HQ852248), grapevine vein-clearing virus (NCBI reference sequence NC_015784), piper yellow mottle virus (GenBank accession KC808712), taro bacilliform virus (NCBI reference sequence NC_004450)] affecting different plant species were downloaded from NCBI. The sequences were aligned to identify conserved regions related to the BSOLV, which is widely distributed in Africa. A relative conserved block of sequence was identified within each ORF of BSOLV strain. The gRNAs were designed for each of CRISPR target from sequences of the BSOLV (GenBank accession KJ013506 and JQ409539) and eBSOLV (GenBank accession HE983609 and AP009334). One gRNA was selected from each ORF (gRNA1-TTGAGCAAAGAAGACGTTG targeting S1 in ORF1; gRNA2-GCTGGAACAACTGGTGACT targeting S2 in ORF2; gRNA3-ATGGAGTTCATATGATCAT targeting S3 in ORF3) based on their specificity to their target site and minimal potential off targets in the Musa genome.
CRISPR/Cas9 plasmid preparation: The gRNA cassette, containing OsU6 promoter followed by two BbsI restriction sites and tracer RNA scaffold, was amplified from pZK_OsU6-gRNA plasmid34 using primers OsU6_Chimera.F1 and OsU6_Chimera.R1 and was cloned into pENTR-D/Topo. The additional BbsI site in pENTR-D/Topo backbone was mutated by amplifying with primers noBbsI.For and noBbsI.Rev (Supplementary Table 1) and self-circularizing using In-Fusion (Takara Bio USA Inc.) according to manufacturer’s instructions to generate pMR185.
One gRNA (Supplementary Table 1) from each ORF of the BSV genomic sequence was synthesized and cloned into pMR185 to generate the gRNA modules as described8. Briefly, for each gRNA sequence, two complementary oligos with 4 bp overhangs were annealed and inserted into BbsI digested pMR185 through T4 ligation.
To allow combining three gRNA modules from the three individual clones (one for each ORF), primers were designed (Supplementary Table 1) to amplify the gRNA modules adding appropriate attB/attBr flanking sites. The amplified fragments were then cloned into the appropriate Multisite Gateway Pro3 pDONR221 vectors by Gateway BP reaction to generate entry clones suitable for MultiSite Gateway LR cloning. Then, all the three gRNAs were combined into one destination vector pDe-CAS9-Kan8 using Gateway LR II Plus clonase reaction (ThermoFisher) to yield the final plasmid pMR215-BSV1. The plasmid has nptII gene as an in planta selection marker and each gRNA was driven by the rice Pol III promoter (OsU6 promoter). The Cas9 gene used in this plasmid was Arabidopsis codon optimized and regulated by parsley ubiquitin promoter (PcUbi).
Delivery of CRISPR/Cas9 plasmid into plantain
The CRISPR/Cas9 construct was introduced into A. tumefaciens strain EHA105 through electroporation. Agrobacterium-mediated transformation of embryogenic cells of Gonja Manjaya was performed with the plasmid construct using the method as described32. The embryogenic cell suspensions of Gonja Manjaya were cocultivated with A. tumefaciens strain EHA105 harboring CRISPR/Cas9 plasmid construct for 3 days in the dark at 22 °C. After 3 days, Agrobacterium-infected cell suspensions was washed with liquid callus induction medium supplemented with 500 mg L−1 cefotaxime, and transferred to selective regeneration medium supplemented with 300 mg L−1 cefotaxime and 100 mg L−1 kanamycin with fortnightly transfer onto fresh medium.
Transformed events regenerated on selective medium containing kanamycin (100 mg L−1) were maintained by routine sub-culturing and multiplied on proliferation medium at 26‒28 °C and a photoperiod of 16 h/8 h (day/night).
PCR analysis to confirm integration of Cas9 gene
Genomic DNA was extracted from putative transgenic events and control non-transgenic plants using the CTAB method as described above. The presence of Cas9 gene was checked by PCR analysis using primers P5 (Supplementary Table 1) amplifying a region between Cas9 and Pea 3 A poly A. PCR was performed in a 20 µL reaction volume containing 1 µL genomic DNA (100 ng µL−1), 2 µL of 10X PCR buffer, 0.4 µL of dNTP, 0.5 µL of 10 µM of each primer, 0.1 µL Taq polymerase, and 15.5 µL nuclease free water. PCR amplification conditions were performed as follows: initial denaturation step at 95 °C for 5 min, followed by 34 cycles of denaturation at 94 °C for 30 s, annealing at 55 °C for 30 s, extension at 72 °C for 1 min, and final extension at 72 °C for 10 min. After amplification, 10 µL of PCR product was resolved on 1% agarose gel stained with gel red.
Molecular analysis of transgenic events
PCR band shift analysis: The mutations in the transgenic events were confirmed by PCR analysis using P1 primer pair (Supplementary Table 1) flanking gRNAs from ORF1 and ORF2 of BSOLV. The primer pair P1 was designed to amplify a DNA fragment of 786 bp surrounding the target sites S1 and S2 (Fig. 1). PCR reaction conditions were: initial denaturation step at 95 °C for 15 min, followed by 34 cycles of denaturation at 94 °C for 30 s, annealing at 57 °C for 30 s, extension at 72 °C for 1 min, and final extension at 72 °C for 10 min. PCR amplicons were resolved on 1% agarose gel stained with gel red and observed for the shift in band size due to the deletions in the target sites.
T7E1 assay: For the T7E1 assay, the DNA fragments containing the targeted sites for ORF1 and ORF2 were amplified using the primer P1 (Fig. 1, Supplementary Table 1) and Phusion high-fidelity DNA polymerase (NEB). The PCR product was then denatured-annealed at 95 °C for 5 min, ramped down to 25 °C at 0.1 °C s−1, and incubated at 25 °C for another 30 min. Annealed PCR products were then digested with 5U of T7E1 for 2 h at 37 °C. The T7E1-digested products were separated by 1% agarose gel electrophoresis and stained with gel red.
Sequence analysis: PCR product obtained using BSOLV-specific primers P1 and P3 flanking target sites were purified with QIAquick PCR purification Kit (Qiagen) according to the instruction manual. The purified products were cloned to pCR™8/GW/TOPO® (Invitrogen) according to the manufacturer’s instruction, transformed to DH5α chemical competent E. coli cells and selected on LB plates containing spectinomycin. Ten clones from each transgenic event were prepared for Sanger sequencing. A total of 20 transgenic events along with four different wild type control plants were sequenced. The sequencing data were analyzed using SnapGene software (www.snapgene.com). All the sequences were aligned with the reference sequence of the BSOLV and wild type non-edited control plant (Supplementary Notes 1 and 2). The wild type control plants were developed along with the genome-edited plants from the same embryogenic cell suspension line.
Analysis of potential off-target mutations
The CRISPR web tool Breaking-Cas35 (http://bioinfogp.cnb.csic.es/tools/breakingcas) was used to check for the potential off target sequences. In addition, a 22 bp-long sequence of each gRNA including the PAM motif was blasted against banana genome A (Musa acuminata) and genome B (M. balbisiana) using the BLASTN program in the Banana Genome Hub36 (http://banana-genome-hub.southgreen.fr) to identify the genes with potential off-target sites. Based on the gene sequences, primers were designed to analyze the potential off-target fragments from seven potential off-target sites, which either showed high sequence similarity with the gRNAs and/or contain the PAM motif. PCR was performed with nine edited events (1, 17, 48, 66, 76, 81, 85, 87, and 97) using these primers. PCR products were purified and cloned to pCR™8/GW/TOPO® (Invitrogen) for sequencing. The sequence chromatograms were analyzed with SnapGene software (www.snapgene.com) and compared with the banana genome and wild type non-edited control plant (Supplementary Note 3).
Evaluation of plants for induction of BSV symptoms
Glasshouse assay under stress conditions: Eight randomly selected genome-edited events (17, 48, 66, 76, 81, 85, 87, and 97) confirmed by sequencing for mutations, and three control wild type plants were acclimatized and grown in pots in the greenhouse. Three replicates per event were used in this assay. Two-month-old potted plants were subjected to water stress for 2 weeks to induce BSV symptoms. At the end of the stress period, the plants were assessed visually for the presence or absence of BSV disease symptoms as chlorosis or yellow streaks. The plants were monitored for 60 days after induction of water stress. The plants were observed for BSV disease symptoms as yellowish streaks and distorted shapes of banana leaves and more severe symptoms including retarded growth, pseudostem split, and leaf necrosis leading to death of infected plants3.
PCR and qPCR analysis for detection of BSOLV: Genomic DNA was extracted from symptomatic (66 and 97) and asymptomatic (17, 48, 76, 81, 85, and 87) edited plants and symptomatic wild type control plants using CTAB33. The leaves showing BSV symptoms from three replicates for each event and three control plants were pooled for DNA extraction. The PCR was performed using P4 primer pair as described above. Primer P4 can detect only the episomal BSOLV and not eBSOLV7. The symptomatic field plants of Agbagba (AAB) and symptomatic wild type control from glasshouse showing BSV symptoms under stress conditions were used as positive controls.
The symptomatic (66 and 97) and asymptomatic (17, 48, 76, 81, 85, and 87) edited plants and symptomatic wild type non-edited control plants from glasshouse were further validated by qPCR for detection of episomal BSOLV using P4 primers. DNA from symptomatic field plant of Agbagba plant and asymptomatic in vitro plantlet of Cavendish Williams were used as positive and negative controls, respectively. qPCR was performed in a 20 µL reaction volume containing 60 ng of genomic DNA, 10 µL of Luna Universal qPCR master mix (New England Biolab), and 0.3 µL of 10 µM of each primer (P4F and P4R, Supplementary Table 1). qPCR was performed using the ABI 7500 real-time machine (Applied Biosystem, USA). Three technical replicates were used for each sample and the experiment was performed twice. The means and standard deviation were calculated for CT values of six technical replicates from two independent experiments, using Minitab 14 statistical software, 2012. Mean CT values of genome-edited plants were compared with wild type non-edited control plant, the positive control, and the negative control.
All data needed to evaluate the conclusions in the paper are present in the paper and/or the supplementary materials. All relevant data are also available from the authors.
Cheng, C., Lockhart, B. E. L. & Olszewski, N. E. The ORF I and II proteins of Commelina yellow mottle virus are virion-associated. Virology 223, 263–271 (1996).
Leclerc, D. et al. The open reading frame III product of cauliflower mosaic virus forms a tetramer through a N-terminal coiled coil. J. Biol. Chem. 273, 29015–29021 (1998).
Fargette, D. et al. Molecular ecology and emergence of tropical plant viruses. Annu. Rev. Phytopathol. 44, 235–260 (2006).
Harper, G., Osuji, J. O., Heslop-Harrison, J. S. & Hull, R. Integration of banana streak virus infection into the Musa genome: molecular and cytogenetic evidence. Virology 255, 207–213 (1999).
Chabannes, M. et al. Three infectious viral species lying in wait in the banana genome. J. Virol. 87, 8624–8637 (2013).
Duroy, P. O., Perrier, X., Laboureau, N., Jacquemoud-Collet, J. P. & Iskra-Caruana, M. L. How endogenous plant pararetroviruses shed light on Musa evolution. Ann. Bot. 117, 625–641 (2016).
Delanoy, M., Salmon, M., Kummert, J., Frison, E. & Lepoivre, P. Development of realtime PCR for the rapid detection of episomal Banana streak virus (BSV). Plant Dis. 87, 33–38 (2003).
Ritter, A. et al. The transcriptional repressor complex FRS7-FRS12 regulates flowering time and growth in Arabidopsis. Nat. Commun. 8, 15235 (2017).
Iskra-Caruana, M. L., Baurens, F. C., Gayral, P. & Chabannes, M. A four-partner plant–virus interaction: enemies can also come from within. Mol. Plant–Microbe Interact. 23, 1394–1402 (2010).
Zaidi, S. S., Tashkandi, M., Mansoor, S. & Mahfouz, M. M. Engineering plant immunity: using CRISPR/Cas9 to generate virus resistance. Front. Plant Sci. 7, 1673 (2016).
Ali, Z., Ali, S., Tashkandi, M., Zaidi, S. S. & Mahfouz, M. M. CRISPR/Cas9-mediated immunity to geminiviruses: differential interference and evasion. Sci. Rep. 6, 26912 (2016).
Baltes, N. J. et al. Conferring resistance to geminiviruses with the CRISPR–Cas prokaryotic immune system. Nat. Plants 1, 15145 (2015).
Gomez, M. A. et al. Simultaneous CRISPR/Cas9‐mediated editing of cassava eIF4E isoforms nCBP‐1and nCBP‐2 reduces cassava brown streak disease symptom severity and incidence. Plant Biotechnol. J. 17, 421–434 (2019).
Zhang, T. et al. Establishing RNA virus resistance in plants by harnessing CRISPR immune system. Plant Biotechnol. J. 16, 1415–1423 (2018).
White, M. K., Hu, W. & Khalili, K. The CRISPR/Cas9 genome editing methodology as a weapon against human viruses. Discov. Med. 19, 255–262 (2015).
Kennedy, E. M. et al. Suppression of hepatitis B virus DNA accumulation in chronically infected cells using a bacterial CRISPR/Cas RNA-guided DNA endonuclease. Virology 476, 196–205 (2014).
Kennedy, E. M. et al. Inactivation of the human papillomavirus E6 or E7 gene in cervical carcinoma cells by using a bacterial CRISPR/Cas RNA-guided endonuclease. J. Virol. 88, 11965–11972 (2014).
Yuen, K. S. et al. CRISPR/Cas9-mediated genome editing of Epstein-Barr virus in human cells. J. Gen. Virol. 96, 626–636 (2015).
Hu, W. et al. RNA-directed gene editing specifically eradicates latent and prevents new HIV-1infection. Proc. Natl. Acad. Sci. USA 111, 11461–11466 (2014).
Ophinni, Y., Inoue, M., Kotaki, T. & Kameoka, M. CRISPR/Cas9 system targeting regulatory genes of HIV-1 inhibits viral replication in infected T-cell cultures. Sci. Rep. 8, 7784 (2018).
Naim, F. et al. Gene editing the phytoene desaturase alleles of Cavendish banana using CRISPR/Cas9. Transgenic Res. 27, 451–460 (2018).
Char, S. N. et al. An Agrobacterium-delivered CRISPR/Cas9 system for high-frequency targeted mutagenesis in maize. Plant Biotechnol. J. 15, 257–268 (2017).
Lowder, L. G. et al. A CRISPR/Cas9 toolbox for multiplexed plant genome editing and transcriptional regulation. Plant Physiol. 169, 971–985 (2015).
Gao, J. et al. CRISPR/Cas9-mediated targeted mutagenesis in Nicotiana tabacum. Plant Mol. Biol. 87, 99–110 (2015).
Brooks, C., Nekrasov, V., Lippman, Z. B. & Van Eck, J. Efficient gene editing in tomato in the first generation using the clustered regularly inter spaced short palindromic repeats/CRISPR-associated 9 systems. Plant Physiol. 166, 1292–1297 (2014).
Kaur, N. et al. CRISPR/Cas9-mediated efficient editing in phytoene desaturase (PDS) demonstrates precise manipulation in banana cv. Rasthali genome. Funct. Integr. Genom. 18, 89–99 (2017).
Tripathi, L. et al. Field trial of Xanthomonas wilt disease-resistant bananas in East Africa. Nat. Biotechnol. 32, 868–870 (2014).
Fu, Y. et al. High-frequency off-target mutagenesis induced by CRISPR-Cas nucleases in human cells. Nat. Biotechnol. 31, 822–826 (2013).
Hohn, T. & Fütterer, J. The proteins and functions of plant pararetroviruses: knowns and unknowns. Crit. Rev. Plant Sci. 16, 133–161 (1997).
Liao, H. K. Use of the CRISPR/Cas9 system as an intracellular defense against HIV-1 infection in human cells. Nat. Commun. 6, 6413 (2015).
Kaminski, R. et al. Elimination of HIV-1 genomes from human T-lymphoid cells by CRISPR/Cas9 gene editing. Sci. Rep. 6, 22555 (2016).
Tripathi, J. N., Muwonge, A. & Tripathi, L. Highly efficient regeneration and transformation protocol for plantain cv. Gonja Manjaya (Musa spp. AAB) using embryogenic cell suspension. In Vitro Cell. Dev. Biol. Plant 48, 216–224 (2012).
Stewart, C. N. & Via, L. E. A rapid CTAB DNA isolation technique useful for RAPD fingerprinting and other PCR applications. Biotechniques 14, 748–750 (1993).
Mikami, M., Toki, S. & Endo, M. Comparison of CRISPR/Cas9 expression constructs for efficient targeted mutagenesis in rice. Plant Mol. Biol. 88, 561–572 (2016).
Oliveros, J. C. Breaking-Cas—interactive design of guide RNAs for CRISPR-Cas experiments for ENSEMBL genomes. Nucleic Acids Res. 44, W267–W271 (2016).
Droc, G. et al. The banana genome hub. Database bat035 (2013). https://doi.org/10.1093/database/bat035.
Authors wish to thank the Roots, Tubers, and Banana program of CGIAR (CRP_RTB) for funding support, Jim Lorenzen, Senior Program Officer at Bill and Melinda Gates Foundation for useful discussion on genome editing of banana for knocking out eBSV, Masaki Endo for providing pZK_OsU6-gRNA plasmid, Lava Kumar, Virologist at IITA for providing DNA of BSV symptomatic plant and discussion on screening of genome-edited events for disease evaluation, and Sarah Macharia, Everlyne Omollo, and Mark Adero for technical assistance.
The authors declare no competing interests.
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Tripathi, J.N., Ntui, V.O., Ron, M. et al. CRISPR/Cas9 editing of endogenous banana streak virus in the B genome of Musa spp. overcomes a major challenge in banana breeding. Commun Biol 2, 46 (2019). https://doi.org/10.1038/s42003-019-0288-7
Frontiers in Plant Science (2020)
Phytopathology Research (2020)
Introducing Chemistry Students to Emerging Technologies in Gene Editing, Their Applications, and Ethical Considerations
Journal of Chemical Education (2020)
Current Plant Biology (2020)
Journal of Basic Microbiology (2020)