Going viral: The role of mobile genetic elements in bacterial immunity.

Bacteriophages and other mobile genetic elements (MGEs) pose a significant threat to bacteria, subjecting them to constant attacks. In response, bacteria have evolved a sophisticated immune system that employs diverse defensive strategies and mechanisms. Remarkably, a growing body of evidence suggests that most of these defenses are encoded by MGEs themselves. This realization challenges our traditional understanding of bacterial immunity and raises intriguing questions about the evolutionary forces at play. Our review provides a comprehensive overview of the latest findings on the main families of MGEs and the defense systems they encode. We also highlight how a vast diversity of defense systems remains to be discovered and their mechanism of mobility understood. Altogether, the composition and distribution of defense systems in bacterial genomes only makes sense in the light of the ecological and evolutionary interactions of a complex network of MGEs.

Systematic interrogation of CRISPR antimicrobials in Klebsiella pneumoniae reveals nuclease-, guide- and strain-dependent features influencing antimicrobial activity.

CRISPR-Cas systems can be utilized as programmable-spectrum antimicrobials to combat bacterial infections. However, how CRISPR nucleases perform as antimicrobials across target sites and strains remains poorly explored. Here, we address this knowledge gap by systematically interrogating the use of CRISPR antimicrobials using multidrug-resistant and hypervirulent strains of Klebsiella pneumoniae as models. Comparing different Cas nucleases, DNA-targeting nucleases outperformed RNA-targeting nucleases based on the tested targets. Focusing on AsCas12a that exhibited robust targeting across different strains, we found that the elucidated modes of escape varied widely, restraining opportunities to enhance killing. We also encountered individual guide RNAs yielding different extents of clearance across strains, which were linked to an interplay between improper gRNA folding and strain-specific DNA repair and survival. To explore features that could improve targeting across strains, we performed a genome-wide screen in different K. pneumoniae strains that yielded guide design rules and trained an algorithm for predicting guide efficiency. Finally, we showed that Cas12a antimicrobials can be exploited to eliminate K. pneumoniae when encoded in phagemids delivered by T7-like phages. Altogether, our results highlight the importance of evaluating antimicrobial activity of CRISPR antimicrobials across relevant strains and define critical parameters for efficient CRISPR-based targeting.

Inference of the Life Cycle of Environmental Phages from Genomic Signature Distances to Their Hosts.

The environmental impact of uncultured phages is shaped by their preferred life cycle (lytic or lysogenic). However, our ability to predict it is very limited. We aimed to discriminate between lytic and lysogenic phages by comparing the similarity of their genomic signatures to those of their hosts, reflecting their co-evolution. We tested two approaches: (1) similarities of tetramer relative frequencies, (2) alignment-free comparisons based on exact k = 14 oligonucleotide matches. First, we explored 5126 reference bacterial host strains and 284 associated phages and found an approximate threshold for distinguishing lysogenic and lytic phages using both oligonucleotide-based methods. The analysis of 6482 plasmids revealed the potential for horizontal gene transfer between different host genera and, in some cases, distant bacterial taxa. Subsequently, we experimentally analyzed combinations of 138 Klebsiella pneumoniae strains and their 41 phages and found that the phages with the largest number of interactions with these strains in the laboratory had the shortest genomic distances to K. pneumoniae. We then applied our methods to 24 single-cells from a hot spring biofilm containing 41 uncultured phage-host pairs, and the results were compatible with the lysogenic life cycle of phages detected in this environment. In conclusion, oligonucleotide-based genome analysis methods can be used for predictions of (1) life cycles of environmental phages, (2) phages with the broadest host range in culture collections, and (3) potential horizontal gene transfer by plasmids.

Human Anelloviruses: Influence of Demographic Factors, Recombination, and Worldwide Diversity.

Anelloviruses represent the major and most diverse component of the healthy human virome, referred to as the anellome. In this study, we determined the anellome of 50 blood donors, forming two sex- and age-matched groups. Anelloviruses were detected in 86% of the donors. The number of detected anelloviruses increased with age and was approximately twice as high in men as in women. A total of 349 complete or nearly complete genomes were classified as belonging to torque teno virus (TTV), torque teno mini virus (TTMV), and torque teno midi virus (TTMDV) anellovirus genera (197, 88, and 64 sequences, respectively). Most donors had intergenus (69.8%) or intragenus (72.1%) coinfections. Despite the limited number of sequences, intradonor recombination analysis showed 6 intragenus recombination events in ORF1. As thousands of anellovirus sequences have been described recently, we finally analyzed the global diversity of human anelloviruses. Species richness and diversity were close to saturation in each anellovirus genus. Recombination was found to be the main factor promoting diversity, although its effect was significantly lower in TTV than in TTMV and TTMDV. Overall, our results suggest that differences in diversity between genera may be caused by variations in the relative contribution of recombination. IMPORTANCE Anelloviruses are the most common human infectious viruses and are considered essentially harmless. Compared to other human viruses, they are characterized by enormous diversity, and recombination is suggested to play an important role in their diversification and evolution. Here, by analyzing the composition of the plasma anellome of 50 blood donors, we find that recombination is also a determinant of viral evolution at the intradonor level. On a larger scale, analysis of anellovirus sequences currently available in databases shows that their diversity is close to saturation and differs among the three human anellovirus genera and that recombination is the main factor explaining this intergenus variability. Global characterization of anellovirus diversity could provide clues about possible associations between certain virus variants and pathologies, as well as facilitate the implementation of unbiased PCR-based detection protocols, which may be relevant for using anelloviruses as endogenous markers of immune status.

Genetic determinants of host tropism in Klebsiella phages.

Bacteriophages play key roles in bacterial ecology and evolution and are potential antimicrobials. However, the determinants of phage-host specificity remain elusive. Here, we isolate 46 phages to challenge 138 representative clinical isolates of Klebsiella pneumoniae, a widespread opportunistic pathogen. Spot tests show a narrow host range for most phages, with <2% of 6,319 phage-host combinations tested yielding detectable interactions. Bacterial capsule diversity is the main factor restricting phage host range. Consequently, phage-encoded depolymerases are key determinants of host tropism, and depolymerase sequence types are associated with the ability to infect specific capsular types across phage families. However, all phages with a broader host range found do not encode canonical depolymerases, suggesting alternative modes of entry. These findings expand our knowledge of the complex interactions between bacteria and their viruses and point out the feasibility of predicting the first steps of phage infection using bacterial and phage genome sequences.

Tracking the Emergence and Dissemination of a blaNDM-23 Gene in a Multidrug Resistance Plasmid of Klebsiella pneumoniae.

Since the discovery of blaNDM-1, NDM ß-lactamases have become one of the most widespread carbapenemases worldwide. To date, 43 different NDM variants have been reported but some, such as blaNDM-23, have not been characterized in detail yet. Here, we describe the emergence of a novel blaNDM-23 allele from a blaNDM-1 ancestor and the multidrug resistance plasmid that has disseminated it through a Klebsiella pneumoniae ST437 clone in several Spanish hospitals. Between 2016 and 2019, 1,972 isolates were collected in an epidemiological survey for extended-spectrum-ß-lactamase (ESBL)-producing Klebsiella pneumoniae in the Comunitat Valenciana (Spain). Three carbapenem-resistant strains failed to be detected by carbapenemase-producing Enterobacteriaceae (CPE) screening tests. These isolates carried a blaNDM-23 gene. To characterize this gene, its emergence, and its dissemination, we performed antimicrobial susceptibility tests, hybrid sequencing with Illumina and Nanopore technologies, and phylogenetic analyses. The MICs of the blaNDM-23 allele were identical to those of the blaNDM-1 allele. The blaNDM-23 allele was found in 14 isolates on a 97-kb nonmobilizable, multidrug-resistant plasmid carrying 19 resistance genes for 9 different antimicrobial families. In this plasmid, the blaNDM-23 gene is in the variable region of a complex class 1 integron with a singular genetic environment. The small genetic distance between blaNDM-23-producing isolates reflects a 5-year-long clonal dispersion involving several hospitals and interregional spread. We have characterized the genomic and epidemiological contexts in the emergence and community spread of a new blaNDM-23 allele in a multidrug resistance (MDR) plasmid of Klebsiella pneumoniae. IMPORTANCE At a time when antimicrobial resistance has become one of the biggest concerns worldwide, the emergence of novel alleles and extremely drug-resistant plasmids is a threat to public health worldwide, especially when they produce carbapenem resistance in one of the most problematic pathogens, such as Klebsiella pneumoniae. We used genomic epidemiology to describe the emergence of a novel NDM-23 allele and identify it in a MDR plasmid that has disseminated through a K. pneumoniae ST437 clone in several hospitals in Spain. Using bioinformatic and phylogenetic analyses, we have traced the evolutionary and epidemiological route of the new allele, the hosting plasmid, and the strain that carried both of them from Pakistan to Spain. A better understanding of the NDM-producing K. pneumoniae populations and plasmids has made evident the spread of this clone through the region, enhancing the importance of genomic surveillance in the control of antimicrobial resistance.

Learning From Mistakes: The Role of Phages in Pandemics.

The misuse of antibiotics is leading to the emergence of multidrug-resistant (MDR) bacteria, and in the absence of available treatments, this has become a major global threat. In the middle of the recent severe acute respiratory coronavirus 2 (SARS-CoV-2) pandemic, which has challenged the whole world, the emergence of MDR bacteria is increasing due to prophylactic administration of antibiotics to intensive care unit patients to prevent secondary bacterial infections. This is just an example underscoring the need to seek alternative treatments against MDR bacteria. To this end, phage therapy has been proposed as a promising tool. However, further research in the field is mandatory to assure safety protocols and to develop appropriate regulations for its use in clinics. This requires investing more in such non-conventional or alternative therapeutic approaches, to develop new treatment regimens capable of reducing the emergence of MDR and preventing future global public health concerns that could lead to incalculable human and economic losses.

One is not enough: On the effects of reference genome for the mapping and subsequent analyses of short-reads.

Mapping of high-throughput sequencing (HTS) reads to a single arbitrary reference genome is a frequently used approach in microbial genomics. However, the choice of a reference may represent a source of errors that may affect subsequent analyses such as the detection of single nucleotide polymorphisms (SNPs) and phylogenetic inference. In this work, we evaluated the effect of reference choice on short-read sequence data from five clinically and epidemiologically relevant bacteria (Klebsiella pneumoniae, Legionella pneumophila, Neisseria gonorrhoeae, Pseudomonas aeruginosa and Serratia marcescens). Publicly available whole-genome assemblies encompassing the genomic diversity of these species were selected as reference sequences, and read alignment statistics, SNP calling, recombination rates, dN/dS ratios, and phylogenetic trees were evaluated depending on the mapping reference. The choice of different reference genomes proved to have an impact on almost all the parameters considered in the five species. In addition, these biases had potential epidemiological implications such as including/excluding isolates of particular clades and the estimation of genetic distances. These findings suggest that the single reference approach might introduce systematic errors during mapping that affect subsequent analyses, particularly for data sets with isolates from genetically diverse backgrounds. In any case, exploring the effects of different references on the final conclusions is highly recommended.

Isolation and Characterization of Two Klebsiella pneumoniae Phages Encoding Divergent Depolymerases.

The emergence of multidrug-resistant bacteria is a major global health concern. The search for new therapies has brought bacteriophages into the spotlight, and new phages are being described as possible therapeutic agents. Among the bacteria that are most extensively resistant to current antibiotics is Klebsiella pneumoniae, whose hypervariable extracellular capsule makes treatment particularly difficult. Here, we describe two new K. pneumoniae phages, πVLC5 and πVLC6, isolated from environmental samples. These phages belong to the genus Drulisvirus within the family Podoviridae. Both phages encode a similar tail spike protein with putative depolymerase activity, which is shared among other related phages and probably determines their ability to specifically infect K. pneumoniae capsular types K22 and K37. In addition, we found that phage πVLC6 also infects capsular type K13 and is capable of striping the capsules of K. pneumoniae KL2 and KL3, although the phage was not infectious in these two strains. Genome sequence analysis suggested that the extended tropism of phage πVLC6 is conferred by a second, divergent depolymerase. Phage πVLC5 encodes yet another putative depolymerase, but we found no activity of this phage against capsular types other than K22 and K37, after testing a panel of 77 reference strains. Overall, our results confirm that most phages productively infected one or few Klebsiella capsular types. This constitutes an important challenge for clinical applications.

Isolation of Four Lytic Phages Infecting Klebsiella pneumoniae K22 Clinical Isolates from Spain.

The emergence of multi-drug-resistant bacteria represents a major public-health threat. Phages constitute a promising alternative to chemical antibiotics due to their high host specificity, abundance in nature, and evolvability. However, phage host specificity means that highly diverse bacterial species are particularly difficult to target for phage therapy. This is the case of Klebsiella pneumoniae, which presents a hypervariable extracellular matrix capsule exhibiting dozens of variants. Here, we report four novel phages infecting K. pneumoniae capsular type K22 which were isolated from environmental samples in Valencia, Spain. Full genome sequencing showed that these phages belong to the Podoviridae family and encode putative depolymerases that allow digestion of specific K22 K. pneumoniae capsules. Our results confirm the capsular type-specificity of K. pneumoniae phages, as indicated by their narrow infectivity in a panel of K. pneumoniae clinical isolates. Nonetheless, this work represents a step forward in the characterization of phage diversity, which may culminate in the future use of large panels of phages for typing and/or for combating multi-drug-resistant K. pneumoniae.

Characterization of New Recombinant Forms of HIV-1 From the Comunitat Valenciana (Spain) by Phylogenetic Incongruence.

Recombination is one of the main processes shaping the evolution of HIV-1, with relevant consequences for its epidemiology. In fact, Circulating and Unique Recombinant Forms (CRFs and URFs) cause 23% of current infections. The routine analyses of antiretroviral resistance yield partial pol gene sequences that can be exploited for molecular epidemiology surveillance but also to study viral diversity and to detect potential recombinant samples. Among the pol sequences derived from a large sample dataset from the Comunitat Valenciana (Spain), we identified nine putative recombinant samples. We aimed at fully characterizing these samples and performing a detailed analysis of the corresponding recombination events. We obtained nearly full-genome sequences and used jpHMM and RDP4 to detect and characterize recombinant fragments. We assessed the confidence of these inferences by likelihood mapping and phylogenetic placement with topology congruence tests. Next, we performed a phylogenetic analysis of each putative recombinant fragment to determine its relationships to previously described recombinant forms. We found that two samples related to CRF44_BF whereas the rest corresponded to new URFs (two URF_AD, one URF_BG that can constitute a new CRF resulting from subtype B and CRF24_BG, and two URF_cpx composed of A, G, K, H, and J subtypes). These URFs have a complex recombination pattern that cannot be determined accurately. They seem to have arisen by successive recombination events among lineages, including other CRFs. Our results highlight the usefulness of routine surveillance analysis for the detection of new HIV-1 recombination forms and, at the same time, the need for full-genome sequencing and recombination detection guidelines to properly characterize this complex process.