Characterization of the microbiota and resistome of bulk tank milk samples from Prince Edward Island dairy farms

Abstract

Bulk tank milk (BTM) is regularly used for surveillance on dairy farms for disease conditions such as mastitis and Johne's disease. In this study, we used 16S rRNA sequencing and bait-capture enrichment to characterize the microbiota and resistome of BTM, and investigate potential differences between the cream or pellet fractions. A total of 12 BTM samples were taken from 12 Prince Edward Island dairy farms, in Atlantic Canada, in duplicates. The DNA was analyzed by high-throughput sequencing of the 16S rRNA gene and a suite of antimicrobial resistance (AMR) genes. Target-capture enrichment of AMR genes was conducted before shotgun sequencing. The bioinformatics pipelines QIIME 2 and AMR++ were used for microbiota and resistome analysis, respectively. Differences between microbiotae were evaluated qualitatively with nonmetric multidimensional scaling and quantitatively with permutational ANOVA of UniFrac distances. A total of 47 phyla were present across the BTM samples. Firmicutes, Bacteroidetes, Proteobacteria, and Actinobacteria were the most abundant phyla. At the genus level, Corynebacterium, Acinetobacter, Lactobacillus, and Turicibacter were the most abundant. There was no significant difference in the Faith's phylogenetic diversity between the cream and pellet fraction. Faith's phylogenetic diversity differed marginally by stall type. There were 10,217 hits across 80 unique AMR genes, with tetracycline resistance being the most common class.