Evolution along the parasitism-mutualism continuum determines the genetic repertoire of prophages

Pseudomonas aeruginosa maintains an inducible array of novel and diverse prophages over lengthy persistence in cystic fibrosis lungs

Pseudomonas aeruginosa has increasing clinical relevance and commonly occupies the cystic fibrosis (CF) airways. Its ability to colonize and persist in diverse niches is attributed to its large accessory genome, where prophages represent a common feature and may contribute to its fitness and persistence. We focused on the CF airways niche and used 197 longitudinal isolates from 12 patients persistently infected by P. aeruginosa. We computationally predicted intact prophages for each longitudinal group and scored their long-term persistence. We then confirmed prophage inducibility and mapped their location in the host chromosome with lysate sequencing. Using comparative genomics, we evaluated prophage genomic diversity, long-term persistence, and level of genomic maintenance. Our findings support previous findings that most P. aeruginosa genomes harbour prophages some of which can self-induce, and that a common CF-treating antibiotic, ciprofloxacin, can induce prophages. Induced prophage genomes displayed high diversity and even genomic novelty. Finally, all induced prophages persisted long-term with their genomes avoiding gene loss and degradation over 4 years of host replication in the stressful CF airways niche. This and our detection of phage genes, which contribute to host competitiveness and adaptation, lends support to our hypothesis that the vast majority of prophages detected as intact and inducible in this study facilitated their host fitness and persistence.

Escherichia coli phage-inducible chromosomal island aids helper phage replication and represses the locus of enterocyte effacement pathogenicity island

In this study, we identify and characterize a novel phage-inducible chromosomal island (PICI) found in commensal Escherichia coli MP1. This novel element, EcCIMP1, is induced and mobilized by the temperate helper phage vB_EcoP_Kapi1. EcCIMP1 contributes to superinfection immunity against its helper phage, impacting bacterial competition outcomes. Genetic analysis of EcCIMP1 led us to uncover a putative transcriptional repressor, which silences virulence gene expression in the murine pathogen Citrobacter rodentium. We also found a putative excisionase encoded by EcCIMP1 which paradoxically does not promote excision of EcCIMP1 but rather supports excision of the helper phage. Another putative excisionase encoded by a presumed integrative conjugative element can also support the excision of vB_EcoP_Kapi1, demonstrating crosstalk between excisionases from multiple classes of mobile genetic elements within the same cell. Although phylogenetically distant from other characterized PICIs, EcCIMP1 and EcCIMP1-like elements are prevalent in both pathogenic and commensal isolates of E. coli from around the world, underscoring the importance of characterizing these abundant genetic elements.

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.

Prophage maintenance is determined by environment-dependent selective sweeps rather than mutational availability

Prophages, viral sequences integrated into bacterial genomes, can be beneficial and costly. Despite the risk of prophage activation and subsequent bacterial death, active prophages are present in most bacterial genomes. However, our understanding of the selective forces that maintain prophages in bacterial populations is limited. Combining experimental evolution with stochastic modeling, we show that prophage maintenance and loss are primarily determined by environmental conditions that alter the net fitness effect of a prophage on its bacterial host. When prophages are too costly, they are rapidly lost through environment-specific sequences of selective sweeps. Conflicting selection pressures that select against the prophage but for a prophage-encoded accessory gene can maintain prophages. The dynamics of prophage maintenance additionally depend on the sociality of this accessory gene. Prophage-encoded genes that exclusively benefit the lysogen maintain prophages at higher frequencies compared with genes that benefit the entire population. That is because the latter can protect phage-free "cheaters," reducing the benefit of maintaining the prophage. Our simulations suggest that environmental variation plays a larger role than mutation rates in determining prophage maintenance. These findings highlight the complexity of selection pressures that act on mobile genetic elements and challenge our understanding of the role of environmental factors relative to random chance events in shaping the evolutionary trajectory of bacterial populations. By shedding light on the key factors that shape microbial populations in the face of environmental changes, our study significantly advances our understanding of the complex dynamics of microbial evolution and diversification.

Endogenous virophages are active and mitigate giant virus infection in the marine protist Cafeteria burkhardae

Endogenous viral elements (EVEs) are common genetic passengers in various protists. Some EVEs represent viral fossils, whereas others are still active. The marine heterotrophic flagellate Cafeteria burkhardae contains several EVE types related to the virophage mavirus, a small DNA virus that parasitizes the lytic giant virus CroV. We hypothesized that endogenous virophages may act as an antiviral defense system in protists, but no protective effect of virophages in wild host populations has been shown so far. Here, we tested the activity of virophage EVEs and studied their impact on giant virus replication. We found that endogenous mavirus-like elements (EMALEs) from globally distributed Cafeteria populations produced infectious virus particles specifically in response to CroV infection. However, reactivation was stochastic, often inefficient, and poorly reproducible. Interestingly, only one of eight EMALE types responded to CroV infection, implying that other EMALEs may be linked to different giant viruses. We isolated and cloned several reactivated virophages and characterized their particles, genomes, and infection dynamics. All tested virophages inhibited the production of CroV during coinfection, thereby preventing lysis of the host cultures in a dose-dependent manner. Comparative genomics of different C. burkhardae strains revealed that inducible EMALEs are common and are not linked to specific geographic locations. We demonstrate that naturally occurring virophage EVEs reactivate upon giant virus infection, thus providing a striking example that eukaryotic EVEs can become active under specific conditions. Moreover, our results support the hypothesis that virophages can act as an adaptive antiviral defense system in protists.

Analysis of CRISPR/Cas Genetic Structure, Spacer Content and Molecular Epidemiology in Brazilian Acinetobacter baumannii Clinical Isolates

CRISPR/Cas is a molecular mechanism to prevent predatory viruses from invading bacteria via the insertion of small viral sequences (spacers) in its repetitive locus. The nature of spacer incorporation and the viral origins of spacers provide an overview of the genetic evolution of bacteria, their natural viral predators, and the mechanisms that prokaryotes may use to protect themselves, or to acquire mobile genetic elements such as plasmids. Here, we report on the CRISPR/Cas genetic structure, its spacer content, and strain epidemiology through MLST and CRISPR typing in Acinetobacter baumannii, an opportunistic pathogen intimately related to hospital infections and antimicrobial resistance. Results show distinct genetic characteristics, such as polymorphisms specific to ancestor direct repeats, a well-defined degenerate repeat, and a conserved leader sequence, as well as showing most spacers as targeting bacteriophages, and several self-targeting spacers, directed at prophages. There was a particular relationship between CRISPR/Cas and CC113 in the study of Brazilian isolates, and CRISPR-related typing techniques are interesting for subtyping strains with the same MLST profile. We want to emphasize the significance of descriptive genetic research on CRISPR loci, and we argue that spacer or CRISPR typing are helpful for small-scale investigations, preferably in conjunction with other molecular typing techniques such as MLST.

The Prophage and Us-Shiga Toxin Phages Revisited

The authors first met in 1998 at the University of Würzburg, Germany, at the Institute of Hygiene and Microbiology, in Helge Karch's lab, where Herbert Schmidt worked as a PostDoc and Maite Muniesa visited the lab for a postdoctoral research stay to work on phages encoding Shiga toxin 2e (Stx2e) [...].

Comparative Genome Analysis of 19 Trueperella pyogenes Strains Originating from Different Animal Species Reveal a Genetically Diverse Open Pan-Genome

Trueperella pyogenes is a Gram-positive opportunistic pathogen that causes severe cases of mastitis, metritis, and pneumonia in a wide range of animals, resulting in significant economic losses. Although little is known about the virulence factors involved in the disease pathogenesis, a comprehensive comparative genome analysis of T. pyogenes genomes has not been performed till date. Hence, present investigation was carried out to characterize and compare 19 T. pyogenes genomes originating in different geographical origins including the draftgenome of the first Indian origin strain T. pyogenes Bu5. Additionally, candidate virulence determinants that could be crucial for their pathogenesis were also detected and analyzed by using various bioinformatics tools. The pan-genome calculations revealed an open pan-genome of T. pyogenes. In addition, an inventory of virulence related genes, 190 genomic islands, 31 prophage sequences, and 40 antibiotic resistance genes that could play a significant role in organism's pathogenicity were detected. The core-genome based phylogeny of T. pyogenes demonstrates a polyphyletic, host-associated group with a high degree of genomic diversity. The identified core-genome can be further used for screening of drug and vaccine targets. The investigation has provided unique insights into pan-genome, virulome, mobiliome, and resistome of T. pyogenes genomes and laid the foundation for future investigations.

Novel lineages of single-stranded DNA phages that coevolved with the symbiotic bacteria Rhizobium

This study describes novel single-stranded DNA phages isolated from common bean agriculture soils by infection of the nitrogen-fixing symbiotic bacteria Rhizobium etli and R. phaseoli. A total of 29 phages analyzed have 4.3-6 kb genomes in size and GC 59-60%. They belong to different clades unrelated to other Microviridae subfamilies. Three-dimensional models of the major capsid protein (MCP) showed a conserved β-barrel structural "jelly-roll" fold. A variable-length loop in the MCPs distinguished three Rhizobium microvirus groups. Microviridae subfamilies were consistent with viral clusters determined by the protein-sharing network. All viral clusters, except for Bullavirinae, included mostly microviruses identified in metagenomes from distinct ecosystems. Two Rhizobium microvirus clusters, chaparroviruses, and chicoviruses, were included within large viral unknown clusters with microvirus genomes identified in diverse metagenomes. A third Rhizobium microvirus cluster belonged to the subfamily Amoyvirinae. Phylogenetic analysis of the MCP confirms the divergence of the Rhizobium microviruses into separate clades. The phylogeny of the bacterial hosts matches the microvirus MCP phylogeny, suggesting a coevolutionary history between the phages and their bacterial host. This study provided essential biological information on cultivated microvirus for understanding the evolution and ecological diversification of the Microviridae family in diverse microbial ecosystems.

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.

Host life-history traits influence the distribution of prophages and the genes they carry

Bacterial strains with a short minimal doubling time-'fast-growing' hosts-are more likely to contain prophages than their slow-growing counterparts. Pathogenic bacterial species are likewise more likely to carry prophages. We develop a bioinformatics pipeline to examine the distribution of prophages in fast- and slow-growing lysogens, and pathogenic and non-pathogenic lysogens, analysing both prophage length and gene content for each class. By fitting these results to a mathematical model of the evolutionary forces acting on prophages, we predict whether the observed differences can be attributed to different rates of lysogeny among the host classes, or other evolutionary pressures. We also test for significant differences in gene content among prophages, identifying genes that are preferentially lost or maintained in each class. We find that fast-growing hosts and pathogens have a greater fraction of full-length prophages, and our analysis predicts that induction rates are significantly reduced in slow-growing hosts and non-pathogenic hosts. Consistent with previous results, we find that several proteins involved in the packaging of new phage particles and lysis are preferentially lost in cryptic prophages. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Streptococcus pneumoniae: a Plethora of Temperate Bacteriophages With a Role in Host Genome Rearrangement

Bacteriophages (phages) are viruses that infect bacteria. They are the most abundant biological entity on Earth (current estimates suggest there to be perhaps 10³¹ particles) and are found nearly everywhere. Temperate phages can integrate into the chromosome of their host, and prophages have been found in abundance in sequenced bacterial genomes. Prophages may modulate the virulence of their host in different ways, e.g., by the secretion of phage-encoded toxins or by mediating bacterial infectivity. Some 70% of Streptococcus pneumoniae (the pneumococcus)—a frequent cause of otitis media, pneumonia, bacteremia and meningitis—isolates harbor one or more prophages. In the present study, over 4000 S. pneumoniae genomes were examined for the presence of prophages, and nearly 90% were found to contain at least one prophage, either defective (47%) or present in full (43%). More than 7000 complete putative integrases, either of the tyrosine (6243) or serine (957) families, and 1210 full-sized endolysins (among them 1180 enzymes corresponding to 318 amino acid-long N-acetylmuramoyl-L-alanine amidases [LytA_PPH]) were found. Based on their integration site, 26 different pneumococcal prophage groups were documented. Prophages coding for tRNAs, putative virulence factors and different methyltransferases were also detected. The members of one group of diverse prophages (PPH090) were found to integrate into the 3’ end of the host lytA_Spn gene encoding the major S. pneumoniae autolysin without disrupting it. The great similarity of the lytA_Spnand lytA_PPH genes (85–92% identity) allowed them to recombine, via an apparent integrase-independent mechanism, to produce different DNA rearrangements within the pneumococcal chromosome. This study provides a complete dataset that can be used to further analyze pneumococcal prophages, their evolutionary relationships, and their role in the pathogenesis of pneumococcal disease.

Prophage Genomics and Ecology in the Family Rhodobacteraceae

Roseobacters are globally abundant bacteria with critical roles in carbon and sulfur biogeochemical cycling. Here, we identified 173 new putative prophages in 79 genomes of Rhodobacteraceae. These prophages represented 1.3 ± 0.15% of the bacterial genomes and had no to low homology with reference and metagenome-assembled viral genomes from aquatic and terrestrial ecosystems. Among the newly identified putative prophages, 35% encoded auxiliary metabolic genes (AMGs), mostly involved in secondary metabolism, amino acid metabolism, and cofactor and vitamin production. The analysis of integration sites and gene homology showed that 22 of the putative prophages were actually gene transfer agents (GTAs) similar to a GTA of Rhodobacter capsulatus. Twenty-three percent of the predicted prophages were observed in the TARA Oceans viromes generated from free viral particles, suggesting that they represent active prophages capable of induction. The distribution of these prophages was significantly associated with latitude and temperature. The prophages most abundant at high latitudes encoded acpP, an auxiliary metabolic gene involved in lipid synthesis and membrane fluidity at low temperatures. Our results show that prophages and gene transfer agents are significant sources of genomic diversity in roseobacter, with potential roles in the ecology of this globally distributed bacterial group.