Comprehensive characterization and RNA-Seq profiling of the HD-Zip transcription factor family in soybean (Glycine max) during dehydration and salt stress

Soybean PHR1-regulated low phosphorus-responsive GmRALF22 promotes phosphate uptake by stimulating the expression of GmPTs

Phosphorus (P) is an essential macronutrient for plant growth and development. Rapid alkalisation factors (RALFs) play crucial roles in plant responses to nutrient stress. However, the functions of Glycine max RALFs (GmRALFs) under low P (LP) stress remain elusive. In this study, we first identified 27 GmRALFs in soybean and then revealed that, under LP conditions, GmRALF10, GmRALF11, and GmRALF22 were induced in both roots and leaves, whereas GmRALF5, GmRALF6, and GmRALF25 were upregulated in leaves. Furthermore, GmRALF22 was found to be the target gene of the transcription factor GmPHR1, which binds to the P1BS cis-element in the promoter of GmRALF22 via electrophoretic mobility shift assay and dual-luciferase experiments. Colonisation with Bacillus subtilis which delivers GmRALF22, increases the expression of the high-affinity phosphate (Pi) transporter genes GmPT2, GmPT11, GmPT13, and GmPT14, thus increasing the total amount of dry matter and soluble Pi in soybeans. RNA sequencing revealed that GmRALF22 alleviates LP stress by regulating the expression of jasmonic acid- (JA-), salicylic acid- (SA-), and immune-related genes. Finally, we verified that GmRALF22 was dependent on FERONIA (FER) to promote Arabidopsis primary root growth under LP conditions. In summary, the GmPHR1-GmRALF22 module positively regulates soybean tolerance to LP.

FaTEDT1L of Octoploid Cultivated Strawberry Functions as a Transcriptional Activator and Enhances Abiotic Stress Tolerance in Transgenic Arabidopsis

Plants may encounter abiotic stresses, such as drought, flooding, salinity, and extreme temperatures, thereby negatively affecting their growth, development, and reproduction. In order to enhance their tolerance to such stresses, plants have developed intricate signaling networks that regulate stress-responsive gene expression. For example, Arabidopsis Enhanced Drought Tolerance1/HOMEODOMAIN GLABROUS 11 (AtEDT1/HDG11), one of the transcription factor genes from the group IV of homeodomain-leucine zipper (HD-ZIP) gene family, has been shown to increase drought tolerance in various transgenic plants. However, the underlying molecular mechanisms of enhanced stress tolerance remain unclear. In this study, we identified a homologous gene related to AtEDT1/HDG11, named FaTEDT1L, from the transcriptome sequencing database of cultivated strawberry. Phylogenetic analysis revealed the close relationship of FaTEDT1L with AtEDT1/HDG11, which is one of the group IV members of the HD-ZIP gene family. Yeast one-hybrid analysis showed that FaTEDT1L functions as a transcriptional activator. Transgenic Arabidopsis plants overexpressing FaTEDT1L under the control of the cauliflower mosaic virus (CaMV) 35S promoter exhibited significantly enhanced tolerance to osmotic stress (both drought and salinity) when compared to the wild-type (WT) plants. Under osmotic stress, the average root length was 3.63 ± 0.83 cm, 4.20 ± 1.03 cm, and 4.60 ± 1.14 cm for WT, 35S::FaTEDT1L T2 #3, and 35S:: FaTEDT1L T2 #5, respectively. Substantially increased root length in 35S::FaTEDT1L T2 #3 and 35S::FaTEDT1L T2 #5 was noted when compared to the WT. In addition, the average water loss rates were 64%, 57.1%, and 55.6% for WT, 35S::FaTEDT1L T2 #3, and 35S::FaTEDT1L T2 #5, respectively, after drought treatment, indicating a significant decrease in water loss rate of 35S:: FaTEDT1L T2 #3 and 35S::FaTEDT1L T2 #5 is a critical factor in enhancing plant drought resistance. These findings thus highlight the crucial role of FaTEDT1L in mitigating drought and salt stresses and regulating plant osmotic stress tolerance. Altogether, FaTEDT1L shows its potential usage as a candidate gene for strawberry breeding in improving crop resilience and increasing agricultural productivity under adverse environmental conditions.

Comprehensive identification and expression analyses of sugar transporter genes reveal the role of GmSTP22 in salt stress resistance in soybean

The transport, compartmentation and allocation of sugar are critical for plant growth and development, as well as for stress resistance, but sugar transporter genes have not been comprehensively characterized in soybean. Here, we performed a genome-wide identification and expression analyses of sugar transporter genes in soybean in order to reveal their putative functions. A total of 122 genes encoding sucrose transporters (SUTs) and monosaccharide transporters (MSTs) were identified in soybean. They were classified into 8 subfamilies according to their phylogenetic relationships and their conserved motifs. Comparative genomics analysis indicated that whole genome duplication/segmental duplication and tandem duplication contributed to the expansion of sugar transporter genes in soybean. Expression analysis by retrieving transcriptome datasets suggested that most of these sugar transporter genes were expressed in various tissues, and a number of genes exhibited tissue-specific expression patterns. Several genes including GmSTP21, GmSFP8, and GmPLT5/6/7/8/9 were predominantly expressed in nodules, and GmPLT8 was significantly induced by rhizobia inoculation in root hairs. Transcript profiling and qRT-PCR analyses suggested that half of these sugar transporter genes were significantly induced or repressed under stresses like salt, drought, and cold. In addition, GmSTP22 was found to be localized in the plasma membrane, and its overexpression promoted plant growth and salt tolerance in transgenic Arabidopsis under the supplement with glucose or sucrose. This study provides insights into the evolutionary expansion, expression pattern and functional divergence of sugar transporter gene family, and will enable further understanding of their biological functions in the regulation of growth, yield formation and stress resistance of soybean.

Identification of GmPT proteins and investigation of their expressions in response to abiotic stress in soybean

MAIN CONCLUSION

Investigation the expression patterns of GmPT genes in response to various abiotic stresses and overexpression of GmPT11 in soybean hairy roots and Arabidopsis exhibited hypersensitivity to salt stress. Soybean is considered to be one of the significant oil crops globally, as it offers a diverse range of essential nutrients that contribute to human health. Salt stress seriously affects the yield of soybean through negative impacts on the growth, nodulation, reproduction, and other agronomy traits. The phosphate transporters 1(PHT1) subfamily, which is a part of the PHTs family in plants, is primarily found in the cell membrane and responsible for the uptake and transport of phosphorus. However, the role of GmPT (GmPT1-GmPT14) genes in response to salt stress has not been comprehensively studied. Here, we conducted a systematic analysis to ascertain the distribution and genomic duplications of GmPT genes, as well as their expression patterns in response to various abiotic stresses. Promoter analysis of GmPT genes revealed that six stress-related cis-elements were enriched in these genes. The overexpression of GmPT11 in soybean hairy roots and Arabidopsis exhibited hypersensitivity to salt stress, while no significant change was observed under low phosphate treatment, suggesting a crucial role in the response to salt stress. These findings provide novel insights into enhancing plant tolerance to salt stress.

Whole-transcriptome characterization and functional analysis of lncRNA-miRNA-mRNA regulatory networks responsive to sugarcane mosaic virus in maize resistant and susceptible inbred lines

Sugarcane mosaic virus (SCMV) is one of the most important pathogens causing maize dwarf mosaic disease, which seriously affects the yield and quality of maize. Currently, the molecular mechanism of non-coding RNAs (ncRNAs) responding to SCMV infection in maize is still uncovered. In this study, a total of 112 differentially expressed (DE)-long non-coding RNAs (lncRNAs), 24 DE-microRNAs (miRNAs), and 1822 DE-messenger RNAs (mRNAs), and 363 DE-lncRNAs, 230 DE-miRNAs, and 4376 DE-mRNAs were identified in maize resistant (Chang7-2) and susceptible (Mo17) inbred lines in response to SCMV infection through whole-transcriptome RNA sequencing, respectively. Moreover, 4874 mRNAs potentially targeted by 635 miRNAs were obtained by degradome sequencing. Subsequently, several crucial SCMV-responsive lncRNA-miRNA-mRNA networks were established, of which the expression levels of lncRNA10865-miR166j-3p-HDZ25/69 (class III homeodomain-leucine zipper 25/69) module, and lncRNA14234-miR394a-5p-SPL11 (squamosal promoter-binding protein-like 11) module were further verified. Additionally, silencing lncRNA10865 increased the accumulations of SCMV and miR166j-3p, while silencing lncRNA14234 decreased the accumulations of SCMV and SPL11 targeted by miR394a-5p. This study revealed the interactions of lncRNAs, miRNAs and mRNAs in maize resistant and susceptible materials, providing novel clues to reveal the mechanism of maize in resistance to SCMV from the perspective of competing endogenous RNA (ceRNA) regulatory networks.

Genome-wide analysis of the SCAMPs gene family of soybean and functional identification of GmSCAMP5 in salt tolerance

The effect of salt damage on plants is mainly caused by the toxic effect of Na+. Studies showed that the secretory carrier membrane proteins were associated with the Na+ transport. However, the salt tolerance mechanism of secretory carrier protein (SCAMP) in soybean was yet to be defined. In this study, ten potential SCAMP genes distributed in seven soybean chromosomes were identified in the soybean genome. The phylogenetic tree of SCAMP domain sequences of several plants can divide SCAMPs into two groups. Most GmSCAMPs genes contained multiple Box4, MYB and MYC cis-elements indicated they may respond to abiotic stresses. We found that GmSCAMP1, GmSCAMP2 and GmSCAMP4 expressed in several tissues and GmSCAMP5 was significantly induced by salt stress. GmSCAMP5 showed the same expression patterns under NaCl treatment in salt-tolerant and salt-sensitive soybean varieties, but the induced time of GmSCAMP5 in salt-tolerant variety was earlier than that of salt-sensitive variety. To further study the effect of GmSCAMP5 on the salt tolerance of soybean plants, compared to GmSCAMP5-RNAi and EV-Control plants, GmSCAMP5-OE had less wilted leave and higher SPAD value. Compared to empty vector control, less trypan blue staining was observed in GmSCAMP5-OE leaves while more staining in GmSCAMP5-RNAi leaves. The Na+ of GmSCAMP5-RNAi plants leaves under NaCl stress were significantly higher than that in EV-Control plants, while significantly lower Na+ in GmSCAMP5-OE plants than in that EV-Control plants. The contents of leaves K+ of GmSCAMP5-RNAi, EV-Control, and GmSCAMP5-OE plants under NaCl stress were opposite to that of leaves Na+ content. Finally, salt stress-related genes NHX1, CLC1, TIP1, SOD1, and SOS1 in transformed hairy root changed significantly compared with the empty control. The research will provide novel information for study the molecular regulation mechanism of soybean salt tolerance.

Genome-Wide Analysis of Cation/Proton Antiporter Family in Soybean (Glycine max) and Functional Analysis of GmCHX20a on Salt Response

Monovalent cation proton antiporters (CPAs) play crucial roles in ion and pH homeostasis, which is essential for plant development and environmental adaptation, including salt tolerance. Here, 68 CPA genes were identified in soybean, phylogenetically dividing into 11 Na+/H+ exchangers (NHXs), 12 K+ efflux antiporters (KEAs), and 45 cation/H+ exchangers (CHXs). The GmCPA genes are unevenly distributed across the 20 chromosomes and might expand largely due to segmental duplication in soybean. The GmCPA family underwent purifying selection rather than neutral or positive selections. The cis-element analysis and the publicly available transcriptome data indicated that GmCPAs are involved in development and various environmental adaptations, especially for salt tolerance. Based on the RNA-seq data, twelve of the chosen GmCPA genes were confirmed for their differentially expression under salt or osmotic stresses using qRT-PCR. Among them, GmCHX20a was selected due to its high induction under salt stress for the exploration of its biological function on salt responses by ectopic expressing in Arabidopsis. The results suggest that the overexpression of GmCHX20a increases the sensitivity to salt stress by altering the redox system. Overall, this study provides comprehensive insights into the CPA family in soybean and has the potential to supply new candidate genes to develop salt-tolerant soybean varieties.

Comprehensive and evolutionary analysis of Spodoptera litura-inducible Cytochrome P450 monooxygenase gene family in Glycine max elucidate their role in defense

Plants being sessile organisms and lacking both circulating phagocytic cells and somatic adaptive immune response, have thrived on various defense mechanisms to fend off insect pests and invasion of pathogens. CYP450s are the versatile enzymes, which thwart plants against insect pests by ubiquitous biosynthesis of phytohormones, antioxidants, and secondary metabolites, utilizing them as feeding deterrents and direct toxins. Therefore, a comprehensive analysis of biotic stress-responsive CYPs from Glycine max was performed to ascertain their function against S. litura-infestation. Phylogenetic analysis and evolutionary studies on conserved domains and motifs disclosed the evolutionary correspondence of these GmCYPs with already characterized members of the CYP450 superfamily and close relatedness to Medicago truncatula. These GmCYPs were mapped on 13 chromosomes; they possess 1-8 exons; they have evolved due to duplication and are localized in endoplasmic reticulumn. Further, identification of methyl-jasmonate, salicylic acid, defense responsive and flavonoid biosynthesis regulating cis-acting elements, their interaction with biotic stress regulating proteins and their differential expression in diverse types of tissues, and during herbivory, depicted their responsiveness to biotic stress. Three-dimensional homology modelling of GmCYPs, docking with heme cofactor required for their catalytic activity and enzyme-substrate interactions were performed to understand the functional mechanism of their action. Moreover, to gain insight into their involvement in plant defense, gene expression analysis was evaluated, which revealed differential expression of 11 GmCYPs upon S. litura-infestation, 12 GmCYPs on wounding while foliar spray of ethylene, methyl-jasmonate and salicylic acid differentially regulated 11 GmCYPs, 6 GmCYPs, and 10 GmCYPs respectively. Our study comprehensively analysed the underlying mechanism of GmCYPs function during S. litura-infestation, which can be further utilized for functional characterization to develop new strategies for enhancing soybean resistance to insect pests.

Twenty years of mining salt tolerance genes in soybean

Current combined challenges of rising food demand, climate change and farmland degradation exert enormous pressure on agricultural production. Worldwide soil salinization, in particular, necessitates the development of salt-tolerant crops. Soybean, being a globally important produce, has its genetic resources increasingly examined to facilitate crop improvement based on functional genomics. In response to the multifaceted physiological challenge that salt stress imposes, soybean has evolved an array of defences against salinity. These include maintaining cell homeostasis by ion transportation, osmoregulation, and restoring oxidative balance. Other adaptations include cell wall alterations, transcriptomic reprogramming, and efficient signal transduction for detecting and responding to salt stress. Here, we reviewed functionally verified genes that underly different salt tolerance mechanisms employed by soybean in the past two decades, and discussed the strategy in selecting salt tolerance genes for crop improvement. Future studies could adopt an integrated multi-omic approach in characterizing soybean salt tolerance adaptations and put our existing knowledge into practice via omic-assisted breeding and gene editing. This review serves as a guide and inspiration for crop developers in enhancing soybean tolerance against abiotic stresses, thereby fulfilling the role of science in solving real-life problems.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01383-3.

Integrated mRNA and miRNA transcriptome analysis of grape in responses to salt stress

Salt stress is an important factor which may negatively affect plant growth and development. High concentrations of Na⁺ ions can destroy the ion balance in plant somatic cells, as well as destroying cell membranes and forming a large number of reactive oxygen species (ROS) and other damage mechanisms. However, plants have evolved numerous defense mechanisms in response to the damages caused by salt stress conditions. Grape (Vitis vinifera L.), a type of economic crop, is widely planted throughout the world. It has been found that salt stress is an important factor affecting the quality and growth of grape crops. In this study, a high-throughput sequencing method was used to identify the differentially expressed miRNAs and mRNAs in grapes as responses to salt stress. A total of 7,856 differentially expressed genes under the salt stress conditions were successfully identified, of which 3,504 genes were observed to have up-regulated expressions and 4,352 genes had down-regulated expressions. In addition, this study also identified 3,027 miRNAs from the sequencing data using bowtie and mireap software. Among those, 174 were found to be highly conserved, and the remaining miRNAs were less conserved. In order to analyze the expression levels of those miRNAs under salt stress conditions, a TPM algorithm and DESeq software were utilized to screen the differentially expressed miRNAs among different treatments. Subsequently, a total of thirty-nine differentially expressed miRNAs were identified, of which fourteen were observed to be up-regulated miRNAs and twenty-five were down-regulated under the salt stress conditions. A regulatory network was built in order to examine the responses of grape plants to salt stress, with the goal of laying a solid foundation for revealing the molecular mechanism of grape in responses to salt stress.

Identification of QTN-by-environment interactions and their candidate genes for soybean seed oil-related traits using 3VmrMLM

Introduction

Although seed oil content and its fatty acid compositions in soybean were affected by environment, QTN-by-environment (QEIs) and gene-by-environment interactions (GEIs) were rarely reported in genome-wide association studies.

Methods

The 3VmrMLM method was used to associate the trait phenotypes, measured in five to seven environments, of 286 soybean accessions with 106,013 SNPs for detecting QTNs and QEIs.

Results

Seven oil metabolism genes (GmSACPD-A, GmSACPD-B, GmbZIP123, GmSWEET39, GmFATB1A, GmDGAT2D, and GmDGAT1B) around 598 QTNs and one oil metabolism gene GmFATB2B around 54 QEIs were verified in previous studies; 76 candidate genes and 66 candidate GEIs were predicted to be associated with these traits, in which 5 genes around QEIs were verified in other species to participate in oil metabolism, and had differential expression across environments. These genes were found to be related to soybean seed oil content in haplotype analysis. In addition, most candidate GEIs were co-expressed with drought response genes in co-expression network, and three KEGG pathways which respond to drought were enriched under drought stress rather than control condition; six candidate genes were hub genes in the co-expression networks under drought stress.

Discussion

The above results indicated that GEIs, together with drought response genes in co-expression network, may respond to drought, and play important roles in regulating seed oil-related traits together with oil metabolism genes. These results provide important information for genetic basis, molecular mechanisms, and soybean breeding for seed oil-related traits.

Overexpression of Pyrus sinkiangensis HAT5 enhances drought and salt tolerance, and low-temperature sensitivity in transgenic tomato

The homeodomain-leucine zipper protein HAT belongs to the homeodomain leucine zipper subfamily (HD-Zip) and is important for regulating plant growth and development and stress tolerance. To investigate the role of HAT5 in tolerance to drought, salt, and low temperature stress, we selected a HAT gene from Pyrus sinkiangensis Yü (Pyrus sinkiangensis T.T. Yu). The sequences were analyzed using ioinformatics, and the overexpressed tomato lines were obtained using molecular biology techniques. The phenotypes, physiological, and biochemical indexes of the wild-type and transgenic tomato lines were observed under different stress conditions. We found that the gene had the highest homology with PbrHAT5. Under drought and NaCl stress, osmotic regulatory substances (especially proline) were significantly accumulated, and antioxidant enzyme activities were enhanced. The malondialdehyde level and relative electrical conductivity of transgenic tomatoes under low temperature (freezing) stress were significantly higher than those of wild-type tomatoes. The reactive oxygen species scavenging system was unbalanced. This study found that PsHAT5 improved the tolerance of tomatoes to drought and salt stress by regulating proline metabolism and oxidative stress ability, reducing the production of reactive oxygen species, and maintaining normal cell metabolism. In conclusion, the PsHAT5 transcription factor has great potential in crop resistance breeding, which lays a theoretical foundation for future excavation of effective resistance genes of the HD-Zip family and experimental field studies.

Molecular Tools and Their Applications in Developing Salt-Tolerant Soybean (Glycine max L.) Cultivars

Abiotic stresses are one of the significant threats to soybean (Glycine max L.) growth and yields worldwide. Soybean has a crucial role in the global food supply chain and food security and contributes the main protein share compared to other crops. Hence, there is a vast scientific saddle on soybean researchers to develop tolerant genotypes to meet the growing need of food for the huge population. A large portion of cultivated land is damaged by salinity stress, and the situation worsens yearly. In past years, many attempts have increased soybean resilience to salinity stress. Different molecular techniques such as quantitative trait loci mapping (QTL), genetic engineering, transcriptome, transcription factor analysis (TFs), CRISPR/Cas9, as well as other conventional methods are used for the breeding of salt-tolerant cultivars of soybean to safeguard its yield under changing environments. These powerful genetic tools ensure sustainable soybean yields, preserving genetic variability for future use. Only a few reports about a detailed overview of soybean salinity tolerance have been published. Therefore, this review focuses on a detailed overview of several molecular techniques for soybean salinity tolerance and draws a future research direction. Thus, the updated review will provide complete guidelines for researchers working on the genetic mechanism of salinity tolerance in soybean.

Systematic Characterization of the OSCA Family Members in Soybean and Validation of Their Functions in Osmotic Stress

Since we discovered OSCA1, a hyperosmolarity-gated calcium-permeable channel that acted as an osmosensor in Arabidopsis, the OSCA family has been identified genome-wide in several crops, but only a few OSCA members' functions have been experimentally demonstrated. Osmotic stress seriously restricts the yield and quality of soybean. Therefore, it is essential to decipher the molecular mechanism of how soybean responds to osmotic stress. Here, we first systematically studied and experimentally demonstrated the role of OSCA family members in the osmotic sensing of soybean. Phylogenetic relationships, gene structures, protein domains and structures analysis revealed that 20 GmOSCA members were divided into four clades, of which members in the same cluster may have more similar functions. In addition, GmOSCA members in clusters III and IV may be functionally redundant and diverged from those in clusters I and II. Based on the spatiotemporal expression patterns, GmOSCA1.6, GmOSCA2.1, GmOSCA2.6, and GmOSCA4.1 were extremely low expressed or possible pseudogenes. The remaining 16 GmOSCA genes were heterologously overexpressed in an Arabidopsis osca1 mutant, to explore their functions. Subcellular localization showed that most GmOSCA members could localize to the plasma membrane (PM). Among 16 GmOSCA genes, only overexpressing GmOSCA1.1, GmOSCA1.2, GmOSCA1.3, GmOSCA1.4, and GmOSCA1.5 in cluster I could fully complement the reduced hyperosmolality-induced [Ca²⁺]_i increase (OICI) in osca1. The expression profiles of GmOSCA genes against osmotic stress demonstrated that most GmOSCA genes, especially GmOSCA1.1, GmOSCA1.2, GmOSCA1.3, GmOSCA1.4, GmOSCA1.5, GmOSCA3.1, and GmOSCA3.2, strongly responded to osmotic stress. Moreover, overexpression of GmOSCA1.1, GmOSCA1.2, GmOSCA1.3, GmOSCA1.4, GmOSCA1.5, GmOSCA3.1, and GmOSCA3.2 rescued the drought-hypersensitive phenotype of osca1. Our findings provide important clues for further studies of GmOSCA-mediated calcium signaling in the osmotic sensing of soybean and contribute to improving soybean drought tolerance through genetic engineering and molecular breeding.

CRISPR/Cas9 mediated gene-editing of GmHdz4 transcription factor enhances drought tolerance in soybean (Glycine max [L.] Merr.)

The HD-Zip transcription factors play a crucial role in plant development, secondary metabolism, and abiotic stress responses, but little is known about HD-Zip I genes in soybean. Here, a homeodomain-leucine zipper gene designated GmHdz4 was isolated. Chimeric soybean plants, GmHdz4 overexpressing (GmHdz4-oe), and gene-editing via CRISPR/Cas9 (gmhdz4) in hairy roots, were generated to examine the GmHdz4 gene response to polyethylene glycol (PEG)-simulated drought stress. Bioinformatic analysis showed GmHdz4 belonged to clade δ, and was closely related to other drought tolerance-related HD-Zip I family genes such as AtHB12, Oshox12, and Gshdz4. The GmHdz4 was located in the plant nucleus and showed transcriptional activation activity by yeast hybrid assay. Quantitative real-time PCR analysis revealed that GmHdz4 expression varied in tissues and was induced by PEG-simulated drought stress. The gmhdz4 showed promoted growth of aboveground parts, and its root system architecture, including the total root length, the root superficial area, and the number of root tips were significantly higher than those of GmHdz4-oe even the non-transgenic line (NT) on root tips number. The better maintenance of turgor pressure by osmolyte accumulation, and the higher activity of antioxidant enzymes to scavenge reactive oxygen species, ultimately suppressed the accumulation of hydrogen peroxide (H₂O₂), superoxide anion (O^2-), and malondialdehyde (MDA), conferring higher drought tolerance in gmhdz4 compared with both GmHdz4-oe and NT. Together, our results provide new insights for future research on the mechanisms by which GmHdz4 gene-editing via CRISPR/Cas9 system could promote drought stress and provide a potential target for molecular breeding in soybean.

The Hidden Potential of High-Throughput RNA-Seq Re-Analysis, a Case Study for DHDPS, Key Enzyme of the Aspartate-Derived Lysine Biosynthesis Pathway and Its Role in Abiotic and Biotic Stress Responses in Soybean

DHDPS is a key enzyme in the aspartate-derived lysine biosynthesis pathway and an evident object of study for biofortification strategies in plants. DHDPS isoforms with novel regulatory properties in Medicago truncatula were demonstrated earlier and hypothesized to be involved in abiotic and biotic stress responses. Here, we present a phylogenetic analysis of the DHPDS gene family in land plants which establishes the existence of a legume-specific class of DHDPS, termed DHDPS B-type, distinguishable from the DHDPS A-type commonly present in all land plants. The G. max genome comprises two A-type DHDPS genes (Gm.DHDPS-A1; Glyma.09G268200, Gm.DHDPS-A2; Glyma.18G221700) and one B-type (Gm.DHDPS-B; Glyma.03G022300). To further investigate the expression pattern of the G. max DHDPS isozymes in different plant tissues and under various stress conditions, 461 RNA-seq experiments were exploited and re-analyzed covering two expression atlases, 13 abiotic and 5 biotic stress studies. Gm.DHDPS-B is seen almost exclusively expressed in roots and nodules in addition to old cotyledons or senescent leaves while both DHDPS A-types are expressed constitutively in all tissues analyzed with the highest expression in mature seeds. Furthermore, Gm.DHDPS-B expression is significantly upregulated in some but not all stress responses including salt stress, flooding, ethylene or infection with Phytophthora sojae and coincides with downregulation of DHDPS A-types. In conclusion, we demonstrate the potential of an in-depth RNA-seq re-analysis for the guidance of future experiments and to expand on current knowledge.

HD-Zip III Gene Family: Identification and Expression Profiles during Leaf Vein Development in Soybean

Leaf veins constitute the transport network for water and photosynthetic assimilates in vascular plants. The class III homeodomain-leucine zipper (HD-Zip III) gene family is central to the regulation of vascular development. In this research, we performed an overall analysis of the HD-Zip III genes in soybean (Glycine max L. Merr.). Our analysis included the phylogeny, conservation domains and cis-elements in the promoters of these genes. We used the quantitative reverse transcription-polymerase chain reaction to characterize the expression patterns of HD-Zip III genes in leaf vein development and analyze the effects of exogenous hormone treatments. In this study, twelve HD-Zip III genes were identified from the soybean genome and named. All soybean HD-Zip III proteins contained four highly conserved domains. GmHB15-L-1 transcripts showed steadily increasing accumulation during all stages of leaf vein development and were highly expressed in cambium cells. GmREV-L-1 and GmHB14-L-2 had nearly identical expression patterns in soybean leaf vein tissues. GmREV-L-1 and GmHB14-L-2 transcripts remained at stable high levels at all xylem developmental stages. GmREV-L-1 and GmHB14-L-2 were expressed at high levels in the vascular cambium and xylem cells. Overall, GmHB15-L-1 may be an essential regulator that is responsible for the formation or maintenance of soybean vein cambial cells. GmREV-L-1 and GmHB14-L-2 were correlated with xylem differentiation in soybean leaf veins. This study will pave the way for identifying the molecular mechanism of leaf vein development.

Genome-Wide Identification and Expression Analysis of Homeodomain Leucine Zipper Subfamily IV (HD-ZIP IV) Gene Family in Cannabis sativa L

The plant-specific homeodomain zipper family (HD-ZIP) of transcription factors plays central roles in regulating plant development and environmental resistance. HD-ZIP transcription factors IV (HDZ IV) have been involved primarily in the regulation of epidermal structure development, such as stomata and trichomes. In our study, we identified nine HDZ IV-encoding genes in Cannabis sativa L. by conducting a computational analysis of cannabis genome resources. Our analysis suggests that these genes putatively encode proteins that have all the conserved domains of HDZ IV transcription factors. The phylogenetic analysis of HDZ IV gene family members of cannabis, rice (Oryza sativa), and Arabidopsis further implies that they might have followed distinct evolutionary paths after divergence from a common ancestor. All the identified cannabis HDZ IV gene promoter sequences have multiple regulation motifs, such as light- and hormone-responsive elements. Furthermore, experimental evidence shows that different HDZ IV genes have different expression patterns in root, stem, leaf, and flower tissues. Four genes were primarily expressed in flowers, and the expression of CsHDG5 (XP_030501222.1) was also correlated with flower maturity. Fifty-nine genes were predicted as targets of HDZ IV transcription factors. Some of these genes play central roles in pathogen response, flower development, and brassinosteroid signaling. A subcellular localization assay indicated that one gene of this family is localized in the Arabidopsis protoplast nucleus. Taken together, our work lays fundamental groundwork to illuminate the function of cannabis HDZ IV genes and their possible future uses in increasing cannabis trichome morphogenesis and secondary metabolite production.

Metabolic adjustment and regulation of gene expression are essential for increased resistance to severe water deficit and resilience post-stress in soybean

Background

Soybean is the main oilseed crop grown in the world; however, drought stress affects its growth and physiology, reducing its yield. The objective of this study was to characterize the physiological, metabolic, and genetic aspects that determine differential resistance to water deficit in soybean genotypes.

Methods

Three soybean genotypes were used in this study, two lineages (L11644 and L13241), and one cultivar (EMBRAPA 48-C48). Plants were grown in pots containing 8 kg of a mixture of soil and sand (2:1) in a greenhouse under sunlight. Soil moisture in the pots was maintained at field capacity until the plants reached the stage of development V4 (third fully expanded leaf). At this time, plants were subjected to three water treatments: Well-Watered (WW) (plants kept under daily irrigation); Water Deficit (WD) (withholding irrigation until plants reached the leaf water potential at predawn of -1.5 ± 0.2 MPa); Rewatered (RW) (plants rehydrated for three days after reached the water deficit). The WW and WD water treatments were evaluated on the eighth day for genotypes L11644 and C48, and on the tenth day for L13241, after interruption of irrigation. For the three genotypes, the treatment RW was evaluated after three days of resumption of irrigation. Physiological, metabolic and gene expression analyses were performed.

Results

Water deficit inhibited growth and gas exchange in all genotypes. The accumulation of osmolytes and the concentrations of chlorophylls and abscisic acid (ABA) were higher in L13241 under stress. The metabolic adjustment of lineages in response to WD occurred in order to accumulate amino acids, carbohydrates, and polyamines in leaves. The expression of genes involved in drought resistance responses was more strongly induced in L13241. In general, rehydration provided recovery of plants to similar conditions of control treatment. Although the C48 and L11644 genotypes have shown some tolerance and resilience responses to severe water deficit, greater efficiency was observed in the L13241 genotype through adjustments in morphological, physiological, genetic and metabolic characteristics that are combined in the same plant. This study contributes to the advancement in the knowledge about the resistance to drought in cultivated plants and provides bases for the genetic improvement of the soybean culture.

Systematic analysis of HD-ZIP transcription factors in sesame genome and gene expression profiling of SiHD-ZIP class I entailing drought stress responses at early seedling stage

BACKGROUND

Sesame is an ancient oilseed crop, known for its high oil content and quality. Its sensitivity to drought at early seedling stage is one of the limiting factors affecting its world-wide growth and productivity. Among plant specific transcription factors, the association of HD-ZIPs with sesame drought responses at early seedling stage is not well-established yet and is very important to develop our molecular understanding on sesame drought tolerance.

METHODS AND RESULTS

In this study, total 61 sesame HD-ZIP proteins were identified, based on their protein sequence homology with Arabidopsis and protein domain(s) architecture prediction, followed by their phylogenetic, conserved domain(s) motifs and gene structure analyses to classify them into four classes (HD-ZIP Class I-IV). HD-ZIP Class I was also subdivided into four subgroups: α (SiHZ25, SiHZ43, SiHZ9 and SiHZ16), β1 (SiHZ10, SiHZ30, SiHZ32 and SiHZ26), β2 (SiHZ42 and SiHZ45) and γ (SiHZ17, SiHZ7 and SiHZ35) by a comparative phylogenetic analysis of sesame with Arabidopsis and maize. Afterwards, twenty-one days old sesame seedlings were exposed to drought stress by withholding water for 7 days (when soil moisture content reduced to ~16%) and gene expression of HD-ZIP Class I (13 members) was performed in well- watered (control) and drought stressed seedlings. The gene expression analysis showed that the expressions of SiHZ7 (6.8 fold) and SiHZ35 (2.6 fold) from γ subgroup were significantly high in drought seedlings.

CONCLUSIONS

This study is useful in demonstrating the role of SiHD-ZIP Class I in sesame drought responses at early seedling stage and to develop its novel drought tolerant varieties.

Identification and characterization of stress responsive homeodomain leucine zipper transcription factors in Medicago truncatula

BACKGROUND

Homeodomain leucine zipper (HD-ZIP) transcription factors play roles in regulating plant development and responses to abiotic stresses; however, how HD-ZIP genes in Medicago truncatula are involved in abiotic stress response remains elusive.

METHODS AND RESULTS

The HD-ZIP I genes in Medicago truncatula were identified and characterized, and their expression patterns in different tissues and under different abiotic stresses were analyzed. A total of 15 Medicago truncatula HD-ZIP I genes were identified and a phylogenetic analysis of HD-ZIP I proteins in Arabidopsis thaliana and Medicago truncatula was conducted. Fifteen HD-ZIP I genes showed diverse tissue preferences. Among them, expressions of MtHB22 and MtHB51 were specially detected in vegetative buds. In addition, they responded to various abiotic stresses, including salinity and osmotic stress and abscisic acid (ABA). For instance, MtHB7 and MtHB12 expression levels were found to be positively associated with salt, osmotic stress and ABA in both shoots and roots, while MtHB13 and MtHB23 were negatively associated with these stresses in Medicago truncatula.

CONCLUSION

The HD-ZIP I genes in Medicago truncatula are evolutionarily conserved, but also exhibit gene duplication and gene loss events. Differential expression analysis of Medicago truncatula HD-ZIP I genes indicated their crucial roles in abiotic stress responses. Our genome-wide analysis of the HD-ZIP I transcription factor family in Medicago truncatula provided a valuable reference for further research.

Genome-wide identification of HD-ZIP transcription factors in maize and their regulatory roles in promoting drought tolerance

Drought is the main limiting factor of maize productivity, therefore improving drought tolerance in maize has potential practical importance. Cloning and functional verification of drought-tolerant genes is of great importance to understand molecular mechanisms under drought stress. Here, we employed a bioinformatic pipeline to identify 42 ZmHDZ drought responsive genes using previously reported maize transcriptomic datasets. The coding sequences, exon-intron structure and domain organization of all the 42 genes were identified. Phylogenetic analysis revealed evolutionary conservation of members of the ZmHDZ genes in maize. Several regulatory elements associated with drought tolerance were identified in the promoter regions of ZmHDZ genes, indicating the implication of these genes in plant response to drought stress. 42 ZmHDZ genes were distributed unevenly on 10 chromosomes, and 24 pairs of gene duplications were the segmental duplication. The expression of several ZmHDZ genes was upregulated under drought stress, and ZmHDZ9 overexpressing transgenic plants exhibited higher SOD and POD activities and higher accumulation of soluble proteins under drought stress which resulted in enhanced developed phenotype and improved resistance. The present study provides evidence for the evolutionary conservation of HD-ZIP transcription factors homologs in maize. The results further provide a comprehensive insight into the roles of ZmHDZ genes in regulating drought stress tolerance in maize.

Nearby transposable elements impact plant stress gene regulatory networks: a meta-analysis in A. thaliana and S. lycopersicum

BACKGROUND

Transposable elements (TE) make up a large portion of many plant genomes and are playing innovative roles in genome evolution. Several TEs can contribute to gene regulation by influencing expression of nearby genes as stress-responsive regulatory motifs. To delineate TE-mediated plant stress regulatory networks, we took a 2-step computational approach consisting of identifying TEs in the proximity of stress-responsive genes, followed by searching for cis-regulatory motifs in these TE sequences and linking them to known regulatory factors. Through a systematic meta-analysis of RNA-seq expression profiles and genome annotations, we investigated the relation between the presence of TE superfamilies upstream, downstream or within introns of nearby genes and the differential expression of these genes in various stress conditions in the TE-poor Arabidopsis thaliana and the TE-rich Solanum lycopersicum.

RESULTS

We found that stress conditions frequently expressed genes having members of various TE superfamilies in their genomic proximity, such as SINE upon proteotoxic stress and Copia and Gypsy upon heat stress in A. thaliana, and EPRV and hAT upon infection, and Harbinger, LINE and Retrotransposon upon light stress in S. lycopersicum. These stress-specific gene-proximal TEs were mostly located within introns and more detected near upregulated than downregulated genes. Similar stress conditions were often related to the same TE superfamily. Additionally, we detected both novel and known motifs in the sequences of those TEs pointing to regulatory cooption of these TEs upon stress. Next, we constructed the regulatory network of TFs that act through binding these TEs to their target genes upon stress and discovered TE-mediated regulons targeted by TFs such as BRB/BPC, HD, HSF, GATA, NAC, DREB/CBF and MYB factors in Arabidopsis and AP2/ERF/B3, NAC, NF-Y, MYB, CXC and HD factors in tomato.

CONCLUSIONS

Overall, we map TE-mediated plant stress regulatory networks using numerous stress expression profile studies for two contrasting plant species to study the regulatory role TEs play in the response to stress. As TE-mediated gene regulation allows plants to adapt more rapidly to new environmental conditions, this study contributes to the future development of climate-resilient plants.

Overexpression of cassava RSZ21b enhances drought tolerance in Arabidopsis

Drought is one of the major environmental constraints affecting crop productivity. Plants have to adjust their developmental and physiological processes to cope with drought. We previously identified 18 cassava serine/arginine-rich (SR) proteins that had a pivotal role in alternative splicing in response to environmental stress. However, functional characterization of SR proteins is rarely explored. Here, we characterized the RSZ subfamily gene MeRSZ21b in cassava. The RSZ21b belongs to the RSZ subfamily, which was widely distributed in major crops and was highly conserved. Quantitative RT-PCR assay showed that the expression of MeRSZ21b was significantly induced by drought. Moreover, overexpression of MeRSZ21b in Arabidopsis was hypersensitive to abscisic acid (ABA) in the phases of seed germination and post-germination seedling growth. Meantime, MeRSZ21b overexpression lines were resistant to sorbitol treatment, and quickly closed the stomata when compared with Col-0 under drought condition. Importantly, overexpression of MeRSZ21b resulted in improved drought tolerance through modulating ABA-dependent signaling. Therefore, our findings refine our knowledge of the SR protein-coding genes and provide novel insights for enhancing plant resistance to environmental stress.

Novel QTL identification and candidate gene analysis for enhancing salt tolerance in soybean (Glycine max (L.) Merr.)

Soybean, a glycophyte that is sensitive to salt stress, is greatly affected by salinity at all growth stages. A mapping population derived from a cross between a salt-sensitive Korean cultivar, Cheongja 3, and a salt-tolerant landrace, IT162669, was used to identify quantitative trait loci (QTLs) conferring salt tolerance in soybean. Following treatment with 120 mM NaCl for 2 weeks, phenotypic traits representing physiological damage, leaf Na+ content, and K+/Na+ ratio were characterized. Among the QTLs mapped on a high-density genetic map harboring 2,630 single nucleotide polymorphism markers, we found two novel major loci, qST6, on chromosome 6, and qST10, on chromosome 10, which controlled traits related to ion toxicity and physiology in response to salinity, respectively. These loci were distinct from the previously known salt tolerance allele on chromosome 3. Other QTLs associated with abiotic stress overlapped with the genomic regions of qST6 and qST10, or with their paralogous regions. Based on the functional annotation and parental expression differences, we identified eight putative candidate genes, two in qST6 and six in qST10, which included a phosphoenolpyruvate carboxylase and an ethylene response factor. This study provides additional genetic resources to breed soybean cultivars with enhanced salt tolerance.

Genome-Wide Characterization of Salt-Responsive miRNAs, circRNAs and Associated ceRNA Networks in Tomatoes

Soil salinization is a major environmental stress that causes crop yield reductions worldwide. Therefore, the cultivation of salt-tolerant crops is an effective way to sustain crop yield. Tomatoes are one of the vegetable crops that are moderately sensitive to salt stress. Global market demand for tomatoes is huge and growing. In recent years, the mechanisms of salt tolerance in tomatoes have been extensively investigated; however, the molecular mechanism through which non-coding RNAs (ncRNAs) respond to salt stress is not well understood. In this study, we utilized small RNA sequencing and whole transcriptome sequencing technology to identify salt-responsive microRNAs (miRNAs), messenger RNAs (mRNAs), and circular RNAs (circRNAs) in roots of M82 cultivated tomato and Solanum pennellii (S. pennellii) wild tomato under salt stress. Based on the theory of competitive endogenous RNA (ceRNA), we also established several salt-responsive ceRNA networks. The results showed that circRNAs could act as miRNA sponges in the regulation of target mRNAs of miRNAs, thus participating in the response to salt stress. This study provides insights into the mechanisms of salt tolerance in tomatoes and serves as an effective reference for improving the salt tolerance of salt-sensitive cultivars.

Genome-Wide Identification, Characterization and Expression Analysis of Soybean CHYR Gene Family

The CHYR (CHY ZINC-FINGER AND RING FINGER PROTEIN) proteins have been functionally characterized in iron regulation and stress response in Arabidopsis, rice and Populus. However, their roles in soybean have not yet been systematically investigated. Here, in this study, 16 GmCHYR genes with conserved Zinc_ribbon, CHY zinc finger and Ring finger domains were obtained and divided into three groups. Moreover, additional 2-3 hemerythrin domains could be found in the N terminus of Group III. Phylogenetic and homology analysis of CHYRs in green plants indicated that three groups might originate from different ancestors. Expectedly, GmCHYR genes shared similar conserved domains/motifs distribution within the same group. Gene expression analysis uncovered their special expression patterns in different soybean tissues/organs and under various abiotic stresses. Group I and II members were mainly involved in salt and alkaline stresses. The expression of Group III members was induced/repressed by dehydration, salt and alkaline stresses, indicating their diverse roles in response to abiotic stress. In conclusion, our work will benefit for further revealing the biological roles of GmCHYRs.

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.

Genome-wide identification of ascorbate-glutathione cycle gene families in soybean (Glycine max) reveals gene duplication events and specificity of gene members linked to development and stress conditions

Ascorbate-glutathione (AsA-GSH) cycle plays an important role in tuning beneficial ROS accumulation for intracellular signals and imparts plant tolerance to oxidative stress by detoxifying excess of ROS. Here, we present genome-wide identification of AsA-GSH cycle genes (APX, MDHAR, DHAR, and GR) in several leguminous species and expression analyses in G. max during stress, germination and tissue development. Our data revealed 24 genes in Glycine genus against the maximum of 15 in other leguminous species, which was due to 9 pars of duplicated genes mostly originated from sub/neofunctionalization. Cytosolic APX and MDHAR genes were highly expressed in different tissues and physiological conditions. Germination induced genes encoding AsA-GSH proteins from different cell compartments, whereas vegetative phase (leaves) stimulated predominantly genes related to chloroplast/mitochondria proteins. Moreover, cytosolic APX-1, 2, MDHAR-1a, 1b and GR genes were the primary genes linked to senescence and biotic stresses, while stAPX-a, b and GR (from organelles) were the most abiotic stress related genes. Biotic and abiotic stress tolerant genotypes generally showed increased MDHAR, DHAR and/or GR mRNA levels compared to susceptible genotypes. Overall, these data clarified evolutionary events in leguminous plants and point to the functional specificity of duplicated genes of the AsA-GSH cycle in G. max.

Genome-wide identification and analysis of WRKY gene family in maize provide insights into regulatory network in response to abiotic stresses

BACKGROUND

The WRKY transcription factor family plays significant roles in biotic and abiotic stress responses, which has been associated with various biological processes in higher plants. However, very little is known regarding the structure and function of WRKY genes in maize.

RESULTS

In this study, a total of 140 ZmWRKY proteins encoded by 125 ZmWRKY genes were eventually identified in maize. On the basis of features of molecular structure and a comparison of phylogenetic relationships of WRKY transcription factor families from Arabidopsis, rice and maize, all 140 ZmWRKY proteins in maize were divided into three main groups (Groups I, II and III) and the Group II was further classified into five subgroups. The characteristics of exon-intron structure of these putative ZmWRKY genes and conserved protein motifs of their encoded ZmWRKY proteins were also presented respectively, which was in accordance with the group classification results. Promoter analysis suggested that ZmWRKY genes shared many abiotic stress-related elements and hormone-related elements. Gene duplication analysis revealed that the segmental duplication and purifying selection might play a significant role during the evolution of the WRKY gene family in maize. Using RNA-seq data, transcriptome analysis indicated that most of ZmWRKY genes displayed differential expression patterns at different developmental stages of maize. Further, by quantitative real-time PCR analysis, twenty-one ZmWRKY genes were confirmed to respond to two different abiotic stress treatments, suggesting their potential roles in various abiotic stress responses. In addition, RNA-seq dataset was used to conduct weighted gene co-expression network analysis (WGCNA) in order to recognize gene subsets possessing similar expression patterns and highly correlated with each other within different metabolic networks. Further, subcellular localization prediction, functional annotation and interaction analysis of ZmWRKY proteins were also performed to predict their interactions and associations involved in potential regulatory network.

CONCLUSIONS

Taken together, the present study will serve to present an important theoretical basis for further exploring function and regulatory mechanism of ZmWRKY genes in the growth, development, and adaptation to abiotic stresses in maize.

Molecular Characterisation of Soybean Osmotins and Their Involvement in Drought Stress Response

Osmotins are multifunctional proteins belonging to the thaumatin-like family related to plant stress responses. To better understand the functions of soybean osmotins in drought stress response, the current study presents the characterisation of four previously described proteins and a novel putative soybean osmotin (GmOLPa-like). Gene and protein structure as well as gene expression analyses were conducted on different tissues and developmental stages of two soybean cultivars with varying dehydration sensitivities (BR16 and EMB48 are highly and slightly sensitive, respectively). The analysed osmotin sequences share the conserved amino acid signature and 3D structure of the thaumatin-like family. Some differences were observed in the conserved regions of protein sequences and in the electrostatic surface potential. P21-like present the most similar electrostatic potential to osmotins previously characterised as promoters of drought tolerance in Nicotiana tabacum and Solanum nigrum. Gene expression analysis indicated that soybean osmotins were differentially expressed in different organs (leaves and roots), developmental stages (R1 and V3), and cultivars in response to dehydration. In addition, under dehydration conditions, the highest level of gene expression was detected for GmOLPa-like and P21-like osmotins in the leaves and roots, respectively, of the less drought sensitive cultivar. Altogether, the results suggest an involvement of these genes in drought stress tolerance.

The underground life of homeodomain-leucine zipper transcription factors

Roots are the anchorage organs of plants, responsible for water and nutrient uptake, exhibiting high plasticity. Root architecture is driven by the interactions of biomolecules, including transcription factors and hormones that are crucial players regulating root plasticity. Multiple transcription factor families are involved in root development; some, such as ARFs and LBDs, have been well characterized, whereas others remain less well investigated. In this review, we synthesize the current knowledge about the involvement of the large family of homeodomain-leucine zipper (HD-Zip) transcription factors in root development. This family is divided into four subfamilies (I-IV), mainly according to structural features, such as additional motifs aside from HD-Zip, as well as their size, gene structure, and expression patterns. We explored and analyzed public databases and the scientific literature regarding HD-Zip transcription factors in Arabidopsis and other species. Most members of the four HD-Zip subfamilies are expressed in specific cell types and several individuals from each group have assigned functions in root development. Notably, a high proportion of the studied proteins are part of intricate regulation pathways involved in primary and lateral root growth and development.

Phospholipases C and D and Their Role in Biotic and Abiotic Stresses

Plants, as sessile organisms, have adapted a fine sensing system to monitor environmental changes, therefore allowing the regulation of their responses. As the interaction between plants and environmental changes begins at the surface, these changes are detected by components in the plasma membrane, where a molecule receptor generates a lipid signaling cascade via enzymes, such as phospholipases (PLs). Phospholipids are the key structural components of plasma membranes and signaling cascades. They exist in a wide range of species and in different proportions, with conversion processes that involve hydrophilic enzymes, such as phospholipase-C (PLC), phospholipase-D (PLD), and phospholipase-A (PLA). Hence, it is suggested that PLC and PLD are highly conserved, compared to their homologous genes, and have formed clusters during their adaptive history. Additionally, they generate responses to different functions in accordance with their protein structure, which should be reflected in specific signal transduction responses to environmental stress conditions, including innate immune responses. This review summarizes the phospholipid systems associated with signaling pathways and the innate immune response.

Plant Acyl-CoA-Binding Proteins-Their Lipid and Protein Interactors in Abiotic and Biotic Stresses

Plants are constantly exposed to environmental stresses during their growth and development. Owing to their immobility, plants possess stress-sensing abilities and adaptive responses to cope with the abiotic and biotic stresses caused by extreme temperatures, drought, flooding, salinity, heavy metals and pathogens. Acyl-CoA-binding proteins (ACBPs), a family of conserved proteins among prokaryotes and eukaryotes, bind to a variety of acyl-CoA esters with different affinities and play a role in the transport and maintenance of subcellular acyl-CoA pools. In plants, studies have revealed ACBP functions in development and stress responses through their interactions with lipids and protein partners. This review summarises the roles of plant ACBPs and their lipid and protein interactors in abiotic and biotic stress responses.

Genome-Wide Identification and Characterization of TALE Superfamily Genes in Soybean (Glycine max L.)

The three-amino-acid-loop-extension (TALE) superfamily genes broadly existed in plants, which played important roles in plant growth, development and abiotic stress responses. In this study, we identified 68 Glycine max TALE (GmTALE) superfamily members. Phylogenetic analysis divided the GmTALE superfamily into the BEL1-like (BLH/BELL homeodomain) and the KNOX (KNOTTED-like homeodomain) subfamilies. Moreover, the KNOX subfamily could be further categorized into three clades (KNOX Class I, KNOX Class II and KNOX Class III). The GmTALE genes showed similarities in the gene structures in the same subfamily or clade, whose coding proteins exhibited analogous motif and conserved domain compositions. Besides, synteny analyses and evolutionary constraint evaluations of the TALE members among soybean and different species provided more clues for GmTALE superfamily evolution. The cis-element analyses in gene promoter regions and relevant gene expression profiling revealed different regulating roles of GmTALE genes during soybean plant development, saline and dehydration stresses. Genome-wide characterization, evolution, and expression profile analyses of GmTALE genes can pave the way for future gene functional research and facilitate their roles for applications in genetic improvement on soybean in saline and dehydration stresses.

In silico Analysis of Acyl-CoA-Binding Protein Expression in Soybean

Plant acyl-CoA-binding proteins (ACBPs) form a highly conserved protein family that binds to acyl-CoA esters as well as other lipid and protein interactors to function in developmental and stress responses. This protein family had been extensively studied in non-leguminous species such as Arabidopsis thaliana (thale cress), Oryza sativa (rice), and Brassica napus (oilseed rape). However, the characterization of soybean (Glycine max) ACBPs, designated GmACBPs, has remained unreported although this legume is a globally important crop cultivated for its high oil and protein content, and plays a significant role in the food and chemical industries. In this study, 11 members of the GmACBP family from four classes, comprising Class I (small), Class II (ankyrin repeats), Class III (large), and Class IV (kelch motif), were identified. For each class, more than one copy occurred and their domain architecture including the acyl-CoA-binding domain was compared with Arabidopsis and rice. The expression profile, tertiary structure and subcellular localization of each GmACBP were predicted, and the similarities and differences between GmACBPs and other plant ACBPs were deduced. A potential role for some Class III GmACBPs in nodulation, not previously encountered in non-leguminous ACBPs, has emerged. Interestingly, the sole member of Class III ACBP in each of non-leguminous Arabidopsis and rice had been previously identified in plant-pathogen interactions. As plant ACBPs are known to play important roles in development and responses to abiotic and biotic stresses, the in silico expression profiles on GmACBPs, gathered from data mining of RNA-sequencing and microarray analyses, will lay the foundation for future studies in their applications in biotechnology.

Genome-wide identification and characterization of GRAS genes in soybean (Glycine max)

BACKGROUND

GRAS proteins are crucial transcription factors, which are plant-specific and participate in various plant biological processes. Thanks to the rapid progress of the whole genome sequencing technologies, the GRAS gene families in different plants have been broadly explored and studied. However, comprehensive research on the soybean (Glycine max) GRAS gene family is relatively lagging.

RESULTS

In this study, 117 Glycine max GRAS genes (GmGRAS) were identified. Further phylogenetic analyses showed that the GmGRAS genes could be categorized into nine gene subfamilies: DELLA, HAM, LAS, LISCL, PAT1, SCL3, SCL4/7, SCR and SHR. Gene structure analyses turned out that the GmGRAS genes lacked introns and were relatively conserved. Conserved domains and motif patterns of the GmGRAS members in the same subfamily or clade exhibited similarities. Notably, the expansion of the GmGRAS gene family was driven both by gene tandem and segmental duplication events. Whereas, segmental duplications took the major role in generating new GmGRAS genes. Moreover, the synteny and evolutionary constraints analyses of the GRAS proteins among soybean and distinct species (two monocots and four dicots) provided more detailed evidence for GmGRAS gene evolution. Cis-element analyses indicated that the GmGRAS genes may be responsive to diverse environmental stresses and regulate distinct biological processes. Besides, the expression patterns of the GmGRAS genes were varied in various tissues, during saline and dehydration stresses and during seed germination processes.

CONCLUSIONS

We conducted a systematic investigation of the GRAS genes in soybean, which may be valuable in paving the way for future GmGRAS gene studies and soybean breeding.

Overexpression of GmUBC9 Gene Enhances Plant Drought Resistance and Affects Flowering Time via Histone H2B Monoubiquitination

Ubiquitylation is a form of post-translational modification of proteins that can alter localization, functionality, degradation, or transcriptional activity within a cell. E2 ubiquitin-conjugating enzyme (UBC) and E3 ubiquitin ligases are the primary determinants of substrate specificity in the context of ubiquitin conjugation. Multiubiquitination modifies target proteins for 26S proteasome degradation, while monoubiquitination controls protein activation and localization. At present, research on the monoubiquitination, especially histone monoubiquitination, has mostly focused on model plants with relatively few on crop species. In this study, we identified 91 UBC-like genes in soybean. The chromosomal localization, phylogenetic relationships, gene structures, and putative cis-acting elements were evaluated. Furthermore, the tissue-specific expression patterns of UBC Class I genes under drought stress were also investigated. Among Class I genes, GmUBC9 induction in response to drought stress was evident, and so this gene was selected for further analysis. GmUBC9 localized to the nucleus and endoplasmic reticulum. The overexpression of GmUBC9 in Arabidopsis led to enhanced tolerance for drought conditions across a range of stages of development, while overexpression in soybean hairy roots similarly led to improvements in tolerance for drought conditions, increased proline content, and reduced MDA content in soybean seedlings compared to wild type plants. HISTONE MONOUBIQUITINATION 2 (HUB2), an E3-like protein involved in histone H2B ubiquitylation (H2Bub1), was found to interact with GmUBC9 through Y2H analysis and BiFC assays in Arabidopsis and soybean. Under drought conditions, the level of H2Bub1 increased, and transcription of drought response genes was activated in GmUBC9 transgenic Arabidopsis and soybean. In addition, GmUBC9 transgenic Arabidopsis and soybean showed a late-flowering phenotype and had increased expression levels of the flowering related genes FLC and MAF4. These findings indicate that GmUBC9 is important for drought stress response and regulation of flowering time in soybean.

Systematic analysis of 1298 RNA-Seq samples and construction of a comprehensive soybean (Glycine max) expression atlas

Soybean (Glycine max [L.] Merr.) is a major crop in animal feed and human nutrition, mainly for its rich protein and oil contents. The remarkable rise in soybean transcriptome studies over the past 5 years generated an enormous amount of RNA-seq data, encompassing various tissues, developmental conditions and genotypes. In this study, we have collected data from 1298 publicly available soybean transcriptome samples, processed the raw sequencing reads and mapped them to the soybean reference genome in a systematic fashion. We found that 94% of the annotated genes (52 737/56 044) had detectable expression in at least one sample. Unsupervised clustering revealed three major groups, comprising samples from aerial, underground and seed/seed-related parts. We found 452 genes with uniform and constant expression levels, supporting their roles as housekeeping genes. On the other hand, 1349 genes showed heavily biased expression patterns towards particular tissues. A transcript-level analysis revealed that 95% (70 963 of 74 490) of the assembled transcripts have intron chains exactly matching those from known transcripts, whereas 3256 assembled transcripts represent potentially novel splicing isoforms. The dataset compiled here constitute a new resource for the community, which can be downloaded or accessed through a user-friendly web interface at http://venanciogroup.uenf.br/resources/. This comprehensive transcriptome atlas will likely accelerate research on soybean genetics and genomics.

Successful field performance in warm and dry environments of soybean expressing the sunflower transcription factor HB4

Soybean yield is limited primarily by abiotic constraints. No transgenic soybean with improved abiotic stress tolerance is commercially available. We transformed soybean plants with genetic constructs able to express the sunflower transcription factor HaHB4, which confers drought tolerance to Arabidopsis and wheat. One line (b10H) carrying the sunflower promoter was chosen among three independent lines because it exhibited the best performance in seed yield, and was evaluated in the greenhouse and in 27 field trials in different environments in Argentina. In greenhouse experiments, transgenic plants showed increased seed yield under stress conditions together with greater epicotyl diameter, larger xylem area, and increased water use efficiency compared with controls. They also exhibited enhanced seed yield in warm and dry field conditions. This response was accompanied by an increase in seed number that was not compensated by a decrease in individual seed weight. Transcriptome analysis of plants from a field trial with maximum difference in seed yield between genotypes indicated the induction of genes encoding redox and heat shock proteins in b10H. Collectively, our results indicate that soybeans transformed with HaHB4 are expected to have a reduced seed yield penalty when cultivated in warm and dry conditions, which constitute the best target environments for this technology.

Transcriptome profile analysis of two Vicia faba cultivars with contrasting salinity tolerance during seed germination

Faba bean (Vicia faba L.) is an important food legume crop. Salinity soils severely constrain the production of faba bean, however, the seed germination of faba bean, which is a vital plant growth stage, is sensitive to salinity. Planting improved varieties of faba bean, which exhibit salt tolerance in seed germination stage, is an optimal strategy for faba bean product. To investigate the genes dynamics during the seed germination stage under salinity, RNA-seq method was used to investigate genome-wide transcription profiles of two faba bean varieties with contrast salt-tolerance during the seed germination. A total of 4,486 differentially expressed genes (DEGs) were identified among the comparison of salt-tolerant variety Y134 and salt-sensitive variety Y078 treated with salinity or not. Of these, 1,410 candidate DEGs were identified as salt-stress response genes. Furthermore, 623 DEGs were identified as variety-specific response gene during seed germination at 16 h or 24 h with salt treatment. Based on the pathway enrichment according to the Kyoto Encyclopedia of Genes and Genomes database (KEGG), these DEGs involving in cell wall loosening (e.g., xyloglucan endotransglucosylase/hydrolase, chitinase, and expansin), hormone metabolism (e.g., LEA genes, genes associated with ABA or ethylene signal pathway), chromatin remodeling (e.g., chromatin structure proteins, LHP1), small interfering RNA pathway, etc., were significantly up-regulated in salt-tolerance variety with salt treatment, indicating that they play critical roles in regulation of seed germination. The results indicated that a clearer mechanism of gene regulation that regulates the seed germination responding to salinity in faba bean. These findings are helpful to increase the understanding of the salt tolerance mechanism of crops during seed germination, and provide valuable genetic resource for the breeding of salt-tolerant faba bean varieties in future.

Molecular adaptation to salinity fluctuation in tropical intertidal environments of a mangrove tree Sonneratia alba

BACKGROUND

Mangroves have adapted to intertidal zones - the interface between terrestrial and marine ecosystems. Various studies have shown adaptive evolution in mangroves at physiological, ecological, and genomic levels. However, these studies paid little attention to gene regulation of salt adaptation by transcriptome profiles.

RESULTS

We sequenced the transcriptomes of Sonneratia alba under low (fresh water), medium (half the seawater salinity), and high salt (seawater salinity) conditions and investigated the underlying transcriptional regulation of salt adaptation. In leaf tissue, 64% potential salinity-related genes were not differentially expressed when salinity increased from freshwater to medium levels, but became up- or down-regulated when salt concentrations further increased to levels found in sea water, indicating that these genes are well adapted to the medium saline condition. We inferred that both maintenance and regulation of cellular environmental homeostasis are important adaptive processes in S. alba. i) The sulfur metabolism as well as flavone and flavonol biosynthesis KEGG pathways were significantly enriched among up-regulated genes in leaves. They are both involved in scavenging ROS or synthesis and accumulation of osmosis-related metabolites in plants. ii) There was a significantly increased percentage of transcription factor-encoding genes among up-regulated transcripts. High expressions of salt tolerance-related TF families were found under high salt conditions. iii) Some genes up-regulated in response to salt treatment showed signs of adaptive evolution at the amino acid level and might contribute to adaptation to fluctuating intertidal environments.

CONCLUSIONS

This study first elucidates the mechanism of high-salt adaptation in mangroves at the whole-transcriptome level by salt gradient experimental treatments. It reveals that several candidate genes (including salt-related genes, TF-encoding genes, and PSGs) and major pathways are involved in adaptation to high-salt environments. Our study also provides a valuable resource for future investigation of adaptive evolution in extreme environments.

Genome-wide characterization and expression analysis of soybean trihelix gene family

Trihelix transcription factors play multiple roles in plant growth, development and various stress responses. In this study, we identified 71 trihelix family genes in the soybean genome. These trihelix genes were located at 19 out of 20 soybean chromosomes unevenly and were classified into six distinct subfamilies: GT-1, GT-2, GTγ, SIP1, SH4 and GTδ. The gene structure and conserved functional domain of these trihelix genes were similar in the same subfamily but diverged between different subfamilies. Thirteen segmental duplicated gene pairs were identified and all of them experienced a strong purifying selective pressure during evolution. Various stress-responsive cis-elements presented in the promoters of soybean trihelix genes, suggesting that the trihelix genes might respond to the environmental stresses in soybean. The expression analysis suggests that trihelix genes are involved in diverse functions during soybean development, flood or salinity tolerance, and plant immunity. Our results provide genomic information of the soybean trihelix genes and a basis for further characterizing their roles in response to environmental stresses.

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.

Genome-wide characterization of NtHD-ZIP IV: different roles in abiotic stress response and glandular Trichome induction

BACKGROUND

The plant-specific homeodomain-leucine zipper class IV (HD-ZIP IV) gene family has been involved in the regulation of epidermal development.

RESULTS

Fifteen genes coding for HD-ZIP IV proteins were identified (NtHD-ZIP-IV-1 to NtHD-ZIP-IV-15) based on the genome of N. tabacum. Four major domains (HD, ZIP, SAD and START) were present in these proteins. Tissue expression pattern analysis indicated that NtHD-ZIP-IV-1, - 2, - 3, - 10, and - 12 may be associated with trichome development; NtHD-ZIP-IV-8 was expressed only in cotyledons; NtHD-ZIP-IV-9 only in the leaf and stem epidermis; NtHD-ZIP-IV-11 only in leaves; and NtHD-ZIP-IV-15 only in the root and stem epidermis. We found that jasmonates may induce the generation of glandular trichomes, and that NtHD-ZIP-IV-1, - 2, - 5, and - 7 were response to MeJA treatment. Dynamic expression under abiotic stress and after application of phytohormones indicated that most NtHD-ZIP IV genes were induced by heat, cold, salt and drought. Furthermore, most of these genes were induced by gibberellic acid, 6-benzylaminopurine, and salicylic acid, but were inhibited by abscisic acid. NtHD-ZIP IV genes were sensitive to heat, but insensitive to osmotic stress.

CONCLUSION

NtHD-ZIP IV genes are implicated in a complex regulatory gene network controlling epidermal development and abiotic stress responses. The present study provides evidence to elucidate the gene functions of NtHD-ZIP IVs during epidermal development and stress response.

Genome-wide characterization and expression analysis of the HD-Zip gene family in response to drought and salinity stresses in sesame

Background

The homeodomain-leucine zipper (HD-Zip) gene family is one of the plant-specific transcription factor families, involved in plant development, growth, and in the response to diverse stresses. However, comprehensive analysis of the HD-Zip genes, especially those involved in response to drought and salinity stresses is lacking in sesame (Sesamum indicum L.), an important oil crop in tropical and subtropical areas.

Results

In this study, 45 HD-Zip genes were identified in sesame, and denominated as SiHDZ01-SiHDZ45. Members of SiHDZ family were classified into four groups (HD-Zip I-IV) based on the phylogenetic relationship of Arabidopsis HD-Zip proteins, which was further supported by the analysis of their conserved motifs and gene structures. Expression analyses of SiHDZ genes based on transcriptome data showed that the expression patterns of these genes were varied in different tissues. Additionally, we showed that at least 75% of the SiHDZ genes were differentially expressed in responses to drought and salinity treatments, and highlighted the important role of HD-Zip I and II genes in stress responses in sesame.

Conclusions

This study provides important information for functional characterization of stress-responsive HD-Zip genes and may contribute to the better understanding of the molecular basis of stress tolerance in sesame.

Differential expression of recently duplicated PTOX genes in Glycine max during plant development and stress conditions

Plastid terminal oxidase (PTOX) is a chloroplast enzyme that catalyzes oxidation of plastoquinol (PQH2) and reduction of molecular oxygen to water. Its function has been associated with carotenoid biosynthesis, chlororespiration and environmental stress responses in plants. In the majority of plant species, a single gene encodes the protein and little is known about events of PTOX gene duplication and their implication to plant metabolism. Previously, two putative PTOX (PTOX1 and 2) genes were identified in Glycine max, but the evolutionary origin and the specific function of each gene was not explored. Phylogenetic analyses revealed that this gene duplication occurred apparently during speciation involving the Glycine genus ancestor, an event absent in all other available plant leguminous genomes. Gene expression evaluated by RT-qPCR and RNA-seq data revealed that both PTOX genes are ubiquitously expressed in G. max tissues, but their mRNA levels varied during development and stress conditions. In development, PTOX1 was predominant in young tissues, while PTOX2 was more expressed in aged tissues. Under stress conditions, the PTOX transcripts varied according to stress severity, i.e., PTOX1 mRNA was prevalent under mild or moderate stresses while PTOX2 was predominant in drastic stresses. Despite the high identity between proteins (97%), molecular docking revealed that PTOX1 has higher affinity to substrate plastoquinol than PTOX2. Overall, our results indicate a functional relevance of this gene duplication in G. max metabolism, whereas PTOX1 could be associated with chloroplast effectiveness and PTOX2 to senescence and/or apoptosis.

Transcriptome profiling of Gerbera hybrida reveals that stem bending is caused by water stress and regulation of abscisic acid

Background

Gerbera hybrida is one of the most popular cut flowers in the world; however, stem bending, which always happens when gerbera flower harvested from the field, greatly limits its vase life. To date the molecular mechanisms underlying stem bending remain poorly understood.

Results

In this study, we performed high-throughput transcriptome sequencing of gerbera during stem bending using the Illumina sequencing technology. Three cDNA libraries constructed from mRNAs of gerbera stem at stem bending stage 0, 2 and 4 were sequenced. More than 300 million high-quality reads were generated and assembled into 96,492 unigenes. Among them, 34,166 unigenes were functionally annotated based on similarity search with known protein. Sequences derived from plants at different stem bending stages were mapped to the assembled transcriptome, and 9,406 differentially expressed genes (DEGs) were identified. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, specific pathways were identified during the stem bending process, such as phenylpropanoid biosynthesis pathway, phenylalanine metabolism pathway, starch and sucrose metabolism pathway, and plant hormone signal transduction pathway. A total of 211 transcription factors (TFs), including TF families involved in plant senescence, such as NAC, MYB, WRKY, and AP2/ERF members, as well as TFs related to water stress tolerance, were shown to be regulated during stem bending. Gene Onotology (GO) functional enrichment analysis indicated that key genes involved in responses to osmotic and oxidative stresses were also varied in expression during this process. Furthermore, analysis of DEGs involved in the hormone signaling pathways and determination of endogenous abscisic acid (ABA) content showed that stem bending may be an ethylene-independent process, but regulated by ABA. In short, our findings suggested that the stem bending of cut gerbera may be caused by the involvement of water stress and regulation of ABA during the postharvest life.

Conclusions

The transcriptome sequences provide a valuable resource in revealing the molecular mechanism underlying stem bending of cut flower and offer novel genes that can be used to guide future studies for ornamental plant breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5961-1) contains supplementary material, which is available to authorized users.

RNA sequencing analysis of salt tolerance in soybean (Glycine max)

Salt stress causes foliar chlorosis and scorch, plant stunting, and eventually yield reduction in soybean. There are differential responses, namely tolerance (excluder) and intolerance (includer), among soybean germplasm. However, the genetic and physiological mechanisms for salt tolerance is complex and not clear yet. Based on the results from the screening of the RA-452 x Osage mapping population, two F4:6 lines with extreme responses, most tolerant and most sensitive, were selected for a time-course gene expression study in which the 250 mM NaCl treatment was initially imposed at the V1 stage and continued for 24 h (hrs). Total RNA was isolated from the leaves harvested at 0, 6, 12, 24 h after the initiation of salt treatment, respectively. The RNA-Seq analysis was conducted to compare the salt tolerant genotype with salt sensitive genotype at each time point using RNA-Seq pipeline method. A total of 2374, 998, 1746, and 630 differentially expressed genes (DEGs) between salt-tolerant line and salt-sensitive line, were found at 0, 6, 12, and 24 h, respectively. The expression patterns of 154 common DEGs among all the time points were investigated, of which, six common DEGs were upregulated and seven common DEGs were downregulated in salt-tolerant line. Moreover, 13 common DEGs were dramatically expressed at all the time points. Based on Log2 (fold change) of expression level of salt-tolerant line to salt-sensitive line and gene annotation, Glyma.02G228100, Glyma.03G226000, Glyma.03G031000, Glyma.03G031400, Glyma.04G180300, Glyma.04G180400, Glyma.05 g204600, Glyma.08G189600, Glyma.13G042200, and Glyma.17G173200, were considered to be the key potential genes involving in the salt-tolerance mechanism in the soybean salt-tolerant line.