Metabolic Profiling of Phloem Exudates as a Tool to Improve Bread-Wheat Cultivars

In a proof of concept study aimed at showing that metabolites in bread wheat (Triticum aestivum L. ssp. aestivum), phloem exudates have potential as biochemical markers for cultivar discrimination, Argentinean cultivars from three quality groups (groups 1, 2, and 3 of high, intermediate, and low quality, respectively) were grown under two nitrogen (N) availabilities and analysed for metabolic profile by electrospray ionisation mass spectrometry. Data as signal strengths of mass/charge (m/z) values binned to a resolution of 0.2 Daltons were subjected to principal component analysis and orthogonal projections to latent structures discriminant analysis. Certain bins were influential in discriminating groups taken in pairs and some were involved in separating all three groups. In high N availability, group 3 cultivars clustered away from the other cultivars, while group 1 cultivars clustered tightly together; group 2 cultivars were more scattered between group 1 and group 3 cultivars. In low N availability, the cultivars were not clustered as tightly; nonetheless, group 1 cultivars tended to cluster together and mainly separated from those of group 2. m/z values also showed potential for discrimination between N availability. In conclusion, phloem exudate metabolic profiles could provide biochemical markers for selection during breeding and for discerning the effects of N fertiliser application.


Introduction
The aim of the work reported in this paper was to examine whether metabolic profiling of wheat phloem exudates could be used to discriminate between cultivars with different bread-making qualities and to discern the effects of nitrogen (N) fertiliser application.Since the development of the wheat grain depends on the supply of nutrients in the phloem [1], such measurements could potentially also provide insights into aspects such as N use efficiency.The identification of biochemical markers useful for discriminatory purposes and for selecting combinations of good quality, high yield, and good tolerance to biotic and abiotic stresses appears feasible [2].Metabolic profiling has been used to distinguish between leaf and fruit extracts of species and cultivars of tomato and examine the effects of N on metabolism in these organs [3][4][5].
In bread wheat (Triticum aestivum L. ssp.aestivum), one important target trait for improvement is grain protein content (GPC), which is one of the determinants of international market price [1].Argentinean wheat cultivars are currently classified into three quality groups for bread-making based upon a range of tests, including grain protein content, wet gluten content, properties in the alveograph and farinograph, and loaf volume.Group 1 is the best quality (including wheats that can be blended with others to correct their visco-elastic properties and that are suitable for industrial bread-making), and group 3, the poorest quality, which tends to include high-yielding cultivars of deficient quality for bread-making (suitable for only short fermentation times, of less than eight hours).Group 2 is of intermediate quality, including cultivars that are not correctors and that are suitable for traditional bread-making and for fermentation of more than eight hours.Groups 1, 2, and 3 are expected to rank from high to low for protein content and other quality traits [6], and phloem composition would be expected to be related to such differences.
Soil or foliar applications of N, applied at different rates and stages of growth, are commonly used to improve GPC.One of the principal contributors to GPC is the process of remobilisation from leaves during organ senescence, but N runoff from soils creates environmental issues.Barneix [1] demonstrated that more than 50-70% of final grain N is accumulated by the plant before anthesis and is remobilised to the grain later.The relationship between the N supplied and that absorbed by the plant is not linear; rather, there is a limit to the potential GPC, which lowers fertiliser use efficiency when high N doses are supplied.The concentration of free amino acids in the phloem acts as a signal to the roots that indicates the N status of the plant, which activates or inhibits NO 3 uptake by the plant.As a consequence, N metabolism in the shoot dictates the rate of NO 3 uptake [1].The concentration of the majority of the amino acids in phloem exudates is proportional to the concentration in the leaves [7], if no pathogens are involved [8].
The concentration of the majority of the amino acids in phloem exudates is proportional to the concentration in the leaves [7], if no pathogens are involved [8].Grain filling is mainly dependent on remobilisation from the flag leaf and the adjacent leaf.Furthermore, the final GPC has been correlated with the amount of free amino acids in the flag leaf during grain filling [9,10].Therefore, the analysis of the metabolites present in the phloem exudates provides a valid indicator of the compounds that will be present in the future grains and of the effect of N fertilisation, given that secondary metabolites serve as a N reserve [11].
These relationships, showing the potential importance for the quality of phloem exudate composition, gave rise to the current work.Since direct injection mass spectrometry (DIMS) through electrospray ionisation time-of-flight (ESI-TOF-MS) provides a rapid method to obtain an initial metabolic profile of samples [3,12], it was chosen as an approach for profile analysis of phloem exudate samples in this proof of concept study.

Experimental Design
The field trial was carried out at the Experimental Field of the Faculty of Agronomy, UNCPBA (Lat.: 36 The trial was a randomised complete block design with three replicate blocks, with a split plot design where the main plot was N fertiliser level and the sub-plot was cultivar.Eight cultivars of contrasting quality were sown: 'ACA 304', 'Klein Proteo', and 'Buck Meteoro' of group 1 (high bread-making quality); 'Bio 3003', 'Bio 1000', and 'Buck Malevo' of group 2 (medium bread-making quality); and 'Klein Gavilán' and 'Klein Guerrero' of group 3 (low bread-making quality).These cultivars have been previously evaluated for quality traits [13][14][15][16][17][18][19], amongst others.
Soil N content before fertilisation was approximately 35 kg N per hectare and the field area was fertilised with triple superphosphate (100 ppm) prior to sowing.Sowing date was adjusted according to the cultivar cycle (short or long) with the aim of ensuring that the cultivars arrived simultaneously at the required physiological stage for extracting phloem exudates.Plots measured 9.5 m × 2.8 m and sowing density was 350 plants/m 2 .Half of each plot, randomly designated, was fertilised to achieve 100 ppm of N availability as urea by hand broadcasting.The herbicide DICAMBA-SPA was applied for weed control.At exudate sampling the plants showed good general health status, with differences in colour intensity between the well-fertilised treatments compared to the rest.

Extraction Method
Phloem exudate extraction was done according to Urquhart and Joy [20] as modified in [21].Samples were taken when the flag leaf was fully expanded, which occurred within a lapse of eight days covering all cultivars, meaning the different sowing dates for short and long cycled cultivars resulted in the desired limited range of dates for the expanded leaf physiological stage, minimising the effect of the cycle.
The exudates should mainly reveal the contents of the phloem; while we cannot rule out the possibility of xylem contents and cellular leakage also occurring in the samples, such contamination would be expected to be slight [20,22].The flag leaf and the adjacent leaf were harvested according to [21].The extremes were cut and complete leaves kept for 15 min in a 20 mm pH 6 ethylenediaminetetraacetic acid (EDTA) solution, which was subsequently discarded to minimise such contamination.Leaves were washed in new EDTA solution for one minute and then left for three hours in a new aliquot of the EDTA solution to allow the phloem to exude.The whole procedure was carried out in the dark to minimise leaf transpiration and absorption of the solution by the xylem.The exudates were stored at −4 • C and lyophilised to −40 • C and 76 mm Hg for 24 h.

Sample Processing
Samples were resuspended in 1 mL of H 2 O and an aliquot of 20 µL was diluted in 100 µL of methanol, 0.1 µL of formic acid, and distilled water to a final volume of 200 µL.The samples were analysed by direct injection mass spectrometry (DIMS) through electrospray ionisation time-of-flight (ESI-TOF-MS) using an Applied Systems/MDS Sciex (Foster City, CA, USA) hybrid quadrupole time-of-flight Q-Star Pulsar-I mass spectrometer.For all samples three technical replicates were analysed.
Of the 48 exudate samples (8 cultivars × 3 blocks × 2 N levels), those from ACA 304 without added N (−N) block II, Bio 3003 −N block I, Bio 1000 with added N (+N) block III and Meteoro +N block III went missing during processing.

Data Analysis
The mass spectra were processed according to Overy et al. [3].Only mass/charge (m/z) values occurring in all three technical replications were included and the data was binned to a resolution of 0.2 Daltons.The metabolic profiles thus obtained were then analysed by unsupervised statistics using principal component analysis (PCA) with SIMCA14 (Umetrics, Umeå, Sweden).When clustering was observed the data was processed by orthogonal projections to latent structures discriminant analysis (OPLS-DA) to resolve the bins contributing towards the clustering of samples and putative metabolites were assigned to these bins from empirical formulae.
Analysis of variance (ANOVA) was carried out using the software Infogen, Córdoba, Argentina [23], for individual bins between quality groups and between N availability.

Resolution of Cultivars
An overall PCA of all data (all samples including all detected masses) showed that the proportion of the total variation accounted for by principal components 1 and 2 was 56.5%.Since the analytical approach used was completely untargeted, the data were examined to reveal any clustering, by using an unsupervised statistical approach and then overlaying the known variables on the plot.Initially the data for the plots with and without added N were analysed separately.From the analysis of plots with added N (Figure 1a), it can be seen that the group 3 cultivars Gavilán and Guerrero cluster away from the other cultivars, while the group 1 cultivars ACA 304, Proteo, and Meteoro cluster tightly together.The group 2 cultivars are more scattered between the group 1 and group 3 cultivars, as might be expected from their intermediate characteristics.For the plots without added N (Figure 1b), the cultivars are not clustered as tightly.The group 3 variety replicates are more scattered, but the group 1 cultivars tend to cluster in the right hand two quartiles and, apart from one Proteo replicate, are separated from the group 2 cultivars.If fewer cultivars are included in the analysis to reduce the total variance, the groups are more clearly separated (data not shown).Since there is evidence of clustering of the data it was examined with an OPLS-DA analysis to identify those bins that contribute most towards the separation.These values are obtained from the loadings plot of the first predictive component of the model and correspond to the covariance of the X-variables and the predictive score vectors.It is interesting to note that many of the bins that separate group 3 from group 2 also separate group 3 from group 1 (Table 1).The cut-off point of the top ten values was an arbitrary decision to illustrate that some bins contribute to the separation of all three groups, such as bin 383; putative metabolites were assigned to these influential bins (Supplementary Table S1).
ANOVAs carried out for individual bins of high influence in the original overall analysis also showed differences between quality groups (Table 2).For example, mean values for bin 203, given in Table 1 as being influential in separating group 3 from the other groups, ranked group 1 < group 2 < group 3 in the ANOVAs, and so may be associated with poor quality.In contrast, bin 370.6 ranked group 1 > group 2 > group 3, and so may favour quality (Table 2).Bin 301.2, influential in separating group 3 from group 1 and group 2 from group 1 (Table 1), gave a ranking consistent with this, of 3 = 2 > 1, and bin 329, influential in separating group 2 from group 1 (Table 1), gave a ranking consistent with this, of 2 > 1 = 3.

Effect of Nitrogen
The effect of adding N was examined by cultivar (Figure 2a (Malevo) and Figure 2b (Guerrero)).Adding N clearly alters the metabolic profiles of the exudates, but an OPLS-DA analysis of these two cultivars do not reveal common masses in the top ten discriminant bins.We hypothesised that the variation in field N could be creating additional variation and, therefore, examined the effect of revealing the blocks.Clearly there is both an effect of block and addition of N on the metabolic profiles of the exudates (Figure 3a,b).Block 3 clusters in the bottom left quartile.This field effect is creating much variation in the samples and, with the small number of replicates used, limits the resolution of the difference in metabolic profiles caused by N.
Effects of N could also be discerned with ANOVA of individual bins, where, in an overall analysis including all cultivars in the three quality groups, nine mass bins were found to show significant differences for N availability (Table 3).Some of these mass bins were also influential in the original overall analysis separating quality groups (bins 301.2 and 329 in Table 1).Values for a certain mass followed by a different letter were significantly different (LSD, p < 0.05).

Discussion
Bread-making quality is determined by genetic and environmental factors.The balance between gliadins and glutenins determines, at least partially, that the dough has properties suitable for baking.Quality is usually determined by conventional and time-costly methods.However, more recently, the genetic control of bread-making quality has been dissected into quantitative trait loci (QTL) in numerous studies, for example, micro-alveograph testing and sub-components [24], sedimentation in sodium dodecyl sulphate (SDS) [24][25][26], GPC [24,[26][27][28][29][30][31][32][33][34], and hardness [24], amongst others.As well as genetic factors, these traits are influenced by environmental factors and management practices, such as N and water availability, temperature and light intensity [35][36][37]; for example, increases in GPC can be achieved by N addition, but after incremental additions of N fertilizer, GPC reaches a maximum and then remains constant, without any increase in N uptake or remobilization by the crop, thus decreasing the efficiency of N fertilizer [1], hence, the importance of efficient tools to assist genetic selection and management practices.
In this work, we propose a novel approach to study bread-making quality by metabolic profiling.Electrospray ionisation mass spectrometry is highly suitable for analysing the wide range of non-volatile compounds present in phloem exudates: it is a soft ionisation method at atmospheric pressure; it is practical for molecules the size of metabolites; it provides good sensitivity; it is adaptable to a wide range of aqueous and organic solvents and, therefore, can be used directly with metabolite extracts prepared from plants; and it can analyse a wide range of different types of molecules, including highly polar molecules, such as peptides, oligonucleotides, and oligosaccharides, as well as small polar molecules, ionic metal complexes, and other soluble inorganic analytes.The analysis of ions by MS-TOF is simple and allows analyses of a virtually unlimited mass range with a resolution of 0.0001 and high sensitivity.
The approach taken in this demonstration of principle experiment was to use a non-targeted analysis of phloem exudates.Non-targeted approaches have been previously used in studies on tomatoes and Arabidopsis [3,12].PCA provides strong evidence that cultivars can be distinguished from each other, as well as between quality groups.For example, in Figure 1, Gavilán germplasm is clearly different to ACA 304.
The overall PCA identified PC 1 and PC 2, accounting for 56.5% of the total variation, which was relatively high compared to values obtained in other work, such as Rogers et al. [38]; presumably these cultivars differed in many aspects of their metabolism responsible for their different agronomic performance.Some mass bins showed large differences between quality groups and could be important for accounting for the differences responsible for belonging to different quality groups.For example, the mass bin 203, which was shown to be the first in the ranking of the top ten bins for separating group 3 from 2 and group 3 from 1 (Table 1), showed a ranking 1< 2 < 3 in the ANOVAs, meaning it may be associated with poor quality.Another mass bin showing large differences between quality groups in the ANOVAs was 370.6, that showed a ranking of 1 > 2 > 3 (Table 2).These mass bins could be markers for the rapid selection of cultivars for quality.
Other mass bins showing different rankings for quality groups in the ANOVAs (Table 2) give similar patterns to those observed in the top ten of Table Differences between N fertiliser levels identified by the PCA were to a certain extent masked by the differences between cultivars, but could still be extracted from the data.For example, mass bins 203 and 305 identified responses to N treatment in the cultivars.
When only three replicates are used, as in this study, the variation between replicates and the large amount of variation between cultivars clouds the separation.When the experimental blocks are overlaid on the data (Figure 3) it is clear that block three samples cluster tightly, whereas blocks one and two are more scattered.Thus, there is a large field effect in this experiment, meaning it may be beneficial in future experiments to increase the number of blocks in order to better take into account such heterogeneity [39].In spite of this, clear effects on the addition of N on the metabolic profiles of the exudates were observed in PCA.Future work will also de aimed at widening the number of cultivars under study and to relate their detailed quality characteristics with the masses analysed here; as mentioned in the introduction, these cultivars have been studied for quality characteristics and we propose to add our own quality data obtained from controlled field trials to these in future studies.
The period over which exudates were collected was short and the environmental conditions changed little over this time; for example, mean maximum, minimum, and mean temperatures were 22.9 ± 3.01, 6.73 ± 2.67, and 15 ± 1.7 • C, respectively, and rainfall was minimal (0.3 and 0.2 mm).Hence, we would not expect phloem exudate composition to be significantly affected.Moreover, Overy et al. [3] and Tetyuk et al. [40] collected exudates from plants differing in age (three days in the former and up to fourteen in the latter) and made no observations on this.The difference in the latter study resulted in collection from different phenological stages; we, in our study, were interested in analysing exudates from plants of the same stage, rather than the same age.
The results presented show that metabolic profiling may be used to extract biochemical markers that may be of potential use in selection, in the discrimination of cultivars of differing quality, and in elucidating the effect of N fertilisation.The challenge will be to identify which metabolites are those associated with that performance and to determine the definitive identification of the metabolites that goes beyond the putative metabolites included in the current study, in order to allow possible reasons to be postulated for those differences, and to generate the potential for biochemical marker selection for the important traits; future work will be directed towards this.Some of the mass bins identified as differing significantly between the cultivars may be involved in grain protein composition and quality characteristics.While these are possible candidates for explaining differences in agronomic performance, further analysis will be needed in order to establish these relationships and, as previously mentioned, the evaluation of a larger number of cultivars will be required.As data becomes available from field trials over several years, designed to explore the consistency of responses for the mass bins, this detailed analysis will be a focus for our resources.

Conclusions
In summary, we believe metabolic profiling has the potential to be developed into a breeding tool to refine plant breeding efficiency.In particular, it would appear to be useful to identify when field effects could be clouding the results of other screening methods.

Table 1 .
m/z bins that separate groups.The top 10 bins showing the strongest contribution in order of their p and q values.

Table 2 .
Influential mass bins showing significant differences between quality groups in ANOVA.

Table 3 .
Mass bins responding to fertiliser levels in ANOVAs applied to quality groups.