# A systems biology study in tomato fruit reveals correlations between the ascorbate pool and genes involved in ribosome biogenesis, translation, and the heat-shock response

Rebecca Stevens, Pierre Baldet, Jean-Paul Bouchet, Mathilde Causse, Catherine Deborde, Claire Deschodt, Mireille Faurobert, Cecile Garchery, Virginie Garcia, Hélène Gautier, et al.

► To cite this version:

Rebecca Stevens, Pierre Baldet, Jean-Paul Bouchet, Mathilde Causse, Catherine Deborde, et al. A systems biology study in tomato fruit reveals correlations between the ascorbate pool and genes involved in ribosome biogenesis, translation, and the heat-shock response. Frontiers in Plant Science, 2018, 9 (137), 16 p. 10.3389/fpls.2018.00137. hal-01721265

HAL Id: hal-01721265

https://hal.science/hal-01721265

Submitted on 1 Mar 2018

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>
<figure><img src="image_2.png" /></figure>

Check for
updates

# A Systems Biology Study in Tomato Fruit Reveals Correlations between the Ascorbate Pool and Genes Involved in Ribosome Biogenesis, Translation, and the Heat-Shock Response

Rebecca G. Stevens¹*, Pierre Baldet², Jean-Paul Bouchet¹, Mathilde Causse¹, Catherine Deborde²,³, Claire Deschodt¹, Mireille Faurobert¹, Cécile Garchery¹, Virginie Garcia², Hélène Gautier⁴, Barbara Gouble⁵, Mickaël Maucourt²,³, Annick Moing²,³, David Page⁵, Johann Petit², Jean-Luc Poëssel¹, Vincent Truffault¹ and Christophe Rothan²

## OPEN ACCESS

**Edited by:**

Atsushi Fukushima,
RIKEN, Japan

**Reviewed by:**

Victoriano Valpuesta,
University of Málaga, Spain
Jedrzej Jakub Szymanski,
Weizmann Institute of Science, Israel

**\*Correspondence:**

Rebecca G. Stevens
Rebecca.stevens@inra.fr

**Specialty section:**

This article was submitted to Plant Systems and Synthetic Biology, a section of the journal Frontiers in Plant Science

**Received:** 09 October 2017

**Accepted:** 24 January 2018

**Published:** 14 February 2018

**Citation:**

Stevens RG, Baldet P, Bouchet J-P, Causse M, Deborde C, Deschodt C, Faurobert M, Garchery C, Garcia V, Gautier H, Gouble B, Maucourt M, Moing A, Page D, Petit J, Poëssel J-L, Truffault V and Rothan C (2018) *A Systems Biology Study in Tomato Fruit Reveals Correlations between the Ascorbate Pool and Genes Involved in Ribosome Biogenesis, Translation, and the Heat-Shock Response*. *Front. Plant Sci.* 9:137. doi: 10.3389/fpls.2018.00137

¹Institut National de la Recherche Agronomique, UR1052, Génétique et Amélioration des Fruits et Légumes, Montfavet, France, ²Institut National de la Recherche Agronomique, Université de Bordeaux, UMR1332, Biologie du Fruit et Pathologie, Villenave d'Ornon, France, ³Plateforme Métabolome du Centre de Génomique Fonctionnelle Bordeaux, Centre Institut National de la Recherche Agronomique de Bordeaux, Villenave d'Ornon, France, ⁴Institut National de la Recherche Agronomique, UR1115, Plantes et Systèmes de culture Horticoles, Avignon, France, ⁵Institut National de la Recherche Agronomique, Université d'Avignon et des Pays du Vaucluse, UMR408 Sécurité et Qualité des Produits d'Origine Végétale, Avignon, France

Changing the balance between ascorbate, monodehydroascorbate, and dehydroascorbate in plant cells by manipulating the activity of enzymes involved in ascorbate synthesis or recycling of oxidized and reduced forms leads to multiple phenotypes. A systems biology approach including network analysis of the transcriptome, proteome and metabolites of RNAi lines for ascorbate oxidase, monodehydroascorbate reductase and galactonolactone dehydrogenase has been carried out in orange fruit pericarp of tomato (*Solanum lycopersicum*). The transcriptome of the RNAi ascorbate oxidase lines is inversed compared to the monodehydroascorbate reductase and galactonolactone dehydrogenase lines. Differentially expressed genes are involved in ribosome biogenesis and translation. This transcriptome inversion is also seen in response to different stresses in Arabidopsis. The transcriptome response is not well correlated with the proteome which, with the metabolites, are correlated to the activity of the ascorbate redox enzymes—ascorbate oxidase and monodehydroascorbate reductase. Differentially accumulated proteins include metacaspase, protein disulphide isomerase, chaperone DnaK and carbonic anhydrase and the metabolites chlorogenic acid, dehydroascorbate and alanine. The hub genes identified from the network analysis are involved in signaling, the heat-shock response and ribosome biogenesis. The results from this study therefore reveal one or several putative signals from the ascorbate pool which modify the transcriptional response and elements downstream.

**Keywords:** ascorbate, cellular signaling, heat-shock response, redox, ribosome biogenesis, translation, tomato (*Solanum lycopersicum*)
INTRODUCTION

Tomato has emerged as a model fruit for studying the relationships between ascorbate and fruit physiology. In addition to its nutritional value, ascorbate is a major antioxidant that contributes to fruit tolerance to biotic and abiotic stresses (Malacrida et al., 2006; Davey et al., 2007; Stevens et al., 2008). Studying the interactions between the ascorbate pool and fruit ripening is of particular interest because the vast majority of fresh market tomatoes produced in Europe are harvested at orange or light-red stages, before being fully ripe. At these stages, fruit undergo major modifications in color, firmness, sweetness, acidity, and aroma. These changes are under developmental and hormonal control and imply profound alterations in the transcriptome, proteome and metabolome (Klee and Giovannoni, 2011). In a previous study, we showed that chilling stress conditions (,<10°C) that can prevail during postharvest transport and storage affect fruit firmness, and that this effect is correlated with the ascorbate oxidation status (Stevens et al., 2008). However, the mechanisms by which ascorbate impinges on fruit physiology remain to be deciphered.

Ripening fruit are non-photosynthetic sink tissues relying on import of sucrose for providing reducing power, ATP, and precursor molecules necessary for synthesis of fruit-specific metabolites or proteins. The conversion of sucrose into hexose-phosphates is the starting point for glycolysis, the TCA cycle and the oxidative pentose phosphate pathway which provide the tissue with its metabolic and energetic requirements. The fruit therefore represents a complex biological system built on multi-layered interactions between metabolites, proteins and genes which form a network characteristic of the tissue. Induced genetic modifications are a good way of perturbing the system by varying the level of a gene and its corresponding protein and/or metabolites (network nodes) to evaluate the new steady state of the network, thus pinpointing the interactions that a particular node will have (Bassel et al., 2012; Krouk et al., 2013). A highly connected node (a hub) will share a number of edges with other nodes and perturbation of a regulatory hub will often give rise to highly pleiotropic phenotypes. This sort of approach is particularly useful to study a biological question as it can generate new hypotheses and identify regulatory genes.

Ascorbate is a major metabolite in plants including fruits with many functions often associated with its role in electron donation and as a cofactor (Gest et al., 2013b). Interestingly, although many studies in plants focus on the photo-protective or antioxidant functions of ascorbate, studies in animals have revealed further roles for ascorbate, for example the molecule functions as a cofactor for the ten-eleven translocation enzymes, which can convert 5-methylcytosine, to 5-hydroxymethylcytosine (5hmC) during demethylation of DNA in mouse embryonic stem cells (Blaschke et al., 2013).

These enzymes are dioxygenases and therefore this is an example of a role for ascorbate as an electron donor to iron-dependent oxidoreductases. Ascorbate is also involved in oxidative protein folding in the endoplasmic reticulum and therefore plays a role in protein maturation (Banhegyi et al., 2003; Zito et al., 2012; Szarka and Lorincz, 2014). Dehydroascorbate has been shown to have a role in neuronal energy metabolism via activation of glucose-6-phosphate dehydrogenase thus increasing flux through the pentose phosphate pathway (Puskas et al., 2000; Cisternas et al., 2014). Furthermore, ascorbate has been identified as a kinase inhibitor (Carcamo et al., 2004) as well as an inhibitor of hexokinase (Fiorani et al., 2000) and its oxidized form can inhibit glyceraldehyde 3-phosphate dehydrogenase leading to an energetic crisis and cell death in cancer cells which uptake more dehydroascorbate via glucose transporters (Yun et al., 2015). These examples have not been shown in plants: the study of ascorbate in non-photosynthetic tissues may be key to revealing additional functions for the molecule; fruit is an ideal tissue, particularly in advanced ripening stages when photosynthetic activity is low.

Ascorbate's functions are profoundly affected by the oxidation state of the ascorbate pool: ascorbate is oxidized to the monodehydroascorbate radical, a relatively stable radical, which can regain an electron to regenerate ascorbate, or lose a second electron to become dehydroascorbate; these three forms are in equilibrium (Bors and Buettner, 1997). Plants also possess ascorbate oxidase enzymes which have roles in oxygen removal and signaling (De Tullio et al., 2007; De Tullio, 2012). Ascorbate oxidase activity has also been shown to have links with sugar metabolism and yield (Garchery et al., 2013), in a similar way to an isoform of monodehydroascorbate reductase (MDHAR), except that the phenotypes are inverted (Truffault et al., 2016): plants with reduced MDHAR activity show reduced growth and yield and lower levels of sugars (sucrose and hexose) particularly in leaves but counter-intuitively contain more ascorbate, particularly under light (Haroldsen et al., 2011; Gest et al., 2013a), whereas plants with reduced ascorbate oxidase (AO) activity have improved yield stability and more sucrose and hexose, particularly in leaves. A mitochondrial enzyme controls the final stage of ascorbate synthesis, galactonolactone dehydrogenase (GLD) (Alhagdow et al., 2007; Leferink et al., 2009), this enzyme is an integral part of complex I and is required for complex I function independently of ascorbate levels (Schimmeyer et al., 2016). Similarly to AO and MDHAR, manipulation of GLD levels in tomato RNAi lines leads to phenotypes including differences in ascorbate redox state, sugar metabolism, plant growth and fruit yield (Alhagdow et al., 2007).

Modification of the ascorbate pool via manipulation of the activity of these enzymes has therefore generated many pleiotropic phenotypes at the physiological level. We hypothesize that the enzymes are affecting the levels of one or several metabolites, closely linked to their activity, which can be seen as a hub, involved in many key cellular processes. We have therefore carried out a large-scale study including phenotypic, metabolic, transcriptomic and proteomic data for orange fruit of RNAi lines for the three enzymes (AO, GLD, and MDHAR) shown in Figure 1 in order to identify their links with cellular

<footnote>**Abbreviations:** RNAi, RNA interference; GLD, galactonolactone dehydrogenase; AO, ascorbate oxidase; MDHAR, monodehydroascorbate reductase; GO, Gene Ontology; FDR, false discovery rate; WGCNA, Weighted gene correlation network analysis; WT, wild-type; TCA, tricarboxylic acid cycle; HSP70, heat-shock protein family; ATP, adenosine triphosphate.</footnote>
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Biological system showing the enzymes studied. Schematic diagram of the enzymes studied in tomato using RNAi. GLD, galactonolactone dehydrogenase, a mitochondrial enzyme; AO, ascorbate oxidase, probably apoplastic; MDHAR, monodehydroascorbate reductase, a cytosolic and peroxisomal isoform. Ascorbate is generated in the mitochondrion and is transported into the cell via an unknown mechanism. Ascorbate oxidase oxidizes ascorbate to monodehydroascorbate which will disproportionate into dehydroascorbate and ascorbate. Monodehydroascorbate is recycled to ascorbate by enzymes with monodehydroascorbate reductase activity.</figcaption></figure>

metabolism and to correlate changes in gene expression, proteins and metabolites with the activity of these enzymes. A link is provided between the redox state of the ascorbate pool, cellular protein synthesis and stability, and ribosomal function.

# MATERIALS AND METHODS

## Greenhouse Plant Culture

West Virginia 106 cherry tomato plants were grown in a multi-span Venlo-type greenhouse in 5l pots (potting compost P3 Tref, Tref EGO substrates BV). Water and nutrients were supplied to the plants using a drip irrigation system to maintain 20-30% drainage. Light intensities of 300-700 photosynthetically active radiation were obtained over the culture period. Flowers were mechanically pollinated three times a week and side shoots removed as they appeared. The 35S-RNAi lines used have been previously described [AO, Solyc04g054690 (Garchery et al., 2013), GLD, Solyc10g079470 (Alhagdow et al., 2007), MDHAR (or MD), Solyc09g009390 (Gest et al., 2013a)]. In each case the specific transcript levels and enzyme activity were checked and shown to decrease. The individual lines chosen were the following: GLD5.3, AO15.1, and MD5 with a wild-type (WT)

untransformed control. Plants homozygous for the transgene from the T2 generation were used. A total of 15 plants per line were used, 10 for fruit sampling and 5 for destructive phenotypic measurements. All 15 plants were used for non-destructive phenotypic measurements.

## Pre- and Post-harvest Phenotypic Measurements

Plant height from the pot to the top of the apical meristem was measured weekly for 7 weeks on 10 plants per line. The growth rate was calculated as cm gained per day. The average leaf area of a fully expanded leaf (one leaf per plant, 10 plants) was calculated by scanning the leaf and calculating the area of the image using ImageJ software. Fruit characteristics (number of days between anthesis and orange maturity stage, and fruit weight) were calculated on at least six fruits per plant from 10 plants. Fruit weight, diameter, firmness measurements (compression and maximal force for skin rupture) and color were measured on a minimum of 20 fruits per line immediately after harvesting (ripe fruits). The weight loss of seven fruits per line, related to conservation properties, was measured after 5 days post-harvest at 20°C. Texture measurements were performed using a multi-purpose texturometer (TaPlus, Lloyd). The global firmness of each fruit was measured by a uniaxial compression test with a 50 mm flat disk (test-speed 10 cm.min⁻¹) and expressed as the pressure (kPa) required for a 3% deformation of fruit height. The maximal force (N) corresponding to skin rupture was recorded using a flat cylinder tip (2 mm diameter) for a 7 mm penetration into the equatorial part of the fruit (test-speed 20 mm.min⁻¹). The color was determined using a CR-400 chromameter (Minolta Co. Ltd., Osaka, Japan) and expressed in the CIE 1976 L*a*b* Color Space coordinates (illuminant D65, 0_view angle, illumination area diameter 8 mm). The chroma a* axis is a good indicator for tomato ripening, varying from green (−60) to red (+60).

Statistical analysis on phenotypic characteristics was carried out using XLstat (Addinsoft, Paris, France) using a non-parametric comparison of means (Kruskal Wallis) with correction (Dunn). Different letters in the tables indicate significant differences (5% significance level).

## Experimental Design for Transcriptome and Metabolome Analyses

Fruit were tagged at anthesis and harvested 20 days later (for ascorbate measurements and proteome analysis only) or when orange. Fruit were harvested during the growing season (from September to November) from trusses 2, 3, and 4 from 10 plants using three fruit/truss/plant to make 3 pools of 30 fruits per line. Fruit pericarp including the epidermis was retained, frozen in liquid nitrogen and stored at −80°C before being ground to a fine powder in liquid nitrogen for the following analyses.

## NMR Metabolomic Profiling

Polar metabolites were extracted from lyophilized powder with a series of hot ethanol/water solutions and quantified by ¹H-NMR as previously described (Mounet et al., 2007). For extract preparation and NMR acquisition parameters, special care was
taken to allow absolute quantification of individual metabolites through addition of 5 mM ethylene diamine tetra-acetic acid disodium salt solution to improve the NMR signal resolution and quantification of organic acids, adequate choice of the NMR acquisition parameters, and use through an electronic reference for quantification (Moing et al., 2004). Unknown compounds were quantified in arbitrary units. The lyophilized extracts were titrated with KOD to apparent pH 6 in 400 mM potassium phosphate buffer solution in D₂O and lyophilized again. The ¹H-NMR spectra were recorded on each dried titrated extract solubilized in 0.5 mL D₂O with the sodium salt of (trimethylsilyl)propionic-2,2,3,3-d4 acid (TSP) at a final concentration of 0.01% for chemical shift calibration, at 500.162 MHz on an Avance spectrometer (Bruker BioSpin, Wissembourg, France) using a 5 mm inverse probe. Sixty-four scans of 32K data points were acquired with a spectral width of 6000 Hz, acquisition time of 2.73 s, 90° pulse angle and recycle delays of 25 s. The raw ¹H NMR spectral profiles have been deposited, with associated metadata, into the Metabolomics Repository of Bordeaux MeRy-B (http://services.cbib.u-bordeaux2.fr/MERYB/public/PublicREF.php?REF=T07001).

## Individual Metabolites

Ascorbate and dehydroascorbate concentrations in fruit powder were analyzed according to previous protocols based on the reduction of an iron-dipyridyl complex leading to a color change detectable spectrophotometrically with a filter at 525 nm (Stevens et al., 2006) with modifications (Garchery et al., 2013).

Individual polyphenols were analyzed using 200 mg of lyophilised powder from the orange fruit (3 pools), to which of 25 µL of taxifolin (1 mg·mL⁻¹, Extrasynthèse, Genay France) was added as an internal standard. Three successive extractions with 70% ethanol were carried out before evaporation to dryness in a vacuum concentrator. The pellet was resuspended in 300 µL water, 700 µL methanol and filtered with a PTFE syringe filter. Ten micro liters were analyzed by reversed-phase HPLC (Shimadzu Prominence system coupled with a Diode Array Detector (HPLC-DAD), Merck LichroCART® 250-4 column; phase Superspher® 100 RP-18 end capped 4 µm) using a gradient elution with ortho-phosphoric acidified water (pH 2.6) (A) and methanol (B) as solvents. The program was as follows: 0 min 15% B, 10 min 15% B, 20 min 21% B, 32 min 21% B, 42 min 23% B, 52 min 23% B, 107 min 60% B; flow rate 0.5 mL·min⁻¹ from 0 to 44 min, 0.5 mL·min⁻¹ to 0.7 mL·min⁻¹ from 44 to 46 min, 0.7 mL·min⁻¹ from 46 to 60 min. Phenolic compounds were monitored at 280 nm. Chlorogenic acid and rutin were quantified using external calibration curves of standard compounds (Extrasynthèse, Genay, France); other compounds were expressed in arbitrary units.

## Proteome Analysis

Protein extraction, two-dimensional gel electrophoresis, and protein identification and classification were carried out as previously described (Xu et al., 2013). Briefly, proteins were extracted using the phenol extraction method developed by Faurobert et al. (2007). Proteins were then separated by two-dimensional gel electrophoresis. After Coomassie colloidal

staining, image analysis was performed with Samespot software (version 4.1), and the normalized spot volumes were obtained. Of the 108 spots revealed, 49 were eliminated as not being variable. Protein identification of the 59 variable spots was performed at the proteome platform of Le Moulon (Gif-sur-Yvette) using a nano-LC-MS/MS method following a previously described procedure (Xu et al., 2013). The database search was run against the tomato unigene database version 4 held by the Solanaceae Genome Network (http://solgenomics.net/; Tomato_200607_build_2) with X!Tandem software (http://www.thegpm.org/TANDEM/, version 2010.12.01.1). Sequences were subsequently converted to Solycs (ITAG version 2.4) by Blast analysis. This conversion was confirmed by re-analyzing the mass spectra or by searching for similarity between the protein peptide sequences and the ITAG2.4 version of the genome by TBLASTN.

## Microarray Analysis and Comparison of Transcriptome Data

Microarray experiments were performed on the three biological replicates (3 pools) and technical duplicates with the dyes reversed (dye swap). The TOM2 cDNA microarray [CGE, Boyce Thompson Institute, Cornell University, http://bti.cornell.edu/CGEP/CGEP.html, (Fei et al., 2011)] containing 12,160 70-mer oligonucleotides corresponding to 11,862 unique and randomly selected transcripts of the tomato genome were used. Sample preparation and microarray hybridization and analysis were carried out as described (Mounet et al., 2009). Briefly, one µg of DNA-free total RNA was amplified with the MessageAmpTM aRNA Ambion Kit (Ambion/Applied Biosystems), labeled using the CyScribe Post-Labeling Kit (GE Healthcare) and the cy3- and cy5-labeled cDNAs were hybridized to TOM2 microarrays using an automatic hybridization station HS 4800 (Tecan) as described. Microarray slides were scanned with a Genepix 4000 B fluorescence reader (Axon Instruments) using Genepix 4.0 Pro image acquisition software. The photomultiplier tube voltage was adjusted to 640 V for Cy3 and 700 V for Cy5. Spot flagging was carried out by Genepix (missing spots) and by visual inspection of the images in order to exclude saturated and heterogeneous spots. The raw data were then submitted to LIMMA R Package v2.12.0 (Ritchie et al., 2015) for data visualization, normalization as described (Alhagdow et al., 2007). Background was first subtracted from the raw median intensity values using the normexp function. Two normalization algorithms were successively applied to the raw data (Yang et al., 2002): print-tip loess (default parameters) as the within-array normalization, followed by a quantile normalization (default parameters) as between-array normalization. The normalized data were then analyzed as described below. To convert to current Solyc identifiers, each of the TOM2 probes was blasted against the tomato genome (ITAG2.3, minimum alignment length 40 bp, minimum identity 84%), which had been sequenced in the meantime, and following this 7633 unique Solycs were retained.

## R/maanova

Differential expression and transcriptome comparisons were carried out using the R/maanova package [https://www.bioconductor.org/packages/release/bioc/html/maanova.html (R
Development-Core-Team, 2015)]. Means and standard errors were calculated for each gene in each line. A simple model taking into account the genotype (line) was used. A matrix of comparisons was used for the pairwise comparison of the transcriptome datasets first of each transgenic line with WT and then with each other. For each comparison, a threshold for the FDR corrected *p*-value of 0.001 was applied. Following initial tests, we added a sub-hypothesis to the comparison matrix to test whether the AO-GLD transcriptome and the AO-MDHAR transcriptomes were inverted i.e., the wild-type expression was the average of the gene expression in the two transgenic lines being compared. Expression differences were expressed as log2 of the Fold Change between the two lines. The linear model was tested using 1,000 iterations and then an FDR adjustment of the *p*-values obtained. A threshold for the corrected *p*-value of 0.001 was set to identify the core gene set for the network visualization.

## Weighted-Gene-Correlation Network Analysis

Weighted gene correlation network analysis (WGCNA) was carried out using the normalized expression data for the 7633 genes with a “trait” file including the 22 metabolites and 40 proteins using the WGCNA package in “R” (Langfelder and Horvath, 2008). No outliers were detected after sample clustering so the whole dataset was used. A soft threshold of 14 and minimum module size of 25 were set. Networks were visualized using Cytoscape version 3.2.0 (Cline et al., 2007) (www.cytoscape.org) using the node and edge output file from the WGCNA package with a threshold set to 0.52 (to fix a number of nodes and edges suitable for visualization). The dataset used contained the 182 differentially expressed genes (differentially expressed in the pairwise comparisons of the three types of transgenic lines), the 22 metabolites and 40 proteins. The number of genes was limited to facilitate network creation and to give weight to the proteins and metabolites. The layout was first designated as degree sorted circle layout and then a degree layout was chosen to aid node identification. Node size is proportional to the degree of connectivity.

## Gene Ontology

Gene ontology (GO) enrichment on the whole transcriptome signatures generated from pairwise comparisons of the lines was carried out using the R package “goseq” (Young et al., 2010) with the corrected *p*-value false detection rate set to 5%. The package requires the log of fold-change and *p*-value from the R/maanova analysis. No length correction was required as the data was from microarrays. Tomato gene ontology for the 182 differentially expressed genes was carried out using the online tools available from the Boyce Thompson Institute, Cornell University (http://bioinfo.bti.cornell.edu/cgi-bin/MetGenMAP/home.cgi, Boyle et al., 2004).

## ROSMETER

The ROSMETER bioinformatic tool allows evaluation of transcriptome signatures linked to reactive oxygen species (ROS) (Rosenwasser et al., 2013). The tool consists of a web-based interface whereby vector-based correlations between stored

transcriptome data, chosen for its ROS signature, and the provided data of interest are calculated. The stored data are microarray-based data from wild type and Arabidopsis mutants for the following genes: ascorbate peroxidase, the *flu* mutant, mitochondrial alternative oxidase, catalase 2 and superoxide dismutase subjected or not to the following stresses: ozone, high light, hydrogen peroxide, drought, rotenone, methyl viologen and 3-aminotriazole as described in the original research paper. The data from the RNAi lines normalized to wild-type (fold change and *p*-value) were compared with this series of transcriptomes arising from different conditions or mutations in *Arabidopsis thaliana*. In order to compare the transcriptomes, the tomato Solycs were first converted to their closest Arabidopsis homolog and submitted to the ROSMETER interface. In the heatmaps, red represents a positive correlation and green a negative one.

## RESULTS

### Transgenic Lines Show Slight Differences in Growth and Development of Plants and Fruit

The three RNAi transgenic lines (AO: ascorbate oxidase; GLD: galactonolactone dehydrogenase and MDHAR: monodehydroascorbate reductase) and wild-type plants grown in the greenhouse were subjected to simple phenotypic analyses to measure plant height, growth rate and final leaf area, the duration of fruit ripening (anthesis to orange fruit) and fruit weight were also measured. The results for the plant phenotypes are shown in **Table 1**: AO and GLD RNAi lines had a slower growth rate and a smaller final plant height compared to wild-type plants. Fruit from the same two transgenic lines also took longer to reach maturity (Table 2). Fruit from the MDHAR silenced lines were significantly smaller than fruit from the other lines, including wild-type (Tables 2, 3), as per previous results (Truffault et al., 2016). Fruits from the GLD RNAi line were also smaller in diameter (Table 3). Harvested fruits were analyzed for fruit firmness, color, and post-harvest weight loss (Table 3). No significant differences were found between the lines in terms of compression or color but the AO and MDHAR lines required less force to penetrate the skin compared to wild type. MDHAR RNAi fruit also showed greater post-harvest weight-loss compared to the WT, AO, and GLD lines.

### Metabolic Differences Are Seen in Phenolic and Amino Acid Metabolism

Metabolome analysis was carried out on the fruit pericarp samples using ¹H-NMR profiling of polar extracts. Ascorbate (Table 4) and major fruit polyphenols (Table 5) were analyzed by specific methods. Ascorbate was significantly increased in orange fruit pericarp of MDHAR lines as previously observed (Gest et al., 2013a) and dehydroascorbate was significantly different in AO and MDHAR lines, although this difference was the opposite of what would be expected given the activities of the enzymes; GLD and WT had intermediate values. Ascorbate data for fruit 20 days after anthesis is presented in **Supplementary Table S1**, unlike
<table><caption>**TABLE 1** | Plant phenotypes of transgenic lines and WT.</caption><thead><tr><th></th><th>Height after 7 days (cm)</th><th>SE</th><th>p</th><th>Height after 7 weeks (cm)</th><th>SE</th><th>p</th><th>Starting growth rate (cm/day)</th><th>SE</th><th>p</th><th>Final growth rate (cm/day)</th><th>SE</th><th>p</th><th>Average leaf area (cm²)</th><th>SE</th><th>p</th></tr></thead><tbody><tr><td>WT</td><td>26.2</td><td>3.7</td><td>a</td><td>193.0</td><td>9.0</td><td>b</td><td>2.7</td><td>0.2</td><td>b</td><td>5.84</td><td>0.32</td><td>b</td><td>313</td><td>46</td><td>a</td></tr><tr><td>AO</td><td>23.7</td><td>1.8</td><td>a</td><td>181.8</td><td>4.9</td><td>ab</td><td>2.5</td><td>0.2</td><td>ab</td><td>5.74</td><td>0.27</td><td>ab</td><td>296</td><td>49</td><td>a</td></tr><tr><td>GLD</td><td>24.5</td><td>2.3</td><td>a</td><td>171.3</td><td>6.3</td><td>a</td><td>2.4</td><td>0.1</td><td>a</td><td>5.22</td><td>0.31</td><td>a</td><td>348</td><td>10</td><td>a</td></tr><tr><td>MDHAR</td><td>25.3</td><td>2.2</td><td>a</td><td>191.8</td><td>14.0</td><td>b</td><td>2.7</td><td>0.3</td><td>b</td><td>6.06</td><td>0.79</td><td>b</td><td>326</td><td>37</td><td>a</td></tr></tbody></table>

Measurements were carried out on greenhouse-grown plants as described in the Materials and methods. Data show means from 10 plants per line with standard error (SE). A comparison of means was carried out using a Kruskal Wallis test with correction (Dunn). Different letters indicate significant differences (5% significance level).

<table><caption>**TABLE 2** | Fruit phenotypes of transgenic lines and WT.</caption><thead><tr><th></th><th>Days to orange</th><th>SE</th><th>p</th><th>Individual fruit weight (g)</th><th>SE</th><th>p</th><th>% Dry weight (orange)</th><th>SE</th><th>p</th></tr></thead><tbody><tr><td>WT</td><td>40.3</td><td>2.1</td><td>a</td><td>6.3</td><td>1.2</td><td>a</td><td>7.1</td><td>0.1</td><td>a</td></tr><tr><td>AO</td><td>41.9</td><td>1.5</td><td>b</td><td>6.1</td><td>0.8</td><td>a</td><td>7.0</td><td>0.1</td><td>a</td></tr><tr><td>GLD</td><td>42.7</td><td>1.4</td><td>b</td><td>6.0</td><td>1.0</td><td>a</td><td>6.8</td><td>0.1</td><td>a</td></tr><tr><td>MDHAR</td><td>40.7</td><td>1.8</td><td>a</td><td>3.8</td><td>0.9</td><td>b</td><td>7.0</td><td>0.1</td><td>a</td></tr></tbody></table>

Measurements were carried out on greenhouse-grown plants as described in the Materials and methods. Data show means from at least 6 fruits per plant of 10 plants per line with standard error (SE). A comparison of means was carried out using a Kruskal Wallis test with correction (Dunn). Different letters indicate significant differences (5% significance level).

orange fruit, no significant differences are seen. The polyphenol that varied most in the transgenic lines was chlorogenic acid, which is also the most prevalent polyphenol, and its stereoisomer cis-chlorogenic acid: as for dehydroascorbate, levels were significantly different in AO and MDHAR lines, AO lines containing less chlorogenic acid than wild-type whereas the MDHAR lines contained more (Table 5 and Figure 2). Levels of these polyphenols in the GLD RNAi fruit were not significantly different to wild-type. The metabolites detected and quantified by ¹H-NMR are shown in Figure 2 (and Supplementary Table S2). Metabolites that differentiate the lines included the amino acids alanine and tyrosine which were highest in MDHAR and lowest in AO lines (significant for alanine) and the hexose sugars glucose and fructose which showed a tendency in contrast to the aforementioned amino acids to be highest in AO lines.

## The Fruit Proteome Is Inversely Affected by MDHAR and AO Expression

The differentially expressed protein spots led to the identification of 34 proteins by mass spectrometry (Table 6). Differences between the lines are presented as a heatmap of the log2 of the ratio with the WT (Figure 3; Table 6, raw data in Supplementary Table S3). The line most affected in terms of its fruit proteome is MDHAR, for which up-regulated proteins included a 60S ribosomal protein, a class I heat-shock protein and a cysteine protease inhibitor. Down-regulated proteins included a heat-shock protein, (hsp70 homolog), lactoylglutathione lyase and an amidase hydantoinase/carbamoylase family protein.

The proteome of the GLD line did not show significant differences to the wild-type, despite visually resembling the ascorbate oxidase profile. A number of proteins showed opposite regulation in AO and MDHAR lines, these included chaperone DnaK (MDHAR down-regulated, AO up-regulated), a probable carbonic anhydrase, a protein disulphide isomerase and a metacaspase 7, a protein with cysteine endopeptidase activity and a positive mediator of cell death which is AO down-regulated and MDHAR up-regulated. One protein was up-regulated in the AO line which is an F1F0-ATPase inhibitor protein. The data for fruit 20 days after anthesis is shown in Supplementary Figure S1: only two proteins showed statistically significant differences in spot volumes between wild-type and the MDHAR line.

## Weighted-Gene Correlation Network Analysis Reveals Clusters of Genes Associated with Both Metabolites and Proteins

The normalized transcriptome data for the four lines was subjected to a weighted-gene correlation network analysis to identify clusters of correlated genes and compare them with different phenotypes (metabolite and protein levels). We identified 38 modules of between 26 and 851 genes (Supplementary Figure S2). The number of genes in each module was too small to give significant differences following GO enrichment analysis but revealed co-regulation of metabolites or proteins with genes in modules or groups of modules. For example, the amino acids were positively correlated with the gene modules "dark red," "blue," "green-yellow" and "yellow." The full data set with module membership is found in Supplementary Table S4.

## Pairwise Comparison of Transcriptome Data Reveals the Same Set of Differentially Regulated Genes in the Comparison AO-MDHAR and AO-GLD

Each of the three transgenic lines was compared with wild-type and then the lines were compared two by two. For the pairwise comparisons of the transgenic lines, the hypothesis that the gene expression in the wild-type was the average of the gene expression in each transgenic line was tested using the maanova model (see Materials and methods). Volcano plots for each pairwise comparison showing the log2 of the fold-change in gene
<table><caption>**TABLE 3** | Fruit physical characteristics of transgenic lines and WT at harvest.</caption><thead><tr><th></th><th>Fruit weight (g)</th><th>SE</th><th>p</th><th>Diameter (mm)</th><th>SE</th><th>p</th><th>Compression (kPa)</th><th>SE</th><th>p</th><th>Max force-skin rupture (N)</th><th>SE</th><th>p</th><th>Color (a*)</th><th>SE</th><th>p</th><th>Postharvest weight loss (%)</th><th>SE</th><th>p</th></tr></thead><tbody><tr><td>WT</td><td>9.3</td><td>0.4</td><td>a</td><td>25.7</td><td>0.4</td><td>a</td><td>49.5</td><td>2.5</td><td>a</td><td>5.7</td><td>0.2</td><td>a</td><td>10.3</td><td>1.0</td><td>a</td><td>0.043</td><td>0.001</td><td>a</td></tr><tr><td>AO</td><td>8.3</td><td>0.4</td><td>a</td><td>24.4</td><td>0.5</td><td>ab</td><td>43.1</td><td>3.0</td><td>a</td><td>4.6</td><td>0.2</td><td>b</td><td>13.8</td><td>1.6</td><td>a</td><td>0.049</td><td>0.003</td><td>ab</td></tr><tr><td>GLD</td><td>7.8</td><td>0.4</td><td>a</td><td>23.6</td><td>0.5</td><td>b</td><td>49.4</td><td>4.0</td><td>a</td><td>5.4</td><td>0.3</td><td>ab</td><td>13.1</td><td>1.7</td><td>a</td><td>0.046</td><td>0.001</td><td>a</td></tr><tr><td>MDHAR</td><td>5.3</td><td>0.3</td><td>b</td><td>21.1</td><td>0.4</td><td>c</td><td>49.9</td><td>2.6</td><td>a</td><td>4.6</td><td>0.3</td><td>b</td><td>15.7</td><td>1.7</td><td>a</td><td>0.063</td><td>0.002</td><td>bc</td></tr></tbody></table>

Fruit physical characteristics were measured on red-ripe fruit of at least 20 fruit per line as described in the Materials and methods. Measurements show means with standard error (SE). A comparison of means was carried out using a Kruskal Wallis test with correction (Dunn). Different letters indicate significant differences (5% significance level). *a* is the name of the color scale.

<table><caption>**TABLE 4** | Ascorbate and dehydroascorbate levels in orange fruit pericarp of transgenic lines and WT.</caption><thead><tr><th></th><th>Reduced AsA (mg/100gfwt)</th><th>SE</th><th>p</th><th>DHA (mg/100gfwt)</th><th>SE</th><th>p</th></tr></thead><tbody><tr><td>WT</td><td>15.46</td><td>1.30</td><td>a</td><td>2.71</td><td>1.39</td><td>ab</td></tr><tr><td>AO</td><td>18.99</td><td>1.97</td><td>ab</td><td>4.73</td><td>1.32</td><td>a</td></tr><tr><td>GLD</td><td>19.19</td><td>0.65</td><td>ab</td><td>2.83</td><td>0.97</td><td>ab</td></tr><tr><td>MDHAR</td><td>21.72</td><td>1.05</td><td>b</td><td>0</td><td>0.68</td><td>b</td></tr></tbody></table>

Ascorbate and dehydroascorbate were assayed in the pericarp tissue (3 biological replicates of 30 fruits per pool) of the lines by a spectrophotometric method. Measurements show means with standard error (SE). A comparison of means was carried out using a Kruskal Wallis test with correction (Dunn). Different letters indicate significant differences (5% significance level).

expression between the lines and the associated *p*-value (−log10 transformed) are shown in **Figure 4**. The AO-MDHAR and AO- GLD comparisons generated the highest number of differentially expressed genes (comparisons 4; 290 genes and comparison 6; 330 genes). The comparisons WT = (transgenic X + transgenic Y)/2 are comparisons 5, 8, and 9. If the hypothesis is true then the *p*-values become non-significant and the −log10 of the *p*- values gives low values which is the case for the comparisons 5 [WT = (AO + MDHAR)/2] and 8 [WT = (AO + GLD)/2] but not the case for the comparison 9 [WT = (GLD + MDHAR)/2]. The whole dataset including means, standard error, *p* and *F*- values and the fold-change with *p*-values for each comparison is available in **Supplementary Table S5**.

## The Differentially Expressed Genes Are Enriched for Functions Involving Translation and Ribosome-Biogenesis

The AO-MDHAR and AO-GLD comparisons generated the highest number of differentially expressed genes and the expression in the WT is always intermediate. Comparison of these two gene sets showed 182 genes common to each comparison, defining a core set of genes for further analysis. Independently of the comparison and the direction of change, differentially expressed genes were part of one “molecular function” category: structural constituent of ribosome (GO:00037) and two “biological process” categories: ribosome

biogenesis and translation (GO:00422 and GO:00064). These genes were up-regulated in the AO line and down-regulated in the MDHAR and GLD RNAi lines. The core set of 182 genes were analyzed using the *Solanaceae* Genome Network GO-enrichment tool and this analysis gave further categories associated with translation, ribosomes, unfolded protein binding, GTPase regulator activity and nucleotide binding. The categories “structural constituent of ribosome,” “small GTPase regulator activity” and “oxidoreductase activity acting on NADH or NADPH” were more than ten-fold enriched in the cluster compared to the genome (Table 7).

## The Subset of 182 Genes Defines a Regulatory Network with Several Hub Genes Including Heat-Shock Transcription Factors, Chaperones, and Ribosomal Proteins

The 182 core gene set was used in network analysis with the metabolites and the protein dataset. A topological overlap matrix using the WGCNA script was produced and an edge file for network visualization using Cytoscape. The network visualized with a degree-sorted circle layout, is shown in **Figure 5** with the top eleven hub genes. The most highly connected genes included two or three ribosomal proteins, an observation which is compatible with those above, and two heat-shock proteins (including DnaJ) and a heat-shock transcription factor.

## The Transcriptomes Show a Signature Common to Stress Responses in Arabidopsis, Particularly Mitochondrial Stress but No Correlation with Hydrogen Peroxide Signaling

Our tomato transcriptome data was compared using ROSMETER, a bioinformatic tool which unites transcriptome data from different stress experiments or mutants in Arabidopsis linked to specific ROS production in specific organelles (see Materials and methods for more information). Its use is based on the hypothesis that the chemical identity and sub-cellular localization of ROS leave a specific imprint on the transcriptomic response (Rosenwasser et al., 2013). We used the RNAi AO and RNAi MDHAR datasets to compare with the Arabidopsis transcriptomes. Our first observation is that independently of the
<table><caption>TABLE 5 | Polyphenols in orange fruit pericarp of transgenic lines and WT.</caption><thead><tr><th></th><th>Chlorogenic acid (µg/gdwt)</th><th>SE</th><th>p</th><th>Caffeic acid glucoside (arbitrary units)</th><th>SE</th><th>p</th><th>Cis-chlorogenic acid (arbitrary units)</th><th>SE</th><th>p</th><th>Quercetin derivative (arbitrary units)</th><th>SE</th><th>p</th><th>Rutin (µg/gdwt)</th><th>SE</th><th>p</th></tr></thead><tbody><tr><td>WT</td><td>3245.2</td><td>31.2</td><td>ab</td><td>3.5</td><td>0.0</td><td>a</td><td>1.3</td><td>0.0</td><td>ab</td><td>0.4</td><td>0.00</td><td>a</td><td>435.5</td><td>2.9</td><td>a</td></tr><tr><td>AO</td><td>2537.5</td><td>81.5</td><td>a</td><td>3.0</td><td>0.2</td><td>a</td><td>1.0</td><td>0.1</td><td>a</td><td>0.4</td><td>0.02</td><td>a</td><td>397.1</td><td>8.1</td><td>a</td></tr><tr><td>GLD</td><td>3029.0</td><td>298.5</td><td>ab</td><td>3.2</td><td>0.5</td><td>a</td><td>1.1</td><td>0.1</td><td>ab</td><td>0.4</td><td>0.01</td><td>a</td><td>406.0</td><td>21.4</td><td>a</td></tr><tr><td>MDHAR</td><td>3554.4</td><td>68.3</td><td>b</td><td>3.7</td><td>0.2</td><td>a</td><td>1.5</td><td>0.1</td><td>b</td><td>0.3</td><td>0.01</td><td>a</td><td>446.7</td><td>21.3</td><td>a</td></tr></tbody></table>

Polyphenolic compounds were assayed in the pericarp tissue (3 biological replicates of 30 fruits per pool) of the lines by HPLC-DAD. When the purified compound was not available for exact quantification, arbitrary units corresponding to peak surface were used. Measurements show means with standard error (SE). A comparison of means was carried out using a Kruskal Wallis test with correction (Dunn). Different letters indicate significant differences (5% significance level).

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

stress and/or organelle, the transcriptome responses are inverted when these two lines are compared (Figure 6). The particular conditions to be highlighted include mitochondrial stress: in this case, the application of rotenone gives a strong negative correlation with the RNAi MDHAR and a positive correlation with RNAi AO 3 h after application. Negative correlations are also found for RNAi MDHAR with Arabidopsis catalase mutants after 3 and 8 h of high light, whereas the correlations are positive for RNAi AO. Positive correlations between RNAi MDHAR exist with 30 min of methyl viologen and generally the knockout ascorbate peroxidase mutant subjected to high light. Again for RNAi AO, in contrast to RNAi MDHAR, a negative correlation is found for 30 min of methyl viologen treatment.

DISCUSSION

This study in fruit tissue has revealed some interesting correlations which might provide clues for further research on ascorbate function but some caveats should be pointed out: here we worked on ripening fruit, a non-photosynthetic tissue, where small differences in ascorbate content are observed between the various transgenic lines. The information provided will help decipher the relationships between ascorbate and fruit physiology, especially in fruit subject to biotic and abiotic stresses. However, it remains to be seen if the correlations are more general and an initial question would be to see if they are transferable to other tissues. As a 35S promoter was used it is
<table><caption>TABLE 6 | Proteins in orange fruit pericarp of the lines as identified by mass spectrometry.</caption><thead><tr><th>Spot ID</th><th>ITAG2.4</th><th>Annotated function</th><th>Significant difference in lines</th><th>Fold change of transcript AO vs. MDHAR</th><th>p-value (FDR)</th><th>Correlation transcript-protein</th></tr></thead><tbody><tr><td>CM148</td><td>Solyc01g099190.2.1</td><td>Lipoxygenase</td><td></td><td>No data</td><td>No data</td><td></td></tr><tr><td>MP25</td><td>Solyc01g103450.2.1</td><td>Chaperone DnaK</td><td>MDHAR down, AO up</td><td></td><td></td><td></td></tr><tr><td>MP26</td><td>Solyc01g104920.2.1</td><td>26S protease regulatory subunit 8 homolog</td><td></td><td>-1.84</td><td>0.002</td><td></td></tr><tr><td>MP14</td><td>Solyc01g111300.2.1</td><td>Cold shock protein-1</td><td></td><td></td><td></td><td></td></tr><tr><td>MP31</td><td>Solyc02g070510.2.1</td><td>Proteasome subunit alpha type</td><td></td><td>No data</td><td>No data</td><td></td></tr><tr><td>MB16</td><td>Solyc02g080630.2.1</td><td>Lactoylglutathione lyase</td><td>MDHAR down</td><td></td><td></td><td></td></tr><tr><td>MP18</td><td>Solyc02g085790.2.1</td><td>T-complex protein 1 subunit zeta</td><td></td><td>5.82</td><td>0.005</td><td></td></tr><tr><td>MP38</td><td>Solyc02g091100.2.1</td><td>Oxalyl-CoA decarboxylase</td><td></td><td></td><td></td><td></td></tr><tr><td>MP42</td><td>Solyc02g092670.1.1</td><td>Subtilisin-like protease</td><td></td><td>No data</td><td>No data</td><td></td></tr><tr><td>JC65</td><td>Solyc03g097270.2.1</td><td>Cysteine proteinase inhibitor</td><td>MDHAR up</td><td>-1.40</td><td>0.017</td><td>Yes</td></tr><tr><td>MP24</td><td>Solyc04g045340.2.1</td><td>Phosphoglucomutase</td><td></td><td></td><td></td><td></td></tr><tr><td>MP29</td><td>Solyc04g045340.2.1</td><td>Phosphoglucomutase</td><td></td><td></td><td></td><td></td></tr><tr><td>CM57</td><td>Solyc05g005490.2.1</td><td>Carbonic anhydrase</td><td>MDHAR up, AO down</td><td>1.98</td><td>0.008</td><td>No</td></tr><tr><td>MP17</td><td>Solyc05g012580.1.1</td><td>Unknown Protein</td><td></td><td></td><td></td><td></td></tr><tr><td>MP11</td><td>Solyc05g054580.2.1/ Solyc09g091520.1.1</td><td>60S acidic ribosomal protein P0</td><td>MDHAR up</td><td>-1.34</td><td>0.023</td><td>Yes</td></tr><tr><td>JC26</td><td>Solyc06g005160.2.1</td><td>Ascorbate peroxidase</td><td></td><td>1.25</td><td>0.036</td><td></td></tr><tr><td>MP13</td><td>Solyc06g005940.2.1</td><td>Protein disulfide isomerase</td><td></td><td></td><td></td><td></td></tr><tr><td>MP33</td><td>Solyc06g005940.2.1</td><td>Protein disulfide isomerase</td><td></td><td></td><td></td><td></td></tr><tr><td>JC60</td><td>Solyc06g059740.2.1</td><td>Alcohol dehydrogenase 2</td><td></td><td>1.49</td><td>0.004</td><td></td></tr><tr><td>TP17</td><td>Solyc06g059740.2.1</td><td>Alcohol dehydrogenase 2</td><td></td><td>1.49</td><td>0.004</td><td></td></tr><tr><td>MP20</td><td>Solyc06g060290.2.1</td><td>Protein disulfide isomerase</td><td>MDHAR up, AO down</td><td>No data</td><td>No data</td><td></td></tr><tr><td>MP36</td><td>Solyc06g065270.2.1</td><td>Adenylate kinase</td><td></td><td></td><td></td><td></td></tr><tr><td>TP55</td><td>Solyc06g076520.1.1</td><td>Class I heat shock protein</td><td>MDHAR up</td><td></td><td></td><td></td></tr><tr><td>CM94</td><td>Solyc07g005820.2.1</td><td>Heat shock protein 70</td><td>MDHAR down</td><td></td><td></td><td></td></tr><tr><td>MP23</td><td>Solyc07g066600.2.1</td><td>Phosphoglycerate kinase</td><td></td><td>-1.36</td><td>0.033</td><td></td></tr><tr><td>MP22</td><td>Solyc08g014130.2.1</td><td>2-isopropylmalate synthase 1</td><td></td><td>2.05</td><td>0.000</td><td></td></tr><tr><td>JC79</td><td>Solyc08g014130.2.1</td><td>2-isopropylmalate synthase 1</td><td></td><td>2.05</td><td>0.000</td><td></td></tr><tr><td>MP21</td><td>Solyc08g062660.2.1</td><td>Ran GTPase binding protein</td><td></td><td>No data</td><td>No data</td><td></td></tr><tr><td>MP41</td><td>Solyc09g082060.2.1</td><td>Cysteine synthase</td><td></td><td></td><td></td><td></td></tr><tr><td>MP37</td><td>Solyc09g083410.2.1</td><td>Amidase hydantoinase/carbamoyl family protein</td><td>MDHAR down</td><td>No data</td><td>No data</td><td></td></tr><tr><td>MP19</td><td>Solyc09g090330.2.1</td><td>Harpin binding protein 1</td><td></td><td>No data</td><td>No data</td><td></td></tr><tr><td>MP12</td><td>Solyc09g098150.2.1</td><td>Metacaspase 7</td><td>MDHAR up, AO down</td><td>2.01</td><td>0.003</td><td>No</td></tr><tr><td>MP16</td><td>Solyc09g098150.2.1</td><td>Metacaspase 7</td><td></td><td>2.01</td><td>0.003</td><td></td></tr><tr><td>MP32</td><td>Solyc09g098150.2.1</td><td>Metacaspase 7</td><td></td><td>2.01</td><td>0.003</td><td></td></tr><tr><td>MP34</td><td>Solyc10g081240.1.1</td><td>Protein grpE</td><td></td><td></td><td></td><td></td></tr><tr><td>MP28</td><td>Solyc10g086580.1.1</td><td>Ribulose-1 5-bisphosphate carboxylase/oxygenase activase 1</td><td></td><td></td><td></td><td></td></tr><tr><td>MP27</td><td>Solyc11g006970.1.1</td><td>Unknown protein DS12 from 2D-PAGE of leaf, chloroplastic</td><td></td><td>No data</td><td>No data</td><td></td></tr><tr><td>MP30</td><td>Solyc11g068510.1.1</td><td>F1F0-ATPase inhibitor protein</td><td>AO up</td><td>No data</td><td>No data</td><td></td></tr><tr><td>MP35</td><td>Solyc11g069000.1.1</td><td>T-complex protein 1 subunit beta</td><td></td><td>-1.46</td><td>0.035</td><td></td></tr><tr><td>MP15</td><td>Solyc12g056230.1.1</td><td>Glutathione peroxidase</td><td></td><td>No data</td><td>No data</td><td></td></tr></tbody></table>

The predicted Solyc (https://solgenomics.net/organism/Solanum_lycopersicum/genome; version 2.4 of the International Tomato Annotation Group, ITAG) is shown with the annotated function and significant difference between lines where applicable (5% significance level). The fold-change and p-value of the corresponding transcript, when present on the microarray and for p-values of < 0.05, is also shown for the AO-MDHAR lines and whether this fold-change is correlated with the difference seen in protein levels.

also possible that the effects we see in fruit are the result of silencing the genes in another part of the plant (for an example of where this is true see our recent results by Truffault et al.,

2017). The use of available fruit-specific promoters (Fernandez et al., 2009) and indeed these genes are expressed in different plant organs and fruit tissues (Supplementary Figures S3, S4)
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Heatmap showing protein levels in orange pericarp fruit (3 biological replicates of 30 fruits) of the transgenic lines and wild-type. All proteins were separated by two-dimensional gel electrophoresis and identified by mass spectrometry. The log2 of the protein ratio with the wild-type for each transgenic line is presented. The scale goes from green (protein decreased compared to wild-type) to purple (protein increased). A comparison of the means was carried out using a Kruskal Wallis test with correction (Dunn). Different letters indicate significant differences (5% significance level). The proteins are identified by their Spot reference ID, the correspondence for which can be found in **Table 6**.</figcaption></figure>

would aid in specifically attributing the effects we observed to the fruit. A final caveat is that the correlations we observed are based on microarray data: the 7600 gene tomato microarray is considerably enriched in genes expressed in ripening fruit (Alba et al., 2005; Fei et al., 2011) and therefore well adapted to the purpose of the study. However, global transcriptome analysis using RNA-seq would additionally allow the detection of weakly expressed genes and thus strengthen our conclusions.

## A Transcriptional Signature Is Observed and Is Correlated to Silencing of the Ascorbate Metabolizing Enzymes

The transcriptomes revealed a switch linked to the activity of ascorbate oxidase on one side of the switch and the enzymes MDHAR/GLD acting in the other direction, always with wild-type gene expression in the middle. There is a tendency for

genes to be up-regulated when ascorbate oxidase is down-regulated. The transcriptome inversion is also seen in datasets from Arabidopsis provided by the ROSMETER tool. The switch seen between AO and MDHAR RNAi lines is capable of distinguishing the conditions found in some of the ROSMETER transcriptome data, although there are limits to this analysis because correlations can be based on a small number of genes (a minimum of 45 according to the original paper). The ROSMETER analysis does not reveal links with hydrogen peroxide metabolism because correlations between our lines and the public data are inversed when comparing ascorbate peroxidase and catalase mutants, both of which in theory should contain more hydrogen peroxide as a result of decreased detoxification. The strong response to rotenone is interesting because rotenone is an inhibitor of mitochondrial complex I, which is also where the galactonolactone dehydrogenase protein is located (Schimmeyer et al., 2016). Also in the paper by Rosenwasser et al. (2013), rotenone-induced stress is also positively correlated with developmental senescence, a process which has parallels with fruit ripening.

The data seem to indicate the existence of a signal affecting transcription closely linked to GLD activity and the redox state of the ascorbate pool. The nature of this signal is unknown but a candidate is monodehydroascorbate, as it is the molecule that links ascorbate oxidase and monodehydroascorbate reductase activity, which could act as a primary signal given (i) the short half-life of the molecule, a prerequisite for a role in signaling and (ii) the proton transfer that occurs when the radical is reduced back to ascorbate (Bielski et al., 1981) which gives potential for protein modification. Most oxidation reactions with ascorbate are indeed single electron reactions (Heber et al., 1996) generating the radical monodehydroascorbate, this radical is considered to be common and a marker for stress. Another possibility is that redox changes in mitochondria act as a retrograde signal (Rhoads and Subbaiah, 2007); this could be a possibility given the involvement of GLD. Interestingly a dehydroascorbate reductase (DHAR) gene (Solyc05g054760) is one of the 182 differentially regulated genes and is down-regulated in GLD and MDHAR lines and up-regulated in the AO line. Compensatory effects, either by enzymes such as DHAR, or other isoforms of the enzymes under study, cannot be ruled out.

## Ascorbate Has Links with Translation and Protein Folding

The transcriptional signature and network analysis suggest a role for ascorbate in control of translation and protein folding, at least in the non-photosynthetic tissue studied. Some evidence from the literature has established links between ascorbate and these cellular mechanisms: a good example is the posttranscriptional negative feedback control that ascorbate exerts over its synthesis via the GDP-I-galactose phosphorylase enzyme (Laing et al., 2015) which appears to affect ribosome density on the 5'UTR of the mRNA. Ribosomes and translation are also regulated by stress and redox processes (Halbeisen and Gerber, 2009; Khandal et al., 2009; Mustroph et al., 2009; Gerashchenko et al., 2012; Gismondi et al., 2014; Moore et al., 2016). Protein folding is
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Pairwise transcriptome comparisons as illustrated using volcano plots (log2 of fold change plotted against -log10 of the *p*-value for the differential expression). The expression of the 7633 genes from the microarray analysis of RNA extracted from the orange fruit pericarp (3 biological replicates of 30 fruit) of the different lines was compared with technical duplicates (dyes swapped). The comparisons were carried out first between the transgenic lines (GalLDH = GLD) and wild-type (AO-WT, comparison 1; GalLDH-WT, comparison 2; MDHAR-WT, comparison 3) before comparing the transgenic lines between themselves (AO-MDHAR, comparison 4, AO-GalLDH, comparison 6, and GalLDH-MDHAR, comparison 7). Afterward the hypothesis was tested that the WT expression was the average of the expression in two given transgenic lines (see Materials and methods; WT = average AO and MD, comparison 5, WT = average AO and GalLDH, comparison 8, WT = average GalLDH and MDHAR, comparison 9). The black line is drawn to represent a threshold *p*-value of 0.001, the red points represent differentially expressed genes.</figcaption></figure>

also affected by environmental challenges: the unfolded protein response results from a period of stress and leads to either up-regulation of proteins required for folding or degradation or in extreme cases cell death if prolonged (Ruberti et al., 2015; Wan and Jiang, 2016). Functional hypotheses such as these would need to be tested. For example, purification of polysomes from the available transgenic lines could lead to information on the RNAs undergoing translation and/or associated proteins and ribosome subunits. Considering the fruit context, this experiment would be especially informative in ripening fruits subjected to high temperature stress or in ripe fruits subjected to stressful

postharvest storage conditions that may result in protein mis-folding such as high CO₂ (Rothan et al., 1997) or chilling temperatures (Stevens et al., 2008).

## Hub Gene Identification Points to Novel Regulatory Roles Linked to Cellular Energy Status and Signaling

Heat-shock proteins also figure amongst the hub genes identified: three heat-shock proteins are identified: SolycHsfB2a (Solyc03g026020), a Class II small heat-shock protein
<table><caption>**TABLE 7** | Gene ontology enrichment on the 182 set of common genes using the tools available at Boyce Thompson Institute, Cornell University (http://bioinfo.bti.cornell.edu/cgi-bin/MetGenMAP/home.cgi).</caption><thead><tr><th>Gene ontology term</th><th>Number of genes</th><th>Cluster frequency (%)</th><th>Genome frequency (%)</th><th>Corrected p-value</th></tr></thead><tbody><tr><td>Unfolded protein binding</td><td>3</td><td>1.6</td><td>0.2</td><td>0.038</td></tr><tr><td>Structural molecule activity</td><td>11</td><td>6.0</td><td>0.7</td><td>0.000</td></tr><tr><td>Structural constituent of ribosome</td><td>10</td><td>5.5</td><td>0.5</td><td>0.000</td></tr><tr><td>Small molecule binding</td><td>16</td><td>8.8</td><td>3.9</td><td>0.028</td></tr><tr><td>Small GTPase regulator activity</td><td>3</td><td>1.6</td><td>0.1</td><td>0.032</td></tr><tr><td>Ribonucleotide binding</td><td>11</td><td>6.0</td><td>2.0</td><td>0.025</td></tr><tr><td>Purine ribonucleotide binding</td><td>11</td><td>6.0</td><td>2.0</td><td>0.023</td></tr><tr><td>Purine ribonucleoside triphosphate binding</td><td>11</td><td>6.0</td><td>2.0</td><td>0.025</td></tr><tr><td>Purine nucleotide binding</td><td>11</td><td>6.0</td><td>2.0</td><td>0.021</td></tr><tr><td>Protein binding</td><td>44</td><td>24.2</td><td>11.9</td><td>0.000</td></tr><tr><td>Oxidoreductase activity acting on NADH or NADPH</td><td>4</td><td>2.2</td><td>0.2</td><td>0.003</td></tr><tr><td>Oxidoreductase activity</td><td>16</td><td>8.8</td><td>2.8</td><td>0.000</td></tr><tr><td>Organic cyclic compound binding</td><td>16</td><td>8.8</td><td>3.3</td><td>0.013</td></tr><tr><td>Nucleotide binding</td><td>16</td><td>8.8</td><td>3.3</td><td>0.015</td></tr><tr><td>Nucleoside phosphate binding</td><td>16</td><td>8.8</td><td>3.3</td><td>0.012</td></tr><tr><td>Catalytic activity</td><td>53</td><td>29.1</td><td>17.9</td><td>0.000</td></tr><tr><td>Binding</td><td>79</td><td>43.4</td><td>25.1</td><td>0.000</td></tr><tr><td>ATP binding</td><td>9</td><td>4.9</td><td>1.4</td><td>0.020</td></tr><tr><td>Adenyl ribonucleotide binding</td><td>9</td><td>4.9</td><td>1.4</td><td>0.019</td></tr><tr><td>Adenyl nucleotide binding</td><td>9</td><td>4.9</td><td>1.4</td><td>0.020</td></tr></tbody></table>

The table shows the cluster frequency (the 182 genes) compared to the genome frequency for each category with the FDR-corrected *p*-value. The three categories (structural constituent of ribosomes, small GTPase regulator activity and oxidoreductase activity acting on NADH or NADPH) are where there is the biggest difference between the genome frequency and the cluster frequency.

Le-HSP17.6 (Solyc08g062340) and a chaperone protein DnaJ (Solyc03g123560). Heat-shock proteins are known to be involved in numerous biological phenomena including stress responses, circadian clock regulation, flowering time, and hypocotyl elongation (Charng et al., 2007; Mittler et al., 2012; Kolmos et al., 2014). The hub gene protein DnaJ (hsp40) stimulates the ATPase activity of DnaK (an hsp70 homolog which is interestingly differentially accumulated in the AO and MDHAR proteomes). These proteins bind to unfolded polypeptides to prevent their aggregation, for example under conditions of heat-shock. The DnaK protein is of interest because it has been detected in the proteome of the Arabidopsis mutant for GDP-L-galactose phosphorylase (vtc2), a major control point for ascorbate synthesis, as being up-regulated after 3 and 5 days of high light compared to the wild-type (Giacomelli et al., 2006). This parallels the situation in the AO line (up-regulation) and is opposite to the MDHAR line phenotype (down-regulation) meaning that ascorbate deficiency mirrors the RNAi AO phenotype. Another hub protein, Arabinogalactan, is involved in many biological processes including germination, cell extension, fertilization, binding, and release of calcium and is linked to the auxin response (Ellis et al., 2010). Only one of the hub genes has already been characterized from tomato (Table 7) which is the tft8 protein, a member of the 14-3-3 protein family involved in defense responses (Roberts and Bowles, 1999).

The major phenotypic differences seen in the transgenic lines are related to growth: in previous studies RNAi AO lines tend

to show improved yield stability (Garchery et al., 2013) whereas the opposite is true for GLD and MDHAR RNAi lines (Alhagdow et al., 2007; Truffault et al., 2016), a phenomenon which is correlated to the changes in the fruit transcriptomes and could be related to fruit growth. Explanations for the phenotypes might be found in the identities of the hub genes which are broadly involved in signaling translation and protein (re)folding. We also noted the presence of many ATP-dependent proteins. The yield phenotypes could therefore reflect differences in cellular energy levels.

## Differences at the Protein and Metabolite Level Do Not Follow Transcriptional Differences and Show a Second Level of Control

As many studies have observed, RNA and protein abundance are often weakly correlated (Halbeisen and Gerber, 2009; Petricka et al., 2012) and our data corroborate this: the quantity of the differentially accumulated proteins is generally not linked to their transcriptional levels (Table 6). For the observable proteins, the GLD proteome is not significantly affected: the GLD proteome actually is remarkably visually similar to the AO proteome which is surprising as the transcriptome data showed contrasted differential expression of genes for these two lines. The metabolic profile is similar to the proteome profile in that several metabolites show opposite changes compared to wild-type between AO and MDHAR RNAi
<figure><img src="image_7.png" /></figure>

<table><thead><tr><th>ITAG2.4</th><th>Degree of association</th><th>Function</th><th>Fold change of transcript AO vs MDHAR</th><th>p value (FDR)</th></tr></thead><tbody><tr><td>Solyc06g084230.2.1</td><td>34</td><td>40S ribosomal protein S24</td><td>2.30</td><td>5.68E-05</td></tr><tr><td>Solyc03g031700.2.1</td><td>18</td><td>ABC transporter FeS assembly protein SufB</td><td>2.79</td><td>2.50E-05</td></tr><tr><td>Solyc12g096220.1.1</td><td>16</td><td>60S ribosomal protein L7</td><td>4.31</td><td>5.68E-05</td></tr><tr><td>Solyc12g010860.1.1</td><td>13</td><td>14-3-3 protein, tft8</td><td>2.19</td><td>6.37E-05</td></tr><tr><td>Solyc02g092790.2.1</td><td>10</td><td>Arabinogalactan</td><td>4.86</td><td>0.00019</td></tr><tr><td>Solyc03g026020.2.1</td><td>9</td><td>Heat stress transcription factor, SolycHsfB2a</td><td>2.45</td><td>0</td></tr><tr><td>Solyc08g091450.2.1</td><td>8</td><td>DNA polymerase delta subunit 4 family</td><td>1.85</td><td>2.50E-05</td></tr><tr><td>Solyc08g062340.2.1</td><td>8</td><td>Class II small heat shock protein Le-HSP17.6</td><td>3.54</td><td>0.00014</td></tr><tr><td>Solyc07g053440.2.1</td><td>7</td><td>Os06g0538200 protein, conserved hypothetical protein</td><td>2.22</td><td>0.00110</td></tr><tr><td>Solyc05g007210.2.1</td><td>7</td><td>Unknown protein? Ribosomal protein L41 interpro domain</td><td>7.34</td><td>2.50E-05</td></tr><tr><td>Solyc03g123560.2.1</td><td>7</td><td>Chaperone protein dnaJ</td><td>2.03</td><td>0.00012</td></tr></tbody></table>

**FIGURE 5** | Network analysis (Cytoscape) showing the most highly connected genes based on their expression correlations. A topological overlap matrix was generated using the WGCNA script for the 182 core genes, the 22 metabolites and the 40 proteins. The edge file produced was imported into Cytoscape for network visualization. The network was visualized with a degree-sorted circle layout. The node sizes and node label sizes are proportional to the degree of connectivity of each node as fixed by Cytoscape. The color scale goes from blue (least connected nodes) to orange (most connected nodes). The hub genes are described in the table which gives the most highly connected genes following visualization of the topology overlap matrix, generated by the WGCNA "R" package, with Cytoscape. The ITAG2.4 identifier (Solyc) with the annotated function is shown with the degree of association as defined by Cytoscape. The transcriptome data (fold change between AO and MDHAR RNAi lines with the FDR corrected *p*-value on the transcript levels) for each gene are also shown.

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

**FIGURE 6** | ROSMETER analysis of the WT vs. AO and WT vs. MDHAR transcriptome data from tomato fruit pericarp. Tomato genes were converted to the Arabidopsis homolog and the fold-changes and *p*-values were compared to the Arabidopsis transcriptomes previously described (Rosenwasser et al., 2013) which covered different stress conditions and/or mutant backgrounds in Arabidopsis as clusters as described in the text to the figure. The heatmap generated shows green for general transcriptome down-regulation (opposite fold-change) and red for up-regulation (similar fold-change) in the tomato lines compared to Arabidopsis.

lines, this includes dehydroascorbate, chlorogenic acid, and alanine. The changes seen in dehydroascorbate are counter-intuitive when compared to what is expected based on enzyme activity. We, with other independent groups, have previously reported increases in ascorbate levels on reduction of MDHAR activity (a peroxisomal and cytosolic isoform; Haroldsen et al., 2011; Gest et al., 2013a) which implies feedback regulation. Links have been shown between polyphenols such as chlorogenic acid or anthocyanins and ascorbate (Page et al., 2012); also oxidation of chlorogenic acid is correlated with oxidation of ascorbate to dehydroascorbate (Takahama, 1998) and ascorbate can regenerate chlorogenic acid from its corresponding radicals (Yamasaki and Grace, 1998). Our results show less chlorogenic acid in lines with lower ascorbate oxidase. We should point out that the genotype used in this study (West Virginia 106) is particularly rich in chlorogenic acid, the wild-type contains 3.25 mg/g dwt of chlorogenic acid whereas reported concentrations range from 0.24 to 3.69 mg/100 g fresh weight (between ~0.024 and 0.37 mg/g dwt; Slimestad and Verheul, 2005; Peng et al., 2008). The differences in alanine are also of interest as this metabolite is produced directly from pyruvate and therefore has close links with glycolysis, the oxidative pentose phosphate pathway and
the TCA cycle and cellular energy production. Alanine is a product of anaerobic metabolism resulting from high rates of amino acid interconversion (Rocha et al., 2010b). Unlike lactate and ethanol, alanine accumulation does not have detrimental side effects but maintains the glycolytic flux (Rocha et al., 2010a).

In conclusion, this study in a non-photosynthetic tissue reveals that the activity of enzymes controlling the oxidation state of the ascorbate pool are correlated to a few metabolic changes and some changes in protein quantity, particularly those involved in protein folding. The relatively low number of metabolic changes compared to gene expression probably reflects buffering of the system. In addition to this post-transcriptional regulation, a transcriptional response is seen which is also linked to the activity of the final biosynthesis enzyme, galactonolactone dehydrogenase. Plants silenced for this enzyme show a transcriptional response similar to that of MDHAR silenced plants and opposite to ascorbate oxidase silenced plants implying signaling from the ascorbate pool. Network analysis reveals hub genes involved in fundamental cellular mechanisms including cellular signaling, translation and protein synthesis as well as the heat-shock response.

## AUTHOR CONTRIBUTIONS

RS: supervised the experiments, carried out the bioinformatic analysis of the data and wrote the article with the contribution of all the authors; MF and J-PB: carried out the proteome study; PB, CD, MM, AM, and J-LP: carried out the metabolic analyses; VG and JP: carried out the microarrays and transcriptome data normalization; ClD: analyzed the ROSMETER results; CG: carried out the greenhouse measurements, sampling and provided technical assistance to RS; BG and DP: carried out the post-harvest study; MC, HG, and VT: helped with the writing; CR: conceived the project and obtained the funding. All authors approved the final version of the manuscript.

## FUNDING

The authors are grateful to the INRA AgroBI programme for financing the study, "VTCFruit" project.

## ACKNOWLEDGMENTS

We acknowledge valuable help from Noé Fernandez from SGN for converting the TOM2 probes to Solyc identifiants and Shilo Rosenwasser for help with application of the ROSMETER tool to tomato data. We would also like to thank Manu Botton for help with the plant culture and pre-harvest phenotyping, Joël Chadoeuf for advice on statistical models in Rmaanova, Camille Croset for polyphenol analysis, Marine Palau for contributing to the proteome analysis and Christopher Sauvage for advice on the R scripts. The metabolome analyses were performed at Bordeaux Metabolome facility-MetaboHUB. We thank Daniel Jacob for developing and maintaining the MeRy-B database. A

version of the abstract was published at http://qualityfruit.inp-toulouse.fr/_attachments/preliminary-program-article/Book%2520of%2520Abstracts%2520COST%2520Porto2016%2520%252814out16%2529.pdf?download=true, as part of the COST Action FA1106 "Quality Fruit."

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00137/full#supplementary-material

**Supplementary Figure S1** | Heatmap showing protein levels in pericarp of fruit 20 days after anthesis (3 pools of 30 fruits) of the transgenic lines and wild-type. All proteins were separated by two-dimensional gel electrophoresis and identified by mass spectrometry. The log2 of the protein ratio with the wild-type for each transgenic line is presented. The scale goes from green (protein decreased compared to wild-type) to purple (protein increased). A comparison of the means was carried out using a Kruskal Wallis test with correction (Dunn). Different letters indicate significant differences (5% significance level). SGN codes are found at https://solgenomics.net/.

**Supplementary Figure S2** | WGCNA module identification. A graphical representation of the data from **Supplementary Table S4**: module eigengenes are labeled by color and the module membership scale from red (high adjacency, positive correlation) to green (low adjacency, negative correlation) shown for each metabolite and protein. The 7633 genes were used to generate the topology overlap matrix. SGN codes are found at https://solgenomics.net/.

**Supplementary Figure S3** | Gene expression profiles for AO, Solyc04g054690; GLD, Solyc10g079470 and MDHAR, Solyc09g009390 as obtained from the TomExpress tool (http://tomexpress.toulouse.inra.fr).

**Supplementary Figure S4** | Gene expression profiles for AO, Solyc04g054690; GLD, Solyc10g079470 and MDHAR, Solyc09g009390 as obtained from the Solgenomics expression atlas tool (http://tea.solgenomics.net).

**Supplementary Table S1** | Ascorbate and dehydroascorbate levels in pericarp fruit 20 days after anthesis of transgenic lines and WT. Ascorbate and dehydroascorbate were assayed in the pericarp tissue (3 biological replicates of 30 fruits per pool) of the lines by a spectrophotometric method. Measurements show means with standard error (SE). A comparison of means was carried out using a Kruskal Wallis test with correction (Dunn). Different letters indicate significant differences (5% significance level).

**Supplementary Table S2** | Fruit metabolite data (orange fruit pericarp). Quantitative data for the following metabolites: galactose, mannose, glutamate, glutamine, alanine, aspartate, tyrosine, fructose, glucose, sucrose, ascorbate, dehydroascorbate, GABA, citramalate, citrate, malate, chlorogenic acid, caffeic acid glucoside, cis-chlorogenic acid, quercetin derivative, and rutin as measured by specific methods (¹H-NMR for polar compounds, HPLC-DAD for the major polyphenols, spectrophotometrically for ascorbate and dehydroascorbate) for each of the three pools for the wild type, AO, GLD and MDHAR RNAi lines.

**Supplementary Table S3** | Proteins. Raw data for quantifiable proteins (see Materials and methods) from orange pericarp with spot identifier code for each of the three pools for the wild type, AO, GLD, and MDHAR RNAi lines.

**Supplementary Table S4** | WGCNA. Solyc identifier, gene function, GO terms, module color, module membership (MM or eigengene-based connectivity) and associated *p*-value for the 38 modules and 7633 genes as generated by the WGCNA R package.

**Supplementary Table S5** | Transcriptome raw data. For each of the 7633 genes mean expression (normalized) in each of the three transgenic lines and wild type with standard error (SE). *F* values, *p*-values, false discovery rate corrected *p*-values and fold change (FC) are also given for each comparison as follows: 1: AO vs. WT, 2: GLD vs. WT, 3: MDHAR vs. WT, 4: AO vs. MDHAR, 5: AO + MDHAR/(2 × WT) = myModel, 6: AO vs. GLD, 7: GLD vs. MDHAR, 8: AO + GLD/(2 × WT) = myModel2, 9: GLD + MDHAR/(2 × WT) = myModel3.
REFERENCES

Alba, R., Payton, P., Fei, Z., McQuinn, R., Debbie, P., Martin, G. B., et al. (2005). Transcriptome and selected metabolite analyses reveal multiple points of ethylene control during tomato fruit development. *Plant Cell* 17, 2954–2965. doi: 10.1105/tpc.105.036053

Alhagdow, M., Mounet, F., Gilbert, L., Nunes-Nesi, A., Garcia, V., Just, D., et al. (2007). Silencing of the mitochondrial ascorbate synthesizing enzyme L-galactono-1,4-lactone dehydrogenase affects plant and fruit development in tomato. *Plant Physiol.* 145, 1408–1422. doi: 10.1104/pp.107.106500

Banhegyi, G., Csala, M., Szarka, A., Varsanyi, M., Benedetti, A., and Mandl, J. (2003). Role of ascorbate in oxidative protein folding. *BioFACTORS* 17, 37–46. doi: 10.1002/biof.5520170105

Bassel, G. W., Gaudinier, A., Brady, S. M., Hennig, L., Rhee, S. Y., and De Smet, I. (2012). Systems analysis of plant functional, transcriptional, physical interaction, and metabolic networks. *Plant Cell* 24, 3859–3875. doi: 10.1105/tpc.112.100776

Bielski, B. H. J., Allen, A. O., and Schwarz, H. A. (1981). Mechanisms of disproportionation of ascorbate radicals. *J. Am. Chem. Soc.* 103, 3516–3518. doi: 10.1021/ja00402a042

Blaschke, K., Ebata, K. T., Karimi, M. M., Zepeda-Martinez, J. A., Goyal, P., Mahapatra, S., et al. (2013). Vitamin C induces Tet-dependent DNA demethylation and a blastocyst-like state in ES cells. *Nature* 500, 222–226. doi: 10.1038/nature12362

Bors, W., and Buettner, G. R. (1997). "The vitamin C radical and its reactions," in *Vitamin C in Health and Disease*, eds L. Packer and J. Fuchs (New York, NY: Marcel Dekker), 75–94.

Boyle, E. I., Weng, S., Gollub, J., Jin, H., Botstein, D., Cherry, J. M., et al. (2004). GO::TermFinder-open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. *Bioinformatics* 20, 3710–3715. doi: 10.1093/bioinformatics/bth456

Carcamo, J. M., Pedraza, A., Borquez-Ojeda, O., Zhang, B., Sanchez, R., and Golde, D. W. (2004). Vitamin C is a kinase inhibitor: dehydroascorbic acid inhibits IkappaBα kinase β. *Mol. Cell. Biol.* 24, 6645–6652. doi: 10.1128/MCB.24.15.6645-6652.2004

Charng, Y. Y., Liu, H. C., Liu, N. Y., Chi, W. T., Wang, C. N., Chang, S. H., et al. (2007). A heat-inducible transcription factor, HsfA2, is required for extension of acquired thermotolerance in Arabidopsis. *Plant Physiol.* 143, 251–262. doi: 10.1104/pp.106.091322

Cisternas, P., Silva-Alvarez, C., Martinez, F., Fernandez, E., Ferrada, L., Oyarce, K., et al. (2014). The oxidized form of vitamin C, dehydroascorbic acid, regulates neuronal energy metabolism. *J. Neurochem.* 129, 663–671. doi: 10.1111/jnc.12663

Cline, M. S., Smoot, M., Cerami, E., Kuchinsky, A., Landys, N., Workman, C., et al. (2007). Integration of biological networks and gene expression data using Cytoscape. *Nat. Protoc.* 2, 2366–2382. doi: 10.1038/nprot.2007.324

Davey, M. W., Auwerkerken, A., and Keulemans, J. (2007). Relationship of apple vitamin C and antioxidant contents to harvest date and postharvest pathogen infection. *J. Sci. Food Agric.* 87, 802–813. doi: 10.1002/jsfa.2777

De Tullio, M. C. (2012). Beyond the antioxidant: the double life of vitamin C. *Subcell. Biochem.* 56, 49–65. doi: 10.1007/978-94-007-2199-9_4

De Tullio, M. C., Ciraci, S., Liso, R., and Arrigoni, O. (2007). Ascorbic acid oxidase is dynamically regulated by light and oxygen. A tool for oxygen management in plants? *J. Plant Physiol.* 164, 39–46. doi: 10.1016/j.jplph.2005.09.016

Ellis, M., Egelund, J., Schultz, C. J., and Bacic, A. (2010). Arabinogalactan-proteins: key regulators at the cell surface? *Plant Physiol.* 153, 403–419. doi: 10.1104/pp.110.156000

Faurobert, M., Mihr, C., Bertin, N., Pawlowski, T., Negroni, L., Sommerer, N., et al. (2007). Major proteome variations associated with cherry tomato pericarp development and ripening. *Plant Physiol.* 143, 1327–1346. doi: 10.1104/pp.106.092817

Fei, Z., Joung, J. G., Tang, X., Zheng, Y., Huang, M., Lee, J. M., et al. (2011). Tomato functional genomics database: a comprehensive resource and analysis package for tomato functional genomics. *Nucleic Acids Res.* 39, D1156–D1163. doi: 10.1093/nar/gkq991

Fernandez, A. I., Viron, N., Alhagdow, M., Karimi, M., Jones, M., Amsellem, Z., et al. (2009). Flexible tools for gene expression and silencing in tomato. *Plant Physiol.* 151, 1729–1740. doi: 10.1104/pp.109.147546

Fiorani, M., De Sanctis, R., Scarlatti, F., Vallorani, L., De Bellis, R., Serafini, G., et al. (2000). Dehydroascorbic acid irreversibly inhibits hexokinase activity. *Mol. Cell. Biochem.* 209, 145–153. doi: 10.1023/A:1007168032289

Garchery, C., Gest, N., Do, P. T., Alhagdow, M., Baldet, P., Menard, G., et al. (2013). A diminution in ascorbate oxidase activity affects carbon allocation and improves yield in tomato under water deficit. *Plant Cell Environ.* 36, 159–175. doi: 10.1111/j.1365-3040.2012.02564.x

Gerashchenko, M. V., Lobanov, A. V., and Gladyshev, V. N. (2012). Genome-wide ribosome profiling reveals complex translational regulation in response to oxidative stress. *Proc. Natl. Acad. Sci. U.S.A.* 109, 17394–17399. doi: 10.1073/pnas.1120799109

Gest, N., Garchery, C., Gautier, H., Jimenez, A., and Stevens, R. (2013a). Light-dependent regulation of ascorbate in tomato by a monodehydroascorbate reductase localized in peroxisomes and the cytosol. *Plant Biotechnol. J.* 11, 344–354. doi: 10.1111/pbi.12020

Gest, N., Gautier, H., and Stevens, R. (2013b). Ascorbate as seen through plant evolution: the rise of a successful molecule? *J. Exp. Bot.* 64, 33–53. doi: 10.1093/jxb/ers297

Giacomelli, L., Rudella, A., and van Wijk, K. J. (2006). High light response of the thylakoid proteome in arabidopsis wild type and the ascorbate-deficient mutant vtc2-2. A comparative proteomics study. *Plant Physiol.* 141, 685–701. doi: 10.1104/pp.106.080150

Gismondi, A., Caldarola, S., Lisi, G., Juli, G., Chellini, L., Iadevaia, V., et al. (2014). Ribosomal stress activates eEF2K-eEF2 pathway causing translation elongation inhibition and recruitment of terminal oligopyrimidine (TOP) mRNAs on polysomes. *Nucleic Acids Res.* 42, 12668–12680. doi: 10.1093/nar/gku996

Halbeisen, R. E., and Gerber, A. P. (2009). Stress-dependent coordination of transcriptome and translatome in yeast. *PLoS Biol.* 7:e1000105. doi: 10.1371/journal.pbio.1000105

Haroldsen, V. M., Chi-Ham, C. L., Kulkarni, S., Lorence, A., and Bennett, A. B. (2011). Constitutively expressed DHAR and MDHAR influence fruit, but not foliar ascorbate levels in tomato. *Plant Physiol. Biochem.* 49, 1244–1249. doi: 10.1016/j.plaphy.2011.08.003

Heber, U., Miyake, C., Mano, J., Ohno, C., and Asada, K. (1996). Monodehydroascorbate radical detected by electron paramagnetic resonance spectrometry is a sensitive probe of oxidative stress in intact leaves. *Plant Cell Physiol.* 37, 1066–1072. doi: 10.1093/oxfordjournals.pcp.a029055

Khandal, D., Samol, I., Buhr, F., Pollmann, S., Schmidt, H., Clemens, S., et al. (2009). Singlet oxygen-dependent translational control in the tigrinad.12 mutant of barley. *Proc. Natl. Acad. Sci. U.S.A.* 106, 13112–13117. doi: 10.1073/pnas.0903522106

Klee, H. J., and Giovannoni, J. J. (2011). Genetics and control of tomato fruit ripening and quality attributes. *Annu. Rev. Genet.* 45, 41–59. doi: 10.1146/annurev-genet-110410-132507

Kolmos, E., Chow, B. Y., Pruneda-Paz, J. L., and Kay, S. A. (2014). HsfB2b-mediated repression of PRR7 directs abiotic stress responses of the circadian clock. *Proc. Natl. Acad. Sci. U.S.A.* 111, 16172–16177. doi: 10.1073/pnas.1418483111

Krouk, G., Lingeman, J., Colon, A. M., Coruzzi, G., and Shasha, D. (2013). Gene regulatory networks in plants: learning causality from time and perturbation. *Genome Biol.* 14:123. doi: 10.1186/gb-2013-14-6-123

Laing, W. A., Martinez-Sanchez, M., Wright, M. A., Bulley, S. M., Brewster, D., Dare, A. P., et al. (2015). An upstream open reading frame is essential for feedback regulation of ascorbate biosynthesis in Arabidopsis. *Plant Cell* 27, 772–786. doi: 10.1105/tpc.114.133777

Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. *BMC Bioinformatics* 9:559. doi: 10.1186/1471-2105-9-559

Leferink, N. G., van Duijn, E., Barendregt, A., Heck, A. J., and van Berkel, W. J. (2009). Galactonolactone dehydrogenase requires a redox-sensitive thiol for optimal production of vitamin C. *Plant Physiol.* 150, 596–605. doi: 10.1104/pp.109.136929

Malacrida, C., Valle, E., and Boggio, S. (2006). Postharvest chilling induces oxidative stress response in the dwarf tomato cultivar Micro-Tom. *Physiol. Plant.* 127, 10–18. doi: 10.1111/j.1399-3054.2005.00636.x
Mittler, R., Finka, A., and Goloubinoff, P. (2012). How do plants feel the heat? *Trends Biochem. Sci.* 37, 118–125. doi: 10.1016/j.tibs.2011.11.007

Moing, A., Maucourt, M., Renaud, C., Gaudillere, M., Brouquisse, R., Leboutellier, B., et al. (2004). Quantitative metabolic profiling through one-dimensional 1H-NMR analyses: application to plant genetics and functional genomics. *Funct. Plant Biol.* 31, 889–902. doi: 10.1071/FP04066

Moore, M., Gossmann, N., and Dietz, K. J. (2016). Redox regulation of cytosolic translation in plants. *Trends Plant Sci.* 21, 388–397. doi: 10.1016/j.tplants.2015.11.004

Mounet, F., Lemaire-Chamley, M., Maucourt, M., Cabasson, C., Giraudel, J. L., Deborde, C., et al. (2007). Quantitative metabolic profiles of tomato flesh and seeds during fruit development: complementary analysis with ANN and PCA. *Metabolomics* 3, 273–288. doi: 10.1007/s11306-007-0059-1

Mounet, F., Moing, A., Garcia, V., Petit, J., Maucourt, M., Deborde, C., et al. (2009). Gene and metabolite regulatory network analysis of early developing fruit tissues highlights new candidate genes for the control of tomato fruit composition and development. *Plant Physiol.* 149, 1505–1528. doi: 10.1104/pp.108.133967

Mustroph, A., Zanetti, M. E., Jang, C. J., Holtan, H. E., Repetti, P. P., Galbraith, D. W., et al. (2009). Profiling translatomes of discrete cell populations resolves altered cellular priorities during hypoxia in Arabidopsis. *Proc. Natl. Acad. Sci. U.S.A.* 106, 18843–18848. doi: 10.1073/pnas.0906131106

Page, M., Sultana, N., Paszkiewicz, K., Florance, H., and Smirnoff, N. (2012). The influence of ascorbate on anthocyanin accumulation during high light acclimation in *Arabidopsis thaliana*: further evidence for redox control of anthocyanin synthesis. *Plant Cell Environ.* 35, 388–404. doi: 10.1111/j.1365-3040.2011.02369.x

Peng, Y., Zhang, Y., and Ye, J. (2008). Determination of phenolic compounds and ascorbic acid in different fractions of tomato by capillary electrophoresis with electrochemical detection. *J. Agric. Food Chem.* 56, 1838–1844. doi: 10.1021/jf0727544

Petricka, J. J., Schauer, M. A., Megraw, M., Breakfield, N. W., Thompson, J. W., Georgiev, S., et al. (2012). The protein expression landscape of the Arabidopsis root. *Proc. Natl. Acad. Sci. U.S.A.* 109, 6811–6818. doi: 10.1073/pnas.1202546109

Puskas, F., Gergely, P. Jr., Banki, K., and Perl, A. (2000). Stimulation of the pentose phosphate pathway and glutathione levels by dehydroascorbate, the oxidized form of vitamin C. *FASEB J.* 14, 1352–1361. doi: 10.1096/fasebj.14.10.1352

R Development-Core-Team (2015). *R: A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing.

Rhoads, D. M., and Subbaiah, C. C. (2007). Mitochondrial retrograde regulation in plants. *Mitochondrion* 7, 177–194. doi: 10.1016/j.mito.2007.01.002

Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. *Nucleic Acids Res.* 43, e47. doi: 10.1093/nar/gkv007

Roberts, M. R., and Bowles, D. J. (1999). Fusicoccin, 14-3-3 proteins, and defense responses in tomato plants. *Plant Physiol.* 119, 1243–1250. doi: 10.1104/pp.119.4.1243

Rocha, M., Licausi, F., Araujo, W. L., Nunes-Nesi, A., Sodek, L., Fernie, A. R., et al. (2010a). Glycolysis and the tricarboxylic acid cycle are linked by alanine aminotransferase during hypoxia induced by waterlogging of *Lotus japonicus*. *Plant Physiol.* 152, 1501–1513. doi: 10.1104/pp.109.150045

Rocha, M., Sodek, L., Licausi, F., Hameed, M. W., Dornelas, M. C., and van Dongen, J. T. (2010b). Analysis of alanine aminotransferase in various organs of soybean (*Glycine max*) and in dependence of different nitrogen fertilisers during hypoxic stress. *Amino Acids* 39, 1043–1053. doi: 10.1007/s00726-010-0596-1

Rosenwasser, S., Fluhr, R., Joshi, J. R., Leviatan, N., Sela, N., Hetzroni, A., et al. (2013). ROSMETER: a bioinformatic tool for the identification of transcriptomic imprints related to reactive oxygen species type and origin provides new insights into stress responses. *Plant Physiol.* 163, 1071–1083. doi: 10.1104/pp.113.218206

Rothan, C., Duret, S., Chevalier, C., and Raymond, P. (1997). Suppression of ripening-associated gene expression in tomato fruits subjected to a high CO₂ concentration. *Plant Physiol.* 114, 255–263. doi: 10.1104/pp.114.1.255

Ruberti, C., Kim, S. J., Stefano, G., and Brandizzi, F. (2015). Unfolded protein response in plants: one master, many questions. *Curr. Opin. Plant Biol.* 27, 59–66. doi: 10.1016/j.pbi.2015.05.016

Schimmeyer, J., Bock, R., and Meyer, E. H. (2016). L-Galactono-1,4-lactone dehydrogenase is an assembly factor of the membrane arm of mitochondrial complex I in Arabidopsis. *Plant Mol. Biol.* 90, 117–126. doi: 10.1007/s11103-015-0400-4

Slimestad, R., and Verheul, M. J. (2005). Seasonal variations in the level of plant constituents in greenhouse production of cherry tomatoes. *J. Agric. Food Chem.* 53, 3114–3119. doi: 10.1021/jf047864e

Stevens, R., Buret, M., Garchery, C., Carretero, Y., and Causse, M. (2006). Technique for rapid, small-scale analysis of vitamin C levels in fruit and application to a tomato mutant collection. *J. Agric. Food Chem.* 54, 6159–6165. doi: 10.1021/jf061241e

Stevens, R., Page, D., Gouble, B., Garchery, C., Zamir, D., and Causse, M. (2008). Tomato fruit ascorbic acid content is linked with monodehydroascorbate reductase activity and tolerance to chilling stress. *Plant Cell Environ.* 31, 1086–1096. doi: 10.1111/j.1365-3040.2008.01824.x

Szarka, A., and Lorincz, T. (2014). The role of ascorbate in protein folding. *Protoplasma* 251, 489–497. doi: 10.1007/s00709-013-0560-5

Takahama, U. (1998). Ascorbic acid-dependent regulation of redox levels of chlorogenic acid and its isomers in the apoplast of leaves of *Nicotiana tabacum* L. *Plant Cell Physiol.* 39, 681–689. doi: 10.1093/oxfordjournals.pcp.a029422

Truffault, V., Gest, N., Garchery, C., Florian, A., Fernie, A. R., Gautier, H., et al. (2016). Reduction of MDHAR activity in cherry tomato suppresses growth and yield and MDHAR activity is correlated with sugar levels under high light. *Plant Cell Environ.* 39, 1279–1292. doi: 10.1111/pce.12663

Truffault, V., Riqueau, G., Garchery, C., Gautier, H., and Stevens, R.G. (2017). Is monodehydroascorbate reductase activity in leaf tissue critical for the maintenance of yield in tomato? *J. Plant Physiol.* 222, 1–8. doi: 10.1016/j.jplph.2017.12.012

Wan, S., and Jiang, L. (2016). Endoplasmic Reticulum (ER) stress and the Unfolded Protein Response (UPR) in plants. *Protoplasma* 253, 753–764. doi: 10.1007/s00709-015-0842-1

Xu, J., Pascual, L., Aurand, R., Bouchet, J. P., Valot, B., Zivy, M., et al. (2013). An extensive proteome map of tomato (*Solanum lycopersicum*) fruit pericarp. *Proteomics* 13, 3059–3063. doi: 10.1002/pmic.201200438

Yamasaki, H., and Grace, S. C. (1998). EPR detection of phytophenoxyl radicals stabilized by zinc ions: evidence for the redox coupling of plant phenolics with ascorbate in the H₂O₂-peroxidase system. *FEBS Lett.* 422, 377–380. doi: 10.1016/S0014-5793(98)00048-9

Yang, Y. H., Dudoit, S., Luu, P., Lin, D. M., Peng, V., Ngai, J., et al. (2002). Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. *Nucleic Acids Res.* 30:e15. doi: 10.1093/nar/30.4.e15

Young, M. D., Wakefield, M. J., Smyth, G. K., and Oshlack, A. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. *Genome Biol.* 11:R14. doi: 10.1186/gb-2010-11-2-r14

Yun, J., Mullarky, E., Lu, C., Bosch, K. N., Kavalier, A., Rivera, K., et al. (2015). Vitamin C selectively kills KRAS and BRAF mutant colorectal cancer cells by targeting GAPDH. *Science* 350, 1391–1396. doi: 10.1126/science.aaa5004

Zito, E., Hansen, H. G., Yeo, G. S., Fujii, J., and Ron, D. (2012). Endoplasmic reticulum thiol oxidase deficiency leads to ascorbic acid depletion and noncanonical scurvy in mice. *Mol. Cell* 48, 39–51. doi: 10.1016/j.molcel.2012.08.010

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Stevens, Baldet, Bouchet, Causse, Deborde, Deschodt, Faurobert, Garchery, Garcia, Gautier, Gouble, Maucourt, Moing, Page, Petit, Poëssel, Truffault and Rothan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.