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.

New Sequence Types and Antimicrobial Drug-Resistant Strains of Streptococcus suis in Diseased Pigs, Italy, 2017-2019.

Streptococcus suis is a pathogen associated with severe diseases in pigs and humans. Human infections have a zoonotic origin in pigs. To assess circulating strains, we characterized the serotypes, sequence types, and antimicrobial susceptibility of 78 S. suis isolates from diseased farmed pigs in Italy during 2017-2019. Almost 60% of infections were caused by serotypes 1/2 and 9. All but 1 of the serotype 2 and 1/2 isolates were confined to a single cluster, and serotype 9 isolates were distributed along the phylogenetic tree. Besides sequence type (ST) 1, the serotype 2 cluster included ST7, which caused severe human infections in China in 1998 and 2005. A large proportion of serotype 9 isolates, assigned to ST123, were resistant to penicillin. The emergence of this clone threatens the successful treatment of S. suis infection. Characterizing S. suis isolates from pigs will promote earlier detection of emerging clones.

SARS-CoV-2 surveillance in Italy through phylogenomic inferences based on Hamming distances derived from pan-SNPs, -MNPs and -InDels.

BACKGROUND: Faced with the ongoing global pandemic of coronavirus disease, the 'National Reference Centre for Whole Genome Sequencing of microbial pathogens: database and bioinformatic analysis' (GENPAT) formally established at the 'Istituto Zooprofilattico Sperimentale dell'Abruzzo e del Molise' (IZSAM) in Teramo (Italy) is in charge of the SARS-CoV-2 surveillance at the genomic scale. In a context of SARS-CoV-2 surveillance requiring correct and fast assessment of epidemiological clusters from substantial amount of samples, the present study proposes an analytical workflow for identifying accurately the PANGO lineages of SARS-CoV-2 samples and building of discriminant minimum spanning trees (MST) bypassing the usual time consuming phylogenomic inferences based on multiple sequence alignment (MSA) and substitution model. RESULTS: GENPAT constituted two collections of SARS-CoV-2 samples. The first collection consisted of SARS-CoV-2 positive swabs collected by IZSAM from the Abruzzo region (Italy), then sequenced by next generation sequencing (NGS) and analyzed in GENPAT (n = 1592), while the second collection included samples from several Italian provinces and retrieved from the reference Global Initiative on Sharing All Influenza Data (GISAID) (n = 17,201). The main results of the present work showed that (i) GENPAT and GISAID detected the same PANGO lineages, (ii) the PANGO lineages B.1.177 (i.e. historical in Italy) and B.1.1.7 (i.e. 'UK variant') are major concerns today in several Italian provinces, and the new MST-based method (iii) clusters most of the PANGO lineages together, (iv) with a higher dicriminatory power than PANGO lineages, (v) and faster that the usual phylogenomic methods based on MSA and substitution model. CONCLUSIONS: The genome sequencing efforts of Italian provinces, combined with a structured national system of NGS data management, provided support for surveillance SARS-CoV-2 in Italy. We propose to build phylogenomic trees of SARS-CoV-2 variants through an accurate, discriminant and fast MST-based method avoiding the typical time consuming steps related to MSA and substitution model-based phylogenomic inference.

Proficiency Testing of Virus Diagnostics Based on Bioinformatics Analysis of Simulated In Silico High-Throughput Sequencing Data Sets.

Quality management and independent assessment of high-throughput sequencing-based virus diagnostics have not yet been established as a mandatory approach for ensuring comparable results. The sensitivity and specificity of viral high-throughput sequence data analysis are highly affected by bioinformatics processing using publicly available and custom tools and databases and thus differ widely between individuals and institutions. Here we present the results of the COMPARE [Collaborative Management Platform for Detection and Analyses of (Re-)emerging and Foodborne Outbreaks in Europe] in silico virus proficiency test. An artificial, simulated in silico data set of Illumina HiSeq sequences was provided to 13 different European institutes for bioinformatics analysis to identify viral pathogens in high-throughput sequence data. Comparison of the participants' analyses shows that the use of different tools, programs, and databases for bioinformatics analyses can impact the correct identification of viral sequences from a simple data set. The identification of slightly mutated and highly divergent virus genomes has been shown to be most challenging. Furthermore, the interpretation of the results, together with a fictitious case report, by the participants showed that in addition to the bioinformatics analysis, the virological evaluation of the results can be important in clinical settings. External quality assessment and proficiency testing should become an important part of validating high-throughput sequencing-based virus diagnostics and could improve the harmonization, comparability, and reproducibility of results. There is a need for the establishment of international proficiency testing, like that established for conventional laboratory tests such as PCR, for bioinformatics pipelines and the interpretation of such results.

Outbreak of unusual Salmonella enterica serovar Typhimurium monophasic variant 1,4 [5],12:i:-, Italy, June 2013 to September 2014.

Monophasic variant of Salmonella enterica subspecies enterica serovar Typhimurium (monophasic S. Typhimurium), with antigenic structure 1,4,[5],12:i:-, appears to be of increasing importance in Europe. In Italy, monophasic S. Typhimurium represented the third most frequent Salmonella serovar isolated from human cases between 2004 and 2008. From June 2013 to October 2014, a total of 206 human cases of salmonellosis were identified in Abruzzo region (Central Italy). Obtained clinical isolates characterised showed S. Typhimurium 1,4,[5],12:i:- with sole resistance to nalidixic acid, which had never been observed in Italy in monophasic S. Typhimurium, neither in humans nor in animals or foods. Epidemiological, microbiological and environmental investigations were conducted to try to identify the outbreak source. Cases were interviewed using a standardised questionnaire and microbiological tests were performed on human as well as environmental samples, including samples from fruit and vegetables, pigs, and surface water. Investigation results did not identify the final vehicle of human infection, although a link between the human cases and the contamination of irrigation water channels was suggested.

webPDBinder: a server for the identification of ligand binding sites on protein structures.

The webPDBinder (http://pdbinder.bio.uniroma2.it/PDBinder) is a web server for the identification of small ligand-binding sites in a protein structure. webPDBinder searches a protein structure against a library of known binding sites and a collection of control non-binding pockets. The number of similarities identified with the residues in the two sets is then used to derive a propensity value for each residue of the query protein associated to the likelihood that the residue is part of a ligand binding site. The predicted binding residues can be further refined using conservation scores derived from the multiple alignment of the PFAM protein family. webPDBinder correctly identifies residues belonging to the binding site in 77% of the cases and is able to identify binding pockets starting from holo or apo structures with comparable performances. This is important for all the real world cases where the query protein has been crystallized without a ligand and is also difficult to obtain clear similarities with bound pockets from holo pocket libraries. The input is either a PDB code or a user-submitted structure. The output is a list of predicted binding pocket residues with propensity and conservation values both in text and graphical format.

Dynamic of SARS-CoV-2 variants circulation in Tunisian pediatric population, during successive waves, from March 2020 to September 2022.

The emergence of SARS-CoV-2 variants has led to several cases among children. However, limited information is available from North African countries. This study describes the SARS-CoV-2 strains circulating in Tunisian pediatric population during successive waves. A total of 447 complete sequences were obtained from individuals aged from 13 days to 18 years, between March 2020 and September 2022: 369 sequences generated during this study and 78 ones, available in GISAID, previously obtained from Tunisian pediatric patients. These sequences were compared with 354 and 274 ones obtained from Tunisian adults and a global dataset, respectively. The variant circulation dynamics of predominant variants were investigated during the study period using maximum-likelihood phylogenetic analysis. Among the studied population, adolescents were the predominant age group, comprising 55.26% of cases. Twenty-three lineages were identified; seven of which were not previously reported in Tunisia. Phylogenetic analysis showed a close relationship between the sequences from Tunisian adults and children. The connections of sequences from other countries were variable according to variants: close relationships were observed for Alpha, B1.160 and Omicron variants, while independent Tunisian clusters were observed for Delta and B.1.177 lineages. These findings highlight the pivotal role of children in virus transmission and underscore the impact of vaccination on virus spread. Vaccination of children, with booster doses, may be considered for better management of future emergences.

Phosfinder: a web server for the identification of phosphate-binding sites on protein structures.

Phosfinder is a web server for the identification of phosphate binding sites in protein structures. Phosfinder uses a structural comparison algorithm to scan a query structure against a set of known 3D phosphate binding motifs. Whenever a structural similarity between the query protein and a phosphate binding motif is detected, the phosphate bound by the known motif is added to the protein structure thus representing a putative phosphate binding site. Predicted binding sites are then evaluated according to (i) their position with respect to the query protein solvent-excluded surface and (ii) the conservation of the binding residues in the protein family. The server accepts as input either the PDB code of the protein to be analyzed or a user-submitted structure in PDB format. All the search parameters are user modifiable. Phosfinder outputs a list of predicted binding sites with detailed information about their structural similarity with known phosphate binding motifs, and the conservation of the residues involved. A graphical applet allows the user to visualize the predicted binding sites on the query protein structure. The results on a set of 52 apo/holo structure pairs show that the performance of our method is largely unaffected by ligand-induced conformational changes. Phosfinder is available at http://phosfinder.bio.uniroma2.it.

One health system supporting surveillance during COVID-19 epidemic in Abruzzo region, southern Italy.

The Istituti Zooprofilattici Sperimentali (IZSs) are public health institutes dealing with the aetiology and pathogenesis of infectious diseases of domestic and wild animals. During Coronavirus Disease 2019 epidemic, the Italian Ministry of Health appointed the IZSs to carry out diagnostic tests for the detection of SARS-CoV-2 in human samples. In particular, the IZS of Abruzzo and Molise (IZS-Teramo) was involved in the diagnosis of SARS-CoV-2 through testing nasopharyngeal swabs by Real Time RT-PCR. Activities and infrastructures were reorganised to the new priorities, in a "One Health" framework, based on interdisciplinary, laboratory promptness, accreditation of the test for the detection of the RNA of SARS-CoV-2 in human samples, and management of confidentiality of sensitive data. The laboratory information system - SILAB - was implemented with a One Health module for managing data of human origin, with tools for the automatic registration of information improving the quality of the data. Moreover, the "National Reference Centre for Whole Genome Sequencing of microbial pathogens - database and bioinformatics analysis" - GENPAT - formally established at the IZS-Teramo, developed bioinformatics workflows and IT dashboard with ad hoc surveillance tools to support the metagenomics-based SARS-CoV-2 surveillance, providing molecular sequencing analysis to quickly intercept the variants circulating in the area. This manuscript describes the One Health system developed by adapting and integrating both SILAB and GENPAT tools for supporting surveillance during COVID-19 epidemic in the Abruzzo region, southern Italy. The developed dashboard permits the health authorities to observe the SARS-CoV-2 spread in the region, and by combining spatio-temporal information with metagenomics provides early evidence for the identification of emerging space-time clusters of variants at the municipality level. The implementation of the One Health module was designed to be easily modelled and adapted for the management of other diseases and future hypothetical events of pandemic nature.

Whole-Genome Sequence and Assembly of Eight Africa Horse Sickness Virus Strains Collected in Namibia and South Africa.

In this report, we describe eight complete genome sequences of African horse sickness virus (AHSV) strains belonging to four different serotypes, namely, AHSV-5, AHSV-6, AHSV-8, and AHSV-9. Samples were collected in Namibia and South Africa from infected horses between 2000 and 2011. As expected, phylogenetic analyses of the variable outer capsid protein VP2 genomic sequences of AHSV-6 and AHSV-8 show higher nucleotide identity between the isolated viruses than that of the relevant reference strains. The full-genome sequence of AHSV will provide useful information on its geographical origin, and it will also be instrumental for comparing the distribution of the Namibian isolate with that of global isolates.

Evaluation of next generation sequencing approaches for SARS-CoV-2.

Within public health control strategies for SARS-CoV-2, whole genome sequencing (WGS) is essential for tracking viral spread and monitoring the emergence of variants which may impair the effectiveness of vaccines, diagnostic methods, and therapeutics. In this manuscript different strategies for SARS-CoV-2 WGS including metagenomic shotgun (SG), library enrichment by myBaits® Expert Virus-SARS-CoV-2 (Arbor Biosciences), nCoV-2019 sequencing protocol, ampliseq approach by Swift Amplicon® SARS-CoV-2 Panel kit (Swift Biosciences), and Illumina COVIDSeq Test (Illumina Inc.), were evaluated in order to identify the best approach in terms of results, labour, and costs. The analysis revealed that Illumina COVIDSeq Test (Illumina Inc.) is the best choice for a cost-effective, time-consuming production of consensus sequences.

Genomic and Antimicrobial Surveillance of Campylobacter Population in Italian Poultry.

Campylobacter is one of the most common foodborne diseases worldwide with increasing rates of antibiotic resistance. Most cases of campylobacteriosis can be traced back to the consumption of poultry meat. Despite many efforts to reduce contamination in farms and in slaughterhouses, the persistence of this pathogen in poultry products remains a problem. This study aimed to evaluate the genetic diversity and antibiotic resistance of 542 C. jejuni and C. coli in Italian poultry, in the framework of two National Monitoring Programs. Genomes were screened for antibiotic resistance, virulence determinants and contextualized within a global collection of C. jejuni. ST2116, ST2863 and ST 832 were the most prevalent and significantly associated with Italian poultry. A worrying increase in resistance to quinolones, fluoroquinolones and tetracycline was observed in C. jejuni, while an increased occurrence of multidrug resistant (MDR) strains and strains resistant to macrolides was detected in C. coli. Low resistance rates were found for aminoglycosides. Molecular resistance determinants were consistent with the phenotypic resistance for tetracycline and quinolones. In silico analysis revealed 119 genes associated with virulence factors, with a notably higher prevalence of some genes in ST2863 genomes. This study highlights the increased resistance to macrolides and the emergence of MDR strains for C. coli, the genetic basis of AMR and the predominance of two genotypes among Campylobacter strains isolated from the Italian poultry farms.

Spatial and temporal dynamics of West Nile virus between Africa and Europe.

It is unclear whether West Nile virus (WNV) circulates between Africa and Europe, despite numerous studies supporting an African origin and high transmission in Europe. We integrated genomic data with geographic observations and phylogenetic and phylogeographic inferences to uncover the spatial and temporal viral dynamics of WNV between these two continents. We focused our analysis towards WNV lineages 1 (L1) and 2 (L2), the most spatially widespread and pathogenic WNV lineages. Our study shows a Northern-Western African origin of L1, with back-and-forth exchanges between West Africa and Southern-Western Europe; and a Southern African origin of L2, with one main introduction from South Africa to Europe, and no back introductions observed. We also noticed a potential overlap between L1 and L2 Eastern and Western phylogeography and two Afro-Palearctic bird migratory flyways. Future studies linking avian and mosquito species susceptibility, migratory connectivity patterns, and phylogeographic inference are suggested to elucidate the dynamics of emerging viruses.

Reemergence of an atypical bluetongue virus strain in goats, Sardinia, Italy.

Bluetongue virus (BTV) is the etiologic agent of bluetongue, a WOAH (founded as Office International des Épizooties, OIE)-notifiable economically important disease of ruminants. BTV is transmitted by Culicoides biting midges and 24 different "classical" serotypes have been reported to date. In recent years, several putative novel BTV serotypes, often referred to as "atypical" BTVs, have been documented. These are characterized by unusual biological characteristics, most notably avirulence and vector-independent transmission. Here, we describe the recurrence of such an atypical virus strain BTV-X ITL2021 detected in goats six years after its first discovery in Sardinia, Italy. Combined serological and genome analysis results clearly suggest that the two strains belong to the same BTV serotype. However, unlike the 2015 strain, BTV-X ITL2021 was successfully isolated in BSR cell-culture allowing further serological characterization. Lastly, seropositivity for BTV-X ITL2021 was detected by virus-neutralization in approximately 74% of animals tested, suggesting that this atypical BTV serotype has been circulating undetected in asymptomatic animals for years.

SARS-CoV-2 in Namibian Dogs.

The pandemic of coronavirus disease 19 (COVID-19) has focused the attention of researchers, and especially public opinion, on the role of the human-animal-environment interface in disease emergence. At the beginning of the COVID-19 pandemic, media reports regarding the role of pets in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused significant concern and social anxiety. Although nowadays proven negligible in developed countries, essentially no studies have been performed in low-income African areas where companion animals are often raised differently from high income countries, and the contact patterns occurring in these scenarios could affect the epidemiological scenario. An extensive molecular biology survey was performed from March 2022 to September 2022 on Namibian dogs residing in urban and rural areas, showing a low but not negligible SARS-CoV-2 prevalence (1%; 95CI: 0.33-2.32%) of 5 out of 500. In only one instance (i.e., a 4-year-old female Labrador) was there a clear association that could be established between the infections of the owner and animal. In all other cases, no evidence of human infection could be obtained and no episodes of COVID-19 were reported by the owners. Although no consistent evidence of pet-to-pet transmission was proven in the present study, a cautionary principle suggests intensive and dedicated investigation into companion animal populations, especially when animal contact is frequent and a particularly susceptible population is present.

Epidemiological and genomic findings of the first documented Italian outbreak of SARS-CoV-2 Alpha variant of concern.

From 24 December 2020 to 8 February 2021, 163 cases of SARS-CoV-2 Alpha variant of concern (VOC) were identified in Chieti province, Abruzzo region. Epidemiological data allowed the identification of 14 epi-clusters. With one exception, all the epi-clusters were linked to the town of Guardiagrele: 149 contacts formed the network, two-thirds of which were referred to the family/friends context. Real data were then used to estimate transmission parameters. According to our method, the calculated Re(t) was higher than 2 before the 12 December 2020. Similar values were obtained from other studies considering Alpha VOC. Italian sequence data were combined with a random subset of sequences obtained from the GISAID database. Genomic analysis showed close identity between the sequences from Guardiagrele, forming one distinct clade. This would suggest one or limited unspecified viral introductions from outside to Abruzzo region in early December 2020, which led to the diffusion of Alpha VOC in Guardiagrele and in neighbouring municipalities, with very limited inter-regional mixing.

Emergence and Spread of SARS-CoV-2 Lineages B.1.1.7 and P.1 in Italy.

Italy's second wave of SARS-CoV-2 has hit hard, with more than three million cases and over 100,000 deaths, representing an almost ten-fold increase in the numbers reported by August 2020. Herein, we present an analysis of 6515 SARS-CoV-2 sequences sampled in Italy between 29 January 2020 and 1 March 2021 and show how different lineages emerged multiple times independently despite lockdown restrictions. Virus lineage B.1.177 became the dominant variant in November 2020, when cases peaked at 40,000 a day, but since January 2021 this is being replaced by the B.1.1.7 'variant of concern'. In addition, we report a sudden increase in another documented variant of concern-lineage P.1-from December 2020 onwards, most likely caused by a single introduction into Italy. We again highlight how international importations drive the emergence of new lineages and that genome sequencing should remain a top priority for ongoing surveillance in Italy.

Characterization and In Vitro Efficacy against Listeria monocytogenes of a Newly Isolated Bacteriophage, ɸIZSAM-1.

Listeria monocytogenes is a bacterial pathogen responsible of listeriosis, a disease that in humans is often related to the contamination of ready-to-eat foods. Phages are candidate biodecontaminants of pathogenic bacteria thanks to their ability to lyse prokaryotes while being safe for eukaryotic cells. In this study, ɸIZSAM-1 was isolated from the drain-waters of an Italian blue cheese plant and showed lytic activity against antimicrobial resistant Listeria monocytogenes strains. This phage was subjected to purification and in vitro efficacy tests. The results showed that at multiplicities of infection (MOIs) ≤ 1, phages were able to keep Listeria monocytogenes at low optical density values up to 8 h, with bacterial counts ranging from 1.02 to 3.96 log10 units lower than the control. Besides, ɸIZSAM-1 was further characterized, showing 25 principal proteins (sodium dodecyl sulfate polyacrylamide gel electrophoresis profile) and a genome of approximately 50 kilo base pairs. Moreover, this study describes a new approach to phage isolation for applications in Listeriamonocytogenes biocontrol in food production. In particular, the authors believe that the selection of phages from the same environments where pathogens live could represent a new approach to successfully integrating the control measures in an innovative, cost effective, safe and environmentally friendly way.

Whole-Genome Sequences of SARS-CoV-2 Lineage B.1.525 Strains (Variant η) Detected from Patients in the Abruzzo Region (Central Italy) during Spring 2021.

Novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants are emerging worldwide. Here, we report the complete genome sequences of 13 severe acute SARS-CoV-2 strains belonging to lineage B.1.525 (variant η).

Multiple detection and spread of novel strains of the SARS-CoV-2 B.1.177 (B.1.177.75) lineage that test negative by a commercially available nucleocapsid gene real-time RT-PCR.

Several lineages of SARS-CoV-2 are currently circulating worldwide. During SARS-CoV-2 diagnostic activities performed in Abruzzo region (central Italy) several strains belonging to the B.1.177.75 lineage tested negative for the N gene but positive for the ORF1ab and S genes (+/+/- pattern) by the TaqPath COVID-19 CE-IVD RT-PCR Kit manufactured by Thermofisher. By sequencing, a unique mutation, synonymous 28948C > T, was found in the N-negative B.1.177.75 strains. Although we do not have any knowledge upon the nucleotide sequences of the primers and probe adopted by this kit, it is likely that N gene dropout only occurs when 28948C > T is coupled with 28932C > T, this latter present, in turn, in all B.1.177.75 sequences available on public databases. Furthermore, epidemiological analysis was also performed. The majority of the N-negative B.1.177.75 cases belonged to two clusters apparently unrelated to each other and both clusters involved young people. However, the phylogeny for sequences containing the +/+/- pattern strongly supports a genetic connection and one common source for both clusters. Though, genetic comparison suggests a connection rather than indicating the independent emergence of the same mutation in two apparently unrelated clusters. This study highlights once more the importance of sharing genomic data to link apparently unrelated epidemiological clusters and to, remarkably, update molecular tests.