Histone H2A variant H2A.B is enriched in transcriptionally active and replicating HSV-1 lytic chromatin

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.

Reanalysis of single-cell RNA sequencing data does not support herpes simplex virus 1 latency in non-neuronal ganglionic cells in mice

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.

Nanopore Guided Annotation of Transcriptome Architectures

High-resolution annotations of transcriptomes from all domains of life are essential for many sequencing-based RNA analyses, including Nanopore direct RNA sequencing (DRS), which would otherwise be hindered by misalignments and other analysis artefacts. DRS allows the capture and full-length sequencing of native RNAs, without recoding or amplification bias, and resulting data 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. Existing software solutions for generating high-resolution transcriptome annotations are poorly suited to small gene dense organisms such as 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 datasets and developed a novel approach to transcriptome. We demonstrate, using a combination of synthetic and original datasets, that our novel approach yields a high level of precision and recall when reconstructing both gene sparse and gene dense transcriptomes from DRS datasets. We further apply this approach to generate a new high resolution transcriptome annotation of the neglected pathogen human adenovirus type F 41 for which we identify 77 distinct transcripts encoding at least 23 different proteins.

tRNA epitranscriptome determines pathogenicity of the opportunistic pathogen Pseudomonas aeruginosa

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 (cmnm5U34) 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.

Histone H2A variant H2A.B is enriched in transcriptionally active HSV-1 lytic chromatin

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 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, ectopic, and endogenous H2B, were assembled into HSV-1 chromatin evenly throughout the genome, but canonical H2A was relatively depleted from the viral chromatin whereas H2A.B was enriched in the most dynamic viral chromatin. When viral transcription was 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

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 mechanisms of regulation of HSV-1 transcription by epigenetics have not been fully characterized and may differ from those regulating cellular transcription. In particular, the lytic HSV-1 chromatin is unusually dynamic, whereas the latent silenced one is not, but the mechanisms resulting in the unique dynamics of the lytic chromatin remain unknown. Here we identify the enrichment on the highly dynamic histone 2A variant H2A in the most dynamic viral chromatin, which provides a mechanistic understanding for 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.

CCR4-NOT differentially controls host versus virus poly(a)-tail length and regulates HCMV infection

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.

Reanalysis of single-cell RNA sequencing data does not support herpes simplex virus 1 latency in non-neuronal ganglionic cells in mice

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 virus is also present in immune cells recovered from ganglia in a mouse model used for studying latency. Here, we reanalyzed the single-cell RNA sequencing (scRNA-Seq) data that formed the basis for this 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 nearexclusively detected in a mixed population of cells undergoing cell death. Specific loss of HSV1 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 by inaccuracies in the data analyses.

Novel viral splicing events and open reading frames revealed by long-read direct RNA sequencing of adenovirus transcripts

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.

DLK-Dependent Biphasic Reactivation of Herpes Simplex Virus Latency Established in the Absence of Antivirals

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.

DRUMMER—Rapid detection of RNA modifications through comparative nanopore sequencing

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.

Single-cell transcriptomics identifies Gadd45b as a regulator of herpesvirus-reactivating neurons

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.

Epstein-Barr Virus-Encoded BILF1 Orthologues From Porcine Lymphotropic Herpesviruses Display Common Molecular Functionality

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.

The architecture of the simian varicella virus transcriptome

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.

Varicella-zoster virus (VZV)–a ubiquitous human pathogen–infects most individuals during childhood, leading to chickenpox, after which the virus persists in the host for decades. Later in life, VZV reactivates to cause shingles, frequently associated with difficult-to-treat chronic pain. Our limited understanding of the viral life-cycle hampers the development of more effective treatment options. Simian varicella virus (SVV) is the non-human primate homologue of VZV and causes a natural disease in Old World monkeys with clinical, pathological, and immunological features resembling human VZV infection. However, it is unclear how similar both viruses are at the molecular level. Here, we have revisited the genome-wide transcriptional activity of SVV during lytic infection of kidney epithelial cells derived from two non-human primate species and validated expression of newly identified viral transcripts in lung tissue from SVV-infected animals. Together, this has led to the identification of numerous alternative RNA isoforms, mostly unique to SVV, and some of which may have functional implications for the virus. Notably, we defined the SVV latency-associated transcript, which is highly similar to its VZV counterpart. In conclusion, our study shows the value of understanding the molecular biology of a given animal model and identifies potentially conserved mechanism of latency.

Mutagenesis of the Varicella-Zoster Virus Genome Demonstrates That VLT and VLT-ORF63 Proteins Are Dispensable for Lytic Infection

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.

ICTV Virus Taxonomy Profile: Herpesviridae 2021

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.

Varicella-Zoster Virus-Genetics, Molecular Evolution and Recombination

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.

Widespread remodeling of the m6A RNA-modification landscape by a viral regulator of RNA processing and export

N6-methyladenosine (m6A) 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 m6A pathway is not known. Here, we show that, in primary fibroblasts, HSV-1 orchestrates a striking redistribution of the nuclear m6A 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 m6A-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 m6A methyltransferase but becomes less impacted as the infection advances. Redistribution of the nuclear m6A machinery is accompanied by a wide-scale reduction in the installation of m6A 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.

Targeting the m6A RNA modification pathway blocks SARS-CoV-2 and HCoV-OC43 replication

N6-methyladenosine (m6A) is an abundant internal RNA modification, influencing transcript fate and function in uninfected and virus-infected cells. Installation of m6A by the nuclear RNA methyltransferase METTL3 occurs cotranscriptionally; however, the genomes of some cytoplasmic RNA viruses are also m6A-modified. How the cellular m6A 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 m6A 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 m6A modification on genomic and subgenomic RNAs of both viruses were mapped by methylated RNA immunoprecipitation sequencing (meRIP-seq). Levels of host factors involved in m6A installation, removal, and recognition were unchanged by HCoV-OC43 infection; however, nuclear localization of METTL3 and cytoplasmic m6A readers YTHDF1 and YTHDF2 increased. This establishes that coronavirus RNAs are m6A-modified and host m6A pathway components control β-coronavirus replication. Moreover, it illustrates the therapeutic potential of targeting the m6A pathway to restrict coronavirus reproduction.

Enteric viruses evoke broad host immune responses resembling those elicited by the bacterial microbiome

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.

Author Correction: The influence of human genetic variation on Epstein-Barr virus sequence diversity

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The influence of human genetic variation on Epstein-Barr virus sequence diversity

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/mm3 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 VLT-ORF63 fusion transcript induces broad viral gene expression during reactivation from neuronal latency

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.

Direct RNA sequencing reveals m6A modifications on adenovirus RNA are necessary for efficient splicing

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 (m6A) to mRNAs. While m6A-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 m6A mapping using short-read sequencing, here we profile m6A within the adenovirus transcriptome using a combination of meRIP-seq and direct RNA long-read sequencing to yield both nucleotide and transcript-resolved m6A detection. Although both early and late viral transcripts contain m6A, depletion of m6A writer METTL3 specifically impacts viral late transcripts by reducing their splicing efficiency. These data showcase a new technique for m6A discovery within individual transcripts at nucleotide resolution, and highlight the role of m6A in regulating splicing of a viral pathogen.

Decoding the Architecture of the Varicella-Zoster Virus Transcriptome

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.

Using Direct RNA Nanopore Sequencing to Deconvolute Viral Transcriptomes

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 N6 -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.

Analysis of the reiteration regions (R1 to R5) of varicella-zoster virus

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.

Recurrent herpes zoster in the Shingles Prevention Study: Are second episodes caused by the same varicella-zoster virus strain?

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.

Whole genome sequencing of Herpes Simplex Virus 1 directly from human cerebrospinal fluid reveals selective constraints in neurotropic viruses

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.

Kallikrein-Mediated Cytokeratin 10 Degradation Is Required for Varicella Zoster Virus Propagation in Skin

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.

High Viral Diversity and Mixed Infections in Cerebral Spinal Fluid From Cases of Varicella Zoster Virus Encephalitis

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.

Co-evolution of sites under immune selection shapes Epstein-Barr Virus population structure

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.

Direct RNA sequencing on nanopore arrays redefines the transcriptional complexity of a viral pathogen

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.

Going the Distance: Optimizing RNA-Seq Strategies for Transcriptomic Analysis of Complex Viral 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.

RNA m6 A modification enzymes shape innate responses to DNA by regulating interferon β

In this study, Rubio et al. researched how the dynamic genome-wide landscape of m⁶A-modified mRNAs impacts virus infection and host immune responses. They show that type I interferon (IFN) production triggered by dsDNA or human cytomegalovirus (HCMV) is controlled by the cellular m⁶A methyltrasferase subunit METTL14 and ALKBH5 demethylase, and their results demonstrate that responses to nonmicrobial dsDNA in uninfected cells are regulated by enzymes controlling m⁶A epitranscriptomic changes.

Modification of mRNA by N⁶-adenosine methylation (m⁶A) on internal bases influences gene expression in eukaryotes. How the dynamic genome-wide landscape of m⁶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⁶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⁶A-modified within both the coding sequence and the 3′ untranslated region (UTR). This establishes that the host m⁶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⁶A epitranscriptomic changes.

2018 Colorado Alphaherpesvirus Latency Society Symposium

Meeting Report on the 8th Annual Symposium of the Colorado Alphaherpesvirus Latency Society (CALS), held on May 16-19, 2018, in Vail, Colorado.

Use of Whole-Genome Sequencing of Adenovirus in Immunocompromised Pediatric Patients to Identify Nosocomial Transmission and Mixed-Genotype Infection

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.

Acute Retinal Necrosis Caused by the Zoster Vaccine Virus

We report acute retinal necrosis caused by the vaccine Oka strain following immunization of a 78-year-old woman with live zoster vaccine. Whole genome sequencing confirmed the ocular vOka strain to be derived from the vaccine and excluded the presence of new mutations or recombination with wild-type Varicella zoster virus.

Molecular Aspects of Varicella-Zoster Virus Latency

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.

Comparative genomic, transcriptomic, and proteomic reannotation of human herpesvirus 6

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4604-2) contains supplementary material, which is available to authorized users.

Viral Genome Sequencing Proves Nosocomial Transmission of Fatal Varicella

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.

Detection of Low Frequency Multi-Drug Resistance and Novel Putative Maribavir Resistance in Immunocompromised Pediatric Patients with Cytomegalovirus

Human cytomegalovirus (HCMV) is a significant pathogen in immunocompromised individuals, with the potential to cause fatal pneumonitis and colitis, as well as increasing the risk of organ rejection in transplant patients. With the advent of new anti-HCMV drugs there is therefore considerable interest in using virus sequence data to monitor emerging resistance to antiviral drugs in HCMV viraemia and disease, including the identification of putative new mutations. We used target-enrichment to deep sequence HCMV DNA from 11 immunosuppressed pediatric patients receiving single or combination anti-HCMV treatment, serially sampled over 1–27 weeks. Changes in consensus sequence and resistance mutations were analyzed for three ORFs targeted by anti-HCMV drugs and the frequencies of drug resistance mutations monitored. Targeted-enriched sequencing of clinical material detected mutations occurring at frequencies of 2%. Seven patients showed no evidence of drug resistance mutations. Four patients developed drug resistance mutations a mean of 16 weeks after starting treatment. In two patients, multiple resistance mutations accumulated at frequencies of 20% or less, including putative maribavir and ganciclovir resistance mutations P522Q (UL54) and C480F (UL97). In one patient, resistance was detected 14 days earlier than by PCR. Phylogenetic analysis suggested recombination or superinfection in one patient. Deep sequencing of HCMV enriched from clinical samples excluded resistance in 7 of 11 subjects and identified resistance mutations earlier than conventional PCR-based resistance testing in 2 patients. Detection of multiple low level resistance mutations was associated with poor outcome.

In vitro system using human neurons demonstrates that varicella-zoster vaccine virus is impaired for reactivation, but not latency

Varicella-zoster virus (VZV) establishes latency in human sensory and cranial nerve ganglia during primary infection (varicella), and the virus can reactivate and cause zoster after primary infection. The mechanism of how the virus establishes and maintains latency and how it reactivates is poorly understood, largely due to the lack of robust models. We found that axonal infection of neurons derived from hESCs in a microfluidic device with cell-free parental Oka (POka) VZV resulted in latent infection with inability to detect several viral mRNAs by reverse transcriptase-quantitative PCR, no production of infectious virus, and maintenance of the viral DNA genome in endless configuration, consistent with an episome configuration. With deep sequencing, however, multiple viral mRNAs were detected. Treatment of the latently infected neurons with Ab to NGF resulted in production of infectious virus in about 25% of the latently infected cultures. Axonal infection of neurons with vaccine Oka (VOka) VZV resulted in a latent infection similar to infection with POka; however, in contrast to POka, VOka-infected neurons were markedly impaired for reactivation after treatment with Ab to NGF. In addition, viral transcription was markedly reduced in neurons latently infected with VOka compared with POka. Our in vitro system recapitulates both VZV latency and reactivation in vivo and may be used to study viral vaccines for their ability to establish latency and reactivate.

Islands of linkage in an ocean of pervasive recombination reveals two-speed evolution of human cytomegalovirus genomes

Human cytomegalovirus (HCMV) infects most of the population worldwide, persisting throughout the host's life in a latent state with periodic episodes of reactivation. While typically asymptomatic, HCMV can cause fatal disease among congenitally infected infants and immunocompromised patients. These clinical issues are compounded by the emergence of antiviral resistance and the absence of an effective vaccine, the development of which is likely complicated by the numerous immune evasins encoded by HCMV to counter the host's adaptive immune responses, a feature that facilitates frequent super-infections. Understanding the evolutionary dynamics of HCMV is essential for the development of effective new drugs and vaccines. By comparing viral genomes from uncultivated or low-passaged clinical samples of diverse origins, we observe evidence of frequent homologous recombination events, both recent and ancient, and no structure of HCMV genetic diversity at the whole-genome scale. Analysis of individual gene-scale loci reveals a striking dichotomy: while most of the genome is highly conserved, recombines essentially freely and has evolved under purifying selection, 21 genes display extreme diversity, structured into distinct genotypes that do not recombine with each other. Most of these hyper-variable genes encode glycoproteins involved in cell entry or escape of host immunity. Evidence that half of them have diverged through episodes of intense positive selection suggests that rapid evolution of hyper-variable loci is likely driven by interactions with host immunity. It appears that this process is enabled by recombination unlinking hyper-variable loci from strongly constrained neighboring sites. It is conceivable that viral mechanisms facilitating super-infection have evolved to promote recombination between diverged genotypes, allowing the virus to continuously diversify at key loci to escape immune detection, while maintaining a genome optimally adapted to its asymptomatic infectious lifecycle.

Rates of vaccine evolution show strong effects of latency: implications for varicella zoster virus epidemiology

Varicella-zoster virus (VZV) causes chickenpox and shingles, and is found in human populations worldwide. The lack of temporal signal in the diversity of VZV makes substitution rate estimates unreliable, which is a barrier to understanding the context of its global spread. Here, we estimate rates of evolution by studying live attenuated vaccines, which evolved in 22 vaccinated patients for known periods of time, sometimes, but not always undergoing latency. We show that the attenuated virus evolves rapidly (∼ 10(-6) substitutions/site/day), but that rates decrease dramatically when the virus undergoes latency. These data are best explained by a model in which viral populations evolve for around 13 days before becoming latent, but then undergo no replication during latency. This implies that rates of viral evolution will depend strongly on transmission patterns. Nevertheless, we show that implausibly long latency periods are required to date the most recent common ancestor of extant VZV to an "out-of-Africa" migration with humans, as has been previously suggested.

Whole-genome enrichment and sequencing of Chlamydia trachomatis directly from clinical samples

Background

Chlamydia trachomatis is a pathogen of worldwide importance, causing more than 100 million cases of sexually transmitted infections annually. Whole-genome sequencing is a powerful high resolution tool that can be used to generate accurate data on bacterial population structure, phylogeography and mutations associated with antimicrobial resistance. The objective of this study was to perform whole-genome enrichment and sequencing of C. trachomatis directly from clinical samples.

Methods

C. trachomatis positive samples comprising seven vaginal swabs and three urine samples were sequenced without prior in vitro culture in addition to nine cultured C. trachomatis samples, representing different serovars. A custom capture RNA bait set, that captures all known diversity amongst C. trachomatis genomes, was used in a whole-genome enrichment step during library preparation to enrich for C. trachomatis DNA. All samples were sequenced on the MiSeq platform.

Results

Full length C. trachomatis genomes (>95-100% coverage of a reference genome) were successfully generated for eight of ten clinical samples and for all cultured samples. The proportion of reads mapping to C. trachomatis and the mean read depth across each genome were strongly linked to the number of bacterial copies within the original sample. Phylogenetic analysis confirmed the known population structure and the data showed potential for identification of minority variants and mutations associated with antimicrobial resistance. The sensitivity of the method was >10-fold higher than other reported methodologies.

Conclusions

The combination of whole-genome enrichment and deep sequencing has proven to be a non-mutagenic approach, capturing all known variation found within C. trachomatis genomes. The method is a consistent and sensitive tool that enables rapid whole-genome sequencing of C. trachomatis directly from clinical samples and has the potential to be adapted to other pathogens with a similar clonal nature.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-014-0591-3) contains supplementary material, which is available to authorized users.

Deep sequencing of viral genomes provides insight into the evolution and pathogenesis of varicella zoster virus and its vaccine in humans

Immunization with the vOka vaccine prevents varicella (chickenpox) in children and susceptible adults. The vOka vaccine strain comprises a mixture of genotypes and, despite attenuation, causes rashes in small numbers of recipients. Like wild-type virus, the vaccine establishes latency in neuronal tissue and can later reactivate to cause Herpes zoster (shingles). Using hybridization-based methodologies, we have purified and sequenced vOka directly from skin lesions. We show that alleles present in the vaccine can be recovered from the lesions and demonstrate the presence of a severe bottleneck between inoculation and lesion formation. Genotypes in any one lesion appear to be descended from one to three vaccine-genotypes with a low frequency of novel mutations. No single vOka haplotype and no novel mutations are consistently present in rashes, indicating that neither new mutations nor recombination with wild type are critical to the evolution of vOka rashes. Instead, alleles arising from attenuation (i.e., not derived from free-living virus) are present at lower frequencies in rash genotypes. We identify 11 loci at which the ancestral allele is selected for in vOka rash formation and show genotypes in rashes that have reactivated from latency cannot be distinguished from rashes occurring immediately after inoculation. We conclude that the vOka vaccine, although heterogeneous, has not evolved to form rashes through positive selection in the mode of a quasispecies, but rather alleles that were essentially neutral during the vaccine production have been selected against in the human subjects, allowing us to identify key loci for rash formation.

Next-generation whole genome sequencing identifies the direction of norovirus transmission in linked patients

BACKGROUND

Noroviruses are a highly transmissible and major cause of nosocomial gastroenteritis resulting in bed and hospital-ward closures. Where hospital outbreaks are suspected, it is important to determine the routes of spread so that appropriate infection-control procedures can be implemented. To investigate a cluster of norovirus cases occurring in children undergoing bone marrow transplant, we undertook norovirus genome sequencing by next-generation methods. Detailed comparison of sequence data from 2 linked cases enabled us to identify the likely direction of spread.

METHODS

Norovirus complementary DNA was amplified by overlapping polymerase chain reaction (PCR) from 13 stool samples from 5 diagnostic real-time PCR-positive patients. The amplicons were sequenced by Roche 454, the genomes assembled by de novo assembly, and the data analyzed phylogenetically.

RESULTS

Phylogenetic analysis indicated that patients were infected by viruses similar to 4 distinct GII.4 subtypes and 2 patients were linked by the same virus. Of the 14 sites at which there were differences between the consensus sequences of the 2 linked viral genomes, 9 had minor variants present within one or the other patient. Further analysis confirmed that minor variants at all 9 sites in patient B w ere present as the consensus sequence in patient A.

CONCLUSIONS

Phylogenetic analysis excluded a common source of infection in this apparent outbreak. Two of 3 patients on the same ward had closely related viruses, raising the possibility of cross-infection despite protective isolation. Analysis of deep sequencing data enabled us to establish the likely direction of nosocomial transmission.

Viral population analysis and minority-variant detection using short read next-generation sequencing

RNA viruses within infected individuals exist as a population of evolutionary-related variants. Owing to evolutionary change affecting the constitution of this population, the frequency and/or occurrence of individual viral variants can show marked or subtle fluctuations. Since the development of massively parallel sequencing platforms, such viral populations can now be investigated to unprecedented resolution. A critical problem with such analyses is the presence of sequencing-related errors that obscure the identification of true biological variants present at low frequency. Here, we report the development and assessment of the Quality Assessment of Short Read (QUASR) Pipeline (http://sourceforge.net/projects/quasr) specific for virus genome short read analysis that minimizes sequencing errors from multiple deep-sequencing platforms, and enables post-mapping analysis of the minority variants within the viral population. QUASR significantly reduces the error-related noise in deep-sequencing datasets, resulting in increased mapping accuracy and reduction of erroneous mutations. Using QUASR, we have determined influenza virus genome dynamics in sequential samples from an in vitro evolution of 2009 pandemic H1N1 (A/H1N1/09) influenza from samples sequenced on both the Roche 454 GSFLX and Illumina GAIIx platforms. Importantly, concordance between the 454 and Illumina sequencing allowed unambiguous minority-variant detection and accurate determination of virus population turnover in vitro.

Specific capture and whole-genome sequencing of viruses from clinical samples

Whole genome sequencing of viruses directly from clinical samples is integral for understanding the genetics of host-virus interactions. Here, we report the use of sample sparing target enrichment (by hybridisation) for viral nucleic acid separation and deep-sequencing of herpesvirus genomes directly from a range of clinical samples including saliva, blood, virus vesicles, cerebrospinal fluid, and tumour cell lines. We demonstrate the effectiveness of the method by deep-sequencing 13 highly cell-associated human herpesvirus genomes and generating full length genome alignments at high read depth. Moreover, we show the specificity of the method enables the study of viral population structures and their diversity within a range of clinical samples types.

Chromosome and gene copy number variation allow major structural change between species and strains of Leishmania

Leishmania parasites cause a spectrum of clinical pathology in humans ranging from disfiguring cutaneous lesions to fatal visceral leishmaniasis. We have generated a reference genome for Leishmania mexicana and refined the reference genomes for Leishmania major, Leishmania infantum, and Leishmania braziliensis. This has allowed the identification of a remarkably low number of genes or paralog groups (2, 14, 19, and 67, respectively) unique to one species. These were found to be conserved in additional isolates of the same species. We have predicted allelic variation and find that in these isolates, L. major and L. infantum have a surprisingly low number of predicted heterozygous SNPs compared with L. braziliensis and L. mexicana. We used short read coverage to infer ploidy and gene copy numbers, identifying large copy number variations between species, with 200 tandem gene arrays in L. major and 132 in L. mexicana. Chromosome copy number also varied significantly between species, with nine supernumerary chromosomes in L. infantum, four in L. mexicana, two in L. braziliensis, and one in L. major. A significant bias against gene arrays on supernumerary chromosomes was shown to exist, indicating that duplication events occur more frequently on disomic chromosomes. Taken together, our data demonstrate that there is little variation in unique gene content across Leishmania species, but large-scale genetic heterogeneity can result through gene amplification on disomic chromosomes and variation in chromosome number. Increased gene copy number due to chromosome amplification may contribute to alterations in gene expression in response to environmental conditions in the host, providing a genetic basis for disease tropism.