Single-nucleotide polymorphism phylotyping of Escherichia coli

We describe a rapid and easily automated phylogenetic grouping technique based on analysis of bacterial genome single-nucleotide polymorphisms (SNPs). We selected 13 SNPs derived from a complete sequence analysis of 11 essential genes previously used for multilocus sequence typing (MLST) of 30 Escherichia coli strains representing the genetic diversity of the species. The 13 SNPs were localized in five genes, trpA, trpB, putP, icdA, and polB, and were selected to allow recovery of the main phylogenetic groups (groups A, B1, E, D, and B2) and subgroups of the species. In the first step, we validated the SNP approach in silico by extracting SNP data from the complete sequences of the five genes for a panel of 65 pathogenic strains belonging to different E. coli pathovars, which were previously analyzed by MLST. In the second step, we determined these SNPs by dideoxy single-base extension of unlabeled oligonucleotide primers for a collection of 183 commensal and extraintestinal clinical E. coli isolates and compared the SNP phylotyping method to previous well-established typing methods. This SNP phylotyping method proved to be consistent with the other methods for assigning phylogenetic groups to the different E. coli strains. In contrast to the other typing methods, such as multilocus enzyme electrophoresis, ribotyping, or PCR phylotyping using the presence/absence of three genomic DNA fragments, the SNP typing method described here is derived from a solid phylogenetic analysis, and the results obtained by this method are more meaningful. Our results indicate that similar approaches may be used for a wide variety of bacterial species.     

Houseflies (Musca domestica) as Vectors for Extended-Spectrum β-Lactamase-Producing Escherichia coli on Spanish Broiler Farms

Flies may act as potential vectors for the spread of resistant bacteria to different environments. This study was intended to evaluate the presence of Escherichia coli strains resistant to cephalosporins in flies captured in the areas surrounding five broiler farms. Phenotypic and molecular characterization of the resistant population was performed by different methods: MIC determination, pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), and phylotyping. The presence of extended-spectrum beta-lactamase (ESBL) genes, their plasmid location, and the mobile genetic elements involved in their mobilization were studied. Additionally, the presence of 35 genes associated with virulence was evaluated. Out of 682 flies captured, 42 yielded ESBL-producing E. coli. Of these isolates, 23 contained bla_CTX-M-1, 18 contained bla_CTX-M-14, and 1 contained bla_CTX-M-9. ESBL genes were associated mainly with the presence of the IncI1 and IncFIB replicons. Additionally, all the strains were multiresistant, and five of them also harbored qnrS. Identical PFGE profiles were found for E. coli isolates obtained from flies at different sampling times, indicating a persistence of the same clones in the farm environment over months. According to their virulence genes, 81% of the isolates were considered avian-pathogenic E. coli (APEC) and 29% were considered extraintestinal pathogenic E. coli (ExPEC). The entrance of flies into broiler houses constitutes a considerable risk for colonization of broilers with multidrug-resistant E. coli. ESBLs in flies reflect the contamination status of the farm environment. Additionally, this study demonstrates the potential contribution of flies to the dissemination of virulence and resistance genes into different ecological niches.     

Multilocus Sequence Typing and Phylogenetic Analyses of Pseudomonas aeruginosa Isolates from the Ocean

Recent isolation of Pseudomonas aeruginosa strains from the open ocean and subsequent pulsed-field gel electrophoresis analyses indicate that these strains have a unique genotype (N. H. Khan, Y. Ishii, N. Kimata-Kino, H. Esaki, T. Nishino, M. Nishimura, and K. Kogure, Microb. Ecol. 53:173-186, 2007). We hypothesized that ocean P. aeruginosa strains have a unique phylogenetic position relative to other strains. The objective of this study was to clarify the intraspecies phylogenetic relationship between marine strains and other strains from various geographical locations. Considering the advantages of using databases, multilocus sequence typing (MLST) was chosen for the typing and discrimination of ocean P. aeruginosa strains. Seven housekeeping genes (acsA, aroE, guaA, mutL, nuoD, ppsA, and trpE) were analyzed, and the results were compared with data on the MLST website. These genes were also used for phylogenetic analysis of P. aeruginosa. Rooted and unrooted phylogenetic trees were generated for each gene locus and the concatenated gene fragments. MLST data showed that all the ocean strains were new. Trees constructed for individual and concatenated genes revealed that ocean P. aeruginosa strains have clusters distinct from those of other P. aeruginosa strains. These clusters roughly reflected the geographical locations of the isolates. These data support our previous findings that P. aeruginosa strains are present in the ocean. It can be concluded that the ocean P. aeruginosa strains have diverged from other isolates and form a distinct cluster based on MLST and phylogenetic analyses of seven housekeeping genes.     

Multicolor Melting Curve Analysis-Based Multilocus Melt Typing of Vibrio parahaemolyticus

Vibrio parahaemolyticus is the leading cause of seafood-borne gastroenteritis outbreaks. To track the source of these diseases in a timely manner, a high throughput typing method is critical. We hereby describe a novel genotyping method for V. parahaemolyticus, termed multilocus melt typing (MLMT), based on multilocus sequence typing (MLST). MLMT utilizes melting curve analysis to interrogate the allelic types of a set of informative single nucleotide polymorphisms (SNPs) derived from the housekeeping genes used in MLST. For each SNP, one allelic type generates distinct T_m values, which are converted into a binary code. Multiple SNPs thus generate a series of binary codes, forming a melt type (MT) corresponding with a sequence type (ST) of MLST. Using a set of 12 SNPs, the MLMT scheme could resolve 218 V.parahaemolyticus isolates into 50 MTs corresponding with 56 STs. The discriminatory power of MLMT and MLST was similar with Simpson’s index of diversity of 0.638 and 0.646, respectively. The global (adjusted Rand index = 0.982) and directional congruence (adjusted Wallace coefficient, MT→ST = 0.965; ST→MT = 1.000) between the two typing approaches was high. The entire procedure of MLMT could be finished within 3 h with negligible hands on time in a real-time PCR machine. We conclude that MLMT provides a reliable and efficient approach for V. parahaemolyticus genotyping and might also find use in other pathogens.     

Amplified Intergenic Locus Polymorphism as a Basis for Bacterial Typing of Listeria spp. and Escherichia coli

DNA-based methods are increasingly important for bacterial typing. The high number of polymorphic sites present among closely related bacterial genomes is the basis for the presented method. The method identifies multilocus genomic polymorphisms in intergenic regions termed AILP (amplified intergenic locus polymorphism). For each locus, a pair of unique PCR primers was designed to amplify an intergenic sequence from one open reading frame (ORF) to the adjacent ORF. Presence, absence, and size variation of the amplification products were identified and used as genetic markers for rapidly differentiating among strains. Polymorphism was evaluated using 18 AILP sites among 28 strains of Listeria monocytogenes and 6 strains of Listeria spp. and 30 AILP markers among 27 strains of Escherichia coli. Up to four alleles per locus were identified among Listeria strains, and up to six were identified among E. coli strains. In both species, more than half of the AILP sites revealed intraspecies polymorphism. The AILP data were applied to phylogenetic analysis among Listeria and E. coli strains. A clear distinction between L. monocytogenes and Listeria spp. was demonstrated. In addition, the method separated L. monocytogenes into the three known lineages and discriminated the most common virulent serotypic group, 4b. In E. coli, AILP analysis separated the known groups as well as the virulent O157:H7 isolates. These findings for both Listeria and E. coli are in agreement with other phylogenetic studies using molecular markers. The AILP method was found to be rapid, simple, reproducible, and a low-cost method for initial bacterial typing that could serve as a basis for epidemiological investigation.     

rhs Genes Are Potential Markers for Multilocus Sequence Typing of Escherichia coli O157:H7 Strains

DNA sequence-based molecular subtyping methods such as multilocus sequence typing (MLST) are commonly used to generate phylogenetic inferences for monomorphic pathogens. The development of an effective MLST scheme for subtyping Escherichia coli O157:H7 has been hindered in the past due to the lack of sequence variation found within analyzed housekeeping and virulence genes. A recent study suggested that rhs genes are under strong positive selection pressure, and therefore in this study we analyzed these genes within a diverse collection of E. coli O157:H7 strains for sequence variability. Eighteen O157:H7 strains from lineages I and II and 15 O157:H7 strains from eight clades were included. Examination of these rhs genes revealed 44 polymorphic loci (PL) and 10 sequence types (STs) among the 18 lineage strains and 280 PL and 12 STs among the 15 clade strains. Phylogenetic analysis using rhs genes generally grouped strains according to their known lineage and clade classifications. These findings also suggested that O157:H7 strains from clades 6 and 8 fall into lineage I/II and that strains of clades 1, 2, 3, and 4 fall into lineage I. Additionally, unique markers were found in rhsA and rhsJ that might be used to define clade 8 and clade 6. Therefore, rhs genes may be useful markers for phylogenetic analysis of E. coli O157:H7.Escherichia coli O157:H7 was first described in 1983 as the causative agent of a food-borne outbreak attributed to contaminated ground beef patties (35), and it has subsequently emerged as a very important food-borne pathogen. Diseases caused by E. coli O157:H7, such as hemorrhagic colitis and hemolytic uremic syndrome, can be very severe or even life-threatening. Cattle are believed to be the main reservoir for E. coli O157:H7 (5, 15, 41), although other animals may also carry this organism (6, 21). Outbreaks are commonly associated with the consumption of beef and fresh produce that come into contact with bovine feces or feces-contaminated environments, such as food contact surfaces, animal hides, or irrigation water (12, 21, 30, 38).It is well-established that strains of E. coli O157:H7 vary in terms of virulence and transmissibility to humans and that strains differing in these characteristics can be distinguished using DNA-based methods (22, 29, 42). For example, octamer-based genome scanning, which is a PCR approach using 8-bp primers, provided the first evidence that there are at least two lineages of O157:H7, termed lineage I and lineage II (22). Strains classified as lineage I are more frequently isolated from humans than are lineage II strains (42). A later refinement of this classification system was coined the lineage-specific polymorphism assay (LSPA), which classified strains based upon the amplicon size obtained using PCRs targeting six chromosomal regions of E. coli O157:H7 and assigned a six-digit code based upon the pattern obtained (42). Most strains of lineage I grouped into LSPA type 111111, while the majority of lineage II strains fell into LSPA types 211111, 212111, and 222222. More recently, it was suggested that LSPA type 211111 strains comprise a separate group called lineage I/II (45).To gain greater insight into the recent evolution of E. coli O157:H7, a method that is more discriminatory than the LSPA method is desirable. Multilocus sequence typing (MLST) is a method that discriminates between strains of a bacterial species by identifying DNA sequence differences in six to eight targeted genes. Satisfactory MLST schemes exist for other bacterial pathogens (28, 43); however, due to the lack of sequence variations in previously targeted gene markers in E. coli O157:H7 (13, 33), MLST approaches for subtyping this pathogen have been more difficult to develop. More recently, high-throughput microarray and sequencing platforms have been used to identify hundreds of single nucleotide polymorphisms (SNPs) that are useful for discriminating between strains of E. coli O157:H7 during epidemiologic investigations and for drawing phylogenetic inferences (11, 20, 29, 44). Particularly noteworthy, Manning et al. (29) developed a subtyping scheme based upon the interrogation of 32 putative SNP loci. This method separated 528 strains into 39 distinct SNP genotypes, which were grouped into nine statistically supported phylogenetic groups called clade 1 through clade 9. By analyzing the rates of hemolytic uremic syndrome observed in patients infected with strains of clades 2, 7, and 8, it was also concluded that clade 8 strains are more virulent to humans than other strains (29).One drawback of current DNA sequence-based subtyping schemes for E. coli O157:H7 is that they require screening of at least 32 SNP loci. We were interested in asking whether a simpler approach that targets a few informative gene markers could be developed for rapid strain discrimination and phylogenetic determination. A recent analysis of E. coli genomes predicted that rearrangement hot spot (rhs) genes are under the strongest positive selection of all coding sequences analyzed (34). Therefore, we hypothesized that these genes would display significant sequence variations for subtyping O157:H7 strains. The rhs genes were first discovered as elements mediating tandem duplication of the glyS locus in E. coli K-12 (26); however, their function remains unknown. There are nine rhs genes within the genome of the prototypical E. coli O157:H7 strain Sakai, and these genes are designated rhsA, -C, -D, -E, -F, -G, -I, -J, and -K (see Table S1 in the supplemental material) (16). Three of these nine rhs genes, rhsF, -J, and -K, were previously studied by Zhang et al. (44), and a number of SNPs were identified among these genes. However, no studies have been conducted to comprehensively investigate rhs genes as markers in an MLST scheme for subtyping E. coli O157:H7.The primary purpose of the present study was to investigate whether there are sufficient DNA sequence variations among rhs genes to develop an MLST approach for subtyping E. coli O157:H7. In this study, a greater level of DNA sequence variation was observed among rhs genes than in gene markers targeted in previous studies (13, 33). Furthermore, phylogenetic analysis using these rhs genes generally agreed with the established lineage and clade classifications of O157:H7 strains defined previously. We also wanted to determine whether there is a correlation between the lineage classification of O157:H7 strains (42) and the recently proposed clade classification (29). The present study reports evidence that O157:H7 strains from clade 8 are classified as lineage I/II, which is a different lineage from well-studied E. coli O157:H7 outbreak strains, such as EDL933 and Sakai. Therefore, we suggest that outbreaks of O157:H7 are caused by two lineages of this pathogen, lineage I and lineage I/II.     

Typing of Campylobacter jejuni and Campylobacter coli isolated from live broilers and retail broiler meat by flaA-RFLP, MLST, PFGE and REP-PCR

We analyzed 100 Campylobacter spp. isolates (C. jejuni and C. coli) from Grenada, Puerto Rico and Alabama, which were collected from live broilers or retail broiler meat. We analyzed these isolates with four molecular typing methods: restriction fragment length polymorphism of the flaA gene (flaA-RFLP), multilocus sequence typing (MLST), pulsed-field gel electrophoresis (PFGE), and automated repetitive extragenic palindromic polymerase chain reaction (REP-PCR) using the DiversiLab system. All methods performed similarly for the typing of C. jejuni and C. coli. The DNA extraction method appears to influence the results obtained with REP-PCR. This method was better for the typing of C. jejuni than C. coli, however both REP-PCR and flaA-RFLP generated types that were indistinguishable between C. jejuni and C. coli and appeared to be random, without any relationship to species, location, or source of isolates. PFGE and MLST generated typing results that had a better correlation with the geographic location of the isolates and showed higher concordance with the Wallace coefficient. The adjusted Rand coefficient did not show higher concordance among the methods, although the PFGE/MLST combination exhibited the highest concordance. PFGE and MLST revealed a better discriminatory power for C. coli isolates than REP-PCR or flaA-RFLP. The use of readily available online tools to calculate the confidence interval of the Simpson's index of diversity and the adjusted Rand and Wallace coefficients helped estimate the discriminatory power of typing methods. Further studies using different C. jejuni and C. coli strains may expand our understanding of the benefits and limitations of each of these typing methods for epidemiological studies of Campylobacter spp.     

Clostridium botulinum Group II Isolate Phylogenomic Profiling Using Whole-Genome Sequence Data

Clostridium botulinum group II isolates (n = 163) from different geographic regions, outbreaks, and neurotoxin types and subtypes were characterized in silico using whole-genome sequence data. Two clusters representing a variety of botulinum neurotoxin (BoNT) types and subtypes were identified by multilocus sequence typing (MLST) and core single nucleotide polymorphism (SNP) analysis. While one cluster included BoNT/B4/F6/E9 and nontoxigenic members, the other comprised a wide variety of different BoNT/E subtype isolates and a nontoxigenic strain. In silico MLST and core SNP methods were consistent in terms of clade-level isolate classification; however, core SNP analysis showed higher resolution capability. Furthermore, core SNP analysis correctly distinguished isolates by outbreak and location. This study illustrated the utility of next-generation sequence-based typing approaches for isolate characterization and source attribution and identified discrete SNP loci and MLST alleles for isolate comparison.     

Molecular serotyping of Salmonella enterica by complete rpoB gene sequencing

Serotyping has been the gold standard for identifying Salmonella, but it requires large amounts of standard antisera. Multilocus sequence typing (MLST) has been applied to identify Salmonella serovars, but the recombination of 4–7 housekeeping genes and multiple analytic steps diminish its applicability. In the present study, we determined the complete sequences of the RNA polymerase beta subunit gene (rpoB) and 7 housekeeping genes (aroC, dnaN, hemD, hisD, purE, sucA, and thrA) for 76 strains of 33 Salmonella enterica serovars and conducted phylogenetic analyses together with the corresponding gene sequences of 24 reference strains registered in the GenBank database. Based on the phylogenetic analyses, 100 strains from 40 serovars and 91 strains from 37 serovars were classified into 60 rpoB (RST) and 49 multilocus sequence types (ST), respectively. The nucleotide similarities were 98.8–100% and 96.9–100% for the complete rpoB gene and the seven concatenated housekeeping genes, respectively. The strains of 35 and 30 serovars formed serovar-specific branches or clusters in the rpoB and housekeeping gene phylogenetic trees, respectively. Therefore, complete rpoB gene sequencing and phylogenetic analysis may be a useful method for identifying Salmonella serovars that is a simpler, more cost-effective, and less time-consuming alternative or complementary method to MLST and conventional serotyping.     

Genetic Affinities within a Large Global Collection of Pathogenic Leptospira: Implications for Strain Identification and Molecular Epidemiology

Leptospirosis is an important zoonosis with widespread human health implications. The non-availability of accurate identification methods for the individualization of different Leptospira for outbreak investigations poses bountiful problems in the disease control arena. We harnessed fluorescent amplified fragment length polymorphism analysis (FAFLP) for Leptospira and investigated its utility in establishing genetic relationships among 271 isolates in the context of species level assignments of our global collection of isolates and strains obtained from a diverse array of hosts. In addition, this method was compared to an in-house multilocus sequence typing (MLST) method based on polymorphisms in three housekeeping genes, the rrs locus and two envelope proteins. Phylogenetic relationships were deduced based on bifurcating Neighbor-joining trees as well as median joining network analyses integrating both the FAFLP data and MLST based haplotypes. The phylogenetic relationships were also reproduced through Bayesian analysis of the multilocus sequence polymorphisms. We found FAFLP to be an important method for outbreak investigation and for clustering of isolates based on their geographical descent rather than by genome species types. The FAFLP method was, however, not able to convey much taxonomical utility sufficient to replace the highly tedious serotyping procedures in vogue. MLST, on the other hand, was found to be highly robust and efficient in identifying ancestral relationships and segregating the outbreak associated strains or otherwise according to their genome species status and, therefore, could unambiguously be applied for investigating phylogenetics of Leptospira in the context of taxonomy as well as gene flow. For instance, MLST was more efficient, as compared to FAFLP method, in clustering strains from the Andaman island of India, with their counterparts from mainland India and Sri Lanka, implying that such strains share genetic relationships and that leptospiral strains might be frequently circulating between the islands and the mainland.     

Development of a Tiered Multilocus Sequence Typing Scheme for Members of the Lactobacillus acidophilus Complex

Members of the Lactobacillus acidophilus complex are associated with functional foods and dietary supplements because of purported health benefits they impart to the consumer. Many characteristics of these microorganisms are reported to be strain specific. Therefore, proper strain typing is essential for safety assessment and product labeling, and also for monitoring strain integrity for industrial production purposes. Fifty-two strains of the L. acidophilus complex (L. acidophilus, L. amylovorus, L. crispatus, L. gallinarum, L. gasseri, and L. johnsonii) were genotyped using two established methods and compared to a novel multilocus sequence typing (MLST) scheme. PCR restriction fragment length polymorphism (PCR-RFLP) analysis of the hsp60 gene with AluI and TaqI successfully clustered 51 of the 52 strains into the six species examined, but it lacked strain-level discrimination. Random amplified polymorphic DNA PCR (RAPD-PCR) targeting the M13 sequence resulted in highly discriminatory profiles but lacked reproducibility. In this study, an MLST scheme was developed using the conserved housekeeping genes fusA, gpmA, gyrA, gyrB, lepA, pyrG, and recA, which identified 40 sequence types that successfully clustered all of the strains into the six species. Analysis of the observed alleles suggests that nucleotide substitutions within five of the seven MLST loci have reached saturation, a finding that emphasizes the highly diverse nature of the L. acidophilus complex and our unconventional application of a typically intraspecies molecular typing tool. Our MLST results indicate that this method could be useful for characterization and strain discrimination of a multispecies complex, with the potential for taxonomic expansion to a broader collection of Lactobacillus species.     

DNA identification and characterization of Campylobacter jejuni and Campylobacter coli isolated from caecal samples of chickens in Grenada

Aims: To speciate Campylobacter strains from the caeca of chickens in Grenada using PCR and to evaluate DNA‐based typing methods for the characterization of these isolates. Methods and Results: Isolates were speciated with two multiplex PCR assays and were typed with flaA‐RFLP, pulsed‐field gel electrophoresis (PFGE) and multilocus sequence typing (MLST). Results confirmed that Campylobacter coli strains were more predominant than Campylobacter jejuni strains. From 56 isolates, 18 were misidentified using biochemical tests. PFGE typing gave the highest discriminatory power among the methods used (Simpson’s index of diversity, D = 0·9061). However, the combination of flaA‐RFLP, PFGE and MLST results gave the highest discrimination for subtyping of these isolates (D = 0·9857). A band position tolerance of 4% in Bio Numerics was the most appropriate for the analysis of this database. MLST profiles were generally concordant with PFGE and/or flaA‐RFLP types. Several isolates exhibited new MLST sequence types (STs), and 43 of the 49 Camp. coli strains belonged to the ST‐828 clonal complex. Conclusions: Campylobacter coli was the most prevalent species isolated from broilers and layers in Grenada, and a combination of restriction and sequence methods was most appropriate for the typing of Camp. coli isolates. Campylobacter coli STs clustered with described poultry‐associated Camp. coli STs by phylogenetic analysis. Significance and Impact of the Study: Further studies to understand the predominance of Camp. coli within Campylobacter spp. from chickens in Grenada may help elucidate the epidemiology of these pathogens in chickens.     

Effect of a glucose impulse on the CcpA regulon in Staphylococcus aureus

Background

A low genetic diversity in Francisella tularensis has been documented. Current DNA based genotyping methods for typing F. tularensis offer a limited and varying degree of subspecies, clade and strain level discrimination power. Whole genome sequencing is the most accurate and reliable method to identify, type and determine phylogenetic relationships among strains of a species. However, lower cost typing schemes are necessary in order to enable typing of hundreds or even thousands of isolates.

Results

We have generated a high-resolution phylogenetic tree from 40 Francisella isolates, including 13 F. tularensis subspecies holarctica (type B) strains, 26 F. tularensis subsp. tularensis (type A) strains and a single F. novicida strain. The tree was generated from global multi-strain single nucleotide polymorphism (SNP) data collected using a set of six Affymetrix GeneChip^® resequencing arrays with the non-repetitive portion of LVS (type B) as the reference sequence complemented with unique sequences of SCHU S4 (type A). Global SNP based phylogenetic clustering was able to resolve all non-related strains. The phylogenetic tree was used to guide the selection of informative SNPs specific to major nodes in the tree for development of a genotyping assay for identification of F. tularensis subspecies and clades. We designed and validated an assay that uses these SNPs to accurately genotype 39 additional F. tularensis strains as type A (A1, A2, A1a or A1b) or type B (B1 or B2).

Conclusion

Whole-genome SNP based clustering was shown to accurately identify SNPs for differentiation of F. tularensis subspecies and clades, emphasizing the potential power and utility of this methodology for selecting SNPs for typing of F. tularensis to the strain level. Additionally, whole genome sequence based SNP information gained from a representative population of strains may be used to perform evolutionary or phylogenetic comparisons of strains, or selection of unique strains for whole-genome sequencing projects.     

Multi-Virulence-Locus Sequence Typing of Listeria monocytogenes

A multi-virulence-locus sequence typing (MVLST) scheme was developed for subtyping Listeria monocytogenes, and the results obtained using this scheme were compared to those of pulsed-field gel electrophoresis (PFGE) and the published results of other typing methods, including ribotyping (RT) and multilocus sequence typing (MLST). A set of 28 strains (eight different serotypes and three known genetic lineages) of L. monocytogenes was selected from a strain collection (n > 1,000 strains) to represent the genetic diversity of this species. Internal fragments (ca. 418 to 469 bp) of three virulence genes (prfA, inlB, and inlC) and three virulence-associated genes (dal, lisR, and clpP) were sequenced and analyzed. Multiple DNA sequence alignment identified 10 (prfA), 19 (inlB), 13 (dal), 10 (lisR), 17 (inlC), and 16 (clpP) allelic types and a total of 28 unique sequence types. Comparison of MVLST with automated EcoRI-RT and PFGE with ApaI enzymatic digestion showed that MVLST was able to differentiate strains that were indistinguishable by RT (13 ribotypes; discrimination index = 0.921) or PFGE (22 profiles; discrimination index = 0.970). Comparison of MVLST with housekeeping-gene-based MLST analysis showed that MVLST provided higher discriminatory power for serotype 1/2a and 4b strains than MLST. Cluster analysis based on the intragenic sequences of the selected virulence genes indicated a strain phylogeny closely related to serotypes and genetic lineages. In conclusion, MVLST may improve the discriminatory power of MLST and provide a convenient tool for studying the local epidemiology of L. monocytogenes.     

Genetic diversity of Staphylococcus equorum isolates from Saeu-jeotgal evaluated by multilocus sequence typing

Staphylococcus equorum, the predominant bacterial species detected in Saeu-jeotgal, a Korean high-salt fermented seafood, is a candidate starter bacterium for Saeu-jeotgal fermentation. A multilocus sequence typing (MLST) scheme was developed to evaluate the genetic diversity and background of S. equorum strains isolated from Saeu-jeotgal. A total of 135 strains, including 117 isolates from Saeu-jeotgal, and others from Myeolchi-jeotgal, sausage, cheese and horse skin, were subjected to MLST, and the internal fragments of seven housekeeping genes, aroE, dnaJ, glpF, gmk, hsp60, mutS, and pta, were compared. This MLST scheme produced 45 sequence types (STs) and the eBURST algorithm clustered the STs into nine clonal groups and seven singletons. Clonal group 1, the major group, consisted of 30 isolates from cheese, Saeu-jeotgal and sausages, which were classified into 12 STs. The predominant ST, ST26, comprised 25 isolates and presented as a singleton. Most of the isolates from Myeolchi-jeotgal and sausages clustered on two different branches of a phylogenetic tree generated with a cluster analysis using the maximum likelihood algorithm. This MLST scheme established the genetic backgrounds of S. equorum strains isolated from different types of food. Among the housekeeping genes used for MLST, gmk had the fewest allele types and fairly low sequence identities (74.0–90.0 %) within the Staphylococcus species. Therefore, sequence analyses of the gmk gene and 16S rRNA gene can be used for the accurate and rapid identification of S. equorum.     

Frequent Combination of Antimicrobial Multiresistance and Extraintestinal Pathogenicity in Escherichia coli Isolates from Urban Rats (Rattus norvegicus) in Berlin,Germany

Urban rats present a global public health concern as they are considered a reservoir and vector of zoonotic pathogens, including Escherichia coli. In view of the increasing emergence of antimicrobial resistant E. coli strains and the on-going discussion about environmental reservoirs, we intended to analyse whether urban rats might be a potential source of putatively zoonotic E. coli combining resistance and virulence. For that, we took fecal samples from 87 brown rats (Rattus norvegicus) and tested at least three E. coli colonies from each animal. Thirty two of these E. coli strains were pre-selected from a total of 211 non-duplicate isolates based on their phenotypic resistance to at least three antimicrobial classes, thus fulfilling the definition of multiresistance. As determined by multilocus sequence typing (MLST), these 32 strains belonged to 24 different sequence types (STs), indicating a high phylogenetic diversity. We identified STs, which frequently occur among extraintestinal pathogenic E. coli (ExPEC), such as STs 95, 131, 70, 428, and 127. Also, the detection of a number of typical virulence genes confirmed that the rats tested carried ExPEC-like strains. In particular, the finding of an Extended-spectrum beta-lactamase (ESBL)-producing strain which belongs to a highly virulent, so far mainly human- and avian-restricted ExPEC lineage (ST95), which expresses a serogroup linked with invasive strains (O18:NM:K1), and finally, which produces an ESBL-type frequently identified among human strains (CTX-M-9), pointed towards the important role, urban rats might play in the transmission of multiresistant and virulent E. coli strains. Indeed, using a chicken infection model, this strain showed a high in vivo pathogenicity. Imagining the high numbers of urban rats living worldwide, the way to the transmission of putatively zoonotic, multiresistant, and virulent strains might not be far ahead. The unforeseeable consequences of such an emerging public health threat need careful consideration in the future.     

Phylogeny and Strain Typing of Escherichia coli, Inferred from Variation at Mononucleotide Repeat Loci

Multilocus sequencing of housekeeping genes has been used previously for bacterial strain typing and for inferring evolutionary relationships among strains of Escherichia coli. In this study, we used shorter intergenic sequences that contained simple sequence repeats (SSRs) of repeating mononucleotide motifs (mononucleotide repeats [MNRs]) to infer the phylogeny of pathogenic and commensal E. coli strains. Seven noncoding loci (four MNRs and three non-SSRs) were sequenced in 27 strains, including enterohemorrhagic (six isolates of O157:H7), enteropathogenic, enterotoxigenic, B, and K-12 strains. The four MNRs were also sequenced in 20 representative strains of the E. coli reference (ECOR) collection. Sequence polymorphism was significantly higher at the MNR loci, including the flanking sequences, indicating a higher mutation rate in the sequences flanking the MNR tracts. The four MNR loci were amplifiable by PCR in the standard ECOR A, B1, and D groups, but only one (yaiN) in the B2 group was amplified, which is consistent with previous studies that suggested that B2 is the most ancient group. High sequence compatibility was found between the four MNR loci, indicating that they are in the same clonal frame. The phylogenetic trees that were constructed from the sequence data were in good agreement with those of previous studies that used multilocus enzyme electrophoresis. The results demonstrate that MNR loci are useful for inferring phylogenetic relationships and provide much higher sequence variation than housekeeping genes. Therefore, the use of MNR loci for multilocus sequence typing should prove efficient for clinical diagnostics, epidemiology, and evolutionary study of bacteria.     

Recombination Blurs Phylogenetic Groups Routine Assignment in Escherichia coli: Setting the Record Straight

The characterization of population structures plays a main role for understanding outbreaks and the dynamics of bacterial spreading. In Escherichia coli, the widely used combination of multiplex-PCR scheme together with goeBURST has some limitations. The purpose of this study is to show that the combination of different phylogenetic approaches based on concatenated sequences of MLST genes results in a more precise assignment of E. coli phylogenetic groups, complete understanding of population structure and reconstruction of ancestral clones. A collection of 80 Escherichia coli strains of different origins was analyzed following the Clermont and Doumith''s multiplex-PCR schemes. Doumith''s multiplex-PCR showed only 1.7% of misassignment, whereas Clermont''s-2000 protocol reached 14.0%, although the discrepancies reached 30% and 38.7% respectively when recombinant C, F and E phylogroups were considered. Therefore, correct phylogroup attribution is highly variable and depends on the clonal composition of the sample. As far as population structure of these E. coli strains, including 48 E. coli genomes from GenBank, goeBURST provides a quite dispersed population structure; whereas NeighborNet approach reveals a complex population structure. MLST-based eBURST can infer different founder genotypes, for instance ST23/ST88 could be detected as the founder genotypes for STC23; however, phylogenetic reconstructions might suggest ST410 as the ancestor clone and several evolutionary trajectories with different founders. To improve our routine understanding of E. coli molecular epidemiology, we propose a strategy based on three successive steps; first, to discriminate three main groups A/B1/C, D/F/E and B2 following Doumith''s protocol; second, visualization of population structure based on MLST genes according to goeBURST, using NeighborNet to establish more complex relationships among STs; and third, to perform, a cost-free characterization of evolutionary trajectories in variants emerging along the clonal expansion using parsimony methods of phylogenetic analysis.     

Correlation between Ureaplasma Subgroup 2 and Genitourinary Tract Disease Outcomes Revealed by an Expanded Multilocus Sequence Typing (eMLST) Scheme

The multilocus sequence typing (MLST) scheme of Ureaplasma based on four housekeeping genes (ftsH, rpL22, valS, and thrS) was described in our previous study; here we introduced an expanded MLST (eMLST) scheme with improved discriminatory power, which was developed by adding two putative virulence genes (ureG and mba-np1) to the original MLST scheme. To evaluate the discriminatory power of eMLST, a total of 14 reference strains of Ureaplasma serovars and 269 clinical strains (134 isolated from symptomatic patients and 135 obtained from asymptomatic persons) were investigated. Our study confirmed that all 14 serotype strains could successfully be differentiated into 14 eMLST STs (eSTs), while some of them could not even be differentiated by the MLST, and a total of 136 eSTs were identified among the clinical isolates we investigated. In addition, phylogenetic analysis indicated that two genetically significantly distant clusters (cluster I and II) were revealed and most clinical isolates were located in cluster I. These findings were in accordance with and further support for the concept of two well-known genetic lineages (Ureaplasma parvum and Ureaplasma urealyticum) in our previous study. Interestingly, although both clusters were associated with clinical manifestation, the sub-group 2 of cluster II had pronounced and adverse effect on patients and might be a potential risk factor for clinical outcomes. In conclusion, the eMLST scheme offers investigators a highly discriminative typing tool that is capable for precise epidemiological investigations and clinical relevance of Ureaplasma.     

A phylogenetic group of Escherichia coli associated with active left-sided Inflammatory Bowel Disease

Background  

Escherichia coli have been found in increased numbers in tissues from patients with Inflammatory Bowel Disease (IBD) and adherent-invasive E. coli have been found in resected ileum from patients with Crohn's disesae. This study aimed to characterize possible differences in phylogenetic group (triplex PCR), extraintestinal pathogenic E. coli (ExPEC) genes and multilocus sequence type (MLST) between E. coli strains isolated from IBD patients with past or present involvement of the left side of the colon and from controls.