Total metagenomes outperform viromes in recovering viral diversity from sulfuric soils

Abstract Recent metagenomic advancements have offered unprecedented insights into soil viral ecology. However, it remains a challenge to select the suitable metagenomic method for investigating soil viruses under different environmental conditions. Here, we assessed the performance of viral size-fraction metagenomes (viromes) and total metagenomes in capturing viral diversity from hypersulfidic soils with neutral pH and sulfuric soils with pH <3.3. Viromes effectively enhanced the sequencing coverage of viral genomes in both soil types. Viomes of hypersulfidic soils outperformed total metagenomes by recovering a significantly higher number of viral operational taxonomic units (vOTUs). However, total metagenomes of sulfuric soils recovered ~4.5 times more vOTUs than viromes on average. Altogether, our findings suggest that the choice between viromes and total metagenomes for studying soil viruses should be carefully considered based on the specific environmental conditions.

Viruses are highly abundant on Earth and contribute significantly to biogeochemical cycles.Total metagenomes, extracting all microbial DNA, have been used to study the soil virosphere [1,2].However, the vast diversity of soil microbiota and "relic" DNA poses obstacles to the reconstruction of viral genomes using total metagenomes [3,4].To address these limitations, researchers have turned to soil viral size-fraction metagenomes (viromes) for reducing cellular contamination and enhancing viral particle concentration [5].Several studies demonstrated that viromes outperformed total metagenomes in uncovering rare virosphere in agriculture and peatland soils [6,7].However, viromes also have drawbacks, including low DNA yields, extraction of only free viral particles, and labour-intensive procedures [5].The broader applicability of viromes remains relatively unexplored due to significant variations in soil biophysiochemical properties across different types [8].Therefore, it becomes imperative to evaluate the effectiveness of viromes versus total metagenomes across diverse soil types for accurately understanding soil viral dynamics.
Here, we collected 36 samples from Adelaide, South Australia, including two types of acid sulfate soils (18 hypersulfidic soil samples with neutral pH from a current mangrove swamp on the Garden Island and 18 sulfuric soil samples with pH <3.3 from adjacent disturbed areas).Since 1935, the disturbed area has experienced reclamation for agriculture and industry, leading to the oxidation of sulfidic materials and generation of substantial acid and sulfuric soils [9,10].We constructed both virome and total metagenome libraries to retrieve the viral diversity, following established protocols [7,11] (see Supplementary materials).
Viromes of hypersulfidic soils revealed the highest number of vOTUs (14094), followed by total metagenomes of hypersulfidic soils (4151) and sulfuric soils (1360), and viromes of sulfuric soils (430) (Fig. 1C).The accumulation curves for vOTUs reached saturation (Fig. S1), suggesting that our sampling effort was adequate to capture the richness of vOTUs.In hypersulfidic soils, each virome yielded 5.1 times more vOTUs compared to total metagenomes.Conversely, in sulfuric soils, the viral richness per total metagenome exceeded viromes by 4.5 times (Fig. 1E).This aligns with observations from human gut viral metadata, where total metagenomes identified more viral contigs than viromes Metagenomic sequencing of viral communities | 3 [12].Furthermore, vOTUs were grouped into viral clusters (VCs) using vConTACT2 (Figs 1F and S2) [13].The number of VCs clustered from the viromes of hypersulfidic was higher than that from total metagenomes, while sulfuric soils showed an opposite trend (Fig. 1F).These results suggest that viromes outperform total metagenomes in capturing higher viral diversity from hypersulfidic soils, whereas sulfuric soils, considered as extreme environments, showed an opposite pattern.
The better performance of total metagenomes over viromes in sulfuric soils may be mainly attributed to the prolonged extreme acidity of the environment, which may directly impact the persistence of most free viral particles [14,15].Additionally, previous studies showed that heavy metals can significantly elevate lysogenic virus abundance [16], and in some extreme environments, one or more viruses were present in nearly every cell [17].In sulfuric soils, harsh conditions may induce the shift from lytic to lysogenic viral strategies for survival adaption.This is partially supported by the significantly higher proportion of lysogenic vOTUs in the total metagenome of sulfuric soils compared to hypersulfidic soils (Fig. 1G).Those factors collectively render viromes less effective in capturing prevalent viral signals within this environment (Table S1).Moreover, technical factors can also inf luence the performance of different metagenomes in capturing the diversity of soil viruses.For instance, the steps involved in virome DNA extraction and library preparation may lead to the reduction of viral diversity in sulfuric soils.Additionally, the lower biological complexity in extreme environments, such as sulfuric soils, may contribute to an enhanced efficiency of vOTU recovery from total metagenomes [5,18].
Viromes have been widely recognized as a powerful tool for studying soil viral ecology.Our results support viromes as a promising approach for studying viruses in natural soil systems.However, in extreme environments such as sulfuric soils, total metagenomes can be the initial consideration due to their efficiency and lower labour requirements in capturing a comprehensive viral profile.Moreover, future studies should encompass a broad array of soil types to refine the application of metagenomic methods in characterizing the soil virosphere.

Figure 1 .
Figure 1.Comparisons of the performance of viromes and total metagenomes in recovering viral diversity from hypersulfidic and sulfuric soils.(A) The number of all contigs assembled and viral contigs identified from viromes and total metagenomes of hypersulfidic and sulfuric soils.(B) The average percentage of contigs belonging to viral genomes, determined by per virome and total metagenome from both soil types.(C) Venn diagram showing the shared and unique vOTUs among the two soil types with the two metagenomic methods.(D) The percentage of reads mapped to vOTUs, assessed by each virome and total metagenome from hypersulfidic and sulfuric soils.(E) The number of vOTUs recovered per virome and total metagenome based on the reads mapping to all vOTUs.(F) The size of viral clusters from viromes or total metagenomes in hypersulfidic and sulfuric soils.(G) Percentage of putative lysogenic viruses in total metagenomes of the two soil types.Different letters above the boxes show the significant difference determined by the Kruskal-Wallis nonparametric test.