phoH-carrying virus communities responded to multiple factors and their correlation network with prokaryotes in sediments along Bohai Sea, Yellow Sea, and East China Sea in China.

Viruses carrying phoH genes are an important functional group that may boost phosphate metabolism of their prokaryote hosts and affect phosphorus cycle in the ocean. However, at present, very little is known about the phoH-carrying viruses' community structure and diversity in marine sediments, as well as their correlation network with prokaryotes and environment. Here, via a large spatial scale investigation along the Bohai Sea, Yellow Sea, and East China Sea, for the first time, diverse unknown benthic phoH-carrying viruses were uncovered, which were mainly affiliated to three clusters. Interestingly, these viruses presented a very distinct community structure compared to those in seawaters. Correlation network analysis implied that these viruses might mainly infect the prokaryotes of Gamm-/Delta-proteobacteria, Thaumarchaeota, and Cyanobacteria in sediments. Distinct virus-prokaryote correlation network modules were shown in different sea areas. These modules' highly nested feature implied their coevolution with prokaryotes during long-term arms race. Their distribution in sediments was influenced by multiple factors including geographic separation and the key environmental variables of total organic carbon and total phosphorus, and responded to terrestrial inputs and coastal aquaculture activities. The results of this study provide novel insights into the benthic virus communities potentially participating in phosphorus cycling in the ocean.