Genomic analysis of 40 prophages located in the genomes of 16 carbapenemase-producing clinical strains of Klebsiella pneumoniae

Complete genome sequence, phenotypic correlation and pangenome analysis of uropathogenic Klebsiella spp

Urinary tract infections (UTI) by antibiotic resistant and virulent K. pneumoniae are a growing concern. Understanding the genome and validating the genomic profile along with pangenome analysis will facilitate surveillance of high-risk clones of K. pneumoniae to underpin management strategies toward early detection. The present study aims to correlate resistome with phenotypic antimicrobial resistance and virulome with pathogenicity in Klebsiella spp. The present study aimed to perform complete genome sequences of Klebsiella spp. and to analyse the correlation of resistome with phenotypic antimicrobial resistance and virulome with pathogenicity. To understand the resistome, pangenome and virulome in the Klebsiella spp, the ResFinder, CARD, IS Finder, PlasmidFinder, PHASTER, Roary, VFDB were used. The phenotypic susceptibility profiling identified the uropathogenic kp3 to exhibit multi drug resistance. The resistome and in vitro antimicrobial profiling showed concordance with all the tested antibiotics against the study strains. Hypermucoviscosity was not observed for any of the test isolates; this phenotypic character matches perfectly with the absence of rmpA and magA genes. To the best of our knowledge, this is the first report on the presence of ste, stf, stc and sti major fimbrial operons of Salmonella enterica serotype Typhimurium in K. pneumoniae genome. The study identifies the discordance of virulome and virulence in Klebsiella spp. The complete genome analysis and phenotypic correlation identify uropathogenic K. pneumoniae kp3 as a carbapenem-resistant and virulent pathogen. The Pangenome of K. pneumoniae was open suggesting high genetic diversity. Diverse K serotypes were observed. Sequence typing reveals the prevalence of K. pneumoniae high-risk clones in UTI catheterised patients. The study also highlights the concordance of resistome and in vitro susceptibility tests. Importantly, the study identifies the necessity of virulome and phenotypic virulence markers for timely diagnosis and immediate treatment for the management of high-risk K. pneumoniae clones.

A universe of human gut-derived bacterial prophages: unveiling the hidden viral players in intestinal microecology

Prophages, which are an existing form of temperate phages that integrate into host genomes, have been found extensively present in diverse bacterial species. The human gut microbiome, characterized by its complexity, dynamism, and interconnectivity among multiple species, remains inadequately understood in terms of the global landscape of bacterial prophages. Here, we analyzed 43,942 human gut-derived bacterial genomes (439 species of 12 phyla) and identified 105,613 prophage regions in ~ 92% of them. 16254 complete prophages were distributed in ~ 24% of bacteria, indicating an extremely uneven prophage distribution across various species within the human gut. Among all identified prophages, ~4% possessed cross-genera (2-20 genera) integration capacity, while ~ 17% displayed broad infection host ranges (targeting 2-35 genera). Functional gene annotation revealed that antibiotic-resistance genes and toxin-related genes were detected in ~ 2.5% and ~ 5.8% of all prophages, respectively. Furthermore, through sequence alignments between our obtained prophages and publicly available human gut phageome contigs, we have observed that ~ 72% of non-redundant prophages are previously unreported genomes, and they illuminate ~ 6.5-9.5% of the individual intestinal "viral dark matter". Our study represents the first comprehensive depiction of human gut-derived prophages, provides a substantial collection of reference sequences for forthcoming human gut phageome-related investigations, and potentially enables better risk assessment of prophage dissemination.

Prophage identification and molecular analysis in the genomes of Pseudomonas aeruginosa strains isolated from critical care patients

Prophages are bacteriophages integrated into the bacterial host's chromosome. This research aims to analyze and characterize the existing prophages within a collection of 53 Pseudomonas aeruginosa strains from intensive care units (ICUs) in Portugal and Spain. A total of 113 prophages were localized in the collection, with 18 of them being present in more than one strain simultaneously. After annotation, five of them were discarded as incomplete, and the 13 remaining prophages were characterized. Of 13, 10 belonged to the siphovirus tail morphology group, 2 to the podovirus tail morphology group, and 1 to the myovirus tail morphology group. All prophages had a length ranging from 20,199 to 63,401 bp and a GC% between 56.2% and 63.6%. The number of open reading frames (ORFs) oscillated between 32 and 88, and in 3/13 prophages, more than 50% of the ORFs had an unknown function. With our findings, we show that prophages are present in the majority of the P. aeruginosa strains isolated from Portuguese and Spanish critically ill patients, many of them found in more than one circulating strain at the same time and following a similar clonal distribution pattern. Although a great sum of ORFs had an unknown function, number of proteins in relation to viral defense (anti-CRISPR proteins, toxin/antitoxin modules, proteins against restriction-modification systems) as well as to prophage interference into their host's quorum sensing system and regulatory cascades were found. This supports the idea that prophages have an influence in bacterial pathogenesis and anti-phage defense. IMPORTANCE Despite being known for decades, prophages remain understudied when compared to the lytic phages employed in phage therapy. This research aims to shed some light into the nature, composition, and role of prophages found within a set of circulating strains of Pseudomas aeruginosa, with special attention to high-risk clones. Given the fact that prophages can effectively influence bacterial pathogenesis, prophage basic research constitutes a topic of growing interest. Furthermore, the abundance of viral defense and regulatory proteins within prophage genomes detected in this study evidences the importance of characterizing the most frequent prophages in circulating clinical strains and in high-risk clones if phage therapy is to be used.

Divergent endophytic viromes and phage genome repertoires among banana (Musa) species

Introduction

Viruses generally cause disease, but some viruses may be beneficial as resident regulators of their hosts or host microbiomes. Plant-associated viruses can help plants survive by increasing stress tolerance or regulating endophytic communities. The goal of this study was to characterize endophytic virus communities in banana and plantain (Musa spp.) genotypes, including cultivated and wild species, to assess virome repertoires and detect novel viruses.

Methods

DNA viral communities were characterized by shotgun sequencing of an enriched endosphere extract from leaves and roots or corm of 7 distinct Musa genotypes (M. balbisiana, Thai Black, M. textilis, M. sikkimensis, Dwarf Cavendish, Williams Hybrid, and FHIA-25 Hybrid).

Results

Results showed abundant virus-like contigs up to 108,191 bp long with higher relative abundance in leaves than roots. Analyses predicted 733 phage species in 51 families, with little overlap in phage communities among plants. Phage diversity was higher in roots and in diploid wild hosts. Ackermanniviridae and Rhizobium phage were generally the most abundant taxa. A Rhizobium RR1-like phage related to a phage of an endophytic tumor-causing rhizobium was found, bearing a holin gene and a partial Shiga-like toxin gene, raising interest in its potential to regulate endophytic Rhizobiaceae. Klebsiella phages were of interest for possible protection against Fusarium wilt, and other phages were predicted with potential to regulate Erwinia, Pectobacterium, and Ralstonia-associated diseases. Although abundant phage-containing contigs were functionally annotated, revealing 1,038 predicted viral protein domains, gene repertoires showed high divergence from database sequences, suggesting novel phages in these banana cultivars. Plant DNA viruses included 56 species of Badnavirus and 26 additional non-Musa plant viruses with distributions that suggested a mixture of resident and transient plant DNA viruses in these samples.

Discussion

Together, the disparate viral communities in these plants from a shared environment suggest hosts drive the composition of these virus communities. This study forms a first step in understanding the endophytic virome in this globally important food crop, which is currently threatened by fungal, bacterial, and viral diseases.

Case report: Analysis of phage therapy failure in a patient with a Pseudomonas aeruginosa prosthetic vascular graft infection

Clinical case of a patient with a Pseudomonas aeruginosa multidrug-resistant prosthetic vascular graft infection which was treated with a cocktail of phages (PT07, 14/01, and PNM) in combination with ceftazidime-avibactam (CZA). After the application of the phage treatment and in absence of antimicrobial therapy, a new P. aeruginosa bloodstream infection (BSI) with a septic residual limb metastasis occurred, now involving a wild-type strain being susceptible to ß-lactams and quinolones. Clinical strains were analyzed by microbiology and whole genome sequencing techniques. In relation with phage administration, the clinical isolates of P. aeruginosa before phage therapy (HE2011471) and post phage therapy (HE2105886) showed a clonal relationship but with important genomic changes which could be involved in the resistance to this therapy. Finally, phenotypic studies showed a decrease in Minimum Inhibitory Concentration (MIC) to ß-lactams and quinolones as well as an increase of the biofilm production and phage resistant mutants in the clinical isolate of P. aeruginosa post phage therapy.

NDM-5-carrying Klebsiella pneumoniae ST437 belonging to high-risk clonal complex (CC11) from an urban river in eastern India

In this study, we described the carbapenem bla NDM-5-carrying extensive drug-resistant (XDR) K. pneumoniae ST437 from an urban river water Kathajodi in Odisha, India. The presence of carbapenem and co-occurrence of other resistance determinants (bla NDM-5, bla CTX-M, bla SHV, and bla TEM), virulence factors (fimH, mrkD, entB, irp-1, and ybtS), and capsular serotype (K54) represent its pathogenic potential. The insertion sequence ISAba125 and the bleomycin resistance gene ble MBL at upstream and downstream, respectively, could play a significant role in the horizontal transmission of the bla NDM-5. Its biofilm formation ability contributes toward environmental protection and its survivability. MLST analysis assigned the isolate to ST437 and clonal lineage to ST11 (CC11) with a single locus variant. The ST437 K. pneumoniae, a global epidemic clone, has been reported in North America, Europe, and Asia. This work contributes in understanding of the mechanisms behind the spread of bla NDM-5 K. pneumoniae ST437 and demands extensive molecular surveillance of river and nearby hospitals for better community health.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03556-5.

Proteomic Study of the Interactions between Phages and the Bacterial Host Klebsiella pneumoniae

Phages and bacteria have acquired resistance mechanisms for protection. In this context, the aims of the present study were to analyze the proteins isolated from 21 novel lytic phages of Klebsiella pneumoniae in search of defense mechanisms against bacteria and also to determine the infective capacity of the phages. A proteomic study was also conducted to investigate the defense mechanisms of two clinical isolates of K. pneumoniae infected by phages. For this purpose, the 21 lytic phages were sequenced and de novo assembled. The host range was determined in a collection of 47 clinical isolates of K. pneumoniae, revealing the variable infective capacity of the phages. Genome sequencing showed that all of the phages were lytic phages belonging to the order Caudovirales. Phage sequence analysis revealed that the proteins were organized in functional modules within the genome. Although most of the proteins have unknown functions, multiple proteins were associated with defense mechanisms against bacteria, including the restriction-modification system, the toxin-antitoxin system, evasion of DNA degradation, blocking of host restriction and modification, the orphan CRISPR-Cas system, and the anti-CRISPR system. Proteomic study of the phage-host interactions (i.e., between isolates K3574 and K3320, which have intact CRISPR-Cas systems, and phages vB_KpnS-VAC35 and vB_KpnM-VAC36, respectively) revealed the presence of several defense mechanisms against phage infection (prophage, defense/virulence/resistance, oxidative stress and plasmid proteins) in the bacteria, and of the Acr candidate (anti-CRISPR protein) in the phages. IMPORTANCE Researchers, including microbiologists and infectious disease specialists, require more knowledge about the interactions between phages and their bacterial hosts and about their defense mechanisms. In this study, we analyzed the molecular mechanisms of viral and bacterial defense in phages infecting clinical isolates of K. pneumoniae. Viral defense mechanisms included restriction-modification system evasion, the toxin-antitoxin (TA) system, DNA degradation evasion, blocking of host restriction and modification, and resistance to the abortive infection system, anti-CRISPR and CRISPR-Cas systems. Regarding bacterial defense mechanisms, proteomic analysis revealed expression of proteins involved in the prophage (FtsH protease modulator), plasmid (cupin phosphomannose isomerase protein), defense/virulence/resistance (porins, efflux pumps, lipopolysaccharide, pilus elements, quorum network proteins, TA systems, and methyltransferases), oxidative stress mechanisms, and Acr candidates (anti-CRISPR protein). The findings reveal some important molecular mechanisms involved in the phage-host bacterial interactions; however, further study in this field is required to improve the efficacy of phage therapy.

Global diversity and distribution of prophages are lineage-specific within the Ralstonia solanacearum species complex

Background

Ralstonia solanacearum species complex (RSSC) strains are destructive plant pathogenic bacteria and the causative agents of bacterial wilt disease, infecting over 200 plant species worldwide. In addition to chromosomal genes, their virulence is mediated by mobile genetic elements including integrated DNA of bacteriophages, i.e., prophages, which may carry fitness-associated auxiliary genes or modulate host gene expression. Although experimental studies have characterised several prophages that shape RSSC virulence, the global diversity, distribution, and wider functional gene content of RSSC prophages are unknown. In this study, prophages were identified in a diverse collection of 192 RSSC draft genome assemblies originating from six continents.

Results

Prophages were identified bioinformatically and their diversity investigated using genetic distance measures, gene content, GC, and total length. Prophage distributions were characterised using metadata on RSSC strain geographic origin and lineage classification (phylotypes), and their functional gene content was assessed by identifying putative prophage-encoded auxiliary genes. In total, 313 intact prophages were identified, forming ten genetically distinct clusters. These included six prophage clusters with similarity to the Inoviridae, Myoviridae, and Siphoviridae phage families, and four uncharacterised clusters, possibly representing novel, previously undescribed phages. The prophages had broad geographical distributions, being present across multiple continents. However, they were generally host phylogenetic lineage-specific, and overall, prophage diversity was proportional to the genetic diversity of their hosts. The prophages contained many auxiliary genes involved in metabolism and virulence of both phage and bacteria.

Conclusions

Our results show that while RSSC prophages are highly diverse globally, they make lineage-specific contributions to the RSSC accessory genome, which could have resulted from shared coevolutionary history.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08909-7.

Analysis of intact prophages in genomes of Paenibacillus larvae: An important pathogen for bees

Paenibacillus larvae is the etiological agent of American Foulbrood (AFB), a highly contagious and worldwide spread bacterial disease that affects honeybee brood. In this study, all complete P. larvae genomes available on the NCBI database were analyzed in order to detect presence of prophages using the PHASTER software. A total of 55 intact prophages were identified in 11 P. larvae genomes (5.0 ± 2.3 per genome) and were further investigated for the presence of genes encoding relevant traits related to P. larvae. A closer look at the prophage genomes revealed the presence of several putative genes such as metabolic and antimicrobial resistance genes, toxins or bacteriocins, potentially influencing host performance. Some of the coding DNA sequences (CDS) were present in all ERIC-genotypes, while others were only found in a specific genotype. While CDS encoding toxins and antitoxins such as HicB and MazE were found in prophages of all bacterial genotypes, others, from the same category, were provided by prophages particularly to ERIC I (enhancin-like toxin), ERIC II (antitoxin SocA) and ERIC V strains (subunit of Panton-Valentine leukocidin system (PVL) LukF-PV). This is the first in-depth analysis of P. larvae prophages. It provides better knowledge on their impact in the evolution of virulence and fitness of P. larvae, by discovering new features assigned by the viruses.

The role of PemIK (PemK/PemI) type II TA system from Klebsiella pneumoniae clinical strains in lytic phage infection

Since their discovery, toxin-antitoxin (TA) systems have captivated the attention of many scientists. Recent studies have demonstrated that TA systems play a key role in phage inhibition. The aim of the present study was to investigate the role of the PemIK (PemK/PemI) type II TA system in phage inhibition by its intrinsic expression in clinical strains of Klebsiella pneumoniae carrying the lncL plasmid, which harbours the carbapenemase OXA-48 and the PemK/PemI TA system. Furthermore, induced expression of the system in an IPTG-inducible plasmid in a reference strain of K. pneumoniae ATCC10031 was also studied. The results showed that induced expression of the whole TA system did not inhibit phage infection, whereas overexpression of the pemK toxin prevented early infection. To investigate the molecular mechanism involved in the PemK toxin-mediated inhibition of phage infection, assays measuring metabolic activity and viability were performed, revealing that overexpression of the PemK toxin led to dormancy of the bacteria. Thus, we demonstrate that the PemK/PemI TA system plays a role in phage infection and that the action of the free toxin induces a dormant state in the cells, resulting in inhibition of phage infections.

Comparative genomics of Chinese and international isolates of Escherichia albertii: population structure and evolution of virulence and antimicrobial resistance

Escherichia albertii is a recently recognized species in the genus Escherichia that causes diarrhoea. The population structure, genetic diversity and genomic features have not been fully examined. Here, 169 E. albertii isolates from different sources and regions in China were sequenced and combined with 312 publicly available genomes (from additional 14 countries) for genomic analyses. The E. albertii population was divided into two clades and eight lineages, with lineage 3 (L3), L5 and L8 more common in China. Clinical isolates were observed in all clades/lineages. Virulence genes were found to be distributed differently among lineages: subtypes of the intimin encoding gene eae and the cytolethal distending toxin gene cdtB were lineage associated, and the second type three secretion system (ETT2) island was truncated in L3 and L6. Seven new eae subtypes and one new cdtB subtype (cdtB-VI) were identified. Alarmingly, 85.9 % of the Chinese E. albertii isolates were predicted to be multidrug-resistant (MDR) with 35.9 % harbouring genes capable of conferring resistance to 10 to 14 different drug classes. The majority of the MDR isolates were of poultry source from China and belonged to four sequence types (STs) [ST4638, ST4479, ST4633 and ST4488]. Thirty-four plasmids with some carrying MDR and virulence genes, and 130 prophages were identified from 17 complete E. albertii genomes. The 130 intact prophages were clustered into five groups, with group five prophages harbouring more virulence genes. We further identified three E. albertii specific genes as markers for the identification of this species. Our findings provided fundamental insights into the population structure, virulence variation and drug resistance of E. albertii.

Genomic Analysis of Prophages from Klebsiella pneumoniae Clinical Isolates

Klebsiella pneumoniae is an increasing threat to public health and represents one of the most concerning pathogens involved in life-threatening infections. The resistant and virulence determinants are coded by mobile genetic elements which can easily spread between bacteria populations and co-evolve with its genomic host. In this study, we present the full genomic sequences, insertion sites and phylogenetic analysis of 150 prophages found in 40 K. pneumoniae clinical isolates obtained from an outbreak in a Portuguese hospital. All strains harbored at least one prophage and we identified 104 intact prophages (69.3%). The prophage size ranges from 29.7 to 50.6 kbp, coding between 32 and 78 putative genes. The prophage GC content is 51.2%, lower than the average GC content of 57.1% in K. pneumoniae. Complete prophages were classified into three families in the order Caudolovirales: Myoviridae (59.6%), Siphoviridae (38.5%) and Podoviridae (1.9%). In addition, an alignment and phylogenetic analysis revealed nine distinct clusters. Evidence of recombination was detected within the genome of some prophages but, in most cases, proteins involved in viral structure, transcription, replication and regulation (lysogenic/lysis) were maintained. These results support the knowledge that prophages are diverse and widely disseminated in K. pneumoniae genomes, contributing to the evolution of this species and conferring additional phenotypes. Moreover, we identified K. pneumoniae prophages in a set of endolysin genes, which were found to code for proteins with lysozyme activity, cleaving the β-1,4 linkages between N-acetylmuramic acid and N-acetyl-D-glucosamine residues in the peptidoglycan network and thus representing genes with the potential for lysin phage therapy.

Comprehensive Scanning of Prophages in Lactobacillus: Distribution, Diversity, Antibiotic Resistance Genes, and Linkages with CRISPR-Cas Systems

Prophage integration, release, and dissemination exert various effects on host bacteria. In the genus Lactobacillus, they may cause bacteriophage contamination during fermentation and even regulate bacterial populations in the gut. However, little is known about their distribution, genetic architecture, and relationships with their hosts. Here, we conducted prophage prediction analysis on 1,472 genomes from 16 different Lactobacillus species and found prophage fragments in almost all lactobacilli (99.8%), with 1,459 predicted intact prophages identified in 64.1% of the strains. We present an uneven prophage distribution among Lactobacillus species; multihabitat species retained more prophages in their genomes than restricted-habitat species. Characterization of the genome features, average nucleotide identity, and landscape visualization presented a high genome diversity of Lactobacillus prophages. We detected antibiotic resistance genes in more than 10% of Lactobacillus prophages and validated that the occurrence of resistance genes conferred by prophage integration was possibly associated with phenotypic resistance in Lactobacillus plantarum. Furthermore, our broad and comprehensive examination of the distribution of CRISPR-Cas systems across the genomes predicted type I and type III systems as potential antagonistic elements of Lactobacillus prophage.

IMPORTANCE Lactobacilli are inherent microorganisms in the human gut and are widely used in the food processing industries due to their probiotic properties. Prophages were reportedly hidden in numerous Lactobacillus genomes and can potentially contaminate entire batches of fermentation or modulate the intestinal microecology once they are released. Therefore, a comprehensive scanning of prophages in Lactobacillus is essential for the safety evaluation and application development of probiotic candidates. We show that prophages are widely distributed among lactobacilli; however, intact prophages are more common in multihabitat species and display wide variations in genome feature, integration site, and genomic organization. Our data of the prophage-mediated antibiotic resistance genes (ARGs) and the resistance phenotype of lactobacilli provide evidence for deciphering the putative role of prophages as vectors of the ARGs. Furthermore, understanding the association between prophages and CRISPR-Cas systems is crucial to appreciate the coevolution of phages and Lactobacillus.

KPC-53, a KPC-3 Variant of Clinical Origin Associated with Reduced Susceptibility to Ceftazidime-Avibactam

This study reports on the characterization of a Klebsiella pneumoniae clinical isolate showing high-level resistance to ceftazidime-avibactam associated with the production of KPC-53, a KPC-3 variant exhibiting a Leu167Glu168 duplication in the Ω-loop and a loss of carbapenemase activity. Whole-genome sequencing (WGS) revealed the presence of two copies of bla _KPC-53, located on a pKpQIL-like plasmid and on a plasmid prophage of the Siphoviridae family, respectively. The present findings provide new insights into the mechanisms of resistance to ceftazidime-avibactam.

Identification, Characterization, and Genomic Analysis of Novel Serratia Temperate Phages from a Gold Mine

Bacteria of the genus Serratia inhabit a variety of ecological niches like water, soil, and the bodies of animals, and have a wide range of lifestyles. Currently, the complete genome sequences of 25 Serratia phages are available in the NCBI database. All of them were isolated from nutrient-rich environments like sewage, with the use of clinical Serratia strains as hosts. In this study, we identified a novel Serratia myovirus named vB_SspM_BZS1. Both the phage and its host Serratia sp. OS31 were isolated from the same oligotrophic environment, namely, an abandoned gold mine (Zloty Stok, Poland). The BZS1 phage was thoroughly characterized here in terms of its genomics, morphology, and infection kinetics. We also demonstrated that Serratia sp. OS31 was lysogenized by mitomycin-inducible siphovirus vB_SspS_OS31. Comparative analyses revealed that vB_SspM_BZS1 and vB_SspS_OS31 were remote from the known Serratia phages. Moreover, vB_SspM_BZS1 was only distantly related to other viruses. However, we discovered similar prophage sequences in genomes of various bacteria here. Additionally, a protein-based similarity network showed a high diversity of Serratia phages in general, as they were scattered across nineteen different clusters. In summary, this work broadened our knowledge on the diverse relationships of Serratia phages.