Next Article in Journal
Understanding the Ecology of Restored Fen Peatlands for Protection and Sustainable Use
Previous Article in Journal
Formation and Fluxes of Soil Trace Gases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Coming-of-Age Characterization of Soil Viruses: A User’s Guide to Virus Isolation, Detection within Metagenomes, and Viromics

1
Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA
2
Department of Biology/Toxicology, Ashland University, Ashland, OH 44805, USA
3
DOE Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
4
Department of Microbiology, The Ohio State University, Mansfield, OH 44906, USA
*
Authors to whom correspondence should be addressed.
Submission received: 25 February 2020 / Revised: 4 April 2020 / Accepted: 17 April 2020 / Published: 21 April 2020

Abstract

:
The study of soil viruses, though not new, has languished relative to the study of marine viruses. This is particularly due to challenges associated with separating virions from harboring soils. Generally, three approaches to analyzing soil viruses have been employed: (1) Isolation, to characterize virus genotypes and phenotypes, the primary method used prior to the start of the 21st century. (2) Metagenomics, which has revealed a vast diversity of viruses while also allowing insights into viral community ecology, although with limitations due to DNA from cellular organisms obscuring viral DNA. (3) Viromics (targeted metagenomics of virus-like-particles), which has provided a more focused development of ‘virus-sequence-to-ecology’ pipelines, a result of separation of presumptive virions from cellular organisms prior to DNA extraction. This separation permits greater sequencing emphasis on virus DNA and thereby more targeted molecular and ecological characterization of viruses. Employing viromics to characterize soil systems presents new challenges, however. Ones that only recently are being addressed. Here we provide a guide to implementing these three approaches to studying environmental viruses, highlighting benefits, difficulties, and potential contamination, all toward fostering greater focus on viruses in the study of soil ecology.

Graphical Abstract

1. Introduction

Viruses are acellular infectious agents that can exist extracellularly as nucleic-acid genomes encapsidated within proteinaceous structures. These infectious agents can be found in environments wherever cellular organisms are present, particularly cellular microorganisms (microbes), such as bacteria, archaea, fungi, and protists [1]. There are an estimated 1031 viruses on Earth with a majority of them infecting microbes, making viruses significant drivers of evolution and essential for life on Earth [2,3]. Soils, as we emphasize here, are among the most virus-rich environments [4,5,6,7].
Soils represent the primary terrestrial habitat of microbes, but research on soil viruses has lagged behind research of viruses in aquatic environments. Our understanding of the soil virosphere consequently consists mostly, at best, of a ‘black box’, particularly in terms of the contributions of viruses to soil ecology. This underdevelopment of soil virus ecology stems mainly from complications in simply gaining access to soil viruses and their nucleic acid.
Though discouraging many from even attempting to analyze the soil virosphere, still there are multiple compelling reasons to do so (Section 2). Here, we present a “user’s guide” to assessing soil viruses, one aimed at better illuminating the soil virus “black box” by describing three prominent approaches, their pros and cons, and providing suggestions on how to reduce methodological difficulties (Section 3 and Section 4). Our objective is to help researchers to better recognize what viruses may be present in a given soil sample, and to do so particularly toward assessment of the ecological impacts of viruses on different soil ecosystems.

Virus-Like Particles (VLPs), Viruses, Microbes, and Other Terms

Virions consist of infectious, encapsidated nucleic acid. As is customary in the field of virus ecology, however, the term ‘virus-like particle’ (VLP) is often used instead of ‘virion’. VLPs specifically are entities of virion size that contain nucleic acid but have not otherwise been identified as viruses (described more in Section 4.1, Section 4.2 and Section 4.3).
In soils, viruses that infect bacteria, also known as bacteriophages (phages), are the most abundant and most studied. The term ‘virus’, however, is used rather than ‘phage’ when a virus’ host has not been identified or when referring to naturally occurring viral communities, even if the most abundant viruses are phages. The term ‘phage’, that is, should not be used to describe the viruses of either domain Archaea or domain Eukarya [8].
To describe microbial cellular organisms, we use the term ‘microbe’. Though a microbe is likely to be a bacterium due to their typically greater abundance, archaea and eukaryotes are also found among soil ‘microbes’. This review thus predominantly covers approaches to determining what viruses of microbes [1,9] are present in soils, keeping in mind that most of those viruses are phages and most of those microbes are bacteria. A glossary of additional terms used in this review is found in the Appendix A.

2. Potential Roles and Forms of Viruses in Soils

The roles that viruses play in soils may be ecologically equivalent to the roles of viruses in oceans. There viral ecogenomics—characterizing the ecology of viruses from their genomes—and the roles of viruses in microbial biogeochemistry have been investigated to a much greater extent [10], albeit with a primary focus on double-stranded DNA (dsDNA) viruses [11]. In this section, we highlight three virus–host interactions that are potentially translatable between marine and soil environments. These include killing and lysing cells (Section 2.2), alteration of host metabolism during infection (Section 2.3), and virus-mediated horizontal gene transfer (i.e., transduction) (Section 2.4). We begin, though, with a more philosophical consideration of the utility of including viruses in considerations of soil ecology (Section 2.1), and then conclude this section with a discussion of the different ways that viruses can exist in soils (Section 2.5).

2.1. Importance of the Soil Virosphere

Why should we care about viruses in soils? For some, the discovery aspect alone is enough to motivate soil virus research, while others want to know to what extent viruses actually matter (e.g., in terms of their impact on soil biogeochemistry). It could be argued, however, that focus should primarily be on characterizing the metabolisms of soil organisms, whereas viruses, at least arguably, do not possess metabolisms. Rather, viruses hijack the already existing metabolisms of their hosts, especially to produce more virions, or at least to produce more virus genomes [12].
This reliance on the metabolisms of their hosts is one of the reasons why viruses technically are often not considered to be “alive” [13,14,15]. One result of these sometimes semantic arguments can be a focus on hosts rather than on the ‘distraction’ of less obviously metabolically contributing viruses. The host-centric dogma is perhaps best summed up by Fierer [16], who suggested that viruses, particularly in their role as predators of cellular organisms, represent the least important factor influencing the composition of soil microbial communities.
Contributing to this viral de-emphasis, little is known about the impacts of viruses in soils, including in terms of their roles as predators of soil microbes. Thus, the following questions should be considered: (Section 2.2) To what extent do soil viruses impact soil microbe numbers and diversity? (Section 2.3) To what degree do viruses directly modify the metabolisms of soil microbes? (Section 2.4) How extensively do viruses horizontally transfer genes between soil microbes? Answering these questions will allow for a better understanding of the importance of the soil virosphere to soil functioning.

2.2. Viral Lysis of Soil Microbes

Viruses in soil ecosystems can kill, that is, serve as ‘predators’ of microbes. Predators, though, are traditionally considered to be organisms that kill and then consume other organisms (prey), assimilating the nutrients associated with the consumed organisms into their own bodies. Predators also tend to consume multiple prey over their lifetimes. Viruses, by contrast, exploit only a single, still metabolizing host per viral generation, and generally consume relatively little of the prey organism, though they can still assimilate 20%–30% of prey mass into new virion particles [17,18,19,20]. This is true even for viruses that are strictly lytic [21], i.e., those which can successfully replicate only in conjunction with host-cell killing. A better metaphor for lytic viruses might be that of parasitoids [22], such as larval wasps, which lethally consume their hosts alive, from within, before emerging in a mature form. Viruses nevertheless still can function across ecosystems in a predator-like manner by driving Lotka–Volterra-like predator–prey dynamics [23,24].
Is the potential of viruses to kill microbes their only aspect that matters in soil? In ocean water columns, viruses are thought to lyse about one-third of microbes each day, and in the process they collectively release about 10 billion tons of nutritious ‘necromass’ into the extracellular environment [25,26]. This nutrient infusion can result in an increase in numbers of existing ‘osmotrophic’ organisms, which in aquatic environments consist mostly of heterotrophic bacteria. The ecological result is known as the ‘viral shunt’ of the ‘microbial loop’, that is, fueling heterotrophic microbial metabolisms in part by lysing autotrophic and heterotrophic microbes [27].
While there is little empirical evidence indicating the extent to which the same microbe lysis-driven biogeochemical processes occur in soils, their occurrence in soils nonetheless seems probable. Process delays likely can exist, however, since the pore size within soils can allow nutrients to become entombed and thereby not immediately available for further microbial utilization, that is, until wetting of soils causes desorption [28,29]. Thus, we can hypothesize that viruses are important and perhaps even the main drivers of nutrient entombment within soils.
Like nutrients, a substantial fraction of virions in soils are thought to also be found adsorbed to abiotic soil materials [30]. Nutrient desorption from soil particles upon wetting therefore might be accompanied also by virion desorption. The resulting release of virions could give rise to additional infection and lysis of soil microbes, ‘pumping’ even more soluble nutrients into soils, analogous to the ‘biological pump’ in oceans [6]. Alternatively, we can speculate that, upon soil wetting, microbial replication and motility could bring cells to virions that have failed to desorb from soil particles, rather than virions desorbing and then diffusively moving toward cells. As a result, in soils some viruses could act as sit-and-wait (ambush) predators of microbes [31], rather than as diffusive ‘pursuit’ predators. Upon subsequent lysis, within wetted soils, substantial numbers of freely diffusing ‘pursuit’ virions may then be released, along with accompanying freely diffusing, lytically released ‘necromass’.

2.3. Viral Modification of Host Metabolism During Lytic Infections

When viruses infect a host cell, often they immediately redirect that cell’s metabolism toward production of virion progeny. This virus-mediated alteration in host biochemistry and physiology can directly impact microbial metabolic outputs. The extent of this impact on microbial metabolism can be higher if a virus carries auxiliary metabolic genes (AMGs). AMGs can represent more efficient versions of genes already used by microbes in their cellular metabolisms, though also can be genes which provide new functions. AMGs, though, are generally thought to be host-derived genes, representing a form of what originally were described as “vegetative viral genes” [32]. Most identified AMGs have been found to impact central carbon metabolism and photosystem II, thus providing an immediate metabolic enhancement over other non-virus-infected cells [10,33].
AMGs are widespread in oceans but seem to be rare in soils [34,35,36]. Glycoside hydrolase AMGs, however, were recently identified in viruses from organic-rich soils where these genes would help break down the complex organic matter present [34,35,36]. AMGs thus may be less rare in soils than previously thought, a notion that can be tested as soil viruses become more sampled and characterized.

2.4. Virus-Mediated Horizontal Gene Transfer

Every second in the oceans there are an estimated 1023 viral infections and these infections are thought to mediate approximately 1016 transduction events between cellular microbes [25,26]. Due to the substantially greater complexity of soil environments (e.g., soil spatial heterogeneity), an equivalent soil calculation is impossible to perform. In principle, though, cellular genetic material should be similarly transferable by soil viruses between microbes, likely greatly expanding soil microbe evolutionary potential [37,38]. Transduced genes thus can provide another mechanism by which viruses can impact ecosystems, one that is in addition to their ability to phenotypically modify cellular organisms (Section 2.3) before killing and lysing them (Section 2.2).
Transduction traditionally has been viewed as a form of accidental horizontal gene transfer [39]. This generally occurs due to virus DNA-handling errors that allow host ‘donor’ genetic material to become encapsidated in a virion. The resulting still structurally functional virions, once released, can then deliver their accidentally packaged genetic cargo to a new ‘recipient’ host. Transduction usually is differentiated into generalized transduction where viruses randomly encapsidate only host DNA vs. specialized transduction, where viruses encapsidate a combination of both host and viral DNA, while usually picking up only a specific portion of host DNA [40]; see also [41].
Specialized transduction [40] along with another virus-associated horizontal gene transfer mechanism known as lysogenic conversion—which is provirus-mediated modification of a cell’s phenotype that occurs during latent virus infections [42]—are both mediated solely by temperate viruses, i.e., ones capable of displaying these latent infections. Lysogenic converting genes, in contrast to donor–host genes being subject to specialized transduction, are considered to be normal components of phage genomes rather than recent accidental acquisitions. These genes can have substantial impacts on ecosystems, such as by encoding bacterial toxin genes [43]. Lysogenic converting genes are also related to, and in many cases even identical to, what are known as phage morons; extra phage genes acquired from hosts that are both stable constituents of virus genomes and expressed during virus infection cycles [44,45,46].

2.5. Many Environmental Viral States

Viruses lately have been conceptualized into two complementary states: free virions (extracellular virus particles) vs. virocells (viruses infecting host cells) [47]. It has long been known, though, that viruses are able to switch back and forth between these two states [48]. From a perspective of environmental virus microbiology, we can consider additional categories of viral states (Figure 1), and specific methods used to characterize environmental viruses will influence the degree to which each state is observed. This section presents this expanded, virus environmental-state framework (Figure 1), which builds on a simpler viewpoint considering proviruses vs. productive infections vs. free virions [34,35].
Virions are part of the encapsidated environmental fraction (category 1). Free virions usually are small in size and virions generally have genomes that are resistant to enzymatic degradation. Virions also are isolatable from unencapsidated materials and rich in viral nucleic acid.
Virocells include latent viral infections (category 2). These can either consist of host genome-integrated proviruses or plasmid proviruses. Integrated proviruses are linked to host-cell genes. Plasmid proviruses are somewhat less coupled with host genes though may be found in many copies both within individual cells and within environments. For both, the viral DNA is unencapsidated.
Virocells also include productive infections (category 3). Like plasmid proviruses, viral genomes undergoing productive infections usually are not physically coupled to host DNA. Unlike plasmid proviruses, productive infections in the near term are highly metabolically active, will typically generate relatively large numbers of newly replicated viral genomes, and also will generate new virions. As a result, category 3 will contain both numerous copies of a given viral genome and encapsidated nucleic acid. For lytic phages the latter will be found within particles (cells) that are much larger than individual virions.
‘Virus-like eDNA’ (vleDNA) (category 4) is extracellular, unencapsidated environmental DNA that has been derived in various ways from virus genomes [49]. This viral nucleic acid often is degradable using DNase and will not be physically linked to host-cell genes unless the vleDNA was derived from integrated proviruses. See Section 4.2 for further discussion of vleDNA.
We suggest an additional, catchall category of virus states that we describe simply as ‘Other’ (category 5). ‘Other’ contains viral genomes that are unencapsidated (contrasting category 1), not physically linked to host genes nor necessarily found in many copies either within cells or across environments (contrasting category 2). They are also not found in numerous copies within cells (contrasting category 3) and not derived from the extracellular environment (contrasting category 4). Examples include restricted virus genomes [50], virus infections that are unsuccessful for other reasons [51,52], viral genomes that are in a stalled pre-replicative state (i.e., ‘pseudolysogenic’) [53,54,55], or viral DNA that is contained within extracellular vesicles [56,57]. In addition, for some virus-like mobile genetic elements of fungi, no encapsidated extracellular states are even known [58,59].
Individual approaches to virus community characterization will tend to result in underassessments of virus presence or impact within environments as (i) not every viral state will be efficiently represented when using a single technique, (ii) not all detected virus-like nucleic acid will be from environmentally propagating viruses, and (iii) not all virions are easily propagated in vitro. Categories 1, 2, 3, and to some degree even 5 can, however, consist of propagatable virus nucleic acid and thereby may in principle be isolated as functional virions in the laboratory (Section 3.1). All five categories can be captured in metagenomes (Section 3.2). Only categories 1, 3, and to a smaller degree also 5 can contribute to viromes (Section 3.3). Thus, metagenomes and viromes will not consist solely of propagating viral nucleic acid, but depending on variations in processing, can permit eDNA to be captured, and not all encapsidated DNA is necessarily of viral origin (Section 4). In contrast, not all virions are easily propagated in vitro, so viromes and metagenomes will tend to capture a greater diversity of potentially propagatable viral nucleic acid than virus isolation can alone.

3. Three Ways to Characterize Soil Viruses

This section describes three different methods used to characterize soil viruses. We specifically consider the pros and cons associated with each approach and how different approaches can complement each other (Figure 2). This is to provide guidance especially to researchers with less expertise on soil viruses. These methods consist of virus isolation and subsequent laboratory propagation (Section 3.1), soil metagenomics (Section 3.2), and the generation and analysis of encapsidated subsets of metagenomes known as viromes, recently dubbed as viromics (Section 3.3).

3.1. Virus Isolation

Until the advent of metagenomics (Section 3.2 but see also Section 3.3), the characterization of environmental viruses first required obtaining a pure virus culture. Though essentially as old as virology itself, at least as a laboratory science, the isolation of viruses remains an important technique that offers a unique lens into understanding host–virus dynamics and can be essential toward fully characterizing a virus’ genotype and phenotype. The functions of most viral genes are unknown, and this is especially so for soil viruses [4]. The primary means of determining the function of a viral gene is by mutationally knocking out that gene and then examining the outcome of virus infection and replication [60]. Further, virus isolation enables measurement of infection metrics as burst size and latent period, and knowledge of those parameters is crucial to understanding the potential of a virus to impact an ecosystem.

3.1.1. Techniques for Isolating Viruses

Methods for isolating viruses from soil and other sources—especially bacteriophages—can be found in previous publications [61,62,63,64]. These include (i) direct isolation, (ii) isolation from soil wash, (iii) isolation following enrichment culture, (iv) isolation following virion concentration, and (v) isolation following induction of proviruses. For all of these approaches, after an appropriate incubation period microbial cells are removed by filtration or centrifugation (or both) and the now clarified fluid is tested for the presence of a virus. Testing typically is by plating to look for host killing as plaques, using a previously isolated indicator host [65,66]. A wealth of information and resources on culturing and characterizing virus isolates is available in the literature [67,68,69,70].
Direct plating, enrichment, and concentration: Virus isolation directly from soils or soil washes can involve simply plating using a previously isolated host strain as an indicator [71,72]. More commonly, especially when viruses are less abundant, a soil sample or wash may be incubated with a broth-cultured microbial isolate, a procedure described as enrichment culturing [61]. Even though enrichment is common practice, multiple attempts at enrichment still may be needed to obtain a single virus isolate. For example, isolation of Mycobacterium phages, one of the most well-cultured category of viruses, can involve 30 parallel enrichment cultures to yield one phage isolate [70,73]. This need for enrichment repetition is likely due to a combination of a high diversity of host species and strains within soils along with most viruses being somewhat host-species or even host-strain specific [51,64]. Polyvalent viruses, however, also exist, meaning that they are able to infect hosts from more than one host genus [74,75,76].
One can also first concentrate viruses after resuspending them from soil. This is then followed by filtration or precipitation, and only then subjecting the resulting concentrate to enrichment culturing [62]. That is, initiating enrichment cultures with potentially more virus particles from soil samples so that the number of enrichment cultures required per successful virus isolation is smaller. The initial virus resuspension step from soils is discussed more fully in Section 3.3.
Latently infecting viruses: A different approach to isolating viruses involves starting with proviruses infecting microbes isolated from soils [38,77]. In soils, approximately 30% of bacteria can harbor one or more prophages that could be induced to produce virions. The number of proviruses present in soils in fact may be even higher than that as not all proviruses are inducible using the typically employed mitomycin C [78,79]. Also, it is important to recognize that temperate phages, upon initial virion adsorption, often infect lytically rather than lysogenically [80], meaning that the same virus from the same environment could potentially be isolated using different isolation methods both as a provirus and as a virion. This ability of a phage to infect other than lysogenically commonly is described as their lytic potential, but it seems to be modifiable in response to how many hosts are present that virions can infect [81,82,83,84,85]. Further information on different proposed methods that bacteriophages use to regulate lysis-lysogeny conversion can be found in the literature (see [84,86,87,88,89,90,91,92]).
Culturing limitations: A major constraint on culturing many viruses is first growing the virus’ host in pure culture, as most microbes have not yet been cultured [93]. This limitation in host availability reduces the types of viruses that can be isolated and thus studied in pure culture. Even for hosts that can be grown in culture, not all can be grown to confluence on an agar plate, i.e., so as to support the growth of virus plaques, and for some viruses, even if their hosts will readily form lawns on agar surfaces, still will not form plaques under standard culture conditions [94,95].
For those viruses that grow poorly as plaques, other approaches may be used, at least for detection, including culture clearing (culture lysis in broth) [96,97] or routine test dilution (meaning culture clearing on a plate as near confluent lysis) [98]. Culture clearing in particular can be performed in multi-well plates in an automated system for high-throughput monitoring [99]. Finally, the original isolation host, especially as it typically will not have been isolated from the same sample as the virus isolate, is not always a primary host of a virus but instead may represent a sub-optimal host, leading to inaccurate estimation of growth parameters [100,101]. These various limitations on growing virus hosts in the laboratory make the isolation, propagation, and also ecological characterization of viral isolates in the laboratory challenging.

3.1.2. Well-Developed Soil-Virus Systems

A few soil virus–host systems have been particularly well developed, especially for phages and bacterial hosts. Buckling and colleagues, for example, used bacteriophage SBW25Φ2 and its soil-living host, Pseudomonas fluorescens SBW25, to study antagonistic coevolution between the host and phage [102] and the role of phages in host diversification [103,104]. Poisot and colleagues also used P. fluorescens SBW25 to isolate a variety of phages from soil washes [105]. They then looked at the range of hosts these bacteriophages could infect using bacteria co-isolated from the same soil washes, so as to examine the role of nutrient resources on the specialization of the phages. From these data, they concluded that soils which are more nutrient limited could contain phages with greater host specialization (narrower host range) than soils which, artificially, have been made more nutrient rich. Vos and colleagues [106] compared phage adaptation to specific hosts using bacteriophages and hosts that were isolated from the same soil samples. They found better adaptation of phages to local hosts—as indicated by infection rates of randomly isolated bacterial hosts—than to hosts isolated a greater distance away.
Among the most well-developed soil-virus systems are those infecting the bacterial genus, Mycobacterium. Mycobacteria include disease-causing along with harmless bacteria commonly found in soils. This makes mycobacteria medically relevant (e.g., Mycobacterium tuberculosis) as it consists of hosts for phage isolation that can be used with few biocontainment precautions (especially Mycobacterium smegmatis). The Science Education Alliance-Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) [107,108] and the Phage Hunters Integrating Research and Education (PHIRE) programs [109,110] successfully integrated phage isolation into mentoring young scientists and providing large collections of phages that infect Actinobacteria (the bacterial phylum that includes genus Mycobacterium). As a consequence, M. smegmatis-infecting phages represent the largest collection in the world of well-characterized virus isolates infecting a single microbial host. Recently from these efforts a patient was successfully treated with a cocktail of three phages able to infect an antibiotic-resistant strain of Mycobacterium abscessus, phages that were isolated using M. smegmatis [111].

3.1.3. Isolation of RNA Fungal Viruses

The methods described in the previous section presuppose that the viruses being isolated have an extracellular phase and are generally lytic to the cell. While this is true for many bacteriophages and archaeal viruses, including those with either a DNA or RNA genome, it is less frequently true for viruses of fungi, also known as mycoviruses. Mycoviruses have been identified in all major taxa of fungi, they predominantly have dsRNA genomes (although both ssRNA and ssDNA genome types exist), and for many no encapsidated extracellular states are known [58,112]. Mycovirus genomes can be isolated or otherwise identified by extracting all of the RNA from growing fungi [113] or instead using RT-PCR to target known mycovirus DNA sequences [114,115].

3.2. Metagenomics

Contrasting the procedures of isolation, which often focus on just a single virus or microbe clone, metagenomics involves extracting all of the DNA from a sample. The DNA is then broken up into many small fragments and sequenced (called shotgun sequencing). The resulting sequence is analyzed en masse to reconstruct the microbe and virus genomes present. As the process does not require culturing, and most microbes cannot be cultured (as noted above), it has greatly expanded our knowledge of microbes across many environments. As it also does not rely on PCR-based detection of a universal marker gene (e.g., the microbial 16S rRNA gene which viruses do not have), it has immensely increased our knowledge of what viruses are present within environments [116,117,118,119,120]. With metagenomics, the composition of environmental viral communities could be described for the first time, substantially accelerating the development of environmental viral ecology.
Contrasting most notably with marine environments, metagenomic studies have not been as successful in analyzing soil viruses. The cause of this deficit has stemmed mainly from low viral DNA-extraction yields, leading to sub-optimal virus genome assemblies. The consequence is poor characterization and detection of meaningful ecological connections between viruses and microbes. As a result, most soil metagenomic studies have disregarded rather than emphasized the viral component [16]. With the advancement of new technologies to amplify lower inputs of DNA, and more sophisticated bioinformatics to analyze the sequencing data, metagenomics for virus ecology in soils is, however, becoming more feasible [34].

3.2.1. Losing Sight of Virus Genes in a ‘Sea’ of Sequence

In this section we consider various challenging aspects to characterizing the viruses found in soil metagenomes. The basis of these challenges is that there are billions of microbes found in a single gram of soil [117,121]. The vastness of these numbers boosts the appeal of soil metagenomics over microbial isolation as it is impossible to isolate all or even most of these organisms. At the same time, the resulting over-abundance of data generated by metagenomic analyses can blur our ability to finely resolve each individual type of organism, and this has especially been an issue for resolving virus genomes. Nonetheless, two general approaches to improve the virus genome-resolving power of metagenomic analyses exist: improved sequencing depth and improved sequence analysis. In addition, there usually will exist biases in terms of what DNA is sequenced or even analyzed.
Sequencing biases: The major benefit of metagenomics stems from the relatively minimal wet lab work needed before sequencing. This is in comparison with virus isolation and subsequent characterization (Section 3.1) or with separation of encapsidated nucleic acid before sequencing for viromic analyses (Section 3.3). Metagenomic analysis of a random sample of DNA nonetheless will result in biases stemming from: (i) differing abundances of community members, (ii) the specific manner in which samples are collected and stored, (iii) the physical and chemical methods used to extract and subsequently amplify DNA (the latter if applicable), and, as considered also in this section, (iv) what bioinformatic tools are used to reconstruct the metagenome [117,122]. Many of these biases, however, can be reduced with implementation especially of more standardized methodologies [123,124,125,126].
Sequencing depth: The enormous diversity of microbes and our inability to physically capture all of the DNA from a soil sample—the latter as resulting, for example, from inefficiencies in microbe lysis and DNA collection—make it impossible to sequence all of the DNA present. The DNA collected is also fragmented, which along with its high diversity makes it difficult to assemble sequenced fragments into complete or even near-complete microbial genomes, where the former, completely assembled genomes, are called metagenome-assembled genomes, or MAGs. In addition, less abundant microbial genomes tend to be not even nearly completely sequenced. The net result is that a metagenome once constructed will not be identical to the actual collection of nucleic acid sequences found in the original sample. In an effort to overcome these issues, the number of sequencing reads for a sample can be increased, which should allow increased recovery of MAGs with lower error (Section 3.2.2). This approach, however, can make assembly of abundant organisms harder due to sequencing errors that can mimic within-species (micro)diversity [127,128].
Similar assembly challenges exist for virus genomes, even though these are generally small (many thousands of base pairs) compared to the genomes of microbes, which are typically much larger, generally several millions of base pairs [129]. In addition, viruses often are not sufficiently abundant in environments to make up for resulting differences in terms of total sequenceable DNA. Target theory [130] thus would predict a lower likelihood that a given sequencing read would ‘hit’ a given viral genome vs. a given microbial genome. For example, as a thought experiment, consider a ‘metagenome’ constructed from only a single sequencing read. The likelihood of that read being of a specific virus genome would be equal, all else held constant, to the total amount of DNA associated with that virus population relative to the total amount of DNA present within a sample; for example, about a trillion base pairs, such as 105 × 107 (virus genome size times number of genomes of a single virus type) vs. quadrillions of base pairs, such as 5 × 106 × 109 (microbe genome size also times number of genomes of a specific microbe type). Even with somewhat more sequencing coverage, these larger genomes still can figuratively act as ‘haystacks’, obscuring virus genome ‘needles’ due to there being many more sequencing reads from microbes than from viruses. The result generally tends to be far less virus sequence and far fewer virus genomes generated in metagenomes than is the case for microbes.
Sequence analysis: Virus detection within metagenomes is further hampered by the often vast diversity of viruses present, which can make de novo assembly of viral contiguous sequences (contigs) challenging [131], that is, assembly without employing already sequenced viral genomes as templates. Specifically, the main and best assembly algorithms are based on overlapping stretches of sequenced nucleotides (i.e., De Bruijn graph assembly [124,131]), and overlapping stretches become rarer the lower the number of copies of specific viral DNA sequences that is originally present in a sample. Indeed, less than 2% of assembled sequences are typically of virus origin [34,132]. The result is decreased virus-sequence detection within metagenomes along with assembly of only partial virus genomes.
Metagenomes also are bioinformatically intensive to assemble and annotate, which can also interfere with virus identification and assembly. In particular, adding more virus detection-and-characterization bioinformatic steps can be unrealistic during metagenome analyses. Furthermore, in attempting to hunt for viruses within a metagenomic ‘sea’, it can quickly become apparent that virus identification itself can be non-trivial and particularly so as often most predicted genes have no annotation (Section 3.2.3) and so consequently can be difficult to assign to viruses. In total, the resulting incomplete bioinformatic ‘snapshot’ of what viruses are present and what specifically their genomes consist of means that virus sequence derived from metagenomes will tend to less readily reveal the functionality of what viruses are present within a sequenced environment.

3.2.2. Vertical Coverage

The concept of sequencing coverage can be used in two ways, horizontal vs. vertical (Figure 3). Horizontal coverage, also known as coverage breadth, refers to what portion of a contig or genome has reads aligning to it at least once, and this is often used to know how complete an assembled genome is relative to a reference genome. For MAGs, which by definition lack a reference genome, researchers rely instead on the identification of universal marker genes to estimate completeness. [133].
Vertical coverage, also known as coverage depth or sequencing depth, is by contrast the average number of reads that align to a base in a contig or an assembled genome. Vertical coverage is often used as a measure of the relative abundance of microbes or viruses within environments and can be used to determine how reliable some analyses are; for example, to assess single nucleotide polymorphisms in a microbial genome you need at least 15× coverage of that base pair [134,135]. Generally speaking, the greater the vertical coverage, the better. For instance, less abundant viruses and less abundant microbes can be missed in studies with too ‘shallow’ vertical coverage because a sequence consensus cannot be reached, and this can impact metagenome diversity estimates.
It unfortunately is difficult to interpret the specifics of ecology from the vertical coverage of sequencing reads. Although the relative abundance of a virus in a metagenome, resulting in the potential for greater vertical sequencing coverage, suggests greater impact by those viruses on an ecosystem, higher abundance does not necessarily determine its impact in more qualitative terms. In addition, metagenomes are only a snapshot of a community and cannot provide information on community dynamics (changes over time) unless generated over a time series, an issue which is not addressed by improving the vertical coverage of only a single sample. It also simply is expensive to increase sequencing depth.

3.2.3. Drawing Information from Bulk Sequence

The primary challenge in metagenomic analysis of bulk DNA to study virus ecology is one of distinguishing viral genomic sequence from background cellular genomic sequence (Section 3.2.1). Major advancements in making these distinctions have been made including identifying viral hallmark genes (VirSorter [136]) or virus-specific motifs (VirFinder/DeepVirFinder [137,138]) to recognize likely viral contigs in metagenomes. Both types of approaches can be performed using publicly available and user-friendly programs found on CyVerse [139] or KBase [140]. These websites provide many advantages to help in the study of viruses. For example, they employ graphical user interfaces (i.e., GUIs as familiarly seen with modern computers and smart phones) rather than command line controls (the latter, e.g., as seen with the original DOS-based personal computers from the 1980s, where ‘DOS’ stands for disc operating system). These GUIs list hundreds of applications (apps) for processing metagenomic data along with the previous versions of those apps, and the user can sort these apps, for example, by topic or function. Additionally, optimal parameters are suggested, making analyses easier to perform, and step-by-step instructions for many of the applications are provided [141,142,143]. Each method has recommended conservative settings but also more-encompassing sensitive settings, determined by how likely the identified sequence represents a genuine virus (greater sensitivity, i.e., results in lower likelihood).
VirSorter, for example, uses multiple lines of evidence to place a contig into a category. Categories include 1 (“most confident”), 2 (“likely”), and 3 (“possible”) for viruses not integrated into a host’s genome or plasmid, with categories 4 to 6 the equivalent for integrated proviruses [123]. VirSorter relies on a database of known viral genes for category prediction. Due to this, it works especially well for marine viruses as they are better characterized genomically, but databases may be improved for virus detection in soils by the addition of new genome sequences of soil viruses, including as following virus isolation (Section 3.1). DeepVirFinder also relies on a virus reference database like VirSorter, but uses a machine learning approach in its database that enables robust detection of virus fragments ≥3 kb, with a conservative approach (likely a virus) selecting contigs with a score ≥0.9 and p-value <0.05 and a sensitive approach (probable virus) ≥0.7 and p-value <0.05. VirSorter and DeepVirFinder can also be used in parallel to optimize viral identification from metagenomic data.
A unique benefit of using metagenomic approaches is the ability to assemble viral and microbial genomes from the same data. MAGs can also be interrogated to identify proviruses. Proviruses found in high-coverage microbial genomes, integrated or not, could have increased coverage (allowing more robust analyses like micro-diversity) over viruses in other states (Section 2.5), simply due to their higher environmental prevalence within highly prevalent microbes [144]. MAGs and identified partial or complete viral genomic sequences (contigs) can be matched using several different approaches (e.g., using spacers in clustered regularly interspaced short palindromic repeats, also known as CRISPRs, and shared nucleotide identity [145]) to identify associations within the same cells and thereby possibly provide information on virus replication lifestyle [146].
Metatranscriptomic datasets can also be obtained through shotgun sequencing of RNA templates and searched for RNA viruses. While often used for assessing gene expression, genomes can be assembled from metatranscriptomes using similar pipelines as for metagenomes, and RNA virus and phage genomes can then be identified in these assemblies, including for soil samples [147]. Importantly, current pipelines including VirSorter and DeepVirFinder are not optimal for RNA virus detection due to (i) a limited number of references for environmental RNA viruses and (ii) fundamental differences in genome structure and gene content for RNA viruses; hence viral sequence mining from metatranscriptomes still requires a substantial amount of manual inspection and curation. One feature that unites all RNA viruses and can aid in their detection and characterization is their RNA-dependent RNA polymerase (RdRp). RdRps are proteins that catalyze the replication of RNA from an RNA template and are essential to RNA viruses. Analyses of RdRps consequently may provide insights into the diversity of RNA viruses and their putative hosts [148].

3.2.4. Outlook

Metagenomics is a powerful approach that has provided numerous insights into the characteristics of uncultured microbes and viruses, along with their possible interactions, and it continues to grow in terms of utilization. The first metagenomics papers only analyzed a small fraction of the microbial data collected and only minimal information about viruses was obtained. With the development of new computational tools and advancements in machine learning [149], however, we are now at a time where virus discovery and exploration can be performed by anyone who generates or has access to a metagenome. Notably, as newer tools become available and metagenomes become routinely generated, their sample collection and analysis needs to be more thorough [120]. This includes in terms of how samples are collected from the environment, how those samples are stored, and then how resulting sequences are documented in terms of meta data. Overall, the metagenomic approach for studying virus ecology is suitable especially for initial characterization of a soil ecosystem, for soil studies aimed at microbial diversity more generally (i.e., beyond ‘just’ viruses), for inferring possible virus–host interactions, and for those just getting their ‘feet wet’ in terms of the study of the omics of soil virus ecology.

3.3. Viromics

A virome is a metagenome that consists, ideally, solely of sequence data obtained from the VLP fraction of environments (encapsidated environmental nucleic acid, also known as a VLP metagenome). Viromes are generated by separating VLPs from microbial cells, lysing those particles, and then sequencing the released nucleic acid. A virome thus can be thought of as a ‘targeted metagenome’, one which focuses on a specific aspect of a metagenome to better describe that fraction’s specific taxonomic content and related characteristics. The first virome, published in 2002 [150], was derived from marine water and since then this approach has become the dominant method for characterizing viruses across many environments [5,10].

3.3.1. Utility and Drawbacks of Viromes

The main advantage viromes have over mining viral signals from less targeted metagenomes is that there is increased coverage specifically of viral genomes. That increased coverage is possible because of prior removal of the DNA of microbes and macroscopic eukaryotes. The latter, as noted (Section 3.2.1), have larger genomes that as a result are represented by a large portion of sequencing reads. The increased vertical coverage afforded by targeting the VLP fraction of biomes for sequencing therefore can yield more complete viral genomes. Consequently, greater horizontal coverage (Section 3.2.2) can increase the diversity of viruses captured, and can reveal micro-diversity within viral populations [151]. All of these benefits accumulate into complete or near-complete viral genomes that can subsequently be used as reference genomes. Reference genomes (i) are (useful for identifying new viruses from metagenomic/viromic data, (ii) can provide viral taxonomic affiliation, and (iii) can better allow for prediction of viral gene functions. Obtaining the virus fraction for targeted sequencing from soils, however, is not without challenges (Section 3.3.2).
Still, this viromics approach has many of the same drawbacks described for metagenomic studies (Section 3.2): biases associated with sample preparation; high expense (due simply to the large number of sequencing reads, though this expense is continually declining); being bioinformatically intensive [120]; and that most predicted genes have no annotation. Unlike untargeted metagenomic approaches, where DNA is extracted en masse from soil, with viromics more wet lab work is required to separate VLPs from various forms of unencapsidated soil DNA before VLP-associated DNA or nucleic acid generally can be extracted.

3.3.2. The Challenge of Separating Virions from Soils

The virome approach has only recently emerged as a viable option in soils, as dramatic differences between aquatic and soil environments, e.g., physical structure of soil, previously have prevented aquatic virome generation protocols from being translatable to soils. Particularly, the problem has been one of virion adsorption to soil matrix and difficulties associated with virion desorption from that matrix, at least in vitro during virome preparation. Separating VLPs from the soil matrix thus is the greatest challenge for characterizing the soil virosphere compared to viromes obtained from less complex environments, and in practice this is a time-consuming and laborious process.
More than 90% of soil viruses are estimated to be adsorbed to the soil matrix [30], and desorbing them can be tricky. There are many forces acting on viruses in soils, but the virion’s isoelectric point—the environmental pH that causes a virion to have no net surface charge—is the primary factor in determining their adsorption to soil matrix [152]. It is currently impossible, however, to determine the isoelectric point of all of the virions in a soil sample. Consequently, to desorb virions, various chemical reagents with different charges and physical methods are employed [5]. Virus desorption methods in particular should be tailored to specific soil types [35,153,154,155]. We therefore first suggest characterizing a soil to understand its anion/cation-exchange capacity (a measure of how many ions can be retained on soil particle surfaces [156]) and, separately, the diversity of the associated microbial community (e.g., via 16S rRNA gene surveys). The latter is because the isoelectric point of viral proteins can be strongly correlated with the isoelectric point associated with their host’s proteins as has been calculated for many microbes [157].

3.3.3. Additional Sources of VLP Losses

VLP desorption from the soil matrix is typically followed by filtration for size fractionation. Viruses range in size from tens of nanometers (nm) to several hundred nm in diameter, making the filtration step a major point of bias, especially since most viromes are generated from viruses that have passed through a 220-nm filter [158]. This specific filter size targets phages which are typically ~50 nm in size. Because 220 nm refers to the maximum pore size, however, it is likely that even virions smaller than that cutoff may not pass through, especially larger virions [159]. In addition, if there is a lot of debris on the filter, viruses can adsorb to that rather than passing through the filter. Larger filters (≥450 nm) have also been used, but less frequently due to fear of microbe contamination (Section 4). After filtration, virions are concentrated and then these steps (i.e., chemical and physical desorption, filtration, and concentration) are repeated on the original soil sample multiple times serially to increase yields.
During these processes there is the additional issue of virions degrading or adsorbing to other surfaces after being desorbed from the soil. Virions often will adhere to every new surface encountered including those associated with the container that holds them, thereby also decreasing yields. Soils that contain a lot of organic matrix material, e.g., humic substances, are particularly difficult in that they contain an array of surface charges and matrix particles, making it hard to both desorb virions and keep them resuspended. Adsorption to organic matrix material not only changes how the virions appear, as adsorption to organic material can make it harder to identify a virion microscopically [160] (more in Section 3.3.4), but also can complicate downstream processing. For example, organic material often can bind DNA, keeping DNA in the organic layer. That DNA, as a result, is then removed from the sample during some DNA extraction methods [161,162].

3.3.4. Efficiency of Virus Resuspension from Soils

The proportion of viruses desorbed from a soil describes a given VLP resuspension method’s virus resuspension efficiency. Different chemical and physical desorption methods can be compared by either enumerating VLPs which are endogenous to a soil sample or instead by recovering a known amount of exogenously added virus particles (the latter, also known as spike-in experiments). In virus ecology, VLP enumeration via microscopic direct counts generally is accomplished via either epifluorescence microscopy (EFM) or transmission electron microscopy (TEM). Direct quantification of VLPs from soils, however, typically is inconsistent between determinations, resulting in high variability between technical replicates and across microscopy techniques [163,164]. In this section we compare these two microscopy techniques as associated and additional approaches to determining virus resuspension efficiency.
Epifluorescence microscopy: EFM is the most widely used environmental-virus direct-count method because it is quick (sample preparation and enumerating accomplished in ~1 h, depending on the number of samples), extremely sensitive (the dyes involved strongly bind to dsDNA and RNA), and is less expensive than TEM [165]. The technique involves a combination of nucleic acid-binding fluorescent dyes and excitatory ultraviolet light that results in visualization of pinpricks of emitted light that individually correspond to VLPs. All the dyes used for EFM, however, will bind to any nucleic acid, although many preferentially bind to dsDNA. This binding promiscuity can mean that many things in a soil sample can ‘light up’ as VLPs during EFM, including DNA contained in extracellular vesicles [56,57] and other ‘fake’ viral particles [166]; see Section 4.1, Section 4.2 and Section 4.3 for other non-virus entities that may be part of the VLP fraction. In addition, the dye fades quickly upon exposure to the ultraviolet light (as known as ‘bleaching’), limiting the time over which a sample can be observed, although this fading can be slowed by using an antifade solution [167].
Transmission electron microscopy: TEM provides higher resolution than EFM, permitting visualization of both virus–cell interactions and virion morphology. In addition, samples can be viewed multiple times, leading to improved precision and accuracy. TEM, however, is much more expensive and time-consuming than EFM (2–3 times longer for the same sample size). It is also not available everywhere and requires considerable expertise. Furthermore, even with the high resolution afforded by TEM, it can be difficult to distinguish true viruses from non-virus particles of similar size, as will typically be found in environmental samples (i.e., non-viral or ‘fake’ VLPs). For an in-depth overview of TEM capabilities for viruses, see [168].
Spiking in functional virions: Both EFM and TEM can be applied to samples with either endogenously or exogenously supplied viruses. Here, we use the term, “spike-in”, to describe exogenously supplied viruses. With spike-ins, virus recovery is typically measured via enumeration of plaque forming units (PFUs). In this case, the recovery of these known viruses acts as a proxy for recovery of all viruses in the soil sample [30,169]. Unfortunately, PFU enumeration is a functional rather than direct measurement, which can be misleading for efficiency determinations as it relies on virus infectivity rather than being a measurement of the absolute quantity of virus particles. In particular, viruses which become inactivated during resuspension without necessarily also losing their viral genomes will not be counted in the course of PFU enumerations, though nevertheless still will contribute to viromes.
Detecting spiked-in encapsidated nucleic acid: To focus efforts on quantification of virus particles rather than their infectivity, that is, rather than detection of PFUs, sequence-specific DNA probes that are tagged with a fluorescent dye can be designed to specifically target virions that have been spiked in, with their abundance measured via qPCR. This approach works well for quantifying known virus pathogens in the environment [170,171], but is not directly representative of native soil viruses due to these spiked-in viruses being added to an environment in which they are not endemic. These nonindigenous viruses will have different adsorption coefficients (how quickly they adsorb to surfaces) and different avidities (overall adsorption strength) for soil constituents. Thus, while quantification of these added viruses is possible, it does not necessarily translate into how well the resuspension process captures the native environmental viruses and as a result soil spiked-in approaches generally are insufficiently quantitative.
Bioinformatic approach: A different measure can be used as a metric to compare estimated efficiencies of virus resuspension after nucleic acid has been sequenced. This involves bioinformatically calculating the amount of sequencing (i.e., number of reads) of identified viruses compared to the total amount of sequencing for the sample, providing a ratio of known virus sequence to total sequenced nucleic acid [35]. A virus resuspension method can be applied to many soils or samples and the ratios determined by this approach can be compared to evaluate how well the virus resuspension method captured viruses vs. contamination. While this approach is not quantitative, since it does not measure the total number of viruses present in a soil sample, it does provide a rough measure of how much non-viral contamination is in each sample and allows comparisons of different resuspension methods and bioinformatic approaches.

3.3.5. Outlook

Characterizing viruses as identified from viromes has become a dominant method in the marine realm, but soil viromic efforts to date have been less rewarding. The relatively limited number of efforts to isolate viruses from soils or characterize viral genomes from soils via metagenomics or viromics has left soil-virus genetic diversity largely unknown. As a result, with each new soil virus study a majority of viruses tend be novel, which is challenging due to the difficulties in assigning functions to otherwise unknown and uncharacterized genes [172]. This usually results in insufficient recognition of virus genomes or of individual virus genes even if previously sequenced, making untargeted metagenomic studies, in particular, less worthwhile for virome characterization. The result is minimal representation of soil viruses in current virus databases, with only ~10% of sequences in the virus RefSeq database [173] (v92) and ~3% in the Integrated Microbial Genome/Virus (IMG/VR) database [174] (v4) arguably representing soil viruses. To overcome the issue of databases mostly lacking in soil-virus sequence and the corresponding large quantity of unknown sequences in soil viromes, reference-sequence independent approaches are emerging that allow comparison of viromes that help to provide insights into spatial and temporal viral diversity [175].
In marine environments, viromics has enabled robust analyses of environmental viruses and their potential impacts on local and global ecosystems [10,176,177]. The hope is that viromics may allow equivalent characterization of virus populations in soils as has been much more readily achieved in non-soil environments. This, however, will likely be achieved only in direct association with improvements in efficiencies of virion desorption from soil matrices. No work as of now—via any of the approaches described here (isolation, metagenomics, viromics)—has characterized every type of virus that may be found in the same sample from any environment, soil or otherwise.

4. Metagenomic Dataset Contaminants

Above we discuss key areas for improvement of de novo assembly (Section 3.2.1), coverage (Section 3.2.2), identification of viral sequences (Section 3.2.3), and virus enrichment (Section 3.3) from metagenomics data (for more on these subjects, see [178]). In this section we consider various forms of ‘contamination’ of metagenomic data, that is, any environmental entities, particularly but not only VLPs, that possess a reasonable likelihood of resembling an active virus within a soil sample. Included, in order of further discussion, are: (i) non-infectious virus-like particles (niVLPs; Section 4.1), (ii) eDNA (Section 4.2), (iii) microbe contamination (Section 4.3), (iv) amplification artifacts (Section 4.4), and (v) ecologically inactive or ‘banked’ virions (Section 4.5).

4.1. Non-Infectious Virus-Like Particles (niVLPs)

Among niVLPs are otherwise intact virions which are no longer capable of successfully infecting a host, should hosts become available (contrast with phage banks; Section 4.5). Included among niVLPs are also VLPs that are not of virus origin, such as gene transfer agents (GTAs; Section 4.3.2). Non-GTA niVLPs are true virion particles which are no longer infectious due to (i) non-wholly catastrophic capsid structural damage (as still allowing inclusion in VLP direct counts but not in virus viable counts), (ii) having faulty genetic material (i.e., lethal mutations or nucleic-acid structural damage), (iii) possessing virion maturation errors (existing as incompletely formed virions) [179], (iv) having become irreversibly attached to soil components in a manner that renders them no longer cell absorbable [180], and (v) which lack genetic material due to injection into a host cell or accidental ejection into the extracellular environment. The latter, now virus capsids lacking in genetic material, would not be detected in a metagenome or virome, but could inflate VLP counts particularly as determined by TEM; dyes for TEM, such as phosphotungstic acid hematoxylin or uranyl acetate, that is, stain the virus capsid material rather than necessarily nucleic acid whereas all EFM dyes would not cause nucleic acid-lacking particles to fluoresce. See Section 3.3.4 for more on virus detection using microscopy.

4.2. Extracellular DNA (eDNA/relic DNA)

The vast majority of eDNA is from microbes and is ubiquitous in soils where it can play a number of ecological roles including serving as a nutrient source, as a component of biofilm matrices, or as a mediator of the horizontal gene transfer mechanism called transformation (i.e., uptake of eDNA such as by microbes). Because eDNA can persist for prolonged periods (then also known as relic DNA), it thereby can obscure our ability to characterize soil ecosystems as they exist in terms of what genomes currently are active. Though eDNA was first thought to come primarily from lysed cells, it was later determined also to be secreted by microbes [181], though it may be released from decaying virions as well; the latter a form of vleDNA (Section 2.5).
The persistence of eDNA is particularly problematic because all of the DNA within a sample is extracted to generate a metagenome. This includes from cellular organisms and viruses but also any eDNA that has persisted. It was recently shown that relic DNA in particular has the potential to inflate microbial richness estimates up to 55% depending on the soil’s geochemical parameters (e.g., pH [182]). It was also recently shown in an aquatic ecosystem that eDNA accounted for about 60% of the total sequenced DNA and that a comparison of eDNA sequences to virome sequences revealed viruses that were only detected in the eDNA samples, implying that vleDNA was present [49].

Removing eDNA from Environmental Samples

In virome generation, virion purification techniques are incorporated to remove non-encapsidated DNA (e.g., DNase treatment to remove eDNA), but nevertheless non-encapsidated DNA is still detected [183,184]. One reason for this is that eDNA, including vleDNA, can be bound to inorganic or organic compounds that can prevent its degradation. Likewise, DNase requires divalent metals for activation and the presence of inorganic (e.g., copper sulfate) and organic compounds in a sample can bind divalent metals, thus partially or completely inhibiting DNA degradation [181,185]. In both cases—eDNA being protected or DNase activity being blocked—the proportion of eDNA contamination persisting past DNase purification depends on the soil composition [186,187,188].
New methods have been proposed to remove eDNA in the laboratory during preparation of metagenomes [182,189] or otherwise predict biases resulting from relic DNA on microbial community structure via modeling [190]. One new method to remove eDNA incorporates propidium monoazide (PMA), which is a photoreactive DNA-binding dye that can enter through pores in cell membranes, binding only to either eDNA or DNA that is found in dead cells. After a short incubation under light, bound PMA modifies DNA, preventing downstream processing, i.e., by blocking amplification and sequencing. PMA has also been used to inactivate DNA associated with damaged viruses [191], though this technique tends to be almost exclusively applied to samples containing viruses that infect humans (see [192] for a comprehensive list of studies). Presumably this technique would not work to remove all niVLPs, because, as noted, a VLP could become non-infectious due to defects in genes or nucleic acid structure rather than due to pores in capsids (Section 4.1). Nevertheless, PMA treatment is still useful, as one environmental metagenomic study, in which samples were collected from a clean-room floor, found that removal of relic DNA allowed detection of microbes and viruses that were not otherwise detected due to their low prevalence relative to that of relic DNA [155]. Once PMA treatment is performed, or any method of eDNA removal, qPCR can be used with 16S or 18S rRNA gene primers [184] to check to see if microbial DNA is still present in a virome, either because microbial relic DNA was not removed during PMA treatment or microbial DNA remained within intact ultrasmall microbial cells (Section 4.3).

4.3. Microbe-Derived Virome Contamination

Removal or even identification of microbial contamination in a virome is not as straightforward as it may seem. VLPs can carry rRNA genes, which is the most common way to detect microbial contamination in a virome [193] and thus genuine VLP DNA may be mistaken for more direct microbial contamination. Alternatively, microbe contamination can be incorrectly inferred if the sequences of actual viruses share similarities to known microbial sequences (e.g., AMGs or specialized transducing particles; Section 2.3 and Section 2.4), which may lead to removal of sequences during bioinformatic processing and thereby loss of legitimate virus data. On the other hand, microbial DNA may represent contamination stemming from the presence of ultrasmall or dormant microbes possessing decreased cell size (Section 4.3.1), GTAs (Section 4.3.2), or even plasmids that may or may not encode virus genomes (Section 4.3.3). Though not discussed further here, note that the converse of small microbes being similar in size to typical VLPs, is large VLPs being similar in size to typical microbes [194].

4.3.1. Ultrasmall Microbes

Ultrasmall microbes, those that can pass through a 0.45-µm filter, and some even 0.2-µm filters [195], are widespread and are found in the Bacteria and Archaea domains. Though not the same, ultramicrocells also exist, which are microbes with reduced cellular size due to dormancy as may be induced for various reasons including starvation [196]. The small cell size of ultrasmall microbes is matched by small genomes that do not include non-essential DNA, resulting in reduced functional potential [197]. Many ultrasmall microbes actually do not have enough metabolic capability to survive in isolation, e.g., as due to some missing complete housekeeping biochemical pathways. Instead, they join with other microbes to form metabolic networks.
Part of bioinformatic virus detection (described in Section 3.2.3) is identifying motifs typically exhibited in viruses, including enrichment of uncharacterized genes or possession of short genes, things which ultrasmall microbes can also exhibit [197]. Ultrasmall microbes thus can be similar to viruses in their genomic properties, which can make them a challenging virome contaminant to remove. Ultrasmall microbes nonetheless are not likely to be present in appreciable quantities in metagenomes for the same reasons that many viruses are also not present in appreciable quantities (i.e., their smaller genomes) (Section 3.2.1). In addition, unlike viruses, ultrasmall microbes can be detected bioinformatically because of their 16S rRNA genes. Nevertheless, due to their small size, ultrasmall microbes can be mistaken for viruses during EFM-based direct counts (Section 3.3.4), thereby inflating perceived VLP numbers.

4.3.2. Gene Transfer Agents (GTAs)

Marine viromes have been rigorously optimized yet still can contain presumptive cellular DNA contamination comprising approximately one third of metagenomic sequencing. This genetic material presumably is of non-viral origin and otherwise is thought to consist mostly of DNA carried by GTAs [198,199]. GTAs are non-viral though nevertheless are VLPs, containing pieces of genetic material obtained from the genome of the microbe they originated from and which they can transfer to other, similar microbes. GTAs, unlike true viruses, cannot however directly create progeny GTAs [200]. GTAs nonetheless have been proposed to be atypical, genetically defective viruses, or viruses that have been otherwise repurposed by a host particularly to horizontally transfer host DNA [201]. In any case, GTAs represent a form of niVLPs. Indeed, to detect GTAs, many studies have focused on the genome sequence characteristics of what microbes are most likely to produce GTAs along with common genotypic characteristics found among identified GTAs, i.e., particularly possession of few if any known viral genes [200].
Currently, more GTAs have been identified from Alphaproteobacteria than any other group of cellular microbes [200]. This consequently presents a potential problem for soil viromics because these bacteria are some of the most abundant microbes found in soils [202]. For example, a recent study proposed that GTA-associated genetic material, based on sequence similarity to Alphaproteobacteria DNA, can represent up to 25% of assembled reads from viromes generated from peat soils [35].

4.3.3. Plasmids

Plasmids are extrachromosomal, semi-autonomous, either circular or linear pieces of DNA, and they are present in most microbes [203]. They regularly encode genes that are non-essential to their cellular hosts, i.e., as known as accessory genes. Most notably from a medical microbiology perspective, this accessory genetic material includes antibiotic resistance genes. Plasmids can move between microbes during conjugation (particularly bacteria connecting via sex pili, effecting DNA movement), via transformation, and by transduction [204] (for more on transduction, see Section 2.4).
Plasmids and viruses can have many similar genes, especially for DNA replication and interaction with host defenses [205,206]. Plasmid DNA sequences, unlike those of viruses generally, are also common in metagenomes and present problems for virus identification as undertaken via automated viral detection, i.e., as due to plasmids encoding virus-like genes [136,207]. For instance, VirSorter (Section 3.2.3) detects viruses based on viral hallmark genes, which can also be picked up by microbial hosts and transferred into plasmids. Discerning between a virus- and a plasmid-encoded virus-like gene within a metagenome also can be difficult because most of the genes in question may be unknown and only genes that are known and previously associated with viral genomes may be described with any certainty as viral genes. Thus, a plasmid with virus-like genes can easily be identified as a virus. New bioinformatic tools, however, are being developed to detect plasmids in metagenomics datasets either for removal or for use in plasmid-focused investigations [208,209,210].
In terms of plasmid inclusion in viromes, it is important to note that plasmids are not encapsidated and are thereby mostly excluded from viromes. Plasmids also can represent a component of eDNA, but like vleDNA, plasmids should be excludable from viromes to the extent that eDNA is removed, for example by DNase treatment, prior to removing encapsidated DNA from viral capsids.

4.4. Amplification Artifacts

Viruses with single-stranded DNA (ssDNA) genomes are diverse, ubiquitous, and infect all domains of life including numerous microbial taxa [211]. The study of ssDNA viruses has arguably benefitted the most from metagenomics due to its greatly expanding the number of known ssDNA viruses, cataloging the hosts they infect, and highlighting their environmental roles, all paving the way for a global analysis of ssDNA viruses and their importance [212]. Even with the advent of metagenomics, however, ssDNA viruses are tricky to both detect and study because they can have segmented genomes, which can appear as separate viruses in metagenomic datasets. Investigations are further impeded because ssDNA viruses undergo rapid mutation while evidence supports widespread horizontal gene transfer [212]. The importance of ssDNA viruses in environments nevertheless may be overstated.
The paradigm in question is that ssDNA viruses are the most abundant virus type in soils. This conclusion, however, appears to have arisen partially because of the need to greatly amplify viral DNA for the sake of generating sufficient quantities for sequencing. Many whole genome amplification methods have been used to overcome the issue of low DNA yield extracted from environmental samples (e.g., random priming-mediated sequence-independent single-primer amplification [159,213]), but multiple displacement amplification (MDA) was the most widespread whole genome amplification method implemented until the 2010s. This technique uses rolling-circle amplification, which has been shown to preferentially amplify circular ssDNA, including that of plasmids, while unevenly amplifying linear genomes [214]. The result is a biased inflation of the abundance of ssDNA viruses in samples, making their actual abundance unknown and quantitative comparisons to other datasets thereby impossible. Thus, while ssDNA viruses are of interest, they are unlikely to be as prevalent as earlier reports suggested.
While it has taken some time, whole genome amplification methods are being replaced with methods that quantitatively capture ssDNA viruses and permit ecological comparisons between ssDNA and dsDNA viruses, providing a more holistic view of the soil virosphere. The first important development was to optimize the first step of traditional high-throughput DNA sequencing protocols (adapter ligation) to allow for PCR amplification and accurate sequencing of both ssDNA and dsDNA [215]. The initial aim of these modifications was to increase the accuracy of sequencing, and because most living things have dsDNA genomes (ssDNA viruses, of course, excepted), adapters were designed to aid in the sequencing of both strands of DNA from a single molecule [215,216]. A library method also was recently developed that included novel adapter attachment chemistry, which permits quantitative amplification and sequencing of ssDNA, dsDNA, and damaged DNA in parallel [217]. To test its fidelity, this method was first applied to mock viral communities [214] and since has been shown to capture both ssDNA and dsDNA viruses in many environments including soil from picogram-level input DNA [155,218]. Library preparation kits and protocols able to generate quantitative metagenomes from nanogram DNA inputs are thus now readily available and should be primarily used, as opposed to non-quantitative amplification.
To aid in the detection of ssDNA viruses, many studies have utilized the fact that the majority of ssDNA viruses are circular and encode known marker genes, such as homologs of genes encoding the rolling-circle replication-associated protein [36,155,206,219]. To date, there has only been one study that quantitatively amplified both ssDNA and dsDNA viruses from the soil samples [155]. Using known ssDNA virus marker genes, it suggested that ssDNA viruses were a small fraction of the microbial viruses observed (~4%). To fully evaluate this ‘ssDNA viruses are dominant in soils’ paradigm, however, additional quantitatively amplified soil viromes are needed that evaluate the relative abundances of ssDNA to dsDNA viruses, with a careful consideration of contaminants, as contaminants can increase perceived abundances of dsDNA viruses, e.g., all known GTAs and ultrasmall microbes carry dsDNA rather than ssDNA [200].

4.5. Ecologically Inactive Viruses

Functionally active but nevertheless ecologically ‘inactive’ viruses can be described as being in a ‘Bank mode’, as equivalent to ‘Seed banks’ for plant populations. This ability vs. inability to potentially cause future infections distinguishes, respectively, banked viruses from niVLPs. The banked mode concept further proposes that only the most abundant viruses within an environment are likely actively replicating [220].
In soils, viruses in banked mode arguably exist as two different subcategories: (i) functionally active viruses that cannot reach a host for a variety of reasons including reversible adsorption to soil matrix, and (ii) functionally active viruses that infect only rare hosts. In the first case, these viruses could become ecologically active when environmental conditions change; for example, when rainfall creates channels in the soil matrix permitting movement (Section 2.2). In the latter case, these viruses are always ecologically active, because their hosts remain, filling a niche. Actually, banked viruses are still ecologically important, at least over longer time frames, because they can help maintain the diversity of viruses, but are nearly impossible to distinguish within metagenomes or viromes from viruses that are more ecologically active. They might be distinguishable instead in metatranscriptomes, time-course experiments, or in experiments where active viruses are labeled (e.g., stable isotope probing [221]).

5. Conclusions

The still young field of soil virus ecology deserves continued and indeed enhanced attention as soils are a central component of many of Earth’s biomes, and viruses are increasingly recognized as important to ecosystem functioning. Different approaches to the study of virus ecology, however, have not been equivalently developed. This can result, in some cases, in intellectual biases where seemingly ‘better’ data come to dominate thinking even if it also underlies different, potentially competing, and not necessarily superior perspectives. Nonetheless, and despite disparate efforts to date, our understanding of the roles of viruses in soils remains meager on nearly all fronts, leading to the functional and ecological importance of viruses in soils to be largely overlooked.
In virus ecology, intellectual biases can perhaps be seen especially in terms of genotypic (sequence-based) characterizations vs. characterizations that are more phenotype-based. This can particularly be the case since sequence-based characterization is often easier to perform and certainly can provide far more data that are more straightforward to analyze using computers. The study of soil virus ecology nevertheless may be relatively unique in this regard in that sequence-based virus analyses, particularly viromics, can also be somewhat difficult to perform with soils, owing to the complexity of the soil environment physically, chemically, and spatially. That is, despite the growing torrent of sequence-based soil viromics data, its role in our understanding of soil virus ecology remains somewhat underdeveloped.
Here we outlined various approaches to undertaking both phenotypic and genotypic characterizations of soil viruses, including the challenges and solutions, with emphasis on improving sequence-based characterizations. Given that soils lag behind other environments in terms of the development of viromics and virus ecology, an important near-term emphasis should be on improving omics approaches in soils and consideration of viruses in all soil microbiome studies.

Author Contributions

G.T. was the primary writer of the manuscript. S.R. and P.H. provided editing and made contributions to some of the writing. S.T.A. was invited by Soil Systems to provide the article, contributed to the writing, and otherwise served as senior author on the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Portions of this work were written under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under contract DE-AC52-07NA27344. The work conducted by the U.S. Department of Energy Joint Genome Institute is supported by the Office of Science of the U.S. Department of Energy under contract no. DE-AC02-05CH11231.

Acknowledgments

We graciously thank Noah Sokol, Rachel Hestrin, Eric Slessarev, Aram Avila Herrera, and Dinesh Adhikari for their feedback on the manuscript. We thank Ryan Goldsberry for making illustrations for the graphical abstract and Figure 2. A special thanks to Ella Sieradzki for doing a courtesy review of the manuscript.

Conflicts of Interest

S.T.A. generated and maintains various virus ecology-emphasizing web pages, including phage.org (the Bacteriophage Ecology Group). The other authors declare no conflicts of interest.

Appendix A. Glossary of Terms

Antagonistic coevolution. Interactions between two species in which evolutionary adaptations in one species negatively affect a second, resulting in evolution of counter adaptations by the second species; for viruses this is typically seen as a coevolutionary arms race where host organisms evolve virus resistance which is then overcome by virus adaptations.
Auxiliary metabolic gene (AMG). Gene encoded by a virus that was acquired from a previous host organism, and that can be expressed during virus infections to alter an infected cell’s metabolic activity over that of uninfected cells.
Bank mode. Refers to virions that are dormant but not inactive, particularly due to current lack of access to absorbable host organisms, but with the dormant state potentially reversible once an absorbable host appear.
Biogeochemistry. Biological, chemical, geological, and physical processes that occur in an environment particularly involving movements of nutrients within and between ecosystems.
Biome. A community of organisms occupying a major habitat.
Chronic release. Virus infections in which virions are released without substantial disruption of host cells, for instance as via virion extrusion or virion budding across or from the host-cell envelope, in contrast to lytic infection.
Community. Multiple species living together in a given area.
Confluent lysis. The inability to delineate where one plaque ends and another begins, making a plate appear to be completely covered by interconnected plaques. This is typically the result of plating too many virus particles that are too numerous to count.
Contiguous sequence (contig). Referring to nucleic acid sequences that are adjacent within the genome of a single organism; sequencing reads that can be assembled into a larger genome fragment are ones which are contiguous.
Coverage. Bioinformatics term that describes the extent to which an assembled genome has sequencing reads that map to it, either across the genome or to a specific region; this can be differentiated into coverage breadth (or horizontal coverage) vs. coverage depth (or vertical coverage).
Coverage breadth. Proportion of a genome to which sequencing reads align, that is, the fraction of a genome that has been successfully sequenced; also known as horizontal coverage.
Coverage depth. Number of sequencing reads that map to a specific region of the genome, that is, the degree of sequencing redundancy achieved; also known as vertical coverage.
De novo assembly. The assembly of a contig using an algorithm, instead of assembly using a reference genome.
Desorption. The release or detachment of a substance or particle from a surface.
Ecogenomics. Methods of determining ecological characteristics and interactions from genome-sequence information.
Enrichment culture. Technique toward amplifying microorganisms of specific phenotype from an environmental sample; consists, for viruses, of adding the sample that might contain a virus to media along with specific host cells to allow amplification of virus numbers.
Environmental DNA (eDNA). DNA that is present in an environment outside of a biological entity, which for most microbes is extracellular.
Epifluorescence microscopy (EFM). Imaging technique that uses a microscope that emits light in ultraviolet wavelengths to cause fluorescence of parts of a specimen.
Extraction. See viral extraction.
Gene transfer agent (GTA). Virus-like particle that is not biased toward packaging the DNA responsible for producing it but rather packages all cellular DNA with roughly equivalent probability.
Generalized transduction. Process by which DNA is moved from one host to a different host due to a virus accidentally, randomly encapsidating host DNA without associated viral DNA; contrast with specialized transduction and gene transfer agents.
Hallmark genes. Genes in a viral genome that are central to virus replication and structure, and are shared by a broad variety of viruses, but are missing from cellular genomes.
Horizontal coverage. Synonymous to coverage breadth.
Horizontal gene transfer. Movement of genetic material between organisms other than in the course of either reciprocal sexual gene exchange or vertically from parent to offspring; virion-mediated horizontal gene transfer generally is called transduction.
Induction. As pertaining to proviruses, the transition from an established latent cycle to a productive infection, including as can be forced, e.g., as in the course of mitomycin C treatment of bacterial lysogens.
Integration. Process of insertion of a provirus’ genome into existing host genomic DNA, within a host cell, as toward establishment of a latent infection; integrated proviruses become physically linked to host genetic material; contrast with plasmid provirus.
Isolation. See viral isolation.
Latent infection. Virus infection during which virion progeny is not produced but viral genome replication occurs.
Library. In the context of metagenomics, a library is a DNA template prepared for sequencing, including following amplification of DNA to adequate levels for sequencing.
Lysogen. Especially a bacterium harboring a prophage; that is, a bacterium hosting a lysogenic cycle.
Lysogenic conversion. Virus-encoded modification of a cell’s phenotype that occurs during latent virus infections and is not a result of normal virus functioning. Lysogenic conversion is not directly associated with retention of the latent-infection state; lysogenic cycle repressor genes, for example, therefore are not also converting genes.
Lysogenic cycle. Ongoing, especially bacteriophage existence as a prophage; a bacteriophage latent infection.
Lytic cycle. Productive viral infection which ends with virion release via host-cell lysis.
Maturation error. Failure of virion components to properly assemble into an infectious virus particle.
Metagenome. Collection of sequences obtained from untargeted sequencing of all nucleic acids extracted from a biome sample.
Metagenomics. Non-culturing set of method to extract, sequence, and analyze a portion of all nucleic acid from a biome sample.
Metagenome-assembled genome (MAG). Near complete to complete genomes of organisms assembled from metagenome sequence information.
Microbial loop. The movement of nutrients, especially carbon, from a dissolved state in an environment up through multiple microbial trophic levels, particularly movement from dissolved organic carbon to heterotrophic bacteria to protozoa.
Micro-diversity. Genetic diversity among individuals within a population (same species).
Mitomycin C. Chemical that alkylates DNA and forms cross-links, causing significant cytotoxicity to cells and resulting in SOS responses and associated induction of proviruses.
Mobile genetic elements. Any entity that moves nucleic acid between loci either within or between cells or organisms.
Moron. Gene acquired from host cells especially as are expressed by viruses during latent cycles but whose function is not necessarily directly related to virus metabolism during latent or productive infections, i.e., as representing ‘more’ DNA.
Multiple displacement amplification (MDA). An amplification technique that uses a polymerase isolated from phi29 bacteriophage to generate sufficient quantities of DNA for sequencing.
Necromass. Total mass associated with dead organisms in an environment or sample.
niVLP. Non-infectious virus-like particles.
Non-infectious virus-like particle (niVLP). VLP that is incapable of infecting a host organism.
Osmotrophic. Referring to organisms obtaining energy and nutrients from dissolved environmental materials, e.g., with heterotrophic bacteria and fungi serving as key osmotrophic organisms in soil environments.
Plasmid provirus. A virus that replicates separate from host chromosomes while latently infecting.
Predator. Organism that kills other (prey) organisms in order to obtain nutrients from that other organism’s now-dead body.
Prophage. Bacteriophage provirus.
Productive infection. Viral infection in which new progeny virions are produced and released, the latter either lytically or chronically depending on the virus.
Provirus. Latently infecting virus genome as present in a host cell.
Pseudolysogeny. Virus infection which has stalled including as due to nutrient limitations but that is capable of restarting toward either a productive or latent infection.
Read. Short for sequencing read, i.e., genotype information of an organism obtained via one individual nucleic acid sequencing process.
Reference genome. A representative example of an organism’s nucleotide sequence.
Relic DNA. DNA that has been preserved in an environment in a non-functional form over extended time periods, e.g., more than seconds, minutes, or hours; see also, for example, niVLPs and vleDNA as well as eDNA.
Restricted infection. Virus infection that cannot be completely executed, thus interfering with virus propagation but not necessarily in which the virus-infected cell is inactivated/killed; for example, as mediated by bacterial restriction-modification systems.
Resuspension. See viral resuspension.
Richness. Number of different populations (different species) found in a given area; short for species richness.
Rolling-circle amplification. In vitro nucleic-acid replication process in which multiple copies of a circular template are generated by a polymerase using one nucleic-acid strand as template while displacing the other strand, i.e., as based on rolling-circle replication.
Sequencing depth. Synonymous to coverage depth.
Shotgun Sequencing. A method where DNA is broken up into many small fragments, which are then sequenced in parallel to obtain multiple overlapping reads to determine the original DNA sequence.
Soil wash. Process where a buffer solution is added to a soil sample, mixed, and the sample then centrifuged, with resulting supernatant recovered.
Specialized transduction. Process where DNA flanking an integrated provirus is encapsidated after an error in provirus excision with it along with virus genomic material then transferred to a new host cell; contrast with generalized transduction.
Strictly lytic. Virus that upon infection is inherently unable to display either latent cycles or chronic release; synonymous with obligately lytic.
Structural damage. Irreversible physical disruption of a virion particle capsid or appendages; contrast with genomic mutation or nucleic-acid damage.
Targeted metagenome. Metagenome generated with specific, biasing steps to focus on a subset of a community; viromes, for example, are targeted metagenomes.
Temperate virus. Virus that can perform both latent and productive replication cycles, though not both at the same time.
Transduction. Process of horizontal gene transfer between cells that is virus effected.
Transmission electron microscopy (TEM). Technique that uses a beam of electrons rather than light to illuminate a specimen and thereby create high resolution images (micrographs); TEM can be used to visualize virus particles within environmental samples.
Viral extraction. A process to lyse viral capsids to release the DNA using a combination of physical and chemical methods.
Viral isolation. A process whereby a single virus is propagated in the laboratory in association with its host.
Viral metagenome. Targeted metagenome focused on viral (or VLP) nucleic acid sequences from a biome sample.
Viral resuspension. Process to desorb virions from soil using a combination of physical and chemical methods.
Vertical coverage. Synonymous with coverage depth.
Viral shunt. Solubilization of cellular organisms, especially microbes, via virus-induced lysis, thereby preventing or delaying energy and organic carbon movement from these organisms to higher trophic levels.
Virion. A complete infectious virus particle including nucleic acid, a capsid, and sometimes an envelope.
Virome. Synonymous with viral metagenome, i.e., a metagenome that has been biased toward sequencing of the VLP portion of biomes.
Viromics. Targeted metagenome method with specific steps to sequence especially viral nucleic acid from a biome.
Virosphere. All of the viruses found in a given area.
Virus-like eDNA (vleDNA). Environmental DNA that is either from or thought to be from a virus, i.e., eDNA of probable virus origin.
Virus-like gene. Genetic material not necessarily explicitly from a virus source that is the best match to a known virus gene and/or which is localized with nearby virus genes to a specific strand of DNA.
Virus-like particle (VLP). Particles of virus size as found in an environmental sample that potentially contains viral nucleic acid, i.e., something that probably is a virion particle but is not necessarily a virion particle; an alternative definition, from the medical virology literature and not used here, is a virus capsid that lacks viral nucleic acid.
Virus-specific motif. Nucleic acid sequence pattern that is indicative of a virus; see also virus-like gene, i.e., the best match to a known virus gene, or otherwise known virus nucleic acid pattern.
VLP metagenome. Synonymous with viral metagenome.

References

  1. Hyman, P.; Abedon, S.T. Viruses of Microorganisms; Caister Academic Press: Norwich, UK, 2018. [Google Scholar]
  2. Witzany, G. Viruses: Essential Agents of Life; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  3. Mushegian, A.R. Are there 1031 virus particles on Earth, or more, or less? J. Bacteriol. 2020. [Google Scholar] [CrossRef]
  4. Williamson, K.E.; Fuhrmann, J.J.; Wommack, K.E.; Radosevich, M. Viruses in soil ecosystems: An unknown quantity within an unexplored territory. Annu. Rev. Virol. 2017, 4, 201–219. [Google Scholar] [CrossRef] [PubMed]
  5. Pratama, A.A.; van Elsas, J.D. The ‘neglected’ soil virome-potential role and impact. Trends Microbiol. 2018, 26, 649–662. [Google Scholar] [CrossRef] [PubMed]
  6. Kuzyakov, Y.; Mason-Jones, K. Viruses in soil: Nano-scale undead drivers of microbial life, biogeochemical turnover and ecosystem functions. Soil Biol. Biochem. 2018, 127, 305–317. [Google Scholar] [CrossRef]
  7. Williamson, K.E. Viruses of microorganisms in soil ecosystems. In Viruses of Microorganisms; Hyman, P., Abedon, S.T., Eds.; Caister Academic Press: Norwich, UK, 2018; pp. 77–93. [Google Scholar]
  8. Abedon, S.T.; Murray, K.L. Archaeal viruses, not archaeal phages: An archaeological dig. Archaea 2013, 2013, 251245. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Hyman, P.; Abedon, S.T. Smaller fleas: Viruses of microorganisms. Scientifica 2012, 2012, 734023. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Breitbart, M.; Bonnain, C.; Malki, K.; Sawaya, N.A. Phage puppet masters of the marine microbial realm. Nat. Microbiol. 2018, 3, 754–766. [Google Scholar] [CrossRef] [PubMed]
  11. Steward, G.F.; Culley, A.I.; Mueller, J.A.; Wood-Charlson, E.M.; Belcaid, M.; Poisson, G. Are we missing half of the viruses in the ocean? ISME J. 2013, 7, 672–679. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Weitz, J.S.; Li, G.; Gulbudak, H.; Cortez, M.H.; Whitaker, R.J. Viral invasion fitness across a continuum from lysis to latency. Virus Evol. 2019, 5, vez006. [Google Scholar] [CrossRef] [Green Version]
  13. Raoult, D.; Forterre, P. Redefining viruses: Lessons from Mimivirus. Nat. Rev. Microbiol. 2008, 6, 315–319. [Google Scholar] [CrossRef]
  14. Dupré, J.; O’Malley, M.A. Varieties of living things: Life at the intersection of lineage and metabolism. Philos. Theor. Biol. 2009, 1, e003. [Google Scholar] [CrossRef] [Green Version]
  15. Forterre, P. Defining life: The virus viewpoint. Orig. Life Evol. Biosph. 2010, 40, 151–160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Fierer, N. Embracing the unknown: Disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 2017, 15, 579–590. [Google Scholar] [CrossRef]
  17. Kozloff, L.M. Biochemical studies of virus reproduction. VII. The appearance of parent nitrogen and phosphorus in the progeny. J. Biol. Chem. 1952, 194, 95–108. [Google Scholar] [PubMed]
  18. Thingstad, T.F.; Bratbak, G.; Heldal, M. Aquatic phage ecology. In Bacteriophage Ecology; Abedon, S.T., Ed.; Cambridge University Press: Cambridge, UK, 2008; pp. 251–280. [Google Scholar]
  19. Mahmoudabadi, G.; Milo, R.; Phillips, R. Energetic cost of building a virus. Proc. Natl. Acad. Sci. USA 2017, 114, E4324–E4333. [Google Scholar] [CrossRef] [Green Version]
  20. Jacquet, S.; Zhong, X.; Peduzzi, P.; Thingstad, T.F.; Parikka, K.J.; Weinbauer, M.G. Virus interactions in the aquatic world. In Viruses of Microorganisms; Hyman, P., Abedon, S.T., Eds.; Caister Academic Press: Norwich, UK, 2018; pp. 115–141. [Google Scholar]
  21. Hobbs, Z.; Abedon, S.T. Diversity of phage infection types and associated terminology: The problem with ‘Lytic or lysogenic’. FEMS Microbiol. Lett. 2016, 363, fnw047. [Google Scholar] [CrossRef] [Green Version]
  22. Forde, S.E.; Thompson, J.N.; Bohannan, B.J.M. Adaptation varies through space and time in a coevolving host–parasitoid interaction. Nature 2004, 431, 841–844. [Google Scholar] [CrossRef] [PubMed]
  23. Chao, L.; Levin, B.R.; Stewart, F.M. A complex community in a simple habitat: An experimental study with bacteria and phage. Ecology 1977, 58, 369–378. [Google Scholar] [CrossRef]
  24. Heilmann, S.; Sneppen, K.; Krishna, S. Sustainability of virulence in a phage-bacterial ecosystem. J. Virol. 2010, 84, 3016–3022. [Google Scholar] [CrossRef] [Green Version]
  25. Suttle, C.A. Marine viruses-major players in the global ecosystem. Nat. Rev. Microbiol. 2007, 5, 801–812. [Google Scholar] [CrossRef]
  26. Keen, E.C. A century of phage research: Bacteriophages and the shaping of modern biology. BioEssays 2015, 37, 6–9. [Google Scholar] [CrossRef] [PubMed]
  27. Wilhelm, S.W.; Suttle, C.A. Viruses and nutrient cycles in the sea: Viruses play critical roles in the structure and function of aquatic food webs. BioScience 1999, 49, 781–788. [Google Scholar] [CrossRef] [Green Version]
  28. Liang, C.; Schimel, J.P.; Jastrow, J.D. The importance of anabolism in microbial control over soil carbon storage. Nat. Microbiol. 2017, 2, 17105. [Google Scholar] [CrossRef] [PubMed]
  29. Barnard, R.L.; Blazewicz, S.J.; Firestone, M.K. Rewetting of soil: Revisiting the origin of soil CO2 emissions. Soil Biol. Biochem. 2020, 107819. [Google Scholar] [CrossRef]
  30. Hurst, C.J.; Gerba, C.P.; Cech, I. Effects of environmental variables and soil characteristics on virus survival in soil. Appl. Environ. Microbiol. 1980, 40, 1067–1079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Abedon, S.T. Ecology of anti-biofilm agents II. bacteriophage exploitation and biocontrol of biofilm bacteria. Pharmaceuticals 2015, 8, 559–589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Barksdale, L.; Arden, S.B. Persisting bacteriophage infections, lysogeny, and phage conversions. Ann. Rev. Microbiol. 1974, 28, 265–299. [Google Scholar] [CrossRef]
  33. Hurwitz, B.L.; U’Ren, J.M. Viral metabolic reprogramming in marine ecosystems. Curr. Opin. Microbiol. 2016, 31, 161–168. [Google Scholar] [CrossRef]
  34. Emerson, J.B.; Roux, S.; Brum, J.R.; Bolduc, B.; Woodcroft, B.J.; Jang, H.B.; Singleton, C.M.; Solden, L.M.; Naas, A.E.; Boyd, J.A.; et al. Host-linked soil viral ecology along a permafrost thaw gradient. Nat. Microbiol. 2018, 3, 870–880. [Google Scholar] [CrossRef]
  35. Trubl, G.; Jang, H.B.; Roux, S.; Emerson, J.B.; Solonenko, N.; Vik, D.R.; Solden, L.; Ellenbogen, J.; Runyon, A.T.; Bolduc, B.; et al. Soil viruses are underexplored players in ecosystem carbon processing. mSystems 2018, 3, e00076-18. [Google Scholar] [CrossRef] [Green Version]
  36. Jin, M.; Guo, X.; Zhang, R.; Qu, W.; Gao, B.; Zeng, R. Diversities and potential biogeochemical impacts of mangrove soil viruses. Microbiome 2019, 7, 58. [Google Scholar] [CrossRef] [PubMed]
  37. Zeph, L.R.; Onaga, M.A.; Stotzky, G. Transduction of Escherichia coli by bacteriophage P1 in soil. Appl. Environ. Microbiol. 1988, 54, 1731–1737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Ghosh, D.; Roy, K.; Williamson, K.E.; White, D.C.; Wommack, K.E.; Sublette, K.L.; Radosevich, M. Prevalence of lysogeny among soil bacteria and presence of 16S rRNA and trzN genes in viral-community DNA. Appl. Environ. Microbiol. 2008, 74, 495–502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Redfield, R.J. Do bacteria have sex? Nat. Rev. Genet. 2001, 2, 634–639. [Google Scholar] [CrossRef] [Green Version]
  40. Schneider, C.L. Bacteriophage-mediated horizontal gene transfer: Transduction. In Bacteriophages: Biology, Technology, Therapy; Harper, D.R., Abedon, S.T., Burrowes, B., McConville, M., Eds.; Springer: New York City, NY, USA, 2017; Available online: https://0-link-springer-com.brum.beds.ac.uk/referenceworkentry/10.1007/978-3-319-40598-8_4-1 (accessed on 18 April 2020).
  41. Chen, J.; Quiles-Puchalt, N.; Chiang, Y.N.; Bacigalupe, R.; Fillol-Salom, A.; Chee, M.S.J.; Fitzgerald, J.R.; Penades, J.R. Genome hypermobility by lateral transduction. Science 2018, 362, 207–212. [Google Scholar] [CrossRef] [Green Version]
  42. Los, M.; Kuzio, J.; McConnell, M.R.; Kropinski, A.M.; Wegrzyn, G.; Christie, G.E. Lysogenic conversion in bacteria of importance to the food industry. In Bacteriophages in the Control of Food- and Waterborne Pathogens; Sabour, P.M., Griffiths, M.W., Eds.; ASM Press: Washington, DC, USA, 2010; pp. 157–198. [Google Scholar]
  43. Christie, G.E.; Allison, H.A.; Kuzio, J.; McShan, M.; Waldor, M.K.; Kropinski, A.M. Prophage-induced changes in cellular cytochemistry and virulence. In Bacteriophages in Health and Disease; Hyman, P., Abedon, S.T., Eds.; CABI Press: Wallingford, UK, 2012; pp. 33–60. [Google Scholar]
  44. Hendrix, R.W.; Lawrence, J.G.; Hatfull, G.F.; Casjens, S. The origins and ongoing evolution of viruses. Trends Microbiol. 2000, 8, 504–508. [Google Scholar] [CrossRef]
  45. Cumby, N.; Davidson, A.R.; Maxwell, K.L. The moron comes of age. Bacteriophage 2012, 2, 225–228. [Google Scholar] [CrossRef]
  46. Taylor, V.L.; Fitzpatrick, A.D.; Islam, Z.; Maxwell, K.L. The diverse impacts of phage morons on bacterial fitness and virulence. Adv. Virus Res. 2019, 103, 1–31. [Google Scholar] [PubMed]
  47. Forterre, P. The virocell concept and environmental microbiology. ISME J. 2013, 7, 233–236. [Google Scholar] [CrossRef] [PubMed]
  48. Delbrück, M. Bacterial viruses or bacteriophages. Biol. Rev. 1946, 21, 30–40. [Google Scholar] [CrossRef]
  49. Mohiuddin, M.; Schellhorn, H.E. Spatial and temporal dynamics of virus occurrence in two freshwater lakes captured through metagenomic analysis. Front. Microbiol. 2015, 6, 960. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Schiffer, J.T.; Aubert, M.; Weber, N.D.; Mintzer, E.; Stone, D.; Jerome, K.R. Targeted DNA mutagenesis for the cure of chronic viral infections. J. Virol. 2012, 86, 8920–8936. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Hyman, P.; Abedon, S.T. Bacteriophage host range and bacterial resistance. Adv. Appl. Microbiol. 2010, 70, 217–248. [Google Scholar]
  52. Labrie, S.J.; Samson, J.E.; Moineau, S. Bacteriophage resistance mechanisms. Nat. Rev. Microbiol. 2010, 8, 317–327. [Google Scholar] [CrossRef] [PubMed]
  53. Miller, R.V.; Day, M. Contribution of lysogeny, pseudolysogeny, and starvation to phage ecology. In Bacteriophage Ecology; Abedon, S.T., Ed.; Cambridge University Press: Cambridge, UK, 2008; pp. 114–143. [Google Scholar]
  54. Abedon, S.T. Disambiguating bacteriophage pseudolysogeny: An historical analysis of lysogeny, pseudolysogeny, and the phage carrier state. In Contemporary Trends in Bacteriophage Research; Adams, H.T., Ed.; Nova Science Publishers: Hauppauge, NY, USA, 2009; pp. 285–307. [Google Scholar]
  55. Los, M.; Wegrzyn, G. Pseudolysogeny. Adv. Virus Res. 2012, 82, 339–349. [Google Scholar] [PubMed]
  56. Soler, N.; Marguet, E.; Verbavatz, J.M.; Forterre, P. Virus-like vesicles and extracellular DNA produced by hyperthermophilic archaea of the order Thermococcales. Res. Microbiol. 2008, 159, 390–399. [Google Scholar] [CrossRef] [PubMed]
  57. Nolte-’t, H.E.; Cremer, T.; Gallo, R.C.; Margolis, L.B. Extracellular vesicles and viruses: Are they close relatives? Proc. Natl. Acad. Sci. USA 2016, 113, 9155–9161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Vainio, E.J.; Hantula, J. Fungal viruses. In Viruses of Microorganisms; Hyman, P., Abedon, S.T., Eds.; Caister Academic Press: Norwich, UK, 2018; pp. 193–209. [Google Scholar]
  59. Sutela, S.; Poimala, A.; Vainio, E.J. Viruses of fungi and oomycetes in the soil environment. FEMS Microbiol. Ecol. 2019, 95, fiz119. [Google Scholar] [CrossRef] [Green Version]
  60. Stahl, F.W. Amber mutants of bacteriophage T4D: Their isolation and genetic characterization. Genetics 2012, 190, 831–832. [Google Scholar]
  61. Van Twest, R.; Kropinski, A.M. Bacteriophage enrichment from water and soil. Meth. Mol. Biol. 2009, 501, 15–21. [Google Scholar]
  62. Wommack, K.E.; Williamson, K.E.; Helton, R.R.; Bench, S.R.; Winget, D.M. Methods for the isolation of viruses from environmental samples. Meth. Mol. Biol. 2009, 501, 3–14. [Google Scholar]
  63. Lobocka, M.; Hejnowicz, M.S.; Gagala, U.; Weber-Dabrowska, B.; Wegrzyn, G.; Dadlez, M. The first step to bacteriophage therapy: How to choose the correct phage. In Phage Therapy: Current Research and Applications; Borysowski, J., Miêdzybrodzki, R., Górski, A., Eds.; Caister Academic Press: Norfolk, UK, 2014; pp. 23–67. [Google Scholar]
  64. Hyman, P. Phages for phage therapy: Isolation, characterization, and host range breadth. Pharmaceuticals 2019, 12, 35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Czajkowski, R.; Ozymko, Z.; Lojkowska, E. Isolation and characterization of novel soilborne lytic bacteriophages infecting Dickeya spp. biovar 3 (‘D. solani’). Plant Pathol. 2013, 63, 758–772. [Google Scholar] [CrossRef]
  66. Anne, J.; Wohlleben, W.; Burkardt, H.J.; Springer, R.; Puhler, A. Morphological and molecular characterization of several actinophages isolated from soil which lyse Streptomyces cattleya or S. venezuelae. J Gen. Microbiol. 1984, 130, 2639–2649. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Clokie, M.R.J.; Kropinski, A.M. Bacteriophages. Methods and Protocols. Volume 1: Isolation, Characterization, and Interactions; Humana Press: New York, NY, USA, 2009; Volume 501. [Google Scholar]
  68. Clokie, M.R.J.; Kropinski, A.M. Bacteriophages. Methods and Protocols. Volume 2: Molecular and Applied Aspects; Humana Press: New York, NY, USA, 2009; Volume 502. [Google Scholar]
  69. Clokie, M.R.J.; Kropinski, A.M.; Lavigne, R. Bacteriophages: Methods and Protocols. Volume 3. Methods and Protocols; Springer protocols (Series); Humana Press; Springer: New York, NY, USA, 2018; Volume 1681. [Google Scholar]
  70. The actinobacteriophage database at PhagesDB.org. 2020. Available online: https://www.phagesdb.org (accessed on 18 April 2020).
  71. Dhar, B.; Singh, B.D.; Singh, R.B.; Srivastava, J.S.; Singh, V.P.; Singh, R.M. Occurrence and distribution of rhizobiophages in Indian soils. Acta Microbiol. Pol. 1979, 28, 319–324. [Google Scholar] [PubMed]
  72. Rombouts, S.; Volckaert, A.; Venneman, S.; Declercq, B.; Vandenheuvel, D.; Allonsius, C.N.; Van Malderghem, C.; Jang, H.B.; Briers, Y.; Noben, J.P. Characterization of novel bacteriophages for biocontrol of bacterial blight in leek caused by Pseudomonas syringae pv. porri. Front. Microbiol. 2016, 7, 279. [Google Scholar] [CrossRef] [Green Version]
  73. Russell, D.A.; Hatfull, G.F. PhagesDB: The actinobacteriophage database. Bioinformatics 2017, 33, 784–786. [Google Scholar] [CrossRef] [Green Version]
  74. Ackermann, H.-W.; DuBow, M.S. Bacteriophage Taxonomy. In Viruses of Prokaryotes. Volume I. General Properties of Bacteriophages; Ackermann, H.-W., DuBow, M.S., Eds.; CRC Press: Boca Raton, FL, USA, 1987; pp. 13–28. [Google Scholar]
  75. Ross, A.; Ward, S.; Hyman, P. More is better: Selecting for broad host range bacteriophages. Front. Microbiol. 2016, 7, 1352. [Google Scholar] [CrossRef] [Green Version]
  76. Yu, P.; Mathieu, J.; Li, M.; Dai, Z.; Alvarez, P.J. Isolation of polyvalent bacteriophages by sequential multiple-host approaches. Appl. Environ. Microbiol. 2016, 82, 808–815. [Google Scholar] [CrossRef] [Green Version]
  77. Williamson, K.E.; Schnitker, J.B.; Radosevich, M.; Smith, D.W.; Wommack, K.E. Cultivation-based assessment of lysogeny among soil bacteria. Microb. Ecol. 2008, 56, 437–447. [Google Scholar] [CrossRef]
  78. Chen, F.; Wang, K.; Stewart, J.; Belas, R. Induction of multiple prophages from a marine bacterium: A genomic approach. Appl. Environ. Microbiol. 2006, 72, 4995–5001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Daly, R.A.; Roux, S.; Borton, M.A.; Morgan, D.M.; Johnston, M.D.; Booker, A.E.; Hoyt, D.W.; Meulia, T.; Wolfe, R.A.; Hanson, A.J.; et al. Viruses control dominant bacteria colonizing the terrestrial deep biosphere after hydraulic fracturing. Nat. Microbiol. 2019, 4, 352–361. [Google Scholar] [CrossRef]
  80. Sinha, V.; Goyal, A.; Svenningsen, S.L.; Semsey, S.; Krishna, S. In silico evolution of lysis-lysogeny strategies reproduces observed lysogeny propensities in temperate bacteriophages. Front. Microbiol. 2017, 8, 1386. [Google Scholar] [CrossRef] [PubMed]
  81. Avlund, M.; Dodd, I.B.; Semsey, S.; Sneppen, K.; Krishna, S. Why does phage play dice? J. Virol. 2009, 83, 11416–11420. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Waller, A.S.; Yamada, T.; Kristensen, D.M.; Kultima, J.R.; Sunagawa, S.; Koonin, E.V.; Bork, P. Classification and quantification of bacteriophage taxa in human gut metagenomes. ISME J. 2014, 8, 1391–1402. [Google Scholar] [CrossRef]
  83. Maslov, S.; Sneppen, K. Well-temperate phage: Optimal bet-hedging against local environmental collapses. Sci. Rep. 2015, 5, 10523. [Google Scholar] [CrossRef] [Green Version]
  84. Erez, Z.; Steinberger-Levy, I.; Shamir, M.; Doron, S.; Stokar-Avihail, A.; Peleg, Y.; Melamed, S.; Leavitt, A.; Savidor, A.; Albeck, S.; et al. Communication between viruses guides lysis-lysogeny decisions. Nature 2017, 541, 488–493. [Google Scholar] [CrossRef]
  85. Abedon, S.T. Look who’s talking: T-even phage lysis inhibition, the granddaddy of virus-virus intercellular communication research. Viruses 2019, 11, 951. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Ghosh, D.; Roy, K.; Williamson, K.E.; Srinivasiah, S.; Wommack, K.E.; Radosevich, M. Acyl-homoserine lactones can induce virus production in lysogenic bacteria: An alternative paradigm for prophage induction. Appl. Environ. Microbiol. 2009, 75, 7142–7152. [Google Scholar] [CrossRef] [Green Version]
  87. Hargreaves, K.R.; Kropinski, A.M.; Clokie, M.R. What does the talking? quorum sensing signalling genes discovered in a bacteriophage genome. PLoS ONE 2014, 9, e85131. [Google Scholar] [CrossRef] [Green Version]
  88. Abedon, S.T. Commentary: Communication between viruses guides lysis-lysogeny decisions. Front. Microbiol. 2017, 8, 983. [Google Scholar] [CrossRef] [PubMed]
  89. Trinh, J.T.; Szekely, T.; Shao, Q.; Balazsi, G.; Zeng, L. Cell fate decisions emerge as phages cooperate or compete inside their host. Nat. Commun. 2017, 8, 14341. [Google Scholar] [CrossRef] [PubMed]
  90. Igler, C.; Abedon, S.T. Commentary: A host-produced quorum-sensing autoinducer controls a phage lysis-lysogeny decision. Front. Microbiol. 2019, 10, 1171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Silpe, J.E.; Bassler, B.L. A host-produced quorum-sensing autoinducer controls a phage lysis-lysogeny decision. Cell 2019, 176, 268–280. [Google Scholar] [CrossRef] [Green Version]
  92. Hynes, A.P.; Moineau, S. Phagebook: The social network. Mol. Cell 2017, 65, 963–964. [Google Scholar] [CrossRef] [Green Version]
  93. Lloyd, K.G.; Steen, A.D.; Ladau, J.; Yin, J.; Crosby, L. Phylogenetically novel uncultured microbial cells dominate Earth microbiomes. mSystems 2018, 3, e00055-18. [Google Scholar] [CrossRef] [Green Version]
  94. Chen, J.; Novick, R.P. Phage-mediated intergeneric transfer of toxin genes. Science 2009, 323, 139–141. [Google Scholar] [CrossRef] [Green Version]
  95. Willner, D.; Hugenholtz, P. From deep sequencing to viral tagging: Recent advances in viral metagenomics. BioEssays 2013, 35, 436–442. [Google Scholar] [CrossRef]
  96. Sullivan, M.B.; Waterbury, J.B.; Chisholm, S.W. Cyanophages infecting the oceanic cyanobacterium Prochlorococcus. Nature 2003, 424, 1047–1051. [Google Scholar] [CrossRef] [PubMed]
  97. Waterbury, J.B.; Valois, F.W. Resistance to co-occurring phages enables marine Synechococcus communities to coexist with cyanophages abundant in seawater. Appl. Environ. Microbiol. 1993, 59, 3393–3399. [Google Scholar] [CrossRef] [Green Version]
  98. Thomas, E.L.; Corbel, M.J. Isolation of a phage lytic for several Brucella species following propagation of Tbilisi phage in the presence of mitomycin C. Arch. Virol. 1977, 54, 259–261. [Google Scholar] [CrossRef] [PubMed]
  99. Henry, M.; Biswas, B.; Vincent, L.; Mokashi, V.; Schuch, R.; Bishop-Lilly, K.A.; Sozhamannan, S. Development of a high throughput assay for indirectly measuring phage growth using the OmniLog(TM) system. Bacteriophage 2012, 2, 159–167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  100. Howard-Varona, C.; Roux, S.; Dore, H.; Solonenko, N.E.; Holmfeldt, K.; Markillie, L.M.; Orr, G.; Sullivan, M.B. Regulation of infection efficiency in a globally abundant marine Bacteriodetes virus. ISME J. 2017, 11, 284–295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Enav, H.; Kirzner, S.; Lindell, D.; Mandel-Gutfreund, Y.; Beja, O. Adaptation to sub-optimal hosts is a driver of viral diversification in the ocean. Nat. Commun. 2018, 9, 4698. [Google Scholar] [CrossRef] [Green Version]
  102. Buckling, A.; Rainey, P.B. Antagonistic coevolution between a bacterium and a bacteriophage. Proc. R. Soc. Lond. B Biol. Sci. 2002, 269, 931–936. [Google Scholar] [CrossRef] [Green Version]
  103. Buckling, A.; Rainey, P.B. The role of parasites in sympatric and allopatric host diversification. Nature 2002, 420, 496–499. [Google Scholar] [CrossRef]
  104. Gomez, P.; Buckling, A. Real-time microbial adaptive diversification in soil. Ecol. Lett. 2013, 16, 650–655. [Google Scholar] [CrossRef]
  105. Poisot, T.; Lepennetier, G.; Martinez, E.; Ramsayer, J.; Hochberg, M.E. Resource availability affects the structure of a natural bacteria-bacteriophage community. Biol. Lett. 2011, 7, 201–204. [Google Scholar] [CrossRef] [Green Version]
  106. Vos, M.; Birkett, P.J.; Birch, E.; Griffiths, R.I.; Buckling, A. Local adaptation of bacteriophages to their bacterial hosts in soil. Science 2009, 325, 833. [Google Scholar] [CrossRef] [Green Version]
  107. Hanauer, D.I.; Graham, M.J.; Betancur, L.; Bobrownicki, A.; Cresawn, S.G.; Garlena, R.A.; Jacobs-Sera, D.; Kaufmann, N.; Pope, W.H.; Russell, D.A.; et al. An inclusive research education community (iREC): Impact of the SEA-PHAGES program on research outcomes and student learning. Proc. Natl. Acad. Sci. USA 2017, 114, 13531–13536. [Google Scholar] [CrossRef] [Green Version]
  108. Sea Phages 2020. Available online: https://www.seaphages.org/ (accessed on 18 April 2020).
  109. Hatfull, G.F. Bacteriophage discovery and genomics. In Bacteriophages: Biology, Technology, Therapy; Harper, D.R., Abedon, S.T., Burrowes, B.H., McConville, M., Eds.; Springer: Berlin/Heidelberg, Germany, 2017; pp. 1–13. [Google Scholar]
  110. Hatfull Lab 2020. Available online: http://www.hatfull.org/courses (accessed on 18 April 2020).
  111. Dedrick, R.M.; Guerrero-Bustamante, C.A.; Garlena, R.A.; Russell, D.A.; Ford, K.; Harris, K.; Gilmour, K.C.; Soothill, J.; Jacobs-Sera, D.; Schooley, R.T.; et al. Engineered bacteriophages for treatment of a patient with a disseminated drug-resistant Mycobacterium abscessus. Nat. Med. 2019, 25, 730–733. [Google Scholar] [CrossRef] [PubMed]
  112. Ghabrial, S.A.; Caston, J.R.; Jiang, D.; Nibert, M.L.; Suzuki, N. 50-plus years of fungal viruses. Virology 2015, 479–480, 356–368. [Google Scholar] [CrossRef] [Green Version]
  113. Liu, C.; Li, M.; Redda, E.T.; Mei, J.; Zhang, J.; Wu, B.; Jiang, X. A novel double-stranded RNA mycovirus isolated from Trichoderma harzianum. Virol. J. 2019, 16, 113. [Google Scholar] [CrossRef]
  114. Arjona-Lopez, J.M.; Telengech, P.; Jamal, A.; Hisano, S.; Kondo, H.; Yelin, M.D.; Arjona-Girona, I.; Kanematsu, S.; Lopez-Herrera, C.J.; Suzuki, N. Novel, diverse RNA viruses from Mediterranean isolates of the phytopathogenic fungus, Rosellinia necatrix: Insights into evolutionary biology of fungal viruses. Environ. Microbiol. 2018, 20, 1464–1483. [Google Scholar] [CrossRef] [PubMed]
  115. Liu, L.; Xie, J.; Cheng, J.; Fu, Y.; Li, G.; Yi, X.; Jiang, D. Fungal negative-stranded RNA virus that is related to bornaviruses and nyaviruses. Proc. Natl. Acad. Sci. USA 2014, 111, 12205–12210. [Google Scholar] [CrossRef] [Green Version]
  116. Handelsman, J. Metagenomics: Application of genomics to uncultured microorganisms. Microbiol. Mol. Biol. Rev. 2004, 68, 669–685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  117. Daniel, R. The metagenomics of soil. Nat. Rev. Microbiol. 2005, 3, 470–478. [Google Scholar] [CrossRef]
  118. Tringe, S.G.; von, M.C.; Kobayashi, A.; Salamov, A.A.; Chen, K.; Chang, H.W.; Podar, M.; Short, J.M.; Mathur, E.J.; Detter, J.C.; et al. Comparative metagenomics of microbial communities. Science 2005, 308, 554–557. [Google Scholar] [CrossRef] [Green Version]
  119. Edwards, R.A.; Rohwer, F. Viral metagenomics. Nat. Rev. Microbiol. 2005, 3, 504–510. [Google Scholar] [CrossRef]
  120. Roux, S. A viral ecogenomics framework to uncover the secrets of nature’s “Microbe whisperers”. mSystems 2019, 4, e00111-19. [Google Scholar] [CrossRef] [Green Version]
  121. Trevors, J.T. One gram of soil: A microbial biochemical gene library. Antonie van Leeuwenhoek J. Microbiol. 2010, 97, 99–106. [Google Scholar] [CrossRef] [PubMed]
  122. Delmont, T.O.; Robe, P.; Clark, I.; Simonet, P.; Vogel, T.M. Metagenomic comparison of direct and indirect soil DNA extraction approaches. J. Microbiol. Methods 2011, 86, 397–400. [Google Scholar] [CrossRef] [PubMed]
  123. Kunin, V.; Copeland, A.; Lapidus, A.; Mavromatis, K.; Hugenholtz, P. A bioinformatician’s guide to metagenomics. Microbiol. Mol. Biol. Rev. 2008, 72, 557–578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  124. Roux, S.; Emerson, J.B.; Eloe-Fadrosh, E.A.; Sullivan, M.B. Benchmarking viromics: An in silico evaluation of metagenome-enabled estimates of viral community composition and diversity. PeerJ 2017, 5, e3817. [Google Scholar] [CrossRef] [Green Version]
  125. Vestergaard, G.; Schulz, S.; Schöler, A.; Schloter, M. Making big data smart—how to use metagenomics to understand soil quality. Biol. Fertil. Soils 2017, 53, 479–484. [Google Scholar] [CrossRef]
  126. McLaren, M.R.; Willis, A.D.; Callahan, B.J. Consistent and correctable bias in metagenomic sequencing experiments. eLife 2019, 8, e46923. [Google Scholar] [CrossRef]
  127. Martinez-Hernandez, F.; Fornas, O.; Lluesma, G.M.; Bolduc, B.; de la Cruz Peña, M.; Martinez, J.M.; Anton, J.; Gasol, J.M.; Rosselli, R.; Rodriguez-Valera, F.; et al. Single-virus genomics reveals hidden cosmopolitan and abundant viruses. Nat. Commun. 2017, 8, 15892. [Google Scholar] [CrossRef] [Green Version]
  128. 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. 2019, 10, 1169. [Google Scholar] [CrossRef]
  129. Rédei, G.P. Encyclopedia of Genetics, Genomics, Proteomics, and Informatics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  130. Nomiya, T. Discussions on target theory: Past and present. J. Radiat. Res. 2013, 54, 1161–1163. [Google Scholar] [CrossRef] [Green Version]
  131. Sutton, T.D.S.; Clooney, A.G.; Ryan, F.J.; Ross, R.P.; Hill, C. Choice of assembly software has a critical impact on virome characterisation. Microbiome 2019, 7, 12. [Google Scholar] [CrossRef] [Green Version]
  132. Goordial, J.; Davila, A.; Greer, C.W.; Cannam, R.; DiRuggiero, J.; McKay, C.P.; Whyte, L.G. Comparative activity and functional ecology of permafrost soils and lithic niches in a hyper-arid polar desert. Environ. Microbiol. 2017, 19, 443–458. [Google Scholar] [CrossRef] [PubMed]
  133. Bowers, R.M.; Kyrpides, N.C.; Stepanauskas, R.; Harmon-Smith, M.; Doud, D.; Reddy, T.B.K.; Schulz, F.; Jarett, J.; Rivers, A.R.; Eloe-Fadrosh, E.A.; et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat. Biotech. 2017, 35, 725–731. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  134. Song, K.; Li, L.; Zhang, G. Coverage recommendation for genotyping analysis of highly heterologous species using next-generation sequencing technology. Sci. Rep. 2016, 6, 35736. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  135. Zaheer, R.; Noyes, N.; Ortega, P.R.; Cook, S.R.; Marinier, E.; Van, D.G.; Belk, K.E.; Morley, P.S.; McAllister, T.A. Impact of sequencing depth on the characterization of the microbiome and resistome. Sci. Rep. 2018, 8, 5890. [Google Scholar] [CrossRef] [PubMed]
  136. Roux, S.; Enault, F.; Hurwitz, B.L.; Sullivan, M.B. VirSorter: Mining viral signal from microbial genomic data. PeerJ 2015, 3, e985. [Google Scholar] [CrossRef]
  137. 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 2017, 5, 69. [Google Scholar] [CrossRef]
  138. Ren, J.; Song, K.; Deng, C.; Ahlgren, N.A.; Fuhrman, J.A.; Li, Y.; Xie, X.; Sun, F. Identifying viruses from metagenomic data by deep learning. arXiv 2018, arXiv:1806.07810. [Google Scholar] [CrossRef] [Green Version]
  139. CyVerse. 2020. Available online: https://www.cyverse.org/ (accessed on 18 April 2020).
  140. KBase Predictive Biology. 2020. Available online: https://www.kbase.us/ (accessed on 18 April 2020).
  141. Bolduc, B.; Youens-Clark, K.; Roux, S.; Hurwitz, B.L.; Sullivan, M.B. iVirus: Facilitating new insights in viral ecology with software and community data sets imbedded in a cyberinfrastructure. ISME J. 2017, 11, 7–14. [Google Scholar] [CrossRef] [Green Version]
  142. Benjamin, B. Welcome to iVirus’s Documentation. 2020. Available online: https://ivirus.readthedocs.io/en/latest/ (accessed on 18 April 2020).
  143. Hurwitz lab. VerveNet. 2020. Available online: https://www.protocols.io/groups/verve-net (accessed on 18 April 2020).
  144. Sangwan, N.; Xia, F.; Gilbert, J.A. Recovering complete and draft population genomes from metagenome datasets. Microbiome 2016, 4, 8. [Google Scholar] [CrossRef] [Green Version]
  145. Roux, S.; Adriaenssens, E.M.; Dutilh, B.E.; Koonin, E.V.; Kropinski, A.M.; Krupovic, M.; Kuhn, J.H.; Lavigne, R.; Brister, J.R.; Varsani, A.; et al. Minimum information about an uncultivated virus genome (MIUViG). Nat. Biotech. 2019, 37, 29–37. [Google Scholar] [CrossRef]
  146. Alrasheed, H.; Jin, R.; Weitz, J.S. Caution in inferring viral strategies from abundance correlations in marine metagenomes. Nat. Commun. 2019, 10, 501. [Google Scholar] [CrossRef] [PubMed]
  147. 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 2019, 116, 25900–25908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  148. de Farias, S.T.; Dos Santos Junior, A.P.; Rego, T.G.; Jose, M.V. Origin and evolution of RNA-dependent RNA polymerase. Front. Genet. 2017, 8, 125. [Google Scholar] [CrossRef]
  149. Ponsero, A.J.; Hurwitz, B.L. The promises and pitfalls of machine learning for detecting viruses in aquatic metagenomes. Front. Microbiol. 2019, 10, 806. [Google Scholar] [CrossRef]
  150. Breitbart, M.; Salamon, P.; Andresen, B.; Mahaffy, J.M.; Segall, A.M.; Mead, D.; Azam, F.; Rohwer, F. Genomic analysis of uncultured marine viral communities. Proc. Natl. Acad. Sci. USA 2002, 99, 14250–14255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  151. Gregory, A.C.; Zayed, A.A.; Conceicao-Neto, N.; Temperton, B.; Bolduc, B.; Alberti, A.; Ardyna, M.; Arkhipova, K.; Carmichael, M.; Cruaud, C.; et al. Marine DNA viral macro- and microdiversity from pole to pole. Cell 2019, 177, 1109–1123. [Google Scholar] [CrossRef] [Green Version]
  152. Michen, B.; Graule, T. Isoelectric points of viruses. J. Appl. Microbiol. 2010, 109, 388–397. [Google Scholar] [CrossRef] [Green Version]
  153. Zablocki, O.; van, Z.L.; Adriaenssens, E.M.; Rubagotti, E.; Tuffin, M.; Cary, C.; Cowan, D. High-level diversity of tailed phages, eukaryote-associated viruses, and virophage-like elements in the metaviromes of Antarctic soils. Appl. Environ. Microbiol. 2014, 80, 6888–6897. [Google Scholar] [CrossRef] [Green Version]
  154. Segobola, J.; Adriaenssens, E.; Tsekoa, T.; Rashamuse, K.; Cowan, D. Exploring viral diversity in a unique South African soil habitat. Sci. Rep. 2018, 8, 111. [Google Scholar] [CrossRef] [Green Version]
  155. Trubl, G.; Roux, S.; Solonenko, N.; Li, Y.F.; Bolduc, B.; Rodriguez-Ramos, J.; Eloe-Fadrosh, E.A.; Rich, V.I.; Sullivan, M.B. Towards optimized viral metagenomes for double-stranded and single-stranded DNA viruses from challenging soils. PeerJ 2019, 7, e7265. [Google Scholar] [CrossRef]
  156. Mehlich, A. Determination of cation-and anion-exchange properties of soils. Soil Sci. 1948, 66, 429–446. [Google Scholar] [CrossRef]
  157. Kozlowski, L.P. Proteome-pI: Proteome isoelectric point database. Nucl. Acids Res. 2017, 45, D1112–D1116. [Google Scholar] [CrossRef]
  158. Thurber, R.V.; Haynes, M.; Breitbart, M.; Wegley, L.; Rohwer, F. Laboratory procedures to generate viral metagenomes. Nat. Protoc. 2009, 4, 470–483. [Google Scholar] [CrossRef] [PubMed]
  159. Parras-Molto, M.; Rodriguez-Galet, A.; Suarez-Rodriguez, P.; Lopez-Bueno, A. Evaluation of bias induced by viral enrichment and random amplification protocols in metagenomic surveys of saliva DNA viruses. Microbiome 2018, 6, 119. [Google Scholar] [CrossRef] [PubMed]
  160. Anesio, A.M.; Hollas, C.; Graneli, W.; Laybourn-Parry, J. Influence of humic substances on bacterial and viral dynamics in freshwaters. Appl. Environ. Microbiol. 2004, 70, 4848–4854. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  161. Tebbe, C.C.; Vahjen, W. Interference of humic acids and DNA extracted directly from soil in detection and transformation of recombinant DNA from bacteria and a yeast. Appl. Environ. Microbiol. 1993, 59, 2657–2665. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  162. Alaeddini, R. Forensic implications of PCR inhibition—A review. Forensic Sci. Int. Genet. 2012, 6, 297–305. [Google Scholar] [CrossRef]
  163. Williamson, K.E.; Corzo, K.A.; Drissi, C.L.; Buckingham, J.M.; Thompson, C.P.; Helton, R.R. Estimates of viral abundance in soils are strongly influenced by extraction and enumeration methods. Biol. Fertil. Soils 2013, 49, 857–869. [Google Scholar] [CrossRef]
  164. Trubl, G.; Solonenko, N.; Chittick, L.; Solonenko, S.A.; Rich, V.I.; Sullivan, M.B. Optimization of viral resuspension methods for carbon-rich soils along a permafrost thaw gradient. PeerJ 2016, 4, e1999. [Google Scholar] [CrossRef] [Green Version]
  165. Patel, A.; Noble, R.T.; Steele, J.A.; Schwalbach, M.S.; Hewson, I.; Fuhrman, J.A. Virus and prokaryote enumeration from planktonic aquatic environments by epifluorescence microscopy with SYBR Green I. Nat. Protoc. 2007, 2, 269–276. [Google Scholar] [CrossRef] [Green Version]
  166. Forterre, P.; Soler, N.; Krupovic, M.; Marguet, E.; Ackermann, H.-W. Fake virus particles generated by fluorescence microscopy. Trends Microbiol. 2013, 21, 1–5. [Google Scholar] [PubMed]
  167. Cunningham, B.R.; Brum, J.R.; Schwenck, S.M.; Sullivan, M.B.; John, S.G. An inexpensive, accurate, and precise wet-mount method for enumerating aquatic viruses. Appl. Environ. Microbiol. 2015, 81, 2995–3000. [Google Scholar] [PubMed] [Green Version]
  168. Zhang, Y.; Hung, T.; Song, J.; He, J. Electron microscopy: Essentials for viral structure, morphogenesis and rapid diagnosis. Sci. China Life Sci. 2013, 56, 421–430. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  169. Burge, W.D.; Enkiri, N.K. Virus adsorption by five soils. J. Environ. Qual. 1978, 7, 73–76. [Google Scholar] [CrossRef]
  170. Rodriguez, R.A.; Pepper, I.L.; Gerba, C.P. Application of PCR-based methods to assess the infectivity of enteric viruses in environmental samples. Appl. Environ. Microbiol. 2009, 75, 297–307. [Google Scholar] [CrossRef] [Green Version]
  171. Farkas, K.; Hassard, F.; McDonald, J.E.; Malham, S.K.; Jones, D.L. Evaluation of molecular methods for the detection and quantification of pathogen-derived nucleic acids in sediment. Front. Microbiol. 2017, 8, 53. [Google Scholar] [CrossRef] [Green Version]
  172. Krishnamurthy, S.R.; Wang, D. Origins and challenges of viral dark matter. Virus Res. 2017, 239, 136–142. [Google Scholar] [CrossRef]
  173. Pruitt, K.D.; Tatusova, T.; Klimke, W.; Maglott, D.R. NCBI reference sequences: Current status, policy and new initiatives. Nucl. Acids Res. 2009, 37, D32–D36. [Google Scholar] [CrossRef] [Green Version]
  174. Paez-Espino, D.; Chen, I.A.; Palaniappan, K.; Ratner, A.; Chu, K.; Szeto, E.; Pillay, M.; Huang, J.; Markowitz, V.M.; Nielsen, T.; et al. IMG/VR: A database of cultured and uncultured DNA Viruses and retroviruses. Nucl. Acids Res. 2017, 45, D457–D465. [Google Scholar] [CrossRef] [PubMed]
  175. Hurwitz, B.L.; U’Ren, J.M.; Youens-Clark, K. Computational prospecting the great viral unknown. FEMS Microbiol. Lett. 2016, 363. [Google Scholar] [CrossRef] [Green Version]
  176. Brum, J.R.; Ignacio-Espinoza, J.C.; Roux, S.; Doulcier, G.; Acinas, S.G.; Alberti, A.; Chaffron, S.; Cruaud, C.; de, V.C.; Gasol, J.M.; et al. Ocean plankton. Patterns and ecological drivers of ocean viral communities. Science 2015, 348, 1261498. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  177. Roux, S.; Brum, J.R.; Dutilh, B.E.; Sunagawa, S.; Duhaime, M.B.; Loy, A.; Poulos, B.T.; Solonenko, N.; Lara, E.; Poulain, J.; et al. Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses. Nature 2016, 537, 689–693. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  178. Rose, R.; Constantinides, B.; Tapinos, A.; Robertson, D.L.; Prosperi, M. Challenges in the analysis of viral metagenomes. Virus Evol. 2016, 2, vew022. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  179. Brooke, C.B. Biological activities of ‘noninfectious’ influenza A virus particles. Future Virol. 2014, 9, 41–51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  180. Wilpiszeski, R.L.; Aufrecht, J.A.; Retterer, S.T.; Sullivan, M.B.; Graham, D.E.; Pierce, E.M.; Zablocki, O.D.; Palumbo, A.V.; Elias, D.A. Soil aggregate microbial communities: Towards understanding microbiome interactions at biologically relevant scales. Appl. Environ. Microbiol. 2019, 85, e00324-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  181. Nagler, M.; Insam, H.; Pietramellara, G.; Ascher-Jenull, J. Extracellular DNA in natural environments: Features, relevance and applications. Appl. Microbiol. Biotechnol. 2018, 102, 6343–6356. [Google Scholar] [CrossRef] [Green Version]
  182. Carini, P.; Marsden, P.J.; Leff, J.W.; Morgan, E.E.; Strickland, M.S.; Fierer, N. Relic DNA is abundant in soil and obscures estimates of soil microbial diversity. Nat. Microbiol. 2016, 2, 16242. [Google Scholar] [CrossRef]
  183. Hurwitz, B.L.; Deng, L.; Poulos, B.T.; Sullivan, M.B. Evaluation of methods to concentrate and purify ocean virus communities through comparative, replicated metagenomics. Environ. Microbiol. 2013, 15, 1428–1440. [Google Scholar] [CrossRef] [Green Version]
  184. Zolfo, M.; Pinto, F.; Asnicar, F.; Manghi, P.; Tett, A.; Bushman, F.D.; Segata, N. Detecting contamination in viromes using ViromeQC. Nat. Biotech. 2019, 37, 1408–1412. [Google Scholar] [CrossRef]
  185. Fox, K.R. DNase I footprinting. In Drug-DNA Interaction Protocols; Springer: Berlin/Heidelberg, Germany, 1997; pp. 1–22. [Google Scholar]
  186. Romanowski, G.; Lorenz, M.G.; Wackernagel, W. Adsorption of plasmid DNA to mineral surfaces and protection against DNase I. Appl. Environ. Microbiol. 1991, 57, 1057–1061. [Google Scholar] [CrossRef] [Green Version]
  187. Crecchio, C.; Stotzky, G. Binding of DNA on humic acids: Effect on transformation of Bacillus subtilis and resistance to DNase. Soil Biol. Biochem. 1998, 30, 1061–1067. [Google Scholar] [CrossRef]
  188. Cai, P.; Huang, Q.Y.; Zhang, X.W. Interactions of DNA with clay minerals and soil colloidal particles and protection against degradation by DNase. Environ. Sci. Technol. 2006, 40, 2971–2976. [Google Scholar] [CrossRef]
  189. Emerson, J.B.; Adams, R.I.; Roman, C.M.B.; Brooks, B.; Coil, D.A.; Dahlhausen, K.; Ganz, H.H.; Hartmann, E.M.; Hsu, T.; Justice, N.B.; et al. Schrödinger’s microbes: Tools for distinguishing the living from the dead in microbial ecosystems. Microbiome 2017, 5, 86. [Google Scholar] [CrossRef]
  190. Lennon, J.T.; Muscarella, M.E.; Placella, S.A.; Lehmkuhl, B.K. How, when, and where relic DNA affects microbial diversity. MBio 2018, 9, e00637-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  191. Fittipaldi, M.; Rodriguez, N.J.; Codony, F.; Adrados, B.; Penuela, G.A.; Morato, J. Discrimination of infectious bacteriophage T4 virus by propidium monoazide real-time PCR. J. Virol. Methods 2010, 168, 228–232. [Google Scholar] [CrossRef] [PubMed]
  192. Biotium. PMA & PMAxx™ Selected References. 2020. Available online: https://biotium.com/wp-content/uploads/2017/10/PMA-PMAxx-References.pdf (accessed on 18 April 2020).
  193. Mizuno, C.M.; Guyomar, C.; Roux, S.; Lavigne, R.; Rodriguez-Valera, F.; Sullivan, M.B.; Gillet, R.; Forterre, P.; Krupovic, M. Numerous cultivated and uncultivated viruses encode ribosomal proteins. Nat. Commun. 2019, 10, 752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  194. Yuan, Y.; Gao, M. Jumbo bacteriophages: An overview. Front. Microbiol. 2017, 8, 403. [Google Scholar] [CrossRef] [Green Version]
  195. Kuhn, E.; Ichimura, A.S.; Peng, V.; Fritsen, C.H.; Trubl, G.; Doran, P.T.; Murray, A.E. Brine assemblages of ultrasmall microbial cells within the ice cover of Lake Vida, Antarctica. Appl. Environ. Microbiol. 2014, 80, 3687–3698. [Google Scholar] [CrossRef] [Green Version]
  196. Velimirov, B. Nanobacteria, ultramicrobacteria and starvation forms: A search for the smallest metabolizing bacterium. Microbes Environ. 2001, 16, 67–77. [Google Scholar] [CrossRef] [Green Version]
  197. Ghuneim, L.J.; Jones, D.L.; Golyshin, P.N.; Golyshina, O.V. Nano-sized and filterable Bacteria and Archaea: Biodiversity and function. Front. Microbiol. 2018, 9, 1971. [Google Scholar] [CrossRef] [Green Version]
  198. Kristensen, D.M.; Mushegian, A.R.; Dolja, V.V.; Koonin, E.V. New dimensions of the virus world discovered through metagenomics. Trends Microbiol. 2010, 18, 11–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  199. Hurwitz, B.L.; Sullivan, M.B. The Pacific Ocean virome (POV): A marine viral metagenomic dataset and associated protein clusters for quantitative viral ecology. PLoS ONE 2013, 8, e57355. [Google Scholar] [CrossRef] [PubMed]
  200. Lang, A.S.; Westbye, A.B.; Beatty, J.T. The distribution, evolution, and roles of gene transfer agents in prokaryotic genetic exchange. Annu. Rev. Virol. 2017, 4, 87–104. [Google Scholar] [CrossRef] [PubMed]
  201. Shakya, M.; Soucy, S.M.; Zhaxybayeva, O. Insights into origin and evolution of alpha-proteobacterial gene transfer agents. Virus Evol. 2017, 3, vex036. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  202. Delgado-Baquerizo, M.; Oliverio, A.M.; Brewer, T.E.; Benavent-Gonzalez, A.; Eldridge, D.J.; Bardgett, R.D.; Maestre, F.T.; Singh, B.K.; Fierer, N. A global atlas of the dominant bacteria found in soil. Science 2018, 359, 320–325. [Google Scholar] [CrossRef] [Green Version]
  203. Shintani, M.; Sanchez, Z.K.; Kimbara, K. Genomics of microbial plasmids: Classification and identification based on replication and transfer systems and host taxonomy. Front. Microbiol. 2015, 6, 242. [Google Scholar] [CrossRef] [PubMed]
  204. Ruhfel, R.E.; Robillard, N.J.; Thorne, C.B. Interspecies transduction of plasmids among Bacillus anthracis, B. cereus, and B. thuringiensis. J. Bacteriol. 1984, 157, 708–711. [Google Scholar] [CrossRef] [Green Version]
  205. Koonin, E.V.; Dolja, V.V. Virus world as an evolutionary network of viruses and capsidless selfish elements. Microbiol. Mol. Biol. Rev. 2014, 78, 278–303. [Google Scholar] [CrossRef] [Green Version]
  206. Kazlauskas, D.; Varsani, A.; Koonin, E.V.; Krupovic, M. Multiple origins of prokaryotic and eukaryotic single-stranded DNA viruses from bacterial and archaeal plasmids. Nat. Commun. 2019, 10, 3425. [Google Scholar] [CrossRef] [Green Version]
  207. Roux, S.; Krupovic, M.; Debroas, D.; Forterre, P.; Enault, F. Assessment of viral community functional potential from viral metagenomes may be hampered by contamination with cellular sequences. Open. Biol. 2013, 3, 130160. [Google Scholar] [CrossRef] [Green Version]
  208. Rozov, R.; Brown, K.A.; Bogumil, D.; Shterzer, N.; Halperin, E.; Mizrahi, I.; Shamir, R. Recycler: An algorithm for detecting plasmids from de novo assembly graphs. Bioinformatics 2017, 33, 475–482. [Google Scholar] [PubMed] [Green Version]
  209. Krawczyk, P.S.; Lipinski, L.; Dziembowski, A. PlasFlow: Predicting plasmid sequences in metagenomic data using genome signatures. Nucl. Acids Res. 2018, 46, e35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  210. Beaulaurier, J.; Zhu, S.; Deikus, G.; Mogno, I.; Zhang, X.S.; Davis-Richardson, A.; Canepa, R.; Triplett, E.W.; Faith, J.J.; Sebra, R.; et al. Metagenomic binning and association of plasmids with bacterial host genomes using DNA methylation. Nat. Biotech. 2018, 36, 61–69. [Google Scholar] [CrossRef] [PubMed]
  211. Shulman, L.M.; Davidson, I. Viruses with circular single-stranded DNA genomes are everywhere! Annu. Rev. Virol. 2017, 4, 159–180. [Google Scholar] [CrossRef] [PubMed]
  212. Malathi, V.G.; Devi, P.R. ssDNA viruses: Key players in global virome. Virusdisease 2019, 30, 3–12. [Google Scholar] [CrossRef]
  213. Karlsson, O.E.; Belak, S.; Granberg, F. The effect of preprocessing by sequence-independent, single-primer amplification (SISPA) on metagenomic detection of viruses. Biosecur. Bioterror. 2013, 11 (Suppl. 1), S227–S234. [Google Scholar] [CrossRef]
  214. Roux, S.; Solonenko, N.E.; Dang, V.T.; Poulos, B.T.; Schwenck, S.M.; Goldsmith, D.B.; Coleman, M.L.; Breitbart, M.; Sullivan, M.B. Towards quantitative viromics for both double-stranded and single-stranded DNA viruses. PeerJ 2016, 4, e2777. [Google Scholar] [CrossRef] [Green Version]
  215. Reuter, J.A.; Spacek, D.V.; Snyder, M.P. High-throughput sequencing technologies. Mol. Cell 2015, 58, 586–597. [Google Scholar] [CrossRef] [Green Version]
  216. Gansauge, M.T.; Meyer, M. Single-stranded DNA library preparation for the sequencing of ancient or damaged DNA. Nat. Protoc. 2013, 8, 737–748. [Google Scholar] [CrossRef]
  217. Aigrain, L.; Gu, Y.; Quail, M.A. Quantitation of next generation sequencing library preparation protocol efficiencies using droplet digital PCR assays - a systematic comparison of DNA library preparation kits for Illumina sequencing. BMC Genom. 2016, 17, 458. [Google Scholar] [CrossRef] [Green Version]
  218. Bekliz, M.; Brandani, J.; Bourquin, M.; Battin, T.J.; Peter, H. Benchmarking protocols for the metagenomic analysis of stream biofilm viromes. PeerJ 2019, 7, e8187. [Google Scholar] [CrossRef] [PubMed]
  219. Reavy, B.; Swanson, M.M.; Cock, P.J.; Dawson, L.; Freitag, T.E.; Singh, B.K.; Torrance, L.; Mushegian, A.R.; Taliansky, M. Distinct circular single-stranded DNA viruses exist in different soil types. Appl. Environ. Microbiol. 2015, 81, 3934–3945. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  220. Breitbart, M.; Rohwer, F. Here a virus, there a virus, everywhere the same virus? Trends Microbiol. 2005, 13, 278–284. [Google Scholar] [CrossRef] [PubMed]
  221. Haig, S.J.; Schirmer, M.; D’Amore, R.; Gibbs, J.; Davies, R.L.; Collins, G.; Quince, C. Stable-isotope probing and metagenomics reveal predation by protozoa drives E. coli removal in slow sand filters. ISME J. 2015, 9, 797–808. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Different categories of sources of environmental viral genomic DNA. Free virions (1) initiate latent infections via cell adsorption (2) or instead initiate productive infections (3), where the latter can be differentiated into lytic vs. chronically releasing virions (not distinguished in the figure). Virus-like eDNA (vleDNA) is a form of extracellular, unencapsidated DNA (4). Viral infections also can take on various forms described here as ‘Other’ (5).
Figure 1. Different categories of sources of environmental viral genomic DNA. Free virions (1) initiate latent infections via cell adsorption (2) or instead initiate productive infections (3), where the latter can be differentiated into lytic vs. chronically releasing virions (not distinguished in the figure). Virus-like eDNA (vleDNA) is a form of extracellular, unencapsidated DNA (4). Viral infections also can take on various forms described here as ‘Other’ (5).
Soilsystems 04 00023 g001
Figure 2. Overview of three common methods of virus characterization. Major methodological steps using virus-isolation (green), metagenomic (blue), and viromic (purple) approaches are shown. Possible contaminants, which are entities that may be described as viruses but which are not actually viruses, are denoted with icons for each approach (as summarized in Section 4). The icons are listed from left to right in order of potential prevalence in each method, although the order will depend on the sample and how it is treated. Finally, the pros and cons are listed for each approach. The pros for viromes are as relative to metagenomes.
Figure 2. Overview of three common methods of virus characterization. Major methodological steps using virus-isolation (green), metagenomic (blue), and viromic (purple) approaches are shown. Possible contaminants, which are entities that may be described as viruses but which are not actually viruses, are denoted with icons for each approach (as summarized in Section 4). The icons are listed from left to right in order of potential prevalence in each method, although the order will depend on the sample and how it is treated. Finally, the pros and cons are listed for each approach. The pros for viromes are as relative to metagenomes.
Soilsystems 04 00023 g002
Figure 3. Illustration of horizontal coverage vs. vertical coverage from sequencing reads, and the impact of increased vertical coverage on horizontal coverage. (A) A hypothetical genome is shown as a double arrow (blue) with sequencing breadth indicated in terms of horizontal coverage of sequencing reads (green bars). (B) Sequencing depth by contrast is indicated in terms of vertical coverage (stacking) of overlapping sequencing reads (also of green bars). (C) Taking the sequencing reads providing increased vertical coverage (from B) and collapsing them into a single layer (bottom) illustrates the potential for greater horizontal coverage resulting from increased vertical coverage compared to decreased vertical coverage. Note that increased vertical coverage can also increase sequencing accuracy in terms of defining consensus sequence, though this potential increase in sequencing accuracy is not indicated in the figure.
Figure 3. Illustration of horizontal coverage vs. vertical coverage from sequencing reads, and the impact of increased vertical coverage on horizontal coverage. (A) A hypothetical genome is shown as a double arrow (blue) with sequencing breadth indicated in terms of horizontal coverage of sequencing reads (green bars). (B) Sequencing depth by contrast is indicated in terms of vertical coverage (stacking) of overlapping sequencing reads (also of green bars). (C) Taking the sequencing reads providing increased vertical coverage (from B) and collapsing them into a single layer (bottom) illustrates the potential for greater horizontal coverage resulting from increased vertical coverage compared to decreased vertical coverage. Note that increased vertical coverage can also increase sequencing accuracy in terms of defining consensus sequence, though this potential increase in sequencing accuracy is not indicated in the figure.
Soilsystems 04 00023 g003

Share and Cite

MDPI and ACS Style

Trubl, G.; Hyman, P.; Roux, S.; Abedon, S.T. Coming-of-Age Characterization of Soil Viruses: A User’s Guide to Virus Isolation, Detection within Metagenomes, and Viromics. Soil Syst. 2020, 4, 23. https://0-doi-org.brum.beds.ac.uk/10.3390/soilsystems4020023

AMA Style

Trubl G, Hyman P, Roux S, Abedon ST. Coming-of-Age Characterization of Soil Viruses: A User’s Guide to Virus Isolation, Detection within Metagenomes, and Viromics. Soil Systems. 2020; 4(2):23. https://0-doi-org.brum.beds.ac.uk/10.3390/soilsystems4020023

Chicago/Turabian Style

Trubl, Gareth, Paul Hyman, Simon Roux, and Stephen T. Abedon. 2020. "Coming-of-Age Characterization of Soil Viruses: A User’s Guide to Virus Isolation, Detection within Metagenomes, and Viromics" Soil Systems 4, no. 2: 23. https://0-doi-org.brum.beds.ac.uk/10.3390/soilsystems4020023

Article Metrics

Back to TopTop