cgMLST Scheme Research Articles

A Whole-Genome-Based Gene-by-Gene Typing System for Standardized High-Resolution Strain Typing of Bacillus anthracis.

ABSTRACTWhole-genome sequencing (WGS) has been established for bacterial subtyping and is regularly used to study pathogen transmission, to investigate outbreaks, and to perform routine surveillance. Core-genome multilocus sequence typing (cgMLST) is a bacterial subtyping method that uses WGS data to provide a high-resolution strain characterization. This study aimed at developing a novel cgMLST scheme for Bacillus anthracis, a notorious pathogen that causes anthrax in livestock and humans worldwide. The scheme comprises 3,803 genes that were conserved in 57 B. anthracis genomes spanning the whole phylogeny. The scheme has been evaluated and applied to 584 genomes from 50 countries. On average, 99.5% of the cgMLST targets were detected. The cgMLST results confirmed the classical canonical single-nucleotide-polymorphism (SNP) grouping of B. anthracis into major clades and subclades. Genetic distances calculated based on cgMLST were comparable to distances from whole-genome-based SNP analysis with similar phylogenetic topology and comparable discriminatory power. Additionally, the application of the cgMLST scheme to anthrax outbreaks from Germany and Italy led to a definition of a cutoff threshold of five allele differences to trace epidemiologically linked strains for cluster typing and transmission analysis. Finally, the association of two clusters of B. anthracis with human cases of injectional anthrax in four European countries was confirmed using cgMLST. In summary, this study presents a novel cgMLST scheme that provides high-resolution strain genotyping for B. anthracis. This scheme can be used in parallel with SNP typing methods to facilitate rapid and harmonized interlaboratory comparisons, essential for global surveillance and outbreak analysis. The scheme is publicly available for application by users, including those with little bioinformatics knowledge.

Expanding of ST11 Carbapenemase-Producing Klebsiella pneumoniae Subclones in a Chinese Hospital, Shenzhen, China.

BackgroundST11 is the most prevalent sequence type of clinical Klebsiella pneumoniae in China.MethodsWe investigated the characteristics of the ST11 subclones using core genome multi-locus sequence typing (cgMLST). Ninety-three carbapenemase-producing K. pneumoniae isolates were collected at Shenzhen People’s Hospital. Then, whole-genome sequencing and cgMLST were used to discriminate apparent subclones within the ST11 group.ResultsWe analyzed the prevalence and genetic relationships of these subclones. ST11 and K. pneumoniae carbapenemase (KPC-2) were the predominant genotype and carbapenemase, respectively, in the clinical carbapenemase-producing K. pneumoniae strains. cgMLST scheme genotyping divided the ST11 group into two clades across seven complex types (CTs). CT1313 was the most prevalent subclone. The deletion of galF and a high frequency of SNPs in genes associated with the stress- and SOS-responses were found in CT1291 and CT2405 over time, respectively.ConclusionOur results indicated that the subclones of the ST11 group had different patterns of prevalence. Highly discriminatory genotyping techniques, such as cgMLST scheme, should be used in further molecular epidemiology investigations.

Development and evaluation of a core genome multilocus sequence typing scheme for Paenibacillus larvae, the deadly American foulbrood pathogen of honeybees

Summary Paenibacillus larvae is the causative agent of the fatal American foulbrood disease in honeybees (Apis mellifera). Strain identification is vital for preventing the spread of the disease. To date, the most accessible and robust scheme to identify strains is the multilocus sequence typing (MLST) method. However, this approach has limited resolution, especially for epidemiological studies. As the cost of whole‐genome sequencing has decreased and as it becomes increasingly available to most laboratories, an extended MLST based on the core genome (cgMLST) presents a valuable tool for high‐resolution investigations. In this study, we present a standardized, robust cgMLST scheme for P. larvae typing using whole‐genome sequencing. A total of 333 genomes were used to identify, validate and evaluate 2419 core genes. The cgMLST allowed fine‐scale differentiation between samples that had the same profile using traditional MLST and allowed for the characterization of strains impossible by MLST. The scheme was successfully used to trace a localized Swedish outbreak, where a cluster of 38 isolates was linked to a country‐wide beekeeping operation. cgMLST greatly enhances the power of a traditional typing scheme, while preserving the same stability and standardization for sharing results and methods across different laboratories.

Using core genome multilocus sequence typing (cgMLST) for vancomycin-resistant Enterococcus faecium isolates to guide infection control interventions and end an outbreak

ObjectivesUntil July 2016, vancomycin-resistantEnterococcus faecium (VREfm) was sporadically detected in Odense University Hospital, Denmark. After July 2016, the number of VREfm cases increased. This study aimed to apply a core genome multilocus sequence typing (cgMLST) scheme for E. faecium to type and analyse VREfm isolates collected at a single Danish hospital and to compare the results with cgMLST data from other regions of Denmark to trace transmission. MethodsA total of 38 VREfm clinical isolates from inpatients at the hospital in the period January 2014 through June 2017 were included in the study and analysed using whole-genome sequencing. Use of SeqSphere + software was initiated from the beginning of June 2017 to obtain MLST, cgMLST and epi curves. Admission histories were incorporated and national surveillance data on cgMLST were used to identify transmission routes. ResultsSix different sequence types (STs) were identified, the most frequent being ST80, ST117 and ST203. cgMLST subdivided the 38 isolates into 18 different complex types (CTs) with 13 isolates (34%) belonging to ST80-CT993. Epi curves indicated transmission of ST80-CT993 in several departments. Transmission from patients transferred from other hospitals was not identifiable. Infection control interventions launched in one department ended the outbreak. ConclusionThe high resolution of cgMLST allowed for detailed interpretation with evidence of nosocomial transmission of specific CTs. cgMLST made it easy to compare our local isolates with national findings, thereby clarifying transmission routes. Supplemented with admission histories, cgMLST targeted the epidemiological investigation and delineated the expensive and time-consuming infection control interventions.

Whole-Genome Sequencing Reveals Diversity of Carbapenem-Resistant Pseudomonas aeruginosa Collected Through the Emerging Infections Program

Background: Carbapenem-resistant Pseudomonas aeruginosa (CRPA) is a frequent cause of healthcare-associated infections (HAIs). The CDC Emerging Infections Program (EIP) conducted population and laboratory-based surveillance of CRPA in selected areas in 8 states from August 1, 2016, through July 31, 2018. We aimed to describe the molecular epidemiology and mechanisms of resistance of CRPA isolates collected through this surveillance. Methods: We defined a case as the first isolate of P. aeruginosa resistant to imipenem, meropenem, or doripenem from the lower respiratory tract, urine, wounds, or normally sterile sites identified from a resident of the EIP catchment area in a 30-day period; EIP sites submitted a systematic random sample of isolates to CDC for further characterization. Of 1,021 CRPA clinical isolates submitted, 707 have been sequenced to date using an Illumina MiSeq. Sequenced genomes were classified using the 7-gene multilocus sequence typing (MLST) scheme, and a core genome MLST (cgMLST) scheme was used to determine phylogeny. Antimicrobial resistance genes were identified using publicly available databases, and chromosomal mechanisms of carbapenem resistance were determined using previously validated genetic markers. Results: There were 189 sequence types (STs) among the 707 sequenced genomes (Fig. 1). The most frequently occurring were high-risk clones ST235 (8.5%) and ST298 (4.7%), which were found across all EIP sites. Carbapenemase genes were identified in 5 (<1%) isolates. Overall, 95.6% of the isolates had chromosomal mutations associated with carbapenem resistance: 93.2% had porinD-associated mutations that decrease membrane permeability to the drugs; 24.8% had mutations associated with overexpression of the multidrug efflux pump MexAB-OprM; and 22.9% had mutations associated with overexpression of the endogenous β-lactamase ampC. More than 1 such chromosomal resistance mutation type was present in 37.8% of the isolates. Conclusions: The diversity of the sequence types demonstrates that HAIs caused by CRPA can arise from a variety of strains and that high-risk clones are broadly disseminated across the EIP sites but are a minority of CRPA strains overall. Carbapenem resistance in P. aeruginosa was predominantly driven by chromosomal mutations rather than acquired mechanisms (ie, carbapenemases). The diversity of the CRPA isolates and the lack of carbapenemase genes suggest that this ubiquitous pathogen can readily evolve chromosomal resistance mechanisms, but unlike carbapenemases, these cannot be easily spread through horizontal transfer.Funding: NoneDisclosures: None

Development of a High-Resolution Single-Nucleotide Polymorphism Strain-Typing Assay Using Whole Genome-Based Analyses for the Lactobacillus acidophilus Probiotic Strain.

Lactobacillus acidophilus is one of the most commonly used industrial products worldwide. Since its probiotic efficacy is strain-specific, the identification of probiotics at both the species and strain levels is necessary. However, neither phenotypic nor conventional genotypic methods have enabled the effective differentiation of L. acidophilus strains. In this study, a whole-genome sequence-based analysis was carried out to establish high-resolution strain typing of 41 L. acidophilus strains (including commercial isolates and reference strains) using the cano-wgMLST_BacCompare analytics platform; consequently, a strain-specific discrimination method for the probiotic strain LA1063 was developed. Using a core-genome multilocus sequence-typing (cgMLST) scheme based on 1390 highly conserved genes, 41 strains could be assigned to 34 sequence types. Subsequently, we screened a set of 92 loci with a discriminatory power equal to that of the 1390 loci cgMLST scheme. A strain-specific polymerase chain reaction combined with a multiplex minisequencing method was developed based on four (phoU, secY, tilS, and uvrA_1) out of 21 loci, which could be discriminated between LA1063 and other L. acidophilus strains using the cgMLST data. We confirmed that the strain-specific single-nucleotide polymorphisms method could be used to quickly and accurately identify the L. acidophilus probiotic strain LA1063 in commercial products.

Development and Application of a Core Genome Multilocus Sequence Typing Scheme for the Health Care-Associated Pathogen Pseudomonas aeruginosa.

Pseudomonas aeruginosa is an opportunistic human pathogen that frequently causes health care-associated infections (HAIs). Due to its metabolic diversity and ability to form biofilms, this Gram-negative nonfermenting bacterium can persist in the health care environment, which can lead to prolonged HAI outbreaks. We describe the creation of a core genome multilocus sequence typing (cgMLST) scheme to provide a stable platform for the rapid comparison of P. aeruginosa isolates using whole-genome sequencing (WGS) data. We used a diverse set of 58 complete P. aeruginosa genomes to curate a set of 4,440 core genes found in each isolate, representing ∼64% of the average genome size. We then expanded the alleles for each gene using 1,991 contig-level genome sequences. The scheme was used to analyze genomes from four historical HAI outbreaks to compare the phylogenies generated using cgMLST to those of other means (traditional MLST, pulsed-field gel electrophoresis [PFGE], and single-nucleotide variant [SNV] analysis). The cgMLST scheme provides sufficient resolution for analyzing individual outbreaks, as well as the stability for comparisons across a variety of isolates encountered in surveillance studies, making it a valuable tool for the rapid analysis of P. aeruginosa genomes.

A Core Genome Multilocus Sequence Typing Scheme for Streptococcus mutans.

Streptococcus mutans is one of the primary pathogens responsible for the development of dental caries. Recent whole-genome sequencing (WGS)-based core genome multilocus sequence typing (cgMLST) approaches have been employed in epidemiological studies of specific human pathogens. However, this approach has not been reported in studies of S. mutans Here, we therefore developed a cgMLST scheme for S. mutans We surveyed 199 available S. mutans genomes as a means of identifying cgMLST targets, developing a scheme that incorporated 594 targets from the S. mutans UA159 reference genome. Sixty-eight sequence types (STs) were identified in this cgMLST scheme (cgSTs) in 80 S. mutans isolates from 40 children that were sequenced in this study, compared to 35 STs identified by multilocus sequence typing (MLST). Fifty-six cgSTs (82.35%) were associated with a single isolate based on our cgMLST scheme, which is significantly higher than in the MLST scheme (11.43%). In addition, 58.06% of all MLST profiles with ≥2 isolates were further differentiated by our cgMLST scheme. Topological analyses of the maximum likelihood phylogenetic trees revealed that our cgMLST scheme was more reliable than the MLST scheme. A minimum spanning tree of 145 S. mutans isolates from 10 countries developed based upon the cgMLST scheme highlighted the diverse population structure of S. mutans This cgMLST scheme thus offers a new molecular typing method suitable for evaluating the epidemiological distribution of this pathogen and has the potential to serve as a benchmark for future global studies of the epidemiological nature of dental caries.IMPORTANCEStreptococcus mutans is regarded as a major pathogen responsible for the onset of dental caries. S. mutans can transmit among people, especially within families. In this study, we established a new epidemiological approach to S. mutans classification. This approach can effectively differentiate among closely related isolates and offers superior reliability relative to that of the traditional MLST molecular typing method. As such, it has the potential to better support effective public health strategies centered around this bacterium that are aimed at preventing and treating dental caries.

Core Genome Multilocus Sequence Typing for Food Animal Source Attribution of Human Campylobacter jejuni Infections.

Campylobacter jejuni is a major foodborne pathogen and common cause of bacterial enteritis worldwide. A total of 622 C. jejuni isolates recovered from food animals and retail meats in the United States through the National Antimicrobial Resistance Monitoring System between 2013 and 2017 were sequenced using an Illumina MiSeq. Sequences were combined with WGS data of 222 human isolates downloaded from NCBI and analyzed by core genome multilocus sequence typing (cgMLST) and traditional MLST. cgMLST allelic difference (AD) thresholds of 0, 5, 10, 25, 50, 100 and 200 identified 828, 734, 652, 543, 422, 298 and 197 cgMLST types among the 844 isolates, respectively, and traditional MLST identified 174 ST. The cgMLST scheme allowing an AD of 200 (cgMLST200) revealed strong correlation with MLST. cgMLST200 showed 40.5% retail chicken isolates, 56.5% swine, 77.4% dairy cattle and 78.9% beef cattle isolates shared cgMLST sequence type with human isolates. All ST-8 had the same cgMLST200 type (cgMLST200-12) and 74.3% of ST-8 and 75% cgMLST200-12 were confirmed as sheep abortion virulence clones by PorA analysis. Twenty-nine acquired resistance genes, including 21 alleles of blaOXA, tetO, aph(3′)-IIIa, ant(6)-Ia, aadE, aad9, aph(2′)-Ig, aph(2′)-Ih, sat4 plus mutations in gyrA, 23SrRNA and L22 were identified. Resistance genotypes were strongly linked with cgMLST200 type for certain groups including 12/12 cgMLST200-510 with the A103V substitution in L22 and 10/11 cgMLST200-608 with the T86I GyrA substitution associated with macrolide and quinolone resistance, respectively. In summary, the cgMLST200 threshold scheme combined with resistance genotype information could provide an excellent subtyping scheme for source attribution of human C. jejuni infections.

A Core Genome Multilocus Sequence Typing Scheme for Pseudomonas aeruginosa.

Pseudomonas aeruginosa is a ubiquitous microorganism and an important opportunistic pathogen responsible for a broad spectrum of infections mainly in immunosuppressed and critically ill patients. Molecular investigations traditionally rely on pulsed field gel electrophoresis (PFGE) and multilocus sequence typing (MLST). In this work we propose a core genome multilocus sequence typing (cgMLST) scheme for P. aeruginosa, a methodology that combines traditional MLST principles with whole genome sequencing data. All publicly available complete P. aeruginosa genomes, representing the diversity of this species, were used to establish a cgMLST scheme targeting 2,653 genes. The scheme was then tested using genomes available at contig, chromosome and scaffold levels. The proposed cgMLST scheme for P. aeruginosa typed over 99% (2,314/2,325) of the genomes available for this study considering at least 95% of the cgMLST target genes present. The absence of a certain number gene targets at the threshold considered for both the creation and validation steps due to low genome sequence quality is possibly the main reason for this result. The cgMLST scheme was compared with previously published whole genome single nucleotide polymorphism analysis for the characterization of the population structure of the epidemic clone ST235 and results were highly similar. In order to evaluate the typing resolution of the proposed scheme, collections of isolates belonging to two important STs associated with cystic fibrosis, ST146 and ST274, were typed using this scheme, and ST235 isolates associated with an outbreak were evaluated. Besides confirming the relatedness of all the isolates, earlier determined by MLST, the higher resolution of cgMLST denotes that it may be suitable for surveillance programs, overcoming possible shortcomings of classical MLST. The proposed scheme is publicly available at: https://github.com/BioinformaticsHIAEMolecularMicrobiology/cgMLST-Pseudomonas-aeruginosa.

Genomic Characterization of Escherichia coli Isolates Belonging to a New Hybrid aEPEC/ExPEC Pathotype O153:H10-A-ST10 eae-beta1 Occurred in Meat, Poultry, Wildlife and Human Diarrheagenic Samples.

Different surveillance studies (2005–2015) in northwest Spain revealed the presence of eae-positive isolates of Escherichia coli O153:H10 in meat for human consumption, poultry farm, wildlife and human diarrheagenic samples. The aim of this study was to explore the genetic and genomic relatedness between human and animal/meat isolates, as well as the mechanism of its persistence. We also wanted to know whether it was a geographically restricted lineage, or whether it was also reported elsewhere. Conventional typing showed that 32 isolates were O153:H10-A-ST10 fimH54, fimAvMT78, traT and eae-beta1. Amongst these, 21 were CTX-M-32 or SHV-12 producers. The PFGE XbaI-macrorestriction comparison showed high similarity (>85%). The plasmidome analysis revealed a stable combination of IncF (F2:A-:B-), IncI1 (STunknown) and IncX1 plasmid types, together with non-conjugative Col-like plasmids. The core genome investigation based on the cgMLST scheme from EnteroBase proved close relatedness between isolates of human and animal origin. Our results demonstrate that a hybrid MDR aEPEC/ExPEC of the clonal group O153:H10-A-ST10 (CH11-54) is circulating in our region within different hosts, including wildlife. It seems implicated in human diarrhea via meat transmission, and in the spreading of ESBL genes (mainly of CTX-M-32 type). We found genomic evidence of a related hybrid aEPEC/ExPEC in at least one other country.

Core Genome Allelic Profiles of Clinical Klebsiella pneumoniae Strains Using a Random Forest Algorithm Based on Multilocus Sequence Typing Scheme for Hypervirulence Analysis.

Hypervirulent Klebsiella pneumoniae (hvKP) infections can have high morbidity and mortality rates owing to their invasiveness and virulence. However, there are no effective tools or biomarkers to discriminate between hvKP and nonhypervirulent K. pneumoniae (nhvKP) strains. We aimed to use a random forest algorithm to predict hvKP based on core-genome data. In total, 272 K. pneumoniae strains were collected from 20 tertiary hospitals in China and divided into hvKP and nhvKP groups according to clinical criteria. Clinical data comparisons, whole-genome sequencing, virulence profile analysis, and core genome multilocus sequence typing (cgMLST) were performed. We then established a random forest predictive model based on the cgMLST scheme to prospectively identify hvKP. The random forest is an ensemble learning method that generates multiple decision trees during the training process and each decision tree will output its own prediction results corresponding to the input. The predictive ability of the model was assessed by means of area under the receiver operating characteristic curve. Patients in the hvKP group were younger than those in the nhvKP group (median age, 58.0 and 68.0 years, respectively; P < .001). More patients in the hvKP group had underlying diabetes mellitus (43.1% vs 20.1%; P < .001). Clinically, carbapenem-resistant K. pneumoniae was less common in the hvKP group (4.1% vs 63.8%; P < .001), whereas the K1/K2 serotype, sequence type (ST) 23, and positive string tests were significantly higher in the hvKP group. A cgMLST-based minimal spanning tree revealed that hvKP strains were scattered sporadically within nhvKP clusters. ST23 showed greater genome diversification than did ST11, according to cgMLST-based allelic differences. Primary virulence factors (rmpA, iucA, positive string test result, and the presence of virulence plasmid pLVPK) were poor predictors of the hypervirulence phenotype. The random forest model based on the core genome allelic profile presented excellent predictive power, both in the training and validating sets (area under receiver operating characteristic curve, 0.987 and 0.999 in the training and validating sets, respectively). A random forest algorithm predictive model based on the core genome allelic profiles of K. pneumoniae was accurate to identify the hypervirulent isolates.

Core/Whole Genome Multilocus Sequence Typing and Core Genome SNP-Based Typing of OXA-48-Producing Klebsiella pneumoniae Clinical Isolates From Spain.

Whole-genome sequencing (WGS)-based typing methods have emerged as promising and highly discriminative epidemiological tools. In this study, we combined gene-by-gene allele calling and core genome single nucleotide polymorphism (cgSNP) approaches to investigate the genetic relatedness of a well-characterized collection of OXA-48-producing Klebsiella pneumoniae isolates. We included isolates from the predominant sequence type ST405 (n = 31) OXA-48-producing K. pneumoniae clone and isolates from ST101 (n = 3), ST14 (n = 1), ST17 (n = 1), and ST1233 (n = 1), obtained from eight Catalan hospitals. Core-genome multilocus sequence typing (cgMLST) schemes from Institut Pasteur’s BIGSdb-Kp (634 genes) and SeqSphere+ (2,365 genes), and a SeqSphere+ whole-genome MLST (wgMLST) scheme (4,891 genes) were used. Allele differences or allelic mismatches and the genetic distance, as the proportion of allele differences, were used to interpret the results from a gene-by-gene approach, whereas the number of SNPs was used for the cgSNP analysis. We observed between 0–10 and 0–14 allele differences among the predominant ST405 using cgMLST and wgMLST from SeqSphere+, respectively, and <2 allelic mismatches when using Institut Pasteur’s BIGSdb-Kp cgMLST scheme. For ST101, we observed 14 and 54 allele differences when using cgMLST and wgMLST SeqSphere+, respectively, and 2–5 allelic mismatches for BIGSdb-Kp cgMLST. A low genetic distance (<0.0035, a previously established threshold for epidemiological link) was generally in concordance with a low number of allele differences (<8) when using the SeqSphere+ cgMLST scheme. The cgSNP analysis showed 6–29 SNPs in isolates with identical allelic SeqSphere+ cgMLST profiles and 16–61 cgSNPs among ST405 isolates. Furthermore, comparison of WGS-based typing results with previously obtained MLST and pulsed-field gel electrophoresis (PFGE) data showed some differences, demonstrating the different molecular principles underlying these techniques. In conclusion, the use of the different WGS-based typing methods that were used to elucidate the genetic relatedness of clonal OXA-48-producing K. pneumoniae all led to the same conclusions. Furthermore, threshold parameters in WGS-based typing methods should be applied with caution and should be used in combination with clinical epidemiological data and population and species characteristics.

Core Genome Multi-locus Sequence Typing Analyses of Leptospira spp. Using the Bacterial Isolate Genome Sequence Database.

With the advent of whole-genome sequencing (WGS), comparative analysis has led to the use of core genome MLST (cgMLST) schemes for the high-resolution reproducible typing of bacterial isolates. In cgMLST, hundreds of loci are used for gene-by-gene comparisons of assembled genomes for studying the genetic diversity of clinically important pathogens. Combination of the cgMLST data and metadata of the isolates is useful for epidemiological investigations.Here we present a cgMLST scheme for the high-resolution typing of isolates from the whole Leptospira genus, enabling identification at the level of species, clades, clonal groups, and sequence types. We show several examples how the cgMLST Leptospira database, which is a publicly available web-based database, can be used for the analyses of WGS data of Leptospira isolates. This effort was undertaken in order to facilitate international collaborations and support the global surveillance of leptospirosis.

Tracing local and regional clusters of carbapenemase-producing Klebsiella pneumoniae ST512 with whole genome sequencing, Finland, 2013 to 2018.

BackgroundTwo epidemiologically-unrelated clusters of Klebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae were detected among several healthcare facilities (HCF) in Finland by routine surveillance using whole genome sequencing (WGS).AimThe objective was to investigate transmission chains to stop further spread of the responsible strain.MethodsIn this observational retrospective study, cases were defined as patients with K. pneumoniae KPC-3 sequence type (ST)512 strain detected in Finland from August 2013 to May 2018. Environmental specimens were obtained from surfaces, sinks and toilets in affected wards. WGS was performed on K. pneumoniae cultures using Illumina MiSeq platform and data were analysed using Ridom SeqShere software K. pneumoniae core genome multilocus sequence typing (cgMLST) scheme. Epidemiological information of the cases was provided by HCFs.ResultsWe identified 20 cases in six HCFs: cluster 1 included 18 cases in five HCFs and cluster 2 two cases in one HCF. In cluster 1, a link with a foreign country was unclear, 6/18 cases without overlapping stay had occupied the same room in one of the five HCFs within > 3 years. In cluster 2, the index case was transferred from abroad, both cases occupied the same room 8 months apart. A strain identical to that of the two cases in cgMLST was isolated from the toilet of the room, suggesting a clonal origin.ConclusionsThe clusters were mostly related to case transfer between facilities and likely involved environmental transmission. We show that CPE surveillance using WGS and collaboration between hospitals are crucial to identify clusters and trace transmission chains.

In vitro characterization and genetic diversity of Bordetella avium field strains

ABSTRACT Bordetella avium (BA) is a respiratory pathogen of particular importance for turkeys. Specific adherence and damage to the respiratory epithelia are crucial steps of the pathogenesis, but knowledge about the mechanisms and the variety of virulence in field strains is limited. We analysed 17 BA field strains regarding their in vitro virulence-associated properties in tracheal organ cultures (TOC) of turkey embryos, and their genetic diversity. The TOC adherence assay indicated that BA field strains differ considerably in their ability to adhere to the tracheal mucosa, while the TOC ciliostasis assay illustrated a high degree of diversity in ciliostatic effects. These two virulence-associated properties were associated with each other in the investigated strains. Three of the investigated strains displayed significantly (P > 0.05) lower in vitro virulence in comparison to other strains. Genetic diversity of BA strains was analysed by core genome multilocus sequence typing (cgMLST). We applied a cgMLST scheme comprising 2667 targets of the reference genome (77.3% of complete genome, BA strain 197N). The results showed a broad genetic diversity in BA field strains but did not demonstrate a correlation between sequence type and virulence-associated properties. The cgMLST analysis revealed that strains with less marked virulence-associated properties had a variety of mutations in the putative filamentous haemagglutinin gene. Likewise, amino acid sequence alignment indicated variations in the protein. The results from our study showed that both adherence and ciliostasis assay can be used for virulence characterization of BA. Variations in the filamentous haemagglutinin protein may be responsible for reduced virulence of BA field strains.

Whole Genome Sequencing Based Surveillance of L. monocytogenes for Early Detection and Investigations of Listeriosis Outbreaks.

In Austria, all laboratories are legally obligated to forward human and food/environmental L. monocytogenes isolates to the National Reference Laboratory/Center (NRL) for Listeria. Two invasive human isolates of L. monocytogenes serotype 1/2a of the same pulsed-field gel electrophoresis (PFGE) pattern, previously unknown in Austria, were cultured for the first time in January 2016. Five further human isolates, obtained from patients with invasive listeriosis between April 2016 and September 2017, showed this PFGE pattern. In Austria the NRL started to use whole-genome sequencing (WGS) based typing in 2016, using a core genome MLST (cgMLST) scheme developed by Ruppitsch et al. 2015, which contains 1701 target genes. Sequence data are submitted to a publicly available nomenclature server (Ridom GmbH, Münster, Germany) for allocation of the core genome complex type (CT). The seven invasive human isolates differed from each other with zero to two alleles and were allocated to CT1234 (declared as outbreak strain). Among the Austrian strain collection of about 6,000 cgMLST-characterized non-human isolates (i.e., food/environmental isolates) 90 isolates shared CT1234. Out of these, 83 isolates were traced back to one meat processing-company. They differed from the outbreak strain by up to seven alleles; one isolate originated from the company's industrial slicer. The remaining seven CT1234-isolates were obtained from food products of four other companies (five fish-products, one ready-to-eat dumpling and one deer-meat) and differed from the outbreak strain by six to eleven alleles. The outbreak described shows the considerable potential of WGS to identify the source of a listeriosis outbreak. Compared to PFGE analysis, WGS-based typing has higher discriminatory power, yields better data accuracy, and allows higher laboratory through-put at lower cost. Utilization of WGS-based typing results of human and food/ environmental L. monocytogenes isolates by appropriate public health analysts and epidemiologists is indispensable to support a successful outbreak investigation.

Genetic diversity and delineation of Salmonella Agona outbreak strains by next generation sequencing, Bavaria, Germany, 1993 to 2018.

BackgroundIn 2017, a food-borne Salmonella Agona outbreak caused by infant milk products from a French supplier occurred in Europe. Simultaneously, S. Agona was detected in animal feed samples in Bavaria.AimUsing next generation sequencing (NGS) and three data analysis methods, this study’s objectives were to verify clonality of the Bavarian feed strains, rule out their connection to the outbreak, explore the genetic diversity of Bavarian S. Agona isolates from 1993 to 2018 and compare the analysis approaches employed, for practicality and ability to delineate outbreaks caused by the genetically monomorphic Agona serovar.MethodsIn this observational retrospective study, three 2017 Bavarian feed isolates were compared to a French outbreak isolate and 48 S. Agona isolates from our strain collections. The later included human, food, feed, veterinary and environmental isolates, of which 28 were epidemiologically outbreak related. All isolates were subjected to NGS and analysed by: (i) a publicly available species-specific core genome multilocus sequence typing (cgMLST) scheme, (ii) single nucleotide polymorphism phylogeny and (iii) an in-house serovar-specific cgMLST scheme. Using additional international S. Agona outbreak NGS data, the cluster resolution capacity of the two cgMLST schemes was assessed.ResultsWe could prove clonality of the feed isolates and exclude their relation to the French outbreak. All approaches confirmed former Bavarian epidemiological clusters.ConclusionEven for S. Agona, species-level cgMLST can produce reasonable resolution, being standardisable by public health laboratories. For single samples or homogeneous sample sets, higher resolution by serovar-specific cgMLST or SNP genotyping can facilitate outbreak investigations.

Lineage-specific evolution in Listeria monocytogenes detected by analysis of a panel of Swiss isolates from food and human origin

Listeria monocytogenes is the causative agent of Listeriosis, a foodborne infection and characterised by outbreaks with high mortality. In the past years, outbreak analysis has shifted from PFGE to core genome MLST (cgMLST) and a difference of ≤10 cgMLST alleles was identified as the most appropriate cut off. The cgMLST scheme contributes significantly to the analysis of outbreaks but is less suitable to provide insight in the overall population structure. In the present work, we characterize the differences in the evolution of lineage I and II and the impact they have on the identification of outbreak isolates starting from the genome sequences of a panel of 166 Swiss L. monocytogenes clinical and food isolates. We included in the analysis recently published genomes from Germany (414 isolates) and Holland (128 isolates). Using data of pairwise cgMLST and SNP differences of these 708 isolates, we can clearly identify the genetic diversity associated to outbreaks (≤10 differences), sublineages (≤150 differences) and lineages. The sublineages identified by SNP and cgMLST match the clustering of the PopPUNK software. When broken down for lineages, data show clearly that lineage II has a lower mutation rate within sublineages but a higher diversity between sublineages. Admixture analysis confirms that the increased diversity between sublineages in lineage II is due to horizontal gene transfer and recombination. While our data show that genome evolution is lineage-specific in L. monocytogenes, the cut off for the identification of outbreaks remains identical between both lineages.

A Core Genome Multilocus Sequence Typing Scheme for Enterococcus faecalis.

Among enterococci, Enterococcus faecalis occurs ubiquitously, with the highest incidence of human and animal infections. The high genetic plasticity of E. faecalis complicates both molecular investigations and phylogenetic analyses. Whole-genome sequencing (WGS) enables unraveling of epidemiological linkages and putative transmission events between humans, animals, and food. Core genome multilocus sequence typing (cgMLST) aims to combine the discriminatory power of classical multilocus sequence typing (MLST) with the extensive genetic data obtained by WGS. By sequencing a representative collection of 146 E. faecalis strains isolated from hospital outbreaks, food, animals, and colonization of healthy human individuals, we established a novel cgMLST scheme with 1,972 gene targets within the Ridom SeqSphere+ software. To test the E. faecalis cgMLST scheme and assess the typing performance, different collections comprising environmental and bacteremia isolates, as well as all publicly available genome sequences from the NCBI and SRA databases, were analyzed. In more than 98.6% of the tested genomes, >95% good cgMLST target genes were detected (mean, 99.2% target genes). Our genotyping results not only corroborate the known epidemiological background of the isolates but exceed previous typing resolution. In conclusion, we have created a powerful typing scheme, hence providing an international standardized nomenclature that is suitable for surveillance approaches in various sectors, linking public health, veterinary public health, and food safety in a true One Health fashion.