Many viruses can invade the central nervous system (CNS), infect the neurons, propagate by exploiting cellular machinery, hijack axonal transport, and cross the synapses to disseminate within neural networks (Taylor and Enquist, 2015). This is also true for the West Nile virus (WNV), which can infect the CNS, use bi-directional axonal transport, and spread trans-synaptically (Maximova et al., 2016). WNV infection has affected an estimated 7 million people in the continental United States since the first outbreak in 1999 (Ronca et al., 2019). In humans, WNV infection can result in WNV neurological disease (WNV-ND) (Davis et al., 2006) with a high in-hospital mortality rate (Yeung et al., 2017), motor and cognitive impairment (Kleinschmidt-DeMasters and Beckham, 2015; Patel et al., 2015), and increased risk of death during the convalescent phase (Philpott et al., 2019). WNV-ND can present as meningoencephalomyelitis with any combination of clinico-pathological features such as aseptic meningitis, encephalitis, and acute flaccid myelitis (Sejvar, 2016). The estimated ratio of the incidence of WNV-ND relative to the number of WNV-RNA-positive blood donors is ≤1% (Betsem et al., 2017), but the current burden of WNV-ND may be underestimated due to underutilization or inaccurate choice of diagnostic tests when viral encephalitis is clinically diagnosed (Vanichanan et al., 2016). The lessons from the past two decades of extensive research using mouse models have provided a better understanding of immune control of WNV infection (Diamond and Gale, 2012). However, current knowledge is yet to be translated into effective treatment or better clinical management of WNV-ND, which remains only supportive. A major question of how virus replication in neurons and ensuing immune responses alter neural functions and result in neuropathological sequelae remains unanswered.

Postulated mechanisms of the neuropathogenic impact of WNV infection center around the loss of neurons due to either direct cytopathic effect of virus replication or damaging effects of the inflammatory environment developing in response to infection. Yet, in order to develop potential neuroprotective therapies, we need to better understand the mechanisms of neuronal dysfunction during acute infection, before neuronal loss has occurred, as well as the characteristics of neurological deficits after recovery. Loss of synapses in the hippocampus has been suggested as one of the potential mechanisms underlying cognitive impairment in patients recovering from WNV encephalitis (Vasek et al., 2016). The impact of WNV infection on neurophysiology is expected to be CNS-region-specific, which would result in the impairment of a function specific to an affected neural network. This necessitates a study of the pathogenic mechanisms of the impairment of motor functions in movement disorders, typically observed in WNV-ND (Lenka et al., 2019).

One unbiased approach to obtain a broad insight into complex pathophysiological changes occurring in the virus-infected CNS is to examine alterations in global gene expression. Although several studies have provided insights into changes in gene expression in the CNS after an experimental WNV infection, they mainly focused on the immune responses and apoptosis (Lim et al., 2017; Kosch et al., 2018; Bourgeois et al., 2011; Clarke et al., 2014a; Venter et al., 2005; Green et al., 2016; Kumar et al., 2016; Clarke et al., 2014b; Vig et al., 2018). Our understanding of the molecular mechanisms underlying changes in the regulation of the CNS during WNV-ND is incomplete. Here, we studied how the CNS is regulated when faced with a challenge to fight virus infection of neurons in a nonhuman primate (NHP) model, using combined analysis of the transcriptome, protein expression, and cell morphology.

Tissue samples for this study were selected from the CNS of NHPs inoculated intracerebrally with wild type WNV (NY99) (Maximova et al., 2014). The rationale for the selection of CNS structures for examination was as follows: (i) high viral burden at the site; (ii) severe histopathological changes at the site; and (iii) presence of neurological signs consistent with the function of the respective structure. Based on these criteria from the intracerebral NHP model (Maximova et al., 2016; Maximova et al., 2014), which are consistent with a body of data on the virology, histopathology, and clinical presentations of WNV-ND in humans (Kleinschmidt-DeMasters and Beckham, 2015; Sejvar, 2016; Lenka et al., 2019; Omalu et al., 2003; Cushing et al., 2004; Guarner et al., 2004; Armah et al., 2007; Hart et al., 2014), the cerebellum and spinal cord samples were selected for analysis. The samples were from the following predetermined time points after WNV infection: (i) 3 days post-infection (dpi; presymptomatic stage; n = 3); (ii) 7 dpi (early-symptomatic stage; moderate incoordination, tremor, and limb weakness; n = 3); and (iii) 9 dpi (advanced-symptomatic stage; severe incoordination, tremor, and limb weakness that required humane euthanasia; n = 3). The mock control cerebellum or spinal cord samples were from NHPs (n = 4) that were inoculated intracerebrally in an identical manner to the virus-inoculated animals, except that the inoculum contained only diluent and not virus (Maximova et al., 2014) (detailed procedure is described by us previously [Maximova et al., 2008]). These four mock-inoculated control animals were used for normalization of gene expression and one animal (euthanized at 10 dpi) was used as a normal control for immunohistochemistry that examined the advanced-symptomatic stage of WNV-ND (9 dpi).

WNV infection of the CNS alters the homeostasis in transcriptional regulation of the immune and neural system domains

We examined gene expression changes in the cerebellar and spinal samples from the CNS of intracerebrally infected NHPs using Affymetrix microarray analysis and validated microarray results by qPCR. A broad examination of changes in global gene expression in response to WNV infection showed transcriptional changes that progressively increased during the neurological symptomatic stages in a CNS-region-specific manner (Figure 1, a and b; top 50 upregulated or downregulated genes are listed in c and d, respectively).

Figure 1

Download asset Open asset

At the presymptomatic stage of WNV-ND (3 dpi), no genes expressed in the cerebellum were significantly differentially regulated compared to mock (false discovery rate, FDR > 0.05). In the spinal cord, however, six genes were significantly upregulated (FDR < 0.05; fold change ≥2 over mock) and none were significantly downregulated. The upregulated genes were: MX1 (12.2-fold change), IFIT3 (7.9-fold change), IFI44L (7.6-fold change), PARP9 (4.0-fold change), STAT1 (3.2-fold change), and TRIM22 (2.1-fold change). Annotation of these genes to gene ontology (GO) terms revealed a functional dichotomy for some of these genes in regulating both the immune and nervous system. While all six genes were annotated to GO biological process (BP) terms such as defense response to virus (GO:0051607) and immune system process (GO:0002376), two of these genes (MX1 and STAT1) were also annotated to the nervous system related terms such as the dendrite (GO:0030425), axon (GO:0030424), synapse (GO:0054202), and regulation of synapse structure or activity (GO: 0050803). These data suggest a dual immune-neural functionality of some early upregulated genes and offer a hint to a transcriptional activation of immune responses and concurrent modulation of neurotransmission, likely associated with sensing the spread of WNV within the CNS, which can occur trans-synaptically (Maximova et al., 2016). Since the small number of the differentially regulated genes precluded the use of statistical analyses for this presymptomatic stage, all functional genomic analyses hereafter will be focused on the early-symptomatic (7 dpi) and/or advanced-symptomatic (9 dpi) stages of WNV-ND.

At the early-symptomatic stage of WNV-ND (7 dpi), there were many more upregulated than downregulated genes in both the cerebellum and spinal cord, with much stronger transcriptional changes in the spinal cord (Figure 1e). At the advanced-symptomatic stage of WNV-ND (9 dpi), compared to early-symptomatic stage, the number of upregulated and downregulated genes increased in both the cerebellum and spinal cord (Figure 1f), with about twice as many differentially expressed genes (DEGs) in the spinal cord, compared to the cerebellum. Notably, the spinal cord also consistently displayed a higher percentage of the downregulated genes, which doubled at the advanced symptomatic stage compared to the early-symptomatic stage (Figure 1e and f). Similarity between the genes that were differentially expressed in the WNV-infected cerebellum and spinal cord was higher for the upregulated genes (40–42%), compared to overlaps in the downregulated genes (9–12%) (Figure 1f and h).

The overrepresentation test (ORT; BP complete; FDR < 0.05), using the lists of genes that were upregulated or downregulated in the spinal cord at the advanced-symptomatic stage, identified the immune and nervous systems as the top biological systems significantly enriched for the upregulated (Figure 1g) and downregulated genes (Figure 1h), respectively.

Since the top largest enriched GO BP term for downregulated genes was 'nervous system development', we further dissected transcriptional regulation of the developmental and repair processes in the CNS that may be altered by WNV infection using the gProfiler (Raudvere et al., 2019). Analysis of the upregulated genes in the cerebellum and spinal cord revealed upregulation of the developmental processes with major ‘spikes’ (i.e., most significant enrichment) related to the immune system development, response to wounding, and wound healing (Figure 2a and c). Both the cerebellum and spinal cord also displayed upregulation of the gliogenesis/astrocyte development and extracellular matrix/tissue remodeling. In contrast, developmental and/or repair processes related to neurons, their cellular components, and respective anatomical structure/network development were significantly enriched for the downregulated genes, with major ‘spikes’ related to neuron/neuron projection development and axon development/guidance (Figure 2b and d). This suggested that there was a downregulation rather than activation of the developmental processes related to neurons. However, given the detected spikes in regulation of axon development/guidance, we asked whether the molecular environment for axonal regeneration was inhibitory or permissive based on the up- or downregulation of expression of established axon growth inhibitory or permissive molecules (Anderson et al., 2016). Of 59 axon-growth-modulating molecules analyzed, we found changes in gene expression for 6 inhibitory and 16 permissive molecules (Figure 2e). Taking into account the fold changes for these 22 axon growth genes, when (i) downregulation of permissive or upregulation of inhibitory molecules would indicate the inhibitory molecular environment for axonal regeneration, while (ii) downregulation of inhibitory or upregulation of permissive molecules would indicate the permissive environment, we found a progressive trend which was skewed to the axon growth permissive environment.

Figure 2

Download asset Open asset

Next, we used the Reactome Knowledgebase (Fabregat et al., 2018) to identify the biological pathways enriched in immune and neural system domains. The Reactome Knowledgebase currently has 27 annotated biological domains, including the immune and neuronal systems, that are visualized as the pathway hierarchy bursts (Figure 3a and b) and can be overlaid by the identified significantly affected pathways to display their coverage.

Figure 3

Download asset Open asset

As expected based on the identification of the immune system as the top biological system significantly enriched for the upregulated genes (Figure 1g), Reactome pathways associated with this system showed significant enrichment only for upregulated genes (Figure 3c, d, g, and h), but not for downregulated genes (at both early-symptomatic and advanced-symptomatic stages of WNV-ND, and in both the cerebellum and spinal cord). Major immune system subdomains such as the innate immune system (including the toll-like receptor cascades), cytokine signaling (including interferon signaling and signaling by interleukins), and the adaptive immune system (including class I MHC mediated antigen processing and presentation, and TCR/BCR receptor signaling) were identified as significantly enriched in both the cerebellum and spinal cord of NHPs at the symptomatic stages of WNV-ND (Figure 3c, d, g, and h). The only difference between the cerebellum and spinal cord in terms of activated immune pathways was that MAP kinase activation reached statistical significance only in the spinal cord at the advanced-symptomatic stage of WNV-ND (compare Figure 3d and h; pathway node #12).

As expected based on the identification of the neuronal system as the top biological system significantly enriched for the downregulated genes (Figure 1h), Reactome pathways associated with this system showed significant enrichment only for downregulated genes (Figure 3e, f, i, and j), but not for upregulated genes (at both early-symptomatic and advanced-symptomatic stages of WNV-ND, and in both the cerebellum and spinal cord). The neuronal system domain was most affected at the advanced-symptomatic stage in both the cerebellum and spinal cord (compare Figure 3e, f, i, and j), coinciding with the progression to severe neurological signs at this stage (Maximova et al., 2014). We noted that the Reactome’s neuronal system domain is limited only to pathways associated with the neurotransmission and does not cover the complexity of the nervous system. Nonetheless, it enabled identification of transmission across chemical synapses as a common significantly affected pathway in the cerebellum and spinal cord (Figure 3f, and j; pathway node #2). The only difference between the cerebellum and spinal cord in terms of affected neuronal pathways was that the protein–protein interactions at synapses reached the statistical significance only in the spinal cord at the advanced-symptomatic stage of WNV-ND (compare Figure 3f and j; pathway node #8).

Together, these results show that WNV infection induces progressive and CNS-structure-dependent gene expression changes that manifest in transcriptional upregulation of biological pathways ascribed to the immune system gene ontologies and downregulation of functions related to the neuronal system biological domain. Therefore, a larger overlap in the upregulated genes compared to downregulated genes between the cerebellum and spinal cord may be explained by a higher similarity in the immune responses, but a site-specific impairment of distinct neural functions carried by these two CNS structures. In addition, we show that WNV infection alters transcriptional regulation of developmental and repair processes in the CNS, which are associated with upregulation of the immune system development/wound healing and downregulation of the neuronal/axonal regenerative processes. However, further dissection of the molecular environment associated with axonal regeneration also indicated a trend that was more skewed to the axon growth permissive rather than inhibitory environment, suggesting initiation of the neuronal network repair programs.

Altered CNS homeostasis in WNV-ND is defined by the large and subtle coordinated transcriptional shifts regulating the immune and nervous system, respectively

To further analyze the magnitude of transcriptional dysregulation within the CNS caused by WNV infection, we employed a statistical functional enrichment method that takes into account the fold changes associated with DEGs (Mi et al., 2019). This method has an advantage of detecting the coordinated changes in expression of groups of genes that are shifted from the overall distribution. Such coordinated transcriptional shifts can be subtle and elude other functional annotation methods (Clark et al., 2003; Mootha et al., 2003). We reasoned that this method could provide a better understanding of the magnitude of transcriptional dysregulation of the nervous system that appeared relatively subtle at the level of pathway analysis (Figure 3). Since the pathway level analysis showed that the immune and neuronal domains were affected the most at the advanced-symptomatic stage of WNV-ND, we analyzed genes that were significantly differentially expressed in WNV-infected cerebellum or spinal cord at this time point, together with their expression values (FC over mock). The results of this test, as expected, revealed coordinated large and subtle transcriptional shifts regulating the immune and nervous system, respectively. The large coordinated shifts (i.e., shifts to the greater fold change values from the overall distribution) were associated with the host defense response to viral infection and innate/adaptive immune responses (Figure 4a–d). In contrast, the subtle coordinated shifts (i.e., shifts to the smaller fold change values from the overall distribution) were associated with the cellular compartments of neurons (i.e., neuronal cell body, dendrite, axon, and glutamatergic and GABA-ergic synapses) and functions related to neurons (i.e., chemical synaptic transmission, neurotransmitter secretion, action potential, ensheathment of neurons, microtubule cytoskeleton organization, and synapse organization) (Figure 4e–h). The complete results are provided in Supplementary file 1 (‘Coordinated transcriptional shifts’).

Figure 4

Download asset Open asset

These results demonstrate that WNV infection induces coordinated shifts in the transcriptional regulation of the CNS homeostasis and that these shifts differ in magnitude and are large when associated with the immune system and subtle when related to the nervous system.

Cell death processes in WNV-ND are regulated bi-directionally and magnitude of neuron cell death regulation does not exceed that of lymphocytes

Since the subtle coordinated transcriptional shifts in WNV-ND were associated with neurons and their functions (Figure 4e–h), but large shifts also included the regulation of cell death (Figure 4b and d), we next used the gProfiler (Raudvere et al., 2019) to simultaneously dissect the GO terms of BP and KEGG pathways (Kanehisa et al., 2019) to gain insight into the magnitude and direction of the regulation of cell death processes at the advanced-symptomatic stage of WNV-ND, with a particular focus on neurons.

Analysis of the genes that were downregulated in the cerebellum or spinal cord did not reveal any significantly enriched terms related to cell death, apoptosis, necroptosis, or necrosis. In contrast, analysis of the genes that were upregulated in the cerebellum or spinal cord showed significant (FDR < 0.05) enrichment of terms associated with regulation of cell death. Further dissection of the transcriptional regulation of cell death processes in WNV-ND revealed that: (i) cell death processes in WNV-ND were regulated bi-directionally (i.e., enriched GO BP terms included both negative and positive regulation) (Figure 5a); (ii) negative regulation of cell death appeared to outweigh the positive regulation based on the larger number of annotated genes in both the cerebellum and spinal cord (Figure 5a), although the FDR-adjusted p-values for the negative and positive regulation were both highly significant and comparable (in the range of 12.2–18.7 negative log10 of adjusted p-value); (iii) regulation of the following types of the cell death was significantly enriched: necroptosis, apoptosis, and necrosis (Figure 5a); (iv) regulation of the cell death processes was significantly enriched in two types of cells: neurons and lymphocytes (Figure 5a); (v) 53% of genes annotated to the ‘regulation of neuron death’ were commonly upregulated in the cerebellum and spinal cord, with only 9% of annotated genes exclusively upregulated in the cerebellum and 38% exclusively upregulated in the spinal cord (Figure 5b); (vi) functional characterization of the common gene set annotated to regulation of neuron death showed enrichment for both negative and positive regulation with a similar level of statistical significance, but genes that were regulated exclusively in the cerebellum contributed more to positive regulation while genes exclusively regulated in the spinal cord more significantly contributed to the negative regulation (Figure 5b); (vii) the percent of genes that became upregulated relative to the total number of genes known to be involved in regulation of the neuron death (i.e., term size) was about twice as a small than that for genes involved in regulation of the apoptotic process in lymphocytes, in both the cerebellum and spinal cord (Figure 5c); and (viii) regulation of the lymphocyte apoptotic processes was a region-specific (in the cerebellum regulation of apoptosis was increased in T cells, while in the spinal cord regulation of apoptosis was increased in B cells) (Figure 5c).

Figure 5

Download asset Open asset

Taken together, these results demonstrate that activation of transcriptional regulation of the cell death in WNV-ND is (i) bi-directional (i.e., concurrently positive and negative), (ii) region-specific (i.e., skewed to a positive regulation in the cerebellum but negative regulation in the spinal cord), and (iii) cell-type-specific (i.e., major cell types with increased regulation of the cell death processes were neurons and lymphocytes). Importantly, the magnitude of cell death regulation in neurons was about twofold less than that in lymphocytes, suggesting that neurons are not a major cell type with increased transcriptional regulation of cell death processes. Instead, it appears that the major cell type with increased regulation of cell death processes was represented by the lymphocytes (T and B cells that are expected to infiltrate the CNS parenchyma as a part of the adaptive immune response to viral infection [Maximova and Pletnev, 2018]). Strikingly, regulation of the apoptotic processes in lymphocytes was region-specific, with increased regulation of T-cell apoptosis in the cerebellum but not in the spinal cord and vice versa – increased regulation of B-cell apoptosis in the spinal cord but not in the cerebellum, suggesting differential lymphocytic responses to WNV infection in these two CNS regions.

WNV-induced transcriptional changes correspond to activation of non-neuronal multicellular responses to infection and disruption of the structural integrity and function of infected neurons

Large coordinated transcriptional shifts correspond to WNV-induced activation of non-neuronal multicellular responses to infection by CNS-intrinsic and CNS-extrinsic cells

Based on the identified large coordinated shifts in transcriptional regulation of the CNS homeostasis in WMV-ND (Figure 4a–d; Supplementary file 1 [‘Coordinated transcriptional shifts’]), we selected three major BPs and examined the infected tissues by immunohistochemistry, using appropriate protein markers. These processes were: (i) reactive astrocytosis; (ii) microglia/macrophage activation and migration; and (iii) lymphocytic infiltration and migration (Table 1).

Table 1

Biological process	Enriched GO terms	Protein marker IR*	Cellular morphology/topology/function	Representative images
Reactive astrocytosis	glial cell activation (GO:0061900); regulation of cell migration (GO:0030334)	GFAP ↑	Astrocytes: hypertrophy of somata and processes; neuron-centripetal migration; perineuronal topology	Figure 6a
Microglia/macrophage activation and migration	innate immune response (GO:0045087); glial cell activation (GO:0061900); regulation of cell migration (GO:0030334); regulation of microglial cell activation (GO:1903978); macrophage activation (GO:0042116); regulation of macrophage chemotaxis (GO:0010758); phagocytosis (GO:0006909); lysosome (GO:0005764)	CD68 ↑	Microglial cells: hypertrophy of somata and processes; neuron-centripetal migration; perineuronal topology; phagocytic activity	Figure 6b
Lymphocytic activation, migration, and infiltration	adaptive immune response (GO:0002250); lymphocyte activation (GO:0046649); lymphocyte migration (GO:0072676); B cell mediated immunity (GO:0019724); regulation of T cell mediated immunity (GO:0002709); regulation of CD4-positive, alpha-beta T cell activation (GO:2000514); regulation of CD8-positive, alpha-beta T cell activation (GO:2001185)	CD20 ↑	B cells: leptomeningeal infiltration in the cerebellar cortex; perivascular infiltration if the spinal cord gray matter; no parenchymal migration	Figure 6c
		CD4 ↑	T helper cells: leptomeningeal infiltration; perivascular infiltration; minimal parenchymal migration	Figure 6c
		CD8 ↑	Cytotoxic T lymphocytes: leptomeningeal infiltration; perivascular infiltration; parenchymal infiltration; neuron-centripetal migration; perineuronal topology	Figure 6c
Virus infection of specific neuronal cell types	negative regulation of viral process (GO:0048525); regulation of viral entry into host cell (GO:0046596)	WNV-Ag positivity	Purkinje cells (PCs; CALB-positive): presence of viral antigens in neuronal perikarya and processes Spinal motor neurons (SMNs; ChAT-positive): presence of viral antigens in neuronal perikarya and processes	Figure 7a and b Figure 7c and d
Disruption of structural integrity and function of infected neurons	neuronal cell body (GO:0043025); somatodendritic compartment (GO:0036477); dendrite (GO:0030425); regulation of cytosolic calcium ion concentration (GO:0051480); cellular calcium ion homeostasis (GO:0006874); calcium ion-regulated exocytosis of neurotransmitter (GO:0048791); neurotransmitter secretion (GO:0007269); synapse (GO:0045202); neuron to neuron synapse (GO:0098984); presynapse (GO:0098793); asymmetric synapse (GO:0032279); distal axon (GO:0150034); glutamatergic synapse (GO:0098978); GABA-ergic synapse (GO:0098982); synapse organization (GO:0050808); synaptic signaling (GO:0099536); chemical synaptic transmission (GO:0007268)	CALB ↓	Purkinje cells: decrease of CALB-IR in somatodendritic compartments of WNV-infected PCs; disturbance of cellular calcium homeostasis in infected PCs	Figure 7b
		ChAT ↓	Spinal motor neurons: decrease of ChAT-IR in SMN cytoplasm and synapses innervating SMNs; disruption of neurotransmitter secretion in SMNs and their afferent innervation	Figure 7d
		SYP ↓	Decrease in IR for structural constituent of synaptic vesicles (SYP-IR) in the cerebellar cortex and spinal cord. Disruption of structural integrity of presynaptic compartments. Cerebellar cortex: putatively affected synapses innervating PCs: CF-PC and PF-PC (asymmetric-glutamatergic) in ML; BC-PC (GABA-ergic) in ML and PCL; and SC-PC (GABA-ergic) in ML Spinal gray matter: putatively affected synapses innervating SMNs: cholinergic C- boutons (ChAT-IR); asymmetric-glutamatergic synapses from descending tracts; and inhibitory synapses from local inhibitory neuron networks	Figure 8a and c Figure 8a Figure 7d and Figure 8c
	neuronal cell body (GO:0043025); somatodendritic compartment (GO:0036477); dendrite (GO:0030425); postsynapse (GO:0098794); microtubule cytoskeleton organization (GO:0000226)	MAP2 ↓	Cerebellar cortex and spinal gray matter: loss of postsynaptic cellular targets for asymmetric (axodendritic) and symmetric (axosomatic) synapses; disruption of integrity of microtubule cytoskeleton	Figure 8b and d

1. ^*IR, immunoreactivity. ↑ - increased compared to mock. ↓ - decreased compared to mock. GFAP, glial fibrillary acidic protein. WNV-Ag, WNV antigens. CALB, calbindin D28k. ChAT, choline acetyltransferase. SYP, synaptophysin. MAP2, microtubule associated protein 2. ML. molecular layer. PCL, Purkinje cell (PC) layer. CF-PC, climbing fiber to PC synapses. PF-PC, parallel fiber to PC synapses. BC-PC basket cell to PC synapses. SC-PC, stellate cell to PC synapses.

Immunoreactivity for glial fibrillary acidic protein (GFAP) was increased in the cerebellum and spinal cord in advanced WNV-ND and showed hypertrophy of somata and processes of astrocytes and their neuron-centripetal migration and perineuronal topology, all consistent with reactive astrocytosis (response of CNS-intrinsic glial cells to damage and disease [Burda and Sofroniew, 2014], including flavivirus infections [Maximova and Pletnev, 2018]; Figure 6a). These changes in astrocyte morphology and topology corresponded to the significantly enriched GO terms related to glial cell activation and regulation of cell migration (Table 1).

Figure 6

Download asset Open asset

Strong immunoreactivity for the lysosomal-associated membrane protein CD68 (undetectable under normal physiological conditions) was observed in the cerebellum and spinal cord in advanced WNV-ND, consistent with microglial/macrophage activation. Activated/reactive microglial cells displayed neuron-centripetal migration and perineuronal topology (Figure 6b). These changes in microglial cell morphology and topology corresponded to the significantly enriched GO terms related to the innate immune response, regulation of microglial cell/macrophage activation, regulation of macrophage chemotaxis, and phagocytosis (Table 1).

Immunoreactivity for the lymphocytic protein markers CD20, CD4, and CD8 was detected in the cerebellum and spinal cord in advanced WNV-ND, consistent with infiltration of the CNS by peripheral CNS-extrinsic immune cells (Figure 6c). CD20 B cells were detected in the leptomeninges but not in the parenchyma of the cerebellum and only occasionally in perivascular sites in the spinal cord, suggesting limited migration of these cells beyond the initial sites of infiltration. CD4 T cells were detected mostly in the leptomeninges and at perivascular sites in both the cerebellum and spinal cord, and minimally in the parenchyma, also suggesting limited parenchymal migration of these cells. In contrast, numerous CD8 T cells were detected in leptomeningeal, perivascular, and parenchymal locations. Moreover, CD8 T cells displayed a neuron-centripetal migration and perineuronal topology. The infiltration of the CNS by these lymphocytic subtypes and differential regulation of their migration patterns corresponded to the significantly enriched GO terms related to the adaptive immune response and lymphocyte activation/migration (Table 1). In addition, the presence of the perivascular infiltration by B cells in the spinal cord gray matter, but not in the cerebellar cortex may be associated with the increased transcriptional regulation of B-cell apoptosis in the spinal cord but not in the cerebellum (Figure 5c), thus providing further support to the notion of differential lymphocytic responses to WNV infection in these two CNS regions.

Subtle coordinated transcriptional shifts correspond to alteration of the structural integrity and function of infected neurons

Based on the identified subtle coordinated shifts in transcriptional regulation of the CNS homeostasis in WMV-ND (Figure 4e–h; Supplementary file 1 [‘Coordinated transcriptional shifts’]), we selected two major BPs and examined WNV-infected tissues by immunohistochemistry, using appropriate protein markers. These processes were: (i) virus infection of specific neuronal cell types; and (ii) disruption of structural integrity and function of infected neurons (Table 1).

To identify the specific neuronal cell types supporting WNV replication in the cerebellum and spinal cord, we used brightfield immunohistochemistry for WNV-antigens followed by double immunofluorescent staining for WNV-antigens and appropriate cell type markers. Colorimetric immunoreactivity for WNV-antigens (brown) was detected in Purkinje cells in the cerebellum (Figure 7a) and motor neurons in the ventral horns of the spinal cord (Figure 7c) in advanced WNV-ND, consistent with typical infection of these neuronal types by flaviviruses. Corresponding to the host response to viral infection, the significantly enriched GO terms were related to the regulation of viral process and viral entry into host cell (Table 1).

Figure 7

Download asset Open asset

Strikingly, double immunofluorescent staining for WNV-antigens and a calcium-binding protein calbindin D28k (CALB, protein marker for Purkinje cells) revealed that immunoreactivity for CALB was diminished in Purkinje cells that were WNV-infected (i.e., WNV+++/CALB+ PCs), compared to highly intense CALB-immunoreactivity in non-infected Purkinje cells (i.e., WNV-/CALB+++ PCs) (Figure 7b). CALB-immunoreactivity was markedly reduced in the somata of WNV-infected Purkinje cells and almost entirely lost from their dendritic trees extending into the molecular layer of the cerebellar cortex (Figure 7b). This finding is consistent with a disturbance of cellular calcium homeostasis and loss of the integrity of somatodendritic compartments of Purkinje cells induced by WNV infection at both protein and gene expression levels (Table 1).

Similar to the staining pattern with WNV-antigens and a neuron-specific marker observed in the Purkinje cells, double immunofluorescent staining for WNV-antigens and spinal motor neuron-specific marker choline acetyltransferase (ChAT) revealed that ChAT-immunoreactivity was markedly reduced in the somata and processes of spinal motor neurons that were WNV-infected (i.e., WNV+++/ChAT+ SMNs), compared to highly intense ChAT-immunoreactivity in non-infected spinal motor neurons (i.e., WNV-/ChAT+++ SMNs) (Figure 7d). Strikingly, normally abundant ChAT+++ cholinergic presynaptic C-boutons innervating the somata and proximal dendrites of spinal motor neurons (see normal SMN example in mock-inoculated control) almost completely disappeared from the SNN somata with a few boutons remaining adjacent to their proximal dendrites (Figure 7d). Interestingly, focal accumulations of WNV+++ granules in the proximal dendrites of SMNs (postsynaptic sites) could be seen immediately adjacent to the remaining ChAT+ cholinergic presynaptic C-boutons, suggesting trans-synaptic spread and/or local postsynaptic replication of WNV (Figure 7d, WNV+++/ChAT+ SMN). These findings are consistent with disruption of the somatodendritic and synaptic structural integrity in WNV-infected SMNs, their neurotransmitter secretion, and afferent innervation (Table 1), which collectively would impair their function.

Immunoreactivity for the synaptic marker protein synaptophysin (SYP) was greatly depleted in advanced WNV-ND compared to mock, encompassing all layers of the cerebellar cortex (i.e., molecular layer, Purkinje cell layer, and granule cell layer) (Figure 8a) and gray matter of the spinal cord (Figure 8c).

Figure 8

Download asset Open asset

Immunoreactivity for microtubule-associated protein 2 (MAP2, a marker for somatodendritic [and thus, postsynaptic] compartments of neurons) was also markedly reduced, confirming the loss of neuronal somatodendritic integrity in Purkinje cells and granule cells in the cerebellum (Figure 8b) and spinal motor neurons in the ventral horns of the spinal cord (Figure 8d). These findings are consistent with the loss of postsynaptic cellular targets for asymmetric (axodendritic) and symmetric (axosomatic) synapses and disruption of integrity of microtubule cytoskeleton in WNV-infected neurons in the cerebellum and spinal cord (Table 1).

Together, these alterations in somatodendritic and synaptic integrity of infected neurons corresponded to the significantly enriched GO terms related to the neuronal cell body, dendrites, microtubule cytoskeleton organization, and synapse (Table 1).

Based on these findings, we were interested in identifying a cellular compartment of neurons that was a major target of differential transcriptional regulation in advanced WNV-ND. Ranking the neuronal CC GO terms based on their enrichment values, we identified the synapse as the most significantly dysregulated neuronal cell compartment (Figure 8e and f; note the smallest FDR value for the GO CC term synapse (GO:0045202)). This phenomenon was observed in both the cerebellum and spinal cord, albeit the spinal cord had about 2.5 times more dysregulated genes annotated to the GO CC term synapse compared to the cerebellum. Within the synapse, significantly enriched child GO terms included the presynapse (GO:0098793), synaptic membrane (GO:0097060), postsynapse (GO:0098794), and ion channel complex (GO:0034702). Further defining the types of synapses affected in WNV-ND, glutamatergic (i.e., asymmetric, axo-dendritic, and excitatory) synapses (GO:0098978) and GABA-ergic (i.e., symmetric, axo-somatic, and inhibitory) synapses (GO:0098982) emerged as significantly transcriptionally dysregulated in both the cerebellum and spinal cord (Figure 8e and f).

Topography of the observed loss of the presynaptic compartments (Figure 7d; Figure 8a and c) and postsynaptic cellular targets (Figure 8b and d) in WNV-ND suggests that (i) putatively affected synapses in the cerebellar cortex may include synapses innervating Purkinje cells (e.g., climbing fiber-PC and parallel fiber-PC [glutamatergic, asymmetric, axo-dendritic, and excitatory] in the molecular layer, Basket cell-PC [GABA-ergic and inhibitory] in the molecular and PC layer, and stellate cell-PC [GABA-ergic and inhibitory] in the molecular layer); and (ii) putatively affected synapses in the spinal cord may include synapses innervating the spinal motor neurons (e.g., cholinergic C-boutons, glutamatergic synapses from descending tracts, and inhibitory synapses from local inhibitory neuron networks) (Table 1).

Taken together, these results reveal numerous alterations of the structural integrity and function of infected neurons and suggest that the synapse is a major target of transcriptional dysregulation in WNV-infected cerebellar and spinal neurons, which is supported by evidence of changes in synaptic integrity at the level of protein expression.

WNV infection disrupts transcriptional regulation of synaptic organization and function

Having identified the synapse as the top transcriptionally dysregulated neuronal cell compartment in advanced WNV-ND, we next asked what specific synaptic subcompartments (hereafter referred as synaptic location) and BP (hereafter referred as function) were affected. We also sought to infer whether there was a loss or gain of synaptic functions based on the synaptic genes that were upregulated or downregulated during the symptomatic stages of WNV-ND. For this, we used the genomic analysis tool SynGO, a knowledge base that accumulates research about the synapse biology and includes about 3000 expert-curated GO annotations for 1112 synaptic genes (Koopmans et al., 2019). SynGO synaptic ontology can be visualized as sunburst diagrams for synaptic location and function (Figure 9a and b, respectively). To determine a precise impact of WNV infection on transcriptional regulation of the synapses, we used high stringency SynGO settings, where synaptic gene annotations are based exclusively on published experimental evidence from neuronal biological systems.

Figure 9

Download asset Open asset

SynGO enrichment analysis of the DEGs annotated to the CC GO term synapse (GO:0045202) (Supplementary file 2, ‘Synaptic genes differentially expressed during WNV-ND’) at the advanced-symptomatic stage of WNV-ND revealed extensive transcriptional dysregulation of synapses in both the cerebellum and spinal cord (Figure 9c–f). In the cerebellum, the postsynapse was the most significantly affected synaptic subcompartment based on the lowest Q-value (Figure 9c; Table 2), followed by the presynapse and synaptic cleft, while the most affected synaptic functions were the synapse organization (regulation of synapse assembly/synapse adhesion between pre- and postsynapse), trans-synaptic signaling (modulation of chemical synaptic transmission), processes in presynapse (synaptic vesicle cycle), and processes in postsynapse (regulation of postsynaptic membrane potential/neurotransmitter receptor levels) (Figure 9d; Table 2). Compared to the cerebellum, the spinal cord had many more synaptic locations (Figure 9e; Table 2) and functions (Figure 9f; Table 2) affected and more enriched (including many sublocations and specific synaptic functions). These included synaptic locations such as the presynapse, synaptic membrane, synaptic cleft, and postsynapse; and synaptic functions such as synapse organization, trans-synaptic signaling, processes in pre- and postsynapse, metabolism, and transport. Full results of the SynGO enrichment analysis for all differentially expressed synaptic genes can be found in Supplementary file 3 (‘SynGO enrichment analysis results’).

Table 2

GO term ID	GO domain	GO term name - hierarchical structure	Q-value
			Spinal cord	Cerebellum
GO:0045202	CC	synapse	7.13E-113	3.56E-35
GO:0097060	CC	├─ synaptic membrane	6.07E-03
GO:0098793	CC	├─ presynapse	4.20E-51	4.19E-17
GO:0099523	CC	│ ├─ presynaptic cytosol	7.07E-05	1.79E-03
GO:0048786	CC	│ ├─ presynaptic active zone	1.52E-15
GO:0048787	CC	│  │ └─ presynaptic active zone membrane	3.25E-11	1.08E-06
GO:0098833	CC	│ ├─ presynaptic endocytic zone	6.87E-03
GO:0008021	CC	│ ├─ synaptic vesicle	6.06E-05
GO:0098992	CC	│ ├─ neuronal dense core vesicle	4.44E-03
GO:0042734	CC	│ ├─ presynaptic membrane	4.46E-21	3.52E-10
GO:0043083	CC	├─synaptic cleft	5.66E-07	2.40E-05
GO:0098794	CC	├─ postsynapse	5.38E-80	4.68E-24
GO:0099524	CC	│ ├─ postsynaptic cytosol	1.61E-02
GO:0099571	CC	│ ├─ postsynaptic cytoskeleton	2.60E-05
GO:0099572	CC	│ ├─ postsynaptic specialization	4.46E-44	4.94E-11
GO:0014069	CC	│  │ ├─ postsynaptic density	1.73E-34	3.60E-08
GO:0045211	CC	│ ├─ postsynaptic membrane	1.12E-21	6.60E-10
SYNGO:postsyn_ser	CC	│  └─ postsynaptic SER	3.30E-04
SYNGO:synprocess	BP	process in the synapse	1.30E-87	2.74E-35
SYNGO:presynprocess	BP	├─ process in the presynapse	4.96E-21	4.11E-07
GO:0099509	BP	│ ├─ regulation of presynaptic cytosolic calcium levels	7.25E-04
GO:0099504	BP	│ ├─ synaptic vesicle cycle	5.53E-21	2.01E-04
GO:0016079	BP	│  │ ├─ synaptic vesicle exocytosis	9.78E-13	3.58E-03
GO:0099502	BP	│  │ │ ├─ calcium-dependent activation of synaptic vesicle fusion	4.74E-04
GO:2000300	BP	│  │ │ ├─ regulation of synaptic vesicle exocytosis	1.82E-04
GO:0016082	BP	│  │ │ └─ synaptic vesicle priming	4.88E-05
GO:0048488	BP	│  │ └─ synaptic vesicle endocytosis	6.15E-06
SYNGO:postsynprocess	BP	├─ process in the postsynapse	5.44E-17	8.45E-06
GO:0099566	BP	│ ├─ regulation of postsynaptic cytosolic calcium levels	2.92E-03
GO:0060078	BP	│ ├─ regulation of postsynaptic membrane potential	2.75E-07	1.63E-03
GO:0099072	BP	│ ├─ regulation of postsynaptic membrane neurotransmitter receptor levels	6.34E-08	1.10E-02
GO:0099645	BP	│  │ ├─ neurotransmitter receptor localization to postsynaptic specialization membrane	1.70E-05
GO:0099149	BP	│  │ └─ regulation of postsynaptic neurotransmitter receptor endocytosis	9.33E-05
GO:0099537	BP	│ ├─ trans-synaptic signaling	1.47E-20	8.11E-13
GO:0007268	BP	│  │ └─ chemical synaptic transmission	1.11E-17	1.43E-10
GO:0050804	BP	│  │ ├─ modulation of chemical synaptic transmission	8.44E-11	1.63E-07
GO:0099171	BP	│  │ │ └─ presynaptic modulation of chemical synaptic transmission	9.43E-05	7.08E-04
GO:0099170	BP	│  │  └─ postsynaptic modulation of chemical synaptic transmission	2.89E-04
GO:0050808	BP	├─ synapse organization	4.53E-45	4.99E-14
GO:0099173	BP	│ ├─ postsynapse organization	5.47E-15	1.91E-02
GO:0099010	BP	│  │ └─ modification of postsynaptic structure	1.92E-05
GO:0099181	BP	│  │ ├─ structural constituent of presynapse	1.46E-06
GO:0098882	BP	│  │ │ ├─ structural constituent of active zone	9.80E-06
GO:0099186	BP	│  │ └─ structural constituent of postsynapse	2.40E-03
GO:0099560	BP	│ ├─ synapse adhesion between pre- and post-synapse	8.52E-06	4.79E-03
GO:0007416	BP	│ ├─ synapse assembly	1.75E-10	2.01E-04
GO:0099054	BP	│  │ ├─ presynapse assembly	8.02E-03
GO:0051963	BP	│  │ ├─ regulation of synapse assembly	4.99E-07	1.19E-04
GO:1905606	BP	│  │ │ ├─ regulation of presynapse assembly	1.37E-04	2.76E-04
GO:0097107	BP	│  │ │ └─ postsynaptic density assembly	1.03E-02
GO:0060074	BP	│ ├─ synapse maturation	6.28E-05
GO:0099188	BP	│  └─ postsynaptic cytoskeleton organization	2.89E-04	1.46E-03
SYNGO:metabolism	BP	├─ metabolism	1.69E-02
SYNGO:catabolic_postsynapse	BP	│  │ └─ protein catabolic process at postsynapse	2.40E-03
SYNGO:transport	BP	└─ transport	1.48E-03
GO:0098887	BP	│  └─ neurotransmitter receptor transport, endosome to postsynaptic membrane	2.40E-03

1. Note: The top-level GO terms 'synapse' and 'process in the synapse' and major successive terms for synaptic subcompartments and processes are highlighted for the cerebellum and spinal cord at the advanced-symptomatic stage of WNV-ND. CC, cellular component; BP, biological process.

Analysis of upregulated versus downregulated synaptic genes in the cerebellum at the early-symptomatic stage of WNV-ND showed that 80% of synaptic genes were upregulated and 20% were downregulated (Table 3). The upregulated synaptic genes were postsynaptic while the downregulated genes were presynaptic. Transcriptional regulation in the spinal cord at this early-symptomatic stage of WNV-ND was characterized by a similar number of upregulated (51.2%) and downregulated (48.8%) synaptic genes, which were annotated to both the presynapse and postsynapse and had functions in synapse organization and chemical synaptic transmission (Table 3).

Table 3

	Number of upregulated genes (percent)	Enriched GO terms	Number of downregulated genes (percent)	Enriched GO terms
Early-symptomatic
Cerebellum	24 (80%)	postsynaptic specialization	6 (20%)	presynapse
Spinal cord	44 (51.2%)	presynapse; postsynapse; postsynaptic specialization; synapse organization; synapse assembly; chemical synaptic transmission	42 (48.8%)	postsynapse; postsynaptic specialization; synapse organization; regulation of postsynaptic membrane neurotransmitter receptor levels; modulation of chemical synaptic transmission
Advanced-symptomatic
Cerebellum	67 (45.9%)	synapse organization; regulation of synapse assembly; presynaptic active zone membrane; presynaptic modulation of chemical synaptic transmission; postsynaptic specialization; metabolism	79 (54.1%)	regulation of synapse organization; regulation of presynapse assembly; presynaptic cytosol; presynaptic active zone membrane; regulation of presynaptic membrane potential; postsynaptic actin cytoskeleton organization; postsynaptic specialization; regulation of postsynaptic membrane neuro-transmitter receptor levels; regulation of postsynaptic membrane potential; synapse adhesion between pre- and postsynapse; regulation of synaptic vesicle cycle; modulation of chemical synaptic transmission
Spinal cord	109 (28.5%)	synapse organization; regulation of synapse assembly; synaptic vesicle membrane; synapse adhesion between pre- and post-synapse; synaptic cleft; presynaptic active zone membrane; presynaptic modulation of chemical synaptic transmission; postsynaptic specialization; synaptic vesicle cycle; modulation of chemical synaptic transmission	274 (71.5%)	regulation of synapse organization; synaptic vesicle membrane; neuronal dense core vesicle; regulation of synaptic vesicle cycle; synaptic vesicle neurotransmitter loading; synaptic vesicle priming; synapse adhesion between pre- and postsynapse; synaptic cleft; structural constituent of presynapse; presynaptic cytosol; presynaptic active zone cytoplasmic component; presynaptic active zone membrane; presynaptic endocytic zone; regulation of presynapse assembly; regulation of presynaptic cytosolic calcium levels; regulation of presynaptic membrane potential; presynaptic modulation of chemical synaptic transmission; structural constituent of postsynapse; regulation of postsynapse organization; postsynaptic cytosol; regulation of modification of postsynaptic actin cytoskeleton; modification of postsynaptic structure; postsynaptic specialization assembly; regulation of postsynaptic cytosolic calcium levels; regulation of calcium-dependent activation of synaptic vesicle fusion; protein catabolic process at postsynapse; regulation of postsynaptic membrane neurotransmitter receptor levels; regulation of postsynaptic neurotransmitter receptor activity; regulation of postsynaptic membrane potential; postsynaptic modulation of chemical synaptic transmission; transport; metabolism

At the advanced-symptomatic stage of WNV-ND, the number of downregulated synaptic genes increased relative to upregulated genes in both the cerebellum and spinal cord, with a higher percentage of downregulated genes in the spinal cord (71.5%). Nonetheless, about one third of the differentially expressed synaptic genes were upregulated at this stage in both the cerebellum and spinal cord (45.9% and 28.5%, respectively). The upregulated and downregulated synaptic genes shared their annotated location and function suggesting a bidirectional dysregulation of transcriptional synaptic homeostasis in WNV-ND, affecting the spinal cord to a higher degree compared to the cerebellum (Table 3). Full results of the SynGO enrichment analysis for the upregulated and downregulated synaptic genes can be found in Supplementary File S4 (‘SynGO enrichment: up- and downregulated synaptic genes’).

Taken together, these findings suggest that WNV infection disrupts transcriptional homeostasis of the CNS synapses and impairs neurotransmission.

WNV infection induces differential expression of pleiotropic genes with a triple immune, neural, and synaptic topology and functionality

Since functional analysis of genes upregulated at the presymptomatic (3 dpi) stage of WNV-ND showed that some of genes (MX1 and STAT1) have a dual immune-neural functionality (i.e., being pleiotropic), we next sought to determine whether the number of such genes increased as infection progressed. With a focus on the advanced-symptomatic stage of WNV-ND, that had the highest number of dysregulated genes, we identified overlaps in DEGs annotated to the immune or nervous system (hereafter referred as immune or neural DEGs, respectively) (Figure 10a).

Figure 10

Download asset Open asset

Relative to the total number of WNV-induced DEGs, the cerebellum had a higher percentage of immune DEGs (33.1%), compared to the spinal cord (18.6%), even though the absolute number of immune DEGs in the spinal cord was higher. Strikingly, despite more than twice the number of the neural DEGs induced in the spinal cord compared to the cerebellum (642 vs. 252), the percentage of neural DEGs relative to the total number of DEGs was almost identical in these two CNS structures (14.6% vs. 14.5%) (Figure 10a). These data suggest that specific CNS regions have differences in immune responses to WNV infection, but more uniform responses related to regulation of neural function.

We next identified 94 and 179 common immune-neural DEGs that were induced in the cerebellum and spinal cord, respectively (Figure 10a). This indicates that these DEGs have a dual immune-neural functionality, suggesting that they are pleiotropic (Radisky et al., 2009; Boulanger, 2009). Common immune-neural as well as exclusively immune or neural DEGs are listed in Supplementary file 5 (‘Venn-diagram results for immune and neural DEGs’).

We next asked what neuronal locations and functions were under transcriptional control by the WNV-induced pleiotropic genes. Somatodendritic compartment, neuron projection, and synapse were identified as the main neuronal targets of regulation by pleiotropic genes (ORT; CC GO domain; FDR < 0.05) (Table 4). Neuron projection development, regulation of the neuron death, microglial/astroglial cell activation/migration, axon guidance, synapse structure/activity, postsynapse to nucleus signaling, and synapse pruning were identified as the main functions regulated by pleiotropic genes (ORT; BP GO domain; FDR < 0.05) (Table 4). The complete functional annotations for the pleiotropic immune-neural genes in WNV-infected cerebellum and spinal cord are provided in Supplementary file 6 (‘ORT: Immune-neural pleiotropic DEGs’).

Table 4

GO domain	GO term	No. of genes	Gene symbols
Cerebellum
CC	somatodendritic compartment (GO:0036477)	19	MBP, HCLS1, ITGA4, PTX3, FCGR2B, ZBP1, CCR2, CYFIP1, CTSL2, CIB1, AIF1, SOS1, ANXA3, CD3E, S100B, FLNA, MAPK1, ITGA1, TGFB2
CC	neuron projection (GO:0043005)	19	MBP, EPHA2, HCLS1, ITGA4, FCGR2B, ZBP1, VIM, CCR2, BCL11B, CYFIP1, CTSL2, CIB1, ANXA3, CD3E, FLNA, MAPK1, ITGA1, TGFB2, ODZ1
CC	postsynapse (GO:0098794)	11	JAK2, HCLS1, C1QA, FCGR2B, ZBP1, CYFIP1, STAT3, SOS1, CD3E, FLNA, MAPK1
BP	neuron projection development (GO:0031175)	25	EPHA2, JAK2, IL6, HCLS1, ITGA4, ADM, RHOG, SHC1, BCL11B, CYFIP1, SDC4, LST1, DOK2, CXCR4, SEMA4A, CSF1R, SOS1, CD3E, S100B, HMGB1, GATA3, MAPK1, CXCL12, ITGA1, NCKAP1L
BP	positive regulation of neuron death (GO:1901216)	13	PTX3, C1QA, FCGR2B, IFNG, BAX, TYROBP, TNF, BCL2L11, GRN, TLR4, ITGA1, TGFB2, CCL3
BP	microglial cell activation (GO:0001774)	10	JAK2, C1QA, IFNG, TLR2, TYROBP, TNF, FPR2, AIF1, GRN, C5AR1
BP	astrocyte activation (GO:0048143)	7	C1QA, IFNG, TNF, FPR2, IL1B, GRN, C5AR1
BP	glial cell migration (GO:0008347)	4	P2RY12, CCL2, TGFB2, CCL3
BP	axon guidance (GO:0007411)	12	EPHA2, RHOG, SHC1, BCL11B, CYFIP1, DOK2, CXCR4, CSF1R, SOS1, GATA3, MAPK1, CXCL12
BP	postsynapse to nucleus signaling pathway (GO:0099527)	3	JAK2, STAT3, RELA
BP	regulation of synapse structure or activity (GO:0050803)	8	HCLS1, FCGR2B, TLR2, SEMA4D, CYFIP1, TNF, SLC7A11, SEMA4A
BP	positive regulation of glutamate receptor signaling pathway (GO:1900451)	3	IFNG, CCR2, CCL2
BP	synapse pruning (GO:0098883)	2	C1QA, C3
Spinal cord
CC	somatodendritic compartment (GO:0036477)	39	REG1A, HSP90AA1, HCLS1, DAB2IP, DHX36, ITGA4, PTX3, FCGR2B, ZBP1, CX3CL1, BCR, HSP90AB1, PTK2B, RARA, FZD3, PSEN1, CCR2, APP, CASP3, PAFAH1B1, DOCK10, ALCAM, ARHGEF7, CTSL2, CIB1, AIF1, PTK2, P2RX4, RTN4, KIF5B, PURA, ADAM10, SNAP25, DNM2, FLNA, MAPK1, NEDD4, BECN1, ITGA1
CC	neuron projection (GO:0043005)	46	REG1A, EPHA2, HSP90AA1, HCLS1, DAB2IP, DHX36, ITGA4, FCGR2B, ZBP1, CX3CL1, BCR, HSP90AB1, PTK2B, RARA, FZD3, PSEN1, VIM, CCR2, APP, LRRC7, BSG, PAFAH1B1, DOCK10, ALCAM, ARHGEF7, CTSL2, CIB1, ARF1, PTK2, ITGA3, CYFIP2, P2RX4, KIF5B, PURA, ADAM10, NCAM1, SNAP25, DNM2, FLNA, MAPK1, NEDD4, BECN1, ITGA1, CEP290, ODZ1, BRAF
CC	axon (GO:0030424)	25	REG1A, HSP90AA1, DAB2IP, DHX36, ITGA4, ZBP1, BCR, HSP90AB1, PTK2B, FZD3, PSEN1, APP, LRRC7, BSG, PAFAH1B1, ALCAM, CIB1, ITGA3, P2RX4, KIF5B, ADAM10, SNAP25, DNM2, FLNA, MAPK1
CC	synapse (GO:0045202)	37	HCLS1, C1QA, FCGR2B, ZBP1, RAB35, BCR, PTK2B, FZD3, PSEN1, FZD3, NCSTN, APP, LRRC7, STAT3, PAFAH1B1, DOCK10, LYN, ARF1, PTK2, ITGA3, CYFIP2, EFNB1, RELA, P2RX4, RTN4, PURA, ADAM10, PTPRZ1, SNAP25, DNM2, FLNA, MAPK1, NEDD4, PTN, FGFR2, PAK2, ANXA1
CC	glutamatergic synapse (GO:0098978)	14	BCR, PTK2B, FZD3, STAT3, LYN, ARF1, RELA, PURA, ADAM10, SNAP25, DNM2, FLNA, NEDD4, PAK2
BP	neuron projection development (GO:0031175)	50	EPHA2, HSP90AA1, HCLS1, IL6, DAB2IP, PIK3R1, ITGA4, SLC11A2, HES1, ADM, JUN, RAB35, HSP90AB1, SEC24B, PTK2B, FZD3, APP, SHC1, PTPN11, CASP3, LST1, BSG, PAFAH1B1, DOK2, DOCK10, ALCAM, LYN, RPS6KA5, PTK2, EIF2AK4, SRF, CYFIP2, SEMA4A, EFNB1, CSF1R, RTN4, KIF5B, PTPRZ1, NCAM1, HMGB1, DNM2, GATA3, MAPK1, CXCL12, NEDD4, PTN, ITGA1, STK4, FGFR2, PAK2
BP	microglial cell activation (GO:0001774)	14	TREM2, C1QA, IFNGR1, JUN, IFNG, CX3CL1, TLR2, APP, TYROBP, TNF, FPR2, AIF1, GRN, C5AR1
BP	astrocyte activation (GO:0048143)	11	TREM2, C1QA, IFNGR1, IFNG, PSEN1, APP, TNF, FPR2, IL1B, GRN, C5AR1
BP	regulation of glial cell migration (GO:1903975)	8	TREM2, CX3CL1, CCR2, GPR183, CSF1, P2RY12, P2RX4, CCL3
BP	regulation of synapse structure or activity (GO:0050803)	17	HCLS1, DAB2IP, DHX36, FCGR2B, TLR2, FZD3, APP, SEMA4D, TNF, PAFAH1B1, SLC7A11, PTK2, SEMA4A, ADAM10, DNM2, NEDD4, PTN
BP	synapse pruning (GO:0098883)	3	TREM2, C1QA, CX3CL1

Since functional characterization of WNV-induced pleiotropic immune-neural DEGs suggested that the synapse is one of their major targets in the neuron, we next determined the overlap between these pleiotropic DEGs and WNV-induced synaptic DEGs identified (Supplementary file 2, ‘Synaptic genes differentially expressed during WNV-ND’). To compare these genes, we focused on the advanced-symptomatic stage of WNV-ND in the spinal cord since it had the highest number of genes in both groups. A Venn diagram comparison identified an overlap containing 36 immune-neural-synaptic pleiotropic DEGs, of which 19 were upregulated and 17 were downregulated (Figure 9b; Supplementary file 7, ‘Immune-neural-synaptic pleiotropic genes’). Functional GO enrichment analysis of the upregulated genes in this overlap indicated that they had multiple roles in the positive regulation of the following functions: (i) immune/defense response to virus; (ii) nervous system development (neurogenesis/gliogenesis); and (iii) neuronal cell body/projections/synapses-specific processes (Figure 10c and d). These processes were also significantly enriched for the downregulated pleiotropic genes, with a trend of less significant (compared to upregulated genes) enrichment of immune processes and more significant enrichment of processes related to the organization, structure, and/or activity of the neuronal compartments (Figure 10e and f). Full results of functional GO enrichment analysis for the upregulated and downregulated pleiotropic immune-neural-synaptic genes are provided in Supplementary file 8 (‘gProfiler: Immune-neural-synaptic pleiotropic DEGs’).

Interestingly, interrogation of specific pleiotropic genes showed that the gene C1QA that encodes a major constituent of the complement system subcomponent C1, was upregulated and mapped to significantly enriched GO terms associated with (i) immune/defense response to virus infection, (ii) glial cell response, (iii) neuron differentiation, and (iv) synapse (Supplementary file 8, gProfiler: Immune-neural-synaptic pleiotropic DEGs’), as well as synaptic pruning (Table 4).

Since the synapse emerged as the neuronal cell compartment most significantly enriched for pleiotropic genes (Figure 10c–f), we applied synaptic GO enrichment analysis (SynGO) and found that both up- and downregulated genes had functionality in synapse organization and topology (i.e., localization of the gene product functions) at the presynapse and postsynapse (Figure 10g). Specifically, the GO term of chemical synaptic transmission was significantly enriched for upregulated pleiotropic genes, while the GO term of protein catabolic process at the postsynapse was significantly enriched for downregulated pleiotropic genes (Figure 10g). A complete list of the enriched synaptic GO terms is provided in Supplementary file 9 (‘Pleiotropic DEGs: Synaptic topology and functionality’).

Taken together, these results indicate that WNV infection of neurons induces differential expression of pleiotropic genes that have multiple functionalities. Strikingly, we found that in addition to their expected role in immune/defense responses during WNV-ND, these genes also regulated distinct functions in neurons and their synapses. This suggests that WNV infection of neurons disrupts transcriptional homeostasis of the immune-neural-synaptic axis with possible off-target effects of virus-induced immune responses on neural function.

Severe disease of the CNS due to infection with flaviviruses is rare, but the consequences for neural function can be devastating (reviewed in Maximova and Pletnev, 2018). Pathogenesis and outcome of CNS diseases caused by flaviviruses are complex and involve an interplay of at least three components: (i) virus properties (i.e., virulence, ability to establish productive replication in neurons, and spread by axonal transport); (ii) host defense responses and/or viral evasion of these responses; and (iii) changes in neural function with ensuing specific neurological impairments. Although we now possess a substantial body of knowledge about flavivirus biology and host immune responses to infection (reviewed in Pierson and Diamond, 2020), including immune responses that develop specifically within the CNS (reviewed in Maximova and Pletnev, 2018), our understanding of how they directly or indirectly influence neural function is incomplete. Here, we aimed to deconvolute the complex changes in CNS physiology that occur during flavivirus infection by examining differential regulation of BPs related to the immune and nervous systems in the NHP model of WNV-ND. We chose to focus on the cerebellum and spinal cord since these CNS structures appear to be among the most affected in both humans with WNV neuroinvasive disease (Kleinschmidt-DeMasters and Beckham, 2015; Sejvar, 2016; Lenka et al., 2019; Omalu et al., 2003; Cushing et al., 2004; Guarner et al., 2004; Armah et al., 2007; Hart et al., 2014) and NHP model of WNV-ND (Maximova et al., 2016; Maximova et al., 2014), causing ataxia and flaccid paralysis, respectively.

CNS neurons are terminally differentiated cells and cannot be replenished, thus virus infection of neurons present a challenge for both the nervous and immune system to coordinate virus clearance while protecting neurons and preserving neural function. While innate immune responses in the CNS are common during virus infection (Maximova and Pletnev, 2018; Griffin, 2011; Griffin, 2003), despite the concept of immune privilege status of the CNS (Ransohoff and Brown, 2012), these responses need to be tightly controlled, since excessive or chronic inflammation can be harmful and detrimental to neural function (Griffin, 2011; Ransohoff and Brown, 2012; Ransohoff and Cardona, 2010). Here we introduce the concept that virus infection of CNS neurons disrupts the immune-neural axis homeostasis due to activation of pleiotropic genes that can regulate both innate immunity and neural function. Pleiotropy describes a concept where one gene and its encoded protein can control disparate, apparently unrelated BPs (Radisky et al., 2009; Boulanger, 2009). This concept is especially intriguing in the context of regulation of critically important neural functions. Many proteins with established functions in the immune system are also expressed in the developing and adult nervous system. Remarkably, some pleiotropic proteins such as proinflammatory cytokines and proteins in the complement system and major histocompatibility complex are essential for the establishment, organization, function, and removal of synapses between neurons (Boulanger, 2009). We found that WNV infection of neurons in the cerebellar and spinal neural networks induces differential regulation (i.e., both up- and downregulation) of close to 200 pleiotropic genes that possess dual immunological and neurological topology and functionality. Strikingly, among the affected neuronal cell compartments, the synapses, both excitatory (glutamatergic) and inhibitory (GABA-ergic), emerged in this study as the most significantly transcriptionally dysregulated compartments. This was supported by the evidence of reduction in expression of the pan-presynaptic protein marker synaptophysin in the cerebellar and spinal motor circuitries. In addition, we identified a set of the pleiotropic genes with a triple topology and functionality in the immune system, nervous system in general (including functions of neuronal cell components and glial cells), and specifically, in neuron-to-neuron synapses (i.e., immune-neural-synaptic pleiotropic genes; n = 36). Focused spatial and functional interrogation of these pleiotropic genes that became possible with a recent release of a comprehensive synaptic biology knowledgebase SynGO (Koopmans et al., 2019), suggested their involvement in the regulation of synapse (including pre- and post-synapses) organization, chemical synaptic transmission, and protein catabolic processes. This builds upon our conclusion that transcriptional responses to WNV infection result in altered regulation of the organization and function of synapses. Many immune molecules in the neuron-microglia-astrocyte-T-cell axes have been implicated in synaptic pruning during development or pathological elimination/loss of synapses during disease, depending on the neuronal network (Perry and O'Connor, 2008; Eroglu and Barres, 2010; Verkhratsky and Nedergaard, 2014; Chung et al., 2013; Tröscher et al., 2019; Kreutzfeldt et al., 2013; Di Liberto et al., 2018; Stephan et al., 2012). Synaptic elimination in the hippocampus by a complement-microglial axis (Vasek et al., 2016) with assistance of T cells (Garber et al., 2019) has also been proposed as an underlying mechanism of flavivirus-associated cognitive dysfunction in a murine model. Our data in NHPs also support a significant role of the neuron-microglia-astrocyte-T-cell axes in the pathogenesis of WNV-ND at the gene and protein expression levels. Furthermore, the identification of the complement system C1QA gene (which was upregulated) as the immune-neural-synaptic pleiotropic gene in this study supports a growing evidence of the multiple roles of its protein product in the immune system and in elimination of synapses during development and disease (Perry and O'Connor, 2008). We propose a scenario in which changes in expression of pleiotropic genes that are intended to activate and maintain immune responses to virus infection would have unintended and potentially devastating off-target effects on neuron-to-neuron synapses and chemical neurotransmission.

Our findings were facilitated by a functional enrichment approach that allows detection of modest but coordinated changes in expression of groups of functionally related genes that may elude other methods (Mi et al., 2019). A similar approach has been used to detecting subtle but important coordinated changes in gene expression associated with complex human disorders that otherwise may be overlooked (Mootha et al., 2003). By using this approach, we were able to detect subtle coordinated shifts in regulation of the nervous system that were accompanied by disruption of the integrity and function of cellular compartments in infected neurons, especially their synapses. While these shifts are subtle, they may have an enormous impact, especially on the regulation of such crucial function as neural function. Consistent with this, we found that WNV infection disrupts transcriptional homeostasis of the cerebellar and spinal synapses, affecting the regulation of synapse organization, chemical synaptic transmission, metabolism, and transport. WNV-infection-induced upregulated and downregulated synaptic genes shared their topology and functionality, suggesting a bidirectional transcriptional dysregulation of synapse biology. These changes offer a mechanistic explanation for WNV-induced impairment of motor functions such as ataxia, tremor, and limb weakness/paralysis in our NHP model (Maximova et al., 2014) and in humans with WNV-ND (Lenka et al., 2019).

In support of the important role of the dysregulation of synapse biology with ensuing impairment of neurotransmission as the leading mechanisms underlying pathogenesis and neurological presentations of WNV-ND in primate hosts (NHPs and humans), we also provide evidence of strict transcriptional regulation of the cell death processes in neurons. In fact, it appears that at the advanced-symptomatic stage of WNV-ND, the magnitude of regulation of cell death processes in neurons does not exceed that in lymphocytes, suggesting that neurons are not a major cell type targeted by increased transcriptional regulation of cell death processes. Instead, it appears that the major cell types with increased regulation of cell death processes were lymphocytes (T and B cells that are expected to infiltrate the CNS parenchyma as a part of the adaptive immune response to viral infection [Maximova and Pletnev, 2018]). Thus, a closer examination of transcriptional regulation of cell death appears to support the lack of extensive activation of cell death pathways in neurons, especially when compared to other cell types in WNV-infected CNS (i.e., infiltrating T and B cells). In addition, a CNS-region-specific activation of the programmed cell death in infiltrated lymphocytes suggests the initiation of the resolution phase of the acute cellular inflammatory responses in WNV-infected CNS. This is consistent with our previously published results with the NHP model for other flavivirus infections of the CNS (Maximova et al., 2009), where we showed that the perivascular infiltrated lymphocytes rather than neurons were undergoing apoptosis.

Since WNV can spread trans-synaptically (Maximova et al., 2016), changes in synaptic homeostasis identified in this study may indicate either a direct damaging impact of infection on neurotransmission or an attempt by the host to arrest virus dissemination and compensate for changes in function. More studies will be required to understand to what extent these synaptic changes represent pathological, protective, and/or compensatory mechanisms, as well as whether they are reversible. This may guide a selection of potential therapeutic targets. For instance, our data implicate a disturbance in transcriptional regulation of glutamatergic (excitatory) and GABA-ergic (inhibitory) synapses, suggesting an imbalance between excitatory and inhibitory neurotransmission in the pathogenesis of WNV-ND. If confirmed, this excitatory-inhibitory imbalance could potentially be pharmacologically targeted with existing drugs as a symptomatic treatment (e.g., with medications increasing GABA-ergic neurotransmission, which are effective in symptomatic treatment of essential tremor, a disease affecting the Purkinje cells in the cerebellum [Louis, 2016], and in this respect, similar to WNV-ND).

Intriguingly, in our intracerebral NHP model of WNV-ND, transcriptional changes were detected earlier in the lumbar spinal cord, the CNS region most remote from the site of virus inoculation (thalamus), compared to the cerebellum. The changes in gene expression also became much stronger in the spinal cord, compared to the cerebellum, as infection progressed. At least three pathophysiological components should be considered to explain this phenomenon: (i) mode of the virus spread to the target neurons in the cerebellum (i.e., Purkinje cells) and spinal cord (i.e., spinal motor neurons) from the site of intracerebral inoculation (i.e., motor thalamus) by axonal transport (i.e., anterograde and/or retrograde); (ii) ability of the virus to establish productive replication in specific neuronal cell types (i.e., PCs and SMNs), and (iii) CNS-region specificity of host responses to infection. Our previous reconstruction of the directionality of trans-synaptic virus spread based on the neuroanatomical connectivity (Maximova et al., 2016) suggested that the virus may have reached the Purkinje cells by 7 dpi only by using the retrograde axonal transport (i.e., motor thalamus → deep cerebellar nuclei → Purkinje cells), while to reach the spinal motor neurons, the virus may have used the anterograde axonal transport (i.e., motor thalamus → corticospinal motor neurons → spinal motor neurons). Therefore, one possibility is that the differences in the mode and speed of axonal transport (Black, 2016) used by the virus to spread to target structures played a role. However, the kinetics of viral replication in the cerebellum and lumbar spinal cord were very similar (1.8 ± 0.2 and <1.7 log₁₀PFU/g [at or below the limit of detection] at 3 dpi; 6.4 ± 0.3 and 5.6 ± 0.4 log₁₀PFU/g at 7 dpi; and 6.9 ± 0.1 and 5.8 ± 0.2 log₁₀PFU/g at 9 dpi, respectively) (Maximova et al., 2016; Maximova et al., 2014), suggesting that the virus had an equal opportunity to infect respective target neurons and to use their cellular machinery to establish productive replication. Indeed, by another measure of virus production, the amount of WNV-antigens was increasing in the Purkinje cells and spinal motor neurons at similar rates (Maximova et al., 2016). Thus, the possibility that these neuronal types have a different ability to support virus replication seems unlikely. Therefore, a more likely explanation of the differences in magnitude and timing of transcriptional regulation of the immune-neural-synaptic axis in WNV-infected cerebellum and spinal cord is that the host responses to infection in these CNS structures are site-specific. CNS-site-specificity in immune responses to viral infections and in mechanisms of the non-cytolytic clearance of the virus from non-renewable neurons by T and B cells has long been recognized (Griffin, 2011; Griffin, 2003; Binder and Griffin, 2001; Cho et al., 2013). Further supporting this concept, we show different patterns of T and B cell infiltration in the cerebellum versus spinal cord, as well as a differential transcriptional regulation of apoptosis of these cells.

This study may have important translational implications associated with the use of NHPs that have a high level of genetic homology to humans, parallels in functioning of the immune system, similar to humans organization of neuroanatomical pathways and skilled motor behavior and hand dexterity, all of which are coupled with their outbred nature (Messaoudi et al., 2011; Lemon and Griffiths, 2005). However, as with any animal model of disease, some limitations need to be considered before translating the findings to humans. Our NHP model of WNV-ND, in which virus is introduced directly into the neural parenchyma (Maximova et al., 2014), may seem artificial since in humans WNV invades the CNS and causes neurological disease following peripheral infection either through mosquito bites or after transfusion/transplantation of infected blood/organs (Sejvar, 2016). However, since viremia is a major neuropathogenic determinant of the flavivirus entry into the CNS (reviewed in Maximova and Pletnev, 2018), virus invasion of the CNS most likely occurs when immune checkpoints controlling peripheral infection and ensuing viremia (Montgomery and Murray, 2015; Tobler et al., 2008; Qian et al., 2015) have failed. Therefore, our intracerebral NHP model may recapitulate the pathogenesis of WNV infection in humans that takes place after virus invasion of the CNS had occurred, by bypassing the peripheral immune control mechanisms that would otherwise limit virus neuroinvasion in the vast majority of infected individuals. Indeed, less than 1% of infected people develop neuroinvasive disease after natural mosquito-borne or iatrogenic blood-borne exposure to WNV (Sejvar, 2016). Clearly, it is not feasible to recapitulate such exposure settings in the NHPs when 100 animals would be needed to potentially induce WNV-ND in only one animal. Nonetheless, we previously showed that the intracerebral NHP model of WNV-ND closely mimics WNV neuroinvasive disease in humans in respect to the (i) gradient in severity of affected CNS structures, (ii) neuropathology, and (iii) ensuing signs of neurological impairment (Maximova et al., 2014), making this model indispensable for studies of disease neuropathogenesis and testing therapeutics and vaccine safety.

Besides an immunocompromised state, age-related alterations in immune responses are the most well-defined risk factors for increased susceptibility to severe WNV-ND (Montgomery, 2017), but how changes associated with aging may influence regulation of the immune-neural-synaptic axis in elderly with WNV-ND warrants further investigation. Peripheral inoculation of immunodeficient or aged NHPs does not result in WNV-ND even with very high doses of WNV (Wertheimer et al., 2010). Therefore, future studies may address the impact of WNV infection on the immune-neural-synaptic axis in older animals using the intracerebral route of infection.

The course of WNV-ND in our animals was abruptly interrupted at the height of neurological signs (9 dpi) due to humane animal care requirements. Therefore, there is no way to know if the disease would lead to a death or recovery, and how the phenotypic features uncovered in this study would change with time. The rapid course of WNV-ND in our NHP model is consistent with limited data showing that acute WNV encephalitis in NHPs leads to death or a moribund state that requires euthanasia 9–15 days after infection (Arroyo et al., 2004; Pogodina et al., 1983). However, a non-fatal encephalitis in NHPs, possibly associated with immunosuppression, can take a subacute course followed by virus persistence (Pogodina et al., 1983). It is important to underscore that even during the acute course of WNV-ND in this study, several developmental processes such as response to wounding, gliogenesis, and extracellular matrix/tissue remodeling became transcriptionally upregulated. However, at the height of neurological signs (i.e., advanced-symptomatic stage of WNV-ND) there was a downregulation rather than activation of neuronal developmental processes such as neuron development and axon guidance. Nevertheless, deeper analysis of the transcriptional regulation of established axon growth inhibitory or permissive molecules (Anderson et al., 2016) revealed a progressive trend which was more skewed to the axon growth permissive rather than inhibitory environment, suggesting that neuronal network repair programs may have been already initiated at the acute stage of neurological disease. Future studies should investigate correlates of recovery versus persistent neurological impairment as outcomes of WNV-ND. This may provide valuable insights and inform on how to harness potentially beneficial processes that lead to resolution of infection and neural repair.

In summary, our findings add a new dimension to understanding of regulation of the immune-neural-synaptic axis and how its homeostasis is altered during virus infection in primates. Induction of pleiotropic genes with distinct functions in each component of the immune-neural-synaptic axis suggests an unintended off-target negative impact of virus-induced immune responses on neurotransmission in the CNS. Since activation of expression of the pleiotropic genes reported here may be a part of conserved host immune responses to many viral infections, our data may serve as a resource in the search for new therapeutic approaches to restore homeostasis in interactions between the nervous and immune system at the time when virus has been cleared from the CNS.

The NHP gene expression data have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE122798 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE122798 ).

1. 1. Maximova OA
   2. Sturdevant DE
   3. Kash JC
   4. Kanakabandi K
   5. Xiao Y
   6. Minai M
   7. Moore IN
   8. Taubenberger J
   9. Martens C
   10. Cohen JI
   11. Pletnev AG

   (2021) NCBI Gene Expression Omnibus

   ID GSE122798. Virus infection of the CNS disrupts the immune-neural-synaptic axis via induction of pleiotropic gene regulation of host responses.

   https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE122798

1. 1. Anderson MA
   2. Burda JE
   3. Ren Y
   4. Ao Y
   5. O'Shea TM
   6. Kawaguchi R
   7. Coppola G
   8. Khakh BS
   9. Deming TJ
   10. Sofroniew MV

   (2016) Astrocyte scar formation aids central nervous system axon regeneration

   Nature 532:195–200.

   https://doi.org/10.1038/nature17623

   • PubMed
   • Google Scholar
2. 1. Armah HB
   2. Wang G
   3. Omalu BI
   4. Tesh RB
   5. Gyure KA
   6. Chute DJ
   7. Smith RD
   8. Dulai P
   9. Vinters HV
   10. Kleinschmidt-DeMasters BK
   11. Wiley CA

   (2007) Systemic distribution of west nile virus infection: postmortem immunohistochemical study of six cases

   Brain Pathology 17:354–362.

   https://doi.org/10.1111/j.1750-3639.2007.00080.x

   • PubMed
   • Google Scholar
3. 1. Arroyo J
   2. Miller C
   3. Catalan J
   4. Myers GA
   5. Ratterree MS
   6. Trent DW
   7. Monath TP

   (2004) ChimeriVax-West nile virus live-attenuated vaccine: preclinical evaluation of safety, immunogenicity, and efficacy

   Journal of Virology 78:12497–12507.

   https://doi.org/10.1128/JVI.78.22.12497-12507.2004

   • PubMed
   • Google Scholar
4. 1. Betsem E
   2. Kaidarova Z
   3. Stramer SL
   4. Shaz B
   5. Sayers M
   6. LeParc G
   7. Custer B
   8. Busch MP
   9. Murphy EL

   (2017) Correlation of west nile virus incidence in donated blood with west nile neuroinvasive disease rates, united states, 2010-2012

   Emerging Infectious Diseases 23:212–219.

   https://doi.org/10.3201/eid2302.161058

   • PubMed
   • Google Scholar
5. 1. Binder GK
   2. Griffin DE

   (2001) Interferon-gamma-mediated site-specific clearance of Alphavirus from CNS neurons

   Science 293:303–306.

   https://doi.org/10.1126/science.1059742

   • PubMed
   • Google Scholar
6. 1. Black MM

   (2016) Axonal transport: the orderly motion of axonal structures

   Methods in Cell Biology 131:1–19.

   https://doi.org/10.1016/bs.mcb.2015.06.001

   • PubMed
   • Google Scholar
7. 1. Boulanger LM

   (2009) Immune proteins in brain development and synaptic plasticity

   Neuron 64:93–109.

   https://doi.org/10.1016/j.neuron.2009.09.001

   • PubMed
   • Google Scholar
8. 1. Bourgeois MA
   2. Denslow ND
   3. Seino KS
   4. Barber DS
   5. Long MT

   (2011) Gene expression analysis in the thalamus and cerebrum of horses experimentally infected with west nile virus

   PLOS ONE 6:e24371.

   https://doi.org/10.1371/journal.pone.0024371

   • PubMed
   • Google Scholar
9. 1. Broom BM
   2. Ryan MC
   3. Brown RE
   4. Ikeda F
   5. Stucky M
   6. Kane DW
   7. Melott J
   8. Wakefield C
   9. Casasent TD
   10. Akbani R
   11. Weinstein JN

   (2017) A galaxy implementation of Next-Generation clustered heatmaps for interactive exploration of molecular profiling data

   Cancer Research 77:e23–e26.

   https://doi.org/10.1158/0008-5472.CAN-17-0318

   • PubMed
   • Google Scholar
10. 1. Burda JE
    2. Sofroniew MV

    (2014) Reactive gliosis and the multicellular response to CNS damage and disease

    Neuron 81:229–248.

    https://doi.org/10.1016/j.neuron.2013.12.034

    • PubMed
    • Google Scholar
11. 1. Cho H
    2. Proll SC
    3. Szretter KJ
    4. Katze MG
    5. Gale M
    6. Diamond MS

    (2013) Differential innate immune response programs in neuronal subtypes determine susceptibility to infection in the brain by positive-stranded RNA viruses

    Nature Medicine 19:458–464.

    https://doi.org/10.1038/nm.3108

    • PubMed
    • Google Scholar
12. 1. Chung WS
    2. Clarke LE
    3. Wang GX
    4. Stafford BK
    5. Sher A
    6. Chakraborty C
    7. Joung J
    8. Foo LC
    9. Thompson A
    10. Chen C
    11. Smith SJ
    12. Barres BA

    (2013) Astrocytes mediate synapse elimination through MEGF10 and MERTK pathways

    Nature 504:394–400.

    https://doi.org/10.1038/nature12776

    • PubMed
    • Google Scholar
13. 1. Clark AG
    2. Glanowski S
    3. Nielsen R
    4. Thomas PD
    5. Kejariwal A
    6. Todd MA
    7. Tanenbaum DM
    8. Civello D
    9. Lu F
    10. Murphy B
    11. Ferriera S
    12. Wang G
    13. Zheng X
    14. White TJ
    15. Sninsky JJ
    16. Adams MD
    17. Cargill M

    (2003) Inferring nonneutral evolution from human-chimp-mouse orthologous gene trios

    Science 302:1960–1963.

    https://doi.org/10.1126/science.1088821

    • PubMed
    • Google Scholar
14. 1. Clarke P
    2. Leser JS
    3. Bowen RA
    4. Tyler KL

    (2014a) Virus-induced transcriptional changes in the brain include the differential expression of genes associated with interferon, apoptosis, interleukin 17 receptor A, and glutamate signaling as well as flavivirus-specific upregulation of tRNA synthetases

    mBio 5:e00902–14.

    https://doi.org/10.1128/mBio.00902-14

    • PubMed
    • Google Scholar
15. 1. Clarke P
    2. Leser JS
    3. Quick ED
    4. Dionne KR
    5. Beckham JD
    6. Tyler KL

    (2014b) Death receptor-mediated apoptotic signaling is activated in the brain following infection with west nile virus in the absence of a peripheral immune response

    Journal of Virology 88:1080–1089.

    https://doi.org/10.1128/JVI.02944-13

    • PubMed
    • Google Scholar
16. 1. Cushing MM
    2. Brat DJ
    3. Mosunjac MI
    4. Hennigar RA
    5. Jernigan DB
    6. Lanciotti R
    7. Petersen LR
    8. Goldsmith C
    9. Rollin PE
    10. Shieh WJ
    11. Guarner J
    12. Zaki SR

    (2004) Fatal west nile virus encephalitis in a renal transplant recipient

    American Journal of Clinical Pathology 121:26–31.

    https://doi.org/10.1309/G23CP54DAR1BCY8L

    • PubMed
    • Google Scholar
17. 1. Davis LE
    2. DeBiasi R
    3. Goade DE
    4. Haaland KY
    5. Harrington JA
    6. Harnar JB
    7. Pergam SA
    8. King MK
    9. DeMasters BK
    10. Tyler KL

    (2006) West nile virus neuroinvasive disease

    Annals of Neurology 60:286–300.

    https://doi.org/10.1002/ana.20959

    • PubMed
    • Google Scholar
18. 1. Di Liberto G
    2. Pantelyushin S
    3. Kreutzfeldt M
    4. Page N
    5. Musardo S
    6. Coras R
    7. Steinbach K
    8. Vincenti I
    9. Klimek B
    10. Lingner T
    11. Salinas G
    12. Lin-Marq N
    13. Staszewski O
    14. Costa Jordão MJ
    15. Wagner I
    16. Egervari K
    17. Mack M
    18. Bellone C
    19. Blümcke I
    20. Prinz M
    21. Pinschewer DD
    22. Merkler D

    (2018) Neurons under T cell attack coordinate Phagocyte-Mediated synaptic stripping

    Cell 175:458–471.

    https://doi.org/10.1016/j.cell.2018.07.049

    • PubMed
    • Google Scholar
19. 1. Diamond MS
    2. Gale M

    (2012) Cell-intrinsic innate immune control of west nile virus infection

    Trends in Immunology 33:522–530.

    https://doi.org/10.1016/j.it.2012.05.008

    • PubMed
    • Google Scholar
20. 1. Edgar R
    2. Domrachev M
    3. Lash AE

    (2002) Gene expression omnibus: ncbi gene expression and hybridization array data repository

    Nucleic Acids Research 30:207–210.

    https://doi.org/10.1093/nar/30.1.207

    • PubMed
    • Google Scholar
21. 1. Eroglu C
    2. Barres BA

    (2010) Regulation of synaptic connectivity by Glia

    Nature 468:223–231.

    https://doi.org/10.1038/nature09612

    • PubMed
    • Google Scholar
22. 1. Fabregat A
    2. Jupe S
    3. Matthews L
    4. Sidiropoulos K
    5. Gillespie M
    6. Garapati P
    7. Haw R
    8. Jassal B
    9. Korninger F
    10. May B
    11. Milacic M
    12. Roca CD
    13. Rothfels K
    14. Sevilla C
    15. Shamovsky V
    16. Shorser S
    17. Varusai T
    18. Viteri G
    19. Weiser J
    20. Wu G
    21. Stein L
    22. Hermjakob H
    23. D'Eustachio P

    (2018) The reactome pathway knowledgebase

    Nucleic Acids Research 46:D649–D655.

    https://doi.org/10.1093/nar/gkx1132

    • PubMed
    • Google Scholar
23. 1. Ganna A
    2. Genovese G
    3. Howrigan DP
    4. Byrnes A
    5. Kurki M
    6. Zekavat SM
    7. Whelan CW
    8. Kals M
    9. Nivard MG
    10. Bloemendal A
    11. Bloom JM
    12. Goldstein JI
    13. Poterba T
    14. Seed C
    15. Handsaker RE
    16. Natarajan P
    17. Mägi R
    18. Gage D
    19. Robinson EB
    20. Metspalu A
    21. Salomaa V
    22. Suvisaari J
    23. Purcell SM
    24. Sklar P
    25. Kathiresan S
    26. Daly MJ
    27. McCarroll SA
    28. Sullivan PF
    29. Palotie A
    30. Esko T
    31. Hultman C
    32. Neale BM

    (2016) Ultra-rare disruptive and damaging mutations influence educational attainment in the general population

    Nature Neuroscience 19:1563–1565.

    https://doi.org/10.1038/nn.4404

    • PubMed
    • Google Scholar
24. 1. Garber C
    2. Soung A
    3. Vollmer LL
    4. Kanmogne M
    5. Last A
    6. Brown J
    7. Klein RS

    (2019) T cells promote microglia-mediated synaptic elimination and cognitive dysfunction during recovery from neuropathogenic flaviviruses

    Nature Neuroscience 22:1276–1288.

    https://doi.org/10.1038/s41593-019-0427-y

    • PubMed
    • Google Scholar
25. 1. Green R
    2. Wilkins C
    3. Thomas S
    4. Sekine A
    5. Ireton RC
    6. Ferris MT
    7. Hendrick DM
    8. Voss K
    9. Pardo-Manuel de Villena F
    10. Baric RS
    11. Heise MT
    12. Gale M

    (2016) Transcriptional profiles of WNV neurovirulence in a genetically diverse collaborative cross population

    Genomics Data 10:137–140.

    https://doi.org/10.1016/j.gdata.2016.10.005

    • PubMed
    • Google Scholar
26. 1. Griffin DE

    (2003) Immune responses to RNA-virus infections of the CNS

    Nature Reviews Immunology 3:493–502.

    https://doi.org/10.1038/nri1105

    • PubMed
    • Google Scholar
27. 1. Griffin DE

    (2011) Viral encephalomyelitis

    PLOS Pathogens 7:e1002004.

    https://doi.org/10.1371/journal.ppat.1002004

    • PubMed
    • Google Scholar
28. 1. Guarner J
    2. Shieh WJ
    3. Hunter S
    4. Paddock CD
    5. Morken T
    6. Campbell GL
    7. Marfin AA
    8. Zaki SR

    (2004) Clinicopathologic study and laboratory diagnosis of 23 cases with west nile virus encephalomyelitis

    Human Pathology 35:983–990.

    https://doi.org/10.1016/j.humpath.2004.04.008

    • PubMed
    • Google Scholar
29. 1. Hart J
    2. Tillman G
    3. Kraut MA
    4. Chiang HS
    5. Strain JF
    6. Li Y
    7. Agrawal AG
    8. Jester P
    9. Gnann JW
    10. Whitley RJ
    11. NIAID Collaborative Antiviral Study Group West Nile Virus 210 Protocol Team

    (2014) West nile virus neuroinvasive disease: neurological manifestations and prospective longitudinal outcomes

    BMC Infectious Diseases 14:248.

    https://doi.org/10.1186/1471-2334-14-248

    • PubMed
    • Google Scholar
30. 1. Kanehisa M
    2. Sato Y
    3. Furumichi M
    4. Morishima K
    5. Tanabe M

    (2019) New approach for understanding genome variations in KEGG

    Nucleic Acids Research 47:D590–D595.

    https://doi.org/10.1093/nar/gky962

    • PubMed
    • Google Scholar
31. 1. Kleinschmidt-DeMasters BK
    2. Beckham JD

    (2015) West nile virus encephalitis 16 years later

    Brain Pathology 25:625–633.

    https://doi.org/10.1111/bpa.12280

    • PubMed
    • Google Scholar
32. 1. Koopmans F
    2. van Nierop P
    3. Andres-Alonso M
    4. Byrnes A
    5. Cijsouw T
    6. Coba MP
    7. Cornelisse LN
    8. Farrell RJ
    9. Goldschmidt HL
    10. Howrigan DP
    11. Hussain NK
    12. Imig C
    13. de Jong APH
    14. Jung H
    15. Kohansalnodehi M
    16. Kramarz B
    17. Lipstein N
    18. Lovering RC
    19. MacGillavry H
    20. Mariano V
    21. Mi H
    22. Ninov M
    23. Osumi-Sutherland D
    24. Pielot R
    25. Smalla K-H
    26. Tang H
    27. Tashman K
    28. Toonen RFG
    29. Verpelli C
    30. Reig-Viader R
    31. Watanabe K
    32. van Weering J
    33. Achsel T
    34. Ashrafi G
    35. Asi N
    36. Brown TC
    37. De Camilli P
    38. Feuermann M
    39. Foulger RE
    40. Gaudet P
    41. Joglekar A
    42. Kanellopoulos A
    43. Malenka R
    44. Nicoll RA
    45. Pulido C
    46. de Juan-Sanz J
    47. Sheng M
    48. Südhof TC
    49. Tilgner HU
    50. Bagni C
    51. Bayés Àlex
    52. Biederer T
    53. Brose N
    54. Chua JJE
    55. Dieterich DC
    56. Gundelfinger ED
    57. Hoogenraad C
    58. Huganir RL
    59. Jahn R
    60. Kaeser PS
    61. Kim E
    62. Kreutz MR
    63. McPherson PS
    64. Neale BM
    65. O’Connor V
    66. Posthuma D
    67. Ryan TA
    68. Sala C
    69. Feng G
    70. Hyman SE
    71. Thomas PD
    72. Smit AB
    73. Verhage M

    (2019) SynGO: an Evidence-Based, Expert-Curated knowledge base for the synapse

    Neuron 103:217–234.

    https://doi.org/10.1016/j.neuron.2019.05.002

    • Google Scholar
33. 1. Kosch R
    2. Delarocque J
    3. Claus P
    4. Becker SC
    5. Jung K

    (2018) Gene expression profiles in neurological tissues during West Nile virus infection: a critical meta-analysis

    BMC Genomics 19:530.

    https://doi.org/10.1186/s12864-018-4914-4

    • PubMed
    • Google Scholar
34. 1. Kreutzfeldt M
    2. Bergthaler A
    3. Fernandez M
    4. Brück W
    5. Steinbach K
    6. Vorm M
    7. Coras R
    8. Blümcke I
    9. Bonilla WV
    10. Fleige A
    11. Forman R
    12. Müller W
    13. Becher B
    14. Misgeld T
    15. Kerschensteiner M
    16. Pinschewer DD
    17. Merkler D

    (2013) Neuroprotective intervention by interferon-γ blockade prevents CD8+ T cell-mediated dendrite and synapse loss

    Journal of Experimental Medicine 210:2087–2103.

    https://doi.org/10.1084/jem.20122143

    • PubMed
    • Google Scholar
35. 1. Kumar M
    2. Belcaid M
    3. Nerurkar VR

    (2016) Identification of host genes leading to west nile virus encephalitis in mice brain using RNA-seq analysis

    Scientific Reports 6:26350.

    https://doi.org/10.1038/srep26350

    • PubMed
    • Google Scholar
36. 1. Lemon RN
    2. Griffiths J

    (2005) Comparing the function of the corticospinal system in different species: organizational differences for motor specialization?

    Muscle & Nerve 32:261–279.

    https://doi.org/10.1002/mus.20333

    • PubMed
    • Google Scholar
37. 1. Lenka A
    2. Kamat A
    3. Mittal SO

    (2019) Spectrum of movement disorders in patients with neuroinvasive west nile virus infection

    Movement Disorders Clinical Practice 6:426–433.

    https://doi.org/10.1002/mdc3.12806

    • PubMed
    • Google Scholar
38. 1. Lim SM
    2. van den Ham HJ
    3. Oduber M
    4. Martina E
    5. Zaaraoui-Boutahar F
    6. Roose JM
    7. van IJcken WFJ
    8. Osterhaus A
    9. Andeweg AC
    10. Koraka P
    11. Martina BEE

    (2017) Transcriptomic analyses reveal differential gene expression of immune and cell death pathways in the brains of mice infected with west nile virus and Chikungunya virus

    Frontiers in Microbiology 8:1556.

    https://doi.org/10.3389/fmicb.2017.01556

    • PubMed
    • Google Scholar
39. 1. Lonsdale J
    2. Thomas J
    3. Salvatore M
    4. Phillips R
    5. Lo E
    6. Shad S
    7. Hasz R
    8. Walters G
    9. Garcia F
    10. Young N
    11. Foster B
    12. Moser M
    13. Karasik E
    14. Gillard B
    15. Ramsey K
    16. Sullivan S
    17. Bridge J
    18. Magazine H
    19. Syron J
    20. Fleming J
    21. Siminoff L
    22. Traino H
    23. Mosavel M
    24. Barker L
    25. Jewell S
    26. Rohrer D
    27. Maxim D
    28. Filkins D
    29. Harbach P
    30. Cortadillo E
    31. Berghuis B
    32. Turner L
    33. Hudson E
    34. Feenstra K
    35. Sobin L
    36. Robb J
    37. Branton P
    38. Korzeniewski G
    39. Shive C
    40. Tabor D
    41. Qi L
    42. Groch K
    43. Nampally S
    44. Buia S
    45. Zimmerman A
    46. Smith A
    47. Burges R
    48. Robinson K
    49. Valentino K
    50. Bradbury D
    51. Cosentino M
    52. Diaz-Mayoral N
    53. Kennedy M
    54. Engel T
    55. Williams P
    56. Erickson K
    57. Ardlie K
    58. Winckler W
    59. Getz G
    60. DeLuca D
    61. MacArthur D
    62. Kellis M
    63. Thomson A
    64. Young T
    65. Gelfand E
    66. Donovan M
    67. Meng Y
    68. Grant G
    69. Mash D
    70. Marcus Y
    71. Basile M
    72. Liu J
    73. Zhu J
    74. Tu Z
    75. Cox NJ
    76. Nicolae DL
    77. Gamazon ER
    78. Im HK
    79. Konkashbaev A
    80. Pritchard J
    81. Stevens M
    82. Flutre T
    83. Wen X
    84. Dermitzakis ET
    85. Lappalainen T
    86. Guigo R
    87. Monlong J
    88. Sammeth M
    89. Koller D
    90. Battle A
    91. Mostafavi S
    92. McCarthy M
    93. Rivas M
    94. Maller J
    95. Rusyn I
    96. Nobel A
    97. Wright F
    98. Shabalin A
    99. Feolo M
    100. Sharopova N
    101. Sturcke A
    102. Paschal J
    103. Anderson JM
    104. Wilder EL
    105. Derr LK
    106. Green ED
    107. Struewing JP
    108. Temple G
    109. Volpi S
    110. Boyer JT
    111. Thomson EJ
    112. Guyer MS
    113. Ng C
    114. Abdallah A
    115. Colantuoni D
    116. Insel TR
    117. Koester SE
    118. Little AR
    119. Bender PK
    120. Lehner T
    121. Yao Y
    122. Compton CC
    123. Vaught JB
    124. Sawyer S
    125. Lockhart NC
    126. Demchok J
    127. Moore HF
    128. GTEx Consortium

    (2013) The Genotype-Tissue expression (GTEx) project

    Nature Genetics 45:580–585.

    https://doi.org/10.1038/ng.2653

    • PubMed
    • Google Scholar
40. 1. Louis ED

    (2016) Essential tremor: a common disorder of purkinje neurons?

    The Neuroscientist : A Review Journal Bringing Neurobiology, Neurology and Psychiatry 22:108–118.

    https://doi.org/10.1177/1073858415590351

    • PubMed
    • Google Scholar
41. 1. Maximova OA
    2. Ward JM
    3. Asher DM
    4. St Claire M
    5. Finneyfrock BW
    6. Speicher JM
    7. Murphy BR
    8. Pletnev AG

    (2008) Comparative neuropathogenesis and neurovirulence of attenuated flaviviruses in nonhuman primates

    Journal of Virology 82:5255–5268.

    https://doi.org/10.1128/JVI.00172-08

    • PubMed
    • Google Scholar
42. 1. Maximova OA
    2. Faucette LJ
    3. Ward JM
    4. Murphy BR
    5. Pletnev AG

    (2009) Cellular inflammatory response to flaviviruses in the central nervous system of a primate host

    Journal of Histochemistry & Cytochemistry 57:973–989.

    https://doi.org/10.1369/jhc.2009.954180

    • PubMed
    • Google Scholar
43. 1. Maximova OA
    2. Speicher JM
    3. Skinner JR
    4. Murphy BR
    5. St Claire MC
    6. Ragland DR
    7. Herbert RL
    8. Pare DR
    9. Moore RM
    10. Pletnev AG

    (2014) Assurance of neuroattenuation of a live vaccine against west nile virus: a comprehensive study of neuropathogenesis after infection with chimeric WN/DEN4Δ30 vaccine in comparison to two parental viruses and a surrogate Flavivirus reference vaccine

    Vaccine 32:3187–3197.

    https://doi.org/10.1016/j.vaccine.2014.04.002

    • PubMed
    • Google Scholar
44. 1. Maximova OA
    2. Bernbaum JG
    3. Pletnev AG

    (2016) West nile virus spreads transsynaptically within the pathways of motor control: anatomical and ultrastructural mapping of neuronal virus infection in the primate central nervous system

    PLOS Neglected Tropical Diseases 10:e0004980.

    https://doi.org/10.1371/journal.pntd.0004980

    • PubMed
    • Google Scholar
45. 1. Maximova OA
    2. Pletnev AG

    (2018) Flaviviruses and the central nervous system: revisiting neuropathological concepts

    Annual Review of Virology 5:255–272.

    https://doi.org/10.1146/annurev-virology-092917-043439

    • PubMed
    • Google Scholar
46. 1. Messaoudi I
    2. Estep R
    3. Robinson B
    4. Wong SW

    (2011) Nonhuman primate models of human immunology

    Antioxidants & Redox Signaling 14:261–273.

    https://doi.org/10.1089/ars.2010.3241

    • PubMed
    • Google Scholar
47. 1. Mi H
    2. Muruganujan A
    3. Huang X
    4. Ebert D
    5. Mills C
    6. Guo X
    7. Thomas PD

    (2019) Protocol update for large-scale genome and gene function analysis with the PANTHER classification system (v.14.0)

    Nature Protocols 14:703–721.

    https://doi.org/10.1038/s41596-019-0128-8

    • PubMed
    • Google Scholar
48. 1. Montgomery RR

    (2017) Age-related alterations in immune responses to west nile virus infection

    Clinical & Experimental Immunology 187:26–34.

    https://doi.org/10.1111/cei.12863

    • PubMed
    • Google Scholar
49. 1. Montgomery RR
    2. Murray KO

    (2015) Risk factors for West Nile virus infection and disease in populations and individuals

    Expert Review of Anti-Infective Therapy 13:317–325.

    https://doi.org/10.1586/14787210.2015.1007043

    • PubMed
    • Google Scholar
50. 1. Mootha VK
    2. Lindgren CM
    3. Eriksson KF
    4. Subramanian A
    5. Sihag S
    6. Lehar J
    7. Puigserver P
    8. Carlsson E
    9. Ridderstråle M
    10. Laurila E
    11. Houstis N
    12. Daly MJ
    13. Patterson N
    14. Mesirov JP
    15. Golub TR
    16. Tamayo P
    17. Spiegelman B
    18. Lander ES
    19. Hirschhorn JN
    20. Altshuler D
    21. Groop LC

    (2003) PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes

    Nature Genetics 34:267–273.

    https://doi.org/10.1038/ng1180

    • PubMed
    • Google Scholar
51. 1. Omalu BI
    2. Shakir AA
    3. Wang G
    4. Lipkin WI
    5. Wiley CA

    (2003) Fatal fulminant pan-meningo-polioencephalitis due to west nile virus

    Brain Pathology 13:465–472.

    https://doi.org/10.1111/j.1750-3639.2003.tb00477.x

    • PubMed
    • Google Scholar
52. 1. Patel H
    2. Sander B
    3. Nelder MP

    (2015) Long-term sequelae of West Nile virus-related illness: a systematic review

    The Lancet Infectious Diseases 15:951–959.

    https://doi.org/10.1016/S1473-3099(15)00134-6

    • PubMed
    • Google Scholar
53. 1. Perry VH
    2. O'Connor V

    (2008) C1q: the perfect complement for a synaptic feast?

    Nature Reviews Neuroscience 9:807–811.

    https://doi.org/10.1038/nrn2394

    • PubMed
    • Google Scholar
54. 1. Philpott DCE
    2. Nolan MS
    3. Evert N
    4. Mayes B
    5. Hesalroad D
    6. Fonken E
    7. Murray KO

    (2019) Acute and delayed deaths after west nile virus infection, Texas, USA, 2002–2012

    Emerging Infectious Diseases 25:256–264.

    https://doi.org/10.3201/eid2502.181250

    • Google Scholar
55. 1. Pierson TC
    2. Diamond MS

    (2020) The continued threat of emerging flaviviruses

    Nature Microbiology 5:796–812.

    https://doi.org/10.1038/s41564-020-0714-0

    • PubMed
    • Google Scholar
56. 1. Pogodina VV
    2. Frolova MP
    3. Malenko GV
    4. Fokina GI
    5. Koreshkova GV
    6. Kiseleva LL
    7. Bochkova NG
    8. Ralph NM

    (1983) Study on West Nile virus persistence in monkeys

    Archives of Virology 75:71–86.

    https://doi.org/10.1007/BF01314128

    • PubMed
    • Google Scholar
57. 1. Qian F
    2. Goel G
    3. Meng H
    4. Wang X
    5. You F
    6. Devine L
    7. Raddassi K
    8. Garcia MN
    9. Murray KO
    10. Bolen CR
    11. Gaujoux R
    12. Shen-Orr SS
    13. Hafler D
    14. Fikrig E
    15. Xavier R
    16. Kleinstein SH
    17. Montgomery RR

    (2015) Systems immunology reveals markers of susceptibility to west nile virus infection

    Clinical and Vaccine Immunology 22:6–16.

    https://doi.org/10.1128/CVI.00508-14

    • PubMed
    • Google Scholar
58. 1. Radisky DC
    2. Stallings-Mann M
    3. Hirai Y
    4. Bissell MJ

    (2009) Single proteins might have dual but related functions in intracellular and extracellular microenvironments

    Nature Reviews Molecular Cell Biology 10:228–234.

    https://doi.org/10.1038/nrm2633

    • PubMed
    • Google Scholar
59. 1. Ransohoff RM
    2. Brown MA

    (2012) Innate immunity in the central nervous system

    Journal of Clinical Investigation 122:1164–1171.

    https://doi.org/10.1172/JCI58644

    • PubMed
    • Google Scholar
60. 1. Ransohoff RM
    2. Cardona AE

    (2010) The myeloid cells of the central nervous system parenchyma

    Nature 468:253–262.

    https://doi.org/10.1038/nature09615

    • PubMed
    • Google Scholar
61. 1. Raudvere U
    2. Kolberg L
    3. Kuzmin I
    4. Arak T
    5. Adler P
    6. Peterson H
    7. Vilo J

    (2019) G:profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update)

    Nucleic Acids Research 47:W191–W198.

    https://doi.org/10.1093/nar/gkz369

    • PubMed
    • Google Scholar
62. 1. Ronca SE
    2. Murray KO
    3. Nolan MS

    (2019) Cumulative incidence of west nile virus infection, continental United States, 1999-2016

    Emerging Infectious Diseases 25:325–327.

    https://doi.org/10.3201/eid2502.180765

    • PubMed
    • Google Scholar
63. 1. Sejvar JJ

    (2016) West nile virus infection

    Microbiology Spectrum 4:EI10-0021-2016.

    https://doi.org/10.1128/microbiolspec.EI10-0021-2016

    • Google Scholar
64. 1. Stephan AH
    2. Barres BA
    3. Stevens B

    (2012) The complement system: an unexpected role in synaptic pruning during development and disease

    Annual Review of Neuroscience 35:369–389.

    https://doi.org/10.1146/annurev-neuro-061010-113810

    • PubMed
    • Google Scholar
65. 1. Taylor MP
    2. Enquist LW

    (2015) Axonal spread of neuroinvasive viral infections

    Trends in Microbiology 23:283–288.

    https://doi.org/10.1016/j.tim.2015.01.002

    • PubMed
    • Google Scholar
66. 1. Tobler LH
    2. Cameron MJ
    3. Lanteri MC
    4. Prince HE
    5. Danesh A
    6. Persad D
    7. Lanciotti RS
    8. Norris PJ
    9. Kelvin DJ
    10. Busch MP

    (2008) Interferon and interferon-induced chemokine expression is associated with control of acute viremia in west nile virus-infected blood donors

    The Journal of Infectious Diseases 198:979–983.

    https://doi.org/10.1086/591466

    • PubMed
    • Google Scholar
67. 1. Tröscher AR
    2. Wimmer I
    3. Quemada-Garrido L
    4. Köck U
    5. Gessl D
    6. Verberk SGS
    7. Martin B
    8. Lassmann H
    9. Bien CG
    10. Bauer J

    (2019) Microglial nodules provide the environment for pathogenic T cells in human encephalitis

    Acta Neuropathologica 137:619–635.

    https://doi.org/10.1007/s00401-019-01958-5

    • PubMed
    • Google Scholar
68. 1. Vanichanan J
    2. Salazar L
    3. Wootton SH
    4. Aguilera E
    5. Garcia MN
    6. Murray KO
    7. Hasbun R

    (2016) Use of testing for west nile virus and other arboviruses

    Emerging Infectious Diseases 22:1587–1593.

    https://doi.org/10.3201/eid2209.152050

    • Google Scholar
69. 1. Vasek MJ
    2. Garber C
    3. Dorsey D
    4. Durrant DM
    5. Bollman B
    6. Soung A
    7. Yu J
    8. Perez-Torres C
    9. Frouin A
    10. Wilton DK
    11. Funk K
    12. DeMasters BK
    13. Jiang X
    14. Bowen JR
    15. Mennerick S
    16. Robinson JK
    17. Garbow JR
    18. Tyler KL
    19. Suthar MS
    20. Schmidt RE
    21. Stevens B
    22. Klein RS

    (2016) A complement-microglial Axis drives synapse loss during virus-induced memory impairment

    Nature 534:538–543.

    https://doi.org/10.1038/nature18283

    • PubMed
    • Google Scholar
70. 1. Venter M
    2. Myers TG
    3. Wilson MA
    4. Kindt TJ
    5. Paweska JT
    6. Burt FJ
    7. Leman PA
    8. Swanepoel R

    (2005) Gene expression in mice infected with West Nile virus strains of different neurovirulence

    Virology 342:119–140.

    https://doi.org/10.1016/j.virol.2005.07.013

    • PubMed
    • Google Scholar
71. 1. Verkhratsky A
    2. Nedergaard M

    (2014) Astroglial cradle in the life of the synapse

    Philosophical Transactions of the Royal Society B: Biological Sciences 369:20130595.

    https://doi.org/10.1098/rstb.2013.0595

    • PubMed
    • Google Scholar
72. 1. Vig PJS
    2. Lu D
    3. Paul AM
    4. Kuwar R
    5. Lopez M
    6. Stokic DS
    7. Leis AA
    8. Garrett MR
    9. Bai F

    (2018) Differential expression of genes related to innate immune responses in ex vivo spinal cord and cerebellar slice cultures infected with west nile virus

    Brain Sciences 9:1.

    https://doi.org/10.3390/brainsci9010001

    • PubMed
    • Google Scholar
73. 1. Wertheimer AM
    2. Uhrlaub JL
    3. Hirsch A
    4. Medigeshi G
    5. Sprague J
    6. Legasse A
    7. Wilk J
    8. Wiley CA
    9. Didier P
    10. Tesh RB
    11. Murray KO
    12. Axthelm MK
    13. Wong SW
    14. Nikolich-Žugich J

    (2010) Immune response to the west nile virus in aged non-human primates

    PLOS ONE 5:e15514.

    https://doi.org/10.1371/journal.pone.0015514

    • PubMed
    • Google Scholar
74. 1. Yeung MW
    2. Shing E
    3. Nelder M
    4. Sander B

    (2017) Epidemiologic and clinical parameters of west nile virus infections in humans: a scoping review

    BMC Infectious Diseases 17:609.

    https://doi.org/10.1186/s12879-017-2637-9

    • PubMed
    • Google Scholar

Author details

1. Olga A Maximova

   Laboratory of Infectious Diseases, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, United States

   Contribution

   Conceptualization, Resources, Data curation, Software, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing - original draft, Project administration, Writing - review and editing

   For correspondence

   maximovao@niaid.nih.gov

   Competing interests

   No competing interests declared

   Additional information

   Senior Author

   "This ORCID iD identifies the author of this article:" 0000-0002-8906-9738

2. Daniel E Sturdevant

   Research Technologies Branch, Genomics Unit, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Hamilton, United States

   Contribution

   Data curation, Formal analysis, Validation, Visualization, Methodology, Writing - review and editing

   Competing interests

   No competing interests declared

3. John C Kash

   Laboratory of Infectious Diseases, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, United States

   Contribution

   Data curation, Formal analysis, Visualization, Methodology

   Competing interests

   No competing interests declared

4. Kishore Kanakabandi

   Research Technologies Branch, Genomics Unit, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Hamilton, United States

   Contribution

   Data curation, Formal analysis, Validation, Investigation

   Competing interests

   No competing interests declared

5. Yongli Xiao

   Laboratory of Infectious Diseases, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, United States

   Contribution

   Data curation

   Competing interests

   No competing interests declared

6. Mahnaz Minai

   Infectious Disease Pathogenesis Section, Comparative Medicine Branch, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, United States

   Contribution

   Methodology, Performed double immunofluorescent staining

   Competing interests

   No competing interests declared

7. Ian N Moore

   Infectious Disease Pathogenesis Section, Comparative Medicine Branch, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, United States

   Contribution

   Supervision, Oversaw double immunofluorescent staining

   Competing interests

   No competing interests declared

8. Jeff Taubenberger

   Laboratory of Infectious Diseases, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, United States

   Contribution

   Resources, Supervision

   Competing interests

   No competing interests declared

9. Craig Martens

   Research Technologies Branch, Genomics Unit, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Hamilton, United States

   Contribution

   Resources, Supervision, Validation, Project administration, Writing - review and editing

   Competing interests

   No competing interests declared

10. Jeffrey I Cohen

    Laboratory of Infectious Diseases, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, United States

    Contribution

    Resources, Supervision, Funding acquisition, Writing - review and editing

    Competing interests

    No competing interests declared

11. Alexander G Pletnev

    Laboratory of Infectious Diseases, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, United States

    Contribution

    Resources, Supervision, Funding acquisition, Writing - review and editing

    Competing interests

    No competing interests declared

• 2,150

  views

• 196

  downloads

• 13

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.