Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mBio
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Host-Microbe Biology

Characteristics of Wetting-Induced Bacteriophage Blooms in Biological Soil Crust

Marc W. Van Goethem, Tami L. Swenson, Gareth Trubl, Simon Roux, Trent R. Northen
Jennifer B. H. Martiny, Editor
Marc W. Van Goethem
aEnvironmental Genomics and Systems Biology Division, Lawrence Berkeley National Laboratory, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tami L. Swenson
aEnvironmental Genomics and Systems Biology Division, Lawrence Berkeley National Laboratory, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gareth Trubl
bPhysical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Gareth Trubl
Simon Roux
cDOE Joint Genome Institute, Walnut Creek, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Simon Roux
Trent R. Northen
aEnvironmental Genomics and Systems Biology Division, Lawrence Berkeley National Laboratory, Berkeley, California, USA
cDOE Joint Genome Institute, Walnut Creek, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Trent R. Northen
Jennifer B. H. Martiny
University of California, Irvine
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mBio.02287-19
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Biological soil crusts (biocrusts) are photosynthetic “hot spots” in deserts and cover ∼12% of the Earth’s terrestrial surface, and yet they face an uncertain future given expected shifts in rainfall events. Laboratory wetting of biocrust communities is known to cause a bloom of Firmicutes which rapidly become dominant community members within 2 days after emerging from a sporulated state. We hypothesized that their bacteriophages (phages) would respond to such a dramatic increase in their host’s abundance. In our experiment, wetting caused Firmicutes to bloom and triggered a significant depletion of cyanobacterial diversity. We used genome-resolved metagenomics to link phage to their hosts and found that the bloom of the genus Bacillus correlated with a dramatic increase in the number of Caudovirales phages targeting these diverse spore-formers (r = 0.762). After 2 days, we observed dramatic reductions in the relative abundances of Bacillus, while the number of Bacillus phages continued to increase, suggestive of a predator-prey relationship. We found predicted auxiliary metabolic genes (AMGs) associated with sporulation in several Caudovirales genomes, suggesting that phages may influence and even benefit from sporulation dynamics in biocrusts. Prophage elements and CRISPR-Cas repeats in Firmicutes metagenome-assembled genomes (MAGs) provide evidence of recent infection events by phages, which were corroborated by mapping viral contigs to their host MAGs. Combined, these findings suggest that the blooming Firmicutes become primary targets for biocrust Caudovirales phages, consistent with the classical “kill-the-winner” hypothesis.

IMPORTANCE This work forms part of an overarching research theme studying the effects of a changing climate on biological soil crust (biocrust) in the Southwestern United States. To our knowledge, this study was the first to characterize bacteriophages in biocrust and offers a view into the ecology of phages in response to a laboratory wetting experiment. The phages identified here represent lineages of Caudovirales, and we found that the dynamics of their interactions with their Firmicutes hosts explain the collapse of a bacterial bloom that was induced by wetting. Moreover, we show that phages carried host-altering metabolic genes and found evidence of proviral infection and CRISPR-Cas repeats within host genomes. Our results suggest that phages exert controls on population density by lysing dominant bacterial hosts and that they further impact biocrust by acquiring host genes for sporulation. Future research should explore how dominant these phages are in other biocrust communities and quantify how much the control and lysis of blooming populations contributes to nutrient cycling in biocrusts.

INTRODUCTION

The Southwestern United States is predicted to experience alterations in the timing of rainfall events with concomitant increases in soil surface temperatures (1). In arid lands, moisture pulses are infrequent, and yet their timing, intensity, and duration dictate biological activity (2). A significant avenue of recent microbial ecological research has been elucidation of the responses of microbial communities to various manifestations of climate change, such as increased rainfall, and of how these responses in turn affect ecosystem processes (3–5).

Biological soil crusts (referred to as biocrusts) are dominant topsoil microbial communities in deserts and serve as the basis for soil surface colonization, stabilization, and activity (6). Biocrusts can be comprised of photosynthetic cyanobacteria, bryophytes (mosses and liverworts), lichens, fungi, a myriad of heterotrophic lineages that cooperate to cycle both soil carbon and nitrogen, and other elements (7). Our previous research in biocrusts dominated by cyanobacteria has shown that biocrust wetting results in significant shifts in microbial community structure (8) and dramatically alters microbial metabolic outputs within the first day of wetting (9). The most dramatic of these observations is that multiple Firmicutes lineages bloom and then collapse following wetting, which we hypothesize is due to bacteriophage (phage) predation.

Recent studies using cultivation-independent metagenomic sequencing have revealed the roles that viruses play across terrestrial biomes, including predation by giant viruses in forests (10) and viral predation of dominant and ecologically relevant microbial lineages along a permafrost thaw gradient (11, 12) and of dominant tailed phages in Antarctic soils (13, 14). Viral metagenomics offers insight into viral diversity within complex samples (15) and unlocks the potential for novel viral diversity to be uncovered without the need for cultivation (16).

We have previously shown that Firmicutes spp. respond positively to a biocrust wetting event (8), which is thought to be due to their emergence from dormancy (17). However, we noted previously that many Bacillus strains declined in abundance soon after the initial wetting event whereas others continued to increase in abundance until a sudden decline occurred (9). We hypothesized that the increase in phage activity upon wetting leads to the predation of blooming Firmicutes populations in biocrust. Several lines of evidence motivate this view: First, phages that infect Bacillus in hot desert soils in both the Mojave Desert and the Namib Desert have been described previously (18, 19). Second, phages have been shown to target blooming bacterial populations in aquatic ecosystems (20), including Prochlorococcus spp. in the ocean (21, 22) and blooms of the cyanobacterium Microcystis aeruginosa in freshwater lakes (23). Finally, the “kill-the-winner” hypothesis posits that viruses act as moderating agents within a community by attacking the most active population of prokaryotes, thus making the most of the limited resources available to them (24, 25). Accordingly, we expected biocrust phages to target the most active hosts following a latent period. To test this concept, we mined microbial and phage sequences from metagenomes obtained from biocrust exposed to a controlled laboratory wetting event and tracked their responses over 2 days. The samples originated from the Moab Desert, a site of long-term research aimed at exploring the effects of climate change in arid regions in Utah in the Southwestern United States.

RESULTS

Our 24 biocrust communities were subjected to a simulated wetting experiment, which we leveraged to explore microbial and phage responses along the 2-day time series. The assembled sequence data comprised 28.4 million contigs (∼20 Gbp) (see Table S1A and B in the supplemental material) from which we reconstructed 40 nearly complete bacterial metagenome-assembled genomes (MAGs). According to the Genomic Standards Consortium standards (26), all of the MAGs were of medium quality, as most lacked an identifiable 16S rRNA gene, and were >75% complete with <10% contamination (Table S2). These genomes encompassed bacterial members of Firmicutes (n = 16), Bacteroidetes (n = 7), Proteobacteria (n = 5), Acidobacteria (n = 4), Chloroflexi (n = 2), and Actinobacteria (n = 2) and included a single representative of Deinococcus-Thermus (see Fig. S1A and B in the supplemental material). The assembled genomes averaged 4.4 Mbp in size, with a range of relative abundances from high abundance (1.33% of assembled metagenome contigs) to rare (0.0067%) community members. We interrogated our MAGs for genes of interest relating to nutrient cycling, sporulation, osmoregulation, carbon starvation, and heat shock tolerances (see Results in Text S1 in the supplemental material).

TEXT S1

Supplementary Results and Discussion pertaining to our co-occurrence network analysis and descriptions of functional genes from the biocrust metagenome-assembled genomes (MAGs). Download Text S1, DOCX file, 0.02 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

FIG S1

(A) Neighbor-Joining tree indicating the placement of gyrB gene sequences retrieved from the biocrust MAGs among 1,000 reference sequences. MAGs are colored according to the biocrust successional stage from which they were reconstituted. Bootstrap values determined after 1,000 iterations of the tree build are presented. (B) Important functional genes identified across all biocrust MAGs. Genes are colored according to abundance level, with red indicating the most common genes and white the less common genes. MAGs are sorted at the phylum level. Download FIG S1, TIF file, 1.6 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

TABLE S1

(A) Raw sequencing statistics for each of the 28 biocrust metagenomes. (B) Assembly quality and 16S rRNA gene descriptions for each of the 28 biocrust metagenomes. Download Table S1, XLSX file, 0.02 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

TABLE S2

Descriptions of each biocrust MAG that was reconstituted from the coassembled metagenomes. Download Table S2, XLSX file, 0.02 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

The most phylogenetically diverse biocrust communities were the prewetting (“dry”) metagenomes and those collected 3 min and 9 h after wetting (∼58% of the community was sequenced on average as estimated using Nonpareil [27]). Biocrust metagenomes from later time points, 18 and 42 h after wetting, captured ∼64% of the community, while the final time point, 49.5 h after wetting, produced communities that were less complex, whereby 71% of the community was sequenced (Fig. S2A). Diversity estimates corroborated these rarefaction estimates and showed significant decreases in α-diversity following wetting (analysis of variance [ANOVA], P < 0.05) (Table 1). The 49.5-h group had significantly lower α-diversity than any of the other groups (Tukey post hoc test, P < 0.05). In contrast, the number of single-copy ribosomal genes (large and small subunits) increased significantly with wetting (P < 0.05). We also observed a significant increase in β-diversity over the course of our wetting experiment (within treatment variation). Note that the decrease in diversity may also have been due to our limited sampling and sequencing of rare community members as their relative abundances decreased dramatically after wetting.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Statistics describing the taxonomic diversity within the biocrust metagenomes grouped according to time after wetting

FIG S2

(A) Nonpareil estimates of sequencing coverage. Samples are colored according to their collection time point after wetting. (B) Nitrogen cycling genes (ORFs) found across our 28 biological soil crust metagenomes indicated as total counts of significant hits. Download FIG S2, TIF file, 2.8 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

We observed significant differences in community composition according to time after wetting (permutational multivariate analysis of variance [PERMANOVA], P < 0.05) (Fig. 1, left panel). The effect of wetting induces changes in composition that are strong enough to supersede those of the maturity stage of the biocrust, which developed over decades (P > 0.05). The increase in β-diversity over time was also observable as individual communities separated within ordination space (nonmetric multidimensional scaling [NMDS] ordinations based on Bray-Curtis dissimilarities of 16S rRNA gene content) (Fig. 1, left panel), which was the result of taxonomic divergence over time and of overall increases in 16S rRNA gene content (Fig. 1, right panel).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

(Left) The changes in taxonomy over the course of the wetting experiment are observable in NMDS ordination space. Biocrust communities grouped according to time after wetting, which was sufficient to explain the differences in community composition (PERMANOVA; P < 0.05). (Right) The communities grew more dissimilar as the copy numbers of 16S rRNA genes increased over time.

Overall, the biocrust communities encompassed 23 bacterial phyla and 1 archaeal phylum (Thaumarchaeota). The metagenomes also provided contigs belonging to Ascomycota (filamentous, saprobic Fusarium and Purpureocillium lineages), metazoa (Tardigrada), and mosses (Physcomitrella patens). The low diversity of fungi might have been due to database bias, and using a dedicated fungal database such as UNITE might have revealed more fungal diversity in these samples. Next, we explored the taxonomic shifts that occurred during biocrust wetting by tracking changes in 16S rRNA gene content in the unassembled reads. On the basis of our previous 16S rRNA gene survey, we anticipated a significant enrichment of Bacillales and other Firmicutes taxa ∼18 h after wetting (8, 9). Notably, these blooms should have been accompanied by substantial reductions in the proportion of filamentous cyanobacterial community members across all samples as nonbloomers are outcompeted by faster growers. However, Cyanobacteria are known to be sensitive to wetting (28, 29), and excessive cellular swelling and lysis are thought to cause their negative response. The dry and 3-min biocrust metagenomes were dominated by photosynthetic Cyanobacteria and by heterotrophic Actinobacteria, Bacteroidetes, and Proteobacteria (Fig. 2, left panel). Overall, the biocrust communities encompassed 23 bacterial phyla and 1 archaeal phylum (Thaumarchaeota). The metagenomes also provided contigs belonging to Ascomycota (filamentous, saprobic Fusarium and Purpureocillium lineages), metazoa (Tardigrada), and mosses (Physcomitrella patens).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Exposure to a sudden wetting event had strong effects on the taxonomy of the biocrust communities across all successional stages. (Left) Relative abundances of microbial phyla in the biocrust communities. Archaea and eukaryotes are classified at the domain level, while Proteobacteria are distinguished at the class level. (Right) The responses of Firmicutes genera are shown as the proportion of 16S rRNA gene copies per metagenome.

The wetting event drove decreases in the relative abundances of all nonblooming phyla. Notably, the relative abundance of Cyanobacteria decreased significantly from ∼64% of the community in the dry samples to just 0.23% 49.5 h after wetting (P < 2.30e−10). We observed striking, significant decreases in the relative abundances of most cyanobacterial lineages over the course of the experiment, culminating in a nearly complete loss of the genera Microcoleus (P < 2.15e−7), unclassified Cyanobacteria (P < 4.78e−5), Tychonema (P < 2.36e−3), and Subsection III, Family I members (P < 0.025) by 49.5 h (Fig. 2, left panel). Cyanobacteria are known to be sensitive to wetting (28, 29) and excessive cellular swelling, and lysis may have caused their negative response here. The reduction in the abundances of cyanobacterial members also had dramatic effects on network topology across the experiment on the basis of co-occurrence analyses of taxonomic information (see Text S1, Fig. S3A to C, and Table S5).

FIG S3

Co-occurrence networks of taxonomic data according to time after wetting. (A) The initial network comprising all dry communities and communities present 3 minutes after wetting. (B) The “Intermediate” network comprising all communities that were present 9 h and 18 h after wetting. (C) The final network comprising all communities that were present 42 h and 49.5 h after wetting. Nodes are colored according to their clusters in the initial network, while rings around nodes indicate their taxonomic group. Nodes are sized according to their relative abundances in each community. Map edges are colored such that positive co-occurrences are represented by blue and solid lines, while negative co-occurrences are represented by both red coloring and dashed lines. Clusters were identified by MCODE and are grouped in circular topologies. Download FIG S3, TIF file, 2.3 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

TABLE S5

Network statistics about each co-occurrence network generated by SparCC. Download Table S5, XLSX file, 0.01 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

Only members of the Firmicutes significantly increased in relative abundance after wetting; their levels increased from 1.07% of 16S rRNA genes in dry samples to 82.3% of 16S rRNA genes at 49.5 h (relative abundance of 16S rRNA genes, ANOVA, Bonferroni corrected, P < 1.77e−6) (Table S3A). The presence of the Firmicutes bloom was pronounced across all biocrust successional stages, but the responses of individual genera were more nuanced. Notably, Brevibacillus showed a steady increase in abundance after wetting (Fig. 2, right panel), whereas Paenibacillus showed slower growth. The population of Bacillus became dominant 9 h after wetting but had collapsed by 49.5 h. We speculated that these diverse responses to wetting would be reflected by increases in the proportions of functional genes affiliated with the changes in the proportions of these taxa.

TABLE S3

(A) Microbial taxa that showed significant changes in relative abundance as a result of wetting. (B) Functional genes that showed significant changes in relative abundance as a result of wetting. Download Table S3, XLSX file, 0.03 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

We tracked how changes in taxonomic composition affected the relative contributions of functional genes. Most of the upregulated genes in our analysis were involved in sporulation (Fig. 3), suggesting that previously dormant Firmicutes had replicated within 9 h after wetting. Specifically, we observed significant increases in the abundance of sporulation cluster stage IIIA genes (AB, AC, AD, and D), stage V genes (spoVS), spore maturation protein A, and genes encoding the spore cortex-lytic enzyme CwlJ (Fig. 3). Consistently, the dramatic reduction in the populations of cyanobacterial taxa over the course of the wetting experiment led to significant reductions in the relative abundances of genes involved in the assembly of photosystem I (ycf4 and psaLK2) and photosystem II (psbADMPVY) and of genes involved in the production of bilin (pcyA) and the phycobilisome (nblB and cpeS). The relative abundances of genes associated with the cyanobacterial circadian clock (input kinase A) and myxoxanthophyll biosynthesis (cruF) and genes encoding products associated with photorespiration (d-glycerate 3-kinase and glycolate oxidase) were also significantly reduced (Table S3B). Strikingly, numerous genes involved in carbon dioxide (CO2) uptake (carboxysome) were also significantly reduced in abundance. These included RuBisCO chaperonin genes (rbcX) as well as CO2-concentrating genes (ccmO), high-affinity carbon uptake genes (hatR), and low-affinity CO2 hydration genes (cphX). Overall, 95 genes decreased in relative abundance over the wetting experiment, while only 22 genes became significantly more abundant following wetting (Table S3B).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Heat map indicating the distribution of major functional genes in biocrust metagenomes following soil wetting. Count data are log normalized to enable comparisons across different sequencing efforts. Genes that are overrepresented are represented by darker green colors; white data represent genes that are underrepresented (i.e., the genes with the lowest relative abundance levels). The x axis indicates biocrust communities following wetting, i.e., ranging from the dry state to 49.5 h after wetting. The y axis data include all functional genes that differed significantly in abundance across the wetting event.

Next, we used VirSorter to mine viral signals from our 28 metagenomes (including the 4 bundle metagenomes), and the mined signals were coassembled into 5 data sets to assist viral sequence recovery. This procedure produced a total of 1,849 unambiguous viral sequences (see Materials and Methods). The removal of short (<10-kb) contigs, together with manual curation (see Materials and Methods), refined our data set to 1,210 viral genomes (ViGs) or genome fragments.

To place our ViGs within the known, relevant diversity of publicly available viral genomes, we created a virus genome-based network based on shared protein content (30, 31) by combining our contigs with the viral genomes and genomic fragments present in the RefSeq viral database (v85) (32). The analysis produced 40,796 virus clusters (VCs) representing approximately genus-level taxonomy, of which 2,294 (6%) contained biocrust viral contigs, representing 504 uncultivated viral genomes (UViGs). Overall, 148 (29%) of the 504 biocrust UViGs that we identified were found to be taxonomically affiliated using these VCs. Clustering the 504 UViGs at 95% average nucleotide identity (ANI) and 85% alignment fraction (AF) (20) produced 433 viral operational taxonomic units (vOTUs), of which 49 could be assigned at the genus level. All classified vOTUs belonged to Caudovirales families (double-stranded DNA [dsDNA] tailed phages), although six vOTUs were broadly categorized as representing dsDNA viruses (Fig. 4, left and top right panels; see also Table S4A).

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

(Left) Network clustering of recovered virus genomes and genome fragments with RefSeq viral genomes (v85) based on shared predicted protein content. Major viral families (nodes) are indicated by unique colors and shapes, with vOTUs found in this study indicated by green circles. Edges (lines between nodes) indicate significant shared protein content. (Top right) Stacked bar chart indicating the number of UViGs that were recovered across the wetting experiment. (Bottom right) Correlation between the relative proportions of Firmicutes and UViGs per time point.

TABLE S4

(A) List of all UViGs (uncultivated viral genomes) derived from this study, including the required statistics prescribed by the consortium behind the minimum information about an uncultivated virus genome (MIUViG) definitions. (B) Sequences of the sporulation genes identified both UViGs and bacterial MAGs. Download Table S4, XLSX file, 0.04 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

The longest circular UViGs belonged to three identical Spbetavirus phages (163,801 bp), which are Siphoviridae phages known to infect soil Bacillus (33). These UViGs were reconstructed independently biocrust metagenomes seen in the early, late-middle (late-mid), and late successional stages, while the early-middle (early-mid) metagenome also contained a smaller Spbetavirus contig of 137,161 bp. the Siphoviridae phages were the most abundant, with 89 vOTUs spanning 17 genera. Podoviridae were represented by 19 vOTUs, including three T7 viruses, while 21 vOTUs assigned to the Myoviridae family were found to include three genera (Abouovirus, Bxz1virus, and P2virus). A further 10 genomes were coarsely assigned to the Caudovirales. Mapping UViGs back to individual metagenomes revealed that all of the phages were detected only after wetting; i.e., no phages were recovered from the bundles or the dry samples (Fig. 4, top right panel), which culminated in the greatest number of UViGs (n = 50) being recovered more than 2 days after wetting. For example, the full-length spbetaviruses emerged at 42 h and 49.5 h after wetting, which correlates with the decline in the abundances of the putative Bacillaceae hosts at those time points (Fig. 2, right panel). Andromedaviruses (length range, 47,888 to 48,100 bp) were found in all successional stages at 18, 42, and 49.5 h after wetting and have previously been shown to infect Bacillus hosts (33). Combined, these UViG data indicated that Firmicutes were predominant targets, particularly at the end of the wetting experiment, which led us to further explore this putative link. The classification of our viral genomes provided primary evidence of phages targeting Firmicutes such as Bacillus, Staphylococcus, and Brevibacillus.

Putative auxiliary metabolic genes (AMGs; genes that originate from microbes and can increase fitness) that code for sporulation proteins (spoIIIE) were identified in 9 UViGs. To confirm that the AMGs were of viral origin, we checked that flanking regions present both upstream and downstream contained common virus genes (integrases, structural genes, or terminases), had significant matches to the gene of interest (E value < 1 × 10−3; bit score > 60), and represented the only AMG per contig. BLASTP queries indicated that all 9 sporulation genes (spoIIIE) were assigned to Bacillus strains. We used the sporulation AMGs to explore the proportion of phages infecting the blooming Firmicutes. As an initial proof of concept, we found that the numbers of UViGs correlated strongly with the relative abundances of Firmicutes across the wetting experiment (r = 0.762). That is, as the proportion of Firmicutes increased throughout the wet-up, so too did the number of recoverable phage genomes (Fig. 4, bottom right panel; see also Fig. S4B).

FIG S4

(A) Flow diagram illustrating our methodological pipeline for analyzing bacterial and viral components of the biocrust communities. (B) Plot of Fig. 4C with a logarithmic fit after omitting the dry sample point containing a zero value (R2 = 0.8935). Download FIG S4, TIF file, 1.6 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

Genes for sporulation were detected in 15 of the 16 Firmicutes MAGs. Genes involved in spore DNA protection (dipicolinate synthase subunits A and B, dpaAB) and small acid-soluble spore proteins (SASP1, SASPS6, and gpr) were detected exclusively in the Firmicutes genomes and are well-known components of endospore production (34, 35). DpaAB genes code for dipicolinic acid, which allows endospores to maintain dormancy under harsh conditions. Sporulation-related genes, including spoIIIE and gerE, were among the most highly overrepresented genes in the metagenomes examined in our study by the end of our wetting experiment (Tables S3B and S4B), reiterating the increase in the relative abundances of spore-formers following wetting. A maximum likelihood tree indicated that the spoIIIE AMGs found in our UViGs showed high similarity to those found in the Firmicutes genomes (Fig. 5). We confirmed the clustering result by predicting protein structures of each spoIIIE gene using I-TASSER and comparing models using TM-align, which showed a similarity score of 0.488.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Maximum likelihood tree indicating the phylogenetic placement of spoIIIE genes retrieved from UViGs with their nearest bacterial neighbors (from Firmicutes MAGs). Protein structures were predicted using I-TASSER, and the first predicted models are presented. Tree lineages are colored blue for bacterial sequences and orange for viral sequences.

Additionally, we searched for CRISPR-Cas repeats and proviral genes within the reconstructed MAGs to link viruses to potential hosts. We found evidence of CRISPR-Cas repeats in 16 MAGs, which links them to previous or ongoing infection events by biocrust viruses (11, 36). Six of the 16 Firmicutes MAGs contained CRISPR elements (Table S2). Aligning the CRISPR-Cas repeats and spacers from the MAGs with our virus contigs revealed perfect homology between a Phietavirus (EM_51679) and the Anoxybacillus MAG EM_40. Both the repeat sequence (5′-CTTCACATCACACATAGTTTCCTCAGAAAC-3′) and three spacers (data not shown) were identical between the virus and its cognate host, which offers a direct link to the microbial host in silico. Next, we searched for prophage elements within each MAG to uncover temperate phages that utilized a lysogenic infection cycle. We found 12 MAGs with viral elements within their genomes, of which 7 were Firmicutes (Table S2).

DISCUSSION

The results presented here provide an initial assessment of the phage diversity in biocrust communities—dominant photosynthetic topsoil communities found around the globe. This description of dominant phage lineages in biocrust further extends our knowledge of viruses in important soil ecosystems (37), which lags considerably behind what is known about marine viruses (20, 38). Our work sought to resolve the contradicting fates of taxa within the Firmicutes after biocrust communities were exposed to a simulated wetting event (8), which we predicted might have been due to phage predation. Here, we provide evidence of a predator-prey relationship between Bacillus and their phages. Comparing the 16S rRNA gene copy numbers of Bacillus throughout wet-up to the numbers of their known UViGs showed that phage abundance peaked as that of their hosts began to decline (Fig. 6).

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Predator-prey relationship between Bacillus organisms and their phages. Bacillus rRNA copy number is provided on the left-hand y axis with error bars calculated across four metagenomes from each time point. The numbers of UViGs linked to Bacillus infection are indicated on the right-hand y axis.

The presence of proviral sequences within Firmicutes genomes suggests that lysogenic replication strategies are common among biocrust viruses. Temperate viruses shift their replication strategy from lysogenic to lytic replication as bacterial host production increases (39). We found prophages in 12 of the 40 MAGs. Seven Paenibacillus MAGs contained proviral sequences, and each of these contained multiple prophages. While the genus Bacillus abundance declined, both Paenibacillus and Brevibacillus continued to increase in relative abundance throughout the wet-up. According to the kill-the-winner model, we would expect phages to attack the fastest-growing host population first (i.e., Bacillus). We anticipate that phages would target members of the Paenibacillus and Brevibacillus after the Bacillus populations have been depleted on the basis of the trajectory of these genera over 2 days (Fig. 2, right panel).

For Bacillus hosts, phage-encoded sporulation genes could play an additional role in these dynamics by providing a mechanism for lysis with concomitant phage and spore release where uninfected host spores would represent a reservoir of future host cells. Note that we cannot differentiate between population reduction mediated by lysis from sporulation and that mediated by some combination of factors with these data. In addition, spores that include integrated phage genomes would also provide additional protection, which may be important in the harsh biocrust environment. Regardless, the lysis of Firmicutes cells probably stimulates the biocrust community as metabolites, such as nitrate, are released into the soil. Conceivably, this viral shunt could provide essential nutrients to mosses or other higher eukaryotes, as predicted in marine ecosystems (40), which has the potential to influence ecological services of the crust.

Our metagenome-assembled genomes also revealed that Firmicutes-infecting phages acquired genes horizontally from their hosts’ genomes. Specifically, nine spoIIIE genes in phage genomes link them to Firmicutes infections. The spoIIIE gene codes for a DNA translocase that plays a critical role during sporulation by assisting the translocation of chromosomal DNA from mother cells into the forespore during the final stages of engulfment (41). It has been shown that marine cyanophages contain photosynthetic genes which maintain essential host functions (22). These spoIIIE genes may have a similar role in these biocrusts if sporulation is an important component of the life history of these Bacillus phages, as endospores can be carriers of phage genomes, thus offering protection from various environmental stresses (42, 43). Alternatively, it is possible that these phage-encoded proteins were originally sporulation related but now perform another function related to DNA export and/or ATP-dependent motor activity. AMGs encoding sporulation machinery (spoVS) in vOTUs from permafrost have been detected previously (12), and it would be intriguing to investigate how widespread these genes are in other soil systems as they may influence host sporulation dynamics (42).

As Firmicutes bloom during precipitation events, they are probably preyed upon by phages, exerting control on their abundance throughout the year. Biocrust saturation after rainfall releases adsorbed phages, which may initiate repeated cycles of lytic replication (37). The drying soil likely draws phages into small pores by shrinking water films (44). Soil microheterogeneity may further limit lytic virus-host interactions as cells become heterogeneously distributed within the dry soil matrix. This paradigm is especially true of desert biocrusts, which experience prolonged periods of dormancy for much of the year and reawaken after moisture events (45, 46). Future research should investigate these dynamics in the field, particularly as our long-term storage of samples (∼11 months in the dark) may have influenced the stress levels that the microbes experienced. Moreover, it is unclear whether blooms of Firmicutes occur as rapidly in naturally occurring biocrust communities or if this response is mitigated by other factors.

Previous surveys of viral populations from desert biotopes, such as cyanobacterium-dominated hypoliths (18) and exposed desert soils (47), revealed a high proportion of Caudovirales (dsDNA tailed phages), which appear to be dominant dsDNA phages in biocrust as well. We found that the family-level hierarchical abundance of Siphoviridae > Myoviridae > Podoviridae was remarkably similar to that described in Antarctic hypoliths (14), potentially indicating that these tailed phages are core constituents of “hot spots” of microbial activity in desert ecosystems. Moreover, the average genome sizes of phages identified here (average length, ∼43 kb) are like those described in Namib soils identified from pulsed-field gel electrophoresis profiles (55 kb). Additionally, those studies provided striking evidence that cyanophages are a potentially low-abundance component of desert viral communities despite the abundance of putative cyanobacterial hosts in hypoliths (47) and biocrust, including the genera Microcoleus, Scytonema, Schizothrix and Nostoc (7, 48). These results should be interpreted with caution, however, as our long-term storage of biocrust samples in the dark for 11 months may have diminished the responsiveness of certain taxa to wetting. Moreover, most knowledge regarding cyanophage diversity stems from oceanic surveys, which may have introduced biases into the RefSeq database that was used here. The number of studies on soil viruses continues to increase rapidly, which certainly offers researchers the potential to detect novel cyanophages in these underexplored soil ecosystems.

Summary.Biocrust wetting induced a bloom of Firmicutes and their phages. Our data suggest that Firmicutes emerge from spores, whereby fewer individuals remain in sporulated states after 2 days, which is reflected in significant enrichment for genes relating to sporulation. In contrast, the reduction in the relative abundance of Cyanobacteria resulted in significant depletion of functions relating to both carbon and nitrogen fixation, and this was observed in collapses in co-occurrence network structures (see Fig. S3C in the supplemental material). Our investigation of phage diversity in biocrust, performed here for the first time to our knowledge, did not link phage predation to cyanobacterial decline but instead revealed that biocrust phages primarily infect blooming Firmicutes, particularly members of the Bacillus genus. The phage-bacterium dynamics observed here appear to follow the kill-the-winner model, in which phages expand on the fastest-growing host population in biocrust. The collapses of population blooms are often mediated by phage infection and lysis and can occur on very short time scales such as in marine systems where phages target blooming Synechococcus species within 30 min (40). Here, we showed that the increase in Firmicutes abundance over 2 days is closely followed by peaks in the frequency of recoverable Caudovirales phage genomes targeting different members of the phylum, followed by a decline of specific Firmicutes lineages as their phages replicate. The presence of putative auxiliary metabolic genes associated with sporulation in these phages suggests that the sporulation process may be an important component of the life history of these phages. Future research could further investigate the contributions of sporulation to phage dispersal and maintenance within endospores. Additional analyses could consolidate the fate of the biocrust communities by investigating how biocrusts return to their initial steady state following drying.

MATERIALS AND METHODS

Sample collection, experimental setup, and DNA extraction.Details of the sample collection procedure, the biocrust wetting experiment, and DNA extraction protocols have been described previously (9) and are illustrated in Fig. S4A in the supplemental material. Briefly, 28 petri dishes (6 cm2 by 1 cm in depth) were used to core biocrust samples from the Green Butte Site near Canyonlands National Park (38°42′54.1′′N, 109°41′27.0′W, Moab, UT, USA) in 2014. The site is an area of long-term ecological research interest aimed at exploring climatic changes in arid regions. Biocrust samples (28 petri dishes) were collected along a maturity gradient of Cyanobacteria-dominated biocrusts ranging from lightly colored, young crust to dark mature crust, which encompassed four maturity stages overall. Dry samples collected in the field were maintained in a dark desiccation chamber in the laboratory for 11 months until required. Biocrust samples underwent a laboratory (12 h light/12 h dark cycle) wetting event by wetting of each soil sample (0.5 g) with liquid chromatography-mass spectrometry (LC-MS)-grade water (1 ml). Soil was removed at five time points following wetting (in addition to dry and bundled samples) for each successional stage (3 min, 9 h, 18 h, 42 h, and 49.5 h after wetting). After pelleting of soil by centrifugation (5,000 × g for 5 min) and removal of porewater, DNA was extracted from soils (0.25 g) using a DNeasy PowerSoil kit (Qiagen, Hilden, Germany), resulting in 100 μl of eluted DNA. Library preparation and sequencing were performed at the QB3 facility at the University of California, Berkeley (UC Berkeley), using an Illumina HiSeq 4000 system.

Metagenomic sequencing and analysis.Metagenomic sequencing provided ∼413 Gbp of Illumina HiSeq data in 2,754,282,840 raw sequences across 28 biocrust samples. The data included four separate metagenomes for microbial communities collected before wetting and those collected 3 min, 9 h, 18 h, 42 h, and 49.5 h after wetting. Sequencing was performed on an Illumina HiSeq 4000 instrument using 2 × 150-bp chemistry. Raw sequences were filtered using Prinseq v0.20.4 (49), which eliminated reads with phred scores of <20 for six consecutive bases (–min_qual_mean), singletons, and those with ambiguous bases (ns_max_n). Estimates of sequence coverage were obtained using Nonpareil v3.30 (27). High-quality sequences were assembled under default conditions using MEGAHIT v1.1.3 (50) as previously suggested for highly complex soil samples (51). These 24 metagenomes were also coassembled according to their maturity stage, i.e., dry, early, early-mid, late-mid, and late successional stage biocrust, as physically identified a priori during the field sampling expedition. Statistics about each assembled metagenome were obtained using MetaQUAST v4.6.3 (52) and are available in Table S1B in the supplemental material. Open reading frames (ORFs) were extracted from each assembled metagenomes with Prodigal v2.6.3 (53). Annotation of contigs was performed using DIAMOND v0.9.14.155 (54) with queries made against the entire NCBI nonredundant nucleotide database. Here, annotations were provided at E values of <0.00001 (-e), while the five best hits were retained (-k). MEGAN software v6.12.5 was used to visualize functional and taxonomic attributes of communities (55). We extracted all 16S rRNA genes from the curated metagenomes using SortMeRNA (56). Finally, we resolved exact sequence variants (ESVs) using DADA2 (Divisive Amplicon Denoising Algorithm [57]), with default parameters except for truncLen (120, 150) and maxEE (2, 5). We assigned taxonomy against the entire SILVA reference database (58) after random rarefaction of each community to the lowest sequencing effort (8,805).

Network construction.We constructed co-occurrence networks at the genus level using our log-transformed taxonomic data. Networks were constructed using Sparse Correlations for Compositional (SparCC) data (59). We chose to divide our 24 metagenomes according to time after wetting. Individual networks were created for the beginning (dry and 3 min after wetting), middle (9 h and 18 h after wetting), and end (42 h and 49.5 h after wetting) of the experiment. We used significant SparCC values (P < 0.05) that were greater than 0.6 or less than −0.6 and included only genera present in at least half of the samples within our OTU table. All networks were visualized using Cytoscape v3.6.1 (60), and highly connected modules were automatically identified using MCODE (61).

Reconstructing genomes from metagenomes.Due to our ultradeep sequencing effort, we were able to reconstruct metagenome-assembled genomes (MAGs) from our assembled biocrust communities using MetaBAT v2.12.1 (62). We increased the sensitivity of the binning process by increasing the minimum contig length to 5,000 bp (–m), the percentage of good contigs to 97 (–maxP), the minimum edge binning score to 85 (–minS), and the maximum number of edges per node to 250 (–maxEdges). We also included the –noAdd option, which prevented small contigs from being binned and prevented the creation of bins that were <1 Mbp in size. MAGs were checked for completeness and contamination using CheckM v1.0.8 (63). Medium- to high-quality MAGs were retained according to the conventional values for minimum information about a MAG (>75% completeness, <10% contamination) (26). MAGs were then functionally annotated online using Prokka v1.12 (64) and the RAST v2.0 server (65), with ORFs validated with UGENE (66). Quality statistics about the MAGs are presented in Table S2.

Phylogenetic tree building.Phylogenetic trees were built using universal bacterial protein gyrase subunit B (gyrB) extracted from MAGs (n = 19) and 2,878 reference sequences. The gyrB sequences were collated and aligned using MAFFT v7 (67). Neighbor-joining trees were built with bootstrap values reported, and visualization was performed in iTOL v3 (68). The phylogenetic placement of the MAGs was also evaluated using the Microbial Genome Atlas (MiGA) online (69).

Statistical analyses.We used STAMP v2.1.3 (70) and R v3.5.0 to make statistical comparisons between the communities according to both time after wetting and successional stage. Within the R statistical environment, we used the vegan package (71) to explore ecologically meaningful differences between communities according to time after wetting and successional stage. Bray-Curtis dissimilarity matrices were used to explore biocrust community structure after Hellinger transformations of species relative abundance data. PERMANOVA was used to test for differences in community structure on the basis of time after wetting (dry biocrust and biocrust at 3 min, 9 h, 18 h, 42 h, and 49.5 h after wetting) and successional maturity (early, early-mid, late-mid, and late) using the adonis function. We used STAMP to test for significant differences in taxonomic and functional composition between time points and successional stages and used ANOVA with Tukey-Kramer post hoc tests and Bonferroni corrections of the P value statistic.

To obtain viral contigs from biocrust metagenomes, we coassembled samples according to their characteristics and community complexity. The coassemblies were produced by combining the 4 filamentous bundles, the 4 dry biocrust samples, the 5 early successional stage biocrust samples, the 5 early-mid succession biocrust samples, the 5 late-mid succession biocrust samples, and the 5 late succession biocrust metagenomes into single, concatenated metagenomes. This produced 6 coassemblies to complement the 28 individual metagenomes.

Recovering and annotating viral contigs.Viral contigs were recovered from the 6 coassembled biocrust metagenomes using VirSorter 1.0.2, which relies on the detection of “hallmark” viral genes, such as those genes involved in the formation of the viral capsid (72). The Virome db database, without additional viral sequences as references, was used to recover viral contigs from all coassembled biocrust metagenomes. All contigs classified into category 1 or category 2 (complete viral contigs) were considered to be of unambiguous viral origin, while contigs in categories 4 and 5 (proviruses within microbial contigs) were modified to include only the predicted proviral regions, thus excluding host genomic signatures (see host analysis below). Only the viral contigs that scored in the highest confidence categories were retained and analyzed. The viral contigs within our metagenomes are likely sequences derived from double-stranded DNA phages and include prophage elements (temperate phages integrated within microbial chromosomes) and plasmids (extrachromosomal DNA) as well as lytic phages. All predicted viral contigs were retained if they were ≥10 kb in length. The most common class was composed of category 2 viral signatures which were “likely” predictions (n = 1,665), followed by category 1 (“most confident” predictions; n = 184). These viral signals were retrieved from our 28 metagenomes (including the 4 bundle metagenomes), which were coassembled into 5 data sets to assist viral sequence recovery. The removal of short (<10-kb) contigs refined our data set to 1,210 viral genomes or genome fragments. The identified viral protein sequences were clustered and taxonomically classified using vConTACT2 (30), which is supported by the International Committee on Taxonomy of Viruses (ICTV). We built a network on the basis of significant gene-sharing similarities (represented as edges) between viral genomes with contigs represented as nodes. We performed all-vs-all BLASTp queries for comparisons between our 2,140 VirSorted sequences and all predicted viral proteins in the NCBI RefSeq database (v85) with an E value of 1 × 10−3. Protein clusters were defined using the Markov Clustering Algorithm (MCL) and then ClusterONE as described previously (73). These viral genomes represent both high-quality draft genomes (single genome fragment greater than 90% of the expected genome size) and viral genome fragments (single genome fragments less than 90% of the expected genome size or of unknown genome size), defined within the Genomics Standards Consortium framework as the minimum information about an uncultivated virus genome (MIUViG) (16). The UViGs identified by this procedure were visualized within the Cytoscape v3.6.1 program. We also mapped our MAGs and UViGs back to our metagenomes to calculate mapping rates and to determine which genomes were recovered from each time point using Bowtie 2 (74). UViGs were clustered into viral operational taxonomic units (vOTUs) at 95% ANI and 85% AF.

Viral host prediction and CRISPR-Cas analyses.To link viruses with their cognate hosts, we explored the distribution of CRISPR-Cas elements within each annotated MAG and metagenome. Clustered regularly interspaced short palindromic repeats (CRISPR) represent phage infection legacy and thus provide a palpable link between phages and their cognate hosts (36). CRISPRs identified with MinCED 0.3.2 were confirmed using BLASTn queries by allowing for 1 mismatch. We also searched for prophage elements within each MAG using VirSorter, within the CyVerse Discovery Environment.

Auxiliary metabolic gene (AMG) analysis.We considered the 148 viral genomes and genomic fragments for exploration of auxiliary metabolic genes (AMGs) on the basis of VirSorter-predicted contigs (only VirSorter categories 1 and 2, thus excluding proviral contigs). This reduced the possibility of including host genes in our analysis. The identified sporulation related AMGs were reannotated to confirm their classifications through BLASTP comparisons whereby those with E values of 1 × 10−3 and bit scores of >60 were retained. Next, the genomic regions both upstream and downstream of the retained genes were explored for common viral genes, including those encoding terminases, integrases, and viral structural genes. Finally, maximum likelihood trees were built using the AMG sequences and those extracted from the host genomes, as described in the “Phylogenetic tree building” section. The protein structures of our AMGs of interest were predicted with I-TASSER (75), after which we tested for correct predictions using TM-align (76).

Data availability.The sequencing data used in this study are available at the NCBI Sequence Read Archive (project accession number PRJNA395099; sample accession numbers SRX3024638 to SRX3024657; https://www.ncbi.nlm.nih.gov/bioproject/395099).

ACKNOWLEDGMENTS

This work was supported by funds provided by the Office of Science Early Career Research Program, Office of Biological and Environmental Research, of the U.S. Department of Energy under contract no. DE-AC02-05CH11231 to Lawrence Berkeley National Laboratory. Work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under contract no. DE-AC52-07NA27344 (to G.T.).

M.W.V.G. conducted the metagenomic and viromics analyses, formulated ideas, and wrote the manuscript. T.R.N. and T.L.S. conceived the study, and T.R.N., T.L.S., S.R., and G.T. edited the manuscript. T.L.S. collected samples and extracted DNA. S.R. and G.T. assisted with the generation and analysis of viral sequence data. All of us read and approved the final manuscript.

We declare that we have no competing financial interests.

FOOTNOTES

    • Received 28 August 2019
    • Accepted 7 November 2019
    • Published 17 December 2019

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

REFERENCES

  1. 1.↵
    1. Schlesinger WH,
    2. Reynolds JF,
    3. Cunningham GL,
    4. Huenneke LF,
    5. Jarrell WM,
    6. Virginia RA,
    7. Whitford WG
    . 1990. Biological feedbacks in global desertification. Science 247:1043–1048. doi:10.1126/science.247.4946.1043.
    OpenUrlAbstract/FREE Full Text
  2. 2.↵
    1. Noy-Meir I
    . 1973. Desert ecosystems: environment and producers. Annu Rev Ecol Syst 4:25–51. doi:10.1146/annurev.es.04.110173.000325.
    OpenUrlCrossRef
  3. 3.↵
    1. Reed SC,
    2. Coe KK,
    3. Sparks JP,
    4. Housman DC,
    5. Zelikova TJ,
    6. Belnap J
    . 2012. Changes to dryland rainfall result in rapid moss mortality and altered soil fertility. Nat Clim Chang 2:752–755. doi:10.1038/nclimate1596.
    OpenUrlCrossRef
  4. 4.↵
    1. Rutherford WA,
    2. Painter TH,
    3. Ferrenberg S,
    4. Belnap J,
    5. Okin GS,
    6. Flagg C,
    7. Reed SC
    . 2017. Albedo feedbacks to future climate via climate change impacts on dryland biocrusts. Sci Rep 7:44188. doi:10.1038/srep44188.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Hultman J,
    2. Waldrop MP,
    3. Mackelprang R,
    4. David MM,
    5. McFarland J,
    6. Blazewicz SJ,
    7. Harden J,
    8. Turetsky MR,
    9. McGuire AD,
    10. Shah MB,
    11. VerBerkmoes NC,
    12. Lee LH,
    13. Mavrommatis K,
    14. Jansson JK
    . 2015. Multi-omics of permafrost, active layer and thermokarst bog soil microbiomes. Nature 521:208. doi:10.1038/nature14238.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Belnap J,
    2. Gillette DA
    . 1998. Vulnerability of desert biological soil crusts to wind erosion: the influences of crust development, soil texture, and disturbance. J Arid Environ 39:133–142. doi:10.1006/jare.1998.0388.
    OpenUrlCrossRefWeb of Science
  7. 7.↵
    1. Elbert W,
    2. Weber B,
    3. Burrows S,
    4. Steinkamp J,
    5. Büdel B,
    6. Andreae MO,
    7. Pöschl U
    . 2012. Contribution of cryptogamic covers to the global cycles of carbon and nitrogen. Nat Geosci 5:459. doi:10.1038/ngeo1486.
    OpenUrlCrossRef
  8. 8.↵
    1. Karaoz U,
    2. Couradeau E,
    3. Da Rocha UN,
    4. Lim H-C,
    5. Northen T,
    6. Garcia-Pichel F,
    7. Brodie EL
    . 2018. Large blooms of Bacillales (Firmicutes) underlie the response to wetting of cyanobacterial biocrusts at various stages of maturity. mBio 9:e01366-16. doi:10.1128/mBio.01366-16.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Swenson TL,
    2. Karaoz U,
    3. Swenson JM,
    4. Bowen BP,
    5. Northen TR
    . 2018. Linking soil biology and chemistry in biological soil crust using isolate exometabolomics. Nat Commun 9:19. doi:10.1038/s41467-017-02356-9.
    OpenUrlCrossRef
  10. 10.↵
    1. Schulz F,
    2. Alteio L,
    3. Goudeau D,
    4. Ryan EM,
    5. Feiqiao BY,
    6. Malmstrom RR,
    7. Blanchard J,
    8. Woyke T
    . 2018. Hidden diversity of soil giant viruses. Nat Commun 9:4881. doi:10.1038/s41467-018-07335-2.
    OpenUrlCrossRef
  11. 11.↵
    1. Emerson JB,
    2. Roux S,
    3. Brum JR,
    4. Bolduc B,
    5. Woodcroft BJ,
    6. Jang HB,
    7. Singleton CM,
    8. Solden LM,
    9. Naas AE,
    10. Boyd JA,
    11. Hodgkins SB,
    12. Wilson RM,
    13. Trubl G,
    14. Li C,
    15. Frolking S,
    16. Pope PB,
    17. Wrighton KC,
    18. Crill PM,
    19. Chanton JP,
    20. Saleska SR,
    21. Tyson GW,
    22. Rich VI,
    23. Sullivan MB
    . 2018. Host-linked soil viral ecology along a permafrost thaw gradient. Nat Microbiol 3:870. doi:10.1038/s41564-018-0190-y.
    OpenUrlCrossRef
  12. 12.↵
    1. Trubl G,
    2. Jang HB,
    3. Roux S,
    4. Emerson JB,
    5. Solonenko N,
    6. Vik DR,
    7. Solden L,
    8. Ellenbogen J,
    9. Runyon AT,
    10. Bolduc B,
    11. Woodcroft BJ,
    12. Saleska SR,
    13. Tyson GW,
    14. Wrighton KC,
    15. Sullivan MB,
    16. Rich VI
    . 2018. Soil viruses are underexplored players in ecosystem carbon processing. mSystems 3:e00076-18. doi:10.1128/mSystems.00076-18.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Adriaenssens EM,
    2. Kramer R,
    3. Van Goethem MW,
    4. Makhalanyane TP,
    5. Hogg I,
    6. Cowan DA
    . 2017. Environmental drivers of viral community composition in Antarctic soils identified by viromics. Microbiome 5:83. doi:10.1186/s40168-017-0301-7.
    OpenUrlCrossRef
  14. 14.↵
    1. Zablocki O,
    2. van Zyl L,
    3. Adriaenssens EM,
    4. Rubagotti E,
    5. Tuffin M,
    6. Cary SC,
    7. Cowan D
    . 2014. High-level diversity of tailed phages, eukaryote-associated viruses, and virophage-like elements in the metaviromes of antarctic soils. Appl Environ Microbiol 80:6888–6897. doi:10.1128/AEM.01525-14.
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    1. Paez-Espino D,
    2. Eloe-Fadrosh EA,
    3. Pavlopoulos GA,
    4. Thomas AD,
    5. Huntemann M,
    6. Mikhailova N,
    7. Rubin E,
    8. Ivanova NN,
    9. Kyrpides NC
    . 2016. Uncovering Earth’s virome. Nature 536:425. doi:10.1038/nature19094.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Roux S,
    2. Adriaenssens EM,
    3. Dutilh BE,
    4. Koonin EV,
    5. Kropinski AM,
    6. Krupovic M,
    7. Kuhn JH,
    8. Lavigne R,
    9. Brister JR,
    10. Varsani A,
    11. Amid C,
    12. Aziz RK,
    13. Bordenstein SR,
    14. Bork P,
    15. Breitbart M,
    16. Cochrane GR,
    17. Daly RA,
    18. Desnues C,
    19. Duhaime MB,
    20. Emerson JB,
    21. Enault F,
    22. Fuhrman JA,
    23. Hingamp P,
    24. Hugenholtz P,
    25. Hurwitz BL,
    26. Ivanova NN,
    27. Labonté JM,
    28. Lee K-B,
    29. Malmstrom RR,
    30. Martinez-Garcia M,
    31. Mizrachi IK,
    32. Ogata H,
    33. Páez-Espino D,
    34. Petit M-A,
    35. Putonti C,
    36. Rattei T,
    37. Reyes A,
    38. Rodriguez-Valera F,
    39. Rosario K,
    40. Schriml L,
    41. Schulz F,
    42. Steward GF,
    43. Sullivan MB,
    44. Sunagawa S,
    45. Suttle CA,
    46. Temperton B,
    47. Tringe SG,
    48. Thurber RV,
    49. Webster NS,
    50. Whiteson KL, et al
    . 2019. Minimum information about an uncultivated virus genome (MIUVIG). Nat Biotechnol 37:29. doi:10.1038/nbt.4306.
    OpenUrlCrossRef
  17. 17.↵
    1. Placella SA,
    2. Brodie EL,
    3. Firestone MK
    . 2012. Rainfall-induced carbon dioxide pulses result from sequential resuscitation of phylogenetically clustered microbial groups. Proc Natl Acad Sci U S A 109:10931–10936. doi:10.1073/pnas.1204306109.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Adriaenssens EM,
    2. Van Zyl L,
    3. De Maayer P,
    4. Rubagotti E,
    5. Rybicki E,
    6. Tuffin M,
    7. Cowan DA
    . 2015. Metagenomic analysis of the viral community in N amib D esert hypoliths. Environ Microbiol 17:480–495. doi:10.1111/1462-2920.12528.
    OpenUrlCrossRef
  19. 19.↵
    1. Prestel E,
    2. Regeard C,
    3. Salamitou S,
    4. Neveu J,
    5. DuBow MS
    . 2013. The bacteria and bacteriophages from a Mesquite Flats site of the Death Valley desert. Antonie Van Leeuwenhoek 103:1329–1341. doi:10.1007/s10482-013-9914-4.
    OpenUrlCrossRef
  20. 20.↵
    1. Suttle CA
    . 2007. Marine viruses—major players in the global ecosystem. Nat Rev Microbiol 5:801. doi:10.1038/nrmicro1750.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    1. Lindell D,
    2. Sullivan MB,
    3. Johnson ZI,
    4. Tolonen AC,
    5. Rohwer F,
    6. Chisholm SW
    . 2004. Transfer of photosynthesis genes to and from Prochlorococcus viruses. Proc Natl Acad Sci U S A 101:11013–11018. doi:10.1073/pnas.0401526101.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Sullivan MB,
    2. Coleman ML,
    3. Weigele P,
    4. Rohwer F,
    5. Chisholm SW
    . 2005. Three Prochlorococcus cyanophage genomes: signature features and ecological interpretations. PLoS Biol 3:e144. doi:10.1371/journal.pbio.0030144.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Yoshida M,
    2. Yoshida T,
    3. Kashima A,
    4. Takashima Y,
    5. Hosoda N,
    6. Nagasaki K,
    7. Hiroishi S
    . 2008. Ecological dynamics of the toxic bloom-forming cyanobacterium Microcystis aeruginosa and its cyanophages in freshwater. Appl Environ Microbiol 74:3269–3273. doi:10.1128/AEM.02240-07.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Winter C,
    2. Bouvier T,
    3. Weinbauer MG,
    4. Thingstad TF
    . 2010. Trade-offs between competition and defense specialists among unicellular planktonic organisms: the “killing the winner” hypothesis revisited. Microbiol Mol Biol Rev 74:42–57. doi:10.1128/MMBR.00034-09.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Thingstad T,
    2. Lignell R
    . 1997. Theoretical models for the control of bacterial growth rate, abundance, diversity and carbon demand. Aquat Microb Ecol 13:19–27. doi:10.3354/ame013019.
    OpenUrlCrossRef
  26. 26.↵
    1. Bowers RM,
    2. Kyrpides NC,
    3. Stepanauskas R,
    4. Harmon-Smith M,
    5. Doud D,
    6. Reddy TBK,
    7. Schulz F,
    8. Jarett J,
    9. Rivers AR,
    10. Eloe-Fadrosh EA,
    11. Tringe SG,
    12. Ivanova NN,
    13. Copeland A,
    14. Clum A,
    15. Becraft ED,
    16. Malmstrom RR,
    17. Birren B,
    18. Podar M,
    19. Bork P,
    20. Weinstock GM,
    21. Garrity GM,
    22. Dodsworth JA,
    23. Yooseph S,
    24. Sutton G,
    25. Glöckner FO,
    26. Gilbert JA,
    27. Nelson WC,
    28. Hallam SJ,
    29. Jungbluth SP,
    30. Ettema TJG,
    31. Tighe S,
    32. Konstantinidis KT,
    33. Liu W-T,
    34. Baker BJ,
    35. Rattei T,
    36. Eisen JA,
    37. Hedlund B,
    38. McMahon KD,
    39. Fierer N,
    40. Knight R,
    41. Finn R,
    42. Cochrane G,
    43. Karsch-Mizrachi I,
    44. Tyson GW,
    45. Rinke C,
    46. Lapidus A,
    47. Meyer F,
    48. Yilmaz P,
    49. Parks DH,
    50. Eren AM, et al
    . 2017. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol 35:725. doi:10.1038/nbt.3893.
    OpenUrlCrossRef
  27. 27.↵
    1. Rodriguez-R LM,
    2. Gunturu S,
    3. Tiedje JM,
    4. Cole JR,
    5. Konstantinidis KT
    . 2018. Nonpareil 3: fast estimation of metagenomic coverage and sequence diversity. mSystems 3:e00039-18. doi:10.1128/mSystems.00039-18.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Angel R,
    2. Conrad R
    . 2013. Elucidating the microbial resuscitation cascade in biological soil crusts following a simulated rain event. Environ Microbiol 15:2799–2815. doi:10.1111/1462-2920.12140.
    OpenUrlCrossRefWeb of Science
  29. 29.↵
    1. Belnap J
    . 2006. The potential roles of biological soil crusts in dryland hydrologic cycles. Hydrol Process 20:3159–3178. doi:10.1002/hyp.6325.
    OpenUrlCrossRefWeb of Science
  30. 30.↵
    1. Jang HB,
    2. Bolduc B,
    3. Zablocki O,
    4. Kuhn J,
    5. Roux S,
    6. Adriaenssens E,
    7. Brister JR,
    8. Kropinski A,
    9. Krupovic M,
    10. Turner D
    . 2019. Gene sharing networks to automate genome-based prokaryotic viral taxonomy. BioRxiv doi:10.1101/533240.
    OpenUrlCrossRef
  31. 31.↵
    1. Bolduc B,
    2. Jang HB,
    3. Doulcier G,
    4. You Z-Q,
    5. Roux S,
    6. Sullivan MB
    . 2017. vConTACT: an iVirus tool to classify double-stranded DNA viruses that infect Archaea and Bacteria. PeerJ 5:e3243. doi:10.7717/peerj.3243.
    OpenUrlCrossRef
  32. 32.↵
    1. Brister JR,
    2. Ako-Adjei D,
    3. Bao Y,
    4. Blinkova O
    . 2015. NCBI viral genomes resource. Nucleic Acids Res 43:D571–D577. doi:10.1093/nar/gku1207.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Mahmud RS,
    2. Garifulina K,
    3. Ulyanova V,
    4. Evtugyn V,
    5. Mindubaeva L,
    6. Khazieva L,
    7. Dudkina E,
    8. Vershinina V,
    9. Kolpakov A,
    10. Ilinskaya O
    . 2017. Bacteriophages of soil bacilli: a new multivalent phage of Bacillus altitudinis. Mol Genet Microbiol Virol 32:87–93. doi:10.3103/S0891416817020082.
    OpenUrlCrossRef
  34. 34.↵
    1. Onyenwoke RU,
    2. Brill JA,
    3. Farahi K,
    4. Wiegel J
    . 2004. Sporulation genes in members of the low G+ C Gram-type-positive phylogenetic branch (Firmicutes). Arch Microbiol 182:182–192. doi:10.1007/s00203-004-0696-y.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Galperin MY
    . 2013. Genome diversity of spore-forming Firmicutes. Microbiol Spectr 1:TBS-0015-2012. doi:10.1128/microbiolspectrum.TBS-0015-2012.
    OpenUrlCrossRef
  36. 36.↵
    1. Skennerton CT,
    2. Imelfort M,
    3. Tyson GW
    . 2013. Crass: identification and reconstruction of CRISPR from unassembled metagenomic data. Nucleic Acids Res 41:e105. doi:10.1093/nar/gkt183.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Williamson KE,
    2. Fuhrmann JJ,
    3. Wommack KE,
    4. Radosevich M
    . 2017. Viruses in soil ecosystems: an unknown quantity within an unexplored territory. Annu Rev Virol 4:201–219. doi:10.1146/annurev-virology-101416-041639.
    OpenUrlCrossRef
  38. 38.↵
    1. Rohwer F,
    2. Thurber RV
    . 2009. Viruses manipulate the marine environment. Nature 459:207. doi:10.1038/nature08060.
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    1. Brum JR,
    2. Hurwitz BL,
    3. Schofield O,
    4. Ducklow HW,
    5. Sullivan MB
    . 2016. Seasonal time bombs: dominant temperate viruses affect Southern Ocean microbial dynamics. ISME J 10:437. doi:10.1038/ismej.2015.125.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Ortmann A,
    2. Lawrence J,
    3. Suttle C
    . 2002. Lysogeny and lytic viral production during a bloom of the cyanobacterium Synechococcus spp. Microb Ecol 43:225–231. doi:10.1007/s00248-001-1058-9.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    1. Sharp MD,
    2. Pogliano K
    . 1999. An in vivo membrane fusion assay implicates SpoIIIE in the final stages of engulfment during Bacillus subtilis sporulation. Proc Natl Acad Sci U S A 96:14553–14558. doi:10.1073/pnas.96.25.14553.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Gabiatti N,
    2. Yu P,
    3. Mathieu J,
    4. Lu GW,
    5. Wang X,
    6. Zhang H,
    7. Soares HM,
    8. Alvarez PJ
    . 2018. Bacterial endospores as phage genome carriers and protective shells. Appl Environ Microbiol 84:e01186-18. doi:10.1128/AEM.01186-18.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. Felde V,
    2. Peth S,
    3. Uteau-Puschmann D,
    4. Drahorad S,
    5. Felix-Henningsen P
    . 2014. Soil microstructure as an under-explored feature of biological soil crust hydrological properties: case study from the NW Negev Desert. Biodivers Conserv 23:1687–1708. doi:10.1007/s10531-014-0693-7.
    OpenUrlCrossRef
  44. 44.↵
    1. Lebre PH,
    2. De Maayer P,
    3. Cowan DA
    . 2017. Xerotolerant bacteria: surviving through a dry spell. Nat Rev Microbiol 15:285. doi:10.1038/nrmicro.2017.16.
    OpenUrlCrossRef
  45. 45.↵
    1. Pointing SB,
    2. Belnap J
    . 2012. Microbial colonization and controls in dryland systems. Nat Rev Microbiol 10:551. doi:10.1038/nrmicro2831.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Griffin JS,
    2. Wells GF
    . 2017. Regional synchrony in full-scale activated sludge bioreactors due to deterministic microbial community assembly. ISME J 11:500. doi:10.1038/ismej.2016.121.
    OpenUrlCrossRef
  47. 47.↵
    1. Zablocki O,
    2. Adriaenssens EM,
    3. Cowan D
    . 2016. Diversity and ecology of viruses in hyperarid desert soils. Appl Environ Microbiol 82:770–777. doi:10.1128/AEM.02651-15.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Garcia-Pichel F,
    2. Loza V,
    3. Marusenko Y,
    4. Mateo P,
    5. Potrafka RM
    . 2013. Temperature drives the continental-scale distribution of key microbes in topsoil communities. Science 340:1574–1577. doi:10.1126/science.1236404.
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    1. Schmieder R,
    2. Edwards R
    . 2011. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27:863–864. doi:10.1093/bioinformatics/btr026.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    1. Li D,
    2. Liu C-M,
    3. Luo R,
    4. Sadakane K,
    5. Lam T-W
    . 2015. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31:1674–1676. doi:10.1093/bioinformatics/btv033.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Van der Walt AJ,
    2. Van Goethem MW,
    3. Ramond J-B,
    4. Makhalanyane TP,
    5. Reva O,
    6. Cowan DA
    . 2017. Assembling metagenomes, one community at a time. BMC Genomics 18:521. doi:10.1186/s12864-017-3918-9.
    OpenUrlCrossRef
  52. 52.↵
    1. Mikheenko A,
    2. Saveliev V,
    3. Gurevich A
    . 2016. MetaQUAST: evaluation of metagenome assemblies. Bioinformatics 32:1088–1090. doi:10.1093/bioinformatics/btv697.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Hyatt D,
    2. Chen G-L,
    3. LoCascio PF,
    4. Land ML,
    5. Larimer FW,
    6. Hauser LJ
    . 2010. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11:119. doi:10.1186/1471-2105-11-119.
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. Buchfink B,
    2. Xie C,
    3. Huson DH
    . 2015. Fast and sensitive protein alignment using DIAMOND. Nat Methods 12:59. doi:10.1038/nmeth.3176.
    OpenUrlCrossRefPubMed
  55. 55.↵
    1. Huson DH,
    2. Auch AF,
    3. Qi J,
    4. Schuster SC
    . 2007. MEGAN analysis of metagenomic data. Genome Res 17:377–386. doi:10.1101/gr.5969107.
    OpenUrlAbstract/FREE Full Text
  56. 56.↵
    1. Kopylova E,
    2. Noé L,
    3. Touzet H
    . 2012. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics 28:3211–3217. doi:10.1093/bioinformatics/bts611.
    OpenUrlCrossRefPubMedWeb of Science
  57. 57.↵
    1. Callahan BJ,
    2. McMurdie PJ,
    3. Rosen MJ,
    4. Han AW,
    5. Johnson AJA,
    6. Holmes SP
    . 2016. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods 13:581. doi:10.1038/nmeth.3869.
    OpenUrlCrossRefPubMed
  58. 58.↵
    1. Quast C,
    2. Pruesse E,
    3. Yilmaz P,
    4. Gerken J,
    5. Schweer T,
    6. Yarza P,
    7. Peplies J,
    8. Glöckner FO
    . 2013. The SILVA ribosomal RNA gene database project: improved data processing and Web-based tools. Nucleic Acids Res 41:D590–D596. doi:10.1093/nar/gks1219.
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.↵
    1. Friedman J,
    2. Alm EJ
    . 2012. Inferring correlation networks from genomic survey data. PLoS Comput Biol 8:e1002687. doi:10.1371/journal.pcbi.1002687.
    OpenUrlCrossRefPubMed
  60. 60.↵
    1. Shannon P,
    2. Markiel A,
    3. Ozier O,
    4. Baliga NS,
    5. Wang JT,
    6. Ramage D,
    7. Amin N,
    8. Schwikowski B,
    9. Ideker T
    . 2003. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 13:2498–2504. doi:10.1101/gr.1239303.
    OpenUrlAbstract/FREE Full Text
  61. 61.↵
    1. Bader GD,
    2. Hogue CW
    . 2003. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics 4:2. doi:10.1186/1471-2105-4-2.
    OpenUrlCrossRefPubMed
  62. 62.↵
    1. Kang D,
    2. Li F,
    3. Kirton ES,
    4. Thomas A,
    5. Egan RS,
    6. An H,
    7. Wang Z
    . 2019. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ Preprints 7:e27522v1. doi:10.7717/peerj.7359.
    OpenUrlCrossRef
  63. 63.↵
    1. Parks DH,
    2. Imelfort M,
    3. Skennerton CT,
    4. Hugenholtz P,
    5. Tyson GW
    . 2015. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res 25:1043–1055. doi:10.1101/gr.186072.114.
    OpenUrlAbstract/FREE Full Text
  64. 64.↵
    1. Seemann T
    . 2014. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30:2068–2069. doi:10.1093/bioinformatics/btu153.
    OpenUrlCrossRefPubMedWeb of Science
  65. 65.↵
    1. Aziz RK,
    2. Bartels D,
    3. Best AA,
    4. DeJongh M,
    5. Disz T,
    6. Edwards RA,
    7. Formsma K,
    8. Gerdes S,
    9. Glass EM,
    10. Kubal M,
    11. Meyer F,
    12. Olsen GJ,
    13. Olson R,
    14. Osterman AL,
    15. Overbeek RA,
    16. McNeil LK,
    17. Paarmann D,
    18. Paczian T,
    19. Parrello B,
    20. Pusch GD,
    21. Reich C,
    22. Stevens R,
    23. Vassieva O,
    24. Vonstein V,
    25. Wilke A,
    26. Zagnitko O
    . 2008. The RAST Server: rapid annotations using subsystems technology. BMC Genomics 9:75. doi:10.1186/1471-2164-9-75.
    OpenUrlCrossRefPubMed
  66. 66.↵
    1. Okonechnikov K,
    2. Golosova O,
    3. Fursov M
    , UGENE team. 2012. Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics 28:1166–1167. doi:10.1093/bioinformatics/bts091.
    OpenUrlCrossRefPubMedWeb of Science
  67. 67.↵
    1. Katoh K,
    2. Standley DM
    . 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30:772–780. doi:10.1093/molbev/mst010.
    OpenUrlCrossRefPubMedWeb of Science
  68. 68.↵
    1. Letunic I,
    2. Bork P
    . 2007. Interactive Tree Of Life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics 23:127–128. doi:10.1093/bioinformatics/btl529.
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.↵
    1. Rodriguez-R LM,
    2. Gunturu S,
    3. Harvey WT,
    4. Rosselló-Mora R,
    5. Tiedje JM,
    6. Cole JR,
    7. Konstantinidis KT
    . 2018. The Microbial Genomes Atlas (MiGA) webserver: taxonomic and gene diversity analysis of Archaea and Bacteria at the whole genome level. Nucleic Acids Res 46:W282–W288. doi:10.1093/nar/gky467.
    OpenUrlCrossRef
  70. 70.↵
    1. Parks DH,
    2. Tyson GW,
    3. Hugenholtz P,
    4. Beiko RG
    . 2014. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics 30:3123–3124. doi:10.1093/bioinformatics/btu494.
    OpenUrlCrossRefPubMed
  71. 71.↵
    1. Oksanen J,
    2. Kindt R,
    3. Legendre P,
    4. O’Hara B,
    5. Stevens MHH,
    6. Oksanen MJ
    . 2007. The vegan package: community ecology package. https://cran.r-project.org/package=vegan.
  72. 72.↵
    1. Roux S,
    2. Enault F,
    3. Hurwitz BL,
    4. Sullivan MB
    . 2015. VirSorter: mining viral signal from microbial genomic data. PeerJ 3:e985. doi:10.7717/peerj.985.
    OpenUrlCrossRef
  73. 73.↵
    1. Roux S,
    2. Brum JR,
    3. Dutilh BE,
    4. Sunagawa S,
    5. Duhaime MB,
    6. Loy A,
    7. Poulos BT,
    8. Solonenko N,
    9. Lara E,
    10. Poulain J,
    11. Pesant S,
    12. Kandels-Lewis S,
    13. Dimier C,
    14. Picheral M,
    15. Searson S,
    16. Cruaud C,
    17. Alberti A,
    18. Duarte CM,
    19. Gasol JM,
    20. Vaqué D,
    21. Bork P,
    22. Acinas SG,
    23. Wincker P,
    24. Sullivan MB
    . 2016. Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses. Nature 537:689. doi:10.1038/nature19366.
    OpenUrlCrossRefPubMed
  74. 74.↵
    1. Langmead B,
    2. Salzberg SL
    . 2012. Fast gapped-read alignment with Bowtie 2. Nat Methods 9:357. doi:10.1038/nmeth.1923.
    OpenUrlCrossRefPubMedWeb of Science
  75. 75.↵
    1. Zhang Y
    . 2008. I-TASSER server for protein 3D structure prediction. BMC Bioinformatics 9:40. doi:10.1186/1471-2105-9-40.
    OpenUrlCrossRefPubMed
  76. 76.↵
    1. Zhang Y,
    2. Skolnick J
    . 2005. TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res 33:2302–2309. doi:10.1093/nar/gki524.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
Characteristics of Wetting-Induced Bacteriophage Blooms in Biological Soil Crust
Marc W. Van Goethem, Tami L. Swenson, Gareth Trubl, Simon Roux, Trent R. Northen
mBio Dec 2019, 10 (6) e02287-19; DOI: 10.1128/mBio.02287-19

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mBio article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Characteristics of Wetting-Induced Bacteriophage Blooms in Biological Soil Crust
(Your Name) has forwarded a page to you from mBio
(Your Name) thought you would be interested in this article in mBio.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Characteristics of Wetting-Induced Bacteriophage Blooms in Biological Soil Crust
Marc W. Van Goethem, Tami L. Swenson, Gareth Trubl, Simon Roux, Trent R. Northen
mBio Dec 2019, 10 (6) e02287-19; DOI: 10.1128/mBio.02287-19
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Firmicutes
bacteriophages
biocrust
soil microbiology
virus-host interactions
wetting

Related Articles

Cited By...

About

  • About mBio
  • Editor in Chief
  • Board of Editors
  • AAM Fellows
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Article Types
  • Ethics
  • Contact Us

Follow #mBio

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2150-7511