Analysis of viral and bacterial communities in groundwater associated with contaminated land

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.

Longitudinal change of gut microbiota in hypertensive disorders in pregnancy: a nested case-control and Mendelian randomization study

Mounting evidence has shown that gut microbiota (GM) is related to hypertensive disorders in pregnancy (HDP), however, most studies only focused on one time point in pregnancy. In this study, we conducted a nested case-control study utilizing a follow-up cohort, resulting in the collection of 47 HDP patients and 30 healthy controls. The GM profiles were explored using 16S rRNA sequencing at three time points during pregnancy. The diversity analysis of GM showed no significant difference between HDP patients and controls, however, we found 21 differential GM during pregnancy. Trend analysis showed that there are statistical differences in the relative abundance of Thermomonas, Xanthomonas, and Phenylobacteriumat during pregnancy in the gestational hypertension group, and of Xanthomonas, Polycyclovorans, and Phenylobacterium in the control group. The correlation study found that six genera of GM are related to blood pressure. Furthermore, the MR analysis identified the causal relationship between Methanobrevibacter and pre-eclampsia (PE). This study first explored the longitudinal change of GM in HDP patients during pregnancy, found the differential GM, and detected the causal association. Our findings may promote the prevention and treatment of HDP from the perspective of GM and provide valuable insights into the pathogenesis of HDP.

A Temperate Sinorhizobium Phage, AP-16-3, Closely Related to Phage 16-3: Mosaic Genome and Prophage Analysis

Soil Sinorhizobium phage AP-16-3, a strain phylogenetically close to Rhizobium phage 16-3, was isolated in a mountainous region of Dagestan, belonging to the origin of cultivated plants in the Caucasus, according to Vavilov N.I. The genome of phage AP-16-3 is 61 kbp in size and contains 62 ORFs, of which 42 ORFs have homologues in the genome of Rhizobium phage 16-3, which was studied in the 1960s-1980s. A search for Rhizobium phage 16-3-related sequences was performed in the genomes of modern strains of root nodule bacteria belonging to different species, genera, and families. A total of 43 prophages of interest were identified out of 437 prophages found in the genomes of 42 strains, of which 31 belonged to Sinorhizobium meliloti species. However, almost all of the mentioned prophages contained single ORFs, and only two prophages contained 51 and 39 ORFs homologous to phages related to 16-3. These prophages were detected in S. meliloti NV1.1.1 and Rh. leguminosarum OyaliB strains belonging to different genera; however, the similarity level of these two prophages did not exceed 14.7%. Analysis of the orphan genes in these prophages showed that they encoded predominantly virion structural elements, but also enzymes and an extensive group of hypothetical proteins belonging to the L, S, and E regions of viral genes of phage 16-3. The data obtained indicate that temperate phages related to 16-3 had high infectivity against nodule bacteria and participated in intragenomic recombination events involving other phages, and in horizontal gene transfer between rhizobia of different genera. According to the data obtained, it is assumed that the repetitive lysogenic cycle of temperate bacteriophages promotes the dissolution of the phage genetic material in the host bacterial genome, and radical updating of phage and host bacterial genomes takes place.

Sample and library preparation approaches for the analysis of the virome of irrigation water

BACKGROUND

The virome (i.e. community of mainly RNA and DNA eukaryotic viruses and bacteriophages) of waters is yet to be extensively explored. In particular, the virome of waters used for irrigation could therefore potentially carry viral pathogens that can contaminate fresh produce. One problem in obtaining viral sequences from irrigation waters is the relatively low amount of virus particles, as well as the presence of human, bacterial and protozoan cells. The present aimed study was to compare different processing, amplification, and sequencing approaches for virome characterization in irrigation waters.

RESULTS

Our analyses considered percentages of viral reads, values for diversity indices and number of families found in sequencing results. The results obtained suggest that enrichment protocols using two (bezonase and microccocal nuclease) or four enzymes at once (bezonase, microccocal nuclease, DNAse and RNase), regardless of an Amicon filtration step, are more appropriate than separated enzymatic treatments for virome characterization in irrigation water. The NetoVIR protocol combined with the ScriptSeq v2 RNA-Seq Library (P0-L20 protocol) showed the highest percentages of RNA viruses and identified the higher number of families.

CONCLUSION

Although virome characterization applied in irrigation waters is an important tool for protecting public health by informing on circulating human and zoonotic infections, optimized and standardized procedures should be followed to reduce the variability of results related to either the sample itself and the downstream bioinformatics analyses. Our results show that virome characterization can be an important tool in the discovery of pathogenic viruses in the environment and can be used to inform and optimize reference-based detection methods provided that appropriate and rigorous controls are included. © 2023 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Virus impacted community adaptation in oligotrophic groundwater environment revealed by Hi-C coupled metagenomic and viromic study

Viruses play a crucial role in microbial mortality, diversity and biogeochemical cycles. Groundwater is the largest global freshwater and one of the most oligotrophic aquatic systems on Earth, but how microbial and viral communities are shaped in this special habitat is largely unexplored. In this study, we collected groundwater samples from 23 to 60 m aquifers at Yinchuan Plain, China. In total, 1920 non-reductant viral contigs were retrieved from metagenomes and viromes constructed by Illumina and Nanopore hybrid sequencing. Only 3% of them could be clustered with known viruses, most of which were Caudoviricetes. Coupling 1.2 Tb Hi-C sequencing with CRISPR matching and homology search, we connected 469 viruses with their hosts while some viral clusters presented a broad-host-range trait. Meanwhile, a large proportion of biosynthesis related auxiliary metabolism genes were identified. Those characteristics might benefit viruses for a better survival in this special oligotrophic environment. Additionally, the groundwater virome showed genomic features distinct from those of the open ocean and wastewater treatment facilities in GC distribution and unannotated gene compositions. This paper expands the current knowledge of the global viromic records and serves as a foundation for a more thorough understanding of viruses in groundwater.

The Deterioration of Agronomical Traits of the Continuous Cropping of Stevia Is Associated With the Dynamics of Soil Bacterial Community

Stevia rebaudiana Bertoni is grown worldwide as an important, natural sweetener resource plant. The yield of steviol glycosides (SVglys) is greatly influenced by continuous cropping. In this study, we collected the roots, rhizosphere soils, and bulk soils from 2 years of continuous cropping (Y2) and 8 years of continuous cropping (Y8). A high-throughput sequencing technology based on Illumina Hiseq 2500 platform was used to study the structure and diversity of bacterial communities in the roots and soils of stevia with different years of continuous cropping. The results demonstrated that although the content of a group of SVglys was significantly increased in stevia of long-term continuous cropping, it inhibited the growth of plants and lowered the leaf dry weight; as a result, the total amount of SVglys was significantly decreased. Meanwhile, continuous cropping changed the physicochemical properties and the bacterial composition communities of soil. The different sampling sources of the root, rhizosphere soil, and bulk soil had no impact on the richness of bacterial communities, while it exhibited obvious effects on the diversity of bacterial communities. Continuous cropping had a stronger effect on the bacterial community composition in rhizosphere soil than in root and bulk soil. Based on linear discriminant analysis effect size (LEfSe), in the rhizosphere soil of Y8, the relative abundance of some beneficial bacterial genera of Sphingomonas, Devosia, Streptomyces, and Flavobacterium decreased significantly, while the relative abundance of Polycyclovorans, Haliangium, and Nitrospira greatly increased. Moreover, the soil pH and nutrient content, especially the soil organic matter, were correlated with the relative abundance of predominant bacteria at the genus level. This study provides a theoretical basis for uncovering the mechanism of obstacles in continuous stevia cropping and provides guidance for the sustainable development of stevia.

Groundwater contaminated with short-chain chlorinated paraffins and microbial responses

The vertical migrations of toxic and persistent short-chain chlorinated paraffins (SCCPs) in soils as well as the microbial responses have been reported, however, there is a paucity of data on the resulting groundwater contamination. Here, we determined the concentration and congener profile of SCCPs in the groundwater beneath a production plant of chlorinated paraffins (CPs) and characterized the microbial community to explore their responses to SCCPs. Results showed that SCCPs ranged from not detected to 70.3 μg/L, with C₁₃-CPs (11.2-65.8%) and Cl₇-CPs (27.2-50.6%), in mass ratio, as the dominant groups. Similar to the distribution pattern in soils, SCCPs in groundwater were distributed in hotspot pattern. CP synthesis was the source of SCCPs in groundwater and the entire contamination plume significantly migrated downgradient, while there was an apparent hysteresis of C₁₃-CP migration. Groundwater microbial community was likely shaped by both hydrogeological condition (pH and depth) and SCCPs. Specifically, the microbial community responded to the contamination by forming a co-occurrence network with "small world" feature, where Desulfobacca, Desulfomonile, Ferritrophicum, Methylomonas, Syntrophobacter, Syntrophorhabdus, Syntrophus, and Thermoanaerobaculum were the keystone taxa. Furthermore, the interrelations between bacterial taxa and SCCPs indicated that the microbial community might cooperate to achieve the dechlorination and mineralization of SCCPs through either anaerobic organohalide respiration mainly functioned by the keystone taxa, or cometabolic degradation processes functioned by Aquabacterium and Hydrogenophaga. Results of this study would provide a better understanding of the environmental behavior and ecological effects of SCCPs in groundwater systems.

Ecogenomics of Groundwater Phages Suggests Niche Differentiation Linked to Specific Environmental Tolerance

Viruses are ubiquitous microbiome components, shaping ecosystems via strain-specific predation, horizontal gene transfer and redistribution of nutrients through host lysis. Viral impacts are important in groundwater ecosystems, where microbes drive many nutrient fluxes and metabolic processes; however, little is known about the diversity of viruses in these environments. We analyzed four groundwater plasmidomes (the entire plasmid content of an environment) and identified 200 viral sequences, which clustered into 41 genus-level viral clusters (approximately equivalent to viral genera) including 9 known and 32 putative new genera. We used publicly available bacterial whole-genome sequences (WGS) and WGS from 261 bacterial isolates from this groundwater environment to identify potential viral hosts. We linked 76 of the 200 viral sequences to a range of bacterial phyla, the majority associated with Proteobacteria, followed by Firmicutes, Bacteroidetes, and Actinobacteria. The publicly available WGS enabled mapping bacterial hosts to several viral sequences. The WGS of groundwater isolates increased the depth of host prediction by allowing host identification at the strain level. The latter included 4 viruses that were almost entirely (>99% query coverage, >99% identity) identified as integrated in the genomes of Pseudomonas, Acidovorax, and Castellaniella strains, resulting in high-confidence host assignments. Lastly, 21 of these viruses carried putative auxiliary metabolite genes for metal and antibiotic resistance, which might drive their infection cycles and/or provide selective advantage to infected hosts. Exploring the groundwater virome provides a necessary foundation for integration of viruses into ecosystem models where they are key players in microbial adaption to environmental stress. IMPORTANCE To our knowledge, this is the first study to identify the bacteriophage distribution in a groundwater ecosystem shedding light on their prevalence and distribution across metal-contaminated and background sites. Our study is uniquely based on selective sequencing of solely the extrachromosomal elements of a microbiome followed by analysis for viral signatures, thus establishing a more focused approach for phage identifications. Using this method, we detected several novel phage genera along with those previously established. Our approach of using the whole-genome sequences of hundreds of bacterial isolates from the same site enabled us to make host assignments with high confidence, several at strain levels. Certain phage genes suggest that they provide an environment-specific selective advantage to their bacterial hosts. Our study lays the foundation for future research on directed phage isolations using specific bacterial host strains to further characterize groundwater phages, their life cycles, and their effects on groundwater microbiome and biogeochemistry.

Ecophysiological Features Shape the Distribution of Prophages and CRISPR in Sulfate Reducing Prokaryotes

Sulfate reducing prokaryotes (SRP) are a phylogenetically and physiologically diverse group of microorganisms that use sulfate as an electron acceptor. SRP have long been recognized as key players of the carbon and sulfur cycles, and more recently, they have been identified to play a relevant role as part of syntrophic and symbiotic relations and the human microbiome. Despite their environmental relevance, there is a poor understanding about the prevalence of prophages and CRISPR arrays and how their distribution and dynamic affect the ecological role of SRP. We addressed this question by analyzing the results of a comprehensive survey of prophages and CRISPR in a total of 91 genomes of SRP with several genotypic, phenotypic, and physiological traits, including genome size, cell volume, minimum doubling time, cell wall, and habitat, among others. Our analysis discovered 81 prophages in 51 strains, representing the 56% of the total evaluated strains. Prophages are non-uniformly distributed across the SRP phylogeny, where prophage-rich lineages belonged to Desulfovibrionaceae and Peptococcaceae. Furthermore, our study found 160 CRISPR arrays in 71 SRP, which is more abundant and widely spread than previously expected. Although there is no correlation between presence and abundance of prophages and CRISPR arrays at the strain level, our analysis showed that there is a directly proportional relation between cellular volumes and number of prophages per cell. This result suggests that there is an additional selective pressure for strains with smaller cells to get rid of foreign DNA, such as prophages, but not CRISPR, due to less availability of cellular resources. Analysis of the prophage genes encoding viral structural proteins reported that 44% of SRP prophages are classified as Myoviridae, and comparative analysis showed high level of homology, but not synteny, among prophages belonging to the Family Desulfovibrionaceae. We further recovered viral-like particles and structures that resemble outer membrane vesicles from D. vulgaris str. Hildenborough. The results of this study improved the current understanding of dynamic interactions between prophages and CRISPR with their hosts in both cultured and hitherto-uncultured SRP strains, and how their distribution affects the microbial community dynamics in several sulfidogenic natural and engineered environments.

Seasonal dynamics in taxonomy and function within bacterial and viral metagenomic assemblages recovered from a freshwater agricultural pond

BACKGROUND

Ponds are important freshwater habitats that support both human and environmental activities. However, relative to their larger counterparts (e.g. rivers, lakes), ponds are understudied, especially with regard to their microbial communities. Our study aimed to fill this knowledge gap by using culture-independent, high-throughput sequencing to assess the dynamics, taxonomy, and functionality of bacterial and viral communities in a freshwater agricultural pond.

RESULTS

Water samples (n = 14) were collected from a Mid-Atlantic agricultural pond between June 2017 and May 2018 and filtered sequentially through 1 and 0.2 μm filter membranes. Total DNA was then extracted from each filter, pooled, and subjected to 16S rRNA gene and shotgun sequencing on the Illumina HiSeq 2500 platform. Additionally, on eight occasions water filtrates were processed for viral metagenomes (viromes) using chemical concentration and then shotgun sequenced. A ubiquitous freshwater phylum, Proteobacteria was abundant at all sampling dates throughout the year. However, environmental characteristics appeared to drive the structure of the community. For instance, the abundance of Cyanobacteria (e.g. Nostoc) increased with rising water temperatures, while a storm event appeared to trigger an increase in overall bacterial diversity, as well as the relative abundance of Bacteroidetes. This event was also associated with an increase in the number of antibiotic resistance genes. The viral fractions were dominated by dsDNA of the order Caudovirales, namely Siphoviridae and Myovirdae.

CONCLUSIONS

Overall, this study provides one of the largest datasets on pond water microbial ecology to date, revealing seasonal trends in the microbial taxonomic composition and functional potential.

Metatranscriptomic Analysis of Oil-Exposed Seawater Bacterial Communities Archived by an Environmental Sample Processor (ESP)

The use of natural marine bacteria as "oil sensors" for the detection of pollution events can be suggested as a novel way of monitoring oil occurrence at sea. Nucleic acid-based devices generically called genosensors are emerging as potentially promising tools for in situ detection of specific microbial marker genes suited for that purpose. Functional marker genes are particularly interesting as targets for oil-related genosensing but their identification remains a challenge. Here, seawater samples, collected in tanks with oil addition mimicking a realistic oil spill scenario, were filtered and archived by the Environmental Sample Processor (ESP), a fully robotized genosensor, and the samples were then used for post-retrieval metatranscriptomic analysis. After extraction, RNA from ESP-archived samples at start, Day 4 and Day 7 of the experiment was used for sequencing. Metatranscriptomics revealed that several KEGG pathways were significantly enriched in samples exposed to oil. However, these pathways were highly expressed also in the non-oil-exposed water samples, most likely as a result of the release of natural organic matter from decaying phytoplankton. Temporary peaks of aliphatic alcohol and aldehyde dehydrogenases and monoaromatic ring-degrading enzymes (e.g., ben, box, and dmp clusters) were observed on Day 4 in both control and oil-exposed and non-exposed tanks. Few alkane 1-monooxygenase genes were upregulated on oil, mostly transcribed by families Porticoccaceae and Rhodobacteraceae, together with aromatic ring-hydroxylating dioxygenases, mostly transcribed by Rhodobacteraceae. Few transcripts from obligate hydrocarbonoclastic genera of Alcanivorax, Oleispira and Cycloclasticus were significantly enriched in the oil-treated exposed tank in comparison to control the non-exposed tank, and these were mostly transporters and genes involved in nitrogen and phosphorous acquisition. This study highlights the importance of seasonality, i.e., phytoplankton occurrence and senescence leading to organic compound release which can be used preferentially by bacteria over oil compounds, delaying the latter process. As a result, such seasonal effect can reduce the sensitivity of genosensing tools employing bacterial functional genes to sense oil. A better understanding of the use of natural organic matter by bacteria involved in oil-biodegradation is needed to develop an array of functional markers enabling the rapid and specific in situ detection of anthropogenic pollution.