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

Tailed double-stranded DNA viruses (order Caudovirales) represent the dominant morphotype among viruses infecting bacteria. Analysis and comparison of complete genome sequences of tailed bacterial viruses provided insights into their origin and evolution. Structural and genomic studies have unexpectedly revealed that tailed bacterial viruses are evolutionarily related to eukaryotic herpesviruses. Organisms from the third domain of life, Archaea, are also infected by viruses that, in their overall morphology, resemble tailed viruses of bacteria. However, high-resolution structural information is currently unavailable for any of these viruses, and only a few complete genomes have been sequenced so far. Here we identified nine proviruses that are clearly related to tailed bacterial viruses and integrated into chromosomes of species belonging to four different taxonomic orders of the Archaea. This more than doubled the number of genome sequences available for comparative studies. Our analyses indicate that highly mosaic tailed archaeal virus genomes evolve by homologous and illegitimate recombination with genomes of other viruses, by diversification, and by acquisition of cellular genes. Comparative genomics of these viruses and related proviruses revealed a set of conserved genes encoding putative proteins similar to virion assembly and maturation, as well as genome packaging proteins of tailed bacterial viruses and herpesviruses. Furthermore, fold prediction and structural modeling experiments suggest that the major capsid proteins of tailed archaeal viruses adopt the same topology as the corresponding proteins of tailed bacterial viruses and eukaryotic herpesviruses. Data presented in this study strongly support the hypothesis that tailed viruses infecting archaea share a common ancestry with tailed bacterial viruses and herpesviruses.     

Genome DNA Sequence Variation, Evolution, and Function in Bacteria and Archaea

Comparative genomics has revealed that variations in bacterial and archaeal genome DNA sequences cannot be explained by only neutral mutations. Virus resistance and plasmid distribution systems have resulted in changes in bacterial and archaeal genome sequences during evolution. The restriction-modification system, a virus resistance system, leads to avoidance of palindromic DNA sequences in genomes. Clustered, regularly interspaced, short palindromic repeats (CRISPRs) found in genomes represent yet another virus resistance system. Comparative genomics has shown that bacteria and archaea have failed to gain any DNA with GC content higher than the GC content of their chromosomes. Thus, horizontally transferred DNA regions have lower GC content than the host chromosomal DNA does. Some nucleoid-associated proteins bind DNA regions with low GC content and inhibit the expression of genes contained in those regions. This form of gene repression is another type of virus resistance system. On the other hand, bacteria and archaea have used plasmids to gain additional genes. Virus resistance systems influence plasmid distribution. Interestingly, the restriction-modification system and nucleoid-associated protein genes have been distributed via plasmids. Thus, GC content and genomic signatures do not reflect bacterial and archaeal evolutionary relationships.     

Global network of specific virus-host interactions in hypersaline environments

Hypersaline environments are dominated by archaea and bacteria and are almost entirely devoid of eukaryotic organisms. In addition, hypersaline environments contain considerable numbers of viruses. Currently, there is only a limited amount of information about these haloviruses. The ones described in detail mostly resemble head-tail bacteriophages, whereas observations based on direct microscopy of the hypersaline environmental samples highlight the abundance of non-tailed virus-like particles. Here we studied nine spatially distant hypersaline environments for the isolation of new halophilic archaea (61 isolates), halophilic bacteria (24 isolates) and their viruses (49 isolates) using a culture-dependent approach. The obtained virus isolates approximately double the number of currently described archaeal viruses. The new isolates could be divided into three tailed and two non-tailed virus morphotypes, suggesting that both types of viruses are widely distributed and characteristic for haloarchaeal viruses. We determined the sensitivity of the hosts against all isolated viruses. It appeared that the host ranges of numerous viruses extend to hosts in distant locations, supporting the idea that there is a global exchange of microbes and their viruses. It suggests that hypersaline environments worldwide function like a single habitat.     

A snapshot of viral evolution from genome analysis of the tectiviridae family

The origin, evolution and relationships of viruses are all fascinating topics. Current thinking in these areas is strongly influenced by the tailed double-stranded (ds) DNA bacteriophages. These viruses have mosaic genomes produced by genetic exchange and so new natural isolates are quite dissimilar to each other, and to laboratory strains. Consequently, they are not amenable to study by current tools for phylogenetic analysis. Less attention has been paid to the Tectiviridae family, which embraces icosahedral dsDNA bacterial viruses with an internal lipid membrane. It includes viruses, such as PRD1, that infect Gram-negative bacteria, as well as viruses like Bam35 with Gram-positive hosts. Although PRD1 and Bam35 have closely related virion morphology and genome organization, they have no detectable sequence similarity. There is strong evidence that the Bam35 coat protein has the "double-barrel trimer" arrangement of PRD1 that was first observed in adenovirus and is predicted to occur in other viruses with large facets. It is very likely that a single ancestral virus gave rise to this very large group of viruses. The unprecedented degree of conservation recently observed for two Bam35-like tectiviruses made it important to investigate those infecting Gram-negative bacteria. The DNA sequences for six PRD1-like isolates (PRD1, PR3, PR4, PR5, L17, PR772) have now been determined. Remarkably, these bacteriophages, isolated at distinctly different dates and global locations, have almost identical genomes. The discovery of almost invariant genomes for the two main Tectiviridae groups contrasts sharply with the situation in the tailed dsDNA bacteriophages. Notably, it permits a sequence analysis of the isolates revealing that the tectiviral proteins can be dissected into a slowly evolving group descended from the ancestor, the viral self, and a more rapidly changing group reflecting interactions with the host.     

Transfection of haloarchaea by the DNAs of spindle and round haloviruses and the use of transposon mutagenesis to identify non‐essential regions

Spindle‐shaped halovirus His2 and spherical halovirus SH1 represent ecologically dominant virus morphotypes in high‐salt environments. Both have linear dsDNA genomes with inverted terminal repeat sequences and terminal proteins, and probably replicate using protein priming. As a first step towards conventional genetic analyses on these viruses, we show that purified viral DNAs can transfect host cells. Intact terminal proteins were essential for this process. Despite the narrow host ranges of these viruses, at least under laboratory conditions, their DNAs were able to transfect a wide range of haloarchaeal species, demonstrating that the cytoplasms of diverse haloarchaea possess all the factors necessary for viral DNA synthesis and virion assembly. Transposon mutagenesis of viral DNAs was then used in conjunction with transfection to produce recombinant viruses, and to then map the insertion sites to identify non‐essential genes. The inserts in 34 His2 mutants were mapped precisely, and most clustered in a few, specific regions, particularly in the inverted terminal repeats and near the ends of ORFs. The results are consistent with the small genome size and densely packed, often overlapping ORFs that are transcribed as long operons. This study is the first demonstration of transfection and transposon mutagenesis in protein‐primed archaeal viruses.     

The first head-tailed virus,MFTV1, infecting hyperthermophilic methanogenic deep-sea archaea

Deep-sea hydrothermal vents are inhabited by complex communities of microbes and their viruses. Despite the importance of viruses in controlling the diversity, adaptation and evolution of their microbial hosts, to date, only eight bacterial and two archaeal viruses isolated from abyssal ecosystems have been described. Thus, our efforts focused on gaining new insights into viruses associated with deep-sea autotrophic archaea. Here, we provide the first evidence of an infection of hyperthermophilic methanogenic archaea by a head-tailed virus, Methanocaldococcus fervens tailed virus 1 (MFTV1). MFTV1 has an isometric head of 50 nm in diameter and a 150 nm-long non-contractile tail. Virions are released continuously without causing a sudden drop in host growth. MFTV1 infects Methanocaldococcus species and is the first hyperthermophilic head-tailed virus described thus far. The viral genome is a double-stranded linear DNA of 31 kb. Interestingly, our results suggest potential strategies adopted by the plasmid pMEFER01, carried by M. fervens, to spread horizontally in hyperthermophilic methanogens. The data presented here open a new window of understanding on how the abyssal mobilome interacts with hyperthermophilic marine archaea.     

An ssDNA virus infecting archaea: a new lineage of viruses with a membrane envelope

Archaeal organisms are generally known as diverse extremophiles, but they play a crucial role also in moderate environments. So far, only about 50 archaeal viruses have been described in some detail. Despite this, unusual viral morphotypes within this group have been reported. Interestingly, all isolated archaeal viruses have a double-stranded DNA (dsDNA) genome. To further characterize the diversity of archaeal viruses, we screened highly saline water samples for archaea and their viruses. Here, we describe a new haloarchaeal virus, Halorubrum pleomorphic virus 1 (HRPV-1) that was isolated from a solar saltern and infects an indigenous host belonging to the genus Halorubrum . Infection does not cause cell lysis, but slightly retards growth of the host and results in high replication of the virus. The sequenced genome (7048 nucleotides) of HRPV-1 is single-stranded DNA (ssDNA), which makes HRPV-1 the first characterized archaeal virus that does not have a dsDNA genome. In spite of this, similarities to another archaeal virus were observed. Two major structural proteins were recognized in protein analyses, and by lipid analyses it was shown that the virion contains a membrane. Electron microscopy studies indicate that the enveloped virion is pleomorphic (approximately 44 × 55 nm). HRPV-1 virion may represent commonly used virion architecture, and it seems that structure-based virus lineages may be extended to non-icosahedral viruses.     

Metagenomic analyses of novel viruses and plasmids from a cultured environmental sample of hyperthermophilic neutrophiles

Two novel viral genomes and four plasmids were assembled from an environmental sample collected from a hot spring at Yellowstone National Park, USA, and maintained anaerobically in a bioreactor at 85°C and pH 6. The double‐stranded DNA viral genomes are linear (22.7 kb) and circular (17.7 kb), and derive apparently from archaeal viruses HAV1 and HAV2. Genomic DNA was obtained from samples enriched in filamentous and tadpole‐shaped virus‐like particles respectively. They yielded few significant matches in public sequence databases reinforcing, further, the wide diversity of archaeal viruses. Several variants of HAV1 exhibit major genomic alterations, presumed to arise from viral adaptation to different hosts. They include insertions up to 350 bp, deletions up to 1.5 kb, and genes with extensively altered sequences. Some result from recombination events occurring at low complexity direct repeats distributed along the genome. In addition, a 33.8 kb archaeal plasmid pHA1 was characterized, encoding a possible conjugative apparatus, as well as three cryptic plasmids of thermophilic bacterial origin, pHB1 of 2.1 kb and two closely related variants pHB2a and pHB2b, of 5.2 and 4.8 kb respectively. Strategies are considered for assembling genomes of smaller genetic elements from complex environmental samples, and for establishing possible host identities on the basis of sequence similarity to host CRISPR immune systems.     

Viruses in soils: morphological diversity and abundance in the rhizosphere

Soil viruses are potentially of great importance as they may influence the ecology and evolution of soil biological communities through both an ability to transfer genes from host to host and as a potential cause of microbial mortality. Despite this importance, the area of soil virology is understudied. Here, we report the isolation and preliminary characterisation of viruses from soils in the Dundee area of Scotland. Different virus morphotypes including tailed, polyhedral (spherical), rod shaped, filamentous and bacilliform particles were detected in the soil samples. An apparent predominance of small spherical and filamentous bacteriophages was observed, whereas tailed bacteriophages were significantly less abundant. In this report, we also present observations and characterisation of viruses from different soil functional domains surrounding wheat roots: rhizosheath, rhizosphere and bulk soil. In spite of the differences in abundance of bacterial communities in these domains, no significant variations in viral population structure in terms of morphology and abundance were found. Typically, there were approximately 1.1–1.2 × 10⁹ virions g⁻¹ dry weight, implicating remarkable differences in virus-to-bacteria ratios in domains close to roots, rhizosphere and rhizosheath (approximately 0.27) and in bulk soil (approximately 4.68).     

Single-cell genomics-based analysis of virus–host interactions in marine surface bacterioplankton

Viral infections dynamically alter the composition and metabolic potential of marine microbial communities and the evolutionary trajectories of host populations with resulting feedback on biogeochemical cycles. It is quite possible that all microbial populations in the ocean are impacted by viral infections. Our knowledge of virus–host relationships, however, has been limited to a minute fraction of cultivated host groups. Here, we utilized single-cell sequencing to obtain genomic blueprints of viruses inside or attached to individual bacterial and archaeal cells captured in their native environment, circumventing the need for host and virus cultivation. A combination of comparative genomics, metagenomic fragment recruitment, sequence anomalies and irregularities in sequence coverage depth and genome recovery were utilized to detect viruses and to decipher modes of virus–host interactions. Members of all three tailed phage families were identified in 20 out of 58 phylogenetically and geographically diverse single amplified genomes (SAGs) of marine bacteria and archaea. At least four phage–host interactions had the characteristics of late lytic infections, all of which were found in metabolically active cells. One virus had genetic potential for lysogeny. Our findings include first known viruses of Thaumarchaeota, Marinimicrobia, Verrucomicrobia and Gammaproteobacteria clusters SAR86 and SAR92. Viruses were also found in SAGs of Alphaproteobacteria and Bacteroidetes. A high fragment recruitment of viral metagenomic reads confirmed that most of the SAG-associated viruses are abundant in the ocean. Our study demonstrates that single-cell genomics, in conjunction with sequence-based computational tools, enable in situ, cultivation-independent insights into host–virus interactions in complex microbial communities.     

Isolation and Genetic Analysis of Haloalkaliphilic Bacteriophages in a North American Soda Lake

Mono Lake is a meromictic, hypersaline, soda lake that harbors a diverse and abundant microbial community. A previous report documented the high viral abundance in Mono Lake, and pulsed-field gel electrophoresis analysis of viral DNA from lake water samples showed a diverse population based on a broad range of viral genome sizes. To better understand the ecology of bacteriophages and their hosts in this unique environment, water samples were collected between February 2001 and July 2004 for isolation of bacteriophages by using four indigenous bacterial hosts. Plaque assay results showed a differential seasonal expression of cultured bacteriophages. To reveal the diversity of uncultured bacteriophages, viral DNA from lake water samples was used to construct clone libraries. Sequence analysis of viral clones revealed homology to viral as well as bacterial proteins. Furthermore, dot blot DNA hybridization analyses showed that the uncultured viruses are more prevalent during most seasons, whereas the viral isolates (Aφ and φ2) were less prevalent, confirming the belief that uncultured viruses represent the dominant members of the community, whereas cultured isolates represent the minority species.     

Novel haloarchaeal viruses from Lake Retba infecting Haloferax and Halorubrum species

The diversity of archaeal viruses is severely undersampled compared with that of viruses infecting bacteria and eukaryotes, limiting our understanding on their evolution and environmental impacts. Here, we describe the isolation and characterization of four new viruses infecting halophilic archaea from the saline Lake Retba, located close to Dakar on the coast of Senegal. Three of the viruses, HRPV10, HRPV11 and HRPV12, have enveloped pleomorphic virions and should belong to the family Pleolipoviridae, whereas the forth virus, HFTV1, has an icosahedral capsid and a long non-contractile tail, typical of bacterial and archaeal members of the order Caudovirales. Comparative genomic and phylogenomic analyses place HRPV10, HRPV11 and HRPV12 into the genus Betapleolipovirus, whereas HFTV1 appears to be most closely related to the unclassified Halorubrum virus HRTV-4. Differently from HRTV-4, HFTV1 encodes host-derived minichromosome maintenance helicase and PCNA homologues, which are likely to orchestrate its genome replication. HFTV1, the first archaeal virus isolated on a Haloferax strain, could also infect Halorubrum sp., albeit with an eightfold lower efficiency, whereas pleolipoviruses nearly exclusively infected autochthonous Halorubrum strains. Mapping of the metagenomic sequences from this environment to the genomes of isolated haloarchaeal viruses showed that these known viruses are underrepresented in the available viromes.     

Constituents of SH1, a novel lipid-containing virus infecting the halophilic euryarchaeon Haloarcula hispanica

Recent studies have indicated that a number of bacterial and eukaryotic viruses that share a common architectural principle are related, leading to the proposal of an early common ancestor. A prediction of this model would be the discovery of similar viruses that infect archaeal hosts. Our main interest lies in icosahedral double-stranded DNA (dsDNA) viruses with an internal membrane, and we now extend our studies to include viruses infecting archaeal hosts. While the number of sequenced archaeal viruses is increasing, very little sequence similarity has been detected between bacterial and eukaryotic viruses. In this investigation we rigorously show that SH1, an icosahedral dsDNA virus infecting Haloarcula hispanica, possesses lipid structural components that are selectively acquired from the host pool. We also determined the sequence of the 31-kb SH1 genome and positively identified genes for 11 structural proteins, with putative identification of three additional proteins. The SH1 genome is unique and, except for a few open reading frames, shows no detectable similarity to other published sequences, but the overall structure of the SH1 virion and its linear genome with inverted terminal repeats is reminiscent of lipid-containing dsDNA bacteriophages like PRD1.     

Bacteriophage genomics

Comparative genomic studies of bacteriophages, especially the tailed phages, together with environmental studies, give a dramatic new picture of the size, genetic structure and dynamics of this population. Sequence comparisons reveal some of the detailed mechanisms by which these viruses evolve and influence the evolution of their bacterial and archaeal hosts. We see rampant horizontal exchange of sequences among genomes, mediated by both homologous and nonhomologous recombination. High frequency exchange among phages occupying similar ecological niches leads to a high degree of mosaic diversity in local populations. Horizontal exchange also takes place at lower frequency across the entire span of phage sequence space.     

High-Level Diversity of Tailed Phages,Eukaryote-Associated Viruses,and Virophage-Like Elements in the Metaviromes of Antarctic Soils

The metaviromes of two distinct Antarctic hyperarid desert soil communities have been characterized. Hypolithic communities, cyanobacterium-dominated assemblages situated on the ventral surfaces of quartz pebbles embedded in the desert pavement, showed higher virus diversity than surface soils, which correlated with previous bacterial community studies. Prokaryotic viruses (i.e., phages) represented the largest viral component (particularly Mycobacterium phages) in both habitats, with an identical hierarchical sequence abundance of families of tailed phages (Siphoviridae > Myoviridae > Podoviridae). No archaeal viruses were found. Unexpectedly, cyanophages were poorly represented in both metaviromes and were phylogenetically distant from currently characterized cyanophages. Putative phage genomes were assembled and showed a high level of unaffiliated genes, mostly from hypolithic viruses. Moreover, unusual gene arrangements in which eukaryotic and prokaryotic virus-derived genes were found within identical genome segments were observed. Phycodnaviridae and Mimiviridae viruses were the second-most-abundant taxa and more numerous within open soil. Novel virophage-like sequences (within the Sputnik clade) were identified. These findings highlight high-level virus diversity and novel species discovery potential within Antarctic hyperarid soils and may serve as a starting point for future studies targeting specific viral groups.     

Diverse viruses of marine archaea discovered using metagenomics

During the past decade, metagenomics became a method of choice for the discovery of novel viruses. However, host assignment for uncultured viruses remains challenging, especially for archaeal viruses, which are grossly undersampled compared to viruses of bacteria and eukaryotes. Here, we assessed the utility of CRISPR spacer targeting, tRNA gene matching and homology searches for viral signature proteins, such as major capsid proteins, for the assignment of archaeal hosts and validated these approaches on metaviromes from Yangshan Harbor (YSH). We report 35 new genomes of viruses which could be confidently assigned to hosts representing diverse lineages of marine archaea. We show that the archaeal YSH virome is highly diverse, with some viruses enriching the previously described virus groups, such as magroviruses of Marine Group II Archaea (Poseidoniales), and others representing novel groups of marine archaeal viruses. Metagenomic recruitment of Tara Oceans datasets on the YSH viral genomes demonstrated the presence of YSH Poseidoniales and Nitrososphaeria viruses in the global oceans, but also revealed the endemic YSH-specific viral lineages. Furthermore, our results highlight the relationship between the soil and marine thaumarchaeal viruses. We propose three new families within the class Caudoviricetes for the classification of the five complete viral genomes predicted to replicate in marine Poseidoniales and Nitrososphaeria, two ecologically important and widespread archaeal groups. This study illustrates the utility of viral metagenomics in exploring the archaeal virome and provides new insights into the diversity, distribution and evolution of marine archaeal viruses.     

The origins and ongoing evolution of viruses

Genome analyses of double strand DNA tailed bacteriophages argue that they evolve by recombinational reassortment of genes and by the acquisition of novel genes as simple genetic elements termed morons. These processes suggest a model for early virus evolution, wherein viruses can be regarded less as having derived from cells and more as being partners in their mutual co-evolution.     

The Vertical Distribution and Diversity of Marine Bacteriophage at a Station off Southern California

Sixty-two bacteriophages were isolated on eight indigenous bacteria from a Pacific Ocean station spanning 887-m vertical depth, on two occasions between 1999 and 2000. On the basis of 16S rRNA sequences, six hosts were tentatively identified to be in the genus Vibrio and the other two were closely related to Altermonas macleodii (W9a) and Pseudoalteromonas spp. (W13a). Restriction fragment length polymorphism (RFLP) analysis of phage genomes using AccI and HapI showed that 16 phages infecting host C4a (Vibrio) displayed 14 unique RFLP patterns. However, identical phages infecting host C4b, C6a, and C6b (all Vibrio) were obtained from both the surface layer and the hypoxic zone at 850 m. Most phage isolates from the second year had a different RFLP pattern but shared genetic similarity to the phages infecting the same host from the previous year based on a hybridization study using phage genome probes. Cluster analysis of RFLP patterns and hybridization results also indicated that phages infecting the same or genetically related hosts, in general, shared higher degrees of homology in spite of the diverse RFLP patterns. Pulsed field gel electrophoresis (PFGE) analysis of native viral genomes indicated a range in genome size from less than 40 to 200 kb, and the dominant band shifted up by about 5-10 kb in the deep samples compared to the shallow ones. Hybridization of phage genome probes with total viral community DNA from various depths suggests these isolates, or at least some of their genes, represent a detectable portion of the natural viral community and were distributed throughout the water column. Thus, the results of this study demonstrated that the genetic diversity of bacteriophage in the ocean is far greater than that of their bacterial hosts. However, host range may have contributed to the evolution of the diverse phage population in the marine environment.     

The genomes of sheeppox and goatpox viruses

Sheeppox virus (SPPV) and goatpox virus (GTPV), members of the Capripoxvirus genus of the Poxviridae, are etiologic agents of important diseases of sheep and goats in northern and central Africa, southwest and central Asia, and the Indian subcontinent. Here we report the genomic sequence and comparative analysis of five SPPV and GTPV isolates, including three pathogenic field isolates and two attenuated vaccine viruses. SPPV and GTPV genomes are approximately 150 kbp and are strikingly similar to each other, exhibiting 96% nucleotide identity over their entire length. Wild-type genomes share at least 147 putative genes, including conserved poxvirus replicative and structural genes and genes likely involved in virulence and host range. SPPV and GTPV genomes are very similar to that of lumpy skin disease virus (LSDV), sharing 97% nucleotide identity. All SPPV and GTPV genes are present in LSDV. Notably in both SPPV and GTPV genomes, nine LSDV genes with likely virulence and host range functions are disrupted, including a gene unique to LSDV (LSDV132) and genes similar to those coding for interleukin-1 receptor, myxoma virus M003.2 and M004.1 genes (two copies each), and vaccinia virus F11L, N2L, and K7L genes. The absence of these genes in SPPV and GTPV suggests a significant role for them in the bovine host range. SPPV and GTPV genomes contain specific nucleotide differences, suggesting they are phylogenetically distinct. Relatively few genomic changes in SPPV and GTPV vaccine viruses account for viral attenuation, because they contain 71 and 7 genomic changes compared to their respective field strains. Notable genetic changes include mutation or disruption of genes with predicted functions involving virulence and host range, including two ankyrin repeat proteins in SPPV and three kelch-like proteins in GTPV. These comparative genomic data indicate the close genetic relationship among capripoxviruses, and they suggest that SPPV and GTPV are distinct and likely derived from an LSDV-like ancestor.