WGS-Based Prediction and Analysis of Antimicrobial Resistance in Campylobacter jejuni Isolates From Israel

Introduction

Whole genome sequencing (WGS) has revolutionized foodborne pathogen surveillance practices (Jagadeesan et al., 2019). This has become achievable due to the rapid developments in high-throughput sequencing technologies, allowing affordable, real-time, large-scale WGS of foodborne pathogens. This technology can produce a large amount of fast, highly accurate, and reliable information permitting species identification and typing, phylogenetic analyses, outbreak investigation, and determination of virulence and resistance traits as a one-stop-shop (Motro and Moran-Gilad, 2017; Ribot et al., 2019). WGS is increasingly introduced as the primary method for foodborne pathogen surveillance in public health laboratories, replacing the preexisting conventional methods for species identification and typing such as serotyping, pulsed-field gel electrophoresis (PFGE), multilocus variable-number tandem-repeat analysis (MLVA), or multilocus sequence typing (MLST) (Inns et al., 2017; Nadon et al., 2017; Revez et al., 2017; Ribot et al., 2019). The superior discriminatory power of WGS compared to the above typing methods has been established in the detection, investigation, and source tracking of local, multi-state, and multi-national foodborne outbreaks (Crowe et al., 2017; Inns et al., 2017; Llarena et al., 2017; Rumore et al., 2018; Gerner-Smidt et al., 2019; Kubota et al., 2019) and is becoming the new gold standard for molecular surveillance of foodborne pathogens (Moran-Gilad, 2017; Gerner-Smidt et al., 2019).

Campylobacter jejuni is a major foodborne pathogen and a leading cause of human gastroenteritis globally, thus it is the target for surveillance in many countries. WGS is rapidly becoming the first line method for Campylobacter subtyping in many public health laboratories, mostly based on whole genome (wg) multilocus sequence typing (MLST) (wgMLST) and core genome MLST (cgMLST) (Nadon et al., 2017; Ribot et al., 2019).

The incidence of C. jejuni infections in Israel is among the highest in industrialized countries, with an alarming burden on young children (Weinberger et al., 2013). wgMLST-based epidemiologic analyses of Israeli C. jejuni isolates suggest that poultry and cattle are likely food sources of infection in Israel, with different strains persisting over years throughout the food chain (Rokney et al., 2018). The analyses also portrayed the spectrum of virulence determinants harbored by the Israeli isolates.

Campylobacter species are becoming increasingly resistant to the antibiotics relevant to the treatment of clinical infection. Resistance to fluoroquinolones is widespread in both C. jejuni and C. coli. In the European Union (EU) countries the overall resistance rate to ciprofloxacin among C. jejuni human isolates was 58% in 2017, with the highest rates reported from Portugal (97%), Lithuania (92%), Spain (87%), Estonia (84%), and Cyprus (80.0%). Resistance among C. coli human isolates was even higher, reaching an overall rate of 63% and highest rates in Estonia (100%), Portugal (100%), Lithuania (97%), and Spain (96%) (European Food Safety Authority and European Centre for Disease Prevention and Control, 2019). Close to 90% resistance rates were also noted in Peru and China (CDC, 2018; Sproston et al., 2018). Resistance to macrolides is also increasing globally, although not at the same pace. It remains relatively low for C. jejuni in Europe (<3% overall) and the United States (<5%), but higher for C. coli (<15% in the United States and ~13% in EU countries) (CDC, 2018, 2019; Florez-Cuadrado et al., 2018; European Food Safety Authority and European Centre for Disease Prevention and Control, 2019). Resistance among isolates from farm and wild animals may be even higher (CDC, 2018, 2019; Florez-Cuadrado et al., 2018; European Food Safety Authority and European Centre for Disease Prevention and Control, 2019). Campylobacter is listed in the most recent report from the Centers for Disease Control and Prevention (CDC) Antibiotic Resistance Threats in the United States (CDC, 2019) and named by the World Health Organization (WHO) as one of the 12 bacteria that pose the greatest threat to human health due to their antibiotic resistance (Sproston et al., 2018).

In Israel, the reported data on Campylobacter resistance is scarce. In a study from Northern Israel, a sample of 66 stool isolates from the years 2015–2016 was tested by both the Sensititre® broth microdilution plate and the E-test methods. Resistance rates to fluoroquinolones were alarmingly high for both C. jejuni and C. coli, reaching >95%, but were low for macrolides (<3%) (Azrad et al., 2018). Similar results were reported in the most recent periodic publication of the Government Central Laboratories, Israeli Ministry of Health, Jerusalem (available at https://www.health.gov.il/PublicationsFiles/LAB_JER2017.pdf, last accessed on February 29, 2020).

The impressive developments in fast and affordable bacterial WGS, the availability of appropriate bioinformatics tools and comprehensive databases, make in silico antimicrobial resistance (AMR) detection a target within reach (Hendriksen et al., 2019). While many hindrances still exist (Ellington et al., 2017), the use of WGS to detect AMR genes and infer phenotypes has shown excellent results in some recent studies (Feldgarden et al., 2019; Hendriksen et al., 2019; Su et al., 2019), including Campylobacter (Zhao et al., 2016; Whitehouse et al., 2018; Painset et al., 2020).

In this study, WGS analysis of representative C. jejuni Israeli isolates retrieved from a national collection was applied to study the genetic basis and dynamics of genotypic AMR and to evaluate the role of WGS in predicting phenotypic AMR.

Materials and Methods

Study Isolates

We studied 263 C. jejuni strains representing the major lineages in Israel as identified by MLST and collected over a 10-year period (2003–2012), as previously described (Weinberger et al., 2016; Rokney et al., 2018). The majority of isolates were recovered from human clinical specimens (n = 239), and the remaining were recovered from poultry and bovine sources (Supplementary Table 1).

Antimicrobial Susceptibility Testing

A sample of 219 strains was tested for antimicrobial susceptibility to erythromycin, ciprofloxacin, tetracycline, gentamicin, nalidixic acid, and streptomycin (1,314 tests). The analysis was performed by broth microdilution at the National Reference Center for Campylobacter. Isolates were cultured on tryptic soy agar plates with 5% defibrinated sheep blood (Hylabs, Israel) and incubated at 42°C under microaerobic conditions (5% O₂, 10% CO₂) for 16–24 h. Broth microdilution was performed using the EUCAMP2 plate (Sensititre^TM, TREK diagnostic systems LTD.), according to manufacturer's instructions. Minimum inhibitory concentrations (MIC) were interpreted using the European Committee on Antimicrobial Susceptibility Testing (EUCAST 2019) breakpoint for erythromycin (R > 4 μg/mL), ciprofloxacin (R > 0.5 μg/mL), and tetracycline (R > 2 μg/mL). The National Antimicrobial Resistance Monitoring System for Enteric Bacteria (NARMS, 2018) breakpoints were applied for gentamicin (R ≥ 4 μg/mL) and for nalidixic acid (≥32 μg/ μg/mL). Streptomycin resistance was determined according to the EUCAST epidemiological cut-off values (ECOFFs) (R > 4 μg/mL). Multidrug resistance was defined as non-susceptibility to at least one agent in three or more antimicrobial classes. Broth microdilution results are detailed in Supplementary Table 2.

WGS Analysis

DNA extraction was performed on the automated QIAsymphony SP platform (Qiagen, Hilden, Germany) and quantified by the Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, CA, USA). DNA libraries were prepared from 1 ng of purified DNA using Nextera XT DNA Sample Prep Kit and Nextera XT Index Kit (Illumina, Inc., San Diego, CA, USA). Libraries were subject to short read sequencing PE 250bpx2, aiming at >100X coverage. The reads were de novo assembled by SPAdes 3.7.1 on the BioNumerics calculation engine (Applied Maths, Belgium). Whole-genome multi locus sequence typing (wgMLST) analysis of assembly-free and assembly-based allele calls for 3,529 loci according to the BioNumerics Campylobacter scheme was performed on the BioNumerics calculation engine (Applied Maths, Belgium). Assemblies are available at the Campylobacter PubMLST database with IDs 76153-76415.

Phylogeny was inferred by calculating a minimal spanning tree based on wgMLST or MLST allelic profiles. Identification of potential clusters was based on a cutoff of ≤ 15 wgMLST allele difference (Rokney et al., 2018; Schürch et al., 2018).

Assemblies were screened by the NCBI AMRFinderPlus tool (version 3.1.1b) for known acquired AMR genes and resistance-conferring point mutations in Campylobacter (available at https://ftp.ncbi.nlm.nih.gov/pathogen/Antimicrobial_resistance/AMRFinder/, last accessed on February 29, 2020) (Feldgarden et al., 2019). The BioNumerics version 7.6.3 gene extraction tool (Applied Maths, Belgium) was used to screen assemblies for the gene cluster that includes cmeR and the multidrug efflux pump cmeABC (Accession AF466820). The G → T transversion in the promoter upstream of bla_OXA−193 was screened with accession NG049489.1 as reference.

Genotypic–Phenotypic Comparisons

The WGS-predicted resistance was compared with the phenotypic resistance for six clinically relevant antimicrobial agents (erythromycin, ciprofloxacin, nalidixic acid, tetracycline, gentamicin, streptomycin). The following calculations were performed: correlation (calculated as the sum of true positive and true negative divided by all tested isolates), sensitivity (calculated by dividing true positive by the sum of true positive and false positive), specificity (calculated by dividing true negative by the sum of true negative and false negative), positive predictive value (PPV) (calculated by dividing true positive by the sum of true positive and false negative), and negative predictive value (NPV) (calculated by dividing true negative by the sum of true negative and false positive).

Statistical Methods

Poisson regression models accounting for overdispersion were used to study the annual trends of the incidence rate for Campylobacter AMR determinants. In these models, the yearly counts of isolates positive for the AMR determinants were the dependent variable; the calendar year was the independent variable, and the total number of tested isolates was the offset parameter.

Ethical Consideration

The study was approved by the local Ethical Committee at Shamir (Assaf Harofeh) Medical Center, Zerifin (#160/11), Israel and the local Ethical Committee of Tel Aviv University, Tel Aviv, Israel.

Results

WGS of Study Isolates

A total of 263 C. jejuni isolates were analyzed by whole genome sequencing (Supplementary Table 1). The strains were collected from human (n = 239) and veterinary (n = 24) sources over a decade. The strains represented the C. jejuni population in Israel and included 18 clonal complexes (CCs) and 51 sequence types (STs), including the major lineages in Israel (Figure 1). Isolate sources, lineages, and the quality metrics of reads and assemblies are detailed in Supplementary Table 1.

FIGURE 1

Figure 1. wgMLST-based phylogeny of 263 C. jejuni isolates. The minimum spanning tree is based on wgMLST analyses of 239 clinical and 24 veterinary isolates. Isolates are represented by circles connected by branches proportional to the allelic distance. The distribution of clonal complexes among the studied population is denoted by color. Partitioned nodes represent closely clustered isolates (≤ 15 allelic distance threshold).

Genomic Prediction of AMR

AMRFinderPlus analysis detected a total of 21 genetic AMR determinants among the 263 C. jejuni studied sequences. These included 18 acquired AMR genes and 3 resistance-conferring point mutations (Table 1). The tet(O) gene, predicting resistance to tetracycline (Feldgarden et al., 2019), was the most prevalent AMR gene occurring in 82% of the isolates. The gene aadE (ant(6)-Ia) predicting resistance to streptomycin (Feldgarden et al., 2019) was detected in 19% of the isolates, while the genes aph(3′)-IIIa (V01547.1) predicting resistance to kanamycin and amikacin (Feldgarden et al., 2019), aad9 (AY701528.1) predicting resistance to spectinomycin, and sat4 (AJ489618.1) predicting resistance to streptothricin (Shen et al., 2017) were only rarely found (0.8, 0.4, and 0.8%, respectively).

TABLE 1

Table 1. Antimicrobial resistance genetic determinants detected by AMRFinderPlus tool in 263 C. jejuni strains.

A variety of 12 known β lactam resistance genes (bla_OXA variants) were detected by AMRFinderPlus in 241 strains (92%), the most prevalent being bla_OXA−193, bla_OXA−461, and bla_OXA−580 (56, 16, and 7%, respectively) (Table 1). The G → T transversion in the −10 promoter region of bla_OXA−193, which confers high-level β lactam resistance in C. jejuni (Zeng et al., 2014; Lopes et al., 2019), was found in 53 (20%) isolates.

The gyrA T86I point mutation predicting reduced susceptibility to fluoroquinolones (Shen et al., 2017) was highly prevalent and was detected in 76% of the isolates by AMRFinderPlus (Table 1). In addition, 19 mutations in GyrA were also detected using the BioNumerics software. The number of mutations in a single isolate ranged from 1 to 6.

The 23S rRNA A2075G point mutation conferring high-level macrolide resistance (Corcoran et al., 2006) was found in two isolates only (0.8%). The 50S rRNA L22-A103V point mutation in the L22 protein was detected in 9.1% of the isolates. This mutation has been reported in both susceptible and resistant Campylobacter strains (Rozynek et al., 2013; Lehtopolku et al., 2020) (Table 1).

Additionally, the operon of the cmeABC multidrug efflux transporter, consisting of cmeA, cmeB, and cmeC genes, was present in all 263 isolates (Table 1). This transporter is the predominant efflux system in Campylobacter and contributes to resistance to structurally diverse agents. Inactivation of cmeABC results in susceptibility to tetracyclines and fluoroquinolones (Shen et al., 2017).

A high prevalence of the arsenic resistance genes acr3 and arsP was observed. The arsP organoarsenical efflux permease (Shen et al., 2017) was found in 182 strains (69%). The presence of arsP is associated with elevated MICs to the arsenic compound roxarsone. The acr3 gene, encoding an arsenite efflux transporter (Shen et al., 2017), was found in 82 strains (31%). Both genes coexisted in 80 isolates.

Phenotypic Determination of AMR

The phenotypic AMR to six antimicrobial agents, determined by broth microdilution, is presented in Table 2. Among the 219 tested isolates, phenotypic antimicrobial resistance rates were as follows: tetracycline, 84.5%; ciprofloxacin, 77.6%; nalidixic acid, 77.2%; streptomycin, 17.8%; and erythromycin, 0.9%. Resistance to gentamicin was not detected. Only 16 (7.3%) isolates were susceptible to all six tested antibiotic agents, while the rest were resistant to between one and five agents (Table 3). More than half of the isolates were resistant to both quinolones and tetracycline (51.6%), and 16.9% were multidrug resistant.

TABLE 2

Table 2. Minimum inhibitory concentration (MIC) distribution among 219 C. jejuni isolates^a.

TABLE 3

Table 3. Phenotypic resistance profiles determined by broth microdilution among 219 C. jejuni isolates.

Comparison Between Genomic Prediction and Phenotypic Determination of AMR

Comparison between WGS-based genotypic predictions of AMR and microdilution-based phenotypic AMR was based on 1,314 susceptibility tests for 219 isolates also screened for the presence of AMR genetic determinants. The overall correlation rate was 98.8%, with a sensitivity of 98.0%, specificity of 99.3%, PPV of 99.1%, and NPV of 98.5% (Table 4). There were 16 (1.2%) discrepancies, mostly associated with an absence of resistance determinants among phenotypically resistant isolates. The highest number of discrepancies was noted for tetracycline, with six (2.7%) of the phenotypically resistant isolates reported as negative for the tet(O) gene. Further analysis of these discrepancies by BLAST indicated that in four strains, tet(O) is present, but it is fragmented among small contigs and in different strands, and therefore is not identified by AMRFinderPlus. Therefore, we attribute 4 discrepancies to assembly fragmentation, while 2 discrepancies could not be accounted for. A combination of the assembly-based approach with read-based mapping may prove useful in minimizing discrepancies due to gene fragmentation.

TABLE 4

Table 4. Comparison between microdilution-based phenotypic resistance and genotype-predicted resistance according to AMR genetic determinant for 219 C. jejuni isolates.

The highest correlation rate was found for erythromycin and gentamicin. The two isolates with the 23S rRNA A2075G point mutation were phenotypically highly resistant to erythromycin (MIC ≥ 128 μg/ML). All isolates were susceptible to gentamicin and no resistance determinant to gentamicin was detected. The correlation rate of ciprofloxacin and nalidixic acid genotypic prediction was also high, reaching 98.2 and 98.6%, respectively. Only the gyrA T86I point mutation showed this high correlation rate, while all the other gyrA point mutations were present in both resistant and susceptible isolates (Figure 2). Notably, the gyrA T86I point mutation was present across the whole range of resistance to ciprofloxacin and nalidixic acid, and the number of the additional point mutations per isolate was not correlated with the level of resistance. Similarly, the range of resistance to tetracycline and streptomycin could not be explained by the single resistance determinant detected for each agent.

FIGURE 2

Figure 2. Correlation between quinolone MIC and the presence of point mutations in the GyrA protein among 219 C. jejuni isolates. The different mutations are denoted by color. Resistance breakpoints: ciprofloxacin (>0.5 μg/ml), nalidixic acid (>16 μg/ml).

Distribution of AMR Determinants Across Genetic Clusters

The distribution of resistance determinants among lineages and within genetic clusters was further investigated. wgMLST-based phylogeny indicated a high level of clonality and clustering among the studied population. Among 263 isolates, 132 isolates (50.2%) formed 38 genetic clusters (single linkage distance ≤ 15 alleles). Cluster sizes ranged between 2 and 13 strains (Figure 1). We investigated the distribution of AMR genes within clades and clusters. Figure 3 depicts the presence of AMR determinants among wgMLST phylogeny of CC-21, the largest clonal complex. A variety of AMR profiles were detected among CC-21 sequences. The predominant STs (ST-21, ST-50, ST-1359, ST-883) were each associated with more than one AMR profile (2–9 profiles).

FIGURE 3

Figure 3. AMR genetic determinants among a wgMLST phylogeny of Clonal Complex 21. A wgMLST-based minimum spanning tree of 104 CC-21 isolates is shown. Each node represents a strain. Partitioned nodes represent closely clustered isolates (≤ 15 allelic distance threshold). Sequence types are denoted in black. The number of allelic differences is shown on the branches connecting the nodes. The AMR gene profile detected by AMRFinderPlus is denoted by color.

WGS-based clustering was predictive of AMR gene content. Closely related isolates that were part of a genetic cluster also shared a homogenous AMR determinant profile. This was observed in 18 of 20 (90.0%) clusters in CC-21. These results may reflect clonal expansion of resistant isolates. Phylogenetic analysis reflected varied patterns of antimicrobial gene distribution across the studied isolates. The prevalent tet(O) gene and the gyrA T86I point mutation were evenly distributed across the wgMLST clusters (Figures 4, 5). The gyrA T86I mutation was less prevalent in some clusters within CC-257 and CC-206. Strong lineage association was noted for the aadE (Figure 6) gene and the various bla_OXA genes (Figure 7). For example, the bla_OXA−580 gene was uniquely found in ST-50 and the bla_OXA−461 was predominantly found in ST-257.

FIGURE 4

Figure 4. wgMLST-based phylogeny of 263 C. jejuni isolates. The distribution of the tet(O) gene associated with tetracycline resistance is shown in color.

FIGURE 5

Figure 5. wgMLST-based phylogeny of 263 C. jejuni isolates. The distribution of the gyrA T86I point mutation associated with quinolone resistance is shown in color.

FIGURE 6

Figure 6. wgMLST-based phylogeny of 263 C. jejuni isolates. The distribution of the aadE (ant(6)-Ia) gene associated with streptomycin resistance is shown in color.

FIGURE 7

Figure 7. The distribution of bla_OXA gene variants across an MLST-based phylogeny of 263 C. jejuni isolates. Each node in the minimum spanning tree represents a sequence type (ST), and the number of allelic differences is denoted on the branches connecting the nodes. The presence of β-lactamase (bla_OXA) gene variants is denoted by color. The G → T promoter transversion associated with high-level ampicillin resistance in bla_OXA−193 is shown in black.

The G → T transversion in bla_OXA−193 was detected mainly in ST-4057 and ST-1359 strains (Figure 7). The aminoglycoside resistance gene cluster aadE (ant(6)-Ia- aph(3')-IIIa-sat4 was detected in two isolates of ST-2337 from bacteremic patients.

Temporal Trends in the Occurrence of AMR Determinants

We investigated trends in the presence of AMR determinants among 239 human C. jejuni studied isolates over a 10-year period (2003–2012). The prevalence of gyrA T86I point mutation increased from 62.5% in 2003 to 91.3% in 2012 (a 1.5-fold increase). This trend differed among the various lineages. Among CC-21 isolates, the prominent clonal complex in Israel, gyrA T86I detection rate was high and stable along the study period, and as high as 100% in several years, while in non-CC-21 lineages detection increased 1.7-fold (from 50.0% in 2003 to 88.2% in 2012) (Figure 8). The detection of tet(O) gene increased 1.4-fold during the study period from 68.8 to 95.7%; no substantial difference was noted between CC-21 and the other CCs (Figure 9). These trends did not reach statistical significance and warrant further investigation on a larger scale.

FIGURE 8

Figure 8. Trends in the prevalence of the gyrA T86I point mutation among 239 human C. jejuni isolates collected over a decade.

FIGURE 9

Figure 9. Trends in the prevalence of the tet(O) gene among 239 human C. jejuni isolates collected over a decade.

Discussion

This study explored resistance trends among Israeli C. jejuni over more than a decade, the genetic basis for resistance phenotypes, and the association of phenotypic and genotypic resistance with phylogeny. The study indicates an alarmingly high rate of quinolone and tetracycline resistance among the Israeli C. jejuni strains emerging over the last decade, and a low resistance rate to macrolides and aminoglycosides. These rates are higher than those reported by NARMS for the United States and the average rates reported by the EU, but comparable with the rates in several EU countries such as Spain, Portugal, Lithuania, or Cyprus (CDC, 2018; European Food Safety Authority and European Centre for Disease Prevention and Control, 2019). Current resistance rates at certain localities in Israel may be even higher (Azrad et al., 2018).

We found a high correlation rate between the presence of AMR determinants detected by the NCBI AMRFinderPlus tool and the phenotypic resistance. The overall correlation rate was 98.8%, with a sensitivity of 98.0%, specificity of 99.3%, PPV of 99.1%, and NPV of 98.5%. Comparable results were reported by Feldgarden et al. for C. jejuni, using the NCBI AMRFinder tool, e.g., 98.9% correlation rate, with a PPV of 97.1% and NPV of 99.2% (Feldgarden et al., 2019). A high correlation rate of 97.5% was also found in a most recent study from England and Wales using another AMR extraction tool, GeneFinder (Painset et al., 2020).

The overall discordance rates between genetic prediction and phenotypic AMR for C. jejuni was low in our study (1.2%), as well as in the above two studies. The highest discordance rate (7 out of 219, 3.2%) in our study was noted for tetracycline, and particularly for resistant isolates lacking the tet(O) gene (6/219, 2.7%). No discordance was noted for gentamicin and erythromycin. The discordance rate for the other antibiotic agents (ciprofloxacin, tetracycline, streptomycin) was <2% (Table 4). Feldgarden et al. and Painset et al. also found 100% consistency for erythromycin based on the presence/absence of mutations in the 23S rRNA, and a low rate of genotype-phenotype discordances for fluoroquinolones, macrolides, and aminoglycosides (Feldgarden et al., 2019; Painset et al., 2020). An exception was resistance to streptomycin in the study by Painset et al., as only 2 of 5 isolates carrying the streptomycin resistant gene (aadE) were phenotypically resistant to the drug. In our study, 36 isolates resistant to streptomycin carried the aadE resistance gene, but it was lacking in two other resistant phenotypes.

There are multiple explanations for possible WGS-based genotypic-phenotypic discrepancies, related to the quality of assembly, bioinformatics tools, the completeness, and curation of WGS-based AMR database, and factors related to the studied isolate, namely gene content, transcriptome and proteome, and plasmid stability during microbiological procedures (Zankari et al., 2017; Chen et al., 2019; Doyle et al., 2020). Low quality of assemblies may lead to defective genes, resulting in incorrect extraction from the AMR databases (Feldgarden et al., 2019; Doyle et al., 2020). Long-read sequencing technology can improve complete genome assembly accuracy with the addition of plasmid sequencing (Chen et al., 2019). Indeed, the current short-read assembly technologies used for AMR detection may miss genetic determinants carried on plasmids, mainly due to repetition or redundancy (Collineau et al., 2019).

The performance of web-based AMR search tools is also important for the prediction of phenotypes. For example, a comparison between AMRFinder and ResFinder (Feldgarden et al., 2019) showed miscalling of aminoglycosides resistance genes due to lack of correspondence between the closest nucleotide hit and actual observed sequence. In other instances ResFinder overspecified, while AMRfinder underspecified the gene symbols, resulting in decreased genotype-phenotype correlation (Feldgarden et al., 2019). Another study took it one step further and compared the results of AMR phenotype prediction of the same bacterial sequences among nine different laboratories (Doyle et al., 2020). It showed a high degree of discordance of genotype-phenotype prediction and substantial variability between the participating laboratories, which stemmed from multiple factors, including the quality of the sequence data, the choice of bioinformatic pipeline, and the interpretation of the results. This study highlighted the understanding that genotypic AMR determination requires standardization of sequencing methods, the AMR databases, and the bioinformatics tools.

At the pathogen level many factors may influence lack of AMR genotype-phenotype correlation, such as the number of alleles, overexpression and under expression of genes (e.g., efflux pumps), mutations in regulatory proteins or their binding sites, protein abundance levels, and contributing factors that are not covered by the searched database (e.g., porin mutations). Finally, novel resistant genes and mutations may be missed by the current databases and search tools (Zankari et al., 2017; Chen et al., 2019; Feldgarden et al., 2019). This may also explain the lack of correlation between the presence of a single resistance determinant and a range of MICs for a specific antibiotic.

The trend of increasing prevalence of resistance determinants to quinolones and tetracycline along the study years is of concern. Similar trends have been reported globally and linked to extensive use of clinically important antibiotics in veterinary medicine (Takahashi et al., 2005; Tang et al., 2017). Monitoring these trends is important for decisions on the choice of appropriate treatment of clinical infections as well as intervention policies. These results also call for a better control of antibiotic-resistant C. jejuni in Israel.

The wgMLST phylogenetic analysis indicates that the C. jejuni population in Israel consists of diverse lineages and genotypes. Notably, half of the strains were part of genetic clusters (≤ 15 allelic distance), indicative of the dominance of several different clones within the population. Clonality may reflect common source or diffuse outbreaks as well as selective expansion of successful clones over time and space. Our analysis shows a link between phylogeny and AMR profile. Among CC-21, closely clustered isolates contained uniform AMR arsenals (18 of 20 clusters). Therefore, the AMR profile could be correctly predicted based on wgMLST clustering in 90% of clusters.

Our study is limited by the small sample size of isolates selected to represent a decade of surveillance. More recent isolates may exhibit a higher rate of resistance, but this is not expected to influence the predictive value of genotypic-phenotypic correlations.

In conclusion, our study illuminates the genetic mechanisms underlying AMR in C. jejuni in Israel. We report a high rate of tetracycline and quinolone resistance among the Israeli C. jejuni strains and a low resistance rate to macrolides and aminoglycosides. Some AMR determinants were associated with distinct genotypes and clones. Our study revealed a high correlation rate between the presence of AMR determinants found by the NCBI AMRFinderPlus tool and phenotypic resistance. However, some resistance phenotypes could not be predicted by the presence of AMR determinants and particularly not the level of resistance expressed as MIC. Further work to improve phenotypic prediction of AMR is central to the envisaged transition toward WGS analysis as the primary approach for characterizing foodborne isolates.

Data Availability Statement

The datasets generated for this study can be found in the Campylobacter PubMLST database, IDs 76153-76415.

Ethics Statement

The studies involving human participants were reviewed and approved by the local Ethical Committee at Shamir (Assaf Harofeh) Medical Center, Zerifin (#160/11), Israel and the local Ethical Committee of Tel Aviv University, Tel Aviv, Israel. Written informed consent from the participants' legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.

Author Contributions

AR, LV, JM-G, VA, and MW conceived the study. AR, LV, and NF performed the laboratory analysis. AR performed the bioinformatic analysis for antibiotic resistance determinants. AR, KV, and MW performed the data analyses. AR and MW drafted the manuscript with help from the other authors. JM-G critically revised the data and the manuscript. All authors read and approved of the submitted version.

Funding

This study was supported by Grant No. 1416/12 by the Israel Science Foundation. The laboratory work of this study was performed at the Public Health Laboratory, Jerusalem, Israel Ministry of Health and supported by it.

Conflict of Interest

KV was employed by Applied Maths NV.

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

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2020.00365/full#supplementary-material

Supplementary Table 1. List of isolates included in the study, metadata and sequence parameters.

Supplementary Table 2. Genotypic antimicrobial resistance determinants and phenotypic broth microdilution results.

References

Azrad, M., Tkhawkho, L., Isakovich, N., Nitzan, O., and Peretz, A. (2018). Antimicrobial susceptibility of Campylobacter jejuni and Campylobacter coli: comparison between etest and a broth dilution method. Ann. Clin. Microbiol. Antimicrob. 17:23. doi: 10.1186/s12941-018-0275-8

PubMed Abstract | CrossRef Full Text | Google Scholar

CDC (2018). National Antimicrobial Resistance Monitoring System for Enteric Bacteria (NARMS): Human Isolates Surveillance Report for 2015 (Final Report). Atlanta, GA: U.S. Department of Health and Human Services, CDC. Available online at: https://www.cdc.gov/narms/pdf/2015-NARMS-Annual-Report-cleared_508.pdf (accessed July 18, 2020).

CDC (2019). Antibiotic Resistance Threats in the United States, 2019. Atlanta, GA: U.S. Department of Health and Human Services, CDC. Available online at: https://www.cdc.gov/drugresistance/biggest-threats.html (accessed July 18, 2020).

Google Scholar

Chen, C. Y., Clark, C. G., Langner, S., Boyd, D. A., Bharat, A., McCorrister, S. J., et al. (2019). Detection of antimicrobial resistance using proteomics and the comprehensive antibiotic resistance database: a case study. Proteomics Clin. Appl. 14:e1800182. doi: 10.1002/prca.201800182

PubMed Abstract | CrossRef Full Text | Google Scholar

Collineau, L., Boerlin, P., Carson, C. A., Chapman, B., Fazil, A., Hetman, B., et al. (2019). Integrating whole-genome sequencing data into quantitative risk assessment of foodborne antimicrobial resistance: a review of opportunities and challenges. Front. Microbiol. 10:1107. doi: 10.3389/fmicb.2019.01107

PubMed Abstract | CrossRef Full Text | Google Scholar

Corcoran, D., Quinn, T., Cotter, L., and Fanning, S. (2006). An investigation of the molecular mechanisms contributing to high-level erythromycin resistance in campylobacter. Int. J. Antimicrob. Agents 27, 40–45. doi: 10.1016/j.ijantimicag.2005.08.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Crowe, S. J., Green, A., Hernandez, K., Peralta, V., Bottichio, L., Defibaugh-Chavez, S., et al. (2017). Utility of combining whole genome sequencing with traditional investigational methods to solve foodborne outbreaks of salmonella infections associated with chicken: a new tool for tackling this challenging food vehicle. J. Food Prot. 80, 654–660. doi: 10.4315/0362-028X.JFP-16-364

PubMed Abstract | CrossRef Full Text | Google Scholar

Doyle, R. M., O'Sullivan, D. M., Aller, S. D., Bruchmann, S., Clark, T., Coello Pelegrin, A., et al. (2020). Discordant bioinformatic predictions of antimicrobial resistance from whole-genome sequencing data of bacterial isolates: an inter-laboratory study. Microb. Genom. 6:e000335. doi: 10.1099/mgen.0.000335

PubMed Abstract | CrossRef Full Text | Google Scholar

Ellington, M. J., Ekelund, O., Aarestrup, F. M., Canton, R., Doumith, M., Giske, C., et al. (2017). The role of whole genome sequencing in antimicrobial susceptibility testing of bacteria: report from the EUCAST Subcommittee. Clin. Microbiol. Infect. 23, 2–22. doi: 10.1016/j.cmi.2016.11.012

PubMed Abstract | CrossRef Full Text | Google Scholar

European Food Safety Authority and European Centre for Disease Prevention and Control (2019). The European Union summary report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2017. EFSA J. 17:e05598. doi: 10.2903/j.efsa.2019.5598

CrossRef Full Text | Google Scholar

Feldgarden, M., Brover, V., Haft, D. H., Prasad, A. B., Slotta, D. J., Tolstoy, I., et al. (2019). Validating the AMRFinder tool and resistance gene database by using antimicrobial resistance genotype-phenotype correlations in a collection of isolates. Antimicrob. Agents Chemother. 63:e00483-19. doi: 10.1128/AAC.00483-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Florez-Cuadrado, D., Moreno, M. A., Ugarte-Ruiz, M., and Dominguez, L. (2018). Antimicrobial resistance in the food chain in the European Union. Adv. Food Nutr. Res. 86, 115–136. doi: 10.1016/bs.afnr.2018.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerner-Smidt, P., Besser, J., Concepcion-Acevedo, J., Folster, J. P., Huffman, J., Joseph, L. A., et al. (2019). Whole genome sequencing: bridging one-health surveillance of foodborne diseases. Front. Public Health 7:172. doi: 10.3389/fpubh.2019.00172

CrossRef Full Text | Google Scholar

Hendriksen, R. S., Bortolaia, V., Tate, H., Tyson, G. H., Aarestrup, F. M., and McDermott, P. F. (2019). Using genomics to track global antimicrobial resistance. Front. Public Health 7:242. doi: 10.3389/fpubh.2019.00242

PubMed Abstract | CrossRef Full Text | Google Scholar

Inns, T., Ashton, P. M., Herrera-Leon, S., Lighthill, J., Foulkes, S., Jombart, T., et al. (2017). Prospective use of whole genome sequencing (WGS) detected a multi-country outbreak of Salmonella Enteritidis. Epidemiol. Infect. 145, 289–298. doi: 10.1017/S0950268816001941

PubMed Abstract | CrossRef Full Text | Google Scholar

Jagadeesan, B., Gerner-Smidt, P., Allard, M. W., Leuillet, S., Winkler, A., Xiao, Y., et al. (2019). The use of next generation sequencing for improving food safety: translation into practice. Food Microbiol. 79, 96–115. doi: 10.1016/j.fm.2018.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubota, K. A., Wolfgang, W. J., Baker, D. J., Boxrud, D., Turner, L., Trees, E., et al. (2019). PulseNet and the changing paradigm of laboratory-based surveillance for foodborne diseases. Public Health Rep. 134(Suppl. 2), 22S−28S. doi: 10.1177/0033354919881650

PubMed Abstract | CrossRef Full Text | Google Scholar

Lehtopolku, M., Kotilainen, P., Haanpera-Heikkinen, M., Nakari, U. M., Hanninen, M. L., Huovinen, P., et al. (2020). Ribosomal mutations as the main cause of macrolide resistance in Campylobacter jejuni and Campylobacter coli. Antimicrob. Agents Chemother. 55, 5939–5941. doi: 10.1128/AAC.00314-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Llarena, A. K., Taboada, E., and Rossi, M. (2017). Whole-genome sequencing in epidemiology of Campylobacter jejuni infections. J. Clin. Microbiol. 55, 1269–1275. doi: 10.1128/JCM.00017-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopes, B. S., Strachan, N. J. C., Ramjee, M., Thomson, A., MacRae, M., Shaw, S., et al. (2019). Nationwide sepwise emergence and evolution of multidrug-resistant Campylobacter jejuni sequence type 5136, United Kingdom. Emerg. Infect. Dis. 25, 1320–1329. doi: 10.3201/eid2507.181572

PubMed Abstract | CrossRef Full Text | Google Scholar

Moran-Gilad, J. (2017). Whole genome sequencing (WGS) for food-borne pathogen surveillance and control-taking the pulse. Euro Surveill. 22:30547. doi: 10.2807/1560-7917.ES.2017.22.23.30547

PubMed Abstract | CrossRef Full Text | Google Scholar

Motro, Y., and Moran-Gilad, J. (2017). Next-generation sequencing applications in clinical bacteriology. Biomol. Detect. Quantif. 14, 1–6. doi: 10.1016/j.bdq.2017.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Nadon, C., Van Walle, I., Gerner-Smidt, P., Campos, J., Chinen, I., Concepcion-Acevedo, J., et al. (2017). PulseNet international: vision for the implementation of whole genome sequencing (WGS) for global food-borne disease surveillance. Euro Surveill. 22:30544. doi: 10.2807/1560-7917.ES.2017.22.23.30544

PubMed Abstract | CrossRef Full Text | Google Scholar

Painset, A., Day, M., Doumith, M., Rigby, J., Jenkins, C., Grant, K., et al. (2020). Comparison of phenotypic and WGS-derived antimicrobial resistance profiles of Campylobacter jejuni and Campylobacter coli isolated from cases of diarrhoeal disease in England and Wales, 2015–16. J. Antimicrob. Chemother. 75, 883–889. doi: 10.1093/jac/dkz539

PubMed Abstract | CrossRef Full Text | Google Scholar

Revez, J., Espinosa, L., Albiger, B., Leitmeyer, K. C., Struelens, M. J., Points, E. N. M. F., et al. (2017). Survey on the use of whole-genome sequencing for infectious diseases surveillance: rapid expansion of European national capacities, 2015–2016. Front. Public Health 5:347. doi: 10.3389/fpubh.2017.00347

PubMed Abstract | CrossRef Full Text | Google Scholar

Ribot, E. M., Freeman, M., Hise, K. B., and Gerner-Smidt, P. (2019). PulseNet: entering the age of next-generation sequencing. Foodborne Pathog. Dis. 16, 451–456. doi: 10.1089/fpd.2019.2634

PubMed Abstract | CrossRef Full Text | Google Scholar

Rokney, A., Valinsky, L., Moran-Gilad, J., Vranckx, K., Agmon, V., and Weinberger, M. (2018). Genomic epidemiology of Campylobacter jejuni transmission in Israel. Front. Microbiol. 9:2432. doi: 10.3389/fmicb.2018.02432

PubMed Abstract | CrossRef Full Text | Google Scholar

Rozynek, E., Mackiw, E., Kaminska, W., Tomczuk, K., Antos-Bielska, M., Dzierzanowska-Fangrat, K., et al. (2013). Emergence of macrolide-resistant Campylobacter strains in chicken meat in Poland and the resistance mechanisms involved. Foodborne Pathog. Dis. 10, 655–660. doi: 10.1089/fpd.2012.1333

PubMed Abstract | CrossRef Full Text | Google Scholar

Rumore, J., Tschetter, L., Kearney, A., Kandar, R., McCormick, R., Walker, M., et al. (2018). Evaluation of whole-genome sequencing for outbreak detection of verotoxigenic Escherichia coli O157:H7 from the Canadian perspective. BMC Genomics 19:870. doi: 10.1186/s12864-018-5243-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Schürch, A. C., Arredondo-Alonso, S., Willems, R. J. L., and Goering, R. V. (2018). Whole genome sequencing options for bacterial strain typing and epidemiologic analysis based on single nucleotide polymorphism versus gene-by-gene-based approaches. Clin. Microbiol. Infect. 24, 350–354. doi: 10.1016/j.cmi.2017.12.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Z., Wang, Y., Zhang, Q., and Shen, J. (2017). Antimicrobial resistance in Campylobacter spp. Microbiol. Spectrum 6:ARBA-0013-2017. doi: 10.1128/microbiolspec.ARBA-0013-2017

CrossRef Full Text | Google Scholar

Sproston, E. L., Wimalarathna, H. M. L., and Sheppard, S. K. (2018). Trends in fluoroquinolone resistance in Campylobacter. Microb. Genom. 4, 1–8. doi: 10.1099/mgen.0.000198

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, M., Satola, S. W., and Read, T. D. (2019). Genome-based prediction of bacterial antibiotic resistance. J. Clin. Microbiol. 57:e01405-18. doi: 10.1128/JCM.01405-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Takahashi, T., Ishihara, K., Kojima, A., Asai, T., Harada, K., and Tamura, Y. (2005). Emergence of fluoroquinolone resistance in Campylobacter jejuni in chickens exposed to enrofloxacin treatment at the inherent dosage licensed in Japan. J. Vet. Med. B Infect. Dis. Vet. Public Health 52, 460–464. doi: 10.1111/j.1439-0450.2005.00890.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, Y., Fang, L., Xu, C., and Zhang, Q. (2017). Antibiotic resistance trends and mechanisms in the foodborne pathogen, Campylobacter. Anim. Health. Res. Rev. 18, 87–98. doi: 10.1017/S1466252317000135

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinberger, M., Lerner, L., Valinsky, L., Moran-Gilad, J., Nissan, I., Agmon, V., et al. (2013). Increased incidence of Campylobacter spp. infection and high rates among children, Israel. Emerg. Infect. Dis. 19, 1828–1831. doi: 10.3201/eid1911.120900

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinberger, M., Moran-Gilad, J., Rokney, A., Davidov, Y., Agmon, V., Peretz, C., et al. (2016). Molecular epidemiology of Campylobacter jejuni infection in Israel-a nationwide study. Clin. Microbiol. Infect. 22:1005.e9–1005.e15. doi: 10.1016/j.cmi.2016.08.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Whitehouse, C. A., Young, S., Li, C., Hsu, C. H., Martin, G., and Zhao, S. (2018). Use of whole-genome sequencing for Campylobacter surveillance from NARMS retail poultry in the United States in 2015. Food Microbiol. 73, 122–128. doi: 10.1016/j.fm.2018.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Zankari, E., Allesoe, R., Joensen, K. G., Cavaco, L. M., Lund, O., and Aarestrup, F. M. (2017). PointFinder: a novel web tool for WGS-based detection of antimicrobial resistance associated with chromosomal point mutations in bacterial pathogens. J. Antimicrob. Chemother. 72, 2764–2768. doi: 10.1093/jac/dkx217

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, X., Brown, S., Gillespie, B., and Lin, J. (2014). A single nucleotide in the promoter region modulates the expression of the beta-lactamase OXA-61 in Campylobacter jejuni. J. Antimicrob. Chemother. 69, 1215–1223. doi: 10.1093/jac/dkt515

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, S., Tyson, G. H., Chen, Y., Li, C., Mukherjee, S., Young, S., et al. (2016). Whole-genome sequencing analysis accurately predicts antimicrobial resistance phenotypes in Campylobacter spp. Appl. Environ. Microbiol. 82, 459–466. doi: 10.1128/AEM.02873-15

PubMed Abstract | CrossRef Full Text | Google Scholar