Anthracnose, caused by Colletotrichum truncatum, leads to large-scale reduction in quality and yield in soybean production. Limited information is available regarding the molecular mechanisms of resistance to anthracnose in soybean. We conducted a transcriptomic and targeted metabolomic analysis of pods from two soybean lines, “Zhechun No. 3” (ZC3) and ZC-2, in response to C. truncatum infection. Factors contributing to the enhanced resistance of ZC-2 to anthracnose compared with that of ZC3, included signal transduction (jasmonic acid, auxin, mitogen-activated protein kinase, and Ca²⁺ signaling), transcription factors (WRKY and bHLH), resistance genes (PTI1, RPP13, RGA2, RPS6, and ULP2B), pathogenesis-related genes (chitinase and lipid transfer protein), and terpenoid metabolism. Targeted metabolomic analysis revealed that terpenoid metabolism responded more promptly and more intensely to C. truncatum infection in ZC-2 than in ZC3. In vitro antifungal activity and resistance induction test confirmed that jasmonic acid, auxin signaling and terpenoids played important roles in soybean resistance to anthracnose. This research is the first study to explore the molecular mechanisms of soybean resistance to anthracnose. The findings are important for in-depth analysis of molecular resistance mechanisms, discovery of resistance genes, and to expedite the breeding of anthracnose-resistant soybean cultivars.

The vegetable soybean (Glycine max L. Merr.) plant is commonly consumed in Southeast Asian countries because of its nutritional value and desirable taste. A “pandan-like” aroma is an important value-added quality trait that is rarely found in commercial vegetable soybean varieties. In this study, three novel aromatic soybean cultivars with a fragrant volatile compound were isolated. We confirmed that the aroma of these cultivars is due to the potent volatile compound 2-acetyl-1-pyrroline (2AP) that was previously identified in soybean. A sequence comparison of GmBADH1/2 (encoding an aminoaldehyde dehydrogenase) between aromatic and non-aromatic soybean varieties revealed a mutation with 10 SNPs and an 11-nucleotide deletion in exon 1 of GmBADH2 in Quxian No. 1 and Xiangdou. Additionally, a 2-bp deletion was detected in exon 10 of GmBADH2 in ZK1754. The mutations resulted in a frame shift and the introduction of premature stop codons. Moreover, genetic analyses indicated that the aromatic trait in these three varieties was inherited according to a single recessive gene model. These results suggested that a mutated GmBADH2 may be responsible for the aroma of these three aromatic soybean cultivars. The expression and function of GmBADH2 in aromatic soybean seeds were confirmed by qRT-PCR and CRISPR/Cas9. A functional marker developed on the basis of the mutated GmBADH2 sequence in Quxian No. 1 and Xiangdou was validated in an F₂ population. A perfect association between the marker genotypes and aroma phenotypes implied that GmBADH2 is a major aroma-conferring gene. The results of this study are potentially useful for an in-depth analysis of the molecular basis of 2-AP formation in soybean and the marker-assisted breeding of aromatic vegetable soybean cultivars.

Homocysteine methyltransferase (HMT) converts homocysteine to methionine using S-methylmethionine (SMM) or S-adenosylmethionine (SAM) as methyl donors in organisms, playing an important role in supplying methionine for the growth and the development of plants. To better understand the functions of the HMT genes in plants, we conducted a wide evolution and expression analysis of these genes. Reconstruction of the phylogenetic relationship showed that the HMT gene family was divided into Class 1 and Class 2. In Class 1, HMTs were only found in seed plants, while Class 2 presented in all land plants, which hinted that the HMT genes might have diverged in seed plants. The analysis of gene structures and selection pressures showed that they were relatively conserved during evolution. However, type I functional divergence had been detected in the HMTs. Furthermore, the expression profiles of HMTs showed their distinct expression patterns in different tissues, in which some HMTs were widely expressed in various organs, whereas the others were highly expressed in some specific organs, such as seeds or leaves. Therefore, according to our results in the evolution, functional divergence, and expression, the HMT genes might have diverged during evolution. Further analysis in the expression patterns of AthHMTs with their methyl donors suggested that the diverged HMTs might be related to supply methionine for the development of plant seeds.

Seed germination is important to soybean (Glycine max) growth and development, ultimately affecting soybean yield. A lower seed field emergence has been the main hindrance for breeding soybeans low in phytate. Although this reduction could be overcome by additional breeding and selection, the mechanisms of seed germination in different low phytate mutants remain unknown. In this study, we performed a comparative transcript analysis of two low phytate soybean mutants (TW-1 and TW-1-M), which have the same mutation, a 2 bp deletion in GmMIPS1, but show a significant difference in seed field emergence, TW-1-M was higher than that of TW-1 . Numerous genes analyzed by RNA-Seq showed markedly different expression levels between TW-1-M and TW-1 mutants. Approximately 30,000–35,000 read-mapped genes and ~21000–25000 expressed genes were identified for each library. There were ~3900–9200 differentially expressed genes (DEGs) in each contrast library, the number of up-regulated genes was similar with down-regulated genes in the mutant TW-1and TW-1-M. Gene ontology functional categories of DEGs indicated that the ethylene-mediated signaling pathway, the abscisic acid-mediated signaling pathway, response to hormone, ethylene biosynthetic process, ethylene metabolic process, regulation of hormone levels, and oxidation-reduction process, regulation of flavonoid biosynthetic process and regulation of abscisic acid-activated signaling pathway had high correlations with seed germination. In total, 2457 DEGs involved in the above functional categories were identified. Twenty-two genes with 20 biological functions were the most highly up/down- regulated (absolute value Log2FC >5) in the high field emergence mutant TW-1-M and were related to metabolic or signaling pathways. Fifty-seven genes with 36 biological functions had the greatest expression abundance (FRPM >100) in germination-related pathways. Seed germination in the soybean low phytate mutants is a very complex process, which involves a series of physiological, morphological and transcriptional changes. Compared with TW-1, TW-1-M had a very different gene expression profile, which included genes related to plant hormones, antioxidation, anti-stress and energy metabolism processes. Our research provides a molecular basis for understanding germination mechanisms, and is also an important resource for the genetic analysis of germination in low phytate crops. Plant hormone- and antioxidation-related genes might strongly contribute to the high germination rate in the TW-1-M mutant.

Unsaturated fatty acids are the main components of vegetable oils. Fatty acid desaturase 2 (FAD2) catalyzes oleic acid (OA) into linoleic acid (LA) transformations, which are essential to the profile of FAs in seeds. To further understand the roles of FAD2s in the synthesis of oil, the evolution and biocatalysis of FAD2s were comprehensively analyzed. The evolution history of the FAD2 gene family showed that most of the FAD2 genes formed monophyletic clades except in eudicots. The FAD2 genes in some eudicots diverged into constitutive and seed-specific expression clades. Notably, the biocatalysis of seed-specific or -abundant expression FAD2s in soybean, perilla, rice, and spruce revealed that their catalytic activity was strongly correlated with the total oil content of their seeds in nature. Additionally, it was found that I and Y in site 143 of GmaFAD2-1 were strictly conserved in the seed-specific and constitutive expression clades of Fabaceae, respectively. Furthermore, the site-directed mutation demonstrated that I and Y are vital to improving and reducing the activity of GmaFAD2s. Therefore, the results indicate that the activity of FAD2s in seeds might be a reference to the total oil content of seeds, and site 143 might have been specifically evolved to be required for the activity of FAD2s in some expression-diverged eudicots, especially in legumes.

Vegetable soybean is a type of value-added specialty soybean, served as a fresh vegetable or snack in China. Due to the difference from other types, it is important to understand the genetic structure and diversity of vegetable soybean for further utilization in breeding programs. The four vegetable cultivars, Taiwan-75, Zhexiandou No. 8, Zhexian No. 9 and Zhexian No. 10 are popular soybean varieties planted in Zhejiang province, and have large pods and intermediate maturity. The clustering showed a close relationship of these four cultivars in simple sequence repeat analysis. To reveal the genome variation of vegetable soybean, these four improved lines were analyzed by whole-genome re-sequencing. The average sequencing depth was 7X and the coverage ratio of each cultivar was at least more than 94%. Compared with the reference genome, a large number of single-nucleotide polymorphisms, insertion/deletions and structure variations were identified with different chromosome distributions. The average heterozygosity rate of the single-nucleotide polymorphisms was 11.99% of these four cultivars. According to the enrichment analysis, there were 23,371 genes identified with putative modifications, and a total of 282 genes were related to carbohydrate metabolic processes. These results provide useful information for genetic research and future breeding, which can facilitate the selection procedures in vegetable soybean breeding.