Dt2 is a gain-of-function MADS-domain factor gene that specifies semideterminacy in soybean

The diversification of the shoot branching system: A quantitative and comparative perspective in meristem determinacy

Reiterative shoot branching largely defines important yield components of crops and is essentially controlled by programs that direct the initiation, dormancy release, and differentiation of meristems in the axils of leaves. Here, we focus on meristem determinacy, defining the number of reiterations that shape the shoot architectures and exhibit enormous diversity in a wide range of species. The meristem determinacy per se is hierarchically complex and context-dependent for the successively emerged meristems, representing a crucial mechanism in shaping the complexity of the shoot branching. In addition, we have highlighted that two key components of axillary meristem developmental programs may have been co-opted in controlling flower/ear number of an axillary inflorescence in legumes/maize, hinting at the diversification of axillary-meristem-patterning programs in different lineages. This begs the question how axillary meristem patterning programs may have diversified during plant evolution and hence helped shape the rich variation in shoot branching systems.

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.

Artificial selection of mutations in two nearby genes gave rise to shattering resistance in soybean

Resistance to pod shattering is a key domestication-related trait selected for seed production in many crops. Here, we show that the transition from shattering in wild soybeans to shattering resistance in cultivated soybeans resulted from selection of mutations within the coding sequences of two nearby genes - Sh1 and Pdh1. Sh1 encodes a C2H2-like zinc finger transcription factor that promotes shattering by repressing SHAT1-5 expression, thereby reducing the secondary wall thickness of fiber cap cells in the abscission layers of pod sutures, while Pdh1 encodes a dirigent protein that orchestrates asymmetric lignin distribution in inner sclerenchyma, creating torsion in pod walls that facilitates shattering. Integration analyses of quantitative trait locus mapping, genome-wide association studies, and allele distribution in representative soybean germplasm suggest that these two genes are primary modulators underlying this domestication trait. Our study thus provides comprehensive understanding regarding the genetic, molecular, and cellular bases of shattering resistance in soybeans.

The MADS-box transcription factor GmFULc promotes GmZTL4 gene transcription to modulate maturity in soybean

Flowering time and maturity are crucial agronomic traits that affect the regional adaptability of soybean plants. The development of soybean cultivars with early maturity adapted to longer days and colder climates of high latitudes is very important for ensuring normal ripening before frost begins. FUL belongs to the MADS-box transcription factor family and has several duplicated members in soybeans. In this study, we observed that overexpression of GmFULc in the Dongnong 50 cultivar promoted soybean maturity, while GmFULc knockout mutants exhibited late maturity. Chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq) revealed that GmFULc could bind to the CArG, bHLH and homeobox motifs. Further investigation revealed that GmFULc could directly bind to the CArG motif in the promoters of the GmZTL3 and GmZTL4 genes. Overexpression of GmZTL4 promoted soybean maturity, whereas the ztl4 mutants exhibited delayed maturity. Moreover, we found that the cis element box 4 motif of the GmZTL4 promoter, a motif of light response elements, played an important role in controlling the growth period. Deletion of this motif shortened the growth period by increasing the expression levels of GmZTL4. Functional investigations revealed that short-day treatment promoted the binding of GmFULc to the promoter of GmZTL4 and inhibited the expression of E1 and E1Lb, ultimately resulting in the promotion of flowering and early maturation. Taken together, these findings suggest a novel photoperiod regulatory pathway in which GmFULc directly activates GmZTL4 to promote earlier maturity in soybean.

Molecular selection of soybean towards adaptation to Central European agroclimatic conditions

Europe is highly dependent on soybean meal imports and anticipates an increase of domestic plant protein production. Ongoing climate change resulted in northward shift of plant hardiness zones, enabling spring-sowing of freezing-sensitive crops, including soybean. However, it requires efficient reselection of germplasm adapted to relatively short growing season and long-day photoperiod. In the present study, a PCR array has been implemented, targeting early maturity (E1-E4, E7, E9, and E10), pod shattering (qPHD1), and growth determination (Dt1) genes. This array was optimized for routine screening of soybean diversity panel (204 accessions), subjected to the 2018-2020 survey of phenology, morphology, and yield-related traits in a potential cultivation region in Poland. High broad-sense heritability (0.84-0.88) was observed for plant height, thousand grain weight, maturity date, and the first pod height. Significant positive correlations were identified between the number of seeds and pods per plant, between these two traits and seed yield per plant as well as between flowering, maturity, plant height, and first pod height. PCR array genotyping revealed high genetic diversity, yielding 98 allelic combinations. The most remarkable correlations were identified between flowering and E7 or E1, between maturity and E4 or E7 and between plant height and Dt1 or E4. The study demonstrated high applicability of this PCR array for molecular selection of soybean towards adaptation to Central Europe, designating recessive qPHD1 and dominant Dt1, E3, and E4 alleles as major targets to align soybean growth season requirements with the length of the frost-free period, improve plant performance, and increase yield.

Deciphering PDH1's role in mung bean domestication: a genomic perspective on pod dehiscence

Mung bean (Vigna radiata) stands as a crucial legume crop in Asia, contributing to food security. However, our understanding of the underlying genetic foundation governing domesticated agronomic traits, especially those linked to pod architecture, remains largely unexplored. In this study, we delved into the genomic divergence between wild and domesticated mung bean varieties, leveraging germplasm obtained from diverse sources. Our findings unveiled pronounced variation in promoter regions (35%) between the two mung bean subpopulations, suggesting substantial changes in gene expression patterns during domestication. Leveraging transcriptome analysis using distinct reproductive stage pods and subpopulations, we identified candidate genes responsible for pod and seed architecture development, along with Genome-Wide Association Studies (GWAS) and Quantitative Trait Locus (QTL) analysis. Notably, our research conclusively confirmed PDH1 as a parallel domesticated gene governing pod dehiscence in legumes. This study imparts valuable insights into the genetic underpinnings of domesticated agronomic traits in mung bean, and simultaneously highlighting the parallel domestication of pivotal traits within the realm of legume crops.

CRISPR/Cas9-mediated editing of GmDWF1 brassinosteroid biosynthetic gene induces dwarfism in soybean

KEY MESSAGE

The study on the GmDWF1-deficient mutant dwf1 showed that GmDWF1 plays a crucial role in determining soybean plant height and yield by influencing the biosynthesis of brassinosteroids. Soybean has not adopted the Green Revolution, such as reduced height for increased planting density, which have proven beneficial for cereal crops. Our research identified the soybean genes GmDWF1a and GmDWF1b, homologous to Arabidopsis AtDWF1, and found that they are widely expressed, especially in leaves, and linked to the cellular transport system, predominantly within the endoplasmic reticulum and intracellular vesicles. These genes are essential for the synthesis of brassinosteroids (BR). Single mutants of GmDWF1a and GmDWF1b, as well as double mutants of both genes generated through CRISPR/Cas9 genome editing, exhibit a dwarf phenotype. The single-gene mutant exhibits moderate dwarfism, while the double mutant shows more pronounced dwarfism. Despite the reduced stature, all types of mutants preserve their node count. Notably, field tests have shown that the single GmDWF1a mutant produced significantly more pods than wild-type plants. Spraying exogenous brassinolide (BL) can compensate for the loss in plant height induced by the decrease in endogenous BRs. Comparing transcriptome analyses of the GmDWF1a mutant and wild-type plants revealed a significant impact on the expression of many genes that influence soybean growth. Identifying the GmDWF1a and GmDWF1b genes could aid in the development of compact, densely planted soybean varieties, potentially boosting productivity.

Mechanisms underlying key agronomic traits and implications for molecular breeding in soybean

Soybean (Glycine max [L.] Merr.) is an important crop that provides protein and vegetable oil for human consumption. As soybean is a photoperiod-sensitive crop, its cultivation and yield are limited by the photoperiodic conditions in the field. In contrast to other major crops, soybean has a special plant architecture and a special symbiotic nitrogen fixation system, representing two unique breeding directions. Thus, flowering time, plant architecture, and symbiotic nitrogen fixation are three critical or unique yield-determining factors. This review summarizes the progress made in our understanding of these three critical yield-determining factors in soybean. Meanwhile, we propose potential research directions to increase soybean production, discuss the application of genomics and genomic-assisted breeding, and explore research directions to address future challenges, particularly those posed by global climate changes.

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.

Innovations in functional genomics and molecular breeding of pea: exploring advances and opportunities

The garden pea (Pisum sativum L.) is a significant cool-season legume, serving as crucial food sources, animal feed, and industrial raw materials. The advancement of functional genomics over the past two decades has provided substantial theoretical foundations and progress to pea breeding. Notably, the release of the pea reference genome has enhanced our understanding of plant architecture, symbiotic nitrogen fixation (SNF), flowering time, floral organ development, seed development, and stress resistance. However, a considerable gap remains between pea functional genomics and molecular breeding. This review summarizes the current advancements in pea functional genomics and breeding while highlighting the future challenges in pea molecular breeding.

Dissection of the E8 locus in two early maturing Canadian soybean populations

Soybean [Glycine max (L.) Merr.] is a short-day crop for which breeders want to expand the cultivation range to more northern agro-environments by introgressing alleles involved in early reproductive traits. To do so, we investigated quantitative trait loci (QTL) and expression quantitative trait loci (eQTL) regions comprised within the E8 locus, a large undeciphered region (~7.0 Mbp to 44.5 Mbp) associated with early maturity located on chromosome GM04. We used a combination of two mapping algorithms, (i) inclusive composite interval mapping (ICIM) and (ii) genome-wide composite interval mapping (GCIM), to identify major and minor regions in two soybean populations (QS15524F2:F3 and QS15544RIL) having fixed E1, E2, E3, and E4 alleles. Using this approach, we identified three main QTL regions with high logarithm of the odds (LODs), phenotypic variation explained (PVE), and additive effects for maturity and pod-filling within the E8 region: GM04:16,974,874-17,152,230 (E8-r1); GM04:35,168,111-37,664,017 (E8-r2); and GM04:41,808,599-42,376,237 (E8-r3). Using a five-step variant analysis pipeline, we identified Protein far-red elongated hypocotyl 3 (Glyma.04G124300; E8-r1), E1-like-a (Glyma.04G156400; E8-r2), Light-harvesting chlorophyll-protein complex I subunit A4 (Glyma.04G167900; E8-r3), and Cycling dof factor 3 (Glyma.04G168300; E8-r3) as the most promising candidate genes for these regions. A combinatorial eQTL mapping approach identified significant regulatory interactions for 13 expression traits (e-traits), including Glyma.04G050200 (Early flowering 3/E6 locus), with the E8-r3 region. Four other important QTL regions close to or encompassing major flowering genes were also detected on chromosomes GM07, GM08, and GM16. In GM07:5,256,305-5,404,971, a missense polymorphism was detected in the candidate gene Glyma.07G058200 (Protein suppressor of PHYA-105). These findings demonstrate that the locus known as E8 is regulated by at least three distinct genomic regions, all of which comprise major flowering genes.

Genome-wide association study reveals GmFulb as candidate gene for maturity time and reproductive length in soybeans (Glycine max)

The ability of soybean [Glycine max (L.) Merr.] to adapt to different latitudes is attributed to genetic variation in major E genes and quantitative trait loci (QTLs) determining flowering time (R1), maturity (R8), and reproductive length (RL). Fully revealing the genetic basis of R1, R8, and RL in soybeans is necessary to enhance genetic gains in soybean yield improvement. Here, we performed a genome-wide association analysis (GWA) with 31,689 single nucleotide polymorphisms (SNPs) to detect novel loci for R1, R8, and RL using a soybean panel of 329 accessions with the same genotype for three major E genes (e1-as/E2/E3). The studied accessions were grown in nine environments and observed for R1, R8 and RL in all environments. This study identified two stable peaks on Chr 4, simultaneously controlling R8 and RL. In addition, we identified a third peak on Chr 10 controlling R1. Association peaks overlap with previously reported QTLs for R1, R8, and RL. Considering the alternative alleles, significant SNPs caused RL to be two days shorter, R1 two days later and R8 two days earlier, respectively. We identified association peaks acting independently over R1 and R8, suggesting that trait-specific minor effect loci are also involved in controlling R1 and R8. From the 111 genes highly associated with the three peaks detected in this study, we selected six candidate genes as the most likely cause of R1, R8, and RL variation. High correspondence was observed between a modifying variant SNP at position 04:39294836 in GmFulb and an association peak on Chr 4. Further studies using map-based cloning and fine mapping are necessary to elucidate the role of the candidates we identified for soybean maturity and adaptation to different latitudes and to be effectively used in the marker-assisted breeding of cultivars with optimal yield-related traits.

A multi-trait GWAS-based genetic association network controlling soybean architecture and seed traits

Ideal plant architecture is essential for enhancing crop yields. Ideal soybean (Glycine max) architecture encompasses an appropriate plant height, increased node number, moderate seed weight, and compact architecture with smaller branch angles for growth under high-density planting. However, the functional genes regulating plant architecture are far not fully understood in soybean. In this study, we investigated the genetic basis of 12 agronomic traits in a panel of 496 soybean accessions with a wide geographical distribution in China. Analysis of phenotypic changes in 148 historical elite soybean varieties indicated that seed-related traits have mainly been improved over the past 60 years, with targeting plant architecture traits having the potential to further improve yields in future soybean breeding programs. In a genome-wide association study (GWAS) of 12 traits, we detected 169 significantly associated loci, of which 61 overlapped with previously reported loci and 108 new loci. By integrating the GWAS loci for different traits, we constructed a genetic association network and identified 90 loci that were associated with a single trait and 79 loci with pleiotropic effects. Of these 79 loci, 7 hub-nodes were strongly linked to at least three related agronomic traits. qHub_5, containing the previously characterized Determinate 1 (Dt1) locus, was associated not only with plant height and node number (as determined previously), but also with internode length and pod range. Furthermore, we identified qHub_7, which controls three branch angle-related traits; the candidate genes in this locus may be beneficial for breeding soybean with compact architecture. These findings provide insights into the genetic relationships among 12 important agronomic traits in soybean. In addition, these studies uncover valuable loci for further functional gene studies and will facilitate molecular design breeding of soybean architecture.

k-mer-based GWAS enhances the discovery of causal variants and candidate genes in soybean

Genome-wide association studies (GWAS) are powerful statistical methods that detect associations between genotype and phenotype at genome scale. Despite their power, GWAS frequently fail to pinpoint the causal variant or the gene controlling a given trait in crop species. Assessing genetic variants other than single-nucleotide polymorphisms (SNPs) could alleviate this problem. In this study, we tested the potential of structural variant (SV)- and k-mer-based GWAS in soybean by applying these methods as well as conventional SNP/indel-based GWAS to 13 traits. We assessed the performance of each GWAS approach based on loci for which the causal genes or variants were known from previous genetic studies. We found that k-mer-based GWAS was the most versatile approach and the best at pinpointing causal variants or candidate genes. Moreover, k-mer-based analyses identified promising candidate genes for loci related to pod color, pubescence form, and resistance to Phytophthora sojae. In our dataset, SV-based GWAS did not add value compared to k-mer-based GWAS and may not be worth the time and computational resources invested. Despite promising results, significant challenges remain regarding the downstream analysis of k-mer-based GWAS. Notably, better methods are needed to associate significant k-mers with sequence variation. Our results suggest that coupling k-mer- and SNP/indel-based GWAS is a powerful approach for discovering candidate genes in crop species.

PH13 improves soybean shade traits and enhances yield for high-density planting at high latitudes

Shading in combination with extended photoperiods can cause exaggerated stem elongation (ESE) in soybean, leading to lodging and reduced yields when planted at high-density in high-latitude regions. However, the genetic basis of plant height in adaptation to these regions remains unclear. Here, through a genome-wide association study, we identify a plant height regulating gene on chromosome 13 (PH13) encoding a WD40 protein with three main haplotypes in natural populations. We find that an insertion of a Ty1/Copia-like retrotransposon in the haplotype 3 leads to a truncated PH13H3 with reduced interaction with GmCOP1s, resulting in accumulation of STF1/2, and reduced plant height. In addition, PH13H3 allele has been strongly selected for genetic improvement at high latitudes. Deletion of both PH13 and its paralogue PHP can prevent shade-induced ESE and allow high-density planting. This study provides insights into the mechanism of shade-resistance and offers potential solutions for breeding high-yielding soybean cultivar for high-latitude regions.

Tapping into the plasticity of plant architecture for increased stress resilience

Plant architecture develops post-embryonically and emerges from a dialogue between the developmental signals and environmental cues. Length and branching of the vegetative and reproductive tissues were the focus of improvement of plant performance from the early days of plant breeding. Current breeding priorities are changing, as we need to prioritize plant productivity under increasingly challenging environmental conditions. While it has been widely recognized that plant architecture changes in response to the environment, its contribution to plant productivity in the changing climate remains to be fully explored. This review will summarize prior discoveries of genetic control of plant architecture traits and their effect on plant performance under environmental stress. We review new tools in phenotyping that will guide future discoveries of genes contributing to plant architecture, its plasticity, and its contributions to stress resilience. Subsequently, we provide a perspective into how integrating the study of new species, modern phenotyping techniques, and modeling can lead to discovering new genetic targets underlying the plasticity of plant architecture and stress resilience. Altogether, this review provides a new perspective on the plasticity of plant architecture and how it can be harnessed for increased performance under environmental stress.

QTL Mapping of Soybean (Glycine max) Vine Growth Habit Trait

The vine growth habit (VGH) is a notable property of wild soybean plants that also holds a high degree of importance in domestication as it can preclude using these wild cultivars for breeding and improving domesticated soybeans. Here, a bulked segregant analysis (BSA) approach was employed to study the genetic etiology of the VGH in soybean plants by integrating linkage mapping and population sequencing approaches. To develop a recombinant inbred line (RIL) population, the cultivated Zhongdou41 (ZD41) soybean cultivar was bred with ZYD02787, a wild soybean accession. The VGH status of each line in the resultant population was assessed, ultimately leading to the identification of six and nine QTLs from the BSA sequencing of the F4 population and F6-F8 population sequence mapping, respectively. One QTL shared across these analyzed generations was detected on chromosome 19. Three other QTLs detected by BSA-seq were validated and localized to the 90.93 kb, 2.9 Mb, and 602.08 kb regions of chromosomes 6 and 13, harboring 14, 53, and 4 genes, respectively. Three consistent VGH-related QTLs located on chromosomes 2 and 19 were detected in a minimum of three environments, while an additional six loci on chromosomes 2, 10, 13, and 18 were detected in at least two environments via ICIM mapping. Of all the detected loci, five had been reported previously whereas seven represent novel QTLs. Together, these data offer new insights into the genetic basis of the VGH in soybean plants, providing a rational basis to inform the use of wild accessions in future breeding efforts.

Ideotype breeding and genome engineering for legume crop improvement

Ideotype breeding is a strategy whereby traits are modeled a priori and then introduced into a model or crop species to assess their impact on yield. Thus, knowledge about the connection between genotype and phenotype is required for ideotype breeding to be deployed successfully. The growing understanding of the genetic basis of yield-related traits, combined with increasingly efficient genome engineering tools, improved transformation efficiency, and high-throughput genotyping of regenerants paves the way for the widespread adoption of ideotype breeding as a complement to conventional breeding. We briefly discuss how ideotype breeding, coupled with such state-of-the-art biotechnological tools, could contribute to knowledge-based legume breeding and accelerate yield gains to ensure food security in the coming decades.

The genetic basis of shoot architecture in soybean

Shoot architecture refers to the three-dimensional body plan of the above ground organs of the plant. The patterning of this body plan results from the tight genetic control of the size and maintenance of meristems, the initiation of axillary growth, and the timing of developmental phase transition. Variation in shoot architecture can result in dramatic differences in plant productivity and/or grain yield due to their effects on light interception, photosynthetic efficiency, response to agronomic inputs, and environmental adaptation. The fine-tuning of shoot architecture has consequently been of great interest to plant breeders, driving the need for deeper understanding of the genes and molecular mechanisms governing these traits. In soybean, the world's most important oil and protein crop, major components of shoot architecture include stem growth habit, plant height, branch angle, branch number, leaf petiole angle, and the size and shape of leaves. Key genes underlying some of these traits have been identified to integrate hormonal, developmental, and environmental signals modulating the growth and orientation of shoot organs. Here we summarize the current knowledge and recent advances in the understanding of the genetic control of these important architectural traits in soybean.

A genome-wide analysis of the USDA Soybean Isoline Collection

The USDA Soybean Isoline Collection has been an invaluable resource for the soybean genetics and breeding community. This collection, established in 1972, consists of 611 near-isogenic lines (NILs) carrying one or multiple genes conferring traits that had been determined to exhibit Mendelian inheritance. It has been used in multiple studies on the genetic basis, physiology, and agronomy of these qualitative traits. Here, we used publicly available genotype (SoySNP50K), phenotype, and pedigree data on this collection to characterize the isogenicity of the NILs and identify chromosomal positions of unmapped genes. A total of 368 NILs had at least 80% identity to their recurrent parent and, thus, were useful for what can be called introgression mapping. Both on-target and off-target introgressions were evaluated. The size of on-target introgressions into individual NILs ranged from 61 kb to 8.4 Mb, whereas off-target introgressions ranged from 2.6 kb to 54.8 Mb. The observed large off-target introgressions indicated that some NILs carry introgressions nearly the size of an entire chromosome. By applying introgression mapping to genes that had never been mapped, we identified the likely chromosomal positions of six such genes: ab, im, lo, Np, pc, and Rpm. The size of mapping intervals was large in some cases (10.28 Mb for im) but small in others (0.21 Mb for Np). The results reported herein will provide future researchers with a resource to help select informative NILs for future studies, and provide a starting point to further fine map, and ultimately clone and functionally characterize these six soybean genes.

Restructuring plant types for developing tailor-made crops

Plants have adapted to different environmental niches by fine-tuning the developmental factors working together to regulate traits. Variations in the developmental factors result in a wide range of quantitative variations in these traits that helped plants survive better. The major developmental pathways affecting plant architecture are also under the control of such pathways. Most notable are the CLAVATA-WUSCHEL pathway regulating shoot apical meristem fate, GID1-DELLA module influencing plant height and tillering, LAZY1-TAC1 module controlling branch/tiller angle and the TFL1-FT determining the floral fate in plants. Allelic variants of these key regulators selected during domestication shaped the crops the way we know them today. There is immense yield potential in the 'ideal plant architecture' of a crop. With the available genome-editing techniques, possibilities are not restricted to naturally occurring variations. Using a transient reprogramming system, one can screen the effect of several developmental gene expressions in novel ecosystems to identify the best targets. We can use the plant's fine-tuning mechanism for customizing crops to specific environments. The process of crop domestication can be accelerated with a proper understanding of these developmental pathways. It is time to step forward towards the next-generation molecular breeding for restructuring plant types in crops ensuring yield stability.

Going north: adaptation of soybean to long-day environments

This article comments on:

Zhu X, Leiser WL, Hahn V, Würschum T. 2023. The genetic architecture of soybean photothermal adaptation to high latitudes. Journal of Experimental Botany 74,2987–3002

In plant breeding, understanding genetic variation in the photoperiodic control of flowering time of crop plants such as soybean is a prerequisite for managing adaptation to new environments. Zhu et al. (2023) analyzed a large diversity panel of >1500 early maturity soybean lines to disclose the genetic architecture behind the timing of flowering and maturity. Their findings confirm known maturity loci and reveal new candidate genes and alleles as well as environmental interactions of individual quantitative trait loci (QTLs) for flowering and maturity time. The results shed light on the complexity of the regulatory network which controls the timing of flowering in soybean. This supports the fine-tuning of plant architectures through the combination of stem termination and flowering genes towards a better adaptation of soybean to high latitudes or other stressful environments.

The genetic architecture of soybean photothermal adaptation to high latitudes

Soybean is a major plant protein source for both human food and animal feed, but to meet global demands as well as a trend towards regional production, soybean cultivation needs to be expanded to higher latitudes. In this study, we developed a large diversity panel consisting of 1503 early-maturing soybean lines and used genome-wide association mapping to dissect the genetic architecture underlying two crucial adaptation traits, flowering time and maturity. This revealed several known maturity loci, E1, E2, E3, and E4, and the growth habit locus Dt2 as causal candidate loci, and also a novel putative causal locus, GmFRL1, encoding a homolog of the vernalization pathway gene FRIGIDA-like 1. In addition, the scan for quantitative trait locus (QTL)-by-environment interactions identified GmAPETALA1d as a candidate gene for a QTL with environment-dependent reversed allelic effects. The polymorphisms of these candidate genes were identified using whole-genome resequencing data of 338 soybeans, which also revealed a novel E4 variant, e4-par, carried by 11 lines, with nine of them originating from Central Europe. Collectively, our results illustrate how combinations of QTL and their interactions with the environment facilitate the photothermal adaptation of soybean to regions far beyond its center of origin.

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.

Genome-wide association study for biomass accumulation traits in soybean

Soybean is one of the most versatile crops for oil production, human diets, and feedstocks. The vegetative biomass of soybean is an important determinant of seed yield and is crucial for the forage usages. However, the genetic control of soybean biomass is not well explained. In this work, we used a soybean germplasm population, including 231 improved cultivars, 207 landraces, and 121 wild soybeans, to investigate the genetic basis of biomass accumulation of soybean plants at the V6 stage. We found that biomass-related traits, including NDW (nodule dry weight), RDW (root dry weight), SDW (shoot dry weight), and TDW (total dry weight), were domesticated during soybean evolution. In total, 10 loci, encompassing 47 putative candidate genes, were detected for all biomass-related traits by a genome-wide association study. Among these loci, seven domestication sweeps and six improvement sweeps were identified. Glyma.05G047900, a purple acid phosphatase, was a strong candidate gene to improve biomass for future soybean breeding. This study provided new insights into the genetic basis of biomass accumulation during soybean evolution.

Supplementary information

The online version contains supplementary material available at 10.1007/s11032-023-01380-6.

Genome-wide survey identified superior and rare haplotypes for plant height in the north-eastern soybean germplasm of China

The proper and efficient utilization of natural genetic diversity can significantly impact crop improvements. Plant height is a quantitative trait governing the plant type as well as the yield and quality of soybean. Here, we used a combined approach including a genome-wide association study (GWAS) and haplotype and candidate gene analyses to explore the genetic basis of plant height in diverse natural soybean populations. For the GWAS analysis, we used the whole-genome resequencing data of 196 diverse soybean cultivars collected from different accumulated temperature zones of north-eastern China to detect the significant single-nucleotide polymorphisms (SNPs) associated with plant height across three environments (E1, E2, and E3). A total of 33 SNPs distributed on four chromosomes, viz., Chr.02, Chr.04, Chr.06, and Chr.19, were identified to be significantly associated with plant height across the three environments. Among them, 23 were consistently detected in two or more environments and the remaining 10 were identified in only one environment. Interestingly, all the significant SNPs detected on the respective chromosomes fell within the physical interval of linkage disequilibrium (LD) decay (± 38.9 kb). Hence, these genomic regions were considered to be four quantitative trait loci (QTLs), viz., qPH2, qPH4, qPH6, and qPH19, regulating plant height. Moreover, the genomic region flanking all significant SNPs on four chromosomes exhibited strong LD. These significant SNPs thus formed four haplotype blocks, viz., Hap-2, Hap-4, Hap-6, and Hap-19. The number of haplotype alleles underlying each block varied from four to six, and these alleles regulate the different phenotypes of plant height ranging from dwarf to extra-tall heights. Nine candidate genes were identified within the four haplotype blocks, and these genes were considered putative candidates regulating soybean plant height. Hence, these stable QTLs, superior haplotypes, and candidate genes (after proper validation) can be deployed for the development of soybean cultivars with desirable plant heights.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01363-7.

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.

A multi-locus linear mixed model methodology for detecting small-effect QTLs for quantitative traits in MAGIC, NAM, and ROAM populations

Although multi-parent populations (MPPs) integrate the advantages of linkage and association mapping populations in the genetic dissection of complex traits and especially combine genetic analysis with crop breeding, it is difficult to detect small-effect quantitative trait loci (QTL) for complex traits in multiparent advanced generation intercross (MAGIC), nested association mapping (NAM), and random-open-parent association mapping (ROAM) populations. To address this issue, here we proposed a multi-locus linear mixed model method, namely mppQTL, to detect QTLs, especially small-effect QTLs, in these MPPs. The new method includes two steps. The first is genome-wide scanning based on a single-locus linear mixed model; the P-values are obtained from likelihood-ratio test, the peaks of negative logarithm P-value curve are selected by group-lasso, and all the selected peaks are regarded as potential QTLs. In the second step, all the potential QTLs are placed on a multi-locus linear mixed model, all the effects are estimated using expectation-maximization empirical Bayes algorithm, and all the non-zero effect vectors are further evaluated via likelihood-ratio test for significant QTLs. In Monte Carlo simulation studies, the new method has higher power in QTL detection, lower false positive rate, lower mean absolute deviation for QTL position estimate, and lower mean squared error for the estimate of QTL size (r²) than existing methods because the new method increases the power of detecting small-effect QTLs. In real dataset analysis, the new method (19) identified five more known genes than the existing three methods (14). This study provides an effective method for detecting small-effect QTLs in any MPPs.

The soybean ubiquitin-proteasome system: Current knowledge and future perspective

Increasing soybean [Glycine max (L.) Merr.] yield has become a worldwide scientific problem in the world. Many studies have shown that ubiquitination plays a key role in stress response and yield formation. In the UniProtKB database, 2,429 ubiquitin-related proteins were predicted in soybean, however, <20 were studied. One key way to address this lack of progress in increasing soybean yield will be a deeper understanding of the ubiquitin-proteasome system (UPS) in soybean. In this review, we summarized the current knowledge about soybean ubiquitin-related proteins and discussed the method of combining phenotype, mutant library, transgenic system, genomics, and proteomics approaches to facilitate the exploration of the soybean UPS. We also proposed the strategy of applying the UPS in soybean improvement based on related studies in model plants. Our review will be helpful for soybean scientists to learn current research progress of the soybean UPS and further lay a theoretical reference for the molecular improvement of soybean in future research by use of this knowledge.

Branch angle and leaflet shape are associated with canopy coverage in soybean

Early canopy coverage is a desirable trait that is a major determinant of yield in soybean (Glycine max). Variation in traits comprising shoot architecture can influence canopy coverage, canopy light interception, canopy-level photosynthesis, and source-sink partitioning efficiency. However, little is known about the extent of phenotypic diversity of shoot architecture traits and their genetic control in soybean. Thus, we sought to understand the contribution of shoot architecture traits to canopy coverage and to determine the genetic control of these traits. We examined the natural variation for shoot architecture traits in a set of 399 diverse maturity group I soybean (SoyMGI) accessions to identify relationships between traits, and to identify loci that are associated with canopy coverage and shoot architecture traits. Canopy coverage was correlated with branch angle, number of branches, plant height, and leaf shape. Using previously collected 50K single nucleotide polymorphism data, we identified quantitative trait locus (QTL) associated with branch angle, number of branches, branch density, leaflet shape, days to flowering, maturity, plant height, number of nodes, and stem termination. In many cases, QTL intervals overlapped with previously described genes or QTL. We also found QTL associated with branch angle and leaflet shape located on chromosomes 19 and 4, respectively, and these QTL overlapped with QTL associated with canopy coverage, suggesting the importance of branch angle and leaflet shape in determining canopy coverage. Our results highlight the role individual architecture traits play in canopy coverage and contribute information on their genetic control that could help facilitate future efforts in their genetic manipulation.

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.

Improvement of Soybean; A Way Forward Transition from Genetic Engineering to New Plant Breeding Technologies

Soybean is considered one of the important crops among legumes. Due to high nutritional contents in seed (proteins, sugars, oil, fatty acids, and amino acids), soybean is used globally for food, feed, and fuel. The primary consumption of soybean is vegetable oil and feed for chickens and livestock. Apart from this, soybean benefits soil fertility by fixing atmospheric nitrogen through root nodular bacteria. While conventional breeding is practiced for soybean improvement, with the advent of new biotechnological methods scientists have also engineered soybean to improve different traits (herbicide, insect, and disease resistance) to fulfill consumer requirements and to meet the global food deficiency. Genetic engineering (GE) techniques such as transgenesis and gene silencing help to minimize the risks and increase the adaptability of soybean. Recently, new plant breeding technologies (NPBTs) emerged such as zinc-finger nucleases, transcription activator-like effector nucleases, and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR/Cas9), which paved the way for enhanced genetic modification of soybean. These NPBTs have the potential to improve soybean via gene functional characterization precision genome engineering for trait improvement. Importantly, these NPBTs address the ethical and public acceptance issues related to genetic modifications and transgenesis in soybean. In the present review, we summarized the improvement of soybean through GE and NPBTs. The valuable traits that have been improved through GE for different constraints have been discussed. Moreover, the traits that have been improved through NPBTs and potential targets for soybean improvements via NPBTs and solutions for ethical and public acceptance are also presented.

Candidate loci for breeding compact plant-type soybean varieties

Plant height and node number are important agronomic traits that influence yield in soybean (Glycine max L.). Here, to better understand the genetic basis of the traits, we used two recombinant inbred line (RIL) populations to detect quantitative trait loci (QTLs) associated with plant height and node number in different environments. This analysis detected 9 and 21 QTLs that control plant height and node number, respectively. Among them, we identified two genomic regions that overlap with Determinate stem 1 (Dt1) and Dt2, which are known to influence both plant height and node number. Furthermore, different combinations of Dt1 and Dt2 alleles were enriched in distinct latitudes. In addition, we determined that the QTLs qPH-13-SE and qPH-13-DW in the two RIL populations overlap with genomic intervals associated with plant height and the QTL qNN-04-DW overlaps with an interval associated with node number. Combining the dwarf allele of qPH-13-SE/qPH-13-DW and the multiple-node allele of qNN-04-DW produced plants with ideal plant architecture, i.e., shorter main stems with more nodes. This plant type may help increase yield at high planting density. This study thus provides candidate loci for breeding elite soybean cultivars for plant height and node number.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01352-2.

CRISPR/Cas9-mediated targeted mutation of the E1 decreases photoperiod sensitivity, alters stem growth habits, and decreases branch number in soybean

The distribution of elite soybean (Glycine max) cultivars is limited due to their highly sensitive to photoperiod, which affects the flowering time and plant architecture. The recent emergence of CRISPR/Cas9 technology has uncovered new opportunities for genetic manipulation of soybean. The major maturity gene E1 of soybean plays a critical role in soybean photoperiod response. Here, we performed CRISPR/Cas9-mediated targeted mutation of E1 gene in soybean cultivar Tianlong1 carrying the dominant E1 to investigate its precise function in photoperiod regulation, especially in plant architecture regulation. Four types of mutations in the E1 coding region were generated. No off-target effects were observed, and homozygous trans-clean mutants without T-DNA were obtained. The photoperiod sensitivity of e1 mutants decreased relative to the wild type plants; however, e1 mutants still responded to photoperiod. Further analysis revealed that the homologs of E1, E1-La, and E1-Lb, were up-regulated in the e1 mutants, indicating a genetic compensation response of E1 and its homologs. The e1 mutants exhibited significant changes in the architecture, including initiation of terminal flowering, formation of determinate stems, and decreased branch numbers. To identify E1-regulated genes related to plant architecture, transcriptome deep sequencing (RNA-seq) was used to compare the gene expression profiles in the stem tip of the wild-type soybean cultivar and the e1 mutants. The expression of shoot identity gene Dt1 was significantly decreased, while Dt2 was significantly upregulated. Also, a set of MADS-box genes was up-regulated in the stem tip of e1 mutants which might contribute to the determinate stem growth habit.

Prospects of Feral Crop De Novo Redomestication

Modern agriculture depends on a narrow variety of crop species, leaving global food and nutritional security highly vulnerable to the adverse effects of climate change and population expansion. Crop improvement using conventional and molecular breeding approaches leveraging plant genetic diversity using crop wild relatives (CWRs) has been one approach to address these issues. However, the rapid pace of the global change requires additional innovative solutions to adapt agriculture to meet global needs. Neodomestication-the rapid and targeted introduction of domestication traits using introgression or genome editing of CWRs-is being explored as a supplementary approach. These methods show promise; however, they have so far been limited in efficiency and applicability. We propose expanding the scope of neodomestication beyond truly wild CWRs to include feral crops as a source of genetic diversity for novel crop development, in this case 'redomestication'. Feral crops are plants that have escaped cultivation and evolved independently, typically adapting to their local environments. Thus, feral crops potentially contain valuable adaptive features while retaining some domestication traits. Due to their genetic proximity to crop species, feral crops may be easier targets for de novo domestication (i.e. neodomestication via genome editing techniques). In this review, we explore the potential of de novo redomestication as an application for novel crop development by genome editing of feral crops. This approach to efficiently exploit plant genetic diversity would access an underutilized reservoir of genetic diversity that could prove important in support of global food insecurity in the face of the climate change.

The Key to the Future Lies in the Past: Insights from Grain Legume Domestication and Improvement Should Inform Future Breeding Strategies

Crop domestication is a co-evolutionary process that has rendered plants and animals significantly dependent on human interventions for survival and propagation. Grain legumes have played an important role in the development of Neolithic agriculture some 12,000 years ago. Despite being early companions of cereals in the origin and evolution of agriculture, the understanding of grain legume domestication has lagged behind that of cereals. Adapting plants for human use has resulted in distinct morpho-physiological changes between the wild ancestors and domesticates, and this distinction has been the focus of several studies aimed at understanding the domestication process and the genetic diversity bottlenecks created. Growing evidence from research on archeological remains, combined with genetic analysis and the geographical distribution of wild forms, has improved the resolution of the process of domestication, diversification and crop improvement. In this review, we summarize the significance of legume wild relatives as reservoirs of novel genetic variation for crop breeding programs. We describe key legume features, which evolved in response to anthropogenic activities. Here, we highlight how whole genome sequencing and incorporation of omics-level data have expanded our capacity to monitor the genetic changes accompanying these processes. Finally, we present our perspective on alternative routes centered on de novo domestication and re-domestication to impart significant agronomic advances of novel crops over existing commodities. A finely resolved domestication history of grain legumes will uncover future breeding targets to develop modern cultivars enriched with alleles that improve yield, quality and stress tolerance.

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.

Natural variation of Dt2 determines branching in soybean

Shoot branching is fundamentally important in determining soybean yield. Here, through genome-wide association study, we identify one predominant association locus on chromosome 18 that confers soybean branch number in the natural population. Further analyses determine that Dt2 is the corresponding gene and the natural variations in Dt2 result in significant differential transcriptional levels between the two major haplotypes. Functional characterization reveals that Dt2 interacts with GmAgl22 and GmSoc1a to physically bind to the promoters of GmAp1a and GmAp1d and to activate their transcription. Population genetic investigation show that the genetic differentiation of Dt2 display significant geographic structure. Our study provides a predominant gene for soybean branch number and may facilitate the breeding of high-yield soybean varieties.

Identification and characterization of putative targets of VEGETATIVE1/FULc, a key regulator of development of the compound inflorescence in pea and related legumes

Inflorescence architecture contributes to essential plant traits. It determines plant shape, contributing to morphological diversity, and also determines the position and number of flowers and fruits produced by the plant, thus influencing seed yield. Most legumes have compound inflorescences, where flowers are produced in secondary inflorescences (I2), formed at the flanks of the main primary inflorescence (I1), in contrast to simple inflorescences of plants like Arabidopsis, in which flowers are directly formed on the I1. The pea VEGETATIVE1/FULc (VEG1) gene, and its homologs in other legumes, specify the formation of the I2 meristem, a function apparently restricted to legumes. To understand the control of I2 development, it is important to identify the genes working downstream of VEG1. In this study, we adopted a novel strategy to identify genes expressed in the I2 meristem, as potential regulatory targets of VEG1. To identify pea I2-meristem genes, we compared the transcriptomes of inflorescence apices from wild-type and mutants affected in I2 development, such as proliferating inflorescence meristems (pim, with more I2 meristems), and veg1 and vegetative2 (both without I2 meristems). Analysis of the differentially expressed genes using Arabidopsis genome databases combined with RT-qPCR expression analysis in pea allowed the selection of genes expressed in the pea inflorescence apex. In situ hybridization of four of these genes showed that all four genes are expressed in the I2 meristem, proving our approach to identify I2-meristem genes was successful. Finally, analysis by VIGS (virus-induced gene silencing) in pea identified one gene, PsDAO1, whose silencing leads to small plants, and another gene, PsHUP54, whose silencing leads to plants with very large stubs, meaning that this gene controls the activity of the I2 meristem. PsHUP54-VIGS plants are also large and, more importantly, produce large pods with almost double the seeds as the control. Our study shows a new useful strategy to isolate I2-meristem genes and identifies a novel gene, PsHUP54, which seems to be a promising tool to improve yield in pea and in other legumes.

An APETALA2/ethylene responsive factor transcription factor GmCRF4a regulates plant height and auxin biosynthesis in soybean

Plant height is one of the key agronomic traits affecting soybean yield. The cytokinin response factors (CRFs), as a branch of the APETALA2/ethylene responsive factor (AP2/ERF) super gene family, have been reported to play important roles in regulating plant growth and development. However, their functions in soybean remain unknown. This study characterized a soybean CRF gene named GmCRF4a by comparing the performance of the homozygous Gmcrf4a-1 mutant, GmCRF4a overexpression (OX) and co-silencing (CS) lines. Phenotypic analysis showed that overexpression of GmCRF4a resulted in taller hypocotyls and epicotyls, more main stem nodes, and higher plant height. While down-regulation of GmCRF4a conferred shorter hypocotyls and epicotyls, as well as a reduction in plant height. The histological analysis results demonstrated that GmCRF4a promotes epicotyl elongation primarily by increasing cell length. Furthermore, GmCRF4a is required for the expression of GmYUCs genes to elevate endogenous auxin levels, which may subsequently enhance stem elongation. Taken together, these observations describe a novel regulatory mechanism in soybean, and provide the basis for elucidating the function of GmCRF4a in auxin biosynthesis pathway and plant heigh regulation in plants.

GmIAA27 Encodes an AUX/IAA Protein Involved in Dwarfing and Multi-Branching in Soybean

Soybean plant height and branching affect plant architecture and yield potential in soybean. In this study, the mutant dmbn was obtained by treating the cultivar Zhongpin 661 with ethylmethane sulfonate. The dmbn mutant plants were shorter and more branched than the wild type. The genetic analysis showed that the mutant trait was controlled by a semi-dominant gene. The candidate gene was fine-mapped to a 91 kb interval on Chromosome 9 by combining BSA-seq and linkage analysis. Sequence analysis revealed that Glyma.09g193000 encoding an Aux/IAA protein (GmIAA27) was mutated from C to T in the second exon of the coding region, resulting to amino acid substitution of proline to leucine. Overexpression of the mutant type of this gene in Arabidopsis thaliana inhibited apical dominance and promoted lateral branch development. Expression analysis of GmIAA27 and auxin response genes revealed that some GH3 genes were induced. GmIAA27 relies on auxin to interact with TIR1, whereas Gmiaa27 cannot interact with TIR1 owing to the mutation in the degron motif. Identification of this unique gene that controls soybean plant height and branch development provides a basis for investigating the mechanisms regulating soybean plant architecture development.

The GmSNAP11 Contributes to Resistance to Soybean Cyst Nematode Race 4 in Glycine max

Soybean cyst nematode (SCN) has devastating effects on soybean production, making it crucial to identify genes conferring SCN resistance. Here we employed next-generation sequencing-based bulked segregant analysis (BSA) to discover genomic regions, candidate genes, and diagnostic markers for resistance to SCN race 4 (SCN4) in soybean. Phenotypic analysis revealed highly significant differences among the reactions of 145 recombinant inbred lines (RILs) to SCN4. In combination with euclidean distance (ED) and Δsingle-nucleotide polymorphism (SNP)-index analyses, we identified a genomic region on Gm11 (designated as rhg1-paralog) associated with SCN4 resistance. Overexpression and RNA interference analyzes of the two candidate genes identified in this region (GmPLAC8 and GmSNAP11) revealed that only GmSNAP11 significantly contributes to SCN4 resistance. We developed a diagnostic marker for GmSNAP11. Using this marker, together with previously developed markers for SCN-resistant loci, rhg1 and Rhg4, we evaluated the relationship between genotypes and SCN4 resistance in 145 RILs and 30 soybean accessions. The results showed that all the SCN4-resistant lines harbored all the three loci, however, some lines harboring the three loci were still susceptible to SCN4. This suggests that these three loci are necessary for the resistance to SCN4, but they alone cannot confer full resistance. The GmSNAP11 and the diagnostic markers developed could be used in genomic-assisted breeding to develop soybean varieties with increased resistance to SCN4.

Identification of superior haplotypes in a diverse natural population for breeding desirable plant height in soybean

KEY MESSAGE

Plant height of soybean is associated with a haplotype block on chromosome 19, which classified 211 soybean accessions into five distinct groups showing significant differences for the target trait. Genetic variation is pivotal for crop improvement. Natural populations are precious genetic resources. However, efficient strategies for the targeted utilization of these resources for quantitative traits, such as plant height (PH), are scarce. Being an important agronomic trait associated with soybean yield and quality, it is imperative to unravel the genetic mechanisms underlying PH in soybean. Here, a genome-wide association study (GWAS) was performed to identify single nucleotide polymorphisms (SNPs) significantly associated with PH in a natural population of 211 cultivated soybeans, which was genotyped with NJAU 355 K Soy SNP Array and evaluated across six environments. A total of 128 SNPs distributed across 17 chromosomes were found to be significantly associated with PH across six environments and a combined environment. Three significant SNPs were consistently identified in at least three environments on Chr.02 (AX-93958260), Chr.17 (AX-94154834), and Chr.19 (AX-93897200). Genomic regions of ~ 130 kb flanking these three consistent SNPs were considered as stable QTLs, which included 169 genes. Of these, 22 genes (including Dt1) were prioritized and defined as putative candidates controlling PH. The genomic region flanking 12 most significant SNPs was in strong linkage disequilibrium (LD). These SNPs formed a single haplotype block containing five haplotypes for PH, namely Hap-A, Hap-B, Hap-C, Hap-D, and Hap-E. Deployment of such superior haplotypes in breeding programs will enable development of improved soybean varieties with desirable plant height.

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.

The control of compound inflorescences: insights from grasses and legumes

A major challenge in biology is to understand how organisms have increased developmental complexity during evolution. Inflorescences, with remarkable variation in branching systems, are a fitting model to understand architectural complexity. Inflorescences bear flowers that may become fruits and/or seeds, impacting crop productivity and species fitness. Great advances have been achieved in understanding the regulation of complex inflorescences, particularly in economically and ecologically important grasses and legumes. Surprisingly, a synthesis is still lacking regarding the common or distinct principles underlying the regulation of inflorescence complexity. Here, we synthesize the similarities and differences in the regulation of compound inflorescences in grasses and legumes, and propose that the emergence of novel higher-order repetitive modules is key to the evolution of inflorescence complexity.

Utilization of Plant Architecture Genes in Soybean to Positively Impact Adaptation to High Yield Environments

Optimization of plant architecture by modifying stem termination and timing of flowering and maturity of soybean is a promising strategy to improve its adaptability to specific production environments. Therefore, it is important to choose a proper stem termination type and to understand morphological differences between each stem termination type under various environmental conditions. Variations in abruptness of stem termination have been generally classified into three classical genetic types, indeterminate (Dt1), determinate (dt1), and semi-determinate (Dt2). However, an additional stem termination type, termed tall determinate, and its genetic symbol, dt1-t, were introduced about 25 years ago. The tall determinate soybean lines show delayed cessation of apical stem growth and about 50% taller plant heights than the typical determinate soybeans, even though the genetic control of the tall determinate phenotype was found to be allelic to dt1. Despite the potential agronomic merits of the alternative stem termination type, knowledge about the tall determinate soybean remains limited. We clarified the molecular basis of the tall determinate stem termination type and examined potential agronomic merits of the alternative stem type under three different production environments in the US. Sequence analysis of the classical tall determinate soybean lines revealed that the dt1-t allele responsible for tall determinate stem architecture is caused by two of the identified independent missense alleles of dt1, dt1-t1 (R130K), and dt1-t2 (R62S). Also, from the comparison among soybean accessions belonging to each of the genotype categories for stem termination types, soybean accessions with tall determinate alleles were found to have a high discrepancy rate in phenotyping. Newly developed tall determinate late-maturing soybean germplasm lines had taller plant heights and a greater number of nodes with a similar stem diameter and similar pod density at the apical stem compared to typical determinate soybeans having dt1 (R166W) alleles in Southern environments in the US. The phenotype of increased pod-bearing nodes with lodging resistance has the potential to improve yield, especially grown in high yield environments. This study suggests an alternative strategy to remodel the shape of soybean plants, which can possibly lead to yield improvement through the modification of soybean plant architecture.

Induced Mutation in GmCOP1b Enhances the Performance of Soybean under Dense Planting Conditions

CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1) is the key photomorphogenic inhibitor that has been extensively studied in higher plants. Nevertheless, its role has not been documented in the economically important soybean. Here we investigated the functions of two COP1 homologous genes, GmCOP1a and GmCOP1b, by analyzing Gmcop1a and Gmcop1b mutants with indels using CRISPR in soybean. We revealed that, although both genes are required for skotomorphogenesis in the dark, the GmCOP1b gene seems to play a more prominent role than GmCOP1a in promoting stem elongation under normal light conditions. Consistently, the bZIP transcriptional factors STF1/2, which repress stem elongation in soybean, accumulated to the highest level in the Gmcop1a1b double mutant, followed by the Gmcop1b and Gmcop1a mutants. Furthermore, the Gmcop1b mutants showed reduced shade response and enhanced performance under high-density conditions in field trials. Taken together, this study provides essential genetic resources for elucidating functional mechanisms of GmCOP1 and breeding of high yield soybean cultivars for future sustainable agriculture.

The miR166 mediated regulatory module controls plant height by regulating gibberellic acid biosynthesis and catabolism in soybean

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that play critical roles in regulating plant growth and development. Here, we used Short Tandem Target Mimic (STTM) technology to generate soybean (Glycine max (L.) Merr.) miRNA knockdown lines and identify miRNAs that regulate plant height, a key agronomic trait that affects yield. STTM166 successfully silenced miR166 in soybean and upregulated the expression of miR166 target genes, such as ATHB14-LIKE. The miR166 knockdown lines (GmSTTM166) displayed a reduced plant height phenotype. Moreover, GmSTTM166 plants contained lower levels of bioactive gibberellic acid (GA3) than wild-type plants, and application of exogenous GA partially rescued the dwarf phenotype of GmSTTM166. Knockdown of miR166 altered the expression of genes involved in GA biosynthesis and catabolism. Further analysis revealed that ATHB14-LIKE directly represses transcription of the GA biosynthesis genes GmGA1 and GmGA2, while activating transcription of the GA catabolic gene GIBBERLLIN 2 OXIDASE 2 (GmGA2ox2). Collectively, these results reveal a pivotal role for miR166 in the genetic control of plant height in soybean, thereby providing invaluable insights for molecular breeding to improve soybean yield.

A functionally divergent SOC1 homolog improves soybean yield and latitudinal adaptation

Soybean (Glycine max) grows in a wide range of latitudes, but it is extremely sensitive to photoperiod, which reduces its yield and ability to adapt to different environments. Therefore, understanding of the genetic basis of soybean adaptation is of great significance for breeding and improvement. Here, we characterized Tof18 (SOC1a) that conditions early flowering and growth habit under both short-day and long-day conditions. Molecular analysis confirmed that the two SOC1 homologs present in soybeans (SOC1a and SOC1b) underwent evolutionary functional divergence, with SOC1a having stronger effects on flowering time and stem node number than SOC1b due to transcriptional differences. soc1a soc1b double mutants showed stronger functional effects than either of the single mutants, perhaps due to the formation of SOC1a and SOC1b homodimers or heterodimers. Additionally, Tof18/SOC1a improves the latitudinal adaptation of cultivated soybeans, highlighting the functional importance of SOC1a. The Tof18G allele facilitates adaptation to high latitudes, whereas Tof18A facilitates adaptation to low latitudes. We demonstrated that SOC1s contribute to floral induction in both leaves and shoot apex through inter-regulation with FTs. The SOC1a-SOC1b-Dt2 complex plays essential roles in stem growth habit by directly binding to the regulatory sequence of Dt1, making the genes encoding these proteins potential targets for genome editing to improve soybean yield via molecular breeding. Since the natural Tof18A allele increases node number, introgressing this allele into modern cultivars could improve yields, which would help optimize land use for food production in the face of population growth and global warming.

A novel Synthetic phenotype association study approach reveals the landscape of association for genomic variants and phenotypes

INTRODUCTION

Genome-Wide Association Studies (GWAS) identify tagging variants in the genome that are statistically associated with the phenotype because of their linkage disequilibrium (LD) relationship with the causative mutation (CM). When both low-density genotyped accession panels with phenotypes and resequenced data accession panels are available, tagging variants can assist with post-GWAS challenges in CM discovery.

OBJECTIVES

Our objective was to identify additional GWAS evaluation criteria to assess correspondence between genomic variants and phenotypes, as well as enable deeper analysis of the localized landscape of association.

METHODS

We used genomic variant positions as Synthetic phenotypes in GWAS that we named "Synthetic phenotype association study" (SPAS). The extreme case of SPAS is what we call an "Inverse GWAS" where we used CM positions of cloned soybean genes. We developed and validated the Accuracy concept as a measure of the correspondence between variant positions and phenotypes.

RESULTS

The SPAS approach demonstrated that the genotype status of an associated variant used as a Synthetic phenotype enabled us to explore the relationships between tagging variants and CMs, and further, that utilizing CMs as Synthetic phenotypes in Inverse GWAS illuminated the landscape of association. We implemented the Accuracy calculation for a curated accession panel to an online Accuracy calculation tool (AccuTool) as a resource for gene identification in soybean. We demonstrated our concepts on three examples of soybean cloned genes. As a result of our findings, we devised an enhanced "GWAS to Genes" analysis (Synthetic phenotype to CM strategy, SP2CM). Using SP2CM, we identified a CM for a novel gene.

CONCLUSION

The SP2CM strategy utilizing Synthetic phenotypes and the Accuracy calculation of correspondence provides crucial information to assist researchers in CM discovery. The impact of this work is a more effective evaluation of landscapes of GWAS associations.