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-m2 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.

Genome-centric view of carbon processing in thawing permafrost.

As global temperatures rise, large amounts of carbon sequestered in permafrost are becoming available for microbial degradation. Accurate prediction of carbon gas emissions from thawing permafrost is limited by our understanding of these microbial communities. Here we use metagenomic sequencing of 214 samples from a permafrost thaw gradient to recover 1,529 metagenome-assembled genomes, including many from phyla with poor genomic representation. These genomes reflect the diversity of this complex ecosystem, with genus-level representatives for more than sixty per cent of the community. Meta-omic analysis revealed key populations involved in the degradation of organic matter, including bacteria whose genomes encode a previously undescribed fungal pathway for xylose degradation. Microbial and geochemical data highlight lineages that correlate with the production of greenhouse gases and indicate novel syntrophic relationships. Our findings link changing biogeochemistry to specific microbial lineages involved in carbon processing, and provide key information for predicting the effects of climate change on permafrost systems.

Evidence for non-methanogenic metabolisms in globally distributed archaeal clades basal to the Methanomassiliicoccales.

Recent discoveries of mcr and mcr-like genes in genomes from diverse archaeal lineages suggest that methane metabolism is an ancient pathway with a complicated evolutionary history. One conventional view is that methanogenesis is an ancestral metabolism of the class Thermoplasmata. Through comparative genomic analysis of 12 Thermoplasmata metagenome-assembled genomes (MAGs) basal to the Methanomassiliicoccales, we show that these microorganisms do not encode the genes required for methanogenesis. Further analysis of 770 Ca. Thermoplasmatota genomes/MAGs found no evidence of mcrA homologues outside of the Methanomassiliicoccales. Together, these results suggest that methanogenesis was laterally acquired by an ancestor of the Methanomassiliicoccales. The 12 analysed MAGs include representatives from four orders basal to the Methanomassiliicoccales, including a high-quality MAG that likely represents a new order, Ca. Lunaplasma lacustris ord. nov. sp. nov. These MAGs are predicted to use diverse energy conservation pathways, including heterotrophy, sulfur and hydrogen metabolism, denitrification, and fermentation. Two lineages are widespread among anoxic, sedimentary environments, whereas Ca. Lunaplasma lacustris has thus far only been detected in alpine caves and subarctic lake sediments. These findings advance our understanding of the metabolic potential, ecology, and global distribution of the Thermoplasmata and provide insight into the evolutionary history of methanogenesis within the Ca. Thermoplasmatota.

Multiplex quantitative PCR for single-reaction genetically modified (GM) plant detection and identification of false-positive GM plants linked to Cauliflower mosaic virus (CaMV) infection.

BACKGROUND: Most genetically modified (GM) plants contain a promoter, P35S, from the plant virus, Cauliflower mosaic virus (CaMV), and many have a terminator, TNOS, derived from the bacterium, Agrobacterium tumefaciens. Assays designed to detect GM plants often target the P35S and/or TNOS DNA sequences. However, because the P35S promoter is derived from CaMV, these detection assays can yield false-positives from non-GM plants infected by this naturally-occurring virus. RESULTS: Here we report the development of an assay designed to distinguish CaMV-infected plants from GM plants in a single multiplexed quantitative PCR (qPCR) reaction. Following initial testing and optimization via PCR and singleplex-to-multiplex qPCR on both plasmid and plant DNA, TaqMan qPCR probes with different fluorescence wavelengths were designed to target actin (a positive-control plant gene), P35S, P3 (a CaMV-specific gene), and TNOS. We tested the specificity of our quadruplex qPCR assay using different DNA extracts from organic watercress and both organic and GM canola, all with and without CaMV infection, and by using commercial and industrial samples. The limit of detection (LOD) of each target was determined to be 1% for actin, 0.001% for P35S, and 0.01% for both P3 and TNOS. CONCLUSIONS: This assay was able to distinguish CaMV-infected plants from GM plants in a single multiplexed qPCR reaction for all samples tested in this study, suggesting that this protocol is broadly applicable and readily transferrable to any interested parties with a qPCR platform.

Genomic resolution of a cold subsurface aquifer community provides metabolic insights for novel microbes adapted to high CO2 concentrations.

As in many deep underground environments, the microbial communities in subsurface high-CO2 ecosystems remain relatively unexplored. Recent investigations based on single-gene assays revealed a remarkable variety of organisms from little studied phyla in Crystal Geyser (Utah, USA), a site where deeply sourced CO2 -saturated fluids are erupted at the surface. To provide genomic resolution of the metabolisms of these organisms, we used a novel metagenomic approach to recover 227 high-quality genomes from 150 microbial species affiliated with 46 different phylum-level lineages. Bacteria from two novel phylum-level lineages have the capacity for CO2 fixation. Analyses of carbon fixation pathways in all studied organisms revealed that the Wood-Ljungdahl pathway and the Calvin-Benson-Bassham Cycle occurred with the highest frequency, whereas the reverse TCA cycle was little used. We infer that this, and selection for form II RuBisCOs, are adaptions to high CO2 -concentrations. However, many autotrophs can also grow mixotrophically, a strategy that confers metabolic versatility. The assignment of 156 hydrogenases to 90 different organisms suggests that H2 is an important inter-species energy currency even under gaseous CO2 -saturation. Overall, metabolic analyses at the organism level provided insight into the biochemical cycles that support subsurface life under the extreme condition of CO2 saturation.

Metagenomic analysis of a high carbon dioxide subsurface microbial community populated by chemolithoautotrophs and bacteria and archaea from candidate phyla.

Research on geologic carbon sequestration raises questions about potential impacts of subsurface microbiota on carbon cycling and biogeochemistry. Subsurface, high-CO2 systems are poorly biologically characterized, partly because of difficulty accessing high-volume, uncontaminated samples. CO2 -driven Crystal Geyser (CG, Utah, USA), an established geologic carbon sequestration analogue, provides high volumes of deep (∼ 200-500 m) subsurface fluids. We explored microbial diversity and metabolic potential in this high-CO2 environment by assembly and analysis of metagenomes recovered from geyser water filtrate. The system is dominated by neutrophilic, iron-oxidizing bacteria, including 'marine' Mariprofundus (Zetaproteobacteria) and 'freshwater' Gallionellales, sulfur-oxidizing Thiomicrospira crunogena and Thiobacillus-like Hydrogenophilales. Near-complete genomes were reconstructed for these bacteria. CG is notably populated by a wide diversity of bacteria and archaea from phyla lacking isolated representatives (candidate phyla) and from as-yet undefined lineages. Many bacteria affiliate with OD1, OP3, OP9, PER, ACD58, WWE3, BD1-5, OP11, TM7 and ZB2. The recovery of nearly 100 genes encoding ribulose-1,5 bisphosphate carboxylase-oxygenase subunit proteins of the Calvin cycle and AMP salvage pathways suggests a strong biological role in high-CO2 subsurface carbon cycling. Overall, we predict microbial impacts on subsurface biogeochemistry via iron, sulfur, and complex carbon oxidation, carbon and nitrogen fixation, fermentation, hydrogen metabolism, and aerobic and anaerobic respiration.

De novo sequences of Haloquadratum walsbyi from Lake Tyrrell, Australia, reveal a variable genomic landscape.

Hypersaline systems near salt saturation levels represent an extreme environment, in which organisms grow and survive near the limits of life. One of the abundant members of the microbial communities in hypersaline systems is the square archaeon, Haloquadratum walsbyi. Utilizing a short-read metagenome from Lake Tyrrell, a hypersaline ecosystem in Victoria, Australia, we performed a comparative genomic analysis of H. walsbyi to better understand the extent of variation between strains/subspecies. Results revealed that previously isolated strains/subspecies do not fully describe the complete repertoire of the genomic landscape present in H. walsbyi. Rearrangements, insertions, and deletions were observed for the Lake Tyrrell derived Haloquadratum genomes and were supported by environmental de novo sequences, including shifts in the dominant genomic landscape of the two most abundant strains. Analysis pertaining to halomucins indicated that homologs for this large protein are not a feature common for all species of Haloquadratum. Further, we analyzed ATP-binding cassette transporters (ABC-type transporters) for evidence of niche partitioning between different strains/subspecies. We were able to identify unique and variable transporter subunits from all five genomes analyzed and the de novo environmental sequences, suggesting that differences in nutrient and carbon source acquisition may play a role in maintaining distinct strains/subspecies.

Impacts of flood damage on airborne bacteria and fungi in homes after the 2013 Colorado Front Range flood.

Flood-damaged homes typically have elevated microbial loads, and their occupants have an increased incidence of allergies, asthma, and other respiratory ailments, yet the microbial communities in these homes remain under-studied. Using culture-independent approaches, we characterized bacterial and fungal communities in homes in Boulder, CO, USA 2-3 months after the historic September, 2013 flooding event. We collected passive air samples from basements in 50 homes (36 flood-damaged, 14 non-flooded), and we sequenced the bacterial 16S rRNA gene (V4-V5 region) and the fungal ITS1 region from these samples for community analyses. Quantitative PCR was used to estimate the abundances of bacteria and fungi in the passive air samples. Results indicate significant differences in bacterial and fungal community composition between flooded and non-flooded homes. Fungal abundances were estimated to be three times higher in flooded, relative to non-flooded homes, but there were no significant differences in bacterial abundances. Penicillium (fungi) and Pseudomonadaceae and Enterobacteriaceae (bacteria) were among the most abundant taxa in flooded homes. Our results suggest that bacterial and fungal communities continue to be affected by flooding, even after relative humidity has returned to baseline levels and remediation has removed any visible evidence of flood damage.

Virus-host and CRISPR dynamics in Archaea-dominated hypersaline Lake Tyrrell, Victoria, Australia.

The study of natural archaeal assemblages requires community context, namely, a concurrent assessment of the dynamics of archaeal, bacterial, and viral populations. Here, we use filter size-resolved metagenomic analyses to report the dynamics of 101 archaeal and bacterial OTUs and 140 viral populations across 17 samples collected over different timescales from 2007-2010 from Australian hypersaline Lake Tyrrell (LT). All samples were dominated by Archaea (75-95%). Archaeal, bacterial, and viral populations were found to be dynamic on timescales of months to years, and different viral assemblages were present in planktonic, relative to host-associated (active and provirus) size fractions. Analyses of clustered regularly interspaced short palindromic repeat (CRISPR) regions indicate that both rare and abundant viruses were targeted, primarily by lower abundance hosts. Although very few spacers had hits to the NCBI nr database or to the 140 LT viral populations, 21% had hits to unassembled LT viral concentrate reads. This suggests local adaptation to LT-specific viruses and/or undersampling of haloviral assemblages in public databases, along with successful CRISPR-mediated maintenance of viral populations at abundances low enough to preclude genomic assembly. This is the first metagenomic report evaluating widespread archaeal dynamics at the population level on short timescales in a hypersaline system.

New approaches indicate constant viral diversity despite shifts in assemblage structure in an Australian hypersaline lake.

It is widely stated that viruses represent the most significant source of biodiversity on Earth, yet characterizing the diversity of viral assemblages in natural systems remains difficult. Viral diversity studies are challenging because viruses lack universally present, phylogenetically informative genes. Here, we developed an approach to estimate viral diversity using a series of functional and novel conserved genes. This approach provides direct estimates of viral assemblage diversity while retaining resolution at the level of individual viral populations in a natural system. We characterized viral assemblages in eight samples from hypersaline Lake Tyrrell (LT), Victoria, Australia, using 39,636 viral contigs. We defined viral operational taxonomic units (OTUs) in two ways. First, we used genes with three different functional predictions that were abundantly represented in the data set. Second, we clustered proteins of unknown function based on sequence similarity, and we chose genes represented by three clusters with numerous members to define OTUs. In combination, diversity metrics indicated between 412 and 735 sampled populations, and the number of populations remained relatively constant across samples. We determined the relative representation of each viral OTU in each sample and found that viral assemblage structures correlate with salinity and solution chemistry. LT viral assemblages were near-replicates from the same site sampled a few days apart but differed significantly on other spatial and temporal scales. The OTU definition approach proposed here paves the way for metagenomics-based analyses of viral assemblages using ecological models previously applied to bacteria and archaea.

Utilizing novel diversity estimators to quantify multiple dimensions of microbial biodiversity across domains.

BACKGROUND: Microbial ecologists often employ methods from classical community ecology to analyze microbial community diversity. However, these methods have limitations because microbial communities differ from macro-organismal communities in key ways. This study sought to quantify microbial diversity using methods that are better suited for data spanning multiple domains of life and dimensions of diversity. Diversity profiles are one novel, promising way to analyze microbial datasets. Diversity profiles encompass many other indices, provide effective numbers of diversity (mathematical generalizations of previous indices that better convey the magnitude of differences in diversity), and can incorporate taxa similarity information. To explore whether these profiles change interpretations of microbial datasets, diversity profiles were calculated for four microbial datasets from different environments spanning all domains of life as well as viruses. Both similarity-based profiles that incorporated phylogenetic relatedness and naïve (not similarity-based) profiles were calculated. Simulated datasets were used to examine the robustness of diversity profiles to varying phylogenetic topology and community composition. RESULTS: Diversity profiles provided insights into microbial datasets that were not detectable with classical univariate diversity metrics. For all datasets analyzed, there were key distinctions between calculations that incorporated phylogenetic diversity as a measure of taxa similarity and naïve calculations. The profiles also provided information about the effects of rare species on diversity calculations. Additionally, diversity profiles were used to examine thousands of simulated microbial communities, showing that similarity-based and naïve diversity profiles only agreed approximately 50% of the time in their classification of which sample was most diverse. This is a strong argument for incorporating similarity information and calculating diversity with a range of emphases on rare and abundant species when quantifying microbial community diversity. CONCLUSIONS: For many datasets, diversity profiles provided a different view of microbial community diversity compared to analyses that did not take into account taxa similarity information, effective diversity, or multiple diversity metrics. These findings are a valuable contribution to data analysis methodology in microbial ecology.

Seasonal variability in bacterial and fungal diversity of the near-surface atmosphere.

Bacteria and fungi are ubiquitous throughout the Earth's lower atmosphere where they often represent an important component of atmospheric aerosols with the potential to impact human health and atmospheric dynamics. However, the diversity, composition, and spatiotemporal dynamics of these airborne microbes remain poorly understood. We performed a comprehensive analysis of airborne microbes across two aerosol size fractions at urban and rural sites in the Colorado Front Range over a 14-month period. Coarse (PM10-2.5) and fine (PM2.5) particulate matter samples were collected at weekly intervals with both bacterial and fungal diversity assessed via high-throughput sequencing. The diversity and composition of the airborne communities varied across the sites, between the two size fractions, and over time. Bacteria were the dominant type of bioaerosol in the collected air samples, while fungi and plants (pollen) made up the remainder, with the relative abundances of fungi peaking during the spring and summer months. As bacteria made up the majority of bioaerosol particles, we analyzed the bacterial communities in greater detail using a bacterial-specific 16S rRNA gene sequencing approach. Overall, bacterial taxonomic richness and the relative abundances of specific bacterial taxa exhibited significant patterns of seasonality. Likewise, airborne bacterial communities varied significantly between sites and across aerosol size fractions. Source-tracking analyses indicate that soils and leaves represented important sources of bacteria to the near-surface atmosphere across all locations with cow fecal bacteria also representing an important source of bioaerosols at the more rural sites during early fall and early spring. Together, these data suggest that a complex set of environmental factors, including changes in atmospheric conditions and shifts in the relative importance of available microbial sources, act to control the composition of microbial bioaerosols in rural and urban environments.

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.

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.

Dynamic viral populations in hypersaline systems as revealed by metagenomic assembly.

Viruses of the Bacteria and Archaea play important roles in microbial evolution and ecology, and yet viral dynamics in natural systems remain poorly understood. Here, we created de novo assemblies from 6.4 Gbp of metagenomic sequence from eight community viral concentrate samples, collected from 12 h to 3 years apart from hypersaline Lake Tyrrell (LT), Victoria, Australia. Through extensive manual assembly curation, we reconstructed 7 complete and 28 partial novel genomes of viruses and virus-like entities (VLEs, which could be viruses or plasmids). We tracked these 35 populations across the eight samples and found that they are generally stable on the timescale of days and transient on the timescale of years, with some exceptions. Cross-detection of the 35 LT populations in three previously described haloviral metagenomes was limited to a few genes, and most previously sequenced haloviruses were not detected in our samples, though 3 were detected upon reducing our detection threshold from 90% to 75% nucleotide identity. Similar results were obtained when we applied our methods to haloviral metagenomic data previously reported from San Diego, CA: 10 contigs that we assembled from that system exhibited a variety of detection patterns on a timescale of weeks to 1 month but were generally not detected in LT. Our results suggest that most haloviral populations have a limited or, possibly, a temporally variable global distribution. This study provides high-resolution insight into viral biogeography and dynamics and it places "snapshot" viral metagenomes, collected at a single time and location, in context.

Gene-centric metagenomics of the fiber-adherent bovine rumen microbiome reveals forage specific glycoside hydrolases.

The complex microbiome of the rumen functions as an effective system for the conversion of plant cell wall biomass to microbial protein, short chain fatty acids, and gases. As such, it provides a unique genetic resource for plant cell wall degrading microbial enzymes that could be used in the production of biofuels. The rumen and gastrointestinal tract harbor a dense and complex microbiome. To gain a greater understanding of the ecology and metabolic potential of this microbiome, we used comparative metagenomics (phylotype analysis and SEED subsystems-based annotations) to examine randomly sampled pyrosequence data from 3 fiber-adherent microbiomes and 1 pooled liquid sample (a mixture of the liquid microbiome fractions from the same bovine rumens). Even though the 3 animals were fed the same diet, the community structure, predicted phylotype, and metabolic potentials in the rumen were markedly different with respect to nutrient utilization. A comparison of the glycoside hydrolase and cellulosome functional genes revealed that in the rumen microbiome, initial colonization of fiber appears to be by organisms possessing enzymes that attack the easily available side chains of complex plant polysaccharides and not the more recalcitrant main chains, especially cellulose. Furthermore, when compared with the termite hindgut microbiome, there are fundamental differences in the glycoside hydrolase content that appear to be diet driven for either the bovine rumen (forages and legumes) or the termite hindgut (wood).

Diversity in the soil virosphere: to infinity and beyond?

Viruses are key members of Earth's microbiomes, shaping microbial community composition and metabolism. Here, we describe recent advances in 'soil viromics', that is, virus-focused metagenome and metatranscriptome analyses that offer unprecedented windows into the soil virosphere. Given the emerging picture of high soil viral activity, diversity, and dynamics over short spatiotemporal scales, we then outline key eco-evolutionary processes that we hypothesize are the major diversity drivers for soil viruses. We argue that a community effort is needed to establish a 'global soil virosphere atlas' that can be used to address the roles of viruses in soil microbiomes and terrestrial biogeochemical cycles across spatiotemporal scales.

Experimental evidence for the impact of soil viruses on carbon cycling during surface plant litter decomposition.

To date, the potential impact of viral communities on biogeochemical cycles in soil has largely been inferred from correlational evidence, such as virus-driven changes in microbial abundances, viral auxiliary metabolic genes, and links with soil physiochemical properties. To more directly test the impact of soil viruses on carbon cycling during plant litter decomposition, we added concentrated viral community suspensions to complex litter decomposer communities in 40-day microcosm experiments. Microbial communities from two New Mexico alpine soils, Pajarito (PJ) and Santa Fe (SF), were inoculated onto grass litter on sand, and three treatments were applied in triplicate to each set of microcosms: addition of buffer (no added virus), live virus (+virus), or killed-virus (+killed-virus) fractions extracted from the same soil. Significant differences in respiration were observed between the +virus and +killed-virus treatments in the PJ, but not the SF microcosms. Bacterial and fungal community composition differed significantly by treatment in both PJ and SF microcosms. Combining data across both soils, viral addition altered links between bacterial and fungal diversity, dissolved organic carbon and total nitrogen. Overall, we demonstrate that increasing viral pressure in complex microbial communities can impact terrestrial biogeochemical cycling but is context-dependent.

Substantial differences in soil viral community composition within and among four Northern California habitats.

Viruses contribute to food web dynamics and nutrient cycles in diverse ecosystems, yet the biogeographical patterns that underlie these viral dynamics are poorly understood, particularly in soil. Here, we identified trends in soil viral community composition in relation to habitat, moisture content, and physical distance. We generated 30 soil viromes from four distinct habitats (wetlands, grasslands, woodlands, and chaparral) by selectively capturing virus-sized particles prior to DNA extraction, and we recovered 3432 unique viral 'species' (dsDNA vOTUs). Viral communities differed significantly by soil moisture content, with viral richness generally higher in wet compared to dry soil habitats. However, vOTUs were rarely shared between viromes, including replicates <10 m apart, suggesting that soil viruses may not disperse well and that future soil viral community sampling strategies may need to account for extreme community differences over small spatial scales. Of the 19% of vOTUs detected in more than one virome, 93% were from the same habitat and site, suggesting greater viral community similarity in closer proximity and under similar environmental conditions. Within-habitat differences indicate that extensive sampling would be required for rigorous cross-habitat comparisons, and results highlight emerging paradigms of high viral activity in wet soils and soil viral community spatial heterogeneity.