Heng-An Lin1
Harrison R. Coker1
Julie A. Howe1
Malak M. Tfaily2
Elek M. Nagy3
Sanjay Antony-Babu3
Steve Hague1
A. Peyton Smith1*- 1Department of Soil and Crop Sciences, Texas A&M University and Texas A&M AgriLife Research, College Station, TX, United States
- 2Department of Environmental Science, University of Arizona, Tucson, AZ, United States
- 3Department of Plant Pathology and Microbiology, Texas A&M University and Texas A&M AgriLife Research, College Station, TX, United States
Root exudates comprise various primary and secondary metabolites that are responsive to plant stressors, including drought. As increasing drought episodes are predicted with climate change, identifying shifts in the metabolome profile of drought-induced root exudation is necessary to understand the molecular interactions that govern the relationships between plants, microbiomes, and the environment, which will ultimately aid in developing strategies for sustainable agriculture management. This study utilized an aeroponic system to simulate progressive drought and recovery while non-destructively collecting cotton (Gossypium hirsutum) root exudates. The molecular composition of the collected root exudates was characterized by untargeted metabolomics using Fourier-Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS) and mapped to the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Over 700 unique drought-induced metabolites were identified throughout the water-deficit phase. Potential KEGG pathways and KEGG modules associated with the biosynthesis of flavonoid compounds, plant hormones (abscisic acid and jasmonic acid), and other secondary metabolites were highly induced under severe drought, but not at the wilting point. Additionally, the associated precursors of these metabolites, such as amino acids (phenylalanine and tyrosine), phenylpropanoids, and carotenoids, were also mapped. The potential biochemical transformations were further calculated using the data generated by FT-ICR MS. Under severe drought stress, the highest number of potential biochemical transformations, including methylation, ethyl addition, and oxidation/hydroxylation, were identified, many of which are known reactions in some of the mapped pathways. With the application of FT-ICR MS, we revealed the dynamics of drought-induced secondary metabolites in root exudates in response to drought, providing valuable information for drought-tolerance strategies in cotton.
1 Introduction
Root systems release a large variety of compounds into the soil environment, commonly referred to as root exudates (Badri and Vivanco, 2009). Root exudates are mainly comprised of various primary and secondary metabolites, such as amino acids, sugars, phenolics, and plant hormones (Oburger and Jones, 2018). Root exudation is an important process that regulates interactions between plants, soil, and soil microorganisms in response to environmental stimuli (Bais et al., 2006; Baetz and Martinoia, 2014; Canarini et al., 2019). Plants may strategically modify the quantity and quality of their exudate profile in response to different abiotic and biotic stresses (Carvalhais et al., 2011; Gao et al., 2018; Li et al., 2019; Zhang et al., 2020), shaping the composition and activity of the rhizosphere microbiome, and promoting ecological feedback to plant hosts. For example, increased release of γ-aminobutyric acid (GABA) and carbohydrates was observed in the root exudates of P-deficient maize (Carvalhais et al., 2011). GABA has been known to regulate anion channels associated with malate efflux (Ramesh et al., 2015), and malate is one of the critical compounds that soybean (Glycine max) exudes to mobilize inorganic P in soil (Krishnapriya and Pandey, 2016). The release of carbohydrates into the rhizosphere has also been linked to enhancing germination and colonization of mycorrhizal fungi, which could improve P acquisition (Bücking et al., 2008; Etesami et al., 2021). The role of root exudates in inducing plant defense response or suppressing root diseases has also been reported in different hosts, such as wheat (Li et al., 2019), soybean (Gao et al., 2018), and tobacco (Zhang et al., 2020). While interest in root exudate research has been primarily focused on plant nutrient acquisition and the signal transductions between plant host and rhizosphere microbiome, there has been recent attention towards the effect of drought on root exudates, rhizosphere microbiomes under drought conditions, and these effects on ecosystem-scale responses (e.g., carbon sequestration, soil respiration, soil aggregation, soil organic matter decomposition, and plant-microbe symbiotic relationships) (Preece and Peñuelas, 2016; Karlowsky et al., 2018; de Vries et al., 2019; Williams and de Vries, 2020).
Several characteristics of root exudates undergo changes during drought. For example, drought-resistant maize (Zea mays) varieties exuded 58.2% more mucilage compared to drought-susceptible varieties (Nazari et al., 2023), and the role of mucilage has been linked to enhanced hydrologic conductance at the root-soil interface (Ahmed et al., 2014). The amino acid proline, an osmotic regulator, is known to accumulate in root exudates under drought in citrus (Citrus sp.) (Vives-Peris et al., 2017) and during recovery in pea (Pisum sativum) (Rubia et al., 2020). However, previous works have shown a range of responses in exudation rate and quantity of targeted compounds in response to drought. For example, increased exudation rates of carbon under drought have been reported in wheatgrass (Agropyron cristatum), mountain grasses, sunflower (Helianthus annuus), and holm oak (Quercus ilex) but not soybean (Glycine max) (Henry et al., 2007; Canarini et al., 2016; Karlowsky et al., 2018; Preece et al., 2018). Moreover, the variability of the drought effect on root exudates or the co-occurrence of multiple stressors [e.g., heat and drought stress (Tiziani et al., 2022)] presents a challenge in predicting how a particular plant species will respond to drought, as plant-soil interactions of different crops lack clear patterns associated with drought treatment or duration (Preece and Peñuelas, 2016). While a recent study hypothesized that root exudation patterns are linked to plant growth strategies and align with the ecosystem response to drought (Williams and de Vries, 2020), the use of various and inconsistent approaches in root exudate collection, different durations and severities of drought treatments, and varying analytical methods have made it challenging to compare the literature and draw conclusive results. In addition, the genetic effects (i.e., species, cultivars), plant growth stage effects, and other environmental factors (e.g., temperature, nutrient status, wind and light) have further complicated the understanding of the relationship between root exudation and drought response. Therefore, there is a need to reveal drought-induced exudation profiles for major crops with advanced technologies to identify potential drought-tolerance regulators.
Analytical techniques such as gas chromatography coupled to mass spectrometry (GC-MS), liquid chromatography-MS (LC-MS), and/or nuclear magnetic resonance spectroscopy (NMR) are often used to characterize the metabolome of plant root exudates (Fan et al., 2001; Luo et al., 2017; Gargallo-Garriga et al., 2018; Salem et al., 2022). Fourier-transform ion cyclotron resonance-MS (FT-ICR MS) is an ultrahigh-resolution analytical technique that has also been proven as a powerful tool for qualitative characterization of the metabolome of complex plant extracts (Maia et al., 2021) and for characterizing the root exudate composition of different plant species (Miao et al., 2020; Lohse et al., 2022; Ulrich et al., 2022). A shift of root exudate’s metabolome profile under different levels of drought has been identified in holm oak (Quercus ilex) using LC-MS (Gargallo-Garriga et al., 2018) and blue grama (Bouteloua gracilis) using GC-MS, NMR, and FT-ICR MS (Ulrich et al., 2022). Both studies demonstrated the dynamics of root exudate compositions and how they corresponded with the severity of drought. While the choice of analytical technology depends on multiple factors, such as high throughput capacity, optimal metabolome coverage, and other factors, FT-ICR MS yields high mass accuracy and resolution power, with a mass-to-charge ratio (m/z) range of 200 to 900 (Hsu et al., 2011; Park et al., 2013), which offers valuable insight into elemental and compound diversity of plant secondary metabolites (Takahashi et al., 2008; Maia et al., 2021). Further, by mapping the molecular formula to databases reveal potential functional pathways and direct the discovery of novel compounds as inherent advantages that are not as achievable with traditional approaches (Salem et al., 2022).
Drought has severe effects on cotton production. For example, in the Southwest region of the United States, a record crop abandonment rate 71% occurred in 2022 due to extreme and unexpected drought throughout the planting and growing season (Meyer et al., 2023). Compared with the previous season’s 12% abandonment rate, this led to a 40% decrease in cotton production in 2022 (Meyer et al., 2023). Therefore, disentangling the molecular feedback between root exudates and plant hosts in response to drought is crucial, as drought is predicted to increase in frequency due to climate change (IPCC, 2021). In this study, we deciphered the qualitative molecular characteristics of root exudates in upland cotton (Gossypium hirsutum), an important fiber crop commonly grown in water-restricted areas. We hypothesized the exudation profile would shift under progressive drought stress. More specifically, unique metabolites would be released during drought, which might be involved in plant stress response or act as signal molecules to interact with beneficial microbes for stress resilience. The objectives were to (1) characterize organic compounds in cotton root exudates under progressive drought stress using untargeted ultra-high resolution mass spectrometric technology (FT-ICR MS), (2) identify key metabolic pathways induced by drought, and (3) identify metabolic profile shifts throughout intensifying drought stress.
2 Materials and methods
2.1 Experimental design
An aeroponic growth chamber trial was performed in 2021. The pots were arranged in a completely randomized design with two treatments, progressive drought stress (drought) and well-watered control (control). Each treatment consisted of 10 pots and each pot had three plants. Treatments were applied for 11 days (water-deficit phase) once plants reached the match-head square stage (Elsner et al., 1979). We specifically targeted this growth stage for cotton because it is particularly vulnerable to drought, which can result in a severe reduction in yield and fiber quality (Zonta et al., 2017). The control treatment received full irrigation (130 mL hr-1)of full-strength Hoagland solution (Hoagland and Arnon, 1950) throughout the entire experiment, while the drought treatment received a 50% reduction in irrigation every 2 to 3 days that ended with 10% full irrigation (water-deficit phase). After 11 days of treatment, irrigation was returned to 100% for all pots for 7 days (recovery phase). Root exudates were non-destructively collected at 0, 2, 4, 7, 9, 11, 14, 16, and 18 days after the experiment was initiated, which includes a baseline (day 0), water-deficit phase (day 2, 4, 7, 9, and 11), and recovery phase (day 14, 16, and 18). Five out of ten pots were randomly selected for sample collection on each sampling day. A total of 90 samples were collected, which consisted of 2 treatments and 5 replicates for 9 sampling days (Figure 1).
Figure 1 Sample collection timeline and sample processing. (A) Root exudates were collected at 0, 2, 4, 7, 9, 11, 14, 16, and 18 days after the experiment was initiated, which includes a baseline (day 0), water-deficit phase (day 2, 4, 7, 9, and 11), and recovery phase (day 14, 16, and 18). (B) Root exudates collected from drought-treated and control-treated plants were filtered through 0.22-μm syringe filter, desalted, and concentrated by solid phase extraction, and analyzed with FT-ICR MS. Each treatment consisted of five pots with three plants per pot. This figure was created with BioRender.com.
2.2 Plant materials, hydroponic, and aeroponic systems
A moderately drought susceptible upland cotton (Gossypium hirsutum, cv. TAM 421) was chosen for this study. Cotton seeds were surface sterilized with 1 min of 70% EtOH, 3 min of 2% NaOCl, and rinsed with sterile distilled water 3 times. After the seeds germinated on the germination papers, three seedlings were placed in individual mesh net cups and supported with hydroponic sponges in each pot (~3.8 L). The hydroponic system was set up with an air pump equipped with two 10-way splitters. Each valve on splitter was connected to an air stone with airline tubing for each pot. Nutrients were supplied with a full-strength Hoagland solution (Hoagland and Arnon, 1950), where the nutrient solution was replaced twice a week, maintaining a pH range of 6.0-6.5. Plants were grown in a hydroponic system with constant aeration with dissolved oxygen levels ranging between 7.7 and 8.0 mg L-1 until 3 to 4 true leaves had developed. Then, the plants were transferred to an aeroponic system with 3 plants per pot. The set up of aeroponic system and the misting cycle was followed by Lin et al. (2022, preprint). In brief, the nutrient solution was pressurized using a diaphragm pump with a built-in pressure shut-off switch set to 827.3 kpa (120 psi). The pressurized solution was then temporarily stored in a 0.75 L pressure accumulator tank before it was delivered to each pot. The misting cycle was controlled by digital timers. For each pot (~7.6 L) in the aeroponic system, a spray ring was installed at the bottom (Figure 1B), composed of 3 nozzles (0.1 mm diameter). Misting was set for 10 sec, and time between misting events ranged between 5 min (well-watered) and 50 min (extreme drought). Environmental conditions in the walk-in growth chamber (EGC, USA) were set as 12-hour light/12-hour dark photoperiod. The air temperature and relative humidity settings were based on the local weather station to mimic the field environment (Supplementary Table 1). A complete description of plant growth conditions, aeroponic system design, and plant performance were reported in Lin et al. (2022, preprint). In brief, the drought symptoms (less than 5% of leaf affected) were first visually observed on day 4 of the drought treatment. The symptoms reached the greatest severity (over 66% of leaf affected) on day 9 and were approaching the wilting point on day 11. The plants in the drought treatment were visually fully recovered by the end of recovery phase.
2.3 Root exudate collection
Root exudates were collected from the whole root system in the aeroponic system at 14:00-16:00 on each collection day (Lin et al., 2022, preprint). The choice of the sampling time was based on our preliminary studies, which revealed a trend of higher concentrations of a targeted compound, abscisic acid, during the afternoon collections. In brief, roots were rinsed with a sampling solution (0.05 mM CaCl2, pH 6.0-6.5) for 2 min using the aeroponic system with a clean mist collection container. The choice of a dilute CaCl2 solution was designed to reduce the drastic change in ionic strength and mimic the soil solution, which is typically dominated by Ca2+ ions. The collected exudates were then filtered through a 22-μm polyethersulfone syringe filter and aliquoted into 1.5 mL centrifuge tubes. Samples were then stored at -80°C until further analysis.
2.4 Sample preparation
The samples were desalted and concentrated by solid phase extraction (SPE) using the method by Dittmar et al. (2008). The SPE cartridge (Bond Elut PPL, 100 mg, Agilent, CA, USA) was rinsed with 3 mL of methanol (one cartridge volume) to activate the column. A 1 mL aliquot of root exudates was first diluted with 14 mL of sterile double-distilled water and acidified with 1 M HCl until the pH was between 2 to 3. The samples were then vacuumed through the SPE cartridge under 170 mbar. Prior to elution, a total of 30 mL of 0.01 M HCl was rinsed through the cartridge to remove possible salt molecules and interferents. After allowing the sorbents to air dry, samples were eluted with 1.5 mL of methanol. Extracted samples were stored at -20°C until analysis.
2.5 Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS)
Extracted samples were analyzed with a 9.4 Tesla Bruker FT-ICR spectrometer located at the University of Arizona. A standard Bruker electrospray ionization (ESI) source was used to generate negatively charged molecular ions, where samples were introduced via direct infusion to the ESI source. 144 scans were averaged for each sample and internally calibrated using an organic matter homologous series separated by 14 Da (CH2 groups). Data Analysis software (BrukerDaltonik version 4.2) was used to convert raw spectra to a list of m/z values applying FTMS peak picker module with a signal-to-noise ratio (S/N) threshold set to 7 and absolute intensity threshold to the default value of 100. Putative chemical formulae were then assigned using Formularity (Tolić et al., 2017) software as previously described in Tfaily et al. (2018). Chemical formulae were assigned based on the following criteria: S/N > 7 and mass measurement error < 1 ppm, taking into consideration the presence of C, H, O, N, S and P and excluding other elements. The data produced by FT-ICR MS (peak masses, peak intensities, and metabolic molecular formula) were then processed through Metabodirect pipeline (Ayala-Ortiz et al., 2023). Peak masses (m/z) < 200 and > 900 were filtered using -m option to target secondary metabolites (Hsu et al., 2011; Park et al., 2013). Data was z-score normalized using –norm_method option. The quality control steps, including 13C isotope filtering and error filtering (0.5 ppm), followed the default setting. Masses must be identified in at least two samples to be included in the analysis. Elemental types (CHO, CHON, CHONP, CHONS, CHONSP, CHOP, CHOS, and CHOSP) were assigned for the filtered masses. KEGG analysis and mass difference analysis were performed using -k and -t options, respectively. Mass differences were used to infer potential biochemical transformations. Potential functional characteristics of the assigned molecular formulas were obtained by mapping to the KEGG database. Despite the molecular formula being the only criterion used for mapping to KEGG, it is important to acknowledge that these assignments are tentative. As noted, there is the possibility of multiple KEGG metabolites, particularly isomers, having the same molecular formula. Therefore, the molecular formula assignments should be interpreted with caution. The drought-induced unique metabolites in water-deficit phase (day 2 to 11) were characterized by comparing with control-treated samples excluding the drought-unique peaks identified at baseline (day 0). The KEGG pathways that were associated with human metabolisms were manually removed from the annotation list to generate the count table. The number of drought-unique metabolites were further visualized with UpSet plot (Lex et al., 2014) to identify the intersections between days (day 2 to 11). The figures were generated with R 4.2.0 (R core team, 2021) or GraphPad Prism v.6.0 (GraphPad Software, La Jolla, CA).
2.6 Statistical analysis for elemental compositions
The abundance matrix of elemental compositions was used to compute Euclidean distance matrices using vegdist function in R package vegan (Oksanen et al., 2020). The data were subset by experimental phases (i.e., baseline, water-deficit, and recovery) to generate each distance matrix. Permutational multivariate analysis of variance (PERMANOVA) analysis was performed for each experimental phase to assess the experiment effect using the adonis2 function in R package vegan (Oksanen et al., 2020). In a PERMANOVA model, the response variable was the distance matrix, and the explanatory variables were treatment (drought and control), days of treatment, and the interaction of treatment × days. Principal component analysis (PCA) was performed to detect the patterns in the elemental compositions using prcomp function in R package stats (R core team, 2021).
3 Results
3.1 Drought-induced unique metabolites were identified in cotton root exudates during water-deficit phase
A total number of 33,870 m/z were identified by FT-ICR MS. After the data was filtered using quality control steps (Ayala-Ortiz et al., 2023), a total of 13,033 m/z remained. After formula assignment, there were 3,985 metabolites with an assigned molecular formula. Overall, the total number of metabolites increased by 78% as drought stress progressed from day 2 to 11 (Figure 2). The highest number of total metabolites and drought-unique metabolites were identified under severe drought (day 9), the time point prior to wilting point (day 11) (Figure 2). Up to 70% of metabolites were shared between drought and control treated plant during water-deficit phase (Figure 2D). During the water-deficit phase, there were two discrete sampling days when the total number of metabolites increased (Figures 2A, C). The first increase occurred on day 4 of the drought treatment when the initial drought symptoms visually appeared. The second increase was observed on day 9, when the plants reached severe drought. Under drought, there was a reduction of the total number of metabolites at day 11 (i.e., near permanent wilting point) and a decrease throughout the recovery phase (day 14-18) (Figures 2A, B). The described trend in number of metabolites corresponded to the drought-unique metabolites. When grouping the data by phase, a greater total number of metabolites and drought-unique metabolites were observed in the water-deficit phase compared to day 0 and recovery phase (Figures 2B, D). Unique metabolites at baseline were identified between drought and control samples and those metabolites were excluded in some of the following analysis (Figure 2).
Figure 2 Stacked bar chart showing the number and percentage of shared and unique metabolites between control and drought treatments (n=5) throughout baseline (day 0), water-deficit (day 2-11) and recovery phases (day 14-18). Number of metabolites (assigned molecular formula) by day (A) and phase (B). Percentage of metabolites by day (C) and phase (D).
3.2 Elemental composition in cotton root exudates
PERMANOVA showed no differences between treatment, days, and treatment × days interaction in elemental composition (i.e., CHO, CHON, CHONP, CHONS, CHONSP, CHOP, CHOS, and CHOSP) when all compounds were used. Treatment had a nearly significant effect on elemental composition (p-value=0.07) during the water-deficit phase (Table 1). PCA analysis of elemental composition showed a minor effect of treatment during water-deficit phase (Figure 3) but not at baseline or during the recovery phase (Supplementary Figure 1). Although high percentages of explained variance in PC1 (72.4%) has been observed, the treatment effect on elemental compositions during water-deficit phase was observed at PC2 with 11.8% explained variance and the variation was mostly driven by CHO and CHON (Figure 3).
Table 1 Permutational multivariate analysis of variance (PERMANOVA) results of elemental composition for baseline (day 0), water deficit (day 2-11), and recovery phase (day 14-18) data (n=5).
Figure 3 Principal Component Analysis (PCA) biplot showing the variation elemental compositions of root exudates during water-deficit phase (day 2-11, n=5). The arrows were labeled by the elemental types.
3.3 Distribution of drought-induced metabolites in the water-deficit phase
The unique peaks identified from the drought-treated samples were further visualized using an UpSet plot (Figure 4). There were 345 metabolites consistently identified throughout the water-deficit phase. The number of distinct metabolites identified at individual days ranged from 118 to 481 with the highest number identified at day 9 in the drought treatment (Figure 4). While looking at the different intersections across the water-deficit phase, there were > 500 unique metabolites identified under severe drought stress (day 9-11) and < 100 metabolites uniquely identified at early drought stress (day 2-4) (Figure 4).
Figure 4 UpSet plot of drought-unique peaks (assigned molecular formula) by each day during water-deficit phase (day 2-11, n=5). UpSet plots visualize set intersections in a matrix layout. The bar plot on the left indicates the total number of unique peaks identified in drought-treated samples (D) at day 2 to 11. The bar chart at the top right and the number on top of each bar indicate the intersection sizes across different days. The connected dots at the bottom right indicate the days used for each intersection.
3.4 Functional characteristics of the drought-induced unique metabolites
KEGG annotation analysis identified a total of 68 pathways and over 17 modules that were induced under water-deficit phase (Tables 2, 3; Supplementary Tables 2, 3). In the drought treatment, 30 out of 68 pathways were consistently observed from day 2 to 11 including flavonoid biosynthesis, carotenoid biosynthesis, phenylpropanoid biosynthesis, and plant hormones biosynthesis (Table 2). A trend of increased counts was observed among those consistently induced pathways and modules (Tables 2, 3). The highest counts for drought-induced unique metabolites across most pathways were primarily observed on day 9 (severe drought) followed by a reduction at day 11 (near wilting point) (Table 2). Several pathways that might be associated with microorganisms, such as phosphotransferase system (PTS), biofilm formation and quorum sensing, and biosynthesis of several antibiotics were also identified. It is important to note that many assigned molecular formulas were mapped to higher KEGG hierarchies “metabolic pathway” and “biosynthesis of secondary metabolites” without further classification (Table 2).
Table 2 Potential KEGG pathways with number of drought-induced unique metabolites in the water-deficit phase (day 2-11, n=5).
Table 3 Potential KEGG modules with number of drought-induced unique metabolites in the water-deficit phase (day 2-11, n=5).
3.5 Higher biochemical transformations were identified under water-deficit plants
A total of 99 potential biochemical transformations were identified in this study (Table 4; Supplementary Table 4). Over 97% of the identified transformations had the highest number of counts in the drought-treated exudates at day 9, including methylation (-H), oxidation/hydroxylation (-H), hydrogenation/dehydrogenation, ethyl addition (-H2O), etc. (Table 4; Supplementary Table 4). The increase of transformation along the progressive drought stress showed a similar trend that was identified in the mass distribution and KEGG annotation results.
Table 4 Potential biochemical transformations identified throughout water-deficit phase (day 2-11, n=5) as the total counts of assigned molecular formula by day.
4 Discussion
The use of aeroponics systems allows for non-destructive root exudate collection under continuous treatments, eliminating the interference from soil particle absorption, root damage, and microbial decomposition (Oburger and Jones, 2018). In aeroponics, the abundance of root hairs and the root structure may differ compared to plants grown in soil (Eldridge et al., 2020). However, it has been reported that in cassava (Manihot esculenta), the genotype response to plant hormones remains consistent between aeroponic and field experiment (Selvaraj et al., 2019). This suggests that the plant’s response in aeroponics can be used to predict its response in the field. Nevertheless, further investigations are necessary to understand the differences in root exudate profiles in various cultivation environments. To the best of our understanding, this was the first time that drought-induced metabolome in cotton root exudates have been characterized, mapped to KEGG databases, and reported.
Previous work has shown that the concentration of primary metabolites in plant roots, such as methanol, glucose, small aminos acids, and other low molecular weight compounds changes under drought (Honeker et al., 2022). Since we focused on compounds with the m/z range of 200 to 900, we were not able to capture most of the low molecular weight compounds, which is also an important component of root exudates. Our research specifically focused on plant secondary metabolites in response to drought. The KEGG and biochemical transformation analysis suggested that changes to the drought-induced metabolome occur in a targetable groups of metabolic pathways (Wang et al., 2019), and the reactions take place intensively when plant experience severe drought. Although PERMANOVA was not able to detect differences in elemental composition between treatments and days of treatment during either the water-deficit or recovery phase (Table 1). This lead to the possibility that the plant alters the structure of metabolites without altering their molecular formula, which are unable to be captured by FT-ICR MS technology (Wang et al., 2019). The high proportion of shared metabolites between control and drought-treated plants (Figure 2D) might be associated with plant development, potentially explained why changes in elemental composition induced by drought were not able to be detected via PERMANOVA analysis. It has been shown that root exudation profile has larger differences between plant species; however, those differences are not only interspecific, but also intraspecific among individuals, and change over time in response to abiotic stresses (Gargallo-Garriga et al., 2018; Preece and Peñuelas, 2020). Therefore, it was not surprising that unique metabolites at baseline were identified between drought-stressed and control samples (Figure 2). To account for variation at baseline, those unique metabolites at baseline were excluded from the KEGG annotation analysis.
4.1 Roles of flavonoid compounds in root exudates under drought
During water-deficit phase, > 100 assigned molecular formula in drought-induced root exudates were mapped to pathways associated with the biosynthesis of flavonoids, including isoflavonoid, flavone, and flavonol (Table 2; Supplementary Table 2). Flavonoids are a group of plant-derived secondary metabolites that consist of > 10,000 compounds (Winkel-shirley, 2001; Sugiyama and Yazaki, 2014). The biosynthesis of flavonoids starts with phenylalanine and malonyl-CoA as the direct precursors. Then it undergoes different chemical reactions, such as hydroxylation, acylation, methylation, malonylation, and prenylation, which results in nine major subgroups (e.g., isoflavonoids, flavones, flavonols, flavandiols, and condensed tannins) having diverse structures and functions (Winkel-shirley, 2001; Sugiyama and Yazaki, 2014). Therefore, it was not surprising that we also observed the pathways associated with the biosynthesis of phenylalanine, tyrosine, and tryptophan; biosynthesis of phenylpropanoids; and KEGG modules involved with biosynthesis of malonyl-CoA under drought conditions (Table 2; Supplementary Table 3). In addition, a large number of methylation and malonylation transformations were also identified in drought-treated samples (Table 4; Supplementary Table 4).
Under drought conditions, a greater concentration of flavonoids have been identified in root exudates in holm oak (Gargallo-Garriga et al., 2018) and pearl millet (Pennisetum glaucum; Ghatak et al., 2022). Flavonoids have been reported as an important player in drought tolerance by directly functioning as antioxidants or serving as signaling molecules for the symbiosis of plants and microbes (Hassan and Mathesius, 2012). In our study, several masses were mapped to the KEGG module related to “daidzein → medicarpin” reactions (Table 3; Supplementary Table 3). Daidzein is an isoflavone compound reported to be secreted by soybean roots, as signaling molecules mediating communication between plants and nitrogen-fixing bacteria (Sugiyama et al., 2007; Okutani et al., 2020). It is unclear if cotton also uses similar strategies to attract nitrogen-fixing bacteria. However, several nitrogen-uptake-related bacteria, such as Mesorhizobium, Sinorhizobium, and Rhizobium, have been identified in the cotton rhizosphere (Qiao et al., 2017), which provide potential directions for future studies. Further investigations are required to characterize and quantify the flavonoid compounds in cotton root exudates and the associated microbes.
4.2 Plant hormones in root exudates
It has been well documented that multiple ubiquitous plant hormones, including abscisic acid (ABA) and jasmonic acid (JA), are coordinated in response to drought stress (Wang et al., 2021; El Sabagh et al., 2022). In this study, drought treatment induced the pathways and modules involved in the biosynthesis of ABA and JA, as well as their associated precursors, carotenoids (Wasilewska et al., 2008) and alpha-linolenic acid (Wasternack and Hause, 2013) (Tables 2, 3). In a previous study, we were able to detect ABA in collected root exudates using GC-MS, although the differences between drought and control treatment were only detected at certain days during water deficit phase (Lin et al., 2022, preprint). We also identified the pathway associated with ATP-binding cassette (ABC) transporters induced under drought, which have been reported to transport phytohormones and root exudation (Badri and Vivanco, 2009; Kang et al., 2011). Because our root exudation collection method was non-destructive compared to other sampling approaches (Oburger and Jones, 2018), we assume that cotton undergoes active transport of those plant hormones into the rhizosphere in order to regulate hormone concentrations inside the root tissues in response to drought, as previously reported in upland rice (Oryzae sativa) (Shi et al., 2014). However, it would be necessary to identify the associated transporters to support this assumption and further understand the mechanisms behind this finding.
4.3 Antibiotics identified in root exudates potentially associated with root endophytes
Interestingly, pathways associated with the biosynthesis of antibiotics (e.g., novobiocin, enediyne antibiotics, and phenazine) were identified under drought, even in soil-less environments. Those antibiotics have been reported to be produced by a variety of saprophytic or endophytic microorganisms (Steffensky et al., 1998; Laursen and Nielsen, 2004; Chen et al., 2010; Li et al., 2015). Streptomyces, which includes both endophytic and saprophytic species, has been known to produce novobiocin and enediyne antibiotics (Steffensky et al., 1998; Chen et al., 2010). Phenazine is also a common antibiotic that can be isolated from Streptomyces and Pseudomonas (Laursen and Nielsen, 2004; Li et al., 2015) that is known to be involved in biofilm formation, which is important process in drought tolerance in plant roots (Mahmoudi et al., 2019). Our growth chamber experiment was not performed in a sterile environment, so by using a soilless aeroponic setup, is it likely that those antibiotics might be produced by host endophytes. Indeed, studies have shown that Streptomyces was enriched in the root under drought and expected to be associated with improving plant fitness (Xu and Coleman-Derr, 2019). Pseudomonas also has been reported as one of the dominant endophytic organisms in cotton roots (Shi et al., 2020). Our current experiment design is not able to differentiate the metabolites that were released by plants or endophytes; however, incorporating isotope labeling techniques could further reveal the mechanisms between plants and microbes under drought.
4.4 Major shift of metabolomic profile during severe drought, but not wilting point
A major shift in the metabolomic profile occurred in our study under severe drought prior to the plant reaching wilting point. Our findings provide temporal responses to progressive drought with activation of various metabolic pathways, which provide important information on relevant sampling time points to further characterize and validate the drought-response mechanisms in cotton root exudates. The transition from severe drought stress to wilting point (day 9 to 11) may be related to cell death, giving insight into the complicated senescence of plants and how root exudates respond to plant death, which has limited attention in recent studies (Prudence et al., 2021). One pathway, ferroptosis, associated with cell death, was identified exclusively on day 9 (Table 2). Ferroptosis is an iron-dependent and highly regulated cell death process (Riyazuddin and Gupta, 2021; Distéfano et al., 2022). In Arabidopsis thaliana, this programmed cell death process was induced by heat stress and was mainly dependent on reactive oxygen species (ROS)-mediated lipid peroxidation (Distéfano et al., 2017). Currently, there are few studies on this newly discovered pathway and no reports indicating whether drought stress will induce this pathway. However, many abiotic stresses, including drought, also induce the accumulation of ROS in plant tissues (Lee and Park, 2012). Shared exudation patterns among heat, drought, and the combination of both stresses have been reported in maize (Tiziani et al., 2022), indicating there are shared responses when plants encounter these types of stresses. The small amount (only 17) of assigned molecular formula that were mapped to the ferroptosis cell death pathway, which were only identified during severe drought (day 9) in our study, suggests that the plant might release “final warning stress signals” prior to wilting point, but it also suggests there might be minor confounding responses between drought and cell death at this collection timepoint.
5 Conclusions
Our study suggested that cotton released unique metabolites in root exudates depending on drought severity. With the use of FT-ICR MS, we were able to characterize all drought-induced metabolites with a m/z range of 200 to 900. The drought-induced metabolome occurred in a small and targetable group of metabolic pathways. To further isolate the specific compounds, targeting the plant under severe drought instead of the wilting point is recommended. By utilizing a non-destructive sampling approach, we were able to capture the temporal changes of exudation profiles of the same plant without the interference of soil particles. Our findings advance the fundamental characteristics of root exudates and plant drought-tolerance mechanisms belowground. Future research should focus on: (1) characterizing the unclassified metabolites, particularly secondary metabolites; (2) utilizing targeted metabolomic approaches (e.g., LC-MS, GC-MS, or HPLC) to quantify the concentration of several drought-induced candidate groups and confirm metabolite annotations, such as flavonoids, to further verified their importance to the plant drought response; and (3) incorporate the results in greenhouse or field trials to further understand the consequence of drought on rhizosphere microbiome and plant performance.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
H-AL contributed to data curation, formal analysis, investigation, methodology, visualization, writing – original draft, writing – review & editing. HRC contributed to data curation, formal analysis, investigation, methodology, writing – review & editing. JAH contributed to conceptualization, funding acquisition, investigation, methodology, resources, supervision, writing – review & editing. MMT contributed to SPE extraction, FT-ICR MS analysis and data curation, writing – review & editing. EMN contributed to methodology, writing – review & editing. SA-B contributed to resources, writing – review & editing. SH contributed to germplasm resources, writing – review & editing. APS contributed to conceptualization, funding acquisition, investigation, methodology, project administration, resources, supervision, writing – review & editing. All authors contributed to the article and approved the submitted version.
Funding
This project was supported by The National Institute of Food and Agriculture within United States Department of Agriculture [grant numbers: 2020-67019-31799, Hatch project 1018999].
Acknowledgments
We thank Chia-Wei Lin and Dr. Youjun Deng for assisting with the processing of biomass samples. We would like to acknowledge Dr. Krishna Parsawar and Dr. Lawrence Walker from the Analytical and Biological Mass Spectrometry (ABMS) Core Facility at the University of Arizona supported in part by the RII (Research, Innovation, and Impact) and TRIF (Technology and Research Initiative Fund) initiative.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1244591/full#supplementary-material
References
Ahmed, M. A., Kroener, E., Holz, M., Zarebanadkouki, M., Carminati, A. (2014). Mucilage exudation facilitates root water uptake in dry soils. Funct. Plant Biol. 41, 1129–1137. doi: 10.1071/FP13330
Ayala-Ortiz, C., Graf-Grachet, N., Freire-Zapata, V., Fudyma, J., Hildebrand, G., AminiTabrizi, R., et al. (2023). MetaboDirect: an analytical pipeline for the processing of FT-ICR MS-based metabolomic data. Microbiome 11, 28. doi: 10.1186/s40168-023-01476-3
Badri, D. V., Quintana, N., El Kassis, E. G., Kim, H. K., Choi, Y. H., Sugiyama, A., et al. (2009). An ABC transporter mutation alters root exudation of phytochemicals that provoke an overhaul of natural soil microbiota. Plant Physiol. 151, 2006–2017. doi: 10.1104/pp.109.147462
Badri, D. V., Vivanco, J. M. (2009). Regulation and function of root exudates. Plant. Cell Environ. 32, 666–681. doi: 10.1111/j.1365-3040.2009.01926.x
Baetz, U., Martinoia, E. (2014). Root exudates: The hidden part of plant defense. Trends Plant Sci. 19, 90–98. doi: 10.1016/j.tplants.2013.11.006
Bais, H. P., Weir, T. L., Perry, L. G., Gilroy, S., Vivanco, J. M. (2006). The role of root exudates in rhizosphere interactions with plants and other organisms. Annu. Rev. Plant Biol. 57, 233–266. doi: 10.1146/annurev.arplant.57.032905.105159
Berg, G., Rybakova, D., Fischer, D., Cernava, T., Vergès, M. C. C., Charles, T., et al. (2020). Microbiome definition re-visited: old concepts and new challenges. Microbiome 8, 1–22. doi: 10.1186/s40168-020-00875-0
Bücking, H., Abubaker, J., Govindarajulu, M., Tala, M., Pfeffer, P. E., Nagahashi, G., et al. (2008). Root exudates stimulate the uptake and metabolism of organic carbon in germinating spores of Glomus intraradices. New Phytol. 180, 684–695. doi: 10.1111/j.1469-8137.2008.02590.x
Canarini, A., Kaiser, C., Merchant, A., Richter, A., Wanek, W. (2019). Root exudation of primary metabolites: mechanisms and their roles in plant responses to environmental stimuli. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.00157
Canarini, A., Merchant, A., Dijkstra, F. A. (2016). Drought effects on Helianthus annuus and Glycine max metabolites: from phloem to root exudates. Rhizosphere 2, 85–97. doi: 10.1016/j.rhisph.2016.06.003
Carvalhais, L. C., Dennis, P. G., Fedoseyenko, D., Hajirezaei, M. R., Borriss, R., von Wirén, N. (2011). Root exudation of sugars, amino acids, and organic acids by maize as affected by nitrogen, phosphorus, potassium, and iron deficiency. J. Plant Nutr. Soil Sci. 174, 3–11. doi: 10.1002/jpln.201000085
Chen, Y., Yin, M., Horsman, G. P., Huang, S., Shen, B. (2010). Manipulation of pathway regulation in Streptomyces globisporus for overproduction of the enediyne antitumor antibiotic C-1027. J. Antibiot. (Tokyo). 63, 482–485. doi: 10.1038/ja.2010.55
Chhaya, Yadav, B., Jogawat, A., Gnanasekaran, P., Kumari, P., Lakra, N., et al. (2021). An overview of recent advancement in phytohormones-mediated stress management and drought tolerance in crop plants. Plant Gene 25, 100264. doi: 10.1016/j.plgene.2020.100264
de Vries, F. T., Williams, A., Stringer, F., Willcocks, R., McEwing, R., Langridge, H., et al. (2019). Changes in root-exudate-induced respiration reveal a novel mechanism through which drought affects ecosystem carbon cycling. New Phytol. 224, 132–145. doi: 10.1111/nph.16001
Distéfano, A. M., López, G. A., Bauer, V., Zabaleta, E., Pagnussat, G. C. (2022). Ferroptosis in plants: regulation of lipid peroxidation and redox status. Biochem. J. 479, 857–866. doi: 10.1042/BCJ20210682
Distéfano, A. M., Martin, M. V., Córdoba, J. P., Bellido, A. M., D’Ippólito, S., Colman, S. L., et al. (2017). Heat stress induces ferroptosis-like cell death in plants. J. Cell Biol. 216, 463–476. doi: 10.1083/jcb.201605110
Dittmar, T., Koch, B., Hertkorn, N., Kattner, G. (2008). A simple and efficient method for the solid-phase extraction of dissolved organic matter (SPE-DOM) from seawater. Limnol. Oceanogr. Methods 6, 230–235. doi: 10.4319/lom.2008.6.230
Eldridge, B. M., Manzoni, L. R., Graham, C. A., Rodgers, B., Farmer, J. R., Dodd, A. N. (2020). Getting to the roots of aeroponic indoor farming. New Phytol. 228, 1183–1192. doi: 10.1111/nph.16780
El Sabagh, A., Islam, M. S., Hossain, A., Iqbal, M. A., Mubeen, M., Waleed, M., et al. (2022). Phytohormones as growth regulators during abiotic stress tolerance in plants. Front. Agron. 4. doi: 10.3389/fagro.2022.765068
Elsner, J. E., Smith, C. W., Owen, D. F. (1979). Uniform stage descriptions in upland cotton. Crop Sci. 19, 361–363. doi: 10.2135/cropsci1979.0011183X001900030021x
Etesami, H., Jeong, B. R., Glick, B. R. (2021). Contribution of arbuscular mycorrhizal fungi, phosphate–solubilizing bacteria, and silicon to P uptake by plant. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.699618
Fan, T. W.-M., Lane, A. N., Shenker, M., Bartley, J. P., Crowley, D., Higashi, R. M. (2001). Comprehensive chemical profiling of gramineous plant root exudates using high-resolution NMR and MS. Phytochemistry 57, 209–221. doi: 10.1016/S0031-9422(01)00007-3
Fraser, C. M., Chapple, C. (2011). The phenylpropanoid pathway in Arabidopsis. Arab. B. 9, e0152. doi: 10.1199/tab.0152
Gao, X., Zhang, S., Zhao, X., Wu, Q. (2018). Potassium-induced plant resistance against soybean cyst nematode via root exudation of phenolic acids and plant pathogen-related genes. PloS One 13, 1–13. doi: 10.1371/journal.pone.0200903
Gargallo-Garriga, A., Preece, C., Sardans, J., Oravec, M., Urban, O., Peñuelas, J. (2018). Root exudate metabolomes change under drought and show limited capacity for recovery. Sci. Rep. 8, 1–15. doi: 10.1038/s41598-018-30150-0
Ghatak, A., Schindler, F., Bachmann, G., Engelmeier, D., Bajaj, P., Brenner, M., et al. (2022). Root exudation of contrasting drought-stressed pearl millet genotypes conveys varying biological nitrification inhibition (BNI) activity. Biol. Fertil. Soils. 58, 291–306. doi: 10.1007/s00374-021-01578-w
Hassan, S., Mathesius, U. (2012). The role of flavonoids in root-rhizosphere signaling: opportunities and challenges for improving plant-microbe interactions. J. Exp. Bot. 63, 3429–3444. doi: 10.1093/jxb/err430
He, M., Ding, N.-Z. (2020). Plant unsaturated fatty acids: multiple roles in stress response. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.562785
Henry, A., Doucette, W., Norton, J., Bugbee, B. (2007). Changes in crested wheatgrass root exudation caused by flood, drought, and nutrient stress. J. Environ. Qual. 36, 904–912. doi: 10.2134/jeq2006.0425sc
Hoagland, D. R., Arnon, D. I. (1950). Preparing the nutrient solution. Water-Culture. Method. Grow. Plants Without. Soil 347, 29–31.
Honeker, L. K., Hildebrand, G. A., Fudyma, J. D., Daber, L. E., Hoyt, D., Flowers, S. E., et al. (2022). Elucidating drought-tolerance mechanisms in plant roots through 1H NMR metabolomics in parallel with MALDI-MS, and NanoSIMS imaging techniques. Environ. Sci. Technol. 56, 2021–2032. doi: 10.1021/acs.est.1c06772
Hsu, C. S., Hendrickson, C. L., Rodgers, R. P., McKenna, A. M., Marshall, A. G. (2011). Petroleomics: advanced molecular probe for petroleum heavy ends. J. Mass. Spectrom. 46, 337–343. doi: 10.1002/jms.1893
IPCC. (2021). Climate change 2021: The physical science basis. Contribution of working group I to the sixth assessment report of the intergovernmental panel on climate change [Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., et al (eds.)]. (Cambridge University Press). In Press
Kang, J., Park, J., Choi, H., Burla, B., Kretzschmar, T., Lee, Y., et al. (2011). Plant ABC transporters. Arab. B. 9, e0153. doi: 10.1199/tab.0153
Karlowsky, S., Augusti, A., Ingrisch, J., Akanda, M. K. U., Bahn, M., Gleixner, G. (2018). Drought-induced accumulation of root exudates supports post-drought recovery of microbes in mountain grassland. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.01593
Krishnapriya, V., Pandey, R. (2016). Root exudation index: screening organic acid exudation and phosphorus acquisition efficiency in soybean genotypes. Crop Pasture Sci. 67, 1096–1109. doi: 10.1071/CP15329
Laursen, J. B., Nielsen, J. (2004). Phenazine natural products: biosynthesis, synthetic analogues, and biological activity. Chem. Rev. 104, 1663–1686. doi: 10.1021/cr020473j
Lee, S., Park, C.-M. (2012). Regulation of reactive oxygen species generation under drought conditions in Arabidopsis. Plant Signal. Behav. 7, 599–601. doi: 10.4161/psb.19940
Lex, A., Gehlenborg, N., Strobelt, H., Vuillemot, R., Pfister, H. (2014). UpSet: visualization of intersecting sets. IEEE Trans. Vis. Comput. Graph. 20, 1983–1992. doi: 10.1109/TVCG.2014.2346248
Li, C.-X., Fu, X.-P., Zhou, X.-G., Liu, S.-W., Xia, Y., Li, N.-H., et al. (2019). Treatment with wheat root exudates and soil microorganisms from wheat/watermelon companion cropping can induce watermelon disease resistance against Fusarium oxysporum f. sp. niveum. Plant Dis. 103, 1693–1702. doi: 10.1094/PDIS-08-18-1387-RE
Li, Y., Han, L., Rong, H., Li, L., Zhao, L., Wu, L., et al. (2015). Diastaphenazine, a new dimeric phenazine from an endophytic Streptomyces diastaticus subsp. ardesiacus. J. Antibiot. (Tokyo). 68, 210–212. doi: 10.1038/ja.2014.124
Lin, H.-A., Coker, H. R., Park, S., Finlayson, S. A., Nagy1, E. M., Hague, S., et al. (2022). Nondestructive root exudate sampling using aeroponics throughout progressive drought stress and recovery. Prepr. (Version 1) available Res. Sq. doi: 10.21203/rs.3.rs-1886373/v1
Lohse, M., Santangeli, M., Steininger-Mairinger, T., Oburger, E., Reemtsma, T., Lechtenfeld, O. J., et al. (2022). The effect of root hairs on exudate composition: a comparative non-targeted metabolomics approach. Anal. Bioanal. Chem 415(5), 823–840. doi: 10.1007/s00216-022-04475-9
Luo, Q., Wang, S., Sun, L., Wang, H. (2017). Metabolic profiling of root exudates from two ecotypes of Sedum alfredii treated with Pb based on GC-MS. Sci. Rep. 7, 39878. doi: 10.1038/srep39878
Maia, M., Figueiredo, A., Cordeiro, C., Sousa Silva, M. (2021). FT-ICR-MS-based metabolomics: A deep dive into plant metabolism. Mass. Spectrom. Rev. 42, 1535–1556. doi: 10.1002/mas.21731
Mahmoudi, T. R., Yu, J. M., Liu, S., Pierson, L. S., Pierson, E. A. (2019). Drought-stress tolerance in wheat seedlings conferred by phenazine-producing rhizobacteria. Front. Microbiol. 10.
Meyer, L., Dew, T., Grace, M., Lanclos, K., MacDonald, S., Soley, G. (2023). The world and United States cotton outlook (United States Department of Agriculture, Agricultural Outlook Forum).
Miao, Y., Lv, J., Huang, H., Cao, D., Zhang, S. (2020). Molecular characterization of root exudates using Fourier transform ion cyclotron resonance mass spectrometry. J. Environ. Sci. (China). 98, 22–30. doi: 10.1016/j.jes.2020.05.011
Nazari, M., Bilyera, N., Banfield, C. C., Mason-Jones, K., Zarebanadkouki, M., Munene, R., et al. (2023). Soil, climate, and variety impact on quantity and quality of maize root mucilage exudation. Plant Soil 482, 25–38. doi: 10.1007/s11104-022-05669-x
Oburger, E., Jones, D. L. (2018). Sampling root exudates – Mission impossible? Rhizosphere 6, 116–133. doi: 10.1016/j.rhisph.2018.06.004
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2020) Vegan: Community ecology package. R package version 2. Available at: https://CRAN.R-project.org/package=vegan.
Okutani, F., Hamamoto, S., Aoki, Y., Nakayasu, M., Nihei, N., Nishimura, T., et al. (2020). Rhizosphere modelling reveals spatiotemporal distribution of daidzein shaping soybean rhizosphere bacterial community. Plant Cell Environ. 43, 1036–1046. doi: 10.1111/pce.13708
Park, K. H., Kim, M. S., Baek, S. J., Bae, I. H., Seo, S.-W., Kim, J., et al. (2013). Simultaneous molecular formula determinations of natural compounds in a plant extract using 15 T Fourier transform ion cyclotron resonance mass spectrometry. Plant Methods 9, 15. doi: 10.1186/1746-4811-9-15
Pichersky, E., Raguso, R. A. (2018). Why do plants produce so many terpenoid compounds? New Phytol. 220, 692–702. doi: 10.1111/nph.14178
Preece, C., Farré-Armengol, G., Llusià, J., Peñuelas, J. (2018). Thirsty tree roots exude more carbon. Tree Physiol. 38, 690–695. doi: 10.1093/treephys/tpx163
Preece, C., Peñuelas, J. (2016). Rhizodeposition under drought and consequences for soil communities and ecosystem resilience. Plant Soil 409, 1–17. doi: 10.1007/s11104-016-3090-z
Preece, C., Peñuelas, J. (2020). A return to the wild: root exudates and food security. Trends Plant Sci. 25, 14–21. doi: 10.1016/j.tplants.2019.09.010
Prudence, S. M. M., Newitt†, J. T., Worsley, S. F., Macey, M. C., Murrell, J. C., Lehtovirta-Morley, L. E., et al. (2021). Soil, senescence and exudate utilisation: characterisation of the Paragon var. spring bread wheat root microbiome. Environ. Microbiome. 16, 12. doi: 10.1186/s40793-021-00381-2
Qiao, Q., Wang, F., Zhang, J., Chen, Y., Zhang, C., Liu, G., et al. (2017). The variation in the rhizosphere microbiome of cotton with soil type, genotype and developmental stage. Sci. Rep. 7, 1–10. doi: 10.1038/s41598-017-04213-7
R Core Team. (2021). R: A language and environment for statistical computing. R Foundation forStatistical Computing Vienna, Austria. https://www.R-project.org/.
Ramesh, S. A., Tyerman, S. D., Xu, B., Bose, J., Kaur, S., Conn, V., et al. (2015). GABA signaling modulates plant growth by directly regulating the activity of plant-specific anion transporters. Nat. Commun. 6, 7879. doi: 10.1038/ncomms8879
Riyazuddin, R., Gupta, R. (2021). Plausible involvement of ethylene in plant ferroptosis: prospects and leads. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.680709
Rubia, M. I., Ramachandran, V. K., Arrese-Igor, C., Larrainzar, E., Poole, P. S. (2020). A novel biosensor to monitor proline in pea root exudates and nodules under osmotic stress and recovery. Plant Soil 452, 413–422. doi: 10.1007/s11104-020-04577-2
Salem, M. A., Wang, J. Y., Al-Babili, S. (2022). Metabolomics of plant root exudates: From sample preparation to data analysis. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.1062982
Selvaraj, M. G., Montoya-P, M. E., Atanbori, J., French, A. P., Pridmore, T. (2019). A low-cost aeroponic phenotyping system for storage root development: unravelling the below-ground secrets of cassava (Manihot esculenta). Plant Methods 15, 131. doi: 10.1186/s13007-019-0517-6
Shi, L., Guo, M., Ye, N., Liu, Y., Liu, R., Xia, Y., et al. (2014). Reduced ABA accumulation in the root system is caused by ABA exudation in upland rice (Oryza sativa L. var. Gaoshan1) and this enhanced drought adaptation. Plant Cell Physiol. 56, 951–964. doi: 10.1093/pcp/pcv022
Shi, Y., Yang, H., Chu, M., Niu, X., Huo, X., Gao, Y., et al. (2020). Endophytic bacterial communities and spatiotemporal variations in cotton roots in Xinjiang, China. Can. J. Microbiol. 67, 506–517. doi: 10.1139/cjm-2020-0249
Steffensky, M., Li, S.-M., Vogler, B., Heide, L. (1998). Novobiocin biosynthesis in Streptomyces spheroides: identification of a dimethylallyl diphosphate:4-hydroxyphenylpyruvate dimethylallyl transferase. FEMS Microbiol. Lett. 161, 69–74. doi: 10.1111/j.1574-6968.1998.tb12930.x
Sugiyama, A., Shitan, N., Yazaki, K. (2007). Involvement of a soybean ATP-binding cassette-type transporter in the secretion of genistein, a signal flavonoid in legume-rhizobium symbiosis. Plant Physiol. 144, 2000–2008. doi: 10.1104/pp.107.096727
Sugiyama, A., Yazaki, K. (2014). Flavonoids in plant rhizospheres: secretion, fate and their effects on biological communication. Plant Biotechnol. 31, 431–443. doi: 10.5511/plantbiotechnology.14.0917a
Takahashi, H., Kai, K., Shinbo, Y., Tanaka, K., Ohta, D., Oshima, T., et al. (2008). Metabolomics approach for determining growth-specific metabolites based on Fourier transform ion cyclotron resonance mass spectrometry. Anal. Bioanal. Chem. 391, 2769–2782. doi: 10.1007/s00216-008-2195-5
Tfaily, M. M., Hess, N. J., Koyama, A., Evans, R. D. (2018). Elevated CO2 changes soil organic matter composition and substrate diversity in an arid ecosystem. Geoderma 330, 1–8. doi: 10.1016/j.geoderma.2018.05.025
Tiziani, R., Miras-Moreno, B., Malacrinò, A., Vescio, R., Lucini, L., Mimmo, T., et al. (2022). Drought, heat, and their combination impact the root exudation patterns and rhizosphere microbiome in maize roots. Environ. Exp. Bot. 203, 105071. doi: 10.1016/j.envexpbot.2022.105071
Tolić, N., Liu, Y., Liyu, A., Shen, Y., Tfaily, M. M., Kujawinski, E. B., et al. (2017). Formularity: software for automated formula assignment of natural and other organic matter from ultrahigh-resolution mass spectra. Anal. Chem. 89, 12659–12665. doi: 10.1021/acs.analchem.7b03318
Ulrich, D. E. M., Clendinen, C. S., Alongi, F., Mueller, R. C., Chu, R. K., Toyoda, J., et al. (2022). Root exudate composition reflects drought severity gradient in blue grama (Bouteloua gracilis). Sci. Rep. 12, 12581. doi: 10.1038/s41598-022-16408-8
Vives-Peris, V., Gómez-Cadenas, A., Pérez-Clemente, R. M. (2017). Citrus plants exude proline and phytohormones under abiotic stress conditions. Plant Cell Rep. 36, 1971–1984. doi: 10.1007/s00299-017-2214-0
Wang, S., Alseekh, S., Fernie, A. R., Luo, J. (2019). The structure and function of major plant metabolite modifications. Mol. Plant 12, 899–919. doi: 10.1016/j.molp.2019.06.001
Wang, N., LI, L., Gao, W., Wu, Y., Yong, H., Weng, J., et al. (2018). Transcriptomes of early developing tassels under drought stress reveal differential expression of genes related to drought tolerance in maize. J. Integr. Agric. 17, 1276–1288. doi: 10.1016/S2095-3119(17)61777-5
Wang, X., Li, Q., Xie, J., Huang, M., Cai, J., Zhou, Q., et al. (2021). Abscisic acid and jasmonic acid are involved in drought priming-induced tolerance to drought in wheat. Crop J. 9, 120–132. doi: 10.1016/j.cj.2020.06.002
Wasilewska, A., Vlad, F., Sirichandra, C., Redko, Y., Jammes, F., Valon, C., et al. (2008). An update on abscisic acid signaling in plants and more. Mol. Plant 1, 198–217. doi: 10.1093/mp/ssm022
Wasternack, C., Hause, B. (2013). Jasmonates: biosynthesis, perception, signal transduction and action in plant stress response, growth and development. An update to the 2007 review in annals of botany. Ann. Bot. 111, 1021–1058. doi: 10.1093/aob/mct067
Williams, A., de Vries, F. T. (2020). Plant root exudation under drought: implications for ecosystem functioning. New Phytol. 225, 1899–1905. doi: 10.1111/nph.16223
Winkel-shirley, B. (2001). Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant phy. 126, 485–493. doi: 10.1104/pp.126.2.485
Xu, L., Coleman-Derr, D. (2019). Causes and consequences of a conserved bacterial root microbiome response to drought stress. Curr. Opin. Microbiol. 49, 1–6. doi: 10.1016/j.mib.2019.07.003
Zhang, C., Feng, C., Zheng, Y., Wang, J., Wang, F. (2020). Root exudates metabolic profiling suggests distinct defense mechanisms between resistant and susceptible tobacco cultivars against black shank disease. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.559775
Keywords: root exudates, untargeted metabolomics, drought, FT-ICR MS, upland cotton, nondestructive sampling, aeroponics
Citation: Lin H-A, Coker HR, Howe JA, Tfaily MM, Nagy EM, Antony-Babu S, Hague S and Smith AP (2023) Progressive drought alters the root exudate metabolome and differentially activates metabolic pathways in cotton (Gossypium hirsutum). Front. Plant Sci. 14:1244591. doi: 10.3389/fpls.2023.1244591
Received: 22 June 2023; Accepted: 14 August 2023;
Published: 30 August 2023.
Edited by:
María José García, University of Cordoba, SpainReviewed by:
Mauro Maver, Free University of Bozen-Bolzano, ItalySanna Sevanto, Los Alamos National Laboratory (DOE), United States
Copyright © 2023 Lin, Coker, Howe, Tfaily, Nagy, Antony-Babu, Hague and Smith. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: A. Peyton Smith, peyton.smith@ag.tamu.edu