Simultaneous changes in seed size, oil content and protein content driven by selection of SWEET homologues during soybean domestication

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.

Unlocking soybean potential: genetic resources and omics for breeding

Soybean (Glycine max) is a vital foundation of global food security, providing a primary source of high-quality protein and oil for human consumption and animal feed. The rising global population has significantly increased the demand for soybeans, emphasizing the urgency of developing high-yield, stress-tolerant, and nutritionally superior cultivars. The extensive collection of soybean germplasm resources-including wild relatives, landraces, and cultivars-represents a valuable reservoir of genetic diversity critical for breeding advancements. Recent breakthroughs in genomic technologies, particularly high-throughput sequencing and multi-omics approaches, have revolutionized the identification of key genes associated with essential agronomic traits within these resources. These innovations enable precise and strategic utilization of genetic diversity, empowering breeders to integrate traits that improve yield potential, resilience to biotic and abiotic stresses, and nutritional quality. This review highlights the critical role of genetic resources and omics-driven innovations in soybean breeding. It also offers insights into strategies for accelerating the development of elite soybean cultivars to meet the growing demands of global soybean production.

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.

Transcription factor ZjABF1 promotes sugar accumulation and abiotic resistance by positively regulating the expression of sugar transport protein ZjSWEET11 and ZjSWEET18 in Chinese jujube

Chinese jujube (Ziziphus jujuba Mill.) exhibits a remarkable resilience to both drought and salinity. Additionally, it is characterized by a high sugar content, with sucrose being the predominant component of its soluble sugars. However, the molecular mechanisms linking robust abiotic stress resistance, sugar accumulation and sugar transport proteins ZjSWEETs remain poorly understood in jujube. In this study, we identified two critical sugar transport proteins, ZjSWEET11 and ZjSWEET18, in Chinese jujube through comprehensive assays and established a positive correlation between sucrose accumulation and the expression of these genes. Furthermore, we discovered that the key transcription factor ZjABF1 within the ABA signaling pathway positively regulated the transcriptional expression of ZjSWEET11 and ZjSWEET18 and increased the sugar contents, consequently improving the drought and salt stress resistance of plants. Basing on these results, we proposed a working module that ZjABF1 promotes sugar accumulation and improves stress resistance by targeting and up-regulating of ZjSWEET11 and ZjSWEET18. Our findings provide valuable insights into the mechanisms underlying sugar accumulation and abiotic stress adaptation in Chinese jujube.

Comparative Proteomic Atlas of Two Soybean Varieties with Contrasting Seed Oil and Protein Content

As complex quantitative traits, soybean seed oil and protein contents are governed by dynamic proteome networks that remain largely unknown. Here, we investigated the dynamic changes of the proteome during seed maturation across two soybean varieties with contrasting seed oil and protein content. Through optimizing the detectability of low-abundance proteins and utilizing library-free data-independent acquisition (directDIA) mass spectrometry, we unprecedentedly identified 7414 proteins and 3975 protein groups (PGs), substantially expanding the soybean seed proteome. Among the PGs, 1391 differentially accumulated between the two varieties. By comparing the abundance of PGs between the two varieties, we identified the core and periphery proteome of soybean seeds and revealed that variations in the oil and protein content are primarily attributed to the peripheral proteome, which significantly fluctuated across seed developmental stages. Our work presents a quantitative proteomic atlas underlying the variation of seed oil and protein content in soybean varieties and provides insight into the mechanisms regulating the seed oil and protein content in soybean.

Characterization of the SWEET Gene Family in Blueberry (Vaccinium corymbosum L.) and the Role of VcSWEET6 Related to Sugar Accumulation in Fruit Development

Sugars will eventually be exported transporters (SWEETs) are essential transmembrane proteins involved in plant growth, stress responses, and plant-pathogen interactions. Despite their importance, systematic studies on SWEETs in blueberries (Vaccinium corymbosum L.) are limited. Blueberries are recognized for their rapid growth and the significant impact of sugar content on fruit flavor, yet the role of the SWEET gene family in sugar accumulation during fruit development remains unclear. In this study, 23 SWEET genes were identified in blueberry, and their phylogenetic relationships, duplication events, gene structures, cis-regulatory elements, and expression profiles were systematically analyzed. The VcSWEET gene family was classified into four clades. Structural and motif analysis revealed conserved exon-intron organization within each clade. RT-qPCR analysis showed widespread expression of VcSWEETs across various tissues and developmental stages, correlating with promoter cis-elements. VcSWEET6a, in particular, was specifically expressed in fruit and showed reduced expression during fruit maturation. Subcellular localization indicated that VcSWEET6a is located in the endoplasmic reticulum. Functional assays in yeast confirmed its role in glucose and fructose uptake, with transport activity inhibited at higher sugar concentrations. Overexpression of VcSWEET6a in blueberries resulted in reduced sugar accumulation. These findings offer valuable insights into the role of VcSWEETs in blueberry sugar metabolism.

A spatially resolved multi-omic single-cell atlas of soybean development

Cis-regulatory elements (CREs) precisely control spatiotemporal gene expression in cells. Using a spatially resolved single-cell atlas of gene expression with chromatin accessibility across ten soybean tissues, we identified 103 distinct cell types and 303,199 accessible chromatin regions (ACRs). Nearly 40% of the ACRs showed cell-type-specific patterns and were enriched for transcription factor (TF) motifs defining diverse cell identities. We identified de novo enriched TF motifs and explored the conservation of gene regulatory networks underpinning legume symbiotic nitrogen fixation. With comprehensive developmental trajectories for endosperm and embryo, we uncovered the functional transition of the three sub-cell types of endosperm, identified 13 sucrose transporters sharing the DNA binding with one finger 11 (DOF11) motif that were co-upregulated in late peripheral endosperm, and identified key embryo cell-type specification regulators during embryogenesis, including a homeobox TF that promotes cotyledon parenchyma identity. This resource provides a valuable foundation for analyzing gene regulatory programs in soybean cell types across tissues and life stages.

Enhanced Analysis of Low-Abundance Proteins in Soybean Seeds Using Advanced Mass Spectrometry

This study presents an advanced approach for the comprehensive analysis of low-abundance proteins in soybean seeds, addressing challenges posed by high-abundance storage proteins. We compared the effectiveness of Data-Dependent Acquisition (DDA), Data-Independent Acquisition (DIA), and BoxCar mass spectrometry techniques to identify low-abundance proteins in two types of soybean seeds: High-Oil and High-Protein seeds. The results indicate that the DIA method, and particularly the BoxCar methods, significantly improve the detection of low-abundance proteins compared to DDA, offering deeper insights into soybean seed biology. Specifically, BoxCar-based analysis revealed distinct proteomic differences between High-Oil and High-Protein seeds, highlighting more active metabolic processes in High-Oil seeds. Additionally, several key proteins were identified and annotated as uniquely expressed in either High-Oil or High-Protein seeds. These findings emphasize the importance of advanced proteomic techniques, such as BoxCar, in deepening our understanding of soybean seed biology and supporting breeding strategies to improve nutritional qualities.

Multi-omics analysis identified the GmUGT88A1 gene, which coordinately regulates soybean resistance to cyst nematode and isoflavone content

Soybean cyst nematode (SCN, Heterodera glycines) is a major pathogen harmful to soybean all over the world, causing huge yield loss every year. Soybean resistance to SCN is a complex quantitative trait controlled by a small number of major genes (rhg1 and Rhg4) and multiple micro-effect genes. Therefore, the continuous identification of new resistant lines and genes is needed for the sustainable development of global soybean production. Here, a novel disease-resistance quantitative trait locus Rscn-16 was identified and fine mapped to an 8.4-kb interval on chromosome 16 using an F2 population. According to transcriptome and metabolome analysis, a UDP-glucosyltransferase encoding gene, GmUGT88A1, was identified as the most likely gene of Rscn-16. Soybean lines overexpressing GmUGT88A1 exhibited increased resistance to SCN, higher isoflavone glycosides and larger seed size while the phenotype of RNA-interference and knockout soybean lines showed sensitivity to SCN and decreased in seed size compared to wild-type plants. GmMYB29 gene could bind to the promoter of GmUGT88A1 and coordinate with GmUGT88A1 to regulate soybean resistance to SCN and isoflavone accumulation. Under SCN infection, GmUGT88A1 participated in the reorientation of isoflavone biosynthetic metabolic flow and the accumulation of isoflavone glycosides, thus protecting soybean from SCN stress. GmUGT88A1 was found to control soybean seed size by affecting transcription abundance of GmSWEET10b and GmFAD3C, which are known to control soybean seed weight. Our findings provide insights into the regulation of SCN resistance, isoflavone content and seed size through metabolic flux redirection, and offer a potential means for soybean improvement.

The comprehensive regulatory network in seed oil biosynthesis

Plant oils play a crucial role in human nutrition, industrial applications and biofuel production. While the enzymes involved in fatty acid (FA) biosynthesis are well-studied, the regulatory networks governing these processes remain largely unexplored. This review explores the intricate regulatory networks modulating seed oil biosynthesis, focusing on key pathways and factors. Seed oil content is determined by the efficiency of de novo FA synthesis as well as influenced by sugar transport, lipid metabolism, FA synthesis inhibitors and fine-tuning mechanisms. At the center of this regulatory network is WRINKLED1 (WRI1), which plays a conserved role in promoting seed oil content across various plant species. WRI1 interacts with multiple proteins, and its expression level is regulated by upstream regulators, including members of the LAFL network. Beyond the LAFL network, we also discuss a potential nuclear factor-Y (NF-Y) regulatory network in soybean with an emphasis on NF-YA and NF-YB and their associated proteins. This NF-Y network represents a promising avenue for future efforts aimed at enhancing oil accumulation and improving stress tolerance in soybean. Additionally, the application of omics-based approaches is of great significance. Advances in omics technologies have greatly facilitated the identification of gene resources, opening new opportunities for genetic improvement. Importantly, several transcription factors involved in oil biosynthesis also participate in stress responses, highlighting a potential link between the two processes. This comprehensive review elucidates the complex mechanisms underlying the regulation of oil biosynthesis, offering insights into potential biotechnological strategies for improving oil production and stress tolerance in oil crops.

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.

FvPHR1 Improves the Quality of Woodland Strawberry Fruit by Up-Regulating the Expression of FvPHT1;7 and FvSWEET9

Phosphorus (P) is vital for plant growth, and continuous P fertiliser application is necessary to increase yield and quality, but it can cause environmental pollution. Plants maintain a steady phosphate (Pi) supply through complex signalling pathways. Phosphate starvation response 1 (PHR1), a key regulator of Pi starvation signals in plants, enables plants to maintain a sufficient Pi level. However, the role of PHR1 in fruit quality remains unclear. In this study, we determined the function of PHR1 in Fragaria vesca (FvPHR1) by overexpressing the FvPHR1 gene. We identified and validated two downstream genes of FvPHR1 by investigating plant phenotypes and analysing RNA-Seq data. FvPHR1 directly enhanced the expression of phosphate transporter 1;7 (FvPHT1;7), increasing Pi uptake and improving photosynthesis efficiency. Additionally, FvPHR1 upregulated the expression of sugar will eventually be exported transporter 9 (FvSWEET9), which encodes a sugar transporter that facilitates sugar transport from leaves to fruit. FvPHR1 can enhance photosynthetic products in a source via the phosphate signalling pathway and facilitate sugar transport to a sink through FvSWEET9. FvPHR1 plays a complicated role in improving fruit quality, providing a molecular foundation for developing strawberry cultivars with highly efficient Pi utilisation processes and high sugar content.

Dissecting the cellular architecture and genetic circuitry of the soybean seed

Seeds are complex structures composed of three regions, embryo, endosperm, and seed coat, with each further divided into subregions that consist of tissues, cell layers, and cell types. Although the seed is well characterized anatomically, much less is known about the genetic circuitry that dictates its spatial complexity. To address this issue, we profiled mRNAs from anatomically distinct seed subregions at several developmental stages. Analyses of these profiles showed that all subregions express similar diverse gene numbers and that the small gene numbers expressed subregion specifically provide information about the biological processes that occur in these seed compartments. In parallel, we profiled RNAs in individual nuclei and identified nuclei clusters representing distinct cell identities. Integrating single-nucleus RNA and subregion mRNA transcriptomes allowed most cell identities to be assigned to specific subregions and cell types and/or cell states. The number of cell identities exceeds the number of anatomically distinguishable cell types, emphasizing the spatial complexity of seeds. We defined gene coexpression networks that underlie distinct biological processes during seed development. We showed that network distribution among subregions and cell identities is highly variable. Some networks operate in single subregions and/or cell identities, and many coexpression networks operate in multiple subregions and/or cell identities. We also showed that single subregions and cell identities possess several networks. Together, our studies provide unique insights into the biological processes and genetic circuitry that underlie the spatial landscape of the seed.

Multiomic analysis of genes related to oil traits in legumes provide insights into lipid metabolism and oil richness in soybean

Soybean (Glycine max) and common bean (Phaseolus vulgaris) diverged approximately 19 million years ago. While these species share a whole-genome duplication (WGD), the Glycine lineage experienced a second, independent WGD. Despite the significance of these WGDs, their impact on gene families related to oil-traits remains poorly understood. Here, we report an in-depth investigation of oil-related gene families in soybean, common bean, and twenty-eight other legume species. We adopted a systematic approach that included 605 RNAseq samples for transcriptome and co-expression analyses, identification of orthologous groups, gene duplication modes and evolutionary rates, and family expansions and contractions. We curated a list of oil candidate genes and found that 91.5% of the families containing these genes expanded in soybean in comparison to common bean. Notably, we observed an expansion of triacylglycerol (TAG) biosynthesis (∼3:1) and an erosion of TAG degradation (∼1.4:1) families in soybean in comparison to common bean. In addition, TAG degradation genes were two-fold more expressed in common bean than in soybean, suggesting that oil degradation is also important for the sharply contrasting seed oil contents in these species. We found 17 transcription factor hub genes that are likely regulators of lipid metabolism. Finally, we inferred expanded and contracted families and correlated these patterns with oil content found in different legume species. In summary, our results do not only shed light on the evolution of oil metabolism genes in soybean, but also present multifactorial evidence supporting the prioritization of promising candidate genes that, if experimentally validated, could accelerate the development of high-oil soybean varieties.

Development and Validation of Kompetitive Allele-Specific Polymerase Chain Reaction Markers for Seed Protein Content in Soybean

Seed protein content is a critical trait in soybean breeding, as it provides a primary source of high-quality protein for both human consumption and animal feed. This study aimed to enhance molecular marker-assisted selection for high-protein soybean varieties by developing Kompetitive Allele-Specific Polymerase Chain Reaction (KASP) markers targeted at loci associated with seed protein content. Nineteen markers with high genotyping efficacy were identified through screening. Utilizing SN76 (a high-protein line) as the male parent and SN49 and DS1 (both low-protein lines) as female parents, 484 F6 generation individuals from these hybrid combinations were selected to validate the predictive accuracy of the 19 KASP markers. Notably, KASP-Pro-1, KASP-Pro-2, and KASP-Pro-3 effectively distinguished genotypes associated with high and low protein content, with prediction accuracies of 68.4%, 75.0%, and 83.3%, respectively. These results underscore the reliability and practical utility of the selected molecular markers, which are located within the genes Glyma.03G219900, Glyma.14G119000, and Glyma.17G074400, respectively. Haplotype analysis and gene pyramiding indicate that these three genes may influence seed protein content. Consequently, these KASP markers can be effectively integrated into genetic and genomic research on soybean seed protein content as well as into marker-assisted breeding.

GmbZIP4a/b Positively Regulate Nodule Number by Affecting Cytokinin Biosynthesis in Glycine max

Legumes have the capability to form nodules that facilitate symbiotic nitrogen fixation (SNF) with rhizobia. Given the substantial energy consumption during the process of SNF, legumes need to optimize nodule number in response to everchanging environmental scenarios. The TGACG BINDING FACTOR1/4 (TGA1/4) are key players in the basal immune response of plants. In this study, both β-glucuronidase staining and quantitative reverse transcription PCR (qRT-PCR) demonstrated that both GmbZIP4a and GmbZIP4b are inducible upon rhizobial inoculation. To investigate their roles further, we constructed gmbzip4a/b double mutants using CRISPR/Cas9 system. Nodulation assessments revealed that these double mutants displayed a reduction in the number of infection threads, which subsequently resulted in a decreased nodule number. However, the processes associated with nodule development including nodule fresh weight, structural characteristics, and nitrogenase activity, remained unaffected in the double mutants. Subsequent transcriptome analyses revealed that zeatin biosynthesis was downregulated in gmbzip4a/b mutants post rhizobial inoculation. Supporting these findings, genes associated with cytokinin (CTK) signaling pathway were upregulated in Williams 82 (Wm82), but this upregulation was not observed in the double mutants after rhizobial treatment. These results suggest that GmbZIP4a/b positively influences nodule formation by promoting the activation of the CTK signaling pathway during the early stages of nodule formation.

Soybean genomics research community strategic plan: A vision for 2024-2028

This strategic plan summarizes the major accomplishments achieved in the last quinquennial by the soybean [Glycine max (L.) Merr.] genetics and genomics research community and outlines key priorities for the next 5 years (2024-2028). This work is the result of deliberations among over 50 soybean researchers during a 2-day workshop in St Louis, MO, USA, at the end of 2022. The plan is divided into seven traditional areas/disciplines: Breeding, Biotic Interactions, Physiology and Abiotic Stress, Functional Genomics, Biotechnology, Genomic Resources and Datasets, and Computational Resources. One additional section was added, Training the Next Generation of Soybean Researchers, when it was identified as a pressing issue during the workshop. This installment of the soybean genomics strategic plan provides a snapshot of recent progress while looking at future goals that will improve resources and enable innovation among the community of basic and applied soybean researchers. We hope that this work will inform our community and increase support for soybean research.

GhDOFD45 promotes sucrose accumulation in cotton seeds by transcriptionally activating GhSWEET10 expression

Cotton seed development and fiber elongation are the inseparable and overlapped development processes requiring the continuous supply of sucrose as the direct carbon source. However, little is known about the molecular mechanism of how sucrose is transported from the source tissues (leaves) into growing cotton seeds. Here, we identify the function of a sucrose transporter gene, Sugars Will Eventually be Exported Transporter 10, GhSWEET10 in cotton seed development. GhSWEET10 encodes a functional sucrose transporter, predominantly expressing in the funiculus, inner seedcoat, and endosperm during fiber elongation. GhSWEET10 RNAi plants (GhSWEET10i) accumulated less sucrose and glucose in growing seeds and that led to shorter fibers and smaller seeds, whereas GhSWEET10 overexpressed plants (GhSWEET10OE) had bigger seeds and longer fibers with more sugar accumulation during fiber elongation. GhSWEET10 gene is transcriptionally controlled by the transcription factor GhDOFD45. GhDOFD45 knockout plants (GhDOFD45-KO) possessed the phenotypes of smaller seeds and shorter fibers like those of GhSWEET10i plants. Furthermore, GhSWEET10 mainly exports the sucrose from the funiculus into developing seeds according to the mimic-analysis of sucrose transporting. Collectively, all these findings show that GhDOFD45 positively regulates GhSWEET10 expression to mainly transport sucrose from leaves into developing cotton seeds. Our findings also imply that the sucrose transport into enlarging seeds benefits fiber development, and thus GhSWEET10 can be selected as a target of breeding novel cotton varieties with larger and more vigorous seeds.

Exploring SoySNP50K and USDA Germplasm Collection Data to Find New QTLs Associated with Protein and Oil Content in Brazilian Genotypes

Genetic diversity within a germplasm collection plays a vital role in the success of breeding programs. However, comprehending this diversity and identifying accessions with desirable traits pose significant challenges. This study utilized publicly available data to investigate SNP markers associated with protein and oil content in Brazilian soybeans. Through this research, twenty-two new QTLs (Quantitative Trait Loci) were identified, and we highlighted the substantial influence of Roanoke, Lee and Bragg ancestor on the genetic makeup of Brazilian soybean varieties. Our findings demonstrate that certain markers are being lost in modern cultivars, while others maintain or even increase their frequency. These observations indicate genomic regions that have undergone selection during soybean introduction in Brazil and could be valuable in breeding programs aimed at enhancing protein or oil content.

Soybean Oil and Protein: Biosynthesis, Regulation and Strategies for Genetic Improvement

Soybean (Glycine max [L.] Merr.) is one of the world's most important sources of oil and vegetable protein. Much of the energy required for germination and early growth of soybean seeds is stored in fatty acids, mainly as triacylglycerols (TAGs), and the main seed storage proteins are β-conglycinin (7S) and glycinin (11S). Recent research advances have deepened our understanding of the biosynthetic pathways and transcriptional regulatory networks that control fatty acid and protein synthesis in organelles such as the plastid, ribosome and endoplasmic reticulum. Here, we review the composition and biosynthetic pathways of soybean oils and proteins, summarizing the key enzymes and transcription factors that have recently been shown to regulate oil and protein synthesis/metabolism. We then discuss the newest genomic strategies for manipulating these genes to increase the food value of soybeans, highlighting important priorities for future research and genetic improvement of this staple crop.

Identification of QTL-allele systems of seed size and oil content for simultaneous genomic improvement in Northeast China soybeans

Northeast China (NEC) is the major production area for soybeans in China, whereas its soybean germplasm has played key roles in world soybean production, especially in the Americas. For plant breeding, genomic selection involves two stages, cross design and progeny selection, with the former determining the latter's potential. In NEC, one of the major breeding purposes is for 100-seed weight (100SW) and seed oil content (SOC). A diverse sample with 361 NEC soybean germplasm accessions was evaluated for their 100SW and SOC in Tieling, Liaoning, China. Both traits exhibited significant phenotypic, genotypic, and G × E variation, with a trait heritability of 82.38% and 86.26%, respectively. A restricted two-stage multi-locus genome-wide association study (RTM-GWAS) with 15,501 SNPLDB (SNP linkage disequilibrium block) markers identified 80 and 92 QTLs with 230 and 299 alleles for 100SW and SOC, respectively. Corresponding to some increase of the two traits, almost all the alleles in the early maturity groups (MG 0 + 00 + 000) were inherited from the late MGs (MG I+II+III), indicating that genetic recombination was the major motivator in addition to a few allele emergence and some allele exclusion fluctuations among early MGs. Using the 95th percentile as indicator, the prediction of recombination potentials showed that 30.43 g 100SW and 27.73% SOC might be achieved, respectively. Three strategies of simultaneous genomic improvement of both traits in designing optimal crosses, namely, 100SW-first, SOC-first, and 100SW-SOC-balance, were proved to be efficient. Thus, the optimal cross design could be extended to multiple traits based on a relatively thorough identification of the QTL-alleles using RTM-GWAS.

Cytosolic fructose-1,6-bisphosphatase isoform mediates metabolic adjustments in bean fruit pericarp to support seed growth

Seed development requires substantial metabolic resources and is influenced by adverse environmental conditions. However, the ability of plants to produce viable seeds under restrictive conditions suggests the existence of mechanisms that make this process less sensitive to environmental stress. Uncovering their regulation could lead to the development of genotypes better adapted to stressful conditions. Plant response to stress is complex, and the contribution of organs such as the fruit pericarp to stress tolerance mechanism may have been underestimated. The bean fruit pericarp, a photosynthetic structure that contributes to seed development, can synthesize starch from surplus sucrose, which is later degraded during the rapid seed growth phase. This metabolic flexibility may be crucial for supporting seed growth when the photosynthate supply is reduced. To explore this possibility, we disrupted phloem continuity at the pedicel level in fruits about to enter the seed reserve accumulation stage. We used the capacity of the pericarp to incorporate 14CO2 to investigate changes in its metabolism. Our findings reveal that, in response to reduced photosynthate availability, the fruit pericarp did not increase 14CO2 fixation. However, the amount of 14C used for starch synthesis decreased, while the proportion used for soluble sugars synthesis increased. This shift resulted in an increase in 14C-products transported to seeds was accompanied by a significant increase in the activity of cytosolic fructose 1,6-bisphosphatase. Our results indicate that photosynthate restriction accelerates the degradation of pericarp storage proteins, and the increase in cFBPase activity could be crucial in converting the carbon produced in carbohydrates.

Metabolic pathways regulated by strigolactones foliar spraying enhance osmoregulation and antioxidant defense in drought-prone soybean

BACKGROUND

Drought stress is a significant abiotic stressor that hinders growth, development, and crop yield in soybeans. Strigolactones (SLs) positively regulate plant resistance to drought stress. However, the impact of foliar application of SLs having different concentrations on soybean growth and metabolic pathways related to osmoregulation remains unknown. Therefore, to clarify the impact of SLs on soybean root growth and cellular osmoregulation under drought stress, we initially identified optimal concentrations and assessed key leaf and root indices. Furthermore, we conducted transcriptomic and metabolic analyses to identify differential metabolites and up-regulated genes.

RESULTS

The results demonstrated that drought stress had a significant impact on soybean biomass, root length, root surface area, water content and photosynthetic parameters. However, when SLs were applied through foliar application at appropriate concentrations, the accumulation of ABA and soluble protein increased, which enhanced drought tolerance of soybean seedlings by regulating osmotic balance, protecting membrane integrity, photosynthesis and activating ROS scavenging system. This also led to an increase in soybean root length, lateral root number and root surface area. Furthermore, the effects of different concentrations of SLs on soybean leaves and roots were found to be time-sensitive. However, the application of 0.5 µM SLs had the greatest beneficial impact on soybean growth and root morphogenesis under drought stress. A total of 368 differential metabolites were screened in drought and drought plus SLs treatments. The up-regulated genes were mainly involved in nitrogen compound utilization, and the down-regulated metabolic pathways were mainly involved in maintaining cellular osmoregulation and antioxidant defenses.

CONCLUSIONS

SLs enhance osmoregulation in soybean plants under drought stress by regulating key metabolic pathways including Arachidonic acid metabolism, Glycerophospholipid metabolism, Linoleic acid metabolism, and Flavone and flavonol biosynthesis. This study contributes to the theoretical understanding of improving soybean adaptability and survival in response to drought stress.

Utility of Arabidopsis KASII Promoter in Development of an Effective CRISPR/Cas9 System for Soybean Genome Editing and Its Application in Engineering of Soybean Seeds Producing Super-High Oleic and Low Saturated Oils

This study reports the use of the Arabidopsis KASII promoter (AtKASII) to develop an efficient CRISPR/Cas9 system for soybean genome editing. When this promoter was paired with Arabidopsis U6 promoters to drive Cas9 and single guide RNA expression, respectively, simultaneous editing of the three fatty acid desaturase genes GmFAD2-1A, GmFAD2-1B, and GmFAD3A occurred in more than 60% of transgenic soybean lines at T2 generation, and all the triple mutants possessed desirable high-oleic traits. In sharp contrast, not a single line underwent simultaneous editing of the three target genes when AtKASII was replaced by the widely used AtEC1.2 promoter. Furthermore, our study showed that the stable and inheritable mutations in the high-oleic lines did not alter the overall contents of oil and protein or amino acid composition while increasing the oleic acid content up to 87.6% from approximately 23.8% for wild-type seeds, concomitant with 34.4- and 3.7-fold reductions in linoleic and linolenic acid, respectively. Collectively, this study demonstrates that the AtKASII promoter is highly promising for optimization of the CRISPR/Cas9 system for genome editing in soybean and possibly beyond.

Soybean hypocotyl elongation is regulated by a MYB33-SWEET11/21-GA2ox8c module involving long-distance sucrose transport

The length of hypocotyl affects the height of soybean and lodging resistance, thus determining the final grain yield. However, research on soybean hypocotyl length is scarce, and the regulatory mechanisms are not fully understood. Here, we identified a module controlling the transport of sucrose, where sucrose acts as a messenger moved from cotyledon to hypocotyl, regulating hypocotyl elongation. This module comprises four key genes, namely MYB33, SWEET11, SWEET21 and GA2ox8c in soybean. In cotyledon, MYB33 is responsive to sucrose and promotes the expression of SWEET11 and SWEET21, thereby facilitating sucrose transport from the cotyledon to the hypocotyl. Subsequently, sucrose transported from the cotyledon up-regulates the expression of GA2ox8c in the hypocotyl, which ultimately affects the length of the hypocotyl. During the domestication and improvement of soybean, an allele of MYB33 with enhanced abilities to promote SWEET11 and SWEET21 has gradually become enriched in landraces and cultivated varieties, SWEET11 and SWEET21 exhibit high conservation and have undergone a strong purified selection and GA2ox8c is under a strong artificial selection. Our findings identify a new molecular pathway in controlling soybean hypocotyl elongation and provide new insights into the molecular mechanism of sugar transport in soybean.

Natural variation in GmSW17 controls seed size in soybean

Seed size/weight plays an important role in determining crop yield, yet only few genes controlling seed size have been characterized in soybean. Here, we perform a genome-wide association study and identify a major quantitative trait locus (QTL), named GmSW17 (Seed Width 17), on chromosome 17 that determine soybean seed width/weight in natural population. GmSW17 encodes a ubiquitin-specific protease, an ortholog to UBP22, belonging to the ubiquitin-specific protease (USPs/UBPs) family. Further functional investigations reveal that GmSW17 interacts with GmSGF11 and GmENY2 to form a deubiquitinase (DUB) module, which influences H2Bub levels and negatively regulates the expression of GmDP-E2F-1, thereby inhibiting the G1-to-S transition. Population analysis demonstrates that GmSW17 undergo artificial selection during soybean domestication but has not been fixed in modern breeding. In summary, our study identifies a predominant gene related to soybean seed weight, providing potential advantages for high-yield breeding in soybean.

Identification of a Major QTL for Seed Protein Content in Cultivated Peanut (Arachis hypogaea L.) Using QTL-Seq

Peanut (Arachis hypogaea L.) is a great plant protein source for human diet since it has high protein content in the kernel. Therefore, seed protein content (SPC) is considered a major agronomic and quality trait in peanut breeding. However, few genetic loci underlying SPC have been identified in peanuts, and the underlying regulatory mechanisms remain unknown, limiting the effectiveness of breeding for high-SPC peanut varieties. In this study, a major QTL (qSPCB10.1) controlling peanut SPC was identified within a 2.3 Mb interval in chromosome B10 by QTL-seq using a recombinant inbred line population derived from parental lines with high and low SPCs, respectively. Sequence comparison, transcriptomic analysis, and annotation analysis of the qSPCB10.1 locus were performed. Six differentially expressed genes with sequence variations between two parents were identified as candidate genes underlying qSPCB10.1. Further locus interaction analysis revealed that qSPCB10.1 could not affect the seed oil accumulation unless qOCA08.1XH13 was present, a high seed oil content (SOC) allele for a major QTL underlying SOC. In summary, our study provides a basis for future investigation of the genetic basis of seed protein accumulation and facilitates marker-assisted selection for developing high-SPC peanut genotypes.

Insights into the ameliorative effect of ZnONPs on arsenic toxicity in soybean mediated by hormonal regulation, transporter modulation, and stress responsive genes

Arsenic (As) contamination of agricultural soils poses a serious threat to crop productivity and food safety. Zinc oxide nanoparticles (ZnONPs) have emerged as a potential amendment for mitigating the adverse effects of As stress in plants. Soybean crop is mostly grown on marginalized land and is known for high accumulation of As in roots than others tissue. Therefore, this study aimed to elucidate the underlying mechanisms of ZnONPs in ameliorating arsenic toxicity in soybean. Our results demonstrated that ZnOB significantly improved the growth performance of soybean plants exposed to arsenic. This improvement was accompanied by a decrease (55%) in As accumulation and an increase in photosynthetic efficiency. ZnOB also modulated hormonal balance, with a significant increase in auxin (149%), abscisic acid (118%), gibberellin (160%) and jasmonic acid content (92%) under As(V) stress assuring that ZnONPs may enhance root growth and development by regulating hormonal signaling. We then conducted a transcriptomic analysis to understand further the molecular mechanisms underlying the NPs-induced As(V) tolerance. This analysis identified genes differentially expressed in response to ZnONPs supplementation, including those involved in auxin, abscisic acid, gibberellin, and jasmonic acid biosynthesis and signaling pathways. Weighted gene co-expression network analysis identified 37 potential hub genes encoding stress responders, transporters, and signal transducers across six modules potentially facilitated the efflux of arsenic from cells, reducing its toxicity. Our study provides valuable insights into the molecular mechanisms associated with metalloid tolerance in soybean and offers new avenues for improving As tolerance in contaminated soils.

A spatially resolved multiomic single-cell atlas of soybean development

Cis-regulatory elements (CREs) precisely control spatiotemporal gene expression in cells. Using a spatially resolved single-cell atlas of gene expression with chromatin accessibility across ten soybean tissues, we identified 103 distinct cell types and 303,199 accessible chromatin regions (ACRs). Nearly 40% of the ACRs showed cell-type-specific patterns and were enriched for transcription factor (TF) motifs defining diverse cell identities. We identified de novo enriched TF motifs and explored conservation of gene regulatory networks underpinning legume symbiotic nitrogen fixation. With comprehensive developmental trajectories for endosperm and embryo, we uncovered the functional transition of the three sub-cell types of endosperm, identified 13 sucrose transporters sharing the DOF11 motif that were co-up-regulated in late peripheral endosperm and identified key embryo cell-type specification regulators during embryogenesis, including a homeobox TF that promotes cotyledon parenchyma identity. This resource provides a valuable foundation for analyzing gene regulatory programs in soybean cell types across tissues and life stages.

Exploring selection signatures in the divergence and evolution of lipid droplet (LD) associated genes in major oilseed crops

BACKGROUND

Oil bodies or lipid droplets (LDs) in the cytosol are the subcellular storage compartments of seeds and the sites of lipid metabolism providing energy to the germinating seeds. Major LD-associated proteins are lipoxygenases, phospholipaseD, oleosins, TAG-lipases, steroleosins, caleosins and SEIPINs; involved in facilitating germination and enhancing peroxidation resulting in off-flavours. However, how natural selection is balancing contradictory processes in lipid-rich seeds remains evasive. The present study was aimed at the prediction of selection signatures among orthologous clades in major oilseeds and the correlation of selection effect with gene expression.

RESULTS

The LD-associated genes from the major oil-bearing crops were analyzed to predict natural selection signatures in phylogenetically close-knit ortholog clusters to understand adaptive evolution. Positive selection was the major force driving the evolution and diversification of orthologs in a lineage-specific manner. Significant positive selection effects were found in 94 genes particularly in oleosin and TAG-lipases, purifying with excess of non-synonymous substitution in 44 genes while 35 genes were neutral to selection effects. No significant selection impact was noticed in Brassicaceae as against LOX genes of oil palm. A heavy load of deleterious mutations affecting selection signatures was detected in T-lineage oleosins and LOX genes of Arachis hypogaea. The T-lineage oleosin genes were involved in mainly anther, tapetum and anther wall morphogenesis. In Ricinus communis and Sesamum indicum > 85% of PLD genes were under selection whereas selection pressures were low in Brassica juncea and Helianthus annuus. Steroleosin, caleosin and SEIPINs with large roles in lipid droplet organization expressed mostly in seeds and were under considerable positive selection pressures. Expression divergence was evident among paralogs and homeologs with one gene attaining functional superiority compared to the other. The LOX gene Glyma.13g347500 associated with off-flavor was not expressed during germination, rather its paralog Glyma.13g347600 showed expression in Glycine max. PLD-α genes were expressed on all the tissues except the seed,δ genes in seed and meristem while β and γ genes expressed in the leaf.

CONCLUSIONS

The genes involved in seed germination and lipid metabolism were under strong positive selection, although species differences were discernable. The present study identifies suitable candidate genes enhancing seed oil content and germination wherein directional selection can become more fruitful.

Co-overexpression of SWEET sucrose transporters modulates sucrose synthesis and defence responses to enhance immunity against bacterial blight in rice

Enhancing carbohydrate export from source to sink tissues is considered to be a realistic approach for improving photosynthetic efficiency and crop yield. The rice sucrose transporters OsSUT1, OsSWEET11a and OsSWEET14 contribute to sucrose phloem loading and seed filling. Crucially, Xanthomonas oryzae pv. oryzae (Xoo) infection in rice enhances the expression of OsSWEET11a and OsSWEET14 genes, and causes leaf blight. Here we show that co-overexpression of OsSUT1, OsSWEET11a and OsSWEET14 in rice reduced sucrose synthesis and transport leading to lower growth and yield but reduced susceptibility to Xoo relative to controls. The immunity-related hypersensitive response (HR) was enhanced in the transformed lines as indicated by the increased expression of defence genes, higher salicylic acid content and presence of HR lesions on the leaves. The results suggest that the increased expression of OsSWEET11a and OsSWEET14 in rice is perceived as a pathogen (Xoo) attack that triggers HR and results in constitutive activation of plant defences that are related to the signalling pathways of pathogen starvation. These findings provide a mechanistic basis for the trade-off between plant growth and immunity because decreased susceptibility against Xoo compromised plant growth and yield.

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.

Genome-wide association study and transcriptome analysis reveal natural variation of key genes regulation flowering time in rapeseed

Flowering time is a crucial determinant of both the yield and oil quality of rapeseed as well as a key indicator of plant maturity. We performed a genome-wide association study and transcriptome analysis to identify key genes/loci regulating flowering time in Brassica napus L. Forty-six haplotype regions harboring candidate genes were determined to be significantly associated with flowering time, and 28 of these haplotype regions overlapped with previously reported quantitative trait loci. A further investigation of these haplotype regions revealed nucleotide variations in the genes BnaFT-A02, BnaFRI-A10, and BnaFPA-A09 that correlated with phenotypic variations in flowering time. Furthermore, the co-expression network analysis indicated that BnaFT-A02 is directly linked to BnaFRI-A10 and BnaFPA-A09, in a subnetwork and also associated with 13 vernalization, 31 photoperiod, 33 autonomous pathway, and 10 gibberellin pathway genes, forming a potential network regulating flowering time in rapeseed. These results provide valuable haplotype markers for the breeding of early maturing rapeseed varieties.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01479-4.

Integrative omics analysis elucidates the genetic basis underlying seed weight and oil content in soybean

Synergistic optimization of key agronomic traits by traditional breeding has dramatically enhanced crop productivity in the past decades. However, the genetic basis underlying coordinated regulation of yield- and quality-related traits remains poorly understood. Here, we dissected the genetic architectures of seed weight and oil content by combining genome-wide association studies (GWAS) and transcriptome-wide association studies (TWAS) using 421 soybean (Glycine max) accessions. We identified 26 and 33 genetic loci significantly associated with seed weight and oil content by GWAS, respectively, and detected 5,276 expression quantitative trait loci (eQTLs) regulating expression of 3,347 genes based on population transcriptomes. Interestingly, a gene module (IC79), regulated by two eQTL hotspots, exhibited significant correlation with both seed weigh and oil content. Twenty-two candidate causal genes for seed traits were further prioritized by TWAS, including Regulator of Weight and Oil of Seed 1 (GmRWOS1), which encodes a sodium pump protein. GmRWOS1 was verified to pleiotropically regulate seed weight and oil content by gene knockout and overexpression. Notably, allelic variations of GmRWOS1 were strongly selected during domestication of soybean. This study uncovers the genetic basis and network underlying regulation of seed weight and oil content in soybean and provides a valuable resource for improving soybean yield and quality by molecular breeding.

Differential selection of yield and quality traits has shaped genomic signatures of cowpea domestication and improvement

Cowpeas (tropical legumes) are important in ensuring food and nutritional security in developing countries, especially in sub-Saharan Africa. Herein, we report two high-quality genome assemblies of grain and vegetable cowpeas and we re-sequenced 344 accessions to characterize the genomic variations landscape. We identified 39 loci for ten important agronomic traits and more than 541 potential loci that underwent selection during cowpea domestication and improvement. In particular, the synchronous selections of the pod-shattering loci and their neighboring stress-relevant loci probably led to the enhancement of pod-shattering resistance and the compromise of stress resistance during the domestication from grain to vegetable cowpeas. Moreover, differential selections on multiple loci associated with pod length, grain number per pod, seed weight, pod and seed soluble sugars, and seed crude proteins shaped the yield and quality diversity in cowpeas. Our findings provide genomic insights into cowpea domestication and improvement footprints, enabling further genome-informed cultivar improvement of cowpeas.

Genome-wide transcriptome analysis reveals key regulatory networks and genes involved in the determination of seed hardness in vegetable soybean

Seed hardness is an important quality trait of vegetable soybean. To determine the factors underlying seed hardness, two landraces with contrasting seed hardness, Niumaohuang (low seed hardness) and Pixiansilicao (high seed hardness), were selected from 216 soybean accessions originating from 26 provinces in China. The contents of the main components in vegetable soybean seeds such as water, soluble sugar, starch, protein and oil were measured, and transcriptome analyses performed during five stages of seed developmental. Transcriptome analysis indicates that during the middle and late stages of seed development, a large number of genes involved in the synthesis or degradation of starch, storage protein, and fatty acids were differentially expressed, leading to differences in the accumulation of stored substances during seed maturation among Niumaohuang and Pixiansilicao. The activity of cell proliferation and the formation of cell walls in the middle and late stages of seed development may also affect the hardness of seeds to a certain extent. In addition, weighted gene co-expression network analysis (WGCNA) was undertaken to identify co-expressed gene modules and hub genes that regulate seed hardness. Overexpression of a candidate seed hardness regulatory hub gene, GmSWEET2, resulted in increased seed hardness. In this study, the important role of GmSWEET2 in regulating the hardness of vegetable soybean seeds was verified and numerous potential key regulators controlling seed hardness and the proportion of seed components were identified, laying the groundwork for improving the texture of vegetable soybean.

Gamma-rays induced genome wide stable mutations in cowpea deciphered through whole genome sequencing

PURPOSE

Gamma rays are the most widely exploited physical mutagen in plant mutation breeding. They are known to be involved in the development of more than 60% of global cowpea (Vigna unguiculata (L.) Walp.) mutant varieties. Nevertheless, the nature and type of genome-wide mutations induced by gamma rays have not been studied in cowpea and therefore, the present investigation was undertaken.

MATERIALS AND METHODS

Genomic DNAs from three stable gamma rays-induced mutants (large seed size, small seed size and disease resistant mutant) of cowpea cultivar 'CPD103' in M6 generation along with its progenitor were used for Illumina-based whole-genome resequencing.

RESULTS

Gamma rays induced a relatively higher frequency (88.9%) of single base substitutions (SBSs) with an average transition to transversion ratio (Ti/Tv) of 3.51 in M6 generation. A > G transitions, including its complementary T > C transitions, predominated the transition mutations, while all four types of transversion mutations were detected with frequencies over 6.5%. Indels (small insertions and deletions) constituted about 11% of the total induced variations, wherein small insertions (6.3%) were relatively more prominent than small deletions (4.8%). Among the indels, single-base indels and, in particular, those involving A/T bases showed a preponderance, albeit indels of up to three bases were detected in low proportions. Distributed across all 11 chromosomes, only a fraction of SBSs (19.45%) and indels (20.2%) potentially altered the encoded amino acids/peptides. The inherent mutation rate induced by gamma rays in cowpea was observed to be in the order of 1.4 × 10-7 per base pair in M6 generation.

CONCLUSION

Gamma-rays with a greater tendency to induce SBSs and, to a lesser extent, indels could be efficiently and effectively exploited in cowpea mutation breeding.

The Idesia polycarpa genome provides insights into its evolution and oil biosynthesis

The deciduous tree Idesia polycarpa can provide premium edible oil with high polyunsaturated fatty acid contents. Here, we generate its high-quality reference genome, which is ∼1.21 Gb, comprising 21 pseudochromosomes and 42,086 protein-coding genes. Phylogenetic and genomic synteny analyses show that it diverged with Populus trichocarpa about 16.28 million years ago. Notably, most fatty acid biosynthesis genes are not only increased in number in its genome but are also highly expressed in the fruits. Moreover, we identify, through genome-wide association analysis and RNA sequencing, the I. polycarpa SUGAR TRANSPORTER 5 (IpSTP5) gene as a positive regulator of high oil accumulation in the fruits. Silencing of IpSTP5 by virus-induced gene silencing causes a significant reduction of oil content in the fruits, suggesting it has the potential to be used as a molecular marker to breed the high-oil-content cultivars. Our results collectively lay the foundation for breeding the elite cultivars of I. polycarpa.

Integrated Transcriptomic and Proteomic Characterization of a Chromosome Segment Substitution Line Reveals the Regulatory Mechanism Controlling the Seed Weight in Soybean

Soybean is the major global source of edible oils and vegetable proteins. Seed size and weight are crucial traits determining the soybean yield. Understanding the molecular regulatory mechanism underlying the seed weight and size is helpful for improving soybean genetic breeding. The molecular regulatory pathways controlling the seed weight and size were investigated in this study. The 100-seed weight, seed length, seed width, and seed weight per plant of a chromosome segment substitution line (CSSL) R217 increased compared with those of its recurrent parent 'Suinong14' (SN14). Transcriptomic and proteomic analyses of R217 and SN14 were performed at the seed developmental stages S15 and S20. In total, 2643 differentially expressed genes (DEGs) and 208 differentially accumulated proteins (DAPs) were detected at S15, and 1943 DEGs and 1248 DAPs were detected at S20. Furthermore, integrated transcriptomic and proteomic analyses revealed that mitogen-activated protein kinase signaling and cell wall biosynthesis and modification were potential pathways associated with seed weight and size control. Finally, 59 candidate genes that might control seed weight and size were identified. Among them, 25 genes were located on the substituted segments of R217. Two critical pathways controlling seed weight were uncovered in our work. These findings provided new insights into the seed weight-related regulatory network in soybean.

Significance of Raffinose Family Oligosaccharides (RFOs) metabolism in plants

Raffinose Family Oligosaccharides (RFOs) are a kind of polysaccharide containing D-galactose, and they widely exist in higher plants. Synthesis of RFOs begins with galactinol synthase (GolS; EC 2.4.1.123) to convert myo-inositol into galactinol. The subsequent formation of raffinose and stachyose are catalyzed by raffinose synthase (RS; EC 2.4.1.82) and stachyose synthase (STS; EC 2.4.1.67) using sucrose and galactinol as substrate, respectively. The hydrolysis of RFOs is finished by α-galactosidase (α-Gal; EC 3.2.1.22) to produce sucrose and galactose. Importance of RFOs metabolism have been summarized, e.g. In RFOs translocating plants, the phloem loading and unloading of RFOs are widely reported in mediating the plant development process. Interference function of RFOs synthesis or hydrolysis enzymes caused growth defect. In addition, the metabolism of RFOs involved in the biotic or abiotic stresses was discussed in this review. Overall, this literature summarizes our current understanding of RFOs metabolism and points out knowledge gaps that need to be filled in future.

Genome-wide identification of SWEET genes reveals their roles during seed development in peanuts

Sugar Will Eventually be Exported Transporter (SWEET) proteins are highly conserved in various organisms and play crucial roles in sugar transport processes. However, SWEET proteins in peanuts, an essential leguminous crop worldwide, remain lacking in systematic characterization. Here, we identified 94 SWEET genes encoding the conservative MtN3/saliva domains in three peanut species, including 47 in Arachis hypogea, 23 in Arachis duranensis, and 24 in Arachis ipaensis. We observed significant variations in the exon-intron structure of these genes, while the motifs and domain structures remained highly conserved. Phylogenetic analysis enabled us to categorize the predicted 286 SWEET proteins from eleven species into seven distinct groups. Whole genome duplication/segment duplication and tandem duplication were the primary mechanisms contributing to the expansion of the total number of SWEET genes. In addition, an investigation of cis-elements in the potential promoter regions and expression profiles across 22 samples uncovered the diverse expression patterns of AhSWEET genes in peanuts. AhSWEET24, with the highest expression level in seeds from A. hypogaea Tifrunner, was observed to be localized on both the plasma membrane and endoplasmic reticulum membrane. Moreover, qRT-PCR results suggested that twelve seed-expressed AhSWEET genes were important in the regulation of seed development across four different peanut varieties. Together, our results provide a foundational basis for future investigations into the functions of SWEET genes in peanuts, especially in the process of seed development.

Dt1 inhibits SWEET-mediated sucrose transport to regulate photoperiod-dependent seed weight in soybean

Soybean is a photoperiod-sensitive short-day crop whose reproductive period and yield are markedly affected by day-length changes. Seed weight is one of the key traits determining the soybean yield; however, the prominent genes that control the final seed weight of soybean and the mechanisms underlying the photoperiod's effect on this trait remain poorly understood. In this study, we identify SW19 as a major locus controlling soybean seed weight by QTL mapping and determine Dt1, an orthologous gene of Arabidopsis TFL1 that is known to govern the soybean growth habit, as the causal gene of the SW19 locus. We showed that Dt1 is highly expressed in developing seeds and regulates photoperiod-dependent seed weight in soybean. Further analyses revealed that the Dt1 protein physically interacts with the sucrose transporter GmSWEET10a to negatively regulate the import of sucrose from seed coat to the embryo, thus modulating seed weight under long days. However, Dt1 does not function in seed development under short days due to its very low expression. Importantly, we discovered a novel natural allelic variant of Dt1 (H4 haplotype) that decouples its pleiotropic effects on seed size and growth habit; i.e., this variant remains functional in seed development but fails to regulate the stem growth habit of soybean. Collectively, our findings provide new insights into how soybean seed development responds to photoperiod at different latitudes, offering an ideal genetic component for improving soybean's yield by manipulating its seed weight and growth habit.

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.

Unveiling the molecular regulatory mechanisms underlying sucrose accumulation and oil reduction in peanut kernels through genetic mapping and transcriptome analysis

Sucrose content is a key factor for the flavor of edible peanut, which determines the sweet taste of fresh peanut and also attribute to pleasant flavor of roasted peanut. To explore the genetic mechanism of the sucrose content in peanut, an F2 population was created by crossing the sweet cultivar Zhonghuatian 1 (ZHT1) with Nanyangbaipi (NYBP). A genomic region spanning 28.26 kb on chromosome A06 was identified for the sucrose content through genetic mapping, elucidating 47.5% phenotypic variance explained. As the sucrose content had a significantly negative correlation with the oil content, this region was also found to be related to the oil content explaining 37.2% of phenotype variation. In this region, Arahy.42CAD1 was characterized as the most likely candidate gene through a comprehensive analysis. The nuclear localization of Arahy.42CAD1 suggests its potential involvement in the regulation of gene expression for sucrose and oil contents in peanut. Transcriptome analysis of the developing seeds in both parents revealed that genes involved in glycolysis and triacylglycerol biosynthesis pathways were not significantly down-regulated in ZHT1, indicating that the sucrose accumulation was not attributed to the suppression of triacylglycerol biosynthesis. Based on the WGCNA analysis, Arahy.42CAD1 was co-expressed with the genes involved in vesicle transport and oil body assembly, suggesting that the sucrose accumulation may be caused by disruptions in TAG transportation or storage mechanisms. These findings offer new insights into the molecular mechanisms governing sucrose accumulation in peanut, and also provide a potential gene target for enhancing peanut flavor.

Identification of the domestication gene GmCYP82C4 underlying the major quantitative trait locus for the seed weight in soybean

KEY MESSAGE

A major quantitative trait locus (QTL) for the hundred-seed weight (HSW) was identified and confirmed in the two distinct soybean populations, and the target gene GmCYP82C4 underlying this locus was identified that significantly associated with soybean seed weight, and it was selected during the soybean domestication and improvement process. Soybean is a major oil crop for human beings and the seed weight is a crucial goal of soybean breeding. However, only a limited number of target genes underlying the quantitative trait loci (QTLs) controlling seed weight in soybean are known so far. In the present study, six loci associated with hundred-seed weight (HSW) were detected in the first population of 573 soybean breeding lines by genome-wide association study (GWAS), and 64 gene models were predicted in these candidate QTL regions. The QTL qHSW_1 exhibits continuous association signals on chromosome four and was also validated by region association study (RAS) in the second soybean population (409 accessions) with wild, landrace, and cultivar soybean accessions. There were seven genes in qHSW_1 candidate region by linkage disequilibrium (LD) block analysis, and only Glyma.04G035500 (GmCYP82C4) showed specifically higher expression in flowers, pods, and seeds, indicating its crucial role in the soybean seed development. Significant differences in HSW trait were detected when the association panels are genotyped by single-nucleotide polymorphisms (SNPs) in putative GmCYP82C4 promoter region. Eight haplotypes were generated by six SNPs in GmCYP82C4 in the second soybean population, and two superior haplotypes (Hap2 and Hap4) of GmCYP82C4 were detected with average HSW of 18.27 g and 18.38 g, respectively. The genetic diversity of GmCYP82C4 was analyzed in the second soybean population, and GmCYP82C4 was most likely selected during the soybean domestication and improvement process, leading to the highest proportion of Hap2 of GmCYP82C4 both in landrace and cultivar subpopulations. The QTLs and GmCYP82C4 identified in this study provide novel genetic resources for soybean seed weight trait, and the GmCYP82C4 could be used for soybean molecular breeding to develop desirable seed weight in the future.

Genome-Wide Association Studies Prioritize Genes Controlling Seed Size and Reproductive Period Length in Soybean

Hundred-seed weight (HSW) and reproductive period length (RPL) are two major agronomic traits critical for soybean production and adaptation. However, both traits are quantitatively controlled by multiple genes that have yet to be comprehensively elucidated due to the lack of major genes; thereby, the genetic basis is largely unknown. In the present study, we conducted comprehensive genome-wide association analyses (GWAS) of HSW and RPL with multiple sets of accessions that were phenotyped across different environments. The large-scale analysis led to the identification of sixty-one and seventy-four significant QTLs for HSW and RPL, respectively. An ortholog-based search analysis prioritized the most promising candidate genes for the QTLs, including nine genes (TTG2, BZR1, BRI1, ANT, KLU, EOD1/BB, GPA1, ABA2, and ABI5) for HSW QTLs and nine genes (such as AGL8, AGL9, TOC1, and COL4) and six known soybean flowering time genes (E2, E3, E4, Tof11, Tof12, and FT2b) for RPL QTLs. We also demonstrated that some QTLs were targeted during domestication to drive the artificial selection of both traits towards human-favored traits. Local adaptation likely contributes to the increased genomic diversity of the QTLs underlying RPL. The results provide additional insight into the genetic basis of HSW and RPL and prioritize a valuable resource of candidate genes that merits further investigation to reveal the complex molecular mechanism and facilitate soybean improvement.

Genome-Wide Identification and Expression Profile Analysis of Sugars Will Eventually Be Exported Transporter (SWEET) Genes in Zantedeschia elliottiana and Their Responsiveness to Pectobacterium carotovora subspecies Carotovora (Pcc) Infection

SWEET, sugars will eventually be exported transporter, is a novel class of sugar transporter proteins that can transport sugars across membranes down a concentration gradient. It plays a key role in plant photosynthetic assimilates, phloem loading, nectar secretion from nectar glands, seed grouting, pollen development, pathogen interactions, and adversity regulation, and has received widespread attention in recent years. To date, systematic analysis of the SWEET family in Zantedeschia has not been documented, although the genome has been reported in Zantedeschia elliottiana. In this study, 19 ZeSWEET genes were genome-wide identified in Z. elliottiana, and unevenly located in 10 chromosomes. They were further clustered into four clades by a phylogenetic tree, and almost every clade has its own unique motifs. Synthetic analysis confirmed two pairs of segmental duplication events of ZeSWEET genes. Heatmaps of tissue-specific and Pectobacterium carotovora subsp. Carotovora (Pcc) infection showed that ZeSWEET genes had different expression patterns, so SWEETs may play widely varying roles in development and stress tolerance in Zantedeschia. Moreover, quantitative reverse transcription-PCR (qRT-PCR) analysis revealed that some of the ZeSWEETs responded to Pcc infection, among which eight genes were significantly upregulated and six genes were significantly downregulated, revealing their potential functions in response to Pcc infection. The promoter sequences of ZeSWEETs contained 51 different types of the 1380 cis-regulatory elements, and each ZeSWEET gene contained at least two phytohormone responsive elements and one stress response element. In addition, a subcellular localization study indicated that ZeSWEET07 and ZeSWEET18 were found to be localized to the plasma membrane. These findings provide insights into the characteristics of SWEET genes and contribute to future studies on the functional characteristics of ZeSWEET genes, and then improve Pcc infection tolerance in Zantedeschia through molecular breeding.

Understanding the Molecular Regulatory Networks of Seed Size in Soybean

Soybean being a major cash crop provides half of the vegetable oil and a quarter of the plant proteins to the global population. Seed size traits are the most important agronomic traits determining the soybean yield. These are complex traits governed by polygenes with low heritability as well as are highly influenced by the environment as well as by genotype x environment interactions. Although, extensive efforts have been made to unravel the genetic basis and molecular mechanism of seed size in soybean. But most of these efforts were majorly limited to QTL identification, and only a few genes for seed size were isolated and their molecular mechanism was elucidated. Hence, elucidating the detailed molecular regulatory networks controlling seed size in soybeans has been an important area of research in soybeans from the past decades. This paper describes the current progress of genetic architecture, molecular mechanisms, and regulatory networks for seed sizes of soybeans. Additionally, the main problems and bottlenecks/challenges soybean researchers currently face in seed size research are also discussed. This review summarizes the comprehensive and systematic information to the soybean researchers regarding the molecular understanding of seed size in soybeans and will help future research work on seed size in soybeans.

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

Advanced multi-omics technologies offer much information that can uncover the regulatory mechanisms from genotype to phenotype. In soybean, numerous multi-omics databases have been published. Although they cover multiple omics, there are still limitations when it comes to the types and scales of omics datasets and analysis methods utilized. This study aims to address these limitations by collecting and integrating a comprehensive set of multi-omics datasets. This includes 38 genomes, transcriptomes from 435 tissue samples, 125 phenotypes from 6686 accessions, epigenome data involving histone modification, transcription factor binding, chromosomal accessibility and chromosomal interaction, as well as genetic variation data from 24 501 soybean accessions. Then, common analysis pipelines and statistical methods were applied to mine information from these multi-omics datasets, resulting in the successful establishment of a user-friendly multi-omics database called SoyMD (https://yanglab.hzau.edu.cn/SoyMD/#/). SoyMD provides researchers with efficient query options and analysis tools, allowing them to swiftly access relevant omics information and conduct comprehensive multi-omics data analyses. Another notable feature of SoyMD is its capability to facilitate the analysis of candidate genes, as demonstrated in the case study on seed oil content. This highlights the immense potential of SoyMD in soybean genetic breeding and functional genomics research.