Dynamic gene expression changes in response to micronutrient, macronutrient, and multiple stress exposures in soybean

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

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 P_i 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.

Coupling VIGS with Short- and Long-Term Stress Exposure to Understand the Fiskeby III Iron Deficiency Stress Response

Yield loss due to abiotic stress is an increasing problem in agriculture. Soybean is a major crop for the upper Midwestern United States and calcareous soils exacerbate iron deficiency for growers, resulting in substantial yield losses. Fiskeby III is a soybean variety uniquely resistant to a variety of abiotic stresses, including iron deficiency. Previous studies identified a MATE transporter (Glyma.05G001700) associated with iron stress tolerance in Fiskeby III. To understand the function of this gene in the Fiskeby III response to iron deficiency, we coupled its silencing using virus-induced gene silencing with RNAseq analyses at two timepoints. Analyses of these data confirm a role for the MATE transporter in Fiskeby III iron stress responses. Further, they reveal that Fiskeby III induces transcriptional reprogramming within 24 h of iron deficiency stress, confirming that like other soybean varieties, Fiskeby III is able to quickly respond to stress. However, Fiskeby III utilizes novel genes and pathways in its iron deficiency response. Identifying and characterizing these genes and pathways in Fiskeby III provides novel targets for improving abiotic stress tolerance in elite soybean lines.

Ionomic Combined with Transcriptomic and Metabolomic Analyses to Explore the Mechanism Underlying the Effect of Melatonin in Relieving Nutrient Stress in Apple

Nutrient stress harms plant growth and yield. Melatonin is a biologically active, multifunctional hormone that relieves abiotic stress in plants. Although previous studies have shown that melatonin plays an important role in improving nutrient-use efficiency, the mechanism of its regulation of nutrient stress remains unclear. In this study, melatonin was applied to apple plants under nutrient stress, and morphological indices, physiological and biochemical indices, and stomatal morphology were evaluated. The response of apple plants to nutrient deficiency and the melatonin mechanism to alleviate nutrient stress were analyzed by combining ionome, transcriptome, and metabolome. The results showed that exogenous melatonin significantly alleviated the inhibitory effect of nutritional stress on the growth of apple plants by regulating stomatal morphology, improving antioxidant enzyme activity, promoting ion absorption, and utilizing and changing the absorption and distribution of minerals throughout the plant. The transcriptome results showed that melatonin alleviated nutrient stress and promoted nutrient absorption and utilization by regulating glutathione metabolism and upregulating some metal ion transport genes. The metabolome results indicated that levels of oxalic acid, L-ascorbic acid, anthocyanins (cyanidin-3-O-galactoside), lignans (lirioresinol A and syringaresinol), and melatonin significantly increased after exogenous melatonin was applied to plants under nutrient stress. These differentially expressed genes and the increase in beneficial metabolites may explain how melatonin alleviates nutrient stress in plants.

Environmental Pollutants and Oxidative Stress in Terrestrial and Aquatic Organisms: Examination of the Total Picture and Implications for Human Health

There is current great international concern about the contribution of environmental pollution to the global burden of disease particularly in the developing, low- and medium-income countries. Industrial activities, urbanization, developmental projects as well as various increased anthropogenic activities involving the improper generation, management and disposal of pollutants have rendered today’s environment highly polluted with various pollutants. These pollutants include toxic metals (lead, cadmium, mercury, arsenic), polycyclic aromatic hydrocarbons, polychlorinated biphenyls, pesticides and diesel exhaust particles most of which appear to be ubiquitous as well as have long-term environmental persistence with a wide range of toxicities such as oxidative stress among others. Oxidative stress, which may arise from increased production of damaging free radicals emanating from increased pollutant burden and depressed bioavailability of antioxidant defenses causes altered biochemical and physiological mechanisms and has been implicated in all known human pathologies most of which are chronic. Oxidative stress also affects both flora and fauna and plants are very important components of the terrestrial environment and significant contributors of nutrients for both man and animals. It is also remarkable that the aquatic environment in which sea animals and creatures are resident is also highly polluted, leading to aquatic stress that may affect the survival of the aquatic animals, sharing in the oxidative stress. These altered terrestrial and aquatic environments have an overarching effect on human health. Antioxidants neutralize the damaging free radicals thus, they play important protective roles in the onset, progression and severity of the unmitigated generation of pollutants that ultimately manifest as oxidative stress. Consequently, human health as well as that of aquatic and terrestrial organisms may be protected from environmental pollution by mitigating oxidative stress and employing the principles of nutritional medicine, essentially based on antioxidants derived mainly from plants, which serve as the panacea of the vicious state of environmental pollutants consequently, the health of the population. Understanding the total picture of oxidative stress and integrating the terrestrial and aquatic effects of environmental pollutants are central to sustainable health of the population and appear to require multi-sectoral collaborations from diverse disciplinary perspectives; basically the environmental, agricultural and health sectors.

Proteome Analysis of the Soybean Nodule Phosphorus Response Mechanism and Characterization of Stress-Induced Ribosome Structural and Protein Expression Changes

In agroecosystems, a plant-usable form of nitrogen is mainly generated by legume-based biological nitrogen fixation, a process that requires phosphorus (P) as an essential nutrient. To investigate the physiological mechanism whereby phosphorus influences soybean nodule nitrogen fixation, soybean root nodules were exposed to four phosphate levels: 1 mg/L (P stress), 11 mg/L (P stress), 31 mg/L (Normal P), and 61 mg/L (High P) then proteome analysis of nodules was conducted to identify phosphorus-associated proteome changes. We found that phosphorus stress-induced ribosomal protein structural changes were associated with altered key root nodule protein synthesis profiles. Importantly, up-regulated expression of peroxidase was observed as an important phosphorus stress-induced nitrogen fixation-associated adaptation that supported two nodule-associated activities: scavenging of reactive oxygen species (ROS) and cell wall growth. In addition, phosphorus transporter (PT) and purple acid phosphatase (PAPs) were up-regulated that regulated phosphorus transport and utilization to maintain phosphorus balance and nitrogen fixation function in phosphorus-stressed root nodules.

Global analysis of switchgrass (Panicum virgatum L.) transcriptomes in response to interactive effects of drought and heat stresses

BACKGROUND

Sustainable production of high-quality feedstock has been of great interest in bioenergy research. Despite the economic importance, high temperatures and water deficit are limiting factors for the successful cultivation of switchgrass in semi-arid areas. There are limited reports on the molecular basis of combined abiotic stress tolerance in switchgrass, particularly the combination of drought and heat stress. We used transcriptomic approaches to elucidate the changes in the response of switchgrass to drought and high temperature simultaneously.

RESULTS

We conducted solely drought treatment in switchgrass plant Alamo AP13 by withholding water after 45 days of growing. For the combination of drought and heat effect, heat treatment (35 °C/25 °C day/night) was imposed after 72 h of the initiation of drought. Samples were collected at 0 h, 72 h, 96 h, 120 h, 144 h, and 168 h after treatment imposition, total RNA was extracted, and RNA-Seq conducted. Out of a total of 32,190 genes, we identified 3912, as drought (DT) responsive genes, 2339 and 4635 as, heat (HT) and drought and heat (DTHT) responsive genes, respectively. There were 209, 106, and 220 transcription factors (TFs) differentially expressed under DT, HT and DTHT respectively. Gene ontology annotation identified the metabolic process as the significant term enriched in DTHT genes. Other biological processes identified in DTHT responsive genes included: response to water, photosynthesis, oxidation-reduction processes, and response to stress. KEGG pathway enrichment analysis on DT and DTHT responsive genes revealed that TFs and genes controlling phenylpropanoid pathways were important for individual as well as combined stress response. For example, hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase (HCT) from the phenylpropanoid pathway was induced by single DT and combinations of DTHT stress.

CONCLUSION

Through RNA-Seq analysis, we have identified unique and overlapping genes in response to DT and combined DTHT stress in switchgrass. The combination of DT and HT stress may affect the photosynthetic machinery and phenylpropanoid pathway of switchgrass which negatively impacts lignin synthesis and biomass production of switchgrass. The biological function of genes identified particularly in response to DTHT stress could further be confirmed by techniques such as single point mutation or RNAi.

New insights into the evolution of SPX gene family from algae to legumes; a focus on soybean

BACKGROUND

SPX-containing proteins have been known as key players in phosphate signaling and homeostasis. In Arabidopsis and rice, functions of some SPXs have been characterized, but little is known about their function in other plants, especially in the legumes.

RESULTS

We analyzed SPX gene family evolution in legumes and in a number of key species from algae to angiosperms. We found that SPX harboring proteins showed fluctuations in domain fusions from algae to the angiosperms with, finally, four classes appearing and being retained in the land plants. Despite these fluctuations, Lysine Surface Cluster (KSC), and the third residue of Phosphate Binding Sites (PBS) showed complete conservation in almost all of SPXs except few proteins in Selaginella moellendorffii and Papaver sumniferum, suggesting they might have different ligand preferences. In addition, we found that the WGD/segmentally or dispersed duplication types were the most frequent contributors to the SPX expansion, and that there is a positive correlation between the amount of WGD contribution to the SPX expansion in individual species and its number of EXS genes. We could also reveal that except SPX class genes, other classes lost the collinearity relationships among Arabidopsis and legume genomes. The sub- or neo-functionalization of the duplicated genes in the legumes makes it difficult to find the functional orthologous genes. Therefore, we used two different methods to identify functional orthologs in soybean and Medicago. High variance in the dynamic and spatial expression pattern of GmSPXs proved the new or sub-functionalization in the paralogs.

CONCLUSION

This comprehensive analysis revealed how SPX gene family evolved from algae to legumes and also discovered several new domains fused to SPX domain in algae. In addition, we hypothesized that there different phosphate sensing mechanisms might occur in S. moellendorffii and P. sumniferum. Finally, we predicted putative functional orthologs of AtSPXs in the legumes, especially, orthologs of AtPHO1, involved in long-distance Pi transportation. These findings help to understand evolution of phosphate signaling and might underpin development of new legume varieties with improved phosphate use efficiency.

Comparative Omics Analysis of Brassica napus Roots Subjected to Six Individual Macronutrient Deprivations Reveals Deficiency-Specific Genes and Metabolomic Profiles

The early and specific diagnosis of a macronutrient deficiency is challenging when seeking to better manage fertilizer inputs in the context of sustainable agriculture. Consequently, this study explored the potential for transcriptomic and metabolomic analysis of Brassica napus roots to characterize the effects of six individual macronutrient deprivations (N, Mg, P, S, K, and Ca). Our results showed that before any visual phenotypic response, all macronutrient deprivations led to a large modulation of the transcriptome and metabolome involved in various metabolic pathways, and some were common to all macronutrient deprivations. Significantly, comparative transcriptomic analysis allowed the definition of a subset of 3282, 2011, 6325, 1384, 439, and 5157 differentially expressed genes (DEGs) specific to N, Mg, P, S, K, and Ca deprivations, respectively. Surprisingly, gene ontology term enrichment analysis performed on this subset of specific DEGs highlighted biological processes that are common to a number of these macronutrient deprivations, illustrating the complexity of nutrient interactions. In addition, a set of 38 biochemical compounds that discriminated the macronutrient deprivations was identified using a metabolic approach. The opportunity to use these specific DEGs and/or biochemical compounds as potential molecular indicators to diagnose macronutrient deficiency is discussed.

Comparing Early Transcriptomic Responses of 18 Soybean (Glycine max) Genotypes to Iron Stress

Iron deficiency chlorosis (IDC) is an abiotic stress that negatively affects soybean (Glycine max [L.] Merr.) production. Much of our knowledge of IDC stress responses is derived from model plant species. Gene expression, quantitative trait loci (QTL) mapping, and genome-wide association studies (GWAS) performed in soybean suggest that stress response differences exist between model and crop species. Our current understanding of the molecular response to IDC in soybeans is largely derived from gene expression studies using near-isogenic lines differing in iron efficiency. To improve iron efficiency in soybeans and other crops, we need to expand gene expression studies to include the diversity present in germplasm collections. Therefore, we collected 216 purified RNA samples (18 genotypes, two tissue types [leaves and roots], two iron treatments [sufficient and deficient], three replicates) and used RNA sequencing to examine the expression differences of 18 diverse soybean genotypes in response to iron deficiency. We found a rapid response to iron deficiency across genotypes, most responding within 60 min of stress. There was little evidence of an overlap of specific differentially expressed genes, and comparisons of gene ontology terms and transcription factor families suggest the utilization of different pathways in the stress response. These initial findings suggest an untapped genetic potential within the soybean germplasm collection that could be used for the continued improvement of iron efficiency in soybean.

Mining Fiskeby III and Mandarin (Ottawa) Expression Profiles to Understand Iron Stress Tolerant Responses in Soybean

The soybean (Glycine max L. merr) genotype Fiskeby III is highly resistant to a multitude of abiotic stresses, including iron deficiency, incurring only mild yield loss during stress conditions. Conversely, Mandarin (Ottawa) is highly susceptible to disease and suffers severe phenotypic damage and yield loss when exposed to abiotic stresses such as iron deficiency, a major challenge to soybean production in the northern Midwestern United States. Using RNA-seq, we characterize the transcriptional response to iron deficiency in both Fiskeby III and Mandarin (Ottawa) to better understand abiotic stress tolerance. Previous work by our group identified a quantitative trait locus (QTL) on chromosome 5 associated with Fiskeby III iron efficiency, indicating Fiskeby III utilizes iron deficiency stress mechanisms not previously characterized in soybean. We targeted 10 of the potential candidate genes in the Williams 82 genome sequence associated with the QTL using virus-induced gene silencing. Coupling virus-induced gene silencing with RNA-seq, we identified a single high priority candidate gene with a significant impact on iron deficiency response pathways. Characterization of the Fiskeby III responses to iron stress and the genes underlying the chromosome 5 QTL provides novel targets for improved abiotic stress tolerance in soybean.

Copper promoting oyster larval growth and settlement: Molecular insights from RNA-seq

As a cofactor of key enzymes, Cu is required in living organisms, although Cu levels in the natural environment are typically low. In this study, the promotion of growth and settlement on the larvae of oyster Crassostrea angulata was observed at an environmentally relevant concentration (10 μg/L Cu). Interestingly, Cu accumulation in the soft tissue of oyster larvae increased during the larval development and exhibited a sharp increase at the late pelagic stage. With the help of RNA-seq, we constructed a high-quality transcriptional database of the oyster C. angulata larvae (24,257 genes with an average length of 1594 bp) via de novo assembly, which provided the basic molecular changes during the larval development. Network analysis of five early developmental stages and differential expression under Cu exposure were integrated to examine the roles of Cu in oyster larvae. Our molecular analysis demonstrated that both ion channels and organic transporters contributed to Cu internalization from the external environment, including zinc transporters and amino acid transporters. The followed distribution of Cu across cells was achieved by ATP7A, the circulatory system, and the Cu transporters (CTRs). Cu exposure enhanced the ribosome and the calcium binding proteins with a higher rate of translation and shell formation, giving rise to faster growth of oyster larvae. Furthermore, Cu facilitated the settling process by upregulating the chitin binding genes and then promoting the formation of the proteinaceous matrix between larvae and substrate. Our study presents the molecular basis for Cu promotion (i.e., hormesis) effects on oyster larval growth and settlement.

Proteome characterization of two contrasting soybean genotypes in response to different phosphorus treatments

Phosphorus (P) is an essential element for the growth and development of plants. Soybean (Glycine max) is an important food crop that is grown worldwide. Soybean yield is significantly affected by P deficiency in the soil. To investigate the molecular factors that determine the response and tolerance at low-P in soybean, we conducted a comparative proteomics study of a genotype with low-P tolerance (Liaodou 13, L13) and a genotype with low-P sensitivity (Tiefeng 3, T3) in a paper culture experiment with three P treatments, i.e. P-free (0 mmol·L-1), low-P (0.05 mmol·L-1) and normal-P (0.5 mmol·L-1). A total of 4126 proteins were identified in roots of the two genotypes. Increased numbers of differentially expressed proteins (DEPs) were obtained from low-P to P-free conditions compared to the normal-P treatment. All DEPs obtained in L13 (660) were upregulated in response to P deficiency, while most DEPs detected in T3 (133) were downregulated under P deficiency. Important metabolic pathways such as oxidative phosphorylation, glutathione metabolism and carbon metabolism were suppressed in T3, which could have affected the survival of the plants in P-limited soil. In contrast, L13 increased the metabolic activity in the 2-oxocarboxylic acid metabolism, carbon metabolism, glycolysis, biosynthesis of amino acids, pentose phosphatase, oxidative phosphorylation, other types of O-glycan biosynthesis and riboflavin metabolic pathways in order to maintain normal plant growth under P deficiency. Three key proteins I1KW20 (prohibitins), I1K3U8 (alpha-amylase inhibitors) and C6SZ93 (alpha-amylase inhibitors) were suggested as potential biomarkers for screening soybean genotypes with low-P tolerance. Overall, this study provides new insights into the response and tolerance to P deficiency in soybean.

Gene Expression Responses to Sequential Nutrient Deficiency Stresses in Soybean

Throughout the growing season, crops experience a multitude of short periods of various abiotic stresses. These stress events have long-term impacts on plant performance and yield. It is imperative to improve our understanding of the genes and biological processes underlying plant stress tolerance to mitigate end of season yield loss. The majority of studies examining transcriptional changes induced by stress focus on single stress events. Few studies have been performed in model or crop species to examine transcriptional responses of plants exposed to repeated or sequential stress exposure, which better reflect field conditions. In this study, we examine the transcriptional profile of soybean plants exposed to iron deficiency stress followed by phosphate deficiency stress (-Fe-Pi). Comparing this response to previous studies, we identified a core suite of genes conserved across all repeated stress exposures (-Fe-Pi, -Fe-Fe, -Pi-Pi). Additionally, we determined transcriptional response to sequential stress exposure (-Fe-Pi) involves genes usually associated with reproduction, not stress responses. These findings highlight the plasticity of the plant transcriptome and the complexity of unraveling stress response pathways.

Examining Short-Term Responses to a Long-Term Problem: RNA-Seq Analyses of Iron Deficiency Chlorosis Tolerant Soybean

Iron deficiency chlorosis (IDC) is a global crop production problem, significantly impacting yield. However, most IDC studies have focused on model species, not agronomically important crops. Soybean is the second largest crop grown in the United States, yet the calcareous soils across most of the upper U.S. Midwest limit soybean growth and profitability. To understand early soybean iron stress responses, we conducted whole genome expression analyses (RNA-sequencing) of leaf and root tissue from the iron efficient soybean (Glycine max) cultivar Clark, at 30, 60 and 120 min after transfer to iron stress conditions. We identified over 10,000 differentially expressed genes (DEGs), with the number of DEGs increasing over time in leaves, but decreasing over time in roots. To investigate these responses, we clustered our expression data across time to identify suites of genes, their biological functions, and the transcription factors (TFs) that regulate their expression. These analyses reveal the hallmarks of the soybean iron stress response (iron uptake and homeostasis, defense, and DNA replication and methylation) can be detected within 30 min. Furthermore, they suggest root to shoot signaling initiates early iron stress responses representing a novel paradigm for crop stress adaptations.