Skip to main content
Advertisement
  • Loading metrics

The Serum Resistome of a Globally Disseminated Multidrug Resistant Uropathogenic Escherichia coli Clone

  • Minh-Duy Phan,

    Affiliations Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

  • Kate M. Peters,

    Affiliations Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

  • Sohinee Sarkar,

    Affiliations Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

  • Samuel W. Lukowski,

    Current address: Department of Genetic Medicine and Development, University of Geneva, Geneva, Switzerland.

    Affiliation Queensland Brain Institute, The University of Queensland, Brisbane, Queensland, Australia

  • Luke P. Allsopp,

    Affiliations Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

  • Danilo Gomes Moriel,

    Affiliations Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

  • Maud E. S. Achard,

    Affiliations Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

  • Makrina Totsika,

    Affiliations Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

  • Vikki M. Marshall,

    Affiliation Queensland Brain Institute, The University of Queensland, Brisbane, Queensland, Australia

  • Mathew Upton,

    Affiliation Faculty of Medicine and Dentistry, University of Plymouth, Plymouth, United Kingdom

  • Scott A. Beatson,

    Affiliations Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

  • Mark A. Schembri

    m.schembri@uq.edu.au

    Affiliations Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

Abstract

Escherichia coli ST131 is a globally disseminated, multidrug resistant clone responsible for a high proportion of urinary tract and bloodstream infections. The rapid emergence and successful spread of E. coli ST131 is strongly associated with antibiotic resistance; however, this phenotype alone is unlikely to explain its dominance amongst multidrug resistant uropathogens circulating worldwide in hospitals and the community. Thus, a greater understanding of the molecular mechanisms that underpin the fitness of E. coli ST131 is required. In this study, we employed hyper-saturated transposon mutagenesis in combination with multiplexed transposon directed insertion-site sequencing to define the essential genes required for in vitro growth and the serum resistome (i.e. genes required for resistance to human serum) of E. coli EC958, a representative of the predominant E. coli ST131 clonal lineage. We identified 315 essential genes in E. coli EC958, 231 (73%) of which were also essential in E. coli K-12. The serum resistome comprised 56 genes, the majority of which encode membrane proteins or factors involved in lipopolysaccharide (LPS) biosynthesis. Targeted mutagenesis confirmed a role in serum resistance for 46 (82%) of these genes. The murein lipoprotein Lpp, along with two lipid A-core biosynthesis enzymes WaaP and WaaG, were most strongly associated with serum resistance. While LPS was the main resistance mechanism defined for E. coli EC958 in serum, the enterobacterial common antigen and colanic acid also impacted on this phenotype. Our analysis also identified a novel function for two genes, hyxA and hyxR, as minor regulators of O-antigen chain length. This study offers novel insight into the genetic make-up of E. coli ST131, and provides a framework for future research on E. coli and other Gram-negative pathogens to define their essential gene repertoire and to dissect the molecular mechanisms that enable them to survive in the bloodstream and cause disease.

Author Summary

The emergence and rapid dissemination of new bacterial pathogens presents multiple challenges to healthcare systems, including the need for rapid detection, precise diagnostics, effective transmission control and effective treatment. E. coli ST131 is an example of a recently emerged multidrug resistant pathogen that is capable of causing urinary tract and bloodstream infections with limited available treatment options. In order to increase our molecular understanding of E. coli ST131, we developed a high-throughput transposon mutagenesis system in combination with next generation sequencing to test every gene for its essential role in growth and for its contribution to serum resistance. We identified 315 essential genes, 270 of which were conserved among all currently available complete E. coli genomes. Fifty-six genes that define the serum resistome of E. coli ST131 were identified, including genes encoding membrane proteins, proteins involved in LPS biosynthesis, regulators and several novel proteins with previously unknown function. This study therefore provides an inventory of essential and serum resistance genes that could form a framework for the future development of targeted therapeutics to prevent disease caused by multidrug-resistant E. coli ST131 strains.

Introduction

Escherichia coli O25b:H4-ST131 (E. coli ST131) is a recently emerged, globally disseminated clone that is often multidrug resistant and is responsible for a high proportion of community- and nosocomially-acquired urinary tract and bloodstream infections [1][6]. E. coli ST131 strains are also capable of causing complicated infections including acute pyelonephritis, osteomyelitis, septic arthritis and septic shock [7], [8]. E. coli ST131 are commonly associated with production of the CTX-M-15 enzyme, currently the most widespread extended spectrum β-lactamase (ESBL) of its type in the world [1], [9]. In addition to resistance against oxyimino-cephalosporins (i.e. cefotaxime, ceftazidime), and monobactams, E. coli ST131 strains are often co-resistant to fluoroquinolones [3], [10]. Indeed, most fluoroquinolone-resistant E. coli strains belong to a recently emerged and dominant subgroup of ST131 strains [11]. Some E. coli ST131 strains have also been reported to produce carbapenemases [12][14], thus severely limiting treatment options that are currently available against this clinically predominant clone [15].

E. coli ST131 strains, like many other uropathogenic E. coli (UPEC) strains, are derived from phylogenetic group B2 [3]. Typically, UPEC strains possess a large and diverse range of virulence factors that contribute to their ability to cause urinary tract and bloodstream infections, including adhesins, toxins, siderophores and protectins [15], [16]. Several studies have demonstrated that E. coli ST131 strains possess a similar suite of virulence factors and cause invasive disease, leading to the hypothesis that the widespread pathogenic success of E. coli ST131 strains may be in part due to enhanced virulence [7], [8], [17]. However, it has become clear from recent studies that E. coli ST131 strains do not possess a heightened virulence potential compared to other UPEC or B2 E. coli strains in causing invasive infections [18] or infections in nematodes and zebrafish embryos [19]. Thus, other factors such as enhanced metabolic capacity have been proposed to contribute to the fitness and pathogenic success of this dominant clone [20], [21].

The genome sequence of one of the best-characterized E. coli ST131 strains, EC958, was recently described [22]. E. coli EC958 is a member of the pulsed-field gel electrophoresis (PFGE) defined UK epidemic strain A, which represents one of the major pathogenic lineages (PFGE strains A–E) of ESBL producing E. coli causing urinary tract infections (UTI) across the UK [23]. E. coli EC958 is resistant to eight antibiotic classes, including oxyimino-cephalosporins, fluoroquinolones and sulphonamides. E. coli EC958 colonizes the bladder of mice in a type 1 fimbriae-dependent manner [22], can invade into bladder epithelial cells and form intracellular bacterial communities, can establish both acute and chronic UTI [24] and can inhibit the contraction of ureters, in vitro [25].

The ability to resist the bactericidal activity of serum, and thus survive in the bloodstream, represents an essential virulence trait for UPEC and other extra-intestinal E. coli (ExPEC) strains, including E. coli ST131 [26][28]. In E. coli, several mechanisms have been shown to contribute to serum resistance. The importance of O-antigens and K capsules in resistance to serum has been recognized since the 1960s and 1980s, respectively [29][31]; and their multiple types, combinations and length contribute differently to serum resistance [32][35]. The major outer membrane protein OmpA [36], plasmid-encoded proteins TraT [37], [38] and Iss [39], and the phage membrane protein Bor [40] have also been reported to contribute to serum resistance in E. coli. Notably, each of these resistance mechanisms has been studied in isolation and in different strain backgrounds. Thus, while serum resistance is clearly a complex phenotype determined by multiple elements, little is known about the combination of factors that contribute to resistance in a single strain.

High-throughput transposon mutagenesis combined with genome-wide targeted sequencing was used recently to study the essential genes in Salmonella enterica serovar Typhi and Caulobacter crescentus [41], [42]. Langridge et al. also used their transposon directed insertion-site sequencing (TraDIS) method to assay every gene for its role in the survival of S. Typhi in the presence of bile salts [42]. Similar approaches (INSeq, HITS, Tn-seq) have also been applied to a range of organisms to study gene requirements for survival in particular niches [43][46]. Here, we adapted TraDIS and designed a multiplexing method to define the essential genes required for in vitro growth (i.e. Luria-Bertani agar media supplemented with 30 µg/ml Cm at 37°C) and the serum resistome (i.e. genes required for resistance to human serum, in E. coli EC958). We show that the essential gene list of E. coli EC958 comprises 315 genes, 231 of which are shared with E. coli K-12. We also define for the first time a comprehensive inventory of genes required for resistance to human serum. Our study provides a molecular blueprint for understanding the mechanisms employed by E. coli ST131 to survive, grow in the bloodstream and cause disease.

Results

Application of multiplexed TraDIS to the E. coli ST131 strain EC958

Approximately 1 million mutants were generated in the E. coli ST131 strain EC958 [22] using an in-house miniTn5 transposon carrying a chloramphenicol (Cm) resistance gene derived from the pKD3 plasmid [47]. A primer comprising four functional regions was designed to facilitate specific sequencing of the transposon insertion sites on the Illumina HiSeq 2000 platform while allowing for intra-lane multiple sample indexing (Figure S1). This primer contained (5′-3′): (i) the P5 sequence to bind to TruSeq flowcells, (ii) the Illumina read 1 sequencing primer binding site, (iii) a 6-bp index sequence for multiple sample barcoding and (iv) a 25-bp transposon specific sequence designed to amplify the last 12 bp of the transposon and its adjacent genomic sequence. Using this custom primer, we successfully sequenced the transposon insertion sites for 6 samples on both TruSeq version 2 and version 3 flowcells (Figure 1A). Each sample yielded from 6.8 million to 15 million reads that were tagged with transposon specific sequence, 71% of which were reliably mapped to EC958 draft chromosome (excluding unscaffolded contigs and plasmids) (Table 1). All experiments were performed in duplicate, with the correlation coefficient for the number of insertions per gene for each pair of samples close to 1 (R2>0.99) and thus demonstrating a high level of reproducibility for each experiment (Figure 1B).

thumbnail
Figure 1. Summary of sequence read data from multiplexed TraDIS.

(A) Number of tagged reads obtained in each lane of TruSeq version 2 and version 3 flowcells. (B) Correlation between the two biological replicates. The number of insertions in each gene from Test A was plotted against that from Test B; the correlation coefficient R2 indicates very low variation between the two replicates.

https://doi.org/10.1371/journal.pgen.1003834.g001

thumbnail
Table 1. Summary of sequencing and mapping results of multiplexed TraDIS runs.

https://doi.org/10.1371/journal.pgen.1003834.t001

Essential genes in E. coli EC958

We initially used our saturated random insertion mutant library to determine the ‘essential genes of EC958’, defined as the set of genes required for growth on LB agar supplemented with Cm 30 µg/ml. We extracted genomic DNA (in duplicates: input A and B) directly from the library pool and sequenced using our multiplexed TraDIS protocol. We combined the reads from input A and B to maximize the coverage resulting in 16 million transposon-tagged reads, of which 11 million uniquely mapped to the EC958 chromosome, resulting in 502,068 unique insertion sites. This equates to an average of one insertion site every 9.92 bp, with a very low probability of having 100 consecutive bp without interruption by chance (P = 4.2×10−5).

The essential gene list was identified using a statistical analysis similar to that described by Langridge et al., which recognized two distinct distributions of insertion indexes (number of insertions divided by gene length) for non-essential genes (gamma) and essential genes (exponential) and called those with insertion indexes less than or equal to the intercept of the two distributions as essential [42]. In our data, an insertion index cut-off of 0.0158, resulted in the identification of 315 genes as essential (Table S1). This cut-off is equivalent to a log2-likelihood-ratio (LLR) of −3.6, which means that our essential genes are at least 12 times more likely to belong to the exponential distribution (essential) than the gamma (non-essential) distribution.

The functional category of each gene was identified based on the COG (Clusters of Orthologous Groups) numbers from the EC958 annotation (accession number PRJEA61443). Figure 2 shows an overview of essential functions in EC958 compared with the total number of genes in each functional category. Genes involved in translation, ribosomal structure and biogenesis account for 25% of the total number of essential genes in EC958, which is 42% of the total number of genes in this category. The second most abundant category in the essential gene list comprised genes involved in cell wall/membrane/envelope biogenesis (12%), followed by genes involved in coenzyme transport and metabolism. There were 23 essential genes with functions not identified in the COG database.

thumbnail
Figure 2. Number of essential genes in each COG functional category.

https://doi.org/10.1371/journal.pgen.1003834.g002

To investigate the conservation of EC958 essential proteins among different E. coli pathotypes, we performed tfastx alignment (FASTA v36) between the essential protein sequences and translated DNA from 50 E. coli complete genomes (Table S2). There were 270 (86%) proteins conserved across all genomes investigated. An additional 17 proteins were also present in more than 90% of the genomes. Only 6 proteins were specific for EC958 (not found in 50 genomes) (Table S3).

The serum resistome of E. coli EC958

Saturated transposon mutagenesis in combination with next-generation sequencing is a powerful tool for whole genome, high-throughput identification of all candidate genes involved in a particular phenotype. Here, we used our transposon mutant library in combination with TraDIS to identify genes from EC958 involved in resistance to human serum, thus enabling us to define the serum resistome of EC958. We designed a mutant selection procedure in which 1 million mutants were exposed to pooled fresh human serum for 90 minutes and then allowed to grow in LB broth for 4 hours before genomic DNA extraction. This procedure permitted the growth of serum resistant mutants while eliminating or inhibiting mutants that were sensitive to serum (Figure 3A). The procedure was performed in parallel with control samples where fresh serum was replaced by inactivated serum that lacked bactericidal activity (data not shown). The genomic DNA from test and control samples were sequenced using our modified Illumina multiplexed TraDIS procedure (Figure 3B) to generate multiple datasets (Figure 1A) that were analysed by the Bioconductor package edgeR after filtering out genes identified as essential [48].

thumbnail
Figure 3. Experimental design to identify serum resistance genes in EC958.

(A) Selection steps employed using fresh serum as test and inactivated serum as control. (B) Schematic illustration of the Illumina sequencing procedure including the use of a custom oligo for indexing and enrichment of insert sites.

https://doi.org/10.1371/journal.pgen.1003834.g003

The serum resistance genes were identified as genes that have significant reduction in read counts in the test samples compared to the control samples (i.e. less mutants survived after serum treatment) (Table S4). A stringent threshold of log2 fold change of read counts (logFC) less than −1 and an adjusted p-value less than 0.001 was used to identify significant genes that are involved in serum resistance (Figure S2). Figure 4 shows the names and genomic locations of the 56 genes that satisfied these stringent criteria.

thumbnail
Figure 4. Overview of the E. coli EC958 serum resistance genes.

The circular diagram depicts the location of 56 serum resistance genes on the E. coli EC958 genome. The two outer rings containing blue and red arrows illustrate the CDS and serum resistance genes, respectively, on the forward and reverse strand of the genome. The inner red ring represents the logFC between the control and the test samples for each CDS. The three insets represent a close-up look at three regions on the genome with the graphs showing the location and relative number of each mutant found in the control (blue) and test (red) samples.

https://doi.org/10.1371/journal.pgen.1003834.g004

Twenty-two (39.3%) of the genes belong to three operons responsible for LPS biosynthesis (including both O-antigen biosynthesis and lipid A-core biosynthesis) and enterobacterial common antigen (ECA) biosynthesis. This result represents the first layer of evidence demonstrating the importance of the O25 antigen as well as ECA in an E. coli ST131 background for protection from the bactericidal activity of human serum. Detailed characterisation of the O25 antigen gene cluster is discussed in subsequent sections. ECA is common to all Enterobacteriaceae and is expressed in both serum resistant (smooth) and sensitive (rough) strains, except for rough strains that are defective in the shared biosynthesis pathway affecting both O-antigen and ECA [49], [50]. Seven out of 12 genes in the ECA operon were required for serum resistance as determined by the TraDIS technique.

The remaining 60.7% of genes in the serum resistome identified by TraDIS included genes encoding lipoprotein, membrane proteins, regulators and hypothetical proteins (Table 2). Some of these also affect LPS such as rfaH (EC958_4322), encoding a known regulator required for LPS biosynthesis [51], [52] and virulence of pathogenic E. coli strains [53], [54], whilst others represent genes that have not previously been shown to be associated with serum resistance. The murein lipoprotein gene lpp (EC958_1897) showed the greatest difference between the test and control samples (logFC of −10), followed by pgm (EC958_0806), encoding phosphoglucomutase. Four hypothetical proteins were also identified, two of which (EC958_0460 and EC958_0461) were further characterized in this study (see below).

thumbnail
Table 2. Genes required for serum resistance identified by multiplexed TraDIS.

https://doi.org/10.1371/journal.pgen.1003834.t002

Validation of serum resistance genes

As mentioned above, we employed a stringent threshold combining fold change and statistical significance to define the set of 56 genes in the EC958 serum resistome. In order to validate these findings, we attempted to test all 56 genes independently for their role in serum resistance. Using a modified lambda red mediated-homologous recombination approach [22], [47] we successfully generated defined knock-out mutants for 54 genes (96.4%) in EC958. We were unable to obtain mutants for the remaining 2 genes (acrA and EC958_2373) despite multiple attempts (Table 2).

The 54 defined mutants were subjected to serum susceptibility testing, whereby the number of surviving colonies after a 90-minute exposure to fresh pooled human serum was compared to the number of colonies prior to treatment (Table 2). A mutant was defined as susceptible to serum when its log difference was at least 1 (i.e. 10 fold reduction after exposure to serum). Forty-one genes (75.9%) contributed to serum resistance in EC958 using this assay. In the case of the remaining 13 mutants, it is possible that the lack of susceptibility to human serum was a reflection of the assay, suggesting that survival in serum within a mixed population of one million mutants may be very different from survival of a pure population carrying the same defective mutation. Therefore, to better mimic the condition of TraDIS library serum selection, a competitive assay was devised where EC958 wild-type was mixed equally with a mutant before exposure to serum and the competitive index of the mutant was measured. Using this competitive assay, five of the twelve mutants were significantly attenuated compared to the wild-type EC958 strain (Table 2). Thus, the overall number of validated susceptible mutants was 46 out of 54 tested (85.2%).

Additional characteristics of serum resistant mutants

We hypothesized that one mechanism associated with enhanced sensitivity to human serum could be due to decreased membrane integrity caused by destabilization of the outer leaflet of the outer membrane. In order to test this, we examined the survival of the 54 mutants in response to outer membrane stresses (i.e. SDS) and osmotic potential (i.e. NaCl). In total, 50.0% (27/54) of the mutants displayed enhanced sensitivity to SDS and 11.1% (6/54) of the mutants displayed enhanced sensitivity to NaCl (Table 2). A comparative analysis of these phenotypes in the context of serum sensitivity is presented below.

O25b antigen biosynthesis genes conferring resistance to serum

It is well established that O-antigen represents the main determinant for serum resistance in E. coli. However, there are more than 180 different O-antigens that have been defined in E. coli and these may contribute differently to serum resistance in individual bacterial strains [32], [33]. Furthermore, direct genetic evidence linking O-antigen biosynthesis to serum resistance is only available for a small number of specific O-antigen types. With this in mind, a detailed characterization of the O25b biosynthesis genes was performed using sequence comparison for function prediction in combination with analysis of LPS composition to deduce the role of each gene in resistance to serum.

Similar to most E. coli strains, the O-antigen biosynthesis operon is located between the galF and gnd genes in EC958. Figure 5A shows a comparison of the EC958 O-antigen operon with the equivalent operon from the K-12 strain MG1655 and the O25 serotype E. coli strain E47a [32], [55].

thumbnail
Figure 5. Characterization of O25b antigen genes in EC958.

(A) Comparison of the O-antigen cluster in EC958 with MG1655 (NC_000913) and E47a (GU014554); genes in green are involved in sugar biosynthesis, genes in purple are involved in O-antigen processing and genes in orange encode glycosyl transferase enzymes; genes with a black outline were shown to be required for serum resistance and genes with a red outline were defined as essential. (B) Predicted function of the four glycosyl transferases in the biosynthesis of the O25b repeat unit [57], [58]. (C) LPS gels showing changes in O-antigen structures corresponding to mutations in each gene.

https://doi.org/10.1371/journal.pgen.1003834.g005

Nucleotide sugar biosynthesis genes.

The four genes encoding dTDP-α-L-rhamnose biosynthesis enzymes were highly conserved in the three strains, consistent with the presence of α-L-rhamnose in the O16 and O25 antigen repeat unit expressed by MG1655 and EC958, respectively [56][58]. The dTDP-α-L-rhamnose biosynthesis pathway starts from α-D-glucose-1-phosphate to produce dTDP-α-L-rhamnose via the catalysis of RmlA (to make dTDP-α-D-glucose), RmlB (to make dTDP-4-dehydro-6-deoxy-α-D-glucose), RmlC (to make dTDP-4-dehydro-6-deoxy-β-mannose) and RmlD (to make dTDP-α-L-rhamnose), respectively. LPS gel analysis of EC958 mutants lacking one of these four enzymes showed that RmlC and RmlD were required for the biosynthesis of O-antigen and mutation in either gene resulted in cells with only the lipid A core (Figure 5C). On the contrary, the LPS patterns of rmlA and rmlB mutants showed only a small change in intensity but retained most of the O-antigen chain length distribution, except for the very long chain length band. This might be explained by the existence of RffH and RffG, two isozymes of RmlA and RmlB, in EC958. RffH and RffG, however, were not able to fully compensate for RmlA or RmlB in resistance to serum (Table 2). The inability to produce the very long chain length O-antigen in rmlA and rmlB mutants might be a crucial factor in determining their susceptibility to human serum. Only the rmlC mutant displayed altered sensitivity to SDS, while mutations in rmlABCD did not affect sensitivity to NaCl (Table 2).

The O-antigen operon in E47a has another set of nucleotide sugar biosynthesis genes (fnlABC) for the biosynthesis of UDP-N-acetyl-α-L-fucosamine (UDP-FucNAc) from UDP-N-acetyl-α-D-glucosamine (UDP-GlcNAc). Curiously, EC958 does not have these three genes in the operon or anywhere else in its genome.

O-antigen processing genes.

The O-antigen flippase (Wzx, EC958_2377) and O-antigen polymerase (Wzy, EC958_2374) from EC958 are highly similar to the corresponding genes from E47a. In contrast, both genes share very low similarity to the corresponding genes from MG1655. Mutation of the wzx gene in MG1655 results in the accumulation of high levels of the UndPP-O unit in the cytoplasm [59] and hypersensitivity to several antibiotics and other agents including nalidixic acid, tetracycline, mitomycin C and hydrogen peroxide [60]. Transcriptional analyses of MG1655 responses to the broad-spectrum biocide polyhexamethylene biguanide suggested that Wzx (also known as RfbX) might be involved in cellular stress responses [61]. In contrast to the previous studies, the wzx gene was defined as essential in this study (Figure 4, Table S1). This might be explained by the fact that Cm was used to select for transposon mutants and the wzx mutants might be hypersensitive to this antibiotic. The wzy O-antigen polymerase gene (EC958_2374) was also defined as essential in EC958 and thus was not further characterized. Mutation of the chain length regulator gene (wzz, EC958_2368) showed a non-modal distribution of O-antigen chain length in its LPS pattern (Figure 5C), an observation consistent with previous study describing the role of this gene [62]. The reduction of long chain length O-antigen in this mutant is likely to account for the serum sensitivity shown in our assay (Table 2).

Glycosyltransferases genes.

Glycosyltransferases (GTs) are required to form the glycosidic bonds between sugars in an O-antigen repeat unit. There are now more than 100,000 GTs within 94 families (http://www.cazy.org). Based on the structure of the O25 repeat unit [63], [64], we predicted a requirement for 4 GTs in EC958. Bioinformatic analysis (blastp) confirmed that there are indeed 4 GTs within the O-antigen operon: EC958_2371, EC958_2372, EC958_2375 and EC958_2376. Surprisingly, two of the predicted GTs in EC958 are very different from those in E47a. Whether this is related to the different source of α-L-FucNAc in EC958 due to the lack of fnlABC remains to be investigated. Since the O25 antigen structure is known, we attempted to predict the glycosidic link formed by each GT using a combination of sequence similarity and O-antigen structure comparison as previously described [64]; the results of which are shown in Table 3 and Figure 5B. We were only able to mutate EC958_2371, while EC958_2372 and EC958_2375 were defined as essential in this study. Analysis of our EC958_2371 mutant confirmed that EC958_2371 is required for O-antigen biosynthesis (Figure 5C), serum resistance and resistance to SDS (Table 2).

thumbnail
Table 3. Summary of E. coli EC958 O25b antigen biosynthesis genes.

https://doi.org/10.1371/journal.pgen.1003834.t003

Other serum resistance mechanisms affecting LPS in EC958

LPS gel analysis was performed on all 54 defined mutants to identify genes that contribute to serum resistance by affecting LPS (Figure S3). The normal LPS pattern of EC958 consists of 12 bands including a thick bottom band representing the lipid A-core and an 11-band laddering pattern of lipid A-core bound O-antigen polymers, followed by approximately 6 thick bands of very long O-antigen chain length (Figure 5C). In addition to the 6 genes involved in O-antigen biosynthesis mentioned above, the LPS patterns of 20 mutants were altered in comparison to wild-type EC958; 6 of these mutants (waaLKYJ, wecA and rfaH) only produced a lipid A-core (Table 2 and Figure S3).

LPS core biosynthesis genes.

The genes involved in biosynthesis of the LPS outer core in EC958 share strong similarity with those from K-12 MG1655 (waaL and waaGPBIJYK) [reviewed in 65], [66], suggesting that the LPS outer core structure in EC958 is the same as that in MG1655. As expected from the functions of waaGBIJK in MG1655, the waaG mutant showed the smallest size of lipid A-core, indicating that its whole outer core was not linked to the lipid A-inner core. Mutations in waaIJYK also produced an expected LPS pattern of lipid A-core only, as these mutations prevent the linkage of O-antigen to the outer core. The LPS of the waaB and waaP mutants produced O-antigen laddering patterns containing an abnormal lipidA-core band, indicating that although the core structures were changed, these changes still allow the linking of O-antigen to the outer core. WaaL is responsible for the ligation of O-antigen to the outer core and mutation of this gene resulted in the synthesis of a lipid A-core without O-antigen. The EC958 waaL and waaGPBIJYK mutants were sensitive to human serum and SDS (except for the waaB and waaY mutant), while only the waaG, waaI and waaY mutants displayed enhanced sensitivity to NaCl (Table 2).

ECA biosynthesis genes.

Seven serum resistant candidate genes within the ECA biosynthesis wec operon were identified by TraDIS (Table 2). Two of these genes (wzxE and wzyE) were not confirmed by serum assays to confer resistance and no changes were observed in the LPS patterns of these two mutants (Table 2, Figure S3). The wzyE mutant (but not the wzxE mutant) was, however, more sensitive to SDS (Table 2). The LPS patterns of the wzzE mutant revealed changes in the lipid A-core and in the modulation of O-antigen chain length (Figure S3). WzxE has previously been shown to preferentially form a protein complex with WzyE and WzzE for the biosynthesis of ECA over Wzy and Wzz [67]. However, our results suggested that WzzE might also contribute to the chain length regulation of O-antigen in EC958. WecA has been shown previously to be involved in the biosynthesis of O7, O18, O75, and O111 antigen [68] and our results indicate that the same is also true for O25b, thus explaining the role of WecA in serum resistance in EC958. The wecD mutant was highly sensitive to serum, possessed an LPS profile that was altered in both the lipid A-core and the amount of O-antigen, and displayed enhanced sensitivity to SDS (Table 2, Figure S3). Further investigation is required to understand the involvement of WecD in lipid A-core and O-antigen biosynthesis. Finally, none of the ECA biosynthesis mutants were altered in sensitivity to NaCl.

Other genes.

The remaining 8 mutants with altered LPS patterns were all confirmed as serum sensitive (Table 2). They included one mutant, rfaH, which produced only lipid A-core, and 7 mutants with different O-antigen patterns: hyxA (EC958_0460), hyxR (EC958_0461), nagA (N-acetylglucosamine-6-phosphate deacetylase), pgm (phosphoglucomutase), galE (UDP-galactose-4-epimerase), EC958_1112 and wcaF (predicted acetyltransferase involved in colanic acid synthesis [69]). Some of these genes were also associated with resistance to additional stresses; EC958 rfaH and wcaF mutants were sensitive to SDS and NaCl, while nagA and pgm mutants were sensitive to SDS (Table 2).

Novel O-antigen chain length regulators

Mutation of the hyxA and hyxR genes in EC958 resulted in the modulation of O-antigen chain length (Figure 6B). The EC958 hyxA mutant exhibited an increased proportion of O-antigen chain of 2 to 6 units with maximum number at 3–5 units and reduction in very high chain length polymer. The EC958 hyxR mutant had an increased proportion of O-antigen polymer of 2–4 units. The hyxA and hyxR genes are located in a pathogenicity island (PAI-X) consisting of 4 genes (fimX and hyxRAB) as previously described in the UPEC strain UTI89 [70]. The hyxB gene (EC958_0459) has also been named upaB due to its function as an autotransporter [71], and we prefer to maintain this nomenclature. This island is present (in several variations) in 24 out of 50 E. coli genomes across all pathotypes (Figure 6A) and current sequence data suggest that it is exclusive to E. coli.

thumbnail
Figure 6. Genetic context of the hyxA and hyxR genes and their role in O-antigen synthesis.

(A) The conserved location of PAI-X containing hyxAR genes in E. coli genomes of various pathotypes. (B) The LPS pattern of the hyxA and hyxR mutants and their complemented strains, demonstrating a role for the hyxA and hyxR genes in O-antigen chain length regulation. Red boxes highlight differences in LPS patterns compared with that of the wild-type.

https://doi.org/10.1371/journal.pgen.1003834.g006

No functional prediction was found for hyxA, and thus this is the first report of hyxA involvement in serum resistance by regulating O-antigen chain length. The hyxR gene encodes a LuxR-like response regulator that suppresses the nitrosative stress response and contributes to intracellular survival in macrophages by regulating hmpA, which encodes a nitric oxide-detoxifying flavohaemoglobin [70]. The expression of the hyxR gene is regulated through bidirectional phase inversion of its promoter region by the upstream gene fimX, which encodes a tyrosine-like recombinase [70]. It is also worth noting that the contribution of hyxA to serum resistance was greater than hyxR, as demonstrated by the 3-log reduction in viability by the hyxA mutant compared to the hyxR mutant. A serum sensitive phenotype for the hyxR mutant was only observed in mixed competition assays, and both mutants did not exhibit altered sensitivity to SDS or NaCl (Table 2).

Novel serum resistance mechanisms not affecting LPS

A major advantage of whole genome approaches such as TraDIS lies in their power of discovery. Out of 46 genes that define the serum resistome of EC958, 21 (46%) genes were confirmed to be required for serum resistance independent of altered LPS patterns (Table 2; Figure S3). The function of these genes ranged across 10 COG functional categories and included ‘Carbohydrate transport and metabolism’ (4 genes), ‘Cell wall/membrane/envelope biogenesis’ (3), ‘Posttranslational modification, protein turnover, chaperones’ (3) and many others (Table 2). Twelve of these genes were also associated with enhanced sensitivity to SDS and one with enhanced NaCl sensitivity (Table 2)

Of the non-LPS genes required for serum resistance, the most notable was lpp (logFC −10, log difference = 6) (Table 2). The lpp gene encodes one of the most abundant proteins in E. coli and is responsible for the stabilisation and integrity of the bacterial cell envelope [72]. Mutation of the lpp gene results in the formation of outer membrane blebs, leakage of periplasmic enzyme ribonuclease, decreased growth rate in media of low ionic strength or low osmolarity and hypersensitive to toxic compounds [73], [74]. Indeed, the EC958 lpp mutant was more sensitive to SDS, suggesting decreased membrane integrity (Table 2). To the best of our knowledge, this study is the first to show a direct link between Lpp and serum resistance.

Another set of genes notable in our TraDIS analysis include tolQAB, which encode three of the six proteins (YbgC-YbgF-TolQ-R-A-B-Pal) that make up the Tol-Pal system of E. coli cell envelope. The Tol-Pal system is responsible for maintaining the integrity of the outer membrane. TolQRA form an inner membrane complex in which TolQR is necessary for its stability [75]. TolB, a periplasmic protein, connects the inner membrane complex with the peptidoglycan-associated lipoprotein, Pal, which is anchored to the outer membrane [76]. In our study, mutation of the tolA and tolQ genes caused sensitivity to human serum and increased sensitivity to SDS, while the tolB mutant did not (Table 2). Our results demonstrate that the Tol-Pal system is important for resistance to human serum, and thus describe a novel function for this important cell wall complex.

BamB is a lipoprotein that is part of the BamABCD complex. Mutation in bamB results in increased outer membrane permeability, thus enhancing sensitivity to rifampin and dramatically reducing growth on SDS and novobiocin [77]. Our data showed that mutation of the bamB gene in EC958 resulted in increased sensitivity to both human serum and SDS (Table 2).

Our TraDIS experiment also indicated that the modification of lipid A with L-Ara4N was important for serum resistance. Three (arnDEF) of the seven genes involved in the biosynthesis and attachment of L-Ara4N to lipid A-core were identified as part of the serum resistome of EC958 and their role was confirmed by mutagenesis (Table 2). This mechanism is known to confer resistance to polymixin B by preventing its binding to lipid A [Reviewed in 78], [79], [80]. ArnD catalyzes a deformylation step to generate UDP-L-Ara4N before it is transported across the inner membrane by ArnEF [79], [81]. The requirement of ArnDEF for serum resistance indicates that EC958 requires L-Ara4N modification to evade the antimicrobial activity of cationic peptides present in human serum. Interestingly, only the arnD mutation conferred sensitivity to SDS (Table 2), which might suggest a role of UDP-L-Ara4N in maintaining membrane integrity. Further investigation is needed to understand why ArnT, the final enzyme required for transferring the L-Ara4N residue to the 4′-phosphate group of lipid A-core, was not identified in our TraDIS-defined serum resistome.

We also identified three genes encoding catabolic enzymes that contributed to the serum resistance phenotype of EC958 (gmm, pgi and fbp) and confirmed their role by mutagenesis (Table 2). Of these genes, only the pgi mutant displayed enhanced sensitivity to SDS (Table 2). Gmm is a GDP-mannose mannosyl hydrolase capable of hydrolyzing both GDP-mannose and GDP-glucose [82]. This enzyme contributes to the biosynthesis of GDP-fucose, a component of colanic acid, possibly by influencing the concentration of GDP-mannose or GDP-glucose in the cell and thus regulating cell wall biosynthesis [82]. Both Pgi and Fbp catalyze the production of D-fructose-6-phosphate from β-D-glucose-6-phosphate and fructose-1,6-bisphosphate, respectively [83], [84]. D-fructose-6-phosphate is a precursor for the biosynthesis of UDP-GlcNAc, which in turn is required for peptidoglycan, lipid A and ECA biosynthesis. Thus, these three enzymes may catalyse key reactions that, if disrupted, could adversely affect the downstream biosynthesis of cell surface components including colanic acid, peptidoglycan, lipid A and ECA.

An ompA mutant is protected from serum killing when present at a low proportion in a mixed bacterial population

Of all the chromosomal genes previously attributed to serum resistance, the only gene that was not identified in our TraDIS screen was ompA. To examine this further we constructed an EC958 ompA mutant and indeed observed it was sensitive to killing by human serum (Figure 7). One way to explain this discrepancy is that the phenotype of an ompA mutant could be complemented in trans by other ompA-intact bacteria in a mixed population such as the mutant library. In fact, OmpA inhibits serum-mediated killing by binding to C4b-binding protein (C4BP) to prevent the activation of C3b via the classical complement pathway [85], and OmpA is known to be released from E. coli cells when treated with serum [86]. In our mutant library, ompA mutants only accounted for approximately 0.02% of the total bacterial cells, and thus we hypothesized that the release of OmpA from 99.98% of the cells, when treated with serum, provided OmpA in trans to complement ompA mutants. We tested this hypothesis by mixing the ompA mutant with wild-type EC958 at various ratios and indeed showed that ompA mutants were protected from serum killing if the proportion of ompA mutants was less than 15% (Figure 7). This result strongly suggests that in trans complementation of OmpA prevents the identification of ompA as a serum resistance gene in our assay.

thumbnail
Figure 7. In trans complementation of wild-type EC958 protects ompA mutants from serum killing.

Wild-type EC958 and an EC958 ompA mutant were mixed at different ratios, exposed to human serum for 90 minutes, and the EC958 ompA mutant was examined for its sensitivity to human serum (expressed as the −log10 CFU difference between time 90 and time 0). The red line indicates the threshold below which a strain was considered resistant to serum killing. The EC958 ompA mutant was protected from killing when present at a percentage of 15% or less in the mixed bacterial suspension.

https://doi.org/10.1371/journal.pgen.1003834.g007

Complementation of selected mutants restores serum resistance

To further demonstrate the function of non-LPS genes in serum resistance, we selected three orphan genes (acnB, greA and fbp) that do not belong to an operon to perform genetic complementation. For these experiments, the selected gene was amplified by PCR, cloned into the low copy number plasmid pSU2718G, and transformed into the respective mutant strain for complementation. In each case, the phenotype of the complemented strain exactly matched that of wild-type EC958 (Table 4). Taken together, these results confirm the role of acnB, greA and fbp in serum resistance and provide a further layer of evidence to support the use of techniques such as TraDIS in functional gene discovery.

Discussion

The rapid advancement of new sequencing technologies has created novel opportunities to interrogate biological systems that were not previously possible. TraDIS was first described as a method that combined high-density mutagenesis with Illumina next generation sequencing technology to study the essential genes of S. Typhi and the conditional essential genes required for survival in bile [42]. The increased data output afforded by next generation DNA sequencing is particularly useful and cost-effective for small bacterial genomes, however it presents technical and bioinformatical challenges for applications such as TraDIS that utilize low complexity DNA libraries. Here we present the application of a modified version of TraDIS that is amenable to multiplexing using the Illumina HiSeq 2000 platform, and we demonstrate its effectiveness by using it to define the essential gene repertoire and the serum resistome of a multidrug resistant strain from the globally disseminated E. coli ST131 lineage.

The multiplexed TraDIS protocol utilizes a newly designed custom oligonucleotide in the library enrichment step of the Illumina library preparation protocol (Figure S1). This oligonucleotide incorporates the Illumina sequencing primer-binding site into transposon specific DNA fragments, enabling the use of the standard Illumina sequencing primer and eliminating the need to design and optimize another sequencing primer for each new transposon sequence. The 6-bp barcode immediately after the sequencing primer-binding site allows 12-sample multiplexing within one lane. The use of 12 barcodes at the first 6 nucleotides of read 1 increased the complexity of the library compared with the original TraDIS protocol, thus reducing data loss due to mis-identification of clusters [87], [88]. However, the number of useable reads from our sequencing runs was still low (15–20% of total reads). We believe non-specific amplification in the enrichment step was the main cause, and further optimization of the enrichment PCR conditions is required. Similar approaches combining transposon mutagenesis with high-throughput sequencing (Tn-seq [45], INSeq [43], HITS [44]) have also been used to address different scientific questions, including the identification of essential genes and genes associated with enhanced fitness in specific growth conditions [41], [42], [45], determination of niche-specific essential genes [43], [44], [89], identification of genes associated with tolerance to various agents/conditions [42], [90] and many other applications as reviewed elsewhere [91], [92]. In terms of insertion density, we achieved 502,068 independent insertion sites with a density of 1/10 bp (i.e. one insertion every 10 bp), which is comparable with the work by Christen et al. (1/8 bp in C. crescentus) [41], Langridge et al. (1/13 bp in Salmonella Typhi) [42] and Barquist et al. (1/9 bp in Salmonella Typhimurium) [46].

The identification of the essential gene set for a single organism is challenging due to several factors, including the presence of transposon insertion cold spots (i.e. regions of low transposon insertion frequency), the difficulty in distinguishing mutations that prevent growth from those that severely reduce growth rate, pre-existing gene duplications and the specific growth conditions used in the experiment [93][95]. In this study, we define essential genes as those genes that, when mutated by transposon insertion, either prevent or severely attenuate growth on LB agar media supplemented with 30 µg/ml Cm at 37°C. The cut-off value to determine whether a gene is essential was defined as the intercept of two distributions of the insertion index of each gene: the exponential distribution representing essential genes and the gamma distribution representing non-essential genes [42]. This means that our essential genes also include those genes that can tolerate insertions but were severely attenuated in the input pool. Out of 315 essential genes, 64 genes had no transposon insertions, 178 genes had 1 to 5 transposon insertions and 73 genes had more than 5 transposon insertions (Table S1).

The high density of insertion sites achieved in our study provided reliable data for the identification of essential genes within the EC958 genomes with a minimal probability of false positive calls due to transposon insertion cold spots. The identification of essential genes has previously been performed using several approaches in E. coli K-12 (strains MG1655 and W3110) [96], [97]. Baba et al. generated null mutations by lambda-red recombination in 3985 E. coli W3110 genes (the Keio library), but were unable to mutate 303 candidate essential genes [97]. This set of essential genes was further consolidated by manual literature review on the EcoGene website (www.ecogene.org), which reduced the set to 289 genes. Of the 315 essential genes identified for EC958 in this study, 231 genes (73%) matched those previously described in the EcoGene list (Table S1). There were 84 essential genes specific for EC958, twenty-four of which do not have homologs in the MG1655 genome. In contrast, 58 genes previously identified as essential for E. coli K-12 were either not present in EC958 or not identified in our analysis. The majority of essential genes in EC958 are conserved with 91% of the genes present in more than 90% of E. coli complete genomes available.

In this study, we provided two layers of evidence for the role of each serum resistance gene: by simultaneously assaying a large mutant library and by generation of defined mutants for independent phenotypic testing. Indeed, using defined mutagenesis we were able to confirm a role for 46 of the 56 genes identified by TraDIS in serum resistance. To the best of our knowledge, this represents the first large scale follow-up of TraDIS data in this manner and highlights the effectiveness of the technique in large-scale functional genomics. Our study also revealed that trans complementation of specific mutants can occur in a large mutant library population, as demonstrated by our findings with an ompA mutant. We also demonstrated complete complementation of mutants containing deletions in the acnB, greA and fbp genes, corroborating their novel role in serum resistance independent of LPS alterations. Finally, we provided further insight into the mechanistic action of the serum resistance genes identified in EC958 by examining the survival of the respective mutants to outer membrane stresses that affect antimicrobial access and osmotic potential.

The search for genetic determinants of serum resistance in bacteria has been ongoing since the 1970s [98]. Our current understanding of the mechanisms that promote bacterial resistance to human serum include a role for surface structures such as O antigens, K antigens, outer membrane proteins (OmpA), and plasmid-encoded proteins (TraT, Iss) [37][39]; notably, however, not all of these mechanisms are required for resistance in a single bacterial strain [32], [33]. Our study represents the first report to simultaneously investigate the entire serum resistome of one strain. Our results demonstrated that both the lipid A-core and O25 antigen are crucial for serum resistance in EC958, while K antigen does not contribute to serum resistance. Out of 54 defined mutants investigated, half had changes in their LPS gel patterns, all of which resulted in serum sensitivity. In contrast, none of the K capsular biosynthesis genes were identified in our TraDIS screen. This result is similar to that reported for the O75:K5 UPEC strain GR-12, where alterations in O75 LPS affected serum resistance more than a K5 null mutation [34]. It is likely, however, that there are strain-specific differences for the role of O antigen and K capsule in serum resistance, and that this reflects differences in the make-up of these structures. For example, previous analysis of an E. coli O4:K54:H5 blood isolate revealed that the K54 antigen contributes more to serum resistance than the O4 antigen [99]. Other K antigens such as the K1 and K2 capsules have also been shown to play an important role in serum resistance [100], [101]. The K antigen expressed by EC958 has not been typed but genomic analysis shows that EC958 has a group 2 capsular gene cluster that conforms to the conserved structure of this group [Reviewed in 102]. However, region 2, which encodes glycosyltransferases specific for each K type, shares such low similarity with available sequences in the GenBank database that deducing its K type in silico was not possible.

The O25 antigen gene cluster was further characterized using sequence analysis, targeted mutation and LPS profiling. All of the dTDP-α-L-rhamnose biosynthesis genes (rmlCADB) were required for serum resistance. However, EC958 lacks the biosynthesis genes for UDP-FucNAc, a component of O25 antigen unit. If, based on the cross reaction of antiserum against the O25 antigen with O25b expressing cells, we assume that EC958 has the same O-antigen repeat unit as the O25 determined from previous studies [57], [58], then EC958 must possess a novel mechanism for the synthesis or uptake of UDP-FucNAc.

Two additional surface antigens that contribute to serum resistance in EC958 are the enterobacterial common antigen and colanic acid (M antigen). Mutations in five ECA biosynthesis genes rendered EC958 susceptible to serum. While mutation of three of these genes (wecA, wzzE and wecD) affected LPS and sensitivity to SDS, mutation of wecE and wecF did not change LPS (although a wecF mutant was more sensitive to SDS), suggesting that the ECA may be involved in serum resistance, perhaps indirectly via its role in membrane integrity. Our data also suggest the involvement of colanic acid in serum resistance. Three genes encoding for colanic acid biosynthesis (wcaI, gmm and wcaF) were identified by TraDIS. Gmm is most likely to be involved in the biosynthesis of colanic acid [82], while mutation of the wcaI gene did not confer serum resistance. The product of wcaF was predicted to be an acetyltransferase [103] required for colanic acid production [69]. The EC958 wcaF mutant possessed an altered LPS pattern with a reduced amount of O-antigen (especially very long chain length O-antigen) and was sensitive to both SDS and high osmolarity. The enhanced sensitivity of the EC958 wcaF mutant could therefore be explained by a number of factors, including altered LPS, altered colonic acid and reduced overall membrane integrity.

A number of other genes were identified that contributed to serum resistance in an LPS-dependent manner. The gene nagA encodes N-acetylglucosamine-6-phosphate deacetylase, an enzyme important for the metabolism of N-acetyl-D-glucosamine [104]. It catalyzes the first step in producing UDP-GlcNAc, a nucleotide sugar required for ECA, lipid A and peptidoglycan biosynthesis, by deacetylating N-acetylglucosamine-6-phosphate to glucosamine-6-phosphate [Reviewed in 105]. However, NagA is not solely responsible for the production of UDP-GlcNAc because glucosamine-6-phosphate can also be obtained via GlmS from fructose-6-phosphate or taken up from the environment by ManXYZ [105]. Indeed, the LPS banding pattern of the nagA mutant was different to that of the parent strain (i.e. it possessed thicker second and third bands from the bottom of the gel; Figure S3), suggesting its enhanced sensitivity phenotype may be associated with a predominantly shorter O antigen.

Pgm is a phosphoglucomutase that catalyses the reversible conversion of glucose-1-phosphate to glucose-6-phosphate, an important step in galactose and maltose catabolism [106]. A pgm mutant has several phenotypes; it is defective in swimming and swarming mobility [107], it possesses an aberrant (shorter and wider) cell morphology, is sensitive to detergents [108] and it stains blue with iodine when grown in the presence of galactose [106]. An EC958 pgm mutant produced little full length O-antigen; the majority of its LPS condensed into a thick band of incomplete lipid A-core and a thin clear band of lipid A-core plus one unit of O-antigen. This feature is consistent with the high serum and SDS sensitivity phenotypes observed for this mutant.

GalE is a well-studied enzyme that catalyzes the interconversion of UDP-galactose and UDP-glucose [Reviewed in 109]. Both nucleotide sugars are required for colanic acid biosynthesis. Furthermore, UDP-glucose is used in three steps to synthesize the LPS outer core (catalyzes by WaaG, WaaI and WaaJ). LPS patterns of the galE mutant exhibited a very thick band of lipid A-core, suggesting that the lipid A-core in this strain has multiple sizes. This may be explained by the limiting effect of UDP-glucose in the three steps involved in its incorporation into the outer core. UDP-glucose can also be synthesized by GalU from glucose-1-phosphate [110], which may explain why an EC958 galE mutant could still make LPS (Figure S3). Despite being able to make LPS, however, the galE mutant was sensitive to human serum. Whether this sensitivity can be attributed to the effect a galE mutation has on LPS or colanic acid remains to be determined.

We have demonstrated the successful application of multiplexed TraDIS for a functional genomics study targeted at E. coli EC958, a prototype strain from the globally disseminated and multidrug resistant E. coli ST131 lineage. This approach enabled the first description of an essential gene set from an ExPEC strain. Our work has also defined the serum resistome in E. coli EC958. This comprehensive inventory of E. coli EC958 genes that contribute to this phenotype provides a framework for the future characterization of virulence genes in ExPEC as well as other Gram-negative pathogens that cause systemic infection.

Materials and Methods

Ethics statement

Approval for the collection of human blood was obtained from the University of Queensland Medical Research Ethics Committee (2008001123). All subjects provided written informed consent.

Bacterial strains and growth conditions

E. coli EC958 was isolated from the urine of a patient presenting with community UTI in the Northwest region of England and is a representative member of the UK epidemic strain A (PFGE type), one of the major pathogenic lineages causing UTI across the United Kingdom [23]. EC958Δlac, which contained a mutation in the lac operon, was used in competitive assays. This strain had an identical growth rate to wild-type EC958. Strains were routinely cultured at 37°C on solid or in liquid Luria Broth (LB) medium supplemented with the appropriate antibiotics (Cm 30 µg/ml or gentamicin 20 µg/ml) unless indicated otherwise.

Generation of miniTn5-Cm mutant library

A miniTn5-Cm transposon containing a Cm cassette flanked by Tn5 mosaic ends (sequence from Epicenter) was PCR amplified from pKD3 plasmid DNA (NotI digested) using primers 2279 5′- CTGTCTCTTATACACATCTcacgtcttgagcgattgtgtagg-3′ and 2280 5′- CTGTCTCTTATACACATCTgacatgggaattagccatggtcc-3′. The PCR reactions were performed using Phusion High-Fidelity DNA polymerase (New England BioLabs). The amplicon was purified using the QIAGEN MinElute PCR purification kit before being phosphorylated using T4 polynucleotide kinase (New England BioLabs) and subjected to the final purification step. A total of at least 800 ng of this miniTn5-Cm transposon DNA was incubated in an 8 µl reaction containing 4 µl of EZ-Tn5 transposase (Epicenter Biotechnologies) at 37°C for 1 h then stored at −20°C.

Bacterial cells were prepared for electroporation as previously described [42]. Briefly, cells were grown in 2×TY broth to an OD600 of 0.3–0.5, then harvested and washed three times in 0.5× volume of 10% cold glycerol before being resuspended in a 1/1000× volume of 10% cold glycerol and kept on ice. A volume of 60 µl cells was mixed with 0.2 µl of transposomes and electroporated in a 2 mm cuvette using a BioRad GenePulser set to 2.5 kV, 25 µF and 200Ω. Cells were resuspended in 1 mL SOC medium and incubated at 37°C for 2 hours, then spread on LB agar plates supplemented with Cm 30 µg/mL. After incubation overnight at 37°C, the total number of colonies was estimated by counting a proportion from multiple plates. Chloramphenicol resistant colonies were resuspended in sterilised LB broth using a bacteriological spreader before adding sterile glycerol to 15% total volume and stored in −80°C. Each batch of mutants contained an estimated 32,000 to 180,000 mutants. The final library of 1 million mutants was created by pooling 11 mutant batches, resulting in a cell suspension of 2×1011 CFU/ml.

Transposon library screening in human serum

Freshly pooled human serum was collected from at least two healthy individuals on the day of the experiment. Ten milliliters of blood was collected from each person and centrifuged at 4000 rpm for 10 minutes to collect the serum. Approximately 2×108 viable mutants were incubated in 1 ml of 50% freshly pooled human serum in LB broth at 37°C for 90 minutes. The control samples were prepared the same way but were incubated with inactivated serum (Millipore) instead of fresh serum. Both control and test experiments were performed in duplicate. The cells were then washed twice with sterile 1×PBS to remove serum, transferred to 100 ml LB broth and allowed to grow at 37°C with 250 rpm shaking for 4 hours. The genomic DNA was then extracted from 5 ml of each culture using Qiagen 100-G genomic tips.

Multiplexed TraDIS

Genomic DNA was standardized to 3.6 µg in a volume of 120 µl before being sheared by Covaris S2 according to the Illumina TruSeq Enrichment gel-free method (TruSeq DNA sample preparation v2 guide). The subsequent steps of DNA end repair, DNA end adenylation and adapter ligation were also done following the Illumina TruSeq v2 instructions. The adapter-ligated fragments containing transposon insertion sites were enriched using a transposon-specific indexing forward primer and the Illumina reverse primer Index 1 (Table S5) at 500 nM each per reaction. This 89 bp forward primer binds specifically to miniTn5-Cm transposon (25 bp) and carries 64 bp overhang which includes 6 bp index sequence, 33 bp binding site for Illumina read 1 sequencing primer and 23 bp P5 sequence for binding to the flowcell. The enrichment step was done using the KAPA Library Amplification kit (KAPA Biosystems) at an annealing temperature of 60°C for 22 cycles. The KAPA Library Quantification kit was used to measure the concentration of DNA fragments in the enriched library. Twelve libraries from 12 samples were pooled to equimolar concentrations to achieve a cluster density of 850 K/mm2 when 10 nM of library pool was loaded onto the flowcell. The 12-plex pool was loaded on 3 lanes of TruSeq v2 and 3 lanes of TruSeq v3 flowcells for sequencing using a 100 cycles, paired-end protocol to access the reproducibility and read quantity among lanes and flowcells. The data from six samples (Figure 1A) were presented in this study. The TraDIS sequence data from this study was deposited on the Sequence Read Archive (SRA) under the BioProject number PRJNA189704.

Analysis of nucleotide sequence data

Sequence reads from the FASTQ files were split according to twelve 6 bp index sequences combined with the 37 bp transposon-specific sequence using fastx_barcode_splitter.pl (total length of 43 bp as barcodes, allowing for 2 mismatches) (FASTX-Toolkit version 0.0.13, http://hannonlab.cshl.edu/fastx_toolkit/index.html). The barcode matching reads were trimmed off the 43 bp barcode at the 5′end and 25 bp of potential low quality at the 3′ end, resulting in high quality sequence reads of 31 bp in length that were used to map to the EC958 chromosome (PRJEA61443) by Maq version 0.7.1 [111]. Subsequent analysis steps were carried out as previously described [42] to calculate the number of sequence reads (raw read counts) and the number of different insertion sites for every gene, which were then used to estimate the threshold to identify essential genes. The read counts and insertion sites were visualized using Artemis version 13.0 [112]. The circular genome diagram was generated by CGView [113] and linear genetic comparison was illustrated using Easyfig version 2.1 [114].

Statistical analyses

We identified genes required for survival in human serum by comparing the differences in read abundance of each gene between the inactivated serum control and active serum test samples using the Bioconductor package edgeR (version 2.6.10) [48]. The raw read counts from two biological replicates of each treatment were loaded into the edgeR package (version 2.6.12) using the R environment (version 2.15.1). Genes that have very low read counts in all the samples (essential genes) were removed from further analysis. The composition bias in each sequence library was normalized using the trimmed mean of M value (TMM) method [115]. We then used the quantile-adjusted conditional maximum likelihood (qCML) for negative binomial models to estimate the dispersions (biological variation between replicates) and to carry out the exact tests for determining genes with significantly lower read counts in the test samples compared to the control samples [116], [117]. Stringent criteria of log fold-change (logFC) ≤−1 and false discovery rate ≤0.001 were chosen to define a list of the most significant genes for further investigation by phenotypic assays.

Molecular methods

Chromosomal DNA purification, PCR and DNA sequencing of PCR products was performed as previously described [118]. Defined mutations were made using the λ-Red recombinase method with some modifications [22], [47]. In brief, the final PCR products were fused and amplified from three fragments containing two 500-bp homologous regions flanking the gene of interest and a Cm cassette from pKD3 plasmid (see Table S5 for list of primers). The fused PCR products were then electroporated into EC958 harbouring a gentamicin resistant plasmid carrying the λ-Red recombinase gene. Mutants were then selected and confirmed by sequencing. Complementation was done by cloning the gene of interest into a gentamicin resistant derivative of pSU2718 [119] at BamHI-XbaI cut sites (primers listed in Table S5). The construct was then transformed into the respective mutant and induced using 1 mM IPTG before and during phenotypic assays.

Serum resistance assay

Overnight bacterial cultures were washed in phosphate buffered saline (PBS) and then standardized to an OD600 of 0.8. Equal volumes (50 µL) of standardized cultures and pooled human sera were mixed and incubated for 90 min at 37°C (in triplicates). Viable counts were performed to estimate the number of bacterial cells prior to serum treatment (t = 0 min) and post serum treatment (t = 90 min). E. coli MG1655 was used as a control as it is completely killed by serum. Serum and PBS only samples served as sterility controls. Competitive serum resistance assays were performed in the same manner, except that a 50∶50 mixture of wild-type (EC958Δlac) and mutant strains were used. Viable counts were performed on MacConkey agar, which allowed the differentiation of EC958Δlac (non-lactose fermenter) and the mutant strains.

SDS and NaCl sensitivity assays (MIC)

The MICs of SDS and NaCl were determined by broth microdilution method as previously described [120]. We used five concentrations for SDS including 0.125%, 0.0625%, 0.031%, 0.016% and 0.008% in LB. For NaCl, the range of concentration was 0.8 M, 0.6 M, 0.5 M, 0.4 M and 0.3 M.

LPS gel assay

LPS was extracted from bacterial strains and LPS patterns were determined by Tricine-SDS Polyacrylamide gel electrophoresis (TSDS-PAGE) and visualized by silver staining as previously described [121], [122].

Supporting Information

Figure S1.

Diagram demonstrating the use of the custom primer with four functional regions for PCR enrichment of transposon derived fragments.

https://doi.org/10.1371/journal.pgen.1003834.s001

(PDF)

Figure S2.

Results of edgeR analysis showing logFC of all genes and highlighting genes that satisfy the stringent cut-off for serum resistome.

https://doi.org/10.1371/journal.pgen.1003834.s002

(PDF)

Figure S3.

LPS patterns of all defined mutants generated in this study. The LPS pattern of wild-type strain EC958 (WT) was used as reference to identify alteration in LPS of other mutants.

https://doi.org/10.1371/journal.pgen.1003834.s003

(PDF)

Table S1.

Summary of transposon insertions in each CDS within the EC958 chromosome.

https://doi.org/10.1371/journal.pgen.1003834.s004

(XLSX)

Table S2.

Sequence comparison of EC958 essential genes and genes from 50 E. coli complete genomes.

https://doi.org/10.1371/journal.pgen.1003834.s005

(XLSX)

Table S3.

Essential genes unique to EC958.

https://doi.org/10.1371/journal.pgen.1003834.s006

(XLSX)

Table S4.

Results of edgeR analysis comparing read counts between control and test samples to identify the serum resistome.

https://doi.org/10.1371/journal.pgen.1003834.s007

(XLSX)

Acknowledgments

We would like to thank Gemma Langridge, Keith Turner, Sabine Eckert, Daniel Turner, Philip Hugenholtz, Nicholas West and Brian Forde for expert advice and Leopold Parts for statistical analysis in R.

Author Contributions

Conceived and designed the experiments: MDP SS MT SAB MAS. Performed the experiments: MDP KMP SS SWL LPA MESA VMM. Analyzed the data: MDP DGM SAB MAS. Contributed reagents/materials/analysis tools: VMM MU SAB. Wrote the paper: MDP SS DGM MT MU SAB MAS.

References

  1. 1. Nicolas-Chanoine MH, Blanco J, Leflon-Guibout V, Demarty R, Alonso MP, et al. (2008) Intercontinental emergence of Escherichia coli clone O25:H4-ST131 producing CTX-M-15. J Antimicrob Chemother 61: 273–281.
  2. 2. Pitout JD, Gregson DB, Campbell L, Laupland KB (2009) Molecular characteristics of extended-spectrum-beta-lactamase-producing Escherichia coli isolates causing bacteremia in the Calgary Health Region from 2000 to 2007: emergence of clone ST131 as a cause of community-acquired infections. Antimicrob Agents Chemother 53: 2846–2851.
  3. 3. Johnson JR, Johnston B, Clabots C, Kuskowski MA, Castanheira M (2010) Escherichia coli sequence type ST131 as the major cause of serious multidrug-resistant E. coli infections in the United States. Clin Infect Dis 51: 286–294.
  4. 4. Peirano G, Pitout JD (2010) Molecular epidemiology of Escherichia coli producing CTX-M beta-lactamases: the worldwide emergence of clone ST131 O25:H4. Int J Antimicrob Agents 35: 316–321.
  5. 5. Sidjabat HE, Derrington P, Nimmo GR, Paterson DL (2010) Escherichia coli ST131 producing CTX-M-15 in Australia. J Antimicrob Chemother 65: 1301–1303.
  6. 6. Rogers BA, Sidjabat HE, Paterson DL (2011) Escherichia coli O25b-ST131: a pandemic, multiresistant, community-associated strain. J Antimicrob Chemother 66: 1–14.
  7. 7. Ender PT, Gajanana D, Johnston B, Clabots C, Tamarkin FJ, et al. (2009) Transmission of an extended-spectrum-beta-lactamase-producing Escherichia coli (sequence type ST131) strain between a father and daughter resulting in septic shock and Emphysematous pyelonephritis. J Clin Microbiol 47: 3780–3782.
  8. 8. Johnson JR, Anderson JT, Clabots C, Johnston B, Cooperstock M (2010) Within-household sharing of a fluoroquinolone-resistant Escherichia coli sequence type ST131 strain causing pediatric osteoarticular infection. Pediatr Infect Dis J 29: 473–475.
  9. 9. Coque TM, Novais A, Carattoli A, Poirel L, Pitout J, et al. (2008) Dissemination of clonally related Escherichia coli strains expressing extended-spectrum beta-lactamase CTX-M-15. Emerg Infect Dis 14: 195–200.
  10. 10. Johnson JR, Menard M, Johnston B, Kuskowski MA, Nichol K, et al. (2009) Epidemic clonal groups of Escherichia coli as a cause of antimicrobial-resistant urinary tract infections in Canada, 2002 to 2004. Antimicrob Agents Chemother 53: 2733–2739.
  11. 11. Johnson JR, Tchesnokova V, Johnston B, Clabots C, Roberts PL, et al. (2013) Abrupt Emergence of a Single Dominant Multidrug-Resistant Strain of Escherichia coli. J Infect Dis 207: 919–928.
  12. 12. Peirano G, Schreckenberger PC, Pitout JD (2011) Characteristics of NDM-1-producing Escherichia coli isolates that belong to the successful and virulent clone ST131. Antimicrob Agents Chemother 55: 2986–2988.
  13. 13. Morris D, Boyle F, Ludden C, Condon I, Hale J, et al. (2011) Production of KPC-2 carbapenemase by an Escherichia coli clinical isolate belonging to the international ST131 clone. Antimicrob Agents Chemother 55: 4935–4936.
  14. 14. Morris D, McGarry E, Cotter M, Passet V, Lynch M, et al. (2012) Detection of OXA-48 carbapenemase in the pandemic clone Escherichia coli O25b:H4-ST131 in the course of investigation of an outbreak of OXA-48-producing Klebsiella pneumoniae. Antimicrob Agents Chemother 56: 4030–4031.
  15. 15. Totsika M, Moriel DG, Idris A, Rogers BA, Wurpel DJ, et al. (2012) Uropathogenic Escherichia coli mediated urinary tract infection. Curr Drug Targets 13: 1386–1399.
  16. 16. Ulett GC, Totsika M, Schaale K, Carey AJ, Sweet MJ, et al. (2013) Uropathogenic Escherichia coli virulence and innate immune responses during urinary tract infection. Curr Opin Microbiol 16: 100–107.
  17. 17. Clermont O, Lavollay M, Vimont S, Deschamps C, Forestier C, et al. (2008) The CTX-M-15-producing Escherichia coli diffusing clone belongs to a highly virulent B2 phylogenetic subgroup. J Antimicrob Chemother 61: 1024–1028.
  18. 18. Johnson JR, Porter SB, Zhanel G, Kuskowski MA, Denamur E (2012) Virulence of Escherichia coli clinical isolates in a murine sepsis model in relation to sequence type ST131 status, fluoroquinolone resistance, and virulence genotype. Infect Immun 80: 1554–1562.
  19. 19. Lavigne JP, Vergunst AC, Goret L, Sotto A, Combescure C, et al. (2012) Virulence potential and genomic mapping of the worldwide clone Escherichia coli ST131. PLoS One 7: e34294.
  20. 20. Gibreel TM, Dodgson AR, Cheesbrough J, Bolton FJ, Fox AJ, et al. (2012) High metabolic potential may contribute to the success of ST131 uropathogenic Escherichia coli. J Clin Microbiol 50: 3202–3207.
  21. 21. Vimont S, Boyd A, Bleibtreu A, Bens M, Goujon JM, et al. (2012) The CTX-M-15-producing Escherichia coli clone O25b: H4-ST131 has high intestine colonization and urinary tract infection abilities. PLoS One 7: e46547.
  22. 22. Totsika M, Beatson SA, Sarkar S, Phan MD, Petty NK, et al. (2011) Insights into a Multidrug Resistant Escherichia coli Pathogen of the Globally Disseminated ST131 Lineage: Genome Analysis and Virulence Mechanisms. PLoS One 6: e26578.
  23. 23. Lau SH, Kaufmann ME, Livermore DM, Woodford N, Willshaw GA, et al. (2008) UK epidemic Escherichia coli strains A-E, with CTX-M-15 beta-lactamase, all belong to the international O25:H4-ST131 clone. J Antimicrob Chemother 62: 1241–1244.
  24. 24. Totsika M, Kostakioti M, Hannan TJ, Upton M, Beatson SA, et al. (2013) A FimH Inhibitor Prevents Acute Bladder Infection and Treats Chronic Cystitis Caused by Multidrug-Resistant Uropathogenic Escherichia coli ST131. J Infect Dis 208(6): 921–8.
  25. 25. Floyd RV, Upton M, Hultgren SJ, Wray S, Burdyga TV, et al. (2012) Escherichia coli-mediated impairment of ureteric contractility is uropathogenic E. coli specific. J Infect Dis 206: 1589–1596.
  26. 26. McCabe WR, Kaijser B, Olling S, Uwaydah M, Hanson LA (1978) Escherichia coli in bacteremia: K and O antigens and serum sensitivity of strains from adults and neonates. J Infect Dis 138: 33–41.
  27. 27. Martínez JL, Cercenado E, Pérez-Díaz JC, Baquero F (1986) Multifactorial determination of systemic invasivity in Escherichia coli. FEMS Microbiology Letters 37: 259–261.
  28. 28. Jacobson SH, ÖStenson C-G, Tullus K, Brauner A (1992) Serum resistance in Escherichia coli strains causing acute pyelonephritis and bacteraemia. APMIS 100: 147–153.
  29. 29. Opal S, Cross A, Gemski P (1982) K antigen and serum sensitivity of rough Escherichia coli. Infect Immun 37: 956–960.
  30. 30. Feingold DS (1969) The serum bactericidal reaction. IV. Phenotypic conversion of Escherichia coli from serum-resistance to serum-sensitivity by diphenylamine. J Infect Dis 120: 437–444.
  31. 31. Gemski P, Cross AS, Sadoff JC (1980) K1 antigen-associated resistance to the bactericidal activity of serum. FEMS Microbiology Letters 9: 193–197.
  32. 32. Stawski G, Nielsen L, ØRskov F, ØRskov IDA (1990) Serum sensitivity of a diversity of Escherichia coli antigenic reference strains. APMIS 98: 828–838.
  33. 33. Cross AS, Kim KS, Wright DC, Sadoff JC, Gemski P (1986) Role of lipopolysaccharide and capsule in the serum resistance of bacteremic strains of Escherichia coli. J Infect Dis 154: 497–503.
  34. 34. Burns SM, Hull SI (1998) Comparison of loss of serum resistance by defined lipopolysaccharide mutants and an acapsular mutant of uropathogenic Escherichia coli O75:K5. Infect Immun 66: 4244–4253.
  35. 35. Porat R, Mosseri R, Kaplan E, Johns MA, Shibolet S (1992) Distribution of polysaccharide side chains of lipopolysaccharide determine resistance of Escherichia coli to the bactericidal activity of serum. J Infect Dis 165: 953–956.
  36. 36. Weiser JN, Gotschlich EC (1991) Outer membrane protein A (OmpA) contributes to serum resistance and pathogenicity of Escherichia coli K-1. Infect Immun 59: 2252–2258.
  37. 37. Moll A, Manning PA, Timmis KN (1980) Plasmid-determined resistance to serum bactericidal activity: a major outer membrane protein, the traT gene product, is responsible for plasmid-specified serum resistance in Escherichia coli. Infect Immun 28: 359–367.
  38. 38. Montenegro MA, Bitter-Suermann D, Timmis JK, Aguero ME, Cabello FC, et al. (1985) traT gene sequences, serum resistance and pathogenicity-related factors in clinical isolates of Escherichia coli and other gram-negative bacteria. J Gen Microbiol 131: 1511–1521.
  39. 39. Smith HW (1974) A search for transmissible pathogenic characters in invasive strains of Escherichia coli: the discovery of a plasmid-controlled toxin and a plasmid-controlled lethal character closely associated, or identical, with colicine V. J Gen Microbiol 83: 95–111.
  40. 40. Barondess JJ, Beckwith J (1995) bor gene of phage lambda, involved in serum resistance, encodes a widely conserved outer membrane lipoprotein. J Bacteriol 177: 1247–1253.
  41. 41. Christen B, Abeliuk E, Collier JM, Kalogeraki VS, Passarelli B, et al. (2011) The essential genome of a bacterium. Mol Syst Biol 7: 528.
  42. 42. Langridge GC, Phan MD, Turner DJ, Perkins TT, Parts L, et al. (2009) Simultaneous assay of every Salmonella Typhi gene using one million transposon mutants. Genome Res 19: 2308–2316.
  43. 43. Goodman AL, McNulty NP, Zhao Y, Leip D, Mitra RD, et al. (2009) Identifying genetic determinants needed to establish a human gut symbiont in its habitat. Cell Host Microbe 6: 279–289.
  44. 44. Gawronski JD, Wong SM, Giannoukos G, Ward DV, Akerley BJ (2009) Tracking insertion mutants within libraries by deep sequencing and a genome-wide screen for Haemophilus genes required in the lung. Proc Natl Acad Sci U S A 106: 16422–16427.
  45. 45. van Opijnen T, Bodi KL, Camilli A (2009) Tn-seq: high-throughput parallel sequencing for fitness and genetic interaction studies in microorganisms. Nat Methods 6: 767–772.
  46. 46. Barquist L, Langridge GC, Turner DJ, Phan MD, Turner AK, et al. (2013) A comparison of dense transposon insertion libraries in the Salmonella serovars Typhi and Typhimurium. Nucleic Acids Res 41(8): 4549–64.
  47. 47. Datsenko KA, Wanner BL (2000) One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci U S A 97: 6640–6645.
  48. 48. Robinson MD, McCarthy DJ, Smyth GK (2010) edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139–140.
  49. 49. Whang HY, Neter E (1973) Production of common enterobacterial antigen by members of the family Enterobacteriaceae. Experientia 29: 96–97.
  50. 50. Makela PH, Mayer H (1976) Enterobacterial common antigen. Bacteriol Rev 40: 591–632.
  51. 51. Beutin L, Manning PA, Achtman M, Willetts N (1981) sfrA and sfrB products of Escherichia coli K-12 are transcriptional control factors. J Bacteriol 145: 840–844.
  52. 52. Bailey MJ, Hughes C, Koronakis V (1996) Increased distal gene transcription by the elongation factor RfaH, a specialized homologue of NusG. Mol Microbiol 22: 729–737.
  53. 53. Nagy G, Dobrindt U, Schneider G, Khan AS, Hacker J, et al. (2002) Loss of regulatory protein RfaH attenuates virulence of uropathogenic Escherichia coli. Infect Immun 70: 4406–4413.
  54. 54. Nagy G, Dobrindt U, Grozdanov L, Hacker J, Emody L (2005) Transcriptional regulation through RfaH contributes to intestinal colonization by Escherichia coli. FEMS Microbiol Lett 244: 173–180.
  55. 55. Wang Q, Wang S, Beutin L, Cao B, Feng L, et al. (2010) Development of a DNA microarray for detection and serotyping of enterotoxigenic Escherichia coli. J Clin Microbiol 48: 2066–2074.
  56. 56. Stevenson G, Neal B, Liu D, Hobbs M, Packer NH, et al. (1994) Structure of the O antigen of Escherichia coli K-12 and the sequence of its rfb gene cluster. J Bacteriol 176: 4144–4156.
  57. 57. Kenne L, Lindberg B, Madden JK, Lindberg AA, Gemski P Jr (1983) Structural studies of the Escherichia coli O-antigen 25. Carbohydr Res 122: 249–256.
  58. 58. Fundin J, Weintraub A, Xu J-G, Widmalm G (2003) NMR analysis of the O-antigen polysaccharide from Escherichia coli strain F171. Magnetic Resonance in Chemistry 41: 202–205.
  59. 59. Liu D, Cole RA, Reeves PR (1996) An O-antigen processing function for Wzx (RfbX): a promising candidate for O-unit flippase. J Bacteriol 178: 2102–2107.
  60. 60. Han X, Dorsey-Oresto A, Malik M, Wang JY, Drlica K, et al. (2010) Escherichia coli genes that reduce the lethal effects of stress. BMC Microbiol 10: 35.
  61. 61. Allen MJ, White GF, Morby AP (2006) The response of Escherichia coli to exposure to the biocide polyhexamethylene biguanide. Microbiology 152: 989–1000.
  62. 62. Franco AV, Liu D, Reeves PR (1998) The wzz (cld) protein in Escherichia coli: amino acid sequence variation determines O-antigen chain length specificity. J Bacteriol 180: 2670–2675.
  63. 63. Stenutz R, Weintraub A, Widmalm G (2006) The structures of Escherichia coli O-polysaccharide antigens. FEMS Microbiol Rev 30: 382–403.
  64. 64. Lundborg M, Modhukur V, Widmalm G (2010) Glycosyltransferase functions of E. coli O-antigens. Glycobiology 20: 366–368.
  65. 65. Schnaitman CA, Klena JD (1993) Genetics of lipopolysaccharide biosynthesis in enteric bacteria. Microbiol Rev 57: 655–682.
  66. 66. Raetz CR, Whitfield C (2002) Lipopolysaccharide endotoxins. Annu Rev Biochem 71: 635–700.
  67. 67. Marolda CL, Tatar LD, Alaimo C, Aebi M, Valvano MA (2006) Interplay of the Wzx translocase and the corresponding polymerase and chain length regulator proteins in the translocation and periplasmic assembly of lipopolysaccharide o antigen. J Bacteriol 188: 5124–5135.
  68. 68. Alexander DC, Valvano MA (1994) Role of the rfe gene in the biosynthesis of the Escherichia coli O7-specific lipopolysaccharide and other O-specific polysaccharides containing N-acetylglucosamine. J Bacteriol 176: 7079–7084.
  69. 69. Danese PN, Pratt LA, Kolter R (2000) Exopolysaccharide production is required for development of Escherichia coli K-12 biofilm architecture. J Bacteriol 182: 3593–3596.
  70. 70. Bateman SL, Seed PC (2012) Epigenetic regulation of the nitrosative stress response and intracellular macrophage survival by extraintestinal pathogenic Escherichia coli. Mol Microbiol 83: 908–925.
  71. 71. Allsopp LP, Beloin C, Ulett GC, Valle J, Totsika M, et al. (2012) Molecular characterization of UpaB and UpaC, two new autotransporter proteins of uropathogenic Escherichia coli CFT073. Infect Immun 80: 321–332.
  72. 72. Keseler IM, Collado-Vides J, Santos-Zavaleta A, Peralta-Gil M, Gama-Castro S, et al. (2011) EcoCyc: a comprehensive database of Escherichia coli biology. Nucleic Acids Res 39: D583–590.
  73. 73. Yem DW, Wu HC (1978) Physiological characterization of an Escherichia coli mutant altered in the structure of murein lipoprotein. J Bacteriol 133: 1419–1426.
  74. 74. Vadillo-Rodriguez V, Schooling SR, Dutcher JR (2009) In situ characterization of differences in the viscoelastic response of individual gram-negative and gram-positive bacterial cells. J Bacteriol 191: 5518–5525.
  75. 75. Cascales E, Lloubes R, Sturgis JN (2001) The TolQ-TolR proteins energize TolA and share homologies with the flagellar motor proteins MotA-MotB. Mol Microbiol 42: 795–807.
  76. 76. Walburger A, Lazdunski C, Corda Y (2002) The Tol/Pal system function requires an interaction between the C-terminal domain of TolA and the N-terminal domain of TolB. Mol Microbiol 44: 695–708.
  77. 77. Onufryk C, Crouch ML, Fang FC, Gross CA (2005) Characterization of six lipoproteins in the sigmaE regulon. J Bacteriol 187: 4552–4561.
  78. 78. Nikaido H, Vaara M (1985) Molecular basis of bacterial outer membrane permeability. Microbiol Rev 49: 1–32.
  79. 79. Yan A, Guan Z, Raetz CR (2007) An undecaprenyl phosphate-aminoarabinose flippase required for polymyxin resistance in Escherichia coli. J Biol Chem 282: 36077–36089.
  80. 80. Trent MS, Ribeiro AA, Doerrler WT, Lin S, Cotter RJ, et al. (2001) Accumulation of a polyisoprene-linked amino sugar in polymyxin-resistant Salmonella typhimurium and Escherichia coli: structural characterization and transfer to lipid A in the periplasm. J Biol Chem 276: 43132–43144.
  81. 81. Raetz CR, Reynolds CM, Trent MS, Bishop RE (2007) Lipid A modification systems in gram-negative bacteria. Annu Rev Biochem 76: 295–329.
  82. 82. Frick DN, Townsend BD, Bessman MJ (1995) A novel GDP-mannose mannosyl hydrolase shares homology with the MutT family of enzymes. J Biol Chem 270: 24086–24091.
  83. 83. Fraenkel DG, Levisohn SR (1967) Glucose and gluconate metabolism in an Escherichia coli mutant lacking phosphoglucose isomerase. J Bacteriol 93: 1571–1578.
  84. 84. Fraenkel DG, Horecker BL (1965) Fructose-1, 6-diphosphatase and acid hexose phosphatase of Escherichia coli. J Bacteriol 90: 837–842.
  85. 85. Prasadarao NV, Blom AM, Villoutreix BO, Linsangan LC (2002) A novel interaction of outer membrane protein A with C4b binding protein mediates serum resistance of Escherichia coli K1. J Immunol 169: 6352–6360.
  86. 86. Hellman J, Loiselle PM, Tehan MM, Allaire JE, Boyle LA, et al. (2000) Outer membrane protein A, peptidoglycan-associated lipoprotein, and murein lipoprotein are released by Escherichia coli bacteria into serum. Infect Immun 68: 2566–2572.
  87. 87. Kircher M, Heyn P, Kelso J (2011) Addressing challenges in the production and analysis of illumina sequencing data. BMC Genomics 12: 382.
  88. 88. Krueger F, Andrews SR, Osborne CS (2011) Large scale loss of data in low-diversity illumina sequencing libraries can be recovered by deferred cluster calling. PLoS One 6: e16607.
  89. 89. Crimmins GT, Mohammadi S, Green ER, Bergman MA, Isberg RR, et al. (2012) Identification of MrtAB, an ABC transporter specifically required for Yersinia pseudotuberculosis to colonize the mesenteric lymph nodes. PLoS Pathog 8: e1002828.
  90. 90. Gallagher LA, Shendure J, Manoil C (2011) Genome-scale identification of resistance functions in Pseudomonas aeruginosa using Tn-seq. MBio 2: e00315–00310.
  91. 91. van Opijnen T, Camilli A (2013) Transposon insertion sequencing: a new tool for systems-level analysis of microorganisms. Nat Rev Microbiol 11: 435–442.
  92. 92. Barquist L, Boinett CJ, Cain AK (2013) Approaches to querying bacterial genomes with transposon-insertion sequencing. RNA Biol 10: 1161–9.
  93. 93. Koonin EV (2003) Comparative genomics, minimal gene-sets and the last universal common ancestor. Nat Rev Microbiol 1: 127–136.
  94. 94. Fang G, Rocha E, Danchin A (2005) How essential are nonessential genes? Mol Biol Evol 22: 2147–2156.
  95. 95. Gerdes S, Edwards R, Kubal M, Fonstein M, Stevens R, et al. (2006) Essential genes on metabolic maps. Curr Opin Biotechnol 17: 448–456.
  96. 96. Gerdes SY, Scholle MD, Campbell JW, Balazsi G, Ravasz E, et al. (2003) Experimental determination and system level analysis of essential genes in Escherichia coli MG1655. J Bacteriol 185: 5673–5684.
  97. 97. Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, et al. (2006) Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol 2: 2006 0008.
  98. 98. Taylor PW (1975) Genetical studies of serum resistance in Escherichia coli. J Gen Microbiol 89: 57–66.
  99. 99. Russo TA, Sharma G, Brown CR, Campagnari AA (1995) Loss of the O4 antigen moiety from the lipopolysaccharide of an extraintestinal isolate of Escherichia coli has only minor effects on serum sensitivity and virulence in vivo. Infect Immun 63: 1263–1269.
  100. 100. Buckles EL, Wang X, Lane MC, Lockatell CV, Johnson DE, et al. (2009) Role of the K2 capsule in Escherichia coli urinary tract infection and serum resistance. J Infect Dis 199: 1689–1697.
  101. 101. Leying H, Suerbaum S, Kroll HP, Stahl D, Opferkuch W (1990) The capsular polysaccharide is a major determinant of serum resistance in K-1-positive blood culture isolates of Escherichia coli. Infect Immun 58: 222–227.
  102. 102. Whitfield C (2006) Biosynthesis and assembly of capsular polysaccharides in Escherichia coli. Annu Rev Biochem 75: 39–68.
  103. 103. Stevenson G, Andrianopoulos K, Hobbs M, Reeves PR (1996) Organization of the Escherichia coli K-12 gene cluster responsible for production of the extracellular polysaccharide colanic acid. J Bacteriol 178: 4885–4893.
  104. 104. White RJ (1968) Control of amino sugar metabolism in Escherichia coli and isolation of mutants unable to degrade amino sugars. Biochem J 106: 847–858.
  105. 105. Park JT, Uehara T (2008) How bacteria consume their own exoskeletons (turnover and recycling of cell wall peptidoglycan). Microbiol Mol Biol Rev 72: 211–227.
  106. 106. Adhya S, Schwartz M (1971) Phosphoglucomutase mutants of Escherichia coli K-12. J Bacteriol 108: 621–626.
  107. 107. Inoue T, Shingaki R, Hirose S, Waki K, Mori H, et al. (2007) Genome-wide screening of genes required for swarming motility in Escherichia coli K-12. J Bacteriol 189: 950–957.
  108. 108. Lu M, Kleckner N (1994) Molecular cloning and characterization of the pgm gene encoding phosphoglucomutase of Escherichia coli. J Bacteriol 176: 5847–5851.
  109. 109. Frey PA (1996) The Leloir pathway: a mechanistic imperative for three enzymes to change the stereochemical configuration of a single carbon in galactose. FASEB J 10: 461–470.
  110. 110. Weissborn AC, Liu Q, Rumley MK, Kennedy EP (1994) UTP: alpha-D-glucose-1-phosphate uridylyltransferase of Escherichia coli: isolation and DNA sequence of the galU gene and purification of the enzyme. J Bacteriol 176: 2611–2618.
  111. 111. Li H, Ruan J, Durbin R (2008) Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Res 18: 1851–1858.
  112. 112. Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, et al. (2000) Artemis: sequence visualization and annotation. Bioinformatics 16: 944–945.
  113. 113. Stothard P, Wishart DS (2005) Circular genome visualization and exploration using CGView. Bioinformatics 21: 537–539.
  114. 114. Sullivan MJ, Petty NK, Beatson SA (2011) Easyfig: a genome comparison visualizer. Bioinformatics 27: 1009–1010.
  115. 115. Robinson MD, Oshlack A (2010) A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol 11: R25.
  116. 116. Robinson MD, Smyth GK (2008) Small-sample estimation of negative binomial dispersion, with applications to SAGE data. Biostatistics 9: 321–332.
  117. 117. Robinson MD, Smyth GK (2007) Moderated statistical tests for assessing differences in tag abundance. Bioinformatics 23: 2881–2887.
  118. 118. Allsopp LP, Totsika M, Tree JJ, Ulett GC, Mabbett AN, et al. (2010) UpaH is a newly identified autotransporter protein that contributes to biofilm formation and bladder colonization by uropathogenic Escherichia coli CFT073. Infect Immun 78: 1659–1669.
  119. 119. Martinez E, Bartolome B, de la Cruz F (1988) pACYC184-derived cloning vectors containing the multiple cloning site and lacZ alpha reporter gene of pUC8/9 and pUC18/19 plasmids. Gene 68: 159–162.
  120. 120. Wiegand I, Hilpert K, Hancock RE (2008) Agar and broth dilution methods to determine the minimal inhibitory concentration (MIC) of antimicrobial substances. Nat Protoc 3: 163–175.
  121. 121. West NP, Jungnitz H, Fitter JT, McArthur JD, Guzman CA, et al. (2000) Role of phosphoglucomutase of Bordetella bronchiseptica in lipopolysaccharide biosynthesis and virulence. Infect Immun 68: 4673–4680.
  122. 122. Tsai CM, Frasch CE (1982) A sensitive silver stain for detecting lipopolysaccharides in polyacrylamide gels. Anal Biochem 119: 115–119.