Whole-genome Multilocus Sequence Typing Research Articles

Dissemination of extensively drug-resistant NDM-producing Providencia stuartii in Europe linked to patients transferred from Ukraine, March 2022 to March 2023

Background The war in Ukraine led to migration of Ukrainian people. Early 2022, several European national surveillance systems detected multidrug-resistant (MDR) bacteria related to Ukrainian patients. Aim To investigate the genomic epidemiology of New Delhi metallo-β-lactamase (NDM)-producing Providencia stuartii from Ukrainian patients among European countries. Methods Whole-genome sequencing of 66 isolates sampled in 2022–2023 in 10 European countries enabled whole-genome multilocus sequence typing (wgMLST), identification of resistance genes, replicons, and plasmid reconstructions. Five bla NDM-1-carrying-P. stuartii isolates underwent antimicrobial susceptibility testing (AST). Transferability to Escherichia coli of a bla NDM-1-carrying plasmid from a patient strain was assessed. Epidemiological characteristics of patients with NDM-producing P. stuartii were gathered by questionnaire. Results wgMLST of the 66 isolates revealed two genetic clusters unrelated to Ukraine and three linked to Ukrainian patients. Of these three, two comprised bla NDM-1-carrying-P. stuartii and the third bla NDM-5-carrying-P. stuartii. The bla NDM-1 clusters (PstCluster-001, n = 22 isolates; PstCluster-002, n = 8 isolates) comprised strains from seven and four countries, respectively. The bla NDM-5 cluster (PstCluster-003) included 13 isolates from six countries. PstCluster-001 and PstCluster-002 isolates carried an MDR plasmid harbouring bla NDM-1, bla OXA-10, bla CMY-16, rmtC and armA, which was transferrable in vitro and, for some Ukrainian patients, shared by other Enterobacterales. AST revealed PstCluster-001 isolates to be extensively drug-resistant (XDR), but susceptible to cefiderocol and aztreonam–avibactam. Patients with data on age (n = 41) were 19–74 years old; of 49 with information on sex, 38 were male. Conclusion XDR P. stuartii were introduced into European countries, requiring increased awareness and precautions when treating patients from conflict-affected areas.

Comparative genomics reveals high genetic similarity among strains of Salmonella enterica serovar Infantis isolated from multiple sources in Brazil.

Salmonella enterica serovar Infantis (Salmonella Infantis) is a zoonotic, ubiquitous and foodborne pathogen of worldwide distribution. Despite Brazil's relevance as a major meat exporter, few studies were conducted to characterize strains of this serovar by genomic analyses in this country. Therefore, this study aimed to assess the diversity of 80 Salmonella Infantis strains isolated from veterinary, food and human sources in Brazil between 2013 and 2018 by comparative genomic analyses. Additional genomes of non-Brazilian countries (n=18) were included for comparison purposes in some analyses. Analyses of whole-genome multi-locus sequence typing (wgMLST), using PGAdb-builder, and of fragmented genomes, using Gegenees, were conducted to compare the 80 Brazilian strains to the 18 non-Brazilian genomes. Pangenome analyses and calculations were performed for all Salmonella Infantis genomes analyzed. The presence of prophages was determined using PHASTER for the 80 Brazilian strains. The genome plasticity using BLAST Ring Image Generator (BRIG) and gene synteny using Mauve were evaluated for 20 selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries. Unique orthologous protein clusters were searched in ten selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries. wgMLST and Gegenees showed a high genomic similarity among some Brazilian Salmonella Infantis genomes, and also the correlation of some clusters with non-Brazilian genomes. Gegenees also showed an overall similarity >91% among all Salmonella Infantis genomes. Pangenome calculations revealed an open pangenome for all Salmonella Infantis subsets analyzed and a high gene content in the core genomes. Fifteen types of prophages were detected among 97.5% of the Brazilian strains. BRIG and Mauve demonstrated a high structural similarity among the Brazilian and non-Brazilian isolates. Unique orthologous protein clusters related to biological processes, molecular functions, and cellular components were detected among Brazilian and non-Brazilian genomes. The results presented using different genomic approaches emphasized the significant genomic similarity among Brazilian Salmonella Infantis genomes analyzed, suggesting wide distribution of closely related genotypes among diverse sources in Brazil. The data generated contributed to novel information regarding the genomic diversity of Brazilian and non-Brazilian Salmonella Infantis in comparison. The different genetically related subtypes of Salmonella Infantis from Brazil can either occur exclusively within the country, or also in other countries, suggesting that some exportation of the Brazilian genotypes may have already occurred.

Flagella Phenotypic Variations of ST34 Type Salmonella Typhimurium and Variants.

Salmonella enterica serovar Typhimurium and its variants are the most common serotypes of human salmonellosis cases. Serotyping Salmonella Typhimurium and its variants has always been challenging. Our previous work found that among 14 Salmonella Typhimurium and variant strains, some different antigenic formulas had 100% pulsed-field gel electrophoresis (PFGE) similarity. The 14 strains were sorted into 3 groups; in each group, the different antigenic formulas had the same PFGE patterns. This phenomenon suggested that different antigenic formula identification might originate from a common ancestor subtyped by PFGE. To assess whether the serotyping method on Salmonella Typhimurium and variant strains reflected the genetic relationship, we improved the discrimination for the phylogenetic relationship among the 14 Salmonella Typhimurium and variant strains using Fourier-transform infrared spectroscopy (FTIR) and whole-genome multilocus sequence typing (wgMLST). We compared the wgMLST assay of 14 Salmonella Typhimurium and variant strains from this study with 50 public ST34 strain data of Salmonella Typhimurium and variant strains. We also compared flagella (H antigen)-related genes based on the whole genome of 14 strains and the other 293 ST34 public database for further understanding of this question. The phylogenetic results (PFGE) showed no regularity between the antigenic formulas and genotypes. The results of the higher discrimination power assays (FTIR and whole-genome multilocus sequence typing) also showed no regularity between the antigenic formulas and genotypes (or phenotypes). The 58 flagella encoding genes of different antigenic formulas were sorted into 13 patterns. However, a similar phenomenon was found: the same flagella encoding gene patterns could express different antigenic formulas. In conclusion, there is no consistency between the antigenic formulas and phylogenetic relationships among ST34 Salmonella Typhimurium and variant strains, even in flagella antigenic formula and flagella encoding genes.

Immersion of debrided diabetic foot ulcer tissue in electrochemically generated pH neutral hypochlorous acid significantly reduces the microbial bioburden: whole-genome sequencing of Staphylococcus aureus, the most prevalent species recovered

Diabetic foot ulcer infections (DFUIs) are the leading cause of lower-limb amputations, mediated predominantly by Staphylococcus aureus. pH-neutral electrochemically generated hypochlorous acid (anolyte) is a non-toxic, microbiocidal agent with significant potential for wound disinfection. To investigate both the effectiveness of anolyte for microbial bioburden reduction in debrided ulcer tissues and the population of resident S.aureus. Fifty-one debrided tissues from 30 people with type II diabetes were aliquoted by wet weight and immersed in 1- or 10-mL volumes of anolyte (200parts per million) or saline for 3min. Microbial loads recovered were determined in colony forming units/g (cfu/g) of tissue following aerobic, anaerobic and staphylococcal-selective culture. Bacterial species were identified and 50S.aureus isolates from 30 tissues underwent whole-genome sequencing (WGS). The ulcers were predominantly superficial, lacking signs of infection (39/51, 76.5%). Of the 42/51 saline-treated tissues yielding ≥105cfu/g, a microbial threshold reported to impede wound-healing, only 4/42 (9.5%) were clinically diagnosed DFUIs. Microbial loads from anolyte-treated tissues were significantly lower than saline-treated tissues using 1mL (1065-fold, 2.0 log) and 10mL (8216-fold, 2.1 log) immersion volumes (P<0.0005). S.aureus was the predominant species recovered (44/51, 86.3%) and 50 isolates underwent WGS. All were meticillin susceptible and comprised 12 sequence types (STs), predominantly ST1, ST5 and ST15. Whole-genome multi-locus sequence typing identified three clusters of closely related isolates from 10 patients indicating inter-patient transmission. Short immersions of debrided ulcer tissue in anolyte significantly reduced microbial bioburden: a potential novel DFUI treatment.

Clostridioides difficile infections were predominantly driven by fluoroquinolone-resistant Clostridioides difficile ribotypes 176 and 001 in Slovakia in 2018–2019

To investigate the epidemiology of Clostridioides difficile infection (CDI) in Slovakian hospitals after the emergence of ribotype 176 (027-like) in 2016. Between 2018 and 2019, European Centre for Disease Control and Prevention CDI surveillance protocol v2.3 was applied to 14 hospitals, with additional data collected on recent antimicrobial use and the characterization of C. difficile isolates. The mean hospital incidence of CDI was 4.1 cases per 10,000 patient bed-days. One hundred and five (27.6%) in-hospital deaths were reported among the 381 cases. Antimicrobial treatment within the previous 4 weeks was recorded in 90.5% (333/368) of cases. Ribotype (RT)176 was detected in 50% (n=185/370, 14 hospitals) and RT001 was detected in 34.6% (n=128/370,13/14 hospitals) of cases with RT data. Overall, 86% (n=318/370) of isolates were resistant to moxifloxacin by Thr82Ile in GyrA (99.7%). Multi-locus variable tandem repeat analysis showed clonal relatedness of predominant RTs within and between hospitals. Seven of 14 sequenced RT176 isolates and five of 13 RT001 isolates showed between zero and three allelic differences by whole-genome multi-locus sequence typing. The majority of sequenced isolates (24/27) carried the erm(B) gene and 16/27 also carried the aac(6')-aph(2'') gene with the corresponding antimicrobial susceptibility phenotypes. Nine RT176 strains carried the cfr(E)gene and one RT001 strain carried the cfr(C) gene, but without linezolid resistance. The newly-predominant RT176 and endemic RT001 are driving the epidemiology of CDI in Slovakia. In addition to fluoroquinolones, the use of macrolide-lincosamide-streptogramin B antibiotics can represent another driving force for the spread of these epidemic lineages. In C. difficile, linezolid resistance should be confirmed phenotypically in strains with detected cfr gene(s).

Efflux Pumps and Different Genetic Contexts of tet(X4) Contribute to High Tigecycline Resistance in Escherichia fergusonii from Pigs.

Tigecycline is a last-resort antibiotic for the treatment of infections caused by multidrug-resistant bacteria. The emergence of plasmid-mediated tigecycline resistance genes is posing a serious threat to food safety and human health and has attracted worldwide attention. In this study, we characterized six tigecycline-resistant Escherichia fergusonii strains from porcine nasal swab samples collected from 50 swine farms in China. All the E. fergusonii isolates were highly resistant to tigecycline with minimal inhibitory concentration (MIC) values of 16-32 mg/L, and all contained the tet(X4) gene. In addition, 13-19 multiple resistance genes were identified in these isolates, revealed by whole-genome sequencing analysis. The tet(X4) gene was identified as being located in two different genetic structures, hp-abh-tet(X4)-ISCR2 in five isolates and hp-abh-tet(X4)-ΔISCR2-ISEc57-IS26 in one isolate. The role of efflux pumps in tigecycline resistance was evaluated by using inhibitor carbonyl cyanide 3-chlorophenylhydrazone (CCCP). The MIC values of tigecycline showed a 2- to 4-fold reduction in the presence of CCCP, indicating the involvement of active efflux pumps in tigecycline resistance in E. fergusonii. The tet(X4) gene was found to be transferable to Escherichia coli J53 by conjugation and resulted in the acquisition of tigcycline resistances in the transconjugants. Whole-genome multilocus sequence typing (wgMLST) and phylogenetic analysis showed a close relationship of five isolates originating from different pig farms, suggesting the transmission of tet(X4)-positive E. fergusonii between farms. In conclusion, our findings suggest that E. fergusonii strains in pigs are reservoirs of a transferable tet(X4) gene and provide insights into the tigecycline resistance mechanism as well as the diversity and complexity of the genetic context of tet(X4) in E. fergusonii.

2048. Molecular Characterization and Resistance Factors of Circulating Acinetobacter baumannii isolates in South-East Michigan

Abstract Background Carbapenem-resistant Acinetobacter baumannii (CRAb) is increasing due to widespread use of antibiotics. Multidrug resistant (MDR) CRAb is a major threat to public health as treatment options are limited. The objective of this study is to elucidate the molecular epidemiology of circulating antibiotic resistance genes causing MDR CRAb infections by using a combination of whole-genome Multi-Locus Sequence Typing (wgMLST) and antibiotic susceptibility phenotyping. Table 1.Molecular characterization of MDR CRAb isolates*The numbers indicate % of sequence identity match for each Beta-lactamase gene. Methods Bacterial isolates were derived from cultures taken from subjects 48 hours following admission as part of routine clinical care for patients between 2017-2020. Isolates were obtained from 16 hospital units (both ICU and non-ICU) across two hospitals in the Detroit area. Whole Genome Sequencing (WGS) was performed using Illumina MiniSeq or Nextseq. WgMLST analysis was performed using BioNumerics software v7.6. ResFinder software was used for analysis of antibiotic resistance genes. Isolates underwent antibiotic susceptibility testing using a broth microdilution method (VITEK2) and Clinical &amp; Laboratory Standards Institute (CLSI) minimum inhibitory concentration (MIC) cut offs. Results Out of the 95 total isolates, 51(54%) were CRAb isolates and of the CRAb isolates, 21(41%) were MDR CRAb. WgMLST identified that majority of the circulating MDR CRAb isolates belonged to ST2Pas (ST195Ox and ST208Ox) based on CDC definitions (Table 1). MDR CRAb isolates were resistant to 3 different classes of antibiotics including aminoglycosides, fluroquinolones and β-lactams. β-lactamase genes present include (blaADC-25, blaOXA-23, blaOXA-66 and blaTEM1D) for both ST195Ox and ST208Ox and (blaADC-25, blaOXA-23 and blaOXA-223) for ST406Pas (ST310Ox). Among the patients with MDR CRAb infections, most were males with respiratory infections in a non-ICU setting. Conclusion The study demonstrated a high proportion of isolates belonged to ST2 Pas carrying multiple beta-lactamase genes including blaOXA-23 gene. ST406Pas might be an emerging lineage carrying the blaOXA-23 gene. In addition to stringent infection control measures, continuous surveillance is recommended in limiting the spread of MDR CRAb isolates in the healthcare settings. Disclosures Chetan Jinadatha, MD, MPH, AHRQ R01 Grant-5R01HS025598: Grant/Research Support|EOS Surfaces: Copper Coupons and materials for testing Keith S. Kaye, MD, MPH, Allecra: Advisor/Consultant|GlaxoSmithKline plc.: Receiving symposia honoraria|GlaxoSmithKline plc.: GlaxoSmithKline plc.-sponsored study 212502|Merck: Advisor/Consultant|qpex: Advisor/Consultant|Shionogi: Grant/Research Support|Spero: Advisor/Consultant Piyali Chatterjee, PhD, AHRQ Grant # 1R03HS027667-01: Grant/Research Support|AHRQ Grant # 1R03HS027667-01: Central Texas Veterans Health Care System.

Comparative Analysis of Salmonella enterica subsp. enterica Serovar Thompson Isolates associated with Outbreaks Using PFGE and wgMLST

The strains associated with foodborne Salmonella enterica Thompson outbreaks in Korea have not been identified. Therefore, we characterized S. Thompson strains isolated from chocolate cakes linked to foodborne outbreaks in Korea. A total of 56 strains were isolated from preserved cake products, products in the supply chain distribution, the manufacturer's apparatus, and egg white liquid products used for cream preparation. Subsequently, serological typing, pathogenic gene-targeted polymerase chain reaction (PCR), pulsed-field gel electrophoresis (PFGE), and whole-genome multi-locus sequence typing (wgMLST) were performed to characterize these isolates. The antigen formula of all isolates was 7:k:1,5, namely Salmonella enterica subsp. enterica Serovar Thompson. All 56 isolates harbored invA, his, hin, and stn, and were negative for sefA and spvC based on gene-targeted PCR analyses. Based on PFGE results, these isolates were classified into one group based on the same SP6X01.011 pattern with 100% similarity. We selected 19 strains based on the region and sample type, which were subjected to wgMLST. Although the examined strains showed 100% similarity, they were classified into seven clusters based on allelic differences. According to our findings, the cause of these outbreaks was chocolate cake manufactured with egg white liquid contaminated with the same Salmonella Thompson. Additionally, comparative analysis of wgMLST on domestic isolates of S. Thompson from the three outbreaks showed genetic similarities of over 99.6%. Based on the results, the PFGE and wgMLST combination can provide highly resolved phylogeny and reliable evidence during Salmonella outbreak investigations.

Transmission Patterns of Clostridioides difficile in a Non-Epidemic Setting Based on WGS Analysis

Clostridioides difficile is a major nosocomial pathogen and has a considerable burden on healthcare systems. Our objective was to determine the transmission patterns of C. difficile in a non-epidemic setting using whole-genome multi-locus sequence typing (wgMLST) and core-genome single-nucleotide polymorphism (cgSNP) analyses. A retrospective study was conducted in a 650-bed university hospital between January 2016 and February 2017. In total, 191 strains isolated from 169 symptomatic C. difficile infection (CDI) patients were analyzed by WGS. Sequences were compared using wgMLST and cgSNP analyses. Genetic data and ward movements were then combined to identify the transmission rate and the type of transmission. The transmission rate varied from 55/169 (19.5%) (wgMLST) to 33/169 (32.5%) (cgSNP). Most transmission was considered cryptic, irrespective of the genetic analysis (38/55 [69.1%] by wgMLST to 25/33 [75.8%] by cgSNP). No transmission within the same ward was observed. In a non-epidemic setting, most C. difficile transmission occurs from sources other than symptomatic CDI patients.

A prospective matched case-control study on the genomic epidemiology of colistin-resistant Enterobacterales from Dutch patients

BackgroundColistin is a last-resort treatment option for infections with multidrug-resistant Gram-negative bacteria. However, colistin resistance is increasing.MethodsA six-month prospective matched case-control study was performed in which 22 Dutch laboratories with 32 associated hospitals participated. Laboratories were invited to send a maximum of five colistin-resistant Escherichia coli or Klebsiella pneumoniae (COLR-EK) isolates and five colistin-susceptible isolates (COLS-EK) to the reference laboratory, matched for patient location, material of origin and bacterial species. Epidemiological/clinical data were collected and included in the analysis. Characteristics of COLR-EK/COLS-EK isolates were compared using logistic regression with correction for variables used for matching. Forty-six ColR-EK/ColS-EK pairs were analysed by next-generation sequencing (NGS) for whole-genome multi-locus sequence typing and identification of resistance genes, including mcr genes. To identify chromosomal mutations potentially leading to colistin resistance, NGS reads were mapped against gene sequences of pmrAB, phoPQ, mgrB and crrB.ResultsIn total, 72 COLR-EK/COLS-EK pairs (75% E. coli and 25% K. pneumoniae) were included. Twenty-one percent of COLR-EK patients had received colistin, in contrast to 3% of COLS-EK patients (OR > 2.9). Of COLR-EK isolates, five contained mcr-1 and two mcr-9. One isolate lost mcr-9 after repeated sub-culturing, but retained colistin resistance. Among 46 sequenced COLR-EK isolates, genetic diversity was large and 19 (41.3%) isolates had chromosomal mutations potentially associated with colistin resistance.ConclusionsColistin resistance is present but uncommon in the Netherlands and caused by the mcr gene in a minority of COLR-EK isolates. There is a need for surveillance of colistin resistance using appropriate susceptibility testing methods.

Whole-genome sequence analysis of Clostridium chauvoei isolated from clinical case of black quarter (BQ) from India.

Black quarter (BQ) is an infectious disease affecting cattle and small ruminants worldwide caused by Gram-positive anaerobic bacterium Clostridium chauvoei. In this study, a draft genome sequence of C. chauvoei NIVEDIBQ1 strain isolated from clinical case of black quarter was analyzed. Sequence analysis indicated that genome had 2653 predicted coding DNA sequences, harbored numerous genes, mobile genetic elements for pathogenesis, and virulence factors. Computational analysis revealed that strain contained 30 virulence-associated genes. An intact genomic region highly similar to the Clostridium phage was present in the genome. Presence of CRISPR systems and the transposon components likely contribute to the genome plasticity. Strain encode diverse spectrum of degradative carbohydrate-active enzymes (CAZymes). Comparative SNP analysis revealed that the genomes of the C. chauvoei strains analyzed were highly conserved. Phylogenetic analysis of strains and available genome (n = 21) based on whole-genome multi-locus sequence typing (wgMLST) and core orthologous genes showed the clustering of strains into two different clusters suggesting geographical links.

Annotated Whole-Genome Multilocus Sequence Typing Schema for Scalable High-Resolution Typing of Streptococcus pyogenes.

Streptococcus pyogenes is a major human pathogen with high genetic diversity, largely created by recombination and horizontal gene transfer, making it difficult to use single nucleotide polymorphism (SNP)-based genome-wide analyses for surveillance. Using a gene-by-gene approach on 208 complete genomes of S. pyogenes, a novel whole-genome multilocus sequence typing (wgMLST) schema was developed, comprising 3,044 target loci. The schema was used for core-genome MLST (cgMLST) analyses of previously published data sets and 265 newly sequenced draft genomes with other molecular and phenotypic typing data. Clustering based on cgMLST data supported the genetic heterogeneity of many emm types and correlated poorly with pulsed-field gel electrophoresis macrorestriction profiling, superantigen gene profiling, and MLST sequence type, highlighting the limitations of older typing methods. While 763 loci were present in all isolates of a data set representative of S. pyogenes genetic diversity, the proposed schema allows scalable cgMLST analysis, which can include more loci for an increased resolution when typing closely related isolates. The cgMLST and PopPUNK clusters were broadly consistent in this diverse population. The cgMLST analyses presented results comparable to those of SNP-based methods in the identification of two recently emerged sublineages of emm1 and emm89 and the clarification of the genetic relatedness among isolates recovered in outbreak contexts. The schema was thoroughly annotated and made publicly available on the chewie-NS online platform (https://chewbbaca.online/species/1/schemas/1), providing a framework for high-resolution typing and analyzing the genetic variability of loci of particular biological interest.

Identifying the Sources of Intestinal Colonization With Extended-Spectrum β-Lactamase-Producing Escherichia coli in Healthy Infants in the Community.

The prevalence of fecal colonization with extended-spectrum β-lactamase-producing Escherichia coli (ESBL-Ec) among children in low- and middle-income countries is alarmingly high. This study aimed to identify the sources of ESBL-Ec colonization in children < 1 year old through comparative analysis of E. coli isolates from child stool, child’s mother stool, and point-of-use drinking water from 46 rural households in Bangladesh. The pairwise similarity in antibiotic susceptibility of E. coli from all three sources was evaluated, followed by phylogenetic clustering using enterobacterial repetitive intergenic consensus polymerase chain reaction and whole-genome sequence analysis of the isolates. Matching antibiotic susceptibility and enterobacterial repetitive intergenic consensus polymerase chain reaction patterns were found among ESBL-Ec isolates from child–mother dyads of 24 and 11 households, respectively, from child–water dyads of 5 and 4 households, respectively, and from child–mother–water triads of 3 and 4 households, respectively. Whole-genome sequence analysis of 30 isolates from 10 households revealed that ESBL-Ec from children in five households (50%) was clonally related to ESBL-Ec either from their mothers (2 households), drinking water sources (2 households), or both mother and drinking-water sources (1 household) based on serotype, phylogroup, sequence type, antibiotic resistance genes, mobile genetic elements, core single-nucleotide polymorphisms, and whole-genome multilocus sequence typing. Overall, this study provides empirical evidence that ESBL-Ec colonization in children is linked to the colonization status of mothers and exposure to the household environments contaminated with ESBL-Ec. Interventions such as improved hygiene practices and a safe drinking water supply may help reduce the transmission of ESBL-Ec at the household level.

Whole-genome sequencing for the characterization of resistance mechanisms and epidemiology of colistin-resistant Acinetobacter baumannii.

BackgroundMultidrug-resistant Acinetobacter baumannii is an important causal pathogen of healthcare-associated infections, and colistin-resistant strains have recently emerged owing to the increased use of colistin. Using next-generation sequencing (NGS), a single whole-genome sequencing (WGS) protocol can identify and type pathogens, analyze genetic relationships among different pathogens, predict pathogenic transmissions, and detect antibiotic resistance genes. However, only a few studies have applied NGS in studying the resistance mechanism and epidemiology of colistin-resistant A. baumannii. This study aimed to elucidate the resistance mechanism of colistin-resistant A. baumannii and analyze its molecular epidemiology through WGS.Materials and methodsThe subjects in this study were patients who visited a university hospital between 2014 and 2018. Thirty colistin-resistant strains with high minimum inhibitory concentrations were selected from various patient samples, and WGS was performed. Comparative genomic analysis was performed for the 27 colistin-resistant A. baumannii strains using a colistin-susceptible strain as the reference genome.ResultsThe WGS analysis found no mutation for lpxA, lpxC, lpx D, pmrA, pmrB, and mcr1, the genes known to be associated with colistin resistance. Fifty-seven coding sequences (CDS) showed differences; they included 13 CDS with known names and functions that contained 21 genes. From the whole-genome multi-locus sequence typing (wgMLST) and single nucleotide polymorphism (SNP) analyses, two major clusters were found for the colistin-resistant A. baumannii strains. However, no differences were observed by the time of detection for each cluster, the samples, the pattern of antibiotic resistance, or the patient characteristics. In the conventional MLST following the Oxford scheme, the typing result showed ST1809, ST451, ST191, ST1837, and ST369 in the global clone 2 (GC2), without any relation with the results of wgMLST and SNP analyses.ConclusionBased on the findings of the resistance gene analysis through WGS and comparative genomic analysis, the potential genes associated with colistin-resistance or CDS were examined. Furthermore, the analysis of molecular epidemiology through WGS regarding colistin-resistant A. baumannii may prove helpful in preventing infection by multidrug-resistant bacteria and controlling healthcare-associated infections.

Using whole-genome sequencing to assess the diversity of Paenibacillus larvae within an outbreak and a beekeeping operation.

The spore-forming bacterium Paenibacillus larvae is the causative agent of American foulbrood (AFB), a devastating disease of honeybees (Apis mellifera). In the present study, we used whole-genome sequencing (WGS) to investigate an extensive outbreak of AFB in northwestern Slovenia in 2019. A total of 59 P . larvae isolates underwent WGS, of which 40 originated from a single beekeeping operation, to assess the diversity of P. larvae within the beekeeping operation, apiary and colony. By applying a case-specific single-linkage threshold of 34 allele differences (AD), whole-genome multilocus sequence typing (wgMLST) identified two outbreak clusters represented by ERIC II-ST11 clones. All isolates from a single beekeeping operation fell within cluster 1 and the median pairwise AD between them was 10 (range=1–22). The median pairwise AD for apiaries of the same beekeeping operation ranged from 8 to 11 (min.=1, max.=22). For colonies of the same apiary and honey samples from these colonies, the median pairwise AD ranged from 8 to 14 (min.=1, max.=20). The maximum within-cluster distance was 33 pairwise AD for cluster 1 and 44 for cluster 2 isolates. The minimum distance between the outbreak-related and non-related isolates was 37 AD, confirming the importance of associated epidemiological data for delineating outbreak clusters. The observed transmission events could be explained by the activities of honeybees and beekeepers. The present study provides insight into the genetic diversity of P. larvae at different levels and thus provides information for future AFB surveillance.

Molecular Identification and Selection of Probiotic Strains Able to Reduce the Serum TMAO Level in Mice Challenged with Choline.

Trimethylamine oxide (TMAO) originates from trimethylamine (TMA), which is oxidized in the liver by hepatic flavin-containing monooxygenases (FMO3). TMA is produced by its dietary precursors such as choline, carnitine, and phosphatidylcholine by gut microbiota. TMAO attracts attention, identified as a novel and independent risk factor for promoting obesity, atherosclerosis and cardiovascular disease (CVD), chronic kidney disease (CKD), insulin tolerance, and colon cancer. Probiotics have been considered as live microorganisms, providing benefits to their host when they are given in sufficient quantities and administered continuously. The objective of this study is to suggest a method to select potential probiotic strains to reduce the serum concentration of TMAO in mice fed with choline. In this work, we chose three lactobacilli with strong adherence capability, and fed multistrain formula (MF) to the mice challenged with choline. On days 7, 14, and day 28, it was found that the MF-containing L. amylovorus LAM1345, Lpb. plantarum LP1145, and Lim. fermentum LF33 showed a significant reduction in serum TMAO and TMA levels. For the single strains, LP1145 reduced TMAO on days 14 and 28, and strain LAM1345 reduced TMAO significantly on days 7 and day 14. For strain LF1143 from strain LF33, it showed no significant effect on TMAO and TMA. Thus, MF showed the best effect, which may be due to the additive and synergetic effect and the contribution of strain LP1145 and LAM1345. Finally, for the LAM1345 and LP1145 strains, we used molecular identification and typing methods to assure that these two strains are unique strains. The methods used for LAM 1345 were leader peptidase A (lepA) gene analysis and phylogenetic analysis, while for strain LP 1145and other strains of Lpb. plantarum subsp. plantarum sequences were compared using the whole-genome multilocus sequence typing (wgMLST) method.

A machine learning-based typing scheme refinement for Listeria monocytogenes core genome multilocus sequence typing with high discriminatory power for common source outbreak tracking.

BackgroundAs whole-genome sequencing for pathogen genomes becomes increasingly popular, the typing methods of gene-by-gene comparison, such as core genome multilocus sequence typing (cgMLST) and whole-genome multilocus sequence typing (wgMLST), are being routinely implemented in molecular epidemiology. However, some intrinsic problems remain. For example, genomic sequences with varying read depths, read lengths, and assemblers influence the genome assemblies, introducing error or missing alleles into the generated allelic profiles. These errors and missing alleles might create “specious discrepancy” among closely related isolates, thus making accurate epidemiological interpretation challenging. In addition, the rapid growth of the cgMLST allelic profile database can cause problems related to storage and maintenance as well as long query search times.MethodsWe attempted to resolve these issues by decreasing the scheme size to reduce the occurrence of error and missing alleles, alleviate the storage burden, and improve the query search time. The challenge in this approach is maintaining the typing resolution when using fewer loci. We achieved this by using a popular artificial intelligence technique, XGBoost, coupled with Shapley additive explanations for feature selection. Finally, 370 loci from the original 1701 cgMLST loci of Listeria monocytogenes were selected.ResultsAlthough the size of the final scheme (LmScheme_370) was approximately 80% lower than that of the original cgMLST scheme, its discriminatory power, tested for 35 outbreaks, was concordant with that of the original cgMLST scheme. Although we used L. monocytogenes as a demonstration in this study, the approach can be applied to other schemes and pathogens. Our findings might help elucidate gene-by-gene–based epidemiology.

Comparison of Whole Genome Sequencing and Repetitive Element PCR for Multidrug-Resistant Pseudomonas aeruginosa Strain Typing

Hospital-acquired infections pose significant costly global challenges to patient care. Rapid and sensitive methods to identify potential outbreaks are integral to infection control measures. Whole-genome sequencing (WGS)-based bacterial strain typing provides higher discriminatory power over standard nucleotide banding pattern-based methods such as repetitive sequence-based PCR (rep-PCR). However, integration of WGS into clinical epidemiology is limited by the lack of consensus in methodology and data analysis/interpretation. In this study, WGS was performed on genomic DNA extracted from 22 multidrug-resistant Pseudomonas aeruginosa (MDR-PA) isolates using next-generation sequencing. Resulting high-quality reads were analyzed for phylogenetic relatedness using a whole-genome multilocus sequence typing (wgMLST)-based software program and single-nucleotide variant phylogenomics (SNVPhyl). WGS-based results were compared with conventional MLST and archived rep-PCR results. Rep-PCR identified three independent clonal clusters of MDR-PA. Only one clonal cluster identified by rep-PCR, an endemic strain within the pediatric cystic fibrosis population at Texas Children's Hospital, was concordantly identified using wgMLST and SNVPhyl. Results were highly consistent between the three sequence-based analyses (conventional MLST, wgMLST, and SNVPhyl), and these results remained consistent with the addition of 74 MDR-PA genomes. These WGS-based methods provided greater resolution for strain discrimination than rep-PCR or standard MLST classification, and the ease of use of wgMLST software renders it clinically viable for analysis, interpretation, and reporting of WGS-based strain typing.

Evaluation of whole-genome sequencing-based subtyping methods for the surveillance of Shigella spp. and the confounding effect of mobile genetic elements in long-term outbreaks.

Many public health laboratories across the world have implemented whole-genome sequencing (WGS) for the surveillance and outbreak detection of foodborne pathogens. PulseNet-affiliated laboratories have determined that most single-strain foodborne outbreaks are contained within 0–10 multi-locus sequence typing (MLST)-based allele differences and/or core genome single-nucleotide variants (SNVs). In addition to being a food- and travel-associated outbreak pathogen, most Shigella spp. cases occur through continuous person-to-person transmission, predominantly involving men who have sex with men (MSM), leading to long-term and recurrent outbreaks. Continuous transmission patterns coupled to genetic evolution under antibiotic treatment pressure require an assessment of existing WGS-based subtyping methods and interpretation criteria for cluster inclusion/exclusion. An evaluation of 4 WGS-based subtyping methods [SNVPhyl, coreMLST, core genome MLST (cgMLST) and whole-genome MLST (wgMLST)] was performed on 9 foodborne-, travel- and MSM-related retrospective outbreaks from a collection of 91 Shigella flexneri and 232 Shigella sonnei isolates to determine the methods’ epidemiological concordance, discriminatory power, robustness and ability to generate stable interpretation criteria. The discriminatory powers were ranked as follows: coreMLST<SNVPhyl<cgMLST<wgMLST (range: 0.970–1.000). The genetic differences observed for non-MSM-related Shigella spp. outbreaks respect the standard 0–10 allele/SNV guideline; however, mobile genetic element (MGE)-encoded loci caused inflated genetic variation and discrepant phylogenies for prolonged MSM-related S. sonnei outbreaks via wgMLST. The S. sonnei correlation coefficients of wgMLST were also the lowest at 0.680, 0.703 and 0.712 for SNVPhyl, coreMLST and cgMLST, respectively. Plasmid maintenance, mobilization and conjugation-associated genes were found to be the main source of genetic distance inflation in addition to prophage-related genes. Duplicated alleles arising from the repeated nature of IS elements were also responsible for many false cg/wgMLST differences. The coreMLST approach was shown to be the most robust, followed by SNVPhyl and wgMLST for inter-laboratory comparability. Our results highlight the need for validating species-specific subtyping methods based on microbial genome plasticity and outbreak dynamics in addition to the importance of filtering confounding MGEs for cluster detection.

Whole-genome sequencing and ad hoc shared genome analysis of Staphylococcus aureus isolates from a New Zealand primary school

Epidemiological studies of communicable diseases increasingly use large whole-genome sequencing (WGS) datasets to explore the transmission of pathogens. It is important to obtain an initial overview of datasets and identify closely related isolates, but this can be challenging with large numbers of isolates and imperfect sequencing. We used an ad hoc whole-genome multi locus sequence typing method to summarise data from a longitudinal study of Staphylococcus aureus in a primary school in New Zealand. Each pair of isolates was compared and the number of genes where alleles differed between isolates was tallied to produce a matrix of “allelic differences”. We plotted histograms of the number of allelic differences between isolates for: all isolate pairs; pairs of isolates from different individuals; and pairs of isolates from the same individual. 340 sequenced isolates were included, and the ad hoc shared genome contained 445 genes. There were between 0 and 420 allelic differences between isolate pairs and the majority of pairs had more than 260 allelic differences. We found many genetically closely related S. aureus isolates from single individuals and a smaller number of closely-related isolates from separate individuals. Multiple S. aureus isolates from the same individual were usually very closely related or identical over the ad hoc shared genome. Siblings carried genetically similar, but not identical isolates. An ad hoc shared genome approach to WGS analysis can accommodate imperfect sequencing of the included isolates, and can provide insights into relationships between isolates in epidemiological studies with large WGS datasets containing diverse isolates.