The evolutionary history of Shigella and enteroinvasive Escherichia coli revised

A novel biochemistry approach combined with MALDI-TOF MS to discriminate Escherichia coli and Shigella species

Escherichia coli and Shigella spp. are closely related, making it crucial to accurately identify them for disease control and prevention. In this study, we utilized MALDI-TOF MS to identify characteristic peaks of decarboxylation products of lysine and ornithine to distinguish between E. coli and Shigella spp. Our findings indicate that the peak at m/z 103.12 ± 0.1 of the product cadaverine from lysine decarboxylase is unique to E. coli, while all Shigella species lack the m/z 103.12 ± 0.1 peak. However, S. sonnei and S. boydii serotype C13 exhibit a specific peak at m/z 89.10 ± 0.1, which is the product of putrescine from ornithine decarboxylase. We were able to correctly identify 97.06% (132 of 136) of E. coli and Shigella isolates and 100% (8 of 8) of S. sonnei isolates using this biochemical-based MALDI-TOF MS detection system. This technology is advantageous for its high-throughput, high quality, and ease of operation, and is of significant value for the diagnosis of E. coli and Shigella-related diseases.

A Tale about Shigella: Evolution, Plasmid, and Virulence

Shigella spp. cause hundreds of millions of intestinal infections each year. They target the mucosa of the human colon and are an important model of intracellular bacterial pathogenesis. Shigella is a pathovar of Escherichia coli that is characterized by the presence of a large invasion plasmid, pINV, which encodes the characteristic type III secretion system and icsA used for cytosol invasion and cell-to-cell spread, respectively. First, we review recent advances in the genetic aspects of Shigella, shedding light on its evolutionary history within the E. coli lineage and its relationship to the acquisition of pINV. We then discuss recent insights into the processes that allow for the maintenance of pINV. Finally, we describe the role of the transcription activators VirF, VirB, and MxiE in the major virulence gene regulatory cascades that control the expression of the type III secretion system and icsA. This provides an opportunity to examine the interplay between these pINV-encoded transcriptional activators and numerous chromosome-encoded factors that modulate their activity. Finally, we discuss novel chromosomal genes icaR, icaT, and yccE that are regulated by MxiE. This review emphasizes the notion that Shigella and E. coli have walked the fine line between commensalism and pathogenesis for much of their history.

The Diversity of Escherichia coli Pathotypes and Vaccination Strategies against This Versatile Bacterial Pathogen

Escherichia coli (E. coli) is a gram-negative bacillus and resident of the normal intestinal microbiota. However, some E. coli strains can cause diseases in humans, other mammals and birds ranging from intestinal infections, for example, diarrhea and dysentery, to extraintestinal infections, such as urinary tract infections, respiratory tract infections, meningitis, and sepsis. In terms of morbidity and mortality, pathogenic E. coli has a great impact on public health, with an economic cost of several billion dollars annually worldwide. Antibiotics are not usually used as first-line treatment for diarrheal illness caused by E. coli and in the case of bloody diarrhea, antibiotics are avoided due to the increased risk of hemolytic uremic syndrome. On the other hand, extraintestinal infections are treated with various antibiotics depending on the site of infection and susceptibility testing. Several alarming papers concerning the rising antibiotic resistance rates in E. coli strains have been published. The silent pandemic of multidrug-resistant bacteria including pathogenic E. coli that have become more difficult to treat favored prophylactic approaches such as E. coli vaccines. This review provides an overview of the pathogenesis of different pathotypes of E. coli, the virulence factors involved and updates on the major aspects of vaccine development against different E. coli pathotypes.

The microbial ecology of Escherichia coli in the vertebrate gut

Escherichia coli has a rich history as biology's 'rock star', driving advances across many fields. In the wild, E. coli resides innocuously in the gut of humans and animals but is also a versatile pathogen commonly associated with intestinal and extraintestinal infections and antimicrobial resistance-including large foodborne outbreaks such as the one that swept across Europe in 2011, killing 54 individuals and causing approximately 4000 infections and 900 cases of haemolytic uraemic syndrome. Given that most E. coli are harmless gut colonizers, an important ecological question plaguing microbiologists is what makes E. coli an occasionally devastating pathogen? To address this question requires an enhanced understanding of the ecology of the organism as a commensal. Here, we review how our knowledge of the ecology and within-host diversity of this organism in the vertebrate gut has progressed in the 137 years since E. coli was first described. We also review current approaches to the study of within-host bacterial diversity. In closing, we discuss some of the outstanding questions yet to be addressed and prospects for future research.

icaR and icaT Are Ancient Chromosome Genes Encoding Substrates of the Type III Secretion Apparatus in Shigella flexneri

Shigella is an Escherichia coli pathovar that colonizes the cytosol of mucosal cells in the human large intestine. To do this, Shigella uses a Type III Secretion Apparatus (T3SA) to translocate several proteins into host cells. The T3SA and its substrates are encoded by genes of the virulence plasmid pINV or by chromosomal genes derived thereof. We recently discovered two chromosomal genes, which seem unrelated to pINV, although they are activated by MxiE and IpgC similarly to some of the canonical substrates of the T3SA. Here, we showed that the production of the corresponding proteins depended on the conservation of a MxiE box in their cognate promoters. Furthermore, both proteins were secreted by the T3SA in a chaperone-independent manner through the recognition of their respective amino-terminal secretion signal. Based on these observations, we named these new genes icaR and icaT, which stand for invasion chromosome antigen with homology for a transcriptional regulator and a transposase, respectively. icaR and icaT have orthologs in commensal and pathogenic E. coli strains belonging mainly to phylogroups A, B1, D and E. Finally, we demonstrated that icaR and icaT orthologs could be activated by the coproduction of IpgC and MxiE in strains MG1655 K-12 (phylogroup A) and O157:H7 ATCC 43888 (phylogroup E). In contrast, the coproduction of EivF and YgeG, which are homologs of MxiE and IpgC in the E. coli T3SS 2 (ETT2), failed to activate icaR and icaT.

IMPORTANCEicaR and icaT are the latest members of the MxiE regulon discovered in the chromosome. The proteins IcaR and IcaT, albeit produced in small amounts, are nonetheless secreted by the T3SA comparably to canonical substrates. The high occurrence of icaR and icaT in phylogroups A, B1, D, and E coupled with their widespread absence in their B2 counterparts agree with the consensus E. coli phylogeny. The widespread conservation of the MxiE box among icaR and icaT orthologs supports the notion that both genes had already undergone coevolution with transcriptional activators ipgC and mxiE- harbored in pINV or a relative- in the last common ancestor of Shigella and of E. coli from phylogroups A, B1, D, and E. The possibility that icaR and icaT may contribute to Shigella pathogenesis cannot be excluded, although some of their characteristics suggest they are fossil genes.

Virulence Factors of Enteric Pathogenic Escherichia coli: A Review

Escherichia coli are remarkably versatile microorganisms and important members of the normal intestinal microbiota of humans and animals. This harmless commensal organism can acquire a mixture of comprehensive mobile genetic elements that contain genes encoding virulence factors, becoming an emerging human pathogen capable of causing a broad spectrum of intestinal and extraintestinal diseases. Nine definite enteric E. coli pathotypes have been well characterized, causing diseases ranging from various gastrointestinal disorders to urinary tract infections. These pathotypes employ many virulence factors and effectors subverting the functions of host cells to mediate their virulence and pathogenesis. This review summarizes new developments in our understanding of diverse virulence factors associated with encoding genes used by different pathotypes of enteric pathogenic E. coli to cause intestinal and extraintestinal diseases in humans.

Occurrence and distribution of fecal indicators and pathogenic bacteria in seawater and Perna perna mussel in the Gulf of Annaba (Southern Mediterranean)

The identification of fecal contamination in coastal marine ecosystems is one of the main requirements for evaluation of potential risks to human health. The objective of this study was to investigate the occurrence and distribution of fecal indicators and pathogenic bacteria in seawaters and mussels collected monthly during a period of 1 year from four different sites in Northeastern Algeria (sites S1 to S4), through biochemical and molecular analyses. Our research is the first to use molecular analysis to unambiguously identify the potentially pathogenic bacteria present in Algerian Perna perna mussels. The obtained results revealed that the levels of fecal indicator bacteria (FIB) from both P. perna and seawater samples largely exceeded the permissible limits at S2 and S3. This is mainly related to their location close to industrial and coastal activity zones, which contain a mixture of urban, agricultural, and industrial pollutants. Besides, P. perna collected from all sites were severalfold more contaminated by FIB than seawater samples, primarily during the warm season of the study period. Biochemical and molecular analyses showed that isolated bacteria from both seawater and mussels were mainly potentially pathogenic species such as E. coli, Salmonella spp., Staphylococcus spp., Klebsiella spp., Pseudomonas spp., and Proteus spp.

Effect of sterilization on the canine vaginal microbiota: a pilot study

BACKGROUND

Surgical sterilization is the most effective method of contraception for dogs. It also prevents pyometra and reduces the risk of mammary tumour development. However, this procedure also has negative effects, such as urinary incontinence. Steroid hormone deprivation following gonadectomy could also affect canine vaginal mucosa conditions and the microbial community colonizing the vaginal tract. This hypothesis was tested by comparing the vaginal cytology and microbial community of two groups of bitches, including 11 in anoestrus and 10 sterilized bitches (post-pubertal sterilization in the last 4 years). Bacteria were identified through metataxonomic analysis, amplifying the V3-V4 regions of 16S rRNA gene, and culturing methods.

RESULTS

Vaginal mucosa cytology was suggestive of dystrophic conditions in sterilized bitches, whereas a typical anoestrus pattern with parabasal and intermediate cells was appreciable in anoestrous animals. Metataxonomic analysis revealed large inter-individual variability. Salmonella, Mycoplasma and Staphylococcus were present in moderate quantities in almost all the samples in both groups. Mollicutes (class level) and Tenericutes (phylum level) were commonly present in moderate quantities in anoestrus samples, whereas these microbes were present at high levels in a single sample from the sterilized group. Based on culturing, a higher number of different species were isolated from the anoestrous bitches, and Mycoplasma canis was exclusively identified in an anoestrous bitch. Staphylococcus spp. was the most frequently isolated genus in both groups, followed by Streptococcus spp., and, among gram-negative bacteria, Escherichia spp. and Haemophilus spp. A comparison of the numbers of the most frequently isolated genera of bacteria from vaginal cultures of bitches revealed that Pasteurella and Proteus were the most frequently identified in sterilized animals based on metataxonomic analysis (p-value = 0.0497 and 0.0382, respectively), whereas Streptococcus was significantly and most frequently isolated from anoestrous bitches using culture methods (p value = 0.0436).

CONCLUSIONS

In this preliminary investigation, no global patterns of the vaginal bacteria community were noted that characterized the condition of the bitches; however, cytology suggested local modifications. Sterilization after puberty caused minimal alterations in the vaginal microbial community of bitches within 4 years after surgery.

Shifting national surveillance of Shigella infections toward geno-serotyping by the development of a tailored Luminex assay and NGS workflow

The phylogenetically closely related Shigella species and enteroinvasive Escherichia coli (EIEC) are responsible for millions of episodes of bacterial dysenteriae worldwide. Given its distinct epidemiology and public health relevance, only Shigellae are subject to mandatory reporting and follow‐up by public health authorities. However, many clinical laboratories struggle to differentiate non‐EIEC, EIEC, and Shigella in their current workflows, leading to inaccuracies in surveillance and rising numbers of misidentified E. coli samples at the National Reference Centre (NRC). In this paper, we describe two novel tools to enhance Shigella surveillance. First, we developed a low‐cost Luminex‐based multiplex assay combining five genetic markers for species identification with 11 markers for serotype prediction for S. sonnei and S. flexneri isolates. Using a test panel of 254 clinical samples, this assay has a sensitivity of 100% in differentiation of EIEC/Shigella pathotype from non‐EIEC strains, and 68.7% success rate in distinction of Shigella and EIEC. A novel, and particularly successful marker was a Shigella‐specific deletion in the spermidine acetyltransferase gene speG, reflecting its metabolic decay. For Shigella serotype prediction, the multiplex assay scored a sensitivity and specificity of 96.6% and 98.4%, respectively. All discrepancies were analyzed with whole‐genome sequencing and shown to be related to causative mutations (stop codons, indels, and promoter mutations) in glycosyltransferase genes. This observation spurred the development of an in silico workflow which extracts the Shigella serotype from Next‐Generation Sequencing (NGS) data, taking into account gene functionality. Both tools will be implemented in the workflow of the NRC, and will play a major role in the shift from phenotypic to genotyping‐based surveillance of shigellosis in Belgium.

Virulence-related genes are associated with clinical and nutritional outcomes of Shigella/Enteroinvasive Escherichia coli pathotype infection in children from Brazilian semiarid region: A community case-control study

Shigella/Enteroinvasive Escherichia coli (EIEC) pathotype is a major enteropathogen associated with diarrhea and malnutrition in children from developing countries. This study aimed to correlate Shigella/EIEC virulence-related genes (VRGs) with clinical symptoms, nutritional status and coenteropathogens in children from the Brazilian semiarid region. We designed a case-control study of community diarrhea in six cities of the Brazil semiarid region with 1200 children aging 2-36 months. Standardized questionnaire was applied for collecting sociodemographic, nutritional status and clinical information of the children. DNA samples were extracted from stools and diagnosed for Shigella/EIEC using PCR-based approaches. Positive samples were tested for 28 VRGs using four multiplex PCRs. Intestinal inflammation was determined by measuring fecal myeloperoxidase (MPO). Shigella/EIEC pathotype was detected in 5% of the children and was significantly associated with diarrhea. The genes sen (encoding Shigella enterotoxin 2), ipgB2, ipgB1 (both encoding type 3 secretion system-T3SS effectors that modulate actin filament), and ospF (encoding a T3SS effector involved in suppression of host responses) were further associated with diarrhea in Shigella/EIEC positive children. Among children presenting diarrhea, virA gene (encoding a T3SS effector that promotes microtubule destabilization) was associated with fever, while virB (encoding a major transcriptional activator) was associated with low height-for-age z-score. In addition, these VRGs were associated with increased fecal MPO, and coinfection with Salmonella spp. was associated with increased abdominal pain. These data reinforce the impact of Shigella/EIEC on diarrhea in children from Brazilian semiarid region and highlighted the contributions of specific virulence genes for its pathobiology.

Taxonomy annotation and guide tree errors in 16S rRNA databases

Sequencing of the 16S ribosomal RNA (rRNA) gene is widely used to survey microbial communities. Specialized 16S rRNA databases have been developed to support this approach including Greengenes, RDP and SILVA. Most taxonomy annotations in these databases are predictions from sequence rather than authoritative assignments based on studies of type strains or isolates. In this work, I investigated the taxonomy annotations and guide trees provided by these databases. Using a blinded test, I estimated that the annotation error rate of the RDP database is ∼10%. The branching orders of the Greengenes and SILVA guide trees were found to disagree at comparable rates with each other and with taxonomy annotations according to the training set (authoritative reference) provided by RDP, indicating that the trees have comparable quality. Pervasive conflicts between tree branching order and type strain taxonomies strongly suggest that the guide trees are unreliable guides to phylogeny. I found 249,490 identical sequences with conflicting annotations in SILVA v128 and Greengenes v13.5 at ranks up to phylum (7,804 conflicts), indicating that the annotation error rate in these databases is ∼17%.

Escherichia coli: an old friend with new tidings

Escherichia coli is one of the most-studied microorganisms worldwide but its characteristics are continually changing. Extraintestinal E. coli infections, such as urinary tract infections and neonatal sepsis, represent a huge public health problem. They are caused mainly by specialized extraintestinal pathogenic E. coli (ExPEC) strains that can innocuously colonize human hosts but can also cause disease upon entering a normally sterile body site. The virulence capability of such strains is determined by a combination of distinctive accessory traits, called virulence factors, in conjunction with their distinctive phylogenetic background. It is conceivable that by developing interventions against the most successful ExPEC lineages or their key virulence/colonization factors the associated burden of disease and health care costs could foreseeably be reduced in the future. On the other hand, one important problem worldwide is the increase of antimicrobial resistance shown by bacteria. As underscored in the last WHO global report, within a wide range of infectious agents including E. coli, antimicrobial resistance has reached an extremely worrisome situation that ‘threatens the achievements of modern medicine’. In the present review, an update of the knowledge about the pathogenicity, antimicrobial resistance and clinical aspects of this ‘old friend’ was presented.

The Intriguing Evolutionary Journey of Enteroinvasive E. coli (EIEC) toward Pathogenicity

Among the intestinal pathogenic Escherichia coli, enteroinvasive E. coli (EIEC) are a group of intracellular pathogens able to enter epithelial cells of colon, multiplicate within them, and move between adjacent cells with a mechanism similar to Shigella, the ethiological agent of bacillary dysentery. Despite EIEC belong to the same pathotype of Shigella, they neither have the full set of traits that define Shigella nor have undergone the extensive gene decay observed in Shigella. Molecular analysis confirms that EIEC are widely distributed among E. coli phylogenetic groups and correspond to bioserotypes found in many E. coli serogroups. Like Shigella, also in EIEC the critical event toward a pathogenic life-style consisted in the acquisition by horizontal gene transfer of a large F-type plasmid (pINV) containing the genes required for invasion, intracellular survival, and spreading through the intestinal mucosa. In Shigella, the ample gain in virulence determinants has been counteracted by a substantial loss of functions that, although important for the survival in the environment, are redundant or deleterious for the life inside the host. The pathoadaptation process that has led Shigella to modify its metabolic profile and increase its pathogenic potential is still in infancy in EIEC, although maintenance of some features typical of E. coli might favor their emerging relevance as intestinal pathogens worldwide, as documented by recent outbreaks in industrialized countries. In this review, we will discuss the evolution of EIEC toward Shigella-like invasive forms going through the epidemiology, including the emergence of new virulent strains, their genome organization, and the complex interactions they establish with the host.

Accurate differentiation of Escherichia coli and Shigella serogroups: challenges and strategies

Shigella spp. and Escherichia coli are closely related; both belong to the family Enterobacteriaceae. Phenotypically, Shigella spp. and E. coli share many common characteristics, yet they have separate entities in epidemiology and clinical disease, which poses a diagnostic challenge. We collated information for the best possible approach to differentiate clinically relevant E. coli from Shigella spp. We found that a molecular approach is required for confirmation. High discriminatory potential is seen with whole genome sequencing analysed for k-mers and single nucleotide polymorphism. Among these, identification using single nucleotide polymorphism is easy to perform and analyse, and it thus appears more promising. Among the nonmolecular methods, matrix-assisted desorption ionization–time of flight mass spectrometry may be applicable when data analysis is assisted with advanced analytic tools.

Diversity of the auxotrophic requirements in natural isolates of Escherichia coli

Isolates of Escherichia coli, except Shigella, are generally prototrophic; they do not require any growth factors to grow in mineral medium. However, a nicotinic acid requirement is common among B2 phylogroup STc95 O18 E. coli clone strains. Nicotinic acid is a precursor of nicotinamide adenine dinucleotide (NAD), an essential molecule that plays central role in cellular metabolism. The defect in NAD synthesis of these strains is due to alterations in de novo biosynthesis pathway nadB gene. Here, by studying growth on minimal medium with glycolytic (glucose) or gluconeogenic (pyruvate or succinate) substrates as the carbon supply in a large panel of E. coli natural isolates representative of the species diversity, we identify an absolute nicotinic acid requirement in non-STc95 strains due in one case to a nadA inactivation. The growth on glucose medium of some extraintestinal pathogenic E. coli strains belonging to various non-O18 B2 phylogroup STc95 clones is restored either by aspartate or nicotinate, demonstrating that the nicotinic acid requirement can also be due to an intracellular aspartate depletion. The auxotrophic requirements depend on the carbon source available in the environment. Moreover, some strains prototrophic in glucose medium become auxotrophic in succinate medium, and conversely, some strainsauxotrophic in glucose medium become prototrophic in succinate medium. Finally, a partial depletion of intracellular aspartate can be observed in some prototrophic strains belonging to various phylogroups. The observed more or less significant depletion according to isolates may be due to differences in tricarboxylic acid cycle enzyme activities. These metabolic defects could be involved in the adaptation of E. coli to its various niches.

Phylogenetic tree construction using trinucleotide usage profile (TUP)

BACKGROUND

It has been a challenging task to build a genome-wide phylogenetic tree for a large group of species containing a large number of genes with long nucleotides sequences. The most popular method, called feature frequency profile (FFP-k), finds the frequency distribution for all words of certain length k over the whole genome sequence using (overlapping) windows of the same length. For a satisfactory result, the recommended word length (k) ranges from 6 to 15 and it may not be a multiple of 3 (codon length). The total number of possible words needed for FFP-k can range from 4⁶=4096 to 4¹⁵.

RESULTS

We propose a simple improvement over the popular FFP method using only a typical word length of 3. A new method, called Trinucleotide Usage Profile (TUP), is proposed based only on the (relative) frequency distribution using non-overlapping windows of length 3. The total number of possible words needed for TUP is 4³=64, which is much less than the total count for the recommended optimal "resolution" for FFP. To build a phylogenetic tree, we propose first representing each of the species by a TUP vector and then using an appropriate distance measure between pairs of the TUP vectors for the tree construction. In particular, we propose summarizing a DNA sequence by a matrix of three rows corresponding to three reading frames, recording the frequency distribution of the non-overlapping words of length 3 in each of the reading frame. We also provide a numerical measure for comparing trees constructed with various methods.

CONCLUSIONS

Compared to the FFP method, our empirical study showed that the proposed TUP method is more capable of building phylogenetic trees with a stronger biological support. We further provide some justifications on this from the information theory viewpoint. Unlike the FFP method, the TUP method takes the advantage that the starting of the first reading frame is (usually) known. Without this information, the FFP method could only rely on the frequency distribution of overlapping words, which is the average (or mixture) of the frequency distributions of three possible reading frames. Consequently, we show (from the entropy viewpoint) that the FFP procedure could dilute important gene information and therefore provides less accurate classification.

Investigating the Relatedness of Enteroinvasive Escherichia coli to Other E. coli and Shigella Isolates by Using Comparative Genomics

Enteroinvasive Escherichia coli (EIEC) is a unique pathovar that has a pathogenic mechanism nearly indistinguishable from that of Shigella species. In contrast to isolates of the four Shigella species, which are widespread and can be frequent causes of human illness, EIEC causes far fewer reported illnesses each year. In this study, we analyzed the genome sequences of 20 EIEC isolates, including 14 first described in this study. Phylogenomic analysis of the EIEC genomes demonstrated that 17 of the isolates are present in three distinct lineages that contained only EIEC genomes, compared to reference genomes from each of the E. coli pathovars and Shigella species. Comparative genomic analysis identified genes that were unique to each of the three identified EIEC lineages. While many of the EIEC lineage-specific genes have unknown functions, those with predicted functions included a colicin and putative proteins involved in transcriptional regulation or carbohydrate metabolism. In silico detection of the Shigella virulence plasmid (pINV), which is essential for the invasion of host cells, demonstrated that a form of pINV was present in nearly all EIEC genomes, but the Mxi-Spa-Ipa region of the plasmid that encodes the invasion-associated proteins was absent from several of the EIEC isolates. The comparative genomic findings in this study support the hypothesis that multiple EIEC lineages have evolved independently from multiple distinct lineages of E. coli via the acquisition of the Shigella virulence plasmid and, in some cases, the Shigella pathogenicity islands.

Phylogenetic Analyses of Shigella and Enteroinvasive Escherichia coli for the Identification of Molecular Epidemiological Markers: Whole-Genome Comparative Analysis Does Not Support Distinct Genera Designation

As a leading cause of bacterial dysentery, Shigella represents a significant threat to public health and food safety. Related, but often overlooked, enteroinvasive Escherichia coli (EIEC) can also cause dysentery. Current typing methods have limited ability to identify and differentiate between these pathogens despite the need for rapid and accurate identification of pathogens for clinical treatment and outbreak response. We present a comprehensive phylogeny of Shigella and EIEC using whole genome sequencing of 169 samples, constituting unparalleled strain diversity, and observe a lack of monophyly between Shigella and EIEC and among Shigella taxonomic groups. The evolutionary relationships in the phylogeny are supported by analyses of population structure and hierarchical clustering patterns of translated gene homolog abundance. Lastly, we identified a panel of 254 single nucleotide polymorphism (SNP) markers specific to each phylogenetic cluster for more accurate identification of Shigella and EIEC. Our findings show that Shigella and EIEC are not distinct evolutionary groups within the E. coli genus and, thus, EIEC as a group is not the ancestor to Shigella. The multiple analyses presented provide evidence for reconsidering the taxonomic placement of Shigella. The SNP markers offer more discriminatory power to molecular epidemiological typing methods involving these bacterial pathogens.

Virulence Gene Regulation in Shigella

Shigella species are the causative agents of bacillary dysentery in humans, an invasive disease in which the bacteria enter the cells of the epithelial layer of the large intestine, causing extensive tissue damage and inflammation. They rely on a plasmid-encoded type III secretion system (TTSS) to cause disease; this system and its regulation have been investigated intensively at the molecular level for decades. The lessons learned have not only deepened our knowledge of Shigella biology but also informed in important ways our understanding of the mechanisms used by other pathogenic bacteria to cause disease and to control virulence gene expression. In addition, the Shigella story has played a central role in the development of our appreciation of the contribution of horizontal DNA transfer to pathogen evolution.A 30-kilobase-pair "Entry Region" of the 230-kb virulence plasmid lies at the heart of the Shigella pathogenesis system. Here are located the virB and mxiE regulatory genes and most of the structural genes involved in the expression of the TTSS and its effector proteins. Expression of the virulence genes occurs in response to an array of environmental signals, including temperature, osmolarity, and pH.At the top of the regulatory hierarchy and lying on the plasmid outside the Entry Region isvirF, encoding an AraC-like transcription factor.Virulence gene expression is also controlled by chromosomal genes,such as those encoding the nucleoid-associated proteins H-NS, IHF, and Fis, the two-component regulators OmpR/EnvZ and CpxR/CpxA, the anaerobic regulator Fnr, the iron-responsive regulator Fur, and the topoisomerases of the cell that modulate DNA supercoiling. Small regulatory RNAs,the RNA chaperone Hfq,and translational modulation also affect the expression of the virulence phenotypetranscriptionally and/orposttranscriptionally.

Distinct mutations led to inactivation of type 1 fimbriae expression in Shigella spp

Shigella spp. are responsible for bacillary dysentery in humans. The acquisition or the modification of the virulence plasmid encoding factors promoting entry of bacteria into and dissemination within epithelial cells was a critical step in the evolution of these bacteria from their Escherichia coli ancestor(s). Incorporation of genomic islands (GI) and gene inactivation also shaped interactions between these pathogens and their human host. Sequence analysis of the GI inserted next to the leuX tRNA gene in S. boydii, S. dysenteriae, S. flexneri, S. sonnei and enteroinvasive E. coli (EIEC) suggests that this region initially carried the fec, yjhATS and fim gene clusters. The fim cluster encoding type I fimbriae is systematically inactivated in both reference strains and clinical isolates and distinct mutations are responsible for this inactivation in at least three phylogenetic groups. To investigate consequences of the presence of fimbriae on the outcome of the interaction of Shigella with host cells, we used a S. flexneri strain harboring a plasmid encoding the E. coli fim operon. Production of fimbriae by this recombinant strain increased the ability of bacteria to adhere to and enter into epithelial cells and had no effect on their ability to disseminate from cell to cell. The observations that production of type I fimbriae increases invasion of epithelial cells and that independent mutations abolish fimbriae production in Shigella suggest that these mutations correspond to pathoadaptive events.

Escherichia fergusonii

Escherichia fergusonii was introduced in the genus Escherichia almost 65 years later than Escherichia coli after which the genus was named. From then (1985) onwards mainly case reports on E. fergusonii associated with disease in individuals of veterinary or human origin have been reported and only very few more extensive studies became available. This has resulted in very fragmented knowledge on this organism. The aim of this manuscript is to give an overview of what is known on E. fergusonii today and to stimulate more research on this organism so that better insight can be obtained in the role that E. fergusonii plays in human and animal infections.

Water quality monitoring and risk assessment by simultaneous multipathogen quantification

Water quality monitoring and microbial risk assessment are important to ensure safe water for drinking, recreational, and agricultural purposes. In this study, we applied a microfluidic quantitative PCR (MFQPCR) approach to simultaneously quantify multiple waterborne pathogens in a natural freshwater lake in Hokkaido, Japan, from April to November, 2012. Tens of thousands of geese stopped over at this lake during their migration in spring and fall. Because lake water is used for irrigation of the surrounding agricultural area, we assessed infection risks through irrigation water usage based on pathogen concentrations directly measured by MFQPCR. We detected various pathogens in the lake water, particularly during the bird migration seasons, suggesting that migratory birds were the main source of the pathogens. However, neither counts of geese nor fecal indicator bacteria were good predictors of pathogen concentrations. On the basis of quantitative microbial risk assessment, concentrations of Campylobacter jejuni and Shigella spp. in water samples were above the concentrations that can potentially cause 10(-4) infections per person per year when water is used to grow fresh vegetables. These results suggest that direct and simultaneous multipathogen quantification can provide more reliable and comprehensive information for risk assessment than the current fecal indicator-based approach.

A DNA microarray for the versatile diagnosis of infectious diarrhea

Several bacteria, viruses, and parasites cause diarrhea as coinfecting pathogens. We designed a DNA microarray comprising 60‐bp probes spotted 194 times for the multiplex detection of 33 enteropathogenic bacteria and seven enteropathogenic viruses, and the archaeon Methanobrevibacter smithii was used as an internal positive control. Nine pathogen‐free stool specimens were used as negative controls. One of these control specimens was further spiked with Salmonella enterica as a positive control. The microarray was then tested with 40 pathological stool specimens, comprising S. enterica (n = 30), Campylobacter jejuni (n = 4), pathogenic Escherichia coli (n = 2), and adenovirus (n = 4). M. smithii was detected in 47/49 (95.9%) specimens, no pathogen was detected in negative controls and S. enterica was identified in the S. enterica‐spiked positive control. The overall specificity was 100% and the overall sensitivity was 97.5% because one S. enterica sample was missed by the microarray. The multiplexed detection of C. jejuni spiked into an adenovirus‐positive stool sample gave positive results, with fluorescence values of 14.3 and 9.1, respectively. These data indicate that using the protocol developed in this article, the DNA array allows for the multiplexed detection of some enteropathogens in stool samples.

Simultaneous quantification of multiple food- and waterborne pathogens by use of microfluidic quantitative PCR

The direct quantification of multiple pathogens has been desired for diagnostic and public health purposes for a long time. In this study, we applied microfluidic quantitative PCR (qPCR) technology to the simultaneous detection and quantification of multiple food- and waterborne pathogens. In this system, multiple singleplex qPCR assays were run under identical detection conditions in nanoliter-volume chambers that are present in high densities on a chip. First, we developed 18 TaqMan qPCR assays that could be run in the same PCR conditions by using prevalidated TaqMan probes. Specific and sensitive quantification was achieved by using these qPCR assays. With the addition of two previously validated TaqMan qPCR assays, we used 20 qPCR assays targeting 10 enteric pathogens, a fecal indicator bacterium (general Escherichia coli), and a process control strain in the microfluidic qPCR system. We preamplified the template DNA to increase the sensitivity of the qPCR assays. Our results suggested that preamplification was effective for quantifying small amounts of the template DNA without any major impact on the sensitivity, efficiency, and quantitative performance of qPCR. This microfluidic qPCR system allowed us to detect and quantify multiple pathogens from fecal samples and environmental water samples spiked with pathogens at levels as low as 100 cells/liter. These results suggest that the routine monitoring of multiple pathogens in food and water samples is now technically feasible. This method may provide more reliable information for risk assessment than the current fecal contamination indicator approach.

Natural genome diversity of AI-2 quorum sensing in Escherichia coli: conserved signal production but labile signal reception

Quorum sensing (QS) regulates the onset of bacterial social responses in function to cell density having an important impact in virulence. Autoinducer-2 (AI-2) is a signal that has the peculiarity of mediating both intra- and interspecies bacterial QS. We analyzed the diversity of all components of AI-2 QS across 44 complete genomes of Escherichia coli and Shigella strains. We used phylogenetic tools to study its evolution and determined the phenotypes of single-deletion mutants to predict phenotypes of natural strains. Our analysis revealed many likely adaptive polymorphisms both in gene content and in nucleotide sequence. We show that all natural strains possess the signal emitter (the luxS gene), but many lack a functional signal receptor (complete lsr operon) and the ability to regulate extracellular signal concentrations. This result is in striking contrast with the canonical species-specific QS systems where one often finds orphan receptors, without a cognate synthase, but not orphan emitters. Our analysis indicates that selection actively maintains a balanced polymorphism for the presence/absence of a functional lsr operon suggesting diversifying selection on the regulation of signal accumulation and recognition. These results can be explained either by niche-specific adaptation or by selection for a coercive behavior where signal-blind emitters benefit from forcing other individuals in the population to haste in cooperative behaviors.

Estimating variation within the genes and inferring the phylogeny of 186 sequenced diverse Escherichia coli genomes

Background

Escherichia coli exists in commensal and pathogenic forms. By measuring the variation of individual genes across more than a hundred sequenced genomes, gene variation can be studied in detail, including the number of mutations found for any given gene. This knowledge will be useful for creating better phylogenies, for determination of molecular clocks and for improved typing techniques.

Results

We find 3,051 gene clusters/families present in at least 95% of the genomes and 1,702 gene clusters present in 100% of the genomes. The former 'soft core' of about 3,000 gene families is perhaps more biologically relevant, especially considering that many of these genome sequences are draft quality. The E. coli pan-genome for this set of isolates contains 16,373 gene clusters.

A core-gene tree, based on alignment and a pan-genome tree based on gene presence/absence, maps the relatedness of the 186 sequenced E. coli genomes. The core-gene tree displays high confidence and divides the E. coli strains into the observed MLST type clades and also separates defined phylotypes.

Conclusion

The results of comparing a large and diverse E. coli dataset support the theory that reliable and good resolution phylogenies can be inferred from the core-genome. The results further suggest that the resolution at the isolate level may, subsequently be improved by targeting more variable genes. The use of whole genome sequencing will make it possible to eliminate, or at least reduce, the need for several typing steps used in traditional epidemiology.

A phylogenomic analysis of Escherichia coli / Shigella group: implications of genomic features associated with pathogenicity and ecological adaptation

Background

The Escherichia coli species contains a variety of commensal and pathogenic strains, and its intraspecific diversity is extraordinarily high. With the availability of an increasing number of E. coli strain genomes, a more comprehensive concept of their evolutionary history and ecological adaptation can be developed using phylogenomic analyses. In this study, we constructed two types of whole-genome phylogenies based on 34 E. coli strains using collinear genomic segments. The first phylogeny was based on the concatenated collinear regions shared by all of the studied genomes, and the second phylogeny was based on the variable collinear regions that are absent from at least one genome. Intuitively, the first phylogeny is likely to reveal the lineal evolutionary history among these strains (i.e., an evolutionary phylogeny), whereas the latter phylogeny is likely to reflect the whole-genome similarities of extant strains (i.e., a similarity phylogeny).

Results

Within the evolutionary phylogeny, the strains were clustered in accordance with known phylogenetic groups and phenotypes. When comparing evolutionary and similarity phylogenies, a concept emerges that Shigella may have originated from at least three distinct ancestors and evolved into a single clade. By scrutinizing the properties that are shared amongst Shigella strains but missing in other E. coli genomes, we found that the common regions of the Shigella genomes were mainly influenced by mobile genetic elements, implying that they may have experienced convergent evolution via horizontal gene transfer. Based on an inspection of certain key branches of interest, we identified several collinear regions that may be associated with the pathogenicity of specific strains. Moreover, by examining the annotated genes within these regions, further detailed evidence associated with pathogenicity was revealed.

Conclusions

Collinear regions are reliable genomic features used for phylogenomic analysis among closely related genomes while linking the genomic diversity with phenotypic differences in a meaningful way. The pathogenicity of a strain may be associated with both the arrival of virulence factors and the modification of genomes via mutations. Such phylogenomic studies that compare collinear regions of whole genomes will help to better understand the evolution and adaptation of closely related microbes and E. coli in particular.

Within-species lateral genetic transfer and the evolution of transcriptional regulation in Escherichia coli and Shigella

BACKGROUND

Changes in transcriptional regulation underlie many of the phenotypic differences observed within and between species of bacteria. Lateral genetic transfer (LGT) can significantly impact the transcription factor (TF) genes which drive these transcriptional changes. Although much emphasis has been placed on LGT of intact genes, the units of transfer and recombination do not necessarily correspond to regions delineated by exact gene boundaries. Here we apply phylogenetic and network-based methods to investigate the relationship between units of lateral transfer and recombination within the Escherichia coli - Shigella clade and the topological properties of genes in the E. coli transcriptional regulatory network (TRN).

RESULTS

We carried out a systematic phylogenetic study of genetic transfer among 5282 sets of putatively orthologous genes from 27 strains belonging to the E. coli - Shigella clade. We then used these results to examine the evolutionary histories of TF genes, as well as the transcriptional regulation of lateral genes. We found evidence of LGT in 2655 (50.3%) gene sets: 678 (12.8%) show evidence of recombination breakpoints within the gene boundaries. Thus, within- and whole- gene lateral transfer is widespread among strains of E. coli and Shigella. We found that unlike global regulators, which have mostly evolved vertically, neighbour regulators (genes which regulate adjacent genes on the chromosome) have frequently been subject to transfer within the E. coli - Shigella clade. At least 56 (62%) of the 90 neighbour regulator gene sets examined show evidence of LGT, 19 (34%) of which have internal recombination breakpoints. Neighbour regulators show no evidence of co-transfer with their nearby target genes. Rather, the frequency of recombination breakpoints, and conflicting evolutionary histories among neighbour regulators and their target genes, suggest that the genomic regions encoding these genes have been constructed through successive layering of LGT events within the clade. We find no difference in the relative complexity of regulation (i.e. the number of regulators) of lateral versus vertical genes.

CONCLUSIONS

Neighbour regulators show higher frequencies of transfer than other types of regulatory genes. This implicates the topological properties of regulatory genes in the TRN, and their physical proximity to targets on the chromosome, as contributing to successful LGT. The prevalence of recombination breakpoints within regulatory and target gene sets indicates that within-gene transfer has had a significant cumulative effect on the evolution of regulatory interactions in E. coli and Shigella.

Whole-genome phylogeny of Escherichia coli/Shigella group by feature frequency profiles (FFPs)

A whole-genome phylogeny of the Escherichia coli/Shigella group was constructed by using the feature frequency profile (FFP) method. This alignment-free approach uses the frequencies of l-mer features of whole genomes to infer phylogenic distances. We present two phylogenies that accentuate different aspects of E. coli/Shigella genomic evolution: (i) one based on the compositions of all possible features of length l = 24 (∼8.4 million features), which are likely to reveal the phenetic grouping and relationship among the organisms and (ii) the other based on the compositions of core features with low frequency and low variability (∼0.56 million features), which account for ∼69% of all commonly shared features among 38 taxa examined and are likely to have genome-wide lineal evolutionary signal. Shigella appears as a single clade when all possible features are used without filtering of noncore features. However, results using core features show that Shigella consists of at least two distantly related subclades, implying that the subclades evolved into a single clade because of a high degree of convergence influenced by mobile genetic elements and niche adaptation. In both FFP trees, the basal group of the E. coli/Shigella phylogeny is the B2 phylogroup, which contains primarily uropathogenic strains, suggesting that the E. coli/Shigella ancestor was likely a facultative or opportunistic pathogen. The extant commensal strains diverged relatively late and appear to be the result of reductive evolution of genomes. We also identify clade distinguishing features and their associated genomic regions within each phylogroup. Such features may provide useful information for understanding evolution of the groups and for quick diagnostic identification of each phylogroup.

The decoupling between genetic structure and metabolic phenotypes in Escherichia coli leads to continuous phenotypic diversity

To assess the extent of intra-species diversity and the links between phylogeny, lifestyle (habitat and pathogenicity) and phenotype, we assayed the growth yield on 95 carbon sources of 168 Escherichia strains. We also correlated the growth capacities of 14 E. coli strains with the presence/absence of enzyme-coding genes. Globally, we found that the genetic distance, based on multilocus sequence typing data, was a weak indicator of the metabolic phenotypic distance. Besides, lifestyle and phylogroup had almost no impact on the growth yield of non-Shigella E. coli strains. In these strains, the presence/absence of the metabolic pathways, which was linked to the phylogeny, explained most of the growth capacities. However, few discrepancies blurred the link between metabolic phenotypic distance and metabolic pathway distance. This study shows that a prokaryotic species structured into well-defined genetic and lifestyle groups can yet exhibit continuous phenotypic diversity, possibly caused by gene regulatory effects.

Representative proteomes: a stable, scalable and unbiased proteome set for sequence analysis and functional annotation

The accelerating growth in the number of protein sequences taxes both the computational and manual resources needed to analyze them. One approach to dealing with this problem is to minimize the number of proteins subjected to such analysis in a way that minimizes loss of information. To this end we have developed a set of Representative Proteomes (RPs), each selected from a Representative Proteome Group (RPG) containing similar proteomes calculated based on co-membership in UniRef50 clusters. A Representative Proteome is the proteome that can best represent all the proteomes in its group in terms of the majority of the sequence space and information. RPs at 75%, 55%, 35% and 15% co-membership threshold (CMT) are provided to allow users to decrease or increase the granularity of the sequence space based on their requirements. We find that a CMT of 55% (RP55) most closely follows standard taxonomic classifications. Further analysis of this set reveals that sequence space is reduced by more than 80% relative to UniProtKB, while retaining both sequence diversity (over 95% of InterPro domains) and annotation information (93% of experimentally characterized proteins). All sets can be browsed and are available for sequence similarity searches and download at http://www.proteininformationresource.org/rps, while the set of 637 RPs determined using a 55% CMT are also available for text searches. Potential applications include sequence similarity searches, protein classification and targeted protein annotation and characterization.

Animal and human pathogenic Escherichia coli strains share common genetic backgrounds

Escherichia coli is a versatile species encompassing both commensals of the digestive tracts of many vertebrates, including humans, and pathogenic strains causing various intra- and extraintestinal infections. Despite extensive gene flow between strains, the E. coli species has a globally clonal population structure, consisting of distinct phylogenetic groups. Little is known about the relationships between phylogenetic groups and host specificity. We therefore used multilocus sequence typing (MLST) to investigate phylogenetic relationships and evaluated the virulence gene content of 35 E. coli strains representative of the diverse diseases encountered in domestic animals. We compared these strains with a panel of 101 human pathogenic and 98 non-human and human commensal strains representative of the phylogenetic and pathovar diversity of this species. A global factorial analysis of correspondence indicated that extraintestinal infections were caused mostly by phylogenetic group B2 strains, whereas intraintestinal infections were caused mostly by phylogenetic group A/B1/E strains, with strains responsible from extraintestinal or intraintestinal infections having specific virulence factors. It was not possible to distinguish between strains of human and animal origin. A detailed phylogenetic analysis of the MLST data showed that numerous pathogenic animal and human strains are very closely related, and had a number of virulence genes in common. However, a set of specific adhesins was identified in animal non-B2 group strains of all pathotypes. In conclusion, human and animal pathogenic strains share common genetic backgrounds, but non-B2 strains of different origins seem to have different sets of adhesins that could be involved in host specificity.

Core and panmetabolism in Escherichia coli

Escherichia coli exhibits a wide range of lifestyles encompassing commensalism and various pathogenic behaviors which its highly dynamic genome contributes to develop. How environmental and host factors shape the genetic structure of E. coli strains remains, however, largely unknown. Following a previous study of E. coli genomic diversity, we investigated its diversity at the metabolic level by building and analyzing the genome-scale metabolic networks of 29 E. coli strains (8 commensal and 21 pathogenic strains, including 6 Shigella strains). Using a tailor-made reconstruction strategy, we significantly improved the completeness and accuracy of the metabolic networks over default automatic reconstruction processes. Among the 1,545 reactions forming E. coli panmetabolism, 885 reactions were common to all strains. This high proportion of core reactions (57%) was found to be in sharp contrast to the low proportion (13%) of core genes in the E. coli pangenome, suggesting less diversity of metabolic functions compared to that of all gene functions. Core reactions were significantly overrepresented among biosynthetic reactions compared to the more variable degradation processes. Differences between metabolic networks were found to follow E. coli phylogeny rather than pathogenic phenotypes, except for Shigella networks, which were significantly more distant from the others. This suggests that most metabolic changes in non-Shigella strains were not driven by their pathogenic phenotypes. Using a supervised method, we were yet able to identify small sets of reactions related to pathogenicity or commensalism. The quality of our reconstructed networks also makes them reliable bases for building metabolic models.

New enterovirulent Escherichia coli serogroup 64474 showing antigenic and genotypic relationships to Shigella boydii 16

Studies based on the analysis of housekeeping genes indicate that Escherichia coli and all Shigella species, except for Shigella boydii type 13, belong to a single species. This study analysed the phenotypic and genotypic characteristics of 23 E. coli strains isolated in different countries from faecal specimens taken from children with diarrhoea. Strains were identified using the VITEK system and typed with rabbit sera obtained against 186 somatic and 53 flagellar E. coli antigens and against 45 Shigella somatic antigens. Biochemical analysis of these strains showed a typical E. coli profile with a defined reaction against both E. coli O179 and S. boydii 16 somatic antisera. Agglutination assays for flagellar antigens showed a response against H2 in 7 (30 %) strains, H10 in 2 (9 %) strains, H32 in 12 (52 %) strains and H34 in 2 (9 %) strains, demonstrating 4 serotypes associated with this new somatic antigen 64474. A serum against one of these E. coli strains (64474) was prepared. Absorption assays of S. boydii 16 and E. coli 64474 antisera with E. coli O179 antigen removed the agglutination response against this O179 antigen completely, while the agglutination titres against both S. boydii 16 and E. coli 64474 remained the same. Four (17 %) E. coli strains showed antimicrobial resistance to piperacillin only, one (4 %) to piperacillin and trimethoprim/sulfamethoxazole, one (4 %) to ciprofloxacin, nitrofurantoin and piperacillin, and two (9 %) strains were resistant to ciprofloxacin, norfloxacin, ofloxacin, piperacillin and trimethoprim/sulfamethoxazole. With regards to PCR assays, one (4 %) of the strains was positive for Shigella gene ipaH, one (4 %) for ipaA, two (9 %) for ipaB, one (4 %) for ipaD, two (9 %) for sepA and three (13 %) for ospF. The rfb gene cluster in the E. coli strains was analysed by RFLP and compared with the gene cluster obtained from S. boydii 16. The rfb-RFLP patterns for all 23 E. coli strains were similar to those obtained for S. boydii 16. The results from PCR tests to detect rfb genes wzx (encoding O unit flippase) and wzy (encoding polymerase) belonging to a cluster related to the biosynthesis of the S. boydii 16-specific O antigen were positive in 21 (91 %) and 22 (96 %) of the strains, respectively. PCR assays to detect E. coli virulence genes were also performed. These assays detected enterotoxigenic E. coli genes ltA1 in 12 of the strains (52 %), st1a in 4 (17 %), cfa1 in 6 (26 %), cs1 in 1 (4 %), cs3 in 3 (13 %), cs13 in 9 (39 %) and cs14 in 5 (22 %) of the strains. Results from the PFGE analyses confirmed the wide geographical distribution of these strains suggesting that 64474 : H2, 64474 : H10, 64474 : H32 and 64474 : H34 are new serotypes of E. coli strains with a defined virulence capacity, and share a common O antigen with S. boydii 16.

Whole-proteome phylogeny of prokaryotes by feature frequency profiles: An alignment-free method with optimal feature resolution

We present a whole-proteome phylogeny of prokaryotes constructed by comparing feature frequency profiles (FFPs) of whole proteomes. Features are l-mers of amino acids, and each organism is represented by a profile of frequencies of all features. The selection of feature length is critical in the FFP method, and we have developed a procedure for identifying the optimal feature lengths for inferring the phylogeny of prokaryotes, strictly speaking, a proteome phylogeny. Our FFP trees are constructed with whole proteomes of 884 prokaryotes, 16 unicellular eukaryotes, and 2 random sequences. To highlight the branching order of major groups, we present a simplified proteome FFP tree of monophyletic class or phylum with branch support. In our whole-proteome FFP trees (i) Archaea, Bacteria, Eukaryota, and a random sequence outgroup are clearly separated; (ii) Archaea and Bacteria form a sister group when rooted with random sequences; (iii) Planctomycetes, which possesses an intracellular membrane compartment, is placed at the basal position of the Bacteria domain; (iv) almost all groups are monophyletic in prokaryotes at most taxonomic levels, but many differences in the branching order of major groups are observed between our proteome FFP tree and trees built with other methods; and (v) previously "unclassified" genomes may be assigned to the most likely taxa. We describe notable similarities and differences between our FFP trees and those based on other methods in grouping and phylogeny of prokaryotes.

Gene and genome trees conflict at many levels

Horizontal gene transfer (HGT) plays a significant role in microbial evolution. It can accelerate the adaptation of an organism, it can generate new metabolic pathways and it can completely remodel an organism's genome. We examine 27 closely related genomes from the YESS group of gamma proteobacteria and a variety of four-taxon datasets from a diverse range of prokaryotes in order to explore the kinds of effects HGT has had on these organisms.

Cryptic lineages of the genus Escherichia

Extended multilocus sequence typing (MLST) analysis of atypical Escherichia isolates was used to identify five novel phylogenetic clades (CI to CV) among isolates from environmental, human, and animal sources. Analysis of individual housekeeping loci showed that E. coli and its sister clade, CI, remain largely indistinguishable and represent nascent evolutionary lineages. Conversely, clades of similar age (CIII and CIV) were found to be phylogenetically distinct. When all Escherichia lineages (named and unnamed) were evaluated, we found evidence that Escherichia fergusonii has evolved at an accelerated rate compared to E. coli, CI, CIII, CIV, and CV, suggesting that this species is younger than estimated by the molecular clock method. Although the five novel clades were phylogenetically distinct, we were unable to identify a discriminating biochemical marker for all but one of them (CIII) with traditional phenotypic profiling. CIII had a statistically different phenotype from E. coli that resulted from the loss of sucrose and sorbitol fermentation and lysine utilization. The lack of phenotypic distinction has likely hindered the ability to differentiate these clades from typical E. coli, and so their ecological significance and importance for applied and clinical microbiology are yet to be determined. However, our sampling suggests that CIII, CIV, and CV represent environmentally adapted Escherichia lineages that may be more abundant outside the host gastrointestinal tract.

Omptin proteins: an expanding family of outer membrane proteases in Gram-negativeEnterobacteriaceae(Review)

The Escherichia coli K-12 outer membrane protein OmpT is a prototype of a unique family of bacterial endopeptidases known as the omptins. This family includes OmpT and OmpP of E. coli, SopA of Shigella flexneri, PgtE of Salmonella enterica, and Pla of Yersinia pestis. Despite their sequence similarities, the omptins vary in their reported functions. The OmpT protease is characterized by narrow cleavage specificity defined by the extracellular loops of the beta-barrel protruding above the lipid bilayer. It employs a distinct proteolytic mechanism that involves a histidine and an aspartate residue. Most of the omptin proteins have been implicated in bacterial pathogenesis. As a result, the omptins are potential targets for antimicrobial drug and vaccine development. This review summarizes recent developments in omptins structure and function, emphasizes their role in pathogenesis, proposes evolutionary relation among the existing omptins, and offers possible directions for future research.

Fate of the H-NS-repressed bgl operon in evolution of Escherichia coli

In the enterobacterial species Escherichia coli and Salmonella enterica, expression of horizontally acquired genes with a higher than average AT content is repressed by the nucleoid-associated protein H-NS. A classical example of an H-NS–repressed locus is the bgl (aryl-β,D-glucoside) operon of E. coli. This locus is “cryptic,” as no laboratory growth conditions are known to relieve repression of bgl by H-NS in E. coli K12. However, repression can be relieved by spontaneous mutations. Here, we investigated the phylogeny of the bgl operon. Typing of bgl in a representative collection of E. coli demonstrated that it evolved clonally and that it is present in strains of the phylogenetic groups A, B1, and B2, while it is presumably replaced by a cluster of ORFans in the phylogenetic group D. Interestingly, the bgl operon is mutated in 20% of the strains of phylogenetic groups A and B1, suggesting erosion of bgl in these groups. However, bgl is functional in almost all B2 isolates and, in approximately 50% of them, it is weakly expressed at laboratory growth conditions. Homologs of bgl genes exist in Klebsiella, Enterobacter, and Erwinia species and also in low GC-content Gram-positive bacteria, while absent in E. albertii and Salmonella sp. This suggests horizontal transfer of bgl genes to an ancestral Enterobacterium. Conservation and weak expression of bgl in isolates of phylogenetic group B2 may indicate a functional role of bgl in extraintestinal pathogenic E. coli.

Horizontal gene transfer, an important mechanism in bacterial adaptation and evolution, requires mechanisms to avoid uncontrolled and possibly disadvantageous expression of the transferred genes. Recently, it was shown that the protein H-NS selectively silences genes gained by horizontal transfer in enteric bacteria. Regulated expression of these genes can then evolve and be integrated into the regulatory network of the new host. Our analysis of the catabolic bgl (aryl-β,D-glucoside) operon, which is silenced by H-NS in E. coli, provides a snapshot on the evolution of such a locus. Genes of the bgl operon were presumably gained by horizontal transfer from Gram-positive bacteria to ancestral enteric bacteria. In E. coli, the bgl operon co-evolved with the diversification of the species into four phylogenetic groups. In one phylogenetic group the bgl operon is functional. However, in two other phylogenetic groups, bgl accumulates disrupting mutations, and it is absent in the fourth group. This indicates that the H-NS–silenced bgl operon evolved differently in E. coli and is presumably positively selected in one phylogenetic group, while it is neutrally or negatively selected in the other groups.

Organised genome dynamics in the Escherichia coli species results in highly diverse adaptive paths

The Escherichia coli species represents one of the best-studied model organisms, but also encompasses a variety of commensal and pathogenic strains that diversify by high rates of genetic change. We uniformly (re-) annotated the genomes of 20 commensal and pathogenic E. coli strains and one strain of E. fergusonii (the closest E. coli related species), including seven that we sequenced to completion. Within the ∼18,000 families of orthologous genes, we found ∼2,000 common to all strains. Although recombination rates are much higher than mutation rates, we show, both theoretically and using phylogenetic inference, that this does not obscure the phylogenetic signal, which places the B2 phylogenetic group and one group D strain at the basal position. Based on this phylogeny, we inferred past evolutionary events of gain and loss of genes, identifying functional classes under opposite selection pressures. We found an important adaptive role for metabolism diversification within group B2 and Shigella strains, but identified few or no extraintestinal virulence-specific genes, which could render difficult the development of a vaccine against extraintestinal infections. Genome flux in E. coli is confined to a small number of conserved positions in the chromosome, which most often are not associated with integrases or tRNA genes. Core genes flanking some of these regions show higher rates of recombination, suggesting that a gene, once acquired by a strain, spreads within the species by homologous recombination at the flanking genes. Finally, the genome's long-scale structure of recombination indicates lower recombination rates, but not higher mutation rates, at the terminus of replication. The ensuing effect of background selection and biased gene conversion may thus explain why this region is A+T-rich and shows high sequence divergence but low sequence polymorphism. Overall, despite a very high gene flow, genes co-exist in an organised genome.

Although abundant knowledge has been accumulated regarding the E. coli laboratory strain K-12, little is known about the evolutionary trajectories that have driven the high diversity observed among natural isolates of the species, which encompass both commensal and highly virulent intestinal and extraintestinal pathogenic strains. We have annotated or re-annotated the genomes of 20 commensal and pathogenic E. coli strains and one strain of E. fergusonii (the closest E. coli related species), including seven that we sequenced to completion. Although recombination rates are much higher than mutation rates, we were able to reconstruct a robust phylogeny based on the ∼2,000 genes common to all strains. Based on this phylogeny, we established the evolutionary scenario of gains and losses of thousands of specific genes, identifying functional classes under opposite selection pressures. This genome flux is confined to very few positions in the chromosome, which are the same for every genome. Notably, we identified few or no extraintestinal virulence-specific genes. We also defined a long-scale structure of recombination in the genome with lower recombination rates at the terminus of replication. These findings demonstrate that, despite a very high gene flow, genes can co-exist in an organised genome.

Analysis of plasmid genes by phylogenetic profiling and visualization of homology relationships using Blast2Network

BACKGROUND

Phylogenetic methods are well-established bioinformatic tools for sequence analysis, allowing to describe the non-independencies of sequences because of their common ancestor. However, the evolutionary profiles of bacterial genes are often complicated by hidden paralogy and extensive and/or (multiple) horizontal gene transfer (HGT) events which make bifurcating trees often inappropriate. In this context, plasmid sequences are paradigms of network-like relationships characterizing the evolution of prokaryotes. Actually, they can be transferred among different organisms allowing the dissemination of novel functions, thus playing a pivotal role in prokaryotic evolution. However, the study of their evolutionary dynamics is complicated by the absence of universally shared genes, a prerequisite for phylogenetic analyses.

RESULTS

To overcome such limitations we developed a bioinformatic package, named Blast2Network (B2N), allowing the automatic phylogenetic profiling and the visualization of homology relationships in a large number of plasmid sequences. The software was applied to the study of 47 completely sequenced plasmids coming from Escherichia, Salmonella and Shigella spps.

CONCLUSION

The tools implemented by B2N allow to describe and visualize in a new way some of the evolutionary features of plasmid molecules of Enterobacteriaceae; in particular it helped to shed some light on the complex history of Escherichia, Salmonella and Shigella plasmids and to focus on possible roles of unannotated proteins.The proposed methodology is general enough to be used for comparative genomic analyses of bacteria.

Safety aspects and implications of regulation of probiotic bacteria in food and food supplements

The application of living bacteria as probiotics in food or food supplements requires a careful safety assessment. This review summarizes key issues concerning the safety aspects of bacteria added to particular products marketed for improvement of general health or treatment of (post)infectious symptoms. The bacteria used in such products should be completely safe; however, it can be challenging to provide evidence for absence of all virulence properties. In some cases, virulence factors have been detected in probiotic bacterial strains, and the implications of these traits for safety assessments are discussed. Horizontal gene transfer can result in acquisition of virulence genes or antimicrobial resistance in probiotic bacteria. Antimicrobial resistance in these bacteria can possibly aid the spread of undesired resistance in intestinal bacterial populations. The relative risk of such gene transfers is considered. The generation of complete bacterial genome sequences can both resolve and create safety issues. Current practices of safety assessment procedures in the United States and the European Union are briefly reviewed and a future outlook is provided.

The pangenome structure of Escherichia coli: comparative genomic analysis of E. coli commensal and pathogenic isolates

Whole-genome sequencing has been skewed toward bacterial pathogens as a consequence of the prioritization of medical and veterinary diseases. However, it is becoming clear that in order to accurately measure genetic variation within and between pathogenic groups, multiple isolates, as well as commensal species, must be sequenced. This study examined the pangenomic content of Escherichia coli. Six distinct E. coli pathovars can be distinguished using molecular or phenotypic markers, but only two of the six pathovars have been subjected to any genome sequencing previously. Thus, this report provides a seminal description of the genomic contents and unique features of three unsequenced pathovars, enterotoxigenic E. coli, enteropathogenic E. coli, and enteroaggregative E. coli. We also determined the first genome sequence of a human commensal E. coli isolate, E. coli HS, which will undoubtedly provide a new baseline from which workers can examine the evolution of pathogenic E. coli. Comparison of 17 E. coli genomes, 8 of which are new, resulted in identification of approximately 2,200 genes conserved in all isolates. We were also able to identify genes that were isolate and pathovar specific. Fewer pathovar-specific genes were identified than anticipated, suggesting that each isolate may have independently developed virulence capabilities. Pangenome calculations indicate that E. coli genomic diversity represents an open pangenome model containing a reservoir of more than 13,000 genes, many of which may be uncharacterized but important virulence factors. This comparative study of the species E. coli, while descriptive, should provide the basis for future functional work on this important group of pathogens.

Prevalence and sequence diversity of a factor required for actin-based motility in natural populations of Burkholderia species

Actin-based motility of the melioidosis pathogen Burkholderia pseudomallei requires BimA. We report a high degree of conservation of bimA in 99 B. pseudomallei isolates from the area of endemicity. A geographically restricted subset of B. pseudomallei isolates harbored a B. mallei-like bimA allele (12.1%), confounding a differential diagnostic test based on amplification of species-specific bimA regions.

Characteristics of human intestinal Escherichia coli with changing environments

To investigate if the characteristics of human intestinal Escherichia coli are changing with the environment of the host, we studied intestinal E. coli from subjects having recently migrated from a temperate to a tropical area. We determined the phylogenetic group, the prevalence of the antibiotic resistance, the presence of integrons and the strain diversity in faecal isolates from 25 subjects originally from metropolitan France and expatriated to French Guyana. These characteristics were compared with those of 25 previously studied Wayampi Amerindian natives of French Guyana and from 25 metropolitan French residents. The three groups of subjects were matched for age and sex, had not taken antibiotics for at least 1 month, nor had been hospitalized within the past year. In all, the characteristics of intestinal E. coli from Expatriates were intermediate between those found in residents from metropolitan France and those found in natives of French Guyana. Prevalence of carriage of resistant Gram-negative bacteria in Expatriates was intermediate between French residents and Wayampi as were the prevalence of integrons in E. coli (12.3% versus 16.3% and 7.8% respectively), and the intra-host diversity of E. coli (2.3 strains/subject versus 1.9 and 3.1, respectively); lastly, in Expatriates, the prevalence of carriage of phylogenetic group B2 strains was lower than in French residents (16% versus 56%, P = 0.005), while carriage of phylogenetic group A strains was lower than in Wayampi (56% versus 88%, P = 0.03). Our results suggest that the composition of the commensal intestinal flora of humans is not static but changes dynamically in response to new environmental conditions.

Molecular pathogenesis of Shigella spp.: controlling host cell signaling, invasion, and death by type III secretion

Shigella spp. are gram-negative pathogenic bacteria that evolved from harmless enterobacterial relatives and may cause devastating diarrhea upon ingestion. Research performed over the last 25 years revealed that a type III secretion system (T3SS) encoded on a large plasmid is a key virulence factor of Shigella flexneri. The T3SS determines the interactions of S. flexneri with intestinal cells by consecutively translocating two sets of effector proteins into the target cells. Thus, S. flexneri controls invasion into EC, intra- and intercellular spread, macrophage cell death, as well as host inflammatory responses. Some of the translocated effector proteins show novel biochemical activities by which they intercept host cell signal transduction pathways. An understanding of the molecular mechanisms underlying Shigella pathogenesis will foster the development of a safe and efficient vaccine, which, in parallel with improved hygiene, should curb infections by this widespread pathogen.

Selection and validation of a multilocus variable-number tandem-repeat analysis panel for typing Shigella spp

The Shigella genus has historically been separated into four species, based on biochemical assays. The classification within each species relies on serotyping. Recently, genome sequencing and DNA assays, in particular the multilocus sequence typing (MLST) approach, greatly improved the current knowledge of the origin and phylogenetic evolution of Shigella spp. The Shigella and Escherichia genera are now considered to belong to a unique genomospecies. Multilocus variable-number tandem-repeat (VNTR) analysis (MLVA) provides valuable polymorphic markers for genotyping and performing phylogenetic analyses of highly homogeneous bacterial pathogens. Here, we assess the capability of MLVA for Shigella typing. Thirty-two potentially polymorphic VNTRs were selected by analyzing in silico five Shigella genomic sequences and subsequently evaluated. Eventually, a panel of 15 VNTRs was selected (i.e., MLVA15 analysis). MLVA15 analysis of 78 strains or genome sequences of Shigella spp. and 11 strains or genome sequences of Escherichia coli distinguished 83 genotypes. Shigella population cluster analysis gave consistent results compared to MLST. MLVA15 analysis showed capabilities for E. coli typing, providing classification among pathogenic and nonpathogenic E. coli strains included in the study. The resulting data can be queried on our genotyping webpage (http://mlva.u-psud.fr). The MLVA15 assay is rapid, highly discriminatory, and reproducible for Shigella and Escherichia strains, suggesting that it could significantly contribute to epidemiological trace-back analysis of Shigella infections and pathogenic Escherichia outbreaks. Typing was performed on strains obtained mostly from collections. Further studies should include strains of much more diverse origins, including all pathogenic E. coli types.

Multilocus sequence typing analysis of Shigella flexneri isolates collected in Asian countries

The multilocus sequence typing scheme used previously for phylogenetic analysis of Escherichia coli was applied to 107 clinical isolates of Shigella flexneri. DNA sequencing of 3423 bp throughout seven housekeeping genes identified eight new allele types and ten new sequence types among the isolates. S. flexneri serotypes 1-5, X and Y were clustered together in a group containing many allelic variants while serotype 6 formed a distinct group, as previously established.

Molecular evolution of pathogenicity-island genes in Pseudomonas viridiflava

The bacterial pathogen Pseudomonas viridiflava possesses two pathogenicity islands (PAIs) that share many gene homologs, but are structurally and phenotypically differentiated (T-PAI and S-PAI). These PAIs are paralogous, but only one is present in each isolate. While this dual presence/absence polymorphism has been shown to be maintained by balancing selection, little is known about the molecular evolution of individual genes on the PAIs. Here we investigate genetic variation of 12 PAI gene loci (7 on T-PAI and 5 on S-PAI) in 96 worldwide isolates of P. viridiflava. These genes include avirulence genes (hopPsyA and avrE), their putative chaperones (shcA and avrF), and genes encoding the type III outer proteins (hrpA, hrpZ, and hrpW). Average nucleotide diversities in these genes (pi = 0.004-0.020) were close to those in the genetic background. Large numbers of recombination events were found within PAIs and a sign of positive selection was detected in avrE. These results suggest that the PAI genes are evolving relatively freely from each other on the PAIs, rather than as a single unit under balancing selection. Evolutionarily stable PAIs may be preferable in this species because preexisting genetic variation enables P. viridiflava to respond rapidly to natural selection.