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.

Soybean is one of the important economic crops, which supplies a great deal of vegetable oil and proteins for human. The content of nutrients in different soybean seeds is different, which is related to the expression of multiple genes, but the mechanisms are complicated and still largely uncertain. In this study, to reveal the possible causes of the nutrients difference in soybeans A7 (containing low oil and high protein) and A35 (containing high oil and low protein), RNA-seq technology was performed to compare and identify the potential differential expressed genes (DEGs) at different seed developmental stages. The results showed that DEGs mainly presented at the early stages of seeds development and more DEGs were up-regulated at the early stage than the late stages. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis showed that the DEGs have diverged in A7 and A35. In A7, the DEGs were mainly involved in cell cycle and stresses, while in A35 were the fatty acids and sugar metabolism. Specifically, when the DEGs contributing to oil and protein metabolic pathways were analyzed, the differences between A7 and A35 mainly presented in fatty acids metabolism and seeds storage proteins (SSPs) synthesis. Furthermore, the enzymes, fatty acid dehydrogenase 2, 3-ketoacyl-CoA synthase and 9S-lipoxygenase, in the synthesis and elongation pathways of fatty acids, were revealed probably to be involved in the oil content difference between A7 and A35, the SSPs content might be due to the transcription factors:Leafy Cotyledon 2and Abscisicacid-intensitive 3, while the sugar transporter, SWEET10a, might contribute to both oil and protein content differences. Finally, six DEGs were selected to analyze their expression using qRT-PCR, and the results were consistent with the RNA-seq results. Generally, the study provided a comprehensive and dynamic expression trends for the seed development processes, and uncovered the potential DEGs for the differences of oil in A7 and A35.

The members of the Old Yellow Enzyme (OYE) family are capable of catalyzing the asymmetric reduction of (E/Z)-citral to (R)-citronellal—a key intermediate in the synthesis of L-menthol. The applications of OYE-mediated biotransformation are usually hampered by its insufficient enantioselectivity and low activity. Here, the (R)-enantioselectivity of Old Yellow Enzyme from Saccharomyces cerevisiae CICC1060 (OYE2y) was enhanced through protein engineering. The single mutations of OYE2y revealed that the sites R330 and P76 could act as the enantioselectivity switch of OYE2y. Site-saturation mutagenesis was conducted to generate all possible replacements for the sites R330 and P76, yielding 17 and five variants with improved (R)-enantioselectivity in the (E/Z)-citral reduction, respectively. Among them, the variants R330H and P76C partly reversed the neral derived enantioselectivity from 32.66% e.e. (S) to 71.92% e.e. (R) and 37.50% e.e. (R), respectively. The docking analysis of OYE2y and its variants revealed that the substitutions R330H and P76C enabled neral to bind with a flipped orientation in the active site and thus reverse the enantioselectivity. Remarkably, the double substitutions of R330H/P76M, P76G/R330H, or P76S/R330H further improved (R)-enantioselectivity to >99% e.e. in the reduction of (E)-citral or (E/Z)-citral. The results demonstrated that it was feasible to alter the enantioselectivity of OYEs through engineering key residue distant from active sites, e.g., R330 in OYE2y.

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.

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.

Possessing aldehyde and carboxyl groups, glyoxylic acid and its ester derivatives serve as platform chemicals for the synthesis of vanillin, (R)-pantolactone, antibiotics or agrochemicals. Methyl glycolate is one of the by-products in the coal-to-glycol industry, and we attempted its value-added use through enzymatic oxidation of methyl glycolate to methyl glyoxylate. The cascade catalysis of glycolate oxidase from Spinacia oleracea (SoGOX), catalase from Helicobacter pylori (HpCAT) and hemoglobin from Vitreoscilla stercoraria (VsHGB) was firstly constructed, despite poor catalytic performance. To enable efficient oxidation of methyl glycolate, eight fusion enzymes of SoGOX, HpCAT and VsHGB were constructed by varying the orientation and the linker length. The fusion enzyme VsHGB-GSG-SoGOX-GGGGS-HpCAT was proved to be best, which reaction yield was 2.9 times higher than that of separated enzymes. The enzyme SoGOX was further subjected to directed evolution and site-saturation mutagenesis. The reaction yield of the resulting variant M267T/S362G was 1.9 times higher than that of the wild type. Then, the double substitution M267T/S362G was integrated with fusion expression to give the fusion enzyme VsHGB-GSG-SoGOX^mut-GGGGS-HpCAT, which crude enzyme was used as biocatalyst. The use of crude enzyme virtually eliminated side reactions and simplified the preparation of biocatalysts. Under the optimized conditions, the crude enzyme VsHGB-GSG-SoGOX^mut-GGGGS-HpCAT catalyzed the oxidation of 200 mM methyl glycolate for 6 h, giving a yield of 95.3%. The development of efficient fusion enzyme and the use of its crude enzyme paved the way for preparative scale application on enzymatic oxidation of methyl glycolate to methyl glyoxylate.