- Open Access
Widespread activation of antisense transcription of the host genome during herpes simplex virus 1 infection
- Emanuel Wyler†1,
- Jennifer Menegatti†2,
- Vedran Franke†1,
- Christine Kocks1,
- Anastasiya Boltengagen1,
- Thomas Hennig3,
- Kathrin Theil1,
- Andrzej Rutkowski4, 5,
- Carmelo Ferrai1,
- Laura Baer2,
- Lisa Kermas1,
- Caroline Friedel6,
- Nikolaus Rajewsky1,
- Altuna Akalin1,
- Lars Dölken3,
- Friedrich Grässer2Email author and
- Markus Landthaler1, 7Email author
© The Author(s). 2017
- Received: 7 June 2017
- Accepted: 29 September 2017
- Published: 31 October 2017
Herpesviruses can infect a wide range of animal species. Herpes simplex virus 1 (HSV-1) is one of the eight herpesviruses that can infect humans and is prevalent worldwide. Herpesviruses have evolved multiple ways to adapt the infected cells to their needs, but knowledge about these transcriptional and post-transcriptional modifications is sparse.
Here, we show that HSV-1 induces the expression of about 1000 antisense transcripts from the human host cell genome. A subset of these is also activated by the closely related varicella zoster virus. Antisense transcripts originate either at gene promoters or within the gene body, and they show different susceptibility to the inhibition of early and immediate early viral gene expression. Overexpression of the major viral transcription factor ICP4 is sufficient to turn on a subset of antisense transcripts. Histone marks around transcription start sites of HSV-1-induced and constitutively transcribed antisense transcripts are highly similar, indicating that the genetic loci are already poised to transcribe these novel RNAs. Furthermore, an antisense transcript overlapping with the BBC3 gene (also known as PUMA) transcriptionally silences this potent inducer of apoptosis in cis.
We show for the first time that a virus induces widespread antisense transcription of the host cell genome. We provide evidence that HSV-1 uses this to downregulate a strong inducer of apoptosis. Our findings open new perspectives on global and specific alterations of host cell transcription by viruses.
Herpesviruses are a family of enveloped DNA viruses which can infect a wide range of animal species. Herpes simplex virus type 1 (HSV-1) is a human pathogen, together with seven other viruses from the alpha, beta, and gammaherpesvirinae. HSV-1 is prevalent worldwide, with more than half of the population latently infected . Beyond its own clinical importance, HSV-1 is used as a model for herpesviruses due to its easy application in the laboratory. A hallmark of HSV-1 is the “host cell shutoff” during lytic infection . This widespread decrease of the expression of host cell genes is achieved at different stages of the RNA life cycle. First, transcription of host cell genes is impaired by reduced levels  and lower chromatin occupancy of RNA Pol II . RNA splicing is then modulated by the viral protein ICP27 . Finally, matured mRNAs are cleaved by the viral vhs/SOX endonuclease, and thus targeted for degradation [6, 7]. Recently, another RNA-related aberrance was observed in HSV-1-infected cells: about half of the transcribed host cell genes experience transcriptional read-through .
Transcription outside of annotated genes, sometimes called “pervasive transcription”, has attracted increasing interest in the past years [9, 10]. First, this term describes transcription from loci not immediately adjacent to known genes, giving rise to long non-coding RNAs (lncRNAs) or enhancer-derived RNAs (eRNAs). Second, antisense transcripts can originate from promoters of annotated genes [11, 12]. Third, transcripts can fully or partially overlap annotated genes, operating in a wide range of mechanisms, from chromatin remodeling to translation [13–16]. Here, we focus on the second and third categories, which we subsume for clarity of presentation as “antisense transcripts” in our text.
While antisense transcription from viral genomes has been described previously, e.g., in HIV-1 , and was also observed in our sequencing data from the HSV-1 genome, to our knowledge, modulation of host cell antisense transcription by virus infections has not been studied so far. Here, we set out to investigate antisense transcription of the host genome after HSV-1 infection, and identified around 1000 antisense transcripts specifically upregulated upon HSV-1 infection. Notably, most of these transcripts are not listed in existing genome annotations. Antisense transcription starts very early in infection. Using replication inhibitor and knockout viruses we show that overexpression of viral ICP4 alone is sufficient to induce expression of a subset of antisense transcripts. Analysis of published RNA-seq data from other herpesviruses revealed that, while induction of a subset of antisense transcripts is conserved within the Alphaherpesvirinae subfamily, infection with more distantly related herpesviruses does not lead to detectable upregulation of antisense transcripts.
Using a reporter assay, we showed that the sequence region upstream of the BBC3as antisense transcript functions as a promoter induced upon infection. Furthermore, we provide evidence that the induced antisense transcript impairs transcription of the BBC3 sense mRNA in cis. We propose that induction of antisense transcripts represents a previously undescribed strategy of how the virus modulates host gene expression.
Detection of antisense transcripts in RNA-seq data from HSV-1-infected cells
For the detection of antisense transcripts, we additionally relied on directional RNA-sequencing data from human foreskin fibroblasts of rRNA-depleted total RNA and 4-thiourdine pulse-labeled and subsequently purified newly transcribed RNA (4sU-RNA)  (Fig. 1a). For maximal sensitivity, the data from 4sU-RNA was analyzed with a newly developed algorithm to detect antisense transcripts. After initial definition of transcribed regions outside of annotated genes, transcriptional read-throughs were filtered out (see “Methods”). Using a running sum algorithm, we defined 12,863 transcribed regions either annotated as antisense transcripts in Ensembl or not annotated at all. Of these, 9765 regions were designated as read-through regions. The remaining 3098 loci were grouped by their expression profiles into ten clusters (Fig. 1b; Additional file 1: Figure S1a; Additional file 2: Table S1). Based on their location, transcripts were categorized as divergent, convergent, or internal antisense transcripts (Fig. 1d), or, if the transcription start site (TSS) of the antisense transcript was not within 10 kb of a TSS of a transcript on the opposite strand, and if it did not overlap a gene body, it was categorized as “undefined”.
Of the 3098 antisense transcripts, 1517 showed at least two-fold upregulation upon infection. Comparing them with a recently published, comprehensive atlas of human non-coding RNAs , we found that only 384 of these begin within less than 100 bp from an annotated TSS, indicating that most of the upregulated transcripts were previously undescribed. As 503 of the 1517 upregulated transcripts were categorized as “undefined”, this leaves around 1014 antisense transcripts upregulated by the HSV-1 infection.
To determine when during infection the antisense transcripts emerge, we clustered the correlation coefficients of their expression values at the different timepoints (Fig. 1c). As expected, biological replicates were highly correlated and the clustering matched the time course of infection. Importantly, sequencing data obtained from RNA labeled 1 to 2 h after infection was separated from the mock and 0–1-h timepoint, indicating that antisense transcripts are already expressed very early in infection.
Transcripts selected for detailed analysis
Convergent, spliced 5′ UTR, asTSS in 5′ UTR intron
Convergent, spliced 5′ UTR, asTSS in 5′ UTR intron
Convergent, spliced 5′ UTR, asTSS in 5′ UTR intron
Overlap or divergent from ALDH4A1
Convergent, spliced 5′ UTR, asTSS in 5′ UTR intron
Convergent, spliced 5′ UTR, asTSS in 5′ UTR intron
Convergent, spliced 5′ UTR, asTSS in 5′ UTR intron
Generally, we observed that the antisense transcripts were not spliced. Noteworthy, antisense transcripts were found in the poly(A)-selected RNA, suggesting that they are polyadenylated.
The identified antisense transcripts can overall be classified into: a) divergent antisense transcripts, where the sense and antisense transcripts likely start from the same promoter region, but do not overlap; b) convergent antisense transcripts, where the 5′ ends of the antisense and the canonical sense transcript overlap; and c) internal antisense transcripts, where several exons of the canonical sense transcript are overlapping with the antisense transcript.
Taken together, we detected 3098 novel antisense transcripts in strand-specific RNA sequencing data, thereby expanding the catalog of lncRNAs . Of these antisense transcripts, 1014 showed increased expression upon HSV-1 infection.
Validation and expression dynamics of antisense transcripts
Confirming the RNA sequencing results, Nanostring counts for antisense transcripts increased during infection (Fig. 2b; Additional file 1: Figure S2c). However, some remarkable differences between cell lines and individual transcripts were observed. The BBC3as transcript (Fig. 2b, top panel) already started to accumulate 1 h after infection in HeLa cells, yet it eventually reached higher levels in the two fibroblast cell lines. Generally, a delayed induction of the antisense transcripts could be monitored in these cells, reflecting the slower progression of infection, which is also observed for the housekeeping genes (Fig. 2a, right). RFX1as and SLC27Aas are the only antisense transcripts of the 12 tested where strong cell line specificities could be observed. RFX1as was only induced in HeLa cells (Fig. 2b, fourth panel), while SLC27A4as was already present in untreated HeLa cells, as shown by comparison of normalized Nanostring counts in the different cell lines (Fig. 2c). Still, there was a relative increase in expression upon infection, although it was considerably lower compared to WI-38 and NHDF. Interestingly, the beginning of the SLC27A4 antisense transcript partially overlapped in some annotations with a tRNA gene. While the coverage profile did not suggest that the antisense transcript is a tRNA, there was no apparent transcription termination at the RNA Pol III terminator. This observation could indicate that the tRNA and the antisense transcript share a common promoter region.
We further validated the antisense transcripts using RT-qPCR from uninfected (0 h) and infected HeLa cells (Additional file 1: Figure S2a). Many of the tested transcripts could not be reliably detected in uninfected cells. Therefore, we normalized qRT-PCR data to a Drosophila melanogaster total RNA spike-in and expression levels at the 2-h timepoint. We included two host MAP kinase target genes as well as pre-miR-183, which are induced during lytic HSV-1 infection [20, 21].
We also considered whether the observed antisense transcripts are a consequence of the reported widespread transcriptional read-through in HSV-1-infected cells . First, we compared the read-through transcripts with the antisense transcripts in the 4sU-seq data (Additional file 1: Figure S2c). There, the two transcript classes showed different expression dynamics. Furthermore, two read-through transcripts were measured using RT-qPCR from the same samples used in Additional file 1: Figure S2a (Additional file 1: Figure S2d). We observed that for these two targets, the qPCR amplicons could only be reproducibly detected at 4 hpi (data not shown), with a strong signal increase from 4 to 6 hpi. Again, this dynamic behavior is different from all antisense transcripts, which also suggests that the antisense transcripts under scrutiny here are not products of the transcriptional read-through observed in HSV-1-infected cells.
To address the question of whether the induced antisense transcripts might also be transcribed in uninfected cells, we used NET-seq data , which is currently one of the most sensitive methods for detection of transcription. RPKM values were calculated from the HeLa total NET-seq control samples, and log2 transformed numbers are plotted in Fig. 2d. A threshold of 3 was applied to separate noise values. Out of 1014 induced antisense transcripts, 762 (75%) have a log2 RPKM in NET-seq larger than this threshold. The majority of induced antisense transcripts are therefore already transcribed in uninfected cells at levels comparable to protein coding genes, but are not detectable in high-throughput sequencing of steady state or newly synthesized RNA (4sU-seq).
Finally, we investigated whether the antisense transcripts could be stabilized upon HSV-1 infection. To this end, we analyzed nucleoplasmic RNA in HeLa cells upon depletion of the exosome component EXOSC3 (also known as RRP40) . Indeed, we could observe at least twofold upregulation of 629 of the 1517 transcripts upregulated in HSV-1-infected cells, among them, e.g., BBC3as (Additional file 1: Figure S2e). However, our antisense detection algorithm identified 12,417 antisense transcripts upregulated in EXOSC3-depleted cells compared to control cells, among them the recently described promoter upstream transcripts [23, 24]. This discrepancy indicates that reduced exosome activity is unlikely a cause of the induction of antisense transcription by HSV-1 infection.
Taken together, we validated the presence of antisense transcripts using two low-throughput methods in three different cell lines. Notably, SLC27A4as was also detected in uninfected HeLa but not NHDF or WI-38 cells. Furthermore, antisense transcripts show expression dynamics distinct from poly(A) read-through and thus represent an independent entity of transcripts.
Viral factors induce antisense transcripts by different mechanisms
To that end, we used published RNA-seq data for the alphaherpesvirus varicella zoster virus (VZV)  as well as the betaherpesvirus human cytomegalovirus (HCMV) , and the gammaherpesvirus Kaposi’s sarcoma-associated herpesvirus (KSHV) . Due to differences in sequencing depths, and since the VZV sequencing is unstranded, a fully quantitative comparison was not possible. For KSHV-infected cells we did not find any indication for the antisense transcripts induced by HSV-1. In VZV-infected differentiated keratinocytes, however, we observed induction of several hundred antisense transcripts. In HCMV-infected cells only a few antisense transcripts were detected, with the important exception of SLC27A4as, which was shown to be ICP4-independent (Fig. 3c). Within the Alphaherpesvirinae family, the Herpes simplex genus includes a second member that infects humans, namely HSV-2. We therefore probed RNA from HSV-2-infected NHDF cells with our Nanostring assay and observed a similar pattern of antisense transcripts compared to infection with HSV-1 (Additional file 1: Figure S3c).
In order to test how the expression of the antisense transcripts depended on the progression of lytic infection, we used the HSV-1 DNA replicase inhibitor phosphonoacetic acid (PAA), which prevents expression of viral late genes. Consequently, levels of viral transcripts were much lower (Fig. 3b, first panel). Due to the attenuated infection, HSV-1 induced host shutoff, as measured by RNA levels of the housekeeping genes POLR1B, SDHA, TBP, and PPIA, which were reduced. While the values with and without PAA at 4 hpi were comparable, at 6.5 and 9 hpi these transcripts were about twofold to threefold higher, respectively (Fig. 3b, second panel). Interestingly, most, but not all, antisense transcripts were more abundant with PAA (Fig. 3b, third panel; Additional file 1: Figure S3a). Some antisense transcripts, such as BBC3as, FOXO3as, IER2as, and ING1as, were several times more abundant in PAA-treated cells (Fig. 3a), while SLC27A4as showed the opposite behavior, with higher levels in the absence of PAA.
The corresponding sense transcripts showed diverse response patterns (Fig. 3b, fourth panel; Additional file 1: Figure S3a). IER2 is induced by the MAP kinase pathway , which is in turn activated by the infection itself . Similarly to IER2, the pro-apoptotic gene ING1  was induced by MAP kinases as well, and, like IER2, showed higher mRNA levels at all timepoints after infection compared to uninfected cells under PAA treatment. This result indicates an induction of these genes by the virus infection, which is, however, counteracted by the host cell shutoff.
In addition, we measured the expression profiles in cells infected with ICP0 or ICP4 deletion viruses. ICP0 functions as a ubiquitin ligase and contributes to an active viral chromatin structure . ICP4 is mainly considered as a regulator of transcription of the viral genome, interacting with several host cell transcription factors [32, 33]. The ICP0 knockout virus (deletion of the full gene) leads to an attenuated infection, i.e., slower progression and viral replication, but otherwise all viral genes are expressed and functional virions are produced. Without functional ICP4, however, the virus does not progress past expression of immediate early genes (the ICP4 deletion virus contains a short deletion in the N-terminus of the RS1 gene, leading to a frameshift and non-functional protein). In our assay, the late genes UL6 and UL56 consequently accumulated more slowly in the ICP0 deletion virus compared to the wild-type virus, and were absent in the ICP4 deletion virus. ICP4/RS1, itself an immediate early gene, was constantly induced and accumulated (Fig. 3c, first panel). The host cell shutoff proceeded in a comparable manner between the two deletion strains, but slower than in the wild-type virus (Fig. 3c, second panel).
The antisense transcripts again showed different profiles (Fig. 3c, third panel; Additional file 1: Figure S3b). BBC3as, C1orf159as, and FOXO3as were more abundant upon infection with the ICP0 deletion virus compared to the wild-type virus; the difference was similar to that with and without PAA (Additional file 1: Figure S3a, b). Again, IER2as and SLC27A4as transcripts showed distinct expression profiles. For both antisense transcripts, induction was stronger in the ICP4 knockout virus compared to the wild-type virus, with the signals from the two knockout viruses being similar to each other. Because only viral immediate early genes are expressed when using the ICP4 knockout virus, and IER2as and SLC27A4as were detectable under this condition, this difference indicates that IER2as and SLC27A4as are induced by immediate early genes, or by a cellular mechanism triggered by the infection. On the other hand, transcription of BBC3as, C1orf159as, and FOXO3as is likely turned on by early genes or by ICP4 itself.
Overexpression of ICP4 is sufficient for induction of some antisense transcripts
To further investigate whether ICP4 directly binds to the promoter region of target genes, we performed ICP4 ChIP-qPCR, as previously performed for RCL1 . The strongest enrichment for the ICP4-ChIP was indeed found in both cases at the beginning of the antisense transcript (Fig. 4b).
Histone marks are similar for induced and constitutively active antisense transcripts
Antisense transcripts appear very early in infection and can be found widespread throughout the host cell genome. To address the question of whether the virus starts transcription from primed regions, we went on to examine ENCODE histone marks  around TSSs. Figure 4c depicts HeLa histone marks of active promoters (H3K4me3, H3K9ac, H3K27ac) and actively transcribed gene bodies (H3K36me3) of a unidirectional gene, POLR1B, and FOXO3, which bears a HSV-1 inducible divergent antisense transcript. While the POLR1B gene showed the expected profile for a unidirectional gene, the island of activating promoter histone marks was much larger around the FOXO3 TSS, reminiscent of bidirectional promoters . We extended this analysis to include all detected antisense transcripts (Fig. 4d; Additional file 1: Figure S4C). In total, five sets of genes were examined. First, 1000 unidirectional control genes were analyzed. For the convergent and divergent antisense transcripts detected here, we subsequently analyzed histone marks around the TSS of induced and uninduced antisense transcripts separately.
Control genes showed a strong single peak at the TSSs for the three active promoter marks, and a strong increase of H3K36me3 occupancy for active gene bodies downstream of the TSS. Genes having divergent antisense transcripts showed the expected bimodal distribution for active promoters, but, remarkably, the induced and constitutive transcripts did not show any difference. This observation was true for H3K36me3 as well, which showed high levels in both directions from the TSS.
Genes with convergent transcripts also showed a single peak for the active promoter marks, but with a shoulder downstream of the TSS, which could reflect the transcription start for the antisense transcript. Again, inducible and constitutive transcripts showed the same profiles. These results indicate that active promoter marks enabling antisense transcription are already present in the uninfected cells.
BBC3as is transcribed from a HSV-1-dependent promoter and inhibits BBC3 sense in cis
To investigate the physiological role of antisense transcripts, we focused on BBC3as, which was among the most highly expressed, and dynamically induced antisense transcripts (Additional file 1: Figure S1b). BBC3as is antisense to the strong apoptosis inducer BBC3, also known as PUMA . Concomitantly, viruses try to counter cell death in different ways to enable maximal replication .
While the signal from the CCNA2 probes likely represented monodispersed single RNA molecules, the bright nuclear foci for NEAT1 and BBC3as likely indicate the accumulation of multiple RNA molecules . Because HeLa cells are triploid for the BBC3 locus , three foci were visible in HeLa cells compared to two in NHDF cells, suggesting BBC3as accumulated at the site of transcription (Fig. 5b). In summary, single molecule RNA FISH data indicated that the BBC3as transcript accumulates at discrete foci in HSV-1-infected cells, likely at the locus which encodes the pro-apoptotic gene BBC3 . In such a scenario, BBC3as might inhibit BBC3 at the level of transcription. Ectopic overexpression of the first 2191 nucleotides of the BBC3as led to a strong signal in multiple foci in the cytoplasm (Additional file 1: Figure S5a), a behavior that was also previously observed for lincRNA-p21 .
Using Nanostring quantification, BBC3 mRNA exhibited a stronger decrease in expression than mRNA derived from housekeeping genes, both with and without PAA treatment (Fig. 3b). Because BBC3as is induced early in infection (Additional file 1: Figure S2a; Fig. 3b), we reasoned that BBC3 is subject to active downregulation during the infection. For closer examination of BBC3 mRNA transcription, we considered the 4sU-seq data sampled at different timepoints after infection. Comparing the sense and antisense transcripts of BBC3 and the four other genes displayed in Fig. 3b, we observed that, after an initial increase at early timepoints, the expression for BBC3 was lower than for housekeeping genes (Fig. 5c). This observation again indicates that an active mechanism could impair BBC3 mRNA transcription, with the antisense transcript as a prime candidate to mediate a possible transcriptional interference.
To define the promoter region of the BBC3 antisense transcript, we performed dual luciferase promoter reporter assays using a region upstream of the BBC3as transcript, which contains two TATA boxes. Different genomic fragments from that region (Fig. 5d) were inserted upstream of an otherwise promoter-less luciferase gene. Luciferase activity was measured in uninfected and HSV-1-infected HeLa cells, showing a luciferase signal increase upon infection (Fig. 5e). This result confirms the existence of a promoter that is activated by the infection and would lead to the transcription of the BBC3as transcript.
Because our attempts to deplete the antisense transcript using siRNAs or gapmers were unsuccessful (data not shown), we went on to deplete BBC3as by excising the promoter region using Cas9 in combination with two different guide RNA pairs (5/6 and 5/9; Fig. 5d) from HeLa cells. As expected, BBC3as was largely undetectable in the two independently generated BBC3as promoter knockout cell lines (termed ∆BBC3as 56 and ∆BBC3as 59), as measured by RT-qPCR from total RNA (Fig. 5f). By contrast, viral UL29 mRNA levels were largely comparable between the HeLa and knockout cell lines. We then measured the expression levels of different transcripts in the newly transcribed RNA population using 4sU-tagging [43, 44]. C1ORF159 and IER2 levels change comparably upon infection in wild-type HeLa and BBC3as knockout cells, as measured by RT-qPCR (Fig. 5g) and Nanostring nCounter assays (Fig. 5h). However, the BBC3 mRNA synthesis showed a marked decrease upon infection in wild type, but not in BBC3as knockout cells. To further investigate the effect of the BBC3 antisense transcript on its sense counterpart, we performed directional high-throughput sequencing using RNA from wild-type and BBC3as knockout cells at 0, 2, 4, and 6 h after infection in two biological replicates. In Fig. 5i, one replicate for all three cell lines at 0 and 6 h after infection is shown. Only in the BBC3as knockout cells was there a remarkable increase of sequencing coverage in the sense direction in the first, slowly processed  half of the BBC3 gene upon infection. This result suggests that, in infected cells, BBC3 pre-mRNA is increasingly transcribed, but only in the absence of antisense transcription.
Taken together, we provide evidence that the BBC3as transcript is expressed upon HSV-1 infection and inhibits transcription of BBC3 in cis, and therefore enables the virus to suppress this potent inducer of apoptosis.
NF-kappaB transcription factor genes show age-dependent induction of antisense transcripts in lipopolysaccharide-treated monocytes
NF-kappaB antisense transcripts after LPS treatment and HSV-1 infection
LPS antisense upregulated
LPS sense upregulated
HSV-1 antisense upregulated
HSV-1 sense upregulated
Adult and older adults yes, newborn no
Adult and older adults yes, newborn no
Older adults weakly, adults and newborn no
In monocytes from newborns, only the sense NFKB2 is induced, to an extent comparable to adults and older adults, but not NFKB2as. The NFKB1/p105 antisense transcript shares the same properties. Interestingly, while in LPS-stimulated monocytes the respective host genes were upregulated together with the antisense transcripts, in HSV-1 infections only the antisense transcripts were upregulated. Among the three other NF-KB factors, REL has an annotated antisense transcript (official gene symbol LINC01185) that was weakly induced in monocytes from older adults, but not in any of the other samples examined here. RELA and RELB finally showed antisense transcripts induced in HSV-1-infected cells, but not in any of the monocyte samples. Again, all three Rel genes were induced in LPS-treated monocytes, but not in HSV-1-infected cells.
In summary, we show age-related antisense transcript induction for genes of the NF-kappaB transcription factor family in LPS-treated monocytes.
An initial analysis of RNA-seq data from HSV-1-infected human cells brought to our attention that, upon infection, antisense transcripts emerge at several loci in the host cell genome. An in-depth analysis of our own and published sequencing data showed indeed that this induction of antisense transcripts occurs throughout the genome, and 1517 such RNA species, most of them novel and not previously annotated, are induced or upregulated upon virus infection.
Validation of the discovered antisense transcripts using RT-qPCR and Nanostring nCounter assays showed different expression dynamics of antisense transcripts and their respective sense transcripts. However, we did not find significant correlations between sense and antisense transcription, indicating that, in general, their mechanisms for transcription initiation are independent of each other.
We investigated two possibilities concerning whether the observed antisense transcripts arise due to the perturbation of transcription-related processes. First, the antisense transcripts could emerge following stabilization of unstable transcripts. Second, they could show up due to impaired transcription termination , which could render short unstable transcripts longer and, concomitantly, also more stable.
Because we relied on sequencing data of pulse-labeled RNA for the bioinformatic detection of antisense transcripts, short-lived transcripts were preferentially detected and therefore could also be observed in uninfected cells.
To address the question of whether the antisense transcripts are in some way a consequence of the infection-induced transcriptional read-through, we looked at the expression dynamics of read-through transcripts and found that they emerge later in infection, but then with a very steep increase (Additional file 1: Figure S2c, d). Nevertheless, antisense transcripts that come up later in infection (e.g., cluster 5 in Fig. 1b) might still be a consequence of transcriptional read-through. We are currently dissecting the causes for this observation, which might provide a more comprehensive answer to the relationship between the two transcriptional phenomena observed in HSV-1-infected cells. In a broader context, the observation of antisense transcription and transcriptional read-through in HSV-1-infected cells could also help to elucidate the mechanisms underlying promoter directionality and transcription termination . Importantly, some of the antisense transcripts appear to be upregulated in nucleoplasmic RNA from EXOSC3-depleted cells . However, many antisense transcripts that are upregulated by exosome depletion become expressed in HSV-1-infected cells. A reduction of exosome components upon infection was not observed, on either the RNA  protein level , indicating a more complex relationship between exosome activity and HSV-1 infection.
The clustering of the expression dynamics of the antisense transcripts (Fig. 1b) already suggested that certain subsets were induced by different mechanisms. We therefore applied modified infection conditions, namely the replication inhibitor PAA and two knockout viruses (Fig. 3). This accentuated differences in expression dynamics of the antisense transcripts, like for SLC27A4as or IER2as. To further explore the different mechanisms causing antisense transcription, we first examined sequencing data from VZV, which belongs to the genus Varicellovirus within the alphaherpesviruses family, HCMV (beta), and KSHV (gamma). While we did not detect HSV-1-induced antisense transcripts in KSHV-infected cells, a subset of antisense transcripts was induced in VZV-infected cells. The SLC27A4as transcript, apparently independent of the viral protein ICP4 (Fig. 3), was one of the few induced also in CMV-infected cells. This is consistent as ICP4 is only present in alphaherpesviruses. Ectopic expression of ICP4 showed that this viral protein is capable of binding to the promoter regions and inducing transcription of some antisense transcripts independently of viral infection (Fig. 4). It should be mentioned, however, that in all assays (RNA-sequencing, RT-qPCR, Nanostring), the levels of antisense transcripts were considerably higher upon infection compared to ICP4 overexpression. While this could be due to technical issues, such as the activity and amount of the ectopically expressed ICP4, it might also indicate that other factors enhance the activity of ICP4 during infection. On a side note, no evidence for the involvment of viral ICP27 was found in recently published  RNA-sequencing data following overexpression of this protein (data not shown).
Taken together, the results indicate that different classes of antisense transcripts are induced by different mechanisms. While some of the antisense transcripts appear to be exclusive for the Herpes simplex genus, others may be more broadly conserved as they are observed in VZV- and partially in CMV-infected cells.
One of the most highly expressed antisense transcripts was BBC3as, which overlaps about half of the BBC3 gene body. Its TSS lies in the last intron of the BBC3 sense transcript. FISH showed strong nuclear foci comparable to the abundant lncRNA NEAT1, with a size and intensity that represents accumulation of multiple molecules. In addition, we found three such foci in infected HeLa cells, which are triploid at this locus , and two foci in most infected NHDF cells (Fig. 5b). Both observations are compatible with BBC3as accumulating at its site of transcription, maybe by binding to nascent BBC3 sense mRNA. To test whether BBC3as transcriptionally antagonizes transcription in the sense direction, we first analyzed gene expression in HSV-1-infected cells treated with the virus replication inhibitor PAA. Upon PAA treatment, infection progresses much more slowly, and the host cell shutoff was attenuated. In such conditions, we observed that the BBC3 mRNA levels decreased quickly, and that the BBC3as was transcribed at greater levels (Fig. 3b). The analysis of the 4sU-seq time course data, which showed how the syntheses of BBC3 sense and antisense transcripts are anti-correlated (Fig. 5c), confirmed this observation. Interestingly, we also observed a slight induction of the BBC3 sense RNA immediately upon infection, as in the RT-qPCR and Nanostring experiments (Figs. 2b and 3c; Additional file 1: Figure S2a). This is not unexpected, as the infection initially triggers apoptosis, which is then actively counteracted by the virus . Furthermore, we provide evidence that the BBC3as downregulates transcription of the BBC3 mRNA on the transcriptional level (Fig. 5g–i). Taken together, this could indicate a completely novel and unexpected mechanism of how HSV-1 impairs apoptosis: by inducing an antisense transcript that transcriptionally downregulates the potent pro-apoptotic BBC3. A similar mechanism was recently shown in the plant Arabidopsis thaliana, where the COOLAIR antisense transcript blocks transcription of the corresponding sense transcript FLC in cis . Notably, the BBC3 gene experiences transcriptional regulation by CTCF, which binds to the first half of the BBC3 gene, thereby inhibiting mRNA processing . The excised fragment in the last intron in the BBC3as knockout cells appears, however, to be outside this region, also judged by the occurrence of CTCF binding sites . This result indicates that the CTCF regulation is not affected in the cell lines used in our study.
Many of the antisense transcripts observed here start in close proximity to the transcription start sites of their respective host genes, giving rise to divergent or convergent transcripts. The cause, importance and function of these so-called bi-directional promoters is currently under discussion [11, 52]. In the course of this study, we asked two questions in this context: to what extent are the antisense transcripts transcribed in uninfected cells? And can the antisense transcripts induced by HSV-1 infection also be found in other biological settings? For an initial investigation of the first question, we analyzed published NET-seq data, nucleoplasmic RNA upon exosome depletion, and histone marks. Although the antisense transcripts are not detected in steady-state or newly synthesized RNA (4sU-seq), 75% of them are strongly transcribed in uninfected HeLa cells according to the NET-seq data (Fig. 2d). Furthermore, upon exosome depletion, e.g., BBC3as becomes also visible in nucleoplasmic RNA. On a similar note, we compared histone marks between induced and non-induced divergent and convergent antisense transcripts and found that the profiles look essentially the same (Fig. 4).
This could indicate that, also in human cells, antisense transcripts are synthesized, and their expression can vary between conditions, as also shown for the antisense transcripts belonging to the NF-kappaB transcription factor family. For example, BBC3as, which might counteract transcription of the BBC3 mRNA, could be induced upon HSV-1 infection or in cancerous cell lines such as HeLa to prevent apoptosis. Still, the question remains why the antisense transcripts becomes detectable in steady-state poly(A) + RNA only upon infection. Investigating this in the context of the virus infection could provide highly interesting insights in RNA processing and degradation of antisense RNA and gene regulation by lncRNAs.
Our results indicate that HSV-1 infections provide a highly valuable and broadly and easily applicable tool to induce a wide range of antisense transcripts. Our study therefore provides a basis for further investigations of this type of transcript in other biological settings.
Cell culture and viruses (HeLa, WI-38, and NHDF cells)
HeLa cells and HeLa BBC3 as promoter knockout cells were maintained in Dulbecco’s modified Eagle medium (DMEM) supplemented with 10% fetal bovine serum (FBS), 100 units/ml penicillin, and 100 μg/ml streptomycin at 37 °C and 5% CO2. Normal human dermal fibroblasts (NHDFs) and WI-38 fibroblasts were cultured in minimal essential medium (MEM) supplemented with 10% FBS, 100 units/ml penicillin, and 100 μg/ml streptomycin. Wild-type HSV-1 was derived from a patient isolate and was obtained from Department of Virology, Saarland University Medical School, Homburg, Germany. The virus was propagated in HeLa cells or NHDF cells. At 72 hpi, viral supernatant was sterile-filtered through a 0.45 μm pore size filter and stored at −80 °C. Viral titers were determined by plaque assay. Briefly, HeLa cells were inoculated with 500 μl of virus-containing supernatant and serial dilutions in DMEM and incubated for 60 minutes at 37 °C and 5% CO2. The supernatant was then replaced by fresh complete DMEM containing 2% methylcellulose. Four to five days post-infection cells were fixed with 4% paraformaldehyde and stained with crystal violet. The PFU was then determined by counting the visible plaques and the experiments were carried out in triplicate. Cells were infected with an MOI of 10.
Plasmids and transfections
Plasmid pEYFP-ICP4 was a kind gift from Chris Boutell . For overexpression of the BBC3as transcript, a genomic region spanning the first 2191 bases of the antisense transcript was amplified using primers BBC3as_pcdna_fwd1/BBC3as_pcdna_rev2 and inserted into pcDNA3.0. Transfections were performed using FuGene HD (Promega). The plasmid expressing the ICP4-3xflag construct was derived from pEYFP-ICP4 by first replacing the EYFP by the oligos EYFP_NheI_SalI_fwd and EYFP_NheI_SalI_rev after NheI/SalI digestion, and then insertion of the ICP4Cterm-3xflag fragment using SgsI/KpnI (Additional file 3: Table S2).
RNA isolation and sequencing
RNA was extracted from cultured cells with Trizol (Thermo Fisher Scientific) and purified from the aqueous phase using the RNA Clean & Concentrator 25 kit (Zymo Research). For sequencing from WI-38 cells, poly(A) RNA was selected using the Dynabeads mRNA DIRECT Kit (Thermo Fisher Scientific). Sequencing libraries were prepared using the NEXTflex Directional RNA-Seq Kit dUTP-based (Bioo Scientific) and sequenced on a Illumina HiSeq 2000 device. After trimming the adapter sequence AGATCGGAAGAGCACACGT, about 25–30 million mappable reads per sample were obtained. For sequencing from EYFP or EYFP-ICP4 transfected HeLa cells as well as HSV-1-infected HeLa and BBC3as knockout cells, poly(A) RNA sequencing libraries were prepared using the TruSeq Stranded mRNA Library Prep Kit (Illumina). For the HeLa infection time-course samples, a 5% D. melanogaster spike-in was added on the level of total RNA. EYFP/EYFP-ICP4 samples were sequenced on an Illumina NextSeq 500 device; infected HeLa samples were sequenced on a Illumina HiSeq 4000 device. After trimming the adapter sequence GATCGGAAGAGCACACGT, about 21–25 million mappable reads per sample were obtained for the EYFP/EYFP-ICP4 samples, and 12–18 million mappable reads for the HeLa/BBC3as knockout samples.
For visualization of RNA-sequencing data, alignment files were normalized by the number of mapped reads, and, after conversion to BigWig files, displayed using IGV . Coverage in the minus direction is displayed as negative values. To compare different samples, the scale was set to the same values for the same strand within a group of samples. Such groups are, e.g., human foreskin fibroblast (HFF) cells on one side and WI-38 cells on the other side in Fig. 1, or WI-38 cells on one side and the HeLa wild-type/BBC3as knockout cells on the other side in Fig. 5.
Bioinformatical identification and definition of antisense transcripts
Data for definition of antisense transcripts were acquired from the herpes-infected 4sU time-course data . For transcript definition, we used all samples except 6–7 and 7–8 hpi. The data were mapped onto the hg19 reference using TopHat2 . Because of the unique transcriptomic profiles in human cells infected with herpes simplex virus, antisense transcripts were defined on the coverage tracks of the HSV-1-infected 4sU samples using a custom running sum algorithm, followed by extensive filtering. For this purpose, all regions that overlapped known Ensembl protein coding, pseudogene, miRNA, lincRNA, and snoRNA loci were first set to have coverage of 0. Regions overlapping lincRNAs which have a TSS within 2 kb upstream of known protein coding genes were regarded as wrongly annotated antisense genes and were not filtered out. To locate the expressed loci, a running sum algorithm was run on the aforementioned tracks (in a strand-specific fashion). The algorithm scans the genome base by base. When it encounters a base that has coverage > 0, it marks a start of the putative region and increments a counter by a predefined value. With each base with coverage > 0 the counter is incremented, and with each base with coverage = 0 the counter is decremented by a predefined value. When the counter reaches 0, it marks the end of a putative region. The transcript definition was done for each timepoint individually. All transcripts that are supported by at least one read per base were kept for further analysis. The regions resulting from individual stages were merged in the following way. 1) The start of the region was set as the TSS that was supported in most of the time-course samples. 2) The width of the region was defined as the median widths of all the putative transcripts that overlapped a certain locus.
Putative antisense regions were filtered in the following way: all regions that were transcribed 1 kb downstream from the transcription termination sites of Ensembl annotated genes were marked as read-through transcripts and were omitted from further analysis.
All regions with length < 100 bp were omitted from further analysis.
Annotation and expression analysis
Divergent—TSS of the antisense transcript is located within 10 kb upstream of the sense gene
Convergent—TSS of the antisense transcript is located 10 kb downstream of the sense gene
Internal—TSS of the antisense transcript is located in the gene body of the antisense transcript, more than 10 kb downstream from the gene TSS.
Annotation with CAGE data and lincRNA base
TSSs of antisense transcripts were overlapped with Fantom5 CAGE data to ascertain how many regions are known to produce capped transcripts. Additionally, the transcripts were overlapped with known lincRNAs from lincRNAdb.
Analysis of antisense expression
To estimate the expression of the antisense transcripts, we first filtered the mapped reads to contain only properly paired read pairs. We then counted the number of reads that overlap each of the antisense transcripts from each of the 4sU and Ribo Zero RNA-seq samples from Rutkowski et al. . The final expression estimates were obtained by normalizing counts for transcript length, and library size, by calculating the size factors from all protein coding genes. The expression for each timepoint was obtained by taking the log2 of timepoint expression/uninfected sample.
For RT-qPCR from total RNA, DNase digestion before the RT was applied to minimize background using DNase I, Amplification Grade (Thermo Fisher Scientific). Per reaction, 1 μg total RNA was used and transcribed using SuperScript III with random hexamer primers (both Thermo Fisher Scientific) according to the manufacturer’s protocol. For the fractioned, 4sU labeled samples, DNase treatment was not applied and 80 ng RNA was used per RT reaction. For qPCR, Power SYBR Green PCR Master Mix (Thermo Fisher Scientific) was used together with 250 nM of each primer and 1:40 final dilution of the RT reaction. Primers used for qPCRs are detailed in Additional file 3: Table S2. Minus RT values and standard curves are shown in Additional file 4: Table S3.
For all experiments, four measurements were done for two biological replicates. Data were normalized to the indicated timepoint/sample using the ∆∆Ct method .
HFF cells, PAA, and knockout viruses
HFF cells were cultured and infection with wild-type virus was performed as described previously . Cells were incubated with wild-type HSV-1 and FXE ICP0 knockout virus for 15’ and n12 ICP4 knockout virus (ICP4KO) for 60’. Where applicable, 300 μg/ml PAA (Sigma-Aldrich, catalog number 284270) was added in conditioned medium when the inoculum was removed. The n12 ICP4 knockout virus  was a kind gift from Neal DeLuca, and the FXE ICP0 knockout  virus was a kind gift from Stacey Efstathiou.
Digital gene expression assays/Nanostring
Gene expression counts in virally infected cells were determined using nCounter Elements TagSets (NanoString) and an nCounter SPRINT profiler. Probe A and B sequences were custom designed with support from NanoString; oligonucleotides were obtained from IDT (Additional file 5: Table S4). Hybridizations were carried out for 18 h. Raw data were exported directly from the SPRINT profiler and normalized using the slope of the linear regression for the positive control values as normalization factor. All normalized values are shown in Additional file 6: Table S5.
Chromatin immunoprecipitation (ChIP) was performed as previously described [59, 60] with the following modifications. Per sample, 40 μg of crosslinked DNA (as determined after reverse crosslinking, protease treatment, phenol extraction, and ethanol precipitation) was used as input. ICP4-3xflag crosslinked to chromatin was precipitated by incubation of sheared chromatin with anti-flag M1 agarose affinity gel (Sigma-Aldrich, catalog number A2220) for 3 h at 4 °C. Yields in the eluates were about 3 ng DNA for control samples (untransfected cells) and 10 ng for ICP4-ChIP samples, as determined by measurements with the Qubit DNA HS assay (Thermo Fisher, catalog number Q32851). For qPCR, 10 ng of DNA template for input and 50 pg for eluate samples was used, as described in the “RT-qPCR analysis” section above.
Single-molecule RNA fluorescence in situ hybridization
Oligonucleotide probes were designed using the Stellaris Probe Designer (https://www.biosearchtech.com/support/tools/design-software/stellaris-probe-designer). Customs probes listed in Additional file 7: Table S6 were ordered either coupled to Quasar 570 (Cy3 replacement) or Quasar 670 (Cy5 replacement) from LGC Biosearch technologies, as well as probes for NEAT1 (catalog number SMF-2036-1). Cells were seeded in 12-well format on gelatin-coated coverslips; 24 h later cells were infected with HSV-1 and fixed with 4% paraformaldehyde at 0, 4, and 6 hpi. Cells were washed once with PBS containing Ca and Mg and fixed in 4% EM-grade formaldehyde in PBS containing 1 mM Ca and Mg for 10 minutes at room temperature, transferred to ice, washed once with PBS, and stored up to several weeks in 70% ethanol at 4 °C.
Hybridizations were performed as described , with minor modifications. Stellaris probes were hybridized overnight at 37 °C in 2 × SSC, 10% dextran sulfate, 0.4% SDS with a stringency of 10% formamide at a concentration of 12.5 to 50 nM. After hybridization, unbound probes were washed away two times, and 5 ng per ml DAPI nuclear stain was included in the second wash. Images were acquired on an inverted Nikon Ti microscope with an Andor iXon Ultra DU-888 camera in CCD mode (1 MHz, conversion gain 2), a 60x NA 1.4 oil objective and Nikon NIS-Elements Ar software (version X), using exposure times of 50 ms for DAPI and 300 to 500 ms in the other fluorescence channels at 50 to 100% intensity of the light source (Sola SEII). Groups of cells for imaging were chosen in the DAPI channel; Z-stacks with 0.3 μm spacing were acquired in all channels including EGFP (to control for non-specific signals generated by autofluorescence). Images were processed using open source software Fiji and merged by maximum intensity projection, as indicated in the figure legends.
Promoter reporter assay
We used a modified pMIR-RNL TK vector (a dual firefly and renilla luciferase vector) described previously . The CMV promoter in this plasmid was replaced by the predicted BBC3 antisense transcript promoter. The sequences inserted into this novel construct pMIR-AS are shown in Additional file 1: Figure S5. The predicted promoter sequence (~700 bp) and shorter sequences (~500 and 300 bp) were cloned in the newly modified pMIR-RNLTK. Typically, HeLa cells were seeded in 24-well format and transfected with 0.2 μg per well reporter vector using PolyFect (Qiagen) and assays were conducted in duplicate. After 48 h cells were infected with HSV-1 or mock infected (0 h) with DMEM. Cells were harvested at 0, 6, and 8 hpi. The ratio of firefly (reporter)/renilla (control) luciferase for each sample was determined (%RLU).
Generation of BBC3as knockout cells
Optimal guide RNA sequences were determined using the CRISPR design website (crispr.mit.edu). The following gRNA sequences were used: TGCGCGGCTGGGTCGCTCCG (gRNA 5), GCTTTTACTCCGTATTTACC (gRNA 6), and TAGAAACACAAGAGCGCACC (gRNA 9). Templates for T7 transcription of candidate RNAs were prepared using oligos BBC3asprom5, BBC3asprom6, BBC3asprom9 as well as T7_tracrRNA (Additional file 3: Table S2) and the Phusion PCR polymerase (Thermo Fisher Scientific) to generate double-stranded DNA. In vitro transcription was performed using T7 polymerase (New England Biolabs, catalog number E2050S). For the in vitro cleavage assay using Cas9 (New England Biolabs, catalog number M0641S), DNA templates were amplified from HeLa genomic DNA using the following primers: BBC3asivtempl12345fw and BBC3asivtempl12345re for gRNA 5, BBC3asivtempl678fwd and BBC3asivtempl678rev for gRNA 6, and ivtemplC_fwd_3 and ivtemplC_rev_3 for gRNA 9. For expression of the gRNAs, we used plasmids pSpCas9(BB)-2A-GFP (pX458, Addgene plasmid 48138), a gift from Feng Zhang , and pU6-(BbsI)_CBh-Cas9-T2A-mCherry (px330-mCherry, Addgene plasmid 64324), a gift from Ralf Kuehn . To prepare gRNA/Cas9 expression plasmids, oligos BBC3aspromg5inpxo1/BBC3aspromg5inpxo2 were inserted into BbsI-digested px330-mCherry and BBC3aspromg6inpxo1/BBC3aspromg6inpxo2 or BBC3aspromg9inpxo1/BBC3aspromg9inpxo2, respectively, into BbsI-digested px458.
Plasmids containing gRNA 5 and gRNA 9 or gRNA 6 and gRNA 9, respectively, were co-transfected into HeLa cells and double transfectants selected and single cells sorted into 96-well plates using FACS. Cells were recovered for 2 to 3 weeks and, after replication of plates, washed with PBS, detached with PBS/1 mM EDTA, and lysed with 0.5% Triton X-100 and 2 mg/ml Proteinase K in 30 mM Tris pH 8. The mixture was incubated for 10 minutes at 70 °C and the Proteinase inactivated by incubation at 95 °C for 10 minutes. The lysate was then used for qPCR using primers BBC3asCRIS_wtonly_fwd2 and BBC3asCRIS_wtonly_rev2, which detect only the wild-type allele, and either BBC3asCRIS_56only_fwd1/BBC3asCRIS_56only_rev1 or BBC3asCRIS_59only_fwd1/BBC3asCRIS_59only_rev1, which detect only mutant alleles. Ploidy of selected clones was then measured using the same primer pairs with BBC3as-2_fwd/BBC3as-2_rev for normalization from purified genomic DNA. For the experiments here, the homozygous knockout clones 56-D2 and 59-B3 were used.
For transcription activity assays, nascent RNA was labeled for 20 minutes with 500 μM 4sU. After MTS-mediated biotinylation  from 80 μg of total RNA, labeled RNA was recovered as previously described .
Analysis of histone modification data
Histone modification profiles were plotted in a region around the promoter of antisense and control genes in the following way. Genes that had a transcription start site within 5 kb of an antisense (novel or annotated) transcript or had expression < 10 RPKM in all time points of the 4sU data were removed from the analysis. From the resulting set, a random set of 1000 genes was chosen.
Antisense transcripts were marked as induced if they had a minimal fold change of 2 in at least one time point of the 4sU data.
Histone modification ChIP-seq signals were extracted in a region of ± 4 kb around the TSS of both the antisense transcripts and selected control genes. The signal for each TSS was scaled to have zero mean and a standard deviation of one, and averaged separately for each of the following groups: divergent antisense transcripts, convergent antisense transcripts, internal antisense transcripts, control genes.
We acknowledge Ouidad Benlasfer and Charlotte Stein for excellent technical assistance, and thank Scott Lacadie for fruitful discussions.
We thank the Deutsche Forschungsgemeinschaft for funding.
Availability of data and materials
The datasets generated in this study are available in the NCBI GEO repository, accession number GSE97009. The following published RNA-sequencing datasets were analyzed for this study: HSV-1-infected fibroblasts, NCBI GEO accession number GSE59717 ; NET-seq data and nucleoplasmic RNA from EXOSC3 depletion in HeLa cells, NCBI GEO accession number GSE81662 . CMV infected fibroblasts, NCBI GEO accession number GSE69906 ; KSHV-infected iSLK cell data obtained via C. Arias ; VSV-infected keratinocytes, EBI ArrayExpress accession number E-MTAB-1717 ; LPS-treated monocytes, NCBI GEO accession number GSE60216 . Histone modification ChIP-seq data was dowloaded from http://hgdownload.cse.ucsc.edu/goldenPath/hg19/encodeDCC/wgEncodeBroadHistone/.
EW performed experiments, analyzed data and wrote the manuscript. JM performed most cell culture and virus work and reporter assays. VF contributed all bioinformatical analysis and wrote parts of the manuscript. CK designed and analyzed FISH and Nanostring assays and ran them together with AB, TH and AR performed the HFF cell culture. KT contributed to RT-qPCR and ChIP-qPCR, CF to ChIP-qPCR, LB to reporter assays, LK to molecular cloning. CCF supported data analysis. NR, AA, and LD, supervised parts of the project and FG and ML supervised the overall project. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Pellet PE, Roizman B. Herpesviridae. In: Knipe DM, Howley PM, editors. Fields Virology. 6th ed. Philadelphia: Wolters Kluwer Health/Lippincott Williams & Wilkins; 2013. p. 1802–2128.Google Scholar
- Rivas HG, Schmaling SK, Gaglia MM. Shutoff of host gene expression in influenza A virus and herpesviruses: similar mechanisms and common themes. Viruses. 2016;8(4):102.View ArticlePubMed CentralPubMedGoogle Scholar
- Dai-Ju JQ, Li L, Johnson LA, Sandri-Goldin RM. ICP27 interacts with the C-terminal domain of RNA polymerase II and facilitates its recruitment to herpes simplex virus 1 transcription sites, where it undergoes proteasomal degradation during infection. J Virol. 2006;80(7):3567–81.View ArticlePubMed CentralPubMedGoogle Scholar
- Abrisch RG, Eidem TM, Yakovchuk P, Kugel JF, Goodrich JA. Infection by herpes simplex virus 1 causes near-complete loss of RNA polymerase II occupancy on the host cell genome. J Virol. 2015;90(5):2503–13.View ArticlePubMedGoogle Scholar
- Sandri-Goldin RM. The many roles of the highly interactive HSV protein ICP27, a key regulator of infection. Future Microbiol. 2011;6(11):1261–77.View ArticlePubMedGoogle Scholar
- Read GS. Virus-encoded endonucleases: expected and novel functions. Wiley Interdiscip Rev RNA. 2013;4(6):693–708.View ArticlePubMedGoogle Scholar
- Sadek J, Read GS. The splicing history of an mRNA affects its level of translation and sensitivity to cleavage by the virion host shutoff endonuclease during herpes simplex virus infections. J Virol. 2016;90(23):10844–56.View ArticlePubMed CentralPubMedGoogle Scholar
- Rutkowski AJ, Erhard F, L'Hernault A, Bonfert T, Schilhabel M, Crump C, Rosenstiel P, Efstathiou S, Zimmer R, Friedel CC, et al. Widespread disruption of host transcription termination in HSV-1 infection. Nat Commun. 2015;6:7126.View ArticlePubMed CentralPubMedGoogle Scholar
- Jensen TH, Jacquier A, Libri D. Dealing with pervasive transcription. Mol Cell. 2013;52(4):473–84.View ArticlePubMedGoogle Scholar
- Cech TR, Steitz JA. The noncoding RNA revolution-trashing old rules to forge new ones. Cell. 2014;157(1):77–94.View ArticlePubMedGoogle Scholar
- Lacadie SA, Ibrahim MM, Gokhale SA, Ohler U. Divergent transcription and epigenetic directionality of human promoters. FEBS J. 2016;283(23):4214–22.View ArticlePubMedGoogle Scholar
- Grzechnik P, Tan-Wong SM, Proudfoot NJ. Terminate and make a loop: regulation of transcriptional directionality. Trends Biochem Sci. 2014;39(7):319–27.View ArticlePubMed CentralPubMedGoogle Scholar
- Villegas VE, Zaphiropoulos PG. Neighboring gene regulation by antisense long non-coding RNAs. Int J Mol Sci. 2015;16(2):3251–66.View ArticlePubMed CentralPubMedGoogle Scholar
- Khorkova O, Myers AJ, Hsiao J, Wahlestedt C. Natural antisense transcripts. Hum Mol Genet. 2014;23(R1):R54–63.View ArticlePubMed CentralPubMedGoogle Scholar
- Nishizawa M, Ikeya Y, Okumura T, Kimura T. Post-transcriptional inducible gene regulation by natural antisense RNA. Front Biosci (Landmark Ed). 2015;20:1–36.View ArticleGoogle Scholar
- Pelechano V, Steinmetz LM. Gene regulation by antisense transcription. Nat Rev Genet. 2013;14(12):880–93.View ArticlePubMedGoogle Scholar
- Saayman S, Ackley A, Turner AM, Famiglietti M, Bosque A, Clemson M, Planelles V, Morris KV. An HIV-encoded antisense long noncoding RNA epigenetically regulates viral transcription. Mol Ther. 2014;22(6):1164–75.View ArticlePubMed CentralPubMedGoogle Scholar
- Hon CC, Ramilowski JA, Harshbarger J, Bertin N, Rackham OJ, Gough J, Denisenko E, Schmeier S, Poulsen TM, Severin J, et al. An atlas of human long non-coding RNAs with accurate 5′ ends. Nature. 2017;543(7644):199–204.View ArticlePubMedGoogle Scholar
- Hu B, Huo Y, Chen G, Yang L, Wu D, Zhou J. Functional prediction of differentially expressed lncRNAs in HSV-1 infected human foreskin fibroblasts. Virol J. 2016;13:137.View ArticlePubMed CentralPubMedGoogle Scholar
- Stark TJ, Arnold JD, Spector DH, Yeo GW. High-resolution profiling and analysis of viral and host small RNAs during human cytomegalovirus infection. J Virol. 2012;86(1):226–35.View ArticlePubMed CentralPubMedGoogle Scholar
- Hargett D, McLean T, Bachenheimer SL. Herpes simplex virus ICP27 activation of stress kinases JNK and p38. J Virol. 2005;79(13):8348–60.View ArticlePubMed CentralPubMedGoogle Scholar
- Schlackow M, Nojima T, Gomes T, Dhir A, Carmo-Fonseca M, Proudfoot NJ. Distinctive patterns of transcription and RNA processing for human lincRNAs. Mol Cell. 2017;65(1):25–38.View ArticlePubMed CentralPubMedGoogle Scholar
- Lloret-Llinares M, Mapendano CK, Martlev LH, Lykke-Andersen S, Jensen TH. Relationships between PROMPT and gene expression. RNA Biol. 2016;13(1):6–14.View ArticlePubMedGoogle Scholar
- Ntini E, Jarvelin AI, Bornholdt J, Chen Y, Boyd M, Jorgensen M, Andersson R, Hoof I, Schein A, Andersen PR, et al. Polyadenylation site-induced decay of upstream transcripts enforces promoter directionality. Nat Struct Mol Biol. 2013;20(8):923–8.View ArticlePubMedGoogle Scholar
- Jones M, Dry IR, Frampton D, Singh M, Kanda RK, Yee MB, Kellam P, Hollinshead M, Kinchington PR, O’Toole EA, et al. RNA-seq analysis of host and viral gene expression highlights interaction between varicella zoster virus and keratinocyte differentiation. PLoS Pathog. 2014;10(1):e1003896.View ArticlePubMed CentralPubMedGoogle Scholar
- Tirosh O, Cohen Y, Shitrit A, Shani O, Le-Trilling VT, Trilling M, Friedlander G, Tanenbaum M, Stern-Ginossar N. The transcription and translation landscapes during human cytomegalovirus infection reveal novel host-pathogen interactions. PLoS Pathog. 2015;11(11):e1005288.View ArticlePubMed CentralPubMedGoogle Scholar
- Arias C, Weisburd B, Stern-Ginossar N, Mercier A, Madrid AS, Bellare P, Holdorf M, Weissman JS, Ganem D. KSHV 2.0: a comprehensive annotation of the Kaposi’s sarcoma-associated herpesvirus genome using next-generation sequencing reveals novel genomic and functional features. PLoS Pathog. 2014;10(1):e1003847.View ArticlePubMed CentralPubMedGoogle Scholar
- Uhlitz F, Sieber A, Wyler E, Fritsche-Guenther R, Meisig J, Landthaler M, Klinger B, Bluthgen N. An immediate-late gene expression module decodes ERK signal duration. Mol Syst Biol. 2017;13(5):928.View ArticlePubMed CentralPubMedGoogle Scholar
- Zachos G, Clements B, Conner J. Herpes simplex virus type 1 infection stimulates p38/c-Jun N-terminal mitogen-activated protein kinase pathways and activates transcription factor AP-1. J Biol Chem. 1999;274(8):5097–103.View ArticlePubMedGoogle Scholar
- Tallen G, Riabowol K. Keep-ING balance: tumor suppression by epigenetic regulation. FEBS Lett. 2014;588(16):2728–42.View ArticlePubMedGoogle Scholar
- Boutell C, Everett RD. Regulation of alphaherpesvirus infections by the ICP0 family of proteins. J Gen Virol. 2013;94(Pt 3):465–81.View ArticlePubMedGoogle Scholar
- Sampath P, Deluca NA. Binding of ICP4, TATA-binding protein, and RNA polymerase II to herpes simplex virus type 1 immediate-early, early, and late promoters in virus-infected cells. J Virol. 2008;82(5):2339–49.View ArticlePubMedGoogle Scholar
- Wagner LM, DeLuca NA. Temporal association of herpes simplex virus ICP4 with cellular complexes functioning at multiple steps in PolII transcription. PLoS One. 2013;8(10):e78242.View ArticlePubMed CentralPubMedGoogle Scholar
- Wang X, Diao C, Yang X, Yang Z, Liu M, Li X, Tang H. ICP4-induced miR-101 attenuates HSV-1 replication. Sci Rep. 2016;6:23205.View ArticlePubMed CentralPubMedGoogle Scholar
- Consortium EP. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489(7414):57–74.View ArticleGoogle Scholar
- Bornelov S, Komorowski J, Wadelius C. Different distribution of histone modifications in genes with unidirectional and bidirectional transcription and a role of CTCF and cohesin in directing transcription. BMC Genomics. 2015;16:300.View ArticlePubMed CentralPubMedGoogle Scholar
- Hikisz P, Kilianska ZM. PUMA, a critical mediator of cell death--one decade on from its discovery. Cell Mol Biol Lett. 2012;17(4):646–69.View ArticlePubMedGoogle Scholar
- Nguyen ML, Blaho JA. Cellular players in the herpes simplex virus dependent apoptosis balancing act. Viruses. 2009;1(3):965–78.View ArticlePubMed CentralPubMedGoogle Scholar
- Dunagin M, Cabili MN, Rinn J, Raj A. Visualization of lncRNA by single-molecule fluorescence in situ hybridization. Methods Mol Biol. 2015;1262:3–19.View ArticlePubMedGoogle Scholar
- Cabili MN, Dunagin MC, McClanahan PD, Biaesch A, Padovan-Merhar O, Regev A, Rinn JL, Raj A. Localization and abundance analysis of human lncRNAs at single-cell and single-molecule resolution. Genome Biol. 2015;16:20.View ArticlePubMed CentralPubMedGoogle Scholar
- Landry JJ, Pyl PT, Rausch T, Zichner T, Tekkedil MM, Stutz AM, Jauch A, Aiyar RS, Pau G, Delhomme N, et al. The genomic and transcriptomic landscape of a HeLa cell line. G3 (Bethesda). 2013;3(8):1213–24.View ArticleGoogle Scholar
- Dimitrova N, Zamudio JR, Jong RM, Soukup D, Resnick R, Sarma K, Ward AJ, Raj A, Lee JT, Sharp PA, et al. LincRNA-p21 activates p21 in cis to promote Polycomb target gene expression and to enforce the G1/S checkpoint. Mol Cell. 2014;54(5):777–90.View ArticlePubMed CentralPubMedGoogle Scholar
- Rabani M, Levin JZ, Fan L, Adiconis X, Raychowdhury R, Garber M, Gnirke A, Nusbaum C, Hacohen N, Friedman N, et al. Metabolic labeling of RNA uncovers principles of RNA production and degradation dynamics in mammalian cells. Nat Biotechnol. 2011;29(5):436–42.View ArticlePubMed CentralPubMedGoogle Scholar
- Marcinowski L, Lidschreiber M, Windhager L, Rieder M, Bosse JB, Radle B, Bonfert T, Gyory I, de Graaf M, Prazeres da Costa O, et al. Real-time transcriptional profiling of cellular and viral gene expression during lytic cytomegalovirus infection. PLoS Pathog. 2012;8(9):e1002908.View ArticlePubMed CentralPubMedGoogle Scholar
- Gomes NP, Espinosa JM. Gene-specific repression of the p53 target gene PUMA via intragenic CTCF-Cohesin binding. Genes Dev. 2010;24(10):1022–34.View ArticlePubMed CentralPubMedGoogle Scholar
- Lu J, Wu X, Hong M, Tobias P, Han J. A potential suppressive effect of natural antisense IL-1beta RNA on lipopolysaccharide-induced IL-1beta expression. J Immunol. 2013;190(12):6570–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Lissner MM, Thomas BJ, Wee K, Tong AJ, Kollmann TR, Smale ST. Age-related gene expression differences in monocytes from human neonates, young adults, and older adults. PLoS One. 2015;10(7):e0132061.View ArticlePubMed CentralPubMedGoogle Scholar
- Kulej K, Avgousti DC, Sidoli S, Herrmann C, Della Fera AN, Kim ET, Garcia BA, Weitzman MD. Time-resolved global and chromatin proteomics during herpes simplex virus type 1 (HSV-1) infection. Mol Cell Proteomics. 2017;16(4 suppl 1):S92–S107.View ArticlePubMedGoogle Scholar
- Tang S, Patel A, Krause PR. Herpes simplex virus ICP27 regulates alternative pre-mRNA polyadenylation and splicing in a sequence-dependent manner. Proc Natl Acad Sci U S A. 2016;113(43):12256–61.View ArticlePubMed CentralPubMedGoogle Scholar
- Rosa S, Duncan S, Dean C. Mutually exclusive sense-antisense transcription at FLC facilitates environmentally induced gene repression. Nat Commun. 2016;7:13031.View ArticlePubMed CentralPubMedGoogle Scholar
- Worsley Hunt R, Mathelier A, Del Peso L, Wasserman WW. Improving analysis of transcription factor binding sites within ChIP-Seq data based on topological motif enrichment. BMC Genomics. 2014;15:472.View ArticlePubMed CentralPubMedGoogle Scholar
- Andersson R. Promoter or enhancer, what’s the difference? Deconstruction of established distinctions and presentation of a unifying model. Bioessays. 2015;37(3):314–23.View ArticlePubMedGoogle Scholar
- Everett RD, Sourvinos G, Orr A. Recruitment of herpes simplex virus type 1 transcriptional regulatory protein ICP4 into foci juxtaposed to ND10 in live, infected cells. J Virol. 2003;77(6):3680–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP. Integrative genomics viewer. Nat Biotechnol. 2011;29(1):24–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.View ArticlePubMed CentralPubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402–8.View ArticlePubMedGoogle Scholar
- DeLuca NA, Schaffer PA. Physical and functional domains of the herpes simplex virus transcriptional regulatory protein ICP4. J Virol. 1988;62(3):732–43.PubMed CentralPubMedGoogle Scholar
- Everett RD. Construction and characterization of herpes simplex virus type 1 mutants with defined lesions in immediate early gene 1. J Gen Virol. 1989;70(Pt 5):1185–202.View ArticlePubMedGoogle Scholar
- Brookes E, de Santiago I, Hebenstreit D, Morris KJ, Carroll T, Xie SQ, Stock JK, Heidemann M, Eick D, Nozaki N, et al. Polycomb associates genome-wide with a specific RNA polymerase II variant, and regulates metabolic genes in ESCs. Cell Stem Cell. 2012;10(2):157–70.View ArticlePubMed CentralPubMedGoogle Scholar
- Stock JK, Giadrossi S, Casanova M, Brookes E, Vidal M, Koseki H, Brockdorff N, Fisher AG, Pombo A. Ring1-mediated ubiquitination of H2A restrains poised RNA polymerase II at bivalent genes in mouse ES cells. Nat Cell Biol. 2007;9(12):1428–35.View ArticlePubMedGoogle Scholar
- Raj A, van den Bogaard P, Rifkin SA, van Oudenaarden A, Tyagi S. Imaging individual mRNA molecules using multiple singly labeled probes. Nat Methods. 2008;5(10):877–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhu JY, Pfuhl T, Motsch N, Barth S, Nicholls J, Grasser F, Meister G. Identification of novel Epstein-Barr virus microRNA genes from nasopharyngeal carcinomas. J Virol. 2009;83(7):3333–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Ran FA, Hsu PD, Wright J, Agarwala V, Scott DA, Zhang F. Genome engineering using the CRISPR-Cas9 system. Nat Protoc. 2013;8(11):2281–308.View ArticlePubMed CentralPubMedGoogle Scholar
- Chu VT, Weber T, Wefers B, Wurst W, Sander S, Rajewsky K, Kuhn R. Increasing the efficiency of homology-directed repair for CRISPR-Cas9-induced precise gene editing in mammalian cells. Nat Biotechnol. 2015;33(5):543–8.View ArticlePubMedGoogle Scholar
- Duffy EE, Simon MD. Enriching s4 U-RNA using methane thiosulfonate (MTS) chemistry. Curr Protoc Chem Biol. 2016;8(4):234–50.View ArticlePubMedGoogle Scholar
- Rädle B, Rutkowski AJ, Ruzsics Z, Friedel CC, Koszinowski UH, Dölken L. Metabolic labeling of newly transcribed RNA for high resolution gene expression profiling of RNA synthesis, processing and decay in cell culture. J Vis Exp. 2013;(78):e50195. doi:10.3791/50195
- Adelson ME, Feola M, Trama J, Tilton RC, Mordechai E. Simultaneous detection of herpes simplex virus types 1 and 2 by real-time PCR and pyrosequencing. J Clin Virol. 2005;33(1):25–34.View ArticlePubMedGoogle Scholar