The Environment and Cyanophage Diversity: Insights From Environmental Sequencing of DNA Polymerase

How might bacteriophages shape biological invasions?

Invasions by eukaryotes dependent on environmentally acquired bacterial mutualists are often limited by the ability of bacterial partners to survive and establish free-living populations. Focusing on the model legume-rhizobium mutualism, we apply invasion biology hypotheses to explain how bacteriophages can impact the competitiveness of introduced bacterial mutualists. Predicting how phage-bacteria interactions affect invading eukaryotic hosts requires knowing the eco-evolutionary constraints of introduced and native microbial communities, as well as their differences in abundance and diversity. By synthesizing research from invasion biology, as well as bacterial, viral, and community ecology, we create a conceptual framework for understanding and predicting how phages can affect biological invasions through their effects on bacterial mutualists.

Diverse Marine T4-like Cyanophage Communities Are Primarily Comprised of Low-Abundance Species Including Species with Distinct Seasonal, Persistent, Occasional, or Sporadic Dynamics

Cyanophages exert important top-down controls on their cyanobacteria hosts; however, concurrent analysis of both phage and host populations is needed to better assess phage-host interaction models. We analyzed picocyanobacteria Prochlorococcus and Synechococcus and T4-like cyanophage communities in Pacific Ocean surface waters using five years of monthly viral and cellular fraction metagenomes. Cyanophage communities contained thousands of mostly low-abundance (<2% relative abundance) species with varying temporal dynamics, categorized as seasonally recurring or non-seasonal and occurring persistently, occasionally, or sporadically (detected in ≥85%, 15-85%, or <15% of samples, respectively). Viromes contained mostly seasonal and persistent phages (~40% each), while cellular fraction metagenomes had mostly sporadic species (~50%), reflecting that these sample sets capture different steps of the infection cycle-virions from prior infections or within currently infected cells, respectively. Two groups of seasonal phages correlated to Synechococcus or Prochlorococcus were abundant in spring/summer or fall/winter, respectively. Cyanophages likely have a strong influence on the host community structure, as their communities explained up to 32% of host community variation. These results support how both seasonally recurrent and apparent stochastic processes, likely determined by host availability and different host-range strategies among phages, are critical to phage-host interactions and dynamics, consistent with both the Kill-the-Winner and the Bank models.

A Review of Cyanophage–Host Relationships: Highlighting Cyanophages as a Potential Cyanobacteria Control Strategy

Harmful algal blooms (HABs) are naturally occurring phenomena, and cyanobacteria are the most commonly occurring HABs in freshwater systems. Cyanobacteria HABs (cyanoHABs) negatively affect ecosystems and drinking water resources through the production of potent toxins. Furthermore, the frequency, duration, and distribution of cyanoHABs are increasing, and conditions that favor cyanobacteria growth are predicted to increase in the coming years. Current methods for mitigating cyanoHABs are generally short-lived and resource-intensive, and have negative impacts on non-target species. Cyanophages (viruses that specifically target cyanobacteria) have the potential to provide a highly specific control strategy with minimal impacts on non-target species and propagation in the environment. A detailed review (primarily up to 2020) of cyanophage lifecycle, diversity, and factors influencing infectivity is provided in this paper, along with a discussion of cyanophage and host cyanobacteria relationships for seven prominent cyanoHAB-forming genera in North America, including: Synechococcus, Microcystis, Dolichospermum, Aphanizomenon, Cylindrospermopsis, Planktothrix, and Lyngbya. Lastly, factors affecting the potential application of cyanophages as a cyanoHAB control strategy are discussed, including efficacy considerations, optimization, and scalability for large-scale applications.

Cyanophage Distribution Across European Lakes of the Temperate-Humid Continental Climate Zone Assessed Using PCR-Based Genetic Markers

Studies of the diversity and distribution of freshwater cyanophages are generally limited to the small geographical areas, in many cases including only one or few lakes. Data from dozens of various lakes distributed at a larger distance are necessary to understand their spatial distribution and sensitivity to biotic and abiotic factors. Thus, the objective of this study was to analyze the diversity and distribution of cyanophages within the infected cells using marker genes (psbA, nblA, and g91) in 21 Polish and Lithuanian lakes. Physicochemical factors that might be related to them were also analyzed. The results demonstrated that genetic markers representing cyanophages were observed in most lakes studied. The frequently detected gene was psbA with 88% of cyanophage-positive samples, while nblA and g91 were found in approximately 50% of lakes. The DNA sequence analyses for each gene demonstrated low variability between them, although the psbA sequences branched within the larger cluster of marine Synechoccocuss counterparts. The principal component analysis allowed to identify significant variation between the lakes that presented high and low cyanobacterial biomass. The lakes with high cyanobacterial biomass were further separated by country and the different diversity of cyanobacteria species, particularly Planktothrix agardhii, was dominant in the Polish lakes and Planktolyngbya limnetica in the Lithuanian lakes. The total phosphorous and the presence of cyanophage genes psbA and nblA were the most important factors that allowed differentiation for the Polish lakes, while the pH and the genes g91 and nblA for the Lithuanian lakes.

Metagenomic community composition and resistome analysis in a full-scale cold climate wastewater treatment plant

BACKGROUND

Wastewater treatment plants are an essential part of maintaining the health and safety of the general public. However, they are also an anthropogenic source of antibiotic resistance genes. In this study, we characterized the resistome, the distribution of classes 1-3 integron-integrase genes (intI1, intI2, and intI3) as mobile genetic element biomarkers, and the bacterial and phage community compositions in the North End Sewage Treatment Plant in Winnipeg, Manitoba. Samples were collected from raw sewage, returned activated sludge, final effluent, and dewatered sludge. A total of 28 bacterial and viral metagenomes were sequenced over two seasons, fall and winter. Integron-integrase genes, the 16S rRNA gene, and the coliform beta-glucuronidase gene were also quantified during this time period.

RESULTS

Bacterial classes observed above 1% relative abundance in all treatments were Actinobacteria (39.24% ± 0.25%), Beta-proteobacteria (23.99% ± 0.16%), Gamma-proteobacteria (11.06% ± 0.09%), and Alpha-proteobacteria (9.18 ± 0.04%). Families within the Caudovirales order: Siphoviridae (48.69% ± 0.10%), Podoviridae (23.99% ± 0.07%), and Myoviridae (19.94% ± 0.09%) were the dominant phage observed throughout the NESTP. The most abundant bacterial genera (in terms of average percent relative abundance) in influent, returned activated sludge, final effluent, and sludge, respectively, includes Mycobacterium (37.4%, 18.3%, 46.1%, and 7.7%), Acidovorax (8.9%, 10.8%, 5.4%, and 1.3%), and Polaromonas (2.5%, 3.3%, 1.4%, and 0.4%). The most abundant class of antibiotic resistance in bacterial samples was tetracycline resistance (17.86% ± 0.03%) followed by peptide antibiotics (14.24% ± 0.03%), and macrolides (10.63% ± 0.02%). Similarly, the phage samples contained a higher prevalence of macrolide (30.12% ± 0.30%), peptide antibiotic (10.78% ± 0.13%), and tetracycline (8.69% ± 0.11%) resistance. In addition, intI1 was the most abundant integron-integrase gene throughout treatment (1.14 × 104 gene copies/mL) followed by intI3 (4.97 × 103 gene copies/mL) while intI2 abundance remained low (6.4 × 101 gene copies/mL).

CONCLUSIONS

Wastewater treatment successfully reduced the abundance of bacteria, DNA phage and antibiotic resistance genes although many antibiotic resistance genes remained in effluent and biosolids. The presence of integron-integrase genes throughout treatment and in effluent suggests that antibiotic resistance genes could be actively disseminating resistance between both environmental and pathogenic bacteria.

Novel Freshwater Cyanophages Provide New Insights into Evolutionary Relationships between Freshwater and Marine Cyanophages

Cyanobacteria and cyanophages are present widely in both freshwater and marine environments. However, freshwater cyanophages remain unknown largely due to the small numbers of cyanophage isolates despite their ecological and environmental significance. In this study, we present the characterization of two novel lytic freshwater cyanophages isolated from a tropical inland lake in Singapore, namely, cyanopodovirus S-SRP01 and cyanomyovirus S-SRM01, infecting two different strains of Synechococcus spp. Functional annotation of S-SRP01 and S-SRM01 genomes revealed a high degree of homology with marine cyanophages. Phylogenetic trees of concatenated genes and whole-genome alignment provided further evidence that S-SRP01 is close evolutionarily to marine cyanopodoviruses, while S-SRM01 is evolutionarily close to marine cyanomyoviruses. Few genetic similarities between freshwater and marine cyanophages have been identified in previous studies. The isolation of S-SRP01 and S-SRM01 expand current knowledge on freshwater cyanophages infecting Synechococcus spp. Their high degree of gene sharing provides new insights into the evolutionary relationships between freshwater and marine cyanophages. This relatedness is further supported by the discovery of similar phenomenon from other freshwater viral metagenomes. IMPORTANCE This study expands the current knowledge on freshwater cyanophage isolates and cyanophage genetic diversity, indicating that freshwater and marine cyanophages infecting Synechococcus spp. may share close genetic similarity and evolutionary relationships.

Spatial patterns in phage-Rhizobium coevolutionary interactions across regions of common bean domestication

Bacteriophages play significant roles in the composition, diversity, and evolution of bacterial communities. Despite their importance, it remains unclear how phage diversity and phage-host interactions are spatially structured. Local adaptation may play a key role. Nitrogen-fixing symbiotic bacteria, known as rhizobia, have been shown to locally adapt to domesticated common bean at its Mesoamerican and Andean sites of origin. This may affect phage-rhizobium interactions. However, knowledge about the diversity and coevolution of phages with their respective Rhizobium populations is lacking. Here, through the study of four phage-Rhizobium communities in Mexico and Argentina, we show that both phage and host diversity is spatially structured. Cross-infection experiments demonstrated that phage infection rates were higher overall in sympatric rhizobia than in allopatric rhizobia except for one Argentinean community, indicating phage local adaptation and host maladaptation. Phage-host interactions were shaped by the genetic identity and geographic origin of both the phage and the host. The phages ranged from specialists to generalists, revealing a nested network of interactions. Our results suggest a key role of local adaptation to resident host bacterial communities in shaping the phage genetic and phenotypic composition, following a similar spatial pattern of diversity and coevolution to that in the host.

Role of Phylogenetic Structure in the Dynamics of Coastal Viral Assemblages

Marine microbes, including viruses, are an essential part of the marine ecosystem, forming the base of the food web and driving biogeochemical cycles. Within this system, the composition of viral assemblages changes markedly with time, and some of these changes are repeatable through time; however, the extent to which these dynamics are reflected within versus among evolutionarily related groups of viruses is largely unexplored. To examine these dynamics, changes in the composition of two groups of ecologically important viruses and communities of their potential hosts were sampled every 2 weeks for 13 months at a coastal site in British Columbia, Canada. We sequenced two marker genes for viruses-the gene encoding the major capsid protein of T4-like phages and their relatives (gp23) and the RNA-dependent RNA polymerase (RdRp) gene of marnavirus-like RNA viruses-as well as marker genes for their bacterial and eukaryotic host communities, the genes encoding 16S rRNA and 18S rRNA. There were strong lagged correlations between viral diversity and community similarity of putative hosts, implying that the viruses influenced the composition of the host communities. The results showed that for both viral assemblages, the dominant clusters of phylogenetically related viruses shifted over time, and this was correlated with environmental changes. Viral clusters contained many ephemeral taxa and few persistent taxa, but within a viral assemblage, the ephemeral and persistent taxa were closely related, implying ecological dynamics within these clusters. Furthermore, these dynamics occurred in both the RNA and DNA viral assemblages surveyed, implying that this structure is common in natural viral assemblages.IMPORTANCE Viruses are major agents of microbial mortality in marine systems, yet little is known about changes in the composition of viral assemblages in relation to those of the microbial communities that they infect. Here, we sampled coastal seawater every 2 weeks for 1 year and used high-throughput sequencing of marker genes to follow changes in the composition of two groups of ecologically important viruses, as well as the communities of bacteria and protists that serve as their respective hosts. Different subsets of genetically related viruses dominated at different times. These results demonstrate that although the genetic composition of viral assemblages is highly dynamic temporally, for the most part the shuffling of genotypes occurs within a few clusters of phylogenetically related viruses. Thus, it appears that even in temperate coastal waters with large seasonal changes, the highly dynamic shuffling of viral genotypes occurs largely within a few subsets of related individuals.

Genomic analysis of Synechococcus phage S-B43 and its adaption to the coastal environment

Synechococcus dominate picocyanobacterial communities in coastal environments. However, only a few Synechococcus phages have been described from the coastal seas of the Northwest Pacific Ocean. Here a new Synechococcus phage, S-B43 was isolated from the Bohai Sea, a semi-closed coastal sea of the Northwest Pacific Ocean. S-B43 is a member of Myoviridae, containing 275 predicted open reading frames. Fourteen auxiliary metabolic genes (AMG) were identified from the genome of S-B43, including five photosynthetic associated genes and several AMGs related to its adaption to the high turbidity and eutrophic coastal environment with a low ratio of phosphorus to nitrogen (HNLP). The occurrences of 31 AMGs among 34 cyanophage genomes indicates that AMGs zwf, gnd, speD, petF and those coding for FECH and thioredoxin were more common in coastal areas than in the open ocean and AMGs pebS and ho1 were more prevalent in the open ocean. The occurrence of cyanophage AMGs in different environments might be a reflection of the environmental adaption of their hosts. This study contributes to our understanding of the interactions between cyanobacteria and cyanophages and their environmental adaption to the coastal environment.

Iroki: automatic customization and visualization of phylogenetic trees

Phylogenetic trees are an important analytical tool for evaluating community diversity and evolutionary history. In the case of microorganisms, the decreasing cost of sequencing has enabled researchers to generate ever-larger sequence datasets, which in turn have begun to fill gaps in the evolutionary history of microbial groups. However, phylogenetic analyses of these types of datasets create complex trees that can be challenging to interpret. Scientific inferences made by visual inspection of phylogenetic trees can be simplified and enhanced by customizing various parts of the tree. Yet, manual customization is time-consuming and error prone, and programs designed to assist in batch tree customization often require programming experience or complicated file formats for annotation. Iroki, a user-friendly web interface for tree visualization, addresses these issues by providing automatic customization of large trees based on metadata contained in tab-separated text files. Iroki’s utility for exploring biological and ecological trends in sequencing data was demonstrated through a variety of microbial ecology applications in which trees with hundreds to thousands of leaf nodes were customized according to extensive collections of metadata. The Iroki web application and documentation are available at https://www.iroki.net or through the VIROME portal http://virome.dbi.udel.edu. Iroki’s source code is released under the MIT license and is available at https://github.com/mooreryan/iroki.