Next Article in Journal
Virus-Like Particle Vaccination Protects Nonhuman Primates from Lethal Aerosol Exposure with Marburgvirus (VLP Vaccination Protects Macaques against Aerosol Challenges)
Next Article in Special Issue
HIV-1 Mutation and Recombination Rates Are Different in Macrophages and T-cells
Previous Article in Journal
Enterovirus Control of Translation and RNA Granule Stress Responses
Previous Article in Special Issue
Characterization of Intra-Type Variants of Oncogenic Human Papillomaviruses by Next-Generation Deep Sequencing of the E6/E7 Region
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Next-Generation Sequencing in the Understanding of Kaposi’s Sarcoma-Associated Herpesvirus (KSHV) Biology

Department of Microbiology and Immunology, School of Medicine, University of Nevada, 1664 N, Virginia Street, MS 320, Reno, NV 89557, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Viruses 2016, 8(4), 92; https://doi.org/10.3390/v8040092
Submission received: 16 December 2015 / Revised: 21 March 2016 / Accepted: 23 March 2016 / Published: 31 March 2016

Abstract

:
Non-Sanger-based novel nucleic acid sequencing techniques, referred to as Next-Generation Sequencing (NGS), provide a rapid, reliable, high-throughput, and massively parallel sequencing methodology that has improved our understanding of human cancers and cancer-related viruses. NGS has become a quintessential research tool for more effective characterization of complex viral and host genomes through its ever-expanding repertoire, which consists of whole-genome sequencing, whole-transcriptome sequencing, and whole-epigenome sequencing. These new NGS platforms provide a comprehensive and systematic genome-wide analysis of genomic sequences and a full transcriptional profile at a single nucleotide resolution. When combined, these techniques help unlock the function of novel genes and the related pathways that contribute to the overall viral pathogenesis. Ongoing research in the field of virology endeavors to identify the role of various underlying mechanisms that control the regulation of the herpesvirus biphasic lifecycle in order to discover potential therapeutic targets and treatment strategies. In this review, we have complied the most recent findings about the application of NGS in Kaposi’s sarcoma-associated herpesvirus (KSHV) biology, including identification of novel genomic features and whole-genome KSHV diversities, global gene regulatory network profiling for intricate transcriptome analyses, and surveying of epigenetic marks (DNA methylation, modified histones, and chromatin remodelers) during de novo, latent, and productive KSHV infections.

1. Introduction

The most recently discovered human tumor virus, Kaposi’s sarcoma-associated herpesvirus (KSHV), or human herpesvirus 8 (HHV-8), is a γ2-lymphotropic-oncogenic virus that causes life-long persistent infection in humans. It is a major cause of AIDS-defining malignancies worldwide [1]. KSHV belongs to the Rhadinovirus genus within the Herpesviridae family, and it is categorized together with Epstein-Barr virus (EBV), murine gammaherpesvirus-68 (MHV-68), and herpesvirus saimiri (HVS) (reviewed in [2]), with a KSHV numbering system to designate open reading frames (ORFs) based on HVS homology. This human pathogen, initially identified nearly two decades ago from Kaposi’s sarcoma (KS) lesions using the representational difference analysis technique [3], is also associated with two distinct lymphoproliferative disorders: primary effusion lymphoma (PEL), or body cavity-based lymphomas (BCBLs), and multicentric Castleman’s disease (MCD)-linked plasmablastic lymphoma [4,5,6]. KSHV infection still remains a major cause of morbidity and mortality for immunosuppressed individuals, particularly organ transplant recipients and patients infected with human immunodeficiency virus (HIV) [7]. KSHV seroprevalence varies according to geographic region and ethnicity; vial infection is widespread in sub-Saharan Africa and the Amazon basin where more than half of the population is infected [8]. Lower levels of KSHV seroprevalence are reported in Northern Europe and North America with seropositivity ranging from 3% to 10% [9].
The KSHV genome is highly conserved, and it has a high degree of sequence identity across the viral strains. However, two major gene regions, including K1/VIP (a variable immunoreceptor tyrosine-based activation motif protein, encoded by the 5′ terminus of the KSHV genome) and K15/LAMP (a latency-associated membrane protein, encoded by the 3′ terminus of the KSHV genome), located at either ends of the viral genome, are highly variable compared to the central genomic region [10,11]. The sequence variability of the K1 gene has led to the identification of five major KSHV subtypes (A, B, C, D, and E), showing up to 35% variability at the amino acid level across the viral strains [10]. The sequence analysis of the K15 gene has led to the characterization of KSHV sequences, with variants designated as P, M, or N alleles, differing by up to 70% at the amino acid level [11,12]. In addition, nine distinct genomic loci (approximately 5.6% of the genome) contain additional variability, including T0.7/K12, K2, K3, ORF18/19, ORF26, K8, ORF73, and two loci within the ORF75 gene regions, within the central, more conserved region of the KSHV genome [12]. Based on the K1/K15 variability, strain classification, and variability of nine ORFs, the known KSHV genomes are classified into 12 principal genotypes.
Advancement in understanding oncovirus genomes is directly associated with the availability of molecular analytical tools, such as Northern blotting, serial analysis of gene expression, custom oligonucleotide microarrays, real-time polymerase chain reaction (PCR), and nucleic acid sequencing. In the late 1970s, two distinct DNA sequencing methods, based on either chemical cleavage of DNA or incorporation of dideoxy-nucleotides, were simultaneously reported by Gibert et al. and Sanger et al., respectively, and are known as first-generation sequencing methods [13,14]. Later, the development of fluorescently labeled dideoxy-nucleotides and automated capillary electrophoresis enabled the use of DNA sequencing or Sanger sequencing technology as part of routine laboratory and medical care [15]. Due to lower error rates, longer-read length (1000 bps obtainable sequence length), and robust performance, Sanger sequencing dominated genomics research and clinical practices until the introduction of new, cost-effective, high-throughput, and high-capacity DNA sequencing platforms were developed to more thoroughly investigate complete genomes and deliver fast, inexpensive, and accurate genetic information.
The novel DNA sequencing methodology, referred to as second-generation sequencing or next-generation sequencing (NGS), has revolutionized the field of cancer genomics and bioinformatics. The major impact of the genome technology revolution is that it has enabled the analysis of cancer genome sequences and the genome structure and the elucidation of mechanisms of cancer pathogenesis, leading to improved tumor diagnosis and treatment. These novel digital NGS methods have provided greater speed, sensitivity, and resolution at a considerably lower per-base cost over other traditional sequencing methods. The arsenal of NGS methods has already been used for whole-genome genotyping, whole-exome mapping, mutation detections, de novo assembly, and reassembly of genome, gene-expression profiling, long non-coding RNA profiling, and protein-DNA or protein-RNA binding site detection. The schematic depicting the different methods for NGS studies is presented in Figure 1. This review will summarize the recent scientific findings from the application of high-throughput technologies, and it will discuss the contributions these studies have made to unravel the intricate KSHV genomes. It will also provide valuable insights into the viral regulation of host and viral gene expression profiles during the different modes of KSHV infection.

2. Unlocking the KSHV Genome and Insights into KSHV Infection from Host Genomics: The Contributions of Whole-Genome Sequencing, Whole-Exome-Sequencing and Targeted-Sequencing

Whole-genome sequencing or DNA sequencing is used to determine the complete DNA make-up of a genome in order to gain pertinent information about genes, genetic variations, and gene functions. Massive-scale sequencing and analysis of genomic sequences provides a comprehensive characterization of complex genomes, structural variants, and a full spectrum of genomic alterations, such as single-nucleotide substitutions, small insertions and deletions, inversions, translocations, chromosomal rearrangements, and copy number alterations [16]. Recent advances in bioinformatic methods for genomic sequencing have resulted in improved variant detection and increased genomic coverage, due to the powerful combination of high-quality, short-read, and paired-end sequencing options.
The first KSHV genome was mapped in 1996, with phage and cosmid libraries from the lymphoid BC-1 cell line (GenBank ID: U75698.1) using the automated Sanger sequencing approach [17]. That study reported an almost-complete nucleotide sequence of the KSHV genome, except for the sequence for a 3-kb region at the end of the genome that could not be cloned into the sequencing vectors. The KSHV genome is reported to be a large (~165 kb) encapsidated double-stranded (ds) DNA genome consisting of a single 140.5 kb long unique coding region (LUR) with 53.5% guanine-cytosine-rich (G+C-rich) 801 bp long terminal repeat (TR) sequences at both ends of the linear viral genome. The LUR sequence has at least 81 viral open reading frames (vORFs) that encode for polypeptides larger than 100 amino acids, including 66 ORFs with homology to HVS ORFs. The study also identified a duplication of an LUR fragment in the TR region, indicative of genomic rearrangements [17]. One year later, Neipel et al. generated a nearly complete sequence of the KSHV genome (GenBank ID: U93872.2) using shotgun sequencing of fragments obtained after partial digestion of DNA extracted from a KS biopsy specimen [18]. The first extensively annotated and fully sequenced KSHV genome, GK18 (GenBank ID: AF148805.2), was obtained from the isolate of KSHV from the biopsy of classic KS lesions from a Greek patient [19].
Since viral mutagenesis is a powerful tool for examining the functions of viral genes, Gao et al. generated and sequenced a full-length recombinant KSHV genome, named KSHV BAC36, by inserting a bacterial artificial chromosome (BAC) cassette between ORF18 and ORF19 of the KSHV genome in the BCBL-1 cell line (GenBank ID: HQ404500.1), and they analyzed the function of a number of viral genes [20]. Furthermore, Yakushko et al. noted a duplication of a 9-kb LUR fragment in the TR region while sequencing the complete KSHV BAC36 genome [21]. Consequently, Jung et al. constructed a new full-length BAC clone of the KSHV genome, called BAC16, derived from the rKSHV.219 virus in JSC-1 PEL cells (GenBank ID: GQ994935.1); complete sequence analysis data indicated a minimal level of sequence variation across the entire viral genome in comparison to the other KSHV strains [22]. All of these KSHV whole-genome sequences were analyzed using the traditional Sanger sequencing. Nevertheless, they are included in this review for completeness.
The most recently sequenced KSHV strain, isolated from cells co-infected with HHV6 and known as DG-1 (GenBank ID: JQ619843.1), represents the first whole-genome KSHV sequence to be completed using Illumina next-generation sequencing of total DNA directly from patient plasma and peripheral blood mononuclear cells (PBMCs), rather than a tissue-cultured or biopsied tumor [23]. A comparison of the whole-genome sequence of this newly discovered KSHV isolate to the KSHV reference sequence (NC_009333) led to the identification of small indels in both the polymorphic regions of K1 and K15 and the TR regions. However, no novel DNA insertions, ORFs, or deletions in the entire coding regions were identified [23].
Very recently, in an effort to correlate KSHV whole-genome diversity and its potential role in KS pathogenesis, genome-scale profiling of KSHV from KS skin lesions of 16 Zambian/sub-Saharan African patients was achieved using Illumina paired-end sequencing [24]. Whole-genome analysis of these KSHV isolates, and multiple alignments of the newly sequenced KSHV genomes, showed a high level of sequence conservation but a low level of variability across the central conserved region. The observed low level variability in the otherwise highly conserved central genomic region indicates that different whole-genome KSHV variants are present in sub-Saharan Africa compared to Western countries, thus resulting in 16 new and distinct Zambian whole-genome phylogenetic structures [24].
In addition, a full genome sequence of the closest homolog of KSHV, retroperitoneal fibromatosis-associated herpesvirus Macaca nemestrina (RFHVMn), was recently analyzed from a clinical KS-like macaque tumor sample using Illumina Genome Analyzer IIx sequencing, indicating a strong genetic and sequence similarity between the two viruses [25]. The RFHVMn genome was annotated based on the KSHV genome, and it was found to have high levels of sequence conservation in specific promoter sequences and within the putative origins of replication. This suggests similarities in the biology and pathology of the two viruses. Another study by Depledge et al. used specific-targeted genomic enrichment and whole-genome profiling of 13 herpesvirus genomes, including varicella-zoster virus (VZV), EBV, and KSHV, directly from a range of low volume clinical samples, including blood, saliva, vesicle fluid, cerebrospinal fluid, and tumor cell lines [26]. That study demonstrated that target enrichment technology was suitable for isolating very low quantities of viral DNA from the pool of complex host DNA.
In addition to KSHV genome sequencing, researchers are currently also interested in identifying the host genetic variants that are important to KSHV infection and disease progression via genome-wide association studies. Because the exome represents only approximately 1% of the genome, whole-exome sequencing yields higher sequence coverage and reduces sequencing cost and time by eliminating the majority of the genome from the analysis. Exome-sequencing also facilitates the identification of possible locations or genes for tumor-causing mutations. Using the whole-exome sequencing approach in a single patient, Byun et al. identified that development of classic KS in childhood might also occur due to an inborn alteration in a gene that is important for immunity [27]. They revealed a homozygous splice-site mutation in the stromal interaction molecule 1 (STIM1) gene that is very important for immunity, and they identified STIM1 deficiency and the associated primary T-cell immunodeficiency as the genetic cause of severe KS in childhood [27]. In another report, they showed that T-cell-mediated immunity to KSHV is greatly impaired by autosomal recessive complete OX40 protein deficiency, due to R65C (Arg65Cys) loss-of-function mutation in a healthy adult with childhood-onset classic KS [28]. OX40 is a co-stimulatory receptor that is present in activated T-cells [29]. The reported data strongly suggest a causal relationship between OX40 deficiency and the development of KS infection. Another research group reported the whole-genome and whole-exome sequencing-based identification of the signal transducer and activator of the transcription factor 4 (STAT4) gene as a potential classic KS predisposing gene, and they suggested that the variation in the STAT4 gene is linked to KS/HHV-8 infection [30]. In addition to whole-exome sequencing, NGS has also been used to target either a handful of host genes or specific genomic regions to detect cancer-causing mutations. Targeted sequencing allows for rapid production of large numbers of low-cost reads on specific genomic locations at much higher coverage levels over conventional methods. In a recent study, Dittmer et al. described the X-chromosome targeted sequencing in 15 PEL cell lines in order to identify tumor-specific single nucleotide variants (SNVs) on the X chromosome [31]. They identified 34 common missense mutations in protein-coding regions in all of the PEL samples, including a Phe196Ser change in the interleukin 1 receptor-associated kinase 1 (IRAK1). This IRAK1 mutant was found to be constitutively phosphorylated and essential for PEL growth and survival. As IRAK1 binding to the MyD88 adapter protein mediates toll-like receptor (TLR) immune signaling, mutation in IRAK1 is considered to be an essential driver for KS-based lymphomas, and inhibition of IRAK1/MyD88 immune signaling might serve as a potential target for drug development [31]. Thus, the combined knowledge of high-throughput genome sequencing and host genetic and genomic variations has provided unprecedented and unbiased insights into viral and host genomes, allowing for the identification of several viral-host factors that influence the progress of KSHV.

3. An Atlas of the KSHV Transcriptome: The Contributions of RNA-Sequencing and Small RNA/Non-Coding RNA-Sequencing

3.1. RNA-Sequencing

The transcriptome is the complete set of transcripts encompassing a multitude of total coding and non-coding RNAs (ncRNAs), particularly small RNAs, including microRNAs (miRNAs) and ncRNAs, to name a few. These are transcribed in a given cell type under given physiological and pathological conditions. An analysis of the complex landscape and dynamics of the transcriptome enables identification of the full set of transcripts, which reveals the differential gene expression patterns, genome annotations, gene rearrangements, ncRNA detection, and sequence variations that are present at the genomic loci [32]. Initial transcriptomics studies using a hybridization-based microarray and traditional low-throughput sequencing-based approaches used a reference genome, and they allowed a limited dynamic range of detection, had a high background signal, and only provided the ability to detect more abundant transcripts. In contrast, high-throughput NGS, referred to as high-throughput transcriptome sequencing or RNA-sequencing (RNA-seq), provides a new method for mapping and reconstructing the entire transcriptome even without any prior knowledge of the genome [32]. The new methodologies for RNA-seq studies provide enormous sequencing depth (base-pair-level resolution), with a very low background signal, a much higher dynamic range of expression levels, and the ability to identify rare transcripts with regulatory functions [33]. Almost all RNA-seq data are generated through three commercially available NGS platforms: Illumina IG, Roche 454 Life Science, and Applied Biosystems SOLiD. Library preparation for RNA-seq requires a sample of purified, total, or fractionated RNA, which is converted to cDNA fragments with adaptors attached to one or both ends, before sequencing on a NGS platform in order to generate short sequences from one end (single-end sequencing) or both ends (pair-end sequencing) [32]. Targeted RNA sequencing, which focuses on sequencing the messenger RNA (mRNA) and small RNA (ncRNA or miRNA) targets, is achieved by including additional isolation or enrichment steps before cDNA synthesis.

3.2. Genome-Wide Analysis of the Latent and Lytic KSHV Transcriptomes

Like all herpesviruses, KSHV undergoes both latent and lytic replication cycles. In a cell culture, primary infection generally results in latency in which only a few promoters that drive the latent RNAs are active, although a small number (1% to 5%) of the latently infected cells spontaneously reactivate the lytic cycle [34,35,36,37]. Viral latent gene expression is critical for inducing and maintaining latency, and it is restricted to the complex multicistronic latency-associated transcript that encodes the viral genes for the latency-associated nuclear antigen (LANA) (ORF73), vCyclin (ORF72), and vFLIP (ORF71), along with a downstream promoter driving the transcription of the viral miRNA cluster and an ORF for the Kaposin (K12) gene [38,39,40,41,42]. After lytic induction using various environmental or chemical inducers [43,44,45,46,47], KSHV transcribes many of the lytic genes that encode the regulatory proteins that are critical for the viral life cycle. Lytic mRNAs of KSHV can be grouped, according to their timing and expression in response to the protein synthesis/DNA replication inhibitors, as immediate early (IE), early (E) and late (L) proteins that encode mRNAs [48,49,50,51]. The IE mRNAs of KSHV include 3.6-kb mRNA encoding RTA/ORF50, 1.7 kb mRNA encoding ORF45, and a 2.0 kb mRNA encoding ORF K4.2 [50]. The latent to lytic switch of KSHV is regulated by the replication and transcription activator (RTA), a 110 kDa protein, encoded by the ORF50 gene, which is capable of inducing the cascade of lytic gene expression. Interestingly, only the over-expression of the RTA protein is necessary and sufficient to disrupt KSHV latency and initiate the lytic replication cascade [52,53,54]. In the lytic cycle, the expression of E transcripts is controlled by the IE transcripts, and it includes polyadenylated nuclear RNA (PAN RNA), Kaposin, ORF57, k-bZIP (K8), K5, K9, K14, K15 ORF6, ORF21, and ORF74 [1], followed by the expression of the L transcripts, including the major capsid protein (MCP) encoded by ORF25 and the small viral capsid (sVCA). Late transcripts are transcribed following DNA replication, and they are the structural genes for virus assembly [48,49].
Transcriptome profiling of KSHV-infected cells, KS lesions, and tissue samples is fundamental to furthering our understanding of the interplay of the large number of underlying mechanisms that coordinate the regulation of KSHV viral gene expression within infected host cells. While the transcription of a few individual genes and promoters of KSHV have been well studied, little is known of the overall landscape of the viral transcriptome on the genome-wide scale [55]. Large-scale gene expression profiling of KSHV transcripts using DNA microarrays, real-time PCR arrays, and Northern blotting have revealed the dynamics and complexity of the viral genome [41,55,56]. Proteomic studies of KSHV-infected cells have assessed the expression of many of the predicted KSHV ORFs [57]. A number of NGS-based research studies have been reported in recent years providing useful insights into the complex KSHV transcriptome, necessary for interpreting the functional elements of the genome and understanding the related KS pathogenesis. During lytic KSHV infection, large segments of the genome that were initially thought to be “non-coding” are actually most likely the stable RNAs [55,58]. Recently, Arias et al. analyzed the transcription and translation profiles of the KSHV genome using a combined mRNA-seq ribosome footprinting (Ribo-seq) and genomic DNA sequencing (DNA-seq) approach to uncover new features about the genomic landscape and peptide-coding capacity of KSHV during the productive (lytic) stage of infection when the progeny virion production is underway [58]. Numerous unknown genomic and functional features were assessed using a tightly controlled and highly inducible epithelial iSLK-219 cell line containing a doxycycline inducible transgene that encodes viral RTA protein. This led to the generation of a novel revised annotation of the KSHV genome: KSHV 2.0, including a more than 45% increase in peptide coding potential with expanded illustration of 49 viral transcripts encoding peptides smaller than 100 aa, 63 new ORFs, ncRNAs, polyadenylation sites, splice junctions, and initiation/termination codons of multiple known ORFs (see Table 1) [58]. In addition, the authors of that study attributed the coding capacity of KSHV to multiple strategies, including splicing, mRNA editing, and the use of alternative transcription start sites leading to multiple upstream ORFs (uORFs) and small ORFs (sORFs). The systematic genome-wide survey of the viral transcriptional and translational activity added to our understanding of the underappreciated complexity of KSHV gene expression during distinct phases of the viral life cycle [58]. However, if the new KSHV genome annotation is conserved among multiple cell types that harbor a stably- and latently-maintained KSHV episome still needs to be elucidated.
Immediately after the de novo infection, an initial burst of lytic transcripts occurs followed by accumulation of the latent transcripts. Additionally, the KSHV viral particles contain a variety of viral mRNAs and several small RNAs, such as viral and cellular miRNAs, ncRNAs, and unusual, small RNAs (usRNAs) [59,60]. Our lab recently performed an in-depth study of the KSHV transcriptome during de novo infection of human PBMCs, CD14+ monocytes, and telomerase-immortalized vascular endothelia (TIVE) cells using Illumina RNA-seq methodology [61]. The transcriptome analysis demonstrated that a significant amount of KSHV latent and lytic transcripts accumulated quickly in all three cell types at 4 h post-induction (hpi); the concomitant silencing of latent and lytic transcripts occurred after 24 hpi for CD14+ cells and after 48 hpi for PBMCs and TIVE cells [61]. These results are indicative of a cell-type-specific KSHV gene expression pattern that occurs early during KSHV infection. During that period, the viral transcripts that were detected in significant amounts included PAN RNA, ORF58/59, kaposin B, K2, K4, K6, ORF11, ORF17, ORF45, ORF27, ORF37, ORF57, ORF64, ORF65, ORF73, and T0.7. RNA-seq analysis showed that the expression of transcripts involved in lytic DNA replication, such as ORF59, ORF9, ORF6, ORF56, ORF40, ORF54, ORF60, ORF61, ORF70, and ORF37, and transcripts encoding viral structural proteins, such as K8.1 and ORF65, occurred as soon as 4 hpi. The detailed analysis of the RNAs composition in the KSHV virions showed that a large number of viral proteins and mRNAs, including PAN RNA, ORF58, ORF59, T0.7, and ORF17, are encapsidated in the infectious virions that are effectively transduced into the target cells [61]. Furthermore, the Click-iT Nascent RNA Capture Kit was used to differentiate between the RNAs present in the virions and those transcribed during the early stage of infection. The RNA-seq of the newly synthesized RNA captured from the 5-Ethynyl-2′-deoxyuridine (EdU)-labeled KSHV-infected PBMCs showed increased expression of several viral genes that were detected at 4 hpi, including ORF2, K3, T1.5, PAN RNA, ORF26, ORF29, ORF36/37, ORF40, ORF50, K8.1, ORF55, vIRF-2, ORF58/59, ORF60, ORF63/64, ORF67, T0.7, ORF72/73, and ORF74. However, at 24 hpi the active transcription was limited to the ORF73 sequence reads. As expected, the analysis of the transcripts necessary for lytic DNA replication indicated a mixture of the transcripts that were actively transcribed and lytic genes that were transduced with the virions [61].
Several research groups have utilized deep RNA-sequencing approaches to demonstrate the role of various viral and host genes in regulating KSHV gene expression during lytic replication and reactivation. LANA, expressed by KSHV ORF73, is a major latency nuclear protein with pleiotropic functions. It plays an important role in repressing lytic genes by binding to and repressing the ORF50/RTA promoter [62]. It was recently found that LANA recruits the host transcriptional repressor, Krüppel-associated box domain-associated protein 1 (KAP1), and represses the RTA promoter [63,64]. However, the location, mechanism, and the dynamics of LANA/KAP1 binding to the RTA promoter have yet to be investigated. Chandran et al. identified the role of the host transcription factor; nuclear respiratory factor 2 (Nrf2) in modulating KSHV gene expression in latently infected PEL-derived cell lines [65]. Cellular Nrf2 plays an important role in the dynamic changes in RTA expression. In the absence of LANA, Nrf2 acts as a transcriptional activator, but it functions as a repressor in the presence of LANA. This switch in the role of Nrf2 on the RTA promoter is mediated by LANA-regulated recruitment of KAP1. Nrf2 inhibition further resulted in increased KSHV lytic gene expression, virion production, and PEL cell death. Collectively, this study demonstrated that KSHV-induced Nrf2 facilitates lytic gene expression during de novo infection, and it later represses this induction via LANA-mediated KAP1 recruitment to the Nrf2 binding site [65]. In a study by Swaminathan et al., the role of KSHV ORF57, a lytic-phase KSHV protein that is essential for lytic replication and virion production, was examined using a combination of deep RNA-seq, quantitative reverse transcription polymerase chain reaction (qRT-PCR), Northern blotting, and rapid amplification of cDNA ends (RACE) analysis [66]. To investigate the mechanism by which the ORF57 gene product enhances KSHV gene expression, a comparison of all of the KSHV lytic transcripts in the presence or absence of ORF57 was conducted. The results indicated that two clusters of the KSHV lytic cycle transcripts are expressed more abundantly in the presence of ORF57 [66]. In that study, eight transcripts that were expressed during early and late lytic replication were found to be particularly ORF57-dependent, and they are involved in DNA replication, DNA packaging, virion assembly, and immune modulation. A 2014 study by Gong et al. sought to examine the role of previously uncharacterized KSHV ORF18 and ORF30 genes during KSHV de novo lytic replication in papilloma-immortalized human oral keratinocyte (HOK16B cells) [67]. Using the KSHV ORF18- and ORF30-deficient virus, high-resolution profiling of transcripts via RNA-seq indicated that these genes are essential for efficient late gene transcription during both de novo KSHV infection and lytic reactivation. Studies of late gene transcription in beta- and gammaherpesviruses provide further evidence that ORF18 and ORF30, along with ORF24, ORF31, ORF34, and ORF66 (gammaherpesvirus nomenclature), belong to the virus-specific pre-initiation complex necessary for the expression of late viral genes [68,69,70]. In comparison to high-resolution global examination of KSHV latent and lytic transcripts, relatively fewer NGS-based studies have been conducted to determine the link between the virus and the host mRNAs. In this context, foremost studies by Glaunsinger et al. revealed that the SOX (host shutoff and alkaline exonuclease) protein encoded by the KSHV ORF37 protein interacts with host translation machinery and shuts off host gene expression in the early stage of infection by accelerating global mRNA turnover [71,72]. KSHV-induced mRNA turnover is governed by a coordinated destruction by both the KSHV exonuclease SOX and the mammalian 5′–3′ mRNA exonuclease XrnI in a two-step process where the mRNAs are initially cleaved internally by the SOX protein, followed by degradation via cellular XrnI [73]. Notably, the cytoplasmic mRNA-targeting SOX protein has been recently shown to cause cellular mRNA decay in mammalian cells by altering RNA polymerase II (RNAPII) transcription in the nucleus via a feedback mechanism that occurs in response to the catalytic activity of XrnI [74]. High-throughput RNA-seq analysis of SOX-expressing cells has indicated selective destabilization of endogenous cellular p53-induced protein with death domain (PIDD), but not Apoptosis Enhancing Nuclease [75] or interleukin-6 (IL-6) transcripts by SOX. This suggests that cellular mRNAs utilize several different mechanisms to escape SOX-mediated shutoff [76,77]. Cleavage of cellular mRNAs by the KSHV SOX protein is both target- and cleavage-site specific, as defined by the sequence of the target RNAs. In order to shed light on the targeting mechanisms of SOX, Gaglia et al. mapped the location of SOX cut sites across the human RNA transcriptome using a degradome sequencing technique, called parallel analysis of RNA ends (PARE), and targeted mutagenesis [78]. The sequences surrounding the SOX cleavage sites in the endogenous mRNA targets were found to have a degenerate sequence pattern that enables SOX to achieve a broad targeting capability while maintaining cut site specificity within the mRNA [78]. Recently, Conrad et al. performed a high-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (HITS-CLIP), to analyze ORF57-cellular RNA interactions in KSHV-infected cells during lytic reactivation [79]. They identified ORF57 interactions with both the virus and the host RNAs throughout the KSHV genome and at the KSHV origins of lytic replication.
In addition to studying viral and cellular gene expression profiles during persistent and/or productive KSHV infection, NGS-based transcriptome studies have been recently extended to discover new virus-tumor associations and the underlying mechanism that drives infection and causes virus-mediated oncogenesis. In a recent large-scale investigation of 50 cancer cell lines, Cao et al. used high-throughput RNA-seq to detect the presence of different viruses that are known to infect humans/mammalian cells and to study virus-host interactions [80]. Along with other significant discoveries, including the detection of a new unknown murine leukemia virus infection and the identification of cellular cytokine regulators disrupted by virus integration, global assessment of the KSHV transcriptome in the PEL-derived BCP-1 cell line indicated that viral immune signaling genes in the left cluster of the genome, such as K2/vIL-6, K4/vIL-8, and K5/E2-ubiquitin ligase 1, showed higher expression values of 840 FPKMs (fragments per kilobase per million mapped reads), 145 FPKMs, and 351 FPKMs, respectively, when compared to a LANA-containing rightward cluster with an expression value of 39 FPKMs; this suggests that secreted cytokines have a significant impact on both latent and lytic KSHV-infected cells [80].

3.3. Small RNA/Non-Coding RNA-Sequencing

3.3.1. Role of miRNAs in KSHV Latency and Lytic Reactivation

Two classes of small (~21–25 nucleotide (nt)) ncRNAs that regulate genes and genomes are: miRNAs and short interfering RNAs (siRNAs) [81]. The miRNAs bind to target messenger RNAs (mRNAs) and cause mRNA translation inhibition or degradation by targeting complimentary sequences in the 3′ untranslated regions (3′ UTR) of the mRNAs [82,83]. These regulatory RNAs play a pivotal role in a variety of biological processes, including viral infection, and dysregulation of miRNAs promotes disease progression and viral cancer pathogenesis [84]. In addition to large numbers of cellular miRNAs, bioinformatics, sequencing, and direct cloning methods have been used to identify and characterize several miRNAs of viral origin [85]. Like all herpesviruses, it is now evident that KSHV encodes many miRNAs in latently and/or productively infected cells (reviewed in detail in [83,84,86,87]). Although only a limited number of protein-coding genes are expressed during latency, KSHV expresses 25 mature miRNAs (at last count) derived from 12 precursor miRNAs (pre-miRNAs) that are located in one cluster in the major latency-associated locus (118–128 kb); it also modulates cellular miRNAs as a means of affecting cell development and cell cycle regulation [88]. The mature KSHV miRNAs are annotated as miR-K1 to miR-K12, depending on their proximity to the Kaposin gene. The KSHV miRNAs are constitutively expressed in the PEL cell lines and de novo infection models during the latent as well as lytic modes of viral infection. The expression levels of KSHV miRNAs in vitro is shown to be controlled by both latent and inducible lytic promoters, although the expression levels of these miRNAs during the natural KSHV infection is yet to be identified [40,89,90]. Furthermore, KSHV miRNAs sequences are highly conserved (approximately 99.6%) among several PEL cell lines and KSHV isolates [91].
Eleven of the 12 currently known distinct KSHV miRNAs, miR-K1 to miR-K11, were initially identified by cDNA cloning and traditional sequencing in latently infected B-cells, while miR-K12 was predicted computationally followed by microarray and Northern blot analysis [92]. Sequence analysis using small RNA high-throughput sequencing on the Solexa/Illumina and ABI SOLiD platforms has offered further insights into the miRNA coding potential of KSHV. The first, in-depth analysis of small viral miRNA-expression in the latently KSHV-infected BC-3 cell line by Umbach et al. demonstrated that KSHV-encoded miRNAs account for approximately 92% of the total miRNA cDNA reads [93]. Significant discrepancies were revealed between the primary obtained viral sequences for miR-K8 and miR-K12, in comparison to the sequences submitted in miRBase [90,92]. The miR-K9 stem-loop located within an approximately 1.7 kb sequence that represents the origin of lytic genome replication, ori-LytB [94,95], appears to be highly mutated in the KSHV-infected BC-3 cell line, thus affecting the miR-K9 expression. The study also reported the first detection of several KSHV-miRNA-offset RNAs (moRNAs) in mammalian somatic cells, for nine of the 12 KSHV miRNAs expressed from the KSHV latency region, although the moRNAs were not abundant [93]. As reported earlier by Pfeffer et al., two distinct seed sequences, miR-K10a and miR-K10b, were detected for miR-K10; these isoforms differ by one nucleotide due to the differential processing at the 5′ end of miR-K10 [89].
In another study, Lin et al. used deep sequencing of small RNA libraries to identify a total of 25 different 5p- (derived from the 5′ arm) and 3p- (derived from the 3′ arm) derivatives of mature KSHV miRNAs obtained from all 12 pre-miRNAs from cells undergoing lytic infection [96]. This study also reported the detection of several novel small viral moRNAs and antisense viral miRNAs (miRNA-AS). KSHV miRNAs, such as miR-K1, miR-K3, miR-K4-3p, miR-K6-3p, and miR-K11 are expressed at detectable levels, unlike the miR-K9 stem-loop, which seemed to have a very low or undetectable expression level; this finding was consistent with earlier reports [86,90,93]. This differential expression of KSHV miRNAs could be attributed to either the differences in the processing and stability of miRNAs or to their incorporation into the RNA-induced-silencing complex (RISC).
Another interesting aspect of KSHV miRNA sequences that might shed light on KSHV virology stems from the observation that several viral and cellular miRNAs as well as usRNAs have been detected in the KSHV virions via Illumina HiSeq 2000 sequencing and in situ hybridization-electron microscopy (ISH-EM) [59]. Numerous KSHV miRNAs, and a subset of poorly expressed host miRNAs, are preferentially packaged into KSHV virions during virus production, and they have been hypothesized to play a role in viral maturation [97,98]. In fact, these virional miRNAs are able to enter host cells where they function to repress gene expression during de novo infection. Given that almost all of the herpesviruses, polyomaviruses, and adenoviruses have been shown to encode their own miRNAs, it is possible that these DNA virus particles also consist of viral-encoded miRNAs. In addition, a significant number of usRNAs have a sequence that is similar to has-miR-K1246, and these are, possibly, processed from a U2 snRNA or are transcribed from their own promoter in the KSHV virions [59]. It is believed that these virional miRNAs and usRNAs might help the virus evade the host immune system and trigger KSHV latency during de novo infection; however, their exact role remains unclear.
Identifying the mRNA targets and functional significance of KSHV miRNAs is critical for understanding KSHV infection and oncogenesis. KSHV miRNAs are known to regulate the viral life cycle, the cell cycle, apoptosis, angiogenesis, and host immune surveillance [99,100]. Recently, ribonomic approaches, such as HITS-CLIP and PAR-CLIP (photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation) have been used to identify the positions of protein-KSHV miRNA interactions with a higher resolution and specificity [101,102]. In addition to encoding KSHV miRNAs, KSHV is also known to regulate the expression levels and functions of numerous cellular miRNAs. Using PAR-CLIP, more than 2000 cellular mRNA targets of KSHV miRNAs have been recently identified, including those that are likely to influence KSHV pathogenesis [103]. Moreover, Ago HITS-CLIP, a technique that combines UV cross-linking, immunoprecipitation of Ago-miRNA-mRNA complexes, and deep sequencing, identified 1170 and 950 putative cellular KSHV miRNA targets in two-KSHV-latently-affected PEL cell lines, BCBL-1 and BC-3, respectively [102]. Significant differences between these cell lines were observed in the Ago-association of KSHV and host miRNAs. This suggests that significant differences might exist in their respective miRNA targetomes.
Small RNA-sequencing and poly-A-enriched mRNA sequencing of viral and cellular miRNAs in latently KSHV-infected cells identified differential expression of miRNA-mRNA target pairs [104]. That study revealed 153 differentially expressed human miRNAs, most of which were significantly downregulated, including the large 14q32 miRNA cluster that is shown to be deregulated in human diseases and cancer. Correlation of miRNA expression levels to their respective mRNA target genes using TaqMan assay indicated up-regulation of miR-708-5p, leading to the down-regulation of caspase-2 and leukemia inhibitory factor (LIF), as well as down-regulation of miR-409-5p, a miRNA associated with an increase in the p53-inhibitor MDM2 [104].
The miRNAs that bind to the 3′UTR of the target transcripts play an important role in the posttranslational regulation of a gene expression program [81,82]. The 3′UTR refers to the portion of an mRNA transcript that extends from the stop codon to the polyadenylated tail, and it affects transcript stability, subcellular localization, and translation efficiency [105,106]. To date, KSHV 3′UTRs have been determined for more than 75% of protein-coding RNAs [58,107,108]. In the past few years, several research groups have performed genome-wide mapping and screening of KSHV transcripts and 3′UTRs to identify numerous bicistronic and polycistronic transcripts in several KSHV genomic loci as novel targets of KSHV miRNAs [58,107,109,110]. The bicistronic and polycistronic transcripts have been identified at numerous KSHV genome loci, including ORF30–33, ORF34-37, ORF50-K8, and ORF71-73 [38,50,111,112]. Knowledge of a likely role for KSHV 3′UTRs in the control of the KSHV life cycle has further increased since Sullivan et al. suggested that some KSHV 3′UTRs might differentially regulate KSHV gene expression during latent vs. lytic infection [107]. Mapping of the 3′UTRome of major KSHV latent and lytic transcripts from TREx BCBL1-RTA cells using Illumina RNA-seq analysis, 3′ rapid amplification of cDNA ends (3′RACE) revealed KSHV transcripts with common polyadenylation sites, leading to shared 3′UTRs and similar 3′UTR-mediated regulation. Approximately half of the KSHV 3′UTRs, including those with higher levels of guanine-cytosine (GC), showed a decrease in gene expression during latent infection. In addition, that study indicated that some KSHV 3′UTRs are sufficient to impart increased gene expression during lytic infection [107]. A study by Gao et al. indicated that KSHV miRNAs frequently target the extended 3′UTR of the 5′-proximal ORFs and 5′-distal ORFs of bicistronic and polycistronic transcripts [109]. Screening of the 3′UTR sequences and polyadenylation sites of 74 KSHV genes has identified that 11 out of 28 KSHV miRNAs targets are bicistronic or polycistronic transcripts that are spread across the entire KSHV genome. Reporter mutagenesis suggests that miR-K3 specifically targets the locus that encodes for lytic genes ORF31–33, thereby resulting in the inhibition of their protein expression. On the other hand, the KSHV isoforms miR-K10a-3p and miR-K10b-3p, and their variants, repressed ORF71, ORF72, and ORF73 transcripts in the latent locus through binding sites in the 5′-distal ORFs and intergenic regions. In addition, a variety of KSHV-encoded miRNAs have been shown to directly target some lytic genes, including ORF50 (miR-K9* and miR-K12-7-5p), ORF56/57 (miR-K5 and miR-K6-3p), and ORF-K2 (miR-K10a-3p), and modulate the latent-to-lytic phase transition [102,113,114,115].

3.3.2. Role of ncRNAs in KSHV Lytic Reactivation

The impact of NGS technologies on ncRNAs discovery and characterization is particularly noteworthy. Because this topic has been reviewed extensively elsewhere [116,117], we only provide a brief summary here.
During lytic infection, one of the most abundant viral transcripts present in KSHV-infected cells is a 1.1 kb, long non-coding PAN RNA [118]. Originally called T1.1 or nut-1 RNA, this early lytic gene product plays a major role in the regulation of viral and host transcriptional programs during all phases of the KSHV life cycle [119,120,121,122]. The PAN RNA locus is located between 28,661 and 29,741 nucleotides and between K6 and ORF16 within the KSHV genome [120]. The high expression of PAN RNA is regulated by the lytic switch K-RTA protein that transcriptionally activates PAN RNA by binding to the closely related response elements, and by ORF57 or mRNA transcript accumulation (Mta), which binds and stabilizes PAN RNA from the cellular RNA decay pathways [123,124,125,126]. The ORF57 interaction and stability of PAN RNA depends on the 9-nucleotide core in the Mta-responsive element (MRE-II) in the 5′ PAN that also binds cellular PABPC1 [121,125]. A recent report by Arias et al. utilized mRNA-seq, along with Ribo-seq and genomic DNA-seq, to provide further knowledge about the coding potential of PAN RNA. By applying these techniques in parallel, that research group demonstrated that PAN RNA accounts for 78%–90% of the KSHV reads as early as 24 hpi [58]. Interestingly, the Ribo-seq analysis of KSHV-infected cells revealed the presence of initiating and elongating ribosomes on the “non-coding” nuclear PAN RNA that protects the PAN RNA during the lytic cycle; this indicates that either a fraction of PAN RNA is cytoplasmic or ribosomes can access the PAN RNA in the nucleus for translation. Thus, ribosome-bound RNA might regulate mRNA encoding for several putative peptides, which, if stable, could be very abundant [58]. In addition, this study identified three PAN-encoded peptides, PAN1.1, PAN1.2, and PAN1.3, which are indicative of the increased coding capacity of PAN. Rossetto et al. used chromatin isolation by RNA purification coupled with next generation sequencing (ChIRP-seq) approach [127], an invaluable tool for mapping specific lncRNA–DNA interactions, and showed that PAN RNA binds and interacts with the KSHV genome to initiate the viral lytic-phase transcription program and occupies much of the KSHV genome during lytic infection, including its own promoter [119]. PAN RNA interacts with specific demethylases, UTX and JMJD3, within the KSHV genome and with the protein components of polycomb repressive complex 2 (PRC2) to globally influence viral and cellular gene expression, respectively [120]. Recently, another non-coding PAN RNA of 3.0 kb (designated as T3.0), transcribed from the opposite strand of the KSHV ORF50, has been identified and reported to encode a small peptide, designated as viral small peptide 1 (vSP-1, 48aa); vSP-1 complexes with RTA proteins and inhibits the proteasome-mediated degradation of RTA, thereby facilitating viral gene expression and lytic replication [128].

4. Whole-KSHV Epigenome-Mapping: The Contributions of DNA Methylation-Sequencing, ChIP-Sequencing and FAIRE-Sequencing

The term “epigenome”, meaning “sequence-independent heritable changes of the genome”, refers to a set of reversible biochemical modifications made to both DNA and associated histone proteins, i.e., the building blocks of nucleosomes that do not involve a change in the underlying nucleotide sequence [129]. Examples of these alterations, generally referred to as epigenetic marks or modifications, include DNA methylation, histone modifications, chromatin remodeling, nucleosome occupancy, and coding and non-coding RNA expression, which collectively govern the gene regulatory network involved in several biological processes [130,131,132,133]. In terms of mapping the epigenetic marks and several DNA–protein interactions, a systematic profiling of the epigenomes in various cell types and stages can advance our understanding of how chromatin modifications are patterned genome-wide, and it can provide a clearer picture of transcriptional gene regulation [134].
Tremendous progress in NGS-based technologies has enabled the reproducible assessment of the viral and cellular epigenetic marks across the genome of multiple cell types. Since NGS allows for the unbiased, accurate, high-resolution, and substantially improved sequencing of global and specialized DNA regions, it is considered to be a cost-effective screening tool to evaluate the epigenetic architecture of the KSHV genome in comparison to traditional sequencing approaches. A more precise mapping of protein-binding sites, transcription factors, enhancers, histone variants, and nucleosome positions in the KSHV epigenome has been achieved using a variety of novel sequencing-based assays for profiling of methylated DNA (MeDIP-sequencing and MBD-sequencing) [135], the characterization of DNA-binding proteins (ChIP-sequencing) [136,137], the sequencing of accessible chromatin, via formaldehyde-assisted isolation of regulatory elements (FAIRE-seq) [138], and the discovery of positioned nucleosomes (MNase-seq), at a single-base resolution with the goal of providing new insights into many KSHV-associated diseases.
DNA methylation refers to the reversible methylation of the fifth carbon atom of cytosine in 5′–C–phosphate–G–3′ (CpG) dinucleotides by methyltransferases [139]. Additionally, 5-hydroxymethylation of cytosine (5hmC), catalyzed by the ten-eleven translocation (TET) family of proteins, has also been reported for mammalian cells [140]. Cytosine methylation is believed to play an important role in epigenetic gene regulation, and it is linked to gene silencing by either interfering with the binding of transcription factors or by modifying the chromatin structure [141]. The following assays are commonly used to identify methylated CpG dinucleotides in a genome: methylated DNA immunoprecipitation coupled with sequencing (MeDIP-seq) and methyl CpG-binding domain (MBD) protein sequencing (MBD-seq). These techniques involve capturing methylated DNA fragments using an antibody raised against 5-methylcytosine, followed by library construction and massive parallel sequencing of those fragments to identify the genomic regions of DNA methylation [142]. Histone marks are frequently investigated by chromatin immunoprecipitation followed by sequencing (ChIP-seq), a technique for characterizing and assaying various DNA-binding proteins [134]. In ChIP-seq, the DNA associated with a specific histone modification is immunoprecipitated (using an antibody specific to the protein of interest), and it is used to generate sequencing libraries for genome-wide profiling by NGS [134]. In addition, the chromosome conformation capture (3C) method, a high-throughput technique used to study the interactions between DNA sequences on the same chromosome or between different chromosomes, has been further developed from its original configuration into the carbon-copy chromosome conformation capture (5C) method, which allows for parallel analysis of physical interactions among different chromosomal loci and higher order chromosomal architecture in a high-throughput manner [143,144]. Another important technique, FAIRE-seq, has been coupled with NGS to identify regions of nucleosome-depleted DNA that are generally associated with active transcription [145]. On the other hand, because nucleosome-associated DNA is not affected by microccocal nuclease (MNase) digestion, it is characterized through MNAse-seq to precisely reveal the position and occupancy of nucleosomes in the functional regulatory regions of chromatin [146]. In recent years, scientific researchers have exploited these genomic assays to generate high-resolution epigenome maps to explain the epigenetics involved in KSHV-associated malignancies. In this section of the review, we briefly describe recent genome-wide elucidations of a broad array of chromatin phenomena that have already begun to change our view of the extent and complexity of KSHV epigenomes and their impact on persistent and productive KSHV infections.

4.1. Sequencing of the KSHV Epigenome during de Novo Infection and Latency

As with other gammaherpesviruses, KSHV virions encapsulate linear double-stranded DNA (dsDNA) in an epigenetically naive configuration without any histones or 5′-methyl cytosine residues. [147] However, upon de novo infection, when the virions enter the target cells, the linear dsDNA is circularized to generate a closed-circular DNA, which is maintained as a circular episome and tightly packed as a nucleosome in the infected host cells. The circular episome is further chromatinized for the regulation of KSHV’s biphasic lifecycle, and it is greatly influenced by the same chromatin modulation mechanisms as cellular chromatin. Thus, the chromatinization of the KSHV genome is the key regulatory step for the establishment of tightly repressed quiescent or latent infection.
Multiple insightful studies have uncovered the different spatiotemporal epigenetic changes that occur on the KSHV genome during pre-latency, latency, and lytic reactivation in distinct cell types (reviewed in [147,148]. During the pre-latency or early stages of KSHV infection in endothelial cells, i.e., within 1–8 hpi, a transcriptionally active chromatin, referred to as euchromatin, is present on the KSHV genome, and it is characterized by enhanced levels of the H3K4me3 and H3K27ac histone marks, accompanied by a temporary induction of a few IE lytic genes, including the RTA promoter [135]. However, at 24–72 hpi, the KSHV genome undergoes a transition to heterochromatin, during which the levels of the activating histone marks and the transient expression of the lytic genes decline, with a concomitant increase in the levels of the H3K27me3/H2AK119ub repressive histone marks on the KSHV genome. This facilitates the establishment of latency [135,149]. Thus, the RTA promoter exhibits a bivalent promoter structure, defined by the simultaneous presence of the H3K4me3 and H3K27me3 regions, and it regulates the different stages of the KSHV life cycle. In addition, during the early stages of de novo KSHV infection of PBMCs, the K8, ORF49, and ORF64 gene promoters have also been shown to have bivalent chromatin domains [150]. In contrast, the KSHV genome displays a transcriptionally permissive chromatin form in oral epithelial cells during primary infection that results in prolonged and robust lytic gene expression [149]. KSHV latency is also governed by other epigenetic alterations, such as cytosine 5-methylation of the CpG islands and positioning of the nucleosomes on the viral episome [151]. Recently, it has been shown that the KSHV genome in latently infected cells is heavily CpG-methylated across the entire genome, except for the LANA promoter region [135].
Dittmer et al. performed FAIRE-Seq to analyze the open chromatin organization and patterns of nucleosome depletion in KSHV latently infected B-lymphoma and epithelial cell lines with a stable and latently maintained KSHV epigenome during latency [138]. They found that approximately 8% of the KSHV genome was associated with the open chromatin regions, while a majority of the viral promoters and genes had closed-chromatin conformation. In general, the quiescent transcription regions or open chromatin regions identified by FAIRE-seq were devoid of nucleosomes and enriched in activating histone marks (H3K4me3 and H3K9/K14-ac). The nucleosome depletion patterns were found to be largely conserved among the latently-infected cells with CTCF-binding sites, coinciding with many nucleosome-depleted loci, excluding the LANA, vIL6, and proximal Kaposin promoter regions, thus, suggesting a common viral genome landscape in all forms of latency [138].
During latency, KSHV maintains a tight control of viral gene transcription, and only a few viral genes that are essential to latent DNA replication and host immune modulations are expressed. Among the KSHV latency genes, the master regulator of latent phase, LANA, which has been reviewed extensively elsewhere [152,153,154], is consistently expressed very early after de novo infection in various KSHV-harboring cells and tissues. LANA is crucial for maintaining the KSHV episome by direct and indirect DNA-binding of LANA to the host and viral chromosome [155,156]. Lu et al. performed two independent LANA-specific ChIP-seq experiments on latently infected BCBL-1 cells, and they identified a total of 256 LANA-binding sites on the host genome [137]. As might be expected, LANA predominantly binds at the TR region of the KSHV genome, especially at the high-affinity LANA-binding sites, i.e., LBS1/2, and indirectly to the LANA promoter. Notably, some of the LANA-binding sites on the cellular genome were identified as the transcription start sites (TSSs) of numerous genes that are associated with the p53 pathway, and the proximal promoter regions of the genes linked to tumor necrosis factor (TNF) network, both of which are pivotal for maintenance of the KSHV genome [137]. Very recently, another research group investigated LANA binding to multiple active viral and cellular genes by performing a detailed ChIP-seq analysis for LANA, polymerase II (Pol II), and histone marks, i.e., H3K4me3 and H3K27me3 occupancy on the KSHV epigenome in BCBL-1 cells [157]. They also proved that LANA binds to chromatin near TSSs for selected viral genes, confirming earlier published data [136,137]. Promoters for LANA, vIRF-3, and vIL-6a-c displayed a correlation between enrichment of activating H3K4me3 marks and Pol II occupancy (but they were depleted for repressive H3K27me3 marks). Instead, H3K27me3 enrichment was primarily observed in heterochromatic regions, such as lytic late gene promoters, and it was comparatively enhanced in latently infected TIVE (TIVE-LTC) cells in comparison to BCBL-1 cells [157]. The coexistence of LANA with H3K4me3 marks on active host gene promoters prompted researchers to investigate the interaction of LANA with H3K4me3 methyltransferases. The results suggest that LANA associates with H3K4me3 methyltransferase hSET1 complex to create activating histone marks, and it prevents H3K27me3 silencing of latency-associated genes to maintain latency [157].
Integrative analysis of ChIP-seq and RNA-seq data provides a deeper understanding of the interplay between the KSHV’s transcriptome and the epigenome. Recently, Mercier et al. used a combined ChIP-seq and RNA-seq approach to study the global relationship between LANA’s occupancy on the cellular and viral genome and its impact on KSHV-driven host gene expression during latency [136]. They found that LANA preferentially binds the host genes near TSSs, and it is highly abundant within transcriptionally active promoters with a sequence identical to the LBS1 motif of KSHV DNA. LANA-bound promoters are enriched in H3K4me3 histone marks, and only a small percentage of the latent genes are differentially regulated in KSHV-infected cells, suggesting that direct LANA-promoter binding is not solely responsible for the altered host gene expression [136]. They also examined the levels of chromatin-bound LANA during lytic induction and found that the binding of LANA to the host and viral chromatin is severely disrupted after 72 hpi in cells undergoing lytic reactivation. This is possibly due to the remodeling of chromatin structure, although the exact mechanism is yet to be identified [136].
As previously mentioned, there is a short period of lytic gene expression during the early stages of KSHV primary infection [158]. Lytic gene expression and the expression of KSHV-encoded LANA protein are quickly shut down, once again suggesting that LANA plays an important role in the establishment of latency. Additionally, LANA is known to down-regulate RTA expression by repressing the RTA promoter [63,159,160]. Two independent research groups recently identified a host protein, KAP1 as a new regulator of KSHV latency that binds to LANA and a significant number of lytic promoters during latency [161,162]. A study by Sun et al. demonstrated the plausible mechanism of transcriptional repression of the RTA promoter by LANA-recruited KAP1, during de novo infection [63]. Genome-wide mapping of the KSHV genome using ChIP-seq detected multiple co-occupation sites for LANA and host KAP1 on the whole KSHV genome, indicating a direct LANA-KAP1 interaction. The results suggest that LANA interacts with both the N- and C-terminal domains of KAP1 in the nucleus, and it recruits KAP1 to the KSHV genome in order to shut down the lytic gene expression and facilitate the establishment of KSHV latency [63].

4.2. Sequencing of the KSHV and Host Epigenomes during Lytic Reactivation and Replication

Several physiological and environmental factors can occasionally trigger the KSHV to transition from the latent to the lytic phase of replication, followed by a complex ordered expression of lytic genes and the release of infectious virus [163]. This dynamic and ordered pattern of viral gene expression is shown to be temporally and epigenetically regulated [135]. The KSHV IE-lytic protein ORF50/RTA functions as the master regulator of KSHV lytic reactivation [53,54,164]. RTA activates the transcription of numerous downstream target genes, either by direct sequence-specific binding to the RTA-responsive elements or by indirect protein–protein interactions with cellular transcriptional factors [165,166,167]. Using genome-wide KSHV genomic microarrays, two independent research groups identified a set of 19 RTA-binding sites located in the promoter, intron, or exon regions in the KSHV genome [168,169]. Two studies have characterized the global presence of histone modifications during reactivation using ChIP-on-chip tiling arrays; they found that activating marks, such as H3K9ac and H3K4me3, are enriched at several loci, including the latency-associated region and the RTA gene promoters, during de novo infection, while repressive H3K9me3 and H3K27me3 marks are developed on the RTA promoter, following primary infection [135,170]. A bivalent chromatin structure of an RTA promoter signifies the complexity of the reactivation process that is controlled by a delicate balance between these activating and repressive histone marks. Upon reactivation, PAN RNA recruits H3K4me3 histone methyltransferase and H3K27 histone demethylases to release transcriptionally repressive marks and up-regulate RTA expression [120,171]. In addition, through association with cellular transcription factor CBF1, RTA binds to its own promoter and recruits histone acetyltransferases to further remodel chromatin, thus favoring reactivation [170].
Host epigenetic regulatory mechanisms are also significant for controlling KSHV lytic reactivation. Thus, high-resolution NGS of a host epigenome is fundamental for studying the epigenetics of virus-host interactions. Several host proteins that regulate chromatin organization and RNA transcription bind KSHV genomes at specific sites and play an important role in controlling viral reactivation [172,173]. CTCF and Rad21 have recently been investigated as potent host cell restriction factors that inhibit KSHV lytic transcription [174]. These proteins are known to bind at distinct sites on herpesvirus genomes during latency, although their exact role in lytic induction is still unknown. In an attempt to unravel the mechanisms of viral gene regulation by host proteins, CTCF, and cohesion during KSHV replication, Li et al. conducted high-resolution ChIP-seq experiments and examined CTCF and cohesin binding to the viral chromatin. CTCF and cohesin act as global negative regulators of transcription, as their dissociation from the viral chromatin during reactivation leads to a 100-fold increase in virus yield [174]. According to the RNA-seq data, these proteins appear to initially exert an activating effect on transcription; however, their prolonged presence significantly inhibits KSHV lytic transcription and the accumulation of viral RNA. Deep sequencing of host epigenome further indicates that KSHV directly or indirectly modulates the host cell’s Small Ubiquitin-like Modifier (SUMO) system in order to overcome the host’s antiviral immune surveillance that limits KSHV lytic reactivation [175]. SUMOylation is a reversible post-translational modification, and three major SUMO protein isoforms—SUMO-1, SUMO-2, and SUMO-3—have been associated with the epigenetic transcriptional regulation of the chromatin states [176]. Initial studies have identified KSHV-encoded K-bZIP as an enzyme with SUMO E3 ligase activity and specificity towards SUMO-2 and SUMO-3 [177]. Further studies by Chang et al. investigated differential genome-wide SUMOylation patterns to identify the role of chromatin SUMOylation during KSHV reactivation [175]. SUMO-1-specific or SUMO-2/3-specific ChIP-seq data has revealed SUMO-2/3 enrichment at the cellular promoter regions of highly transcribed genes and higher K-bZIP-binding on SUMO-2/3 enriched regions. Interestingly, most of the SUMO-2/3-enriched target cellular genes are involved in cell growth, apoptosis, immune response, and cytokine signaling [175]. These studies indicated that SUMO-2/3 marks control the consistent expression of host immune response genes during lytic reactivation. The study was recently extended to demonstrate SUMO-1 and SUMO-2/3 specific chromatin enrichment on KSHV genome during the productive reactivation using chromatin immune-precipitation in conjunction with deep sequencing. An inverse correlation between the SUMO-2/3 marks and the H3K9me3 histone marks was identified. A significant increase in the virus production in SUMO-2/3 knockdown or inactivated viral K-bZIP mutant during reactivation also suggested that SUMO-2/3 modifications on KSHV chromatin could negatively regulate KSHV lytic gene expression and viral reactivation [178].

5. Conclusions

The advent of massively parallel sequencing technologies that are still in the early stages of use have revolutionized the field of KSHV genomics, allowing for new biological insights on genomic, transcriptomic, and epigenomic variations. Although sequencing costs have fallen dramatically, NGS data analysis could be difficult in species with complex genomes that are hundreds of kilobases long and contain direct or inverted repeats. This poses significant challenges during the process of sequence assembly. So far, most herpesvirus genome sequencing data have been collected using large amounts of DNA obtained through culturing the virus in permissive cell lines. In uncultured or unenriched samples, a relatively lower abundance of viral DNA than host DNA leads to a lower number of viral reads; thus, it is difficult to assemble the insufficient viral reads into accurate contiguous sequences (contigs) using de novo algorithms. As a result, only a few reported sequence characterizations of KSHV genomes are available in the public databases.
Though RNA-Seq can give an unbiased view of KSHV transcriptome without the prior knowledge about each sequence of a transcript, the intricate transcriptome is always more complicated than predicted. Thus, analysis of RNA-Seq data is not a straightforward task, as we still do not know the whole potential of the viral transcriptome. For instance, there might be many viral transcripts that have not been annotated yet. Several other issues that complicate the accurate gene expression profiling using RNA-Seq, include, the presence of multiple related sequences within the KSHV genome, overlapping transcripts from different genes and precise estimation of abundance of different transcripts from the same gene. Further complicating the matter, the current existing databases (e.g., Ensembl) used to interpret RNA-Seq data use the outdated genomic assemblies, hg19 and mm9 (for human and mouse, respectively) and gene annotations (general feature format; GFF files) via the popular pipeline for RNA-Seq analysis TopHat and Cufflinks. This is problematic as significantly increased number of total mapped reads can be captured from each RNA-Seq data with the newest genomic assembly hg38, to give a comprehensive view of transcriptome. Nevertheless, thanks to the rapidly decreasing cost of sequencing, higher resolution and unbiased genome-wide characterization of KSHV genomes can be achieved at a single base-pair scale. Thus, in coming years, if the third generation sequencing technologies that are currently being developed deliver as promised, extensive sequencing data collection using user-friendly and robust software tools will continue to provide laboratory scientists and bioinformaticians with a wealth of additional information to unravel the three-dimensional (3D) KSHV genome for multiple applications.

Acknowledgments

We sincerely apologize to the authors whose important work could not be cited here owing to space limitations. The authors also thank the Verma Lab members for critical review of the manuscript. This work was supported from the National Institute of Health (CA174459 and AI105000) and the Research Scholar Grant (124389-RSG-13-230-01-MPC) from the American Cancer Society.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dourmishev, L.A.; Dourmishev, A.L.; Palmeri, D.; Schwartz, R.A.; Lukac, D.M. Molecular genetics of kaposi’s sarcoma-associated herpesvirus (human herpesvirus 8) epidemiology and pathogenesis. Microbiol. Mol. Biol. Rev. 2003, 67, 175–212. [Google Scholar] [CrossRef] [PubMed]
  2. Verma, S.C.; Robertson, E.S. Molecular biology and pathogenesis of kaposi sarcoma-associated herpesvirus. FEMS Microbiol. Lett. 2003, 222, 155–163. [Google Scholar] [CrossRef]
  3. Sarda, S.; Hannenhalli, S. Next-generation sequencing and epigenomics research: A hammer in search of nails. Genom. Inform. 2014, 12, 2–11. [Google Scholar] [CrossRef] [PubMed]
  4. Soulier, J.; Grollet, L.; Oksenhendler, E.; Cacoub, P.; Cazals-Hatem, D.; Babinet, P.; d’Agay, M.F.; Clauvel, J.P.; Raphael, M.; Degos, L.; et al. Kaposi’s sarcoma-associated herpesvirus-like DNA sequences in multicentric castleman’s disease. Blood 1995, 86, 1276–1280. [Google Scholar] [PubMed]
  5. Cesarman, E.; Chang, Y.; Moore, P.S.; Said, J.W.; Knowles, D.M. Kaposi’s sarcoma-associated herpesvirus-like DNA sequences in aids-related body-cavity-based lymphomas. N. Engl. J. Med. 1995, 332, 1186–1191. [Google Scholar] [CrossRef] [PubMed]
  6. Chang, Y.; Cesarman, E.; Pessin, M.S.; Lee, F.; Culpepper, J.; Knowles, D.M.; Moore, P.S. Identification of herpesvirus-like DNA sequences in aids-associated kaposi’s sarcoma. Science 1994, 266, 1865–1869. [Google Scholar] [CrossRef] [PubMed]
  7. Boshoff, C.; Weiss, R. Aids-related malignancies. Nat. Rev. Cancer 2002, 2, 373–382. [Google Scholar] [CrossRef] [PubMed]
  8. Dukers, N.H.; Rezza, G. Human herpesvirus 8 epidemiology: What we do and do not know. Aids 2003, 17, 1717–1730. [Google Scholar] [CrossRef] [PubMed]
  9. Kedes, D.H.; Operskalski, E.; Busch, M.; Kohn, R.; Flood, J.; Ganem, D. The seroepidemiology of human herpesvirus 8 (kaposi’s sarcoma-associated herpesvirus): Distribution of infection in KS risk groups and evidence for sexual transmission. Nat. Med. 1996, 2, 918–924. [Google Scholar] [CrossRef] [PubMed]
  10. Zong, J.C.; Ciufo, D.M.; Alcendor, D.J.; Wan, X.; Nicholas, J.; Browning, P.J.; Rady, P.L.; Tyring, S.K.; Orenstein, J.M.; Rabkin, C.S.; et al. High-level variability in the ORF-K1 membrane protein gene at the left end of the kaposi’s sarcoma-associated herpesvirus genome defines four major virus subtypes and multiple variants or clades in different human populations. J. Virol. 1999, 73, 4156–4170. [Google Scholar] [PubMed]
  11. Poole, L.J.; Zong, J.C.; Ciufo, D.M.; Alcendor, D.J.; Cannon, J.S.; Ambinder, R.; Orenstein, J.M.; Reitz, M.S.; Hayward, G.S. Comparison of genetic variability at multiple loci across the genomes of the major subtypes of kaposi’s sarcoma-associated herpesvirus reveals evidence for recombination and for two distinct types of open reading frame K15 alleles at the right-hand end. J. Virol. 1999, 73, 6646–6660. [Google Scholar] [PubMed]
  12. Hayward, G.S.; Zong, J.C. Modern evolutionary history of the human KSHV genome. Curr. Top. Microbiol. Immunol. 2007, 312, 1–42. [Google Scholar] [PubMed]
  13. Sanger, F.; Coulson, A.R. A rapid method for determining sequences in DNA by primed synthesis with DNA polymerase. J. Mol. Biol. 1975, 94, 441–448. [Google Scholar] [CrossRef]
  14. Maxam, A.M.; Gilbert, W. A new method for sequencing DNA. Proc. Natl. Acad. Sci. USA 1977, 74, 560–564. [Google Scholar] [CrossRef] [PubMed]
  15. Morozova, O.; Marra, M.A. Applications of next-generation sequencing technologies in functional genomics. Genomics 2008, 92, 255–264. [Google Scholar] [CrossRef] [PubMed]
  16. Meyerson, M.; Gabriel, S.; Getz, G. Advances in understanding cancer genomes through second-generation sequencing. Nat. Rev. Genet. 2010, 11, 685–696. [Google Scholar] [CrossRef] [PubMed]
  17. Russo, J.J.; Bohenzky, R.A.; Chien, M.C.; Chen, J.; Yan, M.; Maddalena, D.; Parry, J.P.; Peruzzi, D.; Edelman, I.S.; Chang, Y.; et al. Nucleotide sequence of the kaposi sarcoma-associated herpesvirus (HHV8). Proc. Natl. Acad. Sci. USA 1996, 93, 14862–14867. [Google Scholar] [CrossRef] [PubMed]
  18. Neipel, F.; Albrecht, J.C.; Fleckenstein, B. Cell-homologous genes in the kaposi’s sarcoma-associated rhadinovirus human herpesvirus 8: Determinants of its pathogenicity? J. Virol. 1997, 71, 4187–4192. [Google Scholar] [PubMed]
  19. Rezaee, S.A.; Cunningham, C.; Davison, A.J.; Blackbourn, D.J. Kaposi’s sarcoma-associated herpesvirus immune modulation: An overview. J. Gen. Virol. 2006, 87, 1781–1804. [Google Scholar] [CrossRef] [PubMed]
  20. Zhou, F.C.; Zhang, Y.J.; Deng, J.H.; Wang, X.P.; Pan, H.Y.; Hettler, E.; Gao, S.J. Efficient infection by a recombinant kaposi’s sarcoma-associated herpesvirus cloned in a bacterial artificial chromosome: Application for genetic analysis. J. Virol. 2002, 76, 6185–6196. [Google Scholar] [CrossRef] [PubMed]
  21. Yakushko, Y.; Hackmann, C.; Gunther, T.; Ruckert, J.; Henke, M.; Koste, L.; Alkharsah, K.; Bohne, J.; Grundhoff, A.; Schulz, T.F.; et al. Kaposi’s sarcoma-associated herpesvirus bacterial artificial chromosome contains a duplication of a long unique-region fragment within the terminal repeat region. J. Virol. 2011, 85, 4612–4617. [Google Scholar] [CrossRef] [PubMed]
  22. Brulois, K.F.; Chang, H.; Lee, A.S.; Ensser, A.; Wong, L.Y.; Toth, Z.; Lee, S.H.; Lee, H.R.; Myoung, J.; Ganem, D.; et al. Construction and manipulation of a new kaposi’s sarcoma-associated herpesvirus bacterial artificial chromosome clone. J. Virol. 2012, 86, 9708–9720. [Google Scholar] [CrossRef] [PubMed]
  23. Tamburro, K.M.; Yang, D.; Poisson, J.; Fedoriw, Y.; Roy, D.; Lucas, A.; Sin, S.-H.; Malouf, N.; Moylan, V.; Damania, B.; et al. Vironome of kaposi sarcoma associated herpesvirus-inflammatory cytokine syndrome in an aids patient reveals co-infection of human herpesvirus 8 and human herpesvirus 6a. Virology 2012, 433, 220–225. [Google Scholar] [CrossRef] [PubMed]
  24. Olp, L.N.; Jeanniard, A.; Marimo, C.; West, J.T.; Wood, C. Whole-genome sequencing of kaposi’s sarcoma-associated herpesvirus from zambian kaposi’s sarcoma biopsy specimens reveals unique viral diversity. J. Virol. 2015, 89, 12299–12308. [Google Scholar] [PubMed]
  25. Bruce, A.G.; Ryan, J.T.; Thomas, M.J.; Peng, X.; Grundhoff, A.; Tsai, C.C.; Rose, T.M. Next-generation sequence analysis of the genome of RFHVMn, the macaque homolog of kaposi’s sarcoma (KS)-associated herpesvirus, from a KS-like tumor of a pig-tailed macaque. J. Virol. 2013, 87, 13676–13693. [Google Scholar] [CrossRef] [PubMed]
  26. Depledge, D.P.; Palser, A.L.; Watson, S.J.; Lai, I.Y.; Gray, E.R.; Grant, P.; Kanda, R.K.; Leproust, E.; Kellam, P.; Breuer, J. Specific capture and whole-genome sequencing of viruses from clinical samples. PLoS ONE 2011, 6, e27805. [Google Scholar] [CrossRef] [PubMed]
  27. Byun, M.; Abhyankar, A.; Lelarge, V.; Plancoulaine, S.; Palanduz, A.; Telhan, L.; Boisson, B.; Picard, C.; Dewell, S.; Zhao, C.; et al. Whole-exome sequencing-based discovery of stim1 deficiency in a child with fatal classic kaposi sarcoma. J. Exp. Med. 2010, 207, 2307–2312. [Google Scholar] [CrossRef] [PubMed]
  28. Byun, M.; Ma, C.S.; Akcay, A.; Pedergnana, V.; Palendira, U.; Myoung, J.; Avery, D.T.; Liu, Y.; Abhyankar, A.; Lorenzo, L.; et al. Inherited human OX40 deficiency underlying classic kaposi sarcoma of childhood. J. Exp. Med. 2013, 210, 1743–1759. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Croft, M. Control of immunity by the TNFR-related molecule OX40 (CD134). Annu. Rev. Immunol. 2010, 28, 57–78. [Google Scholar] [CrossRef] [PubMed]
  30. Aavikko, M.; Kaasinen, E.; Nieminen, J.K.; Byun, M.; Donner, I.; Mancuso, R.; Ferrante, P.; Clerici, M.; Brambilla, L.; Tourlaki, A.; et al. Whole-genome sequencing identifies STAT4 as a putative susceptibility gene in classic kaposi sarcoma. J. Infect. Dis. 2015, 211, 1842–1851. [Google Scholar] [CrossRef] [PubMed]
  31. Yang, D.; Chen, W.; Xiong, J.; Sherrod, C.J.; Henry, D.H.; Dittmer, D.P. Interleukin 1 receptor-associated kinase 1 (IRAK1) mutation is a common, essential driver for kaposi sarcoma herpesvirus lymphoma. Proc. Natl. Acad. Sci. USA 2014, 111, E4762–E4768. [Google Scholar] [CrossRef] [PubMed]
  32. Wang, Z.; Gerstein, M.; Snyder, M. RNA-seq: A revolutionary tool for transcriptomics. Nat. Rev. Genet. 2009, 10, 57–63. [Google Scholar] [CrossRef] [PubMed]
  33. Martin, J.A.; Wang, Z. Next-generation transcriptome assembly. Nat. Rev. Genet. 2011, 12, 671–682. [Google Scholar] [CrossRef] [PubMed]
  34. Cannon, J.S.; Ciufo, D.; Hawkins, A.L.; Griffin, C.A.; Borowitz, M.J.; Hayward, G.S.; Ambinder, R.F. A new primary effusion lymphoma-derived cell line yields a highly infectious kaposi’s sarcoma herpesvirus-containing supernatant. J. Virol. 2000, 74, 10187–10193. [Google Scholar] [CrossRef] [PubMed]
  35. Grundhoff, A.; Ganem, D. Inefficient establishment of KSHV latency suggests an additional role for continued lytic replication in kaposi sarcoma pathogenesis. J. Clin. Investig. 2004, 113, 124–136. [Google Scholar] [CrossRef] [PubMed]
  36. Renne, R.; Zhong, W.; Herndier, B.; McGrath, M.; Abbey, N.; Kedes, D.; Ganem, D. Lytic growth of kaposi’s sarcoma-associated herpesvirus (human herpesvirus 8) in culture. Nat. Med. 1996, 2, 342–346. [Google Scholar] [CrossRef] [PubMed]
  37. Wang, C.Y.; Sugden, B. New viruses shake old paradigms. J. Clin. Investig. 2004, 113, 21–23. [Google Scholar] [CrossRef] [PubMed]
  38. Dittmer, D.; Lagunoff, M.; Renne, R.; Staskus, K.; Haase, A.; Ganem, D. A cluster of latently expressed genes in kaposi’s sarcoma-associated herpesvirus. J. Virol. 1998, 72, 8309–8315. [Google Scholar] [PubMed]
  39. Fakhari, F.D.; Dittmer, D.P. Charting latency transcripts in kaposi’s sarcoma-associated herpesvirus by whole-genome real-time quantitative PCR. J. Virol. 2002, 76, 6213–6223. [Google Scholar] [CrossRef] [PubMed]
  40. Samols, M.A.; Hu, J.; Skalsky, R.L.; Renne, R. Cloning and identification of a microrna cluster within the latency-associated region of kaposi’s sarcoma-associated herpesvirus. J. Virol. 2005, 79, 9301–9305. [Google Scholar] [CrossRef] [PubMed]
  41. Sarid, R.; Flore, O.; Bohenzky, R.A.; Chang, Y.; Moore, P.S. Transcription mapping of the kaposi’s sarcoma-associated herpesvirus (human herpesvirus 8) genome in a body cavity-based lymphoma cell line (BC-1). J. Virol. 1998, 72, 1005–1012. [Google Scholar] [PubMed]
  42. Jeong, J.; Papin, J.; Dittmer, D. Differential regulation of the overlapping kaposi’s sarcoma-associated herpesvirus vgcr (ORF74) and lana (ORF73) promoters. J. Virol. 2001, 75, 1798–1807. [Google Scholar] [CrossRef] [PubMed]
  43. Lu, C.; Zeng, Y.; Huang, Z.; Huang, L.; Qian, C.; Tang, G.; Qin, D. Human herpesvirus 6 activates lytic cycle replication of kaposi’s sarcoma-associated herpesvirus. Am. J. Pathol. 2005, 166, 173–183. [Google Scholar] [CrossRef]
  44. Mercader, M.; Taddeo, B.; Panella, J.R.; Chandran, B.; Nickoloff, B.J.; Foreman, K.E. Induction of HHV-8 lytic cycle replication by inflammatory cytokines produced by HIV-1-infected T cells. Am. J. Pathol. 2000, 156, 1961–1971. [Google Scholar] [CrossRef]
  45. Vieira, J.; O’Hearn, P.; Kimball, L.; Chandran, B.; Corey, L. Activation of kaposi’s sarcoma-associated herpesvirus (human herpesvirus 8) lytic replication by human cytomegalovirus. J. Virol. 2001, 75, 1378–1386. [Google Scholar] [CrossRef] [PubMed]
  46. Ye, F.; Zhou, F.; Bedolla, R.G.; Jones, T.; Lei, X.; Kang, T.; Guadalupe, M.; Gao, S.J. Reactive oxygen species hydrogen peroxide mediates kaposi’s sarcoma-associated herpesvirus reactivation from latency. PLoS Pathog. 2011, 7, e1002054. [Google Scholar] [CrossRef] [PubMed]
  47. Davis, D.A.; Rinderknecht, A.S.; Zoeteweij, J.P.; Aoki, Y.; Read-Connole, E.L.; Tosato, G.; Blauvelt, A.; Yarchoan, R. Hypoxia induces lytic replication of kaposi sarcoma-associated herpesvirus. Blood 2001, 97, 3244–3250. [Google Scholar] [CrossRef] [PubMed]
  48. Honess, R.W.; Roizman, B. Regulation of herpesvirus macromolecular synthesis. I. Cascade regulation of the synthesis of three groups of viral proteins. J. Virol. 1974, 14, 8–19. [Google Scholar] [PubMed]
  49. Sun, R.; Lin, S.F.; Staskus, K.; Gradoville, L.; Grogan, E.; Haase, A.; Miller, G. Kinetics of kaposi’s sarcoma-associated herpesvirus gene expression. J. Virol. 1999, 73, 2232–2242. [Google Scholar] [PubMed]
  50. Zhu, F.X.; Cusano, T.; Yuan, Y. Identification of the immediate-early transcripts of kaposi’s sarcoma-associated herpesvirus. J. Virol. 1999, 73, 5556–5567. [Google Scholar] [PubMed]
  51. Jenner, R.G.; Alba, M.M.; Boshoff, C.; Kellam, P. Kaposi’s sarcoma-associated herpesvirus latent and lytic gene expression as revealed by DNA arrays. J. Virol. 2001, 75, 891–902. [Google Scholar] [CrossRef] [PubMed]
  52. Lukac, D.M.; Kirshner, J.R.; Ganem, D. Transcriptional activation by the product of open reading frame 50 of kaposi’s sarcoma-associated herpesvirus is required for lytic viral reactivation in B cells. J. Virol. 1999, 73, 9348–9361. [Google Scholar] [PubMed]
  53. Lukac, D.M.; Renne, R.; Kirshner, J.R.; Ganem, D. Reactivation of kaposi’s sarcoma-associated herpesvirus infection from latency by expression of the ORF 50 transactivator, a homolog of the EBV R protein. Virology 1998, 252, 304–312. [Google Scholar] [CrossRef] [PubMed]
  54. Sun, R.; Lin, S.F.; Gradoville, L.; Yuan, Y.; Zhu, F.; Miller, G. A viral gene that activates lytic cycle expression of kaposi’s sarcoma-associated herpesvirus. Proc. Natl. Acad. Sci. USA 1998, 95, 10866–10871. [Google Scholar] [CrossRef] [PubMed]
  55. Chandriani, S.; Xu, Y.; Ganem, D. The lytic transcriptome of kaposi’s sarcoma-associated herpesvirus reveals extensive transcription of noncoding regions, including regions antisense to important genes. J. Virol. 2010, 84, 7934–7942. [Google Scholar] [CrossRef] [PubMed]
  56. Dittmer, D.P. Transcription profile of kaposi’s sarcoma-associated herpesvirus in primary kaposi’s sarcoma lesions as determined by real-time PCR arrays. Cancer Res. 2003, 63, 2010–2015. [Google Scholar] [PubMed]
  57. Dresang, L.R.; Teuton, J.R.; Feng, H.; Jacobs, J.M.; Camp, D.G., 2nd; Purvine, S.O.; Gritsenko, M.A.; Li, Z.; Smith, R.D.; Sugden, B.; et al. Coupled transcriptome and proteome analysis of human lymphotropic tumor viruses: Insights on the detection and discovery of viral genes. BMC Genom. 2011, 12. [Google Scholar] [CrossRef]
  58. Arias, C.; Weisburd, B.; Stern-Ginossar, N.; Mercier, A.; Madrid, A.S.; Bellare, P.; Holdorf, M.; Weissman, J.S.; Ganem, D. KSHV 2.0: A comprehensive annotation of the kaposi’s sarcoma-associated herpesvirus genome using next-generation sequencing reveals novel genomic and functional features. PLoS Pathog. 2014, 10, e1003847. [Google Scholar] [CrossRef] [PubMed]
  59. Lin, X.; Li, X.; Liang, D.; Lan, K. Micrornas and unusual small rnas discovered in kaposi’s sarcoma-associated herpesvirus virions. J. Virol. 2012, 86, 12717–12730. [Google Scholar] [CrossRef] [PubMed]
  60. Bechtel, J.; Grundhoff, A.; Ganem, D. RNAs in the virion of kaposi’s sarcoma-associated herpesvirus. J. Virol. 2005, 79, 10138–10146. [Google Scholar] [CrossRef] [PubMed]
  61. Purushothaman, P.; Thakker, S.; Verma, S.C. Transcriptome analysis of kaposi’s sarcoma-associated herpesvirus during de novo primary infection of human B and endothelial cells. J. Virol. 2015, 89, 3093–3111. [Google Scholar] [CrossRef] [PubMed]
  62. Lu, F.; Day, L.; Gao, S.J.; Lieberman, P.M. Acetylation of the latency-associated nuclear antigen regulates repression of kaposi’s sarcoma-associated herpesvirus lytic transcription. J. Virol. 2006, 80, 5273–5282. [Google Scholar] [CrossRef] [PubMed]
  63. Sun, R.; Liang, D.; Gao, Y.; Lan, K. Kaposi’s sarcoma-associated herpesvirus-encoded LANA interacts with host KAP1 to facilitate establishment of viral latency. J. Virol. 2014, 88, 7331–7344. [Google Scholar] [CrossRef] [PubMed]
  64. Cai, Q.; Cai, S.; Zhu, C.; Verma, S.C.; Choi, J.Y.; Robertson, E.S. A unique sumo-2-interacting motif within lana is essential for KSHV latency. PLoS Pathog. 2013, 9, e1003750. [Google Scholar] [CrossRef] [PubMed]
  65. Gjyshi, O.; Roy, A.; Dutta, S.; Veettil, M.V.; Dutta, D.; Chandran, B. Activated NRF2 interacts with kaposi’s sarcoma-associated herpesvirus latency protein LANA-1 and host protein KAP1 to mediate global lytic gene repression. J. Virol. 2015, 89, 7874–7892. [Google Scholar] [CrossRef] [PubMed]
  66. Verma, D.; Li, D.J.; Krueger, B.; Renne, R.; Swaminathan, S. Identification of the physiological gene targets of the essential lytic replicative kaposi’s sarcoma-associated herpesvirus ORF57 protein. J. Virol. 2015, 89, 1688–1702. [Google Scholar] [CrossRef] [PubMed]
  67. Gong, D.; Wu, N.C.; Xie, Y.; Feng, J.; Tong, L.; Brulois, K.F.; Luan, H.; Du, Y.; Jung, J.U.; Wang, C.-Y.; et al. Kaposi’s sarcoma-associated herpesvirus ORF18 and ORF30 are essential for late gene expression during lytic replication. J. Virol. 2014, 88, 11369–11382. [Google Scholar] [CrossRef] [PubMed]
  68. Davis, Z.H.; Hesser, C.R.; Park, J.; Glaunsinger, B.A. Interaction between ORF24 and ORF34 in the kaposi’s sarcoma-associated herpesvirus late gene transcription factor complex is essential for viral late gene expression. J. Virol. 2015, 90, 599–604. [Google Scholar] [CrossRef] [PubMed]
  69. Aubry, V.; Mure, F.; Mariame, B.; Deschamps, T.; Wyrwicz, L.S.; Manet, E.; Gruffat, H. Epstein-barr virus late gene transcription depends on the assembly of a virus-specific preinitiation complex. J. Virol. 2014, 88, 12825–12838. [Google Scholar] [CrossRef] [PubMed]
  70. Davis, Z.H.; Verschueren, E.; Jang, G.M.; Kleffman, K.; Johnson, J.R.; Park, J.; von Dollen, J.; Maher, M.C.; Johnson, T.; Newton, W.; et al. Global mapping of herpesvirus-host protein complexes reveals a transcription strategy for late genes. Mol. Cell 2015, 57, 349–360. [Google Scholar] [CrossRef] [PubMed]
  71. Glaunsinger, B.; Ganem, D. Lytic KSHV infection inhibits host gene expression by accelerating global mrna turnover. Mol. Cell 2004, 13, 713–723. [Google Scholar] [CrossRef]
  72. Glaunsinger, B.; Chavez, L.; Ganem, D. The exonuclease and host shutoff functions of the SOX protein of kaposi’s sarcoma-associated herpesvirus are genetically separable. J. Virol. 2005, 79, 7396–7401. [Google Scholar] [CrossRef] [PubMed]
  73. Covarrubias, S.; Gaglia, M.M.; Kumar, G.R.; Wong, W.; Jackson, A.O.; Glaunsinger, B.A. Coordinated destruction of cellular messages in translation complexes by the gammaherpesvirus host shutoff factor and the mammalian exonuclease XRN1. PLoS Pathog. 2011, 7, e1002339. [Google Scholar] [CrossRef] [PubMed]
  74. Abernathy, E.; Gilbertson, S.; Alla, R.; Glaunsinger, B. Viral nucleases induce an mRNA degradation-transcription feedback loop in mammalian cells. Cell Host Microbe 2015, 18, 243–253. [Google Scholar] [CrossRef] [PubMed]
  75. Nishi, J.; Arimura, K.; Utsunomiya, A.; Yonezawa, S.; Kawakami, K.; Maeno, N.; Ijichi, O.; Ikarimoto, N.; Nakata, M.; Kitajima, I.; et al. Expression of vascular endothelial growth factor in sera and lymph nodes of the plasma cell type of castleman’s disease. Br. J. Haematol. 1999, 104, 482–485. [Google Scholar] [CrossRef] [PubMed]
  76. Glaunsinger, B.; Ganem, D. Highly selective escape from KSHV-mediated host mrna shutoff and its implications for viral pathogenesis. J. Exp. Med. 2004, 200, 391–398. [Google Scholar] [CrossRef] [PubMed]
  77. Clyde, K.; Glaunsinger, B.A. Deep sequencing reveals direct targets of gammaherpesvirus-induced mRNA decay and suggests that multiple mechanisms govern cellular transcript escape. PLoS ONE 2011, 6, e19655. [Google Scholar] [CrossRef] [PubMed]
  78. Gaglia, M.M.; Rycroft, C.H.; Glaunsinger, B.A. Transcriptome-wide cleavage site mapping on cellular mrnas reveals features underlying sequence-specific cleavage by the viral ribonuclease SOX. PLoS Pathog. 2015, 11, e1005305. [Google Scholar] [CrossRef] [PubMed]
  79. Sei, E.; Wang, T.; Hunter, O.V.; Xie, Y.; Conrad, N.K. HITS-CLIP analysis uncovers a link between the kaposi’s sarcoma-associated herpesvirus ORF57 protein and host pre-mRNA metabolism. PLoS Pathog. 2015, 11, e1004652. [Google Scholar] [CrossRef] [PubMed]
  80. Cao, S.; Strong, M.J.; Wang, X.; Moss, W.N.; Concha, M.; Lin, Z.; O’Grady, T.; Baddoo, M.; Fewell, C.; Renne, R.; et al. High-throughput RNA sequencing-based virome analysis of 50 lymphoma cell lines from the cancer cell line encyclopedia project. J. Virol. 2015, 89, 713–729. [Google Scholar] [CrossRef] [PubMed]
  81. Carthew, R.W.; Sontheimer, E.J. Origins and mechanisms of miRNAs and siRNAs. Cell 2009, 136, 642–655. [Google Scholar] [CrossRef] [PubMed]
  82. Bartel, D.P. Micrornas: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef]
  83. Lei, X.; Bai, Z.; Ye, F.; Huang, Y.; Gao, S.J. Regulation of herpesvirus lifecycle by viral micrornas. Virulence 2010, 1, 433–435. [Google Scholar] [CrossRef] [PubMed]
  84. Qin, Z.; Jakymiw, A.; Findlay, V.; Parsons, C. KSHV-encoded microRNAs: Lessons for viral cancer pathogenesis and emerging concepts. Int. J. Cell Biol. 2012. [Google Scholar] [CrossRef] [PubMed]
  85. Li, Y.; Kowdley, K.V. MicroRNAs in common human diseases. Genom. Proteom. Bioinf. 2012, 10, 246–253. [Google Scholar] [CrossRef] [PubMed]
  86. Gottwein, E. Kaposi’s sarcoma-associated herpesvirus microRNAs. Front. Microbiol. 2012, 3. [Google Scholar] [CrossRef] [PubMed]
  87. Ramalingam, D.; Kieffer-Kwon, P.; Ziegelbauer, J.M. Emerging themes from EBV and KSHV microRNA targets. Viruses 2012, 4, 1687–1710. [Google Scholar] [CrossRef] [PubMed]
  88. Cai, X.; Cullen, B.R. Transcriptional origin of kaposi’s sarcoma-associated herpesvirus microRNAs. J. Virol. 2006, 80, 2234–2242. [Google Scholar] [CrossRef] [PubMed]
  89. Pfeffer, S.; Sewer, A.; Lagos-Quintana, M.; Sheridan, R.; Sander, C.; Grasser, F.A.; van Dyk, L.F.; Ho, C.K.; Shuman, S.; Chien, M.; et al. Identification of microRNAs of the herpesvirus family. Nat. Methods. 2005, 2, 269–276. [Google Scholar] [CrossRef] [PubMed]
  90. Cai, X.; Lu, S.; Zhang, Z.; Gonzalez, C.M.; Damania, B.; Cullen, B.R. Kaposi’s sarcoma-associated herpesvirus expresses an array of viral microRNAs in latently infected cells. Proc. Natl. Acad. Sci. USA 2005, 102, 5570–5575. [Google Scholar] [CrossRef] [PubMed]
  91. Marshall, V.; Parks, T.; Bagni, R.; Wang, C.D.; Samols, M.A.; Hu, J.; Wyvil, K.M.; Aleman, K.; Little, R.F.; Yarchoan, R.; et al. Conservation of virally encoded microRNAs in kaposi sarcoma-associated herpesvirus in primary effusion lymphoma cell lines and in patients with kaposi sarcoma or multicentric castleman disease. J. Infect. Dis. 2007, 195, 645–659. [Google Scholar] [CrossRef] [PubMed]
  92. Grundhoff, A.; Sullivan, C.S.; Ganem, D. A combined computational and microarray-based approach identifies novel microRNAs encoded by human gamma-herpesviruses. RNA 2006, 12, 733–750. [Google Scholar] [CrossRef] [PubMed]
  93. Umbach, J.L.; Cullen, B.R. In-depth analysis of kaposi’s sarcoma-associated herpesvirus microRNA expression provides insights into the mammalian microRNA-processing machinery. J. Virol. 2010, 84, 695–703. [Google Scholar] [CrossRef] [PubMed]
  94. AuCoin, D.P.; Colletti, K.S.; Xu, Y.; Cei, S.A.; Pari, G.S. Kaposi’s sarcoma-associated herpesvirus (human herpesvirus 8) contains two functional lytic origins of DNA replication. J. Virol. 2002, 76, 7890–7896. [Google Scholar] [CrossRef] [PubMed]
  95. Lin, C.L.; Li, H.; Wang, Y.; Zhu, F.X.; Kudchodkar, S.; Yuan, Y. Kaposi’s sarcoma-associated herpesvirus lytic origin (ori-Lyt)-dependent DNA replication: Identification of the ori-Lyt and association of K8 bZip protein with the origin. J. Virol. 2003, 77, 5578–5588. [Google Scholar] [CrossRef] [PubMed]
  96. Lin, Y.T.; Kincaid, R.P.; Arasappan, D.; Dowd, S.E.; Hunicke-Smith, S.P.; Sullivan, C.S. Small RNA profiling reveals antisense transcription throughout the KSHV genome and novel small RNAs. RNA 2010, 16, 1540–1558. [Google Scholar] [CrossRef] [PubMed]
  97. Qin, Z.; Peruzzi, F.; Reiss, K.; Dai, L. Role of host microRNAs in kaposi’s sarcoma-associated herpesvirus pathogenesis. Viruses 2014, 6, 4571–4580. [Google Scholar] [CrossRef] [PubMed]
  98. Boss, I.W.; Plaisance, K.B.; Renne, R. Role of virus-encoded microRNAs in herpesvirus biology. Trends Microbiol. 2009, 17, 544–553. [Google Scholar] [CrossRef] [PubMed]
  99. Lee, H.R.; Lee, S.; Chaudhary, P.M.; Gill, P.; Jung, J.U. Immune evasion by kaposi’s sarcoma-associated herpesvirus. Future Microbiol. 2010, 5, 1349–1365. [Google Scholar] [CrossRef] [PubMed]
  100. Moody, R.; Zhu, Y.; Huang, Y.; Cui, X.; Jones, T.; Bedolla, R.; Lei, X.; Bai, Z.; Gao, S.J. KSHV microRNAs mediate cellular transformation and tumorigenesis by redundantly targeting cell growth and survival pathways. PLoS Pathog. 2013, 9, e1003857. [Google Scholar] [CrossRef] [PubMed]
  101. Haecker, I.; Renne, R. HITS-CLIP and PAR-CLIP advance viral miRNA targetome analysis. Crit. Rev. Eukaryot. Gene Expr. 2014, 24, 101–116. [Google Scholar] [CrossRef] [PubMed]
  102. Haecker, I.; Gay, L.A.; Yang, Y.; Hu, J.; Morse, A.M.; McIntyre, L.M.; Renne, R. Ago HITS-CLIP expands understanding of kaposi’s sarcoma-associated herpesvirus miRNA function in primary effusion lymphomas. PLoS Pathog. 2012, 8, e1002884. [Google Scholar] [CrossRef] [PubMed]
  103. Gottwein, E.; Corcoran, D.L.; Mukherjee, N.; Skalsky, R.L.; Hafner, M.; Nusbaum, J.D.; Shamulailatpam, P.; Love, C.L.; Dave, S.S.; Tuschl, T.; et al. Viral microRNA targetome of KSHV-infected primary effusion lymphoma cell lines. Cell Host Microbe 2011, 10, 515–526. [Google Scholar] [CrossRef] [PubMed]
  104. Viollet, C.; Davis, D.A.; Reczko, M.; Ziegelbauer, J.M.; Pezzella, F.; Ragoussis, J.; Yarchoan, R. Next-generation sequencing analysis reveals differential expression profiles of miRNA-mRNA target pairs in KSHV-infected cells. PLoS ONE 2015, 10, e0126439. [Google Scholar]
  105. Shaw, G.; Kamen, R. A conserved AU sequence from the 3′ untranslated region of GM-CSF mRNA mediates selective mRNA degradation. Cell 1986, 46, 659–667. [Google Scholar] [CrossRef]
  106. Elkon, R.; Ugalde, A.P.; Agami, R. Alternative cleavage and polyadenylation: Extent, regulation and function. Nat. Rev. Genet. 2013, 14, 496–506. [Google Scholar] [CrossRef] [PubMed]
  107. McClure, L.V.; Kincaid, R.P.; Burke, J.M.; Grundhoff, A.; Sullivan, C.S. Comprehensive mapping and analysis of kaposi’s sarcoma-associated herpesvirus 3′ UTRs identify differential posttranscriptional control of gene expression in lytic vs. latent infection. J. Virol. 2013, 87, 12838–12849. [Google Scholar] [CrossRef] [PubMed]
  108. Majerciak, V.; Ni, T.; Yang, W.; Meng, B.; Zhu, J.; Zheng, Z.M. A viral genome landscape of RNA polyadenylation from KSHV latent to lytic infection. PLoS Pathog. 2013, 9, e1003749. [Google Scholar] [CrossRef] [PubMed]
  109. Bai, Z.; Huang, Y.; Li, W.; Zhu, Y.; Jung, J.U.; Lu, C.; Gao, S.J. Genomewide mapping and screening of kaposi’s sarcoma-associated herpesvirus (KSHV) 3′ untranslated regions identify bicistronic and polycistronic viral transcripts as frequent targets of KSHV microRNAs. J. Virol. 2014, 88, 377–392. [Google Scholar] [CrossRef] [PubMed]
  110. Zhu, Y.; Huang, Y.; Jung, J.U.; Lu, C.; Gao, S.J. Viral miRNA targeting of bicistronic and polycistronic transcripts. Curr. Opin. Virol. 2014, 7, 66–72. [Google Scholar] [CrossRef] [PubMed]
  111. Haque, M.; Wang, V.; Davis, D.A.; Zheng, Z.M.; Yarchoan, R. Genetic organization and hypoxic activation of the kaposi’s sarcoma-associated herpesvirus ORF34-37 gene cluster. J. Virol. 2006, 80, 7037–7051. [Google Scholar] [CrossRef] [PubMed]
  112. Bai, Z.; Zhou, F.; Lei, X.; Ma, X.; Lu, C.; Gao, S.J. A cluster of transcripts encoded by KSHV ORF30-33 gene locus. Virus Genes 2012, 44, 225–236. [Google Scholar] [CrossRef] [PubMed]
  113. Bellare, P.; Ganem, D. Regulation of KSHV lytic switch protein expression by a virus-encoded microRNA: An evolutionary adaptation that fine-tunes lytic reactivation. Cell Host Microbe 2009, 6, 570–575. [Google Scholar] [CrossRef] [PubMed]
  114. Lin, X.; Liang, D.; He, Z.; Deng, Q.; Robertson, E.S.; Lan, K. miR-K12-7-5p encoded by kaposi’s sarcoma-associated herpesvirus stabilizes the latent state by targeting viral ORF50/RTA. PLoS ONE 2011, 6, e16224. [Google Scholar] [CrossRef] [PubMed]
  115. Lin, H.R.; Ganem, D. Viral microRNA target allows insight into the role of translation in governing microrna target accessibility. Proc. Natl. Acad. Sci. USA 2011, 108, 5148–5153. [Google Scholar] [CrossRef] [PubMed]
  116. Rossetto, C.C.; Pari, G.S. PAN’S labyrinth: Molecular biology of kaposi’s sarcoma-associated herpesvirus (KSHV) PAN RNA, a multifunctional long noncoding RNA. Viruses 2014, 6, 4212–4226. [Google Scholar] [CrossRef] [PubMed]
  117. Campbell, M.; Kung, H.J.; Izumiya, Y. Long non-coding RNA and epigenetic gene regulation of KSHV. Viruses 2014, 6, 4165–4177. [Google Scholar] [CrossRef] [PubMed]
  118. Sun, R.; Lin, S.-F.; Gradoville, L.; Miller, G. Polyadenylylated nuclear RNA encoded by kaposi sarcoma-associated herpesvirus. Proc. Natl. Acad. Sci. USA 1996, 93, 11883–11888. [Google Scholar] [CrossRef] [PubMed]
  119. Rossetto, C.C.; Tarrant-Elorza, M.; Verma, S.; Purushothaman, P.; Pari, G.S. Regulation of viral and cellular gene expression by kaposi’s sarcoma-associated herpesvirus polyadenylated nuclear RNA. J. Virol. 2013, 87, 5540–5553. [Google Scholar] [CrossRef] [PubMed]
  120. Rossetto, C.C.; Pari, G. KSHV PAN RNA associates with demethylases UTX and JMJD3 to activate lytic replication through a physical interaction with the virus genome. PLoS Pathog. 2012, 8, e1002680. [Google Scholar] [CrossRef] [PubMed]
  121. Borah, S.; Darricarrere, N.; Darnell, A.; Myoung, J.; Steitz, J.A. A viral nuclear noncoding RNA binds re-localized poly(a) binding protein and is required for late KSHV gene expression. PLoS Pathog. 2011, 7, e1002300. [Google Scholar] [CrossRef] [PubMed]
  122. Campbell, M.; Kim, K.Y.; Chang, P.C.; Huerta, S.; Shevchenko, B.; Wang, D.H.; Izumiya, C.; Kung, H.J.; Izumiya, Y. A lytic viral long noncoding RNA modulates the function of a latent protein. J. Virol. 2014, 88, 1843–1848. [Google Scholar] [CrossRef] [PubMed]
  123. Song, M.J.; Brown, H.J.; Wu, T.T.; Sun, R. Transcription activation of polyadenylated nuclear RNA by Rta in human herpesvirus 8/kaposi’s sarcoma-associated herpesvirus. J. Virol. 2001, 75, 3129–3140. [Google Scholar] [CrossRef] [PubMed]
  124. Chang, P.J.; Shedd, D.; Gradoville, L.; Cho, M.S.; Chen, L.W.; Chang, J.; Miller, G. Open reading frame 50 protein of kaposi’s sarcoma-associated herpesvirus directly activates the viral PAN and K12 genes by binding to related response elements. J. Virol. 2002, 76, 3168–3178. [Google Scholar] [CrossRef] [PubMed]
  125. Massimelli, M.J.; Kang, J.G.; Majerciak, V.; Le, S.Y.; Liewehr, D.J.; Steinberg, S.M.; Zheng, Z.M. Stability of a long noncoding viral RNA depends on a 9-nt core element at the RNA 5′ end to interact with viral ORF57 and cellular PABPC1. Int. J. Biol. Sci. 2011, 7, 1145–1160. [Google Scholar] [CrossRef] [PubMed]
  126. Sahin, B.B.; Patel, D.; Conrad, N.K. Kaposi’s sarcoma-associated herpesvirus ORF57 protein binds and protects a nuclear noncoding RNA from cellular RNA decay pathways. PLoS Pathog. 2010, 6, e1000799. [Google Scholar] [CrossRef] [PubMed]
  127. Chu, C.; Qu, K.; Zhong, F.L.; Artandi, S.E.; Chang, H.Y. Genomic maps of long noncoding RNA occupancy reveal principles of RNA-chromatin interactions. Mol. Cell 2011, 44, 667–678. [Google Scholar] [CrossRef] [PubMed]
  128. Jaber, T.; Yuan, Y. A virally encoded small peptide regulates RTA stability and facilitates kaposi’s sarcoma-associated herpesvirus lytic replication. J. Virol. 2013, 87, 3461–3470. [Google Scholar] [CrossRef] [PubMed]
  129. Goldberg, A.D.; Allis, C.D.; Bernstein, E. Epigenetics: A landscape takes shape. Cell 2007, 128, 635–638. [Google Scholar] [CrossRef] [PubMed]
  130. Berger, S.L. The complex language of chromatin regulation during transcription. Nature 2007, 447, 407–412. [Google Scholar] [CrossRef] [PubMed]
  131. Feinberg, A.P.; Tycko, B. The history of cancer epigenetics. Nat. Rev. Cancer 2004, 4, 143–153. [Google Scholar] [CrossRef] [PubMed]
  132. Mattick, J.S.; Makunin, I.V. Non-coding RNA. Hum. Mol. Genet. 2006, 15, R17–R29. [Google Scholar] [CrossRef] [PubMed]
  133. Li, B.; Carey, M.; Workman, J.L. The role of chromatin during transcription. Cell 2007, 128, 707–719. [Google Scholar] [CrossRef] [PubMed]
  134. Park, P.J. Chip-seq: Advantages and challenges of a maturing technology. Nat. Rev. Genet. 2009, 10, 669–680. [Google Scholar] [CrossRef] [PubMed]
  135. Gunther, T.; Grundhoff, A. The epigenetic landscape of latent kaposi sarcoma-associated herpesvirus genomes. PLoS Pathog. 2010, 6, e1000935. [Google Scholar] [CrossRef] [PubMed]
  136. Mercier, A.; Arias, C.; Madrid, A.S.; Holdorf, M.M.; Ganem, D. Site-specific association with host and viral chromatin by kaposi’s sarcoma-associated herpesvirus LANA and its reversal during lytic reactivation. J. Virol. 2014, 88, 6762–6777. [Google Scholar] [CrossRef] [PubMed]
  137. Lu, F.; Tsai, K.; Chen, H.S.; Wikramasinghe, P.; Davuluri, R.V.; Showe, L.; Domsic, J.; Marmorstein, R.; Lieberman, P.M. Identification of host-chromosome binding sites and candidate gene targets for kaposi’s sarcoma-associated herpesvirus LANA. J. Virol. 2012, 86, 5752–5762. [Google Scholar] [CrossRef] [PubMed]
  138. Hilton, I.B.; Simon, J.M.; Lieb, J.D.; Davis, I.J.; Damania, B.; Dittmer, D.P. The open chromatin landscape of kaposi’s sarcoma-associated herpesvirus. J. Virol. 2013, 87, 11831–11842. [Google Scholar] [CrossRef] [PubMed]
  139. Tucker, K.L. Methylated cytosine and the brain: A new base for neuroscience. Neuron 2001, 30, 649–652. [Google Scholar] [CrossRef]
  140. Tahiliani, M.; Koh, K.P.; Shen, Y.; Pastor, W.A.; Bandukwala, H.; Brudno, Y.; Agarwal, S.; Iyer, L.M.; Liu, D.R.; Aravind, L.; et al. Conversion of 5-methylcytosine to 5-hydroxymethylcytosine in mammalian DNA by MLL partner TET1. Science 2009, 324, 930–935. [Google Scholar] [CrossRef] [PubMed]
  141. Newell-Price, J.; Clark, A.J.; King, P. DNA methylation and silencing of gene expression. Trends Endocrinol. Metab. 2000, 11, 142–148. [Google Scholar] [CrossRef]
  142. Shull, A.; Noonepalle, S.; Lee, E.-J.; Choi, J.-H.; Shi, H. Sequencing the cancer methylome. In Cancer Epigenetics; Verma, M., Ed.; Springer: New York, NY, USA, 2015; pp. 627–651. [Google Scholar]
  143. Gavrilov, A.; Eivazova, E.; Priozhkova, I.; Lipinski, M.; Razin, S.; Vassetzky, Y. Chromosome conformation capture (from 3C to 5C) and its chip-based modification. Methods Mol. Biol. 2009, 567, 171–188. [Google Scholar] [PubMed]
  144. Van Berkum, N.L.; Dekker, J. Determining spatial chromatin organization of large genomic regions using 5C technology. Methods Mol. Biol. 2009, 567, 189–213. [Google Scholar] [PubMed]
  145. Giresi, P.G.; Kim, J.; McDaniell, R.M.; Iyer, V.R.; Lieb, J.D. Faire (formaldehyde-assisted isolation of regulatory elements) isolates active regulatory elements from human chromatin. Genome Res. 2007, 17, 877–885. [Google Scholar] [CrossRef] [PubMed]
  146. Meyer, C.A.; Liu, X.S. Identifying and mitigating bias in next-generation sequencing methods for chromatin biology. Nat. Rev. Genet. 2014, 15, 709–721. [Google Scholar] [CrossRef] [PubMed]
  147. Toth, Z.; Brulois, K.; Jung, J.U. The chromatin landscape of kaposi’s sarcoma-associated herpesvirus. Viruses 2013, 5, 1346–1373. [Google Scholar] [CrossRef] [PubMed]
  148. Uppal, T.; Jha, H.C.; Verma, S.C.; Robertson, E.S. Chromatinization of the KSHV genome during the KSHV life cycle. Cancers 2015, 7, 112–142. [Google Scholar] [CrossRef] [PubMed]
  149. Toth, Z.; Brulois, K.; Lee, H.R.; Izumiya, Y.; Tepper, C.; Kung, H.J.; Jung, J.U. Biphasic euchromatin-to-heterochromatin transition on the KSHV genome following de novo infection. PLoS Pathog. 2013, 9, e1003813. [Google Scholar] [CrossRef] [PubMed]
  150. Jha, H.C.; Lu, J.; Verma, S.C.; Banerjee, S.; Mehta, D.; Robertson, E.S. Kaposi’s sarcoma-associated herpesvirus genome programming during the early stages of primary infection of peripheral blood mononuclear cells. MBio 2014. [Google Scholar] [CrossRef] [PubMed]
  151. Chen, J.; Ueda, K.; Sakakibara, S.; Okuno, T.; Parravicini, C.; Corbellino, M.; Yamanishi, K. Activation of latent kaposi’s sarcoma-associated herpesvirus by demethylation of the promoter of the lytic transactivator. Proc. Natl. Acad. Sci. USA 2001, 98, 4119–4124. [Google Scholar] [CrossRef] [PubMed]
  152. Uppal, T.; Banerjee, S.; Sun, Z.; Verma, S.C.; Robertson, E.S. KSHV LANA—The master regulator of KSHV latency. Viruses 2014, 6, 4961–4998. [Google Scholar] [CrossRef] [PubMed]
  153. Ballestas, M.E.; Kaye, K.M. The latency-associated nuclear antigen, a multifunctional protein central to kaposi’s sarcoma-associated herpesvirus latency. Future Microbiol. 2011, 6, 1399–1413. [Google Scholar] [CrossRef] [PubMed]
  154. Verma, S.C.; Lan, K.; Robertson, E. Structure and function of latency-associated nuclear antigen. Curr. Top. Microbiol. Immunol. 2007, 312, 101–136. [Google Scholar] [PubMed]
  155. Garber, A.C.; Shu, M.A.; Hu, J.; Renne, R. DNA binding and modulation of gene expression by the latency-associated nuclear antigen of kaposi’s sarcoma-associated herpesvirus. J. Virol. 2001, 75, 7882–7892. [Google Scholar] [CrossRef] [PubMed]
  156. Han, S.J.; Hu, J.; Pierce, B.; Weng, Z.; Renne, R. Mutational analysis of the latency-associated nuclear antigen DNA-binding domain of kaposi’s sarcoma-associated herpesvirus reveals structural conservation among gammaherpesvirus origin-binding proteins. J. Gen. Virol. 2010, 91, 2203–2215. [Google Scholar] [CrossRef] [PubMed]
  157. Hu, J.; Yang, Y.; Turner, P.C.; Jain, V.; McIntyre, L.M.; Renne, R. LANA binds to multiple active viral and cellular promoters and associates with the H3K4methyltransferase hSET1 complex. PLoS Pathog. 2014, 10, e1004240. [Google Scholar] [CrossRef] [PubMed]
  158. Krishnan, H.H.; Naranatt, P.P.; Smith, M.S.; Zeng, L.; Bloomer, C.; Chandran, B. Concurrent expression of latent and a limited number of lytic genes with immune modulation and antiapoptotic function by kaposi’s sarcoma-associated herpesvirus early during infection of primary endothelial and fibroblast cells and subsequent decline of lytic gene expression. J. Virol. 2004, 78, 3601–3620. [Google Scholar] [PubMed]
  159. Lan, K.; Kuppers, D.A.; Robertson, E.S. Kaposi’s sarcoma-associated herpesvirus reactivation is regulated by interaction of latency-associated nuclear antigen with recombination signal sequence-binding protein jkappa, the major downstream effector of the notch signaling pathway. J. Virol. 2005, 79, 3468–3478. [Google Scholar] [CrossRef] [PubMed]
  160. Lan, K.; Kuppers, D.A.; Verma, S.C.; Robertson, E.S. Kaposi’s sarcoma-associated herpesvirus-encoded latency-associated nuclear antigen inhibits lytic replication by targeting Rta: A potential mechanism for virus-mediated control of latency. J. Virol. 2004, 78, 6585–6594. [Google Scholar] [CrossRef] [PubMed]
  161. Chang, P.C.; Fitzgerald, L.D.; van Geelen, A.; Izumiya, Y.; Ellison, T.J.; Wang, D.H.; Ann, D.K.; Luciw, P.A.; Kung, H.J. Kruppel-associated box domain-associated protein-1 as a latency regulator for kaposi’s sarcoma-associated herpesvirus and its modulation by the viral protein kinase. Cancer Res. 2009, 69, 5681–5689. [Google Scholar] [CrossRef] [PubMed]
  162. Zhang, L.; Zhu, C.; Guo, Y.; Wei, F.; Lu, J.; Qin, J.; Banerjee, S.; Wang, J.; Shang, H.; Verma, S.C.; et al. Inhibition of KAP1 enhances hypoxia-induced kaposi’s sarcoma-associated herpesvirus reactivation through RBP-Jκ. J. Virol. 2014, 88, 6873–6884. [Google Scholar] [CrossRef] [PubMed]
  163. Yu, Y.; Black, J.B.; Goldsmith, C.S.; Browning, P.J.; Bhalla, K.; Offermann, M.K. Induction of human herpesvirus-8 DNA replication and transcription by butyrate and TPA in BCBL-1 cells. J. Gen. Virol. 1999, 80, 83–90. [Google Scholar] [CrossRef] [PubMed]
  164. Chen, W.Y.; Townes, T.M. Molecular mechanism for silencing virally transduced genes involves histone deacetylation and chromatin condensation. Proc. Natl. Acad. Sci. USA 2000, 97, 377–382. [Google Scholar] [CrossRef] [PubMed]
  165. Bowser, B.S.; DeWire, S.M.; Damania, B. Transcriptional regulation of the K1 gene product of kaposi’s sarcoma-associated herpesvirus. J. Virol. 2002, 76, 12574–12583. [Google Scholar] [CrossRef] [PubMed]
  166. Deng, H.; Song, M.J.; Chu, J.T.; Sun, R. Transcriptional regulation of the interleukin-6 gene of human herpesvirus 8 (kaposi’s sarcoma-associated herpesvirus). J. Virol. 2002, 76, 8252–8264. [Google Scholar] [CrossRef] [PubMed]
  167. Chang, P.J.; Shedd, D.; Miller, G. Two subclasses of kaposi’s sarcoma-associated herpesvirus lytic cycle promoters distinguished by open reading frame 50 mutant proteins that are deficient in binding to DNA. J. Virol. 2005, 79, 8750–8763. [Google Scholar] [CrossRef] [PubMed]
  168. Chen, J.; Ye, F.; Xie, J.; Kuhne, K.; Gao, S.J. Genome-wide identification of binding sites for kaposi’s sarcoma-associated herpesvirus lytic switch protein, Rta. Virology 2009, 386, 290–302. [Google Scholar] [CrossRef] [PubMed]
  169. Ziegelbauer, J.; Grundhoff, A.; Ganem, D. Exploring the DNA binding interactions of the kaposi’s sarcoma-associated herpesvirus lytic switch protein by selective amplification of bound sequences in vitro. J. Virol. 2006, 80, 2958–2967. [Google Scholar] [CrossRef] [PubMed]
  170. Toth, Z.; Maglinte, D.T.; Lee, S.H.; Lee, H.R.; Wong, L.Y.; Brulois, K.F.; Lee, S.; Buckley, J.D.; Laird, P.W.; Marquez, V.E.; et al. Epigenetic analysis of KSHV latent and lytic genomes. PLoS Pathog. 2010, 6, e1001013. [Google Scholar] [CrossRef] [PubMed]
  171. Rossetto, C.C.; Pari, G.S. Kaposi’s sarcoma-associated herpesvirus noncoding polyadenylated nuclear RNA interacts with virus- and host cell-encoded proteins and suppresses expression of genes involved in immune modulation. J. Virol. 2011, 85, 13290–13297. [Google Scholar] [CrossRef] [PubMed]
  172. Kang, H.; Wiedmer, A.; Yuan, Y.; Robertson, E.; Lieberman, P.M. Coordination of KSHV latent and lytic gene control by CTCF-cohesin mediated chromosome conformation. PLoS Pathog. 2011, 7, e1002140. [Google Scholar] [CrossRef] [PubMed]
  173. Tempera, I.; Lieberman, P.M. Chromatin organization of gammaherpesvirus latent genomes. Biochim. Biophys. Acta 2010, 1799, 236–245. [Google Scholar] [CrossRef] [PubMed]
  174. Li, D.J.; Verma, D.; Mosbruger, T.; Swaminathan, S. CTCF and RAD21 act as host cell restriction factors for kaposi’s sarcoma-associated herpesvirus (KSHV) lytic replication by modulating viral gene transcription. PLoS Pathog. 2014, 10, e1003880. [Google Scholar] [CrossRef] [PubMed]
  175. Chang, P.C.; Cheng, C.Y.; Campbell, M.; Yang, Y.C.; Hsu, H.W.; Chang, T.Y.; Chu, C.H.; Lee, Y.W.; Hung, C.L.; Lai, S.M.; et al. The chromatin modification by SUMO-2/3 but not SUMO-1 prevents the epigenetic activation of key immune-related genes during kaposi’s sarcoma associated herpesvirus reactivation. BMC Genom. 2013, 14. [Google Scholar] [CrossRef] [PubMed]
  176. Garcia-Dominguez, M.; Reyes, J.C. SUMO association with repressor complexes, emerging routes for transcriptional control. Biochim. Biophys. Acta 2009, 1789, 451–459. [Google Scholar] [CrossRef] [PubMed]
  177. Chang, P.C.; Izumiya, Y.; Wu, C.Y.; Fitzgerald, L.D.; Campbell, M.; Ellison, T.J.; Lam, K.S.; Luciw, P.A.; Kung, H.J. Kaposi’s sarcoma-associated herpesvirus (KSHV) encodes a SUMO E3 ligase that is SIM-dependent and SUMO-2/3-specific. J. Biol. Chem. 2010, 285, 5266–5273. [Google Scholar] [CrossRef] [PubMed]
  178. Yang, W.S.; Hsu, H.W.; Campbell, M.; Cheng, C.Y.; Chang, P.C. K-bZip mediated SUMO-2/3 specific modification on the KSHV genome negatively regulates lytic gene expression and viral reactivation. PLoS Pathog. 2015, 11, e1005051. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic representation of the different applications of “Next-generation” DNA sequencing technologies.
Figure 1. Schematic representation of the different applications of “Next-generation” DNA sequencing technologies.
Viruses 08 00092 g001
Table 1. New viral features in the KSHV genome including, Upstream and Small ORFs, Alternate start codon usage and Internal ORFs, derived from Next-generation sequencing studies adapted from KSHV 2.0 [57]. All annotation coordinates are based on GQ 994935.1.
Table 1. New viral features in the KSHV genome including, Upstream and Small ORFs, Alternate start codon usage and Internal ORFs, derived from Next-generation sequencing studies adapted from KSHV 2.0 [57]. All annotation coordinates are based on GQ 994935.1.
GeneStrandStart–StopSmall Open Reading Frame (sORF)Upstream Open Reading Frame (uORF)New Features
ORF6.1+3027–3203-uORF-
ORF6.2+3150–3203-uORF-
ORF10.1+14,451–14,531-uORF-
ORF10.2+15,574–15,756--Internal ORF
ORF11.1+15,633–15,722-uORF-
ORF11.2+15,648–15,722-uORF-
ORF11.3+15,693–15,722-uORF-
ORF11.4+15,745–15,756-uORF-
ORF11.5+15,926–15,991--Internal ORF
vIL6.617,873–17,862-uORF-
vIL6.517,915–17,877-uORF-
vIL6.418,003–17,902-uORF-
vIL6.318,047–17,985-uORF-
vIL6.218,086–18,057-uORF-
vIL6.118,116–18,057-uORF-
ORFK3A19,128–18,589--Internal ORF
ORF70A21,099–20,038--Alternate start
ORFK4A21,820–21,743-uORF-
ORFK4.1a22,517–22,416sORF--
ORFK4.1d22,610–22,545sORF--
ORFK4.1e22,653–22,627sORF--
ORFK4.1c22,806–22,723sORF--
ORFK4.1b22,850–22,545sORF--
1.4KbB+24,871–24,915sORF--
1.4KbC+24,921–25,058sORF--
ORFK5.126,569–26,555-uORF-
ORFK6.127,647–27,615-uORF-
ORFK6A27,443–27,087--Alternate start
ORFK6B27,422–27,087--Alternate start
PAN1.1+28,655–28,768sORF--
PAN1.2+28,831–28,965sORF--
PAN1.3+28,888–28,965sORF--
ORF20A35,322–34,429--Internal ORF
ORF20B35,202–34,429--Internal ORF
ORF21.1+35,151–35,177-uORF-
ORF25.1+42,345–42,380-uORF-
ORF28.1+48,758–48,811-uORF-
ORF30.1+50,317–50,358-uORF-
ORF34.1+54,399–54,485-uORF-
ORF35.1+55,419–55,445-uORF-
ORF35.2+55,442–55,474-uORF-
ORF37.1+57,040–57,126-uORF-
ORF38.1+58,251–58,259-uORF-
ORF38.2+58,455–58,589-uORF-
43.1-AS+63,214–63,228sORFs--
ORF43.2-AS+63,254–63,295sORFs--
ORF45.168,447–68,364-uORF-
ORF49.172,425–72,384-uORF-
ORF50AS74,222–74,130sORF--
k8.1 short+75,890–75,971--Internal ORF
ORF54A+77,552–78,439--Alternate start
ORF55.179,501–79,340-uORF-
ORF57A+81,464–83,453--Splice variant
ORF61.2100,071–100,018-uORF-
ORF61.1100,095–100,018-uORF-
ORF62B101,019–100,018--Alternate start
ORF62A101,061–100,018--Alternate start
ORF65.1112,321–112,289-uORF-
ORF69.1+116,043–116,138-uORF-
Kaposin C2119,084–117,738--Alternate start
ORF72.1124,182–124,108-uORF-
ORF75.1134,809–134,729-uORF-
ORF75.2134,894–134,817-uORF-
K15.1135,938–135,846sORF--

Share and Cite

MDPI and ACS Style

Strahan, R.; Uppal, T.; Verma, S.C. Next-Generation Sequencing in the Understanding of Kaposi’s Sarcoma-Associated Herpesvirus (KSHV) Biology. Viruses 2016, 8, 92. https://doi.org/10.3390/v8040092

AMA Style

Strahan R, Uppal T, Verma SC. Next-Generation Sequencing in the Understanding of Kaposi’s Sarcoma-Associated Herpesvirus (KSHV) Biology. Viruses. 2016; 8(4):92. https://doi.org/10.3390/v8040092

Chicago/Turabian Style

Strahan, Roxanne, Timsy Uppal, and Subhash C. Verma. 2016. "Next-Generation Sequencing in the Understanding of Kaposi’s Sarcoma-Associated Herpesvirus (KSHV) Biology" Viruses 8, no. 4: 92. https://doi.org/10.3390/v8040092

APA Style

Strahan, R., Uppal, T., & Verma, S. C. (2016). Next-Generation Sequencing in the Understanding of Kaposi’s Sarcoma-Associated Herpesvirus (KSHV) Biology. Viruses, 8(4), 92. https://doi.org/10.3390/v8040092

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop