Within-host genomic evolution of methicillin-resistant Staphylococcus aureus in long-term carriers

Assessing the genomic evolution of Staphylococcus aureus can help us understand how the bacteria adapt to its environment. In this study, we aimed to assess the mutation rate within 144 methicillin-resistant Staphylococcus aureus (MRSA) carriers with a carriage time from 4 to 11 years, including some carriers who belonged to the same households. We found that 23 of the 144 individuals had completely different MRSA types over time and were therefore not long-term carriers of the same MRSA. From the remaining 121 individuals, we performed whole-genome sequencing (WGS) on 424 isolates and then compared these pairwise using core genome multilocus sequence typing (cgMLST) and single-nucleotide polymorphism (SNP) analyses. We found a median within-host mutation rate in long-term MRSA carriers of 4.9 (3.4-6.9) SNPs/genome/year and 2.7 (1.8-4.2) allelic differences/genome/year, when excluding presumed recombination. Furthermore, we stratified the cohort into subgroups and found no significant difference between the median mutation rate of members of households, individuals with presumed continued exposure, e.g., from travel and persons without known continued exposure. Finally, we found that SNPs occurred at random within the genes in our cohort. KEY POINTS: • Median mutation rate within long-term MRSA carriers of 4.9 (3.4-6.9) SNPs/genome/year • Similar median mutation rates in subgroups (households, travelers) • No hotspots for SNPs within the genome.

A novel cgMLST for genomic surveillance of Yersinia enterocolitica infections in France allowed the detection and investigation of outbreaks in 2017-2021

Enteric yersiniosis, the third most common food-borne zoonosis in Europe, is mainly caused by the pathogen Yersinia enterocolitica. In France, the yersiniosis microbiological surveillance is conducted at the Yersinia National Reference Laboratory (YNRL). Since 2017, isolates have been characterized by whole genome sequencing (WGS) followed by a 500-gene Yersinia-cgMLST. We report here the data of the WGS-based surveillance on Y. enterocolitica isolates for the 2017-2021 period. The YNRL characterized 7,642 Y. enterocolitica strains distributed in 2,497 non-pathogenic isolates from lineages 1Aa and 1Ab, and 5,145 specimens belonging to 8 pathogenic lineages. Among pathogenic isolates, lineage 4 was the most common (87.2%) followed by lineages 2/3-9b (10.6%), 2/3-5a (1.2%), 2/3-9a (0.6%), 3-3b, 3-3c, 1B, and 3-3d (0.1% per each). Importantly, we developed a routine surveillance system based on a new typing method consisting of a 1,727-genes core genome Multilocus Sequence Typing (cgMLST) specific to the species Y. enterocolitica followed by isolate clustering. Thresholds of allelic distances (AD) were determined and fixed for the clustering of isolates: AD ≤ 5 for lineages 4, 2/3-5a, and 2/3-9a, and AD ≤ 3 for lineage 2/3-9b. Clustering programs were implemented in 2019 in routine surveillance to detect genomic clusters of pathogenic isolates. In total, 419 clusters with at least 2 isolates were identified, representing 2,504 of the 3,503 isolates characterized between 2019 and 2021. Most clusters (n = 325) comprised 2 to 5 isolates. The new typing method proved to be useful for the molecular investigation of unusual grouping of cases as well as for the detection of genomic clusters in routine surveillance.

IMPORTANCE

We describe here the new typing method used for molecular surveillance of Yersinia enterocolitica infections in France based on a novel core genome Multilocus Sequence Typing (cgMLST) specific to Y. enterocolitica species. This method can reliably identify the pathogenic Y. enterocolitica subspecies and compare the isolates with a high discriminatory power. Between 2017 and 2021, 5,145 pathogenic isolates belonging to 8 lineages were characterized and lineage 4 was by far the most common followed by lineage 2/3-9b. A clustering program was implemented, and detection thresholds were cross-validated by the molecular and epidemiological investigation of three unusual groups of Y. enterocolitica infections. The routine molecular surveillance system has been able to detect genomic clusters, leading to epidemiological investigations.

Recurrent, ICD-associated L. monocytogenes bacteraemia with multiple septic pulmonary embolisms over a 2-year period

BACKGROUND

Listeria monocytogenes is a bacterial pathogen known for causing listeriosis, a foodborne illness with a wide spectrum of clinical presentations ranging from mild gastroenteritis to severe invasive disease, particularly affecting immunocompromised individuals, pregnant women, newborns, and the elderly. Successful treatment of patients with recurring listeria episodes due to colonised foreign material is often challenging, typically requiring a combination of antimicrobial treatment and surgical removal.

CASE PRESENTATION

Here, we present a particularly complex case of chronic invasive listeriosis with a total of six relapses. After extensive investigations, the patient's ICD device was identified as the focus of infection.

CONCLUSION

The confirmation of relapses through cgMLST analysis highlights the persistence of Listeria monocytogenes and the potential for recurrence even after apparent resolution of symptoms in patients with foreign material. It emphasises the necessity for a comprehensive assessment to identify and mitigate the risk of relapses, thereby ensuring optimal management and outcomes.

The first vanE-type vancomycin resistant Enterococcus faecalis isolates in Norway - phenotypic and molecular characteristics

OBJECTIVES

We aimed to characterize the vanE cluster and its genetic support in the first Norwegian vanE-type isolates and assess genetic relatedness to other vanE isolates.

METHODS

Two vanE-type vancomycin resistant Enterococcus faecalis (vanE-VREfs) isolates (E1 and E2) recovered from the same patient 30 months apart were examined for antimicrobial susceptibility, genome sequence, vancomycin resistance induction, vanE transferability, genome mutation rate, and phylogenetic relationship to E. faecalis closed genomes and two vanE-VREfs from North America.

RESULTS

The ST34 E1 and E2 strains expressed low-level vancomycin resistance and susceptibility to teicoplanin. Their vanE gene clusters were part of a non-transferable Tn6202. The histidine kinase part of vanSE was expressed although a premature stop codon (E1) and insertion of a transposase (E2) truncated their vanSE gene. The vancomycin resistance phenotype in E1 was inducible while constitutive in E2. E1 showed a 125-fold higher mutation rate than E2. Variant calling showed 60 variants but nearly identical chromosomal gene content and synteny between the isolates. Their genomes also showed high similarity to another ST34 vanE-VREfs from Canada.

CONCLUSION

In-depth genomic analyses of the first two vanE-VREfs found in Europe identified a single chromosomal insertion site of two variants of vanE-conferring Tn6202. Single nucleotide polymorphism (SNP) and core genome multilocus sequence type (cgMLST) analyses show the genomes are different. This can be explained by the high mutation rate of E1 and acquisition of different mobile genetic elements; thus, we believe the two isolates from the same patient are genetically related. Genome similarities also suggest relatedness between the Canadian and Norwegian vanE-VREfs.

Introduction and spread of vancomycin-resistant Enterococcus faecium (VREfm) at a German tertiary care medical center from 2004 until 2010: a retrospective whole-genome sequencing (WGS) study of the molecular epidemiology of VREfm

BACKGROUND

In most of Europe and especially in Germany, there is currently a concerning rise in the number of hospital-acquired infections due to vancomycin-resistant Enterococcus faecium (VREfm). Therefore, there is a need to improve our understanding of the way VREfm spreads in hospitals. In this study, we investigated the molecular epidemiology of VREfm isolates from the first appearance at our university hospital in 2004 until 2010. There is only very scarce information about the molecular epidemiology of VREfm from this early time in Germany.

METHODS

Our analysis includes all available first VREfm isolates of each patient at our tertiary care center collected during the years 2004-2010. If available, additional consecutive VREfm isolates from some patients were analyzed. We used multilocus sequence typing (MLST) and core genome multilocus sequence typing (cgMLST) for the analysis and description of nosocomial transmission pathways as well as the detection of outbreaks.

RESULTS

VREfm isolates from 158 patients and 76 additional subsequent patient isolates were included in the analysis. Until 2006, detections of VREfm remained singular cases, followed by a peak in the number of VREfm cases in 2007 and 2008 with a subsequent decline to baseline in 2010. MLST and cgMLST analysis show significant changes in the dominant sequence types (STs) and complex types (CTs) over the study period, with ST192 and ST17 being responsible for the peak in VREfm cases in 2007 and 2008. The four largest clusters detected during the study period are comprised of these two STs. Cluster analysis shows a focus on specific wards and departments for each cluster. In the early years of this study (2004-2006), all analyzed VREfm stemmed from clinical specimens, whereas since 2007, approximately half of the VREfm were detected by screening. Of the 234 VREfm isolates analyzed, 96% had a vanB and only 4% had a vanA resistance genotype.

CONCLUSIONS

This retrospective study contributes significant knowledge about regional VREfm epidemiology from this early VREfm period in Germany. One remarkable finding is the striking dominance of vanB-positive VREfm isolates over the entire study period, which is in contrast with countrywide data. Analysis of cgMLST shows the transition from sporadic VRE cases at our institution to a sharp increase in VRE numbers triggered by oligoclonal spread and specific outbreak clusters with the dominance of ST192 and ST17.

Genotyping of Campylobacter jejuni and prediction tools of its antimicrobial resistance

Although Campylobacter jejuni is the pathogen responsible for the most common foodborne illness, tracing of the infection source remains challenging due to its highly variable genome. Therefore, one of the aim of the study was to compare three genotyping methods (MLST, PFGE, and mP-BIT) to determine the most effective genotyping tool. C. jejuni strains were divided into 4 clusters based on strain similarity in the cgMLST dendrogram. Subsequently, the dendrograms of the 3 tested methods were compared to determine the accuracy of each method compared to the reference cgMLST method. Moreover, a cost-benefit analysis has showed that MLST had the highest inverse discrimination index (97%) and required less workflow, time, fewer consumables, and low bacterial sample quantity. PFGE was shown to be obsolete both because of its low discriminatory power and the complexity of the procedure. Similarly, mP‑BIT showed low separation results, which was compensated by its high availability. Therefore, our data showed that MLST is the optimal tool for genotyping C. jejuni. Another aim was to compare the antimicrobial resistance to ciprofloxacin, erythromycin, and tetracycline in C. jejuni strains isolated from human, water, air, food, and animal samples by two gene sequence-based prediction methods and to compare them with the actual susceptibility of C. jejuni strains using the disc diffusion method. Both tools, ResFinder and RGI, synchronously predict the antimicrobial susceptibility of C. jejuni and either can be used.

Performance evaluation of core genome multilocus sequence typing for genotyping of Mycobacterium tuberculosis strains in China: based on multicenter, population-based collection

PURPOSE

To evaluate the performance of core genome multilocus sequence typing (cgMLST) for genotyping Mycobacterium tuberculosis (M.tuberculosis) Strains in regions where the lineage 2 strains predominate.

METHODS

We compared clustering by whole-genome SNP typing with cgMLST clustering in the analysis of WGS data of 6240 strains from five regions of China. Using both the receiver operating characteristic (ROC) curve and epidemiological investigation to determine the optimal threshold for defining genomic clustering by cgMLST. The performance of cgMLST was evaluated by quantifying the sensitivity, specificity and concordance of clustering between two methods. Logistic regression was used to gauge the impact of strain genetic diversity and lineage on cgMLST clustering.

RESULTS

The optimal threshold for cgMLST to define genomic clustering was determined to be ≤ 10 allelic differences between strains. The overall sensitivity and specificity of cgMLST averaged 99.6% and 96.3%, respectively; the concordance of clustering between two methods averaged 97.1%. Concordance was significantly correlated with strain genetic diversity and was 3.99 times (95% CI, 2.94-5.42) higher in regions with high genetic diversity (π > 1.55 × 10-4) compared to regions with low genetic diversity. The difference missed statistical significance, while concordance for lineage 2 strains (96.8%) was less than that for lineage 4 strains (98.3%). CONCLUSION : cgMLST showed a discriminatory power comparable to whole-genome SNP typing and could be used to genotype clinical M.tuberculosis strains in different regions of China. The discriminative power of cgMLST was significantly correlated with strain genetic diversity and was slightly lower with strains from regions with low genetic diversity.

Antimicrobial Resistance in Enterococcus spp. Isolates from Red Foxes (Vulpes vulpes) in Latvia

Antimicrobial resistance (AMR) is an emerging public health threat and is one of the One Health priorities for humans, animals, and environmental health. Red foxes (Vulpes vulpes) are a widespread predator species with great ecological significance, and they may serve as a sentinel of antimicrobial resistance in the general environment. The present study was carried out to detect antimicrobial resistance, antimicrobial resistance genes, and genetic diversity in faecal isolates of red foxes (Vulpes vulpes). In total, 34 Enterococcus isolates, including E. faecium (n = 17), E. faecalis (n = 12), E. durans (n = 3), and E. hirae (n = 2), were isolated. Antimicrobial resistance to 12 antimicrobial agents was detected with EUVENC panels using the minimum inhibitory concentration (MIC). The presence of antimicrobial resistance genes (ARGs) was determined using whole-genome sequencing (WGS). Resistance to tetracycline (6/34), erythromycin (3/34), ciprofloxacin (2/34), tigecycline (2/34), and daptomycin (2/34) was identified in 44% (15/34) of Enterococcus isolates, while all the isolates were found to be susceptible to ampicillin, chloramphenicol, gentamicin, linezolid, teicoplanin, and vancomycin. No multi-resistant Enterococcus spp. were detected. A total of 12 ARGs were identified in Enterococcus spp., with the presence of at least 1 ARG in every isolate. The identified ARGs encoded resistance to aminoglycosides (aac(6')-I, ant(6)-Ia, aac(6')-Iih and spw), tetracyclines (tet(M), tet(L) and tet(S)), and macrolide-lincosamide-streptogramin AB (lnu(B,G), lsa(A,E), and msr(C)), and their presence was associated with phenotypical resistance. Core genome multilocus sequence typing (cgMLST) revealed the high diversity of E. faecalis and E. faecium isolates, even within the same geographical area. The distribution of resistant Enterococcus spp. in wild foxes in Latvia highlights the importance of a One Health approach in tackling AMR.

Antibiotic resistance, plasmids, and virulence-associated markers in human strains of Campylobacter jejuni and Campylobacter coli isolated in Italy

Campylobacteriosis, a prevalent foodborne gastrointestinal infection in Europe, is primarily caused by Campylobacter jejuni and Campylobacter coli, with rising global concerns over antimicrobial resistance in these species. This study comprehensively investigates 133 human-origin Campylobacter spp. strains (102 C. jejuni and 31 C. coli) collected in Italy from 2013 to 2021. The predominant Multilocus Sequence Typing Clonal complexes (CCs) were ST-21 CC and ST-206 CC in C. jejuni and ST-828 CC in C. coli. Ciprofloxacin and tetracycline resistance, mainly attributed to GyrA (T86I) mutation and tet(O) presence, were prevalent, while erythromycin resistance was associated with 23S rRNA gene mutation (A2075G), particularly in C. coli exhibiting multidrug-resistant pattern CipTE. Notable disparities in virulence factors among strains were observed, with C. jejuni exhibiting a higher abundance compared to C. coli. Notably, specific C. jejuni sequence types, including ST-21, ST-5018, and ST-1263, demonstrated significantly elevated counts of virulence genes. This finding underscores the significance of considering both the species and strain-level variations in virulence factor profiles, shedding light on potential differences in the pathogenicity and clinical outcomes associated with distinct C. jejuni lineages. Campylobacter spp. plasmids were classified into three groups comprising pVir-like and pTet-like plasmids families, exhibiting diversity among Campylobacter spp. The study underscores the importance of early detection through Whole Genome Sequencing to identify potential emergent virulence, resistance/virulence plasmids, and new antimicrobial resistance markers. This approach provides actionable public health data, supporting the development of robust surveillance programs in Italy.

Investigating environmental transmission to resolve a Bacillus cereus group outbreak in a neonatal intensive care unit using core genome multilocus sequence typing

BACKGROUND

We analyzed an outbreak of Bacillus cereus group (Bcg) at a single-center neonatal intensive care unit level IV by conducting comprehensive sampling of both patients and the environment.

METHODS

Between 06/2020 and 10/2021, all Bcg isolates identified by both regular colonization screening and additional sampling of the environment were subjected to whole-genome sequencing, followed by in vitro extraction of MLST ST, resistance genes and virulence factors. Using publicly available genome sequences, we defined an ad hoc core genome multilocus sequence typing (cgMLST) scheme comprising 2759 target genes for Bcg typing, which we applied to the detected isolates. We have compared the results with a stable cgMLST that was published in the meantime and completed the investigation with a SNP analysis.

RESULTS

We analyzed 28 Bcg isolates from patient and environmental samples using MLST and cgMLST. This revealed multiple sequence types, with ST127 being the most common (n = 13). Both cgMLST schemes grouped ten of the 13 ST127 isolates into a cluster, including two invasive isolates from two different patients and several environmental samples. SNP analysis postulated a screen from a ventilation machine as a possible reservoir.

CONCLUSION

In sensitive settings such as neonatal intensive care units, considering the environment in outbreak analyses is crucial, especially when investigating potential transmission routes through shared devices. When dealing with widespread bacteria such as Bcg, high-resolution typing techniques are necessary. In this study, we successfully resolved an outbreak of Bcg infections using a custom cgMLST scheme combined with a SNP analysis.

Molecular characterization of Listeria monocytogenes strains isolated from imported food in China from 14 countries/regions, 2003-2018

Listeria monocytogenes (Lm) is associated with severe foodborne infections and ubiquitous in the nature. Identification of characteristics of Lm transmission through trading of food products is essential for rapidly tracking Lm sources and controlling dissemination of listeriosis. In this study, a total of 44 Lm strains were isolated from food products originating from 14 countries/regions during 2003-2018 at the Shanghai port. The genomes of these Lm strains were sequenced by high-throughput sequencing. Multilocus sequence typing (MLST) analysis showed that 43 isolates were divided into 17 sequence types (STs). The distribution of STs was decentralized, with the dominant ST2 accounting for only 18.18% of the strains. The LM63 strain did not match with any of the existing STs. Core-genome MLST (cgMLST) analysis based on 1748 core genes categorized the 44 strains into 30 cgMLST types (CTs), with CT10153 and CT7892 as the most predominant CTs. Notably, LM63 and LM67 shared the same CT in the cgMLST analysis. The phylogenetic analysis based on single-copy homologous genes revealed that the 44 Lm strains were primarily classified into two lineages. The SNP analysis also indicated that these strains were roughly divided into two clades, with strains in the first clade mainly collected earlier than those in the second clade, which were predominantly collected from 2010 onwards. The analysis using the virulence factor database (VFDB) indicated that the virulence gene inlJ was the most prevalent among these 44 strains. Notably, ddrA, msbA, and sugC were enriched in this dataset, requiring further clarification of their roles in Listeria through future studies. These results might provide a clue for understanding of the global epidemiology and surveillance of Lm and present insights for implementing effective measures to reduce or prevent Listeria contamination outbreaks in imported food products.

Core-genome multilocus sequence typing and core-SNP analysis of Clostridium neonatale strains isolated in different spatio-temporal settings

IMPORTANCE

Clostridium neonatale has been isolated from the fecal samples of asymptomatic neonates and cases of necrotizing enterocolitis (NEC). Taking advantage of a large collection of independent strains isolated from different spatio-temporal settings, we developed and established a cgMLST scheme for the molecular typing of C. neonatale. Both the cgMLST and cgSNP methods demonstrate comparable discrimination power. Results indicate geographic- and temporal- independent clustering of C. neonatale NEC-associated strains. No specific cgMLST clade of C. neonatale was genetically associated with NEC.

Epidemiological trends and antimicrobial resistance in Salmonella enterica serovar Typhimurium clones in Taiwan between 2004 and 2019

OBJECTIVES

We investigated the temporal trends of Salmonella enterica serovar Typhimurium (S. Typhimurium) clones in Taiwan from 2004 to 2019, focusing on antimicrobial resistance (AMR), resistance genetic determinants, and plasmid types.

METHODS

Salmonella isolates were characterized using pulsed-field gel electrophoresis (PFGE), whole-genome sequencing, and antimicrobial susceptibility testing. Clones were defined using PFGE clustering and the hierarchical cgMLST clustering (HierCC) assignments.

RESULTS

Seven major S. Typhimurium clones, HC100_2, 13, 41, 305, 310, 501, and 46261, accounted for 97.6% (8079/8275) of human isolates in Taiwan. Each clone displayed a unique AMR profile, resistance genetic determinants, and plasmid types. Four highly resistant clones (HC100_2, 41, 305, and 310) exhibited multiple resistance in 86.5% to 96.1% of isolates. HC100_305 and HC100_2 were pandemic multidrug-resistant clones, characterized by resistance to ampicillin, chloramphenicol, streptomycin, sulfonamides, and tetracycline (ACSSuT) and ASSuT, respectively. The prevalence of the ACSSuT clone decreased from 68.7% of S. Typhimurium isolates in 2004 to 1.7% in 2019, while the ASSuT clone emerged in 2007 and became the largest clone after 2010. Several plasmids, including IncHI2-IncHI2A, IncC, IncFIB(K), and IncI1-1(α), carried multiple resistance genes or were associated with the carriage of mph(A), blaCMY-2, and blaDHA-1.

CONCLUSIONS

Between 2004 and 2019, Taiwan experienced the emergence, prevalence, and subsequent decline of several highly resistant S. Typhimurium clones. The clones defined using the HierCC approach have global comparability. The increasing resistance to third-generation cephalosporins, cephamycins, ciprofloxacin, and azithromycin in recent years poses a significant medical concern.

Core genome multilocus sequence typing of Clostridioides difficile to investigate transmission in the hospital setting

PURPOSE

Traditional epidemiological investigations of healthcare-associated Clostridioides difficile infection (HA-CDI) are often insufficient. This study aimed to evaluate a procedure that includes secondary isolation and genomic typing of single toxigenic colonies using core genome multilocus sequence typing (cgMLST) for the investigation of C. difficile transmission.

METHODS

We analyzed retrospectively all toxigenic C. difficile-positive stool samples stored at the Lausanne University Hospital over 6 consecutive months. All isolates were initially typed and classified using a modified double-locus sequence typing (DLST) method. Genome comparison of isolates with the same DLST and clustering were subsequently performed using cgMLST. The electronic administrative records of patients with CDI were investigated for spatiotemporal epidemiological links supporting hospital transmission. A comparative descriptive analysis between genomic and epidemiological data was then performed.

RESULTS

From January to June 2021, 86 C. difficile isolates were recovered from thawed samples of 71 patients. Thirteen different DLST types were shared by > 1 patient, and 13 were observed in single patients. A genomic cluster was defined as a set of isolates from different patients with ≤ 3 locus differences, determined by cgMLST. Seven genomic clusters were identified, among which plausible epidemiological links were identified in only 4/7 clusters.

CONCLUSION

Among clusters determined by cgMLST analysis, roughly 40% included unexplained HA-CDI acquisitions, which may be explained by unidentified epidemiological links, asymptomatic colonization, and/or shared common community reservoirs. The use of DLST, followed by whole genome sequencing analysis, is a promising and cost-effective stepwise approach for the investigation of CDI transmission in the hospital setting.

Multidrug-resistant isolates from Ukrainian patients in a German health facility: a genomic surveillance study focusing on antimicrobial resistance and bacterial relatedness

PURPOSE

Antimicrobial resistance is a pressing issue in Ukraine, with healthcare-associated infections caused by multidrug-resistant organisms being a major concern. A recent prospective multicenter study revealed a staggering rate of 48.4% antimicrobial resistance to carbapenems among Enterobacterales causing a healthcare-associated infection. We conducted a systematic survey to investigate the incidence rate and incidence density of carbapenemase-producing Gram-negative bacteria (CPGN) among refugees and war-wounded Ukrainians in connection with the German health system.

METHODS

From the onset of the war until November 2022, seven Ukrainian patients were admitted to our hospital. Upon admission, screening samples and samples from the focus of suspected infection were taken from all seven patients. The incidence rate and the incidence density of CPGN were calculated as a result of the microbiological findings. We sequenced all CPGN using Illumina technology.

RESULTS

The incidence rate of CPGN at our hospital was 0.06 for 2021 and 0.18 for 2022. All seven Ukrainian patients were infected or colonized with at least one CPGN, including K. pneumoniae (14/25), P. aeruginosa (6/25), A. baumannii (1/25), Providencia stutartii (1/25), C. freundii (1/25), and E. coli (2/25). Genomic surveillance revealed that (i) most frequently detected carbapenemases among all sequenced isolates were blaNDM (17/25) and blaOXA-48 (6/25), (ii) most commonly observed plasmid replicons among the K. pneumoniae isolates recovered from Ukrainian patients were Col(pHAD28) (12/14), IncHI1B(pNDM-MAR) (9/14), IncFIB(pNDM-Mar) (12/14), and (iii) clonal relation between the pathogens of the Ukrainian isolates, but not for the isolates from our hospital surveillance system.

CONCLUSION

The rising prevalence of community-acquired colonization and infection with CPGN is having a direct effect on the infection prevention measures, such as higher number of isolations, reprocessing of patient rooms, additional microbiological testing and overall organization within hospitals.

A core-genome multilocus sequence typing scheme for the detection of genetically related Streptococcus pyogenes clusters

The recently observed increase in invasive Streptococcus pyogenes infections causes concern in Europe. However, conventional molecular typing methods lack discriminatory power to aid investigations of outbreaks caused by S. pyogenes. Therefore, there is an urgent need for high-resolution molecular typing methods to assess genetic relatedness between S. pyogenes isolates. In the current study, we aimed to develop a novel high-resolution core-genome multilocus sequence typing (cgMLST) scheme for S. pyogenes and compared its discriminatory power to conventional molecular typing methods. The cgMLST scheme was designed with the commercial Ridom SeqSphere+ software package. To define a cluster threshold, the scheme was evaluated using publicly available data from nine defined S. pyogenes outbreaks in the United Kingdom. The cgMLST scheme was then applied to 23 isolates from a suspected S. pyogenes outbreak and 117 S. pyogenes surveillance isolates both from the Netherlands. MLST and emm-typing results were used for comparison to cgMLST results. The allelic differences between isolates from defined outbreaks ranged between 6 and 31 for isolates with the same emm-type, resulting in a proposed cluster threshold of <5 allelic differences out of 1,095 target loci. Seven out of twenty-three (30%) isolates from the suspected outbreak had an allelic difference of <2, thereby identifying a potential cluster that could not be linked to other isolates. The proposed cgMLST scheme shows a higher discriminatory ability when compared to conventional typing methods. The rapid and simple analysis workflow allows for extended detection of clusters of potential outbreak isolates and surveillance and may facilitate the sharing of sequencing results between (inter)national laboratories.

Whole-Genome Sequencing-Based Profiling of Antimicrobial Resistance Genes and Core-Genome Multilocus Sequence Typing of Campylobacter jejuni from Different Sources in Lithuania

Campylobacter jejuni is known as one of the main causative agents of gastroenteritis in humans worldwide, and the rise of antimicrobial resistance (AMR) in Campylobacter is a growing public health challenge of special concern. Whole-genome sequencing (WGS) was used to characterize genetic determinants of AMR in 53 C. jejuni isolates from dairy cattle, broiler products, wild birds, and humans in Lithuania. The WGS-based study revealed 26 C. jejuni AMR markers that conferred resistance to various antimicrobials. Genetic markers associated with resistance to beta-lactamases, tetracycline, and aminoglycosides were found in 79.3%, 28.3%, and 9.4% of C. jejuni isolates, respectively. Additionally, genetic markers associated with multidrug resistance (MDR) were found in 90.6% of C. jejuni isolates. The WGS data analysis revealed that a common mutation in the quinolone resistance-determining region (QRDR) was R285K (854G > A) at 86.8%, followed by A312T (934G > A) at 83% and T86I (257C > T) at 71.7%. The phenotypic resistance analysis performed with the agar dilution method revealed that ciprofloxacin (CIP) (90.6%), ceftriaxone (CRO) (67.9%), and tetracycline (TET) (45.3%) were the predominant AMR patterns. MDR was detected in 41.5% (22/53) of the isolates tested. Fifty-seven virulence genes were identified in all C. jejuni isolates; most of these genes were associated with motility (n = 28) and chemotaxis (n = 10). Additionally, all C. jejuni isolates harbored virulence genes related to adhesion, invasion, LOS, LPS, CPS, transportation, and CDT. In total, 16 sequence types (STs) and 11 clonal complexes (CC) were identified based on core-genome MLST (cgMLST) analysis. The data analysis revealed distinct diversity depending on phenotypic and genotypic antimicrobial resistance of C. jejuni.

Listeria monocytogenes ST37 Distribution in the Moscow Region and Properties of Clinical and Foodborne Isolates

Listerias of the phylogenetic lineage II (PLII) are common in the European environment and are hypovirulent. Despite this, they caused more than a third of the sporadic cases of listeriosis and multi-country foodborne outbreaks. L. monocytogenes ST37 is one of them. During the COVID-19 pandemic, ST37 appeared in clinical cases and ranked second in occurrence among food isolates in the Moscow region. The aim of this study was to describe the genomic features of ST37 isolates from different sources. All clinical cases of ST37 were in the cohort of male patients (age, 48-81 years) with meningitis-septicemia manifestation and COVID-19 or Influenza in the anamnesis. The core genomes of the fish isolates were closely related. The clinical and meat isolates revealed a large diversity. Prophages (2-4/genome) were the source of the unique genes. Two clinical isolates displayed pseudolysogeny, and excided prophages were A006-like. In the absence of plasmids, the assortment of virulence factors and resistance determinants in the chromosome corresponded to the hypovirulent characteristics. However, all clinical isolates caused severe disease, with deaths in four cases. Thus, these studies allow us to speculate that a previous viral infection increases human susceptibility to listeriosis.

Whole-genome sequencing-based analysis of antimicrobial resistance, virulence factors, and genetic diversity in Yersinia isolated in Wenzhou, China 2020

Yersinia spp. vary significantly in their ability to cause diseases that threaten public health. Their pathogenicity is frequently associated with increasing antimicrobial resistance (AMR) and various virulence factors. The aim of the study was to investigate the AMR genes, virulence factors, and genetic diversity of Yersinia strains isolated from meats and fish in Wenzhou in 2020 by using whole-genome sequencing (WGS). A total of 50 isolates were collected. The phylogenetic relationships among the Yersinia species were also analyzed using multilocus sequence typing (MLST), core genome multi-locus sequence typing (cgMLST), and single nucleotide polymorphism (SNP) analysis. According to the results, all the strains could be classified into five species, with most isolated from beef, followed by poultry, pork, and fish. AMR genes were identified in 23 strains. And the qnrD1 genes were all located in the Col3M plasmid. Virulence genes, such as yaxA, ystB, pla, and yplA, were also found in the 15 Y. enterocolitica strains. And this study also found the presence of icm/dot type IVB-related genes in one Yersinia massiliensis isolate. MLST analysis identified 43 sequence types (STs), 19 of which were newly detected in Yersinia. Moreover, cgMLST analysis revealed that no dense genotype clusters were formed (cgMLST 5341, 5344, 5346-5350, 5353-5390). Instead, the strains appeared to be dispersed over large distances, except when multiple isolates shared the same ST. Isolates Y4 and Y26 were closely related to strains originating from South Korea and Denmark. This study showed considerable diversity in Yersinia spp. isolated from local areas (Wenzhou City). The data generated in our study may enrich the molecular traceability database of Yersinia and provide a basis for the development of more effective antipathogen control strategies.

Introduction and benchmarking of pyMLST: open-source software for assessing bacterial clonality using core genome MLST

Core genome multilocus sequence typing (cgMLST) has gained in popularity for bacterial typing since whole-genome sequencing (WGS) has become affordable. We introduce here pyMLST, a new complete, stand-alone, free and open source pipeline for cgMLST analysis. pyMLST can create or import a core genome database. For each gene, the first allele is aligned against the bacterial genome of interest using BLAT. Incomplete genes are aligned using MAFT. All data are stored in a SQLite database. pyMLST accepts assembly genomes or raw data (with the option pyMLST-KMA) as input. To evaluate our new tool, we selected three genome collections of major bacterial pathogens (Escherichia coli, Pseudomonas aeruginosa and Staphylococcus aureus) and compared them with pyMLST, pyMLST-KMA, ChewBBACA, SeqSphere and the variant calling approach. We compared the sensitivity, precision and false-positive rate for each method with those of the variant calling approach. Minimal spanning trees were generated with each type of software to evaluate their interest in the context of a bacterial outbreak. We found that pyMLST-KMA is a convenient screening method to avoid assembling large bacterial collections. Our data showed that pyMLST (free, open source, available in Galaxy and pipeline ready) performed similarly to the commercial SeqSphere and performed better than ChewBBACA and pyMLST-KMA.

Influence of transition from open bay units to single room units in a neonatal intensive care unit on hospital transmission of multi-drug-resistant Enterobacterales

BACKGROUND

It was shown previously that changing the design of a hospital neonatal intensive care unit (NICU) from open bay units (OBUs) to single room units (SRUs) was not associated with a reduction in Gram-negative multi-drug-resistant organism (MDRO) colonization rates. It was therefore hypothesized that colonization mainly occurs vertically, or through parents and healthcare workers, and not through environmental factors, and that transition to SRUs would not decrease the number of clusters of MDROs with an epidemiological link. To investigate this, core-genome multi-locus sequence typing (cgMLST) was applied on MDROs cultured from infants at the study hospital.

METHODS

This retrospective cohort study included all infants carrying MDROs admitted to the NICU of a tertiary care academic hospital 2 years prior to the transition from OBUs to SRUs in May 2017, and 1.5 years after the transition (2018-2020).

RESULTS

In total, 55 infants were diagnosed with MDRO carriership. Isolates were available from 49 infants for cgMLST. In the OBU period, one cluster involving four of 20 (20%) infants was identified, and in the SRU period, four clusters involving nine of 29 (31%) infants were identified. It was possible to make an epidemiological link in all four SRU MDRO clusters, but this was not possible for the OBU cluster. In the latter case, transmission from an environmental source on the ward seemed likely.

CONCLUSION

After transition to SRUs, there was no decrease in the number of clusters of MDROs with an epidemiological link, suggesting that nursing infants in an NICU with an SRU design is not, in itself, protective against the acquisition of MDROs.

Combining Whole Genome Sequencing Data from Human and Non-Human Sources: Tackling Listeria monocytogenes Outbreaks

Listeria monocytogenes (Lm) is ubiquitous in nature and known for its ability to contaminate foods during production processes. Near real-time monitoring of whole genome sequences from food and human isolates, complemented with epidemiological data, has been used in the Netherlands since 2019 to increase the speed and success rate of source finding in the case of (active) clusters. Nine clusters with 4 to 19 human cases investigated between January 2019 and May 2023 are described. Fish production sites were most often linked to outbreaks of listeriosis (six clusters), though other types of food businesses can face similar Lm problems, as the production processes and procedures determine risk. The results showed that low levels of Lm in food samples can still be linked to disease. Therefore, the investigation of a cluster of cases and deployment of the precautionary principle helps to focus on safe food and to prevent further cases. Good practice of environmental monitoring within a food business allows early detection of potential issues with food safety and helps food businesses to take appropriate measures such as cleaning to prevent regrowth of Lm and thus future outbreaks.

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.

Analysis of transmission-related third-generation cephalosporin-resistant Enterobacterales by electronic data mining and core genome multi-locus sequence typing

BACKGROUND

To contain intra-hospital transmission of third-generation cephalosporin-resistant Enterobacterales (3GCR-E), contact isolation precautions are recommended.

AIM

To quantify transmissions of 3GCR Escherichia coli and 3GCR Klebsiella pneumoniae within a hospital.

METHODS

An automated outbreak detection system (AODS) was used to identify clusters (N≥2) of 3GCR Enterobacterales for the years 2016, 2018 and 2020. Clusters were defined by phenotypic agreement of microbiological results and spatial and temporal relationship. Core genome multi-locus sequence typing (cgMLST) was used to confirm whether the cluster isolates were transmitted between patients.

FINDINGS

A total of 4343 3GCR E. coli and 1377 K. pneumoniae isolates were analysed. Among the 3GCR E. coli isolates, the AODS identified 304 isolates as cluster isolates, the median cluster size was two (range: 2-5). The cgMLST analysis revealed that a total of 23 (7.5%) 3GCR E. coli cluster isolates were transmission-associated, of which 20 isolates (87%) were detected in intensive care patients. Among the 3GCR K. pneumoniae isolates, the AODS identified 73 isolates as cluster isolates, the median cluster size was two (range: 2-4). CgMLST revealed that 35 (48%) 3GCR K. pneumoniae cluster isolates were transmission associated, of which 27 isolates (77%) were detected in intensive care patients.

CONCLUSION

For 3GCR K. pneumoniae, cgMLST confirmed the AODS results more frequently than for 3GCR E. coli. Therefore, contact isolation precautions for 3GCR K. pneumoniae may be appropriate on intensive care units, but only in certain circumstances, such as outbreaks, for Enterobacterales with lower transmissibility, such as E. coli.

Bifidobacterium longum subsp. iuvenis subsp. nov., a novel subspecies isolated from the faeces of weaning infants

The species Bifidobacterium longum currently comprises four subspecies: B. longum subsp. longum , B. longum subsp. infantis , B. longum subsp. suis and B. longum subsp. suillum . Recently, several studies on B. longum suggested the presence of a separate clade containing four strains isolated from infants and one from rhesus macaque. These strains shared a phylogenetic similarity to B. longum subsp. suis DSM 20210T and B. longum subsp. suillum JCM1995T [average nucleotide identity (ANI) of 98.1 %) while showed an ANI of 96.5 % with both B. longum subsp. infantis and B. longum subsp. longum . The current work describes five novel additional B. longum strains isolated from Bangladeshi weaning infants and demonstrates their common phylogenetic origin with those of the previously proposed separated clade. Based on polyphasic taxonomic approach comprising loci multilocus sequence analysis and whole genome multilocus sequence typing, all ten examined strains have been confirmed as a distinct lineage within the species B. longum with B. longum subsp. suis and B. longum subsp. suillum as closest subspecies. Interestingly, these strains are present in weaning infants and primates as opposed to their closest relatives which have been typically isolated from pig and calves. These strains, similarly to B. longum subsp. infantis , show a common capacity to metabolize the human milk oligosaccharide 3-fucosyllactose. Moreover, they harbour a riboflavin synthesis operon, which differentiate them from their closest subspecies, B. longum subsp. suis and B. longum subsp. suillum . Based on the consistent results from genotypical, ecological and phenotypical analyses, a novel subspecies with the name Bifidobacterium longum subsp. iuvenis, with type strain NCC 5000T (=LMG 32752T=CCOS 2034T), is proposed.

Genomic evidence of genetic diversity and functional evolution in Flavobacterium columnare

Flavobacterium columnare is the causative agent of columnaris disease in freshwater fish. Columnaris disease can cause heavy economic losses in aquaculture. In this study, whole-genome sequencing was used to characterize this pathogen. F. columnare isolate AH-01 had a circular chromosome and plasmid that encoded a total of 3,022 genes. Isolate GX-01 only had a circular chromosome and encoded 2,965 genes. Genomic islands, prophage regions, and CRISPR/Cas systems were identified in both genomes. Both genomes presented evidence of gene variation and horizontal transfer, both of which are the essential components of genetic diversity, genome plasticity, and functional evolution. Single-gene phylogeny and comparative genome analyses were performed to investigate the variation and evolution of this pathogen. Genetic analysis of 16S rRNA and housekeeping gene sequences significantly clustered 55 F. columnare isolates into four clades. The intragroup identity of the 16S rRNA gene exceeded 99%, while the intergroup identity was below the species delineation threshold. We discovered significant translocation, inversion, and rearrangement events that influenced local synteny within each group. Notably, the observed alignments varied considerably among all the studied groups. The core genomes of all strains with available sequences comprised 747 genes, corresponding to approximately 25% of the genome. Core genome multilocus sequence typing, genome-wide orthology and phylogenetic analyses, and average nucleotide identity suggested that the currently existing F. columnare was an assemblage of several distinct species, with levels of divergence at least equivalent to those between recognized bacterial species. The present investigation provided genomic evidence of gene variation and horizontal transfer, which were the basis of genetic diversity, genome plasticity, and functional evolution. The findings supported a proposed new taxonomic perspective on F. columnare.

Evaluation of whole and core genome multilocus sequence typing allele schemes for Salmonella enterica outbreak detection in a national surveillance network, PulseNet USA

Salmonella enterica is a leading cause of bacterial foodborne and zoonotic illnesses in the United States. For this study, we applied four different whole genome sequencing (WGS)-based subtyping methods: high quality single-nucleotide polymorphism (hqSNP) analysis, whole genome multilocus sequence typing using either all loci [wgMLST (all loci)] and only chromosome-associated loci [wgMLST (chrom)], and core genome multilocus sequence typing (cgMLST) to a dataset of isolate sequences from 9 well-characterized Salmonella outbreaks. For each outbreak, we evaluated the genomic and epidemiologic concordance between hqSNP and allele-based methods. We first compared pairwise genomic differences using all four methods. We observed discrepancies in allele difference ranges when using wgMLST (all loci), likely caused by inflated genetic variation due to loci found on plasmids and/or other mobile genetic elements in the accessory genome. Therefore, we excluded wgMLST (all loci) results from any further comparisons in the study. Then, we created linear regression models and phylogenetic tanglegrams using the remaining three methods. K-means analysis using the silhouette method was applied to compare the ability of the three methods to partition outbreak and sporadic isolate sequences. Our results showed that pairwise hqSNP differences had high concordance with cgMLST and wgMLST (chrom) allele differences. The slopes of the regressions for hqSNP vs. allele pairwise differences were 0.58 (cgMLST) and 0.74 [wgMLST (chrom)], and the slope of the regression was 0.77 for cgMLST vs. wgMLST (chrom) pairwise differences. Tanglegrams showed high clustering concordance between methods using two statistical measures, the Baker's gamma index (BGI) and cophenetic correlation coefficient (CCC), where 9/9 (100%) of outbreaks yielded BGI values ≥ 0.60 and CCCs were ≥ 0.97 across all nine outbreaks and all three methods. K-means analysis showed separation of outbreak and sporadic isolate groups with average silhouette widths ≥ 0.87 for outbreak groups and ≥ 0.16 for sporadic groups. This study demonstrates that Salmonella isolates clustered in concordance with epidemiologic data using three WGS-based subtyping methods and supports using cgMLST as the primary method for national surveillance of Salmonella outbreak clusters.

Comparative Analysis of Campylobacter jejuni and C. coli Isolated from Livestock Animals to C. jejuni and C. coli Isolated from Surface Water Using DNA Sequencing and MALDI-TOF

This study evaluated the contribution of cattle, sheep, poultry and pigs to the contamination of surface water from rivers by Campylobacter jejuni and C. coli using MLST, cgMLST and considered MALDI-TOF MS as an alternative technique. The 263 strains isolated from cattle (n = 61), sheep (n = 42), poultry (n = 65), pigs (n = 60) and surface water (n = 35) were distributed across 115 sequence types (STs), 49 for C. jejuni and 66 for C. coli. Considering MLST data, 14.2%, 11.4% and 2.8% of the surface water strains could be attributed to cattle, poultry and sheep, respectively, none to pigs, and 85.7% were non-attributed. Analysis of cg-MLST data with STRUCTURE indicated that C. jejuni strains from water were predominantly attributed to poultry (93.5%), weakly to sheep (<1%) and 6.3% non-attributed, and that conversely, C. coli strains from water were predominantly non-attributed (94.3%) and 5.7% attributed to poultry. Considering the protein profiles with a threshold of 94% and 97% of similarity, respectively, strains from surface water could be attributed to poultry (31.4% and 17.1%), and to cattle (17.1% and 5.7%); 54.1% and 77.1% were non-attributed. This study confirmed these livestock animals might contribute to the contamination of surface water, with a level of contribution depending on the typing technique and the method of analysis. MALDI-TOF could potentially be an alternative approach for source attribution.

Genotype diversity of brucellosis agents isolated from humans and animals in Greece based on whole-genome sequencing

BACKGROUND

Brucellosis is a zoonotic disease whose causative agent, Brucella spp., is endemic in many countries of the Mediterranean basin, including Greece. Although the occurrence of brucellosis must be reported to the authorities, it is believed that the disease is under-reported in Greece, and knowledge about the genomic diversity of brucellae is lacking.

METHODS

Thus, 44 Brucella isolates, primarily B. melitensis, collected between 1999 and 2009 from humans and small ruminants in Greece were subjected to whole genome sequencing using short-read technology. The raw reads and assembled genomes were used for in silico genotyping based on single nucleotide substitutions and alleles. Further, specific genomic regions encoding putative virulence genes were screened for characteristic nucleotide changes, which arose in different genotype lineages.

RESULTS

In silico genotyping revealed that the isolates belonged to three of the known sublineages of the East Mediterranean genotype. In addition, a novel subgenotype was identified that was basal to the other East Mediterranean sublineages, comprising two Greek strains. The majority of the isolates can be assumed to be of endemic origin, as they were clustered with strains from the Western Balkans or Turkey, whereas one strain of human origin could be associated with travel to another endemic region, e.g. Portugal. Further, nucleotide substitutions in the housekeeping gene rpoB and virulence-associated genes were detected, which were characteristic of the different subgenotypes. One of the isolates originating from an aborted bovine foetus was identified as B. abortus vaccine strain RB51.

CONCLUSION

The results demonstrate the existence of several distinct persistent Brucella sp. foci in Greece. To detect these and for tracing infection chains, extensive sampling initiatives are required.

Phenotypic and genotypic characterization of Neisseria gonorrhoeae isolates from Yaoundé, Cameroon, 2019 to 2020

This study investigated antimicrobial resistance (AMR) phenotypes and genotypes exhibited by Neisseria gonorrhoeae from Yaoundé, Cameroon. AMR to tetracycline, penicillin and ciprofloxacin was observed although none of the isolates had reduced susceptibility to azithromycin, cefixime or ceftriaxone. Whole genome sequence (WGS) data were obtained and, using a threshold of 300 or fewer locus differences in the N. gonorrhoeae core gene multilocus sequence typing (cgMLST) scheme, four distinct core genome lineages were identified. Publicly available WGS data from 1355 gonococci belonging to these four lineages were retrieved from the PubMLST database, allowing the Cameroonian isolates to be examined in the context of existing data and compared with related gonococci. Examination of AMR genotypes in this dataset found an association between the core genome and AMR with, for example, isolates belonging to the core genome group, Ng_cgc_300 : 21, possessing GyrA and ParC alleles with amino acid substitutions conferring high-level resistance to ciprofloxacin while lineages Ng_cgc_300 : 41 and Ng_cgc_300 : 243 were predicted to be susceptible to several antimicrobials. A core genome lineage, Ng_cgc_300 : 498, was observed which largely consisted of gonococci originating from Africa. Analyses from this study demonstrate the advantages of using the N. gonorrhoeae cgMLST scheme to find related gonococci to carry out genomic analyses that enhance our understanding of the population biology of this important pathogen.

First Report of Campylobacter jejuni Strains Belonging to ST-21 Clonal Complex Isolated from Human, Poultry and Wild Birds in Croatia: Antimicrobial Resistance and Genetic Distance

In the era of growing antimicrobial resistance, a threat affecting humans, endangering animals, as well as livelihoods and food security worldwide, we wanted to find possible explanations for its continuous spread from a new perspective. The ubiquity of resistance genes requires a One Health approach to finding the explanations for continuous AMR spread. The natural transformability of Campylobacter jejuni, its high incidence of infections, and emerging resistance worldwide inspired us to choose C. jejuni ST-21CC to be our pathogen for analyzing its contribution and connection to the cycle of AMR dissemination. ST-21CC is known as a generalist among humans and broilers, the most prevalent lineage worldwide, but it is rarely found in wild birds. Emerging in wild birds, genetic relatedness and similar resistance profiles were expected. We analyzed 23 Croatian C. jejuni strains belonging specifically to ST-21CC from humans, broilers, and wild birds. The genomic data obtained through whole genome sequencing and phenotypic susceptibility data of strains were compared. Our findings suggest high fluoroquinolone resistance in ST-21CC strains, with more diverse genetic backgrounds in wild birds. Intriguing were three isolates of ST-822 (from human and storks), sharing a similar genetic fingerprint.

Molecular characteristics and phylogenetic analysis of Clostridium perfringens from different regions in China, from 2013 to 2021

Clostridium perfringens (C. perfringens) is a significant foodborne pathogen and a common cause of intestinal diseases in both animals and humans. Our study investigated MLST, phenotypic antimicrobial resistance profiles, and resistance genes among isolates from human, animal and food. 186 C. perfringens isolates were obtained from nine provinces in China between 2013 and 2021. Additionally, some specific ST complexes were analyzed by cgMLST and cgSNP to investigate genetic relatedness. MLST indicated the most prevalent STs of C. perfringens of human and animal origin were as follows: ST221 (5/147), ST62 (4/147), ST408 (4/147), and ST493 (4/147) were predominant in humans, while ST479 (5/25) was the major type in animals. Within the same ST complex, genetically unrelated relationships or potential clustering/transmission events were further recognized by cgMLST and cgSNP, illustrating that these two methods are valuable in defining outbreaks and transmission events. All tested isolates were susceptible to vancomycin and meropenem. The rates of resistance to metronidazole, penicillin, cefoxitin, moxifloxacin, and chloramphenicol were low (metronidazole: 1.08%; penicillin: 9.68%; cefoxitin: 0.54%; moxifloxacin: 6.45%; and chloramphenicol: 3.76%). Interestingly, 49.66% of human origin were clindamycin-resistant, and 18.2% were penicillin-insensitive. Importantly, the portion of MDR isolates was significantly lower than in previous reports. The study provides an overview of the epidemiological characteristics of C. perfringens with different origins and hosts in China. C. perfringens demonstrated remarkable genetic diversity and distinct molecular features compared to antibiotic-resistance profiles from other studies.

Core genome multilocus sequence typing scheme for Bacillus cereus group bacteria

Core genome multilocus sequence typing (cgMLST) employs a strategy where the set of orthologous genes common to all members of a group of organisms are used for phylogenetic analysis of the group members. The Bacillus cereus group consists of species with pathogenicity towards insect species as well as warm-blooded animals including humans. While B. cereus is an opportunistic pathogen linked to a range of human disease conditions, including emesis and diarrhoea, Bacillus thuringiensis is an entomopathogenic species with toxicity toward insect larvae, and therefore used as a biological pesticide worldwide. Bacillus anthracis is a classical obligate pathogen causing anthrax, an acute lethal condition in herbivores as well as humans, and which is endemic in many parts of the world. The group also includes a range of additional species, and B. cereus group bacteria have been subject to analysis with a wide variety of phylogenetic typing systems. Here we present, based on analyses of 173 complete genomes from B. cereus group species available in public databases, the identification of a set of 1568 core genes which were used to create a core genome multilocus typing scheme for the group which is implemented in the PubMLST system as an open online database freely available to the community. The new cgMLST system provides unprecedented resolution over existing phylogenetic analysis schemes covering the B. cereus group.

The rare Salmonella enterica serovar Isangi: genomic characterization of the antimicrobial resistance, virulence potential and epidemiology of Brazilian strains in comparison to global isolates

Introduction. Salmonella enterica serovar Isangi (S. Isangi) is a rare non-typhoidal serovar, related to invasive nosocomial infections in various countries and to increasing antimicrobial resistance rates.Gap statement. Despite existing reports on S. Isangi, there is a lack of information of specific traits regarding this serovar, which could be improved through genomic analyses.Aim. Our goals were to characterize the antimicrobial resistance, virulence potential and genomic relatedness of 11 S. Isangi strains from Brazil in comparison to 185 genomes of global isolates using whole-genome sequencing (WGS) data.Methodology. Phenotypic resistance was determined by disc-diffusion. The search for resistance genes, plasmids, prophages, Salmonella pathogenicity islands (SPIs) and virulence genes, plus multi-locus sequence typing (MLST) and core-genome MLST (cgMLST) were performed using WGS.Results. Brazilian S. Isangi strains showed phenotypic resistance to nalidixic acid, ciprofloxacin and streptomycin, and harboured antimicrobial resistance [qnrB19, aac(6')-Iaa, mdsAB] and heavy metal tolerance (arsD, golST) genes. Col(pHAD28) and IncFII(S) plasmids, virulence genes related to adherence, macrophage induction, magnesium uptake, regulation and type III secretion systems, 12 SPIs and eight prophages were detected. The 185 additional global genomes analysed harboured resistance genes against 11 classes of antimicrobial compounds, 22 types of plasmids, 32 prophages, 14 SPIs, and additional virulence genes related to serum resistance, stress adaptation and toxins. Sequence type (ST)216 was assigned to genomes from Brazil and other countries, while ST335 was the most frequent ST, especially among South African genomes. cgMLST showed that Brazilian genomes were more closely related to genomes from European and African countries, the USA and Taiwan, while the majority of South African genomes were more closely related among each other.Conclusion. The presence of S. Isangi strains from Brazil and different countries showing a close genomic correlation, antimicrobial resistance profiles to drugs used in human therapy and a large number of virulence determinants reinforced the need for stronger initiatives to monitor rare non-typhoidal Salmonella serovars such as S. Isangi in order to prevent its dissemination among human and non-human sources.

Comparison of Brucella abortus population structure based on genotyping methods with different levels of resolution

Numerous genotyping techniques based on different principles and with different costs and levels of resolution are currently available for understanding the transmission dynamics of brucellosis worldwide. We aimed to compare the population structure of the genomes of 53 Brazilian Brucella abortus isolates using eight different genotyping methods: multiple-locus variable-number tandem-repeat analysis (MLVA8, MLVA11, MLVA16), multilocus sequence typing (MLST9, MLST21), core genome MLST (cgMLST) and two techniques based on single nucleotide polymorphism (SNP) detection (parSNP and NASP) from whole genomes. The strains were isolated from six different Brazilian states between 1977 and 2008 and had previously been analyzed using MLVA8, MLVA11, and MLVA16. Their whole genomes were sequenced, assembled, and subjected to MSLT9 MLST21, cgMLST, and SNP analyses. All the genotypes were compared by hierarchical grouping method based on the average distances between the correlation matrices of each technique. MLST9 and MLST21 had the lowest level of resolution, both revealing only four genotypes. MLVA8, MLVA11, and MLVA16 had progressively increasing levels of resolution as more loci were analyzed, identifying 6, 16, and 44 genotypes, respectively. cgMLST showed the highest level of resolution, identifying 45 genotypes, followed by the SNP-based methods, both of which had 44 genotypes. In the assessed population, MLVA was more discriminatory than MLST and was easier and cheaper to perform. SNP techniques and cgMLST provided the highest levels of resolution and the results from the two methods were in close agreement. In conclusion, the choice of genotyping technique can strongly affect one's ability to make meaningful epidemiological conclusions but is dependent on available resources: while the VNTR based techniques are more indicated to high prevalence scenarios, the WGS methods are the ones with the best discriminative power and therefore recommended for outbreaks investigation.

Large Multicountry Outbreak of Invasive Listeriosis by a Listeria monocytogenes ST394 Clone Linked to Smoked Rainbow Trout, 2020 to 2021

Whole-genome sequencing (WGS) has revolutionized surveillance of infectious diseases. Disease outbreaks can now be detected with high precision, and correct attribution of infection sources has been improved. Listeriosis, caused by the bacterium Listeria monocytogenes, is a foodborne disease with a high case fatality rate and a large proportion of outbreak-related cases. Timely recognition of listeriosis outbreaks and precise allocation of food sources are important to prevent further infections and to promote public health. We report the WGS-based identification of a large multinational listeriosis outbreak with 55 cases that affected Germany, Austria, Denmark, and Switzerland during 2020 and 2021. Clinical isolates formed a highly clonal cluster (called Ny9) based on core genome multilocus sequence typing (cgMLST). Routine and ad hoc investigations of food samples identified L. monocytogenes isolates from smoked rainbow trout filets from a Danish producer grouping with the Ny9 cluster. Patient interviews confirmed consumption of rainbow trout as the most likely infection source. The Ny9 cluster was caused by a MLST sequence type (ST) ST394 clone belonging to molecular serogroup IIa, forming a distinct clade within molecular serogroup IIa strains. Analysis of the Ny9 genome revealed clpY, dgcB, and recQ inactivating mutations, but phenotypic characterization of several virulence-associated traits of a representative Ny9 isolate showed that the outbreak strain had the same pathogenic potential as other serogroup IIa strains. Our report demonstrates that international food trade can cause multicountry outbreaks that necessitate cross-border outbreak collaboration. It also corroborates the relevance of ready-to-eat smoked fish products as causes for listeriosis. IMPORTANCE Listeriosis is a severe infectious disease in humans and characterized by an exceptionally high case fatality rate. The disease is transmitted through consumption of food contaminated by the bacterium Listeria monocytogenes. Outbreaks of listeriosis often occur but can be recognized and stopped through implementation of whole-genome sequencing-based pathogen surveillance systems. We here describe the detection and management of a large listeriosis outbreak in Germany and three neighboring countries. This outbreak was caused by rainbow trout filet, which was contaminated by a L. monocytogenes clone belonging to sequence type ST394. This work further expands our knowledge on the genetic diversity and transmission routes of an important foodborne pathogen.

Identification of an Outbreak Cluster of Extensively Antibiotic-Resistant GC1 Acinetobacter baumannii Isolates in U.S. Military Treatment Facilities

An outbreak involving an extensively antibiotic-resistant Acinetobacter baumannii strain in three military treatment facilities was identified. Fifty-nine isolates recovered from 30 patients over a 4-year period were found among a large collection of isolates using core genome multilocus sequence typing (MLST). They differed by only 0 to 18 single nucleotide polymorphisms (SNPs) and carried the same resistance determinants except that the aphA6 gene was missing in 25 isolates. They represent a novel sublineage of GC1 lineage 1 that likely originated in Afghanistan. IMPORTANCE A. baumannii is recognized as one of the most important nosocomial pathogens, and carbapenem-resistant strains pose a particularly difficult treatment challenge. Outbreaks linked to this pathogen are reported worldwide, particularly during periods of societal upheaval, such as natural disasters and conflicts. Understanding how this organism enters and establishes itself within the hospital environment is key to interrupting transmission, but few genomic studies have examined these transmissions over a prolonged period. Though historical, this report provides an in-depth analysis of nosocomial transmission of this organism across continents and within and between different hospitals.

Genomic Diversity of Methicillin-Resistant Staphylococcus aureus CC398 Isolates Collected from Diseased Swine in the German National Resistance Monitoring Program GERM-Vet from 2007 to 2019

Livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) clonal complex 398 (CC398) isolates (n = 178) collected in the national resistance monitoring program GERM-Vet from diseased swine in Germany from 2007 to 2019 were investigated for their genomic diversity with a focus on virulence and antimicrobial resistance (AMR) traits. Whole-genome sequencing was followed by molecular typing and sequence analysis. A minimum spanning tree based on core-genome multilocus sequence typing was constructed, and antimicrobial susceptibility testing was performed. Most isolates were assigned to nine clusters. They displayed close phylogenetic relationships but a wide molecular variety, including 13 spa types and 19 known and four novel dru types. Several toxin-encoding genes, including eta, seb, sek, sep, and seq, were detected. The isolates harbored a wide range of AMR properties mirroring the proportions of the classes of antimicrobial agents applied in veterinary medicine in Germany. Multiple novel or rare AMR genes were identified, including the phenicol-lincosamide-oxazolidinone-pleuromutilin-streptogramin A resistance gene cfr, the lincosamide-pleuromutilin-streptogramin A resistance gene vga(C), and the novel macrolide-lincosamide-streptogramin B resistance gene erm(54). Many AMR genes were part of small transposons or plasmids. Clonal and geographical correlations of molecular characteristics and resistance and virulence genes were more frequently observed than temporal relations. In conclusion, this study provides insight into population dynamics of the main epidemic porcine LA-MRSA lineage in Germany over a 13-year-period. The observed comprehensive AMR and virulence properties, most likely resulting from the exchange of genetic material between bacteria, highlighted the importance of LA-MRSA surveillance to prevent further dissemination among swine husbandry facilities and entry into the human community. IMPORTANCE The LA-MRSA-CC398 lineage is known for its low host specificity and frequent multiresistance to antimicrobial agents. Colonized swine and their related surroundings represent a considerable risk of LA-MRSA-CC398 colonization or infection for occupationally exposed people through which such isolates might be further disseminated within the human community. This study provides insight into the diversity of the porcine LA-MRSA-CC398 lineage in Germany. Clonal and geographical correlations of molecular characteristics and resistance and virulence traits were detected and may be associated with the spread of specific isolates through livestock trade, human occupational exposure, or dust emission. The demonstrated genetic variability underlines the lineage's ability to horizontally acquire foreign genetic material. Thus, LA-MRSA-CC398 isolates have the potential to become even more dangerous for various host species, including humans, due to increased virulence and/or limited therapeutic options for infection control. Full-scale LA-MRSA monitoring at the farm, community, and hospital level is therefore essential.

Whole-genome sequencing for genetic diversity analysis of Iranian Brucella spp. isolated from humans and livestock

Brucellosis is one of the most common zoonoses in the Middle East. It is causing economic losses to the livestock industry and has a great public health concern. Little is known about the genetic diversity and distribution of brucellae in Iran. Therefore, forty Brucella spp. strains (B. abortus and B. melitensis) isolated from animals and humans were analyzed by whole genome sequencing (WGS) technology using single nucleotide polymorphism (SNP) analysis and core genome multilocus sequence typing (cgMLST). Brucella isolates were obtained from lymph nodes (cows and camels), milk (cows, camels and sheep), and aborted foetus samples (sheep and goats), as well as cerebrospinal fluid and blood of humans. The isolates were originating from thirteen provinces of Iran and isolated between 2015 and 2020. According to in-silico MLST, ST8 and ST2 were the most frequent sequence types in B. melitensis and B. abortus, respectively. Based on phylogeographic reconstruction using cgSNP analysis, the investigated Iranian B. melitensis strains belonged to the American and Mediterranean lineages of the B. melitensis phylogeny. Furthermore, cgSNP analysis revealed a similarity between Iranian B. abortus isolates and strains from Iraq and Egypt. Therefore, the origin of the Iranian strains can be suggested to be strains from neighboring and Middle East countries. Moreover, cgMLST analysis showed that the Iranian B. melitensis strains were closely relative to strains recovered from sheep and humans in Iraq, Afghanistan, Syria, Turkmenistan, and Pakistan. In the current panel of strains, cgMLST and cgSNP analysis provided an appropriate and accurate tool for effective traceback analyses for Brucella spp. from Iran. The results of cgSNP and cgMLST helped to understand the geographic distribution and interspecies transmission of Iranian strains and highlight the importance of specific brucellosis control measures in Iran with regard to the One-Health approach.

Characterization of Carbapenemase producing Acinetobacter baumannii Isolates from Danish Patients 2014-2021 - Detection of a New International Clone - IC11

The aim of this study was to characterize carbapenemase producing Acinetobacter baumannii isolates from Danish patients using whole genome sequencing (WGS). Furthermore, we compared typing- and epidemiological data for further investigation of the spread and origin of the carbapenemase producing A. baumannii isolates. From January 1^st 2014 through September 30^th 2021, 141 carbapenemase producing A. baumannii isolates, received at the national reference laboratory at Statens Serum Institut, were investigated using WGS. MLST and cgMLST data, obtained by the SeqSphere+ software, were linked to data related to source of isolation, patient age and gender, hospital admission and travel history. The majority of the carbapenemase producing A. baumannii isolates were from males (n=100, 71%). The majority of patients (n=88, 63%) had travelled outside Scandinavia before admission to a Danish hospital. The most prevalent carbapenemase gene was bla_OXA-23 (n=124). Isolates belonging to the dominating international clone IC2 accounted for 78% of all isolates. A new international ST164/OXA-91 clone, proposedly named IC11, was recognized and described. cgMLST analysis revealed 17 clusters, reflecting both sporadic travel to similar geographical areas as well as confirmed outbreaks in Danish hospitals. The occurrence of carbapenemase producing A. baumannii in Denmark was still low; however, isolates belonging to major international clones with a high potential to spread within hospitals, mainly IC2, dominated. OXA-23 was by far the most prevalent carbapenemase detected. Sporadic and travel related introductions to Danish hospitals, but also intra-hospital transmission could be confirmed, emphasizing the need for continuing vigilance.

Evaluation of core genome and whole genome multilocus sequence typing schemes for Campylobacter jejuni and Campylobacter coli outbreak detection in the USA

Campylobacter is a leading causing of bacterial foodborne and zoonotic illnesses in the USA. Pulsed-field gene electrophoresis (PFGE) and 7-gene multilocus sequence typing (MLST) have been historically used to differentiate sporadic from outbreak Campylobacter isolates. Whole genome sequencing (WGS) has been shown to provide superior resolution and concordance with epidemiological data when compared with PFGE and 7-gene MLST during outbreak investigations. In this study, we evaluated epidemiological concordance for high-quality SNP (hqSNP), core genome (cg)MLST and whole genome (wg)MLST to cluster or differentiate outbreak-associated and sporadic Campylobacter jejuni and Campylobacter coli isolates. Phylogenetic hqSNP, cgMLST and wgMLST analyses were also compared using Baker's gamma index (BGI) and cophenetic correlation coefficients. Pairwise distances comparing all three analysis methods were compared using linear regression models. Our results showed that 68/73 sporadic C. jejuni and C. coli isolates were differentiated from outbreak-associated isolates using all three methods. There was a high correlation between cgMLST and wgMLST analyses of the isolates; the BGI, cophenetic correlation coefficient, linear regression model R ² and Pearson correlation coefficients were >0.90. The correlation was sometimes lower comparing hqSNP analysis to the MLST-based methods; the linear regression model R ² and Pearson correlation coefficients were between 0.60 and 0.86, and the BGI and cophenetic correlation coefficient were between 0.63 and 0.86 for some outbreak isolates. We demonstrated that C. jejuni and C. coli isolates clustered in concordance with epidemiological data using WGS-based analysis methods. Discrepancies between allele and SNP-based approaches may reflect the differences between how genomic variation (SNPs and indels) are captured between the two methods. Since cgMLST examines allele differences in genes that are common in most isolates being compared, it is well suited to surveillance: searching large genomic databases for similar isolates is easily and efficiently done using allelic profiles. On the other hand, use of an hqSNP approach is much more computer intensive and not scalable to large sets of genomes. If further resolution between potential outbreak isolates is needed, wgMLST or hqSNP analysis can be used.

Core genome multilocus sequence typing (cgMLST) confirms systemic spread of avian pathogenic Escherichia coli (APEC) in broilers with cellulitis

Broiler cellulitis has emerged as an important cause of economic losses for farmers and slaughter plants from carcass condemnation at processing. Avian pathogenic Escherichia coli (APEC) has been identified as the main causative agent. The aim was to characterize E. coli isolated from cellulitis and organs in broilers at slaughter by whole genome sequencing analysis to study if systemic spread could be confirmed. Isolates were collected post-mortem from 101 carcasses condemned due to dermatitis/cellulitis from five commercial farms and six flocks. Forty-six isolates were characterised to determine serotypes, sequence types and virulence-associated genes. Analysis by cgMLST was performed to study the genetic similarity between isolates from the same broiler, among birds from the same flock and between flocks. Escherichia coli was isolated from 90% of birds from subcutaneous samples. In 20 broilers, E. coli was isolated from organs in pure culture or mixed with sparse growth of other bacteria. In eight of these, there were post-mortem findings suggestive of systemic bacterial spread. The majority of the isolates from the same bird and flock belonged to the same serotype and sequence type and were genetically indistinguishable, but differed when compared between flocks. Common APEC virulence genes, i.e. chuA, fyuA, hlyF, iroN, irp2, iss, ompT, sitA, TerC, TraT, were present in > 87% of the isolates. We conclude that evidence of systemic spread of E. coli from cellulitis was present in some birds at time of slaughter but cannot be reliably detected at meat inspection.

Core-Genome Multilocus Sequence Typing for Epidemiological and Evolutionary Analyses of Phytopathogenic Xanthomonas citri

Xanthomonas citri subsp. citri is the cause of bacterial citrus canker, responsible for major economic losses to the citrus industry. X. citri subspecies and pathovars are responsible for diseases in soybean, common bean, mango, pomegranate, and cashew. X. citri disease has been tracked using several typing methods, but recent studies using genomic sequencing have been key to understanding the evolutionary relationships within the species, including fundamental differences among X. citri subsp. citri pathotypes. Here, we describe a core-genome multilocus sequence typing (cgMLST) scheme for X. citri based on 250 genomes comprising multiple examples of X. citri subsp. citri pathotypes A, A*, and A^w; X. citri subsp. malvacearum; X. citri pv. aurantifolii, pv. fuscans, pv. glycines, pv. mangiferaeindicae, pv. viticola, and pv. vignicola; and single isolates of X. citri pv. dieffenbachiae and pv. punicae. This data set included genomic sequencing of 100 novel X. citri subsp. citri isolates. cgMLST, based on 1,618 core genes across 250 genomes, is implemented at PubMLST (https://pubmlst.org/organisms/xanthomonas-citri/). GrapeTree minimum-spanning tree and Interactive Tree of Life (iTOL) neighbor-joining phylogenies generated from the cgMLST data resolved almost identical groupings of isolates to a core-genome single nucleotide polymorphism (SNP)-based neighbor-joining phylogeny. These resolved identical groupings of X. citri subsp. citri pathotypes and X. citri subspecies and pathovars. X. citri cgMLST should prove to be an increasingly valuable resource for the study of this key species of plant-pathogenic bacteria. Users can submit genomic data and associated metadata for comparison with previously characterized isolates at PubMLST to allow the rapid characterization of the local, national, and global epidemiology of these pathogens and examine evolutionary relationships. IMPORTANCE Xanthomonas citri is a plant pathogen that causes major economic losses to the citrus industry and sweet orange production in particular. Several subspecies and pathogens are recognized, with host ranges including soybean, common bean, mango, pomegranate, and cashew, among others. Recent genomic studies have shown that host-adapted X. citri subspecies and pathovars and X. citri subsp. citri pathotypes form distinct clades. In this study, we describe a core-genome multilocus sequence typing (cgMLST) scheme for this species that can rapidly and robustly discriminate among these ecologically distinct, host-adapted clades. We have established this scheme and associated databases containing genomic sequences and metadata at PubMLST, which users can interrogate with their own genome sequences to determine X. citri subspecies, pathovars, and pathotypes. X. citri cgMLST should prove to be an invaluable tool for the study of the epidemiology and evolution of this major plant pathogen.

Whole genome sequencing reveals two genetically distinct MRSA outbreaks among people who inject drugs and homeless people in Copenhagen

Whole genome sequencing (WGS) has greatly improved the detection of methicillin-resistant Staphylococcus aureus (MRSA) transmission between people. We describe the transmission of two unique MRSA clones among homeless people in Copenhagen using WGS and core genome MLST (cgMLST). In 2014, an accumulation of MRSA bacteremia cases among homeless people admitted to our hospital was recognized, all having the rare MRSA spa t5147/ST88. The European Typology of Homelessness and Housing Exclusion (ETHOS) categories revealed that people who inject drugs (PWID) frequently visiting the milieu but living in private accommodation accounted for most cases. Hoping to terminate the transmission, 161 homeless people were MRSA screened in 2015, but no additional cases were found. From 2009 to 2018 60 patients with genomically related t5147/ST88 isolates were found, of these 70 % were confirmed to come from the homeless setting and 17 % had bacteremia. In 2017 to 2020 cgMLST revealed a smaller MRSA outbreak including 13 PWID with a completely different clone, t1476/ST8, of which 15 % had bacteremia. Our study confirms that WGS and cgMLST is excellent to reveal MRSA outbreaks. The ETHOS categorisation can be useful to find the primary source of spread in the homeless community.

Genome sequencing and comparative genomic analysis of bovine mastitis-associated Staphylococcus aureus strains from India

BACKGROUND

Bovine mastitis accounts for significant economic losses to the dairy industry worldwide. Staphylococcus aureus is the most common causative agent of bovine mastitis. Investigating the prevalence of virulence factors and antimicrobial resistance would provide insight into the molecular epidemiology of mastitis-associated S. aureus strains. The present study is focused on the whole genome sequencing and comparative genomic analysis of 41 mastitis-associated S. aureus strains isolated from India.

RESULTS

The results elucidate explicit knowledge of 15 diverse sequence types (STs) and five clonal complexes (CCs). The clonal complexes CC8 and CC97 were found to be the predominant genotypes comprising 21 and 10 isolates, respectively. The mean genome size was 2.7 Mbp with a 32.7% average GC content. The pan-genome of the Indian strains of mastitis-associated S. aureus is almost closed. The genome-wide SNP-based phylogenetic analysis differentiated 41 strains into six major clades. Sixteen different spa types were identified, and eight isolates were untypeable. The cgMLST analysis of all S. aureus genome sequences reported from India revealed that S. aureus strain MUF256, isolated from wound fluids of a diabetic patient, was the common ancestor. Further, we observed that all the Indian mastitis-associated S. aureus isolates belonging to the CC97 are mastitis-associated. We identified 17 different antimicrobial resistance (AMR) genes among these isolates, and all the isolates used in this study were susceptible to methicillin. We also identified 108 virulence-associated genes and discuss their associations with different genotypes.

CONCLUSION

This is the first study presenting a comprehensive whole genome analysis of bovine mastitis-associated S. aureus isolates from India. Comparative genomic analysis revealed the genome diversity, major genotypes, antimicrobial resistome, and virulome of clinical and subclinical mastitis-associated S. aureus strains.

Countrywide multi-serotype outbreak of Salmonella Bovismorbificans ST142 and monophasic Salmonella Typhimurium ST34 associated with dried pork sausages in France, September 2020* to January 2021

The French National Reference Centre for Escherichia coli, Shigella and Salmonella (FNRC-ESS) detected two human clusters of 33 cases (median age: 10 years; 17 females) infected by Salmonella enterica serotype Bovismorbificans, ST142, HC5_243255 (EnteroBase HierCC‑cgMLST scheme) in September-November 2020 and of 11 cases (median age: 11 years; seven males) infected by S. enterica serotype 4,12:i:-, ST34, HC5_198125 in October-December 2020. Epidemiological investigations conducted by Santé publique France linked these outbreaks to the consumption of dried pork sausages from the same manufacturer. S. Bovismorbificans and S. 4,12:i:- were isolated by the National Reference Laboratory from different food samples, but both strains were identified in a single food sample only by qPCR. Three recalls and withdrawals of dried pork products were issued by the French general directorate of food of the French ministry for agriculture and food in November 2020, affecting eight supermarket chains. A notification on the European Rapid Alert System for Food and Feed and a European urgent enquiry on the Epidemic Intelligence Information System for Food and Waterborne Diseases and Zoonoses (EPIS-FWD) were launched. No cases were reported outside France. Outbreaks caused by multiple serotypes of Salmonella may go undetected by protocols in standard procedures in microbiology laboratories.

Molecular Features and Antimicrobial Susceptibilities of Streptococcus equi ssp. equi Isolates from Strangles Cases in Indonesia

Strangles, caused by Streptococcus equi ssp. equi (S. equi equi), is a highly infectious and frequent disease of equines worldwide. No data are available regarding the molecular epidemiology of strangles in Indonesia. This study aimed to characterize S. equi equi isolates obtained from suspected strangles cases in Indonesia in 2018. Isolates originated from seven diseased horses on four different farms located in three provinces of Indonesia. Whole genome sequences of these isolates were determined and used for seM typing, multilocus sequence typing (MLST), and core genome MLS typing (cgMLST). Genomes were also screened for known antimicrobial resistance genes and genes encoding for the recombinant antigens used in the commercial Strangvac^® subunit vaccine. All seven S. equi equi isolates from Indonesia belonged to ST179 and carried seM allele 166. Isolates differed from each other by only 2 to 14 cgSNPs and built an exclusive sub-cluster within the Bayesian Analysis of Population Structure (BAPS) cluster 2 (BAPS-2) of the S. equi equi cgMLST scheme. All isolates revealed predicted amino acid sequence identity to seven and high similarity to one of the eight antigen fragments contained in Strangvac^®. Furthermore, all isolates were susceptible to beta-lactam antibiotics penicillin G, ampicillin, and ceftiofur. Our data suggest that the horses from this study were affected by strains of the same novel sublineage within globally distributed BAPS-2 of S. equi equi. Nevertheless, penicillin G can be used as a first-choice antibiotic against these strains and Strangvac^® may also be protective against Indonesian strains.

Carbapenemase-Producing Klebsiella pneumoniae in COVID-19 Intensive Care Patients: Identification of IncL-VIM-1 Plasmid in Previously Non-Predominant Sequence Types

During the COVID-19 pandemic, intensive care units (ICUs) operated at or above capacity, and the number of ICU patients coinfected by nosocomial microorganisms increased. Here, we characterize the population structure and resistance mechanisms of carbapenemase-producing Klebsiella pneumoniae (CP-Kpn) from COVID-19 ICU patients and compare them to pre-pandemic populations of CP-Kpn. We analyzed 84 CP-Kpn isolates obtained during the pandemic and 74 CP-Kpn isolates obtained during the pre-pandemic period (2019) by whole genome sequencing, core genome multilocus sequence typing, plasmid reconstruction, and antibiotic susceptibility tests. More CP-Kpn COVID-19 isolates produced OXA-48 (60/84, 71.4%) and VIM-1 (18/84, 21.4%) than KPC (8/84, 9.5%). Fewer pre-pandemic CP-Kpn isolates produced VIM-1 (7/74, 9.5%). Cefiderocol (97.3-100%) and plazomicin (97.5-100%) had the highest antibiotic activity against pandemic and pre-pandemic isolates. Sequence type 307 (ST307) was the most widely distributed ST in both groups. VIM-1-producing isolates belonging to ST307, ST17, ST321 and ST485, (STs infrequently associated to VIM-1) were detected during the COVID-19 period. Class 1 integron Int1-bla_VIM-1-aac(6')-1b-dfrB1-aadAI-catB2-qacEΔ1/sul1, found on an IncL plasmid of approximately 70,000 bp, carried bla_VIM-1 in ST307, ST17, ST485, and ST321 isolates. Thus, CP-Kpn populations from pandemic and pre-pandemic periods have similarities. However, VIM-1 isolates associated with atypical STs increased during the pandemic, which warrants additional monitoring and surveillance.

New insights into the epidemiology of Listeria monocytogenes - A cross-sectoral retrospective genomic analysis in the Netherlands (2010-2020)

Introduction

Listeriosis, caused by infection with Listeria monocytogenes (Lm), is a relatively rare but severe disease with one of the highest mortality rates among bacterial foodborne illnesses. A better understanding on the degree of Lm clustering, the temporal distribution of the clusters, and their association with the various food sources is expected to lead to improved source tracing and risk-based sampling.

Methods

We investigated the genomic epidemiology of Lm in the Netherlands between 2010 and 2020 by analyzing whole-genome-sequencing (WGS) data of isolates from listerioss patients and food sources from nationwide integrated surveillance and monitoring. WGS data of 756 patient and 770 food/environmental isolates was assessed using core-genome multi-locus sequence typing (cgMLST) with Hamming distance as measure for pairwise distances. Associations of genotype with the epidemiological variables such as patient's age and gender, and systematic use of specific drugs were tested by multinomial logistic regressions. Genetic differentiation of the Lm within and between food categories was calculated based on allele frequencies at the 1701 cgMLST loci in each food category.

Results

We confirmed previous results that some clonal complexes (CCs) are overrepresented among clinical isolates but could not identify any epidemiological risk factors. The main findings of this study include the observation of a very weak attribution of Lm types to food categories and a much better attribution to the producer level. In addition, we identified a high degree of temporal persistence of food, patient and mixed clusters, with more than half of the clusters spanning over more than 1 year and up to 10  years.

Discussion

Taken together this would indicate that identifying persistent contamination in food production settings, and producers that process a wide variety of raw food produce, could significantly contribute to lowering the Lm disease burden.

Split k-mer analysis compared to cgMLST and SNP-based core genome analysis for detecting transmission of vancomycin-resistant enterococci: results from routine outbreak analyses across different hospitals and hospitals networks in Berlin, Germany

The increase of Vancomycin-resistant Enterococcus faecium (VREfm) in recent years has been partially attributed to the rise of specific clonal lineages, which have been identified throughout Germany. To date, there is no gold standard for the interpretation of genomic data for outbreak analyses. New genomic approaches such as split k-mer analysis (SKA) could support cluster attribution for routine outbreak investigation. The aim of this project was to investigate frequent clonal lineages of VREfm identified during suspected outbreaks across different hospitals, and to compare genomic approaches including SKA in routine outbreak investigation. We used routine outbreak laboratory data from seven hospitals and three different hospital networks in Berlin, Germany. Short-read libraries were sequenced on the Illumina MiSeq system. We determined clusters using the published Enterococcus faecium-cgMLST scheme (threshold ≤20 alleles), and assigned sequence and complex types (ST, CT), using the Ridom SeqSphere+ software. For each cluster as determined by cgMLST, we used pairwise core-genome SNP-analysis and SKA at thresholds of ten and seven SNPs, respectively, to further distinguish cgMLST clusters. In order to investigate clinical relevance, we analysed to what extent epidemiological linkage backed the clusters determined with different genomic approaches. Between 2014 and 2021, we sequenced 693 VREfm strains, and 644 (93 %) were associated within cgMLST clusters. More than 74 % (n=475) of the strains belonged to the six largest cgMLST clusters, comprising ST117, ST78 and ST80. All six clusters were detected across several years and hospitals without apparent epidemiological links. Core SNP analysis identified 44 clusters with a median cluster size of three isolates (IQR 2-7, min-max 2-63), as well as 197 singletons (41.4 % of 475 isolates). SKA identified 67 clusters with a median cluster size of two isolates (IQR 2-4, min-max 2-19), and 261 singletons (54.9 % of 475 isolates). Of the isolate pairs attributed to clusters, 7 % (n=3064/45 596) of pairs in clusters determined by standard cgMLST, 15 % (n=1222/8500) of pairs in core SNP-clusters and 51 % (n=942/1880) of pairs in SKA-clusters showed epidemiological linkage. The proportion of epidemiological linkage differed between sequence types. For VREfm, the discriminative ability of the widely used cgMLST based approach at ≤20 alleles difference was insufficient to rule out hospital outbreaks without further analytical methods. Cluster assignment guided by core genome SNP analysis and the reference free SKA was more discriminative and correlated better with obvious epidemiological linkage, at least recently published thresholds (ten and seven SNPs, respectively) and for frequent STs. Besides higher overall discriminative power, the whole-genome approach implemented in SKA is also easier and faster to conduct and requires less computational resources.