Development of a versatile resource for post-genomic research through consolidating and characterizing 1500 diverse wild and cultivated soybean genomes

Whole-genome sequencing analysis of soybean diversity across different countries and selection signature of Korean soybean accession

Soybean is an important agricultural crop known for its high protein and oil content, contributing to essential nutritional and health benefits for humans. Domesticated in China over 5,000 years ago, soybean has since adapted to diverse environments and spread worldwide. This study aimed to investigate the genomic characteristics and population structures of 2,317 publicly available soybean whole-genome sequences from diverse geographical regions, including China, Korea, Japan, Europe, North America, and South America. We used large-scale whole-genome sequencing data to perform high-resolution analyses to reveal the genetic characteristics of soybean accessions. Soybean accessions from China and Korea exhibited landrace characteristics, indicating higher genetic diversity and adaptation to local environments. On the other hand, soybean accessions from Japan, the European Union, and South America were found to have low genetic diversity due to artificial selection and breeding for agronomic traits. We also identified key variants and genes associated with the ability to adapt to different environments. In Korean soybean accessions, we observed strong selection signals for isoflavone synthesis, an adaptive trait critical for improving soybean adaptability, survival, and reproductive success by mitigating environmental stress. Identifying specific genomic regions showing unique patterns of selective sweeps for genes such as HIDH, CYP73A11, IFS1, and CYP81E11 associated with isoflavone synthesis provided valuable insights into potential adaptation mechanisms. Our research has significantly improved our understanding of soybean diversity at the genetic level. We have identified key genetic variants and genes influencing adaptability, laying the foundation for future advances in genomics-based breeding programs and crop improvement efforts.

A Novel 10-Base Pair Deletion in the First Exon of GmHY2a Promotes Hypocotyl Elongation, Induces Early Maturation, and Impairs Photosynthetic Performance in Soybean

Plants photoreceptors perceive changes in light quality and intensity and thereby regulate plant vegetative growth and reproductive development. By screening a γ irradiation-induced mutant library of the soybean (Glycine max) cultivar "Dongsheng 7", we identified Gmeny, a mutant with elongated nodes, yellowed leaves, decreased chlorophyll contents, altered photosynthetic performance, and early maturation. An analysis of bulked DNA and RNA data sampled from a population segregating for Gmeny, using the BVF-IGV pipeline established in our laboratory, identified a 10 bp deletion in the first exon of the candidate gene Glyma.02G304700. The causative mutation was verified by a variation analysis of over 500 genes in the candidate gene region and an association analysis, performed using two populations segregating for Gmeny. Glyma.02G304700 (GmHY2a) is a homolog of AtHY2a in Arabidopsis thaliana, which encodes a PΦB synthase involved in the biosynthesis of phytochrome. A transcriptome analysis of Gmeny using the Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed changes in multiple functional pathways, including photosynthesis, gibberellic acid (GA) signaling, and flowering time, which may explain the observed mutant phenotypes. Further studies on the function of GmHY2a and its homologs will help us to understand its profound regulatory effects on photosynthesis, photomorphogenesis, and flowering time.

Soybean type-B response regulator GmRR1 mediates phosphorus uptake and yield by modifying root architecture

Phosphorus (P) plays a pivotal role in plant growth and development. Low P stress can greatly hamper plant growth. Here, we identified a QTL (named QPH-9-1), which is associated with P efficiency across multiple environments through linkage analysis and genome-wide association study. Furthermore, we successfully cloned the underlying soybean (Glycine max) gene GmRR1 (a soybean type-B Response Regulator 1) that encodes a type-B response regulator protein. Knockout of GmRR1 resulted in a substantial increase in plant height, biomass, P uptake efficiency, and yield-related traits due to the modification of root structure. In contrast, overexpression of GmRR1 in plants resulted in a decrease in these phenotypes. Further analysis revealed that knockout of GmRR1 substantially increased the levels of auxin and ethylene in roots, thereby promoting root hair formation and growth by promoting the formation of root hair primordium and lengthening the root apical meristem. Yeast two-hybrid, bimolecular fluorescence complementation, and dual-luciferase assays demonstrated an interaction between GmRR1 and Histidine-containing Phosphotransmitter protein 1. Expression analysis suggested that these proteins coparticipated in response to low P stress. Analysis of genomic sequences showed that GmRR1 underwent a selection during soybean domestication. Taken together, this study provides further insights into how plants respond to low P stress by modifying root architecture through phytohormone pathways.

Anti-Diabetic Effect of Lactobacillus Paracasei Isolated from Malaysian Water Kefir Grains

The prevalence of type 2 diabetes mellitus (T2DM) is alarming because it is always linked to the increase in chronic diseases, mortality, and socioeconomic burden. Water kefir has a wide range of functional and probiotic characteristics attributed to the microorganisms present in the kefir grains. The present study aims to evaluate the in vivo anti-diabetic potential of the isolated Lactobacillus paracasei from Malaysian water kefir grains (MWKG) which was reported to have excellent probiotic properties and high antioxidant activities as reported previously. High-fat diet/streptozotocin (HFD/STZ) induction was used to obtain a T2DM model followed by treatment with the isolated L. paracasei from MWKG. The levels of glucose, insulin, and in vivo liver antioxidants were quantified after 14 weeks. Gene expression analysis of the liver was also carried out using microarray analysis, and several genes were selected for validation using quantitative real-time PCR. Insulin tolerance test demonstrated that the L. paracasei isolated from the MWKG alleviated T2DM by improving the area under the curve of the insulin tolerance test whereby low-dose and high-dose concentrations treated groups showed 2424.50 ± 437.02 mmol/L·min and 2017.50 ± 347.09 mmol/L·min, respectively, compared to untreated diabetic mice which was 3884.50 ± 39.36 mmol/L·min. Additionally, treatment with the isolated L. paracasei from MWKG regulated the expression of several genes related to glucose homeostasis and lipid metabolism in diabetic mice. These results suggested that the isolated L. paracasei from MWKG could be a potential dietary supplement for T2DM.

SoyDBean: a database for SNPs reconciliation by multiple versions of soybean reference genomes

Due to the development of sequence technology and decreased cost, many whole genome sequences have been obtained. As a result, extensive genetic variations have been discovered from many populations and germplasms to understand the genetic diversity of soybean (Glycine max [L.] Merr.). However, assessing the quality of variation is essential because the published variants were collected using different bioinformatic methods and parameters. Furthermore, despite the enhanced genome contiguity and more efficient filling of "N" stretches in the new reference genome, there remains a dearth of endeavors to verify the caliber of variations present in it. The primary goal of this research was to discern a dependable set of SNPs that can withstand reconciliation across multiple reference genomes. Additionally, the investigation aimed to reconfirm the variations through the utilization of numerous whole genome sequencing data obtained from publicly available databases. Based on the result, we created datasets that comprised the thoroughly verified SNP coordinates between the reference assemblies. The resulting "SoyDBean" database is now publicly accessible through the following URL: http://soydbean.plantprofile.net/ .

GmFtsH25 overexpression increases soybean seed yield by enhancing photosynthesis and photosynthates

Increasing plant photosynthetic capacity is a promising approach to boost yields, but it is particularly challenging in C₃ crops, such as soybean (Glycine max (L.) Merr.). Here, we identified GmFtsH25, encoding a member of the filamentation temperature-sensitive protein H protease family, as a major gene involved in soybean photosynthesis, using linkage mapping and a genome-wide association study. Overexpressing GmFtsH25 resulted in more grana thylakoid stacks in chloroplasts and increased photosynthetic efficiency and starch content, while knocking out GmFtsH25 produced the opposite phenotypes. GmFtsH25 interacted with photosystem I light harvesting complex 2 (GmLHCa2), and this interaction may contribute to the observed enhanced photosynthesis. GmFtsH25 overexpression lines had superior yield traits, such as yield per plant, compared to the wild type and knockout lines. Additionally, we identified an elite haplotype of GmFtsH25, generated by natural mutations, which appears to have been selected during soybean domestication. Our study sheds light on the molecular mechanism by which GmFtsH25 modulates photosynthesis and provides a promising strategy for improving the yields of soybean and other crops.

GmEIL4 enhances soybean (Glycine max) phosphorus efficiency by improving root system development

Phosphorus (P) deficiency seriously affects plant growth and development and ultimately limits the quality and yield of crops. Here, a new P efficiency-related major quantitative trait locus gene, GmEIL4 (encoding an ethylene-insensitive 3-like 1 protein), was cloned at qP2, which was identified by linkage analysis and genome-wide association study across four environments. Overexpressing GmEIL4 significantly improved the P uptake efficiency by increasing the number, length and surface area of lateral roots of hairy roots in transgenic soybeans, while interfering with GmEIL4 resulted in poor root phenotypic characteristics compared with the control plants under low P conditions. Interestingly, we found that GmEIL4 interacted with EIN3-binding F box protein 1 (GmEBF1), which may regulate the root response to low P stress. We conclude that the expression of GmEIL4 was induced by low-P stress and that overexpressing GmEIL4 improved P accumulation by regulating root elongation and architecture. Analysis of allele variation of GmEIL4 in 894 soybean accessions suggested that GmEIL4 is undergoing artificial selection during soybean evolution, which will benefit soybean production. Together, this study further elucidates how plants respond to low P stress by modifying root structure and provides insight into the great potential of GmEIL4 in crop P-efficient breeding.

Fitness and Rhizobacteria of F2, F3 Hybrids of Herbicide-Tolerant Transgenic Soybean and Wild Soybean

The introduction of herbicide-tolerant (HT) transgenic soybeans (Glycine max (L.) Merr.) into farming systems raises great concern that transgenes may flow to endemic wild soybeans (Glycine soja Sieb. et Zucc.) via pollen, which may increase the ecological risks by increasing the fitness of hybrids under certain conditions and threaten the genetic diversity of wild soybean populations. In order to demonstrate the potential risk of gene flow from the HT soybean to the wild soybean, the fitness of F2 and F3 hybrids obtained from two wild soybean populations (HLJHRB-1, JSCZ) collected from China and the HT soybean was measured under farmland and wasteland soil conditions, as well as with or without weed competition. Compared with their wild progenitors, the F2 and F3 hybrids of HLJHRB-1 displayed a higher emergence rate, higher aboveground dry biomass, more pods and filled-seed plants, as well as better composite fitness under four planting conditions. The F2 and F3 hybrids of JSCZ also displayed a higher emergence rate, higher aboveground dry biomass, more pods, and more filled seeds per plant under mixed planting, whereas these characteristics were lower under pure planting conditions in wasteland and farmland soil. Therefore, the composite fitness of JSCZ hybrids was higher or lower depending on the planting conditions. Furthermore, the soil microbial communities of the F3 of HLJHRB-1, JSCZ, and the wild soybean were investigated with 16S rDNA sequencing, which showed that low alpha diversity of rhizobacteria was relative to high fitness, and Rhizobium played an important role in promoting F3 plant growth.

POWR1 is a domestication gene pleiotropically regulating seed quality and yield in soybean

Seed protein, oil content and yield are highly correlated agronomically important traits that essentially account for the economic value of soybean. The underlying molecular mechanisms and selection of these correlated seed traits during soybean domestication are, however, less known. Here, we demonstrate that a CCT gene, POWR1, underlies a large-effect protein/oil QTL. A causative TE insertion truncates its CCT domain and substantially increases seed oil content, weight, and yield while decreasing protein content. POWR1 pleiotropically controls these traits likely through regulating seed nutrient transport and lipid metabolism genes. POWR1 is also a domestication gene. We hypothesize that the TE insertion allele is exclusively fixed in cultivated soybean due to selection for larger seeds during domestication, which significantly contributes to shaping soybean with increased yield/seed weight/oil but reduced protein content. This study provides insights into soybean domestication and is significant in improving seed quality and yield in soybean and other crop species.

Multi-Omics Techniques for Soybean Molecular Breeding

Soybean is a major crop that provides essential protein and oil for food and feed. Since its origin in China over 5000 years ago, soybean has spread throughout the world, becoming the second most important vegetable oil crop and the primary source of plant protein for global consumption. From early domestication and artificial selection through hybridization and ultimately molecular breeding, the history of soybean breeding parallels major advances in plant science throughout the centuries. Now, rapid progress in plant omics is ushering in a new era of precision design breeding, exemplified by the engineering of elite soybean varieties with specific oil compositions to meet various end-use targets. The assembly of soybean reference genomes, made possible by the development of genome sequencing technology and bioinformatics over the past 20 years, was a great step forward in soybean research. It facilitated advances in soybean transcriptomics, proteomics, metabolomics, and phenomics, all of which paved the way for an integrated approach to molecular breeding in soybean. In this review, we summarize the latest progress in omics research, highlight novel findings made possible by omics techniques, note current drawbacks and areas for further research, and suggest that an efficient multi-omics approach may accelerate soybean breeding in the future. This review will be of interest not only to soybean breeders but also to researchers interested in the use of cutting-edge omics technologies for crop research and improvement.