Search

Nanopore direct RNA sequencing (DRS) enables the capture and full-length sequencing of native RNAs, without recoding or amplification bias. Resulting data sets may be interrogated to define the identity and location of chemically modified ribonucleotides, as well as the length of poly(A) tails, on individual RNA molecules. The success of these analyses is highly dependent on the provision of high-resolution transcriptome annotations in combination with workflows that minimize misalignments and other analysis artifacts. Existing software solutions for generating high-resolution transcriptome annotations are poorly suited to small gene-dense genomes of viruses due to the challenge of identifying distinct transcript isoforms where alternative splicing and overlapping RNAs are prevalent. To resolve this, we identified key characteristics of DRS data sets that inform resulting read alignments and developed the nanopore guided annotation of transcriptome architectures (NAGATA) software package (https://github.com/DepledgeLab/NAGATA). We demonstrate, using a combination of synthetic and original DRS data sets derived from adenoviruses, herpesviruses, coronaviruses, and human cells, that NAGATA outperforms existing transcriptome annotation software and yields a consistently high level of precision and recall when reconstructing both gene sparse and gene-dense transcriptomes. Finally, we apply NAGATA to generate the first high-resolution transcriptome annotation of the neglected pathogen human adenovirus type F41 (HAdV-41) for which we identify 77 distinct transcripts encoding at least 23 different proteins., Importance: The transcriptome of an organism denotes the full repertoire of encoded RNAs that may be expressed. This is critical to understanding the biology of an organism and for accurate transcriptomic and epitranscriptomic-based analyses. Annotating transcriptomes remains a complex task, particularly in small gene-dense organisms such as viruses which maximize their coding capacity through overlapping RNAs. To resolve this, we have developed a new software nanopore guided annotation of transcriptome architectures (NAGATA) which utilizes nanopore direct RNA sequencing (DRS) datasets to rapidly produce high-resolution transcriptome annotations for diverse viruses and other organisms., Competing Interests: The authors declare no conflict of interest.

During primary varicella zoster virus (VZV) infection, infected lymphocytes drive primary viremia, causing systemic dissemination throughout the host, including the skin. This results in cytokine expression, including interferons (IFNs), which partly limit infection. VZV also spreads from skin keratinocytes to lymphocytes prior to secondary viremia. It is not clear how VZV achieves this while evading the cytokine response. Here, we show that VZV glycoprotein C (gC) binds IFN-γ and modifies its activity, increasing the expression of a subset of IFN-stimulated genes (ISGs), including intercellular adhesion molecule 1 (ICAM1), chemokines and immunomodulatory genes. The higher ICAM1 protein level at the plasma membrane of keratinocytes facilitates lymphocyte function-associated antigen 1-dependent T cell adhesion and expression of gC during infection increases VZV spread to peripheral blood mononuclear cells. This constitutes the discovery of a strategy to modulate IFN-γ activity, upregulating a subset of ISGs, promoting enhanced lymphocyte adhesion and virus spread., (© 2024. The Author(s).)

Competing Interests: The authors declare no conflict of interest.

Herpes simplex virus 1 (HSV-1) transcription is restricted in latently infected neurons and the genomes are in mostly silenced chromatin, whereas all viral genes are transcribed in lytically infected cells, in which the genomes are dynamically chromatinized. Epigenetic regulation modulates HSV-1 transcription during lytic, latent, and reactivating infections but the precise mechanisms are not fully defined. Nucleosomes are dynamic: they slide, breathe, assemble, and disassemble. We and others have proposed that the most dynamic HSV-1 chromatin is transcriptionally competent, whereas the least dynamic is silenced. However, the mechanisms yielding the unusually dynamic viral chromatin remain unknown. Histone variants affect nucleosome dynamics. The dynamics of H2A, H2A.X, and macroH2A were enhanced in infected cells, whereas those of H2A.B were uniquely decreased. We constructed stably transduced cells expressing tagged histone H2A, H2A.B, macroH2A, or H2B, which assembles the H2A/H2B nucleosome dimers with all H2A variants. All H2A variants, as well as ectopic and endogenous H2B were assembled into HSV-1 chromatin evenly throughout the genome but canonical H2A was relatively depleted whereas H2A.B was enriched, particularly in the most dynamic viral chromatin. When viral transcription and DNA replication were restricted, H2A.B became as depleted from the viral chromatin through the entire genome as H2A. We propose that lytic HSV-1 nucleosomes are enriched in the dynamic variant H2A.B/H2B dimers to promote HSV-1 chromatin dynamics and transcriptional competency and conclude that the dynamics of HSV-1 chromatin are determined in part by the H2A variants., Importance: Herpes simplex virus 1 (HSV-1) transcription is epigenetically regulated during latent and lytic infections, and epigenetic inhibitors have been proposed as potential antiviral drugs to modulate latency and reactivation. However, the detailed epigenetic mechanisms of regulation of HSV-1 transcription have not been fully characterized and may differ from those regulating cellular transcription. Whereas lytic HSV-1 chromatin is unusually dynamic, latent silenced HSV-1 chromatin is not. The mechanisms resulting in the unique dynamics of the lytic chromatin remain unknown. Here we identify the enrichment of the highly dynamic histone 2A variant H2A in the most dynamic viral chromatin, which provides a mechanistic understanding of its unique dynamics. Future work to identify the mechanisms of enrichment in H2A.B on the viral chromatin may identify novel druggable epigenetic regulators that modulate HSV-1 latency and reactivation., Competing Interests: The authors declare no conflict of interest.

Most individuals are latently infected with herpes simplex virus type 1 (HSV-1), and it is well-established that HSV-1 establishes latency in sensory neurons of peripheral ganglia. However, it was recently proposed that latent HSV-1 is also present in immune cells recovered from the ganglia of experimentally infected mice. Here, we reanalyzed the single-cell RNA sequencing (scRNA-Seq) data that formed the basis for that conclusion. Unexpectedly, off-target priming in 3' scRNA-Seq experiments enabled the detection of non-polyadenylated HSV-1 latency-associated transcript ( LAT ) intronic RNAs. However, LAT reads were near-exclusively detected in mixed populations of cells undergoing cell death. Specific loss of HSV-1 LAT and neuronal transcripts during quality control filtering indicated widespread destruction of neurons, supporting the presence of contaminating cell-free RNA in other cells following tissue processing. In conclusion, the reported detection of latent HSV-1 in non-neuronal cells is best explained using compromised scRNA-Seq datasets.IMPORTANCEMost people are infected with herpes simplex virus type 1 (HSV-1) during their life. Once infected, the virus generally remains in a latent (silent) state, hiding within the neurons of peripheral ganglia. Periodic reactivation (reawakening) of the virus may cause fresh diseases such as cold sores. A recent study using single-cell RNA sequencing (scRNA-Seq) proposed that HSV-1 can also establish latency in the immune cells of mice, challenging existing dogma. We reanalyzed the data from that study and identified several flaws in the methodologies and analyses performed that invalidate the published conclusions. Specifically, we showed that the methodologies used resulted in widespread destruction of neurons which resulted in the presence of contaminants that confound the data analysis. We thus conclude that there remains little to no evidence for HSV-1 latency in immune cells., Competing Interests: The authors declare no conflict of interest.

The success of bacterial pathogens depends on the coordinated expression of virulence determinants. Regulatory circuits that drive pathogenesis are complex, multilayered, and incompletely understood. Here, we reveal that alterations in tRNA modifications define pathogenic phenotypes in the opportunistic pathogen Pseudomonas aeruginosa . We demonstrate that the enzymatic activity of GidA leads to the introduction of a carboxymethylaminomethyl modification in selected tRNAs. Modifications at the wobble uridine base (cmnm 5 U 34 ) of the anticodon drives translation of transcripts containing rare codons. Specifically, in P. aeruginosa the presence of GidA-dependent tRNA modifications modulates expression of genes encoding virulence regulators, leading to a cellular proteomic shift toward pathogenic and well-adapted physiological states. Our approach of profiling the consequences of chemical tRNA modifications is general in concept. It provides a paradigm of how environmentally driven tRNA modifications govern gene expression programs and regulate phenotypic outcomes responsible for bacterial adaption to challenging habitats prevailing in the host niche., Competing Interests: Competing interests statement:The authors declare no competing interest.

Unlike most RNA and DNA viruses that broadly stimulate mRNA decay and interfere with host gene expression, human cytomegalovirus (HCMV) extensively remodels the host translatome without producing an mRNA decay enzyme. By performing a targeted loss-of-function screen in primary human fibroblasts, we here identify the host CCR4-NOT deadenylase complex members CNOT1 and CNOT3 as unexpected pro-viral host factors that selectively regulate HCMV reproduction. We find that the scaffold subunit CNOT1 is specifically required for late viral gene expression and genome-wide host responses in CCR4-NOT-disrupted cells. By profiling poly(A)-tail lengths of individual HCMV and host mRNAs using nanopore direct RNA sequencing, we reveal poly(A)-tails of viral messages to be markedly longer than those of cellular mRNAs and significantly less sensitive to CCR4-NOT disruption. Our data establish that mRNA deadenylation by host CCR4-NOT is critical for productive HCMV replication and define a new mechanism whereby herpesvirus infection subverts cellular mRNA metabolism to remodel the gene expression landscape of the infected cell. Moreover, we expose an unanticipated host factor with potential to become a therapeutic anti-HCMV target., (© 2023 The Authors. Published under the terms of the CC BY 4.0 license.)

This chapter first details the structure, organization and coding content of the VZV genome to provide a foundation on which the molecular evolution of the virus can be projected. We subsequently describe the evolution of molecular profiling approaches from restriction fragment length polymorphisms to single nucleotide polymorphism profiling to modern day high-throughput sequencing approaches. We describe how the application of these methodologies led to our current model of VZV phylogeograpy including the number and structure of geographic clades and the role of recombination in reshaping these., (© 2021. The Author(s), under exclusive license to Springer Nature Switzerland AG.)

Adenovirus is a common human pathogen that relies on host cell processes for transcription and processing of viral RNA and protein production. Although adenoviral promoters, splice junctions, and polyadenylation sites have been characterized using low-throughput biochemical techniques or short read cDNA-based sequencing, these technologies do not fully capture the complexity of the adenoviral transcriptome. By combining Illumina short-read and nanopore long-read direct RNA sequencing approaches, we mapped transcription start sites and RNA cleavage and polyadenylation sites across the adenovirus genome. In addition to confirming the known canonical viral early and late RNA cassettes, our analysis of splice junctions within long RNA reads revealed an additional 35 novel viral transcripts that meet stringent criteria for expression. These RNAs include fourteen new splice junctions which lead to expression of canonical open reading frames (ORFs), six novel ORF-containing transcripts, and 15 transcripts encoding for messages that could alter protein functions through truncation or fusion of canonical ORFs. In addition, we detect RNAs that bypass canonical cleavage sites and generate potential chimeric proteins by linking distinct gene transcription units. Among these chimeric proteins we detected an evolutionarily conserved protein containing the N-terminus of E4orf6 fused to the downstream DBP/E2A ORF. Loss of this novel protein, E4orf6/DBP, was associated with aberrant viral replication center morphology and poor viral spread. Our work highlights how long-read sequencing technologies combined with mass spectrometry can reveal further complexity within viral transcriptomes and resulting proteomes., Competing Interests: The authors have declared that no competing interests exist.

Attenuation of the live varicella Oka vaccine (vOka) has been attributed to mutations in the genome acquired during cell culture passage of pOka (parent strain); however, the precise mechanisms of attenuation remain unknown. Comparative sequence analyses of several vaccine batches showed that over 100 single-nucleotide polymorphisms (SNPs) are conserved across all vaccine batches; 6 SNPs are nearly fixed, suggesting that these SNPs are responsible for attenuation. By contrast, prior analysis of chimeric vOka and pOka recombinants indicates that loci other than these six SNPs contribute to attenuation. Here, we report that pOka consists of a heterogenous population of virus sequences with two nearly equally represented bases, guanine (G) or adenine (A), at nucleotide 2096 of the ORF31 coding sequence, which encodes glycoprotein B (gB) resulting in arginine (R) or glutamine (Q), respectively, at amino acid 699 of gB. By contrast, 2096A/699Q is dominant in vOka (>99.98%). gB699Q/gH/gL showed significantly less fusion activity than gB699R/gH/gL in a cell-based fusion assay. Recombinant pOka with gB669Q (rpOka_gB699Q) had a similar growth phenotype as vOka during lytic infection in cell culture including human primary skin cells; however, rpOka_gB699R showed a growth phenotype similar to pOka. rpOka_gB699R entered neurons from axonal terminals more efficiently than rpOka_gB699Q in the presence of cell membrane-derived vesicles containing gB. Strikingly, when a mixture of pOka with both alleles equally represented was used to infect human neurons from axon terminals, pOka with gB699R was dominant for virus entry. These results identify a variant allele in gB that contributes to attenuation of vOka. IMPORTANCE The live-attenuated varicella vaccine has reduced the burden of chickenpox. Despite its development in 1974, the molecular basis for its attenuation is still not well understood. Since the live-attenuated varicella vaccine is the only licensed human herpesvirus vaccine that prevents primary disease, it is important to understand the mechanism for its attenuation. Here we identify that a variant allele in glycoprotein B (gB) selected during generation of the varicella vaccine contributes to its attenuation. This variant is impaired for fusion, virus entry into neurons from nerve terminals, and replication in human skin cells. Identification of a variant allele in gB, one of the essential herpesvirus core genes, that contributes to its attenuation may provide insights that assist in the development of other herpesvirus vaccines.

The chemical modification of ribonucleotides plays an integral role in the biology of diverse viruses and their eukaryotic host cells. Mapping the precise identity, location, and abundance of modified ribonucleotides remains a key goal of many studies aimed at characterizing the function and importance of a given modification. While mapping of specific RNA modifications through short-read sequencing approaches has powered a wealth of new discoveries in the past decade, this approach is limited by inherent biases and an absence of linkage information. Moreover, in viral contexts, the challenge is increased due to the compact nature of viral genomes giving rise to many overlapping transcript isoforms that cannot be adequately resolved using short-read sequencing approaches. The recent emergence of nanopore sequencing, specifically the ability to directly sequence native RNAs from virus-infected host cells, provides not just a new methodology for mapping modified ribonucleotides but also a new conceptual framework for what can be derived from the resulting sequencing data. In this minireview, we provide a detailed overview of how nanopore direct RNA sequencing works, the computational approaches applied to identify modified ribonucleotides, and the core concepts underlying both. We further highlight recent studies that have applied this approach to interrogating viral biology and finish by discussing key experimental considerations and how we predict that these methodologies will continue to evolve.

Understanding the molecular mechanisms of herpes simplex virus 1 (HSV-1) latent infection and reactivation in neurons requires the use of in vitro model systems. Establishing a quiescent infection in cultured neurons is problematic, as any infectious virus released can superinfect the cultures. Previous studies have used the viral DNA replication inhibitor acyclovir to prevent superinfection and promote latency establishment. Data from these previous models have shown that reactivation is biphasic, with an initial phase I expression of all classes of lytic genes, which occurs independently of histone demethylase activity and viral DNA replication but is dependent on the cell stress protein DLK. Here, we describe a new model system using HSV-1 Stayput-GFP, a reporter virus that is defective for cell-to-cell spread and establishes latent infections without the need for acyclovir. The establishment of a latent state requires a longer time frame than previous models using DNA replication inhibitors. This results in a decreased ability of the virus to reactivate using established inducers, and as such, a combination of reactivation triggers is required. Using this system, we demonstrate that biphasic reactivation occurs even when latency is established in the absence of acyclovir. Importantly, phase I lytic gene expression still occurs in a histone demethylase and viral DNA replication-independent manner and requires DLK activity. These data demonstrate that the two waves of viral gene expression following HSV-1 reactivation are independent of secondary infection and not unique to systems that require acyclovir to promote latency establishment. IMPORTANCE Herpes simplex virus-1 (HSV-1) enters a latent infection in neurons and periodically reactivates. Reactivation manifests as a variety of clinical symptoms. Studying latency and reactivation in vitro is invaluable, allowing the molecular mechanisms behind both processes to be targeted by therapeutics that reduce the clinical consequences. Here, we describe a novel in vitro model system using a cell-to-cell spread-defective HSV-1, known as Stayput-GFP, which allows for the study of latency and reactivation at the single neuron level. We anticipate this new model system will be an incredibly valuable tool for studying the establishment and reactivation of HSV-1 latent infection in vitro . Using this model, we find that initial reactivation events are dependent on cellular stress kinase DLK but independent of histone demethylase activity and viral DNA replication. Our data therefore further validate the essential role of DLK in mediating a wave of lytic gene expression unique to reactivation.

Infection of immunosuppressed transplant patients with the human γ-herpesvirus Epstein-Barr virus (EBV) is associated with post-transplant lymphoproliferative disease (PTLD), an often fatal complication. Immunosuppressed miniature pigs infected with γ-herpesvirus porcine lymphotropic herpesvirus 1 (PLHV1) develop a similar disease, identifying pigs as a potential preclinical model for PTLD in humans. BILF1 is a G protein-coupled receptor (GPCR) encoded by EBV with constitutive activity linked to tumorigenesis and immunoevasive function downregulating MHC-I. In the present study, we compared BILF1-orthologues encoded by the three known PLHVs (PLHV1-3) with EBV-BILF1 to determine pharmacological suitability of BILF1 orthologues as model system to study EBV-BILF1 druggability. Cell surface localization, constitutive internalization, and MHC-I downregulation as well as membrane proximal constitutive Gα i signaling patterns were conserved across all BILFs. Only subtle differences between the individual BILFs were observed in downstream transcription factor activation. Using Illumina sequencing, PLHV1 was observed in lymphatic tissue from PTLD-diseased, but not non-diseased pigs. Importantly, these tissues showed enhanced expression of PLHV1-BILF1 supporting its involvement in PTLD infection., Competing Interests: Author MJ was employed, in part, by the company the Vaccine Group, Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Mavri, Kubale, Depledge, Zuo, Huang, Breuer, Vrecl, Jarvis, Jovičić, Petan, Ehlers, Rosenkilde and Spiess.)

Motivation: The chemical modification of ribonucleotides regulates the structure, stability and interactions of RNAs. Profiling of these modifications using short-read (Illumina) sequencing techniques provides high sensitivity but low-to-medium resolution i.e. modifications cannot be assigned to specific transcript isoforms in regions of sequence overlap. An alternative strategy uses current fluctuations in nanopore-based long read direct RNA sequencing (DRS) to infer the location and identity of nucleotides that differ between two experimental conditions. While highly sensitive, these signal-level analyses require high-quality transcriptome annotations and thus are best suited to the study of model organisms. By contrast, the detection of RNA modifications in microbial organisms which typically have no or low-quality annotations requires an alternative strategy. Here, we demonstrate that signal fluctuations directly influence error rates during base-calling and thus provides an alternative approach for identifying modified nucleotides., Results: DRUMMER (Detection of Ribonucleic acid Modifications Manifested in Error Rates) (i) utilizes a range of statistical tests and background noise correction to identify modified nucleotides with high confidence, (ii) operates with similar sensitivity to signal-level analysis approaches and (iii) correlates very well with orthogonal approaches. Using well-characterized DRS datasets supported by independent meRIP-Seq and miCLIP-Seq datasets we demonstrate that DRUMMER operates with high sensitivity and specificity., Availability and Implementation: DRUMMER is written in Python 3 and is available as open source in the GitHub repository: https://github.com/DepledgeLab/DRUMMER., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2022. Published by Oxford University Press.)

Regulated loading of eIF3-bound 40S ribosomes on capped mRNA is generally dependent upon the translation initiation factor eIF4E; however, mRNA translation often proceeds during physiological stress, such as virus infection, when eIF4E availability and activity are limiting. It remains poorly understood how translation of virus and host mRNAs are regulated during infection stress. While initially sensitive to mTOR inhibition, which limits eIF4E-dependent translation, we show that protein synthesis in human cytomegalovirus (HCMV)-infected cells unexpectedly becomes progressively reliant upon eIF3d. Targeting eIF3d selectively inhibits HCMV replication, reduces polyribosome abundance, and interferes with expression of essential virus genes and a host gene expression signature indicative of chronic ER stress that fosters HCMV reproduction. This reveals a strategy whereby cellular eIF3d-dependent protein production is hijacked to exploit virus-induced ER stress. Moreover, it establishes how switching between eIF4E and eIF3d-responsive cap-dependent translation can differentially tune virus and host gene expression in infected cells., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

Single-cell RNA sequencing (scRNA-seq) is a powerful technique for dissecting the complexity of normal and diseased tissues, enabling characterization of cell diversity and heterogeneous phenotypic states in unprecedented detail. However, this technology has been underutilized for exploring the interactions between the host cell and viral pathogens in latently infected cells. Herein, we use scRNA-seq and single-molecule sensitivity fluorescent in situ hybridization (smFISH) technologies to investigate host single-cell transcriptome changes upon the reactivation of a human neurotropic virus, herpes simplex virus-1 (HSV-1). We identify the stress sensor growth arrest and DNA damage-inducible 45 beta (Gadd45b) as a critical antiviral host factor that regulates HSV-1 reactivation events in a subpopulation of latently infected primary neurons. We show that distinct subcellular localization of Gadd45b correlates with the viral late gene expression program, as well as the expression of the viral transcription factor, ICP4. We propose that a hallmark of a "successful" or "aborted" HSV-1 reactivation state in primary neurons is determined by a unique subcellular localization signature of the stress sensor Gadd45b., (© 2021 The Authors.)

Primary infection with varicella-zoster virus (VZV) causes varicella and the establishment of lifelong latency in sensory ganglion neurons. In one-third of infected individuals VZV reactivates from latency to cause herpes zoster, often complicated by difficult-to-treat chronic pain. Experimental infection of non-human primates with simian varicella virus (SVV) recapitulates most features of human VZV disease, thereby providing the opportunity to study the pathogenesis of varicella and herpes zoster in vivo. However, compared to VZV, the transcriptome and the full coding potential of SVV remains incompletely understood. Here, we performed nanopore direct RNA sequencing to annotate the SVV transcriptome in lytically SVV-infected African green monkey (AGM) and rhesus macaque (RM) kidney epithelial cells. We refined structures of canonical SVV transcripts and uncovered numerous RNA isoforms, splicing events, fusion transcripts and non-coding RNAs, mostly unique to SVV. We verified the expression of canonical and newly identified SVV transcripts in vivo, using lung samples from acutely SVV-infected cynomolgus macaques. Expression of selected transcript isoforms, including those located in the unique left-end of the SVV genome, was confirmed by reverse transcription PCR. Finally, we performed detailed characterization of the SVV homologue of the VZV latency-associated transcript (VLT), located antisense to ORF61. Analogous to VZV VLT, SVV VLT is multiply spliced and numerous isoforms are generated using alternative transcription start sites and extensive splicing. Conversely, low level expression of a single spliced SVV VLT isoform defines in vivo latency. Notably, the genomic location of VLT core exons is highly conserved between SVV and VZV. This work thus highlights the complexity of lytic SVV gene expression and provides new insights into the molecular biology underlying lytic and latent SVV infection. The identification of the SVV VLT homolog further underlines the value of the SVV non-human primate model to develop new strategies for prevention of herpes zoster., Competing Interests: The authors have declared that no competing interests exist.

Primary varicella-zoster virus (VZV) infection leads to varicella and the establishment of lifelong latency in sensory ganglion neurons. Reactivation of latent VZV causes herpes zoster, which is frequently associated with chronic pain. Latent viral gene expression is restricted to the VZV latency-associated transcript (VLT) and VLT-ORF63 (VLT63) fusion transcripts. Since VLT and VLT63 encode proteins that are expressed during lytic infection, we investigated whether pVLT and pVLT-ORF63 are essential for VZV replication by performing VZV genome mutagenesis using CRISPR/Cas9 and BAC technologies. We first established that CRISPR/Cas9 can efficiently mutate VZV genomes in lytically VZV-infected cells through targeting non-essential genes ORF8 and ORF11 and subsequently show recovery of viable mutant viruses. By contrast, the VLT region was markedly resistant to CRISPR/Cas9 editing. Whereas most mutants expressed wild-type or N-terminally altered versions of pVLT and pVLT-ORF63, only a minority of the resulting mutant viruses lacked pVLT and pVLT-ORF63 coding potential. Growth curve analysis showed that pVLT/pVLT-ORF63 negative viruses were viable, but impaired in growth in epithelial cells. We confirmed this phenotype independently using BAC-derived pVLT/pVLT-ORF63 negative and repaired viruses. Collectively, these data demonstrate that pVLT and/or pVLT-ORF63 are dispensable for lytic VZV replication but promote efficient VZV infection in epithelial cells.

Members of the family Herpesviridae have enveloped, spherical virions with characteristic complex structures consisting of symmetrical and non-symmetrical components. The linear, double-stranded DNA genomes of 125-241 kbp contain 70-170 genes, of which 43 have been inherited from an ancestral herpesvirus. In general, herpesviruses have coevolved with and are highly adapted to their hosts, which comprise many mammalian, avian and reptilian species. Following primary infection, they are able to establish lifelong latent infection, during which there is limited viral gene expression. Severe disease is usually observed only in the foetus, the very young, the immunocompromised or following infection of an alternative host. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Herpesviridae , which is available at ictv.global/report/herpesviridae.

N 6 -methyladenosine (m 6 A) is the most abundant internal messenger RNA (mRNA) modification, contributing to the processing, stability, and function of methylated RNAs. Methylation occurs in the nucleus during pre-mRNA synthesis and requires a core methyltransferase complex consisting of METTL3, METTL14, and WTAP. During herpes simplex virus (HSV-1) infection, cellular gene expression is profoundly suppressed, allowing the virus to monopolize the host transcription and translation apparatus and antagonize antiviral responses. The extent to which HSV-1 uses or manipulates the m 6 A pathway is not known. Here, we show that, in primary fibroblasts, HSV-1 orchestrates a striking redistribution of the nuclear m 6 A machinery that progresses through the infection cycle. METTL3 and METTL14 are dispersed into the cytoplasm, whereas WTAP remains nuclear. Other regulatory subunits of the methyltransferase complex, along with the nuclear m 6 A-modified RNA binding protein YTHDC1 and nuclear demethylase ALKBH5, are similarly redistributed. These changes require ICP27, a viral regulator of host mRNA processing that mediates the nucleocytoplasmic export of viral late mRNAs. Viral gene expression is initially reduced by small interfering RNA (siRNA)-mediated inactivation of the m 6 A methyltransferase but becomes less impacted as the infection advances. Redistribution of the nuclear m 6 A machinery is accompanied by a wide-scale reduction in the installation of m 6 A and other RNA modifications on both host and viral mRNAs. These results reveal a far-reaching mechanism by which HSV-1 subverts host gene expression to favor viral replication., Competing Interests: The authors declare no competing interest.

N 6 -methyladenosine (m 6 A) is an abundant internal RNA modification, influencing transcript fate and function in uninfected and virus-infected cells. Installation of m 6 A by the nuclear RNA methyltransferase METTL3 occurs cotranscriptionally; however, the genomes of some cytoplasmic RNA viruses are also m 6 A-modified. How the cellular m 6 A modification machinery impacts coronavirus replication, which occurs exclusively in the cytoplasm, is unknown. Here we show that replication of SARS-CoV-2, the agent responsible for the COVID-19 pandemic, and a seasonal human β-coronavirus HCoV-OC43, can be suppressed by depletion of METTL3 or cytoplasmic m 6 A reader proteins YTHDF1 and YTHDF3 and by a highly specific small molecule METTL3 inhibitor. Reduction of infectious titer correlates with decreased synthesis of viral RNAs and the essential nucleocapsid (N) protein. Sites of m 6 A modification on genomic and subgenomic RNAs of both viruses were mapped by methylated RNA immunoprecipitation sequencing (meRIP-seq). Levels of host factors involved in m 6 A installation, removal, and recognition were unchanged by HCoV-OC43 infection; however, nuclear localization of METTL3 and cytoplasmic m 6 A readers YTHDF1 and YTHDF2 increased. This establishes that coronavirus RNAs are m 6 A-modified and host m 6 A pathway components control β-coronavirus replication. Moreover, it illustrates the therapeutic potential of targeting the m 6 A pathway to restrict coronavirus reproduction., (© 2021 Burgess et al.; Published by Cold Spring Harbor Laboratory Press.)

The contributions of the viral component of the microbiome-the virome-to the development of innate and adaptive immunity are largely unknown. Here, we systematically defined the host response in mice to a panel of eukaryotic enteric viruses representing six different families. Infections with most of these viruses were asymptomatic in the mice, the magnitude and duration of which was dependent on the microbiota. Flow cytometric and transcriptional profiling of mice mono-associated with these viruses unveiled general adaptations by the host, such as lymphocyte differentiation and IL-22 signatures in the intestine, as well as numerous viral-strain-specific responses that persisted. Comparison with a dataset derived from analogous bacterial mono-association in mice identified bacterial species that evoke an immune response comparable with the viruses we examined. These results expand an understanding of the immune space occupied by the enteric virome and underscore the importance of viral exposure events., Competing Interests: Declaration of interests K.C. has received research support from Pfizer, Takeda, Pacific Biosciences, and Abbvie. K.C. has also consulted for or received an honorarium from Puretech Health, Genentech, and Abbvie. K.C. holds U.S. patent 10,722,600 and provisional patent 62/935,035., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Epstein-Barr virus (EBV) is one of the most common viruses latently infecting humans. Little is known about the impact of human genetic variation on the large inter-individual differences observed in response to EBV infection. To search for a potential imprint of host genomic variation on the EBV sequence, we jointly analyzed paired viral and human genomic data from 268 HIV-coinfected individuals with CD4 + T cell count < 200/mm 3 and elevated EBV viremia. We hypothesized that the reactivated virus circulating in these patients could carry sequence variants acquired during primary EBV infection, thereby providing a snapshot of early adaptation to the pressure exerted on EBV by the individual immune response. We searched for associations between host and pathogen genetic variants, taking into account human and EBV population structure. Our analyses revealed significant associations between human and EBV sequence variation. Three polymorphic regions in the human genome were found to be associated with EBV variation: one at the amino acid level (BRLF1:p.Lys316Glu); and two at the gene level (burden testing of rare variants in BALF5 and BBRF1). Our findings confirm that jointly analyzing host and pathogen genomes can identify sites of genomic interactions, which could help dissect pathogenic mechanisms and suggest new therapeutic avenues.

Varicella-zoster virus (VZV) establishes lifelong neuronal latency in most humans world-wide, reactivating in one-third to cause herpes zoster and occasionally chronic pain. How VZV establishes, maintains and reactivates from latency is largely unknown. VZV transcription during latency is restricted to the latency-associated transcript (VLT) and RNA 63 (encoding ORF63) in naturally VZV-infected human trigeminal ganglia (TG). While significantly more abundant, VLT levels positively correlated with RNA 63 suggesting co-regulated transcription during latency. Here, we identify VLT-ORF63 fusion transcripts and confirm VLT-ORF63, but not RNA 63, expression in human TG neurons. During in vitro latency, VLT is transcribed, whereas VLT-ORF63 expression is induced by reactivation stimuli. One isoform of VLT-ORF63, encoding a fusion protein combining VLT and ORF63 proteins, induces broad viral gene transcription. Collectively, our findings show that VZV expresses a unique set of VLT-ORF63 transcripts, potentially involved in the transition from latency to lytic VZV infection.

Adenovirus is a nuclear replicating DNA virus reliant on host RNA processing machinery. Processing and metabolism of cellular RNAs can be regulated by METTL3, which catalyzes the addition of N6-methyladenosine (m 6 A) to mRNAs. While m 6 A-modified adenoviral RNAs have been previously detected, the location and function of this mark within the infectious cycle is unknown. Since the complex adenovirus transcriptome includes overlapping spliced units that would impede accurate m 6 A mapping using short-read sequencing, here we profile m 6 A within the adenovirus transcriptome using a combination of meRIP-seq and direct RNA long-read sequencing to yield both nucleotide and transcript-resolved m 6 A detection. Although both early and late viral transcripts contain m 6 A, depletion of m 6 A writer METTL3 specifically impacts viral late transcripts by reducing their splicing efficiency. These data showcase a new technique for m 6 A discovery within individual transcripts at nucleotide resolution, and highlight the role of m 6 A in regulating splicing of a viral pathogen.

Varicella-zoster virus (VZV), a double-stranded DNA virus, causes varicella, establishes lifelong latency in ganglionic neurons, and reactivates later in life to cause herpes zoster, commonly associated with chronic pain. The VZV genome is densely packed and produces multitudes of overlapping transcripts deriving from both strands. While 71 distinct open reading frames (ORFs) have thus far been experimentally defined, the full coding potential of VZV remains unknown. Here, we integrated multiple short-read RNA sequencing approaches with long-read direct RNA sequencing on RNA isolated from VZV-infected cells to provide a comprehensive reannotation of the lytic VZV transcriptome architecture. Through precise mapping of transcription start sites, splice junctions, and polyadenylation sites, we identified 136 distinct polyadenylated VZV RNAs that encode canonical ORFs, noncanonical ORFs, and ORF fusions, as well as putative noncoding RNAs (ncRNAs). Furthermore, we determined the kinetic class of all VZV transcripts and observed, unexpectedly, that transcripts encoding the ORF62 protein, previously designated Immediate-Early , were expressed with Late kinetics. Our work showcases the complexity of the VZV transcriptome and provides a comprehensive resource that will facilitate future functional studies of coding RNAs, ncRNAs, and the biological mechanisms underlying the regulation of viral transcription and translation during lytic VZV infection. IMPORTANCE Transcription from herpesviral genomes, executed by the host RNA polymerase II and regulated by viral proteins, results in coordinated viral gene expression to efficiently produce infectious progeny. However, the complete coding potential and regulation of viral gene expression remain ill-defined for the human alphaherpesvirus varicella-zoster virus (VZV), causative agent of both varicella and herpes zoster. Here, we present a comprehensive overview of the VZV transcriptome and the kinetic class of all identified viral transcripts, using two virus strains and two biologically relevant cell types. Additionally, our data provide an overview of how VZV diversifies its transcription from one of the smallest herpesviral genomes. Unexpectedly, the transcript encoding the major viral transactivator protein (pORF62) was expressed with Late kinetics, whereas orthologous transcripts in other alphaherpesviruses are typically expressed during the immediate early phase. Therefore, our work both establishes the architecture of the VZV transcriptome and provides insight into regulation of alphaherpesvirus gene expression., (Copyright © 2020 Braspenning et al.)

The varicella-zoster virus (VZV) genome, comprises both unique and repeated regions. The genome also includes reiteration regions, designated R1 to R5, which are tandemly repeating sequences termed elements. These regions represent an understudied feature of the VZV genome. The R4 region is duplicated, with one copy in the internal repeat short (IRs) which we designated R4A and a second copy in the terminal repeat short (TRs) termed R4B. We developed primers to amplify and Sanger sequence these regions, including independent amplification of both R4 regions. Reiteration regions from >80 cases of PCR-confirmed shingles were sequenced and analyzed. Complete genome sequences for the remaining portions of these viruses were determined using Illumina MiSeq. We identified 28 elements not previously reported, including at least one element for each R region. Length heterogeneity was substantial in R3, R4A and R4B. Length heterogeneity between the two copies of R4 was common., Competing Interests: Declaration of competing interest The authors report no conflicts of interest related to this work., (Copyright © 2020. Published by Elsevier Inc.)

The genomes of DNA viruses encode deceptively complex transcriptomes evolved to maximize coding potential within the confines of a relatively small genome. Defining the full range of viral RNAs produced during an infection is key to understanding the viral replication cycle and its interactions with the host cell. Traditional short-read (Illumina) sequencing approaches are problematic in this setting due to the difficulty of assigning short reads to individual RNAs in regions of transcript overlap and to the biases introduced by the required recoding and amplification steps. Additionally, different methodologies may be required to analyze the 5' and 3' ends of RNAs, which increases both cost and effort. The advent of long-read nanopore sequencing simplifies this approach by providing a single assay that captures and sequences full length RNAs, either in cDNA or native RNA form. The latter is particularly appealing as it reduces known recoding biases whilst allowing more advanced analyses such as estimation of poly(A) tail length and the detection of RNA modifications including N 6 -methyladenosine. Using herpes simplex virus (HSV)-infected primary fibroblasts as a template, we provide a step-by-step guide to the production of direct RNA sequencing libraries suitable for sequencing using Oxford Nanopore Technologies platforms and provide a simple computational approach to deriving a high-quality annotation of the HSV transcriptome from the resulting sequencing data. © 2020 by John Wiley & Sons, Inc. Basic Protocol 1: Productive infection of primary fibroblasts with herpes simplex virus Support Protocol: Cell passage and plating of primary fibroblasts Basic Protocol 2: Preparation and sequencing of dRNA-seq libraries from virus-infected cells Basic Protocol 3: Processing, alignment, and analysis of dRNA-seq datasets., (© 2020 John Wiley & Sons, Inc.)

Varicella zoster virus (VZV) is a skin-tropic virus that infects epidermal keratinocytes and causes chickenpox. Although common, VZV infection can be life-threatening, particularly in the immunocompromized. Therefore, understanding VZV-keratinocyte interactions is important to find new treatments beyond vaccination and antiviral drugs. In VZV-infected skin, kallikrein 6 and the ubiquitin ligase MDM2 are upregulated concomitant with keratin 10 (KRT10) downregulation. MDM2 binds to KRT10, targeting it for degradation via the ubiquitin-proteasome pathway. Preventing KRT10 degradation reduced VZV propagation in culture and prevented epidermal disruption in skin explants. KRT10 knockdown induced expression of NR4A1 and enhanced viral propagation in culture. NR4A1 knockdown prevented viral propagation in culture, reduced LC3 levels, and increased LAMP2 expression. We therefore describe a drug-able pathway whereby MDM2 ubiquitinates and degrades KRT10, increasing NR4A1 expression and allowing VZV replication and propagation., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

The dependence of viruses on the host cell to complete their replicative cycle renders cellular functions potential targets for novel antivirals. We screened a panel of broadly acting cellular ion channel inhibitors for activity against human cytomegalovirus (HCMV) and identified the voltage-gated chloride ion channel inhibitor 4,4'-diisothiocyano-2,2'-stilbenedisulfonic acid (DIDS) as a potent inhibitor of HCMV replication. Time-of-addition studies demonstrated that DIDS inhibited entry via direct interaction with the virion that impeded binding to the plasma membrane. Synthesis and analysis of pharmacological variants of DIDS suggested that intrinsic cysteine, and not lysine, reactivity was important for activity against HCMV. Although sequencing of DIDS-resistant HCMV revealed enrichment of a mutation within UL100 (encoding glycoprotein M) and a specific truncation of glycoprotein RL13, these did not explain the DIDS resistance phenotype. Specifically, only the introduction of the RL13 mutant partially phenocopied the DIDS resistance phenotype. Serendipitously, the entry of DIDS-resistant HCMV also became independent of heparan sulfate proteoglycans (HSPGs), suggesting that evasion of DIDS lowered dependence on an initial interaction with HSPGs. Intriguingly, the DIDS-resistant virus demonstrated increased sensitivity to antibody neutralization, which mapped, in part, to the presence of the gM mutation. Taken together the data characterize the antiviral activity of a novel HCMV inhibitor that drives HCMV infection to occur independently of HSPGs and the generation of increased sensitivity to humoral immunity. The data also demonstrate that compounds with cysteine reactivity have the potential to act as antiviral compounds against HCMV via direct engagement of virions. IMPORTANCE Human cytomegalovirus (HCMV) is major pathogen of nonimmunocompetent individuals that remains in need of new therapeutic options. Here, we identify a potent antiviral compound (4,4'-diisothiocyano-2,2'-stilbenedisulfonic acid [DIDS]), its mechanism of action, and the chemical properties required for its activity. In doing so, the data argue that cysteine-reactive compounds could have the capacity to be developed for anti-HCMV activity. Importantly, the data show that entry of DIDS-resistant virus became independent of heparan sulfate proteoglycans (HSPGs) but, concomitantly, became more sensitive to neutralizing antibody responses. This serendipitous observation suggests that retention of an interaction with HSPGs during the entry process in vivo may be evolutionarily advantageous through better evasion of humoral responses directed against HCMV virions., (Copyright © 2020 Murray et al.)

Herpes Simplex Virus type 1 (HSV-1) chronically infects over 70 per cent of the global population. Clinical manifestations are largely restricted to recurrent epidermal vesicles. However, HSV-1 also leads to encephalitis, the infection of the brain parenchyma, with high associated rates of mortality and morbidity. In this study, we performed target enrichment followed by direct sequencing of HSV-1 genomes, using target enrichment methods on the cerebrospinal fluid (CSF) of clinical encephalitis patients and from skin swabs of epidermal vesicles on non-encephalopathic patients. Phylogenetic analysis revealed high inter-host diversity and little population structure. In contrast, samples from different lesions in the same patient clustered with similar patterns of allelic variants. Comparison of consensus genome sequences shows HSV-1 has been freely recombining, except for distinct islands of linkage disequilibrium (LD). This suggests functional constraints prevent recombination between certain genes, notably those encoding pairs of interacting proteins. Distinct LD patterns characterised subsets of viruses recovered from CSF and skin lesions, which may reflect different evolutionary constraints in different body compartments. Functions of genes under differential constraint related to immunity or tropism and provide new hypotheses on tissue-specific mechanisms of viral infection and latency., (© The Author(s) 2020. Published by Oxford University Press.)

Herpes zoster (HZ) is caused by reactivation of varicella zoster virus (VZV) that established latency in sensory and autonomic neurons during primary infection. In the Shingles Prevention Study (SPS), a large efficacy trial of live attenuated Oka/Merck zoster vaccine (ZVL), PCR-confirmed second episodes of HZ occurred in two of 660 placebo and one of 321 ZVL recipients with documented HZ during a mean follow-up of 3.13 years. An additional two ZVL recipients experienced a second episode of HZ in the Long-Term Persistence Substudy. All episodes of HZ were caused by wild-type VZV. The first and second episodes of HZ occurred in different dermatomes in each of these five participants, with contralateral recurrences in two. Time between first and second episodes ranged from 12 to 28 months. One of the five participants, who was immunocompetent on study enrollment, was immunocompromised at the onset of his first and second episodes of HZ. VZV DNA isolated from rash lesions from the first and second episodes of HZ was used to sequence the full-length VZV genomes. For the unique-sequence regions of the VZV genome, we employed target enrichment of VZV DNA, followed by deep sequencing. For the reiteration regions, we used PCR amplification and Sanger sequencing. Our analysis and comparison of the VZV genomes from the first and second episodes of HZ in each of the five participants indicate that both episodes were caused by the same VZV strain. This is consistent with the extraordinary stability of VZV during the replication phase of varicella and the subsequent establishment of latency in sensory ganglia throughout the body. Our observations also indicate that VZV is stable during the persistence of latency and the subsequent reactivation and replication that results in HZ., (Copyright © 2019 The Authors. Published by Elsevier Ltd.. All rights reserved.)

The Glasgow s17 syn+ strain of herpes simplex virus 1 (HSV1) is arguably the best characterized strain and has provided the reference sequence for HSV1 genetic studies. Here we show that our original s17 syn+ stock was a mixed population from which we have isolated a minor variant that, unlike other strains in the laboratory, fails to be efficiently released from infected cells and spreads predominantly by direct cell-to-cell transmission. Analysis of other s17-derived viruses that had been isolated elsewhere revealed a number with the same release phenotype. Second-generation sequencing of 8 plaque-purified s17-derived viruses revealed sequences that vary by 50 single-nucleotide polymorphisms (SNPs), including approximately 10 coding SNPs. This compared to interstrain variations of around 800 SNPs in strain Sc16, of which a quarter were coding changes. Amongst the variations found within s17, we identified 13 variants of glycoprotein C within the original stock of virus that were predominantly a consequence of altered homopolymeric runs of C residues. Characterization of seven isolates coding for different forms of gC indicated that all were expressed, despite six of them lacking a transmembrane domain. While the release phenotype did not correlate directly with any of these identified gC variations, further demonstration that nine clinical isolates of HSV1 also fail to spread through extracellular release raises the possibility that propagation in tissue culture had altered the HSV1 s17 transmission phenotype. Hence, the s17 intrastrain variation identified here offers an excellent model for understanding both HSV1 transmission and tissue culture adaptation.

During latent infections with herpes simplex virus 1 (HSV-1), viral transcription is restricted and the genomes are mostly maintained in silenced chromatin, whereas in lytically infected cells all viral genes are transcribed and the genomes are dynamically chromatinized. Histones in the viral chromatin bear markers of silenced chromatin at early times in lytic infection or of active transcription at later times. The virion protein VP16 activates transcription of the immediate-early (IE) genes by recruiting transcription activators and chromatin remodelers to their promoters. Two IE proteins, ICP0 and ICP4 which modulate chromatin epigenetics, then activate transcription of early and late genes. Although chromatin is involved in the mechanism of activation of HSV- transcription, its precise role is not entirely understood. In the cellular genome, chromatin dynamics often modulate transcription competence whereas promoter-specific transcription factors determine transcription activity. Here, biophysical fractionation of serially digested HSV-1 chromatin followed by short-read deep sequencing indicates that nuclear HSV-1 DNA has different biophysical properties than protein-free or encapsidated HSV-1 DNA. The entire HSV-1 genomes in infected cells were equally accessible. The accessibility of transcribed or non-transcribed genes under any given condition did not differ, and each gene was entirely sampled in both the most and least accessible chromatin. However, HSV-1 genomes fractionated differently under conditions of generalized or restricted transcription. Approximately 1/3 of the HSV-1 DNA including fully sampled genes resolved to the most accessible chromatin when HSV-1 transcription was active, but such enrichment was reduced to only 3% under conditions of restricted HSV-1 transcription. Short sequences of restricted accessibility separated genes with different transcription levels. Chromatin dynamics thus provide a first level of regulation on HSV-1 transcription, dictating the transcriptional competency of the genomes during lytic infections, whereas the transcription of individual genes is then most likely activated by specific transcription factors. Moreover, genes transcribed to different levels are separated by short sequences with limited accessibility., Competing Interests: The authors have declared that no competing interests exist.

Epstein-Barr virus (EBV) is one of the most common viral infections in humans and persists within its host for life. EBV therefore represents an extremely successful virus that has evolved complex strategies to evade the host's innate and adaptive immune response during both initial and persistent stages of infection. Here, we conducted a comparative genomics analysis on 223 whole genome sequences of worldwide EBV strains. We recover extensive genome-wide linkage disequilibrium (LD) despite pervasive genetic recombination. This pattern is explained by the global EBV population being subdivided into three main subpopulations, one primarily found in East Asia, one in Southeast Asia and Oceania, and the third including most of the other globally distributed genomes we analyzed. Additionally, sites in LD were overrepresented in immunogenic genes. Taken together, our results suggest that host immune selection and local adaptation to different human host populations has shaped the genome-wide patterns of genetic diversity in EBV., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Recent sequencing efforts have led to estimates of human cytomegalovirus (HCMV) genome-wide intrahost diversity that rival those of persistent RNA viruses [Renzette N, Bhattacharjee B, Jensen JD, Gibson L, Kowalik TF (2011) PLoS Pathog 7:e1001344]. Here, we deep sequence HCMV genomes recovered from single and longitudinally collected blood samples from immunocompromised children to show that the observations of high within-host HCMV nucleotide diversity are explained by the frequent occurrence of mixed infections caused by genetically distant strains. To confirm this finding, we reconstructed within-host viral haplotypes from short-read sequence data. We verify that within-host HCMV nucleotide diversity in unmixed infections is no greater than that of other DNA viruses analyzed by the same sequencing and bioinformatic methods and considerably less than that of human immunodeficiency and hepatitis C viruses. By resolving individual viral haplotypes within patients, we reconstruct the timing, likely origins, and natural history of superinfecting strains. We uncover evidence for within-host recombination between genetically distinct HCMV strains, observing the loss of the parental virus containing the nonrecombinant fragment. The data suggest selection for strains containing the recombinant fragment, generating testable hypotheses about HCMV evolution and pathogenesis. These results highlight that high HCMV diversity present in some samples is caused by coinfection with multiple distinct strains and provide reassurance that within the host diversity for single-strain HCMV infections is no greater than for other herpesviruses., (Copyright © 2019 the Author(s). Published by PNAS.)

Characterizing complex viral transcriptomes by conventional RNA sequencing approaches is complicated by high gene density, overlapping reading frames, and complex splicing patterns. Direct RNA sequencing (direct RNA-seq) using nanopore arrays offers an exciting alternative whereby individual polyadenylated RNAs are sequenced directly, without the recoding and amplification biases inherent to other sequencing methodologies. Here we use direct RNA-seq to profile the herpes simplex virus type 1 (HSV-1) transcriptome during productive infection of primary cells. We show how direct RNA-seq data can be used to define transcription initiation and RNA cleavage sites associated with all polyadenylated viral RNAs and demonstrate that low level read-through transcription produces a novel class of chimeric HSV-1 transcripts, including a functional mRNA encoding a fusion of the viral E3 ubiquitin ligase ICP0 and viral membrane glycoprotein L. Thus, direct RNA-seq offers a powerful method to characterize the changing transcriptional landscape of viruses with complex genomes.

Transcriptome profiling has become routine in studies of many biological processes. However, the favored approaches such as short-read Illumina RNA sequencing are giving way to long-read sequencing platforms better suited to interrogating the complex transcriptomes typical of many RNA and DNA viruses. Here, we provide a guide-tailored to molecular virologists-to the ins and outs of viral transcriptome sequencing and discuss the strengths and weaknesses of the major RNA sequencing technologies as tools to analyze the abundance and diversity of the viral transcripts made during infection., (Copyright © 2018 American Society for Microbiology.)

Modification of mRNA by N 6 -adenosine methylation (m 6 A) on internal bases influences gene expression in eukaryotes. How the dynamic genome-wide landscape of m 6 A-modified mRNAs impacts virus infection and host immune responses remains poorly understood. Here, we show that type I interferon (IFN) production triggered by dsDNA or human cytomegalovirus (HCMV) is controlled by the cellular m 6 A methyltrasferase subunit METTL14 and ALKBH5 demethylase. While METTL14 depletion reduced virus reproduction and stimulated dsDNA- or HCMV-induced IFNB1 mRNA accumulation, ALKBH5 depletion had the opposite effect. Depleting METTL14 increased both nascent IFNB1 mRNA production and stability in response to dsDNA. In contrast, ALKBH5 depletion reduced nascent IFNB1 mRNA production without detectably influencing IFN1B mRNA decay. Genome-wide transcriptome profiling following ALKBH5 depletion identified differentially expressed genes regulating antiviral immune responses, while METTL14 depletion altered pathways impacting metabolic reprogramming, stress responses, and aging. Finally, we determined that IFNB1 mRNA was m 6 A-modified within both the coding sequence and the 3' untranslated region (UTR). This establishes that the host m 6 A modification machinery controls IFNβ production triggered by HCMV or dsDNA. Moreover, it demonstrates that responses to nonmicrobial dsDNA in uninfected cells, which shape host immunity and contribute to autoimmune disease, are regulated by enzymes controlling m 6 A epitranscriptomic changes., (© 2018 Rubio et al.; Published by Cold Spring Harbor Laboratory Press.)

Background: Varicella zoster virus (VZV) may cause encephalitis, both with and without rash. Here we investigate whether viruses recovered from the central nervous system (CNS; encephalitis or meningitis) differ genetically from those recovered from non-CNS samples., Methods: Enrichment-based deep sequencing of 45 VZV genomes from cerebral spinal fluid (CSF), plasma, bronchoalveolar lavage (BAL), and vesicles was carried out with samples collected from 34 patients with and without VZV infection of the CNS., Results: Viral sequences from multiple sites in the same patient were identical at the consensus level. Virus from vesicle fluid and CSF in cases of meningitis showed low-level diversity. By contrast, plasma, BAL, and encephalitis had higher numbers of variant alleles. Two CSF-encephalitis samples had high genetic diversity, with variant frequency patterns typical of mixed infections with different clades., Conclusions: Low viral genetic diversity in vesicle fluid is compatible with previous observations that VZV skin lesions arise from single or low numbers of virions. A similar result was observed in VZV from cases of VZV meningitis, a generally self-limiting infection. CSF from cases of encephalitis had higher diversity with evidence for mixed clade infections in 2 cases. We hypothesize that reactivation from multiple neurons may contribute to the pathogenesis of VZV encephalitis.

Background: Adenoviruses are significant pathogens for the immunocompromised, arising from primary infection or reinfection. Serotyping is insufficient to support nosocomial transmission investigations. We investigate whether whole-genome sequencing (WGS) provides clinically relevant information on transmission among patients in a pediatric tertiary hospital., Methods: We developed a target-enriched adenovirus WGS technique for clinical samples and retrospectively sequenced 107 adenovirus-positive residual diagnostic samples, including viremias (>5 × 104 copies/mL), from 37 patients collected January 2011-March 2016. Whole-genome sequencing was used to determine genotype and for phylogenetic analysis., Results: Adenovirus sequences were recovered from 105 of 107 samples. Full genome sequences were recovered from all 20 nonspecies C samples and from 36 of 85 species C viruses, with partial genome sequences recovered from the rest. Whole-genome phylogenetic analysis suggested linkage of 3 genotype A31 cases and uncovered an unsuspected epidemiological link to an A31 infection first detected on the same ward 4 years earlier. In 9 samples from 1 patient who died, we identified a mixed genotype adenovirus infection., Conclusions: Adenovirus WGS from clinical samples is possible and useful for genotyping and molecular epidemiology. Whole-genome sequencing identified likely nosocomial transmission with greater resolution than conventional genotyping and distinguished between adenovirus disease due to single or multiple genotypes.

Primary varicella-zoster virus (VZV) infection causes varicella (chickenpox) and the establishment of a lifelong latent infection in ganglionic neurons. VZV reactivates in about one-third of infected individuals to cause herpes zoster, often accompanied by neurological complications. The restricted host range of VZV and, until recently, a lack of suitable in vitro models have seriously hampered molecular studies of VZV latency. Nevertheless, recent technological advances facilitated a series of exciting studies that resulted in the discovery of a VZV latency-associated transcript (VLT) and provide novel insights into our understanding of VZV latency and factors that may initiate reactivation. Deducing the function(s) of VLT and the molecular mechanisms involved should now be considered a priority to improve our understanding of factors that govern VZV latency and reactivation. In this review, we summarize the implications of recent discoveries in the VZV latency field from both a virus and host perspective and provide a roadmap for future studies.

Varicella-zoster virus (VZV), an alphaherpesvirus, establishes lifelong latent infection in the neurons of >90% humans worldwide, reactivating in one-third to cause shingles, debilitating pain and stroke. How VZV maintains latency remains unclear. Here, using ultra-deep virus-enriched RNA sequencing of latently infected human trigeminal ganglia (TG), we demonstrate the consistent expression of a spliced VZV mRNA, antisense to VZV open reading frame 61 (ORF61). The spliced VZV latency-associated transcript (VLT) is expressed in human TG neurons and encodes a protein with late kinetics in productively infected cells in vitro and in shingles skin lesions. Whereas multiple alternatively spliced VLT isoforms (VLT ly ) are expressed during lytic infection, a single unique VLT isoform, which specifically suppresses ORF61 gene expression in co-transfected cells, predominates in latently VZV-infected human TG. The discovery of VLT links VZV with the other better characterized human and animal neurotropic alphaherpesviruses and provides insights into VZV latency.

Background: Human herpesvirus-6A and -6B (HHV-6) are betaherpesviruses that reach > 90% seroprevalence in the adult population. Unique among human herpesviruses, HHV-6 can integrate into the subtelomeric regions of human chromosomes; when this occurs in germ line cells it causes a condition called inherited chromosomally integrated HHV-6 (iciHHV-6). Only two complete genomes are available for replicating HHV-6B, leading to numerous conflicting annotations and little known about the global genomic diversity of this ubiquitous virus., Results: Using a custom capture panel for HHV-6B, we report complete genomes from 61 isolates of HHV-6B from active infections (20 from Japan, 35 from New York state, and 6 from Uganda), and 64 strains of iciHHV-6B (mostly from North America). HHV-6B sequence clustered by geography and illustrated extensive recombination. Multiple iciHHV-6B sequences from unrelated individuals across the United States were found to be completely identical, consistent with a founder effect. Several iciHHV-6B strains clustered with strains from recent active pediatric infection. Combining our genomic analysis with the first RNA-Seq and shotgun proteomics studies of HHV-6B, we completely reannotated the HHV-6B genome, altering annotations for more than 10% of existing genes, with multiple instances of novel splicing and genes that hitherto had gone unannotated., Conclusion: Our results are consistent with a model of intermittent de novo integration of HHV-6B into host germline cells during active infection with a large contribution of founder effect in iciHHV-6B. Our data provide a significant advance in the genomic annotation of HHV-6B, which will contribute to the detection, diversity, and control of this virus.