Corrigendum: Transcriptomic analysis in tomato fruit reveals divergences in genes involved in cold stress response and fruit ripening.

[This corrects the article DOI: 10.3389/fpls.2023.1227349.].

Class II knotted-like homeodomain protein SlKN5 with BEL1-like homeodomain proteins suppresses fruit greening in tomato fruit.

Knotted1-like homeodomain (KNOX) proteins are essential in regulating plant organ differentiation. Land plants, including tomato (Solanum lycopersicum), have two classes of the KNOX protein family, namely, class I (KNOX I) and class II KNOX (KNOX II). While tomato KNOX I proteins are known to stimulate chloroplast development in fruit, affecting fruit coloration, the role of KNOX II proteins in this context remains unclear. In this study, we employ CRISPR/Cas9 to generate knockout mutants of the KNOX II member, SlKN5. These mutants display increased leaf complexity, a phenotype commonly associated with reduced KNOX II activity, as well as enhanced accumulation of chloroplasts and chlorophylls in smaller cells within young, unripe fruit. RNA-seq data analyses indicate that SlKN5 suppresses the transcriptions of genes involved in chloroplast biogenesis, chlorophyll biosynthesis, and gibberellin catabolism. Furthermore, protein-protein interaction assays reveal that SlKN5 physically interacts with three transcriptional repressors from the BLH1-clade of BEL1-like homeodomain (BLH) protein family, SlBLH4, SlBLH5, and SlBLH7, with SlBLH7 showing the strongest interaction. CRISPR/Cas9-mediated knockout of these SlBLH genes confirmed their overlapping roles in suppressing chloroplast biogenesis, chlorophyll biosynthesis, and lycopene cyclization. Transient assays further demonstrate that the SlKN5-SlBLH7 interaction enhances binding capacity to regulatory regions of key chloroplast- and chlorophyll-related genes, including SlAPRR2-like1, SlCAB-1C, and SlGUN4. Collectively, our findings elucidate that the KNOX II SlKN5-SlBLH regulatory modules serve to inhibit fruit greening and subsequently promote lycopene accumulation, thereby fine-tuning the color transition from immature green fruit to mature red fruit.

Transcriptomic analysis in tomato fruit reveals divergences in genes involved in cold stress response and fruit ripening.

Cold storage is widely used to extend the postharvest life of most horticultural crops, including tomatoes, but this practice triggers cold stress and leads to the development of undesirable chilling injury (CI) symptoms. The underlying mechanisms of cold stress response and CI development in fruits remain unclear as they are often intermingled with fruit ripening changes. To gain insight into cold responses in fruits, we examined the effect of the potent ethylene signaling inhibitor 1-methylcyclopropene (1-MCP) on fruit ripening, CI occurrence and gene expression in mature green tomatoes during storage at 20°C and 5°C. 1-MCP treatments effectively inhibited ethylene production and peel color changes during storage at 20°C. Storage at 5°C also inhibited both ethylene production and peel color change; during rewarming at 20°C, 1-MCP treatments inhibited peel color change but failed to inhibit ethylene production. Furthermore, fruits stored at 5°C for 14 d developed CI symptoms (surface pitting and decay) during the rewarming period at 20°C regardless of 1-MCP treatment. Subsequent RNA-Seq analysis revealed that cold stress triggers a large-scale transcriptomic adjustment, as noticeably more genes were differentially expressed at 5°C (8,406) than at 20°C (4,814). More importantly, we have found some important divergences among genes involved in fruit ripening (up- or down-regulated at 20°C; inhibited by 1-MCP treatment) and those involved in cold stress (up- or down-regulated at 5°C; unaffected by 1-MCP treatment). Transcriptomic adjustments unique to cold stress response were associated with ribosome biogenesis, NcRNA metabolism, DNA methylation, chromatin formation/remodeling, and alternative splicing events. These data should foster further research into cold stress response mechanisms in fruits with the ultimate aim of improving tolerance to low temperature and reduction of CI symptoms during cold storage.

Molecular, hormonal, and metabolic mechanisms of fruit set, the ovary-to-fruit transition, in horticultural crops.

Fruit set is the process by which the ovary develops into a fruit and is an important factor in determining fruit yield. Fruit set is induced by two hormones, auxin and gibberellin, and the activation of their signaling pathways, partly by suppressing various negative regulators. Many studies have investigated the structural changes and gene networks in the ovary during fruit set, revealing the cytological and molecular mechanisms. In tomato (Solanum lycopersicum), SlIAA9 and SlDELLA/PROCERA act as auxin and gibberellin signaling repressors, respectively, and are important regulators of the activity of transcription factors and downstream gene expression involved in fruit set. Upon pollination, SlIAA9 and SlDELLA are degraded, which subsequently activates downstream cascades and mainly contributes to active cell division and cell elongation, respectively, in ovaries during fruit setting. According to current knowledge, the gibberellin pathway functions as the most downstream signal in fruit set induction, and therefore its role in fruit set has been extensively explored. Furthermore, multi-omics analysis has revealed the detailed dynamics of gene expression and metabolites downstream of gibberellins, highlighting the rapid activation of central carbon metabolism. This review will outline the relevant mechanisms at the molecular and metabolic levels during fruit set, particularly focusing on tomato.

Phenotypic Characterization of High Carotenoid Tomato Mutants Generated by the Target-AID Base-Editing Technology.

Our previous study demonstrated that Target-AID which is the modified CRISPR/Cas9 system enabling base-editing is an efficient tool for targeting multiple genes. Three genes, SlDDB1, SlDET1, and SlCYC-B, responsible for carotenoid accumulation were targeted, and allelic variations were previously obtained by Target-AID. In this research, we characterized the effect of new alleles on plant growth and fruit development, as well as carotenoid accumulation, individually in segregating backcross populations or combined in null self-segregant lines. Only lines carrying homozygous substitutions in the three targeted genes and the segregating backcross population of individual mutations were characterized, resulting in the isolation of two allelic versions for SlDDB1, one associated with SlDET1 and the last one with SlCYC-B. All edited lines showed variations in carotenoid accumulation, with an additive effect for each single mutation. These results suggest that Target-AID base-editing technology is an effective tool for creating new allelic variations in target genes to improve carotenoid accumulation in tomato.

Transcriptomic, Hormonomic and Metabolomic Analyses Highlighted the Common Modules Related to Photosynthesis, Sugar Metabolism and Cell Division in Parthenocarpic Tomato Fruits during Early Fruit Set.

Parthenocarpy, the pollination-independent fruit set, can raise the productivity of the fruit set even under adverse factors during the reproductive phase. The application of plant hormones stimulates parthenocarpy, but artificial hormones incur extra financial and labour costs to farmers and can induce the formation of deformed fruit. This study examines the performance of parthenocarpic mutants having no transcription factors of SlIAA9 and SlTAP3 and sldella that do not have the protein-coding gene, SlDELLA, in tomato (cv. Micro-Tom). At 0 day after the flowering (DAF) stage and DAFs after pollination, the sliaa9 mutant demonstrated increased pistil development compared to the other two mutants and wild type (WT). In contrast to WT and the other mutants, the sliaa9 mutant with pollination efficiently stimulated the build-up of auxin and GAs after flowering. Alterations in both transcript and metabolite profiles existed for WT with and without pollination, while the three mutants without pollination demonstrated the comparable metabolomic status of pollinated WT. Network analysis showed key modules linked to photosynthesis, sugar metabolism and cell proliferation. Equivalent modules were noticed in the famous parthenocarpic cultivars 'Severianin', particularly for emasculated samples. Our discovery indicates that controlling the genes and metabolites proffers future breeding policies for tomatoes.

Genome-wide cis-decoding for expression design in tomato using cistrome data and explainable deep learning.

In the evolutionary history of plants, variation in cis-regulatory elements (CREs) resulting in diversification of gene expression has played a central role in driving the evolution of lineage-specific traits. However, it is difficult to predict expression behaviors from CRE patterns to properly harness them, mainly because the biological processes are complex. In this study, we used cistrome datasets and explainable convolutional neural network (CNN) frameworks to predict genome-wide expression patterns in tomato (Solanum lycopersicum) fruit from the DNA sequences in gene regulatory regions. By fixing the effects of trans-acting factors using single cell-type spatiotemporal transcriptome data for the response variables, we developed a prediction model for crucial expression patterns in the initiation of tomato fruit ripening. Feature visualization of the CNNs identified nucleotide residues critical to the objective expression pattern in each gene, and their effects were validated experimentally in ripening tomato fruit. This cis-decoding framework will not only contribute to the understanding of the regulatory networks derived from CREs and transcription factor interactions, but also provides a flexible means of designing alleles for optimized expression.

Modification of tomato breeding traits and plant hormone signaling by target-AID, the genome-editing system inducing efficient nucleotide substitution.

Target activation-induced cytidine deaminase (Target-AID), a novel CRISPR/Cas9-based genome-editing tool, confers the base-editing capability on the Cas9 genome-editing system. It involves the fusion of cytidine deaminase (CDA), which catalyzes cytidine (C) to uridine (U) substitutions, to the mutated nickase-type nCas9 or deactivated-type dCas9. To confirm and extend the applicability of the Target-AID genome-editing system in tomatoes (Solanum lycopersicum L.), we transformed the model tomato cultivar "Micro-Tom" and commercial tomato cultivars using this system by targeting SlDELLA, which encodes a negative regulator of the plant phytohormone gibberellic acid (GA) signaling pathway. We confirmed that the nucleotide substitutions were induced by the Target-AID system, and we isolated mutants showing high GA sensitivity in both "Micro-Tom" and the commercial cultivars. Moreover, by successfully applying this system to ETHYLENE RECEPTOR 1 (SlETR1) with single sgRNA targeting, double sgRNA targeting, as well as dual-targeting of both SlETR1 and SlETR2 with a single sgRNA, we demonstrated that the Target-AID genome-editing system is a promising tool for molecular breeding in tomato crops. This study highlights an important aspect of the scientific and agricultural potential of the combinatorial use of the Target-AID and other base-editing systems.

Organelle genome assembly uncovers the dynamic genome reorganization and cytoplasmic male sterility associated genes in tomato.

To identify cytoplasmic male sterility (CMS)-associated genes in tomato, we determined the genome sequences of mitochondria and chloroplasts in three CMS tomato lines derived from independent asymmetric cell fusions, their nuclear and cytoplasmic donors, and male fertile weedy cultivated tomato and wild relatives. The structures of the CMS mitochondrial genomes were highly divergent from those of the nuclear and cytoplasmic donors, and genes of the donors were mixed up in these genomes. On the other hand, the structures of CMS chloroplast genomes were moderately conserved across the donors, but CMS chloroplast genes were unexpectedly likely derived from the nuclear donors. Comparative analysis of the structures and contents of organelle genes and transcriptome analysis identified three genes that were uniquely present in the CMS lines, but not in the donor or fertile lines. RNA-sequencing analysis indicated that these three genes transcriptionally expressed in anther, and identified different RNA editing levels in one gene, orf265, that was partially similar to ATP synthase subunit 8, between fertile and sterile lines. The orf265 was a highly potential candidate for CMS-associated gene. This study suggests that organelle reorganization mechanisms after cell fusion events differ between mitochondria and chloroplasts, and provides insight into the development of new F1 hybrid breeding programs employing the CMS system in tomato.

Functional disruption of cell wall invertase inhibitor by genome editing increases sugar content of tomato fruit without decrease fruit weight.

Sugar content is one of the most important quality traits of tomato. Cell wall invertase promotes sucrose unloading in the fruit by maintaining a gradient of sucrose concentration between source leaves and fruits, while invertase inhibitor (INVINH) regulates this process. In this study, knock-out of cell wall INVINH in tomato (SlINVINH1) was performed by genome editing using, CRISPR/Cas9 and Target-AID technologies. Most of the genome-edited lines set higher soluble solid content (SSC) fruit than the original cultivar 'Suzukoma', while fruit weight was different among the genome-edited lines. From these genome-edited lines, three lines (193-3, 199-2, and 247-2), whose SSC was significantly higher than 'Suzukoma' and fruit weight were almost the same as the original cultivar, were selected. The fruit weight and overall plant growth of the two lines were comparable to those of the original cultivar. In contrast, the fructose and glucose contents in the mature fruits of the two lines were significantly higher than those of the original cultivar. The mature fruits of genome edited line 193-3 showed the highest sugar content, and the fructose and glucose contents were 29% and 36% higher than that of the original cultivar, respectively. Whole genome sequence data showed no off-target mutations in the genome-edited lines. Non-target metabolome analysis of mature fruits revealed that fructose was the highest loading factor in principal component analysis (PCA) between the genome-edited line and the original cultivar, and no unexpected metabolites appeared in the genome-edited line. In this study, we succeeded in producing tomato lines with high sugar content without a decrease in fruit weight and deterioration of plant growth by knock-out of SlINVINH1 using genome editing technology. This study showed that functional disruption of SlINVINH1 is an effective approach to produce tomato cultivars with high sugar content.

Corrigendum: Loss-of-Function of a Tomato Receptor-Like Kinase Impairs Male Fertility and Induces Parthenocarpic Fruit Set.

[This corrects the article DOI: 10.3389/fpls.2019.00403.].

Overproduction of ascorbic acid impairs pollen fertility in tomato.

Ascorbate is a major antioxidant buffer in plants. Several approaches have been used to increase the ascorbate content of fruits and vegetables. Here, we combined forward genetics with mapping-by-sequencing approaches using an ethyl methanesulfonate (EMS)-mutagenized Micro-Tom population to identify putative regulators underlying a high-ascorbate phenotype in tomato fruits. Among the ascorbate-enriched mutants, the family with the highest fruit ascorbate level (P17C5, up to 5-fold wild-type level) had strongly impaired flower development and produced seedless fruit. Genetic characterization was performed by outcrossing P17C5 with cv. M82. We identified the mutation responsible for the ascorbate-enriched trait in a cis-acting upstream open reading frame (uORF) involved in the downstream regulation of GDP-l-galactose phosphorylase (GGP). Using a specific CRISPR strategy, we generated uORF-GGP1 mutants and confirmed the ascorbate-enriched phenotype. We further investigated the impact of the ascorbate-enriched trait in tomato plants by phenotyping the original P17C5 EMS mutant, the population of outcrossed P17C5 × M82 plants, and the CRISPR-mutated line. These studies revealed that high ascorbate content is linked to impaired floral organ architecture, particularly anther and pollen development, leading to male sterility. RNA-seq analysis suggested that uORF-GGP1 acts as a regulator of ascorbate synthesis that maintains redox homeostasis to allow appropriate plant development.

Efficient base editing in tomato using a highly expressed transient system.

KEY MESSAGE: Base editing in tomatoes was achieved by transient expression. The Solanaceae plants, particularly the tomato (Solanum lycopersicum), is of huge economic value worldwide. The tomato is a unique model plant for studying the functions of genes related to fruit ripening. Deeper understanding of tomatoes is of great importance for both plant research and the economy. Genome editing technology, such as CRISPR/Cas9, has been used for functional genetic research. However, some challenges, such as low transformation efficiency, remain with this technology. Moreover, the foreign Cas9 and gRNA expression cassettes must be removed to obtain null-segregants In this study, we used a high-level transient expression system to improve the base editing technology. A high-level transient expression system has been established previously using geminiviral replication and a double terminator. The pBYR2HS vector was used for this transient expression system. nCas9-CDA and sgRNA-SlHWS were introduced into this vector, and the protein and RNA were then transiently expressed in tomato tissues by agroinfiltration. The homozygous mutant produced by base editing was obtained in the next generation with an efficiency of about 18%. nCas9-free next-generation plants were 71%. All the homozygous base-edited plants in next generation are nCas9-free. These findings show that the high-level transient expression system is useful for base editing in tomatoes.

De novo genome assembly of two tomato ancestors, Solanum pimpinellifolium and Solanum lycopersicum var. cerasiforme, by long-read sequencing.

The ancestral tomato species are known to possess genes that are valuable for improving traits in breeding. Here, we aimed to construct high-quality de novo genome assemblies of Solanum pimpinellifolium 'LA1670' and S. lycopersicum var. cerasiforme 'LA1673', originating from Peru. The Pacific Biosciences (PacBio) long-read sequences with 110× and 104× coverages were assembled and polished to generate 244 and 202 contigs spanning 808.8 Mbp for 'LA1670' and 804.5 Mbp for 'LA1673', respectively. After chromosome-level scaffolding with reference guiding, 14 scaffold sequences corresponding to 12 tomato chromosomes and 2 unassigned sequences were constructed. High-quality genome assemblies were confirmed using the Benchmarking Universal Single-Copy Orthologs and long terminal repeat assembly index. The protein-coding sequences were then predicted, and their transcriptomes were confirmed. The de novo assembled genomes of S. pimpinellifolium and S. lycopersicum var. cerasiforme were predicted to have 71,945 and 75,230 protein-coding genes, including 29,629 and 29,185 non-redundant genes, respectively, as supported by the transcriptome analysis results. The chromosome-level genome assemblies coupled with transcriptome data sets of the two accessions would be valuable for gaining insights into tomato domestication and understanding genome-scale breeding.

Multiple gene substitution by Target-AID base-editing technology in tomato.

The use of Target activation-induced cytidine deaminase (Target-AID) base-editing technology with the CRISPR-Cas 9 system fused with activation-induced cytidine deaminase (AID) resulted in the substitution of a cytidine with a thymine. In previous experiments focusing on a single target gene, this system has been reported to work in several plant species, including tomato (Solanum lycopersicum L.). In this research, we used Target-AID technology to target multiple genes related to carotenoid accumulation in tomato. We selected 3 genes, SlDDB1, SlDET1 and SlCYC-B, for their roles in carotenoid accumulation. Among 12 edited T0 lines, we obtained 10 independent T0 lines carrying nucleotide substitutions in the three targeted genes, with several allelic versions for each targeted gene. The two edited lines showed significant differences in carotenoid accumulation. These results demonstrate that Target-AID technology is a highly efficient tool for targeting multiple genes with several allelic versions.

Fruit setting rewires central metabolism via gibberellin cascades.

Fruit set is the process whereby ovaries develop into fruits after pollination and fertilization. The process is induced by the phytohormone gibberellin (GA) in tomatoes, as determined by the constitutive GA response mutant procera However, the role of GA on the metabolic behavior in fruit-setting ovaries remains largely unknown. This study explored the biochemical mechanisms of fruit set using a network analysis of integrated transcriptome, proteome, metabolome, and enzyme activity data. Our results revealed that fruit set involves the activation of central carbon metabolism, with increased hexoses, hexose phosphates, and downstream metabolites, including intermediates and derivatives of glycolysis, the tricarboxylic acid cycle, and associated organic and amino acids. The network analysis also identified the transcriptional hub gene SlHB15A, that coordinated metabolic activation. Furthermore, a kinetic model of sucrose metabolism predicted that the sucrose cycle had high activity levels in unpollinated ovaries, whereas it was shut down when sugars rapidly accumulated in vacuoles in fruit-setting ovaries, in a time-dependent manner via tonoplastic sugar carriers. Moreover, fruit set at least partly required the activity of fructokinase, which may pull fructose out of the vacuole, and this could feed the downstream pathways. Collectively, our results indicate that GA cascades enhance sink capacities, by up-regulating central metabolic enzyme capacities at both transcriptional and posttranscriptional levels. This leads to increased sucrose uptake and carbon fluxes for the production of the constituents of biomass and energy that are essential for rapid ovary growth during the initiation of fruit set.

Genome editing in PDS genes of tomatoes by non-selection method and of Nicotiana benthamiana by one single guide RNA to edit two orthologs.

The CRISPR/Cas9 system is widely used for targeted mutagenesis in many organisms including plants. For application of this system, tissue culture methods need to be established. In this study, detailed methods for introduction of mutations in tomato and Nicotiana benthamiana plants using the CRISPR/Cas9 system are described. The methods include tissue culture protocols for tomato and N. benthamiana. We also demonstrate the methodology to generate Cas9-free genome edited tomato plants and use of one single guide RNA (sgRNA) to edit two orthologs in N. benthamiana. The examples of editing the PHYTOENE DESATURASE (PDS) genes in these plants are also provided. The Cas9-free tomato line was obtained when tomato plants were cultured on a non-selective medium after transformation with the CRISPR/Cas9 system. Two orthologs of PDS in N. benthamiana were mutated using a sgRNA, because these orthologs contain the same nucleotide sequences with PAM motif. These mutations were inherited to the next generation. The mutations in the PDS genes resulted in an albino phenotype in tomato and N. benthamiana plants. These results demonstrate that the non-selective method is one of the ways to obtain Cas9-free genome editing in tomato plants and that the two orthologs can be edited by one sgRNA in N. benthamiana.

Comparative genomics of muskmelon reveals a potential role for retrotransposons in the modification of gene expression.

Melon exhibits substantial natural variation especially in fruit ripening physiology, including both climacteric (ethylene-producing) and non-climacteric types. However, genomic mechanisms underlying such variation are not yet fully understood. Here, we report an Oxford Nanopore-based high-grade genome reference in the semi-climacteric cultivar Harukei-3 (378 Mb + 33,829 protein-coding genes), with an update of tissue-wide RNA-seq atlas in the Melonet-DB database. Comparison between Harukei-3 and DHL92, the first published melon genome, enabled identification of 24,758 one-to-one orthologue gene pairs, whereas others were candidates of copy number variation or presence/absence polymorphisms (PAPs). Further comparison based on 10 melon genome assemblies identified genome-wide PAPs of 415 retrotransposon Gag-like sequences. Of these, 160 showed fruit ripening-inducible expression, with 59.4% of the neighboring genes showing similar expression patterns (r > 0.8). Our results suggest that retrotransposons contributed to the modification of gene expression during diversification of melon genomes, and may affect fruit ripening-inducible gene expression.

The inhibition of SlIAA9 mimics an increase in endogenous auxin and mediates changes in auxin and gibberellin signalling during parthenocarpic fruit development in tomato.

Parthenocarpic fruit formation can be achieved through the inhibition of SlIAA9, a negative regulator of auxin signalling in tomato plant. During early fruit development under SlIAA9 inhibition, cell division and cell expansion were observed. Bioactive gibberellin (GA) accumulated, but indole-3-acetic acid (IAA) and trans-zeatin did not accumulate substantially. Furthermore, under SlIAA9 inhibition, auxin-responsive genes such as SlIAA2, -3, and -14 were upregulated, and SlARF7 was downregulated. These results indicate that SlIAA9 inhibition mimics an increase in auxin. The auxin biosynthesis genes SlTAR1, ToFZY, and ToFZY5 were stimulated by an increase in auxin and by auxin mimicking under SlIAA9 inhibition. However, SlTAR2 and ToFZY2 were upregulated only by pollination followed by high IAA accumulation. These results suggest that SlTAR2 and ToFZY2 play an important role in IAA synthesis in growing ovaries. GA synthesis was also activated by SlIAA9 inhibition through both the early-13-hydroxylation (for GA1 synthesis) and non-13-hydroxylation (GA4) pathways, indicating that fruit set caused by SlIAA9 inhibition was partially mediated by the GA pathway. SlIAA9 inhibition induced the expression of GA inactivation genes as well as GA biosynthesis genes except SlCPS during early parthenocarpic fruit development in tomato. This result suggests that inactivation genes play a role in fine-tuning the regulation of bioactive GA accumulation.