Whole-genome Multilocus Sequence Typing Research Articles

Whole-Genome Multilocus Sequence Typing of Extended-Spectrum-Beta-Lactamase-Producing Enterobacteriaceae.

Molecular typing has become indispensable in the detection of nosocomial transmission of bacterial pathogens and the identification of sources and routes of transmission in outbreak settings, but current methods are labor-intensive, are difficult to standardize, or have limited resolution. Whole-genome multilocus sequence typing (wgMLST) has emerged as a whole-genome sequencing (WGS)-based gene-by-gene typing method that may overcome these limitations and has been applied successfully for several species in outbreak settings. In this study, genus-, genetic-complex-, and species-specific wgMLST schemes were developed for Citrobacter spp., the Enterobacter cloacae complex, Escherichia coli, Klebsiella oxytoca, and Klebsiella pneumoniae and used to type a national collection of 1,798 extended-spectrum-beta-lactamase-producing Enterobacteriaceae (ESBL-E) isolates obtained from patients in Dutch hospitals. Genus-, genetic-complex-, and species-specific thresholds for genetic distance that accurately distinguish between epidemiologically related and unrelated isolates were defined for Citrobacter spp., the E. cloacae complex, E. coli, and K. pneumoniae wgMLST was shown to have higher discriminatory power and typeability than in silico MLST. In conclusion, the wgMLST schemes developed in this study facilitate high-resolution WGS-based typing of the most prevalent ESBL-producing species in clinical practice and may contribute to further elucidation of the complex epidemiology of antimicrobial-resistant Enterobacteriaceae wgMLST opens up possibilities for the creation of a Web-accessible database for the global surveillance of ESBL-producing bacterial clones.

Virulence Gene Profiles and Clonal Relationships of Escherichia coli O26:H11 Isolates from Feedlot Cattle as Determined by Whole-Genome Sequencing

Escherichia coli O26 is the second most important enterohemorrhagic E. coli (EHEC) serogroup worldwide. Serogroup O26 strains are categorized mainly into two groups: enteropathogenic (EPEC) O26, carrying a locus of enterocyte effacement (LEE) and mostly causing mild diarrhea, and Shiga-toxigenic (STEC) O26, which carries the Shiga toxin (STX) gene (stx), responsible for more severe outcomes. stx-negative O26 strains can be further split into two groups. One O26 group differs significantly from O26 EHEC, while the other O26 EHEC-like group shows all the characteristics of EHEC O26 except production of STX. In order to determine the different populations of O26 E. coli present in U.S. cattle, we sequenced 42 O26:H11 strains isolated from feedlot cattle and compared them to 37 O26:H11 genomes available in GenBank. Phylogenetic analysis by whole-genome multilocus sequence typing (wgMLST) showed that O26:H11/H(-) strains in U.S. cattle were highly diverse. Most strains were sequence type 29 (ST29). By wgMLST, two clear lineages could be distinguished among cattle strains. Lineage 1 consisted of O26:H11 EHEC-like strains (ST29) (4 strains) and O26:H11 EHEC strains (ST21) (2 strains), and lineage 2 (36 strains) consisted of O26:H11 EPEC strains (ST29). Overall, our analysis showed U.S. cattle carried pathogenic (ST21; stx1 (+) ehxA(+) toxB(+)) and also potentially pathogenic (ST29; ehxA(+) toxB(+)) O26:H11 E. coli strains. Furthermore, in silico analysis showed that 70% of the cattle strains carried at least one antimicrobial resistance gene. Our results showed that whole-genome sequence analysis is a robust and valid approach to identify and genetically characterize E. coli O26:H11, which is of importance for food safety and public health. Escherichia coli O26 is the second most important type of enterohemorrhagic E. coli (EHEC) worldwide. Serogroup O26 strains are categorized into two groups: enteropathogenic (EPEC) carrying LEE, causing mild diarrhea, and Shiga toxigenic (STEC) carrying the stx gene, responsible for more severe outcomes. However, there are currently problems in distinguishing one group from the other. Furthermore, several O26 stx-negative strains are consistently misidentified as either EHEC-like or EPEC. The use of whole-genome sequence (WGS) analysis of O26 strains from cattle in the United States (i) allowed identification of O26 strains present in U.S. cattle, (ii) determined O26 strain diversity, (iii) solved the misidentification problem, and (iv) screened for the presence of antimicrobial resistance and virulence genes in the strains. This study provided a framework showing how to easily and rapidly use WGS information to identify and genetically characterize E. coli O26:H11, which is important for food safety and public health.

Implementation of Nationwide Real-time Whole-genome Sequencing to Enhance Listeriosis Outbreak Detection and Investigation.

Listeria monocytogenes (Lm) causes severe foodborne illness (listeriosis). Previous molecular subtyping methods, such as pulsed-field gel electrophoresis (PFGE), were critical in detecting outbreaks that led to food safety improvements and declining incidence, but PFGE provides limited genetic resolution. A multiagency collaboration began performing real-time, whole-genome sequencing (WGS) on all US Lm isolates from patients, food, and the environment in September 2013, posting sequencing data into a public repository. Compared with the year before the project began, WGS, combined with epidemiologic and product trace-back data, detected more listeriosis clusters and solved more outbreaks (2 outbreaks in pre-WGS year, 5 in WGS year 1, and 9 in year 2). Whole-genome multilocus sequence typing and single nucleotide polymorphism analyses provided equivalent phylogenetic relationships relevant to investigations; results were most useful when interpreted in context of epidemiological data. WGS has transformed listeriosis outbreak surveillance and is being implemented for other foodborne pathogens.

Tracing isolates from domestic human Campylobacter jejuni infections to chicken slaughter batches and swimming water using whole-genome multilocus sequence typing

Campylobacter jejuni is the leading cause of bacterial gastroenteritis and chicken is considered a major reservoir and source of human campylobacteriosis. In this study, we investigated temporally related Finnish human (n=95), chicken (n=83) and swimming water (n=20) C. jejuni isolates collected during the seasonal peak in 2012 using multilocus sequence typing (MLST) and whole-genome MLST (wgMLST). Our objective was to trace domestic human C. jejuni infections to C. jejuni isolates from chicken slaughter batches and swimming water. At MLST level, 79% of the sequence types (STs) of the human isolates overlapped with chicken STs suggesting chicken as an important reservoir. Four STs, the ST-45, ST-230, ST-267 and ST-677, covered 75% of the human and 64% of the chicken isolates. In addition, 50% of the swimming water isolates comprised ST-45, ST-230 and ST-677. Further wgMLST analysis of the isolates within STs, accounting their temporal relationship, revealed that 22 of the human isolates (24%) were traceable back to C. jejuni positive chicken slaughter batches. None of the human isolates were traced back to swimming water, which was rather sporadically sampled. The highly discriminatory wgMLST, together with the patient background information and temporal relationship data with possible sources, offers a new, accurate approach to trace back the origin of domestic campylobacteriosis. Our results suggest that potentially a substantial proportion of campylobacteriosis cases during the seasonal peak most probably are due to other sources than chicken meat consumption. These findings warrant further wgMLST-based studies to reassess the role of other reservoirs in the Campylobacter epidemiology both in Finland and elsewhere.

A Nonautochthonous U.S. Strain of Vibrio parahaemolyticus Isolated from Chesapeake Bay Oysters Caused the Outbreak in Maryland in 2010.

In the summer of 2010, Vibrio parahaemolyticus caused an outbreak in Maryland linked to the consumption of oysters. Strains isolated from both stool and oyster samples were indistinguishable by pulsed-field gel electrophoresis (PFGE). However, the oysters contained other potentially pathogenic V. parahaemolyticus strains exhibiting different PFGE patterns. In order to assess the identity, genetic makeup, relatedness, and potential pathogenicity of the V. parahaemolyticus strains, we sequenced 11 such strains (2 clinical strains and 9 oyster strains). We analyzed these genomes by in silico multilocus sequence typing (MLST) and determined their phylogeny using a whole-genome MLST (wgMLST) analysis. Our in silico MLST analysis identified six different sequence types (STs) (ST8, ST676, ST810, ST811, ST34, and ST768), with both of the clinical and four of the oyster strains being identified as belonging to ST8. Using wgMLST, we showed that the ST8 strains from clinical and oyster samples were nearly indistinguishable and belonged to the same outbreak, confirming that local oysters were the source of the infections. The remaining oyster strains were genetically diverse, differing in >3,000 loci from the Maryland ST8 strains. eBURST analysis comparing these strains with strains of other STs available at the V. parahaemolyticus MLST website showed that the Maryland ST8 strains belonged to a clonal complex endemic to Asia. This indicates that the ST8 isolates from clinical and oyster sources were likely not endemic to Maryland. Finally, this study demonstrates the utility of whole-genome sequencing (WGS) and associated analyses for source-tracking investigations. Vibrio parahaemolyticus is an important foodborne pathogen and the leading cause of bacterial infections in the United States associated with the consumption of seafood. In the summer of 2010, Vibrio parahaemolyticus caused an outbreak in Maryland linked to oyster consumption. Strains isolated from stool and oyster samples were indistinguishable by pulsed-field gel electrophoresis (PFGE). The oysters also contained other potentially pathogenic V. parahaemolyticus strains with different PFGE patterns. Since their identity, genetic makeup, relatedness, and potential pathogenicity were unknown, their genomes were determined by using next-generation sequencing. Whole-genome sequencing (WGS) analysis by whole-genome multilocus sequence typing (wgMLST) allowed (i) identification of clinical and oyster strains with matching PFGE profiles as belonging to ST8, (ii) determination of oyster strain diversity, and (iii) identification of the clinical strains as belonging to a clonal complex (CC) described only in Asia. Finally, WGS and associated analyses demonstrated their utility for trace-back investigations.

Comparative genome analysis of Weissella ceti, an emerging pathogen of farm-raised rainbow trout.

BackgroundThe genus Weissella belongs to the lactic acid bacteria and includes 18 currently identified species, predominantly isolated from fermented food but rarely from cases of bacteremia in animals. Recently, a new species, designated Weissella ceti, has been correlated with hemorrhagic illness in farm-raised rainbow trout in China, Brazil, and the USA, with high transmission and mortality rates during outbreaks. Although W. ceti is an important emerging veterinary pathogen, little is known about its genomic features or virulence mechanisms. To better understand these and to characterize the species, we have previously sequenced the genomes of W. ceti strains WS08, WS74, and WS105, isolated from different rainbow trout farms in Brazil and displaying different pulsed-field gel electrophoresis patterns. Here, we present a comparative analysis of the three previously sequenced genomes of W. ceti strains from Brazil along with W. ceti NC36 from the USA and those of other Weissella species.ResultsPhylogenomic and orthology-based analyses both showed a high-similarity in the genetic structure of these W. ceti strains. This structure is corroborated by the highly syntenic order of their genes and the neutral evolution inferred from Tajima’s D. A whole-genome multilocus sequence typing analysis distinguished strains WS08 and NC36 from strains WS74 and WS105. We predicted 10 putative genomic islands (GEI), among which PAIs 3a and 3b are phage sequences that occur only in WS105 and WS74, respectively, whereas PAI 1 is species specific.ConclusionsWe identified several genes putatively involved in the basic processes of bacterial physiology and pathogenesis, including survival in aquatic environment, adherence in the host, spread inside the host, resistance to immune-system-mediated stresses, and antibiotic resistance. These data provide new insights in the molecular epidemiology and host adaptation for this emerging pathogen in aquaculture.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2324-4) contains supplementary material, which is available to authorized users.

Population Genetics and Antimicrobial Susceptibility of Canine Campylobacter Isolates Collected before and after a Raw Feeding Experiment.

In recent years, increasing numbers of consumers have become interested in feeding raw food for their pet dogs as opposed to commercial dry food, in the belief of health advantages. However, raw meat and internal organs, possibly contaminated by pathogens such as Campylobacter spp., may pose a risk of transmission of zoonoses to the pet owners. Campylobacter jejuni is the leading cause of bacterial gastroenteritis in humans but C. upsaliensis has also been associated with human disease. In this study we investigated the effect of different feeding strategies on the prevalence of Campylobacter spp. in Finnish dogs. We further characterized the isolates using multilocus sequence typing (MLST), whole-genome (wg) MLST and antimicrobial susceptibility testing. Dogs were sampled before and after a feeding period consisting of commercial raw feed or dry pellet feed. Altogether 56% (20/36) of the dogs yielded at least one Campylobacter-positive fecal sample. C. upsaliensis was the major species detected from 39% of the dogs before and 30% after the feeding period. Two C. jejuni isolates were recovered, both from raw-fed dogs after the dietary regimen. The isolates represented the same genotype (ST-1326), suggesting a common infection source. However, no statistically significant correlation was found between the feeding strategies and Campylobacter spp. carriage. The global genealogy of MLST types of dog and human C. upsaliensis isolates revealed weakly clonal population structure as most STs were widely dispersed. Major antimicrobial resistance among C. upsaliensis isolates was against streptomycin (STR MIC > 4mg/l). Apart from that, all isolates were highly susceptible against the antimicrobials tested. Mutations were found in the genes rpsL or rpsL and rsmG in streptomycin resistant isolates. In conclusion, increasing trend to feed dogs with raw meat warrants more studies to evaluate the risk associated with raw feeding of pets in transmission of zoonoses to humans.

Refinement of whole-genome multilocus sequence typing analysis by addressing gene paralogy.

We developed a user-friendly program, Genome Profiler (GeP), to refine whole-genome multilocus sequence typing analysis by addressing gene paralogy with conserved gene neighborhoods. In comparison to similar programs, GeP produced overall the best results in terms of accuracy and is thus a useful alternative to resolve relationships of bacterial isolates.

Evidence of land-sea transfer of the zoonotic pathogen Campylobacter to a wildlife marine sentinel species.

Environmental pollution often accompanies the expansion and urbanization of human populations where sewage and wastewaters commonly have an impact on the marine environments. Here, we explored the potential for faecal bacterial pathogens, of anthropic origin, to spread to marine wildlife in coastal areas. The common zoonotic bacterium Campylobacter was isolated from grey seals (Halichoerus grypus), an important sentinel species for environmental pollution, and compared to isolates from wild birds, agricultural sources and clinical samples to characterize possible transmission routes. Campylobacter jejuni was present in half of all grey seal pups sampled (24/50 dead and 46/90 live pups) in the breeding colony on the Isle of May (Scotland), where it was frequently associated with histological evidence of disease. Returning yearling animals (19/19) were negative for C.jejuni suggesting clearance of infection while away from the localized colony infection source. The genomes of 90 isolates from seals were sequenced and characterized using a whole-genome multilocus sequence typing (MLST) approach and compared to 192 published genomes from multiple sources using population genetic approaches and a probabilistic genetic attribution model to infer the source of infection from MLST data. The strong genotype-host association has enabled the application of source attribution models in epidemiological studies of human campylobacteriosis, and here assignment analyses consistently grouped seal isolates with those from human clinical samples. These findings are consistent with either a common infection source or direct transmission of human campylobacter to grey seals, raising concerns about the spread of human pathogens to wildlife marine sentinel species in coastal areas.

Multilocus sequence typing (MLST) and whole-genome MLST of Campylobacter jejuni isolates from human infections in three districts during a seasonal peak in Finland.

A total of 95 human Campylobacter jejuni isolates acquired from domestic infections and collected from three districts in Finland during the seasonal peak (June to September) in 2012 were analyzed by PCR-based multilocus sequence typing (MLST) and by whole-genome sequencing (WGS). Four predominant sequence types (STs) were detected among the isolates: ST-45 (21%) and ST-230 (14%, ST-45 clonal complex [CC]), ST-267 (21%, ST-283 CC), and ST-677 (19%, ST-677 CC). In districts 1 and 3, most of the infections occurred from early July to the middle of August, with a peak at weeks 29 to 31, but in district 2, the infections were dispersed more evenly throughout 3 months (June to August). WGS data were used for further whole-genome MLST (wgMLST) analyses of the isolates representing the four common STs. Shared loci of the isolates within each ST were analyzed as distance matrices of allelic profiles by the neighbor-net algorithm. The highest allelic variations (>400 different alleles) were detected between different clusters of ST-45 isolates (1,121 shared loci), while ST-230 (1,264 shared loci), ST-677 (1,169 shared loci), and ST-267 isolates (1,217 shared loci) were less diverse with the clusters differing by <40 alleles. Closely related isolates showing no allelic variation (subclusters) were detected among all four major STs. In some cases, they originated from different districts, suggesting that isolates can be epidemiologically connected and may have the same infection source despite being originally identified as sporadic infections.

Using multilocus sequence typing to study bacterial variation: prospects in the genomic era.

Multilocus sequence typing (MLST) indexes the sequence variation present in a small number (usually seven) of housekeeping gene fragments located around the bacterial genome. Unique alleles at these loci are assigned arbitrary integer identifiers, which effectively summarizes the variation present in several thousand base pairs of genome sequence information as a series of numbers. Comparing bacterial isolates using allele-based methods efficiently corrects for the effects of lateral gene transfer present in many bacterial populations and is computationally efficient. This 'gene-by-gene' approach can be applied to larger collections of loci, such as the ribosomal protein genes used in ribosomal MLST (rMLST), up to and including the complete set of coding sequences present in a genome, whole-genome MLST (wgMLST), providing scalable, efficient and readily interpreted genome analysis.