Prophage-Driven Genomic Structural Changes Promote Bartonella Vertical Evolution

Exploring the diversity and evolutionary strategies of prophages in Hyphomicrobiales, comparing animal-associated with non-animal-associated bacteria

The Hyphomicrobiales bacterial order (previously Rhizobiales) exhibits a wide range of lifestyle characteristics, including free-living, plant-association, nitrogen-fixing, and association with animals (Bartonella and Brucella). This study explores the diversity and evolutionary strategies of bacteriophages within the Hyphomicrobiales order, comparing animal-associated (AAB) with non-animal-associated bacteria (NAAB). We curated 560 high-quality complete genomes of 58 genera from this order and used the PHASTER server for prophage annotation and classification. For 19 genera with representative genomes, we curated 96 genomes and used the Defense-Finder server to summarize the type of anti-phage systems (APS) found in this order. We analyzed the genetic repertoire and length distributions of prophages, estimating evolutionary rates and comparing intact, questionable, and incomplete prophages in both groups. Analyses of best-fit parameters and bootstrap sensitivity were used to understand the evolutionary processes driving prophage gene content. A total of 1860 prophages distributed in Hyphomicrobiales were found, 695 in AAB and 1165 in the NAAB genera. The results revealed a similar number of prophages per genome in AAB and NAAB and a similar length distribution, suggesting shared mechanisms of genetic acquisition of prophage genes. Changes in the frequency of specific gene classes were observed between incomplete and intact prophages, indicating preferential loss or enrichment in both groups. The analysis of best-fit parameters and bootstrap sensitivity tests indicated a higher selection coefficient, induction rate, and turnover in NAAB genomes. We found 68 types of APS in Hyphomicrobiales; restriction modification (RM) and abortive infection (Abi) were the most frequent APS found for all Hyphomicrobiales, and within the AAB group. This classification of APS showed that NAAB genomes have a greater diversity of defense systems compared to AAB, which could be related to the higher rates of prophage induction and turnover in the latter group. Our study provides insights into the distributions of both prophages and APS in Hyphomicrobiales genomes, demonstrating that NAAB carry more defense systems against phages, while AAB show increased prophage stability and an increased number of incomplete prophages. These results suggest a greater role for domesticated prophages within animal-associated bacteria in Hyphomicrobiales.

Genomic structural plasticity of rodent-associated Bartonella in nature

Rodent‐associated Bartonella species have shown a remarkable genetic diversity and pathogenic potential. To further explore the extent of the natural intraspecific genomic variation and its potential role as an evolutionary driver, we focused on a single genetically diverse Bartonella species, Bartonella krasnovii, which circulates among gerbils and their associated fleas. Twenty genomes from 16 different B. krasnovii genotypes were fully characterized through a genome sequencing assay (using short and long read sequencing), pulse field gel electrophoresis (PFGE), and PCR validation. Genomic analyses were performed in comparison to the B. krasnovii strain OE 1–1. While, single nucleotide polymorphism represented only a 0.3% of the genome variation, structural diversity was identified in these genomes, with an average of 51 ± 24 structural variation (SV) events per genome. Interestingly, a large proportion of the SVs (>40%) was associated with prophages. Further analyses revealed that most of the SVs, and prophage insertions were found at the chromosome replication termination site (ter), suggesting this site as a plastic zone of the B. krasnovii chromosome. Accordingly, six genomes were found to be unbalanced, and essential genes near the ter showed a shift between the leading and lagging strands, revealing the SV effect on these genomes. In summary, our findings demonstrate the extensive genomic diversity harbored by wild B. krasnovii strains and suggests that its diversification is initially promoted by structural changes, probably driven by phages. These events may constantly feed the system with novel genotypes that ultimately lead to inter‐ and intraspecies competition and adaptation.

Deciphering Active Prophages from Metagenomes

Temperate phages (prophages) are ubiquitous in nature and persist as dormant components of host cells (lysogenic stage) before activating and lysing the host (lytic stage). Actively replicating prophages contribute to central community processes, such as enabling bacterial virulence, manipulating biogeochemical cycling, and driving microbial community diversification. Recent advances in sequencing technology have allowed for the identification and characterization of diverse phages, yet no approaches currently exist for identifying if a prophage has activated. Here, we present PropagAtE (Prophage Activity Estimator), an automated software tool for estimating if a prophage is in the lytic or lysogenic stage of infection. PropagAtE uses statistical analyses of prophage-to-host read coverage ratios to decipher actively replicating prophages, irrespective of whether prophages were induced or spontaneously activated. We demonstrate that PropagAtE is fast, accurate, and sensitive, regardless of sequencing depth. Application of PropagAtE to prophages from 348 complex metagenomes from human gut, murine gut, and soil environments identified distinct spatial and temporal prophage activation signatures, with the highest proportion of active prophages in murine gut samples. In infants treated with antibiotics or infants without treatment, we identified active prophage populations correlated with specific treatment groups. Within time series samples from the human gut, 11 prophage populations, some encoding the sulfur metabolism gene cysH or a rhuM-like virulence factor, were consistently present over time but not active. Overall, PropagAtE will facilitate accurate representations of viruses in microbiomes by associating prophages with their active roles in shaping microbial communities in nature.

IMPORTANCE Viruses that infect bacteria are key components of microbiomes and ecosystems. They can kill and manipulate microorganisms, drive planetary-scale processes and biogeochemical cycling, and influence the structures of entire food networks. Prophages are viruses that can exist in a dormant state within the genome of their host (lysogenic stage) before activating in order to replicate and kill the host (lytic stage). Recent advances have allowed for the identification of diverse viruses in nature, but no approaches exist for characterizing prophages and their stages of infection (prophage activity). We develop and benchmark an automated approach, PropagAtE, to identify the stages of infection of prophages from genomic data. We provide evidence that active prophages vary in identity and abundance across multiple environments and scales. Our approach will enable accurate and unbiased analyses of viruses in microbiomes and ecosystems.

Host-Pathogen Adhesion as the Basis of Innovative Diagnostics for Emerging Pathogens

Infectious diseases are an existential health threat, potentiated by emerging and re-emerging viruses and increasing bacterial antibiotic resistance. Targeted treatment of infectious diseases requires precision diagnostics, especially in cases where broad-range therapeutics such as antibiotics fail. There is thus an increasing need for new approaches to develop sensitive and specific in vitro diagnostic (IVD) tests. Basic science and translational research are needed to identify key microbial molecules as diagnostic targets, to identify relevant host counterparts, and to use this knowledge in developing or improving IVD. In this regard, an overlooked feature is the capacity of pathogens to adhere specifically to host cells and tissues. The molecular entities relevant for pathogen–surface interaction are the so-called adhesins. Adhesins vary from protein compounds to (poly-)saccharides or lipid structures that interact with eukaryotic host cell matrix molecules and receptors. Such interactions co-define the specificity and sensitivity of a diagnostic test. Currently, adhesin-receptor binding is typically used in the pre-analytical phase of IVD tests, focusing on pathogen enrichment. Further exploration of adhesin–ligand interaction, supported by present high-throughput “omics” technologies, might stimulate a new generation of broadly applicable pathogen detection and characterization tools. This review describes recent results of novel structure-defining technologies allowing for detailed molecular analysis of adhesins, their receptors and complexes. Since the host ligands evolve slowly, the corresponding adhesin interaction is under selective pressure to maintain a constant receptor binding domain. IVD should exploit such conserved binding sites and, in particular, use the human ligand to enrich the pathogen. We provide an inventory of methods based on adhesion factors and pathogen attachment mechanisms, which can also be of relevance to currently emerging pathogens, including SARS-CoV-2, the causative agent of COVID-19.

Bartonella kosoyi sp. nov. and Bartonella krasnovii sp. nov., two novel species closely related to the zoonotic Bartonella elizabethae, isolated from black rats and wild desert rodent-fleas

The genus Bartonella (Family: Bartonellaceae; Order: Rhizobiales; Class: Alphaproteobacteria) comprises facultative intracellular Gram-negative, haemotropic, slow-growing, vector-borne bacteria. Wild rodents and their fleas harbor a great diversity of species and strains of the genus Bartonella, including several zoonotic ones. This genetic diversity coupled with a fastidious nature of the organism results in a taxonomic challenge that has led to a massive collection of uncharacterized strains. Here, we report the genomic and phenotypic characterization of two strains, members of the genus Bartonella (namely Tel Aviv and OE 1-1), isolated from Rattus rattus rats and Synosternus cleopatrae fleas, respectively. Scanning electron microscopy revealed rod-shaped bacteria with polar pili, lengths ranging from 1.0 to 2.0 µm and widths ranging from 0.3 to 0.6 µm. OE 1-1 and Tel Aviv strains contained one single chromosome of 2.16 and 2.23 Mbp and one plasmid of 29.0 and 41.5 Kbp, with average DNA G+C contents of 38.16 and 38.47 mol%, respectively. These strains presented an average nucleotide identity (ANI) of 89.9 %. Bartonella elizabethae was found to be the closest phylogenetic relative of both strains (ANI=90.9-93.6 %). The major fatty acids identified in both strains were C_18:1ω7c, C_18 : 0 and C_16 : 0. They differ from B. elizabethae in their C_17 : 0 and C_15 : 0 compositions. Both strains are strictly capnophilic and their biochemical profiles resembled those of species of the genus Bartonella with validly published names, whereas differences in arylamidase activities partially assisted in their speciation. Genomic and phenotypic differences demonstrate that OE 1-1 and Tel Aviv strains represent novel individual species, closely related to B. elizabethae, for which we propose the names Bartonella kosoyi sp. nov. and Bartonella krasnovii sp. nov.

Bartonella gene transfer agent: Evolution, function, and proposed role in host adaptation

The processes underlying host adaptation by bacterial pathogens remain a fundamental question with relevant clinical, ecological, and evolutionary implications. Zoonotic pathogens of the genus Bartonella constitute an exceptional model to study these aspects. Bartonellae have undergone a spectacular diversification into multiple species resulting from adaptive radiation. Specific adaptations of a complex facultative intracellular lifestyle have enabled the colonisation of distinct mammalian reservoir hosts. This remarkable host adaptability has a multifactorial basis and is thought to be driven by horizontal gene transfer (HGT) and recombination among a limited genus‐specific pan genome. Recent functional and evolutionary studies revealed that the conserved Bartonella gene transfer agent (BaGTA) mediates highly efficient HGT and could thus drive this evolution. Here, we review the recent progress made towards understanding BaGTA evolution, function, and its role in the evolution and pathogenesis of Bartonella spp. We notably discuss how BaGTA could have contributed to genome diversification through recombination of beneficial traits that underlie host adaptability. We further address how BaGTA may counter the accumulation of deleterious mutations in clonal populations (Muller's ratchet), which are expected to occur through the recurrent transmission bottlenecks during the complex infection cycle of these pathogens in their mammalian reservoir hosts and arthropod vectors.