Complete genome sequence and comparative genomics of Shigella flexneri serotype 2a strain 2457T

Draft genome sequences of the type strains of Shigella flexneri held at Public Health England: comparison of classical phenotypic and novel molecular assays with whole genome sequence

Background

Public Health England (PHE) holds a collection of Shigella flexneri Type strains isolated between 1949 and 1972 representing 15 established serotypes and one provisional type, E1037. In this study, the genomes of all 16 PHE Type strains were sequenced using the Illumina HiSeq platform. The relationship between core genome phylogeny and serotype was examined.

Results

The most common target gene for the detection of Shigella species in clinical PCR assays, ipaH, was detected in all genomes. The type-specific target genes were correctly identified in each genome sequence. In contrast to the S. flexneri in serotype 5 strain described by Sun et al. (2012), the two PHE serotype 5 Type strains possessed an additional oac gene and were differentiated by the presence (serotype 5b) or absence (serotype 5a) of gtrX. The somatic antigen structure and phylogenetic relationship were broadly congruent for strains expressing serotype specific antigens III, IV and V, but not for those expressing I and II. The whole genome phylogenies of the 15 isolates sequenced showed that the serotype 6 Type Strain was phylogenetically distinct from the other S. flexneri serotypes sequenced. The provisional serotype E1037 fell within the serotype 4 clade, being most closely related to the Serotype 4a Type Strain.

Conclusions

The S. flexneri genome sequences were used to evaluate phylogenetic relationships between Type strains and validate genotypic and phenotypic assays. The analysis confirmed that the PHE S. flexneri Type strains are phenotypically and genotypically distinct. Novel variants will continue to be added to this archive.

Safety in numbers: multiple occurrences of highly similar homologs among Azotobacter vinelandii carbohydrate metabolism proteins probably confer adaptive benefits

BACKGROUND

Gene duplication and horizontal gene transfer are common processes in bacterial and archaeal genomes, and are generally assumed to result in either diversification or loss of the redundant gene copies. However, a recent analysis of the genome of the soil bacterium Azotobacter vinelandii DJ revealed an abundance of highly similar homologs among carbohydrate metabolism genes. In many cases these multiple genes did not appear to be the result of recent duplications, or to function only as a means of stimulating expression by increasing gene dosage, as the homologs were located in varying functional genetic contexts. Based on these initial findings we here report in-depth bioinformatic analyses focusing specifically on highly similar intra-genome homologs, or synologs, among carbohydrate metabolism genes, as well as an analysis of the general occurrence of very similar synologs in prokaryotes.

RESULTS

Approximately 900 bacterial and archaeal genomes were analysed for the occurrence of synologs, both in general and among carbohydrate metabolism genes specifically. This showed that large numbers of highly similar synologs among carbohydrate metabolism genes are very rare in bacterial and archaeal genomes, and that the A. vinelandii DJ genome contains an unusually large amount of such synologs. The majority of these synologs were found to be non-tandemly organized and localized in varying but metabolically relevant genomic contexts. The same observation was made for other genomes harbouring high levels of such synologs. It was also shown that highly similar synologs generally constitute a very small fraction of the protein-coding genes in prokaryotic genomes. The overall synolog fraction of the A. vinelandii DJ genome was well above the data set average, but not nearly as remarkable as the levels observed when only carbohydrate metabolism synologs were considered.

CONCLUSIONS

Large numbers of highly similar synologs are rare in bacterial and archaeal genomes, both in general and among carbohydrate metabolism genes. However, A. vinelandii and several other soil bacteria harbour large numbers of highly similar carbohydrate metabolism synologs which seem not to result from recent duplication or transfer events. These genes may confer adaptive benefits with respect to certain lifestyles and environmental factors, most likely due to increased regulatory flexibility and/or increased gene dosage.

Genomic portrait of the evolution and epidemic spread of a recently emerged multidrug-resistant Shigella flexneri clone in China

Shigella flexneri is the major cause of shigellosis in developing countries. A new S. flexneri serotype, Xv, appeared in 2000 and replaced serotype 2a as the most prevalent serotype in China. Serotype Xv is a variant of serotype X, with phosphoethanolamine modification of its O antigen mediated by a plasmid that contained the opt gene. Serotype Xv isolates belong to sequence type 91 (ST91). In this study, whole-genome sequencing of 59 S. flexneri isolates of 14 serotypes (serotypes 1 to 4, Y, Yv, X, and Xv) indicated that ST91 arose around 1993 by acquiring multidrug resistance (MDR) and spread across China within a decade. A comparative analysis of the chromosome and opt-carrying plasmid pSFXv_2 revealed independent origins of 3 serotype Xv clusters in China, with different divergence times. Using 18 cluster-dividing single-nucleotide polymorphisms (SNPs), SNP typing divided 380 isolates from 3 provinces (Henan, Gansu, and Anhui) into 5 SNP genotypes (SGs). One SG predominated in each province, but substantial interregional spread of SGs was also evident. These findings suggest that MDR is the key selective pressure for the emergence of the S. flexneri epidemic clone and that Shigella epidemics in China were caused by a combination of local expansion and interregional spread of serotype Xv.

Genome degeneration and adaptation in a nascent stage of symbiosis

Symbiotic associations between animals and microbes are ubiquitous in nature, with an estimated 15% of all insect species harboring intracellular bacterial symbionts. Most bacterial symbionts share many genomic features including small genomes, nucleotide composition bias, high coding density, and a paucity of mobile DNA, consistent with long-term host association. In this study, we focus on the early stages of genome degeneration in a recently derived insect-bacterial mutualistic intracellular association. We present the complete genome sequence and annotation of Sitophilus oryzae primary endosymbiont (SOPE). We also present the finished genome sequence and annotation of strain HS, a close free-living relative of SOPE and other insect symbionts of the Sodalis-allied clade, whose gene inventory is expected to closely resemble the putative ancestor of this group. Structural, functional, and evolutionary analyses indicate that SOPE has undergone extensive adaptation toward an insect-associated lifestyle in a very short time period. The genome of SOPE is large in size when compared with many ancient bacterial symbionts; however, almost half of the protein-coding genes in SOPE are pseudogenes. There is also evidence for relaxed selection on the remaining intact protein-coding genes. Comparative analyses of the whole-genome sequence of strain HS and SOPE highlight numerous genomic rearrangements, duplications, and deletions facilitated by a recent expansion of insertions sequence elements, some of which appear to have catalyzed adaptive changes. Functional metabolic predictions suggest that SOPE has lost the ability to synthesize several essential amino acids and vitamins. Analyses of the bacterial cell envelope and genes encoding secretion systems suggest that these structures and elements have become simplified in the transition to a mutualistic association.

Lineage-specific distribution of insertion sequence excision enhancer in enterotoxigenic Escherichia coli isolated from swine

Insertion sequences (ISs) are the simplest transposable elements and are widely distributed in bacteria; however, they also play important roles in genome evolution. We recently identified a protein called IS excision enhancer (IEE) in enterohemorrhagic Escherichia coli (EHEC) O157. IEE promotes the excision of IS elements belonging to the IS3 family, such as IS629, as well as several other families. IEE-mediated IS excision generates various genomic deletions that lead to the diversification of the bacterial genome. IEE has been found in a broad range of bacterial species; however, among sequenced E. coli strains, IEE is primarily found in EHEC isolates. In this study, we investigated non-EHEC pathogenic E. coli strains isolated from domestic animals and found that IEE is distributed in specific lineages of enterotoxigenic E. coli (ETEC) strains of serotypes O139 or O149 isolated from swine. The iee gene is located within integrative elements that are similar to SpLE1 of EHEC O157. All iee-positive ETEC lineages also contained multiple copies of IS629, a preferred substrate of IEE, and their genomic locations varied significantly between strains, as observed in O157. These data suggest that IEE may have been transferred among EHEC and ETEC in swine via SpLE1 or SpLE1-like integrative elements. In addition, IS629 is actively moving in the ETEC O139 and O149 genomes and, as in EHEC O157, is promoting the diversification of these genomes in combination with IEE.

Shigella in Brazilian children with acute diarrhoea: prevalence, antimicrobial resistance and virulence genes

Diarrhoeal disease is still considered a major cause of morbidity and mortality among children. Among diarrhoeagenic agents, Shigella should be highlighted due to its prevalence and the severity of the associated disease. Here, we assessed Shigella prevalence, drug susceptibility and virulence factors. Faeces from 157 children with diarrhoea who sought treatment at the Children's Hospital João Paulo II, a reference children´s hospital in Belo Horizonte, state of Minas Gerais, Brazil, were cultured and drug susceptibility of the Shigella isolates was determined by the disk diffusion technique. Shigella virulence markers were identified by polymerase chain reaction. The bacterium was recovered from 10.8% of the children (88.2% Shigella sonnei). The ipaH, iuc, sen and ial genes were detected in strains isolated from all shigellosis patients; set1A was only detected in Shigella flexneri. Additionally, patients were infected by Shigella strains of different ial, sat, sen and set1A genotypes. Compared to previous studies, we observed a marked shift in the distribution of species from S. flexneri to S. sonnei and high rates of trimethoprim/sulfamethoxazole resistance.

Recent advances in understanding enteric pathogenic Escherichia coli

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

Phylogeny rather than ecology or lifestyle biases the construction of Escherichia coli-Shigella genetic exchange communities

Genetic material can be transmitted not only vertically from parent to offspring, but also laterally (horizontally) from one bacterial lineage to another. Lateral genetic transfer is non-uniform; biases in its nature or frequency construct communities of genetic exchange. These biases have been proposed to arise from phylogenetic relatedness, shared ecology and/or common lifestyle. Here, we test these hypotheses using a graph-based abstraction of inferred genetic-exchange relationships among 27 Escherichia coli and Shigella genomes. We show that although barriers to inter-phylogenetic group lateral transfer are low, E. coli and Shigella are more likely to have exchanged genetic material with close relatives. We find little evidence of bias arising from shared environment or lifestyle. More than one-third of donor-recipient pairs in our analysis show some level of fragmentary gene transfer. Thus, within the E. coli-Shigella clade, intact genes and gene fragments have been disseminated non-uniformly and at appreciable frequency, constructing communities that transgress environmental and lifestyle boundaries.

Inferring ancient metabolism using ancestral core metabolic models of enterobacteria

Background

Enterobacteriaceae diversified from an ancestral lineage ~300-500 million years ago (mya) into a wide variety of free-living and host-associated lifestyles. Nutrient availability varies across niches, and evolution of metabolic networks likely played a key role in adaptation.

Results

Here we use a paleo systems biology approach to reconstruct and model metabolic networks of ancestral nodes of the enterobacteria phylogeny to investigate metabolism of ancient microorganisms and evolution of the networks. Specifically, we identified orthologous genes across genomes of 72 free-living enterobacteria (16 genera), and constructed core metabolic networks capturing conserved components for ancestral lineages leading to E. coli/Shigella (~10 mya), E. coli/Shigella/Salmonella (~100 mya), and all enterobacteria (~300-500 mya). Using these models we analyzed the capacity for carbon, nitrogen, phosphorous, sulfur, and iron utilization in aerobic and anaerobic conditions, identified conserved and differentiating catabolic phenotypes, and validated predictions by comparison to experimental data from extant organisms.

Conclusions

This is a novel approach using quantitative ancestral models to study metabolic network evolution and may be useful for identification of new targets to control infectious diseases caused by enterobacteria.

Characterization of intracellular growth regulator icgR by utilizing transcriptomics to identify mediators of pathogenesis in Shigella flexneri

Shigella species Gram-negative bacteria which cause a diarrheal disease, known as shigellosis, by invading and destroying the colonic mucosa and inducing a robust inflammatory response. With no vaccine available, shigellosis annually kills over 600,000 children in developing countries. This study demonstrates the utility of combining high-throughput bioinformatic methods with in vitro and in vivo assays to provide new insights into pathogenesis. Comparisons of in vivo and in vitro gene expression identified genes associated with intracellular growth. Additional bioinformatics analyses identified genes that are present in S. flexneri isolates but not in the three other Shigella species. Comparison of these two analyses revealed nine genes that are differentially expressed during invasion and that are specific to S. flexneri. One gene, a DeoR family transcriptional regulator with decreased expression during invasion, was further characterized and is now designated icgR, for intracellular growth regulator. Deletion of icgR caused no difference in growth in vitro but resulted in increased intracellular replication in HCT-8 cells. Further in vitro and in vivo studies using high-throughput sequencing of RNA transcripts (RNA-seq) of an isogenic ΔicgR mutant identified 34 genes that were upregulated under both growth conditions. This combined informatics and functional approach has allowed the characterization of a gene and pathway previously unknown in Shigella pathogenesis and provides a framework for further identification of novel virulence factors and regulatory pathways.

Small molecule inhibitors of trans-translation have broad-spectrum antibiotic activity

The trans-translation pathway for protein tagging and ribosome release plays a critical role for viability and virulence in a wide range of pathogens but is not found in animals. To explore the use of trans-translation as a target for antibiotic development, a high-throughput screen and secondary screening assays were used to identify small molecule inhibitors of the pathway. Compounds that inhibited protein tagging and proteolysis of tagged proteins were recovered from the screen. One of the most active compounds, KKL-35, inhibited the trans-translation tagging reaction with an IC50 = 0.9 µM. KKL-35 and other compounds identified in the screen exhibited broad-spectrum antibiotic activity, validating trans-translation as a target for drug development. This unique target could play a key role in combating strains of pathogenic bacteria that are resistant to existing antibiotics.

Bacterial serine proteases secreted by the autotransporter pathway: classification, specificity, and role in virulence

Serine proteases exist in eukaryotic and prokaryotic organisms and have emerged during evolution as the most abundant and functionally diverse group. In Gram-negative bacteria, there is a growing family of high molecular weight serine proteases secreted to the external milieu by a fascinating and widely employed bacterial secretion mechanism, known as the autotransporter pathway. They were initially found in Neisseria, Shigella, and pathogenic Escherichia coli, but have now also been identified in Citrobacter rodentium, Salmonella, and Edwardsiella species. Here, we focus on proteins belonging to the serine protease autotransporter of Enterobacteriaceae (SPATEs) family. Recent findings regarding the predilection of serine proteases to host intracellular or extracellular protein-substrates involved in numerous biological functions, such as those implicated in cytoskeleton stability, autophagy or innate and adaptive immunity, have helped provide a better understanding of SPATEs' contributions in pathogenesis. Here, we discuss their classification, substrate specificity, and potential roles in pathogenesis.

Characterization and complete genome sequence of a virulent bacteriophage B4 infecting food-borne pathogenic Bacillus cereus

Bacillus cereus causes food poisoning, resulting in vomiting and diarrhea, due to production of enterotoxins. As a means of controlling this food-borne pathogen, the virulent bacteriophage B4 was isolated and characterized. Bacterial challenge assays showed that phage B4 effectively inhibited growth of members of the B. cereus group as well as B. subtilis, and growth inhibition persisted for over 20 h. One-step growth analysis also revealed the host lysis activity of phage B4, with relatively short eclipse/latent times (10/15 min) and a large burst size (>200 PFU). The complete genome of phage B4, containing a 162-kb DNA with 277 ORFs, was analyzed. The endolysin encoded by the phage B4 genome accounts for the cell lysis activity of this phage. These results suggest that phage B4 has potential as a biological agent to control B. cereus propagation.

Progress and pitfalls in Shigella vaccine research

Renewed awareness of the substantial morbidity and mortality that Shigella infection causes among young children in developing countries, combined with technological innovations in vaccinology, has led to the development of novel vaccine strategies in the past 5 years. Along with advancement of classic vaccines in clinical trials and new sophisticated measurements of immunological responses, much new data has been produced, lending promise to the potential for production of safe and effective Shigella vaccines. Herein, we review the latest progress in Shigella vaccine development within the framework of persistent obstacles.

Cross-protection provided by live Shigella mutants lacking major antigens

The immune response elicited by Shigella infections is dominated by serotype-specific antibodies recognizing the LPS O-antigens. Although a marked antibody response to invasion plasmid antigens (Ipa-s) shared by all virulent strains is also induced, the varying level of immunity elicited by natural infections is serotype-restricted. Previous vaccines have tried to mimic and achieve this serotype-specific, infection-induced immunity. As, however, the four Shigella species can express 50 different types of O-antigens, current approaches with the aim to induce a broad coverage use a mixture of the most common O-antigens combined in single vaccines. In the current study we present data on an alternative approach to generate immunity protective against multiple serotypes. Mutants lacking both major immune-determinant structures (i.e. the Ipa and O-antigens) were not only highly attenuated, but, unlike their avirulent counterparts still expressing these antigens, elicited a protective immune response to heterologous serotypes in a murine model. Evidence is provided that protection was mediated by the enhanced immunogenic potential of minor conserved antigens. Furthermore, the rough, non-invasive double mutants triggered an immune response different from that induced by the smooth, invasive strains regarding the isotype of antibodies generated. These non-invasive, rough mutants may represent promising candidates for further development into live vaccines for the prophylaxis of bacillary dysentery in areas with multiple endemic serotypes.

tmRNA is essential in Shigella flexneri

Nonstop mRNAs pose a challenge for bacteria, because translation cannot terminate efficiently without a stop codon. The trans-translation pathway resolves nonstop translation complexes by removing the nonstop mRNA, the incomplete protein, and the stalled ribosome. P1 co-transduction experiments demonstrated that tmRNA, a key component of the trans-translation pathway, is essential for viability in Shigella flexneri. tmRNA was previously shown to be dispensable in the closely related species Escherichia coli, because E. coli contains a backup system for trans-translation mediated by the alternative release factor ArfA. Genome sequence analysis showed that S. flexneri does not have a gene encoding ArfA. E. coli ArfA could suppress the requirement for tmRNA in S. flexneri, indicating that tmRNA is essential in S. flexneri because there is no functional backup system. These data suggest that resolution of nonstop translation complexes is required for most bacteria.

Isolation and genomic characterization of SfI, a serotype-converting bacteriophage of Shigella flexneri

Background

All Shigella flexneri serotypes except serotype 6 share a common O-antigen tetrasaccharide backbone and nearly all variations between serotypes are due to glucosyl and/or O-acetyl modifications of the common O unit mediated by glycosyltransferases encoded by serotype-converting bacteriophages. Several S. flexneri serotype-converting phages including SfV, SfX, Sf6 and SfII have been isolated and characterized. However, S. flexneri serotype-converting phage SfI which encodes a type I modification of serotype 1 (1a, 1b, 1c and 1d) had not yet been characterized.

Results

The SfI phage was induced and purified from a S. flexneri serotype 1a clinical strain 019. Electron microscopy showed that the SfI phage has a hexagonal head and a long contractile tail, characteristic of the members of Myoviridae family. SfI can convert serotype Y to serotype 1a and serotype X to serotype 1d, but cannot convert 10 other S. flexneri serotypes (1a, 1b, 2a, 2b, 3a, 3b, 4a, 4b, 5a, Xv) tested, suggesting that SfI has a narrow host range. Similar to other S. flexneri serotype-converting phages, SfI integrates into the tRNA-thrW gene adjacent to proA of the host chromosome when lysogenized. The complete sequence of the SfI genome was 38,389 bp, encoding 66 open reading frames and two tRNA genes. Phage SfI shares significant homology with S. flexneri phage SfV, Escherichia coli prophage e14 and lambda, and is classified into the lambdoid phage family. SfI was found to use a cos mechanism for DNA packaging similar to that of phage SfV.

Conclusions

SfI contains features of lambdoid phages and is closely related to S. flexneri phage SfV, E. coli prophage e14 and lambda. The characterization of SfI enhances our understanding of serotype conversion of S. flexneri.

Comparative analysis of the genomes of Shigella dysenteriae type 2 & type 7 isolates

BACKGROUND & OBJECTIVES

The four species of the genus Shigella, namely, S. dysenteriae , S. flexneri, S. boydii and S. sonnei cause a wide spectrum of illness from watery diarrhoea to severe dysentery. Genomes of these four species show great diversity. In this study, NotI, XbaI or I-CeuI restriction enzyme digested genomes of two Shigella dysenteriae isolates belonging to the serotypes 2 and 7 were extensively analyzed to find their relatedness, if any, with the whole genome sequenced strains of S. dysenteriae type 1 and S. flexneri type 2a.

METHODS

Pulsed-field gel electrophoresis (PFGE) technique was used to determine the diversity of Shigella genomes by rapid construction of physical maps. DNA end labelling, Southern hybridization and PCR techniques were also applied for mapping purposes.

RESULTS

The intron-coded enzyme I-CeuI cuts the bacterial genome specifically at its rrn operon. PFGE of I-CeuI digested S. dysenteriae genomes were found to carry seven rrn operons. However, I-CeuI profiles showed distinct restriction fragment polymorphism (RFLP) between the isolates as well as with the whole genome sequenced isolates. Further studies revealed that the genome sizes and I-CeuI linkage maps of the S. dysenteriae type 7 and type 2 isolates were similar to that of S. dysenteriae type 1 and S. flexneri type 2a genomes, respectively.

INTERPRETATION & CONCLUSIONS

Our findings indicate that the type 7 and type 1 isolates of S. dysenteriae were probably evolved from a same precursor, while the type 2 and S. flexneri type 2a were probably evolved and diversified from a common progenitor.

Challenges in enzymatic route of mannitol production

Mannitol is an important biochemical often used as medicine and in food sector, yet its biotechnological is not preffered in Industry for large scale production, which may be due to the multistep mechanism involved in hydrogenation and reduction. This paper is a comparative preview covering present chemical and biotechnological approaches existing today for mannitol production at industrial scale. Biotechnological routes are suitable for adaptation at industrial level for mannitol production, and whatever concerns are there had been discussed in detail, namely, raw materials, broad range of enzymes with high activity at elevated temperature suitable for use in reactor, cofactor limitation, reduced by-product formation, end product inhibition, and reduced utilization of mannitol for enhancing the yield with maximum volumetric productivity.

Myosin-X facilitates Shigella-induced membrane protrusions and cell-to-cell spread

The intracellular pathogen Shigella flexneri forms membrane protrusions to spread from cell to cell. As protrusions form, myosin-X (Myo10) localizes to Shigella. Electron micrographs of immunogold-labelled Shigella-infected HeLa cells reveal that Myo10 concentrates at the bases and along the sides of bacteria within membrane protrusions. Time-lapse video microscopy shows that a full-length Myo10 GFP-construct cycles along the sides of Shigella within the membrane protrusions as these structures progressively lengthen. RNAi knock-down of Myo10 is associated with shorter protrusions with thicker stalks, and causes a >80% decrease in confluent cell plaque formation. Myo10 also concentrates in membrane protrusions formed by another intracellular bacteria, Listeria, and knock-down of Myo10 also impairs Listeria plaque formation. In Cos7 cells (contain low concentrations of Myo10), the expression of full-length Myo10 nearly doubles Shigella-induced protrusion length, and lengthening requires the head domain, as well as the tail-PH domain, but not the FERM domain. The GFP-Myo10-HMM domain localizes to the sides of Shigella within membrane protrusions and the GFP-Myo10-PH domain localizes to host cell membranes. We conclude thatMyo10 generates the force to enhance bacterial-induced protrusions by binding its head region to actin filaments and its PH tail domain to the peripheral membrane.

Shigella flexneri T3SS effector IpaH4.5 modulates the host inflammatory response via interaction with NF-κB p65 protein

Shigella species possess a type III secretion system (T3SS), which is required for human infection and that delivers effector proteins into target host cells. Here, we show that the effector, IpaH4.5 dampens the pro-inflammatory cytokine response. In both the Sereny test and a murine lung infection model, the Shigella ΔipaH4.5 mutant strain caused more severe inflammatory responses and significantly induced higher pro-inflammatory cytokine levels (MIP-2 and TNF-α) in the lung homogenates of wild type-infected mice. Moreover, there was a threefold decrease in bacterial colonization of the mutant compared with the WT and ΔipaH4.5/ipaH4.5-rescued strains. Yeast two-hybrid screening showed that IpaH4.5 specifically interacts with the p65 subunit of NF-κB. Ten truncated versions of IpaH4.5 and p65 spanning different regions were constructed and expressed to further map the IpaH binding sites with p65. The results revealed thatthe p65 region spanning amino acids 1-190 of p65 interacted with the IpaH4.5/1-293 N-terminal region. In vitro, IpaH4.5 displayed ubiquitin ligase activity towards ubiquitin and p65. Furthermore, the transcriptional activity of NF-κB was shown to be inhibited by IpaH4.5 utilizing a dual-luciferase reporter gene detection system containing NF-κB promoter response elements. Thus, we conclude that the IpaH4.5 protein is an E3 ubiquitin ligase capable of directly regulating the host inflammatory response by inhibiting the NF-κB signalling pathway.

Antivirulence genes: insights into pathogen evolution through gene loss

The emergence of new pathogens and the exploitation of novel pathogenic niches by bacteria typically require the horizontal transfer of virulence factors and subsequent adaptation--a "fine-tuning" process--for the successful incorporation of these factors into the microbe's genome. The function of newly acquired virulence factors may be hindered by the expression of genes already present in the bacterium. Occasionally, certain genes must be inactivated or deleted for full expression of the pathogen phenotype to occur. These genes are known as antivirulence genes (AVGs). Originally identified in Shigella, AVGs have improved our understanding of pathogen evolution and provided a novel approach to drug and vaccine development. In this review, we revisit the AVG definition and update the list of known AVGs, which now includes genes from pathogens such as Salmonella, Yersinia pestis, and the virulent Francisella tularensis subspecies. AVGs encompass a wide variety of different roles within the microbe, including genes involved in metabolism, biofilm synthesis, lipopolysaccharide modification, and host vasoconstriction. More recently, the use of one of these AVGs (lpxL) as a potential vaccine candidate highlights the practical application of studying AVG inactivation in microbial pathogens.

Genomic diversity of Escherichia isolates from diverse habitats

Our understanding of the Escherichia genus is heavily biased toward pathogenic or commensal isolates from human or animal hosts. Recent studies have recovered Escherichia isolates that persist, and even grow, outside these hosts. Although the environmental isolates are typically phylogenetically distinct, they are highly related to and phenotypically indistinguishable from their human counterparts, including for the coliform test. To gain insights into the genomic diversity of Escherichia isolates from diverse habitats, including freshwater, soil, animal, and human sources, we carried out comparative DNA-DNA hybridizations using a multi-genome E. coli DNA microarray. The microarray was validated based on hybridizations with selected strains whose genome sequences were available and used to assess the frequency of microarray false positive and negative signals. Our results showed that human fecal isolates share two sets of genes (n>90) that are rarely found among environmental isolates, including genes presumably important for evading host immune mechanisms (e.g., a multi-drug transporter for acids and antimicrobials) and adhering to epithelial cells (e.g., hemolysin E and fimbrial-like adhesin protein). These results imply that environmental isolates are characterized by decreased ability to colonize host cells relative to human isolates. Our study also provides gene markers that can distinguish human isolates from those of warm-blooded animal and environmental origins, and thus can be used to more reliably assess fecal contamination in natural ecosystems.

Shigella flexneri effectors OspE1 and OspE2 mediate induced adherence to the colonic epithelium following bile salts exposure

Shigella flexneri is a Gram-negative pathogen that invades the colonic epithelium. While invasion has been thoroughly investigated, it is unknown how Shigella first attaches to the epithelium. Previous literature suggests that Shigella utilizes adhesins that are induced by environmental signals, including bile salts, encountered in the small intestine prior to invasion. We hypothesized that bile would induce adherence factors to facilitate attachment to colonic epithelial cells. To test our hypothesis, S. flexneri strain 2457T was subcultured in media containing bile salts, and the ability of the bacteria to adhere to the apical surface of polarized T84 epithelial cells was measured. We observed a significant increase in adherence, which was absent in a virulence plasmid-cured strain and a type-III secretion system mutant. Microarray expression analysis indicated that the ospE1/ospE2 genes were induced in the presence of bile, and bile-induced adherence was lost in a ΔospE1/ΔospE2 mutant. Further studies demonstrated that the OspE1/OspE2 proteins were localized to the bacterial outer membrane following exposure to bile salts. The data presented are the first demonstration that the OspE1/OspE2 proteins promote initial adherence to the intestinal epithelium. The adhesins required for Shigella attachment to the colonic epithelium may serve as ideal targets for vaccine development.

Temperature Depended Role of Shigella flexneri Invasion Plasmid on the Interaction with Acanthamoeba castellanii

Shigella flexneri is a Gram-negative bacterium causing the diarrhoeal disease shigellosis in humans. The virulence genes required for invasion are clustered on a large 220 kb plasmid encoding type three secretion system (TTSS) apparatus and virulence factors such as adhesions and invasion plasmid antigens (Ipa). The bacterium is transmitted by contaminated food, water, or from person to person. Acanthamoebae are free-living amoebae (FLA) which are found in diverse environments and isolated from various water sources. Different bacteria interact differently with FLA since Francisella tularensis, Vibrio cholerae, Shigella sonnei, and S. dysenteriae are able to grow inside A. castellanii. In contrast, Pseudomonas aeruginosa induces both necrosis and apoptosis to kill A. castellanii. The aim of this study is to examine the role of invasion plasmid of S. flexneri on the interaction with A. castellanii at two different temperatures. A. castellanii in the absence or presence of wild type, IpaB mutant, or plasmid-cured strain S. flexneri was cultured at 30°C and 37°C and the interaction was analysed by viable count of both bacteria and amoebae, electron microscopy, flow cytometry, and statistical analysis. The outcome of the interaction was depended on the temperature since the growth of A. castellanii was inhibited at 30°C, and A. castellanii was killed by invasion plasmid mediated necrosis at 37°C.

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

Background

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

Methodologies/Principal Findings

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

Conclusions

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

BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons

BACKGROUND

Visualisation of genome comparisons is invaluable for helping to determine genotypic differences between closely related prokaryotes. New visualisation and abstraction methods are required in order to improve the validation, interpretation and communication of genome sequence information; especially with the increasing amount of data arising from next-generation sequencing projects. Visualising a prokaryote genome as a circular image has become a powerful means of displaying informative comparisons of one genome to a number of others. Several programs, imaging libraries and internet resources already exist for this purpose, however, most are either limited in the number of comparisons they can show, are unable to adequately utilise draft genome sequence data, or require a knowledge of command-line scripting for implementation. Currently, there is no freely available desktop application that enables users to rapidly visualise comparisons between hundreds of draft or complete genomes in a single image.

RESULTS

BLAST Ring Image Generator (BRIG) can generate images that show multiple prokaryote genome comparisons, without an arbitrary limit on the number of genomes compared. The output image shows similarity between a central reference sequence and other sequences as a set of concentric rings, where BLAST matches are coloured on a sliding scale indicating a defined percentage identity. Images can also include draft genome assembly information to show read coverage, assembly breakpoints and collapsed repeats. In addition, BRIG supports the mapping of unassembled sequencing reads against one or more central reference sequences. Many types of custom data and annotations can be shown using BRIG, making it a versatile approach for visualising a range of genomic comparison data. BRIG is readily accessible to any user, as it assumes no specialist computational knowledge and will perform all required file parsing and BLAST comparisons automatically.

CONCLUSIONS

There is a clear need for a user-friendly program that can produce genome comparisons for a large number of prokaryote genomes with an emphasis on rapidly utilising unfinished or unassembled genome data. Here we present BRIG, a cross-platform application that enables the interactive generation of comparative genomic images via a simple graphical-user interface. BRIG is freely available for all operating systems at http://sourceforge.net/projects/brig/.

The Vibrio cholerae VctPDGC system transports catechol siderophores and a siderophore-free iron ligand

Vibrio cholerae, the causative agent of cholera, has an absolute requirement for iron. It transports the catechol siderophores vibriobactin, which it synthesizes and secretes, and enterobactin. These siderophores are transported across the inner membrane by one of two periplasmic binding protein-dependent ABC transporters, VctPDGC or ViuPDGC. We show here that one of these inner membrane transport systems, VctPDGC, also promotes iron acquisition in the absence of siderophores. Plasmids carrying the vctPDGC genes stimulated growth in both rich and minimal media of a Shigella flexneri mutant that produces no siderophores. vctPDGC also stimulated the growth of an Escherichia coli enterobactin biosynthetic mutant in low iron medium, and this effect did not require feoB, tonB or aroB. A tyrosine to phenylalanine substitution in the periplasmic binding protein VctP did not alter enterobactin transport, but eliminated growth stimulation in the absence of a siderophore. These data suggest that the VctPDGC system has the capacity to transport both catechol siderophores and a siderophore-free iron ligand. We also show that VctPDGC is the previously unidentified siderophore-independent iron transporter in V. cholerae, and this appears to complete the list of iron transport systems in V. cholerae.

Impact of small repeat sequences on bacterial genome evolution

Intergenic regions of prokaryotic genomes carry multiple copies of terminal inverted repeat (TIR) sequences, the nonautonomous miniature inverted-repeat transposable element (MITE). In addition, there are the repetitive extragenic palindromic (REP) sequences that fold into a small stem loop rich in G–C bonding. And the clustered regularly interspaced short palindromic repeats (CRISPRs) display similar small stem loops but are an integral part of a complex genetic element. Other classes of repeats such as the REP2 element do not have TIRs but show other signatures. With the current availability of a large number of whole-genome sequences, many new repeat elements have been discovered. These sequences display diverse properties. Some show an intimate linkage to integrons, and at least one encodes a small RNA. Many repeats are found fused with chromosomal open reading frames, and some are located within protein coding sequences. Small repeat units appear to work hand in hand with the transcriptional and/or post-transcriptional apparatus of the cell. Functionally, they are multifaceted, and this can range from the control of gene expression, the facilitation of host/pathogen interactions, or stimulation of the mammalian immune system. The CRISPR complex displays dramatic functions such as an acquired immune system that defends against invading viruses and plasmids. Evolutionarily, mobile repeat elements may have influenced a cycle of active versus inactive genes in ancestral organisms, and some repeats are concentrated in regions of the chromosome where there is significant genomic plasticity. Changes in the abundance of genomic repeats during the evolution of an organism may have resulted in a benefit to the cell or posed a disadvantage, and some present day species may reflect a purification process. The diverse structure, eclectic functions, and evolutionary aspects of repeat elements are described.

Emerging vaccine informatics

Vaccine informatics is an emerging research area that focuses on development and applications of bioinformatics methods that can be used to facilitate every aspect of the preclinical, clinical, and postlicensure vaccine enterprises. Many immunoinformatics algorithms and resources have been developed to predict T- and B-cell immune epitopes for epitope vaccine development and protective immunity analysis. Vaccine protein candidates are predictable in silico from genome sequences using reverse vaccinology. Systematic transcriptomics and proteomics gene expression analyses facilitate rational vaccine design and identification of gene responses that are correlates of protection in vivo. Mathematical simulations have been used to model host-pathogen interactions and improve vaccine production and vaccination protocols. Computational methods have also been used for development of immunization registries or immunization information systems, assessment of vaccine safety and efficacy, and immunization modeling. Computational literature mining and databases effectively process, mine, and store large amounts of vaccine literature and data. Vaccine Ontology (VO) has been initiated to integrate various vaccine data and support automated reasoning.

Chromosomal rearrangements in Salmonella enterica serovar Typhi strains isolated from asymptomatic human carriers

Host-specific serovars of Salmonella enterica often have large-scale chromosomal rearrangements that occur by recombination between rrn operons. Two hypotheses have been proposed to explain these rearrangements: (i) replichore imbalance from horizontal gene transfer drives the rearrangements to restore balance, or (ii) the rearrangements are a consequence of the host-specific lifestyle. Although recent evidence has refuted the replichore balance hypothesis, there has been no direct evidence for the lifestyle hypothesis. To test this hypothesis, we determined the rrn arrangement type for 20 Salmonella enterica serovar Typhi strains obtained from human carriers at periodic intervals over multiple years. These strains were also phage typed and analyzed for rearrangements that occurred over long-term storage versus routine culturing. Strains isolated from the same carrier at different time points often exhibited different arrangement types. Furthermore, colonies isolated directly from the Dorset egg slants used to store the strains also had different arrangement types. In contrast, colonies that were repeatedly cultured always had the same arrangement type. Estimated replichore balance of isolated strains did not improve over time, and some of the rearrangements resulted in decreased replicore balance. Our results support the hypothesis that the restricted lifestyle of host-specific Salmonella is responsible for the frequent chromosomal rearrangements in these serovars.

Hypervirulent Chlamydia trachomatis clinical strain is a recombinant between lymphogranuloma venereum (L(2)) and D lineages

Chlamydia trachomatis is an obligate intracellular bacterium that causes a diversity of severe and debilitating diseases worldwide. Sporadic and ongoing outbreaks of lymphogranuloma venereum (LGV) strains among men who have sex with men (MSM) support the need for research on virulence factors associated with these organisms. Previous analyses have been limited to single genes or genomes of laboratory-adapted reference strain L₂/434 and outbreak strain L₂b/UCH-1/proctitis. We characterized an unusual LGV strain, termed L₂c, isolated from an MSM with severe hemorrhagic proctitis. L₂c developed nonfusing, grape-like inclusions and a cytotoxic phenotype in culture, unlike the LGV strains described to date. Deep genome sequencing revealed that L₂c was a recombinant of L₂ and D strains with conserved clustered regions of genetic exchange, including a 78-kb region and a partial, yet functional, toxin gene that was lost with prolonged culture. Indels (insertions/deletions) were discovered in an ftsK gene promoter and in the tarp and hctB genes, which encode key proteins involved in replication, inclusion formation, and histone H1-like protein activity, respectively. Analyses suggest that these indels affect gene and/or protein function, supporting the in vitro and disease phenotypes. While recombination has been known to occur for C. trachomatis based on gene sequence analyses, we provide the first whole-genome evidence for recombination between a virulent, invasive LGV strain and a noninvasive common urogenital strain. Given the lack of a genetic system for producing stable C. trachomatis mutants, identifying naturally occurring recombinants can clarify gene function and provide opportunities for discovering avenues for genomic manipulation.

Lymphogranuloma venereum (LGV) is a prevalent and debilitating sexually transmitted disease in developing countries, although there are significant ongoing outbreaks in Australia, Europe, and the United States among men who have sex with men (MSM). Relatively little is known about LGV virulence factors, and only two LGV genomes have been sequenced to date. We isolated an LGV strain from an MSM with severe hemorrhagic proctitis that was morphologically unique in tissue culture compared with other LGV strains. Bioinformatic and statistical analyses identified the strain as a recombinant of L₂ and D strains with highly conserved clustered regions of genetic exchange. The unique culture morphology and, more importantly, disease phenotype could be traced to the genes involved in recombination. The findings have implications for bacterial species evolution and, in the case of ongoing LGV outbreaks, suggest that recombination is a mechanism for strain emergence that results in significant disease pathology.

An integrated approach for finding overlooked genes in Shigella

BACKGROUND

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

METHODOLOGY/PRINCIPAL FINDINGS

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

CONCLUSIONS/SIGNIFICANCE

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

IS3 profiling identifies the enterohaemorrhagic Escherichia coli O-island 62 in a distinct enteroaggregative E. coli lineage

Background

Enteroaggregative Escherichia coli (EAEC) are important diarrhoeal pathogens that are defined by a HEp-2 adherence assay performed in specialist laboratories. Multilocus sequence typing (MLST) has revealed that aggregative adherence is convergent, providing an explanation for why not all EAEC hybridize with the plasmid-derived probe for this category, designated CVD432. Some EAEC lineages are globally disseminated or more closely associated with disease.

Results

To identify genetic loci conserved within significant EAEC lineages, but absent from non-EAEC, IS3-based PCR profiles were generated for 22 well-characterised EAEC strains. Six bands that were conserved among, or missing from, specific EAEC lineages were cloned and sequenced. One band corresponded to the aggR gene, a plasmid-encoded regulator that has been used as a diagnostic target but predominantly detects EAEC bearing the plasmid already marked by CVD432. The sequence from a second band was homologous to an open-reading frame within the cryptic enterohaemorrhagic E. coli (EHEC) O157 genomic island, designated O-island 62. Screening of an additional 46 EAEC strains revealed that the EHEC O-island 62 was only present in those EAEC strains belonging to the ECOR phylogenetic group D, largely comprised of sequence type (ST) complexes 31, 38 and 394.

Conclusions

The EAEC 042 gene orf1600, which lies within the EAEC equivalent of O-island 62 island, can be used as a marker for EAEC strains belonging to the ECOR phylogenetic group D. The discovery of EHEC O-island 62 in EAEC validates the genetic profiling approach for identifying conserved loci among phylogenetically related strains.

Distribution of genes encoding nucleoid-associated protein homologs in plasmids

Bacterial nucleoid-associated proteins (NAPs) form nucleoprotein complexes and influence the expression of genes. Recent studies have shown that some plasmids carry genes encoding NAP homologs, which play important roles in transcriptional regulation networks between plasmids and host chromosomes. In this study, we determined the distributions of the well-known NAPs Fis, H-NS, HU, IHF, and Lrp and the newly found NAPs MvaT and NdpA among the whole-sequenced 1382 plasmids found in Gram-negative bacteria. Comparisons between NAP distributions and plasmid features (size, G+C content, and putative transferability) were also performed. We found that larger plasmids frequently have NAP gene homologs. Plasmids with H-NS gene homologs had less G+C content. It should be noted that plasmids with the NAP gene homolog also carried the relaxase gene involved in the conjugative transfer of plasmids more frequently than did those without the NAP gene homolog, implying that plasmid-encoded NAP homologs positively contribute to transmissible plasmids.

Genomic characterization of asymptomatic Escherichia coli isolated from the neobladder

The replacement of the bladder with a neobladder made from ileal tissue is the prescribed treatment in some cases of bladder cancer or trauma. Studies have demonstrated that individuals with an ileal neobladder have recurrent colonization by Escherichia coli and other species that are commonly associated with urinary tract infections; however, pyelonephritis and complicated symptomatic infections with ileal neobladders are relatively rare. This study examines the genomic content of two E. coli isolates from individuals with neobladders using comparative genomic hybridization (CGH) with a pan-E. coli/Shigella microarray. Comparisons of the neobladder genome hybridization patterns with reference genomes demonstrate that the neobladder isolates are more similar to the commensal, laboratory-adapted E. coli and a subset of enteroaggregative E. coli than they are to uropathogenic E. coli isolates. Genes identified by CGH as exclusively present in the neobladder isolates among the 30 examined isolates were primarily from large enteric isolate plasmids. Isolations identified a large plasmid in each isolate, and sequencing confirmed similarity to previously identified plasmids of enteric species. Screening, via PCR, of more than 100 isolates of E. coli from environmental, diarrhoeagenic and urinary tract sources did not identify neobladder-specific genes that were widely distributed in these populations. These results taken together demonstrate that the neobladder isolates, while distinct, are genomically more similar to gastrointestinal or commensal E. coli, suggesting why they can colonize the transplanted intestinal tissue but rarely progress to acute pyelonephritis or more severe disease.

The genome sequence of E. coli W (ATCC 9637): comparative genome analysis and an improved genome-scale reconstruction of E. coli

BACKGROUND

Escherichia coli is a model prokaryote, an important pathogen, and a key organism for industrial biotechnology. E. coli W (ATCC 9637), one of four strains designated as safe for laboratory purposes, has not been sequenced. E. coli W is a fast-growing strain and is the only safe strain that can utilize sucrose as a carbon source. Lifecycle analysis has demonstrated that sucrose from sugarcane is a preferred carbon source for industrial bioprocesses.

RESULTS

We have sequenced and annotated the genome of E. coli W. The chromosome is 4,900,968 bp and encodes 4,764 ORFs. Two plasmids, pRK1 (102,536 bp) and pRK2 (5,360 bp), are also present. W has unique features relative to other sequenced laboratory strains (K-12, B and Crooks): it has a larger genome and belongs to phylogroup B1 rather than A. W also grows on a much broader range of carbon sources than does K-12. A genome-scale reconstruction was developed and validated in order to interrogate metabolic properties.

CONCLUSIONS

The genome of W is more similar to commensal and pathogenic B1 strains than phylogroup A strains, and therefore has greater utility for comparative analyses with these strains. W should therefore be the strain of choice, or 'type strain' for group B1 comparative analyses. The genome annotation and tools created here are expected to allow further utilization and development of E. coli W as an industrial organism for sucrose-based bioprocesses. Refinements in our E. coli metabolic reconstruction allow it to more accurately define E. coli metabolism relative to previous models.

Symbionts and pathogens: what is the difference?

The ecological relationships that organisms establish with others can be considered as broad and diverse as the forms of life that inhabit and interact in our planet. Those interactions can be considered as a continuum spectrum, ranging from beneficial to detrimental outcomes. However, this picture has revealed as more complex and dynamic than previously thought, involving not only factors that affect the two or more members that interact, but also external forces, with chance playing a crucial role in this interplay. Thus, defining a particular symbiont as mutualist or pathogen in an exclusive way, based on simple rules of classification is increasingly challenging if not unfeasible, since new methodologies are providing more evidences that depict exceptions, reversions and transitions within either side of this continuum, especially evident at early stages of symbiotic associations. This imposes a wider and more dynamic view of a complex landscape of interactions.

Comparison of 61 sequenced Escherichia coli genomes

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

High-throughput microarray technology in diagnostics of enterobacteria based on genome-wide probe selection and regression analysis

Background

The Enterobacteriaceae comprise a large number of clinically relevant species with several individual subspecies. Overlapping virulence-associated gene pools and the high overall genome plasticity often interferes with correct enterobacterial strain typing and risk assessment. Array technology offers a fast, reproducible and standardisable means for bacterial typing and thus provides many advantages for bacterial diagnostics, risk assessment and surveillance. The development of highly discriminative broad-range microbial diagnostic microarrays remains a challenge, because of marked genome plasticity of many bacterial pathogens.

Results

We developed a DNA microarray for strain typing and detection of major antimicrobial resistance genes of clinically relevant enterobacteria. For this purpose, we applied a global genome-wide probe selection strategy on 32 available complete enterobacterial genomes combined with a regression model for pathogen classification. The discriminative power of the probe set was further tested in silico on 15 additional complete enterobacterial genome sequences. DNA microarrays based on the selected probes were used to type 92 clinical enterobacterial isolates. Phenotypic tests confirmed the array-based typing results and corroborate that the selected probes allowed correct typing and prediction of major antibiotic resistances of clinically relevant Enterobacteriaceae, including the subspecies level, e.g. the reliable distinction of different E. coli pathotypes.

Conclusions

Our results demonstrate that the global probe selection approach based on longest common factor statistics as well as the design of a DNA microarray with a restricted set of discriminative probes enables robust discrimination of different enterobacterial variants and represents a proof of concept that can be adopted for diagnostics of a wide range of microbial pathogens. Our approach circumvents misclassifications arising from the application of virulence markers, which are highly affected by horizontal gene transfer. Moreover, a broad range of pathogens have been covered by an efficient probe set size enabling the design of high-throughput diagnostics.