Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Molecular characterization of Campylobacter spp. recovered from beef, chicken, lamb and pork products at retail in Australia

  • Rhiannon L. Wallace,

    Roles Data curation, Formal analysis, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation National Centre for Epidemiology and Population Health, The Australian National University, Canberra, Australian Capital Territory, Australia

  • Dieter M. Bulach,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Software, Validation, Visualization, Writing – review & editing

    Affiliations Melbourne Bioinformatics, The University of Melbourne, Melbourne, Victoria, Australia, Microbiological Diagnostic Unit Public Health Laboratory, The Peter Doherty Institute, Melbourne, Victoria, Australia

  • Amy V. Jennison,

    Roles Conceptualization, Data curation, Investigation, Methodology, Resources, Software, Validation, Writing – review & editing

    Affiliation Public Health Microbiology, Forensic and Scientific Services, Queensland Health, Brisbane, Queensland, Australia

  • Mary Valcanis,

    Roles Data curation, Investigation, Methodology, Resources, Validation, Writing – review & editing

    Affiliation Microbiological Diagnostic Unit Public Health Laboratory, The Peter Doherty Institute, Melbourne, Victoria, Australia

  • Angus McLure,

    Roles Data curation, Formal analysis, Software, Writing – review & editing

    Affiliation National Centre for Epidemiology and Population Health, The Australian National University, Canberra, Australian Capital Territory, Australia

  • James J. Smith,

    Roles Data curation, Investigation, Methodology, Resources, Writing – review & editing

    Affiliation Food Safety Standards and Regulation, Health Protection Branch, Queensland Health, Brisbane, Queensland, Australia

  • Trudy Graham,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Public Health Microbiology, Forensic and Scientific Services, Queensland Health, Brisbane, Queensland, Australia

  • Themy Saputra,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation New South Wales Food Authority, NSW Government, Sydney, New South Wales, Australia

  • Simon Firestone,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Melbourne Veterinary School, Faculty of Veterinary and Agricultural Sciences, The University of Melbourne, Parkville, Victoria, Australia

  • Sally Symes,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Department of Health and Human Services, Victoria State Government, Melbourne, Victoria, Australia

  • Natasha Waters,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation ACT Government Analytical Laboratory, Australian Capital Territory Health Directorate, Canberra, Australian Capital Territory, Australia

  • Anastasia Stylianopoulos,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Department of Health and Human Services, Victoria State Government, Melbourne, Victoria, Australia

  • Martyn D. Kirk,

    Roles Conceptualization, Funding acquisition, Methodology, Project administration, Supervision, Writing – review & editing

    Affiliation National Centre for Epidemiology and Population Health, The Australian National University, Canberra, Australian Capital Territory, Australia

  • Kathryn Glass

    Roles Conceptualization, Funding acquisition, Methodology, Project administration, Supervision, Writing – review & editing

    kathryn.glass@anu.edu.au

    Affiliation National Centre for Epidemiology and Population Health, The Australian National University, Canberra, Australian Capital Territory, Australia

Abstract

Australian rates of campylobacteriosis are among the highest in developed countries, yet only limited work has been done to characterize Campylobacter spp. in Australian retail products. We performed whole genome sequencing (WGS) on 331 C. coli and 285 C. jejuni from retail chicken meat, as well as beef, chicken, lamb and pork offal (organs). Campylobacter isolates were highly diverse, with 113 sequence types (STs) including 38 novel STs, identified from 616 isolates. Genomic analysis suggests very low levels (2.3–15.3%) of resistance to aminoglycoside, beta-lactam, fluoroquinolone, macrolide and tetracycline antibiotics. A majority (>90%) of isolates (52/56) possessing the fluoroquinolone resistance-associated T86I mutation in the gyrA gene belonged to ST860, ST2083 or ST7323. The 44 pork offal isolates were highly diverse, representing 33 STs (11 novel STs) and harboured genes associated with resistance to aminoglycosides, lincosamides and macrolides not generally found in isolates from other sources. Prevalence of multidrug resistant genotypes was very low (<5%), but ten-fold higher in C. coli than C. jejuni. This study highlights that Campylobacter spp. from retail products in Australia are highly genotypically diverse and important differences in antimicrobial resistance exist between Campylobacter species and animal sources.

Introduction

Thermophilic Campylobacter coli and Campylobacter jejuni are the most common causes of foodborne bacterial infections worldwide. Campylobacteriosis is also the most frequently notified enteric pathogen under surveillance by OzFoodNet, Australia’s government-established foodborne disease surveillance network [1], with 143.5 cases per 100,000 population reported in 2019 [2]. Although only a small fraction of people who become ill from food they have consumed seek medical attention [3], the global burden of foodborne disease is substantial. In Australia, foodborne gastroenteritis is responsible for an estimated 4.1 million cases annually [4]. Globally, the incidence and prevalence of campylobacteriosis have increased dramatically over the past decade [5]. Campylobacteriosis is often self-limiting and generally does not require medical treatment. However, some cases of Campylobacter infection are associated with serious clinical manifestations such as bacteraemia, reactive arthritis, haemolytic uremic syndrome, meningitis, septicaemia and Guillain-Barré syndrome [5, 6].

The global overuse and misuse of antimicrobial agents in humans, animals and plants, has greatly accelerated the development of resistance to antimicrobials by pathogenic microorganisms. In the United States alone, the CDC estimates that 2 million people become infected with antimicrobial(s)-resistant bacterial pathogens, resulting in at least 23,000 deaths annually [7]. The development of antimicrobial resistance (AMR) and emergence of multidrug resistant pathogens are global concerns for both public health agencies and the agri-food industry. Antimicrobial resistant pathogens increase the risk of an infected individual suffering an adverse health effect, such as reduced treatment efficacy, and increased disease severity, hospitalization and mortality than individuals infected with an antimicrobial-susceptible pathogen [810].

Traditionally, phenotypic-methods have been widely used to characterize C. coli and C. jejuni. However, molecular techniques, offering greater accuracy and specificity, have replaced phenotypic methods. These techniques include pulsed-field gel electrophoresis (PFGE), multi-locus sequence typing (MLST) [11], and more recently next generation sequencing [12]. Whole genome sequencing (WGS) provides the highest possible microbial subtyping resolution available to public health authorities [13], enabling faster detection and identification of resistance determinants/mechanisms in microorganisms [14].

Poultry is Australia’s largest meat commodity, with the average Australian consuming 47.4 kg each year [15]. In 2018–2019, 1.24 million tonnes of chicken meat were produced, representing a gross production value of $2.8 billion [16]. The majority of commercial meat chickens in Australia are grown intensively in conventional barns or sheds with controlled environments which provide favourable conditions for the proliferation of pathogenic microorganisms, such as Campylobacter spp. Australia has strict regulations regarding antimicrobial use in livestock. Fluoroquinolones, colistin and 4th generation cephalosporins have never been registered for use in Australian food-producing animals, gentamycin use is banned and 3rd generation cephalosporin usage remains restricted [17].

Evaluation of isolates from meat at the retail level using WGS is an effective way to identify the risk of human exposure to enteric pathogens, particularly microorganisms harbouring AMR genes. Despite the size of the animal agriculture industry in Australia, little is known about the AMR profiles of foodborne pathogens, such as C. coli and C. jejuni in retail products. The aim of this study was to apply WGS to (i) determine the dominant sequence types (STs), (ii) identify the AMR genotypes and (iii) characterize the genetic relatedness of C. coli and C. jejuni from beef, chicken, lamb and pork products at retail in Australia.

Materials and methods

Sample collection and bacterial isolation

Raw meat and offal products were collected from retail outlets in New South Wales (NSW), Queensland (QLD) and Victoria (VIC), between March 2017 and March 2019, as previously described [18]. Additional samples were collected from retail outlets in the Australian Capital Territory (ACT) over a period of five months (May-Sept. 2018). Briefly, fresh and frozen chicken meat and beef-, chicken-, lamb-, and pork-offal were collected from local supermarkets and butcher shops: Canberra in the ACT; Hunter region and metropolitan suburbs around Sydney in NSW; Brisbane, Toowoomba, Rockhampton, Townsville and Cairns in QLD; and Bendigo and Melbourne in VIC. Samples of chicken meat included whole bird, and breast, drumstick, Maryland (thigh and drumstick), thigh, and wing portions. Samples of beef, chicken, lamb and pork offal included giblets, heart, liver, kidney, neck and tongue. Campylobacter spp. prevalence on beef, lamb and pork meat is generally low. Therefore organs (offal) were sampled in order to obtain a suitable number of isolates for WGS.

C. coli and C. jejuni were isolated and identified in their respective jurisdiction according to ISO 10272–1:2017 [19] and AS 5013.06.2015 [20] with minor modifications [18]. Briefly, meat or offal samples were combined with buffered peptone water or enrichment broth (Bolton or Preston) and agitated manually. Samples were incubated at 37°C for 2–4 h followed by microaerobic incubation (85% N2, 5% O2, and 10% CO2) at 41.5°C for 44 ± 4 h and subsequently identified to genus and species level as described previously [18].

Genomic DNA extraction and whole genome sequencing

Genomic DNA was extracted from C. coli (n = 331) and C. jejuni (n = 285) isolated from beef, chicken, lamb and pork using the QiaSymphony® DSP DNA Mini kit (Qiagen) according to the manufacturer’s protocol. The Nextera XT DNA Library Prep kit (Illumina, San Diego, CA, USA) was used to prepare DNA for sequencing. WGS was performed on the Illumina NextSeq500 with 150 base-pair paired-end reads using the NextSeq 500 Mid Output kit (300 cycles) (Illumina). Table 1 summarizes the number of isolates sequenced from each sample type.

thumbnail
Table 1. Summary of whole genome sequenced Campylobacter isolates, collected from retail outlets in Australia over a period of two years (2017–2019), in this study.

https://doi.org/10.1371/journal.pone.0236889.t001

Bioinformatic analysis

Paired-end sequences were analysed using the Nullarbor pipeline v2 (https://github.com/tseemann/nullarbor). Read quality was evaluated to ensure sufficient depth of coverage (min 50x) and isolate purity using Kraken (lack of contaminating reads) (https://github.com/DerrickWood/kraken). Sequences were trimmed using Trimmomatic v0.38 to remove Illumina Nextera adapters and low-quality sequences. Isolates with a genome size (total assembled bases) that differed by >20% from the average genome size in the analysis group were excluded. Reads were aligned to a reference genome NC_022132 (C. coli) or NC_003912 (C. jejuni) using Snippy v4.3.6/BWA-MEM v0.7.17-r1188. Sequences were de novo assembled into contigs using SPAdes v3.13.0 and annotated using Prokka v1.13.3 (https://github.com/tseemann/prokka), as part of the Nullarbor pipeline.

In silico multilocus sequence typing (MLST) was performed on de novo assembled contigs using a BLAST-based tool (https://github.com/tseemann/mlst) using the PubMLST “Campylobacter jejuni/coli” allele database (https://pubmlst.org/campylobacter/) sited at the University of Oxford [21].

Identification of genetic determinants of resistance

Assembled contigs were screened for known AMR genes using the NCBI’s AMRFinderPlus (https://www.ncbi.nlm.nih.gov/pathogens/antimicrobial-resistance/AMRFinder/) and Abricate v.0.8.12 (https://github.com/tseemann/abricate). Quality of hits were filtered using a cut-off set at 95% nucleic acid sequence identity and 50% sequence coverage. Mutations in two housekeeping genes, gyrA (T→I at amino acid at position 86) and 23S rRNA (nucleotides at position 2074 and 2075), associated with quinolone [22] and macrolide/lincosamide/ketolide [23] resistance, respectively, were investigated. A mutation in the promoter region of the blaOXA-61 gene (G→T at position 57), which inhibits transcription resulting in the isolate becoming sensitive to ampicillin, despite the presence of the gene [24] was also investigated. BLAST [25] was used to extract the respective nucleotide sequences associated with the 23S rRNA and blaOXA-61 genes and multiple sequence alignments were generated using Clustal-Omega [26]. The GyrA amino acid sequence was used as the query for a tblastn to extract input sequences for a multiple alignment of the GyrA protein from each isolate.

Temporal analysis of AMR in C. coli and C. jejuni was assessed by classifying isolates as resistant based on the presence of resistance genes or mutations, which have been shown to have a high concordance with phenotypic resistance [14, 27, 28]. Isolates were grouped by quarter and year based on the date the sample was collected from the retail outlet. If an isolate had one or more of the following (i) resistance genes, aad9, aadE-Cc, ant(6)-la, aph(3’)-llla, blaOXA-184, blaOXA-185, lnu(C), tet(O) or tet(44) or (ii) mutations, blaOXA-61 G57T, gyrA T86I or 23S rRNA A2075G, they were classified as resistant. The prevalence of resistant isolates was determined by dividing the number of isolates with at least one resistance marker in each sampling quarter by the total number of isolates collected in that quarter. Ninety-five percent confidence intervals (CI95) for the prevalence of resistant isolates were calculated using the binom package [29] in R [30] with ggplot2 [31] used to plot associated results.

In order to identify resistance-associated factors, a multivariable logistic regression analysis was performed using the glm function in R, and the questionr package [32] was used to calculate the odds ratios (ORs) and CI95. Isolates were classified as resistant as described previously. Isolates that did not possess any resistance genes or mutations were classified as sensitive. Variables in the model included year of sample collection, region (state or territory), food source and Campylobacter species.

Phylogenetic analysis

Phylogenetic trees were inferred from single nucleotide polymorphisms (SNPs) within the core genome of 331 C. coli and 285 C. jejuni. FastTree v.2.1.10 [33] with the Jukes-Cantor model for building approximation of maximum phylogenetic tree based on core genome (regions of reference genomes to which reads mapped from each of the isolates in the species-level analysis groups). Interactive tree of life (iTOL) v4 was used for visualization [34]. PHYLOViZ Online [35] was used to generate goeBURST minimum spanning trees, based on the seven gene MLST allele profiles, in order to examine the population structure of C. coli and C. jejuni. Clonal complexes were determined by grouping multilocus genotypes that shared four or more identical alleles among the seven loci (aspA, ginA, gltA, glyA, pgm, tkt and uncA) with at least one other genotype in the group.

Nucleotide sequence accession numbers

WGS results of the 616 Campylobacter isolates used in this study were submitted to the National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/). The GenBank accession numbers of individual isolates are listed in S1 Table and are available on bioproject PRJNA591966. Note: sequences from 617 isolates are available in this bioproject, including one Campylobacter hyointestinalis isolate not discussed in this publication.

Results

ST prevalence by source

The 331 C. coli isolates differentiated into 60 STs and of these, 25 were observed for the first time in this study (Fig 1). The STs of four of the C. coli isolates could not be determined (S1 Table). Among isolates where the ST was determined, 94.5% belonged to a single clonal complex (ST828 complex), which was dominated by chicken isolates at the core and pork isolates in the periphery. Among C. coli isolates, the most frequently observed STs were ST827 (n = 56; 16.9%), ST1181 (n = 50; 15.1%), ST825 (n = 34; 10.3%) and ST860 (n = 28; 8.4%). The four most common STs included isolates from at least two animal sources, except ST825 which represented isolates from chicken meat or offal only (n = 34). C. coli isolates from pork were widely distributed among STs, with 28 isolates representing 23 STs. A majority of the pork isolates (n = 18; 64.3%) represented 15 STs not found in any other source.

thumbnail
Fig 1. Distribution of C. coli (n = 327) Sequence Types (STs) in Australian retail products.

Chicken meat (dark blue), chicken offal (light blue), lamb offal (yellow), beef offal (orange) and pork offal (green). Increasing circle size represents a larger number of isolates of the respective ST. Connecting lines infer phylogenetic relatedness and represent STs with four or more loci in common. Red numbers indicate the number of alleles differing between the two adjoining STs.

https://doi.org/10.1371/journal.pone.0236889.g001

The 285 C. jejuni isolates differentiated into 53 STs and among these, 13 novel STs were identified (Fig 2). In contrast to the population structure of C. coli, the population structure of C. jejuni was more fragmented, with the majority of isolates (n = 214; 75.1%) belonging to one of seven clonal complexes. Among C. jejuni isolates, ST50 (n = 51; 17.9%) and ST61 (n = 22; 7.7%) were the most common. Only ST50, ST21 and ST42 contained C. jejuni isolates from all four animal sources. ST61 was dominated by isolates from non-chicken sources, namely lamb (n = 11) and beef (n = 9) offal. Although the C. jejuni isolates were predominantly from chicken meat or offal (n = 197; 69.1%), 10 STs were dominated (>50%) by isolates from non-chicken sources.

thumbnail
Fig 2. Distribution of C. jejuni (n = 285) Sequence Types (STs) in Australian retail products.

Chicken meat (dark blue), chicken offal (light blue), lamb offal (yellow), beef offal (orange) and pork offal (green). Increasing circle size represents a larger number of isolates of that respective ST. Connecting lines infer phylogenetic relatedness and represent STs that have four or more loci in common. Red numbers indicate the number of alleles differing between the two adjoining STs.

https://doi.org/10.1371/journal.pone.0236889.g002

AMR prevalence over time

The prevalence of resistance genes and point mutations known to confer phenotypic resistance to aminoglycoside, β-lactam, lincosamide, quinolone and tetracycline antibiotics, in C. coli and C. jejuni collected over two years is shown in Fig 3. Resistance was not analysed by sampling quarter and animal source due to the small number of isolates obtained from non-chicken sources (S1 Fig). Genes and mutations associated with antibiotic resistance were found in all sampling quarters. The prevalence of resistance genes and mutations was lowest in 2017 (6.3–24.3%) and highest in 2018 (27.7–57.6%).

thumbnail
Fig 3. Temporal analysis of antimicrobial resistance determinants.

Prevalence of Campylobacter (n = 616) possessing at least one genetic marker of antimicrobial resistance from retail beef, chicken, lamb and pork products collected over a period of two years. Genetic determinants of resistance used to classify isolates as resistant are described in section 2.4. Error bars indicate 95% confidence intervals.

https://doi.org/10.1371/journal.pone.0236889.g003

Three variables: year, source and species were associated with either an increased or decreased prevalence of AMR at a significance level of P < 0.05 (Table 2). Compared to 2017, isolates collected in 2018 were associated with in an increased risk of resistance (OR 2.41; CI95 1.58–3.68). Campylobacter spp. from pork offal were associated with an increased risk of resistance compared to chicken meat (OR 3.59; CI95 1.79–7.37). Among the two Campylobacter species tested in the present study, C. coli isolates were associated with an increased risk of resistance (OR 1.79; CI95 1.22–2.64) compared to C. jejuni. In a univariable analysis, QLD isolates were associated with a lower prevalence of resistance (OR 0.42; CI95 0.23–0.77), however the extent of this association was reduced in the multivariable analysis (OR 0.70; CI95 0.35–1.40). It should be noted that isolates were not collected in all sampling quarters in each region (S2 Fig).

thumbnail
Table 2. Multivariable logistic regression analysis of factors potentially associated with resistance in Campylobacter coli and Campylobacter jejuni isolated from retail meat and offal products in Australia, 2017–2019.

https://doi.org/10.1371/journal.pone.0236889.t002

AMR prevalence by species

Nine genes and two mutations associated with antibiotic resistance were identified in the 331 C. coli isolates (Table 3; S3S5 Figs). The blaOXA-61 gene encoding resistance to β-lactams was the most commonly detected resistance gene, present on average in 80.4% of isolates (266/331). However, of the isolates possessing the blaOXA-61 gene, none are predicted to have an active promotor, with all isolates having a G at position 57. Among C. coli isolates, only isolates from pork offal (10/28; 35.7%) contained the lnu(C) gene. The prevalence of several resistance genes, namely aad9, lnu(C) and tet(O), as well as the 23S rRNA A → G mutation at position 2075, were highest in isolates recovered from pork offal compared to isolates from other sources. The T86I mutation in the gyrA gene was present in isolates from all sources, except beef.

thumbnail
Table 3. Summary of genetic determinants of antimicrobial resistance (AMR) present in C. coli recovered from beef, chicken, lamb and pork collected at retail outlets in the Australian Capital Territory, New South Wales, Queensland and Victoria, 2017–2019.

https://doi.org/10.1371/journal.pone.0236889.t003

Ten genes and three mutations associated with antibiotic resistance were identified in the 285 C. jejuni isolates (Table 4; S5 Fig). The blaOXA-61 gene was the most commonly detected antimicrobial resistance gene in C. jejuni, present on average in 79.0% of isolates (225/285). However, only 8.9% of these isolates (20/225) had an active blaOXA-61 promoter (G → T mutation at position 57). The tet(O) gene was present on average in 13.3% (38/285) of C. jejuni isolates and was detected in isolates from all sources except beef offal. The T86I mutation in the protein encoded by the gyrA gene was detected infrequently (27/285; 9.5%) in C. jejuni isolates from beef, chicken and pork products. The A → G mutation in the 23S rRNA gene at position 2075 was only detected in one isolate from pork offal.

thumbnail
Table 4. Summary of genetic determinants of antimicrobial resistance (AMR) in C. jejuni recovered from beef, chicken, lamb and pork collected at retail outlets in the Australian Capital Territory, New South Wales, Queensland and Victoria, 2017–2019.

https://doi.org/10.1371/journal.pone.0236889.t004

Phylogeny and AMR

Phylogenetically, C. coli isolates were highly diverse, forming two genetically-distinct clades, with isolates from pork offal generally more closely related to other pork isolates than those from other sources (S3S5 Figs). Half the pork offal C. coli isolates had a multidrug resistance (MDR) genotype and among these, 57.1% (8 isolates) represented novel STs (S3 Fig). The T86I mutation in the protein encoded by the gyrA gene was restricted primarily to three STs: ST2083 (n = 10), ST7323 (n = 14) and ST860 (n = 24) (S3 & S5 Figs). In contrast, the tet(O) and tet(44) genes were identified in over 35 different STs (S3S5 Figs). Fourteen isolates possessed the 23S rRNA A2075G mutation, representing 13 different STs. Multidrug resistance genotypes were found in 26 (7.6%) of C. coli and 2 (0.7%) C. jejuni (S3 & S5 Figs). C. jejuni isolates were highly diverse and generally did not cluster genetically by food source, except ST61 (S5 Fig). Among STs containing more than one isolate, only ST132 (n = 2), ST829 (n = 7), ST4896 (n = 13) and ST10123 (n = 2) had no resistance genes or mutations detected (S3S5 Figs).

Discussion

Antimicrobial resistance in enteric pathogens is considered a serious global public health problem. Among developed countries, several surveillance programs monitoring AMR of Campylobacter to clinically important antimicrobials have been well established, such as the U.S. National Antimicrobial Resistance Monitoring Systems [36], the Canadian Antimicrobial Resistance Surveillance System [37], the European Food Safety Authority [38] and the European Centre for Disease Prevention and Control [39]. In Australia, there is no national surveillance program for AMR monitoring in Campylobacter spp., limiting our knowledge of resistance in Campylobacter. The present study represents the first longitudinal genomic study of Campylobacter in meat at the retail level in Australia, providing insights into the prevalence of AMR genetic markers, as well as the genetic relatedness of isolates from different animal sources. Here we chose to collect samples at retail as this is the closest point before the product comes in contact with consumers, likely providing good representation of Campylobacter spp. that may cause illness in humans.

Although, we did not perform phenotypic antimicrobial susceptibility testing on the isolates in this study, our previous work [40] and reports by others [14, 28, 41], has shown a high concordance between resistance genotype and resistance phenotype in Campylobacter. Genetic evidence for resistance to tetracycline and fluoroquinolone antimicrobials was found in over 11 and 7.5% of C. coli isolates from chicken meat, respectively. These findings are comparable with those in a recent report published by the Australian Chicken Meat Federation (ACMF) that showed 5.2 and 3.1% of C. coli from chicken cecal samples were phenotypically resistant to fluoroquinolone and tetracycline, respectively [42]. Among C. jejuni isolates from retail chicken meat, genetic evidence for resistance to tetracycline and fluoroquinolone was 18.6 and 11.5%, respectively. These findings are slightly lower than those reported by the ACMF where 22.2 and 14.8% of C. jejuni isolates from chicken cecal samples were phenotypically resistant to tetracycline and ciprofloxacin, respectively [42].

Many of our results are understandable in the context of antibiotic usage patterns in Australian food-producing animals and recent literature on the evolution of resistance patterns. Possible explanations for the finding of mutations suggesting quinolone resistance in Australian food-producing animals, in the absence of quinolone usage, are a fitness advantage [43] in the absence of antibiotic selection, or co-selection [44] where found associated with other resistance patterns. The Australian poultry meat industry has suggested the initial introduction of fluoroquinolone resistance could have been from reverse zoonosis [42]. The moderate levels of aminoglycoside resistance could be due to usage of neomycin or spectinomycin rather than gentamycin, which is banned from use in food-producing animals. The range of mutations associated with AMR present in isolates from pork is reflective of the therapeutic usage of first-line antibiotics. Amoxycillin is commonly used to treat endemic respiratory and enteric diseases of pigs, including Actinobacillus pleuroneumonia, colibacillosis, enterotoxic E. coli, diarrhoea and Glasser’s disease (Haemophilus parasuis). Tetracyclines are also commonly used to treat colibacillosis in pigs. Macrolides are commonly used in pigs to treat pneumonia due to mycoplasmosis, swine dysentery (Brachyspira hyodysenteriae) and proliferative enteropathy (Lawsonia intracellularis). It is worth noting that there are no antibiotics used for growth promotion in Australian pigs.

Globally, resistance to clinically important antimicrobials is a serious threat to public health. A recent report by the European Food Safety Authority (EFSA) showed levels of AMR, particularly to tetracycline and ciprofloxacin, are very high in human and animal Campylobacter isolates in some European countries [45]. Prevalence of fluoroquinolone resistance in Campylobacter in Europe from food animals is highly variable, ranging from 1.2% in Norway [46], 44.0% in Belgium [47] and 82.9% in Italy [48] to 90.0% in Spain [49]. In Canada, ciprofloxacin resistance in Campylobacter isolated from retail chicken meat ranged from 6% in Québec to 35% in British Columbia [50]. In China, a recent report found the tet(O) gene in 98% of Campylobacter and the GyrA mutation at codon 86 and the 23S rRNA A2075G point mutation in 99 and 37% of isolates, respectively [51]. In China, antibiotic use is more than five times higher than the international average, primarily due to widespread misuse associated with growth promotion in feed and veterinary use on farms [52]. In Australia, the very low prevalence of AMR genetic markers in Campylobacter spp. from retail products is a testament to good antibiotic stewardship by the Australian agriculture industry. However, our findings suggest the prevalence of AMR in Campylobacter isolates from retail products may be increasing, highlighting the need for routine monitoring of resistance in Campylobacter in the agriculture industry. Genetic mutations conferring AMR may arise after continued exposure to antibiotics, but persist due to neutral selection, e.g. mutation in the gyrA gene conferring fluoroquinolones resistance. Although mutations in the 23S gene conferring erythromycin resistance are detrimental and are rapidly lost in the absence of selection. Future studies should examine resistance in Campylobacter spp. over an extended period of time to determine if the prevalence of resistance is increasing or if our findings can be regarded as background levels of resistance in Australian Campylobacter.

Campylobacteriosis is generally self-limiting and does not require medical treatment. However, people at increased risk of severe complications, such as the immunocompromised or elderly, will likely be prescribed antibiotics to prevent bacteraemia or sequelae [53]. This highlights the importance of monitoring prevalence of antimicrobial resistant Campylobacter in food. We identified five genes; aad9, aadE-Cc, ant(6)-Ia, ant(6)-Ib and aph(3’)-IIIa, associated with resistance to aminoglycosides, as well as lnu(c), associated with lincomycin resistance, in C. coli but not C. jejuni. Additionally, the 23S rRNA mutation associated with resistance to macrolides was rare in both Campylobacter species, but more prevalent in C. coli than C. jejuni. Similarly, prevalence of MDR genotypes was ten-fold higher in C. coli than C. jejuni. In cases of campylobacteriosis in at risk individuals, macrolides or fluoroquinolones may be prescribed, with macrolides preferred due to their low rate of resistance. Our findings highlight that while a particular antimicrobial may be effective against C. jejuni, it may not be as effective for C. coli. In Australia, pathology laboratories do not routinely identify Campylobacter to species level. However, in cases where antibiotics are required, species identification may improve treatment success.

Australia is a large, geographically diverse, island nation supporting a distinct and diverse Campylobacter population evident in the findings from this study. Among the 616 Campylobacter isolates, representing 113 STs, 38 STs had not previously been reported anywhere else in the world. Many of the more common STs found in Australia, such as ST825, ST827 and ST50, have been isolated from various sources globally. C. jejuni isolates were dominated by ST50 which has been found in poultry in Israel [54], Poland [55], Denmark [56] and South Korea [57], among other countries. To date, ST50 has over 3,500 submissions on the PubMLST database. Unlike ST50 which was found in all animal sources, ST61 was dominated by isolates from non-poultry sources, namely beef and lamb offal. This finding agrees with other reports indicating ST61 is commonly recovered from ovine, bovine and human samples, but rarely from poultry [58, 59]. There was no general association between different genetic markers of resistance and ST, with a few exceptions. Most isolates harbouring the GyrA T86I mutation, which is associated with resistance to ciprofloxacin, belonged to ST860, ST2083 or ST7323. Future studies should examine the prevalence of these STs in clinical isolates. The potential for the novel STs identified for the first time in the present study to cause illness in humans remains to be determined.

Of note, the blaOXA-61 gene, shown to confer resistance to β-lactams [60], was present in 79.7% of Campylobacter isolates. However, rather than the presence of blaOXA-61 alone, a single nucleotide mutation (G → T transversion) upstream of blaOXA-61 has been shown to up-regulate expression of the gene, resulting in a high level of β-lactam resistance in C. jejuni [24]. In a previous study, we found 100% of Campylobacter isolates possessing the G → T mutation were resistant to ampicillin, while C. jejuni isolates possessing the blaOXA-61 gene without this mutation, remained susceptible to ampicillin [40]. This information, taken with the other beta-lactamase genes we identified, suggests less than 15% of the C. jejuni isolates in the present study are likely to be resistant to β-lactam antibiotics. Beta-lactams are not recommended for treating campylobacteriosis as Campylobacter is intrinsically resistant to this class of antibiotic [48]. However, our findings suggest this intrinsic resistance may be less prevalent in Campylobacter from retail products in Australia, particularly in C. jejuni.

MDR is defined as resistance to three or more classes of antimicrobials. Among the 616 Campylobacter isolates examined in this study, 26 C. coli and two C. jejuni possessed MDR genotypes. However, among C. coli possessing a MDR genotype, more than 75% of these isolates had the blaOXA-61 gene. Phenotypic testing is necessary to determine if the presence of the blaOXA-61 gene alone in C. coli is in fact associated with phenotypic β-lactam antibiotic resistance. The majority of MDR C. coli represented STs observed for the first time in this study, and the role of these novel STs in human illness is unknown. Overall, our results show that MDR is very rare (< 5%) in Australia. This is comparable to findings in Poland where no MDR Campylobacter were found in poultry [61]. In contrast, resistance to multiple drugs (≥ 4) was found in 31.6% of C. jejuni isolated from poultry meat and related samples in Northern India [62]. Similarly, 8.6% of C. jejuni and 67.6% of C. coli isolates from diarrheal patients and poultry meat in Shanghai, China were identified as MDR [51]. Previous reports have indicated MDR Campylobacter can spread from the food supply chain into the human population [51]. In the present study, we noted differences in the prevalence of genetic determinants of resistance by food source. Generally, genes associated with resistance to aminoglycosides (except aadE-Cc), lincosamides and macrolides were only found in isolates from pork offal. Similarly, the majority (61.5%) of isolates harbouring MDR genotypes were recovered from pork offal. Although the prevalence of AMR is low in Australian Campylobacter spp., Australia should consider establishing integrated surveillance systems to monitor resistance and the transmission of Campylobacter from food to humans.

Our study has some limitations. Although we collected more than 700 samples from non-poultry sources, the majority of isolates that were sequenced (78.6%) were from chicken meat or chicken offal, due to the relatively low prevalence in beef, lamb and pork offal (14–38%) [18]. We may have achieved a higher recovery rate if we had sampled prior to slaughter rather than at retail as chilling significantly reduces the levels of Campylobacter [63]. Although Australians consume more muscle meat than offal, we also chose to collect offal samples of non-poultry products as the prevalence of Campylobacter has been reported to be higher in offal than muscle meat [64]. We found for the more common STs (> 10 isolates/ST) the Campylobacter isolates recovered from chicken meat or chicken offal products belonged to similar STs. As a result, the isolates recovered from offal likely provide a good representation of the STs that may be present on beef, lamb and pork meat.

Conclusions

Our results indicate Australia’s AMR prevalence in Campylobacter spp. from retail products is very low. Our results also suggest prevalence of resistance in Campylobacter spp. from foods of animal origin may be increasing, but ongoing surveillance is needed to confirm such a trend. Although consumption of contaminated poultry is well established as a key risk factor for campylobacteriosis, foods derived from other animals can result in Campylobacter infection. We found that isolates from pork represented a diverse array of STs, many not found among isolates from beef, chicken or lamb. MDR prevalence was also higher in isolates from pork offal, with many of these isolates representing novel STs. As these novel STs were first reported in this study, their significance in human health remains to be determined. Routine surveillance aimed at identifying species and characterizing types and resistance determinants from WGS data of Campylobacter spp. from food and humans will enable early detection of emergent AMR clones and ultimately assist in maintaining the low prevalence of AMR in Australia.

Supporting information

S1 Fig. Summary of Campylobacter isolates, from retail beef, chicken, lamb or pork products, by sampling quarter.

https://doi.org/10.1371/journal.pone.0236889.s001

(TIFF)

S2 Fig. Summary of Campylobacter isolates, collected from retail products in the Australian Capital Territory (ACT), New South Wales (NSW), Queensland (QLD) and Victoria (VIC), by sampling quarter.

https://doi.org/10.1371/journal.pone.0236889.s002

(TIFF)

S3 Fig. Phylogenetic analysis of Campylobacter coli group 1 (n = 317) isolated from retail beef, chicken, lamb or pork products.

The circle lanes from inner to outer order indicate: food source, number of antimicrobial resistance (AMR) genes or mutations, type of AMR gene or mutation (coloured by antibiotic class) and MLST number. Novel STs identified in this study are shown in red text.

https://doi.org/10.1371/journal.pone.0236889.s003

(JPG)

S4 Fig. Phylogenetic analysis of Campylobacter coli group 2 (n = 16) isolated from retail beef, chicken, lamb or pork products.

The lanes from inner to outer order indicate: food source, number of antimicrobial resistance (AMR) genes, type of AMR gene (coloured by antibiotic class) and MLST number. Novel STs identified in this study are shown in red text.

https://doi.org/10.1371/journal.pone.0236889.s004

(JPG)

S5 Fig. Phylogenetic analysis of Campylobacter jejuni (n = 286) isolated from retail beef, chicken, lamb or pork products.

The circle lanes from inner to outer order indicate: food source, number of antimicrobial resistance (AMR) genes or mutations, type of AMR gene or mutation (coloured by antibiotic class) and MLST number. Novel STs identified in this study are shown in red text.

https://doi.org/10.1371/journal.pone.0236889.s005

(JPG)

S1 Table. Accession numbers of Campylobacter coli and Campylobacter jejuni collected from retail beef, chicken, lamb and pork products from Australian retail outlets.

https://doi.org/10.1371/journal.pone.0236889.s006

(PDF)

Acknowledgments

The CampySource Project Team comprises three working groups and a reference panel. The working groups focus on: food and animal sampling, epidemiology and modelling, and genomics. The reference panel includes expert representatives from government and industry. The study includes the following partner organisations: the Australian National University, Massey University, University of Melbourne, Queensland Health, Queensland Health Forensic and Scientific Services, New South Wales Food Authority, New South Wales Health, Hunter New England Health, Victorian Department of Health and Human Services, Food Standards Australia New Zealand, Commonwealth Department of Health and AgriFutures Australia–Chicken Meat Program. CampySource also collaborates with the following organisations: ACT Health, Sullivan Nicolaides Pathology, University of Queensland, Primary Industries and Regions South Australia, Department of Health and Human Services Tasmania, Meat and Livestock Australia, and New Zealand Ministry for Primary Industries. The CampySource Project Team consists of: Nigel P French, Massey University, New Zealand; Mary Valcanis, The University of Melbourne; Dieter Bulach, The University of Melbourne; Emily Fearnley, The Australian National University and South Australian Department for Health and Wellbeing; Russell Stafford, Queensland Health; Amy Jennison, Queensland Health; Trudy Graham, Queensland Health; Keira Glasgow, Health Protection NSW; Kirsty Hope, Health Protection NSW; Themy Saputra, NSW Food Authority; Craig Shadbolt, NSW Food Authority; Arie H Havelaar, The University of Florida, USA; Joy Gregory, Department of Health and Human Services, Victoria; James Flint, Hunter New England Health; Simon Firestone, The University of Melbourne; James Conlan, Food Standards Australia New Zealand; Ben Daughtry, Food Standards Australia New Zealand; James J Smith, Queensland Health; Heather Haines, Department of Health and Human Services, Victoria; Sally Symes, Department of Health and Human Services, Victoria; Barbara Butow, Food Standards Australia New Zealand; Liana Varrone, The University of Queensland; Linda Selvey, The University of Queensland; Tim Sloan-Gardner, ACT Health; Deborah Denehy, ACT Health; Radomir Krsteski, ACT Health; Natasha Waters, ACT Health; Kim Lilly, Hunter New England Health; Julie Collins, Hunter New England Health; Tony Merritt, Hunter New England Health; Rod Givney, Hunter New England Health; Joanne Barfield, Hunter New England Health; Ben Howden, The University of Melbourne; Kylie Hewson, AgriFutures Australia–Chicken Meat Program; Dani Cribb, The Australian National University; Rhiannon Wallace, The Australian National University; Angus McLure, The Australian National University; Ben Polkinghorne, The Australian National University; Cameron Moffatt, The Australian National University; Martyn Kirk, The Australian National University; and Kathryn Glass, The Australian National University. The authors thank Professor Glenn Browning (University of Melbourne) for critically reviewing the manuscript and providing comments on the use of antimicrobials in food-producing animals in Australia.

References

  1. 1. OzFoodNetWorkingGroup. Monitoring the incidence and causes of diseases potentially transmitted by food in Australia Annual report of the OzFoodNet network 2011. Commun Dis Intell Q Rep. 2015;39(2):E236–E64. pmid:26234259
  2. 2. Australian Government Department of Health. National Notifiable Disease Surveillance System- Disease notification rates, Australia, 1991 to 2019 and year-to-date notifications for 2020 2020 [Available from: http://www9.health.gov.au/cda/source/rpt_2.cfm].
  3. 3. Havelaar AH, Kirk MD, Torgerson PR, Gibb HJ, Hald T, Lake RJ, et al. World Health Organization Global Estimates and Regional Comparisons of the Burden of Foodborne Disease in 2010. PLoS medicine. 2015;12(12):e1001923. pmid:26633896
  4. 4. Kirk M, Ford L, Glass K, Hall G. Foodborne Illness, Australia, Circa 2000 and Circa 2010. Emerging infectious diseases. 2014;20(11):1857–64. pmid:25340705
  5. 5. Kaakoush NO, Castano-Rodriguez N, Mitchell HM, Man SIM. Global Epidemiology of Campylobacter Infection. Clin Microbiol Rev. 2015;28(3):687–720. pmid:26062576
  6. 6. McCarthy N, Giesecke J. Incidence of Guillain-Barre syndrome following infection with Campylobacter jejuni. Am J Epidemiol. 2001;153(6):610–4. pmid:11257070
  7. 7. Centers for Disease Control. Antibiotic/ antimicrobial resistance 2016 [Available from: https://www.cdc.gov/drugresistance/about.html].
  8. 8. Helms M, Simonsen J, Olsen KEP, Molbak K. Adverse health events associated with antimicrobial drug resistance in Campylobacter species: A registry-based cohort study. J Infect Dis. 2005;191(7):1050–5. pmid:15747238
  9. 9. Travers K, Barza M. Morbidity of infections caused by antimicrobial-resistant bacteria. Clinical Infectious Diseases. 2002;34:S131–S4. pmid:11988884
  10. 10. Molbak K. Human health consequences of antimicrobial drug-resistant Salmonella and other foodborne pathogens. Clinical Infectious Diseases. 2005;41(11):1613–20. pmid:16267734
  11. 11. Pfaller MA. Molecular epidemiology in the care of patients. Arch Pathol Lab Med. 1999;123(11):1007–10. pmid:10539897
  12. 12. Cantero G, Correa-Fiz F, Ronco T, Strube M, Cerda-Cuellar M, Pedersen K. Characterization of Campylobacter jejuni and Campylobacter coli Broiler Isolates by Whole-Genome Sequencing. Foodborne Pathog Dis. 2018;15(3):145–52. pmid:29256637
  13. 13. World Health Organization. Whole genome sequencing for foodborne disease surveillance: landscape paper. World Health Organization. World Health Organization,; 2018. Report No.: 978-92-4-151386-9.
  14. 14. Whitehouse CA, Young S, Li C, Hsu CH, Martin G, Zhao SH. Use of whole-genome sequencing for Campylobacter surveillance from NARMS retail poultry in the United States in 2015. Food microbiology. 2018;73:122–8. pmid:29526197
  15. 15. Australian Bureau of Agricultural and Resource Economics and Sciences. Pigs and chickens outlook—March 2020. 2020.
  16. 16. Australian Chicken Meat Federation. Australian Industry Facts & Figures 2018 [Available from: https://www.chicken.org.au/facts-and-figures/].
  17. 17. Pesticides Australian and Veterinary Medicines Authority. Antibiotic resistance in animals. A report for the APVMA. 2017.
  18. 18. Walker LJ, Wallace RL, Smith JJ, Graham T, Saputra T, Symes S, et al. Prevalence of Campylobacter coli and Campylobacter jejuni in Retail Chicken, Beef, Lamb, and Pork Products in Three Australian States. J Food Prot. 2019;82(12):2126–34. pmid:31729918
  19. 19. International Organization for Standardization. ISO 10272–1:2017 Microbiology of the food chain—Horizontal method for detection and enumeration of Campylobacter spp.—Part 1: Detection method. International Organization for Standardization. Geneva, Switzerland2017.
  20. 20. Australia Standards. Australian standards for food microbiology—Examination for specific organisms—Campylobacter. AS 5013.6:20152015.
  21. 21. Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018;3:124. pmid:30345391
  22. 22. Hakanen A, Jalava J, Kotilainen P, Jousimies-Somer H, Siitonen A, Huovinen P. gyrA polymorphism in Campylobacter jejuni: detection of gyrA mutations in 162 C. jejuni isolates by single-strand conformation polymorphism and DNA sequencing. Antimicrobial agents and chemotherapy. 2002;46(8):2644–7. pmid:12121947
  23. 23. Ladely SR, Meinersmann RJ, Englen MD, Fedorka-Cray PJ, Harrison MA. 23S rRNA gene mutations contributing to macrolide resistance in Campylobacter jejuni and Campylobacter coli. Foodborne Pathog Dis. 2009;6(1):91–8. pmid:19014274
  24. 24. Zeng X, Brown S, Gillespie B, Lin J. A single nucleotide in the promoter region modulates the expression of the beta-lactamase OXA-61 in Campylobacter jejuni. J Antimicrob Chemother. 2014;69(5):1215–23. pmid:24408987
  25. 25. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402. pmid:9254694
  26. 26. Sievers F, Higgins DG. Clustal Omega, accurate alignment of very large numbers of sequences. Methods in molecular biology (Clifton, NJ). 2014;1079:105–16.
  27. 27. Fabre A, Oleastro M, Nunes A, Santos A, Sifre E, Ducournau A, et al. Whole-Genome Sequence Analysis of Multidrug-Resistant Campylobacter Isolates: a Focus on Aminoglycoside Resistance Determinants. Journal of clinical microbiology. 2018;56(9).
  28. 28. Zhao S, Tyson GH, Chen Y, Li C, Mukherjee S, Young S, et al. Whole-Genome Sequencing Analysis Accurately Predicts Antimicrobial Resistance Phenotypes in Campylobacter spp. Appl Environ Microbiol. 2016;82(2):459–66. pmid:26519386
  29. 29. Dorai-Raj S. binom: Binomial Confidence Intervals For Several Parameterizations. R package version 1.1–1 2014 [Available from: https://CRAN.R-project.org/package=binom].
  30. 30. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria 2019 [Available from: https://www.R-project.org/].
  31. 31. Wickham H. ggplot2: Elegant Graphics for Data Analysis New York, USA: Springer-Verlag; 2016 [Available from: https://ggplot2.tidyverse.org].
  32. 32. Barnier J, Briatte F., Larmarange J. questionr: Functions to Make Surveys Processing Easier. R package version 0.7.0. 2018 [Available from: https://CRAN.R-project.org/package=questionr].
  33. 33. Price MN, Dehal PS, Arkin AP. FastTree 2—approximately maximum-likelihood trees for large alignments. PloS one. 2010;5(3):e9490. pmid:20224823
  34. 34. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. 2019;47(W1):W256–w9. pmid:30931475
  35. 35. Francisco AP, Bugalho M, Ramirez M, Carrico JA. Global optimal eBURST analysis of multilocus typing data using a graphic matroid approach. BMC Bioinformatics. 2009;10:152. pmid:19450271
  36. 36. Centers for Disease Control and Prevention. National Antimicrobial Resistance Monitoring System for Enteric Bacteria (NARMS) 2019 [Available from: https://www.cdc.gov/narms/index.html].
  37. 37. Government of Canada. Canadian Antimicrobial Resistance Surveillance System—Update 2018: Executive Summary 2018 [Available from: https://www.canada.ca/en/public-health/services/publications/drugs-health-products/canadian-antimicrobial-resistance-surveillance-system-2018-report-executive-summary.html].
  38. 38. European Food Safety Authority. Antimicrobial Resistance 2020 [Available from: http://www.efsa.europa.eu/en/topics/topic/antimicrobial-resistance].
  39. 39. European Centre for Disease Prevention and Control. European Antimicrobial Resistance Surveillance Network (EARS-Net) 2020 [Available from: https://www.ecdc.europa.eu/en/about-us/partnerships-and-networks/disease-and-laboratory-networks/ears-net].
  40. 40. Wallace RL, Bulach D, McLure C, Varrone L, Jennison AV, Valcanis M, et al. Antimicrobial Resistance of Campylobacter spp. Causing Human Infection in Australia: An International Comparison. Microb Drug Resist. 2020;(Under revision).
  41. 41. Kayman T, Abay S, Aydin F, Sahin O. Antibiotic resistance of Campylobacter jejuni isolates recovered from humans with diarrhoea in Turkey. J Med Microbiol. 2019;68(2):136–42. pmid:30540246
  42. 42. Australian Chicken Meat Federation. Surveillance for antimicrobial resistance in enteric commensals and pathogens in Australian meat chickens. 2018 October 2018,.
  43. 43. Luo N, Pereira S, Sahin O, Lin J, Huang S, Michel L, et al. Enhanced in vivo fitness of fluoroquinolone-resistant Campylobacter jejuni in the absence of antibiotic selection pressure. Proc Natl Acad Sci U S A. 2005;102(3):541–6. pmid:15634738
  44. 44. Ingram PR, Rogers BA, Sidjabat HE, Gibson JS, Inglis TJJ. Co-selection may explain high rates of ciprofloxacin non-susceptible Escherichia coli from retail poultry reared without prior fluoroquinolone exposure. J Med Microbiol. 2013;62(Pt 11):1743–6. pmid:24136884
  45. 45. EFSA ECDC. The European Union summary report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2017. EFSA Journal. 2019;17(2).
  46. 46. Norstrom M, Hofshagen M, Stavnes T, Schau J, Lassen J, Kruse H. Antimicrobial resistance in Campylobacter jejuni from humans and broilers in Norway. Epidemiol Infect. 2006;134(1):127–30. pmid:16409659
  47. 47. Van Looveren M, Daube G, De Zutter L, Dumont JM, Lammens C, Wijdooghe M, et al. Antimicrobial susceptibilities of Campylobacter strains isolated from food animals in Belgium. J Antimicrob Chemother. 2001;48(2):235–40. pmid:11481294
  48. 48. Marotta F, Garofolo G, di Marcantonio L, Di Serafino G, Neri D, Romantini R, et al. Antimicrobial resistance genotypes and phenotypes of Campylobacter jejuni isolated in Italy from humans, birds from wild and urban habitats, and poultry. PloS one. 2019;14(10):e0223804. pmid:31603950
  49. 49. Saenz Y, Zarazaga M, Lantero M, Gastanares MJ, Baquero F, Torres C. Antibiotic resistance in Campylobacter strains isolated from animals, foods, and humans in Spain in 1997–1998. Antimicrobial agents and chemotherapy. 2000;44(2):267–71. pmid:10639348
  50. 50. Government of Canada. Canadian Integrated Program for Antimicrobial Resistance Surveillance (CIPARS) 2016 Annual Report. Guelph, Ontario, Canada: Public Health Agency of Canada; 2018.
  51. 51. Du YJ, Wang CQ, Ye YL, Liu Y, Wang AM, Li Y, et al. Molecular Identification of Multidrug-Resistant Campylobacter Species From Diarrheal Patients and Poultry Meat in Shanghai, China. Frontiers in microbiology. 2018;9:8.
  52. 52. Wu Z. Antibiotic use and antibiotic resistance in food-producing animals in China. 2019. Contract No.: TAD/CA/APM/WP(2018)19/FINAL.
  53. 53. Yang Y, Feye KM, Shi Z, Pavlidis HO, Kogut M, A JA, et al. A Historical Review on Antibiotic Resistance of Foodborne Campylobacter. Frontiers in microbiology. 2019;10:1509. pmid:31402900
  54. 54. Rokney A, Valinsky L, Moran-Gilad J, Vranckx K, Agmon V, Weinberger M. Genomic Epidemiology of Campylobacter jejuni Transmission in Israel. Frontiers in microbiology. 2018;9(2432).
  55. 55. Wieczorek K, Denis E, Lachtara B, Osek J. Distribution of Campylobacter jejuni multilocus sequence types isolated from chickens in Poland. Poult Sci. 2017;96(3):703–9. pmid:27702925
  56. 56. Joensen KG, Kiil K, Gantzhorn MR, Nauerby B, Engberg J, Holt HM, et al. Whole-Genome Sequencing to Detect Numerous Campylobacter jejuni Outbreaks and Match Patient Isolates to Sources, Denmark, 2015–2017. Emerging infectious diseases. 2020;26(3):523–32. pmid:32091364
  57. 57. Oh JY, Kwon YK, Wei B, Jang HK, Lim SK, Kim CH, et al. Epidemiological relationships of Campylobacter jejuni strains isolated from humans and chickens in South Korea. Journal of microbiology (Seoul, Korea). 2017;55(1):13–20.
  58. 58. Dingle KE, Colles FM, Wareing DR, Ure R, Fox AJ, Bolton FE, et al. Multilocus sequence typing system for Campylobacter jejuni. Journal of clinical microbiology. 2001;39(1):14–23. pmid:11136741
  59. 59. Mullner P, Spencer SE, Wilson DJ, Jones G, Noble AD, Midwinter AC, et al. Assigning the source of human campylobacteriosis in New Zealand: a comparative genetic and epidemiological approach. Infect Genet Evol. 2009;9(6):1311–9. pmid:19778636
  60. 60. Alfredson DA, Korolik V. Isolation and expression of a novel molecular class D beta-lactamase, OXA-61, from Campylobacter jejuni. Antimicrobial agents and chemotherapy. 2005;49(6):2515–8. pmid:15917560
  61. 61. Wozniak-Biel A, Bugla-Ploskonska G, Kielsznia A, Korzekwa K, Tobiasz A, Korzeniowska-Kowal A, et al. High Prevalence of Resistance to Fluoroquinolones and Tetracycline Campylobacter Spp. Isolated from Poultry in Poland. Microb Drug Resist. 2018;24(3):314–22. pmid:28628752
  62. 62. Khan JA, Rathore RS, Abulreesh HH, Qais FA, Ahmad I. Prevalence and Antibiotic Resistance Profiles of Campylobacter jejuni Isolated from Poultry Meat and Related Samples at Retail Shops in Northern India. Foodborne Pathog Dis. 2018;15(4):218–25. pmid:29377719
  63. 63. Di D, Marotta , Nuvoloni , Zilli , Neri , Di S, et al. Prevalence, Population Diversity and Antimicrobial Resistance of Campylobacter coli Isolated in Italian Swine at Slaughterhouse. Microorganisms. 2020;8(2).
  64. 64. Noormohamed A, Fakhr MK. A higher prevalence rate of Campylobacter in retail beef livers compared to other beef and pork meat cuts. Int J Environ Res Public Health. 2013;10(5):2058–68. pmid:23698698