Phylogenetic incongruence in E. coli O104: understanding the evolutionary relationships of emerging pathogens in the face of homologous recombination

Genomic Profiling of Multidrug-Resistant Swine Escherichia coli and Clonal Relationship to Human Isolates in Peru

The misuse of antibiotics is accelerating antimicrobial resistance (AMR) in Escherichia coli isolated from farm animals. The genomes of ten multidrug-resistant (MDR) E. coli isolates from pigs were analyzed to determine their sequence types, serotypes, virulence, and AMR genes (ARGs). Additionally, the relationship was evaluated adding all the available genomes of Peruvian E. coli from humans using the cgMLST + HierCC scheme. Two aEPEC O186:H11-ST29 were identified, of which H11 and ST29 are reported in aEPEC isolates from different sources. An isolate ETEC-O149:H10-ST100 was identified, considered a high-risk clone that is frequently reported in different countries as a cause of diarrhea in piglets. One ExPEC O101:H11-ST167 was identified, for which ST167 is an international high-risk clone related to urinary infections in humans. We identified many ARGs, including extended-spectrum β-lactamase genes, and one ETEC harboring the mcr-1 gene. CgMLST + HierCC analysis differentiated three clusters, and in two, the human isolates were grouped with those of swine in the same cluster. We observed that Peruvian swine MDR E. coli cluster with Peruvian E. coli isolates from healthy humans and from clinical cases, which is of great public health concern and evidence that AMR surveillance should be strengthened based on the One Health approach.

Recombination in Bacterial Genomes: Evolutionary Trends

Bacterial organisms have undergone homologous recombination (HR) and horizontal gene transfer (HGT) multiple times during their history. These processes could increase fitness to new environments, cause specialization, the emergence of new species, and changes in virulence. Therefore, comprehensive knowledge of the impact and intensity of genetic exchanges and the location of recombination hotspots on the genome is necessary for understanding the dynamics of adaptation to various conditions. To this end, we aimed to characterize the functional impact and genomic context of computationally detected recombination events by analyzing genomic studies of any bacterial species, for which events have been detected in the last 30 years. Genomic loci where the transfer of DNA was detected pertained to mobile genetic elements (MGEs) housing genes that code for proteins engaged in distinct cellular processes, such as secretion systems, toxins, infection effectors, biosynthesis enzymes, etc. We found that all inferences fall into three main lifestyle categories, namely, ecological diversification, pathogenesis, and symbiosis. The latter primarily exhibits ancestral events, thus, possibly indicating that adaptation appears to be governed by similar recombination-dependent mechanisms.

Circulation of Shiga Toxin-Producing Escherichia coli Phylogenetic Group B1 Strains Between Calve Stable Manure and Pasture Land With Grazing Heifers

Escherichia coli strains carrying Shiga toxins 1 and 2 (stx₁ and stx₂), intimin (eae), and hemolysin (ehxA) production genes were found in grass shoot, rhizosphere soil, and stable manure samples from a small-scale cattle farm located at the center of Netherlands, using cultivation-dependent and -independent microbiological detection techniques. Pasture land with grazing heifers in the first year of sampling in 2014 and without grazing cattle in 2015 was physically separated from the stable that housed rose calves during both years. Manure from the stable was applied to pasture via injection into soil once per year in early spring. Among a variety of 35 phylogenetic distinctly related E. coli strains, one large group consisting of 21 closely resembling E. coli O150:H2 (18), O98:H21 (2), and O84:H2 (1) strains, all belonging to phylogenetic group B1 and carrying all screened virulence traits, was found present on grass shoots (10), rhizosphere soil (3), and stable manure (8) in 2014, but not anymore in 2015 when grazing heifers were absent. Presence and absence of these strains, obtained via enrichments, were confirmed via molecular detection using PCR-NALFIA in all ecosystems in both years. We propose that this group of Shiga toxin-producing E. coli phylogenetic group B1 strains was originally introduced via stable manure injection into the pasture. Upon grazing, these potential pathogens proliferated in the intestinal track systems of the heifers resulting in defecation with higher loads of the STEC strain onto the grass cover. The STEC strain was further smeared over the field via the hooves of the heifers resulting in augmentation of the potential pathogen in the pasture in 2014, whereas in 2015, in the absence of heifers, no augmentation occurred and only a more diverse group of potentially mild virulent E. coli phylogenetic group A and B1 strains, indigenous to pasture plants, remained present. Via this model, it was postulated that human pathogens can circulate between plants and farm animals, using the plant as an alternative ecosystem. These data indicate that grazed pasture must be considered as a potential carrier of human pathogenic E. coli strains and possibly also of other pathogens.

Differential transcriptome analysis of enterohemorrhagic Escherichia coli strains reveals differences in response to plant-derived compounds

BACKGROUND

Several serious vegetable-associated outbreaks of enterohemorrhagic Escherichia coli (EHEC) infections have occurred during the last decades. In this context, vegetables have been suggested to function as secondary reservoirs for EHEC strains. Increased knowledge about the interaction of EHEC with plants including gene expression patterns in response to plant-derived compounds is required. In the current study, EHEC O157:H7 strain Sakai, EHEC O157:H- strain 3072/96, and the EHEC/enteroaggregative E. coli (EAEC) hybrid O104:H4 strain C227-11φcu were grown in lamb's lettuce medium and in M9 minimal medium to study the differential transcriptional response of these strains to plant-derived compounds with RNA-Seq technology.

RESULTS

Many genes involved in carbohydrate degradation and peptide utilization were similarly upregulated in all three strains, suggesting that the lamb's lettuce medium provides sufficient nutrients for proliferation. In particular, the genes galET and rbsAC involved in galactose metabolism and D-ribose catabolism, respectively, were uniformly upregulated in the investigated strains. The most prominent differences in shared genome transcript levels were observed for genes involved in the expression of flagella. Transcripts of all three classes of the flagellar hierarchy were highly abundant in strain C227-11φcu. Strain Sakai expressed only genes encoding the basal flagellar structure. In addition, both strains showed increased motility in presence of lamb's lettuce extract. Moreover, strain 3072/96 showed increased transcription activity for genes encoding the type III secretion system (T3SS) including effectors, and was identified as a powerful biofilm-producer in M9 minimal medium.

CONCLUSION

The current study provides clear evidence that EHEC and EHEC/EAEC strains are able to adjust their gene expression patterns towards metabolization of plant-derived compounds, demonstrating that they may proliferate well in a plant-associated environment. Moreover, we propose that flagella and other surface structures play a fundamental role in the interaction of EHEC and EHEC/EAEC with plants.

Pandemic fluoroquinolone resistant Escherichia coli clone ST1193 emerged via simultaneous homologous recombinations in 11 gene loci

Global growth in antibiotic resistance is a major social problem. A high level of resistance to fluoroquinolones requires the concurrent presence of at least 3 mutations in the target proteins-2 in DNA gyrase (GyrA) and 1 in topoisomerase IV (ParC), which occur in a stepwise manner. In the Escherichia coli chromosome, the gyrA and parC loci are positioned about 1 Mb away from each other. Here we show that the 3 fluoroquinolone resistance mutations are tightly associated genetically in naturally occurring strains. In the latest pandemic uropathogenic and multidrug-resistant E. coli clonal group ST1193, the mutant variants of gyrA and parC were acquired not by a typical gradual, stepwise evolution but all at once. This happened as part of 11 simultaneous homologous recombination events involving 2 phylogenetically distant strains of E. coli, from an uropathogenic clonal complex ST14 and fluoroquinolone-resistant ST10. The gene exchanges swapped regions between 0.5 and 139 Kb in length (183 Kb total) spread along 976 Kb of chromosomal DNA around and between gyrA and parC loci. As a result, all 3 fluoroquinolone resistance mutations in GyrA and ParC have simultaneously appeared in ST1193. Based on molecular clock estimates, this potentially happened as recently as <12 y ago. Thus, naturally occurring homologous recombination events between 2 strains can involve numerous chromosomal gene locations simultaneously, resulting in the transfer of distant but tightly associated genetic mutations and emergence of a both highly pathogenic and antibiotic-resistant strain with a rapid global spread capability.

Genetic Drift and Indel Mutation in the Evolution of Yeast Mitochondrial Genome Size

Mitochondrial genomes (mitogenomes) are remarkably diverse in genome size and organization, but the origins of dynamic mitogenome architectures are still poorly understood. For instance, the mutational burden hypothesis postulates that the drastic difference between large plant mitogenomes and streamlined animal mitogenomes can be driven by their different mutation rates. However, inconsistent trends between mitogenome sizes and mutation rates have been documented in several lineages. These conflicting results highlight the need of systematic and sophisticated investigations on the evolution and diversity of mitogenome architecture. This study took advantage of the strikingly variable mitogenome size among different yeast species and also among intraspecific strains, examined sequence dynamics of introns, GC-clusters, tandem repeats, mononucleotide repeats (homopolymers) and evaluated their contributions to genome size variation. The contributions of these sequence features to mitogenomic variation are dependent on the timescale, over which extant genomes evolved from their last common ancestor, perhaps due to a combination of different turnover rates of mobile sequences, variable insertion spaces, and functional constraints. We observed a positive correlation between mitogenome size and the level of genetic drift, suggesting that mitogenome expansion in yeast is likely driven by multiple types of sequence insertions in a primarily nonadaptive manner. Although these cannot be explained directly by the mutational burden hypothesis, our results support an important role of genetic drift in the evolution of yeast mitogenomes.

Comparative Genomics and Characterization of the Late Promoter pR' from Shiga Toxin Prophages in Escherichia coli

Shiga-toxin producing Escherichia coli (STEC) causes human illness ranging from mild diarrhea to death. The bacteriophage encoded stx genes are located in the late transcription region, downstream of the antiterminator Q. The transcription of the stx genes is directly under the control of the late promoter pR’, thus the sequence diversity of the region between Q and stx, here termed the pR’ region, may affect Stx toxin production. Here, we compared the gene structure of the pR’ region and the stx subtypes of nineteen STECs. The sequence alignment and phylogenetic analysis suggested that the pR’ region tends to be more heterogeneous than the promoter itself, even if the prophages harbor the same stx subtype. Furthermore, we established and validated transcriptional fusions of the pR’ region to the DsRed reporter gene using mitomycin C (MMC) induction. Finally, these constructs were transformed into native and non-native strains and examined with flow cytometry. The results showed that induction levels changed when pR’ regions were placed under different regulatory systems. Moreover, not every stx gene could be induced in its native host bacteria. In addition to the functional genes, the diversity of the pR’ region plays an important role in determining the level of toxin induction.

Australian porcine clonal complex 10 (CC10) Escherichia coli belong to multiple sublineages of a highly diverse global CC10 phylogeny

We recently identified clonal complex 10 (CC10) Escherichia coli as the predominant clonal group in two populations of healthy Australian food-production pigs. CC10 are highly successful, colonizing humans, food-production animals, fresh produce and environmental niches. Furthermore, E. coli within CC10 are frequently drug resistant and increasingly reported as human and animal extra-intestinal pathogens. In order to develop a high-resolution global phylogeny and determine the repertoire of antimicrobial-resistance genes, virulence-associated genes and plasmid types within this clonal group, we downloaded 228 publicly available CC10 short-read genome sequences for comparison with 20 porcine CC10 we have previously described. Core genome single nucleotide polymorphism phylogeny revealed a highly diverse global phylogeny consisting of multiple lineages that did not cluster by geography or source of the isolates. Australian porcine strains belonged to several of these divergent lineages, indicative that CC10 is present in these animals due to multiple colonization events. Differences in resistance gene and plasmid carriage between porcine strains and the global collection highlighted the role of lateral gene transfer in the evolution of CC10 strains. Virulence profiles typical of extra-intestinal pathogenic E. coli were present in both Australian porcine strains and the broader collection. As both the core phylogeny and accessory gene characteristics appeared unrelated to the geography or source of the isolates, it is likely that the global expansion of CC10 is not a recent event and may be associated with faecal carriage in humans.

Whole-Genome Multilocus Sequence Typing of Extended-Spectrum-Beta-Lactamase-Producing Enterobacteriaceae

Molecular typing has become indispensable in the detection of nosocomial transmission of bacterial pathogens and the identification of sources and routes of transmission in outbreak settings, but current methods are labor-intensive, are difficult to standardize, or have limited resolution. Whole-genome multilocus sequence typing (wgMLST) has emerged as a whole-genome sequencing (WGS)-based gene-by-gene typing method that may overcome these limitations and has been applied successfully for several species in outbreak settings. In this study, genus-, genetic-complex-, and species-specific wgMLST schemes were developed for Citrobacter spp., the Enterobacter cloacae complex, Escherichia coli, Klebsiella oxytoca, and Klebsiella pneumoniae and used to type a national collection of 1,798 extended-spectrum-beta-lactamase-producing Enterobacteriaceae (ESBL-E) isolates obtained from patients in Dutch hospitals. Genus-, genetic-complex-, and species-specific thresholds for genetic distance that accurately distinguish between epidemiologically related and unrelated isolates were defined for Citrobacter spp., the E. cloacae complex, E. coli, and K. pneumoniae wgMLST was shown to have higher discriminatory power and typeability than in silico MLST. In conclusion, the wgMLST schemes developed in this study facilitate high-resolution WGS-based typing of the most prevalent ESBL-producing species in clinical practice and may contribute to further elucidation of the complex epidemiology of antimicrobial-resistant Enterobacteriaceae wgMLST opens up possibilities for the creation of a Web-accessible database for the global surveillance of ESBL-producing bacterial clones.

Escherichia coli O104 in Feedlot Cattle Feces: Prevalence, Isolation and Characterization

Escherichia coli O104:H4, an hybrid pathotype of Shiga toxigenic and enteroaggregative E. coli, involved in a major foodborne outbreak in Germany in 2011, has not been detected in cattle feces. Serogroup O104 with H type other than H4 has been reported to cause human illnesses, but their prevalence and characteristics in cattle have not been reported. Our objectives were to determine the prevalence of E. coli O104 in feces of feedlot cattle, by culture and PCR detection methods, and characterize the isolated strains. Rectal fecal samples from a total of 757 cattle originating from 29 feedlots were collected at a Midwest commercial slaughter plant. Fecal samples, enriched in E. coli broth, were subjected to culture and PCR methods of detection. The culture method involved immunomagnetic separation with O104-specific beads and plating on a selective chromogenic medium, followed by serogroup confirmation of pooled colonies by PCR. If pooled colonies were positive for the wzxO104 gene, then colonies were tested individually to identify wzxO104-positive serogroup and associated genes of the hybrid strains. Extracted DNA from feces were also tested by a multiplex PCR to detect wzxO104-positive serogroup and associated major genes of the O104 hybrid pathotype. Because wzxO104 has been shown to be present in E. coli O8/O9/O9a, wzxO104-positive isolates and extracted DNA from fecal samples were also tested by a PCR targeting wbdDO8/O9/O9a, a gene specific for E. coli O8/O9/O9a serogroups. Model-adjusted prevalence estimates of E. coli O104 (positive for wzxO104 and negative for wbdDO8/O9/O9a) at the feedlot level were 5.7% and 21.2%, and at the sample level were 0.5% and 25.9% by culture and PCR, respectively. The McNemar's test indicated that there was a significant difference (P < 0.01) between the proportions of samples that tested positive for wzxO104 and samples that were positive for wzxO104, but negative for wbdDO8/O9/O9a by PCR and culture methods. A total of 143 isolates, positive for the wzxO104, were obtained in pure culture from 146 positive fecal samples. Ninety-two of the 143 isolates (64.3%) also tested positive for the wbdDO8/O9/O9a, indicating that only 51 (35.7%) isolates truly belonged to the O104 serogroup (positive for wzxO104 and negative for wbdDO8/O9/O9a). All 51 isolates tested negative for eae, and 16 tested positive for stx1 gene of the subtype 1c. Thirteen of the 16 stx1-positive O104 isolates were from one feedlot. The predominant serotype was O104:H7. Pulsed-field gel electrophoresis analysis indicated that stx1-positive O104:H7 isolates had 62.4% homology to the German outbreak strain and 67.9% to 77.5% homology to human diarrheagenic O104:H7 strains. The 13 isolates obtained from the same feedlot were of the same PFGE subtype with 100% Dice similarity. Although cattle do not harbor the O104:H4 pathotype, they do harbor and shed Shiga toxigenic O104 in the feces and the predominant serotype was O104:H7.

Extensive genomic variation within clonal bacterial groups resulted from homologous recombination

Due to divergence, genetic variation is generally believed to be high among distantly related strains, low among closely related ones and little or none within the same classified clonal groups. Several recent genome-wide studies, however, revealed that significant genetic variation resides in a considerable number of genes among strains with identical MLST (Multilocus sequence typing) types and much of the variation was introduced by homologous recombination. Recognizing and understanding genomic variation within clonal bacterial groups could shed new light on the evolutionary path of infectious agents and the emergence of particularly pathogenic or virulent variants. This commentary presents our recent contributions to this line of work.

The epidemic of extended-spectrum-β-lactamase-producing Escherichia coli ST131 is driven by a single highly pathogenic subclone, H30-Rx

The Escherichia coli sequence type 131 (ST131) clone is notorious for extraintestinal infections, fluoroquinolone resistance, and extended-spectrum beta-lactamase (ESBL) production, attributable to a CTX-M-15-encoding mobile element. Here, we applied pulsed-field gel electrophoresis (PFGE) and whole-genome sequencing to reconstruct the evolutionary history of the ST131 clone. PFGE-based cluster analyses suggested that both fluoroquinolone resistance and ESBL production had been acquired by multiple ST131 sublineages through independent genetic events. In contrast, the more robust whole-genome-sequence-based phylogenomic analysis revealed that fluoroquinolone resistance was confined almost entirely to a single, rapidly expanding ST131 subclone, designated H30-R. Strikingly, 91% of the CTX-M-15-producing isolates also belonged to a single, well-defined clade nested within H30-R, which was named H30-Rx due to its more extensive resistance. Despite its tight clonal relationship with H30Rx, the CTX-M-15 mobile element was inserted variably in plasmid and chromosomal locations within the H30-Rx genome. Screening of a large collection of recent clinical E. coli isolates both confirmed the global clonal expansion of H30-Rx and revealed its disproportionate association with sepsis (relative risk, 7.5; P < 0.001). Together, these results suggest that the high prevalence of CTX-M-15 production among ST131 isolates is due primarily to the expansion of a single, highly virulent subclone, H30-Rx.

We applied an advanced genomic approach to study the recent evolutionary history of one of the most important Escherichia coli strains in circulation today. This strain, called sequence type 131 (ST131), causes multidrug-resistant bladder, kidney, and bloodstream infections around the world. The rising prevalence of antibiotic resistance in E. coli is making these infections more difficult to treat and is leading to increased mortality. Past studies suggested that many different ST131 strains gained resistance to extended-spectrum cephalosporins independently. In contrast, our research indicates that most extended-spectrum-cephalosporin-resistant ST131 strains belong to a single highly pathogenic subclone, called H30-Rx. The clonal nature of H30-Rx may provide opportunities for vaccine or transmission prevention-based control strategies, which could gain importance as H30-Rx and other extraintestinal pathogenic E. coli subclones become resistant to our best antibiotics.

When whole-genome alignments just won't work: kSNP v2 software for alignment-free SNP discovery and phylogenetics of hundreds of microbial genomes

Effective use of rapid and inexpensive whole genome sequencing for microbes requires fast, memory efficient bioinformatics tools for sequence comparison. The kSNP v2 software finds single nucleotide polymorphisms (SNPs) in whole genome data. kSNP v2 has numerous improvements over kSNP v1 including SNP gene annotation; better scaling for draft genomes available as assembled contigs or raw, unassembled reads; a tool to identify the optimal value of k; distribution of packages of executables for Linux and Mac OS X for ease of installation and user-friendly use; and a detailed User Guide. SNP discovery is based on k-mer analysis, and requires no multiple sequence alignment or the selection of a single reference genome. Most target sets with hundreds of genomes complete in minutes to hours. SNP phylogenies are built by maximum likelihood, parsimony, and distance, based on all SNPs, only core SNPs, or SNPs present in some intermediate user-specified fraction of targets. The SNP-based trees that result are consistent with known taxonomy. kSNP v2 can handle many gigabases of sequence in a single run, and if one or more annotated genomes are included in the target set, SNPs are annotated with protein coding and other information (UTRs, etc.) from Genbank file(s). We demonstrate application of kSNP v2 on sets of viral and bacterial genomes, and discuss in detail analysis of a set of 68 finished E. coli and Shigella genomes and a set of the same genomes to which have been added 47 assemblies and four “raw read” genomes of H104:H4 strains from the recent European E. coli outbreak that resulted in both bloody diarrhea and hemolytic uremic syndrome (HUS), and caused at least 50 deaths.

Recent advances in understanding enteric pathogenic Escherichia coli

Although Escherichia coli can be an innocuous resident of the gastrointestinal tract, it also has the pathogenic capacity to cause significant diarrheal and extraintestinal diseases. Pathogenic variants of E. coli (pathovars or pathotypes) cause much morbidity and mortality worldwide. Consequently, pathogenic E. coli is widely studied in humans, animals, food, and the environment. While there are many common features that these pathotypes employ to colonize the intestinal mucosa and cause disease, the course, onset, and complications vary significantly. Outbreaks are common in developed and developing countries, and they sometimes have fatal consequences. Many of these pathotypes are a major public health concern as they have low infectious doses and are transmitted through ubiquitous mediums, including food and water. The seriousness of pathogenic E. coli is exemplified by dedicated national and international surveillance programs that monitor and track outbreaks; unfortunately, this surveillance is often lacking in developing countries. While not all pathotypes carry the same public health profile, they all carry an enormous potential to cause disease and continue to present challenges to human health. This comprehensive review highlights recent advances in our understanding of the intestinal pathotypes of E. coli.

Bloody coli: a gene cocktail in Escherichia coli O104:H4

A recent study published in mBio [Y. H. Grad et al., mBio 4(1):e00452-12, 2013] indicates that a rapid introgressive evolution has occurred in Escherichia coli O104:H4 by sequential acquisition of foreign genetic material involving pathogenicity traits. O104 genetic promiscuity cannot be readily explained by high population sizes. However, extensive interactions leading to cumulative assemblies of pathogenicity genes might be assured by small K-strategist populations exploiting particular intestinal niches. Next-generation sequencing technologies will be critical to detect particular “gene cocktails” as potentially pathogenic ensembles and to predict the risk of future outbreaks.

Homologous recombination drives both sequence diversity and gene content variation in Neisseria meningitidis

The study of genetic and phenotypic variation is fundamental for understanding the dynamics of bacterial genome evolution and untangling the evolution and epidemiology of bacterial pathogens. Neisseria meningitidis (Nm) is among the most intriguing bacterial pathogens in genomic studies due to its dynamic population structure and complex forms of pathogenicity. Extensive genomic variation within identical clonal complexes (CCs) in Nm has been recently reported and suggested to be the result of homologous recombination, but the extent to which recombination contributes to genomic variation within identical CCs has remained unclear. In this study, we sequenced two Nm strains of identical serogroup (C) and multi-locus sequence type (ST60), and conducted a systematic analysis with an additional 34 Nm genomes. Our results revealed that all gene content variation between the two ST60 genomes was introduced by homologous recombination at the conserved flanking genes, and 94.25% or more of sequence divergence was caused by homologous recombination. Recombination was found in genes associated with virulence factors, antigenic outer membrane proteins, and vaccine targets, suggesting an important role of homologous recombination in rapidly altering the pathogenicity and antigenicity of Nm. Recombination was also evident in genes of the restriction and modification systems, which may undermine barriers to DNA exchange. In conclusion, homologous recombination can drive both gene content variation and sequence divergence in Nm. These findings shed new light on the understanding of the rapid pathoadaptive evolution of Nm and other recombinogenic bacterial pathogens.