Seasonal Dynamics and Metagenomic Characterization of Marine Viruses in Goseong Bay, Korea

Microbial interactions strengthen deterministic processes during community assembly in a subtropical estuary

Systematic studies on the assembly process and driving mechanisms of microbial communities in estuaries with diverse seasonal and spatial scales are still limited. In this study, high-throughput sequencing, and microbial network analysis were combined to decipher the impact of environmental changes and biological interactions on the maintenance of microbial diversity patterns in the Jiulong River Estuary (JRE). The results showed that overall, stochastic processes dominated the bacterioplankton community assembly in the estuary, accounting for 49.66-74.78 % of the total. Additionally, bacterioplankton community diversity varied significantly across seasons and subzones. Specifically, the concentration of soluble reactive phosphorus (SRP) in the estuary steadily reduced from winter to summer, and the corresponding bacterioplankton community interactions gradually shifted from the weakest interaction in winter to the strongest in summer. The deterministic processes contributed more than half (50.34 %) to microbial assembly in the summer, but only 25.22 % in winter. Deterministic processes prevailed in the seaward with low SRP concentrations and strong bacterioplankton community interactions, while stochastic processes contributed 70.14 % to the assembly of microbial communities riverward. Biotic and abiotic factors, such as nutrients and microbial interactions, jointly drove the seasonal and spatial patterns of bacterioplankton community assembly, but overall, nutrients played a dominant role. Nevertheless, the contributions of nutrients and microbial interactions were equivalent in spatial assembly processes, albeit nutrients were the primary seasonal driver of the bacterioplankton community assembly process. This study emphasizes the significance of microbial interactions in the bacterioplankton community assemblage. These findings provide new and comprehensive insights into the microbial communities' organization in estuaries.

K. Habeebullah S, Ali AK, Al-Zakri WM. Diversity and spatiotemporal variations in bacterial and archaeal communities within Kuwaiti territorial waters of the Northwest Arabian Gulf

Kuwaiti territorial waters of the northwest Arabian Gulf represent a unique aquatic ecosystem prone to various environmental and anthropogenic stressors that pose significant constraints on the resident biota which must withstand extreme temperatures, salinity levels, and reducing conditions, among other factors to survive. Such conditions create the ideal environment for investigations into novel functional genetic adaptations of resident organisms. Firstly, however, it is essential to identify said organisms and understand the dynamic nature of their existence. Thus, this study provides the first comprehensive analysis of bacterial and archaeal community structures in the unique waters of Kuwait located in the Northwest Arabian Gulf and analyzes their variations with respect to depth, season, and location, as well as their susceptibility to changes in abundance with respect to various physicochemical parameters. Importantly, this study is the first of its kind to utilize a shotgun metagenomics approach with sequencing performed at an average depth of 15 million paired end reads per sample, which allows for species-level community profiling and sets the framework for future functional genomic investigations. Results showed an approximately even abundance of both archaeal (42.9%) and bacterial (57.1%) communities, but significantly greater diversity among the bacterial population, which predominantly consisted of members of the Proteobacteria, Cyanobacteria, and Bacteroidetes phyla in decreasing order of abundance. Little to no significant variations as assessed by various metrics including alpha and beta diversity analyses were observed in the abundance of archaeal and bacterial populations with respect to depth down the water column. Furthermore, although variations in differential abundance of key genera were detected at each of the three sampling locations, measurements of species richness and evenness revealed negligible variation (ANOVA p<0.05) and only a moderately defined community structure (ANOSIM r2 = 0.243; p>0.001) between the various locations. Interestingly, abundance of archaeal community members showed a significant increase (log2 median ratio of RA = 2.6) while the bacterial population showed a significant decrease (log2 median ratio = -1.29) in the winter season. These findings were supported by alpha and beta diversity analyses as well (ANOSIM r2 = 0.253; p>0.01). Overall, this study provides the first in-depth analysis of both bacterial and archaeal community structures developed using a shotgun metagenomic approach in the waters of the Northwest Arabian Gulf thus providing a framework for future investigations of functional genetic adaptations developed by resident biota attempting to survive in the uniquely extreme conditions to which they are exposed.

Development and validation of a real-time PCR assay protocol for the specific detection and quantification of pelagiphages in seawater samples

Earth is inhabited by numerous adaptations of cellular forms shaped by the persistent scrutiny of natural selection. Thus, as natural selection has fixed beneficial adaptations of functional traits, cellular life has conquered almost all environmental niches on our planet. However, cellular life succumbs in number and genetic diversity to viruses. Among all viruses, phages are highly prevalent in diverse environments, and due to their vast genetic diversity and abundance, their relevant role as significant players in several ecological processes is now fully recognized. Pelagiphages, bacteriophages infecting bacteria of the SAR11 clade, are the most abundant viruses in the oceans. However, the ecological contribution of pelagiphages on populations of Pelagibacterales remains largely underestimated. An essential aspect of estimating the impact of bacteriophages is their absolute and precise quantification, which provides relevant information about the host-virus interactions and the structure of viral assemblages. Consequently, due to its abundance and claimed influence in the biogeochemical cycling of elements, the accurate quantification of pelagiphages results in an essential task. This study describes the development and validation of a sensitive, specific, accurate and reproducible qPCR platform targeting pelagiphages. Moreover, this method allowed the detection and quantification of pelagiphages in the Gulf of California for the first time.

Prevalence of Viral Frequency-Dependent Infection in Coastal Marine Prokaryotes Revealed Using Monthly Time Series Virome Analysis

Viruses infecting marine prokaryotes have a large impact on the diversity and dynamics of their hosts. Model systems suggest that viral infection is frequency dependent and constrained by the virus-host encounter rate. However, it is unclear whether frequency-dependent infection is pervasive among the abundant prokaryotic populations with different temporal dynamics. To address this question, we performed a comparison of prokaryotic and viral communities using 16S rRNA amplicon and virome sequencing based on samples collected monthly for 2 years at a Japanese coastal site, Osaka Bay. Concurrent seasonal shifts observed in prokaryotic and viral community dynamics indicated that the abundance of viruses correlated with that of their predicted host phyla (or classes). Cooccurrence network analysis between abundant prokaryotes and viruses revealed 6,423 cooccurring pairs, suggesting a tight coupling of host and viral abundances and their "one-to-many" correspondence. Although stable dominant species, such as SAR11, showed few cooccurring viruses, a fast succession of their viruses suggests that viruses infecting these populations changed continuously. Our results suggest that frequency-dependent viral infection prevails in coastal marine prokaryotes regardless of host taxa and temporal dynamics. IMPORTANCE There is little room for doubt that viral infection is prevalent among abundant marine prokaryotes regardless of their taxa or growth strategy. However, comprehensive evaluations of viral infections in natural prokaryotic communities are still technically difficult. In this study, we examined viral infection in abundant prokaryotes by monitoring the monthly dynamics of prokaryotic and viral communities at a eutrophic coastal site, Osaka Bay. We compared the community dynamics of viruses with those of their putative hosts based on genome-based in silico host prediction. We observed frequent cooccurrence among the predicted virus-host pairs, suggesting that viral infection is prevalent in abundant prokaryotes regardless of their taxa or temporal dynamics. This likely indicates that frequent lysis of the abundant prokaryotes via viral infection has a considerable contribution to the biogeochemical cycling and maintenance of prokaryotic community diversity.

Phage strategies facilitate bacterial coexistence under environmental variability

Bacterial communities are often exposed to temporal variations in resource availability, which exceed bacterial generation times and thereby affect bacterial coexistence. Bacterial population dynamics are also shaped by bacteriophages, which are a main cause of bacterial mortality. Several strategies are proposed in the literature to describe infections by phages, such as "Killing the Winner", "Piggyback the loser" (PtL) or "Piggyback the Winner" (PtW). The two temperate phage strategies PtL and PtW are defined by a change from lytic to lysogenic infection when the host density changes, from high to low or from low to high, respectively. To date, the occurrence of different phage strategies and their response to environmental variability is poorly understood. In our study, we developed a microbial trophic network model using ordinary differential equations (ODEs) and performed 'in silico' experiments. To model the switch from the lysogenic to the lytic cycle, we modified the lysis rate of infected bacteria and their growth was turned on or off using a density-dependent switching point. We addressed whether and how the different phage strategies facilitate bacteria coexistence competing for limiting resources. We also studied the impact of a fluctuating resource inflow to evaluate the response of the different phage strategies to environmental variability. Our results show that the viral shunt (i.e. nutrient release after bacterial lysis) leads to an enrichment of the system. This enrichment enables bacterial coexistence at lower resource concentrations. We were able to show that an established, purely lytic model leads to stable bacterial coexistence despite fluctuating resources. Both temperate phage models differ in their coexistence patterns. The model of PtW yields stable bacterial coexistence at a limited range of resource supply and is most sensitive to resource fluctuations. Interestingly, the purely lytic phage strategy and PtW both result in stable bacteria coexistence at oligotrophic conditions. The PtL model facilitates stable bacterial coexistence over a large range of stable and fluctuating resource inflow. An increase in bacterial growth rate results in a higher resilience to resource variability for the PtL and the lytic infection model. We propose that both temperate phage strategies represent different mechanisms of phages coping with environmental variability. Our study demonstrates how phage strategies can maintain bacterial coexistence in constant and fluctuating environments.

Prophage Genomics and Ecology in the Family Rhodobacteraceae

Roseobacters are globally abundant bacteria with critical roles in carbon and sulfur biogeochemical cycling. Here, we identified 173 new putative prophages in 79 genomes of Rhodobacteraceae. These prophages represented 1.3 ± 0.15% of the bacterial genomes and had no to low homology with reference and metagenome-assembled viral genomes from aquatic and terrestrial ecosystems. Among the newly identified putative prophages, 35% encoded auxiliary metabolic genes (AMGs), mostly involved in secondary metabolism, amino acid metabolism, and cofactor and vitamin production. The analysis of integration sites and gene homology showed that 22 of the putative prophages were actually gene transfer agents (GTAs) similar to a GTA of Rhodobacter capsulatus. Twenty-three percent of the predicted prophages were observed in the TARA Oceans viromes generated from free viral particles, suggesting that they represent active prophages capable of induction. The distribution of these prophages was significantly associated with latitude and temperature. The prophages most abundant at high latitudes encoded acpP, an auxiliary metabolic gene involved in lipid synthesis and membrane fluidity at low temperatures. Our results show that prophages and gene transfer agents are significant sources of genomic diversity in roseobacter, with potential roles in the ecology of this globally distributed bacterial group.

The communities and functional profiles of virioplankton along a salinity gradient in a subtropical estuary

Viruses are the major drivers shaping microorganismal communities, and impact marine biogeochemical cycling. They are affected by various environmental parameters, such as salinity. Although the spatiotemporal distribution and dynamics of virioplankton have been extensively studied in saline environments, few detailed studies of community structure and function of viruses along salinity gradients have been conducted. Here, we used the 16S and 18S rRNA gene amplicon and metagenomic sequencing from a subtropical estuary (Pearl River Estuary, PRE; located in Shenzhen, Guangdong Province, China) to explore how viral community composition and function vary along a salinity gradient. Results showed that the detected viruses were mainly bacteriophages. The double-stranded DNA viruses were the most abundant (especially Siphoviridae, Myoviridae, Mimiviridae, Phycodnaviridae, and Podoviridae), followed by a small number of single-stranded DNA (Circoviridae) and RNA (Retroviridae) viruses. Viral biodiversity significantly declined and community structure varied greatly along the salinity gradient. The salinity, ammonium and dissolved oxygen were dominated factors influencing the community composition of viruses. Association network analysis showed that viruses had a negative effect on multiple host taxa (prokaryotic and eukaryotic species). Metagenomic data revealed that the main viral functional potential was involved in organic matter metabolism by carbohydrate-active enzymes (CAZymes). Deeper comparative functional analyses showed that viruses in the low-salinity environment had more carbohydrate-binding module and glycosidase hydrolases activities than those under high-salinity conditions. However, an opposite pattern was observed for carbohydrate esterases. These results suggest that virus-encoded CAZyme genes may alter the bacterial metabolism in estuaries. Overall, our results demonstrate that there is a spatial heterogeneity in the composition and function of virioplankton along a salinity gradient. This study enhances our understanding of viral distribution and their contribution to regulating carbon degradation throughout environments with varying salinities in subtropical estuaries.

The Unexplored Virome of Two Atlantic Coast Fish: Contribution of Next-Generation Sequencing to Fish Virology

Much of the knowledge on viruses is focused on those that can be propagated using cell-cultures or that can cause disease in humans or in economically important animals and plants. However, this only reflects a small portion of the virosphere. Therefore, in this study, we explore by targeted next-generation sequencing, how the virome varies between Atlantic horse mackerels and gilthead seabreams from fisheries and aquaculture from the center and south regions of Portugal. Viral genomes potentially pathogenic to fish and crustaceans, as well as to humans, were identified namelyese included Astroviridae, Nodaviridae, Hepadnaviridae, Birnaviridae, Caliciviridae, and Picornaviridae families. Also bacteriophages sequences were identified corresponding to the majority of sequencese detected, with Myoviridae, Podoviridae, and Siphoviridae, the most widespread families in both fish species. However, these findings can also be due to the presence of bacteria in fish tissues, or even to contamination. Overall, seabreams harbored viruses from a smaller number of families in comparison with mackerels. Therefore, the obtained data show that fish sold for consumption can harbor a high diversity of viruses, many of which are unknown, reflecting the overall uncharacterized virome of fish. While cross-species transmission of bonafide fish viruses to humans is unlikely, the finding of human pathogenic viruses in fish suggest that fish virome can be a potential threat regarding food safety.

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.

Metagenomics revealing new virus species in farm and pet animals and aquaculture

Viral metagenomics is slowly taking over the traditional and widely used molecular techniques for the investigation of pathogenic viruses responsible for illness and inflicting great economic burden on the farm animal industry. Owing to the continued improvements in sequencing technologies and the dramatic reduction of per base costs of sequencing the use of next generation sequencing have been key factors in this progress. Discoveries linked to viral metagenomics are expected to be beneficial to the field of veterinary medicine starting from the development of better diagnostic assays to the design of new subunit vaccines with minimal investments. With these achievements the research has taken a giant leap even toward the better healthcare of animals and, as a result, the animal sector could be growing at an unprecedented pace.

Metagenomics Approaches in Discovery and Development of New Bioactive Compounds from Marine Actinomycetes

Marine actinomycetes are prolific sources of marine drug discovery system contributing for several bioactive compounds of biomedical prominence. Metagenomics, a culture-independent technique through its sequence- and function-based screening has led to the discovery and synthesis of numerous biologically significant compounds like polyketide synthase, Non-ribosomal peptide synthetase, antibiotics, and biocatalyst. While metagenomics offers different advantages over conventional sequencing techniques, they also have certain limitations including bias classification, non-availability of quality DNA samples, heterologous expression, and host selection. The assimilation of advanced amplification and screening methods such as φ29 DNA polymerase, Next-Generation Sequencing, Cosmids, and recent bioinformatics tools like automated genome mining, anti-SMASH have shown promising results to overcome these constrains. Consequently, functional genomics and bioinformatics along with synthetic biology will be crucial for the success of the metagenomic approach and indeed for exploring new possibilities among the microbial consortia for the future drug discovery process.

Seasonal dynamics of bacterial communities in the surface seawater around subtropical Xiamen Island, China, as determined by 16S rRNA gene profiling

This study investigated the microbial structure in the surface seawater from five coastal sites around Xiamen Island, China, over four seasons to evaluate seasonal environmental fluctuations impact on them. This subtropical island is characterized by long, hot, humid summers, and short, mild, dry winters. All sites were dominated by Proteobacteria, Bacteroidetes, Cyanobacteria, and Firmicutes; microbial community composition was similar across four seasons. However, larger proportions of Gammaproteobacteria and Bacillus were observed during the summer than during any other season. The high ratio of Bacillus, Bacteroidetes, and Clostridia richness to Alphaproteobacteria richness in the summer, suggested that the sites we tested were heavily affected by waste water to other seasons. Correlation-based network analyses among the bacterial species and environmental variables indicated important connections between physiochemical variables and specific taxonomic groups. Collectively, our results suggested that seasonal shifts and wastewater pollution together shape the structures of the microbial communities around Xiamen Island.

Antibiotic resistance genes in bacteriophages from diverse marine habitats

Although antibiotic resistance represents a significant and growing threat to human and environmental health worldwide, the contribution of bacteriophages (phages) to the acquisition and spread of antibiotic resistance genes (ARGs) in the environment has not been extensively explored. In this study, a comprehensive analysis of several viromes from diverse marine habitats was performed to investigate whether or not phages carry ARGs. The analysis provides strong evidence that phages from marine habitats are potential reservoirs of ARGs. In fact, genes conferring resistance to aminocoumarin, bacitracin and multidrug resistance (particularly the mexB gene) were found in all analyzed marine viromes. Given this, the role of phages as reservoirs of ARGs should not be underestimated considering their global distribution.

Cobaviruses - a new globally distributed phage group infecting Rhodobacteraceae in marine ecosystems

Bacteriophages are widely considered to influence bacterial communities, however most phages are still unknown or not studied well enough to understand their ecological roles. We have isolated two phages infecting Lentibacter sp. SH36, affiliated with the marine Roseobacter group, and retrieved similar phage genomes from publicly available metagenomics databases. Phylogenetic analysis placed the new phages within the Cobavirus group, in the here newly proposed genus Siovirus and subfamily Riovirinae of the Podoviridae. Gene composition and presence of direct terminal repeats in cultivated cobaviruses point toward a genome replication and packaging strategy similar to the T7 phage. Investigation of the genomes suggests that viral lysis of the cell proceeds via the canonical holin-endolysin pathway. Cobaviral hosts include members of the genera Lentibacter, Sulfitobacter and Celeribacter of the Roseobacter group within the family Rhodobacteraceae (Alphaproteobacteria). Screening more than 5,000 marine metagenomes, we found cobaviruses worldwide from temperate to tropical waters, in the euphotic zone, mainly in bays and estuaries, but also in the open ocean. The presence of cobaviruses in protist metagenomes as well as the phylogenetic neighborhood of cobaviruses in glutaredoxin and ribonucleotide reductase trees suggest that cobaviruses could infect bacteria associated with phototrophic or grazing protists. With this study, we expand the understanding of the phylogeny, classification, genomic organization, biogeography and ecology of this phage group infecting marine Rhodobacteraceae.

The Human Gut Virome in Hypertension

Objectives: Previous studies have reported that the gut microbiome has an important link with the development of hypertension. Though previous researches have focused on the links of gut bacteria with hypertension, little has been known about the linkage of gut viruses to hypertension and the development of hypertension, largely due to the lack of data mining tools for such investigation. In this work, we have analyzed 196 fecal metagenomic data related to hypertension aiming to profile the gut virome and link the gut virome to pre-hypertension and hypertension. Design: Here, we have applied a statistically sound method for mining of gut virome data and linking gut virome to hypertension. We characterized the viral composition and bacterial composition of 196 samples, identified the viral-type of each sample and linked gut virome to hypertension. Results: We stratified these 196 fecal samples into two viral-types and selected 32 viruses as the biomarkers for these groups. We found that viruses could have a superior resolution and discrimination power than bacteria for differentiation of healthy samples and pre-hypertension samples, as well as hypertension samples. Moreover, as to the co-occurrence networks linking viruses and bacteria, we found increasingly pervasive virus-bacteria linkages from healthy people to pre-hypertension people to hypertension patients. Conclusion: Overall, our results have shown ample indications of the link between human gut virome and hypertension, and could help provide microbial solutions toward early diagnoses of hypertension.

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.

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.