Listeria monocytogenes Transmission from Donated Blood to Platelet Transfusion Recipient, Italy.

We report Listeria monocytogenes infection in a patient in Italy who was transfused with pooled platelet concentrate. Genomic analysis revealed that L. monocytogenes isolates from the donor blood unit, the transfused platelets, and the patient's blood culture were genetically closely related, confirming transfusion transmission. Additional surveillance and secondary bacterial screening could improve transfusion safety.

Klebsiella pneumoniae carrying multiple alleles of antigen 43-encoding gene of Escherichia coli associated with biofilm formation.

A clinical strain of Klebsiella pneumoniae typed as sequence type 307 carrying three different alleles of the flu gene encoding the Escherichia coli virulence factor antigen 43 associated with biofilm formation was detected and characterized. The flu alleles are located in the chromosome inside putative integrative conjugative elements. The strain displays the phenotypes associated with Ag43, i.e. bi-phasic colony morphology and enhanced biofilm production. Furthermore, the strain produces low amount of capsule known to affect Ag43 function. Analysis of 1431 worldwide deposited genomes revealed that 3.7% Klebsiella pneumoniae carry one or two flu alleles.

First Italian report of a liver abscess and metastatic endogenous endophthalmitis caused by ST-23 hypervirulent Klebsiella pneumoniae in an immunocompetent individual.

BACKGROUND: Klebsiella pneumoniae is a common species in the gut of mammals and is widely distributed in the environment. However, the environmental source of hvKp that precedes gut colonization is unclear, but once that it reaches the gut there is a possible generalized spread y fecal-oral transmission especially in endemic areas. Liver abscess might develop when the bacteria, using its virulence factors, cross the intestinal barrier and invade the liver by the portal circulation. This syndrome, prevalent mostly in Asian countries, is increasingly reported in Western Countries and leaves open questions about the source of infection. CASE: Here we describe for the first time in Italy, a case of pyogenic liver abscess caused by a hypervirulent Klebsiella pneumoniae (HvKp) complicated by endophthalmitis and other metastatic infections in lung and prostate in an immunocompetent Chinese healthy individual with no recent travel in Asia. CONCLUSION: This case underlines the need for increased awareness of hypervirulent K. pneumoniae, even in settings where it occurs infrequently and where there are not evident epidemiological links.

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.

Reduction trend of mcr-1 circulation in Emilia-Romagna Region, Italy.

This study aims to describe trends of mcr-positive Enterobacterales in humans based on laboratory surveillance with a defined catchment population. The data source is the Micro-RER surveillance system, established in Emilia-Romagna region (Italy), to monitor the trend of mcr resistance. Enterobacterales isolates from human clinical samples with minimum inhibitory concentration (MIC) ≥ 2 mg/L for colistin were sent to the study reference laboratory for the detection of mcr genes. Isolates prospectively collected in the period 2018-2020 were considered for the assessment of population rates and trends; further analyses were carried out for the evaluation of clonality and horizontal mcr gene transfer. Previous isolates from local laboratory collection were also described. In the period 2018-2020, 1164 isolates were sent to the reference laboratory, and 51 (4.4%) were confirmed as mcr-positive: 50 mcr-1 (42 Escherichia coli, 6 Klebsiella pneumoniae, 2 Salmonella enterica) and 1 mcr-4 (Enterobacter cloacae). The number of mcr-positive isolates dropped from 24 in the first half of 2018 to 3 in the whole of 2020 (trend p value < 0.001). Genomic analyses showed the predominant role of the horizontal transfer of mcr genes through plasmids or dissemination of transposable elements compared to clonal dissemination of mcr-positive microorganisms. The study results demonstrate a substantial decrease in the circulation of mcr-1 plasmid genes in Emilia-Romagna Region.

A persistently replicating SARS-CoV-2 variant derived from an asymptomatic individual.

BACKGROUND: Since the first outbreak of SARS-CoV-2, the clinical characteristics of the Coronavirus Disease 2019 (COVID-19) have been progressively changed. Data reporting a viral intra-host and inter-host evolution favouring the appearance of mild SARS-CoV-2 strains are since being accumulating. To better understand the evolution of SARS-CoV-2 pathogenicity and its adaptation to the host, it is therefore crucial to investigate the genetic and phenotypic characteristics of SARS-CoV-2 strains circulating lately in the epidemic. METHODS: Nasopharyngeal swabs have been analyzed for viral load in the early (March 2020) and late (May 2020) phases of epidemic in Brescia, Italy. Isolation of SARS-CoV-2 from 2 high viral load specimens identified on March 9 (AP66) and on May 8 (GZ69) was performed on Vero E6 cells. Amount of virus released was assessed by quantitative PCR. Genotypic characterization of AP66 and GZ69 was performed by next generation sequencing followed by an in-depth in silico analysis of nucleotide mutations. RESULTS: The SARS-CoV-2 GZ69 strain, isolated in May from an asymptomatic healthcare worker, showed an unprecedented capability of replication in Vero E6 cells in the absence of any evident cytopathic effect. Vero E6 subculturing, up to passage 4, showed that SARS-CoV-2 GZ69 infection was as productive as the one sustained by the cytopathic strain AP66. Whole genome sequencing of the persistently replicating SARS-CoV-2 GZ69 has shown that this strain differs from the early AP66 variant in 9 nucleotide positions (C2939T; C3828T; G21784T; T21846C; T24631C; G28881A; G28882A; G28883C; G29810T) which lead to 6 non-synonymous substitutions spanning on ORF1ab (P892S; S1188L), S (K74N; I95T) and N (R203K, G204R) proteins. CONCLUSIONS: Identification of the peculiar SARS-CoV-2 GZ69 strain in the late Italian epidemic highlights the need to better characterize viral variants circulating among asymptomatic or paucisymptomatic individuals. The current approach could unravel the ways for future studies aimed at analyzing the selection process which favours viral mutations in the human host.

In vivo evolution of resistant subpopulations of KPC-producing Klebsiella pneumoniae during ceftazidime/avibactam treatment.

Objectives: KPC-producing Klebsiella pneumoniae (KPC-Kp) represent a serious problem worldwide. Herein, we describe the evolution of ceftazidime/avibactam resistance by sequencing longitudinal clinical isolates from a patient with KPC-Kp bloodstream infection undergoing ceftazidime/avibactam treatment. Methods: WGS was performed on one ceftazidime/avibactam-susceptible KPC-Kp (BOT-CA-S) and two phenotypically different ceftazidime/avibactam-resistant KPC-Kp with low (BOT-CA-R) and high (BOT-EMO) carbapenem MICs. The population diversity was assessed by the frequency of allele mutations and population analysis profiles (PAPs). Results: Phylogenetic analysis demonstrated clonal relatedness of the KPC-Kp isolates, all belonging to the clone ST1519. The D179Y mutation in blaKPC-3 was detected in both of the ceftazidime/avibactam-resistant KPC-Kp, whereas it was absent in the ceftazidime/avibactam-susceptible isolate. The mutation emerged independently in the two ceftazidime/avibactam-resistant isolates and was associated with a significant reduction in carbapenem MICs in BOT-CA-R, but not in BOT-EMO. WGS analysis revealed that the frequency of the D179Y mutation was 96.32% and 51.05% in BOT-CA-R and BOT-EMO, respectively. PAP results demonstrated that carbapenem resistance in BOT-EMO was due to the coexistence of mixed subpopulations harbouring WT and mutated blaKPC-3. A bacterial subpopulation with high ceftazidime/avibactam resistance for BOT-EMO KPC-Kp showed low carbapenem MICs, whereas a subpopulation with high meropenem resistance had a low MIC of ceftazidime/avibactam. Conclusions: Our analysis indicates that mixed subpopulations of ceftazidime/avibactam-resistant KPC-Kp emerge after ceftazidime/avibactam treatment. The evolution of different subpopulations that are highly resistant to ceftazidime/avibactam likely contributes to treatment failure, thereby highlighting the need for combination treatment strategies to limit selection of ceftazidime/avibactam-resistant KPC-Kp subpopulations.

Rise and fall of outbreak-specific clone inside endemic pulsotype of Salmonella 4,[5],12:i:-; insights from high-resolution molecular surveillance in Emilia-Romagna, Italy, 2012 to 2015.

Background and aimEpidemiology of human non-typhoid salmonellosis is characterised by recurrent emergence of new clones of the pathogen over time. Some clonal lines of Salmonella have shaped epidemiology of the disease at global level, as happened for serotype Enteritidis or, more recently, for Salmonella 4,[5],12:i:-, a monophasic variant of serotype Typhimurium. The same clonal behaviour is recognisable at sub-serotype level where single outbreaks or more generalised epidemics are attributable to defined clones. The aim of this study was to understand the dynamics of a clone of Salmonella 4,[5],12:i:- over a 3-year period (2012-15) in a province of Northern Italy where the clone caused a large outbreak in 2013. Furthermore, the role of candidate outbreak sources was investigated and the accuracy of multilocus variable-number tandem repeat analysis (MLVA) was evaluated. Methods: we retrospectively investigated the outbreak through whole genome sequencing (WGS) and further monitored the outbreak clone for 2 years after its conclusion. Results: The study showed the transient nature of the clone in the population, possibly as a consequence of its occasional expansion in a food-processing facility. We demonstrated that important weaknesses characterise conventional typing methods applied to clonal pathogens such as Salmonella 4,[5],12:i:-, namely lack of accuracy for MLVA and inadequate resolution power for PFGE to be reliably used for clone tracking. Conclusions: The study provided evidence for the remarkable prevention potential of whole genome sequencing used as a routine tool in systems that integrate human, food and animal surveillance.

Salmonella Detection and Counting on Pig Carcasses and Cutting Lines in Italian Slaughterhouses.

During 2014-2015, 300 pig carcasses before chilling and 85 food contact surfaces (FCSs) at cutting lines were tested for Salmonella in three slaughterhouses (namely A, B, and C) of northern Italy. In slaughterhouses A and B, four carcass sites of 100 cm2 each (from both the exterior and interior side) were swabbed with a single sponge. In abattoir C, four 100 cm2 sites of the exterior and the interior sides were swabbed with two independent sponges. The population average prevalence of Salmonella-positive carcasses (which takes into account the structure of the study design, with multiple samples collected in a single day) in slaughterhouses A and B was 12.3%, while in slaughterhouse C it was 11.2%. Presence of Salmonella on exterior and interior sides of carcasses showed a low level of concordance (only 3/12 of the contaminated carcasses were positive on both sides). No significant difference was found for FCSs contamination in the three slaughterhouses, with a population average prevalence of Salmonella-positive FCSs of 19.9%. In addition, we found that the clustering due to the day of sampling account for more than 36% and 60% of the overall prevalence variation on carcasses and FCSs, respectively. Eight serovars were identified, with Salmonella Derby as the most common type. The counting of Salmonella on carcasses showed large variability. It was low (<0.0075 most probable number [MPN]/cm2) in 46.6% of the carcasses and as high as 2.7 MPN/cm2 in 4.7%. Specifically, we found that counts on carcasses fit with "heavy tailed" distributions (lognormal and Weibull with a small shape parameter), suggesting not negligible probability of episodes of high Salmonella contamination. The mean values of contamination obtained with the two distributions ranged from 0.235 to 0.435 MPN/cm2.

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.

Characterization of an IncL/M plasmid carrying blaOXA-48 in a Klebsiella pneumoniae strain from Italy.

Here we report the complete nucleotide sequence of a 49.257-bp IncL/M conjugative plasmid (pRAY) carrying the blaOXA-48 gene collected from a Klebsiella pneumoniae clinical strain isolated in Italy. The genetic environment of pRAY plasmid revealed that the blaOXA-48 gene was located within a Tn1999.2 transposon. The pRAY plasmid differed from blaOXA-48-harboring IncL/M plasmids by genetic context and size. Comparative analysis demonstrated that pRAY plasmid lacked a region of ~15 kb carrying genes encoding proteins involved in pilus assembly and plasmid conjugative apparatus.

Processing-Dependent and Clonal Contamination Patterns of Listeria monocytogenes in the Cured Ham Food Chain Revealed by Genetic Analysis.

The quantitative and qualitative patterns of environmental contamination by Listeria monocytogenes were investigated in the production chain of dry-cured Parma ham. Standard arrays of surfaces were sampled in processing facilities during a single visit per plant in the three compartments of the food chain, i.e., ham production (19 plants) and postproduction, which was divided into deboning (43 plants) and slicing (25 plants) steps. The numbers of sampled surfaces were 384 in ham production, with 25 positive for L. monocytogenes, and 1,084 in postproduction, with 83 positives. Statistical analysis of the prevalence of contaminated surfaces showed that in ham production, contamination was higher at the beginning of processing and declined significantly toward the end, while in postproduction, prevalence rose toward the end of processing. Prevalence was higher in the deboning facilities than in slicing facilities and was dependent on the type of surface (floor/drainage > clothing > equipment). The qualitative pattern of contamination was investigated through an analysis of the survey isolates and a set of isolates derived from routine monitoring, including longitudinal isolations. Pulsed-field gel electrophoresis (PFGE) and whole-genome single-nucleotide polymorphism (SNP) analysis revealed a remarkable clonality of L. monocytogenes within plants, with the detection of 16 plant-specific clones out of 17 establishments with multiple isolates. Repeated detections of clonal isolates >6 months apart were also observed. Six was the maximum number of between-isolate differences in core SNPs observed within these clones. Based on the same six-SNP threshold, three clusters of clonal isolates, shared by six establishments, were also identified. The spread of L. monocytogenes within and between plants, as indicated by its clonal behavior, is a matter of concern for the hygienic management of establishments.

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.

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.

The DNA binding mechanism of a SSB protein from Lactococcus lactis siphophage p2.

Virulent lactococcal phages of the Siphoviridae family are responsible for the industrial milk fermentation failures worldwide. Lactococcus lactis, a Gram-positive bacterium widely used for the manufacture of fermented dairy products, is subjected to infections by virulent phages, predominantly those of the 936 group, including phage p2. Among the proteins coded by lactococcal phage genomes, of special interest are those expressed early, which are crucial to efficiently carry out the phage lytic cycle. We previously identified and solved the 3D structure of lactococcal phage p2 ORF34, a single stranded DNA binding protein (SSBp2). Here we investigated the molecular basis of ORF34 binding mechanism to DNA. DNA docking on SSBp2 and Molecular Dynamics simulations of the resulting complex identified R15 as a crucial residue for ssDNA binding. Electrophoretic Mobility Shift Assays (EMSA) and Atomic Force Microscopy (AFM) imaging revealed the inability of the Arg15Ala mutant to bind ssDNA, as compared to the native protein. Since R15 is highly conserved among lactococcal SSBs, we propose that its role in the SSBp2/DNA complex stabilization might be extended to all the members of this protein family.

The porcine corpus luteum as a model for studying the effects of nanoplastics.

Nanoplastics (NPs) affect fertility. We evaluated the effects of NPs treatment on luteal and endothelial cells. We examined crucial markers of growth and redox status. NPs treatment did not induce changes in ATP levels in luteal cells, while it increased (p< 0.05) their proliferation. In endothelial cells, no change in proliferation was detected, while an increase (p<0.05) in ATP levels was observed. The increase of reactive oxygen species, superoxide anion (p<0.05) and nitric oxide (p<0.001) was detected in both cell types, which also showed changes in superoxide dismutase enzyme activity as well as an increase of non-enzymatic antioxidant power (p<0.05). A decrease (p<0.05) in progesterone production as well as an increase of vascular endothelial growth factor A levels were detected (p<0.05). In addition, a dose-dependent accumulation of NPs in endothelial cells was shown, that likely occurred through adhesion and internalization. Results underline potential risk of NPs for corpus luteum functionality.

Development of inhalation powders containing lactic acid bacteria with antimicrobial activity against Pseudomonas aeruginosa.

OBJECTIVES: The aim of the project was to develop and characterise powders containing a probiotic (Lactiplantibacillus plantarum [Lpb. plantarum], Lacticaseibacillus rhamnosus, or Lactobacillus acidophilus) to be administered to the lung for the containment of pathogen growth in patients with lung infections. METHODS: The optimised spray drying process for the powder manufacturing was able to preserve viability of the bacteria, which decreased of only one log unit and was maintained up to 30 days. RESULTS: Probiotic powders showed a high respirability (42%-50% of particles had a size < 5 µm) suitable for lung deposition and were proven safe on A549 and Calu-3 cells up to a concentration of 107 colony-forming units/mL. The Lpb. plantarum adhesion to both cell lines tested was at least 10%. Surprisingly, Lpb. plantarum powder was bactericidal at a concentration of 106 colony-forming units/mL on P. aeruginosa, whereas the other two strains were bacteriostatic. CONCLUSION: This work represents a promising starting point to consider a probiotic inhalation powder a value in keeping the growth of pathogenic microflora in check during the antibiotic inhalation therapy suspension in cystic fibrosis treatment regimen. This approach could also be advantageous for interfering competitively with pathogenic bacteria and promoting the restoration of the healthy microbiota.

Characterization of the genome of the dairy Lactobacillus helveticus bacteriophage {Phi}AQ113.

The complete genomic sequence of the dairy Lactobacillus helveticus bacteriophage ΦAQ113 was determined. Phage ΦAQ113 is a Myoviridae bacteriophage with an isometric capsid and a contractile tail. The final assembled consensus sequence revealed a linear, circularly permuted, double-stranded DNA genome with a size of 36,566 bp and a G+C content of 37%. Fifty-six open reading frames (ORFs) were predicted, and a putative function was assigned to approximately 90% of them. The ΦAQ113 genome shows functionally related genes clustered together in a genome structure composed of modules for DNA replication/regulation, DNA packaging, head and tail morphogenesis, cell lysis, and lysogeny. The identification of genes involved in the establishment of lysogeny indicates that it may have originated as a temperate phage, even if it was isolated from natural cheese whey starters as a virulent phage, because it is able to propagate in a sensitive host strain. Additionally, we discovered that the ΦAQ113 phage genome is closely related to Lactobacillus gasseri phage KC5a and Lactobacillus johnsonii phage Lj771 genomes. The phylogenetic similarities between L. helveticus phage ΦAQ113 and two phages that belong to gut species confirm a possible common ancestral origin and support the increasing consideration of L. helveticus as a health-promoting organism.