Designer viral receptors

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.

Plant genotype is recognized to contribute to variations in microbial community structure in the rhizosphere, soil adherent to roots. However, the extent to which the viral community varies has remained poorly understood and has the potential to contribute to variation in soil microbial communities. Here we cultivated replicates of two Zea mays genotypes, parviglumis and B73, in a greenhouse and harvested the rhizobiome (rhizoplane and rhizosphere) to identify the abundance of cells and viruses as well as rhizobiome microbial and viral community using 16S rRNA gene amplicon sequencing and genome resolved metagenomics. Our results demonstrated that viruses exceeded microbial abundance in the rhizobiome of parviglumis and B73 with a significant variation in both the microbial and viral community between the two genotypes. Of the viral contigs identified only 4.5% (n = 7) of total viral contigs were shared between the two genotypes, demonstrating that plants even at the level of genotype can significantly alter the surrounding soil viral community. An auxiliary metabolic gene associated with glycoside hydrolase (GH5) degradation was identified in one viral metagenome-assembled genome (vOTU) identified in the B73 rhizobiome infecting Propionibacteriaceae (Actinobacteriota) further demonstrating the viral contribution in metabolic potential for carbohydrate degradation and carbon cycling in the rhizosphere. This variation demonstrates the potential of plant genotype to contribute to microbial and viral heterogeneity in soil systems and harbors genes capable of contributing to carbon cycling in the rhizosphere.