ABSTRACT
Variation in disease severity among Escherichia coli O157:H7 infections may result from differential expression of Shiga toxin 2 (Stx2). Eleven strains belonging to four prominent phylogenetic clades, including clade 8 strains representative of the 2006 U.S. spinach outbreak, were examined for stx2 expression by real-time PCR and western blot analysis. Clade 8 strains were shown to overexpress stx2 basally, and following induction with ciprofloxacin when compared to strains from clades 1-3. Differences in stx2 expression generally correlated with Stx2 protein levels. Single-nucleotide polymorphisms identified in regions upstream of stx2AB in clade 8 strains were largely absent in non-clade 8 strains. This study concludes that stx2 overexpression is common to strains from clade 8 associated with hemolytic uremic syndrome, and describes SNPs which may affect stx2 expression and which could be useful in the genetic differentiation of highly-virulent strains.
Subject(s)
Disease Outbreaks , Escherichia coli Infections/microbiology , Escherichia coli O157/pathogenicity , Gene Expression , Hemolytic-Uremic Syndrome/microbiology , Shiga Toxin 2/biosynthesis , Cluster Analysis , Escherichia coli Infections/complications , Escherichia coli Infections/epidemiology , Escherichia coli O157/classification , Escherichia coli O157/genetics , Escherichia coli O157/isolation & purification , Hemolytic-Uremic Syndrome/epidemiology , Humans , Molecular Typing , Polymorphism, Single Nucleotide , Spinacia oleracea/microbiology , United States/epidemiology , VirulenceABSTRACT
Escherichia coli O157:H7 strains often produce Shiga toxins encoded by genes on lambdoid bacteriophages that insert into multiple loci as prophages. O157 strains were classified into distinct clades that vary in virulence. Herein, we used PCR assays to examine Shiga toxin (Stx) prophage occupancy in yehV, argW, wrbA, and sbcB among 346 O157 strains representing nine clades. Overall, yehV was occupied in most strains (n = 334, 96.5%), followed by wrbA (n = 213, 61.6%), argW (n = 103, 29.8%), and sbcB (n = 93, 26.9%). Twelve occupancy profiles were identified that varied in frequency and differed across clades. Strains belonging to clade 8 were more likely to have occupied sbcB and argW sites compared to other clades (p < 0.0001), while clade 2 strains were more likely to have occupied wrbA sites (p < 0.0001). Clade 8 strains also had more than the expected number of occupied sites based on the presence of stx variants (p < 0.0001). Deletion of a 20 kb non-Stx prophage occupying yehV in a clade 8 strain resulted in an ~18-fold decrease in stx2 expression. These data highlight the complexity of Stx prophage integration and demonstrate that clade 8 strains, which were previously linked to hemolytic uremic syndrome, have unique Stx prophage occupancy profiles that can impact stx2 expression.
Subject(s)
Escherichia coli O157/virology , Prophages/physiology , Escherichia coli O157/genetics , Lysogeny , Shiga ToxinABSTRACT
The Escherichia coli O157 : H7 TW14359 strain was implicated in a multi-state outbreak in North America in 2006, which resulted in high rates of severe disease. Similarly, the O157 : H7 RIMD0509952 (Sakai) strain caused the largest O157 : H7 outbreak to date. Both strains were shown to represent divergent phylogenetic lineages. Here we compared global gene expression patterns before and after epithelial cell exposure, as well as the ability to adhere to and invade epithelial cells, between the two outbreak strains. Epithelial cell assays demonstrated a 2.5-fold greater adherence of the TW14359 strain relative to Sakai, while whole-genome microarrays detected significant differential expression of 914 genes, 206 of which had a fold change >/=1.5. Interestingly, most locus of enterocyte effacement (LEE) genes were upregulated in TW14359, whereas flagellar and chemotaxis genes were primarily upregulated in Sakai, suggesting discordant expression of these genes between the two strains. The Shiga toxin 2 genes were also upregulated in the TW14359 strain, as were several pO157-encoded genes that promote adherence, including type II secretion genes and their effectors stcE and adfO. Quantitative RT-PCR confirmed the expression differences detected in the microarray analysis, and expression levels were lower for a subset of LEE genes before versus after exposure to epithelial cells. In all, this study demonstrated the upregulation of major and ancillary virulence genes in TW14359 and of flagellar and chemotaxis genes in Sakai, under conditions that precede intimate bacterial attachment to epithelial cells. Differences in the level of adherence to epithelial cells were also observed, implying that these two phylogenetically divergent O157 : H7 outbreak strains vary in their ability to colonize, or initiate the disease process.
Subject(s)
Disease Outbreaks , Escherichia coli Infections/microbiology , Escherichia coli O157/pathogenicity , Gene Expression Regulation, Bacterial , Animals , Bacterial Adhesion/genetics , Cattle , Cell Line , Epithelial Cells/microbiology , Escherichia coli Infections/epidemiology , Escherichia coli Infections/genetics , Escherichia coli O157/genetics , Escherichia coli O157/metabolism , Escherichia coli Proteins/genetics , Escherichia coli Proteins/metabolism , Gene Expression Profiling , Humans , Molecular Sequence Data , North America/epidemiology , Oligonucleotide Array Sequence Analysis , Reverse Transcriptase Polymerase Chain Reaction , Shiga Toxin 2/genetics , Shiga Toxin 2/metabolism , Species Specificity , Virulence/geneticsABSTRACT
BACKGROUND: Evolutionary analyses of enterohemorrhagic Escherichia coli (EHEC) have identified two distantly related clonal groups: EHEC 1, including serotype O157:H7 and its inferred ancestor O55:H7; and EHEC 2, comprised of several serogroups (O26, O111, O118, etc.). These two clonal groups differ in their virulence and global distribution. Although several fully annotated genomic sequences exist for strains of serotype O157:H7, much less is known about the genomic composition of EHEC 2. In this study, we analyzed a set of 24 clinical EHEC 2 strains representing serotypes O26:H11, O111:H8/H11, O118:H16, O153:H11 and O15:H11 from humans and animals by comparative genomic hybridization (CGH) on an oligoarray based on the O157:H7 Sakai genome. RESULTS: Backbone genes, defined as genes shared by Sakai and K-12, were highly conserved in EHEC 2. The proportion of Sakai phage genes in EHEC 2 was substantially greater than that of Sakai-specific bacterial (non-phage) genes. This proportion was inverted in O55:H7, reiterating that a subset of Sakai bacterial genes is specific to EHEC 1. Split decomposition analysis of gene content revealed that O111:H8 was more genetically uniform and distinct from other EHEC 2 strains, with respect to the Sakai O157:H7 gene distribution. Serotype O26:H11 was the most heterogeneous EHEC 2 subpopulation, comprised of strains with the highest as well as the lowest levels of Sakai gene content conservation. Of the 979 parsimoniously informative genes, 15% were found to be compatible and their distribution in EHEC 2 clustered O111:H8 and O118:H16 strains by serotype. CGH data suggested divergence of the LEE island from the LEE1 to the LEE4 operon, and also between animal and human isolates irrespective of serotype. No correlation was found between gene contents and geographic locations of EHEC 2 strains. CONCLUSION: The gene content variation of phage-related genes in EHEC 2 strains supports the hypothesis that extensive modular shuffling of mobile DNA elements has occurred among EHEC strains. These results suggest that EHEC 2 is a multiform pathogenic clonal complex, characterized by substantial intra-serotype genetic variation. The heterogeneous distribution of mobile elements has impacted the diversification of O26:H11 more than other EHEC 2 serotypes.
Subject(s)
Comparative Genomic Hybridization , Enterohemorrhagic Escherichia coli/genetics , Genome, Bacterial , Animals , Bacterial Typing Techniques , DNA, Bacterial/genetics , Enterohemorrhagic Escherichia coli/classification , Genes, Bacterial , Genetic Variation , Humans , Oligonucleotide Array Sequence Analysis , Phylogeny , Prophages/geneticsABSTRACT
The gut microbiome is intimately related to human health, but it is not yet known which functional activities are driven by specific microorganisms' ecological configurations or transcription. We report a large-scale investigation of 372 human faecal metatranscriptomes and 929 metagenomes from a subset of 308 men in the Health Professionals Follow-Up Study. We identified a metatranscriptomic 'core' universally transcribed over time and across participants, often by different microorganisms. In contrast to the housekeeping functions enriched in this core, a 'variable' metatranscriptome included specialized pathways that were differentially expressed both across participants and among microorganisms. Finally, longitudinal metagenomic profiles allowed ecological interaction network reconstruction, which remained stable over the six-month timespan, as did strain tracking within and between participants. These results provide an initial characterization of human faecal microbial ecology into core, subject-specific, microorganism-specific and temporally variable transcription, and they differentiate metagenomically versus metatranscriptomically informative aspects of the human faecal microbiome.
Subject(s)
Feces/microbiology , Gene Expression Profiling , Metagenome , Microbiota , Aged , Aged, 80 and over , Follow-Up Studies , Gastrointestinal Microbiome , Gene Regulatory Networks , Humans , Longitudinal Studies , Male , Metagenomics , Phylogeny , Prospective StudiesABSTRACT
Characterizing the stability of the gut microbiome is important to exploit it as a therapeutic target and diagnostic biomarker. We metagenomically and metatranscriptomically sequenced the faecal microbiomes of 308 participants in the Health Professionals Follow-Up Study. Participants provided four stool samples-one pair collected 24-72 h apart and a second pair ~6 months later. Within-person taxonomic and functional variation was consistently lower than between-person variation over time. In contrast, metatranscriptomic profiles were comparably variable within and between subjects due to higher within-subject longitudinal variation. Metagenomic instability accounted for ~74% of corresponding metatranscriptomic instability. The rest was probably attributable to sources such as regulation. Among the pathways that were differentially regulated, most were consistently over- or under-transcribed at each time point. Together, these results suggest that a single measurement of the faecal microbiome can provide long-term information regarding organismal composition and functional potential, but repeated or short-term measures may be necessary for dynamic features identified by metatranscriptomics.
Subject(s)
Feces/microbiology , Gastrointestinal Microbiome , Gene Expression , Microbiota , Adult , Aged , Bacteria/classification , Cohort Studies , Follow-Up Studies , Gene Expression Profiling , Health Personnel/statistics & numerical data , High-Throughput Nucleotide Sequencing , Humans , Male , Metagenomics , Middle Aged , Prospective StudiesABSTRACT
Fluoridation of drinking water and dental products prevents dental caries primarily by inhibiting energy harvest in oral cariogenic bacteria (such as Streptococcus mutans and Streptococcus sanguinis), thus leading to their depletion. However, the extent to which oral and gut microbial communities are affected by host fluoride exposure has been underexplored. In this study, we modeled human fluoride exposures to municipal water and dental products by treating mice with low or high levels of fluoride over a 12-week period. We then used 16S rRNA gene amplicon and shotgun metagenomic sequencing to assess fluoride's effects on oral and gut microbiome composition and function. In both the low- and high-fluoride groups, several operational taxonomic units (OTUs) belonging to acidogenic bacterial genera (such as Parabacteroides, Bacteroides, and Bilophila) were depleted in the oral community. In addition, fluoride-associated changes in oral community composition resulted in depletion of gene families involved in central carbon metabolism and energy harvest (2-oxoglutarate ferredoxin oxidoreductase, succinate dehydrogenase, and the glyoxylate cycle). In contrast, fluoride treatment did not induce a significant shift in gut microbial community composition or function in our mouse model, possibly due to absorption in the upper gastrointestinal tract. Fluoride-associated perturbations thus appeared to have a selective effect on the composition of the oral but not gut microbial community in mice. Future studies will be necessary to understand possible implications of fluoride exposure for the human microbiome. IMPORTANCE Fluoride has been added to drinking water and dental products since the 1950s. The beneficial effects of fluoride on oral health are due to its ability to inhibit the growth of bacteria that cause dental caries. Despite widespread human consumption of fluoride, there have been only two studies of humans that considered the effect of fluoride on human-associated microbial communities, which are increasingly understood to play important roles in health and disease. Notably, neither of these studies included a true cross-sectional control lacking fluoride exposure, as study subjects continued baseline fluoride treatment in their daily dental hygiene routines. To our knowledge, this work (in mice) is the first controlled study to assess the independent effects of fluoride exposure on the oral and gut microbial communities. Investigating how fluoride interacts with host-associated microbial communities in this controlled setting represents an effort toward understanding how common environmental exposures may potentially influence health.
ABSTRACT
The emergence of novel pathogens poses a major public health threat causing widespread epidemics in susceptible populations. The Escherichia coli O104:H4 strain implicated in a 2011 outbreak in northern Germany caused the highest frequency of hemolytic uremic syndrome (HUS) and death ever recorded in a single E. coli outbreak. Therefore, it has been suggested that this strain is more virulent than other pathogenic E. coli (e.g., E. coli O157:H7). The E. coli O104:H4 outbreak strain possesses multiple virulence factors from both Shiga toxin (Stx)-producing E. coli (STEC) and enteroaggregative E. coli (EAEC), though the mechanism of pathogenesis is not known. Here, we demonstrate that E. coli O104:H4 produces a stable biofilm in vitro and that in vivo virulence gene expression is highest when E. coli O104:H4 overexpresses genes required for aggregation and exopolysaccharide production, a characteristic of bacterial cells residing within an established biofilm. Interrupting exopolysaccharide production and biofilm formation may therefore represent effective strategies for combating future E. coli O104:H4 infections.
Subject(s)
Biofilms/growth & development , Escherichia coli/genetics , Escherichia coli/physiology , Gene Expression Regulation, Bacterial , Animals , Disease Outbreaks , Escherichia coli/metabolism , Escherichia coli/pathogenicity , Escherichia coli O157/genetics , Escherichia coli O157/metabolism , Escherichia coli O157/pathogenicity , Escherichia coli O157/physiology , Kidney/microbiology , Mice , Polysaccharides, Bacterial/biosynthesis , Species SpecificityABSTRACT
BACKGROUND: Enterohemorrhagic Escherichia coli (EHEC) O157:H7, a food and waterborne pathogen, can be classified into nine phylogenetically distinct lineages, as determined by single nucleotide polymorphism genotyping. One lineage (clade 8) was found to be associated with hemolytic uremic syndrome (HUS), which can lead to kidney failure and death in some cases, particularly young children. Another lineage (clade 2) differs considerably in gene content and is phylogenetically distinct from clade 8, but caused significantly fewer cases of HUS in a prior study. Little is known, however, about how these two lineages vary with regard to phenotypic traits important for disease pathogenesis and in the expression of shared virulence genes. METHODOLOGY/PRINCIPAL FINDINGS: Here, we quantified the level of adherence to and invasion of MAC-T bovine epithelial cells, and examined the transcriptomes of 24 EHEC O157:H7 strains with varying Shiga toxin profiles from two common lineages. Adherence to epithelial cells was >2-fold higher for EHEC O157:H7 strains belonging to clade 8 versus clade 2, while no difference in invasiveness was observed between the two lineages. Whole-genome 70-mer oligo microarrays, which probe for 6088 genes from O157:H7 Sakai, O157:H7 EDL 933, pO157, and K12 MG1655, detected significant differential expression between clades in 604 genes following co-incubation with epithelial cells for 30 min; 186 of the 604 genes had a >1.5 fold change difference. Relative to clade 2, clade 8 strains showed upregulation of major virulence genes, including 29 of the 41 locus of enterocyte effacement (LEE) pathogenicity island genes, which are critical for adherence, as well as Shiga toxin genes and pO157 plasmid-encoded virulence genes. Differences in expression of 16 genes that encode colonization factors, toxins, and regulators were confirmed by qRT-PCR, which revealed a greater magnitude of change than microarrays. CONCLUSIONS/SIGNIFICANCE: These findings demonstrate that the EHEC O157:H7 lineage associated with HUS expresses higher levels of virulence genes and has an enhanced ability to attach to epithelial cells relative to another common lineage.