Genome-wide Transcript Structure Resolution Reveals Abundant Alternate Isoform Usage from Murine Gammaherpesvirus 68

High-density resolution of the Kaposi's sarcoma associated herpesvirus transcriptome identifies novel transcript isoforms generated by long-range transcription and alternative splicing

Kaposi's sarcoma-associated herpesvirus is the etiologic agent of Kaposi's sarcoma and two B-cell malignancies. Recent advancements in sequencing technologies have led to high resolution transcriptomes for several human herpesviruses that densely encode genes on both strands. However, for KSHV progress remained limited due to the overall low percentage of KSHV transcripts, even during lytic replication. To address this challenge, we have developed a target enrichment method to increase the KSHV-specific reads for both short- and long-read sequencing platforms. Furthermore, we combined this approach with the Transcriptome Resolution through Integration of Multi-platform Data (TRIMD) pipeline developed previously to annotate transcript structures. TRIMD first builds a scaffold based on long-read sequencing and validates each transcript feature with supporting evidence from Illumina RNA-Seq and deepCAGE sequencing data. Our stringent innovative approach identified 994 unique KSHV transcripts, thus providing the first high-density KSHV lytic transcriptome. We describe a plethora of novel coding and non-coding KSHV transcript isoforms with alternative untranslated regions, splice junctions and open-reading frames, thus providing deeper insights on gene expression regulation of KSHV. Interestingly, as described for Epstein-Barr virus, we identified transcription start sites that augment long-range transcription and may increase the number of latency-associated genes potentially expressed in KS tumors.

Alternative splicing of CARM1 regulated by LincGET-guided paraspeckles biases the first cell fate in mammalian early embryos

The heterogeneity of CARM1 controls first cell fate bias during early mouse development. However, how this heterogeneity is established is unknown. Here, we show that Carm1 mRNA is of a variety of specific exon-skipping splicing (ESS) isoforms in mouse two-cell to four-cell embryos that contribute to CARM1 heterogeneity. Disruption of paraspeckles promotes the ESS of Carm1 precursor mRNAs (pre-mRNAs). LincGET, but not Neat1, is required for paraspeckle assembly and inhibits the ESS of Carm1 pre-mRNAs in mouse two-cell to four-cell embryos. We further find that LincGET recruits paraspeckles to the Carm1 gene locus through HNRNPU. Interestingly, PCBP1 binds the Carm1 pre-mRNAs and promotes its ESS in the absence of LincGET. Finally, we find that the ESS seen in mouse two-cell to four-cell embryos decreases CARM1 protein levels and leads to trophectoderm fate bias. Our findings demonstrate that alternative splicing of CARM1 has an important role in first cell fate determination.

Open reading frames M12/M13 jointly contribute to MHV-68 latency

Murine gammaherpesvirus 68 (MHV-68), a widely used small-animal model for the analysis of gammaherpesvirus pathogenesis, encodes the MHV-68-specific ORFs M12 and M13. The function of M12 and M13 has not been investigated so far. Therefore, we constructed and analysed recombinant MHV-68 with mutations in either M12, M13 or M12/M13. Both the M12 and M13 mutants did not display any phenotype in vitro or in vivo. However, although the M12/13 double mutant showed similar lytic growth in fibroblasts in vitro and in the lungs of infected mice as wild-type MHV-68, it was significantly attenuated in vivo during latency. This phenotype was completely restored in a revertant of the M12/13 double mutant. Thus, it appears that M12 and M13 might have redundant functions that are only revealed if both genes are lacking. The observation that M12/13 have a function during latency not only contributes to the further understanding of the pathogenesis of MHV-68 infection but might also be of interest considering that M12/13 are located at a genomic position similar to that of LMP2A and K15. The latter are important proteins of their respective human gammaherpesviruses EBV and KSHV that contribute to cellular survival, cell activation and proliferation, which was deduced from in vitro studies.

Lytic Replication and Reactivation from B Cells Is Not Required for Establishing or Maintaining Gammaherpesvirus Latency In Vivo

Gammaherpesviruses (GHVs) are lymphotropic tumor viruses with a biphasic infectious cycle. Lytic replication at the primary site of infection is necessary for GHVs to spread throughout the host and establish latency in distal sites. Dissemination is mediated by infected B cells that traffic hematogenously from draining lymph nodes to peripheral lymphoid organs, such as the spleen. B cells serve as the major reservoir for viral latency, and it is hypothesized that periodic reactivation from latently infected B cells contributes to maintaining long-term chronic infection. While fundamentally important to an understanding of GHV biology, aspects of B cell infection in latency establishment and maintenance are incompletely defined, especially roles for lytic replication and reactivation in this cell type. To address this knowledge gap and overcome limitations of replication-defective viruses, we generated a recombinant murine gammaherpesvirus 68 (MHV68) in which ORF50, the gene that encodes the essential immediate-early replication and transcription activator protein (RTA), was flanked by loxP sites to enable conditional ablation of lytic replication by ORF50 deletion in cells that express Cre recombinase. Following infection of mice that encode Cre in B cells with this virus, splenomegaly and viral reactivation from splenocytes were significantly reduced; however, the number of latently infected splenocytes was equivalent to WT MHV68. Despite ORF50 deletion, MHV68 latency was maintained over time in spleens of mice at levels approximating WT, reactivation-competent MHV68. Treatment of infected mice with lipopolysaccharide (LPS), which promotes B cell activation and MHV68 reactivation ex vivo, yielded equivalent increases in the number of latently infected cells for both ORF50-deleted and WT MHV68, even when mice were simultaneously treated with the antiviral drug cidofovir to prevent reactivation. Together, these data demonstrate that productive viral replication in B cells is not required for MHV68 latency establishment and support the hypothesis that B cell proliferation facilitates latency maintenance in vivo in the absence of reactivation. IMPORTANCE Gammaherpesviruses establish lifelong chronic infections in cells of the immune system and place infected hosts at risk for developing lymphomas and other diseases. It is hypothesized that gammaherpesviruses must initiate acute infection in these cells to establish and maintain long-term infection, but this has not been directly tested. We report here the use of a viral genetic system that allows for cell-type-specific deletion of a viral gene that is essential for replication and reactivation. We employ this system in an in vivo model to reveal that viral replication is not required to initiate or maintain infection within B cells.

G0S2 regulates innate immunity in Kawasaki disease via lncRNA HSD11B1-AS1

BACKGROUND

Kawasaki disease (KD) is a systemic vasculitis that is currently the most common cause of acquired heart disease in children. However, its etiology remains unknown. Long non-coding RNAs (lncRNAs) contribute to the pathophysiology of various diseases. Few studies have reported the role of lncRNAs in KD inflammation; thus, we investigated the role of lncRNA in KD inflammation.

METHODS

A total of 50 patients with KD (median age, 19 months; 29 males and 21 females) were enrolled. We conducted cap analysis gene expression sequencing to determine differentially expressed genes in monocytes of the peripheral blood of the subjects.

RESULTS

About 21 candidate lncRNA transcripts were identified. The analyses of transcriptome and gene ontology revealed that the immune system was involved in KD. Among these genes, G0/G1 switch gene 2 (G0S2) and its antisense lncRNA, HSD11B1-AS1, were upregulated during the acute phase of KD (P < 0.0001 and <0.0001, respectively). Moreover, G0S2 increased when lipopolysaccharides induced inflammation in THP-1 monocytes, and silencing of G0S2 suppressed the expression of HSD11B1-AS1 and tumor necrosis factor-α.

CONCLUSIONS

This study uncovered the crucial role of lncRNAs in innate immunity in acute KD. LncRNA may be a novel target for the diagnosis of KD.

IMPACT

This study revealed the whole aspect of the gene expression profile of monocytes of patients with Kawasaki disease (KD) using cap analysis gene expression sequencing and identified KD-specific molecules: G0/G1 switch gene 2 (G0S2) and long non-coding RNA (lncRNA) HSD11B1-AS1. We demonstrated that G0S2 and its antisense HSD11B1-AS1 were associated with inflammation of innate immunity in KD. lncRNA may be a novel key target for the diagnosis of patients with KD.

The Ends Dictate the Means: Promoter Switching in Herpesvirus Gene Expression

Herpesvirus gene expression is dynamic and complex, with distinct complements of viral genes expressed at specific times in different infection contexts. These complex patterns of viral gene expression arise in part from the integration of multiple cellular and viral signals that affect the transcription of viral genes. The use of alternative promoters provides an increased level of control, allowing different promoters to direct the transcription of the same gene in response to distinct temporal and contextual cues. While once considered rare, herpesvirus alternative promoter usage was recently found to be far more pervasive and impactful than previously thought. Here we review several examples of promoter switching in herpesviruses and discuss the functional consequences on the transcriptional and post-transcriptional regulation of viral gene expression.

Conquering the Host: Determinants of Pathogenesis Learned from Murine Gammaherpesvirus 68

Gammaherpesviruses are an important class of oncogenic pathogens that are exquisitely evolved to their respective hosts. As such, the human gammaherpesviruses Epstein-Barr virus (EBV) and Kaposi sarcoma herpesvirus (KSHV) do not naturally infect nonhuman primates or rodents. There is a clear need to fully explore mechanisms of gammaherpesvirus pathogenesis, host control, and immune evasion in the host. A gammaherpesvirus pathogen isolated from murid rodents was first reported in 1980; 40 years later, murine gammaherpesvirus 68 (MHV68, MuHV-4, γHV68) infection of laboratory mice is a well-established pathogenesis system recognized for its utility in applying state-of-the-art approaches to investigate virus-host interactions ranging from the whole host to the individual cell. Here, we highlight recent advancements in our understanding of the processes by which MHV68 colonizes the host and drives disease. Lessons that inform KSHV and EBV pathogenesis and provide future avenues for novel interventions against infection and virus-associated cancers are emphasized.

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.

Characterizing Expression and Regulation of Gamma-Herpesviral Circular RNAs

Multiple herpesviruses have been recently found to encode viral circular RNAs. Like cellular circular RNAs, these RNAs lack poly-A tails and their 5' and 3' ends have been joined, which confers protection from RNA exonucleases. We examined the expression patterns of circular RNAs from Kaposi's sarcoma herpesvirus (KSHV) in various environments. We performed deep sequencing of circRNA-enriched total RNA from a KSHV-positive patient lymph node for comparison with previous circRNA-Seq results. We found that circvIRF4 is highly expressed in the KSHV-positive patient sample relative to both B cell lines and de novo infected primary vascular and lymphatic endothelial cells (LECs). Overall, this patient sample showed a viral circRNA expression pattern more similar to the pattern from B cell lines, but we also discovered new back-spliced junctions and additional viral circular RNAs in this patient sample. We validated some of these back-spliced junctions as circular RNAs with standard assays. Differential expression patterns of circular RNAs in different cell types led us to investigate what cellular factors might be influencing the ratio of viral linear mRNAs to circular RNAs. We found that repression of certain RNA-binding proteins shifted the balance between viral linear mRNAs and circular RNAs. Taken together, examining viral circular RNA expression patterns may become useful tools for discovering their functions, the regulators of their expression, and determining the stage and cell types of infection in humans.

Th2 Cytokine Modulates Herpesvirus Reactivation in a Cell Type Specific Manner

Gammaherpesviruses, such as Epstein-Barr virus (EBV), Kaposi's sarcoma associated virus (KSHV), and murine γ-herpesvirus 68 (MHV68), establish latent infection in B cells, macrophages, and non-lymphoid cells, and can induce both lymphoid and non-lymphoid cancers. Research on these viruses has relied heavily on immortalized B cell and endothelial cell lines. Therefore, we know very little about the cell type specific regulation of virus infection. We have previously shown that treatment of MHV68-infected macrophages with the cytokine interleukin-4 (IL-4) or challenge of MHV68-infected mice with an IL-4-inducing parasite leads to virus reactivation. However, we do not know if all latent reservoirs of the virus, including B cells, reactivate the virus in response to IL-4. Here we used an in vivo approach to address the question of whether all latently infected cell types reactivate MHV68 in response to a particular stimulus. We found that IL-4 receptor expression on macrophages was required for IL-4 to induce virus reactivation, but that it was dispensable on B cells. We further demonstrated that the transcription factor, STAT6, which is downstream of the IL-4 receptor and binds virus gene 50 N4/N5 promoter in macrophages, did not bind to the virus gene 50 N4/N5 promoter in B cells. These data suggest that stimuli that promote herpesvirus reactivation may only affect latent virus in particular cell types, but not in others.Importance Herpesviruses establish life-long quiescent infections in specific cells in the body, and only reactivate to produce infectious virus when precise signals induce them to do so. The signals that induce herpesvirus reactivation are often studied only in one particular cell type infected with the virus. However, herpesviruses establish latency in multiple cell types in their hosts. Using murine gammaherpesvirus-68 (MHV68) and conditional knockout mice, we examined the cell type specificity of a particular reactivation signal, interleukin-4 (IL-4). We found that IL-4 only induced herpesvirus reactivation from macrophages, but not from B cells. This work indicates that regulation of virus latency and reactivation is cell type specific. This has important implications for therapies aimed at either promoting or inhibiting reactivation for the control or elimination of chronic viral infections.

Empirical Validation of Overlapping Virus lncRNAs and Coding Transcripts by Northern Blot

Viruses, like their metazoan hosts, have evolved to utilize intricate transcriptional mechanisms to generate a vast array of both coding and noncoding RNA transcripts. The resolution of specific noncoding RNA transcripts produced by viruses, particularly herpesviruses, presents a particularly difficult challenge due to their highly dense dsDNA genomes and their complex, overlapping, and context-dependent network of transcripts. While new long read sequencing platforms have facilitated the resolution of some noncoding transcripts from virus genomes, empirical molecular validation of transcripts from individual regions is essential. Herein, we demonstrate that the use of strand specific northern blots is essential for true validation of specific viral noncoding RNAs, and provide here a detailed molecular method for such an approach.

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.

Viral non-coding RNAs: Stealth strategies in the tug-of-war between humans and herpesviruses

Oncogenic DNA viruses establish lifelong infections in humans, and they cause cancers, often in immunocompromised patients, despite anti-viral immune surveillance targeted against viral antigens. High-throughput sequencing techniques allowed the field to identify novel viral non-coding RNAs (ncRNAs). ncRNAs are ideal factors for DNA viruses to exploit; they are non-immunogenic to T cells, thus viral ncRNAs can manipulate host cells without evoking adaptive immune responses. Viral ncRNAs may still trigger the host innate immune response, but many viruses encode decoys/inhibitors to counter-act and evade recognition. In addition, ncRNAs can be secreted to the extracellular space and influence adjacent cells to create a pro-viral microenvironment. In this review, we present recent progress in understanding interactions between oncoviruses and ncRNAs including small and long ncRNAs, microRNAs, and recently identified viral circular RNAs. In addition, potential clinical applications for ncRNA will be discussed. Extracellular ncRNAs are suggested to be diagnostic and prognostic biomarkers and, with the realization of the importance of viral ncRNAs in tumorigenesis, approaches to target critical viral ncRNAs are emerging. Further understanding of viral utilization of ncRNAs will advance anti-viral therapeutics beyond conventional medication and vaccination.

Conserved Herpesvirus Kinase ORF36 Activates B2 Retrotransposons during Murine Gammaherpesvirus Infection

Short interspersed nuclear elements (SINEs) are RNA polymerase III (RNAPIII)-transcribed, retrotransposable noncoding RNA (ncRNA) elements ubiquitously spread throughout mammalian genomes. While normally silenced in healthy somatic tissue, SINEs can be induced during infection with DNA viruses, including the model murine gammaherpesvirus 68 (MHV68). Here, we explored the mechanisms underlying MHV68 activation of SINE ncRNAs. We demonstrate that lytic MHV68 infection of B cells, macrophages, and fibroblasts leads to robust activation of the B2 family of SINEs in a cell-autonomous manner. B2 ncRNA induction requires neither host innate immune signaling factors nor involvement of the RNAPIII master regulator Maf1. However, we identified MHV68 ORF36, the conserved herpesviral kinase, as playing a key role in B2 induction during lytic infection. SINE activation is linked to ORF36 kinase activity and can also be induced by inhibition of histone deacetylases 1 and 2 (HCAC 1/2), which is one of the known ORF36 functions. Collectively, our data suggest that ORF36-mediated changes in chromatin modification contribute to B2 activation during MHV68 infection and that this activity is conserved in other herpesviral protein kinase homologs.IMPORTANCE Viral infection dramatically changes the levels of many types of RNA in a cell. In particular, certain oncogenic viruses activate expression of repetitive genes called retrotransposons, which are normally silenced due to their ability to copy and spread throughout the genome. Here, we established that infection with the gammaherpesvirus MHV68 leads to a dramatic induction of a class of noncoding retrotransposons called B2 SINEs in multiple cell types. We then explored how MHV68 activates B2 SINEs, revealing a role for the conserved herpesviral protein kinase ORF36. Both ORF36 kinase-dependent and kinase-independent functions contribute to B2 induction, perhaps through ORF36 targeting of proteins involved in controlling the accessibility of chromatin surrounding SINE loci. Understanding the features underlying induction of these elements following MHV68 infection should provide insight into core elements of SINE regulation, as well as disregulation of SINE elements associated with disease.

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.

Comprehensive annotations of human herpesvirus 6A and 6B genomes reveal novel and conserved genomic features

Human herpesvirus-6 (HHV-6) A and B are ubiquitous betaherpesviruses, infecting the majority of the human population. They encompass large genomes and our understanding of their protein coding potential is far from complete. Here, we employ ribosome-profiling and systematic transcript-analysis to experimentally define HHV-6 translation products. We identify hundreds of new open reading frames (ORFs), including upstream ORFs (uORFs) and internal ORFs (iORFs), generating a complete unbiased atlas of HHV-6 proteome. By integrating systematic data from the prototypic betaherpesvirus, human cytomegalovirus, we uncover numerous uORFs and iORFs conserved across betaherpesviruses and we show uORFs are enriched in late viral genes. We identified three highly abundant HHV-6 encoded long non-coding RNAs, one of which generates a non-polyadenylated stable intron appearing to be a conserved feature of betaherpesviruses. Overall, our work reveals the complexity of HHV-6 genomes and highlights novel features conserved between betaherpesviruses, providing a rich resource for future functional studies.