Oxygen minimum zones harbour novel viral communities with low diversity

Physiological and metabolic insights into the first cultured anaerobic representative of deep-sea Planctomycetes bacteria

Planctomycetes bacteria are ubiquitously distributed across various biospheres and play key roles in global element cycles. However, few deep-sea Planctomycetes members have been cultivated, limiting our understanding of Planctomycetes in the deep biosphere. Here, we have successfully cultured a novel strain of Planctomycetes (strain ZRK32) from a deep-sea cold seep sediment. Our genomic, physiological, and phylogenetic analyses indicate that strain ZRK32 is a novel species, which we propose be named: Poriferisphaera heterotrophicis. We show that strain ZRK32 replicates using a budding mode of division. Based on the combined results from growth assays and transcriptomic analyses, we found that rich nutrients, or supplementation with NO3- or NH4+ promoted the growth of strain ZRK32 by facilitating energy production through the tricarboxylic acid cycle and the Embden-Meyerhof-Parnas glycolysis pathway. Moreover, supplementation with NO3- or NH4+ induced strain ZRK32 to release a bacteriophage in a chronic manner, without host cell lysis. This bacteriophage then enabled strain ZRK32, and another marine bacterium that we studied, to metabolize nitrogen through the function of auxiliary metabolic genes. Overall, these findings expand our understanding of deep-sea Planctomycetes bacteria, while highlighting their ability to metabolize nitrogen when reprogrammed by chronic viruses.

ViWrap: A modular pipeline to identify, bin, classify, and predict viral-host relationships for viruses from metagenomes

Viruses are increasingly being recognized as important components of human and environmental microbiomes. However, viruses in microbiomes remain difficult to study because of the difficulty in culturing them and the lack of sufficient model systems. As a result, computational methods for identifying and analyzing uncultivated viral genomes from metagenomes have attracted significant attention. Such bioinformatics approaches facilitate the screening of viruses from enormous sequencing datasets originating from various environments. Though many tools and databases have been developed for advancing the study of viruses from metagenomes, there is a lack of integrated tools enabling a comprehensive workflow and analyses platform encompassing all the diverse segments of virus studies. Here, we developed ViWrap, a modular pipeline written in Python. ViWrap combines the power of multiple tools into a single platform to enable various steps of virus analysis, including identification, annotation, genome binning, species- and genus-level clustering, assignment of taxonomy, prediction of hosts, characterization of genome quality, comprehensive summaries, and intuitive visualization of results. Overall, ViWrap enables a standardized and reproducible pipeline for both extensive and stringent characterization of viruses from metagenomes, viromes, and microbial genomes. Our approach has flexibility in using various options for diverse applications and scenarios, and its modular structure can be easily amended with additional functions as necessary. ViWrap is designed to be easily and widely used to study viruses in human and environmental systems. ViWrap is publicly available via GitHub (https://github.com/AnantharamanLab/ViWrap). A detailed description of the software, its usage, and interpretation of results can be found on the website.

ViWrap: A modular pipeline to identify, bin, classify, and predict viral-host relationships for viruses from metagenomes

Viruses are increasingly being recognized as important components of human and environmental microbiomes. However, viruses in microbiomes remain difficult to study because of difficulty in culturing them and the lack of sufficient model systems. As a result, computational methods for identifying and analyzing uncultivated viral genomes from metagenomes have attracted significant attention. Such bioinformatics approaches facilitate screening of viruses from enormous sequencing datasets originating from various environments. Though many tools and databases have been developed for advancing the study of viruses from metagenomes, there is a lack of integrated tools enabling a comprehensive workflow and analyses platform encompassing all the diverse segments of virus studies. Here, we developed ViWrap, a modular pipeline written in Python. ViWrap combines the power of multiple tools into a single platform to enable various steps of virus analysis including identification, annotation, genome binning, species- and genus-level clustering, assignment of taxonomy, prediction of hosts, characterization of genome quality, comprehensive summaries, and intuitive visualization of results. Overall, ViWrap enables a standardized and reproducible pipeline for both extensive and stringent characterization of viruses from metagenomes, viromes, and microbial genomes. Our approach has flexibility in using various options for diverse applications and scenarios, and its modular structure can be easily amended with additional functions as necessary. ViWrap is designed to be easily and widely used to study viruses in human and environmental systems. ViWrap is publicly available via GitHub (https://github.com/AnantharamanLab/ViWrap). A detailed description of the software, its usage, and interpretation of results can be found on the website.

Biogeochemical sulfur cycling of virus auxiliary metabolic genes involved in Napahai plateau wetland

Virus plays important roles in regulating microbial community structure, horizontal gene transfer, and promoting biological evolution, also augmenting host metabolism during infection via the expression of auxiliary metabolic genes (AMGs), and thus affect biogeochemical cycling in the oceans. As the "kidney of the earth," wetlands have rich biodiversity and abundant resources. Based on metagenomic data, 10 AMGs associated with sulfur cycling, i.e., tusA, moaD, dsrE, soxA, soxB, soxC, soxD, soxX, soxY, and soxZ, were analyzed in Napahai plateau wetland. The phylogenetic trees of AMGs involved in sulfur metabolism from different habitats and host origins were constructed. Combined with principal coordinate analysis, it revealed that most AMGs associated with sulfur metabolism clustered separately, indicating the abundance and uniqueness in this region. The sulfur metabolism pathways involved by AMGs were mainly SOX systems, among which sulfur oxidation was associated with moaD and dsrE genes, while sulfur transport was related to tusA genes. It provides an insight into the biogeochemical sulfur cycling in plateau wetlands and lays the foundation for further study on the co-evolution of virus and host.

Virus-to-prokaryote ratio in the Salar de Huasco and different ecosystems of the Southern hemisphere and its relationship with physicochemical and biological parameters

The virus-to-prokaryote ratio (VPR) has been used in many ecosystems to study the relationship between viruses and their hosts. While high VPR values indicate a high rate of prokaryotes' cell lysis, low values are interpreted as a decrease in or absence of viral activity. Salar de Huasco is a high-altitude wetland characterized by a rich microbial diversity associated with aquatic sites like springs, ponds, streams and a lagoon with variable physicochemical conditions. Samples from two ponds, Poza Rosada (PR) and Poza Verde (PV), were analyzed by epifluorescence microscopy to determine variability of viral and prokaryotic abundance and to calculate the VPR in a dry season. In addition, to put Salar de Huasco results into perspective, a compilation of research articles on viral and prokaryotic abundance, VPR, and metadata from various Southern hemisphere ecosystems was revised. The ecosystems were grouped into six categories: high-altitude wetlands, Pacific, Atlantic, Indian, and Southern Oceans and Antarctic lakes. Salar de Huasco ponds recorded similar VPR values (an average of 7.4 and 1.7 at PR and PV, respectively), ranging from 3.22 to 15.99 in PR. The VPR variability was associated with VA and chlorophyll a, when considering all data available for this ecosystem. In general, high-altitude wetlands recorded the highest VPR average (53.22 ± 95.09), followed by the Oceans, Southern (21.91 ± 25.72), Atlantic (19.57 ± 15.77) and Indian (13.43 ± 16.12), then Antarctic lakes (11.37 ± 15.82) and the Pacific Ocean (6.34 ± 3.79). Physicochemical variables, i.e., temperature, conductivity, nutrients (nitrate, ammonium, and phosphate) and chlorophyll a as a biological variable, were found to drive the VPR in the ecosystems analyzed. Thus, the viral activity in the Wetland followed similar trends of previous reports based on larger sets of metadata analyses. In total, this study highlights the importance of including viruses as a biological variable to study microbial temporal dynamics in wetlands considering their crucial role in the carbon budgets of these understudied ecosystems in the southern hemisphere.

A Synthesis of Viral Contribution to Marine Nitrogen Cycling

Nitrogen is an essential component of major cellular macromolecules, such as DNA and proteins. Its bioavailability has a fundamental influence on the primary production of both terrestrial and oceanic ecosystems. Diverse marine microbes consume nitrogen, while only a limited taxon could replenish it, leaving nitrogen one of the most deficient nutrients in the ocean. A variety of microbes are involved in complex biogeochemical transformations of nitrogen compounds, and their ecological functions might be regulated by viruses in different manners. First and foremost, viruses drive marine nitrogen flow via host cell lysis, releasing abundant organic nitrogen into the surrounding environment. Secondly, viruses can also participate in the marine nitrogen cycle by expressing auxiliary metabolic genes (AMGs) to modulate host nitrogen metabolic pathways, such as nitrification, denitrification, anammox, and nitrogen transmembrane transport. Additionally, viruses also serve as a considerable reservoir of nitrogen element. The efficient turnover of viruses fundamentally promotes nitrogen flow in the oceans. In this review, we summarize viral contributions in the marine nitrogen cycling in different aspects and discuss challenges and issues based on recent discoveries of novel viruses involved in different processes of nitrogen biotransformation.

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.

Microbial Ecology of Oxygen Minimum Zones Amidst Ocean Deoxygenation

Oxygen minimum zones (OMZs) have substantial effects on the global ecology and biogeochemical processes of marine microbes. However, the diversity and activity of OMZ microbes and their trophic interactions are only starting to be documented, especially in regard to the potential roles of viruses and protists. OMZs have expanded over the past 60 years and are predicted to expand due to anthropogenic climate change, furthering the need to understand these regions. This review summarizes the current knowledge of OMZ formation, the biotic and abiotic factors involved in OMZ expansion, and the microbial ecology of OMZs, emphasizing the importance of bacteria, archaea, viruses, and protists. We describe the recognized roles of OMZ microbes in carbon, nitrogen, and sulfur cycling, the potential of viruses in altering host metabolisms involved in these cycles, and the control of microbial populations by grazers and viruses. Further, we highlight the microbial community composition and roles of these organisms in oxic and anoxic depths within the water column and how these differences potentially inform how microbial communities will respond to deoxygenation. Additionally, the current literature on the alteration of microbial communities by other key climate change parameters such as temperature and pH are considered regarding how OMZ microbes might respond to these pressures. Finally, we discuss what knowledge gaps are present in understanding OMZ microbial communities and propose directions that will begin to close these gaps.

Ecology of inorganic sulfur auxiliary metabolism in widespread bacteriophages

Microbial sulfur metabolism contributes to biogeochemical cycling on global scales. Sulfur metabolizing microbes are infected by phages that can encode auxiliary metabolic genes (AMGs) to alter sulfur metabolism within host cells but remain poorly characterized. Here we identified 191 phages derived from twelve environments that encoded 227 AMGs for oxidation of sulfur and thiosulfate (dsrA, dsrC/tusE, soxC, soxD and soxYZ). Evidence for retention of AMGs during niche-differentiation of diverse phage populations provided evidence that auxiliary metabolism imparts measurable fitness benefits to phages with ramifications for ecosystem biogeochemistry. Gene abundance and expression profiles of AMGs suggested significant contributions by phages to sulfur and thiosulfate oxidation in freshwater lakes and oceans, and a sensitive response to changing sulfur concentrations in hydrothermal environments. Overall, our study provides fundamental insights on the distribution, diversity, and ecology of phage auxiliary metabolism associated with sulfur and reinforces the necessity of incorporating viral contributions into biogeochemical configurations.

The bacterial sulfur cycle in expanding dysoxic and euxinic marine waters

Dysoxic marine waters (DMW, < 1 μM oxygen) are currently expanding in volume in the oceans, which has biogeochemical, ecological and societal consequences on a global scale. In these environments, distinct bacteria drive an active sulfur cycle, which has only recently been recognized for open-ocean DMW. This review summarizes the current knowledge on these sulfur-cycling bacteria. Critical bottlenecks and questions for future research are specifically addressed. Sulfate-reducing bacteria (SRB) are core members of DMW. However, their roles are not entirely clear, and they remain largely uncultured. We found support for their remarkable diversity and taxonomic novelty by mining metagenome-assembled genomes from the Black Sea as model ecosystem. We highlight recent insights into the metabolism of key sulfur-oxidizing SUP05 and Sulfurimonas bacteria, and discuss the probable involvement of uncultivated SAR324 and BS-GSO2 bacteria in sulfur oxidation. Uncultivated Marinimicrobia bacteria with a presumed organoheterotrophic metabolism are abundant in DMW. Like SRB, they may use specific molybdoenzymes to conserve energy from the oxidation, reduction or disproportionation of sulfur cycle intermediates such as S0 and thiosulfate, produced from the oxidation of sulfide. We expect that tailored sampling methods and a renewed focus on cultivation will yield deeper insight into sulfur-cycling bacteria in DMW.

Potential virus-mediated nitrogen cycling in oxygen-depleted oceanic waters

Viruses play an important role in the ecology and biogeochemistry of marine ecosystems. Beyond mortality and gene transfer, viruses can reprogram microbial metabolism during infection by expressing auxiliary metabolic genes (AMGs) involved in photosynthesis, central carbon metabolism, and nutrient cycling. While previous studies have focused on AMG diversity in the sunlit and dark ocean, less is known about the role of viruses in shaping metabolic networks along redox gradients associated with marine oxygen minimum zones (OMZs). Here, we analyzed relatively quantitative viral metagenomic datasets that profiled the oxygen gradient across Eastern Tropical South Pacific (ETSP) OMZ waters, assessing whether OMZ viruses might impact nitrogen (N) cycling via AMGs. Identified viral genomes encoded six N-cycle AMGs associated with denitrification, nitrification, assimilatory nitrate reduction, and nitrite transport. The majority of these AMGs (80%) were identified in T4-like Myoviridae phages, predicted to infect Cyanobacteria and Proteobacteria, or in unclassified archaeal viruses predicted to infect Thaumarchaeota. Four AMGs were exclusive to anoxic waters and had distributions that paralleled homologous microbial genes. Together, these findings suggest viruses modulate N-cycling processes within the ETSP OMZ and may contribute to nitrogen loss throughout the global oceans thus providing a baseline for their inclusion in the ecosystem and geochemical models.

Genome-resolved viral ecology in a marine oxygen minimum zone

Oxygen minimum zones (OMZs) are critical to marine nitrogen cycling and global climate change. While OMZ microbial communities are relatively well-studied, little is known about their viruses. Here, we assess the viral community ecology of 22 deeply sequenced viral metagenomes along a gradient of oxygenated to anoxic waters (<0.02 μmol/l O₂ ) in the Eastern Tropical South Pacific (ETSP) OMZ. We identified 46 127 viral populations (≥5 kb), which augments the known viruses from ETSP by 10-fold. Viral communities clustered into six groups that correspond to oceanographic features. Oxygen concentration was the predominant environmental feature driving viral community structure. Alpha and beta diversity of viral communities in the anoxic zone were lower than in surface waters, which parallels the low microbial diversity seen in other studies. ETSP viruses were largely endemic, with the majority of shared viruses (87%) also present in other OMZ samples. We detected 543 putative viral-encoded auxiliary metabolic genes (AMGs), of which some have a distribution that reflects physico-chemical characteristics across depth. Together these findings provide an ecological baseline for viral community structure, drivers and population variability in OMZs that will help future studies assess the role of viruses in these climate-critical environments.

Viral elements and their potential influence on microbial processes along the permanently stratified Cariaco Basin redoxcline

Little is known about viruses in oxygen-deficient water columns (ODWCs). In surface ocean waters, viruses are known to act as gene vectors among susceptible hosts. Some of these genes may have metabolic functions and are thus termed auxiliary metabolic genes (AMGs). AMGs introduced to new hosts by viruses can enhance viral replication and/or potentially affect biogeochemical cycles by modulating key microbial pathways. Here we identify 748 viral populations that cluster into 94 genera along a vertical geochemical gradient in the Cariaco Basin, a permanently stratified and euxinic ocean basin. The viral communities in this ODWC appear to be relatively novel as 80 of these viral genera contained no reference viral sequences, likely due to the isolation and unique features of this system. We identify viral elements that encode AMGs implicated in distinctive processes, such as sulfur cycling, acetate fermentation, signal transduction, [Fe–S] formation, and N-glycosylation. These AMG-encoding viruses include two putative Mu-like viruses, and viral-like regions that may constitute degraded prophages that have been modified by transposable elements. Our results provide an insight into the ecological and biogeochemical impact of viruses oxygen-depleted and euxinic habitats.

A Review on Viral Metagenomics in Extreme Environments

Viruses are the most abundant biological entities in the biosphere, and have the ability to infect Bacteria, Archaea, and Eukaryotes. The virome is estimated to be at least ten times more abundant than the microbiome with 10⁷ viruses per milliliter and 10⁹ viral particles per gram in marine waters and sediments or soils, respectively. Viruses represent a largely unexplored genetic diversity, having an important role in the genomic plasticity of their hosts. Moreover, they also play a significant role in the dynamics of microbial populations. In recent years, metagenomic approaches have gained increasing popularity in the study of environmental viromes, offering the possibility of extending our knowledge related to both virus diversity and their functional characterization. Extreme environments represent an interesting source of both microbiota and their virome due to their particular physicochemical conditions, such as very high or very low temperatures and >1 atm hydrostatic pressures, among others. Despite the fact that some progress has been made in our understanding of the ecology of the microbiota in these habitats, few metagenomic studies have described the viromes present in extreme ecosystems. Thus, limited advances have been made in our understanding of the virus community structure in extremophilic ecosystems, as well as in their biotechnological potential. In this review, we critically analyze recent progress in metagenomic based approaches to explore the viromes in extreme environments and we discuss the potential for new discoveries, as well as methodological challenges and perspectives.

Microbial niches in marine oxygen minimum zones

In the ocean's major oxygen minimum zones (OMZs), oxygen is effectively absent from sea water and life is dominated by microorganisms that use chemicals other than oxygen for respiration. Recent studies that combine advanced genomic and chemical detection methods are delineating the different metabolic niches that microorganisms can occupy in OMZs. Understanding these niches, the microorganisms that inhabit them, and their influence on marine biogeochemical cycles is crucial as OMZs expand with increasing seawater temperatures.

Taking Advantage of the Genomics Revolution for Monitoring and Conservation of Chondrichthyan Populations

Chondrichthyes (sharks, rays, skates and chimaeras) are among the oldest extant predators and are vital to top-down regulation of oceanic ecosystems. They are an ecologically diverse group occupying a wide range of habitats and are thus, exploited by coastal, pelagic and deep-water fishing industries. Chondrichthyes are among the most data deficient vertebrate species groups making design and implementation of regulatory and conservation measures challenging. High-throughput sequencing technologies have significantly propelled ecological investigations and understanding of marine and terrestrial species’ populations, but there remains a paucity of NGS based research on chondrichthyan populations. We present a brief review of current methods to access genomic and metagenomic data from Chondrichthyes and discuss applications of these datasets to increase our understanding of chondrichthyan taxonomy, evolution, ecology and population structures. Last, we consider opportunities and challenges offered by genomic studies for conservation and management of chondrichthyan populations.

Metagenomic Characterization of the Viral Community of the South Scotia Ridge

Viruses are the most abundant biological entities in aquatic ecosystems and harbor an enormous amount of genetic diversity. Whereas their influence on marine ecosystems is widely acknowledged, current information about their diversity remains limited. We conducted a viral metagenomic analysis of water samples collected during the austral summer of 2016 from the South Scotia Ridge (SSR), near the Antarctic Peninsula. The taxonomic composition and diversity of the viral communities were investigated, and a functional assessment of the sequences was performed. Phylotypic analysis showed that most viruses belonged to the order Caudovirales, especially the family Podoviridae (41.92⁻48.7%), which is similar to the situation in the Pacific Ocean. Functional analysis revealed a relatively high frequency of phage-associated and metabolism genes. Phylogenetic analyses of phage TerL and Capsid_NCLDV (nucleocytoplasmic large DNA viruses) marker genes indicated that many sequences associated with Caudovirales and NCLDV were novel and distinct from known phage genomes. High Phaeocystis globosa virus virophage (Pgvv) signatures were found and complete and partial Pgvv-like were obtained, which influence host⁻virus interactions. Our study expands existing knowledge of viral communities and their diversities from the Antarctic region and provides basic data for further exploring polar microbiomes.

Genetic and functional diversity of double-stranded DNA viruses in a tropical monsoonal estuary, India

The present study illustrates the genetic diversity of four uncultured viral communities from the surface waters of Cochin Estuary (CE), India. Viral diversity inferred using Illumina HiSeq paired-end sequencing using a linker-amplified shotgun library (LASL) revealed different double-stranded DNA (dsDNA) viral communities. The water samples were collected from four stations PR1, PR2, PR3, and PR4, during the pre-monsoon (PRM) season. Analysis of virus families indicated that the Myoviridae was the most common viral community in the CE followed by Siphoviridae and Podoviridae. There were significant (p < 0.05) spatial variations in the relative abundance of dominant families in response to the salinity regimes. The relative abundance of Myoviridae and Podoviridae were high in the euryhaline region and Siphoviridae in the mesohaline region of the estuary. The predominant phage type in CE was phages that infected Synechococcus. The viral proteins were found to be involved in major functional activities such as ATP binding, DNA binding, and DNA replication. The study highlights the genetic diversity of dsDNA viral communities and their functional protein predictions from a highly productive estuarine system. Further, the metavirome data generated in this study will enhance the repertoire of publicly available dataset and advance our understanding of estuarine viral ecology.

High Incidence of Lysogeny in the Oxygen Minimum Zones of the Arabian Sea (Southwest Coast of India)

Though microbial processes in the oxygen minimum zones (OMZs) of the Arabian Sea (AS) are well documented, prokaryote-virus interactions are less known. The present study was carried out to determine the potential physico-chemical factors influencing viral abundances and their life strategies (lytic and lysogenic) along the vertical gradient in the OMZ of the AS (southwest coast of India). Water samples were collected during the southwest monsoon (SWM) season in two consecutive years (2015 and 2016) from different depths, namely, the surface layer, secondary chlorophyll a maxima (~30–40 m), oxycline (~70–80 m), and hypoxic/suboxic layers (~200–350 m). The high viral abundances observed in oxygenated surface waters (mean ± SD = 6.1 ± 3.4 × 10⁶ viral-like particles (VLPs) mL⁻¹), drastically decreased with depth in the oxycline region (1.2 ± 0.5 × 10⁶ VLPs mL⁻¹) and hypoxic/suboxic waters (0.3 ± 0.3 × 10⁶ VLPs mL⁻¹). Virus to prokaryote ratio fluctuated in the mixed layer (~10) and declined significantly (p < 0.001) to 1 in the hypoxic layer. Viral production (VP) and frequency of virus infected cells (FIC) were maximum in the surface and minimum in the oxycline layer, whereas the viral lysis was undetectable in the suboxic/hypoxic layer. The detection of a high percentage of lysogeny in suboxic (48%) and oxycline zones (9–24%), accompanied by undetectable rates of lytic viral infection support the hypothesis that lysogeny may represent the major survival strategy for viruses in unproductive or harsh nutrient/host conditions in deoxygenated waters.

Contrasting distributions of bacteriophages and eukaryotic viruses from contaminated coastal sediments

Viruses are ubiquitous, abundant and play an important role in all ecosystems. Here, we advance understanding of coastal sediment viruses by exploring links in the composition and abundance of sediment viromes to environmental stressors and sediment bacterial communities. We collected sediment from contaminated and reference sites in Sydney Harbour and used metagenomics to analyse viral community composition. The proportion of phages at contaminated sites was significantly greater than phages at reference sites, whereas eukaryotic viruses were relatively more abundant at reference sites. We observed shifts in viral and bacterial composition between contaminated and reference sites of a similar magnitude. Models based on sediment characteristics revealed that total organic carbon in the sediments explained most of the environmental stress-related variation in the viral dataset. Our results suggest that the presence of anthropogenic contaminants in coastal sediments could be influencing viral community composition with potential consequences for associated hosts and the environment.

Geospatial distribution of viromes in tropical freshwater ecosystems

This study seeks to understand the general distribution of virome abundance and diversity in tropical freshwater ecosystems in Singapore and the geospatial distribution of the virome under different landuse patterns. Correlations between diversity, environmental parameters and land use patterns were analyzed and significant correlations were highlighted. Overall, the majority (65.5%) of the annotated virome belonged to bacteriophages. The percentage of Caudovirales was higher in reservoirs whereas the percentages of Dicistroviridae, Microviridae and Circoviridae were higher in tributaries. Reservoirs showed a higher Shannon-index virome diversity compared to upstream tributaries. Land use (urbanized, agriculture and parkland areas) influenced the characteristics of the virome distribution pattern. Dicistroviridae and Microviridae were enriched in urbanized tributaries while Mimiviridae, Phycodnaviridae, Siphoviridae and Podoviridae were enriched in parkland reservoirs. Several sequences closely related to the emerging zoonotic virus, cyclovirus, and the human-related virus (human picobirnavirus), were also detected. In addition, the relative abundance of PMMoV (pepper mild mottle virus) sequences was significantly correlated with RT-qPCR measurements (0.588 < r < 0.879, p < 0.05). This study shows that spatial factors (e.g., reservoirs/tributaries, land use) are the main drivers of the viral community structure in tropical freshwater ecosystems.

Diversity and Distribution Characteristics of Viruses in Soils of a Marine-Terrestrial Ecotone in East China

A substantial gap remains in our understanding of the abundance, diversity, and ecology of viruses in soil although some advances have been achieved in recent years. In this study, four soil samples according to the salinity gradient from shore to inland in East China have been characterized. Results showed that spherical virus particles represented the largest viral component in all of the four samples. The viromes had remarkably different taxonomic compositions, and most of the sequences were derived from single-stranded DNA viruses, especially from families Microviridae and Circoviridae. Compared with viromes from other aquatic and sediment samples, the community compositions of our four soil viromes resembled each other, meanwhile coastal sample virome closely congregated with sediment and hypersaline viromes, and high salinity paddy soil sample virome was similar with surface sediment virome. Phylogenetic analysis of functional genes showed that four viromes have high diversity of the subfamily Gokushovirinae in family Microviridae and most of Circoviridae replicase protein sequences grouped within the CRESS-DNA viruses. This work provided an initial outline of the viral communities in marine-terrestrial ecotone and will improve our understanding of the ecological functions of soil viruses.

Phage Genome Annotation Using the RAST Pipeline

Phages are complex biomolecular machineries that have to survive in a bacterial world. Phage genomes show many adaptations to their lifestyle such as shorter genes, reduced capacity for redundant DNA sequences, and the inclusion of tRNAs in their genomes. In addition, phages are not free-living, they require a host for replication and survival. These unique adaptations provide challenges for the bioinformatics analysis of phage genomes. In particular, ORF calling, genome annotation, noncoding RNA (ncRNA) identification, and the identification of transposons and insertions are all complicated in phage genome analysis. We provide a road map through the phage genome annotation pipeline, and discuss the challenges and solutions for phage genome annotation as we have implemented in the rapid annotation using subsystems (RAST) pipeline.

Optimizing and evaluating the reconstruction of Metagenome-assembled microbial genomes

Background

Microbiome/host interactions describe characteristics that affect the host's health. Shotgun metagenomics includes sequencing a random subset of the microbiome to analyze its taxonomic and metabolic potential. Reconstruction of DNA fragments into genomes from metagenomes (called metagenome-assembled genomes) assigns unknown fragments to taxa/function and facilitates discovery of novel organisms. Genome reconstruction incorporates sequence assembly and sorting of assembled sequences into bins, characteristic of a genome. However, the microbial community composition, including taxonomic and phylogenetic diversity may influence genome reconstruction. We determine the optimal reconstruction method for four microbiome projects that had variable sequencing platforms (IonTorrent and Illumina), diversity (high or low), and environment (coral reefs and kelp forests), using a set of parameters to select for optimal assembly and binning tools.

Methods

We tested the effects of the assembly and binning processes on population genome reconstruction using 105 marine metagenomes from 4 projects. Reconstructed genomes were obtained from each project using 3 assemblers (IDBA, MetaVelvet, and SPAdes) and 2 binning tools (GroopM and MetaBat). We assessed the efficiency of assemblers using statistics that including contig continuity and contig chimerism and the effectiveness of binning tools using genome completeness and taxonomic identification.

Results

We concluded that SPAdes, assembled more contigs (143,718 ± 124 contigs) of longer length (N50 = 1632 ± 108 bp), and incorporated the most sequences (sequences-assembled = 19.65%). The microbial richness and evenness were maintained across the assembly, suggesting low contig chimeras. SPAdes assembly was responsive to the biological and technological variations within the project, compared with other assemblers. Among binning tools, we conclude that MetaBat produced bins with less variation in GC content (average standard deviation: 1.49), low species richness (4.91 ± 0.66), and higher genome completeness (40.92 ± 1.75) across all projects. MetaBat extracted 115 bins from the 4 projects of which 66 bins were identified as reconstructed metagenome-assembled genomes with sequences belonging to a specific genus. We identified 13 novel genomes, some of which were 100% complete, but show low similarity to genomes within databases.

Conclusions

In conclusion, we present a set of biologically relevant parameters for evaluation to select for optimal assembly and binning tools. For the tools we tested, SPAdes assembler and MetaBat binning tools reconstructed quality metagenome-assembled genomes for the four projects. We also conclude that metagenomes from microbial communities that have high coverage of phylogenetically distinct, and low taxonomic diversity results in highest quality metagenome-assembled genomes.

Electronic supplementary material

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

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.

Assessing Species Diversity Using Metavirome Data: Methods and Challenges

Assessing biodiversity is an important step in the study of microbial ecology associated with a given environment. Multiple indices have been used to quantify species diversity, which is a key biodiversity measure. Measuring species diversity of viruses in different environments remains a challenge relative to measuring the diversity of other microbial communities. Metagenomics has played an important role in elucidating viral diversity by conducting metavirome studies; however, metavirome data are of high complexity requiring robust data preprocessing and analysis methods. In this review, existing bioinformatics methods for measuring species diversity using metavirome data are categorised broadly as either sequence similarity-dependent methods or sequence similarity-independent methods. The former includes a comparison of DNA fragments or assemblies generated in the experiment against reference databases for quantifying species diversity, whereas estimates from the latter are independent of the knowledge of existing sequence data. Current methods and tools are discussed in detail, including their applications and limitations. Drawbacks of the state-of-the-art method are demonstrated through results from a simulation. In addition, alternative approaches are proposed to overcome the challenges in estimating species diversity measures using metavirome data.

The skin microbiome of the common thresher shark (Alopias vulpinus) has low taxonomic and gene function β-diversity

The health of sharks, like all organisms, is linked to their microbiome. At the skin interface, sharks have dermal denticles that protrude above the mucus, which may affect the types of microbes that occur here. We characterized the microbiome from the skin of the common thresher shark (Alopias vulpinus) to investigate the structure and composition of the skin microbiome. On average 618 812 (80.9% ± S.D. 0.44%) reads per metagenomic library contained open reading frames; of those, between 7.6% and 12.8% matched known protein sequences. Genera distinguishing the A. vulpinus microbiome from the water column included, Pseudoalteromonas (12.8% ± 4.7 of sequences), Erythrobacter (5. 3% ± 0.5) and Idiomarina (4.2% ± 1.2) and distinguishing gene pathways included, cobalt, zinc and cadmium resistance (2.2% ± 0.1); iron acquisition (1.2% ± 0.1) and ton/tol transport (1.3% ± 0.08). Taxonomic community overlap (100 - dissimilarity index) was greater in the skin microbiome (77.6), relative to the water column microbiome (70.6) and a reference host-associated microbiome (algae: 71.5). We conclude the A. vulpinus skin microbiome is influenced by filtering processes, including biochemical and biophysical components of the shark skin and result in a structured microbiome.

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.

Putative archaeal viruses from the mesopelagic ocean

Oceanic viruses that infect bacteria, or phages, are known to modulate host diversity, metabolisms, and biogeochemical cycling, while the viruses that infect marine Archaea remain understudied despite the critical ecosystem roles played by their hosts. Here we introduce "MArVD", for Metagenomic Archaeal Virus Detector, an annotation tool designed to identify putative archaeal virus contigs in metagenomic datasets. MArVD is made publicly available through the online iVirus analytical platform. Benchmarking analysis of MArVD showed it to be >99% accurate and 100% sensitive in identifying the 127 known archaeal viruses among the 12,499 viruses in the VirSorter curated dataset. Application of MArVD to 10 viral metagenomes from two depth profiles in the Eastern Tropical North Pacific (ETNP) oxygen minimum zone revealed 43 new putative archaeal virus genomes and large genome fragments ranging in size from 10 to 31 kb. Network-based classifications, which were consistent with marker gene phylogenies where available, suggested that these putative archaeal virus contigs represented six novel candidate genera. Ecological analyses, via fragment recruitment and ordination, revealed that the diversity and relative abundances of these putative archaeal viruses were correlated with oxygen concentration and temperature along two OMZ-spanning depth profiles, presumably due to structuring of the host Archaea community. Peak viral diversity and abundances were found in surface waters, where Thermoplasmata 16S rRNA genes are prevalent, suggesting these archaea as hosts in the surface habitats. Together these findings provide a baseline for identifying archaeal viruses in sequence datasets, and an initial picture of the ecology of such viruses in non-extreme environments.

The Baltic Sea Virome: Diversity and Transcriptional Activity of DNA and RNA Viruses

Inferred virus-host relationships, community structures of ubiquitous ecologically relevant groups, and identification of transcriptionally active populations have been achieved with our Baltic Sea study. Further, these data, highlighting the transcriptional activity of viruses, represent one of the more powerful uses of omics concerning ecosystem health. The use of omics-related data to assess ecosystem health holds great promise for rapid and relatively inexpensive determination of perturbations and risk, explicitly with regard to viral assemblages, as no single marker gene is suitable for widespread taxonomic coverage.

Metagenomic and metatranscriptomic data were generated from size-fractionated samples from 11 sites within the Baltic Sea and adjacent marine waters of Kattegat and freshwater Lake Torneträsk in order to investigate the diversity, distribution, and transcriptional activity of virioplankton. Such a transect, spanning a salinity gradient from freshwater to the open sea, facilitated a broad genome-enabled investigation of natural as well as impacted aspects of Baltic Sea viral communities. Taxonomic signatures representative of phages within the widely distributed order Caudovirales were identified with enrichments in lesser-known families such as Podoviridae and Siphoviridae. The distribution of phage reported to infect diverse and ubiquitous heterotrophic bacteria (SAR11 clades) and cyanobacteria (Synechococcus sp.) displayed population-level shifts in diversity. Samples from higher-salinity conditions (>14 practical salinity units [PSU]) had increased abundances of viruses for picoeukaryotes, i.e., Ostreococcus. These data, combined with host diversity estimates, suggest viral modulation of diversity on the whole-community scale, as well as in specific prokaryotic and eukaryotic lineages. RNA libraries revealed single-stranded DNA (ssDNA) and RNA viral populations throughout the Baltic Sea, with ssDNA phage highly represented in Lake Torneträsk. Further, our data suggest relatively high transcriptional activity of fish viruses within diverse families known to have broad host ranges, such as Nodoviridae (RNA), Iridoviridae (DNA), and predicted zoonotic viruses that can cause ecological and economic damage as well as impact human health.

IMPORTANCE Inferred virus-host relationships, community structures of ubiquitous ecologically relevant groups, and identification of transcriptionally active populations have been achieved with our Baltic Sea study. Further, these data, highlighting the transcriptional activity of viruses, represent one of the more powerful uses of omics concerning ecosystem health. The use of omics-related data to assess ecosystem health holds great promise for rapid and relatively inexpensive determination of perturbations and risk, explicitly with regard to viral assemblages, as no single marker gene is suitable for widespread taxonomic coverage.

Environmental Viral Metagenomics Analyses in Aquaculture: Applications in Epidemiology and Disease Control

Studies on the epidemiology of viral diseases in aquaculture have for a long time depended on isolation of viruses from infected aquatic organisms. The role of aquatic environments in the epidemiology of viral diseases in aquaculture has not been extensively expounded mainly because of the lack of appropriate tools for environmental studies on aquatic viruses. However, the upcoming of metagenomics analyses opens great avenues in which environmental samples can be used to study the epidemiology of viral diseases outside their host species. Hence, in this review I have shown that epidemiological factors that influence the composition of viruses in different aquatic environments include ecological factors, anthropogenic activities and stocking densities of cultured organisms based on environmental metagenomics studies carried out this far. Ballast water transportation and global trade of aquatic organisms are the most common virus dispersal process identified this far. In terms of disease control for outdoor aquaculture systems, baseline data on viruses found in different environments intended for aquaculture use can be obtained to enable the design of effective disease control strategies. And as such, high-risk areas having a high specter of pathogenic viruses can be identified as an early warning system. As for the control of viral diseases for indoor recirculation aquaculture systems (RAS), the most effective disinfection methods able to eliminate pathogenic viruses from water used in RAS can be identified. Overall, the synopsis I have put forth in this review shows that environmental samples can be used to study the epidemiology of viral diseases in aquaculture using viral metagenomics analysis as an overture for the design of rational disease control strategies.

Aquatic viral metagenomics: Lights and shadows

Viruses are the most abundant biological entities on Earth, exceeding bacteria in most of the ecosystems. Specially in oceans, viruses are thought to be the major planktonic predators shaping microorganism communities and controlling ocean biological capacity. Plankton lysis by viruses plays an important role in ocean nutrient and energy cycles. Viral metagenomics has emerged as a powerful tool to uncover viral diversity in aquatic ecosystems through the use of Next Generation Sequencing. However, many of the commonly used viral sample preparation steps have several important biases that must be considered to avoid a misinterpretation of the results. In addition to biases caused by the purification of virus particles, viral DNA/RNA amplification and the preparation of genomic libraries could also introduce biases, and a detailed knowledge about such protocols is required. In this review, the main steps in the viral metagenomic workflow are described paying special attention to the potential biases introduced by each one.

Computational approaches to predict bacteriophage-host relationships

Metagenomics has changed the face of virus discovery by enabling the accurate identification of viral genome sequences without requiring isolation of the viruses. As a result, metagenomic virus discovery leaves the first and most fundamental question about any novel virus unanswered: What host does the virus infect? The diversity of the global virosphere and the volumes of data obtained in metagenomic sequencing projects demand computational tools for virus–host prediction. We focus on bacteriophages (phages, viruses that infect bacteria), the most abundant and diverse group of viruses found in environmental metagenomes. By analyzing 820 phages with annotated hosts, we review and assess the predictive power of in silico phage–host signals. Sequence homology approaches are the most effective at identifying known phage–host pairs. Compositional and abundance-based methods contain significant signal for phage–host classification, providing opportunities for analyzing the unknowns in viral metagenomes. Together, these computational approaches further our knowledge of the interactions between phages and their hosts. Importantly, we find that all reviewed signals significantly link phages to their hosts, illustrating how current knowledge and insights about the interaction mechanisms and ecology of coevolving phages and bacteria can be exploited to predict phage–host relationships, with potential relevance for medical and industrial applications.

New viruses infecting bacteria are increasingly being discovered in many environments through sequence-based explorations. To understand their role in microbial ecosystems, computational tools are indispensable to prioritize and guide experimental efforts. This review assesses and discusses a range of bioinformatic approaches to predict bacteriophage–host relationships when all that is known is their genome sequence.

Biogeography of Viruses in the Sea

Viral ecology is a rapidly progressing area of research, as molecular methods have improved significantly for targeted research on specific populations and whole communities. To interpret and synthesize global viral diversity and distribution, it is feasible to assess whether macroecology concepts can apply to marine viruses. We review how viral and host life history and physical properties can influence viral distribution in light of biogeography and metacommunity ecology paradigms. We highlight analytical approaches that can be applied to emerging global data sets and meta-analyses to identify individual taxa with global influence and drivers of emergent properties that influence microbial community structure by drawing on examples across the spectrum of viral taxa, from RNA to ssDNA and dsDNA viruses.

Combining genomic sequencing methods to explore viral diversity and reveal potential virus-host interactions

Viral diversity and virus-host interactions in oxygen-starved regions of the ocean, also known as oxygen minimum zones (OMZs), remain relatively unexplored. Microbial community metabolism in OMZs alters nutrient and energy flow through marine food webs, resulting in biological nitrogen loss and greenhouse gas production. Thus, viruses infecting OMZ microbes have the potential to modulate community metabolism with resulting feedback on ecosystem function. Here, we describe viral communities inhabiting oxic surface (10 m) and oxygen-starved basin (200 m) waters of Saanich Inlet, a seasonally anoxic fjord on the coast of Vancouver Island, British Columbia using viral metagenomics and complete viral fosmid sequencing on samples collected between April 2007 and April 2010. Of 6459 open reading frames (ORFs) predicted across all 34 viral fosmids, 77.6% (n = 5010) had no homology to reference viral genomes. These fosmids recruited a higher proportion of viral metagenomic sequences from Saanich Inlet than from nearby northeastern subarctic Pacific Ocean (Line P) waters, indicating differences in the viral communities between coastal and open ocean locations. While functional annotations of fosmid ORFs were limited, recruitment to NCBI's non-redundant “nr” database and publicly available single-cell genomes identified putative viruses infecting marine thaumarchaeal and SUP05 proteobacteria to provide potential host linkages with relevance to coupled biogeochemical cycling processes in OMZ waters. Taken together, these results highlight the power of coupled analyses of multiple sequence data types, such as viral metagenomic and fosmid sequence data with prokaryotic single cell genomes, to chart viral diversity, elucidate genomic and ecological contexts for previously unclassifiable viral sequences, and identify novel host interactions in natural and engineered ecosystems.

Ecology and evolution of viruses infecting uncultivated SUP05 bacteria as revealed by single-cell- and meta-genomics

Viruses modulate microbial communities and alter ecosystem functions. However, due to cultivation bottlenecks, specific virus-host interaction dynamics remain cryptic. In this study, we examined 127 single-cell amplified genomes (SAGs) from uncultivated SUP05 bacteria isolated from a model marine oxygen minimum zone (OMZ) to identify 69 viral contigs representing five new genera within dsDNA Caudovirales and ssDNA Microviridae. Infection frequencies suggest that ∼1/3 of SUP05 bacteria is viral-infected, with higher infection frequency where oxygen-deficiency was most severe. Observed Microviridae clonality suggests recovery of bloom-terminating viruses, while systematic co-infection between dsDNA and ssDNA viruses posits previously unrecognized cooperation modes. Analyses of 186 microbial and viral metagenomes revealed that SUP05 viruses persisted for years, but remained endemic to the OMZ. Finally, identification of virus-encoded dissimilatory sulfite reductase suggests SUP05 viruses reprogram their host's energy metabolism. Together, these results demonstrate closely coupled SUP05 virus-host co-evolutionary dynamics with the potential to modulate biogeochemical cycling in climate-critical and expanding OMZs.

Depth-stratified functional and taxonomic niche specialization in the 'core' and 'flexible' Pacific Ocean Virome

Microbes drive myriad ecosystem processes, and their viruses modulate microbial-driven processes through mortality, horizontal gene transfer, and metabolic reprogramming by viral-encoded auxiliary metabolic genes (AMGs). However, our knowledge of viral roles in the oceans is primarily limited to surface waters. Here we assess the depth distribution of protein clusters (PCs) in the first large-scale quantitative viral metagenomic data set that spans much of the pelagic depth continuum (the Pacific Ocean Virome; POV). This established 'core' (180 PCs; one-third new to science) and 'flexible' (423K PCs) community gene sets, including niche-defining genes in the latter (385 and 170 PCs are exclusive and core to the photic and aphotic zones, respectively). Taxonomic annotation suggested that tailed phages are ubiquitous, but not abundant (<5% of PCs) and revealed depth-related taxonomic patterns. Functional annotation, coupled with extensive analyses to document non-viral DNA contamination, uncovered 32 new AMGs (9 core, 20 photic and 3 aphotic) that introduce ways in which viruses manipulate infected host metabolism, and parallel depth-stratified host adaptations (for example, photic zone genes for iron-sulphur cluster modulation for phage production, and aphotic zone genes for high-pressure deep-sea survival). Finally, significant vertical flux of photic zone viruses to the deep sea was detected, which is critical for interpreting depth-related patterns in nature. Beyond the ecological advances outlined here, this catalog of viral core, flexible and niche-defining genes provides a resource for future investigation into the organization, function and evolution of microbial molecular networks to mechanistically understand and model viral roles in the biosphere.

A highly abundant bacteriophage discovered in the unknown sequences of human faecal metagenomes

Metagenomics, or sequencing of the genetic material from a complete microbial community, is a promising tool to discover novel microbes and viruses. Viral metagenomes typically contain many unknown sequences. Here we describe the discovery of a previously unidentified bacteriophage present in the majority of published human faecal metagenomes, which we refer to as crAssphage. Its ~97 kbp genome is six times more abundant in publicly available metagenomes than all other known phages together; it comprises up to 90% and 22% of all reads in virus-like particle (VLP)-derived metagenomes and total community metagenomes, respectively; and it totals 1.68% of all human faecal metagenomic sequencing reads in the public databases. The majority of crAssphage-encoded proteins match no known sequences in the database, which is why it was not detected before. Using a new co-occurrence profiling approach, we predict a Bacteroides host for this phage, consistent with Bacteroides-related protein homologues and a unique carbohydrate-binding domain encoded in the phage genome.

Metagenomic studies of microbial communities often report DNA sequences from unidentified viruses. Here, Dutilh et al. analyse metagenomic data to reveal the complete genome of an abundant, ubiquitous virus from human faeces, and predict that the virus infects bacteria of the Bacteroides group.

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

Long-standing questions in marine viral ecology are centered on understanding how viral assemblages change along gradients in space and time. However, investigating these fundamental ecological questions has been challenging due to incomplete representation of naturally occurring viral diversity in single gene- or morphology-based studies and an inability to identify up to 90% of reads in viral metagenomes (viromes). Although protein clustering techniques provide a significant advance by helping organize this unknown metagenomic sequence space, they typically use only ∼75% of the data and rely on assembly methods not yet tuned for naturally occurring sequence variation. Here, we introduce an annotation- and assembly-free strategy for comparative metagenomics that combines shared k-mer and social network analyses (regression modeling). This robust statistical framework enables visualization of complex sample networks and determination of ecological factors driving community structure. Application to 32 viromes from the Pacific Ocean Virome dataset identified clusters of samples broadly delineated by photic zone and revealed that geographic region, depth, and proximity to shore were significant predictors of community structure. Within subsets of this dataset, depth, season, and oxygen concentration were significant drivers of viral community structure at a single open ocean station, whereas variability along onshore-offshore transects was driven by oxygen concentration in an area with an oxygen minimum zone and not depth or proximity to shore, as might be expected. Together these results demonstrate that this highly scalable approach using complete metagenomic network-based comparisons can both test and generate hypotheses for ecological investigation of viral and microbial communities in nature.

Caught in the middle with multiple displacement amplification: the myth of pooling for avoiding multiple displacement amplification bias in a metagenome

BACKGROUND

Shotgun metagenomics has become an important tool for investigating the ecology of microorganisms. Underlying these investigations is the assumption that metagenome sequence data accurately estimates the census of microbial populations. Multiple displacement amplification (MDA) of microbial community DNA is often used in cases where it is difficult to obtain enough DNA for sequencing; however, MDA can result in amplification biases that may impact subsequent estimates of population census from metagenome data. Some have posited that pooling replicate MDA reactions negates these biases and restores the accuracy of population analyses. This assumption has not been empirically tested.

RESULTS

Using mock viral communities, we examined the influence of pooling on population-scale analyses. In pooled and single reaction MDA treatments, sequence coverage of viral populations was highly variable and coverage patterns across viral genomes were nearly identical, indicating that initial priming biases were reproducible and that pooling did not alleviate biases. In contrast, control unamplified sequence libraries showed relatively even coverage across phage genomes.

CONCLUSIONS

MDA should be avoided for metagenomic investigations that require quantitative estimates of microbial taxa and gene functional groups. While MDA is an indispensable technique in applications such as single-cell genomics, amplification biases cannot be overcome by combining replicate MDA reactions. Alternative library preparation techniques should be utilized for quantitative microbial ecology studies utilizing metagenomic sequencing approaches.

Metagenomic analysis of size-fractionated picoplankton in a marine oxygen minimum zone

Marine oxygen minimum zones (OMZs) support diverse microbial communities with roles in major elemental cycles. It is unclear how the taxonomic composition and metabolism of OMZ microorganisms vary between particle-associated and free-living size fractions. We used amplicon (16S rRNA gene) and shotgun metagenome sequencing to compare microbial communities from large (>1.6 μm) and small (0.2-1.6 μm) filter size fractions along a depth gradient in the OMZ off Chile. Despite steep vertical redox gradients, size fraction was a significantly stronger predictor of community composition compared to depth. Phylogenetic diversity showed contrasting patterns, decreasing towards the anoxic OMZ core in the small size fraction, but exhibiting maximal values at these depths within the larger size fraction. Fraction-specific distributions were evident for key OMZ taxa, including anammox planctomycetes, whose coding sequences were enriched up to threefold in the 0.2-1.6 μm community. Functional gene composition also differed between fractions, with the >1.6 μm community significantly enriched in genes mediating social interactions, including motility, adhesion, cell-to-cell transfer, antibiotic resistance and mobile element activity. Prokaryotic transposase genes were three to six fold more abundant in this fraction, comprising up to 2% of protein-coding sequences, suggesting that particle surfaces may act as hotbeds for transposition-based genome changes in marine microbes. Genes for nitric and nitrous oxide reduction were also more abundant (three to seven fold) in the larger size fraction, suggesting microniche partitioning of key denitrification steps. These results highlight an important role for surface attachment in shaping community metabolic potential and genome content in OMZ microorganisms.

Comparative metagenomics: natural populations of induced prophages demonstrate highly unique, lower diversity viral sequences

To understand the similarities and differences between a free living viral population and its co-occurring temperate population, metagenomes of each type were prepared from the same seawater sample from Tampa Bay, FL. Libraries were prepared from extracted DNA of the ambient viruses and induced prophages from the co-occurring, viral-reduced microbial assemblage. Duplicate libraries were also prepared using the same DNA amplified by multiple displacement amplification. A non-viral-reduced, induced, amplified viral dataset from the same site in 2005 was reanalysed for temporal comparison. The induced viral metagenome was higher in identifiable virus sequences and differed from the other three datasets based on principal component, rarefaction, trinucleotide composition and contig spectrum analyses. This study indicated that induced prophages are unique and have lower overall community diversity than ambient viral populations from the same site. Both of the amplified contemporary metagenomes were enriched in single-stranded DNA (ssDNA) viral sequences. Six and 16 complete, circular ssDNA viral genomes were assembled from the amplified induced and ambient libraries, respectively, mostly similar to circoviruses. The amplified ambient metagenome contained genomes similar to an RNA-DNA hybrid virus recently identified in a hot spring and to an ssDNA virus infecting the diatom Chaetoceros.

Genome of a SAR116 bacteriophage shows the prevalence of this phage type in the oceans

The abundance, genetic diversity, and crucial ecological and evolutionary roles of marine phages have prompted a large number of metagenomic studies. However, obtaining a thorough understanding of marine phages has been hampered by the low number of phage isolates infecting major bacterial groups other than cyanophages and pelagiphages. Therefore, there is an urgent requirement for the isolation of phages that infect abundant marine bacterial groups. In this study, we isolated and characterized HMO-2011, a phage infecting a bacterium of the SAR116 clade, one of the most abundant marine bacterial lineages. HMO-2011, which infects "Candidatus Puniceispirillum marinum" strain IMCC1322, has an ~55-kb dsDNA genome that harbors many genes with novel features rarely found in cultured organisms, including genes encoding a DNA polymerase with a partial DnaJ central domain and an atypical methanesulfonate monooxygenase. Furthermore, homologs of nearly all HMO-2011 genes were predominantly found in marine metagenomes rather than cultured organisms, suggesting the novelty of HMO-2011 and the prevalence of this phage type in the oceans. A significant number of the viral metagenome sequences obtained from the ocean surface were best assigned to the HMO-2011 genome. The number of reads assigned to HMO-2011 accounted for 10.3%-25.3% of the total reads assigned to viruses in seven viromes from the Pacific and Indian Oceans, making the HMO-2011 genome the most or second-most frequently assigned viral genome. Given its ability to infect the abundant SAR116 clade and its widespread distribution, Puniceispirillum phage HMO-2011 could be an important resource for marine virus research.