LC-HRMS Based Non-Targeted Metabolomic Profiling of Wheat ( Triticum aestivum L . ) under Post-Anthesis Drought Stress

Drought stress at the reproductive stage causes severe damage to productivity of wheat. However, little is known about the metabolites associated with drought tolerance. The objectives of this study were to elucidate changes in metabolite levels in wheat under drought, and to identify potential metabolites associated with drought stress through untargeted metabolomic profiling using a liquid chromatography-high resolution mass spectrometry (LC-HRMS)-based technique called Isotopic Ratio Outlier Analysis. Metabolomic analysis was performed on flag leaves of drought-stressed and control (well-watered) plants after 18 days of post-anthesis drought stress at three-hour intervals over a 24-hour period. Out of 723 peaks detected in leaves, 221 were identified as known metabolites. Sixty known metabolites were identified as important metabolites by 3 different methods, PLS-DA, RF and SAM. The most pronounced accumulation due to drought stress was demonstrated by tryptophan, proline, pipecolate and linamarin, whereas the most pronounced decrease was demonstrated by serine, trehalose, N-acetyl-glutamic acid, DIBOA-glucoside etc. Three different patterns of metabolite accumulation were observed over 24-hour period. The increased accumulated metabolites remained higher during all 8 time points in drought stressed leaves. On the contrary, metabolites that showed decreased level remained significantly lower during all or the most time points. However, the levels of some decreased metabolites were lower during the day, but higher during night in drought stressed leaves. Both univariate and multivariate analyses predicted that *Co-first authors. How to cite this paper: Rahman, M.A., Akond, M., Babar, M.A., Beecher, C., Erickson, J., Thomson, K., De Jong, F.A, and Mason, R.E. (2017) LC-HRMS Based Non-Targeted Metabolomic Profiling of Wheat (Triticum aestivum L.) under Post-Anthesis Drought Stress. American Journal of Plant Sciences, 8, 3024-3061. https://doi.org/10.4236/ajps.2017.812205 Received: September 25, 2017 Accepted: November 12, 2017 Published: November 15, 2017 Copyright © 2017 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access


Introduction
Wheat is one of the most widely produced and consumed cereal grains worldwide and global wheat demand is expected to rise in excess of 880 million metric tons by 2050 [1].Production of major crops, including wheat, is declining in some regions of the world due to climate change, especially progressive drought stress, induced by either declining rainfall or higher water stress allied with a hot climate [2] [3] [4].Water stress accompanied by other abiotic stresses can decrease potential crop yield by >50% [5].Water stress during the flowering stage and onwards in wheat typically interrupts photosynthesis and increases translocation of carbohydrates from source to sink, reduces grain number, weight and yield, accelerates phasic development, reduces grain filling duration and enhances senescence [6] [7] [8] [9] [10].Genetic improvement is the most effective and sustainable method to reduce the detrimental effect of water stress.However, the genetic improvement of traits that confer stress tolerance can only be achieved when we have a better understanding of the biochemical mechanism controlling those traits.The analyses of biologically important molecules are necessary to understand which molecules are influencing stress tolerance mechanisms in wheat plants.Under water-stress, plants produce an array of biomolecules including different metabolites [11].Plants can modify their physiology to acclimate to abiotic conditions through metabolic changes [12] [13].Changes in the metabolic level in an organism are likely correlated with the phenotype, as metabolites are the end-products of biological systems [14].Other than the changes in the level of metabolites, little is known about their complex metabolic regulation under water stress or potential role in developing climate resilient wheat.
Crop yield has shown significant correlation with metabolites under drought [15].Metabolic composition is highly variable in a wide variety of plants including major crops like wheat [16], rice [17], maize [15], and cassava [14].These metabolic changes in plants have also occurred due to progressive water deficit conditions accompanied by lower leaf water potentials [18].Water stress induces metabolic changes and accumulation of metabolites including proline, ascorbic acid, glutathione, phenolics, and detoxifying enzymes in cereals [19] [20].Bowne et al. (2012) studied metabolic responses in wheat grown under contrasting water levels (drought and irrigated conditions) and reported that amino acids, most notably proline, tryptophan, and the branched chain amino acids, including leu-American Journal of Plant Sciences cine, isoleucine, and valine exhibited higher levels of expression in the leaves of drought grown plants [21].Additionally, different metabolites including phenolic compounds, amino acids, fatty acids, organic acids, carbohydrates, and sterol-based compounds can be altered depending on the phenology and physiological status of the plant under stress [22].Stress tolerance is a polygenic trait that needs to be incorporated in wheat to develop climate resilient germplasm.
Quantitative trait loci can be identified and tagged on genomic regions in a plant using suitable biomarkers through application of molecular breeding techniques.
If a set of water stress related metabolic biomarkers can be identified for wheat, they could be used as targeted, fast, and low-cost diagnostic tools to select superior performing lines in breeding programs.
Metabolic profiling allows for comprehensive analyses of a range of metabolites that have great value in both phenotyping and plant diagnostics [23].Recent progress in mass spectrometry with advanced data processing technology allows simultaneous measurement of hundreds of chemically different metabolites and investigates more thoroughly the regulation of metabolic networks to study their influence on complex traits.Such metabolic profiling can meaningfully contribute to the study of stress biology in wheat and could offer a set of drought related metabolite biomarkers.To date, there are insufficient biomarkers available to adequately screen cereal crops for drought tolerance or susceptibility.The objectives of this study were: 1) to elucidate the differential metabolite accumulation in wheat leaves under post-anthesis drought stress conditions over 24 hours time period, and the involvement of those metabolites in different pathways in relation to drought tolerance and 2) to identify potential metabolite biomarkers associated with drought stress in wheat.To achieve this objective, we employed a non-targeted LC-HRMS Isotopic Ratio Outlier Analysis (IROA) Global Metabolomics method [24] for identifying metabolites from the leaf tissue of the drought stressed and control wheat plants at 8 time points of 24-hour period.

Plant Material and Growth Conditions
The soft winter wheat variety 'SS8641' (GA-881130/2*GA-881582) is widely grown in the southeastern USA, developed by the University of Georgia Wheat Breeding Program, and was subjected to post-anthesis drought stress and LC-HRMS based untargeted metabolic profiling in this study.The SS8641 is a mid-maturing, high yielding wheat variety and has high test weight with good straw strength.It showed resistance to Hessian fly biotypes B and E, and possesses powdery mildew genes Pm1 and rust resistance genes Lr37/Yr17/Sr38.
Since the wheat flag leaf is a vital source of energy assimilates during grain filling [25] [26], we initially screened a collection of soft wheat germplasm (developed by different US public wheat breeding programs) for thylakoid membrane (indirect assessment of Photosystem II damage) and chlorophyll damage under water stress at flag leaf.The SS8641 wheat genotype showed minimal thylakoid drought treatment (water-deficit, 25% FC).To determine FC, the weights of pot plus dried soil weight were recorded as average 2.24 kg, then pots were irrigated to 100% FC and average weight increased to 3.44 kg (moisture content was assessed as 1.20 kg/pot).Pots in the control treatment were maintained at 100% FC, while drought treatments were maintained up to 25% of the FC.After 18 days of the drought treatment, plants were characterized for photosystem and chlorophyll damages, and flag leaves were collected for metabolic profiling.

Physiological Characterization
Chlorophyll fluorescence (the ratio of variable, Fv, to maximum fluorescence, Fm) and SPAD chlorophyll content were used as indirect methods to assess thylakoid membrane [27] [28] and chlorophyll damages due to stress following the der drought stress and control conditions respectively, and X is assumed to be 100% leakage corresponding to 100% membrane damage.

Measurements of Growth and Yield
Plants were harvested after physiological maturity and oven-dried at 60˚C for five days.After drying, total biomass was measured per plant.The spikes were hand threshed and the chaff was cleaned off.Grain weight was measured and the harvest index was calculated by dividing grain weight/total biomass.Grain yield per spike was calculated by dividing total grain weight by spike number.Grain number per spike was calculated by dividing total grain number by the number of spikes for that plant.Two hundred grains were randomly selected and weighed, and converted to 1000-grain weight.

Leaf Tissue Collection and Sample Preparation for Metabolomics
Flag leaves were collected in 3-hour increments (06:00 AM; 09:00 AM; 12:00 noon; 15:00 PM; 18:00 PM, 21:00 PM, 24:00 PM and 03:00 AM) over a 24 hour time period.Collected flag leaf from each pot was considered a biological replication.A triplicate samples were collected at each sampling time point as our experiment was maintained in a well environmentally controlled greenhouse and experimental treatments showed high accuracy within replication.Sampled leaf tissues were frozen in liquid nitrogen immediately after collection and then stored at −80˚C until processing.Leaf tissue samples were lyophilized and ground using a tissuelyser (24 samples per treatment for a total of 48 samples).
For metabolomic analysis, 5 mg of experimental material was weighed and added to 5 mg of freshly ground wheat internal standard (IS).The IS was an isotopically labeled wheat leaf that had been grown in an atmosphere of 13 C labeled carbon dioxide resulting in a uniform and universal labeling of approximately 97% (IROA Technologies) [24].Next, 500 µl of methanol/10mM aqueous ammonium acetate (50:50) was added to the dried powder and vortexed for 1 min at room temperature.The resulting mixture was further sonicated for 20 min, and centrifuged at 17,000 G for 10 minutes at 4˚C.A 350 µl aliquot of the supernatant was transferred to a clean dry eppendorf tube, and dried under a gentle nitrogen gas stream (at 30˚C).The dried sample was reconstituted in 100 µl of 0.1% aqueous format.

Untargeted Metabolomics
The reconstituted samples were analyzed for untargeted metabolites on a Liquid Chromatography High Resolution Mass Spectrometery (LC-HRMS) platform.
Untargeted and the repeated analysis of a large pooled plasma sample.Injection reproducibility was typically less than 10% even without a ratio to an internal standard.

Metabolic Data Processing
One of the greatest challenges of most metabolic profiling experiments is the ability to differentiate peaks of biological origin from artifact peaks, and to accurately identify and quantitate the peaks of interest.Since we used an IROA-labeled plant material as our internal standard this study followed the IROA "Phenotypic" global labeling and bioinformatics protocols in which the Internal Standard (IS) is labeled at 95% 13 C. Therefore, all biological compounds are paired natural abundance (NA) and IS, and each pair carries distinct molecular signatures.Molecules can be distinguished from each sample set, as they have differing masses [31] [32].
For IROA, control and drought treated samples were analyzed as a single composite sample by LC-MS.Algorithms pair identified biological peaks, and unlabeled NA artifacts were identified and discarded.All biological compounds had two paired peaks; the peak from the 12 C-media is mirrored by a second peak from the 13 C-media.The distance between the monoisotopic peaks readily identified the number of carbons in the compound.The corresponding M +1 and M −1 peaks (and M +2 and M −2 etc. peaks) which are a mass difference of 1.00335 amu (mass difference between a 12 C and 13 C isotope), gave the IROA peaks a characteristic U-shape "smile" pattern.Accurate mass together with the knowledge of the number of carbons in a molecule greatly facilitated metabolite identification.American Journal of Plant Sciences The ClusterFinder program was used to identify, align and quantitate the IROA peaks in the .mzXMLfiles.The files were scanned for IROA peaks down to an intensity level of 1 million in both positive and negative modes with an assumed maximum of 10 ppm mass error.Settings for the experimental samples in the nontargeted analysis assumed a natural abundance (1.1%) isotopic balance while the IS was assumed to have a 97% isotopic balance.The resulting list of IROA peaks was manually curated to define 394 compounds that were seen in both the IS and in the experimental samples.Once the physical attributes for the 494 compounds were identified, a targeted analysis for all of these compounds was imposed on every sample.This resulted in a non-sparse dataset, i.e. there was a value for every compound for every sample and the dataset was exported for analysis.The metabolites were annotated by searching against an in-house metabolite database, Mass Spectrometry Metabolite Library of Standards (MSMLS) (http://iroa.com/page/Mass%20Spectrometry%20Metabolite%20Library%20of%20Standards).

Metabolic Data Analysis
Data tables with metabolite peaks (mz/rt) at 8 time points under both drought and control conditions were formatted as comma separated values (.csv) files and uploaded to the MetaboAnalyst 3.0 server (http://www.metaboanalyst.ca)[33].To shrink any possible variance and to improve the performance for

Physiological Traits
We measured physiological traits to assess the impact of the drought stress 18 days after stress initiation at anthesis.The chlorophyll content (SPAD value), maximum quantum efficiency of PSII (Fv/Fm) and electrolyte leakage (%) were determined on flag leaves for both control (well-watered) and drought conditions (Table 1).Drought stress caused a significant reduction (~11%; p < 0.05) in chlorophyll content of wheat flag leaves compared to control.The photosynthetic apparatus (PSII) was damaged to some extent which was apparent from the reduced Fv/Fm value under drought condition (Table 1), although the reduction was not statistically significant.Electrolyte leakage, in contrast, increased substantially (~18%) under drought condition compared to control resulting in ~30% damage to plasma membranes of the drought stressed plants.

Morphological Traits
Total shoot biomass at maturity, grain weight/plant, grain weight/spike, 1000-grain weight, grains/spike, and harvest index showed significant reduction under drought stress compared to control condition as expected (Table 1).
Table 1.Physiological and morphological traits of SS8641 wheat line under post-anthesis drought stress and control conditions.

Profiling of Leaf Metabolites
Metabolite profiling by LC-HRMS detected a total of 723 peaks from wheat flag leaves.Among the detected peaks, 221 were identified as known metabolites (266 including duplications) and the remaining peaks were unknown metabolites.The identified compounds included amino acids, sugars, organic acids, organic compounds, polyamines, fatty acids, nucleosides/nucleobases and other compounds.
The supervised clustering method, Partial Least Squares-Discriminant Analysis (PLS-DA) was performed for two conditions at 8 different time points.Five PLS-components (PCs) explained 66.5% of the total variation, with the first and second PCs contributing 45.1% and 8.5%, respectively (Figure 1).Drought and control samples were clearly separated across the eight different sampling points in the scores plot between PC1 and PC2 (Figure 1) with minimum overlapping within the group.This separation of drought and control samples clearly indicates the altered state of metabolite levels in the wheat leaves under drought condition.
Significantly different metabolites were analyzed by hierarchical clustering with heat map in order to visualize the effect of drought stress over the control at eight different time points (Figure 2).The heatmap formed two major clusters  A comparison of statistical models was carried out to identify the important metabolites associated with drought condition using: SAM, PLS-DA and RF (Table 2).The SAM plot identified 61 significantly different compounds with Over the course of 24 hours 89 metabolites exhibited significant fluctuations (ANOVA, P ≤ 0.05, Fisher's LSD) across the 8 sampled time-points between the two treatments (Supplemental Table S2).The most significant metabolites including tryptophan, proline, pipecolate, linamarin, indoline and phenylalanine remained significantly higher under drought stress compared to control across all 8 time points with the highest peak at 18 h (3 pm) (Figure 3).On the contrary, there were several sugars, fatty acid, organic acids, polyamines and amino acids which were significantly higher in control compared to drought stress across most or all of time points (Figure 3).Fluctuations were also seen based on day and night cycles.For example, levels of serine, aspartate, raffinose, myo-inositol, spermidine, gluconic acid, DIBOA-glucoside and tribenuron methyl increased beginning at 06 h and remained high during the day time and then decreased at night (21 h) in control leaves, whereas their levels were low during day time but increased at or after 24 h during night in drought stressed leaves.
To better elucidate the biological functions of identified metabolites, a pathway analysis was performed using Oryzasativa and Arabidopsis thaliana as the pathway libraries.As expected, these metabolites were involved in number of different pathways (Table 3).As the rice and Arabidopsis pathway libraries do not have all the compounds, KEGG website was also used to find their metabolic pathways.

Identification of Metabolite Biomarkers
Applying a ROC (Receiver Operating Characteristic)-curve based approach of biomarker analysis [33], MetaboAnalyst 3.0 identified eight metabolites with high AUC values (Area under ROC curve, >0.84) as potential biomarkers for drought tolerance based on classical univariate ROC curve analysis module.

Discussion
Drought This declining photochemical activity and membrane integrity of wheat plants (SS8641) under drought was also reflected in lower shoot dry biomass, grains spike −1 , grain weight plant −1 and harvest index compared to the control.
In addition to physiological and morphological traits, levels of metabolite accumulation during stress at a particular growth stage can provide a more specific and accurate indication of stress tolerance.In the present study, the non-targeted metabolite profiling revealed significant changes in 152 metabolites in wheat flag leaves under drought stress condition compared to control.Of these differentially accumulated metabolites, 60 metabolites were identified as the most important metabolites by three different methods: PLS-DA, RF and SAM, which included an array of different amino acids, sugars, organic acids or compounds, polyamines, fatty acids and their derivatives.Amino acids including tryptophan, proline, phenylalanine, tyrosine and isoleucine accumulated up to 7.1 fold in drought treated flag leaves compared to control, although the levels of some amino acids (e.g.serine, glutathione, glutamine, GABA, threonine) decreased due to drought stress.The accumulation of amino acids under drought stress has been reported previously in several studies [18] [21] [36].An increased level of amino acids is considered to improve stress tolerance in plants through different mechanisms, such as osmotic adjustment, detoxification of reactive oxygen species, and intracellular pH regulation (see review by Krasensky and Jonak 2012 [37]).
Tryptophan acts as an osmolyte in ion transport regulation and in modulating stomatal opening, and maintaining water balance between air and plant [38].
Increased accumulation of tryptophan in wheat under drought is also supported by previous studies [21] [39] [40] [41].Higher levels of tyrosine and phenylalanine under drought stress are also consistent with a study by Bowne et al. ( 2012) [21].The aromatic amino acids, tryptophan, tyrosine and phenylalanine are the precursors of different secondary metabolites including indoleacetate, lipid precursor and lignin in shikimate pathway, which plays a vital role in stress tolerance [42] [43] [44].
Proline level was increased in flag leaves under drought conditions and the accumulation of proline has been correlated with stress tolerance in a wide range of plants by different scientists [45] [46].Proline serves as a compatible solute and helps plants to avoid oxidative stress by preserving osmotic balance or cell turgor, stabilizing membranes, stopping electrolyte leakage, and by keeping reduced concentrations of reactive oxygen species (ROS) [47] [48] [49].Furthermore, the increased level of branched chain amino acids (e.g.isoleucine) under drought stress was reported in tolerant cultivars of wheat [21] [39] [40], barley [50], and Arabidopsis [51], and similar result was found in this study.Taylor et al. (2004) reported that branched chain amino acids are a source of alternate energy when Arabidopsis plants were under sugar starvation, which is a common phenomenon under stress [52].
Sugars and their derivatives were significantly reduced under drought stress.
Lower accumulation of glucose and sucrose is potentially due to reduced photosynthetic capacity, as demonstrated by lower chlorophyll fluorescence and SPAD chlorophyll content, of the drought stressed leaves at the grain filling stage [34] [35].In addition, enhanced maturity under stress condition may have resulted in increased partitioning of carbohydrates from source (leaf) to sink (grain).Depleted starch content in drought-stressed barley flag leaves compared to control has been reported under stress condition [53].Trehalose is reported to be accumulated at a moderate level under abiotic stress and serves as an osmolyte to stabilize proteins and membranes [54] [55].The trehalose of drought stressed leaves in our study may have been utilized as an osmolyte to provide the drought tolerance observed in the genotype SS8641.Myo-inositol acts as precursor of many metabolites especially raffinose that are involved in stress tolerance by providing membrane stability and antioxidative functions [56].Decreased levels of raffinose of the drought stressed leaves can be accounted to the drought tolerance of the wheat plants (SS8641).Obata et al. ( 2015) also reported lower levels of raffinose and myo-inositol in a tolerant maize genotype and suggested that that these could be due to the use of raffinose family oligosacchardises (RFOs) by the tolerant genotypes as their carbon sources [15].
The lower accumulation of sugars was further accompanied by reduced levels of organic acids or compounds including pyruvate, phosphoenolpruvate, α-ketoglutaric acid which are involved in Krebs cycle in leaves of stressed plants.
Drought stress in our study didn't trigger higher accumulation of organic acids or compounds except pipecolate and linamarin.A decreased level of organic acids or compounds was reported in wheat [21] and soybean [36] under drought stress in tolerant genotypes.In contrast, the consistent higher accumulation of pipecolate and linamarin at all time points upon drought stress indicates their potential role in drought tolerance of wheat.
Pipecolate is reported to have a defensive action on proteins and nucleic acid structures by maintaining a stable osmotic status in plants under variable soil water and salt stress [57].Bowne et al. ( 2012) also found significantly changed levels of pipecolate among wheat leaves grown under drought stress compared to control [21].Drought stress led to a significantly higher accumulation of linamarin level in the leaves.Linamarin is a cyanogenic glycoside, which plays an important role in nitrogen transport [58] and can act as an immediate chemical defense against insect and pathogen infestation by discharging toxic cyanide [59].Higher level of linamarin was reported in cassava leaf under drought stress [14].Our study first demonstrated a consistent higher accumulation of linamarin under drought stress at different time points in wheat flag leaves which indicates that linamarin has a strong relationship with drought tolerance in wheat at grain filling stage.The accumulated linamarin under drought stress could further be hydrolysed to release hydrogen cyanide, which can act as precursor of amino acids.This is potentially attributed to higher accumulation of amino acids under drought conditions in our study.
The level of different metabolites changed in wheat flag leaves at 8 different time points across a diel cycle and they changed differently under drought and control conditions (Figure 3, Supplemental Table S2).Metabolites showed mostly 3 different diel patterns of accumulation in flag leaves at 8 different time points.Most of the amino acids (e.g.tryptophan, proline, indoline, phenylalanine) and some organic acids (e.g.pipecolate, linamarin) were accumulated and remained higher in drought stressed leaves over the 24 hrs of day-night period with the maximum accumulation at 3 pm compared to controls.But the levels of sugars, fatty acid, organic acids, amines and some amino acids were lower in drought plants over the most or all of time points compared to controls.Our study also showed that the levels of serine, aspartate, raffinose, myo-inositol, spermidine, gluconic acid, DIBOA-glucoside and tribenuron methyl remained higher during day time but decreased at night under control condition but the diel pattern of metabolite accumulation became completely opposite under drought condition.This pattern of metabolite accumulation was also reported by Benard et al. ( 2015) [60] where sucrose and serine were accumulated during day time and decreased at night under control condition.
The pathway analysis by using MetaboAnalyst was able to link with 33 metabolic and biosynthesis pathways of the KEGG pathway database and Rice Annotation Project database (Table 3).The aromatic amino acids such as phenylalanine, tyrosine and tryptophan are produced via the shikimate pathway having fundamental roles in plant reproduction, development, pest-defense, and environmental stresses [43] [61].Serine, glycine, and threonine all can be converted to the 3-carbon a-ketoacid pyruvate or pyruvate which was downregulated in this study under drought stress.Pyruvate is the end product of glycolysis in the cytosol, and could be transported by mitochondrial pyruvate carriers (MPCs) American Journal of Plant Sciences into mitochondria for consequent cellular substance and energy metabolism [62].Arginine and proline metabolism is one of the unique pathways for the biosynthesis of the amino acids.The regulatory mechanisms controlling proline metabolism, intercellular and intracellular transport, and connections of proline to other metabolic pathways are all imperative to the in vivo functions of proline metabolism [63].Proline accumulated mostly in the cytosol due to water stress reduces cytoplasmic acidosis and conserves NADP+/NADPH ratios well-suited with metabolism [64].Upon relief of stress, the proline metabolic pathway reduces mitochondrial oxidative phosphorylation and produces ATP for repairing of stress-induced damages [49] [64].The increased or decreased levels of metabolites in our study integrated with their linked major metabolic pathways and biosynthetic processes presents a complicated process of drought stress responses in wheat.
We also conducted the analysis to find potential metabolite biomarkers associated with drought tolerance in wheat for the screening of drought tolerant genetic resources of wheat.Obata et al. (2015) opined that promising metabolic traits were stronger in explaining variability in maize grain yield than classical agronomic yield components [15].Six metabolites namely N-acetyl-glutamic acid, tryptophan, proline, linamarin, pipecolate, malonate and DIBOA-glucoside were identified as potential biomarkers by classical univariate and multivariate analysis.The finding of proline and tryptophan accumulation in drought tolerant wheat genotype was not something new.They were already reported as associated with their abiotic stress tolerance by several studies in different crops [37].But the identification of N-acetyl-glutamic acid, linamarin, pipecolate and DIBOA-glucoside in wheat genotype to be associated with drought tolerance was not reported earlier and could be promising metabolite biomarkers for screening the tolerant genotypes.However, further investigation needs to be conducted in a panel of wheat germplasm with different genetics backgrounds to draw a comprehensive conclusion.

Conclusion
This study is the first study in wheat which used a non-targeted LC-HRMS Isotopic Ratio Outlier Analysis (IROA) Global Metabolomics method to identify accumulated metabolites in flag leaves under drought and control conditions at 8 time points during 24 hours day-night period.The study reports increased or decreased level of numerous metabolites simultaneously and consistently at different time points.Some of metabolites received considerable research attention under drought stress, but many of these did not.Metabolites like tryptophan, proline, pipecolate, linamarin, N-acetyl-glutamic acid, DIBOA-glucoside accrued to a greater level after drought stress, which was likely the indication of acclimation in responses to the drought stress.Elevated tryptophan and proline are common in other crops and are not surprising, but increased levels of pipcolate, linamarin, and decreased levels of N-acetyl-glutamic acid and DIBOA-glucoside American Journal of Plant Sciences in response to drought stress have not been widely reported.These metabolites may serve as metabolite biomarkers for screening drought stress tolerance in wheat germplasms to develop climate resilient wheat.The finding of these metabolites could potentially be the effect of genetic variation, as we used soft wheat which is different than other wheat genotypes used in the study before, including in Bowne et al. ( 2012) [21].This difference indicates that global wheat community needs to use diverse set of germplasms for studying metabolite accumulation under drought stress.That will improve our understanding on metabolic control of complex stress conditions including drought in wheat.In addition, we employed LC-HRMS analysis with IROA method which was not also reported in wheat before.This technique is sensitive and potentially more robust than Nuclear magnetic resonance (NMR) and gas chromatography-mass spectrometry (GC-MS) technology as it can identify metabolites in low concentration and has better capacity to control FDR rate.Although some of these newly identified metabolites are promising as biomarkers for improving drought stress tolerance in wheat, their correlation to drought tolerance in wheat and other cereal crops requires further investigation.

M
. A. Rahman et al.DOI: 10.4236/ajps.2017.8122053027 American Journal of Plant Sciences membrane and chlorophyll damage under long-term water stress.This finding led to further investigation of metabolites accumulated during grain filling under drought stress in SS8641.Initially, seeds of SS8641 were exposed to vernalization at 4˚C for 6 weeks to induce flowering.After vernalization, three germinated seeds were planted in each pot and kept in an environmentally controlled greenhouse located at the campus of the University of Florida in Gainesville.The greenhouse was maintained at day/night temperature of 20/15˚C ± 0.5˚C (16 h daylight and 8 h night time) with a relative humidity of 50% ± 2%.Plants were manually irrigated at three day intervals and maintained at 100% field capacity (FC) until exposed to drought stress.Plants were also fertilized with two splits of 12 g of Scott's Osmocote (20-4-8 NPK fertilizer) during the experimental period: one at 7 days and another split at 45 days after transplanting.Once plants reached the flowering stage, pots were randomly assigned to two different watering regimes with three replications as follows: control treatment (well-watered, 100% FC) and methods described by Talukder et al.(2014)  [29].As photosystem II is housed at the thylakoid membrane, damage to thylakoid membrane is considered as an indicator of damage to photosystem II under stress condition.The lower ratio of Fv/Fm indicates more damage to Photosystem II due to stress.Measurements of Fv/Fm were taken on the abaxial surface of intact flag leaves at one third of the length from the ligule after 30 min of dark adaptation on three flag leaves of each plant (biological replication).Fluorescence was measured using a pulse modular fluorometer (Model OS5-FL, Opti-Sciences, Hudson, NH, USA) in both the control and drought stressed leaves at 18 days after flowering.A self-calibrating SPAD chlorophyll meter (Model 502, Spectrum Technologies, Plainfield, IL) was used to measure chlorophyll content on the same flag leaves and leaf blade areas where chlorophyll fluorescence measurements were taken.In each treatment (control and drought), chlorophyll content was measured from three flag leaves (3 different plants) and the average of three measurements was used to represent chlorophyll content.The lower value of SPAD chlorophyll content indicates more damage to chlorophyll due to stress.Plasma membrane damage (PMD) was measured using the method described byRistic and Cass (1993) [30].Leaf disks (diameter = 5 mm) were collected from six individual flag leaves within each biological replication after the drought treatment and placed in de-ionized water (24 ml) in sealed vials.The vials were shaken overnight on a shaker at 5˚C.Electrical conductivity (µS/cm) of the aqueous solution was measured with a Metter Toledo (Seven-Multi S70) conductivity meter.The tissue samples were then autoclaved for 20 min.The conductivity of the solution was again measured after shaking the samples overnight at 5˚C.The percent electrolyte leakage was calculated based on the conductivity before and after autoclaving.The average value of six flag leaves within each biological replication was used to estimate %PMD according to the formula, [(%L D − %L C )/(X − %L C )] × 100, where L D and L C are electrolyte leakages un-

Figure 1 .
Figure 1.Partial least square discriminant analysis (PLS-DA) and 2D Scores plot for the wheat SS8641 leaves over the 24 h time period (8 time points) under control (irrigated) and drought conditions.Samples under control and drought condition at different time points did not overlap with each other (although overlapped within group), indicating an altered state of metabolite levels in the wheat leaves.

Figure 2 .
Figure 2. Heatmap showing levels of top metabolites based on partial least square discriminant analysis (P LS-DA) VIP scores at 8 different time points under control and drought conditions.The heatmap was generated using "Pearson" for distance measure and "Ward" for clustering algorithm.The hierarchical clustering was performed only on the metabolites and the samples were arranged based on their time points in two groups.

Figure 3 .
Figure 3. Significantly different levels of selected metabolites (ANOVA, P ≤ 0.05, Fisher's LSD) in the wheat flag leaves at 8 different time points (03 h, 06 h, 09 h, 12 h, 15 h, 18 h, 21 h, 24 h) during 24 hours of day-night period under drought and control conditions.Dotted line (----) and solid line indicate the metabolite levels under control and drought conditions, respectively.Error bars represent standard errors of the mean at each time point.

Figure 4 .Figure 5 .
Figure 4. Potential biomarkers predicted by classical univariate ROC curve analysis.The left-side plot of the individual metabolite figure is showing ROC (receiver operating characteristics) curve with 95% confidence interval (shadowed) and the solid red dot indicates the optimal cutoff with the associated sensitivity and specificity values.The right-side box-and-whisker plot shows the distribution of abundance values of individual metabolite in control and drought sample with the optimal cutoff as a horizontal dotted red line.

Table 2 .
Important metabolites with their identifier number (KEGG ID/PubChemCID(*)/UNDP ID), molecular formula (MF), identified through Partial Least Square Discrepant Analysis (PLS-DA), Random Forest (RF) and Significant Analysis of Metabolites (SAM).(KEGGKyotoEncyclopedia of Genes and Genomes, UNDP = United Natural Products Database).American Journal of Plant Sciences the delta value of 1.6, FDR of 0.001 and with less than one (0.5) false positive.Similarly the most important metabolites were also identified by PLS-DA method based on the VIP score using five-component model.Random Forest classification ranked the important metabolites in order of decreasing prediction accuracy (Mean Decrease Accuracy) using 5000 trees (permutation) with an overall (OOB, out-of-bag) error of 0.0417.Overall the results were quite similar across all three methods.The top most important 60 metabolites which were identified by at least 2 different methods are shown in Table2.These important metabolites included different amino acids, sugars, organic acids/compounds, amines, fatty acids and others.
**Only variance for component 1 has been shown.

Table 3 .
Selected metabolites and their involvement in different pathways identified by Pathway Analysis of MetaboAnalist 3 using Oryzasativa japonica and Arabidopsis thaliana as the pathway libraries.
N-Acetyl-DL-glutamic acid Arginine and proline metabolism Aspartate Glycine, serine and threonine metabolism, Arginine and proline metabolism, Carbon fixation in photosynthetic organisms, Alanine, aspartate and glutamate metabolism, Aminoacyl-tRNA biosynthesis, Cyanoamino acid metabolism, beta-Alanine metabolism, Cysteine and methionine metabolism, Lysine biosynthesis, Nicotinate and nicotinamide metabolism Pyruvate Glycine, serine and threonine metabolism, Carbon fixation in photosynthetic organisms, Alanine, aspartate and glutamate metabolism, Cysteine and methionine metabolism, Citrate cycle (TCA cycle), Butanoate metabolism, Pyruvate metabolism, C5-Branched dibasic acid metabolism, Glycolysis or Gluconeogenesis, Valine, leucine and isoleucine biosynthesis, Pantothenate and CoA biosynthesis, Terpenoid backbone biosynthesis Others 328.1395 5.55 -2.6999 0.0096736 0.015996 2.02 American Journal of Plant Sciences Supplemental Table2.List of 89 significantly different levels of metabolites (ANOVA, P ≤ 0.05, Fisher's LSD) in the wheat flag leaves at 8 different time points (3 h, 6 h, 9 h, 12 h, 15 h, 18 h, 21 h, 24 h) during 24 hours of day under drought and control conditions with their P-value, false discovery rate (FDR) and significant pairs of time points.