In marine ecosystems, viruses exert control on the composition and metabolism of microbial communities, influencing overall biogeochemical cycling. Deep sea sediments associated with cold seeps are known to host taxonomically diverse microbial communities, but little is known about viruses infecting these microorganisms. Here, we probed metagenomes from seven geographically diverse cold seeps across global oceans to assess viral diversity, virus–host interaction, and virus-encoded auxiliary metabolic genes (AMGs). Gene-sharing network comparisons with viruses inhabiting other ecosystems reveal that cold seep sediments harbour considerable unexplored viral diversity. Most cold seep viruses display high degrees of endemism with seep fluid flux being one of the main drivers of viral community composition. In silico predictions linked 14.2% of the viruses to microbial host populations with many belonging to poorly understood candidate bacterial and archaeal phyla. Lysis was predicted to be a predominant viral lifestyle based on lineage-specific virus/host abundance ratios. Metabolic predictions of prokaryotic host genomes and viral AMGs suggest that viruses influence microbial hydrocarbon biodegradation at cold seeps, as well as other carbon, sulfur and nitrogen cycling via virus-induced mortality and/or metabolic augmentation. Overall, these findings reveal the global diversity and biogeography of cold seep viruses and indicate how viruses may manipulate seep microbial ecology and biogeochemistry.

In marine ecosystems, viruses exert control on the composition and metabolism of microbial communities, thus influencing overall biogeochemical cycling. Deep sea sediments associated with cold seeps are known to host taxonomically diverse microbial communities, but little is known about viruses infecting these microorganisms. Here, we probed metagenomes from seven geographically diverse cold seeps across global oceans, to assess viral diversity, virus-host interaction, and virus-encoded auxiliary metabolic genes (AMGs). Gene-sharing network comparisons with viruses inhabiting other ecosystems reveal that cold seep sediments harbour considerable unexplored viral diversity. Most cold seep viruses display high degrees of endemism with seep fluid flux being one of the main drivers of viral community composition. In silico predictions linked 14.2% of the viruses to microbial host populations, with many belonging to poorly understood candidate bacterial and archaeal phyla. Lysis was predicted to be a predominant viral lifestyle based on lineage-specific virus/host abundance ratios. Metabolic predictions of prokaryotic host genomes and viral AMGs suggest that viruses influence microbial hydrocarbon biodegradation at cold seeps, as well as other carbon, sulfur and nitrogen cycling via virus-induced mortality and/or metabolic augmentation. Overall, these findings reveal the global diversity and biogeography of cold seep viruses and indicate how viruses may manipulate seep microbial ecology and biogeochemistry.

Microbially mediated nitrogen cycling in carbon-dominated cold seep environments remains poorly understood. So far anaerobic methanotrophic archaea (ANME-2) and their sulfate-reducing bacterial partners (SEEP-SRB1 clade) have been identified as diazotrophs in deep sea cold seep sediments. However, it is unclear whether other microbial groups can perform nitrogen fixation in such ecosystems. To fill this gap, we analyzed 61 metagenomes, 1428 metagenome-assembled genomes, and six metatranscriptomes derived from 11 globally distributed cold seeps. These sediments contain phylogenetically diverse nitrogenase genes corresponding to an expanded diversity of diazotrophic lineages. Diverse catabolic pathways were predicted to provide ATP for nitrogen fixation, suggesting diazotrophy in cold seeps is not necessarily associated with sulfate-dependent anaerobic oxidation of methane. Nitrogen fixation genes among various diazotrophic groups in cold seeps were inferred to be genetically mobile and subject to purifying selection. Our findings extend the capacity for diazotrophy to five candidate phyla (Altarchaeia, Omnitrophota, FCPU426, Caldatribacteriota and UBA6262), and suggest that cold seep diazotrophs might contribute substantially to the global nitrogen balance.