Deep sequencing reveals changes in prokaryotic taxonomy and functional diversity of pit muds in different distilleries of China
Annals of Microbiology volume 72, Article number: 12 (2022)
The microbial community in the pit mud correlated closely with the quality of the final product of Chinese strong-flavored Baijiu (CSFB). However, environmental conditions and brewing processes can vary by region and distilleries. This may lead to differences in microbial composition and function in pit mud. Therefore, revealing the features of the pit mud microbial community structure and functions of different distilleries will provide key information for understanding the diversity and difference of microbes in the brewing of CSFB, which will be beneficial for the improvement of the quality of pit mud and CSFB in the future.
Methods and results
Illumina MiSeq sequencing of 16S rRNA gene amplicons was used to analyze the similarities and differences in microbial community structure and function in pit muds of different distilleries located in Shihezi (Xinjiang), Xiangyang (Hubei), and Yibin (Sichuan). At the genus level, Clostridium, Lactobacillus, Aminobacterium, Petrimonas, Syntrophomonas, Methanoculleus, Syntrophaceticus, Sedimentibacter, Caloramator, Ruminococcus, Bacillus, Methanosarcina, and Garciella were the dominated genera of pit muds. There were great differences in the composition of microorganisms in pit muds used by different distilleries. The significantly enriched prokaryotic microbiotas of pit muds collected in the distilleries of Xiangyang were mainly affiliated with Bacillus, Lactobacillus, and Croceifilum, and the relative abundance of methanogens, such as Methanomicrobia and Methanobacteria, were only significantly enriched in the pit mud collected from the distilleries of Yibin (P < 0.05). Functional analysis indicated that the difference of microbial composition in pit mud will further lead to significant differences in various metabolic functions.
The compositions and functions of dominant microorganisms in pit mud used for the production of CSFB by different enterprises across regions in China were greatly different, and there was a close relationship between the compositions and functions of microorganisms in pit mud. Therefore, it may be an effective method to improve CSFB fermentation processes by directionally regulating the microbial community functions of pit mud using specific strains.
Chinese liquor, vodka, rum, tequila, brandy, and whiskey are the six most popular distilled liquors of the world. Of them, Chinese strong-flavored Baijiu (CSFB) is one of the most important traditional distillates, accounting for >70% of the total liquor consumption in China (Liu et al. 2017). The raw materials for making CSFB include grains, typically sorghum or a mixture of corn, rice, millet, and wheat (Zheng and Han 2016). CSFB is fermented under anaerobic conditions in cellars (approximate dimensions: length, 2.0–3.5 m; width, 2.0–3.0 m; and depth, 2.3–2.5 m) lined with pit mud, and the process takes 60–90 days. The prokaryotes present in pit mud break down the macromolecules from the raw material into peptides and monosaccharides, producing aromatic compounds responsible for the flavor of the product (Li et al. 2011; Wang et al. 2017; Xiang et al. 2019). Both recent studies and long-term production practices have revealed that the microbial community in pit mud correlates closely with the quality of the final product (Tao et al. 2017; Liu et al. 2020). However, environmental conditions and brewing processes can vary by region and distilleries. This may lead to differences of microbial composition in pit mud. Therefore, revealing the features of the pit mud microbial community structure and functions of different distilleries will provide key information for understanding the diversity and difference of microbes in CSFB brewing, which will be helpful for the improvement of CSFB production technology in the future.
Classical microbiological methods have provided early insight into the predominant bacteria involved in CSFB fermentation. The development of the next-generation DNA sequencing technologies has revolutionized microbial ecology studies in recent years by allowing high-throughput analysis. Technological progress has allowed more efficient and accurate identification of the microbial communities and functional capacities in different niches, including Daqu (a Chinese traditional fermentation starter) (Xie et al. 2020), soil (Fan et al. 2020), and parts of the body, such as the gut (Hou et al. 2020) and skin (Rainer et al. 2020). At present, this technology has also been applied in some studies related to pit mud microbes (Tao et al. 2017; Chen et al. 2020). However, these studies often collect pit mud samples from specific distilleries located in the same region. As a result, the commonalities and differences of microorganisms in pit mud of different distilleries across regions in China have not been reported yet.
In this study, pit mud samples used for CSFB production were collected from different distilleries located in Shihezi (Xinjiang Province) and Xiangyang (Hubei Province). The community structure and functional characteristics of prokaryotes present in pit muds were analyzed by using Illumina MiSeq sequencing technology. In addition, we integrated and examined sequencing data of pit mud samples taken from Yibin (Sichuan Province) (http://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA597727) (Chen et al. 2020), a region known for the production of high-quality CSFB in China. This study could deepen our understanding of the compositions and functional characteristics of the prokaryotic community in pit mud from different distilleries in China and help further improve the quality of pit mud and CSFB.
Materials and methods
Pit mud samples were collected from Shihezi and Xiangyang. Shihezi is far from the ocean, and air from the ocean is does easily to reach this region. Therefore, there is less precipitation and a dry climate. Xiangyang, in turn, is characterized by the transition from north to south and the consideration of east and west. It belongs to the north subtropical monsoon climate with moderate humidity.
The CSFB distilleries of Shihezi and Xiangyang selected for this study had been established for a long time. However, on the whole, the fermentation cellars of distilleries in Xiangyang collected in this study were younger than those of Shihezi. The collected samples in these fermentation cellars showed characteristics typical of normal pit mud. It does not show degradation of earth color or ester aroma, and it has a moist, soft, and uniform texture. Samples were collected according to the stratified random method (Zheng et al. 2013), and the sampling sites were based on the method adopted by a previous study (Chen et al. 2020). In brief, samples were collected from three depths in each cellar (upper layer, 10 cm below the cellar surface; middle layer, 1.1 m below the cellar surface; and lower layer, collected after removing a 2–3-cm layer of pit mud from the cellar bottom). The sampling procedure is illustrated in Fig. S1. A total of 30 samples were collected from ten distilleries in Shihezi and 27 samples were collected from nine distilleries in Xiangyang. Two subsamples were collected from each layer using a soil core sampler (AMS Inc., USA); 100 g of each subsample was obtained and thoroughly mixed. All samples were kept in ice boxes and transported from the distilleries to the laboratory at low temperatures (transportation time <48 h). In the laboratory, they were stored at −20°C until further analysis.
DNA was extracted from each sample using the FastDNA SPIN kit for soil (Mpbio, USA) following the manufacturer’s instructions. The quality of extracted DNA was checked using 1% agarose gel electrophoresis (Beijing Liuyi Biological Technology Co., Ltd. China) and spectrophotometry (ratio of optical density at 260 and 280 nm) using a NanoDrop™ One spectrophotometer (Thermo Fisher Scientific, USA). All extracted DNA samples were stored at −20°C until further analysis.
The V3–V4 region of the 16S ribosomal RNA (rRNA) gene was amplified, as previously described (Chen et al. 2020), using the forward 341F (5′-CCTAYGGGRBGCASCAG-3′) and reverse 806R (5′-GGACTACHVGGGTWTCTAAT-3′) primers (Yu et al. 2005). The primers contained a set of seven nucleic acid barcodes, which were used to distinguish samples in subsequent analyses. The PCR involved the following steps: 2 min at 95°C; 30 cycles of 20 s at 95°C, 30 s at 55°C, and 30 s at 72°C; final extension at 72°C for 5 min; and cooling at 4°C. The amplicons were sequenced at Majorbio (Shanghai, China) using the Illumina MiSeq system.
Primers and barcode sequences were removed from high-quality reads via in-house Python scripts. Furthermore, based on barcodes, the reads were divided into different sample groups. After merging all sequencing data, including the data of pit mud samples taken from Yibin (Chen et al. 2020), the Quantitative Insights Into Microbial Ecology (QIIME) package (version 1.9.1) (Caporaso et al. 2010) was applied to perform bioinformatics analysis. The specific steps were similar to those from our previous study (Hou et al. 2020). Briefly, UCLUST (Edgar 2010) was employed to classify high-quality sequences into operational taxonomic units (OTUs) at an identity threshold of 97% similarity, while ChimeraSlayer (Haas et al. 2011) was used to remove potential chimeric sequences from the OTU representative set. Singleton OTUs (OTUs with only one sequence) were removed from all datasets. Based on the information extracted from Greengenes (version 13.8) (Desantis et al. 2006), the Ribosomal Database Project (RDP, Release 11.5) (Cole et al. 2007), and SILVA (Version 138) (Quast et al. 2013), each OTU was assigned to the lowest possible taxonomic level on the basis of a minimum bootstrap threshold of 80% (Hou et al. 2015). The OTU table was subsampled correspondingly to adjust the sampling depth using the “multiple_rarefactions.py” program from the QIIME pipeline. Alpha- and beta-diversity were calculated based on the de novo taxonomic tree constructed by the representative chimera-checked OTU set using FastTree (Price et al. 2009). The α-diversity indices, such as the number of observed OTUs and the Shannon diversity index, were calculated in 1000 iterations of rarefied OTU tables in QIIME with the minimum sequencing depth among all study subjects. The functional profiles of prokaryotes in pit mud samples were predicted based on Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) (Langille et al. 2013) software and the COG database.
Statistical analyses were performed using R software version 4.0.2 (https://www.r-project.org/). Nonparametric Mann–Whitney or Kruskal–Wallis tests were used to ascertain the statistical significance of between-group differences in α-diversity indices, prokaryotic composition, and functional categories. False discovery rate values were estimated using the Benjamini-Yekutieli method to control for multiple testing (Benjamini and Yekutieli 2001). P values less than 0.05 were considered statistically significant. To evaluate the prokaryotic community structure in different samples, principal coordinate analysis (PCoA) was carried out based on both the weighted and unweighted UniFrac distances derived from the phylogenetic tree (Lozupone and Knight 2005). The extent of compositional differences between groups was determined using permutational multivariate analysis of variance (PERMANOVA) (Anderson 2001) tests based on UniFrac distances, which examine the significance of inter- and intra-group variations (pseudo F-statistic) via permutations of group assignment. Canonical analysis of principal coordinates (CAP) was conducted using MATLAB software version 2016. The linear discriminant analysis effect size (LEfSe) method was used to determine taxa of prokaryotic biomarkers among the pit mud samples from different regions (LDA score ≥ 3) (Segata et al. 2011). Procrustes analysis was employed to reveal the correlation between the microbial community and functional profiles. Figures were plotted using the R ggplot2 (Wickham 2016) and gplots packages (https://github.com/talgalili/gplots).
Sequence coverage of prokaryotic communities across samples from Shihezi and Xiangyang
We generated a dataset that included 2,127,407 high-quality reads of the V3–V4 regions of prokaryotic 16S rRNA genes present in samples from Shihezi and Xiangyang; there were 37322 ± 7789 high-quality reads (mean ± SD, range 16,519–54,860) per sample. At a high identity cutoff level of 97% sequence similarity, 67,896 OTUs were detected. After removing singleton OTUs, the average number of OTUs per sample was 2484 (range, 712–5496; SD=1243). Based on the information extracted from the RDP, Greengenes, and SILVA databases, each OTU was assigned to the lowest taxonomic level using homologous sequence alignment and clustering. Overall, 5.54% of the sequences could not be assigned to any genus level. Quantification of α-diversity revealed that the Shannon–Wiener diversity curves (rather than the rarefaction curves) reached saturation (Fig. S2), indicating that the obtained sequence depth for all samples was sufficient, although other new phylotypes may be identified by further sequencing.
Prokaryotic composition of all pit mud samples
The prokaryotic community structure of pit mud samples was analyzed at the phylum and genus levels. In total, 68 phyla and 1044 genera were identified from all pit mud samples. Six phyla showed a relative abundance >1%: Firmicutes (76.0%), Bacteroidetes (7.5%), Euryarchaeota (6.3%), Spirochaetes (3.8%), Actinobacteria (2.5%), and Proteobacteria (1.7%) (Fig. 1a). Firmicutes appeared to be the most prevalent phylum in all samples, while Euryarchaeota were mainly present in samples collected from Yibin.
At the genus level, there were obvious differences among the dominant genera in pit muds collected from different regions and distilleries (Fig. 1b). Specifically, there were fifteen and thirteen genera with an average relative abundance >1% in pit muds of Shihezi and Yibin, respectively. However, there were only eight genera with an average relative abundance >1% in pit muds of Xiangyang. The predominant genera in pit mud samples from all distilleries were as follows: Lactobacillus (19.1%), Clostridium (16.1%), Aminobacterium (3.9%), Petrimonas (3.6%), Syntrophomonas (3.2%), Methanoculleus (2.7%), Syntrophaceticus (2.4%), Sedimentibacter (2.3%), Caloramator (2.2%), Ruminococcus (1.8%), Bacillus (1.4%), Methanosarcina (1.1%), and Garciella (1.0%). Interestingly, Methanoculleus and Methanosarcina were mainly present in pit mud samples from Yibin.
Presence of core fermentation-associated prokaryotes in pit mud samples
To identify the core prokaryotes that were most likely to be involved in fermentation, the OTUs that were shared by ≥95% of pit mud samples in at least one region were defined as core OTUs (cOTUs). Overall, 129 OTUs met this threshold (Fig. S3). Among them, 18 cOTUs were found in all samples, of which nine showed annotation results for Clostridium; other cOTUs showed annotation results for Lactobacillus, Ruminococcus, and Bacillus. The cumulative relative abundances of cOTUs in the upper, middle, and lower layers of pit mud from Shihezi were 1.9%, 3.3%, and 19.4%, respectively, while in Xiangyang, these values were much higher (49.3%, 59.9%, and 70.4%, respectively). However, the proportions of cOTUs in the middle and lower layers of pit mud from Yibin were 0.3% and 0.4%, respectively. Notably, many cOTUs that were only found in pit mud samples from Shihezi and Yibin but were missing from Xiangyang samples showed annotations for Sedimentibacter, Caldicoprobacter, Tissierella, Syntrophaceticus, and Petrimonas. Nevertheless, pit mud samples collected from Xiangyang had many unique cOTUs, several of which had annotations for Lactobacillus, Clostridium, Bacillus, Micromonosporaceae_unclassified, Ralstonia, Thermoactinomyces, and Hydrogenispora. Therefore, the community structure of prokaryotes in Xiangyang pit mud collected in this study was more special.
Differences in α-diversity based on distilleries and depth
Prokaryotic species richness and diversity were assessed based on the number of observed OTUs and the Shannon diversity index, respectively (Fig. 2). Pit mud samples collected from Yibin had the highest α-diversity, followed by those from Shihezi and Xiangyang; the α-diversity indices were significantly different among the regions (P<0.05). In addition, the α-diversity increased with increasing depth in pit mud samples from Xiangyang, and the α-diversity was significantly higher in lower-layer pit mud than in upper-layer pit mud (P<0.05). There was no significant difference in the α-diversity between the middle and lower layers of pit mud samples from Yibin.
Comparison of prokaryotic community structure in pit mud samples
To evaluate the variation in prokaryotic community structure in pit mud samples from different distilleries and depths, we used PCoA based on unweighted and weighted UniFrac distances (Fig. 3a, b). When using unweighted UniFrac distances, score plots generated from principal components 1 and 2 showed a clear clustering pattern associated with the collected regions. PCoA based on weighted UniFrac distances showed that the pit mud samples from Xiangyang were clearly distinct from those of Shihezi and Yibin, but that there was some overlap between samples from the latter two regions. There was also a clear overlap between samples obtained from different depths in the same region.
To further compare the influence of distilleries and sampling depth on prokaryotic community structure, the sampling regions and depths were used as grouping factors in PERMANOVA. The results showed that the calculated F value was greater when grouping was performed according to sampling regions (F = 36.75) than when it was performed according to sampling depth (F = 6.35). This indicated that the influence of sampling distilleries on prokaryotic community structure was greater than that of sampling depths. These results were confirmed by calculations performed using CAP (Fig. 3c), as samples from the same region were found to be clustered together.
Differences in the microbial composition of pit mud samples from different distilleries and depths
To compare the differences in microbial composition in distilleries located in different regions, LEfSe analysis was used to identify the taxonomic prokaryotes with significant abundance differences among the three regions (Fig. 4). Due to the lack of upper pit mud samples in Yibin, the analysis was only based on the middle and lower pit mud samples in the region. The significantly enriched prokaryotic microbiotas in Yibin were mainly affiliated with Methanomicrobia and Methanobacteria of Euryarchaeota. The significantly enriched prokaryotic microbiotas in Xiangyang were mainly affiliated with Bacillus, Lactobacillus, and Croceifilum. The significantly enriched prokaryotic microbiotas in Shihezi were mainly affiliated with Clostridium, Ruminococcus, Fermentimonas, Petrimonas, Hazenella, Aminobacterium, Syntrophomonas, Garciella, and Sedimentibacter of the classes Bacteroidia, Bacilli, Clostridia, and Spirochaetia. Based on the Kruskal–Wallis tests, we found that the unclassified prokaryotes in Yibin accounted for 30%, which was significantly higher than those in Shihezi and Xiangyang.
In order to reveal the influence of depth on prokaryotic composition, we compared the relative abundance of groups present at an abundance >0.1% in samples obtained from Shihezi and Xiangyang using paired Kruskal–Wallis tests (Table 1). For samples from Shihezi, the abundance of 21 genera varied significantly across different depths. The genera showing a decrease in relative abundance with an increase in depth included Lactobacillus, Hazenella, Bacillus, Rummeliibacillus, Paenibacillus, Dehalobacterium, and Desulfotomaculum, while those that showed an increase in relative abundance with increasing depth included Aminobacterium, Cloacibacillus, and Gracilibacter. The relative abundance of other prokaryotic genera did not change regularly with depth. In Xiangyang, only the relative abundance of three genera changed significantly with increasing depth. Among these genera, the relative abundance of Lactobacillus decreased, while that of Clostridium and Ruminococcus increased with increasing depth. Using Mann–Whitney tests, we found that the relative abundance of Lactobacillus, Tissierella, Nocardia, Caloribacterium, and Ralstonia in the lower-layer samples from Yibin was significantly higher than that in the corresponding middle-layer samples (P<0.05).
Functional analysis of the microbial community in pit mud samples
Using PICRUSt software and the COG database, 4675 putative COGs were identified from all pit mud samples in the present study. In addition to the poorly characterized categories R and S, COG categories related to cell growth, such as G (carbohydrate transport and metabolism), J (translation, ribosomal structure, and biogenesis), L (replication, recombination, and repair), E (amino acid transport and metabolism), and K (transcription), were found to be predominant across all samples. Comparative analysis of the relative abundance of COG categories in samples from the same depth in each region (Fig. 5a) showed that COG categories Z (cytoskeleton) and U (intracellular trafficking, secretion, and vesicular transport) were significantly more abundant in samples collected from distilleries of Shihezi than in those from the other two regions. In addition, COG category R at all depths and categories L and F (nuclear transport and metabolism) in the upper and middle layers were significantly more abundant in samples from Xiangyang than in those from the other two regions. Yibin samples were significantly more enriched for the COG categories N (cell motility), C (energy production and conversion), H (coenzyme transport and metabolism), B (chromatin structure and dynamics), and S than those from the other two regions. It is worth noting that the relative abundance of COG functional categories N, C, H, and B in samples from Xiangyang was the lowest among the three regions.
Paired Mann–Whitney or Kruskal–Wallis tests were used to compare COG functional categories in samples from different depths in the three regions (Fig. 5b). For Shihezi, functional category U increased significantly in abundance as sample depth increased, and functional category I decreased with increasing depth. For Xiangyang, the relative abundance of functional categories N, H, O (posttranslational modification, protein turnover, and chaperones), T (signal transduction mechanisms), C, and Q (secondary metabolites biosynthesis, transport, and catabolism) significantly increased with increasing sample depth, while that of functional categories L, S, R, G, F, and J decreased with increasing depth. For Yibin, we found that the abundance of functional categories K and I was significantly higher in the lower layer, while that of categories C and S was significantly higher in the middle layer.
Relationship between prokaryotic microbiotas and functions in pit mud samples
To identify co-occurrence patterns among the dominant genera (i.e., those with an average relative abundance >1%) and potential functions of particular genera, we used Spearman’s rank correlation analysis on abundance matrixes of selected prokaryotic genera and COG functional categories (Fig. 6a, b). There were close correlations between prokaryotic genera. In addition, prokaryotic genera also have close correlations with COG functional categories. Notably, the dominant genera in pit mud samples could be divided into two groups: one represented by Lactobacillus and Bacillus and the other containing all the remaining dominant genera. There was a significant negative correlation between the abundance of the two groups, but the abundance of genera within the same group showed a positive correlation with each other (Fig. 6a). There was also a significant negative correlation between the relative abundance of the first group and the COG categories B, U, T, N, O, H, and C, while there was a significant positive correlation with the COG categories L, I, E, Q, V, and K (Fig. 6b). The correlation between the relative abundance of the second group and the abovementioned COG functional units was the complete opposite to the relationship observed for the first group. Procrustes analysis was employed to reveal the consistency of two-dimensional shapes produced by the superimposition of ordination analyses from taxonomic abundance and functional profiles, and it turned out to be extremely significantly related (P = 0.001) (Fig. 6c).
In the present study, we revealed that the dominant phylum in pit mud was Firmicutes, which is consistent with previous reports (Zheng et al. 2013; Liang et al. 2016). Clostridium and Lactobacillus were the dominant genera in all pit mud samples. However, the relative abundance of these two genera was quite different among samples from different distilleries. In addition, half of the cOTUs that were shared by samples from the three regions were also annotated as Clostridium, indicating that the genus exists widely in these regions. Clostridium are strictly anaerobic and play a key role during liquor production (Hu et al. 2016). First, they influence interspecies hydrogen transfer and make a synergistic contribution to metabolism, thereby maintaining the stability of the microbial community (Zheng et al. 2013). Second, they can synthesize caproic acid using sugar, starch, and cellulose as substrates (Dürre 2016; Cavalcante et al. 2017). After reacting with ethanol, caproic acid could form ethyl butyrate and ethyl caproate, the main flavor components of CSFB, through enzyme catalysis and nonenzyme catalysis, respectively (Hu et al. 2015). Therefore, a certain abundance of Clostridium spp. in pit mud is very important for ensuring the excellent flavor of CSFB products.
As another dominant group, Lactobacillus also plays an important role in CSFB production. Lactic acid, the predominant metabolite of Lactobacillus, has been identified as the precursor of ethyl lactate, which is also a major flavoring compound in CSFB (Fan and Qian 2005; Yang et al. 2020). In addition, Lactobacillus reduces the unfavorable characteristics of liquor, as lactic acid reduces the sensory impact of ethanol, enhancing its mellowness and prolonging its aftertaste (Yang et al. 2019). However, an overabundance of Lactobacillus in fermented grains can lead to the production of excess lactic acid and, consequently, excess ethyl lactate. Once the ethyl lactate content surpasses that of ethyl hexanoate, the aroma is reduced, the main body flavor is suppressed, and the flavor loses its balance, with a stuffy sweetness appearing (Li 2010). It is worth noting that Lactobacillus was present in high abundance in all samples from the distilleries of Xiangyang, irrespective of the depth of collection, and this was verified based on the results of cOTUs. This may be related to the fact that the fermentation cellars in Xiangyang collected in this study were relatively young. To improve the quality and flavor of CSFB in Xiangyang, more attention should be given to the richness of Lactobacillus in pit mud.
Bacillus was present in pit mud samples from all three regions. It was the third richest genus in the samples collected from Xiangyang distilleries, and its richness was obviously higher than that of pit mud samples collected from Shihezi and Yibin. Bacillus can produce amylase, glucoamylase, and cellulase (Li et al. 2014; Yang et al. 2017), which degrade cellulose, starch, and protein into substances that can be used for subsequent fermentation. In addition, Bacillus in pit mud can produce specific flavor compounds, such as 2,3-butanediol, 3-hydroxy-2-butanone, 2-methylpropionic acid, and 3-methylbutanoic acid, which are beneficial and improve the flavor of the final product (Zhang et al. 2013; Meng et al. 2015).
The results for α-diversity in this study confirmed that there were significant differences in the community structure of prokaryotes in samples from different distilleries and that samples collected from Yibin and Xiangyang had the highest and lowest α-diversity, respectively. For various macroscale ecosystems, it is well established that high biodiversity increases the essential resource use efficiency of ecological communities and maintains the balance and stability of the ecosystem (Cardinale et al. 2012). Yibin enjoys a mid-subtropical humid monsoon climate (Zhao et al. 2013), which is very suitable for the growth of liquor-making prokaryotes. This is probably why samples from Yibin show the highest prokaryotic α-diversity. In contrast, the relatively low α-diversity in Xiangyang might be closely related to the high relative abundance of Lactobacillus and Bacillus, which were found to be negatively correlated with that of other dominant genera.
The environment (oxygen, moisture, temperature, and pressure) within the pit mud at different depths in fermentation cellars is not the same, which may lead to differences in prokaryotic composition. Therefore, prokaryotic community structures of pit muds from different depths were collected and analyzed separately in this study. By calculating the influence of distilleries and sampling depth on prokaryotic community structure, we showed that the distillery had a greater influence than did the sampling depth. This may be related to the cellar age and brewing technology, and the oxygen content and mixing ratio of grain in different enterprises may be different (Xiang et al. 2015; Xu et al. 2017). Together, these factors could lead to vast differences in the prokaryotic community structure across different enterprises.
We also examined the genera that showed a significantly different abundance across the distilleries. Specifically, among the genera enriched in samples from Shihezi, Sedimentibacter, and Aminobacterium had similar functions, as both ferment amino acids into acetic and butyric acid (Imachi et al. 2016). In addition, Petrimonas produces acetic acid and H2 as the main end products of carbohydrate fermentation, and Fermentimonas produces acetic acid, propionic acid, CO2, traces of isovaleric acid, and H2 (Grabowski et al. 2005; Hahnke et al. 2016). It is worth noting that several methanogenic genera, such as Methanomicrobia and Methanobacteria, were observed in samples from Yibin but were rarely found in samples from Shihezi and Xiangyang. Previous studies have found that specific members of Clostridium, Bacteroidia, Methanobacteria, and Methanomicrobia microorganisms have synergistic effects on anaerobic digestion (Barberan et al. 2012), and more caproate is produced when they are cocultured (Liu et al. 2003), which is beneficial for liquor quality.
In terms of the main metabolic functions of the prokaryotes in pit mud, the COG functional categories with the highest relative abundance across regions, including categories G, J, L, E, K, and R, were related to cell growth (Jeong et al. 2018). The pit mud environment is rich in various microbial species, and prokaryotes increase their relative abundance in the environment through continuous proliferation, to compete for available nutrients. In this process, functional metabolism related to cell growth and reproduction becomes more active in microbes.
There was a close positive correlation between the abundance of Lactobacillus and Bacillus. This may be related to the strong ability of Bacillus to resist acidic conditions (Elshaghabee et al. 2017). It is also worth noting that the existence of Lactobacillus and Bacillus in pit mud is negatively correlated with the existence of almost all other dominant prokaryotes, which indicates growth competition between the two groups. The results of the correlation analysis between the microbial community and functional capacities showed that Lactobacillus and Bacillus had significant negative correlations with COG functional modules such as B, C, H, N, and U, which are necessary for the normal growth of microorganisms, and were significantly positively correlated with functional modules related to the biosynthesis and defense mechanism of microbial secondary metabolites. This suggests that the excessive presence of these two genera in pit mud may interfere with the normal metabolic activities of pit mud microbes and lead to intense competition between pit mud microbes. Procrustes analysis further revealed the close relationship between the taxonomic abundance and functional distribution of microbial communities in pit mud. This suggests that the directional regulation of pit mud prokaryotic community functions using specific strains might become an effective method of improving CSFB fermentation processes and lead to the production of uniformly high-quality products in the future.
In this study, the microbial community and functional capacities of pit mud from distilleries located in three different regions of China were analyzed by Illumina sequencing of 16S rRNA gene amplicons. Clostridium, Lactobacillus, Aminobacterium, Petrimonas, Syntrophomonas, Methanoculleus, Syntrophaceticus, Sedimentibacter, Caloramator, Ruminococcus, Bacillus, Methanosarcina, and Garciella were the dominant genera of pit muds. The relative abundance of those genera in pit mud from different distilleries were quite different. Bacillus, Lactobacillus, and Croceifilum were significantly enriched in pit muds collected in Xiangyang, and the relative abundance of methanogens, such as Methanomicrobia and Methanobacteria, was only significantly enriched in pit mud collected from Yibin. In addition, differences in microbial composition in pit mud will further lead to significant differences in various metabolic functions. Therefore, it may be an effective method to improve CSFB fermentation processes by directionally regulating the microbial community functions of pit mud using specific strains. It should be pointed out that CSFB are widely distributed in various regions of China. Because some trade secrets are involved, it was difficult for us to obtain pit mud samples from a large number of distilleries. Therefore, the results of this study cannot represent the microbial and functional composition of pit mud samples from all distilleries in the regions of Xiangyang, Shihezi, and Yibin.
Availability of data and materials
Raw sequence data are publicly available online through the MG-RAST project number mgp96217 (http://metagenomics.anl.gov/linkin.cgi?project=mgp96217).
Anderson MJ (2001) A new method for non-parametric multivariate analysis of variance. Austral Ecol 26:32–46. https://doi.org/10.1046/j.1442-9993.2001.01070.x
Barberan A, Bates ST, Casamayor EO, Fierer N (2012) Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J 6:343–351. https://doi.org/10.1038/ismej.2011.119
Benjamini Y, Yekutieli D (2001) The control of the false discovery rate in multiple testing under dependency. Ann Stat:1165–1188. https://doi.org/10.1214/aos/1013699998
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, Mcdonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R (2010) QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7:335–336. https://doi.org/10.1038/nmeth.f.303
Cardinale BJ, Duffy JE, Gonzalez A, Hooper DU, Perrings C, Venail P, Narwani A, Mace GM, Tilman D, Wardle DA, Kinzig AP, Daily GC, Loreau M, Grace JB, Larigauderie A, Srivastava DS, Naeem S (2012) Biodiversity loss and its impact on humanity. Nature 486:59–67. https://doi.org/10.1038/nature11148
Cavalcante WDA, Leitao RC, Gehring TA, Angenent LT, Santaella ST (2017) Anaerobic fermentation for n-caproic acid production: a review. Process Biochem 54:106–119. https://doi.org/10.1016/j.procbio.2016.12.024
Chen L, Li Y, Jin L, He L, Ao X, Liu S, Yang Y, Liu A, Chen S, Zou L (2020) Analyzing bacterial community in pit mud of Yibin Baijiu in China using high throughput sequencing. PeerJ 8. https://doi.org/10.7717/peerj.9122
Cole JR, Chai B, Farris RJ, Wang Q, Kulam-Syed-Mohideen AS, Mcgarrell DM, Bandela AM, Cardenas E, Garrity GM, Tiedje JM (2007) The ribosomal database project (RDP-II): introducing myRDP space and quality controlled public data. Nucleic Acids Res 35:D169–D172. https://doi.org/10.1093/nar/gkl889
Desantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL (2006) Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol 72:5069–5072. https://doi.org/10.1128/AEM.03006-05
Dürre P (2016) Physiology and sporulation in Clostridium. The bacterial spore: from molecules to systems, pp 313–329. https://doi.org/10.1128/microbiolspec.TBS-0010-2012
Edgar RC (2010) Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26:2460–2461. https://doi.org/10.1093/bioinformatics/btq461
Elshaghabee FMF, Rokana N, Gulhane RD, Sharma C, Panwar H (2017) Bacillus as potential probiotics: status, concerns, and future perspectives. Front Microbiol 8:1490. https://doi.org/10.3389/fmicb.2017.01490
Fan WL, Qian MC (2005) Headspace solid phase microextraction and gas chromatography-olfactometry dilution analysis of young and aged Chinese “Yanghe Daqu” liquors. J Agric Food Chem 53:7931–7938. https://doi.org/10.1021/jf051011k
Fan X, Chang W, Sui X, Liu Y, Song G, Song F, Feng F (2020) Changes in rhizobacterial community mediating atrazine dissipation by arbuscular mycorrhiza. Chemosphere 256:127046. https://doi.org/10.1016/j.chemosphere.2020.127046
Grabowski A, Tindall BJ, Bardin V, Blanchet D, Jeanthon C (2005) Petrimonas sulfuriphila gen. nov, sp nov, a mesophilic fermentative bacterium isolated from a biodegraded oil reservoir. Int J Syst Evol Microbiol 55:1113–1121. https://doi.org/10.1099/ijs.0.63426-0
Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, Ciulla D, Tabbaa D, Highlander SK, Sodergren E (2011) Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res 21:494–504. https://doi.org/10.1101/gr.112730.110
Hahnke S, Langer T, Koeck DE, Klocke M (2016) Description of Proteiniphilum saccharofermentans sp nov, Petrimonas mucosa sp nov and Fermentimonas caenicola gen. nov, sp nov, isolated from mesophilic laboratory-scale biogas reactors, and emended description of the genus Proteiniphilum. Int J Syst Evol Microbiol 66:1466–1475. https://doi.org/10.1099/ijsem.0.000902
Hou Q, Xu H, Zheng Y, Xi X, Kwok LY, Sun Z, Zhang H, Zhang W (2015) Evaluation of bacterial contamination in raw milk, ultra-high temperature milk and infant formula using single molecule, real-time sequencing technology. J Dairy Sci 98:8464–8472. https://doi.org/10.3168/jds.2015-9886
Hou Q, Zhao F, Liu W, Lv R, Khine WWT, Han J, Sun Z, Lee YK, Zhang H (2020) Probiotic-directed modulation of gut microbiota is basal microbiome dependent. Gut Microbes:1–20. https://doi.org/10.1080/19490976.2020.1736974
Hu XL, Du H, Ren C, Xu Y (2016) Illuminating anaerobic microbial community and cooccurrence patterns across a quality gradient in Chinese Liquor Fermentation Pit Muds. Appl Environ Microbiol 82(8):2506–2515. https://doi.org/10.1128/AEM.03409-15
Hu XL, Du H, Xu Y (2015) Identification and quantification of the caproic acid-producing bacterium Clostridium kluyveri in the fermentation of pit mud used for Chinese strong-aroma type liquor production. Int J Food Microbiol 214:116–122. https://doi.org/10.1016/j.ijfoodmicro.2015.07.032
Imachi H, Sakai S, Kubota T, Miyazaki M, Saito Y, Takai K (2016) Sedimentibacter acidaminivorans sp nov, an anaerobic, amino-acid-utilizing bacterium isolated from marine subsurface sediment. Int J Syst Evol Microbiol 66:1293–1300. https://doi.org/10.1099/ijsem.0.000878
Jeong SE, Chun BH, Kim KH, Park D, Roh SW, Lee SH, Jeon CO (2018) Genomic and metatranscriptomic analyses of Weissella koreensis reveal its metabolic and fermentative features during kimchi fermentation. Food Microbiol 76:1–10. https://doi.org/10.1016/j.fm.2018.04.003
Langille MG, Zaneveld J, Caporaso JG, Mcdonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Vega Thurber RL, Knight R, Beiko RG, Huttenhower C (2013) Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol 31:814–821. https://doi.org/10.1038/nbt.2676
Li WQ (2010) Relationship between luzhou flavor liquor, lactic acid bacteria, lactic acid and ethyl lactate. Liquor Mak 37(03):90–93
Li XR, Ma EB, Yan LZ, Meng H, Du XW, Zhang SW, Quan ZX (2011) Bacterial and fungal diversity in the traditional Chinese liquor fermentation process. Int J Food Microbiol 146:31–37. https://doi.org/10.1016/j.ijfoodmicro.2011.01.030
Li Z, Bai Z, Wang D, Zhang W, Zhang M, Lin F, Gao L, Hui B, Zhang H (2014) Cultivable bacterial diversity and amylase production in three typical Daqus of Chinese spirits. Int J Food Sci Technol 49:776–786. https://doi.org/10.1111/ijfs.12365
Liang H, Luo Q, Zhang A, Wu Z, Zhang W (2016) Comparison of bacterial community in matured and degenerated pit mud from Chinese Luzhou-flavour liquor distillery in different regions. J Inst Brew 122:48–54. https://doi.org/10.1002/jib.296
Liu MK, Tang YM, Guo XJ, Zhao K, Penttinen P, Tian XH, Zhang XY, Ren DQ, Zhang XP (2020) Structural and functional changes in prokaryotic communities in artificial pit mud during Chinese Baijiu Production. msystems 5(2):e00829-19. https://doi.org/10.1128/mSystems.00829-19
Liu MK, Tang YM, Guo XJ, Zhao K, Tian XH, Liu Y, Yao WC, Deng B, Ren DQ, Zhang XP (2017) Deep sequencing reveals high bacterial diversity and phylogenetic novelty in pit mud from Luzhou Laojiao cellars for Chinese strong-flavor Baijiu. Food Res Int 102:68–76. https://doi.org/10.1016/j.foodres.2017.09.075
Liu QS, Li ZB, Yao Z (2003) Main factors influencing the growth and the metaboly of caproic acid bacteria. Liquor Mak Sci Technol:43–45. https://doi.org/10.13746/j.njkj.2003.04.012
Lozupone C, Knight R (2005) UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol 71:8228–8235. https://doi.org/10.1128/AEM.71.12.8228-8235.2005
Meng X, Wu Q, Wang L, Wang D, Chen L, Xu Y (2015) Improving flavor metabolism of Saccharomyces cerevisiae by mixed culture with Bacillus licheniformis for Chinese Maotai-flavor liquor making. J Ind Microbiol Biotechnol 42:1601–1608. https://doi.org/10.1007/s10295-015-1647-0
Price MN, Dehal PS, Arkin AP (2009) FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol 26:1641–1650. https://doi.org/10.1093/molbev/msp077
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glockner FO (2013) The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 41:D590–D596. https://doi.org/10.1093/nar/gks1219
Rainer BM, Thompson KG, Antonescu C, Florea L, Mongodin EF, Bui J, Fischer AH, Pasieka HB, Garza LA, Kang S (2020) Characterization and analysis of the skin microbiota in rosacea: a case–control study. Am J Clin Dermatol 21:139–147. https://doi.org/10.1007/s40257-019-00471-5
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C (2011) Metagenomic biomarker discovery and explanation. Genome Biol 12:R60. https://doi.org/10.1186/gb-2011-12-6-r60
Tao Y, Wang X, Li XZ, Wei N, Jin H, Xu ZC, Tang QL, Zhu XY (2017) The functional potential and active populations of the pit mud microbiome for the production of Chinese strong-flavour liquor. Microb Biotechnol 10:1603–1615. https://doi.org/10.1128/AEM.04070-13
Wang X, Du H, Xu Y (2017) Source tracking of prokaryotic communities in fermented grain of Chinese strong-flavor liquor. Int J Food Microbiol 244:27–35. https://doi.org/10.1016/j.ijfoodmicro.2016.12.018
Wickham H (2016) ggplot2: elegant graphics for data analysis. Springer, New York.
Xiang H, Sun-Waterhouse D, Waterhouse GIN, Cui C, Ruan Z (2019) Fermentation-enabled wellness foods: a fresh perspective. Food Sci Hum Wellness 8:203–243. https://doi.org/10.1016/j.fshw.2019.08.003
Xiang S, Chen J, Zhang Z (2015) The relationship between the fermentation temperature curve and solid state brewing between pre control conditions. Liquor Mak 42:28–31
Xie MW, Lv FX, Ma GX, Farooq A, Li HH, Du Y, Liu Y (2020) High throughput sequencing of the bacterial composition and dynamic succession in Daqu for Chinese sesame flavour liquor. J Inst Brew 126:98–104. https://doi.org/10.1002/jib.592
Xu Y, Sun B, Fan G, Teng C, Xiong K, Zhu Y, Li J, Li X (2017) The brewing process and microbial diversity of strong flavour Chinese spirits: a review. J Inst Brew 123:5–12. https://doi.org/10.1002/jib.404
Yang F, Chen L, Liu Y, Li J, Wang L, Chen J (2019) Identification of microorganisms producing lactic acid during solid-state fermentation of Maotai flavour liquor. J Inst Brew 125:171–177. https://doi.org/10.1002/jib.537
Yang F, Zhang Q, Liu Y, Li J, Wang L, Chen J (2020) Lactic acid biosynthesis pathways and important genes of Lactobacillus panis L7 isolated from the Chinese microbiome. Food Biosci 36:100627. https://doi.org/10.1016/j.fbio.2020.100627
Yang JG, Dou X, Han PJ, Bai FY, Zhou J, Zhang SY, Qin H, Ma YY (2017) Microbial diversity in Daqu during production of luzhou-flavored liquor. J Am Soc Brew Chem 75:136–144. https://doi.org/10.1094/ASBCJ-2017-2879-01
Yu Y, Lee C, Kim J, Hwang S (2005) Group-specific primer and probe sets to detect methanogenic communities using quantitative real-time polymerase chain reaction. Biotechnol Bioeng 89:670–679. https://doi.org/10.1002/bit.20347
Zhang R, Wu Q, Xu Y (2013) Aroma characteristics of Moutai-flavour liquor produced with Bacillus licheniformis by solid-state fermentation. Lett Appl Microbiol 57:11–18. https://doi.org/10.1111/lam.12087
Zhao XL, Zhou YK, Qiao Zong W, Li H, Li WH, Liu BH, Yang QS, Zhang LJ, Du ST, Xu ZH, Yang JZ, Tang Y (2013) A study on unique climate characteristics of Yibin of Sichuan in liquor production. Liquor Mak 40(05):5–10
Zheng J, Liang R, Zhang L, Wu C, Zhou R, Liao X (2013) Characterization of microbial communities in strong aromatic liquor fermentation pit muds of different ages assessed by combined DGGE and PLFA analyses. Food Res Int 54:660–666. https://doi.org/10.1016/j.foodres.2013.07.058
Zheng X, Han B (2016) Baijiu, Chinese liquor: history, classification and manufacture. J Ethn Foods 3:19–25. https://doi.org/10.1016/j.jef.2016.03.001
This work was financially supported by the Key Scientific and Technological Projects in the Field of Achievement Transformation in Xiangyang (2020AAT00469).
Ethics approval and consent to participate
Not applicable: no animals or human subjects were used in this project.
Consent for publication
All of the authors consent to the publication of this manuscript in Annals of Microbiology.
Author WL and SY were employed by the company Hubei Guxiangyang Liquor Industry 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.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Schematic diagram of pit mud sampling sites.
Rarefaction curves and Shannon diversity index curves of different samples.
Core operational taxonomic units (cOTUs) of bacteria in pit mud. The colors of the squares represent the mean relative abundance of each OTU.
About this article
Cite this article
Hou, Q., Wang, Y., Ni, H. et al. Deep sequencing reveals changes in prokaryotic taxonomy and functional diversity of pit muds in different distilleries of China. Ann Microbiol 72, 12 (2022). https://doi.org/10.1186/s13213-022-01671-x