Pervasive transcription of a herpesvirus genome generates functionally important RNAs.

Pervasive transcription is observed in a wide range of organisms, including humans, mice, and viruses, but the functional significance of the resulting transcripts remains uncertain. Current genetic approaches are often limited by their emphasis on protein-coding open reading frames (ORFs). We previously identified extensive pervasive transcription from the murine gammaherpesvirus 68 (MHV68) genome outside known ORFs and antisense to known genes (termed expressed genomic regions [EGRs]). Similar antisense transcripts have been identified in many other herpesviruses, including Kaposi's sarcoma-associated herpesvirus and human and murine cytomegalovirus. Despite their prevalence, whether these RNAs have any functional importance in the viral life cycle is unknown, and one interpretation is that these are merely "noise" generated by functionally unimportant transcriptional events. To determine whether pervasive transcription of a herpesvirus genome generates RNA molecules that are functionally important, we used a strand-specific functional approach to target transcripts from thirteen EGRs in MHV68. We found that targeting transcripts from six EGRs reduced viral protein expression, proving that pervasive transcription can generate functionally important RNAs. We characterized transcripts emanating from EGRs 26 and 27 in detail using several methods, including RNA sequencing, and identified several novel polyadenylated transcripts that were enriched in the nuclei of infected cells. These data provide the first evidence of the functional importance of regions of pervasive transcription emanating from MHV68 EGRs. Therefore, studies utilizing mutation of a herpesvirus genome must account for possible effects on RNAs generated by pervasive transcription. IMPORTANCE The fact that pervasive transcription produces functionally important RNAs has profound implications for design and interpretation of genetic studies in herpesviruses, since such studies often involve mutating both strands of the genome. This is a common potential problem; for example, a conservative estimate is that there are an additional 73,000 nucleotides transcribed antisense to annotated ORFs from the 119,450-bp MHV68 genome. Recognizing the importance of considering the function of each strand of the viral genome independently, we used strand-specific approaches to identify six regions of the genome encoding transcripts that promoted viral protein expression. For two of these regions, we mapped novel transcripts and determined that targeting transcripts from these regions reduced viral replication and the expression of other viral genes. This is the first description of a function for these RNAs and suggests that novel transcripts emanating from regions of pervasive transcription are critical for the viral life cycle.

Identification of alternative transcripts encoding the essential murine gammaherpesvirus lytic transactivator RTA.

UNLABELLED: The essential immediate early transcriptional activator RTA, encoded by gene 50, is conserved among all characterized gammaherpesviruses. Analyses of a recombinant murine gammaherpesvirus 68 (MHV68) lacking both of the known gene 50 promoters (G50DblKo) revealed that this mutant retained the ability to replicate in the simian kidney epithelial cell line Vero but not in permissive murine fibroblasts following low-multiplicity infection. However, G50DblKo replication in permissive fibroblasts was partially rescued by high-multiplicity infection. In addition, replication of the G50DblKo virus was rescued by growth on mouse embryonic fibroblasts (MEFs) isolated from IFN-α/ßR-/- mice, while growth on Vero cells was suppressed by the addition of alpha interferon (IFN-α). 5' rapid amplification of cDNA ends (RACE) analyses of RNAs prepared from G50DblKo and wild-type MHV68-infected murine macrophages identified three novel gene 50 transcripts initiating from 2 transcription initiation sites located upstream of the currently defined proximal and distal gene 50 promoters. In transient promoter assays, neither of the newly identified gene 50 promoters exhibited sensitivity to IFN-α treatment. Furthermore, in a single-step growth analysis RTA levels were higher at early times postinfection with the G50DblKo mutant than with wild-type virus but ultimately fell below the levels of RTA expressed by wild-type virus at later times in infection. Infection of mice with the MHV68 G50DblKo virus demonstrated that this mutant virus was able to establish latency in the spleen and peritoneal exudate cells (PECs) of C57BL/6 mice with about 1/10 the efficiency of wild-type virus or marker rescue virus. However, despite the ability to establish latency, the G50DblKo virus mutant was severely impaired in its ability to reactivate from either latently infected splenocytes or PECs. Consistent with the ability to rescue replication of the G50DblKo mutant by growth on type I interferon receptor null MEFs, infection of IFN-α/ßR-/- mice with the G50DblKo mutant virus demonstrated partial rescue of (i) acute virus replication in the lungs, (ii) establishment of latency, and (iii) reactivation from latency. The identification of additional gene 50/RTA transcripts highlights the complex mechanisms involved in controlling expression of RTA, likely reflecting time-dependent and/or cell-specific roles of different gene 50 promoters in controlling virus replication. Furthermore, the newly identified gene 50 transcripts may also act as negative regulators that modulate RTA expression. IMPORTANCE: The viral transcription factor RTA, encoded by open reading frame 50 (Orf50), is well conserved among all known gammaherpesviruses and is essential for both virus replication and reactivation from latently infected cells. Previous studies have shown that regulation of gene 50 transcription is complex. The studies reported here describe the presence of additional alternatively initiated, spliced transcripts that encode RTA. Understanding how expression of this essential viral gene product is regulated may identify new strategies for interfering with infection in the setting of gammaherpesvirus-induced diseases.

Murine cytomegalovirus protein pM79 is a key regulator for viral late transcription.

Herpesvirus genes are temporally expressed during permissive infections, but how their expression is regulated at late times is poorly understood. Previous studies indicate that the human cytomegalovirus (CMV) gene, UL79, is required for late gene expression. However, the mechanism remains to be fully elucidated, and UL79 homologues in other CMVs have not been studied. Here, we characterized the role of the conserved murine CMV (MCMV) gene M79. We showed that M79 encoded a protein (pM79) which was expressed with early-late kinetics and localized to nuclear viral replication compartments. M79 transcription was significantly decreased in the absence of viral DNA synthesis but markedly stimulated by pM79. To investigate its role, we created the recombinant virus SMin79, in which pM79 expression was disrupted. While marker-rescued virus grew efficiently in fibroblasts, SMin79 failed to produce infectious progeny but was rescued by pM79 expression in trans. During SMin79 infection, representative viral immediate-early and early gene products as well as viral DNA accumulated sufficiently. Formation of viral replication compartments also appeared normal. Pulsed-field gel electrophoresis analysis indicated that the overall structure of replicating viral DNA was indistinguishable between wild-type and SMin79 infection. Viral tiled array and quantitative PCR analysis revealed that many late transcripts sensitive to a viral DNA synthesis inhibitor (phosphonoacetic acid) were markedly reduced by pM79 mutation. This study indicates that cytomegaloviruses use a conserved mechanism to promote transcription at late stages of infection and that pM79 is a critical regulator for at least a subset of viral DNA synthesis-dependent transcripts.

Identification and sequencing of a novel rodent gammaherpesvirus that establishes acute and latent infection in laboratory mice.

Gammaherpesviruses encode numerous immunomodulatory molecules that contribute to their ability to evade the host immune response and establish persistent, lifelong infections. As the human gammaherpesviruses are strictly species specific, small animal models of gammaherpesvirus infection, such as murine gammaherpesvirus 68 (γHV68) infection, are important for studying the roles of gammaherpesvirus immune evasion genes in in vivo infection and pathogenesis. We report here the genome sequence and characterization of a novel rodent gammaherpesvirus, designated rodent herpesvirus Peru (RHVP), that shares conserved genes and genome organization with γHV68 and the primate gammaherpesviruses but is phylogenetically distinct from γHV68. RHVP establishes acute and latent infection in laboratory mice. Additionally, RHVP contains multiple open reading frames (ORFs) not present in γHV68 that have sequence similarity to primate gammaherpesvirus immunomodulatory genes or cellular genes. These include ORFs with similarity to major histocompatibility complex class I (MHC-I), C-type lectins, and the mouse mammary tumor virus and herpesvirus saimiri superantigens. As these ORFs may function as immunomodulatory or virulence factors, RHVP presents new opportunities for the study of mechanisms of immune evasion by gammaherpesviruses.

Hidden Markov model speed heuristic and iterative HMM search procedure.

BACKGROUND: Profile hidden Markov models (profile-HMMs) are sensitive tools for remote protein homology detection, but the main scoring algorithms, Viterbi or Forward, require considerable time to search large sequence databases. RESULTS: We have designed a series of database filtering steps, HMMERHEAD, that are applied prior to the scoring algorithms, as implemented in the HMMER package, in an effort to reduce search time. Using this heuristic, we obtain a 20-fold decrease in Forward and a 6-fold decrease in Viterbi search time with a minimal loss in sensitivity relative to the unfiltered approaches. We then implemented an iterative profile-HMM search method, JackHMMER, which employs the HMMERHEAD heuristic. Due to our search heuristic, we eliminated the subdatabase creation that is common in current iterative profile-HMM approaches. On our benchmark, JackHMMER detects 14% more remote protein homologs than SAM's iterative method T2K. CONCLUSIONS: Our search heuristic, HMMERHEAD, significantly reduces the time needed to score a profile-HMM against large sequence databases. This search heuristic allowed us to implement an iterative profile-HMM search method, JackHMMER, which detects significantly more remote protein homologs than SAM's T2K and NCBI's PSI-BLAST.

Identification of novel microRNA-like molecules generated from herpesvirus and host tRNA transcripts.

We applied deep sequencing technology to small RNA fractions from cells lytically infected with murine gammaherpesvirus 68 (gammaHV68) in order to define in detail small RNAs generated from a cluster of tRNA-related polycistronic structures located at the left end of the viral genome. We detected 10 new candidate microRNAs (miRNAs), six of which were confirmed by Northern blot analysis, leaving four as provisional. In addition, we determined that previously identified and annotated viral miRNA molecules were not necessarily represented as the most abundant sequence originating from a transcript. Based on these new small RNAs and previously reported gammaHV68 miRNAs, we were able to further describe and annotate the distinctive gammaHV68 tRNA-miRNA structures. We used this deep sequencing data and computational analysis to identify similar structures in the mouse genome and validated that these host structures also give rise to small RNAs. This reveals a possible convergent usage of tRNA/polymerase III (pol III) transcripts to generate small RNAs from both mammalian and viral genomes.

Redefining the genetics of murine gammaherpesvirus 68 via transcriptome-based annotation.

Viral genetic studies typically focus on large open reading frames (ORFs) identified during genome annotation (ORF-based annotation). Here we describe tools for examining viral gene expression nucleotide by nucleotide across the genome. Using these tools on the 119,450 base pair (bp) genome of murine gammaherpesvirus 68 (gammaHV68) allowed us to establish that gammaHV68 RNA expression was significantly more complex than predicted from ORF-based annotation, including over 73,000 nucleotides of unexpected transcription within 30 expressed genomic regions (EGRs). Approximately 90% of this RNA expression was antisense to genomic regions containing known large ORFs. We verified the existence of previously undefined transcripts in three EGRs and determined which parts of the transcriptome depend on protein or viral DNA synthesis. This study redefines the genetic map of gammaHV68, indicating that herpesviruses contain significantly more genetic complexity than predicted from ORF-based genome annotations, and provides alternative tools and approaches for viral genetic studies.

Circular box C/D RNAs in Pyrococcus furiosus.

Box C/D RNAs are small, noncoding RNAs that function in RNA modification in eukaryotes and archaea. Here, we report that box C/D RNAs exist in the rare biological form of RNA circles in the hyperthermophilic archaeon Pyrococcus furiosus. Northern analysis of box C/D RNAs reveals two prominent RNA species of different electrophoretic mobilities in total P. furiosus RNA preparations. Together, the results of Northern, ribozyme, RT-PCR, and lariat debranching analyses indicate that the two species are circular and linear RNAs of similar length and abundance. It seems that most, if not all, species of box C/D RNAs exist as circles in P. furiosus. In addition, the circular RNAs are found in complexes with proteins required for box C/D RNA function. Our finding places box C/D RNAs among the extremely few circular RNAs known to exist in nature. Moreover, the unexpected discovery of circular box C/D RNAs points to the existence of a previously unrecognized biogenesis pathway for box C/D RNAs in archaea.