# SIMYB75, an MYB-type transcription factor, promotes anthocyanin accumulation and enhances volatile aroma production in tomato fruits

Wei Jian, Haohao Cao, Shu Yuan, Yudong Liu, Juanfang Lu, Wang Lu, Ning Li, Jianhui Wang, Jian Zou, Ning Tang, et al.

► To cite this version:

Wei Jian, Haohao Cao, Shu Yuan, Yudong Liu, Juanfang Lu, et al.. SIMYB75, an MYB-type transcription factor, promotes anthocyanin accumulation and enhances volatile aroma production in tomato fruits. Horticulture research, 2019, 6 (1), pp.6-22. 10.1038/s41438-018-0098-y. hal-02626965

HAL Id: hal-02626965

https://hal.inrae.fr/hal-02626965

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# SIMYB75, an MYB-type transcription factor, promotes anthocyanin accumulation and enhances volatile aroma production in tomato fruits

Wei Jian<fnref n="1" />, Haohao Cao<fnref n="1" />, Shu Yuan<fnref n="1" />,<fnref n="2" />, Yudong Liu<fnref n="1" />, Juanfang Lu³, Wang Lu³, Ning Li<fnref n="1" />, Jianhui Wang⁴, Jian Zou<fnref n="1" />, Ning Tang⁵, Chan Xu<fnref n="1" />, Yulin Cheng<fnref n="1" />, Yanqiang Gao<fnref n="1" />, Wanpeng Xi³, Mondher Bouzayen<fnref n="1" />,⁶ and Zhengguo Li<fnref n="1" />

## Abstract

Genetic manipulation of genes to upregulate specific branches of metabolic pathways is a method that is commonly used to improve fruit quality. However, the use of a single gene to impact several metabolic pathways is difficult. Here, we show that overexpression of the single gene *SIMYB75* (*SIMYB75*-OE) is effective at improving multiple fruit quality traits. In these engineered fruits, the anthocyanin content reached 1.86 mg g⁻<fnref n="1" /> fresh weight at the red-ripe stage, and these *SIMYB75*-OE tomatoes displayed a series of physiological changes, including delayed ripening and increased ethylene production. In addition to anthocyanin, the total contents of phenolics, flavonoids and soluble solids in *SIMYB75*-OE fruits were enhanced by 2.6, 4, and 1.2 times, respectively, compared to those of wild-type (WT) fruits. Interestingly, a number of aroma volatiles, such as aldehyde, phenylpropanoid-derived and terpene volatiles, were significantly increased in *SIMYB75*-OE fruits, with some terpene volatiles showing more than 10 times higher levels than those in WT fruits. Consistent with the metabolic assessment, transcriptomic profiling indicated that the genes involved in the ethylene signaling, phenylpropanoid and isoprenoid pathways were greatly upregulated in *SIMYB75*-OE fruits. Yeast one-hybrid and transactivation assays revealed that SIMYB75 is able to directly bind to the MYBPLANT and MYBPZM cis-regulatory elements and to activate the promoters of the *LOXC*, *AADC2* and *TPS* genes. The identification of SIMYB75 as a key regulator of fruit quality attributes through the transcriptional regulation of downstream genes involved in several metabolic pathways opens new avenues towards engineering fruits with a higher sensory and nutritional quality.

## Introduction

As an indispensable daily supplier of nutrition to humans, fruit offers abundant nutrients, including vitamins, minerals, antioxidants and fiber<fnref n="1" />. Therefore, improving the nutritional and sensory qualities of fruit is becoming a primary goal to meet consumer requirements

for healthier food. Nevertheless, breeding efforts to date have primarily concentrated on yield and disease resistance traits, whereas sensory and nutritional aspects such as the pigment, aroma, taste and health-promoting compounds have not received sufficient attention. In recent years, tremendous advances in functional genomics and biotechnology tools have presented new prospects for producing high-quality fruits that can better meet consumer expectations.

Anthocyanins are water-soluble pigments that are widely distributed in many plants, where they confer their typical appealing color to various tissues and organs in

Correspondence: Zhengguo Li (zhengguoli@cqu.edu.cn)

<figure><img src="image_2.png" /></figure>

<footnote n="1">¹Key Laboratory of Plant Hormones and Development Regulation of Chongqing, School of Life Sciences, Chongqing University, 401331 Chongqing, China</footnote>

<footnote n="2">²College of Resources, Sichuan Agricultural University, 611130 Chengdu, China
Full list of author information is available at the end of the article.</footnote>
a number of species. With regards to their impact on human health, numerous studies have indicated that anthocyanins have high antioxidant properties, with the scavenging of active oxygen species². In addition, anthocyanins have been reported to play a role in protecting against age-related degenerative diseases, especially by inhibiting tumor cell growth³. Anthocyanins are synthesized via the flavonoid pathway, which is a branch of the phenylalanine pathway⁴. The accumulation of anthocyanins in plants is regulated by a series of environmental conditions and developmental signals⁵,⁶. Two types of genes have been suggested to participate in the anthocyanin pathway, (i) those encoding key enzymes for anthocyanin biosynthesis⁴ and (ii) those coding for the transcription factors (TFs) that regulate the spatiotemporal expression of functional genes⁷. Three types of TFs, R2R3-MYB, bHLH and WD40, are known to control the anthocyanin biosynthetic pathway, and most likely, these TFs work together through the formation of a transcription complex named MBW⁸,⁹. However, several studies have indicated that MYB TFs could also function independently of the MBW complex to regulate the biosynthesis of anthocyanins, thus assigning a crucial role to MYB factors *per se* in controlling anthocyanin biosynthesis⁷,¹⁰,¹¹.

Tomatoes are one of the most widely consumed fruit crops in the world, which motivates their frequent use as a model species to study nutritional metabolism. Moreover, the tomato is an excellent candidate for the transgenic enhancement of its flavonoid content because of its low basal content of these compounds³,¹². Although tomato plants with both *Aft* and *atv* alleles could accumulate small amounts of anthocyanin, the mechanism and key genes controlling this variation have not yet been deciphered¹³,¹⁴. The silencing of the *SIDET1* gene in tomatoes has resulted in altered light-mediated developmental processes along with an enhanced flavonoid content¹⁵. Similarly, overexpression of the chalcone isomerase gene in tomatoes could increase the flavonoid content up to 78-fold exclusively in the peel, but the overall level in the fruits remained low¹². In contrast to the functional genes, manipulating the expression of TFs proved more effective at modulating the secondary metabolism. Lc and C1 are two TFs in maize that control anthocyanin biosynthesis; however, overexpression of these two TFs in tomatoes greatly increased the flavonoid content, but not the anthocyanin content¹⁶. Overexpression of *AtMYB12* under the control of the fruit-specific *E8* promoter resulted in orange tomatoes that contained significantly high contents of flavonoid and caffeoyl quinic acid¹⁷. However, the most significant increase in the anthocyanin content of tomato fruits was first obtained by expressing *Del* and *Ros1*, two snapdragon genes encoding bHLH and MYB TFs, respectively, which were under the control of

the *E8* promoter, resulting in fully purple tomato fruits³. More recently, crossing the *AtMYB12*-overexpression line with the *Del/Ros1* line was reported to enhance anthocyanin accumulation¹⁸. Interestingly, this study indicated that *AtMYB12* not only promoted flavonoid biosynthesis, but it also improved the carbon supply from the primary metabolism, energy and reducing power, which resulted in a larger aromatic amino acid supply for secondary metabolism¹⁸. It is important to mention that the purple tomato fruits obtained in these studies always resulted from the cooperative action of two or more TFs. This result causes difficulty for traditional breeders, since it is challenging to find a wild species that could highly express these two or three key regulatory genes at the same time.

Recently, the *MYB* TF gene *Cs6g17570* was identified to play a critical role in regulating anthocyanin biosynthesis in blood oranges¹⁹. Interestingly, the upregulation of its homologous gene *AtMYB113-like* resulted in the increased accumulation of anthocyanin in *Arabidopsis*, whereas its downregulation led to a reduced amount of this pigment. Moreover, this process was reported to be dependent on the MBW complex²⁰. Through a homology alignment, we identified a tomato gene, named *SlMYB75* (SGN, https://solgenomics.net/search/locus) or *SlAN2*, which showed the highest sequence similarity with *Cs6g17570* (Supplementary Fig. S1). Overexpression of *SlAN2* was previously reported to improve tomato plant resistance to high temperature, cold and oxidative stresses²¹,²². A recent study showed that overexpression of *SlANT1* or *SlAN2* induced accumulation of anthocyanin that was unevenly distributed in tomato fruit, but only *SlAN2* was able to act as a positive regulator of anthocyanin biosynthesis in vegetative tomato tissues under stress conditions²³. At the same time, *SlAN2-overexpression* fruits displayed an orange color, fast softening and higher ethylene content²⁴. In addition, several other MYB TFs are also involved in some fruit quality attributes, such as primary metabolism, secondary metabolism and organic acid metabolism¹⁸,²⁵, and a consumer sensory analysis showed that there is a consumer preference for purple tomato fruits, for their better perceived flavor²⁶. Thus, the potential role of *SlMYB75* in tomato fruit quality remains to be investigated.

In the present study, we attempt to describe the physiological and metabolic changes in purple tomatoes produced by overexpressing a single *SlMYB75* TF (*SlMYB75*-OE). The data indicated that *SlMYB75* can effectively induce the accumulation of anthocyanin in various tissues, and its ectopic expression leads to increased ethylene production and enhanced phenolic, flavonoid and volatile aroma contents. Moreover, transcriptomic profiling showed that many genes involved in the ethylene response, phenylpropanoid and isoprenoid pathways were greatly affected in *SlMYB75*-OE fruits.
Using yeast one-hybrid and dual-luciferase assays, we found that the conserved MYBPLANT (AAACCAACCC) and MYBPZM (ACCTACCC) elements were the core binding sites of SIMYB75, and we revealed the ability of SIMYB75 to activate the promoters of the *LOXC*, *AADC2* and *TPS* genes. This study demonstrates that a single SIMYB75 TF is able to increase the anthocyanin content to 1.86 mg g⁻¹ fresh weight at the red ripening stage and promote volatile aroma accumulation partly through the transcriptional regulation of downstream genes involved in the corresponding metabolic pathways.

# Materials and methods

## Plant materials, growth conditions and chemicals

*Solanum lycopersicum* cv. Micro-Tom was selected as the wild type (WT) in this study. All the tomato seedlings were grown in an intelligent greenhouse under standard conditions (16/8 h and 25 °C/18 °C day/night cycle, 80% humidity and 250 µmol m⁻² s⁻¹ light intensity). For the hormone and stress treatments, the leaves of 1-month-old tomato plants treated with different reagents were harvested. For each tissue/organ type, samples were collected from at least six healthy plants. All the samples were mixed and frozen under liquid nitrogen immediately. The sample experiments were conducted three independent times.

For the chemicals, Folin-Ciocalteu phenol reagent was purchased from Sigma (St. Louis, MO, USA). HPLC-grade methanol was purchased from Merck KgaA (Darmstadt, Germany). The other reagents were all of analytical grade and were purchased from Sangon Biotechnology Co., Ltd. (Shanghai, China).

## Vector construction and plant transformation

The full length of the *SIMYB75* coding sequence was amplified from tomato cDNA and then cloned into the K303 expression vector, which contain two 35s promoters from Gateway technology. The final vector was transferred into GV3101, and *Agrobacterium*-mediated transformation was performed as described previously²⁷,²⁸. The positive transgenic plants were identified by PCR, and homozygous plants from T2 or later generations were used for the experiments. All the primers used in this study are listed in Supplementary Table S2.

## RNA-Seq analysis

Total RNA was extracted from the fruits at the mature green (MG) and breaker (BR + 0) stages using an RNeasy Plant Mini Kit (Tiangen, China) according to the manufacturer's protocol. Samples from both WT and *SIMYB75*-OE fruits were collected with two biological replicates. The concentration and quality of the RNA were assayed using a NanoDrop Lite spectrophotometer

(Thermo Scientific). The cDNA libraries were constructed and then sequenced using a BGISEQ-500 System (BGI Inc.). Clean data were obtained after they were filtered and then mapped to the reference genome of *S. lycopersicum* in the Tomato SGN database (http://solgenomics.net/) with *Bowtie2*. The homogenized data were used to calculate the gene expression levels with *RSEM*. The criteria for defining differentially expressed genes (DEGs) were fold change ≥2.00 and *P* value ≤ 0.05. The raw transcriptome reads of this study have been deposited in the NCBI Short Read Archive under accession number SRP158557.

## Real time quantitative PCR

One microgram of total RNA (RNeasy Plant Mini Kit, Tiangen) was used to synthesize first strand cDNA with PrimeScript™ RT reagent Kit using gDNA Eraser (Perfect Real Time) (TAKARA, Japan). Quantitative real-time PCR was performed with a Bio-Rad CFX system (Bio-Rad, USA) using SYBR® Premix Ex Taq™ (Tli RNaseH Plus) (TAKARA, Japan). Each sample was collected as three independent biological replicates, and the relative fold differences were calculated using a comparative Ct method. *SliUBI* was used as the internal reference over the entire experiment.

## Measurement of metabolites

For the total anthocyanin extraction, 0.1 g of lyophilized powder samples was incubated in 1 mL of pH 1.0 buffer solution (50 mM KCL, 150 mM HCL) and 1 mL of pH 4.5 buffer solution (400 mM CH₃COONa, 240 mM HCL) for 24 h in the dark at 4 °C with gentle shaking, respectively. The samples were centrifuged at 14,000×g for 20 min. The total anthocyanin content in the supernatant was measured spectrophotometrically and expressed in mg of petunidin-3-(p-coumaroyl rutinoside)-5-glucoside (extinction coefficient 17000, molecular weight 934) per gram fresh weight³,²³.

For the ethylene measurement, fruits at different development stages were harvested and placed in an open 50 mL jar for 2 h to minimize the wound ethylene caused by picking. The jars were sealed and incubated at room temperature for 1 h, and 1 mL of headspace gas was collected and then injected into a VARIAN CP3800 gas chromatograph equipped with a flame ionization detector (USA). The samples were compared with reagent-grade ethylene standards of known concentrations and normalized to the fruit weight²⁹. At least 10 individual fruits were measured for each sample.

The total soluble solids and titratable acids were determined using a digital refractometer (Atago PR-101R, Tokyo, Japan) and titration method, respectively³⁰. Each replicate contained 20 fresh tomato fruits and all the determinations were performed in triplicate.
For the total phenolic and flavonoid extractions, 1 g of lyophilized powder samples were incubated with methanol (80%, 24 mL) at 25 °C for 12 h with shaking and then centrifuged at 3000×g for 10 min at 4 °C. The residue was extracted twice more using the same procedure, and then the supernatants were collected and finally diluted to 50 mL with methanol to detect the total phenolics and flavonoids. The total phenolic content was determined using the Folin-Ciocalteu method. The total flavonoid content was measured as described previously³¹.

To measure the primary component of phenolics and flavonoids, 0.5 g of lyophilized powder was extracted with methanol (80%, 12 mL) and dimethyl sulfoxide (1:1, v/v) using the same procedures as the procedure used with the total phenolics and flavonoids. After the samples were filtered through a 0.22 µm syringe filter, the phenolics and flavonoids were measured using the HPLC method described previously³¹.

The concentration of aroma volatiles was determined as described previously, with some modifications³⁰. Three grams of lyophilized powder were homogenized with 5 mL of saturated sodium chloride solution, and then 5 µL of ethyl nonanoate was added as an internal reference. The solution was incubated at 40 °C for 30 min, and a solid-phase microextraction (SPME) needle with a 1 cm-long fiber coated with 65 µm divinylbenzene/carboxen/ polydimethylsiloxane (DVB/CAR/PDMS) fibers (Supelco Co., Bellefonte PA, USA) was used to extract the volatiles. A GCMS-QP2010 gas chromatograph-mass spectrometer system (Shimadzu Corporation, Kyoto, Japan) with an Rtx-5MS (Restek)-fused silica capillary column (5% diphenyl, 95% dimethyl polysiloxane) (0.32 mm, 30 m, 0.5 lm, J&W Scientific, Folsom CA, USA) was used for the compound confirmation. GC-MS Postrun Analysis software (SHIMADZU, GC-MS-QP2010, Japan) was used to evaluate the chromatograms and mass spectra. The compounds were identified by comparing their mass spectra with the data system library (NIST08). The concentrations of the volatile components were expressed as ng g⁻¹ FW h⁻¹.

## Yeast one-hybrid assay

Yeast one-hybrid (Y1H) assays were performed using a Matchmaker Gold Yeast One Hybrid System (Clontech). To construct the prey and bait vectors, the full length of the *SIMYB75* open reading frame (ORF) sequence was cloned and inserted into the pGADT7 vector, and the eight conserved *cis*-elements were cloned into the pAbAi vector. The bait plasmids were transformed into the Y1H Gold strain according to the manufacturer's instructions. Aureobasidin A (AbA) was used to screen the minimal inhibitory concentration for the bait strains. The prey plasmid was transformed into a bait yeast strain to determine the DNA-protein interaction by

screening them on SD medium with AbA and without leucine.

## Dual-luciferase transient expression assay

For the dual-luciferase assay, the full-length ORF of *SIMYB75* was cloned and inserted into a pGreenII 62-SK vector (effector), and the promoter sequence of the specific different expression gene (DEG) was cloned and inserted into pGreenII 0800 LUC vector (reporter). After their transformation into GV3101, the effector and reporter strains were cultured and then resuspended with infiltration buffer (10 mM MES, 10 mM MgCl₂, 200 mM acetosyringone, pH 5.6) to an OD 600 of 1.0–1.5. The mixtures of effector and reporter cultures were infiltrated into *Nicotiana benthamiana* leaves with needleless syringes. Firefly luciferase and *Renilla* luciferase were assayed at 3 days after infiltration using Dual Luciferase Reporter Assay System reagents (Promega). The binding activity of *SIMYB75* to the promoter of each specific DEG was calculated by finding the LUC to REN ratio. At least six biological replicates were conducted for each combination³².

## Statistical analysis

All the experiments were repeated at least three times, and the results were presented with the standard deviations. Student's *t*-tests were used to analyze the data, and a difference was considered to be statistically significant when *P* < 0.05.

## Results

### Overexpression of a single *SIMYB75* induces anthocyanin accumulation

A qRT-PCR investigation of the *SIMYB75* expression pattern indicated that transcripts corresponding to this gene are detected in all tissues, but they display higher expression levels in vegetative tissues, especially in leaf organs (Fig. 1a and b). However, unlike the *Cs6g17570* gene¹⁹, the expression level of *SIMYB75* in tomato fruit organs is very low (Fig. 1a, b). To address the functional significance of *SIMYB75*, a tomato “Micro-Tom” cultivar was transformed with a sense construct of *SIMYB75*, leading to the generation of three independent transgenic homozygous lines. The expression of *SIMYB75* in *SIMYB75*-OE plants was assessed by qRT-PCR, indicating a dramatic upregulation in line #11 and line #19, while the #21 line only showed a slight upregulation (Fig. 1c), and then lines #11 and 19 were selected to perform further experiments. A phenotypic evaluation revealed that the *SIMYB75*-OE lines accumulated abundant amounts of anthocyanins in both vegetative and reproductive organs, especially in the stamen, in which the color turned totally purple (Fig. 1d). Furthermore, the
<figure><img src="image_3.png" /></figure>

SIMYB75-OE lines displayed significantly smaller seeds than the WT (Fig. 1d). To further investigate the expression characteristics of SIMYB75, we checked the expression level of this TF in WT seedlings treated with different hormones or ones that were subjected to stress conditions. The data showed that SIMYB75 is responsive to all the treatments applied here, consistent with the presence of the corresponding cis-elements identified in the promoter of SIMYB75 (Supplementary Fig. S2). These results indicated that overexpressing a single SIMYB75 TF can lead to abundant anthocyanin accumulation in both vegetative and reproductive organs, and this TF could be induced by various hormones or stress conditions.

## Physiological features of the SIMYB75-OE tomato plants

One of the striking features of the SIMYB75-OE lines is the purple color exhibited by several organs, including the seeds, stamens and fruits. An assessment of the anthocyanin accumulation revealed high levels in the SIMYB75-OE fruits, and these compounds were not detected in the WT (Fig. 2a, b). Despite their deep purple color, it was

easy to identify the different ripening stages of the transgenic fruits by checking the color of the low-anthocyanin region (Fig. 2a). Moreover, SIMYB75-OE plants exhibited 3–5-day delays in the occurrence of fruit ripening onset compared to the WT (Fig. 2c, d). Strikingly, ripening-associated ethylene production was significantly higher in SIMYB75-OE tomatoes at all ripening stages (Fig. 2e), which contrasted with the delayed ripening. These data suggested that overexpression of a single SIMYB75 TF can result in a series of physiological changes.

## Transcriptome profiling of the SIMYB75-OE tomato fruits

To gain insight into the extent of transcriptomic reprogramming induced by the ectopic expression of SIMYB75, we performed genome-wide transcriptomic profiling in WT and SIMYB75-OE fruits (line #11) at the MG and BR + 0 stages. Compared to the WT, a total of 1576 and 2843 genes were differently expressed (DEGs) in SIMYB75-OE tomatoes at the MG and BR + 0 stages, respectively (Files S1 and S2). Based on the Gene Ontology (GO) assignment, the common DEGs of the MG and
<figure><img src="image_4.png" /></figure>

BR + 0 stages (675 common DEGs) were assigned to three categories, namely, biological process, cellular component and molecular function. “Biological process” was primarily focused on metabolic and cellular processes. Alterations in the molecular function analysis were primarily involved in the binding and catalytic activities (Supplementary Fig. S3 and File S3). Furthermore, 675 common DEGs were enriched in 73 KEGG pathways. The top five enriched pathways were the biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, protein processing in the endoplasmic reticulum, plant hormone signal transduction and flavonoid biosynthesis (Supplementary Fig. S4 and File S4). These results indicated that SIMYB75 overexpression impacted multiple processes, including transcription, stress responses, secondary metabolism and phytohormone signaling

pathways. Based on the phenotype properties and potential application values of the SIMYB75-OE tomato, we focused on genes involved in metabolic processes, ripening-related processes and TFs (Table 1). To validate the accuracy of the generated transcriptome data, a total of 16 genes related to the above-cited processes were selected for qRT-PCR analysis. The outcome of this targeted expression analysis was highly consistent with the transcriptome data for all the tested genes (Fig. 3).

## Change in the ripening-related metabolic processes in SIMYB75-OE tomatoes

To further investigate the metabolic changes in the SIMYB75-OE tomatoes, we measured the metabolite contents relevant to the DEGs at the BR + 5, BR + 7 and BR + 10 ripening stages. Interestingly, the total soluble
<table><caption>Table 1 List of different expression genes (DEGs) between WT and #11 tomato fruits at the MG and BR + 0 stages. Genes marked with asterisks were validated by qPCR</caption><thead><tr><th>ITAG 2.40 Tomato</th><th>Log2-fold (#11/WT; MG)</th><th>Log2-fold (#11/WT; BR + 0)</th><th>Functional annotation</th></tr></thead><tbody><tr><td>Phenylalanine pathway</td><td></td><td></td><td></td></tr><tr><td>Solyc10g086180.2</td><td>1.66</td><td>1.07</td><td>Phenylalanine ammonia-lyase</td></tr><tr><td>Solyc00g282510.2*</td><td>1.68</td><td>1.03</td><td>Phenylalanine ammonia-lyase-like</td></tr><tr><td>Solyc05g052240.3</td><td>1.34</td><td>1.71</td><td>Probable chalcone-flavonone isomerase 3</td></tr><tr><td>Solyc08g080040.3*</td><td>1.51</td><td>1.68</td><td>Leucoanthocyanidin dioxygenase-like</td></tr><tr><td>Solyc09g091510.3*</td><td>1.87</td><td>1.60</td><td>Chalcone synthase 1</td></tr><tr><td>Solyc12g088170.2</td><td>8.96</td><td>10.01</td><td>Anthocyanin acyltransferase</td></tr><tr><td>Solyc10g083440.1*</td><td>1.85</td><td>1.46</td><td>Anthocyanidin 3-O-glucosyltransferase-like</td></tr><tr><td>Flavor related</td><td></td><td></td><td></td></tr><tr><td>Solyc01g006540.3*</td><td>2.86</td><td>1.19</td><td>LOXC</td></tr><tr><td>Solyc01g099190.3*</td><td>2.49</td><td>1.46</td><td>LOXB</td></tr><tr><td>Solyc01g006560.3</td><td>2.25</td><td>1.38</td><td>LOXF</td></tr><tr><td>Solyc08g006740.3</td><td>2.00</td><td>1.71</td><td>AADC2</td></tr><tr><td>Solyc01g108560.4*</td><td>1.03</td><td>2.16</td><td>CXE1</td></tr><tr><td>Solyc06g060180.2*</td><td>4.40</td><td>1.90</td><td>TPS</td></tr><tr><td>Solyc10g011920.2</td><td>1.10</td><td>1.66</td><td>Aromatic amino acid lyase</td></tr><tr><td>Ripening related pathway</td><td></td><td></td><td></td></tr><tr><td>Solyc01g095080.3*</td><td>2.20</td><td>1.94</td><td>ACS2</td></tr><tr><td>Solyc05g050010.3*</td><td>1.80</td><td>1.98</td><td>ACS4</td></tr><tr><td>Solyc03g111720.3*</td><td>1.91</td><td>3.69</td><td>E4</td></tr><tr><td>Solyc09g089580.3*</td><td>1.74</td><td>1.50</td><td>E8</td></tr><tr><td>Solyc09g075440.3</td><td>1.23</td><td>1.76</td><td>ETR3(NR)</td></tr><tr><td>Solyc05g012020.3*</td><td>3.09</td><td>1.68</td><td>RIN</td></tr><tr><td>Solyc03g044300.3</td><td>2.39</td><td>1.40</td><td>AP2a</td></tr><tr><td>Solyc03g031860.3</td><td>2.60</td><td>1.71</td><td>PSY1</td></tr><tr><td>Solyc10g080210.2</td><td>2.54</td><td>2.33</td><td>PG2a</td></tr><tr><td>Other genes</td><td></td><td></td><td></td></tr><tr><td>Solyc09g065100.2*</td><td>0.62</td><td>1.34</td><td>bHLH150</td></tr><tr><td>Solyc12g007070.2</td><td>2.62</td><td>1.28</td><td>Heat stress transcription factor C-1</td></tr><tr><td>Solyc11g017470.2*</td><td>2.39</td><td>1.02</td><td>NAC domain-containing protein 2-like</td></tr><tr><td>Solyc10g079050.2</td><td>1.46</td><td>1.71</td><td>Transcription factor bHLH130-like</td></tr><tr><td>Solyc10g009550.3</td><td>-2.50</td><td>-1.06</td><td>Probable WRKY transcription factor 30</td></tr><tr><td>Solyc09g014990.3</td><td>-1.93</td><td>-1.18</td><td>Probable WRKY transcription factor 26</td></tr><tr><td>Solyc06g051260.3*</td><td>-1.46</td><td>-1.25</td><td>Transcription factor bHLH51</td></tr><tr><td>Solyc03g115850.3</td><td>-0.99</td><td>-1.63</td><td>NAC domain-containing protein 100</td></tr><tr><td>Solyc10g080030.2</td><td>-1.89</td><td>-1.69</td><td>MADS-box transcription factor 23-like</td></tr><tr><td>Solyc12g014140.2</td><td>-1.93</td><td>-1.81</td><td>TCP transcription factor 3</td></tr><tr><td>Solyc11g028020.2</td><td>-1.35</td><td>-2.09</td><td>TAGL11 transcription factor</td></tr><tr><td>Solyc10g005010.3</td><td>-2.31</td><td>-3.38</td><td>NAC domain-containing protein 43</td></tr></tbody></table>
<table><caption>**Table 1 continued**</caption><thead><tr><th>ITAG 2.40 Tomato</th><th>Log2-fold (#11/WT; MG)</th><th>Log2-fold (#11/WT; BR + 0)</th><th>Functional annotation</th></tr></thead><tbody><tr><td>Solyc08g080490.3</td><td>-1.52</td><td>-1.53</td><td>2S sulfur-rich seed storage protein 2-like</td></tr><tr><td>Solyc09g025210.3</td><td>-1.68</td><td>-2.32</td><td>12S seed storage protein CRA1-like</td></tr><tr><td>Solyc09g072560.3</td><td>-1.68</td><td>-2.55</td><td>11S globulin seed storage protein 2-like</td></tr></tbody></table>

<figure><img src="image_5.png" /></figure>

solids (TSS), which primarily reflect the sugar content, was significantly higher in the *SIMYB75*-OE tomatoes than in the WT (Fig. 4a) while the titratable acids showed no significant difference between WT and *SIMYB75*-OE tomatoes (Fig. 4b). Strikingly, the total contents of phenolics and flavonoids were dramatically increased in

*SIMYB75*-OE tomatoes (Fig. 4c, d). A further characterization of the changes in phenolics and flavonoids contents revealed that the *SIMYB75*-OE fruits exhibited higher levels of chlorogenic acid, neochlorogenic acid, ferulic acid, gallic acid, rutin and quercetin than those in the WT (Fig. 5). These data indicated that in addition to
<figure><img src="image_6.png" /><figcaption>**Fig. 4 Determination of the total soluble solids, titratable acids, total phenolics and total flavonoids in WT and *SIMYB75*-OE fruits at BR + 7.** Total soluble solids (a) and titratable acids (b) were measured using a digital refractometer and a titration method, respectively. The data are the means ± SD of at least 10 individual fruits for each line. The total phenolics (c) and total flavonoids (d) were calculated according to the gallic acid and rutin equivalents, respectively. All data are the means ± SD of three biological replicates, *P* < 0.05; and **P** < 0.01 (Student's t-test)</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 5 Phenolic and flavonoid accumulation in WT and *SIMYB75*-OE fruits at BR + 5, BR + 7, BR + 10.** Chlorogenic acid (a), neochlorogenic acid (b), ferulic acid (c), gallic acid (d), rutin (e), and quercetin (f) were quantified according to their corresponding standards using an HPLC method. The data are the means ± SD of three biological replicates analyzed by Student's t-test; *P* < 0.05; and **P** < 0.01</figcaption></figure>
<table><caption>Table 2 Aroma volatiles emitted by WT and S/MyB75-OE fruits at the BR + 5, BR + 7 and BR + 10 stages. The data are expressed as the means ± standard deviation of triplicate samples</caption><thead><tr><th rowspan="2">Volatile compound</th><th rowspan="2">RT (min)</th><th colspan="3">BR+5 Content (ng g<sup>-1</sup> FW h<sup>-1</sup>)</th><th colspan="3">BR+7 Content (ng g<sup>-1</sup> FW h<sup>-1</sup>)</th><th colspan="3">BR+10 Content (ng g<sup>-1</sup> FW h<sup>-1</sup>)</th></tr><tr><th>WT</th><th>#11</th><th>#19</th><th>WT</th><th>#11</th><th>#19</th><th>WT</th><th>#11</th><th>#19</th></tr></thead><tbody><tr><td>Aldehydes</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>trans-2-Hexenal</td><td>4.90</td><td>259.65 ± 22.48</td><td>421.29 ± 31.71</td><td>331.96 ± 22.38</td><td>193.00 ± 32.87</td><td>368.42 ± 50.84</td><td>257.05 ± 36.45</td><td>219.73 ± 35.18</td><td>352.31 ± 25.68</td><td>303.40 ± 25.90</td></tr><tr><td>Hexanal</td><td>3.66</td><td>77.29 ± 11.74</td><td>116.96 ± 22.90</td><td>97.84 ± 7.12</td><td>44.35 ± 11.82</td><td>89.52 ± 21.54</td><td>77.92 ± 5.11</td><td>74.98 ± 8.37</td><td>196.91 ± 46.12</td><td>128.03 ± 19.14</td></tr><tr><td>3-Hexenal</td><td>3.62</td><td>33.29 ± 2.98</td><td>56.18 ± 3.15</td><td>46.63 ± 2.37</td><td>12.07 ± 0.96</td><td>42.74 ± 2.73</td><td>28.51 ± 2.64</td><td>25.33 ± 2.89</td><td>48.25 ± 2.19</td><td>39.32 ± 3.28</td></tr><tr><td>(Z)-2-Heptenal</td><td>8.01</td><td>4.83 ± 0.54</td><td>11.46 ± 1.66</td><td>7.80 ± 0.61</td><td>6.38 ± 1.00</td><td>11.21 ± 1.80</td><td>8.86 ± 0.35</td><td>6.63 ± 0.52</td><td>14.43 ± 2.11</td><td>10.62 ± 0.74</td></tr><tr><td>(E,E)-2,4-Heptadienal</td><td>9.06</td><td>3.85 ± 0.24</td><td>7.64 ± 0.73</td><td>6.08 ± 0.13</td><td>2.95 ± 0.54</td><td>6.43 ± 1.76</td><td>4.77 ± 0.19</td><td>6.42 ± 1.72</td><td>16.48 ± 1.28</td><td>9.01 ± 1.47</td></tr><tr><td>Heptanal</td><td>6.55</td><td>2.59 ± 0.39</td><td>5.55 ± 0.70</td><td>4.97 ± 0.25</td><td>1.90 ± 0.52</td><td>4.82 ± 1.07</td><td>3.42 ± 0.43</td><td>3.16 ± 0.61</td><td>5.61 ± 1.75</td><td>4.38 ± 1.32</td></tr><tr><td>Geranial</td><td>15.05</td><td>12.27 ± 2.42</td><td>31.05 ± 2.99</td><td>21.22 ± 2.96</td><td>47.15 ± 7.63</td><td>276.86 ± 44.54</td><td>204.11 ± 41.07</td><td>66.92 ± 11.91</td><td>252.34 ± 15.73</td><td>179.60 ± 31.54</td></tr><tr><td>Neral</td><td>14.52</td><td>ND</td><td>6.33 ± 2.33</td><td>4.27 ± 1.03</td><td>ND</td><td>5.45 ± 0.46</td><td>3.18 ± 0.71</td><td>5.70 ± 1.37</td><td>22.40 ± 2.62</td><td>15.37 ± 1.96</td></tr><tr><td>Benzaldehyde</td><td>7.90</td><td>3.49 ± 0.89</td><td>6.38 ± 1.63</td><td>4.91 ± 1.10</td><td>3.73 ± 0.79</td><td>8.70 ± 1.13</td><td>6.70 ± 0.68</td><td>4.23 ± 0.92</td><td>7.63 ± 1.42</td><td>6.46 ± 0.45</td></tr><tr><td>Alcohols</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>6-Methyl-5-hepten-2-ol</td><td>9.30</td><td>9.63 ± 2.62</td><td>11.04 ± 2.40</td><td>10.23 ± 0.43</td><td>19.38 ± 8.69</td><td>23.81 ± 2.63</td><td>21.54 ± 4.15</td><td>16.45 ± 3.52</td><td>27.56 ± 4.53</td><td>24.25 ± 2.73</td></tr><tr><td>2-Propyl-1-pentanol</td><td>10.26</td><td>20.01 ± 5.13</td><td>10.04 ± 0.54</td><td>14.64 ± 2.68</td><td>13.82 ± 2.93</td><td>5.60 ± 3.16</td><td>8.75 ± 1.57</td><td>3.36 ± 1.25</td><td>1.91 ± 0.37</td><td>2.22 ± 0.75</td></tr><tr><td>Eucalyptol</td><td>14.91</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>108.95 ± 23.39</td><td>70.36 ± 9.40</td><td>ND</td><td>85.87 ± 7.43</td><td>53.91 ± 5.61</td></tr><tr><td>1-Octanol</td><td>11.23</td><td>4.23 ± 0.82</td><td>9.57 ± 1.00</td><td>6.20 ± 0.34</td><td>4.71 ± 1.16</td><td>10.07 ± 2.78</td><td>7.21 ± 0.77</td><td>5.16 ± 0.49</td><td>9.60 ± 1.30</td><td>6.73 ± 0.72</td></tr><tr><td>Matsutake alcohol</td><td>8.98</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>21.52 ± 2.68</td><td>14.39 ± 2.82</td><td>ND</td><td>13.66 ± 1.27</td><td>7.89 ± 0.48</td></tr><tr><td>Ketones</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>β-ionone</td><td>18.81</td><td>1.66 ± 0.12</td><td>1.88 ± 0.43</td><td>1.79 ± 0.29</td><td>1.77 ± 0.39</td><td>2.86 ± 0.64</td><td>2.19 ± 0.91</td><td>1.89 ± 0.59</td><td>2.78 ± 0.52</td><td>2.32 ± 0.31</td></tr><tr><td>6-Methyl-5-h eptene-2-one</td><td>8.96</td><td>64.57 ± 14.54</td><td>71.89 ± 18.57</td><td>69.16 ± 15.53</td><td>50.77 ± 7.89</td><td>58.18 ± 13.71</td><td>56.33 ± 12.84</td><td>161.80 ± 12.71</td><td>170.58 ± 36.07</td><td>165.23 ± 13.39</td></tr><tr><td>Geranylacetone</td><td>18.27</td><td>8.46 ± 1.52</td><td>15.84 ± 2.14</td><td>13.73 ± 1.92</td><td>13.82 ± 1.15</td><td>39.76 ± 5.49</td><td>27.06 ± 1.15</td><td>24.37 ± 4.34</td><td>40.88 ± 8.84</td><td>31.89 ± 5.01</td></tr><tr><td>3-Nonanone</td><td>11.49</td><td>2.88 ± 0.73</td><td>13.47 ± 2.38</td><td>8.18 ± 0.97</td><td>91.51 ± 0.40</td><td>202.60 ± 37.01</td><td>151.44 ± 22.51</td><td>140.37 ± 13.72</td><td>279.56 ± 24.59</td><td>216.84 ± 24.57</td></tr><tr><td>Farnesyl acetone</td><td>24.82</td><td>3.27 ± 0.30</td><td>5.35 ± 0.28</td><td>4.42 ± 0.18</td><td>4.35 ± 0.76</td><td>8.42 ± 0.84</td><td>7.09 ± 0.92</td><td>4.32 ± 0.71</td><td>9.78 ± 2.36</td><td>7.47 ± 1.08</td></tr><tr><td>Esters</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Methyl salicylate</td><td>13.59</td><td>4.75 ± 1.33</td><td>13.67 ± 1.07</td><td>11.19 ± 0.61</td><td>4.01 ± 1.69</td><td>20.95 ± 1.66</td><td>12.06 ± 1.06</td><td>2.88 ± 0.45</td><td>8.96 ± 0.89</td><td>6.97 ± 0.59</td></tr><tr><td>Geranyl acetate</td><td>17.15</td><td>ND</td><td>ND</td><td>ND</td><td>6.20 ± 1.04</td><td>31.49 ± 8.72</td><td>21.33 ± 3.07</td><td>6.78 ± 1.73</td><td>33.68 ± 4.80</td><td>21.09 ± 2.35</td></tr><tr><td>Nerol acetate</td><td>16.84</td><td>ND</td><td>ND</td><td>ND</td><td>4.65 ± 0.58</td><td>17.08 ± 2.55</td><td>10.57 ± 0.86</td><td>5.39 ± 1.33</td><td>16.32 ± 2.07</td><td>10.83 ± 1.61</td></tr></tbody></table>
<table><caption>Table 2 continued</caption><thead><tr><th rowspan="2">Volatile compound</th><th rowspan="2">RT (min)</th><th colspan="3">BR+5 Content (ng g<sup>-1</sup> FW h<sup>-1</sup>)</th><th colspan="3">BR+7 Content (ng g<sup>-1</sup> FW h<sup>-1</sup>)</th><th colspan="3">BR+10 Content (ng g<sup>-1</sup> FW h<sup>-1</sup>)</th></tr><tr><th>WT</th><th>#11</th><th>#19</th><th>WT</th><th>#11</th><th>#19</th><th>WT</th><th>#11</th><th>#19</th></tr></thead><tbody><tr><td>α-terpineol acetate</td><td>16.66</td><td>ND</td><td>ND</td><td>ND</td><td>3.82 ± 0.43</td><td>21.86 ± 3.01</td><td>13.74 ± 2.75</td><td>4.95 ± 0.90</td><td>21.79 ± 5.01</td><td>15.64 ± 3.53</td></tr><tr><td>2-Ethylhexyl acrylate</td><td>14.52</td><td>87.61 ± 16.45</td><td>ND</td><td>ND</td><td>61.36 ± 20.21</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>Isobutyl acetate</td><td>4.13</td><td>0.34 ± 0.05</td><td>0.54 ± 0.08</td><td>0.50 ± 0.10</td><td>0.20 ± 0.08</td><td>0.78 ± 0.08</td><td>0.55 ± 0.07</td><td>0.18 ± 0.03</td><td>0.30 ± 0.08</td><td>0.29 ± 0.06</td></tr><tr><td>Terpenoids/Terpenes</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Linalool</td><td>11.88</td><td>28.06 ± 6.20</td><td>134.50 ± 16.79</td><td>99.69 ± 6.14</td><td>34.47 ± 11.09</td><td>313.13 ± 40.92</td><td>234.12 ± 13.87</td><td>59.80 ± 13.92</td><td>370.44 ± 35.57</td><td>260.58 ± 38.45</td></tr><tr><td>Terpineol</td><td>13.68</td><td>4.48 ± 0.96</td><td>13.61 ± 3.12</td><td>11.10 ± 1.87</td><td>6.67 ± 3.57</td><td>27.87 ± 3.05</td><td>17.77 ± 5.06</td><td>8.62 ± 2.00</td><td>19.03 ± 4.04</td><td>12.59 ± 2.01</td></tr><tr><td>()-4-Terpineol</td><td>13.46</td><td>1.02 ± 0.46</td><td>11.08 ± 0.99</td><td>7.22 ± 0.47</td><td>16.37 ± 0.38</td><td>92.76 ± 14.58</td><td>69.05 ± 19.86</td><td>18.34 ± 3.24</td><td>63.52 ± 5.78</td><td>49.16 ± 8.37</td></tr><tr><td>D-Limonene</td><td>10.36</td><td>9.71 ± 1.38</td><td>19.62 ± 4.13</td><td>14.77 ± 1.53</td><td>19.79 ± 4.63</td><td>60.74 ± 14.19</td><td>16.18 ± 5.76</td><td>34.08 ± 6.72</td><td>79.74 ± 12.15</td><td>74.22 ± 9.23</td></tr><tr><td>β-myrcene</td><td>9.46</td><td>1.91 ± 1.05</td><td>4.38 ± 1.02</td><td>3.43 ± 1.19</td><td>12.53 ± 2.52</td><td>44.03 ± 12.17</td><td>25.94 ± 5.84</td><td>26.07 ± 2.44</td><td>50.97 ± 7.38</td><td>44.39 ± 4.12</td></tr><tr><td>Copaene</td><td>17.39</td><td>10.92 ± 2.15</td><td>29.06 ± 5.95</td><td>21.38 ± 2.57</td><td>4.14 ± 0.93</td><td>38.62 ± 9.59</td><td>23.16 ± 4.24</td><td>7.20 ± 1.94</td><td>9.03 ± 1.23</td><td>7.37 ± 0.57</td></tr><tr><td>Ocimene</td><td>10.83</td><td>0.46 ± 0.07</td><td>0.73 ± 0.50</td><td>0.58 ± 0.03</td><td>16.55 ± 2.42</td><td>24.04 ± 4.16</td><td>19.50 ± 3.43</td><td>14.59 ± 1.47</td><td>21.05 ± 3.46</td><td>16.23 ± 4.08</td></tr><tr><td>β-trans-ocimene</td><td>10.57</td><td>0.42 ± 0.07</td><td>1.24 ± 0.37</td><td>0.87 ± 0.13</td><td>5.41 ± 0.47</td><td>18.24 ± 1.73</td><td>12.02 ± 1.00</td><td>6.28 ± 0.90</td><td>15.26 ± 3.26</td><td>11.97 ± 2.47</td></tr><tr><td>γ-terpinen</td><td>11.05</td><td>0.56 ± 0.03</td><td>1.50 ± 0.25</td><td>0.92 ± 0.15</td><td>0.45 ± 0.25</td><td>7.28 ± 1.26</td><td>5.09 ± 0.96</td><td>1.45 ± 1.21</td><td>8.39 ± 0.60</td><td>6.26 ± 0.16</td></tr><tr><td>Terpelene</td><td>10.06</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>5.71 ± 0.50</td><td>3.09 ± 0.35</td><td>ND</td><td>6.18 ± 0.53</td><td>3.62 ± 0.45</td></tr><tr><td>Terpinolen</td><td>11.72</td><td>ND</td><td>1.12 ± 0.15</td><td>ND</td><td>ND</td><td>7.12 ± 2.15</td><td>3.21 ± 0.59</td><td>1.60 ± 0.23</td><td>11.88 ± 1.61</td><td>6.84 ± 1.81</td></tr><tr><td>α-phellandrene</td><td>9.74</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>5.10 ± 0.98</td><td>2.97 ± 0.69</td><td>3.82 ± 0.58</td><td>6.48 ± 0.93</td><td>4.96 ± 0.36</td></tr><tr><td>β-cyclocitral</td><td>14.15</td><td>1.02 ± 0.25</td><td>1.99 ± 0.23</td><td>1.56 ± 0.20</td><td>1.30 ± 0.30</td><td>1.84 ± 0.14</td><td>1.52 ± 0.35</td><td>1.35 ± 0.48</td><td>2.68 ± 0.51</td><td>2.08 ± 0.38</td></tr></tbody></table>

Note: The data are expressed as the means ± standard deviation of triplicate samples. RT retention time, ND not detectable
increasing the anthocyanin content, overexpression of *SIMYB75* also led to elevated levels of sugar, phenolic and flavonoid compounds.

## Overexpression of *SIMYB75* affects volatile aroma metabolism

Considering that some of the DEGs such as *LOX*, *AADC* and *CXE1* have been reported to play critical roles in aroma volatile production³³,³⁴, we performed a comparative analysis of aroma volatiles in WT and *SIMYB75*-OE tomatoes. Nearly 380 volatile compounds were identified in the Micro-Tom tomato fruits analyzed in this experiment, and only 38 characteristic volatiles were screened and quantified (Table 2). Based on the chemical nature of their precursors, tomato volatiles could be categorized into the following three groups: those derived from fatty acids, carotenoids and amino acids³⁴. As shown in Table 2, the aroma compounds derived from fatty acids, especially aldehydes, accounted for a large part of the tomato aroma content. The amount of aroma volatiles in fruits at the BR + 5 stage was much lower than the levels at the BR + 7 and BR + 10 stages, and some of the aroma volatiles could not even be detected at the early ripening stage. Most of the volatiles showed significantly higher levels in *SIMYB75*-OE tomatoes than those in the WT, including trans-2-hexenal, hexanal, 3-hexenal, benzaldehyde, 3-nonanone, methyl salicylate and geranyl acetate. By contrast, some volatiles such as 2-propyl-1-pentanol and 2-ethylhexyl acrylate showed similar or lower contents in *SIMYB75*-OE tomatoes compared to the WT. Noticeably, the contents of terpenoids/terpenes, terpineol, (-)-4-terpineol, copaene, γ-terpinen and β-myrcene in the *SIMYB75*-OE fruits exhibited a dramatic increase, with some of them having more than 10 times higher levels than those in the WT, while the content in the WT stayed very low. Overall, the general trend is that volatiles in *SIMYB75*-OE tomatoes are significantly higher than they are in the WT, in particular those corresponding to aldehydes and terpenes.

## Analysis of the regulatory relationships between *SIMYB75* TF and other genes

To explore whether DEGs could serve as direct targets of *SIMYB75*, we checked the ability of this MYB TF to bind the promoters of some of the DEGs directly using a yeast one-hybrid (Y1H) assay. Since many DEGs were identified in this study, we chose highly conserved elements rather than specific DEGs to verify which ones could be directly bound by *SIMYB75*. By analyzing the promoters with the New Place database tools, we identified eight conserved elements as the primary ones present in the DEG promoters (Supplementary Table S1). Tests on these eight conserved elements in Y1H assays indicated that only MYBPLANT (AAACCAACCC) and

*MYBPZM (ACCTACCC)* could be recognized by *SIMYB75* (Fig. 6a), suggesting they might be direct targets of this MYB gene. To verify whether *SIMYB75* could interact with native promoters of some DEGs, we performed dual-luciferase assays with genes involved in ethylene signaling and volatile aroma metabolism. The data showed that *SIMYB75* could only trans-activate the promoters of the *LOXC*, *AADC2* and *TPS* genes (Fig. 6b).

## Discussion

In the last period, both classical and biotechnology-based breeding methods have made major progress that opened new prospects for crop improvement³,¹⁸,³⁵. Anthocyanins are well known for their strong antioxidant and health-protection properties, but some of the most highly consumed fruits, such as tomatoes, are completely devoid of this compound. In the present study, we obtained remarkable purple tomatoes enriched with this health-promoting compound through the ectopic expression of a single *SIMYB75* TF. Strikingly, *SIMYB75*-OE fruits also exhibited higher phenolics and flavonoids contents and enhanced the production of aroma volatiles. Our results indicated that *SIMYB75* TF plays excellent roles in promoting fruit quality, and it can be directly used in genetic engineering or traditional breeding methods to produce high quality fruit.

The purple phenotype of *SIMYB75*-OE fruits in this study is remarkable, and its anthocyanin content could increase to 1.86 mg g⁻¹ fresh weight (Fig. 2), which exceeded the results of earlier studies through over-expression of a single gene and previously required at least two genes to reach this level³,¹³,¹⁴,¹⁸,²³,²⁴. Given that it is difficult to find the parental lines that could highly express these two or three genes at the same time, together with its easy segregation characteristic, multi-gene strategies are not the best choice for breeders. Thus, expressing a single *SIMYB75* provides a straightforward strategy towards improving the fruit nutritional and sensory quality. In addition, the content of phenolic and flavonoid compounds was also greatly enhanced in *SIMYB75*-OE fruits, which is similar to the results of a study on anthocyanin-producing grapes³⁶. Phenolics, flavonoids, and anthocyains are all derived from phenylalanine, and the transcriptomic data showed that most of the genes involved in the phenylalanine pathway were upregulated in *SIMYB75*-OE fruits. Thus, *SIMYB75* might participate in regulating different branches of the phenylpropanoid pathway.

It has been shown that there was no difference in the seed sizes between the WT and *SLAN2*-expression lines²⁴. However, our results indicated the seeds of *SIMYB75*-OE tomato were significantly smaller than the WT ones (Fig. 1d). Previous studies indicated that plants that produce small seeds usually exhibit superior colonization
<figure><img src="image_8.png" /><figcaption>**Fig. 6 Regulatory roles of SIMYB75 on the conserved cis-elements and promoters of several DEGs. a,** The interactions between SIMYB75 and eight conserved cis-elements were determined by the yeast one-hybrid assay. No auto-activation ability was detected in the eight conserved cis-elements in yeast grown on SD/ura medium containing aureobasidin A (-Ura + AbA). The interactions were determined on SD/leu medium in the presence of AbA (-Leu + AbA). **b,** The transcription activation ability of SIMYB75 on the promoters of several DEGs was tested using dual luciferase assays. The LUC/REN ratio of the empty vector plus the promoter was used as the calibrator (set as 1). At least six biological replicates were conducted for each combination. The values are the means ± SD. \*P < 0.05; and \*\*P < 0.01 (Student's t-test)</figcaption></figure>

abilities³⁷, and thus *SIMYB75*-OE tomatoes may have a better dispersal capacity due to their abundant small seeds. In addition, it has been reported that *SLAN2* over-expression confers enhanced tolerance to abiotic stresses such as cold, high temperature, high light, and oxidative stress²¹⁻²³. This finding is consistent with our data showing that *SIMYB75* is also responsive to NaCl, ABA, and SA (Supplementary Fig. S2). Collectively, these traits might confer a better capacity to survive under harsh environmental conditions to *SIMYB75*-OE tomatoes.

While ethylene is known to regulate fruit ripening positively through inducing ripening-associated genes and autocatalytic ethylene production³⁸, two seemingly conflicting features of *SIMYB75*-OE tomatoes are their high level of ethylene production and delayed ripening (Fig. 2). As an attempt to explain the two seemingly inconsistent features, it may be speculated that the higher anti-oxidative properties of *SIMYB75*-OE tomatoes could mitigate the ethylene effect, since the ripening process comprises a series of physiological and biochemical changes, among which the anti-oxidative characteristics have been suggested to play a fundamental role³⁹,⁴⁰.

Aroma volatiles play a key role in the perception and acceptability of flowers, vegetables and fruits by consumers³⁴. This property leads aroma towards being the focus of horticultural research, but little is known about improving the aroma content through the transcriptional regulation of fruits⁴¹. Consistent with the study in petunia flowers⁴², our results showed that the contents of the phenylalanine-derived aroma volatiles benzaldehyde and methyl salicylate were significantly higher in *SIMYB75*-OE fruits. Furthermore, the transcript level of fatty acid-related genes (*LOXC*, *LOXB* and *LOXF*)⁴³ and the content of most fatty acid-derived aroma volatiles, especially trans-2-hexenal and (z)-2-heptenal, which are important contributors to flavor and consumer acceptance of tomato fruits⁴⁴, were greatly increased in *SIMYB75*-OE fruits. The primary contributors to aroma volatiles in blood oranges ranked from high to low are as follows: terpenes, esters and then aldehydes⁴⁵⁻⁴⁷, but tomato fruits have very low levels of terpenes and contain only minute amounts of accumulated monoterpenes⁴⁸. Interestingly, the terpene levels (Table 2) and transcript level of a terpene synthase gene (*Solyc06g060180.2*) (Table 1) were also greatly upregulated
in *SIMYB75*-OE tomato fruits, suggesting that the isoprenoid pathway was also affected in *SIMYB75*-OE tomatoes. Studies on tomato mutants and rice indicated that there are interactions between the phenylpropanoid and isoprenoid pathways, but the mechanisms underlying the interactions responsible for the production of specific metabolites are still unclear<sup>49,50</sup>. To explore the molecular mechanism of *SIMYB75* in promoting aroma volatiles, yeast one-hybrid and dual-luciferase assays were performed. The yeast one-hybrid assays revealed positive interactions between *SIMYB75* and conserved MYBPLANT (AAAC-CAACCC) and MYBPZM (ACCTACCC) elements, supporting the idea that the CC(T/A)ACC sequence motif might be the core binding site of *SIMYB75* (Fig. 6a), in a similar way with that of AtMYB12 as reported previously<sup>18</sup>. Consistent with this finding, the dual-luciferase assays showed that *LOXC*, *AADC2* and *TPS* genes with promoters containing MYBPLANT (AAACCAACCC) and MYBPZM (ACCTACCC) elements (Supplementary Fig. S5) can be directly trans-activated by *SIMYB75* (Fig. 6b). Thus, *SIMYB75* has a regulatory function and promotes aroma volatile accumulation, probably by targeting downstream aroma volatile-related genes.

Overall, the outcome of the study indicates that over-expression of a single *SIMYB75* TF can improve several tomato quality traits, including the sensory and nutritional aspects. This finding provides new opportunities for innovative breeding strategies aimed at generating tomato lines that better meet consumer requirement in terms of enriching health-promoting metabolites while better coping with a changing environment.

**Acknowledgements**

This work was supported by the National Key Research and Development Program (2016YFD0400101), the National Natural Science Foundation of China (31572175, 31772370), the Fundamental Research Funds for the Central Universities (2018CDXYSM0021) and the Committee of Science and Technology of Chongqing (cstckjcxljrc15).

**Author details**

¹Key Laboratory of Plant Hormones and Development Regulation of Chongqing, School of Life Sciences, Chongqing University, 401331 Chongqing, China. ²College of Resources, Sichuan Agricultural University, 611130 Chengdu, China. ³College of Horticulture and Landscape Architecture, Southwest University, 400716 Chongqing, China. ⁴Horticulture Institute, Sichuan Academy of Agricultural Sciences, 610066 Chengdu, China. ⁵Institute of Special Plants, Chongqing University of Arts and Sciences, 402160 Yongchuan, China. ⁶INRA, Genomique et Biotechnologie des Fruits, Chemin de Borde Rouge, F-31326 Castanet-Tolosan, France

**Conflict of interest**

The authors declare that they have no conflict of interest.

**Publisher's note**

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

**Supplementary Information** accompanies this paper at (https://doi.org/10.1038/s41438-018-0098-y).

Received: 24 May 2018 Revised: 1 September 2018 Accepted: 17 September 2018

Published online: 01 February 2019

**References**

1. Carrari, F. Metabolic regulation underlying tomato fruit development. *J. Exp. Bot.* **57**, 1883–1897 (2006).

2. Yokozawa, T. et al. Study on the inhibitory effect of tannins and flavonoids against the 1,1-diphenyl-2 picrylhydrazyl radical. *Biochem. Pharmacol.* **56**, 213–222 (1998).

3. Butelli, E. et al. Enrichment of tomato fruit with health-promoting anthocyanins by expression of select transcription factors. *Nat. Biotech.* **26**, 1301–1308 (2008).

4. Winkel-Shirley, B. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. *Plant Physiol.* **126**, 485–493 (2001).

5. Solfanelli, C., Poggi, A., Loreti, E., Alpi, A. & Perata, P. Sucrose-Specific Induction of the Anthocyanin Biosynthetic Pathway in Arabidopsis. *Plant Physiol.* **140**, 637–646 (2006).

6. Allan, A. C., Hellens, R. P. & Laing, W. A. MYB transcription factors that colour our fruit. *Trends Plant. Sci.* **13**, 99–102 (2008).

7. Gonzalez, A., Zhao, M., Leavitt, J. M. & Lloyd, A. M. Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings. *Plant J.* **53**, 814–827 (2008).

8. Ramsay, N. A. & Glover, B. J. MYB-bHLH-WD40 protein complex and the evolution of cellular diversity. *Trends Plant. Sci.* **10**, 63–70 (2005).

9. Gao, Y. et al. Tomato SIAN11 regulates flavonoid biosynthesis and seed dormancy by interaction with bHLH proteins but not with MYB proteins. *Hortic. Res.* **5**, https://doi.org/10.1038/s41438-018-0032-3 (2018).

10. Mehtens, F. The Arabidopsis Transcription Factor MYB12 Is a Flavonol-Specific Regulator of Phenylpropanoid Biosynthesis. *Plant Physiol.* **138**, 1083–1096 (2005).

11. Takos, A. M. et al. Light-Induced Expression of a MYB Gene Regulates Anthocyanin Biosynthesis in Red Apples. *Plant Physiol.* **142**, 1216–1232 (2006).

12. Muir, S. R. et al. Overexpression of petunia chalcone isomerase in tomato results in fruit containing increased levels of flavonols. *Nat. Biotech.* **19**, 470–474 (2001).

13. Jones, C., Mes, P. & Myers, J. Characterization and inheritance of the Anthocyanin fruit (Aft) tomato. *J. Hered.* **94**, 449–456 (2003).

14. Sapir, M. et al. Molecular aspects of Anthocyanin fruit tomato in relation to high pigment-1. *J. Hered.* **99**, 292–303 (2008).

15. Davuluri, G. R. et al. Fruit-specific RNAi-mediated suppression of DET1 enhances carotenoid and flavonoid content in tomatoes. *Nat. Biotech.* **23**, 890–895 (2005).

16. Bovy, A. High-flavonol tomatoes resulting from the heterologous expression of the maize transcription factor genes Lc and C1. *Plant Cell* **14**, 2509–2526 (2002).

17. Luo, J. et al. AtMYB12 regulates caffeoyl quinic acid and flavonol synthesis in tomato: expression in fruit results in very high levels of both types of polyphenol. *Plant J.* **56**, 316–326 (2008).

18. Zhang, Y. et al. Multi-level engineering facilitates the production of phenylpropanoid compounds in tomato. *Nat. Commun.* **6**, 8635 (2015).

19. Wang, J. H. et al. Anthocyanin Biosynthesis Regulation in the Fruit of Citrus sinensis cv. Tarocco. *Plant. Mol. Biol. Rep.* **34**, 1–13 (2016).

20. Petroni, K. & Tonelli, C. Recent advances on the regulation of anthocyanin synthesis in reproductive organs. *Plant Sci.* **181**, 219–229 (2011).

21. Meng, X. et al. An R2R3-MYB gene, LeAN2, positively regulated the thermotolerance in transgenic tomato. *J. Plant. Physiol.* **175**, 1–8 (2015).

22. Meng, X. et al. Overexpression of R2R3-MYB gene leads to accumulation of anthocyanin and enhanced resistance to chilling and oxidative stress. *Biol. Plant.* **58**, 121–130 (2014).

23. Kiferle, C. et al. Tomato R2R3-MYB Proteins SIANT1 and SIAN2: Same Protein Activity, Different Roles. *PLoS ONE* **10**, 1–20 (2015).

24. Meng, X. et al. Physiological changes in fruit ripening caused by over-expression of tomato SIAN2, an R2R3-MYB factor. *Plant. Physiol. Biochem.* **89**, 24–30 (2015).

25. Hu, D.-G. et al. MdMYB1 Regulates Anthocyanin and Malate Accumulation by Directly Facilitating Their Transport into Vacuoles in Apples. *Plant Physiol.* **170**, 1315–1330 (2016).
26. Lim, W., Miller, R., Park, J. & Park, S. Consumer sensory analysis of high flavonoid transgenic tomatoes. *J. Food Sci.* **79**, S1212–S1217 (2014).

27. Huang, W. et al. Overexpression of a tomato miR171 target gene SIGRAS24 impacts multiple agronomical traits via regulating gibberellin and auxin homeostasis. *Plant Biotechnol.* **15**, 472–488 (2017).

28. Yang, L. et al. Silencing of SIPL (Solyc03g111690), which encodes a pectate lyase in tomato, confers enhanced fruit firmness, prolonged shelf-life, and reduced susceptibility to gray mold. *Plant Biotechnol.* **15**, 1544–1555 (2017).

29. Liu, M. et al. The chimeric repressor version of an Ethylene Response Factor (ERF) family member, Si-ERF. B3, shows contrasting effects on tomato fruit ripening. *New Phytol.* **203**, 206–218 (2014).

30. Zheng, H., Zhang, Q., Quan, J., Zheng, Q. & Xi, W. Determination of sugars, organic acids, aroma components, and carotenoids in grapefruit pulps. *Food Chem.* **205**, 112–121 (2016).

31. Zhang, Y. et al. Phenolic compositions and antioxidant capacities of Chinese wild mandarin (Citrus reticulata Blanco) fruits. *Food Chem.* **145**, 674–680 (2014).

32. Ye, Y.-J. et al. Banana fruit VQ motif-containing protein5 represses cold-responsive transcription factor MaWRKY26 involved in the regulation of JA biosynthetic genes. *Sci. Rep.* **6**, 23632 (2016).

33. Klee, H. J. & Tieman, D. M. Genetic challenges of flavor improvement in tomato. *Trends Genet.* **29**, 257–262 (2013).

34. Wang, L., Baldwin, E. A. & Bai, J. Recent Advance in Aromatic Volatile Research in Tomato Fruit: The Metabolisms and Regulations. *Food Bioprocess Technol.* **9**, 203–216 (2016).

35. Sagor, G. H. M. et al. A novel strategy to produce sweeter tomato fruits with high sugar contents by fruit-specific expression of a single bZIP transcription factor gene. *Plant. Biotechnol. J.* **14**, 1116–1126 (2015).

36. Deluc, L. et al. Characterization of a grapevine R2R3-MYB transcription factor that regulates the phenylpropanoid pathway. *Plant Physiol.* **140**, 499–511 (2006).

37. Hotta, O. C. & Tanksley, S. D. Natural variation in an ABC transporter gene associated with seed size evolution in tomato species. *PLoS Genet.* **5**, 1–12 (2009).

38. Liu, M., Pirrello, J., Chervin, C., Roustan, J. P. & Bouzayen, M. Ethylene control of fruit ripening: revisiting the complex network of transcriptional regulation. *Plant Physiol.* **169**, 2380–2390 (2015).

39. Conde, C., Delrot, S. & Gerós, H. Physiological, biochemical and molecular changes occurring during olive development and ripening. *J. Plant. Physiol.* **165**, 1545–1562 (2008).

40. Jimenez, A. et al. Changes in oxidative processes and components of the antioxidant system during tomato fruit ripening. *Planta* **214**, 751–758 (2002).

41. Klee, H. J. Improving the flavor of fresh fruits: genomics, biochemistry, and biotechnology. *New Phytol.* **187**, 44–56 (2010).

42. Spitzer-Rimon, B. et al. EOBI1, a Gene Encoding a Flower-Specific Regulator of Phenylpropanoid Volatiles' Biosynthesis in Petunia. *Plant Cell* **22**, 1961–1976 (2010).

43. Shen, J. et al. A 13-lipoxygenase, TomloxC, is essential for synthesis of C5 flavour volatiles in tomato. *J. Exp. Bot.* **65**, 419–428 (2014).

44. Tieman, D. et al. A chemical genetic roadmap to improved tomato flavor. *Science* **355**, 391–394 (2017).

45. Arena, E., Guarrera, N., Campisi, S. & Asmundo, C. N. Comparison of odour active compounds detected by gas-chromatography-olfactometry between hand-squeezed juices from different orange varieties. *Food Chem.* **98**, 59–63 (2006).

46. Maccarone, E., Campisi, S., Fallico, B., Rapisarda, P. & Sgarlata, R. Flavor Components of Italian Orange Juices. *J. Agric. Food Chem.* **46**, 2293–2298 (1998).

47. Moufida, S. d. & Marzouk, B. Biochemical characterization of blood orange, sweet orange, lemon, bergamot and bitter orange. *Phytochemistry* **62**, 1283–1289 (2003).

48. Davidovich-Rikanati, R. et al. Enrichment of tomato flavor by diversion of the early plastidial terpenoid pathway. *Nat. Biotech.* **25**, 899–901 (2007).

49. Enfissi, E. M. A. et al. Integrative transcript and metabolite analysis of nutritionally enhanced DE-ETIOLATED1 downregulated tomato fruit. *Plant Cell* **22**, 1190–1215 (2010).

50. Kim, J. K. et al. Variation and correlation analysis of flavonoids and carotenoids in Korean pigmented rice (Oryza sativa L.) cultivars. *J. Agric. Food Chem.* **58**, 12804–12809 (2010).