- Research article
- Open Access
Comparative genome analysis identifies few traits unique to the Escherichia coli ST131 H30Rx clade and extensive mosaicism at the capsule locus
BMC Genomics volume 15, Article number: 830 (2014)
E.coli ST131 is a globally disseminated clone of multi-drug resistant E. coli responsible for that vast majority of global extra-intestinal E. coli infections. Recent global genomic epidemiological studies have highlighted the highly clonal nature of this group of bacteria, however there appears to be inconsistency in some phenotypes associated with the clone, in particular capsule types as determined by K-antigen testing both biochemically and by PCR.
We performed improved quality assemblies on ten ST131 genomes previously sequenced by our group and compared them to a new reference genome sequence JJ1886 to identify the capsule loci across the drug-resistant clone H30Rx. Our data shows considerable genetic diversity within the capsule locus of H30Rx clone strains which is mirrored by classical K antigen testing. The varying capsule locus types appear to be randomly distributed across the H30Rx phylogeny suggesting multiple recombination events at this locus, but that this capsule heterogeneity has little to no effect on virulence associated phenotypes in vitro.
Our data provides a framework for determining the capsular genetics of E. coli ST131 and further beyond to ExPEC strains, and highlights how capsular mosaicism may be an important strategy in becoming a successful globally disseminated human pathogen.
Extra-intestinal pathogenic Escherichia coli (ExPEC) infections are one of the leading causes of morbidity in the developed world and are particularly associated with infections of the urinary tract (UTI) and with bacteraemia. In recent years one particular clone of ExPEC has emerged to become a globally dominant cause of human infection, E. coli ST131 ( which is also associated with the emergence and spread of multiple-drug resistance in ExPEC infections via the sustained carriage of the CTX-M-15 extended spectrum beta-lactamase enzyme . Recent work has focussed on elucidating the genomic epidemiology of this group of organisms since the report of the genetically homogeneous nature of clinically unrelated isolates in 2012 . Two independent studies identified that all CTX-M-15 positive isolates belonged to a single expanded clone which emerged some time previous to 2000 [4, 5] and which is now referred to as the H30Rx clade of E. coli ST131 . Both studies show this clade to be monomorphic containing a few dozen SNPs difference in data sets spanning geographical and temporal space.
The genetic architecture of the H30Rx clade was also examined , paying particular attention to virulence associated genes of ExPEC and to mobile genetic elements not found in non-ST131 ExPEC. In general these data suggested no ST131 specific virulence gene repertoire as such, though did highlight the seemingly unique nature of the second flagellar cluster Flag-2, which had been previously identified in E. coli ST131 genomes [3, 6]. Additionally the analysis also highlighted the role of intra-ST131 recombination in shaping the lineage  and identified a recombinant fragment common across ST131 within the capsule locus. Classical capsular typing of a collection of E. coli ST131 isolates, many of which were in the H30Rx clade, has shown high diversity in the biochemical profile of capsule antigens  which seems surprising given the monomorphic nature of the H30Rx clade. There were a total of 7 different K capsule types identified within the forty four ST131 isolates tested, which is in contrast to the vast majority of capsule typing which had been performed previously on E. coli ST131 using PCR based methods and which predominantly identified K2 type capsules via kpsMII primers [8, 9]. Indeed none of the strains biochemically tested were identified as K2 but rather as K100 despite testing K2 positive by PCR .
Given that the comparative genomics performed to date on E. coli ST131 have focussed on virulence associated genes, and the confusing data available to date on the diversity of the capsule locus, we sought to investigate loci uniquely associated with the H30Rx clade of E. coli ST131 using previously published genomes [3, 6, 10, 11]. We analysed a pangenome created from our ST131 genomes against reference non-ST131 ExPEC genomes to identify a small number of loci unique to ST131 dominated by lineage unique phages and the Flag-2 locus. Additionally we provide a genetic architecture for the diversity observed in the capsule locus of ST131, and show extensive genetic and biochemical diversity of the capsule region even within the H30Rx lineage of ST131. The random phylogenetic dispersal of these capsule loci suggests recombination occurs frequently at this region within ST131 and concurs with the previous suggestion that the capsule locus may be coming under strong selective pressure in the lifestyle of E. coli ST131 H30Rx .
Results and Discussion
Identification of genetic loci unique to the E. coliST131 H30Rx clade
Given the focus on virulence associated genes in previous gene content studies, we aimed to determine loci unique to ST131 isolates without bias for functionality of the encoding loci. An ExPEC pan genome was constructed using the blast-score ratio method implemented in LS-BSR  containing twelve ST131 genomes and all available non-ST131 reference genome sequences (Table 1). Using the resulting pan-genome matrix we determined the genetic loci uniquely associated with the ST131 group versus the non-ST131 group using the compare_BSR python script implemented in the LS-BSR package. To define the ST131 group we excluded NA114 on the basis that previous work has suggested the methodology used to assemble the genome has resulted in regions missing from that genome that are present in all other H30Rx strains . We also ran the analysis with SE15 as an ST131 but non-H30Rx strain to determine loci unique to the H30Rx clade to which SE15 does not belong [5, 13]. Our resulting data set identified a total of 150 loci unique to ST131 H30Rx strains in comparison to other ExPEC (Additional file 1), dominated by three phages common across the lineage and which have most probably been acquired by the common ancestral H30Rx progenitor and then maintained in the lineage. Hypothetical proteins dominate the functional category of genes (Figure 1) followed by flagellar associated genes and then a small number of metabolic loci. These metabolic loci correlate with data previously published by our group and add further weight to the assertion that ST131 H30Rx is not a metabolically distinct lineage . The most striking locus with respect to potential biology is the confirmation that the Flag-2 accessory flagella locus is unique to ST131 H30Rx amongst ExPEC strains . Again it is likely that this is ancestral to the H30Rx clade but its acquisition within the larger ST131 lineage is suggestive of a possible role in the formation and dissemination of the H30Rx clade and merits a fuller bacterial genetics investigation of its importance and role in the H30Rx clade. Indeed a fuller genetic investigation of all of the H30Rx loci identified as clade associated may be of merit. A saturated transposon mutant library has been constructed in an H30Rx strain and was utilised to determine the essential gene set for serum resistance . Using such a library to test a wider set of environmental and infection conditions would undoubtedly elucidate if the H30Rx unique loci do indeed play a formative role in the success of the lineage.
Genetic architecture of capsule locus variation in the H30Rx clade
Given the reported variability of capsule loci  and capsular antigen type  in E. coli ST131 H30Rx clade strains, we investigated this locus in more detail. We selected the recently released JJ1886 genome  as our reference given it is the only ST131 genome sequenced and assembled to a standard of quality commensurate with being a high quality genome . Using this reference we re-ordered the contigs of the ST131 genomes previously reported by our group  to ensure the genome architecture was as accurate as possible. We then identified the capsule loci of all of the ST131 genomes at our disposal and created separate embl files for each capsule locus of each strain which we then compared using EasyFig . The comparison of the capsule loci (Figure 2) shows a high degree of diversity between the conserved kpsS and kpsTM regions, with no observable similarity between strains in the variable central genes. Blast analysis of each of the variable central genes in each genetic capsule type present returned no significant hits with any reference E. coli sequences. To ascertain how this genetic architecture reflected upon biochemical typing we determined the K antigen type of each of the genetic capsule types for strains which were in our possession by classical capsule typing (Figure 2), and also overlaid any available capsule type information on the other sequenced strains. Our data shows a correlation between the genetic capsular type and the biochemical typing data, and provides a framework for which to contextualise E. coli capsule types from genomic data. More importantly for this study our data clearly shows significant diversity within the capsule locus in E. coli ST131 H30Rx strains suggestive of frequent and targeted recombination in this region .
To examine this in more detail we created a core genome phylogeny for E. coli ST131 using all the ST131 genome sequences in Table 1 and previously published methodology [13, 18] using SE15 as the root of the phylogeny given its phylogenetic position relative to H30Rx . The resulting phylogenetic tree (Figure 3) confirms that all but one of the strains in our analysis, including those previously sequenced by our group prior to the discovery of H30Rx  do indeed belong the H30Rx clade. More importantly when the capsule loci genetics were superimposed on the phylogenetic tree it clearly demonstrates that the capsule loci are randomly distributed across the phylogeny. The only exceptions to this are the small cluster of strains containing UTI18 which have been previously shown by us to essentially be a single clone . Such a random dispersal of the capsule loci across the phylogenetic tree can only be explained by extensive and targeted recombination events at this discreet location on the genome, suggesting there is some pressure acting on the capsule locus resulting in constant switching of capsule genes as the H30Rx clade evolves. Such extensive recombination has been well characterised in Streptococcus pneumonia where capsule locus switching has been shown to play a significant role in vaccine escape  and in the evolutionary dynamics of densely populated infection foci , however such dynamism in capsular recombination in E. coli is hitherto uncharacterised particularly in such a genetically monomorphic clade as ST131 H30Rx.
Capsule diversity has no obvious effect on virulence associated phenotypes in vitro
Given the observation of extensive recombination at the capsule region of the H30Rx clade we sought to determine any obvious phenotypic effects. We compared the ability of our ST131 strains to form capsules at 25°C and 37°C on LB and CLED agar plates over a 14 day period. There was no association with the capsule loci present in the different H30Rx strains and levels of capsulation morphology on agar plates (Table 2). Classical 96 well plate biofilm formation assays also failed to show any significant pattern between different H30Rx capsular variants. We also conducted in vitro cell adhesion and invasion assays on T24 bladder epithelial cells using both the gentamicin protection assay to quantitate invasion, as well as confocal microscopy using strains carrying a medium copy number GFP+ containing plasmid . As with our other virulence associated phenotypes there was no associated difference between different capsular variants of H30Rx. An identical pattern was also observed when the ability of the strains to survive inside cultured U937 macrophage like cell lines was assayed. Finally we determined the levels of serum resistance in our strains using methods previously employed in our lab [22, 23]. We found that the presence of different capsular variants had no effect on serum resistance and that all of our ST131 strains were totally resistant to serum in the 3 hours used for our assay (Table 2). The importance of serum resistance to E. coli ST131 has been documented and functionally characterised  with several glycosylation associated ORFs identified in as playing an essential role in serum resistance. Capsules have classically been considered as important factors in the ability of E. coli to survive human serum however our data suggests that that capsule type may be less important, and that the extensive capsular recombination demonstrated in the ST131 H30Rx clade has no effect on the ability of these pathogens to survive exposure to human serum. It may be that the capsule variability alters phenotypes important for in vivo environments, and there may be merit to future work investigating differences in infection dynamics between the capsule variants using appropriate surrogate infection models.
E. coli ST131 is now the dominant causative agent of extra-intestinal infection by E. coli in the developed world, and is also heavily responsible for the increase in prevalence in multi-drug resistance in E. coli due to extended carriage of the CTX-M-15 ESBL gene . Recent extensive genomic studies have led to a deep understanding of the phylogeography of this lineage of ExPEC [4, 5] and the discovery of a sub-clade of ST131 which is globally dominant and associated with the CTX-M-15 genotype which has been termed the H30Rx clade . Despite these extensive studies the only efforts at comparative genomics of the ST131 lineage have focussed solely on virulence associated genes and large mobile genetic elements unique to the lineage . Here we present an approach where we created an ExPEC pan-genome and then identified loci uniquely associated with the ST131 H30Rx clade. Our data is further suggestive that at a gene content level this clade is rather unremarkable in comparison to other ExPEC, as recently suggested for the clade at a metabolic level , with the secondary flagellar locus Flag-2 the stand-out region unique to ST131 within ExPEC. This region merits further detailed bacterial genetics analysis to uncover its true importance to the emergence and success of the H30Rx clade. Furthermore our analysis shows a surprising level of diversity within the capsule locus of the H30Rx clade with a phylogenetic distribution highly suggestive of frequent recombination at the locus. This recombination has no obvious detectable effect on virulence associated phenotypes in vitro. Given the level of diversity observed at the capsule locus it is tempting to speculate that there is significant selective pressure occurring at this site during the life cycle of the H30Rx clade, and that frequent recombination allows the clade to subvert that pressure. This has been documented to occur in other capsulated pathogens  and also ties in with previous data from our group showing that ST131 strains did not exhibit inter-species recombination across the E. coli species but that rather recombination events were focussed within the ST131 lineage . Temporal studies of ST131 populations from patients and environmental reservoirs may allow us to determine if capsular switching does occur in vivo and if it is an important mechanism in the successful and prolonged dissemination of this important human pathogen.
Strains and genome data
A list of genomes used in our study is provided in Table 1, and of strains used in our study in Table 2. All strains have been previously characterised [3, 10, 13, 23] with the exception of strain JIE186, which is an Australian ST131 CTX-M-15 strain isolated in 2000, and has been submitted to the ENA under our existing ST131 study accession number ERP001095.
Core and pan genome analysis
We created a pan genome for all ExPEC genomes in Table 1 using LS-BSR . We then used the compare_BSR python script implemented in the LS-BSR package to identify loci unique to genomes belonging to the H30Rx clade, with the exception of NA114 which has been shown to have known H30Rx genes missing from its assembly . The resulting 150 loci identified as H30Rx lineage unique were identified by performing BlastX searches against the genome of JJ1886 .
Identification of capsule loci in ST131 genomes
FastQ sequencing data for all of the ST131 genomes produced by our group were re-assembled using Velvet and PAGIT  and using JJ1886 as a reference genome for contig re-ordering. This allowed us to re-order small contigs to the capsule region. The genomes were then annotated using Prokka  and the capsule regions written to new embl files using Artemis. The capsule encoding regions were visually compared using Easyfig  and variable genes were searched against the non-redundant database by BlastX search.
Classical capsule typing
Serotyping was done according to the method of Ørskov and Ørskov. The K antigen was determined by countercurrent immunoelectrophoresis involving K-specific antisera, except for the K1 and K5 antigens, which were detected using K1- and K5-specific phages .
Whole genome phylogeny
All ST131 genomes were aligned using Mugsy  and a core genome extracted as previously described [13, 18]. Maximum likelihood phylogeny was determined using RaxML  implementing the GTR-gamma model. The resulting phylogeny was visualised using Figtree.
Phenotypic characterisation of strains
Biofilm formation was performed at 37°C in static cultures incubated for 5 days in both LB and BHI broth in a 96 well plate, with 5 wells per strain. Assays were performed on three independent occasions and values are representative values of measured levels of crystal violet retention as measured at A600. Capsule production was determined using a scoring system through testing the ability of each strain to form mucoid colonies in LB agar and in CLED agar plates using two incubation temperatures, 37°C and 25°C. Each strain was tested in triplicate. Ability to invade T24 bladder epithelial cells was performed as previously described , and also performed with strains carrying a GFP+ containing plasmid pMN402  which were visualised using confocal microscopy. Serum resistance assays were performed as described previously , as were U937 macrophage cell line survival assays .
Rogers BA, Sidjabat HE, Paterson DL: Escherichia coli O25b-ST131: a pandemic, multiresistant, community-associated strain. J Antimicrob Chemother. 2011, 66 (1): 1-14.
Johnson J, Tchesnokova V, Johnston B, Clabots C, Roberts P, Billig M, Riddell K, Rogers P, Qin X, Butler-Wu S, Price L, Aziz M, Nicolas-Chanoine M, Debroy C, Robicsek A, Hansen G, Urban C, Platell J, Trott D, Zhanel G, Weissman S, Cookson B, Fang F, Limaye A, Scholes D, Chattopadhyay S, Hooper D, Sokurenko E: Abrupt emergence of a single dominant multidrug-resistant strain of Escherichia coli. J Infect Dis. 2013, 207: 919-928.
Clark G, Paszkiewicz K, Hale J, Weston V, Constantinidou C, Penn CW, Achtman M, McNally A: Genomic and molecular epidemiology analysis of clinical Escherichia coli ST131 isolates suggests circulation of a genetically monomorphic but phenotypically heterogeneous ExPEC clone. J Antimicrob Chemother. 2012, 67: 868-877.
Price L, Johnson J, Aziz M, Clabots C, Johnston B, Tchesnokova V, Nordstrom L, Billig M, Chattopadhyay S, Stegger M, Andersen P, Pearson T, Riddell K, Rogers P, Scholes D, Kahl B, Keim P, Sokurenko E: The epidemic of extended-spectrum-β-lactamase-producing Escherichia coli ST131 is driven by a single highly pathogenic subclone, H30-Rx. MBio. 2013, 4: e00377-13.
Petty N, Ben Zakour N, Stanton-Cook M, Skippington E, Totsika M, Forde B, Phan M, Gomes Moriel D, Peters K, Davies M, Rogers B, Dougan G, Rodriguez-Baño J, Pascual A, Pitout J, Upton M, Paterson D, Walsh T, Schembri M, Beatson S: Global dissemination of a multidrug resistant Escherichia coli clone. Proc Natl Acad Sci U S A. 2014, 111 (15): 5694-5699.
Totsika M, Beatson SA, Sarkar S, Phan MD, Petty NK, Bachmann N, Szubert M, Sidjabat HE, Paterson DL, Upton M, Schembri MA: Insights into a multidrug resistant Escherichia coli pathogen of the globally disseminated ST131 lineage: genome analysis and virulence mechanisms. PLoS One. 2011, 6 (10): e26578-
Olesen B, Hansen D, Nilsson F, Frimodt-Møller J, Leihof R, Struve C, Scheutz F, Johnston B, Krogfelt K, Johnson J: Prevalence and characteristics of the epidemic multiresistant Escherichia coli ST131 clonal group among extended-spectrum beta-lactamase-producing E. coli isolates in Copenhagen, Denmark. J Clin Microbiol. 2013, 51: 1779-1785.
Johnson J, O’Bryan T: Detection of the Escherichia coli group 2 polysaccharide capsule synthesis Gene kpsM by a rapid and specific PCR based assay. J Clin Microbiol. 2004, 42: 1773-1776.
Croxall G, Hale J, Weston V, Manning G, Cheetham P, Achtman M, McNally A: Molecular epidemiology of extraintestinal pathogenic Escherichia coli isolates from a regional cohort of elderly patients highlights the prevalence of ST131 strains with increased antimicrobial resistance in both community and hospital care settings. J Antimicrob Chemother. 2011, 66 (11): 2501-2508.
Andersen P, Stegger M, Aziz M, Contente-Cuomo T, Gibbons H, Keim P, Sokurenko E, Johnson J, Price L: Complete genome sequence of the epidemic and highly virulent CTX-M-15-producing H30-Rx Subclone of Escherichia coli ST131. Genome Announc. 2013, 1: e00988-13.
Avasthi TS, Kumar N, Baddam R, Hussain A, Nandanwar N, Jadhav S, Ahmed N: Genome of multidrug-resistant Uropathogenic Escherichia coli Strain NA114 from India. J Bacteriol. 2011, 193 (16): 4272-4273.
Sahl J, Caporaso J, Rasko D, Keim P: The large-scale blast score ratio (LS-BSR) pipeline: a method to rapidly compare genetic content between bacterial genomes. PeerJ. 2014, 2: e332-
McNally A, Cheng L, Harris SR, Corander J: The evolutionary path to extra intestinal pathogenic, drug resistant Escherichia coli is marked by drastic reduction in detectable recombination within the core genome. Genome Biol Evol. 2013, 5: 699-710.
Alqasim A, Emes R, Clark G, Newcombe J, La Ragione R, McNally A: Phenotypic microarrays suggest Escherichia coli ST131 is not a metabolically distinct lineage of extra-intestinal pathogenic E. coli. PLoS One. 2014, 9: e88374-
Phan M, Peters K, Sarkar S, Lukowski S, Allsopp L, Gomes Moriel D, Achard M, Totsika M, Marshall V, Upton M, Beatson S, Schembri M: The serum resistome of a globally disseminated multidrug resistant uropathogenic Escherichia coli clone. PLoS Genet. 2013, 9: e1003834-
Chain PS, Grafham DV, Fulton RS, Fitzgerald MG, Hostetler J, Muzny D, Ali J, Birren B, Bruce DC, Buhay C, Cole JR, Ding Y, Dugan S, Field D, Garrity GM, Gibbs R, Graves T, Han CS, Harrison SH, Highlander S, Hugenholtz P, Khouri HM, Kodira CD, Kolker E, Kyrpides NC, Lang D, Lapidus A, Malfatti SA, Markowitz V, Metha T, et al: Genomics: genome project standards in a new era of sequencing. Science. 2009, 326 (5950): 236-237.
Sullivan M, Petty N, Beatson S: Easyfig: a genome comparison visualizer. Bioinformatics. 2011, 27: 1009-1010.
Sahl JW, Johnson JK, Harris AD, Phillippy AM, Hsiao WW, Thom KA, Rasko DA: Genomic comparison of multi-drug resistant invasive and colonizing Acinetobacter baumannii isolated from diverse human body sites reveals genomic plasticity. BMC Genomics. 2011, 12 (291): doi:10.1186/1471-2164-12-291
Croucher N, Finkelstein J, Pelton S, Mitchell P, Lee G, Parkhill J, Bentley S, Hanage W, Lipsitch M: Population genomics of post-vaccine changes in pneumococcal epidemiology. Nat Genet. 2013, 45: 656-663.
Chewapreecha C, Harris S, Croucher N, Turner C, Marttinen P, Cheng L, Pessia A, Aanensen D, Mather A, Page A, Salter S, Harris D, Nosten F, Goldblatt D, Corander J, Parkhill J, Turner P, Bentley S: Dense genomic sampling identifies highways of pneumococcal recombination. Nat Genet. 2014, 46: 305-309.
McNally A, Dalton T, Ragione RML, Stapleton K, Manning G, Newell DG: Yersinia enterocolitica isolates of differing biotypes from humans and animals are adherent, invasive and persist in macrophages, but differ in cytokine secretion profiles in vitro. J Med Microbiol. 2006, 55 (12): 1725-1734.
Alhashash F, Weston V, Diggle M, McNally A: Multidrug-Resistant Escherichia coli Bacteremia. Emerg Infect Dis. 2013, 19: 1699-1701.
McNally A, Alhashash F, Collins M, Alqasim A, Paszckiewicz K, Weston V, Diggle M: Genomic analysis of Extra-intestinal pathogenic Escherichia coli urosepsis. Clin Microbiol Infect. 2013, 19: 328-334.
Swain MT, Tsai IJ, Assefa SA, Newbold C, Berriman M, Otto TD: A post-assembly genome-improvement toolkit (PAGIT) to obtain annotated genomes from contigs. Nat Protoc. 2012, 7: 1260-1284.
Seemann T: Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014, doi:10.1093/bioinformatics/btu153
Angiuoli SVSS: Mugsy: fast multiple alignment of closely related whole genomes. Bioinformatics. 2011, 27 (3): 334-342.
Stamatakis A, Ludwig T, Meier H: RAxML-III: a fast program for maximum likelihood-based inference of large phylogenetic trees. Bioinformatics. 2005, 21 (4): 456-
Croxall G, Weston V, Joseph S, Manning G, Cheetham P, McNally A: Increased Human Pathogenic Potential of Escherichia coli from Polymicrobial Urinary Tract Infections in Comparison to Isolates from Monomicrobial Culture Samples. J Med Microbiol. 2011, 60: 102-109.
This project was funded by a personal studentship award to AA by King Saud University, and by Royal Society/NSFC international collaboration award IE121459 to AM and ZZ.
There are no ethical considerations relevant to this study, and the authors declare that they have no competing interests.
AA performed all phenotypic lab work. AA and AM analysed genomic data. AM and AZ supplied strains and reagents. FS performed classical capsule typing. AA, AM, and ZZ wrote the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Alqasim, A., Scheutz, F., Zong, Z. et al. Comparative genome analysis identifies few traits unique to the Escherichia coli ST131 H30Rx clade and extensive mosaicism at the capsule locus. BMC Genomics 15, 830 (2014). https://doi.org/10.1186/1471-2164-15-830
- E. coli ST131
- Comparative genomics