Novel Sequence Features of Leucine-rich Repeats in Proteins from Nucleo-cytoplasmic Large DNA Viruses

: Leucine-rich repeats (LRRs) occurring in tandem are 20-29 amino acids long. Eleven LRR types have been recognized. Sequence features of LRRs from viruses were investigated using over 600 LRR proteins from 89 species. Before, metagenome data of nucleo-cytoplasmic large dsDNA viruses (NCLDVs) have been published; the 2,074 NCLDVs encode 199,021 proteins. From the NCLDVs, 549 LRR proteins were identified and analyzed. A comprehensive analysis of TpLRR and FNIP that belong to an LRR class was first performed. The repeating unit lengths (RULs) in five types are 19 residues, which are the shortest among all LRRs. Some RULs are one to five residues shorter than those of the known, corresponding LRR types. The shrinking of RUL is also observed in FNIP. The conserved hydrophobic residues, such as Leu, Val or Ile, in the consensus sequences are frequently substituted by cysteine at one or two positions. Some unique LRR types that are different from those identified previously have been observed. The present study confirms the previous result that the sequence novelty is a general feature of viral LRR proteins.

iVirus 2.0: Cyberinfrastructure-supported tools and data to power DNA virus ecology

Microbes drive myriad ecosystem processes, but under strong influence from viruses. Because studying viruses in complex systems requires different tools than those for microbes, they remain underexplored. To combat this, we previously aggregated double-stranded DNA (dsDNA) virus analysis capabilities and resources into ‘iVirus’ on the CyVerse collaborative cyberinfrastructure. Here we substantially expand iVirus’s functionality and accessibility, to iVirus 2.0, as follows. First, core iVirus apps were integrated into the Department of Energy’s Systems Biology KnowledgeBase (KBase) to provide an additional analytical platform. Second, at CyVerse, 20 software tools (apps) were upgraded or added as new tools and capabilities. Third, nearly 20-fold more sequence reads were aggregated to capture new data and environments. Finally, documentation, as “live” protocols, was updated to maximize user interaction with and contribution to infrastructure development. Together, iVirus 2.0 serves as a uniquely central and accessible analytical platform for studying how viruses, particularly dsDNA viruses, impact diverse microbial ecosystems.

A Two-tiered Functional Screen Identifies Herpesviral Transcriptional Modifiers and their Essential Domains

While traditional methods for studying large DNA viruses allow the creation of individual mutants, CRISPR/Cas9 can be used to rapidly create thousands of mutant dsDNA viruses in parallel. Here, we used this approach to study the human oncogenic Kaposi's sarcoma-associated herpesvirus (KSHV). We designed a sgRNA library containing all possible ~22,000 guides targeting the genome of KSHV - one cut site approximately every 8 base pairs - enabling the pooled screening of the entire genome. We used this tool to phenotype all possible Cas9-targeted viruses for transcription of KSHV late genes, which is required to produce structural components of the viral capsid. By performing targeted deep sequencing of the viral genome to distinguish between knock-out and in-frame alleles created by Cas9, we discovered a novel hit, ORF46 - and more specifically its DNA binding domain - is required for viral DNA replication. Our pooled Cas9 tiling screen followed by targeted deep viral sequencing represents a two-tiered screening paradigm that may be widely applicable to dsDNA viruses.

Polydnavirus Innexins Disrupt Host Cellular Encapsulation and Larval Maturation

Polydnaviruses are dsDNA viruses associated with endoparasitoid wasps. Delivery of the virus during parasitization of a caterpillar and subsequent virus gene expression is required for production of an amenable environment for parasitoid offspring development. Consequently, understanding of Polydnavirus gene function provides insight into mechanisms of host susceptibility and parasitoid wasp host range. Polydnavirus genes predominantly are arranged in multimember gene families, one of which is the vinnexins, which are virus homologues of insect gap junction genes, the innexins. Previous studies of Campoletis sonorensis Ichnovirus Vinnexins using various heterologous systems have suggested the four encoded members may provide different functionality in the infected caterpillar host. Here, we expressed two of the members, vnxG and vnxQ2, using recombinant baculoviruses in susceptible host, the caterpillar Heliothis virescens. Following intrahemocoelic injections, we observed that >90% of hemocytes (blood cells) were infected, producing recombinant protein. Larvae infected with a vinnexin-recombinant baculovirus exhibited significantly reduced molting rates relative to larvae infected with a control recombinant baculovirus and mock-infected larvae. Similarly, larvae infected with vinnexin-recombinant baculoviruses were less likely to survive relative to controls and showed reduced ability to encapsulate chromatography beads in an immune assay. In most assays, the VnxG protein was associated with more severe pathology than VnxQ2. Our findings support a role for Vinnexins in CsIV and more broadly Ichnovirus pathology in infected lepidopteran hosts, particularly in disrupting multicellular developmental and immune physiology.

Extended in vivo transcriptomes of two ascoviruses with different tissue tropisms reveal alternative mechanisms for enhancing virus reproduction in hemolymph

Ascoviruses are large dsDNA viruses characterized by the extraordinary changes they induce in cellular pathogenesis and architecture whereby after nuclear lysis and extensive hypertrophy, each cell is cleaved into numerous vesicles for virion reproduction. However, the level of viral replication and transcription in vesicles compared to other host tissues remains uncertain. Therefore, we applied RNA-Sequencing to compare the temporal transcriptome of Spodoptera frugiperda ascovirus (SfAV) and Trichoplusia ni ascovirus (TnAV) at 7, 14, and 21 days post-infection (dpi). We found most transcription occurred in viral vesicles, not in initial tissues infected, a remarkably novel reproduction mechanism compared to all other viruses and most other intracellular pathogens. Specifically, the highest level of viral gene expression occurred in hemolymph, for TnAV at 7 dpi, and SfAV at 14 dpi. Moreover, we found that host immune genes were partially down-regulated in hemolymph, where most viral replication occurred in highly dense accumulations of vesicles.

Expanding standards in viromics: in silico evaluation of dsDNA viral genome identification, classification, and auxiliary metabolic gene curation

Background Viruses influence global patterns of microbial diversity and nutrient cycles. Though viral metagenomics (viromics), specifically targeting dsDNA viruses, has been critical for revealing viral roles across diverse ecosystems, its analyses differ in many ways from those used for microbes. To date, viromics benchmarking has covered read pre-processing, assembly, relative abundance, read mapping thresholds and diversity estimation, but other steps would benefit from benchmarking and standardization. Here we use in silico-generated datasets and an extensive literature survey to evaluate and highlight how dataset composition (i.e., viromes vs bulk metagenomes) and assembly fragmentation impact (i) viral contig identification tool, (ii) virus taxonomic classification, and (iii) identification and curation of auxiliary metabolic genes (AMGs). Results The in silico benchmarking of five commonly used virus identification tools show that gene-content-based tools consistently performed well for long (≥3 kbp) contigs, while k-mer- and blast-based tools were uniquely able to detect viruses from short (≤3 kbp) contigs. Notably, however, the performance increase of k-mer- and blast-based tools for short contigs was obtained at the cost of increased false positives (sometimes up to ∼5% for virome and ∼75% bulk samples), particularly when eukaryotic or mobile genetic element sequences were included in the test datasets. For viral classification, variously sized genome fragments were assessed using gene-sharing network analytics to quantify drop-offs in taxonomic assignments, which revealed correct assignations ranging from ∼95% (whole genomes) down to ∼80% (3 kbp sized genome fragments). A similar trend was also observed for other viral classification tools such as VPF-class, ViPTree and VIRIDIC, suggesting that caution is warranted when classifying short genome fragments and not full genomes. Finally, we highlight how fragmented assemblies can lead to erroneous identification of AMGs and outline a best-practices workflow to curate candidate AMGs in viral genomes assembled from metagenomes. Conclusion Together, these benchmarking experiments and annotation guidelines should aid researchers seeking to best detect, classify, and characterize the myriad viruses ‘hidden’ in diverse sequence datasets.

Clustering and visualizing the distribution of overlapping reading frames in virus genomes

Gene overlap occurs when two or more genes are encoded by the same nucleotides. This phenomenon is found in all taxonomic domains, but is particularly common in viruses, where it may increase the information content of compact genomes or influence the creation of new genes. Here we report a global comparative study of overlapping reading frames (OvRFs) of 12,609 virus reference genomes in the NCBI database. We retrieved metadata associated with all annotated reading frames in each genome record to calculate the number, length, and frameshift of OvRFs. Our results show that while the number of OvRFs increases with genome length, they tend to be shorter in longer genomes. The majority of overlaps involve +2 frameshifts, predominantly found in dsDNA viruses. However, the longest overlaps involve no shift in reading frame (+0), increasing the selective burden of the same nucleotide positions within codons, instead of exposing additional sites to purifying selection. Next, we develop a new graph-based representation of the distribution of OvRFs among the reading frames of genomes in a given virus family. In the absence of an unambiguous partition of reading frames by homology at this taxonomic level, we used an alignment-free k-mer based approach to cluster protein coding sequences by similarity. We connect these clusters with two types of directed edges to indicate (1) that constituent reading frames are adjacent in one or more genomes, and (2) that the reading frames overlap. These adjacency graphs not only provide a natural visualization scheme, but also a novel statistical framework for analyzing the effects of gene- and genome-level attributes on the frequencies of overlaps.

Tracing dsDNA Virus–Host Coevolution through Correlation of Their G-Quadruplex-Forming Sequences

The importance of gene expression regulation in viruses based upon G-quadruplex may point to its potential utilization in therapeutic targeting. Here, we present analyses as to the occurrence of putative G-quadruplex-forming sequences (PQS) in all reference viral dsDNA genomes and evaluate their dependence on PQS occurrence in host organisms using the G4Hunter tool. PQS frequencies differ across host taxa without regard to GC content. The overlay of PQS with annotated regions reveals the localization of PQS in specific regions. While abundance in some, such as repeat regions, is shared by all groups, others are unique. There is abundance within introns of Eukaryota-infecting viruses, but depletion of PQS in introns of bacteria-infecting viruses. We reveal a significant positive correlation between PQS frequencies in dsDNA viruses and corresponding hosts from archaea, bacteria, and eukaryotes. A strong relationship between PQS in a virus and its host indicates their close coevolution and evolutionarily reciprocal mimicking of genome organization.

Primordial capsid and spooled ssDNA genome structures penetrate ancestral events of eukaryotic viruses

Marine algae viruses are important for controlling microorganism communities in the marine ecosystem, and played a fundamental role during the early events of viral evolution. Here, we have focused on one major group of marine algae viruses, the ssDNA viruses from the Bacilladnaviridae family. We present the capsid structure of the bacilladnavirus, Chaetoceros tenuissimus DNA virus type II (CtenDNAV-II), determined at 2.3 Å resolution. Structural comparison to other viruses supports the previous theory where bacilladnaviruses were proposed to have acquired their capsid protein via horizontal gene transfer from a ssRNA virus. The capsid protein contains the widespread virus jelly-roll fold, but has additional unique features; a third β-sheet and a long C-terminal tail, which are located on the capsid surface and could be important for virus transmission. Further, low-resolution reconstructions of the CtenDNAV-II genome reveal a partially spooled structure, an arrangement previously only described for dsRNA and dsDNA viruses.