Next Article in Journal
Health in Preconception, Pregnancy and Postpartum Global Alliance: International Network Preconception Research Priorities for the Prevention of Maternal Obesity and Related Pregnancy and Long-Term Complications
Previous Article in Journal
FAD/NADH Dependent Oxidoreductases: From Different Amino Acid Sequences to Similar Protein Shapes for Playing an Ancient Function
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phylogeny, Resistome, and Virulome of Escherichia coli Causing Biliary Tract Infections

by
Ángel Rodríguez-Villodres
1,2,
Rémy A. Bonnin
3,4,
José Manuel Ortiz de la Rosa
1,2,
Rocío Álvarez-Marín
1,2,
Thierry Naas
3,4,5,
Javier Aznar
1,2,6,
Jerónimo Pachón
2,7,
José Antonio Lepe
1,2 and
Younes Smani
1,2,*
1
Clinical Unit of Infectious Diseases, Microbiology and Preventive Medicine, University Hospital Virgen del Rocío, 41013 Seville, Spain
2
Institute of Biomedicine of Seville (IBiS), University Hospital Virgen del Rocío, CSIC, University of Seville, 41013 Seville, Spain
3
LabEx Lermit, EA7361, Université Paris-Sud, Université Paris-Saclay, 91190 Saint-Aubin, France
4
The Evolution and Ecology of Resistance to Antibiotics” Unit, Institut Pasteur-APHP-Université Paris Sud, 91400 Orsay, France
5
Bacteriology Hygiene Unit, APHP, Hôpital Bicêtre, Le Kremlin-Bicêtre, 94270 Paris, France
6
Department of Microbiology, University of Seville, 41009 Seville, Spain
7
Department of Medicine, University of Seville, 41009 Seville, Spain
*
Author to whom correspondence should be addressed.
J. Clin. Med. 2019, 8(12), 2118; https://doi.org/10.3390/jcm8122118
Submission received: 4 November 2019 / Revised: 27 November 2019 / Accepted: 29 November 2019 / Published: 2 December 2019
(This article belongs to the Section Infectious Diseases)

Abstract

:
Escherichia coli is the most frequent Gram-negative bacilli involved in intra-abdominal infections. However, despite high mortality rates associated with biliary tract infections due to E. coli, there is no study focusing on this pathogen. In this study, we have characterized a group of 15 E. coli isolates obtained from 12 patients with biliary tract infections. Demographic and clinical data of the patients were recovered. Phylogeny, resistome, and virulome analysis through whole genome sequencing and biofilm formation were investigated. Among the 15 E. coli isolates, no predominant sequence type (ST) was identified, although 3 of them belonged to unknown STs (20%). Resistance to ampicillin, amoxicillin/clavulanic acid, cotrimoxazole, and quinolones was more present in these isolates; whereas, third and fourth generation cephalosporins, carbapenems, amikacin, tigecycline, and colistin were highly active. Moreover, high diversity of virulence factors has been found, with sfa, fimH, and gad the most frequently detected genes. Interestingly, 26.6% of the E. coli isolates were high biofilm-producers. Altogether, our data characterized for the first time E. coli isolates associated with biliary tract infections in terms of genomic relationship, resistome, and virulome.

1. Introduction

Escherichia coli is a Gram-negative bacterium with high clinical relevance. This microorganism may cause severe community- and hospital-acquired infections including bacteraemia, urinary tract, respiratory, and intra-abdominal infections [1]. Among intra-abdominal infections, acute cholangitis, an infection of the biliary system, is usually associated with high morbidity and mortality (5–13%) [2,3], reaching up to 29% in cases of malignant obstruction [3]. In case of biliary tract infection, E. coli was reported to be the most frequent pathogen isolated from bile samples (23%) [4].
E. coli sequence type 131 (ST131) identified by multilocus sequence typing (MLST) has been reported as the most prevalent clonal group worldwide, frequently associated with multidrug-resistance (MDR) and infections [5]. Other pandemic lineages such as ST69, ST95, and ST73 have been also associated with more virulent E. coli isolates [6,7]. Focusing at the abdominal level, the phylogenetic group B2 of E. coli non ST131 has been described in inflammatory bowel disease [8]. However, the epidemiology of E. coli causing biliary tract infection is poorly documented.
E. coli needs to move from the large intestine to the biliary tract to produce a biliary tract infection. This process depends on virulence factors that are usually located in pathogenicity islands and can be divided into five main groups: 1) adhesins, 2) toxins, 3) siderophores, 4) capsular, and 5) protectins and invasins [9]. Thus, some genes coding for virulence factors have been associated with urinary tract infections and bacteraemia such as pap, fimH, sfa, iha, bfp (adhesins), hlyA, cnf1, sat (toxins), fyuA, iutA (siderophores), kpsMTII (capsule) and aer, traT, ompT, usp, and malX (miscellaneous proteins) [10,11,12,13,14]. However, to our knowledge, only limited studies investigated the microbiological features of E. coli causing biliary tract infections [15,16], despite the fact that E. coli is the most frequent etiological agent [17]. Characterizing the virulence factors that facilitate the establishment of biliary tract infections in E. coli may be helpful in order to identify potential targets that could be locked as a therapeutic strategy. To this end, we examined the clonal relationship, the resistome, and virulome in 15 E. coli isolates responsible for biliary tract infections in Spain.

2. Materials and Methods

2.1. Bacterial Strains

Fifteen E. coli isolates obtained from bile and blood of 12 patients with biliary tract infection (BTI), and 21 isolates of 20 patients with non-BTI (bacteremia and non-biliary intraabdominal infections) hospitalized at the University Hospital Virgen del Rocío (Seville, Spain) were included in this study. The isolates were collected at the Clinical Microbiology Service from blood and/or bile samples of these patients and then, stored at −80 °C in Luria Bertani (LB) broth supplemented with 30% glycerol. Identification of the isolates was performed using MALDI-TOF (Bruker Daltonik GmbH, Leipzig, Germany), as described previously [18] and whole genome sequencing (WGS) (HiSeq systems, Illumina, USA). E. coli ATCC 10536 strain was used as the positive control for biofilm formation. E. coli ATCC 25922 was used as the control strain for antibiotic susceptibility testing.

2.2. Demographics, Clinical Data, and Follow Up

The diagnosis of the biliary tract infection was made by members of the Infectious Disease Service and/or the Intensive Care Unit (ICU), according to the defined criteria [19]. The following variables were collected from the 12 BTI and 20 non-BTI patients: age and gender, Charlson score, an index that categorize comorbidity of patients based on the International Classification of Diseases diagnosis code [20] and McCabe score, a score that can obtain comparisons regarding the importance of host factors based on the severity of the underlying disease [21], the acquisition type of the BTI and non-BTI (community, healthcare-associated, or hospital), severity (sepsis or septic shock) [22], antibiotic exposure in the previous 2 months of hospitalization and duration of the antimicrobial treatment. The patients were followed until hospital discharge, death, or 30 days, which ever occurred first. The study was approved by the Ethics Committee of the University Hospital Virgen del Rocío and University Hospital of Virgen Macarena of Seville (approval no. 0023-N-16, 01/10/2016). Written informed consent was signed by all patients before inclusion in the study.

2.3. Antimicrobial Resistant Profile

The minimum inhibitory concentration (MIC) was determined by broth microdilution method using MicroScan Walk Away NM44 panels (Beckman Coulter, Inc., USA) for the following antibiotics: ampicillin, amoxicillin-clavulanic acid, piperacillin-tazobactam, cefuroxime, cefotaxime, ceftazidime, cefepime, imipenem, meropenem, ciprofloxacin, levofloxacin, gentamicin, tobramycin, amikacin, cotrimoxazole, and tigecycline. Colistin was tested by standard broth microdilution as recommended by the European Committee of Antimicrobial Susceptibility Testing (EUCAST) [23] and the interpretation criteria was performed following the EUCAST breakpoints [24]. MDR criteria were established as non-susceptible to ≥1 agent in ≥3 antimicrobial categories [25]. E. coli ATCC 25922 was used as the control strain.

2.4. Whole Genome Sequencing

Genomic DNA was extracted using an UltraClean microbial DNA isolation kit (MO BIO Laboratories, Mo-Bio, Sait-Quentin en Yvelines, France) from overnight cultures in LB agar (Bio-Rad, Marnes-la-Coquette, France). Genomic DNA quantification was performed using a Qubit fluorometer (Life Technologies, Carlsbad, CA, USA) and adjusted to 0.2 ng/μL. The DNA libraries were prepared using the NexteraXT v3 kit (Illumina, San Diego, CA, USA), according to the manufacturer’s instructions, and then run on the HiSeq systems (Illumina, USA) to generate paired-end 150-bp reads. De novo assembly of Illumina reads was performed using CLC genomic workbench 10.1 according the manufacturer’s recommendations (Qiagen, Courtaboeuf, France). The genome was annotated using the Rapid Annotations using Subsystems Technology (RAST) tool. The acquired antimicrobial resistance genes were identified by uploading assembled genomes to the Resfinder server v2.1 (http://cge.cbs.dtu.dk/services/ResFinder-2.1) [26]. Virulence genes were searched using https://cge.cbs.dtu.dk/services/VirulenceFinder/.

2.5. Phylogenetic Analysis

Sequence alignment and phylogenetic trees were performed using the software MEGA7 and Evolview online tool [27,28]. Neighbor-joining trees were built from concatenated sequences of the 7 housekeeping genes (adk, fumC, gyrB, icd, mdh, purA, and recA) used previously in MLST assay [29], and obtained for whole genome sequencing analysis. The evolutionary distances were computed using Kimura’s two-parameter model with gamma-distributed rate variation of 0.8 [30]. A bootstrap consensus tree inferred from 1000 replicates was depicted to represent the evolutionary history of the taxa analyzed [31].

2.6. Biofilm Formation Assay

An abiotic solid surface biofilm formation assay was performed as described previously [32]. In brief, overnight cultures of the clinical E. coli isolates were diluted 1:100 in fresh LB broth in 96 well plates without shaking and incubated at 37 °C for 48 h. Biofilm was stained with crystal violet 0.4% (v/v) and quantified at 580 nm after solubilization with ethanol 95%. E. coli isolates were classified as biofilm-formers if they yielded optical density at 580 nm (OD580nm) values that were at least twice those of the negative control. E. coli ATCC 10536 strain was used as positive control.

2.7. Statistical Analysis

A descriptive analysis was performed for demographics and clinical variables of the patients included in the study, with median and interquartile range for the quantitative variables and frequency distribution (%) for the qualitative variables. Fisher and χ2 tests were used for categorical variables and continuous variables were analyzed using 2-sample t test or Mann Whitney U test. For biofilm formation assay, the group data are presented as mean ± SEM. Differences were considered significant at p < 0.05. All statistical analyses were performed using SPSS software, version 23.0 (IBM Corporation, Somers, New York, USA).

2.8. Nucleotide Sequence Accession Number

The WGS of the E. coli isolates generated in the study were deposited in GenBank under the BioProject accession number PRJNA557044.

3. Results

3.1. Bacterial Isolates

MALDI-TOF identified the 15 bacterial isolates from the 12 patients with bile infections, following the clinical criteria established by Solomkin et al. [19], as E. coli. Among them, 11 (73.3%) were isolated from bile samples and 4 (26.6%) from blood cultures. Three patients had two isolates each: 140-HE and 23-AE from blood and bile cultures, 43-HE and 4-AE from blood and bile cultures, and 60-AE and 61-AE from the bile culture, respectively.

3.2. Demographics and Clinical Data Analysis

Demographic and clinical features of the 12 and 20 patients with biliary tract infections (BTI) and non-BTI, respectively, included in this study are detailed in Table 1. To summarize, for both group of patients, the demographic (age and gender), comorbidities (Charlson and McCabe scores), sepsis or septic shock, and previous treatment and days of treatment were not significantly different. Regarding to the acquisition of the infection, significant difference was observed between both groups (p = 0.035). Importantly, there was no difference in mortality between both groups.

3.3. Antimicrobial Resistance Profile

Antimicrobial susceptibility testing results are shown in Table 2. Non-significant difference has been observed between BTI and non-BTI isolates. The most frequent resistance found in BTI isolates was by far Ampicillin (73.3%), followed by amoxicillin-clavulanic acid (26.6%), cotrimoxazole (26.6%), and fluoroquinolones (13.3%). The rest of the antibiotics (third and fourth generation cephalosporins, carbapenems, amikacin, tigecycline, and colistin) were highly active against these E. coli isolates. For the non-BTI isolates, the most frequent resistance was by far Ampicillin (90.4%), followed by fluorquinolones (33.33%), amoxicillin-clavulanic acid (28.6%), cotrimoxazole (28.5%), and cefuroxime (14.3%). The rest of the antibiotics were highly active.
A MDR (resistance to beta-lactams including or not cephalosporins, beta-lactams/beta-lactamases inhibitors, fluoroquinolones, cotrimoxazole, and/or aminoglycoside) pattern was found in 2 and 5 BTI and non-BTI isolates (13.3% and 23.8%), respectively. Extended spectrum beta-lactamase (ESBL) was only detected in the 2 non-BTI isolates; whereas, carbapenemase production were not detected in any of the BTI and non-BTI isolates.

3.4. Genetic Characterization

3.4.1. Epidemiology

The WGS analysis of the E. coli isolates showed that 12 isolates belonged to 10 different STs and 3 isolates (3-AE, 47-HE, and 61-AE) belonged to novel STs (Table 3). Moreover, the isolate pairs recovered from the same patient 140-HE and 23-AE, and 43-HE and 4-AE belonged to ST131 and ST58, respectively. In contrast, the isolate pair 60-AE and 61-AE corresponded to two different STs, namely ST542 and an unknown ST, respectively.
The phylogenetic analysis of the concatenated sequences of the seven MLST housekeeping genes of each isolate showed a tree consisting of three clusters and one branch formed by only one isolate that belonged to ST69 (66-AE) and unrelated phylogenetically to the other isolates. The cluster 1 contained 5 isolates; two of them (140-HE and 23-AE) belonged to ST131 and 1 (3-AE) to an unknown ST related to ST12 complex. The cluster 2 contained 3 isolates; one of them (61-AE) belonged to an unknown ST related to ST542. Finally, the cluster 3 had 6 isolates; one of them (47-HE) belonged to an unknown ST related to ST3640 (Figure 1).

3.4.2. Virulome

The WGS analysis has identified several genes coding for virulence factors in the E. coli isolates (Table 4). The most frequently detected genes were sfa (S fimbrial adhesin, 93.3%), fimH (type 1 Fimbriae, 93.3%), and gad (glutamate decarboxylase, 86.6%), followed by lpfA (long polar fimbria, 60.0%), iss (increased serum survival, 53.3%), iroN (enterobactin siderophore receptor protein, 53.3%), iutA (ferric aerobactin receptor, 53.3%), mchF (ABC transporter protein, 46.6%), mchC (microcin H47 maduration system, 20.0%), senB (enterotoxin, 20.0%), mchB (microcin H47, 13.3%), astA (heat-stable enterotoxin-1, 13.3%), iha (irgA homologue adhesion, 13.3%), mcmA (microcin M, 13.3%), vat (vacuolating autotransporter toxin, 13.3%), papG (P fimbrial adhesin, 13,3%), tsh (temperature-sensitive hemagglutinin, 6.6%), air (enteroaggregative immunoglobulin repeat protein, 6.6%), eilA (Salmonella hilA homologue, 6.6%), ireA (siderophore receptor, 6.6%), cnf1 (cytotoxic necrotizing factor-1, 6.6%), and ipaH (invasion plasmid antigen, 6.6%). Other relevant virulence factors in E. coli such as hlyA (hemolysin A), stx-1 and stx-2 (Shiga toxins-1 and -2), cdtB (cytolethal distending toxin), eaeA (intimin), and bfp (bundle-forming pilus) were not found in these isolates.

3.4.3. Resistome

The WGS revealed that the blaTEM-1a and blaTEM-1b genes, coding for the broad spectrum beta-lactamase TEM-1, were the most frequently detected enzyme (66.6%). Mutations in gyrA, parC, and parE genes, associated with fluoroquinolone resistance, were detected in 33%, 6.6%, and 20% of the isolates, respectively. Resistance to aminoglycosides (aph, 40%; aadA, 33.3%; and strA, 13.3%), sulphonamides (sul, 46.6%), trimethoprim (dfrA, 33.3%), tetracycline (tet, 20.0%), azithromycin (mph, 13.3%), and chloramphenicol (catA, 13.3%) were also detected in these isolates (Table 5).

3.4.4. Biofilm Formation

The analysis of the biofilm formation in abiotic surface showed that 4 isolates (26.6%) were biofilm-producers with OD580nm higher than 0.2. Two of them belonged to unknown ST (3-AE and 61-AE), one to ST2230 (1-HE), and one to ST542 (60-AE) (Figure 2).

4. Discussion

This study shows, for the first time, the genomic relationship, resistome, and virulome of 15 E. coli isolates obtained from bile and blood cultures of 12 patients with biliary tract infection. Twenty percent of E. coli isolates belonged to unknown ST but with different allele combinations. This percentage is higher than those found in other areas of infections such as urinary tract infection (3.5%), bloodstream infection (10.5%), and left-sided inflammatory bowel disease (0%) [6,8,33]. Of note, twenty-one E. coli isolates recovered from different non-biliary sources in our hospital and over the same period of time were also characterized phylogenetically by MLST and none of them belonged to unknown STs (Table S1).
Although, ST131, ST95, ST73, and ST69 are clearly predominant STs in human E. coli infections such as urinary tract or bloodstream infections [7,34], in our study only 2 E. coli isolates (140-HE and 23-AE, from the same patient) belonged to ST131 and another one (66-AE), to ST-69. The rest of the isolates belonged to other STs, which suggests, even though the sample size of this study is small, that in E. coli biliary tract infections, there is not a predominant ST. It is noteworthy to mention that 140-HE and 23-AE isolates that belonged to ST131 had similar resistance profiles (resistance to ampicillin, amoxicillin-clavulanic acid, and cotrimoxazole) with exception for piperacillin-tazobactam. Indeed, the isolates 140-HE and 23-AE were categorized susceptible to piperacillin-tazobactam by MicroScan and resistant to piperacillin-tazobactam by broth microdilution method. This could be due to the presence of different heteroresistant populations that cannot be detected by the automated Microscan system [35].
Regarding the antimicrobial resistance profiles of the 15 E. coli isolates of this study, ampicillin showed the highest percentage of resistance (73.3%), followed by amoxicillin-clavulanic acid (26.6%), cotrimoxazole (26.6%), and fluoroquinolones (13.3%). These data are in agreement with those of WGS analysis (Table 5). Specifically, blaTEM genes involved in ampicillin resistance have been found in 66.6% of the isolates, dfrA gene involved in trimethoprim resistance has been found in 33.3% of the isolates, whereas gyrA and parE were more prevalent (33% and 20%) than the observed fluoroquinolones resistance rate (13.3%). However, Razaghi et al. reported that 54%, 31.8%, 22.7%, and 0% of E. coli isolated from bile presented resistance to ciprofloxacin, meropenem, ceftazidime, and amoxicillin-clavulanic acid, respectively [15]. These differences are likely the result of differences in local epidemiology. In Spain, the analysis of a total of 1429 E. coli isolates causing intra-abdominal infections in a surveillance study called SMART (Study for Monitoring Antimicrobial Resistance Trends) between 2016 and 2017 revealed that ciprofloxacin and amoxicillin-clavulanic acid presented resistance rates of 27.6% and 18.5%, respectively [36]. However, in Iran, the MDR pattern was detected in 95.5% of the E. coli recovered from bile, reporting different antimicrobial resistance rates in comparison with our results, mainly in the case of third generation cephalosporins (82% vs. 0%), carbapenems (≈70% vs. 0%), gentamicin (36.4% vs. 6.6%), and ciprofloxacin (45.5% vs. 13.3%), respectively [37].
The WGS showed a high diversity in the virulome among the E. coli isolates analyzed in this study, similarly to Fernández-Romero et al. observations [38]. None of the virulence factors genes detected were present in all of the isolates, indicating the absence of a unique gene essential for the development of E. coli-borne biliary tract infections. However, three of the virulence factors were detected in around 90% of the isolates: sfa and fimH that are involved in the adhesion to the cells and facilitate the penetration of bacteria into the tissues [39], and gad, which expression is relevant for the maintenance of E. coli in acid environments such as the biliary tract [40]. These data suggest that adhesion ability and resistance to bile acids seem to be important to produce a biliary tract infection. Transcriptomic confirmatory analyses will help to determine the involvement of these genes in the infections development, especially in the patients with two isolates recovered from bile and blood, respectively, in which the expression of these genes may change during the course of the infection. Similar data have been reported by Wang et al. [16] who studied the role of virulence factors in the development of E. coli bacteraemia in patients with acute cholangitis. They detected the presence of fimH and iutA in 88% and 50% of the isolates, respectively; whereas, sfa was detected only in 8% of the isolates [16]. Another study reported that bfp and astA were the most frequent virulence factors in E. coli isolated from bile [15]; although they did not analyze the presence of fimH, sfa, or gad in their isolates. It is important to mention, that no transcriptomic analyses were performed in this study to confirm the degree of expression of these virulence genes, even though they might be or not be expressed.
Some genetic interconnections have been found between different isolates. It appears that the isolate 3-AE could be the ancestor of the cluster 1. Some virulence genes (iss, ipfA, mchC, and mchB) present in this isolate were lost when compared with the rest of the isolates of the cluster 1. In contrast, other virulence genes (gad, mcmA, vat, and iutA) were acquired in the rest of the isolates of the cluster 1. Similarly, we found that the isolate 60-AE acquired the genes iss, ipfA, iroN, astA, senB, iha, and iutA, when compared with the isolate 61-AE recovered from the same patient. In the same line, the isolate 47-HE, closely related to the isolate 8-AE, acquired the genes astA and fimH, when compared with the isolate 8-AE.
Notably, we report, here, that the E. coli isolates do not only harbor virulence and resistance genes, but also 26.6% of them were high biofilm-formers. However, the relationship between biofilm formation and virulence or resistance pattern was not observed, and needs inclusion of more isolates in the future.

5. Conclusions

This study characterized, for the first time by WGS, the molecular epidemiological traits, virulome, and resistome of a collection of E. coli isolates from patients with biliary tract infections. Low proportion of E. coli ST131 was found and most E. coli isolates belonged to different STs. Antimicrobial susceptibility testing and WGS matched perfectly. Resistance to penicillins, cotrimoxazole, and fluoroquinolones were the most frequently encountered resistances. Furthermore, a high diversity of virulence factors was evidenced, with however, a predominance of genes involved in adhesion and resistance to biliary acids. A better knowledge of the microbiological features of E. coli causing biliary tract infection is important to improve the diagnosis and treatment of these patients.

Supplementary Materials

The following are available online at https://www.mdpi.com/2077-0383/8/12/2118/s1, Table S1: Typing of the 21 non-biliary E. coli isolates by MLST. The results of the typing of the 21 non-biliary E. coli isolates by MLST has been included in the supplementary material as Table S1.

Author Contributions

Y.S. conceived the study and designed the experiments, analyzed the results and wrote the manuscript. Á.R.-V., J.M.O.d.l.R, R.Á.-M., R.A.B., performed the experiments and analyzed the results. T.N., J.A., J.P., and J.A.L. reviewed the manuscript.

Funding

This work was supported by the Miguel Servet Tipo I Project grant, Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Ciencia, Innovación y Universidades (CP15/00132), and by Plan Nacional de I+D+i 2013-2016 and Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Ciencia, Innovación y Universidades, Spanish Network for Research in Infectious Diseases (RD16/0016/0009), cofinanced by the European Development Regional Fund “A Way to Achieve Europe,” Operative Program Intelligent Growth 2014-2020. This study was supported by the Assistance Publique – Hôpitaux de Paris (AP-HP), the University Paris-Sud, the Laboratory of Excellence in Research on Medication and Innovative Therapeutics (LERMIT) supported by a grant from the French National Research Agency (ANR-10-LABX-33). Y.S. is supported by the Subprograma Miguel Servet Tipo I, Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Ciencia, Innovación y Universidades, Spain (CP15/00132). A.R.V. is supported by the Subprograma Río Hortega, Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Ciencia, Innovación y Universidades, Spain (CM18/00122). R.A.M. is supported by the Subprograma Juan Rodés, Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Ciencia, Innovación y Universidades, Spain (JR17/00025).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rogers, B.A.; Sidjabat, H.E.; Paterson, D.L. Escherichia coli O25b-ST131: A pandemic, multiresistant, community-associated strain. J. Antimicrob. Chemother. 2011, 66, 1–14. [Google Scholar] [CrossRef] [PubMed]
  2. Lai, E.C.S.; Mok, F.P.T.; Tan, E.S.Y.; Lo, C.; Fan, S.; You, K.; Wong, J. Endoscopic biliary drainage for severe acute cholangitis. N. Engl. J. Med. 1992, 326, 1582–1586. [Google Scholar] [CrossRef] [PubMed]
  3. Englesbe, M.J.; Dawes, L.G. Resistant pathogens in biliary obstruction: Importance of cultures to guide antibiotic therapy. Hpb 2005, 7, 144–148. [Google Scholar] [CrossRef] [PubMed]
  4. Leung, J.W.; Liu, Y.L.; Lau, G.C.T.; Chan, R.C.Y.; Lai, A.C.W.; Ling, T.K.W.; Cheng, A.F. Bacteriologic analyses of bile and brown pigment stones in patients with acute cholangitis. Gastrointest. Endosc. 2001, 54, 340–345. [Google Scholar] [CrossRef] [PubMed]
  5. Nicolas-Chanoine, M.H.; Bertrand, X.; Madec, J.Y. Escherichia coli ST131, an intriguing clonal group. Clin. Microbiol. Rev. 2014, 27, 543–574. [Google Scholar] [CrossRef] [PubMed]
  6. Adams-Sapper, S.; Diep, B.A.; Perdreau-Remington, F.; Riley, L.W. Clonal composition and community clustering of drug-susceptible and -resistant Escherichia coli isolates from bloodstream infections. Antimicrob. Agents Chemother. 2013, 57, 490–497. [Google Scholar] [CrossRef]
  7. Manges, A.R.; Harel, J.; Masson, L.; Edens, T.J.; Portt, A.; Reid-Smith, R.J.; Zhanel, G.G.; Kropinski, A.M.; Boerlin, P. Multilocus sequence typing and virulence gene profiles associated with Escherichia coli from human and animal sources. Foodborne Pathog. Dis. 2015, 12, 302–310. [Google Scholar] [CrossRef]
  8. Petersen, A.M.; Nielsen, E.M.; Litrup, E.; Brynskov, J.; Mirsepasi, H.; Krogfelt, K.A. A phylogenetic group of Escherichia coli associated with active left-sided Inflammatory Bowel Disease. BMC Microbiol. 2009, 9, 171. [Google Scholar] [CrossRef]
  9. Dale, A.P.; Woodford, N. Extra-intestinal pathogenic Escherichia coli (ExPEC): Disease, carriage and clones. J. Infect. 2015, 71, 615–626. [Google Scholar] [CrossRef]
  10. Ron, E.Z. Distribution and evolution of virulence factors in septicemic Escherichia coli. Int. J. Med. Microbiol. 2010, 300, 367–370. [Google Scholar] [CrossRef]
  11. Hung, W.T.; Cheng, M.F.; Tseng, F.C.; Chen, Y.S.; Lee, S.S.J.; Chang, T.H.; Lin, H.H.; Hung, C.H.; Wang, J.L. Bloodstream infection with extended-spectrum beta-lactamase-producing Escherichia coli: The role of virulence genes. J. Microbiol. Immunol. Infect. 2019. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, M.C.; Tseng, C.C.; Wu, A.B.; Huang, J.J.; Sheu, B.S.; Wu, J.J. Different roles of host and bacterial factors in Escherichia coli extra-intestinal infections. Clin. Microbiol. Infect. 2009, 15, 372–379. [Google Scholar] [CrossRef] [PubMed]
  13. Tarchouna, M.; Ferjani, A.; Ben-Selma, W.; Boukadida, J. Distribution of uropathogenic virulence genes in Escherichia coli isolated from patients with urinary tract infection. Int. J. Infect. Dis. 2013, 17, e450–e453. [Google Scholar] [CrossRef] [PubMed]
  14. Ciesielczuk, H.; Betts, J.; Phee, L.; Doumith, M.; Hope, R.; Woodford, N.; Wareham, D.W. Comparative virulence of urinary and bloodstream isolates of extra-intestinal pathogenic Escherichia coli in a Galleria mellonella model. Virulence 2015, 6, 145–151. [Google Scholar] [CrossRef] [PubMed]
  15. Razaghi, M.; Tajeddin, E.; Ganji, L.; Alebouyeh, M.; Alizadeh, A.H.M.; Sadeghi, A.; Zali, M.R. Colonization, resistance to bile, and virulence properties of Escherichia coli strains: Unusual characteristics associated with biliary tract diseases. Microb. Pathog. 2017, 111, 262–268. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, M.; Tseng, C.; Chen, C.; Wu, J.; Huang, J. The Role of bacterial virulence and host factors in patients with Escherichia coli bacteremia who have acute cholangitis or upper urinary tract infection. Clin. Infect. Dis. 2002, 35, 1161–1166. [Google Scholar] [CrossRef] [PubMed]
  17. Kochar, R.; Banerjee, S. Infections of the Biliary Tract. Gastrointest. Endosc. Clin. N. Am. 2013, 23, 199–218. [Google Scholar] [CrossRef]
  18. Cherkaoui, A.; Hibbs, J.; Emonet, S.; Tangomo, M.; Girard, M.; Francois, P.; Schrenzel, J. Comparison of two matrix-assisted laser desorption ionization-time of flight mass spectrometry methods with conventional phenotypic identification for routine identification of bacteria to the species level. J. Clin. Microbiol. 2010, 48, 1169–1175. [Google Scholar] [CrossRef]
  19. Solomkin, J.S.; Mazuski, J.E.; Bradley, J.S.; Rodvold, K.A.; Goldstein, E.J.; Baron, E.J.; O’Neill, P.J.; Chow, A.W.; Dellinger, E.P.; Eachempati, S.R.; et al. Diagnosis and management of complicated intra-abdominal infection in adults and children: Guidelines by the Surgical Infections Society and the infectious Diseases Society of America. Clin. Infect. Dis. 2010, 50, 133–164. [Google Scholar] [CrossRef]
  20. Charlson, M.E.; Pompei, P.; Ales, K.L.; MacKenzie, C.R. A new method of classifying prognostic comorbidity in longitudinal studies: Development and validation. J. Chronic. Dis. 1987, 40, 373–383. [Google Scholar] [CrossRef]
  21. McCabe, W.R.; Jackson, G.G. Gram-Negative bacteriemia. Arch. Intern. Med. 1962, 110, 847–855. [Google Scholar] [CrossRef]
  22. Singer, M.; Deutschman, C.S.; Seymour, C.W.; Shankar-Hari, M.; Annane, K.; Bauer, M.; Bellomo, R.; Bernard, G.R.; Chiche, J.; Coopersmith, C.M.; et al. The Third International Consensus Definitions for Sepsis and Septic Shock (Sepsis-3). JAMA 2016, 315, 801–810. [Google Scholar] [CrossRef] [PubMed]
  23. European Committee on Antimicrobial Susceptibility Testing. Recommendations for MIC Determination of Colistin (Polymyxin E) As Recommended by the Joint CLSI-EUCAST Polymyxin Breakpoints Working Group; European Committee on Antimicrobial Susceptibility Testing: Växjö, Sweden, 2016. [Google Scholar]
  24. European Committee on Antimicrobial Susceptibility Testing. Breakpoint Tables for Interpretation of MICs and Zone Diameters; European Committee on Antimicrobial Susceptibility Testing: Växjö, Sweden, 2019. [Google Scholar]
  25. Magiorakos, A.P.; Srinivasan, A.; Carey, R.B.; Carmeli, Y.; Falagas, M.E.; Giske, C.G.; Harbarth, S.; Hindler, J.F.; Kahlmeter, G.; Olsson-Liljequist, B.; et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: An international expert proposal for interim standard definitions for acquired resistance. Clin. Microbiol. Infect. 2012, 18, 268–281. [Google Scholar] [CrossRef] [PubMed]
  26. Zankari, E.; Hasman, H.; Cosentino, S.; Vestergaard, M.; Rasmussen, S.; Lund, O.; Aarestrup, F.M.; Larsen, M.V. Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 2012, 67, 2640–2644. [Google Scholar] [CrossRef]
  27. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef]
  28. Zhang, H.; Gao, S.; Lercher, M.J.; Hu, S.; Chen, W. EvolView, an online tool for visualizing, annotating and managing phylogenetic trees. Nucleic Acids Res. 2012, 40, W569–W572. [Google Scholar] [CrossRef]
  29. Rodríguez-Villodres, Á.; Gil-Marqués, M.L.; Álvarez-Marín, R.; Bonnin, R.A.; Pachón-Ibáñez, M.E.; Aguilar-Guisado, M.; Naas, T.; Aznar, J.; Pachón, J.; Lepe, J.A.; et al. Extended-spectrum resistance to β-lactams/β-lactamase inhibitors (ESRI) evolved from low-level resistant Escherichia coli. J. Antimicrob. Chemother. 2019. [Google Scholar] [CrossRef]
  30. Kimura, M. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 1980, 16, 111–120. [Google Scholar] [CrossRef]
  31. Felsenstein, J. Confidence limits on phylogenies: An approach using the bootstrap. Evolution 1985, 39, 783–791. [Google Scholar] [CrossRef]
  32. Vila-Farrés, X.; Parra-Millán, R.; Sánchez-Encinales, V.; Varese, M.; Ayerbe-Algaba, R.; Bayó, N.; Guardiola, S.; Pachón-Ibáñez, M.E.; Kotev, M.; García, J.; et al. Combating virulence of Gram-negative bacilli by OmpA inhibition. Sci. Rep. 2017, 7, 14683. [Google Scholar] [CrossRef]
  33. Hertz, F.B.; Nielsen, J.B.; Schønning, K.; Littauer, P.; Knudsen, J.D.; Løbner-Olesen, A.; Frimodt-Møller, N. Population structure of drug-susceptible, resistant and ESBL-producing Escherichia coli from community-acquired urinary tract infections. BMC Microbiol. 2016, 16, 63. [Google Scholar]
  34. Manges, A.R.; Geum, H.M.; Guo, A.; Edens, T.J.; Fibke, C.D.; Pitout, J.D.D. Global extraintestinal pathogenic Escherichia coli (ExPEC) lineages. Clin. Microbiol. Rev. 2019, 32, e00135-18. [Google Scholar] [CrossRef] [PubMed]
  35. Rodríguez-Villodres, Á.; Ortiz de la Rosa, J.M.; Álvarez-Marín, R.; Pachón, J.; Aznar, J.; Lepe, J.A.; Smani, Y. Heteroresistance to piperacillin-tazobactam in clinical isolates of Escherichia coli sequence type 131. Antimicrob. Agents. Chemother. 2017, 62, e01923-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Cantón, R.; Loza, E.; Aznar, J.; Castillo, F.J.; Cercenado, E.; Fraile-Ribot, P.A.; González-Romo, F.; López-Hontangas, J.L.; Rodríguez-Lozano, J.; Suárez-Barrenechea, A.I.; et al. Monitoring the antimicrobial susceptibility of Gram-negative organisms involved in intraabdominal and urinary tract infections recovered during the SMART study (Spain, 2016 and 2017). Rev. Esp. Quimioter. 2019, 32, 145–155. [Google Scholar] [PubMed]
  37. Tajeddin, E.; Sherafat, S.J.; Majidi, M.R.; Alebouyeh, M.; Alizadeh, A.H.; Zali, M.R. Association of diverse bacterial communities in human bile samples with biliary tract disorders: A survey using culture and polymerase chain reaction-denaturing gradient gel electrophoresis methods. Eur. J. Clin. Microbiol. Infect. Dis. 2016, 35, 1331–1339. [Google Scholar] [CrossRef] [PubMed]
  38. Fernández-Romero, N.; Romero-Gómez, M.P.; Mora-Rillo, M.; Rodríguez-Baño, J.; López-Cerero, L.; Pascual, A.; Mingorance, J. Uncoupling between core genome and virulome in extraintestinal pathogenic Escherichia coli. Can. J. Microbiol. 2015, 61, 647–652. [Google Scholar] [CrossRef]
  39. Sarowska, J.; Futoma-Koloch, B.; Jama-Kmiecik, A.; Frej-Madrzak, M.; Ksiazczyk, M.; Bugla-Ploskonska, G.; Choroszy-Krol, I. Virulence factors, prevalence and potential transmission of extraintestinal pathogenic Escherichia coli isolated from different sources: Recent reports. Gut Pathog. 2019, 11, 10. [Google Scholar] [CrossRef] [Green Version]
  40. Shin, S.; Castanie-Cornet, M.P.; Foster, J.W.; Crawford, J.A.; Brinkley, C.; Kaper, J.B. An activator of glutamate decarboxylase genes regulates the expression of enteropathogenic Escherichia coli virulence genes through control of the plasmid-encoded regulator, Per. Mol. Microb. 2001, 41, 1133–1150. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Phylogenetic tree using neighbor-joining method with 1000 bootstrap replicas of the 15 E. coli isolates from the patients with biliary tract infections.
Figure 1. Phylogenetic tree using neighbor-joining method with 1000 bootstrap replicas of the 15 E. coli isolates from the patients with biliary tract infections.
Jcm 08 02118 g001
Figure 2. Biofilm formation by the 15 E. coli isolates isolated from patients with biliary tract infection. Unknown: unknown ST.
Figure 2. Biofilm formation by the 15 E. coli isolates isolated from patients with biliary tract infection. Unknown: unknown ST.
Jcm 08 02118 g002
Table 1. Demographics and clinical data from patients with biliary tract infections (BTI) and non-BTI by E. coli.
Table 1. Demographics and clinical data from patients with biliary tract infections (BTI) and non-BTI by E. coli.
VariableBTI Patients
(n = 12)
Non-BTI Patients (n = 20)P-Value
Age, median (range)64 (58–72.5)64.5 (55.75–77.75)0.893
Gender (female), n (%)7 (58.3)10 (50.0)0.647
Charlson Score, median (range)1.5 (0–3.75)3 (0.25–6)0.146
McCabe Score ultimately or rapidly fatal, n (%)2 (16.6)11 (55.0)0.062
Acquisition, n (%)
Community9 (75.0)6 (30.0)0.035
Healthcare2 (16.6)5 (25.0)
Nosocomial1 (8.3)9 (45.0)
Sepsis or septic shock, n (%)3 (25.0)10 (50.0)0.267
Previous antibiotic exposure, n (%)7 (58.3)13 (65.0)0.724
Days of treatment, median (range)12.5 (5–21)10 (8–13)0.329
Death, n (%)3 (25.0)3 (15.0)0.647
BTI, biliary tract infection.
Table 2. Antimicrobial resistance in E. coli isolates causing BTI and non-BTI.
Table 2. Antimicrobial resistance in E. coli isolates causing BTI and non-BTI.
BTI E. coli (n = 15)Non-BTI E. coli (n = 21)P-Value
Ampicillin, n (%)11 (73.3)19 (90.4)0.210
Amoxicillin-clavulanic acid, n (%)4 (26.6)6 (28.6)1.000
Piperacillin-tazobactam, n (%)2 (13.3)1 (4.7)0.559
Cefuroxime, n (%)1 (6.6)3 (14.3)0.626
Cefotaxime, n (%)0 (0.0)2 (9.5)0.500
Ceftazidime, n (%)0 (0.0)2 (9.5)0.500
Cefepime, n (%)0 (0.0)2 (9.5)0.500
Imipenem, n (%)0 (0.0)0 (0.0)1.000
Meropenem, n (%)0 (0.0)0 (0.0)1.000
Ciprofloxacin, n (%)2 (13.3)7 (33.3)0.252
Levofloxacin, n (%)2 (13.3)7 (33.3)0.252
Gentamicin, n (%)1 (6.6)2 (9.5)1.000
Tobramycin, n (%)1 (6.6)2 (9.5)1.000
Amikacin, n (%)0 (0.0)0 (0.0)1.000
Cotrimoxazole, n (%)4 (26.6)6 (28.5)1.000
Tigecycline, n (%)0 (0.0)0 (0.0)1.000
Colistin, n (%)0 (0.0)0 (0.0)1.000
ESBL, n (%)0 (0.0)2 (9.5)0.500
MDR, n (%)2 (13.3)5 (23.8)0.674
BTI, biliary tract infection; ESBL, extended spectrum beta-lactamase; MDR, multidrug-resistant.
Table 3. Typing of 15 E. coli isolates by multilocus sequence typing.
Table 3. Typing of 15 E. coli isolates by multilocus sequence typing.
StrainSourceadkfumCgyrBicdmdhpurArecASTST Complex
1-HEBlood634612120137223023
43-HE ^Blood644162481458155
47-HEBlood6951518986Unknown-
140-HE *Blood53404713362829131131
23-AE *bile53404713362829131131
3-AEBile13615131610122Unknown-
4-AE ^Bile644162481458155
8-AEBile6233169863640-
16-AEBile1011488821010
24-AEBile1313913161091212
60-AE ~Bile112115128886542-
61-AE ~Bile112404128886Unknown-
66-AEBile21352765546969
81-AEBile13521014172517141-
82-AEBile64141624814155155
^,* E. coli isolated from blood and bile samples of the same patient. ~ E. coli isolated from bile of the same patient. ST: sequence type.
Table 4. Virulence genes from E. coli causing biliary tract infection.
Table 4. Virulence genes from E. coli causing biliary tract infection.
Strain1-HE43-HE ^47-HE140-HE *23-AE *3-AE4-AE ^8-AE16-AE24-AE60-AE ~61-AE ~66-AE81-AE82-AE
SourceBloodBloodBloodBloodBileBileBileBileBileBileBileBileBileBileBile
iss++ + ++ ++ +
ipfA+++ +++ ++ +
mchC + + +
mchB + +
mchF++ + + +++
iroN++ + + ++++
tsh+
gad+++++ ++++++ ++
astA + ++
senB ++ +
iha + +
air +
ailA +
mcmA + +
vat + +
ireA +
papG+ +
sfa++++++ ++++++++
fimH+++++++ +++++++
iutA++ + ++ ++ +
ipaH +
cnf1 +
hlyA
cdtB
stx-1
stx-2
eaeA
bfp
Total1085438336143101298
iss: increased serum survival, lpfA: long polar fimbria, mchC: microcin H47 maduration system, mchB: microcin H47, mchF: ABC transporter protein, iroN: enterobactin siderophore receptor protein, tsh: temperature-sensitive hemagglutinin, gad: glutamate decarboxylase, astA: heat-stable enterotoxin-1, senB: enterotoxin, iha: irgA homologue adhesion, air: enteroaggregative immunoglobulin repeat protein, eilA: Salmonella hilA homologue, mcmA: microcin M, vat: vacuolating autotransporter toxin, ireA: siderophore receptor, papG: P fimbrial adhesin, Sfa: S fimbrial adhesin, fimH: type 1 Fimbriae, iutA: ferric aerobactin receptor, ipaH: invasion plasmid antigen, cnf1: 4cytotoxic necrotizing factor-1, hlyA: hemolysin A, cdtB: cytolethal distending toxin, stx-1 and stx-2: Shiga toxins-1 and -2, eaeA: intimin, and bfp: bundle-forming pilus. ^,* E. coli isolated from blood and bile samples of the same patient. ~ E. coli isolated from bile of the same patient.
Table 5. Resistance genes from E. coli causing biliary tract infection.
Table 5. Resistance genes from E. coli causing biliary tract infection.
StrainSourcestrAblaTEMcatAsultetdfrAaphgyrAparCparEaadAmph
1-HEBlood+1a +
43-HE^Blood 1b
47-HEBlood 1a + +
140-HE *Blood 1b + + +(I529L)++
23-AE *Bile 1b + + +(I529L)++
3-AEBile+1b++++
4-AE ^Bile +(S83L)
8-AEBile
16-AEBile 1b + +
24-AEBile
60-AE ~Bile 1b +(S83L)
61-AE ~Bile + +++(S83L) +(I529L)+
66-AEBile 1b ++ +
81-AEBile
82-AEBile 1b++++++(S83L) (D87N)+(S80I) +
^,* E. coli isolated from blood and bile samples of the same patient. ~ E. coli isolated from bile of the same patient. 1a and 1b: two different blaTEM genes.

Share and Cite

MDPI and ACS Style

Rodríguez-Villodres, Á.; Bonnin, R.A.; Ortiz de la Rosa, J.M.; Álvarez-Marín, R.; Naas, T.; Aznar, J.; Pachón, J.; Lepe, J.A.; Smani, Y. Phylogeny, Resistome, and Virulome of Escherichia coli Causing Biliary Tract Infections. J. Clin. Med. 2019, 8, 2118. https://doi.org/10.3390/jcm8122118

AMA Style

Rodríguez-Villodres Á, Bonnin RA, Ortiz de la Rosa JM, Álvarez-Marín R, Naas T, Aznar J, Pachón J, Lepe JA, Smani Y. Phylogeny, Resistome, and Virulome of Escherichia coli Causing Biliary Tract Infections. Journal of Clinical Medicine. 2019; 8(12):2118. https://doi.org/10.3390/jcm8122118

Chicago/Turabian Style

Rodríguez-Villodres, Ángel, Rémy A. Bonnin, José Manuel Ortiz de la Rosa, Rocío Álvarez-Marín, Thierry Naas, Javier Aznar, Jerónimo Pachón, José Antonio Lepe, and Younes Smani. 2019. "Phylogeny, Resistome, and Virulome of Escherichia coli Causing Biliary Tract Infections" Journal of Clinical Medicine 8, no. 12: 2118. https://doi.org/10.3390/jcm8122118

APA Style

Rodríguez-Villodres, Á., Bonnin, R. A., Ortiz de la Rosa, J. M., Álvarez-Marín, R., Naas, T., Aznar, J., Pachón, J., Lepe, J. A., & Smani, Y. (2019). Phylogeny, Resistome, and Virulome of Escherichia coli Causing Biliary Tract Infections. Journal of Clinical Medicine, 8(12), 2118. https://doi.org/10.3390/jcm8122118

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop