Monitoring Insect Transposable Elements in Large Double-Stranded DNA Viruses Reveals Host-to-Virus and Virus-to-Virus Transposition

The mechanisms by which transposable elements (TEs) can be horizontally transferred between animals are unknown, but viruses are possible candidate vectors. Here, we surveyed the presence of host-derived TEs in viral genomes in 35 deep sequencing data sets produced from 11 host–virus systems, encompassing nine arthropod host species (five lepidopterans, two dipterans, and two crustaceans) and six different double-stranded (ds) DNA viruses (four baculoviruses and two iridoviruses). We found evidence of viral-borne TEs in 14 data sets, with frequencies of viral genomes carrying a TE ranging from 0.01% to 26.33% for baculoviruses and from 0.45% to 7.36% for iridoviruses. The analysis of viral populations separated by a single replication cycle revealed that viral-borne TEs originating from an initial host species can be retrieved after viral replication in another host species, sometimes at higher frequencies. Furthermore, we detected a strong increase in the number of integrations in a viral population for a TE absent from the hosts’ genomes, indicating that this TE has undergone intense transposition within the viral population. Finally, we provide evidence that many TEs found integrated in viral genomes (15/41) have been horizontally transferred in insects. Altogether, our results indicate that multiple large dsDNA viruses have the capacity to shuttle TEs in insects and they underline the potential of viruses to act as vectors of horizontal transfer of TEs. Furthermore, the finding that TEs can transpose between viral genomes of a viral species sets viruses as possible new niches in which TEs can persist and evolve.     

A sequence-independent in vitro transposon-based strategy for efficient cloning of genomes of large DNA viruses as bacterial artificial chromosomes

Bacterial artificial chromosomes (BACs) derived from genomes of large DNA viruses are powerful tools for functional delineation of viral genes. Current methods for cloning the genomes of large DNA viruses as BACs require prior knowledge of the viral sequences or the cloning of viral DNA fragments, and are tedious because of the laborious process of multiple plaque purifications, which is not feasible for some fastidious viruses. Here, we describe a novel method for cloning the genomes of large DNA viruses as BACs, which entails direct in vitro transposition of viral genomes with a BAC cassette, and subsequent recovery in Escherichia coli. Determination of insertion sites and adjacent viral sequences identify the BAC clones for genetic manipulation and functional characterization. Compared to existing methods, this new approach is highly efficient, and does not require any information on viral sequences or cloning of viral DNA fragments, and plaque purifications. This method could potentially be used for discovering previously unidentified viruses.     

Coenzymes, viruses and the RNA world

The results of a detailed bioinformatic search for ribonucleotidyl coenzyme biosynthetic sequences in DNA- and RNA viral genomes are presented. No RNA viral genome sequence available as of April 2011 appears to encode for sequences involved in coenzyme biosynthesis. In both single- and double-stranded DNA viruses a diverse array of coenzyme biosynthetic genes has been identified, but none of the viral genomes examined here encodes for a complete pathway. Although our conclusions may be constrained by the unexplored diversity of viral genomes and the biases in the construction of viral genome databases, our results do not support the possibility that RNA viruses are direct holdovers from an ancient RNA/protein world. Extrapolation of our results to evolutionary epochs prior to the emergence of DNA genomes suggest that during those early stages living entities may have depended on discontinuous genetic systems consisting of multiple small-size RNA sequences.     

Short prokaryotic DNA fragment binning using a hierarchical classifier based on linear discriminant analysis and principal component analysis

Metagenomics is an emerging field in which the power of genomic analysis is applied to an entire microbial community, bypassing the need to isolate and culture individual microbial species. Assembling of metagenomic DNA fragments is very much like the overlap-layout-consensus procedure for assembling isolated genomes, but is augmented by an additional binning step to differentiate scaffolds, contigs and unassembled reads into various taxonomic groups. In this paper, we employed n-mer oligonucleotide frequencies as the features and developed a hierarchical classifier (PCAHIER) for binning short (≤ 1,000 bps) metagenomic fragments. The principal component analysis was used to reduce the high dimensionality of the feature space. The hierarchical classifier consists of four layers of local classifiers that are implemented based on the linear discriminant analysis. These local classifiers are responsible for binning prokaryotic DNA fragments into superkingdoms, of the same superkingdom into phyla, of the same phylum into genera, and of the same genus into species, respectively. We evaluated the performance of the PCAHIER by using our own simulated data sets as well as the widely used simHC synthetic metagenome data set from the IMG/M system. The effectiveness of the PCAHIER was demonstrated through comparisons against a non-hierarchical classifier, and two existing binning algorithms (TETRA and Phylopythia).     

Engineering large viral DNA genomes using the CRISPR‐Cas9 system

Manipulation of viral genomes is essential for studying viral gene function and utilizing viruses for therapy. Several techniques for viral genome engineering have been developed. Homologous recombination in virus‐infected cells has traditionally been used to edit viral genomes; however, the frequency of the expected recombination is quite low. Alternatively, large viral genomes have been edited using a bacterial artificial chromosome (BAC) plasmid system. However, cloning of large viral genomes into BAC plasmids is both laborious and time‐consuming. In addition, because it is possible for insertion into the viral genome of drug selection markers or parts of BAC plasmids to affect viral function, artificial genes sometimes need to be removed from edited viruses. Herpes simplex virus (HSV), a common DNA virus with a genome length of 152 kbp, causes labialis, genital herpes and encephalitis. Mutant HSV is a candidate for oncotherapy, in which HSV is used to kill tumor cells. In this study, the clustered regularly interspaced short palindromic repeat‐Cas9 system was used to very efficiently engineer HSV without inserting artificial genes into viral genomes. Not only gene‐ablated HSV but also gene knock‐in HSV were generated using this method. Furthermore, selection with phenotypes of edited genes promotes the isolation efficiencies of expectedly mutated viral clones. Because our method can be applied to other DNA viruses such as Epstein–Barr virus, cytomegaloviruses, vaccinia virus and baculovirus, our system will be useful for studying various types of viruses, including clinical isolates.     

MetaPlatanus: a metagenome assembler that combines long-range sequence links and species-specific features

De novo metagenome assembly is effective in assembling multiple draft genomes, including those of uncultured organisms. However, heterogeneity in the metagenome hinders assembly and introduces interspecies misassembly deleterious for downstream analysis. For this purpose, we developed a hybrid metagenome assembler, MetaPlatanus. First, as a characteristic function, it assembles the basic contigs from accurate short reads and then iteratively utilizes long-range sequence links, species-specific sequence compositions, and coverage depth. The binning information was also used to improve contiguity. Benchmarking using mock datasets consisting of known bacteria with long reads or mate pairs revealed the high contiguity MetaPlatanus with a few interspecies misassemblies. For published human gut data with nanopore reads from potable sequencers, MetaPlatanus assembled many biologically important elements, such as coding genes, gene clusters, viral sequences, and over-half bacterial genomes. In the benchmark with published human saliva data with high-throughput nanopore reads, the superiority of MetaPlatanus was considerably more evident. We found that some high-abundance bacterial genomes were assembled only by MetaPlatanus as near-complete. Furthermore, MetaPlatanus can circumvent the limitations of highly fragmented assemblies and frequent interspecies misassembles obtained by the other tools. Overall, the study demonstrates that MetaPlatanus could be an effective approach for exploring large-scale structures in metagenomes.     

Genome-resolved viral and cellular metagenomes revealed potential key virus-host interactions in a deep freshwater lake

Metagenomics has dramatically expanded the known virosphere, but freshwater viral diversity and their ecological interaction with hosts remain poorly understood. Here, we conducted a metagenomic exploration of planktonic dsDNA prokaryotic viruses by sequencing both virion (<0.22 μm) and cellular (0.22–5.0 μm) fractions collected spatiotemporally from a deep freshwater lake (Lake Biwa, Japan). This simultaneously reconstructed 183 complete (i.e., circular) viral genomes and 57 bacterioplankton metagenome-assembled genomes. Analysis of metagenomic read coverage revealed vertical partitioning of the viral community analogous to the vertically stratified bacterioplankton community. The hypolimnetic community was generally stable during stratification, but occasionally shifted abruptly, presumably due to lysogenic induction. Genes involved in assimilatory sulfate reduction were encoded in 20 (10.9%) viral genomes, including those of dominant viruses, and may aid viral propagation in sulfur-limited freshwater systems. Hosts were predicted for 40 (21.9%) viral genomes, encompassing 10 phyla (or classes of Proteobacteria) including ubiquitous freshwater bacterioplankton lineages (e.g., Ca. Fonsibacter and Ca. Nitrosoarchaeum). Comparison with viral genomes derived from published metagenomes revealed viral phylogeographic connectivity in geographically isolated habitats. Notably, analogous to their hosts, actinobacterial viruses were among the most diverse, ubiquitous and abundant viral groups in freshwater systems, with potential high lytic activity in surface waters.     

A method for achieving complete microbial genomes and improving bins from metagenomics data

Metagenomics facilitates the study of the genetic information from uncultured microbes and complex microbial communities. Assembling complete genomes from metagenomics data is difficult because most samples have high organismal complexity and strain diversity. Some studies have attempted to extract complete bacterial, archaeal, and viral genomes and often focus on species with circular genomes so they can help confirm completeness with circularity. However, less than 100 circularized bacterial and archaeal genomes have been assembled and published from metagenomics data despite the thousands of datasets that are available. Circularized genomes are important for (1) building a reference collection as scaffolds for future assemblies, (2) providing complete gene content of a genome, (3) confirming little or no contamination of a genome, (4) studying the genomic context and synteny of genes, and (5) linking protein coding genes to ribosomal RNA genes to aid metabolic inference in 16S rRNA gene sequencing studies. We developed a semi-automated method called Jorg to help circularize small bacterial, archaeal, and viral genomes using iterative assembly, binning, and read mapping. In addition, this method exposes potential misassemblies from k-mer based assemblies. We chose species of the Candidate Phyla Radiation (CPR) to focus our initial efforts because they have small genomes and are only known to have one ribosomal RNA operon. In addition to 34 circular CPR genomes, we present one circular Margulisbacteria genome, one circular Chloroflexi genome, and two circular megaphage genomes from 19 public and published datasets. We demonstrate findings that would likely be difficult without circularizing genomes, including that ribosomal genes are likely not operonic in the majority of CPR, and that some CPR harbor diverged forms of RNase P RNA. Code and a tutorial for this method is available at https://github.com/lmlui/Jorg and is available on the DOE Systems Biology KnowledgeBase as a beta app.     

Presence of replicating virus in recombinant hepadnavirus stocks results from recombination and can be eliminated by the use of a packaging cell line

Mutant hepatitis B viruses are useful tools to study the viral life cycle and viral pathogenesis. Furthermore, recombinant hepatitis B viruses are candidate vectors for liver-directed gene therapy. Because wild-type viruses present in recombinant or mutant virus stocks may falsify experimental results and are detrimental for a viral vector, we investigated whether and to what extent wild-type virus is present in recombinant virus stocks and where it originates from. We took advantage of the duck model of hepatitis B virus infection which allows very sensitive detection of replication-competent viruses by infection of primary duck hepatocytes or of ducklings in vivo. Recombinant hepatitis B virus stocks contained significant amounts of wild-type viruses, which were most probably generated by homologous recombination between plasmids containing homologous viral sequences. In addition, replication-competent viral genomes were reconstituted from plasmids which contained replication-deficient but redundant viral sequences. Using a stable cell line for packaging of deficient viral genomes, no wild-type virus was detected, neither by infection of primary hepatocytes nor in vivo.     

Comparisons of Two Large Phaeoviral Genomes and Evolutionary Implications

The evolution of viral genomes has recently attracted considerable attention. We compare the sequences of two large viral genomes, EsV-1 and FirrV-1, belonging to the family of phaeoviruses which infect different species of marine brown algae. Although their genomes differ substantially in size, these viruses share similar morphologies and similar latent infection cycles. In fact, sequence comparisons show that the viruses have more than 60% of their genes in common. However, the order of genes is completely different in the two genomes, suggesting that extensive recombinational events in addition to several large deletions had occurred during the separate evolutionary routes from a common ancestor. We investigated genes encoding components of signal transduction pathways and genes encoding replicative functions in more detail. We found that the two genomes possess different, although overlapping, sets of genes in both classes, suggesting that different genes from each class were lost, perhaps randomly, after the separate evolution from an ancestral genome. Random loss would also account for the fact that more than one-third of the genes in one viral genome has no counterparts in the other genome. We speculate that the ancestral genome belonged to a cellular organism that had once invaded a primordial brown algal host.     

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.     

Disentangle genus microdiversity within a complex microbial community by using a multi-distance long-read binning method: example of Candidatus Accumulibacter

Complete genomes can be recovered from metagenomes by assembling and binning DNA sequences into metagenome assembled genomes (MAGs). Yet, the presence of microdiversity can hamper the assembly and binning processes, possibly yielding chimeric, highly fragmented and incomplete genomes. Here, the metagenomes of four samples of aerobic granular sludge bioreactors containing Candidatus (Ca.) Accumulibacter, a phosphate-accumulating organism of interest for wastewater treatment, were sequenced with both PacBio and Illumina. Different strategies of genome assembly and binning were investigated, including published protocols and a binning procedure adapted to the binning of long contigs (MuLoBiSC). Multiple criteria were considered to select the best strategy for Ca. Accumulibacter, whose multiple strains in every sample represent a challenging microdiversity. In this case, the best strategy relies on long-read only assembly and a custom binning procedure including MuLoBiSC in metaWRAP. Several high-quality Ca. Accumulibacter MAGs, including a novel species, were obtained independently from different samples. Comparative genomic analysis showed that MAGs retrieved in different samples harbour genomic rearrangements in addition to accumulation of point mutations. The microdiversity of Ca. Accumulibacter, likely driven by mobile genetic elements, causes major difficulties in recovering MAGs, but it is also a hallmark of the panmictic lifestyle of these bacteria.     

G-quadruplexes in genomes of viruses infecting eukaryotes or prokaryotes are under different selection pressures from hosts

G-quadruplexes in viral genomes can be applied as the targets of antiviral therapies, which has attracted wide interest. However, it is still not clear whether the pervasive number of such elements in the viral world is the result of natural selection for functionality. In this study, we identified putative quadruplex-forming sequences (PQSs) across the known viral genomes and analyzed the abundance, structural stability, and conservation of viral PQSs. A Viral Putative G-quadruplex Database (http://jsjds.hzau.edu.cn/MBPC/ViPGD/index.php/home/index) was constructed to collect the details of each viral PQS, which provides guidance for selecting the desirable PQS. The PQS with two putative G-tetrads (G₂-PQS) was significantly enriched in both eukaryotic viruses and prokaryotic viruses, whereas the PQSs with three putative G-tetrads (G₃-PQS) were only enriched in eukaryotic viruses and depleted in prokaryotic viruses. The structural stability of PQSs in prokaryotic viruses was significantly lower than that in eukaryotic viruses. Conservation analysis showed that the G₂-PQS, instead of G₃-PQS, was highly conserved within the genus. This suggested that the G₂-quadruplex might play an important role in viral biology, and the difference in the occurrence of G-quadruplex between eukaryotic viruses and prokaryotic viruses may result from the different selection pressures from hosts.     

Endogenous viral elements in animal genomes

Integration into the nuclear genome of germ line cells can lead to vertical inheritance of retroviral genes as host alleles. For other viruses, germ line integration has only rarely been documented. Nonetheless, we identified endogenous viral elements (EVEs) derived from ten non-retroviral families by systematic in silico screening of animal genomes, including the first endogenous representatives of double-stranded RNA, reverse-transcribing DNA, and segmented RNA viruses, and the first endogenous DNA viruses in mammalian genomes. Phylogenetic and genomic analysis of EVEs across multiple host species revealed novel information about the origin and evolution of diverse virus groups. Furthermore, several of the elements identified here encode intact open reading frames or are expressed as mRNA. For one element in the primate lineage, we provide statistically robust evidence for exaptation. Our findings establish that genetic material derived from all known viral genome types and replication strategies can enter the animal germ line, greatly broadening the scope of paleovirological studies and indicating a more significant evolutionary role for gene flow from virus to animal genomes than has previously been recognized.     

The genome-associated,specific RNA binding proteins of avian and mammalian type C viruses

A structural protein purified from the Rous sarcoma virus (RSV) can specifically bind in vitro to purified avian, but not mammalian, type C viral RNA. Following ultraviolet irradiation of viral particles under conditions which stabilize the polyploid 70S viral RNA, the same polypeptide can be directly purified from the RSV genome. Based on its electrophoretic mobility in polyacrylamide gels containing sodium dodecylsulfate, the RNA binding protein has been identified as the major phosphoprotein (p19) of avian type C viruses. Similar experiments show that the major phosphoproteins of mammalian type C viruses (p12 for murine viruses and p16 for endogenous primate viruses) are also the specific RNA binding proteins and, similarly, are found closely associated with the 70S RNA genomes in the intact viral particles.     

Resolving Prokaryotic Taxonomy without rRNA: Longer Oligonucleotide Word Lengths Improve Genome and Metagenome Taxonomic Classification

Oligonucleotide signatures, especially tetranucleotide signatures, have been used as method for homology binning by exploiting an organism’s inherent biases towards the use of specific oligonucleotide words. Tetranucleotide signatures have been especially useful in environmental metagenomics samples as many of these samples contain organisms from poorly classified phyla which cannot be easily identified using traditional homology methods, including NCBI BLAST. This study examines oligonucleotide signatures across 1,424 completed genomes from across the tree of life, substantially expanding upon previous work. A comprehensive analysis of mononucleotide through nonanucleotide word lengths suggests that longer word lengths substantially improve the classification of DNA fragments across a range of sizes of relevance to high throughput sequencing. We find that, at present, heptanucleotide signatures represent an optimal balance between prediction accuracy and computational time for resolving taxonomy using both genomic and metagenomic fragments. We directly compare the ability of tetranucleotide and heptanucleotide world lengths (tetranucleotide signatures are the current standard for oligonucleotide word usage analyses) for taxonomic binning of metagenome reads. We present evidence that heptanucleotide word lengths consistently provide more taxonomic resolving power, particularly in distinguishing between closely related organisms that are often present in metagenomic samples. This implies that longer oligonucleotide word lengths should replace tetranucleotide signatures for most analyses. Finally, we show that the application of longer word lengths to metagenomic datasets leads to more accurate taxonomic binning of DNA scaffolds and have the potential to substantially improve taxonomic assignment and assembly of metagenomic data.     

Phosphate transporters in marine phytoplankton and their viruses: cross-domain commonalities in viral-host gene exchanges

Phosphate (PO(4)) is an important limiting nutrient in marine environments. Marine cyanobacteria scavenge PO(4) using the high-affinity periplasmic phosphate binding protein PstS. The pstS gene has recently been identified in genomes of cyanobacterial viruses as well. Here, we analyse genes encoding transporters in genomes from viruses that infect eukaryotic phytoplankton. We identified inorganic PO(4) transporter-encoding genes from the PHO4 superfamily in several virus genomes, along with other transporter-encoding genes. Homologues of the viral pho4 genes were also identified in genome sequences from the genera that these viruses infect. Genome sequences were available from host genera of all the phytoplankton viruses analysed except the host genus Bathycoccus. Pho4 was recovered from Bathycoccus by sequencing a targeted metagenome from an uncultured Atlantic Ocean population. Phylogenetic reconstruction showed that pho4 genes from pelagophytes, haptophytes and infecting viruses were more closely related to homologues in prasinophytes than to those in what, at the species level, are considered to be closer relatives (e.g. diatoms). We also identified PHO4 superfamily members in ocean metagenomes, including new metagenomes from the Pacific Ocean. The environmental sequences grouped with pelagophytes, haptophytes, prasinophytes and viruses as well as bacteria. The analyses suggest that multiple independent pho4 gene transfer events have occurred between marine viruses and both eukaryotic and bacterial hosts. Additionally, pho4 genes were identified in available genomes from viruses that infect marine eukaryotes but not those that infect terrestrial hosts. Commonalities in marine host-virus gene exchanges indicate that manipulation of host-PO(4) uptake is an important adaptation for viral proliferation in marine systems. Our findings suggest that PO(4) -availability may not serve as a simple bottom-up control of marine phytoplankton.