A total of 66 Listeria monocytogenes (Lm) isolated from 2013 to 2018 in a small-scale meat processing plant and a dairy facility of Central Italy were studied. Whole Genome Sequencing and bioinformatics analysis were used to assess the genetic relationships between the strains and investigate persistence and virulence abilities. The biofilm forming-ability was assessed in vitro. Cluster analysis grouped the Lm from the meat plant into three main clusters: two of them, both belonging to CC9, persisted for years in the plant and one (CC121) was isolated in the last year of sampling. In the dairy facility, all the strains grouped in a CC2 four-year persistent cluster. All the studied strains carried multidrug efflux-pumps genetic determinants (sugE, mdrl, lde, norM, mepA). CC121 also harbored the Tn6188 specific for tolerance to Benzalkonium Chloride. Only CC9 and CC121 carried a Stress Survival Islet and presented high-level cadmium resistance genes (cadA1C1) carried by different plasmids. They showed a greater biofilm production when compared with CC2. All the CC2 carried a full-length inlA while CC9 and CC121 presented a Premature Stop Codon mutation correlated with less virulence. The hypo-virulent clones CC9 and CC121 appeared the most adapted to food-processing environments; however, even the hyper-virulent clone CC2 warningly persisted for a long time. The identification of the main mechanisms promoting Lm persistence in a specific food processing plant is important to provide recommendations to Food Business Operators (FBOs) in order to remove or reduce resident Lm.

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

Listeria monocytogenes (Lm) is a foodborne pathogen causing listeriosis. Invasive forms of the disease mainly manifest as septicaemia, meningitis and maternal-neonatal infections. Lm-associated respiratory infections are very rare and little known. We reported two Lm respiratory infection cases occurred in Central Italy during the summer of 2020, in the midst of the SARS-CoV2 pandemic. In addition to collect the epidemiological and clinical characteristics of the patients, we used Whole Genome Sequencing to study the genomes of the Lm isolates investigating their virulence and antimicrobial profiles and the presence of genetic mobile elements. Both the strains belonged to hypervirulent MLST clonal complexes (CC). In addition to the Listeria Pathogenicity Island 1 (LIPI-1), the CC1 strain also carried LIPI-3 and the CC4 both LIPI-3 and LIPI-4. Genetic determinants for antimicrobial and disinfectants resistance were found. The CC1 genome presented prophage sequences but they did not interrupt the comK gene, involved in the phagosomal escape of Lm. None of the strains carried plasmids. Lm is an important, although rare, opportunistic pathogen for respiratory tract and lung infections. To avoid dangerous diagnostic delays of these severe clinical forms, it is important to sensitize hospital laboratories to this rare manifestation of listeriosis considering Lm in the differential diagnosis of respiratory infections.

Salmonella enterica serovar Infantis is one of the five main causes of human salmonellosis in the European Union (EU) and in recent years, has been increasingly reported to carry multiple antimicrobial resistance determinants, including extended-spectrum beta-lactamase (ESBL) genes. In our study, we used WGS-based tools to characterize S. Infantis strains circulating in the Abruzzo and Molise regions of Italy between 2017 and 2020 and compared this local dataset to the S. Infantis population present in Italy over the last two decades. Phylogenetic analyses demonstrated that the majority of strains isolated from poultry and turkeys from Abruzzo and Molise were closely related and belonged to one of the two main genetic clusters present in Italy, which were grouped predominantly as ESBL-producing strains that harbored pESI-like plasmid. We showed that 60% of the local strains carried multiple antibiotic resistance genes, including ESBL gene blaCTX–M–1 as well as aadA1, dfrA1, dfrA14, sul1, and tet(A) genes present on the pESI-like megaplasmid. The analysis of strains from Abruzzo and Molise and the publicly available Italian S. Infantis sequences revealed a dramatic increase in the number of identified AMR genes in the strains isolated after 2011. Moreover, the number of strains resistant to five or more antibiotic classes increased from 20–80% in the last decade likely due to the acquisition of the megaplasmid. The persistence of the ESBL-producing and the multidrug-resistant (MDR) clone of S. Infantis in poultry populations in Italy and in Europe requires rapid and efficient intervention strategies to prevent further expansion of the clone.

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.