Zooplankton-associated microbiomes play an important role for host health, and contribute to ecosystem processes such as nutrient cycling. Yet, few studies have assessed how environmental gradients and biotic interactions, including parasitism and diet, may shape the microbiome composition of wild zooplankton. Here, we analyzed the microbiomes of water fleas from the Daphnia longispina species complex using 16S rRNA gene sequencing and a long-term field dataset spanning six sampling events over 13 years. Sampling coincided with outbreaks of the virulent eukaryotic gut parasite Caullerya mesnili. Additionally, we explored how microbiome structure varied in relation to water parameters, phytoplankton density (i.e., Daphnia diet), and zooplankton density and community structure. Daphnia microbiomes displayed strong temporal variation and comparatively small differences based on host infection status. Microbiome beta diversity correlated with phytoplankton density but not with its community composition, including green algae, protists, and cyanobacteria. Environmental conditions, including temperature, dissolved oxygen, and cyanobacterial abundance—previously found to drive Caullerya epidemics—were also associated with distinct microbiome structures. Importantly, microbiome beta diversity co-varied with infection prevalence, suggesting a link between microbiome shifts, epidemic size, and environmental conditions driving large epidemics. Dominant bacterial taxa correlated with Daphnia density, whereas the phylogenetic composition of rare taxa was associated with total zooplankton density. These findings demonstrate the dynamic nature of Daphnia microbiomes and suggest potential mechanisms by which they may mediate disease dynamics, particularly through associations with diet quantity, temperature, and host population density.
