Characterization of the Key Aroma Compounds of Soybean Flavor in Fermented Soybeans with Bacillus subtilis BJ3-2 by Gene Knockout, Gas Chromatography–Olfactometry–Mass Spectrometry, and Aroma Addition Experiments

: Soybean flavor is considered to be essential for the aroma quality of fermented soybeans (FS) with Bacillus subtilis BJ3-2 (BJ3-2) at 37 ◦ C. However, the key aroma compounds of the soybean flavor must be further elucidated. In this study, two candidate genes ( sdaAA and katX ) of BJ3-2 involved in the control of soybean flavor production were screened using prior multi-omics data. FS samples with BJ3-2, BJ3-2 ∆ sdaAA , BJ3-2 ∆ katX , and BJ3-2 ∆ sdaAA ∆ katX were analyzed by quantitative descriptive sensory analysis (QDA), gas chromatography–olfactometry–mass spectrometry (GC-O-MS), relative odor activity values (ROAV), and aroma addition experiments. The QDA revealed that the aroma profile of the soybean flavor in FS consisted of “sweaty”, “smoky”, “beany”, “roasted”, and “sweet” attributes. A total of 20 aroma-active compounds were detected, and 13 of them with ROAV > 1 were identified as key aroma compounds. Moreover, aroma addition experiments were conducted to further confirm the key aroma compounds of soybean flavor. Among them, 2-methylbutyric acid, 2,3,5-trimethylpyrazine, and guaiacol contributed higher aroma intensity values and ROAV, resulting in “sweaty”, “roasted”, and “smoky” attributes of soybean flavor in FS, respectively, while 1-octen-3-ol was associated with the “beany” attribute. These findings provide novel insights into the aroma attributes of soybean flavor in FS and a new strategy for revealing the key aroma compounds in fermented foods.


Introduction
Fermented soybean foods are widely popular in East and Southeast Asian countries due to their unique flavor and nutritional value.Bacillus subtilis (B.subtilis) is a predominant functional bacterium in these foods [1].During the fermentation of soybeans by B. subtilis, the exogenous enzymes secreted by B. subtilis (mainly proteases and lipases) can decompose the proteins and lipids into amino acids, fatty acids, and organic acids, which then form volatile compounds [2,3].The soybean flavor is a special flavor formed by fermenting the raw material by microorganisms (namely Chi-flavor, a special flavor with an ammonia scent distinct from natto) [4].Previous research has reported that the characteristics of B. subtilis BJ3-2 fermented soybean differ at different temperatures.For example, fermented soybeans at 37 •  Chi-flavor-type Baijiu are popular and accepted worldwide due to their prominent soybean flavor [6,7], there are only a few reports on the study of soybean flavor.It is speculated that the synergistic effects of multiple metabolic pathways and multiple genes may produce soybean flavor [5].Meanwhile, the key aroma compounds contributing to soybean flavor in fermented soybeans are poorly known.
Systematic characterization of the aroma-active compounds therein is a prerequisite for improving the flavor quality of fermented soybeans [8].Food aroma compounds are usually studied using an integrated approach encompassing isolation, identification, and quantification of volatile compounds, a deep analysis of aroma compounds, and reconstitution the aroma using the key aroma compounds [9,10].Gas chromatography-olfactometry-mass spectrometry (GC-O-MS), a combination of gas chromatography-olfactometry (GC-O) and gas chromatography-mass spectrometry (GC-MS), separates the volatile compounds in a complex sample through gas chromatography, enables the determination of which of these volatiles are aromatically active through olfactometric detection, and allows the identification of the compounds eluting from the column through mass spectrometry [11].Thus far, the aroma compounds in several B. subtilis fermented soybean foods have been examined using GC-O-MS.Previous studies have demonstrated that pyrazines play a crucial role in the flavor of natto (a traditional Japanese food made from fermented soybeans), among which 2,3,5-trimethylpyrazine, 2,5-dimethylpyrazine, 2,3,5,6-tetramethylpyrazine, and 2,3,5-methyl-6-ethylpyrazine substantially contribute to natto flavor [12].Furthermore, Mathatheeranan et al. identified the dominant aroma compounds of thua nao (indigenous Thai fermented soybeans) fermented by commercial B. subtilis under controlled conditions including 2-methylbutyric acid, phenylethyl alcohol, benzaldehyde, 2,5-dimethylpyrazine, 2,6-dimethylpyrazine, and 2-methoxyphenol, which contribute to the fermented/sour, floral, almond-like, and nutty aromas of fermented soybeans [13].More recently, GC-O-MS combined with relative odor activity values (ROAV) were applied to examine the aroma compounds of sufu [14], Yongchuan Douchi [6], typical fermented soy foods [15], and so on.In these studies, e-nose analysis and quantitative description sensory analysis (QDA) were used to analyze the aroma profiles of different fermented soy products.Differences in the aroma profiles of the various products have been identified and associated with key aroma compounds.
In general, microbial metabolism is crucial for producing flavor in fermented foods [16].To elucidate the mechanism of soybean flavor and soy sauce-like flavor production by B. subtilis BJ3-2, Wu et al. measured and analyzed the transcriptome, proteome, and metabolome of B. subtilis BJ3-2 at three temperatures (37 • C, 45 • C, and 53 • C).The functions of differentially expressed genes (DEGs) have been analyzed by KEGG and GO enrichment analyses.These DEGs are enriched in the relevant metabolic pathways, including the glycine, serine, threonine, cysteine, and methionine metabolite pathways, among others [5].Meanwhile, amino acids are substrates for the formation of aroma compounds [17].The production of soybean flavor may be closely linked to these metabolic pathways.Nevertheless, the candidate genes that may regulate the production of soybean flavor need to be further investigated.
Therefore, this study aimed to (a) screen soybean flavor candidate genes from the analysis of multi-omics data; (b) obtain mutants with significantly altered soybean flavor by knocking out the soybean flavor candidate gene to better analyze the relationship between sensory attributes and aroma compounds of soybean flavor in fermented soybeans (FS); (c) characterize key aroma compounds of soybean flavor in FS using GC-O and ROAV; and (d) verify the correlation between key aroma compounds and sensory attributes by partial least squares regression (PLSR) analysis.In the present study, by employing multi-omics data, gene knockout, and the headspace solid-phase microextraction (HS-SPME) method, combined with GC-O-MS, ROAV, and aroma addition experiments, we have become the first to characterize the aroma of soybean flavor in FS samples with B. subtilis BJ3-2 and its different mutants.These data can help gain a deeper understanding of the key aroma compounds in the production of soybean flavor in fermented soybeans with B. subtilis.
They also provide a new strategy for the screening of soybean flavor candidate genes and the modification of fermentation flavor by genetic engineering.

Reverse Transcription Quantitative Real-Time PCR (RT-qPCR)
To validate the expression of sdaAA and katX in the multi-omics data, the total RNA of BJ3-2 was extracted, and the first strand of cDNA was synthesized using the FastKing method as a template (TIANGEN, Beijing, China).Then, the cDNA was detected by the CFX96 Touch PCR instrument (Bio-Rad, Hercules, CA, USA).The reaction system and the settings applied followed a previously described approach [18].16S rRNA was used as the internal reference gene.The primers used are listed in Table S1.

Construction and Transformation of Homologous Recombinant Knockout Vectors
HLarm is the homologous left arm of the target gene; HRarm is the homologous right arm; cm is the chloramphenicol gene; and erm is the erythromycin gene.DEDP is a double-exchange detection primer (Table S1) used to assess the success of the knockout.The genome of BJ3-2 was used as a template for PCR amplification of the HLarm and HRarm of sdaAA and katX, respectively.The cm was amplified by PCR using pHT01cas9-p43 as the template.The fragment and pUC18 were digested by the corresponding restriction endonuclease EcoR I, Kpn I, BamH I, Sal I, and Pst I and then linked with T4 DNA ligase.The connected vector was then transformed into E. coli DH5α-competent cells.Finally, knockout vectors pUC18-(HLarm-cm-HRarm)/sdaAA and pUC18-(HLarm-cm-HRarm)/katX were constructed.Then, vectors were transformed into BJ3-2 and screened with chloramphenicol; the method used for the transformation was chemical transformation [19,20].PCR validation of positively transformed colonies was performed using DEDP (Figures S1 and S2).Sequencing was then performed to validate the knockout strains (Figures S4 and S5).BJ3-2∆sdaAA and BJ3-2∆katX were obtained.In addition, the genome of BJ3-2 was used as a template to amplify the sdaAA, and the pMarA was used as a template for PCR amplification of the erm.
The aroma profile evaluation of FSs was assessed by quantitative descriptive analysis (QDA), and the sensory evaluation team consisted of 10 panelists from the School of Life Sciences, Guizhou University (4 males and 6 females, 20-27 years old).Each team member underwent a one-week sensory training course on fermented soybeans.Firstly, a descriptive test was conducted to elicit the sensory descriptor.Each FSs (5 g each) was put into 50 mL Teflon vessels, and each member smelled them to freely record the sensory descriptors.The panel then discussed and reached a unanimous agreement based on the frequency and intensity of the sensory attributes of soybean flavor.After that, five sensory attributes were selected to describe the soybean flavor, including "sweaty", "smoky", "beany", "roasted", and "sweet".Finally, for FS samples coded with 3-digit labels, the intensity of the aroma contribution was evaluated using a 6-point scale (0 = no smell, 1 = very weak, 2 = weak, 3 = moderate, 4 = strong, and 5 = very strong with significant irritation) [21], and each sample was analyzed three times.The flavor radar chart was plotted based on the obtained five-dimensional aroma wheel.

Extraction of Volatile Compounds in FS Samples by Using HS-SPME
HS-SPME extraction of FSs was performed according to Peng et al. [22] by putting 2.5 mL of distilled water, 0.625 g of NaCl, and 1.5 g of the fermented soybean sample into a 15 mL headspace vial.Samples were then pre-equilibrated in a 50 • C water bath (Dk-s24, Shanghai, China) for 20 min, after which the aged SPME (50/30 µm DVB/CAR/PDMS (2 cm); Bellefonte, PA, USA) was placed 1 cm above the liquid sample and adsorbed for 40 min.Next, SPME fiber was immediately inserted into the GC injection port (250 • C) of the GC-system and desorbed for 5 min.Each sample was analyzed three times.

Identification of Aroma-Active Compounds in FS Samples by GC-O-MS
GC-O-MS analysis was performed using an Agilent 8890B GC (Agilent Technologies, Santa Clara, CA, USA) equipped with a polar capillary column DB-WAX (60 m × 0.25 mm × 0.25 µm film thickness, Agilent Technologies) and coupled to an olfactory detector (ODP4; Gerstel, Mülheim an der Ruhr, Germany) and an Agilent 5977B MS (Agilent Technologies, Santa Clara, CA, USA).The mass spectrometer and the sniffer detection port received the volatile compounds separated by the capillary column in a 1:1 ratio.Humidified air was pumped into the sniffer arm at 40 mL/min.The following temperature program was applied: the column temperature was initially set at 40 • C and held for 3 min, ramped up to 130 • C at a rate of 3 • C/min and held for 5 min, and finally ramped up to 200 • C at a rate of 5 • C/min and held for 5 min.The GC inlet temperature was 250 • C, and high-purity helium was used as the carrier gas (He > 99.999%) and transported through the capillary column at a constant flow rate of 1.0 L/min.GC-O: The aroma intensity (AI) was evaluated using the approach of Peng et al. with minor adjustments (from a 5-point scale to a 6-point scale) [23].Four panelists (2 males and 2 females, an average of 25 years old) were asked to record the AI values, as well as odor characteristics and retention time.A 6-point scale ranging from 0 (none) to 5 (extreme) was used to evaluate the aroma intensity.A higher number means that the compound has a stronger odor.Analyses was repeated three times by each panelist.Finally, the AI values were the average of four panelists.
For qualitative and quantitative analysis, volatile compounds were identified using mass spectra from the database (NIST 20.0) with reference to retention index (RI) and odor characteristics from the literature.Some important compounds were further identified using standard compounds.Retention indices were calculated for individual compounds based on n-alkane standards from C7 to C30.
For quantification, 8 µL of 2-octanol (0.18 mg/mL) was added to fermented soybeans as an internal standard, and all compounds were analyzed semi-quantitatively, and the relative concentration of each compound was calculated following the method of Zhao et al. [24].The concentration of the volatile compounds is calculated from the peak areas of the volatiles and the concentration of the internal standard.

Calculation of Relative Odor Activity Values (ROAVs)
ROAVs were calculated based on the ratio of the relative content of volatile compounds to their absolute thresholds in water using the following formulae [25]: where Cx is the average concentration of the compounds and OTx is the odor threshold of the compound.The thresholds used in this study were those in water obtained from published literature [26,27].

Aroma Addition Experiment
Aroma addition experiments were conducted to compare the aroma profile changes after the key aroma compounds screened from GC-O and ROAV were added to 1.5 g of FSs.Two FSs with BJ3-2 and BJ3-2∆sdaAA∆katX, which have greater differences in soybean flavor, were selected.The differences in concentrations of key aroma compounds in the two FSs were calculated from the semi-quantitative results (Table S5).All chemical solutions were prepared in water or propylene glycol solution and diluted to the corresponding concentrations [28] (Table S5).Then, 1.5 g of each of the two FSs matrices were placed into 15 mL Teflon vessels; 11 compounds were added to FSs with BJ3-2, and 2 compounds were added to FSs with BJ3-2∆sdaAA∆katX to achieve the same concentration in both FSs (Table S5).The two FSs were then equilibrated at 25 • C for 30 min.The aroma profiles of the two FSs were then evaluated by the ten panelists using QDA, and the flavor radar chart was drawn based on average.

Statistical Analysis
At least three replications of each experiment and sample analysis were conducted, and the results were reported as mean ± standard deviation.Data were statistically analyzed by SPSS software (version 26.0, Inc., Chicago, IL, USA) using a t-test to ascertain a significant difference between BJ3-2, BJ3-2∆sdaAA, BJ3-2∆katX, and BJ3-2∆sdaAA∆katX at p ≤ 0.05.Origin 2022 software (version 9.8.0.200) was used for plotting.The correlation between sensory attributes and key aroma compounds was analyzed by PLSR using the Unscrambler 9.7 version software (CAMO ASA, Oslo, Norway).TBtools v2.042 was used to make the heat map.

Screening of Soybean Flavor Candidate Genes of BJ3-2
To screen genes that may significantly impact the formation of soybean flavor, we conducted further analysis and screening of the prior multi-omics data from BJ3-2 incubated at different temperatures, as described in the literature [29].Briefly, differential genes, proteins, and metabolites with |Log 2 FC (FOLDCHANGE)| ≥ 2, |Log 2 FC| ≥ 1.5, and |Log 2 FC| ≥ 1 were screened according to p-value < 0.05 and from the transcriptome, proteome, and metabolome for KEGG pathway mapping.Co-analyzed genes and metabolites are mainly involved in the glycine, serine, and threonine metabolism; ABC transporter, arginine and proline metabolism; cysteine and methionine metabolism; and tryptophan metabolism metabolic pathways.Meanwhile, differential genes and proteins with |Log 2 FC| ≥ 2.5 were screened from the joint analysis results, and a total of sdaAA and katX were studied as candidate genes.In the transcriptome data, sdaAA was downregulated by 2.50~fold (log 2 FC) after incubating at 53 • C, and katX was down-regulated by 4.40~fold (log 2 FC) after incubating at 53 • C (Table S2).RT-qPCR was further performed to confirm the gene expression of BJ3-2 (Table S3).As shown in Figure 1, the results show that the expression of sdaAA and katX was low at 53 • C and high at 37 • C, indicating the transcriptome data and the RT-qPCR results were consistent.
arginine and proline metabolism; cysteine and methionine metabolism; and tryptop metabolism metabolic pathways.Meanwhile, differential genes and proteins w |Log2FC| ≥ 2.5 were screened from the joint analysis results, and a total of sdaAA and k were studied as candidate genes.In the transcriptome data, sdaAA was down-regula by 2.50~fold (log2FC) after incubating at 53 °C, and katX was down-regulated by 4.40~ (log2FC) after incubating at 53 °C (Table S2).RT-qPCR was further performed to con the gene expression of BJ3-2 (Table S3).As shown in Figure 1, the results show that expression of sdaAA and katX was low at 53 °C and high at 37 °C, indicating the transc tome data and the RT-qPCR results were consistent. .

Knockout of sdaAA and katX Had No Significant Effect on the Morphology and Growth BJ3-2
In this study, we constructed two single-knockout strains and one double-knock strain by homologous recombination.After obtaining the knockout strain, it was de mined whether the knockout affected the colonies phenotypic characteristics and gro ability.The morphology of single colonies of BJ3-2 on solid plates was first identified.onies appeared white, dry, irregular in morphology, and internally wrinkled.The m phology of BJ3-2ΔsdaAA, BJ3-2ΔkatX, and BJ3-2ΔsdaAAΔkatX colonies was comparab that of BJ3-2 (Figure 2A,C,E,G).
Furthermore, the four strains (BJ3-2, BJ3-2ΔsdaAA, BJ3-2ΔkatX, and BJ3-2Δsd AΔkatX) were stained with a blue-purple tint and a short stick form by Gram staining shown in Figure 2B,D,F,H.The growth curves were then used to determine the gro ability of the strains.The results are shown in Figure 2I.The OD values of BJ3-2Δsda BJ3-2ΔkatX, and BJ3-2ΔsdaAAΔkatX were similar to those of BJ3-2 at different time po These findings showed that the morphology and growth of BJ3-2 were not significa affected by the knockout of sdaAA and katX.

Knockout of sdaAA and katX Had No Significant Effect on the Morphology and Growth of BJ3-2
In this study, we constructed two single-knockout strains and one double-knockout strain by homologous recombination.After obtaining the knockout strain, it was determined whether the knockout affected the colonies' phenotypic characteristics and growth ability.The morphology of single colonies of BJ3-2 on solid plates was first identified.Colonies appeared white, dry, irregular in morphology, and internally wrinkled.The morphology of BJ3-2∆sdaAA, BJ3-2∆katX, and BJ3-2∆sdaAA∆katX colonies was comparable to that of BJ3-2 (Figure 2A,C,E,G).
Furthermore, the four strains (BJ3-2, BJ3-2∆sdaAA, BJ3-2∆katX, and BJ3-2∆sdaAA∆katX) were stained with a blue-purple tint and a short stick form by Gram staining, as shown in Figure 2B,D,F,H.The growth curves were then used to determine the growth ability of the strains.The results are shown in Figure 2I.The OD values of BJ3-2∆sdaAA, BJ3-2∆katX, and BJ3-2∆sdaAA∆katX were similar to those of BJ3-2 at different time points.These findings showed that the morphology and growth of BJ3-2 were not significantly affected by the knockout of sdaAA and katX.

Sensory Analysis of Soybean Flavor in FSs with BJ3-2, BJ3-2ΔsdaAA, BJ3-2ΔkatX, and BJ3-2ΔsdaAAΔkatX at 37 °C
Quantitative descriptive analyses (QDAs) were applied to assess the differences in aroma profiles of four FSs after fermentation of soybeans with BJ3-2, BJ3-2ΔsdaAA, BJ3-2ΔkatX, and BJ3-2ΔsdaAAΔkatX at 37 °C (Figure 3A).As shown in Figure 3B, the aroma profiles of soybean flavor from FSs were described by five sensory attributes: "sweaty", "smoky", "beany", "roasted", and "sweet".BJ3-2 was more frequently associated with the "sweaty (4.5 ± 0.45)" attribute, whereas the other attributes were at low levels, probably because acids obscure other attributes.BJ3-2ΔsdaAA and BJ3-2ΔkatX had diminished "sweaty" attributes and enhanced "smoky", "beany", and "sweet" attributes compared to BJ3-2.BJ3-2ΔsdaAAΔkatX presented the lowest intensity of the "sweaty (2.3 ± 0.32)" attribute and the highest intensity of the "beany (2.6 ± 0.12)" and "roasted (3.1 ± 0.33)" attributes, with higher differences from BJ3-2.In addition, the "sweet" attribute had the lowest intensity of soybean flavor in FSs, and knockout of sdaAA and katX resulted in an increase in the intensity of the "sweet" attribute.The QDA results suggest that knockout of sdaAA and katX has a substantial effect on the flavor formation of soybean flavor and also validate the analysis of the mechanism of soybean flavor formation using multi-omics data.Quantitative descriptive analyses (QDAs) were applied to assess the differences in aroma profiles of four FSs after fermentation of soybeans with BJ3-2, BJ3-2∆sdaAA, BJ3-2∆katX, and BJ3-2∆sdaAA∆katX at 37 • C (Figure 3A).As shown in Figure 3B, the aroma profiles of soybean flavor from FSs were described by five sensory attributes: "sweaty", "smoky", "beany", "roasted", and "sweet".BJ3-2 was more frequently associated with the "sweaty (4.5 ± 0.45)" attribute, whereas the other attributes were at low levels, probably because acids obscure other attributes.BJ3-2∆sdaAA and BJ3-2∆katX had diminished "sweaty" attributes and enhanced "smoky", "beany", and "sweet" attributes compared to BJ3-2.BJ3-2∆sdaAA∆katX presented the lowest intensity of the "sweaty (2.3 ± 0.32)" attribute and the highest intensity of the "beany (2.6 ± 0.12)" and "roasted (3.1 ± 0.33)" attributes, with higher differences from BJ3-2.In addition, the "sweet" attribute had the lowest intensity of soybean flavor in FSs, and knockout of sdaAA and katX resulted in an increase in the intensity of the "sweet" attribute.The QDA results suggest that knockout of sdaAA and katX has a substantial effect on the flavor formation of soybean flavor and also validate the analysis of the mechanism of soybean flavor formation using multi-omics data.
The aroma compounds 2,5-dimethylpyrazine, 2,3,5-trimethylpyrazine, and 3-ethyl-2,5-dimethylpyrazine have previously been reported in natto [12].These pyrazines typically provide a "roasted" or "nutty" attribute.In this study, these three pyrazines showed different AI values among the four FSs.Guaiacol also showed a relatively higher AI value in four FSs, and it has a smoky flavor, which should be related to the "smoky" attributes of soybean flavor.1-Octen-3-ol imparts a mushroom flavor, which is here related to the "beany" attribute.In addition, 2,4,5-trimethyloxazole contributed to the "roasted" attribute and 2-pentylfuran to the "beany" attribute.
The aroma compounds 2,5-dimethylpyrazine, 2,3,5-trimethylpyrazine, and 3-ethyl-2,5dimethylpyrazine have previously been reported in natto [12].These pyrazines typically provide a "roasted" or "nutty" attribute.In this study, these three pyrazines showed different AI values among the four FSs.Guaiacol also showed a relatively higher AI value in four FSs, and it has a smoky flavor, which should be related to the "smoky" attributes of soybean flavor.1-Octen-3-ol imparts a mushroom flavor, which is here related to the "beany" attribute.In addition, 2,4,5-trimethyloxazole contributed to the "roasted" attribute and 2-pentylfuran to the "beany" attribute.Several other compounds, including 2-ethylfuran, 2,3-butanedione, methyl 2-methylbutyrate, 2-heptanone, and 3-octanone, were associated with the "sweet" attribute, as these compounds were perceived as having a pleasant flavor, neutralizing the overall flavor of soybean.Concentration differences mainly caused differences in AI values in fermented soybean samples from different strains.Therefore, these compounds need to be quantified and relative odor activity values calculated.
in FSs.The 13 key aroma compounds are listed in Figure 4.The chemical structures of these compounds were further identified by comparing the retention index (RI) and aroma profile with those of the standard compounds and used for further analysis.

Relationship between Key Aroma Compounds and Sensory Attributes of Soybean Flavor in FSs
Partial least squares regression (PLSR) was performed to investigate the relationship between sensory attributes and key aroma compounds in the four FS samples.The X-matrix was set as the thirteen key aroma compounds, and the Y-matrix was set as five sensory attributes.The PLSR consisted of two significant PCs explaining 75% of the X-variance and 92% of the Y-variance (Figure 5), while the inner ellipse and outer ellipse explained 50% and 100% of the variance, respectively.All variables are located outside the small ellipse except 2-ethylfuran (1) and 3-octanone (7), which shows that these variables were well explained by the PLSR model.The low values of RMSEC (0.11) and RMSRP (0.19) also indicate that the regression model performs well.As shown in Figure 5, the closer to the variable Y, the higher the correlation between X and Y [33].The compound 2-Methylbutyric acid ( 12) was highly associated with the "sweaty" attribute.The "smoky" attribute was strongly influenced by guaiacol (13).In addition, 2,3,5-trimethylpyrazine ( 9) and 1octen-3-ol (10) showed robust correlation with "roasted" and "beany" attributes, respectively.In addition, the "sweet" attribute was weakly associated with certain aroma compounds.This might be because more compounds are implicated in the formation of this sensory attribute.These results are consistent with the GC-O and ROAV analyses (Table 1).

Relationship between Key Aroma Compounds and Sensory Attributes of Soybean Flavor in FSs
Partial least squares regression (PLSR) was performed to investigate the relationship between sensory attributes and key aroma compounds in the four FS samples.The X-matrix was set as the thirteen key aroma compounds, and the Y-matrix was set as five sensory attributes.The PLSR consisted of two significant PCs explaining 75% of the X-variance and 92% of the Y-variance (Figure 5), while the inner ellipse and outer ellipse explained 50% and 100% of the variance, respectively.All variables are located outside the small ellipse except 2-ethylfuran (1) and 3-octanone (7), which shows that these variables were well explained by the PLSR model.The low values of RMSEC (0.11) and RMSRP (0.19) also indicate that the regression model performs well.As shown in Figure 5, the closer to the variable Y, the higher the correlation between X and Y [33].The compound 2-Methylbutyric acid (12) was highly associated with the "sweaty" attribute.The "smoky" attribute was strongly influenced by guaiacol (13).In addition, 2,3,5-trimethylpyrazine ( 9) and 1-octen-3-ol (10) showed robust correlation with "roasted" and "beany" attributes, respectively.In addition, the "sweet" attribute was weakly associated with certain aroma compounds.This might be because more compounds are implicated in the formation of this sensory attribute.These results are consistent with the GC-O and ROAV analyses (Table 1).
As shown in Figure 6, in comparison to FSs with BJ3-2, the content of 2-methylbutyric acid in FSs with BJ3-2ΔkatX, BJ3-2ΔsdaAA, and BJ3-2ΔsdaAAΔkatX was significantly lower, with BJ3-2ΔsdaAAΔkatX having the lowest content.The lowest content of 2-methylbutyric acid in BJ3-2ΔsdaAAΔkatX also explains well why it possesses the lowest intensity of "sweaty" attributes (Figure 3B).This compound may have a positive correlation with the intensity of "sweaty" attributes.2,3,5-Trimethylpyrazine was decreased in BJ3-2Δs-daAA and significantly increased in both BJ3-2ΔkatX and BJ3-2ΔsdaAAΔkatX.Meanwhile, in contrast to BJ3-2, the "roasted" attribute of BJ3-2sdaAA was weakened, whereas the "roasted" attribute of BJ3-2ΔkatX and BJ3-2ΔsdaAAΔkatX was enhanced (Figure 3B), indicating that 2,3,5-trimethylpyrazine was associated with the "roasted" attribute.Guaiacol and 1-octen-3-ol were increased in BJ3-2ΔkatX, BJ3-2ΔsdaAA, and BJ3-2ΔsdaAAΔkatX, which is consistent with the QDA results (Figure 3B) of enhanced "smoky" and "beany" attributes of soybean flavor in FSs.Several other compounds also showed some variations in the content of the four FSs.In summary, knockout of sdaAA and katX resulted in a significant increase or decrease in the content of key aroma compounds, which mainly resulted in a significant decrease in 2-methylbutyric acid and a significant increase in 2,5dimethylpyrazine, 2-ethylfuran, 2-pentylfuran, 2-heptanone, 1-octen-3-ol, and guaiacol, thus affecting the soybean flavor.
Additionally, KEGG pathway analysis shows that sdaAA (which encodes serine dehydratase) is involved in metabolic pathways, mainly including glycine, serine, and threonine metabolism and cysteine and methionine metabolism.katX (encodes catalase) is involved in metabolic pathways, mainly tryptophan metabolism, glyoxylate metabolism, and dicarboxylate metabolism.This implies that the knockout of sdaAA and katX may impact the production of aroma compounds by affecting these metabolic pathways, leading to altered sensory attributes in FS with its mutants.However, the roles and specific regulatory mechanisms of sdaAA and katX in the production of soybean flavor still require further investigation.Based on semi-quantitative results, significant difference analyses (p ≤ 0.05) were performed on the content of 13 key aroma compounds between the four FSs with BJ3-2, BJ3-2∆sdaAA, BJ3-2∆katX, and BJ3-2∆sdaAA∆katX.
Additionally, KEGG pathway analysis shows that sdaAA (which encodes serine dehydratase) is involved in metabolic pathways, mainly including glycine, serine, and threonine metabolism and cysteine and methionine metabolism.katX (encodes catalase) is involved in metabolic pathways, mainly tryptophan metabolism, glyoxylate metabolism, and dicarboxylate metabolism.This implies that the knockout of sdaAA and katX may impact the production of aroma compounds by affecting these metabolic pathways, leading to altered sensory attributes in FS with its mutants.However, the roles and specific regulatory mechanisms of sdaAA and katX in the production of soybean flavor still require further investigation.

Aroma Addition Experiment
Based on our analyses and further validation of key aroma compounds, BJ3-2 and BJ3-2ΔsdaAAΔkatX, which have greater differences in aroma profiles, were selected for an aroma addition experiment.Aroma addition experimental models were prepared by adjusting the contents of key aroma compounds in the two samples to be consistent based on the difference in concentrations of the 13 key aroma compounds in the two FSs.After

Aroma Addition Experiment
Based on our analyses and further validation of key aroma compounds, BJ3-2 and BJ3-2∆sdaAA∆katX, which have greater differences in aroma profiles, were selected for an aroma addition experiment.Aroma addition experimental models were prepared by adjusting the contents of key aroma compounds in the two samples to be consistent based on the difference in concentrations of the 13 key aroma compounds in the two FSs.After that, it was determined whether there was a change in the aroma profiles of soybean flavor in the two FSs or whether the aroma profiles converged to verify the effect of these key aroma compounds.As shown in Figure 7, QDA was applied to access the aroma profiles of the two FSs after adding the aroma compounds.The aroma profiles of both samples changed after the addition of the differential aroma compounds, and both samples converged in "sweaty", "smoky", "roasted", and "beany" attributes but still differed in "sweet".This may be due to the fact that the "sweet" attribute of BJ3-2 is masked by the strong intensity of the "sweaty" attribute.The result of the aroma addition experiment suggests that these compounds are closely related to the aroma profiles of soybean flavor and also validates these substances as key aroma compounds of soybean flavor in FS.
Fermentation 2024, 10, x FOR PEER REVIEW 13 of 16 that, it was determined whether there was a change in the aroma profiles of soybean flavor in the two FSs or whether the aroma profiles converged to verify the effect of these key aroma compounds.As shown in Figure 7, QDA was applied to access the aroma profiles of the two FSs after adding the aroma compounds.The aroma profiles of both samples changed after the addition of the differential aroma compounds, and both samples converged in "sweaty", "smoky", "roasted", and "beany" attributes but still differed in "sweet".This may be due to the fact that the "sweet" attribute of BJ3-2 is masked by the strong intensity of the "sweaty" attribute.The result of the aroma addition experiment suggests that these compounds are closely related to the aroma profiles of soybean flavor and also validates these substances as key aroma compounds of soybean flavor in FS.

Source of Key Aroma Compounds of Soybean Flavor in Fermented Soybeans
The results of the key aroma compounds of soybean flavor in FS showed that 2methylbutyric acid produced the main compounds contributing to "sweaty" attributes.Moreover, 2-methylbutyric acid was identified as a short-chain fatty acid.The source of 2methylbutyric acid may be due to the ability of B. subtilis to metabolically convert isoleucine, leucine, and valine to 2-methylbutyric acid during fermentation [34].Meanwhile, 2methylbutyric acid can be converted to methyl 2-methylbutyrate [35].
Guaiacol is also found to be an important aroma compound and can be produced by the degradation of vanillin by the enzyme vanillin dehydrogenase in B. subtilis [38].Moreover, 2,3-butanedione usually has a pleasant, creamy, or buttery taste and can be metabolized by B. subtilis to produce α-acetolactate (which can be converted to acetoin) and 2,3butanedione (which can also be converted to each other under the catalysis of acetoin/butanediol dehydrogenase) [39].Cysteine/alanine has been reported to react with 2,3-butanedione to form 2,4,5-trimethyloxazole [40], and 2-heptanone can also be produced by the Maillard reaction or biosynthesized from acid [41].The metabolic pathways of 3-

Source of Key Aroma Compounds of Soybean Flavor in Fermented Soybeans
The results of the key aroma compounds of soybean flavor in FS showed that 2methylbutyric acid produced the main compounds contributing to "sweaty" attributes.Moreover, 2-methylbutyric acid was identified as a short-chain fatty acid.The source of 2-methylbutyric acid may be due to the ability of B. subtilis to metabolically convert isoleucine, leucine, and valine to 2-methylbutyric acid during fermentation [34].Meanwhile, 2-methylbutyric acid can be converted to methyl 2-methylbutyrate [35].
Guaiacol is also found to be an important aroma compound and can be produced by the degradation of vanillin by the enzyme vanillin dehydrogenase in B. subtilis [38].Moreover, 2,3-butanedione usually has a pleasant, creamy, or buttery taste and can be metabolized by B. subtilis to produce α-acetolactate (which can be converted to acetoin) and 2,3-butanedione (which can also be converted to each other under the catalysis of acetoin/butanediol dehydrogenase) [39].Cysteine/alanine has been reported to react with 2,3-butanedione to form 2,4,5-trimethyloxazole [40], and 2-heptanone can also be produced by the Maillard reaction or biosynthesized from acid [41].The metabolic pathways of 3-octanone and 1-octen-3-ol in fungi are well defined, whereas the biosynthetic pathway in B. subtilis has not been elucidated.The fungal pathway is mainly based on the oxidative degradation of unsaturated fatty acids [42,43].

Conclusions
In this study, the aroma profiles and key aroma compounds of soybean flavor in FS with B. subtilis BJ3-2 and its knockout mutants at 37 • C were analyzed using QDA, GC-O-MS, ROAV calculation, and aroma addition experiments.QDA results revealed that the overall aroma of soybean flavor in FSs consists of "sweaty", "smoky", "roasted", "beany", and "sweet" attributes.In addition, the knockout of soybean flavor candidate genes (sdaAA and katX) did not significantly affect the morphology and growth of BJ3-2 but had a greater effect on the aroma profiles of soybean flavor in FSs with BJ3-2.A total of 20 aroma-active compounds were detected by GC-O-MS from four FSs with BJ3-2, BJ3-2∆sdaAA, BJ3-2∆katX, and BJ3-2∆sdaAA∆katX, and 13 key aroma compounds were identified as key aroma compounds according to their high ROAVs.Furthermore, PLSR revealed a good correlation between the key aroma compounds and sensory attributes.Knockout of sdaAA and katX resulted in a significant increase or decrease in the content of key aroma compounds.Finally, combined with the aroma addition experiments, 2methylbutyric acid, 2,3,5-trimethylpyrazine, guaiacol, and 1-octen-3-ol were identified as the aroma compounds responsible for the main aroma profiles of soybean flavor.The findings of this study provide new insights into the analysis of the soybean flavor of fermented soybeans produced by B. subtilis BJ3-2, which is important for quality control and improvement of the aroma of fermented soybeans.

Institutional Review Board Statement:
As the experimental samples used in the study were consumed in daily life, the study was exempted from ethical review and approval.
Informed Consent Statement: Informed consent was obtained from all panelists involved in the study.Participants agreed to the use of their data/answers.The study strictly adhered to a protocol designed to protect the rights and privacy of all participants.
C enhance soybean flavor compared to those at 45 • C and 53 • C, while fermented soybeans at 45 • C and 53 • C enhance soy sauce-like flavor compared to those at 37 • C [5].Although traditional soybean flavor products such as Douchi and Fermentation 2024, 10, 409 2 of 16

Figure 4 .
Figure 4. Information on thirteen key aroma compounds from GC-O (AI) combined with ROAV. a Thirteen key aroma compounds from GC-O combined with ROAV; b aroma characteristics of key aroma compounds; c odor thresholds of compounds in water reported in the literature.

Figure 4 .
Figure 4. Information on thirteen key aroma compounds from GC-O (AI) combined with ROAV. a Thirteen key aroma compounds from GC-O combined with ROAV; b aroma characteristics of key aroma compounds; c odor thresholds of compounds in water reported in the literature.

Figure 6 .
Figure 6.Changes in content of key aroma compounds among the four fermented soybean samples with BJ3-2, BJ3-2ΔsdaAA, BJ3-2ΔkatX, and BJ3-2ΔsdaAAΔkatX.Significant differences between groups are indicated by letters a, b, c.
Note: Orange color in the table indicates compounds that are common to both GC-O and ROAV > 1.Data are presented as mean ± SD.