Next Article in Journal
Association of Increased Programmed Death Ligand 1 Expression and Regulatory T Cells Infiltration with Higher Hepatocellular Carcinoma Recurrence in Patients with Hepatitis B Virus Pre-S2 Mutant after Curative Surgical Resection
Next Article in Special Issue
Diversity of SIRV-like Viruses from a North American Population
Previous Article in Journal
Repurposing Molnupiravir for COVID-19: The Mechanisms of Antiviral Activity
Previous Article in Special Issue
Halorubrum pleomorphic virus-6 Membrane Fusion Is Triggered by an S-Layer Component of Its Haloarchaeal Host
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Viral Susceptibility of the Haloferax Species

by
Zaloa Aguirre Sourrouille
1,†,
Sabine Schwarzer
1,2,†,
Sebastian Lequime
3,
Hanna M. Oksanen
4,* and
Tessa E. F. Quax
1,2,*
1
Biology of Archaea and Viruses, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, 9747 AG Groningen, The Netherlands
2
Archaeal Virus-Host Interactions, Faculty of Biology, University of Freiburg, 97104 Freiburg im Breisgau, Germany
3
Cluster of Microbial Ecology, Groningen Institute for Evolutionary Life Sciences, University of Groningen, 9747 AG Groningen, The Netherlands
4
Molecular and Integrative Biosciences Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, 00014 Helsinki, Finland
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Viruses 2022, 14(6), 1344; https://doi.org/10.3390/v14061344
Submission received: 16 May 2022 / Revised: 15 June 2022 / Accepted: 17 June 2022 / Published: 20 June 2022
(This article belongs to the Special Issue Archaeal Virology)

Abstract

:
Viruses can infect members of all three domains of life. However, little is known about viruses infecting archaea and the mechanisms that determine their host interactions are poorly understood. Investigations of molecular mechanisms of viral infection rely on genetically accessible virus–host model systems. Euryarchaea belonging to the genus Haloferax are interesting models, as a reliable genetic system and versatile microscopy methods are available. However, only one virus infecting the Haloferax species is currently available. In this study, we tested ~100 haloarchaeal virus isolates for their infectivity on 14 Haloferax strains. From this, we identified 10 virus isolates in total capable of infecting Haloferax strains, which represented myovirus or siphovirus morphotypes. Surprisingly, the only susceptible strain of all 14 tested was Haloferax gibbonsii LR2-5, which serves as an auspicious host for all of these 10 viruses. By applying comparative genomics, we shed light on factors determining the host range of haloarchaeal viruses on Haloferax. We anticipate our study to be a starting point in the study of haloarchaeal virus–host interactions.

1. Introduction

Microbial viruses are widespread and able to infect members of all three domains of life, including archaea. Archaea are ubiquitous microorganisms that can be found in extreme environments, such as salt lakes as well as in mesophilic surroundings such as the oceans and the human body [1,2]. The study of archaeal viruses is essential to understand the origin and the evolution of viruses in general [3]. Archaeal viruses display a high genomic and structural variability, but they also share some common traits with viruses infecting other domains of life [4,5,6,7]. Viruses are divided into different families, currently mainly based on the sequence similarities and the highest taxonomic ranks; realms largely follow groupings based on the characteristics of major virus capsid proteins or genome replication components [5,8,9]. Whereas crenarchaeal viruses come in many different shapes, the majority of viruses infecting euryarchaea display a head–tail morphology and are currently members of 14 families in the class Caudoviricetes [6,9,10,11,12]. The rest of the currently known euryarchaeal viruses are either internal membrane-containing tailless icosahedral (family Sphaerolipoviridae) [13], pleomorphic (family Pleolipoviridae) [14], or spindle-shaped (family Halspiviridae) [15]. Archaeal tailed viruses are the most common isolates infecting halophilic archaea. They are morphologically indistinguishable from tailed double-stranded (ds) DNA bacteriophages that have the myovirus (long and contractile helical tail), podovirus (short tails), or siphovirus (long and non-contractile tails) morphology. At the sequence level, however, archaeal tailed viruses are very diverse, several of them are singletons, and they hardly resemble their bacterial relatives [10,16]. All isolated archaeal viruses have a DNA genome so far, and the majority of their genes encode proteins of unknown function showing limited or nonexistent similarity to tailed bacteriophage proteins, as a result many aspects of archaeal virus–host relationships and virus life cycles remain unknown [6]. However, very recently, comparative genomics and a host range analysis of the tailed archaeal viruses showed the role of the tail fiber adhesin in host recognition [10]. The adhesins from archaeal viruses resemble the adhesins located at the distal tip of the tail fibers of various members of the T-even phage group [17]. The structural core of the adhesins is formed by highly conserved glycine-rich motifs that separate the hypervariable segments [17]. While these conserved glycine-rich domains are used for binding, mutations or shuffling of the hypervariable regions change adhesin receptor specificity and thus primarily determine the host range [10,18]. Previous screenings of haloarchaeal viruses and their hosts highlighted the high abundance of myovirus isolates and their extremely broad host ranges. The Hafunaviridae is the largest family of archaeal tailed viruses and its myoviruses have a broad host range [10,19]. Other haloarchaeal viruses are more specific to a certain host [10,18].
The study of the virus–host relationships and infection mechanisms of archaeal viruses would greatly benefit from the availability of genetically accessible virus–host models, for which molecular biology tools are available. For crenarchaeal virus hosts, several members of the genus Sulfolobales are currently the most common archaeal models with genetic systems available [20,21,22,23,24]. Moreover, there are already some genetic systems for crenarchaeal viruses available such as SSVs, STIV, and STSV1 [25,26,27,28]. The study of euryarchaeal viruses presently relies mainly on haloarchaeal Halorubrum and Haloarcula strains that are infected by a substantial number of known viruses [10,19,29]. These hosts have been used successfully to study archaeal viruses and their structures, entry, and egress mechanisms [30,31,32,33,34]. However, the genetic and the molecular toolset for these hosts is limited. In recent years, Haloferax has become increasingly popular in the archaeal scientific community, and it is the euryarchaeal model for which the most advanced tools for genetic engineering, imaging, and molecular biology are available [35]. The Haloferax tools entail a versatile genetic system for overexpression and genomic knock-out, the availability of several plasmids and different markers, and a CRISPR-based repression system to downregulate gene expression [23,36,37,38,39]. In addition, it is the only archaeon for which several fluorescent fusion proteins are available [40]. It is an excellent organism for light microscopy, and it is also used in microfluidics [40]. The development of these technical advances, and the growing scientific community embracing Haloferax as a central euryarchaeal model, has led to a substantial increase in the understanding of its cell biology [41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56]. This detailed knowledge of euryarchaeal cell biology, obtained using Haloferax as a model, is of great added value in studies of viral infection mechanisms.
In contrast to some other haloarchaea, there are far less viruses known that infect Haloferax strains. Almost 30 years ago, the HF1 virus was reported to infect Haloferax lucentense and Hfx. volcanii, and a defective provirus of Haloferax mediterranei has been identified, both of which are no longer available (M. Dyall-Smith, personal communication) [57,58]. At the moment, Haloferax tailed virus 1 (HFTV1) is only one virus isolated from a Haloferax host, and it was recently isolated together with its host Haloferax gibbonsii LR2-5 from the hypersaline Lake Retba in Senegal [59,60].
Due to the attractiveness of Haloferax for molecular studies, we aim to identify viruses that infect Haloferax strains. In this study, we used the largest available collection of isolated and characterized haloarchaeal virus isolates, and we tested the infectivity of 95 viruses on 14 Haloferax strains. The virus collection contains viruses from all current haloarchaeal virus families, and these viruses represent the majority of haloarchaeal virus isolates isolated to date. This endeavor resulted in an extended virus–host matrix for Haloferax from which we could identify a promising model that is the host to a substantial number of viruses, Haloferax gibbonsii LR2-5. In addition, we used comparative genomics to identify viral factors that allow for the infection of this Haloferax host. With this work, we pave the way to using Haloferax as a model in virus–host interactions.

2. Materials and Methods

2.1. Archaeal Viruses and Strains and Growth Conditions

All Haloferax strains (Supplementary Table S1), virus host strains (Supplementary Table S2), and viruses (Supplementary Table S2) were grown aerobically at 37 °C in a modified growth medium (MGM) [58,61]. The artificial 30% salt water (SW) (240 g NaCl, 30 g MgCl2 ·6H2O, 35 g MgSO4 · 7H2O, 7 g KCl, 5 mL of 1 M CaCl2 · 2H2O, and 80 mL of 1 M Tris-HCl, pH 7.2 per L) was diluted to obtain 18, 20, or 23% SW in the top-agar media, plates, and liquid media, respectively. MGM also contained 5 g of peptone (Oxoid), 1 g of Bacto yeast extract (Becton, Dickinson and Company, Sparks, MD, USA), and Bacto agar (14 g for plates; 4 g per top-layer; Becton, Dickinson and Company) per liter. Viruses were grown on their own host strain listed in Supplementary Table S2 by using a double-layer plaque assay. For plaque assays, the strains were grown over 2–3 nights in liquid media to obtain a dense culture (OD~1) of which 200–300 µL were used per plate for inoculation to obtain an even layer of dense growth of the strain on a soft agar layer. Viruses diluted in MGM broth were added 100 µL per plate. The virus and the strain were combined with melted soft agar (3–4 mL per plate; 50 °C), mixed, and plated. The plates were grown for 2–4 days in a box (in some cases with an additional cup of water) to prevent them from drying out. The plaques observed on plates were used to enumerate the viruses in the samples by taking into account the virus dilution, and plates with confluent or semi-confluent growth of virus were used to produce virus lysates. Virus lysates were prepared by mixing the collected soft agar layer from confluent or semi-confluent plates with 2–3 mL of MGM broth per plate. After 1.5 h of shaking at 37 °C, lysates were cleared from cell debris and agar by centrifugation (Sorvall F14 rotor, 10,000 rpm, 20 min, 4 °C). Lysates were stored at +4 °C for up 2 months before use.

2.2. Sensitivity of Haloferax Strains to Euryarchaeal Viruses

The sensitivity of Haloferax strains (Supplementary Table S1) to viruses (Supplementary Table S2) was determined by placing 10 µL drops of undiluted and 1:100 diluted virus lysates applied to Haloferax strain plated with top-agar on a plate using a double-layer plaque assay. The host strains of the viruses (Supplementary Table S2) were used as positive controls. For virus drops, MGM was used as a negative control. The drops were repeated in duplicate. All cases where cell growth was inhibited on a plate in the presence of virus lysate were further tested by a double-layer plaque assay to confirm viral infection and to obtain the numerical values of efficiency of infection. All titer data was collected at 37 °C.

2.3. Viral Comparative Genomics

The genomic sequences of the viruses were retrieved from NCBI [62]. Phylogenetic trees were constructed using a custom-written R script with the package ggtree [63], using the ANI values calculated with VIRIDIC [64]. Whole genome sequence alignments were visualized with Easyfig [65]. To find putative virus–host determinants, protein sequences found at variable regions of the viral genomes were downloaded from NCBI, aligned using MAFFT v. 7.450 with default settings [66], and visualized by Jalview. Under-represented tetramers were detected using https://www.cmbl.uga.edu//software/signature.html (accessed on 15 December 2021), and for more complex motifs a manual search was carried out in Geneious v 8.1.9.(Auckland, New Zeeland).
The phylogenetic analyses of the adhesin and tail fiber gene sequences were carried out as follows. For both datasets, sequences were aligned using MAFFT v. 7.450 [66], and a maximum likelihood phylogenetic tree constructed with IQ-TREE v. 1.6.12 [67], using ModelFinder for model selection [68]. Additional phylogenetic analyses were conducted in BEAST 1.10.4 [69], using the best nucleotide substitution model as determined in ModelFinder and available in BEAST, namely WAG + F + G4 [70] and BLOSUM62 + F + G4 [71] for the adhesins and tail fibers gene sequences, respectively. All sequences were considered isochronous, and the evolutionary process was reconstructed under a strict molecular clock with a fixed rate of 1 and a constant population size model [72]. The infection capability of the LR2-5 strain was considered as a discrete trait analyzed using a symmetric diffusion model [73], and the number of changes between the two states (Markov jumps) was estimated in the posterior distribution of trees [74]. Proper convergence and mixing (effective sample size > 200) was verified using tracer v. 1.7 [75] and the burn-in (10% of samples) was removed. The presence of a phylogenetic signal linked to the LR2-5 infection capability was tested by comparing the distribution of Markov jumps to a “null distribution” of an estimated number of Markov jumps computed by a set of 10 independent runs, where the states (i.e., LR2-5 infection capability) were randomized. A clear overlap between the 95% highest posterior density (95% HPD) of the estimated number of Markov jumps in the original dataset and the randomized “null” distributions would sign the absence of a phylogenetic signal [73]. If not otherwise stated, default parameters were used.

3. Results and Discussion

3.1. Detection of Viruses Infecting Haloferax Strains

To identify novel virus–host models for Haloferax, a collection of 95 haloarchaeal viruses (Supplementary Table S2) were cross-tested with 14 Haloferax strains (Supplementary Table S1, Figure 1). The viruses were isolated from samples taken from different hypersaline environments, and they are part of the collection at the University of Helsinki. Information on the origins, host strains, virus morphologies, genomes, and taxonomic classification of the viruses can be found in Supplementary Table S2. Some of the 14 tested Haloferax strains are widely used laboratory models, such as Hfx. volcanii H26 and Hfx. mediterranei [76], while several other strains were isolated quite recently from the hypersaline Lake Retba in Senegal (Supplementary Table S2) [59]. A scheme of the experimental pipeline for the screening is depicted in Figure 1.
The screening resulted in the detection of 10 virus isolates that could make plaques on the Haloferax strains tested (Table 1). Curiously, Hfx. gibbonsii LR2-5 was the only susceptible Haloferax strain. The LR2-5 infecting viruses belong to the genera Haloferacalesvirus, Mincapvirus (family Hafunaviridae) or Retbasiphovirus (family Haloferuviridae) (Table 1). It is important to note that during the course of the study, HCTV-6 and HCTV-13 genome sequences were reported to be identical [10], and, as a result, nine unique virus isolates infecting Haloferax strains were identified. Furthermore, the infectivity of the closely related virus isolates assigned to the three genera were re-tested with LR2-5 by plaque assay, but no more interactions were detected (Table 1).
In several cases, however, we observed a growth inhibition zone by spot-on-lawn assay but no plaques during the plaque assay (Table S2). As haloarchaea are known to produce antimicrobial toxins, i.e., halocins [77,78], we assume that halocins produced by the virus host cell are present in some of the virus stocks. These halocins would result in a growth inhibition in spot-on-lawn assays but not result in plaques. Only those virus–host pairs for which plaques were detected were marked as positive and true virus–host pairs. The quantitative plaque assay also allowed for the determination of the viral titer on the Haloferax host, which is an indication of the efficiency of infection. As a summary, we did not find any true virus–host pairs on any of the Haloferax strains, except on strain Haloferax gibbonsii LR2-5, which was isolated a few years ago from Lake Retba in Senegal at the same time as HFTV1 [59,60]. All LR2-5 infecting viruses had either Halorubrum or Haloarcula as their own host and their titers varied between 103 and 1010 PFU/mL on LR2-5 (Table 1). On the other hand, we also calculated the efficiency of plating (EOP) of the viruses that were able to infect LR2-5 to compare them with their original isolation host. The efficiency of plating (EOP) of three mincapviruses HSTV-2, HRTV-2, and HCTV-6 originally grown either on Halorubrum or Haloarcula hosts had a slightly higher or around same EOP on LR2-5 than on their own host. However, the rest of the LR2-5 infecting viruses had three to six magnitudes lower EOPs on LR2-5 than on their own host. These results are in accordance with previous observations, which showed that myoviruses infecting halophilic archaea can have a wide host range [10,29]

3.2. Characteristics of Haloferax Infecting Viruses

To map the determinants of viral host range, we focused on viruses that infect LR2-5 and those viruses belonging to the same genus as LR2-5 infecting viruses (genera Haloferacalesvirus, Mincapvirus, Retbasiphovirus; Table 1). A phylogenetic tree based on the complete genome sequences of the haloferacalesviruses, mincapviruses, and retbasiphovirus (Figure 2) and their EOPs on LR2-5 (Table 1) revealed that all mincapviruses, except HRTV-11, can infect LR2-5.
The LR2-5 infecting viruses have been isolated from hypersaline environmental samples either from Israel, Italy, Senegal, Slovenia, or Thailand (Figure 3A, Supplementary Table S2). All viruses represented tailed morphologies (Table 1). HFTV1 is the only siphovirus and the rest of the LR2-5 infecting viruses are myoviruses (Table 1). A schematic of the morphologies and a typical genome organization of myovirus and siphovirus are shown in Figure 3B.
HRTV-11 is the only mincapvirus isolated from a remote location, Slovenia, where no other haloarchaeal virus were successfully isolated (Figure 3A; [29]). The broad host range of myoviruses has been linked to their ability to exchange their host-specific genetic modules for receptor binding proteins [19]. Thus, we hypothesize that the host range of viruses isolated from locations with a low virus-density might be different from other closely related viruses isolated from virus-dense environments, as the possibility for recombination events might have been limited. We had a closer look into possible host range determining factors to find potential explanations for the few haloferacalesviruses that exceptionally do infect LR2-5.

3.3. Restriction–Modification Systems

Viruses are known to develop strategies to escape host defense mechanisms. Thus, we explored in more detail the antiviral defense mechanisms of the host and the viral escape mechanisms to determine if any of these factors could explain the differences in the ability to infect LR2-5 between viruses from the same family. LR2-5 does not have a CRISPR-cas system but encodes a predicted type I restriction modification (RM) system [60]. These antiviral mechanisms are based on methylation of host DNA (to protect it) and cleavage of unmethylated DNA (viral DNA). Part of the LR2-5 encoded RM system are a Zim methylase (CTAG methyltransferase), a Mrr-like endonuclease, and an RmeRMS (type I restriction enzyme restriction/methylation/specificity subunit).
One strategy of viruses to escape host recognition is the avoidance of certain motifs in their genomes, which are the targets for the RM systems [79]. In line with previous studies [79], the palindromic tetrameric motifs CTAG, GATC, and AGCT are absent in the mincapvirus genomes, whereas haloferacalesviruses lack only CTAG and GATC, except for HRTV-22, which contains a CTAG motif (Supplementary Table S3). Furthermore, mincapviruses and haloferacalesviruses show an under-representation of the TGCA and the CATG motifs. On the other hand, those tetrameric motifs are also under-represented in the genomes of haloferuviruses (Supplementary Table S3). We did not observe any differences between the viruses with respect to under-represented motifs, which would be able to explain the different infectivity of the studied viruses on LR2-5.
Methylated motifs in the genome of LR2-5 were previously predicted [60]. We searched for sequences in viral genomes that might be recognized by the host methyltransferases. We hypothesize that when the viral genome is replicated, the host will recognize those sequences as “self” and methylate the viral DNA, hence, allowing the virus to escape the antiviral-mechanism. Calculation of the frequency of corresponding motifs in the viral genomes revealed that the motif GCGCTG is found more frequently in all mincapviruses than in haloferacalesviruses (Supplementary Table S4). The frequency of the other motifs was similar for all viruses. Therefore, we concluded that host range determinants in case of LR2-5 infecting viruses might rely on another factor.

3.4. Adhesins and Tail Fiber Proteins

Pairwise alignment of the viral genomes of the Hafunaviridae family members showed the same genome organization within their own group (mincapviruses or haloferacalesviruses) (Figure 4). There were only a few variable genomic regions observed. Firstly, at the right-end of the genomes, a region was detected with frequent insertions and inversions (Figure 4). Many of the proteins encoded by genes in this region are designated as hypothetical and they lack homology with proteins found in reference databases. It might be possible that these genes encode proteins involved in viral egress, as several of them have predicted transmembrane domains. However, since egress mechanisms of haloarchaeal viruses are not well understood, and the responsible proteins are not identified, it is not clear if these genes are determinants of differences observed in viral infectivity on LR2-5.
Another highly variable region in the viral genomes is located around genes encoding the viral tail fiber and adhesin (Figure 4 and Figure 5). Adhesins, located at the distal tip of the tail fibers, have been shown to be the host determinants in the tailed dsDNA bacteriophages of the order Caudoviricetes [17]. Long tail fibers of bacteriophages have a modular organization, and they consist of several different proteins [18,80]. Moreover, it has been suggested that haloarchaeal tailed viruses exchange the genes encoding either the tail fiber or adhesin with other viruses by recombination [10]. All haloferacalesviruses and mincapviruses have a gene encoding for a putative tail fiber protein together with a gene encoding a putative adhesin [10]. HFTV1 is a siphovirus with a long non-contractile tail, and it does not contain tail fibers.
A tree of the hafunaviruses (haloferacalesviruses and mincapviruses) tail fiber adhesin proteins showed that they are divided in four clades (Figure 5), as also shown by Liu et al. [10]. Clade 1 and Clade 3 adhesins were previously reported to correlate well with the observed virus host range among the tested Haloarcula, Halobacterium, Halobellus, Halorubrum, and Haloterrigena strains [10]. We aimed to test if this correlation between adhesins and host range could also be seen in our current data set.
We investigated the hypothetical relationship between the susceptibility of the LR2-5 strain to infection and diversification of the hafunavirus tail fiber and adhesin genes by comparing the number of Markov jumps (i.e., infecting or not LR2-5) through the adhesin and the tail fiber gene evolution (Figure 5). A “null” distribution for each gene was generated by randomizing the state’s distribution. We observed a clear overlap between the estimated numbers of Markov jumps (95% highest posterior density; HPD) in the real and the randomized datasets, indicating the lack of a phylogenetic signal linked to the viral infection capability on the LR2-5 strain for both the hafunavirus tail fiber gene and the adhesin gene (Figure 5).
Based on this analysis, we did not observe any clear correlation between the type of adhesin or tail fiber gene and the viral ability to infect Hfx. gibbonsii LR2-5. However, it should be emphasized that our dataset was limited to Haloferax hosts. Such a correlation between host range and adhesin type might only become apparent when a more diverse set of hosts is used, as the Haloferax strains are generally not very susceptible for viral infection. However, we did find a few case examples, which seems to pinpoint the adhesin as a host determining factor. First, HRTV-26, one of the few haloferacalesviruses that infects LR2-5, encodes a Clade 1 adhesin, similar to most mincapviruses. Since the adhesins of haloferacalesviruses that are closely related to HRTV-26 all belong to Clade 3, it is likely that HRTV-26 picked up an adhesin gene via horizontal gene transfer from one of the Clade 1 viruses. This might be the explanation for why HRTV-26, even though it is an haloferacalesvirus, is still capable of infecting LR2-5.
Second, viruses encoding a Clade 3 adhesin do not infect LR2-5, with the exception of HRTV-10. To gain more insight into the HRTV-10 tail fiber adhesin, a multiple sequence alignment was performed. We identified one mutated codon resulting in an amino acid change (T380A) in a conserved motif at the 3′end of the adhesin gene (Supplementary Figure S2). This subtle substitution might be a possible explanation for its increased infectivity on LR2-5, as this could produce alterations in the fold and the topology of adhesins. Liu et al. [10] already discussed how a single amino acid substitution (A217V) in the adhesins of HRTV-19 and HRTV-23 alters the host range. Moreover, Trojet et al. [17] showed how the long tail fiber locus in the T4 superfamily viruses is susceptible to frequent modular shuffling, which results in chimeric adhesins and, thus, viruses acquire new host receptor specificities. However, this hypothesis requires further experiments.

4. Conclusions

We aimed to develop an euryarchaeal virus–host model system with attractive molecular, genetic, and imaging tools to dive deeper into virus–host relationships in archaea. We focused on Haloferax strains, as the most advanced molecular and genetic tools are available for these organisms. Due to the low number of known viruses that infect Haloferax, we set up a large-scale assay to identify viruses capable of infecting Haloferax strains. Our extensive screening of approximately 100 isolated and mostly characterized haloarchaeal viruses showed that only a small subset was capable of infecting Haloferax. Specifically, Hfx. gibbonsii LR2-5 showed itself as an auspicious host, as we identified 10 virus isolates that can infect this host and all other hosts could not be infected. Previously, we studied potential factors that make Hfx. gibbonsii LR2-5 (DSM No. 112399) so much more receptive to virus infection in contrast to related strains. For its virus susceptibility, several possible factors were identified as potential explanations. The most prominent was an absence of CRISPR-cas virus defense systems in LR2-5 [60]. In addition, Hfx. gibbonsii LR2-5 has a different RM system from some other Haloferax strains, and the cell surface might show a different composition than that of other Haloferax strains because LR2-5 encodes divergent surface proteins that might serve as viral receptors, such as pili and S-layer proteins [60]. These latter two proteins are a common component of the archaeal cell-surface, and they are hypothesized to be used by viruses for initial binding and recognition [81,82].
Because we used plaque assays for the host range determination, we cannot distinguish at which step of the infection cycle viruses are not successful in infecting a host: entry, replication, or release. However, the results of the present study might indicate that adhesins, which are required for entry, are an important factor determining host range specificity, whereas other factors such as the viral egress proteins and restriction–modification systems might also play a role. Moreover, it should be noted that one limitation of the approach utilized is the assumption that the successful infection of the host results in plaque formation. It could be that some of the viruses studied here might have integrated their genome when infecting Hfx. gibbonsii LR2-5 and transferred to a lysogenic cycle. Consequently, these types of interactions would not show clear plaques, and they could have been missed by our approach. However, further experimental research on the host range implications of these proteins is needed.
Our analysis has identified Hfx. gibbonsii LR2-5 as the only virus susceptible Haloferax strain. This host can be infected by 10 different viruses with sequenced genomes [10,16], and it thus offers the possibility to compare infection mechanisms between viruses. Moreover, its genome is sequenced and its cell biology is characterized [60]. We are also currently developing a genetic system based on PyrE for Hfx. gibbonsii LR2-5. We conclude that this species is an extremely promising model host to study virus–host interactions in haloarchaea, and we anticipate that this work serves as a stepping stone for future in-depth molecular characterization of archaeal viral infection mechanisms.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14061344/s1, Table S1. Haloferax strains used in the study [19,29,37,59,60,76]; Table S2. Viruses used in the study; Table S3. Underrepresented palindromic motifs in viral genomes; Table S4. Manual counting of modified DNA motifs found in the genomes of viruses infecting Hfx. gibbonsii LR2-5; Figure S1. Multiple sequence alignment of adhesins belonging to the Group 3, Figure S2. VIRIDIC generated heatmap.

Author Contributions

T.E.F.Q. and H.M.O. conceptualized the research. T.E.F.Q., H.M.O., Z.A.S. and S.S. planned experiments and interpreted data. Z.A.S. and S.S. performed experiments. Z.A.S. and S.L. performed comparative genomics. All authors contributed to the writing of the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by funding from the Hector Fellow Academy to T.E.F.Q. and Z.A.S., an Emmy Noether grant (411069969) of the DFG to T.E.F.Q., and by funding from the University of Helsinki and the Academy of Finland for FINStruct and Instruct-FI, part of Biocenter Finland and Instruct-ERIC, respectively.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available within the article and its Supplementary Materials.

Acknowledgments

We would like to thank Sari Korhonen, Xueting Wang, and Zhengqun Li for practical support with the virus-host screen and M. Dyall-Smith for advice on the comparative genomics.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Borrel, G.; Brugère, J.-F.; Gribaldo, S.; Schmitz, R.A.; Moissl-Eichinger, C. The Host-Associated Archaeome. Nat. Rev. Microbiol. 2020, 18, 622–636. [Google Scholar] [CrossRef]
  2. Lloyd, K.G.; May, M.K.; Kevorkian, R.T.; Steen, A.D. Meta-Analysis of Quantification Methods Shows That Archaea and Bacteria Have Similar Abundances in the Subseafloor. Appl. Environ. Microbiol. 2013, 79, 7790–7799. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Forterre, P.; Prangishvili, D. The Origin of Viruses. Res. Microbiol. 2009, 160, 466–472. [Google Scholar] [CrossRef] [PubMed]
  4. Hartman, R.; Munson-McGee, J.; Young, M.J.; Lawrence, C.M. Survey of High-Resolution Archaeal Virus Structures. Curr. Opin. Virol. 2019, 36, 74–83. [Google Scholar] [CrossRef]
  5. Baquero, D.P.; Liu, Y.; Wang, F.; Egelman, E.H.; Prangishvili, D.; Krupovic, M. Structure and Assembly of Archaeal Viruses, 1st ed.; Elsevier Inc.: Amsterdam, The Netherlands, 2020; Volume 108, ISBN 9780128207611. [Google Scholar]
  6. Prangishvili, D.; Bamford, D.H.; Forterre, P.; Iranzo, J.; Koonin, E.V.; Krupovic, M. The Enigmatic Archaeal Virosphere. Nat. Rev. Microbiol. 2017, 15, 724–739. [Google Scholar] [CrossRef]
  7. Pietilä, M.K.; Demina, T.A.; Atanasova, N.S.; Oksanen, H.M.; Bamford, D.H. Archaeal Viruses and Bacteriophages: Comparisons and Contrasts. Trends Microbiol. 2014, 22, 334–344. [Google Scholar] [CrossRef]
  8. Abrescia, N.G.A.; Bamford, D.H.; Grimes, J.M.; Stuart, D.I. Structure Unifies the Viral Universe. Annu. Rev. Biochem. 2012, 81, 795–822. [Google Scholar] [CrossRef]
  9. Koonin, E.v.; Dolja, V.v.; Krupovic, M.; Varsani, A.; Wolf, Y.I.; Yutin, N.; Zerbini, F.M.; Kuhn, J.H. Global Organization and Proposed Megataxonomy of the Virus World. Microbiol. Mol. Biol. Rev. 2020, 84, e00061-19. [Google Scholar] [CrossRef]
  10. Liu, Y.; Demina, T.A.; Roux, S.; Aiewsakun, P.; Kazlauskas, D.; Simmonds, P.; Prangishvili, D.; Oksanen, H.M.; Krupovic, M. Diversity, Taxonomy, and Evolution of Archaeal Viruses of the Class Caudoviricetes. PLoS Biol. 2021, 19, e3001442. [Google Scholar] [CrossRef] [PubMed]
  11. Oren, A. Taxonomy of Halophilic Archaea: Current Status and Future Challenges. Extremophiles 2014, 18, 825–834. [Google Scholar] [CrossRef] [PubMed]
  12. Atanasova, N.S.; Bamford, D.H.; Oksanen, H.M. Haloarchaeal Virus Morphotypes. Biochimie 2015, 118, 333–343. [Google Scholar] [CrossRef] [Green Version]
  13. Demina, T.A.; Pietilä, M.K.; Svirskaitė, J.; Ravantti, J.J.; Atanasova, N.S.; Bamford, D.H.; Oksanen, H.M. HCIV-1 and Other Tailless Icosahedral Internal Membrane-Containing Viruses of the Family Sphaerolipoviridae. Viruses. MDPI AG 2017, 9, 32. [Google Scholar] [CrossRef]
  14. Demina, T.A.; Oksanen, H.M. Pleomorphic Archaeal Viruses: The Family Pleolipoviridae Is Expanding by Seven New Species. Arch. Virol. 2020, 165, 2723–2731. [Google Scholar] [CrossRef]
  15. Bath, C.; Dyall-Smith, M.L. His1, an Archaeal Virus of the Fuselloviridae Family That Infects Haloarcula Hispanica. J. Virol. 1998, 72, 9392–9395. [Google Scholar] [CrossRef]
  16. Senčilo, A.; Jacobs-Sera, D.; Russell, D.A.; Ko, C.C.; Bowman, C.A.; Atanasova, N.S.; Osterlund, E.; Oksanen, H.M.; Bamford, D.H.; Hatfull, G.F.; et al. Snapshot of Haloarchaeal Tailed Virus Genomes. RNA Biol. 2013, 10, 803. [Google Scholar] [CrossRef] [PubMed]
  17. Trojet, S.N.; Caumont-Sarcos, A.; Perrody, E.; Comeau, A.M.; Krisch, H.M. The Gp38 Adhesins of the T4 Superfamily: A Complex Modular Determinant of the Phage’s Host Specificity. Genome Biol. Evol. 2011, 3, 674–686. [Google Scholar] [CrossRef] [Green Version]
  18. Dunne, M.; Denyes, J.M.; Arndt, H.; Loessner, M.J.; Leiman, P.G.; Correspondence, J.K.; Klumpp, J. Salmonella Phage S16 Tail Fiber Adhesin Features a Rare Polyglycine Rich Domain for Host Recognition. Structure 2018, 26, 1573–1582. [Google Scholar] [CrossRef] [Green Version]
  19. Atanasova, N.S.; Demina, T.A.; Buivydas, A.; Bamford, D.H.; Oksanen, H.M. Archaeal Viruses Multiply: Temporal Screening in a Solar Saltern. Viruses 2015, 7, 1902–1926. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Li, Y.; Pan, S.; Zhang, Y.; Ren, M.; Feng, M.; Peng, N.; Chen, L.; Liang, Y.X.; She, Q. Harnessing Type i and Type III CRISPR-Cas Systems for Genome Editing. Nucleic Acids Res. 2015, 44, e34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Albers, S.V.; Driessen, A.J.M. Conditions for Gene Disruption by Homologous Recombination of Exogenous DNA into the Sulfolobus Solfataricus Genome. Archaea 2008, 2, 145–149. [Google Scholar] [CrossRef] [Green Version]
  22. Lewis, A.M.; Recalde, A.; Bräsen, C.; Counts, J.A.; Nussbaum, P.; Bost, J.; Schocke, L.; Shen, L.; Willard, D.J.; Quax, T.E.F.F.; et al. The Biology of Thermoacidophilic Archaea from the Order Sulfolobales. FEMS Microbiol. Rev. 2021, 45, fuaa063. [Google Scholar] [CrossRef]
  23. Leigh, J.A.; Albers, S.V.; Atomi, H.; Allers, T. Model Organisms for Genetics in the Domain Archaea: Methanogens, Halophiles, Thermococcales and Sulfolobales. FEMS Microbiol. Rev. 2011, 35, 577–608. [Google Scholar] [CrossRef] [Green Version]
  24. Wagner, M.; van Wolferen, M.; Wagner, A.; Lassak, K.; Meyer, B.H.; Reimann, J.; Albers, S.V. Versatile Genetic Tool Box for the Crenarchaeote Sulfolobus Acidocaldarius. Front. Microbiol. 2012, 3, 214. [Google Scholar] [CrossRef] [Green Version]
  25. Iverson, E.; Stedman, K. A Genetic Study of SSV1, the Prototypical Fusellovirus. Front. Microbiol. 2012, 3, 200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Wiedenheft, B.; Stedman, K.; Roberto, F.; Willits, D.; Gleske, A.-K.; Zoeller, L.; Snyder, J.; Douglas, T.; Young, M. Comparative Genomic Analysis of Hyperthermophilic Archaeal Fuselloviridae Viruses. J. Virol. 2004, 78, 1954. [Google Scholar] [CrossRef] [Green Version]
  27. Fulton, J.; Bothner, B.; Lawrence, M.; Johnson, J.E.; Douglas, T.; Young, M. Genetics, Biochemistry and Structure of the Archaeal Virus STIV. Biochem. Soc. Trans. 2009, 37, 114–117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Xiang, X.; Chen, L.; Huang, X.; Luo, Y.; She, Q.; Huang, L. Sulfolobus Tengchongensis Spindle-Shaped Virus STSV1: Virus-Host Interactions and Genomic Features. J. Virol. 2005, 79, 8677–8686. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Atanasova, N.S.; Roine, E.; Oren, A.; Bamford, D.H.; Oksanen, H.M. Global Network of Specific Virus-Host Interactions in Hypersaline Environments. Env. Microbiol. 2012, 14, 426–440. [Google Scholar] [CrossRef] [PubMed]
  30. Svirskaitė, J.; Oksanen, H.; Daugelavičius, R.; Bamford, D. Monitoring Physiological Changes in Haloarchaeal Cell during Virus Release. Viruses 2016, 8, 59. [Google Scholar] [CrossRef]
  31. Santos-Pérez, I.; Charro, D.; Gil-Carton, D.; Azkargorta, M.; Elortza, F.; Bamford, D.H.; Oksanen, H.M.; Abrescia, N.G.A. Structural Basis for Assembly of Vertical Single β-Barrel Viruses. Nat. Commun. 2019, 10, 1184. [Google Scholar] [CrossRef] [Green Version]
  32. Pietilä, M.K.; Atanasova, N.S.; Manole, V.; Liljeroos, L.; Butcher, S.J.; Oksanen, H.M.; Bamford, D.H. Virion Architecture Unifies Globally Distributed Pleolipoviruses Infecting Halophilic Archaea. J. Virol. 2012, 86, 5067–5079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. El Omari, K.; Li, S.; Kotecha, A.; Walter, T.S.; Bignon, E.A.; Harlos, K.; Somerharju, P.; De Haas, F.; Clare, D.K.; Molin, M.; et al. The Structure of a Prokaryotic Viral Envelope Protein Expands the Landscape of Membrane Fusion Proteins. Nat. Commun. 2019, 10, 846. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Pietilä, M.K.; Atanasova, N.S.; Oksanen, H.M.; Bamford, D.H. Modified Coat Protein Forms the Flexible Spindle-Shaped Virion of Haloarchaeal Virus His1. Environ. Microbiol. 2013, 15, 1674–1686. [Google Scholar] [CrossRef] [PubMed]
  35. Pohlschroder, M.; Schulze, S. Haloferax Volcanii. Trends Microbiol. 2019, 27, 86–87. [Google Scholar] [CrossRef]
  36. Hartman, A.L.; Norais, C.; Badger, J.H.; Delmas, S.; Haldenby, S.; Madupu, R.; Robinson, J.; Khouri, H.; Ren, Q.; Lowe, T.M.; et al. The Complete Genome Sequence of Haloferax Volcanii DS2, a Model Archaeon. PLoS ONE 2010, 5, e9605. [Google Scholar] [CrossRef] [Green Version]
  37. Allers, T.; Ngo, H.P.; Mevarech, M.; Lloyd, R.G. Development of Additional Selectable Markers for the Halophilic Archaeon Haloferax Volcanii Based on the LeuB and TrpA Genes. Appl. Environ. Microbiol. 2004, 70, 943–953. [Google Scholar] [CrossRef] [Green Version]
  38. Strillinger, E.; Grötzinger, S.W.; Allers, T.; Eppinger, J.; Weuster-Botz, D. Production of Halophilic Proteins Using Haloferax Volcanii H1895 in a Stirred-Tank Bioreactor. Appl. Microbiol. Biotechnol. 2016, 100, 1183–1195. [Google Scholar] [CrossRef] [Green Version]
  39. Stachler, A.E.; Schwarz, T.S.; Schreiber, S.; Marchfelder, A. CRISPRi as an Efficient Tool for Gene Repression in Archaea. Methods 2020, 172, 76–85. [Google Scholar] [CrossRef] [PubMed]
  40. Bisson-Filho, A.W.; Zheng, J.; Garner, E. Archaeal Imaging: Leading the Hunt for New Discoveries. Mol. Biol. Cell 2018, 29, 1675–1681. [Google Scholar] [CrossRef] [PubMed]
  41. Walsh, J.C.; Angstmann, C.N.; Bisson-Filho, A.W.; Garner, E.C.; Duggin, I.G.; Curmi, P.M.G. Division Plane Placement in Pleomorphic Archaea Is Dynamically Coupled to Cell Shape. Mol. Microbiol. 2019, 112, 785–799. [Google Scholar] [CrossRef]
  42. Duggin, I.G.; Aylett, C.H.S.; Walsh, J.C.; Michie, K.A.; Wang, Q.; Turnbull, L.; Dawson, E.M.; Harry, E.J.; Whitchurch, C.B.; Amos, L.A.; et al. CetZ Tubulin-like Proteins Control Archaeal Cell Shape. Nature 2015, 519, 362–365. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. de Silva, R.T.; Abdul-Halim, M.F.; Pittrich, D.A.; Brown, H.J.; Pohlschroder, M.; Duggin, I.G. Improved Growth and Morphological Plasticity of Haloferax Volcanii. Microbiology 2021, 167, 001012. [Google Scholar] [CrossRef] [PubMed]
  44. Kinosita, Y.; Mikami, N.; Li, Z.; Braun, F.; Quax, T.E.F.; van der Does, C.; Ishmukhametov, R.; Albers, S.V.; Berry, R.M. Motile Ghosts of the Halophilic Archaeon, Haloferax Volcanii. Proc. Natl. Acad. Sci. USA 2020, 117, 26766–26772. [Google Scholar] [CrossRef] [PubMed]
  45. Li, Z.; Rodriguez-Franco, M.; Albers, S.-V.V.; Quax, T.E.F.F. The Switch Complex ArlCDE Connects the Chemotaxis System and the Archaellum. Mol. Microbiol. 2020, 114, 468–479. [Google Scholar] [CrossRef]
  46. Quax, T.E.F.; Altegoer, F.; Rossi, F.; Li, Z.; Rodriguez-Franco, M.; Kraus, F.; Bange, G.; Albers, S.V. Structure and Function of the Archaeal Response Regulator CheY. Proc. Natl. Acad. Sci. USA 2018, 115, E1259–E1268. [Google Scholar] [CrossRef] [Green Version]
  47. Esquivel, R.N.; Schulze, S.; Xu, R.; Hippler, M.; Pohlschroder, M. Identification of Haloferax Volcanii Pilin N-Glycans with Diverse Roles in Pilus Biosynthesis, Adhesion, and Microcolony Formation. J. Biol. Chem. 2016, 291, 10602–10614. [Google Scholar] [CrossRef] [Green Version]
  48. Esquivel, R.N.; Xu, R.; Pohlschroder, M. Novel Archaeal Adhesion Pilins with a Conserved N Terminus. J. Bacteriol. 2013, 195, 3808. [Google Scholar] [CrossRef] [Green Version]
  49. Maier, L.K.; Marchfelder, A. It’s All about the T: Transcription Termination in Archaea. Biochem. Soc. Trans. 2019, 47, 461–468. [Google Scholar] [CrossRef]
  50. Humbard, M.A.; Miranda, H.v.; Lim, J.M.; Krause, D.J.; Pritz, J.R.; Zhou, G.; Chen, S.; Wells, L.; Maupin-Furlow, J.A. Ubiquitin-like Small Archaeal Modifier Proteins (SAMPs) in Haloferax Volcanii. Nature 2010, 463, 54–60. [Google Scholar] [CrossRef] [Green Version]
  51. Hawkins, M.; Malla, S.; Blythe, M.J.; Nieduszynski, C.A.; Allers, T. Accelerated Growth in the Absence of DNA Replication Origins. Nature 2013, 503, 544–547. [Google Scholar] [CrossRef]
  52. Haque, R.U.; Paradisi, F.; Allers, T. Haloferax Volcanii for Biotechnology Applications: Challenges, Current State and Perspectives. Appl. Microbiol. Biotechnol. 2020, 104, 1371–1382. [Google Scholar] [CrossRef] [Green Version]
  53. Abdul-Halim, M.F.; Schulze, S.; DiLucido, A.; Pfeiffer, F.; Filho, A.W.B.; Pohlschroder, M. Lipid Anchoring of Archaeosortase Substrates and Midcell Growth in Haloarchaea. mBio 2020, 11, e00349-20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Li, Z.; Kinosita, Y.; Rodriguez-Franco, M.; Nußbaum, P.; Braun, F.; Delpech, F.; Quax, T.E.F.; Albers, S.V. Positioning of the Motility Machinery in Halophilic Archaea. mBio 2019, 10, e00377-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Nußbaum, P.; Ithurbide, S.; Walsh, J.C.; Patro, M.; Delpech, F.; Rodriguez-Franco, M.; Curmi, P.M.G.; Duggin, I.G.; Quax, T.E.F.; Albers, S.V. An Oscillating MinD Protein Determines the Cellular Positioning of the Motility Machinery in Archaea. Curr. Biol. 2020, 30, 4956–4972.e4. [Google Scholar] [CrossRef]
  56. Jarrell, K.F.; Ding, Y.; Meyer, B.H.; Albers, S.-V.; Kaminski, L.; Eichler, J. N-Linked Glycosylation in Archaea: A Structural, Functional, and Genetic Analysis. Microbiol. Mol. Biol. Rev. 2014, 78, 304–341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Li, M.; Liu, H.; Han, J.; Liu, J.; Wang, R.; Zhao, D.; Zhou, J.; Xiang, H. Characterization of CRISPR RNA Biogenesis and Cas6 Cleavage-Mediated Inhibition of a Provirus in the Haloarchaeon Haloferax Mediterranei. J. Bacteriol. 2013, 195, 867–875. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Nuttall, S.D.; Smith, M.L.D. HF1 and HF2: Novel Bacteriophages of Halophilic Archaea. Virology 1993, 197, 678–684. [Google Scholar] [CrossRef]
  59. Mizuno, C.M.; Prajapati, B.; Lucas-Staat, S.; Sime-Ngando, T.; Forterre, P.; Bamford, D.H.; Prangishvili, D.; Krupovic, M.; Oksanen, H.M. Novel Haloarchaeal Viruses from Lake Retba Infecting Haloferax and Halorubrum Species. Environ. Microbiol. 2019, 21, 2129–2147. [Google Scholar] [CrossRef] [Green Version]
  60. Tittes, C.; Schwarzer, S.; Pfeiffer, F.; Dyall-Smith, M.; Rodriguez-Franco, M.; Oksanen, H.M.; Quax, T.E.F. Cellular and Genomic Properties of Haloferax Gibbonsii LR2-5, the Host of Euryarchaeal Virus HFTV1. Front. Microbiol. 2021, 12, 625599. [Google Scholar] [CrossRef]
  61. Dyall-Smith, M. The Halohandbook: Protocols for haloarchaeal genetics. Halohandb. Protoc. haloarchaeal Genet. 2009. Available online: https://haloarchaea.com/wp-content/uploads/2018/10/Halohandbook_2009_v7.3mds.pdf (accessed on 10 January 2022).
  62. Guindon, S.; Dufayard, J.F.; 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, 307–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Yu, G.; Smith, D.K.; Zhu, H.; Guan, Y.; Lam, T.T.Y. Ggtree: An r Package for Visualization and Annotation of Phylogenetic Trees with Their Covariates and Other Associated Data. Methods Ecol. Evol. 2017, 8, 28–36. [Google Scholar] [CrossRef]
  64. Moraru, C.; Varsani, A.; Kropinski, A.M. VIRIDIC—A Novel Tool to Calculate the Intergenomic Similarities of Prokaryote-Infecting Viruses. Viruses 2020, 12, 1268. [Google Scholar] [CrossRef]
  65. Sullivan, M.J.; Petty, N.K.; Beatson, S.A. Easyfig: A Genome Comparison Visualizer. Bioinformatics 2011, 27, 1009–1010. [Google Scholar] [CrossRef] [PubMed]
  66. Katoh, K.; Standley, D.M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  67. Nguyen, L.T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  68. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; Von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast Model Selection for Accurate Phylogenetic Estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [Green Version]
  69. Suchard, M.A.; Lemey, P.; Baele, G.; Ayres, D.L.; Drummond, A.J.; Rambaut, A. Bayesian Phylogenetic and Phylodynamic Data Integration Using BEAST 1.10. Virus Evol. 2018, 4, vey016. [Google Scholar] [CrossRef] [Green Version]
  70. Whelan, S.; Goldman, N. A General Empirical Model of Protein Evolution Derived from Multiple Protein Families Using a Maximum-Likelihood Approach. Mol. Biol. Evol. 2001, 18, 691–699. [Google Scholar] [CrossRef] [Green Version]
  71. Henikoff, S.; Henikoff, J.G. Amino Acid Substitution Matrices from Protein Blocks. Proc. Natl. Acad. Sci. USA 1992, 89, 10915. [Google Scholar] [CrossRef] [Green Version]
  72. Kingman, J.F.C. The Coalescent. Stoch. Proc. Appl. 1982, 13, 235–248. [Google Scholar] [CrossRef] [Green Version]
  73. Lemey, P.; Rambaut, A.; Drummond, A.J.; Suchard, M.A. Bayesian Phylogeography Finds Its Roots. PLOS Comput. Biol. 2009, 5, e1000520. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Minin, V.N.; Suchard, M.A. Counting Labeled Transitions in Continuous-Time Markov Models of Evolution. J. Math Biol. 2008, 56, 391–412. [Google Scholar] [CrossRef] [Green Version]
  75. Rambaut, A.; Drummond, A.J.; Xie, D.; Baele, G.; Suchard, M.A. Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7. Syst. Biol. 2018, 67, 901–904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Rodriguez-Valera, F.; Juez, G.; Kushner, D.J. Halobacterium Mediterranei Spec, Nov., a New 93 Carbohydrate-Utilizing Extreme Halophile. Syst. Appl. Microbiol. 1983, 4, 369–381. [Google Scholar] [CrossRef]
  77. Kumar, V.; Singh, B.; van Belkum, M.J.; Diep, D.B.; Chikindas, M.L.; Ermakov, A.M.; Tiwari, S.K. Halocins, Natural Antimicrobials of Archaea: Exotic or Special or Both? Biotechnol. Adv. 2021, 53, 107834. [Google Scholar] [CrossRef]
  78. Atanasova, N.S.; Pietilä, M.K.; Oksanen, H.M. Diverse Antimicrobial Interactions of Halophilic Archaea and Bacteria Extend over Geographical Distances and Cross the Domain Barrier. Microbiologyopen 2013, 2, 811–825. [Google Scholar] [CrossRef]
  79. Dyall-Smith, M.; Tang, S.L.; Russ, B.; Chiang, P.W.; Pfeiffer, F. Comparative Genomics of Two New HF1-like Haloviruses. Genes 2020, 11, 405. [Google Scholar] [CrossRef] [Green Version]
  80. Cerritelli, M.E.; Wall, J.S.; Simon, M.N.; Conway, J.F.; Steven, A.C. Stoichiometry and Domainal Organization of the Long Tail-Fiber of Bacteriophage T4: A Hinged Viral Adhesin. J. Mol. Biol. 1996, 260, 767–780. [Google Scholar] [CrossRef]
  81. Gambelli, L.; Meyer, B.H.; McLaren, M.; Sanders, K.; Quax, T.E.F.; Gold, V.A.M.; Albers, S.V.; Daum, B. Architecture and Modular Assembly of Sulfolobus S-Layers Revealed by Electron Cryotomography. Proc. Natl. Acad. Sci. USA 2019, 116, 25278–25286. [Google Scholar] [CrossRef] [Green Version]
  82. Zink, I.A.; Pfeifer, K.; Wimmer, E.; Sleytr, U.B.; Schuster, B.; Schleper, C. CRISPR-Mediated Gene Silencing Reveals Involvement of the Archaeal S-Layer in Cell Division and Virus Infection. Nat. Commun. 2019, 10, 4797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. A schematic overview of the virus–host screen. Step 1: Fresh virus stocks made from confluent or semi-confluent plates were prepared on their own host strains and the titers were determined on their own host strains. Step 2: 95 virus stocks (undiluted and 10−2 dilution) were spotted on lawns of 14 Haloferax strains. MGM medium was used as a negative control (CTL). Step 3: After incubation at 37 °C, all virus-Haloferax pairs that resulted in growth inhibition on the spot-on lawn-assay were further tested by plaque assay by making serial dilutions of the virus stock and plating with the Haloferax strains to be tested. Viral plaques observed on Haloferax were counted, the titers were determined, and positive virus-Haloferax pairs were noted in Table 1.
Figure 1. A schematic overview of the virus–host screen. Step 1: Fresh virus stocks made from confluent or semi-confluent plates were prepared on their own host strains and the titers were determined on their own host strains. Step 2: 95 virus stocks (undiluted and 10−2 dilution) were spotted on lawns of 14 Haloferax strains. MGM medium was used as a negative control (CTL). Step 3: After incubation at 37 °C, all virus-Haloferax pairs that resulted in growth inhibition on the spot-on lawn-assay were further tested by plaque assay by making serial dilutions of the virus stock and plating with the Haloferax strains to be tested. Viral plaques observed on Haloferax were counted, the titers were determined, and positive virus-Haloferax pairs were noted in Table 1.
Viruses 14 01344 g001
Figure 2. Phylogenetic tree of viruses belonging to the Haloferuviridae and Hafunaviridae families based on average nucleotide identity (ANI) values calculated using VIRIDIC software. Viruses infecting LR2-5 are surrounded by a box. Scale bar represents the number of substitutions per nucleotide position. Place of isolation: pink Senegal, purple Thailand, green Slovenia, orange Israel, and blue Italy.
Figure 2. Phylogenetic tree of viruses belonging to the Haloferuviridae and Hafunaviridae families based on average nucleotide identity (ANI) values calculated using VIRIDIC software. Viruses infecting LR2-5 are surrounded by a box. Scale bar represents the number of substitutions per nucleotide position. Place of isolation: pink Senegal, purple Thailand, green Slovenia, orange Israel, and blue Italy.
Viruses 14 01344 g002
Figure 3. (A) Isolation sites of the haloferacalesviruses, mincapviruses, and retbasiphovirus (see also Supplementary Table S2). Schemes of the viral morphologies that were observed in each group are also indicated (see also Table 1). LR2-5 infecting viruses are circled with a dash line. (B) Schematic representation of the myovirus and siphovirus virion morphologies (not in scale) and a typical genome organization consisting of different functional modules of tailed archaeal viruses.
Figure 3. (A) Isolation sites of the haloferacalesviruses, mincapviruses, and retbasiphovirus (see also Supplementary Table S2). Schemes of the viral morphologies that were observed in each group are also indicated (see also Table 1). LR2-5 infecting viruses are circled with a dash line. (B) Schematic representation of the myovirus and siphovirus virion morphologies (not in scale) and a typical genome organization consisting of different functional modules of tailed archaeal viruses.
Viruses 14 01344 g003
Figure 4. Schematic genomic alignment of the (A) mincapviruses, and (B) haloferacalesviruses. Grey bars represent homologous genomic regions. The level of nucleotide identity is reflected by the intensity of grey. Genes encoding major capsid protein (orange) and adhesins (pink) are indicated. The LR2-5 infecting viruses are circled with dashed lines. Identical or very similar genomes are shown as one and virus names are separated by /. Figures are prepared with Easyfig.
Figure 4. Schematic genomic alignment of the (A) mincapviruses, and (B) haloferacalesviruses. Grey bars represent homologous genomic regions. The level of nucleotide identity is reflected by the intensity of grey. Genes encoding major capsid protein (orange) and adhesins (pink) are indicated. The LR2-5 infecting viruses are circled with dashed lines. Identical or very similar genomes are shown as one and virus names are separated by /. Figures are prepared with Easyfig.
Viruses 14 01344 g004
Figure 5. Comparison of the adhesin and tail fiber gene sequence phylogenies to the ability of viruses to infect LR2-5. (a,c) Distribution of Markov jumps (95% HPD) for (i) real states (pink) and (ii) randomized states (grey). (b,d) Bayesian maximum clade credibility tree with discrete trait reconstruction based on adhesin or tail fiber gene sequences, respectively. Asterisk indicates node posterior probabilities higher than 0.8, and viruses in the boxes can infect LR2-5. Color code indicates the genus; blue Mincapvirus, brown Haloferacalesvirus.
Figure 5. Comparison of the adhesin and tail fiber gene sequence phylogenies to the ability of viruses to infect LR2-5. (a,c) Distribution of Markov jumps (95% HPD) for (i) real states (pink) and (ii) randomized states (grey). (b,d) Bayesian maximum clade credibility tree with discrete trait reconstruction based on adhesin or tail fiber gene sequences, respectively. Asterisk indicates node posterior probabilities higher than 0.8, and viruses in the boxes can infect LR2-5. Color code indicates the genus; blue Mincapvirus, brown Haloferacalesvirus.
Viruses 14 01344 g005
Table 1. Haloferax gibbonsii LR2-5 infecting viruses and the closely related virus isolates assigned to the same genus.
Table 1. Haloferax gibbonsii LR2-5 infecting viruses and the closely related virus isolates assigned to the same genus.
VirusVirus Morphology (5)Host Strain to Grow the VirusTiter on Own Host, pfu/mLTiter on LR2-5, pfu/mLEOP on LR2-5 (6)
Hafunaviridae (F) (1)
Haloferacalesvirus (G) (2)
HRTV-10MHalorubrum sp. B2-21.2 × 1096.8 × 1036 × 10−6
HRTV-18MHalorubrum sp. SS10-37.5 × 109-
HRTV-20MHalorubrum sp. SS10-97.0 × 1010-
HRTV-22MHalorubrum sp. SS10-91.1 × 1010-
HRTV-26MHalorubrum sp. SS10-91.4 × 1094.3 × 1063 × 10−3
HRTV-5MHalorubrum sp. s5a-32.0 × 1010-
HCTV-7 (3)MHaloarcula californiae4.8 × 1010-
[HCTV-12] (3)MHaloarcula californiae1.3 × 1010-
HCTV-9MHaloarcula californiae2.8 × 1010-
HCTV-11MHaloarcula californiae4.0 × 1010-
HRTV-9MHalorubrum sp. B2-25.3 × 109-
HRTV-16MHaloterrigena sp. SS13-7ndnd
HCTV-8MHaloarcula californiae2.3 × 10103.1 × 1051 × 10−5
HCTV-10MHalorubrum sodomense2.3 × 109-
HJTV-1MHaloarcula japonica1.6 × 109-
HRTV-13MHalorubrum sp. SS8-25.1 × 109-
HRTV-21MHalorubrum sp. SS10-92.8 × 109-
Mincapvirus (G) (2)
HSTV-2MHalorubrum sodomense8.2 × 1092.3 × 10103
HRTV-7MHalorubrum sp. B2-22.1 × 1097.9 × 1054 × 10−4
HRTV-2MHalorubrum sp. s1-22.0 × 10104.6 × 10102
HRTV-11MHalorubrum sp. SL-54.7 × 1010-
HCTV-6 (4)MHaloarcula californiae1.6 × 10103.9 × 1092 × 10−1
[HCTV-13] (4)MHaloarcula californiae2.5 × 109[1.3 × 107][5 × 10−3]
HCTV-15MHalorubrum sp. SS6-21.9 × 10109.6 × 1065 × 10−4
Haloferuviridae (F) (1)
Retbasiphovirus (G) (2)
HFTV1SHaloferax gibbonsii LR2-52.9 × 10122.9 × 10121
(1) Family; (2) Genus; (3) HCTV-7 and HCTV-12 are identical and both were tested here; HCTV-7 will be used in future; (4) HCTV-6 and HCTV-13 are identical and both were tested here; HCTV-6 will be used in future; (5) S, siphovirus morphology; M, myovirus morphology; (6) Relative efficiency of plating (EOP) on LR2-5 compared to the own host with an EOP of 1.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Aguirre Sourrouille, Z.; Schwarzer, S.; Lequime, S.; Oksanen, H.M.; Quax, T.E.F. The Viral Susceptibility of the Haloferax Species. Viruses 2022, 14, 1344. https://doi.org/10.3390/v14061344

AMA Style

Aguirre Sourrouille Z, Schwarzer S, Lequime S, Oksanen HM, Quax TEF. The Viral Susceptibility of the Haloferax Species. Viruses. 2022; 14(6):1344. https://doi.org/10.3390/v14061344

Chicago/Turabian Style

Aguirre Sourrouille, Zaloa, Sabine Schwarzer, Sebastian Lequime, Hanna M. Oksanen, and Tessa E. F. Quax. 2022. "The Viral Susceptibility of the Haloferax Species" Viruses 14, no. 6: 1344. https://doi.org/10.3390/v14061344

APA Style

Aguirre Sourrouille, Z., Schwarzer, S., Lequime, S., Oksanen, H. M., & Quax, T. E. F. (2022). The Viral Susceptibility of the Haloferax Species. Viruses, 14(6), 1344. https://doi.org/10.3390/v14061344

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop