# Metabolic Profile Discriminates and Predicts Arabidopsis Susceptibility to Virus under Field Conditions

Bernadette Rubio, Olivier Fernandez, Patrick Cosson, Thierry Berton, Mélodie Caballero, Roxane Lion, Fabrice Roux, Joy Bergelson, Yves Gibon, Valerie Schurdi-Levraud

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

To cite this version:

Bernadette Rubio, Olivier Fernandez, Patrick Cosson, Thierry Berton, Mélodie Caballero, et al. Metabolic Profile Discriminates and Predicts Arabidopsis Susceptibility to Virus under Field Conditions. Metabolites, 2021, 11 (4), pp.230. 10.3390/metabo11040230. hal-03197197

HAL Id: hal-03197197

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

Submitted on 13 Apr 2021

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

# Metabolic Profile Discriminates and Predicts *Arabidopsis* Susceptibility to Virus under Field Conditions

Bernadette Rubio<fnref n="1" />, Olivier Fernandez <fnref n="1" />,†, Patrick Cosson <fnref n="1" />, Thierry Berton <fnref n="1" />, Mélodie Caballero <fnref n="1" />, Roxane Lion <fnref n="1" />, Fabrice Roux <fnref n="2" />, Joy Bergelson <fnref n="3" />, Yves Gibon <fnref n="1" /> and Valérie Schurdi-Levraud <fnref n="1" />,*

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

check for
updates

**Citation:** Rubio, B.; Fernandez, O.; Cosson, P.; Berton, T.; Caballero, M.; Lion, R.; Roux, F.; Bergelson, J.; Gibon, Y.; Schurdi-Levraud, V. *Metabolic Profile Discriminates and Predicts *Arabidopsis* Susceptibility to Virus under Field Conditions. *Metabolites* **2021**, *11*, 230. https://doi.org/10.3390/metabo11040230

Academic Editors: Victoria Pastor and Pierre Pétriacq

Received: 28 February 2021

Accepted: 2 April 2021

Published: 9 April 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** As obligatory parasites, plant viruses alter host cellular metabolism. There is a lack of information on the variability of virus-induced metabolic responses among genetically diverse plants in a natural context with daily changing conditions. To decipher the metabolic landscape of plant-virus interactions in a natural setting, twenty-six and ten accessions of *Arabidopsis thaliana* were inoculated with *Turnip mosaic virus* (TuMV), in two field experiments over 2 years. The accessions were measured for viral accumulation, above-ground biomass, targeted and untargeted metabolic profiles. The phenotypes of the accessions ranged from susceptibility to resistance. Susceptible and resistant accessions were shown to have different metabolic routes after inoculation. Susceptible genotypes accumulate primary and secondary metabolites upon infection, at the cost of hindered growth. Twenty-one metabolic signatures significantly accumulated in resistant accessions whereas they maintained their growth as mock-inoculated plants without biomass penalty. Metabolic content was demonstrated to discriminate and be highly predictive of the susceptibility of inoculated *Arabidopsis*. This study is the first to describe the metabolic landscape of plant-virus interactions in a natural setting and its predictive link to susceptibility. It provides new insights on plant-virus interactions. In this undomesticated species and in ecologically realistic conditions, growth and resistance are in a permanent conversation.

**Keywords:** *Arabidopsis thaliana*; field conditions; growth; central metabolism; specialized metabolism; trade-off; *Turnip mosaic virus*

## 1. Introduction

Plant health is of primary importance to improve and secure food supply for an increasing human population. Plant viruses represent the major taxonomic group of emergent pathogens of plants [1], and viral infection is one of the most alarming biotic threats due to the impact of climate change on the spatial and temporal distribution of vectors and viruses [2,3]. Compared to other plant pathogens, viruses are particularly unpredictable and difficult to combat. In this context, an understanding of the response of plants to viral infection has great importance in sustainable agricultural solutions.

The genus *Potyvirus*, to which the turnip mosaic virus (TuMV) belongs, is one of the largest genera among plant viruses, causing considerable economic damage in vegetable and fruit crops worldwide [4]. The completion of the viral multiplication and movement cycle results from a complex interplay between virus- and host-encoded factors that can have profound impacts on plant fitness. To invade plants, those obligatory parasites

<footnote n="1">¹ Université de Bordeaux, INRAE, Biologie du Fruit et Pathologie, UMR 1332, F-33140 Villenave d'Ornon, France; bernadette.rubio@inrae.fr (B.R.); olivier.fernandez@univ-reims.fr (O.F.); patrick.cosson@inrae.fr (P.C.); thierry.berton@laposte.net (T.B.); melodie.caballero@inrae.fr (M.C.); roxane.lion@inrae.fr (R.L.); yves.gibon@inrae.fr (Y.G.)</footnote>

<footnote n="2">² CNRS, INRAE, Université de Toulouse, LIPM, F-31320 Castanet-Tolosan, France; fabrice.roux@inrae.fr</footnote>

<footnote n="3">³ Ecology & Evolution, University of Chicago, 1101 E 57th St, Chicago, IL 60637, USA; jbergels@uchicago.edu</footnote>

<footnote>* Correspondence: valerie.schurdi-levraud@inrae.fr</footnote>

<footnote>† Current adress: SFR Condorcet FR CNRS 3417, Université de Reims Champagne-Ardenne, Résistance Induite et Bioprotection des Plantes EA 4707, BP 1039, CEDEX 2, F-51687 Reims, France.</footnote>
have developed tactics to reroute host cellular functions and components for their own benefits. Whatever the outcome of the interaction—compatibility, leading to disease, or incompatibility, leading to resistance—massive reprogramming of metabolism is observed. Indeed, as viruses perturb and exploit the host's carbon [5] and nitrogen [6] metabolism to make their own compounds, demand for energy within the plant increases to sustain viral multiplication, systemic spread and defense responses [7].

To date, viral impacts on plant metabolism have been mainly studied in experiments conducted in controlled conditions and on a limited diversity of host genotypes. In these conditions, respiration [8], photosynthetic efficiency [9] and carbon partitioning [6] were shown to be modified. For example, in a susceptible interaction between tobacco and potato virus Y (PVY), both an increase of soluble carbohydrates and a decrease of photosynthesis were observed 4 days after inoculation (dai) [10]. These observations were confirmed by the demonstration that potato leaves exhibit a decrease in sugar levels one day after PVY infection, with a subsequent increase in both inoculated and systemic leaves a few days later ([5,11]). Similarly, inoculation of susceptible hosts with cucumber mosaic virus (CMV) results in a localized reduction in starch accumulation as a consequence of altered carbohydrate metabolism at viral infection sites [12], although starch accumulates to high levels in systemically infected leaves [10]. Increases in amino acid content have also been demonstrated in terms of total or individual amino acids, as well as in polyamines, in a variety of systems (reviewed in [6,13,14]). Viruses interfere with fatty acids, structural components of intracellular membranes in which replication can take place [15]. They also exploit transport systems in order to invade cells in systemic tissues away from the initial site of infection. This is achieved through an interaction between viral proteins and components of the long distance transport machinery like phloem proteins [16]. Like primary metabolism, specialized metabolism is strongly impacted due to its participation in multiple defense signaling cascades [17]. Numerous compounds, including hormones, are involved in defense ([18,19]). Among them, glucosinolates and phytoalexins play significant roles in defense against a range of pathogens [20].

In a compatible interaction, the outcome of viral colonization can include symptoms such as stunting, chlorosis or necrosis depending on the pathosystem. Incompatible interactions trigger plant resistance and defense signaling that involve the action of antimicrobial components and specific defense proteins [17]. These defense responses are now widely acknowledged to involve a trade-off in model plants, with the cost of resistance generating a negative impact on plant fitness [21–27]. Similarly, in crop plants, high levels of resistance are often associated with yield penalties [28,29].

While extremely informative, these studies performed under controlled optimal conditions reduce environmental effects and increase the likelihood of finding strong relationships between metabolite levels variations [30]. Thereby, it obscures the complexity and the variability of metabolic responses among genetically diverse plants in a natural context where they have to face multiple stresses. Even in main crop species, large-scale metabolic profiling of field-grown populations of genetically diverse accessions remains rare [31–33].

In *Arabidopsis thaliana*, metabolite profiling of large populations of natural accessions or inbred lines has allowed the identification of descriptor sets of metabolites that are predictive of biomass [32,34,35] and physiological traits such as freezing tolerance [36] and herbivore resistance [37]. But, none of these studies addressed the growth and metabolic response to biotic stress in a changing environment.

Here, following the modern standards of ecological genomics [38], we aimed at improving our understanding of the virus-induced reprogramming of metabolism through deciphering variation in the metabolic landscape of the natural *Arabidopsis thaliana/turnip mosaic virus (TuMV) pathosystem* in growth conditions close to environmental reality. To do so, we set up two experiments in the field and explored targeted and untargeted metabolic profiles on 26 accessions and 10 accessions, in 2015 and 2017 respectively. These accessions spanned responses ranging from high susceptibility to full resistance. This work aimed to (i) decipher the trade-offs among responses to TuMV, metabolism and growth
in field conditions in contrasted *Arabidopsis* accessions, (ii) characterize the metabolic disturbances and its kinetics for a contrasted set of accessions, and (iii) identify discriminant metabolic biomarkers of susceptibility/resistance to viruses in *A. thaliana*.

## 2. Results

### 2.1. In the Field, Susceptible Accessions Accumulate PRIMARY METABolites upon Infection, at the Cost of Hindered Growth, Whereas Resistant Accessions Grow with Limited Changes

To analyze the relationship between infection, growth and central metabolism, a set of 26 accessions [39] were mock and TuMV-inoculated in a field experiment in 2015. Thirteen days after inoculation (dai), the viral load of the samples, their above-ground dry biomass and their composition for 10 metabolites features were evaluated.

The quantification of the viral load by ELISA confirmed that eight genotypes out of twenty-six TuMV-inoculated were considered fully resistant as no virus was detected (Table 1). The remaining eighteen accessions fell within a range of susceptibility, OD from 0.258 to 0.66 (Table 1).

<table><caption>**Table 1.** List of the *A. thaliana* accessions with their geographic position and their susceptible status to TuMV infection at 13 dai confirmed by OD values and its SD in the 2015 field experiment.</caption><thead><tr><th>Genotype</th><th>ID <sup>1</sup></th><th>Latitude</th><th>Longitude</th><th>Country</th><th>OD Means</th><th>OD-Standard Deviations</th><th>Susceptibility Groups <sup>2</sup></th></tr></thead><tbody><tr><td>Bay-0</td><td>6899</td><td>49</td><td>11</td><td>GER</td><td>0.004</td><td>0.0155</td><td>R</td></tr><tr><td>Mar2-3</td><td>159</td><td>47.35</td><td>3.93333</td><td>FRA</td><td>0.005</td><td>0.0087</td><td>R</td></tr><tr><td>MIB-39</td><td>190</td><td>47.3833</td><td>5.31667</td><td>FRA</td><td>0.007</td><td>0.0155</td><td>R</td></tr><tr><td>Petergof</td><td>7296</td><td>59</td><td>29</td><td>RUS</td><td>0.011</td><td>0.0252</td><td>R</td></tr><tr><td>Ra-0</td><td>6958</td><td>46</td><td>3.3</td><td>FRA</td><td>0.003</td><td>0.0072</td><td>R</td></tr><tr><td>Rak-2</td><td>8365</td><td>49</td><td>16</td><td>CZE</td><td>0.006</td><td>0.00204</td><td>R</td></tr><tr><td>TOU-A1-69</td><td>335</td><td>46.6667</td><td>4.11667</td><td>FRA</td><td>0.015</td><td>0.0224</td><td>R</td></tr><tr><td>TOU-A1-84</td><td>348</td><td>46.6667</td><td>4.11667</td><td>FRA</td><td>0.007</td><td>0.0204</td><td>R</td></tr><tr><td>Edi-0</td><td>6914</td><td>56</td><td>-3</td><td>UK</td><td>0.258</td><td>0.169</td><td>S</td></tr><tr><td>Bg-2</td><td>6709</td><td>47.6479</td><td>-122.305</td><td>USA</td><td>0.324</td><td>0.211</td><td>S</td></tr><tr><td>MIB-62</td><td>206</td><td>47.3833</td><td>5.31667</td><td>FRA</td><td>0.346</td><td>0.2084</td><td>S</td></tr><tr><td>Gu-0</td><td>6922</td><td>50.3</td><td>8</td><td>GER</td><td>0.389</td><td>0.272</td><td>S</td></tr><tr><td>Col-0</td><td>6909</td><td>38.3</td><td>-92.3</td><td>USA</td><td>0.391</td><td>0.215</td><td>S</td></tr><tr><td>Bu-0</td><td>8271</td><td>50.5</td><td>9.5</td><td>GER</td><td>0.4</td><td>0.203</td><td>S</td></tr><tr><td>TOU-A1-125</td><td>291</td><td>46.6667</td><td>4.11667</td><td>FRA</td><td>0.403</td><td>0.147</td><td>S</td></tr><tr><td>TOU-L-5</td><td>389</td><td>46.6667</td><td>4.11667</td><td>FRA</td><td>0.407</td><td>0.19</td><td>S</td></tr><tr><td>MIB-60</td><td>204</td><td>47.3833</td><td>5.31667</td><td>FRA</td><td>0.425</td><td>0.145</td><td>S</td></tr><tr><td>MIB-28</td><td>178</td><td>47.3833</td><td>5.31667</td><td>FRA</td><td>0.428</td><td>0.242</td><td>S</td></tr><tr><td>Hs-0</td><td>8310</td><td>52.24</td><td>9.44</td><td>GER</td><td>0.46</td><td>0.212</td><td>S</td></tr><tr><td>Mt-0</td><td>6939</td><td>32.34</td><td>22.46</td><td>LIB</td><td>0.467</td><td>0.2</td><td>S</td></tr><tr><td>MIB-20</td><td>171</td><td>47.3833</td><td>5.31667</td><td>FRA</td><td>0.525</td><td>0.294</td><td>S</td></tr><tr><td>MIB-67</td><td>210</td><td>47.3833</td><td>5.31667</td><td>FRA</td><td>0.549</td><td>0.249</td><td>S</td></tr><tr><td>TOU-A1-73</td><td>338</td><td>46.6667</td><td>4.11667</td><td>FRA</td><td>0.55</td><td>0.261</td><td>S</td></tr><tr><td>CUR-10</td><td>79</td><td>45</td><td>1.75</td><td>FRA</td><td>0.576</td><td>0.147</td><td>S</td></tr><tr><td>Mz-0</td><td>6940</td><td>50.3</td><td>8.3</td><td>GER</td><td>0.59</td><td>0.263</td><td>S</td></tr><tr><td>Bs-1</td><td>8270</td><td>47.5</td><td>7.5</td><td>SUI</td><td>0.66</td><td>0.339</td><td>S</td></tr></tbody></table>

Aboveground dry biomass was measured on mock-inoculated and TuMV-inoculated plants. Aboveground dry biomass was significantly lower for the 18 inoculated susceptible accessions compared to those susceptible accessions when treated only with a mock solution, indicating that infection led to a decrease in the growth of susceptible plants. There was no significant difference between the aboveground dry biomass of the eight resistant inoculated accessions compared to their mock-inoculated counterparts (Kruskal-Wallis test; Figure 1).

<footnote n="1">¹ ID according to [40]. ² In 2015, categories have been defined according to the healthy control Col-0 which mean OD value was 0.088 (SD 0.01063). Infected genotypes with mean OD ≤ 0.088 were defined as resistant (R). Infected genotypes with mean OD > 2.5*0.088 were defined as susceptible [41].</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 1.** Above ground Dry biomass for mock-inoculated and TuMV-inoculated susceptible and resistant accessions measured at 13 days after inoculation (dai) on 26 accessions in the '2015' field experiment. Statistical comparisons (Kruskal-Wallis test) on dry biomass were performed between all categories of accessions (N = 21 for Mock and inoculated susceptible categories; N = 8 for mock and inoculated resistant categories). Plots with the same letter are not significantly different at *p* = 0.05.</figcaption></figure>

Ten key metabolic features, amino acids, proteins, glutamate, malate, fumarate, starch, glucose, fructose, sucrose and chlorophyll a, were measured in 2015 for each accession. Thirteen days after inoculation, analysis of the Spearman correlation pattern between dry biomass, OD and the content in the different metabolites features (Table S1 in Supplementary Materials) showed that that viral accumulation was negatively correlated with dry aboveground biomass. In mock-inoculated samples, a significant positive correlation was found between malate, fumarate, starch, glucose, fructose, sucrose and aboveground dry biomass whereas a negative significant correlation was found with amino acids, proteins, glutamate (Table S1). Compared to mock-inoculated plants, the remodeling of metabolic contents in TuMV-inoculated plants was shown as biomass appeared negatively correlated with nine out of 10 of the measured metabolites features (all except chlorophyll a). Whereas, on the contrary, the viral load appears positively correlated to these same nine metabolites' (amino acids, proteins, glutamate, malate, fumarate, starch, glucose, fructose and sucrose) features (Table S1).

A detailed analysis showed that the 18 TuMV-inoculated susceptible accessions exhibited a significant accumulation of eight out of 10 central metabolites (i.e., amino acids, glutamate, malate, fumarate, starch, glucose, fructose and sucrose), when compared to either the eight resistant inoculated accessions or the set of mock inoculated controls (Figure 2). The chlorophyll content in TuMV-inoculated susceptible and resistant accessions was not significantly reduced, confirming that sampling occurred prior to macroscopic symptoms of chlorosis.

TuMV inoculated susceptible accessions are represented by pink boxplot (light pink for Mock-inoculated and dark pink for TuMV inoculated) and resistant accessions by green boxplots (light green for Mock-inoculated and dark green for TuMV inoculated). Statistical analyzes were performed on each of the 10 metabolites features (according to the normality of the data, ANOVA or Kruskall-Wallis test was performed). For each metabolite, plots with the same letter are not significantly different at *p* = 0.05.
<figure><img src="image_7.png" /><figcaption>**Figure 2.** Comparisons of central metabolic content in mock-inoculated and TuMV-inoculated resistant and susceptible *A. thaliana* accessions in the 2015 field experiment. Boxplot sharing the same letter are not significantly different, at the chosen level (default 5%).</figcaption></figure>

The TuMV-inoculated susceptible genotypes therefore accumulate the majority of the measured metabolites features without being used for their growth, whereas, in addition to maintaining their growth, the inoculated resistant accessions exhibited similar amounts of metabolites features relative to those observed for mock-inoculated controls.

## 2.2. *Metabolic Content Discriminates Inoculated A. thaliana Susceptible and Resistant Accessions*

Principal component analysis (PCA) performed with the 10 primary metabolic traits measured at 13 days after inoculation on mock-inoculated and TuMV-inoculated accessions showed that targeted metabolism does not distinguish between susceptible and resistant genotypes when mock- inoculated (Figure 3a).

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Principal component analysis (PCA) performed on 26 *A. thaliana* accessions with 10 primary metabolic traits measured at 13 dai on mock-inoculated plants and TuMV-inoculated plants in the '2015' field experiment. (a) Score plot of individuals. Score plots of resistant accessions are in green, light green for mock-inoculated, dark green for TuMV-inoculated. Score plots of susceptible accessions are in pink, light pink for mock-inoculated, dark pink for TuMV-inoculated. The confidence ellipses around the centroid of individuals are represented. (b) Variable correlation graph. Dry biomass was considered as explanatory traits (in blue on the variable factor map). Chla = Chlorophyll a.</figcaption></figure>
While metabolic change lead to a clear differentiation between susceptible and resistant accessions that have been inoculated with TuMV (Figure 3a), the first axe of the PCA discriminating between susceptible and resistant accessions (46.67% of the explained variance).

The contribution of primary metabolites also differed between mock-inoculated plants and TuMV-inoculated plants with clear contribution of most of the primary metabolites to response in susceptible genotypes (Figure 3b).

## 2.3. *Metabolic Differentiation between Mock-Inoculated and TuMV-Inoculated Genotypes Appears as Early as 5 Days after Inoculation*

In order to specify the kinetics of setting up the metabolic differentiation between susceptible and resistant genotypes, TuMV-inoculated and mock, a field experiment was conducted in 2017. Five resistant genotypes and five susceptible genotypes from the 2015 set were mock-inoculated or inoculated with TuMV. The plants were collected on the day of inoculation (T0) and at 5, 7, 9 and 13 days after inoculation. The viral load was measured at each sampling time. The optical density data taken at 13 dai (Table S2) confirm the classification between susceptible and resistant observed in 2015 for these genotypes.

The PCA carried out using the measurement data of the 10 central metabolites show that at T0, in a non-inoculated situation, the metabolic composition does not allow to distinguish between susceptible and resistant genotypes (Figure 4 T0a). Sugars, fumarate and malate appear positively correlated and clearly participate in the biomass (Figure 4 T0b).

Five days after inoculation, a clear differentiation appears between TuMV-inoculated genotypes and mock-inoculated genotypes, independently of the resistance or susceptibility status (Figure 4 T5a). At 7 and 9 days after inoculation, the pattern appears less clear but the mock-inoculated genotypes, resistant or susceptible, still remain grouped (Figure 4 T7 and T9a) whereas susceptible and resistant genotypes TuMV-inoculated have a different metabolic time course. Thirteen days after inoculation, and as already observed in 2015, the clear differentiation, in terms of metabolic composition, between inoculated resistant genotypes and inoculated susceptible genotypes is in place (Figure 4 T13a).

## 2.4. *The Identification of Metabolic Predictors Reveals 21 Metabolic Signatures Significantly Accumulated in Resistant Accessions*

To analyze the pattern of infection and its relationship to a larger set of metabolites, including specialized metabolites, we conducted an untargeted metabolic analysis by UHPLC-LTQ Orbitrap on the set of 26 accessions showing a contrasting response to TuMV in the 2015 field experiment. This analysis captured a total of 505 workable metabolic signatures (m/z), corresponding primarily to non-polar specialized metabolites.

As previously observed with targeted central metabolism, there was no difference between susceptible and resistant accessions when mock-inoculated (Figure S1a). In contrast, when inoculated with TuMV, principal component analysis on the metabolic variables clearly separated two groups of accessions according to their susceptibility to the TuMV (Figure S1b).

An orthogonal partial least squares-discriminant analysis (OPLS-DA) was then performed to maximize the variation between the two groups of accessions, mock and TuMV-inoculated, and determine the most significant variables contributing to this variation, i.e., variable importance in the projection (VIPs). OPLS-DA analysis was carried out using the central and specialized metabolite data. The quality of the model was validated by the Q² parameter (goodness-of-prediction parameter) with a value of 0.846, thereby showing high predictive capabilities (Figure S2).
<figure><img src="image_9.png" /><figcaption>**Figure 4.** Principal component analysis (PCA) performed with 10 primary metabolic traits measured at 0 (T0), 5 (T5), 7 (T7), 9 (T9) and 13 (T13) days after inoculation on 10 *A. thaliana* accessions mock and TuMV inoculated in the '2017' field experiment. **a** Score plot of individuals. Score plots of resistant accessions are in green, light green for mock-inoculated, dark green for TuMV-inoculated. Score plots of susceptible accessions are in pink, light pink for mock-inoculated, dark pink for TuMV-inoculated. The confidence ellipses around the centroid of individuals are represented. **b** Variable correlation graph. Dry biomass was considered as explanatory traits (in blue on the variable factor map). Chla = Chlorophyll a.</figcaption></figure>

The OPLS-DA analysis performed between the mock- and TuMV-inoculated samples revealed 63 common discriminant metabolic variables (Table S3), most of which (58/63) accumulated in TuMV-inoculated samples and particularly in susceptible accessions (54/58; Table S3). It is worth noting that four VIPs (433;700, 361;491, 512;491 and 64;395) are significantly found in resistant accessions.

To examine which metabolic variables strongly contribute to the OPLS-DA model in TuMV-inoculated susceptible and resistant accessions, variables were ranked according to their VIP values (Table 2).
<table><caption>**Table 2.** List of the variable importance in the projection (VIPs) identify by OPLS-DA analysis performed on TuMV-inoculated resistant and susceptible twenty-six accessions in 2015. For each VIP, comparisons between resistant (R) and susceptible (S) accessions were done. The fold change was calculated for each VIP. Primary metabolites are light-grey highlighted. Metabolites that accumulate significantly more in resistant accessions are at the bottom of the table.</caption><thead><tr><th>VIP OPLS-DA <sup>1</sup></th><th>VIP Values</th><th>m/z <sup>2</sup></th><th>rt <sup>3</sup></th><th>Resistant vs. Susceptible Metabolic Contents</th><th>Fold Change</th></tr></thead><tbody><tr><td>Sucrose</td><td>2.213</td><td>NA</td><td>NA</td><td>S > R *** <sup>4</sup></td><td>2.89</td></tr><tr><td>303;463</td><td>1.973</td><td>303.133</td><td>463.27</td><td>S > R ***</td><td>3.69</td></tr><tr><td>356;452</td><td>1.958</td><td>356.12</td><td>452.352</td><td>S > R ***</td><td>2.17</td></tr><tr><td>219;311</td><td>1.94</td><td>219.101</td><td>311.134</td><td>S > R ***</td><td>4.81</td></tr><tr><td>533;392</td><td>1.926</td><td>533.155</td><td>391.765</td><td>S > R ***</td><td>4.98</td></tr><tr><td>332;430</td><td>1.898</td><td>332.132</td><td>430.481</td><td>S > R ***</td><td>3.78</td></tr><tr><td>116;100</td><td>1.867</td><td>116.07</td><td>100.163</td><td>S > R ***</td><td>6.05</td></tr><tr><td>205;242</td><td>1.837</td><td>205.097</td><td>242.121</td><td>S > R ***</td><td>2.51</td></tr><tr><td>103;242</td><td>1.799</td><td>103.041</td><td>242.089</td><td>S > R ***</td><td>2.86</td></tr><tr><td>302;506</td><td>1.796</td><td>302.102</td><td>506.314</td><td>S > R ***</td><td>5.76</td></tr><tr><td>255;544</td><td>1.769</td><td>255.112</td><td>543.799</td><td>S > R ***</td><td>9.96</td></tr><tr><td>221;216</td><td>1.748</td><td>221.092</td><td>215.744</td><td>S > R ***</td><td>3.13</td></tr><tr><td>385;210</td><td>1.738</td><td>385.106</td><td>210.083</td><td>S > R ***</td><td>45.77</td></tr><tr><td>Glucose</td><td>1.734</td><td>NA</td><td>NA</td><td>S > R ***</td><td>2.49</td></tr><tr><td>175;402</td><td>1.713</td><td>175.148</td><td>401.831</td><td>S > R ***</td><td>3.77</td></tr><tr><td>903;319</td><td>1.71</td><td>903.277</td><td>318.889</td><td>S > R ***</td><td>8.48</td></tr><tr><td>503;391</td><td>1.695</td><td>503.19</td><td>391.457</td><td>S > R ***</td><td>2.51</td></tr><tr><td>343;345</td><td>1.693</td><td>343.117</td><td>345.319</td><td>S > R ***</td><td>2.08</td></tr><tr><td>Fructose</td><td>1.683</td><td>NA</td><td>NA</td><td>S > R ***</td><td>2.49</td></tr><tr><td>474;107</td><td>1.675</td><td>474.218</td><td>107.298</td><td>S > R ***</td><td>2.68</td></tr><tr><td>315;448</td><td>1.674</td><td>315.133</td><td>447.821</td><td>S > R ***</td><td>2.39</td></tr><tr><td>209;488</td><td>1.669</td><td>209.153</td><td>488.35</td><td>S > R ***</td><td>2.7</td></tr><tr><td>209;414</td><td>1.663</td><td>209.153</td><td>413.756</td><td>S > R ***</td><td>2.14</td></tr><tr><td>212;284</td><td>1.657</td><td>211.559</td><td>283.891</td><td>S > R ***</td><td>1.56</td></tr><tr><td>543;99</td><td>1.645</td><td>543.132</td><td>98.5221</td><td>S > R ***</td><td>5.84</td></tr><tr><td>124;346</td><td>1.638</td><td>124.075</td><td>346.354</td><td>S > R ***</td><td>1.92</td></tr><tr><td>370;302</td><td>1.61</td><td>370.149</td><td>301.519</td><td>S > R ***</td><td>28.49</td></tr><tr><td>203;204</td><td>1.591</td><td>203.084</td><td>203.918</td><td>S > R ***</td><td>8.24</td></tr><tr><td>Glutamate</td><td>1.59</td><td>NA</td><td>NA</td><td>S > R ***</td><td>1.47</td></tr><tr><td>430;333</td><td>1.581</td><td>430.17</td><td>332.878</td><td>S > R ***</td><td>3.3</td></tr><tr><td>331;329</td><td>1.578</td><td>331.117</td><td>328.817</td><td>S > R ***</td><td>1.58</td></tr><tr><td>226;406</td><td>1.554</td><td>226.107</td><td>406.369</td><td>S > R ***</td><td>7.97</td></tr><tr><td>151;359</td><td>1.549</td><td>151.075</td><td>359.377</td><td>S > R ***</td><td>1.78</td></tr><tr><td>449;319</td><td>1.545</td><td>449.106</td><td>318.949</td><td>S > R ***</td><td>1.75</td></tr><tr><td>757;319</td><td>1.537</td><td>757.217</td><td>319.13</td><td>S > R ***</td><td>1.69</td></tr><tr><td>270;587</td><td>1.534</td><td>270.133</td><td>587.091</td><td>S > R ***</td><td>43.68</td></tr><tr><td>315;370</td><td>1.524</td><td>315.133</td><td>370.476</td><td>S > R ***</td><td>3.93</td></tr></tbody></table>
<table><caption>Table 2. *Cont.*</caption><thead><tr><th>VIP OPLS-DA <sup>1</sup></th><th>VIP Values</th><th>m/z <sup>2</sup></th><th>rt <sup>3</sup></th><th>Resistant vs. Susceptible Metabolic Contents</th><th>Fold Change</th></tr></thead><tbody><tr><td>162;402</td><td>1.511</td><td>162.055</td><td>402.281</td><td>S > R ***</td><td>1.9</td></tr><tr><td>221;230</td><td>1.503</td><td>221.121</td><td>230.244</td><td>S > R ***</td><td>1.74</td></tr><tr><td>394;517</td><td>1.49</td><td>394.204</td><td>517.072</td><td>S > R ***</td><td>17.28</td></tr><tr><td>302;407</td><td>1.486</td><td>302.101</td><td>406.697</td><td>S > R ***</td><td>7.16</td></tr><tr><td>482;101</td><td>1.475</td><td>482.107</td><td>101.274</td><td>S > R ***</td><td>1.79</td></tr><tr><td>Amino Acids</td><td>1.462</td><td>NA</td><td>NA</td><td>S > R ***</td><td>1.88</td></tr><tr><td>355;306</td><td>1.451</td><td>355.102</td><td>306.288</td><td>S > R ***</td><td>3.29</td></tr><tr><td>191;416</td><td>1.448</td><td>191.143</td><td>416.467</td><td>S > R ***</td><td>2.49</td></tr><tr><td>193;324</td><td>1.445</td><td>193.125</td><td>323.817</td><td>S > R ***</td><td>1.46</td></tr><tr><td>149;360</td><td>1.437</td><td>149.096</td><td>360.18</td><td>S > R ***</td><td>1.66</td></tr><tr><td>642;445</td><td>1.435</td><td>642.254</td><td>445.078</td><td>S > R ***</td><td>2</td></tr><tr><td>321;760</td><td>1.419</td><td>321.114</td><td>760.442</td><td>S > R ***</td><td>3.7</td></tr><tr><td>182;466</td><td>1.415</td><td>182.081</td><td>465.796</td><td>S > R ***</td><td>4.21</td></tr><tr><td>305;210</td><td>1.408</td><td>305.086</td><td>210.166</td><td>S > R ***</td><td>2.24</td></tr><tr><td>348;319</td><td>1.396</td><td>348.274</td><td>318.945</td><td>S > R ***</td><td>3.61</td></tr><tr><td>164;496</td><td>1.375</td><td>164.07</td><td>496.307</td><td>S > R ***</td><td>2.22</td></tr><tr><td>165;424</td><td>1.372</td><td>165.127</td><td>424.14</td><td>S > R ***</td><td>1.54</td></tr><tr><td>317;464</td><td>1.37</td><td>317.101</td><td>463.65</td><td>S > R ***</td><td>1.25</td></tr><tr><td>191;363</td><td>1.36</td><td>191.07</td><td>362.643</td><td>S > R ***</td><td>2.76</td></tr><tr><td>219;486</td><td>1.34</td><td>219.101</td><td>486.477</td><td>S > R ***</td><td>2.7</td></tr><tr><td>105;700</td><td>1.33</td><td>105.069</td><td>700.287</td><td>S > R ***</td><td>1.59</td></tr><tr><td>367;358</td><td>1.323</td><td>367.101</td><td>358.491</td><td>S > R ***</td><td>1.69</td></tr><tr><td>201;344</td><td>1.319</td><td>201.054</td><td>343.832</td><td>S > R ***</td><td>1.71</td></tr><tr><td>189;325</td><td>1.312</td><td>189.127</td><td>324.754</td><td>S > R ***</td><td>1.74</td></tr><tr><td>146;242</td><td>1.311</td><td>146.06</td><td>241.843</td><td>S > R ***</td><td>2.31</td></tr><tr><td>133;296</td><td>1.307</td><td>133.064</td><td>295.68</td><td>S > R ***</td><td>2.65</td></tr><tr><td>179;604</td><td>1.306</td><td>179.106</td><td>603.961</td><td>S > R ***</td><td>1.64</td></tr><tr><td>374;322</td><td>1.305</td><td>374.144</td><td>321.703</td><td>S > R ***</td><td>1.77</td></tr><tr><td>373;453</td><td>1.29</td><td>373.127</td><td>452.6</td><td>S > R ***</td><td>1.81</td></tr><tr><td>302;382</td><td>1.282</td><td>302.041</td><td>381.707</td><td>S > R ***</td><td>1.54</td></tr><tr><td>161;442</td><td>1.274</td><td>161.096</td><td>441.797</td><td>S > R ***</td><td>1.79</td></tr><tr><td>164;375</td><td>1.269</td><td>164.07</td><td>375.012</td><td>S > R ***</td><td>2.95</td></tr><tr><td>109;359</td><td>1.267</td><td>109.064</td><td>359.34</td><td>S > R ***</td><td>1.92</td></tr><tr><td>386;278</td><td>1.256</td><td>386.22</td><td>277.549</td><td>S > R ***</td><td>3.48</td></tr><tr><td>391;774</td><td>1.255</td><td>391.245</td><td>774.477</td><td>S > R ***</td><td>2.29</td></tr><tr><td>409;491</td><td>1.251</td><td>409.169</td><td>490.545</td><td>S > R ***</td><td>27.12</td></tr><tr><td>291;259</td><td>1.241</td><td>291.181</td><td>259.063</td><td>S > R ***</td><td>4.78</td></tr><tr><td>192;462</td><td>1.227</td><td>192.041</td><td>462.37</td><td>S > R ***</td><td>1.27</td></tr><tr><td>107;370</td><td>1.219</td><td>107.085</td><td>370.001</td><td>S > R ***</td><td>1.39</td></tr><tr><td>420;276</td><td>1.217</td><td>419.695</td><td>275.797</td><td>S > R ***</td><td>1.77</td></tr><tr><td>193;373</td><td>1.189</td><td>193.086</td><td>373.266</td><td>S > R ***</td><td>2.04</td></tr></tbody></table>
<table><caption>Table 2. *Cont.*</caption><thead><tr><th>VIP OPLS-DA <sup>1</sup></th><th>VIP Values</th><th>m/z <sup>2</sup></th><th>rt <sup>3</sup></th><th>Resistant vs. Susceptible Metabolic Contents</th><th>Fold Change</th></tr></thead><tbody><tr><td>379;402</td><td>1.179</td><td>379.095</td><td>402.311</td><td>S > R ***</td><td>2.88</td></tr><tr><td>181;464</td><td>1.176</td><td>181.086</td><td>463.916</td><td>S > R ***</td><td>1.42</td></tr><tr><td>181;328</td><td>1.17</td><td>181.086</td><td>327.894</td><td>S > R ***</td><td>1.68</td></tr><tr><td>80;395</td><td>1.168</td><td>80.049</td><td>395.289</td><td>S > R ***</td><td>1.57</td></tr><tr><td>373;344</td><td>1.16</td><td>373.127</td><td>343.531</td><td>S > R ***</td><td>1.54</td></tr><tr><td>627;299</td><td>1.156</td><td>627.155</td><td>299.369</td><td>S > R ***</td><td>2.8</td></tr><tr><td>96;389</td><td>1.155</td><td>96.08</td><td>389.004</td><td>S > R ***</td><td>1.74</td></tr><tr><td>195;453</td><td>1.154</td><td>195.065</td><td>452.614</td><td>S > R ***</td><td>1.55</td></tr><tr><td>79;396</td><td>1.152</td><td>79.041</td><td>395.754</td><td>S > R ***</td><td>1.52</td></tr><tr><td>105;327</td><td>1.15</td><td>105.069</td><td>326.898</td><td>S > R ***</td><td>1.72</td></tr><tr><td>210;465</td><td>1.147</td><td>210.112</td><td>464.799</td><td>S > R ***</td><td>2.44</td></tr><tr><td>335;506</td><td>1.143</td><td>335.127</td><td>506.441</td><td>S > R ***</td><td>3.71</td></tr><tr><td>86;126</td><td>1.114</td><td>86.096</td><td>126.246</td><td>S > R ***</td><td>2.09</td></tr><tr><td>396;257</td><td>1.106</td><td>396.185</td><td>256.51</td><td>S > R ***</td><td>1.43</td></tr><tr><td>611;354</td><td>1.102</td><td>611.158</td><td>353.603</td><td>S > R ***</td><td>1.34</td></tr><tr><td>178;191</td><td>1.102</td><td>178.089</td><td>191.013</td><td>S > R ***</td><td>3.43</td></tr><tr><td>169;496</td><td>1.101</td><td>169.049</td><td>495.5</td><td>S > R ***</td><td>3.67</td></tr><tr><td>162;249</td><td>1.098</td><td>162.055</td><td>249.101</td><td>S > R ***</td><td>1.5</td></tr><tr><td>396;348</td><td>1.092</td><td>396.115</td><td>347.739</td><td>S > R ***</td><td>4.14</td></tr><tr><td>103;389</td><td>1.086</td><td>103.054</td><td>389.444</td><td>S > R ***</td><td>2</td></tr><tr><td>521;325</td><td>1.084</td><td>521.201</td><td>325.328</td><td>S > R ***</td><td>2.52</td></tr><tr><td>133;126</td><td>1.072</td><td>133.105</td><td>126.347</td><td>S > R **</td><td>2.32</td></tr><tr><td>209;426</td><td>1.071</td><td>209.153</td><td>425.939</td><td>S > R ***</td><td>1.98</td></tr><tr><td>162;357</td><td>1.069</td><td>162.055</td><td>356.879</td><td>S > R ***</td><td>1.74</td></tr><tr><td>201;451</td><td>1.063</td><td>201.054</td><td>451.123</td><td>S > R ***</td><td>1.75</td></tr><tr><td>527;460</td><td>1.056</td><td>527.103</td><td>460.095</td><td>S > R ***</td><td>2.02</td></tr><tr><td>Fumarate</td><td>1.054</td><td>NA</td><td>NA</td><td>S > R ***</td><td>1.38</td></tr><tr><td>393;416</td><td>1.052</td><td>393.188</td><td>415.816</td><td>S > R ***</td><td>1.71</td></tr><tr><td>212;774</td><td>1.052</td><td>212.094</td><td>774.047</td><td>S > R ***</td><td>2.76</td></tr><tr><td>402;420</td><td>1.048</td><td>402.162</td><td>419.597</td><td>S > R **</td><td>1.87</td></tr><tr><td>404;388</td><td>1.038</td><td>404.227</td><td>388.155</td><td>S > R **</td><td>1.31</td></tr><tr><td>464;367</td><td>1.035</td><td>464.248</td><td>366.989</td><td>S > R ***</td><td>1.97</td></tr><tr><td>225;371</td><td>1.034</td><td>225.148</td><td>370.798</td><td>S > R ***</td><td>1.57</td></tr><tr><td>103;284</td><td>1.033</td><td>103.054</td><td>283.823</td><td>S > R ***</td><td>3.26</td></tr><tr><td>227;789</td><td>1.023</td><td>227.163</td><td>788.506</td><td>S > R ***</td><td>1.61</td></tr><tr><td>222;216</td><td>1.016</td><td>221.602</td><td>215.822</td><td>S > R **</td><td>1.33</td></tr><tr><td>367;342</td><td>1.011</td><td>367.153</td><td>341.908</td><td>S > R ***</td><td>1.96</td></tr><tr><td>302;463</td><td>1.006</td><td>302.102</td><td>463.078</td><td>S > R ***</td><td>1.49</td></tr><tr><td>309;325</td><td>1.003</td><td>309.116</td><td>324.659</td><td>S > R **</td><td>1.44</td></tr><tr><td>244;325</td><td>1.001</td><td>244.096</td><td>324.677</td><td>S > R **</td><td>1.41</td></tr></tbody></table>
<table><caption>**Table 2. Cont.**</caption><thead><tr><th>VIP OPLS-DA <sup>1</sup></th><th>VIP Values</th><th>m/z <sup>2</sup></th><th>rt <sup>3</sup></th><th>Resistant vs. Susceptible Metabolic Contents</th><th>Fold Change</th></tr></thead><tbody><tr><td>351;184</td><td>2.047</td><td>351.006</td><td>184.222</td><td>R > S ***</td><td>2.88</td></tr><tr><td>432;184</td><td>1.831</td><td>431.971</td><td>183.776</td><td>R > S ***</td><td>3.44</td></tr><tr><td>137;131</td><td>1.618</td><td>136.931</td><td>131.31</td><td>R > S ***</td><td>1.39</td></tr><tr><td>512;491</td><td>1.563</td><td>512.127</td><td>491.367</td><td>R > S ***</td><td>2.54</td></tr><tr><td>394;108</td><td>1.548</td><td>394.2</td><td>108.419</td><td>R > S ***</td><td>1.93</td></tr><tr><td>337;573</td><td>1.51</td><td>337.292</td><td>572.917</td><td>R > S ***</td><td>1.64</td></tr><tr><td>299;181</td><td>1.396</td><td>299.098</td><td>181.363</td><td>R > S ***</td><td>22.3</td></tr><tr><td>181;139</td><td>1.293</td><td>181.053</td><td>138.934</td><td>R > S ***</td><td>1.94</td></tr><tr><td>324;184</td><td>1.263</td><td>323.989</td><td>183.797</td><td>R > S ***</td><td>106.8</td></tr><tr><td>280;184</td><td>1.229</td><td>280.084</td><td>183.768</td><td>R > S ***</td><td>5.9</td></tr><tr><td>361;491</td><td>1.218</td><td>361.092</td><td>490.564</td><td>R > S ***</td><td>2.07</td></tr><tr><td>155;132</td><td>1.141</td><td>154.941</td><td>131.541</td><td>R > S ***</td><td>1.33</td></tr><tr><td>433;482</td><td>1.114</td><td>433.112</td><td>481.588</td><td>R > S ***</td><td>244</td></tr><tr><td>433;700</td><td>1.109</td><td>433.241</td><td>700.261</td><td>R > S ***</td><td>1.82</td></tr><tr><td>64;395</td><td>1.091</td><td>63.934</td><td>395.047</td><td>R > S ***</td><td>3.6</td></tr><tr><td>256;668</td><td>1.089</td><td>256.08</td><td>667.874</td><td>R > S ***</td><td>2.12</td></tr><tr><td>460;740</td><td>1.076</td><td>460.269</td><td>740.005</td><td>R > S **</td><td>1.42</td></tr><tr><td>449;442</td><td>1.069</td><td>449.107</td><td>442.311</td><td>R > S ***</td><td>6.17</td></tr><tr><td>244;130</td><td>1.049</td><td>243.942</td><td>130.167</td><td>R > S ***</td><td>1.35</td></tr><tr><td>347;240</td><td>1.039</td><td>347.159</td><td>239.659</td><td>R > S **</td><td>1.48</td></tr><tr><td>86;184</td><td>1.015</td><td>86.059</td><td>184.278</td><td>R > S ***</td><td>1.87</td></tr></tbody></table>

This ranking confirmed the major role of some metabolites features, including sucrose, glucose, fructose, glutamate, amino acids and fumarate, in the response of susceptible accessions. Of the 140 VIP values identified by OPLS-DA, 119 accumulated to higher levels in susceptible accessions and 21, among which the four VIPs (VIP 433;700, VIP 361;491, VIP 512;491 and VIP 64;395) are found, accumulated to higher levels in resistant accessions (Table 2). The fold change susceptible/resistant reached up to 45 times whereas the fold change resistant/susceptible reached up to 244 times. It is very interesting to highlight two VIPs, VIP 324;184 and VIP 433;482 that are significantly accumulated in resistant accessions when inoculated with a fold-change of 106.8 and 244, respectively, compared to the susceptible accessions.

A complementary analysis, Partial Least Square (PLS), on the dataset obtained with the 26 accessions was performed to test whether viral accumulation can be predicted from metabolic data. Analysis was performed on the optical density values corresponding to the viral accumulation measured by DAS ELISA for each accession (Figure 5). The PLS coefficient estimated in the training data set revealed, after cross-validation, a correlation of 0.61 between predicted and true viral accumulation, confirming the high predictive power of metabolic composition for TuMV susceptibility.

<footnote n="1">¹ When undetermined, VIP are identified through mz/rt values. VIP values are classified in decreasing order,</footnote>

<footnote n="2">² mass to charge ratio, ³ retention time. ⁴ The significance was assessed through a Wilcoxon test at *** p < 0.001,</footnote>

<footnote>** 0.001 < p < 0.01.</footnote>
<figure><img src="image_10.png" /><figcaption>**Figure 5.** Prediction of viral accumulation by the metabolite matrix on 26 *A. thaliana* accessions in the '2015' field experiment.</figcaption></figure>

The metabolite matrix is composed by 10 primary metabolic traits and 505 metabolic signatures (m/z). The replicates of the 26 *A. thaliana* accessions are represented by the blue dots (total data point = 98). The dashed linear red line represents the exact prediction.

Among the best 50 VIP-PLS values, 44 were found common to the VIP values obtained with the OPLS-DA analysis (Table S4). In both analyses, the central metabolites sucrose, glucose and glutamate were found to discriminate susceptible and resistant accessions, ranking among the best VIP values (Table S4). Six common VIPs (VIP 394;108, VIP 351;184, VIP 324;184, VIP 280;184, VIP 432;184 and VIP 86;184) detected by the two analyses are found accumulated more significantly in resistant accessions.

**3. Discussion**

In this study, we characterized the metabolic response of *A. thaliana* to its natural viral pathogen, turnip mosaic virus (TuMV). This study is unusual in studying the metabolic response of a variety of accessions to an important naturally occurring virus on *A. thaliana* [42] and its close relative, *A. halleri* [43], in natural settings.

In neither year of the study, from time 0 to 13 days after inoculation, could we distinguish the metabolic profile of TuMV-susceptible and TuMV-resistant accessions grown in the absence of TuMV. This suggests that there are no constitutive metabolic patterns associated with resistance or susceptibility. Nevertheless, these mock-inoculated samples provided an opportunity to describe elements of central metabolism in *A. thaliana* under common garden field conditions. We found that dry biomass was significantly correlated with fumarate concentrations. Fumarate can accumulate to high levels in *A. thaliana* relative to other plant species, suggesting that it likely constitutes a significant fraction of the fixed carbon in *A. thaliana* rosette leaves [44]. Indeed, the amount of carbon stored in fumarate is similar to that accumulated in starch [44]. This is perhaps not surprising because fast-growing plant species such as *A. thaliana* contain significantly higher concentrations of organic acids, such as fumarate, compared to slow-growing plants [44,45], especially under high light intensity conditions such as in our field experiments [45]. We also found that dry biomass was correlated to glucose. Both fumarate and glucose reached high levels in fast growing accessions, 23 and 19 mM, respectively, suggesting that they might be involved in turgor dynamics and thus growth by cellular expansion [46]. In contrast, we detected a negative correlation between dry biomass and protein. Although speculative, this might follow from the fact that lower protein synthesis contributes to increased efficiency of carbon use, because protein synthesis is a costly process [47,48]. This negative relationship
is in agreement with previous observations among a collection of *A. thaliana* accessions grown under greenhouse/growth chamber and in controlled conditions [49].

Upon viral infection, in the field, the effect of TuMV inoculation could be observed as early as 5 days after inoculation. At 13 dai, susceptible accessions were clearly distinguished from the resistant ones by the accumulation of central metabolites. Eight out of ten metabolic features (amino acids, glutamate, malate, fumarate, starch, glucose, fructose and sucrose), were strongly associated with viral accumulation. Previous efforts to describe the massive reprogramming of the plant central metabolism that takes place in response to pathogens has focused on fungi and bacteria [50]. Here, we add to the limited literature on plant responses to viral infection, which is restricted to controlled laboratory conditions. For example, in a study of the central metabolic response of *Arabidopsis* to tobacco rattle virus (TRV), Fernandez-Calvino and collaborators [6] showed that the susceptible Col-0 accession significantly accumulated sucrose at 8 days' post inoculation. Amino acids were also globally accumulated in infected plants compared to mock plants whereas neither starch nor fumarate were accumulated in infected plants. Differences in the virus species, the number of accessions observed and the fact that theirs was a laboratory experiment could explain these different results. Second, in a controlled multi-stress experiment (including TuMV, drought and heat), accumulation of soluble sugars was again observed in Col-0 plants [51]. The strong metabolic disturbances generated in susceptible plants by TuMV infection suggests that the viral infection stimulated the accumulation of major central metabolites, probably as a result of an imbalance between photosynthesis and growth, which is known to favor oxidative stress [52]. Susceptible accessions were also found to accumulate a large number of specialized compounds. In addition to the fact that such an accumulation was apparently ineffective, the synthesis of specialized compounds is expensive in energy and, even if somewhat controlled [53], accumulation could have toxic effects.

Under the field conditions used in our experiments, inoculated resistant accessions were able to grow at the same rate as controls. As biomass is considered as appropriate proxy for fitness under many circumstances [54], it seems to indicate that resistance was achieved with no penalties on measured traits. Unlike susceptible accessions, resistant accessions accumulate a limited number of specialized metabolites in response to the virus, but in much higher quantities than can be observed in susceptible accessions. This contrasts with a study conducted under controlled conditions on the metabolic response of tomato to tomato yellow leaf curl virus (TYLCV) which found that the defense response of resistant lines is effective through the accumulation of many specialized metabolites [55]. The high predictive power of metabolic composition for TuMV susceptibility when infected might enable this composition to serve as a biomarker [30], as has been done for resistance to *Fusarium graminearum* in wheat [56] and for susceptibility to esca disease in grape [57].

The complex molecular network underlying the balance between growth and immunity has been described [24,25,58] mainly in terms of opposition. Recent reviews [26,59] have emphasized, however, that in a natural context, growth and immunity are in a constant conversation. Our study supports this alternative model and highlights new results on *Arabidopsis*/virus interaction. In particular, we have highlighted the very different metabolic paths between susceptible and resistant genotypes. We found that susceptible accessions experience a large accumulation of central and specialized metabolites with a reduction of growth whereas resistant accessions appear capable of continued growth with a targeted metabolic response. Some compounds as VIP 324;184 and VIP 433;482, that presented high fold-change in resistant accessions compared to susceptible ones are of particular interest. Preliminary putative annotation of major metabolic markers was performed using RT, accurate *m/z* detected by high-resolution MS and MS2 fragments as described previously [60,61]. The resulting predicted molecular formula were screened through chemical databases (HMDB, METLIN, MassBank) to match putative metabolite identification. The first results show that some of them are phenolic compounds including flavonoids and coumarins. The anti-viral properties of some of these compounds such as
quercetagetin or other flavonoids have been demonstrated against Tomato bushy stunt virus [62] and Tobacco mosaic virus [63].

It would be of great interest completing and refining the characterization of these
specialized compounds and their biosynthetic pathways, and then to test their involvement
in resistance by using pharmacological and/or genetic approaches. It is worth mentioning
that anti-phage specialized metabolites molecules, able to block phage replication, have
recently been found in *Streptomyces* [64]. This would allow for better understanding the
complexity of the underlying mechanisms involved in plants' responses to viruses in the
field and to propose new ideotypes.

**4. Materials and Methods**

4.1. *Plant Material*

Genotypes from a worldwide collection of natural accessions of *Arabidopsis thaliana*
were used in our experiments [40]. Twenty-six accessions were challenged with TuMV
during the first field experiment in 2015 and then a subset of 10 accessions were selected for
a second field experiment in 2017. This subset was selected based on their OD value follow-
ing TuMV inoculation (listed in Table 1 and Table S2) to represent extreme phenotypes from
highly susceptible (S) to resistant (R). To prevent misidentifying accessions as resistant due
to inefficient mechanical inoculation, only accessions presenting a resistance phenotype in
all our common garden experiments (see [39]) were selected as resistant accessions.

4.2. Virus Material

Turnip mosaic virus isolate UK1 [65] was routinely propagated by mechanical inocu-
lations on turnip plants *Brassica rapa* L. ssp *rapa* NA FR 490,001 provided by the BraCySol
germplasm center (Ploudaniel, France). To prepare the inoculum, three-week old turnip
plants were mechanically inoculated. Symptoms appeared two weeks later. Young symp-
tomatic leaves of five-week-old turnip were then collected to produce the inoculum.

4.3. *Experimental Design and Growth Conditions*

Two common garden experiments were conducted in 2015 (N = 26 *A. thaliana* acces-
sions) and 2017 (N = 10). Experiments realized in 2015 were organized in a randomized
complete block design (RCBD). The experiment performed in 2015 contained four blocks
with four replicates per accession per block. In 2017, a kinetic analysis was carried out.
Thus, five sampling were performed. The first sampling was done the day of inoculation
(named thereafter 0 days after inoculation (dai)), and at 5, 7, 9 and 13 (dai). The accessions
were arranged in a randomly complete block design of three blocks with three replicates
per accession per sampling date.

As described in [39], seedling trays of 40 wells were used. Seeds were sown on
23 March 2015 and on 20 March 2017 in professional horticultural soil (106Scope, PletrA-
com, Arles, France) under a cold-frame glasshouse without additional light or heating to
ensure homogeneity of germination. At three weeks of age, plants were acclimatized under
an opened tunnel before their transfer to the common garden. Plants were inoculated
during the acclimation step just before their transfer. Soil of the common garden had
been tilled so that seedling trays could be slightly buried. Because the bottoms of the
wells were pierced, roots were able to reach the soil. Climate data were recorded over
the duration of each experiment (Figure S3a,b). The analysis of climatic data across the
three growing environments (cold frame greenhouse, then tunnel, then field) revealed no
significant difference between 2015 and 2017 in the two first steps (cold frame greenhouse
and tunnel) but a significant difference between 2015 and 2017 for temperatures, rainfall
and PAR (photosynthetically active radiation) in common garden (Figure S3c-e).

4.4. *TuMV Inoculation Procedure and Harvest*

One gram of fresh turnip leaves was ground in three volumes of disodium phosphate (Na₂HPO₄·12H₂O) 30 mM and 0.2% diethyldithiocarbamic acid (DIECA). The inoculum
was clarified through 10 min centrifugation at 13,000× g. Supernatant was recovered and maintained at 4 °C until *Arabidopsis*' inoculation, which was performed when plants were 4 weeks old and at a 8–10 leaf stage, corresponding to 1.09, 1.10 Boyes stage [66]. Four young expanded leaves of each plant were mechanically inoculated with 20 µL of inoculum with carborundum added on each leaf. Ten minutes after inoculation, plants were rinsed with water. Mock treatments, in which plants were treated exactly as inoculated plants except for the absence of the virus, were also included. The viral concentration of the TuMV inoculum was quantified after inoculation by quantitative PCR [39]. The viral concentrations of the inocula used was 0.049 g/µL in 2015 and 0.269 ng/µL in 2017 common garden experiments (significantly different through a Wilcoxon-Mann Whitney test *p* = 0.029).

According to previous experiments [39], samples were collected 13 days after inocula-
tion (dpi) in the 2015 experiment at the very beginning of the onset of the first symptoms.
In 2017, the first sampling was done the day of inoculation (named thereafter 0 days after
inoculation (dai)), and at 5, 7, 9 and 13 (dai). The harvest was done in the morning, at
the same time for each year of experiment. All rosette leaves above the inoculated leaves
(systemic leaves) were collected in vials from Zinsser Analytic® (Eschborn, Germany).
Samples were deep-frozen, ground and stored at −80 °C.

4.5. *Quantification of Viral Accumulation*

Viral accumulation was estimated for each individual plant using 100 mg of *A. thaliana*
powder in a semi-quantitative double antibody sandwich assay (DAS-ELISA) with a
commercial anti-potyvirus monoclonal antibody kit (Agdia-Biofords, Evry, France). The
reaction of the substrate (*p*-nitrophenyl phosphate) was followed at 405 nm. Optical
densities (OD) were calculated by removing the mean OD value from the healthy *A. thaliana*
Col-0 control and normalized using a Col-0 positive control deposited on each ELISA plate.
As determined from the ODs of the *A. thaliana* Col-0 positive controls, and in order to
avoid overflow, the 15 min measurements were retained. For each accession, we classified
its degree of resistance or susceptibility based on the average of the ODs obtained across
replicates (Table 1; [39]). Resistant accessions do not accumulate virus and present an OD
equal or lower than the mean of the OD of the healthy *A. thaliana* Col-0 control. Susceptible
accessions have an OD value that exceeds the healthy *A. thaliana* Col-0 control. The category
to which each accession belongs is listed in Table 1 and Table S2, along with its OD values.
Mock-inoculated plants were confirmed to be potyvirus-free.

4.6. *Sample Processing*

A biological replicate was composed by pooling up to four (2015) or three (2017) plants after a successful infection as measured by DAS ELISA. Aliquots of about 17 mg of fresh powder were weighted in 1.1-mL Micronic tubes and used in the targeted metabolite analysis. Samples were then lyophilized and aboveground dry mass determined with an analysis balance. Aliquots of about 10 mg were weighed in 1.1-mL Micronic tubes and used in the untargeted metabolite analysis.

4.7. *Targeted Metabolite Analysis*

Metabolites were extracted in a final volume of 650µL, twice with 80% (v/v) ethanol—
HEPES/KOH 10mM (pH6) and once with 50% (v/v) ethanol—HEPES/KOH 10mM
(ph6) [67]. Chlorophyll content was determined immediately after the extraction [68].
Glucose, fructose, sucrose [69] malate, fumarate ([70], glutamate [71] and amino acids [72]
were determined in the supernatant. Starch [73] and protein [74] contents were determined
on the pellet resuspended in 100 mM NaOH. Analyses were performed in 96-well mi-
croplates using Starlet pipetting robots (Hamilton), and absorbance was read in MP96
microplate readers (SAFAS).
4.7.1. Data Treatments

Metabolite data were normalized by dividing each measure by the value of a biological
standard sample corresponding to a non-inoculated *A. thaliana* Col-0 sample harvested
at the same time as the other samples. The final concentration of each metabolite is the
average of the replicates of each accession.

4.7.2. Statistical Analysis

All principal component analysis (PCA), parametric and nonparametric statistical tests were performed using R version 3.4.0. Statistical significance was set at *p* < 0.05. Partial least squares regression (PLS) and orthogonal partial least-squares regression and discriminant analysis (OPLS-DA) were performed using the R packages *mixOmics* [75] and *PLS* [76].

4.8. Untargeted Metabolic Analysis

Untargeted metabolite measurements were conducted on the 26 accessions in the common garden experiment in 2015. Extraction was realized using a Starlet pipetting robot (Hamilton, Lancaster, PA, USA) with an extraction buffer composed of 80% (v/v) ethanol and 0.1% (v/v) formic acid, using methyl vanillate as an internal standard (50 µg/mL). All samples extracted were filtrated through a Multiscreen Solvinert 96-well filter plate (Merck Millipore, Burlington, MA, USA).

4.8.1. Quality Control (QC)

Twenty five µL of each sample were mixed together to generate a pooled quality control sample (QC). QC samples were analyzed every 10 injections to monitor and correct changes in the instrument response. Solvent blank samples (80% ethanol in water—0.1% formic acid) were also analyzed in-between the other samples.

4.8.2. Liquid Chromatography

Liquid chromatography was performed on a Dionex UHPLC Ultimate 3000 (Thermo Scientific, Waltham, MA, USA). Chromatographic separation was carried out in reverse-phase mode on a Gemini C18 column (2 × 150 mm; 3 µm, Phenomenex, Torrance, CA, USA) equipped with a Gemini C18 guard column (2 × 4 mm, Phenomenex). The mobile phase was composed of milliQ water with 0.1% formic acid (solvent A) and 100% Acetonitrile (solvent B) for a total run time of 18 min. The flow rate was 0.3 mL·min⁻¹ and the column was heated to 30 °C. The autosampler temperature was maintained at 4 °C and the injection volume was 5 µL.

4.8.3. Mass Spectrometry

The UHPLC system was coupled with a LTQ-Orbitrap Elite mass spectrometer (Thermo Scientific). A heated electrospray interface was used and analyses were performed in positive mode. Acquisition was performed in full scan mode with a resolution of 240,000 FWHM in the scan range of *m/z* 50–1000. Data were recorded using the Xcalibur software (Thermo Scientific) and extracted with XCMS.

4.8.4. Data Processing and Statistical Analysis

Data were converted to mzXML file format. Peak picking and alignment were per-
formed using XCMS in R [77]. XCMS-parameters were optimized as described by Patti
and collaborators [78]. XCMS-data processing results in a data matrix which contains
peak intensities that are a unique combination of retention-time and median m/z ratio. To
exclude system-peaks (impurities in the measurement-system, visible in blanks) as well as
poorly detected metabolic features, filter steps were performed based on QCs and blanks as
described in TextS1. Statistical analysis was done in R version 3.4.0 and also using the web
application BioStatFlow version 2.7.7 (http://biostatflow.org/ (accessed on 1 April 2021))
Thus, OPLS-DA (parameters: Kernel type linear, cross validation K-fold with K = 10,
permutations for testing 100) were performed using the 505 features of the untargeted
metabolic analysis combined with the 10 metabolites features assayed by targeted methods.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/metabo11040230/s1, Figure S1: Principal component analysis performed on the 505 metabolic signatures (m/z) measured on 26 *A. thaliana* accessions in the '2015' field experiment, Figure S2: OPLS-DA analysis and its parameters of validation for the TuMV-inoculated resistant and susceptible 26 *A. thaliana* accessions of the '2015' field experiment, Figure S3: Climate raw data and comparison between the 2015 and 2017 field experiment, Table S1: Spearman correlations between dry biomass (DB), viral accumulation (OD) and 10 primary metabolites content measured on the 26 *A. thaliana* accessions of the '2015' field experiment, Table S2: List of the *A. thaliana* accessions with their geographic position and their susceptible status to TuMV infection confirmed by OD values and SD in the '2017' field experiment, Table S3: Common Variable Importance in the Projection (VIPs) obtained with the OPLS-DA performed on mock and TuMV-inoculated and on resistant and susceptible accessions of twenty-six accessions in 2015, Table S4: Comparisons between the Variable Importance in the Projection (VIPs) identify by OPLS-DA analysis and those identify by PLS analysis performed with resistant and susceptible twenty-six accessions in 2015.

**Author Contributions:** Conceptualization, Y.G. and V.S.-L.; Data curation, B.R., O.F. and T.B.; Formal analysis, O.F., P.C., T.B. and V.S.-L.; Funding acquisition, Y.G. and V.S.-L.; Investigation, B.R., M.C., R.L. and V.S.-L.; Methodology, B.R., O.F., P.C. and V.S.-L.; Supervision, Y.G. and V.S.-L.; Writing—original draft, B.R. and V.S.-L.; Writing—review & editing, B.R., O.F., F.R., J.B., Y.G. and V.S.-L. All the authors proofread the final version of the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** The authors thank PHENOME (ANR-11-INBS-0012) and MetaboHUB (ANR-11-INBS-0010) for financing. VS-L thanks the Fulbright Scholarship Program for accommodation and research support in Joy Bergelson's laboratory.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Publicly available datasets were analyzed in this study. This data can be found here: https://data.inrae.fr/privateurl.xhtml?token=7ddc25dc-e76c-4ac4-9acc-61d89e64076f (accessed on 1 April 2021).

**Acknowledgments:** The authors are grateful to Thierry Mauduit and Pascal Briard for their help in the management of field experiments.

**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. Anderson, P.K.; Cunningham, A.A.; Patel, N.G.; Morales, F.J.; Epstein, P.R.; Daszak, P. Emerging Infectious Diseases of Plants: Pathogen Pollution, Climate Change and Agrotechnology Drivers. *Trends Ecol. Evol.* **2004**, *19*, 535–544. [CrossRef]

2. Bebber, D.P.; Ramotowski, M.A.T.; Gurr, S.J. Crop Pests and Pathogens Move Polewards in a Warming World. *Nat. Clim. Chang.* **2013**, *3*, 985–988. [CrossRef]

3. Jones, R.A.C. Future Scenarios for Plant Virus Pathogens as Climate Change Progresses. *Adv. Virus Res.* **2016**, *95*, 87–147. [CrossRef]

4. Adams, M.J.; King, A.M.; Lefkowitz, E.; Carstens, E.B. Part II: The viruses—Family Potyviridae. In *Virus Taxonomy—Ninth Rep. Int. Committee on Taxonomy of Viruses*; Elsevier: San Diego, CA, USA, 2011; pp. 1069–1090.

5. Stare, T.; Ramšak, Ž.; Blejec, A.; Stare, K.; Turnšek, N.; Weckwerth, W.; Wienkoop, S.; Vodnik, D.; Gruden, K. Bimodal Dynamics of Primary Metabolism-Related Responses in Tolerant Potato-Potato Virus Y Interaction. *BMC Genom.* **2015**, *16*, 1–17. [CrossRef]

6. Fernandez-Calvino, L.; Osorio, S.; Hernandez, M.L.; Hamada, I.B.; Del Toro, F.J.; Donaire, L.; Yu, A.; Bustos, R.; Fernie, A.R.; Martinez-Rivas, J.M.; et al. Virus-Induced Alterations in Primary Metabolism Modulate Susceptibility to Tobacco Rattle Virus in Arabidopsis. *Plant Physiol.* **2014**, *166*, 1831–1838. [CrossRef] [PubMed]

7. Biemelt, S.; Sonnewald, U. Plant-Microbe Interactions to Probe Regulation of Plant Carbon Metabolism. *J. Plant Physiol.* **2006**, *163*, 307–318. [CrossRef] [PubMed]
8. Shalitin, D.; Wang, Y.; Omid, A.; Gal-On, A.; Wolf, S. Cucumber Mosaic Virus Movement Protein Affects Sugar Metabolism and Transport in Tobacco and Melon Plants. *Plant Cell Environ.* **2002**, *25*, 989–997. [CrossRef]

9. Kangasjärvi, S.; Neukermans, J.; Li, S.; Aro, E.-M.; Noctor, G. Photosynthesis, Photorespiration, and Light Signalling in Defence Responses. *J. Exp. Bot.* **2012**, *63*, 1619–1636. [CrossRef] [PubMed]

10. Handford, M.G.; Carr, J.P. A Defect in Carbohydrate Metabolism Ameliorates Symptom Severity in Virus-Infected Arabidopsis Thaliana. *J. Gen. Virol.* **2007**, *88*, 337–341. [CrossRef] [PubMed]

11. Kogovsek, P.; Pompe-Novak, M.; Petek, M.; Fragner, L.; Weckwerth, W.; Gruden, K. Primary Metabolism, Phenylpropanoids and Antioxidant Pathways Are Regulated in Potato as a Response to Potato Virus Y Infection. *PLoS ONE* **2016**, *11*, e0146135. [CrossRef]

12. Tecsi, L.I.; Smith, A.M.; Maule, A.J.; Leegood, R.C. A Spatial Analysis of Physiological Changes Associated with Infection of Cotyledons of Marrow Plants with Cucumber Mosaic Virus. *Plant Physiol.* **1996**, *111*, 975–985. [CrossRef] [PubMed]

13. Llave, C. Dynamic Cross-Talk between Host Primary Metabolism and Viruses during Infections in Plants. *Curr. Opin. Virol.* **2016**, *19*, 50–55. [CrossRef] [PubMed]

14. Pilar Lopez-Gresa, M.; Lison, P.; Kim, H.K.; Choi, Y.H.; Verpoorte, R.; Rodrigo, I.; Conejero, V.; Maria Belles, J. Metabolic Fingerprinting of Tomato Mosaic Virus Infected *Solanum lycopersicum*. *J. Plant Physiol.* **2012**, *169*, 1586–1596. [CrossRef]

15. Heaton, N.S.; Randall, G. Multifaceted Roles for Lipids in Viral Infection. *Trends Microbiol.* **2011**, *19*, 368–375. [CrossRef]

16. Wang, A. Dissecting the Molecular Network of Virus-Plant Interactions: The Complex Roles of Host Factors. *Annu. Rev. Phytopathol.* **2015**, *53*, 45–66. [CrossRef]

17. Piasecka, A.; Jedrzejczak-Rey, N.; Bednarek, P. Secondary Metabolites in Plant Innate Immunity: Conserved Function of Divergent Chemicals. *New Phytol.* **2015**, *206*, 948–964. [CrossRef] [PubMed]

18. Verma, V.; Ravindran, P.; Kumar, P.P. Plant Hormone-Mediated Regulation of Stress Responses. *BMC Plant Biol.* **2016**, *16*, 1–10. [CrossRef]

19. Palukaitis, P.; Yoon, J.-Y.; Choi, S.-K.; Carr, J.P. Manipulation of Induced Resistance to Viruses. *Curr. Opin. Virol.* **2017**, *26*, 141–148. [CrossRef] [PubMed]

20. Kliebenstein, D.J. Secondary Metabolites and Plant/Environment Interactions: A View through Arabidopsis Thaliana Tinged Glasses. *Plant Cell Environ.* **2004**, *27*, 675–684. [CrossRef]

21. Tian, D.; Traw, M.B.; Chen, J.Q.; Kreitman, M.; Bergelson, J. Fitness Costs of R-Gene-Mediated Resistance in Arabidopsis Thaliana. *Nature* **2003**, *423*, 74–77. [CrossRef]

22. Denance, N.; Sanchez-Vallet, A.; Goffner, D.; Molina, A. Disease Resistance or Growth: The Role of Plant Hormones in Balancing Immune Responses and Fitness Costs. *Front. Plant Sci.* **2013**, *4*, 155. [CrossRef]

23. Neilson, E.H.; Goodger, J.Q.D.; Woodrow, I.E.; Moller, B.L. Plant Chemical Defense: At What Cost? *Trends Plant Sci.* **2013**, *18*, 250–258. [CrossRef]

24. Huot, B.; Yao, J.; Montgomery, B.L.; He, S.Y. Growth-Defense Tradeoffs in Plants: A Balancing Act to Optimize Fitness. *Mol. Plant* **2014**, *7*, 1267–1287. [CrossRef] [PubMed]

25. Lozano-Duran, R.; Zipfel, C. Trade-off between Growth and Immunity: Role of Brassinosteroids. *Trends Plant Sci.* **2015**, *20*, 12–19. [CrossRef] [PubMed]

26. Karasov, T.L.; Chae, E.; Herman, J.J.; Bergelson, J. Mechanisms to Mitigate the Trade-Off between Growth and Defense. *Plant Cell* **2017**, *29*, 666–680. [CrossRef] [PubMed]

27. Albrecht, T.; Argueso, C.T. Should I Fight or Should I Grow Now? The Role of Cytokinins in Plant Growth and Immunity and in the Growth-Defence Trade-Off. *Ann. Bot.* **2017**, *119*, 725–735. [CrossRef]

28. Bergelson, J.; Purrington, C.B. Surveying Patterns in the Cost of Resistance in Plants. *Am. Nat.* **1996**, *148*, 536–558. [CrossRef]

29. Brown, J.K.M. Yield Penalties of Disease Resistance in Crops. *Curr. Opin. Plant Biol.* **2002**, *5*, 339–344. [CrossRef]

30. Fernandez, O.; Urrutia, M.; Bernillon, S.; Giauffret, C.; Tardieu, F.; Le Gouis, J.; Langlade, N.; Charcosset, A.; Moing, A.; Gibon, Y. Fortune Telling: Metabolic Markers of Plant Performance. *Metabolomics* **2016**, *12*, 158. [CrossRef]

31. Schauer, N.; Semel, Y.; Roessner, U.; Gur, A.; Balbo, I.; Carrari, F.; Pleban, T.; Perez-Melis, A.; Bruedigam, C.; Kopka, J.; et al. Comprehensive Metabolic Profiling and Phenotyping of Interspecific Introgression Lines for Tomato Improvement. *Nat. Biotechnol.* **2006**, *24*, 447–454. [CrossRef]

32. Obata, T.; Witt, S.; Lisec, J.; Palacios-Rojas, N.; Florez-Sarasa, I.; Yousfi, S.; Araus, J.L.; Cairns, J.E.; Fernie, A.R. Metabolite Profiles of Maize Leaves in Drought, Heat, and Combined Stress Field Trials Reveal the Relationship between Metabolism and Grain Yield. *Plant Physiol.* **2015**, *169*, 2665–2683. [CrossRef]

33. Melandri, G.; AbdElgawad, H.; Riewe, D.; Hageman, J.A.; Asard, H.; Beemster, G.T.S.; Kadam, N.; Jagadish, K.; Altmann, T.; Ruyter-Spira, C.; et al. Biomarkers for Grain Yield Stability in Rice under Drought Stress. *J. Exp. Bot.* **2020**, *71*, 669–683. [CrossRef]

34. Meyer, R.C.; Steinfath, M.; Lisec, J.; Becher, M.; Witucka-Wall, H.; Torjek, O.; Fiehn, O.; Eckardt, A.; Willmitzer, L.; Selbig, J.; et al. The Metabolic Signature Related to High Plant Growth Rate in Arabidopsis Thaliana. *Proc. Natl. Acad. Sci. USA* **2007**, *104*, 4759–4764. [CrossRef] [PubMed]

35. Steinfath, M.; Strehmel, N.; Peters, R.; Schauer, N.; Groth, D.; Hummel, J.; Steup, M.; Selbig, J.; Kopka, J.; Geigenberger, P.; et al. Discovering Plant Metabolic Biomarkers for Phenotype Prediction Using an Untargeted Approach. *Plant Biotechnol. J.* **2010**, *8*, 900–911. [CrossRef] [PubMed]
36. Korn, M.; Gärtner, T.; Erban, A.; Kopka, J.; Selbig, J.; Hincha, D.K. Predicting Arabidopsis Freezing Tolerance and Heterosis in Freezing Tolerance from Metabolite Composition. *Mol. Plant* **2010**, *3*, 224–235. [CrossRef]

37. Züst, T.; Heichinger, C.; Grossniklaus, U.; Harrington, R.; Kliebenstein, D.J.; Turnbull, L.A. Natural Enemies Drive Geographic Variation in Plant Defenses. *Science* **2012**, *338*, 116. [CrossRef] [PubMed]

38. Bergelson, J.; Roux, F. Towards Identifying Genes Underlying Ecologically Relevant Traits in Arabidopsis Thaliana. *Nat. Rev. Genet* **2010**, *11*, 867–879. [CrossRef] [PubMed]

39. Rubio, B.; Cosson, P.; Caballero, M.; Revers, F.; Bergelson, J.; Roux, F.; Schurdi-Levraud, V. Genome Wide Association Study Reveals New Loci Involved in Arabidopsis Thaliana and Turnip Mosaic Virus (TuMV) Interactions in the Field. *New Phytologist* **2019**, *221*, 2026–2038. [CrossRef] [PubMed]

40. Horton, M.W.; Hancock, A.M.; Huang, Y.S.; Toomajian, C.; Atwell, S.; Auton, A.; Muliyati, N.W.; Platt, A.; Sperone, F.G.; Vilhjalmsson, B.J.; et al. Genome-Wide Patterns of Genetic Variation in Worldwide Arabidopsis Thaliana Accessions from the RegMap Panel. *Nat. Genet* **2012**, *44*, 212–216. [CrossRef] [PubMed]

41. Sutula, C.L.; Gillett, J.M.; Morrissey, S.M.; Ramsdell, D.C. Interpreting ELISA Data and Establishing the Positive-Negative Threshold. *Plant Dis.* **1986**, *70*, 722. [CrossRef]

42. Pagan, I.; Fraile, A.; Fernandez-Fueyo, E.; Montes, N.; Alonso-Blanco, C.; Garcia-Arenal, F. Arabidopsis Thaliana as a Model for the Study of Plant-Virus Co-Evolution. *Philos. Trans. R. Soc. Lond. B Biol. Sci.* **2010**, *365*, 1983–1995. [CrossRef] [PubMed]

43. Kamitani, M.; Nagano, A.J.; Honjo, M.N.; Kudoh, H. RNA-Seq Reveals Virus–Virus and Virus–Plant Interactions in Nature. *FEMS Microbiol. Ecol.* **2016**, *92*. [CrossRef] [PubMed]

44. Araujo, W.L.; Nunes-Nesi, A.; Fernie, A.R. Fumarate: Multiple Functions of a Simple Metabolite. *Phytochemistry* **2011**, *72*, 838–843. [CrossRef] [PubMed]

45. Chia, D.W.; Yoder, T.J.; Reiter, W.-D.; Gibson, S.I. Fumaric Acid: An Overlooked Form of Fixed Carbon in Arabidopsis and Other Plant Species. *Planta* **2000**, *211*, 743–751. [CrossRef] [PubMed]

46. Fricke, W. Turgor Pressure. In *eLS*; American Cancer Society, John Wiley & Sons, Ltd.: Weinheim, Germany, 2017; pp. 1–6, ISBN 978-0-470-01590-2.

47. Sajitz-Hermstein, M.; Nikoloski, Z. A Novel Approach for Determining Environment-Specific Protein Costs: The Case of Arabidopsis Thaliana. *Bioinformatics* **2010**, *26*, i582–i588. [CrossRef]

48. Kafri, M.; Metzl-Raz, E.; Jona, G.; Barkai, N. The Cost of Protein Production. *Cell Rep.* **2015**, *14*, 22–31. [CrossRef] [PubMed]

49. Sulpice, R.; Nikoloski, Z.; Tschoep, H.; Antonio, C.; Kleessen, S.; Larhlimi, A.; Selbig, J.; Ishihara, H.; Gibon, Y.; Fernie, A.R.; et al. Impact of the Carbon and Nitrogen Supply on Relationships and Connectivity between Metabolism and Biomass in a Broad Panel of Arabidopsis Accessions. *Plant Physiol.* **2013**, *162*, 347–363. [CrossRef] [PubMed]

50. Rojas, C.M.; Senthil-Kumar, M.; Tzin, V.; Mysore, K.S. Regulation of Primary Plant Metabolism during Plant-Pathogen Interactions and Its Contribution to Plant Defense. *Front. Plant Sci.* **2014**, *5*, 17. [CrossRef]

51. Prasch, C.M.; Sonnewald, U. Simultaneous Application of Heat, Drought, and Virus to Arabidopsis Plants Reveals Significant Shifts in Signaling Networks. *Plant Physiol.* **2013**, *162*, 1849–1866. [CrossRef] [PubMed]

52. Foyer, C.H.; Shigeoka, S. Understanding Oxidative Stress and Antioxidant Functions to Enhance Photosynthesis. *Plant Physiol.* **2011**, *155*, 93–100. [CrossRef] [PubMed]

53. Sirikantaramas, S.; Yamazaki, M.; Saito, K. Mechanisms of Resistance to Self-Produced Toxic Secondary Metabolites in Plants. *Phytochem. Rev.* **2007**, *7*, 467. [CrossRef]

54. Younginger, B.S.; Sirová, D.; Cruzan, M.B.; Ballhorn, D.J. Is Biomass a Reliable Estimate of Plant Fitness? *Appl. Plant Sci.* **2017**, *5*, [CrossRef]

55. Sade, D.; Shriki, O.; Cuadros-Inostroza, A.; Tohge, T.; Semel, Y.; Haviv, Y.; Willmitzer, L.; Fernie, A.R.; Czosnek, H.; Brotman, Y. Comparative Metabolomics and Transcriptomics of Plant Response to Tomato Yellow Leaf Curl Virus Infection in Resistant and Susceptible Tomato Cultivars. *Metabolomics* **2015**, *11*, 81–97. [CrossRef]

56. Cuperlovic-Culf, M.; Wang, L.; Forseille, L.; Boyle, K.; Merkley, N.; Burton, I.; Fobert, P.R. Metabolic Biomarker Panels of Response to Fusarium Head Blight Infection in Different Wheat Varieties. *PLoS ONE* **2016**, *11*, e0153642. [CrossRef]

57. Felgueiras, M.L.; Lima, M.R.M.; Dias, A.C.P.; Gil, A.M.; Graça, G.; Rodrigues, J.E.A.; Barros, A. NMR Metabolomics of Esca Disease-Affected Vitis Vinifera Cv. Alvarinho Leaves. *J. Exp. Bot.* **2010**, *61*, 4033–4042. [CrossRef]

58. Chandran, D.; Rickert, J.; Huang, Y.; Steinwand, M.A.; Marr, S.K.; Wildermuth, M.C. Atypical E2F Transcriptional Repressor DEL1 Acts at the Intersection of Plant Growth and Immunity by Controlling the Hormone Salicylic Acid. *Cell Host Microbe* **2014**, *15*, 506–513. [CrossRef] [PubMed]

59. Kliebenstein, D.J. False Idolatry of the Mythical Growth versus Immunity Tradeoff in Molecular Systems Plant Pathology. *Physiol. Mol. Plant Pathol.* **2016**, *95*, 55–59. [CrossRef]

60. Musseau, C.; Jorly, J.; Gadin, S.; Sørensen, I.; Deborde, C.; Bernillon, S.; Mauxion, J.-P.; Atienza, I.; Moing, A.; Lemaire-Chamley, M.; et al. The Tomato Guanylate-Binding Protein SIGBP1 Enables Fruit Tissue Differentiation by Maintaining Endopolyploid Cells in a Non-Proliferative State. *Plant Cell* **2020**, *32*, 3188. [CrossRef] [PubMed]

61. Luna, E.; Flandin, A.; Cassan, C.; Prigent, S.; Chevanne, C.; Kadiri, C.F.; Gibon, Y.; Pétriacq, P. Metabolomics to Exploit the Primed Immune System of Tomato Fruit. *Metabolites* **2020**, *10*, 96. [CrossRef]

62. Rusak, G.; Krajačić, M.; Pleše, N. Inhibition of Tomato Bushy Stunt Virus Infection Using a Quercetagetin Flavonoid Isolated from *Centaurea rupestris* L. *Antivir. Res.* **1997**, *36*, 125–129. [CrossRef]
63. Zhao, W.; Zeng, X.; Zhang, T.; Wang, L.; Yang, G.; Chen, Y.-K.; Hu, Q.; Miao, M. Flavonoids from the Bark and Stems of Cassia Fistula and Their Anti-Tobacco Mosaic Virus Activities. *Phytochem. Lett.* **2013**, *6*, 179–182. [CrossRef]

64. Kronheim, S.; Daniel-Ivad, M.; Duan, Z.; Hwang, S.; Wong, A.I.; Mantel, I.; Nodwell, J.R.; Maxwell, K.L. A Chemical Defence against Phage Infection. *Nature* **2018**, *564*, 283–286. [CrossRef]

65. Jenner, C.E.; Walsh, J.A. Pathotypic Variation in Turnip Mosaic Virus with Special Reference to European Isolates. *Plant Pathol.* **1996**, *45*, 848–856. [CrossRef]

66. Boyes, D.C.; Zayed, A.M.; Ascenzi, R.; McCaskill, A.J.; Hoffman, N.E.; Davis, K.R.; Gorlach, J. Growth Stage-Based Phenotypic Analysis of Arabidopsis: A Model for High Throughput Functional Genomics in Plants. *Plant Cell* **2001**, *13*, 1499–1510. [CrossRef] [PubMed]

67. Geigenberger, P.; Lerchi, J.; Stitt, M.; Sonnewald, U. Phloem-Specific Expression of Pyrophosphatase Inhibits Long Distance Transport of Carbohydrates and Amino Acids in Tobacco Plants. *Plant Cell Environ.* **1996**, *19*, 43–55. [CrossRef]

68. Arnon, D.I. Copper Enzymes in Isolated Chloroplasts. Polyphenoloxidase in Beta Vulgaris. *Plant Physiol.* **1949**, *24*, 1. [CrossRef] [PubMed]

69. Stitt, M.; Quick, P. Photosynthetic Carbon Partitioning: Its Regulation and Possibilities for Manipulation. *Physiol. Plant.* **1989**, *77*, 633–641. [CrossRef]

70. Nunes-Nesi, A.; Carrari, F.; Gibon, Y.; Sulpice, R.; Lytovchenko, A.; Fisahn, J.; Graham, J.; Ratcliffe, R.G.; Sweetlove, L.J.; Fernie, A.R. Deficiency of Mitochondrial Fumarase Activity in Tomato Plants Impairs Photosynthesis via an Effect on Stomatal Function: Inhibition of Fumarase Impairs Photosynthesis. *Plant J.* **2007**, *50*, 1093–1106. [CrossRef]

71. Zhang, N.; Gibon, Y.; Wallace, J.G.; Lepak, N.; Li, P.; Dedow, L.; Chen, C.; So, Y.-S.; Kremling, K.; Bradbury, P.J.; et al. Genome-Wide Association of Carbon and Nitrogen Metabolism in the Maize Nested Association Mapping Population. *Plant Physiol.* **2015**, *168*, 575–583. [CrossRef] [PubMed]

72. Bantan-Polak, T.; Kassai, M.; Grant, K.B. A Comparison of Fluorescamine and Naphthalene-2,3-Dicarboxaldehyde Fluorogenic Reagents for Microplate-Based Detection of Amino Acids. *Anal. Biochem.* **2001**, *297*, 128–136. [CrossRef]

73. Hendriks, J.H.; Kolbe, A.; Gibon, Y.; Stitt, M.; Geigenberger, P. ADP-Glucose Pyrophosphorylase Is Activated by Posttranslational Redox-Modification in Response to Light and to Sugars in Leaves of Arabidopsis and Other Plant Species. *Plant Physiol.* **2003**, *133*, 838–849. [CrossRef] [PubMed]

74. Bradford, M.M. A Rapid and Sensitive Method for the Quantitation of Microgram Quantities of Protein Utilizing the Principle of Protein-Dye Binding. *Anal. Biochem.* **1976**, *72*, 248–254. [CrossRef]

75. Rohart, F.; Gautier, B.; Singh, A.; Le Cao, K.-A. MixOmics: An R Package for omics Feature Selection and Multiple Data Integration. *PLoS Comput. Biol.* **2017**, *13*, e1005752. [CrossRef] [PubMed]

76. Wehrens, R.; Mevik, B.-H. The Pls Package: Principal Component and Partial Least Squares Regression in R. *J. Stat. Softw.* **2007**, *18*, 1–23.

77. Smith, C.A.; Want, E.J.; O'Maille, G.; Abagyan, R.; Siuzdak, G. XCMS: Processing Mass Spectrometry Data for Metabolite Profiling Using Nonlinear Peak Alignment, Matching, and Identification. *Anal. Chem.* **2006**, *78*, 779–787. [CrossRef] [PubMed]

78. Patti, G.J.; Tautenhahn, R.; Siuzdak, G. Meta-Analysis of Untargeted Metabolomic Data from Multiple Profiling Experiments. *Nat. Protoc.* **2012**, *7*, 508–516. [CrossRef]