The Myb73-GDPD2-GA2ox1 transcriptional regulatory module confers phosphate deficiency tolerance in soybean.

Phosphorus is indispensable in agricultural production. An increasing food supply requires more efficient use of phosphate due to limited phosphate resources. However, how crops regulate phosphate efficiency remains largely unknown. Here, we identified a major quantitative trait locus, qPE19, that controls 7 low-phosphate (LP)-related traits in soybean (Glycine max) through linkage mapping and genome-wide association studies. We identified the gene responsible for qPE19 as GLYCEROPHOSPHORYL DIESTER PHOSPHODIESTERASE2 (GmGDPD2), and haplotype 5 represents the optimal allele favoring LP tolerance. Overexpression of GmGDPD2 significantly affects hormone signaling and improves root architecture, phosphate efficiency and yield-related traits; conversely, CRISPR/Cas9-edited plants show decreases in these traits. GmMyb73 negatively regulates GmGDPD2 by directly binding to its promoter; thus, GmMyb73 negatively regulates LP tolerance. GmGDPD2 physically interacts with GA 2-oxidase 1 (GmGA2ox1) in the plasma membrane, and overexpressing GmGA2ox1 enhances LP-associated traits, similar to GmGDPD2 overexpression. Analysis of double mutants for GmGDPD2 and GmGA2ox1 demonstrated that GmGDPD2 regulates LP tolerance likely by influencing auxin and gibberellin dose-associated cell division in the root. These results reveal a regulatory module that plays a major role in regulating LP tolerance in soybeans and is expected to be utilized to develop phosphate-efficient varieties to enhance soybean production, particularly in phosphate-deficient soils.

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

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

GmEID1 modulates light signaling through the Evening Complex to control flowering time and yield in soybean.

Soybean (Glycine max) morphogenesis and flowering time are accurately regulated by photoperiod, which determine the yield potential and limit soybean cultivars to a narrow latitudinal range. The E3 and E4 genes, which encode phytochrome A photoreceptors in soybean, promote the expression of the legume-specific flowering repressor E1 to delay floral transition under long-day (LD) conditions. However, the underlying molecular mechanism remains unclear. Here, we show that the diurnal expression pattern of GmEID1 is opposite to that of E1 and targeted mutations in the GmEID1 gene delay soybean flowering regardless of daylength. GmEID1 interacts with J, a key component of circadian Evening Complex (EC), to inhibit E1 transcription. Photoactivated E3/E4 interacts with GmEID1 to inhibit GmEID1-J interaction, promoting J degradation resulting in a negative correlation between daylength and the level of J protein. Notably, targeted mutations in GmEID1 improved soybean adaptability by enhancing yield per plant up to 55.3% compared to WT in field trials performed in a broad latitudinal span of more than 24°. Together, this study reveals a unique mechanism in which E3/E4-GmEID1-EC module controls flowering time and provides an effective strategy to improve soybean adaptability and production for molecular breeding.

A retrotransposon insertion in the Mao1 promoter results in erect pubescence and higher yield in soybean.

Adaptive changes in crops contribute to the diversity of agronomic traits, which directly or indirectly affect yield. The change of pubescence form from appressed to erect is a notable feature during soybean domestication. However, the biological significance and regulatory mechanism underlying this transformation remain largely unknown. Here, we identified a major-effect locus, PUBESCENCE FORM 1 (PF1), the upstream region of Mao1, that regulates pubescence form in soybean. The insertion of a Ty3/Gypsy retrotransposon in PF1 can recruit the transcription factor GAGA-binding protein to a GA-rich region, which up-regulates Mao1 expression, underpinning soybean pubescence evolution. Interestingly, the proportion of improved cultivars with erect pubescence increases gradually with increasing latitude, and erect-pubescence cultivars have a higher yield possibly through a higher photosynthetic rate and photosynthetic stability. These findings open an avenue for molecular breeding through either natural introgression or genome editing toward yield improvement and productivity.

SoyDNGP: a web-accessible deep learning framework for genomic prediction in soybean breeding.

Soybean is a globally significant crop, playing a vital role in human nutrition and agriculture. Its complex genetic structure and wide trait variation, however, pose challenges for breeders and researchers aiming to optimize its yield and quality. Addressing this biological complexity requires innovative and accurate tools for trait prediction. In response to this challenge, we have developed SoyDNGP, a deep learning-based model that offers significant advancements in the field of soybean trait prediction. Compared to existing methods, such as DeepGS and DNNGP, SoyDNGP boasts a distinct advantage due to its minimal increase in parameter volume and superior predictive accuracy. Through rigorous performance comparison, including prediction accuracy and model complexity, SoyDNGP represents improved performance to its counterparts. Furthermore, it effectively predicted complex traits with remarkable precision, demonstrating robust performance across different sample sizes and trait complexities. We also tested the versatility of SoyDNGP across multiple crop species, including cotton, maize, rice and tomato. Our results showed its consistent and comparable performance, emphasizing SoyDNGP's potential as a versatile tool for genomic prediction across a broad range of crops. To enhance its accessibility to users without extensive programming experience, we designed a user-friendly web server, available at http://xtlab.hzau.edu.cn/SoyDNGP. The server provides two features: 'Trait Lookup', offering users the ability to access pre-existing trait predictions for over 500 soybean accessions, and 'Trait Prediction', allowing for the upload of VCF files for trait estimation. By providing a high-performing, accessible tool for trait prediction, SoyDNGP opens up new possibilities in the quest for optimized soybean breeding.

Novel and multifaceted regulations of photoperiodic flowering by phytochrome A in soybean.

Photoperiod is an important environmental cue. Plants can distinguish the seasons and flower at the right time through sensing the photoperiod. Soybean is a sensitive short-day crop, and the timing of flowering varies greatly at different latitudes, thus affecting yields. Soybean cultivars in high latitudes adapt to the long day by the impairment of two phytochrome genes, PHYA3 and PHYA2, and the legume-specific flowering suppressor, E1. However, the regulating mechanism underlying phyA and E1 in soybean remains largely unknown. Here, we classified the regulation of the E1 family by phyA2 and phyA3 at the transcriptional and posttranscriptional levels, revealing that phyA2 and phyA3 regulate E1 by directly binding to LUX proteins, the critical component of the evening complex, to regulate the stability of LUX proteins. In addition, phyA2 and phyA3 can also directly associate with E1 and its homologs to stabilize the E1 proteins. Therefore, phyA homologs control the core flowering suppressor E1 at both the transcriptional and posttranscriptional levels, to double ensure the E1 activity. Thus, our results disclose a photoperiod flowering mechanism in plants by which the phytochrome A regulates LUX and E1 activity.

Natural variation in Fatty Acid 9 is a determinant of fatty acid and protein content.

Soybean is one of the most economically important crops worldwide and an important source of unsaturated fatty acids and protein for the human diet. Consumer demand for healthy fats and oils is increasing, and the global demand for vegetable oil is expected to double by 2050. Identification of key genes that regulate seed fatty acid content can facilitate molecular breeding of high-quality soybean varieties with enhanced fatty acid profiles. Here, we analysed the genetic architecture underlying variations in soybean seed fatty acid content using 547 accessions, including mainly landraces and cultivars from northeastern China. Through fatty acid profiling, genome re-sequencing, population genomics analyses, and GWAS, we identified a SEIPIN homologue at the FA9 locus as an important contributor to seed fatty acid content. Transgenic and multiomics analyses confirmed that FA9 was a key regulator of seed fatty acid content with pleiotropic effects on seed protein and seed size. We identified two major FA9 haplotypes in 1295 resequenced soybean accessions and assessed their phenotypic effects in a field planting of 424 accessions. Soybean accessions carrying FA9H2 had significantly higher total fatty acid contents and lower protein contents than those carrying FA9H1 . FA9H2 was absent in wild soybeans but present in 13% of landraces and 26% of cultivars, suggesting that it may have been selected during soybean post-domestication improvement. FA9 therefore represents a useful genetic resource for molecular breeding of high-quality soybean varieties with specific seed storage profiles.

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

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

A novel miR160a-GmARF16-GmMYC2 module determines soybean salt tolerance and adaptation.

Salt stress is a major challenge that has a negative impact on soybean growth and productivity. Therefore, it is important to understand the regulatory mechanism of salt response to ensure soybean yield under such conditions. In this study, we identified and characterized a miR160a-GmARF16-GmMYC2 module and its regulation during the salt-stress response in soybean. miR160a promotes salt tolerance by cleaving GmARF16 transcripts, members of the Auxin Response Factor (ARF) family, which negatively regulates salt tolerance. In turn, GmARF16 activates GmMYC2, encoding a bHLH transcription factor that reduces salinity tolerance by down-regulating proline biosynthesis. Genomic analysis among wild and cultivated soybean accessions identified four distinct GmARF16 haplotypes. Among them, the GmARF16H3 haplotype is preferentially enriched in localities with relatively saline soils, suggesting GmARF16H3 was artificially selected to improve salt tolerance. Our findings therefore provide insights into the molecular mechanisms underlying salt response in soybean and provide valuable genetic targets for the molecular breeding of salt tolerance.

EDS1-interacting J protein 1 is an essential negative regulator of plant innate immunity in Arabidopsis.

Plants have evolved precise mechanisms to optimize immune responses against pathogens. ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) plays a vital role in plant innate immunity by regulating basal resistance and effector-triggered immunity. Nucleocytoplasmic trafficking of EDS1 is required for resistance reinforcement, but the molecular mechanism remains elusive. Here, we show that EDS1-INTERACTING J PROTEIN1 (EIJ1), which acts as a DnaJ protein-like chaperone in response to pathogen infection, functions as an essential negative regulator of plant immunity by interacting with EDS1. The loss-of-function mutation of EIJ1 did not affect plant growth but significantly enhanced pathogen resistance. Upon pathogen infection, EIJ1 relocalized from the chloroplast to the cytoplasm, where it interacted with EDS1, thereby restricting pathogen-triggered trafficking of EDS1 to the nucleus and compromising resistance at an early infection stage. During disease development, EIJ1 was gradually degraded, allowing the nuclear accumulation of EDS1 for transcriptional resistance reinforcement. The avirulent strain Pst DC3000 (AvrRps4) abolished the repressive action of EIJ1 by rapidly inducing its degradation in the effector-triggered immunity response. Thus, our findings show that EIJ1 is an essential EDS1-dependent negative regulator of innate plant immunity and provide a mechanistic understanding of how the nuclear versus cytoplasmic distribution of EDS1 is regulated during the immune response.

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

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

Genome-wide identification and expression analysis of the WNK kinase gene family in soybean.

WNK kinases are a unique class of serine/threonine protein kinases that lack a conserved catalytic lysine residue in the kinase domain, hence the name WNK (with no K, i.e., lysine). WNK kinases are involved in various physiological processes in plants, such as circadian rhythm, flowering time, and stress responses. In this study, we identified 26 WNK genes in soybean and analyzed their phylogenetic relationships, gene structures, chromosomal distribution, cis-regulatory elements, expression patterns, and conserved protein motifs. The soybean WNK genes were unevenly distributed on 15 chromosomes and underwent 21 segmental duplication events during evolution. We detected 14 types of cis-regulatory elements in the promoters of the WNK genes, indicating their potential involvement in different signaling pathways. The transcriptome database revealed tissue-specific and salt stress-responsive expression of WNK genes in soybean, the second of which was confirmed by salt treatments and qRT-PCR analysis. We found that most WNK genes were significantly up-regulated by salt stress within 3 h in both roots and leaves, except for WNK5, which showed a distinct expression pattern. Our findings provide valuable insights into the molecular characteristics and evolutionary history of the soybean WNK gene family and lay a foundation for further analysis of WNK gene functions in soybean. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-024-01440-5.

A critical role of the soybean evening complex in the control of photoperiod sensitivity and adaptation.

Photoperiod sensitivity is a key factor in plant adaptation and crop production. In the short-day plant soybean, adaptation to low latitude environments is provided by mutations at the J locus, which confer extended flowering phase and thereby improve yield. The identity of J as an ortholog of Arabidopsis ELF3, a component of the circadian evening complex (EC), implies that orthologs of other EC components may have similar roles. Here we show that the two soybean homeologs of LUX ARRYTHMO interact with J to form a soybean EC. Characterization of mutants reveals that these genes are highly redundant in function but together are critical for flowering under short day, where the lux1 lux2 double mutant shows extremely late flowering and a massively extended flowering phase. This phenotype exceeds that of any soybean flowering mutant reported to date, and is strongly reminiscent of the "Maryland Mammoth" tobacco mutant that featured in the seminal 1920 study of plant photoperiodism by Garner and Allard [W. W. Garner, H. A. Allard, J. Agric. Res. 18, 553-606 (1920)]. We further demonstrate that the J-LUX complex suppresses transcription of the key flowering repressor E1 and its two homologs via LUX binding sites in their promoters. These results indicate that the EC-E1 interaction has a central role in soybean photoperiod sensitivity, a phenomenon also first described by Garner and Allard. EC and E1 family genes may therefore constitute key targets for customized breeding of soybean varieties with precise flowering time adaptation, either by introgression of natural variation or generation of new mutants by gene editing.

The AP2/ERF transcription factor TOE4b regulates photoperiodic flowering and grain yield per plant in soybean.

Photoperiod-mediated flowering determines the phenological adaptability of crops including soybean (Glycine max). A genome-wide association study (GWAS) identified a new flowering time locus, Time of flowering 13 (Tof13), which defined a gene encoding an AP2/ERF transcription factor. This new transcription factor, which we named TOE4b, is localized in the nucleus. TOE4b has been selected for soybean latitude adaptability. The existing natural variant TOE4bH4 was rare in wild soybean accessions but occurred more frequently in landraces and cultivars. Notably, TOE4bH4 improved high-latitude adaptation of soybean to some extent. The gene-edited TOE4b knockout mutant exhibited earlier flowering, conversely, TOE4b overexpression delayed flowering time. TOE4b is directly bound to the promoters and gene bodies of the key flowering integration factor genes FT2a and FT5a to inhibit their transcription. Importantly, TOE4b overexpression lines in field trials not only showed late flowering but also altered plant architecture, including shorter internode length, more internodes, more branches and pod number per plant, and finally boosted grain yield per plant by 60% in Guangzhou and 87% in Shijiazhuang. Our findings therefore identified TOE4b as a pleiotropic gene to increase yield potential per plant in soybean, and these results provide a promising option for breeding a soybean variety with an idealized plant architecture that promotes high yields.

Natural variation of MS2 confers male fertility and drives hybrid breeding in soybean.

A complete and genetically stable male sterile line with high outcrossing rate is a prerequisite for the development of commercial hybrid soybean. It was reported in the last century that the soybean male sterile ms2 mutant has the highest record with seed set. Here we report the cloning and characterization of the MS2 gene in soybean, which encodes a protein that is specifically expressed in the anther. MS2 functions in the tapetum and microspore by directly regulating genes involved in the biosynthesis of secondary metabolites and the lipid metabolism, which is essential for the formation of microspore cell wall. Through comparison of the field performance with the widely used male sterile mutants in the same genetic background, we demonstrated that the ms2 mutant conducts the best in outcrossing rate and makes it an ideal tool in building a cost-effective hybrid system for soybean.

Global dissection of the recombination landscape in soybean using a high-density 600K SoySNP array.

Recombination is crucial for crop breeding because it can break linkage drag and generate novel allele combinations. However, the high-resolution recombination landscape and its driving forces in soybean are largely unknown. Here, we constructed eight recombinant inbred line (RIL) populations and genotyped individual lines using the high-density 600K SoySNP array, which yielded a high-resolution recombination map with 5636 recombination sites at a resolution of 1.37 kb. The recombination rate was negatively correlated with transposable element density and GC content but positively correlated with gene density. Interestingly, we found that meiotic recombination was enriched at the promoters of active genes. Further investigations revealed that chromatin accessibility and active epigenetic modifications promoted recombination. Our findings provide important insights into the control of homologous recombination and thus will increase our ability to accelerate soybean breeding by manipulating meiotic recombination rate.

Diverse flowering responses subjecting to ambient high temperature in soybean under short-day conditions.

Flowering time is one of important agronomic traits determining the crop yield and affected by high temperature. When facing high ambient temperature, plants often initiate early flowering as an adaptive strategy to escape the stress and ensure successful reproduction. However, here we find opposing ways in the short-day crop soybean to respond to different levels of high temperatures, in which flowering accelerates when temperature changes from 25 to 30 °C, but delays when temperature reaches 35 °C under short day. phyA-E1, possibly photoperiodic pathway, is crucial for 35 °C-mediated late flowering, however, does not contribute to promoting flowering at 30 °C. 30 °C-induced up-regulation of FT2a and FT5a leads to early flowering, independent of E1. Therefore, distinct responsive mechanisms are adopted by soybean when facing different levels of high temperatures for successful flowering and reproduction.

Natural variation of FKF1 controls flowering and adaptation during soybean domestication and improvement.

Soybean (Glycine max) is a major source of protein and edible oil world-wide and is cultivated in a wide range of latitudes. However, it is extremely sensitive to photoperiod, which influences flowering time, maturity, and yield, and severely limits soybean latitude adaptation. In this study, a genome-wide association study (GWAS) identified a novel locus in accessions harboring the E1 allele, called Time of flowering 8 (Tof8), which promotes flowering and enhances adaptation to high latitude in cultivated soybean. Gene functional analyses showed that Tof8 is an ortholog of Arabidopsis FKF1. We identified two FKF1 homologs in the soybean genome. Both FKF1 homologs are genetically dependent on E1 by binding to E1 promoter to activate E1 transcription, thus repressing FLOWERING LOCUS T 2a (FT2a) and FT5a transcription, which modulate flowering and maturity through the E1 pathway. We also demonstrate that the natural allele FKF1bH3 facilitated adaptation of soybean to high-latitude environments and was selected during domestication and improvement, leading to its rapid expansion in cultivated soybean. These findings provide novel insights into the roles of FKF1 in controlling flowering time and maturity in soybean and offer new means to fine-tune adaptation to high latitudes and increase grain yield.

CALCIUM-DEPENDENT PROTEIN KINASE38 regulates flowering time and common cutworm resistance in soybean.

Photoperiod-sensitive plants such as soybean (Glycine max) often face threats from herbivorous insects throughout their whole growth period and especially during flowering; however, little is known about the relationship between plant flowering and insect resistance. Here, we used gene editing, multiple omics, genetic diversity and evolutionary analyses to confirm that the calcium-dependent protein kinase GmCDPK38 plays a dual role in coordinating flowering time regulation and insect resistance of soybean. Haplotype 2 (Hap2)-containing soybeans flowered later and were more resistant to the common cutworm (Spodoptera litura Fabricius) than those of Hap3. gmcdpk38 mutants with Hap3 knocked out exhibited similar flowering and resistance phenotypes as Hap2. Knocking out GmCDPK38 altered numerous flowering- and resistance-related phosphorylated proteins, genes, and metabolites. For example, the S-adenosylmethionine synthase GmSAMS1 was post-translationally upregulated in the gmcdpk38 mutants. GmCDPK38 has abundant genetic diversity in wild soybeans and was likely selected during soybean domestication. We found that Hap2 was mostly distributed at low latitudes and had a higher frequency in cultivars than in wild soybeans, while Hap3 was widely selected at high latitudes. Overall, our results elucidated that the two distinct traits (flowering time and insect resistance) are mediated by GmCDPK38.

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.