Genome sequence of Shigella flexneri 2a: insights into pathogenicity through comparison with genomes of Escherichia coli K12 and O157

The two-faced role of cad genes in the virulence of pathogenic Escherichia coli

In enterobacteria, acid stress induces expression of the cad system which is involved in maintaining intracellular pH at levels compatible with cell survival. Despite its crucial role, the cad operon is silenced in Shigella and in other pathogenic Escherichia coli. In the present review, we will address the question of why and how the cad locus has been sacrificed for the sake of optimal expression of virulence traits.

Co-evolution of transcription factors and their targets depends on mode of regulation

Analysis of transcription regulatory networks in γ-proteobacteria reveals that repressors co-evolve tightly with their target genes, whereas activators can be lost independently of their targets.

Background

Differences in the transcription regulation network are at the root of much of the phenotypic variation observed among organisms. These differences may be achieved either by changing the repertoire of regulators and/or their targets, or by rewiring the network. Following these changes and studying their logic is crucial for understanding the evolution of regulatory networks.

Results

We use the well characterized transcription regulatory network of Escherichia coli K12 and follow the evolutionary changes in the repertoire of regulators and their targets across a large number of fully sequenced γ-proteobacteria. By focusing on close relatives of E. coli K12, we study the dynamics of the evolution of transcription regulation across a relatively short evolutionary timescale. We show significant differences in the evolution of repressors and activators. Repressors are only lost from a genome once their targets have themselves been lost, or once the network has significantly rewired. In contrast, activators are often lost even when their targets remain in the genome. As a result, E. coli K12 repressors that regulate many targets are rarely absent from organisms that are closely related to E. coli K12, while activators with a similar number of targets are often absent in these organisms.

Conclusion

We demonstrate that the mode of regulation exerted by transcription factors has a strong effect on their evolution. Repressors co-evolve tightly with their target genes. In contrast, activators can be lost independently of their targets. In fact, loss of an activator can lead to efficient shutdown of an unnecessary pathway.

Polyketide synthases and nonribosomal peptide synthetases: the emerging view from bacterial genomics

A total of 223 complete bacterial genomes are analyzed, with 281 citations, for the presence of genes encoding modular polyketide synthases (PKS) and nonribosomal peptide synthetases (NRPS). We report on the distribution of these systems in different bacterial taxa and, whenever known, the metabolites they synthesize. We also highlight, in the different bacterial lineages, the PKS and NRPS genes and, whenever known, the corresponding products.

Comparative analysis of bacterial genomes: identification of divergent regions in mycobacterial strains using an anchor-based approach

Comparative genomic approaches are useful in identifying molecular differences between organisms. Currently available methods fail to identify small changes in genomes, such as expansion of short repetitive motifs and to analyse divergent sequences. In this report, we describe an anchor-based whole genome comparison (ABWGC) method. ABWGC is based on random sampling of anchor sequences from one genome, followed by analysis of sampled and homologous regions from the target genome. The method was applied to compare two strains of Mycobacterium tuberculosis CDC1551 and H37Rv. ABWGC was able to identify a total of 104 indels including 20 expansion of short repetitive sequences and five recombination events. It included 18 new unidentified genomic differences. ABWGC also identified 188 SNPs including eight new ones. The method was also used to compare M. tuberculosis H37Rv and M. avium genomes. ABWGC was able to correctly pick 1002 additional indels (size >100 nt) between the two organisms in contrast to MUMmer, a popular tool for comparative genomics. ABWGC was able to identify correctly repeat expansion and indels in a set of simulated sequences. The study also revealed important role of small repeat expansion in the evolution of M. tuberculosis strains.

Erosion of interaction networks in reduced and degraded genomes

Unlike eukaryotes, which often recruit duplicated genes into existing protein-protein interaction (PPI) networks, the low levels of gene duplication coupled with the high probability of lateral transfer of novel genes alters the manner in which PPI networks can evolve in bacteria. By inferring the PPIs present in the ancestor to contemporary Gammaproteobacteria, we were able to trace the changes in gene repertoires, and their consequences on PPI network evolution, in several bacterial lineages that have independently undergone reductions in genome size and genome contents. As genomes degrade, virtually all multi-partner proteins have lost interactors; however, the overall average number of connections increases due to the preferential elimination of proteins that interact with only one other protein partner. We also studied the effect of lateral gene transfer on PPI network evolution by analyzing the connectivity of genes that have been gained along the Escherichia coli lineage, as well as those acquired genes subsequently silenced in Shigella flexneri, since diverging from the gammaproteobacterial ancestor. The situation in PPI networks, in which newly acquired genes preferentially attach to the hubs of the network, contrasts that observed in metabolic networks, which evolve by the peripheral gain and loss of genes, and in regulatory networks, in which high connectivity increases the propensity of loss.

Black holes, antivirulence genes, and gene inactivation in the evolution of bacterial pathogens

The evolution of bacterial pathogens from nonpathogenic ancestors is marked principally by the acquisition of virulence gene clusters on plasmids and pathogenicity islands via horizontal gene transfer. The flip side of this evolutionary force is the equally important adaptation of the newly minted pathogen to its new host niche. Pathoadaptive mutations take the form of modification of gene expression such that the pathogen is better fit to survive within the new niche. This mini-review describes the concept of pathoadaptation by loss of gene function. In this process, genes that are no longer compatible with the novel lifestyle of the pathogen are selectively inactivated either by point mutation, insertion, or deletion. These genes are called 'antivirulence genes'. Selective pressure sometimes leads to the deletion of large regions of the genome that contain antivirulence genes generating 'black holes' in the pathogen genome. Inactivation of antivirulence genes leads to a pathogen that is highly adapted to its host niche. Identification of antivirulence genes for a particular pathogen can lead to a better understanding of how it became a pathogen and the types of genetic traits that need to be silenced in order for the pathogen to colonize its new host niche successfully.

Detecting low concentrations of Shigella sonnei in environmental water samples by PCR

Outbreaks of Shigella sonnei associated with contaminated water have been reported and methods for the simultaneous detection of Shigellae and enteroinvasive Escherichia coli in water samples have been developed with detection limits of 10(1)-10(2) CFU mL(-1) of water. Because 10(1)-10(2)Shigellae can cause disease, a more sensitive detection method as an addition to the existing methods for detection of Shigella sonnei in water samples is reported here. Initially, 33 Shigella sonnei and 72 non-Shigella sonnei isolates were tested and one primer pair was found capable of specifically amplifying a 369-bp insertion sequence 1 (IS1) fragment from all 33 Shigella sonnei isolates and one Shigella dysenteriae ATCC isolate by PCR. The detection method was developed, which included filtration of 50 mL of water through a membrane and application of PCR to the membrane using this primer pair. Environmental water samples with total bacterial numbers of 384-2.84 x 10(7) CFU L(-1) were collected and seeded with 13 Shigella sonnei and the Shigella dysenteriae ATCC isolates. Detection limits were determined as 1.7-24.7 and 270-8000 CFU per 50 mL of water, respectively, using this detection method.

L-Sorbose utilization by virulent Escherichia coli and Shigella: different metabolic adaptation of pathotypes

The frequency of L-Sorbose utilization differs significantly between pathotypes of Escherichia coli and Shigella from 93% to 0%. Among 266 strains tested, this frequency increased in the order Shigella, enterotoxigenic E. coli (ETEC), enteroinvasive E. coli (EIEC), Shiga toxin-producing E. coli (STEC), enteroaggregative E. coli, enteropathogenic E. coli (EPEC), and neonatal bacterial meningitis (NBM) E. coli. This suggests an association of pathomechanism with the capability to degrade L-Sorbose. The use of a selective agar, containing L-Sorbose and antibiotics, facilitated the isolation of L-Sorbose-non-utilizing ETEC from stool specimens of patients. The sor operon, comprising seven genes in the order sorCDFBAME, confers L-Sorbose utilization. Surprisingly, L-Sorbose-non-degrading Shigella harbored all genes of the sor operon indicating L-Sorbose-utilizing E. coli as ancestor. Additionally, strains of several EIEC and STEC serotypes harbored an inactivated sor operon. These L-Sorbose-non-utilizing Shigella, EIEC, and STEC showed significantly reduced amounts of transcripts as examined for sorC and sorD. Common surface antigens, types of intimin gene, and hemolysin gene as well as use of L-Sorbose suggested the relatedness of attaching and effacing O26:H11 and O55:H7 EPEC and STEC, respectively. pepE and yibC genes flank the sor operon of E. coli and Shigella strains. Surprisingly, one O7:K1:H- NBM E. coli harbored an aroE-homologous gene between its sor operon and pepE as in Klebsiella pneumoniae suggesting a horizontal gene transfer. In conclusion, L-Sorbose utilization of virulent E. coli and Shigella is characterized by different adaptation that represents a valuable tool for evolutionary and diagnostic analysis of related patho- and serotypes.

Comparative genomics and the evolution of prokaryotes

Although biologists have long recognized the importance of studying evolution to understand the organization of living organisms, only with the development of genomics have evolutionary studies become part of their routine toolkit. Placing genomes into an evolutionary framework has proved useful for understanding the functioning of organisms. It has also substantially increased understanding of the processes by which genomes evolve and led to a re-evaluation of our representation of the diversity and the history of life. In this review, we present some of the most important recent advances and promising leads in the field of microbial evolutionary genomics.

Gene expression profiling of the pH response in Shigella flexneri 2a

The pH response of Shigella flexneri 2a 301 was identified by gene expression profiling. Gene expression profiles of cells grown in pH 4.5 or 8.6 were compared with the profiles of cells grown at pH 7.0. Differential expression was observed for 307 genes: 97 were acid up-regulated, 102 were acid down-regulated, 91 were base up-regulated, and 86 were base down-regulated. Twenty-seven genes were found to be both acid and base up-regulated, and 29 genes were both acid and base down-regulated. This study showed that (1) the most pH-dependent genes regulate energy metabolism; (2) the RpoS-dependent acid-resistance system is induced, while the glutamate-dependent acid resistance system is not; (3) high pH up-regulates some virulence genes, while low pH down-regulates them, consistent with Shigella infection of the low gut; and (4) several cross-stress response genes are induced by pH changes. These results also illustrate that many unknown genes are significantly regulated under acid or basic conditions, providing researchers with important information to characterize their function.

Causes of insertion sequences abundance in prokaryotic genomes

Insertion sequences (ISs) are the smallest and most frequent transposable elements in prokaryotes where they play an important evolutionary role by promoting gene inactivation and genome plasticity. Their genomic abundance varies by several orders of magnitude for reasons largely unknown and widely speculated. The current availability of hundreds of genomes renders testable many of these hypotheses, notably that IS abundance correlates positively with the frequency of horizontal gene transfer (HGT), genome size, pathogenicity, nonobligatory ecological associations, and human association. We thus reannotated ISs in 262 prokaryotic genomes and tested these hypotheses showing that when using appropriate controls, there is no empirical basis for IS family specificity, pathogenicity, or human association to influence IS abundance or density. HGT seems necessary for the presence of ISs, but cannot alone explain the absence of ISs in more than 20% of the organisms, some of which showing high rates of HGT. Gene transfer is also not a significant determinant of the abundance of IS elements in genomes, suggesting that IS abundance is controlled at the level of transposition and ensuing natural selection and not at the level of infection. Prokaryotes engaging in obligatory associations have fewer ISs when controlled for genome size, but this may be caused by some being sexually isolated. Surprisingly, genome size is the only significant predictor of IS numbers and density. Alone, it explains over 40% of the variance of IS abundance. Because we find that genome size and IS abundance correlate negatively with minimal doubling times, we conclude that selection for rapid replication cannot account for the few ISs found in small genomes. Instead, we show evidence that IS numbers are controlled by the frequency of highly deleterious insertion targets. Indeed, IS abundance increases quickly with genome size, which is the exact inverse trend found for the density of genes under strong selection such as essential genes. Hence, for ISs, the bigger the genome the better.

Autotransporter-encoding sequences are phylogenetically distributed among Escherichia coli clinical isolates and reference strains

Autotransporters are secreted bacterial proteins exhibiting diverse virulence functions. Various autotransporters have been identified among Escherichia coli associated with intestinal or extraintestinal infections; however, the specific distribution of autotransporter sequences among a diversity of E. coli strains has not been investigated. We have validated the use of a multiplex PCR assay to screen for the presence of autotransporter sequences. Herein, we determined the presence of 13 autotransporter sequences and five allelic variants of antigen 43 (Ag43) among 491 E. coli isolates from human urinary tract infections, diarrheagenic E. coli, and avian pathogenic E. coli (APEC) and E. coli reference strains belonging to the ECOR collection. Clinical isolates were also classified into established phylogenetic groups. The results indicated that Ag43 alleles were significantly associated with clinical isolates (93%) compared to commensal isolates (56%) and that agn43K12 was the most common and widely distributed allele. agn43 allelic variants were also phylogenetically distributed. Sequences encoding espC, espP, and sepA and agn43 alleles EDL933 and RS218 were significantly associated with diarrheagenic E. coli strains compared to other groups. tsh was highly associated with APEC strains, whereas sat was absent from APEC. vat, sat, and pic were associated with urinary tract isolates and were identified predominantly in isolates belonging to either group B2 or D of the phylogenetic groups based on the ECOR strain collection. Overall, the results indicate that specific autotransporter sequences are associated with the source and/or phylogenetic background of strains and suggest that, in some cases, autotransporter gene profiles may be useful for comparative analysis of E. coli strains from clinical, food, and environmental sources.

A global response to sulfur starvation in Pseudomonas putida and its relationship to the expression of low-sulfur-content proteins

Sulfur is essential for life on Earth, but its availability is limited in many environments. Here the sulfur-starvation response of the model soil bacterium Pseudomonas putida KT2440 is shown to be associated with an approximately fivefold reduction in the total soluble thiol content of the cell. A bioinformatic survey of the P. putida KT2440 genome identified 646 genes encoding proteins with a significantly lower than average sulfur content (low sulfur-content proteins, LSPs), the expression of which may have a role in the global reduction of cellular thiol content during sulfur starvation. Analysis of the genetic organization of the LSP-encoding genes showed that 31% were potentially transcriptionally associated with at least one other gene encoding a protein defined as an LSP. In particular, 55 LSP genes were located in three large clusters, termed low-sulfur islands (LSIs) here. The predicted identities of the proteins encoded by the LSIs strongly suggest that the LSIs have a role in acquiring sulfur from organic sulfur sources during sulfur starvation. This hypothesis was supported by transcription fusion studies on a limited number of LSP promoters under low-sulfur conditions. In a wider survey of bacterial species, LSIs were found to be more prevalent in free-living, Gram-negative bacteria than in Gram-positive or obligately intracellular bacteria.

Revisiting the Molecular Evolutionary History of Shigella 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.

Transcriptional slippage controls production of type III secretion apparatus components in Shigella flexneri

During transcription, series of approximately 9 As or Ts can direct RNA polymerase to incorporate into the mRNA nucleotides not encoded by the DNA, changing the reading frame downstream from the slippage site. We detected series of 9 or 10 As in spa13, spa33 and mxiA encoding type III secretion apparatus components. Analysis of cDNAs indicated that transcriptional slippage occurs in spa13, mxiA and spa33. Changes in the reading frame were confirmed by using plasmids carrying slippage sites in the 5' part of lacZ. Slippage is required for production of Spa13 from two overlapping reading frames and should lead to production of truncated MxiA and Spa33 proteins. Complementation of spa13 and mxiA mutants with plasmids carrying altered sites indicated that slippage in spa13 is required for assembly of the secretion apparatus and that slippage sites in spa13 and mxiA have not been selected to encode Lys residues or to produce two proteins endowed with different activities. The presence of slippage sites decreases production of Spa13 by 70%, of MxiA and Spa33 by 15% and of Spa32 (encoded downstream from spa13) by 50%. These results suggest that transcriptional slippage controls protein production by reducing the proportion of mRNA translated into functional proteins.

Identification of Escherichia coli O157:H7 genomic regions conserved in strains with a genotype associated with human infection

Beta-glucuronidase-negative, sorbitol-nonfermenting isolates of Shiga toxin-producing Escherichia coli O157 comprise part of a clone complex of related enterohemorrhagic E. coli isolates. High-resolution genotyping shows that the O157 populations have diverged into two different lineages that appear to have different ecologies. To identify genomic regions unique to the most common human-associated genotype, suppression subtractive hybridization was used to identify DNA sequences present in two clinical strains representing the human lineage I O157:H7 strains but absent from two bovine-derived lineage II strains. PCR assays were then used to test for the presence of these regions in 10 lineage I strains and 20 lineage II strains. Twelve conserved regions of genomic difference for lineage I (CRD(I)) were identified that were each present in at least seven of the lineage I strains but absent in most of the lineage II strains tested. The boundaries of the lineage I conserved regions were further delimited by PCR. Eleven of these CRD(I) were associated with E. coli Sakai S-loops 14, 16, 69, 72, 78, 82, 83, 91 to 93, 153, and 286, and the final CRD(I) was located on the pO157 virulence plasmid. Several potential virulence factors were identified within these regions, including a putative hemolysin-activating protein, an iron transport system, and several possible regulatory genes. Cluster analysis based on lineage I conserved regions showed that the presence/absence of these regions was congruent with the inferred phylogeny of the strains.

Shigellaas a Foodborne Pathogen and Current Methods for Detection in Food

Shigella, the causative agent of shigellosis or "bacillary dysentery," has been increasingly involved in foodborne outbreaks. According to the Centers for Disease Control and Prevention's Emerging Infections Program, Foodborne Diseases Active Surveillance Network (FoodNet), Shigella was the third most reported foodborne bacterial pathogen in 2002. Foods are most commonly contaminated with Shigella by an infected food handler who practices poor personal hygiene. Shigella is acid resistant, salt tolerant, and can survive at infective levels in many types of foods such as fruits and vegetables, low pH foods, prepared foods, and foods held in modified atmosphere or vacuum packaging. Survival is often increased when food is held at refrigerated temperatures. Detection methods for Shigella include conventional culture methods, immunological methods, and molecular microbiological methods. Conventional culture of Shigella in foods is often problematic due to the lack of appropriate selective media. Immunological methods for Shigella have been researched, yet there is only one commercially available test kit. Molecular microbiological methods such as PCR, oligonucleotide microarrays, and rep-PCR have also been developed for the detection and identification of Shigella. This manuscript reviews the general characteristics, prevalence, growth and survival, and methods for detection of Shigella in food.

The use of comparative genomic hybridization to characterize genome dynamics and diversity among the serotypes of Shigella

Background

Compelling evidence indicates that Shigella species, the etiologic agents of bacillary dysentery, as well as enteroinvasive Escherichia coli, are derived from multiple origins of Escherichia coli and form a single pathovar. To further understand the genome diversity and virulence evolution of Shigella, comparative genomic hybridization microarray analysis was employed to compare the gene content of E. coli K-12 with those of 43 Shigella strains from all lineages.

Results

For the 43 strains subjected to CGH microarray analyses, the common backbone of the Shigella genome was estimated to contain more than 1,900 open reading frames (ORFs), with a mean number of 726 undetectable ORFs. The mosaic distribution of absent regions indicated that insertions and/or deletions have led to the highly diversified genomes of pathogenic strains.

Conclusion

These results support the hypothesis that by gain and loss of functions, Shigella species became successful human pathogens through convergent evolution from diverse genomic backgrounds. Moreover, we also found many specific differences between different lineages, providing a window into understanding bacterial speciation and taxonomic relationships.

No single irreplaceable acidic residues in the Escherichia coli secondary multidrug transporter MdfA

The largest family of solute transporters (major facilitator superfamily [MFS]) includes proton-motive-force-driven secondary transporters. Several characterized MFS transporters utilize essential acidic residues that play a critical role in the energy-coupling mechanism during transport. Surprisingly, we show here that no single acidic residue plays an irreplaceable role in the Escherichia coli secondary multidrug transporter MdfA.

Comparison of the virulence plasmid genomes of two strains of Shigella which lost the ability to bind Congo red

We determined and analyzed the Shigella flexneri serotype 5 (pSF5) and S. dysenteriae serotype 1 (pSD1) virulence plasmid genomes. The total length of pSF5 is 136513 bp, including 165 open reading frames (ORFs). Of these ORFs, 133 were identified and 32 of those had no significant homology to proteins with known functions. The length of pSD1 is 182545 bp, including 224 ORFs, of which we identified 181. The remaining 43 ORFs were not significantly homologous to proteins with known functions. The insertion sequence (IS) elements are 53787 bp in pSF5, and 49616 bp in pSD1, which represents 39.4% and 27.1% of the genome, respectively. There are 22 IS element types in pSF5 and pSD1, among which we report ISEc8 and ISSbo6 for the first time in the Shigella virulence plasmid. Compared to pCP301, there are a large number of deleted genes and gene inversions in both pSF5 and pSD1. The ipa-mxi-spa locus in pSF5 is completely absent, and the genes related to the O-antigen biosynthesis are partially missing. In contrast, the above genes in pSD1 are integral, with the exception of virF. The whole genome analysis of the two plasmids shows that the loss of genes related to gene invasion or regulation also obliterates the ability of pPF5 and pSD1 to bind Congo red (Crb). Whether these genes determine the Crb function requires continued investigation.

Two Frameshift Products Involved in the Transposition of Bacterial Insertion Sequence IS629

IS629 is 1,310 bp in length with a pair of 25-bp imperfect inverted repeats at its termini. Two partially overlapping open reading frames, orfA and orfB, are present in IS629, and two putative translational frameshift signals, TTTTG (T4G) and AAAAT (A4T), are located near the 3'-end of orfA. With the lacZ gene as the reporter, both T4G and A4T motifs are determined to be a -1 frameshift signal. Two peptides representing the two transframe products designated OrfAB' and OrfAB, are identified by a liquid chromatography-tandem mass spectrometric approach. Results of transposition assays show that OrfAB' is the transposase and that OrfAB aids in the transposition of IS629. Pulse-chase experiments and Escherichia coli two-hybrid assays demonstrate that OrfAB binds to and stabilizes OrfAB', thus increasing the transposition activity of IS629. This is the first transposable element in the IS3 family shown to have two functional frameshifted products involved in transposition and to use a transframe product to regulate transposition.

Developing liveShigellavaccines using λ Red recombineering

Live attenuated Shigella vaccines have shown promise in inducing protective immune responses in human clinical trials and as carriers of heterologous antigens from other mucosal pathogens. In the past, construction of Shigella vaccine strains relied on classical allelic exchange systems to genetically engineer the bacterial genome. These systems require extensive in vitro engineering of long homologous sequences to create recombinant replication-defective plasmids or phage. Alternatively, the lambda red recombination system from bacteriophage facilitates recombination with as little as 40 bp of homologous DNA. The process, referred to as recombineering, typically uses an inducible lambda red operon on a temperature-sensitive plasmid and optimal transformation conditions to integrate linear antibiotic resistance cassettes flanked by homologous sequences into a bacterial genome. Recent advances in recombineering have enabled modification of genomic DNA from bacterial pathogens including Salmonella, Yersinia, enteropathogenic Escherichia coli, or enterohemorrhagic E. coli and Shigella. These advances in recombineering have been used to systematically delete virulence-associated genes from Shigella, creating a number of isogenic strains from multiple Shigella serotypes. These strains have been characterized for attenuation using both in vivo and in vitro assays. Based on this data, prototypic Shigella vaccine strains containing multiple deletions in virulence-associated genes have been generated.

Subtractive hybridization and optical mapping of the enterotoxigenic Escherichia coli H10407 chromosome: isolation of unique sequences and demonstration of significant similarity to the chromosome of E. coli K-12

Enterotoxigenic Escherichia coli (ETEC) is a primary cause of diarrhoea in infants in developing countries and in travellers to endemic regions. While several virulence genes have been identified on ETEC plasmids, little is known about the ETEC chromosome, although it is expected to share significant homology in backbone sequences with E. coli K-12. In the absence of genomic sequence information, the subtractive hybridization method and the more recently described optical mapping technique were carried out to determine the degree of genomic variation between virulent ETEC strain H10407 and the non-pathogenic E. coli K-12 strain MG1655. In one round of PCR-based suppression subtractive hybridization, 153 fragments representing sequences unique to strain H10407 were identified. blast searches indicated that few unique sequences showed homology to known pathogenicity island genes identified in related E. coli pathogens. A total of 65 fragments contained sequences that were either linked to hypothetical proteins or showed no homology to any known sequence in the database. The remaining sequences were either phage or prophage related or displayed homology to classifiable genes that function in various aspects of bacterial metabolism. The 153 unique sequences showed variable distribution across different ETEC strains including ETEC strain B7A, which is attenuated in virulence and lacked several H10407-specific sequences. Restriction-enzyme-based optical maps of strain H10407 were compared to in silico restriction maps of strain MG1655 and related E. coli pathogens. The 5.1 Mb ETEC chromosome was approximately 500 kb greater in length than the chromosome of E. coli K-12, collinear with it and indicated several discrete regions where insertions and/or deletions had occurred relative to the chromosome of strain MG1655. No major inversions, transpositions or gross rearrangements were observed on the ETEC chromosome. Based on comparisons with known genomic sequences and related optical-map-based restriction site similarity, the sequence of the H10407 chromosome is expected to demonstrate approximately 96 % identity with that of E. coli K-12.

Tempo and mode of early gene loss in endosymbiotic bacteria from insects

BACKGROUND

Understanding evolutionary processes that drive genome reduction requires determining the tempo (rate) and the mode (size and types of deletions) of gene losses. In this study, we analysed five endosymbiotic genome sequences of the gamma-proteobacteria (three different Buchnera aphidicola strains, Wigglesworthia glossinidia, Blochmannia floridanus) to test if gene loss could be driven by the selective importance of genes. We used a parsimony method to reconstruct a minimal ancestral genome of insect endosymbionts and quantified gene loss along the branches of the phylogenetic tree. To evaluate the selective or functional importance of genes, we used a parameter that measures the level of adaptive codon bias in E. coli (i.e. codon adaptive index, or CAI), and also estimates of evolutionary rates (Ka) between pairs of orthologs either in free-living bacteria or in pairs of symbionts.

RESULTS

Our results demonstrate that genes lost in the early stages of symbiosis were on average less selectively constrained than genes conserved in any of the extant symbiotic strains studied. These results also extend to more recent events of gene losses (i.e. among Buchnera strains) that still tend to concentrate on genes with low adaptive bias in E. coli and high evolutionary rates both in free-living and in symbiotic lineages. In addition, we analyzed the physical organization of gene losses for early steps of symbiosis acquisition under the hypothesis of a common origin of different symbioses. In contrast with previous findings we show that gene losses mostly occurred through loss of rather small blocks and mostly in syntenic regions between at least one of the symbionts and present-day E. coli.

CONCLUSION

At both ancient and recent stages of symbiosis evolution, gene loss was at least partially influenced by selection, highly conserved genes being retained more readily than lowly conserved genes: although losses might result from drift due to the bottlenecking of endosymbiontic populations, we demonstrated that purifying selection also acted by retaining genes of greater selective importance.

Complete genome sequence of Shigella flexneri 5b and comparison with Shigella flexneri 2a

Background

Shigella bacteria cause dysentery, which remains a significant threat to public health. Shigella flexneri is the most common species in both developing and developed countries. Five Shigella genomes have been sequenced, revealing dynamic and diverse features. To investigate the intra-species diversity of S. flexneri genomes further, we have sequenced the complete genome of S. flexneri 5b strain 8401 (abbreviated Sf8401) and compared it with S. flexneri 2a (Sf301).

Results

The Sf8401 chromosome is 4.5-Mb in size, a little smaller than that of Sf301, mainly because the former lacks the SHI-1 pathogenicity island (PAI). Compared with Sf301, there are 6 inversions and one translocation in Sf8401, which are probably mediated by insertion sequences (IS). There are clear differences in the known PAIs between these two genomes. The bacteriophage SfV segment remaining in SHI-O of Sf8401 is clearly larger than the remnants of bacteriophage SfII in Sf301. SHI-1 is absent from Sf8401 but a specific related protein is found next to the pheV locus. SHI-2 is involved in one intra-replichore inversion near the origin of replication, which may change the expression of iut/iuc genes. Moreover, genes related to the glycine-betaine biosynthesis pathway are present only in Sf8401 among the known Shigella genomes.

Conclusion

Our data show that the two S. flexneri genomes are very similar, which suggests a high level of structural and functional conservation between the two serotypes. The differences reflect different selection pressures during evolution. The ancestor of S. flexneri probably acquired SHI-1 and SHI-2 before SHI-O was integrated and the serotypes diverged. SHI-1 was subsequently deleted from the S. flexneri 5b genome by recombination, but stabilized in the S. flexneri 2a genome. These events may have contributed to the differences in pathogenicity and epidemicity between the two serotypes of S. flexneri.

Application of comparative phylogenomics to study the evolution of Yersinia enterocolitica and to identify genetic differences relating to pathogenicity

Yersinia enterocolitica, an important cause of human gastroenteritis generally caused by the consumption of livestock, has traditionally been categorized into three groups with respect to pathogenicity, i.e., nonpathogenic (biotype 1A), low pathogenicity (biotypes 2 to 5), and highly pathogenic (biotype 1B). However, genetic differences that explain variation in pathogenesis and whether different biotypes are associated with specific nonhuman hosts are largely unknown. In this study, we applied comparative phylogenomics (whole-genome comparisons of microbes with DNA microarrays combined with Bayesian phylogenies) to investigate a diverse collection of 94 strains of Y. enterocolitica consisting of 35 human, 35 pig, 15 sheep, and 9 cattle isolates from nonpathogenic, low-pathogenicity, and highly pathogenic biotypes. Analysis confirmed three distinct statistically supported clusters composed of a nonpathogenic clade, a low-pathogenicity clade, and a highly pathogenic clade. Genetic differences revealed 125 predicted coding sequences (CDSs) present in all highly pathogenic strains but absent from the other clades. These included several previously uncharacterized CDSs that may encode novel virulence determinants including a hemolysin, a metalloprotease, and a type III secretion effector protein. Additionally, 27 CDSs were identified which were present in all 47 low-pathogenicity strains and Y. enterocolitica 8081 but absent from all nonpathogenic 1A isolates. Analysis of the core gene set for Y. enterocolitica revealed that 20.8% of the genes were shared by all of the strains, confirming this species as highly heterogeneous, adding to the case for the existence of three subspecies of Y. enterocolitica. Further analysis revealed that Y. enterocolitica does not cluster according to source (host).

A SitABCD homologue from an avian pathogenic Escherichia coli strain mediates transport of iron and manganese and resistance to hydrogen peroxide

An operon encoding a member of the family of ATP-binding cassette (ABC) divalent metal ion transporters, homologous to Salmonella enterica SitABCD, has been identified in the avian pathogenic Escherichia coli (APEC) strain chi7122. The sitABCD genes were located on the virulence plasmid pAPEC-1, and were highly similar at the nucleotide level to the chromosomally encoded sitABCD genes present in Shigella spp. A cloned copy of sitABCD conferred increased growth upon a siderophore-deficient E. coli strain grown in nutrient broth supplemented with the chelator 2,2'-dipyridyl. Ion rescue demonstrated that Sit-mediated growth promotion of this strain was due to the transport of iron. SitABCD mediated increased transport of both iron and manganese as demonstrated by uptake of 55Fe, 59Fe or 54Mn in E. coli K-12 strains deficient for the transport of iron (aroB feoB) and manganese (mntH) respectively. Isotope uptake and transport inhibition studies showed that in the iron transport deficient strain, SitABCD demonstrated a greater affinity for iron than for manganese, and SitABCD-mediated transport was higher for ferrous iron, whereas in the manganese transport deficient strain, SitABCD demonstrated greater affinity for manganese than for iron. Introduction of the APEC sitABCD genes into an E. coli K-12 mntH mutant also conferred increased resistance to the bactericidal effects of hydrogen peroxide. APEC strain chi7122 derivatives lacking either a functional SitABCD or a functional MntH transport system were as resistant to hydrogen peroxide as the wild-type strain, whereas a Deltasit DeltamntH double mutant was more sensitive to hydrogen peroxide. Overall, the results demonstrate that in E. coli SitABCD represents a manganese and iron transporter that, in combination with other ion transport systems, may contribute to acquisition of iron and manganese, and resistance to oxidative stress.

Construction, detection and microarray analysis on Shigella dysenteriae A1 IroN, ShuA single, double mutants

In this study, we constructed single mutants MTS-1, MTS-2 of IroN and ShuA gene and double mutant MTS of them in Shigella dysenteriae A1 strain 51197 by insert and absence. The functional detection of every mutant was performed at the level of culture medium and cell experiment. The gene expression profiles of the mutants and the wild-type strains under iron-enriched and iron-limited conditions were analyzed by the SD51197 whole genomic microarray. The results showed that all the mutants grew obviously less well than the wild-type strains in L broth appending iron chelator DIP. The addition of iron to the cultures can stimulate the growth of mutants back to wild-type levels. In either the experiments on the ability of intracellular multiplication or the cell-to-cell spread in HeLa and U937 cell lines, mutants showed no obvious change in virulence compared with the parental strain SD51197. However when DIP was added to the cultured HeLa cells, the ability of intracellular multiplication of MTS-1, MTS-2, MTS has reduced about 23.4%, 25.2%, 43.6% respectively. The analysis of expression profiles under the iron-limited condition showed that the mutants were more sensitive for the changes of iron deficiency than the wild-type strains, many genes have been altered. Up-regulated genes mainly involved genes of transcription, coenzyme metabolism, amino acid transport and metabolism, and unknown functional genes, while down-regulated genes mainly involved genes of energy and carbohydrate metabolism and unknown function genes; the expression levels of known iron-transport associated genes generally showed up-regulated. The results demonstrated that iron-transport associated genes IroN, ShuA were likely to have some effects on the virulence and growth of S. dysenteriae.

Archaeology and evolution of transfer RNA genes in the Escherichia coli genome

Transfer RNA genes tend to be presented in multiple copies in the genomes of most organisms, from bacteria to eukaryotes. The evolution and genomic structure of tRNA genes has been a somewhat neglected area of molecular evolution. Escherichia coli, the first phylogenetic species for which more than two different strains have been sequenced, provides an invaluable framework to study the evolution of tRNA genes. In this work, a detailed analysis of the tRNA structure of the genomes of Escherichia coli strains K12, CFT073, and O157:H7, Shigella flexneri 2a 301, and Salmonella typhimurium LT2 was carried out. A phylogenetic analysis of these organisms was completed, and an archaeological map depicting the main events in the evolution of tRNA genes was drawn. It is shown that duplications, deletions, and horizontal gene transfers are the main factors driving tRNA evolution in these genomes. On average, 0.64 tRNA insertions/duplications occur every million years (Myr) per genome per lineage, while deletions occur at the slower rate of 0.30 per million years per genome per lineage. This work provides a first genomic glance at the problem of tRNA evolution as a repetitive process, and the relationship of this mechanism to genome evolution and codon usage is discussed.

Genome-wide detection and analysis of homologous recombination among sequenced strains of Escherichia coli

BACKGROUND

Comparisons of complete bacterial genomes reveal evidence of lateral transfer of DNA across otherwise clonally diverging lineages. Some lateral transfer events result in acquisition of novel genomic segments and are easily detected through genome comparison. Other more subtle lateral transfers involve homologous recombination events that result in substitution of alleles within conserved genomic regions. This type of event is observed infrequently among distantly related organisms. It is reported to be more common within species, but the frequency has been difficult to quantify since the sequences under comparison tend to have relatively few polymorphic sites.

RESULTS

Here we report a genome-wide assessment of homologous recombination among a collection of six complete Escherichia coli and Shigella flexneri genome sequences. We construct a whole-genome multiple alignment and identify clusters of polymorphic sites that exhibit atypical patterns of nucleotide substitution using a random walk-based method. The analysis reveals one large segment (approximately 100 kb) and 186 smaller clusters of single base pair differences that suggest lateral exchange between lineages. These clusters include portions of 10% of the 3,100 genes conserved in six genomes. Statistical analysis of the functional roles of these genes reveals that several classes of genes are over-represented, including those involved in recombination, transport and motility.

CONCLUSION

We demonstrate that intraspecific recombination in E. coli is much more common than previously appreciated and may show a bias for certain types of genes. The described method provides high-specificity, conservative inference of past recombination events.

A selC-associated genomic island of the extraintestinal avian pathogenic Escherichia coli strain BEN2908 is involved in carbohydrate uptake and virulence

The complete nucleotide sequence and genetic organization of a new genomic island (AGI-3) isolated from the extraintestinal avian pathogenic Escherichia coli strain BEN2908 is reported. This 49,600-bp island is inserted at the selC locus and contains putative mobile genetic elements such as a phage-related integrase gene, transposase genes, and direct repeats. AGI-3 shows a mosaic structure of five modules. Some of these modules are present in other E. coli strains and in other pathogenic bacterial species. The gene cluster aec-35 to aec-37 of module 1 encodes proteins associated with carbohydrates assimilation such as a major facilitator superfamily transporter (Aec-36), a glycosidase (Aec-37), and a putative transcriptional regulator of the LacI family (Aec-35). The aec-35 to aec-37 cluster was found in 11.6% of the tested pathogenic and nonpathogenic E. coli strains. When present, the aec-35 to aec-37 cluster is strongly associated with the selC locus (97%). Deletion of the aec-35-aec-37 region affects the assimilation of seven carbohydrates, decreases the growth rate of the strain in minimal medium containing galacturonate or trehalose, and attenuates the virulence of E. coli BEN2908 for chickens.

The nature and dynamics of bacterial genomes

Though generally small and gene rich, bacterial genomes are constantly subjected to both mutational and population-level processes that operate to increase amounts of functionless DNA. As a result, the coding potential of bacterial genomes can be substantially lower than originally predicted. Whereas only a single pseudogene was included in the original annotation of the bacterium Escherichia coli, we estimate that this genome harbors hundreds of inactivated and otherwise functionless genes. Such regions will never yield a detectable phenotype, but their identification is vital to efforts to elucidate the biological role of all the proteins within the cell.

Identification of genes subject to positive selection in uropathogenic strains of Escherichia coli: a comparative genomics approach

Escherichia coli is a model laboratory bacterium, a species that is widely distributed in the environment, as well as a mutualist and pathogen in its human hosts. As such, E. coli represents an attractive organism to study how environment impacts microbial genome structure and function. Uropathogenic E. coli (UPEC) must adapt to life in several microbial communities in the human body, and has a complex life cycle in the bladder when it causes acute or recurrent urinary tract infection (UTI). Several studies designed to identify virulence factors have focused on genes that are uniquely represented in UPEC strains, whereas the role of genes that are common to all E. coli has received much less attention. Here we describe the complete 5,065,741-bp genome sequence of a UPEC strain recovered from a patient with an acute bladder infection and compare it with six other finished E. coli genome sequences. We searched 3,470 ortholog sets for genes that are under positive selection only in UPEC strains. Our maximum likelihood-based analysis yielded 29 genes involved in various aspects of cell surface structure, DNA metabolism, nutrient acquisition, and UTI. These results were validated by resequencing a subset of the 29 genes in a panel of 50 urinary, periurethral, and rectal E. coli isolates from patients with UTI. These studies outline a computational approach that may be broadly applicable for studying strain-specific adaptation and pathogenesis in other bacteria.

Differential gene expression of pathogens inside infected hosts

DNA microarray is a useful technology for studying differential gene expression in the context of microbe-host interactions. This review concentrates on recent findings of the survival strategies of three intracellular pathogens: Shigella flexneri, Salmonella enterica serovar Typhimurium and Mycobacterium tuberculosis.

Role of the novel OprD family of porins in nutrient uptake in Pseudomonas aeruginosa

To circumvent the permeability barrier of its outer membrane, Pseudomonas aeruginosa has evolved a series of specific porins. These channels have binding sites for related classes of molecules that facilitate uptake under nutrient-limited conditions. Here, we report on the identification of a 19-member family of porins similar to the basic-amino-acid-specific porin OprD. The members of this family fell into one of two phylogenetically distinct clusters, one bearing high similarity to OprD and the other bearing most similarity to the putative phenylacetic acid uptake porin PhaK of Pseudomonas putida. Analysis of the genome context, operon arrangement, and regulation of the PhaK-like porin OpdK indicated that it might be involved in vanillate uptake. This result was confirmed by demonstrating that an opdK mutant had a deficiency in the ability to grow on vanillate as a carbon source. To extrapolate these data to other paralogues within this family, the substrate specificities of 6 of the 17 remaining OprD homologues were inferred using an approach similar to that used with opdK. The specificities determined were as follows: OpdP, glycine-glutamate; OpdC, histidine; OpdB, proline; OpdT, tyrosine; OpdH, cis-aconitate; and OpdO, pyroglutamate. Thus, members of the OprD subfamily took up amino acids and related molecules, and those characterized members most similar to PhaK were responsible for the uptake of a diverse array of organic acids. These results imply that there is a functional basis for the phylogenetic clustering of these proteins and provide a framework for studying OprD homologues in other organisms.

Genomic comparisons among gamma-proteobacteria

Predicted highly expressed (PHX) genes are compared for 16 gamma-proteobacteria and their similarities and differences are interpreted with respect to known or predicted physiological characteristics of the organisms. Predicted highly expressed genes often reflect the organism's predominant lifestyle, habitat, nutrition sources and metabolic propensities. This technique allows to predict principal metabolic activities of the microorganisms operating in their natural habitats. Among our findings is an unusually high number of PHX enzymes acting in cell wall biosynthesis, amino acid biosynthesis and replication in the ant endosymbiont Blochmannia floridanus. We ascribe the abundance of these PHX genes to specific aspects of the relationship between the bacterium and its host. Xanthomonas campestris is unique with a very high number of PHX genes acting in flagellum biosynthesis, which may play a special role during its pathogenicity. Shewanella oneidensis possesses three protein complexes which all can function as complex I in the respiratory chain but only the Na(+)-transporting NADH:ubiquinone oxidoreductase nqr-2 operon is PHX. The PHX genes of Vibrio parahaemolyticus are consistent with the microorganism's adaptation to extremely fast growth rates. Comparative analysis of PHX genes from complex environmental genomic sequences as well as from uncultured pathogenic microbes can provide a novel, useful tool to predict global flux of matter and key intermediates.

ShiBASE: an integrated database for comparative genomics of Shigella

Among the major enteric bacterial pathogens, Shigella is found to display extreme genome diversity and dynamics, which imposes a challenge in comparative genomic studies. To facilitate further studies in this area, we have constructed an integrated online database, ShiBASE (),which contains Shigella genomic sequences of four species and additional comparative genomic hybridization (CGH) data of 43 serotypes. ShiBASE offers online comparative analysis on DNA sequences, gene orders, metabolic pathways and virulence factors. In addition, ShiBASE has a newly developed online comparative visualization service, Shi-align, which enables the alignment of any query sequence with the reference genome sequences.

Comparative analysis of tandem T7-like promoter containing regions in enterobacterial genomes reveals a novel group of genetic islands

Based on molecular information theory, 10 T7-like promoter models were built for the T7 group of phages and used to scan their host genomes and closely related genomes. 38 genomes were scanned and 12 clusters of tandem promoters were identified in nine enteropathogens. Comparative analysis of these tandem promoter-bearing regions reveals that they are similar to each other, forming prophage-like islands of 4-13 kb. Each island appears to contain two or three tandem T7-like promoters within a stretch of 150-620 bases, but there are no corresponding RNA polymerase (RNAP) genes. The promoters would transcribe two to five putative phage-related proteins, but none of these resemble known phage structural proteins. An integrase belonging to the Int family of site-specific recombinases is encoded upstream of the tandem promoters. A direct repeat of 17-24 bases was found on the ends of all 12 islands. Comparative analysis of the islands shows that these islands appear to have recombined with each other. These results suggest that the islands could encode a group of satellite phages. Activation and function of the islands may depend on transcription by a T7-like RNAP after infection by a T7-like phage or foreign DNA that encodes a T7-like RNAP.

ISPa20 advances the individual evolution of Pseudomonas aeruginosa clone C subclone C13 strains isolated from cystic fibrosis patients by insertional mutagenesis and genomic rearrangements

Pseudomonas aeruginosa clone C strains, which chronically colonize the lungs of cystic fibrosis patients reorganize their genome structure. In this study, a novel member of the IS3 subfamily of IS elements, ISPa20, was detected which was specific for clone C subclone C13 strains. ISPa20, which was present in high copy number, mediated events of genomic reorganization. ISPa20 was inserted into P. aeruginosa backbone genes leading to adaptation to the cystic fibrosis lung habitat and into DNA acquired through horizontal gene transfer. Further on, large chromosomal inversions were mediated by ISPa20. In contrast to strains of other subclonal linages high rates of genomic rearrangements of subclone C13 strains were observed in vitro. The acquisition of mobile elements by P. aeruginosa clone C strains in the lungs of cystic fibrosis patients supports the chronic colonization by insertional mutagenesis and chromosome restructuring leading to microevolution within clone C that reflects macroevolution observed on the species level.

Genomic compositions and phylogenetic analysis of Shigella boydii subgroup

Comparative Genomic Hybridization (CGH) microarray analysis was used to compare the genomic compositions of all eighteen Shigella boydii serotype representative strains. The results indicated the genomic "backbone" of this subgroup contained 2552 ORFs homologous to nonpathogenic E. coli K12. Compared with the genome of K12199 ORFs were found to be absent in all S. boydii serotype representatives, including mainly outer membrane protein genes and O-antigen biosynthesis genes. Yet the specific ORFs of S. boydii subgroup contained basically bacteriophage genes and the function unknown (FUN) genes. Some iron metabolism, transport and type II secretion system related genes were found in most representative strains. According to the CGH phylogenetic analysis, the eighteen S. boydii serotype representatives were divided into four groups, in which serotype C13 strain was remarkably distinguished from the other serotype strains. This grouping result corresponded to the distribution of some metabolism related genes. Furthermore, the analysis of genome backbone genes, specific genes, and the phylogenetic trees allowed us to discover the evolution laws of S. boydii and to find out important clues to pathogenesis research, vaccination and the therapeutic medicine development.

A novel strategy for the identification of genomic islands by comparative analysis of the contents and contexts of tRNA sites in closely related bacteria

We devised software tools to systematically investigate the contents and contexts of bacterial tRNA and tmRNA genes, which are known insertion hotspots for genomic islands (GIs). The strategy, based on MAUVE-facilitated multigenome comparisons, was used to examine 87 Escherichia coli MG1655 tRNA and tmRNA genes and their orthologues in E.coli EDL933, E.coli CFT073 and Shigella flexneri Sf301. Our approach identified 49 GIs occupying ∼1.7 Mb that mapped to 18 tRNA genes, missing 2 but identifying a further 30 GIs as compared with Islander [Y. Mantri and K. P. Williams (2004), Nucleic Acids Res., 32, D55–D58]. All these GIs had many strain-specific CDS, anomalous GC contents and/or significant dinucleotide biases, consistent with foreign origins. Our analysis demonstrated marked conservation of sequences flanking both empty tRNA sites and tRNA-associated GIs across all four genomes. Remarkably, there were only 2 upstream and 5 downstream deletions adjacent to the 328 loci investigated. In silico PCR analysis based on conserved flanking regions was also used to interrogate hotspots in another eight completely or partially sequenced E.coli and Shigella genomes. The tools developed are ideal for the analysis of other bacterial species and will lead to in silico and experimental discovery of new genomic islands.

Fur regulates acid resistance in Shigella flexneri via RyhB and ydeP

Shigella flexneri requires iron for survival, and the genes for iron uptake and homeostasis are regulated by the Fur protein. Microarrays were used to identify genes regulated by Fur and to study the physiological effects of iron availability in S. flexneri. These assays showed that the expression of genes involved in iron acquisition and acid response was induced by low-iron availability and by inactivation of fur. A fur null mutant was acid sensitive in media at pH 2.5, and acid sensitivity was also observed in the wild-type strain grown under iron-limiting conditions. Acid resistance of the fur mutant in minimal medium was restored by addition of glutamate during acid challenge, indicating that the glutamate-dependent acid resistance system was not defective. Inactivation of ryhB, a small regulatory RNA whose expression is repressed by Fur, restored acid resistance in the fur mutant, while overexpressing ryhB increased acid sensitivity in the wild-type strain. RyhB-regulated genes were identified by microarray analysis. The expression of one of the RyhB-repressed genes, ydeP, which encodes a putative oxidoreductase, suppressed acid sensitivity in the fur mutant. Furthermore, an S. flexneri ydeP mutant was defective for both glutamate-independent and glutamate-dependent acid resistance. The repression of ydeP by RyhB may be indirect, as real time polymerase chain reaction (PCR) experiments indicated that RyhB negatively regulates evgA, which encodes an activator of ydeP. These results demonstrate that the acid sensitivity defect of the S. flexneri fur mutant is due to repression of ydeP by RyhB, most likely via repression of evgA.

Host specificity of septicemic Escherichia coli: human and avian pathogens

Extraintestinal pathogenic Escherichia coli (ExPEC) strains are the cause of a diverse spectrum of invasive human and animal infections, often leading to septicemia. ExPEC strains contain virulence factors that enable them to survive in the host blood and tissues. Most of these virulence factors are distributed in ExPEC strains in a host-independent fashion. Genomic analyses of these strains provide evidence for numerous recombinational events and horizontal gene transfer, as well as for a high diversity of virulence factors. In studies of human and avian septicemic strains of serotypes O2 and O78 it appears that there is a positive correlation between virulence, invasiveness and clonal origin. Yet, it is clear that clonal division in these strains, as well as distribution of virulence factors, is independent of the host and closely related clones reside in different hosts. Although the possibility exists that ExPEC strains do have a certain degree of host specificity, which is not obvious from genomic studies, it is clear that the similarity of virulence factors presents a significant zoonotic risk.

(Patho-)Genomics of Escherichia coli

Escherichia coli represents a versatile and diverse enterobacterial species which can be subdivided into (i) nonpathogenic, commensal, (ii) intestinal pathogenic and (iii) extraintestinal pathogenic strains. This classification is mainly based on the presence or absence of DNA regions which are frequently associated with certain pathotypes. In most cases, this genetic information has been horizontally acquired and belongs to the flexible E. coli genome, such as plasmids, bacteriophages and genomic islands. These genomic regions contribute to the rapid evolution of E. coli variants as they are frequently subject to rearrangements, excision and transfer as well as further acquisition of additional DNA thus contributing to the creation of new (pathogenic) variants. Genetic diversity and genome plasticity of E. coli has been underestimated. The accumulating amount of sequence information generated in the era of "genomics" helps to increase our understanding of factors and mechanisms that are involved in diversification of this bacterial species as well as in those that may direct host specificity.

Virulence factors of septicemic Escherichia coli strains

Extraintestinal pathogenic Escherichia coli strains (ExPEC) are the cause of a diverse spectrum of invasive human and animal infections, often leading to septicemia. This review deals with the virulence genes of septicemic ExPEC strains. We discuss the meaning of a virulence gene and survey the genomic, genetic and physiological studies on these strains. Apparently, there are a few virulence factors, which are conserved in the septicemic strains, implying that they are essential for the infection. For the other virulence-related genes a high level of diversity is observed, demonstrating that all stages of the infection can be mediated by a number of alternative virulence factors. The variable profile of virulence genes in septicemic E. coli strains, as well as a prevalence of mobility-related sequences point out the existence of a "mix and match" combinatorial system.

A comparative categorization of protein function encoded in bacterial or archeal genomic islands

Genomes of prokaryotes harbor genomic islands (GIs), which are frequently acquired via horizontal gene transfer (HGT). Here I present an analysis of GIs with respect to gene-encoded functions. GIs were identified by statistical analysis of codon usage and clustering. Genes classified as putatively alien (pA) or putatively native (pN) were categorized according to the COG database. Among pA and pN genes, the distribution of COG functions and classes were studied for different groupings of prokaryotes. Groups were formed according to taxonomical relation or habitats. In all groups, genes related to class L (replication, recombination, and repair) were statistically significantly overrepresented in GIs. GIs of bacteria and archaea showed a distinct pattern of preferences. In archeal GIs, genes belonging to COG class M (cell wall/membrane/envelope biogenesis) or Q (secondary metabolites biosynthesis, transport, and catabolism) were more frequent. In bacterial GIs, genes of classes U (intracellular trafficking, secretion, and vesicular transport), N (cell motility), and V (defense mechanisms) were predominant. Underrepresentation was strongest for genes belonging to class J (translation, ribosomal structure, and biogenesis). Among single COG functions overrepresented in GIs were transferases and transporters. In both superkingdoms, HGT enhances genomic content by meeting demands that are independent of the studied habitats. These findings are in agreement with the complexity theory, which predicts the preferential import of operational genes. However, only specific subsets of operational genes were enriched in GIs. Modification of the cell envelope, cell motility, secretion, and protection of cellular DNA are major issues in HGT.

Genomic comparison of plant pathogenic and nonpathogenic Serratia marcescens strains by suppressive subtractive hybridization

Cucurbit yellow vine disease (CYVD) is caused by disease-associated Serratia marcescens strains that have phenotypes significantly different from those of nonphytopathogenic strains. To identify the genetic differences responsible for pathogenicity-related phenotypes, we used a suppressive subtractive hybridization (SSH) strategy. S. marcescens strain Z01-A, isolated from CYVD-affected zucchini, was used as the tester, whereas rice endophytic S. marcescens strain R02-A (IRBG 502) was used as the driver. SSH revealed 48 sequences, ranging from 200 to 700 bp, that were present in Z01-A but absent in R02-A. Sequence analysis showed that a large proportion of these sequences resembled genes involved in synthesis of surface structures. By construction of a fosmid library, followed by colony hybridization, selection, and DNA sequencing, a phage gene cluster and a genome island containing a fimbrial-gene cluster were identified. Arrayed dot hybridization showed that the conservation of subtracted sequences among CYVD pathogenic and nonpathogenic S. marcescens strains varied. Thirty-four sequences were present only in pathogenic strains. Primers were designed based on one Z01-A-specific sequence, A79, and used in a multiplex PCR to discriminate between S. marcescens strains causing CYVD and those from other ecological niches.

Replication termination in Escherichia coli: structure and antihelicase activity of the Tus-Ter complex

The arrest of DNA replication in Escherichia coli is triggered by the encounter of a replisome with a Tus protein-Ter DNA complex. A replication fork can pass through a Tus-Ter complex when traveling in one direction but not the other, and the chromosomal Ter sites are oriented so replication forks can enter, but not exit, the terminus region. The Tus-Ter complex acts by blocking the action of the replicative DnaB helicase, but details of the mechanism are uncertain. One proposed mechanism involves a specific interaction between Tus-Ter and the helicase that prevents further DNA unwinding, while another is that the Tus-Ter complex itself is sufficient to block the helicase in a polar manner, without the need for specific protein-protein interactions. This review integrates three decades of experimental information on the action of the Tus-Ter complex with information available from the Tus-TerA crystal structure. We conclude that while it is possible to explain polar fork arrest by a mechanism involving only the Tus-Ter interaction, there are also strong indications of a role for specific Tus-DnaB interactions. The evidence suggests, therefore, that the termination system is more subtle and complex than may have been assumed. We describe some further experiments and insights that may assist in unraveling the details of this fascinating process.

Inference and analysis of the relative stability of bacterial chromosomes

The stability of genomes is highly variable, both in terms of gene content and gene order. Here I calibrate the loss of gene order conservation (GOC) through time by fitting a simple probabilistic model on pairwise comparisons involving 126 bacterial genomes. The model computes the probability of separation of pairs of contiguous genes per unit of time and fits the data better than previous ones while allowing a mechanistic interpretation for the loss of GOC with time. Although the information on operons is not used in the model, I observe, as expected, that most highly conserved pairs of genes are indeed within operons. However, even the other pairs are much more conserved than expected given the observed experimental rearrangement rates. After 500 Myr, about 50% of the originally contiguous orthologues remain so in the average genome. Hence, the large majority of rearrangements must be deleterious and random genome rearrangements are unlikely to provide for positively selected structural changes. I then use the deviations from the model to define an intrinsic measure of genome stability that allowed the comparison of distantly related genomes and the inference of ancestral states. This shows that clades differ in genome stability, with cyanobacteria being the least stable and gamma-proteobacteria the most stable. Without correction for phylogeny, free-living bacteria are the least stable group of genomes, followed by pathogens, and then endomutualists. However, after correction for phylogenetic inertia (or the removal of cyanobacteria from the analysis), there is no significant association between genome stability and lifestyle or genome size. Hence, although this method has allowed uncovering some of mechanisms leading to rearrangements, we still ignore the forces that differentially shape selection upon genome stability in different species.