Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Vet. Sci., 23 July 2025

Sec. Veterinary Infectious Diseases

Volume 12 - 2025 | https://doi.org/10.3389/fvets.2025.1649876

Diversity and antibiotic resistance of cultivable bacteria in bulk tank milk from dairy farms in Shandong Province, China

Yijian Qi&#x;Yijian Qi1Zhiyuan Lu&#x;Zhiyuan Lu2Ziru MengZiru Meng3Xiaozhou WangXiaozhou Wang1Huahua ChenHuahua Chen1Muzi LiMuzi Li1Chaonan QuChaonan Qu1Pu Zhang
Pu Zhang4*Yongxia Liu
Yongxia Liu2*Jianzhu Liu
Jianzhu Liu1*
  • 1College of Veterinary Medicine, Shandong Agricultural University, Tai’an, Shandong, China
  • 2Shandong Provincial Key Laboratory of Zoonoses, Shandong Agricultural University, Tai’an, China
  • 3College of Life Sciences, Shandong Agricultural University, Tai'an, China
  • 4The Affiliated Tai’an City Central Hospital of Qingdao University, Tai’an, China

Introduction: This study systematically analyzed bacterial diversity and antimicrobial resistance (AMR) profiles in bulk tank milk from five dairy farms (n = 30) in Shandong Province, China, to assess public health risks associated with microbial contamination and provide critical data for regional quality control and AMR risk assessment in dairy production systems.

Methods: Total bacterial counts were quantified, revealing significant inter-farm variation (P < 0.05) with a range of 3.94–6.68 log CFU/mL. Among 129 bacterial isolates, genus-level dominance and species prevalence were identified. Antimicrobial susceptibility testing (AST) against 10 agents was performed using integrated resistance criteria combining Clinical and Laboratory Standards Institute (CLSI) standards and epidemiological cutoff values (ECOFFs). Nine resistance genes targeting seven antibiotic classes were detected via PCR.

Results: The highest resistance rate was observed for sulfadiazine (53.2%) and the lowest for levofloxacin (6.0%). Multidrug resistance was detected in 23% (20/87) of isolates, with 14 strains meeting ECOFFs-based resistance criteria. PCR analysis showed sul1 (70.5%) and ant(4′)-Ia (54.3%) as the most prevalent resistance genes, while mcr-1, lnu (B), and blaNDM-1 were absent in all isolates. Regional resistance variations correlated significantly with farm management practices.

Discussion: These findings underscore the impact of historical antibiotic use on AMR dissemination. Enhanced AMR surveillance in raw milk, improved antibiotic stewardship, and targeted interventions are crucial to mitigate public health risks from microbial contamination and horizontal gene transfer of resistance determinants.

1 Introduction

Raw milk acts as a critical vector for microbial contamination due to its physicochemical properties, with bacterial loads directly compromising dairy safety (1). The presence of spoilage and pathogenic bacteria in bulk tank milk poses significant risks to public health (2). Cattle ingesting contaminated forage introduce bacteria into their systems, where they proliferate and disseminate into the environment via feces (3). Total bacterial count (TBC) in raw milk is predominantly influenced by sanitation practices. For instance, Sun et al. (4) demonstrated that farms implementing chemical teat disinfection protocols with regular acid washes achieved significantly reduced TBC levels. Contamination by genera such as Pseudomonas, Acinetobacter, Bacillus, Micrococcus, Enterobacter, Enterococcus, and Aerococcus critically threatens dairy quality (5). Notably, thermoduric microorganisms like Bacillus licheniformis frequently detected in dairy powders as well as Bacillus cereus and Bacillus subtilis, are key contributors to milk rancidity (6, 7). Conventional microbial indicators (e.g., Streptococcus, Staphylococcus, coliforms) remain vital for safety evaluation (8). Although the issue of microbial contamination in dairy products has been extensively studied, in the important dairy production base of Shandong Province in China, the relationship between the bacterial diversity in raw milk and the spread of antibiotic resistance, especially the relationship with farm management practices, is still poorly understood.

The escalating use of antimicrobials in livestock management has driven the emergence of multidrug-resistant (MDR) pathogens, critically jeopardizing dairy safety (9). Bastam et al. reported contamination rates of 60% for Staphylococcus aureus (S. aureus) and 53% for Salmonella enterica serovar Typhi (S. typhi) in dairy samples, underscoring the urgency of antimicrobial resistance (AMR) containment (10). Alarmingly, rising resistance to last-resort antibiotics (e.g., carbapenems, cephalosporins) further restricts therapeutic options (11). In China, the National Food Safety Standard (GB 31650.1-2022) establishes maximum residue limits (MRLs) for veterinary antibiotics in raw milk, including β-lactams (penicillin G: 4 μg/kg), tetracyclines (oxytetracycline: 100 μg/kg), and sulfonamides (sulfadiazine: 25 μg/kg). However, small-scale and backyard dairy operations exhibit a 44% non-compliance rate with antibiotic withdrawal periods, largely due to poor medication record-keeping (12). Inadequate milking equipment sanitation further exacerbates cross-contamination risks, leading to antibiotic residues exceeding regulatory thresholds within single milk batches (13). Indiscriminate antibiotic use not only amplifies AMR in livestock pathogens but also facilitates horizontal transfer of antibiotic resistance genes (ARGs) to human-associated bacteria via the food chain, posing dual threats to animal welfare and public health (14, 15).

The persistence of antibiotic resistance genes (ARGs) in environmental bacteria (16), such as Pseudomonas and Acinetobacter, is facilitated by mobile genetic elements that enable horizontal gene transfer (HGT) across phylogenetically distinct species. This dynamic interaction within diverse microbial communities may accelerate the evolution of multidrug-resistant strains, particularly in settings with suboptimal sanitation practices. For decades, antibiotics were routinely added to animal feed not only for therapeutic purposes but also as growth promoters, and such historical practices of overuse have directly contributed to the proliferation of ARGs (17). Although antimicrobial growth promoters have been banned in the European Union since 2006 and restricted in many countries worldwide, with China having significantly reduced their usage in recent years, the legacy of prolonged and extensive antibiotic exposure remains a key driver of persistent ARG dissemination in environmental settings. This abuse leads to ARGs being spread into the ecosystem through husbandry manure, causing environmental pollution and perpetuating the cycle of resistance (18, 19). ARG contamination now poses a critical global health challenge, demanding urgent mitigation strategies.

However, although microbial contamination in dairy products has been extensively studied, in the important dairy production base of Shandong Province, China, the relationship between bacterial diversity in raw milk and the spread of antibiotic resistance, particularly its linkage to regional farm management practices, remains poorly understood. This study aims to systematically analyze the bacterial diversity and antimicrobial resistance profiles in bulk tank milk from Shandong dairy farms. By integrating farm management data, we seek to evaluate the role of environmental bacteria as reservoirs of ARGs and their potential contribution to resistance dissemination, thereby informing region-specific interventions for dairy safety.

2 Materials and methods

2.1 Sample collection

From April to June 2022, 30 raw cow milk samples were systematically collected from 30 bulk milk tanks in five dairy farms (six samples per farm) across three different cities of Shandong Province (one farm from Taian City [TA], two farms from Dongying City [DY], and two farms from Qingdao City [QD]). Each sample was collected through aseptic mixing of raw cow milk obtained from multiple spatially distributed areas within the same bulk milk tank. All the combined samples were placed in 10-ml sterile centrifuge tubes and stored at 4°C and then transported to the laboratory for further analysis. Through the on-site epidemic prevention supervisor and the farm records, the disinfection procedures of the site, the usage of antibiotics, and the operation problems of the equipment were obtained. All experimental procedures complied with the guidelines of the Animal Welfare and Ethics Committee at Shandong Agricultural University and received approval (SDAU2022-00168).

2.2 Enumeration of bacteria

The sample was subjected to ten-fold serial dilutions in Stroke-physiological Saline Solution (SPSS) (Beyotime Biotechnology Co., Ltd., Shanghai, China), achieving a final dilution of 10−2. Subsequently, 5 μL aliquots from each dilution were inoculated onto Brain Heart Infusion agar (Beijing Solarbio Science & Technology Co., Ltd., Qingdao, China) supplemented with 5% sheep blood (BHI-blood agar) and incubated aerobically at 37°C for 24 h. Plates containing from 30 to 300 colonies were counted. Results were expressed as Colony-forming units per milliliter (CFU/ml) and were converted to log (CFU/ml).

2.3 Isolation and identification of bacteria

On BHI blood agar plates, according to differences in colony diameter, shape (round, irregular, filamentous), color (white, cream, yellow), surface texture (smooth, rough, slimy), edge features (smooth, wavy), and hemolysis patterns (β-hemolysis, α-hemolysis, non-hemolytic), all representative colony types showing distinguishable morphological differences were selected. And under the same culture conditions, these were repeatedly streaked onto BHI blood agar plates. This process was continued for 2–3 times, until each isolate obtained a pure culture exhibiting single colony morphology. Using an inoculating loop, pick a single colony and transfer it to a sterile centrifuge tube containing BHI broth, followed by incubation under appropriate conditions for a defined duration. Subsequently, add 20% glycerol to the culture, vortex thoroughly for homogenization, and store at −80°C for future use. Bacterial identification was performed using a matrix-assisted laser desorption/ionization time-of-flight mass spectrometer, (MALDI-TOF MS) (MALDI Biotyper, Bruker, Bremen, Germany) following the manufacturer’s standard operating protocols.

2.4 Antimicrobial susceptibility testing

Drug susceptibility testing was performed according to the broth microdilution, recommended by the Clinical and Laboratory Standard Institute (CLSI) (20). Two-fold serial dilutions of each antibiotic were prepared in 96-well plates (Costa) using appropriate media selected based on bacterial species, with concentration ranges spanning from 0.125 to 512 μg/mL. Incubation periods vary depending on species: 16–20 h for Enterobacteriaceae and Pseudomonas, while 20–24 h for Staphylococcus, Enterococcus, Acinetobacter, and Streptococcus. Escherichia coli (E. coli) ATCC 25922 and S. aureus ATCC 29213 were used as quality control strains. The tested panel included ten antimicrobial agents spanning seven classes of antibiotics: ceftriaxone sodium (CRO), penicillin (PCN), levofloxacin (LVX), sulfadiazine (SDZ), lincomycin (LIN), colistin sulfate (CLS), gentamicin (GM), doxycycline (DOX), neomycin (NEO), and streptomycin (SM). All the above drugs were purchased from the Shanghai MacLean Company (Shanghai, China). The result interpretation time was experimentally observed to vary significantly across bacterial species. The MIC was defined as the lowest drug concentration demonstrating complete inhibition of visible bacterial growth in the test wells. The measured MIC values were compared against the clinical breakpoint (CBP) standards established by the CLSI (CLSI M100 35th ed., 2025) and the European Committee on Antimicrobial Susceptibility Testing (EUCAST, v 12.0 Breakpoint Tables) (20, 21). In cases where CBPs for specific antimicrobial agents are not defined by CLSI or EUCAST, an integrated assessment utilizing epidemiological cutoff values (ECVs [CLSI terminology]/ECOFFs [EUCAST terminology], where resistance was presumptively classified based on non-wild type (NWT) status) and resistance gene detection should be implemented (21, 22). Isolates exhibiting MICs exceeding the ECOFF and demonstrating the presence of a resistance-associated gene known to confer resistance to the specific antimicrobial agent were provisionally classified as candidate resistant strains. MDR bacteria are defined as those that are resistant to at least one antimicrobial agent in each of three or more antimicrobial classes (23).

2.5 DNA extraction

The boiling method was used to extract DNA (24). Single colonies were picked and incubated to the logarithmic phase at 37°C, 200 rpm in Luria-Bertani (LB) broth medium. The bacterial culture was centrifuged at 12,000 × g for 30 s at 4°C to obtain cellular precipitates. The precipitates were mixed with Tris-EDTA (TE) buffer (Thermo Fisher Scientific Co., Ltd., Shanghai, China) and subjected to cellular disruption by boiling at 98°C for 10–15 min. The tubes were immediately placed on ice for ten minutes and then centrifuged at 13,000 rpm for five minutes to collect the clarified supernatant, thereby harvesting DNA from each isolate. The collected samples were stored at −20°C until further processing.

2.6 Screening for antibiotic resistance genes

The 129 isolates were screened for the presence of nine antibiotic resistance genes by PCR, including β-lactam resistance genes (blaKPC and blaNDM-1), sulfonamide resistance gene (sul1), tetracycline resistance gene (tet(M)), aminoglycoside resistance genes (ant(4′)-Ia and aph(2″)-Ic), quinolone resistance gene (qnrS), lincosamide resistance gene (lnu(B)), and polymyxin resistance genes (mcr-1). The PCR was performed with 25 μL volumes composed of 12.5 μL of 2 × Taq Master MIX (Vazyme Biotech Co., Ltd., Nanjing, China), 2 μL of primers, 1 μL of 100 ng of DNA template, and 9.5 μL DEPC-Treated Water (NCM Biotech Co., Ltd., Suzhou, China). All primers are listed in Table 1.

Table 1
www.frontiersin.org

Table 1. Primer sequences used in the study.

2.7 Statistical analysis

Data organization was completed using Microsoft Excel 2021 (Microsoft Corporation, Redmond, WA, United States), and graphs were generated using Origin 2022 (Origin Lab Corporation, Northampton, MA, United States). The data are expressed as means ± SD. Statistical significance was determined based on one-way analysis of variance (ANOVA) in appropriate conditions using GraphPad Prism 8 software. Significance was determined at p < 0.05.

3 Results

3.1 Total bacteria count in raw milk samples

The TBC of raw cow milk samples in five dairy farms ranged from 3.94 to 6.68 log CFU/ml (Figure 1A). Significant inter-farm variation was observed (p < 0.05). QD1 dairy farm had the highest TBC, with values ranging from 6.43 to 6.68 log CFU/ml, while TA dairy farm had the lowest TBC, ranging from 3.94 to 4.27 log CFU/ml. This discrepancy may be attributed to differences in equipment sanitation protocols across farms. Although all farms implemented standardized chemical disinfection protocols before and after daily operations, TA dairy farm uniquely adopted supplementary thermal disinfection (≥85°C for 15 min) followed by thorough drying for heat-resistant components of milk storage tanks. In contrast, DY2 and QD1 dairy farms showed suboptimal cooling system efficiency in their storage tanks, potentially contributing to prolonged microbial proliferation. Overall, TBC in DY2 and QD1 dairy farms were higher than that in other dairy farms.

Figure 1
Image depicting four data visualizations labeled A to D. A shows a box plot with Log CFU/mL measurements for farms TA, DY1, DY2, QD1, and QD2. B is a stacked bar chart illustrating the relative abundance of various bacterial species across the same farms. C is another stacked bar chart with similar data but focusing on specific genera like Psychrobacter, Pantoea, and Bacillus. D is a Venn diagram showing the overlap of bacterial species among the farms, with numbers indicating unique and shared species.

Figure 1. Total bacterial count (TBC) from Taian (TA), Dongying1 (DY1), Dongying2 (DY2), Qingdao1 (QD1) and Qingdao2 (QD2) dairy farms (A). Relative abundance of bacterial species (B) and family (C) from these dairy farms. Venn diagram of correlation of bacterial species in different dairy farms (D).

3.2 Identification of bacteria

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) identified 129 isolates spanning 37 species and 13 families. Notably, pathogenic bacteria associated with bovine mastitis (Staphylococcus xylosus, Staphylococcus chromogenes, Streptococcus parauberis) were isolated, accounting for 8.53% (11/129) of the total isolates. Spoilage-related bacteria (Pseudomonas lundensis, Pseudomonas proteolytica, Pseudomonas gessardii, Pseudomonas fragi) constituted 13.95% (18/129), while opportunistic pathogens (Acinetobacter lwoffii, Acinetobacter johnsonii, E. coli) represented 20.16% (26/129) of the identified microbiota (Table 2). Geographically, TA farm harbored 14 bacterial species, with Bacillus marisflavi as the predominant species. The pathogenic bacteria belonging to the Staphylococcus genus accounted for 20.8% (5/24) (p < 0.05; Figure 1B). DY1 and DY2 farms yielded 10 and 9 species, respectively, dominated by Lactococcus lactis (DY1) and Pseudomonas gessardii (DY2) (p < 0.05). DY1 farm exhibited a notably high prevalence of spoilage-related bacteria 13.6% (3/22), followed by opportunistic pathogens 22.7% (5/22). In contrast, DY2 farm demonstrated comparable microbial patterns, with spoilage-related bacteria constituting 38.1% (8/21) and opportunistic pathogens accounting for 38.1% (8/21). QD1 farm exhibited 13 species, primarily Macrococcus caseolyticus, while QD2 farm isolated 8 species, dominated by Acinetobacter johnsonii (p < 0.05). Notably, spoilage-related bacteria constituted 16.7% (6/36) and opportunistic pathogens represented 25% (9/36) of the microbiota in QD1. In contrast, QD2 showed markedly lower proportions of spoilage-related bacteria 3.8% (1/26) and opportunistic pathogens 15.4% (4/26), but exhibited a higher prevalence of mastitis-associated pathogens 23.1% (6/26). At the family level, TA isolates (n = 24) belonged to 6 families, predominantly Bacillus. DY1/DY2 farms each contained 6 families, with Pseudomonas as the dominant taxon. In QD, QD1 and QD2 farms harbored 7 and 5 families, dominated by Macrococcus and Acinetobacter, respectively (Figure 1C). Notably, A. johnsonii, E. coli, Lactococcus lactis, Pseudomonas fragi, Pseudomonas lundensis, and Rothia endophytica were exclusively detected in DY/QD farms, absent in TA (Figure 1D).

Table 2
www.frontiersin.org

Table 2. Bacterial isolation and identification across sampled dairy farms.

3.3 Antimicrobial resistance

A representative subset of 87 bacterial strains encompassing all identified species/genera was selected from 129 isolates. Antimicrobial susceptibility testing (AST) against ten agents from seven antibiotic classes (with intermediate categorized as non-resistant) revealed significant inter-agent resistance heterogeneity. For isolates lacking established clinical breakpoints for specific antibiotics, data were excluded from final resistance rate calculations (Figures 2AD). SDZ exhibited the highest resistance rate 53.16% (42/79), with marked regional disparities: TA and DY isolates showed high PCN resistance, whereas QD isolates remained susceptible. Region-specific resistance profiling further demonstrated that: TA isolates exhibited prevalent resistance to CRO, SDZ, and PCN, with mastitis-associated S. xylosus showing strong resistance to CRO, SDZ, and PCN, and S. chromogenes displaying high resistance to SDZ, and PCN; DY isolates displayed elevated resistance to SDZ, PCN, and SM, notably with opportunistic A. johnsonii showing SM resistance and E. coli resistant to CRO, SDZ, DOX, PCN, and SM; QD isolates were predominantly resistant to SDZ and SM, where both mastitis-associated S. chromogenes and spoilage-related P. lundensis exhibited strong dual resistance to these agents. Among the tested isolates, only those originating from QD exhibited 50.6% resistance to SM (42/83), as determined by exceeding ECOFFs and confirming the presence of resistance associated genetic determinants in the absence of clinical breakpoints. Within the 23% (20/87) MDR isolates identified (Figure 3A), 70% (14/20) were classified as MDR candidates based on resistance to at least three antimicrobial classes using this integrative assessment. This resistance landscape suggests potential HGT of resistance determinants from environmental commensals to pathogens via mobile genetic elements, particularly in farms with suboptimal hygiene conditions where complex microbial communities provide ecological niches for HGT-driven MDR evolution.

Figure 2
Chart and heatmaps showing antibiotic resistance. Bar chart (A) displays resistance rates for various antibiotics, with SDZ and SM showing higher resistance. Heatmaps (B, C, D) illustrate resistance patterns across bacterial species against different antibiotics. Blue, red, yellow, brown, and green indicate varying resistance levels.

Figure 2. The overall resistance rate of the isolated bacteria to ten antibiotics. Ceftriaxone sodium (CRO), penicillin (PCN), levofloxacin (LVX), sulfadiazine (SDZ), lincomycin (LIN), colistin sulfate (CLS), gentamicin (GM), doxycycline (DOX), neomycin (NEO), and streptomycin (SM) (A). Isolated heat resistance spectrum in Taian (TA), Dongying (DY) and Qingdao (QD) cities. Red: resistant/Blue: susceptible/Yellow: intermediate/Green: non-wild type (NWT)/Brown: clinical breakpoint not established (B–D).

Figure 3
Two graphs display data on antibiotic resistance. Graph A is a bar chart showing the number of multidrug-resistant strains, with 13 presumptive multidrug-resistant and 5 multidrug-resistant for 4 antibiotics. Graph B is a horizontal bar chart depicting total positive numbers of antibiotic resistance genes, highlighting

Figure 3. The number of multidrug resistance (A) and antibiotic resistance genes (B).

3.4 Antimicrobial resistance gene

Genomic analysis of 129 isolates revealed nine ARGs spanning seven drug classes: sulfonamide (sul1), aminoglycosides [ant(4′)-Ia, aph(2″)-Ic], tetracycline (tet(M)), quinolone (qnrS), β-lactam (blaKPC, blaNDM-1), polymyxin (mcr-1), and lincosamide (lnu(B)). The detection rate of sul1 (70.5%) surpassed other ARGs, followed by ant (4′)-Ia (54.3%), while mcr-1, blaNDM-1, and lnu(B) were absent in all isolates (Figure 3B). Geographically, sul1 (70.8%), ant(4′)-Ia (54.2%), and tet(M) (54.2%) dominated TA isolates; DY isolates exhibited higher ant(4′)-Ia (74.4%) and sul1 (60.5%); QD isolates showed sul1 (63.5%) and ant(4′)-Ia (44.5%) prevalence (Table 3). Notably, the ARG distribution demonstrated concordance with regional resistance phenotypes. In TA, the high prevalence of sul1 (70.8%) correlated strongly with observed SDZ resistance (73.7%). In DY, elevated ant(4′)-Ia detection (74.4%) aligned with high SM resistance across tested isolates. For QD, the co-occurrence of sul1 (63.5%) and ant(4′)-Ia (44.5%) paralleled the dual SDZ/SM resistance patterns identified in both pathogenic and spoilage-related bacterial populations.

Table 3
www.frontiersin.org

Table 3. The positive rate of antibiotic resistance genes in different cities.

The congruence between resistance genotypes and phenotypes not only confirms the functional validity of ARGs but also elucidates potential dissemination pathways traversing environmental reservoirs, commensal microbiota, and pathogenic communities. These findings collectively emphasize the imperative for multidimensional containment strategies targeting critical junctions of antimicrobial resistance transmission networks.

4 Discussion

The dairy farm environment harbors diverse bacterial species, including opportunistic pathogens (e.g., Klebsiella pneumoniae (25)), foodborne pathogens (S. aureus (26), Campylobacter (27)), and toxin-producing strains. These bacteria can be directly transmitted through raw milk or enter the dairy product chain through fecal contamination (28). Microbial contamination in bulk tank milk poses significant food safety risks. While enterotoxigenic S. aureus was not detected here, its documented role in dairy-associated poisoning underscores contamination hazards. MDR Campylobacter and Salmonella can spread through raw milk, making post-infection treatment difficult (29, 30). The complex microbial community in raw milk also provides an environment for the transfer of ARGs. For instance, Enterococcus can transfer the tet(M) gene to Lactobacillus spp. through conjugation. Cross-contamination facilitates metabolic exchanges among different bacterial communities and accelerates the spread of ARGs (28).

The microbial diversity in raw milk is profoundly influenced by farm hygiene practices (31, 32). TA farm’s low TBC (3.94 log CFU/mL) and limited MDR isolates correlates with its rigorous thermal disinfection protocols, which reduce bacterial load and potentially limit HGT by disrupting biofilms. Regional variations also stem from environmental adaptation (33): QD farms’ inadequate cooling systems promoted psychrotrophic Pseudomonas growth (Figure 1C), which carried sul1 (60.5%) and acted as ARG reservoirs (34). Optimizing sanitation protocols could disrupt these ecological niches, thereby curbing HGT-driven resistance dissemination (35).

AST revealed alarming resistance rates to SDZ (53.2%) and SM (50.6%) (Figure 2A), consistent with historical overuse of first-line antibiotics in livestock (36). Through interviews with the on-site epidemic prevention supervisor and reviews of farm records, we found that the high SDZ resistance correlated with its routine prophylactic use against mastitis and respiratory infections in sampled farms. In QD farms, inadequate resting space for cows increased the risk of udder injury. Combined with inconsistent feed sourcing and underdeveloped disease prevention systems, this led to reliance on SM and CLS for infection control. This explains the exclusive detection of SM resistance (50.6%) in QD isolates and elevated CLS resistance (24.3%). In contrast, TA and DY farms exhibited no recent SM usage, highlighting regional disparities in antibiotic practices. Integration of ECOFFs and PCR identified 23.0% (20/87) of isolates as MDR/presumptive MDR strains (Figure 3A).

However, certain limitations should be acknowledged. This approach may overlook novel resistance mechanisms and is limited by ECOFF variability across bacterial populations. Conducting a single sampling study only on dairy farms within Shandong Province will limit its general applicability in a larger geographical area and will not be able to accurately assess the dynamics of resistance evolution. Relying on culturable bacteria for screening high-priority genes may underestimate the overall microbial diversity, resulting in an inability to comprehensively evaluate the complete set of resistance genes.

Regional resistance patterns mirrored global trends, with QD’s SM resistance comparable to Turkish dairy isolates (37); while penicillin (PCN) resistance (48%) was lower than Iranian reports (38). Notably, 20.7% of isolates exhibited resistance to 3–4 drug classes, underscoring urgent needs for targeted interventions. Crucially, even partial compliance with antibiotic restrictions can mitigate resistance emergence, but requires synergism with hygiene improvements to block HGT pathways. The dominance of sul1 (70.5%) and tet(M) (54.2%) in environmental bacteria aligns with historical sulfonamide/tetracycline applications in livestock. Selective pressure from prolonged antibiotic use likely enriched these ARGs in bacterial populations (39), facilitating HGT via mobile genetic elements such as plasmids and integrons (40). DY farms’ high ant(4′)-Ia prevalence (74.4%) in Pseudomonas may enable aminoglycoside resistance transfer to co-existing pathogens like E. coli, particularly under suboptimal sanitation conditions that promote mixed microbial communities. Previous studies have monitored the changes in the pool of AMR genes in raw milk and cheese from farm to consumer (41).

In this study, we did not detect the specific genes such as mcr-1, lnu(B) and blaNDM-1. This result may be influenced by various factors, including regulatory restrictions, sampling scope, geographical differences, and the specific bacterial populations present in the sampled milk tanks. However, the persistence of the sulfonamide resistance gene sul1 and the tetracycline resistance gene tet(M) indicates the tenacity of horizontal transferable resistance determinants, emphasizing the importance of continuous vigilance and providing valuable insights for public health monitoring. While our data are specific to cow milk, the detected AMR genes have been similarly reported in small ruminant and buffalo dairy systems under intensive farming conditions (42, 43).

In conclusion, our findings demonstrate the interconnected relationships among microbial diversity, antibiotic, and ARG dissemination within dairy farming systems in Shandong Province. The systematic integration of TBC metrics, AMR profiles, and ARG distribution patterns provides actionable strategies for reducing AMR risks. These evidence-based insights establish a critical foundation for developing targeted intervention protocols and implementing genomic surveillance of multidrug-resistant pathogens, thereby mitigating AMR transmission through dairy products.

5 Conclusion

This study reported the key issues of microbial contamination and AMR characteristics in cow bulk tank milk from dairy farms in Shandong Province, China. The combination of the antibiotic resistance status of the strains and the regional differences observed in farm practices further emphasizes the importance of rational antibiotic use to help alleviate the selection pressure for antibiotic resistance. At the same time, optimizing environmental hygiene management remains crucial for reducing the risk of bacterial transmission and the transfer of antibiotic resistance genes in the farm environment.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The animal studies were approved by the Animal Care and Ethics Committee of Shandong Agricultural University (SDAU2022-00168). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author contributions

YQ: Investigation, Writing – original draft. ZL: Data curation, Writing – original draft. ZM: Writing – review & editing, Software. XW: Validation, Writing – review & editing. HC: Methodology, Writing – review & editing. ML: Writing – review & editing, Validation. CQ: Validation, Writing – review & editing. PZ: Investigation, Writing – review & editing, Supervision, Resources, Project administration. YL: Project administration, Writing – review & editing, Resources, Methodology. JL: Writing – review & editing, Funding acquisition, Supervision, Resources, Project administration, Methodology.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. The project was supported by Shandong Natural Science Foundation of China (ZR2024MC079, ZR2023MC025).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Fusco, V, Chieffi, D, Fanelli, F, Logrieco, AF, Cho, GS, Kabisch, J, et al. Microbial quality and safety of milk and milk products in the 21st century. Compr Rev Food Sci Food Saf. (2020) 19:2013–49. doi: 10.1111/1541-4337.12568

PubMed Abstract | Crossref Full Text | Google Scholar

2. Cao, H, Yan, Y, Wang, L, Dong, L, Pang, X, Tang, S, et al. High-throughput sequencing reveals bacterial diversity in raw Milk production environment and production chain in Tangshan City of China. Food Sci Anim Resour. (2021) 41:452–67. doi: 10.5851/kosfa.2021.e10

PubMed Abstract | Crossref Full Text | Google Scholar

3. Oliver, SP, Jayarao, BM, and Almeida, RA. Foodborne pathogens in milk and the dairy farm environment: food safety and public health implications. Foodborne Pathog Dis. (2005) 2:115–29. doi: 10.1089/fpd.2005.2.115

PubMed Abstract | Crossref Full Text | Google Scholar

4. Sun, L, Lundh, A, Hojer, A, Bernes, G, Nilsson, D, Johansson, M, et al. Milking system and premilking routines have a strong effect on the microbial community in bulk tank milk. J Dairy Sci. (2022) 105:123–39. doi: 10.3168/jds.2021-20661

PubMed Abstract | Crossref Full Text | Google Scholar

5. Straley, BA, Donaldson, SC, Hedge, NV, Sawant, AA, Srinivasan, V, Oliver, SP, et al. Public health significance of antimicrobial-resistant gram-negative bacteria in raw bulk tank milk. Foodborne Pathog Dis. (2006) 3:222–33. doi: 10.1089/fpd.2006.3.222

PubMed Abstract | Crossref Full Text | Google Scholar

6. Ribeiro Junior, JC, Tamanini, R, De Oliveira, ALM, Alfieri, AA, and Beloti, V. Genetic diversity of thermoduric spoilage microorganisms of milk from Brazilian dairy farms. J Dairy Sci. (2018) 101:6927–36. doi: 10.3168/jds.2017-13948

PubMed Abstract | Crossref Full Text | Google Scholar

7. Miller, RA, Kent, DJ, Watterson, MJ, Boor, KJ, Martin, NH, and Wiedmann, M. Spore populations among bulk tank raw milk and dairy powders are significantly different. J Dairy Sci. (2015) 98:8492–504. doi: 10.3168/jds.2015-9943

PubMed Abstract | Crossref Full Text | Google Scholar

8. Gillespie, BE, Lewis, MJ, Boonyayatra, S, Maxwell, ML, Saxton, A, Oliver, SP, et al. Short communication: evaluation of bulk tank milk microbiological quality of nine dairy farms in Tennessee. J Dairy Sci. (2012) 95:4275–9. doi: 10.3168/jds.2011-4881

PubMed Abstract | Crossref Full Text | Google Scholar

9. Sohail, M, Rashid, A, Aslam, B, Waseem, M, Shahid, M, Akram, M, et al. Antimicrobial susceptibility of Acinetobacter clinical isolates and emerging antibiogram trends for nosocomial infection management. Rev Soc Bras Med Trop. (2016) 49:300–4. doi: 10.1590/0037-8682-0111-2016

PubMed Abstract | Crossref Full Text | Google Scholar

10. Bastam, MM, Jalili, M, Pakzad, I, Maleki, A, and Ghafourian, S. Pathogenic bacteria in cheese, raw and pasteurised milk. Vet Med Sci. (2021) 7:2445–9. doi: 10.1002/vms3.604

PubMed Abstract | Crossref Full Text | Google Scholar

11. Rizzo, K, Horwich-Scholefield, S, and Epson, E. Carbapenem and cephalosporin resistance among Enterobacteriaceae in healthcare-associated infections, California, Usa(1). Emerg Infect Dis. (2019) 25:1389–93. doi: 10.3201/eid2507.181938

PubMed Abstract | Crossref Full Text | Google Scholar

12. Dong, J, Wang, L, Liu, L, Zhang, Y, Wu, Q, and Zhao, J. How policy advocacy promotes regulated antibiotic use: evidence from meat duck farmers of China. One Health. (2024) 19:100900. doi: 10.1016/j.onehlt.2024.100900

PubMed Abstract | Crossref Full Text | Google Scholar

13. Zapata-Salas, R, Guarin, JF, and Rios-Osorio, LA. Consumption and informal trade of Milk in the north of Antioquia (Colombia). Vet Med Int. (2024) 2024:6644328. doi: 10.1155/2024/6644328

PubMed Abstract | Crossref Full Text | Google Scholar

14. Wang, W, An, T, Li, G, Xia, D, Zhao, H, Yu, JC, et al. Earth-abundant Ni2P/g-C3N4 lamellar nanohydrids for enhanced photocatalytic hydrogen evolution and bacterial inactivation under visible light irradiation. Appl Catal B Environ. (2017) 217:570–80. doi: 10.1016/j.apcatb.2017.06.027

Crossref Full Text | Google Scholar

15. Van Boeckel, TP, Brower, C, Gilbert, M, Grenfell, BT, Levin, SA, Robinson, TP, et al. Global trends in antimicrobial use in food animals. Proc Natl Acad Sci USA. (2015) 112:5649–54. doi: 10.1073/pnas.1503141112

PubMed Abstract | Crossref Full Text | Google Scholar

16. Jian, Z, Zeng, L, Xu, T, Sun, S, Yan, S, Yang, L, et al. Antibiotic resistance genes in bacteria: occurrence, spread, and control. J Basic Microbiol. (2021) 61:1049–70. doi: 10.1002/jobm.202100201

PubMed Abstract | Crossref Full Text | Google Scholar

17. Sarmah, AK, Meyer, MT, and Boxall, AB. A global perspective on the use, sales, exposure pathways, occurrence, fate and effects of veterinary antibiotics (vas) in the environment. Chemosphere. (2006) 65:725–59. doi: 10.1016/j.chemosphere.2006.03.026

PubMed Abstract | Crossref Full Text | Google Scholar

18. Zhu, YG, Johnson, TA, Su, JQ, Qiao, M, Guo, GX, Stedtfeld, RD, et al. Diverse and abundant antibiotic resistance genes in Chinese swine farms. Proc Natl Acad Sci USA. (2013) 110:3435–40. doi: 10.1073/pnas.1222743110

PubMed Abstract | Crossref Full Text | Google Scholar

19. Zhang, QQ, Ying, GG, Pan, CG, Liu, YS, and Zhao, JL. Comprehensive evaluation of antibiotics emission and fate in the river basins of China: source analysis, multimedia modeling, and linkage to bacterial resistance. Environ Sci Technol. (2015) 49:6772–82. doi: 10.1021/acs.est.5b00729

PubMed Abstract | Crossref Full Text | Google Scholar

20. CLSI ed. Performance standards for antimicrobial susceptibility testing. 35th ed. Wayne, PA: CLSI Supplement M100 Clinical and Laboratory Standards Institute (2025).

Google Scholar

21. Breakpoint Tables for Interpretation of MICs and Zone Diameters Version 12.0. European Committee on Antimicrobial Susceptibility Testing EUCAST clinical breakpoint tables for interpretation of MICs and zone diameters, version 12.0. Sweden: European Committee on Antimicrobial Susceptibility Testing (2022).

Google Scholar

22. Kahlmeter, G, and Turnidge, J. Wild-type distributions of minimum inhibitory concentrations and epidemiological cut-off values-laboratory and clinical utility. Clin Microbiol Rev. (2023) 36:e0010022. doi: 10.1128/cmr.00100-22

PubMed Abstract | Crossref Full Text | Google Scholar

23. Sweeney, MT, Lubbers, BV, Schwarz, S, and Watts, JL. Applying definitions for multidrug resistance, extensive drug resistance and pandrug resistance to clinically significant livestock and companion animal bacterial pathogens. J Antimicrob Chemother. (2018) 73:1460–3. doi: 10.1093/jac/dky043

PubMed Abstract | Crossref Full Text | Google Scholar

24. Peng, X, Yu, KQ, Deng, GH, Jiang, YX, Wang, Y, Zhang, GX, et al. Comparison of direct boiling method with commercial kits for extracting fecal microbiome Dna by Illumina sequencing of 16S rrna tags. J Microbiol Methods. (2013) 95:455–62. doi: 10.1016/j.mimet.2013.07.015

PubMed Abstract | Crossref Full Text | Google Scholar

25. Adhyaru, HJ, Mangroliya, DB, Bhalu, AB, Kabariya, JH, and Ramani, VM. Genomic insights into prevalence of virulence and multi drug resistance genes in milk borne Klebsiella pnuemoniae: face of emerging resistance to last resort antibiotics. Microb Pathog. (2024) 196:106936. doi: 10.1016/j.micpath.2024.106936

PubMed Abstract | Crossref Full Text | Google Scholar

26. Huang, J, Zhang, W, Sun, B, Jiang, Q, Cao, Y, Shang, J, et al. Genetic diversity, antibiotic resistance, and virulence characteristics of Staphylococcus aureus from raw milk over 10 years in Shanghai. Int J Food Microbiol. (2023) 401:110273. doi: 10.1016/j.ijfoodmicro.2023.110273

PubMed Abstract | Crossref Full Text | Google Scholar

27. Ahmadi, Z, Pakbin, B, Kazemi, M, Rahimi, Z, and Mahmoudi, R. Genotyping and antibiotic susceptibility of campylobacter species isolated from raw milk samples in Qazvin, Iran. BMC Res Notes. (2023) 16:314. doi: 10.1186/s13104-023-06576-9

PubMed Abstract | Crossref Full Text | Google Scholar

28. Morandi, S, Silvetti, T, Lopreiato, V, Piccioli-Cappelli, F, Trevisi, E, and Brasca, M. Biodiversity and antibiotic resistance profile provide new evidence for a different origin of enterococci in bovine raw milk and feces. Food Microbiol. (2024) 120:104492. doi: 10.1016/j.fm.2024.104492

PubMed Abstract | Crossref Full Text | Google Scholar

29. Aguiar, RAC, Ferreira, FA, Rubio Cieza, MY, Silva, NCC, Miotto, M, Carvalho, MM, et al. Staphylococcus aureus isolated from traditional artisanal raw Milk cheese from southern Brazil: diversity, virulence, and antimicrobial resistance profile. J Food Prot. (2024) 87:100285. doi: 10.1016/j.jfp.2024.100285

PubMed Abstract | Crossref Full Text | Google Scholar

30. Golas-Pradzynska, M, and Rola, JG. Occurrence and antimicrobial resistance of enterococci isolated from goat's Milk. J Vet Res. (2021) 65:449–55. doi: 10.2478/jvetres-2021-0071

PubMed Abstract | Crossref Full Text | Google Scholar

31. Cui, Y, Song, K, Liu, X, Xu, H, Wang, X, Cheng, G, et al. Research on bacterial diversity and antibiotic resistance in the dairy farm environment in a part of Shandong Province. Animals (Basel). (2024) 14:160. doi: 10.3390/ani14010160

PubMed Abstract | Crossref Full Text | Google Scholar

32. Yan, G, Wang, S, Cui, Y, Xue, K, Liu, Y, and Liu, J. Bacterial diversity and antimicrobial resistance of microorganisms isolated from teat cup liners in dairy farms in Shandong Province, China. Animals (Basel). (2024) 14:2167. doi: 10.3390/ani14152167

PubMed Abstract | Crossref Full Text | Google Scholar

33. Li, Z, Peng, C, Zhang, G, Shen, Y, Zhang, Y, Liu, C, et al. Prevalence and characteristics of multidrug-resistant Proteus mirabilis from broiler farms in Shandong Province, China. Poult Sci. (2022) 101:101710. doi: 10.1016/j.psj.2022.101710

PubMed Abstract | Crossref Full Text | Google Scholar

34. Marchand, S, Heylen, K, Messens, W, Coudijzer, K, de Vos, P, Dewettinck, K, et al. Seasonal influence on heat-resistant proteolytic capacity of Pseudomonas lundensis and Pseudomonas fragi, predominant milk spoilers isolated from Belgian raw milk samples. Environ Microbiol. (2009) 11:467–82. doi: 10.1111/j.1462-2920.2008.01785.x

PubMed Abstract | Crossref Full Text | Google Scholar

35. Munsch-Alatossava, P, Jaaskelainen, S, Alatossava, T, and Gauchi, JP. N(2) gas Flushing limits the rise of antibiotic-resistant bacteria in bovine raw Milk during cold storage. Front Microbiol. (2017) 8:655. doi: 10.3389/fmicb.2017.00655

PubMed Abstract | Crossref Full Text | Google Scholar

36. Yang, Y, Peng, Y, Jiang, J, Gong, Z, Zhu, H, Wang, K, et al. Isolation and characterization of multidrug-resistant Klebsiella pneumoniae from raw cow milk in Jiangsu and Shandong provinces, China. Transbound Emerg Dis. (2021) 68:1033–9. doi: 10.1111/tbed.13787

PubMed Abstract | Crossref Full Text | Google Scholar

37. Ozdemir, R, and Tuncer, Y. Detection of antibiotic resistance profiles and aminoglycoside-modifying enzyme (Ame) genes in high-level aminoglycoside-resistant (Hlar) enterococci isolated from raw milk and traditional cheeses in Turkey. Mol Biol Rep. (2020) 47:1703–12. doi: 10.1007/s11033-020-05262-4

PubMed Abstract | Crossref Full Text | Google Scholar

38. Pahlavanzadeh, S, Khoshbakht, R, Kaboosi, H, and Moazamian, E. Antibiotic resistance and phylogenetic comparison of human, pet animals and raw milk Staphylococcus aureus isolates. Comp Immunol Microbiol Infect Dis. (2021) 79:101717. doi: 10.1016/j.cimid.2021.101717

PubMed Abstract | Crossref Full Text | Google Scholar

39. Pereira, AR, Paranhos, AGO, De Aquino, SF, and Silva, SQ. Distribution of genetic elements associated with antibiotic resistance in treated and untreated animal husbandry waste and wastewater. Environ Sci Pollut Res Int. (2021) 28:26380–403. doi: 10.1007/s11356-021-13784-y

PubMed Abstract | Crossref Full Text | Google Scholar

40. Su, Y, Yu, CY, Tsai, Y, Wang, SH, Lee, C, and Chu, C. Fluoroquinolone-resistant and extended-spectrum β-lactamase-producing Escherichia coli from the milk of cows with clinical mastitis in southern Taiwan. J Microbiol Immunol Infect. (2016) 49:892–901. doi: 10.1016/j.jmii.2014.10.003

PubMed Abstract | Crossref Full Text | Google Scholar

41. Kerek, A, Nemeth, V, Szabo, A, Papp, M, Bányai, K, Kardos, G, et al. Monitoring changes in the antimicrobial-resistance gene set (Arg) of raw Milk and dairy products in a cattle farm, from production to consumption. Vet Sci. (2024) 11:265. doi: 10.3390/vetsci11060265

PubMed Abstract | Crossref Full Text | Google Scholar

42. Veloo, Y, Rajendiran, S, Zakaria, Z, Ismail, R, Rahman, SA, Mansor, R, et al. Prevalence and antimicrobial resistance patterns of Escherichia coli in the environment, cow dung, and Milk of Selangor dairy farms. Antibiotics (Basel). (2025) 14:137. doi: 10.3390/antibiotics14020137

PubMed Abstract | Crossref Full Text | Google Scholar

43. Nelli, A, Voidarou, CC, Venardou, B, Fotou, K, Tsinas, A, Bonos, E, et al. Antimicrobial and methicillin resistance pattern of potential mastitis-inducing Staphylococcus aureus and coagulase-negative staphylococci isolates from the mammary secretion of dairy goats. Biology (Basel). (2022) 11:1591. doi: 10.3390/biology11111591

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: bacterial diversity, antibiotic resistance, bulk tank milk, dairy farm, antibiotic resistance genes

Citation: Qi Y, Lu Z, Meng Z, Wang X, Chen H, Li M, Qu C, Zhang P, Liu Y and Liu J (2025) Diversity and antibiotic resistance of cultivable bacteria in bulk tank milk from dairy farms in Shandong Province, China. Front. Vet. Sci. 12:1649876. doi: 10.3389/fvets.2025.1649876

Received: 19 June 2025; Accepted: 08 July 2025;
Published: 23 July 2025.

Edited by:

Xiangdong Li, Yangzhou University, China

Reviewed by:

Yan-Dong Tang, Chinese Academy of Agricultural Sciences, China
Rui Li, Harvard Medical School, United States

Copyright © 2025 Qi, Lu, Meng, Wang, Chen, Li, Qu, Zhang, Liu and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Pu Zhang, enA4MTk4NDIzQDE2My5jb20=; Yongxia Liu, bGl1eW9uZ3hpYUBzZGF1LmVkdS5jbg==; Jianzhu Liu, bGl1anpAc2RhdS5lZHUuY24=

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.