Analysis of microbial diversity of traditional Chinese starters and screening of special flavor yeast for Chinese baked flatbread

Introduction: Traditional Chinese starters have unclear microbial complexity and community functions.

Methods: To screen functional yeasts for flavor enhancement of Chinese baked flatbread (CBF), Saccharomyces cerevisiae Y4 was isolated and identified from traditional Chinese starters using a combined approach of high-throughput sequencing and culturomics. The traditional starters samples were aseptically collected from Luohe, Xinxiang in Henan province and Rizhao in Shandong province.

Results: The analysis of microbial communities revealed that Luohe starter had the highest fungal (Chao1 = 271.86) and bacterial (Chao1 = 341.76) richness, with dominant genera including Saccharomyces and Lactiplantibacillus. Saccharomyces cerevisiae Y4 showed the following performance: saturated gas production at 48 h, 12% ethanol tolerance, high acid tolerance (OD value below 0.15), and 1975 mL gas production in fermented dough. CBF fermented by Y4 had a high sensory score (82.8) and 40 volatiles, including the unique 1-octen-3-ol and a prominent level of phenethyl acetate content.

Discussion: In this study, a superior yeast strain was identified, providing theoretical insights and valuable microbial resources for optimizing the flavor of CBF.

1 Introduction

Chinese baked flatbread (CBF), a regionally distinctive fermented staple food, is primarily manufactured from wheat flour through dough preparation, proofing, and baking processes. The annual production of CBF in China reached 2.64 billion pieces in 2024 indicating its substantial market presence. Since the flavor quality of CBF directly influences consumers’ purchase intention and repurchase behavior, in-depth research on its flavor profile holds significant practical significance and application value. The type of starter culture employed significantly influences the final flavor and quality attributes of CBF (1). Although traditional starters require a longer fermentation time (6–72 h), staple foods fermented with traditional Chinese starters exhibit superior aroma and taste compared to those produced using pure yeast fermentation, potentially attributed to the microbial diversity inherent in conventional cultures, and microbial communities produce richer flavor substances through synergistic effects (2).

Despite advancements in microbial isolation techniques, the microbial diversity in traditional starters and the mechanisms underlying the formation of volatile flavor substances remain unclear (3). With the deepening of research, the inherent limitations of single technical systems have been overcome by more studies, which combine traditional microbial community analysis methods with modern sequencing technologies to elucidate the microbial diversity inherent in starter microbiota (4). The combination of the two enables cross-validation of data (κ > 0.85) and cross-regional sampling, enhancing the classification accuracy, credibility, and reproducibility of statistical data on microbial diversity and revealing region-specific microbiomes (5).

High-throughput sequencing (HTS), which has the characteristic of high speed, throughput, and accuracy, provides a powerful tool for deciphering microbial diversity in traditional sourdough fermentation (3). Gas chromatography–mass spectrometry (GC–MS) has the characteristic of high sensitivity, high resolution and strong anti-interference ability, which has been used for separating, identifying, and quantifying volatile and semi-volatile compounds in complex mixtures and offers a favorable means for exploring the volatile components of traditional starters (6). Their combination can complement for the limitations of each method individually and enable targeted screening of functional strains in traditional sourdoughs.

In this study, based on these considerations, an integrated approach combining HTS and culturomics was employed to investigate the microorganisms in three traditional Chinese starters previously screened, which were collected from Luohe, Xinxiang, and Rizhao. The microbial community structures and diversity in these starters were revealed. Simultaneously, the impact on the volatile substances of CBF was explored via solid-phase microextraction-gas chromatography–mass spectrometry (SPME-GC–MS).

2 Materials and methods

2.1 Materials

Three traditional Chinese starters were collected from Yuanhui District, Luohe City (LH), Hongqi District, Xinxiang City (XX) in Henan Province, and Donggang District, Rizhao City (RZ) in Shandong Province. Wheat flour was purchased from Henan Jinyuan Grain and Oil Co., Ltd., Zhengzhou, China, and conformed to the national standard GB/T 1355-2021 for first-grade wheat flour. Highly active dry yeast (AQ) was purchased from Angel Yeast Co., Ltd., Yichang, China.

Glucose, yeast extract powder, agar powder, tryptone, magnesium sulfate, manganese sulfate, peptone, beef extract, calcium carbonate, triammonium citrate, sodium acetate, yeast extract, Tween 80, tris, hydrochloric acid, isoamyl alcohol, isopropanol, glacial acetic acid, bromocresol green, ferric chloride, calcium chloride, potassium dihydrogen phosphate, 1,2-dichlorobenzene, potassium chloride, and anhydrous ethanol were supplied by Tianjin Kemiou Chemical Reagent Co., Ltd. The 50 × TAE Buffer, agarose, DNA Marker, 2 × Power Taq PCR Master Mix, and Taq PCR Master Mix were obtained from Sangon Biotech (Shanghai) Co., Ltd. The C₇ - C₃₀ n-alkane mixture standard was acquired from Sigma-Aldrich (Shanghai) Trading Co., Ltd.

2.2 Determination of microbial diversity in traditional starters

Fungal and bacterial DNA were extracted from the starters using the E. Z. N. A™ Mag-Bind Soil DNA Kit (Omega Bio-tek, Norcross, GA, United States, Cat. No. D3370-01), and the concentration and purity of the extracted DNA were analyzed. For fungal DNA, the internal transcribed spacer (ITS) region was amplified with the extracted fungal DNA as the template using primers SSU 0817F (5’-TTAGCATGGAATAATRRAATAGGA-3′) and SSU 1196R (5’-TCTGGACCTGGTGAGTTTCC-3′). For bacterial DNA, the V3-V4 hypervariable region was amplified with the extracted bacterial DNA as the amplification template using primers 338F (5’-ACTCCTACGGGAGGCAGCAG-3′) and 806R9 (5’-GGACTACHVGGGTWTCTAAT-3′). The purified amplicons were quantified, and then sequenced on the Illumina MiSeq PE300 platform. The sequencing results were uploaded to the National Center for Biotechnology Information (NCBI) database. Usearch was used to conduct operational taxonomic unit (OTU) clustering analysis on sequences with > 97% similarity (excluding singletons), and chimeras were removed to obtain representative OTU sequences. The optimal alignment results for the sequences were selected using an alignment threshold of 90%.

2.3 Isolation, purification, and preservation of advantageous yeast and bacteria

2.3.1 Isolation and purification of yeast and lactic acid bacteria

Appropriately diluted solutions (200 μL) were spread on WL medium (for yeast) and MRS medium [for lactic acid bacteria (LAB)], which were then incubated at 30 °C and 36 °C for 2–3 days, respectively. Single colonies were randomly picked based on differences in colony morphology on the media and re-streaked for isolation. After incubation, pure single colonies of yeast and LAB were obtained.

2.3.2 Strain preservation

Single colonies were picked for activation and culture. Then, glycerol and bacterial suspension were mixed at a 1:1 ratio in cryovials, which were stored in a − 80 °C refrigerator for later use.

2.4 Yeast and bacterial DNA extraction

Yeast genomic DNA was extracted using the E. Z. N. A.® Yeast DNA Kit (Omega Bio-tek, D3370-01) following the manufacturer’s protocol: enzymatic pre-lysis with lyticase (30 min, 37 °C), bead-beating with 0.5 mm glass beads (FastPrep-24, 6.0 m·s⁻¹, 30 s × 3), followed by silica-column binding and elution in 50 μL elution buffer.

Bacterial DNA was obtained using the DNeasy UltraClean Microbial DNA Kit (Qiagen, 12,224–50) according to the manufacturer’s instructions: samples were resuspended in MicroBead solution, mechanically lysed in a bead-beater, treated with MD1/MD2 lysis buffers (65 °C, 10 min), and finally purified on silica spin filters with MD4 wash buffer and eluted in 100 μL MD5 buffer.

2.5 Gene amplification

Using the extracted yeast and bacterial DNA as templates, the yeast ITS region was amplified with primers ITS4 (5’-TCCTCCGCTGACTAATATGC-3′) and ITS5 (5’-GGAAGTAAAAGTCGTAACAAGG-3′), and the bacterial 16S rRNA gene was amplified with primers P0 (5’-GAGAGTTTGATCCTGGCTCAG-3′) and P6 (5’-CTACGGCTACCTT-3′).

PCR was conducted using a PTC Genesy 96 T thermal cycler (Xi’an Tianlong Technology Co., Ltd.). All amplifications were performed in 25.5 μL reactions with four replicates. Each reaction volume contained 1 μL of DNA template, 1 μL of Primer ITS4/P0, 1 μL of Primer ITS5/P6, 12.5 μL of 2 × Power Taq PCR Master Mix (Thermo Scientific, United States) and 10 μL of Double-distilled water.

The following cycle parameters were used for yeast: initial denaturation for 5 min at 94 °C; 30 cycles of 30 s at 94 °C (denaturation), 45 s at 52 °C (annealing), and 45 s at 72 °C (extension); and final extension for 10 min at 72 °C.

The following cycle parameters were used for bacteria: initial denaturation for 5 min at 94 °C; 30 cycles of 60 s at 94 °C (denaturation), 60 s at 52 °C (annealing), and 90 s at 72 °C (extension); and final extension for 10 min at 72 °C.

2.6 Electrophoresis detection

Electrophoresis analysis was performed using a 0.8% agarose gel with a DNA Marker as reference, under conditions of 100 V voltage, 100 mA current, and 25 min duration using an electrophoresis apparatus (Model YJ300C, Beijing Junyi Dongfang Electrophoresis Equipment Co., Ltd., Beijing, China). After staining with ethidium bromide solution for 40 min, the PCR amplification fragments were observed under an ultraviolet analyzer (UVP LLC, Upland, CA, United States) to determine their clarity, and the results were photographed and recorded.

2.7 Gene sequencing

The PCR products were sent to Sangon Biotech (Shanghai) Co., Ltd. for sequence determination. Homology search was then performed in the NCBI BLAST database to identify the species of the tested strains.

2.8 Activation and collection of microbial strains

Yeast and LAB were streaked onto YPD medium and MRS medium, separately, and incubated at 28 °C and 36 °C for 48 h. Single colonies were then picked up and inoculated into YPD liquid medium and MRS liquid medium, respectively. After 20 h of cultivation, the supernatant was discarded to obtain fresh wet cells of yeast and LAB.

2.9 Yeast fermentation capacity

The fermentation capacity of yeast was determined using the Durham tube fermentation method. Activated yeast was inoculated into YPD liquid medium containing an inverted Durham tube, followed by incubation at 28 °C for 48 h. The volume of gas in the Durham tube was observed every hour, and the gas production initiation time of the strain and the gas production status at 48 h were recorded.

2.10 Lactate, acetate, and ethanol tolerance of yeast

Yeast seed culture (100 μL) was inoculated into YPD liquid media with different acid and ethanol concentrations, followed by shaking cultivation at 28 °C and 180 r/min for 20 h. The absorbance was then measured at a wavelength of 600 nm. The absorbance value was used to represent the growth density of yeast under different acid and ethanol concentrations, thus reflecting the acid tolerance and ethanol tolerance of the yeast. The volume ratios (v/v) of lactic acid and acetic acid added were as follows: (0 + 0) %, (0.4 + 0.1) %, (0.8 + 0.2) %, (1.2 + 0.3) %, (1.6 + 0.4) %, and (2 + 0.5) %. The addition amounts of ethanol were 0, 5, 10, and 12%, respectively.

2.11 Gas production characteristics of yeast-fermented dough

Fresh wet yeast cells were inoculated into flour at a concentration of 10⁸ CFU/g for dough mixing, and the gas production properties of the dough were determined using a Chopin F3 rheological fermentation instrument (Chopin Technologies, Villeneuve-la-Garenne, France). The test conditions were set as follows: temperature 35 °C; determination duration 3 h; weight of the counterbalance 1.5 kg.

2.12 Preparation of yeast-fermented CBF

Fresh wet yeast cells were added to 200 g of wheat flour at an inoculum size of 10⁷ CFU/g, with 48% water addition for dough mixing. After forming a dough with a smooth and elastic surface, it was placed in a proofing box for the first fermentation (temperature: 30 °C; relative humidity: 85%; duration: 6 h). After proofing, wheat flour (100 g) and water (45 mL) were added for remixing into dough, which was then repeatedly rolled out to fully expel air. The dough was subsequently divided into 100 g portions, which were shaped into round dough blanks using a round dough press, and then placed in the proofing box for the second fermentation (temperature: 30 °C; relative humidity: 85%; duration: 40 min). The proofed dough blanks were then transferred to an oven for baking with the top heat set to 180 °C and the bottom heat to 200 °C, for 16 min. After baking, the products were naturally cooled at room temperature for 30 min, and subsequently stored in a sealed container for later use.

2.13 Sensory evaluation of CBFs

The sensory evaluation was carried out by 20 trained panelists (10 males and 10 females, 20–40 years old) in an isolated sensory booth (room temperature 20–22 °C). The panelists were trained according to ISO standards (ISO 8586: 2012) for sensory analysis and subjected to performance monitoring (ISO 11132: 2012). Compusense software (version 5.2, Compusense Business Avionics B. V., Canada) was used for sample distribution. Panelists finished the sensory evaluation within 10 min and rinsed their mouths with spring water between samples. The appearance, color, palatability, toughness, viscosity, smoothness, and flavor of the CBF was scored by panelists on the basis of Supplementary Table S1.

2.14 Determination of volatiles

Volatiles were extracted from CBF samples using headspace solid-phase microextraction (HS-SPME) and analyzed by Gas chromatography–mass spectrometry (GC–MS) instrument (Model 8,890-5977B, Agilent Technologies (China) Co., Ltd., China) with a DB-5 capillary column (60 m × 0.25 mm × 0.25 μm) (7).

1. Solid-phase microextraction conditions: CBF sample (2 g) was placed into a 20 mL headspace vial and internal standard solution (1 μL, 1,2-dichlorobenzene, 100 mg/L) was immediately added. Sampling was performed using an automatic headspace injection. The extraction head was aged for 30 min, the incubator was preheated at 60 °C for 10 min, the extraction time was 40 min, and the sample was desorbed at 250 °C for 300 s at the injection port.

2. GC–MS conditions: The temperature program was as follows: initial temperature 40 °C held for 3 min, then 3 °C/min increase to 60 °C, followed by 4 °C/min increase to 200 °C, and finally 5 °C/min increase to 260 °C held for 10 min; injection volume: 1 μL; flow rate 1.0 mL/min; no split; full scan mode.

3. Qualitative analysis: The mass spectra obtained using GC–MS data analysis software were compared with the NSIT. L 17 spectral library. The retention index (RI) values of the detected compounds were calculated using the retention times of n-alkanes (C₇–C₃₀) and compared with the database (https://webbook.nist.gov/).

4. Quantitative analysis: Semi-quantification of volatiles was performed using the internal standard method. Following the procedure:

${\mathrm{C}}_{\mathrm{X}}=\frac{{\mathrm{S}}_{\mathrm{X}}}{\mathrm{S}\left(\text{internal}\right)}\times \frac{\mathrm{C}\left(\text{internal}\right)\times \mathrm{V}\left(\text{internal}\right)}{\mathrm{M}\left(\text{sample}\right)}$

In Formula, C_X is the relative mass concentration of compound X in the CBF sample, μg/kg; S_X is the peak area of compound X in the CBF sample; S(internal) is the peak area of the internal standard 1,2-dichlorobenzene; C(internal) is the concentration of the internal standard 1,2-dichlorobenzene; V(internal) is the volume of the internal standard added; M(sample) is the mass of the sample, g.

2.15 Statistical analysis

All experiments were performed in triplicate, and the results were expressed as “mean ± standard deviation.” For significant difference analysis, one-way analysis of variance (ANOVA) and Duncan’s multiple range test (p < 0.05) were used, with statistical analysis conducted using IBM SPSS Statistics 25 (IBM Corporation, Armonk, NY, United States). Graphs were plotted using Origin Pro 2025 (OriginLab Corporation, Northampton, MA, United States).

3 Results and discussion

3.1 Alpha diversity analysis of traditional Chinese starters microbial communities

Alpha diversity analysis, a crucial method for evaluating microbial community complexity, reflects the species richness and evenness of samples (8). As shown in Table 1, for both fungi and bacteria in the traditional starter samples of RZ, XX, and LH, the Coverage indices were all greater than 99%, indicating that the sequencing results were reliable. The LH starter sample exhibited significantly higher species richness in both fungi (Chao1 = 271.86; ACE = 274.36) and bacteria (Chao1 = 341.76; ACE = 351.64) than the other samples. This value range is consistent with the species richness characteristics of traditional dough starters reported in existing studies which detected 490 fungal OTUs and 234 bacterial OTUs in Chinese northwest traditional jiaozi (9). The high species richness of LH may be linked to the combination of its unique raw material composition or fermentation processes and potential cross-kingdom synergistic effects (8).

Table 1

Table 1. Alpha diversity of different traditional Chinese starters.

In contrast, although the fungal richness in the XX starter was moderate, its Shannon index (2.29) and Simpson index (0.22) indicated that the fungal community possessed higher evenness and niche differentiation capacity. This aligns with the findings of which reported that traditional sourdoughs with balanced microbial evenness often establish stable metabolic interaction networks (10). Such high evenness might arise from a more complex metabolic interaction network in its fermentation environment, which is also supported by the findings that high niche differentiation in traditional fermented foods contributes to functional stability (8, 11). Giving that the Coverage indices of all samples exceeding 99%, the Sobs index (271) of the LH fungal community was almost consistent with the Chao1 estimate (271.86). This consistency not only confirms the accuracy of the sequencing data but also suggests that rare microbial groups (abundance < 0.1%) in LH were fully captured. Rare taxa, despite their low abundance, often harbor unique functional genes that may play crucial roles (9). This effective detection provides a valuable dataset to explore uncharacterized microbial functions in traditional starters.

3.2 Fungal and bacterial community structures in traditional Chinese starters

Analysis of the relative abundance of fungal community composition is recognized as a critical method for uncovering the proportional distribution of distinct fungal species within a community (12). The relative abundance of fungal communities in XX, LH, and RZ starters are shown in Figures 1a,b. At the phylum level, three identifiable fungal phyla were detected: Ascomycota, Mucoromycota, and Basidiomycota. Among these, Ascomycota had the highest relative abundance (RA: > 1%), in XX (RA: 91.11%), LH (RA: 96.18%), and RZ (RA: 91.71%) starters, respectively. At the genus level, 15 fungal genera with relative abundances greater than 1% were detected.

Figure 1

Figure 1. Fungal and bacterial community structures in traditional Chinese starters. (a) Species composition at the phylum level of fungi in the traditional Chinese starters. (b) Species composition at the genus level of fungi in the traditional Chinese starters. (c) Species composition at the phylum level of bacteria in the traditional Chinese starters. (d) Species composition at the genus level of bacteria in the traditional Chinese starters. LH, RZ, and XX represent the traditional Chinese starters from Luohe, Rizhao, and Xinxiang.

Saccharomyces, as the dominant genus within Ascomycota, exhibited extremely high relative abundances in samples XX (46.98%), LH (56.83%), and RZ (59.86%). Saccharomyces is the primary microbial group involved in dough fermentation, exerting significant influence on dough expansion height and gluten network structure (13). Furthermore, Saccharomycopsis has an excellent ability to produce β-glucosidase, acid protease, and amylase, which helps promote microbial community succession enriching the flavor of malt (14). Previous studies have shown that bread fermented with exhibits superior flavor properties and lower hardness, as the genus Wickerhamomyces anomalus can convert precursor substances into aromatic substances such as esters, acids, alcohols, and aldehydes via synthetases (15).

Moreover, Saccharomyces exhibited the highest relative abundance in RZ sample, indicating that RZ starter possessed superior fermentation capacity. In addition, the main dominant fungal genera in XX sample also include Saccharomycopsis, Pichia, Rhizomucor, and Thermomyces. In LH sample, the main dominant fungal genera were Hyphopichia (RA: 12.12%) and Wickerhamomyces, (RA: 10.30%). In RZ sample, the main dominant fungal genera were Rhizomucor and Pichia.

Notably, Aspergillus, Monascus Rhizopus, and Rhizomucor were detected in XX, LH, and RZ samples. These are important fermentation strains in mold culture, possessing strong saccharification enzyme activity, which can provide the fermentation agent with multiple enzyme systems, supplying a large amount of small-molecule monosaccharides and free amino acids in the fermentation agent for metabolic utilization by lactic acid bacteria and yeast (16).

Analysis of the relative abundance of bacterial community composition is recognized as a critical method for uncovering the proportional distribution of distinct bacterial species within a community, and the structural characteristics of dominant and rare taxa can be reflected (12). The relative abundances of bacterial communities in XX, LH, and RZ starters are shown in Figures 1c,d. At the phylum level, three bacterial phyla were detected: Firmicutes, Cyanobacteria, and Proteobacteria, and they are commonly found in traditional Chinese starters (17). At the genus level, 12 dominant bacterial genera were detected. The dominant bacterial genera in XX sample were Companilactobacillus (RA: 12.55%) and Lactiplantibacillus (RA: 6.00%). The main dominant bacterial genera in RZ sample were Latilactobacillus (RA: 19.35%), Lactiplantibacillus (RA: 12.97%), and Acetobacter (RA: 8.10%). In LH sample, the main dominant bacterial genera were Companilactobacillus (RA: 18.71%), Latilactobacillus (RA: 19.98%), and Lactiplantibacillus (RA: 8.47%). Companilactobacillus, Latilactobacillus, and Lactiplantibacillus produce lactic acid and other organic acids during dough fermentation, which enrich the flavor of fermented bread products (18). Acetobacter was the only dominant bacterial genus in RZ sample, and it oxidizes ethanol to produce acetic acid during dough fermentation and improves the texture and taste of bread products (14).

3.3 Isolation, purification, and molecular identification of yeast and bacterial strains

Based on morphological differences such as colony shape, color, size, and protrusion degree, a total of 94 yeast strains were selected from XX, LH, and RZ starters, with their colony morphologies shown in Supplementary Table S2, exhibiting seven distinct colony morphologies. Through isolation, purification, analysis, and identification, species composition within microbial communities can be recognized, relative abundances and distribution characteristics of different bacteria in the community determined, and interactions among microorganisms explored (19).

The ITS sequence is extensively applied in species identification and widely recognized as the “barcode” region for fungal classification (20). The gel electrophoresis detection results of partial yeast strains are presented in Supplementary Figure S1a. Target fragments were successfully amplified from all samples. The fragment size was consistent with the expected ITS region length (500–1,000 bp), indicating that the primer design and amplification conditions exhibited favorable specificity and effectiveness.

The molecular identification results of the yeast strains are shown in Table 2. Three yeast species were identified from the three starters (XX, LH, and RZ): Saccharomyces cerevisiae, Wickerhamomyces anomalus, and Hyphopichia burtonni. Strains Y1, Y2, Y3, Y4, and Y5 were all identified as Saccharomyces cerevisiae through molecular identification, and their colony morphology characteristics differed. This is because the colony morphology exhibited by Saccharomyces cerevisiae is related to intraspecific genetic diversity (21).

Table 2

Table 2. Molecular identification of yeast and bacteria in traditional Chinese starters.

During dough fermentation, Saccharomyces cerevisiae primarily metabolizes various sugars to produce CO₂, causing the dough to rise, and can also undergo esterification reactions through metabolic processes, enriching the flavor of fermented foods (13). Moreover, Saccharomyces cerevisiae can impart good flavor to steamed breads and produce substances that inhibit fungal growth during dough fermentation and bread storage (15). Wickerhamomyces anomalus possesses excellent capabilities for producing various glycosidases, including β-D-glucosidase and α-L-arabinofuranosidase, and is widely present in traditional sourdough and mold culture (22).

The results of gel electrophoresis detection for partial bacterial strains are presented in Supplementary Figure S1b. Fragment sizes were mainly concentrated in the range of 1,000–2000 bp. The sizes of the fragments were consistent with the typical length of bacterial 16S rRNA genes, verifying the reliability of the amplification products.

The molecular identification results for the three species of bacteria from XX, LH and RZ are shown in Table 2. A total of four LABs and one acetic acid bacterium (Acetobacter malorum) were identified. The four LABs were Furfurilactobacillus rossiae, Latilactobacillus curvatus, Lactiplantibacillus plantarum, and Companilactobacillus crustorum. Lactiplantibacillus plantarum is a common LAB in traditional Chinese starters, capable of producing large amounts of organic acids during dough fermentation, which positively improves the quality of fermented dough products (18). Acetobacter malorum is an obligate aerobic bacterium that metabolizes ethanol produced by yeast fermentation to generate acetic acid and water, positively influencing the flavor of fermented dough products made with traditional Chinese starters (23).

3.4 Functional characterization of yeast

3.4.1 Yeast fermentation capacity analysis

Glucose, maltose, sucrose, and other carbohydrates in wheat flour serve as the sole carbohydrate sources for yeast, among the most critical characteristics of high-quality yeast is dough fermentation capacity (24). The fermentation capacity of yeast strains Y1–Y7 was determined using the Durham tube fermentation method. As shown in Table 3, within 48 h of fermentation, the Durham tubes of all five Saccharomyces cerevisiae (Y1–Y5) were filled with gas, indicating that these five strains possess strong fermentation capabilities. This performance aligns with the key characteristics of Saccharomyces cerevisiae reported in traditional fermented food studies (23). In contrast, some studies isolated non-Saccharomyces yeasts from traditional Chinese dough starters, which showed rapid initial gas production (within 3 h) but lower total gas yield, highlighting the superior long-term fermentation stability of our Saccharomyces cerevisiae strains (25). Among them, Saccharomyces cerevisiae Y4 began gas production at 6 h of fermentation, with the earliest fermentation initiation time, indicating the best fermentation initiation ability.

Table 3

Table 3. The fermentation ability of yeast and volumetric and sensory scores of CBFs made from different yeasts.

3.4.2 Lactic acid and acetic acid tolerance of yeast

During dough fermentation, the fermentation characteristics of yeast are influenced by extracellular environmental factors such as high temperature, low pH, and high ethanol concentration (21). The lactic acid and acetic acid tolerance of yeast strains Y1–Y7 were determined, with yeast growth density expressed by optical density (OD). As shown in Figure 2a, with the increase of acid content, the OD values of the seven yeast strains gradually decreased, reaching a stable state at an acid content of 20% (v/v). At an acid content of 15% (v/v), the OD values of Saccharomyces cerevisiae Y4 and Y5 were significantly higher than those of Saccharomyces cerevisiae Y1, Y2, Y3, Wickerhamomyces anomalus Y6, and Hyphopichia burtonni Y7, indicating that Saccharomyces cerevisiae Y4 and Y5 had better tolerance to lactic acid and acetic acid. Moreover, differences exist in acid tolerance among Saccharomyces cerevisiae, which is consistent with previous research (23). Furthermore, the pH during dough fermentation typically ranges from 3.8 to 6.7. Therefore, all seven yeast strains can adapt to the acidic environment of dough, and Saccharomyces cerevisiae Y4 and Y5 exhibit superior acid tolerance.

Figure 2

Figure 2. (a) Acid tolerance of yeast. (b) Ethanol tolerance of yeast. (c) Gas release curves of doughs fermented with different yeast strains. (d) Ferment height curves of doughs fermented with different yeast strains. Y1–Y5 represent Saccharomyces cerevisiae Y1–Y5. Y6 represents Wickerhamomyces anomalus Y6. Y7 represents Hyphopichia burtonni Y7.

3.4.3 Ethanol tolerance of yeast

During dough fermentation, yeast produces ethanol through carbohydrate metabolism, and excessively high ethanol concentrations can reduce yeast activity and inhibit its proliferation rate (21). The ethanol tolerance of Saccharomyces cerevisiae Y1, Y2, Y3, Y4, Y5, Wickerhamomyces anomalus Y6, and Hyphopichia burtonni Y7 were determined (Figure 2b). When the ethanol volume fraction exceeded 10%, the OD values of all seven yeast strains began to decrease, indicating that the growth was inhibited. Saccharomyces cerevisiae Y1–Y5 exhibited stronger ethanol tolerance than Wickerhamomyces anomalus Y6 and Hyphopichia burtonni Y7. All seven yeast strains could tolerate an ethanol volume fraction of 10% (v/v), and the ethanol concentration in fermented dough generally did not affect their activity and proliferation rate.

3.4.4 Gas production characteristics of yeast fermentation dough

Yeast aerobic fermentation generates abundant gas, which enables the full stretching of the gluten network and retains CO₂ in the dough to form gas pockets, thereby endowing fermented dough with a soft texture (26). As shown in Figure 2c, significant differences were observed in gas productions and gas retention capacities of doughs fermented by seven yeast strains (Y1–Y7). For doughs fermented by Saccharomyces cerevisiae Y1–Y5, gas production first increased sharply within the initial 40 min of fermentation, then decreased significantly, followed by a slow increase to the maximum gas production. Furthermore, the gas production of the fermented dough of Saccharomyces cerevisiae Y4 was the highest (total gas production of 1975 mL).

As shown in Figure 2d, the expansion heights of the doughs fermented by Saccharomyces cerevisiae Y1–Y5 first increased continuously and then stabilized or decreased slightly. The expansion heights of the doughs fermented by Saccharomyces cerevisiae Y1, Y2, Y4, Y5 showed no significant differences, with maximum expansion heights of 43.2 mm, 41.9 mm, 42.3 mm, and 42.4 mm, respectively, and they were significantly greater than that of Saccharomyces cerevisiae Y3 (35.2 mm). Furthermore, the expansion heights of the doughs fermented with Wickerhamomyces anomalus Y6 and Hyphopichia burtonni Y7 continued to increase slowly. The results showed that Saccharomyces cerevisiae Y4 can quickly initiate the gas production process and significantly shorten the dough fermentation time. Moreover, it steadily maintained a high height with minimal fluctuations, which can ensure the stability of the dough structure and reduce the risk of collapse caused by a sudden drop in gas production. Based on the above results, it can be concluded that the dough fermented with Saccharomyces cerevisiae Y4 has strong gas production and retention capabilities, resulting in better dough fermentation quality.

3.5 Sensory evaluation analysis of yeast fermentation for making CBF

Sensory evaluations of CBFs made with seven yeast strains (Y1–Y7) were conducted, and AQ yeast was as the control. As shown in Table 3 and Supplementary Table S3, CBF fermented individually by the seven yeast strains showed no significant differences in surface color and surface morphology that the sensory scores of these CBFs ranged 79.20–82.80. CBFs fermented with the five Saccharomyces cerevisiae exhibited larger specific volume, better springiness, fine and uniform internal pores, soft texture, crisp taste, and no stickiness to the teeth. Saccharomyces cerevisiae Y4 exhibited remarkable advantages, and the CBF fermented by this strain showed higher volumetric ratios, indicating stronger gas-producing and dough expansion capacities. Moreover, it achieved the highest overall sensory score (82.80 ± 0.5) and performed well in multiple dimensions. In contrast, Y6 and Y7 had lower sensory scores (81.00 ± 0.8 and 79.20 ± 0.7). The sensory score (81.90 ± 0.8) of control AQ performed on par with most Saccharomyces cerevisiae strains lower than Y4.

3.6 Volatiles of CBFs made by Saccharomyces cerevisiae Y4 and yeast (AQ)

The composition of volatiles in baked food serves as the fundamental material basis for reflecting the aroma characteristics of baked food (7). As shown in Table 4, a total of 40 volatiles were detected in the CBF (Y4) sample, including alcohols (12), esters (5), aldehydes (10), acids (6), ketones (2), furan (1), and other compounds (4). CBF (Y4) and CBF (AQ) samples had 38 and 36 volatiles, respectively.

Table 4

Table 4. Compositions of volatiles in CBFs made by Saccharomyces cerevisiae Y4 and Angel yeast (AQ).

During dough fermentation, yeast metabolizes amino acids via the Ehrlich pathway to generate various alcoholic metabolites, which exert a positive impact on the herbal and wine-like aromas of fermented products (27). CBF (Y4) contained high relative contents of phenethyl alcohol, 2-methylbutanol, and isopentanol, consistent with previous research findings. Phenethyl alcohol has a rose-like aroma and is positively correlated with bread aroma, while 2-methylbutanol can add rich, fruity, and floral aromas to fermented dough products (28). Isoamyl alcohol, derived from leucine metabolism, is positively correlated with aroma and exhibits “fragrant, mellow, and malt-like” flavors (29). Compared to CBF (AQ), CBF (Y4) uniquely contained 1-octen-3-ol and trans-5-decenoic alcohol, and the former contributed mushroom and grassy aromas to add a fresh layer of complexity, while the latter further enriched the stereochemical structure of alcohol compounds (7). Moreover, relative content of isobutyl alcohol in CBF (Y4) was significantly lower than in CBF (AQ), avoiding a harsh alcoholic taste (28).

Furthermore, aldehydes are produced via the oxidation of lipids and amino acids. Among the detected aldehyde volatiles, benzaldehyde, nonanal, phenylacetaldehyde, n-hexanal, and trans-2-nonenal are the main ones, characterized by flavor profiles of nutty, floral, and fruity aromas (7). Meanwhile, acidic compounds not only produce acidity during dough fermentation but also serve as important precursors for other aromatic compounds, positively contributing to flavor formation (30). The acidic volatiles detected in CBF (Y4) primarily included caproic acid, benzoic acid, and caprylic acid. An appropriate amount of caproic acid can contribute a cheesy flavor with an oily aroma (27). Caprylic acid exhibits a mild cheesy taste and can react with ethanol to form caprylic acid ethyl ester, imparting a fruity aroma to CBF (Y4).

Esters, as important aromatic compounds, are formed via the esterase-catalyzed reaction between unsaturated fatty acids and alcohols during fermentation (31). Among ester volatiles, those with relatively high contents include ethyl caprylate, phenethyl acetate, γ-nonalactone, and ethyl decanoate. The content of ethyl phenylacetate in CBF (Y4) is 2.46 times that of CBF (AQ), contributing a fruity aroma to CBF (29). In addition, both 3-hydroxy-2-butanone and geranyl propionate were detected in two CBF samples. Ketones are mainly produced through amino acid degradation and β-oxidation of saturated fatty acids (18). 3-Hydroxy-2-butanone contributes a buttery aroma, while geranyl propionate provides a fruity aroma. 2-Pentyl furan imparts a floral aroma to CBF (Y4). The unique indole in CBF (Y4) further enhance the uniqueness of the flavor. The results showed that CBF (Y4) had abundant unique volatiles, synergistically enhancing flavor complexity and uniqueness.

4 Conclusion

Seven yeast strains were isolated and obtained from three traditional Chinese starters, and Saccharomyces cerevisiae Y1–Y5 significantly outperformed Wickerhamomyces anomalus Y6 and Hyphopichia burtonni Y7 in terms of fermentation capacity, alcohol tolerance, and gas production. Moreover, the fermentation characteristics values of Saccharomyces cerevisiae Y4 are the best among the five Saccharomyces cerevisiae strains. CBF made with Y4 exhibited more balanced metabolic characteristics in terms of the types and concentrations of key volatiles (1-octen-3-ol, trans-5-decen-1-ol, 2-phenylethanol, phenethyl acetate, and trans-2-nonenal), and it ensured a rich diversity of volatiles and avoided excessive production of irritating components. In the future, the growth and metabolic dynamics of microorganisms and the complex interactions between them during the dough fermentation process will be further investigated.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors. The original contributions presented in the study of the 16S rRNA gene sequences are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/search/, nucleotide sequence(s): PX775921-PX775937,PX775938-PX775947.

Author contributions

YY: Conceptualization, Data curation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. YW: Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. BL: Investigation, Writing – original draft. CC: Investigation, Writing – original draft. YZ: Formal analysis, Funding acquisition, Writing – original draft. LP: Project administration, Resources, Writing – original draft. RT: Resources, Writing – original draft. GZ: Formal analysis, Methodology, Resources, Software, Supervision, Writing – original draft. QZ: Conceptualization, Formal analysis, Methodology, Project administration, Resources, Supervision, Writing – original draft. HZ: Methodology, Validation, Writing – original draft. JS: Investigation, Validation, Writing – original draft.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the Young Key Teachers Training Program of Henan Province (grant no. 2023GGJS059); the Henan Province Key R&D Special Project (grant no. 231111111800); the Science and Technology Research Project of the Henan Provincial Department of Science and Technology (grant no. 242102110115); the Major Special Projects of Henan Province (grant no. 231100110300); the University Student Innovation and Entrepreneurship Training Program (grant no. 202510463087); and the Science and Technology Research Project of Henan Provincial Department of Science and Technology (grant no. 252102110119).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that Generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2025.1755466/full#supplementary-material

References

1. Wei, G, Chitrakar, B, Wu, J, and Sang, Y. Exploration of microbial profile of traditional starters and its influence on aroma profile and quality of Chinese steamed bread. J Sci Food Agric. (2023) 103:2522–31. doi: 10.1002/jsfa.12428,

PubMed Abstract | Crossref Full Text | Google Scholar

2. Lutter, L, Jõudu, I, and Andreson, H. Volatile organic compounds and their generation in sourdough. Agron Res. (2023) 21:504–36. doi: 10.15159/AR.23.017

Crossref Full Text | Google Scholar

3. Xu, J, Sun, Q, Zheng, O, Wang, Z, Wei, S, Han, Z, et al. Novel insights into the microbial succession and quality variation of golden pompano (Trachinotus ovatus) driven by storage temperatures based on high-throughput sequencing. LWT. (2024) 199:116087. doi: 10.1016/j.lwt.2024.116087

Crossref Full Text | Google Scholar

4. Michel, E, Monfort, C, Deffrasnes, M, Guezenec, S, Lhomme, E, Barret, M, et al. Characterization of relative abundance of lactic acid Bacteria species in French organic sourdough by cultural, Qpcr and Miseq high-throughput sequencing methods. Int J Food Microbiol. (2016) 239:35–43. doi: 10.1016/j.ijfoodmicro.2016.07.034,

PubMed Abstract | Crossref Full Text | Google Scholar

5. Fu, W, Rao, H, Tian, Y, and Xue, W. Bacterial composition in sourdoughs from different regions in China and the microbial potential to reduce wheat allergens. LWT. (2020) 117:108669. doi: 10.1016/j.lwt.2019.108669

Crossref Full Text | Google Scholar

6. Pu, D, Xu, Z, Sun, B, Wang, Y, Xu, J, and Zhang, Y. Advances in food aroma analysis: extraction, separation, and quantification techniques. Foods. (2025) 14:1302. doi: 10.3390/foods14081302,

PubMed Abstract | Crossref Full Text | Google Scholar

7. Wang, Y-H, Yang, Y-Y, Li, H-Q, Zhang, Q-D, Xu, F, and Li, Z-J. Characterization of aroma-active compounds in steamed breads fermented with Chinese traditional sourdough. LWT. (2021) 152:112347. doi: 10.1016/j.lwt.2021.112347

Crossref Full Text | Google Scholar

8. Zhao, Y, Sun, Z, Zhou, C, Ding, Y, Zhang, Y, Liu, Y, et al. Livestock grazing increases soil bacterial alpha-diversity and reduces microbial network complexity in a typical steppe. Agric Ecosyst Environ. (2025) 392:109753. doi: 10.1016/j.agee.2025.109753

Crossref Full Text | Google Scholar

9. Liu, T, Kang, J, Liu, L, Hu, X, Wang, X, Li, X, et al. Microbial community diversity of traditional dough starter (Jiaozi) from two provinces in Northwest China. Ann Microbiol. (2020) 70:18. doi: 10.1186/s13213-020-01544-1

Crossref Full Text | Google Scholar

10. Landis, EA, Oliverio, AM, McKenney, EA, Nichols, LM, Kfoury, N, Biango-Daniels, M, et al. The diversity and function of sourdough starter microbiomes. eLife. (2021) 10:e61644. doi: 10.7554/eLife.61644,

PubMed Abstract | Crossref Full Text | Google Scholar

11. Mandhania Madhvi, H, Paul, D, Suryavanshi Mangesh, V, Sharma, L, Chowdhury, S, Diwanay Sonal, S, et al. Diversity and succession of microbiota during fermentation of the traditional Indian food Idli. Appl Environ Microbiol. (2019) 85:e00368–19. doi: 10.1128/AEM.00368-19

Crossref Full Text | Google Scholar

12. Xiong, R, He, X, Gao, N, Li, Q, Qiu, Z, Hou, Y, et al. Soil Ph amendment alters the abundance, diversity, and composition of microbial communities in two contrasting agricultural soils. Microbiol Spectr. (2024) 12:e0416523. doi: 10.1128/spectrum.04165-23,

PubMed Abstract | Crossref Full Text | Google Scholar

13. Vandenbosch, S, De Bondt, Y, Steensels, J, Verstrepen, K, and Courtin, CM. The Thermotolerance of industrial yeast strains in a dough matrix determines their impact on the oven rise during bread baking. Food Chem. (2025) 491:145151. doi: 10.1016/j.foodchem.2025.145151,

PubMed Abstract | Crossref Full Text | Google Scholar

14. Su, C, Zhang, K-Z, Cao, X-Z, and Yang, J-G. Effects of Saccharomycopsis Fibuligera and Saccharomyces Cerevisiae inoculation on small fermentation starters in Sichuan-style Xiaoqu liquor. Food Res Int. (2020) 137:109425. doi: 10.1016/j.foodres.2020.109425,

PubMed Abstract | Crossref Full Text | Google Scholar

15. Zhou, M, and Li, Z. Characteristics of scalded dough fermented by co-cultures of Saccharomyces cerevisiae Y10, Wickerhamomyces anomalus Y13 and Torulaspora delbrueckii Y22. Int J Food Sci Technol. (2021) 56:5923–30. doi: 10.1111/ijfs.15232

Crossref Full Text | Google Scholar

16. Cao, Y, Xia, Q, Chen, J, Jin, Z, and Aniya,. Understanding of microbial diversity in three representative Qu in China and characterization of the volatile compounds in the corresponding Chinese rice wine. LWT. (2022) 164:113680. doi: 10.1016/j.lwt.2022.113680

Crossref Full Text | Google Scholar

17. Liu, X, Zhou, M, Jiaxin, C, Luo, Y, Ye, F, Jiao, S, et al. Bacterial diversity in traditional sourdough from different regions in China. LWT. (2018) 96:251–9. doi: 10.1016/j.lwt.2018.05.023

Crossref Full Text | Google Scholar

18. Li, Z, Sun, Z, Wu, Q, Huang, H, Tang, Y, Nan, J, et al. Comparative analysis of nutritional and flavor compound alterations in Cordyceps Militaris juice fermented by single or mixed Lactobacillus strains. Food Res Int. (2025) 217:116845. doi: 10.1016/j.foodres.2025.116845,

PubMed Abstract | Crossref Full Text | Google Scholar

19. Daly, SE, Feng, J, Daeschel, D, Kovac, J, and Snyder, AB. The choice of 16s Rrna gene sequence analysis impacted characterization of highly variable surface microbiota in dairy processing environments. mSystems. (2024) 9:e620-24. doi: 10.1128/msystems.00620-24,

PubMed Abstract | Crossref Full Text | Google Scholar

20. Tahir, MF, Ishtiaq, S, Anas, M, Irfan, A, Jaremko, M, Ercişli, S, et al. Integrating biotechnology into diagnostic labs: innovations in PCR and molecular procedures. Curr Res Biotechnol. (2025) 10:100325. doi: 10.1016/j.crbiot.2025.100325

Crossref Full Text | Google Scholar

21. Monnin, L, Nidelet, T, Noble, J, and Galeote, V. Insights into intraspecific diversity of central carbon metabolites in Saccharomyces Cerevisiae during wine fermentation. Food Microbiol. (2024) 121:104513. doi: 10.1016/j.fm.2024.104513,

PubMed Abstract | Crossref Full Text | Google Scholar

22. Sakandar, HA, Hussain, R, Farid Khan, Q, and Zhang, H. Functional microbiota in Chinese traditional baijiu and Mijiu Qu (starters): a review. Food Res Int. (2020) 138:109830. doi: 10.1016/j.foodres.2020.109830,

PubMed Abstract | Crossref Full Text | Google Scholar

23. Rehman, O, Wu, Y, Zhang, Q, Guo, J, Sun, C, Gao, H, et al. Acetic acid- and furfural-based adaptive evolution of Saccharomyces Cerevisiae strains for improving stress tolerance and lignocellulosic ethanol production. Chin J Chem Eng. (2024) 72:26–33. doi: 10.1016/j.cjche.2024.04.002

Crossref Full Text | Google Scholar

24. Yang, S, Zhang, T, Xiong, Y, Zhou, M, and Zhang, D. Effects of different fermentation methods on the Highland barley reconstituted dough and bread. LWT. (2025) 232:118354. doi: 10.1016/j.lwt.2025.118354

Crossref Full Text | Google Scholar

25. Teng, C, Qu, L, Shi, Y, Zhu, Y, Sun, B, and Li, X. Studies on fermentation characteristics of the yeasts screened from traditional starter cultures. J Chinese Institute Food Sci Technol. (2017) 17:74–81. doi: 10.16429/j.1009-7848.2017.11.010

Crossref Full Text | Google Scholar

26. Bu, Z, Wu, S, and Sun, X. Application of enzymes in the quality improvement of whole wheat products: effects on dietary fibre, dough rheology and gas phase. J Cereal Sci. (2025) 124:104222. doi: 10.1016/j.jcs.2025.104222

Crossref Full Text | Google Scholar

27. Fujimoto, A, Ito, K, Narushima, N, and Miyamoto, T. Identification of lactic acid Bacteria and yeasts, and characterization of food components of sourdoughs used in Japanese bakeries. J Biosci Bioeng. (2019) 127:575–81. doi: 10.1016/j.jbiosc.2018.10.014,

PubMed Abstract | Crossref Full Text | Google Scholar

28. Xu, D, Gu, M, Xi, J, Chen, L, Jin, Y, Wu, F, et al. Effect of a commercial peptidase on volatile composition and sensory properties of wheat bread crumb and crust. LWT. (2023) 189:115532. doi: 10.1016/j.lwt.2023.115532

Crossref Full Text | Google Scholar

29. Xu, Y, Zhao, J, Liu, X, Zhang, C, Zhao, Z, Li, X, et al. Flavor mystery of Chinese traditional fermented baijiu: the great contribution of Ester compounds. Food Chem. (2022) 369:130920. doi: 10.1016/j.foodchem.2021.130920,

PubMed Abstract | Crossref Full Text | Google Scholar

30. Pang, X-N, Chen, C, Huang, X-N, Yan, Y-Z, Chen, J-Y, and Han, B-Z. Influence of indigenous lactic acid bacteria on the volatile flavor profile of light-flavor baijiu. LWT. (2021) 147:111540. doi: 10.1016/j.lwt.2021.111540

Crossref Full Text | Google Scholar

31. Yang, H, Qu, T, Liu, J, Zhang, C, and Chen, M. Aroma-producing Saccharomyces Cerevisiae improved the quality of Dangshan pear (Pyrus Spp.) wine through the biosynthesis of key esters. Food Biosci. (2025) 69:106799. doi: 10.1016/j.fbio.2025.106799

Crossref Full Text | Google Scholar