Modeling ecological drivers in marine viral communities using comparative metagenomics and network analyses

Diversity and potential host-interactions of viruses inhabiting deep-sea seamount sediments

Seamounts are globally distributed across the oceans and form one of the major oceanic biomes. Here, we utilized combined analyses of bulk metagenome and virome to study viral communities in seamount sediments in the western Pacific Ocean. Phylogenetic analyses and the protein-sharing network demonstrate extensive diversity and previously unknown viral clades. Inference of virus-host linkages uncovers extensive interactions between viruses and dominant prokaryote lineages, and suggests that viruses play significant roles in carbon, sulfur, and nitrogen cycling by compensating or augmenting host metabolisms. Moreover, temperate viruses are predicted to be prevalent in seamount sediments, which tend to carry auxiliary metabolic genes for host survivability. Intriguingly, the geographical features of seamounts likely compromise the connectivity of viral communities and thus contribute to the high divergence of viral genetic spaces and populations across seamounts. Altogether, these findings provides knowledge essential for understanding the biogeography and ecological roles of viruses in globally widespread seamounts.

Giant Virus Global Proteomics Innovation: Comparative Evaluation of In-Gel and In-Solution Digestion Methods

With their unusually large genome and particle sizes, giant viruses (GVs) defy the conventional definition of viruses. Although most GVs isolated infect unicellular protozoans, such as amoeba, studies in the last decade have established their much wider prevalence infecting most eukaryotic supergroups and some giant viral families with the potential to be human pathogens. Their complexity, almost autonomous life cycle, and enigmatic evolution necessitate the study of GVs. The accurate assessment of GV proteome is a veritable challenge. We have compared the coverage of global protein identification using different methods for GVs isolated in Mumbai, Mimivirus Bombay (MVB), Powai Lake Megavirus (PLMV), and Kurlavirus (KV), along with two previously studied GVs, Acanthamoeba polyphaga Mimivirus (APMV) and Marseillevirus (MV). Our study shows that the simultaneous use of in-gel and in-solution digestion methods can significantly increase the coverage of protein identification in the global proteome analysis of purified GV particles. Combining the two methods of analyses, we identified an additional 72 proteins in APMV and 114 in MV compared with what have been previously reported. Similarly, proteomes of MVB, PLMV, and KV were analyzed, and a total of 242 proteins in MVB, 287 proteins in PLMV, and 174 proteins in KV were identified. Our results suggest that a combined methodology of in-gel and in-solution methods is more efficient and opens up new avenues for innovation in global proteome analysis of GVs. Future planetary health research on GVs can benefit from consideration of a broader range of proteomics methodologies as illustrated by the present study.

Host and geography impact virus diversity in New Zealand's longfin and shortfin eels

The fishing and aquaculture industry is vital for global food security, yet viral diseases can result in mass fish die-off events. Determining the viromes of traditionally understudied species, such as fish, enhances our understanding of the global virosphere and the factors that influence virome composition and disease emergence. Very little is known about the viruses present in New Zealand's native fish species, including the shortfin eel (Anguilla australis) and the longfin eel (Anguilla dieffenbachii), both of which are fished culturally by Māori (the indigenous population of New Zealand) and commercially. Through a total RNA metatranscriptomic analysis of longfin and shortfin eels across three different geographic locations in the South Island of New Zealand, we aimed to determine whether viruses had jumped between the two eel species and whether eel virome composition was impacted by life stage, species, and geographic location. We identified nine viral species spanning eight different families, thereby enhancing our understanding of eel virus diversity in New Zealand and the host range of these viral families. Viruses of the family Flaviviridae (genus Hepacivirus) were widespread and found in both longfin and shortfin eels, indicative of cross-species transmission or virus-host co-divergence. Notably, both host specificity and geographic location appeared to influence eel virome composition, highlighting the complex interaction between viruses, hosts, and their ecosystems. This study broadens our understanding of viromes in aquatic hosts and highlights the importance of gaining baseline knowledge of fish viral abundance and diversity, particularly in aquatic species that are facing population declines.

Gut microbiota-driven metabolic alterations reveal gut-brain communication in Alzheimer's disease model mice

The gut microbiota (GM) and its metabolites affect the host nervous system and are involved in the pathogeneses of various neurological diseases. However, the specific GM alterations under pathogenetic pressure and their contributions to the "microbiota - metabolite - brain axis" in Alzheimer's disease (AD) remain unclear. Here, we investigated the GM and the fecal, serum, cortical metabolomes in APP/PS1 and wild-type (WT) mice, revealing distinct hub bacteria in AD mice within scale-free GM networks shared by both groups. Moreover, we identified diverse peripheral - central metabolic landscapes between AD and WT mice that featured bile acids (e.g. deoxycholic and isodeoxycholic acid) and unsaturated fatty acids (e.g. 11Z-eicosenoic and palmitoleic acid). Machine-learning models revealed the relationships between the differential/hub bacteria and these metabolic signatures from the periphery to the brain. Notably, AD-enriched Dubosiella affected AD occurrence via cortical palmitoleic acid and vice versa. Considering the transgenic background of the AD mice, we propose that Dubosiella enrichment impedes AD progression via the synthesis of palmitoleic acid, which has protective properties against inflammation and metabolic disorders. We identified another association involving fecal deoxycholic acid-mediated interactions between the AD hub bacteria Erysipelatoclostridium and AD occurrence, which was corroborated by the correlation between deoxycholate levels and cognitive scores in humans. Overall, this study elucidated the GM network alterations, contributions of the GM to peripheral - central metabolic landscapes, and mediatory roles of metabolites between the GM and AD occurrence, thus revealing the critical roles of bacteria in AD pathogenesis and gut - brain communications under pathogenetic pressure.

The virome in allergy and asthma: A nascent, ineffable player

Allergic diseases can be affected by virus-host interactions and are increasingly linked with the tissue-specific microbiome. High-throughput metagenomic sequencing has offered the opportunity to study the presence of viruses as an ecologic system, namely, the virome. Even though virome studies are technically challenging conceptually and analytically, they are already producing novel data expanding our understanding of the pathophysiologic mechanisms related to chronic inflammation and allergy. The importance of interspecies and intraspecies interactions is becoming apparent, as they can significantly, directly or indirectly, affect the host's response and antigenic state. Here, we emphasize the challenges and potential insights related to study of the virome in the context of allergy and asthma. We review the limited number of studies that have investigated the virome in these conditions, underlining the need for prospective, repeated sampling designs to unravel the virome's impact on disease development and its interplay with microbiota and immunity. The potential therapeutic use of bacteriophages, which are highly complex components of the virome, is discussed. There is clearly a need for further in-depth investigation of the virome as a system in allergic diseases.

An extended catalog of integrated prophages in the infant and adult fecal microbiome shows high prevalence of lysogeny

Background and aims

The acquisition and gradual maturation of gut microbial communities during early childhood is central to an individual's healthy development. Bacteriophages have the potential to shape the gut bacterial communities. However, the complex ecological interactions between phages and their bacterial host are still poorly characterized. In this study, we investigated the abundance and diversity of integrated prophages in infant and adult gut bacteria by detecting integrated prophages in metagenome assembled genomes (MAGs) of commensal bacteria.

Methods

Our study included 88 infants sampled at 3 weeks, 3 months, 6 months, and 12 months (n = 323 total samples), and their parents around delivery time (n = 138 total samples). Fecal DNA was extracted and characterized by using shotgun metagenomic sequencing, and a collection of prokaryotic MAGs was generated. The MAG collection was screened for the presence of integrated bacteriophage sequences, allowing their taxonomic and functional characterization.

Results

A large collection of 6,186 MAGs from infant and adult gut microbiota was obtained and screened for integrated prophages, allowing the identification of 7,165 prophage sequences longer than 10 kb. Strikingly, more than 70% of the near-complete MAGs were identified as lysogens. The prevalence of prophages in MAGs varied across bacterial families, with a lower prevalence observed among Coriobacteriaceae, Eggerthellaceae, Veillonellaceae and Burkholderiaceae, while a very high prevalence of lysogen MAGs were observed in Oscillospiraceae, Enterococcaceae, and Enterobacteriaceae. Interestingly for several bacterial families such as Bifidobacteriaceae and Bacteroidaceae, the prevalence of prophages in MAGs was higher in early infant time point (3 weeks and 3 months) than in later sampling points (6 and 12 months) and in adults. The prophage sequences were clustered into 5,616 species-like vOTUs, 77% of which were novel. Finally, we explored the functional repertoire of the potential auxiliary metabolic genes carried by these prophages, encoding functions involved in carbohydrate metabolism and degradation, amino acid metabolism and carbon metabolism.

Conclusion

Our study provides an enhanced understanding of the diversity and prevalence of lysogens in infant and adult gut microbiota and suggests a complex interplay between prophages and their bacterial hosts.

Highly diverse and unknown viruses may enhance Antarctic endoliths' adaptability

BACKGROUND

Rock-dwelling microorganisms are key players in ecosystem functioning of Antarctic ice free-areas. Yet, little is known about their diversity and ecology, and further still, viruses in these communities have been largely unexplored despite important roles related to host metabolism and nutrient cycling. To begin to address this, we present a large-scale viral catalog from Antarctic rock microbial communities.

RESULTS

We performed metagenomic analyses on rocks from across Antarctica representing a broad range of environmental and spatial conditions, and which resulted in a predicted viral catalog comprising > 75,000 viral operational taxonomic units (vOTUS). We found largely undescribed, highly diverse and spatially structured virus communities which had predicted auxiliary metabolic genes (AMGs) with functions indicating that they may be potentially influencing bacterial adaptation and biogeochemistry.

CONCLUSION

This catalog lays the foundation for expanding knowledge of virosphere diversity, function, spatial ecology, and dynamics in extreme environments. This work serves as a step towards exploring adaptability of microbial communities in the face of a changing climate. Video Abstract.

Dynamic network modeling of gut microbiota during Alzheimer's disease progression in mice

The intimate association between the gut microbiota (GM) and the central nervous system points to potential intervention strategies for neurological diseases. Nevertheless, there is currently no theoretical framework for selecting the window period and target bacteria for GM interventions owing to the complexity of the gut microecosystem. In this study, we constructed a complex network-based modeling approach to evaluate the topological features of the GM and infer the window period and bacterial candidates for GM interventions. We used Alzheimer's disease (AD) as an example and traced the GM dynamic changes in AD and wild-type mice at one, two, three, six, and nine months of age. The results revealed alterations of the topological features of the GM from a scale-free network into a random network during AD progression, indicating severe GM disequilibrium at the late stage of AD. Through stability and vulnerability assessments of the GM networks, we identified the third month after birth as the optimal window period for GM interventions in AD mice. Further computational simulations and robustness evaluations determined that the hub bacteria were potential candidates for GM interventions. Moreover, our GM functional analysis suggested that Lachnospiraceae UCG-001 - the hub and enriched bacterium in AD mice - was the keystone bacterium for GM interventions owing to its contributions to quinolinic acid synthesis. In conclusion, this study established a complex network-based modeling approach as a practical strategy for disease interventions from the perspective of the gut microecosystem.

Assembly and Annotation of Viral Metagenomes from Short-Read Sequencing Data

Viral metagenomics enables the detection, characterization, and quantification of viral sequences present in shotgun-sequenced datasets of purified virus-like particles and whole metagenomes. Next generation sequencing (Illumina) derived short single or paired-end read runs are a principal platform for metagenomics, and assembly of short reads allows for the identification of distinguishing viral signatures and complex genomic features for taxonomy and functional annotation. Here we describe the identification and characterization of viral genome sequences, bacteriophages, and eukaryotic viruses, from a cohort of human stool samples, using multiple methods. Following the purification of virus-like particles, sequencing, quality refinement, and genome assembly, we begin the protocol with raw short reads deposited in an open-source nucleotide archive. We highlight the use of VIBRANT, an automated computational tool for the characterization of microbial viruses and their viral community function. Finally, we also describe an alternative assembly-free option of mapping reads to established databases of reference genomes and previously characterized metagenome-assembled viral genomes.

Patterns and ecological drivers of viral communities in acid mine drainage sediments across Southern China

Recent advances in environmental genomics have provided unprecedented opportunities for the investigation of viruses in natural settings. Yet, our knowledge of viral biogeographic patterns and the corresponding drivers is still limited. Here, we perform metagenomic deep sequencing on 90 acid mine drainage (AMD) sediments sampled across Southern China and examine the biogeography of viruses in this extreme environment. The results demonstrate that prokaryotic communities dictate viral taxonomic and functional diversity, abundance and structure, whereas other factors especially latitude and mean annual temperature also impact viral populations and functions. In silico predictions highlight lineage-specific virus-host abundance ratios and richness-dependent virus-host interaction structure. Further functional analyses reveal important roles of environmental conditions and horizontal gene transfers in shaping viral auxiliary metabolic genes potentially involved in phosphorus assimilation. Our findings underscore the importance of both abiotic and biotic factors in predicting the taxonomic and functional biogeographic dynamics of viruses in the AMD sediments.

The biogeography of viral communities in extreme environments remains understudied. Here, the authors use metagenomic sequencing on 90 acid mine drainage sediments sampled across Southern China, showing the predominant effects of prokaryotic communities and the influence of environmental variables on viral taxonomy and function.

Viral community analysis in a marine oxygen minimum zone indicates increased potential for viral manipulation of microbial physiological state

Microbial communities in oxygen minimum zones (OMZs) are known to have significant impacts on global biogeochemical cycles, but viral influence on microbial processes in these regions are much less studied. Here we provide baseline ecological patterns using microscopy and viral metagenomics from the Eastern Tropical North Pacific (ETNP) OMZ region that enhance our understanding of viruses in these climate-critical systems. While extracellular viral abundance decreased below the oxycline, viral diversity and lytic infection frequency remained high within the OMZ, demonstrating that viral influences on microbial communities were still substantial without the detectable presence of oxygen. Viral community composition was strongly related to oxygen concentration, with viral populations in low-oxygen portions of the water column being distinct from their surface layer counterparts. However, this divergence was not accompanied by the expected differences in viral-encoded auxiliary metabolic genes (AMGs) relating to nitrogen and sulfur metabolisms that are known to be performed by microbial communities in these low-oxygen and anoxic regions. Instead, several abundant AMGs were identified in the oxycline and OMZ that may modulate host responses to low-oxygen stress. We hypothesize that this is due to selection for viral-encoded genes that influence host survivability rather than modulating host metabolic reactions within the ETNP OMZ. Together, this study shows that viruses are not only diverse throughout the water column in the ETNP, including the OMZ, but their infection of microorganisms has the potential to alter host physiological state within these biogeochemically important regions of the ocean.

iVirus 2.0: Cyberinfrastructure-supported tools and data to power DNA virus ecology

Microbes drive myriad ecosystem processes, but under strong influence from viruses. Because studying viruses in complex systems requires different tools than those for microbes, they remain underexplored. To combat this, we previously aggregated double-stranded DNA (dsDNA) virus analysis capabilities and resources into 'iVirus' on the CyVerse collaborative cyberinfrastructure. Here we substantially expand iVirus's functionality and accessibility, to iVirus 2.0, as follows. First, core iVirus apps were integrated into the Department of Energy's Systems Biology KnowledgeBase (KBase) to provide an additional analytical platform. Second, at CyVerse, 20 software tools (apps) were upgraded or added as new tools and capabilities. Third, nearly 20-fold more sequence reads were aggregated to capture new data and environments. Finally, documentation, as "live" protocols, was updated to maximize user interaction with and contribution to infrastructure development. Together, iVirus 2.0 serves as a uniquely central and accessible analytical platform for studying how viruses, particularly dsDNA viruses, impact diverse microbial ecosystems.

Cyanolichen microbiome contains novel viruses that encode genes to promote microbial metabolism

Lichen thalli are formed through the symbiotic association of a filamentous fungus and photosynthetic green alga and/or cyanobacterium. Recent studies have revealed lichens also host highly diverse communities of secondary fungal and bacterial symbionts, yet few studies have examined the viral component within these complex symbioses. Here, we describe viral biodiversity and functions in cyanolichens collected from across North America and Europe. As current machine-learning viral-detection tools are not trained on complex eukaryotic metagenomes, we first developed efficient methods to remove eukaryotic reads prior to viral detection and a custom pipeline to validate viral contigs predicted with three machine-learning methods. Our resulting high-quality viral data illustrate that every cyanolichen thallus contains diverse viruses that are distinct from viruses in other terrestrial ecosystems. In addition to cyanobacteria, predicted viral hosts include other lichen-associated bacterial lineages and algae, although a large fraction of viral contigs had no host prediction. Functional annotation of cyanolichen viral sequences predicts numerous viral-encoded auxiliary metabolic genes (AMGs) involved in amino acid, nucleotide, and carbohydrate metabolism, including AMGs for secondary metabolism (antibiotics and antimicrobials) and fatty acid biosynthesis. Overall, the diversity of cyanolichen AMGs suggests that viruses may alter microbial interactions within these complex symbiotic assemblages.

Differential microbial assemblages associated with shikonin-producing Borage species in two distinct soil types

Shikonin and its derivatives are the main components of traditional Chinese medicine, Zicao. The pharmacological potential of shikonin and its derivatives have been extensively studied. Yet, less is known about the microbial assemblages associated with shikonin producing Borage plants. We studied microbial profiles of two Borage species, Echium plantagineum (EP) and Lithospermum erythrorhizon (LE), to identify the dynamics of microbial colonization pattern within three rhizo-compatments and two distinct soil types. Results of α and β-diversity via PacBio sequencing revealed significantly higher microbial richness and diversity in the natural soil along with a decreasing microbial gradient across rhizosphere to endosphere. Our results displayed genotype and soil type-dependent fine-tuning of microbial profiles. The host plant was found to exert effects on the physical and chemical properties of soil, resulting in reproducibly different micro-biota. Analysis of differentially abundant microbial OTUs displayed Planctomycetes and Bacteroidetes to be specifically enriched in EP and LE rhizosphere while endosphere was mostly prevailed by Cyanobacteria. Network analysis to unfold co-existing microbial species displayed different types of positive and negative interactions within different communities. The data provided here will help to identify microbes associated with different rhizo-compartments of potential host plants. In the future, this might be helpful for manipulating the keystone microbes for ecosystem functioning.

Epigenetic biosensors for bacteriophage detection and phage receptor discrimination

Environmental monitoring of bacteria using phage-based biosensors has been widely developed for many different species. However, there are only a few available methods to detect specific bacteriophages in raw environmental samples. In this work, we developed a simple and efficient assay to rapidly monitor the phage content of a given sample. The assay is based on the bistable expression of the Salmonella enterica opvAB operon. Under regular growth conditions, opvAB is only expressed by a small fraction of the bacterial subpopulation. In the OpvABON subpopulation, synthesis of the OpvA and OpvB products shortens the O-antigen and confers resistance to phages that use LPS as a receptor. As a consequence, the OpvABON subpopulation is selected in the presence of such phages. Using an opvAB::gfp fusion, we could monitor LPS-binding phages in various media, including raw water samples. To enlarge our phage-biosensor panoply, we also developed biosensors able to detect LPS, as well as protein-binding coliphages. Moreover, the combination of these tools allowed to identify the bacterial receptor triggering phage infection. The epigenetic opvAB::gfp biosensor thus comes in different flavours to detect a wide range of bacteriophages and identify the type of receptor they recognize.

A New Alignment-Free Whole Metagenome Comparison Tool and Its Application on Gut Microbiomes of Wild Giant Pandas

The comparison of metagenomes is crucial for studying the relationship between microbial communities and environmental factors. One recently published alignment-free whole metagenome comparison method based on k-mer frequencies, Libra, showed higher resolutions than the present fastest method, Mash, on whole metagenomic sequencing reads, but it did not perform as well on the assembled contigs. Here, we developed a new alignment-free tool, KmerFreqCalc, for the comparison of the whole metagenomic data, which first calculated the frequencies of both forward and reverse complementary sequences of k-mers like Mash and then computed the cosine distance between the samples based on k-mer frequency vectors like Libra. We applied KmerFreqCalc on the assembled contigs of the gut microbiomes of wild giant pandas and compared the results to Libra and Mash. The results indicated that KmerFreqCalc was able to detect the subtle difference between giant panda samples caused by seasonal diet change, showing better clustering than Libra and Mash. Therefore, KmerFreqCalc has high resolution and accuracy for assembled contigs, being very suitable for comparison of samples with low dissimilarity.

Depth-related variability in viral communities in highly stratified sulfidic mine tailings

BACKGROUND

Recent studies have significantly expanded our knowledge of viral diversity and functions in the environment. Exploring the ecological relationships between viruses, hosts, and the environment is a crucial first step towards a deeper understanding of the complex and dynamic interplays among them.

RESULTS

Here, we obtained extensive 16S rRNA gene amplicon, metagenomics sequencing, and geochemical datasets from different depths of two highly stratified sulfidic mine tailings cores with steep geochemical gradients especially pH, and explored how variations in viral community composition and functions were coupled to the co-existing prokaryotic assemblages and the varying environmental conditions. Our data showed that many viruses in the mine tailings represented novel genera, based on gene-sharing networks. Siphoviridae, Podoviridae, and Myoviridae dominated the classified viruses in the surface tailings and deeper layers. Both viral richness and normalized coverage increased with depth in the tailings cores and were significantly correlated with geochemical properties, for example, pH. Viral richness was also coupled to prokaryotic richness (Pearson's r = 0.65, P = 0.032). The enrichment of prophages in the surface mine tailings suggested a preference of lysogenic viral lifestyle in more acidic conditions. Community-wide comparative analyses clearly showed that viruses in the surface tailings encoded genes mostly with unknown functions while viruses in the deeper layers contained genes mainly annotated as conventional functions related to metabolism and structure. Notably, significantly abundant assimilatory sulfate reduction genes were identified from the deeper tailings layers and they were widespread in viruses predicted to infect diverse bacterial phyla.

CONCLUSIONS

Overall, our results revealed a depth-related distribution of viral populations in the extreme and heterogeneous tailings system. The viruses may interact with diverse hosts and dynamic environmental conditions and likely play a role in the functioning of microbial community and modulate sulfur cycles in situ. Video Abstract.

Genetically similar temperate phages form coalitions with their shared host that lead to niche-specific fitness effects

Temperate phages engage in long-term associations with their hosts that may lead to mutually beneficial interactions, of which the full extent is presently unknown. Here, we describe an environmentally relevant model system with a single host, a species of the Roseobacter clade of marine bacteria, and two genetically similar phages (ɸ-A and ɸ-D). Superinfection of a ɸ-D lysogenized strain (CB-D) with ɸ-A particles resulted in a lytic infection, prophage induction, and conversion of a subset of the host population, leading to isolation of a newly ɸ-A lysogenized strain (CB-A). Phenotypic differences, predicted to result from divergent lysogenic-lytic switch mechanisms, are evident between these lysogens, with CB-A displaying a higher incidence of spontaneous induction. Doubling times of CB-D and CB-A in liquid culture are 75 and 100 min, respectively. As cell cultures enter stationary phase, CB-A viable counts are half of CB-D. Consistent with prior evidence that cell lysis enhances biofilm formation, CB-A produces twice as much biofilm biomass as CB-D. As strains are susceptible to infection by the opposing phage type, co-culture competitions were performed to test fitness effects. When grown planktonically, CB-A outcompeted CB-D three to one. Yet, during biofilm growth, CB-D outcompeted CB-A three to one. These results suggest that genetically similar phages can have divergent influence on the competitiveness of their shared hosts in distinct environmental niches, possibly due to a complex form of phage-mediated allelopathy. These findings have implications for enhanced understanding of the eco-evolutionary dynamics of host-phage interactions that are pervasive in all ecosystems.

Soil microbial diversity drops with land-use change in a high mountain temperate forest: a metagenomics survey

Land-use change has been identified as the most severe threat to biodiversity. Soils are important biodiversity reservoirs, but to what extent conversion of high-altitude temperate forest to arable land affects taxonomic and functional soil biodiversity is still largely unknown. Shotgun metagenomics was used to determine the taxonomic and functional diversity of bacteria, archaea and DNA virus in terms of effective number of species in high-altitude temperate oak and pine-oak forest and arable soils from Mexico. Generally, the soil ecosystem maintained its microbial species richness notwithstanding land-use change. Archaea diversity was not affected by land-use change, but the bacterial diversity decreased with 45-55% when the oak forest was converted to arable land and 65-75% when the pine-oak forest was. Loss in bacterial diversity as a result of land-use change was positively correlated (R² = 0.41) with the 10-25% loss in functional diversity. The archaeal communities were evener than the bacterial ones, which might explain their different response to land-use change. We expected a decrease in DNA viral communities as the bacterial diversity decreased, i.e. their potential hosts. However, a higher viral diversity was found in the arable than in the forest soils. It was found that converting high altitude oak and pine-oak forests to arable land more than halved the bacterial diversity, but did not affect the archaeal and even increased the viral diversity.

Genomic and ecological attributes of marine bacteriophages encoding bacterial virulence genes

Background

Bacteriophages encode genes that modify bacterial functions during infection. The acquisition of phage-encoded virulence genes is a major mechanism for the rise of bacterial pathogens. In coral reefs, high bacterial density and lysogeny has been proposed to exacerbate reef decline through the transfer of phage-encoded virulence genes. However, the functions and distribution of these genes in phage virions on the reef remain unknown.

Results

Here, over 28,000 assembled viral genomes from the free viral community in Atlantic and Pacific Ocean coral reefs were queried against a curated database of virulence genes. The diversity of virulence genes encoded in the viral genomes was tested for relationships with host taxonomy and bacterial density in the environment. These analyses showed that bacterial density predicted the profile of virulence genes encoded by phages. The Shannon diversity of virulence-encoding phages was negatively related with bacterial density, leading to dominance of fewer genes at high bacterial abundances. A statistical learning analysis showed that reefs with high microbial density were enriched in viruses encoding genes enabling bacterial recognition and invasion of metazoan epithelium. Over 60% of phages could not have their hosts identified due to limitations of host prediction tools; for those which hosts were identified, host taxonomy was not an indicator of the presence of virulence genes.

Conclusions

This study described bacterial virulence factors encoded in the genomes of bacteriophages at the community level. The results showed that the increase in microbial densities that occurs during coral reef degradation is associated with a change in the genomic repertoire of bacteriophages, specifically in the diversity and distribution of bacterial virulence genes. This suggests that phages are implicated in the rise of pathogens in disturbed marine ecosystems.

Temporal Variability of Virioplankton during a Gymnodinium catenatum Algal Bloom

Viruses are key biogeochemical engines in the regulation of the dynamics of phytoplankton. However, there has been little research on viral communities in relation to algal blooms. Using the virMine tool, we analyzed viral information from metagenomic data of field dinoflagellate (Gymnodinium catenatum) blooms at different stages. Species identification indicated that phages were the main species. Unifrac analysis showed clear temporal patterns in virioplankton dynamics. The viral community was dominated by Siphoviridae, Podoviridae, and Myoviridae throughout the whole bloom cycle. However, some changes were observed at different phases of the bloom; the relatively abundant Siphoviridae and Myoviridae dominated at pre-bloom and peak bloom stages, while at the post-bloom stage, the members of Phycodnaviridae and Microviridae were more abundant. Temperature and nutrients were the main contributors to the dynamic structure of the viral community. Some obvious correlations were found between dominant viral species and host biomass. Functional analysis indicated some functional genes had dramatic response in algal-associated viral assemblages, especially the CAZyme encoding genes. This work expands the existing knowledge of algal-associated viruses by characterizing viral composition and function across a complete algal bloom cycle. Our data provide supporting evidence that viruses participate in dinoflagellate bloom dynamics under natural conditions.

Viral Communities in the Global Deep Ocean Conveyor Belt Assessed by Targeted Viromics

Viruses are an abundant, diverse and dynamic component of marine and terrestrial ecosystems. In the ocean, viruses play a key role in the biogeochemical cycles and controlling microbial abundance, diversity and evolution. Recent metagenomics studies assessed the structure of the viral community in the upper ocean. However, little is known about the compositional changes in viral communities along the deep ocean conveyor belt. To assess potential changes in the viral community in the global deep-water circulation system, water samples were collected in the core of the North Atlantic Deep Water (NADW) (∼2,500 m) and Pacific Antarctic Bottom Water (∼4,000 m). Microbial and viral abundance were evaluated by flow cytometry. Subsequently, flow cytometry was used to sort virus-like particles and next generation sequencing was applied to build DNA libraries from the sorted virus populations. The viral communities were highly diverse across different oceanic regions with high dissimilarity between samples. Only 18% of the viral protein clusters were shared between the NADW and the Pacific Antarctic Bottom Water. Few viral groups, mainly associated with uncultured environmental and uncultured Mediterranean viruses were ubiquitously distributed along the global deep-water circulation system. Thus, our results point to a few groups of widely distributed abundant viruses in addition to the presence of rare and diverse types of viruses at a local scale.

The Promises and Pitfalls of Machine Learning for Detecting Viruses in Aquatic Metagenomes

Tools allowing for the identification of viral sequences in host-associated and environmental metagenomes allows for a better understanding of the genetics and ecology of viruses and their hosts. Recently, new approaches using machine learning methods to distinguish viral from bacterial signal using k-mer sequence signatures were published for identifying viral contigs in metagenomes. The promise of these content-based approaches is the ability to discover new viruses, with no or few known relatives. In this perspective paper, we examine the use of the content-based machine learning tool VirFinder for the identification of viral sequences in aquatic metagenomes and explore the possibility of using ecosystem-focused models targeted to marine metagenomes. We discuss the impact of the training set composition on the tool performance and the current limitation for the retrieval of low abundance viral sequences in metagenomes. We identify potential biases that could arise from machine learning approaches for viral hunting in real-world datasets and suggest possible avenues to overcome them.

Diversities and potential biogeochemical impacts of mangrove soil viruses

BACKGROUND

Mangroves are ecologically and economically important forests of the tropics. As one of the most carbon-rich biomes, mangroves account for 11% of the total input of terrestrial carbon into oceans. Although viruses are considered to significantly influence local and global biogeochemical cycles, little information is available regarding the community structure, genetic diversity and ecological roles of viruses in mangrove ecosystems.

METHODS

Here, we utilised viral metagenomics sequencing and virome-specific bioinformatics tools to study viral communities in six mangrove soil samples collected from different mangrove habitats in Southern China.

RESULTS

Mangrove soil viruses were found to be largely uncharacterised. Phylogenetic analyses of the major viral groups demonstrated extensive diversity and previously unknown viral clades and suggested that global mangrove viral communities possibly comprise evolutionarily close genotypes. Comparative analysis of viral genotypes revealed that mangrove soil viromes are mainly affected by marine waters, with less influence coming from freshwaters. Notably, we identified abundant auxiliary carbohydrate-active enzyme (CAZyme) genes from mangrove viruses, most of which participate in biolysis of complex polysaccharides, which are abundant in mangrove soils and organism debris. Host prediction results showed that viral CAZyme genes are diverse and probably widespread in mangrove soil phages infecting diverse bacteria of different phyla.

CONCLUSIONS

Our results showed that mangrove viruses are diverse and probably directly manipulate carbon cycling by participating in biomass recycling of complex polysaccharides, providing the knowledge essential in revealing the ecological roles of viruses in mangrove ecosystems.

Meta-network: optimized species-species network analysis for microbial communities

Background

The explosive growth of microbiome data provides ample opportunities to gain a better understanding of the microbes and their interactions in microbial communities. Given these massive data, optimized data mining methods become important and necessary to perform deep and comprehensive analysis. Among the various priorities for microbiome data mining, the examination of species-species co-occurrence patterns becomes one of the key themes in urgent need.

Results

Hence, in this work, we propose the Meta-Network framework to lucubrate the microbial communities. Rooted in loose definitions of network (two species co-exist in a certain samples rather than all samples) as well as association rule mining (mining more complex forms of correlations like indirect correlation and mutual information), this framework outperforms other methods in restoring the microbial communities, based on two cohorts of microbial communities: (a) the loose definition strategy is capable to generate more reasonable relationships among species in the species-species co-occurrence network; (b) important species-species co-occurrence patterns could not be identified by other existing approaches, but could successfully generated by association rule mining.

Conclusions

Results have shown that the species-species co-occurrence network we generated are much more informative than those based on traditional methods. Meta-Network has consistently constructed more meaningful networks with biologically important clusters, hubs, and provides a general approach towards deciphering the species-species co-occurrence networks.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5471-1) contains supplementary material, which is available to authorized users.

Host-hijacking and planktonic piracy: how phages command the microbial high seas

Microbial communities living in the oceans are major drivers of global biogeochemical cycles. With nutrients limited across vast swathes of the ocean, marine microbes eke out a living under constant assault from predatory viruses. Viral concentrations exceed those of their bacterial prey by an order of magnitude in surface water, making these obligate parasites the most abundant biological entities in the ocean. Like the pirates of the 17th and 18th centuries that hounded ships plying major trade and exploration routes, viruses have evolved mechanisms to hijack microbial cells and repurpose their cargo and indeed the vessels themselves to maximise viral propagation. Phenotypic reconfiguration of the host is often achieved through Auxiliary Metabolic Genes - genes originally derived from host genomes but maintained and adapted in viral genomes to redirect energy and substrates towards viral synthesis. In this review, we critically evaluate the literature describing the mechanisms used by bacteriophages to reconfigure host metabolism and to plunder intracellular resources to optimise viral production. We also highlight the mechanisms used when, in challenging environments, a 'batten down the hatches' strategy supersedes that of 'plunder and pillage'. Here, the infecting virus increases host fitness through phenotypic augmentation in order to ride out the metaphorical storm, with a concomitant impact on host substrate uptake and metabolism, and ultimately, their interactions with their wider microbial community. Thus, the traditional view of the virus-host relationship as predator and prey does not fully characterise the variety or significance of the interactions observed. Recent advances in viral metagenomics have provided a tantalising glimpse of novel mechanisms of viral metabolic reprogramming in global oceans. Incorporation of these new findings into global biogeochemical models requires experimental evidence from model systems and major improvements in our ability to accurately predict protein function from sequence data.

Libra: scalable k-mer-based tool for massive all-vs-all metagenome comparisons

Background

Shotgun metagenomics provides powerful insights into microbial community biodiversity and function. Yet, inferences from metagenomic studies are often limited by dataset size and complexity and are restricted by the availability and completeness of existing databases. De novo comparative metagenomics enables the comparison of metagenomes based on their total genetic content.

Results

We developed a tool called Libra that performs an all-vs-all comparison of metagenomes for precise clustering based on their k-mer content. Libra uses a scalable Hadoop framework for massive metagenome comparisons, Cosine Similarity for calculating the distance using sequence composition and abundance while normalizing for sequencing depth, and a web-based implementation in iMicrobe (http://imicrobe.us) that uses the CyVerse advanced cyberinfrastructure to promote broad use of the tool by the scientific community.

Conclusions

A comparison of Libra to equivalent tools using both simulated and real metagenomic datasets, ranging from 80 million to 4.2 billion reads, reveals that methods commonly implemented to reduce compute time for large datasets, such as data reduction, read count normalization, and presence/absence distance metrics, greatly diminish the resolution of large-scale comparative analyses. In contrast, Libra uses all of the reads to calculate k-mer abundance in a Hadoop architecture that can scale to any size dataset to enable global-scale analyses and link microbial signatures to biological processes.

The Promises and Pitfalls of Machine Learning for Detecting Viruses in Aquatic Metagenomes

Tools allowing for the identification of viral sequences in host-associated and environmental metagenomes allows for a better understanding of the genetics and ecology of viruses and their hosts. Recently, new approaches using machine learning methods to distinguish viral from bacterial signal using k-mer sequence signatures were published for identifying viral contigs in metagenomes. The promise of these content-based approaches is the ability to discover new viruses, with no or few known relatives. In this perspective paper, we examine the use of the content-based machine learning tool VirFinder for the identification of viral sequences in aquatic metagenomes and explore the possibility of using ecosystem-focused models targeted to marine metagenomes. We discuss the impact of the training set composition on the tool performance and the current limitation for the retrieval of low abundance viral sequences in metagenomes. We identify potential biases that could arise from machine learning approaches for viral hunting in real-world datasets and suggest possible avenues to overcome them.

Using noninvasive metagenomics to characterize viral communities from wildlife

Microbial communities play an important role in organismal and ecosystem health. While high-throughput metabarcoding has revolutionized the study of bacterial communities, generating comparable viral communities has proven elusive, particularly in wildlife samples where the diversity of viruses and limited quantities of viral nucleic acid present distinctive challenges. Metagenomic sequencing is a promising solution for studying viral communities, but the lack of standardized methods currently precludes comparisons across host taxa or localities. Here, we developed an untargeted shotgun metagenomic sequencing protocol to generate comparable viral communities from noninvasively collected faecal and oropharyngeal swabs. Using samples from common vampire bats (Desmodus rotundus), a key species for virus transmission to humans and domestic animals, we tested how different storage media, nucleic acid extraction procedures and enrichment steps affect viral community detection. Based on finding viral contamination in foetal bovine serum, we recommend storing swabs in RNAlater or another nonbiological medium. We recommend extracting nucleic acid directly from swabs rather than from supernatant or pelleted material, which had undetectable levels of viral RNA. Results from a low-input RNA library preparation protocol suggest that ribosomal RNA depletion and light DNase treatment reduce host and bacterial nucleic acid, and improve virus detection. Finally, applying our approach to twelve pooled samples from seven localities in Peru, we showed that detected viral communities saturated at the attained sequencing depth, allowing unbiased comparisons of viral community composition. Future studies using the methods outlined here will elucidate the determinants of viral communities across host species, environments and time.

Morphological characterization of virus-like particles in coral reef sponges

Marine sponges host complex microbial consortia that vary in their abundance, diversity and stability amongst host species. While our understanding of sponge-microbe interactions has dramatically increased over the past decade, little is known about how sponges and their microbial symbionts interact with viruses, the most abundant entities in the ocean. In this study, we employed three transmission electron microscopy (TEM) preparation methods to provide the first comprehensive morphological assessment of sponge-associated viruses. The combined approaches revealed 50 different morphologies of viral-like particles (VLPs) represented across the different sponge species. VLPs were visualized within sponge cells, within the sponge extracellular mesohyl matrix, on the sponge ectoderm and within sponge-associated microbes. Non-enveloped, non-tailed icosahedral VLPs were the most commonly observed morphotypes, although tailed bacteriophage, brick-shaped, geminate and filamentous VLPs were also detected. Visualization of sponge-associated viruses using TEM has confirmed that sponges harbor not only diverse communities of microorganisms but also diverse communities of viruses.

The North Sea goes viral: Occurrence and distribution of North Sea bacteriophages

Marine viruses are dominated by phages and have an enormous influence on microbial population dynamics, due to lysis and horizontal gene transfer. The aim of this study is to analyze the occurrence and diversity of phages in the North Sea, considering the virus-host interactions and biogeographic factors. The virus community of four sampling stations were described using virus metagenomics (viromes). The results show that the virus community was not evenly distributed throughout the North Sea. The dominant phage members were identified as unclassified phage group, followed by Caudovirales order. Myoviridae was the dominant phage family in the North Sea, which occurrence decreased from the coast to the open sea. In contrast, the occurrence of Podoviridae increased and the occurrence of Siphoviridae was low throughout the North Sea. The occurrence of other groups such as Phycodnaviridae decreased from the coast to the open sea. The coastal virus community was genetically more diverse than the open sea community. The influence of riverine inflow and currents, for instance the English Channel flow affects the genetic virus diversity with the community carrying genes from a variety of metabolic pathways and other functions. The present study offers the first insights in the virus community in the North Sea using viromes and shows the variation in virus diversity and the genetic information moved from coastal to open sea areas.

Soil Viruses Are Underexplored Players in Ecosystem Carbon Processing

Rapidly thawing permafrost harbors ∼30 to 50% of global soil carbon, and the fate of this carbon remains unknown. Microorganisms will play a central role in its fate, and their viruses could modulate that impact via induced mortality and metabolic controls. Because of the challenges of recovering viruses from soils, little is known about soil viruses or their role(s) in microbial biogeochemical cycling. Here, we describe 53 viral populations (viral operational taxonomic units [vOTUs]) recovered from seven quantitatively derived (i.e., not multiple-displacement-amplified) viral-particle metagenomes (viromes) along a permafrost thaw gradient at the Stordalen Mire field site in northern Sweden. Only 15% of these vOTUs had genetic similarity to publicly available viruses in the RefSeq database, and ∼30% of the genes could be annotated, supporting the concept of soils as reservoirs of substantial undescribed viral genetic diversity. The vOTUs exhibited distinct ecology, with different distributions along the thaw gradient habitats, and a shift from soil-virus-like assemblages in the dry palsas to aquatic-virus-like assemblages in the inundated fen. Seventeen vOTUs were linked to microbial hosts (in silico), implicating viruses in infecting abundant microbial lineages from Acidobacteria, Verrucomicrobia, and Deltaproteobacteria, including those encoding key biogeochemical functions such as organic matter degradation. Thirty auxiliary metabolic genes (AMGs) were identified and suggested virus-mediated modulation of central carbon metabolism, soil organic matter degradation, polysaccharide binding, and regulation of sporulation. Together, these findings suggest that these soil viruses have distinct ecology, impact host-mediated biogeochemistry, and likely impact ecosystem function in the rapidly changing Arctic. IMPORTANCE This work is part of a 10-year project to examine thawing permafrost peatlands and is the first virome-particle-based approach to characterize viruses in these systems. This method yielded >2-fold-more viral populations (vOTUs) per gigabase of metagenome than vOTUs derived from bulk-soil metagenomes from the same site (J. B. Emerson, S. Roux, J. R. Brum, B. Bolduc, et al., Nat Microbiol 3:870-880, 2018, https://doi.org/10.1038/s41564-018-0190-y). We compared the ecology of the recovered vOTUs along a permafrost thaw gradient and found (i) habitat specificity, (ii) a shift in viral community identity from soil-like to aquatic-like viruses, (iii) infection of dominant microbial hosts, and (iv) carriage of host metabolic genes. These vOTUs can impact ecosystem carbon processing via top-down (inferred from lysing dominant microbial hosts) and bottom-up (inferred from carriage of auxiliary metabolic genes) controls. This work serves as a foundation which future studies can build upon to increase our understanding of the soil virosphere and how viruses affect soil ecosystem services.

Genetic Diversity and Cooccurrence Patterns of Marine Cyanopodoviruses and Picocyanobacteria

Picocyanobacteria Prochlorococcus and Synechococcus are abundant in the global oceans and subject to active viral infection. In this study, the genetic diversity of picocyanobacteria and the genetic diversity of cyanopodoviruses were synchronously investigated along water columns in the equatorial Indian Ocean and over a seasonal time course in the coastal Sanya Bay, South China Sea. Using the 16S-23S rRNA internal transcribed spacer (ITS)-based clone library and quantitative PCR (qPCR) analyses, the picocyanobacterial community composition and abundance were determined. Sanya Bay was dominated by clade II Synechococcus during all the seasons, and a typical population shift from high-light-adapted Prochlorococcus to low-light-adapted Prochlorococcus was found along the vertical profiles. Strikingly, the DNA polymerase gene sequences of cyanopodoviruses revealed a much greater genetic diversity than we expected. Nearly one-third of the phylogenetic groups were newly described here. No apparent seasonal pattern was observed for the Sanya Bay picocyanobacterial or cyanopodoviral communities. Different dominant cyanopodovirus lineages were identified for the coastal area, upper euphotic zone, and middle-to-lower euphotic zone of the open ocean. Diversity indices of both picocyanobacteria and cyanopodoviruses were highest in the middle euphotic zone and both were lower in the upper euphotic zone, reflecting a host-virus interaction. Cyanopodoviral communities differed significantly between the upper euphotic zone and the middle-to-lower euphotic zone, showing a vertical pattern similar to that of picocyanobacteria. However, in the surface waters of the open ocean, cyanopodoviruses exhibited no apparent biogeographic pattern, differing from picocyanobacteria. This study demonstrates correlated distribution patterns of picocyanobacteria and cyanopodoviruses, as well as the complex biogeography of cyanopodoviruses.IMPORTANCE Picocyanobacteria are highly diverse and abundant in the ocean and display remarkable global biogeography and a vertical distribution pattern. However, how the diversity and distribution of picocyanobacteria affect those of the viruses that infect them remains largely unknown. Here we synchronously analyzed the community structures of cyanopodoviruses and picocyanobacteria at spatial and temporal scales. Both spatial and temporal variations of cyanopodoviral communities can be linked to those of picocyanobacteria. The coastal area, upper euphotic zone, and middle-to-lower euphotic zone of the open ocean have distinct cyanopodoviral communities, showing horizontal and vertical variation patterns closely related to those of picocyanobacteria. These findings emphasize the driving force of host community in shaping the biogeographic structure of viruses. Our work provides important information for future assessments of the ecological roles of viruses and hosts for each other.

Pseudomonas PB1-Like Phages: Whole Genomes from Metagenomes Offer Insight into an Abundant Group of Bacteriophages

Despite the abundance, ubiquity and impact of environmental viruses, their inherent genomic plasticity and extreme diversity pose significant challenges for the examination of bacteriophages on Earth. Viral metagenomic studies have offered insight into broader aspects of phage ecology and repeatedly uncover genes to which we are currently unable to assign function. A combined effort of phage isolation and metagenomic survey of Chicago’s nearshore waters of Lake Michigan revealed the presence of Pbunaviruses, relatives of the Pseudomonas phage PB1. This prompted our expansive investigation of PB1-like phages. Genomic signatures of PB1-like phages and Pbunaviruses were identified, permitting the unambiguous distinction between the presence/absence of these phages in soils, freshwater and wastewater samples, as well as publicly available viral metagenomic datasets. This bioinformatic analysis led to the de novo assembly of nine novel PB1-like phage genomes from a metagenomic survey of samples collected from Lake Michigan. While this study finds that Pbunaviruses are abundant in various environments of Northern Illinois, genomic variation also exists to a considerable extent within individual communities.

Mixing alters the lytic activity of viruses in the dark ocean

In aquatic habitats, viral lysis of prokaryotic cells lowers the overall efficiency of the microbial loop, by which dissolved organic carbon is transfered to higher trophic levels. Mixing of water masses in the dark ocean occurs on a global scale and may have far reaching consequences for the different prokaryotic and virus communities found in these waters by altering the environmental conditions these communities experience. We hypothesize that mixing of deep ocean water masses enhances the lytic activity of viruses infecting prokaryotes. To address this hypothesis, major deep-sea water masses of the Atlantic Ocean such as North Atlantic Deep Water, Mediterranean Sea Overflow Water, Antarctic Intermediate Water, and Antarctic Bottom Water were sampled at five locations. Prokaryotic cells from these samples were collected by filtration and subsequently incubated in virus-reduced water from either the same (control) or a different water mass (transplantation treatment). Additionally, mixtures of prokaryotes obtained from two different water masses were incubated in a mixture of virus-reduced water from the same water masses (control) or in virus-reduced water from the source water masses separately (mixing treatments). Pronounced differences in productivity-related parameters (prokaryotic leucine incorporation, prokaryotic and viral abundance) between water masses caused strong changes in viral lysis of prokaryotes. Often, mixing of water masses increased viral lysis of prokaryotes, indicating that lysogenic viruses were induced into the lytic cycle. Mixing-induced changes in viral lysis had a strong effect on the community composition of prokaryotes and viruses. Our data show that mixing of deep-sea water masses alters levels of viral lysis of prokaryotes and in many cases weakens the efficiency of the microbial loop by enhancing the recycling of organic carbon in the deep ocean.

Dietary energy drives the dynamic response of bovine rumen viral communities

BACKGROUND

Rumen microbes play a greater role in host energy acquisition than that of gut-associated microbes in monogastric animals. Although genome-enabled advancements are providing access to the vast diversity of uncultivated microbes, our understanding of variables shaping rumen microbial communities is in its infancy. Viruses have been shown to impact microbial populations through a myriad of processes, including cell lysis and reprogramming of host metabolism. However, little is known about the processes shaping the distribution of rumen viruses or how viruses may modulate microbial-driven processes in the rumen. To this end, we investigated how rumen bacterial and viral community structure and function responded in five steers fed four randomized dietary treatments in a crossover design.

RESULTS

Total digestible nutrients (TDN), a measure of dietary energy, best explained the variation in bacterial and viral communities. Additional ecological drivers of viral communities included dietary zinc content and microbial functional diversity. Using partial least squares regression, we demonstrate significant associations between the abundances of 267 viral populations and variables driving the variation in rumen viral communities. While rumen viruses were dynamic, 14 near ubiquitous viral populations were identified, suggesting the presence of a core rumen virome largely comprised of novel viruses. Moreover, analysis of virally encoded auxiliary metabolic genes (AMGs) indicates rumen viruses have glycosidic hydrolases to potentially augment the breakdown of complex carbohydrates to increase energy production. Other AMGs identified have a role in redirecting carbon to the pentose phosphate pathway and one carbon pools by folate to boost viral replication.

CONCLUSIONS

We demonstrate that rumen bacteria and viruses have differing responses and ecological drivers to dietary perturbation. Our results show that rumen viruses have implications for understanding the structuring of the previously identified core rumen microbiota and impacting microbial metabolism through a vast array of AMGs. AMGs in the rumen appear to have consequences for microbial metabolism that are largely in congruence with the current paradigm established in marine systems. This study provides a foundation for future hypotheses regarding the dynamics of viral-mediated processes in the rumen.

Phage hunters: Computational strategies for finding phages in large-scale 'omics datasets

A plethora of tools exist for identifying phage sequences in bacterial genomes, single cell amplified genomes, and host-associated and environmental metagenomes. Yet because the genetics of phages and their hosts are closely intertwined, distinguishing viral from bacterial signal remains an ongoing challenge. Further the size, quantity and fragmentary nature of modern 'omics datasets ushers in a new set of computational challenges. Here, we detail the promises and pitfalls of using currently available gene-centric or k-mer based tools for identifying prophage sequences in genomes and prophage and viral contigs in metagenomes. Each of these methods offers a unique piece of the puzzle to elucidating the intriguing signatures of phage-host coevolution.

Diel cycling and long-term persistence of viruses in the ocean's euphotic zone

Marine microbial communities exert a large influence on ocean ecosystem processes, and viruses in these communities play key roles in controlling microbial abundances, nutrient cycling, and productivity. We show here that dominant viruses in the open ocean persist for long time periods and that many appear tightly locked in coordinated diel oscillations with their bacterial hosts. The persistent structure of viral assemblages, as well as synchronized daily oscillations of viruses and hosts, are in part the result of the regular diurnal coupling of viral and host replication cycles. Collectively, our results suggest that viruses, as key components of marine ecosystems, are intrinsically synchronized with the daily rhythms of microbial community processes in the ocean’s photic zone.

Viruses are fundamental components of marine microbial communities that significantly influence oceanic productivity, biogeochemistry, and ecosystem processes. Despite their importance, the temporal activities and dynamics of viral assemblages in natural settings remain largely unexplored. Here we report the transcriptional activities and variability of dominant dsDNA viruses in the open ocean’s euphotic zone over daily and seasonal timescales. While dsDNA viruses exhibited some fluctuation in abundance in both cellular and viral size fractions, the viral assemblage was remarkably stable, with the most abundant viral types persisting over many days. More extended time series indicated that long-term persistence (>1 y) was the rule for most dsDNA viruses observed, suggesting that both core viral genomes as well as viral community structure were conserved over interannual periods. Viral gene transcription in host cell assemblages revealed diel cycling among many different viral types. Most notably, an afternoon peak in cyanophage transcriptional activity coincided with a peak in Prochlorococcus DNA replication, indicating coordinated diurnal coupling of virus and host reproduction. In aggregate, our analyses suggested a tightly synchronized diel coupling of viral and cellular replication cycles in both photoautotrophic and heterotrophic bacterial hosts. A surprising consequence of these findings is that diel cycles in the ocean’s photic zone appear to be universal organizing principles that shape ecosystem dynamics, ecological interactions, and biogeochemical cycling of both cellular and acellular community components.

Soil Viral Communities Vary Temporally and along a Land Use Transect as Revealed by Virus-Like Particle Counting and a Modified Community Fingerprinting Approach (fRAPD)

Environmental surveys on soil viruses are still rare and mostly anecdotal, i. e., they mostly report on viruses at one location or for only a few sampling dates. Detailed time-series analysis with multiple samples can reveal the spatio-temporal dynamics of viral communities and provide important input as to how viruses interact with their potential hosts and the environment. Such surveys, however, require fast, easy-to-apply and reliable methods. In the present study we surveyed monthly across 13 months the abundance of virus-like particles (VLP) and the structure of the viral communities in soils along a land use transect (i.e., forest, pasture, and cropland). We evaluated 32 procedures to extract VLP from soil using different buffers and mechanical methods. The most efficient extraction was achieved with 1× saline magnesium buffer in combination with 20 min vortexing. For community structure analysis we developed an optimized fingerprinting approach (fluorescent RAPD-PCR; fRAPD) by combining RAPD-PCR with fluorescently labeled primers in order to size the obtained fragments on a capillary sequencing machine. With the concomitantly collected data of soil specific factors and weather data, we were able to find correlations of viral abundance and community structure with environmental variables and sampling site. More specifically, we found that soil specific factors such as pH and total nitrogen content played a significant role in shaping both soil viral abundance and community structure. The fRAPD analysis revealed high temporal changes and clustered the viral communities according to sampling sites. In particular we observed that temperature and rainfall shaped soil viral communities in non-forest sites. In summary our findings suggest that sampling site was a key factor for shaping the abundance and community structure of soil viruses, and when site vegetation was reduced, temperature and rainfall were also important factors.

Virioplankton Assemblage Structure in the Lower River and Ocean Continuum of the Amazon

The Amazon River forms a vast plume in the Atlantic Ocean that can extend for more than 1,000 km. Microbial communities promote a globally relevant carbon sink system in the plume. Despite the importance of viruses for the global carbon cycle, the diversity and the possible roles of viruses in the Amazonia are poorly understood. The present work assesses, for the first time, the abundance and diversity of viruses simultaneously in the river and ocean in order to elucidate their possible roles. DNA sequence assembly yielded 29,358 scaffolds, encoding 82,546 viral proteins, with 15 new complete viral genomes from the 12 river and ocean locations. Viral diversity was clearly distinguished by river and ocean. Bacteriophages were the most abundant and occurred throughout the continuum. Viruses that infect eukaryotes were more abundant in the river, whereas phages appeared to have strong control over the host prokaryotic populations in the plume.

The Amazon River watershed and its associated plume comprise a vast continental and oceanic area. The microbial activities along this continuum contribute substantially to global carbon and nutrient cycling, and yet there is a dearth of information on the diversity, abundance, and possible roles of viruses in this globally important river. The aim of this study was to elucidate the diversity and structure of virus assemblages of the Amazon River-ocean continuum. Environmental viral DNA sequences were obtained for 12 locations along the river’s lower reach (n = 5) and plume (n = 7). Sequence assembly yielded 29,358 scaffolds, encoding 82,546 viral proteins, with 15 new complete viral genomes. Despite the spatial connectivity mediated by the river, virome analyses and physical-chemical water parameters clearly distinguished river and plume ecosystems. Bacteriophages were ubiquitous in the continuum and were more abundant in the transition region. Eukaryotic viruses occurred mostly in the river, while the plume had more viruses of autotrophic organisms (Prochlorococcus, Synechococcus) and heterotrophic bacteria (Pelagibacter). The viral families Microviridae and Myoviridae were the most abundant and occurred throughout the continuum. The major functions of the genes in the continuum involved viral structures and life cycles, and viruses from plume locations and Tapajós River showed the highest levels of functional diversity. The distribution patterns of the viral assemblages were defined not only by the occurrence of possible hosts but also by water physical and chemical parameters, especially salinity. The findings presented here help to improve understanding of the possible roles of viruses in the organic matter cycle along the river-ocean continuum.

IMPORTANCE The Amazon River forms a vast plume in the Atlantic Ocean that can extend for more than 1,000 km. Microbial communities promote a globally relevant carbon sink system in the plume. Despite the importance of viruses for the global carbon cycle, the diversity and the possible roles of viruses in the Amazon are poorly understood. The present work assesses, for the first time, the abundance and diversity of viruses simultaneously in the river and ocean in order to elucidate their possible roles. DNA sequence assembly yielded 29,358 scaffolds, encoding 82,546 viral proteins, with 15 new complete viral genomes from the 12 river and ocean locations. Viral diversity was clearly distinguished by river and ocean. Bacteriophages were the most abundant and occurred throughout the continuum. Viruses that infect eukaryotes were more abundant in the river, whereas phages appeared to have strong control over the host prokaryotic populations in the plume.

Drivers of interannual variability in virioplankton abundance at the coastal western Antarctic peninsula and the potential effects of climate change

An 8-year time-series in the Western Antarctic Peninsula (WAP) with an approximately weekly sampling frequency was used to elucidate changes in virioplankton abundance and their drivers in this climatically sensitive region. Virioplankton abundances at the coastal WAP show a pronounced seasonal cycle with interannual variability in the timing and magnitude of the summer maxima. Bacterioplankton abundance is the most influential driving factor of the virioplankton, and exhibit closely coupled dynamics. Sea ice cover and duration predetermine levels of phytoplankton stock and thus, influence virioplankton by dictating the substrates available to the bacterioplankton. However, variations in the composition of the phytoplankton community and particularly the prominence of Diatoms inferred from silicate drawdown, drive interannual differences in the magnitude of the virioplankton bloom; likely again mediated through changes in the bacterioplankton. Their findings suggest that future warming within the WAP will cause changes in sea ice that will influence viruses and their microbial hosts through changes in the timing, magnitude and composition of the phytoplankton bloom. Thus, the flow of matter and energy through the viral shunt may be decreased with consequences for the Antarctic food web and element cycling.

Deep-Sea Hydrothermal Vent Viruses Compensate for Microbial Metabolism in Virus-Host Interactions

Viruses are believed to be responsible for the mortality of host organisms. However, some recent investigations reveal that viruses may be essential for host survival. To date, it remains unclear whether viruses are beneficial or harmful to their hosts. To reveal the roles of viruses in the virus-host interactions, viromes and microbiomes of sediment samples from three deep-sea hydrothermal vents were explored in this study. To exclude the influence of exogenous DNAs on viromes, the virus particles were purified with nuclease (DNase I and RNase A) treatments and cesium chloride density gradient centrifugation. The metagenomic analysis of viromes without exogenous DNA contamination and microbiomes of vent samples indicated that viruses had compensation effects on the metabolisms of their host microorganisms. Viral genes not only participated in most of the microbial metabolic pathways but also formed branched pathways in microbial metabolisms, including pyrimidine metabolism; alanine, aspartate, and glutamate metabolism; nitrogen metabolism and assimilation pathways of the two-component system; selenocompound metabolism; aminoacyl-tRNA biosynthesis; and amino sugar and nucleotide sugar metabolism. As is well known, deep-sea hydrothermal vent ecosystems exist in relatively isolated environments which are barely influenced by other ecosystems. The metabolic compensation of hosts mediated by viruses might represent a very important aspect of virus-host interactions.

Viruses are the most abundant biological entities in the oceans and have very important roles in regulating microbial community structure and biogeochemical cycles. The relationship between virus and host microbes is broadly thought to be that of predator and prey. Viruses can lyse host cells to control microbial population sizes and affect community structures of hosts by killing specific microbes. However, viruses also influence their hosts through manipulation of bacterial metabolism. We found that viral genes not only participated in most microbial metabolic pathways but also formed branched pathways in microbial metabolisms. The metabolic compensation of hosts mediated by viruses may help hosts to adapt to extreme environments and may be essential for host survival.

Marine viruses discovered via metagenomics shed light on viral strategies throughout the oceans

Marine viruses are key drivers of host diversity, population dynamics and biogeochemical cycling and contribute to the daily flux of billions of tons of organic matter. Despite recent advancements in metagenomics, much of their biodiversity remains uncharacterized. Here we report a data set of 27,346 marine virome contigs that includes 44 complete genomes. These outnumber all currently known phage genomes in marine habitats and include members of previously uncharacterized lineages. We designed a new method for host prediction based on co-occurrence associations that reveals these viruses infect dominant members of the marine microbiome such as Prochlorococcus and Pelagibacter. A negative association between host abundance and the virus-to-host ratio supports the recently proposed Piggyback-the-Winner model of reduced phage lysis at higher host densities. An analysis of the abundance patterns of viruses throughout the oceans revealed how marine viral communities adapt to various seasonal, temperature and photic regimes according to targeted hosts and the diversity of auxiliary metabolic genes.

Metagenomic Approaches to Assess Bacteriophages in Various Environmental Niches

Bacteriophages are ubiquitous and numerous parasites of bacteria and play a critical evolutionary role in virtually every ecosystem, yet our understanding of the extent of the diversity and role of phages remains inadequate for many ecological niches, particularly in cases in which the host is unculturable. During the past 15 years, the emergence of the field of viral metagenomics has drastically enhanced our ability to analyse the so-called viral ‘dark matter’ of the biosphere. Here, we review the evolution of viral metagenomic methodologies, as well as providing an overview of some of the most significant applications and findings in this field of research.

Intriguing Interaction of Bacteriophage-Host Association: An Understanding in the Era of Omics

Innovations in next-generation sequencing technology have introduced new avenues in microbial studies through “omics” approaches. This technology has considerably augmented the knowledge of the microbial world without isolation prior to their identification. With an enormous volume of bacterial “omics” data, considerable attempts have been recently invested to improve an insight into virosphere. The interplay between bacteriophages and their host has created a significant influence on the biogeochemical cycles, microbial diversity, and bacterial population regulation. This review highlights various concepts such as genomics, transcriptomics, proteomics, and metabolomics to infer the phylogenetic affiliation and function of bacteriophages and their impact on diverse microbial communities. Omics technologies illuminate the role of bacteriophage in an environment, the influences of phage proteins on the bacterial host and provide information about the genes important for interaction with bacteria. These investigations will reveal some of bio-molecules and biomarkers of the novel phage which demand to be unveiled.