Dispersal, habitat filtering, and eco-evolutionary dynamics as drivers of local and global wetland viral biogeography.

Wetlands store 20-30% of the world's soil carbon, and identifying the microbial controls on these carbon reserves is essential to predicting feedbacks to climate change. Although viral infections likely play important roles in wetland ecosystem dynamics, we lack a basic understanding of wetland viral ecology. Here 63 viral size-fraction metagenomes (viromes) and paired total metagenomes were generated from three time points in 2021 at seven fresh- and saltwater wetlands in the California Bodega Marine Reserve. We recovered 12,826 viral population genomic sequences (vOTUs), only 4.4% of which were detected at the same field site two years prior, indicating a small degree of population stability or recurrence. Viral communities differed most significantly among the seven wetland sites and were also structured by habitat (plant community composition and salinity). Read mapping to a new version of our reference database, PIGEONv2.0 (515,763 vOTUs), revealed 196 vOTUs present over large geographic distances, often reflecting shared habitat characteristics. Wetland vOTU microdiversity was significantly lower locally than globally and lower within than between time points, indicating greater divergence with increasing spatiotemporal distance. Viruses tended to have broad predicted host ranges via CRISPR spacer linkages to metagenome-assembled genomes, and increased SNP frequencies in CRISPR-targeted major tail protein genes suggest potential viral eco-evolutionary dynamics in response to both immune targeting and changes in host cell receptors involved in viral attachment. Together, these results highlight the importance of dispersal, environmental selection, and eco-evolutionary dynamics as drivers of local and global wetland viral biogeography.

High-quality, chromosome-scale genome assemblies: comparisons of three Diaphorina citri (Asian citrus psyllid) geographic populations.

The Asian citrus psyllid, Diaphorina citri, is the insect vector of the causal agent of huanglongbing (HLB), a devastating bacterial disease of commercial citrus. Presently, few genomic resources exist for D. citri. In this study, we utilized PacBio HiFi and chromatin confirmation contact (Hi-C) sequencing to sequence, assemble, and compare three high-quality, chromosome-scale genome assemblies of D. citri collected from California, Taiwan, and Uruguay. Our assemblies had final sizes of 282.67 Mb (California), 282.89 Mb (Taiwan), and 266.67 Mb (Uruguay) assembled into 13 pseudomolecules-a reduction in assembly size of 41-45% compared with previous assemblies which we validated using flow cytometry. We identified the X chromosome in D. citri and annotated each assembly for repetitive elements, protein-coding genes, transfer RNAs, ribosomal RNAs, piwi-interacting RNA clusters, and endogenous viral elements. Between 19,083 and 20,357 protein-coding genes were predicted. Repetitive DNA accounts for 36.87-38.26% of each assembly. Comparative analyses and mitochondrial haplotype networks suggest that Taiwan and Uruguay D. citri are more closely related, while California D. citri are closely related to Florida D. citri. These high-quality, chromosome-scale assemblies provide new genomic resources to researchers to further D. citri and HLB research.

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

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

Minnesota peat viromes reveal terrestrial and aquatic niche partitioning for local and global viral populations.

BACKGROUND: Peatlands are expected to experience sustained yet fluctuating higher temperatures due to climate change, leading to increased microbial activity and greenhouse gas emissions. Despite mounting evidence for viral contributions to these processes in peatlands underlain with permafrost, little is known about viruses in other peatlands. More generally, soil viral biogeography and its potential drivers are poorly understood at both local and global scales. Here, 87 metagenomes and five viral size-fraction metagenomes (viromes) from a boreal peatland in northern Minnesota (the SPRUCE whole-ecosystem warming experiment and surrounding bog) were analyzed for dsDNA viral community ecological patterns, and the recovered viral populations (vOTUs) were compared with our curated PIGEON database of 266,125 vOTUs from diverse ecosystems. RESULTS: Within the SPRUCE experiment, viral community composition was significantly correlated with peat depth, water content, and carbon chemistry, including CH4 and CO2 concentrations, but not with temperature during the first 2 years of warming treatments. Peat vOTUs with aquatic-like signatures (shared predicted protein content with marine and/or freshwater vOTUs) were significantly enriched in more waterlogged surface peat depths. Predicted host ranges for SPRUCE vOTUs were relatively narrow, generally within a single bacterial genus. Of the 4326 SPRUCE vOTUs, 164 were previously detected in other soils, mostly peatlands. None of the previously identified 202,371 marine and freshwater vOTUs in our PIGEON database were detected in SPRUCE peat, but 0.4% of 80,714 viral clusters (VCs, grouped by predicted protein content) were shared between soil and aquatic environments. On a per-sample basis, vOTU recovery was 32 times higher from viromes compared with total metagenomes. CONCLUSIONS: Results suggest strong viral "species" boundaries between terrestrial and aquatic ecosystems and to some extent between peat and other soils, with differences less pronounced at higher taxonomic levels. The significant enrichment of aquatic-like vOTUs in more waterlogged peat suggests that viruses may also exhibit niche partitioning on more local scales. These patterns are presumably driven in part by host ecology, consistent with the predicted narrow host ranges. Although more samples and increased sequencing depth improved vOTU recovery from total metagenomes, the substantially higher per-sample vOTU recovery after viral particle enrichment highlights the utility of soil viromics. Video abstract The importance of Minnesota peat viromes in revealing terrestrial and aquatic niche partitioning for viral populations.

DNase Treatment Improves Viral Enrichment in Agricultural Soil Viromes.

The small genomes of most viruses make it difficult to fully capture viral diversity in metagenomes dominated by DNA from cellular organisms. Viral size fraction metagenomics (viromics) protocols facilitate the enrichment of viral DNA from environmental samples, and these protocols typically include DNase treatment of the post-0.2-µm-filtered viromic fraction to remove contaminating free DNA prior to virion lysis. However, DNase may also remove desirable viral genomic DNA (e.g., contained in virions compromised due to frozen storage or laboratory processing), suggesting that DNase-untreated viromes might be useful in some cases. In order to understand how virome preparation with and without DNase treatment influences the resultant data, here, we compared 15 soil viromes (7 DNase treated and 8 untreated) from 8 samples collected from agricultural fields prior to tomato planting. DNase-treated viromes yielded significantly more assembled viral contigs, contained significantly less nonviral microbial DNA, and recovered more viral populations (viral operational taxonomic units [vOTUs]) through read mapping. However, DNase-treated and untreated viromes were statistically indistinguishable in terms of ecological patterns across viral communities. Although the results suggest that DNase treatment is preferable where possible, in comparison to previously reported total metagenomes from the same samples, both DNase-treated and untreated viromes were significantly enriched in viral signatures by all metrics compared, including a 225-times-higher proportion of viral reads in untreated viromes compared to total metagenomes. Thus, even without DNase treatment, viromics was preferable to total metagenomics for capturing viral diversity in these soils, suggesting that preparation of DNase-untreated viromes can be worthwhile when DNase treatment is not possible. IMPORTANCE Viromics is becoming an increasingly popular method for characterizing soil viral communities. DNase treatment of the viral size fraction prior to DNA extraction is meant to reduce contaminating free DNA and is a common step within viromics protocols to ensure that sequences are of viral origin. However, some samples may not be amenable to DNase treatment due to viral particles being compromised either in storage (i.e., frozen) or during other sample processing steps. To date, the effect of DNase treatment on the recovery of viruses and downstream ecological interpretations of soil viral communities is not thoroughly understood. This work sheds light on these questions and indicates that while DNase treatment of soil viromes improves the recovery of viral populations, this improvement is modest in comparison to the gains made by viromics over total soil metagenomics. Furthermore, DNase treatment may not be necessary to observe the ecological patterns structuring soil viral communities.

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

Viruses are abundant yet understudied members of soil environments that influence terrestrial biogeochemical cycles. Here, we characterized the dsDNA viral diversity in biochar-amended agricultural soils at the preplanting and harvesting stages of a tomato growing season via paired total metagenomes and viral size fraction metagenomes (viromes). Size fractionation prior to DNA extraction reduced sources of nonviral DNA in viromes, enabling the recovery of a vaster richness of viral populations (vOTUs), greater viral taxonomic diversity, broader range of predicted hosts, and better access to the rare virosphere, relative to total metagenomes, which tended to recover only the most persistent and abundant vOTUs. Of 2961 detected vOTUs, 2684 were recovered exclusively from viromes, while only three were recovered from total metagenomes alone. Both viral and microbial communities differed significantly over time, suggesting a coupled response to rhizosphere recruitment processes and/or nitrogen amendments. Viral communities alone were also structured along an 18 m spatial gradient. Overall, our results highlight the utility of soil viromics and reveal similarities between viral and microbial community dynamics throughout the tomato growing season yet suggest a partial decoupling of the processes driving their spatial distributions, potentially due to differences in dispersal, decay rates, and/or sensitivities to soil heterogeneity.

Endogenous Viral Elements Are Widespread in Arthropod Genomes and Commonly Give Rise to PIWI-Interacting RNAs.

Arthropod genomes contain sequences derived from integrations of DNA and nonretroviral RNA viruses. These sequences, known as endogenous viral elements (EVEs), have been acquired over the course of evolution and have been proposed to serve as a record of past viral infections. Recent evidence indicates that EVEs can function as templates for the biogenesis of PIWI-interacting RNAs (piRNAs) in some mosquito species and cell lines, raising the possibility that EVEs may serve as a source of immunological memory in these organisms. However, whether piRNAs are derived from EVEs or serve an antiviral function in other arthropod species is unknown. Here, we used publicly available genome assemblies and small RNA sequencing data sets to characterize the repertoire and function of EVEs across 48 arthropod genomes. We found that EVEs are widespread in arthropod genomes and primarily correspond to unclassified single-stranded RNA (ssRNA) viruses and viruses belonging to the Rhabdoviridae and Parvoviridae families. Additionally, EVEs were enriched in piRNA clusters in a majority of species, and we found that production of primary piRNAs from EVEs is common, particularly for EVEs located within piRNA clusters. While the abundance of EVEs within arthropod genomes and the frequency with which EVEs give rise to primary piRNAs generally support the hypothesis that EVEs contribute to an antiviral response via the piRNA pathway, limited nucleotide identity between currently described viruses and EVEs identified here likely limits the extent to which this process plays a role during infection with known viruses in the arthropod species analyzed.IMPORTANCE Our results greatly expand the knowledge of EVE abundance, diversity, and function in an exceptionally wide range of arthropod species. We found that while previous findings in mosquitoes regarding the potential of EVEs to serve as sources of immunological memory via the piRNA pathway may be generalized to other arthropod species, speculation regarding the antiviral function of EVE-derived piRNAs should take into context the fact that EVEs are, in the vast majority of cases, not similar enough to currently described viruses at the nucleotide level to serve as sources of antiviral piRNAs against them.