Evaluation of Sequencing Library Preparation Protocols for Viral Metagenomic Analysis from Pristine Aquifer Groundwaters

Exploring viral contamination in urban groundwater and runoff

The reliance of the global population on urban aquifers is steadily increasing, and urban aquifers are susceptible to pathogenic contamination through sources such as sewer leakage or urban runoff. However, there is insufficient monitoring of groundwater quality in urban areas. In this study, quantitative polymerase chain reaction (qPCR) was employed to evaluate the presence of human fecal viral indicators and viral pathogens in urban wastewater (n = 13) and groundwater (n = 12) samples from four locations in Barcelona with different degrees of urbanization, as well as in runoff samples (n = 2). Additionally, a target enrichment sequencing (TES) approach was utilized to explore the viral diversity within groundwater and runoff samples, offering insights into viral contamination and potential virus transmission routes in urban areas. Human adenovirus (HAdV) was identified in all wastewater samples, 67 % (8/12) of groundwater samples, and one runoff sample by qPCR indicating human viral fecal contamination. The viral pathogen Norovirus genogroup GI (NoV GI) was detected in wastewater and two winter groundwater samples from highly and medium urbanized areas. NoV genogroup GII (NoV GII), Enterovirus (EV) and SARS-CoV-2 were exclusively detected in wastewater. Human and other vertebrate viruses were detected in groundwater and runoff samples using TES. This study gives insights about the virome present in urban water sources, emphasizing the need for thorough monitoring and deeper understanding to address emerging public health concerns.

Exploring viral diversity and metagenomics in livestock: insights into disease emergence and spillover risks in cattle

Cattle have a significant impact on human societies in terms of both economics and health. Viral infections pose a relevant problem as they directly or indirectly disrupt the balance within cattle populations. This has negative consequences at the economic level for producers and territories, and also jeopardizes human health through the transmission of zoonotic diseases that can escalate into outbreaks or pandemics. To establish prevention strategies and control measures at various levels (animal, farm, region, or global), it is crucial to identify the viral agents present in animals. Various techniques, including virus isolation, serological tests, and molecular techniques like PCR, are typically employed for this purpose. However, these techniques have two major drawbacks: they are ineffective for non-culturable viruses, and they only detect a small fraction of the viruses present. In contrast, metagenomics offers a promising approach by providing a comprehensive and unbiased analysis for detecting all viruses in a given sample. It has the potential to identify rare or novel infectious agents promptly and establish a baseline of healthy animals. Nevertheless, the routine application of viral metagenomics for epidemiological surveillance and diagnostics faces challenges related to socioeconomic variables, such as resource availability and space dedicated to metagenomics, as well as the lack of standardized protocols and resulting heterogeneity in presenting results. This review aims to provide an overview of the current knowledge and prospects for using viral metagenomics to detect and identify viruses in cattle raised for livestock, while discussing the epidemiological and clinical implications.

MuDoGeR: Multi-Domain Genome recovery from metagenomes made easy

Several computational frameworks and workflows that recover genomes from prokaryotes, eukaryotes and viruses from metagenomes exist. Yet, it is difficult for scientists with little bioinformatics experience to evaluate quality, annotate genes, dereplicate, assign taxonomy and calculate relative abundance and coverage of genomes belonging to different domains. MuDoGeR is a user-friendly tool tailored for those familiar with Unix command-line environment that makes it easy to recover genomes of prokaryotes, eukaryotes and viruses from metagenomes, either alone or in combination. We tested MuDoGeR using 24 individual-isolated genomes and 574 metagenomes, demonstrating the applicability for a few samples and high throughput. While MuDoGeR can recover eukaryotic viral sequences, its characterization is predominantly skewed towards bacterial and archaeal viruses, reflecting the field's current state. However, acting as a dynamic wrapper, the MuDoGeR is designed to constantly incorporate updates and integrate new tools, ensuring its ongoing relevance in the rapidly evolving field. MuDoGeR is open-source software available at https://github.com/mdsufz/MuDoGeR. Additionally, MuDoGeR is also available as a Singularity container.

High niche specificity and host genetic diversity of groundwater viruses

Viruses are key members of microbial communities that exert control over host abundance and metabolism, thereby influencing ecosystem processes and biogeochemical cycles. Aquifers are known to host taxonomically diverse microbial life, yet little is known about viruses infecting groundwater microbial communities. Here, we analysed 16 metagenomes from a broad range of groundwater physicochemistries. We recovered 1571 viral genomes that clustered into 468 high-quality viral operational taxonomic units. At least 15% were observed to be transcriptionally active, although lysis was likely constrained by the resource-limited groundwater environment. Most were unclassified (95%), and the remaining 5% were Caudoviricetes. Comparisons with viruses inhabiting other aquifers revealed no shared species, indicating substantial unexplored viral diversity. In silico predictions linked 22.4% of the viruses to microbial host populations, including to ultra-small prokaryotes, such as Patescibacteria and Nanoarchaeota. Many predicted hosts were associated with the biogeochemical cycling of carbon, nitrogen, and sulfur. Metabolic predictions revealed the presence of 205 putative auxiliary metabolic genes, involved in diverse processes associated with the utilization of the host's intracellular resources for biosynthesis and transformation reactions, including those involved in nucleotide sugar, glycan, cofactor, and vitamin metabolism. Viruses, prokaryotes overall, and predicted prokaryotic hosts exhibited narrow spatial distributions, and relative abundance correlations with the same groundwater parameters (e.g. dissolved oxygen, nitrate, and iron), consistent with host control over viral distributions. Results provide insights into underexplored groundwater viruses, and indicate the large extent to which viruses may manipulate microbial communities and biogeochemistry in the terrestrial subsurface.

Evaluation of multiple displacement amplification for metagenomic analysis of low biomass samples

Combining multiple displacement amplification (MDA) with metagenomics enables the analysis of samples with extremely low DNA concentrations, making them suitable for high-throughput sequencing. Although amplification bias and nonspecific amplification have been reported from MDA-amplified samples, the impact of MDA on metagenomic datasets is not well understood. We compared three MDA methods (i.e. bulk MDA, emulsion MDA, and primase MDA) for metagenomic analysis of two DNA template concentrations (approx. 1 and 100 pg) derived from a microbial community standard "mock community" and two low biomass environmental samples (i.e. borehole fluid and groundwater). We assessed the impact of MDA on metagenome-based community composition, assembly quality, functional profiles, and binning. We found amplification bias against high GC content genomes but relatively low nonspecific amplification such as chimeras, artifacts, or contamination for all MDA methods. We observed MDA-associated representational bias for microbial community profiles, especially for low-input DNA and with the primase MDA method. Nevertheless, similar taxa were represented in MDA-amplified libraries to those of unamplified samples. The MDA libraries were highly fragmented, but similar functional profiles to the unamplified libraries were obtained for bulk MDA and emulsion MDA at higher DNA input and across these MDA libraries for the groundwater sample. Medium to low-quality bins were possible for the high input bulk MDA metagenomes for the most simple microbial communities, borehole fluid, and mock community. Although MDA-based amplification should be avoided, it can still reveal meaningful taxonomic and functional information from samples with extremely low DNA concentration where direct metagenomics is otherwise impossible.

The smallest in the deepest: the enigmatic role of viruses in the deep biosphere

It is commonly recognized that viruses control the composition, metabolism, and evolutionary trajectories of prokaryotic communities, with resulting vital feedback on ecosystem functioning and nutrient cycling in a wide range of ecosystems. Although the deep biosphere has been estimated to be the largest reservoir for viruses and their prokaryotic hosts, the biology and ecology of viruses therein remain poorly understood. The deep virosphere is an enigmatic field of study in which many critical questions are still to be answered. Is the deep virosphere simply a repository for deeply preserved, non-functioning virus particles? Or are deep viruses infectious agents that can readily infect suitable hosts and subsequently shape microbial populations and nutrient cycling? Can the cellular content released by viral lysis, and even the organic structures of virions themselves, serve as the source of bioavailable nutrients for microbial activity in the deep biosphere as in other ecosystems? In this review, we synthesize our current knowledge of viruses in the deep biosphere and seek to identify topics with the potential for substantial discoveries in the future.

Combining Flow Cytometry and Metagenomics Improves Recovery of Metagenome-Assembled Genomes in a Cell Culture from Activated Sludge

The recovery of metagenome-assembled genomes is biased towards the most abundant species in a given community. To improve the identification of species, even if only dominant species are recovered, we investigated the integration of flow cytometry cell sorting with bioinformatics tools to recover metagenome-assembled genomes. We used a cell culture of a wastewater microbial community as our model system. Cells were separated based on fluorescence signals via flow cytometry cell sorting into sub-communities: dominant gates, low abundant gates, and outer gates into subsets of the original community. Metagenome sequencing was performed for all groups. The unsorted community was used as control. We recovered a total of 24 metagenome-assembled genomes (MAGs) representing 11 species-level genome operational taxonomic units (gOTUs). In addition, 57 ribosomal operational taxonomic units (rOTUs) affiliated with 29 taxa at species level were reconstructed from metagenomic libraries. Our approach suggests a two-fold increase in the resolution when comparing sorted and unsorted communities. Our results also indicate that species abundance is one determinant of genome recovery from metagenomes as we can recover taxa in the sorted libraries that are not present in the unsorted community. In conclusion, a combination of cell sorting and metagenomics allows the recovery of MAGs undetected without cell sorting.

Novel Virus Identification through Metagenomics: A Systematic Review

Metagenomic Next Generation Sequencing (mNGS) allows the evaluation of complex microbial communities, avoiding isolation and cultivation of each microbial species, and does not require prior knowledge of the microbial sequences present in the sample. Applications of mNGS include virome characterization, new virus discovery and full-length viral genome reconstruction, either from virus preparations enriched in culture or directly from clinical and environmental specimens. Here, we systematically reviewed studies that describe novel virus identification through mNGS from samples of different origin (plant, animal and environment). Without imposing time limits to the search, 379 publications were identified that met the search parameters. Sample types, geographical origin, enrichment and nucleic acid extraction methods, sequencing platforms, bioinformatic analytical steps and identified viral families were described. The review highlights mNGS as a feasible method for novel virus discovery from samples of different origins, describes which kind of heterogeneous experimental and analytical protocols are currently used and provides useful information such as the different commercial kits used for the purification of nucleic acids and bioinformatics analytical pipelines.

Molecular characterization of avian metapneumovirus subtype C detected in wild mallards (Anas platyrhynchos) in The Netherlands

Avian metapneumovirus (AMPV) represents a long-term threat to the poultry industry due to its etiological role in the induction of acute respiratory disease and/or egg drop syndrome in domestic turkeys, chickens, and ducks. Although this disease is commonly referred to as turkey rhinotracheitis, the host range of AMPV encompasses many avian species. We have screened 1323 oropharyngeal- and cloacal swab samples obtained from wild mallards in the Netherlands from 2017 to 2019 by RT-PCR using a degenerate primer pair to detect all members of the Paramyxoviridae and Pneumoviridae or an avian metapneumovirus subtype C (AMPV-C)-specific RT-qPCR assay. We identified a total of seven cases of AMPV-C infections in wild, healthy mallards (Anas platyrhynchos), of which two AMPV-C positive samples were further processed using next-generation sequencing. Phylogenetic analysis of the two complete genomes showed that the newly identified AMPV-C strains share closest sequence identity (97%) with Eurasian lineage AMPV-C strains identified in Muscovy ducks in China that presented with severe respiratory disease and egg production loss in 2011. Further analysis of G protein amino acid sequences showed a high degree of variability between the newly identified AMPV-C variants. PONDR scoring of the G protein has revealed the ectodomain of AMPV-C to be partitioned into a long intrinsically disordered and short ordered region, giving insights into AMPV G protein structural biology. In summary, we provide the first report of full-length AMPV-C genome sequences derived from wild birds in Europe. This emphasizes the need for further surveillance efforts to better characterize the host range, epidemiologic distribution, and pathogenicity of AMPV-C to determine the risk posed by cross-species jumps from wildfowl to domesticated avian species.

Comprehensive Evaluation of RNA and DNA Viromic Methods Based on Species Richness and Abundance Analyses Using Marmot Rectal Samples

Viral metagenomics is the most powerful tool to profile viromic composition for a given sample. Different viromic methods, including amplification-free ones, have been developed, but choosing them for different purposes requires comprehensive benchmarks. Here, we assessed the performance of four routinely used methods, i.e., multiple displacement amplification (MDA), direct metagenomic sequencing (MTG), sequence-independent single-primer amplification (SIA), and metatranscriptomic sequencing (MTT), using marmot rectal samples as the templates spiked with five known viruses of different genome types. The obtained clean data were differently contaminated by host and bacterial genomes, resulting in MDA having the most, with ~72.1%, but MTT had only ~7.5% data, useful for follow-up viromic analysis. MDA showed a broader spectrum with higher efficiency to profile the DNA virome, and MTT captured almost all RNA viruses with extraordinary sensitivity; hence, they are advisable in richness-based viromic studies. MTG was weak in capturing single-stranded DNA viruses, and SIA could detect both RNA and DNA viruses but with high randomness. Due to biases to certain types of viruses, the four methods caused different alterations to species abundance compared to the initial virus composition. SIA and MDA introduced greater stochastic errors to relative abundances of species, genus, and family taxa, whereas the two amplification-free methods were more tolerant toward such errors and thus are recommendable in abundance-based analyses. In addition, genus taxon is a compromising analytic level that ensures technically supported and biologically and/or ecologically meaningful viromic conclusions. IMPORTANCE Viral metagenomics can be roughly divided into species richness-based studies and species abundance-based analyses. Viromic methods with different principles have been developed, but rational selection of these techniques according to different purposes requires comprehensive understanding of their properties. By assessing the four most widely used methods using template samples, we found that multiple displacement amplification (MDA) and metatranscriptomic sequencing (MTT) are advisable for species richness-based viromic studies, as they show excellent efficiency to detect DNA and RNA viruses. Meanwhile, metagenomic sequencing (MTG) and MTT are more compatible with stochastic errors of methods introduced into relative abundance of viromic taxa and hence are rational choices in species abundance-based analyses. This study also highlights that MTG needs to tackle host genome contamination and ameliorate the capacity to detect single-stranded DNA viruses in the future, and the MTT method requires an improvement in bacterial rRNA depletion prior to library preparation.

Assessment of Rapid MinION Nanopore DNA Virus Meta-Genomics Using Calves Experimentally Infected with Bovine Herpes Virus-1

Bovine respiratory disease (BRD), which is the leading cause of morbidity and mortality in cattle, is caused by numerous known and unknown viruses and is responsible for the widespread use of broad-spectrum antibiotics despite the use of polymicrobial BRD vaccines. Viral metagenomics sequencing on the portable, inexpensive Oxford Nanopore Technologies MinION sequencer and sequence analysis with its associated user-friendly point-and-click Epi2ME cloud-based pathogen identification software has the potential for point-of-care/same-day/sample-to-result metagenomic sequence diagnostics of known and unknown BRD pathogens to inform a rapid response and vaccine design. We assessed this potential using in vitro viral cell cultures and nasal swabs taken from calves that were experimentally challenged with a single known BRD-associated DNA virus, namely, bovine herpes virus 1. Extensive optimisation of the standard Oxford Nanopore library preparation protocols, particularly a reduction in the PCR bias of library amplification, was required before BoHV-1 could be identified as the main virus in the in vitro cell cultures and nasal swab samples within approximately 7 h from sample to result. In addition, we observed incorrect assignment of the bovine sequence to bacterial and viral taxa due to the presence of poor-quality bacterial and viral genome assemblies in the RefSeq database used by the EpiME Fastq WIMP pathogen identification software.

Women in the European Virus Bioinformatics Center

Viruses are the cause of a considerable burden to human, animal and plant health, while on the other hand playing an important role in regulating entire ecosystems. The power of new sequencing technologies combined with new tools for processing "Big Data" offers unprecedented opportunities to answer fundamental questions in virology. Virologists have an urgent need for virus-specific bioinformatics tools. These developments have led to the formation of the European Virus Bioinformatics Center, a network of experts in virology and bioinformatics who are joining forces to enable extensive exchange and collaboration between these research areas. The EVBC strives to provide talented researchers with a supportive environment free of gender bias, but the gender gap in science, especially in math-intensive fields such as computer science, persists. To bring more talented women into research and keep them there, we need to highlight role models to spark their interest, and we need to ensure that female scientists are not kept at lower levels but are given the opportunity to lead the field. Here we showcase the work of the EVBC and highlight the achievements of some outstanding women experts in virology and viral bioinformatics.

A snapshot of the global drinking water virome: Diversity and metabolic potential vary with residual disinfectant use

Viruses are important drivers of microbial community ecology and evolution, influencing microbial mortality, metabolism, and horizontal gene transfer. However, the effects of viruses remain largely unknown in many environments, including in drinking water systems. Drinking water metagenomic studies have offered a whole community perspective of bacterial impacts on water quality, but have not yet considered the influences of viruses. In this study, we address this gap by mining viral DNA sequences from publicly available drinking water metagenomes from distribution systems in six countries around the world. These datasets provide a snapshot of the taxonomic diversity and metabolic potential of the global drinking water virome; and provide an opportunity to investigate the effects of geography, climate, and drinking water treatment practices on viral diversity. Both environmental conditions and differences in sample processing were found to influence the viral composition. Using free chlorine as the residual disinfectant was associated with clear differences in viral taxonomic diversity and metabolic potential, with significantly fewer viral populations and less even viral community structures than observed in distribution systems without residual disinfectant. Additionally, drinking water viruses carry antibiotic resistance genes (ARGs), as well as genes to survive oxidative stress and nitrogen limitation. Through this study, we have demonstrated that viral communities are diverse across drinking water systems and vary with the use of residual disinfectant. Our findings offer directions for future research to develop a more robust understanding of how virus-bacteria interactions in drinking water distribution systems affect water quality.

Viruses and Their Interactions With Bacteria and Archaea of Hypersaline Great Salt Lake

Viruses play vital biogeochemical and ecological roles by (a) expressing auxiliary metabolic genes during infection, (b) enhancing the lateral transfer of host genes, and (c) inducing host mortality. Even in harsh and extreme environments, viruses are major players in carbon and nutrient recycling from organic matter. However, there is much that we do not yet understand about viruses and the processes mediated by them in the extreme environments such as hypersaline habitats. The Great Salt Lake (GSL) in Utah, United States is a hypersaline ecosystem where the biogeochemical role of viruses is poorly understood. This study elucidates the diversity of viruses and describes virus–host interactions in GSL sediments along a salinity gradient. The GSL sediment virosphere consisted of Haloviruses (32.07 ± 19.33%) and members of families Siphoviridae (39.12 ± 19.8%), Myoviridae (13.7 ± 6.6%), and Podoviridae (5.43 ± 0.64%). Our results demonstrate that salinity alongside the concentration of organic carbon and inorganic nutrients (nitrogen and phosphorus) governs the viral, bacteria, and archaeal diversity in this habitat. Computational host predictions for the GSL viruses revealed a wide host range with a dominance of viruses that infect Proteobacteria, Actinobacteria, and Firmicutes. Identification of auxiliary metabolic genes for photosynthesis (psbA), carbon fixation (rbcL, cbbL), formaldehyde assimilation (SHMT), and nitric oxide reduction (NorQ) shed light on the roles played by GSL viruses in biogeochemical cycles of global relevance.

A library preparation optimized for metagenomics of RNA viruses

Our understanding of the viral communities associated to animals has not yet reached the level attained on the bacteriome. This situation is due to, among others, technical challenges in adapting metagenomics using high-throughput sequencing to the study of RNA viromes in animals. Although important developments have been achieved in most steps of viral metagenomics, there is yet a key step that has received little attention: the library preparation. This situation differs from bacteriome studies in which developments in library preparation have largely contributed to the democratisation of metagenomics. Here, we present a library preparation optimized for metagenomics of RNA viruses from insect vectors of viral diseases. The library design allows a simple PCR-based preparation, such as those routinely used in bacterial metabarcoding, that is adapted to shotgun sequencing as required in viral metagenomics. We first optimized our library preparation using mock viral communities and then validated a full metagenomic approach incorporating our preparation in two pilot studies with field-caught insect vectors; one including a comparison with a published metagenomic protocol. Our approach provided a fold increase in virus-like sequences compared to other studies, and nearly-full genomes from new virus species. Moreover, our results suggested conserved trends in virome composition within a population of a mosquito species. Finally, the sensitivity of our approach was compared to a commercial diagnostic PCR for the detection of an arbovirus in field-caught insect vectors. Our approach could facilitate studies on viral communities from animals and the democratization of metagenomics in community ecology of viruses.

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.

Deep Impact of Random Amplification and Library Construction Methods on Viral Metagenomics Results

Clinical metagenomics is a broad-range agnostic detection method of pathogens, including novel microorganisms. A major limit is the low pathogen load compared to the high background of host nucleic acids. To overcome this issue, several solutions exist, such as applying a very high depth of sequencing, or performing a relative enrichment of viral genomes associated with capsids. At the end, the quantity of total nucleic acids is often below the concentrations recommended by the manufacturers of library kits, which necessitates to random amplify nucleic acids. Using a pool of 26 viruses representative of viral diversity, we observed a deep impact of the nature of sample (total nucleic acids versus RNA only), the reverse transcription, the random amplification and library construction method on virus recovery. We further optimized the two most promising methods and assessed their performance with fully characterized reference virus stocks. Good genome coverage and limit of detection lower than 100 or 1000 genome copies per mL of plasma, depending on the genome viral type, were obtained from a three million reads dataset. Our study reveals that optimized random amplification is a technique of choice when insufficient amounts of nucleic acid are available for direct libraries constructions.

Metagenomic to evaluate norovirus genomic diversity in oysters: Impact on hexamer selection and targeted capture-based enrichment

Human virus transmission through food consumption has been identified since many years and the international trade increases the risk of dissemination of viral pathogens. The development of metagenomic approach holds many promises for the surveillance of viruses in food and water. This work aimed to analyze norovirus diversity and to evaluate strain-dependent accumulation patterns in three oyster types by using a metagenomic approach. Different hexamer sets to prime cDNA were evaluated before capture-based approach to enhance virus reads recovery during deep sequencing. The study includes the use of technical replicates of artificially contaminated oysters and the analysis of multiple negatives controls. Results showed a clear impact of the hexamer set used for cDNA synthesis. A set of In-house designed (I-HD) hexamers, selected to lower mollusk amplification, gave promising results in terms of viral reads abundancy. However, the best correlation between C_T values, thus concentrations, and number of reads was observed using random hexamers. Random hexamers also provided the highest numbers of reads and allowed the identification of sequence of different human enteric viruses. Regarding human norovirus, different genogroups and genotypes were identified among contigs longer than 500 bp. Two full genomes and six sequences longer than 3600 bases were obtained allowing a precise strain identification. The use of technical triplicates was found valuable to increase the chances to sequence viral strains present at low concentrations. Analyzing viral contamination in shellfish samples is quite challenging, however this work demonstrates that the recovery of full genome or long contigs, allowing clear identification of viral strains is possible.

Prokaryotic and Viral Community Composition of Freshwater Springs in Florida, USA

Aquifers, which are essential underground freshwater reservoirs worldwide, are understudied ecosystems that harbor diverse forms of microbial life. This study investigated the abundance and composition of prokaryotic and viral communities in the outflow of five springs across northern Florida, USA, as a proxy of microbial communities found in one of the most productive aquifers in the world, the Floridan aquifer. The average abundances of virus-like particles and prokaryotic cells were slightly lower than those reported from other groundwater systems, ranging from 9.6 × 103 ml-1 to 1.1 × 105 ml-1 and 2.2 × 103 ml-1 to 3.4 × 104 ml-1, respectively. Despite all of the springs being fed by the Floridan aquifer, sequencing of 16S rRNA genes and viral metagenomes (viromes) revealed unique communities in each spring, suggesting that groundwater microbial communities are influenced by land usage in recharge zones. The prokaryotic communities were dominated by Bacteria, and though the most abundant phyla (Proteobacteria, Cyanobacteria, and Bacteroidetes) were found in relatively high abundance across springs, variation was seen at finer taxonomic resolution. The viral sequences were most similar to those described from other aquatic environments. Sequencing resulted in the completion of 58 novel viral genomes representing members of the order Caudovirales as well as prokaryotic and eukaryotic single-stranded DNA (ssDNA) viruses. Sequences similar to those of ssDNA viruses were detected at all spring sites and dominated the identifiable sequences at one spring site, showing that these small viruses merit further investigation in groundwater systems.IMPORTANCE Aquifer systems may hold up to 40% of the total microbial biomass on Earth. However, little is known about the composition of microbial communities within these critical freshwater ecosystems. Here, we took advantage of Florida's first-magnitude springs (the highest spring classification based on water discharge), each discharging at least 246 million liters of water each day from the Floridan aquifer system (FAS), to investigate prokaryotic and viral communities from the aquifer. The FAS serves as a major source of potable water in the Southeastern United States, providing water for large cities and citizens in three states. Unfortunately, the health of the FAS and its associated springs has declined in the past few decades due to nutrient loading, increased urbanization and agricultural activity in aquifer recharge zones, and saltwater intrusion. This is the first study to describe the prokaryotic and viral communities in Florida's first-magnitude springs, providing a baseline against which to compare future ecosystem change.

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.