Skip to main content

Advertisement

We're creating a new version of this page. See preview

  • Research article
  • Open Access

Clonal diversity, virulence genes content and subclone status of Escherichia coli sequence type 131: comparative analysis of E. coli ST131 and non-ST131 isolates from Iran

BMC Microbiology201919:117

https://doi.org/10.1186/s12866-019-1493-8

  • Received: 15 December 2018
  • Accepted: 14 May 2019
  • Published:

Abstract

Background

The Escherichia coli sequence type 131 (ST131) is a well established clone causing significant extraintestinal infections worldwide. However, no studies have been reported the phenotypic and molecular traits of ST131 isolates in comparison to other clones of E. coli from Iran. So, we determined the differences between 69 ST131 strains collected during a one year surveillance study and 84 non-ST131 isolates, including 56 clinical fluoroquinolone resistant and 28 broiler colibacillosis isolates in terms of clonality and genetic background.

Results

ST131 isolates were associated with phylogroup B2 (68 out of 69 isolates, 98.4%), while clinical non-ST131 and fluoroquinolone resistant broiler isolates mainly belonged to phylogroup A. The highest virulence score was observed in ST131 clone, while they showed less diversity in virulence profiles than other clinical isolates. Almost all of the ST131 isolates (95.6%) were ExPEC and had the highest virulence scores, but their resistance scores were less than clinical non-ST131 isolates. Broiler isolates showed higher prevalence of ExPEC-associated virulence genes and CTX-M-G1/G9 resistance determinants as compared to clinical non-ST131 isolates. While blaOXA-48/NDM carbapenemases were mostly found in ST131 clone, resistance rate against ertapenem was higher among clinical non-ST131 strains. According to ERIC-based fingerprinting, the ST131 strains were more genetically similar, followed by non-ST131 and broiler isolates.

Conclusions

ST131 isolates possess the ability to make a balance between clonality and extent of resistance/virulence genes content, so this phenomenon gives a fitness advantage over other E. coli clones. The broilers E. coli population poses a potential zoonotic risk which could be transmitted to the community through the food chain. A number of factors are involved in the dissemination of and infections due to ST131 clone.

Keywords

  • ST131
  • Carbapenem
  • NDM
  • OXA-48
  • Avian
  • Virulence genes
  • Antibiotic resistance

Background

Extraintestinal infections caused by Escherichia coli, including urinary tract infections, bloodstream infections and neonatal meningitis are common, costly and are major causes of morbidity and mortality. Emerging antimicrobial resistance complicates the management of these infections [1].

E. coli sequence type 131 (ST131) is defined by multilocus sequence typing (MLST) technique and is characterized by resistance to first line antibiotics such as fluoroquinolone and extended spectrum cephalosporins (ESCs). Two major subclones which were identified within the ST131 population are H30 and H30Rx. Except for rare strains, almost all H30 isolates are fluoroquinolone resistant, while H30Rx subclone has been associated with Extended Spectrum β-lactamases (ESBLs), namely CTX-M-15 [2, 3].

Fluoroquinolones are widely used in farm animals, mainly in broilers [4, 5] and fluoroquinolone resistant E. coli strains are frequently isolated from healthy and diseased birds [4]. It has been hypothesized that avian E. coli isolates are reservoirs of virulence and resistance determinants, and strains harboring these elements may be able to cause foodborne diseases in humans. Characterization of overlapping traits between avian pathogenic E. coli (APEC) and urinary pathogenic E. coli (UPEC), such as O-serogroups, virulence markers and assignments to phylogroups, encouraged this hypothesis and subsequent researches comparing APEC and human E. coli isolates including UPEC, neonatal meningitis E. coli (NMEC) and E. coli causing septicemia showed common molecular traits between human and avian strains [4, 6, 7]. While the presence of ST131 clone among E. coli isolates contaminating raw chicken meat wasn’t confirmed in the studies from Spain and Italy [8, 9], another survey from Spain found this specific clone in 7% of retail chicken samples [10]. Vincent et al. also had identified E. coli ST131 from retail chicken samples in Canada but at a significant lower prevalence (0.4%) than in Spain [11].

The presented reasons for the success of the ST131 pandemic clone are the acquisition of antimicrobial resistance and additional virulence factors as well as its predominance in the human gut. There are some reports which suggest that epidemiological differences exist between CTX-M-producing strains of ST131 and non-ST131 clones [12]. Our previous study from Iran reported important data on clinical E. coli ST131 isolates [13]. However, a comparative work describing fluoroquinolone resistant/ESBL producing ST131 E. coli and other non-ST131 isolates or isolates from avian sources has not been reported so far from Iran. Therefore, the objectives of this study were to compare the clonality, virulence characteristics and resistance genotypes of ST131 isolates with that of human non-ST131 and broiler colibacillosis strains in order to determine whether these genotypes were shared or not between these groups of isolates.

Results

Human-sourced E. coli isolates

Sixty-nine strains out of 338 collected human E. coli isolates belonged to ST131 clonal group. The O25b/O16 serogroups, H30/H30Rx subclones and virotypes of 63 isolates have been previously described [13]. The remaining six newly identified ST131 isolates belonged to phylogroup B2, O25b subgroup, virotype C and H30 subclone. Of the remaining 269 E. coli isolates, 56 fluoroquinolone resistant isolates were randomly selected. Thus, a total of 125 unique human E. coli isolates (isolated from 125 admitted patients) were subjected to the thorough investigation described below. Among the ST131 isolates, 39 (56.5%) were from female patients with a mean age of 65.49 years and 30 (43.5%) were from male patients, with a mean age of 73 years, whereas among non-ST131 isolates 39 (69.6%) were from female patients and 17 (30.4%) were from male patients with mean ages of 66.6 and 64.8 years, respectively. Overall, specimen types for ST131 and clinical non-ST131 isolates included blood (0 and 1.8%), wound (5.8% and 1.8%), respiratory (11.6% and 12.5%) and urine (82.6% and 83.9%). The prevalence of ST131 and clinical non-ST131 isolates were higher among women, however, sex distribution did not vary significantly between two groups of clinical isolates.

Prevalence of ESBL-producing E .coli and antimicrobial susceptibility

Fluoroquinolone resistance was detected among 63 ST131 (all belonged to subgroup O25b) and 18 (64.2%) broiler isolates. The highest frequency of resistance was obtained against cefotaxime (95.7%), cefotaxime (100%) and trimethoprim/sulfamethoxazole (71.4%) among ST131, non-ST131 and broiler isolates, respectively. Overall, except for carbapenems, a high proportion of ESBL producing isolates were non-susceptible to antibiotics tested and were categorized as MDR (98.4%, 121/123, P: 0.001).

The resistance rates against studied antibiotics among non-ST131 isolates were higher than ST131 and the differences were significant for trimethoprim/sulfamethoxazole, amoxicillin/clavulanate, piperacillin-tazobactam, ampicillin/sulbactam, fluoroquinolones and tobramycin (Table 1). The highest resistance score was obtained among non-ST-131 (median: 9, IQR: 3) as compared with ST131 (median: 8, IQR: 4, P: 0.01) and broiler isolates (median: 2, IQR: 2, P: 0.001).
Table 1

Antibiotic resistance patterns and molecular determination of resistance genes in studied isolates

ST131 (n = 69)

Non-ST131 (n = 84)

antibiotics

O25b-ST131

n = 63

O16-ST131

n = 6

Total ST131

Human (n = 56)

broiler (n = 28)

P value

ST131 vs. non-ST131

ST131 vs broilers

Antibiotic resistance rates, n(%)

 Imipenem

1 (1.6)

1 (16.7)

2 (2.9)

0

0

  

 Meropenem

1 (1.6)

0

1 (1.4)

1 (1.8)

0

  

 Ertapenem

8 (12.7)

1 (16.7)

9 (13)

12 (21.4)

0

 

0.05

 Ceftazidime

49 (77.8)

5 (83.3)

54 (78.3)

48 (85.7)

1 (3.6)

 

0.001

 Cefepime

37 (58.7)

4 (66.7)

41 (59.4)

37 (66.1)

1 (3.6)

 

0.001

 Cefotaxime

61 (96.8)

5 (83.3)

66 (95,7)

56 (100)

1 (3.6)

 

0.001

 Ciprofloxacin/levofloxacin

63 (100)

0

63 (91.3)

56 (100)

18 (64.3)

0.03

0.002

 Piperacillin/tazobactam

10 (15.9)

2 (33.3)

12 (17.4)

24 (42.8)

1 (3.6)

0.003

 

 Ampicillin/sulbactam

27 (43.5)

5 (83.3)

32 (46.4)

39 (69.6)

0

0.01

0.001

 Amoxicillin/clavulanate

36 (57.1)

5 (83.3)

41 (59.4)

44 (78.6)

4 (14.3)

0.03

0.001

 Aztreonam

59 (93.7)

5 (83.3)

64 (92.8)

52 (93)

1 (3.6)

 

0.001

 Trimethoprim/sulfameth

44 (69.8)

5 (83.3)

49 (71)

48 (85.7)

20 (71.4)

0.05

 

 Tobramycin

32 (50.8)

1 (16.7)

33 (47.8)

38 (67.8)

3 (10.7)

0.03

0.001

 Amikacin

7 (11.1)

0

7 (10.1)

8 (14.3)

0

  

 Gentamicin

29 (46)

2 (33.3)

31 (44.9)

33 (59)

1 (3.6)

 

0.001

 Resistance score (median)

8

8

8

9

2

0.01

0.001

 ESBL phenotype

61 (96.8)

5 (83.3)

66 (95.6)

56 (100)

1 (3.6)

 

0.001

 MDR

60 (95.2)

5 (83.3)

65 (94.2)

56 (100)

7 (25)

 

0.001

Resistance genes, n(%)

blaOXA-48

22 (34.9)

4 (66.7)

26 (37.7)

4 (7.1)

0

0.001

0.001

blaNDM

3

0

3

0

0

  

blaTEM

41 (65)

4 (66.7)

45 (65.2)

45 (80.4)

18 (66.7)

  

blaSHV

29 (46)

2 (33.3)

31 (44.9)

24 (42.1)

5 (18.5)

 

0.01

blaCTX-M-15

51 (81)

6 (100)

57 (82.6)

38 (67.9)

20 (71.4)

  

 CTX-M-G-1

55 (87.3)

6 (100)

61 (88.4)

47 (83.9)

26 (92.9)

  

 CTX-M-G-2

0

0

0

2 (3.6)

2 (7.1)

  

 CTX-M-G-8

0

3 (50)

3 (4.3)

0

0

  

 CTX-M-G-9

20 (31.7)

0

20 (29)

7 (12.5)

13 (46.4)

0.03

 

 CTX-M-G-25

18 (28.6)

3 (50)

21 (30.4)

5 (8.9)

2 (7.1)

0.004

0.01

aac6-Ib-cr

43 (68.3)

2 (33.3)

45 (65.2)

39 (69.6)

7 (25)

 

0.001

qnrB

9 (14.2)

0

9 (13)

8 (14.5)

3 (10.7)

  

qnrS

6 (9.5)

0

6 (8.7)

11 (19.6)

5 (17.9)

  

ArmA

4 (6.3)

0

4 (5.8)

6 (10.7)

0

  

aac6-Ib

42 (66.7)

2 (33.3)

44 (63.8)

36 (64.3)

6 (21.4)

 

0.001

aac3-IIa

40 (63.5)

3 (50)

43 (62.3)

42 (75.1)

9 (32.1)

 

0.01

Characterization of resistance genes

The predominant resistance gene detected among the three groups of isolates was CTX-M-G1, followed by aac6Ib-cr in ST131 (45, 65.2%), blaTEM- in non-ST131 (45, 80.4%) and blaCTX-M-15 in broiler (20, 71.4%) isolates. The prevalence of blaSHV-, blaCTX-M-15, CTX-M-G1, CTX-M-G8, CTX-M-G9, CTX-M-G25 and blaOXA-48 was higher among ST131 as compared with clinical non-ST131 isolates; however the differences were only significant for the last three resistance determinants. The broiler isolates showed higher prevalence of CTX-M-G1 and CTX-M-G9 than both of two clinical groups (Table 1). The main gene contributing to the ESBL phenotype among the ST131 strains was CTX-M-G1 (60/66, 90.9%), followed by blaCTX-M-15 (54/66, 81.8%), while in non-ST131 isolates, it was CTX-M-G1 (47/56, 83.9%) and blaTEM- (45/56, 80.4%). The blaCTX-M-15 was more frequently found among ST131 (57/69, 82.6%) than non-ST131 (38/56, 67.9%) and broiler (20/28, 71.4%) isolates, however, the differences were not significant. Among the clinical isolates carrying blaCTX-M-15, resistance rates were higher against gentamicin (57.9%, P: 0.01), amoxicillin-clavulanate (73.7%, P, 0.02), cefepime (69.5%, P: 0.005), amikacin (15.8%, P: 0.02) and tobramycin (65.3%, P: 0.001). The blaOXA-48 and blaOXA-48/NDM combination were detected in 23 (33.3%) and three (4.3%) of the ST131 isolates respectively, of those 19 isolates including three blaNDM carrying strains were found susceptible to carbapenems. Two out of the four (7.1%) blaOXA-48 carrying non-ST131 isolates were also categorized as susceptible to carbapenems.

Phylogrouping

All the 153 strains formed seven distinct phylogroups (A, B1, B2, C, D, E and F). Table 2 shows the phylogroup distribution for the E. coli isolates stratified by source and by ciprofloxacin susceptibility pattern. Among the total population, group B2 was the most prevalent phylogroup, corresponding to 55.6% (85/153) of the isolates, followed in prevalence by the other groups: A (17%), E (13.1%), F (5.2%), B1 and C (each 3.9%), and D (1.3%). Almost all ST131 isolates were from group B2 (98.6%, 68/69), whereas the non-ST131 and broiler isolates exhibited diverse phylogroups, with group A being predominant in these two groups. Considering the clinical non-ST131 and broiler isolates, fluoroquinolone resistant strains found mainly belonged to phylogroup A, although the differences were not significant with respect to phylogroup distribution of resistant isolates. None of broiler isolates were identified as ST131 clonal group. Multiplex ST-PCR detected one urine cultured non-ST131 isolate as ST95 (phylogroup A), and broiler isolates did not belong to the investigated five major ST.
Table 2

Phylogenetic group distribution for 153 E. coli strains stratified by source (human and broilers) and by ciprofloxacin susceptibility status

Phylogenetic group

Ciprofloxacin susceptible strains, N = 16

Ciprofloxacin resistant isolates N = 137

ST131 (n = 6)

n (%)

Non-ST131 (n = 0)

Broilers (n = 10)

n (%)

ST131 (n = 63)

n (%)

Non-ST131 (n = 56)

n (%)

Broilers (n = 18)

n (%)

A

2 (20)

15 (26.8)

9 (50)

B1

2 (20)

3 (5.4)

1 (5.6)

B2

6 (100)

4 (40)

62 (98.4)

13 (23.2)

0

C

0

5 (8.9)

1 (5.6)

D

1 (10)

1 (1.8)

0

E

0

13 (23.2)

7 (38.9)

F

1 (10)

1 (1.6)

6 (10.7)

0

Virulence genotyping

Based on the molecular definition of ExPEC, all ST131 isolates except one O16 and two O25b strains were attributed with the status of ExPEC, in contrast to only 42.9% (24/56, P: 0.001) and 75% (21/28, P: 0.006) of clinical non-ST131 and broiler isolates, respectively. Considering the 84 non-ST131 isolates, including both of the clinical and broiler strains, fluoroquinolone susceptible isolates showed higher virulence scores as compared with fluoroquinolone resistant strains (median; 11 vs. 6.50, P: 0.045). The conserved pattern of virulence genes (i.e., gene prevalence > 90%) of the ST131 isolates included fyuA, F10 papA, iutA, chuA, sat, malX, usp and traT. Eighty-five different virulence profiles were found (named 1–85). Twenty of them grouped more than one isolate; profile 26 was the most frequently identified (n = 28), while profile 11 showed the highest virulence score (one ST131 isolate, virotype B, virulence score: 16). Clinical non-ST131 and broiler isolates were spread among 42 and 21 different virulence profiles, respectively, whereas ST131 isolates showed 25 different profiles, with the virulence profiles number 26 (n = 26) and 6 (n = 15) containing the majority of isolates (Fig. 1). A single pattern (pattern 84; fimH30+) was detected between human (non-ST131) and broiler isolates.
Fig. 1
Fig. 1

ERIC-PCR based dendogram of the 152 E.coli isolates and virulence gene profile. Cluster analysis of the Dice similarity indices based on the unweighted pair group method using average linkages (UPGMA) was done to generate a dendogram describing the relationship among the ERIC profiles. (One non-ST131 isolate didn’t amplify in ERIC-PCR, so it is not included in the dendogram image). P, poultry (broilers); ST, ST131; NST, non-ST131; VP, virulence profile number. Filed squares indicate the presence of the gene

Among ST131 isolates, the prevalence of K5 and afaDraBC was associated with O25b (fluoroquinolone resistant) and O16 (fluoroquinolone susceptible) subgroups, respectively. The presence of ompT, hlyF, iss, ibeA, iroN and kpsMT II genes was associated with broiler isolates. The median virulence score and the virulence score range values decreased progressively from a high of 11 (IQR: 3, range 7–16) for ST131 isolates to a low of 8.50 (IQR: 8, range, 0–14) for broiler isolates, to a score of 6 (IQR: 5, range, 1–14) for non-ST131 isolates (P < 0.001) (Table 3).
Table 3

Virulence genotypes of studied isolates

ST131 (n = 69)

Non-ST131 (n = 84)

Adhesions

O25b-ST131

n = 63

O16-ST131

n = 6

Total ST131

Human (n = 56)

Broiler (n = 28)

P value

ST131 vs. non-ST131

ST131 vs avian

Virulence markers, n (%)

papA

27 (42.8)

1 (16.7)

28 (40.6)

10 (17.9)

4 (14.3)

0.007

0.01

F10 papA

61 (96.8)

6 (100)

67 (97.1)

17 (30.4)

0

0.001

0.001

papEF

26 (41.3)

1 (16.7)

27 (39.1)

7 (12.5)

3 (10.7)

0.001

0.007

papC

30 (47.6)

1 (16.7)

31 (44.9)

14 (25)

9 (32.1)

0.02

sfa/focDE

1 (1.6)

0

1 (1.4)

1 (1.8)

0

afaFM955459

5 (7.9)

1 (16.7)

6 (8.7)

0

0

0.03

afa/draBC

7 (11.1)

3 (50)

10 (14.5)

5 (8.9)

1 (3.6)

fimH30

62 (98.4)

0

62 (89.9)

27 (48.2)

17 (60.7)

0.001

0.003

Siderophores

fyuA

63 (100)

6 (100)

69 (100)

44 (78.6)

7 (25)

0.001

0.001

iroN

1 (1.6)

0

1 (1.4)

0

17 (60.7)

0.001*

iutA

62 (98.4)

5 (83.3)

67 (97.1)

39 (69.6)

23 (82.1)

0.001

0.02

chuA

62 (98.4)

6 (100)

68 (98.6)

35 (62.5)

12 (43)

0.001

0.001

Toxins

cnf1

18 (28.6)

0

18 (26.1)

2 (3.6)

0

0.001

0.001

sat

62 (98.4)

5 (83.3)

67 (97.1)

18 (32.1)

1 (3.6)

0.001

0.001

hlyA

20 (31.7)

0

20 (29)

4 (7.1)

0

0.003

0.001

hlyF

1 (1.5)

0

1 (1.4)

6 (10.7)

19 (67.9)

0.04

0.001*

Capsules

kpsM T II-K5

34 (54)

0

34 (49.3)

10 (17.9)

13 (46.4)

0.001

kpsMT II-K2

32 (50.8)

2 (33.3)

34 (49.3)

20 (35.7)

20 (71.4)

0.04*

Miscellaneous

ibeA

0

1 (16.7)

1 (1.4)

1 (1.8)

5 (17.9)

0.007*

PAI (malX)

61 (96.8)

5 (83.3)

66 (95.7)

32 (56.1)

5 (17.9)

0.001

0.001

usp

62 (98.4)

6 (100)

68 (98.6)

10 (17.9)

6 (21.4)

0.001

0.001

Protections

traT

60 (95.2)

6 (100)

66 (95.7)

43 (76.8)

20 (71.4)

0.002

0.002

ompT

1 (1.5)

0

1 (1.4)

6 (10.7)

19 (67.9)

0.04

0.001*

iss

1 (1.5)

0

1 (1.4)

7 (12.5)

18 (64.3)

0.02

0.001*

Virulence score (median)

11

8.50

11

6

8

0.001

0.001

Strains with ≥8 virulence genes

62 (98.4)

6 (100)

68 (98.6)

20 (35.7)

16 (57.1)

0.001

0.001

ExPEC status

62 (98.4)

4 (66.7)

66 (95.6)

24 (42.9)

21 (75)

0.001

0.006

*: indicate the positive association with broiler isolates

Regardless of ST131 status, H30 positive isolates were overwhelmingly fluoroquinolone resistant (105/106, 99.1%) (P: 0.001). H30Rx sublineage was mainly detected among ST131 isolates (27/69, 39.1%), whereas one and four strains of broiler and non-ST131 isolates belonged to H30Rx sublineage, respectively (P: 0.001). Resistant phenotypes which were associated with H30Rx subclone were against cefotaxime, amikacin, tobramycin, gentamycin, fluoroquinolones, amoxicillin-clavulanate and aztreonam. Almost all (31/32, 96.8%) H30Rx isolates were qualified as ExPEC and were associated with CTX-M-G1, CTX-M-G9 (P: 0.01 for both) and aac6-Ib/ aac6Ib-cr (P: 0.001, for both). H30Rx isolates showed higher prevalence of CTX-M-15 gene as compared with H30-non-Rx subgroup, but the difference was not significant (87.5% vs. 71.9%). H30Rx sublineage of ST131 clone was significantly associated with the virulence genes of cnf1, hlyA, k5, kpsMT II, papA, papC and papEF.

Genotyping

The highest similarity rates were observed among ST131 isolates, followed by the non-ST131 strains and then the broiler isolates (Fig. 1). Two out the three blaNDM carrying isolates exhibited > 80% similarity in ERIC-PCR based fingerprinting. Several dissimilar isolates belonging to all the three groups harboring blaCTX-M-15 and/or blaOXA-48 were scattered throughout the dendogram, suggesting possible horizontal gene transfer.

Discussion

Alarming increase in the isolation of CTX-M-15 producing E. coli isolates have been reported from different countries and this phenomenon has been linked to the clonal expansion of ST131 [14, 15]. The association of ST131 with higher ExPEC fraction and multidrug resistance phenotype, together with their ability for widespread dissemination, has led to the worldwide predominance of ST131 among other E. coli populations. Past studies have reported the association of selected factors with ST131 and non-ST131 isolates [16]. However, comprehensive analyses to study the association of several determinants with clinical ST131, non-ST131 and E. coli isolates from broiler subjects, particularly from developing countries like Iran, remain scarce.

To our knowledge, this is the first study of phenotypic and molecular traits of ST131 and its main subclones, including H30 and H30Rx and non-H30 strains within a population of systematically collected E. coli clinical isolates from Iran. We found that among fluoroquinolone resistant isolates, ST131 clone shows priority over non-ST131 isolates including more virulence genes profiles, higher proportion of ExPEC and high similarity for phylogroup B2 background (98.4% vs. 23.2%).

We found a predominance of blaCTX-M-15 (82.6%) among ESBL producing ST131 E. coli isolates. Notably, approximately 60% of non-ST131 and broiler isolates also harbored blaCTX-M-15, providing evidence of the extended prevalence of blaCTX-M-15 and a significant association of ST131 with blaCTX-M-15. Moreover, the blaCTX-M-15 was significantly associated with resistance to first line antibiotics, notably amoxicillin/clavulanate, aminoglycosides and cefepime, and also among clinical non-ST131 isolates with greater resistance scores. This event would help the rising expansion of blaCTX-M-15 carrying strains over other ESBL producing E. coli isolates in the presence of selection pressure due to the corresponding antibiotics [17].

Studies have reported a core set of virulence factors among ST131 strains including (fimH, fyuA, iutA, traT, malX), with occasional variation [14, 18]. In our study, these core virulence determinants were detected in > 89% of ST131 isolates screened, along with F10 papA, chuA and usp. Because the present ST131 strains were significantly more likely to be from group B2 than the non-ST131 and broiler E. coli isolates, they may have fitness priority over other E. coli isolates.

The proportions of fimH, fyuA, iutA, traT, and PAI genes were similar to those in previous reports [19, 20], although the proportions of hlyA/cnf1 and ompT/sat were different from those reported in previous Korean/Indian and U. S studies [2022]. This result suggests that differences in virulence profiles among ST131 strains may depend on their geographical origin, whether the study was hospital/ “long-term care facility” or community based, and the source of isolates.

The high diversity of virulence gene profiles that we describe in our ESBL producers has been previously described [23], but we found that ST131 has less profile diversity, which implies to the selective distribution of certain sublineages within ST131. Furthermore, the similarity of broiler isolates with human pathogenic E. coli based on virulence markers was confirmed in this study, since a total of 75% of broiler isolates belonged to the ExPEC fraction. As also suggested by other studies, certain APEC subgroups, especially a large fraction of phylogroup A isolates may be considered potential zoonotic agents [24].

In this study, evaluation of the potential zoonotic risk of the E. coli isolates from broiler relies on the phylogroup distribution associated with ciprofloxacin resistance. Regarding the phylogroups, overall, the proportion of group A strains of clinical non-ST131 isolates was greater than the other phylogroups. Considering that half of the fluoroquinolone resistant broiler E. coli isolates we tested belonged to group A and that chicken meat has previously been found prevailed by A and B1 isolates [25], this finding suggests the hypothesis that avian might be a source for at least some human ciprofloxacin resistant group A strains. The change in phylogroup background of isolates from group B2 in drug-resistant ExPEC has been previously observed [26]; however, the appearance and expansion of the ciprofloxacin-resistant B2 ST131 isolates which was also dominant among this clone seems to disagree with these observations [27].

Interestingly, we didn’t find any document indicative for the carriage of ST131 isolates by broiler species which is in agreement with previous reports of rare isolation of ST131 positive cases among avian sources [27]. However, the main limitation of our study was the relatively small number of broiler isolates which restricted statistical power. Based on our results, the ST131 clone is associated with human microbiota and also predominates among fluoroquinolone-resistant strains, while it can be present among susceptible isolates, even in a low number.

Conclusion

The present study confirms that the H30 allele is prevalent among fluoroquinolone resistant isolates. Our data suggest that ST131 clone possess more virulence traits but less multidrug resistance patterns than other clinical fluoroquinolone resistant/ESBL-producing isolates. Among ESBL-positive isolates, ST131 isolates—mostly representing the H30-Rx subclone— have higher virulence scores than non-ST131 isolates, implying higher virulence power and thereby the possible reason for their high prevalence. Non-ST131 isolates including both of the clinical and broiler strains showed an inverse relationship between the prevalence of virulence markers and the resistance phenotype in which higher resistance is accompanied by lower frequency of virulence traits. Further work is needed to be done to understand the significance of avian E. coli isolates contributing to human disease.

Methods

Isolate collection

In this one year cross-sectional study (March 2015–March 2016), a total of 338 E. coli isolates were cultured from patients with extra-intestinal infections admitted to Kosar university hospital of Semnan, Iran. The clinical samples were taken as part of standard care for the patients. Sixty-three out of the 338 isolates had been previously described [13]. Moreover, a total of 28 tissue swabs of culled-animals affected with colibacillosis arrived to the laboratory at the same time period. Samples were taken from chickens already sacrificed for diagnostic purposes. The colibacillosis samples were collected from cases submitted to the Diagnostic Service of the Veterinary School of the Semnan University (Semnan Province). Swabs were taken from two broiler farms located in Semnan Province. The samples were plated onto MacConkey agar and incubated overnight at 37 °C [28].

Antibiotic susceptibility testing and phenotypic detection of ESBL producers

Antibiotic susceptibility patterns for collected E. coli isolates were determined by using the standard disc diffusion. The results of disc diffusion were interpreted based on the Clinical Laboratory Standard Institute (CLSI) recommendations [29]. We considered intermediate susceptibility patterns as resistant. Isolates with an inhibition zone of < 23 mm to any of the imipenem, meropenem, or ertapenem were categorized as carbapenem-resistant [29]. The resistance score was the number of antibiotics to which an isolate was resistant. Multidrug-resistant (MDR) was defined according to resistance to at least one representatives of ≥3 antimicrobial classes [30]. Phenotypic combined disc test based on the CLSI recommendations was used to detect the Extended Spectrum β-lactamase (ESBL) producing isolates [29].

Phylogenetic analysis and sequence typing PCR

The CTAB method was used for extraction of DNA templates for PCR [31]. Major E. coli phylogroups (A, B1, B2, C, D, E and F) were determined by quadruplex PCR [32]. For detection of major ExPEC STs (ST131, ST95, ST73, ST69 and ST127), multiplex PCR was used as described previously [33].

Screening of ST131clonal group

Single nucleotide polymorphism (SNP) in mdh and gyrB genes associated with ST131 clone was detected by using of the gyrB/mdh duplex PCR for all isolates belonged to phylogroups B2, D and F. The H30 and H30-Rx subclones of ST131 isolates were detected, and the two most common O genotypes, O25b and O16 were screened by PCR among study ST131 isolates [13].

Virulence genotyping

The presence of 24 putative virulence genes was assessed by multiplex PCR [34]. Isolates which were positive for ≥2 of the five markers, including afaDrBC, papAH, and/or papC, sfa, KpsM II, or iutA were considered as ExPEC, and the virulence score was calculated as described previously [35]. The genes described previously by Johnson et al as the minimal predictors of APEC virulence; iroN, ompT, hlyF, iutA and iss were detected by a multiplex PCR [36].

Detection of resistance genes

Isolates carrying resistance markers, including carbapenemases (blaNDM, blaOXA-48, blaIMP-, blaVIM- and blaKPC) [37], (β-lactamases) (blaTEM-, blaSHV-,blaCTX-M groups 1, 2, 8, 9, 25) [38, 39] and plasmid mediated quinolone resistance (PMQR) (qnrA, qnrB, qnrS and aac6Ib-cr) were detected by multiplex PCR according to previously published methods [40]. The presence of the blaCTX-M-15 was detected using a single PCR [41]. Moreover, aminoglycoside resistance determinants (ARD), including aac6-Ib/aac3-IIa and 16S rRNA methylase genes (armA, rmtB, rmtC) were screened by PCR [42, 43].

ERIC-PCR analysis

Clonal relatedness was established by Enterobacterial repetitive intergenic consensus sequence polymerase chain reaction (ERIC-PCR). BioNumerics software, version 6.1, (Applied Maths, Sint-Martens-Laten, Belgium) was used to analyz the ERIC patterns. The similarities in amplicon profiles were compared using a Dice coefficient at 1% tolerance and 0.5% optimization, and by using the unweighted-pair group method with arithmetic mean clustering (UPGMA) method, with a cut-off of 80% similarity a dendrogram was constructed [44].

Statistical analysis

Continuous variables were described as means and standard deviations or as medians and interquartile ranges (IQRs). Dichotomous variables were described using frequencies and percentages, and they were compared using chi-square test, as appropriate. The criterion for statistical significance was P < 0.05. Data were analyzed with software SPSS version 16.

Abbreviations

APEC: 

Avian pathogenic E.coli

E. coli

Escherichia coli

ECOR: 

E. coli reference collection

ERIC-PCR: 

Enterobacterial repetitive intergenic consensus sequence polymerase chain reaction

ESBL: 

Extended Spectrum β-lactamase

ExPEC: 

Extraintestinal Pathogenic E.coli

MDR: 

Multidrug-resistant

MLST: 

Multilocus sequence typing

NMEC: 

Neonatal meningitis E. coli

ST: 

Sequence type

UPEC: 

Urinary pathogenic E. coli

UPGMA: 

Unweighted-pair group method with arithmetic mean clustering

Declarations

Acknowledgements

None.

Funding

This work was supported fully by Semnan University of Medical Sciences (grant No. 1180). This is a report of a database from research projects entitled; “Detection of Virulence Factors Encoding Genes among sequence type 131 of Escherichia coli isolates recovered from hospitalized patients in Kowsar hospital, Iran”, registered in Semnan University of Medical Sciences, Semnan, Iran.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Authors’ contributions

All authors read and approved the final manuscript, ZH: analyzed the dendogram, MA & ND: did laboratory works, OP: wrote the manuscript, FS: analyzed the SPSS.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Microbiology Department, Faculty of Medicine, Semnan University of Medical Sciences, Semnan, Iran
(2)
Student Research Committee, Faculty of Medicine, Semnan University of Medical Sciences, Semnan, Iran
(3)
Microbiology Department, Faculty of medicine, Kerman University of Medical Sciences, Kerman, Iran
(4)
Social Determinants of Health Research Center, Semnan University of Medical Sciences, Semnan, Iran

References

  1. Talan DA, Takhar SS, Krishnadasan A. Fluoroquinolone-resistant and extended-spectrum β-lactamase-producing Escherichia coli infections in patients with pyelonephritis, United States. Emerg Inf Dis. 2016;22:1594–603.View ArticleGoogle Scholar
  2. Nicolas-Chanoine MH, Bertrand X, Madec JY. Escherichia coli ST131, an intriguing clonal group. Clin Microbiol Rev. 2014;27:543–74.View ArticleGoogle Scholar
  3. Price LB, Johnson JR, Aziz M. The epidemic of extended spectrum-b-lactamase-producing Escherichia coli ST131 is driven by a single highly pathogenic subclone, H30-Rx. MBio. 2013;4:e00377–13.PubMedPubMed CentralGoogle Scholar
  4. Jakobsen L, Spangholm DJ, Pedersen K, Pedersen K, Jensen LB, Emborg HD, Agersø Y, et al. Broiler chickens, broiler chicken meat, pigs and pork as sources of ExPEC related virulence genes and resistance in Escherichia coli isolates from community dwelling humans and UTI patients. Int J Food Microbiol. 2010;142:264–72.View ArticleGoogle Scholar
  5. Ansari-Lari M, Shekarforoush S, Mehrshad S, Safari H. Prevalence and risk factors for Salmonella spp. colonization in broiler flocks in shiraz, southern Iran. Vet Res Forum. 2014;5:65–8.PubMedPubMed CentralGoogle Scholar
  6. Ewers C, Li G, Wilking H, Kieling S, Alt K, Antao EM. Avian pathogenic, uropathogenic, and newborn meningitis-causing Escherichia coli: how closely related are they? Int J Med Microbiol. 2007;297:163–76.View ArticleGoogle Scholar
  7. Moulin-Schouleur M, Répérant M, Laurent S, Brée A, Mignon-Grasteau S, Germon P. Extraintestinal pathogenic Escherichia coli strains of avian and human origin: link between phylogenetic relationships and common virulence patterns. J Clin Microbiol. 2007;45:3366–7336.View ArticleGoogle Scholar
  8. Egea P, Lopez-Cerero L, Torres E, Gomez-Sanchez M, del C, Serrano L, et al. Increased raw poultry meat colonization by extended-spectrum β-lactamase-producing Escherichia coli in the south of Spain. Int J Food Microbiol. 2012;159:69–73.View ArticleGoogle Scholar
  9. Caruso G, Giammanco A, Cardamone C, Oliveri G, Mascarella C, Capra G. Extra-intestinal fluoroquinolone resistant Escherichia coli strains isolated from meat. Biomed Res Int. 2019;18:8714975.Google Scholar
  10. Mora A, Herrera A, Mamani R, Lopez C, Alonso MP, et al. Recent emergence of clonal group O25b:K1:H4-B2-ST131 ibeA strains among Escherichia coli poultry isolates, including CTX-M-9-producing strains, and comparison with clinical human isolates. Appl Environ Microbiol. 2010;76:6991–7.View ArticleGoogle Scholar
  11. Vincent C, Boerlin P, Daignault D, Dozois CM, Dutil L, et al. Food reservoir for Escherichia coli causing urinary tract infections. Emerg Infect Dis. 2010;16:88–95.View ArticleGoogle Scholar
  12. Nicolas-Chanoine MH, Robert J, Vigan M, Laouénan C, Brisse S, Mentré F. Different factors associated with CTX-M-producing ST131 and non ST131 E. coli clinical isolates. PLoS One. 2013;8:e72191.View ArticleGoogle Scholar
  13. Hojabri Z, Mirmohammadkhani M, Kamali F, Ghassemi K, Taghavipour S, Pajand O. Molecular epidemiology of Escherichia coli sequence type 131 and its H30/H30-Rx subclones recovered from extra-intestinal infections: first report of OXA-48 producing ST131 clone from Iran. Eur J Clin Microbiol Infect Dis. 2017;36:1859–66.View ArticleGoogle Scholar
  14. Clark G, Paszkiewicz K, Hale J, Weston V, Constantinidou C, Penn C. Genomic analysis uncovers a phenotypically diverse but genetically homogeneous Escherichia coli ST131 clone circulating in unrelated urinary tract infections. J Antimicrob Chemother. 2012;67:868–77.View ArticleGoogle Scholar
  15. Fasciana T, Giordano G, DI CARLO P, Colomba C, Mascarella C, Tricoli M (2017) Virulence factors and antimicrobial resistance of ESCHERICHIA COLI ST131 in community-onset healthcare-associated infections in SICILY, Italy. Pharmacol Online 1: 12–21.Google Scholar
  16. Chung HC, Lai CH, Lin JN, Huang CK, Liang SH, Chen WF. Bacteremia caused by extended-spectrum-beta lactamase- producing Escherichia coli sequence type ST131 and non- ST131 clones: comparison of demographic data, clinical features, and mortality. Antimicrob Agents Chemother. 2012;56:618–22.View ArticleGoogle Scholar
  17. Johnson JR, Urban C, Weissman SJ, Jorgensen JH, Lewis JS 2nd, Hansen G. Molecular epidemiological analysis of Escherichia coli sequence type ST131 (O25:H4) and bla CTX-M-15 among extended-Spectrum-β- lactamase-Producing E. coli from the United States, 2000 to 2009. Antimicrob Agents Chemother. 2012;56:2364–70.View ArticleGoogle Scholar
  18. Calhau V, Ribeiro G, Mendonca N, Da Silva GJ. Prevalent combination of virulence and plasmidic-encoded resistance in ST 131 Escherichia coli strains. Virulence. 2013;4:726–9.View ArticleGoogle Scholar
  19. Dahbi G, Mora A, Lopez C, Alonso MP, Mamani R, Marzoa J. Emergence of new variants of ST131 clonal group among extraintestinal pathogenic Escherichia coli producing extended-spectrum beta-lactamases. Int J Antimicrob Agents. 2013;42:347–51.View ArticleGoogle Scholar
  20. Cha MK, Kang CI, Kim SH, Cho SY, Ha YE, Wi YM. Comparison of the microbiological characteristics and virulence factors of ST131 and non-ST131 clones among extended-spectrum β-lactamase–producing Escherichia coli causing bacteremia. Diagn Microbiol Infect Dis. 2016;84:102–4.View ArticleGoogle Scholar
  21. Hussain A, Ranjan A, Nandanwar N, Babbar A, Jadhav S, Ahmed N. Genotypic and phenotypic profiles of Escherichia coli isolates belonging to clinical sequence type 131 (ST131), clinical non-ST131, and fecal non-ST131 lineages from India. Antimicrob Agents Chemother. 2014;58:7240–9.View ArticleGoogle Scholar
  22. Han JH, Garrigan C, Johnston B, Nachamkin I, Clabots C, Bilker WB. Epidemiology and characteristics of Escherichia coli sequence type 131 (ST131) from long-term care facility residents colonized intestinally with fluoroquinolone-resistant Escherichia coli. Diagn Microbiol Infect Dis. 2017;87:275–80.View ArticleGoogle Scholar
  23. Merino I, Shaw E, Horcajada JP, Cercenado E, Mirelis B, Pallarés MA. CTX-M-15-H30Rx-ST131 subclone is one of the main causes of healthcare-associated ESBL-producing Escherichia coli bacteraemia of urinary origin in Spain. J Antimicrob Chemother. 2016;71:2125–30.View ArticleGoogle Scholar
  24. Solà-Ginés M, Cameron-Veas K, Badiola I, Dolz R, Majó N, Dahbi G. Diversity of multi-drug resistant avian pathogenic Escherichia coli (APEC) causing outbreaks of Colibacillosis in broilers during 2012 in Spain. PLoS One. 2015;10:e0143191.View ArticleGoogle Scholar
  25. Jakobsen L, Kurbasic A, Skjøt-Rasmussen L. Escherichia coli isolates from broiler chicken meat, broiler chickens, pork, and pigs share phylogroups and antimicrobial resistance with community-dwelling humans and patients with urinary tract infection. Foodborne Pathog Dis. 2010;7:537–47.View ArticleGoogle Scholar
  26. Moreno E, Prats G, Sabate´ M. Quinolone, fluoroquinolone and trimethoprim/sulfamethoxazole resistance in relation to virulence determinants and phylogenetic background among uropathogenic Escherichia coli. J Antimicrob Chemother. 2006;57:204–11.View ArticleGoogle Scholar
  27. Giufrè M, Graziani C, Accogli M, Luzzi I, Busani L, Cerquetti M. Escherichia coli of human and avian origin: detection of clonal groups associated with fluoroquinolone and multidrug resistance in Italy. J Antimicrob Chemother. 2012;67:860–7.View ArticleGoogle Scholar
  28. Heininger A, Binder M, Schmidt S, Unertl K, Botzenhart K, Doring G. PCR and blood culture for detection of Escherichia coli bacteremia in rats. J Clin Microbiol. 1999;37:2479–82.PubMedPubMed CentralGoogle Scholar
  29. Clinical and Laboratory Standards Institute. Performance standards for antimicrobial susceptibility testing; twentieth informational supplement M100S-S26. Wayne, PA: CLSI; 2016.Google Scholar
  30. Johnson JR, Johnson B, Clabots C, Kuskowski MA, Pendyala S, DebRoy C. Escherichia coli sequence type ST131 as an emerging fluoroquinolone-resistant uropathogen among renal transplant recipients. Antimicrob Agents Chemother. 2010;54:546–5450.View ArticleGoogle Scholar
  31. Ghotaslou R, Aghazadeh M, Ahangarzadeh Rezaee M, Moshafi MH, Forootanfar H, Hojabri Z. The prevalence of aminoglycoside-modifying enzymes among coagulase negative staphylococci in Iranian pediatric patients. J Infect Chemother. 2014;20:569–73.View ArticleGoogle Scholar
  32. 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:58–65.View ArticleGoogle Scholar
  33. Hojabri Z, Mirmohammadkhani M, Darabi N, Arab M, Pajand O. Characterization of antibiotic-susceptibility patterns and virulence genes of five major sequence types of Escherichia coli isolates cultured from extraintestinal specimens: a 1-year surveillance study from Iran. Infect and Drug Resist. 2019;12:893–903.View ArticleGoogle Scholar
  34. Johnson JR, Stell AL. Extended virulence genotypes of Escherichia coli strains from patients with Urosepsis in relation to phylogeny and host compromise. J Infect Dis. 2000;181:261–72.View ArticleGoogle Scholar
  35. Johnson JR, Porter S, Thuras P, Castanheira M. The pandemic H30 subclone of sequence type 131 (ST131) as the leading cause of multidrug-resistant Escherichia coli infections in the United States (2011-2012). Open Forum Infect Dis. 2017;4:ofx089.PubMedPubMed CentralGoogle Scholar
  36. Johnson TJ, Wannemuehler Y, Doetkott C, Johnson SJ, Rosenberger SC, Nolan LK. Identification of minimal predictors of avian pathogenic Escherichia coli virulence for use as a rapid diagnostic tool. J Clin Microbiol. 2008;46:3987–96.View ArticleGoogle Scholar
  37. Hojabri Z, Arab M, Darabi N, Kia NS, Lopes BS, Pajand O. Evaluation of the commercial combined disk test and minimum inhibitory concentration (MIC) determination for detection of carbapenemase producers among gram-negative bacilli isolated in a region with high prevalence of blaOXA-48 and blaNDM. Int Microbiol. 2019;22:81–9.View ArticleGoogle Scholar
  38. Hojabri Z, Pajand O, Bonura C, Aleo A, Giammanco A, Mammina C. Molecular epidemiology of Acinetobacter baumannii in Iran: endemic and epidemic spread of multiresistant isolates. J Antimicrob Chemother. 2014;69:2383–7.View ArticleGoogle Scholar
  39. Woodford N, Fagan EJ, Ellington MJ. Multiplex PCR for rapid detection of genes encoding CTX-M extended-spectrum (beta)-lactamases. J Antimicrob Chemother. 2006;57:154–5.View ArticleGoogle Scholar
  40. Ciesielczuk H, Hornsey M, Choi V, Woodford N, Wareham DW. Development and evaluation of a multiplex PCR for eight plasmid-mediated quinolone-resistance determinants. J Med Microbiol. 2013;62:1823–7.View ArticleGoogle Scholar
  41. Blanco J, Mora A, Mamani R, López C, Blanco M, Dahbi G, et al. Four Main Virotypes among extended-Spectrum-β-lactamase- producing isolates of Escherichia coli O25b:H4-B2-ST131: bacterial, epidemiological, and clinical characteristics, vol. 51; 2013. p. 3358–67.Google Scholar
  42. Haldorsen BC, Simonsen GS, Sundsfjord A, Samuelsen O, Norwegian Study Group on Aminoglycoside Resistance. Increased prevalence of aminoglycoside resistance in clinical isolates of Escherichia coli and Klebsiella spp. in Norway is associated with the acquisition of AAC (3)-II and AAC (6′) Ib. Diagn Microbiol Infect Dis. 2014;78:66–9.View ArticleGoogle Scholar
  43. Pajand O, Hojabri Z, Nahaei MR, Hajibonabi F, Pirzadeh T, Aghazadeh M. In vitro activities of tetracyclines against different clones of multidrug-resistant Acinetobacter baumannii isolates from two Iranian hospitals. Int J Antimicrob Agents. 2014;43:476–8.View ArticleGoogle Scholar
  44. McLellan SL, Daniels AD, Salmore AK. Genetic characterization of Escherichia coli populations from host sources of fecal pollution by using DNA fingerprinting. Appl Environ Microbiol. 2003;69:2587–94.View ArticleGoogle Scholar

Copyright

© The Author(s). 2019

Advertisement