Correction: Lytic Gene Expression Is Frequent in HSV-1 Latent Infection and Correlates with the Engagement of a Cell-Intrinsic Transcriptional Response

Herpes simplex viruses (HSV) are significant human pathogens that provide one of the best-described examples of viral latency and reactivation. HSV latency occurs in sensory neurons, being characterized by the absence of virus replication and only fragmentary evidence of protein production. In mouse models, HSV latency is especially stable but the detection of some lytic gene transcription and the ongoing presence of activated immune cells in latent ganglia have been used to suggest that this state is not entirely quiescent. Alternatively, these findings can be interpreted as signs of a low, but constant level of abortive reactivation punctuating otherwise silent latency. Using single cell analysis of transcription in mouse dorsal root ganglia, we reveal that HSV-1 latency is highly dynamic in the majority of neurons. Specifically, transcription from areas of the HSV genome associated with at least one viral lytic gene occurs in nearly two thirds of latently-infected neurons and more than half of these have RNA from more than one lytic gene locus. Further, bioinformatics analyses of host transcription showed that progressive appearance of these lytic transcripts correlated with alterations in expression of cellular genes. These data show for the first time that transcription consistent with lytic gene expression is a frequent event, taking place in the majority of HSV latently-infected neurons. Furthermore, this transcription is of biological significance in that it influences host gene expression. We suggest that the maintenance of HSV latency involves an active host response to frequent viral activity. Copyright: ß 2014 Ma et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.


Introduction
Herpes simplex virus (HSV) is a ubiquitous human pathogen that initially replicates in skin and mucosal surfaces prior to establishing an infection in sensory neurons [1]. There, HSV persists life-long with periods of prolonged quiescence punctuated by occasional episodes of reactivation, which may result in recrudescence of disease. As a consequence, HSV has long been viewed as establishing the prototypic latent infection, with minimal antigen expression and little to no active virus production during this phase of infection. While in humans, latency can be short-lived with multiple and frequent bouts of virus reactivation [2], in mice it is especially profound [3,4], making them the ideal model for detailed examination of this feature. Despite this general stability, one part of the HSV genome has long been known to be active in at least a fraction of neurons, giving rise to the expression of noncoding latency-associated transcripts (LAT) [5,6]. These transcripts have been postulated to play various roles in latency and this region is also the source of miRNAs that may also modulate infection [7][8][9][10][11][12][13]. In addition, there is sporadic, but consistent evidence that the transcription of genes normally associated with lytic replication also occurs in sensory ganglia during HSV latency, albeit at far lower levels than LAT [14][15][16][17][18][19][20][21][22]. There are competing explanations for the presence of these lytic transcripts. On the one hand, the lytic products may simply represent sporadic occurrences of complete reactivation, but at frequencies too low to permit ready detection of infectious virus or result in recrudescent disease [23]. An alternative proposal is that lytic transcripts are more widespread, found in a high proportion of infected cells during latency, but they represent abortive reactivation events that fail to progress to productive infection due to additional host control mechanisms [14,24]. Finally it is possible that these transcripts are simply the result of less than total silencing of the HSV genome, but are of no wider biological consequence.
The key to understanding the meaning of non-LAT transcription in HSV latency lies in the ability to examine events at the level of individual neurons. Recent advances in single-cell approaches to transcriptional analysis increase the sensitivity of detecting gene expression and allow the filtering out of samples containing transcripts from irrelevant cell types [25,26]. This means it is now possible to profile gene expression of individual cells using quantitative PCR, as reported in various studies [27][28][29]. One challenge in the application of these technologies to HSV latency is the identification and isolation of latently infected neurons from those that are uninfected and other cell types in sensory ganglia [30]. In this study, we exploited HSV-borne Cre recombinase to indelibly mark infected neurons in order to subject them to transcriptional analysis after laser capture microdissection (LCM) [31]. In doing so, we found that the presence of non-LAT viral transcripts, while at low levels, are detected very frequently in latently-infected neurons. Importantly, infected neurons respond to and match the viral activity with increasing host anti-viral gene expression, demonstrating that this viral transcription has biological consequences. These results demonstrate that the apparent stability of latency in virological terms hides a far more dynamic process in individual neurons.

Presence of low level lytic transcripts during HSV latency in whole ganglia
Mice were inoculated by flank scarification using HSV type 1 (HSV-1) strain KOS, which results in the infection of the skin and innervating dorsal root ganglia (DRG) [32]. Here we define latency as a lack of virus production, and it is in place after the elimination of infectious virus in the DRG from around 7-8 days after inoculation ( Figure 1A). We confirmed the presence of HSV DNA ( Figure 1B) and 2 kb LAT intron ( Figure 1C) in DRG at latency. Expression of lytic genes during latency has been reported by some [14][15][16][17], but not other studies [5,6,13,[33][34][35][36] and so the meaning of the positive reports remains in question. Recent reviews tend to focus on the mechanisms by which lytic genes are silenced during latency, rather than questioning their possible importance [37,38]. Taking a conventional quantitative RT-PCR approach, we found LAT, but no transcripts consistent with lytic gene expression, during latency when whole DRG were examined ( Figure 1C), in contrast to what was seen during the acute period of infection (days 50 and 5 post-inoculation, respectively). However, when a pre-amplification step was added to our protocol (see materials and methods), we were able to detect transcripts from all three tested lytic loci, namely UL54 (ICP27), UL29 (ICP8) and US6 (gD) ( Figure 1D).
In vivo labeled-infected neurons are numerically stable during latency In order to move the analysis to single neurons, we employed a transgenic mouse model in which all infected neurons are marked so that they can be distinguished from other cells [39,40]. Transgenic mice having either yellow fluorescent protein (YFP) or b-galactosidase (LacZ) genes driven by the ubiquitous ROSA26 promoter interrupted with a floxed neomycin cassette were infected with HSVs expressing Cre recombinase under the CMV IE promoter (KOS0152 and KOS/pCMV/eGC, respectively) [40]. Transient expression of Cre occurs in every infected cell, regardless of fate, thereby leaving a legacy of marker gene expression that lasts for the life of the animal [39,40]. In our model, once the acute infection has subsided, DRG neurons marked in this way declined until day 20 and remained numerically stable until day 100, the latest time examined (Figure 2A). Others who have examined similar models in mice with the same methods have extended the observation of stable latency out to 147 days [39]. DRG from KOS0152-infected ROSA-YFP mice, but not mice infected with WT HSV or uninfected mice, contained fluorescence detectable by two-photon microscopy ( Figure 2B) or in microscopy used for laser-dissection ( Figure 2C) consistent with YFP expression in neurons. Samples examined by two-photon microscopy ( Figure 2B) also showed an auto-fluorescent background consisting of scattered small bright particles. Thus, this method showed that once latency was established, the number of marked cells remained constant over a prolonged period and likely for the life of the infected mouse. More importantly, it provided us with a means of determining whether lytic transcripts could be detected in individual infected neurons and if so, determining the frequency and consequence of such expression for host gene expression.

Validation of a method for single cell analysis of transcription in HSV latency
Our approach to obtaining single neuron samples was to use laser capture microscopy (LCM) to identify by YFP expression and cut relevant cells from 8 mm sections of pooled latent DRG ( Figure 3A). Marking of neurons in HSV-infected DRG was confirmed by co-localization of YFP expression with NeuroTrace staining, which detects Nissl bodies ( Figure 3B, arrows). As a negative control, we never found YFP + neurons in DRG from uninfected mice. While the method worked well for the identification of latently-infected neurons, it had limitations; including the likelihood of partial cell capture, variations in capture volume, and inter-cell contamination. We reasoned that the impact of such limitations would be minimized if we examined many neuronal sections (hereafter referred to simply as neurons for brevity) in each analysis. In addition, variations associated with partial cell capture were expected to be consistent across all groups, allowing valid comparisons that would be similar to what could be found if whole neurons were to be used.
To test the reliability of the method, some preliminary studies were done. The first was to extract DNA from neurons and use qPCR to determine levels of HSV DNA in YFP + and YFP 2 neurons. Most individual YFP + cells contained HSV DNA, whereas more than 80% of those with no detectable YFP did not ( Table 1). The detection of HSV in YFP 2 cells was most likely because of incomplete recombination of the ROSA-YFP locus, poor expression from the recombined genes or photo-bleaching of the YFP during capture. A proportion of YFP + cells did not show detectable levels of HSV DNA, probably because we were unable to detect fewer than 10 copies reliably, but also possibly due to

Author Summary
Primary herpes simplex virus (HSV) infections are characterized by acute disease that resolves rapidly, but the virus persists in a latent form in sensory neurons that can be a source of renewed disease. Analyzing gene expression in single mouse neurons harboring latent HSV, we show directly that HSV latency is dynamic and heterogeneous. HSV lytic gene transcripts were frequently detected in latently infected neurons and often in combinations. Expression of selected cellular anti-viral and survival genes showed that transcriptional profiles differed between latently infected and uninfected neurons from the same ganglia. The pattern of host gene expression also differed between latently infected neurons that were and were not experiencing HSV lytic gene expression. Our study suggests that HSV latency is characterized by very frequent switching on of lytic genes and a rapid response by the host, presumably to halt progression to reactivation.
incomplete cell capture. Encouragingly, the distribution of HSV genome copy numbers per neuron above the detection threshold was consistent with previous studies [30,41]. The bulk of YFP + neurons contained fewer than 500 copies of viral DNA, with only a few containing more than 1000 copies ( Figure 3C, Table 2).
To enable analysis of multiple transcripts on a single cell basis, we chose a highly sensitive method that committed each sample in its entirety to a single, random-primed reverse transcription (RT) followed by a multiplexed pre-amplification step before splitting into 61 qPCR (Taqman) reactions. This method has some caveats, for example we cannot verify the strand of transcripts (see discussion). This would also leave no material for controls, such as a sham RT. To ensure our extraction method was not frequently subject to DNA contamination, a set of 16 YFP + neurons was used in a test where each was split into two, enabling a sham RT control. We off-set the reduction in starting nucleic acid by using twice the usual amount of pre-amplified cDNA to quantify LAT. We reasoned that LAT should be detected in at least a proportion of neurons and because viral DNA contamination was the greatest concern. As shown as a heat map in Figure S1, LAT was detected in 15 out of 16 neurons by regular reactions, but none from the sham RT controls. This suggests that DNA contamination is infrequent if not absent in single neuron RNA samples prepared using our chosen method.
For the main analysis of transcription our study samples initially comprised RNA from 205 YFP + and 256 YFP 2 neurons from DRG of mice taken 50 days after infection (combined from four independent infections) and 136 neurons from DRG of uninfected mice (from three independent experiments). Each of these single cell RNA samples was analyzed with 61 primer/probe sets, targeting 51 cellular and 10 viral transcripts. Of the cellular transcripts, 47 of the 51 primers/probes spanned exon-exon boundaries. We included markers that allowed application of stringent selection criteria to neurons for further analyses. These included Gfap and Rbfox3 as markers for satellite cells and neurons with adequate RNA, respectively. To eliminate the impact of possible errors in capturing neurons by LCM and spurious detections, only Gfap 2 , Rbfox3 + neurons in which more than 10% of the host genes could be detected were analysed further. In addition, uninfected neurons (YFP 2 or from uninfected DRG) in which any viral transcripts were detected were also excluded. Numbers of neurons excluded for any of these reasons and those remaining in analyses are noted in Table S1. The following analyses and results were all derived from this main data set.

Ntrk1 + neurons harbor HSV during latency
Within the 51 cellular gene targets, 6 genes (Ntrk1, Ntrk2, Ntrk3, Runx1, Runx3 and Ret) characterized different neuronal lineages. Analysis of these 6 genes showed that infection was heavily biased towards Ntrk1 + neurons ( Figure S2), since Ntrk1 gene (encoding TrkA protein) expression was significantly overrepresented in YFP + neurons: 86% of the YFP + neurons in the infected DRG were Ntrk1 + compared to only 34% YFP 2 neurons ( Figure 3D). Combined, the results demonstrated that YFP only marked infected neurons, predominantly of Ntrk1 + peptidergic nociceptors [42]. Like the HSV DNA copy number, this result is consistent with the literature and suggests that our approach leads  to reliable data [43,44]. Further, this finding led us to narrow the rest of our overall analyses on Ntrk1 + neurons to avoid finding transcriptional signatures associated with neuron subsets and not latency. At the same time, we note that Ntrk1 + neurons are not homogenous and further single cell approaches may reveal details of HSV tropism within this population.

Frequent viral lytic gene expression in the majority of individual infected neurons during latency
Several studies comparing presence of HSV genomes by PCRbased methods and LAT expression in rodent sensory ganglia by conventional in situ hybridization found that only a subset of infected neurons expressed LAT during latency [41,[45][46][47]. However, application of in situ RT-PCR suggests that LAT expression is more widespread [48] and this was also found in our RT controls ( Figure S1). In the main data set, we found RNA consistent with 2 kb LAT in approximately 90% of latentlyinfected neurons ( Figure 4A, 4B and 4D). Comparing this with the frequency of errors associated with our LCM method suggests strongly that all latently-infected neurons express LAT at levels detectable by our method. We therefore excluded YFP + LAT 2 samples from analyses to ensure that we were not including neurons that were not latently-infected.
Differences in sensitivity might also underlie variable accounts of lytic gene expression in latency (e.g. Figure 1C and D). We were able to detect transcripts consistent with lytic gene expression (called lytic transcripts hereafter, see discussion) in single cells and this was very frequent, occurring in 59% of all YFP + neurons ( Figure 4C). When the analysis was narrowed to Ntrk1 + LAT + samples, which we can more confidently state are latently-infected neurons, the frequency rose to 64% ( Figure 4E). Also of note was that many neurons had HSV transcripts from many lytic gene loci. HSV lytic genes are classified into 3 major kinetic classes, namely immediate early (IE), early (E) and late (L). To reduce the complexity in all the possible combinations of the lytic genes in our experiments, we assigned neurons to one of three groups: i) infected (YFP + ) neurons that did not express detectable lytic genes (non-lytic); ii) those that expressed lytic genes from 1 or 2, but not all 3, classes (partial-lytic); iii) those that expressed multiple gene transcripts with at least 1 member of each lytic class (full-lytic). Figure 4E shows that LAT + YFP + neurons were roughly divided into thirds among the 3 categories. So there were roughly as many latently-infected neurons with transcripts consistent with the expression of lytic genes from all kinetic classes as there were that had none of these RNAs. Therefore even at times when latency appears to be highly stable, HSV-infected neurons express non-LAT transcripts at a frequency that has previously only been demonstrated to occur in lytic infection.

Infected neurons respond to latent virus by modulating host gene expression
It has been shown that at the level of the whole ganglia, latent infection alters transcriptional profiles, consistent with an antiviral response working to limit virus reactivation [49]. However, a study of this design cannot dissect changes in neurons from those associated with other cell types, including immune infiltrates and it is possible that non-infected neurons within ganglia respond to infection in a bystander manner. These issues were addressed by an analysis of the anti-viral and cell-survival gene transcripts included in our qPCR panel. The processed gene expression data are shown as a heat map and the differences in the mean expression level between populations have been calculated ( Figure  S3 and Table S2). To reduce the complexity of this data, principal components (PC) analysis was used on 48 cellular genes in these Ntrk1 + neurons. PC1, 2 and 3 captured 44.27%, 6.73% and 5.22% of the variation in the gene expression data, respectively. Most YFP + neurons segregated from those that were YFP 2 and neurons from uninfected DRG on PC2 and PC3, with some overlap captured on PC1 ( Figure 5A). In contrast, YFP 2 neurons from infected mice and neurons from uninfected ganglia were relatively closely bunched together on PC1-3 ( Figure 5B). Another way of looking at the differences between populations was to simply count the number of transcripts that differed both in level of expression and frequency of detection (Tables S3 and S4). In the case of infected and uninfected neurons from infected ganglia, 17 of the 48 genes differed. However, only a single transcript differed

Infected neurons match increasing viral activity by increasing anti-viral and survival gene expression
We reasoned that if the HSV lytic gene expression observed in latently-infected neurons was biologically relevant, it might influence the expression of host genes. For this reason, we compared host gene expression across the three groups of latentlyinfected neurons as identified by expression of HSV lytic genes (non-, partial-and full-lytic). To do this, we visualized the transcriptional data as violin plots (Figures 6A and S4). This representation is essentially a mirrored histogram that combines the transcript levels for each cellular gene (y-axis, represented by Log 2 Ex) and the proportion of individual cells distributed along a given expression value (x-axis). The value of this representation can be seen in that most genes showed a bimodal distribution, indicating that the population contained individual cells with either high or low levels of a given gene transcript. We also calculated differences in mean transcript levels (Table S2) and whether level or frequency of detection was significantly different amongst the three populations of latently-infected neurons (Tables S3 and S4, respectively). We found that more than half (29/48) of the transcripts differed significantly both by level and frequency of detection between full-versus non-lytic groups, with fewer differences between full-versus partial-and partial-versus nonlytic neurons (8/48 and 1/48, respectively). Examples of some of these differences include: 1) Oasl2, Mx1, Ddx58, SerpinB9, Bcl2l11, Pias1 and Pvrl1 were expressed at significantly higher frequencies and levels in the full-lytic subgroup of neurons compared to the non-lytic subset. 2) Transcripts of Tnfrsf10b, Bcl2 and Tnfrsf14 were over-represented in full-lytic when compared to both of the other groups.
3) The single gene that differed between non-lytic and the groups that expressed lytic genes was Pou2f1. These data revealed heterogeneity in host gene expression in latently infected neurons that correlated with the number of classes of viral lytic genes detected. We conclude that transcription of HSV lytic genes during HSV latency is biologically relevant and suggest that even in the absence of overt reactivation, there is a progressive neuronal response against increasing viral activity.

Viral transcription activity alters neuronal transcriptional circuits during latency
Genes involved in the same biological processes or pathways are often expressed with correlated profiles. We wondered whether we could identify such sets of genes that were co-regulated in response to viral lytic gene transcription in latently-infected neurons. Finding these sets would provide evidence of changes to broader transcriptional circuits. To do this, we examined the correlations in expression profiles between pairs of cellular genes from full-and non-lytic neurons. In the full-lytic group, strong gene correlations (t.0.6) were found between 75 pairs of genes (Tables S5 and S7). Importantly, the correlations were weaker and less frequent (only 23 pairs with t.0.6) in the non-lytic group (Tables S6 and S8). Correlations did not only exist between pairs of genes but were also observed between 1 gene and several other genes. To show the co-dependencies between groups of genes, we clustered genes based on their correlation coefficients. First, we identified clusters of genes that were strongly correlated in full-lytic neurons and showed that these same clusters were less prominent in non-lytic neurons ( Figure 6B, Table 3). Separately, we determined the optimal clustering of gene correlations in either full-lytic versus non-lytic neurons and this showed that viral lytic gene expression altered the optimal clustering patterns ( Figure 6C). Thus, both comparisons found that sets of host transcripts and presumably the underlying regulatory mechanisms were altered in neurons where lytic gene loci were transcribed.

Discussion
Long-term HSV infection in mice is generally viewed as a model of profound latency, with few, if any, of the spontaneous reactivation episodes seen in other species and virtually no active production of virus once the lytic phase of infection is extinguished. In the context of this model of stable latency, the frequency with which transcripts consistent with lytic gene expression were found here is striking. Whilst our data does not break ground in terms of the detection of viral transcripts from outside the LAT region during latency, that this might happen so often has previously been only a matter of conjecture [14]. Perhaps more importantly, we also show that the host is equally active in responding to this viral activity, demonstrating that the viral transcription has biological consequences.
Some caveats remain in our study with regard to the nature of the viral transcripts detected. The most important of these is that we were unable to be certain that the transcripts detected were bona fide mRNAs from the various lytic genes. Some regions of the HSV genome are transcribed to produce co-linear messages for a set of genes and further, especially at late times of lytic infection, the termination of transcription is not always accurate. For this reason, in some cases we may not have detected only the HSV transcript intended. In addition, our reverse transcription method used random primers, being chosen to allow multiple transcripts to be detected at high sensitivity, but this also sacrifices strand-specificity. However, the nine HSV genes targeted are spread across the entire genome and different patterns of host gene expression are associated with increasing detection of the virus transcripts when ordered by expected kinetic classes. Furthermore, in some cases others have demonstrated elegantly that transcripts genuinely originating from specific HSV lytic genes are found in latently-infected ganglia [14]. Finally, if we are not detecting lytic gene expression, the alternative is that we have uncovered widespread random transcription of the HSV genome in latency that influences the host, which seems an unlikely proposition. Taking these together, the simplest explanation of our data is that we are detecting expression of HSV lytic genes, if perhaps not necessarily only those stated.
The relative roles of virus and host in maintaining latency remains an area of some contention, as is the extent to which detection of HSV transcripts other than LAT is of any significance. Much work has focused on mechanisms by which transcription from the majority of the genome is inhibited or regulated during latency, either through the silencing effect of heterochromatin or expression of miRNAs [1,7,37,50]. Alternatively, evidence of ongoing detection of HSV proteins by CD8 + T cells [51][52][53][54] and the possibility that these can act in a non-lytic manner on neurons [32,55] suggest adaptive immunity might be an active player in the maintenance of latency. Our data do not bear on the first of these possibilities, except to demonstrate that the mechanisms that impede lytic gene transcription in latency substantially reduce, but do not silence this gene expression in the majority of neurons. The finding of frequent residual HSV lytic gene expression might be used to support a role for the adaptive immune system, but there remains the important caveat that we are only detecting RNA, with no evidence of translation. Such protein expression could be at an extremely low level, given the known sensitivity of T cells [56]. Also intriguingly, we found the expression of Serpinb9, which encodes a protein that can protect cells from CD8 + T cell cytotoxicity, was up-regulated in the neurons showing the most viral transcription [57,58]. At the same time, we note that beyond the presence of activated T cells, evidence for viral protein expression in latency is conflicting. On the one hand, work showing that increased IE gene expression in viruses lacking LAT or miRNAs from the LAT region leads to increased transcription of genes from later classes might be interpreted to suggest that in these cases, IE protein in made [12,13,59]. On the other hand, these were not wild type viruses and a highly sensitive historical analysis based on expression of Cre from a variety of HSV lytic gene promoters that would allow low level or brief protein expression to be detected by Cre-reporter mice failed to detect any lytic gene expression [39].
These difficulties and our results could be used to argue for a third model in which the host acts to maintain latency in a cellintrinsic manner. The viral transcription detected here may represent abortive reactivation events [14,24]. Of interest, we found significant expression of Pou2f1 (Oct1), which is important for viral reactivation [60], in the full-and partial-lytic groups of neurons. We suggest these events fail to progress to productive infection due to detection of virus activity and a host response within neurons. This scenario requires that HSV activity in neurons is recognized by cell-intrinsic sensing mechanisms in the absence of translation of the lytic gene transcripts. Such mechanisms are well documented in cells other than neurons. In this regard, our finding that Zbp1 (also referred to as DAI) is more frequently and significantly expressed in latently-infected, compared with uninfected neurons is of interest because this DNA sensor is associated with the response of fibroblasts to HSV [61]. The identity of the genes and possible mechanisms that might lead to aborting reactivation remains obscure, but our data show significant up-regulation of several genes with established anti-viral roles in at least one relevant comparison (e.g. Mx1, Oasl2, Eif2ak2 (PKR), Ifitm3 and Samhd1). Thus, we propose that there is a barrier to productive infection downstream of the initiation of lytic gene transcription and this is a threshold that is rarely exceeded in mice [14]. Perhaps, understanding why the response to lytic gene transcription is so effective in this species, but not in people who suffer frequent reactivation and recurrence of disease, will highlight pathways of the greatest therapeutic potential. Finally, irrespective of the merits of any of these competing explanations, we would suggest that understanding the consequences of residual lytic gene expression will be important in any complete model of HSV latency.

Ethics statement
All experiments were done according to Australian NHMRC guidelines contained within the Australian Code of Practice for the Care and Use of Animals for Scientific Purposes and under approvals given by The University of Melbourne Animal Ethics Committee (1112345.7) and The Australian National University Animal Experimentation Ethics Committee (A2011/015).

Mice
C57BL/6 mice were bred at the Department of Microbiology and Immunology of the University of Melbourne. The gzmBCreERT2/ROSA26EYFP transgenic mice have been described previously [62]. gzmBCreERT2/ROSA26EYFP mice were not given tamoxifen, and will be referred to as ROSA-YFP. B6.129S4-Gt(ROSA)26Sor tm1So /J, also referred to as ROSA26R mice, were from the Australian Phenomics Facility (ANU, Canberra).

Viruses and infections
The HSV-1 strains KOS0152, KOS/pCMV/eGC and wildtype parental strain KOS (WT HSV) were grown and titrated on Vero cells (CSL, Parkville, Australia) as described in [40]. KOS0152 contains the HCMV IE-Cre expression cassette inserted into the UL43 gene of HSV-1. Similarly, KOS/ pCMV/eGC contains HCMV IE promoter, but driving an EGFP-Cre fusion from the intergenic region beween UL3 and UL4. Mice were infected with 1610 6 PFU via flank scarification as previously described [63]. For single cell gene expression analysis, DRG from infected mice were collected day 50 post inoculation. Other time-points were as noted in the text and figure legends.

Laser capture microdissection
Mice were culled and tissues were fixed in situ by perfusing the mice from the left ventricle with 0.5% paraformaldehyde and 20% sucrose made in nuclease free PBS (fixative). DRG (T8-T12) were harvested, immersed in the fixative, and allowed to equilibrate on ice. Fixed DRG from two mice were pooled and then embedded in OCT (TissueTekIA018; Sakura) and frozen in liquid nitrogen. Frozen DRG were sectioned into 8 mm slices on a cryostat (Leica CM3050S) and consecutive DRG sections were placed on each subsequent slide (LOMB) that has been pre-treated with RNase-Away (Life Technologies). Freshly prepared tissue sections were dehydrated with 50%, 75% and 95% ethanol sequentially and airdried just before laser capture. YFP 2 and YFP + cells were identified from 8-24 non-consecutive tissue sections and laser microdissected using a PALM MicroBeam instrument (Carl Zeiss).
Laser capture was performed at the Biological Optical Microscopy Platform, University of Melbourne, Australia. Cells were then captured into collection tubes containing lysis buffer. We avoided collecting cells from 2-9 consecutive sections to prevent repeated sampling from the same cell. In total, 205 YFP + and 256 YFP 2 cells from HSV infected DRG and 136 cells from uninfected DRG were collected from 3-4 independent experiments with two mice per OCT block.

RNA and DNA extraction
Total RNA from LCM-captured single neurons was made into amplified cDNA using the Single Cell-to-CT kit (Life Technologies) according to manufacturer's instructions. Briefly, LCMcaptured samples were lysed and treated with DNAse; and cell lysates were then reverse transcribed. In some samples, single cell lysates were split into halves whereby one half was reverse transcribed and the other was not. Products of RT (or sham RT) were then pre-amplified in a gene specific manner with a pool of 61 Taqman gene expression assays for 14 cycles; then diluted 1:15 to enable use in the same number of Taqman assays (Tables S9 and S10). Because 2 kb LAT is abundantly expressed in latency [1], we limited the 2 kb LAT primers and probe to a working concentration of 30 nM and 50 nM for pre-amplification, respectively. All other primers and probes were at 180 nM and 50 nM as per manufacturer's instructions, respectively.
Total RNA from whole DRG was extracted and DNAse-treated using RNAqueous microkit (Life Technologies) and reverse transcribed into cDNA using Superscript VILO (Life Technologies) that was then diluted 1:40. Aliquots of some samples were not reverse transcribed and served as 'no RT' controls. In some whole DRG samples, pre-amplification was carried out as described for single cell samples, but these pre-amplified products were not diluted (See Figure 1D).
DNA from LCM-captured neurons and whole DRG were extracted using the Picopure DNA extraction kit (Life Technologies). Briefly, single cells and DRG were incubated in proteinase K at 65uC for 3 h before enzyme inactivation at 95uC for 10 min.

Taqman quantitative PCR
Quantitative PCR assays were performed using Taqman gene expression assays (Tables S9 and S10) and Taqman Fast Advanced mastermix on the StepOnePlus real time PCR machine (Life Technologies). In single cell samples that were split into halves (to allow a matched sham RT control), the amount of pre-amplified cDNA used for quantitative PCR was doubled (See Figure S1). In some whole DRG analysis, viral gene expression was normalized to the housekeeping gene, Pgk1. Relative expression calculated using the 2 2DCt method [64] was used only to represent the absence and presence of viral transcripts.
HSV DNA detection in whole DRG and single neurons were performed as previously described [40]. The sequences, which are specific for the HSV thymidine kinase (UL23) gene, are as follows: forward primer 59-TTGTCTCCTTCCGTGTTTCAGTT-39, reverse primer 59-GGCTCCATACCGACGATCTG-39, and Taqman probe 59-FAM-CCATCTCCCGGGCAAACGTGC-MGB-NFQ-39. DNA extracted from uninfected DRG was used as negative control. Average viral copy number was calculated using a standard curve generated using WT HSV or KOS0152 viral DNA (extracted as described for DRG), which was serially diluted such that 5 ml contained 10 0 copies to 10 7 HSV-1 DNA. Our limit of detection for DNA was 10 copies, as we could not detect below 10 copies consistently. In single cell analysis, a duplex reaction containing the primers/FAM-probe for HSV TK gene, and the primers/VIC-probe for b-actin (Actb) gene (Mm00607939_s1 VIC PL, Life Technologies) that detects genomic DNA was used to filter out LCM samples that did not contain any cells. In all single cell quantitative PCR, mock LCM samples, which contain lysis buffer but no cells, were used as negative controls.
Intravital two-photon microscopy DRG were fixed in situ with 2% paraformaldehyde and 20% sucrose, harvested, secured onto a 10 mm dish with Vetbond tissue adhesive (3M), and immersed in PBS. Images were acquired with an upright LSM710 NLO multiphoton microscope (Carl Zeiss) as described previously [65]. EYFP was excited at 920 nm and collagen was visualized by second harmonic generation. Raw imaging data were processed with Imaris 7.1 (Biplane).
Histology, immunofluorescence and confocal microscopy DRG from HSV infected ROSA26 mice were fixed and stained with X-Gal (BioVectra, Canada) as previously described [40]. DRG were then examined under light microscope (Olympus, Japan) and images of the DRG were taken with an Olympus DP20 camera. The number of b-galactosidase + cells was counted from the images with the aid of ImageJ software.
DRG were fixed in situ with and equilibrated in 2% paraformaldehyde and 20% sucrose before embedding in OCT and frozen in liquid nitrogen. DRG were sectioned at 8 mm thickness. Sections were fixed in 220uC acetone for 5 min, rehydrated with PBS for 5 min, stained with NeuroTrace (N-21482, Life Technologies) at 1:200 v/v in PBS for 20 min, rinsed in PBS as per manufacturer's instructions, incubated with Hoeschst nuclear stain (H33258, 1:3000 vol/vol in PBS) for 3 min, and then mounted with ProLongGold (P36934, Life Technologies). Images were acquired with a Zeiss LSM710 microscope and processed using Imaris 7.1 software (Bitplane).

Single-cell data processing
All Ct values obtained were converted into relative expression threshold (Et) by subtracting the values from the maximum cycle value of 40. When the value of the gene expression is undetermined, we treat this gene as not expressed and assigned an Et value '0' to the gene [66]. No normalization was performed for single cell analyses, which is the standard for this approach, because of large cell-to-cell variances in housekeeping gene expression [67,68] and because of the large biological sample size (e.g. ,200 neurons of each type). Cells that were negative for Rbfox3 (neuron marker), were positive for Gfap (satellite cell marker), and expressed less than 10% of the 51 cellular genes tested were removed from analysis. We consider the remaining Rbfox3 + Gfap 2 cells as neurons. As noted above, transcription of genes follows intermittent bursts of biosynthesis [69], so some neurons may not express Rbfox3 at the time of fixation and these would be wrongly excluded by our approach, but we felt that it was better to set conservative criteria for keeping samples in the final analyses. In a small fraction of YFP 2 neurons from infected DRG (27.4%) and uninfected neurons from uninfected DRG (8.5%), HSV transcripts were detected and these cells were excluded from subsequent analysis.

Single-cell data visualization
Principal components analysis was performed on 48 cellular genes using XLSTAT (Addinsoft SARL, US). Viral genes, Rbfox3, Gfap and Ntrk1 were excluded because single cells were categorized based on these genes. Eigenvectors were then visualized with XLSTAT-3Dplot (Addinsoft SARL, US). Singular analysis toolset (Fluidigm, US) was used to generate violin plots [27]. Gene expression data and correlation coefficients were visualized as heatmaps using TM4 MultiExperimental Viewer [12].

Statistical analysis
Categorical variables were summarized using frequency and percentage. Continuous variables were first assessed for skew using a Shapiro-Wilk test and summarized using mean and standard deviation (SD) or median and inter-quartile range (IQR) as appropriate. Comparison of gene expression levels by neuron and lytic groups were undertaken using a Kruskal-Wallis test with a posthoc Dunn test to correct for multiple simultaneous comparisons. As a sensitivity analysis, gene expression variables were dichotomized into expression recorded/no expression recorded quantities and compared across neuron or lytic groups using a chi-square test or a Fisher's exact test as appropriate. All analyses were conducted using Stata version 12 (StataCorp, College Station, Texas, US).
Correlation coefficients were calculated based on gene expression levels using Kendall rank correlation in XLSTAT (Addinsoft SARL, US). The correlation coefficients were then clustered using k-means or complete-linkage. For k-means clustering, 4 clusters were chosen from the full-lytic YFP + neurons subset based on the elbow method by comparing the sum of squared error (SSE) of a number of clusters. k-means clustering, complete-linkage clustering and SSE were performed in R statistical program.