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.

First imported case of Candida auris infection in Milan, Italy: genomic characterisation.

PURPOSE: Candida auris, an emerging multidrug-resistant yeast, has been reported worldwide. In Italy, the first case was reported in 2019. We describe the first case of C. auris, imported from Greece, in Milan, using whole genome sequencing to characterise mutations associated with antifungal resistance. CASE PRESENTATION: On October 2022 an 80-year-old Italian man was hospitalised in Greece. In the absence of clinical improvement, the patient was transferred to our hospital, in Italy, where blood culture resulted positive for C. auris. Despite therapy, the patient died of septic shock. In a phylogenetic analysis the genome was assigned to Clade I with strains from Kenya, United Arab Emirates and India. D1/D2 region resulted identical to a Greek strain, as for many other strains from different World regions, highlighting the diffusion of this strain. CONCLUSION: Importation of C. auris from abroad has been previously described. We report the first case of C. auris imported into Italy from Greece, according to phylogenetic analysis. This case reinforces the need for monitoring critically ill hospitalised patients also for fungi and addresses the need for the standardisation of susceptibility testing and strategies for diagnosis and therapy.

The evolution of intramitochondriality in Midichloria bacteria.

Midichloria spp. are intracellular bacterial symbionts of ticks. Representatives of this genus colonise mitochondria in the cells of their hosts. To shed light on this unique interaction we evaluated the presence of an intramitochondrial localization for three Midichloria in the respective tick host species and generated eight high-quality draft genomes and one closed genome, showing that this trait is non-monophyletic, either due to losses or multiple acquisitions. Comparative genomics supports the first hypothesis, as the genomes of non-mitochondrial symbionts are reduced subsets of those capable of colonising the organelles. We detect genomic signatures of mitochondrial tropism, including the differential presence of type IV secretion system and flagellum, which could allow the secretion of unique effectors and/or direct interaction with mitochondria. Other genes, including adhesion molecules, proteins involved in actin polymerisation, cell wall and outer membrane proteins, are only present in mitochondrial symbionts. The bacteria could use these to manipulate host structures, including mitochondrial membranes, to fuse with the organelles or manipulate the mitochondrial network.

Comparative analysis of two genomes of Chlamydia pecorum isolates from an Alpine chamois and a water buffalo.

BACKGROUND: To date, whole genome sequencing has been performed mainly for isolates of Chlamydia trachomatis, C. pneumoniae, C. psittaci and C. abortus, but only a few isolates of C. pecorum have been entirely sequenced and this makes it difficult to understand its diversity and population structure. In this study the genome of two C. pecorum strains isolated from the lung of an Alpine chamois affected with pneumonia (isolate PV7855) and the brain of a water buffalo affected with meningoencephalomyelitis (isolate PV6959), were completely sequenced with MiSeq system (Illumina) and analyzed in their most polymorphic regions. RESULTS: The genome length and GC content of the two isolates were found to be consistent with other C. pecorum isolates and the gene content of polymorphic membrane proteins and plasticity zone was found to be very similar. Some differences were observed in the phospholipase genes for both isolates and in the number of genes in the plasticity zone, such as the presence of some hypothetical proteins in PV6959, not present in any other genomes analyzed in this study. Interestingly, PV6959 possesses an extra pmp and has an incomplete tryptophan biosynthesis operon. Plasmids were detected in both isolates. CONCLUSIONS: Genome sequencing of the two C. pecorum strains did not reveal differences in length and GC content despite the origin from different animal species with different clinical disease. In the plasticity zone, the differences in the genes pattern might be related to the onset of specific symptoms or infection of specific hosts. The absence of a tryptophan biosynthesis pathway in PV6959 may suggest a strict relationship between C. pecorum and its host.

The Role of Virulence Factors in Neonatal Sepsis Caused by Enterobacterales: A Systematic Review.

Neonatal sepsis is a life-threatening condition with high mortality. Virulence determinants relevant in causing Gram-negative (GN) neonatal sepsis are still poorly characterized. A better understanding of virulence factors (VFs) associated with GN neonatal sepsis could offer new targets for therapeutic interventions. The aim of this review was to assess the role of GN VFs in neonatal sepsis. We primarily aimed to investigate the main VFs leading to adverse outcome and second to evaluate VFs associated with increased invasiveness/pathogenicity in neonates. MEDLINE, Embase, and Cochrane Library were systematically searched for studies reporting data on the role of virulome/VFs in bloodstream infections caused by Enterobacterales among neonates and infants aged 0-90 days. Twenty studies fulfilled the inclusion criteria. Only 4 studies reported data on the association between pathogen virulence determinants and neonatal mortality, whereas 16 studies were included in the secondary analyses. The quality of reporting was suboptimal in the great majority of the published studies. No consistent association between virulence determinants and GN strains causing neonatal sepsis was identified. Considerable heterogeneity was found in terms of VFs analysed and reported, included population and microbiological methods, with the included studies often showing conflicting data. This variability hampered the comparison of the results. In conclusions, pathogens responsible for neonatal sepsis are widely heterogenous and can use different pathways to develop invasive disease. The recent genome-wide approach needs to include multicentre studies with larger sample sizes, analyses of VF gene profiles instead of single VF genes, alongside a comprehensive collection of clinical information. A better understanding of the roles of virulence genes in neonatal GN bacteraemia may offer new vaccine targets and new markers of highly virulent strains. This information can potentially be used for screening and preventive interventions as well as for new targets for anti-virulence antibiotic-sparing therapies.

MeltingPlot, a user-friendly online tool for epidemiological investigation using High Resolution Melting data.

BACKGROUND: The rapid identification of pathogen clones is pivotal for effective epidemiological control strategies in hospital settings. High Resolution Melting (HRM) is a molecular biology technique suitable for fast and inexpensive pathogen typing protocols. Unfortunately, the mathematical/informatics skills required to analyse HRM data for pathogen typing likely limit the application of this promising technique in hospital settings. RESULTS: MeltingPlot is the first tool specifically designed for epidemiological investigations using HRM data, easing the application of HRM typing to large real-time surveillance and rapid outbreak reconstructions. MeltingPlot implements a graph-based algorithm designed to discriminate pathogen clones on the basis of HRM data, producing portable typing results. The tool also merges typing information with isolates and patients metadata to create graphical and tabular outputs useful in epidemiological investigations and it runs in a few seconds even with hundreds of isolates. AVAILABILITY: https://skynet.unimi.it/index.php/tools/meltingplot/ . CONCLUSIONS: The analysis and result interpretation of HRM typing protocols can be not trivial and this likely limited its application in hospital settings. MeltingPlot is a web tool designed to help the user to reconstruct epidemiological events by combining HRM-based clustering methods and the isolate/patient metadata. The tool can be used for the implementation of HRM based real time large scale surveillance programs in hospital settings.

Geographical reconstruction of the SARS-CoV-2 outbreak in Lombardy (Italy) during the early phase.

The first identification of autochthonous transmission of SARS-CoV-2 in Italy was documented by the Laboratory of Clinical Microbiology, Virology and Bioemergencies of L. Sacco Hospital (Milano, Italy) on 20th February 2020 in a 38 years old male patient, who was found positive for pneumonia at the Codogno Hospital. Thereafter Lombardy has reported the highest prevalence of COVID-19 cases in the country, especially in Milano, Brescia and Bergamo provinces. The aim of this study was to assess the potential presence of different viral clusters belonging to the six main provinces involved in Lombardy COVID-19 cases in order to highlight peculiar province-dependent viral characteristics. A phylogenetic analysis was conducted on 20 full length genomes obtained from patients addressing to several Lombard hospitals from February 20th to April 4th, 2020, aligned with 41 Italian viral genome assemblies available on GISAID database as of 30th March, 2020: two main monophyletic clades, containing 8 and 53 isolates, respectively, were identified. Noteworthy, Bergamo isolates mapped inside the small clade harbouring M gene D3G mutation. The molecular clock analysis estimated a cluster divergence approximately one month before the first patient identification, supporting the hypothesis that different SARS-CoV-2 strains had spread worldwide at different times, but their presence became evident only in late February along with Italian epidemic emergence. Therefore, this epidemiological reconstruction suggests that virus initial circulation in Lombardy was ascribable to multiple introduction. The phylogenetic reconstruction robustness, however, will be improved when more genomic sequences are available, in order to guarantee a complete epidemiological surveillance.

Genomic Characterization of an ST1153 PVL-producing Methicillin Resistant Staphylococcus aureus Clinical Isolate in Italy.

Methicillin-resistant Staphylococcus aureus (MRSA) clones are rapidly increasing beyond the hospital into the community, livestock farming and environmental settings. An Italian man, a professional diver working in Egypt, was admitted to Infectious Diseases Clinic-ASST Fatebenefratelli Sacco for ulcerative skin lesions. An MRSA strain was isolated from the lesions' purulent exudate and the nasal colonization was also ascertained. The strain, characterized by whole genome sequencing, resulted to be Panton-Valentine Leukocidin (PVL) positive, SCCmecI - spa-type t504, and belonging to the sequence type 1153, sporadically described worldwide.

Detection of ST1702 Escherichia coli blaNDM-5 and blaCMY-42 genes positive isolates from a Northern Italian hospital.

We describe two multi drug-resistant (MDR) carbapenemase-producing Escherichia coli clinical isolates from an acute hospital in Milan. Both strains, isolated from a surgical wound sample and a surveillance rectal swab respectively, were positive for a blaNDM-type gene by Xpert Carba-R test. The whole-genome sequence characterization disclosed several resistance determinants: blaNDM-5, blaCMY-42, blaTEM-198, rmtB, mphA. The two isolates belonged to phylogenetic group A, sequence type (ST) 1702 and serotype O89:H9. PCR-based replicon typing and conjugation assay demonstrated an IncI1 plasmid localization for both blaNDM-5 and blaCMY-42 genes. This is the first report of a ST1702 NDM-5 and CMY-42- producing E. coli clone in Italy.

Transcriptome of larvae representing the Rhipicephalus sanguineus complex.

Rhipicephalus sanguineus sensu lato (Ixodida: Ixodidae) is possibly the most widespread tick species worldwide, responsible for transmitting several vector-borne pathogens of medical and veterinary importance. Here, we explore the transcriptome of R. sanguineus s.l. larvae (Putignano strain). We sequenced total RNA from R. sanguineus s.l. larvae. A total of 15,566,986 short paired-end reads were de novo-assembled into 33,396 transcripts and then annotated and analyzed. Particular attention was paid to transcripts putatively encoding ATP-binding proteins, due to their importance as mechanisms of detoxification and acaricide resistance. Additionally, microsatellite loci were investigated, as these are useful markers for population genetic studies. The present data and analyses provide a comprehensive transcriptomic resource for R. sanguineus. The results presented here will aid further genetic and genomic studies of this important tick species.

Three cases of mcr-1-positive colistin-resistant Escherichia coli bloodstream infections in Italy, August 2016 to January 2017.

We describe three cases of bloodstream infection caused by colistin-resistant Escherichia coli in patients in a tertiary hospital in Italy, between August 2016 and January 2017. Whole genome sequencing detected the mcr-1 gene in three isolated strains belonging to different sequence types (STs). This occurrence of three cases with mcr-1-positive E. coli belonging to different STs in six months suggests a widespread problem in settings where high multidrug resistance is endemic such as in Italy.

Genomic epidemiology of Klebsiella pneumoniae in Italy and novel insights into the origin and global evolution of its resistance to carbapenem antibiotics.

Klebsiella pneumoniae is at the forefront of antimicrobial resistance for Gram-negative pathogenic bacteria, as strains resistant to third-generation cephalosporins and carbapenems are widely reported. The worldwide diffusion of these strains is of great concern due to the high morbidity and mortality often associated with K. pneumoniae infections in nosocomial environments. We sequenced the genomes of 89 K. pneumoniae strains isolated in six Italian hospitals. Strains were selected based on antibiotypes, regardless of multilocus sequence type, to obtain a picture of the epidemiology of K. pneumoniae in Italy. Thirty-one strains were carbapenem-resistant K. pneumoniae carbapenemase producers, 29 were resistant to third-generation cephalosporins, and 29 were susceptible to the aforementioned antibiotics. The genomes were compared to all of the sequences available in the databases, obtaining a data set of 319 genomes spanning the known diversity of K. pneumoniae worldwide. Bioinformatic analyses of this global data set allowed us to construct a whole-species phylogeny, to detect patterns of antibiotic resistance distribution, and to date the differentiation between specific clades of interest. Finally, we detected an ∼ 1.3-Mb recombination that characterizes all of the isolates of clonal complex 258, the most widespread carbapenem-resistant group of K. pneumoniae. The evolution of this complex was modeled, dating the newly detected and the previously reported recombination events. The present study contributes to the understanding of K. pneumoniae evolution, providing novel insights into its global genomic characteristics and drawing a dated epidemiological scenario for this pathogen in Italy.

Differential single nucleotide polymorphism-based analysis of an outbreak caused by Salmonella enterica serovar Manhattan reveals epidemiological details missed by standard pulsed-field gel electrophoresis.

We retrospectively analyzed a rare Salmonella enterica serovar Manhattan outbreak that occurred in Italy in 2009 to evaluate the potential of new genomic tools based on differential single nucleotide polymorphism (SNP) analysis in comparison with the gold standard genotyping method, pulsed-field gel electrophoresis. A total of 39 isolates were analyzed from patients (n=15) and food, feed, animal, and environmental sources (n=24), resulting in five different pulsed-field gel electrophoresis (PFGE) profiles. Isolates epidemiologically related to the outbreak clustered within the same pulsotype, SXB_BS.0003, without any further differentiation. Thirty-three isolates were considered for genomic analysis based on different sets of SNPs, core, synonymous, nonsynonymous, as well as SNPs in different codon positions, by Bayesian and maximum likelihood algorithms. Trees generated from core and nonsynonymous SNPs, as well as SNPs at the second and first plus second codon positions detailed four distinct groups of isolates within the outbreak pulsotype, discriminating outbreak-related isolates of human and food origins. Conversely, the trees derived from synonymous and third-codon-position SNPs clustered food and human isolates together, indicating that all outbreak-related isolates constituted a single clone, which was in line with the epidemiological evidence. Further experiments are in place to extend this approach within our regional enteropathogen surveillance system.

Tracking Nosocomial Klebsiella pneumoniae Infections and Outbreaks by Whole-Genome Analysis: Small-Scale Italian Scenario within a Single Hospital.

Multidrug-resistant (MDR) Klebsiella pneumoniae is one of the most important causes of nosocomial infections worldwide. After the spread of strains resistant to beta-lactams at the end of the previous century, the diffusion of isolates resistant to carbapenems and colistin is now reducing treatment options and the containment of infections. Carbapenem-resistant K. pneumoniae strains have spread rapidly among Italian hospitals, with four subclades of pandemic clonal group 258 (CG258). Here we show that a single Italian hospital has been invaded by three of these subclades within 27 months, thus replicating on a small scale the "Italian scenario." We identified a single clone responsible for an epidemic outbreak involving seven patients, and we reconstructed its star-like pattern of diffusion within the intensive care unit. This epidemiological picture was obtained through phylogenomic analysis of 16 carbapenem-resistant K. pneumoniae isolates collected in the hospital during a 27-month period, which were added to a database of 319 genomes representing the available global diversity of K. pneumoniae strains. Phenotypic and molecular assays did not reveal virulence or resistance determinants specific for the outbreak isolates. Other factors, rather than selective advantages, might have caused the outbreak. Finally, analyses allowed us to identify a major subclade of CG258 composed of strains bearing the yersiniabactin virulence factor. Our work demonstrates how the use of combined phenotypic, molecular, and whole-genome sequencing techniques can help to identify quickly and to characterize accurately the spread of MDR pathogens.

Matrix-assisted laser desorption ionization-time of flight and comparative genomic analysis of M-18 group a Streptococcus strains associated with an acute rheumatic fever outbreak in northeast Italy in 2012 and 2013.

Acute rheumatic fever (ARF) is a postsuppurative sequela caused by Streptococcus pyogenes infections affecting school-age children. We describe here the occurrence of an ARF outbreak that occurred in Bologna province, northeastern Italy, between November 2012 and May 2013. Molecular analysis revealed that ARF-related group A Streptococcus (GAS) strains belonged to the M-18 serotype, including subtypes emm18.29 and emm18.32. All M-18 GAS strains shared the same antigenic profile, including SpeA, SpeB, SpeC, SpeL, SpeM, and SmeZ. Matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) analysis revealed that M-18 GAS strains grouped separately from other serotypes, suggesting a different S. pyogenes lineage. Single nucleotide polymorphisms and phylogenetic analysis based on whole-genome sequencing showed that emm18.29 and emm18.32 GAS strains clustered in two distinct groups, highlighting genetic variations between these subtypes. Comparative analysis revealed a similar genome architecture between emm18.29 and emm18.32 strains that differed from noninvasive emm18.0 strains. The major sources of differences between M-18 genomes were attributable to the prophage elements. Prophage regions contained several virulence factors that could have contributed to the pathogenic potential of emm18.29 and emm18.32 strains. Notably, phage ΦSPBO.1 carried erythrogenic toxin A gene (speA1) in six ARF-related M-18 GAS strains but not in emm18.0 strains. In addition, a phage-encoded hyaluronidase gene (hylP.2) presented different variants among M-18 GAS strains by showing internal deletions located in the α-helical and TSßH regions. In conclusion, our study yielded insights into the genome structure of M-18 GAS strains responsible for the ARF outbreak in Italy, thus expanding our knowledge of this serotype.

Genetic barriers more than environmental associations explain Serratia marcescens population structure.

Bacterial species often comprise well-separated lineages, likely emerged and maintained by genetic isolation and/or ecological divergence. How these two evolutionary actors interact in the shaping of bacterial population structure is currently not fully understood. In this study, we investigate the genetic and ecological drivers underlying the evolution of Serratia marcescens, an opportunistic pathogen with high genomic flexibility and able to colonise diverse environments. Comparative genomic analyses reveal a population structure composed of five deeply-demarcated genetic clusters with open pan-genome but limited inter-cluster gene flow, partially explained by Restriction-Modification (R-M) systems incompatibility. Furthermore, a large-scale research on hundred-thousands metagenomic datasets reveals only a partial habitat separation of the clusters. Globally, two clusters only show a separate gene composition coherent with ecological adaptations. These results suggest that genetic isolation has preceded ecological adaptations in the shaping of the species diversity, an evolutionary scenario coherent with the Evolutionary Extended Synthesis.

Cultivation and sequencing-free protocol for Serratia marcescens detection and typing.

Serratia marcescens is an opportunistic pathogen that survives in inhospitable environments causing large outbreaks, particularly in neonatal intensive care units (NICUs). Genomic studies revealed that most S. marcescens nosocomial infections are caused by a specific clone (here "Infectious clone"). Whole genome sequencing (WGS) is the only portable method able to identify this clone, but it requires days to obtain results. We present a cultivation-free hypervariable-locus melting typing (HLMT) protocol for the fast detection and typing of S. marcescens, with 100% detection capability on mixed samples and a limit of detection that can reach the 10 genome copies. The protocol was able to identify the S. marcescens infectious clone with 97% specificity and 96% sensitivity when compared to WGS, yielding typing results portable among laboratories. The protocol is a cost and time saving method for S. marcescens detection and typing for large environmental/clinical surveillance screenings, also in low-middle income countries.

The genome of the heartworm, Dirofilaria immitis, reveals drug and vaccine targets.

The heartworm Dirofilaria immitis is an important parasite of dogs. Transmitted by mosquitoes in warmer climatic zones, it is spreading across southern Europe and the Americas at an alarming pace. There is no vaccine, and chemotherapy is prone to complications. To learn more about this parasite, we have sequenced the genomes of D. immitis and its endosymbiont Wolbachia. We predict 10,179 protein coding genes in the 84.2 Mb of the nuclear genome, and 823 genes in the 0.9-Mb Wolbachia genome. The D. immitis genome harbors neither DNA transposons nor active retrotransposons, and there is very little genetic variation between two sequenced isolates from Europe and the United States. The differential presence of anabolic pathways such as heme and nucleotide biosynthesis hints at the intricate metabolic interrelationship between the heartworm and Wolbachia. Comparing the proteome of D. immitis with other nematodes and with mammalian hosts, we identify families of potential drug targets, immune modulators, and vaccine candidates. This genome sequence will support the development of new tools against dirofilariasis and aid efforts to combat related human pathogens, the causative agents of lymphatic filariasis and river blindness.

Spatial and Temporal Analyses of the Event of Death for 1480 in Milan Using the Data Contained in the Sforza's Registers of the Dead.

Historical death registration was conducted primarily to assess the presence of plague. The Liber Mortuorum of Milan was one of Europe's first registers with many socio-demographical details. In this work, we consider 1480 to make spatial and temporal analyses of the event of death to find possible explanations about the events' distribution and the events' trend over time. The spatial analyses involved Moran's I, the LISA, and the heatmaps; the temporal analysis applied the Durbin-Watson test. All the analyses were conducted separately on all subjects (1813), children (765), and adults (1046). Contrade (districts) were considered for spatial analysis. Moran's I and the Durbin Watson test were significant on all subjects and children's analyses, and the LISA showed the same results for those groups. Children may significantly impact the distribution of death and the trend over time. At least half of the children were 0 years old, and survival in the very first childhood period was closely linked to the family, so that it could be a proxy of the conditions of an area.