Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 12 August 2025

Sec. Industrial Biotechnology

Volume 13 - 2025 | https://doi.org/10.3389/fbioe.2025.1636506

Isolation and evaluation of a Bacillus altitudinis strain to improve cigar tobacco leaves fermentation

Hejun LiuHejun Liu1Xin Fang,Xin Fang1,2Shiping GuoShiping Guo3Bo ZhangBo Zhang1Chongde Wu,Chongde Wu1,2Qiu ZhongQiu Zhong4Ruina ZhangRuina Zhang4Hongzhi ShiHongzhi Shi5Yanqing Qin
Yanqing Qin3*Shuhua Zeng
Shuhua Zeng6*Yao Jin,
Yao Jin1,2*
  • 1College of Biomass Science and Engineering, Sichuan University, Chengdu, China
  • 2Key Laboratory of Leather Chemistry and Engineering, Ministry of Education, Sichuan University, Chengdu, China
  • 3R&D Department, Sichuan Provincial Branch of China National Tobacco Crop Tobacco Science Institute, Chengdu, China
  • 4Deyang Tobacco Company of Sichuan Province, Deyang, Sichuan, China
  • 5College of Tobacco Science, Henan Agricultural University, Zhengzhou, China
  • 6College of Agronomy, Sichuan Agricultural University, Chengdu, China

Introduction: High levels of nitrogenous compounds such as proteins and alkaloids in cigar tobacco leaves can negatively impact combustion, smoke smoothness, and overall sensory quality. At the same time, the formation of aroma-active compounds during fermentation is essential for desirable flavor development. This study aimed to isolate a functional bacterial strain capable of accelerating nitrogen compound degradation and enhancing aroma quality through bioaugmented fermentation.

Methods: A total of 65 bacterial strains were isolated from naturally fermented cigar tobacco using a tobacco extract-based medium. Sixteen showed significant protease activity, and Bacillus altitudinis (CCTCC M20211370) was selected for further study. Laboratory- and industrial-scale fermentation trials were conducted with this strain. Volatile compounds were analyzed using gas chromatography–mass spectrometry (GC-MS), microbial community dynamics were assessed via high-throughput sequencing, and sensory evaluations were performed to assess quality improvements.

Results: Inoculation with Bacillus altitudinis significantly increased desirable aroma compounds such as neophytadiene and β-ionone, while reducing harmful compounds including nicotine and myosmine. The microbial structure was reshaped, with enrichment of beneficial genera like Bacillus and Oceanobacillus. Sensory analysis confirmed improved aroma complexity, featuring enhanced floral, honey-sweet, and resinous notes.

Discussion: Bacillus altitudinis effectively improved cigar tobacco quality by promoting nitrogen compound degradation and modulating microbial metabolism to enhance aroma development. These findings support its potential as a bioaugmentation agent in industrial cigar fermentation. Further genomic and enzymatic studies are warranted to elucidate its functional mechanisms and facilitate large-scale application.

1 Introduction

Tobacco, one of the most widely cultivated and consumed non-food crops globally, is commercially available in diverse forms (Liu T. et al., 2022). Among them, cigars are particularly prized for their distinctive flavor. As a significant commercial tobacco product, cigars enjoy widespread popularity in North America and internationally (Vilela et al., 2020), undergoing drying, fermentation, and aging processes (Jensen and Parmele, 1950).

Cigar tobacco’s upper leaves contain macromolecules such as starch and proteins, with a low sugar-to-nitrogen ratio. High starch content negatively affects the combustion properties of the tobacco, while excessive protein leads to unpleasant odors and irritants during combustion. High pectin content can result in undesirable woody and off-flavors in the smoke (Liu F. et al., 2022). Additionally, the total alkaloid content, total sugars, potassium levels, and the alkali-to-nitrogen ratio all influence the aroma profile, aftertaste, and irritancy of the tobacco. Water-soluble sugars play a crucial role in neutralizing acidic and alkaline smoke compounds during combustion, reducing smoke irritation (Liu F. et al., 2022; Wu Q. et al., 2023) and acting as important precursors for tobacco aroma (Chen et al., 2019). Therefore, fermentation is a key treatment for enhancing both combustion properties and flavor characteristics of the tobacco.

Cigar tobacco fermentation is critical to its final quality and flavor. Through complex microbial metabolic activities, macromolecules like proteins and starch are broken down into smaller molecules, leading to the formation of aromatic compounds (Shibin et al., 2021). However, traditional fermentation processes suffer from long cycles, low efficiency, and inconsistent flavor stability, limiting the large-scale production of high-quality cigars. Recent advancements have focused on the use of functional microorganisms to enhance fermentation and regulate tobacco metabolism, offering promising solutions for improving cigar quality. By selectively screening microorganisms with specific metabolic functions, fermentation cycles can be shortened, and flavor compound production enhanced (Zhang et al., 2024). Therefore, exploring the roles of functional microorganisms in cigar tobacco fermentation and analyzing their metabolic mechanisms hold great theoretical and practical value (Wen et al., 2024).

Previous studies have shown that excessive nitrogen-containing compounds are a major cause of off-flavors in cigar smoke, making their reduction crucial for improving cigar quality (Chunpin, 2015). The microbial community on cigar tobacco leaves (CTLs) is diverse, with genera such as Bacillus, Actinobacteria, and Pseudomonas playing key roles in fermentation (Liu T. et al., 2022; Ning et al., 2023). However, most existing studies focus on microbial diversity, and systematic research on the targeted screening of functional microorganisms and their metabolic regulation mechanisms remains limited (Zhang et al., 2023). The relationships between microorganisms and volatile aroma compounds are not well understood, which hinders the precise application of functional microorganisms in fermentation processes (Liu et al., 2021; Wu X. et al., 2023).

Therefore, this study isolated several core functional microorganisms with protein-degrading abilities from high-quality CTLs. A strain of Bacillus altitudinis (CCTCC NO: M 20211370) was chosen to serve as the bioaugmentation agent for the fermentation process of cigar tobacco. The effect of such strain on volatile compounds, microbial structural of tobacco fermentation were then evaluated.

This study systematically reveals the functional characteristics of Bacillus altitudinis in cigar tobacco fermentation, providing a theoretical and technical foundation for microbial community regulation and fermentation optimization. The findings not only expand the application of functional microorganisms in the tobacco industry but also propose new research directions for developing bioaugmentation strategies for high-quality fermentation industry.

2 Materials and methods

2.1 Tobacco samples and culture conditions

The cigar tobacco samples, consisting of naturally fermented high-grade second-class tobacco leaves including wrapper, binder, and filler types, were from Shifang, provided by the Sichuan Provincial Branch of China National Tobacco Crop Tobacco Science Institute.

Tobacco extract was prepared using fermented cigar tobacco from Shifang as raw material. The tobacco leaves were finely ground and then subjected to a 40-mesh sieve to ensure a uniform particle size. In the preparation process, 15 g of the resulting tobacco powder was combined with 400 mL of distilled water. This mixture was then heated to a boil and allowed to extract for a duration of 30 min. After boiling, the extract was filtered through two layers of gauze to remove any solid residues. To preserve the integrity of the extract for subsequent applications, it was cooled and stored at a temperature of 4 °C.

Protein-hydrolyzing microorganisms were cultured and isolated, and protease activity was assayed under specified conditions: tobacco extract (TE) medium (100 mL L-1 tobacco leaf aqueous extract, 900 mL L-1 distilled water, 4.5 g L-1 yeast extract, 9.0 g L-1 tryptone, 9.0 g L-1 NaCl, with natural pH), fermentation medium (5.0 g L-1 cigar tobacco leaf powder and 20.0 g L-1 glucose), casein medium (10.0 g L-1 casein, 4.0 g L-1 NaCl, 5.0 g L-1 yeast extract and 20.0 g L-1 agar powder, adjusted to pH 7.2). All the chemicals and reagents used were of analytical grade.

2.2 Isolation and purification of microorganisms from tobacco leaf surfaces

Exactly 1 g of high-quality fermented tobacco leaf sample was accurately weighed and finely chopped using sterile scissors under aseptic conditions. The shredded tobacco was immersed in 100 mL of sterile water and allowed to incubate at 37 °C while shaking at 200 rpm for 30 min. Following the incubation period, the mixture underwent filtration and was then centrifuged at 9,000 × g for 15 min. The resulting supernatant was removed, and the pellet was reconstituted in 5 mL of sterile water to obtain the initial bacterial suspension.

The original bacterial suspension was serially diluted (10–1 to 10–5), and 100 μL of every dilution was uniformly spread on TE agar. After incubating the plates upside down at 37 °C for 48 h, bacterial colonies became visible. Single colonies were purified by quadrant streaking until pure isolates were obtained.

2.3 Directed screening of functional microorganisms

Single colonies were inoculated into casein agar and incubated at 37 °C in an inverted position for 48 h. Strains with faster growth and larger casein hydrolysis zones were selected for subsequent purification. The selected isolates were pretreated using the glutaraldehyde fixation method (Feng et al., 2019) for scanning electron microscopy observation with a SU8010 microscope (Hitachi High-Technologies, Japan) at a magnification of ×10,000 , using secondary electron signals with a working distance (WD) of 9.9 mm.

Genomic DNA extraction was carried out utilizing a commercially available DNA extraction kit. The identification of bacteria was achieved through PCR amplification targeting the 16S rDNA region employing universal primers 27F (5′-AGTTTGATCMTGGCTCAG-3′) and 1492R (5′-GGT​TAC​CTT​GTT​ACG​ACT-3′), with products generated using the 1×TSE101 Gold Mix from Qingke Biotechnology Co., Ltd. (China), resulting in amplicons of about 1,500 bp. The assembled amplified sequences were processed through ContigExpress and subsequently uploaded to the NCBI database. For sequence alignment, the Basic Local Alignment Search Tool was utilized, and a phylogenetic tree was created employing the neighbor-joining method in accordance with the Kimura two-parameter model, including 1,000 bootstrap replicates (Shin et al., 2014; Sakulworakan et al., 2021). The reagents and equipment used for sequencing were supplied by Shanghai Meiji Biomedical Technology Co., Ltd. (Shanghai, China).

2.4 Optimization of enhanced fermentation conditions

Functional microorganisms isolated and purified in this study were streaked onto solid media to obtain single colonies, which were subsequently inoculated into liquid fermentation medium and incubated at 37 °C for 48 h to prepare the initial bacterial suspension. Bacterial cells were harvested by centrifuged at 8,000 × g for 10 min, resuspended in sterile physiological saline, and adjusted to a final concentration of approximately 1.0 × 109 CFU mL-1 to prepare the stock inoculum. An equal volume of cryoprotectant solution containing 10% (v/v) glycerol was added, and the mixture was thoroughly homogenized, aliquoted into sterile cryovials, and stored at −80 °C. Prior to use in fermentation experiments, the inoculum was revived from frozen stocks and cultured in liquid medium until reaching a concentration of 3.82 × 108 CFU mL-1 for application. The fermentation process was conducted at 37 °C with a relative humidity of 70% over a period of 10 days, using a seed liquid application rate of 100 mL per kilogram of tobacco material. Three experimental groups were established to evaluate fermentation effects under different treatments: an enhanced fermentation group (GD), in which seed liquid was added; a blank control group (KB), which received the same volume of distilled water instead of seed liquid; and a natural fermentation group (ZR), which was left untreated without any additives or interventions.

2.5 Microbial community structure

Microbial community analysis of CTLs from different fermentation treatments was conducted using triplicate biological replicates. Genomic DNA was extracted from homogenized cigar tobacco samples (40-mesh particle size). DNA purity was assessed using a NanoDrop 2000 spectrophotometer (Shoeb et al., 2017), and DNA integrity was verified via 1% agarose gel electrophoresis.

To analyze bacterial community composition, the V3–V4 region of the 16S rRNA gene was amplified using universal primers 338F (5′-ACT​CCT​ACG​GGA​GGC​AGC​AG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) (Penland et al., 2020). For fungal diversity analysis, the ITS1 region was amplified using primers ITS1F (5′-CTT​GGT​CAT​TTA​GAG​GAA​GTA​A-3′) and ITS2R (5′-GCT​GCG​TTC​TTC​ATC​GAT​GC-3′) (Rimington et al., 2018). The PCR products underwent gel purification and were quantified with a Quantus™ fluorometer (Promega, United States).

Purified amplicons were combined in equal molar concentrations and then processed through paired-end sequencing (2 × 300 bp) on the Illumina MiSeq PE300 platform (Illumina, San Diego, CA, United States). The initial sequencing reads underwent quality filtering with fastp (version 0.20.0) (Chen et al., 2018) and were subsequently merged using FLASH (version 1.2.7) (Magoč and Salzberg, 2011). OTUs were established at 97% sequence similarity (Edgar, 2013) using UPARSE (version 7.1) (Stackebrandt and Goebel, 1994), while chimeric sequences were eliminated. The taxonomic classification was carried out with the RDP Classifier (version 2.2) against the SILVA database for bacteria (16S rRNA) and the UNITE database for fungi (ITS), utilizing a confidence threshold of 70%. The raw sequencing data has been submitted to the NCBI SRA database, with the BioProject accession number PRJNA824254.

2.6 Determination of volatile flavor compounds

The analysis of volatile metabolites from fermented CTLs was conducted after a 7-day fermentation period utilizing a GC-MS system (7890B/5977A, Agilent Technologies, United States) outfitted with a HS-SPME apparatus. An exact weight of 0.5 g of finely ground tobacco was combined with 10 μL of ethyl phenylacetate (chromatographic grade) serving as the internal standard. The samples were subjected to equilibrium at 70°C for 20 min, followed by extraction for 30 min. Finally, qualitative and quantitative analysis of the aroma components was performed using the NIST mass spectral library and Agilent ChemStation.

2.7 Sensory quality evaluation

To further validate the fermentation capability of Bacillus altitudinis, we conducted an industrial-scale fermentation verification experiment with a 1500 kg batch. The CTLs used for sensory evaluation were sourced from this industrial verification experiment. The fermented CTLs was rolled into cigars measuring 120 mm in length and 15 mm in diameter. The cigars were conditioned at 20°C and 60% relative humidity to achieve moisture equilibrium. Subsequently, sensory quality evaluation was conducted by seven trained panelists specializing in cigar production and assessment, following the standardized evaluation form provided by the Changcheng Cigar Factory (Jia et al., 2023). Sensory analysis was performed based on two main aspects: quality characteristics and flavor attributes. For quality characteristics, eleven parameters (e.g., richness, maturity, and irritation) were rated on a 0–9 scale, with higher scores indicating better performance. For flavor attributes (e.g., floral, roasted, and honey-sweet notes) were rated on a 1–5 scale, with higher scores denoted greater intensity. Consensus on evaluation scores for each sample was achieved among all panel members.

2.8 Analytical methods

Proteolytic activity was assessed by casein hydrolysis zone analysis (Shen et al., 2021). The diameter of the bacterial colony (d) and the diameter of the hydrolysis (clear) zone (D) were measured using a ruler. The proteolytic index (I) was calculated to evaluate the strain’s protein degradation ability Equation 1.

I=D/d(1)

Strains with strong protein hydrolysis ability were inoculated into tobacco extract liquid medium and cultured at 37 °C until the optical density at 660 nm (OD660) reached approximately 1.0.

Protease Activity Assay: A crude enzyme solution of 5 mL was subjected to centrifugation at 8,000 × g for 5 min at a temperature of 4 °C. The resulting precipitate was removed, and the supernatant was combined with 50 mL of phosphate buffer, followed by appropriate dilution to create the enzyme test solution. The measurement of protease activity was conducted using UV spectrophotometry at a wavelength of 660 nm with a spectrophotometer (Model L7, INESA Scientific Instrument Co., Ltd., Shanghai, China).

Total nitrogen, protein, starch, reducing sugars, and total sugars were quantified using a continuous flow analysis system (San++Classic, Skalar, Netherlands). The methods for total nitrogen determination followed GB 5009.33–2016, protein content was assessed as per YC/T 166–2003, starch levels were measured according to GB 5009.9–2016, and reducing sugars were evaluated using the GB/T 5513–2019 standard. For each sample, three biological replicates were performed.

2.9 Statistical methods

Principal coordinate analysis (PCoA) was utilized to assess variations in microbial communities among cigar samples originating from various fermentation groups, employing the Bray–Curtis distance matrix for both dimensionality reduction and visualization. The relative quantities of metabolites were determined through a semi-quantitative peak area normalization technique. Compounds that exhibited similarity index (SI) and reverse similarity index (RSI) scores exceeding 750 were chosen for further investigation. The relationships between volatile metabolites and microbial communities post Bacillus altitudinis-enhanced fermentation were illustrated with heatmaps generated using R software (version 4.1.2) (Hoyles et al., 2018). Functional annotation along with pathway analysis for genes and enzymes was executed based on the KEGG database. In examining the correlations between volatile metabolites and microbial communities, a significance level of p < 0.05 (n = 3) was established for evaluating statistical significance.

3 Results

3.1 Strain isolation and identification

Previous studies have shown that the protein content of the macromolecular substances in naturally fermented CTLs is associated with the flavor quality of the tobacco (Wen et al., 2023). A high protein content can lead to poor smoking characteristics, producing unpleasant odors and off-flavors (Song et al., 2024). Therefore, elimination of such off-flavors in CTLs requires a reduction in the protein content, which can improve the smoking quality of cigars. In this study, we focused on selecting microorganisms with the ability to degrade the macromolecular protein substances in tobacco leaves, using protein degradation as a screening principle. The microorganisms were then preserved, identified, and subjected to further fermentation performance testing.

A total of 65 bacterial isolates were obtained through gradient dilution plate method. Given the key regulatory role of protein metabolism in the fermentation process, protein hydrolysis ability was established as the core screening criterion. 16 strains producing protease were selected based on casein hydrolysis circle analysis. Among them, the strain (CCTCC NO: M 20211370) exhibited outstanding degradation efficiency, with a hydrolysis index (D/d value) of 3.955 (Figure 1. A). Further quantitative enzyme activity assays confirmed that its protease activity peak reached 577.15 U mL-1 (Supplementary Figure S1), establishing it as the preferred strain for subsequent fermentation enhancement studies. The colonies of the strain appeared round, opaque, and slightly raised, with rough or irregular edges and a dry to matte surface texture. Gram staining indicated that the strain is Gram-positive (Figure 1. B). Scanning electron microscopy revealed that the cells exhibited a rod-shaped morphology (1.2–2.5 μm) with a rough surface and aggregated arrangement (Figure 1. C). Phylogenetic analysis showed that the strain is a Gram-positive Bacillus, displaying 100% sequence identity with Bacillus altitudinis (Figure 1. D). The strain has been preserved in a −80 °C glycerol stock for further study.

Figure 1
Panel A shows a petri dish with a microbial growth inhibition zone marked D/d=3.955. Panel B features several petri dishes with bacterial colonies at different growth stages. Panel C displays a scanning electron microscope image of rod-shaped bacteria. Panel D presents a phylogenetic tree showing various Bacillus species, highlighting their genetic relationships with branch support values.

Figure 1. Screening and identification of functional microorganisms. (A) The experimental results of casein hydrolysis by Bacillus altitudinis strain; (B) Morphological observation and Gram staining of strains; (C) Scanning electron microscopy results of the strain; (D) Phylogenetic tree of strain Bacillus altitudinis.

3.2 Comparative analysis of volatile compounds

A total of 21 volatile metabolites were identified from cigar tobacco samples, including three tobacco alkaloids, five terpenoids, five esters, six aldehydes and ketones, and two acids (Supplementary Table S1). Compared to the other two groups, the GD group exhibited a significant alteration in the composition of volatile flavor compounds: the content of tobacco alkaloids decreased by 31.6%, aldehydes and ketones increased by 13.4%, terpenoids increased by 64.3%, acids Dramatically increased and small increases were also observed in esters (Figure 2).

Figure 2
Five bar charts labeled A to E compare the concentration of different compounds in micrograms per kilogram across three groups: GD, KB, and ZR. Chart A shows tobacco alkaloids with increasing values from GD to ZR. Chart B displays aldoketones, with GD having the highest concentration. Chart C presents acids, where GD is significantly higher than KB. Chart D indicates that GD has a higher terpene concentration than KB. Chart E shows similar ester levels between GD and KB, slightly reducing in ZR. Asterisks indicate significant differences between groups.

Figure 2. Total volatile substances in cigar tobacco of Shifang binder leaf by Bacillus altitudinis enhanced fermentation. (A) Tobacco alkaloids; (B) Aldoketones; (C) Acids; (D) Terpene; (E) Esters.

In terms of aroma-contributing substances, the GD group presented 10 additional compounds compared to the ZR and KB groups, including plant-derived furans, ethyl palmitate, linalool, geranylgeraniol, (±)-3-hydroxylauric acid, and citrus ketone. Moreover, levels of solanone, geranylacetone, neophytadiene, β-selinene, and cedrene were elevated, while geranylacetone slightly decreased. Furans produced by Maillard reactions are known to contribute caramel-like aromas (Seeman et al., 2003; Song et al., 2024). The increase in aldehydes and ketones is favorable for flavor enhancement due to the presence of carbonyl groups, which are key aroma-active moieties and impart pleasant fragrances (Lei et al., 2015). Solanone contributes to the smoothness of tobacco aroma (Hong-qi, 2005; Teng, 2013), ionone adds sweetness (Wang et al., 2022), and decanal provides a fresh, citrus-like scent (Zhang et al., 2022). Terpenoids are important aromatic components in tobacco, such as neophytadiene and cedrene (Wang et al., 2024).

These results suggest that Bacillus altitudinis not only degraded nitrogen-containing compounds such as tobacco alkaloids, but also activated additional metabolic pathways, thereby increasing both the quantity and diversity of terpenoids and carbonyl compounds. This activity is likely to contribute positively to quality improvement and aroma enhancement during cigar fermentation.

3.3 Comparative analysis on microbial structure

Microbial community diversity was analyzed across the three groups of fermented cigar tobacco leaf samples. A total of 430,590 optimized bacterial sequences and 179,283,491 fungal sequences were obtained, with an average sequence length of 416 bp. Taxonomic classification revealed: 1 domain, 1 kingdom, 11 phyla, 18 classes, 42 orders, 63 families, 95 genera, 121 species, and 142 OTUs. The GD group exhibited a decreasing trend in bacterial diversity and richness, while fungal diversity showed no significant difference (Supplementary Figure S2).

At the phylum level, the top four bacterial phyla were Firmicutes, Actinobacteriota, Proteobacteria, and Chloroflexi. At the genus level, the top four were Bacillus, Rhodococcus, norank_f__Pseudonocardiaceae, and Staphylococcus. To further investigate the impact of the added functional microorganism Bacillus altitudinis on microbial community structure, community composition was analyzed at both the phylum and genus levels (Figure 3).

Figure 3
Bar charts labeled A to D show microbial community compositions across three groups, GD, KB, and ZR. Chart A shows bacterial composition by phylum, with Firmicutes, Actinobacteriota, Proteobacteria, and others. Chart B details bacterial genus composition, including Bacillus and Pseudomonas. Chart C presents fungal composition by phylum, featuring Ascomycota and Basidiomycota. Chart D shows fungal genus composition, including Aspergillus and Penicillium. Each chart uses different colors to represent specific taxa proportions.

Figure 3. Community structure of bacteria and fungi in fermented cigar leaves. (A) Bacteria phylum level; (B) bacteria genus level; (C) fungal phylum level; (D) fungal genus level.

At the phylum level, Firmicutes and Ascomycota were dominant in the GD group, whereas Actinobacteria and Ascomycota dominated in the KB and ZR groups. Compared to the other two groups, the GD group exhibited an increase in Firmicutes and a decrease in Actinobacteria among bacteria, while in fungi, Ascomycota decreased and Basidiomycota showed an increasing trend. At the genus level, the GD group showed elevated abundances of Oceanobacillus and Bacillus, while Rhodococcus, Staphylococcus, and Achromobacter decreased. For fungi, Sampaiozyma, Cladosporium, and Wallemia increased, whereas Aspergillus and Sagenomella declined. Notably, Bacillus became the dominant bacterial genus in the GD group, significantly altering the bacterial community composition of the fermented tobacco leaves (Zhou et al., 2024).

PCoA was conducted to compare the microbial community composition of CTLs across different fermentation groups (Figure 4). For bacterial communities, the variances explained by PC1 and PC2 were 82.87% and 11.18%, respectively, while for fungal communities, PC1 and PC2 accounted for 59.11% and 13.66% of the variance. The results indicated that both bacterial and fungal communities in the GD group were significantly distinct from those in the KB and ZR groups. In contrast, the ZR and KB groups clustered closely together, reflecting similar microbial community structures. These dimensionality reduction results demonstrate that the addition of the functional strain Bacillus altitudinis exerted a significant impact on the composition of microbial communities during cigar tobacco fermentation. Both bacterial and fungal communities in the enhanced fermentation group (GD) exhibited marked differences compared to the other groups.

Figure 4
Panel A shows a PCoA plot of bacterial OTU levels with three groups: GD (red circles), KB (blue triangles), and ZR (green diamonds). Panel B displays a PCoA plot of fungal OTU levels with groups GD (red circles), KB (pink triangles), and ZR (brown diamonds). Each plot includes axis labels indicating the percentage of variation explained.

Figure 4. Principal coordinate analysis of bacteria (A) and fungi (B) in enhanced fermentation cigar leaves.

3.4 Correlation between microoganisms and volatile compounds

As the functional microorganism introduced in this study is a bacterium, the correlation network between bacterial communities and volatile metabolites was specifically analyzed, aiming to facilitate the future identification of beneficial microbes in cigar fermentation. A total of 20 representative volatile metabolites and the top 20 most abundant bacterial genera in the fermentation system were selected for correlation analysis (Figure 5).

Figure 5
Heatmap showing the correlation between microbial genera and chemical compounds. Rows represent bacterial genera, while columns represent various compounds. The scale ranges from -1 (blue) to 1 (red), indicating negative to positive correlations, respectively. Notable clusters of red indicate strong positive correlations, while blue indicates negative correlations.

Figure 5. Heatmap of correlation between volatile metabolites and main microorganisms in CTLs with enhanced fermentation (|ρ | > 0.5, p-value <0.05, *p < 0.05, **p < 0.01, ***p < 0.001).

First, nicotine degradation products such as myosmine and nicotyrine were significantly correlated with Bacillus and Oceanobacillus (Song et al., 2024). Geranylacetone and solanone were strongly associated with Oceanobacillus, Staphylococcus, Ochrobactrum, Delftia, Rhodococcus, and Achromobacter. Metabolites such as (±)-3-hydroxydodecanoic acid, β-ionone, nootkatone, methyl 2,5-octadecadiynoate, ethyl palmitate, (+)-α-cedrene, and 2Z,6E-farnesol exhibited significant correlations with Sphingomonas and Methylobacterium-Methylorubrum. Notably, neophytadiene and cedrene showed significant associations with 13 bacterial genera. Excess alkaloids in cigars not only impart bitterness and harshness during smoking but also pose health risks (Jia et al., 2023). It is widely accepted that alkaloid degradation reduces leaf harshness, resulting in a smoother and milder flavor profile (Xue et al., 2023), thus improving cigar tobacco quality. Correlation analysis confirmed the potential of microbial nicotine degradation. For instance, the nicotine degradation products myosmine and nicotyrine were significantly associated with Bacillus and Oceanobacillus—with the enhanced fermentation system specifically supplemented with Bacillus altitudinis, a member of the Bacillus genus. Additionally, solanone-related compounds are reported to play a key role in reducing impurities and harshness while enhancing the richness, smoothness, aftertaste, sweetness, and cleanliness of cigar aroma (Popova et al., 2019; Zheng et al., 2022). These compounds showed strong correlations with Staphylococcus, Delftia, and Rhodococcus. Terpenoids, which are derived from the MVA pathway and structurally composed of isoprene (C5) units (Tian et al., 2018), are known to be involved in the biosynthesis of plant volatiles (Cai et al., 2014; Wei et al., 2018; Shi et al., 2019; Czajka et al., 2020; Tang et al., 2022) and contribute unique aromas to cigars (Laureati et al., 2020; Song et al., 2020; Chen et al., 2022). These include compounds such as cedrene, neophytadiene, (+)-α-cedrene, and 2Z,6E-farnesol. Significant correlations were observed between terpenoids and genera including Pantoea, Pseudomonas, and Sphingomonas. Of particular note, Sphingomonas was significantly associated with multiple terpenoids and ketone-aldehyde compounds, including β-ionone, cedrene, and neophytadiene.

These findings suggest that further optimization of the bacterial community structure in cigar tobacco could enhance the content of aroma-contributing compounds. This highlights a promising direction for the discovery of beneficial microbes in cigar fermentation.

3.5 Sensory evaluation

Sensory evaluation is generally considered the most effective approach for assessing the fermentation performance of CTLs and the overall quality of cigars. To verify whether the isolated Bacillus altitudinis strain could improve CTLs quality through fermentation, sensory analysis was conducted on CTLs subjected to Bacillus altitudinis-enhanced fermentation, with water-fermented CTLs serving as the control. As shown in Figure 6, the total sensory scores of CTLs in the enhanced fermentation group were higher than those of the control group. Specifically, Figure 6A indicates that Bacillus altitudinis improved the cleanliness and maturity of CTLs while reducing irritation and off-flavors. In terms of aroma intensity, the fermentation process led to a richer, more mellow, and more harmonious aroma profile. Regarding aroma characteristics (Figure 6B), Bacillus altitudinis enhanced floral, honey-sweet, hay-like, and resinous notes in the CTLs. These sensory evaluation results were consistent with the volatile flavor compound analysis, confirming that Bacillus altitudinis-enhanced fermentation can effectively improve the quality of CTLs.

Figure 6
Radar charts labeled A and B compare attributes of GD (orange) and KB (yellow). Chart A evaluates intensity, richness, matureness, irritation, smoothness, sweetness, cleanliness, aroma volume, combustibility, ash color, and ash coagulation. Chart B assesses hay, honey, floral, resinous, baking, and woody characteristics.

Figure 6. Radar plots of sensory evaluation of CTLs subjected to water fermentation and Bacillus altitudinis-enhanced fermentation. (A) Quality characteristics; (B) Flavor attributes. GD represents the enhanced fermentation group, and KB represents the control group.

4 Discussion

Through functional strain-enhanced fermentation, the flavor and smoking characteristics of CTLs can be significantly improved. In order to enhance cigar quality, reduce the nicotine content of CTLs, and optimize cigar flavor, this study isolated a strain of Bacillus altitudinis from high-quality CTLs produced in the Shifang region. This strain was used to enhance the fermentation of CTLs. GC-MS analysis of the flavor compounds after enhanced fermentation revealed a 31.6% reduction in tobacco alkaloid content, a 13.4% increase in aldehyde and ketone levels, a 64.3% increase in terpenoid compounds, significant increases in acid content, and a slight increase in ester compounds. Ten new compounds with desirable flavor characteristics were identified, including linalool, geraniol, (±)-3-hydroxylauric acid, and citrus ketone. The aroma profile analysis showed enhanced floral, honey-like, hay, and resinous notes in CTLs. Regarding smoking characteristics, the effective degradation of nitrogen-containing compounds improved smoothness, maturity, and reduced irritants and off-flavors, making the aroma softer and more refined. Sensory evaluation results aligned with the analysis of volatile flavor compounds, confirming that the Bacillus altitudinis strain could improve the quality of CTLs through fermentation.

The fermentation of CTLs is a synergistic outcome of microbial community interactions. Nicotine degradation is closely associated with the genera Bacillus and Oceanobacillus (Song et al., 2024). The formation of geranylacetone and solanone is significantly correlated with Staphylococcus, Ochrobactrum, and Achromobacter, while (±)-3-hydroxydodecanoic acid, β-ionone, and 2Z,6E-farnesol are strongly associated with Sphingomonas and Methylobacterium-Methylorubrum. In the present study, the addition of Bacillus altitudinis significantly altered the microbial community structure of CTLs, as confirmed by PCoA. Specifically, bacterial community changes were characterized by increased abundances of Oceanobacillus and Bacillus, and decreased abundances of Rhodococcus, Staphylococcus, and Achromobacter. Within the fungal community, the relative abundances of Sampaiozyma, Cladosporium, and Wallemia increased, while those of Aspergillus and Sagenomella decreased; however, these changes in the fungal community were not statistically significant.

To further investigate the enzymatic-level effects of fermentation enhanced by Bacillus altitudinis, based on KEGG database annotations, the relationships among genomic, functional, and enzymatic features were analyzed, and enzyme abundance was quantified accordingly. Following the addition of Bacillus altitudinis, significant shifts in the enzymatic profile of the microbial community in CTLs were observed (Figure 7). Interestingly, although the fungal community showed no significant changes in species abundance or diversity following enhanced fermentation, functional annotation of fungal genes revealed marked differences in enzyme expression compared to the control group. This suggests that the addition of Bacillus altitudinis, while not significantly altering the composition of the fungal community, influenced its gene expression and metabolic pathways, thereby resulting in substantial shifts in the fungal enzymatic expression profile. Specifically, enzymes associated with protein degradation and amino acid metabolism showed markedly increased abundances, including acetyl-CoA C-acetyltransferase, acireductone dioxygenase (Ni2+-requiring), protein-glutamine glutaminase, glycine C-acetyltransferase, metalloendopeptidases, and aspartate kinase. Among them, acetyl-CoA C-acetyltransferase is involved in the degradation of valine, leucine, isoleucine, and lysine, as well as tryptophan metabolism; acireductone dioxygenase (Ni2+-requiring) participates in cysteine and methionine metabolism (Al-Mjeni et al., 2002); protein-glutamine glutaminase is related to glutamine hydrolysis (Kikuchi et al., 1971); glycine C-acetyltransferase and aspartate kinase are associated with glycine (Schmidt et al., 2001) and aspartate metabolism (Chassagnole et al., 2001; Curien et al., 2005), respectively. Metalloendopeptidases catalyze peptide bond hydrolysis in proteins (Bond and Beynon, 1995). These findings further support earlier experimental evidence that this functional strain is capable of degrading nitrogen-containing compounds such as proteins. Additionally, the enhanced fermentation process showed increased abundance of several enzymes involved in carbohydrate metabolism, such as cellulase, protein-Npi-phosphohistidine---sucrose phosphotransferase, levanase, and alpha-galactosidase, which are associated with the degradation of polysaccharides including glucans, fructans, and cellulose. A notable increase was also observed in malate dehydrogenase (oxaloacetate-decarboxylating), an enzyme involved in pyruvate metabolism (de Lorenzo et al., 2024).

Figure 7
Heatmap with hierarchically clustered data, showing enzyme activity levels from three samples: GD, KB, and ZR. Panel A displays enzymes involved in processes like sulfate adenylyltransferase and histidine kinase. Panel B focuses on enzymes such as adenylosuccinate synthase and glutathione peroxidase. Color gradient from blue to red indicates activity from low to high.

Figure 7. Systematic correlation analysis of genomic, functional, and enzymatic information was performed based on the KEGG database. (A) Variation in bacterial enzyme abundance across the GD, KB, and ZR groups; (B) Variation in fungal enzyme abundance across the GD, KB, and ZR groups (Supplementary Table S2 for detailed data).

These enzymatic changes may explain the observed enhancement in the production of volatile flavor compounds during fermentation. The addition of Bacillus altitudinis effectively accelerated the degradation of nitrogenous and carbohydrate substances, resulting in increased availability of small molecules for microbial metabolism. This in turn promoted active microbial metabolism and increased activity of terpene synthases (TPS), which are essential for the biosynthesis of volatile flavor compounds, most of which are derived from carotenoid metabolism (Liu et al., 2024). As illustrated, glucose and fructose are metabolized via glycolysis to yield acetyl-CoA and pyruvate, the starting substrates for the MVA and MEP pathways, respectively. These pathways lead to the formation of isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP), key precursors in carotenoid biosynthesis. Interestingly, the abundance of key enzymes in the MVA pathway decreased markedly after the addition of Bacillus altitudinis, while those in the MEP pathway remained unchanged, suggesting that carotenoid precursor biosynthesis primarily proceeds through the MEP pathway during enhanced fermentation (Figure 8).

Figure 8
Two horizontal bar charts labeled A and B display enzyme activity data. Chart A shows six enzyme codes with bars divided into sections representing ZR (red), GD (blue), and KB (green). Chart B shows six different enzyme codes with sections also representing ZR (blue), GD (red), and KB (green). Each bar represents percentages from zero to one hundred.

Figure 8. Based on KEGG pathway analysis, the variations in key enzymes related to terpenoid biosynthesis were compared among groups. (A) Changes in MVA pathway-related enzyme sys-tems in GD, KB, and ZR groups; (B) Changes in MEP pathway-related enzyme systems in GD, KB, and ZR groups (Supplementary Table S3 for detailed data).

To provide a more comprehensive analysis of the synthesis mechanisms of volatile flavor compounds, the metabolic pathway was visualized based on the KEGG database (Figure 9). The synthesis of carotenoid precursors involves multiple enzymes secreted by microorganisms. Geranylgeranyl diphosphate synthase (EC 2.5.1.1) is a critical enzyme that catalyzes the condensation of IPP and DMAPP to form geranylgeranyl diphosphate (GGPP), a central precursor in carotenoid biosynthesis. In addition to serving as substrates for carotenoids, IPP and DMAPP also lead to the production of mono- and sesquiterpenes such as cedrene (Wang et al., 2024). Carotenoid metabolism further yields various aromatic ketones, including β-ionone, solanone, and geranylacetone (Meléndez-Martínez et al., 2023), all of which contribute significantly to tobacco aroma. Based on the metabolite profiles and metabolic pathway analysis, it is postulated that the enhanced fermentation process activated additional terpene metabolic branches, leading to the biosynthesis of new terpenoid compounds such as (+)-α-cedrol and 2Z,6E-farnesol, which were not detected in the ZR and KB groups. Therefore, precise regulation of IPP and DMAPP metabolic flux may be a key strategy for improving aroma quality (Rasouli et al., 2018).

Figure 9
Diagram illustrating metabolic pathways involved in glucose metabolism, terpenoid backbone biosynthesis, and carotenoid biosynthesis. It shows connections between glycolysis, the TCA cycle, MVA and MEP pathways, and carotenoid synthesis. Key intermediates like Acetyl-CoA, Isopentenyl-PP, Geranyl-PP, and various enzymes (EC numbers provided) are highlighted. End products include alpha-Cedrene, Farnesol, and beta-Ionone.

Figure 9. Metabolic network of major volatile flavor compounds. The green area indicates the glycolysis pathway, the pink area represents the synthesis pathway of terpene skeleton com-pounds, and the blue area corresponds to the carotenoid metabolic pathway. Enzyme codes shown in blue boxes are derived from the KEGG database.

5 Conclusion

In conclusion, Bacillus altitudinis isolated from CTLs plays a key role in the fermentation process of cigar tobacco by regulating the tobacco microbiota, metabolic pathways, and improving aroma. Bacillus altitudinis can enhance the substrate availability on the surface of CTLs by degrading proteins and other macromolecules, activate carotenoid metabolic pathways, and promote the synthesis of flavor compounds. In the regulation of chemical composition, it can reduce the total amount of tobacco alkaloids and promote the synthesis of terpenes and other flavor compounds. In the regulation of the tobacco microbiota, it promotes the accumulation of aldehydes and ketones by selectively modulating key genera such as Oceanobacillus and Staphylococcus. Future research could further elucidate the coupling mechanism between the strain’s genomic functional modules and metabolic pathways, and explore its potential for large-scale industrial fermentation applications.

Data availability statement

The data presented in the study are deposited in NCBI Sequence Read Archive (SRA, https://www.ncbi.nlm.nih.gov/sra/PRJNA824254) database under the BioProject accession number PRJNA824254.

Author contributions

HL: Investigation, Formal Analysis, Writing – original draft. XF: Formal Analysis, Software, Writing – review and editing. SG: Writing – review and editing, Visualization. BZ: Formal Analysis, Writing – original draft. CW: Visualization, Writing – review and editing. QZ: Writing – review and editing, Project administration. RZ: Writing – review and editing, Visualization. HS: Writing – review and editing, Project administration. YQ: Writing – review and editing, Methodology, Funding acquisition. SZ: Methodology, Writing – review and editing. YJ: Funding acquisition, Supervision, Conceptualization, Writing – review and editing, Project administration.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study received support from the Key Research Project of Sichuan Provincial Branch of China National Tobacco Crop Tobacco Science Institute (SCYC202508) and the Key Research Project of China National Tobacco Crop Tobacco Science Institute (110202201036 (XJ-07)).

Conflict of interest

Authors QZ and RZ were employed by Deyang Tobacco Company of Sichuan Province.

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 author(s) declare that no Generative 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.

Supplementary material

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

References

Al-Mjeni, F., Ju, T., Pochapsky, T. C., and Maroney, M. J. (2002). XAS investigation of the structure and function of Ni in acireductone dioxygenase. Biochemistry 41, 6761–6769. doi:10.1021/bi012209a

PubMed Abstract | CrossRef Full Text | Google Scholar

Bond, J. S., and Beynon, R. J. (1995). The astacin family of metalloendopeptidases. Protein Sci. 4, 1247–1261. doi:10.1002/pro.5560040701

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, X., Mai, R., Zou, J., Zhang, H., Zeng, X., Zheng, R., et al. (2014). Analysis of aroma-active compounds in three sweet osmanthus (Osmanthus fragrans) cultivars by GC-olfactometry and GC-MS. J. Zhejiang Univ. Sci. B 15, 638–648. doi:10.1631/jzus.B1400058

PubMed Abstract | CrossRef Full Text | Google Scholar

Chassagnole, C., Raïs, B., Quentin, E., Fell, D. A., and Mazat, J. P. (2001). An integrated study of threonine-pathway enzyme kinetics in Escherichia coli. Biochem. J. 356, 415–423. doi:10.1042/bj3560415

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, S., Lu, Y., Wang, W., Hu, Y., Wang, J., Tang, S., et al. (2022). Efficient production of the β-ionone aroma compound from organic waste hydrolysates using an engineered Yarrowia lipolytica strain. Front. Microbiol. 13, 960558. doi:10.3389/fmicb.2022.960558

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890. doi:10.1093/bioinformatics/bty560

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Zhou, J., Ren, K., Zou, C., Liu, J., Yao, G., et al. (2019). Effects of enzymatic browning reaction on the usability of tobacco leaves and identification of components of reaction products. Sci. Rep. 9, 17850. doi:10.1038/s41598-019-54360-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Chunpin, X. (2015). Influence of complex enzymatic treatment on chemical composition and sensory quality of discarded tobacco leaf. Guizhou Agric. Sci. 1, 35–40. doi:10.3969/j.issn.1001-3601.2015.01.009

CrossRef Full Text | Google Scholar

Curien, G., Ravanel, S., Robert, M., and Dumas, R. (2005). Identification of six novel allosteric effectors of Arabidopsis thaliana aspartate kinase-homoserine dehydrogenase isoforms. J. Biol. Chem. 280, 41178–41183. doi:10.1074/jbc.M509324200

PubMed Abstract | CrossRef Full Text | Google Scholar

Czajka, J. J., Kambhampati, S., Tang, Y. J., Wang, Y., and Allen, D. K. (2020). Application of stable isotope tracing to elucidate metabolic dynamics during Yarrowia lipolytica α-Ionone fermentation. iScience 23, 100854. doi:10.1016/j.isci.2020.100854

PubMed Abstract | CrossRef Full Text | Google Scholar

de Lorenzo, L., Stack, T. M. M., Fox, K. M., and Walstrom, K. M. (2024). Catalytic mechanism and kinetics of malate dehydrogenase. Essays Biochem. 68, 73–82. doi:10.1042/EBC20230086

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi:10.1038/nmeth.2604

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, X., He, P., Pan, C., Xu, J., Xue, B., Yin, W., et al. (2019). Effects of pulsed electric field on the cell wall and infection activity of Rhizoctonia solani. Biol. (Basel) 8, 73. doi:10.3390/biology8040073

PubMed Abstract | CrossRef Full Text | Google Scholar

Hong-qi, Y. (2005). Effect on chiefly Latent fragrant substance on panel test flue-cured tobacco from different producing area. Journal of Hunan Agricultural University 1, 11–14. doi:10.3321/j.issn:1007-1032.2005.01.003

CrossRef Full Text | Google Scholar

Hoyles, L., Fernández-Real, J.-M., Federici, M., Serino, M., Abbott, J., Charpentier, J., et al. (2018). Molecular phenomics and metagenomics of hepatic steatosis in non-diabetic Obese women. Nat. Med. 24, 1070–1080. doi:10.1038/s41591-018-0061-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Jensen, C. O., and Parmele, H. B. (1950). Fermentation of cigar-type tobacco. Ind. Eng. Chem. 42, 519–522. doi:10.1021/ie50483a032

CrossRef Full Text | Google Scholar

Jia, Y., Zhou, W., Yang, Z., Zhou, Q., Wang, Y., Liu, Y., et al. (2023). A critical assessment of the candida strains isolated from cigar tobacco leaves. Front. Bioeng. Biotechnol. 11. doi:10.3389/fbioe.2023.1201957

PubMed Abstract | CrossRef Full Text | Google Scholar

Kikuchi, M., Hayashida, H., Nakano, E., and Sakaguchi, K. (1971). Peptidoglutaminase. Enzymes for selective deamidation of gamma-amide of peptide-bound glutamine. Biochemistry 10, 1222–1229. doi:10.1021/bi00783a019

PubMed Abstract | CrossRef Full Text | Google Scholar

Laureati, M., Cattaneo, C., Tateo, F., and Bononi, M. (2020). Identification of the Volatile Compounds and Sensory Attributes of Long-Term Aging Vin Santo Wine from Malvasia di Candia Aromatic Grapes. Foods 9, 1736. doi:10.3390/foods9121736

PubMed Abstract | CrossRef Full Text | Google Scholar

Lei, G., Wang, L., Liu, X., and Zhang, A. (2015). Fast quantification of phenylethyl alcohol in rose water and chemical profiles of rose water and oil of Rosa damascena and Rosa rugosa from southeast China. J. Liq. Chromatogr. and Relat. Technol. 38, 823–832. doi:10.1080/10826076.2014.976710

CrossRef Full Text | Google Scholar

Liu, F., Fang, X., Li, L., Tan, Z., Ye, M., Zhu, T., et al. (2022). Changes of volatile components in cigar wrapper during deep fermentation. Guangdong Agric. Sci. 49, 158–164. doi:10.16768/j.issn.1004-874X.2022.02.019

CrossRef Full Text | Google Scholar

Liu, F., Wu, Z., Zhang, X., Xi, G., Zhao, Z., Lai, M., et al. (2021). Microbial community and metabolic function analysis of cigar tobacco leaves during fermentation. MicrobiologyOpen 10, e1171. doi:10.1002/mbo3.1171

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, T., Guo, S., Wu, C., Zhang, R., Zhong, Q., Shi, H., et al. (2022). Phyllosphere microbial community of cigar tobacco and its corresponding metabolites. Front. Microbiol. 13. doi:10.3389/fmicb.2022.1025881

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W.-J., Lü, H.-Z., Li, Y., Yao, H., and Luo, H.-M. (2024). The new advance of terpene synthase research in the plant. PPJ 53, 1139–1149. doi:10.13592/j.cnki.ppj.2016.0513

CrossRef Full Text | Google Scholar

Magoč, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi:10.1093/bioinformatics/btr507

PubMed Abstract | CrossRef Full Text | Google Scholar

Meléndez-Martínez, A. J., Esquivel, P., and Rodriguez-Amaya, D. B. (2023). Comprehensive review on carotenoid composition: transformations during processing and storage of foods. Food Res. Int. 169, 112773. doi:10.1016/j.foodres.2023.112773

PubMed Abstract | CrossRef Full Text | Google Scholar

Ning, Y., Zhang, L.-Y., Mai, J., Su, J.-E., Cai, J.-Y., Chen, Y., et al. (2023). Tobacco microbial screening and application in improving the quality of tobacco in different physical states. Bioresour. Bioprocess. 10, 32. doi:10.1186/s40643-023-00651-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Penland, M., Deutsch, S.-M., Falentin, H., Pawtowski, A., Poirier, E., Visenti, G., et al. (2020). Deciphering microbial community dynamics and biochemical changes during nyons black olive natural fermentations. Front. Microbiol. 11, 586614. doi:10.3389/fmicb.2020.586614

PubMed Abstract | CrossRef Full Text | Google Scholar

Popova, V., Ivanova, T., Stoyanova, A., Nikolova, V., Hristeva, T., Gochev, V., et al. (2019). Terpenoids in the essential oil and concentrated aromatic products obtained from Nicotiana glutinosa L. leaves. Molecules 25, 30. doi:10.3390/molecules25010030

PubMed Abstract | CrossRef Full Text | Google Scholar

Rasouli, O., Ahmadi, N., Rashidi Monfared, S., and Sefidkon, F. (2018). Physiological, phytochemicals and molecular analysis of color and scent of different landraces of Rosa damascena during flower development stages. Sci. Hortic. 231, 144–150. doi:10.1016/j.scienta.2017.12.010

CrossRef Full Text | Google Scholar

Rimington, W. R., Pressel, S., Duckett, J. G., Field, K. J., Read, D. J., and Bidartondo, M. I. (2018). Ancient plants with ancient fungi: liverworts associate with early-diverging arbuscular mycorrhizal fungi. Proc. R. Soc. B Biol. Sci. 285, 20181600. doi:10.1098/rspb.2018.1600

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakulworakan, R., Chokmangmeepisarn, P., Dinh-Hung, N., Sivaramasamy, E., Hirono, I., Chuanchuen, R., et al. (2021). Insight into whole genome of Aeromonas veronii isolated from freshwater fish by resistome analysis reveal extensively antibiotic resistant traits. Front. Microbiol. 12, 733668. doi:10.3389/fmicb.2021.733668

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, A., Sivaraman, J., Li, Y., Larocque, R., Barbosa, J. A., Smith, C., et al. (2001). Three-dimensional structure of 2-amino-3-ketobutyrate CoA ligase from Escherichia coli complexed with a PLP-substrate intermediate: inferred reaction mechanism. Biochemistry 40, 5151–5160. doi:10.1021/bi002204y

PubMed Abstract | CrossRef Full Text | Google Scholar

Seeman, J. I., Laffoon, S. W., and Kassman, A. J. (2003). Evaluation of relationships between mainstream smoke acetaldehyde and “tar” and carbon monoxide yields in tobacco smoke and reducing sugars in tobacco blends of U.S. commercial cigarettes. Inhal. Toxicol. 15, 373–395. doi:10.1080/08958370304461

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Y., Cui, F., Wang, D., Li, T., and Li, J. (2021). Quorum quenching enzyme (PF-1240) capable to degrade AHLs as a candidate for inhibiting quorum sensing in food spoilage bacterium Hafnia alvei. Foods 10, 2700. doi:10.3390/foods10112700

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, J., Xie, D., Qi, D., Peng, Q., Chen, Z., Schreiner, M., et al. (2019). Methyl jasmonate-induced changes of flavor profiles during the processing of green, oolong, and black tea. Front. Plant Sci. 10, 781. doi:10.3389/fpls.2019.00781

PubMed Abstract | CrossRef Full Text | Google Scholar

Shibin, R., Jingjing, L., Yuanyuan, Z., Yanqing, Q., Jun, W., Xingyou, Y., et al. (2021). Effect of artificial fermentation temperature and humidity on chemical composition and aroma quality of cigar filler tobacco. Acta Tabacaria Sinica 27, 109–116. doi:10.16472/j.chinatobacco.2020.191

CrossRef Full Text | Google Scholar

Shin, S., Lv, J., Fazio, G., Mazzola, M., and Zhu, Y. (2014). Transcriptional regulation of ethylene and jasmonate mediated defense response in apple (Malus domestica) root during Pythium ultimum infection. Hortic. Res. 1, 14053. doi:10.1038/hortres.2014.53

PubMed Abstract | CrossRef Full Text | Google Scholar

Shoeb, M., Kodali, V. K., Farris, B. Y., Bishop, L. M., Meighan, T. G., Salmen, R., et al. (2017). Oxidative stress, DNA methylation, and telomere length changes in peripheral blood mononuclear cells after pulmonary exposure to metal-rich welding nanoparticles. NanoImpact 5, 61–69. doi:10.1016/j.impact.2017.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, W., Chen, X., Yu, J., Qiao, J., Yang, J., Chen, X., et al. (2024). Effects of Bacillus altitudinis inoculants on cigar tobacco leaf fermentation. Front. Bioeng. Biotechnol. 12, 1417601. doi:10.3389/fbioe.2024.1417601

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, X., Wang, J., Li, N., Yu, J., Meng, F., Wei, C., et al. (2020). Deciphering the high-quality genome sequence of coriander that causes controversial feelings. Plant Biotechnol. J. 18, 1444–1456. doi:10.1111/pbi.13310

PubMed Abstract | CrossRef Full Text | Google Scholar

Stackebrandt, E., and Goebel, B. M. (1994). Taxonomic note: a place for DNA-DNA reassociation and 16S rRNA sequence analysis in the present species definition in bacteriology. Int. J. Syst. Evol. Microbiol. 44, 846–849. doi:10.1099/00207713-44-4-846

CrossRef Full Text | Google Scholar

Tang, D., Shen, Y., Li, F., Yue, R., Duan, J., Ye, Z., et al. (2022). Integrating metabolite and transcriptome analysis revealed the different mechanisms of characteristic compound biosynthesis and transcriptional regulation in tea flowers. Front. Plant Sci. 13, 1016692. doi:10.3389/fpls.2022.1016692

PubMed Abstract | CrossRef Full Text | Google Scholar

Teng, Z. (2013). Canonical correlation between aroma aminophenol contents in flue-cured tobacco leaves and soil physicochemical properties. Soils 45, 277–284. doi:10.3969/j.issn.0253-9829.2013.02.014

CrossRef Full Text | Google Scholar

Tian, N., Liu, F., Wang, P., Yan, X., Gao, H., Zeng, X., et al. (2018). Overexpression of BraLTP2, a lipid transfer protein of Brassica napus, results in increased trichome density and altered concentration of secondary metabolites. Int. J. Mol. Sci. 19, 1733. doi:10.3390/ijms19061733

PubMed Abstract | CrossRef Full Text | Google Scholar

Vilela, A., Ferreira, R., Nunes, F., and Correia, E. (2020). Creation and acceptability of a fragrance with a characteristic tawny port wine-like aroma. Foods 9, 1244. doi:10.3390/foods9091244

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H., Guo, D., Zhang, M., Wu, G., Shi, Y., Zhou, J., et al. (2024). Correlation study on microbial communities and volatile flavor compounds in cigar tobacco leaves of diverse origins. Appl. Microbiol. Biotechnol. 108, 236. doi:10.1007/s00253-024-13032-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Xu, J., and Liu, A. (2022). Identification of the carotenoid cleavage dioxygenase genes and functional analysis reveal DoCCD1 is potentially involved in beta-ionone formation in Dendrobium officinale. Front. Plant Sci. 13, 967819. doi:10.3389/fpls.2022.967819

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, C., Yang, H., Wang, S., Zhao, J., Liu, C., Gao, L., et al. (2018). Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality. Proc. Natl. Acad. Sci. U. S. A. 115, E4151–E4158. doi:10.1073/pnas.1719622115

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, C., Shuanghong, Y., Wanrong, H., Ran, C., Zhishun, C., Danqun, H., et al. (2024). Effects of different adding methods of fermentation medium on the quality of cigar. Front. Bioeng. Biotechnol. 12, 1440961. doi:10.3389/fbioe.2024.1440961

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, C., Zhang, Q., Zhu, P., Hu, W., Jia, Y., Yang, S., et al. (2023). High throughput screening of key functional strains based on improving tobacco quality and mixed fermentation. Front. Bioeng. Biotechnol. 11. doi:10.3389/fbioe.2023.1108766

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Q., Peng, Z., Pan, Y., Liu, L., Li, L., Zhang, J., et al. (2023). Interaction analysis of tobacco leaf microbial community structure and volatiles flavor compounds during cigar stacking fermentation. Front. Microbiol. 14. doi:10.3389/fmicb.2023.1168122

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, X., Hu, Y., Wang, Q., Liu, J., Fang, S., Huang, D., et al. (2023). Study on the correlation between the dominant microflora and the main flavor substances in the fermentation process of cigar tobacco leaves. Front. Microbiol. 14. doi:10.3389/fmicb.2023.1267447

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, F., Yang, J., Luo, C., Li, D., Shi, G., Song, G., et al. (2023). Metagenomic insight into the biodegradation of biomass and alkaloids in the aging process of cigar. Bioresour. Bioprocess. 10, 45–14. doi:10.1186/s40643-023-00667-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Yuan, Y., Zhu, X., Xu, R., Shen, H., Zhang, Q., et al. (2022). The effect of different extraction methods on extraction yield, physicochemical properties, and volatile compounds from field muskmelon seed oil. Foods 11, 721. doi:10.3390/foods11050721

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L.-Y., Mai, J., Shi, J.-F., Ai, K.-B., He, L., Zhu, M.-J., et al. (2024). Study on tobacco quality improvement and bacterial community succession during microbial co-fermentation. Industrial Crops Prod. 208, 117889. doi:10.1016/j.indcrop.2023.117889

CrossRef Full Text | Google Scholar

Zhang, Q., Zheng, T., Yang, Z., Yang, S., Cai, W., Li, P., et al. (2023). Analysis of the structure and metabolic function of microbial community in cigar tobacco leaves in agricultural processing stage. Front. Microbiol. 14. doi:10.3389/fmicb.2023.1230547

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, T., Zhang, Q., Wu, Q., Li, D., Wu, X., Li, P., et al. (2022). Effects of inoculation with acinetobacter on fermentation of cigar tobacco leaves. Front. Microbiol. 13, 911791. doi:10.3389/fmicb.2022.911791

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Q., Yang, J., Feng, Y., Yang, Z., Wang, Y., Zhang, Z., et al. (2024). Analysis of the effects of Bacillus velezensis HJ-16 inoculation on tobacco leaves based on multi-omics methods. Front. Bioeng. Biotechnol. 12, 1493766. doi:10.3389/fbioe.2024.1493766

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: tobacco cigar leaves, fermentation, microbial community, functional strain screening, bioaugmentation

Citation: Liu H, Fang X, Guo S, Zhang B, Wu C, Zhong Q, Zhang R, Shi H, Qin Y, Zeng S and Jin Y (2025) Isolation and evaluation of a Bacillus altitudinis strain to improve cigar tobacco leaves fermentation. Front. Bioeng. Biotechnol. 13:1636506. doi: 10.3389/fbioe.2025.1636506

Received: 28 May 2025; Accepted: 24 July 2025;
Published: 12 August 2025.

Edited by:

Ying-Chien Chung, China University of Science and Technology, Taiwan

Reviewed by:

Digvijay Verma, Babasaheb Bhimrao Ambedkar University, India
Chandrabhan Seniya, Manipal University Jaipur Faculty of Science Technology and Architecture, India

Copyright © 2025 Liu, Fang, Guo, Zhang, Wu, Zhong, Zhang, Shi, Qin, Zeng and Jin. 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: Yanqing Qin, cXF5eXFxMjAwMEAxNjMuY29t; Shuhua Zeng, enNoZ3NwQDE2My5jb20=; Yao Jin, eWFvamluMTJAc2N1LmVkdS5jbg==

Disclaimer: 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.