Front. Microbiol.Frontiers in MicrobiologyFront. Microbiol.1664-302XFrontiers Media S.A.10.3389/fmicb.2020.01699MicrobiologyOriginal ResearchGenomic and Phenotypic Analysis of an ESBL-Producing E. coli ST1159 Clonal Lineage From Wild Birds in MongoliaSchierackPeter1HeidenStefan E.2KhanMuhammad Moman1NikolausLena2KolendaRafal1StubbeMichael3LkhagvasurenDavaa4RödigerStefan15GuentherSebastian2SchauflerKatharina2*1Multiparametric Diagnostics, Brandenburg University of Technology Cottbus – Senftenberg, Senftenberg, Germany2Institute of Pharmacy, University of Greifswald, Greifswald, Germany3Institute of Biology, Martin-Luther-University Halle-Wittenberg, Halle (Saale), Germany4Department of Biology, National University of Mongolia, Ulaanbaatar, Mongolia5Faculty of Health Sciences, Joint Faculty of the Brandenburg University of Technology Cottbus – Senftenberg, The Brandenburg Medical School Theodor Fontane and the University of Potsdam, Senftenberg, Germany
Edited by: Teresa M. Coque, Ramón y Cajal Institute for Health Research, Spain
Reviewed by: Alain Hartmann, Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE), France; Ricardo León-Sampedro, Ramón y Cajal University Hospital, Spain
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Background
In addition to the broad dissemination of pathogenic extended-spectrum beta-lactamase (ESBL)-producing Escherichia (E.) coli in human and veterinary medicine and the community, their occurrence in wildlife and the environment is a growing concern. Wild birds in particular often carry clinically relevant ESBL-producing E. coli.
Objectives
We analyzed ESBL-producing and non-ESBL-producing E. coli obtained from wild birds in Mongolia to identify phylogenetic and functional characteristics that would explain the predominance of a particular E. coli clonal lineage in this area.
Methods
We investigated ESBL-producing E. coli using whole-genome sequencing and phylogenetics to describe the population structure, resistance and virulence features and performed phenotypic experiments like biofilm formation and adhesion to epithelial cells. We compared the phenotypic characteristics to non-ESBL-producing E. coli from the same background (Mongolian wild birds) and genomic results to publicly available genomes.
Results and Conclusion
We found ESBL-producing E. coli sequence type (ST) 1159 among wild birds in Mongolia. This clonal lineage carried virulence features typical for extra-intestinal pathogenic or enterotoxigenic E. coli. Comparative functional experiments suggested no burden of resistance in the ST1159 isolates, which is despite their carriage of ESBL-plasmids. Wild birds will likely disseminate these antibiotic-resistant pathogens further during migration.
ESBL – E. coliST1159Mongoliawild birdsenvironmental epidemiologyIntroduction
Extended-spectrum beta-lactamase (ESBL)-producing Escherichia (E.) coli presents a major threat to public health worldwide (Ferri et al., 2017). Pathogenic representatives cause a range of severe infectious diseases in humans and animals including sepsis, urinary tract and wound infections (Ewers et al., 2012). In addition to their increasing prevalence as infectious agents in human and veterinary medicine, their occurrence in wildlife and the environment is of great concern. It has been previously shown that wild birds carry ESBL-producing E. coli of different sequence types (ST) (Schaufler et al., 2015, 2019; Guenther et al., 2017). However, we know little about factors that influence the spread of these pathogens in the absence of antibiotic selection pressures (Nicolas-Chanoine et al., 2013). In this study, we report the broad occurrence of an ESBL-producing E. coli ST1159 clonal lineage among wild birds in Mongolia, its phylogenetic and functional characteristics and phenotypic differences to other STs and non-ESBL-producing E. coli from the same background.
Materials and MethodsOrigin and Isolation of Bacteria
Cloacal swabs of apparently healthy nestlings and juvenile wild birds were obtained in different remote areas of Mongolia during a period of 4 weeks during July 2017 as part of a “bird ringing expedition” (Figure 1). Samples of cormorants were obtained in three large bird colonies (n > 1000 individuals) during the ringing procedure, samples of birds of prey while ringing birds from individual nests and cranes were caught and sampled individually (without ringing). The sampling locations in this study are among the least densely human-populated areas in the world and large cities were not present (Guenther et al., 2012, 2017; Nicolas-Chanoine et al., 2013). With the exception of moderate ruminant farming, no agricultural activities (e.g., fertilizing fields with manure) were recorded around the sampling locations for hundreds of kilometers. We did not choose any areas with free-ranging livestock for sampling. Overall, 316 individual birds were included (241 cormorants [P. carbo], 37 cranes [A. virgo], and 38 other birds [A. monachus, M. migrans, F. tinnunculus, A. cinerea]). Cloacal swabs (Copan fecal swab, MAST) were shipped to Germany and streaked on agar plates to isolate coliform-appearing bacteria. Figure 1 and Supplementary Table 1 provide details on geographic sampling locations and host species of obtained bacteria. E. coli (both ESBL- and non-ESBL-producers) were isolated from selective chromogenic (MAST CHROMagar) plates with and without 4 μg/mL cefotaxime (the latter for ESBL-producers). One single colony per plate was selected and ESBL-producers were differentiated from AmpC-producers using the Vitek V2 System (Biomerieux, France). Overall, we obtained 49 ESBL-producing E. coli isolates. In addition, we analyzed 42 non-ESBL-producing E. coli from colony 2 (Airag Nur) to phenotypically compare the isolates.
Map of Mongolia depicting the sample locations of the various E. coli isolates. The circle area is proportional to the number of sequence types recovered (ST1159: n = 36 and n = 1; ST38: n = 2 and n = 1; ST711: n = 2; all other STs: n = 1). The new sequence type (ST-) is a single-locus variant (SLV) of ST2179 (1 single-nucleotide variant [SNV] in fumC). Map tiles by Stamen Design, under CC BY 3.0. Data by OpenStreetMap, under ODbL.
Bioinformatics Analysis
All ESBL-producing E. coli (n = 49; [group 1: all ST1159 isolates and group 2: all non-ST1159 ESBL-isolates]) were whole-genome sequenced on an Illumina MiSeq (2× 300 bp PE) in collaboration with Eurofins Genomics (Ebersberg, Germany). Raw reads were quality-trimmed, adapter-trimmed and contaminant-filtered using BBDuk from BBTools v. 38.411. After de novo assembly of trimmed reads into contiguous sequences (contigs) using shovill v. 1.0.4 in combination with SPAdes v. 3.13.12 (Bankevich et al., 2012), draft genomes were polished by mapping trimmed reads back to the contigs with bwa v. 0.7.17-r1188 (Li and Durbin, 2009) and calling SNPs and indels with Pilon v. 1.23 (Walker et al., 2014). Putative plasmid sequences were extracted using mlplasmids (Arredondo-Alonso et al., 2018) and compared to the NCBI database using BLASTN3. Then, we detected sequence types, antibiotic resistance/virulence genes and SNPs using mlst v.2.17.6, abricate v.0.9.8 (including the following databases: Resfinder, CARD, ARG-ANNOT, NCBI ARRGD, EcOH, PlasmidFinder, VFDB and Ecoli_VF) and snippy v.4.4.14,5,6. We inferred a core SNP phylogeny for ST1159. For this, alignments were filtered for recombination using Gubbins v.2.3.4 (Croucher et al., 2015) and core SNPs extracted using snp-sites v.2.4.1 (1675 sites; Page et al., 2016). We inferred a maximum likelihood tree with RAxML-NG v.0.9.0 (Kozlov et al., 2019). The best-scoring maximum likelihood tree was midpoint-rooted and visualized in FigTree7. By applying Roary v. 3.12.0 (Page et al., 2015), we assessed the pan genome of the isolates and visualized it together with the phylogeny and metadata using Phandango (Hadfield et al., 2018). We included twelve publicly available ST1159 genome sequences8, which were processed as mentioned above and one provided (from colleagues at the Institute of Microbiology and Epizootics at Freie Universität Berlin) ST1159 isolate: PBIO1289 (n = 13; [group 3]). Visualization of E. coli virulence factors9 was done with BRIG v.0.95-dev-0004 (Blast Ring Image Generator) (Alikhan et al., 2011). The synteny visualization of ST1159 isolate PBIO845 with UPEC strain UMN026 and ETEC strain UMNK88 was created with tblastx (query and database genetic code: 11; E-value: 1e-10; Camacho et al., 2009) and genoPlotR v. 0.8.9 (minimum bitscore: 50; Guy et al., 2010).
Construction of an ESBL-Plasmid-Cured Variant (PCV)
We constructed an ESBL-plasmid-free variant of one of the ST1159 isolates (PBIO845) by growing isolates at 42°C in BHI broth with daily subcultivations for 4 weeks (Schaufler et al., 2016). To confirm loss of the large ESBL-plasmid, we performed plasmid profile analysis and PCR as described previously (Schaufler et al., 2013).
Phenotypic Experiments
We carried out phenotypic experiments for the 49 ESBL- and 42 non-ESBL-producers and PCV to investigate potential phenotypic differences between and within bacterial groups. Experiments included adhesion (Frommel et al., 2013), long-term colony (Schaufler et al., 2016), and acidity tolerance assays. To test adhesion to LoVo and CHIC-8E11 epithelial cells, we used a 96-well-plate-based assay, measured adherence with a fluorescence microscope (VideoScan system) (Rodiger et al., 2013) and calculated adhesion ratios accordingly. All cells were grown and adhesion assays were carried out as previously described (Frommel et al., 2013). E. coli were grown overnight to an OD600 of 0.8–1.2. Cells were inoculated with an infection dose of about 62,500 bacteria per mm2 of a monolayer using a conversion of 3×108 bacteria/mL/OD600. After 4 h of incubation and washing with 1×PBS, cells and adherent bacteria were fixed with 4% paraformaldehyde. Cells and bacteria were stained with propidium iodide (10 μg/mL in ddH2O) and analyzed. A minimum cell monolayer area of 0.3 mm2 per well was investigated. All 150,000 images were visually analyzed for E. coli pattern formation. Isolates whose images resembled those of the negative control (eukaryotic cells without bacteria) were defined as non-adherent.
Then, we tested the isolates on long-term colony plates to assess their ability to form biofilm-related structures. Three microliters of overnight cultures grown at 37°C were dropped on span agar plates (H. Carroux, Germany) with sodium chloride (5%) and a 0.5% Congo red–0.25% Coomassie brilliant blue solution. Plates were incubated for 5 days at 28°C and screened for curli fibers and cellulose production. A purple long-term colony color indicates the production of curli fibers; elevated, structured surfaces indicate the production of cellulose. Plain white colonies indicate the absence of both (Bokranz et al., 2005). Acidity tolerance was tested in 96-well-plates with different pH values (highest at 7, lowest at 3.5 with 0.5 pH steps). Tolerance to the respective pH value was defined as growth of the isolate in the respective well after 24 h of incubation at 37°C.
All phenotypic experiments were performed in a minimum of three biological replicates and statistically evaluated using the Mann-Whitney U-test and Fisher’s Exact Test with Bonferroni correction to correct for the multiple comparisons and account for non-normally distributed values (p value < 0.0125).
Results and Discussion
Over two thirds (37/49) of ESBL-producing E. coli belonged to ST1159 (group 1; phylogroup B2). Thirty-five ST1159 were obtained from one bird colony and one bird species (Great Cormorant, 46% of the sampled nestlings). The other two ST1159 isolates originated from two species (Black Kite and Gray Heron) suggesting inter-host exchange, however, it is also possible that the isolates occurred independently of each other. The Gray Heron was sampled in the same colony as the cormorants, the Black Kite sample originated from nestlings sampled 300 km away (Figure 1). Other sequence types among the ESBL-producers included ST38 (n = 3; phylogroup D), ST711 (n = 2; phylogroup B1) and others (group 2) (Supplementary Table 1). The new sequence type (ST-) is a single-locus variant (SLV) of ST2179 caused by one single-nucleotide variant (SNV) in fumC allele 65 (C > T; Ala > Val). Previously, we have found ST38 in Mongolian wild birds that contained chromosomally encoded ESBL-genes (Guenther et al., 2017), however, in this study bioinformatics analysis suggested plasmid-encoded bla-genes in ST38.
Figure 2 reveals that all 37 ST1159 from Mongolian wild birds belong to a single clonal lineage (max. 1.4 SNPs/aligned Mbp). The twelve publicly available ST1159 genomes and PBIO1289 (group 3), the latter isolated from a duck coop in Germany in 2005 (Ewers et al., 2009) were further distant, which is not surprising given the fact that they are geographically and timely dispersed and that SNPs accumulate over time. Interestingly, the closest “external” ST1159 to group 1 isolate is a clinical extra-intestinal pathogenic E. coli (ExPEC) isolated from a urinary tract infection (Salipante et al., 2015). The duck sample was identified as avian pathogenic E. coli (APEC) (Ewers et al., 2009). It is well known that ExPEC and APEC share common pathogenic features and often belong to pathogenicity-associated phylogroup B2 (Ewers et al., 2007) which does not only underline the clinical relevance of ESBL-producing ST1159 group 1 but also serves as possible explanation for the success of ST1159 among birds.
E. coli ST1159 core SNP phylogeny. The midpoint-rooted phylogenetic tree was inferred using a maximum likelihood approach with the GTR model of nucleotide substitution and a gamma model of rate heterogeneity. It is depicted with metadata where available and the output of the Roary pangenome analysis. Isolate names in bold designate isolates whose genomes were sequenced as part of this study. Names starting with ERR or SRR designate raw read accession numbers of publicly available ST1159 isolates.
Bioinformatics and plasmid profile analysis suggested that group 1 ST1159 isolates shared identical plasmid contents, namely two large plasmids (∼90–125 kb each), of which one encoded the ESBL-enzyme (CTX-M-14) and the other several virulence features. Determined plasmid incompatibility (Inc.) types were IncFIC(FII) and IncI1. Both IncF- as well as IncI-type plasmids have been previously associated with the spread of ESBL-producing E. coli (Coque et al., 2008; Branger et al., 2018; Montealegre et al., 2020). Group 3 genomes showed only small plasmid analogies with the Mongolian representatives and did not carry any ESBL-genes. Interestingly, the group 1 plasmid sequences combined chromosomally encoded features of uro-pathogenic E. coli (which belong to ExPEC) – the pap-like operon (P pili, PAI [pathogenicity island] I(APEC-O1)) – as well as fae genes (PAI: K88 [F4+] fimbriae), which have previously been described on plasmids from enterotoxigenic E. coli (Figure 3; Shepard et al., 2012).
Synteny plot of E. coli PBIO845 (ST1159) contig00020. The plasmid-derived contig was aligned with the chromosome of UPEC strain UMN026 and plasmid pUMNK88_K88 of ETEC strain UMNK88 using tblastx (E-value 1e-10). The contig combines features commonly found in UPEC strains (P pili encoded by the pap operon) with ETEC-associated features (K88 [F4] fimbriae encoded by the fae operon). Direct comparisons are colored with red hues whereas reverse comparisons are colored with blue hues. Arrows depict coding sequences (CDS) and are colored according to the legend. Arrows with a dotted outline are annotated as pseudo CDS. While the complete sequence of contig00020 is shown (middle), the chromosome (top) and plasmid sequences (bottom) were truncated in the displayed region.
We were surprised by the high number of ST1159 isolates and thus focused on factors that would explain their broad occurrence. First, we investigated the carriage of virulence genes among ST1159 and non-ST1159 genomes. We identified a set of clinically relevant virulence factors (Ewers et al., 2007) i.e., fae (fimbrial proteins), vat (vacuolating autotransporter toxin), malX (maltose transporter), usp (uropathogenic specific protein), and ibeA (invasion protein), which were unique to group 1 and absent to group 2. Some of these were present in the external ST1159 group 3 (Supplementary Figure 1). For example, all external genomes carried vat, malX and ibeA. Note, however, that malX, a phosphotransferase system enzyme, is present in most isolates of E. coli and may be overestimated regarding its importance as virulence factor. Other virulence features were present in all groups (Supplementary Figure 1). We then investigated the potential role of several of these virulence features in functional experiments. As fae genes have been previously described as epithelial adhesins (Xia et al., 2015) we compared groups 1 and 2 as well as non-ESBL-producing isolates (n = 42; [group 4]) in adhesion experiments. We observed no significant differences in bacterial adherence to gut epithelial chicken cells (CHIC) or to those of human origin (LoVo cells). Next, we investigated phenotypes potentially relevant to the success of bacteria in wild birds and the environment. Biofilm formation is an important virulence feature, and associates to increased environmental survival (Schaufler et al., 2019). Interestingly, all group 1 isolates except one formed biofilm-associated curli and cellulose components, which was significantly better than group 4 (p = 0.003). This might point toward an advantage in tolerating environmental challenges when disseminating within and across bird colonies. It is known that bird droppings contain uric acid resulting in low pH values. We performed dilution series to investigate whether ST1159 was more tolerant to acidity than others. The non-ESBL-producers (group 4) showed no significantly better growth in low pH values than ST1159 (group 1). In an earlier study, we report that ESBL-plasmids increased virulence in some E. coli strains. We thus constructed a ST1159 (PBIO845) ESBL-plasmid-cured variant and included it in the above-mentioned phenotypic experiments. We did not observe any virulence benefit in the respective ST1159 wildtype when compared to the PBIO845 mutant.
Our results do not clarify unequivocally why the ESBL-producing ST1159 clonal lineage proliferated in wild birds from Mongolia. ST1159 is a rare sequence type, previously described a dozen times only. Given that we found this ST in 76% of all obtained ESBL-producing isolates from Mongolian wild birds suggests, however, that it has developed unique strategies to adapt to this ecology rather than arising from pure evolutionary coincidence. The intra-host diversity of E. coli is generally high. In wildlife even more than in domestic animals (Escobar-Paramo et al., 2006; Tenaillon et al., 2010). The occurrence of a single ESBL-producing E. coli clonal lineage and its spread within one bird colony and to other individuals is thus surprising, especially when taking the low anthropogenic impact in this area into account. The population density of Mongolia is about two people per square kilometer (for comparison: EU ∼ 250/km2, India ∼ 450/km2, and China ∼ 150/km2), and even lower at the Airag Nur (1/km2). Within 50 kilometers of this lake, only one small permanent settlement is present (Zavkhan ∼ 50–100 inhabitants), and the next village with a medical service is 80 kilometers away (Urgamal ∼1200 inhabitants). Uliastai, which is a city with about 10.000 inhabitants is located 400 kilometers upstream of the river Zavkhan. In addition, as intensive livestock and agricultural farming was completely absent in our sampling locations, the presence of antibiotic residues seems highly unlikely. The available knowledge indicates that different factors are responsible for the success of E. coli clones. In all experiments, ST1159 isolates –despite carrying ESBL-plasmids- were not inferior to the non-ESBL-producing group and showed significantly more often a production of curli and cellulose. These biofilm-associated components might contribute to the environmental success and spread of ST1159 in bird populations. ST1159 proliferated in these wild birds, which suggests that the frequent detection of ESBL-producing isolates in wildlife is not just based on spill-over effects but the expansion of certain clonal lineages able to cope with environmental challenges while simultaneously carrying resistance plasmids.
Conclusion
We speculate that an adult bird has introduced ST1159 to the cormorant colony and/or the Airag Nur lake first. Then, it might have spread among the cormorant nestlings, nearby Heron nestlings at the same lake, and long distance to a 300 km-away bird of prey colony (Figure 1). Cormorant colonies are densely colonized and covered in bird droppings, which results in a situation comparable to industrial poultry production. ST1159’s spread took likely place among nestlings within the colony rather than being driven by repeated introductions through contaminated fish from the lake. Fish are not a natural host for E. coli bacteria and fecal pollution of the lake is unlikely. As cormorants are migrating birds using the Central Asian flyway (Prosser et al., 2011) the origin of the ST1159 clonal lineage might be somewhere in the avian wintering grounds in the southern parts of Central Asia, where anthropogenic impacts are high. Nestlings hatching in this “bio-fermenters-like colony” may further spread ST1159 during their winter migration southwards, similar to what has been shown for cormorants carrying the H5N1 influenza virus (Hu et al., 2011). Keep in mind, however, that the biased sample set with a majority of ESBL-producing isolates and ST1159 obtained from mostly one single cormorant colony does not allow to suggest that this is a generally valid finding for the different locations in Mongolia.
Data Availability Statement
The datasets generated for this study can be found in the European Nucleotide Archive (ENA) at EMBL-EBI with the accession number PRJEB36861 (https://www.ebi.ac.uk/ena/data/view/PRJEB36861).
Author Contributions
PS, KS, MS, and SG conceived and designed the experiments. PS, SG, MS, and DL collected the data and samples. PS, KS, MK, LN, RK, and SR performed laboratory analysis. KS, SG, PS, and SH analyzed the data and wrote the manuscript. All authors have read and approved the final draft of the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The data presents results obtained by the Mongolian–German Biological Expeditions since 1962 (No. 345). We would like to thank all expedition members. We thank Torsten Semmler, Robert Koch Institute Berlin, for his support during the initial sequencing process and Antina Lübke-Becker, Institute of Microbiology and Epizootics, Freie Universität Berlin, for providing one bacterial isolate.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.01699/full#supplementary-material
BLASTN comparison (E-value 1e-10, low complexity filter turned off) of E. coli genomes versus selected E. coli virulence factors. Each concentric ring (from inner to outer) belongs to an isolate as depicted in the legend (top to bottom). Group 1 (ST1159) and 2 (Non-ST1159) contain isolates recovered from Mongolian samples. Group 3 consists of publicly available ST1159 isolates (raw read accession number given) and one isolate originating from a duck coop in 2005.
Overview for ESBL-producing E. coli: sampled avian host species, sample types, sampling location coordinates, sampling dates, identified E. coli multi-locus sequence types (STs), resistance genes and plasmid incompatibility groups based on the whole-genome sequences using mlst, abricate, and snippy.
ReferencesAlikhanN. F.PettyN. K.Ben ZakourN. L.BeatsonS. A. (2011). BLAST ring image generator (BRIG): simple prokaryote genome comparisons.BMC Genomics12:402. 10.1186/1471-2164-12-40221824423Arredondo-AlonsoS.RogersM. R. C.BraatJ. C.VerschuurenT. D.TopJ.CoranderJ. (2018). mlplasmids: a user-friendly tool to predict plasmid- and chromosome-derived sequences for single species.bioRxiv[Preprint]. 10.1101/329045BankevichA.NurkS.AntipovD.GurevichA. A.DvorkinM.KulikovA. S. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing.J. Comput. Biol.19455–477. 10.1089/cmb.2012.002122506599BokranzW.WangX.TschapeH.RomlingU. (2005). Expression of cellulose and curli fimbriae by Escherichia coli isolated from the gastrointestinal tract.J. Med. Microbiol.541171–1182. 10.1099/jmm.0.46064-016278431BrangerC.LeddaA.Billard-PomaresT.DoubletB.FouteauS.BarbeV. (2018). Extended-spectrum beta-lactamase-encoding genes are spreading on a wide range of Escherichia coli plasmids existing prior to the use of third-generation cephalosporins.Microb. Genom.4:e000203. 10.1099/mgen.0.00020330080134CamachoC.CoulourisG.AvagyanV.MaN.PapadopoulosJ.BealerK. (2009). BLAST+: architecture and applications.BMC Bioinform.10:421. 10.1186/1471-2105-10-42120003500CoqueT. M.NovaisA.CarattoliA.PoirelL.PitoutJ.PeixeL. (2008). Dissemination of clonally related Escherichia coli strains expressing extended-spectrum beta-lactamase CTX-M-15.Emerg. Infect. Dis.14195–200.CroucherN. J.PageA. J.ConnorT. R.DelaneyA. J.KeaneJ. A.BentleyS. D. (2015). Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins.Nucleic Acids Res43:e15. 10.1093/nar/gku119625414349Escobar-ParamoP.Le Menac’hA.Le GallT.AmorinC.GouriouS.PicardB. (2006). Identification of forces shaping the commensal Escherichia coli genetic structure by comparing animal and human isolates.Environ. Microbiol.81975–1984. 10.1111/j.1462-2920.2006.01077.x17014496EwersC.AntaoE. M.DiehlI.PhilippH. C.WielerL. H. (2009). Intestine and environment of the chicken as reservoirs for extraintestinal pathogenic Escherichia coli strains with zoonotic potential.Appl. Environ. Microb.75184–192. 10.1128/aem.01324-0818997030EwersC.BetheA.SemmlerT.GuentherS.WielerL. H. (2012). Extended-spectrum beta-lactamase-producing and AmpC-producing Escherichia coli from livestock and companion animals, and their putative impact on public health: a global perspective.Clin. Microbiol. Infect.18646–655. 10.1111/j.1469-0691.2012.03850.x22519858EwersC.LiG.WilkingH.KiesslingS.AltK.AntaoE. M. (2007). Avian pathogenic, uropathogenic, and newborn meningitis-causing Escherichia coli: how closely related are they?Int. J. Med. Microbiol.297163–176. 10.1016/j.ijmm.2007.01.00317374506FerriM.RanucciE.RomagnoliP.GiacconeV. (2017). Antimicrobial resistance: a global emerging threat to public health systems.Crit. Rev. Food Sci.572857–2876.FrommelU.LehmannW.RodigerS.BohmA.NitschkeJ.WeinreichJ. (2013). Adhesion of human and animal Escherichia coli strains in association with their virulence-associated genes and phylogenetic origins.Appl. Environ. Microb.795814–5829. 10.1128/aem.01384-1323872574GuentherS.AschenbrennerK.StammI.BetheA.SemmlerT.StubbeA. (2012). Comparable high rates of extended-spectrum-beta-lactamase-producing Escherichia coli in birds of prey from germany and mongolia.PLoS One7:e53039. 10.1371/journal.pone.005303923300857GuentherS.SemmlerT.StubbeA.StubbeM.WielerL. H.SchauflerK. (2017). Chromosomally encoded ESBL genes in Escherichia coli of ST38 from Mongolian wild birds.J. Antimicrob. Chemother.721310–1313. 10.1093/jac/dkx00628158613GuyL.KultimaJ. R.AnderssonS. G. (2010). genoPlotR: comparative gene and genome visualization in R.Bioinformatics262334–2335. 10.1093/bioinformatics/btq41320624783HadfieldJ.CroucherN. J.GoaterR. J.AbudahabK.AanensenD. M.HarrisS. R. (2018). Phandango: an interactive viewer for bacterial population genomics.Bioinformatics34292–293. 10.1093/bioinformatics/btx61029028899HuX. D.LiuD.WangM. Y.YangL.WangM.ZhuQ. Y. (2011). Clade 2.3.2 avian influenza virus (H5N1), Qinghai Lake region, China, 2009-2010.Emerg. Infect. Dis.17560–562. 10.3201/eid1703.10094821392463KozlovA. M.DarribaD.FlouriT.MorelB.StamatakisA. (2019). RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference.Bioinformatics354453–4455. 10.1093/bioinformatics/btz30531070718LiH.DurbinR. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform.Bioinformatics251754–1760. 10.1093/bioinformatics/btp32419451168MontealegreM. C.RodriguezA. T.RoyS.HossainM. I.IslamM. A.LanzaV. F. (2020). High genomic diversity and heterogenous origins of pathogenic and antibiotic-resistant Escherichia coli in household settings represent a challenge to reducing transmission in low-income settings.mSphere5e704–e719. 10.1128/mSphere.00704-1931941809Nicolas-ChanoineM. H.RobertJ.ViganM.LaouenanC.BrisseS.MentreF. (2013). Different factors associated with CTX-M-producing ST131 and non-ST131 Escherichia coli clinical isolates.PLoS One8:e72191. 10.1371/journal.pone.007219124023730PageA. J.CumminsC. A.HuntM.WongV. K.ReuterS.HoldenM. T. (2015). Roary: rapid large-scale prokaryote pan genome analysis.Bioinformatics313691–3693. 10.1093/bioinformatics/btv42126198102PageA. J.TaylorB.DelaneyA. J.SoaresJ.SeemannT.KeaneJ. A. (2016). SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments.Microb. Genom.2:e000056.ProsserD. J.CuiP.TakekawaJ. Y.TangM. J.HouY. S.CollinsB. M. (2011). Wild bird migration across the qinghai-tibetan plateau: a transmission route for highly pathogenic H5N1.PLoS One6:e17622. 10.1371/journal.pone.001762221408010RodigerS.SchierackP.BohmA.NitschkeJ.BergerI.FrommelU. (2013). A highly versatile microscope imaging technology platform for the multiplex real-time detection of biomolecules and autoimmune antibodies.Adv. Biochem. Eng. Biotech.13335–74. 10.1007/10_2011_132SalipanteS. J.RoachD. J.KitzmanJ. O.SnyderM. W.StackhouseB.Butler-WuS. M. (2015). Large-scale genomic sequencing of extraintestinal pathogenic Escherichia coli strains.Genome Res.25119–128.SchauflerK.SemmlerT.PickardD. J.de ToroM.de la CruzF.WielerL. H. (2016). Carriage of extended-spectrum beta-lactamase-plasmids does not reduce fitness but enhances virulence in some strains of pandemic E. coli Lineages.Front. Microbiol.7:336. 10.3389/fmicb.2016.0033627014251SchauflerK.SemmlerT.WielerL. H.TrottD. J.PitoutJ.PeiranoG. (2019). Genomic and functional analysis of emerging virulent and multi-drug resistant E. coli lineage ST648.Antimicrob. Agents Chemother.63:e00243-19.SchauflerK.SemmlerT.WielerL. H.WohrmannM.BaddamR.AhmedN. (2015). Clonal spread and interspecies transmission of clinically relevant ESBL-producing Escherichia coli of ST410 – another successful pandemic clone?FEMS Microbiol. Ecol.92:fiv155. 10.1093/femsec/fiv15526656065SchauflerK.WielerL. H.SemmlerT.EwersC.GuentherS. (2013). ESBL-plasmids carrying toxin-antitoxin systems can be “cured” of wild-type Escherichia coli using a heat technique.Gut Pathog.5:34. 10.1186/1757-4749-5-3424245987ShepardS. M.DanzeisenJ. L.IsaacsonR. E.SeemannT.AchtmanM.JohnsonT. J. (2012). Genome sequences and phylogenetic analysis of K88-and F18-positive porcine enterotoxigenic Escherichia coli.J. Bacteriol.194395–405. 10.1128/jb.06225-1122081385TenaillonO.SkurnikD.PicardB.DenamurE. (2010). The population genetics of commensal Escherichia coli.Nat. Rev. Microbiol.8207–217. 10.1038/nrmicro229820157339WalkerB. J.AbeelT.SheaT.PriestM.AbouellielA.SakthikumarS. (2014). Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement.PLoS One9:e112963. 10.1371/journal.pone.011296325409509XiaP. P.SongY. J.ZouY. J.YangY.ZhuG. Q. (2015). F4(+) enterotoxigenic Escherichia coli (ETEC) adhesion mediated by the major fimbrial subunit FaeG.J. Basic Microb.551118–1124. 10.1002/jobm.20140090125847483