Next Article in Journal
Genome-Wide Investigation and Expression Analysis of the Catalase Gene Family in Oat Plants (Avena sativa L.)
Next Article in Special Issue
Studies on Improving the Efficiency of Somatic Embryogenesis in Grapevine (Vitis vinifera L.) and Optimising Ethyl Methanesulfonate Treatment for Mutation Induction
Previous Article in Journal
Impact of Chitosan-Based Foliar Application on the Phytochemical Content and the Antioxidant Activity in Hemp (Cannabis sativa L.) Inflorescences
Previous Article in Special Issue
Putative Daucus carota Capsanthin-Capsorubin Synthase (DcCCS) Possesses Lycopene β-Cyclase Activity, Boosts Carotenoid Levels, and Increases Salt Tolerance in Heterologous Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Apple CRISPR-Cas9—A Recipe for Successful Targeting of AGAMOUS-like Genes in Domestic Apple

1
Department of Biology, University of Colorado Colorado Springs, Colorado Springs, CO 80918, USA
2
U.S. Department of Agriculture, Agricultural Research Service (USDA-ARS), The Appalachian Fruit Research Station, 2217 Wiltshire Road, Kearneysville, WV 25430, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Plants 2023, 12(21), 3693; https://doi.org/10.3390/plants12213693
Submission received: 18 September 2023 / Revised: 19 October 2023 / Accepted: 24 October 2023 / Published: 26 October 2023

Abstract

:
Fruit trees and other fruiting hardwood perennials are economically valuable, and there is interest in developing improved varieties. Both conventional breeding and biotechnology approaches are being utilized towards the goal of developing advanced cultivars. Increased knowledge of the effectiveness and efficiency of biotechnology approaches can help guide use of the CRISPR gene-editing technology. Here, we examined CRISPR-Cas9-directed genome editing in the valuable commodity fruit tree Malus x domestica (domestic apple). We transformed two cultivars with dual CRISPR-Cas9 constructs designed to target two AGAMOUS-like genes simultaneously. The main goal was to determine the effectiveness of this approach for achieving target gene changes. We obtained 6 Cas9 control and 38 independent CRISPR-Cas9 events. Of the 38 CRISPR-Cas9 events, 34 (89%) had gene edits and 14 (37%) showed changes to all alleles of both target genes. The most common change was large deletions, which were present in 59% of all changed alleles, followed by small deletions (21%), small insertions (12%), and a combination of small insertions and deletions (8%). Overall, a high rate of successful gene alterations was found. Many of these changes are predicted to cause frameshifts and alterations to the predicted peptides. Future work will include monitoring the floral development and floral form.

Graphical Abstract

1. Introduction

Fruit crops are important sources of food and nutrition for humans. A 2022 estimate of global fruit production was 890 million tons worldwide, with the United States alone accounting for 470 million tons [1]. Major fruit crops include Solanum lycopersicum (tomato), Musa spp. (bananas and plantains), Citrullus lanatus (watermelon), Malus x domestica (domestic apple), and Citrus spp. (citrus) [1]. Given increased consumer demand, both conventional breeding and genetic engineering (GE) approaches are being used to obtain fruit cultivars with altered traits to better address some of the needs of these crops and to meet consumer demand [2,3,4]. Desired changes include improved disease resistance, abiotic stress tolerance, increased nutritional value, and reduced post-harvest waste [5,6,7,8]. Both conventional and GE approaches can sometimes lead to the same desired outcome. Such was the case for red-fleshed domestic apples, obtained either from breeding with red-fleshed cultivars [9] or from overexpression of the apple transcription factor MYB10 for the anthocyanin biosynthesis pathway [10]. Similarly, purple-fruited tomatoes, also with higher anthocyanin content, have been achieved either from trait introgression from a wild relative [11] or by creating a transgenic variety with genes from Antirrhinum (snapdragon), a non-breeding-compatible species [12]. Sometimes, a novel fruit coloration trait achieved via genetic engineering is distinct enough for marketing, as is the case of the pink-fruited Ananas comosus (pineapple) recently released by Del Monte Fresh Produce [13].
CRISPR-Cas9 genome editing has been successfully applied to a wide range of agricultural, ornamental, and horticultural crops [14,15].This approach shows great promise for achieving desired changes in fruit crops of interest. Numerous proof-of-concept studies have shown that it is possible to achieve stable CRISPR-Cas9 genome edits in various fruit and nut species, including apple, banana, Vaccinium corymbosum (blueberry), Theobroma cacao (cacao), Castanea sativa (sweet chestnut), citrus, Coffea canephora (coffee), Vitis vinifera (grape), Pyrus communis (pear), Actinidia deliciosa (kiwi), and Punica granatum (pomegranate) (Table 1, [16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31]). There is interest in developing approaches in many other fruit species, including tropical fruits, but progress is hampered by a need for improved genetic resources, micropropagation, gene transfer abilities, and more [32]. A recent study in apple found that editing the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 6 (MdSPL6) gene led to increased shoot formation in tissue culture; perhaps this approach could be applied to challenging species or cultivars [33]. For most proof-of-concept studies, the gene PHYTOENE DESATURASE (PDS) was targeted, as the inactivation of gene function leads to the development of albino cells, which can be readily visualized during the regeneration process. However, alterations in most genes of interest do not lead to phenotype changes that can be visually assessed in tissue culture. Therefore, some initial studies of CRISPR-Cas9 gene editing of non-PDS targets assessed genetic changes to the targeted gene or genes (Table 1).
Many applied studies of CRISPR-Cas9 gene editing in fruit crops are aimed at achieving improved disease resistance. This biotic stress is a major cause of reductions in crop yield [34,35,36]. Improved resistance to both viral and fungal diseases has been achieved in banana, citrus, grape, and cacao [37,38,39,40,41,42,43,44]. Other traits of interest are related to accelerated breeding and crop production. Early-flowering varieties of apple and pear and were achieved via editing of TFL1, which should shorten the breeding cycle [16]. Bananas with a semi-dwarf form were obtained by simultaneous editing of sites in a gibberellin pathway gene present in multiple copies, which should both assist in ease of harvesting and reduce losses from lodging [45,46]. Bananas with longer fruit shelf life were obtained from edits to a gene in the ethylene synthesis pathway [47]. Other crop-specific changes have been achieved as well, such as grapes with decreased tartaric acid content and apples with lowered foliar phloridzin [48,49].
Table 1. Examples of CRISPR-Cas9 gene editing in perennial fruit and nut species.
Table 1. Examples of CRISPR-Cas9 gene editing in perennial fruit and nut species.
Crop SpeciesGene(s) TargetedOutcomeReferences
Proof-of-concept studies
applePDSAlbino and sectored explants[22]
applePDSAlbino and sectored explants[16]
applePDSTarget gene changes[23]
appleDIPM-1, DIPM-2, DIPM-4Target gene changes[19]
bananaPDS1, PDS2Albino and variegated explants[18]
bananaPDSAlbino and pale green explants[20]
blueberryPDSAlbino and sectored explants[26]
chestnutPDSAlbino explants[24]
citrus—Carrizo citrangePDSAlbino and sectored explants[28]
citrus—grapefruitPDSTarget gene changes[17]
citrus—mini-citrusPDSAlbino and sectored explants[29]
citrus—mini-citrusCCD4Target gene changes[29]
coffeePDSAlbino explants[30]
grapeIdnDHTarget-site changes[23]
grapePDSVariegated explants[21]
grapePDSAlbino explants[25]
grapeMLO-7Target gene changes[19]
kiwiPDSAlbino explants[27]
pomegranatePgUGTsTarget-site changes[31]
Applied work—improved disease resistance
appleMdDIPM4Resistance to fire blight[40]
bananaeBSV (viral gene)Resistance to banana streak virus[41]
cacaoTcNPR3Resistance to Phytophthora tropicalis[37]
citrus—grapefruitCsLOB1Resistance to citrus canker[38]
citrus—Wanjincheng orangeCsLOB1 (promoter)Resistance to citrus canker[39]
citrus—Wanjincheng orangeCsWRKY22Resistance to citrus canker[43]
grapeVvMLO3Resistance to powdery mildew[42]
grapeVvWRKY52Resistance to Botrytis cinerea[44]
Applied work—shortened juvenile period
appleTFL1Early flowering[16]
pearTFL1Early flowering[16]
Applied work—plant form
bananaMaGA20ox2Semi-dwarf size[46]
Applied work—fruit shelf life
bananaMaACO1Increased shelf life[47]
Applied work—altered metabolite content
appleMdPGT1Reduced phloridzin levels[49]
grapeIdnDHReduced tartaric acid levels[48]
Applied work—increased shoot regeneration
appleMdSPL6Increased shoot regeneration[33]
Apple is a valuable fruit crop, with yearly global apple production in 2021 reaching 9.4 million tons [1]. It would be advantageous to improve resistance to key biotic stresses, including fire blight and apple scab, as these factors severely impact tree health and decrease marketable fruit yield. There are resistant wild relatives available for trait introgression into domestic apple [50,51]. However, the traditional breeding approach of using compatible relatives to add traits is challenging in apple, due in part to the long juvenile period, as well as the lack of self-compatibility and the genetic heterozygosity of apple cultivars in general [52]. Apple is amenable to a variety of biotechnology methods, and these approaches have been used to achieve valuable traits. Early-flowering apple, which could be used to accelerate the breeding cycle, has been obtained through diverse approaches, including virus-induced expression of Arabidopsis thaliana (Arabidopsis) FLOWERING LOCUS T (AtFT) [53], overexpression of Betula pendula Roth (silver birch) BpMADS4 [54], stable overexpression of an apple FT homolog MdFT1 [55], RNAi of apple TERMINAL FLOWER 1 (TFL1) homologs [56], and CRISPR targeting of apple TFL1 [16]. The early-flowering BpMADS4 trees are being used in accelerated breeding programs for high-value trait introgression [57,58]. Improved rooting and reduced tree size occurred with expression of the gene ROOTING LOCUS B (rolB) from Agrobacteria rhizogenes in rootstocks grafted to non-transgenic scions [59]. Overexpression of the Lc gene from Zea maize (corn) led to improved resistance to both apple scab and fire blight, but tree form was altered in commercially unacceptable ways [60]. Improved resistance to fire blight has also been achieved by CRISPR-mediated genome editing of the MdDIPM4 gene in apple; however, this study did not report the overall vegetative form [40]. There is a precedent for applying these studies to apple production, as one variety of genetically engineered apple with non-browning fruits obtained from RNAi suppression of POLYPHENOL OXIDASE (PPO) has been commercialized [61]. Therefore, additional knowledge of the efficacy of biotechnology approaches in apple could be applicable to cultivar improvement and commercial production.
A goal of our work was to determine the efficiency and efficacy of dual-gene editing with CRISPR-Cas9 in domestic apple. Prior work with CRISPR gene targeting in apple demonstrated that individual genes can be successfully targeted and edited with gene-specific single-guide RNAs (sgRNAs) [16,19,22,23]. Apple has a diploid genome; however, there is evidence of a historical whole-genome duplication event, and many gene families are greatly duplicated [62]. It would be advantageous to perform editing of two or more genes at the same time, making use of conserved gene regions. For this project, we used CRISPR-Cas9 to target two apple homologs of the key floral development gene AGAMOUS (AG). We previously identified MdMADS15 (hereafter MADS15) and MdMADS221 (hereafter MADS221) as close relatives of AG [63]. Our prior work utilizing RNA interference (RNAi) to suppress MADS15 and MADS221 in the ‘Galaxy’ cultivar reduced both male and female fertility, reduced floral determinacy, and led to extra petals through the homeotic conversion of stamens to petals [63]. An applied goal of this prior study was to reduce fertility and achieve decreased invasive potential while retaining blooms. However, even strongly suppressed RNAi events still showed some target gene expression, as well as reduced expression of the non-target gene MdMADS14. Additionally, the ‘Galaxy’ cultivar previously used for RNAi has the unusual apple trait of fruit set in the absence of pollination, which was an unexpected outcome. Our current CRISPR-Cas9 targeting project utilized a more precise gene-targeting method and two cultivars of domestic apple, ‘Royal Gala’ and ‘M.26’, that should only set fruit after successful pollination and seed set. This work will add to the knowledge of the targeting efficiency of editing two genes simultaneously in a tree fruit crop and could be applied to reduce the invasive potential of flowering ornamental fruiting species such as Pyrus calleryana (Bradford pear), a tree initially popular as an ornamental planting and a source of breeding for fire-blight resistance but now recognized as an invasive species in several states [64].

2. Results

2.1. Analysis of Transformation Efficacy

Prior analysis of the MADS family from domestic apple indicated that two genes, MADS15 and MADS221, appear to be homologous to AG [63], with two alleles per gene. As CRISPR-Cas9 gene targeting relies on precise matches between sgRNAs and target sites, we amplified and sequenced portions of the MADS15 and MADS221 genes from both ‘Royal Gala’ and ‘M.26’ cultivars. Primers were designed based on the published domestic apple genome [62].
We designed a set of sgRNAs to target our two genes simultaneously (Supplementary Figure S1). The sgRNA sites were located in the first exon of each gene to hopefully disrupt protein coding (Figure 1). We created four different sgRNA constructs, each with a different pairing of sgRNAs (pairings and sgRNA sequences are shown in Supplementary Figure S1). We used two sgRNAs in each gene-editing construct to increase target gene editing frequency and potentially create larger deletions (Figure 1). We created a Cas9 control construct to be used as a negative control. Plants transformed with this construct underwent transformation and regeneration with selection. As this construct lacks sgRNAs, no editing was predicted to occur at the target sites.
Four gene-editing sgRNA constructs and a Cas9 control construct were transformed into the ‘M.26’ and ‘Royal Gala’ cultivars via Agrobacterium-mediated genome insertion. We obtained a total of 44 transformed events (Table 2). Here, transformed events were defined as independently isolated shoots after transformation and growth on selective medium, and were confirmed by PCR to have the transgene construct present (workflow and shoot verification in Figure 2). Of these events, 40 were in ‘M.26’ and 4 were in ‘Royal Gala’. We observed 6 total Cas9 control events and 38 total CRISPR events, with the number of events per construct and cultivar ranging from 0 to 13.
We initially used dual-allele amplification of each target gene, meaning amplification directly followed by amplicon sequencing without subcloning or isolating alleles to determine if any events showed evidence of changes to target genes. Primers used are listed in Supplementary Table S1. Analysis of the sequence peaks and alignments with the wild-type sequence showed a variety of outcomes (Supplementary Figure S2). These included events with no changes to target sites (no changes in alignment and single sequence peaks), events with identical bi-allelic changes (changes in alignment and single sequence peaks, indicating identical alleles), and events with allelic differences (changes in alignment and double sequence peaks, indicating non-identical alleles). As these findings indicated there were changes to at least some target sites in some of our events, we decided to pursue allele-specific sequencing.
We then used dual-allele amplification of each target gene followed by subcloning of amplicons to separate alleles from each gene. Subcloned fragments were sent for DNA sequencing and the results aligned with our initial target-site sequences. Cloning of both alleles was verified by the number of amplicons present on the gel and by allele-specific single nucleotide polymorphisms (SNPs).
We observed six Cas9 control events, which had the Cas9 gene added but no sgRNAs to direct Cas9 to target sites. For these events, analysis of the MADS15 and MADS221 genes showed no changes to the target genes. Of our 38 CRISPR-Cas9 events, 34 showed changes to at least one allele (for one or more target genes), while 4 events had no changes to any sites, for an overall event-targeting rate of 89% (Table 3). We found that 14 events had changes to all four alleles (changes in both alleles of MADS15 and both alleles of MADS221), with a rate of 37% for all events. The remaining 20 events were a mix of WT and edited alleles. Complete information on each event and each allele can be found in Supplementary File S1.

2.2. Mutation Types and Frequency in Target Genes

We also analyzed the data in terms of editing by sgRNA. Three of the sgRNAs (sgRNA1, sgRNA2, sgRNA3) were designed to target both MADS15 and MADS221; they were termed dual-gene targeting, while two of the guide RNAs were specific to either MADS15 (sgRNA15) or MADS221 (sgRNA22, Supplementary Figure S1). Examination of targeting by guide RNA showed that target editing by sgRNA ranged from 37–68% of targets changed, with sgRNA3 giving the highest rate of editing (Table 4).
Each of the four editing constructs contained two sgRNAs (Figure 1). Examination of editing frequency by construct showed that one construct led to edits in all events obtained, while the other two constructs gave editing rates of 75–93% (Table 5).
Examination of mutation types showed that most changes were large deletions greater than 16 base pairs (bp); these accounted for 36–82% of edits observed from each construct (Table 6). Small deletions were also common (11–21%), as were small insertions (0–19%). All constructs led to a variety of mutation types, with both large and small deletions observed.
We examined both of our target genes for changes. Alignment of sequences showed a variety of mutations at a given target site for each gene. Small mutations tended to cluster around each sgRNA, while larger deletions spanned the region flanked by two sgRNA sites (Figure 3 and Figure 4).
We used the ExPASy peptide prediction program [65] followed by peptide alignment to illustrate possible peptides produced by our edited alleles. Some mutations were in-frame and led to the alternation of a portion of the predicted peptide, with most of the sequence intact. Other mutations led to predicted frame shifts and severe alterations in predicted peptide sequences, which often led to predicted early stop codons (Figure 5).
No obvious differences were noted between CRISPR-edited lines or within any line in terms of growth or overall appearance. Apple trees require several years to overcome juvenility, during which time they have no capacity to bloom. These trees are now entering maturity, as evidenced by extremely preliminary data: a single bloom from Line T532, which shows the expected phenotype of doubled petals and altered sexual organs (Supplementary Figure S3). In contrast, the M.26 parental cultivar is a single-flowered variety with well-formed stamens.

3. Discussion

Targeted genome editing in fruit crops such as apple has the potential to create high-value traits or remove unwanted features. For apples and other pomme fruits, some desirable traits could include reduced invasive potential, non-browning fruit, improved tolerance to abiotic and biotic stresses, and altered tree form. This sort of editing can also be highly useful in research, allowing for the creation of specific genetic knockouts. Here, we selected the gene targets MADS15 and MADS221, as they appear to be AG homologs in domestic apple. Prior work suppressing these genes via RNAi found that reduced AG-like gene function led to the conversion of anthers to petals and reduced seed set [63,66]. Thus, these genes may be useful targets for engineering floral sterility.
At the time this project was initiated in 2015, the efficiency of CRISPR-Cas9 gene editing in trees was still being discovered. Hence, we created double sgRNA constructs with redundant-gene targeting to try and increase the chances of edits to target sites. Additionally, using two sgRNAs would increase the likelihood of larger deletions. We also created a Cas9 control construct that lacked any guide RNAs to serve as our negative control, as suggested in CRISPR-Cas9 genome editing protocols [67]. We observed no target-site changes in our six Cas9 control events, as would be expected of a lack of genome editing in the absence of sgRNAs. All four editing constructs successfully led to the creation of events with changes to one or more of the target gene alleles. We found that the overall editing rate was high, with most events (34, 89%) showing edits to at least one target allele, and 14 events (37%) showed editing of all four alleles. These findings are similar to initial CRISPR-Cas9 gene-editing work in apple showing efficient target-site changes, with an editing rate of 31.8% for PDS [22] and a rate greater than 85% for PDS and TLF1 [16]. Some of this variation may be due to the different cultivars used, the gene or genes targeted, or to the regeneration and transformations methodologies used.
We found a variety of mutation types in our CRISPR-Cas9 events. Most mutations were large deletions spanning the region between the two sgRNAs, indicating activity at both sgRNA sites. This trend was observed for all four constructs tested. Some alleles had a mixture of small insertions and deletions. The small deletions and insertions we observed were similar to what was found previously for double sgRNA modification in domestic apple [16] and single sgRNA editing in domestic apple [22]. Our finding of many larger deletions spanning the sgRNA sites is different from these prior studies in apple, but is similar to work in Populus, where the use of two sgRNAs led to frequent larger deletions [68]. Unlike prior work using two sgRNAs in apple [16], we did not observe inversions at our target sites. Analysis of editing efficiency by construct and sgRNA showed that all four gene-editing constructs we developed were effective at achieving gene changes. As has been observed before in apple [22], sgRNAs were variable in their effectiveness, but all led to edits.
Analysis of predicted peptides for our events showed a range of predicted outcomes. Some DNA changes had minor impacts on the predicted peptide sequence, such as one deletion followed by one substitution, or the deletion of three amino acids. Other mutations that altered the reading frame led to early stop codons. Assessment of the functional impact of these changes will require waiting until most trees mature and flower. Currently, a single experimental tree has reached maturity. This early bloomer has given a preview of flowers from an event with changes to all four target alleles. Our trees have four total alleles and a variable number of mutated alleles and mutation types, ranging from no mutations to all four alleles changed. It will be interesting to see how this spectrum of mutations influences floral phenotypes.
We did not observe evidence of chimerism in the events characterized in this study, as we did not detect more than two alleles for each gene in any given event. Chimerism after organogenic regeneration is a concern for the study and the usage of CRISPR-Cas9-edited plants [69]. A prior study in apple found that 88% of PDS-targeted events and 100% of TFL1.1-targeted events were chimeric, as were 80% of TFL1.1-targeted events in pear [16]. Indeed, many initial studies targeting PDS in fruit crops observed chimeric explants, as evidenced by the presence of striped, variegated, or sectored regenerated plantlets (Table 1). However, in plants with PDS targeting, loss of gene function is detrimental to plant growth, with albino cells requiring nutritional support from the growth medium or from photosynthetic cells, and plantlets containing at least some green photosynthetic cells would have a growth and health advantage over their fully albino counterparts [70]. Thus, targeting PDS may give an overestimate of the general rate of chimerism that occurs in transformation and regeneration with gene editing. The final outcome of chimerism in a regenerated plant could arise from a genetically variable population of cells forming the final plantlet or could be due to ongoing CRISPR-Cas9 activity over time in the plantlet. An investigation of regenerated events derived from initially edited lines showed the presence of new genetic changes to target sites [69]. Thus, it appears that retained stably-integrated CRISPR editing constructs may have continued activity, or regenerated plants could arise from a subpopulation of cells from a chimeric parental plant. However, studies targeting tree genes other than PDS found a relatively low rate of chimerism. For example, a large-scale study of floral gene targeting in Populus characterized 591 events from a variety of CRISPR-Cas9 constructs and found a global chimerism rate of 1.3%, with most constructs giving 0% chimeric events [68]. Another study in Populus targeting seven tandemly-arrayed Nucleoredoxin1 (NRX1) alleles analyzed 42 individual events and uncovered a wide range of mutation types, including complex fusions and rearrangements, but no evidence of chimerism [71]. In eucalyptus, targeting of the key floral gene LEAFY showed no evidence of chimerism in the 68 events studied [72].
Overall, we found that CRISPR targeting of our two genes simultaneously led to a variety of edited outcomes, with some events fully WT and others edited at all alleles. Most mutations were large deletions, but a wide spectrum of mutation types was observed. As events varied in the number and severity of gene changes, they represent a highly variable population for future study. As our targeted genes are key for overall floral morphology, future work will hopefully include complete assessment of floral form, fertility, and fruits.

4. Materials and Methods

The draft version of the apple genome [62] was used to obtain sequences for MADS15 (Gene ID MDP0000324166, Gene symbol LOC103434929, GenBank cDNA Accession AJ251118, Phytozome Gene Identifier MD07G1022500) and MADS221 (Gene ID MDP0000250080, GenBank cDNA Accession HM122606, Phytozome Gene Identifier MD10G1271000). These data were used to design primers for amplifying and sequencing these genes from the ‘M.26’ and ‘Royal Gala’ cultivars that we selected for transformation and regeneration. Primers used for gene sequencing are found in Supplementary Table S1. The ‘M.26’ cultivar is commonly used as a dwarfing rootstock [73], while ‘Royal Gala’ is used for fruit production [74]. The cultivar-specific gene data were used to select guide RNA sites. Selected sites were 20 bp in length and adjacent to the sequence NRG, based on the cut-site preference for the Cas9 variant used.
  • Construct assembly
The Cas9-only control construct was created previously [68]. Double sgRNA constructs were created as per previously published methods [68]. In brief, two oligos were ordered for the assembly of each sgRNA (oligos for sgRNA construction are listed in Supplementary Table S1). Oligos were phosphorylated, annealed, and ligated into the vector psK-AtU626. Final assembled Cas9 sgRNA cassettes were ligated into the plant transformation vector pK2GW7.
  • Apple transformation and event validation
Clonally propagated apple (Malus x domestica) ‘M.26’ and ‘Royal Gala’ leaves underwent Agrobacterium-mediated transformation as described previously [75]. In brief, young leaves were harvested from axenic cultures, lightly wounded with non-traumatic forceps, inoculated with Agrobacterium tumefaciens strain EHA105 (pCH32) harboring the desired CRISPR constructs, and placed on co-cultivation media in the dark for three days. The leaves were then removed from the media, rinsed, cut into approximately 1 mm wide strips, and placed with the abaxial side down into regeneration media in the dark for three weeks. The regeneration plates were then placed in low light (Figure 2A,B) and regenerants (Figure 2C) transferred to shoot proliferation media (Figure 2D). The resulting regenerants were verified via PCR from tissue-culture plantlet leaves (Figure 2D) using the following primers: Malus x domestica elongation factor 1 (MdEF1, F5′-GACATTGCCCTGTGGAAGTT-3′, R 5′-GGTCTGACCATCCTTGGAAA-3′), used to assess competence of DNA for PCR; VirG (F 5′-GCCGGGGCGAGACCATAGG-3′, R 5′-CGCACGCGCAAGGCAACC-3′), used to assess whether there is residual, endophytic Agrobacterium contamination; and Cas9 and nos terminator specific primers to check for transgene presence (Cas9_end_F2: 5′-CCTACAACAAGCACCGGGAT-3′ and Cas9_tnos_R2: 5′-AACGATCGGGGAAATTCGAG-3′). Explants were propagated in tissue culture at the USDA-ARS-NEA-AFRS facility (Kearneysville, WV, USA), essentially as per published methods [76,77] with root induction [78]. Own-rooted young trees were transferred to soil and grown in a greenhouse with supplemental lighting (Na vapor) to maintain the day length at 16 h, and a maximum/minimum temperature range of 35/20 °C. Trees were watered daily and fertilized regularly (Supplementary Figure S4).
  • Target-site cloning and analysis
Initial amplicons of MADS221 and MADS15 genes were isolated from extracted DNA and amplified using a Q5 polymerase (New England BioLabs® Inc. M0491S, Ipswich, MA, USA) with gene-specific primers (MADS221: Forward 5′ ATGGCCAATGAAAACAAATCC 3′, Reverse 5′ CTGTTGTTGGCATACTCATAGAGG 3′; MADS15: Forward 5′ ATGGCCTATGAAAGCAAATCC 3′, Reverse 5′ CTATTGTTGGCATACTCATAGAGG 3′). PCR was performed with a SimpliAmp™ Thermal Cycler (Applied Biosystems A24812), and annealing temperature was set to 58 °C for 33 cycles. Sterile DI H2O was used for a negative control. The PCR target sequence was 230 bp for each gene. For subcloning, 4 µL of PCR product was set aside, and the remainder was combined with loading dye (40% sucrose, 30% glycerol, 2.0% Orange G, 10× concentration, used at 1X) for visualization on a 1% agarose gel with a 1 Kb DNA Ladder (New England BioLabs® Inc. N3232L, Ipswich, MA, USA) or 1 Kb Plus DNA Ladder (New England BioLabs® Inc. N3232L, Ipswich, MA, USA). Gels were visualized using a GelDoc Go Imaging System (BioRad Laboratories Inc.). Larger amplicons were obtained with the same method using primers MADS15 new F1 (5′ TTTCATTTGTTTCTGCAAGTTTC 3′) and new R1 (5′ CAAGTAAATTGAAGCAATAATATTACCTA 3′), and MADS221 new F2 (5′ GTTGGTGATCAAGAATATAGTAATTG 3′) and new R1 (5′ AGCAAGAAAATTGAAGCAATAATATTATTA 3′). Retained PCR products were subcloned into the vector Zero Blunt® TOPO® PCR Cloning Kit for Sequencing (Invitrogen™ 45-0031, Waltham, MA, USA) according to kit instructions with a 30 min initial incubation time. Cloning reactions were transformed into chemically competent E. coli cells (One Shot™ TOP10 Chemically Competent E. coli, Invitrogen™ C404010, Waltham, MA, USA) as per kit instructions. Transformed cells were divided among two LB + kanamycin (KAN) plates (1.0% tryptone, 0.5% yeast extract, 1.0% NaCl, 1.5% agar powder, 0.05% kanamycin) with 50 µL on one plate and the remainder (~300 µL) on the second plate. Plates were allowed to dry on the benchtop until all liquid was absorbed and were then incubated overnight at 37 °C (VWR™ 1515E, Radnor, PA, USA). Afterwards, the plates were wrapped in parafilm and stored at 4 °C for colony analysis.
  • Colony analysis
Eight colonies were selected from the transformed plates and resuspended in a 1.5 mL centrifuge tube with 100 µL of LB+KAN (1:1000). The colony broth was pipetted or vortexed (Fisher Analog Vortex 02215414, Waltham, MA, USA) and then shaken for 1 h at 37 °C. Colony PCR was performed using the colony broth with EconoTaq® DNA Polymerase (Lucigen 30031-1, Middleton, WI, USA) and the corresponding gene primers (MADS221 or MADS15) to allow for the amplification of cloned fragments. An annealing temperature of 56 °C was used for 33 cycles. Purified TOPO cloning vector with a cloned MADS15 fragment was used as a positive control, and sterile DI H2O was used as a negative control. Loading dye was added to the PCR product and visualized using a 1% agarose gel and 1 Kb DNA Ladder. Generally, 4 colonies with amplified fragments were selected from each sample for each gene for overnight liquid culture and plasmid extraction. Sample selection was based on relative band size of the product(s), and if colonies had different sized bands, then colonies with different sized products were selected for analysis. Next, 30 µL of the colony broth was used to inoculate a 5 mL culture of LB+KAN (1:1000), which was shaken at 37 °C overnight.
  • Plasmid Extraction and Digestion
Cells from the liquid culture were then pelleted using a Sorvall Legend Micro 21 Centrifuge (ThermoFisher Scientific, Waltham, MA, USA), and the DNA was extracted using a Monarch® Plasmid Miniprep Kit (New England BioLabs® Inc. T1010L/S, Ipswich, MA, USA) following kit instructions. DNA was quantified using a NanoDrop 2000 (ThermoFisher Scientific E112352, Waltham, MA, USA). Extracted DNA was then used for an enzyme digest with EcoRI (New England BioLabs® Inc. R0101S/L, Ipswich, MA, USA) to confirm successful gene fragment cloning, as the TOPO vector has EcoRI sites flanking the MCS, and no EcoRI sites are located in the target gene fragments. Solutions were incubated for at least 2 h at 37 °C and then visualized using a 1% agarose gel. Based on visualization results of the enzyme digest, plasmids were selected for sequencing.
  • DNA Sequencing and Peptide Prediction
The selected plasmid samples were mixed to a final DNA quantity of 400 ng in 10 µL with sterile DI H2O and sent out for Sanger Sequencing with GENEWIZ (Azenta Life Sciences, Burlington, WI, USA) using the M13F primer provided by Azenta. Sample sequences were compared with the WT gene form from the proper cultivar using the alignment tools available through the blastn suite [79]. Mutations were catalogued digitally based on size and type (i.e., 1 bp insertion, 15 bp deletion, etc.). Peptide predictions were made using the ExPASy Translate tool [65]. Alignments of predicted peptides were created using Clustal Omega [80].

5. Conclusions

CRISPR-Cas9 gene editing can be successfully applied to obtain dual targeting in domestic apple. The overall efficiency varies both by sgRNA and by gene target, with some trees retaining WT alleles. As the CRISPR-Cas9 system was stably transformed and retained in the trees, a long-term goal is surveying for continued Cas9 activity and allele targeting, along with analysis of vegetative performance and floral form.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12213693/s1, Figure S1: Gene sequences and guide RNA locations; Figure S2: Examples of bulk sequencing outputs; Figure S3: Floral phenotype of line T532; Figure S4: Examples of transformed plants in soil; Table S1: Primer names, sequences, and uses. File S1 raw sequence reads. File S2, predicted peptides.

Author Contributions

Conceptualization A.L.K.; formal analysis S.J., N.B., A.C., L.M., T.A.N., P.W., T.S.A. and A.L.K.; investigation S.J., N.B., A.C., L.M., T.A.N., P.W., T.S.A. and A.L.K.; resources, P.W., T.S.A. and A.L.K.; writing—original draft preparation S.J. and A.L.K.; writing—review and editing A.L.K. and T.S.A.; visualization S.J., T.A.N., P.W. and A.L.K.; supervision T.S.A. and A.L.K.; funding acquisition A.L.K. All authors have read and agreed to the published version of the manuscript.

Funding

Funding for this project was obtained from the College of Letters Arts and Sciences new lab startup funds to AL Klocko, the summer 2022 Undergraduate Research Academy to N Bondarchuk and AL Klocko, a 2022 COVID revitalization grant from UCCS to AL Klocko, and the Schmidt Family Charitable Foundation to AL Klocko. TA and PW were funded under USDA-ARS-CRIS project number 8080-00021-033-00D.

Data Availability Statement

The data presented in this study are available in the article or Supplementary Material.

Acknowledgments

We thank undergraduate research student Jen Yoon for assistance in target-site analysis. We also thank Brent Wallace of UCCS for ordering supplies.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. The findings and conclusions in this publication are those of the author(s) and should not be construed to represent any official USDA or U.S. Government determination or policy. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture. The U.S. Department of Agriculture prohibits discrimination in all its programs and activities on the basis of race, color, national origin, age, disability, and where applicable, sex, marital status, familial status, parental status, religion, sexual orientation, genetic information, political beliefs, reprisal, or because all or part of an individual’s income is derived from any public assistance program. (Not all prohibited bases apply to all programs.) Persons with disabilities who require alternative means for the communication of program information (Braille, large print, audiotape, etc.) should contact the USDA’s TARGET Center at (202) 720-2600 (voice and TDD). To file a complaint of discrimination, write to USDA Director, Office of Civil Rights, 1400 Independence Avenue, S.W., Washington, D.C. 20250-9410, or call (800) 795-3272 (voice) or (202) 720-6382 (TDD). USDA is an equal opportunity provider and employer.

References

  1. FAO. World Food and Agriculture—Statstical Yearbook; FAO: Rome, Italy, 2022. [Google Scholar]
  2. Lobato-Gomez, M.; Hewitt, S.; Capell, T.; Christou, P.; Dhingra, A.; Giron-Calva, P.S. Transgenic and genome-edited fruits: Background, constraints, benefits, and commercial opportunities. Hortic. Res. 2021, 8, 166. [Google Scholar] [CrossRef] [PubMed]
  3. Penna, S.; Jain, S.M. Fruit Crop Improvement with Genome Editing, In Vitro and Transgenic Approaches. Horticulturae 2023, 9, 58. [Google Scholar] [CrossRef]
  4. Sun-Waterhouse, D. The development of fruit-based functional foods targeting the health and wellness market: A review. Int. J. Food Sci. Technol. 2011, 46, 899–920. [Google Scholar] [CrossRef]
  5. Godoy, F.; Olivos-Hernandez, K.; Stange, C.; Handford, M. Abiotic Stress in Crop Species: Improving Tolerance by Applying Plant Metabolites. Plants 2021, 10, 186. [Google Scholar] [CrossRef]
  6. Ofori, K.F.; Antoniello, S.; English, M.M.; Aryee, A.N.A. Improving nutrition through biofortification-A systematic review. Front. Nutr. 2022, 9, 1043655. [Google Scholar] [CrossRef]
  7. Ristaino, J.B.; Anderson, P.K.; Bebber, D.P.; Brauman, K.A.; Cunniffe, N.J.; Fedoroff, N.V.; Finegold, C.; Garrett, K.A.; Gilligan, C.A.; Jones, C.M.; et al. The persistent threat of emerging plant disease pandemics to global food security. Proc. Natl. Acad. Sci. USA 2021, 118, e2022239118. [Google Scholar] [CrossRef] [PubMed]
  8. Yang, L.T.; Zhou, Y.F.; Wang, Y.Y.; Wu, Y.M.; Ye, X.; Guo, J.X.; Chen, L.S. Magnesium Deficiency Induced Global Transcriptome Change in Citrus sinensis Leaves Revealed by RNA-Seq. Int. J. Mol. Sci. 2019, 20, 3129. [Google Scholar] [CrossRef]
  9. Volz, R.K.; Oraguzie, N.C.; Whitworth, C.J.; How, N.; ChagnÃ, D.; Carlisle, C.M.; Gardiner, S.E.; Rikkerink, E.H.A.; Lawrence, T. Breeding for Red Flesh Colour in Apple: Progress and Challenges; International Society for Horticultural Science (ISHS): Leuven, Belgium, 2009; pp. 337–342. [Google Scholar]
  10. Espley, R.V.; Bovy, A.; Bava, C.; Jaeger, S.R.; Tomes, S.; Norling, C.; Crawford, J.; Rowan, D.; McGhie, T.K.; Brendolise, C.; et al. Analysis of genetically modified red-fleshed apples reveals effects on growth and consumer attributes. Plant Biotechnol. J. 2013, 11, 408–419. [Google Scholar] [CrossRef]
  11. Willits, M.G.; Kramer, C.M.; Prata, R.T.; De Luca, V.; Potter, B.G.; Steffens, J.C.; Graser, G. Utilization of the genetic resources of wild species to create a nontransgenic high flavonoid tomato. J. Agric. Food Chem. 2005, 53, 1231–1236. [Google Scholar] [CrossRef]
  12. Butelli, E.; Titta, L.; Giorgio, M.; Mock, H.P.; Matros, A.; Peterek, S.; Schijlen, E.G.; Hall, R.D.; Bovy, A.G.; Luo, J.; et al. Enrichment of tomato fruit with health-promoting anthocyanins by expression of select transcription factors. Nat. Biotechnol. 2008, 26, 1301–1308. [Google Scholar] [CrossRef]
  13. Firoozbady, E.; Young, T.R. Pineapple Plant Named Rosé; Del Monte Fresh Produce: Coral Gables, FL, USA, 2015. [Google Scholar]
  14. Li, Q.; Sapkota, M.; van der Knaap, E. Perspectives of CRISPR/Cas-mediated cis-engineering in horticulture: Unlocking the neglected potential for crop improvement. Hortic. Res. 2020, 7, 36. [Google Scholar] [CrossRef]
  15. Sharma, P.; Pandey, A.; Malviya, R.; Dey, S.; Karmakar, S.; Gayen, D. Genome editing for improving nutritional quality, post-harvest shelf life and stress tolerance of fruits, vegetables, and ornamentals. Front. Genome. Ed. 2023, 5, 1094965. [Google Scholar] [CrossRef]
  16. Charrier, A.; Vergne, E.; Dousset, N.; Richer, A.; Petiteau, A.; Chevreau, E. Efficient Targeted Mutagenesis in Apple and First Time Edition of Pear Using the CRISPR-Cas9 System. Front. Plant. Sci. 2019, 10, 40. [Google Scholar] [CrossRef] [PubMed]
  17. Jia, H.; Wang, N. Targeted genome editing of sweet orange using Cas9/sgRNA. PLoS ONE 2014, 9, e93806. [Google Scholar] [CrossRef] [PubMed]
  18. Kaur, N.; Alok, A.; Shivani; Kaur, N.; Pandey, P.; Awasthi, P.; Tiwari, S. CRISPR/Cas9-mediated efficient editing in phytoene desaturase (PDS) demonstrates precise manipulation in banana cv. Rasthali genome. Funct. Integr. Genom. 2018, 18, 89–99. [Google Scholar] [CrossRef]
  19. Malnoy, M.; Viola, R.; Jung, M.H.; Koo, O.J.; Kim, S.; Kim, J.S.; Velasco, R.; Nagamangala Kanchiswamy, C. DNA-Free Genetically Edited Grapevine and Apple Protoplast Using CRISPR/Cas9 Ribonucleoproteins. Front. Plant Sci. 2016, 7, 1904. [Google Scholar] [CrossRef]
  20. Naim, F.; Dugdale, B.; Kleidon, J.; Brinin, A.; Shand, K.; Waterhouse, P.; Dale, J. Gene editing the phytoene desaturase alleles of Cavendish banana using CRISPR/Cas9. Transgenic Res. 2018, 27, 451–460. [Google Scholar] [CrossRef]
  21. Nakajima, I.; Ban, Y.; Azuma, A.; Onoue, N.; Moriguchi, T.; Yamamoto, T.; Toki, S.; Endo, M. CRISPR/Cas9-mediated targeted mutagenesis in grape. PLoS ONE 2017, 12, e0177966. [Google Scholar] [CrossRef] [PubMed]
  22. Nishitani, C.; Hirai, N.; Komori, S.; Wada, M.; Okada, K.; Osakabe, K.; Yamamoto, T.; Osakabe, Y. Efficient Genome Editing in Apple Using a CRISPR/Cas9 system. Sci. Rep. 2016, 6, 31481. [Google Scholar] [CrossRef]
  23. Osakabe, Y.; Liang, Z.; Ren, C.; Nishitani, C.; Osakabe, K.; Wada, M.; Komori, S.; Malnoy, M.; Velasco, R.; Poli, M.; et al. CRISPR-Cas9-mediated genome editing in apple and grapevine. Nat. Protoc. 2018, 13, 2844–2863. [Google Scholar] [CrossRef]
  24. Pavese, V.; Moglia, A.; Corredoira, E.; Martinez, M.T.; Torello Marinoni, D.; Botta, R. First Report of CRISPR/Cas9 Gene Editing in Castanea sativa Mill. Front. Plant Sci. 2021, 12, 728516. [Google Scholar] [CrossRef]
  25. Ren, F.; Ren, C.; Zhang, Z.; Duan, W.; Lecourieux, D.; Li, S.; Liang, Z. Efficiency Optimization of CRISPR/Cas9-Mediated Targeted Mutagenesis in Grape. Front. Plant Sci. 2019, 10, 612. [Google Scholar] [CrossRef] [PubMed]
  26. Vaia, G.; Pavese, V.; Moglia, A.; Cristofori, V.; Silvestri, C. Knockout of phytoene desaturase gene using CRISPR/Cas9 in highbush blueberry. Front. Plant Sci. 2022, 13, 1074541. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, Z.; Wang, S.; Li, D.; Zhang, Q.; Li, L.; Zhong, C.; Liu, Y.; Huang, H. Optimized paired-sgRNA/Cas9 cloning and expression cassette triggers high-efficiency multiplex genome editing in kiwifruit. Plant Biotechnol. J. 2018, 16, 1424–1433. [Google Scholar] [CrossRef]
  28. Zhang, F.; LeBlanc, C.; Irish, V.F.; Jacob, Y. Rapid and efficient CRISPR/Cas9 gene editing in Citrus using the YAO promoter. Plant Cell Rep. 2017, 36, 1883–1887. [Google Scholar] [CrossRef] [PubMed]
  29. Zhu, C.; Zheng, X.; Huang, Y.; Ye, J.; Chen, P.; Zhang, C.; Zhao, F.; Xie, Z.; Zhang, S.; Wang, N.; et al. Genome sequencing and CRISPR/Cas9 gene editing of an early flowering Mini-Citrus (Fortunella hindsii). Plant Biotechnol. J. 2019, 17, 2199–2210. [Google Scholar] [CrossRef]
  30. Breitler, J.-C.; DeChamp, E.; Campa, C.; Rodriguez, L.A.Z.; Guyot, R.; Marraccini, P.; Etienne, H. CRISPR/Cas9-mediated efficient targeted mutagenesis has the potential to accelerate the domestication of Coffea canephora. Plant Cell Tissue Organ. Cult. 2018, 134, 383–394. [Google Scholar] [CrossRef]
  31. Chang, L.; Wu, S.; Tian, L. Effective genome editing and identification of a regiospecific gallic acid 4-O-glycosyltransferase in pomegranate (Punica granatum L.). Hortic. Res. 2019, 6, 123. [Google Scholar] [CrossRef]
  32. Bapat, V.A.; Jagtap, U.B.; Ghag, S.B.; Ganapathi, T.R. Molecular Approaches for the Improvement of Under-Researched Tropical Fruit Trees: Jackfruit, Guava, and Cusard Apple. Int. J. Fruit Sci. 2020, 20, 233–281. [Google Scholar] [CrossRef]
  33. Li, H.; Sun, H.; Dong, H.; Wang, S.; Fan, X.; Li, Y.; Cheng, L.; Zhang, Z.; Wang, Y.; Zhang, X.; et al. Genome editing of apple SQUAMOSA PROMOTER BINDNG PROTEIN-LIKE 6 enhances adventitious shoot regeneration. Plant Physiol. 2023, 191, 840–843. [Google Scholar] [CrossRef]
  34. Richard, B.; Qi, A.; Fitt, B.D. Control of crop diseases through Integrated Crop Management to deliver climate-smart farming systems for low- and high- input crop production. Plant Pathol. 2021, 71, 187–206. [Google Scholar] [CrossRef]
  35. Savary, S.; Willocquet, L.; Pethybridge, S.J.; Esker, P.; McRoberts, N.; Nelson, A. The global burden of pathogens and pests on major food crops. Nat. Ecol. Evol. 2019, 3, 430–439. [Google Scholar] [CrossRef] [PubMed]
  36. Youssef, K.; Ippolito, A.; Roberto, S.R. Editorial: Post-harvest Diseases of Fruit and Vegetable: Methods and Mechanisms of Action. Front. Microbiol. 2022, 13, 900060. [Google Scholar] [CrossRef]
  37. Fister, A.S.; Landherr, L.; Maximova, S.N.; Guiltinan, M.J. Transient Expression of CRISPR/Cas9 Machinery Targeting TcNPR3 Enhances Defense Response in Theobroma cacao. Front. Plant Sci. 2018, 9, 268. [Google Scholar] [CrossRef]
  38. Jia, H.; Zhang, Y.; Orbovic, V.; Xu, J.; White, F.F.; Jones, J.B.; Wang, N. Genome editing of the disease susceptibility gene CsLOB1 in citrus confers resistance to citrus canker. Plant Biotechnol. J. 2017, 15, 817–823. [Google Scholar] [CrossRef]
  39. Peng, A.; Chen, S.; Lei, T.; Xu, L.; He, Y.; Wu, L.; Yao, L.; Zou, X. Engineering canker-resistant plants through CRISPR/Cas9-targeted editing of the susceptibility gene CsLOB1 promoter in citrus. Plant Biotechnol. J. 2017, 15, 1509–1519. [Google Scholar] [CrossRef]
  40. Pompili, V.; Dalla Costa, L.; Piazza, S.; Pindo, M.; Malnoy, M. Reduced fire blight susceptibility in apple cultivars using a high-efficiency CRISPR/Cas9-FLP/FRT-based gene editing system. Plant Biotechnol. J. 2020, 18, 845–858. [Google Scholar] [CrossRef]
  41. Tripathi, J.N.; Ntui, V.O.; Ron, M.; Muiruri, S.K.; Britt, A.; Tripathi, L. CRISPR/Cas9 editing of endogenous banana streak virus in the B genome of Musa spp. overcomes a major challenge in banana breeding. Commun. Biol. 2019, 2, 46. [Google Scholar] [CrossRef]
  42. Wan, D.Y.; Guo, Y.; Cheng, Y.; Hu, Y.; Xiao, S.; Wang, Y.; Wen, Y.Q. CRISPR/Cas9-mediated mutagenesis of VvMLO3 results in enhanced resistance to powdery mildew in grapevine (Vitis vinifera). Hortic. Res. 2020, 7, 116. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, L.; Chen, S.; Peng, A.; Xie, Z.; He, Y.; Zou, Z. CRISPR/Cas9-mediated editing of CsWRKY22 reduces susceptibility to Xanthomonas citri subsp. citri in Wanjincheng orange (Citrus sinensis (L.) Osbeck). Plant Biotechnol. Rep. 2019, 13, 501–510. [Google Scholar] [CrossRef]
  44. Wang, X.; Tu, M.; Wang, D.; Liu, J.; Li, Y.; Li, Z.; Wang, Y.; Wang, X. CRISPR/Cas9-mediated efficient targeted mutagenesis in grape in the first generation. Plant Biotechnol. J. 2018, 16, 844–855. [Google Scholar] [CrossRef] [PubMed]
  45. Dash, P.K.; Rai, R. Translating the “Banana Genome” to Delineate Stress Resistance, Dwarfing, Parthenocarpy and Mechanisms of Fruit Ripening. Front. Plant Sci. 2016, 7, 1543. [Google Scholar] [CrossRef] [PubMed]
  46. Shao, X.; Wu, S.; Dou, T.; Zhu, H.; Hu, C.; Huo, H.; He, W.; Deng, G.; Sheng, O.; Bi, F.; et al. Using CRISPR/Cas9 genome editing system to create MaGA20ox2 gene-modified semi-dwarf banana. Plant Biotechnol. J. 2020, 18, 17–19. [Google Scholar] [CrossRef]
  47. Hu, C.; Sheng, O.; Deng, G.; He, W.; Dong, T.; Yang, Q.; Dou, T.; Li, C.; Gao, H.; Liu, S.; et al. CRISPR/Cas9-mediated genome editing of MaACO1 (aminocyclopropane-1-carboxylate oxidase 1) promotes the shelf life of banana fruit. Plant Biotechnol. J. 2021, 19, 654–656. [Google Scholar] [CrossRef]
  48. Ren, C.; Liu, X.; Zhang, Z.; Wang, Y.; Duan, W.; Li, S.; Liang, Z. CRISPR/Cas9-mediated efficient targeted mutagenesis in Chardonnay (Vitis vinifera L.). Sci. Rep. 2016, 6, 32289. [Google Scholar] [CrossRef] [PubMed]
  49. Miranda, S.; Piazza, S.; Nuzzo, F.; Li, M.; Lagreze, J.; Mithofer, A.; Cestaro, A.; Tarkowska, D.; Espley, R.; Dare, A.; et al. CRISPR/Cas9 genome-editing applied to MdPGT1 in apple results in reduced foliar phloridzin without impacting plant growth. Plant J. 2023, 113, 92–105. [Google Scholar] [CrossRef] [PubMed]
  50. Kumar, S.; Volz, R.; Alspach, P.; Bus, V. Development of a recurrent apple breeding programme in New Zealand: A synthesis of results, and a proposed revised breeding strategy. Euphytica 2010, 173, 207–222. [Google Scholar] [CrossRef]
  51. Laurens, F.; Durel, C.E.; Patocchia, A.; Peil, A.; Salvi, S.; Tartarini, S.; Velasco, R.; van de Weg, E. Review on apple genetics and breeding programmes and presentation of a new European initiative to increase fruit breeding efficiency. J. Fruit Sci. 2011, 27, 102–107. [Google Scholar]
  52. Cusin, R.; Revers, L.F.; Maraschin, F.S. New biotechnological tools to accelerate scab-resistance trait transfer to apple. Genet. Mol. Biol. 2017, 40 (Suppl. S1), 305–311. [Google Scholar] [CrossRef]
  53. Yamagishi, N.; Sasaki, S.; Yamagata, K.; Komori, S.; Nagase, M.; Wada, M.; Yamamoto, T.; Yoshikawa, N. Promotion of flowering and reduction of a generation time in apple seedlings by ectopical expression of the Arabidopsis thaliana FT gene using the Apple latent spherical virus vector. Plant Mol. Biol. 2011, 75, 193–204. [Google Scholar] [CrossRef]
  54. Flachowsky, H.; Peil, A.; Sopanen, T.; Elo, A.; Hanke, V. Overexpression of BpMADS4 from silver birch (Betula pendula Roth.) induces early-flowering in apple (Malus × domestica Borkh.). Plant Breed. 2007, 126, 137–145. [Google Scholar] [CrossRef]
  55. Trankner, C.; Lehmann, S.; Hoenicka, H.; Hanke, M.V.; Fladung, M.; Lenhardt, D.; Dunemann, F.; Gau, A.; Schlangen, K.; Malnoy, M.; et al. Over-expression of an FT-homologous gene of apple induces early flowering in annual and perennial plants. Planta 2010, 232, 1309–1324. [Google Scholar] [CrossRef]
  56. Flachowsky, H.; Szankowski, I.; Waidmann, S.; Peil, A.; Trankner, C.; Hanke, M.V. The MdTFL1 gene of apple (Malus × domestica Borkh.) reduces vegetative growth and generation time. Tree Physiol. 2012, 32, 1288–1301. [Google Scholar] [CrossRef]
  57. Flachowsky, H.; Le Roux, P.M.; Peil, A.; Patocchi, A.; Richter, K.; Hanke, M.V. Application of a high-speed breeding technology to apple (Malus × domestica) based on transgenic early flowering plants and marker-assisted selection. New Phytol. 2011, 192, 364–377. [Google Scholar] [CrossRef] [PubMed]
  58. Schlatholter, I.; Jansch, M.; Flachowsky, H.; Broggini, G.A.L.; Hanke, M.V.; Patocchi, A. Generation of advanced fire blight-resistant apple (Malus × domestica) selections of the fifth generation within 7 years of applying the early flowering approach. Planta 2018, 247, 1475–1488. [Google Scholar] [CrossRef]
  59. Smolka, A.; Li, X.Y.; Heikelt, C.; Welander, M.; Zhu, L.H. Effects of transgenic rootstocks on growth and development of non-transgenic scion cultivars in apple. Transgenic Res. 2010, 19, 933–948. [Google Scholar] [CrossRef] [PubMed]
  60. Flachowsky, H.; Szankowski, I.; Fischer, T.C.; Richter, K.; Peil, A.; Hofer, M.; Dorschel, C.; Schmoock, S.; Gau, A.E.; Halbwirth, H.; et al. Transgenic apple plants overexpressing the Lc gene of maize show an altered growth habit and increased resistance to apple scab and fire blight. Planta 2010, 231, 623–635. [Google Scholar] [CrossRef] [PubMed]
  61. APHIS. Questions and Answer: Arctic Apple Deregulation; APHIS: Washington, DC, USA, 2015.
  62. Velasco, R.; Zharkikh, A.; Affourtit, J.; Dhingra, A.; Cestaro, A.; Kalyanaraman, A.; Fontana, P.; Bhatnagar, S.K.; Troggio, M.; Pruss, D.; et al. The genome of the domesticated apple (Malus × domestica Borkh.). Nat. Genet. 2010, 42, 833–839. [Google Scholar] [CrossRef]
  63. Klocko, A.L.; Borejsza-Wysocka, E.; Brunner, A.M.; Shevchenko, O.; Aldwinckle, H.S.; Strauss, S.H. Transgenic Suppression of AGAMOUS Genes in Apple Reduces Fertility and Increases Floral Attractiveness. PLoS ONE 2016, 11, e0159421. [Google Scholar] [CrossRef]
  64. Culley, T.M.; Hardiman, N.A. The Beginning of a New Invasive Plant: A History of the Ornamental Callery Pear in the United States. BioScience 2007, 57, 956–964. [Google Scholar] [CrossRef]
  65. Gasteiger, E.; Gattiker, A.; Hoogland, C.; Ivanyi, I.; Appel, R.D.; Bairoch, A. ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003, 31, 3784–3788. [Google Scholar] [CrossRef] [PubMed]
  66. Ireland, H.S.; Tomes, S.; Hallet, I.C.; Karunairetnam, S.; David, K.M.; Yao, J.L.; Schaeffer, R.J. Coreless apples generated by the suppression of carpel genes and hormone-induced fruit set. Fruit Res. 2021, 1, 2. [Google Scholar] [CrossRef]
  67. Ran, F.A.; Hsu, P.D.; Wright, J.; Agarwala, V.; Scott, D.A.; Zhang, F. Genome engineering using the CRISPR-Cas9 system. Nat. Protoc. 2013, 8, 2281–2308. [Google Scholar] [CrossRef] [PubMed]
  68. Elorriaga, E.; Klocko, A.L.; Ma, C.; Strauss, S.H. Variation in Mutation Spectra among CRISPR/Cas9 Mutagenized Poplars. Front. Plant Sci. 2018, 9, 594. [Google Scholar] [CrossRef]
  69. Malabarba, J.; Chevreau, E.; Dousset, N.; Veillet, F.; Moizan, J.; Vergne, E. New Strategies to Overcome Present CRISPR/Cas9 Limitations in Apple and Pear: Efficient Dechimerization and Base Editing. Int. J. Mol. Sci. 2020, 22, 319. [Google Scholar] [CrossRef]
  70. Norris, S.R.; Barrette, T.R.; DellaPenna, D. Genetic dissection of carotenoid synthesis in arabidopsis defines plastoquinone as an essential component of phytoene desaturation. Plant Cell 1995, 7, 2139–2149. [Google Scholar] [CrossRef]
  71. Chen, Y.H.; Sharma, S.; Bewg, W.P.; Xue, L.J.; Gizelbach, C.R.; Tsai, C.J. Multiplex Editing of the Nucleoredoxin1 Tandem Array in Poplar: From Small Indels to Translocations and Complex Inversions. CRISPR J. 2023, 6, 339–349. [Google Scholar] [CrossRef]
  72. Elorriaga, E.; Klocko, A.L.; Ma, C.; du Plessis, M.; An, X.; Myburg, A.A.; Strauss, S.H. Genetic containment in vegetatively propagated forest trees: CRISPR disruption of LEAFY function in Eucalyptus gives sterile indeterminate inflorescences and normal juvenile development. Plant Biotechnol. J. 2021, 19, 1743–1755. [Google Scholar] [CrossRef]
  73. USDA National Institute of Food and Agriculture. Apple Rootstock Info: M.26. Available online: https://apples.extension.org/apple-rootstock-info-m-26/ (accessed on 22 May 2023).
  74. Ten Hove, H.W. Apple Tree—Royal Gala Variety; Stark Brothers Nurseries & Orchards Company: Louisiana, MO, USA, 1977. [Google Scholar]
  75. Borejsza-Wysocka, E.; Norelli, J.L.; Ko, K.; Aldwinckle, H.S. Transformation of authentic M26 for enhanced resistance to fire blight. Acta Hortic. 1999, 489, 259–266. [Google Scholar] [CrossRef]
  76. Ko, K.; Reynoird, J.-P.; Aldwinckle, H.S.; Brown, S.K. T4 lysozyme and attacin genes enhance resistance of transgenic ‘Galaxy’ apple against Erwinia amylorvora. J. Am. Soc. Hortic. Sci. 2002, 127, 515–519. [Google Scholar] [CrossRef]
  77. Norelli, J.L.; Aldwinckle, H.S.; Beer, S.V. Leaf wounding increases efficiency of Agrobacterium-mediated transformation of apple. Phytopath 1988, 78, 1292–1297. [Google Scholar] [CrossRef]
  78. Bolar, J.P.; Hanke, V.; Norelli, J.L.; Aldwinckle, H.S. An efficient method for rooting and acclimation of micropropagated apple cultivars. HortSci 1998, 33, 1251–1252. [Google Scholar] [CrossRef]
  79. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef] [PubMed]
  80. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, R.; McWilliam, H.; Remmert, M.; Soding, J.; et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol. 2011, 7, 539. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Location of sgRNAs in MADS15 and MADS221 and construct diagrams. (A) Schematic of targeting locations in MADS15 and MADS221. Exons are shown with boxes, intron with lines, and scissors represent general locations of sgRNAs. The MADS15 and MADS221 genes are predicted to have the same overall gene structure. As targeting constructs had two sgRNAs, two scissors are shown to indicate the approximate cut sites in exon 1. The total length of MADS15 and MADS221 is ~3102 base pairs. Gene structure is from the alignment of coding and genomic sequences obtained from Phytozome. (B) Construct diagrams. We created four double sgRNA constructs for genome editing (the general structure of which is shown here) and a Cas9-only control construct. LB, T-DNA left border; nptII, neomycin phosphotransferase II gene for kanamycin selection in plants; U6, U6-26 gene promoter from Arabidopsis thaliana; sgRNA, locations of sgRNAs; 2 × 35S, double 35S gene promoter from Cauliflower mosaic virus; Cas9 human, codon-optimized Cas9 from Streptococcus pyogenes; tnos, nopaline synthetase terminator from Agrobacterium tumefaciens; RB, T-DNA right border; specR, spectinomycin resistance gene for selection in bacteria. Promoter directionality is shown with arrows.
Figure 1. Location of sgRNAs in MADS15 and MADS221 and construct diagrams. (A) Schematic of targeting locations in MADS15 and MADS221. Exons are shown with boxes, intron with lines, and scissors represent general locations of sgRNAs. The MADS15 and MADS221 genes are predicted to have the same overall gene structure. As targeting constructs had two sgRNAs, two scissors are shown to indicate the approximate cut sites in exon 1. The total length of MADS15 and MADS221 is ~3102 base pairs. Gene structure is from the alignment of coding and genomic sequences obtained from Phytozome. (B) Construct diagrams. We created four double sgRNA constructs for genome editing (the general structure of which is shown here) and a Cas9-only control construct. LB, T-DNA left border; nptII, neomycin phosphotransferase II gene for kanamycin selection in plants; U6, U6-26 gene promoter from Arabidopsis thaliana; sgRNA, locations of sgRNAs; 2 × 35S, double 35S gene promoter from Cauliflower mosaic virus; Cas9 human, codon-optimized Cas9 from Streptococcus pyogenes; tnos, nopaline synthetase terminator from Agrobacterium tumefaciens; RB, T-DNA right border; specR, spectinomycin resistance gene for selection in bacteria. Promoter directionality is shown with arrows.
Plants 12 03693 g001
Figure 2. Images from transformation and regeneration. (A) Callus forming on regeneration plates after co-cultivation; (B) shoots forming after 4 weeks in the light on regeneration medium; (C) larger shoots were chosen and transferred to shoot proliferation medium; (D) regenerants in shoot proliferation medium; leaves from this stage were used for transformation event confirmation via PCR.
Figure 2. Images from transformation and regeneration. (A) Callus forming on regeneration plates after co-cultivation; (B) shoots forming after 4 weeks in the light on regeneration medium; (C) larger shoots were chosen and transferred to shoot proliferation medium; (D) regenerants in shoot proliferation medium; leaves from this stage were used for transformation event confirmation via PCR.
Plants 12 03693 g002
Figure 3. Examples of edits in MADS15. (A) Examples of small edits in MADS15 from different independent events at the universal sgRNA3 site shown in purple. (B) Examples of large deletion events obtained with sgRNA1 shown in green and those obtained with sgRNA2 shown in yellow. Numbers indicate the amount of inserted (+) or deleted (−) bases.
Figure 3. Examples of edits in MADS15. (A) Examples of small edits in MADS15 from different independent events at the universal sgRNA3 site shown in purple. (B) Examples of large deletion events obtained with sgRNA1 shown in green and those obtained with sgRNA2 shown in yellow. Numbers indicate the amount of inserted (+) or deleted (−) bases.
Plants 12 03693 g003
Figure 4. Examples of edits in MADS221. (A) Examples of small edits in MADS221 from different independent events, with the sgRNA1 site shown in green and the sgRNA3 site shown in purple. (B) Examples of large deletion in events obtained at the sgRNA2 site shown in yellow and at the sgRNA3 site shown in purple. Numbers indicate the amount of inserted (+) or deleted (−) bases.
Figure 4. Examples of edits in MADS221. (A) Examples of small edits in MADS221 from different independent events, with the sgRNA1 site shown in green and the sgRNA3 site shown in purple. (B) Examples of large deletion in events obtained at the sgRNA2 site shown in yellow and at the sgRNA3 site shown in purple. Numbers indicate the amount of inserted (+) or deleted (−) bases.
Plants 12 03693 g004
Figure 5. Predicted peptides of MADS15 and MADS221 from alleles with small indels. (A) Alignment of predicted peptides for the beginning of MADS15 for WT, Cas9, and four different events with small indels. (B) Alignment of predicted peptides for the beginning of MADS221. Dashes indicate absent amino acids, asterisks in alignments indicate stop codons, asterisks under alignments indicate exact sequence matches, colons and periods under alignment indicate partial matches, and numbers on the far right indicate the length of the predicted fragment.
Figure 5. Predicted peptides of MADS15 and MADS221 from alleles with small indels. (A) Alignment of predicted peptides for the beginning of MADS15 for WT, Cas9, and four different events with small indels. (B) Alignment of predicted peptides for the beginning of MADS221. Dashes indicate absent amino acids, asterisks in alignments indicate stop codons, asterisks under alignments indicate exact sequence matches, colons and periods under alignment indicate partial matches, and numbers on the far right indicate the length of the predicted fragment.
Plants 12 03693 g005
Table 2. Transformation events obtained per construct and cultivar. We obtained a total of 44 transgenic events, 6 Cas9 control events and 38 CRISPR events, from our four different CRISPR constructs. Event refers to an independent plant transformation occurrence.
Table 2. Transformation events obtained per construct and cultivar. We obtained a total of 44 transgenic events, 6 Cas9 control events and 38 CRISPR events, from our four different CRISPR constructs. Event refers to an independent plant transformation occurrence.
ConstructCultivarTransformed Events
Cas9‘M.26’5
Cas9‘Royal Gala’1
Total Cas9 eventsBoth6
sgRNA1 sgRNA2‘M.26’8
sgRNA 1 sgRNA2‘Royal Gala’1
sgRNA15 sgRNA3‘M.26’6
sgRNA15 sgRNA3‘Royal Gala’1
sgRNA3 sgRNA2‘M.26’8
sgRNA3 sgRNA2‘Royal Gala’0
sgRNA22 sgRNA3‘M.26’13
sgRNA22 sgRNA3‘Royal Gala’1
Total CRISPR events ‘M.26’‘M.26’35
Total CRISPR events ‘Royal Gala’‘Royal Gala’3
Total CRISPR eventsBoth38
Total transgenic events (Cas9 and CRISPR)Both44
Table 3. Summary of whether editing occurred in each target gene allele and frequencies. We analyzed all 38 CRISPR-Cas9 events for changes in alleles of each target gene. Events with no changes in a given target gene and alleles were designated WT for that gene, events with one WT and one altered allele were designated heterozygous (het) for that gene, and events with alterations to both alleles of a gene were designated dual-edited for that gene. Dual-edited indicates either two different edited alleles, or identical edited alleles.
Table 3. Summary of whether editing occurred in each target gene allele and frequencies. We analyzed all 38 CRISPR-Cas9 events for changes in alleles of each target gene. Events with no changes in a given target gene and alleles were designated WT for that gene, events with one WT and one altered allele were designated heterozygous (het) for that gene, and events with alterations to both alleles of a gene were designated dual-edited for that gene. Dual-edited indicates either two different edited alleles, or identical edited alleles.
Outcome Per GeneEvents
(Number and %)
MADS15MADS221
WTWT4 (11%)
WThet2 (5%)
hetWT3 (8%)
WTDual-edited3 (8%)
Dual-editedWT0 (0%)
hethet4 (11%)
HetDual-edited5 (13%)
Dual-editedhet3 (8%)
Dual-editedDual-edited14 (37%)
Fully WT4 (11%)
Mix of WT and edited20 (53%)
Fully edited14 (37%)
Table 4. Mutation frequency by sgRNA and MADS-like genes. We used three sgRNAs that could target both genes in each cultivar. These dual-gene targeting sgRNAs (sgRNA1, sgRNA2, sgRNA3) could target both alleles of MADS15 and MADS221. The MADS15-specific (sgRNA15) and MADS221-specific (sgRNA22) sgRNAs could each target either MADS15 or MADS221 in both cultivars. As guide RNAs were used in pairs to create each construct and variable events were obtained for each construct, the number of potential targeting sites per guide RNA varies. In total, our sgRNAs had 262 possible sites. Targeting was defined as at least one change from the cultivar-specific WT sequence for a given allele.
Table 4. Mutation frequency by sgRNA and MADS-like genes. We used three sgRNAs that could target both genes in each cultivar. These dual-gene targeting sgRNAs (sgRNA1, sgRNA2, sgRNA3) could target both alleles of MADS15 and MADS221. The MADS15-specific (sgRNA15) and MADS221-specific (sgRNA22) sgRNAs could each target either MADS15 or MADS221 in both cultivars. As guide RNAs were used in pairs to create each construct and variable events were obtained for each construct, the number of potential targeting sites per guide RNA varies. In total, our sgRNAs had 262 possible sites. Targeting was defined as at least one change from the cultivar-specific WT sequence for a given allele.
sgRNATotal SitesTargeting (Number and %)
Dual-gene targeting
sgRNA13617 (47%)
sgRNA26826 (38%)
sgRNA311679 (68%)
MADS15 specific
sgRNA151410 (71%)
MADS221 specific
sgRNA222817 (61%)
All guide RNAs
sgRNA1, sgRNA2, sgRNA3, sgRNA15, sgRNA22262149 (57%)
Table 5. Editing by CRISPR construct. We created one Cas9-only control and four different gene-editing CRISPR-Cas9 constructs, with 7 to 14 events per construct, combined between cultivars. Editing refers to any change from the WT sequence.
Table 5. Editing by CRISPR construct. We created one Cas9-only control and four different gene-editing CRISPR-Cas9 constructs, with 7 to 14 events per construct, combined between cultivars. Editing refers to any change from the WT sequence.
ConstructNumber of Transformation EventsNumber of Edited Events% Editing
Cas9 control600
sgRNA1 sgRNA29889
sgRNA2 sgRNA38675
sgRNA15 sgRNA377100
sgRNA22 sgRNA3141393
Table 6. Frequencies of deletions, insertions, or small indels. Insertion and deletion mutations were classified as small if they involved 15 or fewer bp, and large if they were 16 or more bp. The percentages are calculated as the frequency of that mutation type for all changed alleles for a construct. The most common mutation type per construct is bolded.
Table 6. Frequencies of deletions, insertions, or small indels. Insertion and deletion mutations were classified as small if they involved 15 or fewer bp, and large if they were 16 or more bp. The percentages are calculated as the frequency of that mutation type for all changed alleles for a construct. The most common mutation type per construct is bolded.
ConstructTotal Events (N)Total Alleles (N)Total Edited Alleles (N) %Large
Deletion (N) %
Small
Deletion (N) %
Small
Insertion (N) %
Mix of Small Indels
(N) %
sgRNA1 sgRNA293617 (47%)14 (82%)3 (18%)0 (0%)0 (0%)
sgRNA2 sgRNA383216 (50%)10 (63%)3 (19%)3 (19%)0 (0%)
sgRNA15 sgRNA372822 (79%)14 (50%)3 (11%)5 (18%)0 (0%)
sgRNA22 sgRNA3145643 (77%)20 (36%)12 (21%)4 (7%)7 (13%)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jacobson, S.; Bondarchuk, N.; Nguyen, T.A.; Canada, A.; McCord, L.; Artlip, T.S.; Welser, P.; Klocko, A.L. Apple CRISPR-Cas9—A Recipe for Successful Targeting of AGAMOUS-like Genes in Domestic Apple. Plants 2023, 12, 3693. https://doi.org/10.3390/plants12213693

AMA Style

Jacobson S, Bondarchuk N, Nguyen TA, Canada A, McCord L, Artlip TS, Welser P, Klocko AL. Apple CRISPR-Cas9—A Recipe for Successful Targeting of AGAMOUS-like Genes in Domestic Apple. Plants. 2023; 12(21):3693. https://doi.org/10.3390/plants12213693

Chicago/Turabian Style

Jacobson, Seth, Natalie Bondarchuk, Thy Anh Nguyen, Allison Canada, Logan McCord, Timothy S. Artlip, Philipp Welser, and Amy L. Klocko. 2023. "Apple CRISPR-Cas9—A Recipe for Successful Targeting of AGAMOUS-like Genes in Domestic Apple" Plants 12, no. 21: 3693. https://doi.org/10.3390/plants12213693

APA Style

Jacobson, S., Bondarchuk, N., Nguyen, T. A., Canada, A., McCord, L., Artlip, T. S., Welser, P., & Klocko, A. L. (2023). Apple CRISPR-Cas9—A Recipe for Successful Targeting of AGAMOUS-like Genes in Domestic Apple. Plants, 12(21), 3693. https://doi.org/10.3390/plants12213693

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

Article Metrics

Back to TopTop