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.
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.
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). 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. 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). 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]. 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. 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). 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%). 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. 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]. 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). 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. 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). 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. 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. 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.
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.
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. 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. 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. 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. 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. 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]. 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. 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. 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. 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 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].
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.
The Serum Resistome of a Globally Disseminated Multidrug Resistant Uropathogenic Escherichia coli Clone