Skip to main content

Pathogenicity assessment of Arcobacter butzleri isolated from Canadian agricultural surface water



Water is considered a source for the transmission of Arcobacter species to both humans and animals. This study was conducted to assess the prevalence, distribution, and pathogenicity of A. butzleri strains, which can potentially pose health risks to humans and animals. Cultures were isolated from surface waters of a mixed-use but predominately agricultural watershed in eastern Ontario, Canada. The detection of antimicrobial resistance (AMR) and virulence-associated genes (VAGs), as well as enterobacterial repetitive intergenic consensus-polymerase chain reaction (ERIC-PCR) assays were performed on 913 A. butzleri strains isolated from 11 agricultural sampling sites.


All strains were resistant to one or more antimicrobial agents, with a high rate of resistance to clindamycin (99%) and chloramphenicol (77%), followed by azithromycin (48%) and nalidixic acid (49%). However, isolates showed a significantly (p < 0.05) high rate of susceptibility to tetracycline (1%), gentamycin (2%), ciprofloxacin (4%), and erythromycin (5%). Of the eight VAGs tested, ciaB, mviN, tlyA, and pldA were detected at high frequency (> 85%) compared to irgA (25%), hecB (19%), hecA (15%), and cj1349 (12%) genes. Co-occurrence analysis showed A. butzleri strains resistant to clindamycin, chloramphenicol, nalidixic acid, and azithromycin were positive for ciaB, tlyA, mviN and pldA VAGs. ERIC-PCR fingerprint analysis revealed high genetic similarity among strains isolated from three sites, and the genotypes were significantly associated with AMR and VAGs results, which highlight their potential environmental ubiquity and potential as pathogenic.


The study results show that agricultural activities likely contribute to the contamination of A. butzleri in surface water. The findings underscore the importance of farm management practices in controlling the potential spread of A. butzleri and its associated health risks to humans and animals through contaminated water.


1. Arcobacter butzleri was highly prevalent in agricultural surface water.

2. Strains were highly resistant to antimicrobial agents.

3. Virulence genes were detected at a high frequency.

4. High genetic similarity and association among strains from various sites.

Peer Review reports


Over the past few decades, it has been determined that untreated water is an important source of Arcobacter species which are considered emerging human pathogens causing enteritis [1, 2]. Arcobacter species are gram-negative and are considered the fourth most important pathogenic member of the Campylobacteraceae family [3, 4]. Among these, A. butzleri, A. cryaerophilus, and A. skirrowii have been isolated from humans and animals [1]. A. butzleri is known to cause various human infections making it clinically important. A. butzleri strains isolated from human and animal sources have been reported and classified as emerging human pathogens (EHPs) by the International Commission on Microbiological Specifications for Foods [5].

A. butzleri and A. cryaerophilus can cause abdominal pain, gastroenteritis and acute diarrhea or prolonged watery diarrhea in humans. They can even cause bacteremia, abortion and mastitis in cattle and swine [6,7,8]. The main sources of A. butzleri transmission include ingestion of contaminated water, food contamination, and fecal shedding from livestock such as cattle [1, 7]. Uncooked or minimally processed foods also pose a high risk in the transmission of infection and can consequently be a serious hazard to public health [9]. Since contaminated water is considered one of the possible routes of infection in humans and animals [7], examining the prevalence of Arcobacter species in water is useful to better understand the source and transmission process of these pathogens [10].

Many Arcobacter species are known to acquire resistance against a wide range of commonly used antibiotics due to lateral gene transfer; hence, various antimicrobial agents are proven ineffective in managing A. butzleri infections [7, 11]. The antimicrobial resistance (AMR) profile of A. butzleri previously reported has shown high resistance against erythromycin and ciprofloxacin, which are commonly used antibiotics for treating Campylobacter infections. However, some studies have revealed AMR in A. butzleri against amoxicillin, nalidixic acid, gentamicin, clindamycin, azithromycin, ciprofloxacin, metronidazole, carbenicillin, and cefoperazone, but maximum susceptibility to fluoroquinolones and tetracycline; hence, these two antibiotics can be used to treat infections caused by A. butzleri [12,13,14,15]. Notably, certain Arcobacter isolates have exhibited resistance to multiple drugs (known as multidrug resistance, MDR), including trimethoprim and cephalosporins, as previously reported [7, 16, 17]. Further investigations have revealed that an increase in the degree of antimicrobial resistance in A. butzleri strains may be attributed to their exposure to antibiotics commonly used in veterinary and public health applications as documented by Luangtongkum et al. [18].

To fully comprehend the pathogenicity of A. butzleri, it is important to investigate its virulence factors, such as adhesion, invasion, and cytotoxic capacity. However, current knowledge on these mechanisms and potential virulence factors remains insufficient [19,20,21]. Pathogens with certain virulence genes could potentially lead to infections in humans and other animals [22, 23]. Analysis of the A. butzleri RM 4018 genome has revealed the presence of several putative virulence genes, including cadF, cj1349, ciaB, mviN, pldA, tlyA, irgA, hecA, and hecB. The gene mviN encodes protein MViN, which is essential for peptidoglycan biosynthesis [24]. Similarly, in Campylobacter species, cadF and cj1349 encode fibronectin-binding proteins that enhance bacterial attachment to the host cells [25, 26]. Campylobacter invasive antigen B (ciaB) is involved in host cell invasion [27]. The gene hecA, a member of the filamentous hemagglutinin family, encodes protein HecA (Rojas), while hecB encodes haemolysin activation protein [28], and tlyA encodes hemolysin. Notably, hemolysin proteins are also present in other human infectious agents, such as Mycobacterium tuberculosis and Serpulina [29]. The irgA gene encodes an iron-regulated outer membrane protein (IrgA) that confers pathogenicity to E. coli [30, 31], while pldA encodes an outer membrane phospholipase PldA that mediates hemolysis of erythrocytes in humans [32]. However, it remains unclear whether these putative virulence factors function similarly in other bacterial pathogens.

DNA fingerprinting techniques, such as random amplification of polymorphic DNA (RAPD), restriction fragment length polymorphism (RFLP), pulsed-field gel electrophoresis (PFGE) and enterobacterial repetitive intergenic consensus-polymerase chain reaction (ERIC-PCR), are commonly used to study genotypic diversity of bacterial pathogens [33,34,35]. Of these methods, ERIC-PCR has a higher discriminatory power and can produce more consistent and intricate results than the other techniques [36, 37]. Moreover, ERIC-PCR is more effective and efficient method for strain-level genotyping. ERIC sequences, alternatively known as intergenic repetitive units, have been used in several studies to investigate the level of genetic variation and relatedness among pathogenic bacterial strains or species [36]. The intergenic repetitive units specifically correlated with the level of heterogeneity observed in Arcobacter species was likely the outcome of genetic recombination between progeny of parent genotypes [37]. Bacterial pathogens are known to be continuously evolving, and thus a high level of genetic diversity is expected in the case of A. butzleri.

Therefore, this study was initiated with the aim of investigating the degree of genetic diversity, prevalence of virulence genes, and status of antimicrobial resistance in A. butzleri strains isolated from surface water samples collected from different sites in a mixed-use but predominately agricultural watershed in eastern Ontario, Canada.


Detection and identification of Arcobacter species and A. butzleri by multiplex PCR analysis

Of the total 797 surface water samples, 2040 strains isolated from 11 sites with agricultural-dominated upstream land uses were confirmed as Arcobacter using genus-specific PCR assay. The sites with Arcobacter species isolated include Site 1 (n = 205; 10%), 5 (n = 645; 32%), 6 (n = 357; 18%), 9 (n = 514; 25%), 15 (n = 264; 13%), 17 (n = 1; 0.04%), 18 (n = 6; 0.3%), 19 (n = 12; 0.6%), 20 (n = 28; 1.4%), 22 (n = 2; 0.1%), and 23 (n = 6; 0.3%). Of these 2040 strains, 913 (48%) isolates from Sites 1 (n = 98; 11%), 5 (n = 279; 31%), 6 (n = 154; 17%), 9 (n = 247; 27%), 15 (n = 102; 11%), 17 (n = 1; 0.1%), 18 (n = 2; 0.2%), 19 (n = 5; 0.5%), 20 (n = 17; 2%), 22 (n = 2; 0.2%), and 23 (n = 6; 0.6%) were identified as A. butzleri by species-specific multiplex PCR assay. These A. butzleri isolates were subjected to characterization of Antimicrobial Resistance (AMR), Virulence-associated Genes (VAGs) detection, and genotyping using ERIC-PCR assays.

Antimicrobial resistance analysis

Of the total 913 A butzleri isolates, the majority were resistant to clindamycin (99%), followed by chloramphenicol (77%), nalidixic acid (49%), and azithromycin (48%) (Table 1). However, small proportion of the isolates displayed resistance to erythromycin (5%), ciprofloxacin (4%), gentamycin (2%) and tetracycline (1%). However, a significantly (p < 0.05) high level of resistance to clindamycin was observed in A. butzleri strains isolated from sites 6 (98%), 5 (97%), 15 (93%), 9 (90%) and 1 (88%). Similarly, sites 6 (53%), 5 (51%), and 9 (49%) showed high resistance to azithromycin. The frequency of isolates resistant to chloramphenicol was highest at sites 5 (86%) and 6 (82%) and relatively lower at sites 9 (73%), 15 (73%), and 1 (57%). Resistance against nalidixic acid was found most prevalent in isolates from sites 6 (58%) and 9 (54%), followed by sites 5 (47%) and 15 (45%). Overall, the detection frequency of antimicrobial-resistant isolates was significantly (p < 0.05) high at site 6, where a high level of resistance to all antimicrobial agents was observed. By contrast, many isolates from site 1 showed susceptibility to erythromycin, tetracycline, ciprofloxacin, and gentamycin. Sites 17 to 23, with a very limited number of isolates, showed 100% resistance to clindamycin, chloramphenicol, azithromycin and nalidixic acid. However, all isolates from sites 22 and 23 showed resistance to clindamycin and chloramphenicol as compared to isolates from sites 19 and 20 which showed (100%) resistance to clindamycin.

Table 1 Number (percentage) of A. butzleri isolates resistant to antimicrobial agents

Moreover, A. butzleri isolates (n = 786; 86%) showed significantly (p < 0.05) high multi-drug resistance (MDR), particularly in isolates from sites 5 (93%), 6 (88%), and 15 (86%) as compared to sites 9 (82%) and 1 (72%). Whereas 61% (n = 553) strains showed resistance to four or more antibiotics, while 33% (n = 303) isolates showed resistance to nalidixic acid and one of the following antimicrobials, azithromycin, chloramphenicol, or clindamycin. However, 1% of the isolates from sites 6 and 9 showed resistance to all eight antimicrobials tested. In addition, 4% (n = 37) of the total isolates from all sites were susceptible to all eight antimicrobials tested.

Detection and distribution of VAGs

Overall, all eight VAGs were detected in 913 A. butzleri strains at a variable frequency (Table 2), with ciaB (89%) and tlyA (88%) being most prevalent, followed by pldA (87%) and mviN (86%) genes. However, cj1349 (12%), hecA (15%), hecB (19%) and irgA (25%) genes were detected at relatively low frequency. Among eight VAGS, ciaB, mviN, tlyA and pldA genes were detected significantly (p < 0.05) high across all 11 sites. Whereas no significant (p > 0.05) difference, in the rate of prevalence of cj1349, hecB, hecA and irgA VAGs, across sites was observed. Among these eight VAGs, tlyA (94%) and mviN (90%) were detected in isolates from site 1 as compared to the ciaB (94%) and pldA (89%) genes detected in isolates from site 15. Overall, isolates from site 15 had the highest (> 88%) prevalence of ciaB, mviN, tlyA and pldA genes. On the other hand, a high prevalence of hecB and irgA genes were detected in site 6, as compared to sites 1 and 9 isolates that showed a high prevalence of cj1349 and hecA genes. Although sites 17 to 23 had limited isolates, ciaB, tlyA, mviN, hecB, pldA and irgA genes were detected in site 17 isolates as compared to sites 18 and 19 where all isolates were also positive for ciaB, tlyA, mviN and pldA genes. Whereas, all isolates from site 20 were positive for tlyA and mviN as compared to site 22, positive for tlyA, mviN and pldA genes. Interestingly, only ciaB gene was detected in all isolates from site 23 (Table 2).

Table 2 Number (percent) VAGs detected in A. butzleri isolates

For analyzing the occurrence of multi-VAGs, 889 (97%) A. butzleri strains were detected positive with a significantly high prevalence (p < 0.05) in sites 1 and 6 (99%) than sites 9 and 15 (98%) and 5 (96%), respectively. The majority of isolates (39%) had four or more VAGs with the combination of ciaB, tlyA, mviN and pldA genes. Interestingly, all eight VAGs were harbored only in two of the isolates from sites 9 and 5. Overall, 14 of the isolates, mostly from sites 1 and 5, were detected without any VAGs.

Co-occurrence of AMR and VAGs in A. butzleri strains

Overall, at a significance level of 0.05, A. butzleri strains from different sites that had resistance to clindamycin, chloramphenicol, nalidixic acid, and azithromycin were also detected positive for ciaB, tlyA, mviN and pldA genes. Whereas isolates that showed resistance to clindamycin were also positive for tlyA and mviN genes. A. butzleri isolates from sites 6 and 15 that showed resistance towards nalidixic acid were also resistant to azithromycin also detected positive for ciaB and thyA genes. The strains resistant to clindamycin and chloramphenicol were also positively correlated to the tlyA and mviN genes. Interestingly, mviN and hecA did not show co-occurrence patterns with each other in the tested isolates. On the other hand, mviN showed a co-occurrence pattern with pldA; however, cj1349 gene did not show a co-occurrence pattern with them. Isolates from sites 5 and 9 positive for ciaB and tlyA were frequently detected with mviN and pldA genes. Despite the site of isolation of A. butzleri, the isolates showed a high frequency of resistance to the clindamycin, chloramphenicol, nalidixic acid, and azithromycin were also positive for ciaB, tlyA, mviN and pldA genes, respectively.

ERIC-PCR profiling of A. butzleri strains

The ERIC-PCR fingerprinting pattern of A. butzleri strains from different sites were examined using cluster analysis. Based on theERIC-PCR band patterns (Fig. 1), with a discriminatory index 0.894, a total of 913 A. butzleri genotypes were distinguished. These isolates were grouped into 38 clusters using a 30% similarity cut-off. Overall, a total of 328 (36%) genotypes that were isolated from sites 1, 5, 6, 9 and 15 showed high similarity. Of these isolates, 201 derived from sites 5 and 9 were grouped with strains from site 15. Isolates from site 1 were grouped into 3 clusters, represented by 10 to 39 isolates, of which, 27 genotypes showed > 40% homogeneity. Isolates from site 5 formed 6 clusters; each was represented by 21 to 100 isolates. Of the six clusters, 68 genotypes had > 40% homogeneity. Isolates from site 6 were grouped into eight main clusters and each was represented by 5 to 35 isolates, with 45 genotypes showing > 20% homogeneity. At site 9, a total of eight clusters were observed and each was represented by 6 to 39 isolates, where 52 genotypes showed > 35% homogeneity. Moreover, isolates from site 15 formed five clusters and each was represented by 7 to 35 isolates. At site 15, 16 genotypes showed > 30% homogeneity. Interestingly, the highest similarity of genotypes was detected in isolates from sampling sites 5, 9 and 15, which indicates that the water contamination may occur from the same sources.

Fig. 1
figure 1

Dendrogram analysis showing ERIC-PCR based distinct band patterns of A. butzleri strains isolated from surface water samples collected from an agricultural watershed. * denotes various agricultural sampling sites

Co-occurrence of ERIC-PCR and AMR

Strains that showed similar ERIC-PCR profiles from site 1 exhibited remarkable resistance to clindamycin (94%) and chloramphenicol (85%), in contrast to lower resistance levels to azithromycin (73%) and nalidixic acid (42%), respectively. Similarly, isolates from site 5 were highly resistant to clindamycin (92%) and chloramphenicol (77%), with reduced resistance to azithromycin (66%) and nalidixic acid (50%). Strains from site 6 displayed the highest resistance to chloramphenicol (89%) and clindamycin (84%), as opposed to nalidixic acid (51%) and azithromycin (49%). Moreover, site 9 isolates showed high resistance to clindamycin (94%) and chloramphenicol (88%), compared to nalidixic acid (66%) and azithromycin (49%). Lastly, site 15 isolates showed the highest resistance to each of clindamycin and chloramphenicol (68%), followed by azithromycin (30%) and nalidixic acid (20%) (Table 3).

Table 3 Number (percentage) of isolates showing similar ERIC-PCR and AMR patterns

Conversly, isolates from all sites showed low to no resistance to ciprofloxacin, erythromycin, gentamycin and tetracycline. Isolates from sites 17, 18 and 20 exhibited dissimilar ERIC-PCR profiles, whereas sites 20 and 23 displayed 100% resistance to clindamycin and chloramphenicol. However, site 19 isolates showed only demonstrated 100% resistance to clindamycin alone. Isolates with ERIC-PCR patterns displaying greater than > 80% homogeneity in cluster analysis demonstrated the highest rates of MDR in isolates at sites 1 (79%) and 9 (70%), in contrast to site 5 (58%) and sites 6 and 15 (both at 55%). In contrast, strains from sites 17 to 23 showed no MDR. A notable presence of resistance to a combination of clindamycin, chloramphenicol, and azithromycin was observed in these isolates (Table 3).

Co-occurrence of ERIC-PCR pattern and VAGs

The ciaB gene was found most prevalent at sites 6 (96%) and 5 and 9 (both at 87%) compared to sites 1 (79%) and 15 (67%). Isolates from site 6 exhibited high prevalence of the mviN gene (81%), followed by sites 9 (75%), 5 (74%), 15 (58%) and 1 (38%). Conversly, tlyA was more prevalent at sites 9 (94%) and 6 (85%) compared to sites 5 (78%), 1 (73%) and 15 (64%). The pldA gene was most prevalent at site 6 (74%), followed by sites 9 (70%), 5 (54%), 15 (43%) and 1 (33%). The cj1349, irgA, hecB and hecA genes were present in fewer than 15% of isolates. Moreover, isolates with similar ERIC-PCR patterns showing greater than 80% homogeneity had a high prevalence of multi-VAGs (three or more) at sites 6 and 9 (59%), 1 (56%), followed by sites 15 (47%) and 5 (6%) isolates. However, no isolates with multi-VAGs were found at sites 17 to 23 (Table 4).

Table 4 Number (percentage) of isolates showing similar ERIC-PCR and VAGs patterns

Genotypically, of the 52 isolates from site 1 with similar ERIC-PCR patterns, 37 (71%) showed 100% similarity in cluster analysis, and of those, 5 (10%) isolates exhibited similar AMR and VAGs profiles. From site 5, 98 of 132 isolates (74%) showed 100% similarity, with 19 (14%) similar AMR and VAGs profiles. At site 6, 51 of 80 isolates (64%) had 100% similarity with 7 (9%) isolates demonstrated similar AMR and VAGs profiles. At site 9, 103 of 140 isolates (74%) showed 100% similarity with 21 (15%) showing similar AMR and VAGs profiles. For site 15 isolates, 63 (83%) of 76 displayed 100% homogeneity, and 11 (15%) had similar AMR and VAGs profiles. Sites 17,18 and 22, displayed dissimilar ERIC-PCR patterns, whereas sites 18, 20 and 23 had similar ERIC-PCR patterns but did not achieve 100% homogeneity in cluster analysis (Table 5).

Table 5 Number (percent) of A. butzleri isolates showing similar ERIC-PCR, AMR and VAGs patterns in cluster analysis

Interestingly, seven isolates from site 9 showed high similarity to four isolates from site 15, showed 100% similarity cluster analysis and had similar AMR and VAGs profiles. Moreover, five isolates from site 1 were found 100% similar to three isolates from site 6 in cluster analysis and exhibited similar AMR and VAGs profiles, where these isolates also exhibited MDR and multi-VAGs profiles.


A high number (63%) of human infections are due to the consumption of or close contact with A. butzleri contaminated water [1, 38]. In the present study, 45% of the total isolates from agricultural surface water were identified as A. butzleri. Of eleven stream sites, sites 1, 5, 6, 9 and 15 demonstrated the highest prevalence of A. butzleri. Particularly, site 15 is located in a smaller stream order and possibly impacted by drainages from septic systems of nearby houses and livestock in close proximity to water ways. The other sites, located in medium to large streams, might have additional contamination sources such as wildlife (migratory birds) and urban runoff. Although the rate of isolation of Arcobacter species varies significantly globally, 20.8% of water samples were found positive for Arcobacter species with 90% identified as A. butzleri, as reported by Laishram et al. [37]. Another study showed a relatively lower prevalence of A. butzleri in creeks and streams, with no detection in ponds and drinking water samples [10]. However, in the current study, a moderate (45%) frequency was observed. The isolation rate can be influenced by factors such as sampling and isolation methods, ecological and geographical characteristics of sampling sites, and water pollution mechanisms [39,40,41]. Many studies have suggested that the abundance of A. butzleri in water resources is due to their high viability in water and a potential competitive inhibitory effect on the other species [41]. The presence of Arcobacter species in water resources signals a possible risk of transmission to humans and animals, as well as the potential for contamination of food products.

The phenotypic characterization, including antimicrobial resistance of A. butzleri isolates, demonstrates that isolates from these 11 sites have developed resistance to clindamycin (99%), chloramphenicol (77.7%), azithromycin (48.8%) and nalidixic acid (47.4%). AMR patterns observed were comparable with a previous study where A. butzleri was isolated from different environmental aquatic sources [12]. However, the acquired resistance to clindamycin is particularly concerning because this antibiotic is the first-line drug for treating Campylobacter infection in humans [13, 15]. Similarly, azithromycin and ciprofloxacin are widely used for human infections [13]. Notably, all strains in our study were susceptible to tetracycline (1%), gentamycin (2%), ciprofloxacin (4%), and erythromycin (5%). Erythromycin and gentamicin are also frequently used as alternative treatments for campylobacteriosis [13, 15]. Thus, ciprofloxacin, tetracycline, erythromycin, and gentamycin are recommended as the first-choice antibacterial agents, with erythromycin as a close second for the treatment of Arcobacter-associated infections in humans and animals. A plausible explanation for the emergence and spread of clinical antibiotic resistance found in the isolates could be attributed to the unrestricted usage of these antibiotics in agriculture, as suggested by Chang et al [42]. This study supports the notion that ciprofloxacin, tetracycline, gentamycin, and erythromycin could be the antibiotics of choice for treating other pathogens in agriculture and animal husbandry practices where antimicrobial resistance is prevalent [42].

In our study, a high frequency (33%) of resistance to a combination of four tested antimicrobials (azithromycin, chloramphenicol, clindamycin, and nalidixic acid) was observed, with 1% of the A. butzleri population resistant to all eight antimicrobials tested. Previous studies have reported a variable frequency of MDR (ranging from 20 to 72%) among A. butzleri strains isolated from human and animal sources [15, 38, 43,44,45]. Wastewater from industrial plants, healthcare services, and agriculture are point sources for antimicrobials, antibiotic-resistant bacteria, and antibiotic resistance genes [46, 47]. The presence of A. butzleri in various water sites suggests the role of these isolates in the environment and food chain contamination, posing a potential threat to public health.

In various studies, Arcobacter VAGs have been detected in humans, animals, and food samples [17, 19,20,21,22, 37, 48, 49]. Information on the virulence properties of potentially pathogenic A. butzleri remains scarce, but some studies have elucidated the distribution of virulence markers and the adhesive, invasive and toxicity capacity of arcobacters, highlighting its high invasion ability to the host cells [21]. Although several genes have already been identified in the A. butzleri genome, it is still unclear if the putative virulence factors have similar functions in other microbial species [28]. The eight VAGs in A. butzleri isolates were detected at a variable frequency; four of these, ciaB, mviN, tlyA and pldA were predominantly (> 86%) observed. Similar results were reported by Douidah et al. [19], which showed a 93% prevalence for the ciaB and mviN VAGs. However, low prevalence (< 25%) was noted for cj1349, hecA, hecB, and irgA VAGs across all tested isolates, consistent with previous studies [19, 20]. However, Rathlavath et al. [50] demonstrated a high (97%) prevalence of cj1349 gene in A. butzleri strains isolated from coastal water, fish and shellfish samples, which contradicts to our study results. This discrepancy could be due to the origin of the isolates. The high prevalence of ciaB, tlyA, mviN and pldA VAGs in all isolates also demonstrates a correlation between the presence of these VAGs and infections in humans and animals. Overall, it can be concluded that the common distribution pattern of VAGs among A. butzleri isolates collected from various water sites may have the same pathogenic potential as known for strains isolated from various food and fecal sources. Detailed data analysis demonstrated the prevalence of multiple VAGs (ranging from 2 to 8) in most isolates (97%), with a considerable percent (54%) having four or more VAGs. These results are in accordance with an earlier study by Rathlavath et al. [50]. While research on the roles of these VAGs in A. butzleri continues, the data from this study could be useful in understanding their pathogenic potential and associated health risks to humans and animals.

Since the relationship between resistance and virulence in arcobacters is still unclear, it is often presumed that MDR arcobacters with more antimicrobial drug resistance account for greater virulence [51]. This could be a possible misconception, as it greatly relies on the phylogenetic lineage and resistant determinants of the specific isolate. This study is the first to demonstrate a correlation between AMR and VAGs. Therefore, the association and co-occurrence of AMR and VAGs in A. butzleri strains isolated from various agricultural sites were critically analyzed. Exploring this correlation in other studies could be helpful in managing diseases associated with A. butzleri. The co-occurrence analysis revealed a positive relationship between resistance to azithromycin, chloramphenicol, clindamycin, and nalidixic acid and the presence of ciaB, mviN, tlyA and pldA VAGs. In a previous study, a similar positive correlation between the presence of AMR and VAGs in E. coli and E. faecalis as well as S. aureus and S. pneumoniae were reported [52,53,54]. However, a negative correlation was observed for the presence of tetracycline and hecA VAG in our isolates, in sites 17, 18, 19, 20, 22, and 23 that had fewer isolates but a high prevalence of AMR and VAGs. Further studies are needed to elucidate the epidemiology of these virulence factors in the pathogenesis of A. butzleri.

In previous studies, ERIC-PCR has been used to determine the molecular etiology of arcobacters isolated from various sources such as poultry [6, 37], cattle [55] and food [36]. ERIC-PCR based fingerprinting was used to genotype A. butzleri strains isolated from different agricultural sites. Dendrogram analysis and clustering pattern revealed homogeneity among some of the A. butzleri isolates collected from sites 1, 9 and 15. However, some of the A. butzleri isolates from sites 15 and 23 did not cluster with other isolates, indicating the possibility of heterogeneity among these isolates. High genetic similarity was seen among the isolates from the same site. The results of this study showed a significant genetic similarity among isolates, highlighting a common source of contamination. These results are in congruency with previously reported data where low variability was observed in A. butzleri isolates from a dairy plant [56]. The analysis suggests that the water at these sites could be fecally contaminated [57] supported by previous studies linking the presence of Arcobacter species in water to high levels of fecal contamination [58, 59].

Seven isolates from site 9 and four from site 15 displayed identical ERIC-PCR patterns and 100% homogeneity in cluster analysis, with similar AMR and VAGs profiles. Similarly, five isolates from site 1 showed 100% homogeneity to three isolates from site 6, showing similar AMR and VAGs profiles. These results reveal a close genetic relationship between isolates from sites 9 and 15, and sites 1 and 6, suggesting a common origin. Although site 9 is hydrologically disconnected from site 15 and site 6 has water that could feed site 1, over long distances and under significant dilution, inputs from other sources, including runoff and wildlife (migratory birds), may contribute to the prevalence of A. butzleri strains in these sites. All isolates with similar ERIC-PCR patterns, cluster analysis homogeneity, and similar AMR and VAGs profiles were potentially have a shared origin and possibly more pathogenic or resistant. This study concludes that the presence of AMR and VAGs may increase the resistance and pathogenicity of strains, posing potential health risks to humans and animals.


In conclusion, this study has provided insightes into the prevalence, distribution, and co-occurrence of AMR and VAGs among A. butzleri strains isolated from eleven stream sites within a mixed-use but predominately agricultural river basin in eastern Ontario Canada. The high frequency of MDR in isolates from different sites indicates potential public health risks. Moreover, susceptibility to ciprofloxacin, gentamycin, erythromycin, and tetracycline can be considered drugs of choice in agriculture and animal husbandry. The lack of standardized procedures to determine antibiotic susceptibility patterns underscores the urgent need for developing such testing methods. The majority of the isolates were detected with ciaB, mviN, tlyA and pldA VAGs, necessitating further investigation into their pathogenic roles. A positive correlation between AMR and VAGs in A. butzleri isolates raises concerns about increased resistance and pathogenicity. Furthermore, given the pathogenic potential of these isolates, contaminated water at these sites may serve as a plausible transmission route of virulent strains to humans and animals.


Water sampling, isolation, and culturing

A total number of 797 surface water samples were collected from 11 different water sampling sites (sites designated as 1, 5, 6, 9, 15, 17–20, 22 and 23) located in the South Nation River basin near Ottawa, Ontario, Canada (Fig. 2). A detail watershed description (Table 6) is previously given by Wilkes et al. [57, 60]. The sample collection was carried out on a bi-weekly basis between May to November 2012 to 2019. The samples were collected in sterile bottles, placed on ice in a cooler, transported to the AAFC-Ottawa laboratory, and processed within 24 h of their collection for microbiological analyses. For the isolation of Arcobacter, the samples were processed according to the previously described procedure [61]. Briefly, 1 mL of water sample was resuspended in 9 mL of peptone water (PW) using a single ten-fold serial (ranging from 10–1 to 10–6) dilution. A 100 µL of each sample was plated on Arcobacter selective isolation agar (ASIA) media (Oxoid, Nepean, ON) containing antibiotics supplements; fluorouracil, amphotericin-B, cefoperazone, novobiocin, and trimethoprim. The plates were incubated for three to six days at 30 °C under microaerophilic (85% N2, 10% CO2 and 5% O2) conditions.

Fig. 2
figure 2

Map showing agricultural watershed surface water sampling sites selected for the study

Table 6 Land use description of each site sampled for isolation of A. butzleri in this study

Based on colony morphology, a putative Arcobacter single colony was sub-cultured on modified Agarised Arcobacter Medium (m-AAM) containing selective antibiotic supplements; cefoperazone, amphotericin-B, and teicoplanin. The plates were incubated according to the conditions mentioned above. The putative culture isolates were further confirmed by Gram staining reaction and PCR assays designed specifically for the Arcobacter genus and A. butzleri species, respectively.

For this study, A. butzleri ATCC 49616, A. cryaerophilus NCTC 11885, and A. skirrowii ATCC 51132 reference strains were used as positive and negative controls. The strains were also grown on selective growth media according to the specified culture conditions mentioned above.

Nucleic acid extraction

The DNA extraction from Arcobacter culture isolates, recovered from water samples, were prepared by re-suspending a purified single colony in a sterile 1.5 mL microfuge tube containing 100 μL TE buffer (10 mM Tris–HCl and 1 mM EDTA, pH 8.0). The cells were gently mixed and boiled for 10 min, as described by Khan et al. [62]. The tube was centrifuged at high speed for 60 s., and the supernatant containing purified DNA was used for subsequent examination and quantification by agarose gel electrophoresis and an ND-1000 spectrophotometer using low-range quantitative DNA marker (Fisher Scientific, Ottawa, ON), respectively. The DNA extract was stored at -20 °C for further PCR analyses.

Arcobacter genus and species-specific multiplex PCR assays

To identify Arcobacter isolates, a DNA-based PCR amplification assay was performed using Arcobacter genus-specific oligonucleotide primers and PCR protocol previously described by Harmon and Wesley [63]. The PCR reaction was performed in a 25 µL volume containing 5.0 ng of template DNA, 50 pmol each of Arcobacter 16S rRNA gene primer pair (Table 7), 1.25 U of Taq DNA polymerase 1 × buffer with MgCl2, 200 μM each of the dNTPs (Fisher Scientific, Nepean, ON). The PCR reaction was performed in a Mastercycler Gradient PCR system (Eppendorf, Hauppauge, NY) using an initial denaturation (94 °C for 4 min), followed by 25 amplification cycles consisted of denaturation: 94 °C for 60 s; annealing: 56 °C for 60 s.; extension: 72 °C for 60 s. with final extension cycle at 72 °C for 7 min.

Table 7 List of Arcobacter genus and species-species oligonucleotide PCR primer sequences used in this study

Furthermore, A. butzleri, A. cryaerophilus, A. cibarius, A. faecis, A. lanthieri, and A. skirrowii were identified using species-specific multiplex PCR protocol described by Khan et al. [62]. The PCR amplification reaction was carried out in a 25 μL reaction mixture containing 1 × buffer with MgCl2, 200 μmol L−1 each of the dNTPs (Fisher Scientific, Nepean, ON, Canada), 1 U of Ex-Taq DNA polymerase, 0.1 μmol L−1 of A. butzleri, 0.2 μmol L−1 of A. lanthieri and A. faecis, 0.3 μmol L−1 of A. skirrowii and A. cryaerophilus and 0.4 μmol L−1 of A. cibarius primer pairs (Table 7) and 10 ng of each target DNA template. The PCR protocol was performed with an initial denaturation (94 °C for 3 min) followed by 35 cycles consisted of denaturation (94 °C for 30 s), annealing (58 °C for 30 s), and extension (72 °C for 30 s) with a final extension at 72 °C for 5 min.

The amplified products were electrophoresed on a 1.5% agarose gel and stained with ethidium bromide (0.5 µg mL−1), visualized under an ultraviolet (UV) transilluminator, and photographed using an Alpha Imager (Fisher Scientific) gel documentation system.

Antimicrobial resistance analysis

Antimicrobial susceptibility to eight antibiotics (three concentrations of each), including azithromycin (Azi; 4, 8, 16 μg mL−1;), chloramphenicol (Chl; 16, 32, 64 μg mL−1), ciprofloxacin (Cip; 2, 4, 8 μg mL−1), clindamycin (Cli; 4, 8, 16 μg mL−1), erythromycin (Ery; 16, 32, 64 μg mL−1), gentamycin (Gen; 4, 8, 16 μg mL−1), nalidixic acid (NaI; 32, 64, 128 μg mL−1), and tetracycline (Tet; 8, 16, 32 μg mL−1) obtained from Sigma-Aldrich (ON, Canada) was determined based on the criteria described by CIPARS and standards of the Clinical and Laboratory Standards Institute [64] at breakpoint concentration commonly used for Campylobacter species since no breakpoint values for Arcobacter spp. are available [65]. The agar dilution method, as described by Gaudreau and Gilbert [66], was applied. Briefly, A. butzleri strains were incubated overnight in a 96-well microplate with 200 μL per well of Arcobacter medium (AM) broth at 30 °C under microaerophilic conditions. The cells were then transferred to the surface of rectangular Mueller–Hinton (MH) agar plates by a 96-floating pin replicator (V&P Scientific, San Diego, CA). Agar plates were incubated at 30 ºC under microaerophilic conditions for 3 days, and the growth of A. butzleri isolates on plates with antimicrobials was compared to their growth on control plates without antimicrobial agents. The results were interpreted in accordance with the CLSI [64].

Multiplex PCR assays for detection of virulence-associated genes (VAGs)

A. butzleri culture isolates were analyzed for the presence of VAGs using our previously developed three multiplex PCR-based assays. Each of the assays was designed to detect: 1). ciaB and cj1349, 2). pldA, irgA and hecA, or 3). tlyA, mviN and hecB genes, respectively [67]. Each multiplex PCR (mPCR) amplification assay was carried out in 25 μL of reaction mixtures containing 1 μL (50–70 ng μL−1) of template DNA, 1 U of Ex Taq DNA polymerase (Fisher Scientific), and the compatible PCR reagents, including 1 × buffer with MgCl2, 200 µM each of the dNTPs, and 0.1 µM of each set of forward and reverse primer pair (Table 8). A total volume of 25 µL was adjusted by adding nuclease-free water to the reaction mixture. All three PCR assays were carried out in a Mastercycler Gradient PCR system (Eppendorf, Hauppauge, NY) under the following standardized cycling conditions: initial denaturation at 95 °C for 4 min, 30 cycles of denaturation at 95 °C for 30 s, annealing at 56 °C for 45 s, extension at 72 °C for 45 s, and final elongation at 72 °C for 5 min. Amplified PCR products were analyzed by gel electrophoresis in 2% agarose gel using 100 bp DNA size marker (Life Technologies, Grand Island, NY). The gels were stained, scanned and photographed using the procedure mentioned in the preceding section.

Table 8 List of virulence-associated genes, oligonucleotide PCR primer sequences and amplicon sizes for each multiplex PCR assay used in this study [67]

Enterobacterial repetitive intergenic consensus-polymerase chain reaction (ERIC-PCR) assay

To assess the genetic diversity of A. butzleri isolates, ERIC-PCR assay was carried out using method described by Houf et al. [68]. Briefly, the PCR was performed in a mixture containing a total volume of 25 mL containing 5 mL 5 × PCR buffer, 4 mM MgCl2, 10 pmol uL−1 primers ERIC-F2 (5´-ATG TAA GCT CCT GGG GAT TCA C-3´) and ERIC-R1 (5´-AAG TAA GTG ACT GGG GTG AGC G-3´), 5U Taq polymerase and 10 uL (2 ng) of DNA template. DNA amplification consisted of an initial denaturation at 94 °C for 2.5 min, followed by 40 cycles of amplification (denaturation at 94 ºC for 1 min, annealing at 25 ºC for 1 min and extension at 72 ºC for 2 min) ending with final elongation at 72 ºC for 7 min. The amplified products were analyzed by electrophoresis in a 2% (w/v) 1xTAE agarose gel at 100 V for 2.5 h. The gels were stained, scanned and photographed using the procedure mentioned in the preceding section.

The ERIC-PCR patterns were analyzed by constructing dendrograms using the Gel Compar-II program (Applied Maths BVBA, Kortrijk, Belgium) where all isolates obtained in this study were grouped into different clusters and clades.

Statistical analysis

Statistical analyses, including McNemar Chi-square Contingency and Fisher's Exact Tests, were applied to compare the rate of prevalence, antimicrobial resistance (AMR), virulence-associated genes (VAGs) and genotypic properties in A. butzleri strains across various sampling sites and to determine the pairwise association and co-occurrence of the AMR, VAGs, and genetic determinants in A. butzleri strains using STATISTICA 10.0 [69]. A p-value of < 0.05 was considered statistically significant. The discriminatory index of ERIC-PCR was calculated based on Simpson's Index of Diversity as described by Hunter and Gaston [70].

Availability of data and materials

The material and dataset used and analyzed in this study are added to the manuscript. There is no additional data available from the authors.



Antimicrobial Resistance


Virulence-associated genes


Enterobacterial Repetitive Intergenic Consensus-Polymerase Chain Reaction


Multidrug Resistance


Peptone Water


Modified Agarised Arcobacter Medium


  1. Collado L, Figueras MJ. Taxonomy, epidemiology, and clinical relevance of the genus Arcobacter. Clin Microbiol Rev. 2011;24:174–92.

  2. Prouzet-Mauléon V, Labadi L, Bouges N, Menard A, Megraud F. Arcobacter butzleri: underestimated enteropathogen. Emerg Infect Dis. 2006;12:307.

  3. Vandenberg O, Dediste A, Houf K, Ibekwem S, Souayah H, Cadranel S, Douat N, Zissis G, Butzler JP, Vandamme P. Arcobacter species in humans. Emerg Infect Dis. 2004;10(10):1863–7.

    Article  PubMed Central  Google Scholar 

  4. Yan J-J, Ko W-C, Huang A-H, Chen H-M, Jin Y-T, Wu J-J. Arcobacter butzleri bacteremia in a patient with liver cirrhosis. J Form Med Assoc. 2000;99:166–9.

    CAS  Google Scholar 

  5. International Commission on Microbiological Specifications of Foods (ICMSF): Microorganisms in Foods 7. Microbiological Testing in Food Safety Management. New York: Kluwer Academic/Plenum Publishers; 2002.

    Book  Google Scholar 

  6. Atabay HI, Aydin F. Susceptibility of Arcobacter butzleri isolates to 23 antimicrobial agents. Lett Appl Microbiol. 2001;33:430–3.

  7. Ho HTK, Lipman LJA, Gaastra W. Arcobacter, what is known and unknown about a potential foodborne zoonotic agent! Vet Microbiol. 2006;115:1–13.

  8. Hume ME, Harvey RB, Stanker LH, Droleskey RE, Poole TL, Zhang HB. Genotypic variation among Arcobacter isolates from a farrow-to-finish swine facility. J Food Protect. 2001;64:645–51.

  9. Vicente-Martins S, Oleastro M, Domingues FC, Ferreira S. Arcobacter spp. at retail food from Portugal: Prevalence, genotyping and antibiotics resistance. Food Cont. 2018;85:107–12.

    Article  CAS  Google Scholar 

  10. Çelik E, Ünver A. Isolation and identification of Arcobacter spp. by multiplex PCR from water sources in Kars Region. Curr Microbiol. 2015;71(5):546–50.

    Article  CAS  PubMed  Google Scholar 

  11. Alonso R, Girbau C, Martinez-Malaxetxebarria I, Fernández-Astorga A. Multilocus sequence typing reveals genetic diversity of foodborne Arcobacter butzleri isolates in the North of Spain. Int J Food Microbiol. 2014;191:125–8.

    Article  CAS  PubMed  Google Scholar 

  12. Fera MT, Maugeri TL, Giannone M, Gugliandolo C, La Camera E, Blandino G, Carbone M. In vitro susceptibility of Arcobacter butzleri and Arcobacter cryaerophilus to different antimicrobial agents. Int J Antimicrob Agents. 2003;21:488–91.

  13. Houf K, Devriese LA, Haesebrouck F, Vandenberg O, Butzler JP, van Hoof J, Vandamme P. Antimicrobial susceptibility patterns of Arcobacter butzleri and Arcobacter cryaerophilus strains isolated from humans and broilers. Microb Drug Resist. 2004;10(3):243–7.

    Article  CAS  PubMed  Google Scholar 

  14. Son I, Englen MD, Berrang ME, Fedorka-Cray PJ, Harrison MA. Antimicrobial resistance of Arcobacter and Campylobacter from broiler carcasses. Int J Antimicrob Agents. 2007;29:451–5.

  15. Vandenberg O, Houf K, Douat N, Vlaes L, Retore P, Butzler J-P, Dediste A. Antimicrobial susceptibility of clinical isolates of non-jejuni/coli campylobacters and arcobacters from Belgium. J Antimicrob Chemother. 2006;57:908–13.

  16. Ferreira S, Queiroz JA, Oleastro M, Domingues FC. Insights in the pathogenesis and resistance of Arcobacter: A review. Crit Rev Microbiol. 2016;42:364–83.

  17. Girbau C, Guerra C, Martínez-Malaxetxebarria I, Alonso R, Fernández-Astorga A. Prevalence of ten putative virulence genes in the emerging foodborne pathogen Arcobacter isolated from food products. Food Microbiol. 2015;52:146–9.

  18. Luangtongkum T, Jeon B, Han J, Plummer P, Logue CM, Zhang Q. Antibiotic resistance in Campylobacter: emergence, transmission and persistence. Future Microbiol. 2009;4:189–200.

  19. Douidah L, De Zutter L, Baré J, De Vos P, Vandamme P, Vandenberg O, Van den Abeele A-M, Houf K. Occurrence of putative virulence genes in Arcobacter species isolated from humans and animals. J Clin Microbiol. 2012;50:735–41.

  20. Karadas G, Sharbati S, Hänel I, Messelhäußer U, Glocker E, Alter T, Gölz G. Presence of virulence genes, adhesion and invasion of Arcobacter butzleri. J Appl Microbiol. 2013;115:583–90.

  21. Levican A, Alkeskas A, Günter C, Forsythe SJ, Figueras MJ. Adherence to and invasion of human intestinal cells by Arcobacter species and their virulence genotypes. Appl Environ Microbiol. 2013;79:4951–7.

  22. Ferreira S, Queiroz JA, Oleastro M, Domingues FC. Genotypic and phenotypic features of Arcobacter butzleri pathogenicity. Microb Pathog. 2014;76:19–25.

  23. Levican A, Collado L, Yustes C, Aguilar C, Figueras MJ. Higher water temperature and incubation under aerobic and microaerobic conditions increase the recovery and diversity of Arcobacter spp. from shellfish. Appl Environ Microbiol. 2014;80(1):385–91.

  24. Inoue A, Murata Y, Takahashi H, Tsuji N, Fujisaki S, Kato J-i. Involvement of an essential gene, mviN, in murein synthesis in Escherichia coli. J Bacteriol. 2008;190:7298–301.

  25. Flanagan RC, Neal-McKinney JM, Dhillon AS, Miller WG, Konkel ME. Examination of Campylobacter jejuni putative adhesins leads to the identification of a new protein, designated FlpA, required for chicken colonization. Infect Immun. 2009;77:2399–407.

  26. Monteville MR, Yoon JE, Konkel ME. Maximal adherence and invasion of INT 407 cells by Campylobacter jejuni requires the CadF outer-membrane protein and microfilament reorganization. Microbiol. 2003;149:153–65.

  27. Konkel ME, Kim BJ, Rivera-Amill V, Garvis SG. Identification of proteins required for the internalization of Campylobacter jejuni into cultured mammalian cells. Adv Exp Med Biol. 1999;473:215–24.

    Article  CAS  PubMed  Google Scholar 

  28. Miller WG, Parker CT, Rubenfield M, Mendz GL, Wösten MM, Ussery DW, Stolz JF, Binnewies TT, Hallin PF, Wang G. The complete genome sequence and analysis of the epsilonproteobacterium Arcobacter butzleri. PLoS ONE. 2007;2:e1358.

  29. Sałamaszyńska-Guz A, Klimuszko D. Functional analysis of the Campylobacter jejuni cj0183 and cj0588 genes. Curr Microbiol. 2008;56:592–6.

  30. Johnson JR, Jelacic S, Schoening LM, Clabots C, Shaikh N, Mobley HL, Tarr PI. The IrgA homologue adhesin Iha is an Escherichia coli virulence factor in murine urinary tract infection. Infect Immun. 2005;73:965–71.

  31. Mey AR, Wyckoff EE, Oglesby AG, Rab E, Taylor RK, Payne SM. Identification of the Vibrio cholerae enterobactin receptors VctA and IrgA: IrgA is not required for virulence. Infect Immun. 2002;70:3419–26.

  32. Grant KA, Belandia IU, Dekker N, Richardson PT, Park SF. Molecular characterization of pldA, the structural gene for a phospholipase A from Campylobacter coli, and its contribution to cell-associated hemolysis. Infect Immun. 1997;65:1172–80.

  33. Bagalakote PS, Rathore RS, Ramees TP, Mohan HV, Sumankumar M, Agarwal RK, Kumar A, Dhama K. Molecular characterization of Arcobacter isolates using Randomly Amplified Polymorphic DNA-Polymerase Chain Reaction (RAPD-PCR). Asian J Ani Vet Advan. 2014;9:543–55. doi=ajava.2014.543.555.

  34. de Las RB, Marcobal A, Muñoz R. Allelic diversity and population structure in Oenococcus oeni as determined from sequence analysis of housekeeping genes. Appl Environ Microbiol. 2004;70(12):7210–9.

    Article  CAS  Google Scholar 

  35. Merga JY, Leatherbarrow AJ, Winstanley C, Bennett M, Hart CA, Miller WG, Williams NJ. Comparison of Arcobacter isolation methods, and diversity of Arcobacter spp. in Cheshire United Kingdom. Appl Environ Microbiol. 2011;77(5):1646–50.

  36. Aydin F, Gümüşsoy KS, Atabay HI, Iça T, Abay S. Prevalence and distribution of Arcobacter species in various sources in Turkey and molecular analysis of isolated strains by ERIC-PCR. J Appl Microbiol. 2006;103:27–35.

  37. Laishram M, Rathlavath S, Lekshmi M, Kumar S, Nayak BB. Isolation and characterization of Arcobacter spp. from fresh seafood and the aquatic environment. Int J Food Microbiol. 2016;2(232):87–9.

    Article  CAS  Google Scholar 

  38. Shah AH, Saleha AA, Zunita Z, Murugaiyah M. Arcobacter – An emerging threat to animals and animal origin food products? Trends Food Sci Technol. 2011;22(5):225–36.

    Article  CAS  Google Scholar 

  39. Wilkes G, Brassard J, Edge TA, Gannon V, Gottschall N, Jokinen CC, Jones TH, Khan IU, Marti R, Sunohara MD, Topp E, Lapen DR. Long-term monitoring of waterborne pathogens and microbial source tracking markers in paired agricultural watersheds under controlled and conventional tile drainage management. Appl Environ Microbiol. 2014;80(12):3708–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Wilkes G, Sunohara MD, Topp E, Gottschall N, Craiovan E, Frey SK, Lapen DR. Do reductions in agricultural field drainage during the growing season impact bacterial densities and loads in small tile-fed watersheds? Water Res. 2018;151:423–38.

    Article  CAS  PubMed  Google Scholar 

  41. Talay F, Molva C, Atabay HI. Isolation and identification of Arcobacter species from environmental and drinking water samples. Folia Microbiol. 2016;61:479–84.

    Article  CAS  Google Scholar 

  42. Chang Q, Wang W, Regev G, Lipsitch M, Hanage W. Antibiotics in agriculture and the risk to human health: how worried should we be? Evol Appl. 2014;8(3):240–7.

  43. Abay S, Kayman T, Hizlisoy H, Aydin F. In vitro antibacterial susceptibility of Arcobacter butzleri isolated from different sources. J Vet Med Sci. 2012;74:613–6.

  44. Harrass B, Schwarz S, Wenzel S. Identification and characterization of Arcobacter isolates from broilers by biochemical tests, antimicrobial resistance patterns and plasmid analysis. J Vet Med Ser B. 1998;45:87–94.

  45. Kabeya H, Maruyama S, Morita Y, Ohsuga T, Ozawa S, Kobayashi Y, Abe M, Katsube Y, Mikami T. Prevalence of Arcobacter species in retail meats and antimicrobial susceptibility of the isolates in Japan. Int J Food Microbiol. 2004;90:303–8.

  46. Rizzo L, Manaia C, Merlin C, Schwartz T, Dagot C, Ploy MC, Michael I, Fatta-Kassinos D. Urban wastewater treatment plants as hotspots for antibiotic resistant bacteria and genes spread into the environment: a review. Sci Total Environ. 2013;1(447):345–60.

    Article  CAS  Google Scholar 

  47. Serwecińska L. Antimicrobials and antibiotic-resistant bacteria: a risk to the environment and to public health. Water. 2020;12:3313.

    Article  CAS  Google Scholar 

  48. Tabatabaei M, Aski HS, Shayegh H, Khoshbakht R. Occurrence of six virulence-associated genes in Arcobacter species isolated from various sources in Shiraz, Southern Iran. Microb Pathog. 2014;66:1–4.

  49. Zacharow I, Bystroń J, Wałecka-Zacharska E, Podkowik M, Bania J. Prevalence and antimicrobial resistance of Arcobacter butzleri and Arcobacter cryaerophilus isolates from retail meat in Lower Silesia region. Poland Pol J Vet Sci. 2015;18:63–9.

  50. Rathlavath S, Kohli V, Singh AS, Lekshmi M, Tripathi G, Kumar S, Nayak BB. Virulence genotypes and antimicrobial susceptibility patterns of Arcobacter butzleri isolated from seafood and its environment. Int J Food Microbiol. 2017;18(263):32–7.

    Article  CAS  Google Scholar 

  51. Isidro J, Ferreira S, Pinto M, Domingues F, Oleastro M, Gomes JP, Borges V. Virulence and antibiotic resistance plasticity of Arcobacter butzleri: Insights on the genomic diversity of an emerging human pathogen. Infect Genet Evol. 2020;80:104213.

    Article  CAS  PubMed  Google Scholar 

  52. Badri S, Fassouane A, Bouslikhane M, Filliol I, Hassar M, Cohen N. Relationship between susceptibility to antimicrobials and virulence factors in Escherichia coli isolated from food in Morocco. Int J Food Saf. 2009;11:98–101.

    Google Scholar 

  53. Beceiro A, Tomás M, Bou G. Antimicrobial resistance and virulence: a successful or deleterious association in the bacterial world? Clin Microbiol Rev. 2013;26:185–230.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Zou L-K, Wang H-N, Zeng B, Li J-N, Li X-T, Zhang A-Y, Zhou Y-S, Yang X, Xu C-W, Xia Q-Q. Erythromycin resistance and virulence genes in Enterococcus faecalis from swine in China. New Microbiol. 2011;34:73–80.

    CAS  PubMed  Google Scholar 

  55. van Driessche E, Houf K, Vangroenweghe F, De Zutter L, Hoof JV. Prevalence, enumeration and strain variation of Arcobacter species in the faeces of healthy cattle in Belgium. Vet Microbiol. 2005;105(2):149–54.

    Article  PubMed  Google Scholar 

  56. Ferreira S, Oleastro M, Domingues FC. Occurrence, genetic diversity and antibiotic resistance of Arcobacter sp. in a dairy plant. J Appl Microbiol. 2017;123(4):1019–26.

    Article  CAS  PubMed  Google Scholar 

  57. Wilkes G, Ruecker NJ, Neumann NF, Gannon VP, Jokinen C, Sunohara M, Topp E, Pintar KD, Edge TA, Lapen DR. Spatiotemporal analysis of Cryptosporidium species/genotypes and relationships with other zoonotic pathogens in surface water from mixed-use watersheds. Appl Environ Microbiol. 2013;79(2):434–48.

  58. Collado L, Inza I, Guarro J, Figueras MJ. Presence of Arcobacter spp. in environmental waters correlates with high levels of fecal pollution. Environ Microbiol. 2008;10(6):1635–40.

    Article  PubMed  Google Scholar 

  59. Ramees TS, Suman Kumar M, Dubal ZB, Anjay M, Sivakumar S, Gupta K. Dhama, Javaid M, Milton AAP, Rathore RS, Kumar A. Genotyping of Campylobacter jejuni and C coli from different poultry sources and human by ERIC-PCR. J Vet Pub Health. 2015;13(2):93–8.

    Google Scholar 

  60. Wilkes G, Edge TA, Gannon VP, Jokinen C, Lyautey E, Neumann NF, Ruecker N, Scott A, Sunohara M, Topp E, Lapen DR. Associations among pathogenic bacteria, parasites, and environmental and land use factors in multiple mixed-use watersheds. Water Res. 2011;45(18):5807–25.

    Article  CAS  PubMed  Google Scholar 

  61. Whiteduck-Léveillée K, Whiteduck-Léveillée J, Cloutier M, Tambong JT, Xu R, Topp E, Arts MT, Chao J, Adam Z, André Lévesque C, Lapen DR, Villemur R, Talbot G, Khan IUH. Arcobacter lanthieri sp. nov., isolated from pig and dairy cattle manure. Int J Syst Evol Microbiol. 2015;65(8):2709–16.

    Article  CAS  PubMed  Google Scholar 

  62. Khan IUH, Cloutier M, Libby M, Lapen DR, Wilkes G, Topp E. Enhanced single-tube multiplex PCR assay for detection and identification of six Arcobacter species. J Appl Microbiol. 2017;123(6):1522–32.

    Article  CAS  PubMed  Google Scholar 

  63. Harmon K, Wesley I. Identification of Arcobacter isolates by PCR. Lett Appl Microbiol. 1996;23:241–4.

    Article  CAS  PubMed  Google Scholar 

  64. Clinical Laboratory Standards Institute (CLSI), editor. Performance Standards for Antimicrobial Susceptibility Testing. 21st informational supplement. CLSI M100-S21. Wayne: Clinical and Laboratory Standards Institute; 2011.

  65. Avery BDD, Deckert A, Demczuk W, Dutil D, Finley R, Gow S, Laberge K, Léger D, Mulvey M, Varughese M. Report of the Advisory Committee on Animal Uses of Antimicrobials in Canada: Impact on Resistance and Human Health. Report. Health Canada. 2004.

    Google Scholar 

  66. Gaudreau C, Gilbert H. Comparison of disc diffusion and agar dilution methods for antibiotic susceptibility testing of Campylobacter jejuni subsp. jejuni and Campylobacter coli. J Antimicrob Chemother. 1997;39(6):707–12.

    Article  CAS  PubMed  Google Scholar 

  67. Whiteduck-Léveillée J, Cloutier M, Topp E, Lapen DR, Talbot G, Villemur R, Khan IU. Development and evaluation of multiplex PCR assays for rapid detection of virulence-associated genes in Arcobacter species. J Microbiol Methods. 2016;121:59.

  68. Houf K, De Zutter L, Van Hoof J, Vandamme P. Assessment of the genetic diversity among arcobacters isolated from poultry products by using two PCR-based typing methods. Appl Environ Microbiol. 2002;68(5):2172–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. StatSoft, Inc. Electronic statistics textbook. Tulsa: StatSoft; 2013.

    Google Scholar 

  70. Hunter PR, Gaston MA. Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J Clin Microbiol. 1988;26(11):2465–6.

    Article  CAS  PubMed Central  Google Scholar 

Download references


We thank the field crew and CO-OP students who have assisted with water sample collection and laboratory analysis. We also thank farmers for their cooperation in water sample collection.


This study was funded by Agriculture and Agri-Food Canada (AAFC) under project #s J-001848, J-001012 and J-001367.

Author information

Authors and Affiliations



IUHK, WC and DRL conceived and designed the study. MC, DRL, EC and GW assisted in identifying sites, collected water samples and performed lab analyses. IUHK, WC and DRL drafted and reviewed the manuscript. All co-authors reviewed, edited and approved the final version of the manuscript.

Corresponding author

Correspondence to Izhar U. H. Khan.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors have no competing interest to report.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Khan, I.U.H., Chen, W., Cloutier, M. et al. Pathogenicity assessment of Arcobacter butzleri isolated from Canadian agricultural surface water. BMC Microbiol 24, 17 (2024).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: