# Diversity, taxonomy, and evolution of archaeal viruses of the class Caudoviricetes

Ying Liu, Tatiana A Demina, Simon Roux, Pakorn Aiewsakun, Darius

Kazlauskas, Peter Simmonds, David Prangishvili, Hanna M Oksanen, Mart

Krupovic

► To cite this version:

Ying Liu, Tatiana A Demina, Simon Roux, Pakorn Aiewsakun, Darius Kazlauskas, et al., Diversity, taxonomy, and evolution of archaeal viruses of the class Caudoviricetes. PLoS Biology, 2021, 19 (11), pp.e3001442. 10.1371/journal.pbio.3001442. pasteur-03471084

HAL Id: pasteur-03471084

https://pasteur.hal.science/pasteur-03471084

Submitted on 8 Dec 2021

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 ARTICLE

# Diversity, taxonomy, and evolution of archaeal viruses of the class *Caudoviricetes*

Ying Liu<sup>1,*</sup>, Tatiana A. Demina<sup>2,*</sup>, Simon Roux<sup>3</sup>, Pakorn Aiewsakun<sup>4,5,6</sup>, Darius Kazlauskas<sup>7</sup>, Peter Simmonds<sup>4</sup>, David Prangishvili<sup>1,8</sup>, Hanna M. Oksanen<sup>2,*</sup>, Mart Krupovic<sup>1*</sup>

1 Institut Pasteur, Université de Paris, Archaeal Virology Unit, Paris, France, 2 Molecular and Integrative Biosciences Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland, 3 DOE Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, California, United States of America, 4 Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom, 5 Department of Microbiology, Faculty of Science, Mahidol University, Bangkok, Thailand, 6 Pornchai Matangkasombut Center for Microbial Genomics, Department of Microbiology, Faculty of Science, Mahidol University, Bangkok, Thailand, 7 Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania, 8 Ivane Javakhishvili Tbilisi State University, Tbilisi, Georgia

© These authors contributed equally to this work.

* hanna.oksanen@helsinki.fi (HMO); mart.krupovic@pasteur.fr (MK)

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

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

## OPEN ACCESS

**Citation:** Liu Y, Demina TA, Roux S, Aiewsakun P, Kazlauskas D, Simmonds P, et al. (2021) Diversity, taxonomy, and evolution of archaeal viruses of the class *Caudoviricetes*. PLoS Biol 19(11): e3001442. https://doi.org/10.1371/journal.pbio.3001442

**Academic Editor:** Curtis Suttle, University of British
Colombia, CANADA

**Received:** July 1, 2021

**Accepted:** October 17, 2021

**Published:** November 9, 2021

**Peer Review History:** PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.pbio.3001442

**Copyright:** © 2021 Liu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All new genome sequences of virus isolates were submitted to GenBank (accession numbers MZ334492-MZ334528). All protein sequences used for phylogenetic analyses can be downloaded from the

## Abstract

The archaeal tailed viruses (arTV), evolutionarily related to tailed double-stranded DNA (dsDNA) bacteriophages of the class *Caudoviricetes*, represent the most common isolates infecting halophilic archaea. Only a handful of these viruses have been genomically characterized, limiting our appreciation of their ecological impacts and evolution. Here, we present 37 new genomes of haloarchaeal tailed virus isolates, more than doubling the current number of sequenced arTVs. Analysis of all 63 available complete genomes of arTVs, which we propose to classify into 14 new families and 3 orders, suggests ancient divergence of archaeal and bacterial tailed viruses and points to an extensive sharing of genes involved in DNA metabolism and counterdefense mechanisms, illuminating common strategies of virus-host interactions with tailed bacteriophages. Coupling of the comparative genomics with the host range analysis on a broad panel of haloarchaeal species uncovered 4 distinct groups of viral tail fiber adhesins controlling the host range expansion. The survey of metagenomes using viral hallmark genes suggests that the global architecture of the arTV community is shaped through recurrent transfers between different biomes, including hypersaline, marine, and anoxic environments.

## Introduction

Bacteriophages with helical tails and icosahedral capsids (tailed bacteriophages), classified into the class *Caudoviricetes* [1], represent the most widespread, abundant, and diverse group of viruses on our planet [2] and are likely to infect hosts from most, if not all, known bacterial lineages [3]. Extensive experimental studies conducted over several decades, coupled with comparative analysis of several thousands of complete tailed phage genomes currently available in
IMG/M database at JGI using the accession numbers listed in S2 Data.

**Funding:** This work was supported by l'Agence Nationale de la Recherche grant ANR-20-CE20-0009-02 (to M.K.) and the European Union's Horizon 2020 research and innovation program under grant agreement 685778, project VIRUS X (to D.P.). Y.L. is a recipient of the Pasteur-Roux-Cantarini Fellowship from Institut Pasteur. The Ella and Georg Ehrnrooth Foundation and the Finnish Cultural Foundation are sincerely acknowledged (grants to T.D.). The work conducted by the U.S. Department of Energy Joint Genome Institute (S. R.) is supported by the Office of Science of the U.S. Department of Energy under contract no. DE-AC02-05CH11231. The development and application of GRAViTy analysis was supported by a grant to PS from the Wellcome Trust (WT108418A1A). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

**Abbreviations:** arCOG, archaeal clusters of orthologous genes; arTV, archaeal tailed virus; CGJ, composite generalized Jaccard; CP, carbamoylphosphate; dsDNA, double-stranded DNA; EC, Executive Committee; EOP, efficiency of plating; HMM, hidden Markov model; ICTV, International Committee on Taxonomy of Viruses; IMG/M, Integrated Microbial Genomes and Microbiomes; MCP, major capsid protein; MGM, modified growth medium; MIDAS, metal ion-dependent adhesion site; MTase, methyltransferase; ORF, open reading frame; REase, restriction endonuclease; RM, restriction-modification; RnR, ribonucleotide diphosphate reductase; SW, salt water; TA, toxin-antitoxin; vWA, von Willebrand factor type A.

the public sequence databases, have resulted in detailed understanding of the mechanisms that govern the biology, ecology, and evolution of this group of viruses [2,4–6]. Due to their ubiquity, tailed bacteriophages have a profound impact on the functioning of the biosphere through regulating the structure, composition, and dynamics of bacterial populations in diverse environments, from marine ecosystems to the human gut, and modulate major biogeochemical cycles [7–9]. Archaeal tailed viruses (arTV) are morphologically indistinguishable from tailed bacteriophages [1,4,10–12]. Similar to their bacterial virus relatives, the helical tails of arTVs can be short (podovirus morphology), long noncontractile (siphovirus morphology), or contractile (myovirus morphology) [13–15].

The arTVs have been thus far isolated on halophilic (class Halobacteria) and methanogenic (family Methanobacteriaceae) archaea, both belonging to the phylum Euryarchaeota [16–25]. Related proviruses have also been sighted in other lineages of the Euryarchaeota as well as in ammonia-oxidizing Thaumarchaeota and Aigarchaeota [26–30], whereas recent metagenomics studies revealed novel groups of arTVs putatively infecting marine group II Euryarchaeota, Thaumarchaeota, and Thermoplasmata [31–35]. Ecologically, it has been shown that virus-mediated lysis of archaea in the deep ocean is more rapid than that of bacteria, suggesting an important ecological role of archaeal viruses in marine ecosystems [36]. Evolutionarily, the broad distribution of tailed (pro)viruses in both bacteria and archaea suggests that viruses of this type were part of the virome associated with the last universal cellular ancestor, LUCA [3].

Genomic and structural analyses have shown that archaeal and bacterial tailed viruses have similar genomic organizations, with genes clustered into functional modules, and share homologous virion morphogenesis modules, including the major capsid protein (MCP) with the characteristic HK97 fold and the genome packaging terminase complex, suggesting common principles of virion assembly [11,13,22,30,37,38]. Nevertheless, at the sequence level, arTVs are strikingly diverse showing little similarity to each other and virtually no recognizable similarity to their bacterial virus relatives, indicative of scarce sampling of the archaeal virosphere [38]. Indeed, for several thousands of complete genome sequences of tailed bacteriophages [2], only 25 arTV isolates have been sequenced thus far. The low number of isolates severely limits our appreciation of the ecological impacts of arTVs and obscures the evolutionary history of this important and ancient group of viruses.

Virus discovery in the “age of metagenomics” is increasingly performed by culture-independent methods, whereby viral genomes are sequenced directly from the environment. This is a powerful approach that has already yielded thousands of viral genomes, providing unprecedented insights into virus diversity, environmental distribution, and evolution [39–46]. The limitation of viral metagenomics, however, is that the exact host species for the sequenced viruses typically remain unknown, and many molecular aspects of virus-host interactions cannot be accurately predicted. Here, to further explore the biology and diversity of arTVs, we sequenced the genomes of 37 viruses that infect different species of halophilic archaea and were isolated from hypersaline environments using classical approaches [16,17]. Collectively, our results provide the first global overview of arTV diversity and evolution and establish a taxonomic framework for their classification.

# Results and discussion

## Overview of new haloarchaeal tailed viruses

We sequenced a collection of 37 arTVs (5 siphoviruses and 32 myoviruses) infecting haloarchaeal species belonging to the genera *Halorubrum* and *Haloarcula* [16,17], more than doubling the number of complete genomes of arTVs (S1 Table). The sequenced viruses originate from
geographically remote locations, including Thailand, Israel, Italy, and Slovenia, and in combi-
nation with the previously described isolates, provide a substantially improved genomic insight
into the global distribution of arTVs. The viruses possess double-stranded DNA (dsDNA)
genomes ranging from 35.3 to 104.7 kbp in length. Genomes of several isolates were nearly
identical (<17 nucleotide polymorphisms; S1 Table) and analysis of these genomes, in combi-
nation with host range experiments, was particularly illuminating toward the host range evolu-
tion among halophilic arTVs (see below).

Archaeal tailed viruses represent a distinct group within the prokaryotic
virosphere

To assess the global diversity of arTVs and analyze their relationship to bacterial virus mem-
bers of the class Caudoviricetes, we supplemented the 37 genomes sequenced herein with the
genomes of prokaryotic viruses available in GenBank and analyzed the dataset using GRAViTy
[47,48] and vConTACT v2.0 [49]. The combined dataset included 63 complete arTV genomes
(S1 Table), 3 of which were from viruses infecting methanogenic archaea, whereas all others
were from haloarchaeal viruses. The GRAViTy tool classifies viruses into family-level taxo-
nomic groupings according to homology between viral genes and similarities in genome orga-
nizations, which are expressed using composite generalized Jaccard (CGJ) distances [47,48].
We used a CGJ distance of over 0.8 as the threshold for family-level assignment, consistent
with the family-level classification for eukaryotic viruses and the recently created families of
bacterial viruses [47,48,50]. GRAViTy analysis of the global prokaryotic virome classified
arTVs into 2 large assemblages, which could be further subdivided into 14 family-level group-
ings (CGJ distance ≥ 0.8) (Fig 1). To reveal a finer taxonomic structure within the arTV
assemblage, we relied on the network analytics implemented in vConTACT v2.0, which has
been specifically developed and calibrated to identify genus-level groupings of prokaryotic
viruses [49]. Consistent with the GRAViTy results, the network analysis revealed 2 assem-
blages of arTVs, which were disconnected from all known bacterial and nontailed archaeal
viruses (Fig 2). Notably, only a few genes were shared between arTVs and nontailed archaeal
viruses, including those encoding a methyltransferase (MTase), transposase, adenylyltransfer-
ase, AAA+ domain protein, and proteins of unknown function (see spreadsheet “Fig 2 PCs” in
S2 Data). Viruses within the 2 clades formed 23 genus-level and 14 family-level groups, many
containing just 1 or 2 members, indicating that genetic diversity of archaeal viruses remains
largely undersampled. The family-level and genus-level groupings are described in S1 Text.

Clade I in the GRAViTy analysis forms a sister branch to several groups of bacterial myo-
viruses, including families *Ackermannviridae*, *Herelleviridae*, and *Straboviridae* (T4-like bacte-
riophages), and consists of 4 family-level groups (Fig 1A and 1B, S1 Text), with family (F) 1
being by far the largest, with 39 members, which can be further divided into 4 genus-level (G)
subgroups (Fig 2). Clade I is a cohesive assemblage consisting exclusively of haloarchaeal
viruses and held together by 33 protein clusters, including those involved in virion morpho-
genesis, genome replication and repair, nucleotide metabolism, and several proteins of
unknown function (S2 Table). By contrast, Clade II is less cohesive and comprises 10 family-
level groupings, half of which consist of singletons (Fig 1A and 1C), and includes viruses of
both halophilic and methanogenic archaea. Viruses from different family-level groups in this
clade share an overlapping set of genes (52 protein clusters), including those involved in virion
morphogenesis, genome replication, and repair, DNA metabolism and several other function-
ally diverse proteins. However, unlike in Cluster I, the set of genes holding the Cluster II
together displays sporadic distribution, typically connecting 2 to 3 virus families (S2 Table).
<figure><img src="image_4.png" /></figure>

B Clade I

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

C Clade II

<figure><img src="image_6.png" /><figcaption>**Fig 1. Genome relationships between prokaryotic dsDNA viruses.** (A) Heat map and dendrogram of CGJ distances for classified bacteriophages and archaeal dsDNA viruses. Branches and clusters corresponding to arTVs are shown in red, whereas those of other viruses are in blue. (B) Zoom in on the arTV Clade I. (C) Zoom in on the arTV Clade II. Viruses sequenced in this study are highlighted in bold. CGJ distance of 0.8, chosen as a family-level threshold, is indicated with a broken line, with the family-level groups (F1 to F14) indicated on the left of the heat map. See S1 Data for the underlying dendrogram. arTV, archaeal tailed virus; CGJ, composite generalized Jaccard; dsDNA, double-stranded DNA.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001442.g001
<figure><img src="image_7.png" /><figcaption>**Fig 2.** The network-based analysis of PCs shared among arTVs and the prokaryotic dsDNA viruses. The nodes represent viral genomes, and the edges represent the strength of connectivity between each genome based on shared PCs. Nodes representing genomes of arTVs are in red and other archaeal viruses are in pink, whereas those representing genomes of tailed bacteriophages are in blue and other bacteriophages are in light blue (left panel). The VCs of arTVs are enlarged and labeled in the right panel. See S2 Data for complete VCs and PCs generated by vConTACT v2.0. arTV, archaeal tailed virus; dsDNA, double-stranded DNA; PC, protein cluster; VC, viral cluster.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001442.g002

Based on the results of comprehensive comparative genomics analysis (see S1 Text), we propose classifying all known arTVs into 14 new families. The families *Hafunaviridae*, *Soleiviridae*, *Halomagnusviridae*, and *Pyrstoviridae* include viruses with icosahedral heads and long contractile tails (myovirus morphotype), whereas the families *Druskaviridae*, *Haloferuviridae*, *Graaviviridae*, *Vertoviridae*, *Suolaviridae*, *Saparoviridae*, *Madisaviridae*, *Leisingerviridae*, and *Anaerodiviridae* contain viruses with icosahedral heads and long noncontractile tails (sipho-virus morphotype). The *Shortaselviridae* is the only family of viruses with short tails (podo-virus morphology). The names of the 23 proposed genera are listed in S1 Table. Members of the same proposed genus typically share more than 60% of their proteins, and members of the same family share 20% to 50% of homologous proteins, whereas viruses from different families share less than 10% of proteins (Fig 3, S2 Fig). The 14 families have been approved by the Executive Committee (EC) of the International Committee on Taxonomy of Viruses (ICTV) and await ratification by the whole ICTV membership. Furthermore, the ICTV EC has approved unification of the 4 Clade I families into a new order *Thumleimavirales*; families *Haloferuviridae*, *Pyrstoviridae*, *Shortaselviridae*, and *Graaviviridae* have been included into a new order *Kirjokansivirales*; and families *Leisingerviridae* and *Anaerodiviridae* were included into a new order *Methanobavirales*.

## Gene content of archaeal tailed viruses

All virus genes were functionally annotated using sensitive profile-profile hidden Markov model (HMM) comparisons using HHpred [51] (S3 Table). The virus-encoded proteins were further classified into functional categories based on their affiliation to the archaeal clusters of orthologous genes (arCOG) [52] (S3 Table). Apart from the “Virus-related” proteins, the most frequent functional category assigned to archaeal virus proteins was the “Information storage and processing,” with the “Genome replication, recombination, and repair” subcategory being
<figure><img src="image_8.png" /><figcaption>**Fig 3.** Overview of homologous proteins shared by arTVs. (A) The circus plot displays gene content similarity between representative members of each proposed family. Viral genomes from different families are indicated with distinct colors. The genomic coordinates (kbp) are indicated on genome maps. Putative ORFs distributed on 5 tracks (to avoid overlapping) are represented by green and red tiles on the forward and reverse strands, respectively. Homologous proteins are linked by gray lines. The representative viruses and the proposed virus family names are indicated. (B) The box plot shows the percentage of genes shared by a representative virus from each family with other members of the same genus (G) and family (F) as well as with arTVs from other families (A). Each box represents the middle 50th percentile of the data set and is derived using the lower and upper quartile values. The median value is displayed by a horizontal line. Whiskers represent the maximum and minimum values with the range of 1.5 IQR. Each virus is represented by a dot. In the circus and box plots analyses, proteins with over 30% amino acid sequence identity and E-value < 1 × 10⁻²⁵ in our arTV database are considered to be homologous. Data underlying this figure can be found in S2 Data. arTV, archaeal tailed virus.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001442.g003
<figure><img src="image_9.png" /><figcaption>**Fig 4.** Classification of genes encoded by arTVs. (A) Classification of arTV genes into arCOG functional categories. The homologous gene shared by viruses in the same species is counted as one. The letters represent the following: J, translation, ribosomal structure and biogenesis; K, transcription; L, replication, recombination and repair; D, cell cycle control, cell division, chromosome partitioning; V, defense mechanisms; T, signal transduction mechanisms; M, cell wall/membrane/envelope biogenesis; N, cell motility; O, posttranslational modification, protein turnover, chaperons; G, carbohydrate transport and metabolism; E, amino acid transport and metabolism; F, nucleotide transport and metabolism; H, coenzyme transport and metabolism; I, lipid transport and metabolism; P, inorganic ion transport and metabolism; Q, secondary metabolites biosynthesis, transport and catabolism; R, general function prediction only; S, function unknown; X, virus related. (B) Schematic showing nucleotide-related (partial) metabolic pathways, with a particular enzyme indicated for each reaction. Enzymes that are found encoded by arTVs are in red, otherwise in gray. The number of virus families that harbor the enzyme is indicated next to the enzyme name, with the number of virus species shown in parentheses. NrdA: RnR, MazG: NTP pyrophosphatase, 5'-NCT: 5'-deoxyribonucleotidase, CarA: CP synthase small subunit, Dcd: dCTP deaminase, Tmk: thymidylate kinase, THY1: thymidylate synthase thyX, NDT: nucleoside deoxyribosyltransferase, FolE: GTP cyclohydrolase I, QueC: queuosine biosynthesis protein, QueD: 6-pyruvoyl tetrahydropterin synthase, QueE: 7-carboxy-7-deazaguanine synthase, DpdA: paralog of queuine tRNA-ribosyltransferase, QueFC: NADPH-dependent 7-cyano-7-deazaguanine reductase. See S3 and S5 Tables for individual genes classification. arCOG, archaeal clusters of orthologous gene; arTV, archaeal tailed virus; CP, carbamoylphosphate; RnR, ribonucleotide diphosphate reductase.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001442.g004

most strongly enriched (Fig 4A). Proteins of the “Defense mechanisms” subcategory from the “Cellular processes and signaling” category, primarily including diverse nucleases and DNA MTases, were also abundant. Finally, a substantial fraction of proteins was assigned to the “Metabolism” category, with proteins involved in nucleotide transport and metabolism being most common (Fig 4A). Below we highlight some of the observations with the more detailed description provided in S1 Text.

## Defense and counterdefense factors

Similar to tailed bacteriophages, arTVs are under strong pressure to overcome the host defenses, among which CRISPR-Cas, restriction-modification (RM) and toxin-antitoxin (TA) are the ubiquitous defense systems in both archaea and bacteria. No clear homologs of viral anti-CRISPR proteins were detected in the arTVs genomes, although several members of the *Hafunaviridae* and *Druskaviridae* encode a Cas-4-like nuclease that could potentially function in counterdefense against CRISPR-Cas systems (see S1 Text). By contrast, haloarchaeal tailed viruses (n = 48) from 10 families encode diverse MTases with specificities for N6-adenine, N4-cytosine, and C5-cytosine (S4 Table, S3 Fig). The presence of the *MTase* genes can be linked to the presence of sequence motifs in the genomes of the corresponding viruses. For instance, GATC motif methylated by the phiCh1-like Dam MTase (G<sup>m6</sup>ATC) [53] is present at high frequency (6.47 to 11.04/kbp) in the genomes of viruses that encode this MTase. By contrast, viruses that do not encode the Dam MTase also lack the corresponding motif in their genomes. Thus, arTVs have likely evolved to escape the host RM systems by either self-methylation via different virus-encoded MTases or purging the recognition motifs from their
genomes. In addition to the stand-alone MTases, some arTVs encode accompanying restriction endonucleases (REases), forming apparently functional RM systems, which could be deployed by the viruses for degradation of the host chromosomes or genomes of competing mobile genetic elements (see S1 Text). A similar function could be envisioned for the RNA- or DNA-specific micrococcal nucleases encoded by HRTV-17, phiH1, HCTV-2, and HHTV-2 (S3 Table).

The 7-deazaguanine modifications, produced by the virus-encoded preQ0/G+ pathway, were recently detected in the genomes of diverse viruses, including HVTV-1, and shown to confer viral DNA with resistance to various type II REases [54]. The preQ0/G+ pathway is encoded by all viruses of the *Druskaviridae* (Fig 4B, S3 Table), suggesting that these viruses depend on a similar genome modification for evasion of the host RM systems. Notably, HCTV-1 and HCTV-16 in addition encode a transporter of the queuosine precursor YhhQ. Interestingly, HRTV-29 (family *Haloferuviridae*) does not carry genes for the preQ0/G+ pathway, but encodes a DpdA homolog, a key enzyme mediating replacement of the unmodified guanine base in the DNA (Fig 4B, S3 Table). DNA modification with solely virus-encoded DpdA was demonstrated for the *Salmonella* phage 7-11 [54], suggesting that HRTV-29 also hijacks the host preQ0 pathway for DNA modification through its DpdA protein.

The TA systems are widely distributed in prokaryotes and have been shown to function in bacterial antiviral defense by initiating the programmed cell death, thereby preventing the virus spread [55]. Similar to certain marine bacteriophages [56], viruses of the families *Druskaviridae*, *Vertoviridae*, and *Madisaviridae* encode homologs of the nucleoside pyrophosphohydrolase MazG (S5 Table), which prevents the programmed cell death by degrading the intracellular ppGpp [57]. In addition, several arTVs encode homologs of the VapB antitoxins (arCOG08550), but not the associated toxins, suggesting a function in blocking the VapBC TA systems. Thus, arTVs appear to encode different factors for counteracting TA systems.

## Viral metabolic genes

The pangenome of arTVs includes many predicted metabolic genes with specific or general functions (S5 Table), such as phosphoadenosine phosphosulphate reductase (HRTV-29 and HFTV1), which could be involved in sulfur assimilation pathway; Class II glutamine amidotransferase (phiCh1 and ChaoS9); MaoC-like dehydratase domain protein (HATV-2), which exhibits (R)-specific enoyl-CoA hydratase activity [58]; dual specificity protein phosphatase (HGTV-1 and viruses from the genus *Mincapvirus* of the *Hafunaviridae*), thioredoxin (HGTV-1), peroxide stress response protein YaaA (HRTV-29 and ChaoS9), ADP-ribosyl-transferase (HVTV-1 and HCTV-5), sialidase (HCTV-5) and Hsp90 chaperone protein (HSTV-2 and HGTV-1). However, the most common and widespread metabolic genes in arTVs appear to be linked to the replication, transcription, and translation of the viral genomes.

**Nucleotide biosynthesis.** Most haloarchaeal virus genes involved in DNA metabolism belong to the pyrimidine biosynthesis pathway and encode thymidylate synthase, thymidylate kinase, dCTP deaminase, and nucleoside deoxyribosyltransferase (Fig 4B). These proteins are broadly encoded by viruses from families *Hafunaviridae*, *Druskaviridae*, *Soleiviridae*, *Halo magnusviridae*, *Haloferuviridae*, and *Vertoviridae* (S5 Table). Interestingly, HGTV-1 (family *Halomagnusviridae*) encodes the small subunit, CarA, of carbamoylphosphate (CP) synthetase, which hydrolyzes glutamine to CP, a precursor common for the biosynthesis of pyrimidines and arginine [59] (Fig 4B). Thus, HGTV-1 might be redirecting the host metabolism toward de novo synthesis of pyrimidines. To our knowledge, no other viral isolate thus far has been reported to encode CarA/CarB subunits, although other enzymes in the de novo pyrimidine
synthesis pathway were found encoded by bacterial viruses [60]. Overall, the different genes identified suggest that haloarchaeal tailed viruses likely boost up the pyrimidine metabolism during infection. A similar behavior has been observed in *Pseudomonas* virus infections [61], suggesting that enhanced pyrimidine metabolism is critical for efficient replication of both archaeal and bacterial tailed viruses.

Viruses of the *Hafunaviridae*, *Druskaviridae*, and *Halomagnusviridae* encode class II ribo-nucleotide diphosphate reductases (RnRs), which convert rNTPs into corresponding dNTPs, the essential building blocks for the synthesis of viral DNA (Fig 4B, S5 Table). Viruses from the *Druskaviridae* encode putative 5'(3')-deoxyribonucleotidases, which catalyze the dephosphorylation of nucleoside monophosphates, likely to regulate and maintain the homeostasis of nucleotide and nucleoside pools in the host cells [62] (Fig 4B, S5 Table).

Cobalamin (vitamin B₁₂) is an important cofactor in various metabolic pathways, including DNA biosynthesis (e.g., for class II RnR) [63]. Viruses of the *Druskaviridae* and *Saparoviridae* encode putative cobaltochelatase subunits CobS and CobT [EC Number: 6.6.1.2], which show sequence similarity to the 2 subunits encoded by certain tailed bacteriophages and cellular organisms (S5 Table). CobS is an AAA+ ATPase, whereas CobT contains a von Willebrand factor type A (vWA) domain and a metal ion-dependent adhesion site (MIDAS) domain located in the carboxyl-terminal region [64]. In HHTV-2, *cobT* is split into 2 open reading frames (ORFs) with the vWA domain encoded by a separate ORF. In prokaryotes, CobS interacts with CobT and, together with the third subunit, CobN, forms an active cobaltochelatase complex [65]. The latter catalyzes the insertion of Co²⁺ ion into the corrin ring of hydrogenogenbyrinate a,c-diamide, close to the final step in the aerobic cobalamin biosynthesis pathway [63]. We found that the *cobST* 2-gene cluster is widely encoded in tailed viruses that infect members of 8 bacterial or archaeal orders, including Halobacteriales, as described in this study (S6 Table). In cyanophages, *cobST* gene cluster is part of the core genome [66], although *cobT* is usually mistakenly annotated as a peptidase. Recent bioinformatic study on the genotypes of *cob* subunits in prokaryotic genomes revealed that *cobS* and *cobT* are absent in most prokaryotic genomes, whereas *cobN* is present [65]. Many such prokaryotes, including haloarchaea, were suggested to employ a mosaic cobaltochelatase consisting of CobN and magnesium chelatase subunits [65,67]. Our finding of the broad distribution of the *cobST* gene cluster in tailed viruses, both bacterial and archaeal, suggests that viral CobST likely hijacks the host CobN (from the mosaic cobaltochelatase in the absence of the host encoded CobST), to increase the production of cobalamin, thereby promoting the virus replication.

**Transcription and translation.** It has been previously shown that arTVs encode many proteins involved in RNA metabolism [38], and our current study further expands the complement of viral genes implicated in RNA metabolism and protein translation. For instance, HJTV-2 (*Hafunaviridae*) encodes a Trm112 family protein, which in *Haloferax volcanii* interacts and activates 2 MTases, Trm9 and Mtq2, known to methylate tRNAs and release factors, respectively [68,69]. Remarkably, 4 viruses, HRTV-24, HSTV-3, HJTV-3, and HRTV-15 (*Hafunaviridae*), encode homologs of the ribosomal protein L21e (S5 Fig), which mediates the attachment of the 5S rRNA onto the large ribosomal subunit and stabilizes the orientation of adjacent RNA domains [70]. Ribosomal proteins have been recently shown to be encoded by diverse tailed bacteriophages and at least some of these proteins can be incorporated into host ribosomes [71]. The L21e homologs identified herein are the first examples of ribosomal proteins encoded by archaeal viruses.

*Halomagnusvirus* HGTV-1 encodes a record number of tRNA genes among archaeal viruses: 36 tRNA genes corresponding to all 20 proteinogenic amino acids [38]. In addition, this virus has a number of tRNA metabolism genes: 2 Rnl2-family RNA ligases, which may be involved in the removal of tRNA introns [38]; a tRNA nucleotidyl transferase (CCA-adding
enzyme) involved in tRNA maturation; a Class I lysyl-tRNA synthetase which catalyzes the attachment of lysine to its cognate tRNA; and a tRNA splicing ligase RtcB possibly responsible for the repair of viral tRNAs. RtcB and several tRNAs are also encoded by viruses from the *Hafunaviridae* and *Druskaviridae* (S1 Table).

**Integrases.** All members of *Hafunaviridae*, *Graaviviridae*, *Vertoviridae* and *Leisingerviridae* as well as HCTV-5 from *Druskaviridae* (but not other members of this family) encode site-specific tyrosine recombinases. Consistently, analysis of the archaeal genome sequences available in NCBI database revealed the presence of apparently functional (based on the conservation of all key proteins), integrases-encoding proviruses from the former 4 families, including several previously reported proviruses [20,22] (S7 Table, S6 Fig), suggesting a temperate life cycle for members of the *Hafunaviridae*, *Graaviviridae*, *Vertoviridae* and *Leisingerviridae*. All proviruses were flanked by recognizable direct repeats corresponding to the attachment sites (S7 Table). By contrast, the recombinase of HCTV-5 is homologous to the invertase encoded by viruses of *Vertoviridae*, which is responsible for the inversion of the tail-fiber module of phiCh1 [72]. Thus, HCTV-5 recombinase is likely to be involved in genome rearrangements rather than site-specific integration. Consistently, proviruses related to HCTV-5 or other members of the *Druskaviridae* were not identified. Searches for proviruses related to other arTV families revealed that members of the families *Haloferuviridae* and *Anaerodiviridae* can also have temperate members (S7 Table, S6 Fig), although none of the current isolates encodes an integrase, suggesting that the integration module and hence the integration ability can be occasionally gained by arTVs. In total, 38 proviruses belonging to 6 arTV families (based on the same demarcation criteria which we used to define the families) were identified in diverse species of halophilic and methanogenic archaea (S7 Table), considerably expanding the genetic diversity and distribution of the corresponding virus families. The majority of haloarchaeal proviruses were integrated into diverse tRNA genes, whereas proviruses belonging to the *Anaerodiviridae* were integrated either into intergenic regions or protein-coding genes (S7 Table). No proviruses related to members of the other 8 families could be identified.

## Mutations in tail fiber genes determine the broad host range of hafunaviruses

The arTVs from different families display highly variable host ranges, with some being specific to a single isolate and others infecting multiple species from different haloarchaeal genera [16,17,22,25]. We did not observe obvious relationship between the host range and virus taxonomy, with some representatives of the same family displaying a broad host range and others being able to infect only 1 strain. Thus, host range determinants appear to be virus specific. Unlike in the case of tailed bacteriophages, the factors underlying the host range in haloarchaeal viruses remain obscure. To gain insights into this question, we focused on viruses of the family *Hafunaviridae* which, with 39 myovirus isolates, is currently the largest family of arTVs. Furthermore, most hafunaviruses have broad host ranges, collectively being able to infect a dozen of different strains belonging to 5 haloarchaeal genera [16,17,22]. We assessed the efficiency of plating (EOP) of 24 hafunaviruses (19 isolates from 5 species in the genus *Haloferacalesvirus* and 5 isolates from a single species in the genus *Mincapvirus*) on a panel of 24 haloarchaeal strains from 5 genera [16,17]. The host ranges of viruses within each species varied considerably (S8 Table), consistent with the previous results [16,17].

Availability of complete genome sequences for multiple isolates from the same species allowed us to pinpoint the genes responsible for the observed differences in the host range. The most informative was the comparison of 2 isolates, HRTV-19 and HRTV-23, which differ by 2 nucleotide substitutions, but their EOPs on 2 of the tested strains differ by 3 orders of
<figure><img src="image_10.png" /><figcaption>**Fig 5. Host ranges of hafunaviruses.** The heat map shows the EOP of hafunavirus isolates on different haloarchaeal strains. The EOP of the virus on its original isolation host has been set to 1 and marked with H. Number 1 is equal to the EOP on the original host. Numbers -1, -2, ... refer to $10^{-1}$, $10^{-2}$, ... and +1 refers to $10^1$ when compared to the EOP on the original host, as shown in S8 Table. The isolation host for HRTV-26 is *Halorubrum* sp. SS13-13, which occasionally did not form a lawn and was not included in this analysis. The titer of HRTV-26 on *Halorubrum* sp. SS13-13 was $2 \times 10^8$ plaque-forming units per milliliter, which was considered here as 1 and used for comparison. The hosts were clustered based on the similarity of EOPs of the tested viruses, as represented by the dendrogram on the right side. The upper panel shows the maximum likelihood phylogenetic tree of the tail fiber adhesin proteins encoded by the analyzed viruses. The adhesin of HRTV-2, the amino acid sequence of which is markedly divergent from those of other hafunaviruses, is set as the outgroup. The 4 distinct groups of adhesins are displayed with colored blocks. Bootstrap values greater than 90% are indicated in the nodes by dots. Virus names are colored according to the species to which they belong. EOP, efficiency of plating.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001442.g005

magnitude (Fig 5, S8 Table). The 2 mutations mapped to 2 adjacent genes located at the end of the tail morphogenesis module, suggesting that this locus plays a key role in determining the host specificity. This possibility is further supported by the observation that 5 isolates (HCTV-7, HCTV-8, HCTV-9, HCTV-10, and HCTV-11) in which the 2 genes are identical but sequence divergence in other loci ranges from 1 bp to approximately 1.2 kbp share the same host ranges, with similar EOPs (Fig 5, S8 Table). The first of the 2 genes encodes a glycine-rich tail fiber protein and the second gene encodes a small putative protein (S7 Fig). The glycine-rich protein displays features typical of adhesin proteins located at the distal tip of the tail fibers of diverse T-even phages and has been shown to determine the host specificity in phages [73]. The adhesins encoded by hafunaviruses form 4 distinct phylogenetic groups (Groups 1 to 4), with the corresponding viruses displaying distinct infectivity patterns on the tested haloarchaeal strains (Fig 5, S8 Fig). This is most prominent in the case of viruses encoding the largest adhesin groups, Group 1 and Group 3. For instance, 4 strains, *Halorubrum* sp. SS8-7, *Halorubrum sodomense*, *Haloarcula californiae*, and *Halobellus* sp. SS6-7, are particularly sensitive to
viruses encoding Group 1 adhesins, with the latter 2 strains being exclusively infected by viruses encoding this group of adhesins. Conversely, no virus of this group is able to infect *Halorubrum* strains SS7-4, SS10-3, and SS5-4 or *Haloarcula* sp. SS8-5. Nearly all viruses encoding Group 3 adhesins are able to infect *Halorubrum* sp. SS9-12, whereas *Halorubrum* sp. SS6-1 and *Haloterrigena* sp. S13-7 are completely insensitive to this group of viruses (Fig 5) (see S1 Text for details). In all hafunaviruses, the adhesin gene is located within a hypervariable region (S7 Fig), which has been previously suggested to encode tail fiber proteins [22]. Collectively, our results implicate the adhesin-encoding gene as the key host range determinant in hafunaviruses. Importantly, while the host range pattern can be explained, even if partly, by the adhesin phylogeny, overall intergenomic relationships between hafunaviruses is a poor predictor of the host range. Indeed, viruses belonging to the same species (i.e., >95% genome-wide identity) encode adhesins from different phylogenetic groups and display distinct host ranges.

Continuity between archaeal tailed viruses across biomes

The previous study of viral communities using metagenomics approach revealed abundant presence of tailed viruses in hypersaline environments [74]. To assess the extent of sequence diversity and environmental distribution of halophilic arTVs and to explore the evolutionary relationships between tailed viruses from different environments, we searched the Integrated Microbial Genomes and Microbiomes (IMG/M) database for homologs of MCP, portal and the large terminase (TerL) proteins representing each family of arTVs. The searches collectively yielded 117,049 contigs encoding at least 1 of the 3 viral hallmark proteins. The dataset was supplemented with reference sequences from viruses for which hosts are known or predicted in silico, including tailed (pro)viruses associated with methanogens, MGI and MGII Euryarchaeota, Thaumarchaeota, Thermoplasmata, and Aigarchaeota [26, 27, 29, 31–35]. In addition, sequences of uncultivated haloviruses previously obtained through fosmid sequencing from the Santa Pola saltern, Spain, were also added to the dataset [75]. Among the 3 hallmark proteins, TerL appeared to be the least specific to archaeal viruses, likely due to its relatively high conservation across all tailed bacterial and archaeal viruses. Accordingly, searches queried with haloarchaeal virus TerL against the IMG/M database retrieved a high number of significant hits from diverse environments (S11 Fig). By contrast, MCP and portal proteins were more specific to halophilic arTVs and retrieved a higher ratio of homologs from hypersaline ecosystems compared to those from other habitats. This is consistent with the fact that haloviral MCPs and portal proteins typically display low sequence similarity to their bacteriophage homologs in BLASTP searches. Phylogenetic analysis of the arTV MCP and portal proteins produced largely congruent trees (S12 Fig), suggesting that the 2 proteins coevolve and the corresponding genes rarely undergo nonorthologous replacement. The only possible exceptions to this trend include HGTV-1 and ChaoS9, with the latter being notoriously chimeric [20] (S1G Fig). Consequently, MCP and portal proteins were used as specific markers for arTVs.

The global MCP tree splits into 9 well-supported clades, each represented by haloarchaeal tailed virus isolates from one or more families (Fig 6). Structural comparison of the MCPs from representatives of the 14 arTV families as well as selected uncultured arTVs revealed the conservation of the HK97 fold (Fig 7A) and recapitulated the same 9 clades obtained using sequence-based phylogenetic analysis (Fig 7, S1 Text), testifying to the robustness of these groups. Finally, a similar clustering is also observed in the portal tree (S11 and S12 Figs). In the global MCP phylogeny including sequences from metagenomic databases, almost all clades in the MCP tree contain numerous homologs originating from hypersaline environments (Fig 6),
<figure><img src="image_11.png" /><figcaption>**Fig 6. Phylogenetic analysis of the MCPs encoded by arTVs.** The homologs retrieved from metagenomes were clustered at 0.9 protein sequence identity and the representative of each cluster is indicated as a circle at the tip of the branch with the color reflecting the corresponding biome (see the legend on the left). The ring outside of the tree consists of tiles indicating the biome compositions from the corresponding protein cluster. The dot in the outermost ring indicates the (predicted) virus type in the corresponding branch and is colored based on the host predicted. The names of the arTV isolates of haloarchaea and methanogens as well as the fosmids generated from a saltern of Santa Pola in Spain [75] are shown on the tree. The virus family names are indicated next to the corresponding arTV isolates. The associated data on metagenomic contigs from the IMG/M database can be found in S2 Data. arTV, archaeal tailed virus; IMG/M, Integrated Microbial Genomes and Microbiomes; MCP, major capsid protein.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001442.g006

with the corresponding metagenomes derived from all 7 continents, testifying to the wide geographic distribution of the halophilic arTVs. This result suggests that the haloarchaeal virus isolates described herein adequately represent the overall diversity of haloarchaeal tailed virus communities.

Interestingly, although searches with the hallmark proteins were performed against the global IMG/M database, only a handful of sequences affiliated to bacteria/phages were retrieved (Fig 6, S10 Fig). By contrast, haloarchaeal viruses formed well-supported clades with the previously characterized viruses associated with marine and methanogenic archaea, suggesting that all arTVs are more closely related to each other than to tailed bacteriophages, consistent with the GRAViTy analysis (Fig 1). Indeed, 3 of the clades in the MCP tree contained considerable fractions of taxa originating from nonhypersaline environments (Fig 6). The largest of these clades is dominated by MCP homologs from marine environments, with a relatively small subclade including haloarchaeal viruses of the *Druskaviridae* (e.g., HVTV-1), interspersed between reference viruses associated with MGII Euryarchaeota (magroviruses) [34,35], marine Thaumarchaeota [32], and Thermoplasmata [33] (Fig 6). A similar clustering
<figure><img src="image_12.png" /><figcaption>**Fig 7.** Structural comparison of the arTV MCPs. (A) Structural models of arTV MCPs shown using ribbon representation and colored using rainbow scheme from the N-terminus (blue) to C-terminus (red). The models are shown next to the corresponding tips of the cladogram based on the pairwise structural similarity (Z) scores. (B) All-against-all comparison of the arTV MCP structural models shown as a matrix of pairwise Z scores generated using DALI [93]. The underlying data can be found in S2 Data. arTV, archaeal tailed virus; MagroA/B, group A and B Magroviruses; MCP, major capsid protein; Thaum, a virus of marine Thaumarchaeota; Thermopl, a virus of marine Thermoplasmata.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3001442.g007

is also retrieved in the portal protein tree (S11 and S12 Figs), suggesting that the ancestor of the *Druskaviridae* has evolved from the midst of marine archaeal viruses.

The second MCP tree clade containing sequences from diverse environments is represented by viruses from the *Graaviviridae* and *Vertoviridae* (e.g., BJ1 and phiCh1, respectively; note that the MCP of ChaoS9 from the *Vertoviridae* clusters with that of HHTV-1 from the *Madisaviridae*). This group of halovirus MCPs is embedded within an assemblage of homologs from diverse nonsaline environments, including anoxic habitats (Fig 6). Consistently, the clade includes viruses psiM2 and Drs3 (families *Leisingerviridae* and *Anaerodiviridae*, respectively) infecting methanogenic archaea (Fig 6). The same relationship is recapitulated in the portal protein tree (S11 and S12 Figs). Therefore, viruses from the families *Graaviviridae* and *Vertoviridae* appear to have evolved from viruses inhabiting moderate, nonsaline environments, although the exact evolutionary history remains obscure, given the complicated composition and relatively poor sampling of this clade.

Finally, the third MCP tree clade is represented by viruses from the families *Haloferuviridae* and *Halomagnusviridae* (e.g., HFTV1 and HGTV-1, respectively). Unlike in the other 2 cases, where haloarchaeal viruses were nested among taxa from nonhypersaline environments, in this clade, we observe a reverse situation. Namely, a group of sequences, largely originating from anoxic environments, such as methanogenic microbial communities, anaerobic digesters, wetlands, etc., with at least 1 contig predicted to represent a virus infecting a methanogenic archaeon, forms a sister clade to haloviruses of the *Haloferuviridae* and is further embedded within a clade dominated by sequences from hypersaline environments (Fig 6). A parsimonious explanation for such clustering involves evolution of a group of methanogenic archaeal viruses from within the diversity of haloarchaeal viruses.
Collectively, these results suggest that there is a continuity between arTVs across different biomes. At the same time, however, it becomes clear that haloarchaeal viruses represent a polyphyletic assemblage, with several groups immigrating into hypersaline environments from other habitats with lower salinity. Presumably, the ancestors of these groups gained the ability to infect halophilic archaea by accumulating mutations within the tail fiber genes, akin to those described above for hafunaviruses. We note that natural mixing between marine and hypersaline ecosystems (e.g., through evaporation of coastal marine water) provides an ecological setting for marine viruses to encounter halophilic archaea. Thus, hypersaline environments, once considered as an isolated extreme environment [76,77], emerge as an integral part of the whole ecosystem.

## Concluding remarks

In the present study, we sequenced 37 haloarchaeal tailed virus isolates, which originate from 4 spatially remote locations, more than doubling the current number of sequenced arTVs (S1 Table). Unlike for some hyperthermophilic archaeal viruses [78], for most arTVs, there is no relationship between genetic closeness (and hence taxonomy) and the site of virus isolation. For 6 out of the 7 virus families with more than 1 isolate (*Hafunaviridae, Druskaviridae, Haloferuviridae, Graaviviridae, Vertoviridae, and Leisingerviridae*), members were isolated from 2 to 5 geographically remote locations (S1 Table). In the cases of families *Hafunaviridae* and *Druskaviridae*, which contain members belonging to the same species (i.e., >95% identical at the genome level), nearly identical viruses were isolated from remote sites, indicating global distribution of certain arTV species. Through a systematic analysis of all 63 available arTVs, we establish the first framework for the taxonomy of arTVs and propose 14 new families and 3 orders for their classification to be included within the class *Caudoviricetes* alongside tailed bacteriophages. Half of the arTV families are currently represented by a single virus isolate, indicating the scarce sampling of haloarchaeal virome. While this manuscript was under review, 3 new haloarchaeal viruses and one virus infecting methanogenic archaea have been described [22,79,80]. *Halorubrum* viruses Serpecor1 and Hardycor2 belong to the family *Hafunaviridae* [22], whereas *Halorubrum* virus Hardycor1 [79] and *Methanocaldococcus* virus MFTV1 [80] are only distantly related to other arTVs and will likely represent new virus families.

It has been previously shown that arTVs share virion architecture, assembly principles, and genome organization with their bacterial virus counterparts [11,13,30,38,81]. In this study, we focused on metabolic genes and showed that arTVs carry a rich repertoire of counterdefense and metabolic genes, predominantly those involved in DNA and RNA metabolism. Many of these genes are also commonly encoded by tailed bacteriophages, suggesting common propagation strategies for these bacterial and archaeal viruses. Nevertheless, some of the identified metabolic genes have never been reported in other viruses, e.g., *carA*, which is predicted to promote the de novo synthesis of pyrimidines, RNA MTase activator *trm112* implicated in modulation of protein translation, and the gene encoding a putative ribosomal protein L21e.

Myoviruses of the *Hafunaviridae*, the most populous family of arTVs, have broad host ranges, with strains of the same species displaying considerable variation in the host range. We found that the host range specificity in these viruses is primarily determined by mutations within a gene encoding a tail adhesin protein, which resembles that of T-even bacteriophages [73,82]. High sequence divergence of adhesins encoded by hafunaviruses likely allows these viruses to explore a broad landscape of receptor molecules on the host surface, boosting the competitiveness of this group of viruses in the environments with some of the highest reported virus-to-host ratios. Despite the fundamental differences of archaeal and bacterial cell surface
structures [83], the host range determinant of tailed viruses, the tail adhesins, can be highly similar in sequence features between arTVs and tailed bacteriophages. This highlights the evolutionary potential of the tail fibers adapted to diverse host receptors, which may underlie the overall success of the Caudoviricetes viruses.

Finally, the survey of metagenomic databases uncovered a considerable diversity and global distribution of arTVs. Our results point to the polyphyletic origins of the haloarchaeal tailed viruses, with some of the groups originating from archaeal viruses associated with marine and methanogenic hosts, indicative of virus movement across different biomes. Given that in hypersaline environments the arriving viruses are exposed to extremely high salt concentrations both inside and outside of the cell [84], the adaptation must have entailed a substantial evolution of the virion proteins. Indeed, some of the halophilic arTVs have been shown to be noninfectious under low salt conditions [11,37]. However, the effect was reversible and upon addition of salt, the infectivity was restored, suggesting that virions of haloarchaeal viruses can maintain stability and integrity under a wide range of environmental conditions. Further sampling and isolation of arTVs, especially those infecting nonhalophilic archaea, will be important for further understanding the ecology and evolution of this important group of viruses.

# Materials and methods

## Preparation of virus samples

Viruses sequenced in this study had been isolated previously [16,17]. The host strains used for virus production are listed in S1 Table. Viruses and strains were grown aerobically at 37°C in modified growth medium (MGM) made using artificial salt water (SW) (23% (w/v) in broth, 20% in solid, and 18% in top agar) (http://www.haloarchaea.com/resources/halohandbook/). The double-layer plaque assay method was used for virus growth and quantification, and virus stocks were prepared from semiconfluent plates as previously described [16,17]. Viruses were precipitated from stocks with polyethylene glycol 6000 (final concentration of 10% (w/v)) by magnetic stirring for 1 hour at 4°C, collected by centrifugation (Fiberlite F14 rotor, 9,820 g, 40 minutes, 4°C), and resuspended in 18% SW. Virus samples were purified by rate zonal centrifugation in 5% to 20% (w/v) sucrose (18% SW; Sorvall AH629 rotor, 103,586 g, 40 minutes, 15°C), concentrated by differential centrifugation (Sorvall T647.5 rotor, 113,580 g, 3 hours, 15°C) and resuspended in 18% SW.

## Viral genome isolation and sequencing

Nucleic acids were extracted from purified virus samples either using the PureLink Viral RNA/DNA Mini Kit (Thermo Fisher Scientific, Germany) or by phenol/ether method combined with ethanol precipitation. The viral genome libraries were prepared with the Nextflex PCRFree kit (Bioo Scientific, Austin, Texas) and sequenced by Illumina Miseq (Illumina, San Diego, California, USA) with paired-end 250-bp read length (Biomics Platform, Institut Pasteur, France). Sequenced reads were quality-trimmed using fqCleaner v0.5.0 and assembled using clc_assembler v4.4.2 implemented in Galaxy-Institut Pasteur. All genome sequences were submitted to GenBank, and the accession numbers are provided in S1 Table.

## Viral genome annotation and bioinformatics analysis

The virus genomic termini and packaging mechanism were determined by PhageTerm [85]. ORFs were predicted using RAST v2.0 server [86]. Functional annotation of viral proteins was performed using HHpred [51]. The average nucleotide identity between viruses was calculated by ANI Calculator [87], whereas the global alignment of amino acid sequences was carried out
by EMBOSS Needle tool [88]. Multiple sequence alignments were performed by PROMALS3D [89]. Phylogenetic analyses were carried out using PhyML 3.0, with the automatic model selection [90]. To identify the proviruses related to arTVs, genome sequences of representatives of each proposed genus were used as seeds for blastx searches against the nonredundant protein sequence database at NCBI with the E-value cutoff of $1 \times 10^{-5}$. The exact nucleotide coordinates and attachment sites were identified manually.

## Structural modeling and comparisons

Structural modeling for all but one arTV MCPs was performed using AlphaFold2, with the mmseq2 algorithm for collection and alignment of homologs [91]. In the case of HHTV-1 MCP, AlphaFold2 failed to produce a satisfactory model and hence RoseTTAFold was used instead [92]. The resulting structural models were manually inspected and, when present, the N-terminal extensions corresponding to the scaffolding domains were removed. All structural models were then compared to each other using DALI [93] to produce a matrix and a dendrogram of structural similarity.

## Taxonomy assignment

The viral genomes were subjected to GRAViTy [47,48] (http://gravity.cvr.gla.ac.uk) and vConTACT v2.0 [49] for taxonomy assignment. We incorporated the annotated genomes of arTVs into the genomic datasets that were (i) built for bacterial and archaeal dsDNA viruses from the ICTV 2016 Master Species List 31V1.1 for GRAViTy analysis; and (ii) embedded in vConTACT v2.0 with the NCBI Bacterial and Archaeal Viral RefSeq V88 and default settings for analysis. The virus network generated by vConTACT v2.0 was visualized with Cytoscape software v.3.7.2 (https://cytoscape.org/) using an edge-weighted spring-embedded model that genomes sharing more PCs were positioned closer to each other. Pair-wise genomic comparison was performed by EasyFig with the E-value cutoff of $1 \times 10^{-3}$ and minimum protein length of 15 residues [94]. The overall similarity among arTVs on genus and family levels were analyzed and proteins with over 30% amino acid sequence identity and E-value < $1 \times 10^{-25}$ in our arTV database were counted as homologous. The results were visualized with Circos plot [95], whereas box plots were prepared using OriginLab v9.8 (OriginLab, USA).

## Virus-host interactions

Virus-host pairs tested for hafunaviruses in this study are shown in S8 Table. For initial screening of virus-host interactions, 10-µl drops of undiluted and 100-fold diluted virus stocks were spotted on the lawn of early-stationary growing strains prepared by the double-layer method, and 10-µl drops of MGM broth were spotted as negative control. After 2 to 5 days of incubation, plates were screened for the growth inhibition. When inhibition was observed, plaque assay was used to verify virus-host interactions and quantify the EOP. The titers above $10^3$ plaque-forming units per milliliter could be detected. Data were normalized by comparing the EOP to that on the original host strain (set as 1). The heat map of the virus host ranges was generated based on the EOP results and clustered according to the phylogeny of viral adhesins as well as EOP similarity of the tested hosts using plothmap package in R.

## Metagenomic database screening

The MCP, portal, and TerL protein sequences of the isolated halophilic arTVs were used as queries in a BLAST search of 31,344 public metagenomes available in the IMG/M database at JGI with a threshold of 100 bit score and 80% sequence coverage. Sequences of the 3 markers
from available (pro)viruses of methanogen, marine group I and group II Euryarchaeota, Thau-
marchaeota, Thermoplasmata, Aigarchaeota, and fosmids from a saltern of Santa Pola in
Spain were extracted [26,27,29,31–35,72]. These sequences were subjected to BLAST searches
against our retrieved metagenomic sequences with the same threshold as above and were used
as reference sequences in the phylogenetic analysis. Of the 117,066 retrieved sequences, 20,720
were use-restricted and were not considered further. The 96,345 remaining sequences were
clustered at 90% sequence identity separately and combined into the final datasets. Sequences
in the datasets of the 3 hallmark proteins were aligned using MAFFT v7 [96], followed by
removal of poorly aligned positions by trimAl with gap threshold of 0.2 [97]. The phylogenetic
trees were constructed using a JTT+CAT model by FastTree (version 2.1.11) [98], and visual-
ized using ggtree v2.4.1 [99]. Biome information for IMG/M metagenomes were obtained
from the Gold database [100]. All IMG/M contigs with a hit to an arTV marker were also
cross-referenced with contigs in the IMG/VR database [101]. For contigs detected as viral and
included in IMG/VR, the host prediction was obtained when available at the domain level (i.e.,
bacteria versus archaea) and displayed on the trees.

**Supporting information**

**S1 Text. Supporting information text.**

(PDF)

**S1 Data. The dendrogram of CGJ distances for classified bacterial and archaeal dsDNA viruses.** CGJ, composite generalized Jaccard; dsDNA, double-stranded DNA.

(NWK)

**S2 Data. Raw data used to make graphs presented in this article.**

(XLSX)

**S1 Table. Properties of arTVs. arTV, archaeal tailed virus.**

(XLSX)

**S2 Table. PCs generated for arTVs by vConTACT v2.0.** arTV, archaeal tailed virus; PC, protein cluster.

(XLSX)

**S3 Table. Viral genome annotations.**

(XLSX)

**S4 Table. Virus-encoded DNA MTases.** MTase, methyltransferase.

(XLSX)

**S5 Table. Virus-encoded metabolic genes.**

(XLSX)

**S6 Table. Cobaltochelatase subunits gene cluster encoded by head-tailed viruses.**

(XLSX)

**S7 Table. Proviruses related to arTVs in archaeal genomes.** arTV, archaeal tailed virus.

(XLSX)

**S8 Table. Host ranges of haloarchaeal head-tailed viruses.**

(XLSX)

**S1 Fig. Genome alignment of arTVs.** Comparisons of genomes of (A) representative viruses from the 4 genera in *Hafunaviridae*, (B) viruses from the 2 genera in *Druskaviridae* and viruses
from the 2 genera in *Saparoviridae*, (C) HATV-2 from *Soleiviridae* and HGTV-1 from *Halo-magnusviridae* (singletons), (D) the 3 viruses from *Haloferuviridae* (3 genera), (E) HATV-3 from *Pyrstoviridae*, HFTV1 from *Haloferuviridae*, HSTV-1 from *Shortaselviridae* and HRTV-28 from *Suolaviridae* (all singletons), (F) the 2 members from *Graaviviridae*, (G) the 3 viruses from *Vertoviridae* (2 genera), (H) psiM2 representing *Leisingerviridae* and Drs3 from *Anaero-diviridae*. Putative protein functions are indicated above or below the corresponding ORFs.
Genes encoding virus morphogenesis-related proteins are colored in green, whereas replication-related genes are colored in red. Homologous genes shared between viruses are connected by shadings of different degrees of gray based on the amino acid sequence identity. See S1 Table for more information on families and genera. arTV, archaeal tailed virus.
(PDF)

**S2 Fig. The box plot shows the percentage of genes shared by arTVs.** The percentage of genes of a representative virus from each genus that shared with members of the proposed genus (G) and family (F), as well as arTVs from other families (A) are shown. Each box represents the middle 50th percentile of the data set and is derived using the lower and upper quartile values. The median value is displayed by a horizontal line. Whiskers represent the maximum and minimum values with the range of 1.5 IQR. Each virus is represented by dots. Proteins with over 30% amino acid sequence identity and E-value < 1 × 10⁻²⁵ in our arTV database are counted as homologous proteins. Data underlying this figure can be found in S2 Data. arTV, archaeal tailed virus.
(PDF)

**S3 Fig. Multiple sequence alignments of different categories of MTases encoded by arTVs.** The boxed signature motifs are used to classify the MTases. arTV, archaeal tailed virus; MTase, methyltransferase.
(PDF)

**S4 Fig. Multiple sequence alignment of micrococcal nucleases encoded by arTVs.** The experimentally verified micrococcal nuclease encoded by *Staphylococcus hyicus* (AAA26661.1) is used as a sequence reference in the alignment. The empty and filled stars indicate the calcium binding and catalytic sites, respectively. arTV, archaeal tailed virus.
(PDF)

**S5 Fig. Multiple sequence alignment of putative L21e homologs encoded by arTVs from the family Hafunaviridae.** The L21e encoded by *Halorubrum californiensi*, *Natrinema pallidum*, and *Haloarcula marismortui* are used as references in the alignment. arTV, archaeal tailed virus.
(PDF)

**S6 Fig. Genome comparisons of the representative proviruses with related arTVs from 6 viral families.** Genes encoding virus morphogenesis proteins, genome replication proteins and integrases are colored in green, red and orange, respectively. For visualization purposes, genomes of some proviruses were circularized and reopened at a different position, matching the start site of the reference arTV genome. Homologous genes shared between (pro)viruses are connected by shadings of different degrees of gray based on the amino acid sequence identity. See S7 Table for complete lists of proviruses of *Hafunaviridae*, *Haloferuviridae*, *Graaviviridae*, *Vertoviridae*, *Leisingerviridae*, and *Anaerodiviridae*. arTV, archaeal tailed virus; *att*, attachment site.
(PDF)
**S7 Fig.** The pairwise genome comparisons of viruses from the genera *Haloferacalesvirus* and *Mincapvirus* of the family *Hafunaviridae*. Putative protein functions are indicated above the corresponding ORFs. Genes encoding virus morphogenesis related proteins are colored in green, whereas replication-related genes are in red. The putative tail fiber adhesin-encoding gene is shown in yellow. The genes encoding restriction modification enzymes are in blue. Homologous genes shared between viruses are connected by shadings of different degrees of gray based on the amino acid sequence identity.
(PDF)

**S8 Fig.** The divergence of adhesin proteins encoded by viruses from the family *Hafunaviridae*. (A) Pairwise sequence alignments of adhesins encoded by viruses from the genera *Haloferacalesvirus* and *Mincapvirus*. The HVSs, which are predicted to interact with host receptor, are indicated. (B) Phylogeny of the adhesins. Viruses that were used for EOP assay in this study are indicated in black, otherwise in gray. Nodes with bootstrap support values greater than 90% are indicated with dots. EOP, efficiency of plating; HVS, hypervariable segment.
(PDF)

**S9 Fig.** The effect of mutations in adhesin to the host range changes in viruses of the 4 adhesin groups (see S8 Fig). Pairwise sequence alignments were produced for each of the 4 adhesin groups. The HVSs identified in S8 Fig are indicated in each group. The ratio of shared numbers of sensitive hosts/the total numbers of sensitive hosts of each 2 tested viruses was calculated (see S2 Data), and the heat map of hs versus adhesin pi was generated for each adhesin group. hs, host range similarity; HVS, hypervariable segment; pi, protein identity.
(PDF)

**S10 Fig.** Phylogenetic tree of VP1-like proteins encoded by viruses from the *Hafunaviridae*. The same colors as in S8B Fig were used to indicate the Groups of adhesin proteins encoded by the corresponding viruses. Nodes with the bootstrap support values greater than 90% are indicated with dots.
(PDF)

**S11 Fig.** Phylogenetic analysis of the portal and TerL proteins encoded by arTVs. The homologs retrieved from the metagenomics databases were clustered at 90% identity, and the representative of each cluster is indicated as a dot in the branch with the color denoting the source of its biome. The circle outside of the tree consisting of tiles indicates the biome compositions from the corresponding PC, which is represented by the dot in the branch. The dot in the outermost ring indicates the (predicted) virus type (colored according to its putative host) in the corresponding branch. The names of the tailed virus isolates of haloarchaea and methanogens, as well as the fosmids generated from a saltern of Santa Pola in Spain, are shown on the tree. arTV, archaeal tailed virus; PC, protein cluster.
(PDF)

**S12 Fig.** Maximum likelihood phylogenetic trees of the arTVs MCPs and portal proteins. The best-fitting models identified for the MCP and portal trees were selected by PhyML and were Blosum62+G+I+F and LG+G+I+F, respectively. Same colors are used for clades containing the same group of viral members in both trees. HGTV-1 and related haloarchaeal virus fosmids, as well as ChaoS9, that clustered differently in the 2 trees are indicated by lines. Bootstrap values are shown next to the nodes, with values above 80 shown in black and those below 80 in gray. arTV, archaeal tailed virus; MCP, major capsid protein.
(PDF)
**Acknowledgments**

We thank Sari Korhonen, Matti Ylänne, and Carlos Lapedriza for skillful technical assistance. We are also grateful to Marc Monot (Biomics Platform, Institut Pasteur) for help with the sequencing of viral genomes. The facilities and expertise of the HiLIFE Biocomplex unit at the University of Helsinki, a member of Instruct-ERIC Centre Finland, FINStruct, and Biocenter Finland are gratefully acknowledged.

**Author Contributions**

**Conceptualization:** Ying Liu, Tatiana A. Demina, Hanna M. Oksanen, Mart Krupovic.

**Formal analysis:** Ying Liu, Tatiana A. Demina, Simon Roux, Pakorn Aiewsakun, Darius Kazlauskas.

**Funding acquisition:** David Prangishvili, Mart Krupovic.

**Methodology:** Pakorn Aiewsakun, Peter Simmonds.

**Resources:** Hanna M. Oksanen.

**Supervision:** Hanna M. Oksanen, Mart Krupovic.

**Visualization:** Ying Liu, Simon Roux.

**Writing – original draft:** Ying Liu, Mart Krupovic.

**Writing – review & editing:** Ying Liu, Tatiana A. Demina, Simon Roux, Pakorn Aiewsakun, Darius Kazlauskas, Peter Simmonds, David Prangishvili, Hanna M. Oksanen, Mart Krupovic.

**References**

1. Koonin EV, Dolja VV, Krupovic M, Varsani A, Wolf YI, Yutin N, et al. Global Organization and Proposed Megataonomy of the Virus World. Microbiol Mol Biol Rev. 2020; 84(2):e00061-19. Epub 2020/03/07. https://doi.org/10.1128/MMBR.00061-19 PMID: 32132243.

2. Dion MB, Oechslin F, Moineau S. Phage diversity, genomics and phylogeny. Nat Rev Microbiol 2020; 18(3):125-38. Epub 2020/02/06. https://doi.org/10.1038/s41579-019-0311-5 PMID: 32015529;

3. Krupovic M, Dolja VV, Koonin EV. The LUCA and its complex virome. Nat Rev Microbiol 2020; 18 (11):661-70. Epub 2020/07/16. https://doi.org/10.1038/s41579-020-0408-x PMID: 32665595;

4. Dedeo CL, Cingolani G, Teschke CM. Portal Protein: The Orchestrator of Capsid Assembly for the dsDNA Tailed Bacteriophages and Herpesviruses. Annu Rev Virol. 2019; 6(1):141-60. Epub 2019/07/25. https://doi.org/10.1146/annurev-virology-092818-015819 PMID: 31337287.

5. Nobrega FL, Vlot M, de Jonge PA, Dreesens LL, Beaumont HJE, Lavigne R, et al. Targeting mechanisms of tailed bacteriophages. Nat Rev Microbiol 2018; 16(12):760-73. Epub 2018/08/15. https://doi.org/10.1038/s41579-018-0070-8 PMID: 30104690

6. Casjens SR, Hendrix RW. Bacteriophage lambda: Early pioneer and still relevant. Virology. 2015; 479-480:310-30. Epub 2015/03/07. https://doi.org/10.1016/j.virol.2015.02.010 PMID: 25742714.

7.Breitbart M, Bonnain C, Malki K, Sawaya NA. Phage puppet masters of the marine microbial realm. Nat Microbiol 2018; 3(7):754-66. Epub 2018/06/06. https://doi.org/10.1038/s41564-018-0166-y PMID: 29867096

8. Chow CE, Suttle CA. Biogeography of Viruses in the Sea. Annu Rev Virol 2015; 2(1):41-66. Epub 2016/03/10. https://doi.org/10.1146/annurev-virology-031413-085540 PMID: 26958906

9. Danovaro R, Corinaldesi C, Dell'anno A, Fuhrman JA, Middelburg JJ, Noble RT, et al. Marine viruses and global climate change. FEMS Microbiol Rev 2011; 35(6):993-1034. Epub 2011/01/06. https://doi.org/10.1111/j.1574-6976.2010.00258.x PMID: 21204862

10. Duda RL, Teschke CM. The amazing HK97 fold: versatile results of modest differences. Curr Opin Virol. 2019; 36:9-16. Epub 2019/03/12. https://doi.org/10.1016/j.coviro.2019.02.001 PMID: 30856581.
11. Pietilä MK, Laurinmäki P, Russell DA, Ko CC, Jacobs-Sera D, Hendrix RW, et al. Structure of the archaeal head-tailed virus HSTV-1 completes the HK97 fold story. Proc Natl Acad Sci U S A. 2013; 110(26):10604–9. Epub 2013/06/05. https://doi.org/10.1073/pnas.1303047110 PMID: 23733949.

12. Prangishvili D, Bamford DH, Forterre P, Iranzo J, Koonin EV, Krupovic M. The enigmatic archaeal viro-sphere. Nat Rev Microbiol 2017; 15(12):724–39. Epub 2017/11/11. https://doi.org/10.1038/nrmicro. 2017.125 PMID: 29123227

13. Baquero DP, Liu Y, Wang F, Egelman EH, Prangishvili D, Krupovic M. Structure and assembly of archaeal viruses. Adv Virus Res. 2020; 108:127–64. https://doi.org/10.1016/bs.aivir.2020.09.004 PMID: 33837715

14. Pietilä MK, Demina TA, Atanasova NS, Oksanen HM, Bamford DH. Archaeal viruses and bacterio-phages: comparisons and contrasts. Trends Microbiol 2014; 22(6):334–44. Epub 2014/03/22. https:// doi.org/10.1016/j.tim.2014.02.007 PMID: 24647075

15. Atanasova NS, Oksanen HM, Bamford DH. Haloviruses of archaea, bacteria, and eukaryotes. Curr Opin Microbiol 2015; 25:40–8. Epub 2015/05/02. https://doi.org/10.1016/j.mib.2015.04.001 PMID: 25932531

16. Atanasova NS, Demina TA, Buivydas A, Bamford DH, Oksanen HM. Archaeal viruses multiply: tempo-ral screening in a solar saltern. Viruses. 2015; 7(4):1902–26. Epub 2015/04/14. https://doi.org/10. 3390/v7041902 PMID: 25866903.

17. Atanasova NS, Roine E, Oren A, Bamford DH, Oksanen HM. Global network of specific virus-host interactions in hypersaline environments. Environ Microbiol 2012; 14(2):426–40. Epub 2011/10/19. https://doi.org/10.1111/j.1462-2920.2011.02603.x PMID: 22003883

18. Pfister P, Wasserfallen A, Stettler R, Leisinger T. Molecular analysis of Methanobacterium phage psiM2. Mol Microbiol 1998; 30(2):233–44. Epub 1998/10/29. https://doi.org/10.1046/j.1365-2958. 1998.01073.x PMID: 9791169

19. Wolf S, Fischer MA, Kupczok A, Reetz J, Kern T, Schmitz RA, et al. Characterization of the lytic archaeal virus Drs3 infecting Methanobacterium formicicum. Arch Virol 2019; 164(3):667–74. Epub 2018/12/14. https://doi.org/10.1007/s00705-018-04120-w PMID: 30523430

20. Dyall-Smith M, Palm P, Wanner G, Witte A, Oesterhelt D, Pfeiffer F. Halobacterium salinarum virus ChaoS9, a Novel Halovirus Related to PhiH1 and PhiCh1. Genes (Basel). 2019; 10(3). Epub 2019/03/ 06. https://doi.org/10.3390/genes10030194 PMID: 30832293.

21. Dyall-Smith M, Pfeifer F, Witte A, Oesterhelt D, Pfeiffer F. Complete genome sequence of the model halovirus PhiH1 (PhiH1). Genes (Basel). 2018; 9(10). Epub 2018/10/17. https://doi.org/10.3390/ genes9100493 PMID: 30322017.

22. Dyall-Smith M, Tang SL, Russ B, Chiang PW, Pfeiffer F. Comparative genomics of two new HF1-like haloviruses. Genes (Basel). 2020; 11(4):405. Epub 2020/04/12. https://doi.org/10.3390/ genes11040405 PMID: 32276506.

23. Klein R, Baranyi U, Rossler N, Greineder B, Scholz H, Witte A. Natrialba magadii virus phiCh1: first complete nucleotide sequence and functional organization of a virus infecting a haloalkaliphilic archaeon. Mol Microbiol 2002; 45(3):851–63. Epub 2002/07/26. https://doi.org/10.1046/j.1365-2958. 2002.03064.x PMID: 12139629

24. Pagaling E, Haigh RD, Grant WD, Cowan DA, Jones BE, Ma Y, et al. Sequence analysis of an Archaeal virus isolated from a hypersaline lake in Inner Mongolia, China. BMC Genomics. 2007; 8:410. Epub 2007/11/13. https://doi.org/10.1186/1471-2164-8-410 PMID: 17996081.

25. Mizuno CM, Prajapati B, Lucas-Staat S, Sime-Ngando T, Forterre P, Bamford DH, et al. Novel haloarchaeal viruses from Lake Retba infecting Haloferax and Halorubrum species. Environ Microbiol 2019; 21(6):2129–47. Epub 2019/03/29. https://doi.org/10.1111/1462-2920.14604 PMID: 30920125

26. Krupovic M, Makarova KS, Wolf YI, Medvedeva S, Prangishvili D, Forterre P, et al. Integrated mobile genetic elements in Thaumarchaeota. Environ Microbiol. 2019; 21(6):2056–78. Epub 2019/02/19. https://doi.org/10.1111/1462-2920.14564 PMID: 30773816.

27. Abby SS, Melcher M, Kerou M, Krupovic M, Stieglmeier M, Rossel C, et al. Candidatus Nitrosocaldus cavасurensis, an Ammonia Oxidizing, Extremely Thermophilic Archaeon with a Highly Mobile Genome. Front Microbiol. 2018; 9:28. Epub 2018/02/13. https://doi.org/10.3389/fmicb.2018.00028 PMID: 29434576.

28. Chow CE, Winget DM, White RA III, Hallam SJ, Suttle CA. Combining genomic sequencing methods to explore viral diversity and reveal potential virus-host interactions. Front Microbiol. 2015; 6:265. Epub 2015/04/29. https://doi.org/10.3389/fmicb.2015.00265 PMID: 25914678.

29. Krupovic M, Spang A, Gribaldo S, Forterre P, Schleper C. A thaumarchaeal provirus testifies for an ancient association of tailed viruses with archaea. Biochem Soc Trans 2011; 39(1):82–8. Epub 2011/ 01/27. https://doi.org/10.1042/BST0390082 PMID: 21265751
30. Krupovic M, Forterre P, Bamford DH. Comparative analysis of the mosaic genomes of tailed archaeal viruses and proviruses suggests common themes for virion architecture and assembly with tailed viruses of bacteria. J Mol Biol 2010; 397(1):144-60. Epub 2010/01/30. https://doi.org/10.1016/j.jmb.2010.01.037 PMID: 20109464

31. López-Pérez M, Haro-Moreno JM, de la Torre JR, Rodriguez-Valera F. Novel Caudovirales associated with Marine Group I Thaumarchaeota assembled from metagenomes. Environ Microbiol 2019; 21 (6):1980-8. Epub 2018/10/30. https://doi.org/10.1111/1462-2920.14462 PMID: 30370610

32. Ahlgren NA, Fuchsman CA, Rocab G, Fuhrman JA. Discovery of several novel, widespread, and eco-
logically distinct marine Thaumarchaeota viruses that encode amoC nitrification genes. ISME J. 2019;
13(3):618-31. Epub 2018/10/14. https://doi.org/10.1038/s41396-018-0289-4 PMID: 30315316.

33. Vik DR, Roux S, Brum JR, Bolduc B, Emerson JB, Padilla CC, et al. Putative archaeal viruses from the mesopelagic ocean. PeerJ. 2017; 5:e3428. Epub 2017/06/21. https://doi.org/10.7717/peerj.3428 PMID: 28630803.

34. Philosof A, Yutin N, Flores-Uribe J, Sharon I, Koonin EV, Beja O. Novel Abundant Oceanic Viruses of Uncultured Marine Group II Euryarchaeota. Curr Biol. 2017; 27(9):1362-8. Epub 2017/05/02. https://doi.org/10.1016/j.cub.2017.03.052 PMID: 28457865.

35. Nishimura Y, Watai H, Honda T, Mihara T, Omae K, Roux S, et al. Environmental Viral Genomes Shed New Light on Virus-Host Interactions in the Ocean. mSphere. 2017; 2(2). Epub 2017/03/07. https://doi.org/10.1128/mSphere.00359-16 PMID: 28261669.

36. Danovaro R, Dell'Anno A, Corinaldesi C, Rastelli E, Cavicchioli R, Krupovic M, et al. Virus-mediated archaeal hecatomb in the deep seafloor. Sci Adv. 2016; 2(10):e1600492. Epub 2016/10/21. https://doi.org/10.1126/sciadv.1600492 PMID: 27757416.

37. Pietilä MK, Laurinmäki P, Russell DA, Ko CC, Jacobs-Sera D, Butcher SJ, et al. Insights into head-tailed viruses infecting extremely halophilic archaea. J Virol. 2013; 87(6):3248-60. Epub 2013/01/04. https://doi.org/10.1128/JVI.03397-12 PMID: 23283946.

38. Senčilo A, Jacobs-Sera D, Russell DA, Ko CC, Bowman CA, Atanasova NS, et al. Snapshot of haloarchaeal tailed virus genomes. RNA Biol. 2013; 10(5):803-16. Epub 2013/03/09. https://doi.org/10.4161/rna.24045 PMID: 23470522.

39. Camarillo-Guerrero LF, Almeida A, Rangel-Pineros G, Finn RD, Lawley TD. Massive expansion of human gut bacteriophage diversity. Cell 2021; 184(4):1098-109.e9. Epub 2021/02/20. https://doi.org/10.1016/j.cell.2021.01.029 PMID: 33606979

40. Gregory AC, Zablocki O, Zayed AA, Howell A, Bolduc B, Sullivan MB. The Gut Virome Database Reveals Age-Dependent Patterns of Virome Diversity in the Human Gut. Cell Host Microbe. 2020; 28 (5):724-40.e8. Epub 2020/08/26. https://doi.org/10.1016/j.chom.2020.08.003 PMID: 32841606.

41. Wolf YI, Silas S, Wang Y, Wu S, Bocek M, Kazlauskas D, et al. Doubling of the known set of RNA viruses by metagenomic analysis of an aquatic virome. Nat Microbiol. 2020; 5(10):1262-70. Epub 2020/07/22. https://doi.org/10.1038/s41564-020-0755-4 PMID: 32690954.

42. Gregory AC, Zayed AA, Conceicao-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-23.e14. Epub 2019/04/30. https://doi.org/10.1016/j.cell.2019.03.040 PMID: 31031001.

43. Roux S, Brum JR, Dutilh BE, Sunagawa S, Duhaime MB, Loy A, et al. Ecogenomics and potential bio-
geochemical impacts of globally abundant ocean viruses. Nature 2016; 537(7622):689-93. Epub
2016/09/23. https://doi.org/10.1038/nature19366 PMID: 27654921

44. 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. Epub 2016/08/18. https://doi.org/10.1038/nature19094 PMID: 27533034

45. Callanan J, Stockdale SR, Shkoporov A, Draper LA, Ross RP, Hill C. Expansion of known ssRNA phage genomes: From tens to over a thousand. Sci Adv. 2020; 6(6):eaay5981. Epub 2020/02/23. https://doi.org/10.1126/sciadv.aay5981 PMID: 32083183.

46. Shkoporov AN, Clooney AG, Sutton TDS, Ryan FJ, Daly KM, Nolan JA, et al. The Human Gut Virome Is Highly Diverse, Stable, and Individual Specific. Cell Host Microbe 2019; 26(4):527-41.e5. Epub 2019/10/11. https://doi.org/10.1016/j.chom.2019.09.009 PMID: 31600503

47. Aiewsakun P, Adriaenssens EM, Lavigne R, Kropinski AM, Simmonds P. Evaluation of the genomic diversity of viruses infecting bacteria, archaea and eukaryotes using a common bioinformatic platform: steps towards a unified taxonomy. J Gen Virol. 2018; 99(9):1331-43. Epub 2018/07/18. https://doi.org/10.1099/jgv.0.001110 PMID: 30016225.

48. Aiewsakun P, Simmonds P. The genomic underpinnings of eukaryotic virus taxonomy: creating a sequence-based framework for family-level virus classification. Microbiome. 2018; 6(1):38. Epub 2018/02/21. https://doi.org/10.1186/s40168-018-0422-7 PMID: 29458427.
49. Bin Jang H, Bolduc B, Zablocki O, Kuhn JH, Roux S, Adriaenssens EM, et al. Taxonomic assignment of uncultivated prokaryotic virus genomes is enabled by gene-sharing networks. Nat Biotechnol 2019; 37(6):632-9. Epub 2019/05/08. https://doi.org/10.1038/s41587-019-0100-8 PMID: 31061483

50. Barylski J, Enault F, Dutilh BE, Schuller MB, Edwards RA, Gillis A, et al. Analysis of Spounaviruses as a Case Study for the Overdue Reclassification of Tailed Phages. Syst Biol. 2020; 69(1):110-23. Epub 2019/05/28. https://doi.org/10.1093/sysbio/syz036 PMID: 31127947

51. Gabler F, Nam SZ, Till S, Mirdita M, Steinegger M, Soding J, et al. Protein Sequence Analysis Using the MPI Bioinformatics Toolkit. Curr Protoc Bioinformatics 2020; 72(1):e108. Epub 2020/12/15. https://doi.org/10.1002/cpbi.108 PMID: 33315308

52. Makarova KS, Wolf YI, Koonin EV. Archaeal Clusters of Orthologous Genes (arCOGs): An Update and Application for Analysis of Shared Features between Thermococcales, Methanococcales, and Methanobacteriales. Life (Basel). 2015; 5(1):818-40. Epub 2015/03/13. https://doi.org/10.3390/life5010818 PMID: 25764277

53. Baranyi U, Klein R, Lubitz W, Kruger DH, Witte A. The archaeal halophilic virus-encoded Dam-like methyltransferase M. phiCh1-I methylates adenine residues and complements dam mutants in the low salt environment of Escherichia coli. Mol Microbiol 2000; 35(5):1168-79. Epub 2000/03/11. https://doi.org/10.1046/j.1365-2958.2000.01786.x PMID: 10712697

54. Hutinet G, Kot W, Cui L, Hillebrand R, Balamkundu S, Gnanakalai S, et al. 7-Deazaguanine modifications protect phage DNA from host restriction systems. Nat Commun. 2019; 10(1):5442. Epub 2019/12/01. https://doi.org/10.1038/s41467-019-13384-y PMID: 31784519

55. Lopatina A, Tal N, Sorek R. Abortive Infection: Bacterial Suicide as an Antiviral Immune Strategy. Annu Rev Virol 2020; 7(1):371-84. Epub 2020/06/20. https://doi.org/10.1146/annurev-virology-011620-040628 PMID: 32559405

56. 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. Epub 2019/02/03. https://doi.org/10.1186/s12985-019-1120-1 PMID: 30709355

57. Gross M, Marianovsky I, Glaser G. MazG—a regulator of programmed cell death in Escherichia coli. Mol Microbiol 2006; 59(2):590-601. Epub 2006/01/05. https://doi.org/10.1111/j.1365-2958.2005.04956.x PMID: 16390452

58. Fukui T, Shiomi N, Doi Y. Expression and characterization of (R)-specific enoyl coenzyme A hydratase involved in polyhydroxyalkanoate biosynthesis by Aeromonas caviae. J Bacteriol. 1998; 180(3):667-73. Epub 1998/02/11. https://doi.org/10.1128/JB.180.3.667-673.1998 PMID: 9457873

59. Charlier D, Nguyen Le Minh P, Roovers M. Regulation of carbamoylphosphate synthesis in Escherichia coli: an amazing metabolite at the crossroad of arginine and pyrimidine biosynthesis. Amino Acids. 2018; 50(12):1647-61. Epub 2018/09/22. https://doi.org/10.1007/s00726-018-2654-z PMID: 30238253

60. Crummett LT, Puxty RJ, Weihe C, Marston MF, Martiny JBH. The genomic content and context of auxiliary metabolic genes in marine cyanomyoviruses. Virology 2016; 499:219-29. Epub 2016/10/21. https://doi.org/10.1016/j.virol.2016.09.016 PMID: 27693926

61. Chevallereau A, Blasdel BG, De Smet J, Monot M, Zimmermann M, Kogadeeva M, et al. Next-Generation "-omics" Approaches Reveal a Massive Alteration of Host RNA Metabolism during Bacteriophage Infection of Pseudomonas aeruginosa. PLoS Genet. 2016; 12(7):e1006134. Epub 2016/07/06. https://doi.org/10.1371/journal.pgen.1006134 PMID: 27380413

62. Bianchi V, Pontis E, Reichard P. Interrelations between substrate cycles and de novo synthesis of pyrimidine deoxyribonucleoside triphosphates in 3T6 cells. Proc Natl Acad Sci U S A. 1986; 83(4):986-90. Epub 1986/02/01. https://doi.org/10.1073/pnas.83.4.986 PMID: 3456577

63. Rodionov DA, Vitreschak AG, Mironov AA, Gelfand MS. Comparative genomics of the vitamin B12 metabolism and regulation in prokaryotes. J Biol Chem 2003; 278(42):41148-59. Epub 2003/07/19. https://doi.org/10.1074/jbc.M305837200 PMID: 12869542

64. Lundqvist J, Elmlund D, Heldt D, Deery E, Soderberg CA, Hansson M, et al. The AAA(+) motor complex of subunits CobS and CobT of cobaltochelatase visualized by single particle electron microscopy. J Struct Biol 2009; 167(3):227-34. Epub 2009/06/24. https://doi.org/10.1016/j.jsb.2009.06.013 PMID: 19545636

65. Antonov IV. Two Cobalt Chelatase Subunits Can Be Generated from a Single chID Gene via Programed Frameshifting. Mol Biol Evol 2020; 37(8):2268-78. Epub 2020/03/27. https://doi.org/10.1093/molbev/msaa081 PMID: 32211852

66. Ignacio-Espinoza JC, Sullivan MB. Phylogenomics of T4 cyanophages: lateral gene transfer in the 'core' and origins of host genes. Environ Microbiol 2012; 14(8):2113-26. Epub 2012/02/22. https://doi.org/10.1111/j.1462-2920.2012.02704.x PMID: 22348436
67. Pfeifer F, Dyall-Smith M. Open issues for protein function assignment in Haloferax volcanii and other halophilic archaea. *Genes* (Basel). 2021; 12(7):963. https://doi.org/10.3390/genes12070963 PMID: 34202810

68. van Tran N, Muller L, Ross RL, Lestini R, Letoquart J, Ulryck N, et al. Evolutionary insights into Trm112-methyltransferase holoenzymes involved in translation between archaea and eukaryotes. *Nucleic Acids Res.* 2018; 46(16):8483-99. Epub 2018/07/17. https://doi.org/10.1093/nar/gky638 PMID: 30010922.

69. Letoquart J, van Tran N, Caroline V, Aleksandrov A, Lazar N, van Tilbeurgh H, et al. Insights into molecular plasticity in protein complexes from Trm9-Trm112 tRNA modifying enzyme crystal structure. *Nucleic Acids Res.* 2015; 43(22):10989-1002. Epub 2015/10/07. https://doi.org/10.1093/nar/gkv1009 PMID: 26438534.

70. Klein DJ, Moore PB, Steitz TA. The roles of ribosomal proteins in the structure assembly, and evolution of the large ribosomal subunit. *J Mol Biol* 2004; 340(1):141-77. Epub 2004/06/09. https://doi.org/10.1016/j.jmb.2004.03.076 PMID: 15184028

71. Mizuno CM, Guyomar C, Roux S, Lavigne R, Rodriguez-Valera F, Sullivan MB, et al. Numerous cultivated and uncultivated viruses encode ribosomal proteins. *Nat Commun.* 2019; 10(1):752. Epub 2019/02/16. https://doi.org/10.1038/s41467-019-08672-6 PMID: 30765709.

72. Klein R, Rossler N, Iro M, Scholz H, Witte A. Haloarchaeal myovirus phiCh1 harbours a phase variation system for the production of protein variants with distinct cell surface adhesion specificities. *Mol Microbiol* 2012; 83(1):137-50. Epub 2011/11/25. https://doi.org/10.1111/j.1365-2958.2011.07921.x PMID: 22111759

73. Trojet SN, Caumont-Sarcos A, Perrody E, Comeau AM, Krisch HM. The gp38 adhesins of the T4 superfamily: a complex modular determinant of the phage's host specificity. *Genome Biol Evol.* 2011; 3:674-86. Epub 2011/07/13. https://doi.org/10.1093/gbe/evr059 PMID: 21746838.

74. Roux S, Enault F, Ravet V, Colombet J, Bettarel Y, Auguet JC, et al. Analysis of metagenomic data reveals common features of halophilic viral communities across continents. *Environ Microbiol* 2016; 18 (3):889-903. Epub 2015/10/17. https://doi.org/10.1111/1462-2920.13084 PMID: 26472517

75. Garcia-Heredia I, Martin-Cuadrado AB, Mojica FJ, Santos F, Mira A, Anton J, et al. Reconstructing viral genomes from the environment using fosmid clones: the case of haloviruses. *PLoS ONE.* 2012; 7 (3):e33802. Epub 2012/04/06. https://doi.org/10.1371/journal.pone.0033802 PMID: 22479446.

76. Fendrihan S, Legat A, Pfaffenhuemer M, Gruber C, Weidler G, Gerbl F, et al. Extremely halophilic archaea and the issue of long-term microbial survival. *Rev Environ Sci Biotechnol.* 2006; 5(2-3):203-18. Epub 2006/08/01. https://doi.org/10.1007/s11157-006-0007-y PMID: 21984879.

77. Oren A. The ecology of the extremely halophilic archaea. *FEMS Microbiol Rev.* 1994; 13 (4):415-39.

78. Baquero DP, Contursi P, Piochi M, Bartolucci S, Liu Y, Cvirkaite-Krupovic V, et al. New virus isolates from Italian hydrothermal environments underscore the biogeographic pattern in archaeal virus communities. *ISME J.* 2020; 14(7):1821-33. Epub 2020/04/24. https://doi.org/10.1038/s41396-020-0653-z PMID: 32322010.

79. Dyall-Smith M, Pfeiffer F, Chiang PW, Tang SL. The novel halovirus Hardycor1, and the presence of active (induced) proviruses in four haloarchaea. *Genes* (Basel). 2021; 12(2):149. Epub 2021/01/28. https://doi.org/10.3390/genes12020149 PMID: 33498646.

80. Thiroux S, Dupont S, Nesbo 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. Epub 2020/10/07. https://doi.org/10.1111/1462-2920.15271 PMID: 33022088

81. Senčilo A, Roine E. A Glimpse of the genomic diversity of haloarchaeal tailed viruses. *Front Microbiol.* 2014; 5:84. Epub 2014/03/25. https://doi.org/10.3389/fmicb.2014.00084 PMID: 24659986.

82. Dunne M, Denyes JM, Arndt H, Loessner MJ, Leiman PG, Klumpp J. Salmonella Phage S16 Tail Fiber Adhesin Features a Rare Polyglycine Rich Domain for Host Recognition. *Structure* (London, England: 1993). 2018; 26(12):1573-82.e4. Epub 2018/09/25. https://doi.org/10.1016/j.str.2018.07.017 PMID: 30244968

83. Albers SV, Meyer BH. The archaeal cell envelope. *Nat Rev Microbiol* 2011; 9(6):414-26. Epub 2011/05/17. https://doi.org/10.1038/nrmicro2576 PMID: 21572458

84. Gunde-Cimerman N, Plemenitas A, Oren A. Strategies of adaptation of microorganisms of the three domains of life to high salt concentrations. *FEMS Microbiol Rev* 2018; 42(3):353-75. Epub 2018/03/13. https://doi.org/10.1093/femsre/fuy009 PMID: 29529204

85. Garneau JR, Depardieu F, Fortier LC, Bikard D, Monot M. PhageTerm: a tool for fast and accurate determination of phage termini and packaging mechanism using next-generation sequencing data. *Sci Rep.* 2017; 7(1):8292. Epub 2017/08/16. https://doi.org/10.1038/s41598-017-07910-5 PMID: 28811656.
86. Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014; 42(Data-base issue):D206-14. Epub 2013/12/03. https://doi.org/10.1093/nar/gkt1226 PMID: 24293654.

87. Yoon SH, Ha SM, Lim J, Kwon S, Chun J. A large-scale evaluation of algorithms to calculate average nucleotide identity. Antonie Van Leeuwenhoek 2017; 110(10):1281-6. Epub 2017/02/17. https://doi.org/10.1007/s10482-017-0844-4 PMID: 28204908

88. McWilliam H, Li W, Uludag M, Squizzato S, Park YM, Buso N, et al. Analysis Tool Web Services from the EMBL-EBI. Nucleic Acids Res. 2013; 41(Web Server issue):W597-600. Epub 2013/05/15. https://doi.org/10.1093/nar/gkt376 PMID: 23671338.

89. Pei J, Kim BH, Grishin NV. PROMALS3D: a tool for multiple protein sequence and structure alignments. Nucleic Acids Res. 2008; 36(7):2295-300. Epub 2008/02/22. https://doi.org/10.1093/nar/gkn072 PMID: 18287115.

90. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol 2010; 59(3):307-21. Epub 2010/06/09. https://doi.org/10.1093/sysbio/syq010 PMID: 20525638

91. Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, et al. Highly accurate protein structure prediction with AlphaFold. Nature. 2021; 596(7873):583-9. Epub 2021/07/16. https://doi.org/10.1038/s41586-021-03819-2 PMID: 34265844.

92. Baek M, DiMaio F, Anishchenko I, Dauparas J, Ovchinnikov S, Lee GR, et al. Accurate prediction of protein structures and interactions using a three-track neural network. Science (New York, NY). 2021; 373(6557):871-6. Epub 2021/07/21. https://doi.org/10.1126/science.abj8754 PMID: 34282049

93. Holm L, Laakso LM. Dali server update. Nucleic Acids Res. 2016; 44(W1):W351-5. Epub 2016/05/01. https://doi.org/10.1093/nar/gkw357 PMID: 27131377.

94. Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011; 27(7):1009-10. Epub 2011/02/01. https://doi.org/10.1093/bioinformatics/btr039 PMID: 21278367.

95. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009; 19(9):1639-45. Epub 2009/06/23. https://doi.org/10.1101/gr.092759.109 PMID: 19541911.

96. Yamada KD, Tomii K, Katoh K. Application of the MAFFT sequence alignment program to large data-examination of the usefulness of chained guide trees. Bioinformatics. 2016; 32(21):3246-51. Epub 2016/10/30. https://doi.org/10.1093/bioinformatics/btw412 PMID: 27378296.

97. Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009; 25(15):1972-3. Epub 2009/06/10. https://doi.org/10.1093/bioinformatics/btp348 PMID: 19505945.

98. Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009; 26(7):1641-50. Epub 2009/04/21. https://doi.org/10.1093/molbev/msp077 PMID: 19377059.

99. Yu G, Lam TT, Zhu H, Guan Y. Two Methods for Mapping and Visualizing Associated Data on Phylogeny Using Ggtree. Mol Biol Evol. 2018; 35(12):3041-3. Epub 2018/10/24. https://doi.org/10.1093/molbev/msy194 PMID: 30351396.

100. Mukherjee S, Stamatis D, Bertsch J, Ovchinnikova G, Sundaramurthi JC, Lee J, et al. Genomes OnLine Database (GOLD) v.8: overview and updates. Nucleic Acids Res. 2021; 49(D1):D723-D33. Epub 2020/11/06. https://doi.org/10.1093/nar/gkaa983 PMID: 33152092.

101. Roux S, Paez-Espino D, Chen IA, Palaniappan K, Ratner A, Chu K, et al. IMG/VR v3: an integrated ecological and evolutionary framework for interrogating genomes of uncultivated viruses. Nucleic Acids Res. 2021; 49(D1):D764-D75. Epub 2020/11/03. https://doi.org/10.1093/nar/gkaa946 PMID: 33137183.