Japanese encephalitis virus infection induces changes of mRNA profile of mouse spleen and brain

Background Japanese encephalitis virus (JEV) is a mosquito-borne flavivirus, leading to an acute encephalitis and damage to the central nervous system (CNS). The mechanism of JEV pathogenesis is still unclear. DNA microarray analyses have been recently employed to detect changes in host gene expression, which is helpful to reveal molecular pathways that govern viral pathogenesis. In order to globally identify candidate host genes associated with JEV pathogenesis, a systematic mRNA profiling was performed in spleens and brains of JEV-infected mice. Results The results of microarray analysis showed that 437 genes in spleen and 1119 genes in brain were differentially expressed in response to JEV infection, with obviously upregulated genes like pro-inflammatory chemokines and cytokines, apoptosis-related proteases and IFN inducible transcription factors. And the significant pathways of differentially expressed genes are involved in cytokine-cytokine receptor interaction, natural killer cell mediated cytotoxicity, antigen processing and presentation, MAPK signaling, and toll-like receptor signaling, etc. The differential expression of these genes suggests a strong antiviral response of host but may also contribute to the pathogenesis of JEV resulting in encephalitis. Quantitative RT-PCR (RT-qPCR) assay of some selected genes further confirmed the results of microarray assay. Conclusions Data obtained from mRNA microarray suggests that JEV infection causes significant changes of mRNA expression profiles in mouse spleen and brain. Most of differentially expression genes are associated with antiviral response of host, which may provide important information for investigation of JEV pathogenesis and therapeutic method.


Background
Japanese encephalitis virus (JEV), a mosquito-borne flavivirus belonging to family Flaviviridae, is responsible for an acute encephalitis and damage to the central nervous system (CNS) in wide areas of southern and eastern Asia. And recently, it has been isolated from previously non-affected areas, such as Australia [1]. Japanese encephalitis (JE) has a high fatality rate of 30% and around half of the JE survivors have severe neurological sequelae [2]. Approximately 50,000 JE cases with 10,000 deaths are reported annually [3]. Following entry into the host system through a mosquito bite, JEV may replicates in various organs such as liver and spleen, and then reaches the central nervous system, resulting in a rapid inflammatory response [4]. According to the observations from studies of other flaviviruses, specifically dengue virus, it has been proposed that JEV traverse through a lymphatic route that also involves cells of the monocyte/macrophage lineage. Recently, JEV has been shown to effectively replicate within lymphocytes and macrophages, thereby making these cell types possible carriers of the virus from the periphery to the CNS [5,6]. However, it remains to be elucidated how JEV infects the CNS via these peripheral cells. In addition, although neurological disorders caused by JEV are often characterized by evidence of immune system recognition and the presence of inflammatory components among the neuropathological changes, the mechanisms by which this virus causes neurological disease are not fully understood [7].
Recently, multiple DNA microarray analyses have been employed to detect changes in host gene expression after viral infection, which makes it possible to reveal molecular pathways that govern viral pathogenesis. Genechip analysis of human umbilical vein endothelial cells infected with Dengue Virus (DV) detected the upregulation of 269 genes and downregulation of 126 genes [8]. Gene profiling study of West Nile Virus (WNV) infected human embryonic kidney cells, human glioma cells and mice tissues were also performed [9,10]. Furthermore, identification of gene profiles in JEV-infected neuroblastoma cells and brain tissue have been reported recently, suggesting an increased expression of pro-inflammatory cytokines, chemokines, and anti-viral response genes after JEV infection [11,12]. However, both of studies on JEV were restricted to CNS, and few gene profiling researches about response in peripheral immune system has been carried out.
In present study, to globally identify candidate host genes associated with JEV pathogenesis, DNA microarray technology was utilized to investigate mRNA profile in spleen and brain tissues of mice infected with JEV wild strain P3, and some of the selected genes were further confirmed by quantitative RT-PCR. It was demonstrated that JEV infection resulted in significant changes in the expression of numerous genes in spleen and brain tissues, which could be crucial messages for revealing of JEV pathogenesis.

mRNA expression profile of JEV-infected mice
A mouse whole gene array was used to perform a systematic analysis of mRNA expression profile of spleen and brain tissues of JEV P3-infected mice which were sacrificed at day 3 and day 6 post-inoculation respectively. Genes that had ≥ |2|-fold change were identified as significantly differential expression. Of 41174 genes represented on the chip, 437 genes were differentially expressed in mouse spleens and 1119 genes were differentially regulated in brains in response to JEV infection (change fold ≥ 2.0, p value < 0.05). Unsupervised clustering ( Figure 1) analysis of the expression profiles showed a distinct mRNA signature in both spleens and brains during JEV infection. To elucidate the correlation between gene expression pattern and JEV infectioninduced biological processes, functional classification of mRNA transcripts and pathway analysis were performed. Differentially regulated genes in spleens of JEVinfected mice are involved in the biological processes such as cellular process, biological regulation and immune system process, etc ( Figure 2A). And the significant pathways of differentially expressed genes are known to be involved in cytokine-cytokine receptor interaction, natural killer cell mediated cytotoxicity, antigen processing and presentation, and chemokine signaling, etc (Table 1). While biological processes which showed differentially regulated genes in brains of JEV infected mice are cellular process, metabolic process, and biological regulation, etc ( Figure 2B). And the significant pathways of differentially expressed genes in mouse brain are cytokine-cytokine receptor interaction, MAPK signaling, neuroactive ligand-receptor interaction, and toll-like receptor signaling, etc. (Table 2).

Genes with differential expression in spleens of JEVinfected mice
In spleens of mice with JEV infection, 263 genes were detected to be significantly upregulated and 174 genes were downregulated (Table 3). Genes with increased expression in spleens of JEV infected mice mainly fell into the function of immune response to viral infections. These include pro-inflammatory cytokine IFN-γ, IFN response transcription factor IRF7, IFN-induced proteins like IFIT1, IFITM3 and IFITM7, protein degradation gene ubiquitin-specific protease Usp29 and Usp18, apoptosis related genes granzymA, granzymeB, Porferin, and IMNB2, killer cell lectin-like receptors KLRC1, KLRC2 and KLRC3, and chemokines such as CXCL10, CXCL11, CCL12, CCL2 and CCL9. The marked increase in expression of these genes implies the occurrence of a strong antiviral protective response to JEV infection. The significantly downregulated genes are mainly involved in cell adhesion molecules such as monocyte/macrophage-lineage cell marker CD163, transmembrane cell surface receptor of Langerhans cells CD207, and ligand for myeloid cells receptor CD200. Evidence was also found for decreased expression of interferon transcription factor IRF6 and interleukin 7 receptor, which may contribute to JEV pathogenesis as well.

Genes with differential expression in brains of JEVinfected mice
Compared to the expression after mock infection, 551 genes were detected to be significantly upregulated and 568 genes were downregulated in brains of JEV-infected mice (Table 4). Consistent to the results of spleen, chemokines CCL2, CCL3, CCL4 and CXCL10, IFN response transcription factor IRF7, IFN-induced proteins like Ifit1, Ifit2 and Ifit3, protein degradation gene Usp18 were obviously upregulated, and CD163 and IRF6 showed decreased expression upon JEV infection. In addition, increased expression of IFN-inducible transcription factors STAT1 and STAT2, TNF-α induced protein TNFAIP3, apoptosis-related proteins caspase3 and caspase4, suppressors of cytokine signaling Socs1 and Socs3, pro-inflammatory cytokines IL-1 and IL-6, TLR7, as well as IFN response antiviral genes of OAS family were also observed in microarray analysis. These results suggested the occurrence of a strong inflammatory response in mouse brain.

Confirmation of microarray data by RT-qPCR
To confirm the microarray hybridization results, quantitative RT-PCR was performed on eight selected differentially expressed mRNAs in mouse spleen and brain respectively. As shown in the RT-qPCR result of spleen, granzymA, granzymeB, Porferin, IRF7, IFN-γ, CXCL10, and ILIR2 were significantly upregulated, while CD163 was downregulated ( Table 5). Out of eight tested mRNAs in brains of mice infected with JEV, CCL2, CCL4, CXCL10, Casp3, Casp4, SOCS1 and SOCS3 showed increased expression, and CD163 was also found to have an obviously decreased expression (Table 5).
Although absolute values are not identical due to the different sensitivity between the techniques, all genes showed a well comparative expression pattern with microarray data.

Discussion
Spleen is one of the major peripheral immunity organs and CNS is the ultimate infection target of JEV.
Therefore, identification of the JEV-related host genes in spleens and brains of mice infected with JEV may be helpful for understanding of JEV pathogenesis. To this end, profiles of mRNA expression in both spleen and brain tissues of JEV-infected mice were analyzed systematically in this study.
In mRNA profiling assay, we found chemokines like CCL2 and CXCL10 were significantly up-regulated in both mouse spleen and brain in response to JEV-infection, suggesting a strong inflammatory response of host. Monocyte chemoattractant protein-1 (MCP-1/CCL2) is one of the key chemokines that regulate migration and infiltration of monocytes/macrophages. It was involved in neurological disorders such as encephalitis-related neuronal death, where its levels were elevated in astrocytes leading to neuronal death [13]. Previous study has demonstrated that the WNV-infection stimulated the expression of CCL2 in mice livers, suggesting a consistent result to our study of JEV [10]. CXCL10 is also found to play a crucial role in the host defense response against various viral infection of the CNS by enhancing innate immune responses [14,15], and our result of upregulated CXCL10 mRNA has a good agreement with that was reported by Gupta et al. and Biswas et al. [11,12].
In addition, a strong IFN-pathway-related response was evident in mouse spleen and brain infected with JEV, with increased expression of IFN-γ, IFN response transcription factor IRF7, IFN-induced proteins IFIT1, IFITM3 and IFITM7 in spleen, and IRF7, IFN-inducible transcription factors STAT1 and STAT2, IFN-induced proteins IFIT1, IFIT2 and IFIT3 in brain, implying the occurrence of a protective response of host. Similar results were shown in previous reports as the upregulation of IFN-γ, STAT1 and STAT2 in JEV-infected mice brain [12], and increased expression of STAT1 and IFIT3 in JEV-infected N2A cells [11]. The reason why no increase in expression of IFN-γ was found in brains of JEV-infected mice in our study may be related to the difference of time points with Biswas' study [11]. IRF7 is activated in the presence of double stranded RNA following virus infection, which is functional as one of the regulators of the IFN-α/β gene promoter and the IFNα/β responsive genes to create an antiviral state [16]. The increased expression of IRF7 has also been demonstrated in WNV infected mice, but wasn't found in JEV infected cells [17,18]. This may be due to the different signal pathways between intact host and cell culture.
Pore-forming protein perforin and the family of granzymes have been demonstrated to form an antiviral arsenal central to the function of cytotoxic T lymphocyte (CTL) and natural killer (NK) cells [19][20][21]. After JEV infection, pro-apoptotic genes found to have significant enhancement in mouse spleen, including porferin,   granzyme A and granzymeB, which suggested a strong cytotoxic response against JEV infection. The consistent result as increased expression of granzyme A and B was also shown in the report about WNV infection [17]. Granzyme B is the most characterized granyzme which plays an important role in inducing apoptosis, and it is generally accepted that granzyme A can trigger a distinct nonapoptotic form of cell death [22]. The high expression of granzyme A and B in spleen could help clear virus infection, but may also involve lymphocyte injury. The upregulation of apoptosis-related proteins Caspase3 and 4 were also detected in brains of mice infected with JEV, indicating an inflammation-related neuronal apoptosis. Caspase 3 is an effector caspase that function as a central regulator of apoptosis. It has been reported that JEV infection triggers apoptosis in different cells, such as baby hamster kidney BHK-21 cells, mouse neuroblastoma N18 cells, human neuronal NT-2 cells, and human medulloblastoma cells, resulting in caspase 3 activation [23][24][25]. The function of caspase 4 is not fully known, but it is believed to be an inflammatory caspase, with a role in the immune system [26]. Pro-inflammatory cytokines like IL-1 and IL-6 were found to have increased expression in the brain, which was consistent with the results of the studies on DV and WNV [8,17]. IL-1α and IL-1β has been known to form an important part of the inflammatory response of the body against infection. IL-6 together with IL-1 and TNF-α acts as an endogenous pyrogen by causing fever following viral infections [27]. It's also associated with an unfavorable outcome following yellow fever virus infection [28]. Evidence suggests that circulating IL-6 can activate CNS mechanisms resulting in the development of the febrile response during disease [29]. Upregulation of IL-1 and IL-6 in brain may thus be protective against harmful JEV infection, but also have a pathogenic role in CNS. The 2', 5'-oligoadenylate synthetase (OAS) and its downstream effector RNase L play important roles in host defense against virus infection. The OAS1b protein has been described as a flavivirus resistance factor, and OASl1 as a WNV-resistance factor in wild mice because a truncated version of the protein is expressed in laboratory mice which are susceptible to infection [30][31][32]. Human OAS1 p42/p46 and OAS3 p100 are likely to contribute to host defense against DEN infection and play a role in determining the outcomes of DEN disease severity [33]. Further, the activated expression of OAS2 has been demonstrated in mouse brain in response to JEV infection [34]. In present study, there was increased  expression of various members of the OAS family in the brains of JEV-infected mice, including OAS1a, OAS1b, OAS1c, OAS1e, and OAS2. Therefore, the ability of these proteins to protect against JEV infection should also be further studied. The significant downregulation of CD163 was detected both in mouse spleen and brain in response to JEV infection. CD163 is a novo identified marker for perivascular macrophages in humans, monkeys, and mice. And previously studies have found that perivascular CD163 expression is upregulated and the number of CD163-positive cells increases in HIV and SIV encephalitis (HIVE and SIVE) brains [35,36]. However, CD163 is not a "classical" activation marker, because peripheral blood monocytes and most tissues macrophages of normal uninfected controls all express it and because in vitro pro-inflammatory stimuli largely down-regulate its expression. These findings suggest that CD163 expression is regulated in

Conclusions
In summary, our findings suggested that JEV infection resulted in significant changes in the expression of multiple genes in mouse spleen and brain, including inflammatory cytokines, chemokines, IFN inducible genes, IFN regulators, and apoptosis related genes, etc. These genes may play a critical role on antiviral response of host against JEV infection but could also contribute to the pathogenesis of JEV resulting in encephalitis. The mRNA profile obtained by microarray analysis in this study may provide a foundation for future investigation of JEV pathogenesis and therapeutic method.

Virus infection of mice
Four-weeks-old naive female BALB/c mice were purchased from Wuhan Institute of Biologic Products (Hubei Province, China) and inoculated subcutaneously (s.c.) with 5 × 10 6 PFU of JEV P3 strain. Dilutions were performed in serum-free DMED and experimental controls were mock injected with diluent. A part of JEV-infected mice were sacrificed at day 3 post-inoculation (the day before the neurological symptoms appeared), and spleens were harvested. The left mice were sacrificed at day 6 post-inoculation (the day before the mice started to die), and brains of mice were harvested. Spleen and brain homogenate was made in DMEM for RNA extraction.

Microarrays and bioinformatics
The total RNA was isolated from mouse spleens and brains respectively with trizol reagent (Invitrogen) for mRNA Microarray. mRNA hybridization was performed by shanghaiBio Corporation (shanghai, China) with the use of 4 × 44 K Agilent Whole Mouse Genome Oligo Microarray (total 41,174 oligo probes from 41,174 mouse genes). For each sample pair, the experiments were done with two independent hybridizations (Cy3 and Cy5 interchanging labeling). Hybridized arrays were scanned at 5 μm resolution on an Agilent DNA Microarray Scanner (Model G2565BA). Data extraction from images was done by using Agilent Feature Extraction software. Hierarchical cluster, gene ontology and pathway analysis were analyzed by using SAS (ShanghaiBio Analysis System).
Quantitative real-time RT-PCR (RT-qPCR) For selected mRNA RT-qPCR, total RNA from the same samples used in microarray analysis was tested by using ABI 7500 FAST real-time PCR system. PCR primers were designed with Primer Express 2.0 software (Invitrogen). Results are shown as fold change. For mRNA RT-qPCR, experiments were carried out with the Primer-Script RT reagent Kit (TaKaRa) and SYBR Green Realtim PCR Master Mix (TaKaRa) according to manufacture's instruction. The housekeeping gene GHDAP was used for standardization of the initial RNA content of a sample. Relative changes of gene expression were calculated by the following formula, and the data are represented as fold upregulation/downregulation. fold change = 2 -ΔΔCt , whereΔΔCt =(Ct of gene of interest, treated -Ct of HK gene, treated)-(Ct of gene of interest, control-Ct of HK gene, control), Ct is the threshold cycle number and HK is the house keeping gene.

Statistical analysis
Each gene in each infection group was subjected to a Student's t test to detect large expression differences relative to the mock-infected group. p-values <0.05 were considered to be statistically significant.