Soybean adaption to high-latitude regions is associated with natural variations of GmFT2b, an ortholog of FLOWERING LOCUS T

TWAS facilitates gene-scale trait genetic dissection through gene expression, structural variations, and alternative splicing in soybean

A genome-wide association study (GWAS) identifies trait-associated loci, but identifying the causal genes can be a bottleneck, due in part to slow decay of linkage disequilibrium (LD). A transcriptome-wide association study (TWAS) addresses this issue by identifying gene expression-phenotype associations or integrating gene expression quantitative trait loci with GWAS results. Here, we used self-pollinated soybean (Glycine max [L.] Merr.) as a model to evaluate the application of TWAS to the genetic dissection of traits in plant species with slow LD decay. We generated RNA sequencing data for a soybean diversity panel and identified the genetic expression regulation of 29 286 soybean genes. Different TWAS solutions were less affected by LD and were robust to the source of expression, identifing known genes related to traits from different tissues and developmental stages. The novel pod-color gene L2 was identified via TWAS and functionally validated by genome editing. By introducing a new exon proportion feature, we significantly improved the detection of expression variations that resulted from structural variations and alternative splicing. As a result, the genes identified through our TWAS approach exhibited a diverse range of causal variations, including SNPs, insertions or deletions, gene fusion, copy number variations, and alternative splicing. Using this approach, we identified genes associated with flowering time, including both previously known genes and novel genes that had not previously been linked to this trait, providing insights complementary to those from GWAS. In summary, this study supports the application of TWAS for candidate gene identification in species with low rates of LD decay.

Light recovery after maize harvesting promotes soybean flowering in a maize-soybean relay strip intercropping system

Moving from sole cropping to intercropping is a transformative change in agriculture, contributing to yield. Soybeans adapt to light conditions in intercropping by adjusting the onset of reproduction and the inflorescence architecture to optimize reproductive success. Maize-soybean strip intercropping (MS), maize-soybean relay strip intercropping (IS), and sole soybean (SS) systems are typical soybean planting systems with significant differences in light environments during growth periods. To elucidate the effect of changes in the light environment on soybean flowering processes and provide a theoretical basis for selecting suitable varieties in various planting systems to improve yields, field experiments combining planting systems (IS, MS, and SS) and soybean varieties (GQ8, GX7, ND25, and NN996) were conducted in 2021 and 2022. Results showed that growth recovery in the IS resulted in a balance in the expression of TERMINAL FLOWER 1 (TFL1) and FLOWERING LOCUS T (FT) in the meristematic tissues of soybeans, which promoted the formation of new branches or flowers. IS prolonged the flowering time (2-7 days) and increased the number of forming flowers compared with SS (93.0 and 169%) and MS (67.3 and 103.3%) at the later soybean flowering stage. The higher carbon and nitrogen content in the middle and bottom canopies of soybean contributed to decreased flower abscission by 26.7 and 30.2%, respectively, compared with SS. Canopy light environment recovery promoted branch and flower formation and transformation of flowers into pods with lower flower-pod abscission, which contributed to elevating soybean yields in late-maturing and multibranching varieties (ND25) in IS.

Floral-promoting GmFT homologs trigger photoperiodic after-effects: An important mechanism for early-maturing soybean varieties to regulate reproductive development and adapt to high latitudes

Soybean (Glycine max) is a typical short-day plant, but has been widely cultivated in high-latitude long-day (LD) regions because of the development of early-maturing genotypes which are photoperiod-insensitive. However, some early-maturing varieties exhibit significant responses to maturity under different daylengths but not for flowering, depicting an evident photoperiodic after-effect, a poorly understood mechanism. In this study, we investigated the postflowering responses of 11 early-maturing soybean varieties to various preflowering photoperiodic treatments. We confirmed that preflowering SD conditions greatly promoted maturity and other postflowering developmental stages. Soybean homologs of FLOWERING LOCUS T (FT), including GmFT2a, GmFT3a, GmFT3b and GmFT5a, were highly accumulated in leaves under preflowering SD treatment. More importantly, they maintained a high expression level after flowering even under LD conditions. E1 RNAi and GmFT2a overexpression lines showed extremely early maturity regardless of preflowering SD and LD treatments due to constitutively high levels of floral-promoting GmFT homolog expression throughout their life cycle. Collectively, our data indicate that high and stable expression of floral-promoting GmFT homologs play key roles in the maintenance of photoperiodic induction to promote postflowering reproductive development, which confers early-maturing varieties with appropriate vegetative growth and shortened reproductive growth periods for adaptation to high latitudes.

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

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

The FLOWERING LOCUS T 5b positively regulates photoperiodic flowering and improves the geographical adaptation of soybean

Plants can sense the photoperiod to flower at the right time. As a sensitive short-day crop, soybean (Glycine max) flowering varies greatly depending on photoperiods, affecting yields. Adaptive changes in soybeans rely on variable genetic loci such as E1 and FLOWERING LOCUS T orthologs. However, the precise coordination and control of these molecular components remain largely unknown. In this study, we demonstrate that GmFT5b functions as a crucial factor for soybean flowering. Overexpressed or mutated GmFT5b resulted in significantly early or later flowering, altering expression profiles for several downstream flowering-related genes under a long-day photoperiod. GmFT5b interacts with the transcription factor GmFDL15, suggesting transcriptional tuning of flowering time regulatory genes via the GmFT5b/GmFDL15 complex. Notably, GmFT5a partially compensated for GmFT5b function, as ft5a ft5b double mutants exhibited an enhanced late-flowering phenotype. Association mapping revealed that GmFT5b was associated with flowering time, maturity, and geographical distribution of soybean accessions, all associated with the E1 locus. Therefore, GmFT5b is a valuable target for enhancing regional adaptability. Natural variants or multiple mutants in this region can be utilized to generate optimized soybean varieties with precise flowering times.

Soybean WRINKLED1 protein GmWRI1a promotes flowering under long-day conditions via regulating expressions of flowering-related genes

Flowering time is an important agronomic character that influences the adaptability and yield of soybean [Glycine max (L.) Merrill]. WRINKLED 1 (WRI1) plays an important regulatory role in plant growth and development. In this study, we found that the expression of GmWIR1a could be induced by long days. Compared with the wild type, transgenic soybean overexpressing GmWRI1a showed earlier flowering and maturity under long days but no significant changes under short days. Overexpression of GmWRI1a led to up-regulated expression of genes involved in the regulation of flowering time. The GmWRI1a protein was able to directly bind to the promoter regions of GmAP1, GmFUL1a, GmFUL2 and up-regulated their expression. GmCOL3 was identified by yeast one-hybrid library screening using the GmWRI1a promoter as bait. GmCOL3 was revealed to be a nucleus-localized protein that represses the transcription of GmWRI1a. Expression of GmCOL3 was induced by short days. Taken together, the results show that overexpression of GmWRI1a promotes flowering under long days by promoting the transcriptional activity of flowering-related genes in soybean, and that GmCOL3 binds to the GmWRI1a promoter and directly down-regulates its transcription. This discovery reveals a new function for GmWRI1a, which regulates flowering and maturity in soybean.

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

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

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.

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.

Molecular breeding for improvement of photothermal adaptability in soybean

Soybean (Glycine max (L.) Merr.) is a typical short-day and temperate crop that is sensitive to photoperiod and temperature. Responses of soybean to photothermal conditions determine plant growth and development, which affect its architecture, yield formation, and capacity for geographic adaptation. Flowering time, maturity, and other traits associated with photothermal adaptability are controlled by multiple major-effect and minor-effect genes and genotype-by-environment interactions. Genetic studies have identified at least 11 loci (E1-E4, E6-E11, and J) that participate in photoperiodic regulation of flowering time and maturity in soybean. Molecular cloning and characterization of major-effect flowering genes have clarified the photoperiod-dependent flowering pathway, in which the photoreceptor gene phytochrome A, circadian evening complex (EC) components, central flowering repressor E1, and FLOWERING LOCUS T family genes play key roles in regulation of flowering time, maturity, and adaptability to photothermal conditions. Here, we provide an overview of recent progress in genetic and molecular analysis of traits associated with photothermal adaptability, summarizing advances in molecular breeding practices and tools for improving these traits. Furthermore, we discuss methods for breeding soybean varieties with better adaptability to specific ecological regions, with emphasis on a novel strategy, the Potalaization model, which allows breeding of widely adapted soybean varieties through the use of multiple molecular tools in existing elite widely adapted varieties.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01406-z.

A novel MORN-motif type gene GmMRF2 controls flowering time and plant height of soybean

The flowering time of soybean is a highly important agronomic characteristic, which affects the adaptability and yield. AtMRF1, a MORN-repeat motif gene, acts as a floral promoter in Arabidopsis, its functions in soybean are not yet understood. Here, we employed qRT-PCR to analyze the tissue expression patten of MRF1 homologs in soybean and determined that the GmMRF2 gene, containing a MORN-motif, highly expressed in the shoot and responded to photoperiod. GmMRF2 overexpression soybean lines exhibited earlier flowering time under long-day (LD) conditions, and increased plant height under both LD and short-day (SD) conditions compared to wild-type (WT) plants. The expression levels of gibberellic acid (GA) pathway genes that positively regulate plant height genes and flowering-promoting genes were up-regulated in the GmMRF2 overexpression lines, were up-regulated in the GmMRF2 overexpression lines. Further study revealed that GmMRF2 interacted with GmTCP15 to co-induce the expression of GmSOC1b. Together, our results preliminarily reveal the functions and mechanisms of GmMRF2 in regulating flowering time and plant height, provide a new promising gene for soybean crop improvement.

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.

Haplotype Analysis of GmSGF14 Gene Family Reveals Its Roles in Photoperiodic Flowering and Regional Adaptation of Soybean

Flowering time and photoperiod sensitivity are fundamental traits that determine soybean adaptation to a given region or a wide range of geographic environments. The General Regulatory Factors (GRFs), also known as 14-3-3 family, are involved in protein-protein interactions in a phosphorylation-dependent manner, thus regulating ubiquitous biological processes, such as photoperiodic flowering, plant immunity and stress response. In this study, 20 soybean GmSGF14 genes were identified and divided into two categories according to phylogenetic relationships and structural characteristics. Real-time quantitative PCR analysis revealed that GmSGF14g, GmSGF14i, GmSGF14j, GmSGF14k, GmSGF14m and GmSGF14s were highly expressed in all tissues compared to other GmSGF14 genes. In addition, we found that the transcript levels of GmSGF14 family genes in leaves varied significantly under different photoperiodic conditions, indicating that their expression responds to photoperiod. To explore the role of GmSGF14 in the regulation of soybean flowering, the geographical distribution of major haplotypes and their association with flowering time in six environments among 207 soybean germplasms were studied. Haplotype analysis confirmed that the GmSGF14m^H4 harboring a frameshift mutation in the 14-3-3 domain was associated with later flowering. Geographical distribution analysis demonstrated that the haplotypes related to early flowering were frequently found in high-latitude regions, while the haplotypes associated with late flowering were mostly distributed in low-latitude regions of China. Taken together, our results reveal that the GmSGF14 family genes play essential roles in photoperiodic flowering and geographical adaptation of soybean, providing theoretical support for further exploring the function of specific genes in this family and varietal improvement for wide adaptability.

Transcriptional profiling reveals a critical role of GmFT2a in soybean staygreen syndrome caused by the pest Riptortus pedestris

Soybean staygreen syndrome, characterized by delayed leaf and stem senescence, abnormal pods, and aborted seeds, has recently become a serious and prominent problem in soybean production. Although the pest Riptortus pedestris has received increasing attention as the possible cause of staygreen syndrome, the mechanism remains unknown. Here, we clarify that direct feeding by R. pedestris, not transmission of a pathogen by this pest, is the primary cause of typical soybean staygreen syndrome and that critical feeding damage occurs at the early pod stage. Transcriptome profiling of soybean indicated that many signal transduction pathways, including photoperiod, hormone, defense response, and photosynthesis, respond to R. pedestris infestation. Importantly, we discovered that members of the FLOWERING LOCUS T (FT) gene family were suppressed by R. pedestris infestation, and overexpression of floral inducer GmFT2a attenuates staygreen symptoms by mediating soybean defense response and photosynthesis. Together, our findings systematically illustrate the association between pest infestation and soybean staygreen syndrome and provide the basis for establishing a targeted soybean pest prevention and control system.

Genome-wide signatures of the geographic expansion and breeding of soybean

Soybean is a leguminous crop that provides oil and protein. Exploring the genomic signatures of soybean evolution is crucial for breeding varieties with improved adaptability to environmental extremes. We analyzed the genome sequences of 2,214 soybeans and proposed a soybean evolutionary route, i.e., the expansion of annual wild soybean (Glycine soja Sieb. & Zucc.) from southern China and its domestication in central China, followed by the expansion and local breeding selection of its landraces (G. max (L.) Merr.). We observed that the genetic introgression in soybean landraces was mostly derived from sympatric rather than allopatric wild populations during the geographic expansion. Soybean expansion and breeding were accompanied by the positive selection of flowering time genes, including GmSPA3c. Our study sheds light on the evolutionary history of soybean and provides valuable genetic resources for its future breeding.

CRISPR/Cas9-Mediated Multiple Knockouts in Abscisic Acid Receptor Genes Reduced the Sensitivity to ABA during Soybean Seed Germination

Abscisic acid (ABA) is an important plant hormone that regulates numerous functions in plant growth, development, and stress responses. Several proteins regulate the ABA signal transduction mechanism in response to environmental stress. Among them, the PYR1/PYL/RCAR family act as ABA receptors. This study used the CRISPR/Cas9 gene-editing system with a single gRNA to knock out three soybean PYL genes: GmPYL17, GmPYL18, and GmPYL19. The gRNA may efficiently cause varying degrees of deletion of GmPYL17, GmPYL18, and GmPYL19 gene target sequences, according to the genotyping results of T0 plants. A subset of induced alleles was successfully transferred to progeny. In the T2 generation, we obtained double and triple mutant genotypes. At the seed germination stage, CRISPR/Cas9-created GmPYL gene knockout mutants, particularly gmpyl17/19 double mutants, are less susceptible to ABA than the wild type. RNA-Seq was used to investigate the differentially expressed genes related to the ABA response from germinated seedlings under diverse treatments using three biological replicates. The gmpyl17/19-1 double mutant was less susceptible to ABA during seed germination, and mutant plant height and branch number were higher than the wild type. Under ABA stress, the GO enrichment analysis showed that certain positive germination regulators were activated, which reduced ABA sensitivity and enhanced seed germination. This research gives a theoretical basis for a better understanding of the ABA signaling pathway and the participation of the key component at their molecular level, which helps enhance soybean abiotic stress tolerance. Furthermore, this research will aid breeders in regulating and improving soybean production and quality under various stress conditions.

Using Staphylococcus aureus Cas9 to Expand the Scope of Potential Gene Targets for Genome Editing in Soybean

The CRISPR (clustered regularly interspaced short palindromic repeat)/Cas (CRISPR-associated) is a revolutionary genome editing technology that has been used to achieve site-specific gene knock-out, large fragment deletion, or base editing in many plant species including soybean (Glycinemax). The Streptococcuspyogenes Cas9 (SpCas9) is widely used in plants at present, although there are some reports describing the application of CRISPR/Cpf1 in soybean. Therefore, the selection range of PAM (protospacer adjacent motif) sequences for soybean is currently limited to 5′-NGG-3′ (SpCas9) or 5′-TTTN-3′ (Cpf1), which in turn limits the number of genes that can be mutated. Another Cas9 enzyme from Staphylococcus aureus (SaCas9) recognizes the PAM sequence 5′-NNGRRT-3′ (where R represents A or G), which can provide a wider range of potential target sequences. In this study, we developed a CRISPR/SaCas9 system and used this tool to specifically induce targeted mutations at five target sites in the GmFT2a (Glyma.16G150700) and GmFT5a (Glyma.16G044100) genes in soybean hairy roots. We demonstrated that this tool can recognize the PAM sequences 5′-AAGGGT-3′, 5′-GGGGAT-3′, 5′-TTGAAT-3′, and 5′-TAGGGT-3′ in soybean, and it achieved mutation rates ranging from 34.5% to 73.3%. Our results show that we have established a highly efficient CRISPR/SaCas9 tool that is as suitable as SpCas9 for genome editing in soybean, and it will be useful for expanding the range of target sequences for genome editing.

SELF PRUNING 3C is a flowering repressor that modulates seed germination, root architecture, and drought responses

Allelic variation in the CETS (CENTRORADIALIS, TERMINAL FLOWER 1, SELF PRUNING) gene family controls agronomically important traits in many crops. CETS genes encode phosphatidylethanolamine-binding proteins that have a central role in the timing of flowering as florigenic and anti-florigenic signals. The great expansion of CETS genes in many species suggests that the functions of this family go beyond flowering induction and repression. Here, we characterized the tomato SELF PRUNING 3C (SP3C) gene, and show that besides acting as a flowering repressor it also regulates seed germination and modulates root architecture. We show that loss of SP3C function in CRISPR/Cas9-generated mutant lines increases root length and reduces root side branching relative to the wild type. Higher SP3C expression in transgenic lines promotes the opposite effects in roots, represses seed germination, and also improves tolerance to water stress in seedlings. These discoveries provide new insights into the role of SP paralogs in agronomically relevant traits, and support future exploration of the involvement of CETS genes in abiotic stress responses.

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.

Genomic Dissection and Diurnal Expression Analysis Reveal the Essential Roles of the PRR Gene Family in Geographical Adaptation of Soybean

Pseudo-response regulator (PRR) family members serve as key components of the core clock of the circadian clock, and play important roles in photoperiodic flowering, stress tolerance, growth, and the development of plants. In this study, 14 soybean PRR genes were identified, and classified into three groups according to phylogenetic analysis and structural characteristics. Real-time quantitative PCR analysis revealed that 13 GmPRRs exhibited obvious rhythmic expression under long-day (LD) and short-day (SD) conditions, and the expression of 12 GmPRRs was higher under LD in leaves. To evaluate the effects of natural variations in GmPRR alleles on soybean adaptation, we examined the sequences of GmPRRs among 207 varieties collected across China and the US, investigated the flowering phenotypes in six environments, and analyzed the geographical distributions of the major haplotypes. The results showed that a majority of non-synonymous mutations in the coding region were associated with flowering time, and we found that the nonsense mutations resulting in deletion of the CCT domain were related to early flowering. Haplotype analysis demonstrated that the haplotypes associated with early flowering were mostly distributed in Northeast China, while the haplotypes associated with late flowering were mostly cultivated in the lower latitudes of China. Our study of PRR family genes in soybean provides not only an important guide for characterizing the circadian clock-controlled flowering pathway but also a theoretical basis and opportunities to breed varieties with adaptation to specific regions and farming systems.

GmFT3a fine-tunes flowering time and improves adaptation of soybean to higher latitudes

Onset of flowering of plants is precisely controlled by extensive environmental factors and internal molecular networks, in which FLOWERING LOCUS T (FT) is a key flowering integrator. In soybean, a typical short-day plant, 11 FT homologues are found in its genome, of which several homologues are functionally diversified in flowering pathways and the others including GmFT3a are yet unknown. In the current study, we characterized GmFT3a, which is located on the same chromosome as the flowering promoters GmFT2a and GmFT5a. Overexpression of GmFT3a significantly promoted flowering of Arabidopsis under the inductive long-day (LD) photoperiod. GmFT3a over-expressed soybean also flowered earlier than the control under LD, but they were not significantly different under inductive short-day (SD) conditions, indicating that GmFT3a acts as a flowering promoter in the non-inductive photoperiod in soybean. Compared with other GmFT homologues, GmFT3a exhibited a slighter effect in flowering promotion than GmFT2a, GmFT5a and GmFT2b under LD conditions. GmFT3a promoted flowering by regulating the expression of downstream flowering-related genes and also affected the expression of other GmFTs. According to the re-sequencing data, the regional distributions of two major haplotypes in 176 soybean varieties were analyzed. The varieties with GmFT3a-Hap2 haplotype matured relatively early, and relative higher expression of GmFT3a was detected in early maturing varieties, implying that Hap2 variation may contribute to the adaptation of soybean to higher latitude regions by increasing expression level of genes in metabolism and signaling pathways. The early flowering germplasm generated by overexpression of GmFT3a has potential to be planted at higher latitudes where non-inductive long day is dominant in the growing season, and GmFT3a can be used to fine-tune soybean flowering and maturity time and improve the geographical adaptation.

Functional Redundancy of FLOWERING LOCUS T 3b in Soybean Flowering Time Regulation

Photoperiodic flowering is an important agronomic trait that determines adaptability and yield in soybean and is strongly influenced by FLOWERING LOCUS T (FT) genes. Due to the presence of multiple FT homologs in the genome, their functions in soybean are not fully understood. Here, we show that GmFT3b exhibits functional redundancy in regulating soybean photoperiodic flowering. Bioinformatic analysis revealed that GmFT3b is a typical floral inducer FT homolog and that the protein is localized to the nucleus. Moreover, GmFT3b expression was induced by photoperiod and circadian rhythm and was more responsive to long-day (LD) conditions. We generated a homozygous ft3b knockout and three GmFT3b-overexpressing soybean lines for evaluation under different photoperiods. There were no significant differences in flowering time between the wild-type, the GmFT3b overexpressors, and the ft3b knockouts under natural long-day, short-day, or LD conditions. Although the downstream flowering-related genes GmFUL1 (a, b), GmAP1d, and GmLFY1 were slightly down-regulated in ft3b plants, the floral inducers GmFT5a and GmFT5b were highly expressed, indicating potential compensation for the loss of GmFT3b. We suggest that GmFT3b acts redundantly in flowering time regulation and may be compensated by other FT homologs in soybean.

Progress in soybean functional genomics over the past decade

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

Progress in Soybean Genetic Transformation Over the Last Decade

Soybean is one of the important food, feed, and biofuel crops in the world. Soybean genome modification by genetic transformation has been carried out for trait improvement for more than 4 decades. However, compared to other major crops such as rice, soybean is still recalcitrant to genetic transformation, and transgenic soybean production has been hampered by limitations such as low transformation efficiency and genotype specificity, and prolonged and tedious protocols. The primary goal in soybean transformation over the last decade is to achieve high efficiency and genotype flexibility. Soybean transformation has been improved by modifying tissue culture conditions such as selection of explant types, adjustment of culture medium components and choice of selection reagents, as well as better understanding the transformation mechanisms of specific approaches such as Agrobacterium infection. Transgenesis-based breeding of soybean varieties with new traits is now possible by development of improved protocols. In this review, we summarize the developments in soybean genetic transformation to date, especially focusing on the progress made using Agrobacterium-mediated methods and biolistic methods over the past decade. We also discuss current challenges and future directions.

Enhancing crop diversity for food security in the face of climate uncertainty

Global agriculture is dominated by a handful of species that currently supply a huge proportion of our food and feed. It additionally faces the massive challenge of providing food for 10 billion people by 2050, despite increasing environmental deterioration. One way to better plan production in the face of current and continuing climate change is to better understand how our domestication of these crops included their adaptation to environments that were highly distinct from those of their centre of origin. There are many prominent examples of this, including the development of temperate Zea mays (maize) and the alteration of day-length requirements in Solanum tuberosum (potato). Despite the pre-eminence of some 15 crops, more than 50 000 species are edible, with 7000 of these considered semi-cultivated. Opportunities afforded by next-generation sequencing technologies alongside other methods, including metabolomics and high-throughput phenotyping, are starting to contribute to a better characterization of a handful of these species. Moreover, the first examples of de novo domestication have appeared, whereby key target genes are modified in a wild species in order to confer predictable traits of agronomic value. Here, we review the scale of the challenge, drawing extensively on the characterization of past agriculture to suggest informed strategies upon which the breeding of future climate-resilient crops can be based.

Molecular Genetic Understanding of Photoperiodic Regulation of Flowering Time in Arabidopsis and Soybean

The developmental switch from a vegetative phase to reproduction (flowering) is essential for reproduction success in flowering plants, and the timing of the floral transition is regulated by various environmental factors, among which seasonal day-length changes play a critical role to induce flowering at a season favorable for seed production. The photoperiod pathways are well known to regulate flowering time in diverse plants. Here, we summarize recent progresses on molecular mechanisms underlying the photoperiod control of flowering in the long-day plant Arabidopsis as well as the short-day plant soybean; furthermore, the conservation and diversification of photoperiodic regulation of flowering in these two species are discussed.

Overcoming the genetic compensation response of soybean florigens to improve adaptation and yield at low latitudes

The classical soybean (Glycine max) trait long juvenile (LJ) is essentially a reduction in sensitivity to short-day (SD) conditions for induction and completion of flowering, and has been introduced into soybean cultivars to improve yield in tropical environments. However, only one locus, J, is known to confer LJ in low-latitude varieties. Here, we defined two quantitative trait loci contributing to the LJ trait, LJ16.1 and LJ16.2, and identified them as the florigen (FT) homologs FT2a and FT5a, respectively. The two selected florigen variations both delay flowering time under SD conditions by repressing the floral meristem identity gene GmAPETALA1. Single mutants have a relatively subtle effect on flowering time and displayed a substantial genetic compensation response, but this was absent in ft2a ft5a double mutants, which showed an enhanced LJ phenotype that translated to higher yields under SD conditions. A survey of sequence diversity suggests that FT2a and FT5a variants have diverse origins and have played distinct roles as soybean spread to lower latitudes. Our results show that integration of variants in the florigen genes offers a strategy for customizing flowering time to adjust adaptation and improve crop productivity in tropical regions.

Speed-Breeding System in Soybean: Integrating Off-Site Generation Advancement, Fresh Seeding, and Marker-Assisted Selection

Speed breeding by artificial control of photothermal conditions facilitates generation advancement but was limited in scale and cost. In this study, we demonstrated a cost-saving off-site summer nursery pattern, taking full advantage of shorter daylength and higher temperature with lower latitude compared to the origin of the soybean cultivars used in the study. This substantially reduced the generation cycles under totally natural conditions. Using this approach, two generations of soybean cultivars from Northeastern Spring Planting Region (NE) and Yellow-Huai-Hai Valleys Summer Planting Region (YHH) were successfully obtained in Beijing and Hainan, respectively, compared to one generation in origin. Fresh-seeding method was also used to further shorten the generation duration by 7-10 days, thereby allowing at least four generations per year. Using DNA markers to define haplotypes of maturity genes E1-E4, we proposed a model to predict the optimum adaptation region of the advanced generation lines. Taken together, we present a speed-breeding methodology combining off-site nursery, fresh-seeding method, and marker-assisted selection, aimed at accelerating soybean improvement.

Cotyledons facilitate the adaptation of early-maturing soybean varieties to high-latitude long-day environments

Soybean (Glycine max), a typical short-day plant (SDP) domesticated in temperate regions, has expanded to high latitudes where daylengths are long from soybean emergence to bloom, but rapidly decrease from seed filling to maturity. Cotyledons are well known as the major storage organs in seeds, but it is unclear whether developing cotyledons store flowering substances at filling stage in SD for upcoming seedlings, or instead respond to photoperiod for floral induction after emergence of matured seeds in long-day (LD). Here, we report that cotyledons accelerate flowering of early-maturing varieties not resulting from stored floral stimuli but by perceiving photoperiod after emergence. We found that light signal is indispensable to activate cotyledons for floral induction, and flowering promoting gene GmFT2a is required for cotyledon-dependent floral induction via upregulation of floral identity gene GmAP1. Interestingly, cotyledons are competent to support the entire life cycle of a cotyledon-only plant to produce seeds, underlying a new photoperiod study system in soybean and other dicots. Taken together, these results demonstrate a substantial role for cotyledons in flowering process, whereby we propose a 'cotyledon-based self-reliance' model highlighting floral induction from emergence as a key ecological adaptation for rapid flowering of SDPs grown in LD environments at high latitudes.

Transcriptomal dissection of soybean circadian rhythmicity in two geographically, phenotypically and genetically distinct cultivars

BACKGROUND

In soybean, some circadian clock genes have been identified as loci for maturity traits. However, the effects of these genes on soybean circadian rhythmicity and their impacts on maturity are unclear.

RESULTS

We used two geographically, phenotypically and genetically distinct cultivars, conventional juvenile Zhonghuang 24 (with functional J/GmELF3a, a homolog of the circadian clock indispensable component EARLY FLOWERING 3) and long juvenile Huaxia 3 (with dysfunctional j/Gmelf3a) to dissect the soybean circadian clock with time-series transcriptomal RNA-Seq analysis of unifoliate leaves on a day scale. The results showed that several known circadian clock components, including RVE1, GI, LUX and TOC1, phase differently in soybean than in Arabidopsis, demonstrating that the soybean circadian clock is obviously different from the canonical model in Arabidopsis. In contrast to the observation that ELF3 dysfunction results in clock arrhythmia in Arabidopsis, the circadian clock is conserved in soybean regardless of the functional status of J/GmELF3a. Soybean exhibits a circadian rhythmicity in both gene expression and alternative splicing. Genes can be grouped into six clusters, C1-C6, with different expression profiles. Many more genes are grouped into the night clusters (C4-C6) than in the day cluster (C2), showing that night is essential for gene expression and regulation. Moreover, soybean chromosomes are activated with a circadian rhythmicity, indicating that high-order chromosome structure might impact circadian rhythmicity. Interestingly, night time points were clustered in one group, while day time points were separated into two groups, morning and afternoon, demonstrating that morning and afternoon are representative of different environments for soybean growth and development. However, no genes were consistently differentially expressed over different time-points, indicating that it is necessary to perform a circadian rhythmicity analysis to more thoroughly dissect the function of a gene. Moreover, the analysis of the circadian rhythmicity of the GmFT family showed that GmELF3a might phase- and amplitude-modulate the GmFT family to regulate the juvenility and maturity traits of soybean.

CONCLUSIONS

These results and the resultant RNA-seq data should be helpful in understanding the soybean circadian clock and elucidating the connection between the circadian clock and soybean maturity.

Quantitative Control of Early Flowering in White Lupin (Lupinus albus L.)

White lupin (Lupinus albus L.) is a pulse annual plant cultivated from the tropics to temperate regions for its high-protein grain as well as a cover crop or green manure. Wild populations are typically late flowering and have high vernalization requirements. Nevertheless, some early flowering and thermoneutral accessions were found in the Mediterranean basin. Recently, quantitative trait loci (QTLs) explaining flowering time variance were identified in bi-parental population mapping, however, phenotypic and genotypic diversity in the world collection has not been addressed yet. In this study, a diverse set of white lupin accessions (n = 160) was phenotyped for time to flowering in a controlled environment and genotyped with PCR-based markers (n = 50) tagging major QTLs and selected homologs of photoperiod and vernalization pathway genes. This survey highlighted quantitative control of flowering time in white lupin, providing statistically significant associations for all major QTLs and numerous regulatory genes, including white lupin homologs of CONSTANS, FLOWERING LOCUS T, FY, MOTHER OF FT AND TFL1, PHYTOCHROME INTERACTING FACTOR 4, SKI-INTERACTING PROTEIN 1, and VERNALIZATION INDEPENDENCE 3. This revealed the complexity of flowering control in white lupin, dispersed among numerous loci localized on several chromosomes, provided economic justification for future genome-wide association studies or genomic selection rather than relying on simple marker-assisted selection.

Molecular characterization of a soybean FT homologue, GmFT7

Soybean (Glycine max) is a vital oilseed legume crop that provides protein and oil for humans and feedstock for animals. Flowering is a prerequisite for seed production. Floral transition, from vegetative to reproductive stage, in a plant, is regulated by environmental (light, temperature) and endogenous factors. In Arabidopsis, Flowering Locus T (FT) protein is shown to be a mobile signal that moves from leaf to shoot apical meristem to induce flowering. However, FTs role in soybean is not fully resolved due to the presence of multiple (ten) homologs in the genome. Two of the ten FT homologs (GmFT2a and GmFT5a) have a role in the floral transition while GmFT1a and GmFT4 suppress soybean flowering. Recent deep sequencing data revealed that six FT homologs are expressed in shoot apical meristem and leaves during floral transition. One FT homolog, GmFT7 showed strong expression during soybean floral transition. Though bioinformatic analyses revealed that GmFT7 had high similarity with GmFT2a, ectopic GmFT7 expression in Arabidopsis could not promote flowering or rescue the late-flowering phenotype of Arabidopsis ft-10 mutant.

Soybean AP1 homologs control flowering time and plant height

Flowering time and plant height are key agronomic traits that directly affect soybean (Glycine max) yield. APETALA1 (AP1) functions as a class A gene in the ABCE model for floral organ development, helping to specify carpel, stamen, petal, and sepal identities. There are four AP1 homologs in soybean, all of which are mainly expressed in the shoot apex. Here, we used clustered regularly interspaced short palindromic repeats (CRISPR) - CRISPR-associated protein 9 technology to generate a homozygous quadruple mutant, gmap1, with loss-of-function mutations in all four GmAP1 genes. Under short-day (SD) conditions, the gmap1 quadruple mutant exhibited delayed flowering, changes in flower morphology, and increased node number and internode length, resulting in plants that were taller than the wild type. Conversely, overexpression of GmAP1a resulted in early flowering and reduced plant height compared to the wild type under SD conditions. The gmap1 mutant and the overexpression lines also exhibited altered expression of several genes related to flowering and gibberellic acid metabolism, thereby providing insight into the role of GmAP1 in the regulatory networks controlling flowering time and plant height in soybean. Increased node number is the trait with the most promise for enhancing soybean pod number and grain yield. Therefore, the mutant alleles of the four AP1 homologs described here will be invaluable for molecular breeding of improved soybean yield.

Photoperiod and Vernalization Control of Flowering-Related Genes: A Case Study of the Narrow-Leafed Lupin (Lupinus angustifolius L.)

Narrow-leafed lupin (Lupinus angustifolius L.) is a moderate-yielding legume crop known for its high grain protein content and contribution to soil improvement. It is cultivated under photoperiods ranging from 9 to 17 h, as a spring-sown (in colder locations) or as an autumn-sown crop (in warmer regions). Wild populations require a prolonged cold period, called vernalization, to induce flowering. The key achievement of L. angustifolius domestication was the discovery of two natural mutations (named Ku and Jul) conferring vernalization independence. These mutations are overlapping deletion variants in the promoter of LanFTc1, a homolog of the Arabidopsis thaliana FLOWERING LOCUS T (FT) gene. The third deletion, named here as Pal, was recently found in primitive germplasm. In this study, we genotyped L. angustifolius germplasm that differs in domestication status and geographical origin for LanFTc1 alleles, which we then phenotyped to establish flowering time and vernalization responsiveness. The Ku and Jul lines were vernalization-independent and early flowering, wild (ku) lines were vernalization-dependent and late flowering, whereas the Pal line conferred intermediate phenotype. Three lines representing ku, Pal, and Ku alleles were subjected to gene expression surveys under 8- and 16-h photoperiods. FT homologs (LanFTa1, LanFTa2, LanFTc1, and LanFTc2) and some genes selected by recent expression quantitative trait loci mapping were analyzed. Expression profiles of LanFTc1 and LanAGL8 (AGAMOUS-like 8) matched observed differences in flowering time between genotypes, highlighted by high induction after vernalization in the ku line. Moreover, these genes revealed altered circadian clock control in Pal line under short days. LanFD (FD) and LanCRLK1 (CALCIUM/CALMODULIN-REGULATED RECEPTOR-LIKE KINASE 1) were negatively responsive to vernalization in Ku and Pal lines but positively responsive or variable in ku, whereas LanUGT85A2 (UDP-GLUCOSYL TRANSFERASE 85A2) was significantly suppressed by vernalization in all lines. Such a pattern suggests the opposite regulation of these gene pairs in the vernalization pathway. LanCRLK1 and LanUGT85A2 are homologs of A. thaliana genes involved in the FLOWERING LOCUS C (FLC) vernalization pathway. Lupins, like many other legumes, do not have any FLC homologs. Therefore, candidate genes surveyed in this study, namely LanFTc1, LanAGL8, LanCRLK1, and LanUGT85A2, may constitute anchors for further elucidation of molecular components contributing to vernalization response in legumes.

Progresses, Challenges, and Prospects of Genome Editing in Soybean (Glycine max)

Soybean is grown worldwide for oil and protein source as food, feed and industrial raw material for biofuel. Steady increase in soybean production in the past century mainly attributes to genetic mediation including hybridization, mutagenesis and transgenesis. However, genetic resource limitation and intricate social issues in use of transgenic technology impede soybean improvement to meet rapid increases in global demand for soybean products. New approaches in genomics and development of site-specific nucleases (SSNs) based genome editing technologies have expanded soybean genetic variations in its germplasm and have potential to make precise modification of genes controlling the important agronomic traits in an elite background. ZFNs, TALENS and CRISPR/Cas9 have been adapted in soybean improvement for targeted deletions, additions, replacements and corrections in the genome. The availability of reference genome assembly and genomic resources increases feasibility in using current genome editing technologies and their new development. This review summarizes the status of genome editing in soybean improvement and future directions in this field.

Photoperiod and Vernalization Control of Flowering-Related Genes: A Case Study of the Narrow-Leafed Lupin (Lupinus angustifolius L.)

Narrow-leafed lupin (Lupinus angustifolius L.) is a moderate-yielding legume crop known for its high grain protein content and contribution to soil improvement. It is cultivated under photoperiods ranging from 9 to 17 h, as a spring-sown (in colder locations) or as an autumn-sown crop (in warmer regions). Wild populations require a prolonged cold period, called vernalization, to induce flowering. The key achievement of L. angustifolius domestication was the discovery of two natural mutations (named Ku and Jul) conferring vernalization independence. These mutations are overlapping deletion variants in the promoter of LanFTc1, a homolog of the Arabidopsis thaliana FLOWERING LOCUS T (FT) gene. The third deletion, named here as Pal, was recently found in primitive germplasm. In this study, we genotyped L. angustifolius germplasm that differs in domestication status and geographical origin for LanFTc1 alleles, which we then phenotyped to establish flowering time and vernalization responsiveness. The Ku and Jul lines were vernalization-independent and early flowering, wild (ku) lines were vernalization-dependent and late flowering, whereas the Pal line conferred intermediate phenotype. Three lines representing ku, Pal, and Ku alleles were subjected to gene expression surveys under 8- and 16-h photoperiods. FT homologs (LanFTa1, LanFTa2, LanFTc1, and LanFTc2) and some genes selected by recent expression quantitative trait loci mapping were analyzed. Expression profiles of LanFTc1 and LanAGL8 (AGAMOUS-like 8) matched observed differences in flowering time between genotypes, highlighted by high induction after vernalization in the ku line. Moreover, these genes revealed altered circadian clock control in Pal line under short days. LanFD (FD) and LanCRLK1 (CALCIUM/CALMODULIN-REGULATED RECEPTOR-LIKE KINASE 1) were negatively responsive to vernalization in Ku and Pal lines but positively responsive or variable in ku, whereas LanUGT85A2 (UDP-GLUCOSYL TRANSFERASE 85A2) was significantly suppressed by vernalization in all lines. Such a pattern suggests the opposite regulation of these gene pairs in the vernalization pathway. LanCRLK1 and LanUGT85A2 are homologs of A. thaliana genes involved in the FLOWERING LOCUS C (FLC) vernalization pathway. Lupins, like many other legumes, do not have any FLC homologs. Therefore, candidate genes surveyed in this study, namely LanFTc1, LanAGL8, LanCRLK1, and LanUGT85A2, may constitute anchors for further elucidation of molecular components contributing to vernalization response in legumes.