Skip to main content

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

Abstract

Background

Infections with the West Nile virus (WNV) can attack neurological tissues in the host and alter gene expression levels therein. Several individual studies have analyzed these changes in the transcriptome based on measurements with DNA microarrays. Individual microarray studies produce a high-dimensional data structure with the number of studied genes exceeding the available sample size by far. Therefore, the level of scientific evidence of these studies is rather low and results can remain uncertain. Furthermore, the individual studies concentrate on different types of tissues or different time points after infection. A general statement regarding the transcriptional changes through WNV infection in neurological tissues is therefore hard to make. We screened public databases for transcriptome expression studies related to WNV infections and used different analysis pipelines to perform meta-analyses of these data with the goal of obtaining more stable results and increasing the level of evidence.

Results

We generated new lists of genes differentially expressed between WNV infected neurological tissues and control samples. A comparison with these genes to findings of a meta-analysis of immunological tissues is performed to figure out tissue-specific differences. While 5.879 genes were identified exclusively in the neurological tissues, 15 genes were found exclusively in the immunological tissues, and 44 genes were commonly detected in both tissues. Most findings of the original studies could be confirmed by the meta-analysis with a higher statistical power, but some genes and GO terms related to WNV were newly detected, too. In addition, we identified gene ontology terms related to certain infection processes, which are significantly enriched among the differentially expressed genes. In the neurological tissues, 17 gene ontology terms were found significantly different, and 2 terms in the immunological tissues.

Conclusions

A critical discussion of our findings shows benefits but also limitations of the meta-analytic approach. In summary, the produced gene lists, identified gene ontology terms and network reconstructions appear to be more reliable than the results from the individual studies. Our meta-analysis provides a basis for further research on the transcriptional mechanisms by WNV infections in neurological tissues.

Background

Epidemiology of West Nile fever

West Nile virus [(WNV); family Flaviviridae, genus Flavivirus] is one of the most important emerging virus infections in Europe. Repeated outbreaks in Europe and the dramatic spread of the virus in the United States in the past years [13] have illustrated the high potential of this virus to spread globally. Global trade and travel activities further increase the risk of WNV import to formally unaffected region for example Germany and Central Europe [4, 5]. Since its first isolation in Uganda in 1937 [6], WNV has been isolated from mosquitoes in Eurasia [1] and Australia [7, 8]. Moreover, following a single introduction in New York in 1999, WNV has spread throughout the Americas [2, 3]. In nature West Nile virus is maintained in an enzootic cycle involving ornithophilic mosquitoes and birds but can infect humans, equines and other vertebrates as illustrated by repeated cases of WNV encephalitis in horses and humans [7, 911]. Human infection of WNV usually results in a mild febrile illness including fever, headache and fatigue which last from 1 to a few days (West Nile fever) [12]. However, severe cases with cognitive dysfunction and flaccid paralysis have also been observed [13, 14]. The initial mechanism by which WNV spreads in the body and finally reaches the central nervous system (CNS) following a mosquito bite is not completely clear. It is thought that initial replication of the virus at the bite site takes place in skin epithelial cells and regional lymph nodes. Following primary viremia, WNV spreads to the reticuloendothelial system (RES; e.g kidney, spleen) [15]. This may lead to the onset of unspecific symptoms and a secondary viremia due to replication of the virus within the RES. In rare case of CNS involvement, the virus spreads from the secondary viremia into the CNS. The mechanism by which the virus crosses the blood brain barrier (BBB) is not entirely clear, although some evidence points to an involvement of tumor necrosis factor alpha (TNF- α)-mediated changes in endothelial cell permeability, which enables the virus to enter the CNS [16]. Matrix Metalloproteinase 9 (MMP9) enhances the permeability of the BBB, as well [17]. In contrary to these infection promoting genes, IFN- λ shows protective effects by strengthening the BBB-integrity [18]. Other studies with mice indicate that infection or passive transport through the endothelium or choroid plexus epithelial cells, infection of the olfactory neurons and spread to the olfactory bulb, transport by infected immune cells, or direct axonal retrograde transport from infected peripheral neurons might also allow WNV to cross the BBB [1922]. However, the exact mechanisms as well as the involvement of the immune response mounted by CNS cells towards WNV are still under heavy investigation (reviewed in Winkelmann et al. 2016 in [23]). The analysis of transcriptional responses initiated in different tissues after infection with WNV greatly aids further understanding of diseases development.

High-throughput gene expression studies on West Nile virus infection

Modern technologies for biological research such as DNA microarrays [24] or high-throughput RNA-sequencing (RNA-seq) [25] allow to simultaneously measure gene expression levels of thousands of genes in biological samples. With such methods of gene expression profiling it is possible to identify genes involved in the pathogenesis of diseases and the data can contribute to the understanding of the molecular mechanisms in cells and tissues. In the context of infectious diseases it is often of interest to compare gene expression profiles between infected and non-infected individuals, between different stages after infection or between different tissues. Gene expression studies can help to better understand the role of the transcriptome in immune response [26] or immune dysregulation [27]. In the research on WNV infections many gene expression studies have been performed to correlate expression profiles with different experimental factors such as genetic background [28], infected tissue [29], time after infection [30] or expression of specific genes in the host [31].

Meta-analysis of high-throughput gene expression studies and public databases

Individual scientific experiments or studies that are based on small sample sizes — such as many microarray or RNA-seq studies — usually have only a small statistical power and thus a limited level of scientific evidence. As a consequence, the reproducibility of study results is difficult and methods for research synthesis become more and more important [32, 33]. As one approach of research synthesis, individual studies are often combined by meta-analysis. Employing meta-analysis is widespread in the area of clinical trials and it has also become more important in the field of molecular high-throughput data. Since high-throughput data typically exhibits thousands of features (e.g. genes, mRNAs, proteins, etc.) that are observed on a much smaller number of independent biological samples, results of data analysis are even less robust than analysis of low-dimensional data [34]. Methods for meta-analysis of molecular high-throughput data were for example proposed for DNA microarray data [35, 36] or RNA-seq experiments [37]. The design of individual studies that are aggregated in meta-analyses is usually restricted to standard two-group comparisons (e.g. diseased versus control) [35]. In some rare cases, the correlation of expression data with patient survival is also considered [38]. Most of the proposed methods for meta-analysis of gene expression studies concentrate on either the combination of p-values, the combination of effects sizes (the log2 fold change in the case of gene expression data), or the direct merging of the individual expression data sets [36]. In contrast to meta-analysis of clinical trials, data merging is a possible approach for meta-analysis of transcriptome expression data since many journals ask their contributing authors to submit their gene expression data to public repositories such as ArrayExpress (AE) [39] or Gene Expression Omnibus (GEO) [40].

Meta-analysis of gene expression after West Nile virus infection

In order to identify genes that show an altered expression due to WNV infection and to obtain findings with an increased level of evidence, we performed a meta-analysis of high-throughput gene expression data currently available in public repositories. So far, no such comprehensive analysis has been published. A first approach of a meta-analysis in the context of WNV infections was conducted by Lim et al. (2017) [41], which was however only a comparison of their own transcriptome expression data with external literature findings by means of Venn diagrams. They did not merge the data or individual results by statistical methods. We are further interested in finding gene ontology (GO) terms (categorizing molecular functions, biological processes and cellular components) that are linked to the selected genes which play a role in WNV infections. Therefore, we use gene set enrichment analysis (GSEA) as a further part of our meta-analysis. From the top selected genes, we derive relevance networks that characterize the correlation between these genes. We compare the network derived in the infected samples with that from the control samples. In summary, our analysis provides a holistic approach to study the transcriptome of neurological, but also immunological tissues during WNV infection. Since meta-analysis of high-dimensional expression data is not straightforward, we implement different approaches (e.g. direct merging of the expression data from the individual studies or p-value combination methods). We take the overlap of the results of the different approaches as final result. With this we decided to run a conservative analysis that aims to reduce false positive finding by allowing much more false negative results.

In this article, we first describe the database screening and selection of appropriate data sets, as well as the bioinformatics methods and tools for meta-analysis. The results section details the outcome of different variants of meta-analysis and compares our findings with those of the individual studies. Finally, we critically discuss our findings as well as the benefits and limitations of our meta-analyses.

Methods

In this section, we describe the databases searches and selection process for comparable studies and detail the bioinformatics methods and approaches for meta-analysis.

Database search and grouping of samples

The two databases AE and GEO were searched for gene expression data using the search term 〈“west nile” OR “west nile virus” OR WNV 〉. On 10-11-2017 the GEO query with this search term returned 51 studies, while the AE query returned 36 entries of which 35 were already included in the GEO results. The identified datasets reflect a variety of experimental designs and different study questions. In order to combine comparable studies, we grouped the datasets based on several criteria, including the species from which the samples were taken, the organ system the samples belong to, the type of the experiment (gene expression profiling or other high-throughput screenings) and the sample material (e.g. full organ, isolated cells or cell lines).

Not all samples of an individual study were necessarily used to build a group. We did not consider studies which were provided by the same authors and representing the same experiment. In these cases it could not be clarified whether the same experiment was conducted twice and which version was the correct one. When ‘time after infection’ was an experimental factor, we tried to select samples that were taken at approximately the same time as in other experiments. In the same way, cell or organism line and breed were taken into account were possible. Fig. 1 summarizes the study selection process which resulted in three suitable groups (see Table 1). These groups were all related to samples from mice. No group could be formed for human samples because the identified studies appeared to be replications from the same group of authors. For other species, the number of studies was too small to find comparable experimental settings. Group 3 is the only one using isolated cells or cell cultures. As isolated cells would enhance certain characteristics compared to a more heterogeneous cell population in whole tissues, these studies can not be included in group 2 despite originating from similar tissue types. Focusing on neurological aspects, we decided to omit group 3 from further analysis because of this lack of comparability. In contrast, group 2 was used to compare the expression profiles of the two distinct organ systems and possibly identify genes that would be specific to neurological tissues. All selected studies were based on DNA microarray measurements, none of the studies was conducted by means of RNA-seq.

Fig. 1
figure 1

Flow diagram representing the dataset selection process. 87 records were found in total in GEO and AE using the search term 〈“west nile” OR “west nile virus” OR WNV 〉. Here a group refers to a set of records originating from the same organism and organ system with the same probe type. A subgroup is a set of records belonging to the same group and having similar experimental factors and comparable experiment type. Subgroups are refered to as groups in the paper since groups as defined here are only relevant in the selection process. In the end, 3 subgroups containing 9 records in total were considered for a meta-analysis. Horizontal arrows indicate exclusion relying on the criterion in the box they originate from. SuperSeries are entries of the GEO database that comprise multiple subseries. As the relevant subseries were already found individually by the search terms, the SuperSeries were removed

Table 1 Groups of studies that were built after database searches and selection processes

Throughout this manuscript we use the initials of the person who uploaded the data to AE or GEO. For three out of five studies in group 1, a publication is available. Study ‘PC’ focuses on the comparison of the expression profiles after flavivirus- (WNV and JEV) and reovirus-infections [42]. The most conspicuous pathways were distinguished. In Study ‘YV’ the recovery from WNV-infections and the ensuing neurocognitive deficits were analyzed [43]. Study ‘HH’ depicted the relevance of genes involved in immune and cell death pathways after WNV- and CHIKV-infections [41]. For study ‘H1’ and ‘H2’ no publication could be found. In group 2, only study ‘MK’ was described in a publication outlining the role of RIG-I like receptors (RLR) and type I interferons (IFN) in the restriction of WNV tissue tropism [29].

Data processing

For group 1 and 2, the raw expression data, i.e. before background correction and normalization of all samples was imported from GEO using RStudio [44] with the GEOquery package [45]. Expression data from Affymetrix chips was processed using the RMA method [46]. This method includes background correction with signal and noise close-form transformation, quantile normalization and expression level summarization using medians. We used similar processing strategies for the other types of microarrays included in the groups. Data from Agilent arrays was background-corrected with an implementation of the RMA background correction algorithm in the R-package ‘limma’ [47] and normalized using the quantile method [48]. Illumina BeadArrays data was processed using the ‘neqc’ algorithm which is very similar to RMA, except that there is no summarization step [49]. Summarization is not needed for data from Illumina and Agilent chips since they measure each transcript using a single long probe instead of many short probes targeting the same transcript as on Affymetrix chips [46, 50]. Gene transcripts represented by more than one feature are then aggregated using the median expression value. Only those genes common to all studies in the group were kept in each dataset. Samples that were not retained in the selection process because their experimental factors were not consistent with the rest of the group’s ones were removed.

Besides using the individual data sets, a merged dataset combining all individual datasets was generated for group 1 and 2. To remove batch effects in these cross-platform integrated data, the ComBat algorithm [51] was applied. For exploring grouping of samples and homogeneity between sample groups, principal component analysis (PCA) was used.

Meta-analysis

The meta-analysis was performed in two different ways. In the first variant, differentially expressed genes were detected in the individual studies. Then a p-value combination method was used to combine the individual results. Here, we used the weighted inverse normal p-value combination implemented in the ‘MetaMA’ R-package [35] to account for the varying samples sizes between the studies. This method uses one-sided p-values and therefore guarantees that no genes with contradictory fold changes in the individual studies are selected. In the second analysis variant, differentially expressed genes were selected based on the merged data set. In either variant, the linear models implemented in the R-package ‘limma’ [47] were used for differential testing. Finally, resulting p-values were adjusted to control a false-discovery rate (FDR) of 0.01 using the procedure of Benjamini and Hochberg [52]. We call the analysis variant based on the combination of the results from the individual studies the ‘late merging’ approach, and the variant based on the merged data set the ‘early merging’ approach.

As an additional analysis, GSEA was performed to identify gene sets that play a role in WNV infections. See Subramanian et al. (2005) [53] for a detailed description of the GSEA-algorithm, based on the Kolmogorov-Smirnov test. The gene sets were derived using Gene Ontology (GO) terms [54], assigning to each gene a molecular function (MF), a biological process (BP) or a cellular component (CC). Only those gene sets for which at least two genes were available in the datasets were selected. GSEA was performed on both the list of differentially expressed genes selected in the early merging approach (i.e., the merged data set) and on the list generated in the late merging approach (i.e., merged results from individual differential analyses (DA) by p-value combinations). For GSEA, another analysis variant for the meta-analysis is possible. Besides the early merging of the data and the late combination of the p-values, an ‘intermediate’ variant can be considered. In the intermediate merging approach, the individual data sets were first analyzed to identify the differentially expressed genes and the resulting lists are directly used for GSEA. Finally, the GSEA results itself are combined by the p-value combination method.

In order to study how the WNV infection changes the correlation among the top selected genes, relevance networks based on the infected and control samples were derived using the R-package ‘minet’ [55]. For the network construction, only the early merged data set was used.

Results

In this section, we describe the results of the meta-analyses performed on the studies of group 1 and group 2. First, the meta-analysis of expression data from neurological tissues is presented. Next, we describe the results of meta-analysis of immunological tissues. For both groups, we describe first the group composition and then the findings of the meta-analysis.

Meta-analysis of group 1: neurological tissues

Group composition

Table 2 lists the datasets that were used for a meta-analysis of gene expression in neurological tissues during WNV infection in mice. Each of the experiments followed a two-group design comparing samples of infected subjects with those of uninfected ones. In total, 44 samples were selected for meta-analysis. Microarray platforms from three different manufacturers (Affymetrix, Agilent and Illumina) were used in these experiments. Although we tried to identify similar studies in the selection process, some heterogeneity between the studies remained. In the case of group 1, heterogeneity was incorporated by different mice lines and different times after infection of the samples. Studies H1 and H2 were published by the same authors.

Table 2 Studies included in group 1 (neurological tissues)

Findings of the meta-analysis

The analysis comprised 15.162 genes common in the 44 selected samples from all five studies. On the other hand, between 5.135 and 9.000 genes were excluded, depending on the size of the individual studies. The PCA plot shows a clear separation between the infected and the uninfected samples in the first principal component (Fig. 2). The differential analysis on the merged dataset with the ‘limma’-procedure yielded 6.759 differentially expressed (DE) genes after FDR-correction. The weighted inverse normal p-value combination yielded 7.345 differentially expressed genes (DEG). Figure 3 shows the overlap between the DEGs that could be detected in all of the individual studies and by either merging the data (early strategy) or by combining the p-values (late strategy). The great majority of DEGs, 6.082, could be identified in all three categories. 3.304 genes were found DE in the union of the individual studies, whereas 695 DEGs were only found DE in the meta-analysis variants. Slightly less DEGs (6.759) could be identified in the data merging approach compared to the 7.345 DEGs in the p-value combination approach.

Fig. 2
figure 2

Principal component plot of the five datasets included in group 1 (neurological tissues). The axes show the first two principal components with the amount of variance these components explain from the original data. Studies are represented by different colors, WNV-infected samples are displayed as triangles, mock-infected samples as squares. The plot shows a clear separation of infected and control samples in the first principal component

Fig. 3
figure 3

Overlap of DEGs, detected by the two approaches of meta-analysis and within the individual studies. All p-values have been adjusted to control a FDR of 1%. Approaches for meta-analysis were the direct merging of the study data (early merging) and the combination of individual results by p-value combination methods (late merging). Due to an increased sample size, meta-analyses had a larger power than the individual studies yielding several thousand additional findings

For a more detailed analysis of specific genes, we ranked the early and late merging results by the value of their test statistic resulting from the ‘limma’ procedure, each. Test statistics have been used instead of p-values, because p-value combination can result in p-values of zero or bindings (p-values with the same value). Next, we extracted the top20 genes in each ranking (Table 3) and identified an overlap of eight genes in total. Among the two lists Rsad2 and Cd274 were the top ranked genes. We illustrate the log fold changes of Rsad2 and Cd274 in each individual study and in the merged data set by forest plots (Additional file 1: Figures A5.1 & A5.2). As can be seen in these plots, there is no disagreement regarding the direction of regulation for these genes; both genes are up-regulated in all studies. Furthermore, several of the confidence intervals overlap, indicating a good agreement between the studies. The complete lists of p-values and test statistics is available in Additional file 2.

Table 3 Top20 genes selected by the early merging strategy (first three columns) and late merging approach (last three columns) in group 1 (neurological tissues)

Using the eight overlapping genes found in both top20 lists, a relevance network was calculated for the control samples and for the infected samples, based on the early merged data. In this analysis, we used only the early merging pipeline, since a late merging relevance networks from the individual studies appeared not to be reasonable. The overlap of the two resulting graphs is presented in Fig. 4. A stable correlation seems to exist between the genes Rsad2, Cxcl10 and Oasl2 (black edges). Blue edges indicate gains through the infection while red edges indicate losses. The gene with the largest number of losses is Oasl2 and the gene with the largest number of gains is H2-M3. Thus, the latter one might have a more important regulatory role in infected tissues.

Fig. 4
figure 4

Relevance network with 8 genes selected in group 1 (neurological tissues). The network was built on the merged data set of all five studies. Red edges indicate correlations that go lost after WNF-infection, blue edges indicate new correlations after infection, and black edges indicate stable correlation in both states. Most infection-related losses can be observed for Oasl2 and the most gain are observed for H2-M3

Choi et al. [56] called the proportion of genes that were identified as DE in the meta-analysis but not in any of the individual studies the ‘Integration-driven Discovery Rate’ (IDR). In contrast the ‘Integration-driven Revision Rate’ (IRR) describes the percentage of genes that are declared DE in individual studies but not in meta-analysis. In the early merging analysis, these two quantities were IDR=9.1% and IRR=39.0%, respectively. For the late merging variant, the IDR was 8.8%, whereas the IRR was 33.4% (α=0.01).

GSEA was performed in three different variants as described above, and the intersection of significantly enriched GO terms from the three results was regarded. In total, 12.627 GO terms were studied by the three analysis pipelines. Additional file 1: Figure A5.3 shows the numbers of significantly enriched GO terms for each analysis variant, as well as their overlap. In general, the divergence between the pipelines was rather high: only 17 GO terms were found commonly in each analysis (α=0.01). Detailed result lists from the GSEA are provided in the Additional file 3.

Meta-analysis of group 2: immunological tissues

Group composition

Table 4 lists the studies that were used for meta-analysis of gene expression in immunological tissues during West Nile virus infection in mice. Since only Agilent chips were used in the selected experiments, 20.213 genes were in common between the two studies. 987 genes were excluded from study MK; 3.949 genes from study HT. Mice lines and the sampling times were the same, too. However, the selected samples came from tissues that, despite having immunological functions, are very different from one another.

Table 4 Studies included in group 2 (immunological tissues)

Findings of the meta-analysis

The PCA plot in Fig. 5 shows a separation between the infected and the uninfected tissues along the second principal component. Still, the samples from each organ build separate groups, which probably results from the different nature of the tissues. The differential expression analysis of the early merged dataset yielded 60 DEG after FDR-correction of the p-values (Fig. 6). In contrast, the weighted inverse normal p-value combination gives 1.666 DEGs. The overlap of both methods are 60 DEGs, as well. 683 DEGs could be detected in the union of the individual studies. Therefore, we distinguished IDRs of 18.3% for the early data merging approach and 57.1% for the late p-value combination approach. The IRRs were 96.5% and 48.9%. The top20 genes for the early and the late merging are listed in Table 5, five genes were commonly found in both lists. Among these five genes, Oas1a was the top-ranked gene in both lists. A forest plot for Oas1a is given in Additional file 1: Figure A5.4. Again, there is no disagreement regarding the direction of regulation. The complete lists of p-values and test statistics is available in Additional file 4.

Fig. 5
figure 5

Principal component plot of the three datasets included in group 2 (immunological tissues). The axes show the first two principal components with the amount of variance these components explain from the original data. Studies are represented by different colors, WNV-infected samples are displayed as triangles, mock-infected samples as squares. Although there is a strong variance between the studies on the first principal component, there is also a separation of infected and control samples on the second principal component

Fig. 6
figure 6

Overlap of DEGs, detected by the two approaches of meta-analysis and within the individual studies. All p-values have been adjusted to control a FDR of 1%. Approaches for meta-analysis were the direct merging of the study data (early merging) and the combination of individual results by p-value combination methods (late merging)

Table 5 Top20 genes selected by the early merging strategy (first two columns) and late merging approach (last two columns) in group 2 (immunological tissues)

For the GSEA, we utilized 10.337 GO terms shared in both studies. The three pipelines for GSEA meta-analysis resulted in an overlap of two significantly enriched GO terms, using an α-error of 0.01 (Additional file 1: Figure A5.5). Detailed result lists from the GSEA are provided in the Additional file 5.

Comparison between neurological and immunological tissues

Results of the differential analysis of group 1 and group 2 were compared in Fig. 7. The Venn diagram also shows the overlap of the early and late merging pipeline. Again, we regard the overlap of both pipelines as a conservative results, that gives true positive findings a higher weight by allowing for more false negatives. In the end, 44 genes could be identified as DE in the neurological and immunological tissues. On the other hand, 5.879 DEGs were found specifically in the neurological tissues; 15 DEGs were related only to the immunological tissues. Interestingly, the number of changes in the transcriptome is much larger in the neurological tissues than in the immunological ones, however this could be related to different degrees of homogeneity between the selected studies.

Fig. 7
figure 7

Comparison of the results of neurological and immunological tissues. Each group is further separated by early and late merging. The numbers show the amount of DEGs of each analysis variant, as well as their overlaps

Discussion

In this section, we discuss our findings in the neurological and immunological data and make a short comparison of both. Furthermore, we critically discuss the methodical issues of the meta-analysis.

Findings in neurological tissues

In the group of neurological tissues (group 1), eight overlapping genes were found between the top20 DEGs of the early and late analysis variant were eight DEGs: Rsad2, Cxcl10, Cd274, Oasl2, H2-M3, Trim25, Trim21 and Samd9l (Table 3). The majority of these genes play a role in the innate immune system. Rsad2, Cxcl10, Samd9l and Oasl2 are involved in the type I interferon (IFN) cell signaling pathway and are therefore part of the antiviral immune system [5759]. Rsad2 and Cxcl10 can be induced directly by IRF3 (Interferone regulatory factor) or on detours by IRF5 [57]. Rsad2, also known as Viperin, has been well described for WNV by Szretter et al. [60]. The chemokine Cxcl10 is responsible for the recruitment of CD8+ T-cells after WNV-infection [61]. CD274 (PD-L1) is a costimulatory molecule, involved in the T-cell response [62]. Trim21 and Trim25 are part of the tripartite motif family with ability to mediate ubiquitylation events. They are also induced by interferons. It has been suggested, that Oasl2 has the same effects as the Oasl in humans [58]. Oasl, as well as Trim25 are responsible for the activation, more precisely the ubiquitylation of RIG-I, which is a pattern recognition receptor, i.a. for flavivruses [63]. H2-M3 is a molecule of the MHC class Ib, known for the presentation of CD8+ T-cells. Similar biological functions could have been found for significant genes beyond the top20. We could not find a reference for the sterile alpha motif family protein 91 (Samd9l) in the context of WNV-infections. However, the SAMD9 gene is a downstream target of interferon gamma and the protein is involved in innate immunity and decreases replication of the RNA-virus Japanese Encephalitis virus [59, 64, 65]. The important role of interferons in the context of WNV infections was also described in two of the individual studies [41, 42]. The top20 genes comprise additional interferon-dependent targets, e.g. the Interferon Stimulated Exonuclease Gene 20 (Isg20). The ISG20 protein targets single-stranded RNA and displays antiviral activity towards RNA viruses [66].

The findings in the top20 DEGs match well with the results of the GSEA. Each of the 17 GO terms found in by the three analysis variants (early, intermediate or late analysis) involves a large number of genes ranging between 137 genes for the smallest GO term and 3308 for the largest. Most of these GO terms are in accordance with the results of the differential expression analysis. The GO terms GO:0048583 (regulation of response to stimulus), GO:0010033 (response to organic substance) and GO:0002684 (positive regulation of immune system process) fit to the context of WNV-infection, but are rather unspecific to allow a concrete interpretation. The same is true for the enrichment of GO term GO:0002822 (regulation of adaptive immune response based on somatic recombination of immune receptors built from immunoglobulin superfamily domains).

Findings in immunological tissues

The findings of the meta-analysis are highly compatible with the conclusions of the study from the MK dataset [29].

While group 2 (immunological tissues) had more common DEGs among the individual datasets compared to the findings in group 1, less DEG were found both for the late merging of individual study results and in the early merging analysis strategy. This result hints at a higher heterogeneity among the datasets. However, in contrast to group 1, the IDR was higher than the IRR. Thus, one could argue that the meta-analysis had higher utility in group 2. Regarding the biological interpretation of the DA-results, a clear pattern like in group 1 could not be identified. Many of the DEGs detected in group 2 are affected by IFNs and/or have antiviral functions. Several genes from the OAS- or OASL-family (oligoadenylate synthetase), as well as genes coding for Ifit proteins (Interferon Induced proteins with Tetratricopeptide repeats) can be found [67, 68]. Again, the top20 DEGs were determined for each analysis method (early and late merging). These top20 lists overlapped in the following five genes: Oas1a, Irf7, Oas1f, Oas3 and Oasl2. Oasl2 has also been found as DEG in the neurological tissues and is described above. Except Irf7, all genes are involved in the 2’-5’-oligoadenylate synthetases. It activates the endoribonuclease RNase L, which is responsible for the degradation of viral and cellular RNAs [69]. It has been shown, that Irf7 plays a major role in the regulation of Interferone response [57].

In comparison to group 1, the results of the GSEA differ more, where the finding of the majority of detected GO terms fit to the infection context. GO:0044194 (cytolytic granule) and GO:0001730 (2’-5’-oligoadenylate synthetase activity) were significantly enriched in all three pipelines (early, intermediate, late). The importance of the 2’-5’-oligoadenylate was already depicted. Nevertheless, the overall interpretation of the GSEA-result is limited by this low number of detected gene sets in the first place.

Comparison between neurological and immunological tissues

Regarding the top20 lists (Tables 3 and 5) Oasl2 is the only gene that is differentially expressed in the neurological tissues and in the immunological tissues. Thus, it might be that Oasl2 is a general WNV-affected gene, while the other detected DEGs are rather tissue-specific. Moreover, 43 genes were detected as DE in both tissues, beyond the top20. The much greater amount of 5.879 DEGs in neurological tissues, compared to the 15 DEGs in immunological tissues, enhances the assumption that WNV-infection leads to alterations in the brain.

Methodical and practical issues

Meta-analysis of high-throughput expression data has the clear advantage of an increased sample size and thus an increased power to detect differentially expressed genes and enriched gene sets such as GO terms. In addition, it can be used to clarify contradictory findings regarding the direction of expression fold changes, i.e. up-regulation of a gene in one study and down-regulation in another study. Genes that are regulated in different directions in the individual studies are not selected by the early data merging approach and not by the weighted p-value combination approach.

Besides the benefits a meta-analysis has, it also bears some risks and has some limitations. Like in any meta-analysis it is crucial to select data sets which fit to each other regarding the study question. In the case of gene expression data, studies must focus on approximately the same type of biological samples (i.e. type of tissue), and be based on the same species. This demand may lead to a low number of final studies to be involved in the meta-analysis, even if there is a plethora of studies available before the selection process. Nevertheless, selecting fewer but more suitable studies for a meta-analysis seems to be the better choice for meaningful results.

Since the different approaches for meta-analysis (data merging or result merging) provide clearly different results, the researchers must still be careful with the final biological interpretation. Here, we have chosen a conservative approach and made only a final interpretation about the overlap of findings from different analysis pipelines. This bears, however, the risk of making false negative decisions. From our point of view, false negative conclusion can be acceptable when the goal is to obtain a stable list of top differentially expressed genes. False negatives also occur by each merging step, since genes which do not occur in all studies are omitted from analysis.

Conclusions

Considering the benefits and limitations of meta-analyses, we think that our results provide a contribution to the knowledge about gene expression in neurological tissues after WNV infection that has a higher level of evidence than the individual studies have. In particular, the comparison with immunological tissues shows which genes may play a role after WNV infection in general, and which genes have a tissue-specific regulation.

Abbreviations

AE:

ArrayExpress

BBB:

Blood brain barrier

BP:

Biological process

CC:

Cellular component

CHIKV:

Chikungunya virus

CNS:

Central nervous system

DA:

Differential analysis

DE:

Differentially expressed

DEG:

Differentially expressed gene

FDR:

False-discovery rate

GEO:

Gene expression omnibus

GSEA:

Gene set enrichment analysis

GO:

Gene Ontology

IDR:

Integration-driven discovery rate

IFN:

Interferon

IRR:

Integration-driven revision rate

JEV:

Japanese encephalitis virus

MF:

Molecular function

MHC:

Major histocompatibility complex

PCA:

Principal component analysis

RES:

Reticuloendothelial system

RMA:

Robust multi-array average

WNV:

West Nile virus

References

  1. Calistri P. Epidemiology of west nile in europe and in the mediterranean basin. Open Virol J. 2010; 4(1):29–37. https://0-doi-org.brum.beds.ac.uk/10.2174/1874357901004010029.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Nash D, Mostashari F, Fine A, Miller J, O’Leary D, Murray K, Huang A, Rosenberg A, Greenberg A, Sherman M, Wong S, Layton M, West Nile Outbreak Response Working G. The outbreak of west nile virus infection in the new york city area in 1999. N Engl J Med. 2001; 344(24):1807–14. https://0-doi-org.brum.beds.ac.uk/10.1056/NEJM200106143442401.

    Article  PubMed  CAS  Google Scholar 

  3. Petersen LR, Hayes EB. West nile virus in the americas. Med Clin N Am. 2008; 92(6):1307–22. https://0-doi-org.brum.beds.ac.uk/10.1016/j.mcna.2008.07.004.

    Article  PubMed  Google Scholar 

  4. Gabriel M, Emmerich P, Frank C, Fiedler M, Rashidi-Alavijeh J, Jochum C, Gunther S, Auerhammer K, Rupprecht HJ, Blank RT, Sacher N, Pertzborn L, Stark K, Schrauzer T, Schmidt-Chanasit J. Increase in west nile virus infections imported to germany in 2012. J Clin Virol. 2013; 58(3):587–9. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jcv.2013.08.027.

    Article  PubMed  Google Scholar 

  5. Schultze-Amberger J, Emmerich P, Gunther S, Schmidt-Chanasit J. West nile virus meningoencephalitis imported into germany. Emerg Infect Dis. 2012; 18(10):1698–700. https://0-doi-org.brum.beds.ac.uk/10.3201/eid1810.120204.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Smithburn KC, Hughes TP, Paul JH, Burke AW. A neurotropic virus isolated from the blood of a native of uganda 1. Am J Trop Med Hyg. 1940; s1-20(4):471–92. https://0-doi-org.brum.beds.ac.uk/10.4269/ajtmh.1940.s1-20.471.

    Article  Google Scholar 

  7. Frost MJ, Zhang J, Edmonds JH, Prow NA, Gu X, Davis R, Hornitzky C, Arzey KE, Finlaison D, Hick P, Read A, Hobson-Peters J, May FJ, Doggett SL, Haniotis J, Russell RC, Hall RA, Khromykh AA, Kirkland PD. Characterization of virulent west nile virus kunjin strain, australia, 2011. Emerg Infect Dis. 2012; 18(5):792–800. https://0-doi-org.brum.beds.ac.uk/10.3201/eid1805.111720.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Russell RC, Dwyer DE. Arboviruses associated with human diseasein australia. Microbes Infect. 2000; 2(14):1693–704. https://0-doi-org.brum.beds.ac.uk/10.1016/s1286-4579(00)01324-1.

    Article  PubMed  CAS  Google Scholar 

  9. Tsai TF, Popovici F, Cernescu C, Campbell GL, Nedelcu NI. West nile encephalitis epidemic in southeastern romania. Lancet. 1998; 352(9130):767–71. https://0-doi-org.brum.beds.ac.uk/10.1016/s0140-6736(98)03538-7.

    Article  PubMed  CAS  Google Scholar 

  10. Garcia-Bocanegra I, Jaen-Tellez JA, Napp S, Arenas-Montes A, Fernandez-Morente M, Fernandez-Molera V, Arenas A. West nile fever outbreak in horses and humans, spain, 2010. Emerg Infect Dis. 2011; 17(12):2397–9. https://0-doi-org.brum.beds.ac.uk/10.3201/eid1712.110651.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Platonov AE, Shipulin GA, Shipulina OY, Tyutyunnik EN, Frolochkina TI, Lanciotti RS, Yazyshina S, Platonova OV, Obukhov IL, Zhukov AN, Vengerov YY, Pokrovskii VI. Outbreak of west nile virus infection, volgograd region, russia, 1999. Emerg Infect Dis. 2001; 7(1):128–32. https://0-doi-org.brum.beds.ac.uk/10.3201/eid0701.700128.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Campbell GL, Marfin AA, Lanciotti RS, Gubler DJ. West nile virus. Lancet Infect Dis. 2002; 2(9):519–29. https://0-doi-org.brum.beds.ac.uk/10.1016/s1473-3099(02)00368-7.

    Article  PubMed  Google Scholar 

  13. Ceausu E, Erscoiu S, Calistru P, Ispas D, Dorobat O, Homos M, Barbulescu C, Cojocaru I, Simion CV, Cristea C, Oprea C, Dumitrescu C, Duiculescu D, Marcu I, Mociornita C, Stoicev T, Zolotusca I, Calomfirescu C, Rusu R, Hodrea R, Geamai S, Paun L. Clinical manifestations in the west nile virus outbreak. Rom J Virol. 1997; 48(1-4):3–11.

    PubMed  CAS  Google Scholar 

  14. Sejvar JJ, Haddad MB, Tierney BC, Campbell GL, Marfin AA, Van Gerpen JA, Fleischauer A, Leis AA, Stokic DS, Petersen LR. Neurologic manifestations and outcome of west nile virus infection. JAMA. 2003; 290(4):511–5. https://0-doi-org.brum.beds.ac.uk/10.1001/jama.290.4.511.

    Article  PubMed  Google Scholar 

  15. Deubel V, Fiette L, Gounon P, Drouet MT, Khun H, Huerre M, Banet C, Malkinson M, Despres P. Variations in biological features of west nile viruses. Ann N Y Acad Sci. 2001; 951:195–206.

    Article  PubMed  CAS  Google Scholar 

  16. Wang T, Town T, Alexopoulou L, Anderson JF, Fikrig E, Flavell RA. Toll-like receptor 3 mediates west nile virus entry into the brain causing lethal encephalitis. Nat Med. 2004; 10(12):1366–73. https://0-doi-org.brum.beds.ac.uk/10.1038/nm1140.

    Article  PubMed  CAS  Google Scholar 

  17. Wang P, Dai J, Bai F, Kong KF, Wong SJ, Montgomery RR, Madri JA, Fikrig E. Matrix metalloproteinase 9 facilitates west nile virus entry into the brain. J Virol. 2008; 82(18):8978–85. https://0-doi-org.brum.beds.ac.uk/10.1128/JVI.00314-08.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Lazear HM, Daniels BP, Pinto AK, Huang AC, Vick SC, Doyle SE, Gale JM, Klein RS, Diamond MS. Interferon-lambda restricts west nile virus neuroinvasion by tightening the blood-brain barrier. Sci Transl Med. 2015; 7(284):284–59. https://0-doi-org.brum.beds.ac.uk/10.1126/scitranslmed.aaa4304.

    Article  CAS  Google Scholar 

  19. Kramer-Hammerle S, Rothenaigner I, Wolff H, Bell JE, Brack-Werner R. Cells of the central nervous system as targets and reservoirs of the human immunodeficiency virus. Virus Res. 2005; 111(2):194–213. https://0-doi-org.brum.beds.ac.uk/10.1016/j.virusres.2005.04.009.

    Article  PubMed  CAS  Google Scholar 

  20. Monath TP, Cropp CB, Harrison AK. Mode of entry of a neurotropic arbovirus into the central nervous-system - reinvestigation of an old controversy. Lab Investig. 1983; 48(4):399–410.

    PubMed  CAS  Google Scholar 

  21. Garcia-Tapia D, Loiacono CM, Kleiboeker SB. Replication of west nile virus in equine peripheral blood mononuclear cells. Vet Immunol Immunopathol. 2006; 110(3-4):229–44. https://doi.org/10.1016/j.vetimm.2005.10.003.

    Article  PubMed  CAS  Google Scholar 

  22. Hunsperger EA, Roehrig JT. Temporal analyses of the neuropathogenesis of a west nile virus infection in mice. J Neurovirol. 2006; 12(2):129–39. https://0-doi-org.brum.beds.ac.uk/10.1080/13550280600758341.

    Article  PubMed  CAS  Google Scholar 

  23. Winkelmann ER, Luo H, Wang T. West nile virus infection in the central nervous system. F1000Res. 2016; 5. https://0-doi-org.brum.beds.ac.uk/10.12688/f1000research.7404.1.

  24. Heller MJ. DNA microarray technology: Devices, systems, and applications. Annu Rev Biomed Eng. 2002; 4(1):129–53. https://0-doi-org.brum.beds.ac.uk/10.1146/annurev.bioeng.4.020702.153438.

    Article  PubMed  CAS  Google Scholar 

  25. Wang Z, Gerstein M, Snyder M. RNA-seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009; 10(1):57–63. https://0-doi-org.brum.beds.ac.uk/10.1038/nrg2484.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Mejias A, Dimo B, Suarez NM, Garcia C, Suarez-Arrabal MC, Jartti T, Blankenship D, Jordan-Villegas A, Ardura MI, Xu Z, Banchereau J, Chaussabel D, Ramilo O. Whole blood gene expression profiles to assess pathogenesis and disease severity in infants with respiratory syncytial virus infection. PLoS Med. 2013; 10(11):1001549. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pmed.1001549.

    Article  Google Scholar 

  27. Cameron CM, Cameron MJ, Bermejo-Martin JF, Ran L, Xu L, Turner PV, Ran R, Danesh A, Fang Y, Chan P-KM, Mytle N, Sullivan TJ, Collins TL, Johnson MG, Medina JC, Rowe T, Kelvin DJ. Gene expression analysis of host innate immune responses during lethal h5n1 infection in ferrets. J Virol. 2008; 82(22):11308–17. https://0-doi-org.brum.beds.ac.uk/10.1128/jvi.00691-08.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Green R, Wilkins C, Thomas S, Sekine A, Hendrick DM, Voss K, Ireton RC, Mooney M, Go JT, Choonoo G, Jeng S, de Villena FP-M, Ferris MT, McWeeney S, Gale M. Oas1b-dependent immune transcriptional profiles of west nile virus infection in the collaborative cross. G3 (Bethesda). 2017; 7(6):1665–82. https://0-doi-org.brum.beds.ac.uk/10.1534/g3.117.041624. GGG_041624[PII].

    Google Scholar 

  29. Suthar MS, Brassil MM, Blahnik G, McMillan A, Ramos HJ, Proll SC, Belisle SE, Katze MG, Gale JM. A systems biology approach reveals that tissue tropism to west nile virus is regulated by antiviral genes and innate immune cellular processes. PLoS Pathog. 2013; 9(2):1003168. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.ppat.1003168.

    Article  CAS  Google Scholar 

  30. Graham JB, Swarts JL, Wilkins C, Thomas S, Green R, Sekine A, Voss KM, Ireton RC, Mooney M, Choonoo G, Miller DR, Treuting PM, de Villena FPM, Ferris MT, McWeeney S, Gale M, Lund JM. A mouse model of chronic west nile virus disease. PLoS Pathog. 2016; 12(11):1005996. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.ppat.1005996.

    Article  CAS  Google Scholar 

  31. Zhao J, Vijay R, Zhao J, Gale M, Diamond MS, Perlman S. MAVS expressed by hematopoietic cells is critical for control of west nile virus infection and pathogenesis. J Virol. 2016; 90(16):7098–108. https://0-doi-org.brum.beds.ac.uk/10.1128/jvi.00707-16.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Ioannidis JP, Allison DB, Ball CA, Coulibaly I, Cui X, Culhane AC, Falchi M, Furlanello C, Game L, Jurman G, et al. Repeatability of published microarray gene expression analyses. Nat Genet. 2009; 41(2):149–55.

    Article  PubMed  CAS  Google Scholar 

  33. Altman DG. Some reflections on the evolution of meta-analysis. Res Synth Methods. 2015; 6(3):265–7.

    Article  PubMed  Google Scholar 

  34. Ramasamy A, Mondry A, Holmes CC, Altman DG. Key issues in conducting a meta-analysis of gene expression microarray datasets. PLoS Med. 2008; 5(9):184. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pmed.0050184.

    Article  CAS  Google Scholar 

  35. Marot G, Foulley J-L, Mayer C-D, Jaffrezic F. Moderated effect size and p-value combinations for microarray meta-analyses. Bioinformatics. 2009; 25(20):2692–9. https://0-doi-org.brum.beds.ac.uk/10.1093/bioinformatics/btp444.

    Article  PubMed  CAS  Google Scholar 

  36. Tseng GC, Ghosh D, Feingold E. Comprehensive literature review and statistical considerations for microarray meta-analysis. Nucleic Acids Res. 2012; 40(9):3785–99. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gkr1265.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Rau A, Marot G, Jaffrézic F. Differential meta-analysis of RNA-seq data from multiple studies. BMC Bioinformatics. 2014; 15(1):91. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2105-15-91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Leha A, Beißbarth T, Jung K. Sequential interim analyses of survival data in DNA microarray experiments. BMC Bioinformatics. 2011; 12(1):127. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2105-12-127.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Parkinson H. ArrayExpress–a public repository for microarray gene expression data at the EBI. Nucleic Acids Res. 2004; 33(Database issue):553–5. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gki056.

    Article  Google Scholar 

  40. Edgar R. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002; 30(1):207–10. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/30.1.207.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Lim SM, van den Ham HJ, Oduber M, Martina E, Zaaraoui-Boutahar F, Roose JM, van IWFJ, Osterhaus A, Andeweg AC, Koraka P, Martina BEE. 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. Front Microbiol. 2017; 8:1556. https://0-doi-org.brum.beds.ac.uk/10.3389/fmicb.2017.01556.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Clarke P, Leser JS, Bowen RA, Tyler KL. 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. 2014; 5(2):00902–14. https://0-doi-org.brum.beds.ac.uk/10.1128/mBio.00902-14.

    Article  CAS  Google Scholar 

  43. Vasek MJ, Garber C, Dorsey D, Durrant DM, Bollman B, Soung A, Yu J, Perez-Torres C, Frouin A, Wilton DK, Funk K, DeMasters BK, Jiang X, Bowen JR, Mennerick S, Robinson JK, Garbow JR, Tyler KL, Suthar MS, Schmidt RE, Stevens B, Klein RS. A complement–microglial axis drives synapse loss during virus-induced memory impairment. Nature. 2016; 534(7608):538–43. https://0-doi-org.brum.beds.ac.uk/10.1038/nature18283.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. RStudio Team. RStudio: Integrated Development Environment for R. RStudio, Inc.Boston: RStudio, Inc.; 2016. http://www.rstudio.com/.

    Google Scholar 

  45. Davis S, Meltzer PS. GEOquery: a bridge between the gene expression omnibus (GEO) and BioConductor. Bioinformatics. 2007; 23(14):1846–7. https://0-doi-org.brum.beds.ac.uk/10.1093/bioinformatics/btm254.

    Article  PubMed  CAS  Google Scholar 

  46. Irizarry RA. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003; 4(2):249–64. https://0-doi-org.brum.beds.ac.uk/10.1093/biostatistics/4.2.249.

    Article  PubMed  Google Scholar 

  47. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for rna-sequencing and microarray studies. Nucleic Acids Res. 2015; 43(7):47. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gkv007.

    Article  CAS  Google Scholar 

  48. Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003; 19(2):185–93. https://0-doi-org.brum.beds.ac.uk/10.1093/bioinformatics/19.2.185.

    Article  PubMed  CAS  Google Scholar 

  49. Shi W, Oshlack A, Smyth GK. Optimizing the noise versus bias trade-off for illumina whole genome expression BeadChips. Nucleic Acids Res. 2010; 38(22):204. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gkq871.

    Article  CAS  Google Scholar 

  50. Dalma-Weiszhausz DD, Warrington J, Tanimoto EY, Miyada CG. [1] the affymetrix GeneChip®; platform: An overview. In: Methods in Enzymology, Vol 410. Elsevier: 2006. p. 3–28. https://0-doi-org.brum.beds.ac.uk/10.1016/s0076-6879(06)10001-4.

  51. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical bayes methods. Biostatistics. 2006; 8(1):118–27. https://0-doi-org.brum.beds.ac.uk/10.1093/biostatistics/kxj037.

    Article  PubMed  Google Scholar 

  52. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995:289–300.

  53. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005; 102. https://0-doi-org.brum.beds.ac.uk/10.1073/pnas.0506580102.

  54. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G. Gene ontology: tool for the unification of biology. the gene ontology consortium. Nat Genet. 2000; 25(1):25–9. https://0-doi-org.brum.beds.ac.uk/10.1038/75556.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. Meyer PE, Lafitte F, Bontempi G. minet: Ar/bioconductor package for inferring large transcriptional networks using mutual information. BMC Bioinformatics. 2008; 9(1):461.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Choi JK, Yu U, Kim S, Yoo OJ. Combining multiple microarray studies and modeling interstudy variation. Bioinformatics. 2003; 19 Suppl 1:84–90.

    Article  Google Scholar 

  57. Lazear HM, Lancaster A, Wilkins C, Suthar MS, Huang A, Vick SC, Clepper L, Thackray L, Brassil MM, Virgin HW, Nikolich-Zugich J, Moses AV, Gale JM, Fruh K, Diamond MS. Irf-3, irf-5, and irf-7 coordinately regulate the type i ifn response in myeloid dendritic cells downstream of mavs signaling. PLoS Pathog. 2013; 9(1):1003118. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.ppat.1003118.

    Article  CAS  Google Scholar 

  58. Zhu J, Zhang Y, Ghosh A, Cuevas RA, Forero A, Dhar J, Ibsen MS, Schmid-Burgk JL, Schmidt T, Ganapathiraju MK, Fujita T, Hartmann R, Barik S, Hornung V, Coyne CB, Sarkar SN. Antiviral activity of human oasl protein is mediated by enhancing signaling of the rig-i rna sensor. Immunity. 2014; 40(6):936–48. https://0-doi-org.brum.beds.ac.uk/10.1016/j.immuni.2014.05.007.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  59. Lemos de Matos A, Liu J, McFadden G, Esteves PJ. Evolution and divergence of the mammalian samd9/samd9l gene family. BMC Evol Biol. 2013; 13:121. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2148-13-121.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  60. Szretter KJ, Brien JD, Thackray LB, Virgin HW, Cresswell P, Diamond MS. The interferon-inducible gene viperin restricts west nile virus pathogenesis. J Virol. 2011; 85(22):11557–66.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  61. Klein RS, Lin E, Zhang B, Luster AD, Tollett J, Samuel MA, Engle M, Diamond MS. Neuronal cxcl10 directs cd8+ t-cell recruitment and control of west nile virus encephalitis. J Virol. 2005; 79(17):11457–66. https://0-doi-org.brum.beds.ac.uk/10.1128/JVI.79.17.11457-11466.2005.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  62. Zhu Q, Egelston C, Gagnon S, Sui Y, Belyakov IM, Klinman DM, Berzofsky JA. Using 3 tlr ligands as a combination adjuvant induces qualitative changes in t cell responses needed for antiviral protection in mice. J Clin Invest. 2010; 120(2):607–16. https://0-doi-org.brum.beds.ac.uk/10.1172/JCI39293.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  63. Nazmi A, Dutta K, Basu A. Rig-i mediates innate immune response in mouse neurons following japanese encephalitis virus infection. PLoS One. 2011; 6(6):21761. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0021761.

    Article  CAS  Google Scholar 

  64. Hershkovitz D, Gross Y, Nahum S, Yehezkel S, Sarig O, Uitto J, Sprecher E. Functional characterization of samd9, a protein deficient in normophosphatemic familial tumoral calcinosis. J Invest Dermatol. 2011; 131(3):662–9. https://0-doi-org.brum.beds.ac.uk/10.1038/jid.2010.387.

    Article  PubMed  CAS  Google Scholar 

  65. Zhang LK, Chai F, Li HY, Xiao G, Guo L. Identification of host proteins involved in japanese encephalitis virus infection by quantitative proteomics analysis. J Proteome Res. 2013; 12(6):2666–78. https://0-doi-org.brum.beds.ac.uk/10.1021/pr400011k.

    Article  PubMed  CAS  Google Scholar 

  66. Zheng Z, Wang L, Pan J. Interferon-stimulated gene 20-kda protein (isg20) in infection and disease: Review and outlook. Intractable Rare Dis Res. 2017; 6(1):35–40. https://0-doi-org.brum.beds.ac.uk/10.5582/irdr.2017.01004.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Melchjorsen J, Kristiansen H, Christiansen R, Rintahaka J, Matikainen S, Paludan SR, Hartmann R. Differential regulation of the oasl and oas1 genes in response to viral infections. J Interf Cytokine Res. 2009; 29(4):199–207. https://0-doi-org.brum.beds.ac.uk/10.1089/jir.2008.0050.

    Article  CAS  Google Scholar 

  68. Diamond MS, Farzan M. The broad-spectrum antiviral functions of ifit and ifitm proteins. Nat Rev Immunol. 2013; 13(1):46–57.

    Article  PubMed  CAS  Google Scholar 

  69. Samuel CE. Antiviral actions of interferons. Clin Microbiol Rev. 2001; 14(4):778–809. https://0-doi-org.brum.beds.ac.uk/10.1128/CMR.14.4.778-809.2001.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Funding

This work was supported by the Niedersachsen-Research Network on Neuroinfectiology (N-RENNT) of the Ministry of Science and Culture of Lower Saxony. The funding body was not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

All expression data can be publicly derived from databases NCBI Gene Expression Omnibus or EMBL-EBI ArrayExpress with their corresponding Databank ID, listed in Table 2 and 4.

Author information

Authors and Affiliations

Authors

Contributions

KJ formulated the idea for meta-analysis of WNV related transcriptome expression data and supervised the analyses. RK and JD performed the data analyses. SCB and PC contributed with biological expertise in the interpretation of the results. All authors contributed in writing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Klaus Jung.

Ethics declarations

Ethics approval and consent to participate

No human, animal or plant derived strains have been used directly. This study is instead based on retrospectively analysed, publicly available data.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1

Additional figures A5.1-A5.5. (DOCX 420 kb)

Additional file 2

Result lists of differential analyses in group 1. The table displays result lists of the differential analyses in group 1. Rows representing those genes, common in the five individual studies of group 1. Two identifiers (‘GeneID’ and ‘Gene Symbol’) are used. Columns 3 to 18 are displaying the outcomes for different analysis variants: Columns, starting with the term ‘early’ are results from the ‘early merging’; same applies for the ‘late merging’. Columns of the individual studies begin with their identifier (see Fig. 2). For each, the p-values (‘.p’) and the adjusted p-values (‘.q’) were calculated. For the early and late merging, the absolute value of the test-statistic was also given out (‘.t.abs’). Columns ‘logFC.min’ and ‘logFC.max’ are showing the gene-wise minimum and maximum values of the logFCs from the individual studies. (XLSX 4787 kb)

Additional file 3

Result lists from GSEA analyses in group 1. The table displays result lists of the GSEA in group 1. Rows are representing those GO terms, for which at least two genes were available in each of the analysis variants. The GO-ID, as well as their specific name is given in the first column. The p-values (‘.p’) and adjusted p-values (‘.q’) are shown for each analysis variant (early, late and intermediate merging). The number of genes associated to the GO-term is presented in column ‘nPGenes’. The number of genes associated to the GO-term, which can be found in the data is presented in columns ‘early.nPGenes’ and ‘int.nPGenes’. (XLSX 1322 kb)

Additional file 4

Result lists of differential analyses in group 2. The table displays result lists of the differential analyses in group 2. Rows representing those genes, common in the two individual studies of group 2. Two identifiers (‘GeneID’ and ‘Gene Symbol’) are used. Columns 3 to 14 are displaying the outcomes for different analysis variants: Columns, starting with the term ‘early’ are results from the ‘early merging’; same applies for the ‘late merging’. Columns of the individual studies begin with their identifier (see Fig. 4). For each, the p-values (‘.p’) and the adjusted p-values (‘.q’) were calculated. For the early and late merging, the absolute value of the test-statistic was also given out (‘.t.abs’). Columns ‘logFC.min’ and ‘logFC.max’ are showing the gene-wise minimum and maximum values of the logFCs from the individual studies. (XLSX 5042 kb)

Additional file 5

Result lists from GSEA analyses in group 2. The table displays result lists of the GSEA in group 2. Rows are representing those GO terms, for which at least two genes were available in each of the analysis variants. The GO-ID, as well as their specific name is given in the first column. The p-values (‘.p’) and adjusted p-values (‘.q’) are shown for each analysis variant (early, late and intermediate merging). The number of genes associated to the GO-term is presented in column ‘nPGenes’. The number of genes associated to the GO-term, which can be found in the data is presented in columns ‘early.nPGenes’ and ‘int.nPGenes’. (XLSX 947 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kosch, R., Delarocque, J., Claus, P. et al. Gene expression profiles in neurological tissues during West Nile virus infection: a critical meta-analysis. BMC Genomics 19, 530 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-018-4914-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-018-4914-4

Keywords