Next Article in Journal
Convergent Adaptation of Multiple Herbicide Resistance to Auxin Mimics and ALS- and EPSPS-Inhibitors in Brassica rapa from North and South America
Next Article in Special Issue
Effects of Individual or Combined Deficiency of Phosphorous and Zinc on Phenotypic, Nutrient Uptake, and Molecular Responses of Finger Millet (Eleusine coracana): A Nutri-Rich Cereal Crop
Previous Article in Journal
Transcriptional and Post-Translational Regulation of Plant bHLH Transcription Factors during the Response to Environmental Stresses
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 12
Issue 11
10.3390/plants12112117
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soybean Root Transcriptomics: Insights into Sucrose Signaling at the Crossroads of Nutrient Deficiency and Biotic Stress Responses

by
Leela Chandra Manozna Nidumolu
,
Kristina Mae Lorilla
,
Indrani Chakravarty
and
Claudia Uhde-Stone
*
Department of Biological Sciences, California State University, East Bay, Hayward, CA 94542, USA
*
Author to whom correspondence should be addressed.
Plants 2023, 12(11), 2117; https://doi.org/10.3390/plants12112117
Submission received: 11 April 2023 / Revised: 22 May 2023 / Accepted: 22 May 2023 / Published: 26 May 2023
(This article belongs to the Special Issue Nutrient Signaling Networks in Plants)
Browse Figures
Versions Notes

Abstract

:
Soybean (Glycine max) is an important agricultural crop, but nutrient deficiencies frequently limit soybean production. While research has advanced our understanding of plant responses to long-term nutrient deficiencies, less is known about the signaling pathways and immediate responses to certain nutrient deficiencies, such as Pi and Fe deficiencies. Recent studies have shown that sucrose acts as a long-distance signal that is sent in increased concentrations from the shoot to the root in response to various nutrient deficiencies. Here, we mimicked nutrient deficiency-induced sucrose signaling by adding sucrose directly to the roots. To unravel transcriptomic responses to sucrose acting as a signal, we performed Illumina RNA-sequencing of soybean roots treated with sucrose for 20 min and 40 min, compared to non-sucrose-treated controls. We obtained a total of 260 million paired-end reads, mapping to 61,675 soybean genes, some of which are novel (not yet annotated) transcripts. Of these, 358 genes were upregulated after 20 min, and 2416 were upregulated after 40 min of sucrose exposure. GO (gene ontology) analysis revealed a high proportion of sucrose-induced genes involved in signal transduction, particularly hormone, ROS (reactive oxygen species), and calcium signaling, in addition to regulation of transcription. In addition, GO enrichment analysis indicates that sucrose triggers crosstalk between biotic and abiotic stress responses.

Graphical Abstract

1. Introduction

Soybean (Glycine max) is an important food crop in the US and worldwide, and preserving soybean crop yields is critical for the global economy and food security. Nutrient deficiencies, particularly phosphate (Pi) and iron (Fe) deficiencies, frequently limit crop production [1,2,3]. Iron deficiency chlorosis in soybean alone results in global yield losses of millions of metric tons per year [4,5,6].
Plants need to carefully coordinate growth and development with—often scarce—nutrient availability. Glucose has been identified and thoroughly studied as an important signaling molecule of nutrient deficiency and acts by mechanisms similar to those found in yeast and mammals. Much less is known about the mechanism of sucrose signaling, which appears to be unique to plants [7]. A dual role of sucrose as both a metabolite and a signaling molecule was first suggested many decades ago [8] and has gained momentum in recent years [7,9,10,11,12].
A role for sucrose specifically in nutrient deficiency signaling has been proposed during the last two decades [13,14,15]. Split-root experiments in Arabidopsis [16] and white lupin [17,18] revealed that Pi deficiency responses in the root are triggered, at least in part, by low Pi levels of the shoot. This finding initiated a search for candidates for long-distance signals by which the shoot communicates Pi deficiency to the root. Independent studies identified sucrose, transported in the phloem from shoot to root, as an important long-distance signal of Pi deficiency [14,19]. Further studies implicated sucrose in the signaling of several other nutrient deficiencies, including iron (Fe) [15] and sulfur (S) [20]. Studies of an Arabidopsis mutant that is over-accumulating sucrose (hsp1) indicate that sucrose can act as a global regulator of plant responses to multiple nutrient deficiencies, including deficiencies of Pi, nitrogen (N), and potassium (K) [21]. Plant responses to various nutrient deficiencies tend to overlap, a phenomenon known as crosstalk [22], which may be in part due to a shared sucrose signaling pathway.
It is pertinent to point out that the function of sucrose as an initial signal of various nutrient deficiencies is different from its re-allocation as a carbon source to the root. The latter is an acclimation to Pi and N deficiency that facilitates an increase in the root-to-shoot ratio under these conditions [23] and is not the focus of this study.
While it has been shown that sucrose can regulate the expression of certain genes via the activation of transcription factors such as AtMYB75 and AtWRKY20 in Arabidopsis [7,24], it is still unclear how cells receive and process the sucrose signal. Sucrose non-fermenting-1-related protein kinases (SnRKs) have been proposed to serve both as sugar transporters and sugar sensors [25], though further studies are needed to fully elucidate the potential signaling role of sugar transporters.
Significant research has been devoted to understanding plant responses to nutrient deficiencies. In recent years, transcriptomic approaches, particularly RNA-seq, have made great strides toward the identification of differential gene expression in response to various nutrient deficiencies in plants [26,27,28,29,30,31,32]. However, the underlying mechanism of sucrose signaling is not yet well understood [7,11].
In the current study, RNA-seq analysis has been applied to hydroponically grown soybeans treated with short-term (20 min, 40 min) additions of external sucrose to the roots. Our objective was to dissect the early transcriptomic responses to the signal sucrose in roots to identify potential key players in the sucrose signaling network, specifically as part of nutrient deficiency responses. A better understanding of sucrose signaling in response to various nutrient stresses could aid in developing soybeans with increased stress tolerance, productivity, and more efficient use of fertilizers.

2. Results

2.1. RNA-Seq to Assess Short-Term Responses to Sucrose Resulted in 260 Million Paired-End Reads

To mimic the increased amount of sucrose that is transported from shoot to root in response to various nutrient deficiencies, we added sucrose directly to the roots in a hydroponic system. Previous results in the legume white lupin have shown that cluster root formation, a morphological response to Pi and Fe deficiency, can be mimicked by external sucrose application [33]. Internal sucrose concentrations at the initiation zone of cluster root formation were measured at 3.4 mM sucrose [33]. External sucrose application increased the number of cluster roots in a concentration range of 2.5 mM to 12.5 mM sucrose; a further increase to 25 mM sucrose resulted in unusual root thickening [33]. Based on these data, we decided to add an external sucrose concentration of 10 mM, which is high enough to expect a strong response to the sucrose signal but not too high to expect non-physiological responses. After hydroponic growth for four weeks in nutrient-sufficient conditions, roots were subjected to external sucrose for 0 min (control), 20 min, and 40 min in three biological replications each, for a total of nine samples.
The nine pooled, bar-coded cDNA libraries were sequenced using the Illumina NovaSeq6000 150PE Flow Cell SP platform, resulting in an initial 485.64 million total reads of 150 bp each. FastQC quality checks revealed overall very good sequence quality, with Phred quality scores above 30 even before trimming, except for the first two and last three nucleotides. A quality check after trimming showed Phred quality scores above 30 for all positions, complete adaptor removal, and a total number of paired sequences of 260 million (Table 1).
We used HiSat2 [34] and StringTie [35] to map our paired sequence reads to the G. max reference genome [36], which has been reannotated in 2021 (RefSeq assembly accession: GCF_000004515.6_Glycine_max_v4.0; annotation date 3 October 2021). Our mapping resulted in a total of 62,048 mapped genes, which is more than the currently annotated 59,046 genes in GCF_000004515.6_Glycine_max_v4.0, according to the corresponding gene and feature statistics for this annotation (release 104, https://www.ncbi.nlm.nih.gov/genome/annotation_euk/Glycine_max/104/ (accessed on 14 October 2022)).
Some of these novel genes were found at a very low level, raising the question if they are real. When setting a threshold of at least 3 total reads, we found 61,675 genes, which still includes some novel genes. Applying the same threshold of at least 3 total reads, StringTie identified 99,775 different transcripts, indicating alternative splicing.

2.2. Short-Term Sucrose Exposure Changes the Expression of Hundreds of Genes

MA (mean average) plots (Figure 1a) were used to visualize log2 FC (fold change) against normalized sequence counts at 20 min and 40 min of sucrose exposure, each compared to 0 min (control). The plot reveals some changes in expression (blue dots) at 20 min and more drastic changes after 40 min of sucrose exposure. A PCA (principal component analysis) plot (Figure 1b) displays clear differences between the three time points used in this study (t0, t20, and t40), though it also reveals some variability between the biological replicates.
We extracted normalized expression data (FPKM, fragment per kilobase, and million) from our HiSat2 data using Ballgown [35,37]. A Venn diagram (Figure 1c) represents the total number of up- and down-regulated genes in response to 20 min and 40 min of sucrose treatment. A total of 148 upregulated and 58 downregulated genes are shared between 20 min and 40 min of sucrose responses.

2.3. Genes Involved in Auxin-Mediated Responses and Nutrient Assimilation Are among the Most Highly Upregulated Genes in Response to Sucrose

A heatmap (Figure 2) gives an overview of normalized expression (in FPKM) for the top 40 up-regulated genes after 20 and 40 min of sucrose treatment, revealing low expression levels without sucrose and increased expression at both 20 and 40 min of sucrose.
Table 2 shows the 10 most highly upregulated genes in response to 20 min of sucrose. The description column indicates the involvement of upregulated genes in various functions, including auxin signaling, nitrogen and sulfur assimilation, and transcriptional regulation. A heatmap of the 20 most downregulated genes in response to sucrose is shown in Figure S1.

2.4. Gene Ontology Analysis Reveals High Proportions of Transcription Factors and Internal Signaling in Response to Sucrose

An analysis of GO (gene ontology) terms indicates a high proportion of sucrose-induced genes belong to the categories of signal transduction and stress response (Figure 3a). The largest two groups among molecular functions are protein binding and transcription factors (Figure 3b). Because our GO analysis of molecular functions revealed a high proportion of transcription factors (TFs), we took a closer look at this group of genes. Our RNA-seq analysis revealed a total of 143 transcription factors (TFs) that were differentially expressed (log2FC > 1.5, p-value < 0.5 or <−1.5, p-value < 0.5) in response to sucrose; the majority of these (128) were upregulated, while 15 were downregulated.
At 20 min of sucrose treatment, only nine TFs were significantly upregulated, and only two of these (WRKY 70-like and UNE12) were also upregulated after 40 min of sucrose. In contrast, a total of 121 TFs were upregulated at 40 min of sucrose treatment. In terms of frequency, this means that at t20, 2.5% of upregulated genes were transcription factors, while at t40, this number increased to 5%. Table 3 shows the ten most highly upregulated transcription factors in response to 40 min of sucrose. While only t40 shows significant p-values, the t20 fold changes indicate that the trend of upregulation of these transcription factors starts earlier. The most frequent families among upregulated TFs are WRKY (33 TFs), ethylene response factors (ERFs, 24 TFs), MYB (22 TFs), and basic helix-loop-helix (bLHL, 14 TFs).

2.5. GO Enrichment Indicates Involvement of ROS and Ca2+ Signaling in Responses to Sucrose

Enrichment analysis of GO terms displayed no significantly enriched categories at t20 but highly enriched categories at t40. Table 4 shows the fifteen most highly enriched categories at 40 min of sucrose exposure. The three most highly enriched categories.
The response to chitin, respiratory burst, and hypersensitive response all indicate a strong biotic stress response induced by 40 min of sucrose. Interestingly, ROS production (as in respiratory burst) can also be involved in abiotic stress responses [47].
To further analyze the possible involvement of ROS signaling in response to sucrose, we took a closer look at class III peroxidases (PERs), which can be involved in signaling pathways in response to both biotic and abiotic stresses [48]. We found no PERs significantly induced at 20 min of sucrose treatment. However, at 40 min of sucrose, 18 PERs were significantly upregulated; the ten most highly upregulated peroxidase genes are shown in Table S1. ROS signaling, in turn, can involve calcium-binding proteins and kinases. GO terms related to calcium signaling, such as Ca2+ binding, were enriched among genes upregulated after 40 min of sucrose treatment (corrected p-value 0.0001). We identified 45 genes putatively involved in Ca2+ signaling that were upregulated at t40, six of which were already upregulated at t20; the ten most highly upregulated genes in this category are shown in Table S2.

3. Discussion

3.1. RNA-Seq to Unravel the Regulatory Network in Response to Sucrose

Sucrose can act both as a signal and as a carbon source. Sucrose, transported via phloem from the shoot to the roots, has been identified as a long-distance signal for various nutrient deficiencies, including Pi [13,14], Fe [15], sulfur [20], and nitrogen [21]. It is worth pointing out that sucrose signaling is different from sucrose being reallocated as a carbon source, which is found only under certain nutrient deficiencies (Pi and N) [23] and is not addressed in our study.
RNA-seq proved useful as a high-throughput method to look at the global effect of sucrose on gene expression in order to unravel the regulatory network that is activated by the increase of sucrose in roots. Our results reveal a complex regulatory network, including a high proportion of transcription factors and members of signal transduction pathways, which we will discuss below in more detail.
As we are particularly interested in key regulators, we analyzed our data for differentially expressed transcription factors. Of the transcription factors upregulated in this study, WRKY, ERFs (ethylene response factor), and MYB were the three largest groups. Members of these families of transcription factors are known to play critical roles in regulating plant growth, development, and stress responses. Several studies have shown that members of these families of transcription factors are involved in mediating sucrose signaling in plants. AtWRKY20 and AtMYB75 have been identified as upstream regulators controlling sucrose-responsive genes (reviewed by [7]). Our study identified several transcription factors induced after 20 min and more than hundred transcription factors induced after 40 min of sucrose, providing insights into the regulatory networks that govern plant responses to sucrose signaling. We found that the most frequent families among upregulated TFs are WRKY, ethylene response factors, MYB14, MYB30, WRKY41, and WRKY33, and basic helix-loop-helix, with MYB and WRKY representing the most highly upregulated transcription factors.

3.2. ROS and Ca2+ Signaling Act Downstream of Sucrose

GO (gene ontology) analysis revealed that genes involved in the production of reactive oxygen species and intracellular signaling were highly enriched categories. Several class III peroxidases, calcium-dependent protein kinases (CDPKs), and calcium-binding proteins were upregulated after 20 min, and many more after 40 min of sucrose treatment. ROS signaling is a known response to biotic and abiotic stresses, such as salinity, heavy metals, and nutrient deficiency [47,49]. ROS can be generated by various enzymatic pathways, including NADPH oxidases and class III peroxidases, which can act as ROS producers or consumers [50]. ROS functions in concert with other important second messengers, specifically Ca2+, which in turn is part of Ca2+ signaling via activation of calcium-dependent protein kinases (CDPKs) and calcium-binding proteins, such as calmodulins [51,52]. The fact that some class III peroxidases and proteins involved in Ca2+ signaling were upregulated after 20 min and many more after 40 min of sucrose treatment suggests that both ROS (reactive oxygen species) and Ca2+ signaling pathways are triggered by sucrose and thus act downstream of sucrose signaling.

3.3. Sucrose Activates Plant Hormone Signaling

The two most highly upregulated genes at 20 min of sucrose exposure were auxin response factor 9 (ARF9, an auxin-induced transcriptional regulator) and aminopeptidase M1, which has been shown to be involved in polar auxin transport [39]. Auxin is known as an important mediator acting downstream of sucrose in response to multiple nutrient deficiencies. For example, increased sucrose accumulation has been shown to regulate Fe-deficiency responses in Arabidopsis by promoting auxin signaling [15].
Salicylic acid and jasmonic acid signaling were among the significantly enriched categories at 40 min of sucrose. Both hormones are known to play important roles in ROS signaling, as they are involved in the regulation of ROS production and scavenging in response to biotic and abiotic stresses [53]. These data indicate that plant hormone signaling, particularly auxin, jasmonic acid, and salicylic acid, is involved in the immediate plant responses to sucrose.

3.4. Sucrose May Mediate Crosstalk between Biotic and Abiotic Stresses

Crosstalk is the somewhat surprising phenomenon of plants not only responding to a specific stress but also to other stresses they do not experience at this point [54,55,56]. Because sucrose is a signal of various abiotic stresses, including multiple nutrient deficiencies [21], we speculate that sucrose may be responsible for the reported cross-talk between plant responses to various nutrient deficiencies. Supporting this hypothesis, we found several genes involved in nitrogen and sulfur assimilation among our top upregulated genes in response to 20 min of sucrose. In addition, our findings reveal that sucrose may also mediate crosstalk between biotic and abiotic stresses, as genes involved in both stress types were enriched among sucrose-responsive genes.
Taken together, our results provide insights into transcriptional regulators, ROS- and Ca2+-dependent signaling, as well as hormonal responses, that are all activated by sucrose. In the future, it will be interesting to look at even earlier time points to find the earliest key regulators responding to sucrose signaling. In addition, quantitative phosphoproteomics could be used to identify members of the signal transduction pathways activated in response to sucrose, as some of these are differentially phosphorylated rather than differentially expressed.

4. Materials and Methods

4.1. Seed Germination, Treatments, and Harvest

Glycine max cv. Williams 82 seeds (friendly gift from John Harada and Julie Marie Pelletier, UC Davis) were sterilized by shaking for 4 min in 10% bleach, followed by 6 rinses with autoclaved water. Agar plates were prepared with 0.6% agarose in petri dishes, and seeds were placed on the agar plates, covered with aluminum foil, and incubated at 27 °C. After about 4 days, once the radicle reached 3 to 5 cm in length, germinated seeds were transferred to containers containing about 850 mL of half-strength Hoagland nutrient solution [57] and were grown for 4 weeks; the Hoagland solution was changed every 4 days. The growth chamber conditions were maintained at ~27 °C with a light cycle of 16 h [58]. For sucrose treatment, 8.5 mL of 1 M sucrose (prepared in a half-strength Hoagland solution) was added directly to the hydroponic solution for a final concentration of 10 mM sucrose.
Harvesting occurred at 0 min (control) and after 20 min and 40 min of sucrose addition. To enable statistical data analysis, all-time points were performed in 3 biological replications (3 plants). Per plant, about 100 mg of root tip sections (~2 cm) were harvested in liquid nitrogen and immediately stored at −80 °C.

4.2. RNA Isolation and Quality Check

RNA from frozen root samples was isolated using the RNeasy Plant Mini kit (Qiagen), following the protocol for “Purification of Total RNA from Plant Cells and Tissues, and Filamentous Fungi”. The Qubit 4 Fluorometer (Thermofisher Scientific, Waltham, MA, USA) was used in conjunction with RNA-High Sensitivity and RNA IQ assays to assess both the quantity and quality of each sample. To assess quality, the Qubit RNA IQ assay represents the ratio of large and/or structured RNA to small RNA in the sample. A low score indicates that the sample consists mainly of small RNA, while a high score indicates that the sample consists mainly of large and/or structured RNA. An RNA IQ score of 8 or higher represents good quality RNA; we only used samples with an RNA IQ score of 8 or greater for RNA sequencing.

4.3. cDNA Library Preparation and RNA-Sequencing

Using the Stranded mRNA Kit (Illumina, Foster City, CA, USA), we converted the extracted RNA to cDNA, following the manufacturer’s instructions. Unique dual barcoding for each cDNA library was performed using the IDT for Illumina RNA UD Indexes Set A (Illumina).
Quantity and quality of cDNA libraries were assessed via TapeStation (Agilent, Santa Clara, CA, USA), and remaining adaptors were cleaned up where necessary using AMPure XP magnetic beads (New England Biolabs, Ipswich, MA, USA at a ratio of 1 volume of DNA to 0.8 volumes of beads.
Exact quantification, pooling, and sequencing of the nine barcoded cDNA libraries were performed by QB3 (the Institute for Quantitative Biosciences at UC Berkeley, CA, USA). Sequencing of the pooled cDNA libraries was performed on one flowcell lane on the Illumina Nova-Seq6000 150PE (paired-end) next-generation sequencing platform.

4.4. RNA-Seq Data Analysis

Illumina Conversion Software bcl2fastq2 (v2.20) was used to demultiplex the obtained sequencing data and to convert base call files into FASTQ files.
We transferred demultiplexed FastQ files containing 485.64 million reads from the QB3 server to our storage allocation at the EXPANSE supercomputer housed at SDSC (San Diego Supercomputer Center, CA USA). We checked sequence quality using FastQC (https://qubeshub.org/resources/fastqc, accessed on 1 August 2022). Next, we removed adaptors and any low-quality sequences with TRIMMOMATIC Version 0.32 [59]. The quality of the sequences was again checked by FastQC to ensure TRIMMOMATIC properly removed all adaptors and low-quality regions. Then all paired RNA-seq reads were mapped with HiSat2 [34] to the G. max reference genome [36], which has been reannotated in 2021 (RefSeq assembly accession: GCF_000004515.6_Glycine_max_v4.0, annotation date 3 October 2021). After initial mapping, we used StringTie [35] for transcript assembly.

4.5. Differential Expression Analysis

We used the prepDE.py3 script (http://ccb.jhu.edu/software/stringtie/dl/prepDE.py3, accessed on 20 October 2022) to extract two csv (comma separated value) files with transcript and gene count information. These files were explored further in RStudio. The normalized expression data FPKM (fragment per kilobase and million) was extracted using the R package Ballgown, which allows spliced transcriptome assembly for differential expression analysis [37]. To include DEGs (differentially expressed genes) with expression below detection under some conditions (circumventing the problem of not being able to divide by 0), a small FPKM value of 1 was added to all FPKM data points.
To identify DEGs, we set a threshold of log2FC (fold change) ≥ 1.5 or ≤1.5 and a p-value < 0.05. The corresponding LOC gene IDs were mapped to GLYMA IDs using the rentrez package (Entrez in R, v. 1.2.3; http://cran.nexr.com/web/packages/rentrez/ accessed on 2 November 2022), which is an R-based package provided by NCBI. The terms were mapped by querying the LOC IDs in the NCBI gene database, then parsing the GLYMA IDs from the locus tag of the gene from its summary.
DESeq2 was used to generate MA (mean average) and PCA (principal component analysis) plots [60]. Heatmaps of FPKM data were generated by gplots (https://cran.r-project.org/web/packages/gplots/index.html, accessed on 15 November 2022). Enrichment analysis of GO terms was performed on DEGs using the GO term enrichment tool (https://www.soybase.org/tools.php, accessed on 2 March 2023) on SoyBase [61,62].

5. Conclusions

In conclusion, our RNA-seq analysis revealed that short-term exposure of roots to sucrose activates multiple signaling pathways, including ROS, Ca2+, and hormone signaling. As sucrose induces responses to both biotic and abiotic stresses, we further conclude that sucrose acts as a mediator of crosstalk between biotic and abiotic stresses. Understanding the mechanisms of sucrose signaling in plants is important for developing strategies to improve plant productivity under biotic and abiotic stresses.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12112117/s1, Figure S1: Heatmap of the 40 most down-regulated genes after 20 min and 40 min of sucrose treatment (sorted for t20) across biological replications, Table S1: Ten most highly upregulated class III peroxidases after 40 min of sucrose treatment, Table S2: Ten most highly upregulated genes involved in calcium signaling after 40 min of sucrose treatment.

Author Contributions

Conceptualization, L.C.M.N. and C.U.-S.; Data curation, L.C.M.N., K.M.L. and C.U.-S.; Formal analysis, L.C.M.N., K.M.L. and C.U.-S.; Funding acquisition, C.U.-S.; Methodology, L.C.M.N., K.M.L., I.C. and C.U.-S.; Resources, C.U.-S.; Software, L.C.M.N., K.M.L. and C.U.-S.; Supervision, C.U.-S.; Visualization, L.C.M.N. and K.M.L.; Writing—original draft, L.C.M.N.; Writing—review and editing, L.C.M.N., K.M.L., I.C. and C.U.-S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by an internal RSCA grant for 2021–2022 at California State University, East Bay.

Data Availability Statement

The raw data and processed FPKM data for all samples of this RNA-seq experiment are openly accessible through the NCBI GEO Series accession number GSE228888. (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE228888, submitted 3 April 2023).

Acknowledgments

Soybean seeds (Williams 82) for this experiment were a friendly gift from John Harada and Julie Marie Pelletier at UC Davis.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Marschner, H. Mineral Nutrition of Higher Plants, 2nd ed.; Academic Press Inc.: San Diego, CA, USA, 1995. [Google Scholar]
  2. Mo, X.; Liu, G.; Zhang, Z.; Lu, X.; Liang, C.; Tian, J. Mechanisms Underlying Soybean Response to Phosphorus Deficiency through Integration of Omics Analysis. Int. J. Mol. Sci. 2022, 23, 4592. [Google Scholar] [CrossRef] [PubMed]
  3. Vance, C.P.; Uhde-Stone, C.; Allan, D.L. Phosphorus acquisition and use: Critical adaptations by plants for securing a nonrenewable resource. N. Phytol. 2003, 157, 423–447. [Google Scholar] [CrossRef] [PubMed]
  4. Dobbels, A.A.; Lorenz, A.J. Soybean iron deficiency chlorosis high throughput phenotyping using an unmanned aircraft system. Plant Methods 2019, 15, 97. [Google Scholar] [CrossRef] [PubMed]
  5. Merry, R.; Dobbels, A.A.; Sadok, W.; Naeve, S.; Stupar, R.M.; Lorenz, A.J. Iron deficiency in soybean. Crop Sci. 2022, 62, 36–52. [Google Scholar] [CrossRef]
  6. Naeve, S.L. Iron deficiency chlorosis in soybean. Agron. J. 2006, 98, 1575–1581. [Google Scholar] [CrossRef]
  7. Yoon, J.; Cho, L.H.; Tun, W.; Jeon, J.S.; An, G. Sucrose signaling in higher plants. Plant Sci. 2021, 302, 110703. [Google Scholar] [CrossRef]
  8. Pontis, H.G. On the scent of the riddle of sucrose. Trends Biochem. Sci. 1978, 3, 137–139. [Google Scholar] [CrossRef]
  9. Koch, K. Sucrose metabolism: Regulatory mechanisms and pivotal roles in sugar sensing and plant development. Curr. Opin. Plant Biol. 2004, 7, 235–246. [Google Scholar] [CrossRef]
  10. Ruan, Y.L. Signaling role of sucrose metabolism in development. Mol. Plant 2012, 5, 763–765. [Google Scholar] [CrossRef]
  11. Tognetti, J.A.; Pontis, H.G.; Martinez-Noel, G.M. Sucrose signaling in plants: A world yet to be explored. Plant Signal. Behav. 2013, 8, e23316. [Google Scholar] [CrossRef]
  12. Wind, J.; Smeekens, S.; Hanson, J. Sucrose: Metabolite and signaling molecule. Phytochemistry 2010, 71, 1610–1614. [Google Scholar] [CrossRef]
  13. Hammond, J.P.; White, P.J. Sugar signaling in root responses to low phosphorus availability. Plant Physiol. 2011, 156, 1033–1040. [Google Scholar] [CrossRef]
  14. Hammond, J.; White, P. Sucrose transport in the phloem: Integrating root responses to phosphorus starvation. J. Exp. Bot. 2008, 59, 93–109. [Google Scholar] [CrossRef]
  15. Lin, X.Y.; Ye, Y.Q.; Fan, S.K.; Jin, C.W.; Zheng, S.J. Increased Sucrose Accumulation Regulates Iron-Deficiency Responses by Promoting Auxin Signaling in Arabidopsis Plants. Plant Physiol. 2016, 170, 907–920. [Google Scholar] [CrossRef]
  16. Thibaud, M.C.; Arrighi, J.F.; Bayle, V.; Chiarenza, S.; Creff, A.; Bustos, R.; Paz-Ares, J.; Poirier, Y.; Nussaume, L. Dissection of local and systemic transcriptional responses to phosphate starvation in Arabidopsis. Plant J. 2010, 64, 775–789. [Google Scholar] [CrossRef]
  17. Shen, J.; Li, H.; Neumann, G.; Zhang, F. Nutrient uptake, cluster root formation and exudation of protons and citrate in Lupinus albus as affected by localized supply of phosphorus in a split-root system. Plant Sci. 2005, 168, 837–845. [Google Scholar] [CrossRef]
  18. Shane, M.; De Vos, M.; De Roock, S.; Lambers, H. Shoot P status regulates cluster-root growth and citrate exudation in Lupinus albus grown with a divided root system. Plant Cell Environ. 2003, 26, 265–273. [Google Scholar] [CrossRef]
  19. Liu, J.; Vance, C.P. Crucial roles of sucrose and microRNA399 in systemic signaling of P deficiency: A tale of two team players? Plant Signal. Behav. 2010, 5, 1556–1560. [Google Scholar] [CrossRef]
  20. Wawrzynska, A.; Piotrowska, J.; Apodiakou, A.; Bruckner, F.; Hoefgen, R.; Sirko, A. The SLIM1 transcription factor affects sugar signaling during sulfur deficiency in Arabidopsis. J. Exp. Bot. 2022, 73, 7362–7379. [Google Scholar] [CrossRef]
  21. Lei, M.; Liu, D. Sucrose regulates plant responses to deficiencies in multiple nutrients. Plant Signal. Behav. 2011, 6, 1247–1249. [Google Scholar] [CrossRef]
  22. Ruffel, S. Nutrient-related Long-Distance Signals: Common players and possible crosstalk. Plant Cell Physiol. 2018, 59, 1723–1732. [Google Scholar] [CrossRef] [PubMed]
  23. Hermans, C.; Hammond, J.; White, P.; Verbruggen, N. How do plants respond to nutrient shortage by biomass allocation? Trends Plant Sci. 2006, 11, 610–617. [Google Scholar] [CrossRef] [PubMed]
  24. Teng, S.; Keurentjes, J.; Bentsink, L.; Koornneef, M.; Smeekens, S. Sucrose-specific induction of anthocyanin biosynthesis in Arabidopsis requires the MYB75/PAP1 gene. Plant Physiol. 2005, 139, 1840–1852. [Google Scholar] [CrossRef] [PubMed]
  25. Lalonde, S.; Boles, E.; Hellmann, H.; Barker, L.; Patrick, J.W.; Frommer, W.B.; Ward, J.M. The dual function of sugar carriers. Transport and sugar sensing. Plant Cell 1999, 11, 707–726. [Google Scholar] [CrossRef]
  26. Cao, P.; Zhao, Y.; Wu, F.; Xin, D.; Liu, C.; Wu, X.; Lv, J.; Chen, Q.; Qi, Z. Multi-Omics Techniques for Soybean Molecular Breeding. Int. J. Mol. Sci. 2022, 23, 4994. [Google Scholar] [CrossRef]
  27. Kohlhase, D.R.; McCabe, C.E.; Singh, A.K.; O’Rourke, J.A.; Graham, M.A. Comparing Early Transcriptomic Responses of 18 Soybean (Glycine max) Genotypes to Iron Stress. Int. J. Mol. Sci. 2021, 22, 11643. [Google Scholar] [CrossRef]
  28. Mai, H.J.; Pateyron, S.; Bauer, P. Iron homeostasis in Arabidopsis thaliana: Transcriptomic analyses reveal novel FIT-regulated genes, iron deficiency marker genes and functional gene networks. BMC Plant Biol. 2016, 16, 211. [Google Scholar] [CrossRef]
  29. O’Rourke, J.A.; Graham, M.A. Gene Expression Responses to Sequential Nutrient Deficiency Stresses in Soybean. Int. J. Mol. Sci. 2021, 22, 1252. [Google Scholar] [CrossRef]
  30. Moran Lauter, A.N.; Rutter, L.; Cook, D.; O’Rourke, J.A.; Graham, M.A. Examining Short-Term Responses to a Long-Term Problem: RNA-Seq Analyses of Iron Deficiency Chlorosis Tolerant Soybean. Int. J. Mol. Sci. 2020, 21, 3591. [Google Scholar] [CrossRef]
  31. Wang, S.; Sun, S.; Guo, R.; Liao, W.; Shou, H. Transcriptomic Profiling of Fe-Responsive lncRNAs and Their Regulatory Mechanism in Rice. Genes 2021, 12, 567. [Google Scholar] [CrossRef]
  32. O’Rourke, J.A.; McCabe, C.E.; Graham, M.A. Dynamic gene expression changes in response to micronutrient, macronutrient, and multiple stress exposures in soybean. Funct. Integr. Genom. 2020, 20, 321–341. [Google Scholar] [CrossRef]
  33. Wang, Z.; Shen, J.; Ludewig, U.; Neumann, G. A re-assessment of sucrose signaling involved in cluster-root formation and function in phosphate-deficient white lupin (Lupinus albus). Physiol. Plant. 2015, 154, 407–419. [Google Scholar] [CrossRef]
  34. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef]
  35. Pertea, M.; Kim, D.; Pertea, G.M.; Leek, J.T.; Salzberg, S.L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 2016, 11, 1650–1667. [Google Scholar] [CrossRef]
  36. Schmutz, J.; Cannon, S.B.; Schlueter, J.; Ma, J.; Mitros, T.; Nelson, W.; Hyten, D.L.; Song, Q.; Thelen, J.J.; Cheng, J.; et al. Genome sequence of the palaeopolyploid soybean. Nature 2010, 463, 178–183. [Google Scholar] [CrossRef]
  37. Frazee, A.C.; Pertea, G.; Jaffe, A.E.; Langmead, B.; Salzberg, S.L.; Leek, J.T. Ballgown bridges the gap between transcriptome assembly and expression analysis. Nat. Biotechnol. 2015, 33, 243–246. [Google Scholar] [CrossRef]
  38. Bouzroud, S.; Gouiaa, S.; Hu, N.; Bernadac, A.; Mila, I.; Bendaou, N.; Smouni, A.; Bouzayen, M.; Zouine, M. Auxin Response Factors (ARFs) are potential mediators of auxin action in tomato response to biotic and abiotic stress (Solanum lycopersicum). PLoS ONE 2018, 13, e0193517. [Google Scholar] [CrossRef]
  39. Murphy, A.S.; Hoogner, K.R.; Peer, W.A.; Taiz, L. Identification, purification, and molecular cloning of N-1-naphthylphthalmic acid-binding plasma membrane-associated aminopeptidases from Arabidopsis. Plant Physiol. 2002, 128, 935–950. [Google Scholar] [CrossRef]
  40. Lee, O.R.; Cho, H.T. Cytoplasm localization of aminopeptidase M1 and its functional activity in root hair cells and BY-2 cells. Mol. Biol. Rep. 2012, 39, 10211–10218. [Google Scholar] [CrossRef]
  41. Kherraz, K.; Kherraz, K.; Kameli, A. Homology modeling of Ferredoxin-nitrite reductase from Arabidopsis thaliana. Bioinformation 2011, 6, 115–119. [Google Scholar] [CrossRef]
  42. Mazarei, M.; Lennon, K.A.; Puthoff, D.P.; Rodermel, S.R.; Baum, T.J. Expression of an Arabidopsis phosphoglycerate mutase homologue is localized to apical meristems, regulated by hormones, and induced by sedentary plant-parasitic nematodes. Plant Mol. Biol. 2003, 53, 513–530. [Google Scholar] [CrossRef] [PubMed]
  43. Takano, A.; Kakehi, J.; Takahashi, T. Thermospermine is not a minor polyamine in the plant kingdom. Plant Cell Physiol. 2012, 53, 606–616. [Google Scholar] [CrossRef] [PubMed]
  44. Yariuchi, Y.; Okamoto, T.; Noutoshi, Y.; Takahashi, T. Responses of Polyamine-Metabolic Genes to Polyamines and Plant Stress Hormones in Arabidopsis Seedlings. Cells 2021, 10, 3283. [Google Scholar] [CrossRef] [PubMed]
  45. Martinez-Cuenca, M.R.; Iglesias, D.J.; Talon, M.; Abadia, J.; Lopez-Millan, A.F.; Primo-Millo, E.; Legaz, F. Metabolic responses to iron deficiency in roots of Carrizo citrange [Citrus sinensis (L.) Osbeck. x Poncirus trifoliata (L.) Raf]. Tree Physiol. 2013, 33, 320–329. [Google Scholar] [CrossRef] [PubMed]
  46. Tamadaddi, C.; Verma, A.K.; Zambare, V.; Vairagkar, A.; Diwan, D.; Sahi, C. J-like protein family of Arabidopsis thaliana: The enigmatic cousins of J-domain proteins. Plant Cell Rep. 2022, 41, 1343–1355. [Google Scholar] [CrossRef] [PubMed]
  47. Waszczak, C.; Carmody, M.; Kangasjarvi, J. Reactive Oxygen Species in Plant Signaling. Annu. Rev. Plant Biol. 2018, 69, 209–236. [Google Scholar] [CrossRef]
  48. Kidwai, M.; Ahmad, I.Z.; Chakrabarty, D. Class III peroxidase: An indispensable enzyme for biotic/abiotic stress tolerance and a potent candidate for crop improvement. Plant Cell Rep. 2020, 39, 1381–1393. [Google Scholar] [CrossRef]
  49. Fichman, Y.; Mittler, R. Rapid systemic signaling during abiotic and biotic stresses: Is the ROS wave master of all trades? Plant J. Cell Mol. Biol. 2020, 102, 887–896. [Google Scholar] [CrossRef]
  50. Cosio, C.; Dunand, C. Specific functions of individual class III peroxidase genes. J. Exp. Bot. 2009, 60, 391–408. [Google Scholar] [CrossRef]
  51. Drerup, M.M.; Schlucking, K.; Hashimoto, K.; Manishankar, P.; Steinhorst, L.; Kuchitsu, K.; Kudla, J. The Calcineurin B-like calcium sensors CBL1 and CBL9 together with their interacting protein kinase CIPK26 regulate the Arabidopsis NADPH oxidase RBOHF. Mol. Plant 2013, 6, 559–569. [Google Scholar] [CrossRef]
  52. Steinhorst, L.; Kudla, J. Signaling in cells and organisms—Calcium holds the line. Curr. Opin. Plant Biol. 2014, 22, 14–21. [Google Scholar] [CrossRef]
  53. Myers, R.J.; Fichman, Y.; Zandalinas, S.I.; Mittler, R. Jasmonic acid and salicylic acid modulate systemic reactive oxygen species signaling during stress responses. Plant Physiol. 2023, 191, 862–873. [Google Scholar] [CrossRef]
  54. Bouain, N.; Krouk, G.; Lacombe, B.; Rouached, H. Getting to the Root of Plant Mineral Nutrition: Combinatorial Nutrient Stresses Reveal Emergent Properties. Trends Plant Sci. 2019, 24, 542–552. [Google Scholar] [CrossRef]
  55. Mendoza-Cozatl, D.G.; Gokul, A.; Carelse, M.F.; Jobe, T.O.; Long, T.A.; Keyster, M. Keep talking: Crosstalk between iron and sulfur networks fine-tunes growth and development to promote survival under iron limitation. J. Exp. Bot. 2019, 70, 4197–4210. [Google Scholar] [CrossRef]
  56. Courbet, G.; Gallardo, K.; Vigani, G.; Brunel-Muguet, S.; Trouverie, J.; Salon, C.; Ourry, A. Disentangling the complexity and diversity of crosstalk between sulfur and other mineral nutrients in cultivated plants. J. Exp. Bot. 2019, 70, 4183–4196. [Google Scholar] [CrossRef]
  57. Hoagland, D.R.; Arnon, D.I. The water-culture method for growing plants without soil. Circ. Calif. Agric. Exp. Stn. 1950, 347, 32. [Google Scholar]
  58. Leggett, J.E. Growth and nutrient uptake by soybean plants in nutrient solutions of graded concentrations. Plant Physiol. 1971, 48, 457–460. [Google Scholar] [CrossRef]
  59. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef]
  60. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef]
  61. Brown, A.V.; Conners, S.I.; Huang, W.; Wilkey, A.P.; Grant, D.; Weeks, N.T.; Cannon, S.B.; Graham, M.A.; Nelson, R.T. A new decade and new data at SoyBase, the USDA-ARS soybean genetics and genomics database. Nucleic Acids Res. 2021, 49, D1496–D1501. [Google Scholar] [CrossRef]
  62. Grant, D.; Nelson, R.T.; Cannon, S.B.; Shoemaker, R.C. SoyBase, the USDA-ARS soybean genetics and genomics database. Nucleic Acids Res. 2010, 38, D843–D846. [Google Scholar] [CrossRef] [PubMed]
Plants 12 02117 g001 550
Figure 1. (a) MA (mean average) plot of log2 FC against normalized sequence counts at 20 min (t20) and 40 min (t40) of sucrose treatment, each compared to t0 (no-sucrose control). To reduce background noise, shrinkage of the effect size was applied before visualization. Values of padj (adjusted p-value) ≤0.01 in the DESeq2 gene expression analysis are shown in blue. (b) PCA plot of the three biological replications representing differences between both time points (20 min, 40 min) and control (t0). (c) Venn diagram: significantly up-regulated DEGs in t20 and t40 with a p-value < 0.05 and a log2 FC ≥ 1 are represented on top (green), and significantly down-regulated genes in t20 and t40 with a p-value < 0.05 and a log2 FC ≤ 1 are represented on the bottom (blue).
Figure 1. (a) MA (mean average) plot of log2 FC against normalized sequence counts at 20 min (t20) and 40 min (t40) of sucrose treatment, each compared to t0 (no-sucrose control). To reduce background noise, shrinkage of the effect size was applied before visualization. Values of padj (adjusted p-value) ≤0.01 in the DESeq2 gene expression analysis are shown in blue. (b) PCA plot of the three biological replications representing differences between both time points (20 min, 40 min) and control (t0). (c) Venn diagram: significantly up-regulated DEGs in t20 and t40 with a p-value < 0.05 and a log2 FC ≥ 1 are represented on top (green), and significantly down-regulated genes in t20 and t40 with a p-value < 0.05 and a log2 FC ≤ 1 are represented on the bottom (blue).
Plants 12 02117 g001
Plants 12 02117 g002 550
Figure 2. Heatmap of the top 40 up-regulated genes in both t20 and t40 (sorted for t20) across biological replications. Shown are the relative numbers of normalized reads in FPKM for the three biological replications: t0 (control), t20, and t40.
Figure 2. Heatmap of the top 40 up-regulated genes in both t20 and t40 (sorted for t20) across biological replications. Shown are the relative numbers of normalized reads in FPKM for the three biological replications: t0 (control), t20, and t40.
Plants 12 02117 g002
Plants 12 02117 g003a 550Plants 12 02117 g003b 550
Figure 3. GO analysis of biological processes (a) and molecular functions (b) among highly upregulated genes (Log2FC > 3, p-value < 0.05) at t40.
Figure 3. GO analysis of biological processes (a) and molecular functions (b) among highly upregulated genes (Log2FC > 3, p-value < 0.05) at t40.
Plants 12 02117 g003aPlants 12 02117 g003b
Table
Table 1. Number of paired sequences for each of the nine cDNA libraries.
Table 1. Number of paired sequences for each of the nine cDNA libraries.
cDNA Library,
Biological Replicate (Rep)		Number of Paired
Sequences (in Millions)		Overall Mapping Alignment Rate
Control (t0), rep 124.994.96%
Control (t0), rep 226.589.61%
Control (t0), rep 340.693.34%
20 min sucrose (t20), rep 124.292.29%
20 min sucrose (t20), rep 247.391.00%
20 min sucrose (t20), rep 333.193.34%
40 min sucrose (t40), rep 119.587.91%
40 min sucrose (t40), rep 226.781.77%
40 min sucrose (t40), rep 316.980.44%
Table
Table 2. Ten most highly upregulated genes after 20 min (t20) and 40 min (t40) of sucrose, sorted for t20.
Table 2. Ten most highly upregulated genes after 20 min (t20) and 40 min (t40) of sucrose, sorted for t20.
Annotation/
G. max ID
(if Available)		Descriptiont20t40
Log2FCp-ValueLog2FCp-Value
auxin response factor 9 (ARF9)/transcriptional regulatorpotential mediators of auxin signaling in response to biotic and abiotic stress [38] 4.30.0074.76.92 × 10−5
aminopeptidase M1/
GLYMA.04G053300		metallopeptidase, involved in polar auxin transport [39] and root hair development [40]4.20.0034.00.015
ferredoxin-A/
GLYMA.05G168400		in roots, part of nitrogen assimilation via ferredoxin-nitrite reductase (NiR) [41]4.20.000023.56.76 × 10−6
phosphoglycerate mutase 1enzyme of glycolysis, induced by sucrose and auxin [42]40.000164.40.014
RESPONSE TO LOW SULFUR 3/GLYMA.04G225500function still unknown; possible transcriptional regulator involved in plant responses to environmental challenges [20]3.70.0412.60.023
LOC100305767/
GLYMA.08G158100		uncharacterized3.60.0434.60.014
thermospermine synthase ACAULIS5/
GLYMA.14G099200		involved in the synthesis of thermospermine, which may act as a plant growth regulator [43]; some thermospermine synthases are regulated by plant stress hormones [44]3.60.0203.90.005
Dicarboxylate carrier 1,
transporter of organic acids		may shuttle malate between the cytosol and mitochondria; induced in Fe-deficient roots [45]3.50.0495.10.012
J domain-containing protein/
GLYMA.08G074200		co-chaperones of Hsp70s (heat-shock proteins), likely involved in growth, development, and stress response [46]3.20.0454.30.014
RESPONSE TO LOW SULFUR 3/
GLYMA.06G139300		function still unknown; possible transcriptional regulator involved in plant responses to environmental challenges [20]3.20.0492.00.013
Table
Table 3. Ten most highly upregulated transcription factors (log2FC > 1.5, p-value < 0.5, indicated in bold) at t40. Data for t20 are also shown; while not statistically significant, some trends of upregulation can already be observed at 20 min of sucrose exposure.
Table 3. Ten most highly upregulated transcription factors (log2FC > 1.5, p-value < 0.5, indicated in bold) at t40. Data for t20 are also shown; while not statistically significant, some trends of upregulation can already be observed at 20 min of sucrose exposure.
AnnotationG. max ID
(if Available)		t20t40
log2FCp-Valuelog2FCp-Value
MYB14GLYMA.06G3002003.80.0996.30.018
WRKY41GLYMA.05G2159004.60.0545.70.007
MYB30GLYMA.06G3001002.40.1125.60.034
WRKY33GLYMA.03G0427003.20.0765.10.003
MYB30GLYMA.12G1048002.40.1144.60.015
WRKY SUSIBA2 2.40.1034.50.002
ERF 13GLYMA.03G1117001.20.1164.10.045
MYB13GLYMA.12G1046001.30.1184.10.045
ERF098-likeGLYMA.10G0366001.30.1793.90.00003
WRKY40GLYMA.07G0233002.10.0733.70.017
Table
Table 4. Top fifteen enriched GO categories among genes upregulated at 40 min of sucrose treatment.
Table 4. Top fifteen enriched GO categories among genes upregulated at 40 min of sucrose treatment.
GO TermExpressed
GO		Expected
Expression		Genome
GO Count		Corrected
p-Value *
response to chitin183511135 4 × 10−50
respiratory burst involved in
defense response		91194201.8 × 10−34
regulation of plant-type
hypersensitive response		1224610197 × 10−21
protein targeting the membrane1224610207.6 × 10−21
salicylic acid-mediated
signaling pathway		76214587.9 × 10−21
negative regulation of
programmed cell death		82245258.4 × 10−21
intracellular signal
transduction		76224791.3 × 10−19
response to wounding1194610315.1 × 10−19
negative regulation of the defense response97347665 × 10−18
MAPK cascade78265755.6 × 10−16
ethylene biosynthetic process50122702.7 × 10−15
jasmonic acid-mediated
signaling pathway		94368003 × 10−15
defense response to fungus104429414.7 × 10−15
systemic acquired resistance85316885.5 × 10−15
salicylic acid biosynthetic process79296532.56 × 10−13
* Bonferroni correction (Fisher test p-value multiplied by number of GO terms in gene list of interest).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nidumolu, L.C.M.; Lorilla, K.M.; Chakravarty, I.; Uhde-Stone, C. Soybean Root Transcriptomics: Insights into Sucrose Signaling at the Crossroads of Nutrient Deficiency and Biotic Stress Responses. Plants 2023, 12, 2117. https://doi.org/10.3390/plants12112117

AMA Style

Nidumolu LCM, Lorilla KM, Chakravarty I, Uhde-Stone C. Soybean Root Transcriptomics: Insights into Sucrose Signaling at the Crossroads of Nutrient Deficiency and Biotic Stress Responses. Plants. 2023; 12(11):2117. https://doi.org/10.3390/plants12112117

Chicago/Turabian Style

Nidumolu, Leela Chandra Manozna, Kristina Mae Lorilla, Indrani Chakravarty, and Claudia Uhde-Stone. 2023. "Soybean Root Transcriptomics: Insights into Sucrose Signaling at the Crossroads of Nutrient Deficiency and Biotic Stress Responses" Plants 12, no. 11: 2117. https://doi.org/10.3390/plants12112117

APA Style

Nidumolu, L. C. M., Lorilla, K. M., Chakravarty, I., & Uhde-Stone, C. (2023). Soybean Root Transcriptomics: Insights into Sucrose Signaling at the Crossroads of Nutrient Deficiency and Biotic Stress Responses. Plants, 12(11), 2117. https://doi.org/10.3390/plants12112117

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop