Rapid and high resolution genotyping of all Escherichia coli serotypes using 10 genomic repeat-containing loci

Our laboratory has previously published two multiple-locus variable-number tandem-repeats analysis (MLVA) methods for rapid genotyping of Escherichia coli (E. coli), which are now in routine use for surveillance and outbreak detection. The first assay developed was specific for E. coli O157:H7; however this assay was not suitable for genotyping other E. coli serotypes. A new generic MLVA-assay was then developed with the capability of genotyping all E. coli serotypes. This generic E. coli MLVA (GECM7) was based on polymorphism in seven variable number of tandem repeats (VNTR) loci. GECM7 worked well with the majority of E. coli serotypes; however we wanted to increase the resolution for this method based in part of comparison with PFGE typing of E. coli O26:H11, where PFGE appeared to display higher resolution. The GECM7 method was improved by adding three new repeat-loci to a total of ten (GECM10), and a considerable increase in resolution was observed (from 296 to 507 genotypes on the same set of strains).     

Genotyping of Bifidobacterium longum subsp. longum strains by multilocus variable number of tandem repeat analysis

A multilocus variable number of tandem repeat analysis (MLVA) scheme was developed and 44 isolates of Bifidobacterium longum subsp. longum were typed, including 5 isolates recovered during a clinical trial. The MLVA scheme generated 19 profiles and proved to be a fast, reliable and relatively cheap typing method.     

Multi Locus Variable-Number Tandem Repeat (MLVA) Typing Tools Improved the Surveillance of Salmonella Enteritidis: A 6 Years Retrospective Study

Surveillance of Salmonella enterica subsp. enterica serovar Enteritidis is generally considered to benefit from molecular techniques like multiple-locus variable-number of tandem repeats analysis (MLVA), which allow early detection and confinement of outbreaks. Here, a surveillance study, including phage typing, antimicrobial susceptibility testing and MLVA on 1,535 S. Enteritidis isolates collected between 2007 and 2012, was used to evaluate the added value of MLVA for public health surveillance in Belgium. Phage types PT4, PT8, PT21, PT1, PT6, PT14b, PT28 and PT13 dominate the Belgian S. Enteritidis population. The isolates of S. Enteritidis were most frequently susceptible to all antibiotics tested. 172 different MLVA profiles were detected, of which 9 frequent profiles included 67.2% of the S. Enteritidis population. During a serial passage experiment on selected isolates to investigate the in vitro stability of the 5 MLVA loci, no variations over time were observed indicating that the MLVA profiles were stable. The MLVA profile of isolates originating from different outbreaks in the Democratic Republic of the Congo (DRC) between 2010 and 2011 were distinct from any of the MLVA profiles found in Belgian isolates throughout the six year observational period and demonstrates that MLVA improves public health surveillance of S. Enteritidis. However, MLVA should be complemented with other subtyping methods when investigating outbreaks is caused by the most common MLVA profile.     

A New Methodology for Differentiation and Typing of Closely Related Salmonella enterica serovar Heidelberg Isolates

This study describe the use of a combination of two recently proposed typing approaches, multiple amplification of prophage locus typing (MAPLT) and multiple-locus variable-number tandem-repeat analysis (MLVA) for subdividing within Salmonella enterica serovar Heidelberg (S. Heidelberg). The combined typing method was compared with pulsed-field gel electrophoresis (PFGE) by Simpson’s index of diversity (DI). PFGE was shown to have a DI?=?0.84 and was poor at differentiation of the predominant PT1 (Phage Type 1) phenotype. In comparison, the combined MAPLT/MLVA method comprising 3 MLVA and 9 MAPLT primer pairs provided a higher differentiating ability DI?=?0.92. More importantly, the combined methodology was found to be superior in the differentiation of the predominant PT1 isolates. In conclusion, this study demonstrated the potential of the rapid and simple amalgamated MAPLT/MLVA approach in determining transmission of isolates of clonal phage type groups from various environmental sources to humans.     

Epidemiology of a Salmonella enterica subsp. enterica Serovar Typhimurium Strain Associated with a Songbird Outbreak

Salmonella enterica subsp. enterica serovar Typhimurium is responsible for the majority of salmonellosis cases worldwide. This Salmonella serovar is also responsible for die-offs in songbird populations. In 2009, there was an S. Typhimurium epizootic reported in pine siskins in the eastern United States. At the time, there was also a human outbreak with this serovar that was associated with contaminated peanuts. As peanuts are also used in wild-bird food, it was hypothesized that the pine siskin epizootic was related to this human outbreak. A comparison of songbird and human S. Typhimurium pulsed-field gel electrophoresis (PFGE) patterns revealed that the epizootic was attributed not to the peanut-associated strain but, rather, to a songbird strain first characterized from an American goldfinch in 1998. This same S. Typhimurium strain (PFGE type A3) was also identified in the PulseNet USA database, accounting for 137 of 77,941 total S. Typhimurium PFGE entries. A second molecular typing method, multiple-locus variable-number tandem-repeat analysis (MLVA), confirmed that the same strain was responsible for the pine siskin epizootic in the eastern United States but was distinct from a genetically related strain isolated from pine siskins in Minnesota. The pine siskin A3 strain was first encountered in May 2008 in an American goldfinch and later in a northern cardinal at the start of the pine siskin epizootic. MLVA also confirmed the clonal nature of S. Typhimurium in songbirds and established that the pine siskin epizootic strain was unique to the finch family. For 2009, the distribution of PFGE type A3 in passerines and humans mirrored the highest population density of pine siskins for the East Coast.     

MLVA as a Tool for Public Health Surveillance of Human Salmonella Typhimurium: Prospective Study in Belgium and Evaluation of MLVA Loci Stability

Surveillance of Salmonella enterica subsp. enterica serovar Typhimurium (S. Typhimurium) is generally considered to benefit from molecular techniques like multiple-locus variable-number of tandem repeats analysis (MLVA), which allow earlier detection and confinement of outbreaks. Here, a surveillance study, including phage typing, antimicrobial susceptibility testing and the in Europe most commonly used 5-loci MLVA on 1,420 S. Typhimurium isolates collected between 2010 and 2012 in Belgium, was used to evaluate the added value of MLVA for public health surveillance. Phage types DT193, DT195, DT120, DT104, DT12 and U302 dominate the Belgian S. Typhimurium population. A combined resistance to ampicillin, streptomycin, sulphonamides and tetracycline (ASSuT) with or without additional resistances was observed for 42.5% of the isolates. 414 different MLVA profiles were detected, of which 14 frequent profiles included 44.4% of the S. Typhimurium population. During a serial passage experiment on selected isolates to investigate the in vitro stability of the 5 MLVA loci, variations over time were observed for loci STTR6, STTR10, STTR5 and STTR9. This study demonstrates that MLVA improves public health surveillance of S. Typhimurium. However, the 5-loci MLVA should be complemented with other subtyping methods for investigation of possible outbreaks with frequent MLVA profiles. Also, variability in these MLVA loci should be taken into account when investigating extended outbreaks and studying dynamics over longer periods.     

Detection and Identification of Salmonella enterica,Escherichia coli,and Shigella spp. via PCR-Electrospray Ionization Mass Spectrometry: Isolate Testing and Analysis of Food Samples

An assay to identify the common food-borne pathogens Salmonella, Escherichia coli, Shigella, and Listeria monocytogenes was developed in collaboration with Ibis Biosciences (a division of Abbott Molecular) for the Plex-ID biosensor system, a platform that uses electrospray ionization mass spectroscopy (ESI-MS) to detect the base composition of short PCR amplicons. The new food-borne pathogen (FBP) plate has been experimentally designed using four gene segments for a total of eight amplicon targets. Initial work built a DNA base count database that contains more than 140 Salmonella enterica, 139 E. coli, 11 Shigella, and 36 Listeria patterns and 18 other Enterobacteriaceae organisms. This assay was tested to determine the scope of the assay''s ability to detect and differentiate the enteric pathogens and to improve the reference database associated with the assay. More than 800 bacterial isolates of S. enterica, E. coli, and Shigella species were analyzed. Overall, 100% of S. enterica, 99% of E. coli, and 73% of Shigella spp. were detected using this assay. The assay was also able to identify 30% of the S. enterica serovars to the serovar level. To further characterize the assay, spiked food matrices and food samples collected during regulatory field work were also studied. While analysis of preenrichment media was inconsistent, identification of S. enterica from selective enrichment media resulted in serovar-level identifications for 8 of 10 regulatory samples. The results of this study suggest that this high-throughput method may be useful in clinical and regulatory laboratories testing for these pathogens.     

Comparison of a semi-automated rep-PCR system and multilocus sequence typing for differentiation of Salmonella enterica isolates

The accurate sub-typing of Salmonella enterica isolates is essential for epidemiological investigations and surveillance of Salmonella infections. Salmonella isolates are currently identified using the Kauffman-White serotyping scheme. Multilocus sequence typing (MLST) schemes have been developed for the major bacterial pathogens including Salmonella and have assisted in understanding the molecular epidemiology and population biology of these organisms. Recently, the DiversiLab rep-PCR system has been developed using micro-fluidic chips to provide standardized, semi-automated fingerprinting for pathogens including S. enterica. In the current study, 71 isolates of S. enterica, representing 21 serovars, were analyzed using MLST and the DiversiLab rep-PCR system. MLST was able to identify 31 sequence types (STs), while the DiversiLab system revealed 38 DiversiLab types (DTs). The rep-PCR distinguished isolates of different serovars and showed greater discriminatory power (0.95) than MLST typing (0.89). Rep-PCR exhibited 92% concordance with MLST and 90% with serotyping, while the concordance level of MLST typing with serotyping was 96%, representing a strong association. Comparison of rep-PCR profiles with those held in an online library database led to the accurate prediction of serovar in 63% of cases and resulted in inaccurate predictions for 10% of profiles. MLST and the rep-PCR system may provide useful additional informative techniques for the molecular identification of S. enterica. We conclude that the DiversiLab rep-PCR system may provide a rapid (less than 4 h) and standardized method for sub-typing isolates of S. enterica.     

Genotyping of Bifidobacterium longum subsp. longum strains by multilocus variable number of tandem repeat analysis

A multilocus variable number of tandem repeat analysis (MLVA) scheme was developed and 44 isolates of Bifidobacterium longum subsp. longum were typed, including 5 isolates recovered during a clinical trial. The MLVA scheme generated 19 profiles and proved to be a fast, reliable and relatively cheap typing method.     

Evolution of Pan-Genomes of Escherichia coli,Shigella spp., and Salmonella enterica

Multiple sequencing of genomes belonging to a bacterial species allows one to analyze and compare statistics and dynamics of the gene complements of species, their pan-genomes. Here, we analyzed multiple genomes of Escherichia coli, Shigella spp., and Salmonella enterica. We demonstrate that the distribution of the number of genomes harboring a gene is well approximated by a sum of two power functions, describing frequent genes (present in many strains) and rare genes (present in few strains). The virtual absence of Shigella-specific genes not present in E. coli genomes confirms previous observations that Shigella is not an independent genus. While the pan-genome size is increasing with each new strain, the number of genes present in a fixed fraction of strains stabilizes quickly. For instance, slightly fewer than 4,000 genes are present in at least half of any group of E. coli genomes. Comparison of S. enterica and E. coli pan-genomes revealed the existence of a common periphery, that is, genes present in some but not all strains of both species. Analysis of phylogenetic trees demonstrates that rare genes from the periphery likely evolve under horizontal transfer, whereas frequent periphery genes may have been inherited from the periphery genome of the common ancestor.     

Development of Multiple-Locus Variable-Number Tandem-Repeat Analysis for Molecular Subtyping of Campylobacter jejuni by Using Capillary Electrophoresis

Campylobacter jejuni is a common cause of the frequently reported food-borne diseases in developed and developing nations. This study describes the development of multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) using capillary electrophoresis as a novel typing method for microbial source tracking and epidemiological investigation of C. jejuni. Among 36 tandem repeat loci detected by the Tandem Repeat Finder program, 7 VNTR loci were selected and used for characterizing 60 isolates recovered from chicken meat samples from retail shops, samples from chicken meat processing factory, and stool samples. The discrimination ability of MLVA was compared with that of multilocus sequence typing (MLST). MLVA (diversity index of 0.97 with 31 MLVA types) provided slightly higher discrimination than MLST (diversity index of 0.95 with 25 MLST types). The overall concordance between MLVA and MLST was estimated at 63% by adjusted Rand coefficient. MLVA predicted MLST type better than MLST predicted MLVA type, as reflected by Wallace coefficient (Wallace coefficient for MLVA to MLST versus MLST to MLVA, 86% versus 51%). MLVA is a useful tool and can be used for effective monitoring of C. jejuni and investigation of epidemics caused by C. jejuni.     

Evaluation of Multiple-Locus Variable Number of Tandem Repeats Analysis for Typing Livestock-Associated Methicillin-Resistant Staphylococcus aureus

Background

The increasing occurrence of livestock-associated (LA) methicillin-resistant Staphylococcus aureus (MRSA) associated with the clonal complex (CC) 398 within the past years shows the importance of standardized and comparable typing methods for the purposes of molecular surveillance and outbreak detection. Multiple-locus variable number of tandem repeats analysis (MLVA) has recently been described as an alternative and highly discriminative tool for S. aureus. However, until now the applicability of MLVA for the typing of LA-MRSA isolates from different geographic origin has not been investigated in detail. We therefore compared MLVA and S. aureus protein A (spa) typing for characterizing porcine MRSA from distinct Dutch and German farms.

Methodology/Principal Findings

Overall, 134 MRSA isolates originating from 21 different pig-farms in the Netherlands and 36 farms in Germany comprising 21 different spa types were subjected to MLVA-typing. Amplification and subsequent automated fragment sizing of the tandem repeat loci on a capillary sequencer differentiated these 134 isolates into 20 distinct MLVA types. Whereas overall MLVA and spa typing showed the same discriminatory power to type LA-MRSA (p  = 0.102), MLVA was more discriminatory than spa typing for isolates associated with the prevalent spa types t011 and t034 (Simpson’s Index of Diversity 0.564 vs. 0.429, respectively; p<0.001).

Conclusion

Although the applied MLVA scheme was not more discriminatory than spa typing in general, it added valuable information to spa typing results for specific spa types (t011, t034) which are highly prevalent in the study area, i.e. Dutch-German border area. Thus, both methods may complement each other to increase the discriminatory power to resolute highly conserved clones such as CC398 (spa types t011, t034) for the detection of outbreaks and molecular surveillance of zoonotic MRSA.     

Subtyping of a Large Collection of Historical Listeria monocytogenes Strains from Ontario,Canada, by an Improved Multilocus Variable-Number Tandem-Repeat Analysis (MLVA)

Listeria monocytogenes is responsible for severe and often fatal food-borne infections in humans. A collection of 2,421 L. monocytogenes isolates originating from Ontario''s food chain between 1993 and 2010, along with Ontario clinical isolates collected from 2004 to 2010, was characterized using an improved multilocus variable-number tandem-repeat analysis (MLVA). The MLVA method was established based on eight primer pairs targeting seven variable-number tandem-repeat (VNTR) loci in two 4-plex fluorescent PCRs. Diversity indices and amplification rates of the individual VNTR loci ranged from 0.38 to 0.92 and from 0.64 to 0.99, respectively. MLVA types and pulsed-field gel electrophoresis (PFGE) patterns were compared using Comparative Partitions analysis involving 336 clinical and 99 food and environmental isolates. The analysis yielded Simpson''s diversity index values of 0.998 and 0.992 for MLVA and PFGE, respectively, and adjusted Wallace coefficients of 0.318 when MLVA was used as a primary subtyping method and 0.088 when PFGE was a primary typing method. Statistical data analysis using BioNumerics allowed for identification of at least 8 predominant and persistent L. monocytogenes MLVA types in Ontario''s food chain. The MLVA method correctly clustered epidemiologically related outbreak strains and separated unrelated strains in a subset analysis. An MLVA database was established for the 2,421 L. monocytogenes isolates, which allows for comparison of data among historical and new isolates of different sources. The subtyping method coupled with the MLVA database will help in effective monitoring/prevention approaches to identify environmental contamination by pathogenic strains of L. monocytogenes and investigation of outbreaks.     

Fine-Scale Genetic Diversity among Burkholderia pseudomallei Soil Isolates in Northeast Thailand

Burkholderia pseudomallei soil isolates from northeast Thailand were genotyped using multiple-locus variable-number tandem repeat (VNTR) analysis (MLVA) and multilocus sequence typing (MLST). MLVA identified 19 genotypes within three clades, while MLST revealed two genotypes. These close genetic relationships imply a recent colonization followed by localized expansion, similar to what occurs in an outbreak situation.     

An evaluation of the species and subspecies of the genus Salmonella with whole genome sequence data: Proposal of type strains and epithets for novel S. enterica subspecies VII,VIII, IX,X and XI

Species and subspecies within the Salmonella genus have been defined for public health purposes by biochemical properties; however, reference laboratories have increasingly adopted sequence-based, and especially whole genome sequence (WGS), methods for surveillance and routine identification. This leads to potential disparities in subspecies definitions, routine typing, and the ability to detect novel subspecies. A large-scale analysis of WGS data from the routine sequencing of clinical isolates was employed to define and characterise Salmonella subspecies population structure, demonstrating that the Salmonella species and subspecies were genetically distinct, including those previously identified through phylogenetic approaches, namely: S. enterica subspecies londinensis (VII), subspecies brasiliensis (VIII), subspecies hibernicus (IX) and subspecies essexiensis (X). The analysis also identified an additional novel subspecies, reptilium (XI). Further, these analyses indicated that S. enterica subspecies arizonae (IIIa) isolates were divergent from the other S. enterica subspecies, which clustered together and, on the basis of ANI analysis, subspecies IIIa was sufficiently distinct to be classified as a separate species, S. arizonae. Multiple phylogenetic and statistical approaches generated congruent results, suggesting that the proposed species and subspecies structure was sufficiently biologically robust for routine application. Biochemical analyses demonstrated that not all subspecies were distinguishable by these means and that biochemical approaches did not capture the genomic diversity of the genus. We recommend the adoption of standardised genomic definitions of species and subspecies and a genome sequence-based approach to routine typing for the identification and definition of novel subspecies.     

MLVA Genotyping of Brucella melitensis and Brucella abortus Isolates from Different Animal Species and Humans and Identification of Brucella suis Vaccine Strain S2 from Cattle in China

In China, brucellosis is an endemic disease and the main sources of brucellosis in animals and humans are infected sheep, cattle and swine. Brucella melitensis (biovars 1 and 3) is the predominant species, associated with sporadic cases and outbreak in humans. Isolates of B. abortus, primarily biovars 1 and 3, and B. suis biovars 1 and 3 are also associated with sporadic human brucellosis. In this study, the genetic profiles of B. melitensis and B. abortus isolates from humans and animals were analyzed and compared by multi-locus variable-number tandem-repeat analysis (MLVA). Among the B. melitensis isolates, the majority (74/82) belonged to MLVA8 genotype 42, clustering in the ‘East Mediterranean’ group. Two B. melitensis biovar 1 genotype 47 isolates, belonging to the ‘Americas’ group, were recovered; both were from the Himalayan blue sheep (Pseudois nayaur, a wild animal). The majority of B. abortus isolates (51/70) were biovar 3, genotype 36. Ten B. suis biovar 1 field isolates, including seven outbreak isolates recovered from a cattle farm in Inner Mongolia, were genetically indistinguishable from the vaccine strain S2, based on MLVA cluster analysis. MLVA analysis provided important information for epidemiological trace-back. To the best of our knowledge, this is the first report to associate Brucella cross-infection with the vaccine strain S2 based on molecular comparison of recovered isolates to the vaccine strain. MLVA typing could be an essential assay to improve brucellosis surveillance and control programs.     

A new typing technique for the Rickettsiales Ehrlichia ruminantium: multiple-locus variable number tandem repeat analysis

Ehrlichia ruminantium (ER) is a member of the order Rickettsiales transmitted by Amblyomma ticks. This obligatory intracellular bacterium is the causative agent of a fatal disease in ruminants, named heartwater. It represents a constraint on breeding development in sub-Saharan Africa and in the Caribbean. The genetic diversity of the strains of ER, which could be a limiting factor to obtain effective vaccines, needs to be better characterized. For this purpose, we developed a molecular typing technique based on the polymorphism of variable number tandem repeat (VNTR) sequences, MLVA (multiple locus VNTR analysis).Eight (out of 21) VNTR candidates were validated using 17 samples representing a panel of ER strains from different geographical origins from West, South Africa, and Caribbean areas and in ER infected ticks and goat tissues. This result demonstrated the ability of these VNTRs to type a wide range of strains. The stability of the selected VNTR markers was very good, at the time scale needed for epidemiological purposes: in particular, no difference in the VNTR profiles was observed between virulent and attenuated strains (for Gardel and Senegal strains) and between strains (Gardel and Blonde strains) isolated in the same area 19 years apart. We validated the strong discriminatory power of MLVA for ER and found a high level of polymorphism between the available strains, with 10 different profiles out of 13 ER strains.The MLVA scheme described in this study is a rapid and efficient molecular typing tool for ER, which allows rapid and direct typing of this intracellular pathogen without preliminary culture and gives reliable results that can be used for further epidemiological studies.     

Genotypic and PFGE/MLVA analyses of Vibrio cholerae O1: geographical spread and temporal changes during the 2007-2010 cholera outbreaks in Thailand

Background

Vibrio cholerae O1 El Tor dominated the seventh cholera pandemic which occurred in the 1960s. For two decades, variants of V. cholerae O1 El Tor that produce classical cholera toxin have emerged and spread globally, replacing the prototypic El Tor biotype. This study aims to characterize V. cholerae O1 isolates from outbreaks in Thailand with special reference to genotypic variations over time.

Methods/Findings

A total of 343 isolates of V. cholerae O1 from cholera outbreaks from 2007 to 2010 were investigated, and 99.4% were found to carry the classical cholera toxin B subunit (ctxB) and El Tor rstR genes. Pulsed-field gel electrophoresis (PFGE) differentiated the isolates into 10 distinct pulsotypes, clustered into two major groups, A and B, with an overall similarity of 88%. Ribotyping, multiple-locus variable-number tandem-repeat analysis (MLVA), and PCR to detect Vibrio seventh pandemic island II (VSP-II) related genes of randomly selected isolates from each pulsotype corresponded to the results obtained by PFGE. Epidemiological investigations revealed that MLVA type 2 was strongly associated with a cholera outbreak in northeastern Thailand in 2007, while MLVA type 7 dominated the outbreaks of the southern Gulf areas in 2009 and MLVA type 4 dominated the outbreaks of the central Gulf areas during 2009–2010. Only MLVA type 16 isolates were found in a Thai-Myanmar border area in 2010, whereas those of MLVA types 26, 39, and 41 predominated this border area in 2008. Type 39 then disappeared 1–2 years later as MLVA type 41 became prevalent. Type 41 was also found to infect an outbreak area.

Conclusions

MLVA provided a high-throughput genetic typing tool for understanding the in-depth epidemiology of cholera outbreaks. Our epidemiological surveys suggest that some clones of V. cholerae O1 with similar but distinctive genetic traits circulate in outbreak sites, while others disappear over time.     

Streptococcus suis outbreak investigation using multiple‐locus variable tandem repeat number analysis

Two outbreaks of Streptococcus suis ST7 occurred in humans in 1998 and 2005 in China. PFGE of chromosome restriction fragments found all ST7 isolates to be indistinguishable. Due to the genetic homogeneity of ST7 isolates, development of a rapid sub‐typing method with high discriminatory power for ST7 isolates is required. In this study, a novel method, MLVA, was developed to type S. suis serotype 2 strains. Further, this method was used to analyze outbreak‐associated ST7 strains in China. A total of 144 ST7 S. suis isolates were sub‐typed into 34 MLVA types. Among these, eight isolates from the 1998 outbreak were sub‐typed into five MLVA types, of which four MLVA types were also detected in Sichuan in 2005. These data indicate that the pathogens responsible for the two outbreaks had the same origin. In addition, some observations also provided molecular evidence for the transmission route, possibly indicating that the MLVA method has usefulness in epidemiology. The developed MLVA scheme for S. suis has greater discriminative power than PFGE. The method described here may be useful for identifying the source of S. suis infection and monitoring its spread.