Co-Occurrence Patterns of Bacterial Communities and Resistance Genes: A Comprehensive Multi-Pen Fecal Microbiome and Resistome Study in Dairy Farms

Antimicrobial resistance (AMR) poses a critical public health threat, with rising multidrug resistance cases compromising treatment effectiveness. Knowledge about the resistome in dairy production systems remains limited, particularly regarding lactating cows. This study investigated the microbiome and resistome across the hospital, fresh, and mid-lactation pens on 18 conventional dairy farms in California and Ohio using shotgun metagenomic sequencing of pooled fecal samples. Pooled fecal pat samples were collected as part of a larger field study using a quasi-experimental design that assigned farms to the training intervention group (six per state) or the control group (three per state). For the training intervention group, farm worker(s), identified as having the task of diagnosing and treating adult cows on the farm, participated in a training program on antimicrobial stewardship practices. Pooled fecal samples (<i>n</i> = 7) were collected at enrollment and 3 months after the intervention was completed on each participating farm (<i>n</i> = 18). A total of 10,221 bacterial species and 345 AMR genes conferring resistance to 22 antimicrobial classes were identified. The hospital pen exhibited a higher AMR gene diversity compared to fresh and mid-lactation pens (<i>p</i> < 0.05). Several AMR genes showed bimodal distribution, suggesting complex transmission mechanisms. Network analysis revealed distinct gene correlation profiles across pens, with the hospital pen showing fewer gene interactions. Our findings suggest that farm-level antimicrobial drug use may not be the sole or primary driver of resistome composition in pooled fecal samples from dairy cattle, highlighting the need to investigate other factors influencing AMR dynamics in livestock systems.