Zinc Excess Triggered Polyamines Accumulation in Lettuce Root Metabolome, As Compared to Osmotic Stress under High Salinity

Abiotic stresses such as salinity and metal contaminations are the major environmental stresses that adversely affect crop productivity worldwide. Crop responses and tolerance to abiotic stress are complex processes for which “-omic” approaches such as metabolomics is giving us a newest view of biological systems. The aim of the current research was to assess metabolic changes in lettuce (Lactuca sativa L.), by specifically probing the root metabolome of plants exposed to elevated isomolar concentrations of NaCl and ZnSO4. Most of the metabolites that were differentially accumulated in roots were identified for stress conditions, however the response was more intense in plants exposed to NaCl. Compounds identified in either NaCl or ZnSO4 conditions were: carbohydrates, phenolics, hormones, glucosinolates, and lipids. Our findings suggest that osmotic stress and the consequent redox imbalance play a major role in determining lettuce root metabolic response. In addition, it was identified that polyamines and polyamine conjugates were triggered as a specific response to ZnSO4. These findings help improve understanding of how plants cope with abiotic stresses. This information can be used to assist decision-making in breeding programs for improving crop tolerance to salinity and heavy metal contaminations.


INTRODUCTION
Farmers around the world only harvest an average of 50% of the potential yield they would obtain under optimal conditions (FAO, 2014). The gap between the actual yield and yield potential, known as "yield gap" is attributable to abiotic factors, especially salinity, drought, and heavy metals contamination (Moshelion and Altman, 2015). Excessive sodium chloride (NaCl) salinity in soil or water (Tester and Davenport, 2003) is a serious threat for agriculture. High concentrations of NaCl can inhibit plant growth and result in decline of productivity. This is due to (i) ionic and osmotic effects on nutritional balance and metabolic process, such as photosynthetic machinery (Parida and Das, 2005), nitrogen assimilation , protein synthesis (Giridara Kumar et al., 2003); (ii) induction of the over synthesis of reactive oxygen species (ROS) (Colla et al., 2010).
In addition, heavy metal toxicity is also considered a major abiotic stress affecting worldwide agricultural production (Kumar et al., 2015a,b). Unlike several heavy metals, such as cadmium (Cd), chromium (Cr), lead (Pb), and mercury (Hg), zinc (Zn) is considered an essential microelement for higher plants. Zinc is required at optimal concentration of 20 µg g −1 dry weight for the normal functioning of cell metabolism as well as plant growth and development (Broadley et al., 2007;Marschner, 2012). The normal range of Zn in the nutrient solutions for satisfying the crop requirement of Zn is between 0.05 and 0.50 mg/L (Jones, 2005). It is also involved in protein synthesis and many metabolic processes as one of the major cofactors of numerous enzymes (e.g., carbonic anhydrase, Cu/Zn superoxide dismutase, and matrix metalloproteinases; Cakmak, 2000). For most crop species, the critical toxicity concentration in leaf tissue ranges from 100 to >300 µg g −1 dw, the latter being more typical (Ruano et al., 1988;Marschner, 2012). However, anthropogenic activities such as mining, burning of fossil fuels, and agricultural practices have lead to Zn accumulation in soil (Nagajyoti et al., 2010). High Zn concentrations in soil and water can disturb physiological, biochemical, and metabolic processes leading to stunted plant growth by altering carbohydrate metabolism (Marschner, 2012) and photosynthesis, lowering the concentration of essential nutrients such as magnesium and iron (Sagardoy et al., 2010), causing oxidative damage to membranes, and interfering with DNA replication (Broadley et al., 2007;Vassilev et al., 2007).
Recent research into the development of "omics, " including genomics, ionomics, metabolomics, proteomics, and transcriptomics, has boosted "plant science" and helped clarify the functions of many key genes, proteins, and metabolite networks involved in plant responses under unfavorable soil and environmental conditions (Rodziewicz et al., 2014). The field of metabolomics has grown substantially over the past decade and has proven to be an important and efficient tool in plant responses to abiotic stresses (Rodziewicz et al., 2014) that allows the identification of potential biomarkers linked to improved stress tolerance. This can then lead to efficient genetic improvement programs (Weckwerth and Kahl, 2013). Knowledge of root molecular changes in response to stress is fundamental since roots are highly sensitive to several types of abiotic stresses (Jiang et al., 2007). Although plant response to salt stress has been only investigated in vegetables at leaf metabolic level , to our knowledge no published data is available related to root metabolomic responses to Zn exposure. As a large amount of information has been published concerning soil salinity, this study considered NaCl as a positive control, able to exert an iso-osmotic stress similar to ZnSO 4 .
The objective of the present study was to elucidate the metabolic changes in lettuce (Lactuca sativa L.), by probing the root metabolome of plants exposed to elevated concentrations of NaCl and ZnSO 4 .

Plant Material, Growth Conditions, and Treatments
Lettuce (L. sativa L. var. longifolia cv "Bionda degli Ortolani") seedlings purchased from a local nursery were transplanted at the three true-leaves stage, into plastic pots (diameter 16 cm) filled with commercial peat-based substrate having pH 6.2 in water, electrical conductivity 0.25 dS/m, porosity 92% v/v, CSC 98 meq/100 g, 140 mg/L N, 43 mg/L P, 150 mg/L K, 150 mg/L Ca, 50 mg/L Mg, 2.5 mg/L Fe, 1.3 mg/L Mn, 0.1 mg/L Zn. Lettuce plants were cultivated under controlled environment (average air temperature of 20 • C). Treatments were arranged in a randomized complete block design with 10 pots per treatment, each of them having four plants. The pot experiment consisted of three nutrient solution treatments: a non-salt control, 100 mM NaCl and 100 mM ZnSO 4 as described by Lucini and Bernardo (2015). The study was performed in terms of equimolar concentrations of the two different salts in order to evaluate the ion effects of Na and Zn on root metabolic profiling. The NaCl and ZnSO 4 nutrient solutions were applied 1 week after transplanting. Pots were watered as required by the crop (3-5 irrigations per week); all treatments received the same watering regime. At the end of the experiment (37 days after transplanting) leaves were collected, fresh tissue biomass measured and samples were then quenched in liquid nitrogen and grounded into a fine powder using mortar and pestle.

Mineral Analysis
The root tissues were analyzed for the following elements: K, Mg, Fe, Na, and Zn. Potassium and Mg were selected for their antagonism with Na for root uptake while Fe was selected for its antagonistic interaction with Zn. A portion (250 mg) of each sample was mineralized in 0.5 mL of hydrogen peroxide (H 2 O 2 ) and 1 mL of nitric acid (HNO 3 ) for 6 h, diluted in 5% nitric acid and analyzed by inductively coupled plasma atomic emission spectroscopy (ICP Iris; Thermo Optek, Milano, Italy; Karla, 1998). A multi-element source was used for calibration purposes and a certified reference material was analyzed prior to sample analyses. Each sample was read in triplicate.

Metabolomic Analysis
The untargeted screening of plant metabolites was carried out through UHPLC chromatographic coupled to a hybrid quadrupole-time-of-flight mass spectrometer (UHPLC/Q-TOF). An Agilent 1290 UHPLC liquid chromatograph was coupled to an Agilent G6550 mass spectrometer detector (Agilent Technologies Santa Clara, CA, USA) through a Dual Electrospray JetStream ionization system.
Reverse phase chromatography was carried out on an Agilent Zorbax Eclipse-plus column (75 × 2.1 mm i.d., 1.8 µm). The LC mobile phase consisted of water (A) and methanol (B), flowing at 220 µL min −1 and 35 • C. Formic acid 0.1% (v/v) and ammonium formate (5 mM) (both from Sigma Aldrich, St. Louis, MO, USA) were added to A and B. The gradient started with 5% B and increased to 90% B within 30 min, to be held for 5 min. The injection volume was 3 µL and mass spectrometer was operated in the positive scan mode to acquire spectra in the range of 100-1600 m/z. Regarding source conditions, sheath gas was nitrogen at 10 L min −1 and 350 • C, drying gas was nitrogen at 10 L min −1 and 280 • C, nebulizer pressure was 60 psig, nozzle voltage was 300 V and capillary voltage was 3.5 kV. Lock masses (m/z 121.0509 and 922.0098) were continuously infused in the ionization source by a quaternary pump through a dedicated electrospray.
Raw data were processed through a recursive analysis workflow, using Profinder B.05 (from Agilent Technologies) for feature alignment and filtering after initial deconvolution. Compounds identification was carried out using the whole isotopic pattern of deconvoluted features, hence taking into account accurate mass, isotope accurate spacing, and ratio.
Those compounds that were not present in 100% of the replicates within at least one treatment were discarded. Filtered features were next identified using the database exported from PlantCyc 9.5 (Plant Metabolic Network, http://www.plantcyc.org; released November 2014), using a mass accuracy tolerance threshold of 5 ppm. The peak volume of each compound identified with a score above 80/100 was extracted from the total ions current and exported for statistics and data interpretation. Additionally, the database was integrated using the software Pathway-to-PCDL (Agilent technologies), to create a compounds database (defined as "Personal Compounds Database Library" in the software tool) from the plant pathways available via KEGG (Kanehisa et al., 2016). In this way, annotated compounds could be identified, thus allowing pathway analysis from differential metabolites by using their KEGG ID.

Statistical Analysis
Interpretation of metabolomic results was carried out after recursive identification, using Mass Profiler Professional B.12.05 (Agilent technologies). The identified compounds were log 2 normalized and abundances baselined against the median of each compound in all samples. Statistics and chemometrics included analysis of variance (ANOVA at P = 0.001, Bonferroni multiple testing correction) and fold-change analysis (cut-off = 5), which were combined into volcano plots. Differential metabolites gained from the KEGG-derived database and selected through Volcano analysis were then used to perform pathway analysis using KEGG annotations to link metabolites and metabolic pathways.
Multivariate statistics were also carried out: Unsupervised cluster analysis (hierarchical cluster algorithm; similarity measure set to Euclidean) and Partial Least Squares Discriminant Analysis (PLS-DA; N-fold validation, N = 4) were performed. Concerning PLS-DA, the loadings (compounds) used to build the class prediction model were plotted according to their weight within the latent vectors, and those being most relevant in class prediction (i.e., those having a score of above +0.05 rather than below −0.05) were exported to support or complement information from Volcano analysis.

Plant Biomass and Mineral Composition
The lettuce biomass decreased in response to an increase in sodium chloride and zinc sulfate concentration in the nutrient solution with more detrimental effects resulting from the NaCl treatment. The percentage reduction in biomass production in comparison to the non-salt control was significantly lower in plants treated with ZnSO 4 (−25%) than in plants treated with NaCl (−53%, data not shown).
The macro and microelements concentration determined by ICP spectroscopy as a function of the nutrient solution treatments are displayed in Table 1. Sodium and zinc concentrations sharply increased in the corresponding treatments (NaCl and ZnSO 4 , respectively), and the effect was more marked for zinc, which increased more than 16 fold, compared to a <5-fold increase for sodium. Moreover, NaCl treatment decreased the lettuce root content of K while Fe concentration in root was reduced by ZnSO 4 treatment.

Metabolic Profile
Metabolomics in crop science are supposed to be prone to false positives, because databases are poorer than for model plants, and because they lack of tandem MS information. The combination of comprehensive databases like PlantCyc and KEGG in combination with proper data handling (i.e., the use of whole isotopic profile for identification, compounds alignments and filtering, and then recursive analysis) was effective in providing a wide metabolomic dataset and limiting false positives. All compounds were detected with low mass error (nominally below 5 ppm, and in the sub-ppm range in several cases). Starting from a total of more than 2400 compounds detected, recursive analysis and the following filters in Mass Profiler Professional dramatically reduced the compounds in the dataset to a number slightly above a thousand. The whole dataset, as gained after data processing, is provided as Supplementary Material.
Differential compounds were selected through Volcano analysis, where analysis of variance and fold-change were combined (using a fold-change cut off = 5 and a p-value of 0.001). Unsupervised cluster analysis (Figure 1), carried out from fold-change patterns (heat maps), and resulted in the definition of three well distinct clusters, one per treatment. Consistently, PLS-DA class prediction models gave accurate results (prediction  model overall accuracy = 100% after N-fold validation with N = 4), in which replications clustered each other within a treatment (Figure 2). Considering that PLS-DA could clearly differentiate the treatments, those compounds being differential (i.e., hose having the highest weight in the class prediction model) were exported from the compound loading view of PLS-DA and pooled with other differential metabolites. The majority of these differential metabolites had already been selected through Volcano analysis thus strengthening differences at metabolome level also through multivariate statistics. The compounds selected from Volcano analysis and PLS-DA scores are summarized in Tables 2, 3 (ZnSO 4 treatment against control, and NaCl treatment against control, respectively).
NaCl stress, despite being osmolar to Zn stress, triggered a more pronounced metabolic response, both in terms of fold-change values and as number of differential metabolites (Tables 2, 3). The volcano analysis for Zn-stressed vs. NaClstressed plants is given in detail as Supplementary Table. Ascorbate (2-carboxy-L-xylonolactone and 4-O-oxalyl-Lthreonate) and photosynthetic pigments (porphyrin ring and hydroxymethylbilane) degradation products were increased by NaCl stress, as well as homoglutathione. Analogously, phytoalexins (camalexin, medicarpin, zealexins, and TRIBOAglucoside) were mainly related to NaCl stress conditions (Table 3).

DISCUSSION
Root is the first organ feeling abiotic stresses in soil such as salinity and heavy metal contamination (Jiang et al., 2007), thus conditioning stress sensitivity and limiting plant biomass production (Steppuhn and Raney, 2005;Paradisone et al., 2015). Moreover, high Na concentrations in the soil or pore water may depress nutrient-ion activities and produce an extreme ratio of Na/K (Grattan and Grieve, 1999), whereas high zinc concentrations may induce deficiency of Mg and Fe because of the similar ion radius of Zn 2+ and Fe 2+ (Sagardoy et al., 2010), and Zn 2+ and Mg 2+ (Boardman and McGuire, 1990). Consistently, biomass reduction and imbalance of mineral nutrients were observed in our study. Indeed, the lowest K concentration was observed with the NaCl treatment, whereas the lowest Mg and Fe concentrations were recorded with the ZnSO 4 treatment ( Table 1).
As the pots were regularly watered with NaCl or ZnSO 4 enriched water solutions during the whole test period, we expect to have gained a steady state regarding the ions available in the substrate, after an initial transition period. Drainage of at least 50% of irrigation volume was achieved to avoid accumulation of salts in the substrate. At this steady state, the Zn 2+ and Na + concentrations in substrate solution were expected to be comparable and the plants should have been exposed to a similar osmotic stress.
Taking into account that metabolic changes induced by salinity have been widely reported in literature, our interpretation was mainly focused on biochemical changes related to zinc stress. In this view, salinity has been used as a positive control that exerted a similar osmotic stress to zinc sulfate treatment. As a general consideration for the differential metabolites, several compounds were identified in both ZnSO 4 and NaCl treated plants, with the latter inducing a more pronounced response to stress. The main classes of differential compounds were: carbohydrates, phenolics, hormones (gibberellins a cytokinin), glucosinolates and lipids. Salinity also induced changes in ascorbate and other compounds related to oxidative stress metabolism, together with phytoalexins and terpenes. However, zinc treatment also affected alkaloids and the polyamine profile in the roots ( Table 2).
The reduction of carbohydrates is in agreement with other studies, where a decrease in the carbohydrate route 2 | Output of Volcano analysis: metabolites resulting differential in lettuce root exposed to zinc stress against control, grouped according to functional classes.

Metabolites p (Corr) FC (abs) Regulation
Carbohydrates 2-amino-3,7-dideoxy-D-threo-hept-6-ulosonate 5.  of photosynthetic metabolism was associated to Zn stress (Chernyad'ev and Monakhova, 2006). In addition the increase in glucosinolates, traditionally linked to biotic stress, has been reported in heavy metals stressed plants (Poschenrieder et al., 2006) and even more in hyperaccumulator plant species (Asad et al., 2013). Phytoalexins, a class of compounds usually linked to biotic stress, have also been to abiotic stress (Jeandet et al., 2014). Therefore, their involvement in NaCl stress response, although less expected, can be justified. Osmotic stress is known to induce an oxidative secondary stress that can explain the increase in ascorbate breakdown products and in homoglutathione, which was observed in NaCl stressed roots. Sakuraba et al. (2014) have also reported that chlorophyll breakdown occurs during stress-induced yellowing, according to our findings in NaCl stressed roots. In previous work, on the use of a protein hydrolysate biostimulant on NaCl stressed lettuce,  achieved comparable results, reporting the involvement of glucosinolates, terpenes, homoglutathione, hormones, and carbohydrates in response to salinity stress.
As a large number of differential metabolites are common to both stresses (Figure 3, Venn analysis from differential metabolites under either Zn excess or salinity stress), it can be postulated that osmotic/oxidative stress is a process influencing lettuce roots following exposure to ZnSO 4 , in agreement with previous results achieved through proteomics using the same experimental design (Lucini and Bernardo, 2015). Furthermore, a larger number of differential metabolites could be evidenced under NaCl stress, suggesting this treatment altered lettuce metabolism in a way stronger than Zn.
The differential metabolites triggered by Zn stress were linked to KEGG metabolic pathways using their KEGG ID (where available via the Pathway-to-PCDL software), in order to drive and support the involvement of specific physiological responses during results interpretation. This analysis highlighted the involvement of the super pathway of polyamine biosynthesis (two matched entities: spermidine and 5-mehylthioadenosine; p < 0.001), spermidine hydroxycinnamic acid conjugate biosynthesis (two matched entities: spermidine and S-adenosyl-L-homocysteine; p < 0.001), chlorophyll a degradation (one matched entity: pheophorbide a; p < 0.001), and 2,3-flavanols biosynthesis (one matched entity: catechin; p < 0.001). Therefore, together with chlorophyll degradation and flavonols biosynthesis, the pathway analysis clearly confirmed the specific involvement of polyamine metabolism, and spermidine in particular, to ameliorate disorders from zinc stress. To provide further confidence regarding the involvement of spermidine in lettuce response to excess of Zn, two additional investigations were carried out. As first, a Venn diagram (Figure 3) was generated starting from the differential metabolites resulting from Volcano analyses (Zn vs. control and NaCl vs. control, respectively). Indeed, Venn diagram evidenced that spermidine and S-adenosyl-L-homocysteine were exclusive of Zn-related differential metabolites. In the second confirmatory approach, the mathematical model underlying PLS-DA analysis was exported, and the model formula coefficient for each compound was recorded. Spermidine showed a coefficient of +0.146 and −0.512 on first and second vector, respectively. These values, despite confirming that polyamine response under Zn excess was not the univocal metabolic response to stress, still are in agreement with the distribution of Zn-related samples in the PLS-DA hyperspace as provided in Figure 2. Our findings are in agreement with previous work on the role of polyamines in alleviating the adverse effects of heavy metals stress (Luo et al.,FIGURE 3 | Venn analysis from metabolites resulting differential after Volcano analysis under either ZnSO 4 or NaCl stressed plants. 2009; Aldesuquy et al., 2014). In addition to playing an important role in cell elongation and division in several plant species, polyamines are involved in the response to stress and their levels in plant increase under a number of environmental stress conditions such as drought, nutrient deficiency and phytotoxic metals (Aldesuquy et al., 2014). In several plants, polyamines can exist as free amines rather than being conjugated with hydroxycinnamic acids (hydroxycinnamic acid amides; Luo et al., 2009). Luo et al. (2009) reported that conjugated polyamines are accumulated in seeds and sometimes in roots, as a response to both biotic and abiotic stresses. A specific polyamine glucosyl transferase has been isolated (Luo et al., 2009), however its role in stress response is not fully elucidated and it has been postulated on the basis of correlative data. Recently, lettuce root pre-treatment with spermidine has been reported to enhance CO 2 assimilation efficiency and to affect stomatal regulation (Sun et al., 2016). In the current experiment, polyamines have been specifically induced under Zn conditions, presumably to cope with the stress this metal induces to roots.

CONCLUSIONS
Despite emerging as a diffuse soil contamination, not only linked to well-known point sources such as mining or steel industry, the stress from Zn excess to plant has not been deeply investigated yet. On the contrary, information on the effects of salinity to plants at molecular level is pretty available. In our work, lettuce root response to NaCl and zinc-induced stress was investigated at metabolome level under equimolar concentrations, to use elevate salinity as a "positive control." The majority of differential roots metabolites were shared between the two stresses, although fold-changes were higher under NaCl stress condition. Most of the identified compounds could be ascribed to carbohydrates, phenolics, hormones, glucosinolates, and lipids. Our findings suggest that osmotic stress and redox imbalance play a major role in determining lettuce root response to NaCl and ZnSO 4 . However, accumulation of polyamines and polyamine conjugates has been specifically triggered by zinc.
The results provide insight into plant function under different environmental stimuli and are important in revealing the physiological/molecular mechanisms determining plant response to abiotic stresses, and to zinc exposure in particular. These findings at molecular level are useful to understand the biochemical processes that plant uses to cope with zinc stress. Our results might be also useful to assist the selection of tolerant cultivars and, possibly, be translated to other horticultural crops.

AUTHOR CONTRIBUTIONS
LL designed the study, in cooperation with YR and MT. LL also coordinated the experimental work, developed the mass spectrometric method, performed statistics, and helped to draft the manuscript. LB carried out the plant experiments and metabolomic work, contributed to interpretation of data and drafted the manuscript. YR, GC MT, and DK participated in the interpretation of data, helped to draft the manuscript and to critically revise it for intellectual content. All authors read and approved the final manuscript, whereas LL was responsible of final approval of the version to be published.