2023
Molecular Composition of Soybean Seeds from Drought Susceptible and Tolerant Cultivars
Category:
Sustainable Production
Keywords:
Abiotic stressAgricultureLand Use Water supply
Lead Principal Investigator:
Troy Wood , Research Foundation for SUNY on behalf of University at Buffalo
Co-Principal Investigators:
Project Code:
SYBN 23 010
Brief Project Summary:
A researcher at the State University of New York at Buffalo has been collaborating with one of the world’s foremost soybean geneticists to compare two soybean varieties at the metabolic level. One of the varieties has been crossbred with soybean plants that aren’t native to North America to become more drought resistant—and the research has demonstrated that resistance. The variety, called “731” by researchers, inhibits the process that allows water molecules to escape plant leaves and has a higher chlorophyll content, among other characteristics.
Key Beneficiaries:
#agronomists, #Soybean producers
Unique Keywords:
#chlorophyll , #drought, #drought resistant, #environmental stress, #genetics, #metabolic, #water
Information And Results
Project Summary

Agricultural crops can endure a matrix of stress resultant from a variety of sources including biotic or abiotic stressors such as drought, flooding, salinity, or nutrient availability. Among the different sources of stress plants can undergo, varying levels of water deficiency and drought have the most prolific and detrimental effect to agricultural farms on the global and national scale. Numerous plant traits have been identified for the potential of improving the performance of drought-affected crops, mainly through conservation of water, with more recent works identifying the importance of slow canopy wilting (SW) phenotypes for their potential stress tolerance in water deficient environments. Legumes have a particular intolerance to water deficiency in the early stages of growth and flowering, where a decrease in water availability by half can result in up to a loss of half the expected yields.

An exotic soybean germplasm, plant introduction (PI) 567731 in maturity group III (MG III), was identified to consistently express the SW phenotype in the field compared to the drought sensitive cultivar Pana; PI 567731 showed lower yield loss than Pana under drought stress with greater than 13% more yield index (yield under rain-fed/ yield under irrigation). PI 567731 uses significantly less water under drought, and this water conservation strategy was identified to be associated with limited-maximum transpiration rates. The transpiration of PI 567731 was found to be sensitive to an aquaporin inhibitor (silver-nitrate) indicating the independence of a limited-maximum transpiration to a lack of silver-sensitive aquaporins in these SW genotypes. In efforts to understand many findings from field trials and further mapping of QTLs, many researchers have adopted use of proteomic and metabolomic techniques and platforms for data analysis to further reinforce and probe the mechanisms of plant stress responses. As responses are characteristic to either acute or prolonged effects to drought stress, the initial impacts primarily effect net photosynthesis and photosynthetic performance of the plants. Under drought, stomatal closures and hormonal signaling through abscisic acid have been identified as the key reductants to net photosynthesis, with increased efficiency of the photosystem (PS) II denoted in stress tolerance. Even though water deficiencies do not directly impact the primary components of C3 plants PSI or PSII directly, these secondary impacts are well known in a variety of crops to reversibly impact photosynthesis, prior to photosynthetic decay. This emphasizes the need for targeted approaches of profiling phytochemicals as a reliable means of screening for stress tolerances.

Targeted and non-targeted approaches for determining metabolic profiles of agronomical crops have been entailed with instrumental approaches ranging from gas or liquid chromatography (GC/LC) coupled with mass spectrometry (MS), nuclear magnetic resonance (NMR) and a variety of spectroscopic techniques. No all-inclusive method for simultaneous detection of all metabolites is available. With a broad array of expression in a variety of primary and secondary metabolites in model plants and agricultural crops, methods either prove to be either moderate throughput with high specificity in extracts, or lack specificity with high-throughput analysis. High-resolution accurate mass MS platforms allows for the determination of molecular formulas; when combined with tandem mass spectrometry (MS/MS) molecular structures can be proposed.

Research supported by the sponsors in the previous grant period led to a number of important metabolomics findings for the drought-susceptible cultivar Pana and the drought-tolerant cultivar PI 567731. The plant introduction (PI) PI 567731, which demonstrates a slow wilting canopy phenotype in maturity group III, was profiled in drought stress field trials against a drought susceptible check cultivar, Pana. Relative phytochemical content of chlorophyll (chl) a/b, and pheophytin (pheo) was profiled by direct infusion electrospray Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry. High-throughput detection of metabolic profiles in twenty-four experimental groups occurred in triplicate within a few hours, without chromatographic separation. Multivariate analysis was able to form predictive models, encompassing the variance of growth and drought stress, within the experimental groupings at two physiological ages. Statistically significant increases within the Chl content in control conditions were detected, and an expanded photosynthetic antenna within the drought affected treatment condition could account for increased photosynthetic content; in particular, the distinct enhancement of chl b is noted from PI 567731. Moreover, the existence of unique chl-related metabolites (m/z >900) were confirmed through tandem mass spectrometry. The chlorophyll-related metabolite at m/z 1073.71 is a metabolic product of chlorophyll a with an additional C12H20O- group.

We completed the metabolite libraries of the methanolic leaf extracts of Pana and PI 567731 by filtering the high mass resolution data obtained from each through the SoyCyc and Human Metabolome Database and performing Kendrick mass defect analysis (KMD) to identify ionic formulas. Direct infusion electrospray ionization FT-ICR mass spectra for both cultivars is shown in Figure 1 of the attachment.

A total of 60 ionic formulas certified to be present in soybeans are shared by both Pana and PI 567731 from the over 460 found for Pana and 350 for PI 567731. Prominent amongst these are mono- and diacylglycerols, pheophytin a and chlorophyll a, monosaccharides, disaccharides, xanthins, and vicenin-2 (a flavonoid diglucosylation product). Notable also is the simultaneous presence of plastoquinone, detected with products echinone and plastoquinol, essential components of photosynthetic electron transfer. Likewise, ubiquinol-8 and -9 are detected along with 3-demethylubiquinol-9 and demethylmenaquinol-8, key components of aerobic respiration and photosynthethic electron transfer. The metabolite cycloeucalenone is involved in phytosterol biosynethesis.

In addition, 23 species are unique to Pana in controls, and are shown in Table 1 of the attachmebt, while five are unique to PI 567731 in controls, and are shown in Table 2 of the attachment.

Pana is a drought-sensitive soybean cultivar. Using the known soybean metabolites identified in Table 1, there are several carboxylic acid molecules present in Pana that were not detected in PI 567731; these are essential precursors to lipids. Carlactone is an oxidation product of cartenal, possibly indicating oxidative stress in Pana even in the control which has not experienced drought. This is further supported by the presence of glutathione disulfide, the oxidized dimer of glutathione. Galactopinitols are required substrates and products of galactosylcyclitol biosynthesis. The compound 15-cis-phytoene is needed for production of plastoquinol and carotenes. Likewise, the substance menoquinol-8 is a polyprenyl quinone required for electron transport. A richer complement of pheophytins and chlorophylls are detected in Pana in comparison to PI 567731 (e.g., chlorophyll b was only detected in Pana). However, our earlier work showed that PI 567731 maintains greater levels of pheophytins and chlorophylls during drought.

In contrast, PI 567731 is a drought-tolerant soybean cultivar. The metabolites uniquely detected in the methanolic extract of PI 567731 are shown in Table 2. The galactosyl glycerol compound is 3-ß-D-galactosyl-sn-glycerol, formed from the degradation of diacyl glycerols. Soyasapogenol B is a key precursor in the formation of its glucuronide. There are many possible structures for the trisaccharides, so anabolism of more complex saccharides from mono- and disaccharides might explain the appearance of trisaccharides here. Plastoquinones are electron carriers that are necessary building blocks for plastoquinol, and are found in chloroplasts, thus playing a central role in the photosynthetic electron transport chain. Therefore, PI 567731 may adapt better to drought conditions because of how it processes sugar molecules and builds a reservoir of electron transport carriers.

Having established a methodology to catalog the leaf metabolites, the next phase of this research is to apply high resolution mass spectrometry to identify the novel molecular structures of the metabolites detected in soybean seeds from the drought susceptible cultivar Pana and drought tolerant cultivar PI 567731. In order to accomplish this, these metabolites will be separated from extracts using liquid chromatography (LC) and further be characterized using tandem mass spectrometry (MS/MS). LC is essential here, as noted in Tables 1 and 2, some of the unique ionic formulas detected in methanolic leaf extracts have multiple possible isomeric structures, and indeed a multiplicity of these structures may be present! A previous study of different methanolic, acetone, and methylene chloride extracts of soybean seeds using LC and MS/MS showed that key metabolites such as phytochemicals, flavones and isoflavones, and soyasaponins could be distinguished using this method. Interestingly, the authors found that use of low polarity solvents such as methylene chloride and acetone did not impact the overall species detected in the seed extracts but did lower their abundances; thus, we will exclusively utilize methanol as the extract solvent for the Pana and PI 567731 seeds.
Based on the metabolite identities established by LC-MS/MS from the methanolic extracts of the Pana and PI 567731 seeds, we will compare the known identifications for correlations linking these metabolites to various metabolic pathways; in essence, we are using molecular profiling as an approach to distinguish soybean phenotypes, particularly for the drought-tolerant cultivar so that plant biologists could develop soybean cultivars that possess enhanced drought tolerance. We will compare and contrast the molecular composition of the seeds of the drought tolerant and drought susceptible cultivars. An emphasis in the analysis of this data is to determine whether there may be key nutritional differences between the seeds.


Project Objectives

There are two specific research objectives in this proposal:
1. Having previously built leaf metabolite databases for two cultivars of soybeans (Glycine max), Pana (drought susceptible) and PI 567731 (drought tolerant), perform a similar analysis for the soybean seeds themselves, now using liquid chromatography in combination with high resolution tandem mass spectrometry (MS/MS).
2. Having built seed metabolite databases for the two cultivars, compare and contrast these data to determine how significantly different their molecular compositions are. In particular, a focus of this data is to determine whether there may be key nutritional differences between the seeds.

Project Deliverables

1. A list of metabolites identified in the extracts of seeds of the drought-susceptible cultivar Pana.
2. A list of metabolites identified in the extracts of seeds of the drought-tolerant cultivar PI 567731.
3. Compare the metabolic maps of each cultivar to determine differences in molecular composition between the seeds, with key emphasis on those species which provide nutritional value to humans and livestock.

Progress Of Work

Updated July 30, 2023:
As funding was just recently established, this report will largely detail the plan over the next quarter. As this grant is investigating the nutritional molecular composition of soybean seeds. In our study two cultivars will be compared, Pana (drought susceptible) and PI 567731 (drought tolerant). There are two goals. The first goal is to build a metabolite database for each seed type using the combination of liquid chromatography and tandem mass spectrometry (MS/MS). Much of the development of the method will be done using low resolution mass spectrometry before transfer to the high resolution platform. The second goal is to compare and contrast the metabolite databases of these two cultivars to determine how significantly different their molecular compositions are. In particular, a focus of this data is to determine whether there may be key nutritional differences between the seeds.

Work has begun on milling seeds from each group using mortar and pestle. Once powders of each cultivar have been produced, liquid extraction must be performed. However, no singular extraction method is ideal to obtain the known rich diversity in molecular composition of soybean seeds. Each cultivar powder will be divided into two subsets. Subset A will undergo liquid extraction using ethyl acetate to extract nonpolar analytes from the seeds; this subset should be rich in lipids, a major nutrient in soybeans. Subset B will be extracted by a methanol/water/acetic acid to obtain the polar analytes, including amino acids, soyasaponins, and oligosaccharides. Of these, amino acids and soyasaponins will be analyzed using reversed-phase liquid chromatography, while oligosaccharides will be analyzed using hydrophobic interaction liquid chromatography (HILIC).

Updated October 30, 2023:
This grant is investigating the nutritional molecular composition of soybean seeds, and we are comparing two different cultivars, Pana (drought susceptible) and PI 567731 (drought tolerant). There are two goals. The first goal is to build a metabolite database for each seed type using the combination of liquid chromatography and tandem mass spectrometry (MS/MS). Much of the development of the method will be done using low resolution mass spectrometry before transfer to the high resolution platform. The second goal is to compare and contrast the metabolite databases of these two cultivars to determine how significantly different their molecular compositions are. In particular, this data focuses on key nutritional differences between the seeds.

Due to problems on our high-resolution LC-MS/MS system, identification of products will not be made until this final quarter. However, work has progressed with a lower resolution LC-MS/MS system during the current reporting period. Two different approaches have been developed to obtain powder of soybean seeds. The first approach uses a hard granite mortar and pestle to mill the seeds into a fine powder; this approach works fairly well except for the seed coat, which tends not be ground well with this approach. A far superior approach has been to use a household pepper mill to grind the seeds, which also works for the seed coat as well.

To conduct a metabolic analysis of the seed powder, it is divided into two subsets. The first subset is targeted toward recovery and detection of nonpolar analytes. For this the seed powder is extracted with ethyl acetate. A comparison of this extract to USP soybean oil by low-resolution LC-MS/MS shows numerous matches to lipids and fatty acids known in the USP soybean oil; high-resolution LC-MS/MS will be used to confirm these matches, as well as attempts to identify compounds which are found in the nonpolar fraction but do not match those in USP soybean oil standard.

The second subset is targeted toward polar analytes. Here, the seed powder is extracted using a mixture of methanol/water/acetic acid. To date, reversed-phase liquid chromatography has been used for the analysis of this extract. Numerous of the essential amino acids have been identified in the powders using the low-resolution LC-MS/MS system, which need to be confirmed by high-resolution LC-MS/MS. Later eluting compounds which may indeed be from the soyasaponin family are detectable, but will require the high-resolution LC-MS/MS to obtain formulas of fragment ions needed to assess identity. We are also interested in the composition of oligosaccharides in the polar solvent extracts of the powders. Unsurprisingly, our reversed-phase LC approach is unable to resolve oligosaccharides, so during the current quarter we will be testing a method known as hydrophobic interaction liquid chromatography (HILIC). HILIC is sometimes promising where reversed-phase LC fails because compounds like oligosaccharides, which elute with the void volume in reversed-phase LC, often have long elution times by HILIC. The reverse is also true. This is because the stationary phase in HILIC is also hydrophilic, which means polar analytes have stronger interactions with the stationary phase and hence longer eluting times.

Once the high-resolution LC-MS/MS is back online and identities of metabolites are finalized, we will begin a comparison between the metabolites detected in Pana and PI 567731 to determine nutritional differences between the seeds.

Final Project Results

Updated January 31, 2024:

View uploaded report PDF file

This grant investigated the nutritional molecular composition of soybean seeds. In our study the drought susceptible cultivar Pana (drought susceptible) was investigated. There were two goals. The first goal was to build a metabolite database for each seed type using the combination of liquid chromatography and tandem mass spectrometry (MS/MS). The second goal was to compare and contrast the metabolite databases of the Pana cultivar with the PI 567731 (drought tolerant) cultivar to determine how different their molecular compositions are. While ultimately we were unable to develop a working LC-MS/MS method to separate the components extracted from the seeds, and therefore unable to achieve the second goal, we were able to glean significant information on the molecular composition of Pana seeds even in the absence of chromatography. Whether analyzing soybean powders extracted from solvents compatible with water (hydrophilic or polar) or those with little affinity for water (hydrophobic or nonpolar), from 160 to a little over 200 different peaks attributed to different molecular formulas were detected. The most abundant molecule, detected in both positive and negative ion modes, corresponds to the formula C7H16N2O4. Amino acids, small peptides, and small nucleic acids were detected in the polar and nonpolar powders of soybean seeds. Reversed-phase liquid chromatography, a technique used to separate components, was unsuccessful, suggesting the need for implementation of more polar column chemistries for separation of these compounds.

Benefit To Soybean Farmers

Drought can significantly decrease soybean (Glycine max) yields, and current climate models suggest that drought will become increasingly common in all regions of the United States, including the Northeast. The development of drought-resistant soybean cultivars depends on an improved understanding of stress-response mechanisms. While it is understood that drought stress impacts photosynthesis in soybeans, how this is manifested in the actual metabolic products observed in the plant is largely unknown. Previous support from the New York Corn & Soybean Growers Association allowed us to compile catalogs of two plant introductions (PIs) of soybeans, one which is drought susceptible, Pana (PI 597387), and one which is drought tolerant, (PI 567731). Catalogs of leaf metabolites for both cultivars have been generated, both for young leaves (1 week) and old (2-3 weeks) and abundances of these metabolites have been obtained. Leaf extracts from the two cultivars grown under control conditions share 60 ionic formulas which are matched to the SoyCyc database. Prominent amongst these are mono- and diacylglycerols, pheophytin a, chlorophyll a, monosaccharides, disaccharides, xanthins, and vicenin-2 (a flavonoid diglucosylation product). Notable also is the simultaneous presence of plastoquinone, detected with products echinone and plastoquinol, essential components of photosynthetic electron transfer. Likewise, ubiquinol-8 and -9 are detected along with 3-demethylubiquinol-9 and demethylmenaquinol-8, key components of aerobic respiration and photosynthetic electron transfer. The metabolite cycloeucalenone is involved in phytosterol biosynthesis.

In the drought-susceptible cultivar Pana, drought had a substantial effect on a large number of metabolites within the young plant leaf extracts, almost exclusively expressed at higher level in leaves exposed to drought. Of those that showed 2x or more abundance in drought-treated leaves include monosaccharides, 1-18:3,-2-18:3-monogalactosyldiacylglycerol, pheophytin a, chlorophyll a and a number of chlorophyll-related metabolites. It appears in the Pana cultivar, as water starts to become scarce, most metabolites, including the chlorophylls, are elevated in early stages of drought. However, to our surprise, drought had a relatively small effect on most metabolites as it progresses; in the old Pana leaf extracts, the only identified metabolites which changed in relative abundance by more than a factor of 2, were metabolites with m/z 349.178, 367.188, 635.439, 787.520, and 813.491, with the first four decreasing upon drought. After filtering these m/z values through databases and subsequent Kendrick mass defect (KMD) analysis, the first two metabolites having decreased upon drought were identified as C18H30O4 and C18H32O5. However, the species at m/z 813.491 is glycinolprenol 11, which increased in abundance (2.6x) in the old leaves upon drought-treatment. This metabolite is involved in biosynthesis of polyisoprenoids and suggests a primary mechanism by which the Pana cultivar uses to adapt to drought conditions is to increase its production of isoprenoid lipids, possibly to minimize loss of water by producing more wax on the leaf surface.

For PI 567731, the metabolomics data reveals a far more complex situation. Drought had a significant impact on many metabolites within the young plant leaf extracts; a 2x greater increase in abundance was observed for over 40 metabolites. In addition, a few metabolites are uniquely detected in the methanolic extract of PI 567731 control vs. Pana. One is 3-ß-D-galactosyl-sn-glycerol, formed from the degradation of diacyl glycerols. Soyasapogenol B is a key precursor in the formation of its glucuronide. There are many possible structures for the trisaccharides, so anabolism of more complex saccharides from mono- and disaccharides might explain the appearance of trisaccharides here. Plastoquinones are electron carriers that are necessary building blocks for plastoquinol, and are found in chloroplasts, thus playing a central role in the photosynthetic electron transport chain. Therefore, PI 567731 may adapt better to drought conditions because of how it processes sugar molecules and builds a reservoir of electron transport carriers. In summary, it appears that in the PI 567731 cultivar, much like the Pana cultivar, as water starts to become scarce, most metabolites, including the chlorophylls, are elevated in early stages of drought. Unlike Pana, for the PI 567731 cultivar, drought had a substantial effect on many metabolites within the old plant leaf extracts as well. Overall, the drought-tolerant cultivar, PI 567731, sees a greater change in metabolite levels as drought progresses, so PI 567731 does not adapt to drought in the same way as Pana; indeed, chlorophyll and its related metabolites appear to shift in favor of the metabolite found at m/z 1073.71, which is related to pheophytin a with an additional C12H20O- group (C67H94N4O6).

The United Soybean Research Retention policy will display final reports with the project once completed but working files will be purged after three years. And financial information after seven years. All pertinent information is in the final report or if you want more information, please contact the project lead at your state soybean organization or principal investigator listed on the project.