Moisture modulates soil reservoirs of active DNA and RNA viruses

Wu, Ruonan; Davison, Michelle R.; Gao, Yuqian; Nicora, Carrie D.; Mcdermott, Jason E.; Burnum-Johnson, Kristin E.; Hofmockel, Kirsten S.; Jansson, Janet K.

doi:10.1038/s42003-021-02514-2

Download PDF

Article
Open access
Published: 26 August 2021

Moisture modulates soil reservoirs of active DNA and RNA viruses

Communications Biology volume 4, Article number: 992 (2021) Cite this article

4180 Accesses
26 Citations
66 Altmetric
Metrics details

Subjects

Abstract

Soil is known to harbor viruses, but the majority are uncharacterized and their responses to environmental changes are unknown. Here, we used a multi-omics approach (metagenomics, metatranscriptomics and metaproteomics) to detect active DNA viruses and RNA viruses in a native prairie soil and to determine their responses to extremes in soil moisture. The majority of transcribed DNA viruses were bacteriophage, but some were assigned to eukaryotic hosts, mainly insects. We also demonstrated that higher soil moisture increased transcription of a subset of DNA viruses. Metaproteome data validated that the specific viral transcripts were translated into proteins, including chaperonins known to be essential for virion replication and assembly. The soil viral chaperonins were phylogenetically distinct from previously described marine viral chaperonins. The soil also had a high abundance of RNA viruses, with highest representation of Reoviridae. Leviviridae were the most diverse RNA viruses in the samples, with higher amounts in wet soil. This study demonstrates that extreme shifts in soil moisture have dramatic impacts on the composition, activity and potential functions of both DNA and RNA soil viruses.

RNA-viromics reveals diverse communities of soil RNA viruses with the potential to affect grassland ecosystems across multiple trophic levels

Article Open access 08 April 2022

Luke S. Hillary, Evelien M. Adriaenssens, … James E. McDonald

Viromes outperform total metagenomes in revealing the spatiotemporal patterns of agricultural soil viral communities

Article Open access 21 February 2021

Christian Santos-Medellin, Laura A. Zinke, … Joanne B. Emerson

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

Article 28 September 2023

Christian Santos-Medellín, Steven J. Blazewicz, … Joanne B. Emerson

Introduction

Viruses are ubiquitous and abundant in nature¹. However, they are underexplored in soil environments due to the experimental and computational difficulties in extracting and decoding their genetic complexity^2,3,4,5,6. With increasing access to large sequence datasets, i.e., metagenomes and metatranscriptomes, it has become possible to perform more high-throughput and detailed surveys of viruses in soil environments^2,4.

Metagenomes, while valuable for assembling DNA viral genomes, do not reveal whether the viruses detected are actively contributing to ecosystem functions at the time of collection. By contrast, study of viral community expression through metatranscriptomics and metaproteomics has the potential to reveal which viruses are active and what genes are expressed under given environmental conditions. In addition, metatranscriptome data can also be used to reassemble RNA viruses. Our knowledge of soil RNA viruses is currently even more limited than for DNA viruses, due to greater challenges in extraction of viral RNA from natural environments⁷. However, in the time of COVID-19, the knowledge of environmental RNA viruses is of utmost importance, given that soil and other natural ecosystems could function as environmental reservoirs for potentially pathogenic RNA viruses. To date, there has been only one survey of metatranscriptomes for RNA viruses in soil microcosms³, but their types and abundances in natural field sites and in other soil types are yet to be investigated. Metaproteomics provides an additional validation of viral activity, including which viral proteins are produced and evidence of lytic viral lifestyles. Metaproteomics is, however, challenging to perform in soil and there have been few successful examples to date^8,9,10. However, improving techniques and new methodologies in proteomics have recently made it possible to delve deeper into soil proteomes for the detection of less-abundant proteins^11,12. Therefore, here we used a multi-omics approach, combining DNA and RNA sequencing with metaproteomics, leveraging the strength of each omics method, to better understand how soil viruses respond to environmental changes.

In this study, we aimed to address how differences in soil moisture influence the activities of soil DNA and RNA viruses. We focused our study on a native prairie soil at the Konza Experimental Station in Kansas, which sits at the crossroads for predicted shifts in precipitation with climate change, either increasing drought towards the southwest, or increasing rainfall to the northeast¹³. To investigate both possible scenarios, we hypothesized that wet and dry soils would harbor different soil viruses and that their respective activities would be influenced by soil moisture. To address this hypothesis, we leveraged existing data from a simple experiment that compared Konza prairie soil that was saturated to soil that was desiccated^8,13,14,15. Briefly, three replicate soil samples from two independent field locations in the Konza prairie were either wet to saturation or air-dried to represent wet and dry soil treatments, respectively¹³, resulting in a total of 12 soil microcosms. After a 15-day incubation at 21 °C under wet or dry conditions, DNA, RNA, and protein were extracted from each of the soil microcosms and the expression data (metatranscriptome and metaproteome) were mapped to metagenomes that were constructed from the original soil prior to incubation. The DNA was used for 16S rRNA gene sequencing as previously described¹³. The resultant multi-omics data were then used to screen for soil viruses and to assess their activities in response to differences in soil moisture.

Results and discussion

A diverse and active DNA virosphere

We first leveraged two existing metagenomes that were constructed from the Konza native prairie soil^14,15 to screen for viral sequences at the site. Each of the metagenomes was obtained from a composite of all the replicate soils collected at ambient field moisture conditions. One of the metagenomes was de novo assembled from deep sequence data (1.1 Tb)¹⁴ and the second was a hybrid assembly of short and long reads (267.0 Gb)¹⁶. The combination of the two metagenomes was used to maximize the coverage of viral sequences from the Konza prairie site. To balance between the detection limits of the viral detection tools and the wide range of viral genome size, the viral contigs > 2.5 kb in length were combined with those obtained from screening of the two largest public viral databases (i.e., IMG/VR¹⁷ and NCBI Virus¹⁶) to further increase the coverage of DNA viral sequences. We acknowledge that the length cutoff of 2.5 kb would preclude detection of some ssDNA viruses with small segmented genome sizes (e.g., Nanoviridae¹⁸). As a result, a DNA viral database for the site was curated that included 726,108 de-replicated viral contigs. The DNA viral database then served as a scaffold for mapping of metatranscriptome and metaproteome datasets to determine the activities of soil DNA viruses and their responses to differences in soil moisture. This approach was also recently applied to detect the transcriptional activity of marine prokaryotic and eukaryotic viruses^19,20,21,22 and giant viruses in soil⁵.

The metatranscriptome reads from both wet and dry treatments were mapped to a total of 416 unique DNA viral contigs using stringent criteria (% sequence identity > 95% and % sequence coverage > 80%). The 416 DNA viral contigs with an average sequence length of 19 kb were highly diverse and grouped into 139 clusters, with 111 of the clusters being singletons (Supplementary Data 1).

We aimed to assign putative host taxa to the viral clusters by combining several approaches: CRISPR spacer matching, and screening for host and viral sequence similarities to respective databases (details in ‘Methods’). As a result, we assigned putative viral host taxa to 160 out of the 416 transcribed DNA viral contigs. Some of these were assigned to more than one host (Supplementary Data 1), resulting in a total of 181 virus–host pairings (Fig. 1a). Of these, 79 host–virus pairs were detected only in the dry soil treatment, 51 were only in the wet soil treatment, and an additional 51 were found in both dry and wet treatments (Fig. 1a). Consistent with previous reports⁴, the majority of the transcribed DNA viral contigs were annotated as bacteriophage sequences. Different sets of transcribed DNA viral contigs were unique to wet or dry soils and assigned to specific hosts at the phylum level, whereas others were shared (Fig. 1a). However, the dominant soil taxa, i.e., Proteobacteria and Actinobacteria that were previously identified by 16S rRNA gene sequencing in this soil environment, were predicted as hosts under both wet and dry conditions (Supplementary Fig. 1a). Eukaryotic DNA viruses, such as Bracovirus and Ichnovirus belonging to a family of insect viruses within the Polydnaviridae family, were also transcribed in the soils (Fig. 1a and Supplementary Data 1). Most of these insect viruses were only detected in dry soil conditions. These differences in virus–host pairings suggest that some of the respective hosts were impacted differently by the dry and wet incubation conditions.

**Fig. 1: Transcribed DNA viral communities and their responses to wet and dry soil conditions.**

There were 21 DNA viral contigs that were assigned to hosts across multiple bacterial phyla suggesting the presence of viral generalists^1,23 (Supplementary Data 1). We recognize that host assignment based on CRISPR spacer matching, however, is limited to detection of recent or historical virus–host interactions that were captured at the time of sampling²⁴. As bioinformatics assignment of virus–host linkages only suggests possible pairings based on sequence features, there are also chances of introducing false positives. However, we applied the most stringent criteria possible to provide confident host assignments.

Increased activity of a subset of DNA viruses in wet soil

Soil moisture has a strong influence on the community structures of transcribed DNA viruses. The majority of the transcriptionally active DNA viral contigs were unique to wet or dry conditions, with only 111 viral contigs (~ 26.7%) detected in both wet and dry soils, suggesting that the different soil moisture conditions may shape the activity of the DNA viral community differently (Fig. 1b). Interestingly, although a significantly higher number of transcribed DNA viral contigs were detected in dry soils (Fig. 1b, c), the levels of transcriptional activity were significantly higher (based on the normalized abundance of RNA reads that mapped to the viral contigs) for DNA viruses in wet soils irrespective of sampling site location (Fig. 1d). DNA viral contigs with mapped transcripts could represent either prophages that are passively replicated along with their host genomes, or (lytic) viruses that are actively regulating early/middle/late expression of viral gene clusters²⁵. In soil, a lysogenic lifestyle is considered to be an adaptive strategy for viruses to cope with long periods of low host activity^26,27. Therefore, the 1.5-fold increase in the number of transcribed DNA viral contigs representing transcriptionally active DNA viruses, but with lower levels of overall transcription, in dry soil suggests that the increase was due to a higher prevalence of lysogeny in dry conditions. This hypothesis is strengthened by our finding of a 20-fold increase in transcripts for lysogenic markers (i.e., integrase and excisionase) in one of our replicates (A-2) in dry compared to wet conditions (Supplementary Data 2). High number of lysogenic phages were also previously reported in dry Antarctic soils using a cultivation-independent induction assay²⁸. By contrast, under wet soil conditions we found a 2-fold increase in transcription of fewer viral contigs representing a subset of DNA viruses, suggesting that those viruses were more transcriptionally active in response to higher soil moisture. In addition, there was a higher correlation between prokaryotic abundances, as estimated by 16S rRNA gene sequencing, with DNA viral transcript counts in wet soils (R² = 0.593, Supplementary Fig. 1d) in comparison to dry soils (R² = 0.069, Supplementary Fig. 1d), supporting this hypothesis.

We then identified which soil DNA viruses were most transcriptionally active and how they responded to the differences in soil moisture. As the majority of the transcribed DNA viral contigs (97%) were environmental viruses with unclassified taxonomy assignment, we were not able to calculate the taxonomic abundance of each and instead compared the differential abundances of the transcribed viral contigs. There were four DNA viral contigs with significantly different levels of transcription under wet and dry conditions (VC_1, VC_19, VC_282, VC_412; Fig. 2a). Contigs VC_1 and VC_19 correspond to unclassified viral contigs deposited in IMG/VR (identifiers of ‘REF:2547132004_2547132004’ and ‘3300010038_Ga0126315_10000854’) that were previously detected in metagenomes from the Rifle site²⁹ and from serpentine soil in the UC McLaughlin Reserve³⁰, respectively. Contigs VC_282 and VC_412 were extracted from our Kansas metagenomes. Contigs VC_1 and VC_19 had significantly higher levels of transcriptional activity in wet soils compared to dry soils (p < 0.01, Fig. 2b), whereas VC_412 had significantly more transcripts in dry soils (p < 0.01, Fig. 2b). However, VC_282 was detected with higher transcript levels only in the dry soil replicates from site C. Five specific regions of contig VC_1 had the highest transcriptional activity in both wet and dry soil conditions, with lengths of 546, 338, 175, 420, and 663 bp, respectively (Fig. 2c). The finding that the same specific regions had differential transcript mapping frequencies is suggestive of active/lytic viruses with highly regulated transcription of early and/or late genes^25,31, in comparison to lysogenic viruses that are passively transcribed along with their host genomes. VC_1 was originally detected from the Rifle site, an aquifer environment. Therefore, this virus–host pair may be better adapted to wet conditions, reflecting our finding of higher transcription levels for VC-1 in wet soil (p < 0.05, Fig. 2b). A similar increase in activity of DNA viruses together with a bloom of their respective hosts has also previously been observed following laboratory wetting of soil biocrusts³², suggesting that this may be a common phenomenon to soil wetting.

**Fig. 2: DNA viral contigs with differential transcription in wet and dry soil treatments.**

It is interesting to note that the quality-filtered transcripts were mapped to both protein coding and noncoding regions of the DNA viral contigs. Noncoding RNAs with a phage origin have previously been reported to regulate viral replication at the early stage of infection and to maintain a lysogenic state by silencing the expression of late structural genes^33,34. We observed a trend towards a higher percentage of viral noncoding RNAs in drier soils from site A, along with higher transcriptional levels of lysogenic markers (i.e., integrase and excisionase) in these samples (Fig. 3a and Supplementary Data 2). These findings suggest higher levels of lysogenic phages in these samples, but this hypothesis needs further experimental validation. Interestingly, recent studies identified prophage-encoded noncoding RNAs that can also contribute to the virulence of bacterial hosts³⁵ and protective functions such as superinfection-immunity³⁶. Future studies are therefore needed to better characterize the functions of viral noncoding RNAs.

**Fig. 3: Functional characterization of viral transcripts and proteins.**

After removal of noncoding RNAs, the remaining transcripts were mapped to 314 viral genes. Only 149 of those transcribed viral genes were annotated to functional gene categories (Supplementary Data 3), reflecting that a large proportion of viral genes remain uncharacterized. The annotated viral genes with transcripts detected are shown in Fig. 3b with the majority having low transcript read depths (< 20 per sample). A range of viral structural genes (e.g., phage tail, head, capsid), genes encoding DNA/RNA polymerase, and genes related to DNA recombination or re-arrangement (resolvase, Rhs element vgr protein) were transcribed, indicative of active viruses in the soil incubations. After verifying the gene positions on the respective viral contigs, one auxiliary metabolic gene (AMG), acyl-CoA dehydrogenase, was also found to be transcribed (Fig. 3b and Supplementary Data 3). This gene encodes an enzyme that is involved in initial steps of fatty acid metabolism. Similarly, virus-encoded AMGs involved in the fatty acid oxidation were previously reported in viromes from the Pacific Ocean³⁷ and have been reported as a conserved strategy for a wide range of viruses in ocean³⁸. Our findings suggest that viruses may similarly be involved in host metabolism in soil.

In addition to transcripts, we detected several viral peptides in our soil metaproteome data that were indicative of viruses with a lytic lifestyle. To our knowledge this is the first report of using an untargeted, metaproteomic approach to detect viral proteins in soil, although it has previously been used with success for aquatic samples³⁹. For highly confident assignment of viral peptides with low abundances in the metaproteome, quality searching criteria (i.e., 5% false discovery rate (FDR), < 5 ppm mass error) were applied with manual inspection (Supplementary Data 4 and 6). Confident peptides for phage tail proteins, a virus-encoded protease, an RNA polymerase, and numerous chaperonins are captured in Fig. 3b. Inferring viral activity from metaproteomic data is supported by the fact that substantial components of viral structures (e.g., capsids and tails) are comprised of proteins³⁹. Due to the low abundances of viral proteins detected, we refrained from statistical comparisons of the impact of soil wetting and drying on the proteome data.

Interestingly, we detected a collection of chaperonin-like genes (K01802, K03554, K04043, and K04077) that were expressed by soil viruses at both the transcript and peptide level. Examples of confidently identified viral chaperonin peptide sequences (PSM < 1.19E-13, PPM < 2.83) are shown in Fig. 3d. These soil viral chaperonins shared similar sequence regions with bacterial, eukaryotic, and marine viral GroELs (Group I chaperonin), and also contained conserved features that clearly distinguished them as a novel group of soil viral chaperonins, as shown in the alignments (Supplementary Data 6). One example of a conserved region is displayed as an outer circle of the GroEL phylogenetic tree (Fig. 3c). Chaperonins have previously been expressed from a phage⁴⁰ and a plant virus (Closteroviruses)⁴¹. A recent study also reported a high prevalence of viral chaperonins including ones related to bacterial (GroELs) and archaeal (thermosome) chaperonins in an aquatic system⁴². Our phylogenetic analysis of bacterial, viral and eukaryotic GroEL genes (Fig. 3c) suggests that although chaperonins from the marine ecosystem are more phylogenetically similar to bacterial chaperonins, the majority of the soil viral chaperonins that we detected were most likely derived from eukaryotes, with some having bacterial origins. Such an apparent evolutionary separation of viral chaperonins in marine and soil ecosystems warrants further study. Because of the recognized functional importance of viral chaperonins in viral assembly and a lysogenic-lytic lifestyle transition⁴², the detection of soil viral GroEL at all levels, genomic, transcriptional, and translational, in our study suggests that the respective DNA viruses were actively infecting their hosts.

A diverse soil RNA virosphere

Our second aim was to examine the soil metatranscriptome data for RNA viruses. To date there have been few reports of RNA viruses in soil. Recently, RNA viruses in a marine study were found to be more abundant than DNA viruses implying important ecological functions in marine ecosystems⁷. The current knowledge of RNA viruses in soil is fragmentary, mainly focusing on culturable viruses and crop pathogens^43,44. Recently, Starr et al. (2019) reported detection of RNA viruses from metatranscriptomes of a California annual grassland soil³. Here we reassembled RNA viruses from the quality-filtered metatranscriptomic reads and uncovered a diverse RNA viral community in the Kansas native prairie soil.

The taxonomies of the identified soil RNA viral contigs are shown in rooted phylogenetic trees based on the alignments of the marker gene, RNA-dependent RNA polymerase (RdRP) for double-stranded RNA viruses, negative single-stranded RNA viruses, and positive single-stranded RNA viruses⁴⁵ (Fig. 4). Most transcripts from double-stranded RNA viruses were mapped to viral contigs annotated as Reoviridae, which are known to infect a wide variety of eukaryotic hosts⁴⁶ (Fig. 4a). One Reoviridae member that is closely related to a Bluetongue virus (MH559812, max. E-value of 1.86E-09, min. % identity of 88%) had the highest cumulative genomic coverage (Fig. 4a and Supplementary Data 7). This could be explained as either due to a bloom of the respective host during incubation, or due to their highly segmented genomes, i.e., ~10 dsRNA segments per genome, making them more resilient to sample processing steps⁴⁷. This is in contrast to Starr et al. (2019) who did not detect Reoviridae in the California grassland study³. The negative single-stranded RNA viruses had highest representation of Nairoviridae, Peribunyaviridae, and Hantaviridae across both wet and dry soils and Paramyxovirade in dry soil only (Fig. 4b). The positive single-stranded RNA viruses had highest representation of Secoviridae, Bromoviridae, Closteroviridae, Ifaviridae, Piconaviridae, Hypoviridae, and Leviviridae (Fig. 4c). The Picornavirales-like RNA viruses detected in Kansas soil were not previously found in California soil². More study is therefore needed to further extend our understanding of soil RNA viruses and their potential ecological functions.

**Fig. 4: Phylogenetic diversity of RNA viruses and their estimated abundances.**

More than half of the identified RNA viral sequences (124 out of 209) were assigned to Leviviridae, a family of negative-sense single-stranded RNA bacteriophages and the most diverse family found in Kansas soils (Fig. 4). The Leviviridae-like RNA viral sequences were clustered into two deeply branched major clades (Fig. 4d, e), supporting the recent hypothesis that the current Leviviridae family is evolving into two distinct lineages^45,48. Similar to our study, Starr et al. (2019) also found Leviviridae phages were the most diverse in California grassland soil. Leviviridae are only known to infect Proteobacteria⁴⁹ (Virus–Host DB, https://www.genome.jp/virushostdb/, accessed on 11 August 2020), one of the most abundant and diverse bacteria phyla in the Kansas soil (898 out of the total 4419 OTUs at 97% 16S rRNA gene sequence identity, Supplementary Fig. 1a). Similar observations of environmental viruses assigned to the most diverse and abundant taxa (e.g., Proteobacteria) have also been reported previously^1,2,4.

RNA viruses are responsive to soil moisture

Soil wetting and drying treatments impacted the RNA viral communities. Soil moisture treatments shaped the RNA viral communities from each of the Kansas soil locations (A and C) differently (Fig. 5a). The RNA viral communities generally grouped together according to sampling locations, except for one outlier from site A (Fig. 5a). The total RNA viral abundances were strongly correlated with the abundance of active eukaryotic species (based on 18S rRNA gene transcript counts from the metatranscriptomes, R² = 0.829 in Supplementary Fig. 1e). The correlation was higher in soils under wet conditions (R² = 0.856 in Supplementary Fig. 1f) compared to dry conditions (R² = 0.811 in Supplementary Fig. 1f). Among the detected RNA viruses, there were two viral families that responded significantly to the experimental soil wetting and drying. The estimated abundances of Leviviridae were higher in wet soils (p < 0.01, Fig. 5b). By contrast, a family of eukaryotic viruses (Paramyxoviridae) was more abundant in dry soils (p < 0.01, Fig. 5b). Due to the predominantly lytic lifestyle ascribed to Leviviridae³, higher abundances in wet soil compared to dry soil may reflect a higher degree of host lysis in wet soils. Proteobacteria, the host of Leviviridae⁵⁰, were also more abundant in wet soil (Supplementary Fig. 1a). Sampling across multiple time-points is needed to specifically investigate the continuous population dynamics of the virus–host communities, to better resolve interactions and transcriptional regulation at the community scale, and to elucidate infection dynamics.

**Fig. 5: Composition profiling of the RNA viral community and its abundance shift in response to wet and dry soil treatments.**

Conclusion

In summary, the highly diverse DNA and RNA viral composition, the increased viral activity in wet soil, and the potential virus–host interactions discovered here bring new knowledge about the ecological importance of the soil virosphere. As we demonstrate in this study, integration of information from multi-omics data is a promising approach for identifying both DNA and RNA viruses and assessing their potential activities in response to environmental changes. It has not escaped our notice during the global COVID-19 pandemic that the approaches described here can also be applied to detect novel pathogenic viral reservoirs in the environment. Because the soil we analyzed was a native prairie soil, with low levels of anthropogenic impacts, the majority of the active DNA viruses that we detected were bacteriophage or eukaryotic viruses known to have insect hosts. We found a high diversity of RNA viruses, including a group of relatively undefined eukaryotic viruses. Therefore, it would be interesting to compare this study to other soil environments, including those that are impacted by anthropogenic disturbances to further characterize the soil virosphere, and investigate their potential as reservoirs of viral pathogens.

Methods

Sampling site and soil treatment

Soil samples (top 0–15 cm) were collected in three locations at least 10 m apart from each of two sites (noted as sites A: 39°06'11ʺN, 96°36'48ʺW, 339 meters above sea level (MASL) and site C: 39°04'20ʺN, 96°34'33ʺW, 415 MASL in our previous study¹³) from native grassland locations at the Konza Experimental Field Station in Kansas. The soil samples from sites A and C were processed separately. For each site, the three soil samples were combined, homogenized, frozen and shipped in 1 gal Ziploc bags on dry ice to PNNL. Upon receipt at the laboratory, the samples were quickly processed through a 2 mm sieve. The soils from the two field locations (sites A and C) were separately aliquoted into ~ 50 g portions in Falcon tubes, flash frozen in liquid nitrogen and stored at −80 °C until use.

DNA was extracted from the original composite soil that was collected from the Konza field site for construction of two metagenomes as previously described¹³. One of the metagenomes was constructed for deep sequence coverage, with 1.1 Tb of sequence data¹⁴. The second metagenome used a hybrid assembly approach with long reads and short reads for improved assembly¹⁵.

Soil microcosms were prepared for each site (A and C) in three replicates and were either air-dried (dry treatment) or wet to saturation (wet treatment) and incubated at 21 °C for 15 days, as previously described¹³. DNA, RNA, and protein were extracted from each of the 12 microcosms: 2 field locations × 2 treatments × 3 replicates. The DNA was used for 16S rRNA gene sequencing as previously described¹³. RNA was extracted from 2 g of soil from each microcosm using the PowerSoil RNA extraction kit (MoBio Laboratories Inc., Carlsbad, CA, USA), according to the manufacturers’ instructions. Protein was extracted using the MPLex approach, as previously described¹¹.

Mining viral sequences from metagenome assemblies

DNA viral sequences (> 2500 bp) were mined from assembled contigs of the two existing metagenomes^14,15 described above. The de novo assembled contigs were screened for viral sequences by three independent methods including: (1) a probabilistic approach (taxonomy assignment based on multiple open reading frames, modified from Paez-Espino et al. (2017)); (2) public database searches (VirSorter⁵¹ and IMG/VR¹⁷); and (3) machine learning (VirFinder⁵²). The three methods were integrated and viral contigs were specified using four criteria: (1) contigs had more genes with viral origins than ‘bacterial or archaeal-like’ origins; (2) contigs were classified as viral contigs by VirSorter; (3) contigs had E-value cutoffs ≤1.0e⁻⁰⁵ after searching the IMG/VR database; and (4) contigs had p-values less than 0.05 and a VirFinder score greater than 0.90. Contigs that fulfilled at least three out of these four criteria were classified as viral contigs. Details of the different methods are provided below:

(1)
Probabilistic method

Protein sequences of the open reading frames were predicted and translated using Prodigal⁵³ and searched against viral, bacterial, and archaeal genome databases using hmmsearch (Hmmer v3.1b2⁵⁴, E-value cutoff of 1.0e⁻⁰⁵). The five databases searched included three viral databases: (1) a collection of 25,218 viral protein Hidden Markov Models (HMMs) built upon viral protein-coding genes from NCBI viral genomes⁵⁵; (2) a custom database comprised of 1147 curated viral protein families (Pfam⁵⁶, Supplementary Data 8); and (3) the existing Nucleo-Cytoplasmic Viruses Orthologous Groups (NCVOG)⁵⁷. In addition, the EggNOG bacterial and archaeal database⁵⁸ were used to annotate the proteins that are more likely from bacterial or archaeal genomes. The bit scores from searching against the five databases per protein were ranked and taxon annotation was assigned based on the search with the highest bit score. Modified from the previous study with arbitrary cutoffs⁵⁵, we retained the contigs with more ‘viral-like’ than ‘bacterial or archaeal-like’ genes on check and considered them as likely viral candidates, thus fulfilling our first criterion.
(2)
Public database searching

Similar to the ‘Probabilistic method’, we also screened for viral contigs using VirSorter⁵¹, one of the most widely used viral classification tools. VirSorter searches against profile HMMs and protein sets curated from public databases of NCBI RefSeqVirus (v69, January 2014) and virome sequences from selected habitats, which complement the databases searched in the ‘Probabilistic method’. The contigs classified as Category 1–3 by the VirSorter scoring matrix were retained as filling our second criterion above.

The contigs were also searched against IMG/VR (released on 1 July 2018), the largest publicly available database that includes both cultured and uncultured viral sequences¹⁷. Contigs having E-value cutoffs ≤1.0e⁻⁰⁵ fulfilled our third criterion and were retained.
(3)
Machine learning method

To access as many viral features as possible, contigs were also screened by VirFinder⁵², a machine learning-based method of modeling viral tetranucleotide patterns. Contigs with p-values less than 0.05 and a VirFinder score greater than 0.90 passed our fourth criterion.

After examination by the four criteria, the length-filtered assemblies (> 2500 bp) satisfying at least three out of the four criteria were identified as the viral contigs used for the following analyses.

DNA viral database curation

To expand the DNA viral reference dataset for transcript searching, the resulting viral contigs identified from the two metagenomes constructed from the source soil as described above, were combined with the two largest public viral databases, IMG/VR (released on 1 July 2018) and NCBI Virus¹⁶. The combined DNA viral sequences were then de-replicated to remove the exact and sub-strings of the sequences (VSEARCH v2.13.4, derep_fulllength) and curated into a larger DNA viral database including 726,108 unique DNA viral sequences. The newly curated DNA viral database was searched against the metatranscriptomes sequenced from wet and dry treatments for a comparison of active DNA viral communities.

Quality control of metatranscriptomic data

The existing metatranscriptomic data¹³ were screened both to determine which DNA viruses were transcriptionally active and to determine the RNA virosphere in wet and dry soil treatments. As described above, the metatranscriptomes were obtained from three replicate samples for each of the two independent locations after the wet and dry soil incubations, resulting in a total of 12 metatranscriptomes. The raw cDNA sequence reads were trimmed to remove the Illumina adapters and to retain high-quality reads (score > 30 and length > 36 bases), as recommended by Trimmomatic (v0.33). The metatranscriptomic reads were mapped to PhiX genome using Burrows-Wheeler Aligner (BWA, v0.7.17) and the exact matches were removed. The quality-filtered reads were then separated into reads corresponding to ribosomal and non-ribosomal RNA using SortMeRNA (version 2.1b). The ribosomal RNA reads were then filtered for 18S rRNA transcripts by mapping to a SILVA database including a set of well-curated 18S rRNA reference genes (SortMeRNA version 2.1b, silva-euk-18s-id95 database) to determine the transcriptionally active eukaryotic community composition (more details in Supplementary Fig. 1b). The remaining reads, corresponding to non-ribosomal RNA, were used to determine which DNA viruses were transcriptionally active and also to assemble RNA viruses.

Transcriptional activity and diversity of DNA viruses

The non-ribosomal RNA reads were mapped to our self-curated DNA viral database using BamM (v1.7.3, bamm make, https://github.com/Ecogenomics/BamM) and the mapping was filtered using stringent cutoffs, i.e., % identity > 0.95 and % alignment > 0.80 (BamM v1.7.3, bamm filter). DNA viral sequences that were mapped by the quality-filtered metatranscriptomic reads were considered as transcribed DNA viruses. Transcriptional activity was estimated by the average base coverage of the viral sequences (samtools v1.9, samtools depth, http://www.htslib.org/doc/) normalized by the total counts of reads per sample and the length of the searched DNA viral contigs.

In addition, transcripts of lysogenic markers, integrase and excisionase, were recruited from the metatranscriptomes. A total of 20,712 unique viral integrases and 329 unique excisionases were collected from NCBI Virus (https://www.ncbi.nlm.nih.gov/labs/virus/vssi/#/, accessed on 16 December 2020; sequence de-replication using VSEARCH v2.13.4, derep_fulllength). The same transcript read mapping strategies and quality filtering criteria as mentioned above were applied. The transcript counts of the lysogenic markers were normalized by the total number of reads per sample.

The diversity of the transcriptionally active DNA viruses was evaluated by clustering the viral sequences based on their average amino acid identities (AAIs) using the Ruby script in enveomics toolbox⁵⁹ and converting the output into a distance matrix. The single linkage clustering in scikit-learn (v0.23.1) was applied to the matrix to cluster the viral sequences at the cutoff of 70% AAI. Because the viral genome assemblies from the soil metagenomes were incomplete, fragments from different regions of the same viral genome were not clustered together based on AAI, when no significant amounts of proteins/genes were shared. We also incorporated a complementary method for clustering the viral contigs according to their tetranucleotide patterns. The tetranucleotide frequency-based method was previously shown to make linkages between nonoverlapping but related viral contigs⁶⁰ and is the main principle underlying bioinformatics tools like VirFinder⁵² and PhaMers⁶¹. Pair-wise Pearson correlations were calculated between z-score distributions for every tetranucleotide in the screened viral contigs (pyani v0.2.9, https://github.com/widdowquinn/pyani). The viral contigs that were not clustered by the AAI method were assigned to the clusters containing the paired contigs with the highest correlation coefficients of the tetranucleotide frequency z-scores.

Transcription of coding and noncoding regions of DNA viruses

The quality-filtered transcripts that mapped to the DNA viral contigs as mentioned above (% identity > 0.95 and % alignment > 0.80) were further separated into those that mapped to gene-coding and noncoding regions of the DNA viral contigs (predicted by Prodigal, v2.6.3) and counted individually. Potential AMGs were screened from the transcribed genes of the DNA viral contigs. The most confident AMG candidates were recognized as those that were annotated with auxiliary metabolic functions (E-value cutoff of 1 × 10⁻⁵) and located on viral contigs with ‘viral-like’ genes both up and downstream.

Assembly of RNA viruses from metatranscriptome data

The quality-filtered metatranscriptomic reads were also de novo assembled by Trinity (v2.8.5). RNA viruses were screened from the de novo assemblies in two ways: (1) searching for the marker gene, RdRP, using a set of curated HMMs, and (2) the whole-sequence similarity searches against NCBI RNA viral genomes¹⁶ (E-value cutoff of 1 × 10⁻⁵). Taxonomy annotation of the identified RNA viral contigs was based on RdRP gene similarity, which has been used as a phylogenetic marker for RNA viruses³. Due to the incomplete assembly of the viral genomes, some identified RNA viral contigs were not detected by RdRP gene-screening using HMMs, but were instead annotated by their significant hits in NCBI RNA viral database (E-value cutoff of 1 × 10⁻⁵). Multiple Alignment with the Fast Fourier Transform (MAFFT v7) program was used to align the RdRP genes and a phylogenetic tree was constructed using the maximum-likelihood method via FastTree (v2.1.10) to display the RNA virosphere. The quality-filtered metatranscriptomic reads were mapped to the detected RNA viral assemblies and filtered with stringent criteria (95% identity, 80% coverage) by BamM (v1.7.3, bamm make and bamm filter, https://github.com/Ecogenomics/BamM). The abundances of the RNA contigs were estimated by the average base coverage of the mined RNA viral sequences after filtering the mapping at % identity > 0.95 and % alignment > 0.80 (samtools v1.9, samtools depth, http://www.htslib.org/doc/). The relative abundance of each taxon was measured by taking the average of the estimated coverage of RNA viral contigs assigned to the taxon. The genome coverage of each RNA viral taxon was used to represent the composition of the detected RNA viral community under each treatment.

Mining viral peptides from MS/MS spectra

Metaproteomes were obtained from the same soil samples as for the metatranscriptomes, resulting in 12 matched metaproteomes (three replicates each for wet and dry incubations and two locations). The LC-MS/MS were performed using two different approaches as described previously⁸. Briefly, all of the samples were analyzed individually with one-dimensional LC-MS/MS. In addition, the pooled replicate samples went through off-line liquid-chromatography-based fractionation followed by LC-MS/MS analysis⁸.

We first curated a viral protein database including the proteins predicted from contigs in DNA viral database mentioned above (Prodigal, v2.6.3), plus a set of unique reference viral structural proteins collected from NCBI Virus (a total of 63,088 capsid proteins, 25 envelope and 15,817 tail proteins). The LC-MS/MS spectra were searched against the metaproteome database using MS-GF+ search engine⁶². The spectrum level peptide confidence score of the peptide-spectrum match (PSM; i.e., MSGFDB_SpecProb in MS-GF+) and mass error (in ppm) of the precursor peptide ion (i.e., DelM_PPM in MS-GF+) were optimized to achieve the highest number of peptide identification within each dataset while keeping the target-decoy-based FDR of peptide identification below 5%. In order to obtain the most confident viral identification, any redundant peptide identifications between the viral metaproteome and the rest of the metaproteome were excluded. The remaining peptide identification were from the viral metaproteome only. The data also went through a manual quality control process to remove any PSM with low confidence based on MS/MS fragmentation coverage of the peptide and MS peak quality. Final spectra are shown in Supplementary Data 5, and their qualification are shown in Supplementary Data 4.

Host prediction of identified DNA viruses

An integrated approach was applied to assign host taxa to the DNA viral contigs. DNA viral contigs were first searched for exact matches in the spacers detected from the clustered regulatory interspaced short palindromic repeat (CRISPR) regions on the non-viral contigs from the two metagenomes (Method 1). Spacers were detected using MinCED (https://github.com/ctSkennerton/minced) and searched against the viral contigs using optimized parameters as previously described^63,64. The viral contigs were assigned to hosts based on the taxonomy annotations of the respective non-viral contigs extracted with the spacer hits using the Contig Annotation Tool (CAT, v4.6, https://github.com/dutilh/CAT). Sequence similarity searches to non-viral contigs (potential hosts) and NCBI reference viruses were then applied to the remaining viral contigs that were not assigned by the spacer hits. A non-viral contig was assigned to a putative viral host when they shared a high sequence similarity (Method 2, BLASTN, E-value threshold = 10⁻³ and bit score = 50). Alternatively, when a viral contig had a close hit in the NCBI viral genome database (BLASTN, E-value threshold = 10⁻³ and bit score = 50), we adopted the host information of the closest reference virus (Method 3). In addition, we trained homogeneous Markov models (Method 4) using 1853 NCBI RefSeq genomes and applied a probabilistic approach coded in WIsh⁶⁵, known for accurate host predictions of short viral contigs. We only retained the virus–host pairings that were detected using our nested workflow (Methods 1-2-3) and WIsh⁶⁵ (Method 4) for highly confident host assignments.

Phylogenetic placement of soil viral chaperonin

First, 4821 bacterial and eukaryotic Group I chaperonin (GroEL) protein sequences were extracted from the NCBI Reference Sequences (https://www.ncbi.nlm.nih.gov/refseq/, accessed July 2020) and de-replicated using VSEARCH (v2.13.4). A total of 235 bacterial and eukaryotic GroEL protein sequences from unique Orders were selected to represent the phylogenetic diversity. An additional 108 viral chaperonins from an intensive study of marine viruses possessing a bacterial-derived GroEL⁴² were added. Protein sequences predicted from the Kansas viral contigs that were annotated as chaperonins and detected in either metatranscriptome or metaproteome were aligned together with the bacterial, eukaryotic and marine viral GroELs using MAFFT (v7). The alignment with the estimated sequence conservation was visualized in Jalview (v2.11.1.0). A phylogenetic tree of GroEL across kingdoms was constructed using maximum-likelihood method via FastTree (v2.1.10).

Statistics and reproducibility

t-Tests were conducted to assess the significance of the different counts of viral sequences and transcripts mapped under wet and dry treatments. The DESeq2 package was used for differential gene expression analysis. Differential abundances with p-values that were less than 0.05 were considered as significant. Sample size and number of replicates are described in ‘Sampling site and soil treatment’ and provided in the results and figure legends. The experiment was run using replicate samples for each treatment and at least reproducible results were found when comparing between replicates.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

The two published metagenomes can be downloaded from the PNNL Data Hub (https://data.pnnl.gov/project/12620) and White et al.¹⁵. The metatranscriptome sequence data and the 16S rRNA sequence data are available in Chowdhury et al.¹³. The MS proteomics data have been deposited to the ProteomeXchange Consortium via the MASSIVE partner repository with the accession number of MSV000086144 (https://massive.ucsd.edu/ProteoSAFe/dataset.jsp?task=500afc8e54e1481a9d90fcd35b556511). The DNA and RNA viral sequences identified from de novo assemblies are deposited under NCBI BioProject accession number PRJNA744880. All of the curated databases, including the reference viral structural proteins, viral lysogenic gene markers and the curated DNA viral database, and source data that were used to plot the figures are packaged at the public repository Zenodo⁶⁶.

Code availability

The R codes used to plot the main and supplementary figures are available at Zenodo with no restriction to access⁶⁶.

References

Paez-Espino, D. et al. Uncovering Earth’s virome. Nature 536, 425 (2016).
Article CAS PubMed Google Scholar
Emerson, J. B. et al. Host-linked soil viral ecology along a permafrost thaw gradient. Nat. Microbiol. 3, 870 (2018).
Article CAS PubMed PubMed Central Google Scholar
Starr, E. P., Nuccio, E. E., Pett-Ridge, J., Banfield, J. F. & Firestone, M. K. Metatranscriptomic reconstruction reveals RNA viruses with the potential to shape carbon cycling in soil. Proc. Natl Acad. Sci. USA 116, 25900–25908 (2019).
Article CAS PubMed PubMed Central Google Scholar
Trubl, G. et al. Soil viruses are underexplored players in ecosystem carbon processing. mSystems 3, e00076–00018 (2018).
Article CAS PubMed PubMed Central Google Scholar
Schulz, F. et al. Hidden diversity of soil giant viruses. Nat. Commun. 9, 4881 (2018).
Article PubMed PubMed Central CAS Google Scholar
Dias, R. S. et al. Viral concentration methods for diversity studies in soil samples. Appl. Soil Ecol. 155, 103666 (2020).
Article Google Scholar
Steward, G. F. et al. Are we missing half of the viruses in the ocean? ISME J. 7, 672–679 (2013).
Article CAS PubMed Google Scholar
Callister, S. J. et al. Addressing the challenge of soil metaproteome complexity by improving metaproteome depth of coverage through two-dimensional liquid chromatography. Soil Biol. Biochem. 125, 290–299 (2018).
Article CAS Google Scholar
Hultman, J. et al. Multi-omics of permafrost, active layer and thermokarst bog soil microbiomes. Nature 521, 208–212 (2015).
Article CAS PubMed Google Scholar
Woodcroft, B. J. et al. Genome-centric view of carbon processing in thawing permafrost. Nature 560, 49–54 (2018).
Article CAS PubMed Google Scholar
Nicora, C. D. et al. The MPLEx protocol for multi-omic analyses of soil samples. J. Vis. Exp. 30, e57343 (2018).
Burnum-Johnson, K. E. et al. Ion mobility spectrometry and the omics: distinguishing isomers, molecular classes and contaminant ions in complex samples. Trends Anal. Chem. 116, 292–299 (2019).
Article CAS Google Scholar
Chowdhury, T. R. et al. Metaphenomic responses of a native prairie soil microbiome to moisture perturbations. mSystems 4, 00061-19 (2019).
Nelson, W. C. et al. Terabase metagenome sequencing of grassland soil microbiomes. Microbiol. Resour. Announc. 9, e00718-20. (2020).
White, R. A. et al. Moleculo long-read sequencing facilitates assembly and genomic binning from complex soil metagenomes. mSystems 1, e00045-16 (2016).
Brister, J. R., Ako-Adjei, D., Bao, Y. & Blinkova, O. NCBI viral genomes resource. Nucleic Acids Res. 43, D571–D577 (2015).
Article CAS PubMed Google Scholar
Paez-Espino, D. et al. IMG/VR v. 2.0: an integrated data management and analysis system for cultivated and environmental viral genomes. Nucleic Acids Res. 47, D678–D686 (2018).
Article PubMed Central CAS Google Scholar
Krupovic, M. et al. Cressdnaviricota: a virus phylum unifying 7 families of Rep-encoding viruses with single-stranded, circular DNA genomes. J. Virol. 94, e00582-20 (2020).
Aylward, F. O. et al. Diel cycling and long-term persistence of viruses in the ocean’s euphotic zone. Proc. Natl Acad. Sci. USA 114, 11446–11451 (2017).
Article CAS PubMed PubMed Central Google Scholar
Allen, L. Z. et al. The Baltic Sea virome: diversity and transcriptional activity of DNA and RNA viruses. mSystems 2, e00125-16 (2017).
Dupont, C. L. et al. Genomes and gene expression across light and productivity gradients in eastern subtropical Pacific microbial communities. ISME J. 9, 1076–1092 (2015).
Article CAS PubMed Google Scholar
Sieradzki, E. T., Ignacio-Espinoza, J. C., Needham, D. M., Fichot, E. B. & Fuhrman, J. A. Dynamic marine viral infections and major contribution to photosynthetic processes shown by spatiotemporal picoplankton metatranscriptomes. Nat. Commun. 10, 1–9 (2019).
Article CAS Google Scholar
Malki, K. et al. Bacteriophages isolated from Lake Michigan demonstrate broad host-range across several bacterial phyla. Virol. J. 12, 164 (2015).
Article PubMed PubMed Central CAS Google Scholar
Puschnik, A. S., Majzoub, K., Ooi, Y. S. & Carette, J. E. A CRISPR toolbox to study virus–host interactions. Nat. Rev. Microbiol. 15, 351–364 (2017).
Article CAS PubMed PubMed Central Google Scholar
Neufing, P. J., Shearwin, K. E. & Egan, J. B. Establishing lysogenic transcription in the temperate coliphage 186. J. Bacteriol. 183, 2376–2379 (2001).
Article CAS PubMed PubMed Central Google Scholar
Marsh, P. & Wellington, E. Phage–host interactions in soil. FEMS Microbiol. Ecol. 15, 99–107 (1994).
Article CAS Google Scholar
Stewart, F. M. & Levin, B. R. The population biology of bacterial viruses: why be temperate. Theor. Popul. Biol. 26, 93–117 (1984).
Article CAS PubMed Google Scholar
Williamson, K. E., Radosevich, M., Smith, D. W. & Wommack, K. E. Incidence of lysogeny within temperate and extreme soil environments. Environ. Microbiol. 9, 2563–2574 (2007).
Article CAS PubMed Google Scholar
Hug, L. A. et al. Aquifer environment selects for microbial species cohorts in sediment and groundwater. ISME J. 9, 1846–1856 (2015).
Article PubMed PubMed Central CAS Google Scholar
Gravuer, K. & Eskelinen, A. Nutrient and rainfall additions shift phylogenetically estimated traits of soil microbial communities. Front. Microbiol. 8, 1271 (2017).
Article PubMed PubMed Central Google Scholar
Yang, H. et al. Transcription regulation mechanisms of bacteriophages: recent advances and future prospects. Bioengineered 5, 300–304 (2014).
Article PubMed PubMed Central Google Scholar
Van Goethem, M. W., Swenson, T. L., Trubl, G., Roux, S. & Northen, T. R. Characteristics of wetting-induced bacteriophage blooms in biological soil crust. mBio 10, e02287-19 (2019).
Bloch, S., Lewandowska, N., Węgrzyn, G. & Nejman-Faleńczyk, B. Bacteriophages as sources of small non-coding RNA molecules. Plasmid 113, 102527 (2020).
Tycowski, K. T. et al. Viral noncoding RNAs: more surprises. Genes Dev. 29, 567–584 (2015).
Article CAS PubMed PubMed Central Google Scholar
Owen, S. V. et al. A window into lysogeny: revealing temperate phage biology with transcriptomics. Microb. Genom. 6, e000330 (2020).
Ranade, K. & Poteete, A. R. A switch in translation mediated by an antisense RNA. Genes Dev. 7, 1498–1507 (1993).
Article CAS PubMed Google Scholar
Hurwitz, B. L., Hallam, S. J. & Sullivan, M. B. Metabolic reprogramming by viruses in the sunlit and dark ocean. Genome Biol. 14, R123 (2013).
Article PubMed PubMed Central Google Scholar
Rosenwasser, S., Ziv, C., Van Creveld, S. G. & Vardi, A. Virocell metabolism: metabolic innovations during host–virus interactions in the ocean. Trends Microbiol. 24, 821–832 (2016).
Article CAS PubMed Google Scholar
Brum, J. R. et al. Illuminating structural proteins in viral “dark matter” with metaproteomics. Proc. Natl Acad. Sci. USA 113, 2436–2441 (2016).
Article CAS PubMed PubMed Central Google Scholar
Kurochkina, L. P. et al. Expression and functional characterization of the first bacteriophage-encoded chaperonin. J. Virol. 86, 10103–10111 (2012).
Article CAS PubMed PubMed Central Google Scholar
Satyanarayana, T. et al. Closterovirus encoded HSP70 homolog and p61 in addition to both coat proteins function in efficient virion assembly. Virology 278, 253–265 (2000).
Article CAS PubMed Google Scholar
Marine, R. L., Nasko, D. J., Wray, J., Polson, S. W. & Wommack, K. E. Novel chaperonins are prevalent in the virioplankton and demonstrate links to viral biology and ecology. ISME J. 11, 2479–2491 (2017).
Article PubMed PubMed Central Google Scholar
Andika, I. B., Kondo, H. & Sun, L. Interplays between soil-borne plant viruses and RNA silencing-mediated antiviral defense in roots. Front. Microbiol. 7, 1458 (2016).
Article PubMed PubMed Central Google Scholar
Zhang, R. et al. A novel single-stranded RNA virus isolated from a phytopathogenic filamentous fungus, Rosellinia necatrix, with similarity to hypo-like viruses. Front. Microbiol. 5, 360 (2014).
Article PubMed PubMed Central Google Scholar
Wolf, Y. I. et al. Origins and evolution of the global RNA virome. mBio 9, e02329-18 (2018).
Hillman, B. I., Supyani, S., Kondo, H. & Suzuki, N. A reovirus of the fungus Cryphonectria parasitica that is infectious as particles and related to the Coltivirus genus of animal pathogens. J. Virol. 78, 892–898 (2004).
Article CAS PubMed PubMed Central Google Scholar
McDonald, S. M., Nelson, M. I., Turner, P. E. & Patton, J. T. Reassortment in segmented RNA viruses: mechanisms and outcomes. Nat. Rev. Microbiol. 14, 448 (2016).
Article CAS PubMed PubMed Central Google Scholar
Callanan, J. et al. Expansion of known ssRNA phage genomes: from tens to over a thousand. Sci. Adv. 6, eaay5981 (2020).
Article CAS PubMed PubMed Central Google Scholar
Bollback, J. P. & Huelsenbeck, J. P. Phylogeny, genome evolution, and host specificity of single-stranded RNA bacteriophage (family Leviviridae). J. Mol. evolution 52, 117–128 (2001).
Article CAS Google Scholar
Krishnamurthy, S. R., Janowski, A. B., Zhao, G., Barouch, D. & Wang, D. Hyperexpansion of RNA bacteriophage diversity. PLoS Biol. 14, e1002409 (2016).
Article PubMed PubMed Central CAS Google Scholar
Roux, S., Enault, F., Hurwitz, B. L. & Sullivan, M. B. VirSorter: mining viral signal from microbial genomic data. PeerJ 3, e985 (2015).
Ren, J., Ahlgren, N. A., Lu, Y. Y., Fuhrman, J. A. & Sun, F. VirFinder: a novel k-mer based tool for identifying viral sequences from assembled metagenomic data. Microbiome 5, 69 (2017).
Article PubMed PubMed Central Google Scholar
Hyatt, D. et al. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11, 119 (2010).
Eddy, S. R. Accelerated profile HMM searches. PLoS Comput. Biol. 7, e1002195 (2011).
Article CAS PubMed PubMed Central Google Scholar
Paez-Espino, D., Pavlopoulos, G. A., Ivanova, N. N. & Kyrpides, N. C. Nontargeted virus sequence discovery pipeline and virus clustering for metagenomic data. Nat. Protoc.12, 1673–1682 (2017).
Finn, R. D. et al. Pfam: the protein families database. Nucleic Acids Res. 42, D222–D230 (2013).
Article PubMed PubMed Central CAS Google Scholar
Yutin, N., Wolf, Y. I., Raoult, D. & Koonin, E. V. Eukaryotic large nucleo-cytoplasmic DNA viruses: clusters of orthologous genes and reconstruction of viral genome evolution. Virol. J. 6, 223 (2009).
Huerta-Cepas, J. et al. eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucleic Acids Res. 44, D286–D293 (2016).
Rodriguez-R, L. M. & Konstantinidis, K. T. The enveomics collection: a toolbox for specialized analyses of microbial genomes and metagenomes. PeerJ Preprints 4, e1900v1; Report No. 2167-9843, https://doi.org/10.7287/peerj.preprints.1900v1 (2016).
Trifonov, V. & Rabadan, R. Frequency analysis techniques for identification of viral genetic data. mBio 1, e00156–e00110 (2010).
Article PubMed PubMed Central CAS Google Scholar
Deaton, J., Yu, F. B. & Quake, S. R. Mini‐Metagenomics and Nucleotide Composition Aid the Identification and Host Association of Novel Bacteriophage Sequences. Advanced biosystems. 3, 1900108 (2019).
Kim, S. & Pevzner, P. A. MS-GF+ makes progress towards a universal database search tool for proteomics. Nat. Commun. 5, 5277 (2014).
Article CAS PubMed Google Scholar
Biswas, A., Gagnon, J. N., Brouns, S. J., Fineran, P. C. & Brown, C. M. CRISPRTarget: bioinformatic prediction and analysis of crRNA targets. RNA Biol. 10, 817–827 (2013).
Article CAS PubMed PubMed Central Google Scholar
Edwards, R. A., McNair, K., Faust, K., Raes, J. & Dutilh, B. E. Computational approaches to predict bacteriophage–host relationships. FEMS Microbiol. Rev. 40, 258–272 (2015).
Article PubMed PubMed Central CAS Google Scholar
Galiez, C., Siebert, M., Enault, F., Vincent, J. & Söding, J. WIsH: who is the host? Predicting prokaryotic hosts from metagenomic phage contigs. Bioinformatics 33, 3113–3114 (2017).
Article CAS PubMed PubMed Central Google Scholar
Wu, R. et al. Moisture modulates soil reservoirs of active DNA and RNA viruses. zenodo https://doi.org/10.5281/zenodo.5081148 (2021).

Download references

Acknowledgements

This work was supported by the U.S. Department of Energy (DOE), Office of Biological and Environmental Research (BER) and is a contribution of the Scientific Focus Area ‘Phenotypic response of the soil microbiome to environmental perturbations’. Pacific Northwest National Laboratory (PNNL) is operated for the DOE by Battelle Memorial Institute under Contract DE-AC05-76RL01830. Portions of this work was performed in the W. R. Wiley Environmental Molecular Sciences Laboratory (EMSL) (grid.436923.9), a DOE Office of Science User Facility sponsored by the Office of BER and located at PNNL. The metaproteomic work was funded by the U.S. DOE, Office of Science, Office of BER, Early Career Research Program to K.E.B.J. We also thank Joon Yong Lee for his valuable insights into analysis of the metaproteomic data.

Author information

Authors and Affiliations

Earth and Biological Sciences Directorate, Pacific Northwest National Laboratory, Richland, WA, USA
Ruonan Wu, Michelle R. Davison, Yuqian Gao, Carrie D. Nicora, Jason E. Mcdermott, Kristin E. Burnum-Johnson, Kirsten S. Hofmockel & Janet K. Jansson
Department of Agronomy, Iowa State University, Ames, IA, USA
Kirsten S. Hofmockel

Authors

Ruonan Wu
View author publications
You can also search for this author in PubMed Google Scholar
Michelle R. Davison
View author publications
You can also search for this author in PubMed Google Scholar
Yuqian Gao
View author publications
You can also search for this author in PubMed Google Scholar
Carrie D. Nicora
View author publications
You can also search for this author in PubMed Google Scholar
Jason E. Mcdermott
View author publications
You can also search for this author in PubMed Google Scholar
Kristin E. Burnum-Johnson
View author publications
You can also search for this author in PubMed Google Scholar
Kirsten S. Hofmockel
View author publications
You can also search for this author in PubMed Google Scholar
Janet K. Jansson
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

J.K.J. and K.S.H. obtained funding for the study. J.K.J. and R.W. designed the experiment. R.W. performed the majority of the analyses. Y.G and K.E.B.J. performed peptide prediction from MS/MS spectra data and conducted similarity search to viral proteins. C.D.N. extracted soil proteins and prepared samples for MS/MS spectra detection. J.S.M. and M.R.D. provided expertise in virology and helped to interpret the data. All authors contributed to writing of the manuscript.

Corresponding author

Correspondence to Janet K. Jansson.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Communications Biology thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editor: Caitlin Karniski.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary information.

Description of Supplementary Files.

Supplementary Data 1. Summary of the 416 transcribed DNA viral contigs.

Supplementary Data 2. Transcript reads of lysogenic marker genes encoding viral integrase and excisionase.

Supplementary Data 3. Annotations and normalized transcript counts of DNA viral genes.

Supplementary Data 4. Viral peptides detected in metaproteomes.

Supplementary Data 5. Raw spectrum graphs of the 60 highly confident mass spectrometry spectral hits.

Supplementary Data 6. The alignments of bacterial, eukaryotic, marine and soil viral chaperonin protein sequences.

Supplementary Data 7 RNA viral contigs identified from the Kansas grassland soil.

Supplementary Data 8. Viral Pfam list for viral gene annotation.

Reporting summary.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Wu, R., Davison, M.R., Gao, Y. et al. Moisture modulates soil reservoirs of active DNA and RNA viruses. Commun Biol 4, 992 (2021). https://doi.org/10.1038/s42003-021-02514-2

Download citation

Received: 21 October 2020
Accepted: 18 July 2021
Published: 26 August 2021
DOI: https://doi.org/10.1038/s42003-021-02514-2

This article is cited by

Virus diversity and activity is driven by snowmelt and host dynamics in a high-altitude watershed soil ecosystem
- Clement Coclet
- Patrick O. Sorensen
- Simon Roux
Microbiome (2023)
Propagation of viral genomes by replicating ammonia-oxidising archaea during soil nitrification
- Sungeun Lee
- Ella T Sieradzki
- Christina Hazard
The ISME Journal (2023)
Extreme summers impact cropland and grassland soil microbiomes
- Qicheng Bei
- Thomas Reitz
- Anna Heintz-Buschart
The ISME Journal (2023)
Soil viral diversity, ecology and climate change
- Janet K. Jansson
- Ruonan Wu
Nature Reviews Microbiology (2023)
Geochemical constraints on bacteriophage infectivity in terrestrial environments
- Hans K Carlson
- Denish Piya
- Vivek K Mutalik
ISME Communications (2023)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.