Basic Sequencing Data and Alpha Diversity Analysis
Table 1
Basic Sequencing Data of the Pit Mud from F-S, G-Z, H-X and I-D
Sample | 16S | ITS |
Effective sequence | Shannon index | Chao 1 index | Coverage rate | Effective sequence | Shannon index | Chao 1 index | Coverage rate |
F-S-1 | 112576 | 5.2948 | 826 | 0.9976 | 128083 | 2.4599 | 86 | 0.9999 |
F-S-2 | 122754 | 5.3320 | 842 | 0.9978 | 120498 | 2.4471 | 91 | 0.9999 |
F-S-3 | 114750 | 5.3160 | 812 | 0.9979 | 120712 | 2.6153 | 92 | 0.9998 |
G-Z-1 | 127751 | 5.4721 | 785 | 0.9983 | 117690 | 3.1094 | 89 | 0.9998 |
G-Z-2 | 119701 | 5.4386 | 773 | 0.9982 | 117107 | 2.5827 | 87 | 0.9999 |
G-Z-3 | 119528 | 5.3970 | 762 | 0.9981 | 117972 | 2.2831 | 61 | 0.9999 |
H-X-1 | 117698 | 5.2033 | 815 | 0.9978 | 126473 | 3.6627 | 101 | 0.9998 |
H-X-2 | 116475 | 5.2284 | 754 | 0.9978 | 123248 | 3.5273 | 93 | 0.9999 |
H-X-3 | 113759 | 5.1491 | 768 | 0.9976 | 121489 | 3.4006 | 98 | 0.9999 |
I-D-1 | 125031 | 4.3747 | 815 | 0.9979 | 126599 | 3.1670 | 93 | 0.9999 |
I-D-2 | 120550 | 4.5248 | 800 | 0.9978 | 117278 | 3.2523 | 102 | 0.9999 |
I-D-3 | 120392 | 4.4931 | 804 | 0.9979 | 126217 | 3.1045 | 100 | 0.9999 |
Note: F-S: upper layer; G-Z: middle layer; H-X: lower layer; I-D: bottom layer. |
As can be seen from Table 1, effective sequences are obtained by the quality control of original sequences and the removal of chimeras. The average remaining effective sequences of pit mud samples in F-S, G-Z, H-X and I-D are 116693, 122327, 115977 and 121991, respectively. The coverage rate is more than 0.99, thus indicating that the sequencing depth is sufficient, the sequences in the samples are basically completely detected, and the results are true and reliable, which can be used for the subsequent analysis.
According to the Alpha diversity analysis of pit mud samples, Chao 1 index is mainly related to the abundance of samples, and the larger the index, the higher the abundance; Shannon index is mainly related to the diversity information of samples, which not only reflects the abundance of species, but also reflects the evenness of species. The larger the index, the higher the diversity. As can be seen from Table 1, in terms of the sequencing of 16S rDNA sequence, Chao 1 index can be ranked as F-S༞I-D༞H-X༞G-Z, and Shannon index can be ranked as G-Z༞F-S༞H-X༞I-D; In terms of the sequencing of ITS sequence, Chao 1 index can be ranked as I-D༞H-X༞F-S༞G-Z, and Shannon index can be ranked as H-X༞I-D༞G-Z༞F-S.
Venn Diagram Analysis
Figure 1 presents a Veen diagram based on OTUs (Operational Taxonomic Units). The overlapping part between graphs in the Venn diagram represents the number of common species, while the non-overlapping part represents the number of unique species of samples, which can clearly show the similarities and differences among samples. As can be seen from Fig. 3, the OTU number of the pit mud from F-S, G-Z, H-X and I-D are 691, 651, 629 and 662, respectively. The total number of common OTUs shared by four pit mud samples is 334, which indicates that 334 bacterial species exist simultaneously in the pit mud from F-S, G-Z, H-X and I-D; The number of unique OTUs in the pit mud from F-S, G-Z, H-X and I-D is 153, 120, 100 and 128, respectively.
Analysis of Bacterial Community Structure at the Phylum Level
As can be seen from Fig. 2, at the phylum level, with the relative abundance >0.5% as the threshold, there are 5 dominant bacterial phyla in the pit mud from F-S, including Firmicutes (52.8%), Bacteroidetes (29.6%), Synergistetes (6.4%), Chloroflexi (2.1%) and Spirochaetes (0.6%); There are 3 dominant bacterial phyla in the pit mud from G-Z, including Firmicutes (67.4%), Bacteroidetes (25.5%) and Synergistetes (0.9%); There are 5 dominant bacterial phyla in the pit mud from H-X, including Bacteroidetes (48.2%), Firmicutes (34.4%), Synergistetes (8.0%), Chloroflexi (2.7%) and Cloacimonetes (2.6%); There are 5 dominant bacterial phyla in the pit mud from I-D, including Firmicutes (66.7%), Bacteroidetes (14.7%), Synergistetes (10.3%), Kiritimatiellaeota (0.9%) and Chloroflexi (0.9%).
Analysis of Fungal Community Structure at the Phylum Level
As can be seen from Fig. 3, at the phylum level, with the relative abundance >0.5% as the threshold, the unclassified fungi (Unclassified) in the pit mud from F-S account for 0.06%, and there are 3 dominant fungal phyla, including Ascomycota (70.5%), Mucoromycota (28.9%) and Basidiomycota (0.6%); The unclassified fungi (Unclassified) in the pit mud from G-Z account for 0.2%, and there are 3 dominant fungal phyla, including Ascomycota (52.2%), Mucoromycota (44.6%) and Basidiomycota (3.0%). The unclassified fungi (Unclassified) in the pit mud from H-X account for 0.9%, and there are 3 dominant fungal phyla, including Ascomycota (53.2%), Mucoromycota (44.1%) and Basidiomycota (1.8%); The unclassified fungi (Unclassified) in the pit mud from I-D account for 2.0%, and there are 3 dominant fungal phyla, including Mucoromycota (58.5%), Ascomycota (37.4%) and Basidiomycota (2.0%).
Analysis of Bacterial Community Structure at the Genus Level
As can be seen from Fig. 4, at the genus level, with the relative abundance >1.0% as the threshold, the unclassified bacteria (Unclassified) in the pit mud from F-S account for 14.02%, and there are 11 dominant bacterial genera, including Hydrogenispora (23.67%), Petrimonas (12.75%), Caproiciproducens (7.32%), Proteiniphilum (6.61%), Ruminofilibacter (5.15%), Aminobacterium (4.06%), Lentimicrobium (3.27%), Christensenellaceae R-7 group (2.60%), Syner-01 (2.16%), Sedimentibacter (1.53%) and Syntrophomonas (1.49%); The unclassified bacteria (Unclassified) in the pit mud from G-Z account for 17.43%, and there are 11 dominant bacterial genera, including Caproiciproducens (28.00%), Lactobacillus (10.87%), Lentimicrobium (8.02%), Petrimonas (7.46%), Proteiniphilum (6.58%), Fermentimonas (2.13%), Hydrogenispora (1.64%), Herbinix (1.46%), Caldicoprobacter (1.22%), Sedimentibacter (1.12%), and Syntrophomonas (1.07%); The unclassified bacteria (Unclassified) in the pit mud from H-X account for 19.44%, and there are 9 dominant bacterial genera, including Proteiniphilum (16.10%), Blvii28_wastewater-sludge group (14.27%), Petrimonas (10.21%), Aminobacterium (7.78%), Hydrogenispora (7.63%), Caproiciproducens (3.92%), LNR_A2-18 (2.67%), Fermentimonas (1.76%) and Syntrophomonas (1.62%). The unclassified bacteria (Unclassified) in the pit mud from I-D account for 21.54%, and there are 8 dominant bacterial genera, including Hydrogenispora (36.92%), Petrimonas (10.33%), Aminobacterium (10.04%), Proteiniphilum (3.22%), Sedimentibacter (1.98%), Sporosarcina (1.68%), Syntrophomonas (1.56%) and Caproiciproducens (1.31%).
Analysis of Fungal Community Structure at the Genus Level
As can be seen from Fig. 5, at the genus level, with the relative abundance >1.0% as the threshold, the unclassified fungi (Unclassified) in the pit mud from F-S account for 1.5%, and there are 4 dominant fungal genera, including Thermomyces (42.7%), Rhizopus (28.8%), Aspergillus (19.6%) and Thermoascus (3.7%); The unclassified fungi (Unclassified) in the pit mud from G-Z account for 18.7%, and there are 7 dominant fungal genera, including Rhizopus (44.0%), Acremonium (8.5%), Cyphellophora (7.4%), Thermomyces (5.3%), Aspergillus (4.6%), Trichosporon (3.3%) and Thermoascus (1.4%); The unclassified fungi (Unclassified) in the pit mud from H-X account for 8.3%, and there are 7 dominant fungal genera, including Rhizopus (43.8%), Aspergillus (31.0%), Thermoascus (4.1%), Cladosporium (1.6%), Thermomyces (1.6%), Pseudeurotium (1.3%) and Penicillium (1.2%); The unclassified fungi (Unclassified) in the pit mud from I-D account for 18.5%, and there are 5 dominant fungal genera, including Rhizopus (57.7%), Aspergillus (10.2%), Thermoascus (3.3%), Penicillium (1.2%) and Hyphopichia (1.0%).
Composition and Content of Volatile Compounds in Pit Mud
The composition and content of volatile compounds in pit mud are listed in Table 2.
Table 2
Composition and Content of Volatile Compounds in Pit Mud
Compound Type | No. | Compound Name | Relative Content/% |
Pit mud from F-S | Pit mud from G-Z | Pit mud from H-X | Pit mud from I-D |
Esters | 1 | Ethyl 2-hydroxypropionate | 5.43 | 5.65 | 1.05 | 0.35 |
2 | Butyl caproate | 0.55 | 0.12 | 0.05 | 0.7 |
3 | Hexyl butyrate | 0.17 | - | - | 0.21 |
4 | Ethyl caprylate | 1.47 | 0.64 | 0.16 | 3.07 |
5 | Isoamyl caproate | 0.1 | - | - | 0.63 |
6 | Butyl lactate | 0.46 | - | - | - |
7 | Hexyl hexanoate | 0.3 | 0.27 | 0.27 | 10.68 |
8 | Ethyl caprate | 0.14 | 0.33 | 0.1 | 0.44 |
9 | Diethyl succinate | 0.34 | 0.33 | 0.08 | 0.02 |
10 | Ethyl methoxyacetate | 0.28 | 0.24 | 0.04 | 0.07 |
11 | Ethyl phenylacetate | 0.67 | 0.46 | 0.3 | 0.26 |
12 | Ethyl 3-Phenylpropionate | 0.75 | 0.33 | 0.29 | 0.3 |
13 | Ethyl methyl-4-pentenoate | 0.27 | 0.04 | - | - |
14 | Ethyl tetradecanoate | 0.16 | 0.18 | 0.14 | 0.16 |
15 | Ethyl pentadecanoate | 0.24 | 0.09 | 0.1 | 0.04 |
16 | Ethyl hexadecanoate | 0.91 | 0.77 | 0.49 | 0.78 |
17 | Ethyl 9-octadecenoic acid | 0.26 | 0.1 | - | 0.13 |
18 | Ethyl 9, 12-octadecadienoic acid | 0.34 | 0.21 | - | - |
19 | Ethylene acetate | 0.27 | - | - | - |
20 | Dibutyl phthalate | 0.23 | - | - | - |
21 | Ethyl caproate | 0.67 | 3.94 | 0.06 | 9.23 |
22 | Ethyl heptanoate | - | 0.14 | 0.06 | 0.65 |
23 | Butyl 2-hydroxypropanoate | - | 0.17 | - | - |
24 | Caproic acid 2,2-Dimethylhexanoic acid | - | 0.1 | - | - |
25 | Ethyl 13-methyl-tetradecanoate | - | 0.13 | - | - |
26 | Caproic acid 4-octyl ester | - | - | 0.04 | - |
27 | Allyl 2-ethylbutyrate | - | - | 0.2 | - |
28 | Ethyl tridecanoate | - | - | 0.11 | - |
29 | Amyl hexanoate | - | - | - | 0.66 |
30 | Propyl caprylate | - | - | - | 0.06 |
31 | Ethyl nonanoate | - | - | - | 0.12 |
32 | 2-ethylbutyric acid, 4-heptyl ester | - | - | - | 0.05 |
33 | Heptyl acetate | - | - | - | 1.04 |
34 | N-hexyl caprylate | - | - | - | 1.75 |
35 | Furfuryl acetate | - | - | - | 0.05 |
36 | Allyl 2-ethylbutyrate | - | - | - | 0.08 |
37 | Octyl pelargonate | - | - | - | 0.1 |
38 | Octyl octanoate | - | - | - | 0.07 |
39 | 9-Hexadecenoicacid, ethylester | - | - | - | 0.04 |
Alcohols | 1 | 1-hexanol | 2.56 | 1.83 | 0.41 | 0.95 |
2 | 1-octanol | 0.31 | 0.15 | - | 0.18 |
3 | 6-hendecanol | 0.09 | - | - | - |
4 | Benzyl alcohol | 0.15 | - | - | - |
5 | Phenethyl alcohol | 1.31 | 1.75 | 0.64 | 0.18 |
6 | Triethylene glycol | 0.09 | - | - | - |
7 | 1-butanol | - | 0.63 | - | - |
8 | 3-methyl-1-butanol | - | 0.37 | - | - |
9 | 1-decyl alcohol | - | - | - | 0.04 |
Acids | 1 | 2-methylpropionic acid | 1.45 | 0.69 | 0.54 | 0.06 |
2 | Acetic acid | 3.76 | 3.65 | 2.15 | 1.1 |
3 | Butyric acid | 16.83 | 12.11 | 12.53 | 1.77 |
4 | 3-methylbutyric acid | 0.06 | 1.9 | 2.45 | 0.28 |
5 | 3-methylpentanoic acid | 1.94 | - | - | - |
6 | Valeric acid | 2.78 | 3.06 | 11.28 | 0.93 |
7 | α-methyl phenylpropionic acid | 0.05 | - | - | - |
8 | 4-methylpentanoic acid | 0.14 | 0.18 | 0.19 | - |
9 | Caproic acid | 43.82 | 50.43 | 50.63 | 39.99 |
10 | Heptanoic acid | 1.04 | 1.52 | 4.25 | 4.72 |
11 | Caprylic acid | 1.86 | 2.22 | 3.57 | 14.36 |
12 | Nonanoic acid | 0.1 | 0.07 | 0.15 | 0.44 |
13 | Capric acid | 0.09 | 0.15 | 0.18 | 0.79 |
14 | Benzoic acid | 0.33 | 0.11 | 0.11 | - |
15 | Tetradeconic acid | 0.12 | - | - | - |
16 | Palmitic acid | 2.36 | 0.14 | 0.11 | - |
17 | Propionic acid | - | 0.06 | 0.58 | 0.07 |
18 | Phenylacetic acid | - | 0.07 | 0.08 | - |
19 | 2-methylbutyric acid | - | | 0.01 | - |
20 | 2-methylpentanoic acid | - | - | 0.04 | - |
21 | 2-phenethyl hexanoic acid | - | - | - | 0.19 |
22 | Cis10-heptadecenoic acid | - | - | - | 0.14 |
Others | 1 | Phenol | 0.19 | 0.32 | 0.3 | 0.32 |
2 | P-methylphenol | 1.49 | 2.17 | 4.77 | 0.82 |
3 | Cresol | - | 0.05 | - | - |
4 | 3-methylphenol | - | - | - | 0.03 |
5 | 2,4-di-tert-butylphenol | - | - | - | 0.16 |
6 | Tetramethylpyrazine | - | 0.38 | 0.14 | - |
7 | D-limonene | - | 0.16 | - | - |
A total of 78 volatile compounds are detected from the pit mud used for manufacturing Taorong-type Baijiu, including 39 esters, 9 alcohols, 22 acids and 7 others. Esters and acids are two dominant components in pit mud. There are significant differences in esters and acids, as well as their content in different layers of pit mud; while, the differences in alcohols are not significant. The total content of volatile compounds in pit mud shows an upward-downward-upward trend with the pit depth. There are 46 kinds of volatile compounds in the pit mud from F-S, 45 kinds from G-Z, 39 kinds from H-X and 49 kinds from I-D. Besides, there are also differences in the relative contents of various components in the pit mud from F-S, G-Z, H-X and I-D. Ester compounds are the most volatile compounds with the highest content and variety in pit mud, 39 of which are the main contributors to the aroma of Taorong-type Baijiu, and there are 21, 20, 17 and 28 kinds in the pit mud from F-S, G-Z, H-X and I-D, respectively. Ethyl caproate ranks first with respect to the content of esters, with its content in the pit mud from F-S, G-Z, H-X and I-D being 0.67%, 3.94%, 0.06% and 9.23% respectively. It is mainly generated under the synergistic action of various bacteria and enzymes[11]. Ethanol and acetic acid combine to form butyric acid, followed by the synthesis of caproic acid through esterase. Subsequently, caproic acid is synthesized through ethanol[27]. Ethyl caproate is considered as the key component contributing to the flavor and quality of Baijiu[28]. Ethyl caprylate and ethyl heptanoate ranks second and third with respect to the content of esters. It can be found that the content of ethyl esters in pit mud is the highest, and there are 13, 15, 12 and 14 kinds of ethyl esters in the pit mud from F-S, G-Z, H-X and I-D, respectively, which are the main esters in Taorong-type Baijiu. There is no significant difference in ethanol compounds in pit mud, and 9 kinds of ethanol compounds are detected, including 6, 5, 2 and 4 kinds in the pit mud from F-S, G-Z, H-X and I-D, respectively. 1-Hexanol is the ethanol compound with the highest content in pit mud, with its relative contents in the pit mud from F-S, G-Z, H-X and I-D being 2.56%, 1.83%, 0.41% and 0.95%, respectively. There are abundant acid compounds in pit mud, and there are significant differences in different layers of pit mud. A total of 22 kinds of acid compounds are detected, including 16, 25, 17 and 13 kinds in the pit mud from F-S, G-Z, H-X and I-D, respectively. Among them, caproic acid, butyric acid, acetic acid, valeric acid, caprylic acid and heptanoic acid have a higher content. The content of caproic acid is extremely high, with its relative contents in the pit mud from F-S, G-Z, H-X and I-D being 43.82%, 50.43%, 50.63% and 39.99%, respectively. Under the coupling action of Caproiciproducens and Methanogen, acetic acid is produced by ethanol oxidation, and then ethanol reacts with butyric acid to produce caproic acid[29]. Caproic acid and ethyl caproate produced by pit mud fermentation are the main aroma components of Taorong-type Baijiu[11]. A large amount of caproic acid is produced during fermentation, and then reacts with ethanol to produce ethyl caproate, the main aroma component of Baijiu[30]. Acetic acid, butyric acid, heptanoic acid and caprylic acid have a second higher content, and they are the main organic acid components of Taorong-type Baijiu.
Correlation between Microbes and Main Volatile Compounds in Pit Mud
Note
Z-A: ethyl 2-hydroxypropionate; Z-B: butyl caproate; Z-C: ethyl caprylate; Z-D: hexyl hexanoate; Z-E: ethyl caprate; Z-F: ethyl phenylacetate; Z-G: ethyl 3-phenylpropionate; Z-H: ethyl hexadecanoate; Z-I: ethyl caproate; Z-J: ethyl heptanoate; Z-K: heptyl acetate; Z-L: N-hexyl caprylate; C-A: 1-hexanol; C-B: 1-octanol; C-C: 6-hendecanol; C-D: benzyl alcohol; C-E: phenethyl alcohol; C-F: triethylene glycol; C-G: 1-butanol C-H: 3-methyl-1-butanol C-I: 1-decyl alcohol; S-A: 2-methylpropionic acid; S-B: acetic acid; S-C: butyric acid; S-D: 3-methylbutyric acid; S-E: 3-methylpentanoic acid; S-F: valeric acid; S-G: caproic acid; S-H: heptanoic acid; S-I: caprylic acid; S-J: nonanoic acid; S-K: capric acid; S-L: palmitic acid; Q-A: phenol; Q-B: p-methylphenol; Q-C: cresol; Q-D: 3-methylphenol; Q-E: 2,4-di-tert-butylphenol
The correlation analysis between the main volatile compounds and bacteria and fungi in pit mud was conducted respectively, and the correlation heat map was obtained. As can be seen in Fig. 6, the bacterial genera that are closely correlated with the main volatile compounds in pit mud include Hydrogenispora, Aminobacterium, Lentimicrobium, Sedimentibacter, Ruminofilibacter, Christensenellaceae_R-7_group and Syner-01. As can be seen from Fig. 7, the fungal genera that are closely correlated with the main volatile compounds in pit mud include Rhizopus, Thermomyces, Monascus and Penicillium. Ethyl caproate is the main aroma substance in Baijiu, which has the highest correlation with Sedimentibacter and Monascus, followed by Hydrogenispora and Rhizopus. Sedimentibacter can synthesize caproic acid, butyric acid, acetic acid, hexanol, ethanol and butanol with carbon source and protein as substrates, and can generate ethyl caproate. Acetic acid is positively correlated with Syner-01, Ruminofilibacter, Lentimicrobium, Caproiciproducens and Thermomyces, and negatively correlated with Aminobacterium, Monascus and Penicillium. Hexanol, octanol, 6-hendecanol and benzyl alcohol are positively correlated with Syner-01, Christensenellaceae_R-7_group, Ruminofilibacter and Thermomyces. Therefore, the microbial community structure in pit mud has a certain influence on the flavor and quality of Baijiu.