# Virus diversity and interactions with hosts in deep-sea hydrothermal vents

Ruolin Cheng, Xiaofeng Li, Lijing Jiang, Linfeng Gong, Claire Geslin, Zongze

Shao

► To cite this version:

Ruolin Cheng, Xiaofeng Li, Lijing Jiang, Linfeng Gong, Claire Geslin, et al. Virus diversity and interactions with hosts in deep-sea hydrothermal vents. Microbiome, 2022, 10 (1), pp.235. 10.1186/s40168-022-01441-6. hal-04729795

HAL Id: hal-04729795

https://hal.science/hal-04729795v1

Submitted on 15 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# Virus diversity and interactions with hosts in deep-sea hydrothermal vents

Ruolin Cheng<sup>1,2</sup>, Xiaofeng Li<sup>1,3</sup>, Lijing Jiang<sup>1</sup>, Linfeng Gong<sup>1</sup>, Claire Geslin<sup>4,5</sup> and Zongze Shao<sup>1,5,6*</sup>

## Abstract

**Background:** The deep sea harbors many viruses, yet their diversity and interactions with hosts in hydrothermal ecosystems are largely unknown. Here, we analyzed the viral composition, distribution, host preference, and metabolic potential in different habitats of global hydrothermal vents, including vent plumes, background seawater, diffuse fluids, and sediments.

**Results:** From 34 samples collected at eight vent sites, a total of 4662 viral populations (vOTUs) were recovered from the metagenome assemblies, encompassing diverse phylogenetic groups and defining many novel lineages. Apart from the abundant unclassified viruses, tailed phages are most predominant across the global hydrothermal vents, while single-stranded DNA viruses, including Microviridae and small eukaryotic viruses, also constitute a significant part of the viromes. As revealed by protein-sharing network analysis, hydrothermal vent viruses formed many novel genus-level viral clusters and are highly endemic to specific vent sites and habitat types. Only 11% of the vOTUs can be linked to hosts, which are the key microbial taxa of hydrothermal habitats, such as *Gammaproteobacteria* and *Campylobacterota*. Intriguingly, vent viromes share some common metabolic features in that they encode auxiliary genes that are extensively involved in the metabolism of carbohydrates, amino acids, cofactors, and vitamins. Specifically, in plume viruses, various auxiliary genes related to methane, nitrogen, and sulfur metabolism were observed, indicating their contribution to host energy conservation. Moreover, the prevalence of sulfur-relay pathway genes indicated the significant role of vent viruses in stabilizing the tRNA structure, which promotes host adaptation to steep environmental gradients.

**Conclusions:** The deep-sea hydrothermal systems hold untapped viral diversity with novelty. They may affect both vent prokaryotic and eukaryotic communities and modulate host metabolism related to vent adaptability. More explorations are needed to depict global vent virus diversity and its roles in this unique ecosystem.

**Keywords:** Deep-sea hydrothermal vent, Viruses, Metagenomes, Viral ecology

## Background

Deep-sea hydrothermal vents are one of the most extreme and dynamic environments on Earth [1]. In this dark world, mixing between anoxic hydrothermal fluids and oxic cold seawater results in wide chemical and thermal gradients, providing energy sources for the vent

ecosystems [2]. Unlike most ecosystems that are fueled by photosynthesis, biological productivity is primarily driven by chemoautotrophs in deep-sea hydrothermal vents [3]. Chemoautotrophs use the energy produced by the oxidation of sulfur, hydrogen, methane, ammonia, or iron to fix carbon [2, 4], converting dissolved inorganic carbon into the organic phase within the biota. Diffuse vent fluids are hot spots of primary productivity in the deep ocean and provide a window into the subseafloor microbial habitat [5, 6]. Hydrothermal fluids are highly diluted in plumes, which can rise hundreds of meters and

<footnote>*Correspondence: shaozz@163.com</footnote>

<footnote n="1">¹ Key Laboratory of Marine Genetic Resources, Third Institute of Oceanography, Ministry of Natural Resources, Xiamen, China
Full list of author information is available at the end of the article</footnote>
disperse hundreds of kilometers away from their source and impact broader deep-sea microbial communities and biogeochemistry [7]. Over the past decade, significant efforts have been made to explore the source, diversity, and function of the microbes inhabiting hydrothermal vents [3, 5–12]. These studies have suggested that the prokaryotic communities in hydrothermal plumes are distinct from those in diffuse fluids [13] and hydrothermal sediments [7, 10]. As important components of hydrothermal vent microbiomes, viral communities have received less attention.

Viruses are the most abundant, pervasive, and genetically diverse biological entities in the biosphere [14]. In the ocean, the total estimated number of viruses is approximately $10^{30}$, making up the second largest relative biomass (but the most abundant) in comparison with prokaryotes and protists, despite their small size [15, 16]. They play a pivotal role in marine ecosystems not only by lysing their hosts but also through horizontal gene transfer and manipulating host metabolism via the expression of viral-encoded auxiliary metabolic genes (AMGs) [17]. Each day, viruses in surface waters kill 20–40% of prokaryotes and release up to $10^9$ tons of carbon and other nutrients, which has a significant influence on ocean biogeochemical cycles [16]. Additionally, it is estimated that marine viruses transduce approximately $10^{14}$–$10^{17}$ Gbp of DNA per day [18], affecting host diversity and function. Comparatively, our knowledge of viral diversity and processes in the deep sea is quite limited, partially due to the difficulties in obtaining and processing samples from the deep sea.

In deep-sea hydrothermal vent ecosystems, virus-like particles are more abundant than prokaryotes and are believed to have a profound impact on microbial communities [19]. The viral abundances of hydrothermal plume samples were reported to be $10^5–10^6$ VLPs ml⁻¹, higher than in the surrounding seawater samples [20, 21]. Moreover, it has been suggested that hydrothermal vent microbes harbor substantial populations of temperate viruses [22–24], which may improve host fitness and facilitate horizontal gene transfer. A well-known example of phage AMGs is the gene encoding the reverse dissipative sulfite reductase (rdsr) [25]. The alpha (rdsrA) and gamma (rdsrC) subunit genes of this enzyme were identified in hydrothermal plume phages that putatively infect sulfur-oxidizing bacteria, suggesting that viruses play a direct role in the sulfur cycle. The AMGs of hydrothermal vent viruses are also involved in many other metabolic pathways, including nitrogen, methane metabolism, and amino acid biosynthesis [24, 26], or they may even compensate for novel metabolic pathways for their host microorganisms [19]. These findings suggest that hydrothermal vent viruses are a large reservoir of genetic

diversity and have complex interactions with their hosts and habitats, which remain to be fully elucidated.

To date, microbes identified in hydrothermal vent habitats have largely remained uncultured, and few virus isolates have been reported [27, 28]. Instead, metagenomics technology is commonly applied to characterize the diversity, ecology, and evolution of hydrothermal vent viruses [19, 22–25, 29–31]. Some of these studies have tried to identify viruses in the cellular fractions (0.22 µm filtered), especially in recent years [29–31], because many bioinformatic tools now enable us to recover and analyze viral sequences from complex metagenomes that are generated without viral particle enrichment [32]. The cellular metagenomes may contain sequences of integrated/extrachromosomal proviruses, viruses undergoing the lytic cycle, large virions that are retained on filters, and small virions that are adsorbed onto filters, providing new insights into viral communities and virus-host interactions [29, 33–35].

In this study, we sought to obtain a comprehensive view of viral populations and their interactions with hosts in deep-sea hydrothermal vents. For this purpose, metagenomic data were generated from plume and background seawater samples collected at two hydrothermal fields on the Carlsberg Ridge (CR), northwest Indian Ocean. In addition, 29 publicly available metagenomes from different hydrothermal vent habitats around the world were compiled for analysis. Here, we provide a characterization of the community structure, the virus-host associations, and the potential ecological roles of viruses in deep-sea hydrothermal vents across the global oceans. A collection of hydrothermal vent viral genomes was recovered from microbial metagenomic datasets, and the unique features of the vent viral communities were revealed. To date, this is the largest survey of viruses inhabiting deep-sea hydrothermal vents, and the results will expand our understanding of viral diversity and functions in extreme marine ecosystems.

# Results and discussion

## Overview of deep-sea hydrothermal vent microbial communities

To explore the community composition of microorganisms in the CR hydrothermal vents, metagenomes were sequenced from five plume and background seawater samples collected at the “Wocan” and “Tianxiu” hydrothermal fields (Supplementary Table 1). For comparison, metagenomes of 29 other samples derived from seven different hydrothermal vent sites in the Pacific Ocean, Atlantic Ocean, and southwest Indian Ocean were retrieved from the public database (Supplementary Table 1). These sites are abbreviated as follows: Axial Seamount (Axial), Southwest Indian Ridge (SWIR), Menez
Gwen (Menez), Eastern Lau Spreading Center (Lau), Guaymas Basin (Guaymas), Mid-Cayman Rise (Cayman), and Southern Mariana Trough (Mariana). The final datasets include 34 metagenomes representing four different types of habitats in hydrothermal environments, i.e., vent plumes (n = 16), background seawater (n = 5), diffuse fluids (n = 6), and sediments (n = 7). The sampling locations are shown in Fig. 1.

The 16S rRNA gene fragments (16S miTags) were extracted from clean reads of these metagenomes for taxonomic profiling [36]. Classification of the 16S miTags at the phylum level (the class level for Proteobacteria) revealed dominance by Gammaproteobacteria in almost all of the metagenomes, followed by Deltaproteobacteria, Alphaproteobacteria, and Bacteroidetes (Supplementary Fig. 1A). Campylobacterota (previously named Epsilonproteobacteria [37, 38]) showed high relative abundances in diffuse fluid and some plume metagenomes. The dominant archaeal lineage was Thaumarchaeota, which occurred primarily in seawater and hydrothermal plumes. Among the samples, taxonomic composition at the phylum level was more similar between samples from the same hydrothermal vent site, and samples from the same habitat type tended to be clustered together. It should be noted that the plume samples from CR (this study) and Lau [11] were processed by multiple displacement amplification

(MDA). To assess the bias introduced by the MDA process, we compared the taxonomic composition inferred from CR metagenomic data to those from 16S rRNA gene pyrosequencing (Jiang et al., unpublished). Similar to previous observations in Lau plume samples [11], no significant differences were found in the abundance patterns, suggesting that MDA bias did not obscure the major trends.

De novo assembly and binning of these metagenome sequencing data resulted in the reconstruction of 581 high- or medium-quality prokaryotic metagenome-assembled genomes (MAGs) [39] with ≥ 50% completeness and ≤ 10% contamination. These MAGs were clustered at 95% average nucleotide identity (ANI), representing species-level groups spanning 54 phyla, including 515 bacterial and 66 archaeal MAGs (Supplementary Fig. 1B). Most of the bacterial MAGs belong to dominant lineages, such as Gammaproteobacteria (n = 138), Bacteroidetes (n = 54), Deltaproteobacteria (n = 53), and Campylobacterota (n = 41), while archaeal MAGs are primarily affiliated with Candidatus Thermoplasmatota (n = 25) and Thaumarchaeota (n = 18). These MAGs, in addition to 440 high-quality single-cell amplified genomes (SAGs) retrieved from the CR hydrothermal vents and the publicly available genomes of hydrothermal vent microbial isolates, provide a good basis for investigating the connections between viruses and prokaryotes.

<figure><img src="image_3.png" /></figure>
## Diversity and phylogeny of hydrothermal vent viruses

VirSorter v1.0.6 [40] and VIBRANT v1.2.0 [41] were used to identify viral contigs in the hydrothermal vent metagenome assemblies, followed by manual curation. Contigs ≥ 5 kb or ≥ 2 kb and circular were pooled together, resulting in 8847 putative viral sequences. We retained the small circular contigs because these may represent small single-stranded DNA (ssDNA) virus genomes (such as phages in the family Microviridae and eukaryotic viruses in the families Circoviridae, Geminiviridae, and Smacoviridae). The 8847 candidate viral contigs were then clustered at 95% ANI over 80% of the sequence length, producing 4662 species-level viral populations (viral operational taxonomic units (vOTUs)) [10] (Supplementary Table 2). These vOTUs ranged from 2000 to 226,341 bp in size (total length of the vOTU contigs = 37,751,900 bp and N50 of the vOTU contigs = 10,916 bp). The largest contig has a length of more than 200 kb and might be classified as a “huge phage” [42] rather than a nucleocytoplasmic large DNA virus (NCLDV) [43], according to the presence of phage-specific genes [42].

As determined by CheckV [44], almost half of the vOTUs (47%) represented viral genomes of medium quality and above (Supplementary Fig. 2A). Based on the presence of terminal repeats and provirus integration sites, 1731 vOTUs were identified as complete viral genomes (37%). Of these, 1727 vOTUs were circular and were probably free viruses, while 4 vOTUs were predicted to be intact proviruses.

A large proportion of the hydrothermal vent vOTUs were classified as double-stranded DNA (dsDNA) viruses of the order *Caudovirales* (45%), of which the family *Myoviridae* was predominant (Supplementary Fig. 2B). Because the terminase large subunit (TerL) gene is conserved in all head-tail phages [45], a phylogenetic analysis of the TerL gene was performed to assess the diversity and genetic distance of *Caudovirales* in hydrothermal vents. Within the viral contigs, a total of 638 complete ORFs encoding the TerL genes were identified and used to construct the phylogenetic tree (Fig. 2). Most of the sequences from hydrothermal vent metagenomes fell into 14 known lineages, which were defined by different

<figure><img src="image_4.png" /></figure>
DNA packaging strategies, while the rest of them formed
several novel branches, indicating the remarkable diver-
sity of head-tail phages in our datasets.

ssDNA viruses accounted for approximately 36% of
the total set of vOTUs, with primary assignment to the
family Microviridae. Phages in the family Microviridae
are among the smallest DNA viruses [46], with circular
genomes ranging from 4.2 to 6.5 kb [47, 48]. A total of
872 complete or near-complete genomes of Microviridae
were recovered from hydrothermal vent metagenomes.
The amino acid sequences of the well-conserved major
capsid protein VP1 were used as a phylogenetic marker
for the classification of these viruses. The majority of the
VP1 sequences showed < 70% shared identities relative
to their best matches in NCBI's NR database and showed
< 70% shared identities compared to each other, reflect-
ing high levels of divergence. The phylogenetic analysis
showed that most (694 genomes) of the hydrothermal
vent-derived Microviridae belonged to the subfamily
Gokushovirinae, followed by group D (143 genomes). A
few sequences were clustered within the clades Pequeno-
virus (8 genomes), Pichovirinae (4 genomes), Bullavirinae
(1 genome), and Alpavirinae (1 genome), and the other
22 genomes were not clustered with any known subfami-
lies (Fig. 3A). Interestingly, the new clade seemed to con-
sist of viruses with small genomes (approximately or less
than 4 kb). One of these has a genome size of 3559 bp
and encodes only three putative ORFs, including a capsid

protein, a replication initiator, and a protein of unknown
function. It represents the smallest microvirus recovered
in our datasets and the smallest ssDNA phage with the
least ORFs reported to date.

CRESS-DNA (circular rep. encoding single-stranded
DNA) viruses, including Circoviridae and its related fam-
ilies, were also highly represented (15% of the vOTUs,
Supplementary Fig. 2B, Supplementary Table 2). This
group of eukaryotic viruses has small circular genomes
and commonly encodes only 2 proteins, of which the
replication initiation protein (Rep) is the only univer-
sally conserved gene [49]. Based on the sequences of
Rep proteins, 694 CRESS-DNA viruses identified from
hydrothermal vent metagenomes were clustered within
13 known families, while the other 145 sequences defined
several potentially new clades (Fig. 3B).

The NCLDVs are another group of eukaryotic viruses [43], including the families *Poxviridae*, *Iridoviridae*, *Ascoviridae*, *Asfarviridae*, *Marseilleviridae*, *Mimiviridae*, and *Phycodnaviridae*, as well as several lineages of unclassified viruses. Viruses in this group were also present in our hydrothermal vent metagenomics datasets (Supplementary Fig. 2B, Supplementary Table 2). As shown in Supplementary Table 2, 27 vOTUs were classified as viruses from the *Mimiviridae*, *Phycodnaviridae*, or other NCLDV families. However, the sizes of these vOTUs ranged from 6112 to 24,611 bp and only represented small genome fragments of the viruses.

<figure><img src="image_5.png" /></figure>
Approximately 19% of the total vOTUs did not show any significant sequence similarity to any known viral families and could not be taxonomically classified for the time (Supplementary Fig. 2B).

It should be noted that all the hydrothermal plume, diffuse fluid, and background seawater samples used in this study were passed through 0.22 µm filters (Supplementary Table 1 and references therein). Giant viruses, integrated proviruses, and actively infecting viruses within the cells would be retained on the membrane, while most free virus particles of small size would be lost during this step, such as the non-tailed ssDNA viruses. Thus, we analyzed the metagenomes of the cellular fraction and the virus-like particle (VLP) fraction of the sediment samples from the southwest India Ridge [19] to evaluate the extent to which they reflect viral diversity. The results showed that the cellular fraction was comparable to the VLP fraction with respect to viral recovery (data not shown). The number of recovered CRESS-DNA viruses in VLP fractions doubled those recovered in cellular fractions, but the number of microviruses recovered in cellular fractions was greater than those recovered in VLP fractions, despite their small size. One explanation for these differences is that the intracellular microviruses were captured on filters, while most of the CRESS-DNA viruses replicating in eukaryotic hosts were excluded

from sampling. No significant difference was observed for other viral groups. For hydrothermal plume and diffuse fluid samples involved in this study, metagenomic data on the viral fractions were not available. However, a comparative analysis of the cellular and viral metagenomes derived from a fluid sample was performed in a previous study [24], and high enrichment of mobile elements and proviruses was observed in the cellular fraction. Given that a high proportion of viruses in hydrothermal ecosystems are lysogenic [22–24], we supposed that the viral sequences identified in the cellular metagenomes could represent the diversity of hydrothermal viruses to a large extent. However, to fully characterize the viromes in deep-sea hydrothermal vents, the metagenomes of both the VLP fraction and cellular fraction are still needed, preferably with the addition of RNA-seq data, which will enable us to discover and analyze RNA viruses.

## Viral communities across different zones of hydrothermal vents

To investigate the viral community structure in hydrothermal ecosystems, the relative abundances of vOTUs in each metagenome were calculated and normalized (Fig. 4, Supplementary Fig. 3). The 34 samples used in this study were collected from 8 different hydrothermal vents across various geographical zones, including

<figure><img src="image_6.png" /><figcaption>**Fig. 4** Bubble plot showing the relative abundance of viral groups in different hydrothermal samples. The metagenomic samples from this study are indicated in bold</figcaption></figure>
hydrothermal plume, background water, diffuse fluid, and sediment samples. As a result, the vOTU abundance patterns were primarily clustered by sample types and secondarily by hydrothermal vent sites (Supplementary Fig. 3). The viral community composition of hydrothermal vent sediments was significantly different from that of other habitats. Hydrothermal plumes and the surrounding deep-sea water samples showed similar vOTU abundance patterns, as plume and water samples from the same hydrothermal vent field were always clustered together. This result is consistent with previous studies suggesting that plume microbial communities resemble those from background seawater samples [7], indicating that hydrothermal plumes are strongly influenced by ambient seawater. Overall, these results showed that the virome structures varied across different hydrothermal vent habitats and different hydrothermal vent fields. The hierarchical clustering of samples based on vOTU and 16S miTags showed similar patterns (Supplementary Fig. 1A, Supplementary Fig. 3), implying a close link between viral and prokaryotic communities.

Similar to other studies of marine viromes, the majority of the hydrothermal vent viral communities were composed of unclassified viruses (up to 69% of the total viral reads, Fig. 4). At the family level, the tailed phage *Myoviridae* (on average 9.2%, 7.8%, and 1.5% of plume, fluid, and sediment samples, respectively) was the most dominant group of dsDNA viruses in most samples, followed by Siphoviridae (9.1%, 4.6%, and 2.4%) and Podoviridae (7.6%, 4.8%, and 0.7%).

The ssDNA viruses also accounted for a large fraction of viral communities. In all of the sediment metagenomes, the relative abundance of Microviridae was higher than that of any other dsDNA virus family. The CRESS-DNA virus group was present in most plume samples and three of the sediment samples, accounting for approximately 2.2% of the total viral reads (Fig. 4). It is possible that the enrichment of ssDNA viruses in these datasets was caused by the MDA process, which used phi29 DNA polymerase and preferentially amplified small circular ssDNA molecules [50, 51]. However, the quantification of viral DNA without amplification also revealed the dominance of ssDNA viruses in the total DNA viral assemblages of deep-sea sediments [52]. Thus, we suggest that ssDNA viruses are abundant and play an important role in hydrothermal vent environments.

The eukaryotic NCLDVs accounted for 0.4% of the total viral reads in hydrothermal plume samples, on average. Of these families, Mimiviridae and Phycodnaviridae are the most abundant. In contrast, samples from diffuse fluid and sediment contained fewer NCLDVs. Large DNA viruses in this group have been frequently detected in marine metagenomes, including those derived from

hydrothermal vents [30, 53–55]. The NCLDV sequences in metagenomic datasets may come from marine unicellular eukaryotes or free giant viruses, but their roles in these ecosystems are largely unknown.

## Hydrothermal vent viruses are novel and endemic

To gain further insight into this viral diversity and distribution, we then used an extensively validated, network-based method [56] to investigate the relationship among hydrothermal vent vOTUs and viral sequences identified from other marine ecosystems. The 4662 vOTUs recovered from hydrothermal vents were compared to NCBI Viral RefSeq v97 and viral contigs in other marine metagenomic datasets: GOV 2.0 seawater [57] and cold seeps [58]. vConTACT2 was used to de novo predict genus-level groups (viral clusters, VCs) from viral population data [56]. Lastly, a total of 16,618 clusters were generated, reflecting the huge and unexplored diversity of marine viruses (Fig. 5A). As the largest marine virus database to date, the GOV 2.0 datasets contributed the largest number of VCs (14,968 clusters), while taxonomically known viruses from NCBI RefSeq only formed 482 clusters. Consistent with previous studies [58], the viral compositions in different habitats varied considerably, with only 30 clusters shared by all three marine metagenomic datasets (Fig. 5A).

The 4662 hydrothermal vent vOTUs were grouped into 1138 genus-level clusters, of which only 32 VCs contained genomes of known viruses from NCBI RefSeq. Moreover, 584 VCs (~51%) were exclusively composed of hydrothermal vent viruses, which may represent completely new genus candidates (Fig. 5A). These VCs contained 1107 vOTUs, most of which were ssDNA viruses within Microviridae (548 vOTUs) and the CRESS-DNA virus families (122). Approximately one-third of the vOTUs (314) could be classified as tailed phages, probably of novel genera within *Myoviridae* (144), Podoviridae (65), Siphoviridae (57), and unclassified *Caudovirales* (48). The remaining 123 vOTUs could not be taxonomically classified at the family or even higher levels.

Within hydrothermal vent habitats, a high proportion of the viruses seemed to be endemic, given that a majority (818 VCs, 71.9%) of the VCs only occurred in a specific vent site (Fig. 5B). This observation was consistent with previous work showing that most viruses in hydrothermal vent fluids have limited distributions [31]. Generally, the VC richness in hydrothermal plumes was higher than those in sediments or diffuse fluids, suggesting a greater viral diversity in hydrothermal plumes. Although the viral communities of different hydrothermal plume samples were somewhat similar at the family level, they shared a small fraction of clusters. The number
<figure><img src="image_7.png" /></figure>

of clusters shared between sediments and fluid was even lower, and only one cluster was detected across all sample types (Fig. 5B). The cluster turned out to be viruses in the family Microviridae, reinforcing the ubiquity of this group.

## Virus-host connections in hydrothermal vents

The interactions between viruses and their hosts exert a strong influence on microbial diversity and are essential for understanding the ecology and functioning of microbial communities [16]. We sought to link the hydrothermal vent vOTUs to their potential hosts by using a combination of four in silico methods based on their CRISPR spacer match, tRNA match, sequence similarity, and k-mer frequencies [35]. As a result, putative targeted hosts were predicted for a small fraction (494 vOTUs, ~11%) of the hydrothermal vent vOTUs (Fig. 6, Supplementary Table 3). Specifically, 230 vOTUs were linked to MAGs (Supplementary Fig. 1B) and SAGs derived from hydrothermal vents. Most connections were predicted by WIsH and CRISPR spacer matching (252 viral-host pairs each), 114 pairs by sequence homology, and 119 pairs by tRNA matching. Among them, the linkages of 78 viral-host pairs were supported by two or more prediction strategies. The majority (~89%) of these vOTUs were predicted to infect a specific host, and only 31 OTUs were linked to hosts from different prokaryotic phyla. This result is consistent with the common perceptions and previous findings that most viruses have a narrow host range [35, 40, 58].

The predicted hosts of hydrothermal vent phages include bacterial and archaeal species from 39 different phyla (Fig. 6, Supplementary Table 3). Seventy-nine vOTUs were linked to archaea, of which the phylum *Candidatus* Thermoplasmatota was the most frequently predicted (33 associated vOTUs). This newly proposed phylum contains the Marine group II (MGII) and Marine group III (MGIII) archaea [59]. MGII dominate ocean surface waters and may play important roles in the marine carbon cycle [60], whereas members of MGIII live in deep mesopelagic and bathypelagic environments at relatively low abundance [61]. Both of these groups have been found in deep-sea hydrothermal vents and are thought to contribute to organic compound degradation [62].

*Gammaproteobacteria* were the most frequently predicted bacterial hosts, with 144 associated vOTUs, followed by *Actinobacteria* (39 vOTUs), *Alphaproteobacteria* (37 vOTUs), *Bacteroidota* (34 vOTUs), *Firmicutes* (27 vOTUs), and *Campylobacterota* (26 vOTUs). These groups were among the most abundant and active bacterial lineages in the hydrothermal vent ecosystems, as previously reported [9, 12] and as revealed in this study (Supplementary Fig. 1A). For example, *Gammaproteobacteria* is a large bacterial class with metabolic versatility and is observed in almost all habitats surrounding hydrothermal vents [63]. We found that the most frequently predicted hosts within *Gammaproteobacteria* were the genera *Acinetobacter* (with 15 associated vOTUs), *Alteromonas* (15 vOTUs), *Pseudomonas* (14 vOTUs), and *Alcanivorax* (8 vOTUs), which were
<figure><img src="image_8.png" /></figure>

dominant in most of the samples involved in this study. According to the well-known kill-the-winner hypothesis [16], abundant microbes are more likely to be infected and lysed by viruses because a high population density will increase the host-virus encounter rate [64, 65]. Thus, it is not surprising that many viruses target *Gammaproteobacteria* in hydrothermal vents. These viruses showed high abundances (Supplementary Fig. 4) and might play important roles in regulating the vent microbial communities.

As another abundant and ubiquitous group that inhabits hydrothermal vent environments, *Campylobacterota* is a group of chemolithotrophic primary producers that primarily use sulfur compounds and hydrogen as electron donors [66] and are regarded as indicators of hydrothermal activity and passive tracers of vent fluids [9, 63, 67]. However, only a few potential prophage regions have been reported in the complete genomes of deep-sea *Campylobacterota* [68, 69], and only one of them has been isolated [70] to date. In this study, 26 vOTUs were shown to be potentially able to infect members of the phylum *Campylobacterota*, particularly the genus *Sulfurimonas* (14 associated vOTUs). These findings included viruses from the families Myoviridae, Podoviridae, Siphoviridae, Herelleviridae, and Microviridae,

while 8 vOTUs remained unclassified at the family level, indicating the unrevealed diversity of viruses infecting hydrothermal vent Campylobacterota. Further analysis of these viral genomes will provide new insight into the interactions of this ecologically important group and their phages.

CRESS-DNA viruses from the existing families were reported to infect hosts across the eukaryotic domain, including plants, fungi, and animals [49], but a recent study based on CRISPR analysis suggested that viruses from the CRESS-DNA family Smacoviridae infected methanogenic archaea instead of humans [71]. Interestingly, our results also revealed some potential connections between CRESS-DNA viruses and bacterial or archaeal hosts (Fig. 6, Supplementary Table 3), suggesting a broader host range for this group. To date, CRESS-DNA virus isolates with definitive hosts include members from five CRESS-DNA families [71]. However, the number of CRESS-DNA viruses discovered in metagenomics surveys now far exceeds the number of biologically characterized viral isolates [49]. Considering the diversity of CRESS-DNA viruses [72, 73] and their origin from bacterial rolling circle-replicating plasmids [74, 75], it is possible that some CRESS-DNA viruses infect hosts beyond eukaryotes.
<figure><img src="image_9.png" /></figure>

## AMGs of vent viruses are involved in various metabolic pathways

Viral infections can affect host metabolism via the expression of viral AMGs. To better understand the ecological impact of viruses in deep-sea hydrothermal ecosystems, we searched the AMGs in hydrothermal vent viral genomes and calculated their relative abundances. Based on the comprehensive annotation of viral ORFs, a total of 608 genes were considered to be putative AMGs (Fig. 7, Supplementary Table 4). Sequence homology searches against the NCBI NR database showed that a large portion of these AMGs were probably acquired from *Proteobacteria*, especially the *Gammaproteobacteria* and *Alphaproteobacteria*, while ~21% of the AMGs came from unclassified source species (Fig. 7B). The origins of AMGs reflect virus-host connections, because phages generally acquire AMGs from their hosts [76]. According to the KEGG annotation, the identified AMGs of hydrothermal vent viruses were involved in a variety of metabolic pathways, including those related to carbohydrate metabolism, amino acid metabolism, and the metabolism of cofactors and vitamins (Fig. 7A). This trend is consistent with the viral metabolic profiles revealed by analyzing 6 hydrothermal vent metagenomes [41], suggesting that there are some common features in the patterns of metabolic capabilities of hydrothermal vent viromes.

The AMG composition and abundance profiles across the hydrothermal vent samples indicated that viruses in hydrothermal plumes encoded a larger number of AMGs with diverse functions (Fig. 7C). Most of the AMGs had higher abundances in the plume samples compared with hydrothermal fluid or sediments. This observation may reflect the role of viruses in facilitating host adaptation to the dynamic nature of hydrothermal plumes and is in congruence with the metabolic versatility of microbes therein [10]. Specifically, AMGs involved in energy metabolism were found only in plume samples, including those related to methane, nitrogen, and sulfur metabolism (Fig. 7C). Several AMGs were associated with sulfur metabolism pathways, such as genes coding for phosphoadenosine phosphosulfate reductase (*cysH*), adenylylsulfate kinase (*cysC*), sulfate adenylyltransferase (*sat*), and dissimilatory sulfite reductase subunits A (*dsrA* or *rdsrA*, as referred to in a previous study [25]). The *cysH* and *cysC* genes are involved in assimilatory sulfate reduction, whereas the *sat* and *dsrA* genes are related to dissimilatory sulfur reduction/oxidation [77]. Sulfur oxidation and sulfate reduction are both important parts of sulfur cycling in hydrothermal vent ecosystems [78], and the presence of these AMGs suggested that phages participate extensively in these pathways.
The most abundant AMG identified in hydrothermal vent metagenomes was DNA cytosine methyltransferase (DNMT1, *dcm*). However, viral *dcm* and another 13 AMGs were present in metagenomes derived from several diverse environments and were thus considered to perform central functions in the viral life cycle [38]. Except for these globally conserved AMGs, the thiouridine synthase subunit E (*tusE*, a homolog of *dsrC*) gene related to the sulfur relay system was identified in high abundance. This gene encodes a sulfur transfer protein for tRNA thiol modifications, which is required for protein synthesis machinery [79]. The expression of tRNA thiolation genes has been reported to increase the stability of tRNA structure and is essential for bacterial survival at high temperatures [80, 81]. In addition, it has been suggested that the sulfur relay system is involved in microbial tolerance against acid stress [82, 83], heavy metals [84], and organic solvents [85]. Therefore, this AMG may benefit the hosts by improving their adaptability to various stress conditions and thus provide them with a growth advantage in hydrothermal vent environments in which the temperature and chemistry fluctuate substantially.

It is worth noting that some of the predicted AMGs may not be bona fide AMGs. First, it is difficult to avoid false-positive predictions of viral contigs completely since these metagenomes were derived from cellular fractions. Although the candidate AMGs were co-localized with at least one viral hallmark gene, they might belong to the host regions that were retained due to the miscall of a prophage boundary [32]. Another concern is the true function of candidate AMGs. As described above, some metabolic genes are more likely to be involved in the viral life cycle rather than the host metabolism, such as methyltransferases, glycosyl transferases, glycoside hydrolases, and adenylyltransferases [86–88]. Thus, further investigations, such as genome context assessments and functional analyses of putative AMGs, are required for a better understanding of the viral impacts on hydrothermal vents.

## Conclusions

In this study, we explored the viral community of CR hydrothermal vents and additional hydrothermal vent sites across a wide geographical area based on metagenomic data. We found that deep-sea hydrothermal systems are large reservoirs of novel viruses. Both vent prokaryotic and eukaryotic communities are affected by viruses that target diverse hosts and modulate their metabolisms. These interactions shape the structure and function of microbiomes in deep-sea hydrothermal vents and may profoundly influence the broader oceans via vent fluid circulation and plume drift. Our exploration highlights

global vent viral diversity and suggests the significant roles of viruses in this unique ecosystem. With improvements in deep-sea sampling and culture-dependent and culture-independent technologies, a more comprehensive understanding of virus-host interactions in hydrothermal ecosystems is soon to be expected.

## Methods

### Sample collection

Plume and background seawater samples were collected from two different deep-sea hydrothermal vents, “Wocan” and “Tianxiu,” at the Carlsberg Ridge of the northwest Indian Ocean during the COMRA cruise DY 38 in March 2017 (Supplementary Table 1). The human-operated vehicle “Jiaolong” was used to collect 1.5-L water samples in individual Niskin bottles. The collected water was filtered through 0.22 µm polycarbonate membranes (diameter 45 mm; Whatman, Clifton, NJ, USA) and frozen at –80 °C on board for DNA extraction. For the single-cell sequencing, plume samples from the “Wocan” vent were fixed with glycerol-Tris-EDTA buffer [89] and frozen at –80 °C until further processing.

### DNA extraction and metagenomic sequencing

The total DNA was extracted from filtration membranes as described previously [7]. Multiple displacement amplification of genomic DNA was performed using the illustra Ready-To-Go GenomiPhi V3 DNA Amplification Kit (GE Healthcare, Piscataway, NJ, USA). Paired-end library was constructed using NEXTFLEX Rapid DNA-Seq (Bioo Scientific, Austin, TX, USA). Adapters containing the full complement of sequencing primer hybridization sites were ligated to the blunt end of fragments. Shotgun sequencing was performed on Illumina HiSeq PE150 platform (Illumina Inc., San Diego, CA, USA) at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China) according to the manufacturer’s instructions (www.illumina.com).

### Metagenome assembly and annotation

The raw reads obtained by Illumina paired-end sequencing were trimmed and quality filtered using fastp software [90]. Clean reads were then assembled using MEGAHIT [91] with default options. Metagene [92] was used to predict protein coding sequences (CDS) from the metagenomic assemblies. Nonredundant genes generated by CD-HIT clustering (with 95% shared sequence identity and 90% coverage) were aligned to the NCBI NR database using BLASTp (with an e-value cut-off of 1e-5) for taxonomic classification. Functional annotations were conducted based on comparisons with KEGG [93], eggNOG v5.0 [94], and the Carbohydrate-Active enZymes (CAZy) databases [95] using the BLASTp [96] program
(with an *e*-value cut-off of 1e-5). Additionally, 29 publicly available metagenomic datasets generated from hydrothermal vent samples (Supplementary Table 1) were downloaded from the NCBI Sequence Read Archive (SRA) database and were quality controlled, assembled, and annotated as described above.

## Metagenomic binning and metagenome-assembled genome (MAG) classification

Contigs larger than 1500 bp from the final assemblies were included for metagenomic binning using MetaBAT2 v2.12.1 [97] with default parameters. The original bins were then run through MetaWRAP's reassemble_ bins module [98] to improve their quality, and the completeness and contamination of the resulting bins were evaluated by CheckM [99]. The high- and medium-quality bins (completeness ≥ 50% and contamination ≤ 10%) were then dereplicated at 95% average nucleotide identity (ANI) using dRep v2.3.2 [100], resulting in 581 species-level MAGs. The taxonomic assignment of the MAGs was performed using the GTDB-Tk package v0.3.2 [101] and was ultimately converted into the corresponding NCBI taxonomy. The phylogenomic relationships of MAGs were inferred using IQ-TREE 2 [102] based on a concatenation of 120 bacterial or 122 archaeal marker genes identified by GTDB-Tk. Support for nodes in the ML trees was evaluated with 1000 ultrafast bootstrap replicates [103], and the generated tree was visualized using iTOL v4 [104].

## SAG library construction, sequencing, and analysis

The sequencing of plume-derived SAGs was performed in the Bigelow Laboratory Single Cell Genomics Center (https://scgc.bigelow.org/). Fluorescence-activated cell sorting, cell lysis, multiple displacement amplification, Illumina sequencing, and de novo genome assembly were conducted as previously described [89]. The quality assessment and taxonomic classification of SAGs were performed as described above for MAGs.

## Identification of viral contigs

Contigs ≥ 2 kb from metagenome assemblies were used to recover viral sequences. VirSorter analysis [40] was run with the parameter "--db 2 (viromes database)," and only the highest confidence contig categories 1, 2, 4, and 5 were included in this study, with categories 4 and 5 being manually curated. Contigs containing at least one of the viral hallmark genes (such as "virion structure," "capsid," "portal," "head," "tail," "baseplate," or "terminase") were retained. VIBRANT v1.2.1 [41] was also used to identify viral contigs using the default parameters, and only the complete circular, high- and medium-quality drafts were kept for further analysis. The contigs identified by

VirSorter and VIBRANT were then compiled and clustered at 95% shared nucleotide identity and 80% coverage [105], yielding 4662 viral populations, or viral operational taxonomic unit (vOTUs). Lastly, the CheckV pipeline was used to estimate the completeness of the viral genomes and to predict viral lifestyles [44].

## Abundance profiling in metagenomics data

For taxonomic profiling of prokaryotic communities, 16S miTags were recovered from the metagenomic reads using phyloFlash [36]. The extracted 16S miTAGS were mapped against the SILVA SSU Ref. database (v132) [106] for taxonomic assignment. To calculate the relative abundances of vOTUs and host microorganisms in each sample, clean reads from metagenomes were mapped to the viral contigs or microbial genomes using the CoverM package (https://github.com/wwood/CoverM) with contig mode and genome mode, respectively. RPKM (reads per kilobase per million mapped reads) values were selected to represent the relative abundances of the viral and host populations. Pearson correlation was used to calculate the distances between samples for hierarchical clustering. Heatmaps were generated using the pheatmap R package and TBtools [107].

## Host prediction

Four computational host prediction strategies were used to identify virus-host interactions [35]. (i) CRISPR spacers match: a clustered regularly interspaced short palindromic repeats (CRISPRs) spacer database was created for a set of microbial genomes using the MinCED tool [108]. For metagenomics reads, Crass v1.0.1 [109] was used with the default parameters to recover CRISPR spacers and repeat elements. The identified spacers were queried for exact sequence matches against all viral contigs using the BLASTn-short mode in the BLAST+ package [110]. Match requirements were at least 95% identity over 95% spacer length, and only ≤ 1 mismatch was allowed. The corresponding CRISPR direct repeat types were connected to microbial genomes via BLASTn (with an *e*-value cut-off of 1e-10, 100% nucleotide identity) [111]. (ii) Transfer RNA (tRNA) match: tRNAs were recovered from the microbial genomes and viral contigs using ARAGORN with the "-t" option [112]. The identified tRNA sequences were compared using BLASTn, and only a perfect match (100% coverage and 100% identity) was considered indicative of putative host-virus pairs. (iii) Nucleotide sequence homology search [113]: To link prophages with hosts, viral contigs were searched against microbial genomes using BLASTn with the following thresholds: 75% minimum coverage of the viral
contig length, 70% minimum nucleotide identity, 50 minimum bit score, and 0.001 maximum e-value. (iv) k-mer frequencies: WISH v1.0 [114] was run with the default parameters against the host database. Connections were inferred when *p* < 0.001. If multiple hosts were predicted for a vOTU, the one supported by different approaches was chosen as the one with the most confidence. The host database employed for these prediction methods was composed of (i) all reference genomes from the Genome Taxonomy Database (GTDB), (ii) all MAGs (≥ 50% completeness, and ≤ 10% contamination) recovered from hydrothermal vent metagenomes (*n* = 581), (iii) all SAGs obtained from the CR hydrothermal vent (*n* = 440), and (iv) a custom collection of marine microbial genomes from the Marine Culture Collection of China (*n* = 1452).

## Viral taxonomic assignment and network analysis

The predicted open reading frames (ORFs) of the viral contigs were mapped against the NR protein database using DIAMOND v0.9.21 [115], and their taxonomic affiliations were determined using the CAT v5.0.3 package [116] based on the last common ancestor (LCA) algorithm. Contig classification is based on a voting approach of all classified ORFs by summing up all the bit scores from ORFs supporting a specific classification. Protein-sharing network analysis of the hydrothermal vent vOTUs, the reference phage genomes (from NCBI Viral RefSeq version 97), the Global Oceans Viromes 2 (GOV 2.0) datasets [57], and viral contigs from cold seeps [58] was performed using vConTACT v2.0 [56]. Briefly, Prodigal v2.6.3 [117] was used for ORF prediction from the vOTUs. The predicted protein sequences were then subjected to all-to-all BLASTp using DIAMOND, and the BLAST result file was used as input for vConTACT2. The similarity score between vOTUs was calculated based on the number of shared protein clusters, and related vOTUs with a similarity score of ≥ 1 were grouped into viral clusters.

## Construction of phylogenetic trees

For phylogenetic trees, the deduced amino acid sequences of selected marker genes were aligned using the MUSCLE program [118], and the multiple alignments were trimmed with TrimAl v1.2 [119]. IQ-TREE2 [102] was used to infer the maximum-likelihood (ML) tree with the best substitution model selected by Modelfinder [120], and support for nodes in the ML trees was evaluated with 1000 ultrafast bootstrap replicates [103]. The resulting trees were visualized and annotated using FigTree v1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/) and the iTOL v4 online tool [104].

## Identification of auxiliary metabolic genes

Functional annotations on the ORFs in the viral contigs were conducted based on comparisons with the eggNOG v5.0 [94] database using eggNOG-mapper v2 [121]. Genes with a KEGG annotation falling under the “metabolic pathways” category or “sulfur relay system” were considered to be putative vAMGs, as previously defined in VIBRANT [41]. Lastly, manual curation was performed to remove the metabolic genes involved in common viral functions. To generate the abundance profiles for AMGs, clean reads were mapped to the metagenomes using bowtie2 [122], and the RPKM values of each gene were calculated. The sum of the RPKM values of genes with the same KO annotations was used to represent the relative abundance of each gene category.

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s40168-022-01441-6.

**Additional file 1: Supplementary Table 1.** Metagenome datasets used for analysis in this study.

**Additional file 2: Supplementary Figure 1.** Microbial communities in deep-sea hydrothermal vents. (A) Relative abundance of 16S miTag in 34 hydrothermal vent samples. The top 30 most abundant phyla (class level for Proteobacteria) among the metagenomes are shown. (B) Phylogenetic tree of high-quality metagenome-assembled genomes (MAGs) recovered from 34 hydrothermal vent metagenomes. Maximum-likelihood phylogenetic trees of bacterial and archaeal MAGs at the phylum level (class level for Proteobacteria) were inferred from 120 bacterial or 122 archaeal single-copy marker genes, respectively. Support for nodes in the ML trees was evaluated with 1000 ultrafast bootstrap replicates, and bootstrap scores >70% are flagged with dots. The number of MAGs related to viruses and the total number of recovered MAGs in the clade are shown in brackets.

**Additional file 3: Supplementary Table 2.** Characteristics of vOTUs recovered from hydrothermal vent metagenomes. The taxonomic affiliation of vOTUs was inferred using the last common ancestor algorithm in CAT v5.0.3. The completeness of the viral genomes and viral lifestyles were determined using the CheckV pipeline.

**Additional file 4: Supplementary Figure 2.** Genome quality and taxonomic composition of hydrothermal vent vOTUs. (A) Proportion of genome quality categories assessed by CheckV. (B) Taxonomic classification of vOTUs at the family level.

**Additional file 5: Supplementary Figure 3.** Distribution patterns of all hydrothermal vent vOTUs. The relative abundances of vOTUs (y-axis) in each sample (x-axis) were calculated as reads per kilobase per million mapped reads (RPKM values) and were normalized on the log2 scale. The vOTUs and the samples were hierarchically clustered.

**Additional file 6: Supplementary Table 3.** Predicted virus-host connections for vOTUs. The highest confidence host for each vOTU is highlighted.

**Additional file 7: Supplementary Figure 4.** Distribution patterns of viruses and their predicted hosts in deep-sea hydrothermal vents. The relative abundances of vOTUs (top left triangle) and their predicted hosts (bottom right triangle) were grouped by host taxonomy and were normalized on the log2 scale.

**Additional file 8: Supplementary Table 4.** Annotation of putative auxiliary metabolic genes.
## Acknowledgements

The authors are grateful to the whole R/V Xiang-yang-hong 9 team of the scientific cruise DY 381. The authors also thank Prof. Rui Zhang and Min Jin for their helpful suggestions and comments on improving the manuscript.

## Authors' contributions

Conceptualization, RC and ZS; methodology, RC and XL; investigation, LJ; data curation, XL, LJ, and LG; writing—original draft preparation, RC; writing—review and editing, CG and ZS; supervision, ZS; and funding acquisition, ZS and RC. The authors read and approved the final manuscript.

## Funding

This work was funded by the National Key Research and Development Program of China (No. 2018YFC0310705; No. 2018YFC0310701); the grant of Laboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology (Qingdao) (No. OF2019NO05); Scientific Research Foundation of Third Institute of Oceanography, MNR (No. 2018022); Natural Science Foundation of China (No. 42006088); and the China Ocean Mineral Resources R&D Association (COMRA) program (No. DY135-B2-01).

## Availability of data and materials

The metagenomic datasets used for analysis in this study are publicly available in the NCBI SRA repository at https://www.ncbi.nlm.nih.gov/sra, and accession numbers are listed in Supplementary Table 1. Metagenomic data from CR hydrothermal vent were also deposited in the SRA database (accession number PRJNA685608). The sequences of the vOTUs generated from the current study have been deposited in the National Omics Data Encyclopedia (NODE) database at https://www.biosino.org/, accession number OEP003015.

## Declarations

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Competing interests

The authors declare that they have no competing interests.

### Author details

¹Key Laboratory of Marine Genetic Resources, Third Institute of Oceanography, Ministry of Natural Resources, Xiamen, China. ²State Key Laboratory Breeding Base of Marine Genetic Resource, Fujian Key Laboratory of Marine Genetic Resources, Xiamen 361005, China. ³State Key Laboratory for Managing Biotic and Chemical Threats to the Quality and Safety of Agro-Products, Key Laboratory of Biotechnology in Plant Protection of Ministry of Agriculture and Zhejiang Province, Institute of Plant Virology, Ningbo University, Ningbo 315211, China. ⁴Univ Brest, CNRS, IFREMER, IRP 1211 MicrobSea, Laboratoire de Microbiologie des Environnements Extrêmes LM2E, IUEM, Rue Dumont d'Urville, F-29280 Plouzané, France. ⁵Sino-French Laboratory of Deep-Sea Microbiology (MICROBSEA-LIA), Plouzané, France. ⁶Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519000, China.

Received: 16 February 2022 Accepted: 30 November 2022

Published online: 24 December 2022

## References

1. Thornburg CC, Zabriskie TM, McPhail KL. Deep-sea hydrothermal vents: potential hot spots for natural products discovery? *J. Nat. Prod.* 2010;73(3):489–99.

2. Dick GJ. The microbiomes of deep-sea hydrothermal vents: distributed globally, shaped locally. *Nat. Rev. Microbiol.* 2019;17(5):271–83.

3. Lesniewski RA, Jain S, Anantharaman K, Schloss PD, Dick GJ. The metatranscriptome of a deep-sea hydrothermal plume is dominated by water column methanotrophs and lithotrophs. *ISME J.* 2012;6(12):2257–68.

4. Anantharaman K, Breier JA, Sheik CS, Dick GJ. Evidence for hydrogen oxidation and metabolic plasticity in widespread deep-sea sulfur-oxidizing bacteria. *Proc. Natl. Acad. Sci.* 2013;110(1):330.

5. Galambos D, Anderson RE, Reveillaud J, Huber JA. Genome-resolved metagenomics and metatranscriptomics reveal niche differentiation in functionally redundant microbial communities at deep-sea hydrothermal vents. *Environ. Microbiol.* 2019;21(11):4395–410.

6. Trembath-Reichert E, Butterfield DA, Huber JA. Active subseafloor microbial communities from Mariana back-arc venting fluids share metabolic strategies across different thermal niches and taxa. *ISME J.* 2019;13(9):2264–79.

7. Dick GJ, Tebo BM. Microbial diversity and biogeochemistry of the Guaymas Basin deep-sea hydrothermal plume. *Environ. Microbiol.* 2010;12(5):1334–47.

8. Sheik CS, Anantharaman K, Breier JA, Sylvan JB, Edwards KJ, Dick GJ. Spatially resolved sampling reveals dynamic microbial communities in rising hydrothermal plumes across a back-arc basin. *ISME J.* 2015;9(6):1434–45.

9. Djurhuus A, Mikalsen S-O, Giebel H-A, Rogers AD. Cutting through the smoke: the diversity of microorganisms in deep-sea hydrothermal plumes. *R. Soc. Open Sci.* 2017;4(4):160829.

10. Dick G, Anantharaman K, Baker B, Li M, Reed D, Sheik C. The microbiology of deep-sea hydrothermal vent plumes: ecological and biogeochemical linkages to seafloor and water column habitats. *Front. Microbiol.* 2013;4(124).

11. Anantharaman K, Breier JA, Dick GJ. Metagenomic resolution of microbial functions in deep-sea hydrothermal plumes across the Eastern Lau Spreading Center. *ISME J.* 2016;10(1):225–39.

12. Li M, Jain S, Dick GJ. Genomic and transcriptomic resolution of organic matter utilization among deep-sea bacteria in Guaymas Basin hydrothermal plumes. *Front. Microbiol.* 2016;7(1125).

13. Anderson RE, Beltrán MT, Hallam SJ, Baross JA. Microbial community structure across fluid gradients in the Juan de Fuca Ridge hydrothermal system. *FEMS Microbiol. Ecol.* 2013;83(2):324–39.

14. Edwards RA, Rohwer F. Viral metagenomics. *Nat. Rev. Microbiol.* 2005;3(6):504–10.

15. Suttle CA. Viruses in the sea. *Nature*. 2005;437(7057):356–61.

16. Suttle CA. Marine viruses — major players in the global ecosystem. *Nat. Rev. Microbiol.* 2007;5(10):801–12.

17. Breitbart M. Marine viruses: truth or dare. *Annu. Rev. Mar. Sci.* 2012;4(1):425–48.

18. Rohwer F, Prangishvili D, Lindell D. Roles of viruses in the environment. *Environ. Microbiol.* 2009;11(11):2771–4.

19. He T, Li H, Zhang X, Bailey MJ. Deep-sea hydrothermal vent viruses compensate for microbial metabolism in virus-host interactions. *mBio.* 2017;8(4):e00893–17.

20. Ortmann AC, Suttle CA. High abundances of viruses in a deep-sea hydrothermal vent system indicates viral mediated microbial mortality. *Deep-Sea Res. I Oceanogr. Res. Pap.* 2005;52(8):1515–27.

21. Ray J, Dondrup M, Modha S, Steen IH, Sandaa R-A, Clokie M. Finding a needle in the virus metagenome haystack—micro-metagenome analysis captures a snapshot of the diversity of a bacteriophage armoire. *PLoS One.* 2012;7(4):e34238.

22. Williamson SJ, Cary SC, Williamson KE, Helton RR, Bench SR, Winget D, et al. Lysogenic virus-host interactions predominate at deep-sea diffuse-flow hydrothermal vents. *ISME J.* 2008;2(11):1112–21.

23. Anderson RE, Brazelton WJ, Baross JA. Is the genetic landscape of the deep subsurface biosphere affected by viruses? *Frontiers in microbiology.* 2011;2:219.

24. Anderson RE, Sogin ML, Baross JA. Evolutionary strategies of viruses, bacteria and archaea in hydrothermal vent ecosystems revealed through metagenomics. *PLoS one.* 2014;9:e109696.

25. Anantharaman K, Duhaime MB, Breier JA, Wendt KA, Toner BM, Dick GJ. Sulfur oxidation genes in diverse deep-sea viruses. *Science.* 2014;344(6185):757.

26. Ahlgren NA, Fuchsman CA, Rocab G, Fuhrman JA. Discovery of several novel, widespread, and ecologically distinct marine Thaumarchaeota viruses that encode amoC nitrification genes. *ISME J.* 2019;13(3):618–31.

27. Lossouarn J, Dupont S, Gorlas A, Mercier C, Bienvenu N, Marguet E, et al. An abyssal mobilome: viruses, plasmids and vesicles from deep-sea hydrothermal vents. *Res. Microbiol.* 2015;166(10):742–52.
28. Thiroux S, Dupont S, Nesbø CL, Bienvenu N, Krupovic M, L'Haridon S, et al. The first head-tailed virus, MFTV1, infecting hyperthermophilic methanogenic deep-sea archaea. *Environ. Microbiol.* 2021;23(7):3614–26.

29. Nigro OD, Jungbluth SP, Lin H-T, Hsieh C-C, Miranda JA, Schvarcz CR, et al. Viruses in the oceanic basement. *mBio*. 2017;8.

30. Castelán-Sánchez HG, Lopéz-Rosas I, García-Suastegui WA, Peralta R, Dobson ADW, Batista-García RA, et al. Extremophile deep-sea viral communities from hydrothermal vents: structural and functional analysis. *Mar. Genomics*. 2019;46:16–28.

31. Thomas E, Anderson RE, Li V, Rogan LJ, Huber JA, Petersen JM, et al. Diverse viruses in deep-sea hydrothermal vent fluids have restricted dispersal across ocean basins. *mSystems*. 2021;6(3):e00068–21.

32. Pratama AA, Bolduc B, Zayed AA, Zhong ZP, Guo J, Vik DR, et al. Expanding standards in viromics: in silico evaluation of dsDNA viral genome identification, classification, and auxiliary metabolic gene curation. *PeerJ*. 2021;9:e11447.

33. Williamson SJ, Allen LZ, Lorenzi HA, Fadrosh DW, Brami D, Thiagarajan M, et al. Metagenomic exploration of viruses throughout the Indian Ocean. *PLoS One*. 2012;7(10):e42047.

34. Roux S, Hallam SJ, Woyke T, Sullivan MB. Viral dark matter and virus-host interactions resolved from publicly available microbial genomes. *eLife*. 2015;4:e08490.

35. Paez-Espino D, Eloe-Fadrosh EA, Pavlopoulos GA, Thomas AD, Huntemann M, Mikhailova N, et al. Uncovering Earth's virome. *Nature*. 2016;536(7617):425–30.

36. Gruber-Vodicka HR, BKB S, Pruesse E. phyloFlash: rapid small-subunit rRNA profiling and targeted assembly from metagenomes. *mSystems*. 2020;5(5):e00920.

37. Waite DW, Vanwonterghem I, Rinke C, Parks DH, Zhang Y, Takai K, et al. Comparative genomic analysis of the class Epsilonproteobacteria and proposed reclassification to Epsilonbacteraea (phyl. nov.). *Front. Microbiol*. 2017;8.

38. Waite DW, Vanwonterghem I, Rinke C, Parks DH, Zhang Y, Takai K, et al. Erratum: addendum: comparative genomic analysis of the class Epsilonproteobacteria and proposed reclassification to Epsilonbacteraea (phyl. nov.). *Front. Microbiol*. 2018;9:772.

39. Bowers RM, Kyrpides NC, Stepanauskas R, Harmon-Smith M, Doud D, Reddy TBK, et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. *Nat. Biotechnol*. 2017;35(8):725–31.

40. Roux S, Enault F, Hurwitz BL, Sullivan MB. VirSorter: mining viral signal from microbial genomic data. *PeerJ*. 2015;3:e985.

41. Kieft K, Zhou Z, Anantharaman K. VIBRANT: automated recovery, annotation and curation of microbial viruses, and evaluation of viral community function from genomic sequences. *Microbiome*. 2020;8(1):90.

42. Al-Shayeb B, Sachdeva R, Chen L-X, Ward F, Munk P, Devoto A, et al. Clades of huge phages from across Earth's ecosystems. *Nature*. 2020;578(7795):425–31.

43. Koonin EV, Yutin N. Origin and evolution of eukaryotic large nucleocytoplasmic DNA viruses. *Intervirology*. 2010;53(5):284–92.

44. Nayfach S, Camargo AP, Schulz F, Eloe-Fadrosh E, Roux S, Kyrpides NC. CheckV assesses the quality and completeness of metagenome-assembled viral genomes. *Nat. Biotechnol*. 2021;39(5):578–85.

45. Maniloff J, Ackermann HW. Taxonomy of bacterial viruses: establishment of tailed virus genera and the other Caudovirales. *Arch. Virol*. 1998;143(10):2051–63.

46. Krupovic M. Networks of evolutionary interactions underlying the polyphyletic origin of ssDNA viruses. *Curr. Opin. Virol*. 2013;3(5):578–86.

47. King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ. Virus taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses. San Diego: Elsevier Academic Press; 2012.

48. Zhan Y, Chen F. The smallest ssDNA phage infecting a marine bacterium. *Environ. Microbiol*. 2019;21(6):1916–28.

49. Zhao L, Rosario K, Breitbart M, Duffy S. Chapter three - Eukaryotic circular rep-encoding single-stranded DNA (CRESS DNA) viruses: ubiquitous viruses with small genomes and a diverse host range. In: Kielian M, Mettenleiter TC, Roossinck MJ, editors. Advances in Virus Research, vol. 103: Academic Press; 2019. p. 71–133.

50. Dean FB, Nelson JR, Giesler TL, Lasken RS. Rapid amplification of plasmid and phage DNA using Phi29 DNA polymerase and multiply-primed rolling circle amplification. *Genome Res.* 2001;11(6):1095–9.

51. Binga EK, Lasken RS, Neufeld JD. Something from (almost) nothing: the impact of multiple displacement amplification on microbial ecology. *ISME J*. 2008;2(3):233–41.

52. Yoshida M, Mochizuki T, Urayama S-I, Yoshida-Takashima Y, Nishi S, Hirai M, et al. Quantitative viral community DNA analysis reveals the dominance of single-stranded DNA viruses in offshore Upper Bathyal sediment from Tohoku, Japan. *Front. Microbiol*. 2018;9(75).

53. Endo H, Blanc-Mathieu R, Li Y, Salazar G, Henry N, Labadie K, et al. Biogeography of marine giant viruses reveals their interplay with eukaryotes and ecological functions. *Nat. Ecol. Evol.* 2020;4(12):1639–49.

54. Kristensen DM, Mushegian AR, Dolja VV, Koonin EV. New dimensions of the virus world discovered through metagenomics. *Trends Microbiol*. 2010;18(11):1–9.

55. Jian H, Yi Y, Wang J, Hao Y, Zhang M, Wang S, et al. Diversity and distribution of viruses inhabiting the deepest ocean on Earth. *ISME J*. 2021;15(10):3094–110.

56. Gregory AC, Zayed AA, Conceição-Neto N, Temperton B, Bolduc B, Alberti A, et al. Marine DNA viral macro- and microdiversity from pole to pole. *Cell*. 2019;177(5):1109–1123.e1114.

57. Li Z, Pan D, Wei G, Pi W, Zhang C, Wang J-H, et al. Deep sea sediments associated with cold seeps are a subsurface reservoir of viral diversity. *ISME J*. 2021.

58. Rinke C, Rubino F, Messer LF, Youssef N, Parks DH, Chuvochina M, et al. A phylogenomic and ecological analysis of the globally abundant Marine group II archaea (Ca. Poseidonia ord. nov.). *ISME J*. 2019;13(3):663–75.

59. Rinke C, Rubino F, Messer LF, Youssef N, Parks DH, Chuvochina M, et al. A phylogenomic and ecological analysis of the globally abundant Marine group II archaea (Ca. Poseidonia ord. nov.). *ISME J*. 2019;13(3):663–75.

60. Zhang CL, Xie W, Martin-Cuadrado A-B, Rodriguez-Valera F. Marine group II archaea, potentially important players in the global ocean carbon cycle. *Frontiers in microbiology*. 2015;6:1108.

61. Haro-Moreno JM, Rodriguez-Valera F, López-García P, Moreira D, Martin-Cuadrado A-B. New insights into marine group III Euryarchaeota, from dark to light. *ISME J*. 2017;11(5):1102–17.

62. Li M, Baker BJ, Anantharaman K, Jain S, Breier JA, Dick GJ. Genomic and transcriptomic evidence for scavenging of diverse organic compounds by widespread deep-sea archaea. *Nat. Commun.* 2015;6:8933.

63. Ding J, Zhang Y, Wang H, Jian H, Leng H, Xiao X. Microbial community structure of deep-sea hydrothermal vents on the ultraslow spreading Southwest Indian Ridge. *Front. Microbiol*. 2017;8(1012).

64. Thingstad TF, Lignell R. Theoretical models for the control of bacterial growth rate, abundance, diversity and carbon demand. *Aquat. Microb. Ecol.* 1997;13(1):19–27.

65. Zhao Y, Temperton B, Thrash JC, Schwalbach MS, Vergin KL, Landry ZC, et al. Abundant SAR11 viruses in the ocean. *Nature*. 2013;494(7437):357–60.

66. Fortunato CS, Larson B, Butterfield DA, Huber JA. Spatially distinct, temporally stable microbial populations mediate biogeochemical cycling at and below the seafloor in hydrothermal vent fluids. *Environ. Microbiol*. 2018;20(2):769–84.

67. Huber JA, Cantin HV, Huse SM, Mark Welch DB, Sogin ML, Butterfield DA. Isolated communities of Epsilonproteobacteria in hydrothermal vent fluids of the Mariana arc seamounts. *FEMS Microbiol. Ecol.* 2010;73(3):538–49.

68. Campbell BJ, Smith JL, Hanson TE, Klotz MG, Stein LY, Lee CK, et al. Adaptations to submarine hydrothermal environments exemplified by the genome of Nautilia profundicola. *PLoS Genet.* 2009;5(2):e1000362.

69. Nakagawa S, Takaki Y, Shimamura S, Reysenbach A-L, Takai K, Horikoshi K. Deep-sea vent e-proteobacterial genomes provide insights into emergence of pathogens. *Proc. Natl. Acad. Sci.* 2007;104(29):12146.

70. Yoshida-Takashima Y, Takaki Y, Shimamura S, Nunoura T, Takai K. Genome sequence of a novel deep-sea vent epsilonproteobacterial phage provides new insight into the co-evolution of Epsilonproteobacteria and their phages. *Extremophiles: life under extreme conditions*. 2013;17(3):405–19.
71. Díez-Villaseñor C, Rodriguez-Valera F. CRISPR analysis suggests that small circular single-stranded DNA smacoviruses infect archaea instead of humans. *Nat. Commun.* 2019;10(1):294.

72. Yoshida M, Takaki Y, Eitoku M, Nunoura T, Takai K. Metagenomic analysis of viral communities in (hado)pelagic sediments. *PLoS One*. 2013;8(2):e57271.

73. Martin DP, Biagini P, Lefeuvre P, Golden M, Roumagnac P, Varsani A. Recombination in eukaryotic single stranded DNA viruses. *Viruses*. 2011;3(9):1699–738.

74. Krupovic M, Ravantti JJ, Bamford DH. Geminiviruses: a tale of a plasmid becoming a virus. *BMC Evol. Biol.* 2009;9(1):112.

75. Kazlauskas D, Varsani A, Koonin EV, Krupovic M. Multiple origins of prokaryotic and eukaryotic single-stranded DNA viruses from bacterial and archaeal plasmids. *Nat. Commun.* 2019;10(1):3425.

76. Warwick-Dugdale J, Buchholz HH, Allen MJ, Temperton B. Host-hijacking and planktonic piracy: how phages command the microbial high seas. *Virol. J.* 2019;16(1):15.

77. Kieft K, Zhou Z, Anderson RE, Buchan A, Campbell BJ, Hallam SJ, et al. Ecology of inorganic sulfur auxiliary metabolism in widespread bacteriophages. *Nat. Commun.* 2021;12(1):3503.

78. Frank KL, Rogers DR, Olins HC, Vidoudez C, Girguis PR. Characterizing the distribution and rates of microbial sulfate reduction at Middle Valley hydrothermal vents. *ISME J.* 2013;7(7):1391–401.

79. Ikeuchi Y, Shigi N, Kato J-i, Nishimura A, Suzuki T: Mechanistic insights into sulfur relay by multiple sulfur mediators involved in thiouridine biosynthesis at tRNA wobble positions. *Mol. Cell.* 2006;21(1):97–108.

80. Murata M, Fujimoto H, Nishimura K, Charoensuk K, Nagamitsu H, Raina S, et al. Molecular strategy for survival at a critical high temperature in Eschierichia coli. *PLoS One*. 2011;6(6):e20063.

81. Shigi N, Sakaguchi Y, Suzuki T, Watanabe K. Identification of two tRNA thiolation genes required for cell growth at extremely high temperatures. *J. Biol. Chem.* 2006;281(20):14296–306.

82. Alvarez-Ordóñez A, Cummins C, Deasy T, Clifford T, Begley M, Hill C. Acid stress management by Cronobacter sakazakii. *Int. J. Food Microbiol.* 2014;178:21–8.

83. Cao L, Wang J, Sun L, Kong Z, Wu Q, Wang Z. Transcriptional analysis reveals the relativity of acid tolerance and antimicrobial peptide resistance of Salmonella. *Microb. Pathog.* 2019;136:103701.

84. Gang H, Xiao C, Xiao Y, Yan W, Bai R, Ding R, et al. Proteomic analysis of the reduction and resistance mechanisms of Shewanella oneidensis MR-1 under long-term hexavalent chromium stress. *Environ. Int.* 2019;127:94–102.

85. Liu J, Chen L, Wang J, Qiao J, Zhang W. Proteomic analysis reveals resistance mechanism against biofuel hexane in Synechocystis sp. PCC 6803. *Biotechnol. Biofuels.* 2012;5(1):68.

86. Shaffer M, Borton MA, McGivern BB, Zayed AA, La Rosa SL, Solden LM, et al. DRAM for distilling microbial metabolism to automate the curation of microbiome function. *Nucleic Acids Res.* 2020;48(16):8883–900.

87. James MN. The peptidases from fungi and viruses. *Biol. Chem.* 2006;387(8):1023–9.

88. Jeudy S, Rigou S, Alempic J-M, Claverie J-M, Abergel C, Legendre M. The DNA methylation landscape of giant viruses. *Nat. Commun.* 2020;11(1):2657.

89. Stepanauskas R, Fergusson EA, Brown J, Poulton NJ, Tupper B, Labonté JM, et al. Improved genome recovery and integrated cell-size analyses of individual uncultured microbial cells and viral particles. *Nat. Commun.* 2017;8(1):84.

90. Chen S, Zhou Y, Chen Y. Gu J: fastp: an ultra-fast all-in-one FASTQ preprocessor. *Bioinformatics*. 2018;34(17):i884–90.

91. Li D, Liu C-M, Luo R, Sadakane K, Lam T-W. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. *Bioinformatics*. 2015;31(10):1674–6.

92. Hideki N, Jungho P, Toshihisa T. MetaGene: prokaryotic gene finding from environmental genome shotgun sequences. *Nucleic Acids Research.* 2006;34(19):5623–30.

93. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto Encyclopedia of Genes and Genomes. *Nucleic Acids Res.* 1999;27(1):29–34.

94. Huerta-Cepas J, Szklarczyk D, Heller D, Hernández-Plaza A, Forslund SK, Cook H, et al. eggNOG 5.0: a hierarchical, functionally and

phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. *Nucleic Acids Res.* 2018;47(D1):D309–14.

95. Vincent L, Hemalatha GR, Elodie D, Coutinho PM, Bernard H. The carbohydrate-active enzymes database (CAZy) in 2013. *Nucleic Acids Research*. 2014;D1:D490.

96. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. *J. Mol. Biol.* 1990;215(3):403–10.

97. Kang DD, Li F, Kirton E, Thomas A, Egan R, An H, et al. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. *PeerJ*. 2019;7:e7359.

98. Uritskiy GV, DiRuggiero J, Taylor J. MetaWRAP-a flexible pipeline for genome-resolved metagenomic data analysis. *Microbiome*. 2018;6(1):158.

99. Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. *Genome Res.* 2015;25(7):1043–55.

100. Olm MR, Brown CT, Brooks B, Banfield JF. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. *ISME J.* 2017;11(12):2864–8.

101. Chaumeil P-A, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. *Bioinformatics*. 2019;36(6):1925–7.

102. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. *Mol. Biol. Evol.* 2020;37(5):1530–4.

103. Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS. UFBoot2: improving the ultrafast bootstrap approximation. *Mol. Biol. Evol.* 2017;35(2):518–22.

104. Ivica L. Peer, Bork: Interactive Tree of Life (iTOL) v4: recent updates and new developments. *Nucleic Acids Res.* 2019.

105. Clustering viral genomes in iVirus. [https://www.protocols.io/view/clustering-viral-genomes-in-ivirus-gwebxbe].

106. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. *Nucleic Acids Res.* 2012;41(D1):D590–6.

107. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. *Mol. Plant.* 2020;13(8):1194–202.

108. Bland C, Ramsey TL, Sabree F, Lowe M, Brown K, Kyrpides NC, et al. CRISPR Recognition Tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats. *BMC Bioinformatics.* 2007;8(1):209.

109. Skennerton CT, Imelfort M, Tyson GW. Crass: identification and reconstruction of CRISPR from unassembled metagenomic data. *Nucleic Acids Res.* 2013;41(10):e105.

110. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. *BMC Bioinformatics.* 2009;10(1):421.

111. Emerson JB, Roux S, Brum JR, Bolduc B, Woodcroft BJ, Jang HB, et al. Host-linked soil viral ecology along a permafrost thaw gradient. *Nat. Microbiol.* 2018;3(8):870–80.

112. Laslett D, Canback B. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. *Nucleic Acids Res.* 2004;32(11):11–6.

113. Tominaga K, Morimoto D, Nishimura Y, Ogata H, Yoshida T. In silico prediction of virus-host interactions for marine bacteroidetes with the use of metagenome-assembled genomes. *Front. Microbiol.* 2020;11(738).

114. Galiez C, Siebert M, Enault F, Vincent J, Söding J. WISH: who is the host? Predicting prokaryotic hosts from metagenomic phage contigs. *Bioinformatics.* 2017;33(19):3113–4.

115. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. *Nat. Methods.* 2015;12(1):59–60.

116. von Meijenfeldt FAB, Arkhipova K, Cambuy DD, Coutinho FH, Dutilh BE. Robust taxonomic classification of uncharted microbial sequences and bins with CAT and BAT. *Genome Biol.* 2019;20(1):217.

117. Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. *BMC Bioinformatics.* 2010;11(1):119.

118. Robert C. Edgar: MUSCLE: multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Res.* 2004.
119. Capella-Gutiérrez S, Silla-Martínez JM. Gabaldón T: trimAI: a tool for automated alignment trimming in large-scale phylogenetic analyses. *Bioinformatics*. 2009;25(15):1972-3.

120. Kalyaanamoorthy S, Minh BQ, Wong TKF, Haeseler AV, Jermiin LS. Modelfinder: fast model selection for accurate phylogenetic estimates. *Nat. Methods*. 2017;14(6).

121. Cantalapiedra CP, Hernández-Plaza A, Letunic I, Bork P, Huerta-Cepas J. eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale. *bioRxiv*. 2021.

122. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. *Nat. Methods*. 2012;9(4):357-9.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_10.png" /></figure>

BMC