Spatial turnover of soil viral populations and genotypes overlain by cohesive responses to moisture in grasslands

Targeted viromes and total metagenomes capture distinct components of bee gut phage communities

BACKGROUND

Despite being among the most abundant biological entities on earth, bacteriophage (phage) remain an understudied component of host-associated systems. One limitation to studying host-associated phage is the lack of consensus on methods for sampling phage communities. Here, we compare paired total metagenomes and viral size fraction metagenomes (viromes) as methods for investigating the dsDNA viral communities associated with the GI tract of two bee species: the European honey bee Apis mellifera and the eastern bumble bee Bombus impatiens.

RESULTS

We find that viromes successfully enriched for phage, thereby increasing phage recovery, but only in honey bees. In contrast, for bumble bees, total metagenomes recovered greater phage diversity. Across both bee species, viromes better sampled low occupancy phage, while total metagenomes were biased towards sampling temperate phage. Additionally, many of the phage captured by total metagenomes were absent altogether from viromes. Comparing between bees, we show that phage communities in commercially reared bumble bees are significantly reduced in diversity compared to honey bees, likely reflecting differences in bacterial titer and diversity. In a broader context, these results highlight the complementary nature of total metagenomes and targeted viromes, especially when applied to host-associated environments.

CONCLUSIONS

Overall, we suggest that studies interested in assessing total communities of host-associated phage should consider using both approaches. However, given the constraints of virome sampling, total metagenomes may serve to sample phage communities with the understanding that they will preferentially sample dominant and temperate phage. Video Abstract.

The direct and indirect drivers shaping RNA viral communities in grassland soils

Recent studies have revealed diverse RNA viral communities in soils. Yet, how environmental factors influence soil RNA viruses remains largely unknown. Here, we recovered RNA viral communities from bulk metatranscriptomes sequenced from grassland soils managed for 5 years under multiple environmental conditions including water content, plant presence, cultivar type, and soil depth. More than half of the unique RNA viral contigs (64.6%) were assigned with putative hosts. About 74.7% of these classified RNA viral contigs are known as eukaryotic RNA viruses suggesting eukaryotic RNA viruses may outnumber prokaryotic RNA viruses by nearly three times in this grassland. Of the identified eukaryotic RNA viruses and the associated eukaryotic species, the most dominant taxa were Mitoviridae with an average relative abundance of 72.4%, and their natural hosts, Fungi with an average relative abundance of 56.6%. Network analysis and structural equation modeling support that soil water content, plant presence, and type of cultivar individually demonstrate a significant positive impact on eukaryotic RNA viral richness directly as well as indirectly on eukaryotic RNA viral abundance via influencing the co-existing eukaryotic members. A significant negative influence of soil depth on soil eukaryotic richness and abundance indirectly impacts soil eukaryotic RNA viral communities. These results provide new insights into the collective influence of multiple environmental and community factors that shape soil RNA viral communities and offer a structured perspective of how RNA virus diversity and ecology respond to environmental changes.

IMPORTANCE

Climate change has been reshaping the soil environment as well as the residing microbiome. This study provides field-relevant information on how environmental and community factors collectively shape soil RNA communities and contribute to ecological understanding of RNA viral survival under various environmental conditions and virus-host interactions in soil. This knowledge is critical for predicting the viral responses to climate change and the potential emergence of biothreats.

Effects of vegetation cover and aquaculture pollution on viral assemblages in mangroves sediments

Mangrove forests, a critical coastal ecosystem, face numerous anthropogenic threats, particularly from aquaculture activities. Despite the acknowledged significance of viruses in local and global biogeochemical cycles, there is limited knowledge regarding the community structure, genomic diversity, and ecological roles of viruses in mangrove forests ecosystems, especially regarding their responses to aquaculture. In this study, we identified 17,755 viral operational taxonomic units (vOTUs) from nine sediments viromes across three distinct ecological regions of the mangrove forests ecosystem: mangrove, bare flat, and aquaculture regions. Viral assemblages varied among three regions, and the pathogenic viruses associated with marine animals, such as the white spot syndrome virus (WSSV) from Nimaviridae, were identified in this study. The relative abundance of Nimaviridae in the bare flat region was higher than in other regions. Furthermore, viruses in distinct mangrove forests sediments regions have adapted to their environments by adopting distinct survival strategies and encoding various auxiliary metabolic genes involved in carbon metabolism and antibiotic resistance. These adaptations may have profound impacts on biogeochemical cycles. This study provides the first insights into the effects of vegetation cover and aquaculture on the community structure and ecological roles of viruses in mangrove forests sediments. These findings are crucial for understanding the risks posed by anthropogenic threats to mangrove forests ecosystems and informing effective management strategies.

A global atlas of soil viruses reveals unexplored biodiversity and potential biogeochemical impacts

Historically neglected by microbial ecologists, soil viruses are now thought to be critical to global biogeochemical cycles. However, our understanding of their global distribution, activities and interactions with the soil microbiome remains limited. Here we present the Global Soil Virus Atlas, a comprehensive dataset compiled from 2,953 previously sequenced soil metagenomes and composed of 616,935 uncultivated viral genomes and 38,508 unique viral operational taxonomic units. Rarefaction curves from the Global Soil Virus Atlas indicate that most soil viral diversity remains unexplored, further underscored by high spatial turnover and low rates of shared viral operational taxonomic units across samples. By examining genes associated with biogeochemical functions, we also demonstrate the viral potential to impact soil carbon and nutrient cycling. This study represents an extensive characterization of soil viral diversity and provides a foundation for developing testable hypotheses regarding the role of the virosphere in the soil microbiome and global biogeochemistry.

Habitat- and lifestyle-dependent structural and functional characteristics of viruses in mangrove wetlands of different functional zonings

Mangrove wetlands, as one of the natural ecosystems with the most ecological services, have garnered widespread attention about their microbial driven biogeochemical cycling. Urbanization have led to different spatial patterns of environmental conditions and microbial communities in mangroves. However, viruses, as the pivotal drivers of biogeochemical cycling in mangroves, remain inadequately explored in terms of how their ecological potential and complex interactions with host respond to functional zonings. To address this knowledge gap, we conducted a comprehensive investigation on the structural and functional properties of temperate and lytic viruses in mangrove wetlands from different functional zonings by jointly using high-throughput sequencing, prokaryotic and viral metagenomics. Multiple environmental factors were found to significantly influence the taxonomic and functional composition, as well as lysogen-lysis decision-making of mangrove viruses. Furthermore, enriched auxiliary metabolic genes (AMGs) involved in methane, nitrogen and sulfur metabolism, and heavy metal resistance were unveiled in mangrove viruses, whose community composition was closely related to lifestyle and host. The virus-host pairs with different lifestyles were also discovered to react to environmental changes in different ways, which provided an empirical evidence for how virus and bacteria dynamics were specific to viral lifestyles in nature. This study expands our comprehension of the intricate interactions among virus, prokaryotic host and the environment in mangrove wetlands from multiple perspectives, including viral lifestyles, virus-host interactions, and habitat dependence. Importantly, it provides a new ecological perspective on how mangrove viruses are adapted to the stress posed by urbanization.

Host and nonhost bacteria support bacteriophage dissemination along mycelia and abiotic dispersal networks

Bacteriophages play a crucial role in shaping bacterial communities, yet the mechanisms by which nonmotile bacteriophages interact with their hosts remain poorly understood. This knowledge gap is especially pronounced in structured environments like soil, where spatial constraints and air-filled zones hinder aqueous diffusion. In soil, hyphae of filamentous microorganisms form a network of 'fungal highways' (FHs) that facilitate the dispersal of other microorganisms. We propose that FHs also promote bacteriophage dissemination. Viral particles can diffuse in liquid films surrounding hyphae or be transported by infectable (host) or uninfectable (nonhost) bacterial carriers coexisting on FH networks. To test this, two bacteriophages that infect Pseudomonas putida DSM291 (host) but not KT2440 (nonhost) were used. In the absence of carriers, bacteriophages showed limited diffusion on 3D-printed abiotic networks, but diffusion was significantly improved in Pythium ultimum-formed FHs when the number of connecting hyphae exceeded 20. Transport by both host and nonhost carriers enhanced bacteriophage dissemination. Host carriers were five times more effective in transporting bacteriophages, particularly in FHs with over 30 connecting hyphae. This study enhances our understanding of bacteriophage dissemination in nonsaturated environments like soils, highlighting the importance of biotic networks and bacterial hosts in facilitating this process.

Investigating eco-evolutionary processes of microbial community assembly in the wild using a model leaf litter system

Microbial communities are not the easiest to manipulate experimentally in natural ecosystems. However, leaf litter-topmost layer of surface soil-is uniquely suitable to investigate the complexities of community assembly. Here, we reflect on over a decade of collaborative work to address this topic using leaf litter as a model system in Southern California ecosystems. By leveraging a number of methodological advantages of the system, we have worked to demonstrate how four processes-selection, dispersal, drift, and diversification-contribute to bacterial and fungal community assembly and ultimately impact community functioning. Although many dimensions remain to be investigated, our initial results demonstrate that both ecological and evolutionary processes occur simultaneously to influence microbial community assembly. We propose that the development of additional and experimentally tractable microbial systems will be enormously valuable to test the role of eco-evolutionary processes in natural settings and their implications in the face of rapid global change.

Viral metagenomics reveals diverse virus-host interactions throughout the soil depth profile

IMPORTANCE

Soil viruses can moderate the roles that their host microbes play in global carbon cycling. However, given that most studies investigate the surface layer (i.e., top 20 cm) of soil, the extent to which this occurs in subsurface soil (i.e., below 20 cm) is unknown. Here, we leveraged public sequencing data to investigate the interactions between viruses and their hosts at soil depth intervals, down to 115 cm. While most viruses were detected throughout the soil depth profile, their adaptation to host microbes varied. Nonetheless, we uncovered evidence for the potential of soil viruses to encourage their hosts to recycle plant-derived carbon in both surface and subsurface soils. This work reasons that our understanding of soil viral functions requires us to continue to dig deeper and compare viruses existing throughout soil ecosystems.

Viromes vs. mixed community metagenomes: choice of method dictates interpretation of viral community ecology

Background

Viruses, the majority of which are uncultivated, are among the most abundant biological entities on Earth. From altering microbial physiology to driving community dynamics, viruses are fundamental members of microbiomes. While the number of studies leveraging viral metagenomics (viromics) for studying uncultivated viruses is growing, standards for viromics research are lacking. Viromics can utilize computational discovery of viruses from total metagenomes of all community members (hereafter metagenomes) or use physical separation of virus-specific fractions (hereafter viromes). However, differences in the recovery and interpretation of viruses from metagenomes and viromes obtained from the same samples remain understudied.

Results

Here, we compare viral communities from paired viromes and metagenomes obtained from 60 diverse samples across human gut, soil, freshwater, and marine ecosystems. Overall, viral communities obtained from viromes were more abundant and species rich than those obtained from metagenomes, although there were some exceptions. Despite this, metagenomes still contained many viral genomes not detected in viromes. We also found notable differences in the predicted lytic state of viruses detected in viromes vs metagenomes at the time of sequencing. Other forms of variation observed include genome presence/absence, genome quality, and encoded protein content between viromes and metagenomes, but the magnitude of these differences varied by environment.

Conclusions

Overall, our results show that the choice of method can lead to differing interpretations of viral community ecology. We suggest that the choice of whether to target a metagenome or virome to study viral communities should be dependent on the environmental context and ecological questions being asked. However, our overall recommendation to researchers investigating viral ecology and evolution is to pair both approaches to maximize their respective benefits.

Zea mays genotype influences microbial and viral rhizobiome community structure

Plant genotype is recognized to contribute to variations in microbial community structure in the rhizosphere, soil adherent to roots. However, the extent to which the viral community varies has remained poorly understood and has the potential to contribute to variation in soil microbial communities. Here we cultivated replicates of two Zea mays genotypes, parviglumis and B73, in a greenhouse and harvested the rhizobiome (rhizoplane and rhizosphere) to identify the abundance of cells and viruses as well as rhizobiome microbial and viral community using 16S rRNA gene amplicon sequencing and genome resolved metagenomics. Our results demonstrated that viruses exceeded microbial abundance in the rhizobiome of parviglumis and B73 with a significant variation in both the microbial and viral community between the two genotypes. Of the viral contigs identified only 4.5% (n = 7) of total viral contigs were shared between the two genotypes, demonstrating that plants even at the level of genotype can significantly alter the surrounding soil viral community. An auxiliary metabolic gene associated with glycoside hydrolase (GH5) degradation was identified in one viral metagenome-assembled genome (vOTU) identified in the B73 rhizobiome infecting Propionibacteriaceae (Actinobacteriota) further demonstrating the viral contribution in metabolic potential for carbohydrate degradation and carbon cycling in the rhizosphere. This variation demonstrates the potential of plant genotype to contribute to microbial and viral heterogeneity in soil systems and harbors genes capable of contributing to carbon cycling in the rhizosphere.

Dispersal, habitat filtering, and eco-evolutionary dynamics as drivers of local and global wetland viral biogeography

Wetlands store 20-30% of the world's soil carbon, and identifying the microbial controls on these carbon reserves is essential to predicting feedbacks to climate change. Although viral infections likely play important roles in wetland ecosystem dynamics, we lack a basic understanding of wetland viral ecology. Here 63 viral size-fraction metagenomes (viromes) and paired total metagenomes were generated from three time points in 2021 at seven fresh- and saltwater wetlands in the California Bodega Marine Reserve. We recovered 12,826 viral population genomic sequences (vOTUs), only 4.4% of which were detected at the same field site two years prior, indicating a small degree of population stability or recurrence. Viral communities differed most significantly among the seven wetland sites and were also structured by habitat (plant community composition and salinity). Read mapping to a new version of our reference database, PIGEONv2.0 (515,763 vOTUs), revealed 196 vOTUs present over large geographic distances, often reflecting shared habitat characteristics. Wetland vOTU microdiversity was significantly lower locally than globally and lower within than between time points, indicating greater divergence with increasing spatiotemporal distance. Viruses tended to have broad predicted host ranges via CRISPR spacer linkages to metagenome-assembled genomes, and increased SNP frequencies in CRISPR-targeted major tail protein genes suggest potential viral eco-evolutionary dynamics in response to both immune targeting and changes in host cell receptors involved in viral attachment. Together, these results highlight the importance of dispersal, environmental selection, and eco-evolutionary dynamics as drivers of local and global wetland viral biogeography.

Viral but not bacterial community successional patterns reflect extreme turnover shortly after rewetting dry soils

As central members of soil trophic networks, viruses have the potential to drive substantial microbial mortality and nutrient turnover. Pinpointing viral contributions to terrestrial ecosystem processes remains a challenge, as temporal dynamics are difficult to unravel in the spatially and physicochemically heterogeneous soil environment. In Mediterranean grasslands, the first rainfall after seasonal drought provides an ecosystem reset, triggering microbial activity during a tractable window for capturing short-term dynamics. Here, we simulated precipitation in microcosms from four distinct dry grassland soils and generated 144 viromes, 84 metagenomes and 84 16S ribosomal RNA gene amplicon datasets to characterize viral, prokaryotic and relic DNA dynamics over 10 days. Vastly different viral communities in each soil followed remarkably similar successional trajectories. Wet-up triggered a significant increase in viral richness, followed by extensive compositional turnover. Temporal succession in prokaryotic communities was much less pronounced, perhaps suggesting differences in the scales of activity captured by viromes (representing recently produced, ephemeral viral particles) and total DNA. Still, differences in the relative abundances of Actinobacteria (enriched in dry soils) and Proteobacteria (enriched in wetted soils) matched those of their predicted phages, indicating viral predation of dominant bacterial taxa. Rewetting also rapidly depleted relic DNA, which subsequently reaccumulated, indicating substantial new microbial mortality in the days after wet-up, particularly of the taxa putatively under phage predation. Production of abundant, diverse viral particles via microbial host cell lysis appears to be a conserved feature of the early response to soil rewetting, and results suggest the potential for 'Cull-the-Winner' dynamics, whereby viruses infect and cull but do not decimate dominant host populations.

How might bacteriophages shape biological invasions?

Invasions by eukaryotes dependent on environmentally acquired bacterial mutualists are often limited by the ability of bacterial partners to survive and establish free-living populations. Focusing on the model legume-rhizobium mutualism, we apply invasion biology hypotheses to explain how bacteriophages can impact the competitiveness of introduced bacterial mutualists. Predicting how phage-bacteria interactions affect invading eukaryotic hosts requires knowing the eco-evolutionary constraints of introduced and native microbial communities, as well as their differences in abundance and diversity. By synthesizing research from invasion biology, as well as bacterial, viral, and community ecology, we create a conceptual framework for understanding and predicting how phages can affect biological invasions through their effects on bacterial mutualists.

A subset of viruses thrives following microbial resuscitation during rewetting of a seasonally dry California grassland soil

Viruses are abundant, ubiquitous members of soil communities that kill microbial cells, but how they respond to perturbation of soil ecosystems is essentially unknown. Here, we investigate lineage-specific virus-host dynamics in grassland soil following "wet-up", when resident microbes are both resuscitated and lysed after a prolonged dry period. Quantitative isotope tracing, time-resolved metagenomics and viromic analyses indicate that dry soil holds a diverse but low biomass reservoir of virions, of which only a subset thrives following wet-up. Viral richness decreases by 50% within 24 h post wet-up, while viral biomass increases four-fold within one week. Though recent hypotheses suggest lysogeny predominates in soil, our evidence indicates that viruses in lytic cycles dominate the response to wet-up. We estimate that viruses drive a measurable and continuous rate of cell lysis, with up to 46% of microbial death driven by viral lysis one week following wet-up. Thus, viruses contribute to turnover of soil microbial biomass and the widely reported CO2 efflux following wet-up of seasonally dry soils.

Enumerating soil biodiversity

Soil is an immense habitat for diverse organisms across the tree of life, but just how many organisms live in soil is surprisingly unknown. Previous efforts to enumerate soil biodiversity consider only certain types of organisms (e.g., animals) or report values for diverse groups without partitioning species that live in soil versus other habitats. Here, we reviewed the biodiversity literature to show that soil is likely home to 59 ± 15% of the species on Earth. We therefore estimate an approximately two times greater soil biodiversity than previous estimates, and we include representatives from the simplest (microbial) to most complex (mammals) organisms. Enchytraeidae have the greatest percentage of species in soil (98.6%), followed by fungi (90%), Plantae (85.5%), and Isoptera (84.2%). Our results demonstrate that soil is the most biodiverse singular habitat. By using this estimate of soil biodiversity, we can more accurately and quantitatively advocate for soil organismal conservation and restoration as a central goal of the Anthropocene.

Meeting report: The first soil viral workshop 2022

Soil viral ecology is a growing research field; however, the state of knowledge still lags behind that of aquatic systems. Therefore, to facilitate progress, the first Soil Viral Workshop was held to encourage international scientific discussion and collaboration, suggest guidelines for future research, and establish soil viral research as a concrete research area. The workshop took place at Søminestationen, Denmark, between 15 and 17th of June 2022. The meeting was primarily held in person, but the sessions were also streamed online. The workshop was attended by 23 researchers from ten different countries and from a wide range of subfields and career stages. Eleven talks were presented, followed by discussions revolving around three major topics: viral genomics, virus-host interactions, and viruses in the soil food web. The main take-home messages and suggestions from the discussions are summarized in this report.

Highly diverse and unknown viruses may enhance Antarctic endoliths' adaptability

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Propagation of viral genomes by replicating ammonia-oxidising archaea during soil nitrification

Ammonia-oxidising archaea (AOA) are a ubiquitous component of microbial communities and dominate the first stage of nitrification in some soils. While we are beginning to understand soil virus dynamics, we have no knowledge of the composition or activity of those infecting nitrifiers or their potential to influence processes. This study aimed to characterise viruses having infected autotrophic AOA in two nitrifying soils of contrasting pH by following transfer of assimilated CO₂-derived ¹³C from host to virus via DNA stable-isotope probing and metagenomic analysis. Incorporation of ¹³C into low GC mol% AOA and virus genomes increased DNA buoyant density in CsCl gradients but resulted in co-migration with dominant non-enriched high GC mol% genomes, reducing sequencing depth and contig assembly. We therefore developed a hybrid approach where AOA and virus genomes were assembled from low buoyant density DNA with subsequent mapping of ¹³C isotopically enriched high buoyant density DNA reads to identify activity of AOA. Metagenome-assembled genomes were different between the two soils and represented a broad diversity of active populations. Sixty-four AOA-infecting viral operational taxonomic units (vOTUs) were identified with no clear relatedness to previously characterised prokaryote viruses. These vOTUs were also distinct between soils, with 42% enriched in ¹³C derived from hosts. The majority were predicted as capable of lysogeny and auxiliary metabolic genes included an AOA-specific multicopper oxidase suggesting infection may augment copper uptake essential for central metabolic functioning. These findings indicate virus infection of AOA may be a frequent process during nitrification with potential to influence host physiology and activity.

Spatial turnover of soil viral populations and genotypes overlain by cohesive responses to moisture in grasslands

Viruses shape microbial communities, food web dynamics, and carbon and nutrient cycling in diverse ecosystems. However, little is known about the patterns and drivers of viral community composition, particularly in soil, precluding a predictive understanding of viral impacts on terrestrial habitats. To investigate soil viral community assembly processes, here we analyzed 43 soil viromes from a rainfall manipulation experiment in a Mediterranean grassland in California. We identified 5,315 viral populations (viral operational taxonomic units [vOTUs] with a representative sequence ≥10 kbp) and found that viral community composition exhibited a highly significant distance-decay relationship within the 200-m² field site. This pattern was recapitulated by the intrapopulation microheterogeneity trends of prevalent vOTUs (detected in ≥90% of the viromes), which tended to exhibit negative correlations between spatial distance and the genomic similarity of their predominant allelic variants. Although significant spatial structuring was also observed in the bacterial and archaeal communities, the signal was dampened relative to the viromes, suggesting differences in local assembly drivers for viruses and prokaryotes and/or differences in the temporal scales captured by viromes and total DNA. Despite the overwhelming spatial signal, evidence for environmental filtering was revealed in a protein-sharing network analysis, wherein a group of related vOTUs predicted to infect actinobacteria was shown to be significantly enriched in low-moisture samples distributed throughout the field. Overall, our results indicate a highly diverse, dynamic, active, and spatially structured soil virosphere capable of rapid responses to changing environmental conditions.