1. Introduction
Antimicrobial resistance (AMR) threatens the successful treatment of bacterial infections not only in human but also in animal health [
1,
2]. The use of antimicrobials (AB) in humans and animals is a driver in the increase of AMR in bacterial populations, even following guidelines for prudent use of AB [
3,
4]. This risk significantly increases with the misuse of these drugs [
5]. Moreover, the AMR reservoir of bacteria from livestock has been increasingly investigated for its potential to transfer AMR to humans via direct contact, the environment and contaminated food [
6]. Nevertheless, the extent of this transmission remains uncertain due to the enormous complexity of the AMR epidemiology involving animals, environment, and humans [
7,
8,
9]. Nevertheless, policy makers, in the European Union (EU), have developed legislation to monitor and regulate the antimicrobial use in animals with the goal to decrease AMR burden in humans in the long run [
10,
11]. However, the global effect of these actions, regarding the reduction of AMR at the human-animal-environment interface, is still under investigation, and very few scientific studies have shown encouraging results, limited to some antimicrobials such as colistin [
12,
13] probable due that AMU is one key driver for AMR but other socio-economic factors should be also taken into account in AMR epidemiology as recently assessed [
14]. On the other hand, this long-term reduction of AB consumption in veterinary medicine could seriously hamper the care of animals and generate severe welfare issues if animals are not treated with the right antimicrobial when it is really needed.
The current EU legislation regarding antimicrobials [
10] have focused special attention to restrict as much as possible the use of last resource antimicrobials (3rd and 4th generation cephalosporins, polymyxins and quinolones) in animals following the recommendations addressed by the European Medicine Agency in 2019 [
15]. Thus, these last resource AB can only be used when no other options belonging to less risky categories (C and D) for AMR are available to treat animals [
15]. However, up to date, most of the long-term surveillance data available are only from healthy animals that may not reflect the situation in veterinary bacterial pathogens [
16]. Thus, the European Food Safety Authority (EFSA) coordinates a mandatory active monitoring of AMR in zoonotic (
Salmonella spp and
Campylobacter spp), indicator bacteria (
Escherichia coli) and extended-spectrum-cephalosporin-resistant and carbapenemase-producing
E. coli from healthy food-producing animals (cattle, poultry, pigs) at slaughter and meat following European directives [
17,
18]. On the other hand, a coordinated and harmonized strategy for AMR monitoring in diseased animals has just started at European level [
19] to fill the gap for AMR data in pathogens from diseased animals. Thus, updated information will be generated to guide antimicrobial stewardship initiatives such as treatment guidelines, and to guide policymakers in regulating veterinary antibiotic use [
20].
The use of antimicrobials with therapeutic or metaphylactic purpose in pigs may be necessary to control the relevant pathogens involved in respiratory and enteric disorders, contributing to most of the pig antimicrobial consumption [
21,
22,
23]. Thus, porcine respiratory disease complex (PRDC) and post-weaning diarrhoea (PWD) are some of the most challenging diseases affecting the pig industry worldwide [
24,
25]. PRDC is a syndrome that results from a combination of infectious (bacteria and viruses) and non-infectious factors.
Actinobacillus pleuropneumoniae (APP),
Pasteurella multocida,
Mycoplasma hyopneumoniae, and
Bordetella bronchiseptica are the most common bacterial agents involved [
26]. On the other hand,
Escherichia coli is the main causative agent of PWD, affecting piglets after weaning. PWD is characterized by a profuse diarrhea, dehydration, significant mortality, and loss of body weight in surviving pigs [
27,
28,
29]. When clinical signs appear, prescription of antimicrobials is in many cases the only solution to control the spread of the PRDC and PWD within the herd [
21,
22,
29,
30,
31]. Thus, it may be necessary to use last resource antimicrobials if no other option is available according to an antimicrobial stewardship program [
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32]. It must be highlighted that during the last four years, the sales of last resource antimicrobials in European´s livestock are between 0.2 and 2.8% of the total sales of antimicrobials [
33], suggesting that the bacterial populations are hardly exposed to these family of drugs across Europe.
An important aspect of dealing with the AMR crisis is surveillance [
34], which provides susceptibility data allowing to act more effectively when necessary. Another goal of AMR surveillance is to analyze the temporal trends of AMR patterns for early warning of potential threats and decipher the impact of policies in animals regarding the use of antimicrobials in the long term. Unfortunately, there is scarce knowledge on the antimicrobial susceptibility profiles of veterinary bacterial pathogens in Europe due to a lack of coordinated strategy between member states [
35]. The objective of this study was to describe and analyze the temporal trends during the last four years of last resource antimicrobials in Spanish porcine pathogens as a suitable model for other countries, considering the low consumption of these drugs in Spain compared with the total antimicrobial consumption (3-4,1%) and the consistent decrease in the total antimicrobial use in livestock [
33].
3. Discussion
Antimicrobial susceptibility is usually measured by the minimum inhibitory concentration (MIC), which is the lowest concentration that stops in vitro growth of the targeted bacteria using microdilution methods in veterinary laboratories. Modelling the MIC values is challenging since these types of data are interval-censored and ordinal [
36,
37]. One approach to deal with these data is to dichotomize the MIC values into two categories, resistant (R) and susceptible (S) using established clinical breakpoints or epidemiological cut-off values (ECOFF), followed by logistic regression [
38,
39]. However, this is not an ideal approach since there is a loss of quantitative information from the MIC values when they are dichotomized [
36,
40]. Other critical point to dichotomize the MIC values into R and S categories is the existence of accepted clinical breakpoints to obtain comparable results between different studies. In the case of pig respiratory pathogens, there are a reasonable amount of internationally accepted clinical breakpoints, but this is not the case for pig enteric pathogens. Moreover, EUCAST ECOFFs are missing for 45.3% (MIC) and 76.9% (disk diffusion) of bacterial species in the veterinary field [
41]. Since we work with clinical cases, it was decided to interpret our MIC results using clinical breakpoints instead of ECOFFs. Therefore, we can monitor the antimicrobial susceptibility pattern for different antibiotics, but we cannot monitor resistance in bacterial populations as suggested by the EARS-VET surveillance network [
41]. Moreover, our study is based on clinical cases (passive collection) whose representativeness of the general animal population is unknown [
42]. Considering the limited information available for some antibiotic-microorganism pair, we have extrapolated clinical breakpoints available for quinolones and cephalosporins and respiratory pathogens [
43,
44] to enteric ones, and we have used the clinical breakpoint for colistin and
E. coli from humans [
45]. This approach seems reasonable to study the antimicrobial susceptibility temporal trends for all the porcine pathogens, but it has not allowed extrapolating directly these findings to clinical efficacy in pigs, especially for digestive pathogens. Despite these limitations, we consider that our data provide robust information about the evolution of the antimicrobial susceptibility pattern of the main pig pathogens in Spain during the study period.
The qualitative categorization into S and R, does not allow to determine dynamics of bacterial population, in particular wild type populations approaching the clinical breakpoint. This is especially important for cases of decrease susceptibility to antimicrobials associated to punctual mutations, like fluoroquinolones and
E. coli, where increase in the MIC is associated with chromosomic mutations in the quinolone resistance determining regions [
46]. MIC outcome data could be more appropriately modelled using statistical models other than logistic regression such as Cox proportional hazards, multinomial logistic, ordinal logistic, linear and tobit regression models [
36,
37,
38,
40,
47]. In this case, we have used a multinominal logistic model based on distributing the range of MIC values into four categories (from 1 to 4) that include two MIC values in each category, being category 1 the most susceptible (lowest MIC value) and 4 the less susceptible (highest MIC value), as suggested by other authors with a similar database for
E. coli [
48]. Finally, the antimicrobial panel was selected to represent commonly used compounds for the treatment of pig diseases in practice [
31,
32], and not focused on monitoring antimicrobial resistance in surveillance programs. This is a clear limitation of this study since antimicrobials tested herein were not the same for all the porcine pathogens.
Our data clearly showed a different pattern in the evolution of antimicrobial susceptibility for each combination of drug and microorganism. However, in the case of both fluoroquinolones, marbofloxacin and enrofloxacin, in combination with
A. pleuropneumoniae, the proportion of isolates susceptible to each of the antimicrobials was practically the same. Similarly occurred for
P. multocida, indicating that testing one of those fluoroquinolones in these two pathogens would be sufficient to test for this antimicrobial family [
32]. Contrarily, data on susceptibility obtained for
E. coli in combination with ceftiofur, could not be extrapolated to cefquinome as it has been also previously suggested by other authors [
49]. This is not surprising as cefquinome has been reported not useful in separating isolates with extended spectrum betalactamases or plasmidic AmpC from the cephalosporin-susceptible isolates [
50]. These results reinforced that the evolution of antimicrobial susceptibility must be studied in a case-by-case situation where generalization for drug families and bacteria is not possible as described previously [
32]. Finally, one interesting line of research could be studying the evolution mechanisms shaping the maintenance of antibiotic resistance in pig pathogens as carried out by Durao et al [
51] but it is out of the scope of this paper.
In general terms, pig pathogens involved in respiratory diseases analyzed herein appeared to remain susceptible or tended to increase susceptibility to critical antimicrobials over the study period. For
E. coli, there was also a tendency to increase susceptibility for most antimicrobials, except for ceftiofur, where there was a significant decrease in susceptibility for MIC category 1 from 2019 to 2020. Taken together, results obtained using dichotomized versus categorized MIC data were generally similar for all the pairs of drug/microorganism combinations with some exceptions, where categorized MIC was more sensitive detecting slight changes in antimicrobial susceptibility patterns (i.e. cefquinome and marbofloxacin in combination with
E. coli). Finally, for the combination colistin with
E. coli, by using dichotomized MIC data, a dramatic increase in susceptibility to colistin from 2019 to 2021 was observed, with slight decrease in 2022. This is interesting since there was a voluntary reduction in the sales of colistin in pig production in Spain from 34.9 mg/PCU to 3 mg/PCU between 2015 and 2018, which could explain these results, but we do not have figures of colistin consumption by farm and a sound study linking consumption with antimicrobial susceptibility cannot be carried out with our database. Still, by using dichotomized MIC data (S and R,
Figure 8) this decrease in susceptibility observed for the year 2022 was not detected, suggesting that categorized MIC data may be more sensible that dichotomized to detect slight changes in antimicrobial susceptibility pattern. Despite not achieving enough sample size to have robust data, the evolution of antimicrobial susceptibility for
Salmonella spp and colistin is very close to the observed tendency for
E. coli and colistin (additional
Figure S1). It must be highlighted that both bacteria are in the same ecological niche.
In Spain, the antimicrobial susceptibility for last resource antimicrobials in pig pathogens remained stable or increased in the last four years. These are sound results in terms of preserving the efficacy of critical important antimicrobials and minimizing the burden and spread of resistance from farm to fork.
4. Materials and Methods
4.1. Clinical samples
Between January 2019 and December 2022, samples were taken from diseased or recently deceased pigs from farms across Spain showing acute clinical signs of respiratory tract infections or pigs showing diarrhea. None of these animals had been exposed to antimicrobial treatment for, at least, 15 days prior sampling. Thus, the sampled animals were between 3 and 24 weeks old showing overt respiratory symptoms with or without depression and/or hyperthermia (>39.8ºC). For each clinical case, samples of lungs of two recently deceased pigs (<12 hours) were submitted under refrigeration to the laboratory. If no recently dead pigs were suitable for sampling, at least, two animals with acute respiratory signs were humanely sacrificed and lung samples were drawn. On the other hand, for piglets showing PWD, the sampled animals were between 3 and 12 weeks old showing clinical symptoms of the disease. Intestinal content obtained from humanly euthanized animals or watery diarrhea from sick pigs were obtained. In both cases, the samples were submitted under refrigeration to the laboratory and processed during the following 24 hours after collection. Only one isolate was included by farm across the study to avoid redundancy and overrepresentation of bacterial clones.
4.2. Bacterial isolation and identification
Clinical specimens were cultured aseptically onto blood agar (Columbia agar with 5% Sheep blood, 254005 BD), chocolate agar (GC II agar with IsoVitaleX, 254060, BD or blood Agar No. 2 Base, 257011, BD) and MacConkey agar (4016702, Biolife Italiana Srl) and incubated at 35–37°C in aerobic conditions with 5-10% CO2 for 24–48 hours to address the isolation of respiratory bacterial pathogens. Finally, for the isolation of digestive pathogens, specimens were cultured aseptically onto blood agar, MacConkey agar and Xylose-Lysine-Desoxycholate Agar (XLD, CM0469, Oxoid). The plates were incubated at 35–37°C in aerobic conditions for 24 hours.
Identification of isolates for respiratory pathogens and enteric pathogens was carried out by matrix assisted laser desorption ionization-time of flight (MALDI-TOF Biotyper System, Bruker Daltonics, Bremen, Germany) as previously described (25). Individual isolates were stored at -80°C in brain heart infusion (CM1135, Oxoid) with 30% of glycerol (G9012, Sigma-aldrich).
4.3. Antimicrobial susceptibility testing
Antimicrobial susceptibility testing was determined using minimum inhibitory concentration (MIC) value for each combination of bacterial species and antimicrobial tested. Thus, MIC was performed in accordance with the recommendations presented by the Clinical and Laboratory Standards Institute [
31,
32] in a customized 96-well microtitre plate (Sensititre, Trek diagnostic Systems Inc., East Grinstead, UK) containing a total of 12 and 8 antibiotics/concentrations for respiratory and digestive pathogens, respectively. The antimicrobials tested for swine respiratory pathogens belong to category D [
15]: Sulfamethoxazole/trimethoprim, doxycycline, oxytetracycline and amoxicillin; Category C: Florfenicol, tiamulin, tulathromycin, tildipirosin and tilmicosin and category B: Ceftiofur, enrofloxacin and marbofloxacin. On the other hand, the antimicrobials tested for swine enteric pathogens belong to category D: Sulfamethoxazole/trimethoprim and spectinomycin; Category C: florfenicol, apramycin, gentamycin, neomycin and amoxicillin/clavulanic acid and category B: ceftiofur, cefquinome enrofloxacin, marbofloxacin and colistin.
Bacteria were thawed, cultured on chocolate agar or blood agar, and incubated at 35-37ºC in aerobiosis (or with 5-10% CO2 for APP) for 18-24h. Three to five colonies were picked and emulsified in demineralized water (or Cation Adjusted Muëller-Hinton Broth (CAMHB) for APP) to obtain a turbidity of 0.5 McFarland standard (Sensititre™ nephelometer V3011). Suspensions were further diluted in CAMHB for E. coli, CAMHB or CAMHB with 2.5-5% Lysed Horse Blood for P. multocida and Veterinary Fastidious Medium (VFM) or Mueller Hinton Fastidious broth with Yeast (MHF-Y) for APP to reach a final inoculum concentration of 5x105 cfu/ml. Then, the Sensititre panel was reconstituted by adding 100μl/well of the inoculum. Plates containing E. coli isolates were incubated at 35 ± 2ºC for 16-20h, P. multocida isolates were incubated at 35 ± 2ºC for 18-24h. In the case of APP isolates, plates were covered with a perforated seal and incubated at 35 ± 2ºC with 5-10% CO2 for 20-24h.
The antibiotic panels were read manually using Sensititre™ Vizion (V2021) and the MIC value was established as the lowest drug concentration inhibiting visible growth. For each isolate tested, a colony count and a purity check were performed following CLSI and manufacturer recommendations. Moreover, quality control strains were also included. Thus,
Actinobacillus pleuropneumoniae (ATCC 27090™),
Escherichia coli (ATCC 25922™),
Streptococcus pneumoniae (ATCC 49619™) and
Enterococcus faecalis (ATCC 29212™) were included as quality control following CLSI recommendations [
31,
32]. The MICs of the quality control strains had to be within acceptable CLSI ranges to accept the results obtained in the laboratory.
4.5. Statistical methods
All the data analysis was carried out with JMP®, Version 13 (SAS Institute Inc., Cary, NC, USA, 1989–2019). Descriptive statistics (MIC range, MIC
50 and MIC
90) were performed to summarize the distribution of the isolates within each MIC category. The number of categories was based on distributing the range of MIC values in four categories (from one to four) that include two MIC values for category, being category one the most susceptible (lowest MIC value) and category four the less susceptible (highest MIC value). The range of concentrations tested were 0,06-8, 0,03-4, 0,25-32 g/mL for 3rd and 4th cephalosporins, quinolones and polymyxins, respectively. Moreover, clinical susceptibility (susceptible/resistant for each isolate) was determined according to CLSI clinical breakpoints for APP,
P. multocida, and
E. coli for quinolones and cephalosporins and EUCAST guidelines for colistin in the case of E. coli, respectively [
43,
44,
45] (
Table 12).
A logistic (susceptible/resistant for each isolate) and multinomial logistic regression model (four MIC categories) was used to analyze the susceptibility data for the antimicrobials from year 2019 to 2022, only for those pairs of antimicrobial/microorganisms if at least 100 isolates were available for each year, as recommended by De Jong et al (2022) [
16]. Susceptible/resistant and categorized MIC data (MIC category 1, 2, 3 and 4) were used for logistic and multinominal logistic regression model, respectively as dependent variables, and the year as independent one. Thus, year of sampling was categorized by individual years and modelled as a hierarchical indicator variable, where for each year the preceding year was used as the referent [
52]. The final multinomial model was executed with outcome category 1 as the base referent category (the most susceptible one). The model assumptions and goodness-of-fit were evaluated as appropriate for these models [
52]. Thus, the level of significance used to reject the null hypothesis was p ≤ 0.05.