Ecological drivers and potential functions of viral communities in flooded arsenic-contaminated paddy soils

This work revealed the profile of viral communities in paddy soils with different levels of arsenic (As) contamination during the flooded period. The structure of viral communities differed significantly in highly and moderately As-contaminated soils. The diversity of soil viral communities under high As contamination decreased. Siphoviridae, Podoviridae, Myoviridae, and Microviridae were the dominant viral families in all samples, and the relative abundances of five of the top 20 viral genera were significantly different between highly and moderately As-contaminated groups. Seventeen dissimilatory As(V)-reducing bacteria were predicted to host 161 viral operational taxonomic units (vOTUs), mainly affiliated with the genera of Sulfurospirillum, Deferribacter, Bacillus and Fusibacter. Among them, 28 vOTUs were also associated with Fe(III)-reducing bacteria, which belonged to different species of the genus Shewanella. Procrustes analysis showed that the community structure of soil viruses was strongly correlated with both prokaryotic community structure and geochemical properties. Random forest analyses revealed that the Total-Fe, DCB-Fe and oxalate-Fe were the most significant variables on viral community richness, while the total-As concentration was an important factor on the Shannon index. Furthermore, As resistance genes (ArsC, ArsR and ArsD), As methylation genes (arsM) and As transporter genes (Pst and Pit) were identified among the auxiliary metabolic genes (AMGs) of the virome. This work revealed that the viruses might influence microbial adaptation in response to As-induced stress, and provided a perspective on the potential virus-mediated biogeochemical cycling of As.