Identification of host genes leading to West Nile virus encephalitis in mice brain using RNA-seq analysis

Differential host responses may be critical determinants of distinct pathologies of West Nile virus (WNV) NY99 (pathogenic) and WNV Eg101 (non-pathogenic) strains. We employed RNA-seq technology to analyze global differential gene expression in WNV-infected mice brain and to identify the host cellular factors leading to lethal encephalitis. We identified 1,400 and 278 transcripts, which were differentially expressed after WNV NY99 and WNV Eg101 infections, respectively, and 147 genes were common to infection with both the viruses. Genes that were up-regulated in infection with both the viruses were mainly associated with interferon signaling. Genes associated with inflammation and cell death/apoptosis were only expressed after WNV NY99 infection. We demonstrate that differences in the activation of key pattern recognition receptors resulted in the induction of unique innate immune profiles, which corresponded with the induction of interferon and inflammatory responses. Pathway analysis of differentially expressed genes indicated that after WNV NY99 infection, TREM-1 mediated activation of toll-like receptors leads to the high inflammatory response. In conclusion, we have identified both common and specific responses to WNV NY99 and WNV Eg101 infections as well as genes linked to potential resistance to infection that may be targets for therapeutics.

Scientific RepoRts | 6:26350 | DOI: 10.1038/srep26350 WNV. Very recently, high-throughput RNA sequencing (RNA-seq) technology, which is a powerful way to profile the transcriptome with great efficiency and higher accuracy, has been employed in various viral infections and diseases 8,9 . RNA-seq technology is highly sensitive, quantitative, and provides unbiased measurement of gene transcripts in a large dynamic range of gene abundance. In this study, we used RNA-seq technology to analyze host transcriptomes in the brains of mice infected with WNV NY99 and WNV Eg101 to gain insight into the underlying mechanisms of the differential pathogenicity of these viruses. We have characterized and compared the host responses to the two WNV strains in terms of canonical pathways and biological functions, as well as in terms of individual genes.

Results and Discussions
Distinct and dynamic changes in the expression of cellular genes in WNV NY99-and WNV Eg101-infected mice brain. To identify the dynamics of host gene expression associated with WNVE, RNA-seq analysis was conducted to explore the transcriptomes from mice brains infected with lethal, WNV NY99, and non-lethal, WNV Eg101 viruses, as well as mock-infected mice. After peripheral inoculation, WNV is first detected in the brain between days 4 and 6 after footpad inoculation and peak virus load is observed at day 8 after inoculation 10,11 . WNV NY99-infected mice demonstrate various neurological symptoms such as paresis, tremors, and ataxic gait at day 8 after inoculation. Therefore, to understand the transcriptome changes occurring at the time of high WNV replication in brain and when animals display clinical symptoms of encephalitis, we examined gene expression at day 8 after infection. Virus was not detected in the mock-infected brains. WNV was detected in the brains of both WNV NY99-and WNV Eg101-infected mice by plaque assay. However, virus titers were significantly higher in the brains of WNV NY99-infected mice when compared to WNV Eg101infected mice (Fig. 1A). RNA harvested from the brains of mock-, WNV NY99-and WNV Eg101-infected mice (four mice per group) were subjected to high-throughput sequencing on an Illumina HiSeq 2500 sequencer as described previously 8,9 .
Gene expression profiles for WNV NY99-and WNV Eg101-infected brains were compared with mock-infected brains. Figure 1B demonstrates total numbers of up-and down-regulated differentially expressed genes for WNV NY99-and WNV Eg101-infected mice. Number of differentially expressed genes was higher in the brains of WNV NY99-infected mice, which correlate with significantly higher brain viral load observed in these mice compared to WNV Eg101-infected mice (Fig. 1A). In WNV NY99-infected mice, 1,036 genes were up-regulated with log 2 fold change values ranging from 0.48 to 9.6 and 364 genes were down-regulated with log 2 fold change values ranging from − 0.48 to − 9.6. In WNV Eg101-infected mice, 184 genes were up-regulated with log 2 fold change values ranging from 0.48 to 7.4 and 94 genes were down-regulated with log 2 fold change values ranging from − 0.48 to − 2.6 ( Fig. 1B). The number of up-regulated genes was larger than the number of down-regulated genes in both groups of virus-infected mice. Genes involved in biological pathways related to innate antiviral immunity and inflammation responses were among the top differentially expressed genes after WNV NY99 infection. These include well-characterized responses such as type I interferon (IFN) signaling and chemokine signaling. Similar to WNV NY99, genes associated with IFN signaling were among the top differentially expressed genes after WNV Eg101 infection. Tables 1 and 2 depict the top 10 up-regulated and down-regulated genes in WNV NY99-and WNV Eg101-infected mice brains, respectively. The complete list of differentially expressed genes can be found in Supplementary tables S1 and S2.
Furthermore, Venn diagram was generated to examine the unique and overlapping genes for both groups. As depicted in overlapping circle in Fig. 1C, 147 genes were common to infection with WNV NY99 and WNV Eg101. Nighty two up-regulated genes and 48 down-regulated genes were common to infection with both viruses, indicating that 50% (92/184) of genes that are up-regulated in the brain following infection with WNV Eg101 are also up-regulated following infection with WNV NY99 and that 51% (48/94) of genes that are down-regulated in the brain following infection with WNV Eg101 are also down-regulated following infection with WNV NY99. Interestingly, seven genes common to infection with both the viruses demonstrated inverse trend ( Table 3). Expression of HOXB5, GFAP, HBA-A1, CYR61 and α 2M genes was up-regulated in WNV NY99 infection and down-regulated in WNV Eg101 infection. In contrast, expression of SYNPO2 and PTPN3 genes was down-regulated in WNV NY99 infection and up-regulated in WNV Eg101 infection. GFAP and α 2M genes are particularly of interest. GFAP is an astrocyte marker and its up-regulation has been demonstrated in WNV infection. Others and we have demonstrated that activation of astrocytes is associated with enhanced production of pro-inflammatory mediators after WNV infection 12,13 . α 2M is a protease inhibitor and cytokine transporter with important role in inflammation. α 2M is synthesized in the brain primarily by astrocytes, and expression of its receptor has been identified in neurons and astrocytes 14 . Several studies demonstrated that α 2M is a specific cytokine carrier, which binds major pro-and anti-inflammatory cytokines, such as interleukin (IL)-1β , IL-6, TNF-α , and others 15 . It has been demonstrated that α 2M −/− mice are resistant to the lethal effects of lipopolysaccharides 16 . α 2M is also implicated in Alzheimer disease (AD) due to its ability to mediate the clearance and degradation of α -beta, the major component of beta-amyloid deposits 17 . However, role of α 2M is yet to be determined in WNV infection and other inflammatory infectious diseases. Up-regulation of both GFAP and α 2M genes in brain cells especially astrocytes may contribute to enhanced inflammation observed in the brains of WNV NY99-infected mice.
Genes that were up-regulated in infection with both the viruses are mainly associated with IFN signaling. Interferon stimulated genes (ISG) such as MX1, IFIT1, IFIT2, IFIT3, and RSAD2 (viperin) were significantly up-regulated after infection with both the viruses ( Table 4). The IFN response is central to the innate defense mechanisms of the host against WNV infection 18 . The paracrine and autocrine secretion of IFN renders cells "antiviral" by inducing several genes including ISG, the dsRNA-dependent protein kinase R (PKR) and 2′ -5′ -oligoadenylate synthase (OAS)/RNase L. These ISG confer an antiviral state by blocking virus replication at different levels such as early-stage virus infection, inhibition of post-transcriptional modification and virus maturation 18 . It has been demonstrated that viperin can inhibit WNV infection in cells and in vivo 19 . Members of the IFN-inducible tetratricopeptide (IFIT) family of genes (IFIT1, IFIT2, and IFIT3) also contribute to the control of WNV infection 18 . IFIT2 is important for controlling virus replication in the brain 20 . The absence of PKR signaling or a deficient OAS/RNaseL pathway also leads to increased susceptibility to WNV infection 21,22 . Genetic variation in OAS1 also has been suggested as a risk factor for initial infection with WNV in humans 23 .
Interestingly, genes associated with inflammation such as cytokines, chemokines and their receptors were exclusively expressed only after WNV NY99 infection. WNV NY99 infection induced a strong up-regulation of multiple cytokines and chemokines in mice brains (Table 4). Chemokine pathway-associated genes, such as C-X-C motif chemokine 10 (CXCL10), chemokine (C-C motif) ligand 5 (CCL5), and CCL2, were among the 10 most up-regulated genes after WNV NY99 infection. CXCL10 was the topmost up-regulated gene in WNV NY99-infected mice. Most of the chemokine, cytokines and their receptors were significantly up-regulated after WNV NY99 infection, including CXCL1, CXCL13, CXCL16, CCL7, CCL4, CCL3, CCL8, CCL19, IL-1β , IL-1α , and IL-7. In contrast, WNV Eg101 infection did not induce the expression of these pro-inflammatory mediators ( Table 4). Only CXCL12 gene was up-regulated in WNV Eg101 infection. Similar to IFN response, WNV-induced pro-inflammatory mediators are also known to protect mice from lethal WNV disease. CXCL10 promotes trafficking of WNV-specific CD8 + T cells via binding to its cognate receptor CXCR3 24 . Enhanced expression of CCL3, CCL4, and CCL5 by WNV infection leads to CCR5-dependent trafficking of CD4 + and CD8 + T cells, NK cells, and macrophages into the brain 25 . CCL2 is important in the trafficking of inflammatory monocyte into the brain after WNV infection 26 . IL-1β promotes immune cell trafficking into the brain 27 . Interestingly, we observed an enhanced production of key pro-inflammatory cytokines and chemokines in only WNV NY99-infected mice  Table 3. Commonly expressed genes demonstrating inverse trend in WNV NY99 and WNV Eg101 infections.
despite increased disease severity observed in these mice. This could be due to significantly increased antigen load in WNV NY99-infected mice compared to WNV Eg101-infected mice (Fig. 1A). Collectively, we demonstrate both predicted and novel gene changes in response to infection with lethal, WNV NY99 and non-lethal, WNV Eg101 viruses. Previous studies have not reported the entire transcriptome profile of WNV-infected mice brain using RNA-seq. Qian and colleagues reported RNA-seq analysis of WNV-infected human macrophages in vitro 28 . Similar to our data, this study also identified cytokine/chemokine signaling and IFN signaling as key signaling pathways modulated by WNV infection in human macrophages. Further, Qian and colleagues demonstrated the important roles of IFI27, AIM2, CCR3 (receptor of chemokines CCL5 and CCL8) and CXCR3 (receptor of chemokines CXCL10) in resistance to infection with WNV. Our data also demonstrates up-regulation of CCL5, CCL8, and CXCL10 in WNV NY99-infected mice (Table 4). However, we did not observe up-regulation of IFI27 and AIM2. This could be due to the differential antiviral response to WNV infection in peripheral macrophages compared to resident brain cells.
The major set of genes/pathways with increased expression observed in the present study with WNV are similar to those reported for other neuroinvasive viruses like Japanese encephalitis virus (JEV). Clarke and colleagues conducted microarray analysis of WNV NY99-and JEV-infected mice brains and reported that flaviviruses-induced transcriptional changes in the brain include the differential expression of genes associated with inflammation, apoptosis, interleukin 17 receptor A, IFN signaling, glutamate signaling, and tRNA charging pathway 29 . Similarly, our data also demonstrate induction of aforementioned genes in the brain after WNV NY99 infection; however, we did not observe the differential expression of genes associated with tRNA charging pathway. This could be due to the difference in the mice strain (C57BL/6J vs. Swiss Webster) and route of inoculation (subcutaneous vs. intracranial) used in the study by Clarke and colleagues. We inoculated mice by footpad route, which mimics the natural route of WNV infection. Other studies have also reported large-scale gene expression changes in the brain following subcutaneous JEV infection. Similar to our study, genes that are up-regulated in subcutaneous JEV infection are mainly associated with IFN signaling, pattern recognition receptors, chemotactic genes, activation of inflammasome, and cell death/survival signaling 30,31 . These results illustrate the pathways that are commonly activated during neurotropic viral infections of the central nervous system (CNS).
Confirmation of RNA-seq data using qRT-PCR. qRT-PCR was conducted on additional virus-infected (n = 7 in each group) and mock-infected (n = 4) animals to confirm the gene expression changes of a selected number of differentially expressed genes. Similar to RNA-seq data, MX1, IFIT1, and IFIT2 genes were significantly up-regulated in the brains of both WNV NY99 and WNV Eg101-infected mice compared to mock-infected mice (Fig. 2). The up-regulation of inflammatory mediators following WNV NY99 infection, but not following WNV Eg101 infection, was also confirmed by qRT-PCR. Pro-inflammatory mediators, IL-1β and CCL3, expression was only up-regulated in WNV NY99 infection. Others and we have previously demonstrated dramatic increase in the protein levels of multiple cytokines and chemokines in WNV NY99-infected mice brains at day 8 after inoculation, but this increase was not observed in the brains of WNV Eg101-infected mice 7 . Moreover, expression of GFAP and α 2M genes was up-regulated in WNV NY99 infection and down-regulated after WNV Eg101 infection (Fig. 2).

Functional analysis of WNV-modulated genes.
To investigate the biological interactions of differentially expressed genes and to identify important functional networks, significantly modulated differentially expressed genes derived from RNA-seq analyses were imported into the Ingenuity Pathway Analysis (IPA) tool as described previously 8,9,[32][33][34] . The highest activated networks (high z-score) were identified using IPA as described in the materials and methods. Figure 3A,B depict the top ten activated canonical pathways after WNV NY99 and  (Fig. 3A). Moreover, pathological pathways such as 'Death Receptor Signaling' and 'Retinoic acid-Mediated Apoptosis Signaling' were also activated after WNV NY99 infection. Two of the top ten activated canonical pathways after WNV NY99 infection were shared with those affected by WNV Eg101 infection; 'IFN Signaling' and the 'Role of PRR in Recognition of Bacteria and Viruses' (Fig. 3B). Other highly activated pathways related to WNV Eg101 infection included 'Growth Hormone Signaling' , 'Cardiac β -adrenergic Signaling' , 'Calcium Signaling' , and ' Activation of IRF by Cytosolic Pattern Recognition Receptors' . We next generated heat map to compare the activated canonical pathways after WNV NY99 and WNV Eg101 infections. As depicted in Fig. 4, most of the canonical pathways were activated by infection with both the strains, although to different extents. Activation was more enhanced in WNV NY99-infected mice than in WNV Eg101-infected mice, which is consistent with higher brain viral titers observed in WNV NY99-infected mice. The three important canonical pathways activated exclusively in WNV NY99 infection were the 'TREM-1 Signaling' , 'Death Receptor Signaling' , and 'TLR Signaling' . Pathways related to cell death were only activated after WNV NY99 infection. To further understand the role of these differentially activated canonical pathways after WNV NY99 and WNV Eg101 infection in disease pathogenesis, we generated the network maps of 'PRR in Recognition of Viruses and Bacteria' , 'TREM-1 Signaling' , and 'Death Receptor Signaling' .
Pattern recognition receptors. The innate immune system encodes a series of PRR that, upon recognition of a viral pathogen-associated molecular pattern (PAMP), induce a potent antiviral host immune response 35 . The activation of these PRR is very tightly regulated at multiple steps to prevent tissue damage due to uncontrolled production of cytokines. The major PRR are the retinoic-acid inducible gene-I (RIG-I)-like receptors (RLR), TLR, and the nucleotide oligomerization domain (Nod)-like receptors (NLR). Members of the RLR (RIG-I and MDA5) and TLR (TLR3 and TLR7) families are the dominant PRR that detect WNV infection 36,37 . These receptors drive complementary antiviral defense programs by activation of downstream signaling resulting in the production of antiviral IFN, ISG and pro-inflammatory cytokines 36,37 . Our data demonstrate that WNV NY99 infection induces the expression of all three major PRR; RLR, NLR and TLR (Fig. 5A). Interestingly only RLR were induced after WNV Eg101 infection (Fig. 5B). The RLR are a family of cytosolic RNA helicase proteins comprised of three members: RIG-I, myeloma differentiation antigen 5 (MDA5), and LGP2 36 . The importance of the RLR signaling pathway in protection against WNV has been demonstrated by several studies in vivo. RLRs detect and respond to WNV by inducing a potent antiviral defense program, characterized by production of type I IFN and expression of antiviral effector genes including ISG. Both RIG-I and MDA5 are required for optimal antiviral responses against WNV 38,39 . RIG-I or MDA5 knockout mice show increased mortality compared to wild type mice after WNV infection 40,41 . Toll-like receptors. While RLR pathway is known to be protective after WNV infection, the role of TLR pathway has been controversial. The TLR family is composed of more than 10 members, with each acting as a sensor of conserved microbial component, that drive the induction of immune responses 42 . Recognition of WNV by TLR occurs in endosomes and is largely mediated by TLR7 and TLR3, which bind ssRNA and dsRNA, respectively. The role of TLR3 in WNV infection remains controversial, as two independent studies using the qRT-PCR was conducted on RNA extracted from mock and WNV NY99-and WNV Eg101-infected mice brains at day 8 after infection to determine fold-change in MX1, IFIT1, IFIT3, CCL3, IL-1β , GFAP, and α 2M gene expression. Changes in the levels of each gene were first normalized to the β -actin gene and then the fold-change in WNV-infected mice brain was calculated in comparison to corresponding mock-infected mice brain. Data represents the mean ± SEM, representing two independent experiments (n = 7 per group). *P < 0.05. same mouse model demonstrated opposing phenotypes of enhanced susceptibility and resistance to WNV infection 43,44 . Similarly, the function of TLR7 in WNV infection is less clear, as one study found no difference in susceptibility to WNV between wild-type (WT) and TLR7 −/− mice after intradermal WNV infection 45 , and in other study TLR7 −/− mice were more susceptible to intraperitoneal WNV infection than WT mice 46 . Protection in TLR3 knockout mice is attributed to the reduced production of pro-inflammatory mediators from brain cells, which in turn reduce primary neuronal damage and promote blood-brain barrier stability 43 . Knockout of both TLR3 and TLR7 have little systemic impact on IFN production after WNV infection. Our data demonstrate that WNV NY99 infection induces the expression of seven members of TLR family including TLR1, TLR2, TLR3, TLR4, TLR7, TLR8, and TLR9 ( Fig. 5A and Table 5). In contrast, WNV Eg101 infection does not induce expression of any TLR (Fig. 5B and Table 5). Furthermore, we observed increased levels of several pro-inflammatory mediators and their cognate receptors only in the brains of WNV NY99-infected mice, which correlate with TLR activation in these mice.
Although TLR are important for viral recognition, excessive TLR activation may result in a pathogenic rather than a protective outcome in infections of the CNS. In this regard, antagonizing TLR signaling has been suggested as a therapeutic approach to control brain inflammation and overall disease pathology 43,47,48 . Other than WNVE, activation of TLR pathway is also demonstrated to be associated with encephalitis in the SIV model of neuroAIDS 49 . Recent study demonstrated that macrophages from young individuals can down-regulate TLR3 following infection with WNV, whereas macrophages of the elderly cannot, resulting in production of high levels of pro-inflammatory and vasculogenic cytokines in this group 50 . This is an interesting observation as elderly individuals are far more susceptible to develop WNVE than young individuals 51 . Moreover, TLR4 signaling during infection with highly pathogenic H5N1 influenza A virus has been reported to contribute towards lung pathology 52 . Similar to our study, it has been recently demonstrated that peripheral inoculation of WNV Eg101 in WT mice at doses lower than 10 6 PFU did not produce any clinical sign of disease or mortality, and no TLR activation 53 . However, inoculation of WNV Eg101 at a higher dose (10 7 PFU) induced significant mortality and TLR activation, suggesting that activation of TLR pathway is a result of high virus replication 53 . Similarly, in this study, high virus load in WNV NY99-infected mice results in TLR activation and subsequently enhanced inflammation observed in these mice.
Nod-like receptors. NLR are soluble or cytosolic receptors in the mammalian cell cytoplasm. Our data demonstrate activation of NOD1, NOD2, NLRP3, NLRP6, and NLRC5 following WNV NY99 infection (Figs 5A and 6A). In contrast, only NLRP6 was induced after WNV Eg101 infection (Figs 5B and 6B). NLRC5, and NLRP3 can form inflammasomes by recruiting apoptosis-associated speck-like protein containing a CARD (ASC) to activate caspase-1, leading to the cleavage of pro-IL-1β to mature IL-1β . Activation of the NLRP3 inflammasome and subsequent secretion of IL-1β family pro-inflammatory cytokines have been implicated in the protection against WNV infection 54 . NLRP3 inflammasome has also been demonstrated to play role in the pathogenesis of Japanese encephalitis 55 . We have previously demonstrated that ASC is critical for the immune response and survival in WNVE 56 . However, role of NLRC5 and NLRP6 is yet to be determined in WNV pathogenesis. In addition to RLR, TLR, and NLR, PKR and OAS are classes of IFN-inducible PRR that can recognize dsRNA and restrict a number of viruses, including WNV. Our data demonstrate activation of PKR and OAS after infection with both the viruses (Fig. 5A,B).
Collectively, this analysis demonstrated that differences in the activation of key PRR resulted in the induction of unique innate immune profiles, which corresponded with the induction of IFN and inflammatory responses.  TREM-1 signaling. TREM1 is known to activate neutrophils and monocytes/macrophages by signaling through the 12-kDa-adapter protein DNAX activation protein (DAP12), and it activates TLR and NLR signaling pathways via secretion of pro-inflammatory cytokines and chemokines in response to microbial infection 57 . TREM-1 signaling has been indicated to play role in the pathogenesis of several inflammatory disorders of both infectious and non-infectious etiology. However, the significance of TREM signaling in innate immunity to virus infections and the underlying mechanisms remain largely unclear. Our data demonstrates that activation of TREM1 signaling is specific to pathogenic WNV NY99 infection. Majority of the genes involved in TREM-1 pathway including TLR, NLR, MCP-1, JAK2, IL-1β , cell adhesion and co-stimulator molecules, were up-regulated after WNV NY99 infection (Fig. 6A). In contrast, only MCP-1 and NLRP6 were up-regulated after WNV Eg101 infection (Fig. 6B). This finding could partially explain the high inflammatory response observed following WNV NY99 infection, as TREM-1 is directly implicated in the enhancement of inflammation by activating TLR and NLR signaling pathways. A substantial amount of data supports the link between TREM-1 signaling and several diseases, such as polymicrobial septic shock and inflammatory bowel disease, and in animal models of pneumonia and asthma 57 . TREM-1 signaling contributes to the cytokine storm associated with lethal filovirus disease 58,59 . It has been recently demonstrated that TREM-1 deficient mice were protected from severe influenza disease 60 . Another recent study by Suthar et al. used transcriptional profiling and pathway modeling to show that TREM-1 signaling was enriched in the liver following infection with WNV in mice 61 . Although this study implies that TREM-1 signaling might be one of the pathways responsible for restricting tissue tropism, the precise role of TREM-1 in WNV pathogenesis has not been explored. Based on our data, it is likely that TREM-1 signaling may promote inflammation and cause substantial tissue damage thereby contributing to disease pathogenesis of WNV NY99 infection.
Death receptor signaling. Many viruses encode proteins that can inhibit apoptosis and the triggering of apoptosis has been noted in WNV NY99-infected neuronal tissue 62 . We next generated a network map of death receptor signaling for WNV NY99-and WNV Eg101-infected mice. FasL, TNF-α , and Apo2L (TRAIL)-mediated downstream signaling leading to cell apoptosis were only activated in WNV NY99-infected mice, which also correlate with more severe pathogenicity observed in WNV NY99-infected mice (Fig. 7A,B). Downstream genes involved in cell death (DAXX, FLIP, TBID, PARP, CASP7, and CASP8) were activated upon WNV NY99 infection, but not following WNV Eg101 infection.
In Summary, we have provided a comprehensive transcriptome analysis of mice brains that were infected with lethal and non-lethal WNV strains. One of the limitations of this study is the lack of over-expression or knockdown studies to delineate the role of the genes involved in the observed signaling pathways. Differential expression of specific host genes observed between virulent and non-virulent WNV strains could explain the variation in disease pathogenicity. In addition, we have identified novel genes and pathways that are not previously associated with WNV infection and might provide novel therapeutic targets for the management of WNVE.

Materials and Methods
Animal experiments. Nine-week old C57BL/6J mice were purchased from The Jackson Laboratory.
Animals were housed four per cage and allowed to eat and drink ad libitum. The animal suite was maintained at 72 °F, 45% humidity and on 12/12-light/dark cycles. Sawdust bedding was provided along with paper towel. This study was approved by the University of Hawaii Institutional Animal Care and Use Committee (IACUC) (protocol number 13-1759), and conducted in strict accordance with guidelines established by the National Institutes of Health and the University of Hawaii IACUC. All animal experiments were conducted in consultation with veterinary and animal care staff at the University of Hawaii in the animal biosafety level-3 laboratory.
Mice (four animals per group) were inoculated via the footpad route with 100 plaque-forming units (PFU) of WNV NY99 or WNV Eg101 or PBS (mock). A power calculation assuming a type I error of 0.05 gave a power of 92% using four mice per group. At day 8 after infection, mice were anesthetized using isoflurane and perfused with cold PBS as described previously 10 . Brains were harvested and flash frozen in 2-methylbutane (Sigma) and stored at − 80 °C until further use.
Plaque assay. One-half of the frozen brain was weighed and homogenized in a bullet blender (Next Advance) using glass or zirconium beads, and plaque assay was conducted using Vero cells for analysis of WNV titer as described previously 10,63 . Gene expression analysis using RNA-seq. One-half of the frozen brain from four mice was powdered over dry ice to obtain a homogenous sampling. An aliquot of the frozen powder was used to isolate total RNA using RNeasy Mini Kit (Qiagen, Cat #217004) as described previously 63 . Genomic DNA contamination was eliminated by digesting the RNA with RNase-free DNase (Ambion, Austin, TX). The quality of all RNA samples was assessed using an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). RNA-seq libraries were constructed using the Illumina TruSeq RNA sample preparation v2 as per manufacturer's protocol (Illumina, San Diego, CA, USA), to generate 100-bp paired-end reads. Libraries were sequenced using an Illumina HiSeq 2500 with v3 chemistry at the University of California Irvine's Genomic High Throughput facility. A total of 1,227,480,292 reads were inspected using the FastQC (http://www.bioinformatics.babraham.ac.uk/ projects/fastqc/) program and minimally trimmed to remove low quality bases. The processed reads were then mapped to the mouse genome and the its gene annotations (mm_ref_GRCm38.p3) using TopHat (ver. 2.0.9), allowing for two mismatches 64 . The Cufflinks program (version 2.2) was used with the multi-read-correct option to estimate the normalized transcription profiles from the mapping results 65 . The Cuffdiff program from Cufflinks was used with a Benjamini-Hochberg False Discovery Rate (FDR) cut-off of 5% to assess statistically significant differential expression changes between the control and each of the infected groups.
Pathway analysis. Data were analyzed using IPA (QIAGEN Redwood City, USA) as described previously 29 .
The most significant canonical pathways and functional processes of biological importance were selected using the list of differentially expressed genes identified by RNA-seq and the IPA Knowledge Base as described previously 8,9 . Pathway enrichment p-values (Fisher's exact test) and activation z-scores were calculated by IPA. Fisher's exact test, using IPA, was used to calculate cut-off point of significance. p < 0.05 is considered significant.  Validation by qRT-PCR. qRT-PCR was conducted on additional virus-infected (n = 7 in each group) and mock-infected (n = 4) brains. Brain tissues were powdered over dry ice to obtain a homogenous sampling and an aliquot of the frozen brain powder was used to extract total RNA as described previously. The mRNA levels of multiple host genes were determined using qRT-PCR and the fold-change in infected brains compared to mock brains was calculated after normalizing to the β -actin gene as described previously 63 . The primer sequences used for qRT-PCR are listed in Table 6.
Statistical Analysis. Differentially expressed genes with an FDR adjusted p-value of the cuffdiff test statistics < 0.05 were considered significant. For plaque assay and qRT-PCR, unpaired student t-test using Graphpad was used to calculate p values as described previously 10,66 . p < 0.05 is considered significant.