GmKIX8-1 regulates organ size in soybean and is the causative gene for the major seed weight QTL qSw17-1

Loss of phytochromobilin synthase activity leads to larger seeds with higher protein content in soybean

Seed weight is an important agronomic trait that is related to seed size and determines yield in soybean (Glycine max). We previously identified a spontaneous soybean mutant with light green leaves called ygl2. Here, we cloned YGL2, which encodes a phytochromobilin (PΦB) synthase involved in synthesizing the chromophore of the photoreceptor phytochrome. The lesion in ygl2 is a 10-bp deletion, causing a frameshift mutation and a premature stop codon that truncates the encoded protein. In contrast to the wild type, ygl2 lacks PΦB synthase activity and function. This appears to promote cell expansion, thus increasing seed weight. Surprisingly, the ygl2 mutant also exhibits excellent traits including early maturity and high protein content. Moreover, under the condition of dense planting (3 cm), the yield of YGL2 mutant was significantly increased. Mutants harboring ygl2 mutations that we generated via gene editing had enlarged seeds with high protein content. Moreover, the expression levels of the photoperiod sensitive genes (E1, FT2a, FT5a) were lower in the ygl2 mutant than in the wild type. Mutating the YGL2 gene resulted in increased biliverdin content and decreased heme content. We determined that Lhcb4, a chlorophyll a/b binding protein in photosystem II, interacts with YGL2 but not with the mutant version of the protein. We thus identified a mutation in a PΦB synthase gene that enhances seed weight in soybean, providing a promising breeding target for this important crop.

Key Determinants of Seed Size for Enhancing Genetic Gain in Legumes

Legumes play a pivotal role in human nutrition due to their high nutritional value, especially protein content. Therefore, enhancing the productivity of grain legumes is desirable for ensuring food and nutritional security. Seed size and seed weight are key factors influencing productivity. This article consolidates the substantial amount of research conducted to uncover the molecular signatures associated with seed size into a structured format, providing a one-stop platform of available resources for enhancing genetic gains in legumes. The advent of NGS technologies enabled the decryption of genomes and transcriptomes of important grain legumes. Moreover, molecular signatures such as SSRs, SNPs, transcription factors, methylation patterns and so forth scanned from phenotypically and genotypically well-characterized natural and mapping populations helped identify the QTLs, MTAs and candidate genes associated with seed size. Many of these QTLs and candidate genes have been utilized in marker-assisted breeding for achieving larger seeds and enhanced yield in legumes. Besides, the characterization of legume orthologs of candidate genes from other crops using different omics approaches helped in understanding the regulatory pathways involved in seed size determination in legumes. This review provides a direction for the effective utilization of available resources to enhance legume productivity.

A R2R3-MYB Transcription Factor of GmMYB62 Regulates Seed-Coat Color and Seed Size in Arabidopsis

The seed-coat color and seed size have an impact on both the evolutionary fitness and the grain yield of crops. Soybean is a major oil crop, and the seed-coat color and seed size exhibit natural diversity among the different soybean varieties. Here, we found an R2R3-MYB transcription factor of GmMYB62, which shows a significant increase in expression as the seed-coat color changes from yellow to black in different soybean varieties. The GmMYB62 was specifically highly expressed in reproductive organs, especially in floral organs in soybeans. The GmMYB62 encodes a nuclear protein that contains two MYB domains. In the phylogenetic analysis, the GmMYB62 was relatively conserved after the divergence of the monocots and dicots, and it also grouped with transcriptional repressors of MYBs in anthocyanin synthesis. The GmMYB62 was overexpressed in Arabidopsis and the seeds displayed a pale-brown coat in GmMYB62 overexpression lines, in contrast to the dark-brown seed coat observed in wild-type of Col-0. The anthocyanin content in the GmMYB62 overexpression lines was dramatically reduced when compared to Col-0. Additionally, the seeds in overexpression lines showed shorter lengths, larger widths, and lower thousand-seed weights than those in Col-0. Furthermore, the genes related to anthocyanin synthesis and seed size regulation were investigated, and expression of eight genes that involved in anthocyanin synthesis pathway, like chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), and anthocyanidin synthase (ANS) were severely inhibited in the GmMYB62 overexpression lines when compared to Col-0. In addition, the ARGOS-LIKE (ARL), B-Type Cyclin 1 (CYCB1), and enhancer of DA1-1 (EOD3), which govern cell expansion and proliferation, were highly expressed in GmMYB62 overexpression lines when compared to Col-0. Overall, this study sheds new light on the control of seed-coat color and seed size by GmMYB62 and provides potentially valuable targets for improving crop seed quality.

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.

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

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

Supplementary Information

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

Genome-Wide Identification, Expression and Interaction Analysis of GLN Gene Family in Soybean

As a globally significant economic crop, the seed size of soybean (Glycine max [L.] Merr.) is jointly regulated by internal genetic factors and external environmental signals. This study discovered that the GLN family proteins in soybean are similar to the KIX-PPD-MYC transcriptional repressor complex in Arabidopsis, potentially influencing seed size by regulating the expression of the downstream gene GIF1. Additionally, β-1,3-glucanase (βGlu) plays a crucial role in antifungal activity, cell composition, flower development, pollen development, abiotic resistance, seed germination, and maturation in soybean. Through a detailed analysis of the structure, chromosomal localization, phylogenetic relationships, and expression situations in different tissues at different stages of the soybean GLN gene family members, this research certifies a theoretical foundation for subsequent research on the biological functions of GLN genes in soybean. This research incorporated a comprehensive genomic identification and expression analysis of the GLN gene family in soybean. The results indicate that the 109 soybean GLN genes are unevenly distributed across soybean chromosomes and exhibit diverse expression patterns in different tissues, suggesting they may have distinct functions in soybean morphogenesis. GO enrichment analysis shows that the GLN gene family may participate in a variety of biological activities, cellular components, and molecular biological processes, particularly in catalytic activity, cellular components, and metabolic processes. These findings provide important information for comprehending the role of the GLN gene family in soybean and offer potential targets for molecular breeding of soybean.

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.

Application of genome editing in plant reproductive biology: recent advances and challenges

KEY MESSAGE

This comprehensive review underscores the application of genome editing in plant reproductive biology, including recent advances and challenges associated with it. Genome editing (GE) is a powerful technology that has the potential to accelerate crop improvement by enabling efficient, precise, and rapid engineering of plant genomes. Over the last decade, this technology has rapidly evolved from the use of meganucleases (homing endonucleases), zinc-finger nucleases, transcription activator-like effector nucleases to the use of clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (CRISPR/Cas), which has emerged as a popular GE tool in recent times and has been extensively used in several organisms, including plants. GE has been successfully employed in several crops to improve plant reproductive traits. Improving crop reproductive traits is essential for crop yields and securing the world's food supplies. In this review, we discuss the application of GE in various aspects of plant reproductive biology, including its potential application in haploid induction, apomixis, parthenocarpy, development of male sterile lines, and the regulation of self-incompatibility. We also discuss current challenges and future prospects of this technology for crop improvement, focusing on plant reproduction.

Comparative Analysis of Two Soybean Cultivars Revealed Tolerance Mechanisms Underlying Soybean Adaptation to Flooding

Flooding stress poses a significant challenge to soybean cultivation, impacting plant growth, development, and ultimately yield. In this study, we investigated the responses of two distinct soybean cultivars: flooding-tolerant Nanxiadou 38 (ND38) and flooding-sensitive Nanxiadou 45 (ND45). To achieve this, healthy seedlings were cultivated with the water surface consistently maintained at 5 cm above the soil surface. Our objective was to elucidate the physiological and molecular adaptations of the two cultivars. Under flooding stress, seedlings of both cultivars exhibited significant dwarfing and a notable decrease in root length. While there were no significant differences in the dry weight of aboveground shoots, the dry weight of underground shoots in ND38 was strikingly decreased following flooding. Additionally, total chlorophyll content decreased significantly following flooding stress, indicating impaired photosynthetic performance of the cultivars. Moreover, malondialdehyde (MDA) levels increased significantly after flooding, particularly in the ND45 cultivar, suggesting heightened oxidative stress. Expression analysis of methylation and demethylation genes indicated that MET1 and DME play crucial roles in response to flooding stress in soybeans. Meanwhile, analysis of the hemoglobin family (GLBs), aquaporin family (AQPs), glycolytic pathway-related genes, and NAC transcription factor-related genes identified GLB1-1 and GLB1-2, GLB2-2, PIP2-6, PIP2-7, TIP2-2, TIP4-1, TIP5-1, Gm02G222400 (fructose-bisphosphate aldolase), Gm19G017200 (glucose-6-phosphate isomerase), and Gm04G213900 (alcohol dehydrogenase 1) as key contributors to flooding tolerance in both soybean cultivars. These findings provide crucial insights into the physiological and molecular mechanisms underlying flooding tolerance in soybeans, which could guide future molecular breeding strategies for the development of flooding-tolerant soybean cultivars.

Advances in Soybean Genetic Improvement

The soybean (Glycine max) is a globally important crop due to its high protein and oil content, which serves as a key resource for human and animal nutrition, as well as bioenergy production. This review assesses recent advancements in soybean genetic improvement by conducting an extensive literature analysis focusing on enhancing resistance to biotic and abiotic stresses, improving nutritional profiles, and optimizing yield. We also describe the progress in breeding techniques, including traditional approaches, marker-assisted selection, and biotechnological innovations such as genetic engineering and genome editing. The development of transgenic soybean cultivars through Agrobacterium-mediated transformation and biolistic methods aims to introduce traits such as herbicide resistance, pest tolerance, and improved oil composition. However, challenges remain, particularly with respect to genotype recalcitrance to transformation, plant regeneration, and regulatory hurdles. In addition, we examined how wild soybean germplasm and polyploidy contribute to expanding genetic diversity as well as the influence of epigenetic processes and microbiome on stress tolerance. These genetic innovations are crucial for addressing the increasing global demand for soybeans, while mitigating the effects of climate change and environmental stressors. The integration of molecular breeding strategies with sustainable agricultural practices offers a pathway for developing more resilient and productive soybean varieties, thereby contributing to global food security and agricultural sustainability.

Molecular Regulation of Shoot Architecture in Soybean

Soybean (Glycine max [L.] Merr.) serves as a major source of protein and oil for humans and animals. Shoot architecture, the spatial arrangement of a plant's above-ground organs, strongly affects crop yield and is therefore a critical agronomic trait. Unlike wheat and rice crops that have greatly benefitted from the Green Revolution, soybean yield has not changed significantly in the past six decades owing to its unique shoot architecture. Soybean is a pod-bearing crop with pods adhered to the nodes, and variation in shoot architecture traits, such as plant height, node number, branch number and number of seeds per pod, directly affects the number of pods and seeds per plant, thereby determining yield. In this review, we summarize the relationship between soybean yield and these major components of shoot architecture. We also describe the latest advances in identifying the genes and molecular mechanisms underlying soybean shoot architecture and discuss possible directions and approaches for breeding new soybean varieties with ideal shoot architecture and improved yield.

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.

CRISPR/Cas genome editing in soybean: challenges and new insights to overcome existing bottlenecks

BACKGROUND

Soybean is a worldwide-cultivated crop due to its applications in the food, feed, and biodiesel industries. Genome editing in soybean began with ZFN and TALEN technologies; however, CRISPR/Cas has emerged and shortly became the preferable approach for soybean genome manipulation since it is more precise, easy to handle, and cost-effective. Recent reports have focused on the conventional Cas9 nuclease, Cas9 nickase (nCas9) derived base editors, and Cas12a (formally Cpf1) as the most commonly used genome editors in soybean. Nonetheless, several challenges in the complex plant genetic engineering pipeline need to be overcome to effectively edit the genome of an elite soybean cultivar. These challenges include (1) optimizing CRISPR cassette design (i.e., gRNA and Cas promoters, gRNA design and testing, number of gRNAs, and binary vector), (2) improving transformation frequency, (3) increasing the editing efficiency ratio of targeted plant cells, and (4) improving soybean crop production.

AIM OF REVIEW

This review provides an overview of soybean genome editing using CRISPR/Cas technology, discusses current challenges, and highlights theoretical (insights) and practical suggestions to overcome the existing bottlenecks.

KEY SCIENTIFIC CONCEPTS OF REVIEW

The CRISPR/Cas system was discovered as part of the bacterial innate immune system. It has been used as a biotechnological tool for genome editing and efficiently applied in soybean to unveil gene function, improve agronomic traits such as yield and nutritional grain quality, and enhance biotic and abiotic stress tolerance. To date, the efficiency of gRNAs has been validated using protoplasts and hairy root assays, while stable plant transformation relies on Agrobacterium-mediated and particle bombardment methods. Nevertheless, most steps of the CRISPR/Cas workflow require optimizations to achieve a more effective genome editing in soybean plants.

Morphological characterization and transcriptome analysis of rolled and narrow leaf mutant in soybean

BACKGROUND

In plants, the leaf functions as a solar panel, where photosynthesis converts carbon dioxide and water into carbohydrates and oxygen. In soybean, leaf type traits, including leaf shape, leaf area, leaf width, and leaf width so on, are considered to be associated with yield. In this study, we performed morphological characterization, transcriptome analysis, and endogenous hormone analysis of a rolled and narrow leaf mutant line (rl) in soybean.

RESULTS

Compared with wild type HX3, mutant line rl showed rolled and narrower leaflet, and smaller leaf, meanwhile rl also performed narrower pod and narrower seed. Anatomical analysis of leaflet demonstrated that cell area of upper epidermis was bigger than the cell area of lower epidermis in rl, which may lead rolled and narrow leaf. Transcriptome analysis revealed that several cytokinin oxidase/dehydrogenase (CKX) genes (Glyma.06G028900, Glyma.09G225400, Glyma.13G104700, Glyma.14G099000, and Glyma.17G054500) were up-regulation dramatically, which may cause lower cytokinin level in rl. Endogenous hormone analysis verified that cytokinin content of rl was lower. Hormone treatment results indicated that 6-BA rescued rolled leaf enough, rescued partly narrow leaf. And after 6-BA treatment, the cell area was similar between upper epidermis and lower epidermis in rl. Although IAA content and ABA content were reduced in rl, but exogenous IAA and ABA didn't affect leaf type of HX3 and rl.

CONCLUSIONS

Our results suggest abnormal cytokinin metabolism caused rolled and narrow leaf in rl, and provide valuable clues for further understanding the mechanisms underlying leaf development in soybean.

KINASE-INDUCIBLE DOMAIN INTERACTING 8 regulates helical pod morphology in Medicago truncatula

Leguminosae exhibits a wide diversity of legume forms with varying degrees of spiral morphologies, serving as an ideal clade for studying the growth and development of spiral organs. While soybean (Glycine max) develops straight pods, the pod of the model legume Medicago truncatula is a helix structure. Despite the fascinating structures and intensive description of the pods in legumes, little is known regarding the genetic mechanism underlying the highly varied spirality of the legume pods. In this study, we found that KINASE-INDUCIBLE DOMAIN INTERACTING 8 (MtKIX8) plays a key role in regulating the pod structure and spirality in M. truncatula. Unlike the coiled and barrel-shaped helix pods of the wild type, the pods of the mtkix8 mutant are loose and deformed and lose the topologic structure as observed in the wild-type pods. In the pods of the mtkix8 mutant, the cells proliferate more actively and overly expand, particularly in the ventral suture, resulting in uncoordinated growth along the dorsal and ventral sutures of pods. The core cell cycle genes CYCLIN D3s are upregulated in the mtkix8 pods, leading to the prolonged growth of the ventral suture region of the pods. Our study revealed the key role of MtKIX8 in regulating seed pod development in M. truncatula and demonstrates a genetic regulatory model underlying the establishment of the helical pod in legumes.

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.

Genome-wide identification of the TIFY gene family in tobacco and expression analysis in response to Ralstonia solanacearum infection

The TIFY gene family plays an essential role in plant development and abiotic and biotic stress responses. In this study, genome-wide identification of TIFY members in tobacco and their expression pattern analysis in response to Ralstonia solanacearum infection were performed. A total of 33 TIFY genes were identified, including the TIFY, PPD, ZIM&ZML and JAZ subfamilies. Promoter analysis results indicated that a quantity of light-response, drought-response, SA-response and JA-response cis-elements exist in promoter regions. The TIFY gene family exhibited expansion and possessed gene redundancy resulting from tobacco ploidy change. In addition, most NtTIFYs equivalently expressed in roots, stems and leaves, while NtTIFY1, NtTIFY4, NtTIFY18 and NtTIFY30 preferentially expressed in roots. The JAZ III clade showed significant expression changes after inoculation with R. solanacearum, and the expression of NtTIFY7 in resistant varieties, compared with susceptible varieties, was more stably induced. Furthermore, NtTIFY7-silenced plants, compared with the control plants, were more susceptible to bacterial wilt. These results lay a foundation for exploring the evolutionary history of TIFY gene family and revealing gene function of NtTIFYs in tobacco bacterial wilt resistance.

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.

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.

Leaf growth - complex regulation of a seemingly simple process

Understanding the underlying mechanisms of plant development is crucial to successfully steer or manipulate plant growth in a targeted manner. Leaves, the primary sites of photosynthesis, are vital organs for many plant species, and leaf growth is controlled by a tight temporal and spatial regulatory network. In this review, we focus on the genetic networks governing leaf cell proliferation, one major contributor to final leaf size. First, we provide an overview of six regulator families of leaf growth in Arabidopsis: DA1, PEAPODs, KLU, GRFs, the SWI/SNF complexes, and DELLAs, together with their surrounding genetic networks. Next, we discuss their evolutionary conservation to highlight similarities and differences among species, because knowledge transfer between species remains a big challenge. Finally, we focus on the increase in knowledge of the interconnectedness between these genetic pathways, the function of the cell cycle machinery as their central convergence point, and other internal and environmental cues.

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.

Exploring the role of FBXL fbxl gene family in Soybean: Implications for plant height and seed size regulation

F-box proteins constitute a significant family in eukaryotes and, as a component of the Skp1p-cullin-F-box complex, are considered critical for cellular protein degradation and other biological processes in plants. Despite their importance, the functions of F-box proteins, particularly those with C-terminal leucine-rich repeat (LRR) domains, remain largely unknown in plants. Therefore, the present study conducted genome-wide identification and in silico characterization of F-BOX proteins with C-terminal LRR domains in soybean (Glycine max L.) (GmFBXLs). A total of 45 GmFBXLs were identified. The phylogenetic analysis showed that GmFBXLs could be subdivided into ten subgroups and exhibited a close relationship with those from Arabidopsis thaliana, Cicer aretineum, and Medicago trunculata. It was observed that most cis-regulatory elements in the promoter regions of GmFBXLs are involved in hormone signalling, stress responses, and developmental stages. In silico transcriptome data illustrated diverse expression patterns of the identified GmFBXLs across various tissues, such as shoot apical meristem, flower, green pods, leaves, nodules, and roots. Overexpressing (OE) GmFBXL12 in Tianlong No.1 cultivar resulted in a significant difference in seed size, number of pods, and number of seeds per plant, indicated a potential increase in yield compared to wild type. This study offers valuable perspectives into the role of FBXLs in soybean, serving as a foundation for future research. Additionally, the identified OE lines represent valuable genetic resources for enhancing seed-related traits in soybean.

Regulation of seed traits in soybean

Soybean (Glycine max) is an essential economic crop that provides vegetative oil and protein for humans, worldwide. Increasing soybean yield as well as improving seed quality is of great importance. Seed weight/size, oil and protein content are the three major traits determining seed quality, and seed weight also influences soybean yield. In recent years, the availability of soybean omics data and the development of related techniques have paved the way for better research on soybean functional genomics, providing a comprehensive understanding of gene functions. This review summarizes the regulatory genes that influence seed size/weight, oil content and protein content in soybean. We also provided a general overview of the pleiotropic effect for the genes in controlling seed traits and environmental stresses. Ultimately, it is expected that this review will be beneficial in breeding improved traits in soybean.

Global analysis of seed transcriptomes reveals a novel PLATZ regulator for seed size and weight control in soybean

Seed size and weight are important factors that influence soybean yield. Combining the weighted gene co-expression network analysis (WGCNA) of 45 soybean accessions and gene dynamic changes in seeds at seven developmental stages, we identified candidate genes that may control the seed size/weight. Among these, a PLATZ-type regulator overlapping with 10 seed weight QTLs was further investigated. This zinc-finger transcriptional regulator, named as GmPLATZ, is required for the promotion of seed size and weight in soybean. The GmPLATZ may exert its functions through direct binding to the promoters and activation of the expression of cyclin genes and GmGA20OX for cell proliferation. Overexpression of the GmGA20OX enhanced seed size/weight in soybean. We further found that the GmPLATZ binds to a 32-bp sequence containing a core palindromic element AATGCGCATT. Spacing of the flanking sequences beyond the core element facilitated GmPLATZ binding. An elite haplotype Hap3 was also identified to have higher promoter activity and correlated with higher gene expression and higher seed weight. Orthologues of the GmPLATZ from rice and Arabidopsis play similar roles in seeds. Our study reveals a novel module of GmPLATZ-GmGA20OX/cyclins in regulating seed size and weight and provides valuable targets for breeding of crops with desirable agronomic traits.

Molecular, genetic, and genomic basis of seed size and yield characteristics in soybean

Soybean (Glycine max L. Merr.) is a crucial oilseed cash crop grown worldwide and consumed as oil, protein, and food by humans and feed by animals. Comparatively, soybean seed yield is lower than cereal crops, such as maize, rice, and wheat, and the demand for soybean production does not keep up with the increasing consumption level. Therefore, increasing soybean yield per unit area is the most crucial breeding objective and is challenging for the scientific community. Moreover, yield and associated traits are extensively researched in cereal crops, but little is known about soybeans' genetics, genomics, and molecular regulation of yield traits. Soybean seed yield is a complex quantitative trait governed by multiple genes. Understanding the genetic and molecular processes governing closely related attributes to seed yield is crucial to increasing soybean yield. Advances in sequencing technologies have made it possible to conduct functional genomic research to understand yield traits' genetic and molecular underpinnings. Here, we provide an overview of recent progress in the genetic regulation of seed size in soybean, molecular, genetics, and genomic bases of yield, and related key seed yield traits. In addition, phytohormones, such as auxin, gibberellins, cytokinins, and abscisic acid, regulate seed size and yield. Hence, we also highlight the implications of these factors, challenges in soybean yield, and seed trait improvement. The information reviewed in this study will help expand the knowledge base and may provide the way forward for developing high-yielding soybean cultivars for future food demands.

A loss-of-function mutant allele of a glycosyl hydrolase gene has been co-opted for seed weight control during soybean domestication

The resultant DNA from loss-of-function mutation can be recruited in biological evolution and development. Here, we present such a rare and potential case of "to gain by loss" as a neomorphic mutation during soybean domestication for increasing seed weight. Using a population derived from a chromosome segment substitution line of Glycine max (SN14) and Glycine soja (ZYD06), a quantitative trait locus (QTL) of 100-seed weight (qHSW) was mapped on chromosome 11, corresponding to a truncated β-1, 3-glucosidase (βGlu) gene. The novel gene hsw results from a 14-bp deletion, causing a frameshift mutation and a premature stop codon in the βGlu. In contrast to HSW, the hsw completely lost βGlu activity and function but acquired a novel function to promote cell expansion, thus increasing seed weight. Overexpressing hsw instead of HSW produced large soybean seeds, and surprisingly, truncating hsw via gene editing further increased the seed size. We further found that the core 21-aa peptide of hsw and its variants acted as a promoter of seed size. Transcriptomic variation in these transgenic soybean lines substantiated the integration hsw into cell and seed size control. Moreover, the hsw allele underwent selection and expansion during soybean domestication and improvement. Our work cloned a likely domesticated QTL controlling soybean seed weight, revealed a novel genetic variation and mechanism in soybean domestication, and provided new insight into crop domestication and breeding, and plant evolution.

Genome-wide analysis of KIX gene family for organ size regulation in soybean (Glycine max L.)

The KIX domain, conserved among various nuclear and co-activator factors, acts as a binding site that interacts with other transcriptional activators and co-activators, playing a crucial role in gene expression regulation. In plants, the KIX domain is involved in plant hormone signaling, stress response regulation, cell cycle control, and differentiation, indicating its potential relevance to crop productivity. This study aims to identify and characterize KIX domains within the soybean (Glycine max L.) genome to predict their potential role in improving crop productivity. The conservation and evolutionary history of the KIX domains were explored in 59 plant species, confirming the presence of the KIX domains in diverse plants. Specifically, 13 KIX domains were identified within the soybean genome and classified into four main groups, namely GmKIX8/9, GmMED15, GmHAC, and GmRECQL, through sequence alignment, structural analysis, and phylogenetic tree construction. Association analysis was performed between KIX domain haplotypes and soybean seed-related agronomic traits using re-sequencing data from a core collection of 422 accessions. The results revealed correlations between SNP variations observed in GmKIX8-3 and GmMED15-4 and soybean seed phenotypic traits. Additionally, transcriptome analysis confirmed significant expression of the KIX domains during the early stages of soybean seed development. This study provides the first characterization of the structural, expression, genomic haplotype, and molecular features of the KIX domain in soybean, offering a foundation for functional analysis of the KIX domain in soybean and other plants.

Recent advances in the improvement of soybean seed traits by genome editing

Genetic improvement of soybean seed traits is important for developing new varieties that meet the demand for soybean as a food, forage crop, and industrial products. A large number of soybean genome sequences are currently publicly available. This genome sequence information provides a significant opportunity to design genomic approaches to improve soybean traits. Genome editing represents a major advancement in biotechnology. The production of soybean mutants through genome editing is commonly achieved with either an Agrobacterium-mediated or biolistic transformation platform, which have been optimized for various soybean genotypes. Currently, the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated endonuclease 9 (Cas9) system, which represents a major advance in genome editing, is used to improve soybean traits, such as fatty acid composition, protein content and composition, flavor, digestibility, size, and seed-coat color. In this review, we summarize the recent advances in the improvement of soybean seed traits through genome editing. We also discuss the characteristics of genome editing using the CRISPR/Cas9 system with transformation platforms.

The MIO1-MtKIX8 module regulates the organ size in Medicago truncatula

Plant organ size is an important agronomic trait tightly related to crop yield. However, the molecular mechanisms underlying organ size regulation remain largely unexplored in legumes. We previously characterized a key regulator F-box protein MINI ORGAN1 (MIO1)/SMALL LEAF AND BUSHY1 (SLB1), which controls plant organ size in the model legume Medicago truncatula. In order to further dissect the molecular mechanism, MIO1 was used as the bait to screen its interacting proteins from a yeast library. Subsequently, a KIX protein, designated MtKIX8, was identified from the candidate list. The interaction between MIO1 and MtKIX8 was confirmed further by Y2H, BiFC, split-luciferase complementation and pull-down assays. Phylogenetic analyses indicated that MtKIX8 is highly homologous to Arabidopsis KIX8, which negatively regulates organ size. Moreover, loss-of-function of MtKIX8 led to enlarged leaves and seeds, while ectopic expression of MtKIX8 in Arabidopsis resulted in decreased cotyledon area and seed weight. Quantitative reverse-transcription PCR and in situ hybridization showed that MtKIX8 is expressed in most developing organs. We also found that MtKIX8 serves as a crucial molecular adaptor, facilitating interactions with BIG SEEDS1 (BS1) and MtTOPLESS (MtTPL) proteins in M. truncatula. Overall, our results suggest that the MIO1-MtKIX8 module plays a significant and conserved role in the regulation of plant organ size. This module could be a good target for molecular breeding in legume crops and forages.

Critical role for uricase and xanthine dehydrogenase in soybean nitrogen fixation and nodule development

De novo purine biosynthesis is required for the incorporation of fixed nitrogen in ureide exporting nodules, as formed on soybean [Glycine max (L.) Merr.] roots. However, in many cases, the enzymes involved in this pathway have been deduced strictly from genome annotations with little direct genetic evidence, such as mutant studies, to confirm their biochemical function or importance to nodule development. While efforts to develop large mutant collections of soybean are underway, research on this plant is still hampered by the inability to obtain mutations in any specific gene of interest. Using a forward genetic approach, as well as CRISPR/Cas9 gene editing via Agrobacterium rhizogenes-mediated hairy root transformation, we identified and characterized the role of GmUOX (Uricase) and GmXDH (Xanthine Dehydrogenase) in nitrogen fixation and nodule development in soybean. The gmuox knockout soybean mutants displayed nitrogen deficiency chlorosis and early nodule senescence, as exemplified by the reduced nitrogenase (acetylene reduction) activity in nodules, the internal greenish-white internal appearance of nodules, and diminished leghemoglobin production. In addition, gmuox1 nodules showed collapsed infected cells with degraded cytoplasm, aggregated bacteroids with no discernable symbiosome membranes, and increased formation of poly-β-hydroxybutyrate granules. Similarly, knockout gmxdh mutant nodules, generated with the CRISPR/Cas9 system, also exhibited early nodule senescence. These genetic studies confirm the critical role of the de novo purine metabolisms pathway not only in the incorporation of fixed nitrogen but also in the successful development of a functional, nitrogen-fixing nodule. Furthermore, these studies demonstrate the great utility of the CRISPR/Cas9 system for studying root-associated gene traits when coupled with hairy root transformation.

The elite variations in germplasms for soybean breeding

The genetic base of soybean cultivars (Glycine max (L.) Merr.) has been narrowed through selective domestication and specific breeding improvement, similar to other crops. This presents challenges in breeding new cultivars with improved yield and quality, reduced adaptability to climate change, and increased susceptibility to diseases. On the other hand, the vast collection of soybean germplasms offers a potential source of genetic variations to address those challenges, but it has yet to be fully leveraged. In recent decades, rapidly improved high-throughput genotyping technologies have accelerated the harness of elite variations in soybean germplasm and provided the important information for solving the problem of a narrowed genetic base in breeding. In this review, we will overview the situation of maintenance and utilization of soybean germplasms, various solutions provided for different needs in terms of the number of molecular markers, and the omics-based high-throughput strategies that have been used or can be used to identify elite alleles. We will also provide an overall genetic information generated from soybean germplasms in yield, quality traits, and pest resistance for molecular breeding.

Genetic regulatory networks of soybean seed size, oil and protein contents

As a leading oilseed crop that supplies plant oil and protein for daily human life, increasing yield and improving nutritional quality (high oil or protein) are the top two fundamental goals of soybean breeding. Seed size is one of the most critical factors determining soybean yield. Seed size, oil and protein contents are complex quantitative traits governed by genetic and environmental factors during seed development. The composition and quantity of seed storage reserves directly affect seed size. In general, oil and protein make up almost 60% of the total storage of soybean seed. Therefore, soybean's seed size, oil, or protein content are highly correlated agronomical traits. Increasing seed size helps increase soybean yield and probably improves seed quality. Similarly, rising oil and protein contents improves the soybean's nutritional quality and will likely increase soybean yield. Due to the importance of these three seed traits in soybean breeding, extensive studies have been conducted on their underlying quantitative trait locus (QTLs) or genes and the dissection of their molecular regulatory pathways. This review summarized the progress in functional genome controlling soybean seed size, oil and protein contents in recent decades, and presented the challenges and prospects for developing high-yield soybean cultivars with high oil or protein content. In the end, we hope this review will be helpful to the improvement of soybean yield and quality in the future breeding process.

Novel Seed Size: A Novel Seed-Developing Gene in Glycine max

Soybean-seed development is controlled in multiple ways, as in many known regulating genes. Here, we identify a novel gene, Novel Seed Size (NSS), involved in seed development, by analyzing a T-DNA mutant (S006). The S006 mutant is a random mutant of the GmFTL4pro:GUS transgenic line, with phenotypes with small and brown seed coats. An analysis of the metabolomics and transcriptome combined with RT-qPCR in the S006 seeds revealed that the brown coat may result from the increased expression of chalcone synthase 7/8 genes, while the down-regulated expression of NSS leads to small seed size. The seed phenotypes and a microscopic observation of the seed-coat integument cells in a CRISPR/Cas9-edited mutant nss1 confirmed that the NSS gene conferred small phenotypes of the S006 seeds. As mentioned in an annotation on the Phytozome website, NSS encodes a potential DNA helicase RuvA subunit, and no such genes were previously reported to be involved in seed development. Therefore, we identify a novel gene in a new pathway controlling seed development in soybeans.

Understandings and future challenges in soybean functional genomics and molecular breeding

Soybean (Glycine max) is a major source of plant protein and oil. Soybean breeding has benefited from advances in functional genomics. In particular, the release of soybean reference genomes has advanced our understanding of soybean adaptation to soil nutrient deficiencies, the molecular mechanism of symbiotic nitrogen (N) fixation, biotic and abiotic stress tolerance, and the roles of flowering time in regional adaptation, plant architecture, and seed yield and quality. Nevertheless, many challenges remain for soybean functional genomics and molecular breeding, mainly related to improving grain yield through high-density planting, maize-soybean intercropping, taking advantage of wild resources, utilization of heterosis, genomic prediction and selection breeding, and precise breeding through genome editing. This review summarizes the current progress in soybean functional genomics and directs future challenges for molecular breeding of soybean.

Identification of major quantitative trait loci and candidate genes for seed weight in soybean

Four major quantitative trait loci for 100-seed weight were identified in a soybean RIL population under five environments, and the most likely candidate genes underlying these loci were identified. Seed weight is an important target of soybean breeding. However, the genes underlying the major quantitative trait loci (QTL) controlling seed weight remain largely unknown. In this study, a soybean population of 300 recombinant inbred lines (RILs) derived from a cross between PI595843 (PI) and WH was used to map the QTL and identify candidate genes for seed weight. The RIL population was genotyped through whole genome resequencing, and phenotyped for 100-seed weight under five environments. A total of 38 QTL were detected, and four major QTL, each explained at least 10% of the variation in 100-seed weight, were identified. Six candidate genes within these four major QTL regions were identified by analyses of their tissue expression patterns, gene annotations, and differential gene expression levels in soybean seeds during four developmental stages between two parental lines. Further sequence variation analyses revealed a C to T substitution in the first exon of the Glyma.19G143300, resulting in an amino acid change between PI and WH, and thus leading to a different predicted kinase domain, which might affect its protein function. Glyma.19G143300 is highly expressed in soybean seeds and encodes a leucine-rich repeat receptor-like protein kinase (LRR-RLK). Its predicted protein has typical domains of LRR-RLK family, and phylogenetic analyses reveled its similarity with the known LRR-RLK protein XIAO (LOC_Os04g48760), which is involved in controlling seed size. The major QTL and candidate genes identified in this study provide useful information for molecular breeding of new soybean cultivars with desirable seed weight.

Systematic trait dissection in oilseed rape provides a comprehensive view, further insight, and exact roadmap for yield determination

Background

Yield is the most important and complex trait that is influenced by numerous relevant traits with very complicated interrelations. While there are a large number of studies on the phenotypic relationship and genetic basis of yield traits, systematic studies with further dissection focusing on yield are limited. Therefore, there is still lack of a comprehensive and in-depth understanding of the determination of yield.

Results

In this study, yield was systematically dissected at the phenotypic, genetic to molecular levels in oilseed rape (Brassica napus L.). The analysis of correlation, network, and principal component for 21 traits in BnaZN-RIL population showed that yield was determined by a complex trait network with key contributors. The analysis of the constructed high-density single nucleotide polymorphism (SNP) linkage map revealed the concentrated distribution of distorted and heterozygous markers, likely due to selection on genes controlling the growth period and yield heterosis. A total of 134 consensus quantitative trait loci (QTL) were identified for 21 traits, of which all were incorporated into an interconnecting QTL network with dozens of hub-QTL. Four representative hub-QTL were further dissected to the target or candidate genes that governed the causal relationships between the relevant traits.

Conclusions

The highly consistent results at the phenotypic, genetic, and molecular dissecting demonstrated that yield was determined by a multilayer composite network that involved numerous traits and genes showing complex up/down-stream and positive/negative regulation. This provides a systematic view, further insight, and exact roadmap for yield determination, which represents a significant advance toward the understanding and dissection of complex traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02134-w.

An enhancing effect attributed to a nonsynonymous mutation in SOYBEAN SEED SIZE 1, a SPINDLY-like gene, is exploited in soybean domestication and improvement

Soybean (Glycine max) was domesticated from its wild relative Glycine soja. One-hundred-seed weight is one of the most important domesticated traits determining soybean yield; however, its underlying genetic basis remains elusive. We characterized a soybean seed size 1 (sss1) mutant featuring large seeds compared to its wild-type background. Positional cloning revealed that the candidate gene GmSSS1 encoded a SPINDLY homolog and was co-located in a well-identified quantitative trait locus (QTL)-rich region on chromosome 19. Knocking out GmSSS1 resulted in small seeds, while overexpressing GmSSS1/Gmsss1 induced large seeds. Modulating GmSSS1/Gmsss1 in transgenic plants can positively influence cell expansion and cell division. Relative to GmSSS1, one mutation leading to an E to Q substitution at the 182^nd residue in Gmsss1 conferred an enhancing effect on seed weight. GmSSS1 underwent diversification in wild-type and cultivated soybean, and the alleles encoding the Gmsss1-type substitution of 182^nd -Q, which originated along the central and downstream parts of the Yellow River, were selected and expanded during soybean domestication and improvement. We cloned the causative gene for the sss1 mutant, which is linked with a seed weight QTL, identified an elite allele of this gene for increasing seed weight, and provided new insights into soybean domestication and breeding.

Developing Genetic Engineering Techniques for Control of Seed Size and Yield

Many signaling pathways regulate seed size through the development of endosperm and maternal tissues, which ultimately results in a range of variations in seed size or weight. Seed size can be determined through the development of zygotic tissues (endosperm and embryo) and maternal ovules. In addition, in some species such as rice, seed size is largely determined by husk growth. Transcription regulator factors are responsible for enhancing cell growth in the maternal ovule, resulting in seed growth. Phytohormones induce significant effects on entire features of growth and development of plants and also regulate seed size. Moreover, the vegetative parts are the major source of nutrients, including the majority of carbon and nitrogen-containing molecules for the reproductive part to control seed size. There is a need to increase the size of seeds without affecting the number of seeds in plants through conventional breeding programs to improve grain yield. In the past decades, many important genetic factors affecting seed size and yield have been identified and studied. These important factors constitute dynamic regulatory networks governing the seed size in response to environmental stimuli. In this review, we summarized recent advances regarding the molecular factors regulating seed size in Arabidopsis and other crops, followed by discussions on strategies to comprehend crops' genetic and molecular aspects in balancing seed size and yield.

Genome-wide association studies provide genetic insights into natural variation of seed-size-related traits in mungbean

Although mungbean (Vigna radiata (L.) R. Wilczek) is an important legume crop, its seed yield is relatively low. To address this issue, here 196 accessions with 3,607,508 SNP markers were used to identify quantitative trait nucleotides (QTNs), QTN-by-environment interactions (QEIs), and their candidate genes for seed length (SL), seed width, and 100-seed weight (HSW) in two environments. As a result, 98 QTNs and 20 QEIs were identified using 3VmrMLM, while 95, >10,000, and 15 QTNs were identified using EMMAX, GEMMA, and CMLM, respectively. Among 809 genes around these QTNs, 12 were homologous to known seed-development genes in rice and Arabidopsis thaliana, in which 10, 2, 1, and 0 genes were found, respectively, by the above four methods to be associated with the three traits, such as VrEmp24/25 for SL and VrKIX8 for HSW. Eight of the 12 genes were significantly differentially expressed between two large-seed and two small-seed accessions, and VrKIX8, VrPAT14, VrEmp24/25, VrIAR1, VrBEE3, VrSUC4, and Vrflo2 were further verified by RT-qPCR. Among 65 genes around these QEIs, VrFATB, VrGSO1, VrLACS2, and VrPAT14 were homologous to known seed-development genes in A. thaliana, although new experiments are necessary to explore these novel GEI-trait associations. In addition, 54 genes were identified in comparative genomics analysis to be associated with seed development pathway, in which VrKIX8, VrABA2, VrABI5, VrSHB1, and VrIKU2 were also identified in genome-wide association studies. This result provided a reliable approach for identifying seed-size-related genes in mungbean and a solid foundation for further molecular biology research on seed-size-related genes.

Natural allelic variation of GmST05 controlling seed size and quality in soybean

Seed size is one of the most important agronomic traits determining the yield of crops. Cloning the key genes controlling seed size and pyramiding their elite alleles will facilitate yield improvement. To date, few genes controlling seed size have been identified in soybean, a major crop that provides half of the plant oil and one quarter of the plant protein globally. Here, through a genome-wide association study of over 1800 soybean accessions, we determined that natural allelic variation at GmST05 (Seed Thickness 05) predominantly controlled seed thickness and size in soybean germplasm. Further analyses suggested that the two major haplotypes of GmST05 differed significantly at the transcriptional level. Transgenic experiments demonstrated that GmST05 positively regulated seed size and influenced oil and protein contents, possibly by regulating the transcription of GmSWEET10a. Population genetic diversity analysis suggested that allelic variations of GmST05 were selected during geographical differentiation but have not been fixed. In summary, natural variation in GmST05 determines transcription levels and influences seed size and quality in soybean, making it an important gene resource for soybean molecular breeding.

PEAPOD repressors modulate and coordinate developmental responses to light intensity in Arabidopsis

Higher plants adapt to different light intensities by altering hypocotyl elongation, stomatal density, seed size, and flowering time. Despite the importance of this developmental plasticity, knowledge of the underlying genetic and molecular mechanisms modulating and coordinating responses to light intensity remains incomplete. Here, I report that in Arabidopsis the PEAPOD (PPD) repressors PPD1 and PPD2 prevent exaggerated responses to light intensity. Genetic and transcriptome analyses, of a ppd deletion mutant and a PPD1 overexpression genotype, were used to identify how PPD repressors modulate the light signalling network. A ppd1/ppd2 deletion mutant has elongated hypocotyls, elevated stomatal density, enlarged seed, and delayed flowering, whereas overexpression of PPD1 results in the reverse. Transcription of both PPD1 and PPD2, upregulated in low light and downregulated in higher light, is activated by PHYTOCHROME INTERACTING FACTOR 4. I found PPDs modulate light signalling by negative regulation of SUPPRESSOR OF phyA-105 (SPA1) transcription. Whereas PPDs coordinate many of the responses to light intensity - hypocotyl elongation, flowering time, and stomatal density - by repression/de-repression of SPA1, PPD regulation of seed size occurs independent of SPA1. In conclusion PPD repressors modulate and coordinate developmental responses to light intensity by altering light signal transduction.

Experimental Evidence for Seed Metabolic Allometry in Barrel Medic (Medicago truncatula Gaertn.)

Seed size is often considered to be an important trait for seed quality, i.e., vigour and germination performance. It is believed that seed size reflects the quantity of reserve material and thus the C and N sources available for post-germinative processes. However, mechanisms linking seed size and quality are poorly documented. In particular, specific metabolic changes when seed size varies are not well-known. To gain insight into this aspect, we examined seed size and composition across different accessions of barrel medic (Medicago truncatula Gaertn.) from the genetic core collection. We conducted multi-elemental analyses and isotope measurements, as well as exact mass GC–MS metabolomics. There was a systematic increase in N content (+0.17% N mg−1) and a decrease in H content (–0.14% H mg−1) with seed size, reflecting lower lipid and higher S-poor protein quantity. There was also a decrease in 2H natural abundance (δ2H), due to the lower prevalence of 2H-enriched lipid hydrogen atoms that underwent isotopic exchange with water during seed development. Metabolomics showed that seed size correlates with free amino acid and hexoses content, and anticorrelates with amino acid degradation products, disaccharides, malic acid and free fatty acids. All accessions followed the same trend, with insignificant differences in metabolic properties between them. Our results show that there is no general, proportional increase in metabolite pools with seed size. Seed size appears to be determined by metabolic balance (between sugar and amino acid degradation vs. utilisation for storage), which is in turn likely determined by phloem source metabolite delivery during seed development.

Development of Ethyl Methanesulfonate Mutant Edamame Soybean (Glycine max (L.) Merr.) Populations and Forward and Reverse Genetic Screening for Early-Flowering Mutants

Induced mutation is a viable breeding strategy that is widely utilized in the development of elite plant varieties. We aimed to improve a variety of edamame by constructing novel mutant populations using the ethyl methanesulfonate in soybeans (Glycine max (L.) Merr.). In the M₂ population, the flowering stage showed a considerable standard deviation compared to the wild type, confirming that the mutant populations had the expected DNA mutations. To identify the DNA mutations in the mutant populations, we used the targeting induced local lesions in genomes (TILLING) method, which is a reverse genetic method, to search for soybean flowering-related gene mutants. A total of 30 mutants from E1, E3, E4, and PhyA1 genes, which are known to be highly effective genes, or their homologous gene for flowering and maturation found in soybean quantitative trait locus analyses were isolated from our TILLING screening. Among these mutants, there were eleven nonsynonymous substitution mutants, one nonsense mutant, and two single nucleotide deletion mutants that could be expected to reduce or eliminate gene function. The e1, e3, and e4 mutants obtained in this study flowered considerably earlier than the wild type. In particular, the e1 mutant with a nonsynonymous substitution flowered approximately 1 month after sowing regardless of the sowing date, and its harvest date was approximately 1 month earlier than that of the wild type. Mutations identified using the TILLING method could not only be used as gel-based DNA markers with the same manipulation method, but the mutations could also be detected as DNA markers by the high-resolution melting method. These results indicate that mutations achieved without chromosome modification by crossbreeding are effective for early and practical improvement of superior varieties and that efficient selection of mutants by reverse genetics is an effective method for the identification of genetic modifications. The edamame mutant populations developed in this study are believed to possess various useful alleles which may be applicable in the search for mutations that lead to improved edamame yield and eating quality beyond the flowering stage.

Downregulation of a gibberellin 3β-hydroxylase enhances photosynthesis and increases seed yield in soybean

Seed yield, determined mainly by seed numbers and seed weight, is the primary target of soybean breeding. Identifying the genes underlying yield-related traits is of great significance. Through joint linkage mapping and a genome-wide association study for 100-seed weight, we cloned GmGA3ox1, a gene encoding gibberellin 3β-hydroxylase, which is the key enzyme in the gibberellin synthesis pathway. Genome resequencing identified a beneficial GmGA3ox1 haplotype contributing to high seed weight, which was further confirmed by soybean transformants. CRISPR/Cas9-generated gmga3ox1 mutants showed lower seed weight, but promoted seed yield by increasing seed numbers. The gmga3ox1 mutants reduced gibberellin biosynthesis while enhancing photosynthesis. Knockout of GmGA3ox1 resulted in the upregulation of numerous photosynthesis-related genes, particularly the GmRCA family encoding ribulose-1,5-bispho-sphate carboxylase-oxygenase (Rubisco) activases. The basic leucine zipper transcription factors GmbZIP97 and GmbZIP159, which were both upregulated in the gmga3ox1 mutants and induced by the gibberellin synthesis inhibitor uniconazole, could bind to the promoter of GmRCAβ and activate its expression. Analysis of genomic sequences with over 2700 soybean accessions suggested that GmGA3ox1 is being gradually utilized in modern breeding. Our results elucidated the important role of GmGA3ox1 in soybean yield. These findings reveal important clues for future high-yield breeding in soybean and other crops.

Identification of closely associated SNPs and candidate genes with seed size and shape via deep re-sequencing GWAS in soybean

KEY MESSAGE

A soybean natural population was genotyped by deep re-sequencing and phenotyped for six seed size- and shape-related traits under six environments to identify closely associated SNPs and candidate genes. Seed size and shape are important determining factors for soybean yield formation, while their genetic basis and molecular mechanism are still largely unknown, which seriously constrains the increasing of soybean yield at present. In view of this, a natural population was genotyped via the deep re-sequencing technique (~ 20 ×) and phenotyped for six related traits under six environments. In total, 154 SNPs were closely associated with seed length across diverse environments, and 323, 483, 565, 394 and 2038 SNPs were closely associated with seed width, seed diameter, seed circumference, seed area and ratio of length to width under multiple environments. Moreover, 98.70%, 96.28%, 48.24%, 85.13%, 97.21% and 98.58% of them were further demonstrated by the BLUP and mean values of the related traits. Furthermore, 218 genes flanking the associated SNPs on chromosomes 6 and 10 were analyzed for DNA mutations and RNA expressions through SNP alleles and transcriptome data, simultaneously. The candidate genes, Glyma.10G035200 (Sn1-specific diacylglycerol lipase), Glyma.10G035400 (transcription factor) and Glyma.10G058200 (phenylalanine ammonia-lyase), were discovered to relate with the seed size and shape for their different DNA sequences or differential RNA expressions among soybean varieties at five seed developmental stages. Thus, these closely associated SNPs and related genes provide novel insights and useful information for the seed size and shape genetic basis dissection and breeding improvement in soybean.

Identification of major genomic regions for soybean seed weight by genome-wide association study

The hundred-seed weight (HSW) is an important yield component and one of the principal breeding traits in soybean. More than 250 quantitative trait loci (QTL) for soybean HSW have been identified. However, most of them have a large genomic region or are environmentally sensitive, which provide limited information for improving the phenotype in marker-assisted selection (MAS) and identifying the candidate genes. Here, we utilized 281 soybean accessions with 58,112 single nucleotide polymorphisms (SNPs) to dissect the genetic basis of HSW in across years in the northern Shaanxi province of China through one single-locus (SL) and three multi-locus (ML) genome-wide association study (GWAS) models. As a result, one hundred and fifty-four SNPs were detected to be significantly associated with HSW in at least one environment via SL-GWAS model, and 27 of these 154 SNPs were detected in all (three) environments and located within 7 linkage disequilibrium (LD) block regions with the distance of each block ranged from 40 to 610 Kb. A total of 15 quantitative trait nucleotides (QTNs) were identified by three ML-GWAS models. Combined with the results of different GWAS models, the 7 LD block regions associated with HSW detected by SL-GWAS model could be verified directly or indirectly by the results of ML-GWAS models. Eleven candidate genes underlying the stable loci that may regulate seed weight in soybean were predicted. The significantly associated SNPs and the stable loci as well as predicted candidate genes may be of great importance for marker-assisted breeding, polymerization breeding, and gene discovery for HSW in soybean.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01310-y.

The Organ Size and Morphological Change During the Domestication Process of Soybean

Soybean is one of the most important legume crops that can provide the rich source of protein and oil for human beings and livestock. In the twenty-one century, the total production of soybean is seriously behind the needs of a growing world population. Cultivated soybean [Glycine max (L.) Merr.] was domesticated from wild soybean (G. soja Sieb. and Zucc.) with the significant morphology and organ size changes in China around 5,000 years ago, including twisted stems to erect stems, small seeds to large seeds. Then it was spread worldwide to become one of the most popular and important crops. The release of the reference soybean genome and omics data provides powerful tools for researchers and breeders to dissect the functional genes and apply the germplasm in their work. Here, we summarized the function genes related to yield traits and organ size in soybean, including stem growth habit, leaf size and shape, seed size and weight. In addition, we also summarized the selection of organ traits during soybean domestication. In the end, we also discussed the application of new technology including the gene editing on the basic research and breeding of soybean, and the challenges and research hotspots in the future.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

SlKIX8 and SlKIX9 are negative regulators of leaf and fruit growth in tomato

Plant organ size and shape are major agronomic traits that depend on cell division and expansion, which are both regulated by complex gene networks. In several eudicot species belonging to the rosid clade, organ growth is controlled by a repressor complex consisting of PEAPOD (PPD) and KINASE-INDUCIBLE DOMAIN INTERACTING (KIX) proteins. The role of these proteins in asterids, which together with the rosids constitute most of the core eudicot species, is unknown. We used Clustered Regularly Interspaced Short Palindromic Repeats-CRISPR-associated protein 9 genome editing to target SlKIX8 and SlKIX9 in the asterid model species tomato (Solanum lycopersicum) and analyzed loss-of-function phenotypes. Loss-of-function of SlKIX8 and SlKIX9 led to the production of enlarged, dome-shaped leaves and these leaves exhibited increased expression of putative Solanum lycopersicum PPD (SlPPD target genes. Unexpectedly, kix8 kix9 mutants carried enlarged fruits with increased pericarp thickness due to cell expansion. At the molecular level, protein interaction assays indicated that SlKIX8 and SlKIX9 act as adaptors between the SlPPD and SlTOPLESS co-repressor proteins. Our results show that KIX8 and KIX9 are regulators of organ growth in asterids and can be used in strategies to improve important traits in produce such as thickness of the fruit flesh.