Next Article in Journal
Assessing Silicon-Mediated Growth Performances in Contrasting Rice Cultivars under Salt Stress
Next Article in Special Issue
UV-B Irradiation Effect on Microalgae Performance in the Remediation of Effluent Derived from the Cigarette Butt Cleaning Process
Previous Article in Journal
Grapevine Pinot Gris Virus Is Present in Different Non-Vitis Hosts
Previous Article in Special Issue
Growth and Biosynthesis of Phenolic Compounds of Canola (Brassica napus L.) to Different Ultraviolet (UV)-B Wavelengths in a Plant Factory with Artificial Light
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Foliar and Root Comparative Metabolomics and Phenolic Profiling of Micro-Tom Tomato (Solanum lycopersicum L.) Plants Associated with a Gene Expression Analysis in Response to Short Daily UV Treatments

by
Alessia Mannucci
1,2,†,
Marco Santin
1,*,†,
Lucas Vanhaelewyn
2,3,4,†,
Maria Calogera Sciampagna
1,
Maria Begoña Miras-Moreno
5,
Leilei Zhang
5,
Luigi Lucini
5,
Mike Frank Quartacci
1,
Dominique Van Der Straeten
2,
Antonella Castagna
1 and
Annamaria Ranieri
1
1
Department of Agriculture, Food and Environment, University of Pisa, 56124 Pisa, Italy
2
Laboratory of Functional Plant Biology, Department of Biology, Faculty of Sciences, Ghent University, 9000 Ghent, Belgium
3
Department of Agricultural Economics, Ghent University, Coupure Links 653, 9000 Ghent, Belgium
4
Deroose Plants NV, Weststraat 129A, 9940 Evergem, Belgium
5
Department for Sustainable Food Process, Università Cattolica del Sacro Cuore, 29122 Piacenza, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Plants 2022, 11(14), 1829; https://doi.org/10.3390/plants11141829
Submission received: 26 May 2022 / Revised: 6 July 2022 / Accepted: 9 July 2022 / Published: 12 July 2022

Abstract

:
Tomato (Solanum lycopersicum L.) is globally recognised as a high-value crop both for commercial profit and nutritional benefits. In contrast to the extensive data regarding the changes in the metabolism of tomato fruit exposed to UV radiation, less is known about the foliar and root metabolome. Using an untargeted metabolomic approach through UHPLC-ESI-QTOF-MS analysis, we detected thousands of metabolites in the leaves (3000) and roots (2800) of Micro-Tom tomato plants exposed to 11 days of short daily UV radiation, applied only on the aboveground organs. Multivariate statistical analysis, such as OPLS-DA and volcano, were performed to allow a better understanding of the modifications caused by the treatment. Based on the unexpected modulation to the secondary metabolism, especially the phenylpropanoid pathway, of which compounds were down and up accumulated respectively in leaves and roots of treated plants, a phenolic profiling was carried out for both organs. The phenolic profile was associated with a gene expression analysis to check the transcription trend of genes involved in the UVR8 signalling pathway and the early steps of the phenolic biosynthesis. The retention of the modifications at metabolic and phenolic levels was also investigated 3 days after the UV treatment, showing a prolonged effect on the modulation once the UV treatment had ceased.

1. Introduction

UV-A and part of UV-B radiation (UV, ultraviolet radiation; 290–400 nm) reach the Earth’s surface, triggering the photomorphogenic responses required for the acclimation of plants to this environmental stimulus [1,2,3]. Several studies have elucidated the role of UV radiation as a plant growth regulator; however, if exposed at a high dose, UV can negatively affect many physiological processes, causing the onset of oxidative stress that can no longer be counteracted by the boost of the antioxidant production, which is also triggered by UV light [4]. However, plants have developed several ways to cope with UV radiation, such as modifications of the metabolic pathways, along with an increase in the phenylpropanoid biosynthesis, which is one of the main topics under this environmental condition. Several other pathways are known to be responsive to UV radiation, including amino acid biosynthesis [5], as they also are precursors of specific secondary metabolites needed in the UV-B responses.
UV RESISTANCE LOCUS 8 (UVR8) has been identified as the specific photoreceptor able to perceive UV-B and trigger downstream signalling [6]. This dimeric protein is localized in both the nucleus and cytoplasm under normal light conditions; however, after UV-B exposure, it splits into monomers and binds to several light-regulating proteins such as CONSTITUTIVELY PHOTOMORPHOGENIC 1 (COP1) [7,8], leading to its accumulation in the nucleus compartment. The sequestration of COP1 inhibits the proteolysis of target proteins [1,8], among which is the transcription factor HY5 (ELONGATED HYPOCOTYL 5), enabling the transcription of downstream genes such as phenolic-related genes [4].
It has been widely demonstrated that metabolites such as phenolics and flavonoids are fundamental as protectants against oxidative damages caused by reactive oxygen species (ROSs) produced during stress responses, and they also act as UV-B-absorbing compounds in the epidermal cell layers, reducing the transmission of this wavelength to the mesophyll, and hence protecting the photosystems [2,9,10]. UV-B induces the transcription of genes encoding the enzymes involved in the first steps of the phenylpropanoid pathway, such as PHENYLALANINE AMMONIA LYASE (PAL), CHALCONE SYNTHASE (CHS), and FLAVONOL SYNTHASE (FLS) [11,12]. Therefore, the activation of the UVR8 signalling pathway enables plant morphologic modifications and the accumulation of UV screening compounds that enhance the plant’s antioxidant defence [6,13]. Notably, roots also express the UVR8 photoreceptor and several regulators required for UV-B signal transduction with some root specific gene induction [6,14,15]. Moreover, this wavelength is known to influence root development and direction of growth [16,17]; however, less is known regarding the biochemical changes occurring in this organ when only the shoot is exposed to UV-B.
This study aimed at having a collective view of the foliar and root metabolic modifications, with a focus on the phenolic profiling, during the acclimation of Micro-Tom tomato plants to short daily UV treatments for 11 days exclusively on the aboveground organs. Subsequently, the retention of the modifications at metabolic and phenolic levels was investigated at 3 days after the removal of UV treatment (recovery phase). Moreover, the quantification of phenolics with an “-omic” approach was associated with a gene expression analysis to check the transcription of genes involved in the UVR8 signalling and the early steps of the phenolic biosynthesis. Literature in which a metabolomic approach was applied to investigate UV and particularly UV-B effects on plants is still scarce, and, to the best of our knowledge, they were mainly addressed to fruit and foliar metabolome [18,19,20,21]. However, we know that roots are also equipped with the whole signalling molecules needed for UV responses [15,22], and that several signalling factors can generally move from shoot-to-root and vice versa [23]. Although light piping mechanism has been observed regarding far red light [24], similar evidence has not been reported for UV-B radiation yet. While the potential of UV radiation in increasing the plant antioxidant defence is well documented, here, we provide new insights regarding how UV exposure in the aerial part of the plant can affect below-ground tissues and how phenolics are differently modulated in the shoots, even after several days of recovery.

2. Results and Discussion

2.1. Effects of Short Daily UV Treatment of the Phyllosphere on the Leaf and Root Metabolism of Tomato Plants

The UHPLCESI/QTOF-MS analysis, coupled with a comprehensive database for primary and secondary metabolites identification (PlantCyc), allowed for the detection of about 3000 metabolites in leaves and 2800 metabolites in roots. The whole list of metabolites detected is provided as supplementary material, together with abundances and composite mass spectra (Table S2). Given the broad chemical diversity in the metabolome, several multivariate statistical analyses were performed to allow a better understanding of the modifications caused by the UV treatment.
The output of the unsupervised fold-change-based hierarchical clustering (HCA) performed on the foliar metabolomic profile (Figure 1) showed a separation between UV and CTR samples, with the main clustering between UV-rec and the other groups. However, in the case of roots, the output of the HCA (Figure 1) showed the main clustering between CTR-rec and the other samples, with UV groups being clustered together irrespective of the time point considered. This preliminary analysis highlighted a different behavior of the two organs’ metabolism.
Then, an OPLS-DA supervised model (Figure 2a,b for leaves, Figure 2c,d for roots) was used to enhance the interpretation of the variables that maximize the differences among sample groups for each investigated organ.
Figure 2a,c clearly empathized the difference existing between CTR and UV samples (independently from the time point considered) in both leaves and roots. Consequently, Variable Importance for Prediction (VIP) analysis was carried out to show which metabolites weighed the most in the loading plot of these supervised models. Table 1 lists the metabolic markers in leaves, with 122 compounds in total, related to Figure 2a, with a VIP score higher than 1.1.
Univariate statistical analysis was performed to identify the VIP markers displaying a significant fold-change (FC) between UV and CTR foliar samples. About half of the compounds listed displayed a significant p value. Some compounds such as 2-oxoadipate, oxalosuccinate, scopolamine, 9-[6(S),9-diamino-5,6,7,8,9-pentadeoxy-β-D-ribo-nonafuranosyl]-9H-purin-6-amine, 3-hydroxy-9-apo-δ-caroten-9-one, and tricoumaroyl spermidine showed a remarkable and statistically significant enhanced FC in respect to CTR samples. On the contrary, N-P-tosyl-L-phenylalanyl chloromethyl ketone, 1-[18-hydroxyoeoyl]-2-[18-hydroxy-linoleoyl]-sn-glycerol, (22R,23R)-28-homocastasterone, and cyclo-dopa 5-O-glucoside reported a marked and significantly reduced FC in respect to CTR samples.
Roots VIP metabolic markers related to the OPLS-DA plot are shown in Figure 2c and listed in Table 2.
As for leaves, univariate statistical analysis led to identifying the compounds with a significant FC in respect to CTR root samples. Among 169 markers, 62 compounds displayed a significant FC. In particular, 10-deacetyl-2-debenzoylbaccatin III, a pentose, N-methylanthranilate, L-α-(methylenecyclopropyl)-glycine, 5,10-methylenetetrahydropteroyl mono-L-glutamate, all-trans-hexaprenyl diphosphate, 4-(β-D-glucosyloxy) benzoate, L-nicotianamine, dTMP, and N-(4-aminobenzoyl)-L-glutamate showed a remarkably enhanced or reduced (cycloheptadienyl/sinapate) FC in respect to CTR samples.
Finally, among the hundreds of VIP markers resulting from the OPLSDA plots, only eight markers were found to be in common between leaves and roots, and these were compounds related to the secondary metabolism such as O-sinapoylglucarolactone, cyanidin 3-O-β-D-p-coumaroylglucoside/cyanidin 3-(p-coumaroyl) glucoside/pelargonidin 3-O-β-D-caffeoylglucoside, scopolamine, aminoacids such as cyclogutamate and N-methylanthranilate, and 5-[[4-methoxy-3-(phenylmethoxy)phenyl]methyl]-2,4-pyrimidinediamine.
The OPLS-DA in Figure 2b,d shows the effects of UV and time factors in leaves and roots, respectively. A good separation among the different groups is visible in both organs; however, in leaves, the treatment played the main role for the distribution along the first vector, with UV samples placed in the negative half of the plot and CTR samples in the positive half, while the sampling time did not have a great effect. Surprisingly roots exhibited not only the effect of the UV treatment, with UV samples and CTR samples in the positive and negative half of the score plot considering the second vector, but also a time effect, with a good separation of the samples along the first vector according to this factor.

2.2. Outputs of the Volcano Analysis on the Metabolome

Differential compounds, with individual p values and fold-change resulting from the volcano plot, are fully provided as supplementary material (Table S3). To visualise how the metabolic pathways were modulated by short daily UV treatments, the differential compounds resulting from the volcano analysis were interpreted by the Omics Viewer Dashboard of the PlantCyc Pathway Tool software. Figure 3 shows several metabolic pathways in leaves and roots, with the main focus on biosynthesis and degradation processes.
Leaves of UV-11d plants displayed an increased fold-change in the amino acid and hormone biosynthesis (+28 and +23 log2FC, respectively) in respect to the relative CTR, while the fatty acid/lipid and secondary metabolites showed a decrease in respect to CTR leaves on the same day (−24 and −100 log2FC, respectively, Figure 3). Especially the latter biosynthetic process was still reduced after the recovery, while fatty acid/lipid biosynthesis increased in respect to CTR (+47 log2FC). Fatty acids, as components of the cellular membrane, also play a role in stress signalling responses leading to the alteration in the fluidity of the membrane. As an example, the release of unsaturated fatty acids is able to activate the defensive genes in tomato leaves under herbivore wounding [25]. Yang and collaborators [26], studying the effects of high UV-B irradiation (120.8 μW cm−2 for 5 h followed by incubation with dark) on Clematis terniflora, found an accumulation of linolenic and linoleic acids. In our experiment, the fluctuation in the FA/lipid synthesis reported in UV-treated leaves and the increase in their biosynthesis after 3 days of recovery suggest a remodelling of the membrane in respect to untreated plants and once the treatment was removed. Specifically, the lipids assigned to the biosynthetic process comprised, among others, two galactolipids, one phospholipid, three sterols, and three polyunsaturated fatty acids. Among these latter, (12Z,15Z)-9,10-epoxyoctadeca-12,15-dienoate and (9Z,12Z)-15,16-dihydroxyoctadeca-9,12-dienoate, which are involved in the production of hydroxy fatty acids from α-linolenic acid, were oppositely affected by the UV-B exposure, being increased or decreased, respectively. Since some hydroxylated fatty acids act as antifungal agents [27], the influence of UV-B radiation on this pathway may lead to changes in plant–pathogen interactions.
Two of the three sterols identified as differential compounds according to volcano analysis (4α-formyl-4β-methyl-5α-cholesta-8,24-dien-3β-ol and 4α-carboxy-5α-cholesta-8,24-dien-3β-ol) markedly accumulated in UV-B-treated leaves after three days of recovery, similarly to the triterpene precursor, presqualene diphosphate (Table S3). The two sterols identified in the present research are cholesterol precursors. Differently from mammals, who synthesize cholesterol, plants produce a complex mix of sterols, among which sitosterol, stigmasterol, and campesterol are the most representative compounds, though, in some plants, as in the Solanaceae family, cholesterol is present in significant amounts [28]. Sterols are intrinsic components of the cell membrane, therefore their increase after the end of UV-B treatment may represent a mechanism of repairing possible damage induced by UV-B radiation to restore membrane fluidity and permeability. A slightly increased lipid peroxidation was indeed observed in UV-B-irradiated Micro-Tom leaves on day 11, while, at the end of the recovery period, this parameter was similar between control and treated leaves [29]. It is well known that flavonoids can protect against lipid peroxidation caused by stressing conditions such as photoinhibition [30,31], and different studies have indicated that UV-B can increase lipid peroxidation levels, usually linked with a decrease in chlorophylls concentration and photosynthesis inhibition [32]. In particular, Liu et al. [33] found a significant correlation between anthocyanidins and lipid peroxidation, suggesting that the content of those antioxidant compounds may be linked to the extent of lipid oxidation. In this study, the FC of anthocyanin in leaves was slightly decreased in respect to CTR plants (paragraph 2.4), which is in line with the quantification reported in our recent paper [34] after 11 days. One could assume that flavonoids might have been consumed in reactions aimed to interrupt the oxidative process occurring at the expense of cellular components such as membranes. However, as we reported in our previous papers [29,34], the UV-B treatment applied in this study did not cause severe damages to the photosynthetic apparatus, as indicated by the higher total chlorophyll concentration detected after 11 days in treated plants, and no differences in respect to CTR, considering the PSII efficiency, were detected. Moreover, as showed in [34], we found a reduced expression of a gene related to anthocyanin biosynthesis on previous days, suggesting a likely reduced biosynthesis of these compounds in the leaves of treated plants.
An increase in sterol (sitosterol and stigmasterol) accumulation is reported in grapevine leaves following a ‘‘field-simulating’’ dose of UV-B radiation (4.75 kJ m−2 per day) administered at low intensity (16 h at 8.25 μW cm−2) and, less intensely, also at high intensity (4 h at 33 μW cm−2). This latter treatment induced an increase in antioxidant compounds, such as mono- and diterpenes and tocopherol and phytol [35]. In our study, the accumulation of some monoterpenes, diterpenes, and sesquiterpenes, among which were some phytoalexins, was stimulated in UV-exposed tomato leaves, indicating the activation of both MEP and mevalonate biosynthetic pathways. Such a response was generally more pronounced after 11 days of treatment (Table S3). Beside changes in sterol levels, in our research, UV-B treated leaves underwent a decrease in 1-18:3-2-16:2-monogalactosyldiacylglycerol and an increase in the amount of 1-16:0-2-18:3-digalactosyldiacylglycerol, particularly after 3 days of recovery, indicating that UV-B radiation modifies the lipidic composition of membranes.
After 3 days of recovery, hormone degradation was strongly increased (+43 log2FC), suggesting the necessity of modulating the hormone homeostasis. Interestingly, in UV-11d leaves, the reduced secondary metabolites synthesis was parallel to a lower degradation process (Figure 3). A more detailed representation of the secondary metabolites pathway, which is the typical process involved in the UV response and acclimation, is shown in Figure S1, where, in particular, the phenylpropanoid biosynthetic and degradation pathways demonstrated a strong reduction in treated leaves.
Roots displayed a different behaviour from leaves, in particular, in the case of the secondary metabolism biosynthesis, which appeared to be enhanced in UV-11d in respect to its CTR (+191 log2FC, Figure 3). After examining the secondary metabolites in detail, pathway biosynthesis led us to identify an increase in terpenoid (+60 log2FC) and phenylpropanoid derivative (+16 log2FC) biosynthesis (Figure S1) in irradiated plants with respect to CTR at the 11-day treatment. However, phenylpropanoid biosynthesis was greatly reduced after the recovery (−47 log2FC). Degradative reactions of the secondary metabolites (Figure 3) showed a reduction in the fold-change after the recovery in respect to CTR roots (−13 log2FC). Concerning the carbohydrate metabolism, degradative processes were strongly enhanced on both day 11 of treatment and after the recovery (about +13 log2FC), with respect to roots of control plants, likely suggesting a mobilization of storage forms. The degradation of root carbohydrates might have been the source of precursors needed for the accumulation of secondary metabolites in this organ such as the flavonoid class, which will be discussed later in this article.
As observed in leaves, two monogalactosyldiacylglycerols showed an opposite behaviour following UV-B treatment also in roots. This finding, together with the decreased 4a-carboxy-5a-cholesta-7,24-dien-3b-ol root content detected in the roots of UV-B-treated plants, in accordance with the report of decreased cholesterol and accumulation of two steryl esters in roots of Withania somnifera plants irradiated with supplemental UV-B (3.6 kJ m−2 day−1 above ambient level [36], suggests a rearrangement of membrane lipid composition due to UV-B irradiation also in the non-irradiated organ.
Based on the unexpected modulation of the biochemical processes particularly related to the phenylpropanoid pathway, phenolic compounds were profiled in tomato leaves and roots, since it is highly recognised that phenolic compounds are involved in many plant processes and play a fundamental role in the acclimation towards UV, especially UV-B radiation.

2.3. Effects of Short Daily UV Treatment on the Phenolic Profile of Tomato Leaves and Roots

Since the UHPLC-ESI/QTOF-MS system coupled with a comprehensive phenolic database (Phenol-Explorer) allowed for the detection and identification of hundreds of compounds (232 and 183 compounds in leaves and roots, respectively), an OPLS-DA model was built for both organs, and for the metabolomics, to extrapolate the UV-induced phenolic modifications.
In leaves, the OPLS-DA analysis (Figure 4a) revealed a clear distinction between CTR and UV-treated groups, indicating that, regardless of the sampling time, the irradiation did result in effectively modulating the phenolic pattern.
Variable Importance for Prediction (VIP) analysis was carried out to show which phenolic compounds weighed the most in the loading plot of the OPLS-DA model reported in Figure 4a. Table 3 shows the 33 phenolic markers with a VIP score higher than 1.1, classified according to their class and subclass.
The most abundant class was the flavonoid one (15 compounds in total) followed by the phenolic acids (12), while anthocyanins (5), flavonols (5), and hydroxycinnamic acids (10) were the most represented subclasses. At the same time, to understand which of the VIP markers displayed a significant fold-change between UV and CTR foliar samples, a univariate statistical analysis was performed and highlighted three compounds with a significant p value. These are rosmarinic acid (−4 log2FC), caffeoyl aspartic acid (+0.52 log2FC), and cirsimaritin (+0.72 log2FC). Besides the UV protection of phenolic acids and flavonoids, they also play an important role as antifungal and antimicrobial compounds.
Furthermore, both sampling time and UV treatment played a role in the group separation along with the vectors (Figure 4b). Indeed, the first vector provided a clear separation of all groups, with CTR-11d being the only one located in the positive half of the plot, while, according to the second vector, only UV-11d was positioned in the positive half. However, upon obtaining a general overview of the OPLS-DA output, control, and UV-treated samples, they are shown to be well clustered within each group and well separated from the others, indicating good homogeneity among the replicates and a highly different phenolic profile according to both the treatment and the sampling time.
The OPLS-DA analysis on root tissue revealed a clear distinction between control and UV-treated plants (Figure 4c). This result suggests that the UV irradiation of the plant’s aerial part did indirectly alter the phenolic profile also within the below-ground tissues. The details of the VIP analysis performed on root tissue data are indicated in Table 4, where 25 phenolics with a VIP score higher than 1.1 are classified according to their class and subclass.
The most abundant class was the flavonoid one (12 compounds in total) followed by the phenolic acids (7), while anthocyanins (6) and hydroxycinnamic acids (5) represent the large subclasses. As for leaves, a univariate statistical analysis of the FC was performed, and p-Coumaroyl glycolic acid, a hydroxycinnamic acid, was the only compound with a significant p value. This phenolic acid was significantly reduced by UV treatment (−10 log2FC) and, as reported by Kadam et al. [37], it is known for its antioxidant and anti-inflammatory properties.
Moreover, considering both time and UV treatment factors (Figure 4d), a good separation among all groups occurred within the score plot. In detail, the first vector highlights the effect of the sampling time factor, since both CTR-11d and UV-11d plants were overlapped and placed in the negative half of the plot, while CTR-rec and UV-rec were in the positive half, although well distinct from each other. According to the second vector, the UV-rec group resulted in being the only one located in the positive half of the plot, distant from all the other groups.

2.4. Fold-Change Analysis on the Comprehensive Dataset of Leaves and Roots Phenolic Compounds

To investigate the accumulation of the most UV-responsive phenolic classes and subclasses in detail, a fold change (FC) analysis of UV in respect to controls was carried out on the base of the whole phenolics dataset, considering both the sampling times (11 days of UV exposure and 3 days after the end of the treatment) and both tissues (leaves and roots) (Figure 5). A detailed list of the phenolic molecules and their fold change are listed in Tables S4 and S5.
A general decrease in the phenolic class abundance was detected in the UV-11d leaves, especially for phenolic acids, which showed a −0.57 log2FC in respect to control leaves (Figure 5a). After the 3 days of recovery, most classes were still down accumulated in UV-rec, with flavonoids and phenolic acids showing a −0.13 and −0.18 log2FC and lignans being the most affected phenolic subclass, with −0.68 log2FC in respect to control (Figure 5a). UV-induced responses typically trigger the accumulation of phenylpropanoids due to their ROS scavenging and screening activity; although, recently, several studies have showed the complexity of the phenolic response to this wavelength, depending on the exposure time, plant species, and phenolic class considered [38]. However, after 11 days of treatments, Micro-Tom tomato plants have already acclimated to the supplemental UV irradiation, as showed in our previous study [29], where many physiological parameters were checked to ensure that the UV dose applied did not induce stress-related responses. Indeed, we detected a reduced content of total phenols on day 11 of UV (according to the Folin–Ciocalteau method), and no differences were found in terms of flavonoid concentration. In contrast, a slightly increased phenols and flavonoids content was detected in previous days. According to these data, we speculate that the well-known and widely described induction of phenolics biosynthesis could be an early response under UV-B conditions, with fluctuations depending on the antioxidant status of the plant. After the plant acclimated to low UV-B doses and the morphogenic responses took place, e.g., after an increase in leaf thickness, lignification, trichome density [39,40,41], it is likely that the plant energy and resources are re-addressed mainly to primary rather than secondary metabolism, so that physiological growth and development can be restored. Indeed, as reported by Coffey et al. [42], the inhibiting effect of UV-B radiation on the leaf growth of A. thaliana, grown in outdoor conditions for 7 months, ten days a month, was detected only under summer intensities, and it was not related to the activation of the UVR8-signalling pathway. The authors also found no changes in the concentration of UV-absorbing pigments; however, nothing was reported about the composition of the pool. Robson et al. [16] also argued that the UV effects on leaf growth are complex and that, once the plants acclimated to the UV stimulus, normal development can be restored, leading to compensatory effects on the enlargement of the leaf area. According to our hypothesis of a UV-triggered modification in the morphology of the leaves, the reduced levels of lignans detected in UV-rec leaves suggest that these compounds were cross-linked for lignin biosynthesis [43].
Very few studies have focused on UV-B perception in roots [15,44], and most of them have applied high UV-B intensities, which likely induced the onset of stress responses overlapping the activation of the UVR8 pathway [45,46]. To the best of our knowledge, no previous works have investigated the changes in the phenolic profile of roots under short daily UV conditions. The fold-change analysis on the root phenolic abundance dataset (Figure 5b) revealed a different behaviour between the two plant organs. The roots of treated plants showed an increase in flavonoids (+0.35 log2FC) and lignans (+0.26 log2FC) in the UV-11d group. After the recovery, an opposite behaviour was shown by lignans (−0.27 log2FC) and phenolic acids (−0.60-fold) whose concentrations were both reduced, while flavonoids were still increased (+0.64-fold).
Consequently, to better understand how the treatment influenced the flavonoids (Figure 5c,d) and phenolic acids (Figure 5e,f), the fold change analysis was also carried out at subclass levels in both organs.
Flavonoids in UV-11d leaves showed an increased log2FC of dihydrochalcones, flavanols, flavanones, and flavones (+0.25, +0.91, +1.47, +0.93, respectively), while anthocyanins, dihydroflavonols, and isoflavonoids were reduced (−0.16, −0.81 and −0.43, respectively) compared to the correspondent control. After the 3-day recovery, leaves of treated plants still exhibited an increase in flavanols and flavones (+0.44 and +0.22 log2FC, respectively) and a persistent decrease in anthocyanins and dihydroflavonols (−0.13 and −0.68, respectively). Conversely, flavanones (−0.22 log2FC), flavonols (−0.30), and isoflavonoids (+1.45) demonstrated an opposite behaviour in respect to what was observed after 11 days. In root tissue, all the flavonoid subclasses generally increased their abundance at both sampling times, with flavones showing the highest increase (+0.74 and +1.52 log2FC in UV-11d and UV-rec, respectively). Only isoflavonoids were reduced after the recovery period, with a −1.69-fold decrease in respect to controls.
Besides their function as antioxidants, flavonoids can modulate the phytohormone signalling as well. Polar auxin (IAA) transport was increased in studies with Arabidopsis mutant, with reduced levels of flavonoids aglycones [47,48,49]. In our study flavonols and flavones, which have been demonstrated to be able to inhibit polar IAA transport and signalling [50,51], were found to be enhanced in UV-11d leaves, likely causing a reduced basipetal movement of this hormone towards roots. Indeed, a reduced level of IAA concentrations was detected in the roots of UV-treated plants after 3 days of recovery [29], although a reduced biosynthesis in this organ could be another factor affecting its amount. Although less concentrated in respect to the aerial part of the plant, root flavonoids have significant roles in the regulation of root growth, within the nitrogen cycle, and among allelopathic interactions [52]. Flavonoids can exert signalling actions in lower concentrations than those necessary for ROS scavenging and UV absorption [50,53]. The increased levels of flavonols and flavones detected in roots of UV-11d may play a role in creating an IAA gradient at cellular and tissue levels, thus regulating root growth. However, since the flavonoid aglycones are located in tissue involved in auxin translocation, it is reasonable to consider that 8 flavonol glycosides out of 11 flavonols were found, and only 1 flavone glycoside out of 6 flavones in leaves, while, in roots, there were 6 flavonol glycosides out of 9 flavonols and 1 flavone glycoside out of 5 flavones. Further investigations are needed to verify the link between the IAA metabolism and flavonoids under short daily UV treatments, considering both leaves and roots, to gain a complete understanding of this phenomenon. The presence of flavonoids has been reported in the root exudates released into the rhizosphere, where these compounds act as regulators of nutrient availability, thanks to their reducing potential and chelating ability and promotion of root-rhizobia symbiosis [54]. Though we did not study the flavonoid content and composition of the root exudated, the increased flavonoid content (FC) in the roots of treated plants, detected in our experiment, suggests a potential benefit for irradiated plants in minerals uptake.
Concerning the phenolic acids, a general decrease in hydroxybenzoic (HB) and hydroxycinnamic (HC) acids occurred in leaves at both UV-11d and UV-rec. Notably, the strongest decrease was detected at 11 days for HC (−0.60 log2FC) and after the 3-day recovery period for HB (−1.43-fold). In a study on the effects of UV-modulated resistance to herbivory in tomato plants (0.34 kJ m−2 day−1 UV dose), it was found that chlorogenic acid and rutin, the two major leaf phenolic compounds in this species, but generally the whole leaf metabolome, did not change after 14-day UV exposure [55]. This might suggest that the secondary metabolism could have been affected at earlier days of the UV treatment. HC and HB in roots displayed an opposite trend, with −1.29 and +0.61 log2FC for HB and HC at 11 days, respectively, and +0.87 and −0.48 log2FC for HB and HC after the 3-day recovery period. As already mentioned, phenolic acids are involved in plant growth, antioxidant defence, and plant/microbe interaction [56,57]. Indeed, these compounds can be secreted by roots and act as allelochemicals by suppressing the proliferation of pathogens [58].

2.5. Outputs of the Volcano Analysis on the Phenolic Profile

A volcano analysis was performed on the phenolic dataset, combining moderated t-test (p < 0.05) and fold-change analysis (cut-off ≥ 1.2), and the differential compounds according to each sampling time were grouped in their respective class and subclass. The output revealed a similar trend between leaves and roots in the number of differentially modulated compounds (Table 5 and Table 6).
Indeed, a higher number of up and down accumulated compounds was found after the recovery, indicating that the short daily UV treatment applied for 11 days affected the plants not only during the daily irradiation, but also after 3 days of recovery.
Moreover, the effect on the UV-driven modulation of the phenolics differed between the two organs. Indeed, in leaves, 14 phenolics increased, while 10 phenolics decreased, following the UV treatment (Table 5). Considering the roots, only three phenolics increased, while 13 phenolics decreased, by UV exposure (Table 6).

2.6. Effects of Short UV Treatment on the Gene Expression of Tomato Leaves and Roots

To integrate the timing of response between the UV-mediated gene modulation and phenolics accumulation, leaves and roots were sampled at 8 d and 11 d of UV treatment and 3 days after the end of the UV irradiation (rec). For this purpose, some candidate genes related to the UVR8 specific pathway, HY5 and COP1, were investigated due to their well-known implication in the responses to low UV-B doses [4]. Simultaneously, genes related to the phenylpropanoid biosynthetic pathway, PAL, CHS, and FLS, were tested to find a link with the changes in the phenolic profile.

2.6.1. Differences in the Gene Expression of UVR8 Pathway-Related Genes

Leaf HY5 and COP1 gene expression is reported in Figure 6a,c. As expected, the HY5 transcript level increased at both 8 d and 11 d days (+2.1- and +3.88-fold, respectively), while no differences compared to control plant levels were found three days after the end of the UV exposure. This trend confirmed the activation of the specific UVR8 pathway in Micro-Tom tomato plants by the application of short UV treatment. Indeed, HY5 is known to play a pivotal role in light signalling-mediated responses, and it is also known to be induced by UV-B radiation [4].
Considering the root tissue, a decrease in HY5 expression was observed in treated plants (Figure 6b), being statistically significant at 8 d of UV and after the recovery (−0.53 and −0.18-fold, respectively). To the best of our knowledge, the expression of HY5 in roots not directly exposed to short daily UV doses has never been investigated in tomatoes during a long-term experiment. Moreover, in general, studies on the roles and functions of HY5 in roots are scarce. HY5 is known to take part in many biological processes, such as nitrogen assimilation and other activities related to light and hormone control, as in the case of auxin [4,59]. In a study on Arabidopsis, Cluis et al. [59] found that two negative regulators of auxin signalling, AUXIN RESISTANT 2/INDOLE ACETIC ACID 7 and SOLITARY ROOT/IAA14, were downregulated in hy5 mutants, suggesting that HY5 might regulate their transcription by binding their promoters. In addition, the same authors stated that hy5 mutants have higher IAA signalling due to the HY5-regulating effect on the expression of those IAA inhibitors. From our data, we speculate that the reduced expression of the HY5 gene in the roots of UV-treated plants could potentially lead to changes in the IAA signalling and alterations in the morphology of this organ. Further studies regarding this aspect are encouraged.
No modification in COP1 gene expression was detected in both leaves and roots (Figure 6c,d), suggesting that its activation could have occurred in the early days of the treatment, or that it might require some recovery time after the treatment to be detectable. However, in accordance with our results, a study by Liu et al. [12] on the gene regulation of flavonol synthase by high and low UV-B in Vitis vinifera L. grapes found no significant responses COP1 to UV-B.

2.6.2. Differences in the Gene Expression of Some Phenylpropanoid Biosynthetic Genes

The UV effects on the phenolic biosynthetic pathways were checked by measuring the gene expression of PAL, CHS, and FLS, whose enzymes are involved in the early steps of phenolic acid, flavonoids, and flavonols biosynthesis, respectively.
PAL, an enzyme catalysing the cinnamic acid biosynthesis, which is a precursor for lignin and flavonoids, plays a pivotal role in regulating compounds resulting from the phenylpropanoid pathway and as a branch point between aromatic amino acids and secondary metabolites pathway [60]. The transcripts levels of PAL were increased (+1.94-fold) in leaves of treated plants at 11 d (Figure 6e) in respect to controls, while no differences were found on the other days. This is in line with the differences in the fold-change of phenolic acids’ content, presented in Figure 5a, and which, despite being less concentrated than in the control leaves and being less decreased in UV-rec than in UV11d, indicate a trend to recover the control values, likely due to an upregulation of the related biosynthetic enzymes. In roots (Figure 6f), an increase in the PAL transcript level of treated plants was detected at 8 d (+2.76-fold), followed by a reduction in the latter days (−0.44 and −0.73-fold at 11 d and rec, respectively). The changes in PAL transcript level occurring in roots reflect the changes in the phenolic acids’ content (Figure 5b). Remarkably, the increase in PAL gene expression at 8 d led to an increase in the phenolic acids’ content in UV-11d, while the reduced concentration of these compounds after the 3 days of recovery matched the reduced PAL levels detected at 11 d. Moreover, a persistently reduced amount of root phenolic acids was likely to be detected in the later days, since the downregulation of PAL transcription was still visible after the end of the treatment.
CHS is an enzyme involved in the early steps of flavonoid biosynthesis, catalysing the conversion from coumaroyl-CoA to naringenin chalcone, while FLS catalyses the conversion from dihydroflavonoids to flavonoids. The literature states that the increase in CHS expression represents a UV-B-specific response, upregulated by very short UV-B exposures [61]. Moreover, Brown and Jenkins [4] have indicated that changes in CHS expression do not occur when non-specific UV-B responses take place. Foliar CHS transcript levels (Figure 6g) were reduced by UV treatment at both 8 d and 11 d (−0.65 and −0.69-fold, respectively) and, although the flavonoid content in UV-11d was comparable to controls (Figure 5a), lower CHS levels might be the cause of the reduction in the flavonoid amount detected in UV-rec (Figure 5a). After the recovery, CHS gene expression slightly increased (+0.19-fold) compared to control samples, suggesting the activation of the flavonoid biosynthesis and their accumulation in the following days. No differences were found in foliar FLS gene expression (Figure 6i) according to the days considered, even if a decreased flavonol accumulation was evident in UV-rec (Figure 5a).
In roots, CHS transcript levels were reduced (−0.65-fold) at the 8 d (Figure 6h), followed by a 3-fold increase at 11 d and rec, along with the increase in flavonoid concentration (Figure 5b). Root FLS gene expression (Figure 6j) exhibited an increase in UV-treated plants at 8 d (+0.58-fold), possibly leading to a slight increase in the flavonols detected in the UV-11d root samples (Figure 5b).
Generally, a good relationship between gene expression and biochemical analysis was provided. One of the main results worth highlighting is that that the UV irradiation applied on the plant’s aerial part caused the activation of some genes related to the phenolic biosynthetic pathway of the root. As mentioned before, in the case of UV, no light piping has been observed so far, contrary to far red light [24]; therefore, one might hypothesise that some kind of biochemical signals moved to the root organ to provide a coordination between the upper and lower ground parts of the plants.
Besides a local accumulation mediated by enhanced gene expression, a long-distance transport from leaves to root cannot be excluded. As demonstrated by Buer et al. [62] by the local application of a fluorescent flavonoid and its detection in distal tissues in Arabidopsis, flavonoids accumulated in light-shielded roots of plants with light-exposed shoots thanks to a shoot-to-root transport.
However, it should be remembered that the metabolites level results from both biosynthetic and degradation processes, and the overall decrease in phenolic acids and flavonoids in leaves could be also due to their consumption in ROS scavenging activity. Additional time points on the earlier days are highly encouraged to have a comprehensive view of the acclimation process of Micro-Tom tomato plants under short daily UV treatment.

3. Materials and Methods

3.1. Plant Cultivation and UV Treatment

Surface sterilized seeds of Solanum lycopersicum L. cultivar Micro-Tom (JustSeed Ltd., Wrexham, United Kingdom) were germinated on water-soaked paper and then transferred in pots with perlite for 7 days, after which the plantlets were placed in a hydroponic system and grown in a Hoagland solution (pH~6), refreshed every week. The climate chamber was set at 24 °C, with 16 h light/8 h dark photoperiod and a photosynthetic photon flux density (PPFD) of 228 μmol m−2 s−1 supplied by an LED illumination system (blue/red/green 3/6/1 ratio, C-LED, Imola, Italy). The UV treatment was applied after 24 days and lasted 11 days. Plantlets were divided into a control group (CTR), receiving PAR radiation only, and a UV-treated group (UV), grown with PAR + UV radiation (15 min a day corresponding to a daily unweighted UV dose of 1.19 kJ m−2) provided by Philips Ultraviolet-B Narrowband lamps (TL 20 W/01-RS, Koninklijke Philips Electronics, Eindhoven, The Netherlands). The UV bulbs emission, measured by a JAZ EL-XR1 spectroradiometer (OCEAN OPTICS, Dunedin, FL, USA), provided a UV-B biologically effective intensity of 0.132 W m−2 and UV-A intensity of 0.015 W m−2, calculated according to Flint and Caldwell [63], with the latter corresponding to only 36% contribution of the total UV emission of the lamps. Different shading between CTR and UV plants was avoided through the application of empty lamp holders inside the CTR chamber, and all plants were rotated within the chambers to obtain uniform irradiation.
For phenolic profiling, samples of leaves and roots of both control and treated groups were collected after 11 days of short daily UV treatment (CTR-11d and UV-11d) and after 3 days from the end of the irradiation (CTR-rec and UV-rec). For gene expression analysis, samples were collected after 8 (8 d) and 11 (11 d) days of UV exposure and after 3 days from the end of the irradiation (rec). A pool of leaves and the whole root were immediately frozen in liquid nitrogen and stored at −80 °C until use.

3.2. Extraction and UHPLC-ESI-QTOF-MS Analysis of Leaves and Roots Metabolites

Freeze-dried samples (100 mg) of leaves and roots were extracted by homogenizer-assisted solvent extraction (Ultra-Turrax; Polytron PT, Switzerland) in 4 mL of an acidified methanol/water solution (80:20, v/v), centrifuged and filtered through a 0.22 µm membrane. For the analysis, 4 biological replicates per treatment and time, each obtained from a pool of 3 different plants, were used. The analysis of the metabolomics and phenolic profile was performed as reported by Senizza et al. [64] by using an Agilent 6550 iFunnel quadrupole-time-of-flight mass spectrometer and an Agilent 1200 series ultra-high-pressure liquid chromatographic system (UHPLC-ESI/QTOF-MS). Reverse phase chromatographic separation was achieved under a water-acetonitrile gradient elution, starting from 6% acetonitrile to 94% in 33 min. The mass spectrometer worked in SCAN mode with a range from 100 to 1200 m/z, positive and negative polarity, and extended dynamic range mode (nominal mass resolution = 30,000 FWHM).
Metabolite annotation was performed by the Profinder B.07 software tool (Agilent Technologies, Santa Clara, CA, USA) by their monoisotopic accurate mass and isotopes patterns (accurate spacing and isotopes ratio) according to [64], using the databases Phenol-Explorer 3.6 and PlantCyc to annotate phenolic compounds and plant metabolome, respectively [65,66]. Identification was carried out in compliance with Level 2 (putatively annotated compounds), as set out by the COSMOS Metabolomics Standards Initiative [67].

3.3. RNA Extraction and qRT-PCR Analysis

RNA from frozen homogenized plant tissue (3 biological replicates per treatment and time, each obtained from a pool of 3 different plants per treatment and time) was extracted from leaves and roots separately and prepared using Qiagen RNeasy Plant Mini kits (QIAGEN S.r.l., Milan, Italy) according to the manufacturer’s instructions and followed by cDNA synthesis using the cDNA Verso kit (Thermo Fisher Scientific-Life Technologies, Carlsbad, CA, USA). The SyGreen qPCR kit (qPCRbio; PCR Biosystems, London, UK) was used to perform qRT-PCR in an iCycler thermal cycler with an iQ5 optical system (Bio-Rad, Hercules, CA, USA) as an optical module for signal detection. The data were analysed with QBase software (Biogazelle, Ghent, Belgium) [68], where expression values were normalised with reference genes ACTIN7, UBIQUITIN, and EF1. The relative expression was automatically calculated in Qbase+. Primer sequences for HY5, COP1, PAL, CHS, FLS, and reference genes are listed in Table S1.

3.4. Data Elaboration and Statistics

Elaboration of metabolomics data was carried out as previously reported by Ghisoni et al. [69]. The software Mass Profiler Professional 12.6 (Agilent Technologies, Santa Clara, CA, USA) was used for log2 transformation, normalization at the 75th percentile, and baselining against the median of each compound. Thereafter, the unsupervised hierarchical cluster analysis was used to underline the relatedness across the different treatments, according to Euclidean distance and Ward’s linkage rule. The volcano plot analysis was also performed for each treatment by combining a moderated t-test (p < 0.05; Benjamini–Hochberg multiple testing correction) and fold-change analysis (cut-off ≥ 1.2). The SIMCA 16 (Umetrics, Malmo, Sweden) software was then used to build the orthogonal projection to latent structures discriminant analysis (OPLS-DA) supervised modelling. Outliers were investigated according to Hotelling’s T2 (95% and 99% confidence limits for suspect and strong outliers), CV-ANOVA (p < 0.05) cross-validation and permutation testing (for overfitting, N = 100) were also performed from the OPLS-DA model. Thereafter, model fitness parameters (goodness-of-fit, R2Y, and goodness-of-prediction, Q2Y) were recorded, and the variable importance in projection (VIP) analysis was used to identify the most discriminant compounds (VIP score > 1.1). Differential compounds were finally imported into the PlantCyc Omic Viewer Pathway Tool (Stanford, CA, USA) software for biochemical interpretations [70].
JMP software (SAS Institute, Inc., Cary, NC, USA) was used for gene expression data to perform Student t-test statistical analysis. According to each day (8 d, 11 d and rec), asterisks (*) indicate a significant difference between CTR and UV group (* p ≤ 0.05, ** p < 0.01, *** p < 0.001).

4. Conclusions

To the best of our knowledge, this is the first time that a comparative metabolomic and phenolic profiling analysis has been performed in leaves and roots following UV irradiation. In the present work, we focused on the late responses of tomato plants acclimated to short daily UV treatment. The untargeted approach revealed an extensive modulation of the metabolome in leaves. Nonetheless, the root metabolome was also significantly affected by UV, despite not being directly exposed to the radiation, thus giving a new insight on the indirect effects in this organ, indicative of shoot-to-root signalling. Our results demonstrate that both biosynthesis and degradation processes were affected, involving AA and FA/lipids, in addition to a modulation of the hormone profile in leaves as compared to non-irradiated plants. Moreover, roots displayed a marked increase in the carbohydrate degradative pathway. However, a particular trend related to the secondary metabolism was noticed in both organs. Indeed, the biosynthesis of secondary metabolites appeared reduced in the leaves of treated plants. This result may indicate that, once the plant had acclimated to short daily UV treatment, the compounds produced first as part of the antioxidant defence could be used for other purposes, e.g., cross-linked for lignin synthesis in the case of lignans, thus causing changes in the leaf morphology and development. Roots showed an opposite behavior, with an enhancement of the secondary metabolism, particularly after the 11 days of treatment, where an increase in flavonoid content was found. These changes might affect not only the nutrient uptake, but also the architecture of this organ. Further studies are encouraged to verify this hypothesis. Our data also suggest a prolonged effect on the modulation of phenolic compounds in both organs once the UV treatment had ceased.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants11141829/s1. Figure S1: Fold-change (FC) of leaves and roots secondary metabolites biosynthesis and degradation in leaves and roots of UV treated plants, at 11 days of treatment and after 3 days of recovery, with respect to control plants. FC values were elaborated by the Omics Viewer Dashboard of the PlantCyc Pathway Tool software. Data for each biochemical pathway represent the sum of Log2FC of each detected compound. Syn, biosynthesis; Deg, degradation; Table S1: Primers list for qRT-PCR in Micro-Tom tomato plants (Solanum lycopersicum L.); Table S2: List of metabolites detected with UHPLC-ESI-QTOF in leaves and roots; Table S3: Differential compounds in leaves and roots resulting from volcano plot; Table S4: Dataset of foliar phenolic compounds using the database exported from Phenol Explorer. The average of the compound intensity is provided (where 1 indicated not detected), together with annotations (raw formula and ID score); Table S5: Dataset of root phenolic compounds using the database exported from Phenol Explorer. The average of the compound Intensity is provided (where 1 indicated not detected), together with annotations (raw formula and ID score).

Author Contributions

A.M.: investigation, formal analysis, data curation, visualization, writing—original draft and review and editing. M.S.: formal analysis, data curation, visualization, writing review, and editing. L.V.: supervision, writing, review, and editing. M.C.S.: investigation, review, and editing. M.B.M.-M.: methodology, investigation, formal analysis, data curation, visualization. L.Z.: methodology, investigation, formal analysis, data curation, visualization. L.L.: methodology, supervision, writing, review, and editing. M.F.Q.: conceptualization, supervision, writing, review, and editing. D.V.D.S.: supervision, writing, review, and editing. A.C.: conceptualization, supervision, writing, review, and editing. A.R.: conceptualization, project administration, supervision, writing, review, and editing. All authors have read and agreed to the published version of the manuscript.

Funding

A.M. was financially supported by the PhD program of Agriculture, Food and Environment fellowship, University of Pisa. D.V.D.S. acknowledges Ghent University for financial support.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets presented in this study are available in Supplementary Materials.

Acknowledgments

We thank C-LED (Imola, Italy) for providing us with the LEDs equipment.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Jenkins, G.I. Signal transduction in responses to UV-B radiation. Annu. Rev. Plant Biol. 2009, 60, 407–431. [Google Scholar] [CrossRef] [PubMed]
  2. Hideg, É.; Jansen, M.A.; Strid, Å. UV-B exposure, ROS, and stress: Inseparable companions or loosely linked associates? Trends Plant Sci. 2013, 18, 107–115. [Google Scholar] [CrossRef] [Green Version]
  3. Ulm, R.; Baumann, A.; Oravecz, A.; Máté, Z.; Ádám, É.; Oakeley, E.J.; Schäfer, E.; Nagy, F. Genome-wide analysis of gene expression reveals function of the bZIP transcription factor HY5 in the UV-B response of Arabidopsis. Proc. Natl. Acad. Sci. USA 2004, 101, 1397–1402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Brown, B.A.; Jenkins, G.I. UV-B signaling pathways with different fluence-rate response profiles are distinguished in mature Arabidopsis leaf tissue by requirement for UVR8, HY5, and HYH. Plant Physiol. 2008, 146, 576–588. [Google Scholar] [CrossRef] [Green Version]
  5. Kusano, M.; Tohge, T.; Fukushima, A.; Kobayashi, M.; Hayashi, N.; Otsuki, H.; Kondou, Y.; Goto, H.; Kawashima, M.; Matsuda, F.; et al. Metabolomics reveals comprehensive reprogramming involving two independent metabolic responses of Arabidopsis to UV-B light. Plant J. 2011, 67, 354–369. [Google Scholar] [CrossRef] [PubMed]
  6. Rizzini, L.; Favory, J.J.; Cloix, C.; Faggionato, D.; O’Hara, A.; Kaiserli, E.; Baumeister, R.; Schäfer, E.; Nagy, F.; Jenkins, G.I.; et al. Perception of UV-B by the Arabidopsis UVR8 protein. Science 2011, 332, 103–106. [Google Scholar] [CrossRef] [Green Version]
  7. Kaiserli, E.; Jenkins, G.I. UV-B promotes rapid nuclear translocation of the Arabidopsis UV-B–specific signaling component UVR8 and activates its function in the nucleus. Plant Cell 2007, 19, 2662–2673. [Google Scholar] [CrossRef] [Green Version]
  8. Favory, J.J.; Stec, A.; Gruber, H.; Rizzini, L.; Oravecz, A.; Funk, M.; Albert, A.; Cloix, C.; Jenkins, G.I.; Oakeley, E.J.; et al. Interaction of COP1 and UVR8 regulates UV-B-induced photomorphogenesis and stress acclimation in Arabidopsis. EMBO J. 2009, 28, 591–601. [Google Scholar] [CrossRef]
  9. Kaspar, S.; Matros, A.; Mock, H.P. Proteome and flavonoid analysis reveals distinct responses of epidermal tissue and whole leaves upon UV−B Radiation of barley (Hordeum vulgare L.) seedlings. J. Proteome Res. 2010, 9, 2402–2411. [Google Scholar] [CrossRef] [PubMed]
  10. Burchard, P.; Bilger, W.; Weissenböck, G. Contribution of hydroxycinnamates and flavonoids to epidermal shielding of UV-A and UV-B radiation in developing rye primary leaves as assessed by ultraviolet-induced chlorophyll fluorescence measurements. Plant Cell Environ. 2000, 23, 1373–1380. [Google Scholar] [CrossRef]
  11. Kalbin, G.; Ohlsson, A.B.; Berglund, T.; Rydström, J.; Strid, A. Ultraviolet-B-radiation-induced changes in nicotinamide and glutathione metabolism and gene expression in plants. Eur. J. Biochem. 1997, 249, 465–472. [Google Scholar] [CrossRef] [Green Version]
  12. Liu, L.; Gregan, S.; Winefield, C.; Jordan, B. From UVR 8 to flavonol synthase: UV-B-induced gene expression in Sauvignon blanc grape berry. Plant Cell Env. 2015, 38, 905–919. [Google Scholar] [CrossRef] [PubMed]
  13. Heijde, M.; Ulm, R. UV-B photoreceptor-mediated signalling in plants. Trends Plant Sci. 2012, 17, 230–237. [Google Scholar] [CrossRef]
  14. Tong, H.; Leasure, C.D.; Hou, X.; Yuen, G.; Briggs, W.; He, Z.H. Role of root UV-B sensing in Arabidopsis early seedling development. Proc. Natl. Acad. Sci. USA 2008, 105, 21039–21044. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Vanhaelewyn, L.; Bernula, P.; Van Der Straeten, D.; Vandenbussche, F.; Viczián, A. UVR8-dependent reporters reveal spatial characteristics of signal spreading in plant tissues. Photochem. Photobiol. Sci. 2019, 18, 1030–1045. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Robson, T.M.; Klem, K.; Urban, O.; Jansen, M.A. Re-interpreting plant morphological responses to UV-B radiation. Plant Cell Environ. 2015, 38, 856–866. [Google Scholar] [CrossRef]
  17. Wan, J.; Zhang, P.; Wang, R.; Sun, L.; Wang, W.; Zhou, H.; Xu, J. UV-B radiation induces root bending through the flavonoid-mediated auxin pathway in Arabidopsis. Front. Plant Sci. 2018, 9, 618. [Google Scholar] [CrossRef]
  18. Santin, M.; Ranieri, A.; Hauser, M.T.; Miras-Moreno, B.; Rocchetti, G.; Lucini, L.; Stride, Å.; Castagna, A. The outer influences the inner: Postharvest UV-B irradiation modulates peach flesh metabolome although shielded by the skin. Food Chem. 2021, 338, 127782. [Google Scholar] [CrossRef]
  19. Santin, M.; Lucini, L.; Castagna, A.; Rocchetti, G.; Hauser, M.T.; Ranieri, A. Comparative “phenol-omics” and gene expression analyses in peach (Prunus persica) skin in response to different postharvest UV-B treatments. Plant Physiol. Biochem. 2019, 135, 511–519. [Google Scholar] [CrossRef]
  20. Liu, Y.; Liu, J.; Wang, H.Z.; Wu, K.X.; Guo, X.R.; Mu, L.Q.; Tang, Z.H. Comparison of the global metabolic responses to UV-B radiation between two medicinal Astragalus species: An integrated metabolomics strategy. Environ. Exp. Bot. 2020, 176, 104094. [Google Scholar] [CrossRef]
  21. Kaling, M.; Kanawati, B.; Ghirardo, A.; Albert, A.; Winkler, J.B.; Heller, W.; Barta, C.; Loreto, F.; Schmitt-Kopplin, P.; Schnitzler, J.P. UV-B mediated metabolic rearrangements in poplar revealed by non-targeted metabolomics. Plant Cell Environ. 2015, 38, 892–904. [Google Scholar] [CrossRef] [PubMed]
  22. Tossi, V.E.; Regalado, J.J.; Iannicelli, J.; Laino, L.E.; Burrieza, H.P.; Escandón, A.S.; Pitta-Álvarez, S.I. Beyond Arabidopsis: Differential UV-B response mediated by UVR8 in diverse species. Front. Plant Sci. 2019, 10, 780. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, X.; Yao, Q.; Gao, X.; Jiang, C.; Harberd, N.P.; Fu, X. Shoot-to-root mobile transcription factor HY5 coordinates plant carbon and nitrogen acquisition. Curr. Biol. 2016, 26, 640–646. [Google Scholar] [CrossRef] [Green Version]
  24. Lee, H.J.; Ha, J.H.; Park, C.M. Underground roots monitor aboveground environment by sensing stem-piped light. Commun. Integr. Biol. 2016, 9, e1261769. [Google Scholar] [CrossRef]
  25. Conconi, A.; Miquel, M.; Browse, J.A.; Ryan, C.A. Intracellular levels of free linolenic and linoleic acids increase in tomato leaves in response to wounding. Plant Physiol. 1996, 111, 797–803. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Yang, B.; Wang, X.; Gao, C.; Chen, M.; Guan, Q.; Tian, J.; Komatsu, S. Proteomic and metabolomic analyses of leaf from Clematis terniflora DC. exposed to high-Level ultraviolet-B irradiation with dark treatment. J. Proteome Res. 2016, 15, 2643–2657. [Google Scholar] [CrossRef] [PubMed]
  27. Guimarães, A.; Venâncio, A. The Potential of Fatty Acids and Their Derivatives as Antifungal Agents: A Review. Toxins 2022, 14, 188. [Google Scholar] [CrossRef] [PubMed]
  28. Rogowska, A.; Szakiel, A. The role of sterols in plant response to abiotic stress. Phytochem. Rev. 2020, 19, 1525–1538. [Google Scholar] [CrossRef]
  29. Mannucci, A.; Mariotti, L.; Castagna, A.; Santin, M.; Trivellini, A.; Reyes, T.H.; Mensuali-Sodi, A.; Ranieri, A.; Quartacci, M.F. Hormone profile changes occur in roots and leaves of Micro-Tom tomato plants when exposing the aerial part to low doses of UV-B radiation. Plant Physiol. Biochem. 2020, 148, 291–301. [Google Scholar] [CrossRef]
  30. Chauhan, N.P.; Fatma, T.; Mishra, R.K. Protection of wheat chloroplasts from lipid peroxidation and loss of photosynthetic pigments by quercetin under strong illumination. J. Plant Physiol. 1992, 140, 409–413. [Google Scholar] [CrossRef]
  31. Liu, M.; Li, X.; Liu, Y.; Cao, B. Regulation of flavanone 3-hydroxylase gene involved in the flavonoid biosynthesis pathway in response to UV-B radiation and drought stress in the desert plant, Reaumuria soongorica. Plant Physiol. Biochem. 2013, 73, 161–167. [Google Scholar] [CrossRef] [PubMed]
  32. Tuo, C.; Lizhe, A.; Huyuan, F.; Jinghong, Y.; Xunling, W. The effect of UV-B radiation on membrane lipid peroxidation and mechanisms in broad bean leaves. Acta Ecol. Sin. 2001, 21, 579–583. [Google Scholar]
  33. Liu, M.; Cao, B.; Zhou, S.; Liu, Y. Responses of the flavonoid pathway to UV-B radiation stress and the correlation with the lipid antioxidant characteristics in the desert plant Caryopteris mongolica. Acta Ecol. Sin. 2012, 32, 150–155. [Google Scholar] [CrossRef]
  34. Mannucci, A.; Scartazza, A.; Santaniello, A.; Castagna, A.; Santin, M.; Quartacci, M.F.; Ranieri, A. Short daily ultraviolet exposure enhances intrinsic water-use efficiency and delays senescence in Micro-Tom tomato plants. Funct. Plant Biol. 2022. [Google Scholar] [CrossRef]
  35. Gil, M.; Pontin, M.; Berli, F.; Bottini, R.; Piccoli, P. Metabolism of terpenes in the response of grape (Vitis vinifera L.) leaf tissues to UV-B radiation. Phytochemistry 2012, 77, 89–98. [Google Scholar] [CrossRef]
  36. Takshak, S.; Agrawal, S.B. Alterations in metabolite profile and free radical scavenging activities of Withania somnifera leaf and root extracts under supplemental ultraviolet-B radiation. Acta Physiol. Plant. 2015, 37, 1–12. [Google Scholar] [CrossRef]
  37. Kadam, D.; Palamthodi, S.; Lele, S.S. LC–ESI-Q-TOF–MS/MS profiling and antioxidant activity of phenolics from L. Sativum seedcake. Food Sci. Technol. 2018, 55, 1154–1163. [Google Scholar] [CrossRef]
  38. Santin, M.; Ranieri, A.; Castagna, A. Anything New under the Sun? An Update on Modulation of Bioactive Compounds by Different Wavelengths in Agricultural Plants. Plants 2021, 10, 1485. [Google Scholar] [CrossRef]
  39. Jansen, M.A.; van den Noort, R.E.; Tan, M.A.; Prinsen, E.; Lagrimini, L.M.; Thorneley, R.N. Phenol-oxidizing peroxidases contribute to the protection of plants from ultraviolet radiation stress. Plant Physiol. 2001, 126, 1012–1023. [Google Scholar] [CrossRef] [Green Version]
  40. Yan, A.; Pan, J.; An, L.; Gan, Y.; Feng, H. The responses of trichome mutants to enhanced ultraviolet-B radiation in Arabidopsis thaliana. J. Photochem. Photobiol. B: Biol. 2012, 113, 29–35. [Google Scholar] [CrossRef]
  41. Karabourniotis, G.; Papadopoulos, K.; Papamarkou, M.; Manetas, Y. Ultraviolet-B radiation absorbing capacity of leaf hairs. Physiol. Plant. 1992, 86, 414–418. [Google Scholar] [CrossRef]
  42. Coffey, A.; Prinsen, E.; Jansen, M.A.K.; Conway, J. The UVB photoreceptor UVR8 mediates accumulation of UV-absorbing pigments, but not changes in plant morphology, under outdoor conditions. Plant Cell Environ. 2017, 40, 2250–2260. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, Y.; Chantreau, M.; Sibout, R.; Hawkins, S. Plant cell wall lignification and monolignol metabolism. Front. Plant Sci. 2013, 4, 220. [Google Scholar] [CrossRef] [Green Version]
  44. Yokawa, K.; Baluška, F. Pectins, ROS homeostasis and UV-B responses in plant roots. Phytochemistry 2015, 112, 80–83. [Google Scholar] [CrossRef] [PubMed]
  45. Ghisi, R.; Trentin, A.R.; Masi, A.; Ferretti, M. Carbon and nitrogen metabolism in barley plants exposed to UV-B radiation. Physiol. Plant. 2002, 116, 200–205. [Google Scholar] [CrossRef]
  46. Zhang, R.; Huang, G.; Wang, L.; Zhou, Q.; Huang, X. Effects of elevated ultraviolet-B radiation on root growth and chemical signaling molecules in plants. Ecotoxicol. Environ. Saf. 2019, 171, 683–690. [Google Scholar] [CrossRef] [PubMed]
  47. Murphy, A.; Peer, W.A.; Taiz, L. Regulation of auxin transport by aminopeptidases and endogenous flavonoids. Planta 2000, 211, 315–324. [Google Scholar] [CrossRef] [PubMed]
  48. Brown, D.E.; Rashotte, A.M.; Murphy, A.S.; Normanly, J.; Tague, B.W.; Peer, W.A.; Taiz, L.; Muday, G.K. Flavonoids act as negative regulators of auxin transport in vivo in Arabidopsis. Plant Physiol. 2001, 126, 524–535. [Google Scholar] [CrossRef] [Green Version]
  49. Jansen, M.A. Ultraviolet-B radiation effects on plants: Induction of morphogenic responses. Physiol. Plant. 2002, 116, 423–429. [Google Scholar] [CrossRef]
  50. Peer, W.A.; Blakeslee, J.J.; Yang, H.; Murphy, A.S. Seven things we think we know about auxin transport. Mol. Plant 2011, 4, 487–504. [Google Scholar] [CrossRef]
  51. Jacobs, M.; Rubery, P.H. Naturally occurring auxin transport regulators. Science 1988, 241, 346–349. [Google Scholar] [CrossRef] [PubMed]
  52. Rao, A.S. Root flavonoids. Bot. Rev. 1990, 56, 1–84. [Google Scholar] [CrossRef]
  53. Agati, G.; Tattini, M. Multiple functional roles of flavonoids in photoprotection. New Phytol. 2010, 186, 786–793. [Google Scholar] [CrossRef] [PubMed]
  54. Cesco, S.; Neumann, G.; Tomasi, N.; Pinton, R.; Weisskopf, L. Release of plant-borne flavonoids into the rhizosphere and their role in plant nutrition. Plant Soil 2010, 329, 1–25. [Google Scholar] [CrossRef]
  55. Escobar-Bravo, R.; Chen, G.; Kim, H.K.; Grosser, K.; van Dam, N.M.; Leiss, K.A.; Klinkhamer, P.G. Ultraviolet radiation exposure time and intensity modulate tomato resistance to herbivory through activation of jasmonic acid signaling. J. Exp. Bot. 2019, 70, 315–327. [Google Scholar] [CrossRef] [PubMed]
  56. Ghasemzadeh, A.; Ghasemzadeh, N. Flavonoids and phenolic acids: Role and biochemical activity in plants and human. J. Med. Plant Res. 2011, 5, 6697–6703. [Google Scholar] [CrossRef]
  57. Mandal, S.M.; Chakraborty, D.; Dey, S. Phenolic acids act as signaling molecules in plant-microbe symbioses. Plant Signal. Behav. 2010, 5, 359–368. [Google Scholar] [CrossRef] [Green Version]
  58. Gao, X.; Zhang, S.; Zhao, X.; Wu, Q. Potassium-induced plant resistance against soybean cyst nematode via root exudation of phenolic acids and plant pathogen-related genes. PLoS ONE 2018, 13, e0200903. [Google Scholar] [CrossRef]
  59. Cluis, C.P.; Mouchel, C.F.; Hardtke, C.S. The Arabidopsis transcription factor HY5 integrates light and hormone signaling pathways. Plant J. 2004, 38, 332–347. [Google Scholar] [CrossRef]
  60. Tohge, T.; Watanabe, M.; Hoefgen, R.; Fernie, A.R. The evolution of phenylpropanoid metabolism in the green lineage. Crit. Rev. Biochem. Mol. Biol. 2013, 48, 123–152. [Google Scholar] [CrossRef]
  61. Frohnmeyer, H.; Loyall, L.; Blatt, M.R.; Grabov, A. Millisecond UV-B irradiation evokes prolonged elevation of cytosolic-free Ca2+ and stimulates gene expression in transgenic parsley cell cultures. Plant J. 1999, 20, 109–117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Buer, C.S.; Muday, G.K.; Djordjevic, M.A. Flavonoids are differentially taken up and transported long distances in Arabidopsis. Plant Physiol. 2007, 145, 478–490. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Flint, S.D.; Caldwell, M.M. A biological spectral weighting function for ozone depletion research with higher plants. Physiol. Plant. 2003, 117, 137–144. [Google Scholar] [CrossRef]
  64. Senizza, B.; Zhang, L.; Miras-Moreno, B.; Righetti, L.; Zengin, G.; Ak, G.; Bruni, R.; Lucini, L.; Sifola, M.I.; El-Nakhel, C.; et al. The Strength of the Nutrient Solution Modulates the Functional Profile of Hydroponically Grown Lettuce in a Genotype-Dependent Manner. Foods 2020, 9, 1156. [Google Scholar] [CrossRef]
  65. Rothwell, J.A.; Pérez-Jiménez, J.; Neveu, V.; Medina-Ramon, A.; M’Hiri, N.; Garcia Lobato, P.; Manach, C.; Knox, K.; Eisner, R.; Wishart, D.; et al. Phenol-Explorer 3.0: A major update of the Phenol-Explorer database to incorporate data on the effects of food processing on polyphenol content. Database 2013, 2013, bat070. [Google Scholar] [CrossRef] [PubMed]
  66. Zhang, P.; Dreher, K.; Karthikeyan, A.; Chi, A.; Pujar, A.; Caspi, R.; Karp, P.; Kirkup, V.; Latendresse, M.; Lee, C.; et al. Creation of a Genome-Wide Metabolic Pathway Database for Populus trichocarpa Using a New Approach for Reconstruction and Curation of Metabolic Pathways for Plants. Plant Physiol. 2010, 153, 1479–1491. [Google Scholar] [CrossRef] [Green Version]
  67. Salek, R.M.; Steinbeck, C.; Viant, M.R.; Goodacre, R.; Dunn, W.B. The role of reporting standards for metabolite annotation and identification in metabolomic studies. Gigascience 2013, 2, 2047-217X. [Google Scholar] [CrossRef]
  68. Hellemans, J.; Mortier, G.; De Paepe, A.; Speleman, F.; Vandesompele, J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007, 8, 1–14. [Google Scholar] [CrossRef] [Green Version]
  69. Ghisoni, S.; Lucini, L.; Angilletta, F.; Rocchetti, G.; Farinelli, D.; Tombesi, S.; Trevisan, M. Discrimination of extra-virgin-olive oils from different cultivars and geographical origins by untargeted metabolomics. Food Res. J. 2019, 121, 746–753. [Google Scholar] [CrossRef]
  70. Caspi, R.; Dreher, K.; Karp, P.D. The Challenge of Constructing, Classifying, and Representing Metabolic Pathways. FEMS Microbiol. Lett. 2013, 345, 85–93. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Unsupervised hierarchical clustering analysis of the Micro-Tom tomato metabolome in leaves and roots clustered considering the treatment (control, CTR, or UV-treated, UV) and time (11 days of UV, 11 d, or 3 days after the end of the UV exposure, rec). A fold-change-based heat map was used to create dendrograms. Clustering and dendrograms were produced by choosing the Euclidean distance and Ward’s linkage rule.
Figure 1. Unsupervised hierarchical clustering analysis of the Micro-Tom tomato metabolome in leaves and roots clustered considering the treatment (control, CTR, or UV-treated, UV) and time (11 days of UV, 11 d, or 3 days after the end of the UV exposure, rec). A fold-change-based heat map was used to create dendrograms. Clustering and dendrograms were produced by choosing the Euclidean distance and Ward’s linkage rule.
Plants 11 01829 g001
Figure 2. OPLS-DA carried out from the UHPLC-ESI-QTOF metabolome of leaves and roots of Micro-Tom tomato plants, considering (a,c) the UV treatment factor; (b,d) UV treatment and time factors. Individual replications are given in the class prediction model score plot.
Figure 2. OPLS-DA carried out from the UHPLC-ESI-QTOF metabolome of leaves and roots of Micro-Tom tomato plants, considering (a,c) the UV treatment factor; (b,d) UV treatment and time factors. Individual replications are given in the class prediction model score plot.
Plants 11 01829 g002
Figure 3. Fold-change (FC) of metabolic processes in leaves and roots of UV treated plants, at 11 days of treatment and after 3 days of recovery, with respect to control plants. FC values were elaborated by the Omics Viewer Dashboard of the PlantCyc Pathway Tool software. Data for each biochemical pathway represent the sum of Log2FC of each detected compound. Syn, biosynthesis; Deg, degradation; FA/Lip, fatty acids and lipids; Carbo, carbohydrates; Sec Metab, secondary metabolites; AA, amino acids.
Figure 3. Fold-change (FC) of metabolic processes in leaves and roots of UV treated plants, at 11 days of treatment and after 3 days of recovery, with respect to control plants. FC values were elaborated by the Omics Viewer Dashboard of the PlantCyc Pathway Tool software. Data for each biochemical pathway represent the sum of Log2FC of each detected compound. Syn, biosynthesis; Deg, degradation; FA/Lip, fatty acids and lipids; Carbo, carbohydrates; Sec Metab, secondary metabolites; AA, amino acids.
Plants 11 01829 g003
Figure 4. OPLS-DA carried out from the UHPLC-ESI-QTOF phenolic profile of leaves and roots of Micro-Tom tomato plants, considering (a,c) UV treatment (b,d) UV treatment and time factors. Individual replications are given in the class prediction model score plot.
Figure 4. OPLS-DA carried out from the UHPLC-ESI-QTOF phenolic profile of leaves and roots of Micro-Tom tomato plants, considering (a,c) UV treatment (b,d) UV treatment and time factors. Individual replications are given in the class prediction model score plot.
Plants 11 01829 g004
Figure 5. Fold-change (FC) of phenolic accumulation in leaves and roots of UV treated plants at 11 days of treatment and after 3 days of recovery with respect to control plants. Compounds were grouped on the basis of phenolic class (a,b), flavonoid (c,d), and phenolic acids (e,f) subclasses. Data for each biochemical pathway represent the sum of the Log2FC of each detected compound.
Figure 5. Fold-change (FC) of phenolic accumulation in leaves and roots of UV treated plants at 11 days of treatment and after 3 days of recovery with respect to control plants. Compounds were grouped on the basis of phenolic class (a,b), flavonoid (c,d), and phenolic acids (e,f) subclasses. Data for each biochemical pathway represent the sum of the Log2FC of each detected compound.
Plants 11 01829 g005
Figure 6. Expression levels of genes involved and in the UVR8 pathway (HY5, COP1) and in phenolic biosynthesis (PAL, CHS, FLS) in (a,c,e,g,i) leaves and (b,d,f,h,j) roots of control (CTR) and UV-treated (UV) Micro-Tom tomato plants after 8 and 11 days of UV treatment (8 d and 11 d) and after 3 days from the end of the irradiation (rec). Data are mean ± SE of 3 biological replicates. According to each day, asterisks (*) indicate a significant difference between CTR and UV group (* p ≤ 0.05, ** p < 0.01, *** p < 0.001).
Figure 6. Expression levels of genes involved and in the UVR8 pathway (HY5, COP1) and in phenolic biosynthesis (PAL, CHS, FLS) in (a,c,e,g,i) leaves and (b,d,f,h,j) roots of control (CTR) and UV-treated (UV) Micro-Tom tomato plants after 8 and 11 days of UV treatment (8 d and 11 d) and after 3 days from the end of the irradiation (rec). Data are mean ± SE of 3 biological replicates. According to each day, asterisks (*) indicate a significant difference between CTR and UV group (* p ≤ 0.05, ** p < 0.01, *** p < 0.001).
Plants 11 01829 g006
Table 1. Leaf metabolic markers, according to their VIP score ± standard deviation (STD), related to the OPLS-DA model in Figure 2a (effect of UV treatment). Compounds are provided with p value and their Log2FC calculated with respect to control samples. * p ≤ 0.05, ** p < 0.01, *** p < 0.001, n.s., not significant.
Table 1. Leaf metabolic markers, according to their VIP score ± standard deviation (STD), related to the OPLS-DA model in Figure 2a (effect of UV treatment). Compounds are provided with p value and their Log2FC calculated with respect to control samples. * p ≤ 0.05, ** p < 0.01, *** p < 0.001, n.s., not significant.
CompoundVIP Score ± STDp ValueLog2 FC
2-oxoadipate2.23 ± 0.51***14.39
2,4-pyridinedicarboxylate/quinolinate2.09 ± 0.41***−0.58
rosmarinate2.07 ± 0.54n.s.0.34
desmethylxanthohumol2.01 ± 0.43n.s.−10.09
gibberellin A1.98 ± 0.40***−0.62
phlormethylbutanophenone1.96 ± 0.54n.s.4.00
O-sinapoylglucarolactone1.95 ± 0.89***−0.64
(-)-epipodophyllotoxin/β-peltatin/podophyllotoxin/5′-demethoxy-6-methoxypodophyllotoxin1.94 ± 0.55***0.52
L-α-amino-ε-keto-pimelate1.92 ± 0.54n.s.−8.37
ubiquinone-21.92 ± 0.50n.s.5.06
α-difluoromethylarginine1.91 ± 0.76n.s.−6.03
oxalosuccinate1.90 ± 0.86***12.36
1-O-(4-coumaroyl)-β-D-glucose1.88 ± 0.37n.s.13.85
(E)-4-hydroxy-3-methylbut-2-en-1-yl diphosphate1.87 ± 0.65***−0.31
3-deoxy-D-manno-octulosonate 8-phosphate1.86 ± 0.88***−0.57
pelargonidin/genistein/emodin/apigenin1.84 ± 0.54n.s.−10.19
4,5-seco-dopa1.83 ± 0.51n.s.−10.18
sn-glycero-3-phosphocholine1.83 ± 0.61n.s.13.75
N-methylanthranilate1.82 ± 0.36***−0.67
(+)-pisatin/1,7-dihydroxy-6,8-dimethoxy-2-methylanthraquinone/robustaquinone B/dimethylkaempferol/2′,7-dihydroxy-4′,5′-dimethoxyisoflavone/cirsimaritin/ladanein1.82 ± 0.83***0.51
solasodine1.80 ± 0.91**0.78
S-adenosyl-4-methylsulfanyl-2-oxobutanoate1.79 ± 0.76n.s.n.s.
2′-deoxymugineate1.78 ± 0.37*−0.26
2-hydroxy-4-carboxypyrimidine1.75 ± 0.77n.s.−10.39
3-C-glucosyl-2,4,4′,6-tetrahydroxydibenzoylmethane/6-C-glucosyl-2-hydroxynaringenin1.74 ± 0.70n.s.n.s.
A pentose1.72 ± 0.50***−0.72
acetonedicarboxylate/2-oxoglutarate1.72 ± 0.89*−0.46
5-methylcytidine/N4-aminocytidine1.71 ± 0.81n.s.14.58
2-deoxyglucose 6-phosphate/β-L-fucose 1-phosphate/1.69 ± 0.46n.s.−12.70
methylsuccinate/glutarate1.69 ± 0.98**−1.62
pyridoxal1.69 ± 0.70n.s.−10.19
(indol-3-yl)pyruvate1.69 ± 0.99**−0.62
tetraketide pyrone1.69 ± 0.74n.s.n.s.
capsidiol/lubimin1.68 ± 0.58*0.29
coniferaldehyde1.68 ± 0.37n.s.11.10
arabidopate1.67 ± 0.64n.s.0.34
dihydropapaverine/(S)-tetrahydrocolumbamine1.67 ± 0.64**−0.67
furcatin1.66 ± 0.55*−0.40
1-[18-hydroxyoeoyl]-2-[18-hydroxy-linoleoyl]-sn-glycerol1.66 ± 0.51*−9.12
1,2-dioleoylglycerol1.66 ± 0.54n.s.−9.53
(R)-S-lactoylglutathione1.65 ± 0.59n.s.−4.00
(3Z,6Z)-nonadienal1.65 ± 0.64n.s.−2.09
Fe(II)-nicotianamine1.64 ± 0.56*−0.77
(-)-4’-demethyl-deoxypodophyllotoxin/S-adenosyl-L-homocysteine1.64 ± 0.67**−0.59
quercetin-3-rhamnoside-7-rhamnoside/vicenin-2/isovitexin 2″-O-β-D-glucoside/vitexin 2″-O-β-D-glucoside/rhamnosylisoorientin1.64 ± 0.59*−0.48
A phenylpropanoid1.64 ± 0.71n.s.11.57
A phenylpropanoid conjugate1.63 ± 0.60*−0.48
L-dopachrome1.63 ± 0.97**1.97
protoheme1.63 ± 0.62**−1.21
pelargonidin 3-O-β-D-p-coumaroylglucoside1.62 ± 0.62*−0.47
D-tryptophan1.62 ± 0.73n.s.9.09
(4S)-2,3-dehydroleucopelargonidin/1,3,8-trihydroxy-2-methoxyanthraquinone/scutellarein/aureusidin1.61 ± 0.71*−0.43
cyanidin 3-O-β-D-p-coumaroylglucoside/cyanidin 3-(p-coumaroyl)-glucoside/pelargonidin 3-O-β-D-caffeoylglucoside1.60 ± 0.78*−1.30
N-vanillate-L-glutamate1.60 ± 0.71**−1.56
(2S)-eriodictyol/2-hydroxynaringenin/1-(4-hydroxyphenyl)-3-(2,4,6-trihydroxyphenyl)propane-1,3-dione/dalbergioidin1.60 ± 0.71*−0.43
7-(methylsulfanyl)heptyl-desulfoglucosinolate1.60 ± 0.71**0.49
1-aminocyclopropane-1-carboxylate/azetidine-2-carboxylate1.60 ± 0.77n.s.n.s.
(-)-isopiperitenone1.60 ± 0.56*2.08
(-)-yatein1.60 ± 0.73**0.47
1,2-di-O-sinapoyl-β-D-glucose/7-O-methylvitexin 2″-O-β-L-rhamnoside1.60 ± 0.61n.s.−11.11
TRIBOA-β-D-glucoside1.59 ± 0.85**−0.72
(4S)-2,3-dehydro-leucocyanidin/2-hydroxyeriodictyol/cis-12-sulfojasmonate1.58 ± 0.71*−0.42
N-hydroxyhomomethionine1.58 ± 1.19*−0.46
galactopinitol1.58 ± 0.76n.s.−8.41
scopolamine1.57 ± 0.49*9.38
cyanidin-3-O-rutinoside-5-O-β-D-glucoside1.57 ± 1.02*−0.68
(22R,23R)-28-homocastasterone1.57 ± 0.33*−8.88
acenaphthenequinone1.56 ± 0.91*−9.78
4-coumaroylhexanoylmethane1.56 ± 1.13n.s.−10.56
2-hydroxyferulate/2-hydroxycaffeate1.56 ± 0.52n.s.−9.25
hypoxanthine1.56 ± 0.35n.s.4.00
7,8-dihydromonapterin1.56 ± 0.35n.s.13.90
(R)-prunasin1.56 ± 0.35n.s.23.85
9-[6(S),9-diamino-5,6,7,8,9-pentadeoxy-β-D-ribo-nonafuranosyl]-9H-purin-6-amine1.55 ± 0.35*8.78
hydroxyechinenone1.55 ± 0.34n.s.21.57
luteoforol/(+)-catechin/leucopelargonidin1.54 ± 0.90n.s.−9.26
(5Z,8Z,11Z,14Z,17Z)-icosapentaenoate1.54 ± 0.47n.s.−7.60
(E)-2-(1H-indol-3-yl)-1-thioacetohydroximate1.54 ± 0.72n.s.n.s.
solavetivone/artemisinic aldehyde/germacra-1(10),4,11(13)-trien-12-al/(4S)-4-(5,5-dimethylcyclohex-1-en-1-yl)cyclohex-1-ene-1-carbaldehyde/zerumbone1.53 ± 0.54n.s.0.24
cysteine1.53 ± 0.99**−1.27
olivetolate1.53 ± 1.06n.s.n.s.
17-O-acetylajmaline1.53 ± 0.43n.s0.24
phenylmethanesulfenate1.53 ± 1.00n.s−0.24
3-carboxy-8-(methylsulfanyl)-2-oxooctanoate1.53 ± 0.99n.s.−8.76
(E)-phenylacetaldehyde oxime/2-phenylacetamide/N-benzylformamide/(Z)-phenylacetaldehyde oxime1.52 ± 0.54n.s.−0.35
solasodine 3-O-β-D-glucoside1.52 ± 0.66*−2.57
6,8-dihydroxypurine1.51 ± 0.53n.s.12.65
D,L-α-methylphosphinothricin1.50 ± 0.49n.s.−4.12
5-[[4-methoxy-3-(phenylmethoxy)phenyl]methyl]-2,4-pyrimidinediamine1.50 ± 1.04*0.88
cyclo-dopa 5-O-glucoside1.50 ± 0.52**−9.12
(E)-8-(methylsulfanyl)octanal oxime1.50 ± 0.52n.s.2.33
6-amino-2-oxohexanoate1.50 ± 0.58n.s.−9.18
adenosine 3′,5′-bisphosphate/GDP group1.50 ± 0.53n.s.n.s.
avenastenone1.48 ± 0.76n.s.n.s.
3-chlorodiaminopimelate1.48 ± 0.83*−0.36
3-hydroxylubimin1.48 ± 0.92n.s.n.s.
(S)-4-hydroxymandelonitrile/(R)-4-hydroxymandelonitrile/2-formylaminobenzaldehyde/3-hydroxyindolin-2-one/5,6-dihydroxyindole/p-hydroxymandelonitrile/1.48 ± 0.87n.s.−7.84
dihydroconiferyl alcohol/iridotrial1.48 ± 0.58*0.75
1-18:3-2-18:3-digalactosyldiacylglycerol1.47 ± 0.48n.s.−0.17
Se-methyl-Se-L-methionine1.46 ± 0.70n.s.−6.46
(12Z,15Z)-9,10-dihydroxyoctadeca-12,15-dienoate/HPODE/(9Z,12Z)-15,16-dihydroxyoctadeca-9,12-dienoate/(9Z,15Z)-12,13-dihydroxyoctadeca-9,15-dienoate/9,10-12,13-diepoxyoctadecanoate1.46 ± 0.53n.s.−3.01
benzyl alcohol 6-O-β-D-xylopyranosyl-β-D-glucopyranoside1.46 ± 0.60n.s.−3.59
4-hydroxy-5-methyl-2-propyl-3(2H)-furanone1.46 ± 0.63n.s.6.86
S-adenosyl 3-(methylsulfanyl)propylamine1.45 ± 0.22n.s.−5.26
gibberellin A251.45 ± 1.23n.s.n.s.
2-hydroxysuccinamate/aspartate1.45 ± 0.48n.s.n.s.
1-18:3-2-16:3-monogalactosyldiacylglycerol1.44 ± 0.39n.s.−0.17
3-hydroxy-9-apo-δ-caroten-9-one1.44 ± 0.52*8.23
(S)-2-amino-6-oxohexanoate1.43 ± 0.28n.s.−0.10
(3S,5R,6R)-3,5-dihydroxy-6,7-didehydro-5,6-dihydro-12′-apo-β-caroten-12′-al1.43 ± 0.65*0.45
tricoumaroyl spermidine1.43 ± 0.81*7.62
CMP group1.43 ± 0.63n.s.−6.18
2-C-methyl-D-erythritol-2,4-cyclodiphosphate1.43 ± 0.69n.s.n.s.
D-galactosylononitol1.42 ± 1.02n.s.−9.05
(indol-3-yl)acetyl-L-glutamine1.42 ± 0.58*−2.57
2-carboxy-1,4-naphthoquinone1.42 ± 1.06*−0.43
artemisinin1.42 ± 0.56n.s.12.07
enol-oxaloacetate/oxaloacetate1.42 ± 0.99n.s.−0.37
hemigossypol-6-methyl ether1.42 ± 0.78n.s.n.s.
A terpenoid1.41 ± 0.65n.s.6.49
chelirubine1.41 ± 0.75n.s.−2.27
tartrate1.41 ± 1.03n.s.−0.37
palmitoleate1.41 ± 0.63n.s.6.00
1-heptanal1.40 ± 0.76n.s.−6.31
16,17-dihydro-16α,17-dihydroxy gibberellin A41.40 ± 0.56n.s.4.88
Table 2. Root metabolic markers, according to their VIP score ± standard deviation (STD), related to the OPLS-DA model in Figure 2c (effect of UV treatment). Compounds are provided with p value and their Log2FC calculated in respect to control samples. * p ≤ 0.05, ** p < 0.01, *** p < 0.001, n.s., not significant.
Table 2. Root metabolic markers, according to their VIP score ± standard deviation (STD), related to the OPLS-DA model in Figure 2c (effect of UV treatment). Compounds are provided with p value and their Log2FC calculated in respect to control samples. * p ≤ 0.05, ** p < 0.01, *** p < 0.001, n.s., not significant.
CompoundVIP Score ± STDp ValueLog2 FC
vanillylamine2.49 ± 0.60n.s.−0.08
10-deacetyl-2-debenzoylbaccatin III2.41 ± 0.26***12.13
A pentose2.40 ± 0.56***13.43
gardenin B/dalnigrein/3,7,3′,4′-tetramethylquercetin2.39 ± 0.60n.s.4.00
heliocide B12.28 ± 0.31***2.46
N-methylanthranilate2.16 ± 0.71***12.27
L-α-(methylenecyclopropyl)-glycine2.10 ± 0.83***12.89
5,10-methylenetetrahydropteroyl mono-L-glutamate2.10 ± 0.58***14.55
1,2-dehydroreticuline/(S)-corytuberine2.06 ± 0.32n.s.11.07
1-hydroxycumene2.03 ± 0.62n.s.0.98
cyclogutamate2.01 ± 0.41n.s.14.83
all-trans-hexaprenyl diphosphate1.96 ± 1.19***13.26
syn-copalyl diphosphate1.93 ±0.56n.s.n.s.
salutaridine1.90 ± 0.69n.s.3.66
quercetin 3-O-(4\-O-p-coumaroyl)-glucoside″1.90 ± 0.85n.s.n.s.
geranate/(1R)-trans-chrysanthemate/N-hydroxy-L-valine1.86 ± 0.35*0.56
(+)-carvone/(-)-isopiperitenone/(-)-carvone/4-isopropylbenzyl alcohol/(4S)-perillyl aldehyde/(+)-menthofuran1.86 ± 0.37*0.56
(+)-bornane-2,5-dione1.86 ± 0.46n.s.2.77
O-sinapoylglucarolactone1.83 ± 1.15*0.49
butyl propanoate1.82 ± 0.52n.s.−9.14
quercetin 3-O-(-O-p-coumaroyl)-glucoside1.82 ± 0.90n.s.14.77
2,4-diamino-6-methyl-5,3′-(3-nitrophenoxy)prop-1′-yloxypyrimidine1.78 ± 1.15n.s.n.s.
robustaquinone F1.77 ± 1.06**0.76
(3-hydroxy-2-oxindol-3-yl)acetyl-L-aspartate1.77 ± 1.07**0.77
cyanidin 3-O-β-D-caffeoylglucoside1.76 ± 1.04**0.77
kaempferol 3-O-β-D-glucosylgalactoside/kaempferol 3-O-β-D-glucosyl-(1->2)-β-D-glucoside/1.76 ± 1.04**0.77
cyanidin 3-O-β-D-p-coumaroylglucoside/cyanidin 3-(p-coumaroyl)-glucoside/pelargonidin 3-O-β-D-caffeoylglucoside1.76 ± 1.05**0.79
4-(β-D-glucosyloxy)benzoate1.76 ± 0.85**10.55
α-hydroxyheme1.76 ± 0.98n.s.2.90
scopolamine1.73 ± 0.64n.s.−2.15
5-oxooctanal1.72 ± 0.94*−0.47
22-oxo-docosanoyl-CoA/3-oxobehenoyl-CoA1.71 ± 0.67n.s.−9.97
L-nicotianamine1.70 ± 1.00***15.39
soyasaponin I1.68 ± 0.54n.s.−0.40
epoxypheophorbide a1.68 ± 0.94n.s.9.89
heliannuol C1.66 ± 1.09n.s.−3.30
1-α-linolenoyl-2-α-linolenoyl-phosphatidylcholine1.66 ± 0.81n.s.−0.54
presqualene diphosphate1.65 ± 1.05n.s.12.64
iridotrial/dihydroconiferyl alcohol1.65 ± 0.67n.s.2.67
ricinine1.65 ± 0.67n.s.2.67
conhydrinone/hygrine/1.64 ± 0.58n.s.−1.92
5-hydroxy-γ-coniceine1.61 ± 0.58n.s.−3.51
2′-hydroxypseudobaptigenin1.61 ± 1.00n.s.−0.10
rhamnetin1.60 ± 1.01n.s.−0.11
grasshopper ketone1.59 ± 1.25n.s.−0.58
(E)-1-(L-cysteinylglycin-S-yl)-N-hydroxy-ω-(methylsulfanyl)pentan-1-imine1.58 ± 0.51n.s.8.29
hydroxy-β-zeacarotene1.58 ± 1.32n.s.−6.89
β-D-glucosyl crocetin1.58 ± 0.49n.s.9.40
N-acetyl-α-D-galactosamine 1-phosphate/N-acetyl-α-D-glucosamine 1-phosphate1.58 ± 0.89n.s.−0.08
curcumin 4′-O-β-D-gentiotrioside/curcumin 4′-O-β-D-gentiobiosyl 4″-O-β-D-glucoside1.57 ± 1.25n.s.7.27
1-palmitoyl-2-vernoloyl-phosphatidylcholine1.56 ± 0.70n.s.−1.99
4α-carboxy-5α-cholesta-7,24-dien-3β-ol1.56 ± 0.74*−1.24
1-18:1-2-16:0-monogalactosyldiacylglycerol1.55 ± 0.86n.s.−0.39
munjistin1.55 ± 1.14n.s.−0.14
salicylate 2-O-β-D-glucoside1.55 ± 0.95*−1.75
glutathione disulfide1.54 ± 1.06n.s.1.51
1-18:0-2-18:3-phosphatidylethanolamine1.54 ± 1.05n.s.−1.67
1-18:3-2-16:0-digalactosyldiacylglycerol/1-16:0-2-18:3-digalactosyldiacylglycerol1.54 ± 1.16n.s.−0.90
N-(4-aminobenzoyl)-L-glutamate1.54 ± 1.03*8.41
aloesone1.54 ± 0.49n.s.13.60
A phenylpropanoid1.53 ± 0.98n.s.−0.08
1-18:2-2-18:2-monogalactosyldiacylglycerol1.52 ± 0.99*1.79
ginsenoside Rg11.52 ± 0.26n.s.−1.74
7,8-diaminopelargonate1.52 ± 0.85n.s.−2.56
3β-hydroxy-β-cyclocitral/(6E)-8-hydroxygeranial/(6E)-8-oxogeraniol1.52 ± 0.71*1.22
2,6-diaminopurine/heptanoate1.51 ± 0.37n.s.−0.19
stevioside/rebaudioside B1.51 ± 0.88n.s.−7.24
ajmaline1.51 ± 0.94n.s.−8.37
pheophorbide b1.51 ± 0.86n.s.−6.15
1-palmitoyl-2-linoleoyl-phosphatidylcholine1.50 ± 0.94n.s.−0.31
22-hydroxydocosanoate1.49 ± 1.10n.s.−0.39
(E)-1-(L-cysteinylglycin-S-yl)-N-hydroxy-ω-(methylsulfanyl)heptan-1-imine/sinapoyl-(S)-malate1.47 ± 0.30n.s.−6.58
eupatolin1.47 ± 1.01n.s.0.42
deoxyhumulone1.47 ± 1.29n.s.−0.54
(S)-scoulerine/(S)-coreximine1.47 ± 0.67n.s.0.28
phosphopantetheine1.47 ± 0.30n.s.−6.63
delphinidin 3-O-rutinoside-7-O-glucoside/quercetin 3-O-rhamnosyl(1->2)glucoside-7-O-rhamnoside/quercetin 3-O-gentiobioside-7-O-rhamnoside/kaempferol 3-O-β-D-glucosyl-(1->2)-glucosyl-(1->2)-β-D-glucoside1.47 ± 0.56n.s.4.44
cycloheptadienyl/sinapate1.47 ± 0.96**−10.73
4α-carboxy-5α-cholesta-8,24-dien-3β-ol1.46 ± 0.78n.s.−3.30
1-deoxy-2,3-hexodiulose-6-phosphate1.46 ± 1.02n.s.−0.07
A saccharide1.46 ± 0.98n.s.−0.07
(indol-3-yl)acetyl-L-proline1.46 ± 0.57n.s.−6.48
catharanthine1.46 ± 0.95n.s.−0.44
geranyl β-D-glucopyranoside1.46 ± 0.89*0.63
ajmaline-x1.46 ± 1.17n.s.n.s.
6-methoxypodophyllotoxin1.45 ± 1.14n.s.3.41
p-coumaroylserotonin1.45 ± 0.49n.s.8.25
(2E,4E,6E)-4-methylocta-2,4,6-trienedial1.45 ± 0.39n.s.−0.18
(-)-medicarpin-3-O-glucoside1.45 ± 0.26n.s.5.54
3-oxocerotoyl-CoA1.45 ± 1.05n.s.−7.81
1,3,5-trimethoxybenzene1.45 ± 0.39n.s.−0.18
1-linoleoyl-2-palmitoyl-phosphatidylglycerol1.45 ± 0.90*−1.29
1-stearoyl-sn-glycerol 3-phosphate1.44 ± 0.64n.s.−0.17
β-citraurin1.44 ± 0.62n.s.−0.19
myristate1.44 ± 0.95n.s.−0.32
phosphocholine1.43 ± 1.16n.s.−5.00
tetramethylmyricetin/tetramethylquercetagetin1.43 ± 0.31n.s.−5.43
norbelladine1.43 ± 0.98n.s.7.39
cyanidin O-O-[6-O-(6-O-4-hydroxycinnamoyl-β-D-glucosyl)-2-O-β-D-xylosyl-β-D-galactoside]1.43 ± 1.49n.s.0.48
cytochrome c1.43 ± 0.64n.s.n.s.
gibberellin A51-catabolite1.42 ± 0.52n.s.0.40
26-hydroxybrassinolide1.42 ± 1.22n.s.−0.47
2-O-caffeoylglucarate1.42 ± 0.61n.s.−7.41
allosamidin/N,N′,N″-triacetylchitotriose1.41 ± 0.73n.s.−2.50
5-phospho-D-arabinonohydroxamate1.41 ± 1.07n.s.−0.06
1,2-dipalmitoyl-phosphatidylglycerol-phosphate1.41 ± 1.26n.s.−1.36
DIMBOA-β-D-glucoside1.41 ± 0.36n.s.−0.76
coumarinate1.41 ± 0.40n.s.−0.50
(R)-3-(4-hydroxyphenyl)lactate1.40 ± 0.39n.s.−0.99
scopolin1.40 ± 0.38n.s.−0.72
protochlorophyllide a1.40 ± 0.84n.s.−1.85
(E)-1-(L-cysteinylglycin-S-yl)-N-hydroxy-2-(1H-indol-3-yl)ethan-1-imine1.40 ± 0.39n.s.−0.67
chlorophyllide b1.40 ± 1.29n.s.0.40
anthraniloyl-O-glucopyranose1.40 ± 0.51n.s.−0.97
5-(methylsulfanyl)pentyl-desulfoglucosinolate1.40 ± 0.31n.s.−1.30
gibberellin A12-aldehyde1.40 ± 0.52n.s.8.42
isatin1.39 ± 1.43n.s.−0.05
tetradecan-1-ol1.39 ± 0.55n.s.0.00
1-18:3-2-16:4-monogalactosyldiacylglycerol1.38 ± 0.98n.s.−2.10
(indol-3-yl)acetyl-myo-inositol L-arabinoside1.38 ± 0.93n.s.0.36
cis-tuberonic acid1.38 ± 1.31n.s.−6.50
(1,4)-β-xylobiose/7,8-dimethoxy-flavone/3-methylinosine1.38 ± 0.92n.s.7.30
p-aminobenzoate-β-D-glucopyranosyl ester1.38 ± 0.48n.s.−7.95
1-(3,4-dihydroxyphenyl)-5-hydroxy-3-decanone/3″-hydroxy-geranylhydroquinone1.38 ± 1.50n.s.−0.74
phenylarsine oxide1.37 ± 1.02n.s.−0.05
tyramine1.37 ± 0.63n.s.−0.16
trans-tuberonic acid1.37 ± 1.31n.s.−1.18
CDP-1-18:1(9Z)-2-18:1(9Z)-glycerol1.36 ± 0.89*0.22
phlorisobutanophenone1.36 ± 1.34n.s.5.51
16-epivellosimine1.36 ± 1.23n.s.0.00
5-mercaptodeoxyuridine1.36 ± 0.77n.s.−0.08
indole-3-carbinol1.35 ± 1.06n.s.−6.80
capsanthin/antheraxanthin/4-hydroxyzeaxanthin1.35 ± 0.85n.s.−0.64
3,4-dihydrocoumarin/4-coumaraldehyde/1-phenylpropane-1,2-dione1.35 ± 0.68n.s.−0.15
dimethylsulfoniopropanoate-amine1.34 ± 0.96n.s.−1.39
kaempferide 3-O-β-D-glucopyranosyl-(1->2)-O-α-L-rhamnoside1.34 ± 1.49n.s.0.41
5-methyl-DL-tryptophan1.34 ± 0.51n.s.n.s.
3,8-divinyl protochlorophyllide a1.33 ± 1.50n.s.0.409
1-18:1-2-16:0-phosphatidate1.33 ± 0.91n.s.−7.65
N-demethylvindolidine1.32 ± 1.20n.s.−2.73
dihydrogeranylgeranyl diphosphate1.32 ± 0.38n.s.n.s.
gibberellin A91.32 ± 0.99n.s.−0.29
18-hydroxyoleate/9,10-epoxystearate/ricinoleate1.32 ± 0.73n.s.−0.36
(7Z,10Z,13Z,16Z,19Z)-docosa-7,10,13,16,19-pentenoate1.32 ± 0.60n.s.−0.47
pheophytin b1.32 ± 1.08n.s.−7.60
3-oxo-2-(cis-2′-pentenyl)-cyclopentane-1-butanoyl-CoA1.32 ± 1.51n.s.0.29
kaempferol 3-O-(\-O-p-coumaroyl)-glucoside″/kaempferol-3-O-β-D-glucopyranoside-7-O-α-L-rhamnopyranoside1.32 ± 1.41n.s.n.s.
linoleate1.32 ± 0.71n.s.−0.36
α,ω-9Z-octadecenedioate1.32 ± 0.68n.s.−0.51
prenyl diphosphate/isopentenyl diphosphate1.31 ± 1.00n.s.−0.05
(+)-jasmonate1.31 ± 0.91n.s.4.85
6-isobutyl-4-hydroxy-2-pyrone1.31 ± 0.84n.s.−2.59
α-chaconine1.31 ± 1.14n.s.−7.90
HPOTE/(9Z,11E,14Z)-(13S)-hydroperoxyoctadeca-(9,11,14)-trienoate1.31 ± 0.71n.s.−6.35
palmitoleoyl-CoA1.31 ± 1.60n.s.0.25
pheophytin a1.31 ± 1.14n.s.−6.55
imidazole acetol-phosphate1.31 ± 1.01n.s.−0.05
N-hydroxypentahomomethionine1.30 ± 0.83n.s.−0.75
all-trans-nonaprenyl diphosphate1.30 ± 0.99n.s.0.08
kauralexin B31.30 ± 0.70n.s.0.72
Table 3. Discriminant phenolic compounds identified in leaves according to VIP associated with the OPLS-DA model reported in Figure 4a (effect of UV treatment). Phenolic compounds, classified according to their class and subclass, are provided with p value and their Log2FC calculated in respect to control samples. ** p < 0.01, *** p < 0.001, n.s., not significant.
Table 3. Discriminant phenolic compounds identified in leaves according to VIP associated with the OPLS-DA model reported in Figure 4a (effect of UV treatment). Phenolic compounds, classified according to their class and subclass, are provided with p value and their Log2FC calculated in respect to control samples. ** p < 0.01, *** p < 0.001, n.s., not significant.
CompoundClassSubclassVIP Score ± STDp ValueLog2 FC
Cyanidin 3-O-(6″-p-coumaroyl-glucoside)/Petunidin 3-O-rutinoside/Pelargonidin 3-O-sophoroside/Cyanidin 3-O-rutinosideFlavonoidsAnthocyanins1.40 ± 0.42n.s.−0.26
CyanidinFlavonoidsAnthocyanins1.38 ± 0.49n.s.−0.20
Cyanidin 3-O-glucosyl-rutinosideFlavonoidsAnthocyanins1.35 ± 0.64n.s.−0.49
Cyanidin 3-O-glucoside/Cyanidin 3-O-galactoside/Petunidin 3-O-arabinoside/Peonidin 3-O-arabinoside/p-Coumaric acidFlavonoidsAnthocyanins1.24 ± 0.90n.s.−6.47
Delphinidin 3-O-glucoside/Delphinidin 3-O-galactosideFlavonoidsAnthocyanins1.13 ± 0.40n.s.−6.09
DihydroquercetinFlavonoidsDihydroflavonols1.38 ± 0.50n.s.−0.20
(-)-Epicatechin 3-O-gallate/(+)-Catechin 3-O-gallateFlavonoidsFlavanols1.22 ± 0.97n.s.−6.54
Eriocitrin/NeoeriocitrinFlavonoidsFlavanones1.38 ± 0.48n.s.−0.26
CirsimaritinFlavonoidsFlavones1.59 ± 0.76***0.72
Apigenin 7-O-(6″-malonyl-apiosyl-glucoside)FlavonoidsFlavones1.13 ± 0.39n.s.−0.54
Kaempferol 7-O-glucosideFlavonoidsFlavonols1.48 ± 0.72n.s.−2.51
Kaempferol 3-O-rutinoside/Chrysoeriol 7-O-apiosyl-glucoside/Apigenin 6.8-di-C-glucoside/Kaempferol 3-O-galactoside 7-O-rhamnoside/Luteolin 7-O-rutinosideFlavonoidsFlavonols1.41 ± 0.40n.s.−0.26
Kaempferol/Scutellarein/LuteolinFlavonoidsFlavonols1.38 ± 0.50n.s.−0.19
Kaempferol 3-O-glucosyl-rhamnosyl-galactoside/Quercetin 3-O-rhamnosyl-rhamnosyl-glucoside/Kaempferol 3-O-glucosyl-rhamnosyl-glucosideFlavonoidsFlavonols1.35 ± 0.64n.s.−0.49
3-Methoxysinensetin/NobiletinFlavonoidsFlavonols1.14 ± 0.61n.s.−1.65
3-Methylcatechol/Guaiacol/4-MethylcatecholOther polyphenolsAlkylphenols1.79 ± 0.34n.s.4.94
5-NonadecenylresorcinolOther polyphenolsAlkylphenols1.25 ± 0.75n.s.7.75
5-HeneicosenylresorcinolOther polyphenolsAlkylphenols1.10 ± 0.81n.s.−0.56
Umbelliferone/4-HydroxycoumarinOther polyphenolsHydroxycoumarins1.10 ± 1.00n.s.0.10
p-HPEA-EA/Ligstroside-aglyconeOther polyphenolsTyrosols1.53 ± 0.58n.s.−0.36
LigstrosideOther polyphenolsTyrosols1.14 ± 0.52n.s.0.10
Galloyl glucosePhenolic acidsHydroxybenzoic acids1.21 ± 1.28n.s.1.78
Gallic acid 4-O-glucosidePhenolic acidsHydroxybenzoic acids1.18 ± 0.82n.s.0.08
Rosmarinic acidPhenolic acidsHydroxycinnamic acids1.69 ± 0.64***−4.00
Caffeoyl aspartic acidPhenolic acidsHydroxycinnamic acids1.52 ± 0.72**0.52
p-Coumaroyl tartaric acidPhenolic acidsHydroxycinnamic acids1.44 ± 0.89n.s.1.85
4-Feruloylquinic acid/3-Feruloylquinic acid/5-Feruloylquinic acidPhenolic acidsHydroxycinnamic acids1.38 ± 0.77n.s.−3.71
Ferulic acid 4-O-glucoside/Feruloyl glucosePhenolic acidsHydroxycinnamic acids1.32 ± 0.65n.s.−0.44
Isoferulic acid/Ferulic acidPhenolic acidsHydroxycinnamic acids1.30 ± 0.62n.s.−0.54
Caffeoyl tartaric acidPhenolic acidsHydroxycinnamic acids1.28 ± 1.15n.s.−6.84
m-Coumaric acid/o-Coumaric acidPhenolic acidsHydroxycinnamic acids1.24 ± 1.25n.s.−8.72
Avenanthramide 2c/Avenanthramide KPhenolic acidsHydroxycinnamic acids1.18 ± 0.91n.s.0.82
Dihydrocaffeic acid/Syringaldehyde/Homovanillic acidPhenolic acidsHydroxyphenylpropanoic acids1.20 ± 1.18n.s.−7.57
Table 4. Discriminant phenolic compounds identified in roots according to VIP associated with the OPLS-DA model reported in Figure 4c (effect of UV treatment). Phenolic compounds, classified according to their class and subclass, are provided with p value and their Log2FC calculated in respect to control samples. * p ≤ 0.05, n.s., not significant.
Table 4. Discriminant phenolic compounds identified in roots according to VIP associated with the OPLS-DA model reported in Figure 4c (effect of UV treatment). Phenolic compounds, classified according to their class and subclass, are provided with p value and their Log2FC calculated in respect to control samples. * p ≤ 0.05, n.s., not significant.
CompoundClassSubclassVIP Score ± STDp ValueLog2 FC
Delphinidin 3,5-O-diglucoside/Delphinidin 3-O-glucosyl-glucosideFlavonoidsAnthocyanins1.70 ± 0.59n.s.0.60
Delphinidin 3-O-(6″-p-coumaroyl-glucoside)FlavonoidsAnthocyanins1.70 ± 0.62n.s.0.64
Cyanidin 3-O-(6″-caffeoyl-glucoside)FlavonoidsAnthocyanins1.70 ± 0.62n.s.0.64
Pelargonidin 3,5-O-diglucosideFlavonoidsAnthocyanins1.68 ± 0.48n.s.8.60
Pigment A/Peonidin 3-O-(6″-p-coumaroyl-glucoside)FlavonoidsAnthocyanins1.51 ± 0.63n.s.−0.58
Peonidin 3-O-rutinosideFlavonoidsAnthocyanins1.17 ± 1.05n.s.−7.18
Prodelphinidin dimer B3FlavonoidsFlavanols1.70 ± 0.60n.s.0.61
Gardenin BFlavonoidsFlavones2.19 ± 0.86n.s.4.00
Nepetin/Isorhamnetin/RhamnetinFlavonoidsFlavones1.41 ± 0.75n.s.−0.28
Apigenin 7-O-(6″-malonyl-apiosyl-glucoside)FlavonoidsFlavones1.31 ± 0.79n.s.−0.99
Quercetin 3-O-galactoside 7-O-rhamnoside/Kaempferol 3-O-sophoroside/Quercetin 3-O-rutinoside/Kaempferol 3,7-O-diglucoside/Quercetin 3-O-rhamnosyl-galactosideFlavonoidsFlavonols1.72 ± 0.56n.s.3.03
Quercetin/Morin/p-Coumaroyl malic acid/6-HydroxyluteolinFlavonoidsFlavonols1.70 ± 0.62n.s.0.58
4-Vinylguaiacol/3-MethoxyacetophenoneOther polyphenolsAlkylmethoxyphenols1.40 ± 0.45n.s.−0.36
p-AnisaldehydeOther polyphenolsHydroxybenzaldehydes1.14 ± 1.06n.s.−0.83
Thymol/CarvacrolOther polyphenolsPhenolic terpenes1.67 ± 0.56n.s.0.39
Hydroxytyrosol 4-O-glucosideOther polyphenolsTyrosols1.24 ± 0.73n.s.−1.00
3,4-DHPEA-ACOther polyphenolsTyrosols1.10 ± 1.09n.s.5.34
4-Hydroxybenzoic acid 4-O-glucosidePhenolic acidsHydroxybenzoic acids1.16 ± 0.53n.s.0.018
Sinapic acidPhenolic acidsHydroxycinnamic acids1.54 ± 1.26n.s.8.77
p-Coumaroyl glycolic acidPhenolic acidsHydroxycinnamic acids1.44 ± 0.63*−10.91
5-Caffeoylquinic acid/3-Caffeoylquinic acid/4-Caffeoylquinic acidPhenolic acidsHydroxycinnamic acids1.25 ± 0.19n.s.9.47
p-Coumaric acidPhenolic acidsHydroxycinnamic acids1.21 ± 0.25n.s.5.25
Caffeic acidPhenolic acidsHydroxycinnamic acids1.11 ± 0.98n.s.5.34
Homoveratric acidPhenolic acidsHydroxyphenylacetic acids1.10 ± 1.09n.s.−0.36
d-Viniferin/e-Viniferin/PallidolStilbenesStilbenes1.19 ± 0.28n.s.0.57
Table 5. List of the phenolic compounds in UV-treated leaves, selected by combining analysis of variance and fold-change into volcano plot (Benjamini–Hochberg multiple testing correction, p < 0.05; fold-change cut-off = 1.2). The symbols “↑” and “↓” mean increased and decreased accumulation, respectively.
Table 5. List of the phenolic compounds in UV-treated leaves, selected by combining analysis of variance and fold-change into volcano plot (Benjamini–Hochberg multiple testing correction, p < 0.05; fold-change cut-off = 1.2). The symbols “↑” and “↓” mean increased and decreased accumulation, respectively.
SamplePhenolic ClassPhenolic SubclassMarkerLog2 FCAccumulation
UV-11dOther polyphenolsTyrosolsOleoside 11-methylester17.42
Phenolic acidsHydroxycinnamic acidsRosmarinic acid−19.09
UV-recFlavonoidsFlavonesCirsimaritin17.68
Flavonols6,8-Dihydroxykaempferol/Myricetin13.85
Flavonols/FlavonolsNepetin/Isorhamnetin/Rhamnetin16.44
FlavonesTetramethylscutellarein−15.90
Lignans-Dimethylmatairesinol19.12
Lignans/Hydroxybenzoic acidsProtocatechuic aldehyde−18.90
Lignans/Hydroxybenzoic acidsSesamol/4-Hydroxybenzoic acid/2-Hydroxybenzoic acid−18.92
Phenolic acidsHydroxybenzoic acids3-Hydroxybenzoic acid−18.90
Gallic acid 4-O-glucoside−8.90
Galloyl glucose−8.90
Hydroxycinnamic acidsCaffeoyl aspartic acid18.02
Ferulic acid 4-O-glucoside/Feruloyl glucose−17.44
m-Coumaric acid/o-Coumaric acid−15.02
p-Coumaric acid 15.95
p-Coumaroyl tartaric acid18.59
Rosmarinic acid−18.83
Hydroxyphenylpropanoic acidsDihydro-p-coumaric acid/Methoxyphenylacetic acid−12.05
Other polyphenolsAlkylmethoxyphenols/Hydroxybenzoketones4-Vinylguaiacol/3-Methoxyacetophenone12.67
Alkylphenols4-Vinylphenol20.91
Hydroxybenzaldehydesp-Anisaldehyde−12.43
HydroxycoumarinsCoumarin14.96
Tyrosol3,4-DHPEA-EDA−15.68
Tyrosols/AlkylphenolsTyrosol/4-Ethylcatechol20.92
Table 6. List of the phenolic compounds in UV-treated roots, selected by combining analysis of variance and fold-change into volcano plot (Benjamini–Hochberg multiple testing correction, p < 0.05; fold-change cut-off = 1.2). The symbols “↑” and “↓” mean increased and decreased accumulation, respectively.
Table 6. List of the phenolic compounds in UV-treated roots, selected by combining analysis of variance and fold-change into volcano plot (Benjamini–Hochberg multiple testing correction, p < 0.05; fold-change cut-off = 1.2). The symbols “↑” and “↓” mean increased and decreased accumulation, respectively.
SamplePhenolic ClassPhenolic SubclassMarkerLog2 FCAccumulation
UV-11dFlavonoidsFlavonesGardenin B20.62
UV-recFlavonoidsFlavonesGardenin B21.69
Jaceidin 4′-O-glucuronide−15.19
Phenolic acidsHydroxycinnamic acids5-Caffeoylquinic acid/3-Caffeoylquinic acid/4-Caffeoylquinic acid−21.25
p-Coumaroyl glycolic acid−13.02
p-Coumaric acid−19.58
Sinapic acid4.45
Lignans-Pinoresinol/Matairesinol−8.43
Phenolic acids/Other polyphenols/Phenolic acidsHydroxyphenylpropanoic acids/Hydroxybenzaldehydes/Hydroxyphenylacetic acidsDihydrocaffeic acid/Syringaldehyde/Homovanillic acid16.32
Other polyphenolsAlkylmethoxyphenols/Hydroxybenzoketones4-Vinylguaiacol/3-Methoxyacetophenone−19.39
Alkylphenols5-Pentadecylresorcinol−18.23
4-Vinylphenol−21.38
HydroxyphenylpropenesAcetyl eugenol−19.07
Phenolic terpenesCarnosic acid−17.60
StilbeneResveratrol−18.66
TyrosolsHydroxytyrosol 4-O-glucoside19.83
Tyrosols/AlkylphenolsTyrosol/4-Ethylcatechol−21.37
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mannucci, A.; Santin, M.; Vanhaelewyn, L.; Sciampagna, M.C.; Miras-Moreno, M.B.; Zhang, L.; Lucini, L.; Quartacci, M.F.; Van Der Straeten, D.; Castagna, A.; et al. Foliar and Root Comparative Metabolomics and Phenolic Profiling of Micro-Tom Tomato (Solanum lycopersicum L.) Plants Associated with a Gene Expression Analysis in Response to Short Daily UV Treatments. Plants 2022, 11, 1829. https://doi.org/10.3390/plants11141829

AMA Style

Mannucci A, Santin M, Vanhaelewyn L, Sciampagna MC, Miras-Moreno MB, Zhang L, Lucini L, Quartacci MF, Van Der Straeten D, Castagna A, et al. Foliar and Root Comparative Metabolomics and Phenolic Profiling of Micro-Tom Tomato (Solanum lycopersicum L.) Plants Associated with a Gene Expression Analysis in Response to Short Daily UV Treatments. Plants. 2022; 11(14):1829. https://doi.org/10.3390/plants11141829

Chicago/Turabian Style

Mannucci, Alessia, Marco Santin, Lucas Vanhaelewyn, Maria Calogera Sciampagna, Maria Begoña Miras-Moreno, Leilei Zhang, Luigi Lucini, Mike Frank Quartacci, Dominique Van Der Straeten, Antonella Castagna, and et al. 2022. "Foliar and Root Comparative Metabolomics and Phenolic Profiling of Micro-Tom Tomato (Solanum lycopersicum L.) Plants Associated with a Gene Expression Analysis in Response to Short Daily UV Treatments" Plants 11, no. 14: 1829. https://doi.org/10.3390/plants11141829

APA Style

Mannucci, A., Santin, M., Vanhaelewyn, L., Sciampagna, M. C., Miras-Moreno, M. B., Zhang, L., Lucini, L., Quartacci, M. F., Van Der Straeten, D., Castagna, A., & Ranieri, A. (2022). Foliar and Root Comparative Metabolomics and Phenolic Profiling of Micro-Tom Tomato (Solanum lycopersicum L.) Plants Associated with a Gene Expression Analysis in Response to Short Daily UV Treatments. Plants, 11(14), 1829. https://doi.org/10.3390/plants11141829

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

Article Metrics

Back to TopTop