Genomic analysis of Escherichia coli strains isolated from diseased chicken in the Czech Republic

doi:10.21203/rs.2.13990/v3

Download PDF

Research article

Genomic analysis of Escherichia coli strains isolated from diseased chicken in the Czech Republic

https://doi.org/10.21203/rs.2.13990/v3

This work is licensed under a CC BY 4.0 License

Journal Publication

published 09 Jun, 2020

Read the published version in BMC Veterinary Research →

You are reading this latest preprint version

Background: Avian pathogenic Escherichia coli (APEC) can cause various extraintestinal infections in poultry, resulting in massive economic losses in poultry industry. In addition, some avian E. coli strains may have zoonotic potential, making poultry a possible source of infection for humans. Due to its extreme genetic diversity, this pathotype remains poorly defined. This study aimed to investigate the diversity of colibacillosis-associated E. coli isolates from Central European countries with a focus on the Czech Republic.

Results: Of 95 clinical isolates subjected to preliminary characterization, 32 were selected for whole-genome sequencing. A multi resistant phenotype was detected in a majority of the sequenced strains with the predominant resistance to β-lactams and quinolones being associated with TEM-type beta-lactamase genes and chromosomal gyrA mutations respectively. The phylogenetic analysis confirmed a great diversity of isolates, that were derived from nearly all phylogenetic groups, with predominace of B2, B1 and C phylogroups. Clusters of closely related isolates within ST23 (phylogroup C) and ST429 (phylogroup B2) indicated a possible local spread of these clones. Besides, the ST429 cluster carried bla_{CMY-2, -59}genes for AmpC beta-lactamase and isolates of both clusters were generally well-equipped with virulence-associated genes, with considerable differences in distribution of certain virulence-associated genes between phylogenetically distant lineages. Other important and potentially zoonotic APEC STs were detected, incl. ST117, ST354 and ST95, showing several molecular features typical for human ExPEC.

Conclusions: The results support the concept of local spread of virulent APEC clones, as well as of zoonotic potential of specific poultry-associated lineages, and highlight the need to investigate the possible source of these pathogenic strains.

Small Animal Medicine

Large Animal Medicine

Avian pathogenic E. coli

Extraintestinal pathogenic E. coli

virulence-associated genes

avian colibacillosis

whole-genome sequencing

Avian colibacillosis is a complex of several localized or systemic syndromes, affecting poultry of all age and production categories. It comprises yolk sac infection and omphalitis, leading to increased mortality rates in newly hatched chicks, cellulitis in broilers or reproductive tract infections in laying hens. Other forms of manifestation include swollen head syndrome (SHS), respiratory infections and septicemia which frequently result in death or chronic forms of infection. Avian colibacillosis thus represents a great economic burden for the poultry industry (1). Despite its importance as a significant cause of disease, the pathogenesis of these infections is not utterly understood. For a long time APEC strains were considered merely opportunistic pathogens, predominantly, but not exclusively associated with O1, O2, O8, O78 and several other serogroups (2). It has been demonstrated that disease-associated E. coli strains encode multiple putative virulence genes and significantly differ from commensals, particularly in the carriage of the ColV plasmid-associated genes, which are considered markers of poultry-adapted pathogenic strains (3-5).

The ability to cause colibacillosis in chicken defines the APEC (avian-pathogenic E. coli) pathotype. However, not every strain isolated from diseased chicken carries typical virulence-associated genes, underlining an opportunistic character of some types of E. coli infections (6). On the other hand, APEC-like strains (carrying APEC-associated virulence traits) can be found also in the gut of healthy chicken (7,8). It has been suggested by Maturana et al. (9) that the APEC population composes of distinct subpathotypes associated with different syndromes, similar to the human extraintestinal pathogenic E. coli (ExPEC). Interestingly, the authors showed that SHS and omphalitis isolates formed two distinct groups differing in virulence, suggesting primary and opportunistic character of those infections, respectively. Similarly, chronic salpingitis-peritonitis syndrome resulting from an ascending infection and an acute peritonitis without salpingitis, probably originating from respiratory infection or gut translocation after a stress insult, can be distinguished in layers (10-12).

There is a close genetic relationship between APEC and human ExPEC. Zoonotic potential of poultry strains has been implicated. ExPEC are the main cause of urinary-tract infections (as so called uropathogenic E. coli, UPEC) in humans and meningitis in neonates (neonatal-meningitis E. coli, NMEC), and are also associated with bacteremia, sepsis, cellulitis and other potentially fatal infections (13). Similar to APEC, these strains are characterized by the presence of various virulence-associated genes, participating in adhesion and colonization of different tissues, invasion of internal organs, iron acquisition and avoiding host’s immune responses. ExPEC are typically associated with the phylogenetic groups B2 and D, in contrast to commensal and intestinal pathogenic strains derived from groups A and B1 (14) and to APEC, which are highly variable in distribution to various phylogenetic groups (15). Although there is no specific set of genes to define the subpathotypes (16,17), APEC, UPEC and NMEC generally form genetically distinct groups. There is, however, a substantial overlap especially within the B2 phylogenetic group, which comprises strains isolated from both humans and chickens, showing high virulence in the chicken infection model and in the neonatal rat meningitis model with low or no host specifity (18,19,17). Moreover, an isolate showing high virulence in the rat meningitis model has been found in faeces of a healthy chicken (5), another finding that suggests an association between poultry and poultry products as a potential source of human pathogens.

Recently, several highly virulent and resistant ExPEC lineages with worldwide distribution have emerged (e.g. ST131, ST95 etc.) (20). Whereas some of them are associated exclusively with human infections, others are frequently isolated from diseased poultry or poultry products (21-26). It is however difficult to assess the real importance of poultry as a source of human infections. Mechanisms of transmission of pathogenic clones through the production chain to humans are very complex and not quite elucidated, as well as the relationships between genetic „arsenal“ of virulence and resistance-associated genes and pathogenesis of the disease. Whole-genome sequencing (WGS) represents a revolutionary tool to study these mechanisms in their complexity (27). Moreover, an immense variability of APEC pathotype and differences in geographic distribution of specific clones underlines the importance of mapping the local situation. While several papers have reported occurrence of highly pathogenic APEC clones in different counries, the information for the Central Europe have been sporadic or is lacking (28).

2.1 Samples collection and preliminary characterization

A total of 95 isolates were subjected to preliminary characterization including serogrouping, antimicrobial resistance (AMR) testing and PCR detection of virulence and antibiotic resistance genes. The disc diffusion test showed that 69.5% were resistant to three or more groups of antimicrobials, which we considered as a criterion of multiresistance. Resistance to ampicillin was recorded in 78 isolates (82.0%), followed by resistance to nalidixic acid (62 isolates; 65.3%), sulphonamides (45; 47.4%) and sulphonamides-trimethoprim (28; 29.5%). Nineteen isolates (20.0%) showed reduced susceptibility to ciprofloxacin (additional file 1 – Figure 1). Using four antisera (O1, O8, O18 and O78), 49 isolates (52%) were typeable, with predominant serogroups O1 (30; 32%) and O8 (13; 14%). Four isolates reacted to O78 and two against O18 antisera.

The bla_TEMgene was detected in 45 isolates (47.4%). Other prevalent genes detected by PCR included tet(A) (26 isolates; 27.4%) and sul2 (28.4%). As for virulence genes, most isolates carried typical APEC plasmid-associated genes iss (75; 78.9%), iroN (73; 76.8%), iut(A) (68; 71.6%), whereas others, plasmid- or chromosome-associated genes, e.g. cvaC (49; 51.6%), frz_orf4 (44; 46.3%), tsh (32; 33.7%) and felA (7; 7.4%) were less prevalent.

2.2 In silico serotyping, MLST and phylogenetic analysis of 32 selected isolates

We found a diversity of serogroups in the collection of 32 isolates subjected to whole genome sequencing. Overall 14 different O types and 16 H types were identified and 10 isolates failed to be typed by WGS (the results are summarized in Table 1, additional file 1). Except for serogroup O8 (7 isolates), the remaining serogroups were only represented by one or two isolates. The predominance of O8 serogroup appeared as a selection bias since only Czech isolates were selected for sequencing. Of the O8 serogroup, six isolates belonged to the O8:H9 serotype, most of them to ST23 type.

The isolates were derived from all phylogroups according to the Clermont scheme (29) except for the group E: group F (3 isolates), B2 (9 isolates), D (2 isolates), clade I (1 isolate), A (4 isolates), C (6 isolates), B1 (7 isolates).

The MLST analysis identified 22 distinct sequence types (please see Table 1 in the supplementary material and Figure 1), most of them represented only by a single isolate (ST352, ST95, ST140, ST354, ST93, ST4110, ST1249, ST1914, ST770, ST2223, ST746, ST1249, ST162, ST1157, ST602, ST1841, ST533, ST7104). Two isolates were typed as ST117 belonging to the phylogenetic group F; ST429 of the B2 group and ST23 of the C group were detected in 4 and 6 isolates, respectively.

The core genome consisted of 2763 genes (55.28 kbp). The phylogenetic tree based on the core SNPs corresponded to the structure of E. coli phylogeny. Groups F, D and clade I were represented by only a few isolates and did not form any distinct clusters except the minor subcluster of the two group D isolates; two ST117 isolates from the F phylogroup were unrelated to other isolate of F phylogroup (ST354) and formed their own distinct clade. In the B2 cluster two subclusters (B2a, B2b) were found; B2b subcluster was formed by four closely related ST429 isolates and one ST4110. Another cluster included isolates from phylogroups A, C and B1. Interestingly, all isolates of the C group belonged to ST23, O8:H9 serotype (with one exception of O78:H9 serotype).

2.3 Identification of resistance genes

bla_TEM-1(9/32; 28.1%) and a combination of bla_TEM-106, ₁₃₅ (6/32; 18.8%) belonged amongst the most prevalent resistance genes. A combination of plasmid-mediated β-lactamase genes bla_{CMY-2, -59} was detected in four isolates (12.5%), three of them belonging to the ST429, the remainig one to ST354. PMQR (plasmid-mediated quinolone resistance) gene qnrS1 was carried by seven isolates (21.8%) within sequence types 23 and 429. Other identified AMR genes were sul1 (7/32; 21.8%), sul2 (8/32; 25.0% ), dfrA14 (1/32; 3.0%), dfrA15 (4/32; 12.5%), dfrA5 (1/32; 3.0%), tet(A) (12/32; 37.5%), tet(B) (2/32; 6.3%), aadA (6/32; 18.8%), aac(3)-VIa (3/32; 9.4%), ant(2)-1a (1/32; 3.0%), aph(3)-1b (6/32; 18.8%), aph(3)-1a (2/32; 6.3%), aph(6)-1d (5/32; 15.6%), catA1 (2/32; 6.3%), floR (1/32; 3.0%) and bla_TEM-30(1/32; 3,0%). In addition, all isolates showed the presence of genes encoding components of various multidrug efflux pumps, participating in resistance to aminoglycosides, macrolides and fluoroquinolones. Except for qnrS1, which is associated with partial resistance to fluoroquinolones, no other PMQR gene was detected. Reduced susceptibility to quinolones in most isolates appeared to be due to chromosomal mutations, especially in the gyrA gene (21; 65.6%), to lesser extent also in parC (5; 15.6%) and parE (1; 3%). In five ST23 isolates (15.6%), a mutation in the ampC promoter was detected. (For overview of resistance genes, please see the Table 1, additional file 1.)

2.4 Identification of virulence genes

The genomic analysis confirmed a great diversity of selected isolates (see figure 2 and supplementary material, file 3). Overall, factors associated with adhesion and invasion, as well as siderophores were found in most isolates; more than 90 % of isolates encoded F1 fimbriae, curli, E. coli common pilus and enterobactin. All but one isolate carried ibeB gene, while ibeA was present mostly in B2 and F phylogenetic groups, but not in isolates from other groups. A siderophore system salmochelin (81%), haemolysin F (90.6%) and serum-resistance associated proteins, Iss (87.5%) and TraT (78%) were present in most isolates with generally equal distribution in all phylogenetic groups. Full SitABCD iron transport system was detected in 78% isolates, outer membrane protease (OmpT) and colicin V synthesis protein (CvaC) in 68.8% and 59% isolates, respectively.

Several virulence- genes were associated with particular branches of the phylogenetic tree. For example, Stg fimbriae, Ycb fimbriae, CFA/1 fimbriae and genes associated with ETT2 (E. coli type III secretion system 2) were common in B1 and C groups, but almost or entirely absent in B2 phylogenetic group. In contrast, yersiniabactin, aerobactin and hemin receptor (chuA), as well as OmpA, capsular antigens (kpsD, kpsT, kpsM), pathogenicity-island marker malX, uropathopathogenic-specific protein (usp) and afore-mentioned brain-endothelium invasin (ibeA) appeared to be B2 group-associated.

Lastly, some well known virulence-associated genes were detected uniquely amongst our isolates. For example, complete pap operon was present only in the one ST95 isolate, as well as the intimin-like adhesin (fdeC). Similarly, neuC gene was found in ST95, ST140 and two ST429 isolates (all B2 phylogroup) and K99 (F5) fimbriae only in ST354 isolate of the F phylogroup. For overview of all virulence-associated genes detected please see the Table 2, additional file 2.

2.5 Identification of plasmid replicons

All but one isolate harboured a replicon of the F incompatibility group, FIB replicon being the most commonly detected (31/32 isolates; 96.9%). At the same time, Col replicons were detected in most (22/32; 68.8%) isolates. Groups IncB/O/K/Z (10; 31.3%) and IncX1 (6; 18.8%) also appeared relatively frequently, while others were identified only in individual isolates. The IncHI1B replicon was identified in two isolates of the ST429 and ST23 cluster, respectively. Overall, types and number of replicons varied greatly even within the two closely related clusters. (For overview of replicons, please see the Table 1, additional file 1.)

The aim of current study was to evaluate diversity of colibacillosis-associated isolates from the Czech Republic. Indeed, analysis showed an immense phenotypic and genotypic variability, the isolates differing greatly in their antimicrobial-resistance phenotype, virulence genes profile and plasmid content, together having little in common. As generally acknowledged, there is no specific combination of virulence genes that would accurately define the APEC pathotype (16). The most prevalent APEC genes are also frequently present in commensal strains. There is an abundance of adhesins and iron-transporting systems, which may be considered essential prerequsites of extraintestinal pathogenicity in all types of avian and mammalian disease, but also fitness factors enabling asymptomatic colonization of healthy hosts and effective transmission. Presence of Col-V-associated genes such as iroN, iss, iutA, ompT etc. is characteristic for most APEC, more than UPEC and NMEC (15), nevertheless, their exact role in pathogenesis remains unclear or controversial (30,31). Col-V-like plasmids are, however, acknowledged as markers of poultry-adapted pathogenic strains (5,22).

As expected, the phylogenetic analysis also revealed a substantial diversity of isolates, that originated from all phylogenetic groups with the exception of group E, the most prevalent was B2 phylogroup, which is, along with D, considered typical group for human ExPEC (14). However, the second most prevalent phylogroup was B1 (7 isolates), a group commonly associated with intestinal pathogenic or commensal fecal strains (32). Interestingly, there were notable differences in virulence trait distribution among phylogenetic groups, although the isolates had been collected from the same types of infection. The idea of pathogenic strains with quite a different combination of virulence genes with alternative functions causing the same clinical disease has been proposed by Mokady et al. (33) and points out the importance of horizontal gene transfer enabling rapid adaptation to new niches by expression of certain genes in a different genetic background (34). Notably, it was the presence of typical Col-V plasmid-associated genes such as ompT, iss, cvaC, iro and sit (but suprisingly not iut, iuc for aerobactin) that were equally distributed among isolates from all phylogenetic groups.

Despite the overall diversity, the phylogenetic analysis revealed two clusters (ST429, group B2, and ST23, groups C), both containing four similar isolates that were obtained from different farms in Northern Moravia. Two ST23 isolates identical according to the core genome analysis were collected at the same day on two different farms, indicating a possible clonal spread in the locality. Both were isolates of yolk sac infection of one-day-old chicken, however, coming from different hatcheries. Colibacillosis outbreaks caused by a specific pathogenic clone have been repeatedly reported (e.g. 12,35). On the other hand, a closely related isolate (25 SNPs difference) had been collected on an unrelated farm approximately half a year before. Similar situation was observed in the ST429 cluster – the most similar isolates were from the same date and were separated from the other isolates of this cluster (with 26 – 61 SNPs difference) by a span of several months. One may speculate these isolates could have a common origin, however, the question, whether these clones may become established somehow in the production chain and circulate between flocks or farms for a long time period or whether a repeated introduction occurs from a specific source, remains unanswered. The isolates of the ST429 cluster were obtained from one-day-old chicken coming from four different hatcheries. Although an evidence for „pseudo-vertical“ spread through the production pyramid has been proposed recently (36), this fact suggests that the hatchery is probably not the source. The problem of possible reservoir of pathogenic strains for Northern Moravian farms should be addressed more closely in the future.

Both ST429 and ST23 are considered as predominant APEC lineages that are frequently isolated from poultry with clinically manifested disease (35,37), but also poultry products (26). Although representing quite unrelated APEC clades, they both appear to be poultry-specific, with little pathogenic potential for humans (7,17). In fact, an APEC strain χ7122 (ST23) has been shown to be phylogenetically closer to human ETEC (without any enterotoxin production) than to ExPEC (38). Therefore, in our collection, one may consider the two clusters, ST429 and ST23, representatives of phylogenetically distant lineages presumably associated with the same disease, again underlining the importance of accessory genome in virulence potential of APEC. The ST429 isolates had a slightly higher average number of virulence-asociated genes than ST23 isolates (172 vs. 154) including genes encoding capsule production (kpsM, T, D, neuC), invasins (ibeA, ompA) and iron-binding systems (aerobactin, yersiniabactin, chu) that the ST23 cluster (not all ST23 isolates) lacked. In contrast, ST23 isolates were characterized by presence of Stg fimbriae and ETT2-related genes. This transport system, even in degenerate state, has been reported to enhance virulence in APEC (39). Both sequence types coded for curli, F1 fimbriae, salmochelin, OmpT, TraT, Iss, however, only Iro, OmpT nad Iss have been reported to occur in significantly higher prevalence in APEC than avian-faecal E. coli (AFEC) (4). Nevertheless, it probably supports the idea of feasibility and usefulness of PCR typing targeting such potential markers of APEC derived from distant phylogenetic groups (e.g. 3).

Two isolates were assigned to ST117 (phylogenetic group F). Recent studies indicate that this sequence type comprises important APEC lineages that are repeatedly reported from colibacillosis outbreaks in different countries (25,37,40-42), but are also highlighted as potential zoonotic pathogens for containing ExPEC-related virulence genes and being isolated from both retail poultry meat and human clinical urinary tract infections (43). The remaining phylogroup F isolate was ST354, another potentially zoonotic ST, reported particularly from human and animal healthcare facilities and characterized by common resistance to antimicrobials including fluoroquinolones (44,45). This isolate carried bla_CMY-2,-59and encoded multiple adhesins including K99/F5 fimbriae, which were not found anywhere else. Both ST117 and ST354 were highly prevalent among ESBL/AmpC positive chicken isolates and it has been proposed that these lineages exhibit particularly effective host colonization and persistence in the environment (41,45).

ST95 is probably the most important pandemic ExPEC lineage that is frequently isolated from chickens (23,25,46). In fact, it may represent, along with closely related ST140, that part of B2 phylogroup where human ExPEC and APEC form a single „subpathotype“ of genetically indistingushable strains (17-19). In humans, ST95 was associated with bloodstream infections, UTIs and meningitis, often characterized by serogroups O1, O2, O45, flagellar antigen H7 and K1 capsule (typical feature of NMEC) and, in contrast to other pandemic lineages, relatively low tendency to acquire antimicrobial resistence (20,47). Indeed, not every ST95 seems to be zoonotic, as was shown with APEC O1 in a murine infectious model (48). On the other hand, our ST95 isolate fulfilled the molecular criteria for UPEC as defined by Johnson et al. (49).

Antimicrobial-resistance profile ranged from full susceptibility to all antimicrobials tested to multidrug resistance, with dominating resistance to β-lactams (ampicillin) and first generation quinolones (nalidixic acid). Resistance to β-lactams was associated largely with TEM-type β-lactamase production. No selection procedure to obtain ESBL/AmpC producing isolates had been used and we did not detect any bla_CTX-Mgene, while four isolates carried bla_CMY-2gene. This gene, along with bla_CTX-M-1, is the most common ESBL/AmpC β-lactamase in poultry E. coli isolates (50). While in most quinolone-resistant isolates a chromosomal mutation in gyrA gene was detected, seven both susceptible or resistant isolates carried qnrS-1. Co-occurence of bla-_{CMY-2, -59} and qnrS-1 was observed in two isolates from the ST429 cluster and all but one ST23 isolates carried the remaining qnrS-1 genes. One may assume that the afore-mentioned fact that these STs are not commonly associated with human disease does not make them epidemiologically irrelevant, for they may still serve as a source of resistance or virulence determinants in horizontal gene transfer. Indeed, the importance of horizontal gene transfer may be assumed from the detection of a multitude of replicons previously associated with both resistance and virulence gene spread (51-55).

Despite its limitations due to relatively small number of isolates of completely sequenced isolates, this study could be considered a basic overview of the diversity of colibacillosis-associated E. coli isolates and a delineation of paths that are to be followed in more extensive monitoring of virulent clones occuring in central Europe, as well as more elaborate analysis of their phylogenetic background and accessory genome and the role they play in adaptation of different APEC lineages to different hosts, infection types and routes of transmission. Genomic analysis of a collection of poultry colibacillosis-associated isolates revealed two clusters of phylogenetically distant lineages (ST429 and ST23) alongside a great diversity of other sequence types. In general, the collection showed a split into isolates from phylogroups F, B2 and D on one side and A, C and B1 on the other, distinctly differing in distribution of several virulence-associated genes. Clearly more research is needed to assess whether they differ also in their virulence potential and other features.

5.1 Strains isolation and preliminary characterisation

Samples have been collected since 2014 at various Czech, Slovakian and Romanian farms with increased mortality due to colibacillosis: mostly internal organ swabs from one-day chicken with yolk sac infection and septicemia or from broilers and layers suffering from colisepticemia and polyserositis (peritonitis, perihepatitis, pericarditis, airsaculitis or haemorrhagic septicemia). The origin of isolates is shown in the Table 4 (additional file 2). The chicken, from which the sequenced isolates were obtained, originated from 4 different hatcheries in the Czech Republic, Slovakia and Hungary. For 3 isolates, the hatchery has not been traced back. For the codes indicating the farm and hatchery of origin please see Table 4 (additional file 2). The samples were cultivated on McConkey agar (37°C for 18 hours aerobically) (Oxoid, UK), subcultivated on Columbia blood agar (the same conditions) (Oxoid, UK) and identified by MALDI-TOF MS (Bruker-Daltonics, Germany). Antibiotic susceptibility to selected antimicrobials – ampicillin (10 μg), amoxicillin-clavulanic acid (20/10 μg), cephalotin (30 μg), sulphonamide compounds (250 – 300 μg), gentamicin (10 μg), nalidixic acid (30 μg), sulphamethoxasol-trimethoprim (1,25/23,75 μg), tetracycline (30 μg), chloramphenicol (30 μg), ciprofloxacin (5 μg) – was tested by disc diffusion method (Oxoid, UK) and interpreted according to the Clinical and Laboratory Standard Institute (56).

Within preliminary characterization of 95 isolates, a slide agglutination test with four commercial antisera (O1, O8, O18 and O78) was performed according to the manufacturer‘s instructions (Denka Seiken, Japan). Presence of several selected resistance and virulence-associated genes (additional file 3 – table 3) was detected by PCR. After this preliminary characterization 32 strains were selected for whole-genome sequencing. To encompass the greatest possible diversity, we exluded isolates from the same individual, the same farm or the same date of isolation, if they showed the identical resistance phenotype and gene profile.

5.2 DNA extraction and whole-genome sequencing

NucleoSpin Tissue DNA extraction kit (Macherey-Nagel, Germany) following manufacturer’s instructions was used to obtain pure DNA. The DNA libraries were prepared with Nextera XT Library preparation kit (Illumina, USA). Finally, Illumina Next-Seq and Mi-Seq platforms were used for the whole-genome sequencing to obtain 2 x 150-bp or 2 x 300-bp paired-end reads, respectively.

5.3 Data processing

Adaptor residues and low quality (Q ≤ 20) ends were removed from the reads using Trimmomatic v0.36. (57). De novo assembly was performed using SPAdes assembler v3.12.0 (58). Contigs were submitted to online typing tools (Centre for Genomic Epidemiology, Technical University of Denmark; http://www.genomicepidemiology.org/): ResFinder 3.1 (59), PlasmidFinder 2.0 (60), SeroTypeFinder 2.0 (61), MLST 2.0 (62). One isolate was assigned as a novel ST8874 by Enterobase. Presence of resistance and virulence genes was predicted using CARD (Comprehensive Antibiotic Resistance Database) and VFDB (Virulence Factor Database) (63), respectively, in the ABRicate v0.8.13 programme (https://card.mcmaster.ca/, https://github.com/tseemann/abricate). The threshold for gene identity was set to 95%. The Clermont typing tool was used to classify isolates into phylogenetic groups (64)(http://clermontyping.iame-research.center/). In order to investigate the genetic relationships between isolates, genomes were annotated using Prokka v1.13 (65) and the core genome alignment was performed using Roary (66). The core genome alignment was used to determine the single nucleotide polymorphism (SNP) distance using snp-dist (https://github.com/tseemann/snp-dists). Phylogenetic tree was constructed using RAxML v8.2.10 using GTR+GAMMA+I model (67). Phylogenetic tree was then visualised via iTOL (68)(https://itol.embl.de/). The raw sequencing data were deposited to GenBank under BioProject PRJNA553636 and corresponding accession numbers to SRA for each sample can be found in Table 4 (additional file 4).

3.1 Ethic approval and consent to participate

No animal care committee approval was necessary for the purposes of this study, as no animal work was required. The samples were collected from diseased or dead animals by practical veterinarians in cooperation with the farm owners and with their consent. The authors declare that all applicable ethical guidelines were followed.

3.2 Consent for publication

Not applicable.

3.3 Availability of data and materials

The genome sequences of the 32 isolates are available in the GenBank under the bioproject ID PRJNA553636. The corresponding accession numbers are stated in the Additional file 4. https://www.ncbi.nlm.nih.gov/bioproject/PRJNA553636/.

3.4 Competing interests

The authors declare they have no conflict of interest.

3.5 Funding

The project was funded by Internal Grant Agenture of the University of Veterinary and Pharmaceutical Sciences, Brno, project 107/2017/FVL, and project NAZV QK1810462. The funding body did not have any role in the study design and collection, analysis, and interpretation of data and in writing the manuscript.

3.6 Authors’ contributions

DS collected the samples and provided anamnestic information. MM performed the laboratory analysis and interpreted the data. DC and EJ performed the whole-genome sequencing. AV processed and interpreted the genomic data. AP performed the laboratory analysis and data interpretation and is the major contributor in writing. AC provided supervision and consultaion and was the major author of the project. All authors have read and approved the manuscript.

3.7 Acknowledgements

Field vets are thanked for participation in the project, for collecting samples and data. We would like to thank to Mrs. Jana Hofirkova, Mrs. Jarmila Kovarova and MVDr. Ludmila Kohoutova for their excellent help during precessing of the samples in the laboratory. We also thank to Iva Jamborova for her assistance with DNA sequencing.

AMR – antimicrobial resistance; AFEC – avian fecal Escherichia coli; APEC – avian pathogenic Escherichia coli; ESBL – extended-spectrum beta-lactamases; ETEC – enterotoxinogenic Escherichia coli; ExPEC – extraintestinal pathogenic Escherichia coli; MLST – multilocus sequence typing; NMEC – neonatal meningitis Escherichia coli; PCR – polymerase chain reaction; PMQR – plasmid-mediated quinolone resistance; SePEC – septicemia-associated pathogenic Escherichia coli; SHS – swollen head syndrome; SNP – single nucleotide polymorphism; ST – sequence type; UPEC – uropathogenic Escherichia coli; UTI – urinary tract infection; WHS – whole genome sequencing.

Nolan L K, Barnes HJ, Vaillancourt JP, Abdul-Aziz T, Logue CM. Colibacillosis. In: Swayne D. E., Glisson J. R., et al. (eds). Diseases of Poultry. Ames, Iowa, USA: John Wiley & Sons, Inc., 2013. p. 751-805.
La Ragione RM, Woodward MJ. Virulence factors of Escherichia coli serotypes associated with avian colisepticaemia. Res Vet Sci. 2002; 73: 27 – 35.
Johnson TJ, Wannemuehler Y, Doetkott C, Johnson SJ, Rosenberger SC, Nolan LK. Identification of Minimal Predictors of Avian Pathogenic Escherichiacoli Virulence for Use as a Rapid Diagnostic Tool. J Clin Microbiol. 2008; 46: 3987–3996.
Rodriguez-Siek KE, Giddings CW, Doetkott C, Johnson TJ, Nolan LK. Characterizing the APEC pathotype. Vet Res. 2005; 36: 241–256.
Stromberg ZR, Johnson JR, Fairbrother JM, Kilbourne J, Van Goor A, Curtiss R 3^rd, Mellata M. Evaluation of Escherichia coli isolates from healthy chickens to determine their potential risk to poultry and human health. PLoS ONE. 2017; 12: e0180599.
Collingwood C, Kemmett K, Williams N, Wigley P. Is the Concept of Avian Pathogenic Escherichia coli as a Single Pathotype Fundamentally Flawed? Front Vet Sci. 2014; 1, 5: 1 – 4.
Ewers C, Antão EM, Diehl I, Philipp HC, Wieler LH. Intestine and Environment of the Chicken as Reservoirs for Extraintestinal Pathogenic Escherichia coli Strains with Zoonotic Potential. Appl Environ Microbiol. 2009; 75: 184–192.
Kemmett K, Humphrey T, Rushton S, Close A, Wigley P, Williams NJ. A Longitudinal Study Simultaneously Exploring the Carriage of APEC Virulence Associated Genes and the Molecular Epidemiology of Faecal and Systemic coli in Commercial Broiler Chickens. PLoS ONE. 2013; 8: e67749.
Maturana VG, de Pace F, Carlos C, Mistretta Pires M, Amabile de Campos T, Nakazato G, Guedes Stheling E, Logue CM, Nolan LK, Dias da Silveira. Subpathotypes of Avian Pathogenic Escherichia coli (APEC) Exist as Defined by their Syndromes and Virulence Traits. Open Microbiol J. 2011; 5: 55–64.
Jordan FTW, Williams NJ, Wattret A, Jones T. Observations on salpingitis, peritonitis and salpingoperitonitis in a layer breeder flock. Vet Rec. 2005; 157: 573–577.
Olsen RH, Thøfner ICN, Pors SE, Pires dos Santos T, Christensen JP. Experimental induced avian E. coli salpingitis: Significant impact of strain and host factors on the clinical and pathological outcome. Vet Microbiol. 2016; 188: 59–66.
Pires-dos-Santos T, Bisgaard M, Christensen H. Genetic diversity and virulence profiles of Escherichia coli causing salpingitis and peritonitis in broiler breeders. Vet Microbiol. 2013; 162: 873–880.
Mellata M. Human and Avian Extraintestinal Pathogenic Escherichia coli: Infections, Zoonotic Risks, and Antibiotic Resistance Trends. Foodborne Pathog Dis. 2013; 10: 916–932.
Johnson JR, Russo TA. Extraintestinal pathogenic Escherichia coli : “The other bad E coli .” J. Lab Clin Med. 2002; 139: 155–162.
Ewers C, Antão EM, Diehl I, Philipp HC, Wieler LH. Intestine and Environment of the Chicken as Reservoirs for Extraintestinal Pathogenic Escherichia coli Strains with Zoonotic Potential. Appl Environ Microbiol. 2009; 75: 184–192.
Guabiraba R, Schouler C. Avian colibacillosis: still many black holes. FEMS Microbiol Lett. 2015; 362: fnv118.
Zhu Ge X, Jiang J, Pan Z, Hu L, Wang S, Wang H, Leung FC, Dai J, Fan H. Comparative genomic analysis shows that avian pathogenic Escherichia coli isolate IMT5155 (O2:K1:H5; ST complex 95, ST140) shares close relationship with ST95 APEC O1:K1 and human ExPEC O18:K1 strains. PLoS ONE. 2014; 9: e112048–e112048.
Moulin-Schouleur M, Répérant M, Laurent S, Brée A, Mignon-Grasteau S, Germon P, Rasschaert D, Schouler C. Extraintestinal Pathogenic Escherichia coli Strains of Avian and Human Origin: Link between Phylogenetic Relationships and Common Virulence Patterns. J Clin Microbiol. 2007; 45: 3366–3376.
Tivendale KA, Logue CM, Kariyawasam S, Jordan D, Hussein A, Li G, Wannemuehler Y, Nolan LK. Avian-pathogenic Escherichia coli strains are similar to neonatal meningitis coli strains and are able to cause meningitis in the rat model of human disease. Infect Immun. 2010; 78: 3412–3419.
Riley LW. Pandemic lineages of extraintestinal pathogenic Escherichia coli. Clin Microbiol Infect. 2014; 20: 380–390.
Cunha MPV, Saidenberg AB., Moreno AM, Ferreira AJP, Vieira MAM., Gomes TAT, Knöbl T. Pandemic extra-intestinal pathogenic Escherichia coli (ExPEC) clonal group O6-B2-ST73 as a cause of avian colibacillosis in Brazil. PLoS One. 2017; 12: e0178970–e0178970.
Liu CM, Stegger M, Aziz M, Johnson TJ, Waits K, Nordstrom L, Gauld L, Weaver B, Rolland D, Statham S, Horwinski J, Sariya S, Davis GS, Sokurenko E, Keim P, Johnson JR, Price LB. Escherichia coli ST131-H22 as a Foodborne Uropathogen. mBio. 2018; 9, 4: e00470-18.
Maluta RP, Logue CM, Casas MRT, Meng T, Guastalli EAL, Rojas TCG, Montelli AC, Sadatsune T, de Carvalho Ramos M, Nolan LK, da Silveira WD. Overlapped sequence types (STs) and serogroups of avian pathogenic (APEC) and human extra-intestinal pathogenic (ExPEC) Escherichia coli isolated in Brazil. PLoS ONE. 2014; 9: e105016.
Manges AR. Escherichia coli and urinary tract infections: the role of poultry-meat. J Clin Microbiol Infect. 2016; 22: 122–129.
Mora A, Viso S, López, C, Alonso MP, García-Garrote F, Dabhi G, Mamani R, Herrera A, Marzoa J, Blanco M, Blanco JE, Moulin-Schouleur M, Schouler C, Blanco J. Poultry as reservoir for extraintestinal pathogenic Escherichia coli O45:K1:H7-B2-ST95 in humans. Vet Microbiol. 2013; 167: 506–512.
Pietsch M, Irrgang A, Roschanski N, Brenner Michael G, Hamprecht A, Rieber H, Käsbohrer A, Schwarz S, Rösler U, Kreienbrock L, Pfeifer Y, Fuchs S, Werner G, RESET Study Group. Whole genome analyses of CMY-2-producing Escherichia coli isolates from humans, animals and food in Germany. BMC Genomics. 2018; 19: 601.
Cordoni G, Woodward MJ, Wu H, Alanazi M, Wallis T, La Ragione RM. Comparative genomics of European avian pathogenic coli (APEC). BMC Genomics. 2016; 17: 960.
Paudel S, Stessl B, Hess C, Zloch A, Hess M. High genetic diversity among extraintestinal Escherichia coli isolates in pullets and layers revealed by a longitudinal study. BMC Vet Res. 2016;12(1):221.
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.
Johnson TJ, Siek KE, Johnson SJ, Nolan LK. DNA Sequence of a ColV Plasmid and Prevalence of Selected Plasmid-Encoded Virulence Genes among Avian Escherichia coli J Bacteriol. 2006; 188: 745.
Mellata M, Dho-Moulin M, Dozois CM, Curtiss R 3rd, Brown PK, Arné P, Brée A, Desautels C, Fairbrother JM. Role of virulence factors in resistance of avian pathogenic Escherichia coli to serum and in pathogenicity. Infect Immun. 2003; 71: 536–540.
Tenaillon O, Skurnik D, Picard B, Denamur E. The population genetics of commensal Escherichia coli. Nat Rev Microbiol. 2010; 8: 207–217.
Mokady D, Gophna U, Ron EZ. Extensive gene diversity in septicemic Escherichia coli strains. J Clin Microbiol. 2005; 43: 66–73.
Escobar-Parámo P, Blanc-Potard AB, Bui H, Le Bouguénec C, Denamur E. A Specific Genetic Background Is Required for Acquisition and Expression of Virulence Factors in Escherichia coli. Mol Biol Evol. 2004; 21: 1085–1094.
Solà-Ginés M, Cameron-Veas K, Badiola I, Dolz R, Majó N, Dahbi G, Viso S, Mora A, Blanco J, Piedra-Carrasco N, González-López JJ, Migura-Garcia L. 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.
Projahn M, Daehre K, Roesler U, Friese A. Extended-Spectrum-Beta-Lactamase- and Plasmid-Encoded Cephamycinase-Producing Enterobacteria in the Broiler Hatchery as a Potential Mode of Pseudo-Vertical Transmission. J Appl Environ Microbiol. 2016; 83: e02364-16.
Cummins ML, Reid CJ, Roy Chowdhury P, Bushell RN, Esbert N, Tivendale KA, Noormohammadi AH, Islam S, Marenda MS, Browning GF, Markham PF, Djordjevic SP. Whole genome sequence analysis of Australian avian pathogenic Escherichia coli that carry the class 1 integrase gene. Microb Genom. 2019; 5: e000250.
Dziva F, Hauser H, Connor TR, van Diemen PM, Prescott G, Langridge GC, Eckert S, Chaudhuri RR, Ewers C, Mellata M, Mukhopadhyay S, Curtiss R 3rd, Dougan G, Wieler LH, Thomson NR, Pickard DJ, Stevens MP. Sequencing and functional annotation of avian pathogenic Escherichia coli serogroup O78 strains reveal the evolution of coli lineages pathogenic for poultry via distinct mechanisms. Infect Immun. 2013; 81: 838–849.
Ideses D, Gophna U, Paitan Y, Chaudhuri RR, Pallen MJ, Ron EZ. A Degenerate Type III Secretion System from Septicemic Escherichia coli Contributes to Pathogenesis. J Bacteriol. 2005; 187: 8164 - 8171.
Dissanayake DRA, Octavia S, Lan R. Population structure and virulence content of avian pathogenic Escherichia coli isolated from outbreaks in Sri Lanka. Vet Microbiol. 2014; 168: 403–412.
Projahn M, Daehre K, Semmler T, Guenther S, Roesler U, Friese A. Environmental adaptation and vertical dissemination of ESBL-/pAmpC-producing Escherichia coli in an integrated broiler production chain in the absence of an antibiotic treatment. Microb Biotechnol. 2018; 11: 1017–1026.
Ronco T, Stegger M, Olsen RH, Sekse C, Nordstoga AB, Pohjanvirta T, Lilje B, Lyhs U, Andersen PS, Pedersen K. Spread of avian pathogenic Escherichia coli ST117 O78:H4 in Nordic broiler production. BMC Genomics 2017; 18: 13.
Vincent C, Boerlin P, Daignault D, Dozois CM, Dutil L, Galanakis C, Reid-Smith RJ, Tellier PP, Tellis PA, Ziebell K, Manges AR. Food reservoir for Escherichia coli causing urinary tract infections. Emerg Infect Dis. 2010; 16: 88–95.
Guo S, Wakeham D, Brouwers HJM, Cobbold RN, Abraham S, Mollinger JL, Johnson JR, Chapman TA, Gordon DM, Barrs VR, Trott DJ. Human-associated fluoroquinolone-resistant Escherichia coli clonal lineages, including ST354, isolated from canine feces and extraintestinal infections in Australia. Microbes Infect. 2015; 17: 266–274.
Vangchhia B, Abraham S, Bell JM, Collignon P, Gibson JS, Ingram PR, Johnson JR, Kennedy K, Trott DJ, Turnidge JD, Gordon DM. Phylogenetic diversity, antimicrobial susceptibility and virulence characteristics of phylogroup F Escherichia coli in Australia. Microbiol. 2016; 162: 1904–1912.
Jørgensen SL, Stegger M, Kudirkiene E, Lilje B, Poulsen LL, Ronco T, Pires Dos Santos T, Kiil K, Bisgaard M, Pedersen K, Nolan LK, Price LB, Olsen RH, Andersen PS, Christensen H. Diversity and Population Overlap between Avian and Human Escherichia coli Belonging to Sequence Type 95. 2019; mSphere 4, e00333-18.
Mora A, López C, Dabhi G, Blanco M, Blanco JE, Alonso MP, Herrera A, Mamani R, Bonacorsi S, Moulin-Schouleur M, Blanco J. Extraintestinal pathogenic Escherichia coli O1:K1:H7/NM from human and avian origin: detection of clonal groups B2 ST95 and D ST59 with different host distribution. BMC Microbiol. 2009; 9: 132.
Johnson TJ, Wannemuehler Y, Kariyawasam S, Johnson JR, Logue CM, Nolan LK. Prevalence of Avian-Pathogenic Escherichia coli Strain O1 Genomic Islands among Extraintestinal and Commensal coli Isolates. J Bacteriol. 2012; 194: 2846–2853.
Johnson JR, Murray AC, Gajewski A, Sullivan M, Snippes P, Kuskowski MA, Smith KE. Isolation and Molecular Characterization of Nalidixic Acid-Resistant Extraintestinal Pathogenic Escherichia coli from Retail Chicken Products. Antimicrob Agents Chemother. 2003; 47: 2161–2168.
van Hoek AHAM, Stalenhoef JE, van Duijkeren E, Franz E. Comparative virulotyping of extended-spectrum cephalosporin-resistant coli isolated from broilers, humans on broiler farms and in the general population and UTI patients. Vet Microbiol. 2016; 194: 55–61.
Dolejska M, Villa L, Poirel L, Nordmann P, Carattoli A. Complete sequencing of an IncHI1 plasmid encoding the carbapenemase NDM-1, the ArmA 16S RNA methylase and a resistance–nodulation–cell division/multidrug efflux pump. J Antimicrob Chemother. 2013; 68: 34–39.
Dolejska M, Duskova E, Rybarikova J, Janoszowska D, Roubalova E, Dibdakova K, Maceckova G, Kohoutova L, Literak I, Smola J, Cizek A. Plasmids carrying blaCTX-M-1 and qnr genes in Escherichia coli isolates from an equine clinic and a horseback riding centre. J Antimicrob Chemother. 2011; 66: 757–764.
Johnson TJ, Nolan LK. Pathogenomics of the Virulence Plasmids of Escherichiacoli. Microbiol. Mol Biol Rev. 2009; 73: 750–774.
Kao CY, Chen JW, Liu TL, Yan JJ, Wu JJ. Comparative Genomics of Escherichia coli Sequence Type 219 Clones From the Same Patient: Evolution of the IncI1 bla_CMY-Carrying Plasmid in Vivo. Front Microbiol. 2018; 9: 1518.
Royden A, Ormandy E, Pinchbeck G, Pascoe B, Hitchings MD, Sheppard SK, Williams NJ. Prevalence of faecal carriage of extended-spectrum β-lactamase (ESBL)-producing Escherichia coli in veterinary hospital staff and students. Vet Rec Open. 2019; 6: e000307.
Clinical and Laboratory Standards Institute (CLSI). Performance Standards for Antimicrobial Susceptibility Testing; Twenty-fifth Informational Supplement. CLSI document M100-S25. Wayne, PA, USA, Clinical and Laboratory Standards Institute; 2015.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 2014; 30: 2114–2120.
Bankevich A, Nurk, S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012; 19: 455–477.
Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, Aarestrup FM, Larsen MV. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012; 67: 2640–2644.
Carattoli A, Zankari E, García-Fernández A, Voldby Larsen M, Lund O, Villa L, Møller Aarestrup F, Hasman H. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob. Agents Chemother. 2014; 58: 3895–3903.
Joensen K G, Tetzschner AM, Iguchi A, Aarestrup FM, Scheutz F. Rapid and easy in silico serotyping of Escherichia coliusing whole genome sequencing (WGS) data. J Clin Microbiol. 2015; 53(8): 2410-2426.
Larsen MV, Cosentino S, Rasmussen S, Friis C, Hasman H, Marvig RL, Jelsbak L, Sicheritz-Pontén T, Ussery DW, Aarestrup FM, Lund O. Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol. 2012; 50: 1355–1361.
Liu B, Zheng D, Jin Q, Chen L, Yang J. VFDB 2019: a comparative pathogenomic platform with an interactive web interface. Nucleic Acids Res. 2019; 47(D1): D687–D692.
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: e000192.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014; 30(14): 2068-2069.
Page JA, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, Fookes M, Falush D, Keane JA, Parkhill J. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015; 15, 31: 3691-3693.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014; 30(9): 1312-1313.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016; 44: W242–W245.

Download PDF

Journal Publication

published 09 Jun, 2020

Read the published version in BMC Veterinary Research →

Submission checks completed at journal
03 Jun, 2020
Editorial decision: Accept
02 Jun, 2020

You are reading this latest preprint version

Genomic analysis of Escherichia coli strains isolated from diseased chicken in the Czech Republic

Status:

Journal Publication

Version 3

Abstract

Figures

Background

Results

2.1 Samples collection and preliminary characterization

2.2 In silico serotyping, MLST and phylogenetic analysis of 32 selected isolates

2.3 Identification of resistance genes

2.4 Identification of virulence genes

2.5 Identification of plasmid replicons

Discussion

Conclusions

Methods

5.1 Strains isolation and preliminary characterisation

5.2 DNA extraction and whole-genome sequencing

5.3 Data processing

Declarations

3.1 Ethic approval and consent to participate

3.2 Consent for publication

3.3 Availability of data and materials

3.4 Competing interests

3.5 Funding

3.6 Authors’ contributions

3.7 Acknowledgements

Abbreviations

References

Supplementary Files

Status:

Journal Publication

Version 3