Phage Population Research Articles

Leveraging mathematical modeling framework to guide regimen strategy for phage therapy

Bacteriophage (phage) cocktail therapy has been relied upon more and more to treat antibiotic-resistant infections. Understanding of the complex kinetics between phages, target bacteria, and the emergence of phage resistance remain hurdles to successful clinical outcomes. Building upon previous mathematical concepts, we develop biologically-motivated nonlinear ordinary differential equation models to explore single, cocktail, and sequential phage treatment modalities. While the optimal pairwise phage treatment strategy was the double simultaneous administration of two highly potent and asymmetrically binding phage strains, it appears unable to prevent the evolution of resistance. This treatment regimen did have a greater lysis efficiency, promoted higher phage population sizes, reduced bacterial density the most, and suppressed the evolution of resistance the longest compared to all other treatments strategies tested. Conversely, the combination of phages with polar potencies allows the more efficiently replicating phages to monopolize susceptible host cells, thereby quickly negating the intended compounding effect of cocktails. Together, we demonstrate that a biologically-motivated modeling-based framework can be leveraged to quantify the effects of each phage’s properties to more precisely predict treatment responses.

Genetic variation in individuals from a population of the minimalist bacteriophage Merri-merri-uth nyilam marra-natj driving evolution of the virus.

Bacteriophages (phages) are viruses that prey on bacteria. This study sampled natural phage populations to test the hypothesis that untapped genetic variation within a population can be the basis for the selection of phages to diversify their host-range. Sampling of a freshwater site revealed two populations of the phage Merri-merri-uth nyilam marra-natj (phage MMNM), differing by a variant residue (Val134Ala) in the baseplate protein MMNM_26. This sequence variation modulated bacterial killing in plaques, and further evolution of the phages on a semi-permissive bacterial host led to a new generation of phages with more diverse phenotypes in killing the bacterium Klebsiella pneumoniae.

The phageome of patients with ulcerative colitis treated with donor fecal microbiota reveals markers associated with disease remission

Bacteriophages are influential within the human gut microbiota, yet they remain understudied relative to bacteria. This is a limitation of studies on fecal microbiota transplantation (FMT) where bacteriophages likely influence outcome. Here, using metagenomics, we profile phage populations - the phageome - in individuals recruited into two double-blind randomized trials of FMT in ulcerative colitis. We leverage the trial designs to observe that phage populations behave similarly to bacterial populations, showing temporal stability in health, dysbiosis in active disease, modulation by antibiotic treatment and by FMT. We identify a donor bacteriophage putatively associated with disease remission, which on genomic analysis was found integrated in a bacterium classified to Oscillospiraceae, previously isolated from a centenarian and predicted to produce vitamin B complex except B12. Our study provides an in-depth assessment of phage populations during different states and suggests that bacteriophage tracking has utility in identifying determinants of disease activity and resolution.

Alterations of the gut microbiome in HIV infection highlight human anelloviruses as potential predictors of immune recovery

BackgroundHIV-1 infection is characterized by a massive depletion of mucosal CD4 T cells that triggers a cascade of events ultimately linking gut microbial dysbiosis to HIV-1 disease progression and pathogenesis. The association between HIV infection and the enteric virome composition is less characterized, although viruses are an essential component of the gut ecosystem. Here, we performed a cross-sectional analysis of the fecal viral (eukaryotic viruses and bacteriophages) and bacterial microbiome in people with HIV (PWH) and in HIV-negative individuals. To gain further insight into the association between the gut microbiome composition, HIV-associated immunodeficiency, and immune recovery, we carried out a longitudinal study including 14 PWH who initiated antiretroviral therapy (ART) and were followed for 24 months with samplings performed at baseline (before ART) and at 2, 6, 12, and 24 months post-ART initiation.ResultsOur data revealed a striking expansion in the abundance and prevalence of several human virus genomic sequences (Anelloviridae, Adenoviridae, and Papillomaviridae) in stool samples of PWH with severe immunodeficiency (CD4 < 200). We also noted a decreased abundance of sequences belonging to two plant viruses from the Tobamovirus genus, a reduction in bacterial alpha diversity, and a decrease in Inoviridae bacteriophage sequences. Short-term ART (24 months) was linked to a significant decrease in human Anelloviridae sequences. Remarkably, the detection of Anellovirus sequences at baseline independently predicted poor immune recovery, as did low CD4 T cell counts. The bacterial and bacteriophage populations were unique to each PWH with individualized trajectories; we found no discernable pattern of clustering after 24 months on ART.ConclusionAdvanced HIV-1 infection was associated with marked alterations in the virome composition, in particular a remarkable expansion of human anelloviruses, with a gradual restoration after ART initiation. In addition to CD4 T cell counts, anellovirus sequence detection might be useful to predict and monitor immune recovery. This study confirms data on the bacteriome and expands our knowledge on the viral component of the gut microbiome in HIV-1 infection.3oR8qDfXoWY9s_P5ES648YVideo

Improving viral annotation with artificial intelligence.

Viruses of bacteria, "phages," are fundamental, poorly understood components of microbial community structure and function. Additionally, their dependence on hosts for replication positions phages as unique sensors of ecosystem features and environmental pressures. High-throughput sequencing approaches have begun to give us access to the diversity and range of phage populations in complex microbial community samples, and metagenomics is currently the primary tool with which we study phage populations. The study of phages by metagenomic sequencing, however, is fundamentally limited by viral diversity, which results in the vast majority of viral genomes and metagenome-annotated genomes lacking annotation. To harness bacteriophages for applications in human and environmental health and disease, we need new methods to organize and annotate viral sequence diversity. We recently demonstrated that methods that leverage self-supervised representation learning can supplement statistical sequence representations for remote viral protein homology detection in the ocean virome and propose that consideration of the functional content of viral sequences allows for the identification of similarity in otherwise sequence-diverse viruses and viral-like elements for biological discovery. In this review, we describe the potential and pitfalls of large language models for viral annotation. We describe the need for new approaches to annotate viral sequences in metagenomes, the fundamentals of what protein language models are and how one can use them for sequence annotation, the strengths and weaknesses of these models, and future directions toward developing better models for viral annotation more broadly.

Phage-Based Therapy in Combination with Antibiotics: A Promising Alternative against Multidrug-Resistant Gram-Negative Pathogens.

The continued rise in antimicrobial resistance poses a serious threat to public health worldwide. The use of phages that can have bactericidal activity without disrupting the normal flora represents a promising alternative treatment method. This practice has been successfully applied for decades, mainly in Eastern Europe, and has recently been used as an emergency therapy for compassionate care in the United States. Here, we provide a comprehensive review of the pre-clinical and clinical applications of phage therapy concerning three major Gram-negative pathogens: Klebsiella pneumoniae, Pseudomonas aeruginosa, and Acinetobacter baumannii. The advantages and the challenges of expanding the usage of phages as an alternative or adjunctive treatment for antimicrobial-resistant bacterial infections are discussed. We emphasize the virologic complexities of using the highly adaptable phage populations as molecular tools, along with antibiotic chemical compounds, to effectively combat rapidly coevolving pathogenic bacteria in the host microenvironment. Pre-clinical studies, isolated clinical reports and a few randomized clinical trials have shown that bacteriophages can be effective in treating multidrug-resistant bacterial infections. The ability of some phages to revert the resistance against antibiotics, and possibly also against the human complement and other phages, appears to be a great advantage of phage therapy despite the inevitable emergence of phage-resistant strains. Bacteriophages (or specific phage-derived products) can enhance antimicrobial efficacy by reducing bacterial virulence via the alteration of basic bacterial structures, primarily of the cellular wall and membrane. Although several issues remain open regarding their effective clinical application, it appears that phage-based therapeutics in combination with antibiotics can provide an effective solution to the spread of antimicrobial resistance.

Specific aspects of bacteriophage regulation: international practices and future developments

INTRODUCTION. The increasing prevalence of multidrug-resistant strains of pathogens determines the need for fundamentally new antibacterial agents, including bacteriophage preparations. The consistent implementation of phage therapy is hindered by the lack of generally accepted standardised regulatory documents governing the legal and methodological aspects of the production and preclinical and clinical studies of bacteriophage preparations.AIM. This study aimed to analyse the international experience with the production and lifecycle management of bacteriophage preparations, as well as the main regulatory requirements for the control of their quality, safety, and efficacy.DISCUSSION. It is difficult to develop virulent bacteriophage preparations in accordance with the existing requirements for other medicinal products because of the biological characteristics of bacteriophages, the wide variety of bacteriophage strains, and the potential for rapid changes both in the bacteriophage population and in the pathogen population. Therefore, it is reasonable to develop streamlined marketing authorisation routes for phage therapies and methods for the assessment of their safety and efficacy. As part of these efforts, it is necessary to assess the adverse events specific to this group of medicinal products, such as the risks of lysogeny, resistance to bacteriophages, and antibiotic resistance gene transfer between bacterial strains. The pharmaceutical development of bacteriophage preparations can be based on several approaches. Many countries worldwide, including the United States, are implementing the concept of Quality by Design, considering approaches based on the Biological Master File, and conducting Expanded Access programmes. The Active Substance Master File procedure allows the submission of a separate document package covering only part of the registration dossier for regulatory approval. Expanded Access programmes provide individual patients with access to innovative medicinal products without approved treatment protocols. In the Russian Federation, the commercial production of bacteriophage medicinal products is conducted in accordance with the quality standards specified in the State Pharmacopoeia of the Russian Federation.CONCLUSIONS. There are fundamental differences in the approaches to phage therapy and its regulation around the world and in the Russian Federation. It is reasonable to supplement the current national guidelines for the safety and efficacy evaluation of bacteriophage preparations, in particular, to specify the requirements for conducting preclinical studies.

Abundance measurements reveal the balance between lysis and lysogeny in the human gut microbiome.

The human gut contains diverse communities of bacteriophage, whose interactions with the broader microbiome and potential roles in human health are only beginning to be uncovered. Here, we combine multiple types of data to quantitatively estimate gut phage population dynamics and lifestyle characteristics in human subjects. Unifying results from previous studies, we show that an average human gut contains a low ratio of phage particles to bacterial cells (~1:100), but a much larger ratio of phage genomes to bacterial genomes (~4:1), implying that most gut phage are effectively temperate (e.g., integrated prophage, phage-plasmids, etc.). By integrating imaging and sequencing data with a generalized model of temperate phage dynamics, we estimate that phage induction and lysis occurs at a low average rate (~0.001-0.01 per bacterium per day), imposing only a modest fitness burden on their bacterial hosts. Consistent with these estimates, we find that the phage composition of a diverse synthetic community in gnotobiotic mice can be quantitatively predicted from bacterial abundances alone, while still exhibiting phage diversity comparable to native human microbiomes. These results provide a foundation for interpreting existing and future studies on links between the gut virome and human health.

Viral adaptations to alternative hosts in the honey bee pathogen Paenibacillus larvae.

Bacteriophages (phages) that are intended to be used to treat bacterial infections are often improved using genetic engineering or experimental evolution. A protocol called "Appelmans" utilizes evolution in microtiter plates to promote the evolution of phages that can infect nonpermissive hosts. We tested a modification of the Appelmans protocol using the honey bee pathogen, Paenibacillus larvae. Three phages evolved together on four P. larvae strains following the standard Appelmans protocol and a modified version to ensure high phage diversity throughout ten rounds of passaging. The host range of 360 plaques were characterized and six new phage lysis patterns were identified. These new phage lysis patterns included plaque formation on previously nonpermissive, phage-resistant isolates that were used to identify phage types. The modified protocol did not drastically change the rate or number of new phage types observed but did prevent the phage population from being dominated by one phage that tended to rapidly raise in frequency. These findings showed how a minor modification of the Appelmans protocol influenced the development of phages for phage therapy. The method also provided improved phages for the treatment of bacterial infections in honey bees.

Design, optimization, and inference of biphasic decay of infectious virus particles.

Virus population dynamics are driven by counter-balancing forces of production and loss. Whereas viral production arises from complex interactions with susceptible hosts, the loss of infectious virus particles is often approximated as a first-order kinetic process. As such, experimental protocols to measure infectious virus loss are not typically designed to identify non-exponential decay processes. Here, we propose methods to evaluate if an experimental design is adequate to identify multiphasic virus particle decay and to optimize the sampling times of decay experiments, accounting for uncertainties in viral kinetics. First, we evaluate synthetic scenarios of biphasic decays, with varying decay rates and initial proportions of subpopulations. We show that robust inference of multiphasic decay is more likely when the faster decaying subpopulation predominates insofar as early samples are taken to resolve the faster decay rate. Moreover, design optimization involving non-equal spacing between observations increases the precision of estimation while reducing the number of samples. We then apply these methods to infer multiple decay rates associated with the decay of bacteriophage ('phage') , an evolved isolate derived from phage . A pilot experiment confirmed that decay is multiphasic, but was unable to resolve the rate or proportion of the fast decaying subpopulation(s). We then applied a Fisher information matrix-based design optimization method to propose non-equally spaced sampling times. Using this strategy, we were able to robustly estimate multiple decay rates and the size of the respective subpopulations. Notably, we conclude that the vast majority (94%) of the phage population decays at a rate 16-fold higher than the slow decaying population. Altogether, these results provide both a rationale and a practical approach to quantitatively estimate heterogeneity in viral decay.

Genetically similar Xanthomonas arboricola pv. pruni strains and associated phage display phenotypic and genotypic variation across 35 years

Bacteriophages (phages), viruses that infect bacteria, have key ecological and evolutionary functions in the phytobiome. Despite the importance of phages as primary drivers for bacterial evolution, phage-bacteria interactions across spatiotemporal scales in natural, agricultural settings are underexplored. With increased interest in phage-based therapies to manage bacterial pathogens, an enhanced understanding of phage genetic and functional diversity at the population level, and how this, in turn, impacts bacterial evolution and virulence, is necessary. This study presents data on the genetic similarity among Xanthomonas arboricola pv. pruni (Xap) strains isolated from different geographic locations that display different lytic phenotypes when challenged with a panel of six phage isolates collected in the same region over four decades. The minor yet significant genetic variation among this small population of Xap strains is structured by both geographic location and response to phage infection. Phage genomes are also highly similar, with conserved and diverse genomic loci that correspond to isolation year. The six phages characterized here cluster into the Kantovirinae subfamily and possibly form a new genus. Only future studies will elucidate the role of Xap and Xapφ phage genes identified here in the virulence and lysis of Xap and how these, in turn, impact bacterial spot disease outcomes. The research and tripartite pathosystem presented here provides a unique opportunity to investigate the coevolution of phage-phytobacterial pathogen-plant host in depth in an agricultural setting with the potential to monitor the rate at which phage populations contribute to bacterial genetic diversity across geographic and temporal scales.

Distributions, interactions, and dynamics of prokaryotes and phages in a hybrid biological wastewater treatment system

BackgroundUnderstanding the interactions and dynamics of microbiotas within biological wastewater treatment systems is essential for ensuring their stability and long-term sustainability. In this study, we developed a systematic framework employing multi-omics and Hi-C sequencing to extensively investigate prokaryotic and phage communities within a hybrid biofilm and activated sludge system.ResultsWe uncovered distinct distribution patterns, metabolic capabilities, and activities of functional prokaryotes through the analysis of 454 reconstructed prokaryotic genomes. Additionally, we reconstructed a phage catalog comprising 18,645 viral operational taxonomic units (vOTUs) with high length and contiguity using hybrid assembly, and a distinct distribution of phages was depicted between activated sludge (AS) and biofilm. Importantly, 1340 host-phage pairs were established using Hi-C and conventional in silico methods, unveiling the host-determined phage prevalence. The majority of predicted hosts were found to be involved in various crucial metabolic processes, highlighting the potential vital roles of phages in influencing substance metabolism within this system. Moreover, auxiliary metabolic genes (AMGs) related to various categories (e.g., carbohydrate degradation, sulfur metabolism, transporter) were predicted. Subsequent activity analysis emphasized their potential ability to mediate host metabolism during infection. We also profiled the temporal dynamics of phages and their associated hosts using 13-month time-series metagenomic data, further demonstrating their tight interactions. Notably, we observed lineage-specific infection patterns, such as potentially host abundance- or phage/host ratio-driven phage population changes.ConclusionsThe insights gained from this research contribute to the growing body of knowledge surrounding interactions and dynamics of host-phage and pave the way for further exploration and potential applications in the field of microbial ecology.7mTZubWt51_LqpKWPutADaVideo

Fecal microbiota transplantation alters gut phage communities in a clinical trial for obesity

BackgroundFecal microbiota transplantation (FMT) is a therapeutic intervention used to treat diseases associated with the gut microbiome. In the human gut microbiome, phages have been implicated in influencing human health, with successful engraftment of donor phages correlated with FMT treatment efficacy. The impact that gastrointestinal phages exert on human health has primarily been connected to their ability to modulate the bacterial communities in the gut. Nonetheless, how FMT affects recipients’ phage populations, and in turn, how this influences the gut environment, is not yet fully understood. In this study, we investigated the effects of FMT on the phageome composition of participants within the Gut Bugs Trial (GBT), a double-blind, randomized, placebo-controlled trial that investigated the efficacy of FMT in treating obesity and comorbidities in adolescents. Stool samples collected from donors at the time of treatment and recipients at four time points (i.e., baseline and 6 weeks, 12 weeks, and 26 weeks post-intervention), underwent shotgun metagenomic sequencing. Phage sequences were identified and characterized in silico to examine evidence of phage engraftment and to assess the extent of FMT-induced alterations in the recipients’ phageome composition.ResultsDonor phages engrafted stably in recipients following FMT, composing a significant proportion of their phageome for the entire course of the study (33.8 ± 1.2% in females and 33.9 ± 3.7% in males). Phage engraftment varied between donors and donor engraftment efficacy was positively correlated with their phageome alpha diversity. FMT caused a shift in recipients’ phageome toward the donors’ composition and increased phageome alpha diversity and variability over time.ConclusionsFMT significantly altered recipients' phage and, overall, microbial populations. The increase in microbial diversity and variability is consistent with a shift in microbial population dynamics. This proposes that phages play a critical role in modulating the gut environment and suggests novel approaches to understanding the efficacy of FMT in altering the recipient’s microbiome.Trial registrationThe Gut Bugs Trial was registered with the Australian New Zealand Clinical Trials Registry (ACTR N12615001351505). Trial protocol: the trial protocol is available at https://bmjopen.bmj.com/content/9/4/e026174.A4HmAXiKAUhp3nsywzn7_QVideo

Phage-induced efflux down-regulation boosts antibiotic efficacy.

The interactions between a virus and its host vary in space and time and are affected by the presence of molecules that alter the physiology of either the host or the virus. Determining the molecular mechanisms at the basis of these interactions is paramount for predicting the fate of bacterial and phage populations and for designing rational phage-antibiotic therapies. We study the interactions between stationary phase Burkholderia thailandensis and the phage ΦBp-AMP1. Although heterogeneous genetic resistance to phage rapidly emerges in B. thailandensis, the presence of phage enhances the efficacy of three major antibiotic classes, the quinolones, the beta-lactams and the tetracyclines, but antagonizes tetrahydrofolate synthesis inhibitors. We discovered that enhanced antibiotic efficacy is facilitated by reduced antibiotic efflux in the presence of phage. This new phage-antibiotic therapy allows for eradication of stationary phase bacteria, whilst requiring reduced antibiotic concentrations, which is crucial for treating infections in sites where it is difficult to achieve high antibiotic concentrations.

Long-read powered viral metagenomics in the oligotrophic Sargasso Sea

Dominant microorganisms of the Sargasso Sea are key drivers of the global carbon cycle. However, associated viruses that shape microbial community structure and function are not well characterised. Here, we combined short and long read sequencing to survey Sargasso Sea phage communities in virus- and cellular fractions at viral maximum (80 m) and mesopelagic (200 m) depths. We identified 2,301 Sargasso Sea phage populations from 186 genera. Over half of the phage populations identified here lacked representation in global ocean viral metagenomes, whilst 177 of the 186 identified genera lacked representation in genomic databases of phage isolates. Viral fraction and cell-associated viral communities were decoupled, indicating viral turnover occurred across periods longer than the sampling period of three days. Inclusion of long-read data was critical for capturing the breadth of viral diversity. Phage isolates that infect the dominant bacterial taxa Prochlorococcus and Pelagibacter, usually regarded as cosmopolitan and abundant, were poorly represented.

SEA-PHAGES and SEA-GENES: Advancing Virology and Science Education.

Research opportunities for undergraduate students are strongly advantageous, but implementation at a large scale presents numerous challenges. The enormous diversity of the bacteriophage population and a supportive programmatic structure provide opportunities to engage early-career undergraduates in phage discovery, genomics, and genetics. The Science Education Alliance (SEA) is an inclusive Research-Education Community (iREC) providing centralized programmatic support for students and faculty without prior experience in virology at institutions from community colleges to research-active universities to participate in two course-based projects, SEA-PHAGES (SEA Phage Hunters Advancing Genomic and Evolutionary Science) and SEA-GENES (SEA Gene-function Exploration by a Network of Emerging Scientists). Since 2008, the SEA has supported more than 50,000 undergraduate researchers who have isolated more than 23,000 bacteriophages of which more than 4,500 are fully sequenced and annotated. Students have functionally characterized hundreds of phage genes, and the phage collection has fueled the therapeutic use of phages for treatment of Mycobacterium infections. Participation in the SEA promotes student persistence in science education, and its inclusivity promotes a more equitable scientific community.

Characterisation of PVL-Positive Staphylococcus argenteus from the United Arab Emirates.

Staphylococcus argenteus is a recently described staphylococcal species that is related to Staphylococcus aureus but lacks the staphyloxanthin operon. It is able to acquire both resistance markers such as the SCCmec elements and mobile genetic elements carrying virulence-associated genes from S. aureus. This includes those encoding the Panton-Valentine leukocidin (PVL), which is associated mainly with severe and/or recurrent staphylococcal skin and soft tissue infections. Here, we describe the genome sequences of two PVL-positive, mecA-negative S. argenteus sequence type (ST) 2250 isolates from the United Arab Emirates in detail. The isolates were found in a dental clinic in the United Arab Emirates (UAE). Both were sequenced using Oxford Nanopore Technology (ONT). This demonstrated the presence of temperate bacteriophages in the staphylococcal genomes, including a PVL prophage. It was essentially identical to the published sequence of phiSa2wa_st78 (GenBank NC_055048), a PVL phage from an Australian S. aureus clonal complex (CC) 88 isolate. Besides the PVL prophage, one isolate carried another prophage and the second isolate carried two additional prophages, whereby the region between these two prophages was inverted. This "flipped" region comprised about 1,083,000 bp, or more than a third of the strain's genome, and it included the PVL prophage. Prophages were induced by Mitomycin C treatment and subjected to transmission electron microscopy (TEM). This yielded, in accordance to the sequencing results, one or, respectively, two distinct populations of icosahedral phages. It also showed prolate phages which presumptively might be identified as the PVL phage. This observation highlights the significance bacteriophages have as agents of horizontal gene transfer as well as the need for monitoring emerging staphylococcal strains, especially in cosmopolitan settings such as the UAE.

Assessing phage-host population dynamics by reintroducing virulent viruses to synthetic microbiomes

Microbiomes feature complex interactions between diverse bacteria and bacteriophages. Synthetic microbiomes offer a powerful way to study these interactions; however, a major challenge is obtaining a representative bacteriophage population during the bacterial isolation process. We demonstrate that colony isolation reliably excludes virulent viruses from sample sources with low virion-to-bacteria ratios such as feces, creating "virulent virus-free" controls. When the virulent dsDNA virome is reintroduced to a 73-strain synthetic gut microbiome in a bioreactor model of the human colon, virulent viruses target susceptible strains without significantly altering community structure or metabolism. In addition, we detected signals of prophage induction that associate with virulent predation. Overall, our findings indicate that dilution-based isolation methods generate synthetic gut microbiomes that are heavily depleted, if not devoid, of virulent viruses and that such viruses, if reintroduced, have a targeted effect on community assembly, metabolism, and prophage replication.

Viruses of Nitrogen-Fixing Mesorhizobium Bacteria in Globally Distributed Chickpea Root Nodules

Legume nodules are specialized environments on plant roots that are induced and dominated by nitrogen-fixing bacteria. Bacteriophages (phages) in these nodules could potentially provide top-down controls on the population size and, therefore, the function of nitrogen-fixing symbionts. Here we sought to characterize the diversity and biogeographical patterns of phages that infect nitrogen-fixing Mesorhizobium symbionts isolated from root nodules, leveraging 266 genomes of Mesorhizobium isolated from nodules and 648 nodule metagenomes collected from three species of chickpea plants ( Cicer spp.) under different agricultural management practices, spanning eight countries on five continents. We identified 106 phage populations (viral operational taxonomic units [vOTUs]) in Mesorhizobium draft genomes, 37% of which were confirmed as likely prophages. These vOTUs were detected in 64% of the Mesorhizobium-dominated nodule metagenomes and 58% of the Mesorhizobium isolates. Per metagenome, 1 to 16 putative Mesorhizobium vOTUs were detected, with more than half of the nodules containing only one such vOTU. The majority of vOTUs were detected exclusively in Ethiopia, followed by India and Morocco, with the lowest richness of putative Mesorhizobium phages in countries that applied industrial Mesorhizobium inoculants to crops. Two vOTUs were identified in five or more countries and in nodules dominated by different strains of Mesorhizobium, suggesting infection of diverse Mesorhizobium hosts and long-term interactions. Beta-diversity of these Mesorhizobium phage assemblages was significantly correlated with the dominant Mesorhizobium strain, but not with measured environmental parameters. Our findings indicate that nitrogen-fixing nodules in chickpea plants can contain distinct viral assemblages, with potential impacts on the nodule microbiome that bear further exploration.

Bacteriophage specificity is impacted by interactions between bacteria.

There is significant natural diversity in how many different types of bacteria a bacteriophage can infect, but the mechanisms driving this diversity are unclear. This study uses a combination of mathematical modeling and an in vitro system consisting of Escherichia coli, Salmonella enterica, a T5-like generalist phage, and the specialist phage P22vir to highlight the connection between bacteriophage specificity and interactions between their potential microbial prey. Mathematical modeling suggests that competing bacteria tend to favor generalist bacteriophage, while bacteria that benefit each other tend to favor specialist bacteriophage. Experimental results support this general finding. The experiments also show that the optimal phage strategy is impacted by phage degradation and bacterial physiology. These findings enhance our understanding of how complex microbial communities shape selection on bacteriophage specificity, which may improve our ability to use phage to manage antibiotic-resistant microbial infections.