Phage-bacteria infection networks

Application of machine learning in bacteriophage research

Phages are one of the key components in the structure, dynamics, and interactions of microbial communities in different bins. It has a clear impact on human health and the food industry. Bacteriophage characterization using in vitro approaches are time/cost consuming and laborious tasks. On the other hand, with the advent of new high-throughput sequencing technology, the development of a powerful computational framework to characterize the newly identified bacteriophages is inevitable for future research. Machine learning includes powerful techniques that enable the analysis of complex datasets for knowledge discovery and pattern recognition. In this study, we have conducted a comprehensive review of machine learning methods application using different types of features were applied in various aspects of bacteriophage research including, automated curation, identification, classification, host species recognition, virion protein identification, and life cycle prediction. Moreover, potential limitations and advantages of the developed frameworks were discussed.

Plant virus evolution under strong drought conditions results in a transition from parasitism to mutualism

Environmental conditions are an important factor driving pathogens' evolution. Here, we explore the effects of drought stress in plant virus evolution. We evolved turnip mosaic potyvirus in well-watered and drought conditions in Arabidopsis thaliana accessions that differ in their response to virus infection. Virus adaptation occurred in all accessions independently of watering status. Drought-evolved viruses conferred a significantly higher drought tolerance to infected plants. By contrast, nonsignificant increases in tolerance were observed in plants infected with viruses evolved under standard watering. The magnitude of this effect was dependent on the plant accessions. Differences in tolerance were correlated to alterations in the expression of host genes, some involved in regulation of the circadian clock, as well as in deep changes in the balance of phytohormones regulating defense and growth signaling pathways. Our results show that viruses can promote host survival in situations of abiotic stress, with the magnitude of such benefit being a selectable trait.

Characteristics on host specificity, infection, and temperature stability of Weissella phages from watery kimchi

Three bacteriophages (phage), ΦWC53, ΦWC54, and ΦWC56, of Weissella were isolated from watery kimchi and characterized. ΦWC53 belonged to Siphoviridae and ΦWC54 and ΦWC56 belonged to Myoviridae family. By one-step growth, the burst sizes were 5-260 particles/infected cells and the latent periods were 20-45 min. The phages infected Weissella spp., Leuconostoc mesenteroides, and Lactobacillus spp., differently by showing clear or turbid plaques. The phage adsorption rates on lactic acid bacteria were high on Weissella and low on Leuconostoc and Lactobacillus. However, the adsorption of ΦWC53 occurred variously among Weissella spp. and Weissella host grew well in the liquid culture without lysis after challenging by ΦWC53. Tolerances of these phages to temperature showed more various than those to pH. ΦWC53 was stable at 7 °C and 30 °C, but ΦWC54 and ΦWC56 were stable only at 7 °C. Therefore, three Weissella phages belonged to the different families and indicated diverse infection patterns on Weissella, Leuconostoc, and Lactobacillus with various stabilities for pH and temperature.

Campylobacter phages use hypermutable polyG tracts to create phenotypic diversity and evade bacterial resistance

Phase variation is a common mechanism for creating phenotypic heterogeneity of surface structures in bacteria important for niche adaptation. In Campylobacter, phase variation occurs by random variation in hypermutable homonucleotide 7-11 G (polyG) tracts. To elucidate how phages adapt to phase-variable hosts, we study Fletchervirus phages infecting Campylobacter dependent on a phase-variable receptor. Our data demonstrate that Fletcherviruses mimic their host and encode hypermutable polyG tracts, leading to phase-variable expression of two of four receptor-binding proteins. This creates phenotypically diverse phage populations, including a sub-population that infects the bacterial host when the phase-variable receptor is not expressed. Such population dynamics of both phage and host promote co-existence in a shared niche. Strikingly, we identify polyG tracts in more than 100 phage genera, infecting more than 70 bacterial species. Future experimental work may confirm phase variation as a widespread strategy for creating phenotypically diverse phage populations.

Genomic Analysis of a Novel Phage Infecting the Turkey Pathogen Escherichia coli APEC O78 and Its Endolysin Activity

Due to the increasing spread of multidrug-resistant (MDR) bacteria, phage therapy is considered one of the most promising methods for addressing MDR bacteria. Escherichia coli lives symbiotically in the intestines of humans and some animals, and most strains are beneficial in terms of maintaining a healthy digestive tract. However, some E. coli strains can cause serious zoonotic diseases, including diarrhea, pneumonia, urinary tract infections, and hemolytic uremic syndrome. In this study, we characterized a newly isolated Myoviridae phage, vB_EcoM_APEC. The phage vB_EcoM_APEC was able to infect E. coli APEC O78, which is the most common MDR E. coli serotype in turkeys. Additionally, the phage's host range included Klebsiella pneumoniae and other E. coli strains. The genome of phage vB_EcoM_APEC (GenBank accession number MT664721) was 35,832 bp in length, with 52 putative open reading frames (ORFs) and a GC content of 41.3%. The genome of vB_EcoM_APEC exhibited low similarity (79.1% identity and 4.0% coverage) to the genome of Acinetobacter phage vB_AbaM_IME284 (GenBank no. MH853787.1) according to the nucleotide Basic Local Alignment Search Tool (BLASTn). Phylogenetic analysis revealed that vB_EcoM_APEC was a novel phage, and its genome sequence showed low similarity to other available phage genomes. Gene annotation indicated that the protein encoded by orf11 was an endolysin designated as LysO78, which exhibited 64.7% identity (91.0% coverage) with the putative endolysin of Acinetobacter baumannii phage vB_AbaM_B9. The LysO78 protein belongs to glycoside hydrolase family 19, and was described as being a chitinase class I protein. LysO78 is a helical protein with 12 α-helices containing a large domain and a small domain in terms of the predicted three-dimensional structure. The results of site-directed mutagenesis indicated that LysO78 contained the catalytic residues E54 and E64. The purified endolysin exhibited broad-spectrum bacteriolytic activity against Gram-negative strains, including the genera Klebsiella, Salmonella, Shigella, Burkholderia, Yersinia, and Pseudomonas, as well as the species Chitinimonas arctica, E. coli, Ralstonia solanacearum, and A. baumannii. An enzymatic assay showed that LysO78 had highly lytic peptidoglycan hydrolases activity (64,620,000 units/mg) against E. coli APEC O78, and that LysO78 had lytic activity in the temperature range of 4-85 °C, with an optimal temperature of 28 °C and optimal pH of 8.0, and was active at pH 3.0-12.0. Overall, the results suggested that LysO78 might be a promising therapeutic agent for controlling MDR E. coli APEC O78 and nosocomial infections caused by multidrug-resistant bacteria.

Characteristics of Bacteriophage Isolates and Expression of Shiga Toxin Genes Transferred to Non Shiga Toxin-Producing E. coli by Transduction

A risk analysis of Shiga toxin (Stx)-encoding bacteriophage was carried out by confirming the transduction phage to non-Stx-producing Escherichia coli (STEC) and subsequent expression of the Shiga toxin genes. The virulence factor stx1 was identified in five phages, and both stx1 and stx2 were found in four phages from a total of 19 phage isolates with seven non-O157 STEC strains. The four phages, designated as φNOEC41, φNOEC46, φNOEC47, and φNOEC49, belonged morphologically to the Myoviridae family. The stabilities of these phages to temperature, pH, ethanol, and NaClO were high with some variabilities among the phages. The infection of five non-STEC strains by nine Stx-encoding phages occurred at a rate of approximately 40%. Non-STEC strains were transduced by Stx-encoding phage to become lysogenic strains, and seven convertant strains had stx1 and/or stx2 genes. Only the stx1 gene was transferred to the receptor strains without any deletion. Gene expression of a convertant having both stx1 and stx2 genes was confirmed to be up to 32 times higher for Stx1 in 6% NaCl osmotic media and twice for Stx2 in 4% NaCl media, compared with expression in low-salt environments. Therefore, a new risk might arise from the transfer of pathogenic genes from Stx-encoding phages to otherwise harmless hosts. Without adequate sterilization of food exposed to various environments, there is a possibility that the toxicity of the phages might increase.

Sewage and sewage-contaminated environments are the most prominent sources to isolate phages against Pseudomonas aeruginosa

Background

P. aeruginosa is the primary source of hospital-acquired infections. Unfortunately, antibiotic resistance is growing to precariously high levels, making the infections by this pathogen life-threatening and hard to cure. One possible alternative to antibiotics is to use phages. However, the isolation of phages suitable for phage therapy— be lytic, be efficient, and have a broad host range —against some target bacteria has proven difficult. To identify the best places to look for these phages against P. aeruginosa we screened hospital sewages, soils, and rivers in two cities.

Results

We isolated eighteen different phages, determined their host range, infection property, and plaque morphology. We found that the sewage and sewage-contaminated environments are the most reliable sources for the isolation of Pseudomonas phages. In addition, phages isolated from hospital sewage showed the highest efficiency in lysing the bacteria used for host range determination. In contrast, phages from the river had larger plaque size and lysed bacteria with higher levels of antibiotic resistance.

Conclusions

Our findings provided additional support for the importance of sewage as the source of phage isolation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02197-z.

A New Pipeline for Designing Phage Cocktails Based on Phage-Bacteria Infection Networks

In recent years, the spread of antibiotic-resistant bacteria and efforts to preserve food microbiota have induced renewed interest in phage therapy. Phage cocktails, instead of a single phage, are commonly used as antibacterial agents since the hosts are unlikely to become resistant to several phages simultaneously. While the spectrum of activity might increase with cocktail complexity, excessive phages could produce side effects, such as the horizontal transfer of genes that augment the fitness of host strains, dysbiosis or high manufacturing costs. Therefore, cocktail formulation represents a compromise between achieving substantial reduction in the bacterial loads and restricting its complexity. Despite the abovementioned points, the observed bacterial load reduction does not increase significantly with the size of phage cocktails, indicating the requirement for a systematic approach to their design. In this work, the information provided by host range matrices was analyzed after building phage-bacteria infection networks (PBINs). To this end, we conducted a meta-analysis of 35 host range matrices, including recently published studies and new datasets comprising Escherichia coli strains isolated during ripening of artisanal raw milk cheese and virulent coliphages from ewes' feces. The nestedness temperature, which reflects the host range hierarchy of the phages, was determined from bipartite host range matrices using heuristic (Nestedness Temperature Calculator) and genetic (BinMatNest) algorithms. The latter optimizes matrix packing, leading to lower temperatures, i.e., it simplifies the identification of the phages with the broadest host range. The structure of infection networks suggests that generalist phages (and not specialist phages) tend to succeed in infecting less susceptible bacteria. A new metric (Φ), which considers some properties of the host range matrices (fill, temperature, and number of bacteria), is proposed as an estimator of phage cocktail size. To identify the best candidates, agglomerative hierarchical clustering using Ward's method was implemented. Finally, a cocktail was formulated for the biocontrol of cheese-isolated E. coli, reducing bacterial counts by five orders of magnitude.

A Spontaneous rapZ Mutant Impairs Infectivity of Lytic Bacteriophage vB_EcoM_JS09 against Enterotoxigenic Escherichia coli

Our understanding of the mechanisms underlying phage-bacterium interactions remains limited. In Escherichia coli, RapZ regulates glucosamine-6-phosphate (GlcN6P) metabolism, the formation of which initiates synthesis of the bacterial cell envelope, including lipopolysaccharides (LPS). However, the role of RapZ, if any, on phage infectivity remains to be investigated. Here, we isolated strains of enterotoxigenic E. coli (ETEC) resistant to its specific lytic bacteriophage vB_EcoM_JS09 (JS09) in a phage aerosol spray experiment. Whole-genome analysis of phage-resistant bacteria revealed the rapZ gene acquired a premature stop mutation at amino acid 227. Here, we report that the mutation in the rapZ gene confers resistance by inhibiting 93.5% phage adsorption. Furthermore, this mutation changes the morphology of phage plaques, reduces efficiency of plating and phage propagation efficiency, and impairs the infectivity of phage JS09 against ETEC. Using scanning electron microscopy assays, we attribute the inability of the phage to adsorb to the loss of receptors in strains with defective RapZ. Analysis of the LPS profile shows that strains with defective RapZ inhibit phage infection by changing the LPS profile in E. coli. Preincubation of phage JS09 with LPS extracted from a wild-type (WT) strain blocked infection, suggesting LPS is the host receptor for phage JS09 adsorption. Our data uncover the mechanism by which ETEC resists infection of phage JS09 by mutating the rapZ gene and then increasing the expression of glmS and changing the phage receptor-LPS profile. These findings provide insight into the function of the rapZ gene for efficient infection of phage JS09.

IMPORTANCE The development of phage-resistant bacteria is a challenging problem for phage therapy. However, our knowledge of phage resistance mechanisms is still limited. RapZ is an RNase adaptor protein encoded by the rapZ gene and plays an important function in Gram-positive and Gram-negative bacteria. Here, we report the whole-genome analysis of a phage-resistant enterotoxigenic Escherichia coli (ETEC) strain, which revealed that the rapZ gene acquired a premature stop mutation (E227Stop). We show that the premature stop mutation of rapZ impairs the infectivity of phage JS09 in ETEC. Furthermore, our findings indicate that ETEC becomes resistant against the adsorption and infection of phage JS09 by mutating the rapZ gene, increasing the expression of glmS, and changing the phage receptor-LPS profile. It is also first reported here that RapZ is essential for efficient infection of phage JS09.

Balance between promiscuity and specificity in phage λ host range

As hosts acquire resistance to viruses, viruses must overcome that resistance to re-establish infectivity, or go extinct. Despite the significant hurdles associated with adapting to a resistant host, viruses are evolutionarily successful and maintain stable coevolutionary relationships with their hosts. To investigate the factors underlying how pathogens adapt to their hosts, we performed a deep mutational scan of the region of the λ tail fiber tip protein that mediates contact with the receptor on λ's host, Escherichia coli. Phages harboring amino acid substitutions were subjected to selection for infectivity on wild type E. coli, revealing a highly restrictive fitness landscape, in which most substitutions completely abrogate function. A subset of positions that are tolerant of mutation in this assay, but diverse over evolutionary time, are associated with host range expansion. Imposing selection for phage infectivity on three λ-resistant hosts, each harboring a different missense mutation in the λ receptor, reveals hundreds of adaptive variants in λ. We distinguish λ variants that confer promiscuity, a general ability to overcome host resistance, from those that drive host-specific infectivity. Both processes may be important in driving adaptation to a novel host.

Predicting bacteriophage hosts based on sequences of annotated receptor-binding proteins

Nowadays, bacteriophages are increasingly considered as an alternative treatment for a variety of bacterial infections in cases where classical antibiotics have become ineffective. However, characterizing the host specificity of phages remains a labor- and time-intensive process. In order to alleviate this burden, we have developed a new machine-learning-based pipeline to predict bacteriophage hosts based on annotated receptor-binding protein (RBP) sequence data. We focus on predicting bacterial hosts from the ESKAPE group, Escherichia coli, Salmonella enterica and Clostridium difficile. We compare the performance of our predictive model with that of the widely used Basic Local Alignment Search Tool (BLAST). Our best-performing predictive model reaches Precision-Recall Area Under the Curve (PR-AUC) scores between 73.6 and 93.8% for different levels of sequence similarity in the collected data. Our model reaches a performance comparable to that of BLASTp when sequence similarity in the data is high and starts outperforming BLASTp when sequence similarity drops below 75%. Therefore, our machine learning methods can be especially useful in settings in which sequence similarity to other known sequences is low. Predicting the hosts of novel metagenomic RBP sequences could extend our toolbox to tune the host spectrum of phages or phage tail-like bacteriocins by swapping RBPs.

Sustainable Microbiome: a symphony orchestrated by synthetic phages

We are surrounded by microbes, mostly bacteria and their viruses or phages, on the inside and outside of our bodies. These bacteria in constant interactions with phages are regulating multiple functions critical to our health. Luckily, they are amenable, but we need precise tools for their safe manipulation and improving human health. Here, we argue that recent advances in single-cell technologies, culturomics and synthetic biology offer exciting opportunities to create these tools as well as revealing specific phages-bacteria interactions in the body.

The network structure and eco-evolutionary dynamics of CRISPR-induced immune diversification

As a heritable sequence-specific adaptive immune system, CRISPR-Cas is a powerful force shaping strain diversity in host-virus systems. While the diversity of CRISPR alleles has been explored, the associated structure and dynamics of host-virus interactions have not. We explore the role of CRISPR in mediating the interplay between host-virus interaction structure and eco-evolutionary dynamics in a computational model and compare the results with three empirical datasets from natural systems. We show that the structure of the networks describing who infects whom and the degree to which strains are immune, are respectively modular (containing groups of hosts and viruses that interact strongly) and weighted-nested (specialist hosts are more susceptible to subsets of viruses that in turn also infect the more generalist hosts with many spacers matching many viruses). The dynamic interplay between these networks influences transitions between dynamical regimes of virus diversification and host control. The three empirical systems exhibit weighted-nested immunity networks, a pattern our theory shows is indicative of hosts able to suppress virus diversification. Previously missing from studies of microbial host-pathogen systems, the immunity network plays a key role in the coevolutionary dynamics.

Diversity and Local Coadaptation of Escherichia coli and Coliphages From Small Ruminants

Bacteriophages are highly specific predators that drive bacterial diversity through coevolution while striking tradeoffs among preserving host populations for long-term exploitation and increasing their virulence, structural stability, or host range. Escherichia coli and other coliform bacteria present in the microbiota of milk and during early ripening of raw milk cheeses have been linked to the production of gas, manifested by the appearance of eyes, and the development of off-flavors; thus, they might cause early blowing and cheese spoilage. Here, we report the characterization of coliphages isolated from manure from small ruminant farms and E. coli strains isolated from goat and sheep raw milk cheese. Additionally, the virulence and host range of locally isolated and laboratory collection phages were determined by comparing the susceptibility of E. coli strains from different sources. In agreement with the high genetic diversity found within the species E. coli, clustering analysis of whole-cell protein revealed a total of 13 distinct profiles but none of the raw milk cheese isolates showed inhibition of growth by reference or water-isolated coliphages. Conversely, 10 newly isolated phages had a broad host range (i.e., able to lyse ≥50% of bacterial hosts tested), thus exhibiting utility for biocontrol and only one cheese-isolated E. coli strain was resistant to all the phages. Whereas there was a high positive correlation between bacterial susceptibility range and lysis intensity, the phages virulence decreased as range increased until reaching a plateau. These results suggest local gene-for-gene coevolution between hosts and phages with selective tradeoffs for both resistance and competitive ability of the bacteria and host-range extension and virulence of the phage populations. Hence, different phage cocktail formulations might be required when devising long-term and short-term biocontrol strategies.

High-throughput mapping of the phage resistance landscape in E. coli

Bacteriophages (phages) are critical players in the dynamics and function of microbial communities and drive processes as diverse as global biogeochemical cycles and human health. Phages tend to be predators finely tuned to attack specific hosts, even down to the strain level, which in turn defend themselves using an array of mechanisms. However, to date, efforts to rapidly and comprehensively identify bacterial host factors important in phage infection and resistance have yet to be fully realized. Here, we globally map the host genetic determinants involved in resistance to 14 phylogenetically diverse double-stranded DNA phages using two model Escherichia coli strains (K-12 and BL21) with known sequence divergence to demonstrate strain-specific differences. Using genome-wide loss-of-function and gain-of-function genetic technologies, we are able to confirm previously described phage receptors as well as uncover a number of previously unknown host factors that confer resistance to one or more of these phages. We uncover differences in resistance factors that strongly align with the susceptibility of K-12 and BL21 to specific phage. We also identify both phage-specific mechanisms, such as the unexpected role of cyclic-di-GMP in host sensitivity to phage N4, and more generic defenses, such as the overproduction of colanic acid capsular polysaccharide that defends against a wide array of phages. Our results indicate that host responses to phages can occur via diverse cellular mechanisms. Our systematic and high-throughput genetic workflow to characterize phage-host interaction determinants can be extended to diverse bacteria to generate datasets that allow predictive models of how phage-mediated selection will shape bacterial phenotype and evolution. The results of this study and future efforts to map the phage resistance landscape will lead to new insights into the coevolution of hosts and their phage, which can ultimately be used to design better phage therapeutic treatments and tools for precision microbiome engineering.

A methodological framework to analyse determinants of host-microbiota networks, with an application to the relationships between Daphnia magna's gut microbiota and bacterioplankton

The past 30 years have seen both a surge of interest in assessing ecological interactions using tools borrowed from network theory and an explosion of data on the occurrence of microbial symbionts thanks to next-generation sequencing. Given that classic network methods cannot currently measure the respective effects of different environmental and biological drivers on network structure, we here present two methods to elucidate the determinants of bipartite interaction networks. The first method is based on classifications and compares communities within networks to the grouping of nodes by treatment or similar controlling groups. The second method assesses the link between multivariate explanatory variables and network structure using redundancy analyses after singular value decomposition. In both methods, the significance of effects can be gauged through two randomizations. Our methods were applied to experimental data on Daphnia magna and its interactions with gut microbiota and bacterioplankton. The whole network was affected by Daphnia's diet (algae and/or cyanobacteria) and sample type, but not by Daphnia genotype. At coarse grains, bacterioplankton and gut microbiota communities were different. At this scale, the structure of the gut microbiota-based network was not linked to any explanatory factors, while the bacterioplankton-based network was related to both Daphnia's diet and genotype. At finer grains, Daphnia's diet and genotype affected both microbial networks, but the effect of diet on gut microbiota network structure was mediated solely by differences in microbial richness. While no reciprocal effect between the microbial communities could be found, fine-grained analyses presented a more nuanced picture, with bacterioplankton likely affecting the composition of the gut microbiota. Our methods are widely applicable to bipartite networks, can elucidate both controlled and environmental effects in experimental setting using a large amount of sequencing data and can tease apart reciprocal effects of networks on one another. The twofold approach we propose has the advantage of being able to tease apart effects at different scales of network structure, thus allowing for detailed assessment of reciprocal effects of linked networks on one another. As such, our network methods can help ecologists understand huge datasets reporting microbial co-occurrences within different hosts.

Modular prophage interactions driven by capsule serotype select for capsule loss under phage predation

Klebsiella species are able to colonize a wide range of environments and include worrisome nosocomial pathogens. Here, we sought to determine the abundance and infectivity of prophages of Klebsiella to understand how the interactions between induced prophages and bacteria affect population dynamics and evolution. We identified many prophages in the species, placing these taxa among the top 5% of the most polylysogenic bacteria. We selected 35 representative strains of the Klebsiella pneumoniae species complex to establish a network of induced phage-bacteria interactions. This revealed that many prophages are able to enter the lytic cycle, and subsequently kill or lysogenize closely related Klebsiella strains. Although 60% of the tested strains could produce phages that infect at least one other strain, the interaction network of all pairwise cross-infections is very sparse and mostly organized in modules corresponding to the strains' capsule serotypes. Accordingly, capsule mutants remain uninfected showing that the capsule is a key factor for successful infections. Surprisingly, experiments in which bacteria are predated by their own prophages result in accelerated loss of the capsule. Our results show that phage infectiousness defines interaction modules between small subsets of phages and bacteria in function of capsule serotype. This limits the role of prophages as competitive weapons because they can infect very few strains of the species complex. This should also restrict phage-driven gene flow across the species. Finally, the accelerated loss of the capsule in bacteria being predated by their own phages, suggests that phages drive serotype switch in nature.

Biofilm Structure Promotes Coexistence of Phage-Resistant and Phage-Susceptible Bacteria

Encounters among bacteria and their viral predators (bacteriophages) are among the most common ecological interactions on Earth. These encounters are likely to occur with regularity inside surface-bound communities that microbes most often occupy in natural environments. Such communities, termed biofilms, are spatially constrained: interactions become limited to near neighbors, diffusion of solutes and particulates can be reduced, and there is pronounced heterogeneity in nutrient access and physiological state. It is appreciated from prior theoretical work that phage-bacteria interactions are fundamentally different in spatially structured contexts, as opposed to well-mixed liquid culture. Spatially structured communities are predicted to promote the protection of susceptible host cells from phage exposure, and thus weaken selection for phage resistance. The details and generality of this prediction in realistic biofilm environments, however, are not known. Here, we explore phage-host interactions using experiments and simulations that are tuned to represent the essential elements of biofilm communities. Our simulations show that in biofilms, phage-resistant cells-as their relative abundance increases-can protect clusters of susceptible cells from phage exposure, promoting the coexistence of susceptible and phage-resistant bacteria under a large array of conditions. We characterize the population dynamics underlying this coexistence, and we show that coexistence is recapitulated in an experimental model of biofilm growth measured with confocal microscopy. Our results provide a clear view into the dynamics of phage resistance in biofilms with single-cell resolution of the underlying cell-virion interactions, linking the predictions of canonical theory to realistic models and in vitro experiments of biofilm growth.IMPORTANCE In the natural environment, bacteria most often live in communities bound to one another by secreted adhesives. These communities, or biofilms, play a central role in biogeochemical cycling, microbiome functioning, wastewater treatment, and disease. Wherever there are bacteria, there are also viruses that attack them, called phages. Interactions between bacteria and phages are likely to occur ubiquitously in biofilms. We show here, using simulations and experiments, that biofilms will in most conditions allow phage-susceptible bacteria to be protected from phage exposure, if they are growing alongside other cells that are phage resistant. This result has implications for the fundamental ecology of phage-bacteria interactions, as well as the development of phage-based antimicrobial therapeutics.

Initial Mapping of the New York City Wastewater Virome

Wastewater is a rich source of microbial life and contains bacteria, viruses, and other microbes found in human waste as well as environmental runoff sources. As part of an effort to characterize the New York City wastewater metagenome, we profiled the viral community of sewage samples across all five boroughs of NYC and found that local sampling sites have unique sets of viruses. We focused on bacteriophages, or viruses of bacteria, to understand how they may influence the microbial ecology of this system. We identified several new clusters of phages and successfully associated them with bacterial hosts, providing insight into virus-host interactions in urban wastewater. This study provides a first look into the viral communities present across the wastewater system in NYC and points to their functional importance in this environment.

Bacteriophages are abundant members of all microbiomes studied to date, influencing microbial communities through interactions with their bacterial hosts. Despite their functional importance and ubiquity, phages have been underexplored in urban environments compared to their bacterial counterparts. We profiled the viral communities in New York City (NYC) wastewater using metagenomic data collected in November 2014 from 14 wastewater treatment plants. We show that phages accounted for the largest viral component of the sewage samples and that specific virus communities were associated with local environmental conditions within boroughs. The vast majority of the virus sequences had no homology matches in public databases, forming an average of 1,700 unique virus clusters (putative genera). These new clusters contribute to elucidating the overwhelming proportion of data that frequently goes unidentified in viral metagenomic studies. We assigned potential hosts to these phages, which appear to infect a wide range of bacterial genera, often outside their presumed host. We determined that infection networks form a modular-nested pattern, indicating that phages include a range of host specificities, from generalists to specialists, with most interactions organized into distinct groups. We identified genes in viral contigs involved in carbon and sulfur cycling, suggesting functional importance of viruses in circulating pathways and gene functions in the wastewater environment. In addition, we identified virophage genes as well as a nearly complete novel virophage genome. These findings provide an understanding of phage abundance and diversity in NYC wastewater, previously uncharacterized, and further examine geographic patterns of phage-host association in urban environments.

IMPORTANCE Wastewater is a rich source of microbial life and contains bacteria, viruses, and other microbes found in human waste as well as environmental runoff sources. As part of an effort to characterize the New York City wastewater metagenome, we profiled the viral community of sewage samples across all five boroughs of NYC and found that local sampling sites have unique sets of viruses. We focused on bacteriophages, or viruses of bacteria, to understand how they may influence the microbial ecology of this system. We identified several new clusters of phages and successfully associated them with bacterial hosts, providing insight into virus-host interactions in urban wastewater. This study provides a first look into the viral communities present across the wastewater system in NYC and points to their functional importance in this environment.

Stabilization of extensive fine-scale diversity by ecologically driven spatiotemporal chaos

It has recently become apparent that the diversity of microbial life extends far below the species level to the finest scales of genetic differences. Remarkably, extensive fine-scale diversity can coexist spatially. How is this diversity stable on long timescales, despite selective or ecological differences and other evolutionary processes? Most work has focused on stable coexistence or assumed ecological neutrality. We present an alternative: extensive diversity maintained by ecologically driven spatiotemporal chaos, with no assumptions about niches or other specialist differences between strains. We study generalized Lotka-Volterra models with antisymmetric correlations in the interactions inspired by multiple pathogen strains infecting multiple host strains. Generally, these exhibit chaos with increasingly wild population fluctuations driving extinctions. But the simplest spatial structure, many identical islands with migration between them, stabilizes a diverse chaotic state. Some strains (subspecies) go globally extinct, but many persist for times exponentially long in the number of islands. All persistent strains have episodic local blooms to high abundance, crucial for their persistence as, for many, their average population growth rate is negative. Snapshots of the abundance distribution show a power law at intermediate abundances that is essentially indistinguishable from the neutral theory of ecology. But the dynamics of the large populations are much faster than birth-death fluctuations. We argue that this spatiotemporally chaotic "phase" should exist in a wide range of models, and that even in rapidly mixed systems, longer-lived spores could similarly stabilize a diverse chaotic phase.

Mapping co-ancestry connections between the genome of a Medieval individual and modern Europeans

Historical genetic links among similar populations can be difficult to establish. Identity by descent (IBD) analyses find genomic blocks that represent direct genealogical relationships among individuals. However, this method has rarely been applied to ancient genomes because IBD stretches are progressively fragmented by recombination and thus not recognizable after few tens of generations. To explore such genealogical relationships, we estimated long IBD blocks among modern Europeans, generating networks to uncover the genetic structures. We found that Basques, Sardinians, Icelanders and Orcadians form, each of them, highly intraconnected sub-clusters in a European network, indicating dense genealogical links within small, isolated populations. We also exposed individual genealogical links -such as the connection between one Basque and one Icelandic individual- that cannot be uncovered with other, widely used population genetics methods such as PCA or ADMIXTURE. Moreover, using ancient DNA technology we sequenced a Late Medieval individual (Barcelona, Spain) to high genomic coverage and identified IBD blocks shared between her and modern Europeans. The Medieval IBD blocks are statistically overrepresented only in modern Spaniards, which is the geographically closest population. This approach can be used to produce a fine-scale reflection of shared ancestry across different populations of the world, offering a direct genetic link from the past to the present.

Impact of phages on soil bacterial communities and nitrogen availability under different assembly scenarios

BACKGROUND

Bacteriophages, the viruses infecting bacteria, are biological entities that can control their host populations. The ecological relevance of phages for microbial systems has been widely explored in aquatic environments, but the current understanding of the role of phages in terrestrial ecosystems remains limited. Here, our objective was to quantify the extent to which phages drive the assembly and functioning of soil bacterial communities. We performed a reciprocal transplant experiment using natural and sterilized soil incubated with different combinations of two soil microbial communities, challenged against native and non-native phage suspensions as well as against a cocktail of phage isolates. We tested three different community assembly scenarios by adding phages: (a) during soil colonization, (b) after colonization, and (c) in natural soil communities. One month after inoculation with phage suspensions, bacterial communities were assessed by 16S rRNA amplicon gene sequencing.

RESULTS

By comparing the treatments inoculated with active versus autoclaved phages, our results show that changes in phage pressure have the potential to impact soil bacterial community composition and diversity. We also found a positive effect of active phages on the soil ammonium concentration in a few treatments, which indicates that increased phage pressure may also be important for soil functions.

CONCLUSIONS

Overall, the present work contributes to expand the current knowledge about soil phages and provide some empirical evidence supporting their relevance for soil bacterial community assembly and functioning. Video Abstract.

Ecological Assembly Processes Are Coordinated between Bacterial and Viral Communities in Fractured Shale Ecosystems

Interactions between viral communities and their microbial hosts have been the subject of many recent studies in a wide range of ecosystems. The degree of coordination between ecological assembly processes influencing viral and microbial communities, however, has been explored to a much lesser degree. By using a combined null modeling approach, this study investigated the ecological assembly processes influencing both viral and microbial community structure within hydraulically fractured shale environments. Among other results, significant relationships between the structuring processes affecting both the viral and microbial community were observed, indicating that ecological assembly might be coordinated between these communities despite differing selective pressures. Within this deep subsurface ecosystem, these results reveal a potentially important balance of ecological dynamics that must be maintained to enable long-term microbial community persistence. More broadly, this relationship begins to provide insight into the development of communities across trophic levels.

The ecological drivers that concurrently act upon both a virus and its host and that drive community assembly are poorly understood despite known interactions between viral populations and their microbial hosts. Hydraulically fractured shale environments provide access to a closed ecosystem in the deep subsurface where constrained microbial and viral community assembly processes can be examined. Here, we used metagenomic analyses of time-resolved-produced fluid samples from two wells in the Appalachian Basin to track viral and host dynamics and to investigate community assembly processes. Hypersaline conditions within these ecosystems should drive microbial community structure to a similar configuration through time in response to common osmotic stress. However, viral predation appears to counterbalance this potentially strong homogeneous selection and pushes the microbial community toward undominated assembly. In comparison, while the viral community was also influenced by substantial undominated processes, it assembled, in part, due to homogeneous selection. When the overall assembly processes acting upon both these communities were directly compared with each other, a significant relationship was revealed, suggesting an association between microbial and viral community development despite differing selective pressures. These results reveal a potentially important balance of ecological dynamics that must be in maintained within this deep subsurface ecosystem in order for the microbial community to persist over extended time periods. More broadly, this relationship begins to provide knowledge underlying metacommunity development across trophic levels.

IMPORTANCE Interactions between viral communities and their microbial hosts have been the subject of many recent studies in a wide range of ecosystems. The degree of coordination between ecological assembly processes influencing viral and microbial communities, however, has been explored to a much lesser degree. By using a combined null modeling approach, this study investigated the ecological assembly processes influencing both viral and microbial community structure within hydraulically fractured shale environments. Among other results, significant relationships between the structuring processes affecting both the viral and microbial community were observed, indicating that ecological assembly might be coordinated between these communities despite differing selective pressures. Within this deep subsurface ecosystem, these results reveal a potentially important balance of ecological dynamics that must be maintained to enable long-term microbial community persistence. More broadly, this relationship begins to provide insight into the development of communities across trophic levels.

Coexistence of nestedness and modularity in host-pathogen infection networks

The long-term coevolution of hosts and pathogens in their environment forms a complex web of multi-scale interactions. Understanding how environmental heterogeneity affects the structure of host-pathogen networks is a prerequisite for predicting disease dynamics and emergence. Although nestedness is common in ecological networks, and theory suggests that nested ecosystems are less prone to dynamic instability, why nestedness varies in time and space is not fully understood. Many studies have been limited by a focus on single habitats and the absence of a link between spatial variation and structural heterogeneity such as nestedness and modularity. Here we propose a neutral model for the evolution of host-pathogen networks in multiple habitats. In contrast to previous studies, our study proposes that local modularity can coexist with global nestedness, and shows that real ecosystems are found in a continuum between nested-modular and nested networks driven by intraspecific competition. Nestedness depends on neutral mechanisms of community assembly, whereas modularity is contingent on local adaptation and competition. The structural pattern may change spatially and temporally but remains stable over evolutionary timescales. We validate our theoretical predictions with a longitudinal study of plant-virus interactions in a heterogeneous agricultural landscape.

The Role of Evolution in Shaping Ecological Networks

The structure of ecological networks reflects the evolutionary history of their biotic components, and their dynamics are strongly driven by ecoevolutionary processes. Here, we present an appraisal of recent relevant research, in which the pervasive role of evolution within ecological networks is manifest. Although evolutionary processes are most evident at macroevolutionary scales, they are also important drivers of local network structure and dynamics. We propose components of a blueprint for further research, emphasising process-based models, experimental evolution, and phenotypic variation, across a range of distinct spatial and temporal scales. Evolutionary dimensions are required to advance our understanding of foundational properties of community assembly and to enhance our capability of predicting how networks will respond to impending changes.

Diversity and Host Interactions Among Virulent and Temperate Baltic Sea Flavobacterium Phages

Viruses in aquatic environments play a key role in microbial population dynamics and nutrient cycling. In particular, bacteria of the phylum Bacteriodetes are known to participate in recycling algal blooms. Studies of phage–host interactions involving this phylum are hence important to understand the processes shaping bacterial and viral communities in the ocean as well as nutrient cycling. In this study, we isolated and sequenced three strains of flavobacteria—LMO6, LMO9, LMO8—and 38 virulent phages infecting them. These phages represent 15 species, occupying three novel genera. Additionally, one temperate phage was induced from LMO6 and was found to be competent at infecting LMO9. Functions could be predicted for a limited number of phage genes, mainly representing roles in DNA replication and virus particle formation. No metabolic genes were detected. While the phages isolated on LMO8 could infect all three bacterial strains, the LMO6 and LMO9 phages could not infect LMO8. Of the phages isolated on LMO9, several showed a host-derived reduced efficiency of plating on LMO6, potentially due to differences in DNA methyltransferase genes. Overall, these phage–host systems contribute novel genetic information to our sequence databases and present valuable tools for the study of both virulent and temperate phages.

Filamentous phages reduce bacterial growth in low salinities

Being non-lytic, filamentous phages can replicate at high frequencies and often carry virulence factors, which are important in the evolution and emergence of novel pathogens. However, their net effect on bacterial fitness remains unknown. To understand the ecology and evolution between filamentous phages and their hosts, it is important to assess (i) fitness effects of filamentous phages on their hosts and (ii) how these effects depend on the environment. To determine how the net effect on bacterial fitness by filamentous phages changes across environments, we constructed phage-bacteria infection networks at ambient 15 practical salinity units (PSU) and stressful salinities (11 and 7 PSU) using the marine bacterium, Vibrio alginolyticus and its derived filamentous phages as model system. We observed no significant difference in network structure at 15 and 11 PSU. However, at 7 PSU phages significantly reduced bacterial growth changing network structure. This pattern was mainly driven by a significant increase in bacterial susceptibility. Our findings suggest that filamentous phages decrease bacterial growth, an indirect measure of fitness in stressful environmental conditions, which might impact bacterial communities, alter horizontal gene transfer events and possibly favour the emergence of novel pathogens in environmental Vibrios.

Genomic Signatures for Sedimentary Microbial Utilization of Phytoplankton Detritus in a Fast-Flowing Estuary

In fast-flowing, river-dominated estuaries, “hotspots” of microbial biogeochemical cycling can be found within areas of extended water retention. Lateral bays located off of the North and South channels of the Columbia River estuary are proposed to be such hotspots. Previous metagenomic studies on water samples indicated that these regions function both as sources and sinks of biogenic particles, with potential to impact organic matter fluxes in the estuary. To extend this work, we analyzed 11 sediment metagenomes from three disparate bays: the freshwater Cathlamet Bay, and the brackish Youngs Bay and more saline Baker Bay located nearer the mouth to the south and north of the main channel, respectively. Samples were collected from upper layers of sediments in August of 2011 and 2013 for DNA extraction and metagenome sequencing. All metagenomes were dominated by bacterial sequences, although diatom sequences as high as 26% of the total annotated sequences were observed in the higher salinity samples. Unsupervised 2D hierarchical clustering analysis resulted in the eleven metagenome samples clustered into four groups by microbial taxonomic composition, with Bacteroides, diatom, and phage levels driving most of the grouping. Results of functional gene clustering further indicated that diatom bloom degradation stage (early vs. late) was an important factor. While the Flavobacteriia and Cytophagia classes were well represented in metagenomes containing abundant diatoms, taxa from the Bacteroidia class, along with certain members of the Sphingobacteriia class, were particularly abundant in metagenomes representing later stages of diatom decomposition. In contrast, the sediment metagenomes with low relative abundance of diatom and Bacteroidetes sequences appeared to have a metabolic potential biased toward microbial growth under nutrient limitation. While differences in water salinity clearly also influenced the microbial community composition and metabolic potential, our results highlight a central role for allochthonous labile organic matter (i.e., diatom detritus), in shaping bacterial taxonomic and functional properties in the Columbia River estuary lateral bay sediments. These results suggest that in fast-flowing, river-dominated estuaries, sediment microbial communities in areas of extended water retention, such as the lateral bays, may contribute disproportionately to estuarine organic matter degradation and recycling.

Complete genome analysis of a novel temperate bacteriophage induced from Corynebacterium striatum

A temperate bacteriophage, IME1320_01, was induced by mitomycin C treatment from Corynebacterium striatum. This phage possesses a double-stranded DNA genome of 40,086 bp with a G+C content of 58%. A total of 53 putative open reading frames (ORFs) were identified in its genome. BLASTn analysis revealed that IME1320_01 had the highest sequence similarity to Corynebacterium striatum strain 216, with a genome homology coverage of 44% and highest sequence identity of 95%. The termini of the phage genome was non-redundant, with a 13-nt 3'-protruding cohesive end. To the best of our knowledge, phage IME1320_01 is the first inducible phage to be identified in Corynebacterium striatum.

Revisiting the Concept of Host Range of Plant Pathogens

Strategies to manage plant disease-from use of resistant varieties to crop rotation, elimination of reservoirs, landscape planning, surveillance, quarantine, risk modeling, and anticipation of disease emergences-all rely on knowledge of pathogen host range. However, awareness of the multitude of factors that influence the outcome of plant-microorganism interactions, the spatial and temporal dynamics of these factors, and the diversity of any given pathogen makes it increasingly challenging to define simple, all-purpose rules to circumscribe the host range of a pathogen. For bacteria, fungi, oomycetes, and viruses, we illustrate that host range is often an overlapping continuum-more so than the separation of discrete pathotypes-and that host jumps are common. By setting the mechanisms of plant-pathogen interactions into the scales of contemporary land use and Earth history, we propose a framework to assess the frontiers of host range for practical applications and research on pathogen evolution.

Role of host genetic diversity for susceptibility-to-infection in the evolution of virulence of a plant virus

Predicting viral emergence is difficult due to the stochastic nature of the underlying processes and the many factors that govern pathogen evolution. Environmental factors affecting the host, the pathogen and the interaction between both are key in emergence. In particular, infectious disease dynamics are affected by spatiotemporal heterogeneity in their environments. A broad knowledge of these factors will allow better estimating where and when viral emergence is more likely to occur. Here, we investigate how the population structure for susceptibility-to-infection genes of the plant Arabidopsis thaliana shapes the evolution of Turnip mosaic virus (TuMV). For doing so we have evolved TuMV lineages in two radically different host population structures: (1) a metapopulation subdivided into six demes (subpopulations); each one being composed of individuals from only one of six possible A. thaliana ecotypes and (2) a well-mixed population constituted by equal number of plants from the same six A. thaliana ecotypes. These two populations were evolved for twelve serial passages. At the end of the experimental evolution, we found faster adaptation of TuMV to each ecotype in the metapopulation than in the well-mixed heterogeneous host populations. However, viruses evolved in well-mixed populations were more pathogenic and infectious than viruses evolved in the metapopulation. Furthermore, the viruses evolved in the demes showed stronger signatures of local specialization than viruses evolved in the well-mixed populations. These results illustrate how the genetic diversity of hosts in an experimental ecosystem favors the evolution of virulence of a pathogen.

Understanding and Exploiting Phage-Host Interactions

Initially described a century ago by William Twort and Felix d’Herelle, bacteriophages are bacterial viruses found ubiquitously in nature, located wherever their host cells are present. Translated literally, bacteriophage (phage) means ‘bacteria eater’. Phages interact and infect specific bacteria while not affecting other bacteria or cell lines of other organisms. Due to the specificity of these phage–host interactions, the relationship between phages and their host cells has been the topic of much research. The advances in phage biology research have led to the exploitation of these phage–host interactions and the application of phages in the agricultural and food industry. Phages may provide an alternative to the use of antibiotics, as it is well known that the emergence of antibiotic-resistant bacterial infections has become an epidemic in clinical settings. In agriculture, pre-harvest and/or post-harvest application of phages to crops may prevent the colonisation of bacteria that are detrimental to plant or human health. In addition, the abundance of data generated from genome sequencing has allowed the development of phage-derived bacterial detection systems of foodborne pathogens. This review aims to outline the specific interactions between phages and their host and how these interactions may be exploited and applied in the food industry.

Coevolutionary dynamics shape the structure of bacteria-phage infection networks

Coevolution-reciprocal evolutionary change among interacting species driven by natural selection-is thought to be an important force in shaping biodiversity. This ongoing process takes place within tangled networks of species interactions. In microbial communities, evolutionary change between hosts and parasites occurs at the same time scale as ecological change. Yet, we still lack experimental evidence of the role of coevolution in driving changes in the structure of such species interaction networks. Filling this gap is important because network structure influences community persistence through indirect effects. Here, we quantified experimentally to what extent coevolutionary dynamics lead to contrasting patterns in the architecture of bacteria-phage infection networks. Specifically, we look at the tendency of these networks to be organized in a nested pattern by which the more specialist phages tend to infect only a proper subset of those bacteria infected by the most generalist phages. We found that interactions between coevolving bacteria and phages become less nested over time under fluctuating dynamics, and more nested under arms race dynamics. Moreover, when coevolution results in high average infectivity, phages and bacteria differ more from each other over time under arms race dynamics than under fluctuating dynamics. The tradeoff between the fitness benefits of evolving resistance/infectivity traits and the costs of maintaining them might explain these differences in network structure. Our study shows that the interaction pattern between bacteria and phages at the community level depends on the way coevolution unfolds.

The genera of bacteriophages and their receptors are the major determinants of host range

The host range of phages is a key to understand their impact on bacterial ecology and evolution. Because of the complexity of phage-host interactions, the variables that determine the breadth of a phage host range remain poorly understood. Here, we propose a novel holistic approach to identify the host range determinants of a new collection of phages infecting Salmonella, isolated from animal, environmental and wastewater samples that were able to infect 58 of the 71 Salmonella strains in our collection. By using a set of statistic approaches (non-metric dimensional scaling, Bray-Curtis distance, PERMANOVA), we analysed phenotypic (host range on wild-type and receptor mutants) and genetic data (taxonomic assignment and receptor binding proteins) to evaluate the impact of isolation strain and niche, phage receptor and genus on the host range. Statistical analysis revealed that two phage characteristics influence the host range by explaining the most variance: the receptor by 45% and the genus by 51%. Interestingly, phage genus and receptor in combination explained 79% of the variance, establishing these characteristics as the major determinants of the host range. This study demonstrates the power and the novelty of applying statistical approaches to phenotypic and genetic data to investigate the ecology of phage-host interactions.

The first known virus isolates from Antarctic sea ice have complex infection patterns

Viruses are recognized as important actors in ocean ecology and biogeochemical cycles, but many details are not yet understood. We participated in a winter expedition to the Weddell Sea, Antarctica, to isolate viruses and to measure virus-like particle abundance (flow cytometry) in sea ice. We isolated 59 bacterial strains and the first four Antarctic sea-ice viruses known (PANV1, PANV2, OANV1 and OANV2), which grow in bacterial hosts belonging to the typical sea-ice genera Paraglaciecola and Octadecabacter. The viruses were specific for bacteria at the strain level, although OANV1 was able to infect strains from two different classes. Both PANV1 and PANV2 infected 11/15 isolated Paraglaciecola strains that had almost identical 16S rRNA gene sequences, but the plating efficiencies differed among the strains, whereas OANV1 infected 3/7 Octadecabacter and 1/15 Paraglaciecola strains and OANV2 1/7 Octadecabacter strains. All the phages were cold-active and able to infect their original host at 0°C and 4°C, but not at higher temperatures. The results showed that virus-host interactions can be very complex and that the viral community can also be dynamic in the winter-sea ice.

Analysing ecological networks of species interactions

Network approaches to ecological questions have been increasingly used, particularly in recent decades. The abstraction of ecological systems - such as communities - through networks of interactions between their components indeed provides a way to summarize this information with single objects. The methodological framework derived from graph theory also provides numerous approaches and measures to analyze these objects and can offer new perspectives on established ecological theories as well as tools to address new challenges. However, prior to using these methods to test ecological hypotheses, it is necessary that we understand, adapt, and use them in ways that both allow us to deliver their full potential and account for their limitations. Here, we attempt to increase the accessibility of network approaches by providing a review of the tools that have been developed so far, with - what we believe to be - their appropriate uses and potential limitations. This is not an exhaustive review of all methods and metrics, but rather, an overview of tools that are robust, informative, and ecologically sound. After providing a brief presentation of species interaction networks and how to build them in order to summarize ecological information of different types, we then classify methods and metrics by the types of ecological questions that they can be used to answer from global to local scales, including methods for hypothesis testing and future perspectives. Specifically, we show how the organization of species interactions in a community yields different network structures (e.g., more or less dense, modular or nested), how different measures can be used to describe and quantify these emerging structures, and how to compare communities based on these differences in structures. Within networks, we illustrate metrics that can be used to describe and compare the functional and dynamic roles of species based on their position in the network and the organization of their interactions as well as associated new methods to test the significance of these results. Lastly, we describe potential fruitful avenues for new methodological developments to address novel ecological questions.

Characterization of Bacteriophage vB-EcoS-95, Isolated From Urban Sewage and Revealing Extremely Rapid Lytic Development

Morphological, biological, and genetic characteristics of a virulent Siphoviridae phage, named vB-EcoS-95, is reported. This phage was isolated from urban sewage. It was found to infect some Escherichia coli strains giving clear plaques. The genome of this phage is composed of 50,910 bp and contains 89 ORFs. Importantly, none of the predicted ORFs shows any similarity with known pathogenic factors that would prevent its use in medicine. Genome sequence analysis of vB-EcoS-95 revealed 74% similarity to genomic sequence of Shigella phage pSf-1. Compared to pSf-1, phage vb-EcoS-95 does not infect Shigella strains and has an efficient bacteriolytic activity against some E. coli strains. One-step growth analysis revealed that this phage has a very short latent period (4 min), and average burst size of 115 plaque forming units per cell, which points to its high infectivity of host cells and strong lytic activity. The bacteriolytic effect of vB-EcoS-95 was tested also on biofilm-producing strains. These results indicate that vB-EcoS-95 is a newly discovered E. coli phage that may be potentially used to control the formation of biofilms.

Evolution of plant-virus interactions: host range and virus emergence

Changes in host range are central to virus emergence. Host range, together with its evolution, is determined by virus intrinsic factors, such as genetic traits determining its fitness in different hosts. Experimental analyses have shown the relevance in host range evolution of across-host fitness trade-offs. Host range is also determined by ecological factors extrinsic to the virus such as the distribution, abundance, and interaction of species, and understanding their role in host range evolution is a current challenge. Indeed, intrinsic and extrinsic factors, and the complexity of biotic and abiotic interactions, must be considered in order to provide generalisations on patterns of transmission, host range evolution, and disease emergence. This exciting new field of research is still in its infancy.

Modeling Multispecies Gene Flow Dynamics Reveals the Unique Roles of Different Horizontal Gene Transfer Mechanisms

Horizontal gene transfer within diverse bacterial populations occurs through multiple mechanisms of exchange. The most established routes of gene transfer, transduction, transformation, and conjugation, have been characterized in detail, revealing the advantages and limitations of each mechanism. More recently, interspecies gene exchange via extracellular vesicles has been reported and characterized, making vesicle-mediated exchange a fourth, general mechanism of gene transfer. Despite an understanding of each individual pathway, how all of these mechanisms act in concert has not been explored. Here we develop a model of gene exchange in a multispecies bacterial community that takes into account the rates and limitations of all four gene transfer mechanisms. Our results reveal unique roles for each gene exchange mechanism, and highlight how multiple pathways working together are required for widespread gene exchange within diverse bacterial populations.

Oil reservoirs, an exceptional habitat for microorganisms

•

Water-containing parts within oil reservoirs extend the zone of biodegradation.

•

Biodegradation is controlled by environmental factors.

•

Proteobacteria and Euryarchaeota are ubiquitous in oil reservoirs over all temperature ranges.

•

Biofilms as microbial adaption in oil reservoirs.

•

Viruses as potential control for microbial activity and function.

Microorganisms are present in oil reservoirs around the world where they degrade oil and lead to changes in oil quality. Unfortunately, our knowledge about processes in deep oil reservoirs is limited due to the lack of undisturbed samples. In this review, we discuss the distribution of microorganisms at the oil-water transition zone as well as in water saturated parts of the oil leg and their possible physiological adaptations to abiotic and biotic ecological factors such as temperature, salinity and viruses. We show the importance of studying the water phase within the oil, because small water inclusions and pockets within the oil leg provide an exceptional habitat for microorganisms within a natural oil reservoir and concurrently enlarge the zone of oil biodegradation. Environmental factors such as temperature and salinity control oil biodegradation. Temperature determines the type of microorganisms which are able to inhabit the reservoir. Proteobacteria and Euryarchaeota, are ubiquitous in oil reservoirs over all temperature ranges, whereas some others are tied to specific temperatures. It is proposed that biofilm formation is the dominant way of life within oil reservoirs, enhancing nutrient uptake, syntrophic interactions and protection against environmental stress. Literature shows that viruses are abundant in oil reservoirs and the possible impact on microbial community composition due to control of microbial activity and function is discussed.

Computational prediction of inter-species relationships through omics data analysis and machine learning

BACKGROUND

Antibiotic resistance and its rapid dissemination around the world threaten the efficacy of currently-used medical treatments and call for novel, innovative approaches to manage multi-drug resistant infections. Phage therapy, i.e., the use of viruses (phages) to specifically infect and kill bacteria during their life cycle, is one of the most promising alternatives to antibiotics. It is based on the correct matching between a target pathogenic bacteria and the therapeutic phage. Nevertheless, correctly matching them is a major challenge. Currently, there is no systematic method to efficiently predict whether phage-bacterium interactions exist and these pairs must be empirically tested in laboratory. Herein, we present our approach for developing a computational model able to predict whether a given phage-bacterium pair can interact based on their genome.

RESULTS

Based on public data from GenBank and phagesDB.org, we collected more than a thousand positive phage-bacterium interactions with their complete genomes. In addition, we generated putative negative (i.e., non-interacting) pairs. We extracted, from the collected genomes, a set of informative features based on the distribution of predictive protein-protein interactions and on their primary structure (e.g. amino-acid frequency, molecular weight and chemical composition of each protein). With these features, we generated multiple candidate datasets to train our algorithms. On this base, we built predictive models exhibiting predictive performance of around 90% in terms of F1-score, sensitivity, specificity, and accuracy, obtained on the test set with 10-fold cross-validation.

CONCLUSION

These promising results reinforce the hypothesis that machine learning techniques may produce highly-predictive models accelerating the search of interacting phage-bacteria pairs.

Cross-resistance is modular in bacteria-phage interactions

Phages shape the structure of natural bacterial communities and can be effective therapeutic agents. Bacterial resistance to phage infection, however, limits the usefulness of phage therapies and could destabilise community structures, especially if individual resistance mutations provide cross-resistance against multiple phages. We currently understand very little about the evolution of cross-resistance in bacteria–phage interactions. Here we show that the network structure of cross-resistance among spontaneous resistance mutants of Pseudomonas aeruginosa evolved against each of 27 phages is highly modular. The cross-resistance network contained both symmetric (reciprocal) and asymmetric (nonreciprocal) cross-resistance, forming two cross-resistance modules defined by high within- but low between-module cross-resistance. Mutations conferring cross-resistance within modules targeted either lipopolysaccharide or type IV pilus biosynthesis, suggesting that the modularity of cross-resistance was structured by distinct phage receptors. In contrast, between-module cross-resistance was provided by mutations affecting the alternative sigma factor, RpoN, which controls many lifestyle-associated functions, including motility, biofilm formation, and quorum sensing. Broader cross-resistance range was not associated with higher fitness costs or weaker resistance against the focal phage used to select resistance. However, mutations in rpoN, providing between-module cross-resistance, were associated with higher fitness costs than mutations associated with within-module cross-resistance, i.e., in genes encoding either lipopolysaccharide or type IV pilus biosynthesis. The observed structure of cross-resistance predicted both the frequency of resistance mutations and the ability of phage combinations to suppress bacterial growth. These findings suggest that the evolution of cross-resistance is common, is likely to play an important role in the dynamic structure of bacteria–phage communities, and could inform the design principles for phage therapy treatments.

Phage therapy is a promising alternative to antibiotics for treating bacterial infections. Yet as with antibiotics, bacteria readily evolve resistance to phage attack, including cross-resistance that protects against multiple phages at once and so limits the usefulness of phage cocktails. Here we show, using laboratory experimental evolution of resistance against 27 phages in P. aeruginosa, that cross-resistance is common and determines the ability of phage combinations to suppress bacterial growth. Using whole-genome sequencing, we show that cross-resistance is most common against multiple phages that use the same receptor but that global regulator mutations provide generalist resistance, probably by simultaneously affecting the expression of multiple different phage receptors. Future trials should test if these features of cross-resistance evolution translate to more complex in vivo environments and can therefore be exploited to design more effective phage therapies for the clinic.