Cellular responses to hsv-1 infection are linked to specific types of alterations in the host transcriptome

ABSTRACT Pathogen invasion triggers a number of cellular responses and alters the host transcriptome. Here we report that the type of changes to cellular transcriptome is related to the type

of cellular functions affected by lytic infection of Herpes Simplex Virus type I in Human primary fibroblasts. Specifically, genes involved in stress responses and nuclear transport

exhibited mostly changes in alternative polyadenylation (APA), cell cycle genes showed mostly alternative splicing (AS) changes, while genes in neurogenesis, rarely underwent these changes.

Transcriptome wide, the infection resulted in 1,032 cases of AS, 161 incidences of APA, 1,827 events of isoform changes and up regulation of 596 genes and down regulations of 61 genes

compared to uninfected cells. Thus, these findings provided important and specific links between cellular responses to HSV-1 infection and the type of alterations to the host transcriptome,

highlighting important roles of RNA processing in virus-host interactions. SIMILAR CONTENT BEING VIEWED BY OTHERS INFLUENZA VIRUS REPURPOSES THE ANTIVIRAL PROTEIN IFIT2 TO PROMOTE

TRANSLATION OF VIRAL MRNAS Article 24 August 2020 THE RNA QUALITY CONTROL PATHWAY NONSENSE-MEDIATED MRNA DECAY TARGETS CELLULAR AND VIRAL RNAS TO RESTRICT KSHV Article Open access 03 July

2020 NON-AUG HIV-1 UORF TRANSLATION ELICITS SPECIFIC T CELL IMMUNE RESPONSE AND REGULATES VIRAL TRANSCRIPT EXPRESSION Article Open access 18 February 2025 INTRODUCTION Herpes Simplex Virus

type I (HSV-1) is a 152 kb double-strand DNA virus containing around 80 genes that infects about 80% of the human population1,2. Following a primary infection HSV-1 normally enters a latent

infection in sensory neurons of periphery sensory ganglions3. Stress signals and weakened immunity cause the reactivation of HSV-1 from sensory neurons and lytic infection ensues in

epithelial cells which these neurons innervate, resulting in cold sores or Herpes keratitis4. HSV-1 lytic infection in cultured cells unfolds rapidly, with the expression of immediate early

(IE) genes, including ICP0 and ICP4, the virus quickly recruits host RNA polymerase, transcription co-regulators and chromatin modifying complexes to facilitate viral early gene expression

and prepare for viral DNA synthesis followed by late gene expression, occurring at approximately 6 hours post infection (hpi). The incoming virus triggers a number of host responses

including the activation of the interferon pathway5,6, the DNA damage response7,8,9, apoptosis10,11 and other host defense mechanisms limiting viral growth12. In return, many of the viral

genes are designed to modulate these responses to ensure viral transcription, genome synthesis and assembly. ICP34.5, for example, is a key viral factor interfering with the interferon β

(IFN-β)pathway13,14,15. ICP27, ICP4 and ICP22, on the other hand, are negative regulators of the host apoptotic response, while ICP8 inhibits the host DNA damage response by inactivating the

ATR kinase16,17,18,19. ICP0 inhibits host transcription silencing activity by displacing host CoREST silencing complex, while it also degrades RNF8 and RNF168, two ubiquitin ligases in the

DDR pathway8,20,21 and components of the PML body22. The ICP27 factor also inhibits host RNA splicing23,24, while the viral host shutoff protein (vhs) degrades host, as well as viral

RNAs25,26,27,28,29. These viral factors, together with host responses, result in complex viral host interactions, the outcome of which determines whether the virus enters a lytic infection,

or becomes suppressed and enters latency. Many details of these processes are reflected in the alterations in the transcriptome of infected host cells. Thus how host cells respond to viral

infection at transcriptomic level is an important but under explored question. Studies in HSV-1 infected mouse embryonic fibroblast cells (MEF), mouse cornea and trigeminal ganglion have led

to the discovery of new genes and pathways of virus-host interactions14,30,31. However, these studies were done using DNA microarrays with limited genome coverage, thus likely missing many

important genes, furthermore these studies are unable to analyze other types of changes such as alternative splicing (AS), alternative polyadenylation (APA) and gene isoform composition.

Transcriptome-sequencing (RNA-seq) has revealed that approximately 94% of human genes are alternatively spliced, generating a much larger diversity of functional variants from a fixed number

of genes in the genome32. AS is widely known to participate in a wide range of biological processes including virus-host interactions. The HSV-1 encoded ICP27 protein, which exports

unspliced mRNA to the cytosol, is also reported to alter PML protein23 and glycoprotein C24 isoform composition via alternative splicing, while the SM factor from EBV changes isoform

expression of STAT for the benefit of EBV infection33. At genome wide level, profound changes in AS occur as cells adapt to external stimulus. For example, dendritic cells produced wide

spread changes in alternative splicing when encountering bacteria34. Also, when compared to normal cells, genome wide changes in splice isoforms were also observed in carcinoma35. Although

many splicing factors are described, the mechanisms regulating the splicing process in response to external signals are not well understood. The CTCF factor is one of the few known host

regulators reported to interact near many intron-exon junctions and alters the rate of RNA Pol II elongation and favors the splicing of downstream exons. CTCF binding sites are subject to

imprinting, thus providing means of epigenetic regulation of alternative splicing36. In addition to AS, more than 50% of human genes are also subject to APA37, which changes in response to

different physiological conditions or during cell differentiation. For example, increased proliferation, dedifferentiation and disease conditions are associated with proximal polyadenylation

sites, such as C2C12 myoblast compared to C2C12 differentiated myotubes38, resting B cell compared to activated B cells39, or MCF7 breast cancer cells compared to normal breast epithelial

cells40. In contrast, when cells exit the cell cycle and become differentiated, distal polyadenylation sites are preferred. Terminally differentiated neurons41 and activated B cells are

examples39. The longer 3′ UTRs tend to contain multiple miRNA targets42 and other cis elements43, providing fine tuned regulation of gene transcript levels. Although the mechanism is not

well understood, the choice of APA sites is affected by the interaction between the C-terminal domain of the largest subunit RNA polymerase II and a processing complex CPSF (cleavage and

polyadenylation specificity factor), CstF (cleavage stimulation factor) and the canonical poly(A) signal AAUAAA44. It is also affected by the usage of alternative promoters and by epigenetic

modification of local chromatin44. To date there are no reported studies analyzing whether and how viral infection affects host cell RNA splicing or polyadenylation site choices at the

transcriptomic level. Here we report an RNA-seq analysis of the HSV-1 infected host transcriptome of human primary fibroblast BJ cells at 6hpi and found that the infection up regulated 596

genes and down regulated 61 genes, significantly more than previous reports. This analysis identified many new important viral induced genes including epigenetic regulators BRD2 and CBX4,

regulator of RNA metabolism such as ZNF36. But, importantly, the infection resulted in 1,032 alternative splicing events, 161 APA and 1,827 cases of gene isoform expression changes. GO

analyses revealed that stress response, cell cycle and nuclear transport genes are almost exclusively regulated by AS or APA, genes related to neurogenesis are mostly regulated by

differential expression. This pattern suggests important regulatory and mechanistic differences among host responses and highlights important roles of RNA processing in virus-host

interactions. RESULTS Human BJ Skin Fibroblasts cells were used to study human HSV-1 lytic infection. The infection was done with the 17+ strain of HSV-1 at a multiplicity of infection (MOI)

of 5 to ensure close to 100% of the cells become infected and the infection is more synchronized among different cells. The infected cells were processed at 6 hpi and three independent

biological repeat samples were pooled for library construction and sequencing. 6 hpi samples reflect an important time point, when viral transcription and replication is active and the host

nucleus and host chromatin are still mostly intact, but most genes that are up regulated by HSV-1 infection could be readily detected. HSV-1 INFECTION CAUSED DIFFERENTIAL GENE EXPRESSION IN

HUMAN BJ CELLS As the differential expression analysis from RNA-seq data allowed the identification of important new candidate genes, including coding and noncoding genes mediating

virus-host interactions. The RNA-seq data consisted of 26 million (26,693,203) 90 bp paired-end reads for infected and similar amount for uninfected control (26,797,247), which was processed

by TopHat45,46 and Cufflinks47,48,49,50. Approximately 84% of reads from the control group and 40% from infected sample were aligned to the human genome. The infected sample contained a

large amount of paired-end reads aligned to the HSV-1 genome (30%), thus, at 6hpi and 5 MOI, viral RNAs comprise about 40% of the total mapped RNA reads in the infected human BJ cells.

RNA-seq result showed that the infection led to increased expression of 596 annotated genes and down-regulation of just 61. Differentially expressed genes are plotted in Fig. 1a. The number

of down-regulated genes accounts for only 9.28% of differentially expressed genes, a ratio similar to that of a previous report in infected MEF cells14. Based on gene ontology (GO) analysis

differentially expressed genes performed by DAVID51, top 20 GO terms from these genes ranked by _p_ value are listed in Fig. 1b. Notably, these GO terms could be grouped into three broad

categories based on their function, with transcriptional regulation comprised of 6 GO terms (positive regulation of transcription from RNA polymerase II promoter, positive regulation of

transcription, regulation of transcription, positive regulation of gene expression, regulation of transcription from RNA polymerase II promoter and regulation of transcription,

DNA-dependent), 131 genes, metabolic (mostly nucleic acid, RNA and DNA metabolic) processes has 6 GO terms (positive regulation of biosynthetic process, positive regulation of nitrogen

compound metabolic process, positive regulation of cellular biosynthetic process, positive regulation of nucleobase, nucleoside, nucleotide and nucleic acid metabolic process, positive

regulation of macromolecule biosynthetic process and regulation of RNA metabolic process), 120 genes and developmental genes has 8 GO terms (sensory organ development, embryonic skeletal

system development, embryonic organ development, regionalization, cell fate commitment, embryonic morphogenesis, pattern specification process and neuron differentiation), 101 genes. These

results suggest that transcriptional regulation, metabolic processes and developmental genes are the three most significantly affected classes of genes. Viral modification of cellular

metabolic processes, especially nucleic acid metabolism is crucial for viral transcription, DNA synthesis and alternation of host RNA processing. In contrast, the biological function of

activating neurogenic genes by HSV-1 infection is not clear, but this is likely due to the viral ICP0 protein, which displaces the host CoREST repressor complex from silencing viral

transcription, in this process ICP0 also derepresses many host neurogenic genes 8,20,21,52. As shown in the heatmap in Fig. 1c, we extracted the four categories of genes known to be affected

by HSV-1 infection from GO analyses result: regulation of transcription (123 genes), neurogenesis (49 genes), apoptosis (46 genes) and immune system development (17 genes) to group the

differentially expressed genes (listed in Supplementary Table S1). The infection affected genes in these categories are similar to a previous analyses in MEF cells14, but containing many

newly identified genes (listed in Supplementary Table S2), suggesting that the RNA-seq used in the present analysis revealed changes in the virus infected transcriptome in much greater

details. Transcriptional regulation may play a key role in driving many cellular processes in response to viral infection. For example, two key epigenetic regulators, BRD253,54 and

CBX455,56,57 are highly induced genes (Fig. 1d,e), possibly to reprogram the cellular epigenetic landscape in respond to the infection, or to facilitate viral transcription. Apoptosis is

known to be induced and modified by HSV infection14. BCL2L11, an important factor in apoptosis pathway58,59,60, was up regulated 103 fold (Fig. 1d,e), suggesting increased apoptotic

signaling as a result of cellular stress. TP73 is another pro-apoptotic gene activated by the DNA damage response61. It was up regulated by 81-fold. In contrast, as the virus takes over the

control of the host cells, it also inhibits the apoptotic process to allow the virus to reproduce, or to establish latent infection. For example, PIM3, a Serine/Threonine kinase family of

proto-oncogene with an anti-apoptotic property62,63, was up regulated 8-fold by HSV-1 infection (Supplementary Table S1). HSV-1 infection is known to induce DNA damage response (DDR) and

recruits DNA repair factors to the viral replication centers64. Here we found three genes H2AFX (Fig. 1d, e), XAB2 (Supplementary Fig. S1) and PAPD7 were up regulated by the infection. H2AFX

is the precursor of modified histone γH2A.X, a key chromatin mark of DNA damage response (DDR)65,66, XAB2 is involved in transcription-coupled repair (TCR)67,68, while PAPD7 is a DNA

polymerase involved in DNA repair and sister chromatid cohesion69. The induction of these genes suggests that HSV-1 infection activated DDR may also lead to gene expression level changes in

DDR associated factors. Like all viral infections, host inflammatory reaction (cytokines and chemokines) and intrinsic antiviral response genes (the interferon pathway) are activated. A

table of HSV-1 infection induced inflammatory genes and interferon β (IFN β) regulated genes are listed in Supplementary Table S1, Table S2. During infection, HSV-1 also modulates the host

interferon response through ICP34.5 and many commonly activated immunity genes are attenuated over time14. Several examples of newly identified innate immunity genes include DDX58 (DEAD Box

Polypeptide 58), OAS1 (2′-5′-Oligoadenylate Synthetase 1) and TLR4 (Toll-Like Receptor 4) (Fig. 1d,e). DDX58 (also known as RIG-I), a cellular viral sensor70,71, was activated 6.6 fold by

the infection. OAS1, an interferon activated 2-5A synthetase72,73, which synthesizes 2′,5′-oligoadenylates to activate RNase L to degrade viral (as well as cellular) RNAs, was activated 2.2

fold by HSV-1 infection. In contrast, TLR4, a member of the Toll-Like Receptor family of cellular membrane viral sensors74,75, was down regulated 2.6 fold by the infection. These examples

highlighted the activation or modulation of cellular antiviral response by the incoming virus. In addition to protein coding genes, we also discovered 51 Long noncoding RNAs which were

significantly regulated after HSV-1 infection, including H19 and MEG9 (Fig. 1d,e). We subsequently selected several newly identified host genes that are either up or down regulated following

infection for qRT-PCR validation. qRT-PCR results exhibited the similar gene expression level changes as the RNA-seq analysis, with a correlation coefficient, _r_ of 0.85 (_p_ value < 

0.01, Pearson test), thus validating the analysis (Fig. 1d,e, see primer information in Supplementary Table S3). HSV-1 INFECTION INDUCED HOST ALTERNATIVE SPLICING Alternative splicing (AS)

events32,76 include skipped exon (SE), retained intron (RI), alternative to 5′ splicing site (AS5), alternative to 3′ splicing site (AS3), mutually exclusive exon (MXE), alternative start

(altstart), alternative end (altend) and skip multiple exons (skip_multi_exon) (see Fig. 2a for a graphic explanation). We used ASD software76 to analyze the RNA-seq data and found 1,032

significant AS events occurred in 883 genes after HSV-1 infection. While most AS events (771) happened once per gene, they occurred twice in 80 genes, three times in 25 genes, four times in

4 genes (ARNTL, ATXN2L, MAP3K12 and U2AF1L4) and five times in 2 genes (B4GALT4 and FHOD1). The occurrence of different types of AS is summarized in Fig. 2b. Notably, the SE type of AS (274

cases) comprising the largest share, over a quarter, of all AS events. The second largest group, 21.7% of all AS events, is intron retention (RI). To confirm this analysis, we validated 4 AS

events in 4 genes, HCFC1R1 (host cell factor C1 regulator 1), STK11 (Serine/Threonine Kinase 11), c-Fos and NUFIP2 (Nuclear fragile X mental retardation-interacting protein 2), using

primers designed to detect the changes (Supplementary Table S3). In Fig. 2c,d, both HCFC1R1 and STK11 showed more retained introns after infection, as RT-PCR using primers specifically

amplifying the introns in question detected an increase in amount of fragment size of 358 bp and 742 bp for these two genes, respectively, confirming the increase of unspliced intron (arrows

in Fig. 2c,d). In contrast, c-Fos, an important immediate early cellular gene, which senses stress77, showed increased splicing after infection (Fig. 2e). Normally, the c-Fos transcript is

highly unstable with a half-life of only 15 minutes. Its RNA is only partially spliced with the unspliced precursor RNA quickly degraded due to a destabilizing element present in the third

intron77,78. Stress signals promote the splicing and produce stable c-Fos mRNA. From the RNA-seq result, the third intron became completely spliced after the infection, along with a

significant increase of the more stable c-Fos mRNA. Using primers that could distinguish spliced (686 bp: arrow in Fig. 2e) and unspliced intron number 3 (116 bp) we detected two fragments

of 802 bp and 686 bp with similar intensity before infection, but detected mostly the 686 bp following infection (Fig. 2e). This result confirmed the increased splicing of intron number 3 in

the c-Fos gene, suggesting that HSV-1 infection activated the cellular stress response by increasing the splicing of c-Fos pre mRNA. An example of SE type of alternative splicing is seen in

the NUFIP2 gene. To verify this, we used primers specific to detect the skipped exon and the RT-PCR result showed a drop in the amount of 424 bp fragment after infection, indicating exon

skipping (Fig. 2f). HSV-1 INFECTION LED TO CHANGES IN POLYADENYLATION IN THE HOST TRANSCRIPTOME The 3′ UTRs of genes often contain miRNA targets42 and other cis elements43, which could

modulate the turnover of mRNA. For genes with multiple polyadenylation sites, the selection of a proximal or a distal site could be an important regulatory mechanism to control gene

expression. Thus we analyzed the changes in host APA pattern using DaPars79 software to identify de novo dynamic APAs from RNA-seq data and found a total of 161 significant cases of APA

changes, with 61 cases displaying a shift from proximal polyadenylation sites to distal ones, resulting in longer 3′ UTRs and 100 cases showed switching from distal polyadenylation sites to

proximal sites, leading to shorter 3′ UTRs (Fig. 3a, Supplementary Table S4). As shortening of the 3′UTRs is masked to a certain degree by existing, longer UTRs from the same genes prior to

the infection, the 100 cases are likely an underestimate of the actual 3′ UTR shortening events. Thus, the overall changes of APA pattern are biased toward shortening of the 3′UTR, which,

together with the increase of skipped exons (Fig. 2b), suggests that HSV-1 infection led to an overall shortening of affected gene transcripts. We then chose 6 genes that exhibited

shortening of 3′ UTR after the infection and conducted qRT-PCR and RT-PCR (primer information in Supplementary Table S3) experiments to validate the analysis. SREK180 is a RNA processing

factor, CSRNP281 and TFAP2A82 are transcription activators, PAPD769,83 and NIT1 are DNA damage response factor84 and RUNX1 is a cell immunity factor85. As shown in Fig. 3b, qRT-PCR

validation of the 3′UTR of these genes showed reduction of regions of the 3′UTRs specific to the longer form, changes consistent with the APA analysis using the DaPars software79. The RT-PCR

validation of these genes is shown in Fig. 3c. Here, the PAPD7 gene showed an increase of reads in the last exon, but a drop in the 3′ UTR (boxed area), indicating a shortening. RT-PCR

result showed a reduction of a predicted 1,540bp product. Similar results were seen in the remainder 5 genes, while HSV-1 infection slightly reduced the reads of the exons of these genes,

the 3′ UTR regions were significantly and specifically reduced compared to the rest of these genes (see Fig. 3c boxed areas). Thus, both qRT-PCR and RT-PCR analyses showed changes consistent

with RNA-seq analysis, confirming the increased usage of proximal APA sites. In contrast, two examples of the lengthening of 3′UTR are shown in Fig. 3d. Here the 3′ reads corresponding to

the longer form of the UTR of both ABCC5 and IGF2 genes showed increases (see arrows in Fig. 3d and qRT-PCR validation in Fig. 3b). HSV-1 INFECTION INDUCED WIDE SPREAD CHANGES IN ISOFORM

COMPOSITION IN THE HOST TRANSCRIPTOME Isoform composition is an important aspect of the cellular transcriptome and reflects different cellular states and cell identity. Using TopHat45,46 and

Cufflinks47,48,49,50 softwares, we found 1,239 cases of increase in isoform expression and 588 cases of descrease in isoform expression after HSV-1 infection (plotted in Fig. 4a). These

changes occurred in 1,674 genes, with 1,544 genes containing changes in one isoform, 114 genes showing changes in two isoforms, 11 genes with three isoform changes, 4 genes with 4 isoform

changes and 1 gene with changes in 6 isoforms. Examples of isoform-specific expression changes are shown in Fig. 4b–f. SNHG9, a long intergenic noncoding RNA, has 2 isoforms, 00087626 and

00087627, the infection resulted in the significant expression of isoform 00087627 (Fig. 4b,d, note the boxed area). As the 00087626 isoform is a spliced form of 00087627, this infection

induced isoform switch is also a result of retained intron. The NUFIP2 gene contains 4 isoforms, with the infection reduced 00110828 specific reads (Fig. 4c, red boxed area). KLHL21 has 9

isoforms, HSV-1 infection increased the expression of transcript 00012543, which could be seen from the vertical boxed area where reads specific to isoform 00012543 showed a robust increase

(Fig. 4b,e). Finally, DGAT1 also has 11 isoforms and HSV-1 infection induced the expression of 00241783, as evidenced by the appearance of reads inside the red vertical box (Fig. 4b,f). To

determine how changes in isoform composition is related to that in differentiall gene expression, we compared 1,255 annotated genes that showed isoform composition changes with 657

annotated, differentially expressed genes and found 294 cases (or 23.43%) overlap, indicating a significantly large propotion of isoform composition changes (76.57%) did not result in gene

expression changes (Fig. 5a). In most of the 294 cases examined, the changes of isoform level were in the same direction of the changes of gene expression, however, we found 12 cases where

isoform expression level changes happened in the opposite directions with gene expression level changes, suggesting that these genes were subject to mutiple levels of regulation. This result

revealed that gene isoform changes were more prevalent than differential gene expression and in most cases, these changes did not lead to significant changes in overall gene expression

levels in viral infected host transcriptome. When genes with changes in AS and isoform compositon were compared, we found only 110 genes with AS changes also had changes in isoform

compositions (Fig. 5a). This is far fewer than expected, as alternative splicing is believed to be a major cause of isoform composition change. However, similar results have been observed by

others86,87, suggesting that other mechanisms, for example, alternative promoter usage or RNA turnover, may contribute more to changes in the composition of gene isoforms. We also compared

genes with isoform changes and genes with APA changes and found only 20 cases of overlap, a small fraction of the total isoform changes, as APA is not significantly related to gene isoform

composition changes (Fig. 5a). Finally, comparison of AS with APA genes revealed 22 cases of overlap, suggesting again that APA is an independently regulated event and does not contribute

significantly to AS (Fig. 5a). GO ANALYSIS LINKS CELLULAR RESPONSES TO VIRAL INFECTION TO THE TYPES OF ALTERATION IN THE HOST TRANSCRIPTOME As described earlier in Fig. 1b, we performed GO

analysis to reveal the type of biological processes exhibiting changes in AS, APA and isoform expression after HSV-1 infection. The top 20 significant GO terms in differential gene

expression, AS, APA and isoform changes are listed in Fig. 1b and Fig. 5b–d. The GO terms for alternative splicing changes are concentrated in cell cycle (4 terms, 67 genes), metabolic and

catabolic processes (12 terms, 187 genes) and transcription (2 GO terms, 40 genes) (Fig. 5b, Supplementary Table S5), suggesting that AS is an important mechanism regulating these processes,

which may quickly generate functionally different variants in response to HSV-1 infection. In the APA GO terms (Fig. 5c and Supplementary Table S4), chromatin and transcription regulation

occupy the largest shares with 8 GO terms, 32 genes, among these most displayed 3′UTR shortening. The remainder, nuclear transport has 2 terms, 13 genes, metabolic process has 4 terms, 13

genes as did stress responses with 13 genes, suggesting that processing the 3′ UTR played an important role in these cellular processes in response to viral infection. In the isoform change

category (Fig. 5d), 12 of the 20 GO terms including 235 genes are involved in metabolic, especially nucleic acid metabolic pathways, while 7 GO terms with 247 genes are involved in

transcription, suggesting that isoform expression switch is a key regulatory mechanism for these two processes in the host cells after viral infection. Finally, as described earlier (Fig.

1b) in differential gene expression analysis, we found that the largest group of GO terms belongs to neurogenesis and development process (8 GO terms, 101 genes), but these terms are notably

missing from AS, APA and isoform analyses, again suggesting that these genes are primarily regulated by changing the amount of transcript, not at the level of splicing. We next combined the

similar GO terms into six main cellular processes, 1) transcription and chromatin, 2) cell cycle, 3) metabolic processes, 4) nuclear transport, 5) neurogenesis and 6) stress response and

compared the changes in AS, APA, isoform composition and differential gene expression in each of these processes (Table 1). We found that genes in transcription and chromatin and metabolic

process employ all 4 types of regulation. However, the genes in the remaining GO terms are affected mostly through one specific type. For example, HSV-1 infection induced cell cycle genes

underwent only in AS changes (see a list in Supplementary Table S5), while nuclear transport genes and stress responses genes were mostly regulated by APA (also see Supplementary Table S4).

Two main host responses, immunity and apoptosis are also significant but are not within the top 20 for each type of transcriptomic changes. We found 25 immunity genes underwent isoform

switch and 17 genes exhibited differential gene expression after viral infection (Supplementary Table S6). None of immunity genes has AS or APA changes. In contrast, apoptotic genes have

undergone all four types of changes with 81 genes having changes in gene isoform compositions, 30 in differential gene expression, 52 in AS and 16 in APA, respectively (see list in

Supplementary Table S6). Taken together, these results strongly suggest mechanistic links between viral infection activated cellular pathways and specific types of alterations in the host

transcriptome (Fig. 6). DISCUSSION HSV-1 lytic infection triggers a cascade of cellular events and responses, its effects on cellular gene expression and function are the product of complex

virus-host interactions. Although many details of HSV productive infection are known, our understanding of its effects on the host transcriptome is very limited. In an attempt to obtain

insights into the mechanism of virus-host interactions, Kent _et al_.30 examined latently infected mouse trigeminal ganglia by microarray and found that HSV-1 induced 56 differentially

expressed genes, including neuronal-specific genes and several signaling molecules that could promote the initiation of HSV-1 infection. Pasieka _et al_.14 detected the change of host cell

gene expression also through microarray in HSV-1 infected mouse MEF cells and found 347 up regulated and 19 down regulated genes. Pasieka _et al_.31 performed array analysis of HSV-1

infected mouse corneas and observed roles for viral vhs and the host STAT signaling pathway in viral host interactions. These analyses, while providing valuable and important mechanistic

clues to host responses and virus-host interactions, are of limited scope due to the coverage of DNA arrays used. To gain a more comprehensive picture of the host transcriptomic changes in

HSV-1 infected human cells, we performed RNA-seq analysis of infected human primary fibroblast BJ cells and found a significantly larger number of differentially expressed genes in infected

cells when compared to previous DNA array studies. The newly identified genes include epigenetic regulators, components of the DDR, regulators of apoptotic pathway and genes mediating the

immune response. For example, CBX4 and BRD2 are previously unreported genes that were highly up regulated, suggesting that these epigenetic regulators may play important roles in the viral

infected transcriptome. Overall, the RNA-seq analysis revealed many more cases of increase than decrease of gene expression following HSV-1 infection, an observation similar to Pasieka _et

al_.14, which is likely due to the fact that existing RNAs prior to the infection masked the decrease. In a recent study, Rutkowski _et al_., analyzed nascent transcripts in HSV-1 infected

human epithelial cells and found that the infection induced wide spread read through of RNA polymerase resulting in fortuitous expression of hundreds of downstream genes88. Thus, to a

certain degree, a portion of the infection-induced transcriptome changes noted in our RNA-seq analyses of thousands of genes may not serve relevant biological functions. Importantly, our

RNA-seq analysis reveals detailed changes in host RNA splicing, 3′UTR length and gene isoform composition after the infection. Together we found 1,032 cases of AS changes, encompassing

skipped exons, retained introns, alternative to 5′ or 3′ splicing site, mutually exclusive exon, alternative start, alternative end and skip multiple exons (Fig. 2b). We suspect that many of

the SE cases observed here were due to HSV-1 infection induced activation of DNA damage response, as skipped exons were seen in cells under genotoxic stress77. Indeed, many of the cell

cycle genes underwent AS are of the SE type, CENPE and ZC3HC1 are two examples. A notable case highlighting the importance of AS is the c-Fos gene, which becomes more spliced as a result of

infection, producing a more stable mRNA and thus serves as a mechanism of activating the stress response77, How viral infection led to changes of AS is an important yet little understood

question. The viral factor ICP27 is known to disrupt host splicing, resulting in more unspliced mRNAs23,24. Indeed, we found 224 cases, or 21.7% of the changes in AS belong to retained

introns and presumably, many of the increases in intron retention were due to the activity of ICP27. However, the remaining AS events are most likely the result of regulatory activities

other than that of ICP27. For example, our analysis detected 304 cases (29.5%) of skipped exons or multiple skipped exons following HSV-1 infection (Fig. 2b). Part of these is likely due to

an HSV-1 infection induced DDR, which is known to cause skipped exons89. HSV-1 infection also elicited a change in the 3′ UTR through APA. In a total of 161 cases, we found 100 resulted in

shortening, suggesting that the majority of affected genes shifted towards making shorter transcripts under viral infection. For example, two stress response genes, SSRP1, PCBP4, (and 6

other validated examples of proximal polyadenylation site usage (Fig. 3b,c) after the infection produced shorter transcripts. In addition to AS and APA, the HSV-1 infected cell transcriptome

showed profound changes in gene isoform composition, including 1,239 cases of increases in isoform expression and 588 cases of descrease in isoform expression occurring in 1,674 genes (Fig.

4a). Taken together, the infection resulted in 657 differential expressed genes, but 2,149 genes with AS, APA and isoform changes combined (Fig. 5a), suggesting that the regulation of the

host transcriptome at the level of RNA processing is more extensive and plays at least as important a role in virus-host interactions, as regulation at the transcription level. More

importantly, our analysis revealed a pattern suggesting that HSV-1 induced cellular responses or pathways that are linked to specific types of changes to the host transcriptome (Table 1 and

Fig. 6). For example, cell cycle related genes are almost exclusively regulated at the level of AS after infection (listed in Supplementary Table S5), stress response genes are regulated at

the level of APA (also see a list in Supplementary Table S4), neural specific genes are regulated mostly by differential expression, while immunity genes are regulated at the level of

differential expression and isoform composition changes and nuclear transport genes are regulated both by AS and APA. In contrast, genes involved in transcription and chromatin and genes in

metabolic, mostly nucleic acid, RNA and DNA metabolic processes displayed all four types of changes after the infection. Although, the finding that viral ICP0 protein induced derepression of

neuronal genes and genotoxic stress induced changes of AS in cell cycle genes offer some clues as to why cellular responses to viral infection are linked to specific patterns of

transcriptome changes, the underlying mechanisms are likely complex and are mostly unknown. Here we observed that genes involved in RNA metabolism underwent the most profound changes, which

include differential gene expression, AS, APA and isoform composition changes. From differentially expressed genes alone, we recovered 108 genes involved in RNA metabolism, which may offer

candidates connecting cellular pathways to transcriptomic changes. Particularly, we found up-regulation of several genes that are specifically involved in splicing regulation (Supplementary

Fig. S1). For example, PABPC1 is known to participate in AS and APA regulation90. YBX1 mediates pre-mRNA alternative splicing regulation91,92. XAB2 is involved in transcription-coupled

repair (TCR), transcription and pre-mRNA splicing67,68, while, ZFP3693,94 is also an RNA splicing regulator. Thus, this analysis offers an interesting opportunity to further investigate

mechanism of HSV-1 infection induced regulation of RNA processing and host transcriptomic changes. Future studies in this direction should shed light on how cellular processes affect the

transcriptome. MATERIALS AND METHODS The methods were carried out in accordance with the approved guidelines. CELLS AND VIRUS Human BJ Skin Fibroblasts cells (ATCC CRL-2522) were purchased

from ATCC (Manassas, VA) and grown in complete Dulbecco’s modified Eagle’s medium (DMEM, Life technology Gibco, USA) containing 10% fetal bovine serum (FBS, Gibco, USA), 1%

penicillin/streptomycin (Gibco, USA) in a humidified 5% CO2 atmosphere at 37 °C. HSV-1 stain 17+ was used in this study. Virus was grown and titrated on Vero cells (ATCC, USA). Viral

infections were done according to standard protocols. Briefly, cultured cells were replaced with serum free DMEM, followed by adding the virus and incubating for 1 hour with occasional

rotation to get an even spread, then the culture medium was replaced by regular DMEM with 10% FBS and 1% antibiotics. All experiments were carried out in accordance with the approved

guidelines of ethics committee of Kunming Institute of Zoology and all experimental protocol were approved by ethics committee of Kunming Institute of Zoology, Chinese Academy of Sciences.

RNA EXTRACTION AND SEQUENCING BJ cells were infected with HSV-1 at an MOI of 5 and harvested at 6 hours post infection. The total RNA was extracted with TRIzol reagent (Ambion, 15596-018)

following the manufacturer’s instructions. As HSV-1 infection is a rapid process, sequencing individual infections would lead to high level of variation. To circumvent this issue, we pooled

the RNA samples from three biological repeats. After the DNase I treatment, magnetic beads with Oligo (dT) were used to isolate mRNA and then mixed with a fragmentation buffer to degrade the

mRNA into shorter fragments. Next, cDNA is synthesized using the mRNA fragments as templates. Short fragments are purified and resolved with EB buffer for sticky ends repair and single

nucleotide A addition. After that, the short fragments are connected with adapters. After agarose gel electrophoresis, the fragments averaging 200 bp) are selected for the PCR amplification

as templates. During the quality control steps, an Agilent 2100 Bioanaylzer and ABI StepOnePlus Real-Time PCR System were used in quantification and qualification of the sample library.

Finally, the library was sequenced using Illumina HiSeqTM 2000. RT-PCR VALIDATION BJ cells were infected with HSV-1 at an MOI of 5 and harvested at 6 hours post infection. The total RNA was

extracted with TRIzol reagent (Ambion, 15596-018) following the manufacturer’s instructions. One microgram of RNA was subjected to a DNase treatment with RQ1 RNase-Free DNase (Promega). cDNA

derived from this RNA was synthesized using RevertAid H Minus First Strand cDNA Synthesis Kit and Random Hexamer Primer. The synthesized cDNAs were amplified with TaKaRa PCR Amplification

Kit. Primers used to amplify each gene are listed in supplementary Table S3b. The products were analyzed using gel electrophoresis. RT-QPCR VALIDATION 1 μg RNA was reverse transcribed using

a primeScript RT and DNA Eraser (TaKaRa, DRR047A) Reagent Kits and stored at −80 °C. Real time PCR was run in triplicate with 50 ng cDNA using FastStart Universal SYBR Green Master (Roche,

04913914001) and ABI7900HT. Relative differences were determined using the ΔΔCt approach. ΔΔCt = (Ctinfection−Ct18S rRNA)−(Ctunfection−Ct18S rRNA). The fold enrichment value is 2−ΔΔCt. DATA

ANALYSIS FOR RNA-SEQ DATA: GENE AND ISOFORM DIFFERENTIAL EXPRESSION ANALYSIS Clean paired-end RNA-seq reads with 90 bp in each end were aligned to the human genome (Homo sapiens (release

37.72), Homo_sapiens.GRCh37.72.gtf) using the TopHat program (v2.0.9) with default parameters. Cufflinks (v2.1.1) was used to calculate Gene and Isoform expression level. We used cutoff

(_p__Value≤ 0.05, FPKM (fragments per kilobaseof exon per million fragments mapped) ≥1 at least one group; FPKM fold ≥2) to select significant, differentially expressed genes and isoforms at

6hpi. We used Tophat (v2.0.9) to align the remainder of unmapped paired-end reads to the HSV-1 strain 17 plus genome (GenBank: JN555585.1). ALTERNATIVE SPLICING ANALYSIS The

“accepted_hits.bam” files generated by Tophat (v2.0.9) were used to detect alternative splicing using the Java program, ASD (AS detector) (v1.2) with the annotation file

Homo_sapiens.GRCh37.72.gtf. We selected significant alternatively splicing cases with adjusted__p__Value ≤ 0.05. ALTERNATIVE POLYADENYLATION ANALYSIS Clean paired-end RNA-seq reads were

aligned to the human genome (Hg19) using TopHat (v2.0.9). We used DaPars (Dynamic analysis of Alternative PolyAdenylation from RNA-seq) to identify APA with default parameters (hg19 bed

file). PDUI is a unit measure of distal polyadenylation site usage (dPAS)79, a large PDUI value indicates higher distal site usage. To quantify the relative PAS (poly A site) usage, Dapars

defined the percentage of dPAS usage for each sample as PDUI index79. The greater the PDUI is, the more the dPAS of a transcript is used and vice versa. We selected significant APA with

following criterion: FDR ≤ 0.05, |ΔPDUI| = |PDUIinfected − PDUIcontrol| ≥ 0.2, |log2(PDUIinfected/PDUIcontrol)| ≥ 1.5. GENE ONTOLOGY ANALYSIS We separately uploaded differential expression

genes/AS genes/APA genes/Isoform genes into DAVID. DAVID calculated a _p_ value for gene enrichment with a modified Fisher’s exact test and a Benjamin-Hochberg multiple test correction. We

selected significant GO terms with _p__value ≤ 0.05. STATISTICAL ANALYSIS We used R relative packages, such as pheatmap (pheatmap: Pretty Heatmaps, Raivo Kolde, 2015) and VennDiagram

(VennDiagram: Generate High-Resolution Venn and Euler Plots, Hanbo Chen, 2015) and functions, such as cor.test() to analyze data and draw figures. ADDITIONAL INFORMATION HOW TO CITE THIS

ARTICLE: Hu, B. _et al_. Cellular responses to HSV-1 infection are linked to specific types of alterations in the host transcriptome. _Sci. Rep._ 6, 28075; doi: 10.1038/srep28075 (2016).

REFERENCES * Whitley, R. J. & Roizman, B. Herpes simplex virus infections. The Lancet 357, 1513–1518 (2001). Article  CAS  Google Scholar  * Roizman, B. & Whitley, R. J. An inquiry

into the molecular basis of HSV latency and reactivation. Annu Rev Microbiol 67, 355–74 (2013). Article  CAS  PubMed  Google Scholar  * Koelle, D. M. & Corey, L. Herpes simplex: insights

on pathogenesis and possible vaccines. Annu Rev Med 59, 381–95 (2008). Article  CAS  PubMed  Google Scholar  * Padgett, D. A. et al. Social stress and the reactivation of latent herpes

simplex virus type 1. Proc Natl Acad Sci USA 95, 5 (1998). Article  Google Scholar  * Lafaille, F. G. et al. Impaired intrinsic immunity to HSV-1 in human iPSC-derived TLR3-deficient CNS

cells. Nature 491, 769–73 (2012). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Griffiths, S. J. et al. A systematic analysis of host factors reveals a Med23-interferon-lambda

regulatory axis against herpes simplex virus type 1 replication. PLos Pathog 9, e1003514 (2013). Article  CAS  PubMed  PubMed Central  Google Scholar  * Smith, S., Reuven, N., Mohni, K. N.,

Schumacher, A. J. & Weller, S. K. Structure of the Herpes Simplex Virus 1 Genome: Manipulation of Nicks and Gaps Can Abrogate Infectivity and Alter the Cellular DNA Damage Response. J

Virol 88, 10146–10156 (2014). Article  PubMed  PubMed Central  CAS  Google Scholar  * Lilley, C. E., Chaurushiya, M. S., Boutell, C., Everett, R. D. & Weitzman, M. D. The intrinsic

antiviral defense to incoming HSV-1 genomes includes specific DNA repair proteins and is counteracted by the viral protein ICP0. PLos Pathog 7, e1002084 (2011). Article  CAS  PubMed  PubMed

Central  Google Scholar  * Volcy, K. & Fraser, N. W. DNA damage promotes herpes simplex virus-1 protein expression in a neuroblastoma cell line. J Neurovirol 19, 57–64 (2013). Article 

CAS  PubMed  PubMed Central  Google Scholar  * Prasad, A., Remick, J. & Zeichner, S. L. Activation of Human Herpesvirus Replication by Apoptosis. J Virol 87, 10641–10650 (2013). Article

  CAS  PubMed  PubMed Central  Google Scholar  * Wang, X., Patenode, C. & Roizman, B. US3 protein kinase of HSV-1 cycles between the cytoplasm and nucleus and interacts with programmed

cell death protein 4 (PDCD4) to block apoptosis. Proc Natl Acad Sci USA 108, 14632–6 (2011). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Negorev, D. G., Vladimirova, O. V.,

Ivanov, A., Rauscher, F. 3rd & Maul, G. G. Differential role of Sp100 isoforms in interferon-mediated repression of herpes simplex virus type 1 immediate-early protein expression. J

Virol 80, 8019–29 (2006). Article  CAS  PubMed  PubMed Central  Google Scholar  * Wang, Y. et al. p32 Is a Novel Target for Viral Protein ICP34.5 of Herpes Simplex Virus Type 1 and

Facilitates Viral Nuclear Egress. J Biol Chem 289, 35795–35805 (2014). Article  CAS  PubMed  PubMed Central  Google Scholar  * Pasieka, T. J. et al. Functional genomic analysis of herpes

simplex virus type 1 counteraction of the host innate response. J Virol 80, 7600–12 (2006). Article  CAS  PubMed  PubMed Central  Google Scholar  * Rasmussen, S. B. et al. Activation of

Autophagy by alpha-Herpesviruses in Myeloid Cells Is Mediated by Cytoplasmic Viral DNA through a Mechanism Dependent on Stimulator of IFN Genes. J Immunol 187, 5268–5276 (2011). Article  CAS

  PubMed  Google Scholar  * Gillis, P. A., Okagaki, L. H. & Rice, S. A. Herpes simplex virus type 1 ICP27 induces p38 mitogen-activated protein kinase signaling and apoptosis in HeLa

cells. J Virol 83, 1767–77 (2009). Article  CAS  PubMed  Google Scholar  * Jackson, B. R., Noerenberg, M. & Whitehouse, A. A novel mechanism inducing genome instability in Kaposi’s

sarcoma-associated herpesvirus infected cells. PLos Pathog 10, e1004098 (2014). Article  PubMed  PubMed Central  CAS  Google Scholar  * Aubert, M. & Blaho, J. A. Modulation of apoptosis

during herpes simplex virus infection in human cells. Microbes Infect 3, 859–66 (2001). Article  CAS  PubMed  Google Scholar  * Aubert, M. et al. The antiapoptotic herpes simplex virus

glycoprotein J localizes to multiple cellular organelles and induces reactive oxygen species formation. J Virol 82, 617–29 (2008). Article  CAS  PubMed  Google Scholar  * Ferenczy, M. W.,

Ranayhossaini, D. J. & DeLuca, N. A. Activities of ICP0 involved in the reversal of silencing of quiescent herpes simplex virus 1. J Virol 85, 4993–5002 (2011). Article  CAS  PubMed 

PubMed Central  Google Scholar  * Roizman, B. The checkpoints of viral gene expression in productive and latent infection: the role of the HDAC/CoREST/LSD1/REST repressor complex. J Virol

85, 7474–7482 (2011). Article  CAS  PubMed  PubMed Central  Google Scholar  * Everett, R. D. The spatial organization of DNA virus genomes in the nucleus. PLos Pathog 9, e1003386 (2013).

Article  CAS  PubMed  PubMed Central  Google Scholar  * Nojima, T. et al. Herpesvirus protein ICP27 switches PML isoform by altering mRNA splicing. Nucleic Acids Res 37, 6515–27 (2009).

Article  CAS  PubMed  PubMed Central  Google Scholar  * Sedlackova, L. et al. Herpes simplex virus type 1 ICP27 regulates expression of a variant, secreted form of glycoprotein C by an

intron retention mechanism. J Virol 82, 7443–55 (2008). Article  CAS  PubMed  PubMed Central  Google Scholar  * Barzilai, A., Zivony-Elbom, I., Sarid, R., Noah, E. & Frenkel, N. The

herpes simplex virus type 1 vhs-UL41 gene secures viral replication by temporarily evading apoptotic cellular response to infection: Vhs-UL41 activity might require interactions with

elements of cellular mRNA degradation machinery. J Virol 80, 505–513 (2006). Article  CAS  PubMed  PubMed Central  Google Scholar  * Esclatine, A., Taddeo, B., Evans, L. & Roizman, B.

The herpes simplex virus 1 UL41 gene-dependent destabilization of cellular RNAs is selective and may be sequence-specific. Proc Natl Acad Sci USA 101, 3603–3608 (2004). Article  ADS  CAS 

PubMed  PubMed Central  Google Scholar  * Taddeo, B., Zhang, W. & Roizman, B. The UL41 protein of herpes simplex virus 1 degrades RNA by endonucleolytic cleavage in absence of other

cellular or viral proteins. Proc Natl Acad Sci USA 103, 2827–2832 (2006). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Taddeo, B., Zhang, W. & Roizman, B. The herpes

simplex virus host shutoff RNase degrades cellular and viral mRNAs made before infection but not viral mRNA made after infection. J Virol 87, 4516–4522 (2013). Article  CAS  PubMed  PubMed

Central  Google Scholar  * Esclatine, A., Taddeo, B. & Roizman, B. The UL41 protein of herpes simplex virus mediates selective stabilization or degradation of cellular mRNAs. Proc Natl

Acad Sci USA 101, 18165–70 (2004). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Kent, J. R. & Fraser, N. W. The cellular response to herpes simplex virus type 1 (HSV-1)

during latency and reactivation. J Neurovirol 11, 376–383 (2005). Article  CAS  PubMed  Google Scholar  * Pasieka, T. J. et al. Host responses to wild-type and attenuated herpes simplex

virus infection in the absence of Stat1. J Virol 83, 2075–87 (2009). Article  CAS  PubMed  Google Scholar  * Wang, E. T. et al. Alternative isoform regulation in human tissue transcriptomes.

Nature 456, 470–6 (2008). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Verma, D. & Swaminathan, S. Epstein-Barr virus SM protein functions as an alternative splicing

factor. J Virol 82, 7180–8 (2008). Article  CAS  PubMed  PubMed Central  Google Scholar  * Rodrigues, R., Grosso, A. R. & Moita, L. Genome-wide analysis of alternative splicing during

dendritic cell response to a bacterial challenge. PLoS One 8, e61975 (2013). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Valletti, A. et al. Genome-wide analysis of

differentially expressed genes and splicing isoforms in clear cell renal cell carcinoma. PLos One 8, e78452 (2013). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Shukla, S. et

al. CTCF-promoted RNA polymerase II pausing links DNA methylation to splicing. Nature 479, 74–9 (2011). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Tian, B., Hu, J., Zhang,

H. & Lutz, C. S. A large-scale analysis of mRNA polyadenylation of human and mouse genes. Nucleic Acids Res 33, 201–212 (2005). Article  CAS  PubMed  PubMed Central  Google Scholar  *

Ji, Z. & Tian, B. Reprogramming of 3′untranslated regions of mRNAs by alternative polyadenylation in generation of pluripotent stem cells from different cell types. PLos One 4, e8419

(2009). Article  ADS  PubMed  PubMed Central  CAS  Google Scholar  * Takagaki, Y., Seipelt, R. L., Peterson, M. L. & Manley, J. L. The polyadenylation factor CstF-64 regulates

alternative processing of IgM heavy chain pre-mRNA during B cell differentiation. Cell 87, 941–952 (1996). Article  CAS  PubMed  Google Scholar  * Fu, Y. et al. Differential genome-wide

profiling of tandem 3′UTRs among human breast cancer and normal cells by high-throughput sequencing. Genome Res 21, 741–747 (2011). Article  CAS  PubMed  PubMed Central  Google Scholar  *

An, J. J. et al. Distinct role of long 3′UTR BDNF mRNA in spine morphology and synaptic plasticity in hippocampal neurons. Cell 134, 175–187 (2008). Article  CAS  PubMed  PubMed Central 

Google Scholar  * Sandberg, R., Neilson, J. R., Sarma, A., Sharp, P. A. & Burge, C. B. Proliferating cells express mRNAs with shortened 3′untranslated regions and fewer microRNA target

sites. Science 320, 1643–1647 (2008). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Chen, C.-Y. A. & Shyu, A.-B. AU-rich elements: characterization and importance in mRNA

degradation. Trends Biochem Sci 20, 465–470 (1995). Article  CAS  PubMed  Google Scholar  * Di Giammartino, D. C., Nishida, K. & Manley, J. L. Mechanisms and consequences of alternative

polyadenylation. Mol Cell 43, 853–66 (2011). Article  CAS  PubMed  PubMed Central  Google Scholar  * Langmead, B., Trapnell, C., Pop, M. & Salzberg, S. L. Ultrafast and memory-efficient

alignment of short DNA sequences to the human genome. Genome Biol 10, R25 (2009). Article  PubMed  PubMed Central  CAS  Google Scholar  * Trapnell, C., Pachter, L. & Salzberg, S. L.

TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25, 1105–11 (2009). Article  CAS  PubMed  PubMed Central  Google Scholar  * Roberts, A., Pimentel, H., Trapnell, C. &

Pachter, L. Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics 27, 2325–9 (2011). Article  CAS  PubMed  Google Scholar  * Roberts, A., Trapnell, C.,

Donaghey, J., Rinn, J. L. & Pachter, L. Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol 12, R22 (2011). Article  CAS  PubMed  PubMed Central  Google

Scholar  * Trapnell, C. et al. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol 31, 46–53 (2013). Article  CAS  PubMed  Google Scholar  *

Trapnell, C. et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 28, 511–5 (2010).

Article  CAS  PubMed  PubMed Central  Google Scholar  * Huang da, W., Sherman, B. T. & Lempicki, R. A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics

resources. Nat Protoc 4, 44–57 (2009). Article  PubMed  CAS  Google Scholar  * Lilley, C. E. et al. A viral E3 ligase targets RNF8 and RNF168 to control histone ubiquitination and DNA damage

responses. EMBO J 29, 943–955 (2010). Article  CAS  PubMed  PubMed Central  Google Scholar  * Hnilicova, J. et al. The C-terminal domain of Brd2 is important for chromatin interaction and

regulation of transcription and alternative splicing. Mol Biol Cell 24, 3557–68 (2013). Article  CAS  PubMed  PubMed Central  Google Scholar  * LeRoy, G., Rickards, B. & Flint, S. J. The

double bromodomain proteins Brd2 and Brd3 couple histone acetylation to transcription. Mol Cell 30, 51–60 (2008). Article  CAS  PubMed  PubMed Central  Google Scholar  * Doerks, T., Copley,

R. R., Schultz, J., Ponting, C. P. & Bork, P. Systematic identification of novel protein domain families associated with nuclear functions. Genome Res 12, 47–56 (2002). Article  CAS 

PubMed  PubMed Central  Google Scholar  * Li, J. et al. Cbx4 governs HIF-1alpha to potentiate angiogenesis of hepatocellular carcinoma by its SUMO E3 ligase activity. Cancer Cell 25, 118–31

(2014). Article  CAS  PubMed  Google Scholar  * Luis, N. M. et al. Regulation of human epidermal stem cell proliferation and senescence requires polycomb- dependent and -independent

functions of Cbx4. Cell Stem Cell 9, 233–46 (2011). Article  CAS  PubMed  Google Scholar  * U. M., Miyashita, T., Shikama, Y., Tadokoro, K. & Yamada, M. Molecular cloning and

characterization of six novel isoforms of human Bim, a member of the proapoptotic Bcl-2 family. FEBS Lett 509, 7 (2001). Google Scholar  * O’Connor, L. et al. Bim a novel member of the Bcl-2

family that promotes apoptosis. EMBO J 17, 12 (1998). Google Scholar  * Weber, K., Harper, N., Schwabe, J. & Cohen, G. M. BIM-mediated membrane insertion of the BAK pore domain is an

essential requirement for apoptosis. Cell Rep 5, 409–20 (2013). Article  CAS  PubMed  PubMed Central  Google Scholar  * Jones, E. V., Dickman, M. J. & Whitmarsh, A. J. Regulation of

p73-mediated apoptosis by c-Jun N-terminal kinase. Biochem J 405, 617–623 (2007). Article  CAS  PubMed  PubMed Central  Google Scholar  * Li, Y. Y. et al. Pim-3, a proto-oncogene with

serine/threonine kinase activity, is aberrantly expressed in human pancreatic cancer and phosphorylates bad to block bad-mediated apoptosis in human pancreatic cancer cell lines. Cancer Res

66, 6741–7 (2006). Article  CAS  PubMed  Google Scholar  * Popivanova, B. K. et al. Proto-oncogene, Pim-3 with serine/threonine kinase activity, is aberrantly expressed in human colon cancer

cells and can prevent Bad-mediated apoptosis. Cancer Sci 98, 321–8 (2007). Article  CAS  PubMed  Google Scholar  * Wilkinson, D. E. & Weller, S. K. Recruitment of cellular recombination

and repair proteins to sites of herpes simplex virus type 1 DNA replication is dependent on the composition of viral proteins within prereplicative sites and correlates with the induction

of the DNA damage response. J Virol 78, 4783–4796 (2004). Article  CAS  PubMed  PubMed Central  Google Scholar  * Sharma, A., Singh, K. & Almasan, A. Histone H2AX phosphorylation: a

marker for DNA damage. Methods Mol Biol 920, 613–26 (2012). Article  CAS  PubMed  Google Scholar  * Tanaka, T., Halicka, D., Traganos, F. & Darzynkiewicz, Z. Cytometric analysis of DNA

damage: phosphorylation of histone H2AX as a marker of DNA double-strand breaks (DSBs). Methods Mol Biol 523, 161–8 (2009). Article  CAS  PubMed  Google Scholar  * Yonemasu, R. et al.

Disruption of mouse XAB2 gene involved in pre-mRNA splicing, transcription and transcription-coupled DNA repair results in preimplantation lethality. DNA Repair (Amst) 4, 479–91 (2005).

Article  CAS  Google Scholar  * Kuraoka, I. et al. Isolation of XAB2 complex involved in pre-mRNA splicing, transcription and transcription-coupled repair. J Biol Chem 283, 940–50 (2008).

Article  CAS  PubMed  Google Scholar  * Wang, Z. H., Castano, I. B., De Las Penas, A., Adams, C. & Christman, M. F. Pol kappa: A DNA polymerase required for sister chromatid cohesion.

Science 289, 774–779 (2000). Article  ADS  CAS  PubMed  Google Scholar  * Takahasi, K. et al. Nonself RNA-sensing mechanism of RIG-I helicase and activation of antiviral immune responses.

Mol Cell 29, 428–40 (2008). Article  CAS  PubMed  Google Scholar  * Yoneyama, M. et al. The RNA helicase RIG-I has an essential function in double-stranded RNA-induced innate antiviral

responses. Nat Immunol 5, 730–7 (2004). Article  CAS  PubMed  Google Scholar  * Li, C. Z. et al. Polymorphism of OAS-1 determines liver fibrosis progression in hepatitis C by reduced ability

to inhibit viral replication. Liver Int 29, 1413–21 (2009). Article  CAS  PubMed  PubMed Central  Google Scholar  * Knapp, S. et al. Polymorphisms in interferon-induced genes and the

outcome of hepatitis C virus infection: roles of MxA, OAS-1 and PKR. Genes Immun 4, 411–9 (2003). Article  CAS  PubMed  Google Scholar  * Medzhitov, R., Preston-Hurlburt, P. & Janeway

Jr., C. A. A human homologue of the Drosophila Toll protein signals activation of adaptive immunity. Nature 388, 4 (1997). Article  CAS  Google Scholar  * Chowdhury, P., Sacks, S. H. &

Sheerin, N. S. Toll-like receptors TLR2 and TLR4 initiate the innate immune response of the renal tubular epithelium to bacterial products. Clin Exp Immunol 145, 346–56 (2006). Article  CAS

  PubMed  PubMed Central  Google Scholar  * Zhou, X. et al. Transcriptome analysis of alternative splicing events regulated by SRSF10 reveals position-dependent splicing modulation. Nucleic

Acids Res 42, 4019–30 (2014). Article  CAS  PubMed  PubMed Central  Google Scholar  * Blattner, C. et al. UV-Induced Stabilization of c-fos and Other Short-Lived mRNAs. Mol Cell Biol 20, 10

(2000). Article  Google Scholar  * Jurado, J., Fuentes-Almagro, C. A., Prieto-Alamo, M. J. & Pueyo, C. Alternative splicing of c-fos pre-mRNA: contribution of the rates of synthesis and

degradation to the copy number of each transcript isoform and detection of a truncated c-Fos immunoreactive species. BMC Mol Biol 8, 83 (2007). Article  PubMed  PubMed Central  CAS  Google

Scholar  * Xia, Z. et al. Dynamic analyses of alternative polyadenylation from RNA-seq reveal a 3′-UTR landscape across seven tumour types. Nat Commun 5, 5274 (2014). Article  ADS  CAS 

PubMed  Google Scholar  * Baltz, A. G. et al. The mRNA-bound proteome and its global occupancy profile on protein-coding transcripts. Molecular Cell 46, 674–690 (2012). Article  CAS  PubMed

  Google Scholar  * Hendrickx, A. et al. Docking motif-guided mapping of the interactome of protein phosphatase-1. Chem Biol 16, 365–371 (2009). Article  CAS  PubMed  Google Scholar  *

Eckert, D., Buhl, S., Weber, S., Jager, R. & Schorle, H. The AP-2 family of transcription factors. Genome Biol 6, 246 (2005). Article  PubMed  PubMed Central  CAS  Google Scholar  *

Ogami, K., Cho, R. & Hoshino, S.-i. Molecular cloning and characterization of a novel isoform of the non-canonical poly (A) polymerase PAPD7. Biochem Biophys Res Commun 432, 135–140

(2013). Article  CAS  PubMed  Google Scholar  * Semba, S. et al. Biological functions of mammalian Nit1, the counterpart of the invertebrate NitFhit Rosetta stone protein, a possible tumor

suppressor. J Biol Chem 281, 28244–28253 (2006). Article  CAS  PubMed  Google Scholar  * Ono, M. et al. Foxp3 controls regulatory T-cell function by interacting with AML1/Runx1. Nature 446,

685–689 (2007). Article  ADS  CAS  PubMed  Google Scholar  * Katz, Y., Wang, E. T., Airoldi, E. M. & Burge, C. B. Analysis and design of RNA sequencing experiments for identifying

isoform regulation. Nat Methods 7, 1009–U101 (2010). Article  CAS  PubMed  PubMed Central  Google Scholar  * Wang, E. T. et al. Transcriptome-wide regulation of pre-mRNA splicing and mRNA

localization by muscleblind proteins. Cell 150, 710–24 (2012). Article  CAS  PubMed  PubMed Central  Google Scholar  * Rutkowski, A. J. et al. Widespread disruption of host transcription

termination in HSV-1 infection. Nat Commun 6, 7126 (2015). Article  ADS  PubMed  Google Scholar  * Lenzken, S. C., Loffreda, A. & Barabino, S. M. RNA splicing: a new player in the DNA

damage response. Int J Cell Biol 2013, 153634 (2013). Article  PubMed  PubMed Central  CAS  Google Scholar  * Hosoda, N., Lejeune, F. & Maquat, L. E. Evidence that poly(A) binding

protein C1 binds nuclear pre-mRNA poly(A) tails. Mol Cell Biol 26, 3085–97 (2006). Article  CAS  PubMed  PubMed Central  Google Scholar  * Chansky, H. A., Hu, M., Hickstein, D. D. &

Yang, L. Oncogenic TLS/ERG and EWS/Fli-1 fusion proteins inhibit RNA splicing mediated by YB-1 protein. Cancer research 61, 3586–3590 (2001). CAS  PubMed  Google Scholar  * Stickeler, E. et

al. The RNA binding protein YB‐1 binds A/C‐rich exon enhancers and stimulates splicing of the CD44 alternative exon v4. EMBO J 20, 3821–3830 (2001). Article  CAS  PubMed  PubMed Central 

Google Scholar  * Stoecklin, G. et al. MK2-induced tristetraprolin14-3-3 complexes prevent stress granule association and ARE-mRNA decay. EMBO J 23, 12 (2004). Article  CAS  Google Scholar 

* Brooks, S. A., Connolly, J. E. & Rigby, W. F. C. The Role of mRNA Turnover in the Regulation of Tristetraprolin Expression: Evidence for an Extracellular Signal-Regulated

Kinase-Specific, AU-Rich Element-Dependent, Autoregulatory Pathway. The Journal of Immunology 172, 7263–7271 (2004). Article  CAS  PubMed  Google Scholar  Download references

ACKNOWLEDGEMENTS We thank the excellent technical assistance of Ms. Liping Yang. This work was in part supported by grants from Chinese Academy of Sciences (KSCXZ-EW-BR-6), Startup fund from

Kunming Institute of Zoology, Chinese Academy of Sciences (Y102421081), Grants from Yunnan Provincial Government (2013FA051), National Science Foundation of China (NSFC 81471966), a Common

Project of Panzhihua, Science and Technology Bureau of China (2012CY-S-22(9)) and a Visiting professorship for senior international scientist from CAS to NWF (2012T1S0001). AUTHOR

INFORMATION Author notes * Hu Benxia and Li Xin contributed equally to this work. AUTHORS AND AFFILIATIONS * Key Laboratory of Animal Models and Human Disease Mechanisms of the Chinese

Academy of Sciences & Yunnan Province, Kunming Institute of Zoology, Kunming, 650223, Yunnan, China Benxia Hu, Xin Li, Yongxia Huo, Yafen Yu, Qiuping Zhang, Guijun Chen & Jumin Zhou

* Kunming College of Life Science, University of Chinese Academy of Sciences, Kunming, 650204, Yunnan, China Benxia Hu & Xin Li * School of Life Sciences, Anhui University, Hefei,

230601, Anhui, China Benxia Hu & Yafen Yu * State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, Yunnan,

China Yaping Zhang & Dongdong Wu * Department of Microbiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, 19104, PA, USA Nigel W. Fraser Authors * Benxia Hu

View author publications You can also search for this author inPubMed Google Scholar * Xin Li View author publications You can also search for this author inPubMed Google Scholar * Yongxia

Huo View author publications You can also search for this author inPubMed Google Scholar * Yafen Yu View author publications You can also search for this author inPubMed Google Scholar *

Qiuping Zhang View author publications You can also search for this author inPubMed Google Scholar * Guijun Chen View author publications You can also search for this author inPubMed Google

Scholar * Yaping Zhang View author publications You can also search for this author inPubMed Google Scholar * Nigel W. Fraser View author publications You can also search for this author

inPubMed Google Scholar * Dongdong Wu View author publications You can also search for this author inPubMed Google Scholar * Jumin Zhou View author publications You can also search for this

author inPubMed Google Scholar CONTRIBUTIONS B.H. Data analysis, RT-PCR and qRT-PCR validation and manuscript writing. X.L. RNA-sequencing, qRT-PCR validation and manuscript writing. Y.H.

Data analysis and analysis method. Y.Y. RT-qPCR validation. Q.Z. Writed the AS and APA part of the paper. G.C. Cell culture related experiments. Y.Z. Analyses methods. D.W. Data analyses and

manuscript. N.W.F. Writing manuscript. J.Z. Experiment design, data analyses and manuscript writing. All authors read and approved the final manuscript. ETHICS DECLARATIONS COMPETING

INTERESTS The authors declare no competing financial interests. ELECTRONIC SUPPLEMENTARY MATERIAL SUPPLEMENTARY INFORMATION RIGHTS AND PERMISSIONS This work is licensed under a Creative

Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in

the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of

this license, visit http://creativecommons.org/licenses/by/4.0/ Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Hu, B., Li, X., Huo, Y. _et al._ Cellular responses to HSV-1

infection are linked to specific types of alterations in the host transcriptome. _Sci Rep_ 6, 28075 (2016). https://doi.org/10.1038/srep28075 Download citation * Received: 01 March 2016 *

Accepted: 26 May 2016 * Published: 29 June 2016 * DOI: https://doi.org/10.1038/srep28075 SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get

shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative