The soybean stem growth habit gene Dt1 is an ortholog of Arabidopsis TERMINAL FLOWER1

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.

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.

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.

MsTFL1A delays flowering and regulates shoot architecture and root development in Medicago sativa

KEY MESSAGE

MsTFL1A is an important gene involved in flowering repression in alfalfa (Medicago sativa) which conditions not only above-ground plant shoot architecture but also root development and growth. Delayed flowering is an important trait for forage species, as it allows harvesting of high-quality forage for a longer time before nutritional values decline due to plant architecture changes related to flowering onset. Despite the relevance of delayed flowering, this trait has not yet been thoroughly exploited in alfalfa. This is mainly due to its complex genetics, sensitivity to inbreeding and to the fact that delayed flowering would be only advantageous if it allowed increased forage quality without compromising seed production. To develop new delayed-flowering varieties, we have characterized the three TERMINAL FLOWERING 1 (TFL1) family of genes in alfalfa: MsTFL1A, MsTFL1B and MsTFL1C. Constitutive expression of MsTFL1A in Arabidopsis caused late flowering and changes in inflorescence architecture, indicating that MsTFL1A is the ortholog of Arabidopsis TFL1. Overexpression of MsTFL1A in alfalfa consistently led to delayed flowering in both controlled and natural field conditions, coupled to an increase in leaf/stem ratio, a common indicator of forage quality. Additionally, overexpression of MsTFL1A reduced root development, reinforcing the role of MsTFL1A not only as a flowering repressor but also as a regulator of root development.We conclude that the precise manipulation of MsTFL1A gene expression may represent a powerful tool to improve alfalfa forage quality.

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.

Common signatures of selection reveal target loci for breeding across soybean populations

Understanding the underlying genetic bases of yield-related selection and distinguishing these changes from genetic drift are critical for both improved understanding and future success of plant breeding. Soybean [Glycine max (L.) Merr.] is a key species for world food security, yet knowledge of the mechanism of selective breeding in soybean, such as the century-long program of artificial selection in U.S. soybean germplasm, is currently limited to certain genes and loci. Here, we identify genome-wide signatures of selection in separate populations of soybean subjected to artificial selection for increased yield by multiple breeding programs in the United States. We compared the alternative soybean breeding population (AGP) created by USDA-ARS to the conventional public soybean lines (CGP) developed at three different stages of breeding (ancestral, intermediate, and elite) to identify shared signatures of selection and differentiate these from drift. The results showed a strong selection for specific haplotypes identified by single site frequency and haplotype homozygosity methods. A set of common selection signatures was identified in both AGP and CGP that supports the hypothesis that separate breeding programs within similar environments coalesce on the fixation of the same key haplotypes. Signatures unique to each breeding program were observed. These results raise the possibility that selection analysis can allow the identification of favorable alleles to enhance directed breeding approaches.

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.

Natural and artificial selection of multiple alleles revealed through genomic analyses

Genome-to-phenome research in agriculture aims to improve crops through in silico predictions. Genome-wide association study (GWAS) is potent in identifying genomic loci that underlie important traits. As a statistical method, increasing the sample quantity, data quality, or diversity of the GWAS dataset positively impacts GWAS power. For more precise breeding, concrete candidate genes with exact functional variants must be discovered. Many post-GWAS methods have been developed to narrow down the associated genomic regions and, ideally, to predict candidate genes and causative mutations (CMs). Historical natural selection and breeding-related artificial selection both act to change the frequencies of different alleles of genes that control phenotypes. With higher diversity and more extensive GWAS datasets, there is an increased chance of multiple alleles with independent CMs in a single causal gene. This can be caused by the presence of samples from geographically isolated regions that arose during natural or artificial selection. This simple fact is a complicating factor in GWAS-driven discoveries. Currently, none of the existing association methods address this issue and need to identify multiple alleles and, more specifically, the actual CMs. Therefore, we developed a tool that computes a score for a combination of variant positions in a single candidate gene and, based on the highest score, identifies the best number and combination of CMs. The tool is publicly available as a Python package on GitHub, and we further created a web-based Multiple Alleles discovery (MADis) tool that supports soybean and is hosted in SoyKB (https://soykb.org/SoybeanMADisTool/). We tested and validated the algorithm and presented the utilization of MADis in a pod pigmentation L1 gene case study with multiple CMs from natural or artificial selection. Finally, we identified a candidate gene for the pod color L2 locus and predicted the existence of multiple alleles that potentially cause loss of pod pigmentation. In this work, we show how a genomic analysis can be employed to explore the natural and artificial selection of multiple alleles and, thus, improve and accelerate crop breeding in agriculture.

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.

Characterization of terminal flowering cowpea (Vigna unguiculata (L.) Walp.) mutants obtained by induced mutagenesis digs out the loss-of-function of phosphatidylethanolamine-binding protein

Cowpea (Vigna unguiculata (L.) Walp) is one of the major food legume crops grown extensively in arid and semi-arid regions of the world. The determinate habit of cowpea has many advantages over the indeterminate and is well adapted to modern farming systems. Mutation breeding is an active research area to develop the determinate habit of cowpea. The present study aimed to develop new determinate habit mutants with terminal flowering (TFL) in locally well-adapted genetic backgrounds. Consequently, the seeds of popular cowpea cv P152 were irradiated with doses of gamma rays (200, 250, and, 300 Gy), and the M1 populations were grown. The M2 populations were produced from the M1 progenies and selected determinate mutants (TFLCM-1 and TFLCM-2) from the M2 generation (200 Gy) were forwarded up to the M5 generation to characterize the mutants and simultaneously they were crossed with P152 to develop a MutMap population. In the M5 generation, determinate mutants (80-81 days) were characterized by evaluating the TFL growth habit, longer peduncles (30.75-31.45 cm), erect pods (160°- 200°), number of pods per cluster (4-5 nos.), and early maturity. Further, sequencing analysis of the VuTFL1 gene in the determinate mutants and MutMap population revealed a single nucleotide transversion (A-T at 1196 bp) in the fourth exon and asparagine (N) to tyrosine (Y) amino acid change at the 143rd position of phosphatidylethanolamine-binding protein (PEBP). Notably, the loss of function PEPB with a higher confidence level modification of anti-parallel beta-sheets and destabilization of the protein secondary structure was observed in the mutant lines. Quantitative real-time PCR (qRT-PCR) analysis showed that the VuTFL1 gene was downregulated at the flowering stage in TFL mutants. Collectively, the insights garnered from this study affirm the effectiveness of induced mutation in modifying the plant's ideotype. The TFL mutants developed during this investigation have the potential to serve as a valuable resource for fostering determinate traits in future cowpea breeding programs and pave the way for mechanical harvesting.

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.

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.

UV-B irradiation-activated E3 ligase GmILPA1 modulates gibberellin catabolism to increase plant height in soybean

Plant height is a key agronomic trait that affects yield and is controlled by both phytohormone gibberellin (GA) and ultraviolet-B (UV-B) irradiation. However, whether and how plant height is modulated by UV-B-mediated changes in GA metabolism are not well understood. It has not been reported that the E3 ubiquitin ligase Anaphase Promoting Complex/Cyclosome (APC/C) is involved in the regulation of plant growth in response to environmental factors. We perform a forward genetic screen in soybean and find that a mutation in Glycine max Increased Leaf Petiole Angle1 (GmILPA1), encoding a subunit of the APC/C, lead to dwarfism under UV-B irradiation. UV-B promotes the accumulation of GmILPA1, which ubiquitinate the GA catabolic enzyme GA2 OXIDASE-like (GmGA2ox-like), resulting in its degradation in a UV-B-dependent manner. Another E3 ligase, GmUBL1, also ubiquitinate GmGA2ox-like and enhance the GmILPA1-mediated degradation of GmGA2ox-like, which suggest that GmILPA1-GmGA2ox-like module counteract the UV-B-mediated reduction of bioactive GAs. We also determine that GmILPA1 is a target of selection during soybean domestication and breeding. The deletion (Indel-665) in the promoter might facilitate the adaptation of soybean to high UV-B irradiation. This study indicates that an evolutionary GmILPA1 variant has the capability to develop ideal plant architecture with soybean cultivars.

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.

Genome-Wide Characterization of PEBP Gene Family and Functional Analysis of TERMINAL FLOWER 1 Homologs in Macadamia integrifolia

Edible Macadamia is one of the most important commercial nut trees cultivated in many countries, but its large tree size and long juvenile period pose barriers to commercial cultivation. The short domestication period and well-annotated genome of Macadamia integrifolia create great opportunities to breed commercial varieties with superior traits. Recent studies have shown that members of the phosphatidylethanolamine binding protein (PEBP) family play pivotal roles in regulating plant architecture and flowering time in various plants. In this study, thirteen members of MiPEBP were identified in the genome of M. integrifolia, and they are highly similarity in both motif and gene structure. A phylogenetic analysis divided the MiPEBP genes into three subfamilies: MFT-like, FT-like and TFL1-like. We subsequently identified two TERMINAL FLOWER 1 homologues from the TFL1-like subfamily, MiTFL1 and MiTFL1-like, both of which were highly expressed in stems and vegetative shoots, while MiTFL1-like was highly expressed in young leaves and early flowers. A subcellular location analysis revealed that both MiTFL1 and MiTFL1-like are localized in the cytoplasm and nucleus. The ectopic expression of MiTFL1 can rescue the early-flowering and terminal-flower phenotypes in the tfl1-14 mutant of Arabidopsis thaliana, and it indicates the conserved functions in controlling the inflorescence architecture and flowering time. This study will provide insight into the isolation of PEBP family members and the key targets for breeding M. integrifolia with improved traits in plant architecture and flowering time.

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.

An Improved Genome-Wide Association Procedure Explores Gene-Allele Constitutions and Evolutionary Drives of Growth Period Traits in the Global Soybean Germplasm Population

In soybeans (Glycine max (L.) Merr.), their growth periods, DSF (days of sowing-to-flowering), and DFM (days of flowering-to-maturity) are determined by their required accumulative day-length (ADL) and active temperature (AAT). A sample of 354 soybean varieties from five world eco-regions was tested in four seasons in Nanjing, China. The ADL and AAT of DSF and DFM were calculated from daily day-lengths and temperatures provided by the Nanjing Meteorological Bureau. The improved restricted two-stage multi-locus genome-wide association study using gene-allele sequences as markers (coded GASM-RTM-GWAS) was performed. (i) For DSF and its related ADL_DSF and AAT_DSF, 130-141 genes with 384-406 alleles were explored, and for DFM and its related ADL_DFM and AAT_DFM, 124-135 genes with 362-384 alleles were explored, in a total of six gene-allele systems. DSF shared more ADL and AAT contributions than DFM. (ii) Comparisons between the eco-region gene-allele submatrices indicated that the genetic adaptation from the origin to the geographic sub-regions was characterized by allele emergence (mutation), while genetic expansion from primary maturity group (MG)-sets to early/late MG-sets featured allele exclusion (selection) without allele emergence in addition to inheritance (migration). (iii) Optimal crosses with transgressive segregations in both directions were predicted and recommended for breeding purposes, indicating that allele recombination in soybean is an important evolutionary drive. (iv) Genes of the six traits were mostly trait-specific involved in four categories of 10 groups of biological functions. GASM-RTM-GWAS showed potential in detecting directly causal genes with their alleles, identifying differential trait evolutionary drives, predicting recombination breeding potentials, and revealing population gene networks.

Characterization of CcTFL1 Governing Plant Architecture in Pigeon pea (Cajanus cajan (L.) Millsp.)

Growth habits are among the essential adaptive traits acted upon by evolution during plant speciation. They have brought remarkable changes in the morphology and physiology of plants. Inflorescence architecture varies greatly between wild relatives and cultivars of pigeon pea. The present study isolated the CcTFL1 (Terminal Flowering Locus 1) locus using six varieties showing determinate (DT) and indeterminate (IDT) growth habits. Multiple alignments of CcTFL1 sequences revealed the presence of InDel, which describes a 10 bp deletion in DT varieties. At the same time, IDT varieties showed no deletion. InDel altered the translation start point in DT varieties, resulting in the shortening of exon 1. This InDel was validated in ten varieties of cultivated species and three wild relatives differing in growth habits. The predicted protein structure showed the absence of 27 amino acids in DT varieties, which was reflected in mutant CcTFL1 by the absence of two α-helices, a connecting loop, and shortened β-sheet. By subsequent motif analysis, it was found that the wild-type protein had a phosphorylation site for protein kinase C, but the mutant protein did not. In silico analysis revealed that the InDel-driven deletion of amino acids spans, containing a phosphorylation site for kinase protein, may have resulted in the non-functionality of the CcTFL1 protein, rendering the determinate growth habit. This characterization of the CcTFL1 locus could be used to modulate growth habits through genome editing.

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.

Epistatic interaction between CsCEN and CsSHBY in regulating indeterminate/determinate growth of lateral branch in cucumber

KEY MESSAGE

Two genetic loci, det-ma (CsCEN) and det-lb, showed epistatic interaction on indeterminate/determinate growth of LB in cucumber. CsSHBY was identified as the candidate gene for det-lb locus. Plant architecture depends on the spatial regulation of meristems from both main axis (MA) and lateral branches (LBs). Fate (indeterminate or determinate) of these meristems is a crucial source of architectural diversity determining crop productivity and management. CENTRORADIALIS/TERMINAL FLOWER 1/SELF-PRUNING (CETS) gene family have been well known as pivotal regulators for indeterminate/determinate growth of MA. Nevertheless, genes that regulate LB indeterminacy/determinacy remained unclear. Cucumber (Cucumis sativus L.) has typical monopodial growth and multiple lateral branches. Both MA and LBs had indeterminate or determinate growth, and indeterminate/determinate growth of LB was controlled by two distinct loci, det-ma (CsCEN) and det-lb. In our study, based on bulked segregant analysis (BSA) method, the det-lb locus was mapped on a 60.6 kb region on chromosome 1 harboring only one gene CsaV3_1G044330, which encoded a putative vacuolar-sorting protein (designated as CsSHBY). Multipoint mutations in CsSHBY were identified in D082 and D226, compared with CCMC, including nonsynonymous SNP mutations and a 6-bp deletion in exons. Further, qPCR showed that CsSHBY was highly expressed in lateral bud of CCMC, suggesting that CsSHBY might play an active role in regulating indeterminate/determinate growth of LB. Genetic analyses showed that det-ma (CsCEN) had an epistatic effect on det-lb (CsSHBY), and CsCEN could activate CsSHBY promoter by Dual luciferase and GUS activity assays. Meanwhile, Cscen or Csshby was found to influence auxin contents and CsYUCs and CsPINs expression levels. These findings provided new insights into precisely optimizing plant architecture for yield improvements.

The developmental dynamics in cool season legumes with focus on chickpea

Chickpea is one of the most widely consumed grain legume world-wide. Advances in next-generation sequencing and genomics tools have led to genetic dissection and identification of potential candidate genes regulating agronomic traits in chickpea. However, the developmental particularities and its potential in reforming the yield and nutritional value remain largely unexplored. Studies in crops such as rice, maize, tomato and pea have highlighted the contribution of key regulator of developmental events in yield related traits. A comprehensive knowledge on the development aspects of a crop can pave way for new vistas to explore. Pea and Medicago are the close relatives of genus Cicer and the basic developmental events in these legumes are similar. However, there are some distinct developmental features in chickpea which hold potential for future crop improvement endeavours. The global chickpea germplasm encompasses wide range of diversities in terms of morphology at both vegetative and reproductive stages. There is an immediate need for understanding the genetic and molecular basis of this diversity and utilizing them for the yield contributing trait improvement. The review discusses some of the key developmental events which have potential in yield enhancement and the lessons which can be learnt from model legumes in this regard.

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.

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.

Time series canopy phenotyping enables the identification of genetic variants controlling dynamic phenotypes in soybean

Advances in plant phenotyping technologies are dramatically reducing the marginal costs of collecting multiple phenotypic measurements across several time points. Yet, most current approaches and best statistical practices implemented to link genetic and phenotypic variation in plants have been developed in an era of single-time-point data. Here, we used time-series phenotypic data collected with an unmanned aircraft system for a large panel of soybean (Glycine max (L.) Merr.) varieties to identify previously uncharacterized loci. Specifically, we focused on the dissection of canopy coverage (CC) variation from this rich data set. We also inferred the speed of canopy closure, an additional dimension of CC, from the time-series data, as it may represent an important trait for weed control. Genome-wide association studies (GWASs) identified 35 loci exhibiting dynamic associations with CC across developmental stages. The time-series data enabled the identification of 10 known flowering time and plant height quantitative trait loci (QTLs) detected in previous studies of adult plants and the identification of novel QTLs influencing CC. These novel QTLs were disproportionately likely to act earlier in development, which may explain why they were missed in previous single-time-point studies. Moreover, this time-series data set contributed to the high accuracy of the GWASs, which we evaluated by permutation tests, as evidenced by the repeated identification of loci across multiple time points. Two novel loci showed evidence of adaptive selection during domestication, with different genotypes/haplotypes favored in different geographic regions. In summary, the time-series data, with soybean CC as an example, improved the accuracy and statistical power to dissect the genetic basis of traits and offered a promising opportunity for crop breeding with quantitative growth curves.

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.

Divergence of functions and expression patterns of soybean bZIP transcription factors

Soybean (Glycine max) is a major protein and oil crop. Soybean basic region/leucine zipper (bZIP) transcription factors are involved in many regulatory pathways, including yield, stress responses, environmental signaling, and carbon-nitrogen balance. Here, we discuss the members of the soybean bZIP family and their classification: 161 members have been identified and clustered into 13 groups. Our review of the transcriptional regulation and functions of soybean bZIP members provides important information for future study of bZIP transcription factors and genetic resources for soybean breeding.

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.

Phenotypic analysis of Longya-10 × pale flax hybrid progeny and identification of candidate genes regulating prostrate/erect growth in flax plants

Flax is a dual-purpose crop that is important for oil and fiber production. The growth habit is one of the crucial targets of selection during flax domestication. Wild hybridization between cultivated flax and wild flax can produce superior germplasms for flax breeding and facilitate the study of the genetic mechanism underlying agronomically important traits. In this study, we used pale flax, Linum grandiflorum, and L. perenne to pollinate Longya-10. Only pale flax interspecific hybrids were obtained, and the trait analysis of the F₁ and F₂ generations showed that the traits analyzed in this study exhibited disparate genetic characteristics. In the F₁ generation, only one trait, i.e., the number of capsules per plant (140) showed significant heterosis, while the characteristics of other traits were closely associated with those of the parents or a decline in hybrid phenotypes. The traits of the F₂ generation were widely separated, and the variation coefficient ranged from 9.96% to 146.15%. The quantitative trait locus underlying growth habit was preliminarily found to be situated on chromosome 2 through Bulked-segregant analysis sequencing. Then linkage mapping analysis was performed to fine-map GH2.1 to a 23.5-kb interval containing 4 genes. Among them, L.us.o.m.scaffold22.109 and L.us.o.m.scaffold22.112 contained nonsynonymous SNPs with Δindex=1. Combined with the qRT-PCR results, the two genes might be possible candidate genes for GH2.1. This study will contribute to the development of important germplasms for flax breeding, which would facilitate the elucidation of the genetic mechanisms regulating the growth habit and development of an ideal architecture for the flax plant.

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.

Genomic analyses of rice bean landraces reveal adaptation and yield related loci to accelerate breeding

Rice bean (Vigna umbellata) is an underexploited domesticated legume crop consumed for dietary protein in Asia, yet little is known about the genetic diversity of this species. Here, we present a high-quality reference genome for a rice bean landrace (FF25) built using PacBio long-read data and a Hi-C chromatin interaction map, and assess the phylogenetic position and speciation time of rice bean within the Vigna genus. We sequence 440 landraces (two core collections), and GWAS based on data for growth sites at three widely divergent latitudes reveal loci associated with flowering and yield. Loci harboring orthologs of FUL (FRUITFULL), FT (FLOWERING LOCUS T), and PRR3 (PSEUDO-RESPONSE REGULATOR 3) contribute to the adaptation of rice bean from its low latitude center of origin towards higher latitudes, and the landraces which pyramid early-flowering alleles for these loci display maximally short flowering times. We also demonstrate that copy-number-variation for VumCYP78A6 can regulate seed-yield traits. Intriguingly, 32 landraces collected from a mountainous region in South-Central China harbor a recently acquired InDel in TFL1 (TERMINAL FLOWER1) affecting stem determinacy; these materials also have exceptionally high values for multiple human-desired traits and could therefore substantially advance breeding efforts to improve rice bean.

Rice bean is an underexploited legume crop that has many desirable properties against bio and abiotic stresses. Here, the authors report the genome assembly of this species, conduct population genetics studies and reveal the genetic variations associated with adaptation and yield traits.

QTL and candidate genes for heterophylly in soybean based on two populations of recombinant inbred lines

Heterophylly, the existence of different leaf shapes and sizes on the same plant, has been observed in many flowering plant species. Yet, the genetic characteristics and genetic basis of heterophylly in soybean remain unknown. Here, two populations of recombinant inbred lines (RILs) with distinctly different leaf shapes were used to identify loci controlling heterophylly in two environments. The ratio of apical leaf shape (LSUP) to basal leaf shape (LSDOWN) at the reproductive growth stage (RLS) was used as a parameter for classifying heterophylly. A total of eight QTL were detected for RLS between the two populations and four of them were stably identified in both environments. Among them, qRLS20 had the largest effect in the JS population, with a maximum LOD value of 46.9 explaining up to 47.2% of phenotypic variance. This locus was located in the same genomic region as the basal leaf shape QTL qLSDOWN20 on chromosome 20. The locus qRLS19 had the largest effect in the JJ population, with a maximum LOD value of 15.2 explaining up to 27.0% of phenotypic variance. This locus was located in the same genomic region as the apical leaf shape QTL qLSUP19 on chromosome 19. Four candidate genes for heterophylly were identified based on sequence differences among the three parents of the two mapping populations, RT-qPCR analysis, and gene functional annotation analysis. The QTL and candidate genes detected in this study lay a foundation for further understanding the genetic mechanism of heterophylly and are invaluable in marker-assisted breeding.

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.

Thirty Years of Mungbean Genome Research: Where Do We Stand and What Have We Learned?

Mungbean is a socioeconomically important legume crop in Asia that is currently in high demand by consumers and industries both as dried beans and in plant-based protein foods. Marker-assisted and genomics-assisted breeding are promising approaches to efficiently and rapidly develop new cultivars with improved yield, quality, and resistance to biotic and abiotic stresses. Although mungbean was at the forefront of research at the dawn of the plant genomics era 30 years ago, the crop is a "slow runner" in genome research due to limited genomic resources, especially DNA markers. Significant progress in mungbean genome research was achieved only within the last 10 years, notably after the release of the VC1973A draft reference genome constructed using next-generation sequencing technology, which enabled fast and efficient DNA marker development, gene mapping, and identification of candidate genes for complex traits. Resistance to biotic stresses has dominated mungbean genome research to date; however, research is on the rise. In this study, we provide an overview of the past progress and current status of mungbean genomics research. We also discuss and evaluate some research results to provide a better understanding of mungbean genomics.

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.

GmMDE genes bridge the maturity gene E1 and florigens in photoperiodic regulation of flowering in soybean

Soybean (Glycine max) is highly sensitive to photoperiod, which affects flowering time and plant architecture and thus limits the distribution range of elite soybean cultivars. The major maturity gene E1 confers the most prominent effect on photoperiod sensitivity, but its downstream signaling pathway remains largely unknown. Here, we confirm that the encoded E1 protein is a transcriptional repressor. The expression of seven GmMDE genes (Glycine max MADS-box genes downregulated by E1) was suppressed when E1 was overexpressed and promoted when E1 was knocked out through clustered regularly-interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 (Cas9)-mediated mutagenesis. These GmMDEs exhibited similar tissue specificity and expression patterns, including in response to photoperiod, E1 expression, and E1 genotype. E1 repressed GmMDE promoter activity. Results for two GmMDEs showed that E1 epigenetically silences their expression by directly binding to their promoters to increase H3K27me3 levels. The overexpression of GmMDE06 promoted flowering and post-flowering termination of stem growth. The late flowering phenotype of E1-overexpressing soybean lines was reversed by the overexpression of GmMDE06, placing GmMDE06 downstream of E1. The overexpression of GmMDE06 increased the expression of the soybean FLOWERING LOCUS T orthologs GmFT2a and GmFT5a, leading to feedback upregulation of GmMDE, indicating that GmMDE and GmFT2a/GmFT5a form a positive regulatory feedback loop promoting flowering. GmMDE06 also promoted post-flowering termination of stem growth by repressing the expression of the shoot identity gene Dt1. The E1-GmMDEs-GmFT2a/5a-Dt1 signaling pathway illustrates how soybean responds to photoperiod by modulating flowering time and post-flowering stem termination.

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.

Hard versus soft selective sweeps during domestication and improvement in soybean

Genome scans for selection can provide an efficient way to dissect the genetic basis of domestication traits and understand mechanisms of adaptation during crop evolution. Selection involving soft sweeps (simultaneous selection for multiple alleles) is probably common in plant genomes but is under-studied, and few if any studies have systematically scanned for soft sweeps in the context of crop domestication. Using genome resequencing data from 302 wild and domesticated soybean accessions, we conducted selection scans using five widely employed statistics to identify selection candidates under classical (hard) and soft sweeps. Across the genome, inferred hard sweeps are predominant in domesticated soybean landraces and improved varieties, whereas soft sweeps are more prevalent in a representative subpopulation of the wild ancestor. Six domestication-related genes, representing both hard and soft sweeps and different stages of domestication, were used as positive controls to assess the detectability of domestication-associated sweeps. Performance of various test statistics suggests that differentiation-based (FST ) methods are robust for detecting complete hard sweeps, and that LD-based strategies perform well for identifying recent/ongoing sweeps; however, none of the test statistics detected a known soft sweep we previously documented at the domestication gene Dt1. Genome scans yielded a set of 66 candidate loci that were identified by both differentiation-based and LD-based (iHH) methods; notably, this shared set overlaps with many previously identified QTLs for soybean domestication/improvement traits. Collectively, our results will help to advance genetic characterizations of soybean domestication traits and shed light on selection modes involved in adaptation in domesticated plant species.

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.

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.