Relationship between the Dynamics of Flavor Compounds and Microbial Succession in the Natural Fermentation of Zhalajiao, a Popular Traditional Chinese Fermented Chili Paste

Zhalajiao, a traditional Chinese fermented food, is popular due to its unique flavor. Traditional Zhalajiao fermentation is closely related to flavor compounds production. However, the mechanisms underlying the formation of these crucial flavor components in Zhalajiao remain unclear. Here, we explored the dynamic changes in physical and chemical properties, microbial diversity, and flavor components of Zhalajiao at various fermentation times. In total, 6 organic acids, 17 amino acids, and 21 key volatile compounds were determined as flavor components. In Zhalajiao, Lactobacillus and Cyanobacterium were the main bacteria that were involved in the formation of crucial flavor compounds. Candida showed a significant correlation with 14 key flavor compounds during fermentation (p < 0.05) and was the main fungal genus associated with flavor formation in Zhalajiao. This research offers a theoretical foundation for the flavor regulation and quality assurance of Zhalajiao.


Introduction
Zhalajiao is a traditional Chinese fermented food with distinctive characteristics that has become increasingly popular among people in the central and south-western regions of China due to its vibrant color, appetizing sourness, and mellow aftertaste [1].Traditionally, it is produced through spontaneous fermentation by using cereals (mainly corn or rice) and fresh chilies (Xiaomila, Erjingtiao, and Niujiaojiao) as raw materials under the action of naturally occurring microorganisms [2].The entire fermentation cycle lasts 30-90 days at room temperature.With its characteristic flavors of sourness, aroma, spiciness, and freshness, Zhalajiao can be served as a standalone dish after being fried or used as a spice to improve the taste of other dishes [3].During fermentation, carbohydrates, lipids, proteins, and other substances in Zhalajiao are decomposed, resulting in the formation of various flavor components, such as organic acids, free amino acids (FAAs), and volatile compounds (esters, acids, terpenes), which have been shown to significantly contribute to the flavor characteristics of fermented Zhalajiao [4].Additionally, microorganisms are crucial to the creation of flavor components during fermentation.The diversity of microbial composition places fermented foods in a complex micro-ecological environment, leading to the diversity of the flavors of these products.Given the relatively uncontrolled fermentation of Zhalajiao, the presence of various microorganisms in the air, water, and fermentation tanks and the use of various raw materials result in the formation of a microbial community with a complex structure in the final product.In traditional fermented foods, this complex natural microbial community is related to physicochemical properties and affects flavor characteristics [5].
Some research has been performed on microbial diversity in Zhalajiao products.Guo et al. [3] used high-throughput sequencing and culture-dependent methods to collect and study corn-based Zhalajiao samples from two regions.They found remarkable discrepancies in the bacterial community composition of corn-based Zhalajiao samples from two regions.Companilactobacillus, Lactiplantibacillus, and Corynebacterium were the common dominant genera in the two samples.Cai et al. [1,2,6] studied the bacterial community and flavor characteristics of rice-based Zhalajiao produced in different regions from differing raw materials.Rice-based Zhalajiao made in different regions from varying ingredients had significantly distinct microbial communities and a sensory quality, with Lactobacillus being the predominant bacterial genus.However, the above research merely discussed the bacterial community at the culmination of rice-based Zhalajiao fermentation.In a previous research paper, corn-based Zhalajiao with a fermentation time of 0-60 days was collected to study the dynamic changes in microflora and flavor components during its fermentation [4].Lactobacillus and Kazachistania were considered to be the most preponderant bacterial and fungal genera, respectively.Furthermore, Lactobacillus, Acetobacter, Kazachistania, Wickerhamyces, Kluyveromyces, and Cyberlindnera were the major microbes participating in the generate of flavor components.Nevertheless, to our knowledge, no research on the composition and succession of bacteria and fungi in the fermentation of rice-based Zhalajiao is available, and the relationship among the major flavor compounds, physicochemical properties, and associated microorganisms of this product is unclear.Therefore, studying the dynamic succession of physicochemical properties, flavor compounds, and microflora of rice-based Zhalajiao during fermentation and the relationship among these factors is crucial for the production and quality control of rice-based Zhalajiao.
For this study, the changes in physicochemical indices, organic acids, FAAs, and volatile components during Zhalajiao fermentation were measured.High-throughput sequencing was applied to analyze the combination of and change in microbial communities during Zhalajiao fermentation.In addition, the relationship between flavor components and microflora was explored.The aim of this research is to lay the foundation for the identification of the core functional microbial and improvement in the quality of traditional Zhalajiao.

Zhalajiao Preparation and Collection
Fresh red chilies were cleaned and cut into pieces and then mixed well with an equal weight of rice flour and 6% (w/w) salt.After being crushed and mixed, the raw materials were placed in a clean fermentation tank.The mouth of the tank was sealed, and the tank was placed upside down in a water basin.The chilies were naturally fermented for 90 days in Tong zi, Guizhou Province, China.Eight samples were collected on days 0, 7, 15, 22, 30, 45, 60, and 90 of fermentation and denoted as D0, D7, D15, D22, D30, D45, D60, and D90, respectively.The samples were kept at 4 • C for the measurement of physicochemical indices, FAAs, organic acids, and volatile compounds.For DNA extraction and Illumina MiSeq sequencing, the samples (20 g) were promptly kept at 80 • C after being collected at different fermentation cycles.

Physicochemical Properties
The moisture content in Zhalajiao was measured through direct drying [7].The pH was determined with a pH meter (PHS-3C, Shanghai Instrument and Electronics Science Instrument Co., Ltd., Shanghai, China) in reference to Xu et al. [8].Total acid was determined in accordance with [9].Reducing sugars was measured by 3,5-dinitrosalicylic acid colorimetry [10].Total soluble sugars were measured using anthrone colorimetry in accordance with [11].

Determination of Organic Acids and FAAs
Organic acids were examined using high-performance liquid chromatography.The chromatographic conditions were as described by Xiao et al. [12] with slight modifications.Separations were performed on Agilent AAAC18 (4.6 mm × 150 mm i.d., 5 µm).The column temperature was set at 28 • C. The mobile phase comprised 0.01 mol/L diammonium hydrogen phosphate (pH 2.5) and methanol (98:2, v:v), and the flow rate was 0.6 mL/min.The injection volume was 10 µL, and the detection wavelength was 220 nm.
The FAAs in Zhalajiao were detected with an automatic amino acid analyzer in accordance with a previous report with appropriate modifications [13].The sample of Zhalajiao was settled with 2.0 mL of 10% sulfosalicylic acid for 2 h at 4 • C and then centrifuged at 8000 r/min for 1 min.Then, 1.0 mL of the supernatant fluid was filtrated using a 0.22 µm microporous membrane and determined on the L-8900 automatic amino acid analyzer (Hitachi, Tokyo, Japan).Various FAA standards were used for quantification.

Determination of Volatiles
Volatile compounds were determined in reference to the method of Chen et al. [14] with slight modifications.In total, 5 g of samples was placed in a 20 mL headspace injection vial containing 10 µL of 2-octanol standard with a mass concentration of 500 mg/L as an internal standard.Then, the vial was equilibrated for 30 min at 50 • C in a thermostat water bath.Next, an HS-SPME fiber (100 µm PDMS) was inserted into the vial to absorb the volatile compounds for 40 min.Finally, the head of the fiber was desorbed for 5 min at 250 • C in a gas chromatographic mass spectrometer (GCMS-GP2010, Shimadzu, Kyoto, Japan).The chromatographic column used was a DB-5MS capillary column with a size of 30 m × 0.25 mm × 0.25 µm, and the column carrier gas was helium at a flow rate of 1.0 mL/min.The temperature program was 40 The mass spectrum was handled in the electron ionization mode at 70 eV with an ion-source temperature of 250 • C and a scan range from 40 m/z to 400 m/z.Volatile compounds were qualitatively analyzed using a combination of spectral library searches such as NIST05 and NIST08 (compounds with matches > 80) and retention indices.The retention index (RI) of each volatile component was calculated from the retention time of C7~C40 n-alkanes [13].The quantification of each volatile component in the samples was performed using the internal standard method.The calculation formula was as follows: where Ci is the concentration of the identified compound (µg/kg), ρ is the mass concentration of the internal standard (µg/µL); V is the volume of the internal standard (µL); Ai is the peak area of each volatile component; A is the peak area of the internal standard substance; and m is the mass of the sample (kg) [8].

Calculation of Odor Activity Value
Odor activity value (OAV) was applied to characterize the contribution of various aroma compounds to the main aroma components in Zhalajiao.It was determined by dividing the concentration of each compound in water by its detection threshold [4].OAV > 1 indicated that the compound is conductive to the aroma of the sample and is a key volatile component.A high OAV is indicative of the high individual contribution of a compound [13].
2.6.Microbiological Analysis 2.6.1.DNA Extraction and PCR Amplification Total DNA was extracted from each Zhalajiao sample by using an E.Z.N.A Soil DNA Kit (Omega Bio-Tek, Norcross, GA, USA) in accordance with the specification.The concentration and pureness of extracted DNA were tested using Nanodrop 2000 (Thermo Fisher Scientific, Waltham, MA, USA) and 1% agarose gels.In reference to a previous report with appropriate modifications [8], the V3-V4 variable regions of the bacterial 16S rDNA genes were amplified using primers 338F and 806R following the parameters: pre-denaturation at 95 • C for 3 min; 27 cycles at 95 • C for 30 s, 55 • C for 30 s, and 72 • C for 45 s; and final extension at 72 • C for 10 min.Each amplification was supplemented with ddH2O to a volume of 20 µL and included 4 µL of 5× FastPfu Buffer, 2 µL of 2.5 mM dNTPs, 0.8 µL of primers (5 µM), 4 µL of FastPlu Polymerase, and 10 ng of template DNA.The fungal ITS1-ITS2 gene regions were amplified with ITS1F and ITS2R primers in accordance with Chen et al. [14].The PCR was implemented as follows: pre-denaturation at 95 • C for 3 min; 35 cycles at 95 • C for 30 s, 55 • C for 30 s, and 72 • C for 45 s; and final extension at 72 • C for 10 min.Each amplification was supplemented with ddH2O to a volume of 20 µL and included 2 µL of 10× buffer, 2 µL of 2.5 mM dNTPs, 0.8 µL of primers (5 µM), 0.2 µL of rTaq polymerase, 0.2 µL of BSA, and 10 ng of template DNA.

Sequencing and Data Analysis
The amplified PCR products were examined by electrophoresis with 2% agarose gels and purified using an AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA).After purification and building of the paired-end library, sequencing was performed using an Illumina HiSeq 2500 platform (Majorbio Bio-pharm Technology Co., Ltd., Shanghai, China).The sequence data were analyzed and quality filtered using QIIME (version 1.17).After being merged by FLASH and quality filtered, the sequences were clustered into operational taxonomic units (OTUs) at a 97% similarity threshold by UPARSE, and UCHIME was employed to identify chimeras in accordance with Wang et al. [15].

Statistical Analysis
All the experiments were repeated three times in parallel, and values were expressed as mean ± sd.SPSS statistical software 20.0 (IBM, Chicago, IL, USA) was applied for statistical analysis, and Pearson correlation coefficients between microorganisms and flavor substances were calculated.Significant differences among groups were determined by a one-way analysis of variance with Tukey's honest significant difference post hoc test (p < 0.05).Figures were generated with Origin 2021 8.6 software.The correlation clustering heat map and networks were performed using the Lianchuan Biological Cloud Platform (https://www.omicstudio.cn/toolaccessed on 10 May 2023).

Changes in Physicochemical Properties and Organic Acids during Fermentation
The physicochemical indices and organic acid determination results of all samples (D0, D7, D15, D22, D30, D45, D60, and D90) during Zhalajiao fermentation are shown in Table 1.The results showed that during fermentation, the moisture content of Zhalajiao increased from 45.03% (D0) to 51.31% (D90).The pH and total acid content of the samples also changed significantly (p < 0.05).The raw chili had an original pH of 4.66 ± 0.01, which decreased to an ultimate pH of 4.21 ± 0.03 after 90 days of fermentation.The total acid content showed the lowest value on day 0 (4.32 ± 0.06 g/kg), increased conspicuously to its highest value on day 45 (10.26 ± 0.06 g/kg), then dropped to 9.29 ± 0.03 g/kg on day 90.As critical indicators of fermented vegetables, pH and total acids endow Zhalajiao with a distinctive sour taste that reflects and is influenced by the growth and metabolism of microorganisms during fermentation [16].Furthermore, the growth of harmful microorganisms can be inhibited and shelf life can be extended by using the appropriate pH and acidity levels [17].Total soluble and reducing sugars not only impart sweetness to foods but are also important nutrients for the survival of microorganisms [18].In this study, during the first 15 days of fermentation, a large number of microorganisms utilized sugar fermentation to produce acid, resulting in a sharp decrement in the content of reducing acids from 1.25 ± 0.02 g/kg (D0) to 0.57 ± 0.02 g/kg (D15).However, the reducing sugar content did not change significantly after 45 days of fermentation, likely due to the inhibition of microbial activity with the increase in acidity.Furthermore, during days 0-15 of fermentation, the total soluble sugar content in Zhalajiao did not change significantly and remained at approximately 3.46 g/kg because microorganisms mainly use reducing sugars as carbon and energy sources.In the later stage of fermentation, the content of reducing sugars in Zhalajiao cannot meet the growth and metabolism requirements of microorganisms, leading to the decomposition and utilization of polysaccharides by microorganisms.The total soluble sugar content decreased from 3.54 ± 0.10 g/kg (D22) to 3.20 ± 0.05 g/kg (D45) and remained stable thereafter.Note: Different letters in the same row represent significant differences (p < 0.05).nd: not detected.
Acidity is a very important taste quality of Zhalajiao.Although Zhalajiao fermentation produces a variety of acidic substances, organic acids, such as lactic, acetic, malic, and succinic acids, are the main contributors to acidity.The proportion of organic acids also helps improve the flavor and balance of fermented chili paste [13].The total organic acids in Zhalajiao first increased and then decreased during fermentation and reached the maximum value of 13.92 g/kg on day 15 of fermentation (Table 1).The contents of oxalic and malic acids decreased continuously during the first 22 days of fermentation, increased and stabilized after 30 days of fermentation, and significantly decreased after 60 days of fermentation (p < 0.05).The content of tartaric acid increased continuously during fermentation and remained at approximately 1.1 g/kg after 30 days of fermentation.Succinic acid had an initial content of 0.78 ± 0.041 g/kg and peaked at 3.60 ± 0.064 g/kg on day 15 of fermentation but could not be detected after day 30.The lactic acid content reached a maximum of 3.85 ± 0.073 g/kg on day 15 of fermentation and then declined slightly to 3.21 ± 0.024 g/kg.On day 90 of Zhalajiao fermentation, lactic acid accounted for 33.72% of the total organic acid content, whereas acetic acid accounted for 29.62%.Consistent with our results, the findings of Xu et al. [8] showed that lactic and acetic acids are the most important organic acids in chili paste.Lactic acid is an important metabolite of microbial cells, and lactic acid produced by microbial fermentation is used as a souring agent and preservative in food because of its mellow taste.Two types of lactic acid fermentation exist: homolactic fermentation and heterolactic fermentation.In homolactic fermentation, LAB uses glucose to produce lactic acid through glycolysis.In heterolactic fermentation, ethanol or acetic acid are produced in addition to lactic acid [19].Lactic acid accounted for the largest proportion of organic acids in Zhalajiao.As Zhalajiao fermentation proceeded, the lactic acid content continued to increase, whereas the acetic acid content decreased, indicating that Zhalajiao fermentation is mainly homolactic fermentation.

Changes in FAAs during Fermentation
FAAs not only have nutritional functions, but they are also important flavoring substances.Asp and Glu are umami amino acids; His, Arg, Val, Phe, Leu, Met, Ile, and Tyr are slightly bitter amino acids; and Ser, Thr, Gly, and Ala are classified as sweet amino acids.These FAAs are composed of different taste-related peptides, and changes in amino acid composition may have a certain effect on aromatic components, such as alcohols and esters [13].Table 2 shows that 17 FAAs were detected during Zhalajiao fermentation.The FAA content had an initial value of 10.88 g/kg, which peaked at 16.86 g/kg on day 7, and then decreased.Similar results were reported for fermented fish-chili paste [20] and may be attributed to lactic acid fermentation, in which soluble proteins were degraded into FAAs by endogenous enzymes secreted by microorganisms in Zhalajiao during fermentation [21].The inhibition of the growth and metabolism of microorganisms with the progress of fermentation reduced the decomposition rate of proteins.At the same time, FAAs participated in the reaction for the formation of flavor substances, resulting in a reduction in the FAA content from 16.86 g/kg (D7) to 3.74 g/kg (D45).Subsequently, the FAAs content remained relatively stable at approximately 4 g/kg from days 45 to 90 of fermentation.At the beginning of fermentation, sweet amino acids dominated among the amino acids, contributing to the flavor of fresh chili.As fermentation progressed, sweet amino acids were consumed and decreased, whereas bitter, umami, and tasteless amino acids increased.Ser is the main sweet amino acid.It increased first and then decreased during fermentation.Ser was present at the highest content in the first 30 days of fermentation, accounting for 70.42-86.18% of the total content of amino acids.Different from our study, a previous work showed that Thr, followed by Asp and Ser, was the most abundant amino acid in Paojiao [13].The difference in the amino acid contents may be ascribed to the different chili varieties and fermentation environments used.After 45 days of fermentation, the content of Ser was only approximately 9%, and the content of various FAAs was relatively balanced.Additionally, the content of His in Zhalajiao before fermentation was 43.93 ± 0.02 mg/kg and decreased to almost undetectable levels after fermentation.Given that FAAs are important components in the growth of LAB, the decrease in the His content during the late fermentation period may be due to utilization by LAB [14].Note: Different letters in the same row represent significant differences (p < 0.05).nd: not detected.

Changes in Volatile Flavor Compounds during Fermentation
In this study, the volatile compounds were analyzed, and the results are shown in Supplementary Table S1.Volatile compounds play an important role as flavor determinants in the quality evaluation of fermented chili foods, which are affected by microbial communities and metabolic reactions during fermentation [22].Approximately 66 different volatile components were detected throughout the whole fermentation process of Zhalajiao.These components included 28 esters, 18 terpenes, 12 alkanes, 5 acids, 1 alcohol, 1 aldehyde, and 1 pyrazine.The content of total volatile components decreased from an original value of 36.11mg/kg to 13.30 mg/kg on day 7 and then increased continuously, reaching the maximum value of 82.04 mg/kg on day 60.Furthermore, the analysis of the OAV of the above volatiles showed that a total of 21 aroma substances, including 13 esters, 7 terpenes, and 1 pyrazine, had an OAV > 1 during Zhalajiao fermentation (Table 3).Although these substances did not have the highest content in the whole flavor system, they had a significant effect on the flavor quality of Zhalajiao [20].Ethyl caproate, (Z)-4-decenoic acid, ethyl ester, (+)-Limonene, Ocimene, Linalool, and β-Caryophyllene had an OAV > 1 on day 0, indicating that these six substances contributed to the aroma of fresh chili.Among these substances, (+)-Limonene had an OAV > 1000, which is evidence showing that it is the major aroma component of fresh chili.Zhalajiao contained 21 compounds with an OAV > 1 after 45 days of fermentation.Of these compounds, Ethyl butyrate, Ethyl caproate, Hexyl acetate, Hexanoic acid, hexyl ester, Linalool, β-Ionone, and 2-isobutyl-3-methoxypyrazine had an OAV > 100, indicating that these seven substances were the main aroma components in Zhalajiao after fermentation.B-Ionone was not detected in the first 22 days of fermentation.However, its content reached 965.86 µg/kg on day 45 of fermentation, likely due to the degradation of carotene during Zhalajiao fermentation, resulting in β-Ionone formation [23].Furthermore, given its extremely low threshold, β-Ionone contributed greatly to the aroma of Zhalajiao after fermentation.Remarkably, although alkane compounds were present at high amounts in Zhalajiao, they were not the main aroma components of Zhalajiao due to their high flavor threshold and small contribution to the flavor of Zhalajiao [22].The dynamic changes in volatile compounds in samples at different fermentation times were further explored by cluster heat mapping and principal component analysis (PCA).Cluster analysis (Figure 1A) showed that the volatile flavor compounds in Zhalajiao fermentation mainly accumulated during the 30-60-day fermentation period.During Zhalajiao fermentation, the content of esters, alcohols, aldehydes, acids, and pyrazines increased first then decreased.Esters are typically produced through the esterification of short-chain acids and alcohols and are considered to be an important source of pleasant sweet and fruity flavors in fermented foods [10].In this study, only eight types of ester compounds were detected with a content of 5.22 mg/kg on day 0. With the progress of fermentation, the types and contents of ester compounds increased continuously, reaching the maximum content of 42.89 mg/kg on day 60 of fermentation.Terpenoids are typically present in various plant cells in the form of glycosides and are one of the secondary metabolites of plants [24].Seven terpenoids were detected in Zhalajiao with a total content of 24.88 mg/kg at day 0 of fermentation.Although other terpenoids were detected after 30 days of fermentation, the total content of terpenoids at this time point was lower than that on day 0. Alcohols, aldehydes, acids, and pyrazines were detected after 30 days of fermentation at relatively low contents.Notably, the composition of volatile components in Zhalajiao differed at different fermentation times.Terpenoids accounted for 68.9% and esters and alkanes accounted for approximately 15% of the total content of volatile components on day 0 of fermentation.Among terpenoids, (+)-limonene was present at the highest content, accounting for approximately 40% of the total content.With the extension of fermentation time, the proportion of terpenoids decreased, whereas that of esters and alkanes increased.This change trend has also been reported previously [16].After 30 days of fermentation, ester compounds accounted for approximately 50% and terpenoids and alkanes accounted for approximately 20% of the total content, likely because during fermentation, alcohols and acids can be used as precursors for the synthesis of ester compounds and esterification yields esters [15].Nevertheless, fermented Zha-chili is acidic, and esters may be hydrolyzed to acids and alcohols under acidic conditions.In other words, the esterification is reversible.Such reversibility may account for the decrement in the ester content at the later stage of fermentation.
PCA was performed on the basis of the average content of all volatile compounds to visualize the dynamics of volatile compounds during fermentation.The PCA revealed that the variance contribution rate of PC1 (70.6%) and PC2 (11.9%) was 82.5%, which could well represent the overall sample information (Figure 1B).Each group had excellent repeatability, and significant discrepancies were found in the aroma characteristics of the samples, indicating that the increase in the fermentation time could dramatically affect the overall aroma characteristics of Zhalajiao during natural fermentation [25].The PCA divided the Zhalajiao samples taken at different fermentation stages into three independent groups, among which the D0, D7, D15, D22, D30, and D90 samples showed high similarity and clustered into the same category.These samples were distant from the D45 and D60 samples in terms of distribution area, suggesting significant differences in the volatile compound content between these samples and the Zhalajiao samples collected on days 45 and 60 [26].Additionally, the 21 key volatile compounds with an OAV > 1 in Zhalajiao fermentation were subjected to cluster analysis, as shown in Figure 1C.More types and higher contents of key volatiles, particularly esters and terpenes, were found on day 45 of fermentation than at previous stages (red areas of clusters in Figure 1C).Consequently, this stage can be deduced to be the key stage of flavor formation.Interestingly, (+)-limonene, ocimene, and β-caryophyllene have fruity, fragrant, and spicy properties, respectively [9], and are the major aroma compounds in erijingtiao chili before fermentation.Ethyl butyrate, ethyl caproate, hexyl acetate, hexanoic acid, hexyl ester, linalool, and β-ionone had significant contributions to the flavor of Zhalajiao on day 45 of fermentation, endowing fruity and grassy aromas [22].Furthermore, the aroma of 2-isobutyl-3-methoxypyrazine, which appeared after fermentation, is similar to that of pepper and coffee [3].Although 2osobutyl-3-methoxypyrazine was present at low contents in Zhalajiao, it had a threshold of only 0.016 µg/kg and an OAV > 1000 and was therefore the characteristic aroma component of Zhalajiao after 30 days of fermentation.

Alpha Diversity of Microorganisms
This study analyzed the alpha diversity of the microbial communities in Zhalajiao at different fermentation stages by using Illumina Miseq high-throughput sequencing technology to gain a good understanding of the characteristics and dynamic succession of microbial communities throughout Zhalajiao fermentation.A total of 940,098 valid bacterial sequences were obtained after sequence screening.The optimized sequences had an average length of 442 bp.Moreover, 892,596 valid fungal sequences, with an average length of 274 bp for the optimized sequences, were obtained.The microbial diversity in Zhajiao was analyzed at the level of 3% dissimilarity (Figure 2A).The Ace and Chao indices represent bacterial abundance.They are used to estimate the number of species in a sample and are positively correlated with microbial abundance [27].As shown in Figure 2A, during Zhalajiao fermentation, the Ace and Chao indices of bacteria were considerably higher than those of fungi, indicating that bacterial flora were far more abundant than fungal flora.The bacterial Ace and Chao indices were high at the beginning of Zhalajiao fermentation, and the abundance of bacterial flora fluctuated during fermentation.The fungal Ace and Chao indices first decreased, then increased, and continued to decrease throughout the entire fermentation process.Overall, the abundance of fungal communities in Zhalajiao decreased over 90 days of fermentation.The Shannon and Simpson indices reflect community diversity, with the Shannon index being positively related and the Simpson index being negatively related to the complexity of microbial flora [28].During the fermentation of Zhalajiao, the variation range of the bacterial Shannon index was significantly greater than that of the fungal Shannon index.At the initiation of fermentation (D0), the difference between the Shannon indices of bacteria and fungi was small.However, as fermentation progressed, the bacterial Shannon index tended to increase and then decreased and remained stable after 22 days, whereas the fungal Shannon index decreased significantly and then remained stable.The trend of the Simpson index was opposite to that of the Shannon index.These data indicated that large numbers of bacteria and fungi were present at the initial stage of Zhalajiao fermentation.Nevertheless, as a result of the development of dominant flora, such as Lactobacillus, Cyanobacteria, and Candida during fermentation, some species decreased, resulting in a reduction in microbial diversity.In general, bacterial communities showed greater diversity and richness than fungal communities.

Changes in Bacterial Community Diversity and Composition during Fermentation
A total of 60 different bacterial genera were identified during Zhalajiao fermentation.The taxonomic group with a relative abundance of more than 1% was defined as the dominant group [3] and included six species, namely, Lactobacillus, Cyanobacteria, Staphylococcus, Mitochondria, Enterobacter, and Xanthomonas.The relative abundance of bacterial communities at the genus level in Zhalajiao at different fermentation times is shown in Figure 2B.In the beginning, the abundance of Lactobacillus in Zhalajiao was only 1.99%, whereas the abundances of Cyanobacteria, Staphylococcus, and Mitochondria were 46.11%, 45.82%, and 1.95%, respectively.After 7 days of fermentation, the abundance of Lactobacillus increased rapidly to 64.23%, that of Cyanobacteria decreased to 29.65%, and that of Staphylococcus was only 2.21%.Subsequently, the abundance of Lactobacillus ranged from 64% to 72%, whereas that of Cyanobacteria ranged from 16% to 33%.As a result of the increase in oxygen consumption and acid production, Lactobacillus occupied a dominant position throughout the entire fermentation process.This result was consistent with the finding reported by Dong et al. [4].Notably, Lactobacillus is the dominant microorganism in many fermented foods, including Koumiss [29], Chinese fermented vegetables [5], and vinegar [30].Staphylococcus has been found to be the predominant genus in fermented foods, including Chinese strong-flavored liquor [31] and traditional Chinese low-salt fermented fish [32].It has been demonstrated that Staphylococci from raw milk and cheese exhibits lipolytic and proteolytic activities and creates significant amounts of flavor compounds, such as FAAs, aldehydes, amines, and free fatty acids [33].

Changes in Fungal Community Diversity and Composition during Fermentation
A total of 61 different fungal genera were identified during Zhalajiao fermentation.Five species had a relative abundance greater than 1%: Candida, Meyerozyma, Aspergillus, Gibberella, and Fusarium.The relative abundance of fungi at the genus level in Zhalajiao at different fermentation times is shown in Figure 2C.On day 0 of fermentation, the abundance of Candida was 87.46%, that of Meyerozyma was 6.9%, and that of Gibberella and Fusarium was 1%.After 7 days of fermentation, the abundance of Candida remained above 98%.Candida was the dominant fungus in the fermentation of Zhalajiao.Candida and Meyerozyma have been reported to be critical for the flavor formation of fermented chili paste [8] and pickled radishes [34].In addition, Meyerozyma is one of the main yeasts found in fermented soybean paste and koji from Japan and China [35].It can produce volatile phenolic compounds during the fermentation of Thai soy sauce and improve the flavor characteristics of Thai soy sauce [36].

Pearson Correlation Analysis between Core Microorganisms and Flavor Components during Fermentation
In this study, 17 amino acids, 6 organic acids, and 21 key volatile compounds with a crucial role in the taste and flavor of Zhalajiao were detected during fermentation.Remarkably, given that fermentation is a process based on microbial metabolism, microorganisms are closely related to the production of typical fermented flavor substances [8].The potential relationship between microbial community succession and significant changes in key flavor components during Zhalajiao fermentation was revealed by using Pearson correlation analysis (|r| > 0.70, p < 0.05).
Figure 3A,B show the correlation between the succession of dominant bacterial communities and the changes in key flavor compounds during Zhalajiao fermentation.The diagram shows that Lactobacillus was mainly positively correlated with most organic acids and FAAs but was weakly correlated with volatile compounds.Interestingly, Lactobacillus was significantly associated with lactic acid (p < 0.05).Lactic acid provides sourness and can interact with alcohols, aldehydes, and ketones during the fermentation of many vegetables to produce a variety of new flavor substances [16].The results of this study are consistent with those of Fang et al. [30], who reported that Lactobacillus was significantly positively correlated with the various flavor compounds (including lactic acid, acetic acid, ethyl ester, dodecanoic acid, and 2,4-di-tert-butylphenol) produced during the fermentation of Zhejiang rose vinegar.Furthermore, Lactobacillus has been reported to produce D-amino acids, which can enhance the taste of fermented foods [37].Enterobacter exhibited a significant positive correlation with some esters with high odor activity.By contrast, Staphylococcus was closely correlated with histidine (bitterness) (p < 0.01) likely because it can secrete bitter amino acids and peptides [38].In addition, Staphylococcus showed a significantly positive correlation with (+)-limonene and ocimene.Xanthomonas presented a significantly positive correlation with succinic acid.However, different from our study, the work of Xu et al. [8] found that Xanthomonas was positively correlated with phenol and styrene.These results indicated that Lactobacillus, Enterobacter, and Staphylococcus are the core bacteria that may play an important role in the production of flavor compounds.Rao et al. [34] also proposed that Lactobacillus and Enterobacter contributed significantly to the flavor of pickled radish.
As shown in Figure 3C,D, a significant correlation (p < 0.05) was found between 4 fungi and 17 flavor substances in Zhalajiao fermentation.The network graph shows that fungi mainly exhibited significant correlations with changes in amino and organic acids during Zhalajiao fermentation.In the fungal community, Candida was conspicuously correlated with 14 flavor substances in fermentation, suggesting that it is the main fungal genus responsible for the formation of Zhalajiao flavor.Specifically, Candida showed a high correlation with various amino acids, such as Val, Phe, Leu, Lys, Ile, Tyr, Thr, and Gly (p < 0.05).This result has also been reported by Xiao et al. [12], who pointed out that Candida plays an important role in the production of amino acids (Asp, Gly, Leu, Ile, His, Ser, and Glu) in fermented foods.Candida etchellsii have been reported to enhance the flavor of soy sauce [39].Furthermore, our results demonstrated that Meyerozyma, Gibberella, and Fusarium were highly positively correlated with His (bitter FAAs) and (+)-limonene (with citrus and lemon aroma) (p < 0.01), positively correlated with Arg, acetic acid, and ocimene, but were negatively correlated with other FAAs and organic acids.These results indicate that Candida is mainly the core fungus that produces the key flavor of Zhalajiao.

Conclusions
The correlation between the flavor components and microbial composition during Zhalajiao fermented for 90 days was elucidated.During Zhalajiao fermentation, pH decreased, the total acid content continuously increased, and the organic acid content first increased then decreased.Lactic and acetic acids were the main organic acids.The content of total FAAs decreased after fermentation.Bitter, umami, and tasteless FAAs increased, whereas sweet FAAs decreased after fermentation.In addition, 21 volatile compounds (OAV > 1), mainly esters and terpenoids, were identified as the key volatile compounds of Zhalajiao.Finally, high-throughput sequencing and correlation analysis revealed that in Zhalajiao, Lactobacillus and Candida were the most dominant bacterial and fungal genera, respectively.These genera were significantly positively correlated with most of the key flavor components and had important contributions to the flavor formation of Zhalajiao during natural fermentation.These results contribute to improving the quality of Zhalajiao and provide theoretical guidance for stabilizing the flavor of Zhalajiao.

Figure 1 .
Figure 1.Dynamics in volatile compounds of Zhalajiao at different fermentation stages.Heatmap of the volatile compounds (A), principal component analysis (B), heatmap cluster analysis of the 21 key volatile compounds (C).The heatmap shows the relative abundance of volatiles at different stages of fermentation, which is directly correlated with the color intensity.

Figure 1 .
Figure 1.Dynamics in volatile compounds of Zhalajiao at different fermentation stages.Heatmap of the volatile compounds (A), principal component analysis (B), heatmap cluster analysis of the 21 key volatile compounds (C).The heatmap shows the relative abundance of volatiles at different stages of fermentation, which is directly correlated with the color intensity.

Figure 2 .
Figure 2. Dynamics in microbial community of Zhalajiao at different fermentation stages.Alpha diversity (A); the relative abundance of bacterial community (B) and fungal community (C) at genus level.Some classes with relative abundance < 1% are regarded as "others".The Ace and Chao indices represent bacterial abundance.They are used to estimate the number of species in a sample and are positively correlated with microbial abundance [27].As shown in Figure 2A, during Zhalajiao fermentation, the Ace and Chao indices of bacteria were considerably higher than those of fungi, indicating that bacterial flora were far more abundant than fungal flora.The bacterial Ace and Chao indices were high at the beginning of Zhalajiao fermentation, and the abundance of bacterial flora fluctuated during fermentation.The fungal Ace and Chao indices first decreased, then increased, and continued to decrease throughout the entire fermentation process.Overall, the abundance

Figure 2 .
Figure 2. Dynamics in microbial community of Zhalajiao at different fermentation stages.Alpha diversity (A); the relative abundance of bacterial community (B) and fungal community (C) at genus level.Some classes with relative abundance < 1% are regarded as "others".

Figure 3 . 3 .
Figure 3.The clustering heatmap (bacteria: (A); fungi: (C)) and correlation network (bacteria: (B); fungi: (D)) between microbial community at genus level and main flavor compounds.The circular Figure 3.The clustering heatmap (bacteria: (A); fungi: (C)) and correlation network (bacteria: (B); fungi: (D)) between microbial community at genus level and main flavor compounds.The circular and square symbols indicate the preponderant genus and the key flavor components, respectively.Orange and grey lines show a positive and negative correlation, respectively.

Table 1 .
Changes in physiochemical properties during the fermentation of Zhalajiao.

Table 2 .
Changes in free amino acids (FAAs, mg/kg) during the fermentation of Zhalajiao.

Table 3 .
The average OAV values of major volatile compounds in different fermentation stages.
Note: nd: not detected.Odor thresholds were mainly obtained from an online database (http://www.leffingwell.com (accessed on 19 August 2023)) with water applied as the matrix.