Skip to main content

Improved CRISPR/Cas9 gene editing by fluorescence activated cell sorting of green fluorescence protein tagged protoplasts



CRISPR/Cas9 is widely used for precise genetic editing in various organisms. CRISPR/Cas9 editing may in many plants be hampered by the presence of complex and high ploidy genomes and inefficient or poorly controlled delivery of the CRISPR/Cas9 components to gamete cells or cells with regenerative potential. Optimized strategies and methods to overcome these challenges are therefore in demand.


In this study we investigated the feasibility of improving CRISPR/Cas9 editing efficiency by Fluorescence Activated Cell Sorting (FACS) of protoplasts. We used Agrobacterium infiltration in leaves of Nicotiana benthamiana for delivery of viral replicons for high level expression of gRNAs designed to target two loci in the genome, NbPDS and NbRRA, together with the Cas9 nuclease in fusion with the 2A self-splicing sequence and GFP (Cas9-2A-GFP). Protoplasts isolated from the infiltrated leaves were then subjected to FACS for selection of GFP enriched protoplast populations. This procedure resulted in a 3–5 fold (from 20 to 30% in unsorted to more than 80% in sorted) increase in mutation frequencies as evidenced by restriction enzyme analysis and the Indel Detection by Amplicon Analysis, which allows for high throughput profiling and quantification of the generated mutations.


FACS of protoplasts expressing GFP tagged CRISPR/Cas9, delivered through A. tumefaciens leaf infiltration, facilitated clear CRISPR/Cas9 mediated mutation enrichment in selected protoplast populations.


CRISPR/Cas has emerged as a powerful tool for precise genetic editing (PGE) in a wide range of organisms [1], including plants [2]. CRISPR/Cas relies on the Cas DNA nuclease being guided by the small guide RNA (gRNA), to make a double stranded break (DSB) at the desired place in the genome (reviewed in [3]) leading to activation of inherent repair mechanisms (Non-Homologous End Joining (NHEJ) or Homologous Recombination (HR) if a DNA molecule with identical flanking sequences is co-delivered. CRISPR/Cas mediated PGE in plants may be complicated by the presence of complex and high ploidy genomes or by inefficient or poorly controlled delivery of PGE components to gamete cells or cells with regenerative potential. Moreover, subsequent regeneration and tissue culturing post PGE is often lengthy, labor-intensive and prone to produce random somatic mutations and targeted insertion mediated mutagenesis through homologous recombination is still a main challenge within PGE [2]. There is therefore a demand for optimizing PGE in plants towards efficient generation and propagation of stable heritable editing at the organism level.

Nucleic acids may be introduced into plant cells/tissues by biolistic particle bombardment [4], which, however, often results in insertion of multiple copies at multiple sites in the genome [5]. Other strategies include transformation of protoplast by chemical means using polyethylene glycol (PEG) in combination with calcium ions or by electroporation (reviewed in [5]), where the latter requires elaborate tissue culturing for regeneration to fertile plants and may introduce genetic instability and resulting somaclonal variation. PEG-mediated transformation, in particular, has been used to deliver constructs encoding the PGE components, incl. Zinc Finger-Nucleases (ZFNs) [6], Transcription activator-like effector nucleases (TALENs) [7, 8] and CRISPR/Cas9 [8, 9] and lately also for delivery of the Cas9 enzyme and associated gRNA into plant cell protoplasts in vitro [10]. Excess DNA is regularly used for PEG mediated transformation of protoplasts (typically in molar ratios of 1: 1–2 × 107 (protoplast: plasmid DNA) [11]) and has been reported to confer unintended random integrations in the recipient genomes [12]. Agrobacterium-mediated transformation on the other hand is generally perceived to be an efficient and a more controlled way of delivering transgenes [13] and the use of strains, with putatively downregulated integration capacity [14], in combination with down-regulation of host factor integration genes may facilitate alternative ways of non-integrative delivery of PGE components. Also, Agrobacterium may in some cases be the only viable option for delivering of transgenes. In recent years, Agrobacterium-mediated delivered viral constructs has attracted increasing interest because of their high copy number and resulting expression capabilities [15, 16]. Deconstructed viral vectors (replicons) have proved extremely effective for rapid, high-yield production of a number of pharmaceutical proteins, of which some are currently undergoing clinical evaluation [16]. As efficient gene editing relies on PGE component expression, virus replicons have likewise attracted attention as delivery vehicles [17]. Deconstructed geminivirus type replicons (as delivery vehicles) have been shown to generate mutations in the solanaceous species Nicotiana benthamiana [17] and Solanum lycopersicum (tomato) [18] and recently in Triticum aestivum (wheat) [19]. N. benthamiana can be grown in high density and still produce large amounts of biomass in a matter of weeks [16], and has a track record for production of therapeutic glycoproteins in mg scale ([20,21,22,23,24]) through the use of leaf or leaf disc infiltration [25]. In addition, N. benthamiana may readily be subjected to protoplast transformation [26] and explant/protoplast regeneration [27, 28]. Several approaches have been reported to confer enrichment of PGE mutations in cells. Fluorescence Activated Cell Sorting (FACS) of edited cells, for example, is regularly used as means of PGE mutation enrichment in mammalian cell systems [29], and the present study addresses the feasibility of applying this strategy to plant cells.

So far reports on FACS and post FACS cultivation of plant protoplasts are relatively scarce [30], due to the removal of the rigid and structure providing cell wall, which otherwise stabilize the cell integrity [31,32,33]. The present study explores the combined use of Agrobacterium-mediated delivery of viral replicons for expression of GFP tagged gRNA/Cas9 in leaves of N. benthamiana with FACS in order to obtain protoplast populations with significantly increased gene editing.


The overall strategy for Agrobacterium-mediated delivery of deconstructed replicons expressing gRNA/Cas9-2A-GFP in leaves of N. benthamiana combined with FACS of GFP expressing protoplasts is outlined in Fig. 1.

Fig. 1
figure 1

Scheme for Agrobacterium-mediated in-leaf GFP tagged CRISPR/Cas9 mutation generation combined with FACS enrichment of GFP expressing protoplasts. a Guide RNA (gRNA) target sequence may be selected on the basis of in silico prediction analysis and the presence of a Restriction Enzyme (RE) recognition motif spanning the SpCas9 cleavage site (− 3 bp upstream of the protospacer adjacent motif (PAM) [34]) for fast RE-mediated mutation screening. Primers flanking the gRNA target for PCR mediated-mutation scoring are indicated. The deconstructed bean yellow dwarf virus (BeYDV) replicon is produced from the Agrobacterium tumefaciens delivered T-DNA, that contains the viral cis-acting Long (LIR) and Short Intergenic Regions (SIR) in a Long-Short-Long region (pLSL) arrangement, which together with the co-expressed trans acting Rep/RepA replication initiation proteins facilitate replicational release and Gemini Virus Replicon (GVR) circularization allowing joining of the two BeYDV replicon LIR elements within plant cell nuclei [17]. Abbreviations: Left and Right T-DNA border, LB & RB, Cauliflower mosaic virus 35S promoter, CMV35S, Arabidopsis thaliana U6 promotor, AtU6-Pro [35, 36], hygromycin phosphotransferase, HPT, Streptococcus pyogenes Cas9, SpCas9, nopaline synthase terminator, NOS, Nucleus Localization Signal, NLS, 2A self-cleaving sequence of foot-and-mouth disease virus (FMDV), 2A [37, 38], Agrobacterium tumefaciens, A. tumefaciens, Nicotiana benthamiana, N. benthamiana. The replicon constructs (a) are transformed into A. tumefaciens by electroporation, grown under selection overnight and re-suspended in infiltration buffer to a final total OD600 of ca. 0.2 where after the abaxial side of young expanding leaves of 3–4 week old N. benthamiana plants are infiltrated with the agrobacterium strain carrying the construct of interest using a syringe and left for 2–4 days (b). Protoplasts are isolated (c) and subjected to florescence microscopy (for estimation of protoplast isolation and transformation efficiencies) and to Fluorescence Activated Cell Sorting (FACS) (d) of GFP (SpCas9-2A-GFP) expressing protoplasts for mutation enrichment. The target region on the genome is amplified by PCR (e) with mutations scored by the high throughput screening technique Indel Detection by Amplicon Analysis (IDAA) [39] (f), which allows for detection of down to 1 bp deletions and insertions (indels) and by restriction enzyme (RE) analysis (g), which monitors resistant mutated RE recognition/cleavage sites. Optionally, explants with stable PGE editing can be obtained by embedding the protoplasts in alginate, followed by callus induction and shoot regeneration as outlined in [40]. Protoplasts shown in (c) are presented as light-, fluorescent micrographs and overlay hereof

gRNA and replicon construct design

In the present study we targeted the Nicotiana benthamiana PHYTOENE DESATURASE (NbPDS) and REDUCED RESIDUAL ARABINOSE arabinosyl transferase (NbRRA) loci, orthologous to the Arabidopsis thaliana arabinosyltransferase encoding genes involved arabinosylation of plant cell wall extensins (AtRRA1–3) [41, 42], which have a proven [43] and an untested CRISPR/Cas9 editing record, respectively (Fig. 2a). gRNA target sequences were confined to early exons and identified on the basis of in silico prediction analysis (, [45]), and the presence of a Restriction Enzyme (RE) recognition sequence spanning the predicted cut site of SpCas9–3 bp upstream of the protospacer adjacent motif (PAM) [34] for RE-mediated mutation screening.

Fig. 2
figure 2

NbRRAall1/NbPDS2-gRNA generated indels. a gRNA targets of the N. benthamiana loci, REDUCED RESIDUAL ARABINOSE arabinosyl transferase (NbRRA) and PHYTOENE DESATURASE (NbPDS), were a BtgI and a AvrII site is situated 2 and 0 bp upstream of the protospacer adjacent motif (PAM), respectively. Given the predicted cut site of SpCas9, 3 bp upstream of the PAM sequence [44], all of the NbPDS2-gRNA derived mutation combinations will destroy the AvrII site in the NbPDS target site and only insertions starting with ‘G’ at the cut site, i.e. less than one fourth the insertions possible, will restore the BtgI site in the NbRRAall1 target site. Primers, flanking the gRNA targets, are indicated by arrows. b Western blot analysis of day 4 post infiltration leaves using anti Flag and anti GFP mAbs, cross-reacting to SpCas9 (154 kDa) and to a faint protein band corresponding to the un-cleaved fusion protein (SpCas9-2A-GFP, ca 180 kDa), respectively. c, d DNA from 4 days post infiltration leaf samples of NbRRAall1- and NbPDS2-gRNA/Cas9 infiltrations were isolated, PCR amplified and subjected to restriction enzyme digestions using BtgI (NbRRAall1) and AvrII (NbPDS2), respectively, with the resistant bands (indicated by arrow) isolated, cloned into pJet and 12 clones of each target sequenced revealing the resulting indels depicted

A deconstructed immobilized mild strain of the bean yellow dwarf virus (BeYDV), allowing for a high replicon copy number in the nucleus, has recently been used to construct an Agrobacterium T-DNA that integrates into the host cell chromosome and delivers a geminivirus replicon (GVR) [17, 46]. The minimal immobilized replicons are delivered by Agrobacterium infiltration (here to N. benthamiana leaves) along with co-infiltrated constructs for expression of replicon trans-acting replication initiation proteins (Rep or RepA) [47] (Fig. 1a). While the replicons are non-integrative and transiently expressed the initial Agrobacterium T-DNA (LB-RB) delivery of the replicon is integrative [17]. Lately, GVRs were constructed and used to propagate and express PGE components, such as ZFNs and TALENs and CRISPR/Cas9 [17]. In the present study, we inserted the Streptococcus pyogenes Cas9 enzyme (SpCas9) [48] in translational fusion with the 2A self-splicing sequence of the foot-and-mouth disease virus [37, 38] and GFP [49] (SpCas9-2A-GFP) under control of the CMV 35S promoter and the gRNAs under control of the AtU6 promoter [35, 36] in the BeYDV GVR replicon [17] as depicted in Fig. 1a and detailed in the Methods section.

In-leaf gRNA/Cas9 generated mutations

The SpCas9-2A-GFP/gRNA expressing GVR replicons (Fig. 1a) targeting the NbPDS and NbRRA loci (Fig. 2a) were electroporated into Agrobacterium tumefaciens and grown under selection overnight and re-suspended in infiltration buffer to a final total OD600 of 0.2 where after the abaxial sides of young expanding leaves of N. benthamiana were subjected to Agrobacterium infiltration. The infiltrated plants were left for 2–4 days allowing for gRNA/Cas9 expression and mutation generation within the intact leaves. Western blot analysis of total proteinacious extracts, using anti-Flag and anti-GFP mAbs as primary antibodies against the Flag- and GFP-tagged SpCas9 revealed the presence of a distinct band at the expected MW (154 kDa) of mature SpCas9 with a faint band corresponding to the un-cleaved fusion protein (SpCas9-2A-GFP, ca 180 kDa) in infiltrated leaves demonstrating expression and efficient 2A mediated auto cleavage of SpCas9-2A-GFP (Fig. 2b).

RE-mediated mutation analysis of PCR fragments, using primers flanking the gRNA target sites, revealed the presence of non-digestible bands indicative of mutated RE recognition/cleavage sequence for the two target sites (Fig. 2c and d). The RE resistant band of each locus was isolated, sub-cloned and sequenced with the presence of insertions or deletions (indels) demonstrated (Fig. 2c and d).

Protoplast isolation and FACS-mediated mutation enrichment

Protoplasts of WT and infiltrated leaves were essentially obtained using the protocol devised by Dovzhenko et al. 1998 [27] with minor modifications as outlined in the Methods section. Protoplast quality and yield varied significantly apparently influenced by growth conditions pre- and post-infiltration. Here a temperature of 22–24 °C, a 16 h/8 h (light/dark) regime of moderate sunlight (See ‘Growth conditions’, Methods section) generally conferred a high amount of intact protoplasts. Protoplast integrity and transformation were assessed by comparative bright field and fluorescent microscopy frequently with varying estimated transformation rates of 20 - > 80% (Additional file 1: Figure S1). GFP fluorescence accumulated in particular in the cytoplasmic strands and the contours of the cell (Additional file 1: Figure S1), which is in agreement with a cytoplasmic 2A-mediated release of GFP. This was corroborated by the western blot analysis (Fig. 2b) showing presence of the mature SpCas9 with only traces un-cleaved product. Also, in agreement with soluble non-tagged GFP being able to pass into and accumulate in nuclei [50], some accumulation in nuclei structures was observed (Additional file 1: Figure S1).

FACS of fluorescent protoplasts were done using a FACSAria III (BD Biosciences) apparatus with settings to accommodate for the approximate size of N. benthamiana protoplasts [51] as described in the Methods section. Two fluorescent enriched populations, protoplasts with medium GFP intensity (P4) and with high intensity (P5), were selected for sorting corresponding to 17% & 10 and 14% & 5% of the total population for the NbRRAall1-gRNA/SpCas9-2A-GFP and NbPDS2-gRNA/SpCas9-2A-GFP infiltrations, respectively (Fig. 3b). RE analysis of PCR amplicons suggested an estimated indel frequency of unsorted, P4 and P5 sorted populations of 20–30, 50% and 70–80% for the NbRRAall1-gRNA and 40, 50 and > 80% for the NbPDS2-gRNA (Fig. 3c). This was corroborated by Indel Detection by Amplicon Analysis (IDAA) (Fig. 3d) and sequence analysis of the cloned PCR fragments of the two P5 populations (10 clones of each), which showed indel to WT ratio of 60 and 70%, respectively. The indel distributions obtained for the NbRRAall1- and NbPDS2-gRNA infiltrations − 3(1), − 1(2) & + 1(4) and − 1(3) & + 1(3) (Fig. 3e), respectively, are in agreement with earlier findings for SpCas9 mediated mutations in plants [52].

Fig. 3
figure 3

FACS mediated enrichment of gRNA/SpCas9 expressing protoplast cells and resulting mutations. 3–5 N. benthamiana leaves were infiltrated with Agrobacterium-delivered replicons expressing SpCas9-2A-GFP together with NbRRAall-gRNA or NbPDS2-gRNA (a), respectively, and left for 2–4 days. b WT protoplasts and protoplasts expressing SpCas9-2A-GFP and NbRRAall1- or NbPDS2-gRNA were subjected to GFP mediated FACS. The DAPI and FITC intensities for protoplasts were recorded and three populations, P3, P4 & P5, with the P3 population corresponding to non-transformed cells and the P4 and P5 populations representing intermediary and high stringently sorted cell populations, respectively, were selected from the Dot Scattering Chromatograms. Transfected protoplasts were defined as FITC-positive events and gates were set to separate WT and GFP enriched protoplast populations, using the WT sample to define non-transfected wild type populations (P3) in the transfected samples. P4 and P5 (GFP enriched populations) were gated with medium and high FITC signal intensity. c RE analysis of PCR amplified target regions using BtgI and AvrII for NbRRAall1- and NbPDS2 gRNAs, respectively, demonstrating indel formation in unsorted and indel enrichment in FACS sorted (P4 and P5) populations. Indel enrichment in P5 populations were cooperated by the IDAA technique (d) where the additional restriction enzyme digest allows for visualization of the mutated population without the presence of non-mutated PCR amplicons (‘mutated / RE resistant’ designates the RE site were mutated rendering it resistant for digestion while ‘WT/cut’ designates WT sites that were cut and moved downstream in the chromatogram). e Sequence analysis of RE-resistant PCR fragments of the two P5 populations. f Post viability of protoplasts was assessed in WT protoplasts (dark circular objects) without detectable GFP signal and GFP fluorescence in Cas9-2A-GFP sorted protoplasts (presented as an overlay of light and fluorescent micrographs). FACS was carried out using a FACSAria III (BD Biosciences) apparatus with procedure and parameters as outlined in the Methods section and IDDA as described in [39]. For the viability test shown in f crude protoplasts were prepared and sorted on a Sony Cell sorter SH800S with sorting gating parameters similar to those used on the BD FACSAria III sorter and with the W5 buffer as recipient buffer

Bright field microscopy suggested that 10–20% of a WT protoplast population was intact after FACS when using PBS buffer as sheath fluid and MMM550 as recipient buffer (Additional file 2: Figure S2). Viability of post FACS GFP positive protoplasts was assessed by bright field and fluorescence microscopy (Fig. 3f) and confirmed by propidium iodide exclusion assays (Additional file 2: Figure S2). Sorting into PBS as recipient buffer resulted in instant lysis as evidenced by bright field microscopy (data not shown). In agreement with ribonucleoprotein, i.e. in vitro transcribed gRNA and heterologous expressed Cas9, conferring nuclease activity in-vitro [10, 53], we tested (post FACS) for PBS lysis mediated editing activity and found an 2–3 fold increased editing when the PBS lysed protoplasts were left in PBS for 2 h at room temperature (Additional file 3: Figure S3). All post FACS protoplast samples were immediately incubated on ice accordingly. Potential continued editing in the timespan from FACS to further processing may on the other hand likewise increase the ‘in cell’ editing.

Embedment of GFP transformed protoplasts in alginate with initial callus formation (Additional file 4: Figure S4) demonstrated the feasibility of obtaining gene edited lines using an explant shoot regeneration systems as described in [40].


The use of PGE in plants may be complicated by the presence of complex genomes and inefficient or poorly controlled delivery of PGE components to gamete cells or recipient pluripotent cells. DNA encoding PGE components may be delivered to the plant cell either directly, i.e. by biolistic transformation or protoplast transformation (reviewed by [5]), or indirectly, mainly via bacteria, usually Agrobacterium tumefaciens or (less commonly) Agrobacterium rhizogenes [54], which is generally perceived to be a controlled way of delivering transgenes [55]. Virus replicons provide high copy number of expression units and thus a means of significantly boosting PGE component expression levels [46, 56] and methods for increasing identification/selection of PGE-edited cells have been introduced and applied successfully e.g. for mammal cells [29].

In the present study we combined Agrobacterium-mediated delivery of a viral replicon expressing GFP labeled gRNA/SpCas9 for generation of in-leaf mutations with the use of FACS of GFP-fluorescent protoplasts for enrichment of mutated protoplast populations. BeYDV GVR replicons, expressing gRNAs targeting the NbPDS and the NbRRA locus in N. benthamiana, respectively, together with the SpCas9 nuclease, fused to the 2A self-splicing sequence and GFP (SpCas9-2A-GFP), were introduced into leaves of N. benthamiana by Agrobacterium-mediated infiltration and left for expression and mutation generation within the intact leaf. In leaf expression of the GFP- and Flag-tagged SpCas9 enzyme was readily verified by western blot analysis and generated mutations as evidenced by the presence of restriction enzyme (RE) resistant bands of PCR amplicons comprising the mutated recognition site were cooperated by cloning and sequence analysis of the RE resistant bands. Indel distribution was found to be in accordance with earlier studies for SpCas9 mediated genome editing in plants [52]. With the aim of selecting and concentrating edited cells GFP-expressing protoplasts were isolated and subjected to FACS. The two fluorescence-enriched populations were selected for FACS with the most stringently sorted population yielding a 3–5 fold enrichment in mutations as evidenced by RE-mediated mutation analysis of PCR amplicons and sequence analysis.

The IDAA method allows for fast and direct assessment of indel prevalence and distribution [39]. In the current study, IDAA was combined with RE analysis for visualization of the isolated mutation population, where potential single nucleotide substitutions within the RE recognition site will otherwise co-migrate with the WT peak. While the observed > 50% reduction of FAM-fluorescence signal in IDAA analyses of overnight RE digestions may complicate absolute peak quantification between samples, quantification of the WT peak and indel peak(s) within single samples, provides a means of estimating relative mutation efficacies between samples. The combined use of the RE analysis and the IDAA technique adds an extra analytical layer to the versatile IDAA technique. 10–20% of the protoplasts appeared to be intact in post FACS populations, when PBS and the MMM550 buffer were used as sheath fluid and recipient buffer, respectively. This ratio may, however, be increased by replacing the sheath fluid PBS buffer with a more osmotically favorable buffer and, if feasible on the FACS apparatus used, decrease shearing forces by lowering the psi. Here FACS on a Sony SH3800S cell sorter yielded ample intact protoplasts post FACS probably due to the available 130 μm sorting chip with accordingly lower psi. Isolation of non-ruptured protoplasts through the use of a sucrose gradient significantly aids the identification of protoplast populations with and without GFP-expression. Once an initial delineation of the protoplast populations on the cell sorter has been established, this step may potentially be omitted.

Extracellular gRNA/Cas9 activity from lysed protoplasts, e.g. mediated by FACS sorting, was significant and this residual activity, which may lead to an over estimated indel frequency, was abolished by incubation on ice or FACS sorting into an RNAse containing or protein denaturing buffer.

Also, in this study protoplast yields were generally highly variable. A recent study on Agrobacterium infiltration-mediated expression of a reporter in leaves of N. benthamiana recommended infiltration of more plants but less leaves and sample more positions on the leaf as opposed to run a high number of technical replicates [57]. In addition, it is conceivable that Agrobacterium infection/pathogenesis may affect intact protoplasts yield.

Recently, an unexpected high level of integrations in the recipient genomes associated with PEG-mediated plasmid transformation of protoplasts was reported [12]. Further optimization of the here devised PGE approach may include exploring the use of integration deficient Agrobacterium strains [58] or Virus Induced Gene Silencing (VIGS) mediated down-regulation of host plant factors [59] also important for T-DNA integration, as means of non-integrative delivery of the PGE components [60]. The obtained mutation enrichment may facilitate mutation detection e.g. in situations where the activity of a particular gRNA is weak and reduce laborious explant generation and screening steps. Alternatively, the protoplast based PGE system may be used e.g. in promoter-reporter editing test-screens.


The present study outlines a strategy for enrichment of CRISPR/Cas9 editing in leaf protoplasts. GFP tagged gRNA/Cas9 (gRNA/Cas9-2A-GFP) was delivered by Agrobacterium-infiltration to leaves of N. benthamiana and protoplasts isolated. Subsequent FACS of GFP expressing protoplasts resulted in several fold mutation enrichment in the selected fluorescence enriched protoplast populations.


Growth conditions

Seeds of wild-type Nicotiana benthamiana were sown and grown in soil (Pindstrup substrate number 2) for 4 weeks in greenhouse with a 16/8 h light/dark cycle, app. 70% relative humidity and a day/night temperature cycle of 24 and 17 °C.

2 days prior to infiltration plants were subjected to regular sunlight at a photosynthetic flux of 20–40 μmol photons m− 2s− 1, Photosynthetic Active Radiation (PAR): 20.5 μE.m− 2s− 1, Red – Far Red ratio (R:FR): 1,69), 22–24 °C temperature, an app. 16 h/8 h (light/dark) diurnal rhythm and 70% relative humidity, which were also imposed in the post-infiltration period.

Vectors and construct designs

Descriptive naming of vectors, constructs, primers and primer sequences are provided in Additional file 5: Table S1. The vector pLSLGFP-R (V82), described in [17], containing GFP insert in front of the CMV35S promoter and Gateway destination site in front of CMV 35S promoter-LIR, respectively, was kindly provided by Nicholas Baltes, Michigan University, US. The gRNA Gateway entry vector V26 (pUC57_attL1-AtU6:BbsI-BbsI-tracr-TT_AttL2) was synthesized by Genscript. To obtain insertion of NbRRAall1- or NbPDS2 gRNAs V26 was linearized with BbsI and gRNA targets NbRRAall1 and NbPDS2 were inserted by ligation of the annealed oligonucleotides P042 & P043 and P149 & P150, respectively, yielding V207 (attL1-AtU6: NbPDS2-tracr-TT_AttL2) and V208 (attL1-AtU6: NbRRAall1-tracr-TT_AttL2). V207 and V208 were linearized using EcoRI and cloned together with the Streptococcus pyogenes Cas9 (SpCas9) fragment [11], which was PCR amplified from HBT-Cas9 (Gift from Jen Sheen, Harvard Medical School) using the primer-set P077 & P212, the GFP-Nos fragment amplified from pLSLGFP. R using the primerset L1 & L2, all together using the In-fusion cloning kit (Clontech), yielding V197 (pUC57_AttL1-SpCas9-2A-GFP-Nos; AtU6-NbRRAall1-gRNA-TT.AttL2) and V198 (pUC57_AttL1-SpCas9-2A-GFP-Nos; AtU6-NbPDS2-gRNA-TT-AttL2). V197 & V198 was gateway cloned using pLSL_v2 as destination vector yielding V199 (pLSL_V2_ LIR-AttB1-SpCas9-2A-GFP-Nos; AtU6-NbRRAall1-gRNA-TT-AttB2 SIR-35S-LIR) and V200 (pLSL_V2_ LIR-AttB1-SpCas9-2A-GFP-Nos; AtU6-NbPDS2-gRNA-TT-AttB2 SIR-35S-LIR), respectively. V199 & V200 will, when co-expressed with pREP, express SpCas9 in fusion with the 2A self-splicing sequence of the foot-and-mouth disease virus (FMDV) [37, 38] and GFP [49] (SpCas9-2A-GFP) under control of the CMV35S promoter.

For GFP expression only V82 (pLSLGFP-R_v2) was used.

PDS (NbPDS) and RRA (NbRRA) target loci in the N. benthamiana chromosome

N. benthamiana genes were obtained from [61] based on homology with the Arabidopsis thaliana genes. As N. benthamiana is allotetraploid both chromosome variations of a gene in the given locus are obtained. In contrast to e.g. the presence of 1 and 3 isogenes of AtPDS [62] and AtRRA [41, 42] in diploid Arabidopsis, respectively, NbPDS and NbRRA appear to be single gene loci in allotetraploid N. benthamiana.

The NbRRA gene SolGenomics: Niben101Scf18348 with exons (33526..33687, 35895..36708 & 36767..37113) and Niben101Scf09172 with exons (260530..260692, 261438..262553) with the NbRRAall1-gRNA situated in exon 2 (35905..35924, 261512..261531).

The NbPDS gene SolGenomics: Niben101Scf14708 with exons (13814..14036, 14118..14251, 15346..15435, 16328..16386, 16604..16760, 17017.. 17166, 17412..17532, 17695..17909 & 18003..18104) and Niben101Scf01283 with exons (198006..198228, 198317..198449, 199413..199501, 200074..200127, 200369..200501, 200792..200940, 201104..201223, 201388..201601 201694..201796, 202066..202113 & 202983..203028) and with the NbPDS2-gRNA situated in Exon 3 (15409..15428, 199476..199495).

Agrobacterium mediated leaf infiltration and expression in Nicotiana benthamiana

Agrobacterium tumefaciens pGV3850, harboring constructs (pREP, p19 and (pLSL_V2_ LIR-AttB1-SpCas9-2A-GFP-Nos; AtU6-NbRRAall1-gRNA-TT- AttB2 SIR-35S-LIR (V199) or pLSL_V2_LIR-AttB1-SpCas9-2A-GFP-Nos; AtU6-NbPDS2-gRNA-TT-AttB2 SIR-35S-LIR (V200) and empty vector control were inoculated in 5 mL YEP media with kanamycin (50 mg/L) and rifampicillin (50 mg/L) and incubated at 28 °C, 250 rpm for 24 h. Cells were harvested by centrifugation for 20 min at 4000×g and re-suspended in infiltration buffer (10 mM MES (Sigma-Aldrich), 10 mM MgCl2 and 10 μM acetosyringone (3′,5′-Dimethoxy-4′-Hydroxyacetophenone, Sigma-Aldrich) to a final OD600 of ~ 0.2 and incubated for 3 h at room temperature.

The abaxial side of 3–5 young expanding leaves (4–6 × 6–8 cm (Width, Length)) of N. benthamiana was infiltrated with A. tumesfaciens pGV3850 containing the various constructs and co-infiltrated with the p19 construct [63] (Final OD600 = 0.2) essentially as described by Sainsbury and Lomonossoff (2008) [64], and left for 2–4 days depending on the experimental setting.

Protoplast isolation

Protoplasts were obtained using the protocol devised by Dovzhenko et al. 1998 [27]. Inoculated N. benthamiana leaves for subsequent protoplast-alginate embedment were sterilized by dipping in 96% ethanol and floating in 1.5% hypochlorite solution for 15 min. 3–5 leaves were cut into 0.5–1 mm strips with a scalpel and submerged in 10 ml enzyme solution (400 mM mannitol, 20 mM MES-KOH, pH 5.7, 20 mM KCl, supplemented with 1% Cellulase R10 (w/v) (Duchefa Biochemie, C8001), 0.25% Macerozyme (Duchefa Biochemie, C8002), heated to 55 °C, 10 min, then supplemented with 10 mM CaCl2 and 0.1% BSA) and incubated 2–5 h at 26 °C, 100 rpm, then filtered through a 100 μm filter into a 50 ml Falcon tube, centrifuged for 5 min at 100×g, where after the supernatant was poured off and the protoplast-containing pellet was re-suspended in 3 ml of 10 mM MgSO4, 10 mM MgCl2, 10 mM MES-KOH, pH 5.8, buffer, 0.5 M mannitol (MMM550) which was carefully layered on top on 8 ml 0.6 M sucrose cushion and spun down at 100×g, 2 min, at room temperature. Intact protoplasts at the interface were collected and spun down at 100×g for 2 mins then re-suspended in MMM550 –– if used for alginate imbedding this step was repeated three times.

For viability test the protoplast-containing pellet was re-suspended in 5 ml 2.5 mM MES-KOH, pH 5.7, 125 mM CaCl2, 154 mM NaCl, 5 mM KCl, 0.5 mM glucose (W5), centrifuged for 5 min at 100×g, the supernatant poured off, and the pellet re suspended in 0.5 ml W5 and placed on ice until FACS, which was initiated immediately after the wash step.

Embedment of GFP-fluorescent protoplast in alginate

Protoplast embedment in alginate was essentially done as described in [27] except the thin alginate layer was formed using the ‘droplet on Ca-agar’ method as described in [65]. Briefly, protoplasts re-suspended in 200 μl MMM550 were mixed with 200 μl alginate solution (MMM550 + 2.8% alginate (low viscosity)). A 300 μl droplet was left on a Ca-Agar plate (0.4 M mannitol, 50 mM CaCl2, 1% plant agar (Duchefa 1001.5000)) which was tilted to spread out the droplet, and after 30 min a floating solution (0.4 M mannitol, 50 mM CaCl2) was added to the plates to allow for movement of the layer. The layer was taken up by a spatula and moved to small Petri dishes containing F-PCN (described in [8]) .

gDNA extraction

A single fully infiltrated leaf was thoroughly ground in liquid nitrogen and DNA was extracted using DNeasy Plant Mini Kit (Qiagen).

PCR of genome target NbRRA and NbPDS loci

PCR-amplicons containing the NbRRA & NbPDS targets were amplified using nested PCR: First 5 μl of protoplast suspension (obtained as described in ‘Protoplast isolation’) was used in a 50 μl PCR reaction using Phire Plant Direct PCR Master Mix (ThermoFisher F160S) with the cycle parameters: 5 min at 98 °C followed by 40 cycles of 10 s at 98 °C, 10 s at (65 °C for RRA and 62 °C for PDS) and 40 s at 72 °C followed by 7 min at 72 °C using the primers P348 & P232 and P346 & P342 for NbPDS2 and NbRRAall1, respectively. Nested NbRRA PCR was performed in a 50 μl reaction using X7 polymerase [66] with 1:100 diluted 1’th PCR reaction as template and the cycle parameters: 5 min at 94 °C followed by 25 cycles of 30 s at 94 °C, 30 s at 58 °C and 30 s at 72 °C followed by 7 min at 72 °C and the primers P319 and P320.

Nested PCR of NbPDS2 was done in a 25 μl reaction using ClonAMP HiFi master mix 2x (Takara 639,298) with the cycle parameters: 5 min at 98 °C followed by 20 cycles of 30 s at 98 °C, 30 s at 65 °C and temperature dropping 0.5 °C per cycle and 30 s at 72 °C followed by 20 cycles of 30 s at 98 °C, 30 s at 58 °C and 30 s at 72 °C followed by 3 min at 72 °C and the primers P321 and P322.

Primers for scoring in leaf mutations were P321 & P322 (NbRRAall1) and P232 & P233 (NbPDS2).

Cloning in pJet and sequencing

10 μl of PCR product was digested ON in a 50 μl reaction with BtgI (NbRRAall1 amplicon) and AvrII (NbPDS amplicon). Enzyme Resistant bands were isolated from agarose gels using NucleoSpin® Gel and Monarch® DNA Gel Extraction Kit (New England Biolabs) and cloned into pJet1.2 using CloneJET PCR Cloning Kit #K1232. Sequences were aligned using CLC Workbench.

Indel detection by amplicon analysis (IDAA) and semi quantification of IDAA peaks

Indel Detection by Amplicon Analysis (IDAA) was done essentially as described in and outlined in [39] and in the Method section ‘PCR of genome target NbRRA and NbPDS loci’. Briefly A tri-primer PCR setup which relies of the incorporation of a florescent universal 6-FAM 5′-labelled primer (FamF), with the corresponding non-labeled primer in a 1:10 diluted concentration, was used for FAM labeling of PCR amplicons. PCR amplification of the NbRRAall1 and NbPDS2 regions were done using the ClonAMP HiFi master mix 2x (Takara 639,298) in a 25 μl reaction with the cycle parameters: 5 min at 95 °C followed by 30 cycles of 30 s at 95 °C, 30 s at 58 °C and 30 s at 72 °C followed by 3 min at 72 °C. Primers were P230 & P231 (NbRRAall1) and P232 & P233 (NbPDS2), where bold designates FAM primer overhang (Additional file 5: Table S1).

Mutation frequencies, as identified by quantification of peak area in IDAA chromatograms, were estimated using the Open Source Software program ImageJ ( from and with areas identified as described (

Fluorescence microscopy

Fluorescence imaging (presence GFP) was carried out with an epifluorescence microscope Olympus BX41 equipped with a CCD camera (FITC filter for GFP fluorescence and DAPI filter for FDA staining) or a laser scanning confocal microscope Leica SP5 equipped with an Argon (448 nm) and a Argon laser (448 nm).

Western blot analysis

App. 50 μl of seedling powder, crushed in liquid N2, was boiled in 50 μl 2 × SDS-PAGE loading buffer (280 mM SDS, 400 mM Tris, 40% glycerol, 1.4 M mercaptoethanol, 0.6 mM Bromophenol Blue) for 15 min and separated (200 V, 50 min) on 12% Criterion XT Bis-Tris gels (Bio-rad). Proteins were electrotransferred onto polyvinylidene difluoride (PVDF) membranes (Bio-rad) using a Trans-Blot® TurboTM Blotting instrument (Bio-rad). The membrane was blocked in blocking solution (PBS pH 7.5, 5% non-fat dry milk) overnight at 4 °C under mild shaking. The membrane was probed with anti-GFP mouse IgG (Roche) and Anti-Flag M2 mouse IgG (Sigma) at 11000 dilution in blocking solution overnight at 4 °C, followed by 3 × 5′ wash in PBS buffer (PBS pH 7.5). The membrane was then incubated with goat anti-mouse IgG conjugated to Alkaline Phosphatase (AP) (Sigma) (11000 dilution in blocking solution) for 1 h at room temperature, and rinsed 3 × 5′ with PBST. Pre-mixed NBT/BCIP AP solution (UCPH, DK) was added to the blot and incubated for color development.

Post FACS residual gRNA activity of lysed protoplasts

20 μl of protoplasts expressing gRNA-NbPDS2/SpCas9 were added to 80 μl PBS, briefly vortexed and left at room temperature for 2 h; 20 μl of protoplasts were added to PBS buffer with 5 μl RNAseA/T1 (Thermo fisher #EN0551), briefly vortexed and left at room temperature for 2 h; and 20 μl of protoplasts was flash frozen with immediately addition of 80 μl PBS which were then heated for 3 min at 95 °C. Flanking primers used were P233 and P232.

Flow cytometry and fluorescence activated cell sorting (FACS) of N. benthamiana protoplasts

The protoplast solution was first passed through 50-μm filcons (BD Biosciences) to achieve a single-cell suspension. Protoplast suspensions were cytometrically analyzed and sorted with a FACSAria III (BD Biosciences) fitted with a 100-μm nozzle and using phosphate-buffered saline (PBS) as a sheath fluid. The procedure and setting used were as described in [29] with a large nozzle size (100 μm) to provide optimal survival for most cell types and sorting based on ~ 10,000 events. Briefly, the sheath pressure was set at 20 psi, and the defection plate voltage was set at 5000 V (default “low” setting). A 488 nm Coherent Sapphire Solid State laser was used for excitation, and emission was measured at 530 nm for GFP. The photomultiplier tube voltage was set at 183 V for forward scatter, 286 V for side scatter, 308 V for GFP, and 518 V for Allophycocyanin. The threshold value for event detection was set at 8835 on forward scattering. The drop drive frequency was set to approximately 30 kHz, and the amplitude was set to approximately 45 V; the drop delay value was approximately 26 (these settings will vary slightly with day-to-day operation of the FACSAria III). Identification of viable, single protoplasts through the use of forward scatter (FSC) and side scatter (SSC) as a first gating strategy, which is routinely used for gating mammalian cells, was not attempted due to the high variability of protoplast size. Instead, the FITC and DAPI intensities were recorded as represented in dot plots. 10,000 events are displayed in each plot. Gates were set to separate and thus enable enrichment of WT and GFP transfected protoplasts, using the WT sample to define non-transfected wild type populations in the transfected samples. Transfected protoplasts were defined as FITC-positive events. Data were processed using the FACSDiva 8.0.1 software (BD Biosciences).

Viability test was done on a Sony SH800S Cell sorter, with automated setup for 130 μm microfluidics sorting chips, psi 9. PBS was used as sheath fluid, with samples sorted into flat-bottomed 96 well microtiter-plates containing 200 μl W5 buffer. For visualization purposes protoplasts were layered at the bottom of the microtiter plate by a brief centrifugation step, 100×g, 5 min.

Gating strategy on Sony SH800S cell sorter was similar to those used on the BD FACSAria III sorter.

Availability of data and materials

All constructs used in the present study are listed in Additional file 5: Table S1 and will be available upon request. Basic vector constructs (51491, 51493, 51494, 52255) were from and are available from Addgene (


A. tumefaciens :

Agrobacterium tumefaciens


Arabidopsis thaliana U6 promoter


Bean yellow dwarf virus


Cauliflower mosaic virus 35S promoter


Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated systems (Cas))


Fluorescence Activated Cell Sorting


Self-cleaving sequence of foot-and-mouth disease virus, 2A


guide RNA


Gemini Virus Replicon


Hygromycin phosphotransferase


Homologous Recombination


Indel Detection by Amplicon Analysis


Deletions and insertions

LB & RB:

Left and Right T-DNA border

N. benthamiana :

Nicotiana benthamiana


Phytoene desaturase


REDUCED RESIDUAL ARABINOSE arabinosyl transferase


Non-Homologous End Joining


Nucleus Localisation Signal


Nopaline synthase terminator


Protospacer adjacent motif


Precise genome editing


Restriction Enzyme


Streptococcus pyogenes Cas9


Transcription activator-like effector nucleases


Zinc Finger-Nucleases


  1. Tsai SQ, Joung JK. Defining and improving the genome-wide specificities of CRISPR-Cas9 nucleases. Nat Rev Genet. 2016;17(5):300–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Gao C. The future of CRISPR technologies in agriculture. Nat Rev Mol Cell Biol. 2018;19(5):275–6.

    Article  CAS  PubMed  Google Scholar 

  3. Ding Y, Li H, Chen LL, Xie K. Recent advances in genome editing using CRISPR/Cas9. Front Plant Sci. 2016;7:703.

    PubMed  PubMed Central  Google Scholar 

  4. Sanford JC, Smith FD, Russell JA. Optimizing the biolistic process for different biological applications. Methods Enzymol. 1993;217:483–509.

    Article  CAS  PubMed  Google Scholar 

  5. Rivera AL, Gomez-Lim M, Fernandez F, Loske AM. Physical methods for genetic plant transformation. Phys Life Rev. 2012;9(3):308–45.

    Article  PubMed  Google Scholar 

  6. Townsend JA, Wright DA, Winfrey RJ, Fu F, Maeder ML, Joung JK, Voytas DF. High-frequency modification of plant genes using engineered zinc-finger nucleases. Nature. 2009;459(7245):442–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Zhang Y, Zhang F, Li X, Baller JA, Qi Y, Starker CG, Bogdanove AJ, Voytas DF. Transcription activator-like effector nucleases enable efficient plant genome engineering. Plant Physiol. 2013;161(1):20–7.

    Article  CAS  PubMed  Google Scholar 

  8. Mao Y, Zhang H, Xu N, Zhang B, Gou F, Zhu JK. Application of the CRISPR-Cas system for efficient genome engineering in plants. Mol Plant. 2013;6(6):2008–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Li JF, Norville JE, Aach J, McCormack M, Zhang DD, Bush J, Church GM, Sheen J. Multiplex and homologous recombination-mediated genome editing in Arabidopsis and Nicotiana benthamiana using guide RNA and Cas9. Nat Biotechnol. 2013;31(8):688–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Woo JW, Kim J, Il Kwon S, Corvalan C, Cho SW, Kim H, Kim SG, Kim ST, Choe S, Kim JS. DNA-free genome editing in plants with preassembled CRISPR-Cas9 ribonucleoproteins. Nat Biotechnol. 2015;33(11):1162–U1156.

    Article  CAS  PubMed  Google Scholar 

  11. Li JF, Norville JE, Aach J, McCormack M, Zhang D, Bush J, Church GM, Sheen J. Multiplex and homologous recombination-mediated genome editing in Arabidopsis and Nicotiana benthamiana using guide RNA and Cas9. Nat Biotechnol. 2013;31(8):688–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Clasen BM, Stoddard TJ, Luo S, Demorest ZL, Li J, Cedrone F, Tibebu R, Davison S, Ray EE, Daulhac A, et al. Improving cold storage and processing traits in potato through targeted gene knockout. Plant Biotechnol J. 2016;14(1):169–76.

    Article  CAS  PubMed  Google Scholar 

  13. Gelvin SB. Agrobacterium-mediated plant transformation: the biology behind the “gene-jockeying” tool. Microbiol Mol Biol Rev. 2003;67(1):16–37 table of contents.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Nester EW. Agrobacterium: nature’s genetic engineer. Front Plant Sci. 2014;5:730.

    PubMed  Google Scholar 

  15. Gleba YY, Tuse D, Giritch A. Plant viral vectors for delivery by agrobacterium. Curr Top Microbiol Immunol. 2014;375:155–92.

    CAS  PubMed  Google Scholar 

  16. Peyret H, Lomonossoff GP. When plant virology met agrobacterium: the rise of the deconstructed clones. Plant Biotechnol J. 2015;13(8):1121–35.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Baltes NJ, Gil-Humanes J, Cermak T, Atkins PA, Voytas DF. DNA replicons for plant genome engineering. Plant Cell. 2014;26(1):151–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Cermak T, Baltes NJ, Cegan R, Zhang Y, Voytas DF. High-frequency, precise modification of the tomato genome. Genome Biol. 2015;16:232.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Gil-Humanes J, Wang Y, Liang Z, Shan Q, Ozuna CV, Sanchez-Leon S, Baltes NJ, Starker C, Barro F, Gao C, et al. High-efficiency gene targeting in hexaploid wheat using DNA replicons and CRISPR/Cas9. Plant J. 2017;89(6):1251–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Schneider JD, Marillonnet S, Castilho A, Gruber C, Werner S, Mach L, Klimyuk V, Mor TS, Steinkellner H. Oligomerization status influences subcellular deposition and glycosylation of recombinant butyrylcholinesterase in Nicotiana benthamiana. Plant Biotechnol J. 2014;12(7):832–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Alkanaimsh S, Karuppanan K, Guerrero A, Tu AM, Hashimoto B, Hwang MS, Phu ML, Arzola L, Lebrilla CB, Dandekar AM, et al. Transient expression of tetrameric recombinant human Butyrylcholinesterase in Nicotiana benthamiana. Front Plant Sci. 2016;7:743.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Dirnberger D, Steinkellner H, Abdennebi L, Remy JJ, van de Wiel D. Secretion of biologically active glycoforms of bovine follicle stimulating hormone in plants. Eur J Biochem. 2001;268(16):4570–9.

    Article  CAS  PubMed  Google Scholar 

  23. Le Mauff F, Mercier G, Chan P, Burel C, Vaudry D, Bardor M, Vezina LP, Couture M, Lerouge P, Landry N. Biochemical composition of haemagglutinin-based influenza virus-like particle vaccine produced by transient expression in tobacco plants. Plant Biotechnol J. 2015;13(5):717–25.

    Article  PubMed  CAS  Google Scholar 

  24. Dicker M, Tschofen M, Maresch D, Konig J, Juarez P, Orzaez D, Altmann F, Steinkellner H, Strasser R. Transient Glyco-engineering to produce recombinant IgA1 with defined N- and O-Glycans in plants. Front Plant Sci. 2016;7:18.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Yang Y, Li R, Qi M. In vivo analysis of plant promoters and transcription factors by agroinfiltration of tobacco leaves. Plant J. 2000;22(6):543–51.

    Article  CAS  PubMed  Google Scholar 

  26. Shen J, Fu J, Ma J, Wang X, Gao C, Zhuang C, Wan J, Jiang L. Isolation, culture, and transient transformation of plant protoplasts. Curr Protoc Cell Biol. 2014;63:2.8.1–17.

    Article  Google Scholar 

  27. Dovzhenko A, Bergen U, Koop HU. Thin-alginate-layer technique for protoplast culture of tobacco leaf protoplasts: shoot formation in less than two weeks. Protoplasma. 1998;204(1–2):114–8.

    Article  Google Scholar 

  28. Shepard JF, Totten RE. Isolation and regeneration of tobacco Mesophyll cell protoplasts under low osmotic conditions. Plant Physiol. 1975;55(4):689–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Lonowski LA, Narimatsu Y, Riaz A, Delay CE, Yang Z, Niola F, Duda K, Ober EA, Clausen H, Wandall HH, et al. Genome editing using FACS enrichment of nuclease-expressing cells and indel detection by amplicon analysis. Nat Protoc. 2017;12(3):581–603.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Bargmann BO, Birnbaum KD. Fluorescence activated cell sorting of plant protoplasts. J Vis Exp. 2010;(36):1673.

  31. Birnbaum K, Jung JW, Wang JY, Lambert GM, Hirst JA, Galbraith DW, Benfey PN. Cell type-specific expression profiting in plants via cell sorting of protoplasts from fluorescent reporter lines. Nat Methods. 2005;2(8):615–9.

    Article  CAS  PubMed  Google Scholar 

  32. Kirchhoff J, Raven N, Boes A, Roberts JL, Russell S, Treffenfeldt W, Fischer R, Schinkel H, Schiermeyer A, Schillberg S. Monoclonal tobacco cell lines with enhanced recombinant protein yields can be generated from heterogeneous cell suspension cultures by flow sorting. Plant Biotechnol J. 2012;10(8):936–44.

    Article  CAS  PubMed  Google Scholar 

  33. Carqueijeiro I, Guimaraes AL, Bettencourt S, Martinez-Cortes T, Guedes JG, Gardner R, Lopes T, Andrade C, Bispo C, Martins NP, et al. Isolation of cells specialized in anticancer alkaloid metabolism by fluorescence-activated cell sorting. Plant Physiol. 2016;171(4):2371–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Nishimasu H, Ran FA, Hsu PD, Konermann S, Shehata SI, Dohmae N, Ishitani R, Zhang F, Nureki O. Crystal structure of Cas9 in complex with guide RNA and target DNA. Cell. 2014;156(5):935–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Waibel F, Filipowicz W. U6 snRNA genes of Arabidopsis are transcribed by RNA polymerase III but contain the same two upstream promoter elements as RNA polymerase II-transcribed U-snRNA genes. Nucleic Acids Res. 1990;18(12):3451–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Li X, Jiang D-H, Yong K, Zhang D-B. Varied transcriptional efficiencies of multiple Arabidopsis U6 small nuclear RNA genes. J Integr Plant Biol. 2007;49(2):222–9.

    Article  CAS  Google Scholar 

  37. El Amrani A, Barakate A, Askari BM, Li X, Roberts AG, Ryan MD, Halpin C. Coordinate expression and independent subcellular targeting of multiple proteins from a single transgene. Plant Physiol. 2004;135(1):16–24.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Szymczak AL, Workman CJ, Wang Y, Vignali KM, Dilioglou S, Vanin EF, Vignali DA. Correction of multi-gene deficiency in vivo using a single ‘self-cleaving’ 2A peptide-based retroviral vector. Nat Biotechnol. 2004;22(5):589–94.

    Article  CAS  PubMed  Google Scholar 

  39. Yang Z, Steentoft C, Hauge C, Hansen L, Thomsen AL, Niola F, Vester-Christensen MB, Frodin M, Clausen H, Wandall HH, et al. Fast and sensitive detection of indels induced by precise gene targeting. Nucleic Acids Res. 2015;43(9):e59.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Vasil V, Vasil IK. Regeneration of tobacco and petunia plants from protoplasts and culture of corn protoplasts. In Vitro. 1974;10:83–96.

    Article  CAS  PubMed  Google Scholar 

  41. Egelund J, Obel N, Ulvskov P, Geshi N, Pauly M, Bacic A, Petersen BL. Molecular characterization of two Arabidopsis thaliana glycosyltransferase mutants, rra1 and rra2, which have a reduced residual arabinose content in a polymer tightly associated with the cellulosic wall residue. Plant Mol Biol. 2007;64(4):439–51.

    Article  CAS  PubMed  Google Scholar 

  42. Velasquez SM, Ricardi MM, Dorosz JG, Fernandez PV, Nadra AD, Pol-Fachin L, Egelund J, Gille S, Harholt J, Ciancia M, et al. O-glycosylated cell wall proteins are essential in root hair growth. Science. 2011;332(6036):1401–3.

    Article  CAS  PubMed  Google Scholar 

  43. Doench JG, Hartenian E, Graham DB, Tothova Z, Hegde M, Smith I, Sullender M, Ebert BL, Xavier RJ, Root DE. Rational design of highly active sgRNAs for CRISPR-Cas9-mediated gene inactivation. Nat Biotechnol. 2014;32(12):1262–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Garneau JE, Dupuis ME, Villion M, Romero DA, Barrangou R, Boyaval P, Fremaux C, Horvath P, Magadan AH, Moineau S. The CRISPR/Cas bacterial immune system cleaves bacteriophage and plasmid DNA. Nature. 2010;468(7320):67–71.

    Article  CAS  PubMed  Google Scholar 

  45. Gao Y, Zhang Y, Zhang D, Dai X, Estelle M, Zhao Y. Auxin binding protein 1 (ABP1) is not required for either auxin signaling or Arabidopsis development. Proc Natl Acad Sci. 2015;112(7):2275–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Regnard GL, Halley-Stott RP, Tanzer FL, Hitzeroth II, Rybicki EP. High level protein expression in plants through the use of a novel autonomously replicating geminivirus shuttle vector. Plant Biotechnol J. 2010;8(1):38–46.

    Article  CAS  PubMed  Google Scholar 

  47. Hefferon KL, Dugdale B. Independent expression of rep and RepA and their roles in regulating bean yellow dwarf virus replication. J Gen Virol. 2003;84(12):3465–72.

    Article  CAS  PubMed  Google Scholar 

  48. Jinek M, Chylinski K, Fonfara I, Hauer M, Doudna JA, Charpentier E. A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. Science. 2012;337(6096):816–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Siemering KR, Golbik R, Sever R, Haseloff J. Mutations that suppress the thermosensitivity of green fluorescent protein. Curr Biol. 1996;6(12):1653–63.

    Article  CAS  PubMed  Google Scholar 

  50. Chiu W, Niwa Y, Zeng W, Hirano T, Kobayashi H, Sheen J. Engineered GFP as a vital reporter in plants. Curr Biol. 1996;6(3):325–30.

    Article  CAS  PubMed  Google Scholar 

  51. Takebe I, Otsuki Y, Aoki S. Isolation of tobacco mesophyll cells in intact and active state. Plant Cell Physiol. 1968;9(1):115–24.

    Google Scholar 

  52. Feng Z, Mao Y, Xu N, Zhang B, Wei P, Yang DL, Wang Z, Zhang Z, Zheng R, Yang L, et al. Multigeneration analysis reveals the inheritance, specificity, and patterns of CRISPR/Cas-induced gene modifications in Arabidopsis. Proc Natl Acad Sci U S A. 2014;111(12):4632–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Liu Y, Tao W, Wen S, Li Z, Yang A, Deng Z, Sun Y. In vitro CRISPR/Cas9 system for efficient targeted DNA editing. mBio. 2015;6(6):e01714–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Tzfira T, Yarnitzky O, Vainstein A, Altman A. Agrobacterium rhizogenes-mediated DNA transfer inPinus halepensis mill. Plant Cell Rep. 1996;16(1–2):26–31.

    Article  CAS  PubMed  Google Scholar 

  55. Ziemienowicz A, Shim YS, Matsuoka A, Eudes F, Kovalchuk I. A novel method of transgene delivery into triticale plants using the agrobacterium transferred DNA-derived nano-complex. Plant Physiol. 2012;158(4):1503–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Yin K, Han T, Liu G, Chen T, Wang Y, Yu AYL, Liu Y. A geminivirus-based guide RNA delivery system for CRISPR/Cas9 mediated plant genome editing. Sci Rep-Uk. 2015;5:14926.

    Article  CAS  Google Scholar 

  57. Bashandy H, Jalkanen S, Teeri TH. Within leaf variation is the largest source of variation in agroinfiltration of Nicotiana benthamiana. Plant Methods. 2015;11:47.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. Mysore KS, Bassuner B, Deng XB, Darbinian NS, Motchoulski A, Ream W, Gelvin SB. Role of the agrobacterium tumefaciens VirD2 protein in T-DNA transfer and integration. Mol Plant-Microbe Interact. 1998;11(7):668–83.

    Article  CAS  PubMed  Google Scholar 

  59. Anand A, Krichevsky A, Schornack S, Lahaye T, Tzfira T, Tang Y, Citovsky V, Mysore KS. Arabidopsis VIRE2 INTERACTING PROTEIN2 is required for agrobacterium T-DNA integration in plants. Plant Cell. 2007;19(5):1695–708.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Waltz E. Tiptoeing around transgenics. Nat Biotechnol. 2012;30(3):215–7.

    Article  CAS  PubMed  Google Scholar 

  61. Bombarely A, Rosli HG, Vrebalov J, Moffett P, Mueller LA, Martin GB. A draft genome sequence of Nicotiana benthamiana to enhance molecular plant-microbe biology research. Mol Plant-Microbe Interact. 2012;25(12):1523–30.

    Article  CAS  PubMed  Google Scholar 

  62. Qin GJ, Gu HY, Ma LG, Peng YB, Deng XW, Chen ZL, Qu LJ. Disruption of phytoene desaturase gene results in albino and dwarf phenotypes in Arabidopsis by impairing chlorophyll, carotenoid, and gibberellin biosynthesis. Cell Res. 2007;17(5):471–82.

    Article  CAS  PubMed  Google Scholar 

  63. Voinnet O, Rivas S, Mestre P, Baulcombe D. An enhanced transient expression system in plants based on suppression of gene silencing by the p19 protein of tomato bushy stunt virus. Plant J. 2003;33(5):949–56.

    Article  CAS  PubMed  Google Scholar 

  64. Sainsbury F, Lomonossoff GP. Extremely high-level and rapid transient protein production in plants without the use of viral replication. Plant Physiol. 2008;148(3):1212–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Maćkowska K, Jarosz A, Grzebelus E. Plant regeneration from leaf-derived protoplasts within the Daucus genus: effect of different conditions in alginate embedding and phytosulfokine application. Plant Cell, Tissue Organ Cult (PCTOC). 2014;117(2):241–52.

    Article  CAS  Google Scholar 

  66. Nørholm MH. A mutant Pfu DNA polymerase designed for advanced uracil-excision DNA engineering. BMC Biotechnol. 2010;10(1):21.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references


Nicholas Baltes, Minnesota University, U.S.A, is acknowledged for providing the Gemini Virus Replicon vector constructs. Jen Sheen, Harvard Medical School, U.S.A, is acknowledged for providing the HBT-Cas9 construct.


This work was supported by The Danish Councils for Strategic and Independent Research (12–125709, 12–131859) (PGE for controlled yeast and plant and glycosylation), The Danish National Research Foundation (DNRF107), The Copenhagen University Excellence Program for Interdisciplinary Research (CDO2016) PGE for controlled plant glycosylation, Villum Foundation (00017489) (Confocal imaging and microscopy).

Author information

Authors and Affiliations



All authors have seen and approved the manuscript. RSM, BLP, JM, MC, YL and ZY conducted the experiments and BLP, RSM, BJ, HHW, EPB and ZY designed the experiments. BLP wrote the paper. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Bent Larsen Petersen.

Ethics declarations

Ethics approval and consent to participate

Nothing to declare.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Localized GFP-fluorescence of SpCas9-2A-GFP/NbPDS2-gRNA. GFP-fluorescence of SpCas9-2A-GFP/NbPDS2-gRNA from Agrobacterium infiltrated in leaves of N. benthamiana 3 days post infiltration was evident in the contours epidermis cells of intact leaves (A). Overlay of bright field and fluorescence (FITC filter) microscopy of isolated SpCas9-2A-GFP/NbPDS2-gRNA transformed protoplasts regularly showed 60 - > 80% transformation efficiency (B). GFP fluorescence was seen in cytoplasmic strands with some nuclei accumulation (A and C), which both are in accordance with a primarily cytoplasmic localization of the GFP. (PDF 9243 kb)

Additional file 2:

Figure S2.WT N. benthamiana protoplasts pre and post FACS. Protoplasts were isolated as described in the Methods section, stored in buffer MMM550 on ice and immediately FACS sorted (total population sorted) into the MMM550 buffer and stored on ice. An estimated survival rate of ca 10–20% (concentric intact protoplast) was observed as evidenced by bright field (A) microscopy. (B, C) GFP expression analysis using confocal microscopy and viability test using propidium iodide (PI). Left panels are scan of the protoplasts expressing SpCas9-2A-GFP construct. Right panels non-transformed control. (B) Distinguishable GFP signal can be observed in transformed protoplasts (arrowhead). (C) Viability analysis using propidium iodide (PI). Non PI stained protoplast expressing GFP were observed. (PDF 7907 kb)

Additional file 3:

Figure S3. Post FACS residual gRNA/Cas9 activity of lysed protoplasts. Ribonucleoprotein, i.e. in vitro transcribed gRNA mixed with heterologous expressed SpCas9 enzyme, delivered by PEG transformation, have been shown to confer efficient nuclease activity in Arabidopsis thaliana, tobacco, lettuce and rice protoplasts [10, 53]. We tested whether PBS mediated protoplast lysis could mediate additional extra-cellular derived indel formation resulting in an over-estimated gRNA/SpCas9 activity. Incubation 2 h at room temperature in PBS buffer resulted in a 2–3 fold increased indel formation, compared to immediate activity abolishment through flash freezing/boiling or RNAse addition, as judged by resistant RE band intensities. Lanes: Pos Ctrl (NbPDS2-gRNA/SpCas9 positive from leaves), Neg ctrl (WT without NbPDS2-gRNA/SpCas9), RT 2 h (PBS mediated lysis followed by 2 h incubation at room temperature), flash freezing (flash freezing in liquid N2 followed by boiling), +RNAase (RNAase addition). For experimental setup see Method section (PDF 191 kb)

Additional file 4:

Figure S4. GFP-fluorescent protoplasts embedded in alginate. Single fluorescent protoplasts are visible as evidenced by fluorescent (FITC filter) microscopy before (A) and after alginate embedment (B). Calli formation (C) of a single protoplast as evidenced by bright field microscopy. Protoplast embedment in alginate is described in the Methods section. (PDF 4381 kb)

Additional file 5:

Table S1. Vector construct and primer list (DOCX 18 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Petersen, B.L., Möller, S.R., Mravec, J. et al. Improved CRISPR/Cas9 gene editing by fluorescence activated cell sorting of green fluorescence protein tagged protoplasts. BMC Biotechnol 19, 36 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: