Comparison of Whole Genome (wg-) and Core Genome (cg-) MLST (BioNumericsTM) Versus SNP Variant Calling for Epidemiological Investigation of Pseudomonas aeruginosa

Recombinase-aided amplification assay for rapid detection of imipenem-resistant Pseudomonas aeruginosa and rifampin-resistant Pseudomonas aeruginosa

The indiscriminate use of antibiotics has resulted in a growing resistance to drugs in Pseudomonas aeruginosa. The identification of antibiotic resistance genes holds considerable clinical significance for prompt diagnosis. In this study, we established and optimized a Recombinase-Aided Amplification (RAA) assay to detect two genes associated with drug resistance, oprD and arr, in 101 clinically collected P. aeruginosa isolates. Through screening for the detection or absence of oprD and arr, the results showed that there were 52 Imipenem-resistant P. aeruginosa (IRPA) strains and 23 Rifampin-resistant P. aeruginosa (RRPA) strains. This method demonstrated excellent detection performance even when the sample concentration is 10 copies/μL at isothermal conditions and the results could be obtained within 20 minutes. The detection results were in accordance with the results of conventional PCR and Real-time PCR. The detection outcomes of the arr gene were consistently with the resistance spectrum. However, the antimicrobial susceptibility results revealed that 65 strains were resistant to imipenem, while 49 strains sensitive to imipenem with oprD were identified. This discrepancy could be attributed to genetic mutations. In summary, the RAA has higher sensitivity, shorter time, and lower-cost instrument requirements than traditional detection methods. In addition, to analyze the epidemiological characteristics of the aforementioned drug-resistant strains, we conducted Multilocus Sequence Typing (MLST), virulence gene, and antimicrobial susceptibility testing. MLST analysis showed a strong correlation between the sequence types ST-1639, ST-639, ST-184 and IRPA, while ST-261 was the main subtype of RRPA. It was observed that these drug-resistant strains all possess five or more virulence genes, among which exoS and exoU do not coexist, and they are all multidrug-resistant strains. The non-coexistence of exoU and exoS in P.aeruginosa is related to various factors including bacterial regulatory mechanisms and pathogenic mechanisms. This indicates that the relationship between the presence of virulence genes and the severity of patient infection is worthy of attention. In conclusion, we have developed a rapid and efficient RAA (Recombinase-Aided Amplification) detection method that offers significant advantages in terms of speed, simplicity, and cost-effectiveness (especially in time and equipment aspect). This novel approach is designed to meet the demands of clinical diagnostics.

Five-year VIM-producing Pseudomonas aeruginosa outbreak in four Belgian ICUs, an investigation report (2019-2023)

BACKGROUND

Verona integron-encoded metallo-β-lactamase-producing Pseudomonas aeruginosa (VIM-PA) outbreaks are frequently linked to contaminated sink-drains in the intensive care unit (ICU). This study aims to investigate a VIM-PA outbreak occurring at 4 ICUs in a Belgian university center.

METHODS

Between 01/01/2019 and 30/07/2023, data were retrospectively retrieved. Whole-genome sequencing of VIM-PA was carried out for available isolates and the core genome multilocus sequencing typing (cgMLST) was used to confirm clonality. New case incidence was estimated by analyzing the weekly data of at-risk and VIM-PA-colonized patients, fitting a regression model.

RESULTS

Fifty-one patients were colonized, among them, 32 (63%) were infected by VIM-PA, which contributed to 7 deaths. The outbreak investigation showed that 19 (47%) of the examined sink-drains grew at least once a VIM-PA. Two major clusters were observed by cgMLST: ST111 (59 clones with 40 clinical isolates), and ST17 (8 clones with 6 clinical isolates). The estimated incidence rate of new cases was significantly higher in one unit.

CONCLUSIONS

A 5-year prolonged outbreak at the UZ Brussel ICUs was caused by only 2 VIM-PA clones, both linked to sink-drains, with minimal mutations occurring throughout the years. Statistical modeling found different incidence rates between units. Tailored interventions were hence prioritized.

refMLST: reference-based multilocus sequence typing enables universal bacterial typing

BACKGROUND

Commonly used approaches for genomic investigation of bacterial outbreaks, including SNP and gene-by-gene approaches, are limited by the requirement for background genomes and curated allele schemes, respectively. As a result, they only work on a select subset of known organisms, and fail on novel or less studied pathogens. We introduce refMLST, a gene-by-gene approach using the reference genome of a bacterium to form a scalable, reproducible and robust method to perform outbreak investigation.

RESULTS

When applied to multiple outbreak causing bacteria including 1263 Salmonella enterica, 331 Yersinia enterocolitica and 6526 Campylobacter jejuni genomes, refMLST enabled consistent clustering, improved resolution, and faster processing in comparison to commonly used tools like chewieSnake.

CONCLUSIONS

refMLST is a novel multilocus sequence typing approach that is applicable to any bacterial species with a public reference genome, does not require a curated scheme, and automatically accounts for genetic recombination.

AVAILABILITY AND IMPLEMENTATION

refMLST is freely available for academic use at https://bugseq.com/academic .

Tracing the origin of NDM-1-producing and extensively drug-resistant Pseudomonas aeruginosa ST357 in the Netherlands

BACKGROUND

In the hospital environment, carbapenemase-producing Pseudomonas aeruginosa (CPPA) may lead to fatal patient infections. However, the transmission routes of CPPA often remain unknown. Therefore, this case study aimed to trace the origin of CPPA ST357, which caused a hospital-acquired pneumonia in a repatriated critically ill patient suffering from Guillain-Barré Syndrome in 2023.

METHODS

Antimicrobial susceptibility of the CPPA isolate for 30 single and combination therapies was determined by disk-diffusion, Etest or broth microdilution. Whole-genome sequencing was performed for three case CPPA isolates (one patient and two sinks) and four distinct CPPA ST357 patient isolates received in the Dutch CPPA surveillance program. Furthermore, 193 international P. aeruginosa ST357 assemblies were collected via three genome repositories and analyzed using whole-genome multi-locus sequence typing in combination with antimicrobial resistance gene (ARG) characterization.

RESULTS

A Dutch patient who carried NDM-1-producing CPPA was transferred from Kenya to the Netherlands, with subsequent dissemination of CPPA isolates to the local sinks within a month after admission. The CPPA case isolates presented an extensively drug-resistant phenotype, with susceptibility only for colistin and cefiderocol-fosfomycin. Phylogenetic analysis showed considerable variation in allelic distances (mean = 150, max = 527 alleles) among the ST357 isolates from Asia (n = 92), Europe (n = 58), Africa (n = 21), America (n = 16), Oceania (n = 2) and unregistered regions (n = 4). However, the case isolates (n = 3) and additional Dutch patient surveillance program isolates (n = 2) were located in a sub-clade of isolates from Kenya (n = 17; varying 15-49 alleles), the United States (n = 7; 21-115 alleles) and other countries (n = 6; 14-121 alleles). This was consistent with previous hospitalization in Kenya of 2/3 Dutch patients. Additionally, over half of the isolates (20/35) in this sub-clade presented an identical resistome with 9/17 Kenyan, 5/5 Dutch, 4/7 United States and 2/6 other countries, which were characterized by the blaNDM-1, aph(3')-VI, ARR-3 and cmlA1 ARGs.

CONCLUSION

This study presents an extensively-drug resistant subclone of NDM-producing P. aeruginosa ST357 with a unique resistome which was introduced to the Netherlands via repatriation of critically ill patients from Kenya. Therefore, the monitoring of repatriated patients for CPPA in conjunction with vigilance for the risk of environmental contamination is advisable to detect and prevent further dissemination.

Comparison of gene-by-gene and genome-wide short nucleotide sequence-based approaches to define the global population structure of Streptococcus pneumoniae

Defining the population structure of a pathogen is a key part of epidemiology, as genomically related isolates are likely to share key clinical features such as antimicrobial resistance profiles and invasiveness. Multiple different methods are currently used to cluster together closely related genomes, potentially leading to inconsistency between studies. Here, we use a global dataset of 26 306 Streptococcus pneumoniae genomes to compare four clustering methods: gene-by-gene seven-locus MLST, core genome MLST (cgMLST)-based hierarchical clustering (HierCC) assignments, life identification number (LIN) barcoding and k-mer-based PopPUNK clustering (known as GPSCs in this species). We compare the clustering results with phylogenetic and pan-genome analyses to assess their relationship with genome diversity and evolution, as we would expect a good clustering method to form a single monophyletic cluster that has high within-cluster similarity of genomic content. We show that the four methods are generally able to accurately reflect the population structure based on these metrics and that the methods were broadly consistent with each other. We investigated further to study the discrepancies in clusters. The greatest concordance was seen between LIN barcoding and HierCC (adjusted mutual information score=0.950), which was expected given that both methods utilize cgMLST, but have different methods for defining an individual cluster and different core genome schema. However, the existence of differences between the two methods shows that the selection of a core genome schema can introduce inconsistencies between studies. GPSC and HierCC assignments were also highly concordant (AMI=0.946), showing that k-mer-based methods which use the whole genome and do not require the careful selection of a core genome schema are just as effective at representing the population structure. Additionally, where there were differences in clustering between these methods, this could be explained by differences in the accessory genome that were not identified in cgMLST. We conclude that for S. pneumoniae, standardized and stable nomenclature is important as the number of genomes available expands. Furthermore, the research community should transition away from seven-locus MLST, whilst cgMLST, GPSC and LIN assignments should be used more widely. However, to allow for easy comparison between studies and to make previous literature relevant, the reporting of multiple clustering names should be standardized within the research.

Development of the Pneumococcal Genome Library, a core genome multilocus sequence typing scheme, and a taxonomic life identification number barcoding system to investigate and define pneumococcal population structure

Investigating the genomic epidemiology of major bacterial pathogens is integral to understanding transmission, evolution, colonization, disease, antimicrobial resistance and vaccine impact. Furthermore, the recent accumulation of large numbers of whole genome sequences for many bacterial species enhances the development of robust genome-wide typing schemes to define the overall bacterial population structure and lineages within it. Using the previously published data, we developed the Pneumococcal Genome Library (PGL), a curated dataset of 30 976 genomes and contextual data for carriage and disease pneumococci recovered between 1916 and 2018 in 82 countries. We leveraged the size and diversity of the PGL to develop a core genome multilocus sequence typing (cgMLST) scheme comprised of 1222 loci. Finally, using multilevel single-linkage clustering, we stratified pneumococci into hierarchical clusters based on allelic similarity thresholds and defined these with a taxonomic life identification number (LIN) barcoding system. The PGL, cgMLST scheme and LIN barcodes represent a high-quality genomic resource and fine-scale clustering approaches for the analysis of pneumococcal populations, which support the genomic epidemiology and surveillance of this leading global pathogen.

Development and implementation of a core genome multilocus sequence typing scheme for Haemophilus influenzae

Haemophilus influenzae is part of the human nasopharyngeal microbiota and a pathogen causing invasive disease. The extensive genetic diversity observed in H. influenzae necessitates discriminatory analytical approaches to evaluate its population structure. This study developed a core genome multilocus sequence typing (cgMLST) scheme for H. influenzae using pangenome analysis tools and validated the cgMLST scheme using datasets consisting of complete reference genomes (N = 14) and high-quality draft H. influenzae genomes (N = 2297). The draft genome dataset was divided into a development dataset (N = 921) and a validation dataset (N = 1376). The development dataset was used to identify potential core genes, and the validation dataset was used to refine the final core gene list to ensure the reliability of the proposed cgMLST scheme. Functional classifications were made for all the resulting core genes. Phylogenetic analyses were performed using both allelic profiles and nucleotide sequence alignments of the core genome to test congruence, as assessed by Spearman's correlation and ordinary least square linear regression tests. Preliminary analyses using the development dataset identified 1067 core genes, which were refined to 1037 with the validation dataset. More than 70% of core genes were predicted to encode proteins essential for metabolism or genetic information processing. Phylogenetic and statistical analyses indicated that the core genome allelic profile accurately represented phylogenetic relatedness among the isolates (R 2 = 0.945). We used this cgMLST scheme to define a high-resolution population structure for H. influenzae, which enhances the genomic analysis of this clinically relevant human pathogen.

Epidemiology of healthcare-associated Pseudomonas aeruginosa in intensive care units: are sink drains to blame?

BACKGROUND

Pseudomonas aeruginosa (PA) is a common cause of healthcare-associated infection (PA-HAI) in the intensive care unit (ICU).

AIM

To describe the epidemiology of PA-HAI in ICUs in Ontario, Canada, and to identify episodes of sink-to-patient PA transmission.

METHODS

This was a prospective cohort study of patients in six ICUs from 2018 to 2019, with retrieval of PA clinical isolates, and PA-screening of antimicrobial-resistant organism surveillance rectal swabs, and of sink drain, air, and faucet samples. All PA isolates underwent whole-genome sequencing. PA-HAI was defined using US National Healthcare Safety Network criteria. ICU-acquired PA was defined as PA isolated from specimens obtained ≥48 h after ICU admission in those with prior negative rectal swabs. Sink-to-patient PA transmission was defined as ICU-acquired PA with close genomic relationship to isolate(s) previously recovered from sinks in a room/bedspace occupied 3-14 days prior to collection date of the relevant patient specimen.

FINDINGS

Over ten months, 72 PA-HAIs occurred among 60/4263 admissions. The rate of PA-HAI was 2.40 per 1000 patient-ICU-days; higher in patients who were PA-colonized on admission. PA-HAI was associated with longer stay (median: 26 vs 3 days uninfected; P < 0.001) and contributed to death in 22/60 cases (36.7%). Fifty-eight admissions with ICU-acquired PA were identified, contributing 35/72 (48.6%) PA-HAIs. Four patients with five PA-HAIs (6.9%) had closely related isolates previously recovered from their room/bedspace sinks.

CONCLUSION

Nearly half of PA causing HAI appeared to be acquired in ICUs, and 7% of PA-HAIs were associated with sink-to-patient transmission. Sinks may be an under-recognized reservoir for HAIs.

Genomic insights of Salmonella isolated from dry fermented sausage production chains in Spain and France

The presence of Salmonella in dry fermented sausages is source of recalls and outbreaks. The genomic diversity of 173 Salmonella isolates from the dry fermented sausage production chains (pig carcasses, pork, and sausages) from France and Spain were investigated through their core phylogenomic relationships and accessory genome profiles. Ten different serovars and thirteen sequence type profiles were identified. The most frequent serovar from sausages was the monophasic variant of S. Typhimurium (1,4,[5],12:i:-, 72%) while S. Derby was in pig carcasses (51%). Phylogenomic clusters found in S. 1,4,[5],12:i:-, S. Derby, S. Rissen and S. Typhimurium serovars identified closely related isolates, with less than 10 alleles and 20 SNPs of difference, displaying Salmonella persistence along the pork production chain. Most of the S. 1,4,[5],12:i:- contained the Salmonella genomic island-4 (SGI-4), Tn21 and IncFIB plasmid. More than half of S. Derby strains contained the SGI-1 and Tn7. S. 1,4,[5],12:i:- genomes carried the most multidrug resistance genes (91% of the strains), whereas extended-spectrum β-lactamase genes were found in Typhimurium and Derby serovars. Salmonella monitoring and characterization in the pork production chains, specially S. 1,4,[5],12:i:- serovar, is of special importance due to its multidrug resistance capacity and persistence in dry fermented sausages.

Phenotypic and molecular characterization of multi-drug resistant Klebsiella spp. isolates recovered from clinical settings

Klebsiella pneumoniae is a Gram-negative bacterium that colonizes the gastrointestinal tract and nasopharynx with many being linked to nosocomial infections. Extended-spectrum β-lactamases (ESBL)-producing and carbapenem-resistant K. pneumoniae is recognized by the World Health Organization (WHO) as a critical public health concern. In this study, whole-genome sequencing (WGS) - based analysis was performed to understand the molecular epidemiology of multi-drug resistant Klebsiella spp. clinical isolates. Genome comparison, multi-locus sequence typing (MLST), pulsed-field gel electrophoresis (PFGE), and whole-genome-SNP-based phylogenetic analysis (wg-SNP) were used for in-depth molecular characterization. in silico typing was used to determine the resistance genes, virulence factors, Inc. groups, and capsular types. All except one isolate were non-susceptible to meropenem and 89% were non-susceptible to ertapenem and imipenem. blaNDM, blaOXA-48, and blaKPC were the detected carbapenemases with blaNDM-1 found in half of the sequenced genomes. Resistance to colistin was detected in one isolate and was linked to several genetic alterations in crrB, pmrB, and pmrC genes. The most common plasmid type was IncFIB followed by IncR, and the Type 3 fimbriae, encoded by the mrkABCDF operon, was conserved among all isolates. The most common sequence- (ST) and K-type detected were ST147 and K64. The prevelance and the genomic relatedness of ST147 isolates, as shown by the data from SNP-based phylogenetic analysis, PFGE, and genomic clustering, may be an outbreak marker. However, this can only be validated through a more comprehensive study encompassing a wider sampling scheme and over an extended timeframe.

gcPathogen: a comprehensive genomic resource of human pathogens for public health

Here, we present the manually curated Global Catalogue of Pathogens (gcPathogen), an extensive genomic resource designed to facilitate rapid and accurate pathogen analysis, epidemiological exploration and monitoring of antibiotic resistance features and virulence factors. The catalogue seamlessly integrates and analyzes genomic data and associated metadata for human pathogens isolated from infected patients, animal hosts, food and the environment. The pathogen list is supported by evidence from medical or government pathogenic lists and publications. The current version of gcPathogen boasts an impressive collection of 1 164 974 assemblies comprising 986 044 strains from 497 bacterial taxa, 4794 assemblies encompassing 4319 strains from 265 fungal taxa, 89 965 assemblies featuring 13 687 strains from 222 viral taxa, and 646 assemblies including 387 strains from 159 parasitic taxa. Through this database, researchers gain access to a comprehensive 'one-stop shop' that facilitates global, long-term public health surveillance while enabling in-depth analysis of genomes, sequence types, antibiotic resistance genes, virulence factors and mobile genetic elements across different countries, diseases and hosts. To access and explore the data and statistics, an interactive web interface has been developed, which can be accessed at https://nmdc.cn/gcpathogen/. This user-friendly platform allows seamless querying and exploration of the extensive information housed within the gcPathogen database.

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.

New insights into the pathogenesis and transmission of Brucella pinnipedialis: systemic infection in two bottlenose dolphins (Tursiops truncatus)

IMPORTANCE

Brucella spp. are zoonotic pathogens that can affect both terrestrial and marine mammals. Brucella ceti has been identified in various cetacean species, but only one sequence type (ST27) has been reported in humans. However, it is important to conduct surveillance studies to better understand the impact of marine Brucella species on marine mammals, a typically understudied host group. Here, we describe a systemic infection by two related strains of Brucella pinnipedialis (ST25) in a couple of live-stranded bottlenose dolphins, with more severe lesions in the younger animal. Furthermore, B. pinnipedialis was first detected in milk from a female cetacean that stranded with its offspring. Our study reveals novel insights into the epidemiology and pathological consequences of B. pinnipedialis infections in cetaceans, emphasizing the crucial importance of ongoing surveillance and accurate diagnosis to understand the impact of this pathogen on marine mammal populations.

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.

Assessment of plasmids for relating the 2020 Salmonella enterica serovar Newport onion outbreak to farms implicated by the outbreak investigation

BACKGROUND

The Salmonella enterica serovar Newport red onion outbreak of 2020 was the largest foodborne outbreak of Salmonella in over a decade. The epidemiological investigation suggested two farms as the likely source of contamination. However, single nucleotide polymorphism (SNP) analysis of the whole genome sequencing data showed that none of the Salmonella isolates collected from the farm regions were linked to the clinical isolates-preventing the use of phylogenetics in source identification. Here, we explored an alternative method for analyzing the whole genome sequencing data driven by the hypothesis that if the outbreak strain had come from the farm regions, then the clinical isolates would disproportionately contain plasmids found in isolates from the farm regions due to horizontal transfer.

RESULTS

SNP analysis confirmed that the clinical isolates formed a single, nearly-clonal clade with evidence for ancestry in California going back a decade. The clinical clade had a large core genome (4,399 genes) and a large and sparsely distributed accessory genome (2,577 genes, at least 64% on plasmids). At least 20 plasmid types occurred in the clinical clade, more than were found in the literature for Salmonella Newport. A small number of plasmids, 14 from 13 clinical isolates and 17 from 8 farm isolates, were found to be highly similar (> 95% identical)-indicating they might be related by horizontal transfer. Phylogenetic analysis was unable to determine the geographic origin, isolation source, or time of transfer of the plasmids, likely due to their promiscuous and transient nature. However, our resampling analysis suggested that observing a similar number and combination of highly similar plasmids in random samples of environmental Salmonella enterica within the NCBI Pathogen Detection database was unlikely, supporting a connection between the outbreak strain and the farms implicated by the epidemiological investigation.

CONCLUSION

Horizontally transferred plasmids provided evidence for a connection between clinical isolates and the farms implicated as the source of the outbreak. Our case study suggests that such analyses might add a new dimension to source tracking investigations, but highlights the need for detailed and accurate metadata, more extensive environmental sampling, and a better understanding of plasmid molecular evolution.

Genomic Epidemiology and Phenotypic Characterization of Staphylococcus aureus from a Tertiary Hospital in Tianjin Municipality, Northern China

Staphylococcus aureus remains a dangerous pathogen and poses a great threat to public health worldwide. The prevalence of the S. aureus clonotype is temporally and geographically variable. The genomic and phenotypic characteristics of S. aureus isolates in Tianjin, which is among the four big municipalities in China, are unclear. In the present study, 201 nonduplicate S. aureus isolates, including 70 methicillin-resistant S. aureus (MRSA) and 131 methicillin-susceptible S. aureus (MSSA), were collected from 2015 to 2021 in a tertiary hospital in Tianjin. Whole-genome sequencing of S. aureus isolates was carried out to investigate bacterial molecular characteristics, genomic phylogeny, antimicrobial resistance (AMR) gene carriage, and virulence factor gene distribution. The antibiotic resistance profiles, hemolytic activities, and biofilm formation abilities of the S. aureus isolates were also determined. In total, 31 distinct sequence types (STs) and 68 spa types were identified. ST59 (15.9%, 32/201) was the predominant clonotype, followed by ST398 (14.9%, 30/201) and several other major STs (ST1, ST5, ST6, ST22, ST25, ST188, and the newly emerging ST5527). ST59 and ST5527 mainly included MRSA isolates, while ST398 and the other major STs mainly included MSSA isolates. The unique characteristics of the S. aureus isolates belonging to the major STs were determined. ST59 isolates exhibited strong hemolytic activity, and ST398 strains had high biofilm formation capacity, while ST5527 isolates presented the greatest AMR. The genomic epidemiology and phenotypic characteristics of S. aureus isolates determined in this study will help in disease control in nosocomial environments. IMPORTANCE Staphylococcus aureus is an important bacterium pathogen in tertiary hospitals, which provide rich medical resources. Tianjin is one of the four municipalities in China with a population of more than 13 million. However, the epidemiology and molecular characteristics of S. aureus isolates in Tianjin are unknown. In this study, the genomic and phenotypic analyses were performed to investigate 201 S. aureus isolates collected from a tertiary hospital in Tianjin over a time span of 6 years. The refined analysis of predominant clones ST59, ST398, the newly emerging clone ST5527, as well as other major clones, will undoubtedly aid in the control and prevention of infections caused by S. aureus in tertiary hospitals.

Whole-genome sequencing reveals genomic characterization of Listeria monocytogenes from food in China

Introduction

Listeria monocytogenes is a foodborne bacterium that could persist in food and food processing environments for a long time. Understanding the population structure and genomic characterization of foodborne L. monocytogenes is essential for the prevention and control of listeriosis.

Methods

A total of 322 foodborne L. monocytogenes isolates from 13 geographical locations and four food sources in China between 2000 and 2018 were selected for whole-genome sequencing.

Results

In silico subtyping divided the 322 isolates into five serogroups, 35 sequence types (STs), 26 clonal complexes (CCs) and four lineages. Serogroup IIa was the most prevalent serogroup and ST9 was the most prevalent ST of foodborne L. monocytogenes strains isolated in China. The in-depth phylogenetic analysis on CC9 revealed that ST122 clone might be original from ST9 clone. Furthermore, 23 potentially relevant clusters were identified by pair-wised whole-genome single nucleotide polymorphism analysis, indicating that persistent- and/or cross-contamination had occurred in markets in China. ST8 and ST121 were the second and third top STs of L. monocytogenes in China, which had heterogeneity with that of L. monocytogenes isolates from other countries. The antibiotic resistance genes aacA4, tetM, tetS, dfrG carried by different mobile elements were found in L. monocytogenes strains. One lineage II strain carrying Listeria Pathogenicity Island 3 was first reported. In addition, a novel type of premature stop codon in inlA gene was identified in this study.

Discussion

These findings revealed the genomic characteristics and evolutionary relationship of foodborne L. monocytogenes in China on a scale larger than previous studies, which further confirmed that whole-genome sequencing analysis would be a helpful tool for routine surveillance and source-tracing investigation.

Pseudomonas aeruginosa left ventricular assist device (LVAD) driveline infection acquired from the bathroom at home

We describe a patient with a left ventricular assist device (LVAD) infection by Pseudomonas aeruginosa acquired at home. The Pseudomonas from the driveline was similar to several surface cultures of the patient's home shower. This case illustrates the potential and importance of infection prevention measures at home.

Development and Evaluation of a Core Genome Multilocus Sequencing Typing (cgMLST) Scheme for Serratia marcescens Molecular Surveillance and Outbreak Investigations

Serratia marcescens can cause a range of severe infections and contributes to nosocomial outbreaks. Although whole-genome sequencing (WGS)-based typing is the standard method for molecular surveillance and outbreak investigation, there is no standardized analytic scheme for S. marcescens core genome multilocus sequence typing (cgMLST). Here, the development and evaluation of a S. marcescens cgMLST scheme is reported with the goal of enabling a standardized methodology and typing nomenclature. Four hundred ninety-one high-quality S. marcescens WGS data sets were extracted from public databases and-using the genomic sequence of NCBI reference strain S. marcescens Db11 (NZ_HG326223.1) as a starting point-all Db11 genes present in ≥97% data sets used to create a cgMLST scheme. The novel scheme was evaluated using WGS data from 24 outbreak investigations (n = 175 isolates) distributed over three continents. Analysis of Db11 genes within the 491 data sets identified 2,692 target genes present in ≥97% of genomes (mean, 99.1%; median, 99.9%). These genes formed the novel cgMLST scheme, covering 47.8% of nucleotides in the Db11 genome. Analyzing 175 isolates from 24 outbreaks using the novel scheme gave comparable results to previous typing efforts for both general groupings and allelic distances within clusters. In summary, a novel cgMLST scheme for S. marcescens was developed and evaluated. The scheme and its associated nomenclature will improve standardization of typing efforts for molecular surveillance and outbreak investigation, allowing better understanding of S. marcescens genomic epidemiology and facilitating interlaboratory comparisons.

Whole Genome Multi-Locus Sequence Typing and Genomic Single Nucleotide Polymorphism Analysis for Epidemiological Typing of Pseudomonas aeruginosa From Indonesian Intensive Care Units

We have previously studied carbapenem non-susceptible Pseudomonas aeruginosa (CNPA) strains from intensive care units (ICUs) in a referral hospital in Jakarta, Indonesia (Pelegrin et al., 2019). We documented that CNPA transmissions and acquisitions among patients were variable over time and that these were not significantly reduced by a set of infection control measures. Three high risk international CNPA clones (sequence type (ST)235, ST823, ST357) dominated, and carbapenem resistance was due to carbapenemase-encoding genes and mutations in the porin OprD. Pelegrin et al. (2019) reported core genome analysis of these strains. We present a more refined and detailed whole genome-based analysis of major clones represented in the same dataset. As per our knowledge, this is the first study reporting Single Nucleotide Polymorphisms (wgSNP) analysis of Pseudomonas strains. With whole genome-based Multi Locus Sequence Typing (wgMLST) of the 3 CNPA clones (ST235, ST357 and ST823), three to eleven subgroups with up to 200 allelic variants were observed for each of the CNPA clones. Furthermore, we analyzed these CNPA clone clusters for the presence of wgSNP to redefine CNPA transmission events during hospitalization. A maximum number 35350 SNPs (including non-informative wgSNPs) and 398 SNPs (ST-specific_informative-wgSNPs) were found in ST235, 34,570 SNPs (including non-informative wgSNPs) and 111 SNPs (ST-specific_informative-wgSNPs) in ST357 and 26,443 SNPs (including non-informative SNPs) and 61 SNPs (ST-specific_informative-wgSNPs) in ST823. ST-specific_Informative-wgSNPs were commonly noticed in sensor-response regulator genes. However, the majority of non-informative wgSNPs was found in conserved hypothetical proteins or in uncharacterized proteins. Of note, antibiotic resistance and virulence genes segregated according to the wgSNP analyses. A total of 8 transmission chains for ST235 strains followed by 9 and 4 possible transmission chains for ST357 and ST823 were traceable on the basis of pairwise distances of informative-wgSNPs (0 to 4 SNPs) among the strains. The present study demonstrates the value of detailed whole genome sequence analysis for highly refined epidemiological analysis of P. aeruginosa.

Performance Evaluation of the IR Biotyper® System for Clinical Microbiology: Application for Detection of Staphylococcus aureus Sequence Type 8 Strains

Background: Methicillin-resistant S. aureus (MRSA) clonal lineages have been classified based on sequence type (ST) and pulsotype associated with human infection. Providing rapid and accurate epidemiological insight is important to address proper infection control in both community-acquired and nosocomial hospital settings. In this regard, this study was performed to evaluate the IR Biotyper^® (IRBT^®) for strain typing of S. aureus clinical isolates on three media. Methods: A total of 24 S. aureus clinical isolates comprising 15 MRSA isolates (six ST5, three ST72, three ST8, and three ST188 isolates) and nine methicillin-susceptible S. aureus (MSSA) isolates (three ST5, three ST72, and three ST8 isolates) were included for evaluating the IRBT^®. Molecular characterization of all S. aureus isolates was performed by conventional PCR and sequencing methods. The IRBT^® was evaluated according to manufacturer instructions and a modified sample procedure on commonly used BAP, MHA, and TSA media. Subsequently, the spectra obtained by IRBT^® software were compared with dendrograms of PFGE analysis. Results: In this study, the modified sample procedure for reducing the amount of bacteria and bacterial concentration improved the acquisition quality pass rate of the IRBT^®. Each spectrum of S. aureus ST5, ST72, and ST188 isolates on all three media could not be clustered by IRBT^®. However, the dendrogram obtained from the spectra of S. aureus ST8 isolates on TSA medium were in concordance with that obtained by PFGE analysis. In addition, the visual distribution of S. aureus ST8 isolates on TSA medium in a 2D scatter plot appeared as separated point set from those of S. aureus ST5, ST72, and ST188 isolates. Conclusions: The IRBT^® system is a rapid strain typing tool using the FTIR spectroscopic method. This system demonstrated the possibility of discriminating the strain types of S. aureus clinical isolates. Indeed, S. aureus ST8 isolates on TSA medium were successfully differentiated from other strain type isolates.

Comparative genome analysis of Pasteurella multocida serogroup B:2 strains causing haemorrhagic septicaemia (HS) in bovines

Pasteurella multocida, a Gram-negative bacterium with ubiquitous nature, is known to affect wide range of host species worldwide with varied clinical manifestations including haemorrhagic septicaemia (HS) in bovines. Although, HS causing P. multocida strains were identified and characterized by conventional tools and PCR assays, diverse strains are indistinguishable by these tools in the face of disease outbreaks. In this study, draft genomes of three virulent P. multocida serotype B:2 strains (NIVEDIPm32, NIVEDIPm34 and NIVEDIPm35) were analyzed following whole genome sequencing, assembly, annotation and compared them with existing global genomes (n = 43) of bovine origin in the database. Three draft genomes of NIVEDIPm strains consisted of 40-52 contigs with GC content of ∼40.4%. The genome size and predicted genes content was ∼2.3 Mb and 2181-2189, respectively. Besides, the presence of various mobile genetic elements, antimicrobial resistance genes and biofilm related genes suggested their vital roles in virulence; further, adaptation to the host immune system as well as host pathogen interaction. Multi locus sequence analysis based on RIRDC scheme showed the presence of ST122 in all the three strains. wgMLST based phylogenic analysis suggested that HS causing Indian virulent field strains differed geographically and showed diversity from existing HS vaccine strain P52. The phylogenetic tree revealed that North Indian strains share high similarity with strains of Pakistan than South Indian Strain. Notably, a high divergence of SNPs between the HS causing circulating virulent strains of India and current HS vaccine strain P52 suggested an imminent need for relook in to HS vaccination strategy for livestock in India.

Majorbio Cloud: A one-stop, comprehensive bioinformatic platform for multiomics analyses

The platform consists of three modules, which are pre-configured bioinformatic pipelines, cloud toolsets, and online omics' courses. The pre-configured bioinformatic pipelines not only combine analytic tools for metagenomics, genomes, transcriptome, proteomics and metabolomics, but also provide users with powerful and convenient interactive analysis reports, which allow them to analyze and mine data independently. As a useful supplement to the bioinformatics pipelines, a wide range of cloud toolsets can further meet the needs of users for daily biological data processing, statistics, and visualization. The rich online courses of multi-omics also provide a state-of-art platform to researchers in interactive communication and knowledge sharing.

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

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

In vitro and in silico parameters for precise cgMLST typing of Listeria monocytogenes

Background

Whole genome sequencing analyzed by core genome multi-locus sequence typing (cgMLST) is widely used in surveillance of the pathogenic bacteria Listeria monocytogenes. Given the heterogeneity of available bioinformatics tools to define cgMLST alleles, our aim was to identify parameters influencing the precision of cgMLST profiles.

Methods

We used three L. monocytogenes reference genomes from different phylogenetic lineages and assessed the impact of in vitro (i.e. tested genomes, successive platings, replicates of DNA extraction and sequencing) and in silico parameters (i.e. targeted depth of coverage, depth of coverage, breadth of coverage, assembly metrics, cgMLST workflows, cgMLST completeness) on cgMLST precision made of 1748 core loci. Six cgMLST workflows were tested, comprising assembly-based (BIGSdb, INNUENDO, GENPAT, SeqSphere and BioNumerics) and assembly-free (i.e. kmer-based MentaLiST) allele callers. Principal component analyses and generalized linear models were used to identify the most impactful parameters on cgMLST precision.

Results

The isolate’s genetic background, cgMLST workflows, cgMLST completeness, as well as depth and breadth of coverage were the parameters that impacted most on cgMLST precision (i.e. identical alleles against reference circular genomes). All workflows performed well at ≥40X of depth of coverage, with high loci detection (> 99.54% for all, except for BioNumerics with 97.78%) and showed consistent cluster definitions using the reference cut-off of ≤7 allele differences.

Conclusions

This highlights that bioinformatics workflows dedicated to cgMLST allele calling are largely robust when paired-end reads are of high quality and when the sequencing depth is ≥40X.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08437-4.

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

Background

Multidrug-resistant Acinetobacter baumannii is an important causal pathogen of healthcare-associated infections, and colistin-resistant strains have recently emerged owing to the increased use of colistin. Using next-generation sequencing (NGS), a single whole-genome sequencing (WGS) protocol can identify and type pathogens, analyze genetic relationships among different pathogens, predict pathogenic transmissions, and detect antibiotic resistance genes. However, only a few studies have applied NGS in studying the resistance mechanism and epidemiology of colistin-resistant A. baumannii. This study aimed to elucidate the resistance mechanism of colistin-resistant A. baumannii and analyze its molecular epidemiology through WGS.

Materials and methods

The subjects in this study were patients who visited a university hospital between 2014 and 2018. Thirty colistin-resistant strains with high minimum inhibitory concentrations were selected from various patient samples, and WGS was performed. Comparative genomic analysis was performed for the 27 colistin-resistant A. baumannii strains using a colistin-susceptible strain as the reference genome.

Results

The WGS analysis found no mutation for lpxA, lpxC, lpx D, pmrA, pmrB, and mcr1, the genes known to be associated with colistin resistance. Fifty-seven coding sequences (CDS) showed differences; they included 13 CDS with known names and functions that contained 21 genes. From the whole-genome multi-locus sequence typing (wgMLST) and single nucleotide polymorphism (SNP) analyses, two major clusters were found for the colistin-resistant A. baumannii strains. However, no differences were observed by the time of detection for each cluster, the samples, the pattern of antibiotic resistance, or the patient characteristics. In the conventional MLST following the Oxford scheme, the typing result showed ST1809, ST451, ST191, ST1837, and ST369 in the global clone 2 (GC2), without any relation with the results of wgMLST and SNP analyses.

Conclusion

Based on the findings of the resistance gene analysis through WGS and comparative genomic analysis, the potential genes associated with colistin-resistance or CDS were examined. Furthermore, the analysis of molecular epidemiology through WGS regarding colistin-resistant A. baumannii may prove helpful in preventing infection by multidrug-resistant bacteria and controlling healthcare-associated infections.

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

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

Staphylococcus aureus and Cystic Fibrosis-A Close Relationship. What Can We Learn from Sequencing Studies?

Staphylococcus aureus is next to Pseudomonas aeruginosa the most isolated pathogen from the airways of cystic fibrosis (CF) patients, who are often infected by a dominant S. aureus clone for extended periods. To be able to persist, the pathogen has to adapt to the hostile niche of the airways to counteract host defence, antibiotic therapy and the competition with coinfecting pathogens. S. aureus is equipped with many virulence factors including adhesins, toxins that are localized on the chromosome, on plasmids or are phage-related. S. aureus is especially versatile and adaptation and evolution of the pathogen occurs by the acquisition of new genes by horizontal gene transfer (HGT), changes in nucleotides (single nucleotide variations, SNVs) that can cause a selective advantage for the bacteria and become fixed in subpopulations. Methicillin-resistant S. aureus are a special threat to CF patients due to the more severe lung disease occurring in infected patients. Today, with decreasing costs for sequencing, more and more studies using S. aureus isolates cultured from CF patients are being published, which use whole genome sequencing (WGS), multilocus sequence typing (MLST) or spa-sequence typing (spa-typing) to follow the population dynamics of S. aureus, elucidate the underlying mechanisms of phenotypic variants, newly acquired resistance or adaptation to the host response in this particular niche. In the first part of this review, an introduction to the genetic make-up and the pathogenesis of S. aureus with respect to CF is provided. The second part presents an overview of recent studies and their findings using genotypic methods such as single or multilocus sequencing and whole genome sequencing, which identify factors contributing to the adaptation of S. aureus and its evolution in the airways of individuals with CF.

Outbreak of Pseudomonas aeruginosa producing VIM carbapenemase in an intensive care unit and its termination by implementation of waterless patient care

BACKGROUND

Long-term outbreaks of multidrug-resistant Gram-negative bacilli related to hospital-building water systems have been described. However, successful mitigation strategies have rarely been reported. In particular, environmental disinfection or replacement of contaminated equipment usually failed to eradicate environmental sources of Pseudomonas aeruginosa.

METHODS

We report the investigation and termination of an outbreak of P. aeruginosa producing VIM carbapenemase (PA-VIM) in the adult intensive care unit (ICU) of a Swiss tertiary care hospital with active case finding, environmental sampling and whole genome sequencing (WGS) of patient and environmental strains. We also describe the implemented control strategies and their effectiveness on eradication of the environmental reservoir.

RESULTS

Between April 2018 and September 2020, 21 patients became either infected or colonized with a PA-VIM strain. For 16 of them, an acquisition in the ICU was suspected. Among 131 environmental samples collected in the ICU, 13 grew PA-VIM in sink traps and drains. WGS confirmed the epidemiological link between clinical and environmental strains and the monoclonal pattern of the outbreak. After removing sinks from patient rooms and implementation of waterless patient care, no new acquisition was detected in the ICU within 8 months after the intervention.

DISCUSSION

Implementation of waterless patient care with removal of the sinks in patient rooms was successful for termination of a PA-VIM ICU outbreak linked to multiple environmental water sources. WGS provides highly discriminatory accuracy to investigate environment-related outbreaks.

Techniques in bacterial strain typing: past, present, and future

PURPOSE OF REVIEW

The advancement of molecular techniques such as whole-genome sequencing (WGS) has revolutionized the field of bacterial strain typing, with important implications for epidemiological surveillance and outbreak investigations. This review summarizes state-of-the-art techniques in strain typing and examines barriers faced by clinical and public health laboratories in implementing these new methodologies.

RECENT FINDINGS

WGS-based methodologies are on track to become the new 'gold standards' in bacterial strain typing, replacing traditional methods like pulsed-field gel electrophoresis and multilocus sequence typing. These new techniques have an improved ability to identify genetic relationships among organisms of interest. Further, advances in long-read sequencing approaches will likely provide a highly discriminatory tool to perform pangenome analyses and characterize relevant accessory genome elements, including mobile genetic elements carrying antibiotic resistance determinants in real time. Barriers to widespread integration of these approaches include a lack of standardized workflows and technical training.

SUMMARY

Genomic bacterial strain typing has facilitated a paradigm shift in clinical and molecular epidemiology. The increased resolution that these new techniques provide, along with epidemiological data, will facilitate the rapid identification of transmission routes with high confidence, leading to timely and effective deployment of infection control and public health interventions in outbreak settings.

Whole Genome-Based Characterization of Listeria monocytogenes Isolates Recovered From the Food Chain in South Africa

Listeria monocytogenes is an important foodborne pathogen which has the ability to adapt and survive in food and food processing facilities where it can persist for years. In this study, a total of 143 L. monocytogenes isolates in South Africa (SA) were characterized for their strain’s genetic relatedness, virulence profiles, stress tolerance and resistance genes associated with L. monocytogenes. The Core Genome Multilocus Sequence Typing (cgMLST) analysis revealed that the most frequent serogroups were IVb and IIa; Sequence Types (ST) were ST204, ST2, and ST1; and Clonal Complexes (CC) were CC204, CC1, and CC2. Examination of genes involved in adaptation and survival of L. monocytogenes in SA showed that ST1, ST2, ST121, ST204, and ST321 are well adapted in food processing environments due to the significant over-representation of Benzalkonium chloride (BC) resistance genes (bcrABC cassette, ermC, mdrL and Ide), stress tolerance genes (SSI-1 and SSI-2), Prophage (φ) profiles (LP_101, vB LmoS 188, vB_LmoS_293, and B054 phage), plasmids profiles (N1-011A, J1776, and pLM5578) and biofilm formation associated genes. Furthermore, the L. monocytogenes strains that showed hyper-virulent potential were ST1, ST2 and ST204, and hypo-virulent were ST121 and ST321 because of the presence and absence of major virulence factors such as LIPI-1, LIPI-3, LIPI-4 and the internalin gene family members including inlABCEFJ. The information provided in this study revealed that hyper-virulent strains ST1, ST2, and ST204 could present a major public health risk due to their association with meat products and food processing environments in SA.

Whole Genome Sequencing Applied to Pathogen Source Tracking in Food Industry: Key Considerations for Robust Bioinformatics Data Analysis and Reliable Results Interpretation

Whole genome sequencing (WGS) has arisen as a powerful tool to perform pathogen source tracking in the food industry thanks to several developments in recent years. However, the cost associated to this technology and the degree of expertise required to accurately process and understand the data has limited its adoption at a wider scale. Additionally, the time needed to obtain actionable information is often seen as an impairment for the application and use of the information generated via WGS. Ongoing work towards standardization of wet lab including sequencing protocols, following guidelines from the regulatory authorities and international standardization efforts make the technology more and more accessible. However, data analysis and results interpretation guidelines are still subject to initiatives coming from distinct groups and institutions. There are multiple bioinformatics software and pipelines developed to handle such information. Nevertheless, little consensus exists on a standard way to process the data and interpret the results. Here, we want to present the constraints we face in an industrial setting and the steps we consider necessary to obtain high quality data, reproducible results and a robust interpretation of the obtained information. All of this, in a time frame allowing for data-driven actions supporting factories and their needs.