The molecular chaperone binding protein BiP prevents leaf dehydration-induced cellular homeostasis disruption

In-Depth Characterization of bZIP Genes in the Context of Endoplasmic Reticulum (ER) Stress in Brassica campestris ssp. chinensis

Numerous studies have been conducted to investigate the genomic characterization of bZIP genes and their involvement in the cellular response to endoplasmic reticulum (ER) stress. These studies have provided valuable insights into the coordinated cellular response to ER stress, which is mediated by bZIP transcription factors (TFs). However, a comprehensive and systematic investigations regarding the role of bZIP genes and their involvement in ER stress response in pak choi is currently lacking in the existing literature. To address this knowledge gap, the current study was initiated to elucidate the genomic characteristics of bZIP genes, gain insight into their expression patterns during ER stress in pak choi, and investigate the protein-to-protein interaction of bZIP genes with the ER chaperone BiP. In total, 112 members of the BcbZIP genes were identified through a comprehensive genome-wide analysis. Based on an analysis of sequence similarity, gene structure, conserved domains, and responsive motifs, the identified BcbZIP genes were categorized into 10 distinct subfamilies through phylogenetic analysis. Chromosomal location and duplication events provided insight into their genomic context and evolutionary history. Divergence analysis estimated their evolutionary history with a predicted divergence time ranging from 0.73 to 80.71 million years ago (MYA). Promoter regions of the BcbZIP genes were discovered to exhibit a wide variety of cis-elements, including light, hormone, and stress-responsive elements. GO enrichment analysis further confirmed their roles in the ER unfolded protein response (UPR), while co-expression network analysis showed a strong relationship of BcbZIP genes with ER-stress-responsive genes. Moreover, gene expression profiles and protein-protein interaction with ER chaperone BiP further confirmed their roles and capacity to respond to ER stress in pak choi.

Physiological and transcriptomic analyses reveal that phytohormone pathways and glutathione metabolism are involved in the arsenite toxicity response in tomatoes

The main forms of inorganic arsenic (As) in soil are arsenate [As(V)] and arsenite [As(III)]. Both forms inhibit plant growth. Here, we investigate the effects of As(III) toxicity on the growth of tomatoes by integrating physiological and transcriptomic analyses. As(III) toxicity induces oxidative damage, inhibits photosynthetic efficiency, and reduces soluble sugar levels. As(III) toxicity leads to reductions in auxin, cytokinin and jasmonic acid contents by 29 %, 39 % and 55 %, respectively, but leads to increases in the ethylene precursor 1-amino-cyclopropane carboxylic acid, abscisic acid and salicylic acid contents in roots, by 116 %, 79 % and 39 %, respectively, thereby altering phytohormone signalling pathways. The total glutathione, reduced glutathione (GSH) and oxidized glutathione (GSSG) contents are reduced by 59 %, 49 % and 94 % in roots; moreover, a high GSH/GSSG ratio is maintained through increased glutathione reductase activity (increased by 214 %) and decreased glutathione peroxidase activity (decreased by 40 %) in the roots of As(III)-treated tomato seedlings. In addition, As(III) toxicity affects the expression of genes related to the endoplasmic reticulum stress response. The altered expression of aquaporins and ABCC transporters changes the level of As(III) accumulation in plants. A set of hub genes involved in modulating As(III) toxicity responses in tomatoes was identified via a weighted gene coexpression network analysis. Taken together, these results elucidate the physiological and molecular regulatory mechanism underlying As(III) toxicity and provide a theoretical basis for selecting and breeding tomato varieties with low As(III) accumulation. Therefore, these findings are expected to be helpful in improving food safety and to developing sustainable agricultural.

The miRNA-mRNA Regulatory Modules of Pinus massoniana Lamb. in Response to Drought Stress

Masson pine (Pinus massoniana Lamb.) is a major fast-growing woody tree species and pioneer species for afforestation in barren sites in southern China. However, the regulatory mechanism of gene expression in P. massoniana under drought remains unclear. To uncover candidate microRNAs, their expression profiles, and microRNA-mRNA interactions, small RNA-seq was used to investigate the transcriptome from seedling roots under drought and rewatering in P. massoniana. A total of 421 plant microRNAs were identified. Pairwise differential expression analysis between treatment and control groups unveiled 134, 156, and 96 differential expressed microRNAs at three stages. These constitute 248 unique microRNAs, which were subsequently categorized into six clusters based on their expression profiles. Degradome sequencing revealed that these 248 differentially expressed microRNAs targeted 2069 genes. Gene Ontology enrichment analysis suggested that these target genes were related to translational and posttranslational regulation, cell wall modification, and reactive oxygen species scavenging. miRNAs such as miR482, miR398, miR11571, miR396, miR166, miRN88, and miRN74, along with their target genes annotated as F-box/kelch-repeat protein, 60S ribosomal protein, copper-zinc superoxide dismutase, luminal-binding protein, S-adenosylmethionine synthase, and Early Responsive to Dehydration Stress may play critical roles in drought response. This study provides insights into microRNA responsive to drought and rewatering in Masson pine and advances the understanding of drought tolerance mechanisms in Pinus.

Regulation of the unfolded protein response during dehydration stress in African clawed frogs, Xenopus laevis

The unfolded protein response (UPR) is a wide-ranging cellular response to accumulation of malfolded proteins in the endoplasmic reticulum (ER) and acts as a quality control mechanism to halt protein processing and repair/destroy malfolded proteins under stress conditions of many kinds. Among vertebrate species, amphibians experience the greatest challenges in maintaining water and osmotic balance, the high permeability of their skin making them very susceptible to dehydration and challenging their ability to maintain cellular homeostasis. The present study evaluates the involvement of the UPR in dealing with dehydration-mediated disruption of protein processing in the tissues of African clawed frogs, Xenopus laevis. This primarily aquatic frog must deal with seasonal drought conditions in its native southern Africa environment. Key markers of cellular stress that impact protein processing were identified in six tissues of frogs that had lost 28% of total body water, as compared with fully hydrated controls. This included upregulation of glucose-regulated proteins (GRPs) that are resident chaperones in the ER, particularly 2-ninefold increases in GRP58, GRP75, and/or GRP94 in the lung and skin. Activating transcription factors (ATF3, ATF4, ATF6) that mediate UPR responses also responded to dehydration stress, particularly in skeletal muscle where both ATF3 and ATF4 rose strongly in the nucleus. Other protein markers of the UPR including GADD34, GADD153, EDEM, and XBP-1 also showed selective upregulation in frog tissues in response to dehydration and nuclear levels of the transcription factors XBP-1 and P-CREB rose indicating up-regulation of genes under their control.

Novel molecules and target genes for vegetative heat tolerance in wheat

To prevent yield losses caused by climate change, it is important to identify naturally tolerant genotypes with traits and related pathways that can be targeted for crop improvement. Here we report on the characterization of contrasting vegetative heat tolerance in two UK bread wheat varieties. Under chronic heat stress, the heat-tolerant cultivar Cadenza produced an excessive number of tillers which translated into more spikes and higher grain yield compared to heat-sensitive Paragon. RNAseq and metabolomics analyses revealed that over 5000 genotype-specific genes were differentially expressed, including photosynthesis-related genes, which might explain the observed ability of Cadenza to maintain photosynthetic rate under heat stress. Around 400 genes showed a similar heat-response in both genotypes. Only 71 genes showed a genotype × temperature interaction. As well as known heat-responsive genes such as heat shock proteins (HSPs), several genes that have not been previously linked to the heat response, particularly in wheat, have been identified, including dehydrins, ankyrin-repeat protein-encoding genes, and lipases. Contrary to primary metabolites, secondary metabolites showed a highly differentiated heat response and genotypic differences. These included benzoxazinoid (DIBOA, DIMBOA), and phenylpropanoids and flavonoids with known radical scavenging capacity, which was assessed via the DPPH assay. The most highly heat-induced metabolite was (glycosylated) propanediol, which is widely used in industry as an anti-freeze. To our knowledge, this is the first report on its response to stress in plants. The identified metabolites and candidate genes provide novel targets for the development of heat-tolerant wheat.

Cadmium-mediated toxicity in plant cells is associated with the DCD/NRP-mediated cell death response

Cadmium (Cd2+ ) is highly harmful to plant growth. Although Cd2+ induces programmed cell death (PCD) in plant cells, Cd2+ stress in whole plants during later developmental stages and the mechanism underlying Cd2+ -mediated toxicity are poorly understood. Here, we showed that Cd2+ limits plant growth, causes intense redness in leaf vein, leaf yellowing, and chlorosis during the R1 reproductive stage of soybean (Glycine max). These symptoms were associated with Cd2+ -induced PCD, as Cd2+ -stressed soybean leaves displayed decreased number of nuclei, enhanced cell death, DNA damage, and caspase 1 activity compared to unstressed leaves. Accordingly, Cd2+ -induced NRPs, GmNAC81, GmNAC30 and VPE, the DCD/NRP-mediated cell death signalling components, which execute PCD via caspase 1-like VPE activity. Furthermore, overexpression of the positive regulator of this cell death signalling GmNAC81 enhanced sensitivity to Cd2+ stress and intensified the hallmarks of Cd2+ -mediated PCD. GmNAC81 overexpression enhanced Cd2+ -induced H2 O2 production, cell death, DNA damage, and caspase-1-like VPE expression. Conversely, BiP overexpression negatively regulated the NRPs/GmNACs/VPE signalling module, conferred tolerance to Cd2+ stress and reduced Cd2+ -mediated cell death. Collectively, our data indicate that Cd2+ induces PCD in plants via activation of the NRP/GmNAC/VPE regulatory circuit that links developmentally and stress-induced cell death.

Endoplasmic Reticulum Stress and Unfolded Protein Response Signaling in Plants

Plants are sensitive to a variety of stresses that cause various diseases throughout their life cycle. However, they have the ability to cope with these stresses using different defense mechanisms. The endoplasmic reticulum (ER) is an important subcellular organelle, primarily recognized as a checkpoint for protein folding. It plays an essential role in ensuring the proper folding and maturation of newly secreted and transmembrane proteins. Different processes are activated when around one-third of newly synthesized proteins enter the ER in the eukaryote cells, such as glycosylation, folding, and/or the assembling of these proteins into protein complexes. However, protein folding in the ER is an error-prone process whereby various stresses easily interfere, leading to the accumulation of unfolded/misfolded proteins and causing ER stress. The unfolded protein response (UPR) is a process that involves sensing ER stress. Many strategies have been developed to reduce ER stress, such as UPR, ER-associated degradation (ERAD), and autophagy. Here, we discuss the ER, ER stress, UPR signaling and various strategies for reducing ER stress in plants. In addition, the UPR signaling in plant development and different stresses have been discussed.

Contrasting roles of GmNAC065 and GmNAC085 in natural senescence, plant development, multiple stresses and cell death responses

NACs are plant-specific transcription factors involved in controlling plant development, stress responses, and senescence. As senescence-associated genes (SAGs), NACs integrate age- and stress-dependent pathways that converge to programmed cell death (PCD). In Arabidopsis, NAC-SAGs belong to well-characterized regulatory networks, poorly understood in soybean. Here, we interrogated the soybean genome and provided a comprehensive analysis of senescence-associated Glycine max (Gm) NACs. To functionally examine GmNAC-SAGs, we selected GmNAC065, a putative ortholog of Arabidopsis ANAC083/VNI2 SAG, and the cell death-promoting GmNAC085, an ANAC072 SAG putative ortholog, for analyses. Expression analysis of GmNAC065 and GmNAC085 in soybean demonstrated (i) these cell death-promoting GmNACs display contrasting expression changes during age- and stress-induced senescence; (ii) they are co-expressed with functionally different gene sets involved in stress and PCD, and (iii) are differentially induced by PCD inducers. Furthermore, we demonstrated GmNAC065 expression delays senescence in Arabidopsis, a phenotype associated with enhanced oxidative performance under multiple stresses, higher chlorophyll, carotenoid and sugar contents, and lower stress-induced PCD compared to wild-type. In contrast, GmNAC085 accelerated stress-induced senescence, causing enhanced chlorophyll loss, ROS accumulation and cell death, decreased antioxidative system expression and activity. Accordingly, GmNAC065 and GmNAC085 targeted functionally contrasting sets of downstream AtSAGs, further indicating that GmNAC85 and GmNAC065 regulators function inversely in developmental and environmental PCD.

BiP-overexpressing soybean plants display accelerated hypersensitivity response (HR) affecting the SA-dependent sphingolipid and flavonoid pathways

Biotic and abiotic environmental stresses have limited the increase in soybean productivity. Overexpression of the molecular chaperone BiP in transgenic plants has been associated with the response to osmotic stress and drought tolerance by maintaining cellular homeostasis and delaying hypersensitive cell death. Here, we evaluated the metabolic changes in response to the hypersensitivity response (HR) caused by the non-compatible bacteria Pseudomonas syringae pv. tomato in BiP-overexpressing plants. The HR-modified metabolic profiles in BiP-overexpressing plants were significantly distinct from the wild-type untransformed. The transgenic plants displayed a lower abundance of HR-responsive metabolites as amino acids, sugars, carboxylic acids and signal molecules, including p-aminobenzoic acid (PABA) and dihydrosphingosine (DHS), when compared to infected wild-type plants. In contrast, salicylic acid (SA) biosynthetic and signaling pathways were more stimulated in transgenic plants, and both pathogenesis-related genes (PRs) and transcriptional factors controlling the SA pathway were more induced in the BiP-overexpressing lines. Furthermore, the long-chain bases (LCBs) and ceramide biosynthetic pathways showed alterations in gene expression and metabolite abundance. Thus, as a protective pathway against pathogens, HR regulation by sphingolipids and SA may account at least in part by the enhanced resistance of transgenic plants. GmNAC32 transcriptional factor was more induced in the transgenic plants and it has also been reported to regulate flavonoid synthesis in response to SA. In fact, the BiP-overexpressing plants showed an increase in flavonoids, mainly prenylated isoflavones, as precursors for phytoalexins. Our results indicate that the BiP-mediated acceleration in the hypersensitive response may be a target for metabolic engineering of plant resistance against pathogens.

Transcriptomics of Biostimulation of Plants Under Abiotic Stress

Plant biostimulants are compounds, living microorganisms, or their constituent parts that alter plant development programs. The impact of biostimulants is manifested in several ways: via morphological, physiological, biochemical, epigenomic, proteomic, and transcriptomic changes. For each of these, a response and alteration occur, and these alterations in turn improve metabolic and adaptive performance in the environment. Many studies have been conducted on the effects of different biotic and abiotic stimulants on plants, including many crop species. However, as far as we know, there are no reviews available that describe the impact of biostimulants for a specific field such as transcriptomics, which is the objective of this review. For the commercial registration process of products for agricultural use, it is necessary to distinguish the specific impact of biostimulants from that of other legal categories of products used in agriculture, such as fertilizers and plant hormones. For the chemical or biological classification of biostimulants, the classification is seen as a complex issue, given the great diversity of compounds and organisms that cause biostimulation. However, with an approach focused on the impact on a particular field such as transcriptomics, it is perhaps possible to obtain a criterion that allows biostimulants to be grouped considering their effects on living systems, as well as the overlap of the impact on metabolism, physiology, and morphology occurring between fertilizers, hormones, and biostimulants.

GmNAC81 Inversely Modulates Leaf Senescence and Drought Tolerance

Glycine max NAC81 (GmNAC81) is a downstream effector of the DCD/NRP-mediated cell death signaling, which interacts with GmNAC30 to fully induce the caspase 1-like vacuolar processing enzyme (VPE) expression, the executioner of the cell death program. GmNAC81 has been previously shown to positively modulate leaf senescence via the NRP/GmNAC81/VPE signaling module. Here, we examined the transcriptome induced by GmNAC81 overexpression and leaf senescence and showed that GmNAC81 further modulates leaf senescence by regulating an extensive repertoire of functionally characterized senescence-associated genes (SAGs). Because the NRP/GmNAC81/VPE signaling circuit also relays stress-induced cell death signals, we examined the effect of GmNAC81 overexpression in drought responses. Enhanced GmNAC81 expression in the transgenic lines increased sensitivity to water deprivation. Under progressive drought, the GmNAC81-overexpressing lines displayed severe leaf wilting, a larger and faster decline in leaf Ψw, relative water content (RWC), photosynthesis rate, stomatal conductance, and transpiration rate, in addition to higher Ci/Ca and lower Fm/Fv ratios compared to the BR16 control line. Collectively, these results indicate that the photosynthetic activity and apparatus were more affected by drought in the transgenic lines. Consistent with hypersensitivity to drought, chlorophyll loss, and lipid peroxidation were higher in the GmNAC81-overexpressing lines than in BR16 under dehydration. In addition to inducing VPE expression, GmNAC81 overexpression uncovered the regulation of typical drought-responsive genes. In particular, key regulators and effectors of ABA signaling were suppressed by GmNAC81 overexpression. These results suggest that GmNAC81 may negatively control drought tolerance not only via VPE activation but also via suppression of ABA signaling.

Genome-wide identification and characterization of Solanum tuberosum BiP genes reveal the role of the promoter architecture in BiP gene diversity

The endoplasmic reticulum (ER) immunoglobulin binding proteins (BiPs) are molecular chaperones involved in normal protein maturation and refolding malformed proteins through the unfolded protein response (UPR). Plant BiPs belong to a multi-gene family contributing to development, immunity, and responses to environmental stresses. This study identified three BiP homologs in the Solanum tuberosum (potato) genome using phylogenetic, amino acid sequence, 3-D protein modeling, and gene structure analysis. These analyses revealed that StBiP1 and StBiP2 grouped with AtBiP2, whereas StBiP3 grouped with AtBiP3. While the protein sequences and folding structures are highly similar, these StBiPs are distinguishable by their expression patterns in different tissues and in response to environmental stressors such as treatment with heat, chemicals, or virus elicitors of UPR. Ab initio promoter analysis revealed that potato and Arabidopsis BiP1 and BiP2 promoters were highly enriched with cis-regulatory elements (CREs) linked to developmental processes, whereas BiP3 promoters were enriched with stress related CREs. The frequency and linear distribution of these CREs produced two phylogenetic branches that further resolve the groups identified through gene phylogeny and exon/intron phase analysis. These data reveal that the CRE architecture of BiP promoters potentially define their spatio-temporal expression patterns under developmental and stress related cues.

Genome-wide identification and characterization of Solanum tuberosum BiP genes reveal the role of the promoter architecture in BiP gene diversity

The endoplasmic reticulum (ER) immunoglobulin binding proteins (BiPs) are molecular chaperones involved in normal protein maturation and refolding malformed proteins through the unfolded protein response (UPR). Plant BiPs belong to a multi-gene family contributing to development, immunity, and responses to environmental stresses. This study identified three BiP homologs in the Solanum tuberosum (potato) genome using phylogenetic, amino acid sequence, 3-D protein modeling, and gene structure analysis. These analyses revealed that StBiP1 and StBiP2 grouped with AtBiP2, whereas StBiP3 grouped with AtBiP3. While the protein sequences and folding structures are highly similar, these StBiPs are distinguishable by their expression patterns in different tissues and in response to environmental stressors such as treatment with heat, chemicals, or virus elicitors of UPR. Ab initio promoter analysis revealed that potato and Arabidopsis BiP1 and BiP2 promoters were highly enriched with cis-regulatory elements (CREs) linked to developmental processes, whereas BiP3 promoters were enriched with stress related CREs. The frequency and linear distribution of these CREs produced two phylogenetic branches that further resolve the groups identified through gene phylogeny and exon/intron phase analysis. These data reveal that the CRE architecture of BiP promoters potentially define their spatio-temporal expression patterns under developmental and stress related cues.

Comprehensive mapping of abiotic stress inputs into the soybean circadian clock

The plant circadian clock evolved to increase fitness by synchronizing physiological processes with environmental oscillations. Crop fitness was artificially selected through domestication and breeding, and the circadian clock was identified by both natural and artificial selections as a key to improved fitness. Despite progress in Arabidopsis, our understanding of the crop circadian clock is still limited, impeding its rational improvement for enhanced fitness. To unveil the interactions between the crop circadian clock and various environmental cues, we comprehensively mapped abiotic stress inputs to the soybean circadian clock using a 2-module discovery pipeline. Using the "molecular timetable" method, we computationally surveyed publicly available abiotic stress-related soybean transcriptomes to identify stresses that have strong impacts on the global rhythm. These findings were then experimentally confirmed using a multiplexed RNA sequencing technology. Specific clock components modulated by each stress were further identified. This comprehensive mapping uncovered inputs to the plant circadian clock such as alkaline stress. Moreover, short-term iron deficiency targeted different clock components in soybean and Arabidopsis and thus had opposite effects on the clocks of these 2 species. Comparing soybean varieties with different iron uptake efficiencies suggests that phase modulation might be a mechanism to alleviate iron deficiency symptoms in soybean. These unique responses in soybean demonstrate the need to directly study crop circadian clocks. Our discovery pipeline may serve as a broadly applicable tool to facilitate these explorations.

Mechanism of the drought tolerance of a transgenic soybean overexpressing the molecular chaperone BiP

Drought is one of major constraints that limits agricultural productivity. Some factors, including climate changes and acreage expansion, indicates towards the need for developing drought tolerant genotypes. In addition to its protective role against endoplasmic reticulum (ER) stress, we have previously shown that the molecular chaperone binding protein (BiP) is involved in the response to osmotic stress and promotes drought tolerance. Here, we analyzed the proteomic and metabolic profiles of BiP-overexpressing transgenic soybean plants and the corresponding untransformed line under drought conditions by 2DE-MS and GC/MS. The transgenic plant showed lower levels of the abscisic acid and jasmonic acid as compared to untransformed plants both in irrigated and non-irrigated conditions. In contrast, the level of salicylic acid was higher in transgenic lines than in untransformed line, which was consistent with the antagonistic responses mediated by these phytohormones. The transgenic plants displayed a higher abundance of photosynthesis-related proteins, which gave credence to the hypothesis that these transgenic plants could survive under drought conditions due to their genetic modification and altered physiology. The proteins involved in pathways related to respiration, glycolysis and oxidative stress were not signifcantly changed in transgenic plants as compared to untransformed genotype, which indicate a lower metabolic perturbation under drought of the engineered genotype. The transgenic plants may have adopted a mechanism of drought tolerance by accumulating osmotically active solutes in the cell. As evidenced by the metabolic profiles, the accumulation of nine primary amino acids by protein degradation maintained the cellular turgor in the transgenic genotype under drought conditions. Thus, this mechanism of protection may cause the physiological activities including photosynthesis to be active under drought conditions.

Revisiting the Soybean GmNAC Superfamily

The NAC (NAM, ATAF, and CUC) genes encode transcription factors involved with the control of plant morph-physiology and stress responses. The release of the last soybean (Glycine max) genome assembly (Wm82.a2.v1) raised the possibility that new NAC genes would be present in the soybean genome. Here, we interrogated the last version of the soybean genome against a conserved NAC domain structure. Our analysis identified 32 putative novel NAC genes, updating the superfamily to 180 gene members. We also organized the genes in 15 phylogenetic subfamilies, which showed a perfect correlation among sequence conservation, expression profile, and function of orthologous Arabidopsis thaliana genes and NAC soybean genes. To validate our in silico analyses, we monitored the stress-mediated gene expression profiles of eight new NAC-genes by qRT-PCR and monitored the GmNAC senescence-associated genes by RNA-seq. Among ER stress, osmotic stress and salicylic acid treatment, all the novel tested GmNAC genes responded to at least one type of stress, displaying a complex expression profile under different kinetics and extension of the response. Furthermore, we showed that 40% of the GmNACs were differentially regulated by natural leaf senescence, including eight (8) newly identified GmNACs. The developmental and stress-responsive expression profiles of the novel NAC genes fitted perfectly with their phylogenetic subfamily. Finally, we examined two uncharacterized senescence-associated proteins, GmNAC065 and GmNAC085, and a novel, previously unidentified, NAC protein, GmNAC177, and showed that they are nuclear localized, and except for GmNAC065, they display transactivation activity in yeast. Consistent with a role in leaf senescence, transient expression of GmNAC065 and GmNAC085 induces the appearance of hallmarks of leaf senescence, including chlorophyll loss, leaf yellowing, lipid peroxidation and accumulation of H₂O₂. GmNAC177 was clustered to an uncharacterized subfamily but in close proximity to the TIP subfamily. Accordingly, it was rapidly induced by ER stress and by salicylic acid under late kinetic response and promoted cell death in planta. Collectively, our data further substantiated the notion that the GmNAC genes display functional and expression profiles consistent with their phylogenetic relatedness and established a complete framework of the soybean NAC superfamily as a foundation for future analyses.

Insights into soybean transcriptome reconfiguration under hypoxic stress: Functional, regulatory, structural, and compositional characterization

Soybean (Glycine max) is one of the major crops worldwide and flooding stress affects the production and expansion of cultivated areas. Oxygen is essential for mitochondrial aerobic respiration to supply the energy demand of plant cells. Because oxygen diffusion in water is 10,000 times lower than in air, partial (hypoxic) or total (anoxic) oxygen deficiency is important component of flooding. Even when oxygen is externally available, oxygen deficiency frequently occurs in bulky, dense or metabolically active tissues such as phloem, meristems, seeds, and fruits. In this study, we analyzed conserved and divergent root transcriptional responses between flood-tolerant Embrapa 45 and flood-sensitive BR 4 soybean cultivars under hypoxic stress conditions with RNA-seq. To understand how soybean genes evolve and respond to hypoxia, stable and differentially expressed genes were characterized structurally and compositionally comparing its mechanistic relationship. Between cultivars, Embrapa 45 showed less up- and more down-regulated genes, and stronger induction of phosphoglucomutase (Glyma05g34790), unknown protein related to N-terminal protein myristoylation (Glyma06g03430), protein suppressor of phyA-105 (Glyma06g37080), and fibrillin (Glyma10g32620). RNA-seq and qRT-PCR analysis of non-symbiotic hemoglobin (Glyma11g12980) indicated divergence in gene structure between cultivars. Transcriptional changes for genes in amino acids and derivative metabolic process suggest involvement of amino acids metabolism in tRNA modifications, translation accuracy/efficiency, and endoplasmic reticulum stress in both cultivars under hypoxia. Gene groups differed in promoter TATA box, ABREs (ABA-responsive elements), and CRT/DREs (C-repeat/dehydration-responsive elements) frequency. Gene groups also differed in structure, composition, and codon usage, indicating biological significances. Additional data suggests that cis-acting ABRE elements can mediate gene expression independent of ABA in soybean roots under hypoxia.

Characterization of BiP Genes from Pepper (Capsicum annuum L.) and the Role of CaBiP1 in Response to Endoplasmic Reticulum and Multiple Abiotic Stresses

Adverse environmental conditions have a detrimental impact on crop growth and development, and cause protein denaturation or misfolding. The binding protein (BiP) plays an important protective role by alleviating endoplasmic reticulum (ER) stress induced by misfolded proteins. In this study, we characterized three BiP genes (CaBiP1, CaBiP2, and CaBiP3) in pepper, an economically important vegetable and spice species. The role of CaBiP1 in plant tolerance to ER stress and adverse environmental conditions (including heat, salinity, osmotic and drought stress) were investigated. All the expected functional and signaling domains were detected in three BiP proteins, but the motifs and exon-intron distribution differed slightly in CaBiP3. CaBiP1 and CaBiP2 were constitutively expressed in all the tested tissues under both normal and stressed conditions, whereas CaBiP3 was mainly expressed following stress. Silencing of CaBiP1 reduced pepper tolerance to ER stress and various environment stresses, and was accompanied by increased H₂O₂ accumulation, MDA content, relative electric leakage (REL), water loss rate, and a reduction in soluble protein content and relative water content (RWC) in the leaves. Conversely, overexpression of CaBiP1 in Arabidopsis enhanced tolerance to ER stress and multiple environment stresses, as demonstrated by an increase in germination rate, root length, survival rate, RWC, the unfolded protein response (UPR) pathway, and a decrease in water loss rate. Our results suggest that CaBiP1 may contribute to plant tolerance to abiotic stresses by reducing ROS accumulation, increasing the water-retention ability, and stimulating UPR pathways and expression of stress-related genes.

Functional and regulatory conservation of the soybean ER stress-induced DCD/NRP-mediated cell death signaling in plants

BACKGROUND

The developmental and cell death domain (DCD)-containing asparagine-rich proteins (NRPs) were first identified in soybean (Glycine max) as transducers of a cell death signal derived from prolonged endoplasmic reticulum (ER) stress, osmotic stress, drought or developmentally-programmed leaf senescence via the GmNAC81/GmNAC30/GmVPE signaling module. In spite of the relevance of the DCD/NRP-mediated signaling as a versatile adaptive response to multiple stresses, mechanistic knowledge of the pathway is lacking and the extent to which this pathway may operate in the plant kingdom has not been investigated.

RESULTS

Here, we demonstrated that the DCD/NRP-mediated signaling also propagates a stress-induced cell death signal in other plant species with features of a programmed cell death (PCD) response. In silico analysis revealed that several plant genomes harbor conserved sequences of the pathway components, which share functional analogy with their soybean counterparts. We showed that GmNRPs, GmNAC81and VPE orthologs from Arabidopsis, designated as AtNRP-1, AtNRP-2, ANAC036 and gVPE, respectively, induced cell death when transiently expressed in N. benthamiana leaves. In addition, loss of AtNRP1 and AtNRP2 function attenuated ER stress-induced cell death in Arabidopsis, which was in marked contrast with the enhanced cell death phenotype displayed by overexpressing lines as compared to Col-0. Furthermore, atnrp-1 knockout mutants displayed enhanced sensitivity to PEG-induced osmotic stress, a phenotype that could be complemented with ectopic expression of either GmNRP-A or GmNRP-B. In addition, AtNRPs, ANAC036 and gVPE were induced by osmotic and ER stress to an extent that was modulated by the ER-resident molecular chaperone binding protein (BiP) similarly as in soybean. Finally, as putative downstream components of the NRP-mediated cell death signaling, the stress induction of AtNRP2, ANAC036 and gVPE was dependent on the AtNRP1 function. BiP overexpression also conferred tolerance to water stress in Arabidopsis, most likely due to modulation of the drought-induced NRP-mediated cell death response.

CONCLUSION

Our results indicated that the NRP-mediated cell death signaling operates in the plant kingdom with conserved regulatory mechanisms and hence may be target for engineering stress tolerance and adaptation in crops.

Deciduous forest responses to temperature, precipitation, and drought imply complex climate change impacts

Changes in spring and autumn phenology of temperate plants in recent decades have become iconic bio-indicators of rapid climate change. These changes have substantial ecological and economic impacts. However, autumn phenology remains surprisingly little studied. Although the effects of unfavorable environmental conditions (e.g., frost, heat, wetness, and drought) on autumn phenology have been observed for over 60 y, how these factors interact to influence autumn phenological events remain poorly understood. Using remotely sensed phenology data from 2001 to 2012, this study identified and quantified significant effects of a suite of environmental factors on the timing of fall dormancy of deciduous forest communities in New England, United States. Cold, frost, and wet conditions, and high heat-stress tended to induce earlier dormancy of deciduous forests, whereas moderate heat- and drought-stress delayed dormancy. Deciduous forests in two eco-regions showed contrasting, nonlinear responses to variation in these explanatory factors. Based on future climate projection over two periods (2041-2050 and 2090-2099), later dormancy dates were predicted in northern areas. However, in coastal areas earlier dormancy dates were predicted. Our models suggest that besides warming in climate change, changes in frost and moisture conditions as well as extreme weather events (e.g., drought- and heat-stress, and flooding), should also be considered in future predictions of autumn phenology in temperate deciduous forests. This study improves our understanding of how multiple environmental variables interact to affect autumn phenology in temperate deciduous forest ecosystems, and points the way to building more mechanistic and predictive models.

Comprehensive analysis of the endoplasmic reticulum stress response in the soybean genome: conserved and plant-specific features

BACKGROUND

Despite the relevance of the eukaryotic endoplasmic reticulum (ER)-stress response as an integrator of multiple stress signals into an adaptive response, knowledge about these ER-mediated cytoprotective pathways in soybean (Glycine max) is lacking. Here, we searched for genes involved in the highly conserved unfolded protein response (UPR) and ER stress-induced plant-specific cell death signaling pathways in the soybean genome.

METHODS

Previously characterized Arabidopsis UPR genes were used as prototypes for the identification of the soybean orthologs and the in silico assembly of the UPR in soybean, using eggNOG v4.0 software. Functional studies were also conducted by analyzing the transcriptional activity of soybean UPR transducers.

RESULTS

As a result of this search, we have provided a complete profile of soybean UPR genes with significant predicted protein similarities to A. thaliana UPR-associated proteins. Both arms of the plant UPR were further examined functionally, and evidence is presented that the soybean counterparts are true orthologs of previously characterized UPR transducers in Arabidopsis. The bZIP17/bZI28 orthologs (GmbZIP37 and GmbZIP38) and ZIP60 ortholog (GmbZIP68) from soybean have similar structural organizations as their Arabidopsis counterparts, were induced by ER stress and activated an ERSE- and UPRE-containing BiP promoter. Furthermore, the transcript of the putative substrate of GmIREs, GmbZIP68, harbors a canonical site for IRE1 endonuclease activity and was efficiently spliced under ER stress conditions. In a reverse approach, we also examined the Arabidopsis genome for components of a previously characterized ER stress-induced cell death signaling response in soybean. With the exception of GmERD15, which apparently does not possess an Arabidopsis ortholog, the Arabidopsis genome harbors conserved GmNRP, GmNAC81, GmNAC30 and GmVPE sequences that share significant structural and sequence similarities with their soybean counterparts. These results suggest that the NRP/GmNAC81 + GmNAC30/VPE regulatory circuit may transduce cell death signals in plant species other than soybean.

CONCLUSIONS

Our in silico analyses, along with current and previous functional data, permitted generation of a comprehensive overview of the ER stress response in soybean as a framework for functional prediction of ER stress signaling components and their possible connections with multiple stress responses.