Longitudinal transcriptomic characterization of viral genes in HSV-1 infected tree shrew trigeminal ganglia
Virology Journal volume 17, Article number: 95 (2020)
Following acute infection, Herpes Simplex virus-1 (HSV-1) establishes lifelong latency and recurrent reactivation in the sensory neurons of trigeminal ganglia (TG). Infected tree shrew differs from mouse and show characteristics similar to human infection. A detailed transcriptomic analysis of the tree shrew model could provide mechanistic insights into HSV-1 infection in humans.
We sequenced the transcriptome of infected TGs from tree shrews and mice, and 4 human donors, then examined viral genes expression up to 58 days in infected TGs from mouse and tree shrew, and compare the latency data with that in human TGs.
Here, we found that all HSV-1 genes could be detected in mouse TGs during acute infection, but 22 viral genes necessary for viral transcription, replication and viral maturation were not expressed in tree shrew TGs during this stage. Importantly, during latency, we found that LAT could be detected both in mouse and tree shrew, but the latter also has an ICP0 transcript signal absent in mouse but present in human samples. Importantly, we observed that infected human and tree shrew TGs have a more similar LAT region transcription peak. More importantly, we observed that HSV-1 spontaneously reactivates from latently infected tree shrews with relatively high efficiency.
These results represent the first longitudinal transcriptomic characterization of HSV-1 infection in during acute, latency and recurrent phases, and revealed that tree shrew infection has important similar features with human infection.
HSV-1 is a ubiquitous but important human pathogen carried by over half of the world’s population; HSV-1 infection starts with primary infection at the periphery and subsequent lifelong latency in the peripheral nervous system . In experimental animals such as mouse, acute infection develops following cornea inoculation, the virus replicates in the epithelial cells on the corneal surface and is later transported into trigeminal ganglia, where it establishes latency [2,3,4]. The acute stage of infection involves lytic, or productive infection of HSV-1 at the site of inoculation, and frequently in the TGs of infected mice; in the tree shrew TGs, however, signs of lytic infected were not seen and no infectious virus could be detected . During the lytic phase of infection in cultured cells, all viral genes are believed to be expressed in a cascade-dependent manner [6, 7], but during latency, most viral genes are silenced with the exception of the latency-associated transcript (LAT), multiple miRNAs [8,9,10] and two small RNAs . Although latent, the virus reactivates from individual neurons periodically, and could cause more serious diseases including herpes keratitis or herpes encephalitis [12, 13]. Recently, HSV-1 has been suggested to play a role in Alzheimer’s disease [14, 15].
Multiple animal species including mice, rabbits and tree shrews [5, 16, 17] have been used to model human HSV-1 infection, with mice being the most widely used. However, we have recently studied HSV-1 infection in tree shrews, which are more closely related to primates than rodents , and observed differences between tree shrews and mice during acute infection, latency and reactivation [5, 19]. Here, we performed an in-depth transcriptional profiling of the infected mouse, tree shrew and human TGs, and reveal the differences in viral gene expression patterns. Importantly, we show that latent tree shrew TGs can express many viral genes, including UL6, UL8 and ICP8, which are consistent with spontaneous reactivation. More importantly, during latency tree shrew and human TGs had more similar transcription peaks in the LAT region and possibly more ICP0 transcription, supporting that the tree shrew model better mimics HSV-1 latent infection in human.
Viral infection dynamics vary among animal models
To analyze the viral transcriptional patterns during acute and latent stages of HSV-1 infection, mice and tree shrews were infected with HSV-1 strain 17+ by corneal scarification, and infected TGs collected over a period of 58 days (Fig. 1a, b). For each time point, three biological replicates were generated. Human TG samples were also collected for comparison with latently infected mouse and tree shrew TGs (Fig. 1a).
In total, we sequenced 52 samples (24 mouse TGs, 24 tree shrew TGs and 4 human TGs) using an Illumina HiSeq Platform and generated about 5.9 million reads per sample. After mapping to the genome of the hosts and HSV-1 strain 17+, we performed differentially expressed gene (DEG) detection using DESeq2 (differential expression analysis for sequence count data2) , in which we compared viral and host gene expression between samples with different infection times and uninfected control samples or each pair of consecutive time points to determine what viral or host genes are expressed significantly. The host genes were then filtered out and are not discussed further here. When using read count > 10 as the detection threshold of viral genes, we found that the number of viral genes were quite different between mice and tree shrews (Fig. 1c). Almost all of the 74 unique viral genes (LAT, ICP0, ICP34.5 and ICP4 all have a duplicated in HSV-1 genome) were expressed in mouse samples from 3 to 7 days post infection (dpi), while only a total of 52 viral genes were detected during the acute phase (3 to 10 dpi) in tree shrew TGs. The differences extended to latency, only LAT was detected in mouse TGs. However, infected tree shrew TGs also showed a high level of ICP0 transcript signal at both 28 and 58 dpi (Fig. 2a).
To provide a dynamic change in viral genes throughout the infection period, we plotted the distribution of absolute fold changes in viral genes expression for each model separately. For both models, the largest differences were observed between the 0 and 3dpi. After 3dpi, the tree shrew model showed an abrupt decline until 10dpi, then increased between 10 and 14dpi, and then slowly declined again until 58dpi. The differences were that mouse model displayed an ascending change during 5 to 10 dpi and showed no significant change in the range of 14 to 28dpi (p > 0.05, ANOVA test, Fig. 1d). Next, we asked whether there were differences in the up- or down-regulation of viral genes during this period between the two models. The number of differentially expressed genes for each pairwise comparison and each model separately was shown in Fig. 1e, which was highest in mouse model between 0 and 3dpi (up-regulated), and showed down-regulated genes within 5 to 10dpi. This demonstrates that the infected moue TG showed the largest changes in viral genes expression levels during acute phase, while in tree shrew the changes are mild in acute stage, but more significant changes occur during latency. Therefore, it appears that strong acute infections tend to lead to a calm latency, while the milder acute stage results in a latency with greater change in viral transcription and a higher frequency of reactivation.
The acute infection in mouse TG is more consistent with productive infection based on the expression of essentially all HSV-1 genes, therefore we compared mouse TG viral gene expression with that in human primary fibroblasts at 2, 4 and 8 h post infection (hpi), and found that the transcriptional pattern of lytic infection was similar to that in the mouse model (Supplemental Figure 1). When considering that mice have a higher mortality than tree shrews during the acute stage  and immunohistochemistry failed to detect any HSV-1 antigens in acutely infected tree shrew TGs while positive signals were detected in mouse TGs , these results support the conclusion that HSV-1 can build a highly active productive infection in mouse, but only mount a limited or abortive infection in tree shrew TGs.
Viral gene expression patterns differ between infected mouse and tree shrew TGs
To determine the temporal pattern of viral gene expression following inoculation, we generated a heat map of the viral DEGs expression data. As shown in Fig. 2a, in mouse TGs essentially all HSV-1 genes were detected with higher transcription levels during the acute stage, but in tree shrew TGs, only some of these genes were detected. During latency, from day 28–58, tree shrew and mouse differed in that tree shrew TGs expressed ICP0 transcript, in addition to LAT, while mouse TGs express only LAT with little or no ICP0 signal.
Next, we focused on the viral genes that are not expressed or are expressed at extremely low levels in acutely infected tree shrew TGs, but are expressed in acutely infected mouse TGs. We found that there are at least 22 such genes (Fig. 2b). Using UniProt database and HSV-1 17+ genome annotation information in NCBI (Supplemental Table 1), we sorted these 22 viral genes into 6 clusters (Fig. 2c). Four of these clusters contain at least 3 viral genes, and these are tegument proteins, envelope glycoproteins, DNA packaging proteins and nuclear proteins, while the other clusters included ICP4 and UL9. Selected RT-qPCR validation of these transcripts is shown in Fig. 2d.
Importantly, of the 22 undetected genes, 6 are envelope proteins, functioning in viral entry and membrane fusion. Thus HSV-1 infected cells from tree shrew TGs are likely missing these essential gene products for viral maturation and subsequent infection. With most of the important late structural genes not being made, HSV-1 in acutely infected tree shrew TGs was likely undergoing an abortive lytic infection, which explains the low toxicity and absence of infectious HSV-1 in tree shrews. Furthermore, two of the viral genes not detectable in infected tree shrew TGs are ICP4 and UL9. ICP4 mainly acts as a transactivator of viral transcription, which is essential for expression of almost all early and late genes; UL9 encodes the viral origin binding protein, which acts as an ATPase and helicase and is required for initiation of DNA synthesis . The lack of detectable expression of these genes could account for the fact that viral transcription is limited in tree shrew TGs.
Many HSV-1 genes that are below the detection threshold by immunohistochemistry or immunofluorescence were expressed at extremely low levels in tree shrew TGs
In our previous study, we could detect HSV-1 antigens by immunohistochemistry and in particular ICP4 by immunofluorescence in acutely infected mouse TGs, but we could not detect any HSV-1 proteins in acutely infected tree shrew TGs by immunohistochemistry using the same anti-HSV-1 polyclonal antibodies, or by immunofluorescence using monoclonal antibody against viral essential genes ICP0 and ICP4 . Consistent with this data, ICP4 transcript was noticeably nearly absent by Real Time PCR analysis (Fig. 2d). When examining the sequencing data, we noted that ICP4 reads were extremely low. We selected ICP4 and unique short region (US) genes for further analyses, as this region also includes genes US2, US3 and US5, which are also below the detection threshold (log2 fold change > 1, and adjusted P value < 0.01). The sequence coverage profiles are shown in Fig. 3a. The maximal mean ICP4 reads in mouse TGs were very high, at nearly 12,000, contrasting with tree shrew which were only 200. A more common phenomenon in the data from acutely infected tree shrew TGs is that many viral genes show incomplete coverage in coding regions. To determine whether this incomplete coverage reflects the difference in length of the transcript, we designed primers to amplify about 1 kb fragments of selected viral genes (Fig. 3b and c). We also sequenced these PCR products and found that they correctly aligned to the viral genome. These results support the existence of very low levels of intact transcripts in areas with little or no reads in RNAseq.
The uneven coverage is seen in many HSV-1 genes and is likely due to the low efficiency in amplification during sequencing as a result of the high GC-content of the viral genome . Consistent with this possibility, the detected reads correspond to reads that are also abundant in mouse samples (see red boxes in Fig. 3a), while the region with little or no reads correspond to relatively low reads in the mouse samples (see green boxes in Fig. 3a). We assessed the relationship between GC-content and genome coverage in tree shrew acute stage samples (Fig. 3d and Supplemental Figure 2A). These result support that the “coverage” or “non-coverage” in viral genome is not related with GC-content in the case of a large amount of virus genes transcription (tree shrew 5dpi, 7dpi and 10dpi samples); however, if HSV-1 genome has a low transcriptional activity, GC-content has a significant effect on genome coverage. The non-covered region in tree shrew samples both has low GC-content sites and high GC-content sites (Fig. 3e and Supplemental Figure 2B), which indicated that high GC-content is not the only reason for those undetectable viral genes in tree shrew samples.
Since ICP4 is essential for viral transcription, it’s extremely low level of expression contributes to the lack of a productive infection in tree shrew TGs [22,23,24,25]. The lack of ICP4 is likely due to transcriptional inhibition and not a result of degradation of the ICP4 protein. In our previous study, we had reported that HSV-1 could become latent and later spontaneously reactivate from tree shrew TGs, and that the virus could mount a full productive infection in tree shrew CNS in vivo , and in vitro in cultured ganglion neurons from tree shrew TGs (Supplemental Figure 3), thus, ICP4 protein can function in tree shrew neurons, and the low expression is more likely due to the specific inhibition of the ICP4 transcription in TG neurons in vivo. Nonetheless, we believe, that the extremely low level of ICP4 transcripts might have been sufficient for the virus to sustain a low level of genome replication as reported before  and subsequently achieve latency in the tree shrew without triggering a full blown infection and immune response.
Spontaneous reactivation is detected during latency in tree shrew TGs
RNA sequencing data suggested that low levels of viral transcripts were widely present in infected tree shrew samples as could be seen in samples beyond day 14, contrasting that of mouse TGs, where the average reads are close to zero (Fig. 4a). To understand the distribution of low levels of viral genes in tree shrew TGs, we separated the three replicates of latency samples and analyzed them individually. We found that the low levels of viral transcripts were from individual tree shrew TGs and not a result of uniform low level of expression in all samples. As seen in Fig. 4b, sample #3 from day 28 and #1 from day 58 appear to express many more viral genes than the remaining samples.
We then calculated the significantly transcribed viral genes in these sample using an R package, DEGseq, which allowed us to identify differentially expressed genes in each latent sample individually . Several viral genes, such as ICP8, UL6, UL8 were found to be differentially expressed in more than one latent tree shrew samples (Fig. 4c), and therefore we validated some selected viral genes using RT-qPCR from 58dpi samples (Fig. 4d). The results indicate that tree shrew #3 from day 28 sample and #1 from day 58 sample are likely undergoing spontaneously reactivation as predicted from our earlier study analyzing tree shrew eye swabbing .
Human and tree shrew are similar in HSV-1 ICP0 and LAT expression pattern in latently infected TGs
Studies using the mouse and other models have established that during HSV-1 latent infection, almost all viral genes are silenced with the exception of the LAT region, which is expressed and produces a stable intron [28, 29], small and microRNAs [8, 30], and helps to maintain latency or reactivation. This is indeed the case in mouse with HSV-1 transcriptome profiling over the course of 2 months post infection. However, viral gene expression in infected tree shrew TGs, with the additional expression of ICP0 transcript during latency, differs significantly from that in mouse TGs. This is evident in Fig. 2a, where both the LAT and the ICP0 region were transcribed at robust levels in tree shrews during latency, whereas only LAT signal was found in mouse.
To determine what HSV-1 genes are expressed in human TGs and which model is more similar to human in terms of ICP0 and LAT expression patterns during latency, we analyzed human TG sequencing data, and found that the LAT region and ICP0 region were also activated in human TGs, and no other viral genes were detected. To confirm that the observed ICP0 transcription was not an artifact of the specific analytic tool, and also to account for the fact that the region of the genome encoding ICP0 overlaps the LAT region, we used individual strand specific primers for ICP0 and LAT and validated the expression of ICP0 (Fig. 5a), and the strand specificity of this approach had been demonstrated in supplemental Figure 4. When this assay was applied to amplify a longer ICP0 transcript (657pb), we only detected the ICP0 signal in the latently infected tree shrew samples (supplemental Figure 4C). In our previous study, ICP0 transcript was also detected by in situ hybridization in 58dpi tree shrew TGs . Human TGs were harvested a considerable time (1.5–2.5 days) after death, which may result in HSV-1 reactivation, but no viral genes other than LAT and ICP0 were detected in the RNAseq data, and whether this data can really reflect the latent state in the living human is a question we cannot answer.
Although most viral genes were at lower levels in tree shrew than that in mouse during acute stage, the level of LAT was higher in tree shrew TGs than in mouse TGs. This trend continued into latency where LAT is expressed at a much higher level in tree shrew TGs than in mouse (Fig. 5b, c). To reveal the difference of viral gene expression during latency, we compared the LAT region transcription patterns between human and animal models. Two of four human samples we collected (Fig. 5d) and 2 of 16 GEO data (Supplemental Figure 5) provided significant LAT signals. Moreover, the LAT transcription patterns were more similar to those of tree shrews, and they both had a transcription peak in the 120,758–121,214 region, which was not observed in the mouse latency samples (Fig. 5c, d). We also compared the coverage of the peak region (120,758-121,214), LAT intron (remove peak region) and LAT exon (Fig. 5e), and tree shrew and human have higher readings in the peak region. Considering that the animal model is artificially infected, and the human sample is a natural infection, which may result in a low LAT reading in human samples.
HSV-1 infection in tree shrew TGs differs from that in mouse, and viral transcriptional level in tree shrews and mice are affected by many factors. In our previous relative HSV-1 genome quantitative data, acutely infected mouse TGs contained higher levels of viral genome than tree shrew TGs , which also contribute to the differences of viral transcription. Another factor that cannot be ignored is whether HSV-1 replicates effectively in the eye of the tree shrew model, because surface replication is the source of the virus transported into the ganglion. In fact, we had tested HSV-1 titer and genome copy number in tree shrew eyes within a period of 46 days after ocular infection (Li et al., in press), and these results support that HSV-1 replication in tree shrew eyes could be comparable with that in mouse eyes.
An important finding from this study is the apparent spontaneous reactivation from latency from two of the tree shrew TGs sequenced, a result supporting our previous study where live virus could be recovered from eye tears in latently infected animals---a feature that mimics human infections seen in tree shrew but less efficient in mouse. From genes detected during spontaneous reactivation, we found immediate early genes ICP0, ICP22 and ICP8, which are important for viral early transcription and genome replication. However, a number of studies support the spontaneous reactivation of infectious HSV-1 in murine sensory ganglia [31, 32], and ICP0 signal was detected in individual latently infected mice in our RNAseq data, but other lytic viral genes reads are too low to tell if they are positive.
What caused the virus to spontaneously reactivate is currently unknown. LAT is the only known transcript in latently infected tissues [10, 33, 34], but it is not just a latency marker, it maintains viral latency [35, 36] by inhibiting apoptosis [37,38,39,40], repressing lytic gene expression , etc. There have been many studies supporting LAT is essential for the highly induced reactivation phenotype in the mouse model [42,43,44], and for the high induced or spontaneous reactivation in the rabbit model [45,46,47,48,49,50,51,52,53]. If true, then the high level of LAT transcript in tree shrew TGs compared to mouse could partially explain why the virus spontaneously reactivates in tree shrew but not in mouse.
In our analyses comparing latently infected mouse, tree shrew and human TGs, we observed that LAT could be detected from all three, but tree shrew TG was more similar to human TG in that both have high reads in a LAT intron region. Considering that the oligo (dT) method cannot enrich non-polyadenylated LATs, the high level of LAT signal are mainly supported by polyadenylated LATs [54,55,56,57], the signal in the LAT intron region is in fact ICP0, but not LAT intron, which is not polyadenylated. This is supported by the strand-specific RT-qPCR experiments shown in Fig. 5a.
A more important similarity between infected tree shrew and human TGs is that both express transcript from the ICP0 region. In contrast, there are no ICP0 transcripts detected in infected mouse TGs, which is inconsistent with some previous reports that the ICP0 signal they detected was located near the first exon at the 5′ end [58, 59]. In light of this, this comparison could also include infected rabbits TGs in future analysis to examine whether this ICP0 transcripts exist in reactivating samples.
Both in tree shrew and human TGs, most of the ICP0 reads are concentrated at the 3′ end of the ICP0 gene (Fig. 5c & d), while the remaining region of ICP0 has very few reads, raising the possibility that this region is independently transcribed. However, it is also possible that the uneven reads were a result of the high GC content of the HSV genome. Since we could not amplify the full-length ICP0 transcript, we could not rule out the possibility of partial transcripts arising from the 3′ of ICP0. In addition, there is a 0.8 kb LAT insulator and a CTCF binding site within LAT intron region , and these elements are both located up-stream of the LAT/ICP0 transcription peak, making it possible that the CTCF binding site could serve as a promoter to for this transcript, as reporter by others .
In this study, we compared the viral transcriptome of infected mouse and tree shrew TGs during the course of infection, and we compared the pattern of viral transcripts in latently infected mouse and tree shrew TGs with that of human TGs. We found that HSV-1 transcription in acutely infected TGs differs dramatically between mouse and tree shrew, with HSV-1 in mouse TG undergoing productive infection, while that in tree shrew TGs appears to be going through an abortive infectious cycle, missing keys genes needed for viral transcription, replication and maturation. During the latent phase of infection, LAT was detected in mouse, tree shrew and human TGs, but we also detected an ICP0 transcript fragment from tree shrew and human TGs, making tree shrew latent infection more similar to human than mouse. When we analyzed tree shrew TGs individually, we found samples that appeared to be undergoing spontaneous reactivation. Together these analyses support the tree shrew as a better model of human HSV-1 infection in the peripheral nervous system, offering the possibility of a better understand HSV-1 latency and reactivation, and the discovery of potential novel targets for therapeutic interventions. Taken together, the transcriptome data reveals that tree shrews and humans have a more similar transcription pattern in the LAT region during latency than that of mice and humans, supporting tree shrews as a more accurate animal model for research on HSV-1 latency and reactivation.
Materials and methods
Chinese tree shrews, 6-month-old female, were obtained from the experimental animal core facility of the Kunming Institute of Zoology, Chinese Academy of Sciences. During the experiment, tree shrews were kept in experimental cages of 54 cm × 45 cm × 50 cm, and no more than two animals per cage. The temperature in the room was controlled at 15 ~ 28 °C, the relative humidity was 40 ~ 70%, the daily light was 12 h, and the noise should not exceed 60 dB. BALB/C mice, 6-week-old female, were obtained from the Kunming Medical University. All experimental procedures and animal care were carried out in accordance with the protocols approved by the Institutional Animal Care and Use Committee of the Kunming Institute of Zoology, Chinese Academy of Sciences. The research program was reviewed and approved by the Institutional Animal Care and Use Committee of Kunming Institute of Zoology, Chinese Academy of Sciences.
Virus and cells
HSV-1 strain 17+ was used to infect mouse and tree shrews. Viral culture was performed in Vero cells and titrated by plaque forming assay on RS1 cells, and these two cells were both obtained from Conservation Genetics Academy of Science (CAS) Kunming Cell Bank. The infected cells were cultured in DMEM supplemented with 2% fetal bovine serum (GibcoTM). All cells were maintained at 37 °C with 5% CO2. Experiments involving infectious virus were conducted in a Biosafety level 2 laboratory.
Animal infections and collection of samples
Animal anesthesia, corneal scarification and inoculation with HSV-1 17+ virus was as previously described . 1 × 104 PFU of HSV-1 were used to inoculate each mouse eye, and 1 × 106 PFU were used on tree shrews, which was consistent with our previous reports [5, 19]. The control (mock infected) animals were also scratched and treated with Vero cell supernatant. To generate biological replicates, three groups of animals were infected independently on different days. Samples of infected trigeminal ganglions were collected at 3, 5, 7, 10, 14, 28, 58dpi, and mock infected samples were collected after 24 h of treatment, and then were ground to a fine powder in liquid nitrogen using 50 ml grinding beakers and 20-mm grinding ball for RNA extraction.
Human trigeminal ganglions collection and preparation
Four human trigeminal ganglia were obtained at autopsy and provided by Academy of Forensic Science (Shanghai, China) and Second Affiliated Hospital of Kunming Medical University. According to the record, the material was taken after 2–3 days of death, then stored in liquid nitrogen and transported with dry ice, and demographics are shown in supplemental Table 3. This is about as early as it is possible to legally obtain human tissue for our studies. Tissues were stored at − 80 °C until processing. To disrupt the tissue, frozen ganglia were wrapped in aluminum foil and mechanically broken apart on dry ice using a metal anvil and hammer, both also frozen on dry ice. Small fragments of tissue from a trigeminal ganglion from were used for RNA extraction.
RNA extraction and sequencing
Mouse, tree shrew or human TGs were individually ground to a fine powder in liquid nitrogen before RNA extraction. For each sample, approximately 1 g of powder (from two TGs of one animal) was resuspended in 1 ml TRIzol reagent (Life Technologies) and total RNA was extracted according to the manufacturer’s recommendation. Afterwards, RNA samples were sent to BGI (Wuhan, China) for purification, library preparation and sequencing. The Ambion Trubo DNA-free kit (Life Technologies) was used to eliminate genomic DNA contamination, and an Agilent 2100 Bioanalyzer (Agilent RNA 6000 Nano Kit) to perform the total RNA sample quality control. mRNAs were isolated from total RNA using the oligo (dT) method, and then purified and fragmented using divalent cations under elevated temperature. The oligo (dT) method could not extract, and thus excluded non-polyadenylated transcripts, such as non-polyadenylated LAT introns, and retains only polyadenylated LATs [54,55,56,57]. First strand cDNA was synthesized using random primers, and the second strand cDNA was synthesised with Polymerase I and RNase H. cDNA fragments were purified and resolved with EB buffer for end reparation and single nucleotide A (adenine) addition. After that, the cDNA fragments were linked with adapters. Those cDNA fragments with suitable size (300 bp) were selected for PCR amplification. An Agilent 2100 Bioanaylzer and ABI StepOnePlus Real-Time PCR System were used in quantification and qualification of the libraries. Finally, all of the transcriptome libraries were sequenced using an Illumina HiSeq X Ten sequencer with a paired-end protocol.
Read mapping, normalization, and statistical analysis of differential gene expression
According to the BGI’s instructions, the low quality reads (More than 20% of the bases qualities were lower than 10), reads with adaptors and reads with unknown bases (N bases more than 5%) were filtered using an internal software, SOAPnuke, to get the clean reads. After filtering, the remaining reads were stored in FASTQ format. All reads were aligned to the host genome or HSV-1 17+ using HISAT2 software with default parameters . In viral transcriptome analysis, all reads aligned to host genomes were filtered out, and all unmapped reads were then mapped against HSV-1 17+ genome. The quantification of transcript abundance (raw count) was conducted using featureCounts  software supported by the Subread package , and multi-mapping and multi-overlapping reads were excluded from expression analysis. We chose to normalize our data set using the relative log expression method implemented in the DESeq2 package (v1.22.2)  in R (www.r-project.org), and only viral genes that had a read count of at least 5 were considered. After DESeq2 analyses, genes with a log2 fold change > 1 and FDR- adjusted p-value <− 0.01 were considered differentially expressed.
Reverse transcription PCR (RT-PCR) and quantitative real-time RT-PCR (RT-qPCR)
In Fig. 5a and b, reverse transcription used strand-specific primers for ICP0 and LAT transcript respectively and was performed with a high fidelity RT-PCR kit (TaKaRa), according to the manufacturer’s instructions. In a 20 μl reaction mixtures contained RNase inhibitor, RTase, PrimeScript buffer, dNTP mixture, 1 μg template RNA and 5pM a strand-specific RT primer. After RT reaction, 1 μL cDNA was used for PCR reaction flowing by 30 cycles of 10 s at 98 °C, then 5 s at 60 °C, and 20s at 72 °C. The RT-PCR product was electrophoresed on a 2% agarose gel for 20 min at 200 V. For gene expression analysis via RT-qPCR, isolated total RNA was reverse-transcribed using a PrimeScriptTM RT reagent Kit with gDNA Eraser (TaKaRa). The RT-qPCR was performed in an Applied Biosystems 7900HT using FastStart Universal SYBR Green Master (ROX, Roche). Per RT-qPCR reaction, cDNA derived from 50 ng RNA was deployed. Cycling conditions were 10 min 95 °C, followed by 45 cycles of 10s at 95 °C, 30s at 60 °C, and 30s at 72 °C. Primers are listed in supplemental Table 2. For the expression of fold changes of RT-qPCR data, the 2 − ΔΔCt method was used .
Availability of data and materials
The HSV-1 reference genome and annotation file were downloaded from NCBI (ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/859/985/GCF_000859985.2_ViralProj15217), GCF_000859985.2_ViralProj15217_genomic.fna.gz and GCF_000859985.2_ViralProj15217_genomic.gff.gz). The RNA-seq data generated for this paper have been submitted to the Genome Sequence Archive  in BIG Data Center , Beijing Institute of Genomics (BIG), Chinese Academy of Sciences (GSA, http://bigd.big.ac.cn/gsa), under accession number CRA001750.
The datasets used and analysed during the current study are available from the corresponding author on reasonable request.
Herpes simplex virus
Herpes simplex virus-1
days post infection
Conservation Genetics Academy of Science
Reverse transcription PCR
Quantitative real-time RT-PCR
Differentially expressed gene
Roizman B, Knipe DM, Whitely R: Herpes Simplex Viruses. In Fields Virology (Knipe, Fields Virology)-2 Volume Set). Edited by David M K, Peter H; 2013: 1823-1897.
Valyi-Nagy T, Deshmane S, Dillner A, Fraser NW. Induction of cellular transcription factors in trigeminal ganglia of mice by corneal scarification, herpes simplex virus type 1 infection, and explantation of trigeminal ganglia. J Virol. 1991;65:4142–52.
Stevens JG, Cook ML. Latent herpes simplex virus in spinal ganglia of mice. Science. 1971;173:843–5.
Cook ML, Bastone VB, Stevens JG. Evidence that neurons harbor latent herpes simplex virus. Infect Immun. 1974;9:946–51.
Li L, Li Z, Wang E, Yang R, Xiao Y, Han H, Lang F, Li X, Xia Y, Gao F, et al. Herpes simplex virus 1 infection of tree shrews differs from that of mice in the severity of acute infection and viral transcription in the peripheral nervous system. J Virol. 2015;90:790–804.
Honess RW, Roizman B. Regulation of herpesvirus macromolecular synthesis: sequential transition of polypeptide synthesis requires functional viral polypeptides. Proc Natl Acad Sci U S A. 1975;72:1276–80.
Honess RW, Roizman B. Regulation of herpesvirus macromolecular synthesis. I. Cascade regulation of the synthesis of three groups of viral proteins. J Virol. 1974;14:8–19.
Umbach JL, Kramer MF, Jurak I, Karnowski HW, Coen DM, Cullen BR. MicroRNAs expressed by herpes simplex virus 1 during latent infection regulate viral mRNAs. Nature. 2008;454:780–3.
Kramer MF, Jurak I, Pesola JM, Boissel S, Knipe DM, Coen DM. Herpes simplex virus 1 microRNAs expressed abundantly during latent infection are not essential for latency in mouse trigeminal ganglia. Virology. 2011;417:239–47.
Stevens JG, Wagner EK, Devi-Rao GB, Cook ML, Feldman LT. RNA complementary to a herpesvirus alpha gene mRNA is prominent in latently infected neurons. Science. 1987;235:1056–9.
Shen W, Sa e Silva M, Jaber T, Vitvitskaia O, Li S, Henderson G, Jones C: Two small RNAs encoded within the first 1.5 kilobases of the herpes simplex virus type 1 latency-associated transcript can inhibit productive infection and cooperate to inhibit apoptosis. J Virol 2009, 83:9131–9139.
Rowe AM, St Leger AJ, Jeon S, Dhaliwal DK, Knickelbein JE, Hendricks RL. Herpes keratitis. Prog Retin Eye Res. 2013;32:88–101.
Fatahzadeh M, Schwartz RA. Human herpes simplex virus infections: epidemiology, pathogenesis, symptomatology, diagnosis, and management. J Am Acad Dermatol. 2007;57:737–63 quiz 764-736.
Eimer WA, Vijaya Kumar DK, Navalpur Shanmugam NK, Rodriguez AS, Mitchell T, Washicosky KJ, Gyorgy B, Breakefield XO, Tanzi RE, Moir RD. Alzheimer's disease-associated beta-amyloid is rapidly seeded by Herpesviridae to protect against brain infection. Neuron. 2018;99:56–63 e53.
Ezzat K, Pernemalm M, Palsson S, Roberts TC, Jarver P, Dondalska A, Bestas B, Sobkowiak MJ, Levanen B, Skold M, et al. The viral protein corona directs viral pathogenesis and amyloid aggregation. Nat Commun. 2019;10:2331.
Darai G, Schwaier A, Komitowski D, Munk K. Experimental infection of Tupaia belangeri (tree shrews) with herpes simplex virus types 1 and 2. J Infect Dis. 1978;137:221–6.
Darai G, Koch HG, Flugel RM, Gelderblom H. Tree shrew (Tupaia) herpesviruses. Dev Biol Stand. 1982;52:39–51.
Fan Y, Huang ZY, Cao CC, Chen CS, Chen YX, Fan DD, He J, Hou HL, Hu L, Hu XT, et al. Genome of the Chinese tree shrew. Nat Commun. 2013;4:1426.
Li L, Li Z, Li X, Wang E, Lang F, Xia Y, Fraser NW, Gao F, Zhou J. Reactivation of HSV-1 following explant of tree shrew brain. J Neuro-Oncol. 2016;22:293–306.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Risso D, Schwartz K, Sherlock G, Dudoit S. GC-content normalization for RNA-Seq data. BMC Bioinformatics. 2011;12:480.
Watson RJ, Clements JB. A herpes simplex virus type 1 function continuously required for early and late virus RNA synthesis. Nature. 1980;285:329–30.
DeLuca NA, Schaffer PA. Activation of immediate-early, early, and late promoters by temperature-sensitive and wild-type forms of herpes simplex virus type 1 protein ICP4. Mol Cell Biol. 1985;5:1997–2008.
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:2339–49.
Dremel SE, DeLuca NA. Genome replication affects transcription factor binding mediating the cascade of herpes simplex virus transcription. Proc Natl Acad Sci U S A. 2019;116:3734–9.
DeLuca NA, Courtney MA, Schaffer PA. Temperature-sensitive mutants in herpes simplex virus type 1 ICP4 permissive for early gene expression. J Virol. 1984;52:767–76.
Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26:136–8.
Farrell MJ, Dobson AT, Feldman LT. Herpes simplex virus latency-associated transcript is a stable intron. Proc Natl Acad Sci U S A. 1991;88:790–4.
Krummenacher C, Zabolotny JM, Fraser NW. Selection of a nonconsensus branch point is influenced by an RNA stem-loop structure and is important to confer stability to the herpes simplex virus 2-kilobase latency-associated transcript. J Virol. 1997;71:5849–60.
Peng W, Vitvitskaia O, Carpenter D, Wechsler SL, Jones C. Identification of two small RNAs within the first 1.5-kb of the herpes simplex virus type 1-encoded latency-associated transcript. J Neuro-Oncol. 2008;14:41–52.
Feldman LT, Ellison AR, Voytek CC, Yang L, Krause P, Margolis TP. Spontaneous molecular reactivation of herpes simplex virus type 1 latency in mice. Proc Natl Acad Sci U S A. 2002;99:978–83.
Margolis TP, Elfman FL, Leib D, Pakpour N, Apakupakul K, Imai Y, Voytek C. Spontaneous reactivation of herpes simplex virus type 1 in latently infected murine sensory ganglia. J Virol. 2007;81:11069–74.
Wagner EK, Flanagan WM, Devi-Rao G, Zhang YF, Hill JM, Anderson KP, Stevens JG. The herpes simplex virus latency-associated transcript is spliced during the latent phase of infection. J Virol. 1988;62:4577–85.
Wagner EK, Devi-Rao G, Feldman LT, Dobson AT, Zhang YF, Flanagan WM, Stevens JG. Physical characterization of the herpes simplex virus latency-associated transcript in neurons. J Virol. 1988;62:1194–202.
Thompson RL, Sawtell NM. The herpes simplex virus type 1 latency-associated transcript gene regulates the establishment of latency. J Virol. 1997;71:5432–40.
Perng GC, Slanina SM, Yukht A, Ghiasi H, Nesburn AB, Wechsler SL. The latency-associated transcript gene enhances establishment of herpes simplex virus type 1 latency in rabbits. J Virol. 2000;74:1885–91.
Ahmed M, Lock M, Miller CG, Fraser NW. Regions of the herpes simplex virus type 1 latency-associated transcript that protect cells from apoptosis in vitro and protect neuronal cells in vivo. J Virol. 2002;76:717–29.
Carpenter D, Hsiang C, Jiang X, Osorio N, BenMohamed L, Jones C, Wechsler SL. The herpes simplex virus type 1 (HSV-1) latency-associated transcript (LAT) protects cells against cold-shock-induced apoptosis by maintaining phosphorylation of protein kinase B (AKT). J Neuro-Oncol. 2015;21:568–75.
Tormanen K, Allen S, Mott KR, Ghiasi H. The Latency-Associated Transcript Inhibits Apoptosis via Downregulation of Components of the Type I Interferon Pathway during Latent Herpes Simplex Virus 1 Ocular Infection. J Virol. 2019;93(10):e00103–19.
Kent JR, Kang W, Miller CG, Fraser NW. Herpes simplex virus latency-associated transcript gene function. J Neuro-Oncol. 2003;9:285–90.
Nicoll MP, Hann W, Shivkumar M, Harman LE, Connor V, Coleman HM, Proenca JT, Efstathiou S. The HSV-1 latency-associated transcript functions to repress latent phase lytic gene expression and suppress virus reactivation from latently infected neurons. PLoS Pathog. 2016;12:e1005539.
Leib DA, Bogard CL, Kosz-Vnenchak M, Hicks KA, Coen DM, Knipe DM, Schaffer PA. A deletion mutant of the latency-associated transcript of herpes simplex virus type 1 reactivates from the latent state with reduced frequency. J Virol. 1989;63:2893–900.
Leib DA, Nadeau KC, Rundle SA, Schaffer PA. The promoter of the latency-associated transcripts of herpes simplex virus type 1 contains a functional cAMP-response element: role of the latency-associated transcripts and cAMP in reactivation of viral latency. Proc Natl Acad Sci U S A. 1991;88:48–52.
Mador N, Braun E, Haim H, Ariel I, Panet A, Steiner I. Transgenic mouse with the herpes simplex virus type 1 latency-associated gene: expression and function of the transgene. J Virol. 2003;77:12421–9.
Trousdale MD, Steiner I, Spivack JG, Deshmane SL, Brown SM, MacLean AR, Subak-Sharpe JH, Fraser NW. In vivo and in vitro reactivation impairment of a herpes simplex virus type 1 latency-associated transcript variant in a rabbit eye model. J Virol. 1991;65:6989–93.
Bloom DC, Devi-Rao GB, Hill JM, Stevens JG, Wagner EK. Molecular analysis of herpes simplex virus type 1 during epinephrine-induced reactivation of latently infected rabbits in vivo. J Virol. 1994;68:1283–92.
Perng GC, Dunkel EC, Geary PA, Slanina SM, Ghiasi H, Kaiwar R, Nesburn AB, Wechsler SL. The latency-associated transcript gene of herpes simplex virus type 1 (HSV-1) is required for efficient in vivo spontaneous reactivation of HSV-1 from latency. J Virol. 1994;68:8045–55.
Perng GC, Ghiasi H, Slanina SM, Nesburn AB, Wechsler SL. The spontaneous reactivation function of the herpes simplex virus type 1 LAT gene resides completely within the first 1.5 kilobases of the 8.3-kilobase primary transcript. J Virol. 1996;70:976–84.
Watson ZL, Washington SD, Phelan DM, Lewin AS, Tuli SS, Schultz GS, Neumann DM, Bloom DC. In Vivo Knockdown of the Herpes Simplex Virus 1 Latency-Associated Transcript Reduces Reactivation from Latency. J Virol. 2018;92(16):e00812–18.
Hill JM, Nolan NM, McFerrin HE, Clement C, Foster TP, Halford WP, Kousoulas KG, Lukiw WJ, Thompson HW, Stern EM, Bhattacharjee PS. HSV-1 latent rabbits shed viral DNA into their saliva. Virol J. 2012;9:221.
Hill JM, Sedarati F, Javier RT, Wagner EK, Stevens JG. Herpes simplex virus latent phase transcription facilitates in vivo reactivation. Virology. 1990;174:117–25.
Jarman RG, Loutsch JM, Devi-Rao GB, Marquart ME, Banaszak MP, Zheng X, Hill JM, Wagner EK, Bloom DC. The region of the HSV-1 latency-associated transcript required for epinephrine-induced reactivation in the rabbit does not include the 2.0-kb intron. Virology. 2002;292:59–69.
Bloom DC, Hill JM, Devi-Rao G, Wagner EK, Feldman LT, Stevens JG. A 348-base-pair region in the latency-associated transcript facilitates herpes simplex virus type 1 reactivation. J Virol. 1996;70:2449–59.
Kang W, Mukerjee R, Gartner JJ, Hatzigeorgiou AG, Sandri-Goldin RM, Fraser NW. Characterization of a spliced exon product of herpes simplex type-1 latency-associated transcript in productively infected cells. Virology. 2006;356:106–14.
Devi-Rao GB, Goodart SA, Hecht LM, Rochford R, Rice MK, Wagner EK. Relationship between polyadenylated and nonpolyadenylated herpes simplex virus type 1 latency-associated transcripts. J Virol. 1991;65:2179–90.
Zwaagstra JC, Ghiasi H, Slanina SM, Nesburn AB, Wheatley SC, Lillycrop K, Wood J, Latchman DS, Patel K, Wechsler SL. Activity of herpes simplex virus type 1 latency-associated transcript (LAT) promoter in neuron-derived cells: evidence for neuron specificity and for a large LAT transcript. J Virol. 1990;64:5019–28.
Mitchell WJ, Lirette RP, Fraser NW. Mapping of low abundance latency-associated RNA in the trigeminal ganglia of mice latently infected with herpes simplex virus type 1. J Gen Virol. 1990;71(Pt 1):125–32.
Chen SH, Lee LY, Garber DA, Schaffer PA, Knipe DM, Coen DM. Neither LAT nor open reading frame P mutations increase expression of spliced or intron-containing ICP0 transcripts in mouse ganglia latently infected with herpes simplex virus. J Virol. 2002;76:4764–72.
Maillet S, Naas T, Crepin S, Roque-Afonso AM, Lafay F, Efstathiou S, Labetoulle M. Herpes simplex virus type 1 latently infected neurons differentially express latency-associated and ICP0 transcripts. J Virol. 2006;80:9310–21.
Chen Q, Lin L, Smith S, Huang J, Berger SL, Zhou J. CTCF-dependent chromatin boundary element between the latency-associated transcript and ICP0 promoters in the herpes simplex virus type 1 genome. J Virol. 2007;81:5192–201.
Beishline K, Vladimirova O, Tutton S, Wang Z, Deng Z, Lieberman PM. CTCF driven TERRA transcription facilitates completion of telomere DNA replication. Nat Commun. 2017;8:2114.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.
Liao Y, Smyth GK, Shi W. The subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 2013;41:e108.
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:402–8.
Wang Y, Song F, Zhu J, Zhang S, Yang Y, Chen T, Tang B, Dong L, Ding N, Zhang Q, et al. GSA: genome sequence archive<sup/>. Genomics Proteomics Bioinformatics. 2017;15:14–8.
Members BIGDC. The BIG data center: from deposition to integration to translation. Nucleic Acids Res. 2017;45:D18–24.
Council. NR. Guide for the care and use of laboratory animals. Washington, DC: National Academy Press; 1996.
The authors would like to thank Yangyang YU at BGI-Shenzhen for guidance of bioinformatics analysis.
This study was supported by grants from National Natural Science Foundation of China-Yunnan Joint Found (NSFC, U1602226), National Natural Science Foundation of China (NSFC, 81672040) and Ministry of Science and Technology of China (MOST, SQ2018YFE020419) and a Thousand Foreign Talent scholarship from Yunnan province to Jumin Zhou, Yunnan applied basic research project (2017FE467) and International Science and Technology Cooperation Project (2017IB011) to Xia Cao and Chinese Academy of Sciences President’s International Fellowship Initiative (PIFI, 2019VBA0045) to Nigel W. Fraser.
Ethics approval and consent to participate
All tree shrew and mouse experiments were performed according to the procedures and policies of the Guide for the Care and Use of Laboratory Animals . The Study protocol of ocular scarification and ocular HSV-1 infection was reviewed and approved by the Institutional Animal Care and Use Committee of the Kunming Institute of Zoology of the Chinese Academy of Sciences. The identification number is SMKX-20150317-65. All of the experimental protocols involving human trigeminal ganglia were approved by the Human Ethical Commission at Kunming Medical University and the Ethical Commission of the Academy of Forensic Science (Shanghai, China). All of the donors’ parental or legal guardians provided written informed consent.
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, E., Ye, Y., Zhang, K. et al. Longitudinal transcriptomic characterization of viral genes in HSV-1 infected tree shrew trigeminal ganglia. Virol J 17, 95 (2020). https://doi.org/10.1186/s12985-020-01344-8