Molecular and Functional Profiling of the Polyamine Content in Enteroinvasive E. coli?: Looking into the Gap between Commensal E. coli and Harmful Shigella

Polyamines are small molecules associated with a wide variety of physiological functions. Bacterial pathogens have developed subtle strategies to exploit polyamines or manipulate polyamine-related processes to optimize fitness within the host. During the transition from its innocuous E. coli ancestor, Shigella, the aetiological agent of bacillary dysentery, has undergone drastic genomic rearrangements affecting the polyamine profile. A pathoadaptation process involving the speG gene and the cad operon has led to spermidine accumulation and loss of cadaverine. While a higher spermidine content promotes the survival of Shigella within infected macrophages, the lack of cadaverine boosts the pathogenic potential of the bacterium in host tissues. Enteroinvasive E. coli (EIEC) display the same pathogenicity process as Shigella, but have a higher infectious dose and a higher metabolic activity. Pathoadaption events affecting the cad locus have occurred also in EIEC, silencing cadaverine production. Since EIEC are commonly regarded as evolutionary intermediates between E. coli and Shigella, we investigated on their polyamine profile in order to better understand which changes have occurred along the path to pathogenicity. By functional and molecular analyses carried out in EIEC strains belonging to different serotypes, we show that speG has been silenced in one strain only, favouring resistance to oxidative stress conditions and survival within macrophages. At the same time, we observe that the content of spermidine and putrescine, a relevant intermediate in the synthesis of spermidine, is higher in all strains as compared to E. coli. This may represent an evolutionary response to the lack of cadaverine. Indeed, restoring cadaverine synthesis decreases the expression of the speC gene, whose product affects putrescine production. In the light of these results, we discuss the possible impact of pathoadaptation events on the evolutionary emergence of a polyamine profile favouring to the pathogenic lifestyle of Shigella and EIEC.     

Isolation of Shigella dysenteriae Type 1 and S. flexneri Strains from Surface Waters in Bangladesh: Comparative Molecular Analysis of Environmental Shigella Isolates versus Clinical Strains

Bacillary dysentery caused by Shigella species is a public health problem in developing countries including Bangladesh. Although, shigellae-contaminated food and drinks are often the source of the epidemic's spread, the possible presence of the pathogen and transmission of it through environmental waters have not been adequately examined. We analyzed surface waters collected in Dhaka, Bangladesh, for the presence of shigellae by a combination of PCR assays followed by concentration and culturing of PCR-positive samples. Analysis of 128 water samples by PCR assays for Shigella-specific virulence genes including ipaBCD, ipaH, and stx1 identified 14 (10.9%) samples which were positive for one or more of these virulence genes. Concentration of the PCR-positive samples by filtration followed by culturing identified live Shigella species in 11 of the 14 PCR-positive samples. Analysis of rRNA gene restriction patterns (ribotype) showed that the environmental isolates shared ribotypes with a collection of clinical isolates, but in contrast to the clinical isolates, 10 of the 11 environmental isolates were either negative or carried deletions in the plasmid-encoded invasion-associated genes ipaB, ipaC, and ipaD. However, all environmental Shigella isolates were positive for the chromosomal multicopy invasion-associated gene ipaH and all Shigella dysenteriae type 1 isolates were positive for the stx1 gene in addition to ipaH. This study demonstrated the presence of Shigella in the aquatic environment and dispersion of different virulence genes among these isolates which appear to constitute an environmental reservoir of Shigella-specific virulence genes. Since critical virulence genes in Shigella are carried by plasmids or mobile genetic elements, the environmental gene pool may contribute to an optimum combination of genes, causing the emergence of virulent Shigella strains which is facilitated in particular by close contact of the population with surface waters in Bangladesh.     

Gene decay in Shigella as an incipient stage of host-adaptation

Background

Many facultative bacterial pathogens have undergone extensive gene decay processes, possibly due to lack of selection pressure during evolutionary conversion from free-living to intracellular lifestyle. Shigella, the causative agents of human shigellosis, have arisen from different E. coli-like ancestors independently by convergent paths. As these bacteria all have lost large numbers of genes by mutation or deletion, they can be used as ideal models for systematically studying the process of gene function loss in different bacteria living under similar selection pressures.

Methodologies/Principal Findings

We compared the sequenced Shigella genomes and re-defined decayed genes (pseudogenes plus deleted genes) in these bacteria. Altogether, 85 genes are commonly decayed in the five analyzed Shigella strains and 1456 genes are decayed in at least one Shigella strain. Genes coding for carbon utilization, cell motility, transporter or membrane proteins are prone to be inactivated. Decayed genes tend to concentrate in certain operons rather than distribute averagely across the whole genome. Genes in the decayed operon accumulated more non-synonymous mutations than the rest genes and meanwhile have lower expression levels.

Conclusions

Different Shigella lineages underwent convergent gene decay processes, and inactivation of one gene would lead to a lesser selection pressure for the other genes in the same operon. The pool of superfluous genes for Shigella may contain at least two thousand genes and the gene decay processes may still continue in Shigella until a minimum genome harboring only essential genes is reached.     

Evolution of Pan-Genomes of Escherichia coli,Shigella spp., and Salmonella enterica

Multiple sequencing of genomes belonging to a bacterial species allows one to analyze and compare statistics and dynamics of the gene complements of species, their pan-genomes. Here, we analyzed multiple genomes of Escherichia coli, Shigella spp., and Salmonella enterica. We demonstrate that the distribution of the number of genomes harboring a gene is well approximated by a sum of two power functions, describing frequent genes (present in many strains) and rare genes (present in few strains). The virtual absence of Shigella-specific genes not present in E. coli genomes confirms previous observations that Shigella is not an independent genus. While the pan-genome size is increasing with each new strain, the number of genes present in a fixed fraction of strains stabilizes quickly. For instance, slightly fewer than 4,000 genes are present in at least half of any group of E. coli genomes. Comparison of S. enterica and E. coli pan-genomes revealed the existence of a common periphery, that is, genes present in some but not all strains of both species. Analysis of phylogenetic trees demonstrates that rare genes from the periphery likely evolve under horizontal transfer, whereas frequent periphery genes may have been inherited from the periphery genome of the common ancestor.     

Revisiting the Molecular Evolutionary History of <Emphasis Type="Italic">Shigella</Emphasis> spp.

The theory that Shigella is derived from multiple independent origins of Escherichia coli (Pupo et al. 2000) has been challenged by recent findings that the virulence plasmids (VPs) and the chromosomes share a similar evolutionary history (Escobar-Paramo et al. 2003), which suggests that an ancestral VP entered an E. coli strain only once, which gave rise to Shigella spp. In an attempt to resolve these conflicting theories, we constructed three phylogenetic trees in this study: a robust chromosomal tree using 23 housekeeping genes from 46 strains of Shigella and enteroinvasive E. coli (EIEC), a chromosomal tree using 4 housekeeping genes from 19 EcoR strains and 46 Shigella/EIEC strains, and a VP tree using 5 genes outside of the VP cell-entry region from 38 Shigella/EIEC strains. Both chromosomal trees group Shigella into three main clusters and five outliers, and strongly suggest that Shigella has multiple origins within E. coli. Most strikingly, the VP tree shows that the VPs from two main Shigella clusters, C1 and C2, are more closely related, which contradicts the chromosomal trees that place C2 and C3 next to each other but C1 at a distance. Additionally, we have identified a complete tra operon of the F-plasmid in the genome sequence of an EIEC strain and found that two other EIEC strains are also likely to possess a complete tra operon. All lines of evidence support an alternative multiorigin theory that transferable diverse ancestral VPs entered diverse origins of E. coli multiple times during a prolonged period of time, resulting in Shigella species with diverse genomes but similar pathogenic properties. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Martin Kreitman] Jian Yang and Huan Nie contributed equally to this work.     

Comparison of 61 Sequenced Escherichia coli Genomes

Escherichia coli is an important component of the biosphere and is an ideal model for studies of processes involved in bacterial genome evolution. Sixty-one publically available E. coli and Shigella spp. sequenced genomes are compared, using basic methods to produce phylogenetic and proteomics trees, and to identify the pan- and core genomes of this set of sequenced strains. A hierarchical clustering of variable genes allowed clear separation of the strains into clusters, including known pathotypes; clinically relevant serotypes can also be resolved in this way. In contrast, when in silico MLST was performed, many of the various strains appear jumbled and less well resolved. The predicted pan-genome comprises 15,741 gene families, and only 993 (6%) of the families are represented in every genome, comprising the core genome. The variable or ‘accessory’ genes thus make up more than 90% of the pan-genome and about 80% of a typical genome; some of these variable genes tend to be co-localized on genomic islands. The diversity within the species E. coli, and the overlap in gene content between this and related species, suggests a continuum rather than sharp species borders in this group of Enterobacteriaceae.     

Surface display of the HPV L1 capsid protein by the autotransporter Shigella IcsA

Autotransporters have become attractive tools for surface expression of foreign proteins in Gram-negative bacteria. In this study, the Shigella autotransporter IcsA, has been exploited to express the human papillomavirus (HPV) type 16 L1 capsid protein in Shigella sonnei and Escherichia coli. The L1 gene was fused in-frame to replace the coding sequence of the IcsA passenger domain that is responsible for actin-based motility. The resultant hybrid protein could be detected by an anti-L1 antibody on the surface of S. sonnei and E. coli. In E. coli, the protein was expressed on the entire surface of the bacterium. In contrast, the protein was detected mainly at one pole of the Shigella bacterium. However, the protein became evenly distributed on the surface of the Shigella bacterium when the icsP gene was removed. Our study demonstrated the possibility of exploiting autotransporters for surface expression of large, heterologous viral proteins, which may be a useful strategy for vaccine development.     

Comparative Genomics and Characterization of Hybrid Shigatoxigenic and Enterotoxigenic Escherichia coli (STEC/ETEC) Strains

Background

Shigatoxigenic Escherichia coli (STEC) and enterotoxigenic E. coli (ETEC) cause serious foodborne infections in humans. These two pathogroups are defined based on the pathogroup-associated virulence genes: stx encoding Shiga toxin (Stx) for STEC and elt encoding heat-labile and/or est encoding heat-stable enterotoxin (ST) for ETEC. The study investigated the genomics of STEC/ETEC hybrid strains to determine their phylogenetic position among E. coli and to define the virulence genes they harbor.

Methods

The whole genomes of three STEC/ETEC strains possessing both stx and est genes were sequenced using PacBio RS sequencer. Two of the strains were isolated from the patients, one with hemolytic uremic syndrome, and one with diarrhea. The third strain was of bovine origin. Core genome analysis of the shared chromosomal genes and comparison with E. coli and Shigella spp. reference genomes was performed to determine the phylogenetic position of the STEC/ETEC strains. In addition, a set of virulence genes and ETEC colonization factors were extracted from the genomes. The production of Stx and ST were studied.

Results

The human STEC/ETEC strains clustered with strains representing ETEC, STEC, enteroaggregative E. coli, and commensal and laboratory-adapted E. coli. However, the bovine STEC/ETEC strain formed a remote cluster with two STECs of bovine origin. All three STEC/ETEC strains harbored several other virulence genes, apart from stx and est, and lacked ETEC colonization factors. Two STEC/ETEC strains produced both toxins and one strain Stx only.

Conclusions

This study shows that pathogroup-associated virulence genes of different E. coli can co-exist in strains originating from different phylogenetic lineages. The possibility of virulence genes to be associated with several E. coli pathogroups should be taken into account in strain typing and in epidemiological surveillance. Development of novel hybrid E. coli strains may cause a new public health risk, which challenges the traditional diagnostics of E. coli infections.     

DNA Microarray-Based Identification of Serogroups and Virulence Gene Patterns of Escherichia coli Isolates Associated with Porcine Postweaning Diarrhea and Edema Disease

Escherichia coli strains causing postweaning diarrhea (PWD) and edema disease (ED) in pigs are limited to a number of serogroups, with O8, O45, O138, O139, O141, O147, O149, and O157 being the most commonly reported worldwide. In this study, a DNA microarray based on the O-antigen-specific genes of all 8 E. coli serogroups, as well as 11 genes encoding adhesion factors and exotoxins associated with PWD and ED, was developed for the identification of related serogroups and virulence gene patterns. The microarray method was tested against 186 E. coli and Shigella O-serogroup reference strains, 13 E. coli reference strains for virulence markers, 43 E. coli clinical isolates, and 12 strains of other bacterial species and shown to be highly specific with reproducible results. The detection sensitivity was 0.1 ng of genomic DNA or 10³ CFU per 0.3 g of porcine feces in mock samples. Seventeen porcine feces samples from local hoggeries were examined using the microarray, and the result for one sample was verified by the conventional serotyping methods. This microarray can be readily used to screen for the presence of PWD- and ED-associated E. coli in porcine feces samples.     

The Evolutionary History of <Emphasis Type="Italic">Shigella</Emphasis> and Enteroinvasive <Emphasis Type="Italic">Escherichia coli</Emphasis> Revised

In Shigella and enteroinvasive Escherichia coli (EIEC), the etiologic agents of shigellosis in humans, the determinants responsible for entry of bacteria into and dissemination within epithelial cells are encoded by a virulence plasmid. To understand the evolution of the association between the virulence plasmid and the chromosome, we performed a phylogenetic analysis using the sequences of four chromosomal genes (trpA, trpB, pabB, and putP) and three virulence plasmid genes (ipaB, ipaD, and icsA) of a collection of 51 Shigella and EIEC strains. The phylogenetic tree derived from chromosomal genes showed a typical star phylogeny, indicating a fast diversification of Shigella and EIEC groups. Phylogenetic groups obtained from the chromosomal and plasmidic genes were similar, suggesting that the virulence plasmid and the chromosome share similar evolutionary histories. The few incongruences between the trees could be attributed to exchanges of fragments of different plasmids and not to the transfer of an entire plasmid. This indicates that the virulence plasmid was not transferred between the different Shigella and EIEC groups. These data support a model of evolution in which the acquisition of the virulence plasmid in an ancestral E. coli strain preceded the diversification by radiation of all Shigella and EIEC groups, which led to highly diversified but highly specialized pathogenic groups.     

An Investigation of the Diversity of Strains of Enteroaggregative Escherichia coli Isolated from Cases Associated with a Large Multi-Pathogen Foodborne Outbreak in the UK

Following a large outbreak of foodborne gastrointestinal (GI) disease, a multiplex PCR approach was used retrospectively to investigate faecal specimens from 88 of the 413 reported cases. Gene targets from a range of bacterial GI pathogens were detected, including Salmonella species, Shigella species and Shiga toxin-producing Escherichia coli, with the majority (75%) of faecal specimens being PCR positive for aggR associated with the Enteroaggregative E. coli (EAEC) group. The 20 isolates of EAEC recovered from the outbreak specimens exhibited a range of serotypes, the most frequent being O104:H4 and O131:H27. None of the EAEC isolates had the Shiga toxin (stx) genes. Multilocus sequence typing and single nucleotide polymorphism analysis of the core genome confirmed the diverse phylogeny of the strains. The analysis also revealed a close phylogenetic relationship between the EAEC O104:H4 strains in this outbreak and the strain of E. coli O104:H4 associated with a large outbreak of haemolytic ureamic syndrome in Germany in 2011. Further analysis of the EAEC plasmids, encoding the key enteroaggregative virulence genes, showed diversity with respect to FIB/FII type, gene content and genomic architecture. Known EAEC virulence genes, such as aggR, aat and aap, were present in all but one of the strains. A variety of fimbrial genes were observed, including genes encoding all five known fimbrial types, AAF/1 to AAF/V. The AAI operon was present in its entirety in 15 of the EAEC strains, absent in three and present, but incomplete, in two isolates. EAEC is known to be a diverse pathotype and this study demonstrates that a high level of diversity in strains recovered from cases associated with a single outbreak. Although the EAEC in this study did not carry the stx genes, this outbreak provides further evidence of the pathogenic potential of the EAEC O104:H4 serotype.     

Comparative genomics of enterohemorrhagic Escherichia coli O145:H28 demonstrates a common evolutionary lineage with Escherichia coli O157:H7

Background

Although serotype O157:H7 is the predominant enterohemorrhagic Escherichia coli (EHEC), outbreaks of non-O157 EHEC that cause severe foodborne illness, including hemolytic uremic syndrome have increased worldwide. In fact, non-O157 serotypes are now estimated to cause over half of all the Shiga toxin-producing Escherichia coli (STEC) cases, and outbreaks of non-O157 EHEC infections are frequently associated with serotypes O26, O45, O103, O111, O121, and O145. Currently, there are no complete genomes for O145 in public databases.

Results

We determined the complete genome sequences of two O145 strains (EcO145), one linked to a US lettuce-associated outbreak (RM13514) and one to a Belgium ice-cream-associated outbreak (RM13516). Both strains contain one chromosome and two large plasmids, with genome sizes of 5,737,294 bp for RM13514 and 5,559,008 bp for RM13516. Comparative analysis of the two EcO145 genomes revealed a large core (5,173 genes) and a considerable amount of strain-specific genes. Additionally, the two EcO145 genomes display distinct chromosomal architecture, virulence gene profile, phylogenetic origin of Stx2a prophage, and methylation profile (methylome). Comparative analysis of EcO145 genomes to other completely sequenced STEC and other E. coli and Shigella genomes revealed that, unlike any other known non-O157 EHEC strain, EcO145 ascended from a common lineage with EcO157/EcO55. This evolutionary relationship was further supported by the pangenome analysis of the 10 EHEC str ains. Of the 4,192 EHEC core genes, EcO145 shares more genes with EcO157 than with the any other non-O157 EHEC strains.

Conclusions

Our data provide evidence that EcO145 and EcO157 evolved from a common lineage, but ultimately each serotype evolves via a lineage-independent nature to EHEC by acquisition of the core set of EHEC virulence factors, including the genes encoding Shiga toxin and the large virulence plasmid. The large variation between the two EcO145 genomes suggests a distinctive evolutionary path between the two outbreak strains. The distinct methylome between the two EcO145 strains is likely due to the presence of a BsuBI/PstI methyltransferase gene cassette in the Stx2a prophage of the strain RM13514, suggesting a role of horizontal gene transfer-mediated epigenetic alteration in the evolution of individual EHEC strains.     

An integrated approach for finding overlooked genes in Shigella

Background

The completion of numerous genome sequences introduced an era of whole-genome study. However, many genes are missed during genome annotation, including small RNAs (sRNAs) and small open reading frames (sORFs). In order to improve genome annotation, we aimed to identify novel sRNAs and sORFs in Shigella, the principal etiologic agents of bacillary dysentery.

Methodology/Principal Findings

We identified 64 sRNAs in Shigella, which were experimentally validated in other bacteria based on sequence conservation. We employed computer-based and tiling array-based methods to search for sRNAs, followed by RT-PCR and northern blots, to identify nine sRNAs in Shigella flexneri strain 301 (Sf301) and 256 regions containing possible sRNA genes. We found 29 candidate sORFs using bioinformatic prediction, array hybridization and RT-PCR verification. We experimentally validated 557 (57.9%) DOOR operon predictions in the chromosomes of Sf301 and 46 (76.7%) in virulence plasmid.We found 40 additional co-expressed gene pairs that were not predicted by DOOR.

Conclusions/Significance

We provide an updated and comprehensive annotation of the Shigella genome. Our study increased the expected numbers of sORFs and sRNAs, which will impact on future functional genomics and proteomics studies. Our method can be used for large scale reannotation of sRNAs and sORFs in any microbe with a known genome sequence.