Soybean steroids improve crop abiotic stress tolerance and increase yield.

Sterols have long been associated with diverse fields, such as cancer treatment, drug development, and plant growth; however, their underlying mechanisms and functions remain enigmatic. Here, we unveil a critical role played by a GmNF-YC9-mediated CCAAT-box transcription complex in modulating the steroid metabolism pathway within soybeans. Specifically, this complex directly activates squalene monooxygenase (GmSQE1), which is a rate-limiting enzyme in steroid synthesis. Our findings demonstrate that overexpression of either GmNF-YC9 or GmSQE1 significantly enhances soybean stress tolerance, while the inhibition of SQE weakens this tolerance. Field experiments conducted over two seasons further reveal increased yields per plant in both GmNF-YC9 and GmSQE1 overexpressing plants under drought stress conditions. This enhanced stress tolerance is attributed to the reduction of abiotic stress-induced cell oxidative damage. Transcriptome and metabolome analyses shed light on the upregulation of multiple sterol compounds, including fucosterol and soyasaponin II, in GmNF-YC9 and GmSQE1 overexpressing soybean plants under stress conditions. Intriguingly, the application of soybean steroids, including fucosterol and soyasaponin II, significantly improves drought tolerance in soybean, wheat, foxtail millet, and maize. These findings underscore the pivotal role of soybean steroids in countering oxidative stress in plants and offer a new research strategy for enhancing crop stress tolerance and quality from gene regulation to chemical intervention.

GmDof41 regulated by the DREB1-type protein improves drought and salt tolerance by regulating the DREB2-type protein in soybean.

Despite their essential and multiple roles in biological processes, the molecular mechanism of Dof transcription factors (TFs) for responding to abiotic stresses is rarely reported in plants. We identified a soybean Dof gene GmDof41 which was involved in the responses to drought, salt, and exogenous ABA stresses. Overexpression of GmDof41 in soybean transgenic hairy roots attenuated H2O2 accumulation and regulated proline homeostasis, resulting in the drought and salt tolerance. Yeast one-hybrid and electrophoretic mobility shift assay (EMSA) illustrated that GmDof41 was regulated by the DREB1-type protein GmDREB1B;1 that could improve drought and salt tolerance in plants. Further studies illustrated GmDof41 can directly bind to the promoter of GmDREB2A which encodes a DREB2-type protein and affects abiotic stress tolerance in plants. Collectively, our results suggested that GmDof41 positively regulated drought and salt tolerance by correlating with GmDREB1B;1 and GmDREB2A. This study provides an important basis for further exploring the abiotic stress-tolerance mechanism of Dof TFs in soybean.

Genome-Wide Analysis of the Soybean TIFY Family and Identification of GmTIFY10e and GmTIFY10g Response to Salt Stress.

TIFY proteins play crucial roles in plant abiotic and biotic stress responses. Our transcriptome data revealed several TIFY family genes with significantly upregulated expression under drought, salt, and ABA treatments. However, the functions of the GmTIFY family genes are still unknown in abiotic stresses. We identified 38 GmTIFY genes and found that TIFY10 homologous genes have the most duplication events, higher selection pressure, and more obvious response to abiotic stresses compared with other homologous genes. Expression pattern analysis showed that GmTIFY10e and GmTIFY10g genes were significantly induced by salt stress. Under salt stress, GmTIFY10e and GmTIFY10g transgenic Arabidopsis plants showed higher root lengths and fresh weights and had significantly better growth than the wild type (WT). In addition, overexpression of GmTIFY10e and GmTIFY10g genes in soybean improved salt tolerance by increasing the PRO, POD, and CAT contents and decreasing the MDA content; on the contrary, RNA interference plants showed sensitivity to salt stress. Overexpression of GmTIFY10e and GmTIFY10g in Arabidopsis and soybean could improve the salt tolerance of plants, while the RNAi of GmTIFY10e and GmTIFY10g significantly increased sensitivity to salt stress in soybean. Further analysis demonstrated that GmTIFY10e and GmTIFY10g genes changed the expression levels of genes related to the ABA signal pathway, including GmSnRK2, GmPP2C, GmMYC2, GmCAT1, and GmPOD. This study provides a basis for comprehensive analysis of the role of soybean TIFY genes in stress response in the future.

Proteomics analysis reveals that foreign cp4-epsps gene regulates the levels of shikimate and branched pathways in genetically modified soybean line H06-698.

Although genetically modified (GM) glyphosate-resistant soybeans with cp4-epsps gene have been widely planted all over the world, their proteomic characteristics are not very clear. In this study, the soybean seeds of a GM soybean line H06-698 (H) with cp4-epsps gene and its non-transgenic counterpart Mengdou12 (M), which were collected from two experiment fields in two years and used as 4 sample groups, were analyzed with label-free proteomics technique. A total of 1706 proteins were identified quantitatively by label-free quantification, and a total of 293 proteins were detected as common differential abundance proteins (DAPs, FC is not less than 1.5) both in two groups or more. Functional enrichment analysis of common DAPs identified from four groups, shows that most up-regulated proteins were clustered into stress response, carbon and energy metabolism, and genetic information processing. Further documentary analysis shows that 15 proteins play important roles in shikimate pathways, reactive oxygen species (ROS) and stress response. These results indicated that the change of protein abundance in different samples were affected by various factors, but except shikimate and branched pathways related proteins, only ROS and stress-related proteins were found to be stably regulated by cp4-epsps gene, and no unexpected and safety-related proteins such as antinutritional factors, allergenic proteins, and toxic proteins were found as DAPs. The influence of foreign genes in genetically modified plants is worthy of attention and this work provides new clues for exploring the regulated proteins and pathways in GM plants.

Genome-Wide Characterization and Expression Analysis of Soybean TGA Transcription Factors Identified a Novel TGA Gene Involved in Drought and Salt Tolerance.

The TGA transcription factors, a subfamily of bZIP group D, play crucial roles in various biological processes, including the regulation of growth and development as well as responses to pathogens and abiotic stress. In this study, 27 TGA genes were identified in the soybean genome. The expression patterns of GmTGA genes showed that several GmTGA genes are differentially expressed under drought and salt stress conditions. Among them, GmTGA17 was strongly induced by both stress, which were verificated by the promoter-GUS fusion assay. GmTGA17 encodes a nuclear-localized protein with transcriptional activation activity. Heterologous and homologous overexpression of GmTGA17 enhanced tolerance to drought and salt stress in both transgeinc Arabidopsis plants and soybean hairy roots. However, RNAi hairy roots silenced for GmTGA17 exhibited an increased sensitivity to drought and salt stress. In response to drought or salt stress, transgenic Arabidopsis plants had an increased chlorophyll and proline contents, a higher ABA content, a decreased MDA content, a reduced water loss rate, and an altered expression of ABA- responsive marker genes compared with WT plants. In addition, transgenic Arabidopsis plants were more sensitive to ABA in stomatal closure. Similarly, measurement of physiological parameters showed an increase in chlorophyll and proline contents, with a decrease in MDA content in soybean seedlings with overexpression hairy roots after drought and salt stress treatments. The opposite results for each measurement were observed in RNAi lines. This study provides new insights for functional analysis of soybean TGA transcription factors in abiotic stress.

The WRKY Transcription Factor GmWRKY12 Confers Drought and Salt Tolerance in Soybean.

WRKYs are important regulators in plant development and stress responses. However, knowledge of this superfamily in soybean is limited. In this study, we characterized the drought- and salt-induced gene GmWRKY12 based on RNA-Seq and qRT-PCR. GmWRKY12, which is 714 bp in length, encoded 237 amino acids and grouped into WRKY II. The promoter region of GmWRKY12 included ABER4, MYB, MYC, GT-1, W-box and DPBF cis-elements, which possibly participate in abscisic acid (ABA), drought and salt stress responses. GmWRKY12 was minimally expressed in different tissues under normal conditions but highly expressed under drought and salt treatments. As a nucleus protein, GmWRKY12 was responsive to drought, salt, ABA and salicylic acid (SA) stresses. Using a transgenic hairy root assay, we further characterized the roles of GmWRKY12 in abiotic stress tolerance. Compared with control (Williams 82), overexpression of GmWRKY12 enhanced drought and salt tolerance, increased proline (Pro) content and decreased malondialdehyde (MDA) content under drought and salt treatment in transgenic soybean seedlings. These results may provide a basis to understand the functions of GmWRKY12 in abiotic stress responses in soybean.

Development and utilization of a new chemically-induced soybean library with a high mutation density .

Mutagenized populations have provided important materials for introducing variation and identifying gene function in plants. In this study, an ethyl methanesulfonate (EMS)-induced soybean (Glycine max) population, consisting of 21,600 independent M2 lines, was developed. Over 1,000 M4 (5) families, with diverse abnormal phenotypes for seed composition, seed shape, plant morphology and maturity that are stably expressed across different environments and generations were identified. Phenotypic analysis of the population led to the identification of a yellow pigmentation mutant, gyl, that displayed significantly decreased chlorophyll (Chl) content and abnormal chloroplast development. Sequence analysis showed that gyl is allelic to MinnGold, where a different single nucleotide polymorphism variation in the Mg-chelatase subunit gene (ChlI1a) results in golden yellow leaves. A cleaved amplified polymorphic sequence marker was developed and may be applied to marker-assisted selection for the golden yellow phenotype in soybean breeding. We show that the newly developed soybean EMS mutant population has potential for functional genomics research and genetic improvement in soybean.

Co-expression of G2-EPSPS and glyphosate acetyltransferase GAT genes conferring high tolerance to glyphosate in soybean.

Glyphosate is a widely used non-selective herbicide with broad spectrum of weed control around the world. At present, most of the commercial glyphosate tolerant soybeans utilize glyphosate tolerant gene CP4-EPSPS or glyphosate acetyltransferase gene GAT separately. In this study, both glyphosate tolerant gene G2-EPSPS and glyphosate degraded gene GAT were co-transferred into soybean and transgenic plants showed high tolerance to glyphosate. Molecular analysis including PCR, Sothern blot, qRT-PCR, and Western blot revealed that target genes have been integrated into genome and expressed effectively at both mRNA and protein levels. Furthermore, the glyphosate tolerance analysis showed that no typical symptom was observed when compared with a glyphosate tolerant line HJ06-698 derived from GR1 transgenic soybean even at fourfold labeled rate of Roundup. Chlorophyll and shikimic acid content analysis of transgenic plant also revealed that these two indexes were not significantly altered after glyphosate application. These results indicated that co-expression of G2-EPSPS and GAT conferred high tolerance to the herbicide glyphosate in soybean. Therefore, combination of tolerant and degraded genes provides a new strategy for developing glyphosate tolerant transgenic crops.

Functional marker development of miR1511-InDel and allelic diversity within the genus Glycine.

BACKGROUND: Single-stranded non-protein coding small RNAs, 18-25 nucleotides in length, are ubiquitous throughout plants genomes and are involved in post-transcriptional gene regulation. Several types of DNA markers have been reported for the detection of genetic diversity or sequence variation in soybean, one of the most important legume crops in worldwide for seed protein and oil content. Recently, with the available of public genomic databases, there has been a shift from the labor-intensive development of PCR-based markers to sequence-based genotyping and the development of functional markers within genes, often coupled with the use of RNA information. But thus far miRNA-based markers have been only developed in rice and tobacco. Here we report the first functional molecular miRNA marker, miR1511-InDel, in soybean for a specific single copy locus used to assess genetic variation in domesticated soybean (Glycine max [L.] Merr) and its wild progenitor (Glycine soja Sieb. & Zucc.). RESULTS: We genotyped a total of 1,669 accessions of domesticated soybean (G. max) and its wild progenitor G. soja which are native throughout the China and parts of Korea, Japan and Russia. The results indicate that the miR1511 locus is distributed in cultivated soybean and has three alleles in annual wild soybean. Based on this result, we proposed that miR-InDel marker technology can be used to assess genetic variation. The inclusion of geo-reference data with miR1511-InDel marker data corroborated that accessions from the Yellow River basin (Huanghuai) exhibited high genetic diversity which provides more molecular evidence for gene diversity in annual wild soybean and domestication of soybean. CONCLUSIONS: These results provide evidence for the use of RNA marker, miRNA1511-InDel, as a soybean-specific functional maker for the study of genetic diversity, genotyping of germplasm and evolution studies. This is also the first report of functional marker developed from soybean miRNA located within the functional region of pre-miRNA1511.

Salinity tolerance in soybean is modulated by natural variation in GmSALT3.

The identification of genes that improve the salt tolerance of crops is essential for the effective utilization of saline soils for agriculture. Here, we use fine mapping in a soybean (Glycine max (L.) Merr.) population derived from the commercial cultivars Tiefeng 8 and 85-140 to identify GmSALT3 (salt tolerance-associated gene on chromosome 3), a dominant gene associated with limiting the accumulation of sodium ions (Na+) in shoots and a substantial enhancement in salt tolerance in soybean. GmSALT3 encodes a protein from the cation/H+ exchanger family that we localized to the endoplasmic reticulum and which is preferentially expressed in the salt-tolerant parent Tiefeng 8 within root cells associated with phloem and xylem. We identified in the salt-sensitive parent, 85-140, a 3.78-kb copia retrotransposon insertion in exon 3 of Gmsalt3 that truncates the transcript. By sequencing 31 soybean landraces and 22 wild soybean (Glycine soja) a total of nine haplotypes including two salt-tolerant haplotypes and seven salt-sensitive haplotypes were identified. By analysing the distribution of haplotypes among 172 Chinese soybean landraces and 57 wild soybean we found that haplotype 1 (H1, found in Tiefeng 8) was strongly associated with salt tolerance and is likely to be the ancestral allele. Alleles H2-H6, H8 and H9, which do not confer salinity tolerance, were acquired more recently. H1, unlike other alleles, has a wide geographical range including saline areas, which indicates it is maintained when required but its potent stress tolerance can be lost during natural selection and domestication. GmSALT3 is a gene associated with salt tolerance with great potential for soybean improvement.

De novo assembly of soybean wild relatives for pan-genome analysis of diversity and agronomic traits.

Wild relatives of crops are an important source of genetic diversity for agriculture, but their gene repertoire remains largely unexplored. We report the establishment and analysis of a pan-genome of Glycine soja, the wild relative of cultivated soybean Glycine max, by sequencing and de novo assembly of seven phylogenetically and geographically representative accessions. Intergenomic comparisons identified lineage-specific genes and genes with copy number variation or large-effect mutations, some of which show evidence of positive selection and may contribute to variation of agronomic traits such as biotic resistance, seed composition, flowering and maturity time, organ size and final biomass. Approximately 80% of the pan-genome was present in all seven accessions (core), whereas the rest was dispensable and exhibited greater variation than the core genome, perhaps reflecting a role in adaptation to diverse environments. This work will facilitate the harnessing of untapped genetic diversity from wild soybean for enhancement of elite cultivars.

Role of physiological mechanisms and EPSPS gene expression in glyphosate resistance in wild soybeans (Glycine soja).

The physiological mechanisms underlying glyphosate resistance in wild soybean germplasm and relevant EPSPS gene expression were evaluated. These germplasms were selected by gradually increasing glyphosate selection pressure started from 2010. As indicated by a whole-plant dose response bioassay, ZYD-254 plants were resistant to glyphosate at concentrations of 1230gaeha(-1), but the susceptible plants (ZYD-16) were unable to survive in the presence of 300gaeha(-1) glyphosate. The ED50 values of resistant germplasm were approximately 8.8 times of the susceptible germplasm. Chlorophyll content was significantly decreased in ZYD-16 plants in comparison with ZYD-254 plants. ZYD-16 plants accumulated 10.1 times more shikimate in leaves at 5days after glyphosate treatment at 1230gaeha(-1) than ZYD-254 did. GST activity differed between ZYD-254 and ZYD-16 in three tissues. It was highest in leaves. There were no significant differences in EPSPS1 or EPSPS3 expression between two germplasms before exposure to glyphosate treatment. After glyphosate treatment, there was a 2- to 4-fold increase in EPSPS1 mRNA levels in ZYD-254, but there was no change in EPSPS3 mRNA levels in ZYD-254 or ZYD-16.

Molecular footprints of domestication and improvement in soybean revealed by whole genome re-sequencing.

BACKGROUND: Artificial selection played an important role in the origin of modern Glycine max cultivars from the wild soybean Glycine soja. To elucidate the consequences of artificial selection accompanying the domestication and modern improvement of soybean, 25 new and 30 published whole-genome re-sequencing accessions, which represent wild, domesticated landrace, and Chinese elite soybean populations were analyzed. RESULTS: A total of 5,102,244 single nucleotide polymorphisms (SNPs) and 707,969 insertion/deletions were identified. Among the SNPs detected, 25.5% were not described previously. We found that artificial selection during domestication led to more pronounced reduction in the genetic diversity of soybean than the switch from landraces to elite cultivars. Only a small proportion (2.99%) of the whole genomic regions appear to be affected by artificial selection for preferred agricultural traits. The selection regions were not distributed randomly or uniformly throughout the genome. Instead, clusters of selection hotspots in certain genomic regions were observed. Moreover, a set of candidate genes (4.38% of the total annotated genes) significantly affected by selection underlying soybean domestication and genetic improvement were identified. CONCLUSIONS: Given the uniqueness of the soybean germplasm sequenced, this study drew a clear picture of human-mediated evolution of the soybean genomes. The genomic resources and information provided by this study would also facilitate the discovery of genes/loci underlying agronomically important traits.

Glyphosate effects on the gene expression of the apical bud in soybean (Glycine max).

Glyphosate is a broad spectrum, non-selective herbicide which has been widely used for weed control. Much work has focused on elucidating the high accumulation of glyphosate in shoot apical bud (shoot apex). However, to date little is known about the molecular mechanisms of the sensitivity of shoot apical bud to glyphosate. Global gene expression profiling of the soybean apical bud response to glyphosate treatment was performed in this study. The results revealed that the glyphosate inhibited tryptophan biosynthesis of the shikimic acid pathway in the soybean apical bud, which was the target site of glyphosate. Glyphosate inhibited the expression of most of the target herbicide site genes. The promoter sequence analysis of key target genes revealed that light responsive elements were important regulators in glyphosate induction. These results will facilitate further studies of cloning genes and molecular mechanisms of glyphosate on soybean shoot apical bud.

A comparative miRNAome analysis reveals seven fiber initiation-related and 36 novel miRNAs in developing cotton ovules.

An increasing number of microRNAs (miRNAs) have been shown to play crucial regulatory roles in the process of plant development. Here, we used high-throughput sequencing combined with computational analysis to characterize miRNAomes from the ovules of wild-type upland cotton and a fiberless mutant during fiber initiation. Comparative miRNAome analysis combined with northern blotting and RACE-PCR revealed seven fiber initiation-related miRNAs expressed in cotton ovules and experimentally validated targets of these miRNAs are involved in different cellular responses and metabolic processes, including transcriptional regulation, auxin and gibberellin signal transduction, actin bundles, and lignin biosynthesis. This paper describes a complex regulatory network consisting of these miRNAs expressed in cotton ovules to coordinate fiber initiation responses. In addition, 36 novel miRNAs and two conserved miRNAs were newly identified, nearly doubling the number of known cotton miRNA families to a total of 78. Furthermore, a chromatin remodeling complex subunit and a pre-mRNA splicing factor are shown for the first time to be miRNA targets. To our knowledge, this study is the first systematic investigation of fiber initiation-related miRNAs and their targets in the developing cotton ovule, deepening our understanding of the important regulatory functions of miRNAs in cotton fiber initiation.

Molecular characterization of three ethylene responsive element binding factor genes from cotton.

Ethylene-responsive factors (ERFs) are important regulators of plant gene expression. In this study, three novel ERF genes, GhERF2, GhERF3 and GhERF6, were isolated from cotton (Gossypium hirstum) using rapid amplification of cDNA ends-polymerase chain reaction. Transient expression analysis using GhERF-green fluorescent protein fusions showed that these three proteins were targeted to the nucleus. Fusion proteins consisting of GhERF2, GhERF3 or GhERF6 coupled to the GAL4 DNA binding domain strongly activated transcription in yeast. Furthermore, GhERF6 was shown to be able to bind specifically to GCC boxes using a particle bombardment assay in tobacco cells. Semi-quantitative reverse transcription-polymerase chain reaction revealed that GhERF2 and GhERF3 are constitutively expressed in all organs, while GhERF6 is only constitutively expressed in vegetative organs. When plants were treated with ethylene, abscisic acid, salt, cold and drought, the transcripts of GhERF2, GhERF3 and GhERF6 were rapidly induced to high levels. Promoter analysis also indicated that the 5' upstream regions of the three genes possess elements induced by these physiological and environmental factors. Collectively, our data suggest that GhERF2, GhERF3 and GhERF6 might function as positive trans-acting factors in the plant responses to ethylene, abscisic acid and other stresses and provide useful clues for further research into the mechanism of them in regulating cotton multiple stress responses.

Identification and characterization of the novel gene GhDBP2 encoding a DRE-binding protein from cotton (Gossypium hirsutum).

A cDNA encoding one novel DRE-binding protein, GhDBP2, was isolated from cotton seedlings. It is classified into the A-6 group of DREB subfamily based on multiple sequence alignment and phylogenetic characterization. Using semi-quantitative RT-PCR, we found that the GhDBP2 transcripts were greatly induced by drought, NaCl, low temperature and ABA treatments in cotton cotyledons. The DNA-binding properties of GhDBP2 were analyzed by electrophoretic mobility shift assay (EMSA), showing that GhDBP2 successfully binds to the previously characterized DRE cis-element as well as the promoter region of the LEA D113 gene. Consistent with its role as a DNA-binding protein, GhDBP2 is preferentially localized to the nucleus of onion epidermal cells. In addition, when GhDBP2 is transiently expressed in tobacco cells, it activates reporter gene expression driven by the LEA D113 promoter. Taken together, our results indicate that GhDBP2 is a DRE-binding transcriptional activator involved in activation of down-stream genes such as LEA D113 expression through interaction with the DRE element, in response to environmental stresses as well as ABA treatment.

Molecular cloning, expression profile and promoter analysis of a novel ethylene responsive transcription factor gene GhERF4 from cotton (Gossypium hirstum).

Ethylene-responsive element binding factors (ERFs) are plant-specific transcription factors, many of which have been linked to stress responses. A novel ERF gene, designated GhERF4, was isolated by RACE-PCR from Gossypium hirsutum. The GhERF4 cDNA has a total length of 1061bp with an open reading frame of 669bp, encoding a protein of 222 amino acids with a molecular weight of 23.5kDa and a calculated pI of 9.03. Sequence alignment shows that GhERF4 contains a 58 amino acid long AP2/ERF domain and a RKRP nuclear localization signal, and belongs to a group II protein in the ERF subfamily as typified by the C-terminal ERF-associated Amphiphilic Repression (EAR) motif. Southern blot analysis indicates that GhERF4 is a single copy gene in cotton genome. Using green fluorescent protein fusion, we demonstrate that GhERF4 accumulates specifically in the nucleus of onion epidermis cells. Semi-quantitative RT-PCR reveals that GhERF4 is constitutively expressed in true leaves, roots, seeds and stems. The transcripts of GhERF4 accumulate highly and rapidly when plants are treated with exogenous ethylene, salt, cold, drought stresses and exogenous abscisic acid (ABA) treatment, suggesting that GhERF4 is regulated by certain components of the stress signaling pathway. Promoter analysis indicates that the 5' upstream region of GhERF4 possesses some elements induced by physiological and environmental factors. These results indicate that GhERF4 may play an important role in response to ethylene, ABA and environmental stresses.

Molecular cloning and functional characterization of a DREB1/CBF-like gene (GhDREB1L) from cotton.

The transcription factors DREB1s/CBFs play important roles in the regulation of plant resistance to environmental stresses and are quite useful for generating transgenic plants tolerant to these stresses. In the present work, a cDNA encoding DREB1/CBF-like protein (GhDREB1L) from cotton was isolated, and its sequence features, DNA binding preference, and expression patterns of the transcripts were also characterized. GhDREB1L contained one conserved AP2/ERF domain and its amino acid sequence was similar to the DREB1/CBF group of the DREB family from other plants. The DNA-binding domain of GhDREB1L was successfully expressed as a fusion protein in Escherichia coli BL21 (DE3) and purified by Ni-NTA affinity chromatography. Electrophoretic mobility shift assay revealed that the purified GhDREB1L fusion protein had a specific binding activity with the previously characterized DRE element (core sequence, ACCGAC) and also with the DRE-like sequence (core sequence, GCCGAC) in the promoter of the dehydration-responsive late embryogenesis-abundant gene LEA D113. Semi-quantitative RT-PCR showed that GhDREB1L was induced in the cotton cotyledons by low temperature, as well as drought and NaCl treatments. These results suggested that the novel cotton GhDREB1L might play an important role in response to low temperature as well as drought and high salinity through binding to the DRE cis-element.