- Research
- Open access
- Published:
Survey in ruminants from Rwanda revealed high diversity and prevalence of extended-spectrum cephalosporin-resistant Enterobacterales
BMC Veterinary Research volume 20, Article number: 523 (2024)
Abstract
Background
Antimicrobial resistance (AMR) in Enterobacterales constitutes a significant threat to the health of both humans and animals and a socioeconomic problem. Enterobacterales, mainly Escherichia coli, carrying β-lactamases has become one of the main indicators to estimate the burden of AMR in animals within “One Health” approach.
Objectives
To assess the presence of extended-spectrum cephalosporin-resistant Enterobacterales associated with ruminants (cattle, sheep, goats) habituated in all five provinces of Rwanda and to perform in depth characterization of isolates.
Methods
We screened 454 rectal swabs from 203 cows, 170 goats, and 81 sheep and selective isolation of extended-spectrum cephalosporin-resistant Enterobacterales was conducted. Isolates were identified as a members of the order Enterobacterales by MALDI-TOF MS and further characterized by susceptibility testing and by whole-genome sequencing.
Results
Out of the 454 samples, 64 extended-spectrum cephalosporin-resistant Enterobacterales were isolated from 58 animals. Isolates belonged to seven bacterial species and were identified as Escherichia coli (n = 54), Enterobacter bugandensis (n = 4), Enterobacter mori (n = 2), Klebsiella pneumoniae (n = 2), Enterobacter dykesii (n = 1), and Citrobacter freundii (n = 1). All isolates displayed an Extended-spectrum β-lactamases (ESBL) phenotype, with exception of Citrobacter freundii isolate displayed both an ESBL and AmpC phenotype. In addition, all Enterobacter isolates were identified as stably de-repressed AmpC-producers. ESBLs genes, blaCTX−M−15 was predominant. Resistance to tetracycline and tet(A) was most frequently observed among non-β-lactam resistance. Forty-eight isolates displayed multidrug-resistance phenotypes. A shiga toxin-producing E. coli and an enterotoxigenic E. coli isolate were observed. Genome comparisons revealed thirty-five E. coli sequence types (ST) (ST10, ST307 being predominate).
Conclusions
Considering the high proximity between ruminants and humans in Rwanda, the dissemination of antimicrobial drug resistance highlights the public health threats and requires the joint and multisectoral action of human and veterinary medicine, at human-animal-environment interfaces. Therefore, it is important to establish national and global “One Health” surveillance programs of AMR to tackle the antibiotic-resistant crisis in human and veterinary medicine.
Background
The different β-lactam antibiotics are most extensively used for the treatment of bacterial infections in human and veterinary medicine [1, 2]. Antimicrobial resistance (AMR) in bacteria constitutes a significant threat to the health of both humans and animals and a socioeconomic problem [3]. For multidrug-resistant members of Enterobacterales – and primarily β-lactamase producers –limited therapeutic options exist [4]. The most common β-lactamases are SHV, TEM and CTX-M variants of extended-spectrum β-lactamases (ESBLs) [5, 6]. ESBLs confer resistance to penicillins, I., II., and III. generation cephalosporins and monobactams [7], but can inhibited by β-lactamase inhibitors [5]. ESBL-producing Enterobacterales are associated with a high burden in terms of morbidity and mortality [8]. Therefore, the World Health Organization (WHO) listed ESBL-producing Enterobacterales among others in the most critical group according to the urgency of the need for the development of new antibiotics [9]. Enterobacterales harboring ESBL genes have been isolated across continents from humans, domestic and wild animals, and the environment [3, 7, 10].
Rwanda recognized the challenging situation of AMR as a significant public health problem and a threat to its unique population of mountain gorillas (Gorilla beringei ssp. beringei) and other endangered wildlife. In consequence, based on five objectives of the Global Action Plan of the WHO [11], Rwanda developed the National Action Plan on Antimicrobial Resistance [12]. Nevertheless, information on AMR in bacteria associated with human infections and domestic animals remains scarce.
The present study aimed to determine the presence of extended-spectrum cephalosporin-resistant Enterobacterales associated with ruminants (cattle, sheep, goats) located in all five provinces of Rwanda and to characterize these isolates.
The Rwandan economy relies mainly on the agriculture sector, where subsistence farming is still predominant. Cattle farming has a long tradition and is anchored in the history of the country [13]. Farmers have free access (without prescription) to veterinary drugs, like antibiotics and most medicines, while on the other hand, microbiological diagnostic options in veterinary medicine are limited [12].
Materials and methods
Sample collection and identification of Enterobacterales
To determine the presence of enteric colonization by extended-spectrum cephalosporin-resistant Enterobacterales, a total of 454 rectal swabs from 203 cows, 170 goats, and 81 sheep, were sampled between September until December 2021. Animals were selected based on owner´s willingness to participate to the study and only if they did not receive any antimicrobial treatment in the past six months.
The study was discussed, and the rectal swabbing was approved by the Research Screening and Ethical Clearance Committee of the College of Agriculture, Animal Sciences and Veterinary Medicine, University of Rwanda (003/2021/DRI).
In the microbiological laboratory at the New Vision Veterinary Hospital in Rwanda (https://nvvh.rw/) each rectal sample was incubated at 37 °C overnight in buffered peptone water (Merck), with cefotaxime (1 mg/L) and then cultivated at 37 °C overnight on MacConkey agar (Rapid Labs, Colchester, United Kingdom) supplemented with cefotaxime (1 mg/L) (MacCTX). After incubation on MacCTX, one colony representing a distinct colony morphotype (e.g. shape, size, color, mucoid yes/no, lactose fermentation was subcultured on the same medium and cryo-conserved. For species identification and further in-depth characterization, isolates were sent to the Institute of Microbiology, University of Veterinary Medicine Vienna, Austria. Isolates were identified to the species level by matrix-assisted laser desorption ionization-time-of-flight mass spectrometry (MALDI-TOF MS) (Bruker Daltonik, Bremen, Germany). Only isolates that belonged to the order Enterobacterales and that could grow on MacCTX were selected for further characterization.
Susceptibility testing
Antimicrobial susceptibility testing was performed by agar disk diffusion according to the Clinical & Laboratory Standards Institute (CLSI) standards [14]. Escherichia coli ATCC® 25922 served as the quality control. Disks containing the following antimicrobial agents were used: cefotaxime (30 µg), ceftazidime (30 µg), cefoxitin (30 µg), meropenem (10 µg), gentamicin (10 µg), tobramycin (10 µg), amikacin (30 µg), ciprofloxacin (5 µg), trimethoprim-sulfamethoxazole (1.25/23.75 µg), tetracycline (30 µg), chloramphenicol (30 µg), fosfomycin (200 µg), and nitrofurantoin (300 µg) (Becton Dickinson, Heidelberg, Germany). AmpC-hyperproducing isolates (stably de-repressed) of Enterobacter spp. and Citrobacter spp. were defined as previously described [15].
Molecular analyses
Isolates were analyzed by whole-genome sequencing (WGS). Isolation of bacterial DNA, library preparation, and paired-end-sequencing using NextSeq2000 platform were performed as described previously [16]. De novo assembly of the raw reads was performed using SPAdes v.3.9.0 [17], and typing was performed with SeqSphere + software (Ridom, Münster, Germany). The species allocation was confirmed by JSpecies workspace using the ANIb (average nucleotide identity via BLAST) analysis tool [18].
E. coli phylotypes were determined based on WGS by the Clermont Typing tool (http://clermontyping.iame-research.center/, accessed on 2023 May 30) [19, 20]. Clonal relatedness of E. coli isolates was assessed by two-locus sequence typing of combined data of fumC and fimH sequences using CHTyper [21] hosted at the Center for Genomic Epidemiology (https://cge.cbs.dtu.dk/services/chtyper/, accessed on 2023 May 30). Sequence types (STs) of E. coli were determined using the Escherichia/Shigella database hosted on EnteroBase [22]. SeqSphere + software (Ridom, Münster, Germany) was used for core-genome Multilocus sequence typing (cgMLST) analysis of E. coli and Klebsiella pneumoniae isolates. The same software was used to predict serogenotypes of E. coli. STs of K. pneumoniae were determined using the K. pneumoniae species complex database hosted at the Institute Pasteur (https://bigsdb.pasteur.fr/klebsiella/). The same database was used to examine the presence of virulence-associated genes (VAGs). STs of Enterobacter (En.) cloacae complex and Citrobacter (C.) freundii were assigned using databases hosted at PubMLST (Enterobacter cloacae complex: https://pubmlst.org/organisms/enterobacter-cloacae, Citrobacter sp.: https://pubmlst.org/organisms/citrobacter-spp [23]. ABRicate [24] v1.0.0 was used to identify antimicrobial resistance genes using ResFinder 4.1 [25], VAGs were identified using Virulence Factor Database [26], and virulence-associated genes in E. coli using VirulenceFinder 2.0 [27]. The presence of plasmids was assessed using PlasmidFinder 2.1 [28]. The beta-lactamase database (BLDB) was used to classify some of β-lactamase genes [29]. Probability prediction of the location of a given bla resistance gene in Enterobacterales was achieved by applying mlplasmids trained on E. coli or Klebsiella [30]. Posterior probability scores (ppp) > 0.7 and a minimum contig length of 700 bp indicate that a given contig is plasmid-derived. For bla genes with ppp scores below 0.7, Basic Local Alignment Search Tool (BLAST) searches with the respective contig sequence were performed against all Enterobacterales plasmids at National Center for Biotechnology Information (NCBI). The most promising plasmid entries were then used as backbones in mapping all plasmid contigs of a given isolat using GENEious Prime software 2022.2.2. A plasmid location was assumed only if the combined BLAST/GENEious analysis was able to identify a published plasmid which could be reconstructed with > 90% identity. Additional gene prediction and annotation were performed using the Bioinformatics Resource Center PATRIC [31].
The genomes of isolates subjected to WGS were deposited under PRJNA950594 in the NCBI BioProject database.
Statistical analyses and data visualization
Statistical analyses and visualizations were computed in R 4.2.3 [32]. Descriptive statistics were used to characterize the study population and describe the results. The 95% confidence interval (95%CI) of the prevalence was calculated using the R package epiR [33] A generalized linear model (GLM) with binomial distribution and logit function was applied to identify variables that may explain the presence/absence of extended-spectrum cephalosporin-resistant Enterobacterales in animals (presence was defined as the carriage of at least one resistant isolate). The impact of the following explanatory variables was assessed: province (Northern, Eastern, Western, Southern, Kigali City), animal species (bovine, ovine, caprine,, with ovine used as reference category), and sex (female, male, with female used as reference category). The best-fitted model (outlining the most important predictors and their effects on the outcome variable) was identified using a stepwise backward selection based on the Akaike Information Criterion. Estimates were transformed to odds ratios (OR). Moreover, difference between prevalences was analyzed by using the 2-tailed 2-proportion z-test.
Results
Isolation of cefotaxime resistant isolates, phenotypic resistance and species identification
A large majority of the sampled animals were female (95.4%). Animals´ age ranged between 0.3 and 11 months (29% of the data on age was unavailable). Overall, 166 swabs were collected in the Northern Province (71 bovines/51 caprines/44 ovines), 127 in the Eastern Province (66/57/4), 103 in the Western Province (38/49/16), 55 in the Southern Province (25/13/17), and three in Kigali City (3/0/0). Out of the 454 samples, 64 extended-spectrum cephalosporin-resistant Enterobacterales were isolated from 58 animals (prevalence = 12.8%; 95%CI = [9.8–16.2]), including 33 bovines (16.3%; 95%CI = [11.5–22.1]), 20 caprines (11.8%; 95%CI = [7.3–17.6]), and 5 ovines (6.2%; 95%CI = [2.0–13.8]). The prevalence in females was 12.9% (95%CI = [9.9–16.5]) and was statistically not significant (p = 0.9) from 9.5% in males (95%CI = [1.2–30.4]).
Isolates belonged to seven bacterial species and were identified as E. coli (n = 54), En. bugandensis (n = 4), En. mori (n = 2), K. pneumoniae (n = 2), En. dykesii (n = 1), and C. freundii (n = 1) (Additional file 1: Table S1). All E. coli and K. pneumoniae isolates displayed an ESBL phenotype, whereas the C. freundii isolate displayed both an ESBL and AmpC phenotype. In addition, all Enterobacter isolates were identified as stably de-repressed AmpC-producers.
Molecular characterization
Among genes associated with resistance to β-lactams, blaCTX−M−15 was identified in 47 isolates, followed by blaTEM−1 (n = 35). These two β-lactamase genes were observed together in 33 isolates. Among other blaCTX−M alleles, blaCTX−M−14 (n = 3), blaCTX−M−137 (n = 2), and single blaCTX−M−121 and blaCTX−M−27 were identified in E. coli isolates. The blaSHV genes were less frequently detected, blaSHV−208 genes were identified together with blaCTX−M−15 and blaTEM−1 in two K. pneumoniae isolates. The ESBL gene blaSHV−12 was detected in an E. coli isolate. The oxacillinase gene blaOXA−1 was observed in three E. coli isolates. The genes belonging to blaEC family belonging to class C β-lactamases (cephalosporinase) were observed in 53 isolates. In addition to one not further specified blaEC gene (designated blaEC−NT) and the following blaEC genes were identified blaEC−1 (n = 3), blaEC−4 (n = 2), blaEC−8 (n = 10), blaEC−13 (n = 8), blaEC−15 (n = 9), blaEC−18 (n = 15), blaEC−1502 (n = 2), blaEC−2021 (n = 1), and blaEC−2045 (n = 2) (Additional file 1: Table S1).
Among non-β-lactam resistance, resistance to tetracycline (n = 45; 70.3%), trimethoprim-sulfamethoxazole (n = 41; 64.1%), fosfomycin (n = 13; 20.3%), chloramphenicol (n = 11, 17. 2%), tobramycin (n = 8, 12.5%), ciprofloxacin (n = 6, 9.4%), gentamicin (n = 3, 4.7%), and nitrofurantoin (n = 1, 1.6%) was detected. Visualization of genotypes per isolated bacterial species is shown in Fig. 1 The most prevalent non-β-lactamase genes detected were tet(A) (n = 41), sul2 (n = 38), strB (also known as aph [6]-Id) (n = 38), aph(3’’)-Ib (n = 33), dfrA14 (n = 29), and qnrS1 (n = 28). Forty-eight isolates displayed multidrug-resistance phenotypes [34].
Among the ciprofloxacin-resistant isolates, significant mutations in the quinolone resistance-determining regions (QRDR) of gyrA, parC and parE were observed. The mutations led to amino acid substitutions at positions 83 (Ser→Leu) and 87 (Asp→Asn) in GyrA in five isolates and to one exchange at position 83 (Ser→Leu) in one isolate. ParC was altered at position 80 (Ser→Ile) in all six ciprofloxacin-resistant isolates as well as at position 458 (Ser→Ala) in two isolates, at position 458 (Ser→Thr) in one isolate, and at position 416 (Leu→Phe) in another isolate (Additional file 1: Table S1).
Genes associated with biocide resistance were detected in 23 isolates, namely qacE (n = 13), sitABCD (n = 6), oqxA (n = 2), and oqxB (n = 2), which are associated with the decreased susceptibility to benzalkonium chloride, chlorhexidine, and hydrogen peroxide (Additional file 1: Table S1).
Different virulence-associated genes were observed in E. coli isolates, with the tellurite resistance gene terC (n = 54) and the fimH gene encoding the mannose-binding type 1 pilus tip protein FimH (n = 49) being predominant. Among E. coli pathotypes, one caprine Shiga toxin-producing E. coli (STEC) carrying stx1c and stx2b, one bovine enterotoxigenic E. coli (ETEC) isolate carrying astA (coding for heat-stable enterotoxin EAST-1) and eltIAB (coding for heat-labile enterotoxin LTIh) were observed. Among K. pneumoniae virulence factors, type 3 fimbriae encoded by the mrk operon and genes coding for yersiniabactin (ybt) were detected in all K. pneumoniae isolates (Additional file 1: Table S1).
E. coli phylotyping revealed that the most predominant phylogroups were A and B1 (each n = 20), followed by D (n = 9), C (n = 2), B1 mash E (n = 2) and cladeI (n = 1). The fumC and fimH (CH) typing divided the E. coli isolates into 38 distinct CH clonotypes, wherefrom CH clonotype 29–87 was predominant (n = 4). Two novel fimH alleles were detected: fimH3351 (isolate 511398-22) and fimH3352 (isolate 511421-22).
Overall, 44 STs were evidenced that showed a heterogeneous spatial distribution over the Rwandan provinces. Thirty-five E. coli sequence types (ST) (based on seven loci MLST extracted from WGS data) were detected: ST10, ST3076 (each n = 4), ST9849, ST3489 (each n = 3), ST69, ST5614, ST38, ST224, ST5909, ST196, ST2614, ST14131, ST349 (each n = 2) and single isolates of ST2334, ST14130, ST90, ST297, ST23, ST155, ST7809, ST1788, ST2465, ST13718, ST5708, ST361, ST6636, ST2178, ST75, ST398, ST847, ST448, ST9347, ST48, ST1303, and ST5204, wherefrom ST13718, ST14130, and ST14131 are novel STs. The analysis of WGS data revealed that ten isolates could be assigned to WGS-predicted serotypes (somatic O and flagellar H antigens) O1:H9, O5:H11, O8:H9 (n = 3), O8:H53, O15:H18, O101:H10 (n = 2), and O131:H45. In the remaining isolates, the O or H type could not be assigned (Additional file 1: Table S1). WGS-based cgMLST revealed eight clusters comprising four isolates (Cluster 1, ST3076), three isolates (Cluster 2, ST9849), two isolates (Cluster 3, ST14131; Cluster 4, ST38; Cluster 5, ST196; Cluster 6, ST224; Cluster 7, ST5614, and Cluster 8, ST3489) (Fig. 2). All other isolates were singletons. Out of seven Enterobacter spp. isolates, six were novel, ST2013 (En. dykesii), ST2014, ST2015, ST2017 (En. bugadensis), ST2016 and ST2088 (En. mori). ST733 (En. bugadensis) was an already established ST. Other STs obtained were ST610 (both K. pneumoniae isolates) and C. freundii ST163 (Additional file 1: Table S1).
Minimum spanning tree comparing E. coli isolates. Easch circle represents isolates with an allelic profile based on the cgMLST which consists out of 2513 alleles. Numbers correspond to the allelic differences between isolates; isolates with closely related genotypes are in bolded circles and marked as clusters. Isolates were labeled according to isolate ID (/) classical MLST
The presence of the plasmid replicons among β-lactamase gene carrying isolates was analyzed by PlasmidFinder and revealed that ten different plasmid replicons could be detected, with IncB/O/K/Z being predominant. The mlplasmids analyses and BLAST searches predicted that the majority of the β-lactamase genes might be plasmid-encoded and the majority of the other resistance genes appeared to be chromosomally encoded. In the present study, it was observed that some bla genes were carried on plasmids, including blaCTX−M−15 by IncB/O/K/Z, IncFIB(AP001918), IncFIB (H89-PhagePlasmid), IncFII, IncFII(K), IncHI2A, and IncQ1IncY plasmids, blaCTX−M−27 by IncFII (pRSB107) plasmid, blaTEM−1B by IncFII (pHN7A8), IncI1-I (Alpha), and IncY plasmids, blaSHV−12 on IncY plasmid. The full list of AMR genes and their predicted plasmid probability are shown in Additional file 2: Table S2. Genome data are shown in additional file 3: Table S3.
Occurrence of extended-spectrum cephalosporin-resistant Enterobacterales among different animals
The best GLM model selected according to the AIC values included only animal species as explanatory variable. The odds of extended-spectrum cephalosporin-resistant Enterobacterales occurrence was nearly threefold greater for bovines compared to ovines (OR = 2.95, 95% IC = [1.20–8.88]). No significant difference was observed between ovines and caprines (OR = 2.03, 95% CI = [0.79–6.28]).
Discussion
The present study aimed at investigating phenotypic and genotypic characteristics of extended-spectrum cephalosporin-resistant Enterobacterales isolates obtained from rectal swabs of ruminants (cattle, sheep, goats) from all five provinces in Rwanda. To the best of authors knowledge, the present study is the first study dealing with in-depth analyses of extended-spectrum cephalosporin resistance in Enterobacterales associated with ruminants in Rwanda.
In this study, the extended-spectrum cephalosporin resistance in Enterobacterales was detected at a frequency of 12.8% (58/454) in rectal samples of Rwandan ruminant livestock species. This is of interest because third- or fourth-generation cephalosporins are either not used for the treatment of domestic animals in Rwanda yet [35] and rarely used for the treatment of wild animals [36]. In addition, none of the animals examined during the present study was associated with known risk factors for the infection and/or colonization with AMR, such as recent veterinary health care (during the last six months) or treatment with antimicrobial substances in the preceding six months.
Antimicrobials are usually not used to treat goats and sheep in Rwanda. Even though the occurrence of extended-spectrum cephalosporin-resistant Enterobacterales in caprine and ovine samples was less frequent than in bovine, this is of particular concern. In Rwanda, small ruminants usually live in close proximity with their owners, often sharing the same premises with humans during the night. Interestingly, caprine ESBL-producing E. coli isolates were obtained solely from the Western provinces. Five out of 13 STs observed among caprine E. coli isolates have already been isolated from humans in Rwanda (ST10, ST155, ST361, ST38, and ST69) and from an animal of unknown species (ST38) [37]. Among them, members of ST10, ST38, ST69, ST155 belong to important pandemic clonal lineages, which are often multidrug-resistant and usually associated with numerous ESBL genes [38, 39]. With the exception of one E. coli harboring blaCTX−M−27 belonging to the ST90, and an E. coli isolate harboring blaSHV−12, all other caprine isolates carried blaCTX−M−15. CTX-M-15, a subtype in the CTX-M-1 group, is a globally disseminated β-lactamase and was detected in E. coli isolates associated with humans and animals [40]. The blaCTX−M−15 gene was also the most predominant β-lactamase gene among ESBL E. coli isolates in Rwanda from hospitalized patients, caregivers, and the community [37] as well as from flies from a tertiary hospital [41]. CTX-M-27, together with CTX-M-15 and CTX-M-14, are dominant among Enterobacterales worldwide [42]. While blaCTX−M−15 and blaCTX−M−27 positive E. coli have already been reported [37], to our knowledge, this is the first report of E. coli harboring blaSHV−12, in Rwanda. SHV-12 is an important ESBL associated with nosocomial and community-acquired infections in several countries and has also been reported in E. coli associated with domestic and wild animals [7, 43, 44]. Besides resistance-associated genes, genes coding for virulence factors were observed in all isolates. In the blaCTX−M−15-carrying isolate belonging to ST75, the Shiga toxin genes were observed. ST75 has already been reported in numerous countries in Asia, Europe, North America, and Oceania. Still, there are no ST75 entries from Africa neither in the Enterobase Escherichia/Shigella Database (https://enterobase.warwick.ac.uk/species/index/ecoli, last accessed 2023 May 30) nor in the E. coli PubMLST Database (https://pubmlst.org/organisms/escherichia-spp), last accessed 2023 May 30). In a multi-resistant blaCTX−M−15-E. coli ST297, genes coding for F17 fimbriae, the heat-labile enterotoxin LTII (eltIIAB), and other VAGs were observed. The fimbrial adhesin F17 is associated with ETEC in calves [45].
Among the characterized ovine isolates, an atypical enteropathogenic E. coli (aEPEC) belonging to ST2178, a common clone in Africa (https://enterobase.warwick.ac.uk/species/index/ecoli, last accessed 2023 May 30) was observed. This isolate carried the astA gene (coding for heat-stable enterotoxin EAST-1) and other virulence-associated genes (VAGs). EAST-1 has been associated with diarrheal disease in humans and numerous animal species, mainly cattle and swine [46].
Since cattle play such an important role in the lives of the people in Rwanda, the detection of multidrug-resistant E. coli isolates is of great importance, especially the very close proximity between humans and cattle induces a risk of animal-to-human transmission of AMR. Six out of 22 different STs detected among bovine ESBL E. coli isolates (ST10, ST196, ST224, ST48, ST448, and ST69) have previously been observed among human isolates [37], whereas ST10 and ST69 are well-known pandemic lineages [38, 39]. The blaCTX−M−15 was the most predominant ESBL gene among bovine isolates. Furthermore, blaCTX−M−14, blaCTX−M−121, and blaCTX−M−137 were observed. CTX-M-137, a hybrid of CTX-M-14 [47] and CTX-M-121 [48], both CTX-M-9 group-like β-lactamases, is rarely reported among E. coli isolates, and to our knowledge, there is no published data reporting the presence of E. coli harboring this β-lactamase gene in Rwanda. An astA gene accompanied by other VAGs was observed in three STs, ST10, ST349, and ST847. Even though EAST-1 was associated with diseased humans and animals, the virulence capacity and pathogenesis of the EAST-1 toxin remains uncertain and disputed [46].
The predominant E. coli phylotypes in the present study were groups A and B1. Both phylogroups are described as being widespread among human commensal E. coli isolates [49]. Phylogroup D was common in our study (n = 9). It was previously described to be associated with E. coli isolates associated with extra-intestinal infections in humans [50].
Based on cgMLST analysis eight clusters and 35 singletons were observed. Two clusters (cluster 3 (ST14131) and cluster 6 (ST196)) comprise two isolates from the same farm and originated from the same species (bovine or caprine). All other clusters contain isolates that originated (a) from the same sector and the same animal species but different farms (cluster 7 (ST38)), (b) from the same province and the same animal species but different sectors (cluster 1 (ST3076)), (c) from different provinces but the same species (clusters 5 (ST224) and 8 (ST3489)), and finally (d) from different provinces and different animals species (cluster 4 (ST5614). Interestingly, from one rectal sample two E. coli belonging to different clusters could be detected (strain IDs 511405-22 and 511415-22). These findings suggested clonal distribution and interspecies transmission of extended-spectrum cephalosporin-resistant E. coli in examined animals throughout Rwanda.
Two CTX-M-15 positive K. pneumoniae isolates have all characteristics in common (i.e. resistance pheno- and genotype, biocide resistance genes, identical set of virulence associated genes, sequence type) including the same cgMLST type. They originated from different provinces and animals. K. pneumoniae is a known nosocomial pathogen with increasing multi-drug resistant rates and global dissemination [51]. Reports about K. pneumoniae ST610 are scarce. ST610 has two entries in the K. pneumoniae species complex database hosted at the Institute Pasteur (https://bigsdb.pasteur.fr/klebsiella/, last accessed 2023 May 30), whose corresponding isolates originated from a bovine mastitis in the USA and a human in Ireland. K. pneumoniae harboring blaCTX−M−15 has been observed in neighboring countries [52,53,54] but to our knowledge, there are no reports from Rwanda.
The blaCTX−M−15 gene was also carried by an ST163 C. freundii isolate, which is not a common clone. Until present, ST163 has three entries in PubMLST, from Canada (environment) and two human patients from China and Israel, respectively. Recently in Rwanda multi-drug resistant C. freundii harboring blaCTX−M−15, blaOXA−10, blaOXA−1, blaTEM−1b, and several non-β-lactamase resistance genes were described [41]. These isolates belonged to the ST328, a rarely observed lineage among Citrobacter spp [55, 56].
The majority of β-lactamase genes detected during the present study were predicted to be plasmid-borne. Plasmids play a significant role in the dissemination of ESBL genes [5]. Therefore, these findings could be of particular interest for public health in Rwanda. Several recent studies reported the occurrence of chromosomally located β-lactamase genes. Therefore, the differentiation between plasmid and chromosome localization of β-lactamase genes could be challenging, which is a limitation of the method.
Conclusion
The present study contributes to the growing evidence that the presence of extended-spectrum cephalosporin-resistant Enterobacterales are currently part of the microbiome of ruminants in Rwanda. β-Lactamase genes frequently encoded on plasmids can easily be transferred horizontally. In addition, our data suggested that the spread of resistant clones could already have happened vertically through Rwanda. Considering the close contact between ruminants and humans in Rwanda, the dissemination of pandemic clones requires the joint and multisectoral action of human and veterinary medicine, at human-animal-environment interfaces. It is nowadays critical to establish national and global “One Health” surveillance programs of AMR to tackle the antibiotic-resistant crisis in human and veterinary medicine to ultimately protect the unique Rwandan wildlife, especially the population of endangered mountain gorillas (Gorilla beringei ssp. beringei) that plays a critical role in the lives of the people in Rwanda and in the resilience of the Rwandan natural ecosystem. Main limitations of this study rely in the challenges related to field work and data collection in resource-limited settings, difficult weather conditions, accessibility to farms, and communication with animal owners. Despite these inherent challenges, our results highlight the importance of AMR research in Africa.
Data availability
All data generated or analysed during this study are included in this published article [and its supplementary information files].
Abbreviations
- AMR:
-
Antimicrobial resistance
- aEPEC:
-
Atypical enteropathogenic E. coli
- BLAST:
-
Basic Local Alignment Search Tool
- BLDB:
-
Beta-lactamase database
- cgMLST:
-
Core-genome Multilocus sequence typing
- CH:
-
fumC and fimH
- ESBL:
-
Extended-spectrum β-lactamases
- ETEC:
-
Enterotoxigenic E. coli
- GLM:
-
Generalized linear model
- MLST:
-
Multilocus sequence typing
- NCBI:
-
National Center for Biotechnology Information
- QRDR:
-
Quinolone resistance-determining regions
- ST:
-
Sequence type
- STEC:
-
Shiga toxin-producing E. coli
- VAG:
-
Virulence-associated genes
- WGS:
-
Whole-genome sequencing
- WHO:
-
World Health Organization
References
Bush K, Bradford PA. β-lactams and β-lactamase inhibitors: An overview. Cold Spring Harb Perspect Med. 2016;6(8).
Ogunrinu OJ, Norman KN, Vinasco J, Levent G, Lawhon SD, Fajt VR, et al. Can the use of older-generation beta-lactam antibiotics in livestock production over-select for beta-lactamases of greatest consequence for human medicine? An in vitro experimental model. PLoS ONE. 2020;15(11):e0242195.
Murray CJ, Ikuta KS, Sharara F, Swetschinski L, Robles Aguilar G, Gray A et al. Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. Lancet. 2022;399(10325).
World Health Organization. Global antimicrobial resistance and use surveillance system (GLASS) report 2021. World Health Organization; 2021.
Bush K. Past and present perspectives on β-lactamases. Volume 62. Antimicrobial Agents and Chemotherapy; 2018.
Zamudio R, Boerlin P, Beyrouthy R, Madec JY, Schwarz S, Mulvey MR et al. Dynamics of extended-spectrum cephalosporin resistance genes in Escherichia coli from Europe and North America. Nat Commun. 2022;13(1).
Poirel L, Madec JY, Lupo A, Schink AK, Kieffer N, Nordmann P et al. Antimicrobial Resistance in Escherichia coli. Microbiol Spectr. 2018;6(4).
Ling W, Furuya-Kanamori L, Ezure Y, Harris PNA, Paterson DL. Adverse clinical outcomes associated with infections by Enterobacterales producing ESBL (ESBL-E): A systematic review and meta-analysis. JAC Antimicrob Resist. 2021;3(2).
Butler MS, Gigante V, Sati H, Paulin S, Al-Sulaiman L, Rex JH, et al. Analysis of the Clinical Pipeline of Treatments for Drug-Resistant Bacterial Infections: Despite Progress, More Action Is Needed. Volume 66. Antimicrobial Agents and Chemotherapy; 2022.
Dolejska M, Papagiannitsis CC. Plasmid-mediated resistance is going wild. Plasmid [Internet]. 2018;99:99–111. https://linkinghub.elsevier.com/retrieve/pii/S0147619X18300349
Global Action Plan on Antimicrobial Resistance. Microbe Magazine. 2015;10(9).
Ministry of Health NAP and M and, Rwanda E. National action plan on antimicrobial resistance 2020–2024, NAPAMR (2021). [cited 2023 May 30]; https://www.who.int/publications/m/item/rwanda-national-action-plan-on-antimicrobial-resistance-2020-2024
Manishimwe R, Nishimwe K, Ojok L. Assessment of antibiotic use in farm animals in Rwanda. Trop Anim Health Prod. 2017;49(6):1101–6.
CLSI. M100 Performance Standards for Antimicrobial Susceptibility Testing, 32nd Edition. Clinical Laboratory Standard Institute. 2021.
Grünzweil OM, Palmer L, Cabal A, Szostak MP, Ruppitsch W, Kornschober C et al. Presence of β-lactamase-producing enterobacterales and salmonella isolates in marine mammals. Int J Mol Sci. 2021;22(11).
Lepuschitz S, Huhulescu S, Hyden P, Springer B, Rattei T, Allerberger F et al. Characterization of a community-acquired-MRSA USA300 isolate from a river sample in Austria and whole genome sequence based comparison to a diverse collection of USA300 isolates. Sci Rep. 2018;8(1).
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5).
Richter M, Rosselló-Móra R, Oliver Glöckner F, Peplies J. JSpeciesWS: A web server for prokaryotic species circumscription based on pairwise genome comparison. Bioinformatics. 2016;32(6).
Beghain J, Bridier-Nahmias A, Le Nagard H, Denamur E, Clermont O. ClermonTyping: an easy-to-use and accurate in silico method for Escherichia genus strain phylotyping. Microb Genom. 2018;4(7).
Clermont O, Dixit OVA, Vangchhia B, Condamine B, Dion S, Bridier-Nahmias A, et al. Characterization and rapid identification of phylogroup G in Escherichia coli, a lineage with high virulence and antibiotic resistance potential. Environ Microbiol. 2019;21(8):3107–17.
Roer L, Johannesen TB, Hansen F, Stegger M, Tchesnokova V, Sokurenko E et al. CHTyper, a Web Tool for Subtyping of Extraintestinal Pathogenic Escherichia coli Based on the fumC and fimH Alleles. J Clin Microbiol. 2018;56(4).
Zhou Z, Alikhan NF, Mohamed K, Fan Y, Achtman M. The EnteroBase user’s guide, with case studies on Salmonella transmissions, Yersinia pestis phylogeny, and Escherichia core genomic diversity. Genome Res. 2020;30(1):138–52.
Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications [version 1; referees: 2 approved]. Wellcome Open Res. 2018;3.
Seemann T. ABRicate. 2020.
Bortolaia V, Kaas RS, Ruppe E, Roberts MC, Schwarz S, Cattoir V et al. ResFinder 4.0 for predictions of phenotypes from genotypes. J Antimicrob Chemother. 2020;75(12).
Liu B, Zheng D, Zhou S, Chen L, Yang J. VFDB 2022: a general classification scheme for bacterial virulence factors. Nucleic Acids Res. 2022;50(D1):D912–7.
Tetzschner AMM, Johnson JR, Johnston BD, Lund O, Scheutz F. In Silico genotyping of Escherichia coli isolates for extraintestinal virulence genes by use of whole-genome sequencing data. J Clin Microbiol. 2020;58(10).
Carattoli A, Zankari E, Garciá-Fernández A, Larsen MV, Lund O, Villa L et al. In Silico detection and typing of plasmids using plasmidfinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58(7).
Naas T, Oueslati S, Bonnin RA, Dabos ML, Zavala A, Dortet L, et al. Beta-lactamase database (BLDB) - structure and function. J Enzyme Inhib Med Chem. 2017;32(1):917–9.
Arredondo-Alonso S, Rogers MRC, Braat JC, Verschuuren TD, Top J, Corander J et al. Mlplasmids: A user-friendly tool to predict plasmid- and chromosome-derived sequences for single species. Microb Genom. 2018;4(11).
Wattam AR, Abraham D, Dalay O, Disz TL, Driscoll T, Gabbard JL, et al. PATRIC, the bacterial bioinformatics database and analysis resource. Nucleic Acids Res. 2014;42:D1.
R Core Team. R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria}, 2023. 2023 [cited 2023 May 30]; https://www.R-project.org/
Stevenson MNEHCMJSJTRRJRCJSPSPYKJGPSFSKRPJJMRCCASNSARA. _epiR: Tools for the Analysis of Epidemiological Data_. R package version 2.0.58. 2023 [cited 2023 May 30]; https://CRAN.R-project.org/package=epiR
Sweeney MT, Lubbers BV, Schwarz S, Watts JL. Applying definitions for multidrug resistance, extensive drug resistance and pandrug resistance to clinically significant livestock and companion animal bacterial pathogens. J Antimicrob Chemother. 2018;73(6).
Rosine Manishimwe BNRH et al. Consumption Trends in Antibiotics for Veterinary Use in Rwanda: A Retrospective Study between 2019 and 2021, 21 July 2022, PREPRINT (Version 1) available at Research Square. 2022 [cited 2023 May 30]; https://www.researchsquare.com/article/rs-1807218/v1
Nsengimana O, Habarugira G, Ojok L, Ruhagazi D, Kayitare A, Shyaka A. Infectious coryza in a grey crowned crane (Balearica regulorum) recovered from captivity. Vet Med Sci. 2022;8(2).
Eger E, Heiden SE, Korolew K, Bayingana C, Ndoli JM, Sendegeya A et al. Circulation of Extended-Spectrum Beta-Lactamase-Producing Escherichia coli of Pandemic Sequence Types 131, 648, and 410 Among Hospitalized Patients, Caregivers, and the Community in Rwanda. Front Microbiol. 2021;12.
Riley LW. Pandemic lineages of extraintestinal pathogenic Escherichia coli. Clin Microbiol Infect. 2014;20(5):380–90.
Mathers AJ, Peirano G, Pitout JDD. The role of epidemic resistance plasmids and international high- risk clones in the spread of multidrug-resistant Enterobacteriaceae. Clin Microbiol Rev. 2015;28(3).
Castanheira M, Simner PJ, Bradford PA. Extended-spectrum β-lactamases: an update on their characteristics, epidemiology and detection. JAC Antimicrob Resist. 2021;3(3):dlab092.
Heiden SE, Kurz MSE, Bohnert J, Bayingana C, Ndoli JM, Sendegeya A et al. Flies from a tertiary hospital in Rwanda carry multidrug-resistant Gram-negative pathogens including extended-spectrum beta-lactamase-producing E. coli sequence type 131. Antimicrob Resist Infect Control. 2020;9(1).
Yoon EJ, Choi YJ, Kim D, Won D, Choi JR, Jeong SH. Amplification of the Chromosomal bla CTX-M-14 Gene in Escherichia coli Expanding the Spectrum of Resistance under Antimicrobial Pressure. Microbiol Spectr. 2022;10(3).
Liakopoulos A, Mevius D, Ceccarelli D. A review of SHV extended-spectrum β-lactamases: Neglected yet ubiquitous. 7, Front Microbiol. 2016.
Alonso CA, Michael GB, Li J, Somalo S, Simón C, Wang Y, et al. Analysis of blaSHV-12-carrying Escherichia coli clones and plasmids from human, animal and food sources. J Antimicrob Chemother. 2017;72(6):1589–96.
Bihannic M, Ghanbarpour R, Auvray F, Cavalié L, Châtre P, Boury M et al. Identification and detection of three new F17 fimbrial variants in Escherichia coli strains isolated from cattle. Vet Res. 2014;45(1).
Dubreuil JD. EAST1 toxin: An enigmatic molecule associated with sporadic episodes of diarrhea in humans and animals. 57, J Microbiol. 2019.
Tian GB, Huang YM, Fang ZL, Qing Y, Zhang XF, Huang X. CTX-M-137, a hybrid of CTX-M-14-like and CTX-M-15-like β-lactamases identified in an Escherichia coli clinical isolate. J Antimicrob Chemother. 2014;69(8).
Rao L, Lv L, Zeng Z, Chen S, He D, Chen X, et al. Increasing prevalence of extended-spectrum cephalosporin-resistant Escherichia coli in food animals and the diversity of CTX-M genotypes during 2003–2012. Vet Microbiol. 2014;172:3–4.
Clermont O, Christenson JK, Denamur E, Gordon DM. The Clermont Escherichia coli phylo-typing method revisited: Improvement of specificity and detection of new phylo-groups. Environ Microbiol Rep. 2013;5(1).
Johnson JR, Delavari P, Kuskowski M, Stell AL. Phylogenetic distribution of extraintestinal virulence-associated traits in Escherichia coli. J Infect Dis. 2001;183(1).
Navon-Venezia S, Kondratyeva K, Carattoli A. Klebsiella pneumoniae: A major worldwide source and shuttle for antibiotic resistance. Volume 41. FEMS Microbiology Reviews; 2017.
Mshana SE, Hain T, Domann E, Lyamuya EF, Chakraborty T, Imirzalioglu C. Predominance of Klebsiella pneumoniae ST14 carrying CTX-M-15 causing neonatal sepsis in Tanzania. BMC Infect Dis. 2013;13(1).
Decano AG, Pettigrew K, Sabiiti W, Sloan DJ, Neema S, Bazira J et al. Pan-resistome characterization of uropathogenic escherichia coli and klebsiella pneumoniae strains circulating in uganda and kenya, isolated from 2017–2018. Antibiotics. 2021;10(12).
Irenge LM, Ambroise J, Bearzatto B, Durant JF, Bonjean M, Gala JL. Genomic Characterization of Multidrug-Resistant Extended Spectrum β-Lactamase-Producing Klebsiella pneumoniae from Clinical Samples of a Tertiary Hospital in South Kivu Province, Eastern Democratic Republic of Congo. Microorganisms. 2023;11(2).
Liu L, Qin L, Hao S, Lan R, Xu B, Guo Y et al. Lineage, antimicrobial resistance and virulence of citrobacter spp. Pathogens. 2020;9(3).
Dong H, Li Y, Cheng J, Xia Z, Liu W, Yan T et al. Genomic Epidemiology Insights on NDM-Producing Pathogens Revealed the Pivotal Role of Plasmids on bla NDM Transmission. Microbiol Spectr. 2022;10(2).
Acknowledgements
We are grateful to all participating farmers and to the team of the NVVH in Musanze (Flora, Beatrice, Gaudiose, Grace, Jophride, Antoine, Jean Paul, Sylvie, Geraldine, Vital, Adelphine, Pierre, Aron, and Felicien). We would also like to express our thanks to Anna Stöger, Katrin Berghammer, and Michael Steinbrecher for technical assistance.
Funding
This work was partially supported by “Veterinaer Projekt Ruanda – united vets friendship group”.
Author information
Authors and Affiliations
Contributions
Authors’ contributionsEI, HK, OWF and IL conceived and designed the experiments. EI, HK, ACR, CN, EM, WR and IL performed the experiments. EI, HK, MPS, ACR, BS, ADL, TW, ATF, SDB, SS, SM, RE, JS and IL analyzed the data. EI, HK, MPS, ACR, ADL and IL drafted the manuscript. MPS, ADR, OWF, ATF, SDB, SS, SM, RE, JS, WR and IL revised the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The study was discussed, and the rectal swabbing was approved by the Research Screening and Ethical Clearance Committee of the College of Agriculture, Animal Sciences and Veterinary Medicine, University of Rwanda (003/2021/DRI). The study was performed in compliance with institutional guidelines for research on animals and an “informed consent named “Animal Health Survey Informed Consent Form” was obtained and signed from animal owners. All methods were performed in accordance with the relevant intitutaion guidelines and regulations.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic Supplementary Material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Irimaso, E., Keinprecht, H., Szostak, M.P. et al. Survey in ruminants from Rwanda revealed high diversity and prevalence of extended-spectrum cephalosporin-resistant Enterobacterales. BMC Vet Res 20, 523 (2024). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12917-024-04359-3
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12917-024-04359-3