The MADS-box transcription factors play fundamental roles in reproductive developmental control. Although the roles of many plant MADS-box proteins have been extensively studied, there are almost no functional studies of them in soybean, an important protein and oil crop in the world. In addition, the MADS-box protein orthologs may have species-specific functions. Controlling male fertility is an important goal in plant hybrid breeding but is difficult in some crops like soybean. The morphological structure of soybean flowers prevents the cross-pollination. Understanding the molecular mechanisms for floral development will aid in engineering new sterile materials that could be applied in hybrid breeding programs 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.

The MADS-box transcription factors are an ancient family of genes that regulate numerous physiological and biochemical processes in plants and facilitate the development of floral organs. However, the functions of most of these transcription factors in soybean remain unknown. In this work, a MADS-box gene, GmAGL1, was overexpressed in soybean. Phenotypic analysis showed that GmAGL1 overexpression not only resulted in early maturation but also promoted flowering and affected petal development. Furthermore, the GmAGL1 was much more effective at promoting flowering under long-day conditions than under short-day conditions. Transcriptome sequencing analysis showed that before flowering, the photoperiod pathway photoreceptor CRY2 and several circadian rhythm genes, such as SPA1, were significantly down-regulated, while some other flowering-promoting circadian genes, such as GI and LHY, and downstream genes related to flower development, such as FT, LEAFY, SEP1, SEP3, FUL, and AP1, were up-regulated compared with the control. Other genes related to the flowering pathway were not noticeably affected. The findings reported herein indicate that GmAGL1 may promote flowering mainly through the photoperiod pathway. Interestingly, while overexpression of GmAGL1 promoted plant maturity, no reduction in seed production or oil and protein contents was observed.

Pod dehiscence (shattering) is the main cause of serious yield loss during the soybean mechanical harvesting process. A better understanding of the genetic architecture and molecular mechanisms of pod dehiscence is of great significance for soybean breeding. In this study, genome-wide association analysis (GWAS) with NJAU 355K SoySNP array was performed to detect single nucleotide polymorphisms (SNPs) associated with pod dehiscence in an association panel containing 211 accessions across five environments. A total of 163 SNPs were identified as significantly associated with pod dehiscence. Among these markers, 136 SNPs identified on chromosome 16 were located in the known QTL qPDH1. One, one, three, eleven, three, one, three, three and one SNPs were distributed on chromosome 1, 4, 6, 8, 9, 11, 17, 18 and 20, respectively. Favorable SNPs and six haplotypes were identified based on ten functional SNPs; among those Hap2 and Hap3 were considered as optimal haplotypes. In addition, based on GWAS results, the candidate gene Glyma09g06290 was identified. Quantitative real-time PCR (qRT-PCR) results and polymorphism analysis suggested that Glyma09g06290 might be involved in pod dehiscence. Furthermore, a derived cleaved amplified polymorphic sequences (dCAPS) marker for Glyma09g06290 was developed. Overall, the loci and genes identified in this study will be helpful in breeding soybean accessions resistant to pod dehiscence.