Stochastic Processes Drive the Assembly and Metabolite Profiles of Keystone Taxa during Chinese Strong-Flavor Baijiu Fermentation

Multispecies communities participate in the fermentation of Chinese strong-flavor Baijiu (CSFB), and the metabolic activity of the dominant and keystone taxa is key to the flavor quality of the final product. However, their roles in metabolic function and assembly processes are still not fully understood. Here, we identified the variations in the metabolic profiles of dominant and keystone taxa and characterized their community assembly using 16S rRNA and internal transcribed spacer (ITS) gene amplicon and metatranscriptome sequencing. We demonstrate that CSFB fermentations with distinct metabolic profiles display distinct microbial community compositions and microbial network complexities and stabilities. We then identified the dominant taxa (Limosilactobacillus fermentum, Kazachstania africana, Saccharomyces cerevisiae, and Pichia kudriavzevii) and the keystone ecological cluster (module 0, affiliated mainly with Thermoascus aurantiacus, Weissella confusa, and Aspergillus amstelodami) that cause changes in metabolic profiles. Moreover, we highlight that the alpha diversity of keystone taxa contributes to changes in metabolic profiles, whereas dominant taxa exert their influence on metabolic profiles by virtue of their relative abundance. Additionally, our results based on the normalized stochasticity ratio (NST) index and the neutral model revealed that stochastic and deterministic processes together shaped CSFB microbial community assemblies. Stochasticity and environmental selection structure the keystone and dominant taxa differently. This study provides new insights into understanding the relationships between microbial communities and their metabolic functions. IMPORTANCE From an ecological perspective, keystone taxa in microbial networks with high connectivity have crucial roles in community assembly and function. We used CSFB fermentation as a model system to study the ecological functions of dominant and keystone taxa at the metabolic level. We show that both dominant taxa (e.g., those taxa that have the highest relative abundances) and keystone taxa (e.g., those taxa with the most cooccurrences) affected the resulting flavor profiles. Moreover, our findings established that stochastic processes were dominant in shaping the communities of keystone taxa during CSFB fermentation. This result is striking as it suggests that although the controlled conditions in the fermentor can determine the dominant taxa, the uncontrolled rare keystone taxa in the microbial community can alter the resulting flavor profiles. This important insight is vital for the development of potential manipulation strategies to improve the quality of CSFB through the regulation of keystone species.