The fold changes associated with the differentially expressed genes at day 14 post-infection were superimposed on the Chemokine signaling pathway and visualized using Cytoscape (Figure 4). Chemokine signaling clearly see more contributes to the upregulation of ISGs since the following signaling cascade is upregulated at the transcriptional level: Chemokine → Chemokine receptor (R) → JAK2/3 → STAT → ISG expression (Figure 4). Figure 4 Chemokine Signaling Pathway from the KEGG database (ID: mmu04062) overlaid

with log 2 fold change values for genes differentially expressed between DBA/2 and C57BL/6 at day 14. The scale for log2 fold change values is indicated at the bottom of the pathway diagram, where red shading indicates greater expression in DBA/2 compared to C57BL/6 mice and blue shading represents lesser expression. Genes not differentially expressed, i.e., with buy FK506 a fold change between −2 and +2 (log2 fold change between −1 and +1) are depicted in white. The identification of gene ontology (GO) terms significantly over-represented in the set of 1334 differentially expressed genes was performed using the Biological Networks Gene Ontology (BiNGO) tool [16], which preserves the hierarchical relationship among ontology

terms (Figure 5). Using an FDR corrected p-value cut-off <0.001 the three most significant FRAX597 GO terms were: immune system process, immune response, and defense response. Therefore, the immune related terms revealed by GO analysis agree with the results obtained from pathway analysis. The entire list of GO terms that were significantly enriched for differentially expressed genes at an FDR corrected p-value <0.05 are available in Additional file 2: Table S2. Figure 5 Hierarchical depiction of GO terms significantly

over-represented in the set of genes that were differentially expressed with a fold change ≥ 2 or ≤ -2 (log 2 fold change ≥ 1 or ≤ -1, respectively) between DBA/2 and C57BL/6 mice at any time point (N = 1334). The size of the node associated with each GO term is relative to the number of differentially Tyrosine-protein kinase BLK expressed genes belonging to that term. The color scale indicates the level of significance associated with each node with red being the most significant. For display purposes only GO terms with an FDR corrected p-value <0.001 are depicted. The full list of significant GO terms using an FDR corrected p-value cut off <0.05 is available in Additional file 2: Table S2. Protein network analysis Protein-protein and protein-DNA interactions between 416 genes that were differentially expressed between mice strains at day 14 were identified using MetaCore (GeneGo, St. Joseph, MI). The resulting protein interaction network depicted in Figure 6 consists of four major hubs: hypoxia inducible factor 1A (HIF1A), interferon regulatory factor 1 (IRF1), STAT1, and Yin Yang 1 (YY1).