Introduction {#Sec1}
============

As an important branch of agriculture, horticulture originated thousands of years ago and has developed greatly during the course of human history. Horticultural crops are generally considered to include vegetable and fruit crops as well as floricultural and ornamental plants, which are cultivated for food, for nutritional and medical use, and for esthetic enjoyment^[@CR1]^. Vegetable and fruit crops are low in calories but contain high levels of vitamins and minerals^[@CR2]^, making them indispensable for balancing our daily diet. Although the supply of horticultural products is increasing, the diversity and nutritional value of the products are decreasing^[@CR3]^. These decreases can be partially attributed to the narrow genetic diversity of horticultural crops resulting from domestication and breeding as well as reproductive barriers that inhibit genetic introgression from wild relatives. Therefore, the generation of genetic resources with diverse and desirable characteristics will be of great value for improving horticultural products.

Thousands of years ago, humans began to improve crops by introducing new traits from crossable relatives. The essential goal of this process was the transfer of desirable genetic variations. As late as 1930s, the available variations were generated solely through natural or spontaneous processes. Breeders subsequently learned to produce mutants by using chemical mutagens or radiation^[@CR4]^. Both spontaneous and induced mutations have significantly increased crop yield and quality^[@CR5]^. Given the rareness and randomness of these mutations, however, obtaining suitable materials for crop improvement has proven to be laborious and time consuming^[@CR4]^.

With the rapid progress in molecular biology, DNA sequence-specific manipulation has become a powerful tool. In 1987, several animal scientists invented gene-targeting technology that relies on homologous recombination (HR). This innovative technology enabled researchers to precisely edit (though with a low frequency) an endogenous gene after introducing a donor template into mouse embryonic stem cells^[@CR6],[@CR7]^. Similar progress was subsequently reported by plant researchers, but with an extremely low editing frequency of 0.5--7.2 × 10^−4 [@CR8],[@CR9]^. DNA double-stranded breaks (DSBs), which commonly result in HR in meiotic chromosomes^[@CR10]^, were later used to increase the HR frequency in gene targeting^[@CR11]^. In addition to HR, DSBs can be repaired through the error-prone nonhomologous end-joining (NHEJ) pathway in somatic cells, which can generate mutations via the small deletions or insertions that occur at a break site^[@CR12]^. Scientists have used the following kinds of engineered endonucleases to introduce site-specific DSBs: meganucleases (MNs), zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated 9 (Cas9), and CRISPR from *Prevotella* and *Francisella* 1 (CRISPR/Cpf1). These engineered endonucleases have enabled genome editing in various biological systems^[@CR13]--[@CR16]^.

With the advent of CRISPR/Cas9, the application of genome editing to horticultural crops has greatly advanced. In this review, we first introduce and compare the engineered nucleases that are used for genome editing. We then consider their current applications in horticulture. Finally, we discuss the implications and challenges of genome editing for the improvement of horticultural crops.

Genome-editing systems {#Sec2}
======================

Sequence-specific DNA binding, such as the interaction between a transcription factor and a promoter, is a common phenomenon. For genome editing, the previously mentioned nucleases can target specific sequences to generate DSBs under the guidance of protein--DNA interaction (for MNs, ZFNs, and TALENs) or RNA--DNA base-pairing (for CRISPR/Cas9 and CRISPR/Cpf1)^[@CR16],[@CR17]^.

Meganucleases or homing nucleases {#Sec3}
---------------------------------

The first class of nucleases for genome editing, MNs or homing endonucleases, was discovered in the genomes of microorganisms or organelles. By recognizing DNA sequence elements ranging from 12 to 40 bp, these nucleases cut both strands of DNA in a site-specific manner (Fig. [1a](#Fig1){ref-type="fig"})^[@CR18]^. Among MNs, the I-*Cre*I protein has received the most research attention and has been reported to be effective in maize^[@CR19]^, but the rare occurrence of recognizable sites limits the ability of I-*Cre*I and other MNs to edit desired target sites^[@CR17]^. To broaden the application of MNs, researchers have used mutagenesis or combinatorial assembly to produce MN variants that target the desired DNA sequence^[@CR20],[@CR21]^. Nevertheless, the overlapping recognition and catalytic domains of modified MNs cause difficulties and often compromise their catalytic activity^[@CR15]^. For these reasons, MNs have not been widely used by plant scientists.Fig. 1Schematic models of genome-editing systems.**a** A meganuclease can recognize a DNA sequence element of 12--40 bp and cut both strands at specific sites, forming sticky double-stranded breaks (DSBs). **b** In ZFNs, each zinc finger recognizes a 3-bp DNA sequence. Target specificity is achieved by arrays of several zinc fingers. Each DNA strand is bound by one zinc finger array linked with FokI, which in dimer form cuts DNA strands. **c** In TALENs, the central binding domain of each TALE consists of 13--28 repeats. Each repeat (a highly conserved sequence of 34 amino acids) can recognize and bind one nucleotide through the variable di-residues at the 12th and 13th positions. Paired TALENs lead to the dimerization of FokI, and the dimers cut the DNA stands, forming sticky DSBs at the target site. **d** In the CRISPR/Cas9 system, a single guide RNA (sgRNA) pairs with the target sequence upstream of a 5′-NGG-3′ PAM motif (N=A, T, C or G). The Cas9 endonuclease cuts the noncomplementary and complementary DNA strands at a location 3 nucleotides upstream of the PAM motif with RuvC and HNH domains, respectively. The cutting forms a blunt end DSB. **e** In the CRISPR/Cpf1 system, target specificity is achieved by the pairing of crRNA with the DNA strand downstream of a 5′-TTN-3′ PAM motif. The Cpf1 endonuclease uses the RuvC and Nuc domains to cut noncomplementary and complementary DNA strands at different positions, producing DSBs with sticky ends

ZFNs and TALENs {#Sec4}
---------------

As suggested by their names, ZFNs or TALENs are generated by fusing the DNA cleavage domain of the endonuclease FokI with zinc fingers (ZFs) or with transcriptional activator-like effectors (TALEs). The FokI endonuclease domain mediates independent and nonspecific DNA cleavage upon dimerization and is not involved in any sequence recognition^[@CR22]^. Therefore, a pair of ZFs or TALEs, each fused with a FokI endonuclease domain, is designed to achieve site-specific cleavage^[@CR23]--[@CR25]^. ZFs are found in transcription factors, with each finger domain recognizing three specific nucleotides. ZFNs typically exhibit an array of 3 or 4 finger domains, which can recognize 18--24 bp sequences when a ZFN occurs as a dimer^[@CR23],[@CR25]^. Many studies have been conducted to improve ZFN applicability, efficiency, and precision^[@CR26],[@CR27]^, but there are still concerns about interference from neighboring finger domains and the limited number of recognition sites (Fig. [1b](#Fig1){ref-type="fig"})^[@CR15]^.

In contrast to ZFNs, TALENs achieve sequence specificity via the customizable DNA-binding domains of TALEs, which are proteins excreted by the common bacterial plant pathogen *Xanthomonas*^[@CR28]^. During pathogenesis, TALEs bind to a specific sequence of plant promoters to activate gene expression to facilitate infection^[@CR28]^. The central binding domain of TALEs consists of 13--28 repeat sequences. Each repeat, which encodes a highly conserved sequence of 34 amino acids, can recognize and bind to one nucleotide through the variable di-residues at the 12th and 13th positions^[@CR29]--[@CR31]^. Such one-to-one pairing, together with the negligible context dependency on neighboring repeats, enables TALENs to target desired sequences (Fig. [1c](#Fig1){ref-type="fig"})^[@CR32],[@CR33]^. In general, TALENs outperform ZFNs in terms of precision and accessibility.

CRISPR/Cas9 and CRISPR/Cpf1 {#Sec5}
---------------------------

Unlike ZFN and TALEN systems, which depend on protein--DNA binding specificity, the CRISPR system relies on RNA--DNA binding to achieve sequence specificity. During the functional elucidation of the CRISPR/Cas system, its involvement in bacterial resistance to viruses was experimentally demonstrated^[@CR34]^, and several components, including crRNA, PAM motif, and tracrRNA, were discovered to be necessary for this system^[@CR35]--[@CR37]^. More interestingly, reconstructed key components of the CRISPR/Cas9 system can introduce DSBs in a site-specific way, suggesting the potential use of this programmable RNA-guided CRISPR/Cas9 system for genome editing in organisms other than bacteria^[@CR38],[@CR39]^. This possibility was soon demonstrated in human and mouse cells^[@CR40]--[@CR42]^, zebrafish^[@CR43]^, and plants^[@CR44]--[@CR48]^. In the system, site-specific binding to the target is achieved via RNA-DNA pairing of a 20-nt sequence in the chimeric single-guide RNA (sgRNA) with the target. The other crRNA- and tracrRNA-derived sequences also interact with the target to form an RNA:DNA heteroduplex that is recognized by the collective interactions of several Cas9 domains: PI, REC1, RuvC, and NUC. Thereafter, the RuvC and HNH domains cut the noncomplementary and complementary DNA strands at a location 3 nucleotides upstream of the PAM motif, respectively (Fig. [1d](#Fig1){ref-type="fig"}). The recognizable PAM motif of Cas9 is 5′-NGG-3′ (N=A, T, C, or G), and this G-rich feature prevents the design of sgRNAs in T-rich regions^[@CR49]^.

Cpf1, another endonuclease in the class 2 Type V CRISPR system, has also been found to be efficient in plant genome editing^[@CR50]^ and to present unique features^[@CR51]^. First, Cpf1 does not require an additional tracrRNA to form a mature crRNA. Second, unlike Cas9, which recognizes G-rich PAM sequences, Cpf1 recognizes T-rich PAM sequences. Finally, whereas cutting by the Cas9 endonuclease produces blunt ends, cutting by the Cpf1 endonuclease produces cohesive ends (Fig. [1e](#Fig1){ref-type="fig"}). In addition to causing site-specific mutations, CRISPR genome-editing systems can be used to achieve gene regulation^[@CR52],[@CR53]^ through the manipulation of the nuclease-inactivated Cas9 (dCas9).

Each of the endonucleases used for genome editing has unique properties because of differences in their underlying mechanisms (Fig. [1](#Fig1){ref-type="fig"} and Table [1](#Tab1){ref-type="table"}, Zhang et al.^[@CR16],[@CR54]^; Knott and Doudna^[@CR55]^). In addition to generating indel mutations at target sequences, CRISPR/Cas systems have been adapted for precise base editing^[@CR56]--[@CR59]^. Base editors usually consist of an sgRNA-guided Cas9 nickase (nCas9) fused with a deaminase that causes C to T or A to G base conversions. These resources greatly increase the versatility of the tools that can be used for precise manipulation of horticultural crops.Table 1Comparison of genome-editing systems\*PropertyMNsZFNsTALENsCRISPR/Cas9 or CRISPR/Cpf1Site-recognition domainMN binding domainZinc fingersTranscription activator-like effectorssgRNA or crRNAInteraction patternProtein--DNAProtein--DNAProtein--DNARNA--DNA pairingDNA cleavageMNsFokIFokICas9 or Cpf1Available sites\*\*1/1000 bp1/140 bpAny site (in principle)1/13 bpPrecision++++++++++ or ++++Efficiency++++++ or ++Ease of design++++++++++ or +++++Specificity++++++++ or +++Multiplex editing+++++++ or ++++\*This table is based on Boglioli and Richard^[@CR60]^, Rocha-Martins et al.^[@CR17]^, and Zhang et al.^[@CR16]^. "+" indicates the level\*\*This information is based on human genome data

Current status of genome editing in horticultural crops {#Sec6}
=======================================================

To obtain genetic resources with diverse characteristics for breeding, both spontaneous and induced mutations have been commonly used^[@CR60]^. The rareness and uncertainty of these mutations have motivated scientists to find ways to introduce precise mutations at target sites^[@CR15],[@CR17]^. Recently, most genome-editing studies on plants have been carried out in model systems and staple crops^[@CR44]--[@CR46]^, but the application of genome editing to horticultural crops is rapidly increasing^[@CR61]^. In 2013, the first example of genome editing in a horticultural crop was achieved via a TALEN in *Brassica* *oleracea*^[@CR62]^. In the following years, the number of studies involving genome editing in horticulture has exponentially increased (Fig. [2a](#Fig2){ref-type="fig"}, Table [2](#Tab2){ref-type="table"}), and CRISPR-based systems now dominate. The functions of genes targeted by genome editing are very diverse, but researchers have focused most on targets affecting development, followed by targets affecting metabolism and stress responses. In addition, studies that focus on the improvement of the CRISPR/Cas9 system in horticultural crops frequently use marker/reporter genes as targets such as *phytoene desaturase* (*PDS*), whose mutation results in an albino phenotype (Fig. [2b](#Fig2){ref-type="fig"}). Among horticultural crops, tomato has received much more attention regarding genome editing than other crops: \~42% of genome-editing studies have involved tomato, whereas \~13% have involved potato. Although most (72%) genome editing with horticultural crops is performed in vegetables (Fig. [2c](#Fig2){ref-type="fig"}), some floral and medicinal plants have also been successfully manipulated by genome editing (Fig. [2c](#Fig2){ref-type="fig"}).Fig. 2Number of research articles involving gene editing.The information used in this figure was retrieved through May 31 of 2019. According to the information from https://aps.dac.gov.in/Public/Crops.pdf, horticultural crops include vegetables, fruits, florals, and medicinal plants. **a** The number of research articles involving the editing of horticultural crops with ZFNs, TALENs, and CRISPR/Cas9 from 2013 to 2019 (only the first 5 months). **b** The number of research articles in which the edited genes were mainly associated with development, metabolism, stress tolerance and other functions. **c** The number of research articles involving gene editing of different kinds of horticultural cropsTable 2A list of publications on genome editing in horticultural cropsSpeciesCrop typeGenome editing toolTargeted geneGene function or phenotypeClassification of targeted geneReference*Solanum lycopersicum*VegetableCRISPR*SlALS1*Enhanced herbicide resistanceStress response^[@CR103]^*Solanum lycopersicum*VegetableCRISPR*SlJAZ2*Resistance to bacterial speckStress response^[@CR104]^*Solanum lycopersicum*VegetableCRISPR*APETALA2a (AP2a), NON-RIPENING (NOR) and FRUITFULL (FUL1/TDR4 and FUL2/MBP7)*Fruit development and ripeningDevelopment^[@CR105]^*Solanum lycopersicum*VegetableCRISPR*Pectate lyase (PL), polygalacturonase 2a (PG2a), and beta-galactanase (TBG4)*Cell wall gene, altered fruit color and firmnessDevelopment^[@CR106]^Solanum lycopersicumVegetableCRISPR*SlNPR1*Reduced drought toleranceStress response^[@CR107]^*Solanum lycopersicum*VegetableCRISPR*SlALS1, SlALS2*Enhanced herbicide resistanceStress response^[@CR101]^*Solanum lycopersicum*VegetableCRISPR*SlGAI*Gibberellin response and dwarfismDevelopment^[@CR108]^*Solanum lycopersicum*VegetableCRISPR*SlEIN2, SlERFE1, SlARF2B, SlGRAS8, SlACS2, SlACS4*Ethylene response and fruit developmentDevelopment^[@CR97]^*Solanum lycopersicum*VegetableCRISPR*SBPase*Leaf senescence (SBPase in primary metabolism)Metabolism^[@CR109]^*Solanum lycopersicum*VegetableCRISPR*CBF1*Chilling toleranceStress response^[@CR110]^*Solanum lycopersicum*VegetableCRISPR*POLYGALACTURONASE (PG) and PECTATE LYASE (PL)*Cell wall geneDevelopment^[@CR111]^*Solanum lycopersicum*VegetableCRISPR*NPTII*N.A.Others^[@CR112]^*Solanum lycopersicum*VegetableCRISPR*Psy1 and CrtR-b2*Carotenoid metabolismMetabolism^[@CR113]^*Solanum lycopersicum*VegetableCRISPR*NADK2A, IAA9*NAD Kinase 2A; IAA9Development^[@CR114]^*Solanum lycopersicum*VegetableCRISPR*DDM1a, b*Decrease in DNA methylationDevelopment^[@CR115]^*Solanum lycopersicum*VegetableCRISPR*SlMAPK20*Aborted pollen developmentDevelopment^[@CR116]^*Solanum lycopersicum*VegetableCRISPR*Carotenoid isomerase and Psy1*Carotenoid metabolismMetabolism^[@CR117]^*Solanum lycopersicum*VegetableCRISPR*Solyc08g075770*Fusarium wilt susceptibilityStress response^[@CR118]^*Solanum lycopersicum*VegetableCRISPR*TypeII GRX 14, 15, 16, 17*Redox regulationMetabolism^[@CR119]^*Solanum lycopersicum*VegetableCRISPR*lncRNA1459*Repressed fruit ripening, lycopene, ethylene and carotenoid biosynthesisMetabolism^[@CR120]^*Solanum lycopersicum*VegetableCRISPR*SGR1, Blc, LCY-E, LCY-B1, LCY-B2*Increased lycopene contentMetabolism^[@CR121]^*Solanum lycopersicum*VegetableCRISPR*PDS*Albino phenotypeReporter^[@CR122]^*Solanum lycopersicum*VegetableCRISPR*SlDML2*DNA methylation and fruit ripeningReporter^[@CR66]^*Solanum lycopersicum*VegetableCRISPR*PDS and GABA-TP1, GABA-TP2, GABA-TP3, CAT9 and SSADH*γ-aminobutyric acid metabolismMetabolism^[@CR123]^*Solanum lycopersicum*VegetableCRISPR*SlMYB12*Pink tomato fruit colorMetabolism^[@CR124]^*Solanum lycopersicum*VegetableCRISPR*Coat protein, Replicase from TYLCV*Obtained resistance to tomato yellow leaf curl virusStress response^[@CR125]^*Solanum lycopersicum*VegetableCRISPR*RIN*Ethylene production and fruit ripeningMetabolism^[@CR126]^*Solanum pimpinellifolium*VegetableCRISPR*SP, MULT, FAS, CyCb, OVUTE and FW2.2*Plant and inflorescence architecture, fruit shape and lycopene biosynthesisDevelopment, metabolism^[@CR69]^*Solanum pimpinellifolium*VegetableCRISPR*SP, SP5, CLV3 and WUS, GGP1*plant architecture, day-length insensitivity, enlarged fruit size and vitamin CDevelopment, metabolism^[@CR70]^*Solanum lycopersicum*VegetableCRISPR*RIN*Ethylene production and fruit ripeningDevelopment^[@CR68]^*Solanum lycopersicum*VegetableCRISPR*SlORRM4*RNA editing and fruit ripeningDevelopment^[@CR67]^*Solanum lycopersicum*VegetableCRISPR*ALC*Shelf lifeMetabolism^[@CR127]^*Solanum lycopersicum*VegetableCRISPR*CLAVATA-WUSCHEL*Altered locule numberDevelopment^[@CR65]^*Solanum lycopersicum*VegetableCRISPR*SlMAPK3*Drought stressStress response^[@CR128]^*Solanum lycopersicum*VegetableCRISPR*Glutamate decarboxylase (GAD)*γ-aminobutyric acid metabolismMetabolism^[@CR129]^*Solanum lycopersicum*VegetableCRISPR*Solyc12g038510*Jointless mutant, abscissionDevelopment^[@CR130]^*Solanum lycopersicum*VegetableCRISPR*Multiple genes*Generate a pool of mutantsOthers^[@CR131]^*Solanum lycopersicum*VegetableCRISPR*PSY*Fruit colorDevelopment^[@CR132]^*Solanum lycopersicum*VegetableCRISPR*Solyc12g038510*Jointless and branchingDevelopment^[@CR133]^*Solanum lycopersicum*VegetableCRISPR*L1L4*Involved in fruit metabolism during ripeningMetabolism^[@CR134]^*Solanum lycopersicum*VegetableCRISPR*DELLA and ETR*Hormone responseDevelopment^[@CR135]^*Solanum lycopersicum*VegetableCRISPR*SlMlo1*Powdery mildew resistanceStress response^[@CR136]^*Solanum lycopersicum*VegetableCRISPR*SlIAA9*Parthenocarpic tomato plantsDevelopment^[@CR137]^*Solanum lycopersicum*VegetableCRISPR*SP5G*More rapid floweringDevelopment^[@CR64]^*Solanum lycopersicum*VegetableCRISPR*Genes involved tomato domestication*Development and plant architectureDevelopment^[@CR138]^*Solanum lycopersicum*VegetableCRISPR*SlAGL6*Production of parthenocarpic fruit under high temperatureDevelopment^[@CR139]^*Solanum lycopersicum*VegetableCRISPR*N.A*N.AOthers^[@CR140]^*Solanum lycopersicum*VegetableCRISPR*SlBOP*Inflorescence structureDevelopment^[@CR63]^*Solanum lycopersicum*VegetableZFN*L1L4*Heterochronic phenotype, plant architectureDevelopment^[@CR141]^*Solanum lycopersicum*VegetableCRISPR*PDS and PIF*Albino phenotypeReporter^[@CR142]^*Solanum lycopersicum*VegetableCRISPR*N.A*.N.A.Others^[@CR143]^*Solanum lycopersicum*VegetableTALEN, CRISPR*ANT1*Anthocyanin biosynthesisMetabolism^[@CR144]^*Solanum lycopersicum*VegetableCRISPR*RIN*Fruit ripeningDevelopment^[@CR145]^*Solanum lycopersicum*VegetableTALEN*PROCERA*GA response and taller plantDevelopment^[@CR146]^*Solanum lycopersicum*VegetableCRISPR*AGO7*Leaf morphologyDevelopment^[@CR147]^*Solanum* *tuberosum*VegetableCRISPR*St16DOX*Steroidal glycoalkaloids metabolismMetabolism^[@CR148]^*Solanum* *tuberosum*VegetableCRISPR*GBSS genes*Starch biosynthesisMetabolism^[@CR149]^*Solanum* *tuberosum*VegetableCRISPR*S-RNase*Self-incompatibilityDevelopment^[@CR150]^*Solanum* *tuberosum*VegetableCRISPR*Coilin gene*Enhanced resistance to biotic and abiotic agentsStress response^[@CR151]^*Solanum* *tuberosum*VegetableCRISPR*StALS1, StALS2*Enhanced herbicide resistanceStress response^[@CR101]^*Solanum* *tuberosum*VegetableCRISPR*GBSS1*Starch biosynthesisMetabolism^[@CR152]^*Solanum* *tuberosum*VegetableCRISPR*S-Rnase*Self-incompatibilityDevelopment^[@CR72]^*Solanum* *tuberosum*VegetableCRISPR*Coilin gene*Enhanced resistance to biotic and abiotic agentsStress response^[@CR153]^*Solanum* *tuberosum*VegetableTALEN*SBE1 and StvacINV22*Sugar metabolismMetabolism^[@CR154]^*Solanum* *tuberosum*VegetableCRISPR*StMYB44*Phosphorus homeostasisStress response^[@CR155]^*Solanum* *tuberosum*VegetableCRISPR*GBSS*Starch metabolism and tuber qualityMetabolism^[@CR156]^*Solanum* *tuberosum*VegetableTALEN*StALS1*Enhanced herbicide resistanceStress response^[@CR157]^*Solanum* *tuberosum*VegetableTALEN*StALS1*Enhanced herbicide resistanceStress response^[@CR158]^*Solanum* *tuberosum*VegetableTALEN*vINV*Postharvest cold storage and processingMetabolism^[@CR71]^*Solanum* *tuberosum*VegetableCRISPR*StALS1*Enhanced herbicide resistanceMetabolism^[@CR159]^*Solanum* *tuberosum*VegetableTALEN*StALS1*Enhanced herbicide resistanceMetabolism^[@CR160]^*Solanum* *tuberosum*VegetableCRISPR*StIAA2*Aux/IAA protein, shoot morphogenesisDevelopment^[@CR161]^*Brassica oleracea*VegetableCRISPR*BolC.GA4.a*GA response and dwarfismDevelopment^[@CR162]^*Brassica oleracea*VegetableCRISPR*BoPDS, BoSRK3, BoMS1*Albino phenotype, self-incompatibility, male sterilityDevelopment^[@CR163]^*Brassica napus*VegetableCRISPR*LMI1*Leaf lobe developmentDevelopment^[@CR164]^*Brassica oleracea, rapa*VegetableCRISPR*PDS and FRI*Albino phenotype and floweringReporter, development^[@CR165]^*Brassica napus*VegetableCRISPR*FAD2*Fatty acid metabolismMetabolism^[@CR76]^*Brassica carinata*VegetableCRISPR*Fascilin-like arabinogalactan protein*Regulation of root hairs under phosphorus stressDevelopment, stress response^[@CR166]^*Brassica napus*VegetableCRISPR*WRKY11 and WRKY70*Enhanced biotic resistanceStress response^[@CR167]^*Brassica napus*VegetableCRISPR*SDG8*Histone lysine methyltransferaseDevelopment^[@CR168]^*Brassica napus*VegetableCRISPR*CLV3 and CLV1, CLV2*Regulate multilocular seedsDevelopment^[@CR169]^*Brassica rapa and napus*VegetableCRISPR*AP2a, AP2b*Sepal to carpal modificationDevelopment^[@CR170]^*Brassica napus*VegetableCRISPR*BnaRGA, BnaDA1, BnaDA2, BnaFUL*Multiple genes involved in plant developmentDevelopment^[@CR171]^*Brassica carinata*VegetableCRISPR*Fascilin-like arabinogalactan protein*Root hair developmentDevelopment^[@CR172]^*Brassica napus*VegetableCRISPR*ALC*Valve margin development, seed shatteringDevelopment^[@CR173]^*Brassica oleracea*VegetableTALEN*FRIGIDA*Early flowering phenotypeDevelopment^[@CR62]^*Dendrobium officinale*FlowerCRISPR*C3H, C4H, 4CL, CCR, and IRX*Lignocellulose biosynthesisMetabolism^[@CR174]^*Lettuce sativa*VegetableCRISPR*LsBIN2*Impaired brassinosteroid responseDevelopment^[@CR83]^*Lettuce sativa*VegetableCRISPR*LsNCED4*Thermo-inhibition of seed germinationDevelopment^[@CR175]^*Cucumis sativus*VegetableCRISPR*eIF4E*Enhanced viral resistanceStress response^[@CR176]^*Cucumis sativus*VegetableCRISPR*CmWIP1*Gynoecious phenotypeDevelopment^[@CR177]^*Musa balbisiana*FruitCRISPR*eBSV*Control of virus pathogenesisStress response^[@CR178]^*Musa acuminata*FruitCRISPR*PDS*Albino phenotypeReporter^[@CR179]^*Musa acuminata*FruitCRISPR*PDS*Albino phenotypeReporter^[@CR180]^*Actinidia deliciosa*FruitCRISPR*PDS*Albino phenotypeReporter^[@CR181]^*Vitis vinifera*FruitCRISPR*VvPDS*Albino phenotypeReporter^[@CR182]^*Vitis vinifera*FruitCRISPR*IdnDH*Biosynthesis of tartaric acidMetabolism^[@CR183]^*Vitis vinifera*FruitCRISPR*VvWRKY52*Increased the resistance to *Botrytis cinerea*Stress response^[@CR75]^*Vitis vinifera*FruitCRISPR*VvPDS*Albino phenotypeReporter^[@CR184]^*Vitis vinifera*FruitCRISPR*MLO-7*Powdery mildew resistanceStress response^[@CR185]^*Vitis vinifera*FruitCRISPR*IdnDH*Biosynthesis of tartaric acidMetabolism^[@CR186]^*Citrus sinensis*FruitCRISPR*DMR6*Huanglongbin resistanceStress response^[@CR187]^*Citrus sinensis*FruitCRISPR*PDS*Albino phenotypeReporter^[@CR188]^*Citrus paradisi*FruitCRISPR*CsPDS, Cs2g12470 and Cs7g03360*Albino phenotypeReporter^[@CR189]^*Citrus sinensis*FruitCRISPR*PDS*Albino phenotypeReporter^[@CR190]^*Citrus sinensis*FruitCRISPR*CsLOB1*Canker resistanceStress response^[@CR73]^*Citrus paradisi*FruitCRISPR*CsLOB1*Canker resistanceStress response^[@CR74]^*Citrus sinensis*FruitCRISPR*CsPDS*Albino phenotypeReporter^[@CR191]^*Chrysanthemum morifolium*FlowerCRISPR*CpYGFP*Targeted editing of the YGFP reporter geneOthers^[@CR192]^*Ipomoea nil*FlowerCRISPR*InDFR-B*Anthocyanin biosynthesis and white flowersMetabolism^[@CR193]^*Ipomoea nil*FlowerCRISPR*InCCD4*Altered petal colorDevelopment^[@CR194]^*Petunia inflata*FlowerCRISPR*PiSSK1*Self-incompatibilityDevelopment^[@CR195]^*Petunia hybrid*FlowerCRISPR*PDS*Albino phenotypeReporter^[@CR196]^*Citrullus lanatus*FruitCRISPR*ALS*Increased herbicide resistanceStress response^[@CR197]^*Citrullus lanatus*FruitCRISPR*PDS*Albino phenotypeReporter^[@CR198]^*Salvia miltiorrhiza*Medicinal plantCRISPR*SmCPS1*Tanshinone biosynthesisMetabolism^[@CR199]^*Camelina sativa*VegetableCRISPR*FAE1*Reduced long-chain fatty acidsMetabolism^[@CR77]^*Camelina sativa*VegetableCRISPR*CsDGAT1 or CsPDAT1*Altered fatty acid composition and reduced oil contentMetabolism^[@CR200]^*Camelina sativa*VegetableCRISPR*FAD2*Reduced levels of polyunsaturated fatty acidsMetabolism^[@CR78]^*Camelina sativa*VegetableCRISPR*FAD2*Decreased polyunsaturated fatty acidsMetabolism^[@CR79]^Malus pumilaFruitCRISPR*PDS, TFL1.1*Albino phenotype, early floweringDevelopment^[@CR201]^*Malus* *pumila*FruitCRISPR*PDS*Albino phenotypeReporter^[@CR183]^*Malus* *pumila*FruitCRISPR*PDS*Albino phenotypeReporter^[@CR202]^*Malus* *pumila*FruitCRISPR*DIPM*Blight resistanceStress response^[@CR185]^*Malus* *pumila*FruitZFN*udiA*Edited reporter geneOthers^[@CR203]^*Pyrus communis*FruitCRISPR*TFL1.1*Early floweringDevelopment^[@CR201]^*Daucus carota*VegetableCRISPR*PDS, MYB113-like*Albino phenotypeReporter^[@CR204]^*Daucus carota*VegetableCRISPR*F3H*Altered anthocyanin biosynthesisMetabolism^[@CR205]^*Torenia fournieri*FlowerCRISPR*F3H*Altered flower pigmentationMetabolism^[@CR206]^*Fragaria vesca*FruitCRISPR*FveTAA1, FveARF8*Auxin signaling, plant developmentDevelopment^[@CR207]^*Fragaria vesca, Fragaria x Ananassa*FruitCRISPR*FvMYB10, FvCHS*Anthocyanin biosynthesisMetabolism^[@CR208]^*Fragaria x Ananassa*FruitCRISPR*FaTM6*Anther developmentDevelopment^[@CR209]^*Fragaria vesca, Fragaria x Ananassa*FruitCRISPR*PDS*Albino phenotypeReporter^[@CR210],\ [@CR211]^

In tomato, development-related genes have been edited to manipulate flowering patterns and fruit development. The tomato *BLADE-ON-PETIOLE* (*BOP*) genes, which encode transcriptional cofactors, can regulate inflorescence structure, and knock-out of *SlBOP* genes by gene editing reduces the number of flowers per inflorescence^[@CR63]^. CRISPR/Cas9-induced mutations in the flowering repressor *self-pruning 5G* lead to rapid flowering and early harvest^[@CR64]^. In addition, editing of the cis-regulatory region of *SlCLV3*^[@CR65]^ or the coding regions of *SlDML2*^[@CR66]^, *SlORRM4*^[@CR67]^ and the *RIN* locus^[@CR68]^ alters fruit development and ripening. Interestingly, multiplex targeting of several genes that are important for tomato domestication was found to greatly alter the properties of the wild tomato relative *Solanum pimpinellifolium* such that the generated mutants were similar to cultivated tomato^[@CR69],[@CR70]^. In potato, when the vacuolar invertase gene was disrupted by TALEN, the cold storage and processing of tubers were improved^[@CR71]^. Another recent study in potato showed the possibility of overcoming self-incompatibility by editing the *S-RNase* gene, which would provide an alternative method of propagation through seeds^[@CR72]^. In addition to tomato and potato, other horticultural crops have also been edited to obtain desirable traits. Genes related to resistance to plant pathogens such as *Xanthomonas citri*^[@CR73],[@CR74]^ and *Botrytis cinerea*^[@CR75]^ have been manipulated in citrus, apple, and grape. In oilseed crops, genes involved in fatty acid metabolism have been frequently targeted to improve oil quality^[@CR76]--[@CR79]^. The application of genome editing to improve crops is based on knowledge of the association between genes and their controlled traits. In the future, functional characterization of genes in different crops will help to identify valuable targets that could be edited for potential horticultural improvement, such as increased productivity, marketing quality, and nutritional value.

Possible implications of genome editing in horticulture {#Sec7}
=======================================================

The goal of breeding is to harness genetic variations to introduce desirable traits. These genetic variations can arise in various ways, such as by spontaneous mutation, chemical mutagenesis, and physical mutagenesis. Gene editing could be regarded as biological mutagenesis. In comparison with other approaches, genome-editing technology is superior in terms of versatility, efficiency, and specificity. For instance, CRISPR-based genome editing can cause many types of mutations in target sequences, including small insertions/deletions, deletions of large fragments, gene replacement, and precise base substitutions^[@CR16]^. In addition, genome-editing technology is continuously advancing: the endonuclease Cpf1^[@CR51]^ and newly discovered or designed Cas9 variants^[@CR80],[@CR81]^ can recognize different PAM sequences, thereby broadening the genome-wide sites that can be targeted for editing.

Genome-edited plants are not considered genetically modified organisms (GMOs) in countries such as the U.S. and Japan but are still under strict GMO regulation in Europe. The largest difference between genome-edited plants and GMOs is that the genomes of edited plants can be free of exogenous DNA sequences. The exogenous DNA of the editing tools can be removed through genetic segregation^[@CR82]^ or may never have to be introduced if CRISPR reagents are delivered as ribonucleoproteins^[@CR83],[@CR84]^.

Mutants generated via genome editing can be directly used for crop production or as prebreeding materials. Through genome editing, desirable traits can be directly introgressed into elite or heirloom lines without compromising other properties, and the resulting lines with targeted improvement will be ready for use in production. The wild relatives of cultivated varieties are also potential materials for genome editing because they generally present unique features in many important traits. For instance, wild species of cultivated tomato are more resistant to unfavorable environments than commercial cultivars^[@CR85]^. Wild *Solanum pimpinellifolium* was recently domesticated by the editing of several important genes affecting plant architecture and fruit development, resulting in new tomato varieties with the desirable properties of cultivated tomato combined with the favorable traits of the wild species^[@CR69],[@CR70]^. Mutations can generally be introduced in either the coding region or the cis-regulatory region of the targeted gene, and mutations in the cis-regulatory region could be used to generate quantitative variation for breeding selection. In tomato, for example, fruit locule number is determined by several naturally occurring mutations in the cis-regulatory regions of *CLAVATA-WUSCHEL*^[@CR65]^. This finding motivated researchers to design a multiplexed CRISPR/Cas9 system targeting the *CLAVATA-WUSCHEL* promoters to generate tomato lines with a wide range of locule numbers. Quantitative variations have also been observed when the genes responsible for inflorescence and plant architecture are engineered^[@CR65]^. In addition to regulating gene activity by editing the DNA sequence of the cis-regulatory region, gene activity can be regulated by the its epigenetic status of this region. By integrating genome editing (CRISPR/Cas9) with epigenetic regulation, researchers are able to target a gene of interest and modify its epigenetic status. For instance, an sgRNA-guided fusion protein between the dead Cas9 (dCas9) variant and the catalytic domain of the TEN-ELEVEN TRANSLOCATION1 (TET1cd) demethylase can remove 5mC at specific sites, thereby increasing gene expression^[@CR86]^. An epigenetic mutant can also be crossed with the corresponding wild type to generate epigenetic recombinant inbred lines (epiRILs). Individuals from these populations are genetically identical but epigenetically distinct. Such populations have been constructed in Arabidopsis and exhibit considerable phenotypic variations^[@CR87]--[@CR90]^. These examples demonstrate that genome editing is an excellent tool for producing new alleles and epialleles, which are important sources of phenotypic variation for crop improvement.

Challenges and future perspectives for the improvement of horticultural crops through genome editing {#Sec8}
====================================================================================================

Although genome editing has many advantages over conventional crop breeding, some challenges remain for its application to horticultural crops. In horticultural crops, molecular and genetic studies are difficult, which hinders the identification of genes responsible for desirable traits. Sequencing the genomes of horticultural crops of interest will be important for identifying genes associated with desirable traits. For crops lacking a reference genome, the target sequence could be cloned by using degenerate primers designed for conserved protein motifs with putative functions related to desirable traits. A good example is the *mildew-resistance locus* (*MLO*), which has been characterized in detail in barley^[@CR91]^; the phylogenetically conservative nature of the *MLO* has facilitated the generation of powdery mildew-resistant plants in wheat, tomato, and strawberry^[@CR92],[@CR93]^.

Once a gene to be edited has been identified, researchers must take into account the methods used to deliver editing reagents and the procedure for regenerating the edited mutants. To date, more than 25 horticultural plant species have been successfully edited (Table [2](#Tab2){ref-type="table"}), usually with editing reagents delivered via *Agrobacteria* or virus systems, and the edited plants are regenerated via in vitro tissue culture. Although tissue culture-based transformation and regeneration is most widely used for genome editing, no well-established protocol for transformation and regeneration from tissue culture is available for many horticultural crops. *In planta* transformation, which is an alternative to in vitro tissue culture-based *Agrobacterium* transformation, refers to the infection of in vivo explants in which the targeted tissues are apical or auxiliary meristems, stigmas, pollens, or inflorescences^[@CR94]^. This method has been successfully used to transform tomato^[@CR95]^ and *Brassica* species^[@CR96]^ and should be further explored for use in horticultural crops that are recalcitrant to traditional genetic transformation. Additionally, successful genetic transformation of horticultural crops requires the consideration of editing efficiency, which is affected by many factors, such as sgRNA number and GC content, the expression levels of sgRNA and Cas9, and the secondary structure of the paired sgRNA and target sequence^[@CR97],[@CR98]^. In the future, the editing system should be further optimized in different crop species.

The elimination of foreign DNA fragments (transferred T-DNAs) to obtain transgene-free edited plants remains difficult in some highly heterozygous and clonally propagated horticultural species^[@CR99]^, such as potato, sweet potato, and banana. One possibility is to generate many transformants, followed by high-throughput screening of transgene-free mutants^[@CR100]^. This approach has been used to generate \~10% of mutants without foreign DNA^[@CR100],[@CR101]^. Another approach for transgene-free genome editing is to deliver editing reagents as in vitro transcripts^[@CR102]^ or ribonucleoproteins^[@CR83],[@CR84]^.

In conclusion, mutagenesis via genome editing outperforms spontaneous and induced mutations in terms of precision and efficiency. Although this technology is being increasingly used in many crops, its widespread use in the breeding of horticultural crops will require three challenges to be surmounted. First, clear breeding traits of the horticultural crop in question should be identified via communication among consumers, breeders, and biologists. Second and third, suitable methods must be developed for delivering editing reagents and for subsequently regenerating mutants. Given the great potential of genome editing and the importance of horticultural crops, we expect that these challenges will be overcome in the near future.

The authors declare that they have no conflict of interest.
