Multi-omics reveal microbial functional traits and antifungal metabolites associated with lower Pseudogymnoascus destructans loads in bat cave soils.
White-nose syndrome, caused by Pseudogymnoascus destructans (Pd), is a major fungal disease threatening hibernating bats. Cave soils can serve as environmental reservoirs for Pd, yet the microbial and biochemical mechanisms underlying naturally low Pd burdens in some cave environments remain poorly understood. Here, we integrated soil microbiome profiling, metagenomics, metabolomics, multi-omics network analysis, and in vitro validation to investigate the ecological and functional basis of differential Pd loads in hibernating bat caves in Northeast China. The three caves shared cold, humid, and weakly acidic microenvironments, but differed significantly in electrical conductivity, soil water content, nutrient availability, and extracellular enzyme activities. Soil microbial communities showed significant inter-cave variation in composition, diversity, and niche breadth, with stochastic processes contributing substantially to community assembly. Environmental variables, particularly pH and Pd load, were important predictors of microbial community structure. Functional analyses revealed that the low-Pd Gezi Cave was enriched in genes associated with organic carbon degradation, nitrogen input and retention, and secondary metabolism. Metabolomic profiling further identified cave-specific metabolite signatures, among which Biochanin A, 4-Hydroxybenzaldehyde, Vanillin, and Arachidonic acid were negatively correlated with Pd loads. Integrated pathway and network analyses showed that differential genes and metabolites jointly mapped to secondary metabolite biosynthesis, aminobenzoate degradation, and flavonoid degradation pathways, forming a microbe-metabolite-functional gene coupling network involving key taxa such as Rhodococcus, Pseudorhodoplanes, and Rhodoplanes. In vitro assays confirmed that 4-Hydroxybenzaldehyde, Coumarin, and Vanillin inhibited Pd growth. Structural equation modelling further indicated that environmental heterogeneity was associated with variation in Pd loads through microbial functional attributes and metabolite profiles. These findings suggest that naturally low-Pd cave soils are associated with coordinated environmental filtering, microbial functional specialization, and antifungal metabolite production, providing mechanistic insight into microbial and biochemical constraints on Pd persistence in cave reservoirs.