Pan-pathogen deep sequencing of nosocomial bacterial pathogens in Italy in spring 2020: a prospective cohort study.

BACKGROUND: Nosocomial infections pose a considerable risk to patients who are susceptible, and this is particularly acute in intensive care units when hospital-associated bacteria are endemic. During the first wave of the COVID-19 pandemic, the surge of patients presented a significant obstacle to the effectiveness of infection control measures. We aimed to assess the risks and extent of nosocomial pathogen transmission under a high patient burden by designing a novel bacterial pan-pathogen deep-sequencing approach that could be integrated with standard clinical surveillance and diagnostics workflows. METHODS: We did a prospective cohort study in a region of northern Italy that was severely affected by the first wave of the COVID-19 pandemic. Inpatients on both ordinary and intensive care unit (ICU) wards at the San Matteo hospital, Pavia were sampled on multiple occasions to identify bacterial pathogens from respiratory, nasal, and rectal samples. Diagnostic samples collected between April 7 and May 10, 2020 were cultured on six different selective media designed to enrich for Acinetobacter baumannii, Escherichia coli, Enterococcus faecium, Enterococcus faecalis, Klebsiella pneumoniae, Pseudomonas aeruginosa, Staphylococcus aureus, and Streptococcus pneumoniae, and DNA from each plate with positive growth was deep sequenced en masse. We used mSWEEP and mGEMS to bin sequencing reads by sequence cluster for each species, followed by mapping with snippy to generate high quality alignments. Antimicrobial resistance genes were detected by use of ARIBA and CARD. Estimates of hospital transmission were obtained from pairwise bacterial single nucleotide polymorphism distances, partitioned by within-patient and between-patient samples. Finally, we compared the accuracy of our binned Acinetobacter baumannii genomes with those obtained by single colony whole-genome sequencing of isolates from the same hospital. FINDINGS: We recruited patients from March 1 to May 7, 2020. The pathogen population among the patients was large and diverse, with 2148 species detections overall among the 2418 sequenced samples from the 256 patients. In total, 55 sequence clusters from key pathogen species were detected at least five times. The antimicrobial resistance gene prevalence was correspondingly high, with key carbapenemase and extended spectrum ß-lactamase genes detected in at least 50 (40%) of 125 patients in ICUs. Using high-resolution mapping to infer transmission, we established that hospital transmission was likely to be a significant mode of acquisition for each of the pathogen species. Finally, comparison with single colony Acinetobacter baumannii genomes showed that the resolution offered by deep sequencing was equivalent to single-colony sequencing, with the additional benefit of detection of co-colonisation of highly similar strains. INTERPRETATION: Our study shows that a culture-based deep-sequencing approach is a possible route towards improving future pathogen surveillance and infection control at hospitals. Future studies should be designed to directly compare the accuracy, cost, and feasibility of culture-based deep sequencing with single colony whole-genome sequencing on a range of bacterial species. FUNDING: Wellcome Trust, European Research Council, Academy of Finland Flagship program, Trond Mohn Foundation, and Research Council of Norway.

Molecular characterization of emerging Echovirus 11 (E11) shed light on the recombinant origin of a variant associated with severe hepatitis in neonates.

Echovirus 11 (E11) has gained attention owing to its association with severe neonatal infections. Due to the limited data available, the World Health Organization (WHO) considers public health risk to the general population to be low. The present study investigated the genetic variation and molecular evolution of E11 genomes collected from May to December 2023. Whole genome sequencing (WGS) was performed for 16 E11 strains. Phylogenetic analysis on WG showed how all Italian strains belonged to genogroup D5, similarly to other E11 strains recently reported in France and Germany all together aggregated into separate clusters. A cluster-specific recombination pattern was also identified using phylogenetic analysis of different genome regions. Echovirus 6 was identified as the major recombinant virus in 3Cpro and 3Dpol regions. The molecular clock analysis revealed that the recombination event probably occurred in June 2018 (95% HPD interval: Jan 2016-Jan 2020). Shannon entropy analyses, within P1 region, showed how 11 amino acids exhibited relatively high entropy. Five of them were exposed on the canyon region which is responsible for receptor binding with the neonatal Fc receptor. The present study showed the recombinant origin of a new lineage of E11 associated with severe neonatal infections.

Bloodstream Infection Caused by Erysipelothrix rhusiopathiae in an Immunocompetent Patient.

Erysipelothrix rhusiopathiae is a facultative anaerobe Gram-positive bacillus, which is considered a zoonotic pathogen. E. rhusiopathiae causes erysipeloid, mainly in occupational groups such as veterinarians, slaughterhouse workers, farmers, and fishermen. Two cutaneous forms (localised and generalised) and a septicaemic form have been described. Here, we report the isolation of a strain of E. rhusiopathiae from a 56-year-old immunocompetent obese male admitted to Fondazione IRCCS Policlinico San Matteo Pavia (Italy). Blood cultures were collected and Gram-positive bacilli were observed. E. rhusiopathiae grew and was identified. Antimicrobial susceptibility tests were performed and interpreted with EUCAST breakpoints (PK-PD). The strain was susceptible to all the antibiotics tested, while it was intrinsically resistant to vancomycin. The clinical diagnosis of E. rhusiopathiae can be challenging, due to the broad spectrum of symptoms and potential side effects, including serious systemic infections such as heart diseases. In the case described, bacteraemia caused by E. rhusiopathiae was detected in a immunocompetent patient. Bacteraemia caused by E. rhusiopathiae is rare in immunocompetent people and blood cultures were proven to be essential for the diagnosis and underdiagnosis of this pathogen, which is possible due to its resemblance to other clinical manifestations.

Impact of Whole Genome Sequencing to investigate transmission of Serratia marcescens in Neonatal Intensive Care Unit.

Newborns admitted to neonatal intensive care units (NICU) are at increased risk of health care-associated infections. Serratia marcescens represent the third most common pathogen in NICU outbreaks. Here we present an outbreak investigation performed using Whole Genome Sequencing (WGS) analyses and the control measures implemented to limit the spread of S. marcescens in the NICU of an Italian hospital. In February 2023 S. marcescens was isolated from six newborns, when in 2022 this pathogen was isolated only from two samples in the same ward. Measures for infection prevention were adopted. Routinary surveillance screening, performed with rectal swabs collected at admission and weekly thereafter, was implemented to search for S. marcescens presence. Environmental samples were collected. All the isolates, obtained from the conjunctival swab of six newborns, from rectal swab of two newborns who did not develop infections, as well as from the aerators of two faucets, were sequenced. WGS analyses showed no correlation between the isolates from newborns and environmental isolates. The implementation of the measures for infection prevention and control had enabled us to successfully control the outbreak within a short period. WGS analyses proved to be crucial in outbreak investigation to limit the spreading of the pathogens.

Enterobacter asburiae ST229: an emerging carbapenemases producer.

Enterobacter asburiae, member of the Enterobacter cloacae complex (ECC) group, shows an increasing clinical relevance being responsible for infections like pneumonia, urinary tract infections and septicemia. The aim of the present study was the investigation of the genomic features of two XDR E. asburiae ST229 clinical strains co-carrying blaNDM-1 and blaVIM-1 determinants, collected in October 2021 and in June 2022, respectively. Two E. asburiae strains were collected from rectal swabs of as many patients admitted to the cardiopulmonary intensive care unit of Fondazione I.R.C.C.S. "Policlinico San Matteo" in Pavia, Italy. Based on the antibiotic susceptibility profile results, both isolates showed an XDR phenotype, retaining susceptibility only to fluoroquinolones. Both isolates shared identical resistome, virulome, plasmid content, and belonged to ST229, a rarely reported sequence type. They co-harbored blaNDM-1 and blaVIM-1 genes, that resulted located on transferable plasmids by conjugation and transformation. Moreover, both strains differed in 24 SNPs and showed genetic relatedness with E. asburiae ST709 and ST27. We described the first case of ST229 E. asburiae co-harboring blaNDM-1 and blaVIM-1 in Italy. This study points out the emergence of carbapenemases in low-risk pathogens, representing a novel challenge for public health, that should include such types of strains in dedicated surveillance programs. Antimicrobial susceptibility testing was carried out using Thermo Scientific™ Sensititre™ Gram Negative MIC Plates DKMGN. Both strains underwent whole-genome sequencing (WGS) using Illumina Miseq platform. Resistome, plasmidome, virulome, MLST, plasmid MLST and a SNPs-based phylogenetic tree were in silico determined.

Dynamics of viral DNA shedding and culture viral DNA positivity in different clinical samples collected during the 2022 mpox outbreak in Lombardy, Italy.

BACKGROUND: Mpox virus (MPXV) has recently spread outside of sub-Saharan Africa. This large multicentre study was conducted in Lombardy, the most densely populated Italian region accounting for more than 40% of Italian cases. The present study aims to: i) evaluate the presence and the shedding duration of MPXV DNA in different body compartments correlating the MPXV viability with the time to onset of symptoms; ii) provide evidence of MPXV persistence in different body compartment as a source of infection and iii) characterize the MPXV evolution by whole genome sequencing (WGS) during the outbreak occurred in Italy. MATERIAL AND METHODS: The study included 353 patients with a laboratory-confirmed diagnosis of MPXV infection screened in several clinical specimens in the period May 24th - September 1st, 2022. Viral isolation was attempted from different biological matrices and complete genome sequencing was performed for 61 MPXV strains. RESULTS: MPXV DNA detection was more frequent in the skin (94.4%) with the longest median time of viral clearance (16 days). The actively-replicating virus in cell culture was obtained for 123/377 (32.6%) samples with a significant higher viral quantity on isolation positive samples (20 vs 31, p < 0.001). The phylogenetic analysis highlighted the high genetic identity of the MPXV strains collected, both globally and within the Lombardy region. CONCLUSION: Skin lesion is gold standard material and the high viral load and the actively-replicating virus observed in genital sites confirms that sexual contact plays a key role in the viral transmission.

Clinical isolates of ST131 blaOXA-244-positive Escherichia coli, Italy, December 2022 to July 2023.

The dissemination of carbapenemase-producing Escherichia coli, although still at low level, should be continuously monitored. OXA-244 is emerging in Europe, mainly in E. coli. In Italy, this carbapenemase was reported from an environmental river sample in 2019. We report clinical isolates of OXA-244-producing ST131 E. coli in four patients admitted to an acute care hospital in Pavia, Italy. The association of this difficult-to-detect determinant with a globally circulating high-risk clone, ST131 E. coli, is of clinical relevance.

Concomitant Resistance to Cefiderocol and Ceftazidime/Avibactam in Two Carbapenemase-Producing Klebsiella pneumoniae Isolates from Two Lung Transplant Patients.

In this study, we present two cases of Klebsiella pneumoniae, one KPC-33- and one NDM-1-producing, showing resistance to cefiderocol and ceftazidime/avibactam, collected in the intensive care unit of a hospital in Northern Italy from two patients who had recently undergone lung transplantation. Whole-genome sequencing was performed to investigate the molecular features of these strains.

Rapid spread of a novel NDM-producing clone of Klebsiella pneumoniae CC147, Northern Italy, February to August 2023.

New Delhi metallo-beta-lactamase (NDM)-producing Klebsiella pneumoniae (Kp) ST147 caused a large multi-hospital outbreak in Italy from 2018 to 2021. We describe a new ST6668 NDM-producing Kp clone, belonging to CC147, which rapidly spread across hospitals in the Pavia province (Northern Italy) from February to August 2023. Genomic analyses revealed that ST6668 is different from ST147 and fast evolving. As shown here, genomic surveillance programmes are useful for tracking the spread of new clones with reduced susceptibility to most antibiotics.

P-DOR, an easy-to-use pipeline to reconstruct bacterial outbreaks using genomics.

SUMMARY: Bacterial Healthcare-Associated Infections (HAIs) are a major threat worldwide, which can be counteracted by establishing effective infection control measures, guided by constant surveillance and timely epidemiological investigations. Genomics is crucial in modern epidemiology but lacks standard methods and user-friendly software, accessible to users without a strong bioinformatics proficiency. To overcome these issues we developed P-DOR, a novel tool for rapid bacterial outbreak characterization. P-DOR accepts genome assemblies as input, it automatically selects a background of publicly available genomes using k-mer distances and adds it to the analysis dataset before inferring a Single-Nucleotide Polymorphism (SNP)-based phylogeny. Epidemiological clusters are identified considering the phylogenetic tree topology and SNP distances. By analyzing the SNP-distance distribution, the user can gauge the correct threshold. Patient metadata can be inputted as well, to provide a spatio-temporal representation of the outbreak. The entire pipeline is fast and scalable and can be also run on low-end computers. AVAILABILITY AND IMPLEMENTATION: P-DOR is implemented in Python3 and R and can be installed using conda environments. It is available from GitHub https://github.com/SteMIDIfactory/P-DOR under the GPL-3.0 license.

The COVID-19 Pandemic Sparked Off a Large-Scale Outbreak of Carbapenem-Resistant Acinetobacter baumannii from the Endemic Strains at an Italian Hospital.

Acinetobacter baumannii is a nosocomial pathogen that poses a serious threat due to the rise of incidence of multidrug-resistant (MDR) strains. During the COVID-19 pandemic, MDR A. baumannii clones have caused several outbreaks worldwide. Here, we describe a detailed investigation of an MDR A. baumannii outbreak that occurred at Policlinico San Matteo (Pavia, Italy). A total of 96 A. baumannii strains, isolated between January and July 2020 from 41 inpatients (both SARS-CoV-2 positive and negative) in different wards, were characterized by phenotypic and genomic analyses combining Illumina and Nanopore sequencing. Antibiotic susceptibility testing revealed that all isolates were resistant to carbapenems, and the sequence analysis attributed this to the carbapenemase gene blaOXA-23. Virulence factor screening unveiled that all strains carried determinants for biofilm formation, while plasmid analysis revealed the presence of two plasmids, one of which was ~100 kbp long and encoded a phage sequence. A core genome-based phylogeny was inferred to integrate outbreak strain genomes with background genomes from public databases and the local surveillance program. All strains belonged to the globally disseminated sequence type 2 (ST2) clone and were mainly divided into two clades. Isolates from the outbreak clustered with surveillance isolates from 2019, suggesting that the outbreak was caused by two strains that were already circulating in the hospital before the start of the pandemic. The intensive spread of A. baumannii in the hospital was enhanced by the extreme emergency situation of the first COVID-19 pandemic wave that resulted in reduced attention to infection prevention and control practices. IMPORTANCE The COVID-19 pandemic, especially during the first wave, posed a great challenge to the hospital management and generally promoted nosocomial pathogen dissemination. MDR A. baumannii can easily spread and persist for a long time on surfaces, causing outbreaks in health care settings. Infection prevention and control practices, epidemiological surveillance, and microbiological screening are fundamental in order to control such outbreaks. Here, we sequenced the genomes of 96 isolates from an outbreak of MDR A. baumannii strains using both short- and long-read technology in order to reconstruct the outbreak events in fine detail. The sequence data demonstrated that two endemic clones of MDR A. baumannii were the source of this large hospital outbreak during the first COVID-19 pandemic wave, confirming the effect of COVID-19 emergency disrupting the protection provided by the use of the standard prevention procedures.

Surveillance in a Neonatal Intensive Care Unit Allowed the Isolation of a Strain of VIM-Producing Pantoea brenneri.

Here, we describe the isolation of a strain of the genus Pantoea encoding a VIM carbapenemase, the first to our knowledge. The strain, isolated from a rectal swab of a 10-day-old newborn admitted to a neonatal intensive care unit (NICU), was identified through whole-genome sequencing analyses as Pantoea brenneri. The strain harbored the carbapenemases gene blaVIM-1. The prompt application of contact measures and the isolation of the newborn prevented the dissemination of VIM-producing P. brenneri and of the plasmid carrying the VIM-1 gene to other newborns.

Comparative analysis of SARS-CoV-2 quasispecies in the upper and lower respiratory tract shows an ongoing evolution in the spike cleavage site.

Studies are needed to better understand the genomic evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This study aimed to describe viral quasispecies population of upper and lower respiratory tract by next-generation sequencing in patients admitted to intensive care unit. A deep sequencing of the S gene of SARS-CoV-2 from 109 clinical specimens, sampled from the upper respiratory tract (URT) and lower respiratory tract (LRT) of 77 patients was performed. A higher incidence of non-synonymous mutations and indels was observed in the LRT among minority variants. This might be explained by the ability of the virus to invade cells without interacting with ACE2 (e.g. exploiting macrophage phagocytosis). Minority variants are highly concentrated around the gene portion encoding for the Spike cleavage site, with a higher incidence in the URT; four mutations are highly recurring among samples and were found associated with the URT. Interestingly, 55.8% of minority variants detected in this locus were T>G and G>T transversions. Results from this study evidenced the presence of selective pressure and suggest that an evolutionary process is still ongoing in one of the crucial sites of spike protein associated with the spillover to humans.

Isolation of a Colistin-Susceptible MDR Pantoea calida Harboring the mcr-9 Gene Suggests the Silent Spread of the Resistance Factor.

Pantoea spp. are bacteria that are often detected in the environment and as symbionts of arthropods. They sporadically cause infections in humans and recently extended spectrum beta-lactamases (ESBL)- and carbapenemase-producing strains have started to emerge. In this study, we report the isolation and the complete genome sequence of a strain of Pantoea calida encoding the colistin-resistance gene mcr-9. The strain was isolated from a preterm newborn in a neonatal pathology ward. On clinical examination, his vital signs were normal and blood culture was negative. Rectal swab screening for ESBL-producing Enterobacterales allowed to isolate the bacterium, and a complete genome was obtained using both short and long read sequencing. The mcr-9 gene was found to be encoded on a IncHI2 superplasmid, which confers resistance to six classes of antibiotics, including beta lactams (ESBL). Despite the presence of mcr-9, the isolate retains susceptibility to colistin, which could be explained by the absence of compatible regulatory genes (qseBC) from the genome. The presence of the resistance gene is undetectable with the routine clinical procedures, that is, phenotypic tests. This suggests that a silent spread might be ongoing in the ward. To our knowledge, this is the first description of an MDR P. calida and of a Pantoea spp. encoding any mobile colistin resistance gene.

Spread of multiple SARS-CoV-2 lineages April-August 2020 anticipated the second pandemic wave in Lombardy (Italy).

During the early phase of the pandemic (20 February-4 April 2020), we have investigated the temporal and geographical evolution of the virus in Lombardy showing the circulation of at least seven lineages distributed differently in the Region. In the present study, the molecular epidemiology of SARS-CoV-2 was monitored in a period between two pandemic waves in order to track the circulation of new variants (April-August 2020). A great majority of SARS-CoV-2 strains (70.8%) belonged to lineages B, B.1, B.1.1 and B.1.1.1, and five strains belonging to four lineages were already reported in Italy (B.1.1.148, B.1.1.162, B.1.1.71, and B.1.425). In addition, 21 SARS-CoV-2 strains belonged to six lineages not previously observed in Italy were detected. No variants of concern were observed. A total of 152/1274 (11.3%) amino acid changes were observed among spike gene sequences and only 26/152 (17.1%) occurred in the receptor-binding domain region of the spike protein. Results of this study are indicative of ongoing transmission throughout the lockdown period, rather than re-introduction of novel lineages past lockdown. The use of molecular epidemiology in Italy should be promoted in order to provide additional understanding of the transmission of the disease and to have major effect on controlling the spread of disease.

Bispecific IgG neutralizes SARS-CoV-2 variants and prevents escape in mice.

Neutralizing antibodies that target the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein are among the most promising approaches against COVID-191,2. A bispecific IgG1-like molecule (CoV-X2) has been developed on the basis of C121 and C135, two antibodies derived from donors who had recovered from COVID-193. Here we show that CoV-X2 simultaneously binds two independent sites on the RBD and, unlike its parental antibodies, prevents detectable spike binding to the cellular receptor of the virus, angiotensin-converting enzyme 2 (ACE2). Furthermore, CoV-X2 neutralizes wild-type SARS-CoV-2 and its variants of concern, as well as escape mutants generated by the parental monoclonal antibodies. We also found that in a mouse model of SARS-CoV-2 infection with lung inflammation, CoV-X2 protects mice from disease and suppresses viral escape. Thus, the simultaneous targeting of non-overlapping RBD epitopes by IgG-like bispecific antibodies is feasible and effective, and combines the advantages of antibody cocktails with those of single-molecule approaches.

Epidemiological Characterization of Listeria monocytogenes Infections in Pavia Province in 2017 Reveals the Presence of Multiple Concurrently Circulating Strains.

Consumption of raw food, especially smoked fish, meat, soft cheeses, and vegetables, contaminated with Listeria monocytogenes can cause listeriosis, which can be invasive in pregnant women, elderly, and immunocompromised and diabetic patients. Through June to November of 2017, 11 patients developed invasive listeriosis in a small area of northern Italy. In the same period, 15 food samples (ready-to-eat seafood, raw vegetables, cheese samples, and salami) collected during the routine screening programs in the same area were found to be contaminated with L. monocytogenes. We characterized the isolates to determine the relatedness of L. monocytogenes strains isolated from patients and isolates from food samples and food-processing plants. Whole genome sequencing analysis showed that multiple L. monocytogenes strains were circulating in the area and no association was found between clinical and food isolates.

Genomic epidemiology of SARS-CoV-2 reveals multiple lineages and early spread of SARS-CoV-2 infections in Lombardy, Italy.

From February to April 2020, Lombardy (Italy) reported the highest numbers of SARS-CoV-2 cases worldwide. By analyzing 346 whole SARS-CoV-2 genomes, we demonstrate the presence of seven viral lineages in Lombardy, frequently sustained by local transmission chains and at least two likely to have originated in Italy. Six single nucleotide polymorphisms (five of them non-synonymous) characterized the SARS-CoV-2 sequences, none of them affecting N-glycosylation sites. The seven lineages, and the presence of local transmission clusters within three of them, revealed that sustained community transmission was underway before the first COVID-19 case had been detected in Lombardy.