KSHV 2.0: A Comprehensive Annotation of the Kaposi's Sarcoma-Associated Herpesvirus Genome Using Next-Generation Sequencing Reveals Novel Genomic and Functional Features: e1003847

Abstract

Productive herpesvirus infection requires a profound, time-controlled remodeling of the viral transcriptome and proteome. To gain insights into the genomic architecture and gene expression control in Kaposi's sarcoma-associated herpesvirus (KSHV), we performed a systematic genome-wide survey of viral transcriptional and translational activity throughout the lytic cycle. Using mRNA-sequencing and ribosome profiling, we found that transcripts encoding lytic genes are promptly bound by ribosomes upon lytic reactivation, suggesting their regulation is mainly transcriptional. Our approach also uncovered new genomic features such as ribosome occupancy of viral non-coding RNAs, numerous upstream and small open reading frames (ORFs), and unusual strategies to expand the virus coding repertoire that include alternative splicing, dynamic viral mRNA editing, and the use of alternative translation initiation codons. Furthermore, we provide a refined and expanded annotation of transcription start sites, polyadenylation sites, splice junctions, and initiation/termination codons of known and new viral features in the KSHV genomic space which we have termed KSHV 2.0. Our results represent a comprehensive genome-scale image of gene regulation during lytic KSHV infection that substantially expands our understanding of the genomic architecture and coding capacity of the virus.

Abstract

Productive herpesvirus infection requires a profound, time-controlled remodeling of the viral transcriptome and proteome. To gain insights into the genomic architecture and gene expression control in Kaposi's sarcoma-associated herpesvirus (KSHV), we performed a systematic genome-wide survey of viral transcriptional and translational activity throughout the lytic cycle. Using mRNA-sequencing and ribosome profiling, we found that transcripts encoding lytic genes are promptly bound by ribosomes upon lytic reactivation, suggesting their regulation is mainly transcriptional. Our approach also uncovered new genomic features such as ribosome occupancy of viral non-coding RNAs, numerous upstream and small open reading frames (ORFs), and unusual strategies to expand the virus coding repertoire that include alternative splicing, dynamic viral mRNA editing, and the use of alternative translation initiation codons. Furthermore, we provide a refined and expanded annotation of transcription start sites, polyadenylation sites, splice junctions, and initiation/termination codons of known and new viral features in the KSHV genomic space which we have termed KSHV 2.0. Our results represent a comprehensive genome-scale image of gene regulation during lytic KSHV infection that substantially expands our understanding of the genomic architecture and coding capacity of the virus.

Author Summary

Kaposi's sarcoma-associated herpesvirus (KSHV) is a cancer-causing agent in immunocompromised patients that establishes long-lasting infections in its hosts. Initially described in 1994 and extensively studied ever since, KSHV molecular biology is understood in broad outline, but many detailed questions are still to be resolved. After almost two decades, specific aspects pertaining to the organization of the KSHV genome as well as the fate of the viral transcripts during the productive stages of infection remain unexplored. Here we use a systematic genome-wide approach to investigate changes in gene and protein expression during the productive stage of infection known as the lytic cycle. We found that the viral genome has a large coding capacity, capable of generating at least 45% more products than initially anticipated by bioinformatic analyses alone, and that it uses multiple strategies to expand its coding capacity well beyond what is determined solely by the DNA sequence of its genome. We also provide an expanded and highly detailed annotation of known and new genomic features in KSHV. We have termed this new architectural and functional annotation KSHV 2.0. Our results indicate that viral genomes are more complex than anticipated, and that they are subject to tight mechanisms of regulation to ensure correct gene expression.

Citation: Arias C, Weisburd B, Stern-Ginossar N, Mercier A, Madrid AS, et al. (2014) 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 10(1): e1003847. doi:10.1371/journal.ppat.1003847

Editor: Dirk P. Dittmer, University of North Carolina at Chapel Hill, United States of America

Received: June 21, 2013; Accepted: October 20, 2013; Published: January 16, 2014

Copyright: © 2014 Arias et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: CA, ASM and AM were supported by a Novartis Presidential Postdoctoral fellowship. NSG was supported by a human frontier science program postdoctoral fellowship. JSW and DG were supported by the Howard Hughes Medical Institute. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Kaposi's sarcoma-associated herpesvirus (KSHV) is a member of the gamma-herpesvirus family and the etiologic agent of Kaposi's sarcoma, primary effusion lymphoma (PEL), and multicentric Castleman's disease [1], [2]. This human pathogen, initially identified in Kaposi's sarcoma lesions from AIDS patients, has been extensively studied since its discovery and isolation in 1994 [3]. Shortly thereafter, the KSHV genome, a dsDNA molecule of ~165 kb, was sequenced from the lymphoid cell line BC-1, allowing the in silico annotation of open reading frames (ORFs) that fit the following criteria: (1) they start with a canonical initiator AUG codon and (2) they encode polypeptides larger than 100 amino acids (aa). Many of these ORFs had functional homologues in herpesvirus saimiri (HVS), a gamma-herpesvirus related to KSHV [4]. This study identified a total of 81 such viral ORFs, and except for the more recent addition of microRNAs, non-coding RNAs, and a few small ORFs [5]-[8], the genome map of KSHV has changed little ever since. Gene expression profiling of KSHV transcripts using northern blots, custom oligonucleotide microarrays and real time PCR arrays have demonstrated extensive transcription of the viral genome, hinting at a complex transcriptional profile [6], [9], [10] (unpublished data). More recently, proteomic studies of KSHV-infected cells have assessed the expression of many of the predicted ORFs [11]. However, and in spite of all the aforementioned efforts, a detailed understanding of the genomic architecture, translational state, and biological functions of KSHV gene products remains incomplete.

In an attempt to extend our current knowledge of the coding capacity of the KSHV genome during the productive stage of infection, we employed an unbiased functional genomics approach to study the transcription and translation profiles of lytic KSHV using mRNA-sequencing (mRNA-Seq), ribosome footprinting (Ribo-Seq), and genomic DNA sequencing (DNA-Seq). When combined, these methods provide a comprehensive, high-resolution view of gene regulation and expression dynamics [12]-[14].

By employing these techniques in parallel, we have generated a state-of-the-art annotation of the KSHV genome. Our approach confirms the presence and timing of expression of the majority of previously annotated ORFs, while revealing several novel and, in some cases, unexpected genomic features including ribosome protection of non-coding RNAs, new splice variants, and a plethora of upstream and small ORFs. In addition, we have confirmed and expanded the annotation of transcription start sites, polyadenylation sites, and initiation/termination codons of multiple known ORFs. Our analyses have also uncovered new instances of viral mRNA editing, strongly hinting at a new layer of viral gene regulation during reactivation. The wealth of information generated by integrating the data obtained from our combined methods has expanded our understanding of the viral genome architecture and dynamics, revealing a surprising coding capacity of KSHV that goes well beyond what was initially described based on its genome sequence alone.

Results

mRNA-Seq and Ribo-Seq reveal the architecture of the coding and non-coding viral transcriptome at a single-nucleotide resolution

The life cycle of KSHV can be separated in two very distinct stages: the dormant state known as latency and the productive state referred to as the lytic cycle [2]. While viral gene expression in latency is limited and most of the genome is silent, the lytic cycle is a transcriptionally dynamic state where the timing of gene expression is tightly regulated to ensure the ordered synthesis of viral products [6], [15], [16]. We sought to study the kinetics of latent and lytic viral transcription in detail, as well as the translational fate of newly synthesized mRNAs. To this end, we employed a system developed by our group that allows the study of the KSHV lytic cycle in a tightly-controlled manner [17]. This system comprises the epithelial iSLK-219 cell line, which is latently infected with a heterologous KSHV strain (see below) and harbors a doxycycline (Dox) inducible transgene encoding the viral transcription factor RTA (replication and transcriptional activator). The exogenous expression of RTA by Dox treatment in iSLK-219 cells is sufficient to induce the lytic reactivation of latent KSHV. Notably, latency in iSLK-219 cells is very strict with less than 0.1% of the cells showing lytic markers in the absence of induction [17]. This is the principal experimental advantage of SLK cells, allowing the study of KSHV latency in the near total absence of contaminating lytically infected cells. The viral strain in iSLK-219 is the recombinant KSHV.219, which encodes a constitutive GFP reporter as well as an RTA-inducible RFP reporter in the viral genome, thereby facilitating the monitoring of viral reactivation [18].

To finely resolve the transcriptional profile and the ribosome occupancy of viral mRNAs, we induced iSLK-219 cells with Dox for 0, 8, 24, 48 and 72 hr (Figure 1A). We evaluated KSHV lytic reactivation by epifluorescence microscopy analysis of GFP and RFP expression, as well as by immunodetection of viral products and quantification of viral DNA replication (Figure 1B, Figure S1). In iSLK-219, lytic DNA replication, the traditional border between early and late times, commences at ~48 hr post induction (hpi) (Figure S1B). The selected time points represent the different stages of the lytic cycle, known as latent (0 hr), immediate early- (8 hr), delayed early- (24 hr), and late-lytic (48 and 72 hr) (Figure 1A). At each time point, we recovered polyadenylated RNA (mRNA) and 3 sets of ribosome footprints (described in Materials and Methods). To map actively elongating ribosomes on viral transcripts, we isolated ribosome footprints from cells treated with cycloheximide (CHX), a translation inhibitor that binds the ribosomal-E-site and arrests elongating ribosomes [19]. In the same manner, we mapped initiating ribosomes by treating cells with harringtonine (Harr), a translational inhibitor that binds the 60S subunit and hinders the progression of the initiating ribosome, causing ribosomes to stall at translation start sites [20]. Finally we mapped releasing ribosomes accumulating at the stop codon in samples that were not treated with any drug as previously described [14]. We then constructed Illumina-compatible libraries from fragmented and size-selected mRNA segments (40-100 nt), or ribosome protected RNA (ribosome footprints ~30 nt in length, Figure S2A) following the standard ribosome profiling protocol previously described [14], [21]. The libraries were deep-sequenced and the resulting reads aligned to the KSHV genome (GQ994935). As expected, the number of reads aligning to the KSHV genome increased as the lytic cycle progressed (Figure S2B). To annotate viral splice junctions, we used two splice junction mapping tools; TopHat and HMMSplicer [22], [23]. With these tools, we detected the majority of the known splice junctions and discovered 7 new events including one at the 3' end of ORF57. Lastly, we annotated putative ORFs by training a support vector machine (SVM) to identify translation initiation sites throughout the KSHV genome based on characteristic peaks within the harringtonine Ribo-Seq data. The list of ORFs produced by extending each of the putative initiation sites to the next in-frame stop codon (taking into account any intervening splice junctions) was then finalized through manual curation.

[Figure omitted, see PDF]

Figure 1. High-resolution mapping of KSHV genomic features.

(A-B) Strategy for identification of transcripts and ribosome footprints in iSLK-219 cells. Latent iSLK-219 cells were induced to enter the lytic cycle by expression of the doxycycline-inducible transcription factor RTA for 8, 24, 48 and 72 hr. mRNA and ribosome footprints were isolated from cells and the purified RNA was deep sequenced. (C) Ribosome occupancy and mRNA profiles of the late lytic gene K8.1 at 48 hr post induction (hpi).

Figure 1C shows an example of the data obtained with our combined approach. In this case the read coverage from mRNA-Seq and Ribo-Seq (CHX and Harr) libraries for the late-lytic viral gene K8.1, one of the best-mapped genes in KSHV. The data clearly illustrates the single-nucleotide resolution and high-coverage of our methods which here allow the delineation of transcription start and end sites, splice junctions, and coding region boundaries. Notably, the coordinates derived from our combined approach correspond precisely to those previously reported for this gene [24], [25], providing strong validation of our methodology. Together, the data obtained using our multipronged approach generated a high-resolution map of the viral genome architecture.

KSHV 2.0: A revised, high-resolution functional and temporal annotation of the KSHV genome

mRNA-Seq and Ribo-Seq allowed us to perform an unabridged temporal analysis of viral gene expression coupled to a blueprint of the viral episome architecture, granting the opportunity to develop a revised version of the KSHV genome annotation, which we have designated KSHV 2.0 (Figure 2 and Tables 1, 2 and 3). In KSHV 2.0 we annotate the coordinates for 49 viral transcripts and 70 ORFs, as well as those for non-coding RNAs, polyadenylation signals, and splice junctions. In addition, KSHV 2.0 incorporates information pertaining to the timing of expression of the aforementioned elements. Remarkable novel features of KSHV 2.0 include a set of 50 ribosome-loaded segments not previously annotated as bona fide ORFs because of their small size (3-100 aa) and the use of non-canonical start sites (Figure 2B). Together with novel peptide isoforms, and splice variants, these short and upstream ORFs (sORFs and uORFs, respectively) increase the coding repertoire of KSHV by more than 45% and add a new level of potential gene regulation to an already complex landscape. The novel features annotated in this study are summarized in file S1 and can be visualized using the mochiview database in file S2 [26]. In spite of the comprehensive annotation generated for KSHV 2.0, some known features of the viral genome were not detected or could not be rigorously assigned (Figure 2B), due to ambiguities generated by regions of low sequencing coverage, overlapping transcription and translation, and cell line specific patterns of gene expression.

[Figure omitted, see PDF]

Figure 2. KSHV 2.0: A high-resolution functional genome map of KSHV.

(A) Circular map of the KSHV genome. The outer circle represents the genomic coordinates. The concentric circles illustrate, from the outermost to the innermost: viral ORFs (Red sense, blue antisense), sORFs (purple) and uORFs (orange), transcripts (maroon sense, dark blue antisense) and timing of expression (from yellow to dark green: latent, early lytic, delayed early lytic, late lytic). (B) Linear map of the KSHV genome showing the details of genomic features illustrated in A. ORF29, ORF32 and ORF36 are depicted as annotated ORFs as we detected the presence of ribosomes in their coding sequences. The exact boundaries of translation for these ORFs could not be determined. See also Tables 1, 2 and 3.

Table 1. Transcription and translation coordinates, timing of expression, and function of latent, early lytic and delayed early lytic ORFs.

Table 2. Transcription and translation coordinates, timing of expression, and function of late lytic ORFs.

Table 3. Transcription and translation coordinates, and function of primary ORFs without determined time of expression.

Inspection of KSHV 2.0 reveals three prominent features: (1) coding and non-coding elements are densely packed in the episome, (2) multiple strategies are used to increase its polypeptide repertoire, including splicing, mRNA editing, and alternative start codon use, and (3) sORFs and uORFs populate many regions of the viral genome. The specific transcriptional and translational features of KSHV 2.0 are discussed in detail below.

mRNA-Seq reveals a dense viral transcriptome encoded in the KSHV genome

The transcriptional capacity of the KSHV genome has been traditionally studied using northern blotting and gene expression profiling with oligonucleotide microarrays [6], [9], [27], [28]. While these studies have exposed many features of the viral transcriptome, the limitations of these methods prevent the fine mapping of transcripts, which can require single-nucleotide resolution. For this reason, we performed mRNA-Seq in cells lytically infected with KSHV, to explore the transcriptional landscape of the KSHV genome, resolve the boundaries of viral messages, and uncover novel cis-regulatory elements including transcription start sites (TSS), polyadenylation signals (PAS), and splice junctions.

Taking advantage of the peaks visible at the 5' ends of transcripts in mRNA profiles (Figure 1C), which are a natural consequence of visualizing the 5' ends of fragments produced via random fragmentation of multiple mRNA copies of any given transcript, we mapped 49 TSS upstream of 54 out of the 85 officially annotated genes (See Materials and Methods, Tables 1 and 2). The annotation of the TSS coordinates for the remaining 31 viral genes was impeded by low coverage or the presence of overlapping transcripts. Among the mapped viral genes, the discrepancy between the number of TSS and genes, stems from the existence of bi- and poly-cistronic mRNAs (Figure 3A). Of the 49 TSS mapped, 28 are novel while 21 correspond to annotated transcripts whose TSS were previously characterized. Of the 21 previously documented TSS described in our study, 13 are mapped exactly as in the literature and 8 are located within 50 nucleotides of their reported coordinates as previously resolved by 5' rapid amplification of cDNA ends (5'RACE) (Tables 1 and 2). Interestingly, sequence alignment of the promoter regions corresponding to the TSS unveiled in our analyses, shows the presence of a TATA-box 30 nucleotides upstream of 41 TSS, which remarkably corresponds to the same location of this cis-regulatory element in human promoters. The remaining eight TSS are TATA-less (Table S1). Our observations clearly reflect the strict evolutionary dependence of the pathogen on the host's transcriptional machinery (Figure 3B) [29]-[31].

[Figure omitted, see PDF]

Figure 3. Host cis-regulatory elements are conserved in viral transcripts.

(A) mRNA-Seq of the locus encoding the ORF8/9/10 and 11 genes at 72 hr post reactivation. Three TSS sharing a common PAS were mapped within this polycistronic locus. (B) Sequence analysis of 50 nucleotides flanking the predicted TSS site for 49 viral genes. Note the presence of the TATA box 30 bp upstream of the TSS. (C) Sequence analysis of 30 nt flanking the 42 polyA sites mapped in KSHV. Note the conserved AA/UUAAA motif 5-20 nt upstream (top panel), and the GU rich region (middle panel) 5-20 nt downstream of the cleavage site. See also Tables S1 and S2.

Like cellular messages, KSHV transcripts are protected by a 5' 7-methylguanosine cap and a 3' poly-adenylate (polyA) tail [32], [33]. To map functional polyadenylation signals (PAS) in the viral genome, we selected the RNA-Seq reads that contained a stretch of 5 or more adenosines at their 3' end and, after trimming this poly-A sequence, aligned the reads to the KSHV genome. The 3' positions of aligned reads were then marked as polyadenylation cleavage sites, except when the genome contained a complementary poly-T stretch at the same location as the poly-A stretch. Using this approach, we mapped 94 putative cleavage sites, corresponding to 42 transcripts and 74 genes (Tables 1, 2 and 3, Table S2). Our data, and recent studies mapping the 3'UTRs of KSHV genes [34], [35], support the existence of bi- and poly-cistronic messages, as well as transcript clusters with distinct transcription start sites (TSS) that end in a common PAS, suggesting the existence of alternative nested promoters upstream of such PAS's (Figure 3A, Table S2).

Interestingly, sequence analysis of a 60-nucleotide window centered on the predicted cleavage site for the polyadenylation machinery confirmed the presence of the canonical AAUAAA motif in 83% of mRNAs (35 out of 42), the alternative AUUAAA motif in 17% (7 out of 42), and the accompanying downstream GU rich element in all of these transcripts (Figure 3C and Table S2). Similar observations were recently reported in genome-wide analyses of polyadenylation sites in PEL cells infected with KSHV [36]. The presence of these conserved elements highlights once more the strict dependence of the virus on host factors that control the RNA processing of Pol II transcripts [37].

Identification of viral, ribosome-protected non-coding RNAs by mRNA-Seq and Ribo-Seq

Besides the densely packed coding regions and regulatory features we annotated, our mRNA-Seq data also show the massive accumulation of sequence reads that map outside of previously annotated coding regions, thus indicating highly permissive transcription of most of the viral genome late in the lytic cycle (Figure S3). Two of these regions correspond to two long transcripts recently discovered by our group, the 10 kb antisense RNA to the latent transcripts (ALT) and the 17 kb K1-ORF11 antisense (K1/11-AS) [6]. Intriguingly, these long RNAs show short regions modestly populated by ribosomes, suggesting they may have a coding potential (Figure S4A and S4B).

It is noteworthy that this observation was not restricted to ALT and K1/11-AS. Surprisingly, our Ribo-Seq data also revealed the presence of ribosomes on the "non-coding" RNA PAN (polyadenylated nuclear RNA). PAN is the most abundant viral transcript during the lytic cycle and is required for viral gene expression and virion production [33], [38], [39]. Interestingly, and in spite of PAN's reported nuclear localization, we observed initiating ribosomes accumulating at the start codon in the harringtonine treated samples, elongating ribosomes throughout the body of the transcript in the CHX treated samples, and an accumulation of releasing ribosomes at the stop codon in the samples not treated with any translation inhibitor, starting at 8 hr following reactivation and throughout the lytic cycle (Figure 4A, S5A and S5B). Taking in consideration this pattern of ribosome protection, classical of coding regions [14], we identified three predominant sORFs hosted within the PAN transcript: PAN1.1 (37 aa, 28655), PAN1.2 (44 aa, 28831) and PAN1.3 (25 aa, 28888) (Figure 4B). Besides these, we also identified 3 minor sORFs at the 3' end of PAN, with very low ribosome occupancy (data not shown). To evaluate the coding capacity of PAN, we calculated the ribosome release score (RSS) for the three main putative ORFs, PAN1.1, 1.2 and 1.3. The RRS is a metric that takes in consideration that ribosome protection within a coding region ends after the stop codon and that no ribosomes should be present at the 3'UTR of the transcript following an ORF. A recent report by Guttman et al. indicates that the RSS provides an indirect measure of translation that allows the reliable differentiation between coding and non-coding transcripts [40]. The RSS calculated for PAN1.1, 1.2 and 1.3 are comparable to those of known coding RNAs and are similar to scores previously determined for small ORFs within mammalian transcripts [40], further supporting the translation potential of PAN (Figure S5D, E). Notwithstanding the low translation efficiency of the three major PAN sORFs (0.05 at 8 hpi to 0.2 at 72 hpi), our mRNA-Seq and Ribo-Seq data suggest that, owing to the significant accumulation of the PAN transcript during the lytic cycle, the putative peptides encoded in these sORFs could be quite abundant. In fact, PAN RNA represented up to 92% of the total viral mRNA-Seq reads and the ribosome-protected RNA corresponding to the small PAN peptides represented up to 1.7% of the total cycloheximide Ribo-Seq reads (Table S3). Thus, our data strongly indicate that this transcript is available for ribosome binding, and that in addition to its documented functions as a "non-coding" RNA, PAN may also be a presumptive coding RNA. It is important to note that the putative coding regions for PAN are overlapping with the ORFK7 transcript. However, close inspection of the ribosome accumulation at the start codon of ORFK7 indicates that in order for the PAN peptides to be encoded by the K7 transcript, the translation efficiency of the internal peptides would need to be 1000 to 10000 times more efficient than that of the main ORF K7 (Figure S5C). Based on this observation and the vast number of reads seen for the mRNA and ribosome protected fragments in PAN1.1, 1.2 and 1.3, we conclude that these putative coding regions are within the PAN transcript.

[Figure omitted, see PDF]

Figure 4. mRNA-Seq reveals ribosome-protected non-coding RNAs in KSHV.

(A) The PAN transcript is protected by ribosomes during the lytic cycle. mRNA-Seq and Ribo-Seq of the PAN gene at 48 hr post reactivation. Ribo-Seq in cells treated with CHX (blue) or Harr (teal) show the accumulation of ribosomes on three distinct coding regions of this transcript. (B) PAN codes for three putative sORFs, PAN1.1 (37aa), PAN1.2 (44aa), and PAN1.3 (25aa), which can be identified as ribosome-populated regions (in capital letters). The predicted signal peptide in PAN1.1 is underlined. (C) The K5/K6 antisense transcript is devoid of nucleosomes following lytic reactivation. mRNA-Seq and Ribo-Seq of the genomic region from 24200 bp to 29700 bp at 48 hr post induction. The light blue solid line corresponds to the K5/K6 antisense RNA. Solid black arrows indicate transcripts and grey arrowheads indicate coding regions. (D) Northern blot for the same region depicted in (C). The A and B probes are indicated as short broken lines in (C). PFA: phosphonoformate. See also Figure S5 and Table S3.

We hypothesize that despite their minuscule size, the putative peptides encoded by PAN may be functional. Indeed, the number and characterized functions of such small peptides are continuously increasing, and there is overwhelming evidence in other viral systems, as well as in eukaryotic cells, for the abundance and relevance of small peptides [14], [21], [41]. Encouraged by these findings, we started to look for a possible function of the predicted peptides encoded by PAN1.1, 1.2 and 1.3. To such end, we used bioinformatics tools that included diverse motif finding and peptide-function prediction engines [42], [43]. Surprisingly, multiple independent analyses predicted a putative signal peptide in PAN1.1 (Figure 4B), thus suggesting that this peptide may traverse the secretory apparatus. Importantly, the sequence of PAN1.1 and the other small peptides predicted within PAN show 100% conservation at the nucleotide level between different isolates of KSHV (data not shown). Further studies of these putative gene products are underway.

Three additional regions show extensive mRNA-Seq coverage, particularly at 72 hr post reactivation. These are the antisense transcripts corresponding to ORFK5/K6, as well as the antisense transcripts for ORFK9-ORF62 and ORFK2 (vIL6)-ORFK4.2 (Figure 4C, data not shown). We confirmed the existence of a 6 kb RNA antisense to ORFK5/K6 by northern blotting (Figure 4D). This antisense RNA, which we have denominated K5/K6-AS, corresponds to the T6.1 RNA previously described by Taylor et al. [28]. Strikingly, K5/K6-AS is devoid of any initiating or elongating ribosomes, and therefore may represent a bona-fide long non-coding RNA in KSHV that is inherently distinct from PAN, ALT and K1-11AS (Figure 4).

mRNA-Seq and Ribo-Seq unveil strategies used by the virus to expand its peptide repertoire

Like other herpesviruses, KSHV makes widespread use of the cellular mRNA splicing machinery [14], [44]. To confirm known splice junctions and discover new ones, we annotated the possible viral splice junctions in a genome-wide fashion by employing HMMSplicer and TopHat on the mRNA sequences that did not align to either the viral or the human genome [22], [23]. Our results confirmed the presence of 27 splice junctions, corresponding to one or more introns in 17 viral genes (20% of genes). These included the well-characterized splice variants observed in ORF50 and ORF57, as well as the multiple splice variants of K8 and K8.1 (Figure 1 and Table S4). The coordinates of the splice junctions annotated using our experimental data confirm those from previous reports (Table S4), affirming the reliability of our combined methods.

Notably, our data not only correctly annotated known splice junctions but revealed 7 novel ones, thus increasing the number of splicing events from 20 to 27 (Table S4). One such splice junction is located at the 3' end of the ORF57 transcript. ORF57 is a well-characterized KSHV protein thought to be an activator of mRNA maturation and transport, enhancing viral gene expression [45], [46]. Our data support the existence of canonical splice donor and acceptor sites in the new predicted junction (Figure S6A), which give rise to a novel splice variant of ORF57 in which the truncation of a second exon results in the accumulation of ribosomes on a previously uncharacterized third exon (Figure 5A,B). We confirmed the second splicing event of the ORF57 transcript by end-point PCR in iSLK (iSLK-219) and lymphatic endothelial cells (LEC-219) infected with recombinant KSHV.219, but not in infected B cells (BCBL-1) (Figure 5C and S6B).

[Figure omitted, see PDF]

Figure 5. KSHV employs multiple strategies to expand and regulate its coding repertoire.

(A-C) Novel splice variant of ORF57 detected by mRNA-Seq. Tracks for mRNA-Seq and Ribo-Seq in CHX treated cells (48 hpi) show the new, shorter-splice variant of ORF57, ORF57A (red box). Splicing of the second intron in ORF57 removes the UAA stop codon, resulting in the generation of a 33 aa C-terminal exon. (B) Zoom in of the region highlighted in (A). Note that ribosomes accumulate downstream of the ORF57 stop codon, indicating that the new exon may be translated. (C) Verification of ORF57A splicing by end-point PCR in lytic iSLK-219 (72 hpi) and LEC-219 cells. The primers used for amplification flank the intron boundaries and are shown in panel (A) as solid black lines. The high molecular weight product (~800 bp) corresponds to the unspliced form of ORF57. The low molecular weight product (165 bp), corresponds to ORF57A, from which the second intron has been removed. (D-E) mRNA editing in the RTA and Kaposin transcripts. Sanger sequencing of RTA (D) and Kaposin (E) genomic DNA (top panel) and cDNAs (48 hpi) (bottom panel) from iSLK-219 cells. The edited nucleotide is bold, italic font. (F-G) ORF54 translation starts from two alternative initiation codons. (F) Ribo-Seq tracks for cycloheximide (CHX blue) and harringtonine (Harr teal) treated cells at 72 hpi. The dotted lines indicate the two translation initiation sites used for ORF54 translation. (G) Immunoblot for ORF54 in latent and lytic (48 hpi) iSLK cells, infected with Wt KSHV or an ORF54 knock-out virus (ΔORF54). *Non-specific product. See also Figures S6, S7 and S8.

The second splicing event in the ORF57 transcript results in the removal of a 571 nucleotide fragment encoding amino acids 266-455 within the second exon and leads to the loss of the canonical UAA stop codon of ORF57, resulting in the generation of a novel isoform, here named ORF57A, with a different C-terminus that contains 33 amino acids (Figure 5B). Interestingly, the stop codon of ORF57A (position 83464) is located downstream of the canonical polyA cleavage site in ORF57 (position 83453), suggesting the presence of a transcript with an extended 3' end. We confirmed the existence of such longer mRNAs in iSLK-219 and BCBL-1 cells by PCR using a primer set annealing within the ORF57 and ORF57A coding regions, and downstream of the annotated polyA cleavage site (Figure S6C). ORF57A, the new alternative splice variant of ORF57, is 299 amino acids in length and lacks the C-terminal leucine zipper (aa 343-364), the second arginine-glycine-glycine rich domain (RGG) (aa 372-374), the zinc finger domain (aa 423-432) and the glycine-leucine-phenylalanine-phenylalanine(GLFF) domain (aa 447-450) (Figure 5A). While the expression of a truncated form of ORF57 could have functional implications, it is important to note that the ORF57A splice variant was detected only in cells infected with the recombinant KSHV.219 virus (iSLK-219 and LEC-219), and could reflect a secondary effect of the insertion of the GFP/RFP reporter cassette downstream of ORF57 [18], leading to the activation of this cryptic splice site within the ORF57 transcript. This observation demonstrates that elements inserted within the viral genome, even in regions that are seemingly devoid of regulatory/functional elements, may not be inert and could have repercussions on viral gene expression and/or function.

A second posttranscriptional mechanism employed by KSHV to expand its coding capacity is mRNA editing. The post-transcriptional recoding of RNA results in single nucleotide discrepancies between the genomic and transcript sequences [47], [48]. By comparing our mRNA-Seq and DNA-Seq data sets, we found 6 instances of mRNA editing in KSHV in two or more time points (Table S5). Three of such editing events include the previously reported A-to-G transition in genomic position 117,809 within the transcript encoding Kaposin [49] and two novel G-to-A transitions at genomic positions 72,795 and 74,281 of the mRNA encoding RTA, leading to amino acid substitutions in the corresponding encoded polypeptides (Table S5). We confirmed these editing events by end-point PCR amplification of cDNA followed by Sanger sequencing (Figure 5D, E, data not shown). Interestingly, we note that the Kaposin message is edited starting at 24 hr following reactivation, and that the relative amount of edited transcript increases dramatically as the lytic cycle progresses, leading to a highly penetrant A-to-G transition at 72 hr (Figure 5E, Figure S7A and B). The surge in the levels of Kaposin mRNA editing is concomitant with the up-regulation of all isoforms of the adenosine deaminase acting on RNA (ADAR), the enzyme implicated in A-to-I editing (Figure S7C). Our observations indicate that ADAR is at least partially insensitive to the generalized host shutoff mediated by the viral endonuclease SOX (Figure S7D) [50]. The consequence(s) of higher levels of ADAR on host mRNA and other viral transcripts, if any, remain to be determined. In addition to investigating the mechanisms of regulation and the activity of ADAR during lytic infection, it would be of great interest to ascertain the biological impact of the A638T substitution in RTA. The post-translational modification prediction tool NetPhosK [51] suggests that such a mutation improves the sequence context for S634 and S636 phosphorylation in the C-terminus of the protein [52]. The biochemical and functional consequences of this mRNA editing in RTA are yet to be determined and will be the focus of future studies. Two of the predicted events, a U-to-G transversion in position 6144 within ORF6 and a G-to-U transversion in position 96434 within ORF59, were not affirmed by Sanger sequencing, highlighting the importance of validation of putative mRNA editing sites identified through next-generation sequencing (data not shown).

A third mechanism used by KSHV to increase the coding capacity of its viral genome is independent of transcriptional control and involves the manipulation of translation. Ribo-Seq allowed us to accurately map most of the annotated viral ORFs while affording the opportunity to discover several dozen undocumented peptides and putative protein isoforms (Figure 2B). Our data show five coding regions that are of particular interest in that regard. ORF70, ORF K6, ORF54, ORF62 and Kaposin exhibit a remarkable accumulation of initiating ribosomes on multiple in-frame translation start sites, strongly arguing in favor of the presence of at least two protein variants for each one of these ORFs (Figure 5F, Tables 1 and 2, Table S6, Figure S8, File S1 and data not shown). We examined the expression of one of these proteins, ORF54 by immunoblotting in lysates from latent and lytic iSLK cells infected with Wt KSHV or an ORF54 deleted virus [53]. In perfect agreement with our Ribo-Seq data, we detected two isoforms of ORF54 using an antibody directed against the C-terminus of the protein. These migrate at ~ kDa (318 aa -ORF54) and 32 kDa (291 aa-ORF54A) in denaturing SDS-PAGE gels, consistent with our finding of 2 polypeptides that share a common C-terminal domain but possess distinct amino-termini owing to the usage of alternative translation initiation sites. Furthermore, our sequencing data indicate that the previously uncharacterized short form of ORF54 is the most abundant one (Figure 5F), which is also in exact agreement with our immunoblot analysis (Figure 5G). The ORF54 and ORF54A products detected during lytic infection of iSLK cells are also clearly seen in HEK293 cells transfected with C-terminally tagged versions of the gene (data not shown). Taken together, our data affirm that KSHV can selectively use alternative start codons to amplify the peptide repertoire synthesized during the lytic cycle.

Ribo-Seq reveals a cryptic translational regulatory network

The peptide coding capacity of KSHV has been defined globally employing in silico approaches and proteomics studies [4], [11] and at a single-gene level by mutagenesis, epitope-tagging and immunodetection. We sought to obtain a unifying and comprehensive understanding of the viral peptide coding capacity using ribosome profiling. Using Ribo-Seq, we mapped most of the previously annotated viral ORFs with precision and, remarkably, discovered 63 new ORFs, representing a higher than 45% increase of the annotated coding capacity of KSHV to date. The vast majority of these new ORFs encode peptides smaller than 100 amino acids and, in 44% of the cases, peptides that are translated from initiation codons with consensus or near-consensus Kozak sequences [54]. Thus, we have reclassified the coding regions of KSHV into primary ORFs, alternative splice variants, internal ORFs, ORFs with alternative start codons, small (sORFs) and upstream ORFS (uORFs) (Figure 2, Tables 1 and 2, Figure S9B-C, File S1).

We defined sORFs as all of those regions encoding peptides of ~100 aa or less that are not found at the 5' of an annotated viral gene [55]. In total we found 14 sORFs within 6 transcripts (Table S7). Among these, we clearly detected ribosomes populating the 5' end of the ORF50-antisense (50-AS) transcript at 24-48 hr post reactivation (Figure S9B), confirming recent reports that indicate that this mRNA is indeed present in polysomal fractions [7]. While previous transfection-based studies from our lab and others have characterized peptides ranging from 17 to 48 amino acids starting from multiple AUG initiation codons, our Ribo-Seq data indicate that the accumulation of initiating ribosomes in an authentic viral infection involves at least three non-canonical start codons giving rise to small peptides from 8 to 76 aa [7], [8].

The second class of small coding regions revealed by our Ribo-Seq data consists of a group of 36 upstream ORFs (uORFs). These uORFs are present in the leader sequence of annotated ORFs and encode peptides of ~100 aa or less [56]. We noted that uORFs are very numerous and widely distributed across the whole genome (Table S7). In total, 24 genes have between 1 and 6 uORFs that are either in-frame or out-of-frame with the main ORF (Figure S9C). Interestingly, and as has been previously reported for HCMV and mammalian cells [14], [21], 44% of uORFs are translated from a non-canonical start codon and are highly detected at late times during reactivation (Figure S9A, Table S7).

An example of the regulatory capacity of uORFs in KSHV was recently documented by Kronstad and colleagues, who described the functions of two uORFs identified in our Ribo-Seq data as uORF35.1 and uORF35.2 (Figure 6A) [57]. uORF35.1 and uORF35.2 have opposing regulatory functions on the translation of the downstream ORFs ORF35 and ORF36, in a mechanism akin to that described for eukaryotic uORFs regulating cell-stress response genes [58]. These uORFs are located in the 5' leader sequence of the ORF35-36 bicistronic transcript. The uORF35.1 small peptide (8 aa) is in-frame with respect to ORF35, while the uORF35.2 small peptide (10 aa) is an out-of-frame overlapping ORF with respect to ORF35 (Figure 6B). Both of these uORFs inhibit the expression of ORF35, as their deletion promotes accumulation of this protein. However, uORF35.2 has stimulatory effects on the translation of the most 3' gene, ORF36, via a continuous scanning mechanism [57]. These data affirm the existence and functional significance of two of the uORFs identified by Ribo-Seq, and support the reliability of this method for identifying such elements.

[Figure omitted, see PDF]

Figure 6. KSHV encodes functional uORFs.

(A) mRNA-Seq and Ribo-Seq tracks (CHX and Harr) show the accumulation of ribosomes in the regions of uORF35.1 and uORF35.2. These uORFs regulate the expression of ORF35 and ORF36 (adapted from [57]). (B) Start and coding sequences for uORF35.1, uORF35.2 and ORF35. See also Table S7.

The time-dependent phase switch controlling gene expression in KSHV

The phase switch from latency to the lytic cycle is a highly regulated process that requires the temporally controlled expression of genes. Our strategy of mapping transcripts and coding sequences across different stages of the lytic cycle in an RTA-regulatable expression system revealed a built-in timer for viral reactivation that relies on the use of specific TSS during the distinct stages of the viral life cycle. Fine temporal mapping of viral transcription include 4 latent messages, 13 early genes that are expressed starting at 8 hr, 19 genes expressed between 24 hr and 48 hr after reactivation and 38 genes at 48-72 hr following DNA replication (Tables 1, 2 and 3, Figure 7A, Figure 8A and 8B). We took advantage of the restricted latency and protracted lytic cycle observed in iSLK-219 cells when compared to cells of lymphoid origin (unpublished observations) to study the kinetics of viral transcription in much finer detail. As expected, our results show that only a handful of transcripts are expressed during latency, namely the K1-ORF4 bicistronic message, vIL6, Kaposin and the LANA-vCyclin-vFLIP tricistronic transcript. Interestingly, the ribosome profiling of latent cells shows that only vIL6 and LANA are protected by ribosomes (Figure 7A and 7B). Furthermore, we confirmed by immunoblot the presence of LANA and the absence of Kaposin and vCyclin proteins in latently infected cells (Figure 7D and 7E). These observations raised questions about whether the K1, vCyclin and vFLIP proteins might be importantly regulated at the level of translation. Consistent with this, the Kaposin, K1, vCyclin and vFLIP transcripts are abundantly protected by ribosomes upon induction of the lytic cycle, and their cognate proteins can be detected by immunoblotting after such induction (Figure 7C, 7D and 7E). However, subsequent Northern blot analysis revealed the pattern of accumulation of vCyclin and vFLIP transcripts in iSLK differs from that previously observed in B cells (Figure 7F and 7G) [59]. To our surprise, we could detect only the tricistronic (LANA-vCyclin-vFLIP), but not the bicistronic (vCyclin-vFLIP) message, in latent iSLK-219 cells (Figure 7F). However, the bicistronic transcript, which has previously been proposed to be the mRNA for these 2 proteins, is abundantly expressed in lytic iSLK-219 cells (Figure 7F) [60]. This observation, in combination with our mRNA-seq and Ribo-seq data, suggests that vCyclin and vFLIP proteins are indeed mainly expressed from the bicistronic message and that their expression is primarily regulated at the RNA level during latency in iSLK-219 cells. It remains possible that translational control governs the latent expression of ORF K1.

[Figure omitted, see PDF]

Figure 7. Gene expression is limited during latency in iSLK-219 cells.

(A) Log2 of normalized mRNA-Seq, Ribo-Seq counts (reads per kilobase million, rpkM) and translation efficiency (TE = Ribo-seq rpkM/mRNA-seq rpkM) for the latent transcripts and two late lytic ORFs (ORF25 and ORF62). Note the low TE for the latent genes at 0 h and 8 h. (B-C) Ribosome occupancy of latent transcripts changes with time. mRNA-seq (red) and Ribo-seq (CHX blue) for (B) the latency locus and (C) the K1-ORF4 region. (D-E) Immunoblot for the products of the latent transcripts (D) LANA and vCyclin, and (E) Kaposin. (F) Northern blot for the tricistronic (LANA-vCyclin-vFLIP) and the bicistronic (vCyclin-vFLIP) illustrated in (G). mRNA (100 ng-iSLK-219) or total RNA (10 ug-BCBL-1) were probed for vCyclin and vFLIP (see black bar in G). (G) Schematic of the tricistronic and bicistronic transcripts in the major latency locus.

[Figure omitted, see PDF]

Figure 8. Lytic reactivation is subject to a transcriptional timer.

(A-B) Log2 of normalized mRNA-Seq (A) and Ribo-Seq (B) counts (reads per kilobase million, rpkM) was calculated for the primary ORFs. The log2 scores reflect mRNA-expression (mRNA-seq) or ribosome occupancy (Ribo-seq) for each gene. The ORFs are organized by genomic position from ORFK1 to ORF75. The ORFs that were not annotated in our study have not been included in the plots (C) mRNA-Seq for the ORF58-62 locus. Three TSS sharing a single PAS were mapped in this region. Note the timing of expression for each of the transcripts.

As to the lytic cycle, our data support that most gene expression during this phase is controlled through transcriptional regulation, as the vast majority of the newly synthesized mRNAs are protected by ribosomes without delay (Figure 8A and 8B). A notable example of the temporal selection of distinct TSS can be seen in the transcription of the ORF58-62 locus which encodes the EBV-BMRF2 homologue (ORF58), the DNA polymerase processivity factor (ORF59), the small (ORF60) and large (ORF61) subunits of the viral ribonucleotide reductase, and a small capsid protein (ORF62) (Figure 8C) [61], [62]. In this case, our data clearly distinguish at least three independent transcripts with different expression kinetics: the RNA of ORF58-59 is expressed first, followed by the delayed early ORF60-61, and the late ORF62 transcript. The differential expression of these mRNAs correlates with their biological function, as genes required for DNA replication (ORF59, 60 and 61) are expressed before structural proteins (ORF62). Through the detection of these three distinct cistrons, our results unequivocally support the existence of three independent promoters that are integral part of the aforementioned temporally-regulated gene activation. Indeed, our TATA-box analyses, as well as previous reports, have identified at least three temporally-regulated promoters in this region (Table S1), supporting the differential expression control of ORF58-62. The time regulated selection of TSS has also been reported in HCMV, where the expression of transcripts from alternative TSS results in the translation of different protein products at specific times during infection, and represents a conserved mechanism of gene regulation in herpesviruses [14]. Taken together, our results, alongside previous studies, suggest that the transcription of KSHV genes is tightly regulated by promoter availability and the dynamic interplay of host and viral transcription factors in a time-dependent phase switch operating in the transition from latency into the lytic cycle [63], [64].

Discussion

The wealth of information we generated by combining our DNA-Seq, mRNA-Seq and Ribo-Seq datasets allowed us to build a comprehensive, high-resolution map of the viral genome over the KSHV life cycle. Our approach showcases the great analytical power of next-generation sequencing technologies, as we were able to pan-genomically map coding, non-coding, and regulatory features of the KSHV episome in iSLK-219 cells. Most importantly, this study provides tangible evidence derived from experimental data, as opposed to in silico prediction approaches, of the mechanisms employed by herpesviruses to widen the coding capacity of their genome through the use of diverse strategies including splicing, mRNA recoding, and alternative start codon usage. Furthermore, we demonstrate that the viral genome is not a conventional source of coding sequences as traditionally defined by in silico predictions and homology analyses, but a rather rich collection of diverse coding sequences that include numerous viral sORFs and uORFs. Some of those may exert translational modulation of other viral ORFs, thereby enriching and adding complexity to the viral gene regulatory profile. It is noteworthy that several dozen viral features uncovered by our group in this study were not documented earlier because of the limitations imposed by conventional methods employed to study gene expression.

Our results confirm the striking conservation of eukaryotic cis-regulatory elements in the KSHV genome, including TATA boxes, polyadenylation signals, and splice junctions. The conservation of both, the sequence and position of these features reflects the strict dependence of the virus on host factors and importantly, imposes a major constraint for the pathogen because of the need to compete with host factors for transcription and translation of virally-encoded products. To circumvent this competition, the virus causes the gradual but massive inhibition of host protein expression by increasing mRNA turnover. This mechanism, known as host shutoff, is orchestrated by the concerted action of the viral endonuclease SOX and the cellular exonuclease Xrn1 [65]. The elimination of competing cellular mRNAs renders the translational machinery components available to the virus, relieving the restrictions with respect to the expression of viral gene products. In addition to this upsurge in the availability of translation factors caused by the decreasing amounts of host transcripts, KSHV also hijacks essential cellular pathways that directly impinge on translation. Notable examples are the virus-orchestrated activation of mTOR and MAP kinases, which promote translation during the lytic cycle and result in the expression of viral ORFs required for the progression of infection [66]-[68]. It thus follows that the tug of war between host and virus would result in the commandeering of essential cellular factors by the virus to promote its own replication. Based on our analyses, we posit that the increased accessibility to ribosomes, the creation of an environment that is conducive to high translation, and the accumulation of viral transcripts results in an extremely favorable scenario for pervasive translation of viral encoded ORFs herein annotated, including sORFs and uORFs, and the short ORFs in PAN. Future studies in cells infected with viruses defective for host shutoff will help clarify the contribution of this mechanism to the viral translational output.

Another striking yet puzzling finding derived from our data is the protection of PAN by initiating and elongating ribosomes. This well-characterized and abundant viral RNA is expressed during the lytic cycle and has a predominant nuclear localization [33], [38]. Surprisingly, this nuclear RNA is clearly bound by ribosomes, indicating that either a fraction of PAN is cytoplasmic and available for translation, or (less likely) that ribosomes can access PAN in the nuclear compartment [69], [70]. Even more remarkable is the fact that the relative abundance of ribosome-bound PAN fragments suggests that the encoded peptides, if stable, could be abundant. Irrespective of their abundance, the biological functions of such putative peptides remain to be investigated. Bioinformatics analyses failed to identify discernible domains or particular motifs in 2 out of the 3 PAN-encoded peptides identified in our study. We did however find recognizable motifs within the primary structure of PAN1.1. This peptide harbors a putative signal sequence but lacks a discernible transmembrane domain or an ER retention signal, thus suggesting that PAN1.1 may be a secreted product. Studies aimed at the identification and characterization of PAN1.1 products are currently ongoing.

An interesting observation from our analyses is the detection of mRNA editing in two viral transcripts. This posttranscriptional mechanism of coding and non-coding RNA editing is conserved in eukaryotes. In humans extensive RNA editing has been reported, the majority of the events corresponding to A-to-I transitions mediated by the family of deaminases ADAR. Indeed, alterations in the activity of these enzymes are associated with disease [48], [71]. RNA editing is not restricted to eukaryotic messages and has been observed on viral transcripts in cells infected with RNA and DNA viruses [72]. The effects of RNA editing of viral transcripts antagonize (HCV, MV), or promote (HIV, HDV) viral activity and may affect the function of particular non-coding RNAs (EBV) or viral proteins (KSHV). As in humans, most of the RNA editing events detected on viral transcripts are ADAR dependent and correspond to A-to-I transitions, as the one we observe in Kaposin (position 117809). Interestingly, we see the protein levels of ADAR1 increasing throughout the viral lytic cycle, indicating the message for this enzyme partially escapes the widespread host shutoff caused by SOX. It is not clear however what would be the effects of such up-regulation of ADAR on RNA editing of host messages, given that a large number of transcripts are degraded during KSHV lytic infection. In addition to the nucleotide change in Kaposin, we also detect the non-canonical G-to-A editing in RTA (positions 72795 and 74801). While rare and less frequent than the A-to-I change, the G-to-A recoding has been previously identified in HIV and prostate and colorectal cancer [73], [74]. The enzyme responsible for this editing event has not been characterized to date, and KSHV lytic infection may provide a useful system for its identification.

Remarkably, our results strongly suggest that the peptide coding capacity of KSHV has been previously underestimated. This is illustrated in the overall high representation of sORFs and uORFs in the viral genome, often translated from near canonical start codons. The functional implications of such translation initiation events have been documented in several eukaryotes (yeast and metazoans) where translation initiation from near canonical codons under physiological conditions occurs more often than anticipated [21], [75]-[77]. In the case of KSHV, it is tempting to speculate that the increased use of non-canonical start codons, particularly late during infection (48 and 72 hr), is a probable consequence of the high availability of ribosomes and the translational permissiveness observed during the lytic cycle. A plausible explanation for the non-canonical start codon selection could be the abundance of eukaryotic translation initiation factors (eIFs) such as eIF1 and eIF5. It has been shown that these two proteins have opposite effects in the selection of the start codon; eIF1 increases the stringency of codon selection favoring initiation from AUG initiator codons whereas eIF5 favors translation from non-AUG codons and AUG codons nested within a poor sequence context [78], [79]. The relative abundance and stability of these factors has not been characterized in KSHV-infected cells. One could propose that the ratio of these factors may change during reactivation, favoring the translation from non-AUG codons. Future studies exploring the sequence context of the start codons of viral sORFs, uORFs, and alternative variants of main ORFs, as well as studies aimed at uncovering the interplay between eIF1 and eIF5 during infection will undoubtedly shed light on the mechanisms governing the intricate translation patterns we observed during the KSHV lytic cycle.

Another surprise revealed in our studies was the abundance of uORFs encoded in the 5' ends of viral transcripts. Like sORFs, uORFs have been found in eukaryotes, where they serve as modulators of gene expression during cellular stresses [80]. Their role has been extensively characterized in the face of amino acid starvation in yeast, where the accumulation of the transcription factor GCN4 depends on regulatory uORFs [81], [82]. This regulatory mechanism is evolutionarily conserved. In metazoans, uORFs regulate the abundance of the stress-responsive transcription factors C/EBPa/b, ATF4 and CHOP as well as the regulatory subunit of protein phosphatase 1, encoded by the GADD34 gene [58], [83], [84]. Although the role of uORFs in viruses remains largely unexplored, these regulatory elements permeate many viral families [14], [85]-[87] suggesting they may also control viral gene expression in instances where cellular stress pathways are engaged. In KSHV, the translational regulatory function of uORFs controlling the expression of ORF35 and ORF36, has been recently described [57]. The existence of a plethora of uORFs throughout the viral genome strongly indicates that this mechanism may be more widely used by KSHV than previously suspected.

Taken together, our results illustrate the dynamics of gene regulation during the different stages of the KSHV life cycle, while they also reveal that the coding capacity of its genome goes well beyond what was anticipated by in silico analyses. The striking conservation of the mechanisms employed by host and virus to regulate transcription, translation, and the diversity of the peptide repertoire, elevates KSHV as a valuable model system to explore the mechanistic underpinnings of the host-virus interactions in herpesviruses at large, as well as those of fundamental cellular processes, including the control of translation initiation in response to cellular stress.

Materials and Methods

Cell lines

iSLK and iSLK-219 cells (kindly provided by JinJong Myoung) were maintained in DMEM supplemented with 10% fetal bovine serum, L-glutamine (2 mM, Invitrogen), penicillin (100 IU/ml, Gibco) and streptomycin (100 ug/ml, Gibco) at 37°C under a 5% CO2 atmosphere. iSLK-219 cells were grown in the presence of puromycin (10 mg/ml, Invivogen) to maintain selection for the viral episome. LEC-219 cells were maintained in EBM-2 (Lonza cc3156) media supplemented with the EGM2-MV kit (Lonza cc3203) in presence of 0.25 ug/ml puromycin to maintain selection for the viral episome. BCBL-1 cells were maintained in RPMI supplemented with 10% fetal bovine serum, L-glutamine (2 mM, Invitrogen), penicillin (100 IU/ml, Gibco) and streptomycin (100 ug/ml, Gibco) at 37°C under a 5% CO2 atmosphere.

Induction and assessment of lytic replication in iSLK-219 cells

To induce entry to the lytic cycle, iSLK-219 cells were seeded at 1-2.5×105 cells/ml and 24 h after seeding (~70-80% confluent) cells were induced with doxycycline (1 µg/ml, BD Bioscience cat 631311,). To prevent viral DNA replication in the cells collected at 24 hpi, we induced these cells with Dox in the presence of phosphonoformate (500 uM). All other time points were treated with Dox alone. At the indicated times after induction viral reactivation was evaluated by microscopy detection of the PAN-RFP reporter. To determine the timing of KSHV DNA replication, DNA was isolated at the indicated times following reactivation using the DNeasy blood and tissue kit following manufacturer guidelines (Qiagen cat 69581). The DNA (20 ng) was used for qPCR using primers for the LANA promoter (Fwd: 5' AGGATGGAGATCGCAGACAC 3'; Rev: 5' CCAGCAAACCCACTTTAACC 3') or GAPDH as a normalization control (Fwd: 5' AGCCACATCGCTCAGACAC 3'; Rev: 5' TGGAAGATGGTGATGGGATT 3').

Immunoblotting and antibodies

Cells were washed and collected in RIPA buffer (10 mM Tris pH 8; 1 mM EDTA; 150 mM NaCl; 5% glycerol; 0.1% sodium deoxycholate; 0.1%SDS; 1%Triton x-100) at the indicated time points. Cell lysates were clarified and protein concentration was quantified using the Bio-Rad DC protein assay following manufacturer guidelines. 10 ug of protein per sample were fractionated by Bis-Tris PAGE (4-20% gels in MES SDS-running buffer) and transferred to nitrocellulose membranes. Immunoblots were processed, incubated with primary antibody overnight and developed using ECL reagent according to manufacturer recommendations.

The following antibodies were purchased from commercial suppliers: K8.1 (Advanced biotechnologies cat. 13-213-100) Actin (Sigma-Aldrich cat. A2228), ADAR1 (Abcam cat. ab126745), vCyclin (SCBT cat. 19415) and K-bZIP (SCBT cat. F33P1). The LANA rabbit polyclonal antibody was raised against a synthetic peptide from the acidic domain of LANA (Polson and Ganem, unpublished). The Kaposin B and C antibody was raised against the DR1/DR2 regions of the protein (Bisson and Ganem, unpublished). The ORF54 rabbit polyclonal antibody was raised against a synthetic peptide from the C-terminal region of the protein (aa 280-298 EDTNSVRKHTNEDNPVHEP) (Covance).

Preparation of ribosome profiling samples

Ribosome profiling was performed as previously described [13], [21] with some modifications. In brief, iSLK and iSLK-219 cells were left untreated, or treated with cycloheximide (100 ug/ml, 1 minute) or harringtonine (1 ug/ml, 90 seconds) followed by cycloheximide (100 ug/ml, 1 minute). After treatment cells were washed in cold PBS twice and lysed in lysis buffer (20 mM Tris, 1% triton, 220 mM NaCl, 15 mM MgCl2, 1 mM DTT, 8% glycerol) supplemented with Turbo DNAse (Ambion cat. AM2238). The lysate was clarified, treated with RNAse I (Ambion cat. 2295) and overlaid on a 34% sucrose cushion. Monosomes were isolated by centrifugation at 69000 rpm for 4 hr in a TLA 110 rotor. Ribosome protected RNA fragments were isolated from the monosomal fraction by acidic-phenol extraction.

Total RNA extraction and polyadenylated RNA selection for mRNA-Seq

Total RNA was extracted from 1×107 cells using RNABee (AMSBIO cat. Cs-104B) following the manufacturer guidelines. Polyadenylated RNA was isolated from the total fraction using Oligotex mRNA kit (Qiagen cat. 70022). The resulting mRNA was partially fragmented by alkaline hydrolysis with sodium carbonate to ~150 nt segments on average, and then fragments between 40-100 nt were isolated from gel. The laboratory of Jonathan Weissman has previously documented that this partial fragmentation results in the preferential accumulation of mRNA 5' terminal fragments for most transcripts with a non-overlapping transcription start site [14]. We determined that for the transcripts where we were able to annotate a transcription start site, there is a 4-5 fold enrichment (average 5-fold, median 4-fold) of the number of reads for the first 10 nt of the transcript compared to three 10 nt windows within the gene body (nt 20-30, 30-40, 50-60 data not shown). The mRNA profiles in all figures show counts of the 5'-most bases of sequencing reads.

Library generation, sequencing, sequence alignments, normalization, and splice junction discovery

Strand specific libraries were generated as in Ingolia et al., 2012, with the modifications described in Stern-Ginossar et al., 2012. Samples were sequenced on either the Illumina Genome Analyzer II or the HiSeq 2000 using the truseq sbs kit v3 50 cycles (Illumina cat. FC-401-3002). Sequence analysis was done as described in Stern-Ginossar et al., 2012. Briefly, linker and polyA sequences were removed from the 3' end of the reads preceding the alignment. Sequencing reads were aligned with Bowtie2 [88] allowing for 2 mismatches. Sequences aligning to rRNA were discarded and the remaining reads were aligned to KSHV (GQ994935.1) and human (hg19) genomes. For normalization, uniquely mapped reads were used to calculate the mRNA and ribosome footprint reads per kilobase per million (rpkM), and regions containing multi-mapped reads were masked out. The rpkM/gene for mRNA-Seq and Ribo-Seq were visualized using TreeView (1.1.6).

The sequences that did not align to the viral and host genomes were analyzed with TopHat and HMMSplicer for splice junction discovery using default options [22], [23]. We annotated splice junctions present in at least two time points, with an HMMSplicer score>900 and a TopHat score>7. We determined these thresholds based on previously characterized splice junctions. Two previously reported and two novel splice junctions had low or no TopHat score, but were included in the annotations due to their detection in multiple samples and high HMMSplicer score (ORFK1/ORF4, ORF70, ORF46/47, and Kaposin).

While other strong and reliable bioinformatics tools, such as Cufflinks and Scripture, are available for transcript reconstruction, we were not able to use them for viral mRNAs annotation. These tools are optimized for the analysis of transcripts in genomes where transcriptional units are well spaced and well defined. In the case of the compact KSHV genome, most transcripts use common regulatory features, are overlapping, or are very close to each other. Furthermore, and as shown in our study and previous studies by our lab and others, virtually the entire viral genome is transcribed late during the lytic cycle [6], [9]. These conditions make it difficult for these bioinformatics tools to parse out, identify and predict viral transcripts.

Prediction of translation initiation sites

The identification of translation initiation sites was done using a machine learning approach as previously described Stern-Ginossar et al., 2012. In total our approach successfully predicted the translation initiation site of 64% (56 of 87) of the previously annotated ORFs (Tables 1, 2 and 3). The sites that were not predicted correspond to regions of low read coverage or overlapping ORFs.

Annotation of 3' termini

Stretches of 5 or more consecutive adenosines (polyA), allowing one non-A base for every 5, were removed from the 3' ends of mRNA sequencing reads before alignment. These reads were aligned to the viral genome and the 3' end was determined as the last nucleotide before the start of the polyA stretch. To prevent false positives, polyA sequences were only used as evidence of 3' termini when they mismatched the underlying reference genome sequence.

Validation of ORF57 splice junction

For the validation of the second splicing event within the 3' end of the ORF57 transcript we prepared cDNA from 1 ug of total RNA using the qScript cDNA-Supermix (Quanta cat. 95048-025) or SuperScript III First Strand Synthesis System for RT-PCR (Invitrogen cat. 18080-051) with a combination of oligo(dT) and random hexamers, following manufacturers recommendations. PCR was done using as a template 1% of the resulting cDNA. The following primers flanking the 5' and 3' ends of the second splice junction were used for 30 cycles of amplification (Fwd: 5' GGCAAAGACGACGAACTCAT 3' Rev: 5' GAGAAGAGACCACGCCTGACT 3'). The resulting products were separated in a 1.2% agarose gel and stained with ethidium bromide for visualization.

For the validation of the 3' end extension of ORF57, PCR was done as described above using the following primers (A-Fwd: 5' GGGTGGTTTGATGAGAAGGA, B-Fwd: 5' TGGCAGAGTGTCTCCCGTAT, C-Rev: 5' GAGAAGAGACCACGCCTGACT, D-Rev: 5' ATAATGCCGAAGCCGTTATG)

Validation of mRNA editing by Sanger sequencing

Total RNA was extracted from cells at different times following induction of reactivation by doxycycline treatment, and cDNA was prepared as described above. Gene specific amplification was done using Phusion High-Fidelity DNA Polymerase (NEB cat. M0530S) for 32 cycles following manufacturer guidelines. The following primers were used for amplification: RTA-74281: Fwd-T7 5' taatacgactcactatagggACGCGCTGTTGTCCAGTATT 3', Rev-T3 5' aattaaccctcactaaagggGTACAGTGTGCCGGACTCCT 3'; RTA-72795: Fwd-T7 5' taatacgactcactatagggCCTCTCGAATGAGGACCAAA 3', Rev-T3 5' aattaaccctcactaaagggGTAGACCGGTTGGAAAACCA 3'; Kaposin-117809: Fwd-T7 5' taatacgactcactatagggGTTGCAACTCGTGTCCTGAA 3', Rev-T3 5' aattaaccctcactaaagggAGGCTTAACGGTGTTTGTGG 3' ; ORF6-6144: Fwd-T7 5' taatacgactcactatagggGGGATACTTCTCGGGGAGAG 3', Rev-T3 5' aattaaccctcactaaagggGGCCCACTGTGCTCAGTAAT 3'; ORF63-102377: Fwd-T7 5' taatacgactcactatagggGTTGGAAAATATCGCGTGCT 3', Rev-T3 5' aattaaccctcactaaagggTTGTGTGTTCGGTCCTGTGT 3'; ORF59-96434 Fwd-T7 5' taatacgactcactatagggGGACGTGACCCTCCTGTCTA 3', Rev-T3 5' aattaaccctcactaaagggTAACGTCTCCACTGCCTTCC3'.

Northern blotting

Total RNA was extracted from cells using RNABee (AMSBIO cat. Cs-104B) following manufacturer instructions. Northern blotting was done for 10 ug of RNA, or 100 ng of mRNA per lane using the Ambion NorthernMax system (Invitrogen cat. AM1940). Gene and sense specific riboprobes were synthesized from PCR products using the Ambion MAXIscript T7-T3 Transcription Kit (Ambion cat. AM1326) according to manufacturer guidelines. The primers used for the PCR amplification of gene specific probes are: Probe A (24151-25437): Fwd-T7 5' taatacgactcactatagggagaCAGTCACAAGCACACAACCC3', Rev-T3 5' aattaaccctcactaaagggagaTTCGGGTGATTAAGCAAAGG3'; Probe B (25437-25929) : Fwd-T7 taatacgactcactatagggagaCCTTTGCTTAATCACCCGAA, Rev-T3 5' aattaaccctcactaaagggagaGGTGACCGTACTGCCATACC3'; Probe C (123441-122984): Fwd-T7 5' taatacgactcactatagggCGCTAACAGGGGAAACGTTAACCTGC3', Rev-T3 5' aattaaccctcactaaagggCTCATTGCCCGCCTCTATTA 3'

Mochiview files

File S2 contains the updated KSHV 2.0 annotations, mRNA profiles, and ribosome occupancy (CHX) plots. The database can be opened using the Mochiview file software [26], free for download at http://johnsonlab.ucsf.edu/sj/mochiview-software/

To open the file, import and activate the database. Restart the program and the new database will be available for viewing. The following files are included: GQ994935 sequence, KSHV2.0 location, mRNA_dox8h_minus, mRNA_dox8h_plus, mRNA_dox24h_minus, mRNA_dox24h_plus, mRNA_dox48h_minus, mRNA_dox48h_plus, mRNA_dox72h_minus, mRNA_dox72h_plus, fp_dox8h_minus, fp_dox8h_plus, fp_dox24h_minus, fp_dox24h_plus, fp_dox48h_minus, fp_dox48h_plus, fp_dox72h_minus, fp_dox72h_plus, fp_harr_dox48h_minus, fp_harr_dox48h_plus, fp_harr_dox72h_minus, fp_harr_dox72h_plus

Supporting Information

Figure S1.

Lytic reactivation of KSHV is induced by exogenous expression of RTA in iSLK-219 cells. (A) Immunoblot of latent (LANA) and late lytic (K8.1) products in iSLK-219 cells (B) Quantitative PCR of viral DNA, shows DNA replication starting at 48 hpi and increasing with time.

(TIFF)

Figure S2.

Read length distribution and sequencing-coverage for mRNA-seq and Ribo-seq in iSLK-219. (A) Length of fragmented, size selected mRNA (40-100 nt) and monosome-protected footprints (~30 nt) in lysates from cells in the lytic cycle (72 hr). (B) Number of reads for mRNA-seq and Ribo-seq of all samples included in this study. Note that the number of mRNA-seq and Ribo-seq reads that align to the viral cycle increase as the lytic cycle progresses.

(TIF)

Figure S3.

Transcription of viral genes is highly permissive during the lytic cycle. mRNA-seq (red) and Ribo-seq (CHX blue) profiles for the entire KSHV genome at 8 and 72 hr following Dox induction. Note the change in the read density as the viral lytic cycle progresses. The Y axis (number of reads) was cut at 50 (plus strand) and -50 (minus strand) for ease of visualization.

(TIFF)

Figure S4.

ALT and K1-11-AS are modestly bound by ribosomes. (A and B) mRNA-seq and Ribo-seq (CHX and Harringtonine) profiles for the lincRNAs K1/11-Antisense (A) and ALT (B) at 72 hr post lytic reactivation. Only the strand of RNA coding for the lincRNA is shown. Solid blue arrows represent transcripts, light blue arrow heads coding regions and red thick arrows represent the lincRNA.

(TIFF)

Figure S5.

The highly abundant viral transcript PAN, codes for three putative small peptides. (A) Timecourse of ribosome accumulation on the putative peptides within PAN. Notice the difference in the scale of the number of reads (y axis). The double head arrow indicates that the PAN transcript continues after the region shown in this figure (B) Accumulation of releasing ribosomes at the stop codon of PAN1.1 (top panel) and K8.1 (bottom panel). Where indicated the mRNA-seq and Ribo-seq (CHX and no drug) profiles for PAN (8 hr) and K8.1 (48 hr) coding regions are shown. (C) mRNA-seq and Ribo-seq of the ORFK7 and putative PAN peptides at 72 hr post reactivation. The right panel is a zoom in of the start codon of ORFK7. Notice the difference in the number of reads on the y axis. (D and E) Ribosome Release Score (RRS) was calculated as RSS = [(footprint reads coding region/footprint reads 3'UTR)/(mRNA reads coding region/mRNA reads 3'UTR)] using the read values from the 48 and 72 hr timepoints. RRS above 14 have been previously calculated for coding transcripts by Guttman et al., thus strongly supporting the coding capacity of PAN [40].

(TIFF)

Figure S6.

Identification of a novel splice junction within the 3' end of ORF57. (A) Presence of a novel splice site within the 3' of ORF57. The ORF57 sequence from 81886 to 83490 (GQ994935) includes the previously annotated 5' splice and the novel 3' splice annotated in this study. The splice donor and acceptor sites are highlighted in italic bold. Underlined is the sequence of the 3' and 5' introns. (B) Splicing of the second intron of ORF57 was evaluated by end-point PCR on cDNA from latent and lytic BCBL-1 and iSLK-219 cells, using primers flanking the intron boundaries (black solid lines). (C) The presence of a 3' extension of the ORF57 transcript was confirmed in lytic iSLK-219 and BCBL-1 cells by end-point PCR, using three combinations of primers (A-forward, B-forward, C-reverse and D-reverse) flanking the annotated polyA cleavage site for this transcript. 28S was used as an internal loading control.

(TIF)

Figure S7.

Higher levels of Kaposin mRNA-editing correlate with an increase in ADAR during the lytic cycle. (A) PCR amplification followed by Sanger sequencing of the region of Kaposin cDNA where mRNA editing occurs. Highlighted by the red box is the edited position, nucleotide 117809. Note the gradual and dramatic change from A to G as the lytic cycle progresses. (B) Quantification of the results obtained in (A). (C) Changes in the levels of ADAR isoforms (p150, p110 and p80) during the lytic cycle. Actin: Loading control, KbZip: Viral reactivation control, NS: not specific. (D) Fold changes in the mRNA and ribosome occupancy levels for two cellular genes that escape the host shutoff response, interleukin 6 (IL6) and the Hypoxia inducible factor 1 (HIF1), one gene sensitive to host shutoff protein kinase, DNA-activated, catalytic polypeptide (PRKDC) and ADAR. The fold change was calculated as the Log2 of the ratio of rpkM at 48 hr over 0 hr, for the mRNA-seq (left panel) or Ribo-Seq (right panel) samples. The black and grey bars represent two independent biological replicates. Note that the changes in ADAR mRNA and ribosome occupancy are very similar to the host shutoff escapee HIF1.

(TIFF)

Figure S8.

KSHV uses alternative start codons to increase peptide diversity. mRNA-Seq and Ribo-Seq (CHX, Harringtonine and No drug) profiles for ORF K6 (A) and ORF62 at 72 h post reactivation (B). The blue arrowheads represent the annotated ORF and the red arrowheads represent the ORFs expressed from the alternative start codons. The dashed lines mark the multiple translation initiation start sites.

(TIF)

Figure S9.

uORFs and sORFs populate the viral genome in KSHV. (A) Percentage of uORFs and sORFs starting from canonical and near cognate start codons. Other codons are ATT, ATC, ATA, TTG and ACG (see Table S7). (B) The ORF50-AS transcript is protected by ribosomes. mRNA-seq and Ribo-seq profiles (CHX and Harringtonine) of the ORF50-AS/ORF49 region. The blue and grey arrows represent transcripts, the blue and grey arrowheads represent the coding regions of ORF49 and ORF50 and the yellow solid line marks the region of ORF50-AS where ribosomes bind. (C) Numerous sORFs are found within the K4.1/K4.2 transcript. mRNA-seq and Ribo-seq (CHX and Harringtonine) for the K4/K4.1/K4.2 region. The blue arrow and arrowhead represent the K4 transcript, coding region and uORF. The red arrow and arrowheads represent the K4.1/K4.2 transcript and the numerous uORFs encoded by this transcript. The grey arrow and arrowhead represents the K4.1 and K4.2 genes previously annotated, not detected in this study.

(TIFF)

Table S1.

Transcription start sites in KSHV. Annotations for the transcription start sites described in this study. The second column contains the sequence of 50 nt 5' and 10 nt 3' from the annotated TSS. The TATA box is in bold, italic font. The third column indicates the distance from the start of the TATA box to the TSS. The eight TATA-less genes are indicated as No TATA.

(DOCX)

Table S2.

PolyA sites in KSHV. Polyadenylation sites were annotated in the KSHV genome using mRNA-seq. Cleavage sites were determined as the last nucleotide before a stretch of 5 or more consecutive adenosines. The sequence in column 7 corresponds to 30 nt flanking the polyA cleavage site. The AA/UUAA and GU sites are in italic, bold font and the cleavage site(s) are in bold font.

(DOCX)

Table S3.

PAN represent a large percentage of the KSHV mRNA-seq. The PAN mRNA or ribosome footprint (FP) reads were calculated for the region between 28661 to 29690. The percentage of KSHV reads was calculated using the total number of viral reads for mRNA-seq or Ribo-seq for each timepoint. 1. The PAN RNA reads in this latent sample are likely coming from the 0.1% of spontaneous lytic cells or from minor cross-contamination during library preparation. 2. Technical replicates.

(DOCX)

Table S4.

Splice junctions in the KSHV transcriptome. Splice junctions were identified and annotated using TopHat and HMMsplicer. The novel splice junctions identified in this study are in bold font. Donor and acceptor sites are in italic, bold font (column 6). Columns 7 and 8 contain the TopHat and HMM generated scores for two replicates at 48 hr (rep1 and rep2) and one at 72 hr post reactivation. These scores reflect the strength of the alignment (HMM splicer) and the number of reads that mapped to the region (TopHat and HMM splicer).

(DOCX)

Table S5.

mRNA editing in KSHV. mRNA editing events were predicted in silico and were annotated by position and gene. We include in this table a DNA mutation identified in the GQ994935.1 genome in ORF63 at position 102377 (A-to-C change). Three editing events, Kaposin 117809, RTA 72841 and RTA 72795 have been confirmed by Sanger sequencing.

(DOCX)

Table S6.

Alternative start codon usage and internal ORFs in KSHV. From our Ribo-seq data we identified five viral ORFs (ORF70, K6, ORF54, ORF62 and Kaposin) with multiple in frame initiation codons where ribosomes accumulate. The size of the predicted protein products originating from the alternative start codons is indicated in the column labeled amino acids. In addition, five viral ORFs (ORF10, ORF11, K3, ORF20 and K8.1) contain internal initiation codons, in or out of frame with the primary ORF, where we detect initiating ribosomes. The size of the predicted products is indicated in the column labeled amino acids (Nuc.: nucleotides, M: Manual annotation, P: Predicted in silico).

(DOCX)

Table S7.

Upstream and small ORFs are widely distributed in KSHV. The translation boundaries for uORFs and sORFs were predicted by the SVM based on the Ribo-seq data from harringtonine treated cells at 72 h post reactivation. The ORF coordinates were then manually curated. The translation efficiency (TE) for the ORFs larger than 15 amino acids was determined as Footprint-rpkM/mRNA-rpkM. For the calculation we excluded the first 45 nucleotides from the start codon. (M: Manual annotation, P: Predicted in silico, N/A: not available, *: Start codons are ambiguous).

(DOCX)

File S1.

Summary of novel features annotated for KSHV. The novel genomic features annotated in our study include the ORFs derived from alternative start codons, internal, small and upstream ORFs, and polyadenylation sites. All coordinates are annotated based on the human herpesvirus 8 accession number GQ994935.1.

(XLSX)

File S2.

Mochiview files for updated annotation and mRNA-seq/Ribo-seq profiles. Supplementary file 1 is a Mochiview file [26] containing the updated KSHV 2.0 annotations, mRNA profiles, and ribosome occupancy (CHX) plots for all of the time points analyzed (0, 8, 24, 48, 72 hr). Ribosome occupancy plots for Harr-treated cells are included for the 48 and 72 hr time points. See materials and methods for uploading instructions.

(ZIP)

Acknowledgments

We thank JinJong Myoung, Calvin Jan, Joshua Dunn, Eugene Oh, Clement Chu, Amy Kistler and Sanjay Chandriani for providing reagents and technical assistance. We thank Matthew Holderfield, Claudia Rubio, Lisa Kronstad, Britt Glaunsinger, Joe DeRisi and the members of the DeRisi and Ganem labs for insightful scientific discussions. We thank Diego Acosta-Alvear for valuable technical advice and critical reading of the manuscript.

Author Contributions

Conceived and designed the experiments: CA JSW DG. Performed the experiments: CA NSG ASM MH PB. Analyzed the data: CA BW AM. Contributed reagents/materials/analysis tools: AM ASM NSG. Wrote the paper: CA DG.

References

Cesarman E, Knowles DM (1999) The role of Kaposi's sarcoma-associated herpesvirus (KSHV/HHV-8) in lymphoproliferative diseases. Seminars in cancer biology 9: 165-174. doi: 10.1006/scbi.1998.0118

Ganem D (2010) KSHV and the pathogenesis of Kaposi sarcoma: listening to human biology and medicine. The Journal of clinical investigation 120: 939-949. doi: 10.1172/jci40567

Chang Y, Cesarman E, Pessin MS, Lee F, Culpepper J, et al. (1994) Identification of herpesvirus-like DNA sequences in AIDS-associated Kaposi's sarcoma. Science (New York, NY) 1865-1869. doi: 10.1126/science.7997879

Russo JJ (1996) Nucleotide sequence of the Kaposi sarcoma-associated herpesvirus (HHV8). Proceedings of the National Academy of Sciences 93: 14862-14867. doi: 10.1073/pnas.93.25.14862

Gottwein E (2012) Kaposi's Sarcoma-Associated Herpesvirus microRNAs. Frontiers in microbiology 3: 165. doi: 10.3389/fmicb.2012.00165

Chandriani S, Xu Y, Ganem D (2010) The lytic transcriptome of Kaposi's sarcoma-associated herpesvirus reveals extensive transcription of noncoding regions, including regions antisense to important genes. Journal of virology 84: 7934-7942. doi: 10.1128/jvi.00645-10

Xu Y, Ganem D (2010) Making sense of antisense: seemingly noncoding RNAs antisense to the master regulator of Kaposi's sarcoma-associated herpesvirus lytic replication do not regulate that transcript but serve as mRNAs encoding small peptides. Journal of virology 84: 5465-5475. doi: 10.1128/jvi.02705-09

Jaber T, Yuan Y (2013) A virally encoded small peptide regulates RTA stability and facilitates Kaposi's sarcoma-associated herpesvirus lytic replication. Journal of virology 87: 3461-3470. doi: 10.1128/jvi.02746-12

Sarid R, Flore O, Bohenzky RA, Chang Y, Moore PS (1998) 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 72: 1005-1012. doi: 10.1097/00042560-199804010-00059

Dittmer DP (2003) Transcription Profile of Kaposi's Sarcoma-associated Herpesvirus in Primary Kaposi's Sarcoma Lesions as Determined by Real-Time PCR Arrays. Cancer Res 63: 2010-2015

Dresang LR, Teuton JR, Feng H, Jacobs JM, Camp DG, et al. (2011) Coupled transcriptome and proteome analysis of human lymphotropic tumor viruses: insights on the detection and discovery of viral genes. BMC genomics 12: 625. doi: 10.1186/1471-2164-12-625

Ingolia NT, Ghaemmaghami S, Newman JRS, Weissman JS (2009) Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling. Science (New York, NY) 324: 218-223. doi: 10.1126/science.1168978

Ingolia NT, Brar GA, Rouskin S, McGeachy AM, Weissman JS (2012) The ribosome profiling strategy for monitoring translation in vivo by deep sequencing of ribosome-protected mRNA fragments. Nature protocols 7: 1534-1550. doi: 10.1038/nprot.2012.086

Stern-Ginossar N, Weisburd B, Michalski A, Le VTK, Hein MY, et al. (2012) Decoding human cytomegalovirus. Science (New York, NY) 338: 1088-1093. doi: 10.1126/science.1227919

Jenner RG, Albà MM, Boshoff C, Kellam P (2001) Kaposi's sarcoma-associated herpesvirus latent and lytic gene expression as revealed by DNA arrays. Journal of virology 75: 891-902. doi: 10.1128/jvi.75.2.891-902.2001

Paulose-Murphy M, Ha NK, Xiang C, Chen Y, Gillim L, et al. (2001) Transcription program of human herpesvirus 8 (kaposi's sarcoma-associated herpesvirus). Journal of virology 75: 4843-4853. doi: 10.1128/jvi.75.10.4843-4853.2001

Myoung J, Ganem D (2011) Generation of a doxycycline-inducible KSHV producer cell line of endothelial origin: maintenance of tight latency with efficient reactivation upon induction. Journal of virological methods 174: 12-21. doi: 10.1016/j.jviromet.2011.03.012

Vieira J, O'Hearn PM (2004) Use of the red fluorescent protein as a marker of Kaposi's sarcoma-associated herpesvirus lytic gene expression. Virology 325: 225-240. doi: 10.1016/j.virol.2004.03.049

Schneider-Poetsch T, Ju J, Eyler DE, Dang Y, Bhat S, et al. (2010) Inhibition of eukaryotic translation elongation by cycloheximide and lactimidomycin. Nature chemical biology 6: 209-217. doi: 10.1038/nchembio.304

Fresno M, Jimenez A, Vazquez D (1977) Inhibition of Translation in Eukaryotic Systems by Harringtonine. European Journal of Biochemistry 72: 323-330. doi: 10.1111/j.1432-1033.1977.tb11256.x

Ingolia NT, Lareau LF, Weissman JS (2011) Ribosome Profiling of Mouse Embryonic Stem Cells Reveals the Complexity and Dynamics of Mammalian Proteomes. Cell 147 (4) 789-802. doi: 10.1016/j.cell.2011.10.002

Trapnell C, Pachter L, Salzberg SL (2009) TopHat: discovering splice junctions with RNA-Seq. Bioinformatics (Oxford, England) 25: 1105-1111. doi: 10.1093/bioinformatics/btp120

Dimon MT, Sorber K, DeRisi JL (2010) HMMSplicer: a tool for efficient and sensitive discovery of known and novel splice junctions in RNA-Seq data. PloS one 5: e13875. doi: 10.1371/journal.pone.0013875

Tang S, Zheng Z-M (2002) Kaposi's sarcoma-associated herpesvirus K8 exon 3 contains three 5'-splice sites and harbors a K8.1 transcription start site. The Journal of biological chemistry 277: 14547-14556. doi: 10.1074/jbc.m111308200

Chandran B, Bloomer C, Chan SR, Zhu L, Goldstein E, et al. (1998) Human herpesvirus-8 ORF K8.1 gene encodes immunogenic glycoproteins generated by spliced transcripts. Virology 249: 140-149. doi: 10.1006/viro.1998.9316

Homann OR, Johnson AD (2010) MochiView: versatile software for genome browsing and DNA motif analysis. BMC biology 8: 49. doi: 10.1186/1741-7007-8-49

Zhu FX, Cusano T, Yuan Y (1999) Identification of the Immediate-Early Transcripts of Kaposi's Sarcoma-Associated Herpesvirus. Journal of virology 73: 5556-5567

Taylor JL, Bennett HN, Snyder BA, Moore PS, Chang Y (2005) Transcriptional analysis of latent and inducible Kaposi's sarcoma-associated herpesvirus transcripts in the K4 to K7 region. Journal of virology 79: 15099-15106. doi: 10.1128/jvi.79.24.15099-15106.2005

Sandelin A, Carninci P, Lenhard B, Ponjavic J, Hayashizaki Y, et al. (2007) Mammalian RNA polymerase II core promoters: insights from genome-wide studies. Nature reviews Genetics 8: 424-436. doi: 10.1038/nrg2026

Carninci P, Sandelin A, Lenhard B, Katayama S, Shimokawa K, et al. (2006) Genome-wide analysis of mammalian promoter architecture and evolution. Nature genetics 38: 626-635. doi: 10.1038/ng1789

Smale ST, Kadonaga JT (2003) The RNA polymerase II core promoter. Annual review of biochemistry 72: 449-479. doi: 10.1146/annurev.biochem.72.121801.161520

Bartkoski M, Roizman B (1976) RNA synthesis in cells infected with herpes simplex virus. XIII. Differences in the methylation patterns of viral RNA during the reproductive cycle. Journal of Virology 20: 583-588

Sun R, Lin S-F, Gradoville L, Miller G (1996) Polyadenylated nuclear RNA encoded by Kaposi sarcoma-associated herpesvirus. PNAS 93: 11883-11888. doi: 10.1073/pnas.93.21.11883

McClure LV, Kincaid RP, Burke JM, Grundhoff A, Sullivan CS (2013) Comprehensive Mapping and Analysis of Kaposi's Sarcoma-Associated Herpesvirus 3' UTRs Identify Differential Posttranscriptional Control of Gene Expression in Lytic versus Latent Infection. Journal of virology 87: 12838-12849. doi: 10.1128/jvi.02374-13

Bai Z, Huang Y, Li W, Zhu Y, Jung JU, et al. (2013) Genome-wide mapping and screening of KSHV 3' UTRs identify bicistronic and polycistronic viral transcripts as frequent targets of KSHV microRNAs. Journal of virology doi: 10.1128/jvi.02689-13

Majerciak V, Ni T, Yang W, Meng B, Zhu J, et al. (2013) A Viral Genome Landscape of RNA Polyadenylation from KSHV Latent to Lytic Infection. PLoS pathogens 9: e1003749. doi: 10.1371/journal.ppat.1003749

Proudfoot NJ, Brownlee GG (1976) 3' Non-coding region sequences in eukaryotic messenger RNA. Nature 263: 211-214. doi: 10.1038/263211a0

Zhong W, Ganem D (1997) Characterization of ribonucleoprotein complexes containing an abundant polyadenylated nuclear RNA encoded by Kaposi's sarcoma-associated herpesvirus (human herpesvirus 8). J Virol 71: 1207-1212

Rossetto CC, Pari G (2012) KSHV PAN RNA associates with demethylases UTX and JMJD3 to activate lytic replication through a physical interaction with the virus genome. PLoS pathogens 8: e1002680. doi: 10.1371/journal.ppat.1002680

Guttman M, Russell P, Ingolia NT, Weissman JS, Lander ES (2013) Ribosome Profiling Provides Evidence that Large Noncoding RNAs Do Not Encode Proteins. Cell 154: 240-251. doi: 10.1016/j.cell.2013.06.009

Wilson BA, Masel J (2011) Putatively Noncoding Transcripts Show Extensive Association with Ribosomes. Genome biology and evolution 3: 1245-1252. doi: 10.1093/gbe/evr099

Petersen TN, Brunak S, von Heijne G, Nielsen H (2011) SignalP 4.0: discriminating signal peptides from transmembrane regions. Nature methods 8: 785-786. doi: 10.1038/nmeth.1701

Hiller K, Grote A, Scheer M, Münch R, Jahn D (2004) PrediSi: prediction of signal peptides and their cleavage positions. Nucleic Acids Research 32: 375-379. doi: 10.1093/nar/gkh378

Conrad NK (2009) Posttranscriptional gene regulation in Kaposi's sarcoma-associated herpesvirus. Advances in applied microbiology 68: 241-261. doi: 10.1016/s0065-2164(09)01206-4

Jackson BR, Noerenberg M, Whitehouse A (2012) The Kaposi's Sarcoma-Associated Herpesvirus ORF57 Protein and Its Multiple Roles in mRNA Biogenesis. Frontiers in microbiology 3: 59. doi: 10.3389/fmicb.2012.00059

Majerciak V, Zheng Z-M (2009) Kaposi's sarcoma-associated herpesvirus ORF57 in viral RNA processing. Frontiers in Bioscience 14: 1516-1528. doi: 10.2741/3322

Mallela A, Nishikura K (2012) A-to-I editing of protein coding and noncoding RNAs. Critical reviews in biochemistry and molecular biology 47: 493-501. doi: 10.3109/10409238.2012.714350

Nishikura K (2010) Functions and regulation of RNA editing by ADAR deaminases. Annual review of biochemistry 79: 321-349. doi: 10.1146/annurev-biochem-060208-105251

Gandy SZ, Linnstaedt SD, Muralidhar S, Cashman KA, Rosenthal LJ, et al. (2007) RNA editing of the human herpesvirus 8 kaposin transcript eliminates its transforming activity and is induced during lytic replication. Journal of virology 81: 13544-13551. doi: 10.1128/jvi.01521-07

Glaunsinger B, Ganem D (2004) Lytic KSHV Infection Inhibits Host Gene Expression by Accelerating Global mRNA Turnover. Molecular Cell 13: 713-723. doi: 10.1016/s1097-2765(04)00091-7

Blom N, Sicheritz-Pontén T, Gupta R, Gammeltoft S, Brunak S (2004) Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence. Proteomics 4: 1633-1649. doi: 10.1002/pmic.200300771

Tsai W-H, Wang P-W, Lin S-Y, Wu I-L, Ko Y-C, et al. (2012) Ser-634 and Ser-636 of Kaposi's Sarcoma-Associated Herpesvirus RTA are Involved in Transactivation and are Potential Cdk9 Phosphorylation Sites. Frontiers in microbiology 3: 60. doi: 10.3389/fmicb.2012.00060

Madrid AS, Ganem D (2012) Kaposi's sarcoma-associated herpesvirus ORF54/dUTPase downregulates a ligand for the NK activating receptor NKp44. Journal of virology 86: 8693-8704. doi: 10.1128/jvi.00252-12

Kozak M (1981) Possible role of flanking nucleotides in recognition of the AUG initiator codon by eukaryotic ribosomes. Nucleic acids research 5233-5252. doi: 10.1093/nar/9.20.5233

Basrai MA, Hieter P, Boeke JD (1997) Small Open Reading Frames: Beautiful Needles in the Haystack. Genome Res 7: 768-771

Calvo SE, Pagliarini DJ, Mootha VK (2009) Upstream open reading frames cause widespread reduction of protein expression and are polymorphic among humans. Proceedings of the National Academy of Sciences of the United States of America 106: 7507-7512. doi: 10.1073/pnas.0810916106

Kronstad LM, Brulois KF, Jung JU, Glaunsinger BA (2013) Dual short upstream open reading frames control translation of a herpesviral polycistronic mRNA. PLoS pathogens 9: e1003156. doi: 10.1371/journal.ppat.1003156

Vattem KM, Wek RC (2004) Reinitiation involving upstream ORFs regulates ATF4 mRNA translation in mammalian cells. Proceedings of the National Academy of Sciences of the United States of America 101: 11269-11274. doi: 10.1073/pnas.0400541101

Dittmer D, Lagunoff M, Renne R, Staskus K, Haase A, et al. (1998) A Cluster of Latently Expressed Genes in Kaposi's Sarcoma-Associated Herpesvirus. J Virol 72: 8309-8315

Grundhoff A, Ganem D (2001) Mechanisms governing expression of the v-FLIP gene of Kaposi's sarcoma-associated herpesvirus. Journal of virology 75: 1857-1863. doi: 10.1128/jvi.75.4.1857-1863.2001

Majerciak V, Yamanegi K, Zheng Z-M (2006) Gene structure and expression of Kaposi's sarcoma-associated herpesvirus ORF56, ORF57, ORF58, and ORF59. Journal of virology 80: 11968-11981. doi: 10.1128/jvi.01394-06

Chan SR, Chandran B (2000) Characterization of Human Herpesvirus 8 ORF59 Protein (PF-8) and Mapping of the Processivity and Viral DNA Polymerase-Interacting Domains. Journal of Virology 74: 10920-10929. doi: 10.1128/jvi.74.23.10920-10929.2000

Günther T, Grundhoff A (2010) The epigenetic landscape of latent Kaposi sarcoma-associated herpesvirus genomes. PLoS pathogens 6: e1000935. doi: 10.1371/journal.ppat.1000935

Toth Z, Maglinte DT, Lee SH, Lee H-R, Wong L-Y, et al. (2010) Epigenetic analysis of KSHV latent and lytic genomes. PLoS pathogens 6: e1001013. doi: 10.1371/journal.ppat.1001013

Covarrubias S, Gaglia MM, Kumar GR, Wong W, Jackson AO, et al. (2011) Coordinated destruction of cellular messages in translation complexes by the gammaherpesvirus host shutoff factor and the mammalian exonuclease Xrn1. PLoS pathogens 7: e1002339. doi: 10.1371/journal.ppat.1002339

Arias C, Walsh D, Harbell J, Wilson AC, Mohr I (2009) Activation of host translational control pathways by a viral developmental switch. PLoS pathogens 5: e1000334. doi: 10.1371/journal.ppat.1000334

Cohen A, Brodie C, Sarid R (2006) An essential role of ERK signalling in TPA-induced reactivation of Kaposi's sarcoma-associated herpesvirus. The Journal of general virology 87: 795-802. doi: 10.1099/vir.0.81619-0

Kuang E, Wu F, Zhu F (2009) Mechanism of sustained activation of ribosomal S6 kinase (RSK) and ERK by kaposi sarcoma-associated herpesvirus ORF45: multiprotein complexes retain active phosphorylated ERK AND RSK and protect them from dephosphorylation. The Journal of biological chemistry 284: 13958-13968. doi: 10.1074/jbc.m900025200

Reid DW, Nicchitta CV (2012) The enduring enigma of nuclear translation. The Journal of cell biology 197: 7-9. doi: 10.1083/jcb.201202140

David A, Dolan BP, Hickman HD, Knowlton JJ, Clavarino G, et al. (2012) Nuclear translation visualized by ribosome-bound nascent chain puromycylation. The Journal of cell biology 197: 45-57. doi: 10.1083/jcb.201112145

Gallo A, Locatelli F (2012) ADARs: allies or enemies? The importance of A-to-I RNA editing in human disease: from cancer to HIV-1. Biological reviews of the Cambridge Philosophical Society 87: 95-110. doi: 10.1111/j.1469-185x.2011.00186.x

Samuel C (2012) ADARs: viruses and innate immunity. Current Topics in Microbiology and Immunology 353: 163-195. doi: 10.1007/82_2011_148

Bourara K (2000) Generation of G-to-A and C-to-U Changes in HIV-1 Transcripts by RNA Editing. Science 289: 1564-1566. doi: 10.1126/science.289.5484.1564

Klimek-Tomczak K, Mikula M, Dzwonek A, Paziewska A, Karczmarski J, et al. (2006) Editing of hnRNP K protein mRNA in colorectal adenocarcinoma and surrounding mucosa. British journal of cancer 94: 586-592

Zhang F, Hinnebusch AG (2011) An upstream ORF with non-AUG start codon is translated in vivo but dispensable for translational control of GCN4 mRNA. Nucleic acids research 39: 3128-3140. doi: 10.1093/nar/gkq1251

Menschaert G, Van Criekinge W, Notelaers T, Koch A, Crappe J, et al.. (2013) Running title: Molecular and cellular proteomics mcp.M113.0: 1-41

Fritsch C, Herrmann A, Nothnagel M, Szafranski K, Huse K, et al. (2012) Genome-wide search for novel human uORFs and N-terminal protein extensions using ribosomal footprinting. Genome research 22: 2208-2218. doi: 10.1101/gr.139568.112

Ivanov IP, Loughran G, Sachs MS, Atkins JF (2010) Initiation context modulates autoregulation of eukaryotic translation initiation factor 1 (eIF1). Proceedings of the National Academy of Sciences 107: 18056-18060. doi: 10.1073/pnas.1009269107

Loughran G, Sachs MS, Atkins JF, Ivanov IP (2012) Stringency of start codon selection modulates autoregulation of translation initiation factor eIF5. Nucleic acids research 40: 2898-2906. doi: 10.1093/nar/gkr1192

Morris DR, Geballe AP (2000) Upstream Open Reading Frames as Regulators of mRNA Translation. Molecular and Cellular Biology 20: 8635-8642. doi: 10.1128/mcb.20.23.8635-8642.2000

Gaba A, Wang Z, Krishnamoorthy T, Hinnebusch AG, Sachs MS (2001) Physical evidence for distinct mechanisms of translational control by upstream open reading frames. The EMBO journal 20: 6453-6463. doi: 10.1093/emboj/20.22.6453

Hinnebusch AG (1997) Translational Regulation of Yeast GCN4. A window on factors that control initiator-tRNA binding to the ribosome. Journal of Biological Chemistry 272: 21661-21664. doi: 10.1074/jbc.272.35.21661

Palam LR, Baird TD, Wek RC (2011) Phosphorylation of eIF2 facilitates ribosomal bypass of an inhibitory upstream ORF to enhance CHOP translation. The Journal of biological chemistry 286: 10939-10949. doi: 10.1074/jbc.m110.216093

Lincoln AJ (1998) Inhibition of CCAAT/Enhancer-binding Protein alpha and beta Translation by Upstream Open Reading Frames. Journal of Biological Chemistry 273: 9552-9560. doi: 10.1074/jbc.273.16.9552

Cao J, Geballe A (1995) Translational inhibition by a human cytomegalovirus upstream open reading frame despite inefficient utilization of its AUG codon. J Virol 69: 1030-1036

Shabman RS, Hoenen T, Groseth A, Jabado O, Binning JM, et al. (2013) An upstream open reading frame modulates ebola virus polymerase translation and virus replication. PLoS pathogens 9: e1003147. doi: 10.1371/journal.ppat.1003147

Chen A, Kao YF, Brown CM (2005) Translation of the first upstream ORF in the hepatitis B virus pregenomic RNA modulates translation at the core and polymerase initiation codons. Nucleic acids research 33: 1169-1181. doi: 10.1093/nar/gki251

Langmead B, Salzberg SL (2012) Fast gapped-read alignment with Bowtie 2. Nature methods 9: 357-359. doi: 10.1038/nmeth.1923

Bowser BS, DeWire SM, Damania B (2002) Transcriptional Regulation of the K1 Gene Product of Kaposi's Sarcoma-Associated Herpesvirus. Journal of Virology 76: 12574-12583. doi: 10.1128/jvi.76.24.12574-12583.2002

Spiller OB, Robinson M, O'Donnell E, Milligan S, Morgan BP, et al. (2003) Complement Regulation by Kaposi's Sarcoma-Associated Herpesvirus ORF4 Protein. Journal of Virology 77: 592-599. doi: 10.1128/jvi.77.1.592-599.2003

Ozgur S, Damania B, Griffith J (2011) The Kaposi's sarcoma-associated herpesvirus ORF6 DNA binding protein forms long DNA-free helical protein filaments. Journal of structural biology 174: 37-43. doi: 10.1016/j.jsb.2010.10.015

Davison AJ, Stow ND (2005) New genes from old: redeployment of dUTPase by herpesviruses. Journal of virology 79: 12880-12892. doi: 10.1128/jvi.79.20.12880-12892.2005

Neipel F, Albrecth JC, Ensser A, Huang Y, Li JJ, et al. (1997) Human herpesvirus 8 encodes a homolog of interleukin-6. Journal of virology 71: 839-842

Coscoy L, Ganem D (2000) Kaposi's sarcoma-associated herpesvirus encodes two proteins that block cell surface display of MHC class I chains by enhancing their endocytosis. Proceedings of the National Academy of Sciences of the United States of America 97: 8051-8056. doi: 10.1073/pnas.140129797

Rimessi P, Bonaccorsi A, Stürzl M, Fabris M, Brocca-Cofano E, et al. (2001) Transcription pattern of human herpesvirus 8 open reading frame K3 in primary effusion lymphoma and Kaposi's sarcoma. Journal of virology 75: 7161-7174. doi: 10.1128/jvi.75.15.7161-7174.2001

Persson LM, Wilson AC (2010) Wide-scale use of Notch signaling factor CSL/RBP-Jkappa in RTA-mediated activation of Kaposi's sarcoma-associated herpesvirus lytic genes. Journal of virology 84: 1334-1347. doi: 10.1128/jvi.01301-09

Moore PS, Boshoff C, Weiss RA, Chang Y (1996) Molecular Mimicry of Human Cytokine and Cytokine Response Pathway Genes by KSHV. Science 274: 1739-1744. doi: 10.1126/science.274.5293.1739

Cheng EH-Y (1997) A Bcl-2 homolog encoded by Kaposi sarcoma-associated virus, human herpesvirus 8, inhibits apoptosis but does not heterodimerize with Bax or Bak. Proceedings of the National Academy of Sciences 94: 690-694. doi: 10.1073/pnas.94.2.690

Sarid R, Sato T, Bohenzky RA, Russo JJ, Chang Y (1997) Kaposi's sarcoma-associated herpesvirus encodes a functional Bcl-2 homologue. Nature Medicine 3: 293-298. doi: 10.1038/nm0397-293

Chang J, Ganem D (2000) On the control of late gene expression in Kaposi's sarcoma-associated herpesvirus (human herpesvirus-8). J Gen Virol 81: 2039-2047

Arumugaswami V, Wu T-T, Martinez-Guzman D, Jia Q, Deng H, et al. (2006) ORF18 is a transfactor that is essential for late gene transcription of a gammaherpesvirus. Journal of virology 80: 9730-9740. doi: 10.1128/jvi.00246-06

Haque M, Wang V, Davis DA, Zheng Z-M, Yarchoan R (2006) Genetic organization and hypoxic activation of the Kaposi's sarcoma-associated herpesvirus ORF34-37 gene cluster. Journal of virology 80: 7037-7051. doi: 10.1128/jvi.00553-06

Masa S-R, Lando R, Sarid R (2008) Transcriptional regulation of the open reading frame 35 encoded by Kaposi's sarcoma-associated herpesvirus. Virology 371: 14-31. doi: 10.1016/j.virol.2007.08.023

Hamza MS, Reyes RA, Izumiya Y, Wisdom R, Kung H-J, et al. (2004) ORF36 protein kinase of Kaposi's sarcoma herpesvirus activates the c-Jun N-terminal kinase signaling pathway. The Journal of biological chemistry 279: 38325-38330. doi: 10.1074/jbc.m400964200

Koyano S (2003) Glycoproteins M and N of human herpesvirus 8 form a complex and inhibit cell fusion. Journal of General Virology 84: 1485-1491. doi: 10.1099/vir.0.18941-0

Kuang E, Tang Q, Maul GG, Zhu F (2008) Activation of p90 ribosomal S6 kinase by ORF45 of Kaposi's sarcoma-associated herpesvirus and its role in viral lytic replication. Journal of virology 82: 1838-1850. doi: 10.1128/jvi.02119-07

Wang S-S, Chang P-J, Chen L-W, Chen L-Y, Hung C-H, et al. (2012) Positive and negative regulation in the promoter of the ORF46 gene of Kaposi's sarcoma-associated herpesvirus. Virus research 165: 157-169. doi: 10.1016/j.virusres.2012.02.010

Naranatt PP, Akula SM, Chandran B (2002) Characterization of gamma2-human herpesvirus-8 glycoproteins gH and gL. Archives of virology 147: 1349-1370. doi: 10.1007/s00705-002-0813-7

Lukac DM, Renne R, Kirshner JR, Ganem D (1998) 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 252: 304-312. doi: 10.1006/viro.1998.9486

Gruffat H, Portes-Sentis S, Sergeant A, Manet E (1999) Kaposi's sarcoma-associated herpesvirus (human herpesvirus-8) encodes a homologue of the Epstein-Barr virus bZip protein EB1. J Gen Virol 80: 557-561

Li H, Komatsu T, Dezube BJ, Kaye KM (2002) The Kaposi's Sarcoma-Associated Herpesvirus K12 Transcript from a Primary Effusion Lymphoma Contains Complex Repeat Elements, Is Spliced, and Initiates from a Novel Promoter. Journal of Virology 76: 11880-11888. doi: 10.1128/jvi.76.23.11880-11888.2002

Matsumura S, Fujita Y, Gomez E, Tanese N, Wilson AC (2005) Activation of the Kaposi's sarcoma-associated herpesvirus major latency locus by the lytic switch protein RTA (ORF50). Journal of virology 79: 8493-8505. doi: 10.1128/jvi.79.13.8493-8505.2005

Pertel PE, Spear PG, Longnecker R (1998) Human herpesvirus-8 glycoprotein B interacts with Epstein-Barr virus (EBV) glycoprotein 110 but fails to complement the infectivity of EBV mutants. Virology 251: 402-413. doi: 10.1006/viro.1998.9412

Wu FY, Ahn JH, Alcendor DJ, Jang WJ, Xiao J, et al. (2001) Origin-independent assembly of Kaposi's sarcoma-associated herpesvirus DNA replication compartments in transient cotransfection assays and association with the ORF-K8 protein and cellular PML. Journal of virology 75: 1487-1506. doi: 10.1128/jvi.75.3.1487-1506.2001

Bisson SA, Page A-L, Ganem D (2009) A Kaposi's Sarcoma-Associated Herpesvirus Protein That Forms Inhibitory Complexes with Type I Interferon Receptor Subunits, Jak and STAT Proteins, and Blocks Interferon-Mediated Signal Transduction. Journal of Virology 83: 5056-5066. doi: 10.1128/jvi.02516-08

Unal A, Pray T, Lagunoff M, Pennington M, Ganem D, et al. (1997) The protease and the assembly protein of Kaposi's sarcoma-associated herpesvirus (human herpesvirus 8). J Virol 71: 7030-7038

Cannon JS, Hamzeh F, Moore S, Nicholas J, Ambinder RF (1999) Human Herpesvirus 8-Encoded Thymidine Kinase and Phosphotransferase Homologues Confer Sensitivity to Ganciclovir. J Virol 73: 4786-4793

Ohno S, Steer B, Sattler C, Adler H (2012) ORF23 of murine gammaherpesvirus 68 is non-essential for in vitro and in vivo infection. The Journal of general virology 93: 1076-1080. doi: 10.1099/vir.0.041129-0

Wong E, Wu T-T, Reyes N, Deng H, Sun R (2007) Murine gammaherpesvirus 68 open reading frame 24 is required for late gene expression after DNA replication. Journal of virology 81: 6761-6764. doi: 10.1128/jvi.02726-06

Nealon K, Newcomb WW, Pray TR, Craik CS, Brown JC, et al. (2001) Lytic replication of Kaposi's sarcoma-associated herpesvirus results in the formation of multiple capsid species: isolation and molecular characterization of A, B, and C capsids from a gammaherpesvirus. Journal of virology 75: 2866-2878. doi: 10.1128/jvi.75.6.2866-2878.2001

May JS, Walker J, Colaco S, Stevenson PG (2005) The murine gammaherpesvirus 68 ORF27 gene product contributes to intercellular viral spread. Journal of virology 79: 5059-5068. doi: 10.1128/jvi.79.8.5059-5068.2005

May JS, Coleman HM, Boname JM, Stevenson PG (2005) Murine gammaherpesvirus-68 ORF28 encodes a non-essential virion glycoprotein. The Journal of general virology 86: 919-928. doi: 10.1099/vir.0.80661-0

Renne R, Blackbourn D, Whitby D, Levy J, Ganem D (1998) Limited Transmission of Kaposi's Sarcoma-Associated Herpesvirus in Cultured Cells. J Virol 72: 5182-5188

Bai Z, Zhou F, Lei X, Ma X, Lu C, et al. (2012) A cluster of transcripts encoded by KSHV ORF30-33 gene locus. Virus genes 44: 225-236. doi: 10.1007/s11262-011-0698-1

Guo H, Wang L, Peng L, Zhou ZH, Deng H (2009) Open reading frame 33 of a gammaherpesvirus encodes a tegument protein essential for virion morphogenesis and egress. Journal of virology 83: 10582-10595. doi: 10.1128/jvi.00497-09

AuCoin DP, Pari GS (2002) The human herpesvirus-8 (Kaposi's sarcoma-associated herpesvirus) ORF 40/41 region encodes two distinct transcripts. J Gen Virol 83: 189-193

Wang L, Guo H, Reyes N, Lee S, Bortz E, et al. (2012) Distinct domains in ORF52 tegument protein mediate essential functions in murine gammaherpesvirus 68 virion tegumentation and secondary envelopment. Journal of virology 86: 1348-1357. doi: 10.1128/jvi.05497-11

Deng B, O'Connor CM, Kedes DH, Zhou ZH (2007) Direct visualization of the putative portal in the Kaposi's sarcoma-associated herpesvirus capsid by cryoelectron tomography. Journal of virology 81: 3640-3644. doi: 10.1128/jvi.02254-06

Bortz E, Wang L, Jia Q, Wu T-T, Whitelegge JP, et al. (2007) Murine gammaherpesvirus 68 ORF52 encodes a tegument protein required for virion morphogenesis in the cytoplasm. Journal of virology 81: 10137-10150. doi: 10.1128/jvi.01233-06

Cunningham C (2003) Transcription mapping of human herpesvirus 8 genes encoding viral interferon regulatory factors. Journal of General Virology 84: 1471-1483. doi: 10.1099/vir.0.19015-0

Lin S, Sun R, Heston L, Gradoville L, Shedd D, et al. (1997) Identification, expression, and immunogenicity of Kaposi's sarcoma- associated herpesvirus-encoded small viral capsid antigen. J Virol 71: 3069-3076

Yang T-C, Chang C-P, Lan Y-C, Liu C-L, Shih M-C, et al. (2009) Recombinant ORF66 and ORFK12 antigens for the detection of human herpesvirus 8 antibodies in HIV-positive and -negative patients. Biotechnology letters 31: 629-637. doi: 10.1007/s10529-009-9918-3

Desai PJ, Pryce EN, Henson BW, Luitweiler EM, Cothran J (2012) Reconstitution of the Kaposi's sarcoma-associated herpesvirus nuclear egress complex and formation of nuclear membrane vesicles by coexpression of ORF67 and ORF69 gene products. Journal of virology 86: 594-598. doi: 10.1128/jvi.05988-11

Santarelli R, Farina A, Granato M, Gonnella R, Raffa S, et al. (2008) Identification and characterization of the product encoded by ORF69 of Kaposi's sarcoma-associated herpesvirus. Journal of virology 82: 4562-4572. doi: 10.1128/jvi.02400-07

Nador RG, Milligan LL, Flore O, Wang X, Arvanitakis L, et al. (2001) Expression of Kaposi's sarcoma-associated herpesvirus G protein-coupled receptor monocistronic and bicistronic transcripts in primary effusion lymphomas. Virology 287: 62-70. doi: 10.1006/viro.2001.1016

Zhu FX, Chong JM, Wu L, Yuan Y (2005) Virion proteins of Kaposi's sarcoma-associated herpesvirus. Journal of virology 79: 800-811. doi: 10.1128/jvi.79.2.800-811.2005

Cinquina CC, Grogan E, Sun R, Lin SF, Beardsley GP, et al. (2000) Dihydrofolate reductase from Kaposi's sarcoma-associated herpesvirus. Virology 268: 201-217. doi: 10.1006/viro.1999.0165

Gáspár G, De Clercq E, Neyts J (2002) Gammaherpesviruses encode functional dihydrofolate reductase activity. Biochemical and Biophysical Research Communications 297: 756-759. doi: 10.1016/s0006-291x(02)02286-6

González CM, Wong EL, Bowser BS, Hong GK, Kenney S, et al. (2006) Identification and characterization of the Orf49 protein of Kaposi's sarcoma-associated herpesvirus. Journal of virology 80: 3062-3070. doi: 10.1128/jvi.80.6.3062-3070.2006

Gregory SM, Davis BK, West JA, Taxman DJ, Matsuzawa S, et al. (2011) Discovery of a viral NLR homolog that inhibits the inflammasome. Science (New York, NY) 331: 330-334. doi: 10.1126/science.1199478

Wong EL, Damania B (2006) Transcriptional regulation of the Kaposi's sarcoma-associated herpesvirus K15 gene. Journal of virology 80: 1385-1392. doi: 10.1128/jvi.80.3.1385-1392.2006