Although next-generation sequencing (NGS) technology provides a comprehensive means with which to identify potential pathogens from clinical specimens, simple and user-friendly bioinformatics pipelines are expected to obtain the entire viral genome sequence, subsequently providing traceability, based on extensive molecular phylogenetic analyses. We have developed a web-based integrated NGS analysis tool for the viral genome (virus genome-targeted assembly pipeline: VirusTAP), which includes extensive sequence subtraction of host- or bacteria-related NGS reads prior to de novo assembly, leading to the prompt and accurate assembly of viral genome sequences from metagenomic NGS reads. The VirusTAP web site is at https://gph.niid.go.jp/cgi-bin/virustap/index.cgi/.

Next-generation sequencing technologies have allowed the rapid determination of the complete genomes of many organisms. Although shotgun sequences from large genome organisms are still difficult to reconstruct perfect contigs each of which represents a full chromosome, those from small genomes have been assembled successfully into a very small number of contigs. In this study, we show that shotgun reads from phage genomes can be reconstructed into a single contig by controlling the number of read sequences used in de novo assembly. We have developed a pipeline to assemble small viral genomes with good reliability using a resampling method from shotgun data. This pipeline, named V-GAP (Viral Genome Assembly Pipeline), will contribute to the rapid genome typing of viruses, which are highly divergent, and thus will meet the increasing need for viral genome comparisons in metagenomic studies.

Papillomaviruses infect humans and animals, most often causing benign proliferations on skin or mucosal surfaces. Rarely, these infections persist and progress to cancer. In humans, this transformation most often occurs with high-risk papillomaviruses, where viral integration is a critical event in carcinogenesis. The first aim of this study was to sequence the viral genome of canine papillomavirus (CPV) 16 from a pigmented viral plaque that progressed to metastatic squamous cell carcinoma in a dog. The second aim was to characterize multiple viral genomic deletions and translocations as well as host integration sites. The full viral genome was identified using a combination of PCR and high throughput sequencing. CPV16 is most closely related to chipapillomaviruses CPV4, CPV9, and CPV12 and we propose CPV16 be classified as a chipapillomavirus. Assembly of the full viral genome enabled identification of deletion of portions of the E1 and E2/E4 genes and two viral translocations within the squamous cell carcinoma. Genome walking was performed which identified four sites of viral integration into the host genome. This is the first description of integration of a canine papillomavirus into the host genome, raising the possibility that CPV16 may be a potential canine high-risk papillomavirus type.

We describe a new genome alignment-based model for understanding the diversity of viruses based on evolutionary genetic relationships. This approach uses information theory and a physical model to determine the information shared by the genes in two genomes. Pairwise comparisons of genes from the viruses are created from alignments using NCBI BLAST, and their match scores are combined to produce a metric between genomes, which is in turn used to determine a global classification using the 5,817 viruses on RefSeq. In cases where there is no measurable alignment between any genes, the method falls back to a coarser measure of genome relationship: the mutual information of 4-mer frequency. This results in a principled model which depends only on the genome sequence, which captures many interesting relationships between viral families, and which creates clusters which correlate well with both the Baltimore and ICTV classifications. The incremental computational cost of classifying a novel virus is low and therefore newly discovered viruses can be quickly identified and classified. The model goes beyond alignment-free classifications by producing a full phylogeny similar to those constructed by virologists using qualitative features, while relying only on objective genes. These results bolster the case for mathematical models in microbiology which can characterize organisms using only their genetic material and provide an independent check for phylogenies constructed by humans, considerably faster and more cheaply than less modern approaches.

It is well appreciated that the evolutionary divergence of genes and genomes from a common ancestor ultimately leads to incompatibilities if those genomes are hybridized. Far less is known about the ability and nature of compensatory evolution to yield the recovery of function in hybrid genomes. Here the major capsid gene of the bacteriophage T7 (40-kb dsDNA) was replaced with the homologous gene of either T3 or K11, each 22% different at the protein level from the T7 homolog. Initial fitness was moderately impaired for the T3 exchange, but the K11 exchange was not viable without a compensatory change in the T7 scaffolding protein. Subsequent adaptation of the transgenic phages led to nearly complete fitness recoveries. Compensatory changes were few, mostly in the transgene and its main interacting partner, the scaffolding protein gene. The large magnitude of fitness recovery with relatively few mutations suggests that the fitness costs of hybridizations and horizontal gene exchanges between moderately diverged genomes can potentially be short-lived through compensatory evolution.

No abstract available

Advances in cloning and sequencing technology are yielding a massive number of viral genomes. The classification and annotation of these genomes constitute important assets in the discovery of genomic variability, taxonomic characteristics and disease mechanisms. Existing classification methods are often designed for specific well-studied family of viruses. Thus, the viral comparative genomic studies could benefit from more generic, fast and accurate tools for classifying and typing newly sequenced strains of diverse virus families.

Advances in next generation sequencing make it possible to obtain high-coverage sequence data for large numbers of viral strains in a short time. However, since most bioinformatics tools are developed for command line use, the selection and accessibility of computational tools for genome assembly and variation analysis limits the ability of individual labs to perform further bioinformatics analysis.

Kaposi's sarcoma associated herpesvirus (KSHV) persists in a highly-ordered chromatin structure inside latently infected cells with the majority of the viral genome having repressive marks. However, upon reactivation the viral chromatin landscape changes into 'open' chromatin through the involvement of lysine demethylases and methyltransferases. Besides methylation of lysine residues of histone H3, arginine methylation of histone H4 plays an important role in controlling the compactness of the chromatin. Symmetric methylation of histone H4 at arginine 3 (H4R3me2s) negatively affects the methylation of histone H3 at lysine 4 (H3K4me3), an active epigenetic mark deposited on the viral chromatin during reactivation. We identified a novel binding partner to KSHV viral DNA processivity factor, ORF59-a protein arginine methyl transferase 5 (PRMT5). PRMT5 is an arginine methyltransferase that dimethylates arginine 3 (R3) of histone H4 in a symmetric manner, one hallmark of condensed chromatin. Our ChIP-seq data of symmetrically methylated H4 arginine 3 showed a significant decrease in H4R3me2s on the viral genome of reactivated cells as compared to the latent cells. Reduction in arginine methylation correlated with the binding of ORF59 on the viral chromatin and disruption of PRMT5 from its adapter protein, COPR5 (cooperator of PRMT5). Binding of PRMT5 through COPR5 is important for symmetric methylation of H4R3 and the expression of ORF59 competitively reduces the association of PRMT5 with COPR5, leading to a reduction in PRMT5 mediated arginine methylation. This ultimately resulted in a reduced level of symmetrically methylated H4R3 and increased levels of H3K4me3 marks, contributing to the formation of an open chromatin for transcription and DNA replication. Depletion of PRMT5 levels led to a decrease in symmetric methylation and increase in viral gene transcription confirming the role of PRMT5 in viral reactivation. In conclusion, ORF59 modulates histone-modifying enzymes to alter the chromatin structure during lytic reactivation.

Viral infectious diseases are significant threats to the welfare of world populations. Besides the widespread acute viral infections (e.g., dengue fever) and chronic infections [e.g., those by the human immunodeficiency virus (HIV) and hepatitis B virus (HBV)], emerging viruses, such as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), pose great challenges to the world. Genome editing technologies, including clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated (Cas) proteins, zinc-finger nucleases (ZFNs), and transcription activator-like effector nucleases (TALENs), have played essential roles in the study of new treatment for viral infectious diseases in cell lines, animal models, and clinical trials. Genome editing tools have been used to eliminate latent infections and provide resistance to new infections. Increasing evidence has shown that genome editing-based antiviral strategy is simple to design and can be quickly adapted to combat infections by a wide spectrum of viral pathogens, including the emerging coronaviruses. Here we review the development and applications of genome editing technologies for preventing or eliminating infections caused by HIV, HBV, HPV, HSV, and SARS-CoV-2, and discuss how the latest advances could enlighten further development of genome editing into a novel therapy for viral infectious diseases.

The genomes of positive-strand RNA viruses serve as a template for both protein translation and genome replication. In enteroviruses, a cloverleaf RNA structure at the 5' end of the genome functions as a switch to transition from viral translation to replication by interacting with host poly(C)-binding protein 2 (PCBP2) and the viral 3CDpro protein. We determined the structures of cloverleaf RNA from coxsackievirus and poliovirus. Cloverleaf RNA folds into an H-type four-way junction and is stabilized by a unique adenosine-cytidine-uridine (A•C-U) base triple involving the conserved pyrimidine mismatch region. The two PCBP2 binding sites are spatially proximal and are located on the opposite end from the 3CDpro binding site on cloverleaf. We determined that the A•C-U base triple restricts the flexibility of the cloverleaf stem-loops resulting in partial occlusion of the PCBP2 binding site, and elimination of the A•C-U base triple increases the binding affinity of PCBP2 to the cloverleaf RNA. Based on the cloverleaf structures and biophysical assays, we propose a new mechanistic model by which enteroviruses use the cloverleaf structure as a molecular switch to transition from viral protein translation to genome replication.

Understanding the relationships between genomic sequences is essential to the classification and characterization of living beings. The classes and characteristics of an organism can be identified in the corresponding genome space. In the genome space, the natural metric is important to describe the distribution of genomes. Therefore, the similarity of two biological sequences can be measured. Here, we report that all of the viral genomes are in 32-dimensional Euclidean space, in which the natural metric is the weighted summation of Euclidean distance of k-mer natural vectors. The classification of viral genomes in the constructed genome space further proves the convex hull principle of taxonomy, which states that convex hulls of different families are mutually disjoint. This study provides a novel geometric perspective to describe the genome sequences.

Influenza viral particles are assembled at the plasma membrane concomitantly with Rab11a-mediated endocytic transport of viral ribonucleoprotein complexes (vRNPs). The mechanism of spatiotemporal regulation of viral budozone formation and its regulatory molecules on the endocytic vesicles remain unclear. Here, we performed a proximity-based proteomics approach for Rab11a and found that ARHGAP1, a Rho GTPase-activating protein, is transported through the Rab11a-mediated apical transport of vRNP. ARHGAP1 stabilized actin filaments in infected cells for the lateral clustering of hemagglutinin (HA) molecules, a viral surface membrane protein, to the budozone. Disruption of the HA clustering results in the production of virions with low HA content, and such virions were less resistant to protease and had enhanced antigenicity, presumably because reduced clustering of viral membrane proteins exposes hidden surfaces. Collectively, these results demonstrate that Rab11a-mediated endocytic transport of ARHGAP1 with vRNPs stimulates budozone formation to ensure the integrity of virion surface required for viral survival. IMPORTANCE The endocytic transport of the influenza viral genome triggers the clustering of viral membrane proteins at the plasma membrane to form the viral budozone. However, host factors that promote viral budozone formation in concert with viral genome transport have not been identified. Here, we found that ARHGAP1, a negative regulator of the Rho family protein, is transported with the viral genome and stabilizes actin filaments to promote budozone formation. We have shown that ARHGAP1-mediated efficient formation of viral budozone was crucial for the clustering of viral HA protein to the progeny viral particles. The clustering of HA proteins on the virions is responsible for the structural integrity of the viral particles, which promotes viral stability and viral immune evasion. This study highlights the molecular mechanism that works in concert with viral genome packaging to ensure the structural integrity of viral particles.

We report the first use of whole viral genome sequencing to identify nosocomial transmission of varicella-zoster virus with fatal outcome. The index case patient, nursed in source isolation, developed disseminated zoster with rash present for 1 day before being transferred to the intensive care unit (ICU). Two patients who had received renal transplants while inpatients in an adjacent ward developed chickenpox and 1 died; neither patient had direct contact with the index patient.

Despite the accelerating number of uncultivated virus sequences discovered in metagenomics and their apparent importance for health and disease, the human gut virome and its interactions with bacteria in the gastrointestinal tract are not well understood. This is partly due to a paucity of whole-virome datasets and limitations in current approaches for identifying viral sequences in metagenomics data. Here, combining a deep-learning based metagenomics binning algorithm with paired metagenome and metavirome datasets, we develop Phages from Metagenomics Binning (PHAMB), an approach that allows the binning of thousands of viral genomes directly from bulk metagenomics data, while simultaneously enabling clustering of viral genomes into accurate taxonomic viral populations. When applied on the Human Microbiome Project 2 (HMP2) dataset, PHAMB recovered 6,077 high-quality genomes from 1,024 viral populations, and identified viral-microbial host interactions. PHAMB can be advantageously applied to existing and future metagenomes to illuminate viral ecological dynamics with other microbiome constituents.

Human cells encode up to 15 DNA polymerases with specialized functions in chromosomal DNA synthesis and damage repair. In contrast, complex DNA viruses, such as those of the herpesviridae family, encode a single B-family DNA polymerase. This disparity raises the possibility that DNA viruses may rely on host polymerases for synthesis through complex DNA geometries. We tested the importance of error-prone Y-family polymerases involved in translesion synthesis (TLS) to human cytomegalovirus (HCMV) infection. We find most Y-family polymerases involved in the nucleotide insertion and bypass of lesions restrict HCMV genome synthesis and replication. In contrast, other TLS polymerases, such as the polymerase ζ complex, which extends past lesions, was required for optimal genome synthesis and replication. Depletion of either the polζ complex or the suite of insertion polymerases demonstrate that TLS polymerases suppress the frequency of viral genome rearrangements, particularly at GC-rich sites and repeat sequences. Moreover, while distinct from HCMV, replication of the related herpes simplex virus type 1 is impacted by host TLS polymerases, suggesting a broader requirement for host polymerases for DNA virus replication. These findings reveal an unexpected role for host DNA polymerases in ensuring viral genome stability.

DUX4 is a germline transcription factor and a master regulator of zygotic genome activation. During early embryogenesis, DUX4 is crucial for maternal to zygotic transition at the 2-8-cell stage in order to overcome silencing of genes and enable transcription from the zygotic genome. In adult somatic cells, DUX4 expression is silenced and its activation in adult muscle cells causes the genetic disorder Facioscapulohumeral Muscular Dystrophy (FSHD). Here we show that herpesviruses from alpha-, beta- and gamma-herpesvirus subfamilies as well as papillomaviruses actively induce DUX4 expression to promote viral transcription and replication. We demonstrate that HSV-1 immediate early proteins directly induce expression of DUX4 and its target genes including endogenous retroelements, which mimics zygotic genome activation. We further show that DUX4 directly binds to the viral genome and promotes viral transcription. DUX4 is functionally required for herpesvirus infection, since genetic depletion of DUX4 by CRISPR/Cas9 abrogates viral replication. Our results show that herpesviruses induce DUX4 expression and its downstream germline-specific genes and retroelements, thus mimicking an early embryonic-like transcriptional program that prevents epigenetic silencing of the viral genome and facilitates herpesviral gene expression.

An increasing amount of evidence emphasizes the role of metabolic reprogramming in immune cells to fight infections. However, little is known about the regulation of metabolite transporters that facilitate and support metabolic demands. In this study, we found that the expression of equilibrative nucleoside transporter 3 (ENT3, encoded by solute carrier family 29 member 3, Slc29a3) is part of the innate immune response, which is rapidly upregulated upon pathogen invasion. The transcription of Slc29a3 is directly regulated by type I interferon-induced signaling, demonstrating that this metabolite transporter is an interferon-stimulated gene (ISG). Suprisingly, we unveil that several viruses, including SARS-CoV-2, require ENT3 to facilitate their entry into the cytoplasm. The removal or suppression of Slc29a3 expression is sufficient to significantly decrease viral replication in vitro and in vivo. Our study reveals that ENT3 is a pro-viral ISG co-opted by some viruses to gain a survival advantage.

Hepatitis C virus (HCV) is a positive single-stranded RNA virus of enormous global health importance, with direct-acting antiviral therapies replacing an immunostimulatory interferon-based regimen. The dynamics of HCV positive and negative-strand viral RNAs (vRNAs) under antiviral perturbations have not been studied at the single-cell level, leaving a gap in our understanding of antiviral kinetics and host-virus interactions. Here, we demonstrate quantitative imaging of HCV genomes in multiple infection models, and multiplexing of positive and negative strand vRNAs and host antiviral RNAs. We capture the varying kinetics with which antiviral drugs with different mechanisms of action clear HCV infection, finding the NS5A inhibitor daclatasvir to induce a rapid decline in negative-strand viral RNAs. We also find that the induction of host antiviral genes upon interferon treatment is positively correlated with viral load in single cells. This study adds smFISH to the toolbox available for analyzing the treatment of RNA virus infections.

In mammals, approximately 10% of genome sequences correspond to endogenous viral elements (EVEs), which are derived from ancient viral infections of germ cells. Although most EVEs have been inactivated, some open reading frames (ORFs) of EVEs obtained functions in the hosts. However, EVE ORFs usually remain unannotated in the genomes, and no databases are available for EVE ORFs. To investigate the function and evolution of EVEs in mammalian genomes, we developed EVE ORF databases for 20 genomes of 19 mammalian species. A total of 736,771 non-overlapping EVE ORFs were identified and archived in a database named gEVE (http://geve.med.u-tokai.ac.jp). The gEVE database provides nucleotide and amino acid sequences, genomic loci and functional annotations of EVE ORFs for all 20 genomes. In analyzing RNA-seq data with the gEVE database, we successfully identified the expressed EVE genes, suggesting that the gEVE database facilitates studies of the genomic analyses of various mammalian species.Database URL: http://geve.med.u-tokai.ac.jp.