ESBL/AmpC-Producing Escherichia coli in Wild Boar: Epidemiology and Risk Factors

Simple Summary Antimicrobial resistance (AMR) represents a complex global issue due to the many factors involved. Extended-spectrum β-lactamase and AmpC (ESBL/AmpC)-producing Escherichia coli deserves attention for its broad repercussions on public health. Moreover, wild host species are of interest, particularly wild boar. Indeed, the constantly increasing population densities and the limited data on AMR in this species lead to health risks where spatial overlap with humans and domestic animals occurs. Therefore, 1504 wild boar fecal samples were analyzed to investigate ESBL/AmpC-producing E. coli and the effects of host-related factors and of human population density on their spread. A high prevalence of ESBL/AmpC-producing E. coli emerged in wild boar, species not treated with antibiotics, supporting that infection may be acquired through environmental contamination, whether of human or animal origin. Young animals were more colonized than older ones, demonstrating higher susceptibility as seen in domestic animals. Moreover, a positive association recorded between frequency of the TEM resistance gene and human population density suggests that spatial overlap may influence the infection in wild boar. Further analyses would be desirable to investigate the origin of the recorded environmental contamination, although a role of wild boar as a maintenance host of AMR strains emerged. Abstract The complex health problem of antimicrobial resistance (AMR) involves many host species, numerous bacteria and several routes of transmission. Extended-spectrum β-lactamase and AmpC (ESBL/AmpC)-producing Escherichia coli are among the most important strains. Moreover, wildlife hosts are of interest as they are likely antibiotics free and are assumed as environmental indicators of AMR contamination. Particularly, wild boar (Sus scrofa) deserves attention because of its increased population densities, with consequent health risks at the wildlife–domestic–human interface, and the limited data available on AMR. Here, 1504 wild boar fecal samples were microbiologically and molecularly analyzed to investigate ESBL/AmpC-producing E. coli and, through generalized linear models, the effects of host-related factors and of human population density on their spread. A prevalence of 15.96% of ESBL/AmpC-producing E. coli, supported by blaCTX-M (12.3%), blaTEM (6.98%), blaCMY (0.86%) and blaSHV (0.47%) gene detection, emerged. Young animals were more colonized by ESBL/AmpC strains than older subjects, as observed in domestic animals. Increased human population density leads to increased blaTEM prevalence in wild boar, suggesting that spatial overlap may favor this transmission. Our results show a high level of AMR contamination in the study area that should be further investigated. However, a role of wild boar as a maintenance host of AMR strains emerged.


Introduction
Antimicrobial resistance (AMR) represents a global health problem that involves humans, animal species and ecosystems [1][2][3]. Animals can play a role as reservoir of pathogenic and non-pathogenic antimicrobial-resistant microorganisms that, through direct contact or toward the food chain, can contribute to the spread and/or maintenance of AMR [4].
Escherichia coli is one of the most important bacteria that contributes to the complexity of AMR, and among its strains, those producing β-lactamases, extended-spectrum β-lactamases (ESBLs) and other β-lactamases such as AmpCs encoded by plasmid-located genes are of particular interest. Indeed, resistance to most beta-lactam antibiotics, including third and fourth-generation cephalosporins, can limit treatment options in case of infection [5] and lead to major public health concerns [6]. Authors reported 8750 deaths ascribed to E. coli resistant to third-generation cephalosporins in the European Union and European Economic Area during 2015, and comparing these data with others from 2007, an approximately 4-fold increase was detected [7]. The presence of E. coli ESBL and AmpC has been reported in several animal species [7][8][9], but, among them, wildlife is of particular epidemiological interest [10].
Wild species are generally not treated with antimicrobials and may acquire AMR bacteria or antimicrobial residues just through food and water in environments contaminated by domestic animals or humans [11]. Indeed, wildlife in anthropized environments or in close contact to humans or agricultural areas showed higher levels of AMR [12]. Thus, wildlife can be assumed as an important indicator for assessing the environmental spread of antimicrobial bacteria or AMR genes [4].
Among free-range species, wild boar (Sus scrofa) has been of high interest because of its role of an 'invasive species', the very close phylogenetical similarity to pigs (Sus scrofa domesticus) and the great increase of its population densities [13]. The continuous expansion of the habitat ranges of this species in close proximity to farms and pastures (anthropized environment) leads to opportunities for spatial overlap and contact with other species favoring health risks of sharing infections across the wildlife-domestic animal-human interface [14], such as antimicrobial resistant E. coli. Indeed, the use of antibiotics in both veterinary and human medicine and the known relationship between usage of antimicrobials and the occurrence of resistance in bacterial isolates in manure or manure contaminated surface water leads to AMR risks for wild boars in anthropized environments and highlights the role of drug resistance "conductor" of the environment [11,15]. E. coli is constantly present both in the digestive tract and in the environment, e.g., in water and soil [15], and this can contribute to the transfer of resistant genes between strains and to the increase in the drug resistance in environmental bacteria [15], an issue that is not yet sufficiently researched. Moreover, wild boar can act as a potential, but still to be defined, reservoir of AMR [16], maintaining strains/AMR genes through direct or indirect contact within its populations. This issue of the likely transmission of these pathogenic strains to humans through handling or consumption of contaminated game meat [7,17] should be taken into account. Nevertheless, data on ESBL and AmpC-producing E. coli in wild boar are still limited, especially with respect to factors that may influence its spread or the potential health risks related to the contact of these populations with humans or domestic animals. Therefore, we carried out an epidemiological investigation of ESBL/AmpC-producing E. coli in free-ranging wild boar from Northwest Italy and aimed to evaluate (i) the prevalence within the population and (ii) the effects of host-related factors and of human population density on the spread of these pathogenic bacteria.

Study Area
The study area was in Northern Italy within the easternmost part of the Lombardy region (Province of Brescia; 45 • 32 20 N, 10 • 13 10 E). The territory includes four wild boar hunting areas (Figure 1). From the official data on hunting activities provided by the local hunting office, consistent hunting behavior (dependent on the number of hunters and hunting days) was assumed among hunting areas and years. For this reasons, Chiari et al. [18] used the total number of wild boar hunted per year as an approximation of wild boar abundance and calculated a relative abundance index to take into account the different sizes of hunting areas, scaling abundance to district areas in km 2 ( Table 1). For each hunting area, the human population density was calculated based on data from the Italian Statistical Institute (ISTAT). Briefly, the number of people inhabiting the municipalities in the hunting areas was collected. When only a part of a municipality was included in a hunting area, the population was divided based on the percentage of the area in the municipality assuming that people were evenly distributed. To calculate the human population density within each hunting area, the total population present was scaled to district areas in km 2 ( Table 1).

Study Area
The study area was in Northern Italy within the easternmost part of the Lom region (Province of Brescia; 45°32′20″ N, 10°13′10″ E). The territory includes four wil hunting areas (Figure 1). From the official data on hunting activities provided by th hunting office, consistent hunting behavior (dependent on the number of hunte hunting days) was assumed among hunting areas and years. For this reasons, Chia [18] used the total number of wild boar hunted per year as an approximation of wil abundance and calculated a relative abundance index to take into account the di sizes of hunting areas, scaling abundance to district areas in km 2 ( Table 1). For each ing area, the human population density was calculated based on data from the Statistical Institute (ISTAT). Briefly, the number of people inhabiting the municipal the hunting areas was collected. When only a part of a municipality was include hunting area, the population was divided based on the percentage of the area in th nicipality assuming that people were evenly distributed. To calculate the human p tion density within each hunting area, the total population present was scaled to d areas in km 2 ( Table 1).

Isolation and Identification of ESBL/AmpC E. coli
The identification of β-lactamase-producing E. coli was performed through a double synergy diagnostic method. Specifically, 1 g of feces was diluted in 9 mL (1:10 dilution) of brain-heart infusion (BHI) broth supplemented with 1 mg/L cefotaxime for a preenrichment phase. After an overnight incubation, a drop of the BHI broth was used to inoculate MacConkey agar supplemented with 1 mg/L cefotaxime [19]. Positive growths were identified as pink to dark-pink colonies (lactose +), tested for activity of indole (+) and citrate (−), and one of these was selected for further molecular characterization. On the basis of staining and morphological characteristics, suspect colonies were selected and subjected to subsequent molecular characterization. A single bacterial colony from each phenotype-positive sample was resuspended in 250 µL of DNase-RNase free water, and DNA was extracted by lysis boiling (98 • C for 10 min). Identification of E. coli was conducted by PCR phylogenetic group analysis according to Clermont et al. [20]. This is a quadruplex method that can detect the eight E. coli phylogroups A, B1, B2, C, D, E, F and clade I and can exclude other strains with typical Escherichia phenotype, i.e., E. albertii and E. fergusonii.

Analysis of Resistance Genes
The detection of resistance genes was performed with a panel of reactions. Specifically, multiple PCR was used for the identification of bla CTX-M [21]; single PCRs were used for bla SHV [22], bla TEM [23] and bla CMY [24]. All the PCR gene targets, thermal profiles and primer sequences used for the detection of ESBL/AmpC E. coli are described in Table S1. In Figure S1 the agarose gel electrophoresis analysis of each specific genes, bla CTX-M , bla SHV , bla TEM and bla CMY , is shown. All the amplicons found to be positive for bla TEM and bla SHV were sequenced [25].

Statistical Analysis
The analysis can be split into three components: descriptive analysis of hunted wild boar, estimation of the apparent prevalence (p) with the corresponding confidence intervals estimated using an exact method and analysis of factors associated with prevalence. We used logistic regression models (GLM) [26] to define the effects of host-related explanatory variables (gender, age class, hunting season, hunting area, wild boar abundance) and environment-related explanatory variables (human population density) on the probability of wild boar to be positive for ESBL-/AmpC-producing E. coli. To investigate the effects of the same explanatory variables on AMR genes spread, we modeled each gene of microbiological positive subjects using a Binomial GLM. The models were compared using the likelihood ratio test (LRT). A post-hoc Tukey test was used to perform pairwise comparisons when statistically significant factors were detected in the models. Statistical analyses were carried out using R software version 4.0.3 [27], specifically the binom and emmeans packages. Results were considered statistically significant if the p-value (p) was less than 0.05.

Results
The phylogenetic analysis showed the presence of 240 E. coli distributed among all the eight phylo-groups.

Discussion
The present study revealed a high prevalence of ESBL/AmpC-producing E. coli in free-range wild boar. Young animals were more colonized by ESBL/AmpC strains, and particularly by the bla CTX-M resistance gene, than older subjects. Wild boar hunted during hunting season 2017-2018 showed a higher prevalence of ESBL/AmpC strains, bla CTX-M and bla TEM than animals of 2018-2019 and 2019-2020. The increasing of human population density leads to an increase of prevalence of bla TEM in wild boar.
In this study the most frequent ESBL/AmpC-producing E. coli phylogenetic groups were B1 (26.25%, 63/240) and A (24.58%, 59/240), in accordance with Holtmann et al. [7] (B1 = 54.55%, 12/22; A = 4/22, 18.18%). The results of both ESBL/AmpC-producing E. coli and AMR genes recorded are consistent with previous reports for this species [28][29][30]. Although our prevalence of ESBL/AmpC-producing E. coli (15.96%, 240/1504) appears higher than those previously recorded in wild boars by Wasyl et al. [1] (2.7%, 9/332), Literak et al. [31] (2%, 5/293), Bonardi et al. [32] (0.9%, 1/108) and Holtmann et al. [7] (5.9%, 22/375), the lack of a standardized diagnostic method makes difficult the comparison of results of different studies [12]. However, bla CTX-M (12.3%, 185/1504) was confirmed as the most frequent among AMR genes as reported by Literak et al. [31] (80%, 4/5) and Holtmann et al. [7] (68.18%, 15/22). The findings of TEM-1 (92/98) and SHV-12 (7/7) as the most frequent variants of bla TEM and bla SHV are consistent with the results of Holtmann et al. [7] (TEM-1 = 5/7; SHV-12 = 4/4). Considering that free-range species should be free of antibiotics/AMR, the recorded prevalence suggests that the acquisition of ESBL/AmpCproducing E. coli or of AMR genes occurs through direct contact with the environment, soils or surface water (feeding and drinking) contaminated with livestock/pig manure [7,33]. A possible explanation for the differences in prevalence recorded between our study and the previous ones could be related to a higher AMR environmental contamination of our study area. Indeed, all samples were collected in an area with a high density of farms, and Italy is one of the largest consumers of veterinary antimicrobials in Europe [34], even in age groups such as finisher pigs, where the use of antimicrobials should be low [35]. In this regard, the intrinsic characteristics of the wild boar's diet, an opportunistic omnivore that can change its diet drastically depending on the availability of food, should also be taken into account. Although currently no data on the prevalence of AMR in domestic animals living in spatial overlap with wild boar are available to be included in the analysis, further investigations will address this aspect.
Statistical models showed a higher prevalence of E. coli ESBL/AmpC (20.58%) and bla CTX-M (16.52%) in young animals compared to sub-adults (ESBL/AmpC = 13.32%, bla CTX-M = 10.05%) and adults (ESBL/AmpC = 15.17%, bla CTX-M = 11.5%). To the best of our knowledge, in the only available reference in which an age comparison of prevalence of these strains was performed for wild boar, no significant differences were registered [7]. However, the fact that in pig farms piglets have an increased susceptibility to E coli strains, which can cause severe enteritis especially when they are only a few weeks old [36], suggests that young wild boar may also have an increased susceptibility to this infection. Furthermore, similarly to domestic animals, young individuals may have a greater susceptibility to these strains due to a not yet fully developed intestinal microflora and may be more susceptible to new bacterial colonization [37]. In order to better understand the transmission of E. coli that occurred in young wild boars, an in-depth analyses of virulence factors, i.e., adhesins (K88, K99, F41, 987P and F18), of isolated E. coli would be useful to determine whether these strains were part of the normal microflora or could potentially cause diseases in domestic pigs.
The , 2019/2020 = 5.85%) recorded between hunting seasons could be associated with both the role of the environmental indicator of wild boar and that of the "conductor" in the transmission of pathogenic strains/resistance genes of the environment [2]. The suggested hypothesis is that until 2017, the natural ecosystem of wild boar was more contaminated by resistant bacteria than in the subsequent years. However, with our data we cannot discriminate which mechanism actually occurs and what origin this contamination had, whether human or animal, in light of the surprising rapidity with which this change in trend of prevalence between hunting seasons occurred. However, the occurrence of bla TEM in E. coli of wild boar being positively related to human population density suggests that this spatial overlap, particularly due to feeding with anthropogenic food sources and waste, can expose wild boar to traces of antimicrobial substances or directly to resistant bacteria of human origin [7]. In the late 1990s, when most ESBLs were mutants of the classical TEM-1, -2 and SHV-1 enzymes [38], bla TEM was the most frequent isolated among nosocomial strains [39], being also detected in commensal strains from healthy adults [40]. From the 2000s the bla CTX-M genes became the dominant enzymes in human populations [41]. Indeed, Livermore et al. [38] reported that in Italy, ESBL prevalence increased steadily from 0.2% in 1999 to 1.6% in 2003, with this change largely attributable to CTX-M-β-lactamase-positive isolates, which comprised 38% of all ESBL producers in 2003 compared with 12% in 1999. Our results may appear surprising considering the findings of Holtmann et al. [7] who showed a positive relationship between ESBL strains and human population density. Nevertheless, as bla TEM still occurs in the most clinically relevant enterobacterial species [39] with the emergence of its new variants [42], and given that no data are currently available on the prevalence of bla CTX-M in the human population of our study area, the hypothesis that bla CTX-M may not yet be as widespread in the human population and therefore does not have a direct effect on wild boar could be suggested.
The high prevalence of ESBL/AmpC-producing E. coli recorded in wild boar in this study leads to strategies for reducing and preventing the spread of resistant bacteria in the environment and between species. As the main concern of AMR is related to the overuse of antimicrobials [43], the reduction of antibiotic administration in domestic animals should be the starting point. In addition, farm biosecurity standards and infrastructure should be improved in order to minimize the impact of animal husbandry manure in the environment and to prevent the related potential spreading of AMR bacteria. Attention should also be paid to preventing wild boar feeding with anthropogenic waste or food sources that can further expose them to traces of antimicrobials or directly to resistant bacteria of human origin.

Conclusions
The high prevalence of ESBL/AmpC-producing E. coli recorded in this study shows the potential role of wild boar as a maintenance host of AMR strains. However, the creation of historical data series would be desirable in order to monitor and better evaluate these temporal trends in wild boar. In addition, the role of this species as indicator of environmental contamination, supported by the recorded positive association between bla TEM and human population density, leads to an in-depth analysis that also includes domestic animals. Indeed, through a "One Health" approach, further specific sampling should be carried out at both farm and human/environmental levels within hunting areas to assess prevalence of these strains/AMR genes and to evaluate whether their presence actually influences their spread in wild boar.