Soybean Gene Research Articles

Integrated Transcriptome and Metabolome Analysis Reveals Molecular Mechanisms Underlying Resistance to Phytophthora Root Rot.

Soybean production is significantly impacted by Phytophthora root rot (PRR), which is caused by Phytophthora sojae. The nucleotide-binding leucine-rich repeat (NLR) gene family plays a crucial role in plant disease resistance. However, current understanding of the function of soybean NLR genes in resistance to PRR is limited. To address this knowledge gap, transgenic soybean plants overexpressing the NLR gene (Glyma.18g283200) were generated to elucidate the molecular mechanism of resistance. Here, transcript changes and metabolic differences were investigated at three time points (12, 24, and 36 h) after P. sojae infection in hypocotyls of two soybean lines, Dongnong 50 (susceptible line, WT) and Glyma.18g283200 overexpression line (resistant line, OE). Based on the changes in differentially expressed genes (DEGs) in response to P. sojae infection in different lines and at different time points, it was speculated that HOPZ-ACTIVATED RESISTANCE 1 (ZAR1), valine, leucine, and isoleucine degradation, and phytohormone signaling may be involved in the defense response of soybean to P. sojae at the transcriptome level by GO term and KEGG pathway enrichment analysis. Differentially accumulated metabolites (DAMs) analysis revealed that a total of 223 and 210 differential metabolites were identified in the positive ion (POS) and negative ion (NEG) modes, respectively. An integrated pathway-level analysis of transcriptomics (obtained by RNA-seq) and metabolomics data revealed that isoflavone biosynthesis was associated with disease resistance. This work provides valuable insights that can be used in breeding programs aiming to enhance soybean resistance against PRR.

Genome-Wide Analysis of Serine Carboxypeptidase-like Genes in Soybean and Their Roles in Stress Resistance.

The serine carboxypeptidase-like (SCPL) gene family plays a crucial role in the regulation of plant growth, development, and stress response through activities such as acyltransferases in plant secondary metabolism pathways. Although SCPL genes have been identified in various plant species, their specific functions and characteristics in soybean (Glycine max) have not yet been studied. We identified and characterized 73 SCPL genes, grouped into three subgroups based on gene structure and phylogenetic relationships. These genes are distributed unevenly across 20 soybean chromosomes and show varied codon usage patterns influenced by both mutation and selection pressures. Gene ontology (GO) enrichment suggests these genes are involved in plant cell wall regulation and stress responses. Expression analysis in various tissues and under stress conditions, including the presence of numerous stress-related cis-acting elements, indicated that these genes have varied expression patterns. This suggests that they play specialized roles such as modulating plant defense mechanisms against nematode infections, enhancing tolerance to drought and high salinity, and responding to cold stress, thereby helping soybean adapt to environmental stresses. Moreover, the expression of specific GmSCPLs was significantly affected following exposure to nematode infection, drought, high salt (NaCl), and cold stresses. Our findings underscore the potential of SCPL genes in enhancing stress resistance in soybean, providing a valuable resource for future genetic improvement and breeding strategies.

Network pharmacology study on fermented soybeans for the prevention of Alzheimer's disease in older individuals.

Alzheimer's disease (AD) is a progressive neurodegenerative disorder characterized by the disruption of synaptic communication among millions of neurons. Recent research has highlighted the potential therapeutic effectiveness of natural polyphenolic compounds in addressing AD. Soybeans are abundant in polyphenols, and their polyphenolic composition undergoes significant alteration through fermentation by Eurotium cristatum. Through comprehensive database searches, we identified active components within fermented soybean polyphenols and genes associated with AD. Subsequently, we utilized Venn diagrams to analyze the overlap between AD-related genes and these components. Furthermore, we visualized the network between intersecting targets and proteins using Cytoscape software. The anti-AD effects of soybeans were further explored through comprehensive analysis, including protein-protein interaction analysis, pathway enrichment analysis, and molecular docking studies. Our investigation unveiled 6-hydroxydaidzein as a major component of fermented soybean polyphenols, shedding light on its potential therapeutic significance in combating AD. The intersection between target proteins of fermented soybeans and disease-related targets in AD comprised 34 genes. Protein-protein interaction analysis highlighted key potential targets, including glyceraldehyde-3-phosphate dehydrogenase (GAPDH), glycogen synthase kinase 3 beta (GSK3B), amyloid precursor protein (APP), cyclin-dependent kinase 5 (CDK5), and beta-site APP cleaving enzyme 1 (BACE1). Molecular docking results demonstrated a robust binding effect between major components from fermented soybeans and the aforesaid key targets implicated in AD treatment. These findings suggest that fermented soybeans demonstrate a degree of efficacy and present promising prospects in the prevention of AD.

TWAS facilitates gene-scale trait genetic dissection through gene expression, structural variations, and alternative splicing in soybean

A genome-wide association study (GWAS) identifies trait-associated loci, but identifying the causal genes can be a bottleneck, due in part to slow decay of linkage disequilibrium (LD). A transcriptome-wide association study (TWAS) addresses this issue by identifying gene expression–phenotype associations or integrating gene expression quantitative trait loci with GWAS results. Here, we used self-pollinated soybean (Glycine max [L.] Merr.) as a model to evaluate the application of TWAS to the genetic dissection of traits in plant species with slow LD decay. We generated RNA sequencing data for a soybean diversity panel and identified the genetic expression regulation of 29 286 soybean genes. Different TWAS solutions were less affected by LD and were robust to the source of expression, identifing known genes related to traits from different tissues and developmental stages. The novel pod-color gene L2 was identified via TWAS and functionally validated by genome editing. By introducing a new exon proportion feature, we significantly improved the detection of expression variations that resulted from structural variations and alternative splicing. As a result, the genes identified through our TWAS approach exhibited a diverse range of causal variations, including SNPs, insertions or deletions, gene fusion, copy number variations, and alternative splicing. Using this approach, we identified genes associated with flowering time, including both previously known genes and novel genes that had not previously been linked to this trait, providing insights complementary to those from GWAS. In summary, this study supports the application of TWAS for candidate gene identification in species with low rates of LD decay.

Genome-wide identification and expression analysis of the SPL transcription factor family and its response to abiotic stress in Pisum sativum L

Squamous promoter binding protein-like (SPL) genes encode plant-specific transcription factors (TFs) that play essential roles in modulating plant growth, development, and stress response. Pea (Pisum sativum L.) is a coarse grain crop of great importance in food production, biodiversity conservation and molecular genetic research, providing genetic information and nutritional resources for improving agricultural production and promoting human health. However, only limited researches on the structure and functions of SPL genes exist in pea (PsSPLs). In this study, we identified 22 PsSPLs and conducted a genome-wide analysis of their physical characteristics, chromosome distribution, gene structure, phylogenetic evolution and gene expression patterns. As a result, the PsSPLs were unevenly distributed on the seven chromosomes of pea and harbored the SBP domain, which is composed of approximately 76 amino acid residues. The phylogenetic analysis revealed that the PsSPLs clustered into eight subfamilies and showed high homology with SPL genes in soybean. Further analysis showed the presence of segmental duplications in the PsSPLs. The expression patterns of 22 PsSPLs at different tissues, developmental stages and under various stimulus conditions were evaluated by qRT-PCR method. It was found that the expression patterns of PsSPLs from the same subfamily were similar in different tissues, the transcripts of most PsSPLs reached the maximum peak value at 14 days after anthesis in the pod. Abiotic stresses can cause significantly up-regulated PsSPL19 expression with spatiotemporal specificity, in addition, four plant hormones can cause the up-regulated expression of most PsSPLs including PsSPL19 in a time-dependent manner. Therefore, PsSPL19 could be a key candidate gene for signal transduction during pea growth and development, pod formation, abiotic stress and plant hormone response. Our findings should provide insights for the elucidating of development regulation mechanism and breeding for resistance to abiotic stress pea.

Genome-wide analysis of filamentous temperature-sensitive H protease (ftsH) gene family in soybean.

The filamentous temperature-sensitive H protease (ftsH) gene family belongs to the ATP-dependent zinc metalloproteins, and ftsH genes play critical roles in plant chloroplast development and photosynthesis. In this study, we performed genome-wide identification and a systematic analysis of soybean ftsH genes. A total of 18 GmftsH genes were identified. The subcellular localization was predicted to be mainly in cell membranes and chloroplasts, and the gene structures, conserved motifs, evolutionary relationships, and expression patterns were comprehensively analyzed. Phylogenetic analysis of the ftsH gene family from soybean and various other species revealed six distinct clades, all of which showed a close relationship to Arabidopsis thaliana. The members of the GmftsH gene family were distributed on 13 soybean chromosomes, with intron numbers ranging from 3 to 15, 13 pairs of repetitive segment. The covariance between these genes and related genes in different species of Oryza sativa, Zea mays, and Arabidopsis thaliana was further investigated. The transcript expression data revealed that the genes of this family showed different expression patterns in three parts, the root, stem, and leaf, and most of the genes were highly expressed in the leaves of the soybean plants. Fluorescence-based real-time quantitative PCR (qRT-PCR) showed that the expression level of GmftsH genes varied under different stress treatments. Specifically, the genes within this family exhibited various induction levels in response to stress conditions of 4℃, 20% PEG-6000, and 100 mmol/L NaCl. These findings suggest that the GmftsH gene family may play a crucial role in the abiotic stress response in soybeans. It was also found that the GmftsH7 gene was localized on the cell membrane, and its expression was significantly upregulated under 4 ℃ treatment. In summary, by conducting a genome-wide analysis of the GmftsH gene family, a strong theoretical basis is established for future studies on the functionality of GmftsH genes. This research can potentially serve as a guide for enhancing the stress tolerance characteristics of soybean.

Genome-Wide Association Analysis of Yield-Related Traits and Candidate Genes in Vegetable Soybean.

Owing to the rising demand for vegetable soybean products, there is an increasing need for high-yield soybean varieties. However, the complex correlation patterns among quantitative traits with genetic architecture pose a challenge for improving vegetable soybean through breeding. Herein, a genome-wide association study (GWAS) was applied to 6 yield-related traits in 188 vegetable soybean accessions. Using a BLINK model, a total of 116 single nucleotide polymorphisms (SNPs) were identified for plant height, pod length, pod number, pod thickness, pod width, and fresh pod weight. Furthermore, a total of 220 genes were found in the 200 kb upstream and downstream regions of significant SNPs, including 11 genes encoding functional proteins. Among them, four candidate genes, Glyma.13G109100, Glyma.03G183200, Glyma.09G102200, and Glyma.09G102300 were analyzed for significant haplotype variations and to be in LD block, which encode MYB-related transcription factor, auxin-responsive protein, F-box protein, and CYP450, respectively. The relative expression of candidate genes in V030 and V071 vegetable soybean (for the plant height, pod number, and fresh pod weight of V030 were lower than those of the V071 strains) was significantly different, and these genes could be involved in plant growth and development via various pathways. Altogether, we identified four candidate genes for pod yield and plant height from vegetable soybean germplasm. This study provides insights into the genomic basis for improving soybean and crucial genomic resources that can facilitate genome-assisted high-yielding vegetable soybean breeding.

Bowman-Birk Inhibitor Mutants of Soybean Generated by CRISPR-Cas9 Reveal Drastic Reductions in Trypsin and Chymotrypsin Inhibitor Activities.

Despite the high quality of soybean protein, raw soybeans and soybean meal cannot be directly included in animal feed mixtures due to the presence of Kunitz (KTi) and Bowman-Birk protease inhibitors (BBis), which reduces animal productivity. Heat treatment can substantially inactivate trypsin and chymotrypsin inhibitors (BBis), but such treatment is energy-intensive, adds expense, and negatively impacts the quality of seed proteins. As an alternative approach, we have employed CRISPR/Cas9 gene editing to create mutations in BBi genes to drastically lower the protease inhibitor content in soybean seed. Agrobacterium-mediated transformation was used to generate several stable transgenic soybean events. These independent CRISPR/Cas9 events were examined in comparison to wild-type plants using Sanger sequencing, proteomic analysis, trypsin/chymotrypsin inhibitor activity assays, and qRT-PCR. Collectively, our results demonstrate the creation of an allelic series of loss-of-function mutations affecting the major BBi gene in soybean. Mutations in two of the highly expressed seed-specific BBi genes lead to substantial reductions in both trypsin and chymotrypsin inhibitor activities.

GmSABP2-1 encodes methyl salicylate esterase and functions in soybean defense against soybean cyst nematode.

The soybean gene GmSABP2-1 encodes methyl salicylate esterase and its overexpression led to significant reduction in development of pathogenic soybean cyst nematode. Soybean cyst nematode (SCN, Heterodera glycines) is one of the most devastating pests of soybean (Glycine max L. Merr.). In searching for SCN-defense genes, a soybean gene of the methylesterase (MES) family was found to be upregulated in an SCN-resistant soybean line and downregulated in an SCN-susceptible line upon SCN infection. This gene was designated as GmSABP2-1. Here, we report on biochemical and overexpression studies of GmSABP2-1 to examine its possible function in SCN resistance. The protein encoded by GmSABP2-1 is closely related to known methyl salicylate esterases. To determine the biochemical function of GmSABP2-1, a full-length cDNA of GmSABP2-1 was cloned into a protein expression vector and expressed in Escherichia coli. The resulting recombinant GmSABP2-1 was demonstrated to catalyze the demethylation of methyl salicylate. The biochemical properties of GmSABP2-1 were determined. Its apparent Km value was 46.2 ± 2.2μM for methyl salicylate, comparable to those of the known methyl salicylate esterases. To explore the biological significance of GmSABP2-1 in soybean defense against SCN, we first overexpressed GmSABP2-1 in transgenic hairy roots of an SCN-susceptible soybean line. When infected with SCN, GmSABP2-1-overexpressing hairy roots showed 84.5% reduction in the development of SCN beyond J2 stage. To provide further genetic evidence for the role of GmSABP2-1 in SCN resistance, stable transgenic soybean plants overexpressing GmSABP2-1 were produced. Analysis of the GmSABP2-1-overexpressing lines showed a significant reduction in SCN development compared to non-transgenic plants. In conclusion, we demonstrated that GmSABP2-1 encodes methyl salicylate esterase and functions as a resistance-related gene against SCN.

Genome-Wide Identification of Members of the Soybean CBL Gene Family and Characterization of the Functional Role of GmCBL1 in Responses to Saline and Alkaline Stress.

Calcium ions function as key messengers in the context of intracellular signal transduction. The ability of plants to respond to biotic and abiotic stressors is highly dependent on the calcineurin B-like protein (CBL) and CBL-interacting protein kinase (CIPK) signaling network. Here, a comprehensive effort was made to identify all members of the soybean CBL gene family, leading to the identification of 15 total genes distributed randomly across nine chromosomes, including 13 segmental duplicates. All the GmCBL gene subfamilies presented with similar gene structures and conserved motifs. Analyses of the expression of these genes in different tissues revealed that the majority of these GmCBLs were predominantly expressed in the roots. Significant GmCBL expression and activity increases were also observed in response to a range of stress-related treatments, including salt stress, alkaline stress, osmotic stress, or exposure to salicylic acid, brassinosteroids, or abscisic acid. Striking increases in GmCBL1 expression were observed in response to alkaline and salt stress. Subsequent analyses revealed that GmCBL1 was capable of enhancing soybean salt and alkali tolerance through the regulation of redox reactions. These results offer new insight into the complex mechanisms through which the soybean CBL gene family regulates the responses of these plants to environmental stressors, highlighting promising targets for efforts aimed at enhancing soybean stress tolerance.

Silencing GmATG5 genes accelerated senescence and enhanced disease resistance in soybean

Autophagy plays an essential role in recycling/re-utilizing nutrients and in adaptions to numerous stresses. However, the roles of autophagy in soybean have not been investigated extensively. In this study, a virus-induced gene silencing approach mediated by bean pod mottle virus (BPMV) was used to silence autophagy-related gene 5 (ATG5) genes in soybean (referred to as GmATG5). Our results showed that ATG8 proteins were massively accumulated in the dark-treated leaves of the GmATG5-silenced plants relative to the vector control plants (BPMV-0), indicating that autophagy pathway is impaired in the GmATG5-silenced plants. Consistent with the impaired autophagy, an accelerated senescence phenotype was observed on the leaves of the dark-treated GmATG5-silenced plants, which was not shown on the leaves of the dark-treated BPMV-0 plants. In addition, the accumulation levels of both reactive oxygen species (ROS) and salicylic acid (SA) were significantly induced in the GmATG5-silenced plants compared with that of the vector control plants (BPMV-0), indicating an activated immunity. Accordingly, the GmATG5-silenced plants exhibited significantly enhanced resistance against Pseudomonas syringae pv. glycinea (Psg) in comparison with the BPMV-0 plants. Nevertheless, the activated immunity observed in the GmATG5-silenced plant was independent of the activation of mitogen-activated protein kinase (MAPK).

Functional study of the soybean stamen-preferentially expressed gene GmFLA22a in regulating male fertility.

China has a high dependence on soybean imports, yield increase at a faster rate is an urgent problem that need to be solved at present. The application of heterosis is one of the effective ways to significantly increase crop yield. In recent years, the development of an intelligent male sterility system based on recessive nuclear sterile genes has provided a potential solution for rapidly harnessing the heterosis in soybean. However, research on male sterility genes in soybean has been lagged behind. Based on transcriptome data of soybean floral organs in our research group, a soybean stamen-preferentially expressed gene GmFLA22a was identified. It encodes a fasciclin-like arabinogalactan protein with the FAS1 domain, and subcellular localization studies revealed that it may play roles in the endoplasmic reticulum. Take advantage of the gene editing technology, the Gmfla22a mutant was generated in this study. However, there was a significant reduction in the seed-setting rate in the mutant plants at the reproductive growth stage. The pollen viability and germination rate of Gmfla22a mutant plants showed no apparent abnormalities. Histological staining demonstrated that the release of pollen grains in the mutant plants was delayed and incomplete, which may due to the locule wall thickening in the anther development. This could be the reason of the reduced seed-setting rate in Gmfla22a mutants. In summary, our study has preliminarily revealed that GmFLA22a may be involved in regulating soybean male fertility. It provides crucial genetic materials for further uncovering its molecular function and gene resources and theoretical basis for the utilization of heterosis in soybean.

Integrating Genome-Wide Association Study, Transcriptome and Metabolome Reveal Novel QTL and Candidate Genes That Control Protein Content in Soybean.

Protein content (PC) is crucial to the nutritional quality of soybean [Glycine max (L.) Merrill]. In this study, a total of 266 accessions were used to perform a genome-wide association study (GWAS) in three tested environments. A total of 23,131 high-quality SNP markers (MAF ≥ 0.02, missing data ≤ 10%) were identified. A total of 40 association signals were significantly associated with PC. Among them, five novel quantitative trait nucleotides (QTNs) were discovered, and another 32 QTNs were found to be overlapping with the genomic regions of known quantitative trait loci (QTL) related to soybean PC. Combined with GWAS, metabolome and transcriptome sequencing, 59 differentially expressed genes (DEGs) that might control the change in protein content were identified. Meantime, four commonly upregulated differentially abundant metabolites (DAMs) and 29 commonly downregulated DAMs were found. Remarkably, the soybean gene Glyma.08G136900, which is homologous with Arabidopsis hydroxyproline-rich glycoproteins (HRGPs), may play an important role in improving the PC. Additionally, Glyma.08G136900 was divided into two main haplotype in the tested accessions. The PC of haplotype 1 was significantly lower than that of haplotype 2. The results of this study provided insights into the genetic mechanisms regulating protein content in soybean.

CRISPR/Cas9-mediated editing of GmDWF1 brassinosteroid biosynthetic gene induces dwarfism in soybean.

The study on the GmDWF1-deficient mutant dwf1 showed that GmDWF1 plays a crucial role in determining soybean plant height and yield by influencing the biosynthesis of brassinosteroids. Soybean has not adopted the Green Revolution, such as reduced height for increased planting density, which have proven beneficial for cereal crops. Our research identified the soybean genes GmDWF1a and GmDWF1b, homologous to Arabidopsis AtDWF1, and found that they are widely expressed, especially in leaves, and linked to the cellular transport system, predominantly within the endoplasmic reticulum and intracellular vesicles. These genes are essential for the synthesis of brassinosteroids (BR). Single mutants of GmDWF1a and GmDWF1b, as well as double mutants of both genes generated through CRISPR/Cas9 genome editing, exhibit a dwarf phenotype. The single-gene mutant exhibits moderate dwarfism, while the double mutant shows more pronounced dwarfism. Despite the reduced stature, all types of mutants preserve their node count. Notably, field tests have shown that the single GmDWF1a mutant produced significantly more pods than wild-type plants. Spraying exogenous brassinolide (BL) can compensate for the loss in plant height induced by the decrease in endogenous BRs. Comparing transcriptome analyses of the GmDWF1a mutant and wild-type plants revealed a significant impact on the expression of many genes that influence soybean growth. Identifying the GmDWF1a and GmDWF1b genes could aid in the development of compact, densely planted soybean varieties, potentially boosting productivity.

Knockout of miR396 genes increases seed size and yield in soybean.

Yield improvement has long been an important task for soybean breeding in the world in order to meet the increasing demand for food and animal feed. miR396 genes have been shown to negatively regulate grain size in rice, but whether miR396 family members may function in a similar manner in soybean is unknown. Here, we generated eight soybean mutants harboring different combinations of homozygous mutations in the six soybean miR396 genes through genome editing with clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated nuclease (Cas)12SF01 in the elite soybean cultivar Zhonghuang 302 (ZH302). Four triple mutants (mir396aci, mir396acd, mir396adf, and mir396cdf), two quadruple mutants (mir396abcd and mir396acfi), and two quintuple mutants (mir396abcdf and mir396bcdfi) were characterized. We found that plants of all the mir396 mutants produced larger seeds compared to ZH302 plants. Field tests showed that mir396adf and mir396cdf plants have significantly increased yield in growth zones with relatively high latitude which are suited for ZH302 and moderately increased yield in lower latitude. In contrast, mir396abcdf and mir396bcdfi plants have increased plant height and decreased yield in growth zones with relatively high latitude due to lodging issues, but they are suited for low latitude growth zones with increased yield without lodging problems. Taken together, our study demonstrated that loss-of-function of miR396 genes leads to significantly enlarged seed size and increased yield in soybean, providing valuable germplasms for breeding high-yield soybean.

An enzyme-free and ratiometric biosensor for single nucleotide polymorphism with portable electrochemical system

Portable, low-cost and efficient biosensing systems are valuable for single nucleotide polymorphisms (SNPs) analysis, particularly in resource-limited fields. Herein, a rapid, enzyme-free and ratiometric SNPs biosensor with portable electrochemical system enabling simultaneous recognition and signal amplification is proposed. To demonstrate its versatility and effectiveness, we have selected the human cancer-related KRAS gene and the soybean mosaic virus resistance Rsv3 gene as our analysis models. Our biosensor is constructed based on a combination of strand displacement reactions and hybridization chain reactions (HCR). This unique design enables to differentiate between mutant (MT) and wild (WD) sequences based on the differential release of HCR products on the electrode surfaces. Importantly, this enzyme-free SNPs biosensor concurrently achieves both SNPs recognition and signal amplification, significantly reducing the reaction time to just 40 minutes. Furthermore, we have developed a miniaturized biosensor system that incorporates screen printed gold electrodes (SPGE) and smartphones. This system is integrated with a ratiometric strategy with an limit of detection as low as 1.2 fM, enhancing accuracy and eliminating external interference. This simple SNPs biosensing system possesses high reliability, positioning it as a versatile biosensor for the determination of various SNPs mutations.

Overexpression of soybean GmDHN9 gene enhances drought resistance of transgenic Arabidopsis

ABSTRACT Soybean is one of the important oil crops and a major source of protein and lipids. Drought can cause severe soybean yields. Dehydrin protein (DHN) is a subfamily of LEA proteins that play an important role in plant responses to abiotic stresses. In this study, the soybean GmDHN9 gene was cloned and induced under a variety of abiotic stresses. Results showed that the GmDHN9 gene response was more pronounced under drought induction. Subcellular localization results indicated that the protein was localized in the cytoplasm. The role of transgenic Arabidopsis plants in drought stress response was further studied. Under drought stress, the germination rate, root length, chlorophyll, proline, relative water content, and antioxidant enzyme content of transgenic Arabidopsis thaliana transgenic genes were higher than those of wild-type plants, and transgenic plants contained less O2 −, H2O2 and MDA contents. In short, the GmDHN9 gene can regulate the homeostasis of ROS and enhance the drought resistance of plants.

Overexpression of Chalcone Isomerase-like Genes, GmCHI4A and GmCHI4B, Enhances Salt Tolerance of Cotyledon Hairy Roots and Composite Plant in Soybean (Glycine max (L.) Merr.)

Chalcone isomerase (CHI) is an important enzyme involved in the biosynthesis of flavonoids, one that is crucial in both plant defense and human health. Although many CHI genes have been previously identified, the function of CHI-like genes in soybean remains unclear. In this study, we cloned the CHI-like genes GmCHI4A and GmCHI4B (GmCHI4s) in soybean. The real-time quantitative polymerase chain reaction showed that GmCHI4s were expressed primarily in soybean root, but were also present in other tissues, including the stem, leaf, and seed with a low expression level. Overexpression of GmCHI4s was able to significantly improve some beneficial traits of the transformed hair roots of cotyledon or composite plants under salt stress conditions. Root length, root wet weight, and the underground biomass was increased, and the elevation of MDA content was inhibited under 100 mmol L−1 or 150 mmol L−1 NaCl treatment. Leaf chlorophyll content was elevated in overexpressed GmCHI4A composite plants under 150 mmol L−1 NaCl treatment. The expression levels of salt-stress-related genes GmSOD1, GmAPX1, GmSOS1, and GmNHX1 were significantly upregulated in overexpressed GmCHI4 hairy roots compared to that in empty-vector-expressed hairy roots. The above results indicated GmCHI4s’ potential action against salt stress. Furthermore, overexpression of GmCHI4A and GmCHI4B increased the total isoflavone content by six times and three times, respectively. Glycitin and glycitein levels were significantly elevated in the overexpressed GmCHI4A hairy roots, while glycitin, genistin, daidzein, and genistein were significantly increased in overexpressed GmCHI4B hairy roots. This study identified a new function of the CHI-like gene, as well as providing a new selected gene for salt tolerance and isoflavone improvement using biotechnological approaches in soybean.

Multi-omics analysis reveals novel loci and a candidate regulatory gene of unsaturated fatty acids in soybean (Glycine max (L.) Merr)

BackgroundSoybean is a major oil crop; the nutritional components of soybean oil are mainly controlled by unsaturated fatty acids (FA). Unsaturated FAs mainly include oleic acid (OA, 18:1), linoleic acid (LLA, 18:2), and linolenic acid (LNA, 18:3). The genetic architecture of unsaturated FAs in soybean seeds has not been fully elucidated, although many independent studies have been conducted. A 3 V multi-locus random single nucleotide polymorphism (SNP)-effect mixed linear model (3VmrMLM) was established to identify quantitative trait loci (QTLs) and QTL-by-environment interactions (QEIs) for complex traits.ResultsIn this study, 194 soybean accessions with 36,981 SNPs were calculated using the 3VmrMLM model. As a result, 94 quantitative trait nucleotides (QTNs) and 19 QEIs were detected using single-environment (QTN) and multi-environment (QEI) methods. Three significant QEIs, namely rs4633292, rs39216169, and rs14264702, overlapped with a significant single-environment QTN.ConclusionsFor QTNs and QEIs, further haplotype analysis of candidate genes revealed that the Glyma.03G040400 and Glyma.17G236700 genes were beneficial haplotypes that may be associated with unsaturated FAs. This result provides ideas for the identification of soybean lipid-related genes and provides insights for breeding high oil soybean varieties in the future.

Development of an NLR-ID Toolkit and Identification of Novel Disease-Resistance Genes in Soybean.

The recognition of pathogen effectors through the nucleotide-binding leucine-rich repeat receptor (NLR) family is an important component of plant immunity. In addition to typical domains such as TIR, CC, NBS, and LRR, NLR proteins also contain some atypical integrated domains (IDs), the roles of which are rarely investigated. Here, we carefully screened the soybean (Glycine max) genome and identified the IDs that appeared in the soybean TNL-like proteins. Our results show that multiple IDs (36) are widely present in soybean TNL-like proteins. A total of 27 Gm-TNL-ID genes (soybean TNL-like gene encoding ID) were cloned and their antiviral activity towards the soybean mosaic virus (SMV)/tobacco mosaic virus (TMV) was verified. Two resistance (R) genes, SRA2 (SMV resistance gene contains AAA_22 domain) and SRZ4 (SMV resistance gene contains zf-RVT domain), were identified to possess broad-spectrum resistance characteristics towards six viruses including SMV, TMV, plum pox virus (PPV), cabbage leaf curl virus (CaLCuV), barley stripe mosaic virus (BSMV), and tobacco rattle virus (TRV). The effects of Gm-TNL-IDX (the domain of the Gm-TNL-ID gene after the TN domain) on the antiviral activity of a R protein SRC7TN (we previously reported the TN domain of the soybean broad-spectrum resistance gene SRC7) were validated, and most of Gm-TNL-IDX inhibits antiviral activity mediated by SRC7TN, possibly through intramolecular interactions. Yeast-two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays showed that seven Gm-TNL-IDX interacted with SMV-component proteins. Truncation analysis on a broad-spectrum antiviral protein SRZ4 indicated that SRZ4TIR is sufficient to mediate antiviral activity against SMV. Soybean cDNA library screening on SRZ4 identified 48 interacting proteins. In summary, our results indicate that the integration of IDs in soybean is widespread and frequent. The NLR-ID toolkit we provide is expected to be valuable for elucidating the functions of atypical NLR proteins in the plant immune system and lay the foundation for the development of engineering NLR for plant-disease control in the future.