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

Mini-Multilocus Sequence Typing Scheme for the Global Population of Neisseria gonorrhoeae

The increasing problem of antimicrobial resistance in N. gonorrhoeae necessitates the development of molecular typing schemes that are suitable for rapid and mass screening. The objective of this study was to design and validate a mini-MLST scheme for N. gonorrhoeae based on global pathogen population data. Using sequences of seven housekeeping genes of 21,402 isolates with known MLSTs from the PubMLST database, we identified eighteen informative polymorphisms and obtained mini-MLST nucleotide profiles to predict MLSTs of isolates. We proposed a new MLST grouping system for N. gonorrhoeae based on mini-MLST profiles. Phylogenetic analysis revealed that MLST genogroups are a stable characteristic of the N. gonorrhoeae global population. The proposed grouping system has been shown to bring together isolates with similar antimicrobial susceptibility, as demonstrated by the characteristics of major genogroups. Established MLST prediction algorithms based on nucleotide profiles are now publicly available. The mini-MLST scheme was evaluated using a MLST detection/prediction method based on the original hydrogel DNA microarray. The results confirmed a high predictive ability up to the MLST genogroup. The proposed holistic approach to gonococcal population analysis can be used for the continuous surveillance of known and emerging resistant N. gonorrhoeae isolates.

Rapid identification of lactic acid bacteria at species/subspecies level via ensemble learning of Ramanomes

Rapid and accurate identification of lactic acid bacteria (LAB) species would greatly improve the screening rate for functional LAB. Although many conventional and molecular methods have proven efficient and reliable, LAB identification using these methods has generally been slow and tedious. Single-cell Raman spectroscopy (SCRS) provides the phenotypic profile of a single cell and can be performed by Raman spectroscopy (which directly detects vibrations of chemical bonds through inelastic scattering by a laser light) using an individual live cell. Recently, owing to its affordability, non-invasiveness, and label-free features, the Ramanome has emerged as a potential technique for fast bacterial detection. Here, we established a reference Ramanome database consisting of SCRS data from 1,650 cells from nine LAB species/subspecies and conducted further analysis using machine learning approaches, which have high efficiency and accuracy. We chose the ensemble meta-classifier (EMC), which is suitable for solving multi-classification problems, to perform in-depth mining and analysis of the Ramanome data. To optimize the accuracy and efficiency of the machine learning algorithm, we compared nine classifiers: LDA, SVM, RF, XGBoost, KNN, PLS-DA, CNN, LSTM, and EMC. EMC achieved the highest average prediction accuracy of 97.3% for recognizing LAB at the species/subspecies level. In summary, Ramanomes, with the integration of EMC, have promising potential for fast LAB species/subspecies identification in laboratories and may thus be further developed and sharpened for the direct identification and prediction of LAB species from fermented food.

Inoviridae prophage and bacterial host dynamics during diversification, succession, and Atlantic invasion of Pacific-native Vibrio parahaemolyticus

IMPORTANCE

An understanding of the processes that contribute to the emergence of pathogens from environmental reservoirs is critical as changing climate precipitates pathogen evolution and population expansion. Phylogeographic analysis of Vibrio parahaemolyticus hosts combined with the analysis of their Inoviridae phage resolved ambiguities of diversification dynamics which preceded successful Atlantic invasion by the epidemiologically predominant ST36 lineage. It has been established experimentally that filamentous phage can limit host recombination, but here, we show that phage loss is linked to rapid bacterial host diversification during epidemic spread in natural ecosystems alluding to a potential role for ubiquitous inoviruses in the adaptability of pathogens. This work paves the way for functional analyses to define the contribution of inoviruses in the evolutionary dynamics of environmentally transmitted pathogens.

Core genome multilocus sequence typing analysis of Finnish Taylorella equigenitalis isolates

In Finland, Taylorella equigenitalis, the causative agent of contagious equine metritis (CEM), was first detected in 1992. The aim of this study was to genotype Finnish T. equigenitalis isolates to investigate the epidemiology of the infection in the Finnish horse population. A total of 34 T. equigenitalis isolates from 24 horses obtained during 1992-2021 were subjected to whole genome sequencing (WGS) and subsequent local ad hoc core genome multi-locus sequence typing (cgMLST) targeting 1259 loci. Classical MLST profiles were extracted from the whole-genome sequence data. Three novel MLST types, ST81, ST82 and ST83, and four previously described sequence types, ST16, ST17, ST50 and ST63 were detected among the isolates. cgMLST minimum spanning tree analysis using 12 allele difference as threshold, resulted in five clusters and three singletons. cgMLST clusters were congruent with the MLST-defined groups, except for the ST83 isolates which were divided into two clusters. However, the high discriminatory power cgMLST allowed differentiation between isolates of the same MLST type as each isolate had a unique core genome ST. Our study suggests that cgMLST has the prospective for being a standardised typing method for T. equigenitalis in the future, and further contributes to worldwide phylogenetic and spatio-temporal analyses needed to better understand the epidemiology of the bacterium.

Single-cell rapid identification, in situ viability and vitality profiling, and genome-based source-tracking for probiotics products

Rapid expansion of the probiotics industry demands fast, sensitive, comprehensive, and low-cost strategies for quality assessment. Here, we introduce a culture-free, one-cell-resolution, phenome-genome-combined strategy called Single-Cell Identification, Viability and Vitality tests, and Source-tracking (SCIVVS). For each cell directly extracted from the product, the fingerprint region of D2O-probed single-cell Raman spectrum (SCRS) enables species-level identification with 93% accuracy, based on a reference SCRS database from 21 statutory probiotic species, whereas the C-D band accurately quantifies viability, metabolic vitality plus their intercellular heterogeneity. For source-tracking, single-cell Raman-activated Cell Sorting and Sequencing can proceed, producing indexed, precisely one-cell-based genome assemblies that can reach ~99.40% genome-wide coverage. Finally, we validated an integrated SCIVVS workflow with automated SCRS acquisition where the whole process except sequencing takes just 5 h. As it is >20-fold faster, >10-time cheaper, vitality-revealing, heterogeneity-resolving, and automation-prone, SCIVVS is a new technological and data framework for quality assessment of live-cell products.

ReporTree: a surveillance-oriented tool to strengthen the linkage between pathogen genetic clusters and epidemiological data

BACKGROUND

Genomics-informed pathogen surveillance strengthens public health decision-making, playing an important role in infectious diseases' prevention and control. A pivotal outcome of genomics surveillance is the identification of pathogen genetic clusters and their characterization in terms of geotemporal spread or linkage to clinical and demographic data. This task often consists of the visual exploration of (large) phylogenetic trees and associated metadata, being time-consuming and difficult to reproduce.

RESULTS

We developed ReporTree, a flexible bioinformatics pipeline that allows diving into the complexity of pathogen diversity to rapidly identify genetic clusters at any (or all) distance threshold(s) or cluster stability regions and to generate surveillance-oriented reports based on the available metadata, such as timespan, geography, or vaccination/clinical status. ReporTree is able to maintain cluster nomenclature in subsequent analyses and to generate a nomenclature code combining cluster information at different hierarchical levels, thus facilitating the active surveillance of clusters of interest. By handling several input formats and clustering methods, ReporTree is applicable to multiple pathogens, constituting a flexible resource that can be smoothly deployed in routine surveillance bioinformatics workflows with negligible computational and time costs. This is demonstrated through a comprehensive benchmarking of (i) the cg/wgMLST workflow with large datasets of four foodborne bacterial pathogens and (ii) the alignment-based SNP workflow with a large dataset of Mycobacterium tuberculosis. To further validate this tool, we reproduced a previous large-scale study on Neisseria gonorrhoeae, demonstrating how ReporTree is able to rapidly identify the main species genogroups and characterize them with key surveillance metadata, such as antibiotic resistance data. By providing examples for SARS-CoV-2 and the foodborne bacterial pathogen Listeria monocytogenes, we show how this tool is currently a useful asset in genomics-informed routine surveillance and outbreak detection of a wide variety of species.

CONCLUSIONS

In summary, ReporTree is a pan-pathogen tool for automated and reproducible identification and characterization of genetic clusters that contributes to a sustainable and efficient public health genomics-informed pathogen surveillance. ReporTree is implemented in python 3.8 and is freely available at https://github.com/insapathogenomics/ReporTree .

Clonal diversity of Staphylococcus aureus isolates in clinical specimens from selected health facilities in Ethiopia

Staphylococcus aureus is among the top three causative agents of nosocomial infection in Ethiopia. The majority of studies in Ethiopia have focused on the epidemiology of S. aureus in hospital settings, with limited molecular genotyping results. Molecular characterization of S. aureus is essential for identification of strains, and contributes to the control and prevention of S. aureus infection. The aim of the current study was to determine the molecular epidemiology of methicillin-susceptible S. aureus (MSSA) and methicillin-resistant S. aureus (MRSA) isolates recovered from clinical specimens in Ethiopia. A total of 161 MSSA and 9 MRSA isolates were characterized using pulsed-field gel electrophoresis (PFGE) and staphylococcal protein A (spa) typing. Based on the PFGE analysis, MSSA isolates were grouped into eight pulso-types groups (from A to I), while MRSA isolates clustered into three (A, B and C) pulso-types with more than 80% similarity. The spa typing analysis showed diversity of S. aureus with 56 distinct spa types. Spa type t355 was most prevalent (56/170, 32.9%), while eleven new spa types were detected including t20038, t20039, and t20042. The identified spa types were clustered into 15 spa-clonal complexes (spa-CCs) using BURP analysis; novel/unknown spa types were further subjected to MLST analysis. The majority of isolates belonged to spa-CC 152 (62/170, 36.4%), followed by spa-CC 121 (19/170, 11.2%), and spa-CC 005 (18 /170, 10.6%). Of the nine MRSA isolates, 2 (22.2%) were spa-CC 239 with staphylococcal cassette chromosome (SCC)mec III. These findings highlight the diversity of S. aureus strains in Ethiopia, as well as the presence of potentially epidemic strains circulating in the country necessitating further characterization of S. aureus for antimicrobial resistance detection and infection prevention purposes.

Uncovering the Biotechnological Importance of Geotrichum candidum

Fungi make a fundamental contribution to several biotechnological processes, including brewing, winemaking, and the production of enzymes, organic acids, alcohols, antibiotics, and pharmaceuticals. The present review explores the biotechnological importance of the filamentous yeast-like fungus Geotrichum candidum, a ubiquitous species known for its use as a starter in the dairy industry. To uncover G. candidum's biotechnological role, we performed a search for related work through the scientific indexing internet services, Web of Science and Google Scholar. The following query was used: Geotrichum candidum, producing about 6500 scientific papers from 2017 to 2022. From these, approximately 150 that were associated with industrial applications of G. candidum were selected. Our analysis revealed that apart from its role as a starter in the dairy and brewing industries, this species has been administered as a probiotic nutritional supplement in fish, indicating improvements in developmental and immunological parameters. Strains of this species produce a plethora of biotechnologically important enzymes, including cellulases, β-glucanases, xylanases, lipases, proteases, and α-amylases. Moreover, strains that produce antimicrobial compounds and that are capable of bioremediation were identified. The findings of the present review demonstrate the importance of G. candidum for agrifood- and bio-industries and provide further insights into its potential future biotechnological roles.

Advanced Molecular-Genetic Methods and Prospects for Their Application for the Indication and Identification of <i>Yersinia pestis</i> Strains

The review provides an analysis of the literature data on the use of various modern molecular-genetic methods for the indication and identification of Yersinia pestis strains with different properties and degree of virulence, which is caused by the diverse natural conditions in which they circulate. The methods are also considered from the perspective of their promising application at three levels (territorial, regional and federal) of the system for laboratory diagnosis of infectious diseases at the premises of Rospotrebnadzor organizations to solve the problem of maintaining the sanitary and epidemiological well-being of the country’s population. The main groups of methods considered are as follows: based on the analysis of the lengths of restriction fragments (ribo- and IS-typing, pulse gel electrophoresis); based on the analysis of specific fragments (DFR typing, VNTR typing); based on sequencing (MLST, CRISPR analysis, SNP analysis); PCR methods (including IPCR, SPA); isothermal amplification methods (LAMP, HDA, RPA, SEA, PCA, SHERLOCK); DNA-microarray; methods using aptamer technology; bio- and nano-sensors; DNA origami; methods based on neural networks. We can conclude that the rapid development of molecular diagnostics and genetics is aimed at increasing efficiency, multi-factorial approaches and simplifying the application of techniques with no need for expensive equipment and highly qualified personnel for analysis. At all levels of the system for laboratory diagnosis of infectious diseases at the Rospotrebnadzor organizations, it is possible to use methods based on PCR, isothermal amplification, SHERLOCK, biosensors, and small-sized sequencing devices. At the territorial level, at plague control stations, the use of immuno-PCR and SPA for the indication of Y. pestis is viable. At the regional level, introduction of the technologies based on the use of aptamers and DNA chips looks promising. For the federal level, the use of DNA origami methods and new technologies of whole genome sequencing is a prospect within the framework of advanced identification, molecular typing and sequencing of the genomes of plague agent strains.

Listeria monocytogenes in foods-From culture identification to whole-genome characteristics

Listeria monocytogenes is an important foodborne pathogen, which is able to persist in the food production environments. The presence of these bacteria in different niches makes them a potential threat for public health. In the present review, the current information on the classical and alternative methods used for isolation and identification of L. monocytogenes in food have been described. Although these techniques are usually simple, standardized, inexpensive, and are routinely used in many food testing laboratories, several alternative molecular-based approaches for the bacteria detection in food and food production environments have been developed. They are characterized by the high sample throughput, a short time of analysis, and cost-effectiveness. However, these methods are important for the routine testing toward the presence and number of L. monocytogenes, but are not suitable for characteristics and typing of the bacterial isolates, which are crucial in the study of listeriosis infections. For these purposes, novel approaches, with a high discriminatory power to genetically distinguish the strains during epidemiological studies, have been developed, e.g., whole-genome sequence-based techniques such as NGS which provide an opportunity to perform comparison between strains of the same species. In the present review, we have shown a short description of the principles of microbiological, alternative, and modern methods of detection of L. monocytogenes in foods and characterization of the isolates for epidemiological purposes. According to our knowledge, similar comprehensive papers on such subject have not been recently published, and we hope that the current review may be interesting for research communities.

Prevalence, distribution, enterotoxin profiles, antimicrobial resistance, and genetic diversity of Bacillus cereus group isolates from lettuce farms in Korea

Lettuce wraps are popular in Korean cuisine for their high nutritional value and versatility as healthy additions to multiple dishes. Microbial contamination of lettuce is a major concern, as lettuce is consumed fresh without cooking. Among foodborne pathogens, the spore-forming, facultative anaerobic bacterium, Bacillus cereus is one of the frequently detected pathogen in lettuce in Korea. In this study, we investigated the prevalence and distribution of Bacillus cereus strains in lettuce production farms and further evaluated the enterotoxin gene profiles, antibiotic susceptibility, multidrug resistance pattern, and genetic differences among the B. cereus group isolates. Of the 140 samples isolated from 10 lettuce production farms, 30 samples (21.42%) were positive for B. cereus in which 19 (31.6%) and 10 (23.25%) were from soil and lettuce, respectively. The enterotoxin patterns A (hblCDA, nheABC, entFM, and cytK genes) and B (hblCDA, nheABC, and entFM genes) accounted for 50% and 20% of all the isolates, whereas the emetic gene cesB was not detected in any of the B. cereus group isolates. Antibiotic susceptibility testing of the B. cereus group isolates revealed that all the strains were predominantly resistant to β-lactam antibiotics except imipenem and generally susceptible to most of the non β-lactam antibiotics, including gentamycin, streptomycin, chloramphenicol, and tetracycline. ERIC-PCR and MLST analysis revealed high genetic diversity among the 30 B. cereus group isolates, which belonged to 26 different sequence types (STs) and seven new STs. Moreover, isolates with identical STs exhibited similar patterns of antibiotic resistance and enterotoxin profiles. Results of this study indicate a high prevalence of B. cereus group isolates in lettuce production farms in the Republic of Korea.

Genome Sequence of Escherichia coli Isolated from an Adult in Kibera, an Urban Informal Settlement in Nairobi, Kenya

An Escherichia coli strain (sequence type 636) was isolated from an adult residing in an urban informal settlement in Nairobi, Kenya, and was sequenced using the Illumina MiSeq platform. The draft genome was 5,075,726 bp, with a Col(BS512) plasmid plus aph(6)-Id, bla_TEM-1B, and dfrA7 genes, which encode kanamycin, ampicillin, and trimethoprim resistance proteins, respectively.

Multilocus Sequence Typing of Non-JP2 Serotype b Aggregatibacter actinomycetemcomitans Strains of Ghanaian and Swedish Origin

Objective and Methods

The Gram-negative bacterium, Aggregatibacter actinomycetemcomitans is associated with periodontitis affecting young individuals. The geographic dissemination of the highly leukotoxic JP2 genotype of serotype b of this species was previously studied by multilocus sequence typing (MLST). Here, we have used MLST to genetically characterize non-JP2 genotype strains of serotype b, isolated from individuals living in Ghana (n=41), and in Sweden (n=13), respectively.

Results

The MLST analysis revealed a total of nine sequence types (ST). Both Ghanaian and Swedish isolates were distributed in ST 1-3. ST 5 and 6 were only identified among the Ghanaian strains, whereas ST 4, 7, 8 and 9 were uniquely represented among the Swedish strains. Previously, we characterized these non-JP2 genotype strains of A. actinomycetemcomitans serotype b by arbitrarily-primed (AP)-PCR, which distributed them into three groups, AP-PCR type 1, 2, and 3, respectively. AP-PCR type 1 strains are generally highly leukotoxic, and are associated with progression of periodontal attachment loss. As AP-PCR type 1 includes both JP2 genotype strains and a proportion of non-JP2 genotype strains of serotype b, a straightforward diagnostic procedure has been sought. This has revealed a gene, cagE, which appears to be conserved only in this AP-PCR type. According to our results, MLST was not a highly discriminatory method to identify AP-PCR type 1, as strains of this AP-PCR type could be found within three different ST: ST 2, ST 3 and ST 8.

Conclusion

According to MLST, a geographic dissemination of non-JP2 genotype A. actinomycetemcomitans serotype b appears to exist. However, aiming to identify carriers of AP-PCR type 1, non-JP2 genotype serotype b, PCR with cagE-specific primers is likely the most efficient diagnostic procedure known today.

Molecular subtyping for source tracking of Escherichia coli using core genome multilocus sequence typing at a food manufacturing plant

When harmful bacteria are detected in the final product at a food manufacturing plant, it is necessary to identify and eliminate the source of contamination so that it does not occur again. In the current study, the source of contamination was tracked using core genome multilocus sequence typing (cgMLST) analysis in cases where Escherichia coli was detected in the final product at a food manufacturing plant. cgMLST analysis was performed on 40 strains of E. coli collected from the environment [floor (26 strains), drainage ditch (5 strains), container (4 strains), post-heating production line (1 strain)] and products [final product (3 strains) and intermediate product (1 strain)]. In total, 40 E. coli isolates were classified into 17 genogroups by cgMLST analysis. The 4 E. coli strains isolated from the intermediate and final products were classified into two genogroups (I and II). Certain isolates collected from the environment also belonged to those genogroups, it was possible to estimate the transmission of E. coli in the manufacturing plant. Thus, the dynamics of E. coli in the food manufacturing location were clarified by using cgMLST analysis. In conclusion, our results indicate that cgMLST analysis can be effectively used for hygiene management at food manufacturing locations.

Molecular Investigation of Recurrent Streptococcus iniae Epizootics Affecting Coral Reef Fish on an Oceanic Island Suggests at Least Two Distinct Emergence Events

Streptococcus iniae is an emerging zoonotic pathogen of increasing concern for aquaculture and has caused several epizootics in reef fishes from the Caribbean, the Red Sea and the Indian Ocean. To study the population structure, introduction pathways and evolution of S. iniae over recurring epizootics on Reunion Island, we developed and validated a Multi Locus Sequence Typing (MLST) panel using genomic data obtained from 89 isolates sampled during epizootics occurring over the past 40years in Australia, Asia, the United States, Israel and Reunion Island. We selected eight housekeeping loci, which resulted in the greatest variation across the main S. iniae phylogenetic clades highlighted by the whole genomic dataset. We then applied the developed MLST to investigate the origin of S. iniae responsible for four epizootics on Reunion Island, first in inland aquaculture and then on the reefs from 1996 to 2014. Results suggest at least two independent S. iniae emergence events occurred on the island. Molecular data support that the first epizootic resulted from an introduction, with inland freshwater aquaculture facilities acting as a stepping-stone. Such an event may have been facilitated by the ecological flexibility of S. iniae, able to survive in both fresh and marine waters and the ability of the pathogen to infect multiple host species. By contrast, the second epizootic was associated with a distinct ST of cosmopolitan distribution that may have emerged as a result of environment disturbance. This novel tool will be effective at investigating recurrent epizootics occurring within a given environment or country that is despite the fact that S. iniae appears to have low genetic diversity within its lineage.

Genomic Surveillance and Improved Molecular Typing of Bordetella pertussis Using wgMLST

Multilocus sequence typing (MLST) provides allele-based characterization of bacterial pathogens in a standardized framework. However, classical MLST schemes for Bordetella pertussis, the causative agent of whooping cough, seldom reveal diversity among the small number of gene targets and thereby fail to delineate population structure. To improve the discriminatory power of allele-based molecular typing of B. pertussis, we have developed a whole-genome MLST (wgMLST) scheme from 225 reference-quality genome assemblies. Iterative refinement and allele curation resulted in a scheme of 3,506 coding sequences and covering 81.4% of the B. pertussis genome. This wgMLST scheme was further evaluated with data from a convenience sample of 2,389 B. pertussis isolates sequenced on Illumina instruments, including isolates from known outbreaks and epidemics previously characterized by existing molecular assays, as well as replicates collected from individual patients. wgMLST demonstrated concordance with whole-genome single nucleotide polymorphism (SNP) profiles, accurately resolved outbreak and sporadic cases in a retrospective comparison, and clustered replicate isolates collected from individual patients during diagnostic confirmation. Additionally, a reanalysis of isolates from two statewide epidemics using wgMLST reconstructed the population structures of circulating strains with increased resolution, revealing new clusters of related cases. Comparison with an existing core genome (cgMLST) scheme highlights the stable gene content of this bacterium and forms the initial foundation for necessary standardization. These results demonstrate the utility of wgMLST for improving B. pertussis characterization and genomic surveillance during the current pertussis disease resurgence.

Antibiotic Resistance and Genotypes of Nosocomial Strains of Acinetobacter baumannii in Kazakhstan

The aim of this study was to determine the prevalence of A. baumannii antibiotic-resistant strains in Kazakhstan and to characterize genotypes related to epidemic “high-risk” clones. Two hundred and twenty four A. baumannii isolates from four cities of Kazakhstan in 2011–2019 were studied. Antibiotic susceptibility testing was performed by using broth microdilutions method according to EUCAST (v 11.0) recommendations. The presence of bla_OXA-23-like, bla_{OXA-24/40-like,}bla_OXA-58-like,bla_VIM,bla_IMP, and bla_NDM genes was determined by PCR. Genotyping was performed using high-throughput real-time PCR detection of 21 SNPs at 10 chromosomal loci used in existing MLST schemes. Resistance rates to imipenem, meropenem, amikacin, gentamicin, and ciprofloxacin were 81.3%, 78.6%, 79.9%, 65.2%, and 89.3%, respectively. No colistin resistant isolates were detected. The values of the MIC 50% and the MIC 90% of tigecycline were 0.125 mg/L, only four isolates (1.8%) had the ECOFF value >0.5 mg/L. The presence of acquired carbapenemase genes was found in 82.2% strains, including bla_OXA-23-like (78.6%) or bla_OXA-58-like (3.6%) genes. The spreading of carbapenem resistant A. baumannii strains in Kazakhstan was associated with epidemic “high-risk” clonal groups, predominantly, CG208(92)^OXF/CG2^PAS (80.8%) and less often CG231(109)^OXF/CG1^PAS (1.8%).

Phenotypic and Genetic Characterization of the Cheese Ripening Yeast Geotrichum candidum

The yeast Geotrichum candidum (teleomorph Galactomyces candidus) is inoculated onto mold- and smear-ripened cheeses and plays several roles during cheese ripening. Its ability to metabolize proteins, lipids, and organic acids enables its growth on the cheese surface and promotes the development of organoleptic properties. Recent multilocus sequence typing (MLST) and phylogenetic analyses of G. candidum isolates revealed substantial genetic diversity, which may explain its strain-dependant technological capabilities. Here, we aimed to shed light on the phenotypic and genetic diversity among eight G. candidum and three Galactomyces spp. strains of environmental and dairy origin. Phenotypic tests such as carbon assimilation profiles, the ability to grow at 35°C and morphological traits on agar plates allowed us to discriminate G. candidum from Galactomyces spp. The genomes of these isolates were sequenced and assembled; whole genome comparison clustered the G. candidum strains into three subgroups and provided a reliable reference for MLST scheme optimization. Using the whole genome sequence as a reference, we optimized an MLST scheme using six loci that were proposed in two previous MLST schemes. This new MLST scheme allowed us to identify 15 sequence types (STs) out of 41 strains and revealed three major complexes named GeoA, GeoB, and GeoC. The population structure of these 41 strains was evaluated with STRUCTURE and a NeighborNet analysis of the combined six loci, which revealed recombination events between and within the complexes. These results hint that the allele variation conferring the different STs arose from recombination events. Recombination occurred for the six housekeeping genes studied, but most likely occurred throughout the genome. These recombination events may have induced an adaptive divergence between the wild strains and the cheesemaking strains, as observed for other cheese ripening fungi. Further comparative genomic studies are needed to confirm this phenomenon in G. candidum. In conclusion, the draft assembly of 11 G. candidum/Galactomyces spp. genomes allowed us to optimize a genotyping MLST scheme and, combined with the assessment of their ability to grow under different conditions, provides a reliable tool to cluster and eventually improves the selection of G. candidum strains.

Hospital outbreak caused by linezolid resistant Enterococcus faecium in Upper Austria

Background

Enterococcus faecium is part of the human gastrointestinal flora but may act as opportunistic pathogen. Environmental persistence, high colonization capability and diverse intrinsic and acquired resistance mechanisms make it especially successful in nosocomial high-risk settings. In March 2014, an outbreak of Linezolid resistant Enterococcus faecium (LREfm) was observed at the hematooncology department of a tertiary care center in Upper Austria.

Methods

We report on the outbreak investigation together with the whole genome sequencing (WGS)-based typing results including also non-outbreak LREfm and susceptible isolates.

Results

The 54 investigated isolates could be divided in six clusters based on cgMLST. Cluster one comprised LREfm isolates of genotype ST117 and CT24, which was identified as the causative clone of the outbreak. In addition, the detection of four other clusters comprising isolates originating from hematooncology patients but also at other hospitals, pointed to LREfm transmission between local healthcare facilities. LREfm patients (n = 36) were typically at risk for acquisition of nosocomial pathogens because of immunosuppression, frequent hospitalization and antibiotic therapies. Seven of these 36 patients developed LREfm infection but were successfully treated. After termination of the initial outbreak, sporadic cases occurred despite a bundle of applied outbreak control interventions.

Conclusions

WGS proved to be an effective tool to differentiate several LREfm clusters in an outbreak. Active screening for LREfm is important in a high-risk setting such as hematooncology, where multiple introductions are possible and occur despite intensified infection control measures.

Electronic supplementary material

The online version of this article (10.1186/s13756-019-0598-z) contains supplementary material, which is available to authorized users.

Fourier-Transform InfraRed Spectroscopy Can Quickly Type Gram-Negative Bacilli Responsible for Hospital Outbreaks

The typing of epidemic bacterial pathogens in hospitals relies on DNA-based, expensive, and time-consuming techniques, that are often limited to retrospective studies. However, the quick identification of epidemic pathogens in the routine of the microbiology laboratories would expedite infection control procedures that limit the contamination of new patients. IR Biotyper (Bruker Daltonics GmbH) is a new typing machine based on Fourier-transform infrared (FTIR) spectroscopy which generates spectra, aiming at typing the micro-organisms within 3 h. This technique discriminates the isolates by exploring the differences of the surface cell polysaccharides. In this work, we evaluated the ability of the FTIR spectroscopy to recognize Gram-negative bacilli clones responsible for hospital outbreaks. Isolates of Pseudomonas aeruginosa (n = 100), Klebsiella pneumoniae (n = 16), Enterobacter cloacae (n = 23), and Acinetobacter baumannii (n = 20) were typed by the reference methods Multi-Locus Sequence Typing (defining sequence types – STs) along with or without pulsed field gel electrophoresis (PFGE) (defining pulsotypes), and by FTIR spectroscopy. The congruence of FTIR spectroscopy clustering was compared to those of MLST and PFGE by Adjusted Rand index and Adjusted Wallace coefficient. We found that FTIR spectroscopy accurately clustered P. aeruginosa, K. pneumoniae, and E. cloacae isolates belonging to the same ST. The performance of the FTIR spectroscopy was slightly lower for A. baumannii. Furthermore, FTIR spectroscopy also correctly clustered P. aeruginosa isolates having a similar pulsotype. Overall, the IR Biotyper can quickly (in less than 3 h) detect the spread of clones of P. aeruginosa, K. pneumoniae, E. cloacae, and A. baumannii. The use of this technique by clinical microbiology laboratories may help to tackle the spread of epidemic clones by the quick implementation of infection control measures.

Parallel sequencing of porA reveals a complex pattern of Campylobacter genotypes that differs between broiler and broiler breeder chickens

Chicken meat represents an important source of Campylobacter infections of humans world-wide. A better understanding of Campylobacter epidemiology in commercial chicken flocks will facilitate the development of more effective intervention strategies. We developed a gene-specific parallel sequencing approach that efficiently indicated genetic diversity in farm-derived samples and revealed Campylobacter genotypes that would not be detected using microbiological culture. Parallel sequencing of the porA nucleotide fragment identified a different pattern of diversity in broiler flocks compared with broiler-breeder flocks at both individual bird and flock levels. Amongst the flocks tested, broiler flocks and individual birds were dominated by one or two porA fragment types whereas co-dominance with up to six porA fragment types was evident in breeder birds. A high proportion (83.6-93.3%) of porA variants were shared between broiler and breeder flocks. The porA-based diversity profiling could be a useful addition to the repertoire of tools employed to attribute potential sources of contamination for broiler flocks, including the environment, wild animals or other chickens. This approach can be extended to include other loci within Campylobacter and developed for molecular epidemiology studies of other bacterial species.

Epidemiological Typing of Serratia marcescens Isolates by Whole-Genome Multilocus Sequence Typing

Serratia marcescens is an opportunistic bacterial pathogen. It is notorious for its increasing antimicrobial resistance and its potential to cause outbreaks of colonization and infections, predominantly in neonatal intensive care units (NICUs). There, its spread requires rapid infection control response. To understand its spread, detailed molecular typing is key. We present a whole-genome multilocus sequence typing (wgMLST) method for S. marcescens Using a set of 299 publicly available whole-genome sequences (WGS), we developed an initial wgMLST system consisting of 9,377 gene loci. This included 1,455 loci occurring in all reference genomes and 7,922 accessory loci. This closed system was validated using three geographically diverse collections of S. marcescens consisting of 111 clinical isolates implicated in nosocomial dissemination events in three hospitals. The validation procedure showed a full match between epidemiological data and the wgMLST analyses. We set the cutoff value for epidemiological (non)relatedness at 20 different alleles, though for the majority of outbreak-clustered isolates, this difference was limited to 4 alleles. This shows that the wgMLST system for S. marcescens provides prospects for successful future monitoring for the epidemiological containment of this opportunistic pathogen.

Genomic Surveillance of 4CMenB Vaccine Antigenic Variants among Disease-Causing Neisseria meningitidis Isolates, United Kingdom, 2010-2016

In September 2015, 4CMenB meningococcal vaccine was introduced into the United Kingdom infant immunization program without phase 3 trial information. Understanding the effect of this program requires enhanced surveillance of invasive meningococcal disease (IMD) Neisseria meningitidis isolates and comparison with prevaccination isolates. Bexsero Antigen Sequence Types (BASTs) were used to analyze whole-genome sequences of 3,073 prevaccine IMD N. meningitidis isolates obtained during 2010−2016. Isolates exhibited 803 BASTs among 31 clonal complexes. Frequencies of antigen peptide variants were factor H binding protein 1, 13.4%; Neisserial heparin-binding antigen 2, 13.8%; Neisseria adhesin A 8, 0.8%; and Porin A-VR2:P1.4,10.9%. In 2015−16, serogroup B isolates showed the highest proportion (35.7%) of exact matches to >1 Bexsero components. Serogroup W isolates showed the highest proportion (93.9%) of putatively cross-reactive variants of Bexsero antigens. Results highlighted the likely role of cross-reactive antigens. BAST surveillance of meningococcal whole-genome sequence data is rapid, scalable, and portable and enables international comparisons of isolates.

Whole genome sequencing for investigations of meningococcal outbreaks in the United States: a retrospective analysis

Although rare in the U.S., outbreaks due to Neisseria meningitidis do occur. Rapid, early outbreak detection is important for timely public health response. In this study, we characterized U.S. meningococcal isolates (N = 201) from 15 epidemiologically defined outbreaks (2009-2015) along with temporally and geographically matched sporadic isolates using multilocus sequence typing, pulsed-field gel electrophoresis (PFGE), and six whole genome sequencing (WGS) based methods. Recombination-corrected maximum likelihood (ML) and Bayesian phylogenies were reconstructed to identify genetically related outbreak isolates. All WGS analysis methods showed high degree of agreement and distinguished isolates with similar or indistinguishable PFGE patterns, or the same strain genotype. Ten outbreaks were caused by a single strain; 5 were due to multiple strains. Five sporadic isolates were phylogenetically related to 2 outbreaks. Analysis of 9 outbreaks using timed phylogenies identified the possible origin and estimated the approximate time that the most recent common ancestor emerged for outbreaks analyzed. U.S. meningococcal outbreaks were caused by single- or multiple-strain introduction, with organizational outbreaks mainly caused by a clonal strain and community outbreaks by divergent strains. WGS can infer linkage of meningococcal cases when epidemiological links are uncertain. Accurate identification of outbreak-associated cases requires both WGS typing and epidemiological data.

Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications

The PubMLST.org website hosts a collection of open-access, curated databases that integrate population sequence data with provenance and phenotype information for over 100 different microbial species and genera.  Although the PubMLST website was conceived as part of the development of the first multi-locus sequence typing (MLST) scheme in 1998 the software it uses, the Bacterial Isolate Genome Sequence database (BIGSdb, published in 2010), enables PubMLST to include all levels of sequence data, from single gene sequences up to and including complete, finished genomes.  Here we describe developments in the BIGSdb software made from publication to June 2018 and show how the platform realises microbial population genomics for a wide range of applications.  The system is based on the gene-by-gene analysis of microbial genomes, with each deposited sequence annotated and curated to identify the genes present and systematically catalogue their variation.  Originally intended as a means of characterising isolates with typing schemes, the synthesis of sequences and records of genetic variation with provenance and phenotype data permits highly scalable (whole genome sequence data for tens of thousands of isolates) means of addressing a wide range of functional questions, including: the prediction of antimicrobial resistance; likely cross-reactivity with vaccine antigens; and the functional activities of different variants that lead to key phenotypes.  There are no limitations to the number of sequences, genetic loci, allelic variants or schemes (combinations of loci) that can be included, enabling each database to represent an expanding catalogue of the genetic variation of the population in question.  In addition to providing web-accessible analyses and links to third-party analysis and visualisation tools, the BIGSdb software includes a RESTful application programming interface (API) that enables access to all the underlying data for third-party applications and data analysis pipelines.

Microbial sequence typing in the genomic era

Next-generation sequencing (NGS), also known as high-throughput sequencing, is changing the field of microbial genomics research. NGS allows for a more comprehensive analysis of the diversity, structure and composition of microbial genes and genomes compared to the traditional automated Sanger capillary sequencing at a lower cost. NGS strategies have expanded the versatility of standard and widely used typing approaches based on nucleotide variation in several hundred DNA sequences and a few gene fragments (MLST, MLVA, rMLST and cgMLST). NGS can now accommodate variation in thousands or millions of sequences from selected amplicons to full genomes (WGS, NGMLST and HiMLST). To extract signals from high-dimensional NGS data and make valid statistical inferences, novel analytic and statistical techniques are needed. In this review, we describe standard and new approaches for microbial sequence typing at gene and genome levels and guidelines for subsequent analysis, including methods and computational frameworks. We also present several applications of these approaches to some disciplines, namely genotyping, phylogenetics and molecular epidemiology.

GrapeTree: visualization of core genomic relationships among 100,000 bacterial pathogens

Current methods struggle to reconstruct and visualize the genomic relationships of large numbers of bacterial genomes. GrapeTree facilitates the analyses of large numbers of allelic profiles by a static “GrapeTree Layout” algorithm that supports interactive visualizations of large trees within a web browser window. GrapeTree also implements a novel minimum spanning tree algorithm (MSTree V2) to reconstruct genetic relationships despite high levels of missing data. GrapeTree is a stand-alone package for investigating phylogenetic trees plus associated metadata and is also integrated into EnteroBase to facilitate cutting edge navigation of genomic relationships among bacterial pathogens.

Systematic characterization of Bacillus Genetic Stock Center Bacillus thuringiensis strains using Multi-Locus Sequence Typing

The goal of this work was to perform a systematic characterization of Bacillus thuringiensis (Bt) strains from the Bacillus Genetic Stock Center (BGSC) collection using Multi-Locus Sequence Typing (MLST). Different genetic markers of 158 Bacillus thuringiensis (Bt) strains from 73 different serovars stored in the BGSC, that represented 92% of the different Bt serovars of the BGSC were analyzed, the 8% that were not analyzed were not available. In addition, we analyzed 72 Bt strains from 18 serovars available at the pubMLST bcereus database, and Bt strains G03, HBF18 and Bt185, with no H serovars provided by our laboratory. We performed a systematic MLST analysis using seven housekeeping genes (glpF, gmK, ilvD, pta, pur, pycA and tpi) and analyzed correlation of the results of this analysis with strain serovars. The 233 Bt strains analyzed were assigned to 119 STs from which 19 STs were new. Genetic relationships were established by phylogenetic analysis and showed that STs could be grouped in two major Clusters containing 21 sub-groups. We found that a significant number of STs (101 in total) correlated with specific serovars, such as ST13 that corresponded to nine Bt isolates from B. thuringiensis serovar kenyae. However, other serovars showed high genetic variability and correlated with multiple STs; for example, B. thuringiensis serovar morrisoni correlated with 11 different STs. In addition, we found that 16 different STs correlated with multiple serovars (2-4 different serovars); for example, ST12 correlated with B. thuringiensis serovar alesti, dakota, palmanyolensis and sotto/dendrolimus. These data indicated that only partial correspondence between MLST and serotyping can be established.

Draft Genome Sequences of Four Clinical Legionella pneumophila Isolates from Ontario, Canada

Legionella pneumophila outbreak investigations require the development of reliable typing methods to better understand the genetic relationships of the isolates involved. Here, we report the draft genome sequences of four clinical Legionella pneumophila isolates obtained between 2000 and 2012 in Ontario, Canada.

The Applied Development of a Tiered Multilocus Sequence Typing (MLST) Scheme for Dichelobacter nodosus

Dichelobacter nodosus (D. nodosus) is the causative pathogen of ovine footrot, a disease that has a significant welfare and financial impact on the global sheep industry. Previous studies into the phylogenetics of D. nodosus have focused on Australia and Scandinavia, meaning the current diversity in the United Kingdom (U.K.) population and its relationship globally, is poorly understood. Numerous epidemiological methods are available for bacterial typing; however, few account for whole genome diversity or provide the opportunity for future application of new computational techniques. Multilocus sequence typing (MLST) measures nucleotide variations within several loci with slow accumulation of variation to enable the designation of allele numbers to determine a sequence type. The usage of whole genome sequence data enables the application of MLST, but also core and whole genome MLST for higher levels of strain discrimination with a negligible increase in experimental cost. An MLST database was developed alongside a seven loci scheme using publically available whole genome data from the sequence read archive. Sequence type designation and strain discrimination was compared to previously published data to ensure reproducibility. Multiple D. nodosus isolates from U.K. farms were directly compared to populations from other countries. The U.K. isolates define new clades within the global population of D. nodosus and predominantly consist of serogroups A, B and H, however serogroups C, D, E, and I were also found. The scheme is publically available at https://pubmlst.org/dnodosus/.

Large-scale genomic analyses reveal the population structure and evolutionary trends of Streptococcus agalactiae strains in Brazilian fish farms

Streptococcus agalactiae is a major pathogen and a hindrance on tilapia farming worldwide. The aims of this work were to analyze the genomic evolution of Brazilian strains of S. agalactiae and to establish spatial and temporal relations between strains isolated from different outbreaks of streptococcosis. A total of 39 strains were obtained from outbreaks and their whole genomes were sequenced and annotated for comparative analysis of multilocus sequence typing, genomic similarity and whole genome multilocus sequence typing (wgMLST). The Brazilian strains presented two sequence types, including a newly described ST, and a non-typeable lineage. The use of wgMLST could differentiate each strain in a single clone and was used to establish temporal and geographical correlations among strains. Bayesian phylogenomic analysis suggests that the studied Brazilian population was co-introduced in the country with their host, approximately 60 years ago. Brazilian strains of S. agalactiae were shown to be heterogeneous in their genome sequences and were distributed in different regions of the country according to their genotype, which allowed the use of wgMLST analysis to track each outbreak event individually.

DNA-Sequence Based Typing of the Cronobacter Genus Using MLST, CRISPR-cas Array and Capsular Profiling

The Cronobacter genus is composed of seven species, within which a number of pathovars have been described. The most notable infections by Cronobacter spp. are of infants through the consumption of contaminated infant formula. The description of the genus has greatly improved in recent years through DNA sequencing techniques, and this has led to a robust means of identification. However some species are highly clonal and this limits the ability to discriminate between unrelated strains by some methods of genotyping. This article updates the application of three genotyping methods across the Cronobacter genus. The three genotyping methods were multilocus sequence typing (MLST), capsular profiling of the K-antigen and colanic acid (CA) biosynthesis regions, and CRISPR-cas array profiling. A total of 1654 MLST profiled and 286 whole genome sequenced strains, available by open access at the PubMLST Cronobacter database, were used this analysis. The predominance of C. sakazakii and C. malonaticus in clinical infections was confirmed. The majority of clinical strains being in the C. sakazakii clonal complexes (CC) 1 and 4, sequence types (ST) 8 and 12 and C. malonaticus ST7. The capsular profile K2:CA2, previously proposed as being strongly associated with C. sakazakii and C. malonaticus isolates from severe neonatal infections, was also found in C. turicensis, C. dublinensis and C. universalis. The majority of CRISPR-cas types across the genus was the I-E (Ecoli) type. Some strains of C. dublinensis and C. muytjensii encoded the I-F (Ypseudo) type, and others lacked the cas gene loci. The significance of the expanding profiling will be of benefit to researchers as well as governmental and industrial risk assessors.

Will the emergence of core genome MLST end the role of in silico MLST?

The technological advancement of molecular epidemiological analysis using next-generation sequencing (NGS) for foodborne pathogens has a groundbreaking impact over the past three years. In particular, the emergence of cg (core genome) multilocus sequence typing(MLST) has a significant impact. This is because this technology made it possible for many researchers to carry out molecular epidemiological analysis on foodborne pathogens in a common language, using common definitions. The resolution of core genome MLST (cgMLST) far surpasses that of MLST, which only uses seven (usually, in some cases five) housekeeping genes. Therefore, cgMLST would in no doubt terminate the role of conventional MLST as the molecular epidemiological tool. However, the role of MLST would probably not end all together. Rather, the sequence type (ST) of the conventional MLST is expected to be used as in silico MLST by a wider range of researchers than ever in the next 10 years. This is because, with the arrival of the NGS era, we have come to be able to obtain ST of conventional MLST by simply entering the NGS text file into one's own PC. In other words, acquisition of ST data is no longer limited to researchers aiming to conduct MLST for the first place. The impact of such a change is large. In silico MLST will continue to be used as a tool for understanding the broad characteristics of bacterial strains. This review aimed to summarize the main information on STs that have been accumulated for representative foodborne pathogens, in particular for potential NGS users in this new era who have been not familiar with MLST until now.

Defining a Core Genome Multilocus Sequence Typing Scheme for the Global Epidemiology of Vibrio parahaemolyticus

Vibrio parahaemolyticus is an important human foodborne pathogen whose transmission is associated with the consumption of contaminated seafood, with a growing number of infections reported over recent years worldwide. A multilocus sequence typing (MLST) database for V. parahaemolyticus was created in 2008, and a large number of clones have been identified, causing severe outbreaks worldwide (sequence type 3 [ST3]), recurrent outbreaks in certain regions (e.g., ST36), or spreading to other regions where they are nonendemic (e.g., ST88 or ST189). The current MLST scheme uses sequences of 7 genes to generate an ST, which results in a powerful tool for inferring the population structure of this pathogen, although with limited resolution, especially compared to pulsed-field gel electrophoresis (PFGE). The application of whole-genome sequencing (WGS) has become routine for trace back investigations, with core genome MLST (cgMLST) analysis as one of the most straightforward ways to explore complex genomic data in an epidemiological context. Therefore, there is a need to generate a new, portable, standardized, and more advanced system that provides higher resolution and discriminatory power among V. parahaemolyticus strains using WGS data. We sequenced 92 V. parahaemolyticus genomes and used the genome of strain RIMD 2210633 as a reference (with a total of 4,832 genes) to determine which genes were suitable for establishing a V. parahaemolyticus cgMLST scheme. This analysis resulted in the identification of 2,254 suitable core genes for use in the cgMLST scheme. To evaluate the performance of this scheme, we performed a cgMLST analysis of 92 newly sequenced genomes, plus an additional 142 strains with genomes available at NCBI. cgMLST analysis was able to distinguish related and unrelated strains, including those with the same ST, clearly showing its enhanced resolution over conventional MLST analysis. It also distinguished outbreak-related from non-outbreak-related strains within the same ST. The sequences obtained from this work were deposited and are available in the public database (http://pubmlst.org/vparahaemolyticus). The application of this cgMLST scheme to the characterization of V. parahaemolyticus strains provided by different laboratories from around the world will reveal the global picture of the epidemiology, spread, and evolution of this pathogen and will become a powerful tool for outbreak investigations, allowing for the unambiguous comparison of strains with global coverage.

Genetic and genomic diversity of NheABC locus from Bacillus strains

Non-hemolytic enterotoxin (NHE), a tri-partite, proteinaceous toxin encoded by contiguous nheA, nheB and nheC genes of Bacillus cereus sensu lato (B. cereus s.l.), is considered to be associated with the foodborne diarrheic syndrome. However, B. cereus s.l. includes a number of closely related strains, and the occurrence of NHE among them, and other members of Bacillus is unclear. Consequently, we aimed to determine the distribution and evolution of NHE within Bacillus by confirming the presence of the nheA, B and C sequences and variation within them using published data, and to analyze the genomic and genetic diversity. The phylogenetic tree of NHE proteins (NheA, NheB and NheC) from 81 different B. cereus s.l. strains was constructed. And on the genetic determinants of the NHE toxin did not bring any obvious link between the nheABC genes sequence of a strain and its virulence in the diarrhoeal pathogenesis. Analysis of the genomic diversity of the nheA, B and C loci revealed that their upstream regions were more conserved than the downstream sequences. Multilocus sequence typing schemes (MLST) based on seven concatenated housekeeping genes and nheA, B and C genes of the 75 strains were developed. The neighbor joining phylogenetic tree based on seven housekeeping genes together with nheA, B and C genes was similiar to published Bacillus phylogenetic trees. And on the genetic determinants of the NHE toxin did not bring any obvious link between the nheABC genes sequence of a strain and its virulence in the diarrhoeal pathogenesis.The results indicate that nheA, B and C genes do not affect the diversity of housekeeping genes, and this specific NHE protein does not participate in the diarrheic syndrome.

Molecular Assay Validation Using Genomic Sequence Databases

Whole-genome sequence databases offer new in silico approaches for designing and validating PCR assays in the clinical microbiology laboratory. An article in this issue of the Journal of Clinical Microbiology (M. J. Jansen van Rensburg, C. Swift, A. J. Cody, C. Jenkins, and M. C. J. Maiden, J Clin Microbiol, 54:2882-2890, 2016, http://dx.doi.org/10.1128/JCM.01522-16) demonstrates the use of publicly available genomic sequence data to evaluate a PCR assay for distinguishing Campylobacter species.

Next-Generation Epidemiology: Using Real-Time Core Genome Multilocus Sequence Typing To Support Infection Control Policy

Multidrug-resistant bacteria are responsible for substantial morbidity and mortality worldwide. Tracking the nosocomial spread of resistant bacteria is critical to infection control. Mellmann et al. (J. Clin. Microbiol. 54:2874-2881, 2016, http://dx.doi.org/10.1128/JCM.00790-16) have described prospective whole-genome sequencing with core genome multilocus sequencing typing (cgMLST) analysis for real-time surveillance and have addressed the practical aspects of implementing this type of operation in the hospital setting.

Population Genetic Structure of Listeria monocytogenes Strains as Determined by Pulsed-Field Gel Electrophoresis and Multilocus Sequence Typing

Listeria monocytogenes is a ubiquitous bacterium that may cause the foodborne illness listeriosis. Only a small amount of data about the population genetic structure of strains isolated from food is available. This study aimed to provide an accurate view of the L. monocytogenes food strain population in France. From 1999 to 2014, 1,894 L. monocytogenes strains were isolated from food at the French National Reference Laboratory for L. monocytogenes and classified according to the five risk food matrices defined by the European Food Safety Authority (EFSA). A total of 396 strains were selected on the basis of different pulsed-field gel electrophoresis (PFGE) clusters, serotypes, and strain origins and typed by multilocus sequence typing (MLST), and the MLST results were supplemented with MLST data available from Institut Pasteur, representing human and additional food strains from France. The distribution of sequence types (STs) was compared between food and clinical strains on a panel of 675 strains. High congruence between PFGE and MLST was found. Out of 73 PFGE clusters, the two most prevalent corresponded to ST9 and ST121. Using original statistical analysis, we demonstrated that (i) there was not a clear association between ST9 and ST121 and the food matrices, (ii) serotype IIc, ST8, and ST4 were associated with meat products, and (iii) ST13 was associated with dairy products. Of the two major STs, ST121 was the ST that included the fewest clinical strains, which might indicate lower virulence. This observation may be directly relevant for refining risk analysis models for the better management of food safety.

IMPORTANCE This study showed a very useful backward compatibility between PFGE and MLST for surveillance. The results enabled better understanding of the population structure of L. monocytogenes strains isolated from food and management of the health risks associated with L. monocytogenes food strains. Moreover, this work provided an accurate view of L. monocytogenes strain populations associated with specific food matrices. We clearly showed that some STs were associated with food matrices, such as meat, meat products, and dairy products. We opened the way to source attribution modeling in order to quantify the relative importance of the main food matrices.

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.

IMPORTANCE

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.

Does typing of Chlamydia trachomatis using housekeeping multilocus sequence typing reveal different sexual networks among heterosexuals and men who have sex with men?

Background

Chlamydia trachomatis infections remain the most common bacterial sexually transmitted infection worldwide. To gain more insight into the epidemiology and transmission of C. trachomatis, several schemes of multilocus sequence typing (MLST) have been developed. We investigated the clustering of C. trachomatis strains derived from men who have sex with men (MSM) and heterosexuals using the MLST scheme based on 7 housekeeping genes (MLST-7) adapted for clinical specimens and a high-resolution MLST scheme based on 6 polymorphic genes, including ompA (hr-MLST-6).

Methods

Specimens from 100 C. trachomatis infected men who have sex with men (MSM) and 100 heterosexual women were randomly selected from previous studies and sequenced. We adapted the MLST-7 scheme to a nested assay to be suitable for direct typing of clinical specimens. All selected specimens were typed using both the adapted MLST-7 scheme and the hr-MLST-6 scheme. Clustering of C. trachomatis strains derived from MSM and heterosexuals was assessed using minimum spanning tree analysis.

Results

Sufficient chlamydial DNA was present in 188 of the 200 (94 %) selected samples. Using the adapted MLST-7 scheme, full MLST profiles were obtained for 187 of 188 tested specimens resulting in a high success rate of 99.5 %. Of these 187 specimens, 91 (48.7 %) were from MSM and 96 (51.3 %) from heterosexuals. We detected 21 sequence types (STs) using the adapted MLST-7 and 79 STs using the hr-MLST-6 scheme. Minimum spanning tree analyses was used to examine the clustering of MLST-7 data, which showed no reflection of separate transmission in MSM and heterosexual hosts. Moreover, typing using the hr-MLST-6 scheme identified genetically related clusters within each of clusters that were identified by using the MLST-7 scheme.

Conclusion

No distinct transmission of C. trachomatis could be observed in MSM and heterosexuals using the adapted MLST-7 scheme in contrast to using the hr-MLST-6. In addition, we compared clustering of both MLST schemes and demonstrated that typing using the hr-MLST-6 scheme is able to identify genetically related clusters of C. trachomatis strains within each of the clusters that were identified by using the MLST-7 scheme.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-016-1486-2) contains supplementary material, which is available to authorized users.

Evaluation of two multi-locus sequence typing schemes for commensal Escherichia coli from dairy cattle in Washington State

Multi-locus sequence typing (MLST) is a useful system for phylogenetic and epidemiological studies of multidrug-resistant Escherichiacoli. Most studies utilize a seven-locus MLST, but an alternate two-locus typing method (fumC and fimH; CH typing) has been proposed that may offer a similar degree of discrimination at lower cost. Herein, we compare CH typing to the standard seven-locus method for typing commensal E. coli isolates from dairy cattle. In addition, we evaluated alternative combinations of eight loci to identify combinations that maximize discrimination and congruence with standard seven-locus MLST among commensal E. coli while minimizing the cost. We also compared both methods when used for typing uropathogenic E. coli (UPEC). CH typing was less discriminatory for commensal E. coli than the standard seven-locus method (Simpson's Index of Diversity=0.933 [0.902-0.964] and 0.97 [0.96-0.979], respectively). Combining fimH with housekeeping gene loci improved discriminatory power for commensal E. coli from cattle but resulted in poor congruence with MLST. We found that a four-locus typing method including the housekeeping genes adk, purA, gyrB and recA could be used to minimize cost without sacrificing discriminatory power or congruence with Achtman seven-locus MLST when typing commensal E. coli.

Genomic Analysis of Serogroup Y Neisseria meningitidis Isolates Reveals Extensive Similarities Between Carriage-Associated and Disease-Associated Organisms

BACKGROUND

Neisseria meningitidis is a frequent colonizer of the human nasopharynx, with asymptomatic carriage providing the reservoir for invasive, disease-causing strains. Serogroup Y (MenY) strains are a major cause of meningococcal disease. High-resolution genetic analyses of carriage and disease isolates can establish epidemiological relationships and identify potential virulence factors.

METHODS

Whole-genome sequence data were obtained for 99 MenY carriage isolates recovered in the United Kingdom during 1997-2010. Sequences were compared to those of 73 MenY invasive isolates recovered during 2010-2011, using a gene-by-gene approach.

RESULTS

Comparisons across 1605 core genes resolved 91% of isolates into one of 8 clusters containing closely related disease and carriage isolates. Six clusters contained carried meningococci isolated during 1997-2001, suggesting temporal stability. One cluster of isolates, predominately sharing the designation Y: P1.5-1,10-1: F4-1: ST-1655 (cc23), was resolved into one subcluster with 86% carriage isolates and a second with 90% invasive isolates. These subclusters were defined by specific allelic differences in 5 core genes encoding glycerate kinase (glxK), valine-pyruvate transaminase (avtA), superoxide dismutase (sodB), and 2 hypothetical proteins.

CONCLUSIONS

High-resolution genetic analyses detected long-term temporal stability and temporally overlapping carriage and disease populations for MenY clones but also evidence of a disease-associated clone.

Whole genome multilocus sequence typing as an epidemiologic tool for Yersinia pestis

Human plague is a severe and often fatal zoonotic disease caused by Yersinia pestis. For public health investigations of human cases, nonintensive whole genome molecular typing tools, capable of defining epidemiologic relationships, are advantageous. Whole genome multilocus sequence typing (wgMLST) is a recently developed methodology that simplifies genomic analyses by transforming millions of base pairs of sequence into character data for each gene. We sequenced 13 US Y. pestis isolates with known epidemiologic relationships. Sequences were assembled de novo, and multilocus sequence typing alleles were assigned by comparison against 3979 open reading frames from the reference strain CO92. Allele-based cluster analysis accurately grouped the 13 isolates, as well as 9 publicly available Y. pestis isolates, by their epidemiologic relationships. Our findings indicate wgMLST is a simplified, sensitive, and scalable tool for epidemiologic analysis of Y. pestis strains.

Clonal Clustering Using 10-Gene Multilocus Sequence Typing Reveals an Association Between Genotype and Listeria monocytogenes Maximum Growth Rate in Defined Medium

We used a 10-gene (10G) multilocus sequence typing scheme to investigate the diversity and phylogenetic distribution of 124 Listeria monocytogenes strains across major lineages, major serotypes, and seven epidemic clones that have been previously associated with outbreaks. The 124 isolates proved to be diverse, with a total of 81 sequence types (10G-STs) belonging to 13 clonal complexes (CCs), where all STs of the same CC differ from one another in up to 3 of the 10 alleles (named as 10G-triple-locus-variant-clonal-complexes [10G-TLV-CCs]). Phenotypic characterization for 105 of the 124 strains showed that L. monocytogenes had variable maximum growth rate (μ(max)) in a defined medium at 16°C, and classification by lineage or serotype was not able to reflect the genetic basis for the difference of this phenotype. Among the six major 10G-TLV-CCs, 10G-TLV-CC4 that included lineage I strains had significantly lower μ(max) (Tukey honestly significant difference adjusted [adj.] p < 0.05) compared to 10G-TLV-CC1 and 10G-TLV-CC3 that both comprised lineage II strains, indicating a distinct difference in growth of these L. monocytogenes isolates under nutrient-limited conditions among some of the CCs. However, the other three (10G-TLV-CC2, 6, and 10) of the six major 10G-TLV-CCs containing either lineage I or lineage II strains did not show significantly different μ(max) compared to the others (adj. p < 0.05). Our findings highlighted the importance of using molecular typing methods that can be used in evolutionary analyses as a framework for further understanding the phenotypic characteristics of subgroups of L. monocytogenes.

BacWGSTdb, a database for genotyping and source tracking bacterial pathogens

Whole genome sequencing has become one of the routine methods in molecular epidemiological practice. In this study, we present BacWGSTdb (http://bacdb.org/BacWGSTdb), a bacterial whole genome sequence typing database which is designed for clinicians, clinical microbiologists and hospital epidemiologists. This database borrows the population structure from the current multi-locus sequence typing (MLST) scheme and adopts a hierarchical data structure: species, clonal complex and isolates. When users upload the pre-assembled genome sequences to BacWGSTdb, it offers the functionality of bacterial genotyping at both traditional MLST and whole-genome levels. More importantly, users are told which isolates in the public database are phylogenetically close to the query isolate, along with their clinical information such as host, isolation source, disease, collection time and geographical location. In this way, BacWGSTdb offers a rapid and convenient platform for worldwide users to address a variety of clinical microbiological issues such as source tracking bacterial pathogens.