GmNAC039 and GmNAC018 activate the expression of cysteine protease genes to promote soybean nodule senescence.

Root nodules are major sources of nitrogen for soybean (Glycine max (L.) Merr.) growth, development, production, and seed quality. Symbiotic nitrogen fixation is time-limited, as the root nodule senesces during the reproductive stage of plant development, specifically during seed development. Nodule senescence is characterized by the induction of senescence-related genes, such as papain-like cysteine proteases (CYPs), which ultimately leads to the degradation of both bacteroids and plant cells. However, how nodule senescence-related genes are activated in soybean is unknown. Here, we identified 2 paralogous NAC transcription factors, GmNAC039 and GmNAC018, as master regulators of nodule senescence. Overexpression of either gene induced soybean nodule senescence with increased cell death as detected using a TUNEL assay, whereas their knockout delayed senescence and increased nitrogenase activity. Transcriptome analysis and nCUT&Tag-qPCR assays revealed that GmNAC039 directly binds to the core motif CAC(A)A and activates the expression of 4 GmCYP genes (GmCYP35, GmCYP37, GmCYP39, and GmCYP45). Similar to GmNAC039 and GmNAC018, overexpression or knockout of GmCYP genes in nodules resulted in precocious or delayed senescence, respectively. These data provide essential insights into the regulatory mechanisms of nodule senescence, in which GmNAC039 and GmNAC018 directly activate the expression of GmCYP genes to promote nodule senescence.

Dynamics of cis-regulatory sequences and transcriptional divergence of duplicated genes in soybean.

Transcriptional divergence of duplicated genes after whole genome duplication (WGD) has been described in many plant lineages and is often associated with subgenome dominance, a genome-wide mechanism. However, it is unknown what underlies the transcriptional divergence of duplicated genes in polyploid species that lack subgenome dominance. Soybean is a paleotetraploid with a WGD that occurred 5 to 13 Mya. Approximately 50% of the duplicated genes retained from this WGD exhibit transcriptional divergence. We developed accessible chromatin region (ACR) datasets from leaf, flower, and seed tissues using MNase-hypersensitivity sequencing. We validated enhancer function of several ACRs associated with known genes using CRISPR/Cas9-mediated genome editing. The ACR datasets were used to examine and correlate the transcriptional patterns of 17,111 pairs of duplicated genes in different tissues. We demonstrate that ACR dynamics are correlated with divergence of both expression level and tissue specificity of individual gene pairs. Gain or loss of flanking ACRs and mutation of cis-regulatory elements (CREs) within the ACRs can change the balance of the expression level and/or tissue specificity of the duplicated genes. Analysis of DNA sequences associated with ACRs revealed that the extensive sequence rearrangement after the WGD reshaped the CRE landscape, which appears to play a key role in the transcriptional divergence of duplicated genes in soybean. This may represent a general mechanism for transcriptional divergence of duplicated genes in polyploids that lack subgenome dominance.

Isolation and Identification of Salinity-Tolerant Rhizobia and Nodulation Phenotype Analysis in Different Soybean Germplasms.

Increasing the soybean-planting area and increasing the soybean yield per unit area are two effective solutions to improve the overall soybean yield. Northeast China has a large saline soil area, and if soybeans could be grown there with the help of isolated saline-tolerant rhizobia, the soybean cultivation area in China could be effectively expanded. In this study, soybeans were planted in soils at different latitudes in China, and four strains of rhizobia were isolated and identified from the soybean nodules. According to the latitudes of the soil-sampling sites from high to low, the four isolated strains were identified as HLNEAU1, HLNEAU2, HLNEAU3, and HLNEAU4. In this study, the isolated strains were identified for their resistances, and their acid and saline tolerances and nitrogen fixation capacities were preliminarily identified. Ten representative soybean germplasm resources in Northeast China were inoculated with these four strains, and the compatibilities of these four rhizobium strains with the soybean germplasm resources were analyzed. All four isolates were able to establish different extents of compatibility with 10 soybean resources. Hefeng 50 had good compatibility with the four isolated strains, while Suinong 14 showed the best compatibility with HLNEAU2. The isolated rhizobacteria could successfully establish symbiosis with the soybeans, but host specificity was also present. This study was a preliminary exploration of the use of salinity-tolerant rhizobacteria to help the soybean nitrogen fixation in saline soils in order to increase the soybean acreage, and it provides a valuable theoretical basis for the application of saline-tolerant rhizobia.

Natural variation in Fatty Acid 9 is a determinant of fatty acid and protein content.

Soybean is one of the most economically important crops worldwide and an important source of unsaturated fatty acids and protein for the human diet. Consumer demand for healthy fats and oils is increasing, and the global demand for vegetable oil is expected to double by 2050. Identification of key genes that regulate seed fatty acid content can facilitate molecular breeding of high-quality soybean varieties with enhanced fatty acid profiles. Here, we analysed the genetic architecture underlying variations in soybean seed fatty acid content using 547 accessions, including mainly landraces and cultivars from northeastern China. Through fatty acid profiling, genome re-sequencing, population genomics analyses, and GWAS, we identified a SEIPIN homologue at the FA9 locus as an important contributor to seed fatty acid content. Transgenic and multiomics analyses confirmed that FA9 was a key regulator of seed fatty acid content with pleiotropic effects on seed protein and seed size. We identified two major FA9 haplotypes in 1295 resequenced soybean accessions and assessed their phenotypic effects in a field planting of 424 accessions. Soybean accessions carrying FA9H2 had significantly higher total fatty acid contents and lower protein contents than those carrying FA9H1 . FA9H2 was absent in wild soybeans but present in 13% of landraces and 26% of cultivars, suggesting that it may have been selected during soybean post-domestication improvement. FA9 therefore represents a useful genetic resource for molecular breeding of high-quality soybean varieties with specific seed storage profiles.

Single-Cell RNA-Sequencing of Soybean Reveals Transcriptional Changes and Antiviral Functions of GmGSTU23 and GmGSTU24 in Response to Soybean Mosaic Virus.

Soybean mosaic virus (SMV) stands as a prominent and widespread threat to soybean (Glycine max L. Merr.), the foremost legume crop globally. Attaining a thorough comprehension of the alterations in the transcriptional network of soybeans in response to SMV infection is imperative for a profound insight into the mechanisms of viral pathogenicity and host resistance. In this investigation, we isolated 50 294 protoplasts from the newly developed leaves of soybean plants subjected to both SMV infection and mock inoculation. Subsequently, we utilized single-cell RNA sequencing (scRNA-seq) to construct the transcriptional landscape at a single-cell resolution. Nineteen distinct cell clusters were identified based on the transcriptomic profiles of scRNA-seq. The annotation of three cell types-epidermal cells, mesophyll cells, and vascular cells-was established based on the expression of orthologs to reported marker genes in Arabidopsis thaliana. The differentially expressed genes between the SMV- and mock-inoculated samples were analyzed for different cell types. Our investigation delved deeper into the tau class of glutathione S-transferases (GSTUs), known for their significant contributions to plant responses against abiotic and biotic stress. A total of 57 GSTU genes were identified by a thorough genome-wide investigation in the soybean genome G. max Wm82.a4.v1. Two specific candidates, GmGSTU23 and GmGSTU24, exhibited distinct upregulation in all three cell types in response to SMV infection, prompting their selection for further research. The transient overexpression of GmGSTU23 or GmGSTU24 in Nicotiana benthamiana resulted in the inhibition of SMV infection, indicating the antiviral function of soybean GSTU proteins.

Tandem mass tag-based proteomics reveals the antiepileptic mechanism of steroidal saponins from Anemarrhena asphodeloides in Kainic acid induced epileptic rat model.

Epilepsy (EP) is one of the most common neurological diseases in the world. Anemarrhena asphodeloides Bunge. (AA), as a typical heat-cleaning medicine, has been proven to possess the antiepileptic effect in clinical and experimental studies. Anemarrhena asphodeloides steroidal saponins (AAS) are main components. However, the therapeutic effects and underlying mechanisms of AAS against EP are not been fully elucidated. In this study, 63 steroidal saponins were discovered in AAS by UPLC-Q-TOF/MS analysis. Pharmacological and behavioral analysis demonstrated that AAS could significantly lower the Racine classification and reduce the frequency of generalized spike rhythm the rate of tetanic seizures in kainic acid-induced epileptic rats. Hematoxylin and eosin and Nissl staining-indicated AAS could significantly improve hippocampal injury and neuron loss in epileptic rats. TMT proteomic analysis discovered 26 different expressed proteins (DEPs), which were identified as the rescue proteins. After bioinformatic analysis, Heat Shock Protein 90 Alpha Family Class B Member 1 (Hsp90ab1) and Tyrosine 3-Monooxygenase (Ywhab) were screened as key DEPs and verified by western blotting. AAS could significantly inhibited the up-regulation of Hsp90ab1 and Ywhab in EP rats; these two proteins might be the key targets of AAS in treating EP.

Genetic Analysis of Soybean Flower Size Phenotypes Based on Computer Vision and Genome-Wide Association Studies.

The dimensions of organs such as flowers, leaves, and seeds are governed by processes of cellular proliferation and expansion. In soybeans, the dimensions of these organs exhibit a strong correlation with crop yield, quality, and other phenotypic traits. Nevertheless, there exists a scarcity of research concerning the regulatory genes influencing flower size, particularly within the soybean species. In this study, 309 samples of 3 soybean types (123 cultivar, 90 landrace, and 96 wild) were re-sequenced. The microscopic phenotype of soybean flower organs was photographed using a three-eye microscope, and the phenotypic data were extracted by means of computer vision. Pearson correlation analysis was employed to assess the relationship between petal and seed phenotypes, revealing a strong correlation between the sizes of these two organs. Through GWASs, SNP loci significantly associated with flower organ size were identified. Subsequently, haplotype analysis was conducted to screen for upstream and downstream genes of these loci, thereby identifying potential candidate genes. In total, 77 significant SNPs associated with vexil petals, 562 significant SNPs associated with wing petals, and 34 significant SNPs associated with keel petals were found. Candidate genes were screened by candidate sites, and haplotype analysis was performed on the candidate genes. Finally, the present investigation yielded 25 and 10 genes of notable significance through haplotype analysis in the vexil and wing regions, respectively. Notably, Glyma.07G234200, previously documented for its high expression across various plant organs, including flowers, pods, leaves, roots, and seeds, was among these identified genes. The research contributes novel insights to soybean breeding endeavors, particularly in the exploration of genes governing organ development, the selection of field materials, and the enhancement of crop yield. It played a role in the process of material selection during the growth period and further accelerated the process of soybean breeding material selection.

QTL Mapping and Data Mining to Identify Genes Associated with Soybean Epicotyl Length Using Cultivated Soybean and Wild Soybean.

Soybean (Glycine max) plants first emerged in China, and they have since been established as an economically important oil crop and a major source of daily protein for individuals throughout the world. Seed emergence height is the first factor that ensures seedling adaptability to field management practices, and it is closely related to epicotyl length. In the present study, the Suinong 14 and ZYD00006 soybean lines were used as parents to construct chromosome segment substitution lines (CSSLs) for quantitative trait loci (QTL) identification. Seven QTLs were identified using two years of epicotyl length measurement data. The insertion region of the ZYD00006 fragment was identified through whole genome resequencing, with candidate gene screening and validation being performed through RNA-Seq and qPCR, and Glyma.08G142400 was ultimately selected as an epicotyl length-related gene. Through combined analyses of phenotypic data from the study population, Glyma.08G142400 expression was found to be elevated in those varieties exhibiting longer epicotyl length. Haplotype data analyses revealed that epicotyl data were consistent with haplotype typing. In summary, the QTLs found to be associated with the epicotyl length identified herein provide a valuable foundation for future molecular marker-assisted breeding efforts aimed at improving soybean emergence height in the field, with the Glyma.08G142400 gene serving as a regulator of epicotyl length, offering new insight into the mechanisms that govern epicotyl development.

Multi-Omics Analysis of a Chromosome Segment Substitution Line Reveals a New Regulation Network for Soybean Seed Storage Profile.

Soybean, a major source of oil and protein, has seen an annual increase in consumption when used in soybean-derived products and the broadening of its cultivation range. The demand for soybean necessitates a better understanding of the regulatory networks driving storage protein accumulation and oil biosynthesis to broaden its positive impact on human health. In this study, we selected a chromosome segment substitution line (CSSL) with high protein and low oil contents to investigate the underlying effect of donor introgression on seed storage through multi-omics analysis. In total, 1479 differentially expressed genes (DEGs), 82 differentially expressed proteins (DEPs), and 34 differentially expressed metabolites (DEMs) were identified in the CSSL compared to the recurrent parent. Based on Gene Ontology (GO) term analysis and the Kyoto Encyclopedia of Genes and Genomes enrichment (KEGG), integrated analysis indicated that 31 DEGs, 24 DEPs, and 13 DEMs were related to seed storage functionality. Integrated analysis further showed a significant decrease in the contents of the seed storage lipids LysoPG 16:0 and LysoPC 18:4 as well as an increase in the contents of organic acids such as L-malic acid. Taken together, these results offer new insights into the molecular mechanisms of seed storage and provide guidance for the molecular breeding of new favorable soybean varieties.

A highly efficient soybean transformation system using GRF3-GIF1 chimeric protein.

Expression of GRF3-GIF1 chimera significantly enhanced regeneration and transformation efficiency in soybean, increasing the number of transformable cultivars. Moreover, GmGRF3-GIF1 can be combined with CRISPR/Cas9 for highly effective gene editing.