- Research article
- Open Access
Phylogenetic and genomic diversity of human bacteremic Escherichia coli strains
BMC Genomicsvolume 9, Article number: 560 (2008)
Extraintestinal pathogenic Escherichia coli (ExPEC) strains represent a huge public health burden. Knowledge of their clonal diversity and of the association of clones with genomic content and clinical features is a prerequisite to recognize strains with a high invasive potential. In order to provide an unbiased view of the diversity of E. coli strains responsible for bacteremia, we studied 161 consecutive isolates from patients with positive blood culture obtained during one year in two French university hospitals. We collected precise clinical information, multilocus sequence typing (MLST) data and virulence gene content for all isolates. A subset representative of the clonal diversity was subjected to comparative genomic hybridization (CGH) using 2,324 amplicons from the flexible gene pool of E. coli.
Recombination-insensitive phylogenetic analysis of MLST data in combination with the ECOR collection revealed that bacteremic E. coli isolates were highly diverse and distributed into five major lineages, corresponding to the classical E. coli phylogroups (A+B1, B2, D and E) and group F, which comprises strains previously assigned to D. Compared to other strains of phylogenetic group B2, strains belonging to MLST-derived clonal complexes (CCs) CC1 and CC4 were associated (P < 0.05) with a urinary origin. In contrast, no CC appeared associated with severe sepsis or unfavorable outcome of the bacteremia. CGH analysis revealed genomic characteristics of the distinct CCs and identified genomic regions associated with CC1 and/or CC4.
Our results demonstrate that human bacteremia strains distribute over the entire span of E. coli phylogenetic diversity and that CCs represent important phylogenetic units for pathogenesis and comparative genomics.
Escherichia coli is the most abundant facultative anaerobic bacteria of the human intestinal flora. Whereas E. coli usually appears to be a harmless commensal, in other circumstances, E. coli strains can be pathogenic to humans and were grouped into various pathotypes . Among these, extraintestinal pathogenic Escherichia coli (ExPEC) are responsible for urinary tract, intra-abdominal and soft tissue infections, meningitis, pneumonia and osteomyelitis often associated to bacteremia . Bacteremia represents the tenth major cause of death in developed countries and among Gram-negative bacteria, Escherichia coli represents the first cause of bacteremia, with 30% of the total number of bacteremias being due to this pathogen .
Among bacteremic isolates, those that are characterized by specific virulence factors (VFs) such as adhesins, capsule, cytotoxins and siderophores are considered as extraintestinal pathogenic E. coli (ExPEC) , as these VFs are classically described as being necessary to overcome host defenses, invade host tissues and trigger a local inflammatory response . Phylogenetic analyses [5, 6] suggested that E. coli can be divided into four major groups (A, B1, B2 and D). Generally, human ExPEC strains belong to group B2, and to a lesser extent, group D, whereas commensal strains and less virulent strains belong to A or B1 groups [2, 7–9].
Whether E. coli virulence factors influence the severity of sepsis, in humans, is not clearly established. In fact, no association between the severity of sepsis and the four phylogenetic groups or other bacterial determinants was evidenced . The difficulty to establish a link between phylogenetic groups on the one hand, and severity of infection in humans on the other hand, can have several causes. In particular, clonal structure within phylogenetic groups and heterogeneity among clonal groups in terms of severity could constitute an important confounding factor. Therefore, it appears important to characterize the clonal structure within phylogenetic groups and its association with clinical features.
Multilocus Sequence Typing (MLST), in which internal portions of multiple housekeeping genes are sequenced to define clonal diversity, has emerged as a powerful tool to describe the genetic structure of bacterial populations . In addition, MLST data allow determining the phylogenetic relationships among deep lineages, providing a complementary view of the population structure .
Several sequence-based studies were used to characterize clones or phylogenetic subgroups within E. coli [13–21]. Previous analyses of ExPEC were carried out on selected strains [21–25]. However, the identity and distribution of clonal groups in unselected bacteremia episodes have not been investigated systematically, and thus, no unbiased view of the overall diversity and distribution of human bacteremic strains is currently available.
The aims of this work were (i) To define the clonal diversity of all E. coli strains isolated from consecutive bacteremia during one year in two major university hospitals in Paris, (ii) To establish the possible association between clonal groups and clinical determinants, and (iii) To characterize the genomic content of the disclosed clonal groups.
Nucleotide polymorphism and phylogenetic diversity
The internal portions of eight selected housekeeping genes were sequenced in the 161 isolates, with only two exceptions: three isolates (AVC062, NCK061 and NCK062) were PCR-negative for trpA, and three isolates (AVC025, AVC041, NCK011) were PCR-negative for uidA. All sequences were aligned without insertion or deletion, with the exception of two isolates (NCK018 and NCK043), which showed a 5-nucleotide deletion in uidA, resulting in a frameshift. The proportion of variable sites ranged from 9.3% (uidA) to 15.5% (trpA) (Table 1). Considering the 4,095 nucleotides of the eight gene portions together, there were 471 (11.5%) variable sites. The nucleotide diversity (π, average number of nucleotide differences per site between two randomly-selected isolates) ranged from 0.0174 (putP) to 0.0393 (trpA) (Table 1).
E. coli strains are generally considered to be grouped into the four major phylogenetic groups A, B1, B2 and D [5, 6], although some strains may belong to additional groups [19, 26]. We characterized all isolates by the triplex PCR method, designed to typify strains into four categories equated to the four major groups . Based on this PCR scheme, the number of bacteremic isolates classified into groups A, B1, B2 and D were 33 (20%), 11 (7%), 80 (50%) and 37 (23%), respectively. To determine precisely the phylogenetic diversity of the 161 bacteremic isolates, we sequenced the eight gene portions in strains of the ECOR collection and performed a joint phylogenetic analysis, together with seven genome reference strains (Figure 1). The results highlighted the strong phylogenetic clustering of E. coli strains into five sharply separated branches. Based on triplex PCR results and ECOR strains, these branches could be equated to groups B2, A+B1, D, E and a group, which we refer to as group F, comprising strains assigned to group D based on triplex PCR. Of note, when genetic distances among strains were considered, the B2 branch was more heterogeneous than the branch comprising strains of groups A and B1 together. Group A appeared paraphyletic, as B1 strains were nested within the diversity of A strains. A few confirmed inconsistencies (for isolates NCK23, NCK45 and AVC003, Figure 2) were found between triplex-PCR and sequence-based phylogenetic relationships, as also reported by others [20, 26]. Remarkably, the bacteremic isolates were distributed into the five major branches, clearly demonstrating that they do not represent a restricted subset of E. coli strains.
Genotype diversity and identification of clonal complexes
The 161 clinical isolates were distributed into 87 sequence types (ST), seven of which comprised five or more strains, representing 38% of the strains in total (Additional file 1). The most frequent ST encountered was ST1 (ST78 or ST79 in Mark Achtman's MLST scheme; please see Additional file 2 for known correspondence of STs in our scheme and Mark Achtman's one), with 15 isolates. Sixty-nine STs corresponded to a single isolate, and Simpson's diversity index was 97.6%, underlying the diversity of bacteremic E. coli isolates revealed by MLST. Of note, frequent STs mostly belonged to groups A and B2, resulting in a lower Simpson's index of diversity of groups A (87.7%) and B2 (92.9%), compared in particular to group B1, were all sampled isolates were distinct (Simpson's index, 100%). There was no clear difference in genotype diversity between the two hospital centers, and frequent STs were recovered from each of them.
Clonal families within bacterial populations in which homologous recombination occurs, as is the case for E. coli [19, 28], are best identified based on allelic profile comparisons, as the collapsing of nucleotide polymorphisms into allelic numbers is less sensitive to the disturbing impact of recombination . As determined using eBURST and a minimum spanning tree (MStree) analysis (Figure 2), the 87 STs were grouped into 19 clonal complexes (CC), i.e. groups of closely related STs differing by no more than one allele from another member of the group. Eleven major CCs (CC1 to CC4, CC6, CC10, CC32, CC36, CC40, CC43 and CC66) comprised five or more isolates. Based on available complete genome sequences and ECOR strains, correspondence of CCs with a previous classification  was established (Figure 2). Interestingly, no member of subgroup VIII (ST149)  was found among the bacteremic clinical isolates. Strains of this clone, represented by the ED1a genome strain (Figure 2), have been reported as human-specific and strictly commensal . Each individual CC probably has a single ancestor for most parts of its genome, and CCs can therefore be equated to clones, as it is very likely that strains sharing seven out of eight alleles have a common descent. Although large allelic profile distances separated most CCs, the links among clones disclosed by MStree analysis (Figure 2A) were consistent with phylogenetic group assignments. For example, CCs of group B2 generally shared more alleles in common than with STs of other groups.
Association of clonal complexes with clinical determinants
In a previous work, differences were found among phylogenetic groups with respect to the primary source of bacteremia . B2 was significantly associated with urosepsis and immunocompetent hosts, whereas non-B2 isolates were associated with non-urinary tract origins and immunocompromised hosts. Here, we tested the hypothesis of a possible association of particular clones (the above described 11 major clones), even within a given phylogenetic group, with clinical determinants (Additional file 1).
Based on a logistic regression model analysis, we did not find statistical evidence for an association of CCs with a severe sepsis or a defavorable outcome. However, the two major clones of group B2 (CC1 and CC4) were significantly associated with urinary tract as a source of bacteremia (Table 2), in comparison to other B2 genotypes. There were 12 (70%) urosepsis cases out of 17 isolates of CC1 (p = 0.017), and 11 (65%) urosepsis cases out of 17 isolates of CC4 (p = 0.041). In contrast, in the 46 remaining B2 isolates, there were only 14 (30%) urosepsis cases. Conversely, B2 clone CC32 was significantly associated with a non urinary tract origin as primary source of infection (p = 0.023). The other clones of group B2 showed mixed origins with few urosepsis cases (Table 2). A number of other characteristics of particular clones compared to other B2 isolates were statistically supported: CC36 was associated with female gender (p = 0.005), although it was not associated with urinary origin; CC32 was found only in males (p = 0.029), CC40 was associated with community-acquired infection (p = 0.048) and CC6 tended to be associated with diabetes mellitus (p = 0.059). These results may indicate heterogeneity of biological characteristics among CCs of phylogenetic group B2. Similarly, among phylogenetic group D isolates, CC3 was significantly associated with a digestive tract origin of the infection (p = 0.038). Among group A isolates, isolates of CC66 were significantly associated with renal failure (p = 0.049), and CC2 isolates were weakly associated with neoplasia (p = 0.069).
Virulence factor content of clonal complexes
Virulence factor (VF) distribution is known to differ among the four phylogenetic groups [5–7]. In our previous study, isolates of group B2 were characterized by a high number of VFs, whereas those belonging to groups A and D exhibited lower numbers . Here, we compared at the clone level the presence of nine genes associated with ExPEC virulence, and of the svg ORF, recently described as specific for a highly virulent B2 subgroup  (Additional file 1). As expected, we found a strong association between CCs and VF content, particularly among B2 isolates (Table 3). For example, clones CC1 and CC4 differed significantly by the number of isolates containing hlyC, cnf and sfa (low in CC1 but high in CC4). In addition, iucC was significantly associated with CC1 and CC4 isolates, compared to the remaining B2 isolates. The absence of the adhesin papGII was statistically associated with B2-CC32, B2-CC10 and B2-CC40, compared to the other B2 isolates. In addition, B2-CC43 and B2-CC6 were almost devoid of the investigated VFs, with statistical support for papC, sfa, hlyC, cnf1 and iucC (p < 0.05). Finally, all isolates of ST1 (the central genotype of CC1) except one were positive upon svg PCR, consistent with previous results . Only one non-ST1 isolate (ST41, member of CC40) was also positive for svg.
Flexible gene pool content of bacteremic isolates
In order to determine the genomic characteristics of clones and to identify possible associations between genetic factors and clinical features, we designed a macroarray composed of 2,324 probes and analyzed the genome content of a subset of 60 bacteremic isolates selected to be representative of the diversity of genotypes (Additional file 1). Generally, the DNA array hybridization results were in accordance with housekeeping genes-based phylogeny and genotypic (ST) classification. Cluster analysis of hybridization data grouped all B2 isolates in a single major group and other phylogroups were also recovered with few exceptions (Figure 2B). Interestingly, subgroups recovered by cluster analysis of DNA array data corresponded to clonal complexes with a single exception (NCK065-CC1, Figure 2B), showing that CCs are characterized by a specific genomic content. In addition, DNA array data were fully consistent with the above-described virulence factor PCR results.
Because the strongest link of clonal diversity with clinical features within major phylogenetic groups was the association of two B2 clones (CC1 and CC4) with urinary tract origin, we concentrated our analysis on the genomic particularities of CC1 and CC4. We identified 145 (6%) probes that were distributed differently in CC1 and/or CC4 when compared to other B2 isolates (Figure 3, Additional file 3), as determined by the template matching clustering approach .
Among probes that reacted positively with both CC1 and CC4 but not with other B2 isolates, were two clusters of six and eight sequences, respectively (c1445-c1453 and c1569-c1577) (Fisher exact test, p < 0.05) (Figure 3). These prophage sequences are included in previously described genomics islands [32, 33] in the CFT073 chromosome, respectively φ-CFT073-potB and PAI-CFT073-icdA. Consistently, these sequences were also present in the O45:K1 strain S88 (ST1, CC1), a meningitis-causing strain that is closely related  to E. coli K1 RS218 strain (ST1, CC1). The two groups of phage-associated sequences were not restricted to CC1-CC4 B2 isolates (present in strain 536 and B2-CC32, Figure 3) and also detected in some isolates of groups A, B1 and D.
Genomic regions that were specific for CC4 were also identified. Isolates of CC4 were characterized, compared to other B2 isolates, by the presence of three operons composed respectively of three (c1884-c1886), six (c2893-c2898) and seven genes (c5372-c5381) (p < 0.001, Fisher exact test), encoding putative proteins of unknown function. Consistently, these genes were present in the genome of CFT073 (ST4, CC4), but not in the genomes of strains 536 (ST261, CC32) and S88 (ST1, CC1). Moreover, these operons were detected in none of the non-B2 isolates. In addition, a large DNA region from CFT073 (c0358-c0365) was detected more significantly in CC4 isolates than in other B2 isolates (p = 0.01 Fisher exact test, present also in B2-CC36), and all isolates of groups A, B and D were negative for this DNA region (Figure 3 and Additional file 3). This region comprised three genes (c0361-c0363) encoding two putative membrane export proteins and a putative RTX exotoxin A. Regarding CC1-associated genes, only two bacteriophage proteins described in RS218 (ST1, CC1) genome, one Q homolog phage protein and a K-12 MG1655 hypothetical protein (ydfC), were found. Conversely, two groups of sequences were always absent from CC1 and CC4 but present in most other B2 isolates (Figure 3; p < 0.05, Fisher exact test).
We analyzed the phylogenetic diversity of an unselected set of 161 consecutive bacteremic E. coli isolates, and compared the genomic content of a representative subset of genotypes in order to investigate the link between clonal diversity, gene content and clinical correlates of bacteremia.
Bacteremic isolates were distributed across the entire span of E. coli phylogenetic diversity. The amount of polymorphism found herein is consistent with values reported in E. coli, based on other housekeeping genes [18, 19, 35, 36]. Amounts of polymorphism among bacteremia isolates were only slightly lower than a recent study of 185 isolates from freshwater beaches . For example, the authors found 49 uidA alleles with 12% polymorphic sites, whereas we found 27 uidA alleles and 9.3% polymorphic sites (note that our collection of bacteremic isolates included only 7% of isolates of group B1, whereas this group represented 56% of isolates from the environment ). Consequently, our results clearly indicate that human E. coli bacteremia strains are genetically highly diverse. Of note, avian pathogenic E. coli strains (APEC) overlap only partly with the human ExPEC population, with some potentially zoonotic clonal groups being found both in humans and birds, while other APEC clonal groups are rare or absent among human ExPEC [37–39]. Thus, total ExPEC diversity may exceed the large diversity of human bacteremic isolates alone.
Although it was not the main purpose of this study, we noted that our combined analysis of bacteremic isolates and of ECOR reference strains using ClonalFrame revealed the existence of sharp demarcations among deep branches, suggesting a much stronger structure within the E. coli species than previously disclosed based on recombination-sensitive phylogenetic analyses, such as the neighbor-joining method. So far, methods that detect homologous recombination events and account for these in phylogenetic reconstruction have not been used widely in E. coli, and it is likely that phylogenetic group demarcation may have been problematic until now given the effects of homologous recombination, which tends to homogenize E. coli diversity and blur the neat borders that may delineate independently evolving phylogenetic lineages [19, 26]. In addition, classification into only four major groups appears to be an oversimplification of a more complex reality. The five major branches disclosed herein do not correspond totally to previous group definitions, as groups A and B1 were not well separated, and as branch F could not be equated to any previously described group. These results challenge the classical view of the internal phylogenetic structure of E. coli, in agreement with recent studies [18, 19].
Our initial hypothesis was that clones, rather than entire phylogenetic groups, may be more relevant natural entities to establish an association of genotype with phenotype, including clinical correlates of bacteremic isolates. In order to define groups of closely related genotypes, we used allelic profile-based comparisons, rather than nucleotide sequence-based analysis [16, 21, 23]. Because the later approach is sensitive to homologous recombination, analysis of allelic profiles is preferred for clone definition in many bacterial species [12, 40]. Among E. coli bacteremic isolates, a number of CCs, which can be equated to clones, could thus be defined. Interestingly, the central genotype of CCs, inferred to represent the founder of the clone, generally had a higher frequency than its variants, consistent with currently ongoing or recent expansion of the clones . For the purposes of classification, it is noteworthy that most CCs were separated by large allelic profile distances (four or five mismatches), highlighting the neat demarcation among them. Phylogenetic group B2 appeared to be the most strongly structured group, as most isolates were grouped into nine CCs (Figure 2).
We explored the hypothesis that clinical correlates and distribution of particular genes would be associated with specific clonal backgrounds, rather than being distributed widely across E. coli diversity. Because horizontal gene transfer occurs among E. coli strains, virulence genes that confer specific advantages e.g. for urinary tract infection can be distributed in various genomic backgrounds. However, association of specific genes with particular genomic backgrounds (clonal groups) can be retained, at least in the short term, by predominantly vertical transmission e.g. for genes that are not harbored on mobile elements. In this case, knowledge of the clonal background would have predictive value regarding gene content and corresponding phenotypes .
Consistent with our initial hypothesis that particular clones may exhibit clinically relevant features, several correlations were established between some CCs and clinical data. Most notably, CC1 and CC4 were clearly associated with isolates responsible for urosepsis; these CCs correspond to previously-described subgroups II (strain CFT073) and IX (strain RS218), corresponding to strains responsible of urosepsis and of neonatal meningitis, respectively . Hence, the previously held association of B2 as a whole with urosepsis may be an oversimplification due to the strong contribution of CC1 and CC4 to isolates of this group, and may thus not be valid for all B2 strains. Interestingly, we detected by DNA array analysis, two specific gene clusters associated with CC1 and CC4. These sequences, which correspond to cryptic prophages, were previously described as being located in two consecutive genomic islands . Consistent with our results, one of these clusters, PAI-CFT073-icdA, was described as being more frequent among pyelonephritis isolates . Likewise, a specific DNA region encoding to a putative RTX protein was significantly associated to CC4 isolates. A first RTX-toxin described in E. coli, hemolysin A, has been clearly associated with pyelonephritis and implicated in inflammation during urinary tract infection [42, 43]. Some UPEC strains including CFT073 harbor two operons of this toxin, and the specific role of each copy remains unclear [33, 44]. Whether the new putative RTX toxin constitutes an advantage to CC4 isolates and could explain their specific urovirulence remains to be determined.
As assessed by PCR and DNA array hybridization, clonal complexes were characterized by specific gene content. For example, isolates belonging to CC4 exhibited significantly more VFs (sfa, cnf1, hlyC) than CC1 isolates, and ORF svg was associated with CC1. In addition, the pattern of gene content variation (Figure 3) was highly concordant with clonal complexes, clearly illustrating that genomic background as assessed by MLST and gene content are strongly correlated. This pattern is consistent with the well-established mechanisms of deletion or acquisition of entire PAIs [45–50]. In order to get clues into the possible biological and clinical significance of gene content differences among clones, it will be necessary to combine functional studies with the analysis of E. coli isolates from other sources. For example, comparing the gene content of bacteremic and commensal isolates within a single clone should provide insights into microevolutionary events leading to increased or decreased pathogenic potential.
Our results confirmed that most of the usually recognized extra-intestinal VFs (e.g. pap, sfa, hly) were concentrated in the virulent CCs, particularly those belonging to phylogenetic group B2. In contrast, others VFs, for example those related to the plasmid-encoded iron uptake system (e.g. iucC, iroN), were more broadly distributed [45, 51]. Therefore, for highly mobile genetic elements (e.g. plasmids), the association with clonal background may not be retained.
In conclusion, we characterized the clonal diversity of all consecutive isolates responsible for bacteremia during one year in two hospitals, revealing for the first time their distribution into phylogenetic groups and clonal groups in an unbiased way. We demonstrated the existence of differences in gene content or clinical parameters among some of the major clonal families identified. Although they should be confirmed on larger populations to be firmly established, these results confirm that particular STs, CCs or phylogenetic subgroups, rather than the more inclusive higher categories known as major phylogenetic groups, will be relevant units of pathogenicity. However, the heterogeneity of gene content within clones, although limited, may be highly significant clinically, which in turn suggests that finer phylogenetic subdivisions within CCs may also differ in clinically important characteristics. The failure to associate CCs with severity of clinical outcome may also indicate that the characteristics of the host play a dominant role.
We studied 161 well-characterized E. coli isolates from bacteremia; this set of isolates (Additional file 1) has been previously described . Briefly, isolates were collected between December 2002 and December 2003 and correspond to all consecutive episodes of E. coli bacteremia in two major university hospitals in Paris. Epidemiological (age, sex gender), clinical (community/nosocomial acquired infection, immune status, underlying diseases, primary source of infection, severity sepsis scoring and outcome), and laboratory data for each episode were recorded in an anonymous computer database in accordance with French law. We determined the primary source of bacteremia by clinical and radiological presentation and/or by evidence, based on antibiogram and/or serotyping, of an identical strain isolated from blood and other body sites culture . When the primary source of infection was not found, the origin of infection was assigned to the digestive tract. All the above 161 isolates are listed in supplementary Table 1. The ECOR collection of 72 reference strains  was included for phylogenetic comparisons; five ECOR strains could not be sequenced at all eight genes (see below) and were removed from the analysis: ECOR5 (A), ECOR6 (A), ECOR9 (A), ECOR13 (A) and ECOR64 (B2). In addition, seven strains for which the complete genome sequence is available (K-12 MG1655, O157:H7-EDL933, 536, CFT073, UTI89, RS218) or underway (ED1a, ColiScope project, http://www.genoscope.cns.fr/spip/Escherichia-fergusonii.html) were used.
Triplex-PCR for tentative group assignment
A triplex-PCR method  was used to tentatively assign the isolates to the four major E. coli phylogenetic groups A, B1, B2 and D.
Multilocus Sequence Typing
Primer pairs were designed for PCR amplification and sequencing of internal portions of eight housekeeping genes (Table 4). Selected genes included dinB, icdA, pabB, polB, putP, trpA and trpB, previously used for phylogenetic analysis of E. coli/Shigella strains [53–55]. New PCR primers were designed in internal portions of the genes based on previously obtained sequences, in order to amplify target regions of approximately 500 – 600 bp (Table 4). These seven genes represent six distinct loci on the E. coli chromosome, as trpA and trpB are located in the same operon. To increase the number of loci to seven, we added gene uidA, which is used for E. coli MLST by the group of Tom Whittam . Universal sequencing primer sequences were added to the 5' end of the PCR primers. All PCR products were thus sequenced using the same two sequencing primers (Table 4). Further details on this MLST scheme can be found at http://www.pasteur.fr/mlst.
Nucleotide sequences were obtained using Big Dye version 3.1 chemistry on ABI 3100 or 3730 apparatuses. In order to eliminate the risk of sample mix-up, PCR and sequencing were performed using a molecular biology robot (RoboAmp 4200-PE; MWG Biotech, Courtaboeuf, France). Sequence chromatograms were edited and stored using BioNumerics version 4.5 (Applied-Maths, St. Maartens-Latem, Belgium). All nucleotides within the consensus sequence template were supported by at least two sequence chromatograms. A different allele number was given to each distinct sequence within a locus, and a distinct sequence type (ST) number was attributed to each distinct combination of alleles. Null alleles corresponding to negative PCR amplification were considered as alignment gaps in phylogenetic analyses and as allele '999' in profile-based analyses. Isolates were grouped into clonal complexes (CCs) by eBURST, if they differed at no more than 1 locus from at least one other member of the group . Founder genotypes of CCs were defined as the ST of the CC with the highest number of neighboring STs (single locus variants). Nucleotide diversity was calculated using DNAsp version 4 . Minimum spanning tree analysis was performed using BioNumerics version 5.10. MEGA  was used to draw the consensus phylogenetic tree obtained using ClonalFrame  after 100,000 iterations, including 50,000 burn-in.
Allele sequences and STs are available on Institut Pasteur's MLST web site at http://www.pasteur.fr/mlst.
A DNA macroarray was adapted from a previously described array . This E. coli pathoarray was developed with the aim to contain a large E. coli flexible gene pool. DNA probes were selected using the genome sequence (GenBank accession CP000247) of the O6:K15:H31 uropathogenic E. coli strain 536 , the sequence (GenBank accession AE014075) of the O6:K2:H1 uropathogenic E. coli strain CFT073 , and E. coli K-12 (MG1655, GenBank accession U00096) . DNA sequences showing an homology above 95% within the three genomes were excluded. In addition, specific sequences of the meningitis-associated strain RS218 , i.e. sequences with less that 95% homology in the above three strains, were included. The genome sequence of strain RS218 was only partially available at http://www.genome.wisc.edu at the time of membrane conception.
Amplicons were designed to correspond to annotated ORFs. Additional file 4 describes the composition of the membrane. When the size of an ORF was less than 800 bp the entire sequence corresponding to the ORF was selected. ORF above 1 kb were represented on the membrane by at least 2 amplicons. In addition sequences of intergenic regions which were not present in all three genomes were amplified and represented on the membrane if above 500 bp. 2,324 specific genomic DNA fragments corresponding to 2,196 ORFs and 91 intergenic regions (IRs) were spotted on the membrane. Fourteen additional sequences were added in the DNA array. These sequences corresponded to specific virulence-associated ExPEC genes characterized in recent studies but not represented after comparative genomic analysis of the above strains.
The amplicons were spotted in duplicate by a robot (Eurogentec, Seraing, Belgium) on nylon membranes and fixed by alkali treatment. Positive (186 spots containing 16S rDNA and genomic DNA of CFT073 strain) and negative (15 spots containing both human and mouse gene DNA and 14 spots with absence of DNA) control spots were used for normalization.
Genomic DNA extraction, labeling and macroarray hybridization
Bacteria were grown overnight at 37°C in aerobic condition on agar plate. Genomic DNA was isolated using the DNA Wizard Genomic DNA Purification Kit (Promega, Madison, USA). Genomic DNA extract (100 ng to 500 ng) were labelled with 50 mCi of [α33P] dCTP (Amersham Biosciences, Orsay, France), using the Rediprime II random prime labelling system (GE Healthcare, Amersham, UK) according to manufacturer's instructions. The membranes were soaked in 6× SSC and pre-hybridized in Shake 'n' Stack hybridization over (Hybaid, Thermo Scientific, Canada) for 4 hours at 65°C with 6 ml Church & Gilbert hybridization buffer (0.5 M NaPi, 1 mM EDTA, 7% SDS) containing 100 μg/ml heat-denatured salmon sperm DNA (Invitrogen, Life Technologies, France). Labeled genomic DNA was denatured at 100°C for 5 min and hybridization was performed during 15–18 h at 65°C in the Church & Gilbert buffer. After hybridization, each nylon membrane was washed three times in the same solution (40 mM NaPi, 1 mM EDTA, 1% SDS) for 45 min at 65°C. The membranes were then exposed to a PhosphorImager screen for 48 h. The PhosphorImager screens were scanned on a Storm 860 PhosphorImager device (GE Healthcare, Life Science, UK) at a pixel size of 50 μm. Before being re-hybridized, filters were stripped in 125 ml of a bowled buffer (10 mM Tris-HCl pH 7.6, 1 mM EDTA pH 8.0, 1% SDS) for 25 min at 100°C. A unique hybridization experiment was performed for each isolate investigated in this study and at least two independently experiments were performed for E. coli control strains.
Macroarray quality control
This macroarray was first tested using the three sequenced strains, i.e. E. coli strain 536, strain CFT073 and E. coli K-12. This experiment was performed twice for strain 536, four times for strain CFT073 and three times with E. coli K-12 MG1655. False negative results were obtained in less than 0.01%. The average number of false positive was 8%, 9% and 12.5% using strain 536, CFT073 and K-12 MG1655 E. coli, respectively. Thirty-five amplicons, giving repeatedly false positive or negative results, were removed from the final analysis and are labeled as such in Additional file 1. The remaining sequences on the array were annotated as coding for (i) cell structure membrane proteins (127), (ii) putative functional enzymes involved in metabolic pathways (493), (iii) putative transcriptional regulators (63), (iv) known virulence factors (204), (v) hypothetical proteins (880), and (vi) mobile genetic elements such as transposases, IS or phages (219). To get insight into the specific function of RS218-derived sequences, we used the complete genome sequence (ColiScope project) of the phylogenetically related strain S88 .
Macroarray data submission
Macroarray data were submitted to the ArrayExpress database under accession number: A-MEXP-1451.
Macroarray data analysis
The macroarray data were analysed using the ArrayVision software (Imaging Reseach, St Catharines, Canada) for signal quantification. To avoid empirical determination of cut-offs, we developed R scripts (based on the Package Mclust for Normal Mixture modeling) adapted from a program developed for comparative genomic hybridization analysis of Enterococcus faecalis . This analysis was based on bimodal distribution of hybridization signals, thus highlighting the existence of two different gene populations corresponding to high or low intensities. Each population of gene fits a Gaussian distribution model and allows the use of an algorithm for gene classification based on the probability that each value of probe belonged to the population of high or low hybridization signals (i.e. present, class 1 or absent genes, class -1). Parameters of the model are estimated by maximum likelihood using an expectation maximization (EM) algorithm running from the median of negative controls to the median of positive controls. Using this mathematic model, genes were classified to be either present (probability of belonging to class 1 > 80%) or absent (probability of belonging to class 1 < 20%). Some genes were considered ambiguous (probability of belonging to class 1 between 20% and 80%). R analysis was done for all the hybridized arrays and a matrix of 1 (gene present in this strain), 0 (ambiguous gene) and -1 (gene absent is this strain) was obtained. TMeV software was used with this matrix for data clustering and template matching analysis .
Comparisons were based on Chi-square test for categorical variables and Fisher exact test when numbers were below five. All tests were two-tailed and p < 0.05 was considered significant. Multiple logistic regressions were used to assess the relationships between two clonal complexes, the binary variables, and explanatory variables: host characteristics and bacterial determinants. Variables with a p value less than 0.25 in univariate analyses were included into the model and interaction between significant variable was estimated. Only significant variables in the model were conserved to determine the statistically significant relationships between CC and clinical determinants. STATA version 8 was the statistical software used for Chi-square and multiple logistic regression models.
Multilocus sequence typing
Escherichia coli reference collection
Extraintestinal pathogenic Escherichia coli
comparative genomic hybridization
minimum spanning tree.
Dobrindt U: (Patho-)Genomics of Escherichia coli. Int J Med Microbiol. 2005, 295: 357-371. 10.1016/j.ijmm.2005.07.009.
Russo TA, Johnson JR: Proposal for a new inclusive designation for extraintestinal pathogenic isolates of Escherichia coli: ExPEC. J Infect Dis. 2000, 181: 1753-1754. 10.1086/315418.
Martin GS, Mannino DM, Eaton S, Moss M: The epidemiology of sepsis in the United States from 1979 through 2000. N Engl J Med. 2003, 348: 1546-1554. 10.1056/NEJMoa022139.
Johnson JR, Gajewski A, Lesse AJ, Russo TA: Extraintestinal pathogenic Escherichia coli as a cause of invasive nonurinary infections. J Clin Microbiol. 2003, 41: 5798-5802. 10.1128/JCM.41.12.5798-5802.2003.
Selander RK, Caugant DA, Whittam TS: Genetic structure and variation in natural populations of Escherichia coli. Escherichia coli and Salmonella typhimurium: cellular and molecular biology. Edited by: Neidhardt FC, Ingraham JL, Low KB, Magasanik B, Schaechter M, Umbarger HE. 1987, Washington D.C.: ASM, 1625-1648.
Herzer PJ, Inouye S, Inouye M, Whittam TS: Phylogenetic distribution of branched RNA-linked multicopy single-stranded DNA among natural isolates of Escherichia coli. J Bacteriol. 1990, 172: 6175-6181.
Picard B, Garcia JS, Gouriou S, Duriez P, Brahimi N, Bingen E, Elion J, Denamur E: The link between phylogeny and virulence in Escherichia coli extraintestinal infection. Infect Immun. 1999, 67: 546-553.
Johnson JR, Stell AL: Extended virulence genotypes of Escherichia coli strains from patients with urosepsis in relation to phylogeny and host compromise. J Infect Dis. 2000, 181: 261-272. 10.1086/315217.
Johnson JR, Russo TA: Extraintestinal pathogenic Escherichia coli: "the other bad E coli". J Lab Clin Med. 2002, 139: 155-162. 10.1067/mlc.2002.121550.
Martinez JA, Soto S, Fabrega A, Almela M, Mensa J, Soriano A, Marco F, Jimenez de Anta MT, Vila J: Relationship of phylogenetic background, biofilm production, and time to detection of growth in blood culture vials with clinical variables and prognosis associated with Escherichia coli bacteremia. J Clin Microbiol. 2006, 44: 1468-1474. 10.1128/JCM.44.4.1468-1474.2006.
Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, Zhang Q, Zhou J, Zurth K, Caugant DA: Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci USA. 1998, 95: 3140-3145. 10.1073/pnas.95.6.3140.
Feil EJ: Small change: keeping pace with microevolution. Nat Rev Microbiol. 2004, 2: 483-495. 10.1038/nrmicro904.
Goldberg TL, Gillespie TR, Singer RS: Optimization of analytical parameters for inferring relationships among Escherichia coli isolates from repetitive-element PCR by maximizing correspondence with multilocus sequence typing data. Appl Environ Microbiol. 2006, 72: 6049-6052. 10.1128/AEM.00355-06.
Nemoy LL, Kotetishvili M, Tigno J, Keefer-Norris A, Harris AD, Perencevich EN, Johnson JA, Torpey D, Sulakvelidze A, Morris JG, Stine OC: Multilocus sequence typing versus pulsed-field gel electrophoresis for characterization of extended-spectrum beta-lactamase-producing Escherichia coli isolates. J Clin Microbiol. 2005, 43: 1776-1781. 10.1128/JCM.43.4.1776-1781.2005.
Tartof SY, Solberg OD, Manges AR, Riley LW: Analysis of a uropathogenic Escherichia coli clonal group by multilocus sequence typing. J Clin Microbiol. 2005, 43: 5860-5864. 10.1128/JCM.43.12.5860-5864.2005.
Johnson JR, Owens KL, Clabots CR, Weissman SJ, Cannon SB: Phylogenetic relationships among clonal groups of extraintestinal pathogenic Escherichia coli as assessed by multi-locus sequence analysis. Microbes Infect. 2006, 8: 1702-1713. 10.1016/j.micinf.2006.02.007.
Bidet P, Metais A, Mahjoub-Messai F, Durand L, Dehem M, Aujard Y, Bingen E, Nassif X, Bonacorsi S: Detection and identification by PCR of a highly virulent phylogenetic subgroup among extraintestinal pathogenic Escherichia coli B2 strains. Appl Environ Microbiol. 2007, 73: 2373-2377. 10.1128/AEM.02341-06.
Walk ST, Alm EW, Calhoun LM, Mladonicky JM, Whittam TS: Genetic diversity and population structure of Escherichia coli isolated from freshwater beaches. Environ Microbiol. 2007, 9: 2274-2288. 10.1111/j.1462-2920.2007.01341.x.
Wirth T, Falush D, Lan R, Colles F, Mensa P, Wieler LH, Karch H, Reeves PR, Maiden MC, Ochman H, Achtman M: Sex and virulence in Escherichia coli: an evolutionary perspective. Mol Microbiol. 2006, 60: 1136-1151. 10.1111/j.1365-2958.2006.05172.x.
Lacher DW, Steinsland H, Blank TE, Donnenberg MS, Whittam TS: Molecular evolution of typical enteropathogenic Escherichia coli: clonal analysis by multilocus sequence typing and virulence gene allelic profiling. J Bacteriol. 2007, 189: 342-350. 10.1128/JB.01472-06.
Le Gall T, Clermont O, Gouriou S, Picard B, Nassif X, Denamur E, Tenaillon O: Extraintestinal virulence is a coincidental by-product of commensalism in B2 phylogenetic group Escherichia coli strains. Mol Biol Evol. 2007, 11: 2373-2384. 10.1093/molbev/msm172.
Johnson JR, Delavari P, Stell AL, Whittam TS, Carlino U, Russo TA: Molecular comparison of extraintestinal Escherichia coli isolates of the same electrophoretic lineages from humans and domestic animals. J Infect Dis. 2001, 183: 154-159. 10.1086/317662.
Bidet P, Mahjoub-Messai F, Blanco J, Blanco J, Dehem M, Aujard Y, Bingen E, Bonacorsi S: Combined multilocus sequence typing and O serogrouping distinguishes Escherichia coli subtypes associated with infant urosepsis and/or meningitis. J Infect Dis. 2007, 196: 297-303. 10.1086/518897.
Johnson JR, Stell AL, O'Bryan TT, Kuskowski M, Nowicki B, Johnson C, Maslow JN, Kaul A, Kavle J, Prats G: Global molecular epidemiology of the O15:K52:H1 extraintestinal pathogenic Escherichia coli clonal group: evidence of distribution beyond Europe. J Clin Microbiol. 2002, 40: 1913-1923. 10.1128/JCM.40.6.1913-1923.2002.
Manges AR, Johnson JR, Foxman B, O'Bryan TT, Fullerton KE, Riley LW: Widespread distribution of urinary tract infections caused by a multidrug-resistant Escherichia coli clonal group. N Engl J Med. 2001, 345: 1007-1013. 10.1056/NEJMoa011265.
Gordon DM, Clermont O, Tolley H, Denamur E: Assigning Escherichia coli strains to phylogenetic groups: multi-locus sequence typing versus the PCR triplex method. Environ Microbiol. 2008, 10: 2484-2496. 10.1111/j.1462-2920.2008.01669.x.
Clermont O, Bonacorsi S, Bingen E: Rapid and simple determination of the Escherichia coli phylogenetic group. Appl Environ Microbiol. 2000, 66: 4555-4558. 10.1128/AEM.66.10.4555-4558.2000.
Dykhuizen DE, Green L: Recombination in Escherichia coli and the definition of biological species. J Bacteriol. 1991, 173: 7257-7268.
Clermont O, Lescat M, O'Brien CL, Gordon DM, Tenaillon O, Denamur E: Evidence for a human-specific Escherichia coli clone. Environ Microbiol. 2008, 10: 1000-1006. 10.1111/j.1462-2920.2007.01520.x.
Jaureguy F, Carbonnelle E, Bonacorsi S, Clec'h C, Casassus P, Bingen E, Picard B, Nassif X, Lortholary O: Host and bacterial determinants of initial severity and outcome of Escherichia coli sepsis. Clin Microbiol Infect. 2007, 13: 854-862. 10.1111/j.1469-0691.2007.01775.x.
Pavlidis P, Noble WS: Analysis of strain and regional variation in gene expression in mouse brain. Genome Biol. 2001, 2: RESEARCH0042-10.1186/gb-2001-2-10-research0042.
Lloyd AL, Rasko DA, Mobley HL: Defining genomic islands and uropathogen-specific genes in uropathogenic Escherichia coli. J Bacteriol. 2007, 189: 3532-3546. 10.1128/JB.01744-06.
Welch RA, Burland V, Plunkett G, Redford P, Roesch P, Rasko D, Buckles EL, Liou SR, Boutin A, Hackett J: Extensive mosaic structure revealed by the complete genome sequence of uropathogenic Escherichia coli. Proc Natl Acad Sci USA. 2002, 99: 17020-17024. 10.1073/pnas.252529799.
Bonacorsi S, Clermont O, Houdouin V, Cordevant C, Brahimi N, Marecat A, Tinsley C, Nassif X, Lange M, Bingen E: Molecular analysis and experimental virulence of French and North American Escherichia coli neonatal meningitis isolates: identification of a new virulent clone. J Infect Dis. 2003, 187: 1895-1906. 10.1086/375347.
Feil EJ, Holmes EC, Bessen DE, Chan MS, Day NP, Enright MC, Goldstein R, Hood DW, Kalia A, Moore CE: Recombination within natural populations of pathogenic bacteria: short-term empirical estimates and long-term phylogenetic consequences. Proc Natl Acad Sci USA. 2001, 98: 182-187. 10.1073/pnas.98.1.182.
Reid SD, Herbelin CJ, Bumbaugh AC, Selander RK, Whittam TS: Parallel evolution of virulence in pathogenic Escherichia coli. Nature. 2000, 406: 64-67. 10.1038/35017546.
Rodriguez-Siek KE, Giddings CW, Doetkott C, Johnson TJ, Fakhr MK, Nolan LK: Comparison of Escherichia coli isolates implicated in human urinary tract infection and avian colibacillosis. Microbiology. 2005, 151: 2097-2110. 10.1099/mic.0.27499-0.
Johnson TJ, Wannemuehler Y, Johnson SJ, Stell AL, Doetkott C, Johnson JR, Kim KS, Spanjaard L, Nolan LK: Comparison of Extraintestinal Pathogenic Escherichia coli from Human and Avian Sources Reveals a Mixed Subset Representing Potential Zoonotic Pathogens. Appl Environ Microbiol. 2008
Moulin-Schouleur M, Reperant M, Laurent S, Bree A, Mignon-Grasteau S, Germon P, Rasschaert D, Schouler C: Extraintestinalpathogenic Escherichia coli strains of avian and human origin: link between phylogenetic relationships and common virulence patterns. J Clin Microbiol. 2007, 45: 3366-3376. 10.1128/JCM.00037-07.
Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG: eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol. 2004, 186: 1518-1530. 10.1128/JB.186.5.1518-1530.2004.
Turner KM, Feil EJ: The secret life of the multilocus sequence type. Int J Antimicrob Agents. 2007, 29: 129-135. 10.1016/j.ijantimicag.2006.11.002.
Blanco M, Blanco JE, Alonso MP, Blanco J: Virulence factors and O groups of Escherichia coli isolates from patients with acute pyelonephritis, cystitis and asymptomatic bacteriuria. Eur J Epidemiol. 1996, 12: 191-198. 10.1007/BF00145506.
Uhlen P, Laestadius A, Jahnukainen T, Soderblom T, Backhed F, Celsi G, Brismar H, Normark S, Aperia A, Richter-Dahlfors A: Alpha-haemolysin of uropathogenic E. coli induces Ca2+ oscillations in renal epithelial cells. Nature. 2000, 405: 694-697. 10.1038/35015091.
Nagy G, Altenhoefer A, Knapp O, Maier E, Dobrindt U, Blum-Oehler G, Benz R, Emody L, Hacker J: Both alpha-haemolysin determinants contribute to full virulence of uropathogenic Escherichia coli strain 536. Microbes Infect. 2006, 8: 2006-2012. 10.1016/j.micinf.2006.02.029.
Johnson JR, Delavari P, Kuskowski M, Stell AL: Phylogenetic distribution of extraintestinal virulence-associated traits in Escherichia coli. J Infect Dis. 2001, 183: 78-88. 10.1086/317656.
Bingen-Bidois M, Clermont O, Bonacorsi S, Terki M, Brahimi N, Loukil C, Barraud D, Bingen E: Phylogenetic analysis and prevalence of urosepsis strains of Escherichia coli bearing pathogenicity island-like domains. Infect Immun. 2002, 70: 3216-3226. 10.1128/IAI.70.6.3216-3226.2002.
Blanco J, Blanco M, Alonso MP, Blanco JE, Garabal JI, Gonzalez EA: Serogroups of Escherichia coli strains producing cytotoxic necrotizing factors CNF1 and CNF2. FEMS Microbiol Lett. 1992, 75: 155-159. 10.1111/j.1574-6968.1992.tb05409.x.
Landraud L, Gauthier M, Fosse T, Boquet P: Frequency of Escherichia coli strains producing the cytotoxic necrotizing factor (CNF1) in nosocomial urinary tract infections. Lett Appl Microbiol. 2000, 30: 213-216. 10.1046/j.1472-765x.2000.00698.x.
Blum G, Falbo V, Caprioli A, Hacker J: Gene clusters encoding the cytotoxic necrotizing factor type 1, Prs-fimbriae and alpha-hemolysin form the pathogenicity island II of the uropathogenic Escherichia coli strain J96. FEMS Microbiol Lett. 1995, 126: 189-195.
Swenson DL, Bukanov NO, Berg DE, Welch RA: Two pathogenicity islands in uropathogenic Escherichia coli J96: cosmid cloning and sample sequencing. Infect Immun. 1996, 64: 3736-3743.
Boyd EF, Hartl DL: Diversifying selection governs sequence polymorphism in the major adhesin proteins fimA, papA, and sfaA of Escherichia coli. J Mol Evol. 1998, 47: 258-267. 10.1007/PL00006383.
Ochman H, Selander RK: Standard reference strains of Escherichia coli from natural populations. J Bacteriol. 1984, 157: 690-693.
Escobar-Paramo P, Sabbagh A, Darlu P, Pradillon O, Vaury C, Denamur E, Lecointre G: Decreasing the effects of horizontal gene transfer on bacterial phylogeny: the Escherichia coli case study. Mol Phylogenet Evol. 2004, 30: 243-250. 10.1016/S1055-7903(03)00181-7.
Bjedov I, Lecointre G, Tenaillon O, Vaury C, Radman M, Taddei F, Denamur E, Matic I: Polymorphism of genes encoding SOS polymerases in natural populations of Escherichia coli. DNA Repair (Amst). 2003, 2: 417-426. 10.1016/S1568-7864(02)00241-0.
Lecointre G, Rachdi L, Darlu P, Denamur E: Escherichia coli molecular phylogeny using the incongruence length difference test. Mol Biol Evol. 1998, 15: 1685-1695.
Jolley KA, Feil EJ, Chan MS, Maiden MC: Sequence type analysis and recombinational tests (START). Bioinformatics. 2001, 17: 1230-1231. 10.1093/bioinformatics/17.12.1230.
Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.
Didelot X, Falush D: Inference of bacterial microevolution using multilocus sequence data. Genetics. 2007, 175: 1251-1266. 10.1534/genetics.106.063305.
Hochhut B, Wilde C, Balling G, Middendorf B, Dobrindt U, Brzuszkiewicz E, Gottschalk G, Carniel E, Hacker J: Role of pathogenicity island-associated integrases in the genome plasticity of uropathogenic Escherichia coli strain 536. Mol Microbiol. 2006, 61: 584-595. 10.1111/j.1365-2958.2006.05255.x.
Blattner FR, Plunkett G, Bloch CA, Perna NT, Burland V, Riley M, Collado-Vides J, Glasner JD, Rode CK, Mayhew GF: The complete genome sequence of Escherichia coli K-12. Science. 1997, 277: 1453-1474. 10.1126/science.277.5331.1453.
Xie Y, Kolisnychenko V, Paul-Satyaseela M, Elliott S, Parthasarathy G, Yao Y, Plunkett G, Blattner FR, Kim KS: Identification and characterization of Escherichia coli RS218-derived islandsin the pathogenesis of E. coli meningitis. J Infect Dis. 2006, 194: 358-364. 10.1086/505429.
Lepage E, Brinster S, Caron C, Ducroix-Crepy C, Rigottier-Gois L, Dunny G, Hennequet-Antier C, Serror P: Comparative genomic hybridization analysis of Enterococcus faecalis: identification of genes absent from food strains. J Bacteriol. 2006, 188: 6858-6868. 10.1128/JB.00421-06.
Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34: 374-378.
Platform Genotyping of Pathogens and Public Health acknowledges financial support from the Institut Pasteur and the Institut de Veille Sanitaire (Saint-Maurice, France). XN was supported by INSERM and University Paris Descartes. ED was partially supported by the "Fondation pour la Recherche Médicale".
FJ, LL, OL, BP, XN and SB designed the study. VP, LD, EC and OC conceived the experiments, and gathered data and isolates. FJ, LL, VP, LD, EF, GG, BP, XN and SB analyzed the data. FJ, LL, LD, GG, OL, ED, BP, XN and SB wrote the manuscript. All authors read and approved the final manuscript.
Françoise Jaureguy, Luce Landraud contributed equally to this work.