- 1Technology Center, China Tobacco Henan Industrial Co., Ltd, Zhengzhou, China
- 2State Key Laboratory of Crop Stress Adaptation and Improvement, College of Agriculture, Henan University, Kaifeng, China
- 3College of Tobacco Science and Engineering, Zhengzhou University of Light Industry, Zhengzhou, China
Flavors and fragrances are essential for product quality, yet they are highly susceptible to contamination due to high moisture content and rich nutrients. This study investigates microbial growth, pH changes, volatile compound dynamics, and microbial community changes during the storage of flavors and fragrances. Results indicate that total viable counts (TVC) remained stable for the first three days but increased rapidly afterward, exceeding the acceptable limit of 5 log CFU/mL by day 7. The pH levels initially rose slightly, followed by a steady decline, which indicates spoilage progression. Gas chromatography–mass spectrometry (GC–MS) analysis revealed significant degradation of key aromatic compounds, such as 5-hydroxymethylfurfural (5-HMF), vanillin, and its derivative ethyl vanillin. Whole genome shotgun (WGS) sequencing demonstrated a marked increase in microbial community richness and diversity as storage progressed, with a notable shift in composition. Early storage stages were dominated by fungal species from the Ascomycota phylum, while later stages saw a rise in spoilage-associated bacteria, particularly from the Firmicutes and Proteobacteria phyla. Throughout the storage process, Zygosaccharomyces and its dominant species, Zygosaccharomyces bailii, remained prevalent, though their average relative abundance decreased from 81.26 to 32.29%. In addition, the bacterial species Oceanobacillus sojae and Niallia nealsonii showed significant increases in relative abundance, suggesting that bacteria were one of the key contributors to the spoilage of flavors and fragrances. Functional analysis based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database indicated a shift in metabolic pathways within the microbial community, with heightened metabolic activity correlating with spoilage. These findings provide valuable insights for improving storage methods and quality control of flavors and fragrances.
1 Introduction
Flavors and fragrances are intricate mixtures of aromatic compounds that provide distinctive and pleasant scents. In 2024, the global market for flavors and fragrances reached approximately $28.72 billion, with widespread applications across food, beverages, tobacco, perfumes, and cosmetics fields (Flavors and Fragrances Market, 2024; Adeboye et al., 2014; Zhu and Xiao, 2022). The flavor and aroma of food are critical determinants of consumer preferences, and even small additions of flavorings can significantly enhance product quality and customer satisfaction, ultimately driving consumption (Shukla et al., 2019).
Most food products naturally harbor bacteria unless they undergo specific sterilization (Bloomfield et al., 2023). Irradiation has been demonstrated to effectively inhibit the growth of harmful microorganisms, such as Staphylococcus aureus (Nortjé et al., 2006), Salmonella (Calle et al., 2021), and Bacillus cereus (Hong et al., 2008), without compromising food quality. According to Chinese national standards (GB/T 18526.4-2001), the total bacterial count in irradiated products should not exceed 1.5 × 104 CFU/mL. However, due to the high moisture content and organic composition, flavors and fragrances are particularly susceptible to microbial proliferation, which poses storage challenges and leads to spoilage, resulting in economic losses and resource waste. Spoilage in flavors and fragrances typically results from oxidation from air exposure, the development of off-flavors and off-odors, as well as microbiological contamination (Koutsoumanis et al., 2021). Of these, microbiological spoilage is the most prevalent, with bacteria, yeasts, and molds playing major roles. Currently, flavor quality evaluation heavily relies on subjective sensory assessments, lacking objective microbial monitoring techniques.
In recent years, microbial spoilage has emerged as a major concern in the industry, with particular attention to the dynamics of spoilage bacteria during storage (Xia et al., 2023). Microorganisms are the primary contributors to spoilage in most seafood products. However, only a few specific microbial communities can proliferate extensively, ultimately leading to spoilage, which are referred to as specific spoilage organisms (SSOs) (Gram and Dalgaard, 2002). Understanding these spoilage microorganisms across various foods is crucial for controlling microbial growth and extending shelf life. Traditional molecular biology methods, such as polymerase chain reaction (PCR) and 16S rRNA gene analysis, have proven insufficient for comprehensive studies of microbial dynamics in food (Bekaert et al., 2015; Ezeokoli et al., 2016; Olofsson et al., 2007). High-throughput sequencing (HTS) technologies, such as the Illumina NovaSeq platform, have emerged as a powerful tool for accurately characterizing SSOs in food products. For example, Shewanella spp. is recognized as a dominant spoilage organism due to its significant amine metabolism capabilities (Zhang et al., 2022; Zhu et al., 2016). Furthermore, HTS has also been widely applied to examine the predominant spoilage microorganisms under various storage conditions. Studies have shown that Staphylococcus and Serratia are key spoilage genera in modified atmosphere and vacuum-packaged braised chicken, respectively (Lei et al., 2022). In Daokou braised chickens, predominant spoilage bacteria include Enterococcus, Serratia, Vagococcus, and Lactobacillus (Zhang et al., 2023), and Pseudomonas is recognized as the primary spoilage organism in Spanish mackerel (Zheng et al., 2020). HTS enables detailed analysis of dominant microbial characteristics, such as their aerobic or anaerobic nature, which can provide insights into strategies for extending shelf life (Zhu et al., 2021). Additionally, functional predictions based on the KEGG database can further elucidate the potential metabolic functions of microbial communities, highlighting their impact on spoilage mechanisms and the stability of flavors and fragrances.
Despite the application of HTS in studying SSOs across various products, research on the microorganisms responsible for spoilage in flavors and fragrances remains limited. This study aims to investigate the total viable count (TVC) and pH changes of flavors and fragrances during storage. Gas chromatography–mass spectrometry (GC–MS) is employed to analyze the effects of storage on volatile compounds. Furthermore, whole-genome sequencing (WGS) is utilized to examine microbial community diversity, composition, and function throughout storage, ultimately providing a theoretical foundation for extending the shelf life of flavors and fragrances.
2 Materials and methods
2.1 Sample collection and storage
Flavors and fragrances samples were provided by China Tobacco Henan Industrial Co., Ltd. To simulate typical production process, we opened fully packaged samples, took portions, and then resealed and transported them to the laboratory. In the lab, the samples were divided into two groups: three were stored at −80°C as the initial storage group (Group A), and the remaining three were stored at room temperature (25°C, 60% relative humidity) for 14 days as the end of storage group (Group B). Additionally, samples were collected on days 0, 3, 7, 10, 14, and 17 for total viable count (TVC) and pH values measurement. Volatile compounds analysis was conducted on samples collected on days 0, 7, and 14.
2.2 Total viable count (TVC) and pH measurement
TVC was measured according to the China National Food Safety Standard Method for Food Microbiological Examination (GB 4789.2-2022). Briefly, 25 mL of sample was mixed with 225 mL of sterile saline solution to prepare a 1:10 dilution. This mixture was then serially diluted in ten-fold increments. For each dilution, 1 mL was transferred to a sterile Petri dish, followed by the addition of 15 mL of plate count agar. After thorough mixing, the plates were incubated at 37°C for 48 h, and colony counts were expressed as CFU/mL.
For pH values measurements, a Mettler Toledo Seven2GO™ pH meter was used, following the manufacturer’s instructions. Measurements were taken at regular intervals to monitor pH changes throughout the storage period at room temperature.
2.3 Volatile compounds analysis
Volatile compounds in flavors and fragrances were analyzed using gas chromatography–mass spectrometry (GC–MS). The 2 mL sample of flavors and fragrances was placed in a centrifuge tube, followed by the addition of 10 mL PBS. The mixture was vortexed for 5 min at 2,500 rpm. Next, 10 mL of HPLC-grade acetonitrile and 50 μL of internal standard (2,3-dichlorotoluene, 0.2096 g/L) were added, and the mixture was vortexed for 20 min at 2500 rpm. After freezing the sample at −20°C for 30 min, 1.5 g NaCl and 6 g MgSO₄ were added followed by quick shaking. Then, 5 mL of HPLC-grade dichloromethane was added, vortexed for 10 min at 2,500 rpm, and centrifuged for 3 min at 8,000 rpm. The supernatant was collected, concentrated, filtered, and subsequently analyzed using GC–MS.
The chromatographic conditions were as follows: an HP-5MS capillary column (60 m × 250 μm × 0.25 μm) was used. The inlet temperature was set to 240°C, with a 1.0 μL injection volume and helium as the carrier gas at a flow rate of 1.0 mL/min in splitless mode. The temperature program was as follows: hold at 50°C for 4 min, increase by 3°C/min to 70°C and hold for 5 min, then increase by 2°C/min to 100°C and hold for 17 min, followed by an increase to 120°C (hold for 10 min), and finally increase by 6°C/min to 280°C.
For mass spectrometry, electron ionization (EI) was used with an electron energy of 70 eV. The solvent delay was set to 10 min, and full-scan detection was performed with a mass range of m/z 35–500. Compound identification was achieved by comparing the retention times with known flavors and fragrances components and using the NIST20 spectral library. Quantification was performed using the internal standard method for compounds with available standards, while those without standards were quantified based on relative peak areas (compound peak area to internal standard peak area).
2.4 Whole genome shotgun (WGS) analysis
Microbial genomic DNA was extracted using the OMEGA Mag-Bind Soil DNA Kit (M5635-02) following the manufacturer’s instructions. The extracted DNA was stored at −20°C. DNA quality and quantity were measured using agarose gel electrophoresis and a Qubit™ 4 Fluorometer. Metagenome shotgun sequencing libraries were prepared using the Illumina TruSeq Nano DNA LT Library Preparation Kit with an insert size of 400 bp. Sequencing was conducted on the Illumina NovaSeq platform (PE150) at Personal Biotechnology Co., Ltd. (Shanghai, China).
Raw sequencing reads were processed to generate high-quality reads. First, sequencing adapters were removed using Cutadapt (v1.2.1), followed by trimming of low-quality reads using the sliding-window algorithm in fastp. Taxonomic classifications of reads were performed using Kraken2 against a RefSeq-derived database (including archaea, bacteria, viruses, fungi, protozoans, metazoans, and plants). Reads assigned to metazoans and plants were excluded from downstream analysis. Alternatively, Kaiju with greedy-5 mode was used against an NR-derived protein database including similar taxa.
Assembly of quality-filtered reads was conducted with Megahit (v1.1.2) using meta-large preset parameters. Contigs longer than 300 bp were clustered using mmseqs2 with an identity threshold of 95 and 90% coverage of the shorter contig. Taxonomic assignment was performed by aligning contigs against the NCBI-NT database using mmseqs2 with the taxonomy mode. Contigs assigned to plants or metazoans were excluded. Gene prediction was performed using MetaGeneMark, and coding sequences (CDS) were clustered using mmseqs2 with 90% identity and coverage.
To assess gene abundances, high-quality reads were mapped to predicted genes using Salmon in quasi-mapping mode. Abundance was normalized using CPM (copies per kilobase per million mapped reads). Functional annotation was performed using mmseqs2 against Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Linear discriminant analysis effect size (LEfSe) was applied to detect differentially abundant taxa and functions across groups. Beta diversity analysis was conducted using Bray-Curtis distance metrics and visualized through principal coordinate analysis (PCoA).
3 Results and discussion
3.1 Total viable counts and pH changes during storage
Microbial activity significantly influences the shelf life of flavors and fragrances, with total viable counts (TVC) serving as a direct indicator of microbial growth during storage. As shown in Figure 1A, the initial TVC was 2 log CFU/mL, indicating minimal contamination during production and packaging. For the first 3 days, TVC remained stable at 2–3 log CFU/mL, with no significant changes observed. However, after day 3, a rapid increase occurred, with TVC reaching 11 log CFU/mL by day 14. Beyond this point, microbial growth plateaued, likely due to the depletion of available organic nutrients, which limited further proliferation. According to the National Standard of China (GB 30616-2020), the acceptable TVC level is 5 log CFU/mL. Based on this standard, the sample remained within acceptable limits for up to 7 days.
Figure 1. Changes in total viable count (TVC) (A) and pH values (B) during the storage of flavors and fragrances.
In addition to TVC, monitoring pH values changes also provides valuable insights into spoilage progression. As shown in Figure 1B, the pH rose slightly from 3.69 to 3.74 early in storage, followed by a steady decline to 3.39 by day 14. This initial rise may be due to the production of alkaline compounds during early microbial metabolism. As storage continued, the accumulation of acidic by-products, particularly from acid-producing bacteria, became more pronounced, leading to a decline in pH as spoilage advanced. This trend is consistent with findings by Zhang et al. (2023), who observed a similar pH pattern in the spoilage of braised chicken, with an initial increase followed by a decline.
3.2 Volatile components analysis
Aroma loss significantly impacts the quality of flavors and fragrances during storage (Obenland et al., 2011). In this study, gas chromatography–mass spectrometry (GC–MS) was employed to analyze dynamic changes in volatile compounds over time. Volatile compounds with initial concentrations of no less than 1 μg/mL were listed in Table 1. The internal standard, 2,3-dichlorotoluene, showed consistent concentrations across three measurements, confirming the instrument’s reproducibility.
A total of 3 aldehydes, 12 esters, 3 phenols, 3 acids, 4 alcohols, and 4 ketones were identified. Among the aldehydes, 5-hydroxymethylfurfural (5-HMF), vanillin, and its derivative ethyl vanillin are notable aroma contributors. However, these compounds were not detected in the later stages of storage. Vanillin can undergo biotransformation into vanillic acid and vanillyl alcohol, facilitated by microorganisms like Saccharomyces or Cystobasidium larynges (Rönnander et al., 2018). By day 7, vanillin had completely degraded, while vanillyl alcohol, initially undetectable, increased to 0.0102 mg/mL by day 14. This suggested that vanillyl alcohol might be the degradation product of vanillin. Both vanillin and ethyl vanillin are widely used flavoring agents with a characteristic rich vanilla scent, with ethyl vanillin providing a more intense and long-lasting aroma. Notably, in addition to their aromatic properties, these compounds have demonstrated antimicrobial activity (Fitzgerald et al., 2003). Vanillin, for instance, has been shown to inhibit the growth of bacteria and fungi such as Escherichia coli and Saccharomyces cerevisiae, and to exert cytotoxic effects on mammalian cells (Adeboye et al., 2014; Bezerra et al., 2016; Cheng et al., 2007). Therefore, the presence of vanillin in samples contributes not only to aroma but also to the preservation of flavors and fragrances. Its degradation diminishes both the aroma profile and the preservative potential.
Ester compounds play a pivotal role in contributing fruity, floral, or sweet aromas, making them essential for enhancing the complexity and longevity of fragrances. For example, phenethyl isobutyrate and phenethyl acetate are widely used in fine perfumes, shampoos, soaps, and non-cosmetic products such as detergents and household cleaners (McGinty et al., 2012a; McGinty et al., 2012b). Sotolon is another significant compound, recognized for contributing to the characteristic roasted chicory aroma (Wu and Cadwallader, 2019). Additionally, palmarosa essential oil contains several key aromatic components, with benzyl benzoate being one of the most significant (Kou et al., 2023). In the later stages of storage, a noticeable decline in the concentrations of several esters, including phenethyl isobutyrate, DL-malic acid diethyl ester, phenethyl acetate, ethyl glycolate, and benzyldimethylcarbinyl butyrate, was observed. Notably, the concentration of phenethyl isobutyrate decreased by 77.97%. Conversely, some esters remained relatively stable or even showed slight increases, such as benzyl benzoate, ethyl L(-)-lactate, sotolon, bis (2-ethylhexyl) phthalate, ethyl linoleate, dibutyl terephthalate, and 2-hydroxypropylacetate. This variation is likely due to complex biochemical reactions that occur during storage, leading to the production of these compounds.
Sorbic acid, along with its sodium, potassium, and calcium salts, is often used in food preservation for its antifungal and antibacterial properties, working by disrupting microbial metabolism (Brul and Coote, 1999; de Jesus et al., 2021). When combined with vanillin, potassium sorbate has been shown to have a synergistic preservative effect (Kimani et al., 2023; Matamoros-León et al., 1999). In this study, a decline in sorbic acid concentration was also observed in the later stages of storage, further reducing its preservative efficacy. Whether sorbic acid and vanillin exhibit a similar synergistic effect in this study remains to be confirmed through further experimentation. Additionally, with prolonged storage, fluoroacetone levels were detected to gradually increase, raising concerns about the safety and quality of the products.
3.3 Whole genome shotgun (WGS) analysis
3.3.1 Diversity analysis
The microbial diversity in flavors and fragrances during storage was analyzed using high-throughput sequencing. Alpha diversity, which assesses microbial richness and diversity within individual samples, was evaluated. As shown in the boxplot in Figure 2A, the Goods coverage index for all six samples exceeded 0.999, indicating that the sequencing depth was sufficient to capture the microbial diversity in each sample. Species richness was assessed using the Chao1, ACE, and observed species indices. By the end of storage (Group B), these indices were significantly higher compared to those at the early storage stage (Group A) (p < 0.05), demonstrating a substantial increase in microbial richness over time. Community diversity, measured by the Shannon index, also significantly increased in Group B compared to Group A, suggesting that the microbial community had become more diverse as storage progressed. The rarefaction curve based on the observed species showed a clear trend toward a saturation plateau, confirming that the sequencing adequately reflected community richness, with Group B showing greater diversity (Figure 2B).
Figure 2. Analysis of microbial community diversity. (A) Alpha diversity of microbial communities. (B) Rarefaction curves of microbial communities. (C) Principal coordinate analysis (PCoA) of microbial communities based on Bray-Curtis distances. The horizontal and vertical axes represent relative distances without specific practical meaning. (D) Hierarchical cluster analysis (HCA) of microbial communities. In panels B, C, and D, Group A represents the early stage of storage, and Group B represents the end of storage.
Further analysis of beta diversity revealed significant differences in microbial community structure between the two groups. Principal Coordinates Analysis (PCoA), based on Bray-Curtis distances, indicated that the first principal component (PCo1) accounted for 58.6% of the total variation, while the second component (PCo2) explained 34.5%. The two groups were separated along the PCo1 axis, illustrating a distinct variation in community structure between them (Figure 2C). Hierarchical clustering analysis (HCA), based on Bray-Curtis distances and the top 10 genera by average relative abundance, depicted the similarity among the samples as a dendrogram, with Group A and Group B forming distinct clusters (Figure 2D). The PCoA and HCA results confirmed that the microbial community structures in Group A and Group B were distinctly different.
3.3.2 Taxonomic composition of microbial communities
To investigate the impact of storage on microbial communities within flavors and fragrances, we conducted species annotation analysis across various taxonomic levels for both sample groups. The results show the top 15 microbial taxa by average relative abundance at the phylum and genus levels (Figures 3A,B).
Figure 3. Composition analysis of microbial communities. (A) Relative abundance of the top 15 microbial phyla. (B) Relative abundance of the top 15 microbial genera. (C) Relative abundance of the top 15 microbial species. (D) Cladogram showing biomarkers identified through linear discriminant analysis (LDA) with a threshold score of 5. (E) LDA of significantly different microbial taxa with a threshold score of 5.
At the phylum level (Figure 3A), Ascomycota initially dominated the microbial community in the early storage stage (Group A), with an average relative abundance of 97.86%, while other phyla had abundances below 1%. By the end of storage (Group B), microbial diversity increased significantly. Although Ascomycota remained the most prevalent phylum, its relative abundance decreased to 40.60%. Meanwhile, Firmicutes increased from 0.098% to 35.14%, Proteobacteria from 0.01% to 7.54%, and Firmicutes A from 0.03% to 4.72%. These shifts indicate a compositional change in the microbial community over time, with bacteria such as Firmicutes and Proteobacteria emerging as dominant by the end of storage. Firmicutes and Proteobacteria are commonly associated with spoilage in various foods, including fermented sausages (Zhang et al., 2021), braised chicken (Lei et al., 2022) and ready-to-eat chicken meat (Qiu et al., 2022).
At the genus level (Figure 3B), Zygosaccharomyces and Saccharomyces were predominant in Group A, with average relative abundances of 81.26% and 16.31%, respectively. In Group B, genera with average relative abundances exceeding 5% included Zygosaccharomyces (32.29%), Oceanobacillus (7.59%), Niallia (6.89%), and Paenibacillus (5.21%). Among the top 15 genera, fungal genera including Zygosaccharomyces, Saccharomyces, Monascus, and Rachicladosporium showed a marked decrease from 97.57% in the early stage to 37.06% by the end of storage. Both Zygosaccharomyces and Saccharomyces belong to the Ascomycota phylum and the Saccharomycetaceae family, and are well-known for their ability to thrive in extreme conditions, such as high sugar and salt concentrations. This resilience makes them common contaminants in products with high osmotic pressure, like syrups and soy sauce. The rise of bacterial genera such as Oceanobacillus, Niallia, and Paenibacillus suggests that these bacteria are significant contributors to spoilage during storage.
At the species level (Figure 3C), Zygosaccharomyces bailii, Saccharomyces cerevisiae, Zygosaccharomyces parabailii, and Brachyspira catarrhinii accounted for nearly all microbial populations in Group A. Among them, the fungal species Zygosaccharomyces bailii was dominant, with an average relative abundance of 76.40%. As storage progressed, the relative abundance of Zygosaccharomyces bailii decreased significantly to 10.57%, although it remained the predominant species. In contrast, the bacterial species Oceanobacillus sojae and Niallia nealsonii showed substantial increases, with their average relative abundances rising from 0.05% to 7.28% and from 0.03% to 5.00%, respectively.
We used Linear Discriminant Analysis Effect Size (LEfSe) to identify microbial biomarkers that differed significantly between the two groups (Zhu et al., 2021). A higher Linear Discriminant Analysis (LDA) score indicates a greater contribution of a species to the difference between groups. Using a Wilcoxon test and an LDA score threshold of 5, five phyla were found with higher relative abundances in Group B: Cyanobacteria, Firmicutes, Firmicutes A, Proteobacteria, and Ascomycota (Figures 3D,E). Among these, only Ascomycota is fungal, while the others are bacterial.
3.3.3 Functional analysis of microbial communities
We analyzed the functional differences between samples based on Bray-Curtis distance. Principal Coordinate Analysis (PCoA) indicated that the first axis (PCo1) explained 88.4% of the total variance, while the second axis (PCo2) accounted for 8.3%, separating the two sample groups along the PCo1 axis (Figure 4A). Hierarchical clustering analysis (HCA), using Bray-Curtis distances and KEGG Level 2 categories, also confirmed distinct functional differences between the groups, with Group A and Group B forming separate clusters (Figure 4B).
Figure 4. Functional analysis of microbial communities. (A) Beta diversity of microbial community functions, presented through PCoA. (B) HCA of microbial functions based on Bray-Curtis distances. (C) LDA of significantly different functions with a threshold score of 3.8.
Further functional analysis via LEfSe, using an LDA score threshold of 3.8, identified significant KEGG pathway biomarkers between the groups. Group A was enriched in pathways related to genetic information processing, organismal systems, and human diseases at KEGG Level 1. In contrast, Group B was enriched in metabolism, suggesting that microbial metabolic activity is a primary driver of spoilage in flavors and fragrances during storage. Detailed analysis of Group B’s KEGG Level 2 pathways identified significant biomarkers in amino acid metabolism, metabolism of cofactors and vitamins, and xenobiotics biodegradation and metabolism. Additionally, cell motility emerged as a key functional biomarker for Group B (Figure 4C). These findings emphasize the functional shifts in microbial communities during storage, with metabolic activity playing a crucial role in the spoilage of flavors and fragrances.
4 Conclusion
This study investigated the relationship between microbial activity, volatile compound degradation, and the quality of flavors and fragrances during storage. The findings indicated that microbial growth, evidenced by increasing total viable counts (TVC) and a consistent decline in pH, significantly affected product quality. By day 7, TVC levels exceeded the acceptable limit, indicating spoilage. Simultaneously, key volatile compounds, such as ethyl vanillin, 5-hydroxymethylfurfural (5-HMF), and vanillin, degraded over time, resulting in diminished aroma and reduced preservative efficacy.
Metagenomic analysis revealed a notable increase in microbial richness and diversity, with a significant shift in community structure as storage progressed. While fungi dominated the microbial community during the early stages of storage, bacterial relative abundance increased substantially in the later stages. Functional analysis further confirmed that microbial metabolic activity played a critical role in driving spoilage, highlighting the importance of controlling microbial metabolism to extend the shelf life of flavors and fragrances. Although yeasts are generally considered beneficial microorganisms, certain spoilage yeasts, particularly those belonging to the Zygosaccharomyces genus, are major contaminants in the food and beverage industry (Fleet, 2011). In this study, Zygosaccharomyces bailii remained the dominant species throughout storage. This yeast is recognized as one of the most aggressive spoilage microorganisms, frequently contaminating products such as wine, fruit juices, carbonated soft drinks, and canned foods, especially those with high sugar content (Kuanyshev et al., 2017). Its high resistance to preservatives further complicates contamination control (Karaman et al., 2016). Interestingly, studies have shown that mint essential oil exhibits significant anti-yeast activity against Zygosaccharomyces bailii in apple juice. However, the combination of sodium benzoate and mint essential oil unexpectedly promoted Zygosaccharomyce bailii growth, demonstrating an antagonistic effect (Karaman et al., 2016). Additionally, some commercially available fermentates have shown inhibitory effects on Zygosaccharomyce bailii. For instance, a 1% concentration of Polyfence™ 35TD, a blend of cultured citrus extracts, significantly inhibited Zygosaccharomyce bailii and could serve as a potential replacement for potassium sorbate in dressings without affecting sensory quality (Samapundo et al., 2017).
In summary, this study provides valuable insights into the biological mechanisms underlying spoilage in flavors and fragrances during storage. It offers a theoretical foundation for developing strategies to prevent spoilage and extend product shelf life. Future research should focus on inhibiting the growth of dominant microbial species identified through high-throughput sequencing to improve product stability and longevity.
Data availability statement
The raw sequencing data of microbial communities related to flavors and fragrances have been deposited in the NCBI Sequence Read Archive (SRA) under BioProject accession number: PRJNA1191206.
Author contributions
YF: Conceptualization, Formal analysis, Funding acquisition, Methodology, Visualization, Writing – original draft, Writing – review & editing. TZ: Formal analysis, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. JY: Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. WL: Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. YY: Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. JH: Data curation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. SH: Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. ZY: Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. QL: Data curation, Formal analysis, Writing – original draft, Writing – review & editing. WZ: Data curation, Formal analysis, Writing – original draft, Writing – review & editing. QZ: Conceptualization, Methodology, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Henan China Tobacco Industry Co., Ltd. Science and Technology Project [grant number: AW2022017].
Conflict of interest
YF, TZ, JY, WL, YY, ZY, QL, WZ, and QZ were employed by China Tobacco Henan Industrial Co., Ltd.
The remaining authors declare that the research 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 authors declare that no Gen AI was used in the creation of this manuscript.
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.
References
Adeboye, P. T., Bettiga, M., and Olsson, L. (2014). The chemical nature of phenolic compounds determines their toxicity and induces distinct physiological responses in Saccharomyces cerevisiae in lignocellulose hydrolysates. AMB Express 4:46. doi: 10.1186/s13568-014-0046-7
Bekaert, K., Devriese, L., Maes, S., and Robbens, J. (2015). Characterization of the dominant bacterial communities during storage of Norway lobster and Norway lobster tails (Nephrops norvegicus) based on 16S rDNA analysis by PCR-DGGE. Food Microbiol. 46, 132–138. doi: 10.1016/j.fm.2014.06.022
Bezerra, D. P., Soares, A. K. N., and de Sousa, D. P. (2016). Overview of the role of vanillin on redox status and Cancer development. Oxidative Med. Cell. Longev. 2016:9734816. doi: 10.1155/2016/9734816
Bloomfield, S. J., Zomer, A. L., O'Grady, J., Kay, G. L., Wain, J., Janecko, N., et al. (2023). Determination and quantification of microbial communities and antimicrobial resistance on food through host DNA-depleted metagenomics. Food Microbiol. 110:104162. doi: 10.1016/j.fm.2022.104162
Brul, S., and Coote, P. (1999). Preservative agents in foods: mode of action and microbial resistance mechanisms. Int. J. Food Microbiol. 50, 1–17. doi: 10.1016/S0168-1605(99)00072-0
Calle, A., Fernandez, M., Montoya, B., Schmidt, M., and Thompson, J. (2021). UV-C LED irradiation reduces Salmonella on chicken and food contact surfaces. Food Secur. 10:1459. doi: 10.3390/foods10071459
Cheng, W.-Y., Hsiang, C.-Y., Bau, D.-T., Chen, J.-C., Shen, W.-S., Li, C.-C., et al. (2007). Microarray analysis of vanillin-regulated gene expression profile in human hepatocarcinoma cells. Pharmacol. Res. 56, 474–482. doi: 10.1016/j.phrs.2007.09.009
de Jesus, J. H. F., Szilágyi, I. M., Regdon, G., and Cavalheiro, E. T. G. (2021). Thermal behavior of food preservative sorbic acid and its derivates. Food Chem. 337:127770. doi: 10.1016/j.foodchem.2020.127770
Ezeokoli, O. T., Gupta, A. K., Mienie, C., Popoola, T. O. S., and Bezuidenhout, C. C. (2016). PCR-denaturing gradient gel electrophoresis analysis of microbial community in soy-daddawa, a Nigerian fermented soybean (Glycine max (L.) Merr.) condiment. Int. J. Food Microbiol. 220, 58–62. doi: 10.1016/j.ijfoodmicro.2016.01.003
Fitzgerald, D. J., Stratford, M., and Narbad, A. (2003). Analysis of the inhibition of food spoilage yeasts by vanillin. Int. J. Food Microbiol. 86, 113–122. doi: 10.1016/S0168-1605(03)00059-X
Flavors and Fragrances Market. (2024) Flavors and fragrances market: global industry trends, share, size, growth, opportunity and forecast 2024–2032.
Fleet, G. H. (2011). “Yeast spoilage of foods and beverages” in The yeasts. eds. C. P. Kurtzman, J. W. Fell, and T. Boekhout. Fifth ed (London: Elsevier), 53–63.
Gram, L., and Dalgaard, P. (2002). Fish spoilage bacteria – problems and solutions. Curr. Opin. Biotechnol. 13, 262–266. doi: 10.1016/S0958-1669(02)00309-9
Hong, Y.-H., Park, J.-Y., Park, J.-H., Chung, M.-S., Kwon, K.-S., Chung, K., et al. (2008). Inactivation of Enterobacter sakazakii, Bacillus cereus, and Salmonella typhimurium in powdered weaning food by electron-beam irradiation. Radiat. Phys. Chem. 77, 1097–1100. doi: 10.1016/j.radphyschem.2008.05.004
Karaman, K., Sagdic, O., and Yilmaz, M. T. (2016). Multiple response surface optimization for effects of processing parameters on physicochemical and bioactive properties of apple juice inoculated with Zygosaccharomyces rouxii and Zygosaccharomyces bailii. LWT Food Sci. Technol. 69, 258–272. doi: 10.1016/j.lwt.2016.01.006
Kimani, B. G., Takó, M., Veres, C., Krisch, J., Papp, T., Kerekes, E. B., et al. (2023). Activity of binary combinations of natural phenolics and synthetic food preservatives against food spoilage yeasts. Food Secur. 12:1338. doi: 10.3390/foods12061338
Kou, Z., Zhang, J., Lan, Q., Liu, L., Su, X., Islam, R., et al. (2023). Antifungal activity and mechanism of palmarosa essential oil against pathogen Botrytis cinerea in the postharvest onions. J. Appl. Microbiol. 134. doi: 10.1093/jambio/lxad290
Koutsoumanis, K., Tsaloumi, S., Aspridou, Z., Tassou, C., and Gougouli, M. (2021). Application of quantitative microbiological risk assessment (QMRA) to food spoilage: principles and methodology. Trends Food Sci. Technol. 114, 189–197. doi: 10.1016/j.tifs.2021.05.011
Kuanyshev, N., Adamo, G. M., Porro, D., and Branduardi, P. (2017). The spoilage yeast Zygosaccharomyces bailii: foe or friend? Yeast 34, 359–370. doi: 10.1002/yea.3238
Lei, Y., Huang, J., Cheng, Y., Zhang, Y., Huang, T., and Huang, M. (2022). Changes in bacterial communities and the volatilome of braised chicken with different packaging stored at 4 °C. Food Res. Int. 155:111056. doi: 10.1016/j.foodres.2022.111056
Matamoros-León, B., Argaiz, A., and López-Malo, A. (1999). Individual and combined effects of vanillin and potassium sorbate on Penicillium digitatum, Penicillium glabrum, and Penicillium italicum growth. J. Food Prot. 62, 540–542. doi: 10.4315/0362-028X-62.5.541
McGinty, D., Letizia, C. S., and Api, A. M. (2012a). Fragrance material review on phenethyl isobutyrate. Food Chem. Toxicol. 50, S416–S419. doi: 10.1016/j.fct.2012.02.065
McGinty, D., Vitale, D., Letizia, C. S., and Api, A. M. (2012b). Fragrance material review on phenethyl acetate. Food Chem. Toxicol. 50, S491–S497. doi: 10.1016/j.fct.2012.02.089
Nortjé, K., Buys, E. M., and Minnaar, A. (2006). Use of γ-irradiation to reduce high levels of Staphylococcus aureus on casein–whey protein coated moist beef biltong. Food Microbiol. 23, 729–737. doi: 10.1016/j.fm.2006.02.001
Obenland, D., Collin, S., Mackey, B., Sievert, J., and Arpaia, M. L. (2011). Storage temperature and time influences sensory quality of mandarins by altering soluble solids, acidity and aroma volatile composition. Postharvest Biol. Technol. 59, 187–193. doi: 10.1016/j.postharvbio.2010.09.011
Olofsson, T. C., Ahrné, S., and Molin, G. (2007). The bacterial flora of vacuum-packed cold-smoked salmon stored at 7°C, identified by direct 16S rRNA gene analysis and pure culture technique. J. Appl. Microbiol. 103, 109–119. doi: 10.1111/j.1365-2672.2006.03216.x
Qiu, M., Xiao, X., Xiao, Y., Ma, J., Yang, H., Jiang, H., et al. (2022). Dynamic changes of bacterial communities and microbial association networks in ready-to-eat chicken meat during storage. Food Secur. 11:3733. doi: 10.3390/foods11223733
Rönnander, J., Ljunggren, J., Hedenström, E., and Wright, S. A. I. (2018). Biotransformation of vanillin into vanillyl alcohol by a novel strain of Cystobasidium laryngis isolated from decaying wood. AMB Express 8:137. doi: 10.1186/s13568-018-0666-4
Samapundo, S., de Baenst, I., Eeckhout, M., and Devlieghere, F. (2017). Inhibitory activity of fermentates towards Zygosaccharomyces bailii and their potential to replace potassium sorbate in dressings. LWT Food Sci. Technol. 79, 309–315. doi: 10.1016/j.lwt.2017.01.055
Shukla, R., Chauhan, N., Rajak, C., and Flora, S. (2019). “Flavors and fragrances” in Flavor development for functional foods and nutraceuticals, vol. 141. eds. M. Selvamuthukumaran and Y. Pathak (Boca Raton, FL: CRC Press LLC).
Wu, T., and Cadwallader, K. R. (2019). Identification of characterizing aroma components of roasted chicory "coffee" brews. J. Agric. Food Chem. 67, 13848–13859. doi: 10.1021/acs.jafc.9b00776
Xia, J., Jiang, N., Zhang, B., Sun, R., Zhu, Y., Xu, W., et al. (2023). Bacterial changes in boiled crayfish between different storage periods and characterizations of the specific spoilage Bacteria. Food Secur. 12:3006. doi: 10.3390/foods12163006
Zhang, Q. Q., Li, D., Zhang, W., Jiang, M., Chen, X. H., and Dong, M. S. (2021). Comparative analysis of the bacterial diversity of Chinese fermented sausages using high-throughput sequencing. LWT 150:111975. doi: 10.1016/j.lwt.2021.111975
Zhang, L., Wang, Y., Ma, Y., Deng, J., Ni, Y., and Xu, B. (2023). Monitoring and traceability of spoilage bacteria from three representative braised chickens based on culture-dependent and high-throughput sequencing. LWT 182:114824. doi: 10.1016/j.lwt.2023.114824
Zhang, W., Yu, Y., He, H., Lv, X., Liu, Z., and Ni, L. (2022). The adhesion and spoilage of Shewanella putrefaciens in Tilapia. Food Secur. 11:1913. doi: 10.3390/foods11131913
Zheng, R., Xu, X., Xing, J., Cheng, H., Zhang, S., Shen, J., et al. (2020). Quality evaluation and characterization of specific spoilage organisms of Spanish mackerel by high-throughput sequencing during 0 °C cold chain logistics. Food Secur. 9:312. doi: 10.3390/foods9030312
Zhu, Y., Guo, Y., Yang, F., Zhou, C., Tang, C., and Zhou, G. (2021). Combined application of high-throughput sequencing and UHPLC-Q/TOF-MS-based metabolomics in the evaluation of microorganisms and metabolites of dry-cured ham of different origins. Int. J. Food Microbiol. 359:109422. doi: 10.1016/j.ijfoodmicro.2021.109422
Zhu, G., and Xiao, Z. (2022). “Flavors and fragrances: structure of various flavors with food ingredients” in Flavors and fragrances in food processing: Preparation and characterization methods, vol. 1433. eds. P. Balakrishnan and S. Gopi (Washington, DC: American Chemical Society), 21–188.
Keywords: flavors and fragrances, whole genome shotgun (WGS), high-throughput sequencing, spoilage, GC–MS
Citation: Feng Y, Zhang T, Yang J, Liu W, Yang Y, Huang J, Huang S, Yang Z, Liu Q, Zheng W and Zhou Q (2025) Characterization of microbial communities in flavors and fragrances during storage. Front. Microbiol. 16:1516594. doi: 10.3389/fmicb.2025.1516594
Edited by:
Giovanna Suzzi, University of Teramo, ItalyReviewed by:
Bowen Li, Southwest University, ChinaDilhun Keriman Arserim Uçar, Bingöl Univesity, Türkiye
Copyright © 2025 Feng, Zhang, Yang, Liu, Yang, Huang, Huang, Yang, Liu, Zheng and Zhou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qing Zhou, emhvdXFpbmcwMTI0QDE2My5jb20=