Membrane vesicles from antibiotic-resistant Staphylococcus aureus transfer antibiotic-resistance to antibiotic-susceptible Escherichia coli.

AIM: Bacteria naturally produce membrane vesicles (MVs), which have been shown to contribute to the spread of multi-drug resistant bacteria (MDR) by delivering antibiotic-resistant substances to antibiotic-susceptible bacteria. Here, we aim to show that MVs from Gram-positive bacteria are capable of transferring ß-lactam antibiotic-resistant substances to antibiotic-sensitive Gram-negative bacteria. MATERIALS AND METHODS: MVs were collected from a methicillin-resistant strain of Staphylococcus aureus (MRSA) and vesicle-mediated fusion with antimicrobial-sensitive Escherichia coli (RC85). It was performed by exposing the bacteria to the MVs to develop antimicrobial-resistant E. coli (RC85-T). RESULTS: The RC85-T exhibited a higher resistance to ß-lactam antibiotics compared to the parent strain. Although the secretion rates of the MVs from RC85-T and the parent strain were nearly equal, the ß-lactamase activity of the MVs from RC85-T was 12-times higher than that of MVs from the parent strain, based on equivalent protein concentrations. Moreover, MVs secreted by RC85-T were able to protect ß-lactam-susceptible E. coli from ß-lactam antibiotic-induced growth inhibition in a dose-dependent manner. CONCLUSION: MVs play a role in transferring substances from Gram-positive to Gram-negative bacteria, shown by the release of MVs from RC85-T that were able to protect ß-lactam-susceptible bacteria from ß-lactam antibiotics. SIGNIFICANCE AND IMPACT OF STUDY: MVs are involved in the emergence of antibiotic-resistant strains in a mixed bacterial culture, helping us to understand how the spread of multidrug-resistant bacteria could be reduced.

Analysis of Spounaviruses as a Case Study for the Overdue Reclassification of Tailed Phages.

Tailed bacteriophages are the most abundant and diverse viruses in the world, with genome sizes ranging from 10 kbp to over 500 kbp. Yet, due to historical reasons, all this diversity is confined to a single virus order-Caudovirales, composed of just four families: Myoviridae, Siphoviridae, Podoviridae, and the newly created Ackermannviridae family. In recent years, this morphology-based classification scheme has started to crumble under the constant flood of phage sequences, revealing that tailed phages are even more genetically diverse than once thought. This prompted us, the Bacterial and Archaeal Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV), to consider overall reorganization of phage taxonomy. In this study, we used a wide range of complementary methods-including comparative genomics, core genome analysis, and marker gene phylogenetics-to show that the group of Bacillus phage SPO1-related viruses previously classified into the Spounavirinae subfamily, is clearly distinct from other members of the family Myoviridae and its diversity deserves the rank of an autonomous family. Thus, we removed this group from the Myoviridae family and created the family Herelleviridae-a new taxon of the same rank. In the process of the taxon evaluation, we explored the feasibility of different demarcation criteria and critically evaluated the usefulness of our methods for phage classification. The convergence of results, drawing a consistent and comprehensive picture of a new family with associated subfamilies, regardless of method, demonstrates that the tools applied here are particularly useful in phage taxonomy. We are convinced that creation of this novel family is a crucial milestone toward much-needed reclassification in the Caudovirales order.

Integrative omics analysis of Pseudomonas aeruginosa virus PA5oct highlights the molecular complexity of jumbo phages.

Pseudomonas virus vB_PaeM_PA5oct is proposed as a model jumbo bacteriophage to investigate phage-bacteria interactions and is a candidate for phage therapy applications. Combining hybrid sequencing, RNA-Seq and mass spectrometry allowed us to accurately annotate its 286,783 bp genome with 461 coding regions including four non-coding RNAs (ncRNAs) and 93 virion-associated proteins. PA5oct relies on the host RNA polymerase for the infection cycle and RNA-Seq revealed a gradual take-over of the total cell transcriptome from 21% in early infection to 93% in late infection. PA5oct is not organized into strictly contiguous regions of temporal transcription, but some genomic regions transcribed in early, middle and late phases of infection can be discriminated. Interestingly, we observe regions showing limited transcription activity throughout the infection cycle. We show that PA5oct upregulates specific bacterial operons during infection including operons pncA-pncB1-nadE involved in NAD biosynthesis, psl for exopolysaccharide biosynthesis and nap for periplasmic nitrate reductase production. We also observe a downregulation of T4P gene products suggesting mechanisms of superinfection exclusion. We used the proteome of PA5oct to position our isolate amongst other phages using a gene-sharing network. This integrative omics study illustrates the molecular diversity of jumbo viruses and raises new questions towards cellular regulation and phage-encoded hijacking mechanisms.

Taxonomy of prokaryotic viruses: 2018-2019 update from the ICTV Bacterial and Archaeal Viruses Subcommittee.

This article is a summary of the activities of the ICTV's Bacterial and Archaeal Viruses Subcommittee for the years 2018 and 2019. Highlights include the creation of a new order, 10 families, 22 subfamilies, 424 genera and 964 species. Some of our concerns about the ICTV's ability to adjust to and incorporate new DNA- and protein-based taxonomic tools are discussed.

Closely related viruses of the marine picoeukaryotic alga Ostreococcus lucimarinus exhibit different ecological strategies.

In marine ecosystems, viruses are major disrupters of the direct flow of carbon and nutrients to higher trophic levels. Although the genetic diversity of several eukaryotic phytoplankton virus groups has been characterized, their infection dynamics are less understood, such that the physiological and ecological implications of their diversity remain unclear. We compared genomes and infection phenotypes of the two most closely related cultured phycodnaviruses infecting the widespread picoprasinophyte Ostreococcus lucimarinus under standard- (1.3 divisions per day) and limited-light (0.41 divisions per day) nutrient replete conditions. OlV7 infection caused early arrest of the host cell cycle, coinciding with a significantly higher proportion of infected cells than OlV1-amended treatments, regardless of host growth rate. OlV7 treatments showed a near-50-fold increase of progeny virions at the higher host growth rate, contrasting with OlV1's 16-fold increase. However, production of OlV7 virions was more sensitive than OlV1 production to reduced host growth rate, suggesting fitness trade-offs between infection efficiency and resilience to host physiology. Moreover, although organic matter released from OlV1- and OlV7-infected hosts had broadly similar chemical composition, some distinct molecular signatures were observed. Collectively, these results suggest that current views on viral relatedness through marker and core gene analyses underplay operational divergence and consequences for host ecology.

Characterization of a specific monoclonal antibody against immunoglobulin light kappa/L1 chain in olive flounder (Paralichthys olivaceus).

Immunoglobulins (Ig) are heterodimeric proteins that play critical roles in the adaptive immune system of vertebrates. Because of their plasticity, teleostean Igs are more diverse, and thus do not conform to mammalian classifications. Because of this, mammalian-based Ig cell markers cannot be used successfully to study immune responses in fish. There is therefore a need to produce Ig-specific cell markers for fish. Here, we attempted to identify the specific isotype detected by an Ig light chain-specific monoclonal antibody (anti-olive flounder IgL-mAb: M7C3-4) that we had previously produced [11]. Three newly identified sequences of the Ig light chain from olive flounder were classified according to their isotypes. Subsequent analyses revealed that M7C3-4 was able to specifically detect lymphocytes expressing one of the κ chains (Igκ-a) in olive flounder. Interestingly, Igκ-a+ B cells were more abundant in spleen and trunk-kidney than in peripheral blood, indicating a distribution different from that of IgM+ B cells. Our work reveals interesting aspects of B cell distribution and differentiation, and may aid in the production of suitable and effective cell markers for olive flounder.

Phylogenomic network and comparative genomics reveal a diverged member of the ΦKZ-related group, marine vibrio phage ΦJM-2012.

Bacteriophages are the largest reservoir of genetic diversity. Here we describe the novel phage ΦJM-2012. This natural isolate from marine Vibrio cyclitrophicus possesses very few gene contents relevant to other well-studied marine Vibrio phages. To better understand its evolutionary history, we built a mathematical model of pairwise relationships among 1,221 phage genomes, in which the genomes (nodes) are linked by edges representing the normalized number of shared orthologous protein families. This weighted network revealed that ΦJM-2012 was connected to only five members of the Pseudomonas ΦKZ-like phage family in an isolated network, strongly indicating that it belongs to this phage group. However, comparative genomic analyses highlighted an almost complete loss of colinearity with the ΦKZ-related genomes and little conservation of gene order, probably reflecting the action of distinct evolutionary forces on the genome of ΦJM-2012. In this phage, typical conserved core genes, including six RNA polymerase genes, were frequently displaced and the hyperplastic regions were rich in both unique genes and predicted unidirectional promoters with highly correlated orientations. Further, analysis of the ΦJM-2012 genome showed that segments of the conserved N-terminal parts of ΦKZ tail fiber paralogs exhibited evidence of combinatorial assortment, having switched transcriptional orientation, and there was recruitment and/or structural changes among phage endolysins and tail spike protein. Thus, this naturally occurring phage appears to have branched from a common ancestor of the ΦKZ-related groups, showing a distinct genomic architecture and unique genes that most likely reflect adaptation to its chosen host and environment.

Combination treatment against scuticociliatosis by reducing the inhibitor effect of mucus in olive flounder, Paralichthys olivaceus.

The olive flounder, Paralichthys olivaceus, is an economically important food fish in Japan and Korea. Scuticociliatosis is a major parasitic disease, and fatal infection with scuticociliates, or mixed infections with scuticociliates and other pathogenic agents (e.g., Vibrio spp.) cause severe mortalities in farmed olive flounders. To date, however, effective chemotherapeutic treatment of scuticociliatosis has only been reported at the in vitro level. In this study, we employed combination treatment, using benzalkonium chloride (to remove excess mucus from the body surface) and bronopol (to kill the parasites), to overcome the protective effect of mucus by some medicine to the scuticociliates. In the presence of the mucus mixture, the higher dose of bronopol (156 ppm) yielded morphologies and motilities similar to those of ciliates treated with the lower dose of bronopol (80 ppm) in the absence of mucus. We also investigated the in vivo effects of this treatment in field trials involving a total of 15,025 naturally infected flounders. We observed that short-term bath treatments with benzalkonium chloride (50 ppm) followed by bronopol (500 ppm) were effective, assessed by the relative percentage mortality (RPS) value. Thus, this study provides a notable therapeutic strategy by removing the mucus to treat scuticociliatosis in olive flounders at the aquaculture field level.

Comparative sequence analysis of a multidrug-resistant plasmid from Aeromonas hydrophila.

Aeromonas hydrophila is a pathogenic bacterium that has been implicated in fish, animal, and human disease. Recently, a multidrug resistance (MDR) plasmid, pR148, was isolated from A. hydrophila obtained from a tilapia (Oreochromis niloticus) farm in Thailand. pR148 is a 165,906-bp circular plasmid containing 147 coding regions showing highest similarity to pNDM-1_Dok1, an MDR plasmid isolated from a human pathogen. pR148 was also very similar to other IncA/C plasmids isolated from humans, animals, food, and fish. pR148 contains a mercuric resistance operon and encodes the complete set of genes for the type 4 secretion system. pR148 encodes a Tn21 type transposon. This transposon contains the drug resistance genes qacH, bla(OXA-10), aadA1, and sul1 in a class 1 integron; tetA and tetR in transposon Tn1721; and catA2 and a duplicate sul1 in a locus showing 100% similarity to IncU plasmids isolated from fish. The bla(OXA-10) and aadA1 genes showed 100% similarity to those from the Acinetobacter baumannii AYE genome. The similarity of pR148 to a human pathogen-derived plasmid indicates that the plasmids were either transferred between different genera or that they are derived from a common origin. Previous studies have shown that IncA/C plasmids retain a conserved backbone, while the accessory region points to lateral gene transfer. These observations point out the dangers of indiscriminate use of antibiotics in humans and in animals and the necessity of understanding how drug resistance determinants are disseminated and transferred.

Complete genome sequence of the bacteriophages ECBP1 and ECBP2 isolated from two different Escherichia coli strains.

Escherichia coli is recognized as one of the most abundant avian bacterial pathogens. In this study, we report the sequencing by the traditional Sanger method of ECBP1 and ECBP2: bacteriophages that infected two different E. coli strains which might be used as therapeutic agents in combination with alternative antibiotics.

Heat shock protein profiles on the protein and gene expression levels in olive flounder kidney infected with Streptococcus parauberis.

Heat shock proteins (HSPs) have been observed in cells exposed to a variety of stresses, including infectious pathogens. This study used a label-free, quantitative proteomic approach and transcriptional gene expression analysis to investigate infection-related HSP proteins and their encoding genes in whole kidneys from olive flounder (Paralichthys olivaceus). During Streptococcus parauberis infection in the flounder, the genes encoding Hsp10, Hsp40A4, Hsp40B6, Hsp40B11, Hsp60, Hsp70, glucose regulated protein 78 (Grp78), Hsp90α, Hsp90ß and Grp94 were induced, and the protein levels of Hsp60, Hsp70, Hsp90α, Hsp90ß and Grp94 were differentially regulated over time. Subsequent results also revealed that Hsp60, Hsp70, Hsp90α, Hsp90ß and Grp94 appear to be the dominant and critical HSPs in olive flounder during bacterial infection. This is the first estimation of the differential involvement of HSPs in the immune response of olive flounder exposed to bacterial infection.

Recombinant interferon-γ activates immune responses against Edwardsiella tarda infection in the olive flounder, Paralichthys olivaceus.

Interferon gamma (IFN-γ) is a cytokine that plays a very important role in defining Th1 immune response in all vertebrates. In this study, recombinant IFN-γ (rIFN-γ) from the olive flounder (Paralichthys olivaceus) was produced in an Escherichia coli system using a pET expression vector. Stimulation of whole kidney leukocytes (immune-related cells) in vitro with the resulting rIFN-γ significantly induced the gene expression of interleukin-1ß (IL-1ß), signal transducer and activator of transcription 1 (STAT1), CXCL13-like chemokine (CXCL13), and IFN-γ. rIFN-γ also weakly induced the expression of IL-1ß, tumor necrosis factor-α (TNF-α), CXCL13, and IFN-γ in olive flounder-derived HINAE (non-immune) cells. The effects of rIFN-γ against Edwardsiella tarda infection in vivo were assessed by intraperitoneally injecting a mixture of rIFN-γ (100 ng) and E. tarda (1 × 10(5) CFU/ml) into the olive flounder. The survival rate in the rIFN-γ-injected group was 60% compared to 0% in the group treated with E. tarda only, demonstrating that olive flounder IFN-γ is effective in reinforcing immune responses and preventing against edwardsiellosis.

Viral tag and grow: a scalable approach to capture and characterize infectious virus-host pairs.

Viral metagenomics (viromics) has reshaped our understanding of DNA viral diversity, ecology, and evolution across Earth's ecosystems. However, viromics now needs approaches to link newly discovered viruses to their host cells and characterize them at scale. This study adapts one such method, sequencing-enabled viral tagging (VT), to establish "Viral Tag and Grow" (VT + Grow) to rapidly capture and characterize viruses that infect a cultivated target bacterium, Pseudoalteromonas. First, baseline cytometric and microscopy data improved understanding of how infection conditions and host physiology impact populations in VT flow cytograms. Next, we extensively evaluated "and grow" capability to assess where VT signals reflect adsorption alone or wholly successful infections that lead to lysis. Third, we applied VT + Grow to a clonal virus stock, which, coupled to traditional plaque assays, revealed significant variability in burst size-findings that hint at a viral "individuality" parallel to the microbial phenotypic heterogeneity literature. Finally, we established a live protocol for public comment and improvement via protocols.io to maximally empower the research community. Together these efforts provide a robust foundation for VT researchers, and establish VT + Grow as a promising scalable technology to capture and characterize viruses from mixed community source samples that infect cultivable bacteria.

Complete genome sequence and immunoproteomic analyses of the bacterial fish pathogen Streptococcus parauberis.

Although Streptococcus parauberis is known as a bacterial pathogen associated with bovine udder mastitis, it has recently become one of the major causative agents of olive flounder (Paralichthys olivaceus) streptococcosis in northeast Asia, causing massive mortality resulting in severe economic losses. S. parauberis contains two serotypes, and it is likely that capsular polysaccharide antigens serve to differentiate the serotypes. In the present study, the complete genome sequence of S. parauberis (serotype I) was determined using the GS-FLX system to investigate its phylogeny, virulence factors, and antigenic proteins. S. parauberis possesses a single chromosome of 2,143,887 bp containing 1,868 predicted coding sequences (CDSs), with an average GC content of 35.6%. Whole-genome dot plot analysis and phylogenetic analysis of a 60-kDa chaperonin-encoding gene and the glyceraldehyde-3-phosphate dehydrogenase (GAPDH)-encoding gene showed that the strain was evolutionarily closely related to Streptococcus uberis. S. parauberis antigenic proteins were analyzed using an immunoproteomic technique. Twenty-one antigenic protein spots were identified in S. parauberis, by reaction with an antiserum obtained from S. parauberis-challenged olive flounder. This work provides the foundation needed to understand more clearly the relationship between pathogen and host and develops new approaches toward prophylactic and therapeutic strategies to deal with streptococcosis in fish. The work also provides a better understanding of the physiology and evolution of a significant representative of the Streptococcaceae.

Enhanced reliability of avian influenza virus (AIV) and Newcastle disease virus (NDV) identification using matrix-assisted laser desorption/ionization-mass spectrometry (MALDI-MS).

In-solution enzymatic and nonenzymatic digestion methods have been successfully implemented in matrix-assisted laser desorption/ionization-mass spectrometry (MALDI-MS)-based virus identification, extending to typing/subtyping of deadly influenza viruses. However, these methods are inefficient in obtaining more precise information on surface proteins of myxovirus particles, not only the hemagglutinin and neuraminidase of influenza virus but also the hemagglutinin-neuraminidase of Newcastle disease virus (NDV). Imbalances in viral protein composition cause ion suppression of tryptic fragments from low-abundant target proteins (surface proteins), adversely affecting reproducibility of mass spectra. Additionally, the coexistence of tryptic peptides from several proteins requires sophisticated statistical solutions for precise result interpretations. To circumvent these, we apply detergent-based (gel-free) partitioning of whole viruses into soluble surface proteins and insoluble virus materials, using differential centrifugation. MALDI-TOF or MALDI-TOF/TOF MS was applied to analyze tryptic peptides from separated viral proteins. In this study, we achieved type/subtype of avian influenza virus (AIV) within 5 h, based on 4 major proteins, by significantly reducing ion suppression and signal overlap from various protein sources. Hence, our approach can both yield dependable results and allow Web-based search engines to be directly employed, obviating the need for additional statistical strategy. Additionally, we demonstrate the utility of the method using NDV.

Viral Ecogenomics of Arctic Cryopeg Brine and Sea Ice.

Arctic regions, which are changing rapidly as they warm 2 to 3 times faster than the global average, still retain microbial habitats that serve as natural laboratories for understanding mechanisms of microbial adaptation to extreme conditions. Seawater-derived brines within both sea ice (sea-ice brine) and ancient layers of permafrost (cryopeg brine) support diverse microbes adapted to subzero temperatures and high salinities, yet little is known about viruses in these extreme environments, which, if analogous to other systems, could play important evolutionary and ecosystem roles. Here, we characterized viral communities and their functions in samples of cryopeg brine, sea-ice brine, and melted sea ice. Viral abundance was high in cryopeg brine (1.2 × 108 ml-1) and much lower in sea-ice brine (1.3 × 105 to 2.1 × 105 ml-1), which roughly paralleled the differences in cell concentrations in these samples. Five low-input, quantitative viral metagenomes were sequenced to yield 476 viral populations (i.e., species level; ≥10 kb), only 12% of which could be assigned taxonomy by traditional database approaches, indicating a high degree of novelty. Additional analyses revealed that these viruses: (i) formed communities that differed between sample type and vertically with sea-ice depth; (ii) infected hosts that dominated these extreme ecosystems, including Marinobacter, Glaciecola, and Colwellia; and (iii) encoded fatty acid desaturase (FAD) genes that likely helped their hosts overcome cold and salt stress during infection, as well as mediated horizontal gene transfer of FAD genes between microbes. Together, these findings contribute to understanding viral abundances and communities and how viruses impact their microbial hosts in subzero brines and sea ice.IMPORTANCE This study explores viral community structure and function in remote and extreme Arctic environments, including subzero brines within marine layers of permafrost and sea ice, using a modern viral ecogenomics toolkit for the first time. In addition to providing foundational data sets for these climate-threatened habitats, we found evidence that the viruses had habitat specificity, infected dominant microbial hosts, encoded host-derived metabolic genes, and mediated horizontal gene transfer among hosts. These results advance our understanding of the virosphere and how viruses influence extreme ecosystems. More broadly, the evidence that virally mediated gene transfers may be limited by host range in these extreme habitats contributes to a mechanistic understanding of genetic exchange among microbes under stressful conditions in other systems.

Antibiotic susceptibility and resistance of Streptococcus iniae and Streptococcus parauberis isolated from olive flounder (Paralichthys olivaceus).

The rates of antibiotic susceptibility and resistance were investigated in Streptococcus iniae and Streptococcus parauberis isolates obtained from diseased olive flounders (Paralichthys olivaceus) collected from fish farms in Jeju Island, Korea. Isolates of S. iniae (n=65) were susceptible to cefotaxime, erythromycin, ofloxacin, penicillin, tetracycline and vancomycin, as demonstrated by the minimum inhibitory concentration (MIC) test. Isolates of S. parauberis (n=86) were highly resistant to erythromycin (58% of the 86 isolates tested) and tetracycline (63% of the 86 isolates tested). Fifty-four isolates of tetracycline-resistant S. parauberis contained the tet(M/O/S) genes, of which 39 and 12 isolates contained the tet(M) and tet(S) genes, respectively, whereas 3 isolates contained both the tet(M) and tet(S) genes. Among the erythromycin-resistant isolates of S. parauberis (n=50) only 14 contained the erm(B) gene. These results suggest that the tet(S) and erm(B) genes of S. parauberis are involved in the acquisition of high-level resistance to erythromycin and tetracycline. Our findings reveal a high rate of antibiotic resistance among strains of S. parauberis and emphasize the need to develop an appropriate vaccine to reduce the use of antibiotics.

Soil Viruses Are Underexplored Players in Ecosystem Carbon Processing.

Rapidly thawing permafrost harbors ∼30 to 50% of global soil carbon, and the fate of this carbon remains unknown. Microorganisms will play a central role in its fate, and their viruses could modulate that impact via induced mortality and metabolic controls. Because of the challenges of recovering viruses from soils, little is known about soil viruses or their role(s) in microbial biogeochemical cycling. Here, we describe 53 viral populations (viral operational taxonomic units [vOTUs]) recovered from seven quantitatively derived (i.e., not multiple-displacement-amplified) viral-particle metagenomes (viromes) along a permafrost thaw gradient at the Stordalen Mire field site in northern Sweden. Only 15% of these vOTUs had genetic similarity to publicly available viruses in the RefSeq database, and ∼30% of the genes could be annotated, supporting the concept of soils as reservoirs of substantial undescribed viral genetic diversity. The vOTUs exhibited distinct ecology, with different distributions along the thaw gradient habitats, and a shift from soil-virus-like assemblages in the dry palsas to aquatic-virus-like assemblages in the inundated fen. Seventeen vOTUs were linked to microbial hosts (in silico), implicating viruses in infecting abundant microbial lineages from Acidobacteria, Verrucomicrobia, and Deltaproteobacteria, including those encoding key biogeochemical functions such as organic matter degradation. Thirty auxiliary metabolic genes (AMGs) were identified and suggested virus-mediated modulation of central carbon metabolism, soil organic matter degradation, polysaccharide binding, and regulation of sporulation. Together, these findings suggest that these soil viruses have distinct ecology, impact host-mediated biogeochemistry, and likely impact ecosystem function in the rapidly changing Arctic. IMPORTANCE This work is part of a 10-year project to examine thawing permafrost peatlands and is the first virome-particle-based approach to characterize viruses in these systems. This method yielded >2-fold-more viral populations (vOTUs) per gigabase of metagenome than vOTUs derived from bulk-soil metagenomes from the same site (J. B. Emerson, S. Roux, J. R. Brum, B. Bolduc, et al., Nat Microbiol 3:870-880, 2018, https://doi.org/10.1038/s41564-018-0190-y). We compared the ecology of the recovered vOTUs along a permafrost thaw gradient and found (i) habitat specificity, (ii) a shift in viral community identity from soil-like to aquatic-like viruses, (iii) infection of dominant microbial hosts, and (iv) carriage of host metabolic genes. These vOTUs can impact ecosystem carbon processing via top-down (inferred from lysing dominant microbial hosts) and bottom-up (inferred from carriage of auxiliary metabolic genes) controls. This work serves as a foundation which future studies can build upon to increase our understanding of the soil virosphere and how viruses affect soil ecosystem services.