Monitoring longitudinal antimicrobial resistance trends of Staphylococcus aureus strains worldwide over the past 100 years to decipher its evolution and transmission.

Staphylococcus aureus inhabits diverse habitats including food waste and wastewater treatment plants. Cases of S. aureus-induced infection are commonly reported worldwide. The emergence of antimicrobial resistance (AMR) of S. aureus is a growing public health threat worldwide. Here, we longitudinally monitored global trends in antibiotic resistance genes (ARGs) of 586 S. aureus strains, isolated between 1884 and 2022. The ARGs in S. aureus exhibited a significant increase over time (P < 0.0001). Mobile genetic elements play a crucial role in the transfer of ARGs in S. aureus strains. The structural equation model results revealed a significant correlation between the human development index and rising antibiotic consumption, which subsequently leads to an indirect escalation of AMR in S. aureus strains. Lastly, a machine learning algorithm successfully predicted the AMR risk of global terrestrial S. aureus with over 70% accuracy. Overall, these findings provided valuable insights for managing AMR in S. aureus.

Horizontal gene transfer in activated sludge enhances microbial antimicrobial resistance and virulence.

Activated sludge (AS) plays a vital role in removing organic pollutants and nutrients from wastewater. However, the risks posed by horizontal gene transfer (HGT) between bacteria in AS are still unclear. Here, a total of 478 high-quality non-redundant metagenome-assembled genomes (MAGs) were obtained. >50 % and 5 % of MAGs were involved in at least one HGT and recent HGT, respectively. Most of the transfers (82.4 %) of antimicrobial resistance genes (ARGs) occurred among the classes of Alphaproteobacteria and Gammaproteobacteria. The bacteria involved in the transfers of virulence factor genes (VFGs) mainly include Alphaproteobacteria (42.3 %), Bacteroidia (19.2 %), and Gammaproteobacteria (11.5 %). Moreover, the number of ARGs and VFGs in the classes of Alphaproteobacteria and Gammaproteobacteria was higher than that in other bacteria (P < 0.001). Mobile genetic elements were important contributors to ARGs and VFGs in AS bacteria. These results have implications for the management of antimicrobial resistance and virulence in activated sludge microorganisms.

Distinct increase in antimicrobial resistance genes among Vibrio parahaemolyticus in recent decades worldwide.

Vibrio parahaemolyticus is a common pathogen, and has emerged with multiple antimicrobial resistance (AMR). However, few studies have conducted large-scale investigations of AMR and virulence trends of V. parahaemolyticus worldwide. This study longitudinally monitored antibiotic resistance genes (ARGs) and virulence factor genes (VFGs) trends of 1540 V. parahaemolyticus isolates isolated from 1951 to 2021. The number of ARGs in V. parahaemolyticus isolates distinctly increased over the years (P = 5.9e-10), while the number of VFGs decreased significantly (P < 2.2e-16). However, the number of VFGs of isolates isolated from humans has not changed significantly over the years (R = 0.013, P = 0.74), suggesting that the pathogenic risk to humans has not been reduced. Besides, mobile genetic elements are important contributors to ARGs in V. parahaemolyticus (R = 0.34, P < 2.2e-16), but have no promoting effect on VFGs (P = 0.50). The structural equation model illustrated that the human development index promoted the consumption of antibiotics, thereby indirectly promoting an increase in the AMR of the V. parahaemolyticus isolates. Finally, the random forest was performed to predict the ARG and VFG risks of global terrestrial V. parahaemolyticus isolates, and successfully map these threats with over 80% accuracy. This study aimed to evaluate the global risks posed by AMR and virulence, which helps to develop methods specifically targeting V. parahaemolyticus to mitigate these threats.

Monitoring Longitudinal Trends and Assessment of the Health Risk of Shigella flexneri Antimicrobial Resistance.

Shigella flexneri infection is the main cause of diarrhea in humans worldwide. The emergence of antimicrobial resistance (AMR) of S. flexneri is a growing public health threat worldwide, while large-scale studies monitoring the longitudinal AMR trends of isolates remain scarce. Here, the AMR gene (ARG) profiles of 717 S. flexneri isolates from 1920 to 2020 worldwide were determined. The results showed that the average number of ARGs in isolates has increased significantly, from 19.2 ± 2.4 before 1970 to 29.6 ± 5.3 after 2010. In addition, mobile genetic elements were important contributors to ARGs in S. flexneri isolates. The results of the structural equation model showed that the human development index drove the consumption of antibiotics and indirectly promoted the antibiotic resistance. Finally, a machine learning algorithm was used to predict the antibiotic resistance risk of global terrestrial S. flexneri isolates and successfully map the antibiotic resistance threats in global land habitats with over 80% accuracy. Collectively, this study monitored the longitudinal AMR trends, quantitatively surveilled the health risk of S. flexneri AMR, and provided a theoretical basis for mitigating the threat of antibiotic resistance.

Metabolite-Based Mutualistic Interaction between Two Novel Clostridial Species from Pit Mud Enhances Butyrate and Caproate Production.

Pit mud microbial consortia play crucial roles in the formation of Chinese strong-flavor baijiu's key flavor-active compounds, especially butyric and caproic acids. Clostridia, one of the abundant bacterial groups in pit mud, were recognized as important butyric and caproic acid producers. Research on the interactions of the pit mud microbial community mainly depends on correlation analysis at present. Interaction between Clostridium and other microorganisms and its involvement in short/medium-chain fatty acid (S/MCFA) metabolism are still unclear. We previously found coculture of two clostridial strains isolated from pit mud, Clostridium fermenticellae JN500901 (C.901) and Novisyntrophococcus fermenticellae JN500902 (N.902), could enhance S/MCFA accumulation. Here, we investigated their underlying interaction mechanism through the combined analysis of phenotype, genome, and transcriptome. Compared to monocultures, coculture of C.901 and N.902 obviously promoted their growth, including shortening the growth lag phase and increasing biomass, and the accumulation of butyric acid and caproic acid. The slight effects of inoculation ratio and continuous passage on the growth and metabolism of coculture indicated the relative stability of their interaction. Transwell coculture and transcriptome analysis showed the interaction between C.901 and N.902 was accomplished by metabolite exchange, i.e., formic acid produced by C.901 activated the Wood-Ljungdahl pathway of N.902, thereby enhancing its production of acetic acid, which was further converted to butyric acid and caproic acid by C.901 through reverse ß-oxidation. This work demonstrates the potential roles of mutually beneficial interspecies interactions in the accumulation of key flavor compounds in pit mud. IMPORTANCE Microbial interactions played crucial roles in influencing the assembly, stability, and function of the microbial community. The metabolites of pit mud microbiota are the key to flavor formation of Chinese strong-flavor baijiu. So far, researches on the interactions of the pit mud microbial community have been mainly based on the correlation analysis of sequencing data, and more work needs to be performed to unveil the complicated interaction patterns. Here, we identified a material exchange-based mutualistic interaction system involving two fatty acid-producing clostridial strains (Clostridium fermenticellae JN500901 and Novisyntrophococcus fermenticellae JN500902) isolated from pit mud and systematically elucidated their interaction mechanism for promoting the production of butyric acid and caproic acid, the key flavor-active compounds of baijiu. Our findings provide a new perspective for understanding the complicated interactions of pit mud microorganisms.

Spatial heterogeneity of the microbiome and metabolome profiles of high-temperature Daqu in the same workshop.

High-temperature Daqu, usually used as a fermentation starter for sauce-flavor Baijiu production, plays an essential role in the yield and flavor quality of Baijiu. The environmental heterogeneity of different locations in the workshop during fermentation led to the final production of Daqu with three different types (i.e., white, yellow, and black Daqu). How to use these three types of Daqu in Baijiu production mainly depends on the workers' experience so far. Here, we aimed to reveal the potential functions of different types of Daqu by comparing enzyme activity, volatile metabolites, and microbiota characteristics. White_Qu exhibited the highest liquefaction and saccharification enzyme activities, while the highest neutral proteinase and cellulase enzyme activities were detected in black_Qu. The total volatile content of yellow_Qu and black_Qu was roughly double that of white_Qu, and multivariate analysis revealed distinct volatile dissimilarities across different types of Daqu. Significant differences in bacterial and fungal community structures, assembly patterns, and potential functional profiles were discovered among different types of Daqu. At the genus level, Oceanobacillus and Thermomyces dominated the white_Qu microbiota, and the abundant microbes in yellow_Qu and black_Qu were scattered in Kroppenstedtia and Thermoascus. Bacterial and fungal communities were dominated by deterministic and stochastic assembly processes, respectively, suggesting that bacteria may be more affected by abiotic environmental factors and species interaction than fungi. Co-occurrence network analysis showed positive correlations characterized Daqu microbial networks, and network topological features indicated stronger interactions between bacterial taxa compared with fungal community. The Spearman correlation analysis revealed that four bacterial genera (Kroppenstedtia, Virgibacillus, Scopulibacillus, and Staphylococcus) and two fungal genera (Thermoascus and Aspergillus) exhibited positive correlations with almost all of the abundant volatiles. This work reveals that spatially varying environments lead to the microbiome and metabolome heterogeneity of high-temperature Daqu in the same workshop.

Comparative Genomics Unveils the Habitat Adaptation and Metabolic Profiles of Clostridium in an Artificial Ecosystem for Liquor Production.

Clostridium inhabiting pit mud (PM) is one of the important bacterial populations for synthesizing flavor compounds of Chinese strong-flavor baijiu. The long-term cereal fermentation with sorghum as the main raw material creates an environment rich in starch, ethanol, and organic acids (mainly lactic acid). However, the genetic factors underpinning Clostridium's adaptation to PM remain poorly understood. Here, we performed comparative genomic analysis between 30 pit mud-associated (PMA) and 100 non-pit mud-associated (NPMA) Clostridium strains. Comparison analysis of the enrichment of KEGG pathways between PMA and NPMA Clostridium strains showed two-component system, flagellar assembly, and bacterial chemotaxis pathways related to environmental adaptation were enriched in PMA strains. The number of genes encoding alcohol dehydrogenase and l-lactate dehydrogenase in PMA Clostridium strains was significantly higher than that in NPMA, which is helpful for them to adapt to the ethanol- and lactic acid-rich environment. The analysis of carbohydrate-active enzymes demonstrated that glycoside hydrolases (GHs) was the most abundant family in all Clostridium strains, and genes encoding GH4 and GH13, involved in starch and sucrose metabolism, were enriched in PMA Clostridium. Horizontal gene transfer analysis revealed that multiple genes encoding the enzymes involved in carbohydrate and amino acid metabolism were transferred from Bacillus to Clostridium in pit mud. Most of the PMA Clostridium strains had good potential for butyric acid synthesis from ethanol, lactic acid, and starch. Collectively, this study furthers our understanding of the habitat adaptation and metabolic potential of PMA Clostridium strains. IMPORTANCE Pit mud is a typical artificial ecosystem for Chinese liquor production. Clostridium inhabiting pit mud plays essential roles in the flavor formation of strong-flavor baijiu. The relative abundance of Clostridium increased with pit mud quality, further influencing the quality of baijiu. So far, the ecological adaptation of Clostridium to a pit mud-associated lifestyle is largely unknown. Here, comparative genomic analysis of pit mud-associated (PMA) and non-pit mud-associated (NPMA) Clostridium strains was performed. We found genes related to the metabolism of starch, ethanol, and lactic acid were enriched in PMA Clostridium strains, which facilitated their adaptation to the unique brewing environment. In addition, horizontal gene transfer contributed to the adaptation of Clostridium to pit mud. Our findings provide genetic insights on PMA Clostridium strains' ecological adaptation and metabolic characteristics.

Novisyntrophococcus fermenticellae gen. nov., sp. nov., isolated from an anaerobic fermentation cellar of Chinese strong-flavour baijiu.

A Gram-stain-negative, coccus-shaped, obligately anaerobic, non-motile and non-spore-forming bacterium, designated strain JN500902T, was isolated from the mud in a fermentation cellar used continuously over 30 years for Chinese strong-flavour baijiu production. Colonies were white, circular, convex and smooth-edged. Growth was observed at 20-40 °C (optimum, 37 °C), at pH 5.0-10 (optimum, pH 7.5), with 0-2 % (w/v) NaCl and with 0-4 % (v/v) ethanol. The Biolog assay demonstrated positive reactions of strain JN500902T in the metabolism of l-fucose and pyruvate. The predominant cellular fatty acids (>10 %) consisted of C16 : 0 and C14 : 0. The major end metabolites of strain JN500902T were acetic acid and ethanol when incubated anaerobically in liquid reinforced clostridial medium. Acetate was the major organic acid end product. The complete genome size of strain JN500902T was 3 420 321 bp with 3327 identified genes. The G+C content was 43.5 mol%. Phylogenetic analysis based on 16S rRNA gene sequences affiliated strain JN500902T with the family Lachnospiraceae, having low sequence similarity (92.8 %) to the nearest type strain, Syntrophococcus sucromutans DSM 3224T and forming a clearly distinct branch. Core genome phylogenetic analysis of the isolate and 134 strains belonging to the family Lachnospiraceae also revealed that strain JN500902T was well-separated from other genera of this family as a monophyletic clade. The average nucleotide identity and amino acid identity values between strain JN500902T and 134 Lachnospiraceae strains were less than 74 and 65 %, respectively. Considering its polyphasic characteristics, strain JN500902T represents a novel genus and species within the family Lachnospiraceae, for which the name Novisyntrophococcus fermenticellae gen. nov., sp. nov. is proposed. The type strain is JN500902T (=CICC 24502T=JCM 33939T).

Deciphering the succession patterns of bacterial community and their correlations with environmental factors and flavor compounds during the fermentation of Zhejiang rosy vinegar.

Zhejiang Rosy Vinegar (ZRV) is a traditional condiment in Southeast China, produced using semi-solid-state fermentation techniques under an open environment, yet little is known about the functional microbiota involved in the flavor formation of ZRV. In this study, 43 kinds of volatile flavor substances were identified by HS-SPME/GC-MS, mainly including ethyl acetate (relative content at the end of fermentation: 1104.1 mg/L), phenylethyl alcohol (417.6 mg/L) and acetoin (605.2 mg/L). The most abundant organic acid was acetic acid (59.6 g/L), which kept rising during the fermentation, followed by lactic acid (7.0 g/L), which showed a continuously downward trend. Amplicon sequencing analysis revealed that the richness and diversity of bacterial community were the highest at the beginning and then maintained decreasing during the fermentation. The predominant bacteria were scattered in Acetobacter (average relative abundance: 63.7%) and Lactobacillus (19.8%). Both sequencing and culture-dependent analysis showed Lactobacillus dominated the early stage (day 10 to 30), and Acetobacter kept highly abundant from day 40 to the end. Spearman correlation analysis displayed that the potential major groups involved in the formation of flavor compounds were Acetobacter and Lactobacillus, which were also showed strong relationships with other bacteria through co-occurrence network analysis (edges attached to Acetobacter: 61.7%; Lactobacillus: 14.0%). Moreover, structural equation model showed that the contents of ethanol, titratable acid and reducing sugar were the major environmental factors playing essential roles in influencing the succession of bacterial community and their metabolism during the fermentation. Overall, these findings illuminated the dynamic profiles of bacterial community and flavor compounds and the potential functional microbes, which were expected to help us understand the formation of flavor substances in ZRV.