SoyMD: a platform combining multi-omics data with various tools for soybean research and breeding

Comprehensive genome-wide analysis of the GmFRIGIDA gene family in soybean: identification, characterization, and expression dynamics

Background

Frigida (FRI) genes are crucial for regulating flowering time in plants. While the biological importance of the Frigida-like (FRL) gene family has been recognized in Arabidopsis, a systematic analysis of these genes in soybean is lacking. Characterizing FRL genes in soybean will help uncover their roles in flowering regulation, offering valuable insights for improving soybean adaptation.

Results

In this study, we identified 16 Frigida genes in soybean, naming them based on their relationship to the FRL genes in Arabidopsis thaliana. These genes are unevenly distributed across thirteen chromosomes. Phylogenetic analysis categorizes Frigida-like proteins from Arabidopsis, soybean, and rice into four distinct subfamilies (I-IV). Our findings indicate that eight GmFRLs arose from whole-genome duplication (WGD) events, alongside two tandem duplication events. Gene structure analysis confirmed that all GmFRL members contain Frigida domains. Additionally, promoter analysis revealed numerous cis-acting elements related to photoperiodic response, suggesting their significant role in soybean's light response mechanisms. RNA-seq data demonstrated variable expression levels of GmFRL genes across tissues, including flower, leaf, pod, and seed, and other tissues, while subcellular localization and qPCR analyses further support their vital role in light responsiveness in soybean.

Conclusion

In summary, our comprehensive analysis offers valuable insights into the evolution and potential functions of GmFRL genes, emphasizing their significance in photoperiodic responses and establishing a foundation for further research on the GmFRL family.

Seed-specific expression of phosphatidate phosphohydrolases increases soybean oil content and seed weight

BACKGROUND

Soybean is a major oil crop and a primary protein source for livestock, and soybean oil is the most common input for biodiesel. Identifying genes that enhance soybean yield and oil content is crucial for breeding programs. Phosphatidic acid (PA) phosphohydrolase (PAH), which dephosphorylates PA to diacylglycerol (DAG), plays a critical role in lipid synthesis, and yet their potential in improving agronomic traits of oil crops remains unexplored.

RESULTS

This study shows that seed-specific expression of AtPAH1/2 enhances PA turnover into DAG and triacylglycerol (TAG) accumulation in soybean seeds. PAH overexpression upregulated the expression of DAG acyltransferase (DGAT) but suppressed phospholipid: DAG acyltransferase (PDAT). In addition, seed-specific expression of AtPAH1/2 increases soybean seed size and weight. Furthermore, analysis of the variation of the soybean PAHs in 4414 soybean accessions indicated that the advantageous effects of GmPAHs on oil content and seed weight were selected during domestication.

CONCLUSION

These findings suggest that targeting PAHs represents a promising strategy for enhancing soybean seed oil content and yield in current cultivars and landraces soybeans.

The transcription factors GmVOZ1A and GmWRI1a synergistically regulate oil biosynthesis in soybean

Soybean (Glycine max [L.] Merr.) is a major oil-producing crop worldwide. Although several related proteins regulating soybean oil accumulation have been reported, little is known about the regulatory mechanisms. In this study, we characterized vascular plant one-zinc-finger 1A (GmVOZ1A) that interacts with WRINKLED 1a (GmWRI1a) using yeast 2-hybrid library screening. The GmVOZ1A-GmWRI1a interaction was further verified by protein-protein interaction assays in vivo and in vitro. GmVOZ1A enhanced the seed fatty acid and oil contents by regulating genes involved in lipid biosynthesis. Conversely, a loss-of-function mutation in GmVOZ1A resulted in a reduction in triacylglycerol (TAG) content in soybean. Protein-DNA interaction assays revealed that GmVOZ1A and GmWRI1a cooperate to upregulate the expression level of acyl-coenzyme A-binding protein 6a (GmACBP6a) and promote the accumulation of TAG. In addition, GmACBP6a overexpression promoted seed fatty acid and oil contents, as well as increased seed size and 100-seed weight. Taken together, these findings indicate that the transcription factor GmVOZ1A regulates soybean oil synthesis and cooperates with GmWRI1a to upregulate GmACBP6a expression and oil biosynthesis in soybean. The results lay a foundation for a comprehensive understanding of the regulatory mechanisms underlying soybean oil biosynthesis and will contribute to improving soybean oil production through molecular breeding approaches.

Soybean2035: A decadal vision for soybean functional genomics and breeding

Soybean, the fourth most important crop in the world, uniquely serves as a source of both plant oil and plant protein for the world's food and animal feed. Although soybean production has increased approximately 13-fold over the past 60 years, the continually growing global population necessitates further increases in soybean production. In the past, especially in the last decade, significant progress has been made in both functional genomics and molecular breeding. However, many more challenges should be overcome to meet the anticipated future demand. Here, we summarize past achievements in the areas of soybean omics, functional genomics, and molecular breeding. Furthermore, we analyze trends in these areas, including shortages and challenges, and propose new directions, potential approaches, and possible outputs toward 2035. Our views and perspectives provide insight into accelerating the development of elite soybean varieties to meet the increasing demands of soybean production.

SoyOD: An Integrated Soybean Multi-omics Database for Mining Genes and Biological Research

Soybean is a globally important crop for food, feed, oil, and nitrogen fixation. A variety of multi-omics studies have been carried out, generating datasets ranging from genotype to phenotype. In order to efficiently utilize these data for basic and applied research, a soybean multi-omics database with extensive data coverage and comprehensive data analysis tools was established. The Soybean Omics Database (SoyOD) integrates important new datasets with existing public datasets to form the most comprehensive collection of soybean multi-omics information. Compared to existing soybean databases, SoyOD incorporates an extensive collection of novel data derived from the deep-sequencing of 984 germplasms, 162 novel transcriptomic datasets from seeds at different developmental stages, 53 phenotypic datasets, and more than 2500 phenotypic images. In addition, SoyOD integrates existing data resources, including 59 assembled genomes, genetic variation data from 3904 soybean accessions, 225 sets of phenotypic data, and 1097 transcriptomic sequences covering 507 different tissues and treatment conditions. Moreover, SoyOD can be used to mine candidate genes for important agronomic traits, as shown in a case study on plant height. Additionally, powerful analytical and easy-to-use toolkits enable users to easily access the available multi-omics datasets, and to rapidly search genotypic and phenotypic data in a particular germplasm. The novelty, comprehensiveness, and user-friendly features of SoyOD make it a valuable resource for soybean molecular breeding and biological research. SoyOD is publicly accessible at https://bis.zju.edu.cn/soyod.

Identification of superior haplotypes and candidate gene for seed size-related traits in soybean (Glycine max L.)

Seed size is an economically important trait that directly determines the seed yield in soybean. In the current investigation, we used an integrated strategy of linkage mapping, association mapping, haplotype analysis and candidate gene analysis to determine the genetic makeup of four seed size-related traits viz., 100-seed weight (HSW), seed area (SA), seed length (SL), and seed width (SW) in soybean. Linkage mapping identified a total of 23 quantitative trait loci (QTL) associated with four seed size-related traits in the F2 population; among them, 17 were detected as novel QTLs, whereas the remaining six viz., qHSW3-1, qHSW4-1, qHSW18-1, qHSW19-1, qSL4-1 and qSW6-1 have been previously identified. Six out of 23 QTLs were major possessing phenotypic variation explained (PVE) ≥ 10%. Besides, the four QTL Clusters/QTL Hotspots harboring multiple QTLs for different seed size-related traits were identified on Chr.04, Chr.16, Chr.19 and Chr.20. Genome-wide association study (GWAS) identified a total of 62 SNPs significantly associated with the four seed size-related traits. Interestingly, the QTL viz., qHSW18-1 was identified by both linkage mapping and GWAS, and was regarded as the most stable loci regulating HSW in soybean. In-silico, sequencing and qRT-PCR analysis identified the Glyma.18G242400 as the most potential candidate gene underlying the qHSW18-1 for regulating HSW. Moreover, three haplotype blocks viz., Hap2, Hap6A and Hap6B were identified for the SW trait, and one haplotype was identified within the Glyma.18G242400 for the HSW. These four haplotypes harbor three to seven haplotype alleles across the association mapping panel of 350 soybean accessions, regulating the seed size from lowest to highest through intermediate phenotypes. Hence, the outcome of the current investigation can be utilized as a potential genetic and genomic resource for breeding the improved seed size in soybean.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01525-1.

Genome-wide association study of sucrose content in vegetable soybean

BACKGROUND

Vegetable soybean is an important legume vegetable. High sucrose content is a significant quality characteristic of vegetable soybean that influences consumers' taste. However, the genetic basis of sucrose content in vegetable soybean is currently unclear.

RESULTS

In this study, the genome-wide association study (GWAS) of sucrose content in vegetable soybean was performed using Chinese soybean mini-core collection. The results showed a wide genetic variation for the sucrose content in the mini-core collection. The sucrose content of genotypes from HHR (Huanghuai region) and SR (Southern region) was higher than that of genotypes from NER (Northeast region) and NR (Northern region). Furthermore, 82,187 high quality SNPs (Single nucleotide polymorphism) were used for GWAS of sucrose content. Based on SNPs detected in multiple environments, the chromosome 8 19,496,314-19,698,413 bp interval was identified as the candidate interval. And Glyma.08g234100 was most likely to affect the sucrose content of vegetable soybean seeds.

CONCLUSIONS

This study has created new details to be used for breeding for high sucrose content in vegetable soybean.

Genome-wide association study revealed some new candidate genes associated with flowering and maturity time of soybean in Central and West Siberian regions of Russia

The duration of flowering and maturity is an important agricultural trait determining the suitability of a variety for cultivation in the target region. In the present study, we used genome-wide association analysis (GWAS) to search for loci associated with soybean flowering and maturity in the Central and West Siberian regions of Russia. A field experiment was conducted in 2021/2022 at two locations (Orel and Novosibirsk). A germplasm collection of 180 accessions was genotyped using SoySNP50K Illumina Infinium Bead-Chip. From the initial collection, we selected 129 unrelated accessions and conducted GWAS on this dataset using two multi-locus models: FarmCPU and BLINK. As a result, we identified 13 loci previously reported to be associated with duration of soybean development, and 17 new loci. 33 candidate genes were detected in these loci using analysis of co-expression, gene ontology, and literature data, with the best candidates being Glyma.03G177500, Glyma.13G177400, and Glyma.06G213100. These candidate genes code the Arabidopis orthologs TOE1 (TARGET OF EAT 1), SPL3 (SQUAMOSA PROMOTER BINDING PROTEIN LIKE 3), the DELLA protein, respectively. In these three genes, we found haplotypes which may be associated with the length of soybean flowering and maturity, providing soybean adaptation to a northern latitudes.

Harnessing Multi-Omics Strategies and Bioinformatics Innovations for Advancing Soybean Improvement: A Comprehensive Review

Soybean improvement has entered a new era with the advent of multi-omics strategies and bioinformatics innovations, enabling more precise and efficient breeding practices. This comprehensive review examines the application of multi-omics approaches in soybean-encompassing genomics, transcriptomics, proteomics, metabolomics, epigenomics, and phenomics. We first explore pre-breeding and genomic selection as tools that have laid the groundwork for advanced trait improvement. Subsequently, we dig into the specific contributions of each -omics field, highlighting how bioinformatics tools and resources have facilitated the generation and integration of multifaceted data. The review emphasizes the power of integrating multi-omics datasets to elucidate complex traits and drive the development of superior soybean cultivars. Emerging trends, including novel computational techniques and high-throughput technologies, are discussed in the context of their potential to revolutionize soybean breeding. Finally, we address the challenges associated with multi-omics integration and propose future directions to overcome these hurdles, aiming to accelerate the pace of soybean improvement. This review serves as a crucial resource for researchers and breeders seeking to leverage multi-omics strategies for enhanced soybean productivity and resilience.

Bioinformatics Identification and Expression Analysis of Acetyl-CoA Carboxylase Reveal Its Role in Isoflavone Accumulation during Soybean Seed Development

Isoflavones belong to the class of flavonoid compounds, which are important secondary metabolites that play a crucial role in plant development and defense. Acetyl-CoA carboxylase (ACCase) is a biotin-dependent enzyme that catalyzes the conversion of Acetyl-CoA into Malonyl-CoA in plants. It is a key enzyme in fatty acid synthesis and also catalyzes the production of various secondary metabolites. However, information on the ACC gene family in the soybean (Glycine max L. Merr.) genome and the specific members involved in isoflavone biosynthesis is still lacking. In this study, we identified 20 ACC family genes (GmACCs) from the soybean genome and further characterized their evolutionary relationships and expression patterns. Phylogenetic analysis showed that the GmACCs could be divided into five groups, and the gene structures within the same groups were highly conserved, indicating that they had similar functions. The GmACCs were randomly distributed across 12 chromosomes, and collinearity analysis suggested that many GmACCs originated from tandem and segmental duplications, with these genes being under purifying selection. In addition, gene expression pattern analysis indicated that there was functional divergence among GmACCs in different tissues. The GmACCs reached their peak expression levels during the early or middle stages of seed development. Based on the transcriptome and isoflavone content data, a weighted gene co-expression network was constructed, and three candidate genes (Glyma.06G105900, Glyma.13G363500, and Glyma.13G057400) that may positively regulate isoflavone content were identified. These results provide valuable information for the further functional characterization and application of GmACCs in isoflavone biosynthesis in soybean.

MS2/GmAMS1 encodes a bHLH transcription factor important for tapetum degeneration in soybean

KEY MESSAGE

GmAMS1 is the only functional AMS and works with GmTDF1-1 and GmMS3 to orchestrate the tapetum degeneration in soybean. Heterosis could significantly increase the production of major crops as well as soybean [Glycine max (L.) Merr.]. Stable male-sterile/female-fertile mutants including ms2 are useful resources to apply in soybean hybrid production. Here, we identified the detailed mutated sites of two classic mutants ms2 (Eldorado) and ms2 (Ames) in MS2/GmAMS1 via the high-throughput sequencing method. Subsequently, we verified that GmAMS1, a bHLH transcription factor, is the only functional AMS member in soybean through the complementary experiment in Arabidopsis; and elucidated the dysfunction of its homolog GmAMS2 is caused by the premature stop codon in the gene's coding sequence. Further qRT-PCR analysis and protein-protein interaction assays indicated GmAMS1 is required for expressing downstream members in the putative DYT1-TDF1-AMS-MYB80/MYB103/MS188-MS1 cascade module, and might regulate the upstream members in a feedback mechanism. GmAMS1 could interact with GmTDF1-1 and GmMS3 via different region, which contributes to dissect the mechanism in the tapetum degeneration process. Additionally, as a core member in the conserved cascade module controlling the tapetum development and degeneration, AMS is conservatively present in all land plant lineages, implying that AMS-mediated signaling pathway has been established before land plants diverged, which provides further insight into the tapetal evolution.

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.

Genome-Wide Studies of FH Family Members in Soybean (Glycine max) and Their Responses under Abiotic Stresses

Formins or formin homology 2 (FH2) proteins, evolutionarily conserved multi-domain proteins in eukaryotes, serve as pivotal actin organizers, orchestrating the structure and dynamics of the actin cytoskeleton. However, a comprehensive investigation into the formin family and their plausible involvement in abiotic stress remains undocumented in soybean (Glycine max). In the current study, 34 soybean FH (GmFH)family members were discerned, their genomic distribution spanning the twenty chromosomes in a non-uniform pattern. Evolutionary analysis of the FH gene family across plant species delineated five discernible groups (Group I to V) and displayed a closer evolutionary relationship within Glycine soja, Glycine max, and Arabidopsis thaliana. Analysis of the gene structure of GmFH unveiled variable sequence lengths and substantial diversity in conserved motifs. Structural prediction in the promoter regions of GmFH gene suggested a large set of cis-acting elements associated with hormone signaling, plant growth and development, and stress responses. The investigation of the syntenic relationship revealed a greater convergence of GmFH genes with dicots, indicating a close evolutionary affinity. Transcriptome data unveiled distinctive expression patterns of several GmFH genes across diverse plant tissues and developmental stages, underscoring a spatiotemporal regulatory framework governing the transcriptional dynamics of GmFH gene. Gene expression and qRT-PCR analysis identified many GmFH genes with a dynamic pattern in response to abiotic stresses, revealing their potential roles in regulating plant stress adaptation. Additionally, protein interaction analysis highlighted an intricate web of interactions among diverse GmFH proteins. These findings collectively underscore a novel biological function of GmFH proteins in facilitating stress adaptation in soybeans.

The 2024 Nucleic Acids Research database issue and the online molecular biology database collection

The 2024 Nucleic Acids Research database issue contains 180 papers from across biology and neighbouring disciplines. There are 90 papers reporting on new databases and 83 updates from resources previously published in the Issue. Updates from databases most recently published elsewhere account for a further seven. Nucleic acid databases include the new NAKB for structural information and updates from Genbank, ENA, GEO, Tarbase and JASPAR. The Issue's Breakthrough Article concerns NMPFamsDB for novel prokaryotic protein families and the AlphaFold Protein Structure Database has an important update. Metabolism is covered by updates from Reactome, Wikipathways and Metabolights. Microbes are covered by RefSeq, UNITE, SPIRE and P10K; viruses by ViralZone and PhageScope. Medically-oriented databases include the familiar COSMIC, Drugbank and TTD. Genomics-related resources include Ensembl, UCSC Genome Browser and Monarch. New arrivals cover plant imaging (OPIA and PlantPAD) and crop plants (SoyMD, TCOD and CropGS-Hub). The entire Database Issue is freely available online on the Nucleic Acids Research website (https://academic.oup.com/nar). Over the last year the NAR online Molecular Biology Database Collection has been updated, reviewing 1060 entries, adding 97 new resources and eliminating 388 discontinued URLs bringing the current total to 1959 databases. It is available at http://www.oxfordjournals.org/nar/database/c/.