- Research article
- Open Access
Influenza H3N2 infection of the collaborative cross founder strains reveals highly divergent host responses and identifies a unique phenotype in CAST/EiJ mice
BMC Genomics volume 17, Article number: 143 (2016)
Influenza A virus is a zoonotic pathogen that poses a major threat to human and animal health. The severe course of influenza infection is not only influenced by viral virulence factors but also by individual differences in the host response. To determine the extent to which the genetic background can modulate severity of an infection, we studied the host responses to influenza infections in the eight genetically highly diverse Collaborative Cross (CC) founder mouse strains.
We observed highly divergent host responses between the CC founder strains with respect to survival, body weight loss, hematological parameters in the blood, relative lung weight and viral load. Mouse strain was the main factor with highest effect size on body weight loss after infection, demonstrating that this phenotype was highly heritable. Sex represented another significant main effect, although it was less strong. Analysis of survival rates and mean time to death suggested three groups of susceptibility phenotypes: highly susceptible (A/J, CAST/EiJ, WSB/EiJ), intermediate susceptible (C57BL/6J, 129S1/SvImJ, NOD/ShiLtJ) and highly resistant strains (NZO/HlLtJ, PWK/PhJ). These three susceptibility groups were significantly different with respect to death/survival counts. Viral load was significantly different between susceptible and resistant strains but not between intermediate and highly susceptible strains. CAST/EiJ mice showed a unique phenotype. Despite high viral loads in their lungs, CAST/EiJ mice exhibited low counts of infiltrating granulocytes and showed increased numbers of macrophages in the lung. Histological studies of infected lungs and transcriptome analyses of peripheral blood cells and lungs confirmed an abnormal response in the leukocyte recruitment in CAST/EiJ mice.
The eight CC founder strains exhibited a large diversity in their response to influenza infections. Therefore, the CC will represent an ideal mouse genetic reference population to study the influence of genetic variation on the susceptibility and resistance to influenza infections which will be important to understand individual variations of disease severity in humans. The unique phenotype combination in the CAST/EiJ strain resembles human leukocyte adhesion deficiency and may thus represent a new mouse model to understand this and related abnormal immune responses to infections in humans.
Every year, about 500,000 people die world-wide from influenza virus infections . In recent history, the emergence of new influenza subtypes has caused several pandemics [2–4]. The most severe pandemic in 1918 resulted in about 30–50 million deaths worldwide  and a new variant of a seasonal H1N1 virus, pH1N1, caused a world-wide pandemic in 2009 [6, 7]. The course and outcome of an infection is determined by pathogen virulence, environmental and host factors. Genetic diversity of the host has been found as an important factor that strongly influences the host response and severity to influenza infections in humans [8–12]. Environmental and genetic factors can be well controlled in experimental animal models. It has been shown in mouse models that the genetic make-up of the host strongly influences the disease outcome after influenza infections (reviewed in [11, 13]). Subsequently, several genes have been described also in humans that contribute to disease severity .
Mouse genetic reference populations (GRPs) have become an important experimental system that allows investigating the influence of genetic variation on human diseases (reviewed in [15–21]. In addition, GRPs have been used to study the influence of genetic variation on the host response to infectious pathogens (reviewed in [13, 22–24]. The Collaborative Cross (CC) is a recently created population that has been generated from eight inbred mouse strains, referred to as founder . The addition of three wild-derived strains (CAST/EiJ, PWK/PhJ and WSB/EiJ) strongly enhanced the genetic diversity of the CC . Several studies in not fully inbred CC strains (pre-CC strains) and CC strains themselves led to the identification of quantitative trait loci (QTLs) influencing human diseases [27–31]. More recently, CC F1 hybrids have been identified that exhibited severe hemorrhagic fever, similar to humans, after infection with Ebola virus . Another CC strain has been shown to represent a new model for spontaneous colitis .
Here, we investigated the pathology, disease severity and the host response in the eight CC founder strains to H3N2 influenza A virus infections. We observed a broad range of phenotypic outcomes for many pathological and molecular parameters. Also, we found a specific defect in highly susceptible CAST/EiJ mice that resembled leukocyte recruitment deficiency in humans.
Three different groups of susceptibility phenotypes are observed after infection of CC founder strains with H3N2 virus
Eight to twelve weeks old female and male mice of the eight founder strains were infected intra-nasally with a dose of 1 × 101 FFU of the mouse-adapted influenza virus subtype H3N2 (A/HK/01/68). Survival rates of infected mice were monitored for a period of 14 days (Fig. 1). Our results suggest three different susceptibility groups based on survival rates and mean time to death: Highly susceptible, intermediate susceptible and highly resistant strains. Pair-wise chi-square tests between the three groups revealed significant differences for the death/survival counts (Additional file 1: Table S1). The strains A/J, CAST/EiJ and WSB/EiJ belonged to the highly susceptible group in which all infected mice died. NZO/HlLtJ and PWK/PhJ were found to be highly resistant; all infected mice survived the infection. Three strains, C57BL/6J, 129S1/SvImJ and NOD/ShiLtJ, showed intermediate susceptibility with survival rates ranging from 10 % to 90 %. No significant differences were observed in the survival curves of male and female mice for all strains except for NOD/ShiLtJ, whereby females were more susceptible than males. At an increased infection dose of 2 × 103 FFU all mice from highly and intermediate susceptible strains (A/J, C57BL/6J, 129S1/SvImJ, NOD/ShiLtJ, CAST/EiJ and WSB/EiJ) succumbed to the infection (Fig. 1). In contrast, the two highly resistant strains (NZO/HlLtJ and PWK/PhJ) were also resistant to infections with higher infection doses of 2 × 103 FFU and 2 × 105 FFU (Fig. 1). For all strains, except the highly resistant ones, survival rates in females and males at the lower dose were significantly higher than survival rates at higher doses (Fig. 1).
Next, we analyzed mean time to death (MTTD) to investigate how rapidly susceptible mice succumbed to infection. At an infection dose of 1 × 101 FFU, the highly susceptible strains exhibited the lowest MTTD. The MTTD was 7 and 9 days p.i. whereas the intermediate susceptible strains died at day 11 to 14 p.i. (Fig. 2a). At higher infection dose (2 × 103 FFU) the highly susceptible (day 5 to 6) and the intermediate susceptible group (day 7 to 8) had a lower MTTD (Fig. 2b). When comparing male and female mice, a significant difference of MTTD was found for NOD/ShiLtJ and WSB/EiJ mice at the low virus dose whereas at the high infection dose, only C57BL/6J mice exhibited a sex-dependent difference in MTTD (Fig. 2).
Infection of CC founder strains leads to highly variable body weight changes
In a previous study , body weight changes in CC mice were recorded daily until day 4 after infection with H1N1 (A/PR/8/34 (PR8)) virus. In our studies, we aimed to provide a more complete phenotype description after H3N2 infection that covers all phases of the host response (innate, adaptive and recovery). Therefore, we measured body weight changes at each day after infection until recovery of body weight after 14 days p.i.. We also included both sexes in our study. The analysis of body weight changes over 14 days confirmed a large phenotypic diversity between susceptible and resistant strains and also within susceptible and resistant groups (Fig. 3). In addition, the resistant strains were infected with a dose of 2 × 105 FFU. In this case, NZO/HlLtJ mice showed a drop in body weight (10 %) at day 2 which stayed low until day 14 after infection whereas PWK/PhJ showed higher body weight loss (10-20 %) at early days after infection and rapidly regained weight from day 4 on (Fig. 3).
A multi-factorial analysis of variance (ANOVA) was performed to investigate variation of body weight change as response to infection for the factors strains, sex and days. The analysis was limited to days 1 to 7 p.i. because for some strains (A/J, CAST/EiJ, and WSB/EiJ), more than 60 % of infected mice died by day 8. Thus, the surviving mice were not any more representative for the group. ANOVA (model: body weight loss ~ strain * sex * day and deleting the non-significant interaction strain:day:sex) revealed that all main effects (strain, sex and day) and all pair-wise interactions were significant (Additional file 2: Table S2). Subsequent Tukey-HSD post-hoc test was used to determine significance of pairs (strain by day p.i.; Additional file 3: Table S3). The changes in body weight were most pronounced on day 7 p.i.. Strains A/J, CAST/EiJ, and WSB/EiJ exhibited the highest loss of body weight whereas NZO/HlLtJ was least affected. The broad-sense heritability was determined for each day separately and ranged from 0.28 to 0.65 (Table 1) indicating high heritability of this trait. The sex-effect (using the optimized AVOVA model) was significant from days 3 to 7 p.i. but with a lower effect size (maximum of 8 %, Table 1).
Comparisons for all strains over 14 days p.i. (Fig. 3) confirmed that female mice from the highly and intermediate susceptible group exhibited a higher weight loss than male mice at the lower infection dose of 1 × 101 FFU (Fig. 3). For the resistant strains NZO/HlLtJ and PWK/PhJ no differences between sexes were observed (Fig. 3). In contrast, at the higher infection dose (2 × 103 FFU) no significant differences were observed between males and females except for some single days for C57BL/6J on day 3 and 4, PWK/PhJ on days 2, 11, 12 and 14 and WSB/EiJ on day 2 (Fig. 3). Male mice generally showed less severe changes in body weight than females. Only for the PWK/PhJ strain female mice exhibited significantly less body weight loss than male mice. However, this was only observed late after infection for days 11, 12 and 14.
Selection of infection dose and gender for further experiments
For all following studies, we decided to perform infections with a dose of 1 × 101 FFU because at this dose, variation between the founder strains was highest for survival, MTTD and body weight changes. Also, we limited the following analyses to female mice. Females were used extensively in our and other laboratories studies thus allowing better comparability to existing data than males. In addition, the sex effect was relatively small (ANOVA) compared to the strain effect. Also, in this study, we mainly focused on effects caused by genetic variation which is represented by the strain effect (heritability). Therefore, the additional information gained by including males would have been rather limited in our studies and, with respect to the 3R principle, was omitted.
Increase in relative lung weight differs among CC founder strains
The relative lung weight (wet lung weight relative to body weight on day of preparation) was investigated at different days after infection as measurement for the degree of immune cell infiltration and fluid influx into the lung (Fig. 4). For highly susceptible strains the relative lung weight increased significantly on day 5 p.i.. Lung weight increased for intermediate susceptible strains and reached a peak at day 8 p.i.. Resistant strains showed only a slight increase in relative lung weight on day 5 p.i. and 8 p.i.. To confirm that the increase in relative lung weight was not simply due to the loss in body weight at the respective days p.i., we also calculated relative lung weight for the different days p.i. with reference to the starting body weight. Also in this case, we still observed a significant increase in relative lung weight in susceptible strains and no significant increase in the two resistant strains (data not shown).
Kinetics of viral load in lungs revealed large differences between resistant and susceptible CC founder strains
To investigate if the severity of infection was related to viral replication in the lungs of mice, we determined viral load at different time points (day 3, 5 and 8) post infection. All highly and intermediate susceptible strains (A/J, CAST/EiJ, WSB/EiJ, C57BL/6J, 129S1/SvlmJ, NOD/ShiLtJ) showed similar high viral loads of about 105 FFU at all time points analyzed (Fig. 5). In contrast, the highly resistant strains NZO/HlLtJ and PWK/PhJ had much lower viral titers in the lung (Fig. 5). PWK/PhJ exhibited lowest viral loads of all strains (103 FFU/lung) and no increase between day 3 and 5 whereas NZO/HlLtJ mice showed an increase in viral load from day 3 to day 5 p.i. (Fig. 5). Both strains had cleared the virus from their lungs by day 8 p.i.. Viral load was significantly different on day 3 and day 5 p.i. between susceptible and resistant mice but not different between highly and intermediate susceptible mice (ANOVA, model: lg.viral.ld ~ suscept; Additional file 4: Table S4 and Additional file 5: Table S5).
Each CC founder strain exhibits a characteristic leukocyte profile in peripheral blood
The changes in cellular composition of the peripheral blood after infection was monitored via VetScan hematological assay (Fig. 6). After infection, the percentage of lymphocytes in all strains decreased, whereby the percentage of granulocytes increased. Only the highly resistant strains NZO/HlLtJ and PWK/PhJ showed nearly constant levels of both cell types. The percentage of monocytes was similar for all strains. Comparison between the eight founder strains revealed a high percentage of granulocytes in A/J, 129S1/SvImJ, NOD/ShiLtJ and WSB/EiJ on day 3 and 5 after infection, resulting in a ratio of about 35 % of the leukocytes. C57BL/6J, NZO/HlLtJ, CAST/EiJ and PWK/PhJ exhibited low levels of granulocytes at day 3, 5 and 8 after infection (20 % or below). CAST/EiJ was the only susceptible strain showing a decrease from day 3 to day 5. Notably, granulocyte amounts in mock infected NOD/ShiLtJ mice were at 18 % of total WBCs in mock infection. In conclusion, CAST/EiJ, although carrying high viral loads in the lung, exhibited only very low granulocyte counts in the periphery.
Correlation analysis of phenotypic traits
In our study, we measured changes of several phenotypic traits (body weight, hematology, relative lung weight, viral load) in single mice after infection. This allowed us to correlate these different traits within each subset. Fig. 7a shows the correlation analysis of body weight loss, viral load and percentages of cell populations in the peripheral blood for setup 2 (Table 2). Fig. 7b illustrates the correlation between body weight loss, relative lung weight and percentages of cell populations in the peripheral blood for setup 3 (Table 2). The analysis showed a strong negative correlation (coefficients of –0.98 and –0.97) between percentages of lymphocytes and granulocytes for both subsets which suggests that the main change in peripheral blood cell composition was mainly due to an increase in granulocytes (Fig. 7a, b). Furthermore, we found a strong negative correlation (coefficient of −0.87) between body weight loss and relative lung weight which suggests that the loss in body weight may serve as a surrogate for pathological changes in the lung (edema and immune cell infiltrates) (Fig. 7b). Body weight loss and percentage of lymphocytes appeared only weakly correlated with viral load (Fig. 7a, b). Similarly, body weight loss and percentage of lymphocytes were only moderately or weekly correlated (Fig. 7a, b).
Detailed analysis of unique CAST/EiJ phenotype by flow cytometry
CAST/EiJ mice were highly susceptible to H3N2 infection and showed a unique combination of high viral load in the lungs (Fig. 5) but low granulocyte counts in the periphery (Fig. 6). Therefore, a more detailed investigation of the CAST/EiJ mice in comparison to two intermediate susceptible mouse strains was performed by flow cytometry of blood and lung cells at day 3 and 5 p.i.. We chose C57BL/6J and 129S1/SvImJ mice for this comparison because 129S1/SvImJ mice exhibited a very strong response whereas C57BL/6J showed an intermediate peripheral response.
In blood, all three strains had comparable numbers of granulocytes (Ly6G+CD11b+ cells) before infection (Fig. 8a). At day 3 p.i., the absolute cell numbers of granulocytes increased in all strains, but much less in CAST/EiJ (1.5-fold) compared to C57BL/6J and 129S1/SvImJ mice (5-and 3-fold, respectively) (Fig. 8a). Furthermore, granulocyte counts significantly decreased from day 3 to day 5 in CAST/EiJ (Fig. 8a) but were similar to day 3 in the other two strains (Fig. 8a). After infection, monocyte (CD115+ cells) numbers remained high in CAST/EiJ mice at day 5 p.i., whereas they decreased in the other two strains from day 3 to day 5 p.i. (although not statistically significant) (Fig. 8a). Furthermore, NK cell (NKP46+ cells) numbers were lower in non-infected and at day 3 p.i. in CAST/EiJ compared to C57BL/6J and 129S1/SvImJ mice and at similar levels in CAST/EiJ mice compared to the other two strains on day 5 p.i. (Fig. 8a). Also, CAST/EiJ mice exhibited the highest absolute numbers of B cells (CD19+) in non-infected animals which strongly decreased after infection (Fig. 8a). CAST/EiJ mice also had low numbers of T-helper (CD4+ cells) and cytotoxic cells (CD8+ cells) compared to C57BL/6 J and 129S1/SvImJ before and after infection (Fig. 8a). Both T-cell populations diminished significantly in CAST/EiJ mice on day 3. In summary, the numbers of granulocytes in the blood were lower after infection in CAST/EiJ mice compared to the other two strains whereas the number of monocytes stayed high in the blood of CAST/EiJ mice.
The most remarkable difference that we observed in flow cytometry analysis of the lungs from infected CAST/EiJ mice (Fig. 8b) was the very low number of all investigated immune cell populations compared to C57BL/6 J and 129S1/SvImJ mice: granulocytes (Ly6G+ cells), dendritic cells (CD11c+cells), T cells (CD4+ and CD8+ cells), and macrophages (F4/80+ cells). Granulocytes in CAST/EiJ mice increased by 4-fold from day 0 to 3 p.i. with no further increase from day 3 to 5 p.i.. In contrast, granulocyte numbers in the lungs of 129S1/SvImJ and C57BL/6 J increased considerably from day 3 to day 5 p.i. and were 30-times higher on day 5 p.i. compared to CAST/EiJ. The number of macrophages was generally higher in 129S1/SvImJ and C57BL/6 J lungs compared to CAST/EiJ mice. However, the increase in macrophage numbers was much higher in CAST/EiJ (25-fold from day 0 to day 5) in comparison to 129S1/SvImJ and C57BL/6 J lungs (4-and 8-fold, respectively). The amount of dendritic cells (CD11c+) decreased similarly in all three strains after infection. CD4+ (T helper) cells as well as CD8+ (cytotoxic T) cells decreased until day 3 p.i. and then significantly increased until day 5 p.i. in all three strains. In summary, CAST/EiJ mice exhibited in general much lower numbers of immune cells in the lungs before and after infection. After infection, granulocyte numbers increased much less in CAST/EiJ lungs compared to the other two strains whereas the relative increase in macrophages was higher in CAST/EiJ.
Histopathological and immunohistochemical studies of lungs revealed differences in lung pathology of CAST/EiJ mice in comparison to C57BL/6J and 129S1/SvImJ
We investigated the lungs of CAST/EiJ mice in comparison to C57BL/6J and 129S1/SvImJ mice by histopathology and immunohistochemistry at different days after infection. C57BL/6J, 129S1/SvImJ and CAST/EiJ mice developed a multifocal mild to severe necrotizing and/or granulocytic broncho-interstitial pneumonia (Fig. 9, Table 3). In C57BL/6J mice at 3 day p.i., the alveolar septa were mildly thickened with infiltration of few granulocytes and there was mild to moderate necrosis of alveolar, bronchiolar and bronchial epithelial cells. In the alveolar, bronchiolar and bronchial lumina there were moderate inflammatory infiltrates consisting of granulocytes, increased numbers of alveolar macrophages, small amounts of proteinaceous material (edema) and fibrin, and mild hemorrhage (Table 3). Peribronchiolarly, peribronchially and perivascularly there was moderate edema and infiltration of moderate numbers of granulocytes and fewer macrophages, lymphocytes and plasma cells. At day 5 p.i., there was an increase in numbers of inflammatory cells in the alveolar lumina and amounts of intraalveolar edema and hemorrhage. Also, there was more alveolar and bronchiolar epithelial necrosis and pneumocyte type II hyperplasia. At day 8 p.i., in addition to previously described changes, there were intra-alveolar hyaline membranes and flattening of bronchiolar and bronchial epithelium. In C57BL/6J mock infected mice no significant lesions were observed. Lungs of 129S1/SvImJ mice showed less granulocyte infiltration and type II pneumocyte hyperplasia in the alveoli, less peribronchiolar, peribronchial, and perivascular infiltration, but slightly more necrosis of alveolar, bronchiolar and bronchial epithelium (Table 3). In lungs of CAST/EiJ mice we observed lower numbers of granulocytes and mild or no peribronchiolar, peribronchial and perivascular infiltrates. There was a trend for fewer intra-alveolar granulocytes and more intra-alveolar macrophages in CAST/EiJ mice compared to C57BL/6J and 129S1/SvImJ mice at day 5 p.i. (Additional file 6: Table S6). These observations were consistent with the flow cytometry data (Fig. 8) showing that CAST/EiJ mice exhibited an increase in granulocytes, macrophages and T lymphocytes whereas the absolute amount of leukocytes remained very low compared to the other susceptible strains. Immunohistochemistry showed that influenza virus antigen was present in respiratory epithelial cells of all strains (Fig. 9, Table 3). The number of positive alveolar epithelial cells peaked at 5 days p.i. in all strains, with about twice as many positive cells in the CAST/EiJ mice as in C57BL/6J or 129S1/SvImJ mice. The number of positive bronchiolar and bronchial epithelial cells was higher in CAST/EiJ mice (day 3 and 5 p.i.) than in C57BL/6J and 129S1/SvImJ mice. This observation revealed that cellular distribution of virus and number of infected cells was different in CAST/EiJ lungs compared to the other highly susceptible strains although the total viral load was similar (Fig. 5). This can be explained by the fact that CAST/EiJ lungs are smaller in comparison to the other analyzed strains. Thus, similar total amounts of infectious virus (viral load) in the lung of CAST/EiJ mice would result in a higher number of infected cells. In mock infected mice no virus antigen was detected in any tissue examined.
Analysis of global gene expression changes in lung and blood of CAST/EiJ mice revealed abnormal immune response after infection
To further investigate the high susceptibility of CAST/EiJ mice and to gain insights into the gene regulatory networks involved, we infected female CAST/EiJ, C57BL/6J, 129S1/SvlmJ and PWK/PhJ mice with 1 × 101 FFU H3N2 virus and performed transcriptome analysis of the lungs and peripheral blood at 3 and 5 days p.i. and mock day 3. Principal component analysis (PCA) of all individual samples from the lungs revealed good grouping of biological replicates and a clear separation of infected versus mock-infected mice (Fig. 10a). The first principal component reflected expression changes associated with disease whereas the second principal component was driven by mouse subspecies. C57BL/6J and 129S1/SvlmJ exhibited a change in gene expression from mock-infected to infected lungs at day 3 to day 5 p.i.. Also, CAST/EiJ infected mice did show a strong change from mock-to influenza infected samples. However, samples from days 3 and 5 p.i. were not separated from each other in CAST/EiJ mice in contrast to the changes in C57BL/6J and 129S1/SvlmJ. These observations indicate a strong peak inflammatory response already at day 3 p.i. in infected CAST/EiJ mice. PWK/PhJ mice exhibited a shift from mock-infected mice to day 3 and day 5 p.i. whereby the distance to mock-treated mice was larger at day 3 than at day 5 suggesting that the low replication rate of influenza virus in the lungs of PWK/PhJ mice induced a much lower inflammatory response which was already partly resolved at day 5 p.i.. Identification of differentially expressed probe sets (DEPS, showing a 2-fold or more absolute difference in expression compared to mock-treated mice with an adjusted p-value of < 0.05) revealed highest number of DEPS for infected versus mock-infected CAST/EiJ mice on day 3 (2054) and day 5 p.i. (2513), respectively (Fig. 11). C57BL/6J infected mice had intermediate numbers of DEPS that were at similar levels on day 3 p.i. (966 and 975, respectively) and higher in 129S1/SvlmJ (1627) than in C57BL/6J (1108), respectively, at day 5 p.i.. Lowest number of DEPS were observed for PWK/PhJ for both days, exhibiting higher numbers on day 3 (604) compared to day 5 p.i. (306). These results corroborate the conclusions above from the PCA analysis. Analysis of 13 inflammatory genes that are strongly regulated in influenza infected mice  revealed that all infected mice exhibited a similar level of up-regulation confirming that all mice were successfully infected. In addition, these results showed that CAST/EiJ mice were able to initiate a proper inflammatory host response (Additional file 7: Figure S1).
For the PCA of the blood transcriptomes, we first performed an ANOVA analysis and selected probe sets that were significantly regulated (p < 0.001). The resulting PCA of these 2079 probe sets revealed good grouping of infected and non-infected groups as well as for the different strains (Fig. 10b). Analysis of 13 inflammatory genes  showed similar up-regulation of the inflammatory response in CAST/EiJ mice compared to 129S1/SvlmJ and C57BL/6J. PWK/PhJ mice exhibited only weak up-regulation or even down-regulation of these genes (Additional file 8: Figure S2). The results again show that CAST/EiJ mice were capable to initiate a proper inflammatory response at the transcriptional level in the blood.
We also performed a Digital Cell Quantification (DCQ, ) to investigate possible differences in immune cell populations of CAST/EiJ mice. DCQ of blood from infected CAST/EiJ mice revealed lower numbers of granulocytes and a higher number of monocytes and macrophages compared to infected 129S1/SvlmJ mice (Fig. 12a). DCQ analysis of the lung transcriptome showed equal numbers of monocytes and macrophages in CAST/EiJ and 129S1/SvlmJ mice and a much stronger reduction in dendritic cells in CAST/EiJ compared to 129S1/SvlmJ mice after infection (Fig. 12b).
To identify genes that may be involved in the high susceptibility in CAST/EiJ mice, we determined probe sets that were up-regulated in the lungs of C57BL/6J and 129S1/SvlmJ at days 3 and 5 p.i. compared to the mock control but not in CAST/EiJ. However, it should be noted that differences in expression signals from array analysis of the genetically diverse CC mice may also be caused by SNPs in probes. We therefore confirmed expression changes obtained from an array data set with expression changes observed by RNAseq . Only those genes were retained that were differentially regulated in both, our arrays and the RNAseq study. The filtering revealed eight genes that failed to be up-regulated in CAST/EiJ mice (Fig. 13): Prm1 (protamine 1), Gdf3 (growth differentiation factor 3), Nts (neurotensin), Plekhs1 (pleckstrin homology domain containing, family S member 1), Hpse (heparanase), F830016B08Rik (RIKEN cDNA F830016B08 gene), Insrr (insulin receptor-related receptor), Wfdc17 (WAP four-disulfide core domain 17). In addition, the analysis of the chemokine/cytokine KEGG pathway showed a strong down-regulation of Ccr6 (chemokine (C-C motif) receptor 6) specifically in CAST/EiJ mice (Fig. 14).
We subsequently investigated other genes involved in the regulatory networks of Hpse and Ccr6 (see discussion) and found specific gene expression changes of Ccr7 (chemokine (C-C motif) receptor 7), Cxcr5 (chemokine (C-X-C motif) receptor 5), Cxcr1 (chemokine (C-X-C motif) receptor 1), and Hspg2 (perlecan, heparan sulfate proteoglycan 2) in CAST/EiJ mice after infection (Fig. 14).
We also filtered for genes that were up-regulated in the blood of C57BL/6J and 129S1/SvlmJ at days 3 and 5 p.i. but not in CAST/EiJ and found four DEPS: Il1b (interleukin 1 beta), Pi16 (peptidase inhibitor 16) and I830012O16Rik (RIKEN cDNA I830012O16 gene). The fourth probe set belonged to a gene from which a second probe set was not lower in CAST/EiJ mice and it was therefore omitted (Additional file 9: Figure S3).
Here, we have analyzed various phenotypic traits after influence A (H3N2) infection in the CC founder strains. We found a strong influence of genetic variation on all traits and a large variation of the host response and pathology. In particular, the response in CAST/EiJ mice revealed a unique phenotype which may represent a new model for the human leucocyte recruitment deficiency disease.
Susceptibility and resistance to influenza infection and disease is strongly influenced by genetic variation
Our comprehensive analysis of the CC founder strains revealed highly variable phenotypes for all investigated traits (Fig. 15). On the basis of mortality rate, MTTD, and body weight changes, we could distinguish three different susceptibility groups: highly susceptible strains, intermediate susceptible strains, and highly resistant strains. Amongst the resistant strains, NZO/HlLtJ only marginally lost weight (less than 10 %) at the low infection dose but exhibited weight loss at the higher doses whereas PWK/PhJ mice also exhibited a larger loss of body weight early after infection.
In previous studies, CC founders and incipient lines of the CC had been studied after infection with H1N1 (A/PR/8/34 (PR8)) virus and the authors described a large variation on phenotypes [34, 37]. Here, we investigated the CC founder strains after infection with another contemporary influenza virus subtype, H3N2 (A/HK/01/68), and we also observed high diversity of phenotypic responses between the founder strains. The study by Ferris et al., phenotyped infected mice only until day 4 p.i. . They identified a highly susceptible and a highly resistant group. Similarly, Bottomly and colleagues categorize a severe and mild response to infection based on global transcriptome profiling . In our studies, we monitored the infected mice over a much larger period, 14 days p.i., and also performed infections with different doses of the same virus, and we distinguish three susceptibility groups. Analysis of variance (ANOVA) for body weight changes revealed that the factor strain had the strongest effect on body weight loss indicating high heritability for this trait. These results thus demonstrate that CC founder strains and CC strains will represent a highly valuable model system to identify genetic factors that regulate susceptibility to infections and other host responses. Indeed, in an analysis of incipient pre-CC strains, several chromosomal regions were identified that influence gene expression of the host response as well as pathology in the lung [34, 37].
Influence of sex on disease progression and outcome
Using ANOVA, we also found that sex had also a significant effect on body weight loss. At a dose of 1 × 101 FFU female mice exhibited higher body weight loss compared to males. However, the effect of sex was weaker than the strain effect. The influence of sex had not been reported in previous studies which were performed only in females [34, 37]. In pairwise comparisons of strains, we found significant sex-specific differences in body weight loss in the highly and intermediate susceptible group after infection with 1 × 101 FFU, but not in the highly resistant group. Differential responses to H3N2 influenza A infections between female and male mice have been reported previously [38, 39]. Also sex-specific differences to influenza A virus infections in humans have been observed. However, the present knowledge is rather limited because epidemiological data often lack information on sex. Furthermore, tendencies showing that either females or males are more prone to severe disease progression are influenced by additional factors like general life style, social behavior, coinfections or viral exposure (reviewed in: ). Our results suggest that influenza infections in the CC population may represent a highly valuable resource to obtain a better understanding of sex-specific responses to pathogens in mammals, including humans.
We did not observe significant differences in survival between male and female mice in the founder stains, except for NOD/ShiLtJ. NOD/ShiLtJ mice spontaneously develop insulin-dependent type 1 diabetes mellitus (T1D) which is characterized by the destruction of pancreatic cells [41, 42]. They also show several abnormalities in immune cells, like defective antigen presentation , abnormalities in T  and NK cells  as well as impaired cytokine signaling in macrophages . Female NOD/ShiLtJ mice exhibit a strong decrease in insulin production already at the age of 12 weeks, which starts several weeks later in males. The incidence of spontaneous T1D is much higher in females (60-80 %) compared to males (20-30 %) . Thus, the increased susceptibility to influenza A infections in females which we observed may be directly related to their pre-diabetic/diabetic state or to an abnormal immune response. It has been shown, that there was an increased risk of complications among patients with diabetes after pH1N1 infection . Therefore, NOD/ShiLtJ mice may provide a valuable mouse model for investigating the role of diabetic state on influenza A infections.
Different Mx1 alleles influence survival rates and viral load but are modulated by host genetic background
Many studies on influenza pathogenesis in mice were performed in classical laboratory strains which carry a mutation in the influenza resistance gene Mx1 first described in the A2G mouse strain . The CC founder strains carry five different Mx1 alleles (, Fig. 16). The amino acid sequences of the protein encoded by the Mx1 alleles in NZO/HlLtJ and A2G were identical to the consensus sequence from the multiple alignments. The PWK/PhJ allele encodes one and the CAST/EiJ allele two amino acid changes compared to the consensus sequence. The WSB/EiJ transcript of Mx1 harbors a stop codon in its open reading frame leading to a truncated protein. All classical laboratory strains (A/J, C57BL/6J, 129S1/SvImJ and NOD/ShiLtJ) carry a deletion from the 8th to the 12th exon, resulting in a protein that has previously been shown to be non-functional [34, 49]. In our study, the presence or absence of a functional Mx1 allele correlated strongly with resistance and susceptibility. PWK/PhJ and NZO/HlLtJ were highly resistant and survived a virus challenge with increasing doses of infection, 1 × 101 FFU, 2 × 103 FFU and 2 × 105 FFU. Thus, the single amino acid difference had no influence on the protective functions of the MX1 protein. In contrast, the amino acid exchanges in the CAST/EiJ protein seemed to be detrimental leading to a highly susceptible phenotype. Similar results were obtained after infection with H1N1 (A/PR/8/34 (PR8)) virus . We recently found that the protective function of MX1 is dependent on genetic back ground . Thus, further experiments are needed to elucidate the function of the CAST/EiJ MX1 protein and the potential contribution of the CAST/EiJ genetic background.
MX1 is known to inhibit IAV replication [51, 52]. Therefore, we investigated viral load at different time points p.i. in the CC founder stains and related to the absence or presence of a functional Mx1 allele. Viral load was lowest in PWK/PhJ mice, was only detected at early times p.i. and was cleared by day 8 p.i.. NZO/HlLtJ mice were also able to clear the virus, but viral titer increased until day 5. Thus, there is a clear difference in viral load even in strains that carry a functional Mx1 allele which may be caused by genes that modulate the activity of MX1 or other immune response genes. In contrast, all intermediate and highly susceptible mice strains that carried the non-functional Mx1 alleles exhibited high viral loads already on day 3, and virus was not cleared on day 8 or at the time of death. The CC population now offers the possibility of investigating the same Mx1 allele on different genetic backgrounds and to identify potential modifying loci that interact with MX1.
Correlation of viral load with other phenotypes
We did not observe a general correlation between viral load and survival in the susceptible strains. All susceptible and intermediate susceptible strains had similar viral loads. For other mouse strains it had been shown that high viral loads in the lung correlate well with increased susceptibility to IAV. For example, the highly susceptible DBA/2 J strain exhibited a higher viral load compared to resistant C57BL/6J mice [53, 54]. Boon et al.  studied viral load in 21 inbred mouse laboratory strains infected with H5N1 (A/HK/213/03). They found significant association of viral load and mortality whereby an increased host inflammatory response was caused by a higher viral load. However, other studies demonstrate no differences in viral load between resistant and susceptible strains . Ferris et al. concluded from studies of pre-CC mice that both viral replication and virus-induced inflammation were predictors for disease severity being independent of each other. Nevertheless, when both were combined they were better predictors than either of them alone . Thus, in conclusion, either increased viral load causing extensive tissue damage or a hyper-inflammatory host response may result in a lethal outcome of an influenza infection. Therefore, the analysis of the divergent host responses in CC mice will allow investigating the biological mechanisms and the cellular and molecular networks in detail in a controlled experimental system and thus substantially contribute to better understand individual host factors that are involved in causing severe influenza disease outcomes in humans.
CAST/EiJ mice exhibit a distinct deficiency in their immune response to influenza infections
Hematological analysis of the peripheral blood in infected CAST/EiJ mice showed only a slight increase in granulocytes in the peripheral blood compared to other strains that had high viral loads in their lungs. Thus, we hypothesized that CAST/EiJ may exhibit a specific deficiency in their response to influenza infection. We therefore studied the immune response in CAST/EiJ mice in more detail by analyzing immune cells in blood and lung in comparison to other strains. Flow cytometry studies revealed that the increase of granulocytes in the blood was much lower after infection in CAST/EiJ mice compared to the other two strains whereas the number of monocytes stayed high in the blood of CAST/EiJ mice. In the lung, CAST/EiJ mice exhibited an increase of granulocytes, macrophages and T lymphocytes. However, the absolute amount of leukocytes remained very low compared to the other susceptible strains. After infection, granulocyte numbers increased much less in CAST/EiJ lungs compared to the other two strains whereas the relative increase in macrophages was higher in CAST/EiJ. Histological analysis confirmed the low level of neutrophil infiltrations in the lung and increased macrophage numbers. Also, DCQ analysis of global gene expression changes confirmed the low neutrophil response in the blood, reduced neutrophil infiltration into the lung and an increase in macrophage numbers in the lung. In line with our results, Ferris et al. showed that CAST/EiJ mice developed high viral loads in lungs after influenza infection but showed less infiltrating immune cells in their lungs . They and others observed low levels of T and NK cells in the blood of CAST/EiJ mice [34, 57]. CAST/EiJ mice are also highly susceptible to monkeypox virus (MPXV) infection lacking interferon gamma (IFN-γ) induction which might be related to a defect in activation and recruitment of IFN-γ producing immune cells . CAST/EiJ mice were rescued from lethal challenge with MPXV by intra-nasal administration of IFNγ . However, IFNγ is not crucial for the host response to influenza infection  whereas, interferons type I and III are essential [60–62]. Therefore, we investigated if pre-treatment with IFN-I was able to rescue CAST/EiJ mice from lethal influenza A infection. We observed a shift of about two to three days in the survival curve for IFN-treated versus non-treated CAST/EiJ mice. However, IFN pre-treatment did not rescue CAST/EiJ mice from lethal infection (data not shown). In addition, we did not observe reduced levels of Ifng transcripts in CAST/EiJ compared to the other mouse strains (Additional file 7: Figure S1 and Additional file 8: Figure S2).
Our observations suggest that granulocytes have a deficiency to enter the site of infection in CAST/EiJ mice. Members of the CXC chemokine family, notably CCL4 and CCL5, as well as selectins, play a major role in the recruitment of leukocytes, especially of neutrophils . We showed earlier that genes from these two families, notably CXC chemokines, are upregulated very early after influenza infection in lungs of C57BL/6J mice . However, we observed that after infection the transcript levels of Ccl4 and Ccl5 were similarly up-regulated in the lungs of CAST/EiJ as for 129S1/SvImJ and C57BL/6J (data not shown). Also, deficiencies in several CXC chemokine ligands and receptors (Cxcr2, Cxcl10, Cxcr3, Cxcl5) diminish the leukocytes influx into lungs [64–67]. The expression of Cxcl10, Cxcr3 and Cxcl5 was not differently regulated in CAST/EiJ compared to the two other susceptible mouse strains (data not shown). However, expression of Cxcr2 dropped to base line levels at day 5 p.i. only in CAST/EiJ mice (data not shown). It was reported for various inflammation models that in mice deficient in E-, L-, or P-selectin genes (Sele, Sell, Selp), the levels of neutrophils, T lymphocytes and lung inflammation are reduced [68–72]. All selectin genes were similarly up-regulated in all three susceptible genes in the lung after infection. However, Sele was strongly down-regulated in CAST/EiJ on day 5 p.i. (Additional file 10: Figure S4). Furthermore, the continuous recruitment of leukocytes is regulated by a feed-forward loop in which infiltrated immune cells in the lung enhance the release of chemokines leading to the recruitment of more neutrophils . The depletion of neutrophils resulted in uncontrolled virus growth and increased mortality in mice . Thus, reduced granulocyte infiltration in CAST/EiJ may secondarily also affect the subsequent recruitment of more granulocytes and other immune cells. However, we do observe an increase in macrophages. Further analyses will be necessary to unravel the details of these interactions phenomenon.
In summary, the amount of granulocytes in the blood of CAST/EiJ mice increased after infection which suggests a normal primary signaling of infected cells to the hematopoietic system and the release of leukocytes into the blood. However, the number of granulocytes remained constant in the lungs of CAST/EiJ from day 3 on and did not further increase. In contrast, granulocyte numbers strongly increased in lungs of 129S1/SvImJ and C57BL/6J mice. The strong specific down-regulation of Sele at day 5 supports the hypothesis of a leukocyte recruitment deficiency in CAST/EiJ mice. Together, these results suggest that the lack of sustained neutrophil recruitment may be the primary defect in CAST/EiJ mice.
We also performed a more detailed analysis of gene expression changes at the level of gene interaction pathways and individual genes to obtain more insights into the underlying mechanism of the abnormal CAST/EiJ immune cell phenotype. The canonical inflammatory genes, interferons and interferon response genes as well as many inflammatory chemokines such as Cxcl10, are all highly up-regulated in the lungs of CAST/EiJ mice. In addition, CAST/EiJ mice showed the highest number of DEPS after infection in blood and lungs. Thus, CAST/EiJ mice exhibited not a general deficiency in the overall activation of an inflammatory response but rather a specific defect in the recruitment of cells to the lung. Seven genes were specifically not up-regulated in CAST/EiJ after infection which have not been described in the context of the host response to influenza infections, but may be also important host factors for virus defense in humans . Furthermore, the Hpse gene was not as strongly up-regulated after infection in CAST/EiJ compared to the other mouse strains. Hpse encodes heparanase, a mammalian endoglycosidase which degrades heparan sulfate glycosaminoglycan in extracellular space and within the cells. It was postulated that heparanase may play an important role in inflammatory and autoimmune processes including leukocyte recruitment, immune cell extravasation and migration, release of cytokines and chemokines and activation of innate immune cells (reviewed in [75, 76]). Recently, the role of heparanase in regulating the migration of DCs was demonstrated in heparanase-deficient mice . When we analyzed the expression profiles of genes in the heparanase network, we found a specific down-regulation of the Hspg2 gene in CAST/EiJ mice. Hspg2 encodes perlecan [78–80] also known as basement membrane-specific heparan sulfate proteoglycan 2 (HSPG2). Perlecan is a large multi-domain proteoglycan that binds to and cross-links many extracellular matrix (ECM) components and cell-surface molecules. It is a key component of the vascular extracellular matrix, where it interacts with a variety of other matrix components helping to maintain the endothelial barrier function. Three long chains of glycosaminoglycans are attached to the perlecan core protein that plays an important role in the recruitment of leukocytes via binding to chemokines [81, 82]. Also, the Wfdc17 gene, encoding the AMWAP protein (activated microglia/macrophage whey acid protein domain protein), was not expressed in CAST/EiJ mice before and after infection. AMWAP has been postulated to be a negative regulator of pro-inflammatory microglia/macrophage activation and a potential modulator of innate immunity in neurodegeneration . Interestingly, Wfdc17 was up-regulated in lungs of C57BL/6J, 129S1/SvlmJ and PWK/PhJ after infection, but not in CAST/EiJ which is consistent with the observed increase of infiltrating macrophages in CAST/EiJ lungs. In addition, analysis of the KEGG chemokine/cytokine pathway revealed a strong down-regulation of Ccr6 (chemokine (C-C motif) receptor 6) in the lung that was unique for CAST/EiJ mice. CCR6 is an important chemokine receptor that mediates recruitment of immature and mature DCs and antigen presenting cells (APCs) to the sites of epithelia inflammation in the lung (reviewed in ). Epithelial cells in the lung increase their production of CCL20, which is the only chemokine ligand for CCR6, which in turn causes a further recruitment of immature CCR6-positive DCs into the lung. Ccr6 expression is down-regulated during the maturation process and accordingly Ccr7 is upregulated leading to the migration of DCs to the lymph node to fulfill its APC function. The expression of Ccr6 was down-regulated in all analyzed mouse strains, but very strongly in CAST/EiJ lungs, especially at day 5 after infection. Interestingly, Ccr7 expression was up-regulated in all strains, but not in CAST/EiJ, where Ccr7 was down-regulated. After infection a large number of immature DC continuously enter the lung tissue and then differentiate into mature DCs . This process was reflected by the up-regulation of Ccr7 in C57BL/6J and 129S1/SvlmJ. However, in CAST/EiJ mice less immature DC cells enter the lung which may be reflected by the decrease of Ccr6 and Ccr7 expression. Another down-regulated gene in CAST/EiJ was Cxcr1, a member of the G-protein-coupled receptor family. CXCR1 is a receptor for the CXC chemokine interleukin IL-8, which is a potent neutrophil recruiting  and activating factor [86, 87]. IL-8 is an important modulator of monocyte-endothelial interactions under flow conditions, since it can rapidly cause rolling monocytes to adhere firmly onto monolayers expressing E-selectin . Furthermore, Cxcr5 was down-regulated in CAST/EiJ lungs but not in the other strains. CXCR5 is a cytokine receptor that binds to B-lymphocyte attractant (BLC) and is involved in B-cell migration [89, 90].
Filtering for genes that were up-regulated in the blood of C57BL/6J and 129S1/SvlmJ but not in CAST/EiJ identified three genes. One of those Il1b (interleukin 1 beta) encodes IL1B which is an important inflammatory cytokine that mediates the inflammation and initiate the immune response against Influenza A virus infection . Pro-IL-1B is cleaved by caspase-1 which is activated through the formation of the NLRP3 inflammasome  and induces the expression of a variety of inflammatory mediators, which initiate the cascade of inflammatory responses . It was postulated that the synergy between CD40L and IL-1beta represents a potent, T cell-independent mechanism for DC activation during the earliest stages of inflammatory responses . Most interestingly, in humans it has been described that genetic variants in IL1A and IL1B contribute to the susceptibility to 2009 H1N1 infections .
In summary, the specific changes of gene expression patters that were unique in CAST/EiJ lungs belong to gene networks that play a crucial role in migration and recruiting of leukocytes. Together, these observations suggest a deficiency in CAST/EiJ mice in leukocyte adhesion and/or recruitment. The recruitment of leukocytes to the site of infection is a complex process of different steps including leucocyte rolling, activation, arrest and transmigration into the tissue which involves many proteins and interactions. Three types of leukocyte adhesion deficiencies (LADI-III) have been described in humans (reviewed in ). Mutations in ITGB2 (integrin beta 2), SLC35C1 (solute carrier family 35) and FERMT3 (fermitin family homolog 3) have been identified as the genetic cause of disease in LADI, II and III, respectively . However, we did not observe differences in the expression of these LAD genes in CAST/EiJ mice (data not shown). Thus the specific CAST/EiJ phenotype may be a valuable model system to identify and study additional genes related to human leukocyte adhesion and recruitment deficiencies. Of course, further functional analyses under different conditions will be necessary to describe the leukocyte phenotype in CAST/EiJ mice in more detail.
In conclusion, our studies on the CC founder strains emphasize the importance of host genetic backgrounds for the course and outcome of IAV infections. The CC founder strains as well as the CC strains will thus represent an extremely valuable experimental model system to understand the biological mechanisms as well as the genetic factors that cause severe influenza disease in humans. In most cases, high mortality was associated with high viral loads and a strong immune response. However, we also observed unexpected phenotype combinations such as high viral load but a weak immune response. In particular, the CAST/EiJ strain may provide a novel model for human leukocyte adhesion and recruitment deficiencies.
All experiments in mice were approved by an external committee according to the national guidelines of the animal welfare law in Germany (BGBl. I S. 1206, 1313 and BGBl. I S. 1934). The protocol used in these experiments has been reviewed by an ethics committee (equivalent to US IACUC, according to the German law) and approved by the ‘Niedersächsisches Landesamt für Verbraucherschutz und Lebensmittelsicherheit, Oldenburg, Germany’ (Permit Numbers: 33.9.42502-04-051/09 and 3392 42502-04-13/1234).
Virus, mouse strains, infection
The mouse-adapted virus strain influenza A/HK/01/68,obtained from Otto Haller, University of Freiburg  was produced in the allantoic cavity of 10-day-old embryonated hen eggs for 48 h at 37 °C. The CC founder (A/J, C57BL/6J, 129S1/SvlmJ, NOD/ShiLtJ, NZO/HILtJ, CAST/EiJ, PWK/PhJ, WSB/EiJ) were purchased from Jackson Laboratories (Bar Harbor, ME) and bred in our animal facility in Braunschweig for two to six generations depending on the respective strain. All mice were maintained under specific pathogen free conditions and according to the German animal welfare law. The entire age range of mice was 7 to 13 weeks whereby 95.5 % of the mice had an age range of 8–12 weeks. Female and male mice were anesthetized by intra-peritoneal injection with a mixture of Ketamine/Xylazine (100 mg/ml Ketamine and 20 mg/ml Xylazine) in sterile natrium chloride solution. The doses were adjusted to the individual body weight using 200 μl/20 g body weight. Mice were then intra-nasally infected with 20 μl virus solution of sterile phosphate-buffered saline using different virus concentrations (1 × 101, 2 × 103 and 2 × 105 FFU). Body weight was determined daily as % of initial weight at day 0. Mice that lost more than 30 % of their body weight had to be killed for ethical reasons and were also scored as dead. Three different experimental setups with combination of different phenotyping parameters were used (Table 2).
Relative lung weight
Mock-infected and infected female mice (n = 3-6 per strain and time point) were anesthetized using 160 μl Isofluoran (cp-pharma, Burgdorf, Germany). After exsanguinating via puncture of the retro-bulbar vein plexus, lungs were collected and weighted directly. Relative lung weight was calculated as followed: (lung weight on day of preparation/body weight on day of preparation) × 100 and displayed as percentage of total body weight.
Viral load in lungs
Viral load in infected lungs was determined on MDCK II (Madin-Darby Canine Kidney II) cells using the FFU assay as described previously . Briefly, lungs of mice were homogenized and debris removed by centrifugation. Samples were stored in aliquots at −70 °C. Serial dilutions of lung homogenates were prepared and viral titers determined by the FFU assay.
RNA isolation for arrays
At different time points after virus infection (day 3 and 5) and mock infection (day 3 post sterile phosphate-buffered saline infection) lung and blood were taken for RNA isolation. Body weight was measured until the day of sample collection. Eye blood was collected in RNA protect animal blood tubes (Qiagen) and immediately stored at −20 °C. Lung RNA was isolated using Qiagen Midi Kit as described previously (104). Blood RNA was isolated with the RNeasy protect animal blood kit (Qiagen). RNA concentration was measured with the NanoDrop (Thermo Scientific).
For DNA microarray hybridization and analysis, the quality and integrity of the total RNA was controlled on a 2100 Bioanalyzer (Agilent Technologies; Waldbronn, Germany). 100 ng of total RNA were applied for Cy3-labelling reaction using the one color Quick Amp Labeling protocol (Agilent Technologies; Waldbronn, Germany). Labeled cRNA was hybridized to Agilent’s mouse 4 × 44 k microarrays for 16 h at 68 °C and scanned using the Agilent DNA Microarray Scanner. Expression values were calculated by the software package Feature Extraction 10.5.1.1 (Agilent Technologies; Waldbronn, Germany).
Data were analyzed using GraphPad Prism version 5.04 for Windows (GraphPad Software, San Diego, California). Mean and standard error of the mean (SEM) were calculated for all groups. Dunn’s post test with multiple testing correction was used for significance testing between individual groups (e.g. different gender, virus concentration and time points). Non-parametric Mann Whitney U test was used to determine p-values for the significance of differences between two groups. Log rank test was used to calculate significances for survival rates. Analysis of variance (ANOVA) was performed using the lm and aov functions of the R software package . A multi-factorial ANOVA model (weight-loss ~ strain * sex * day) with all interaction terms and with days (from day 1 to day 7 p.i.) as a categorical factor was used first. In this analysis, the three-way interaction term strain:day:sex was not significant and was therefore removed to generate an optimized ANOVA model (Additional file 2: Table S2). Subsequent Tukey HSD post-hoc test was used to identify pairwise significant differences. A one-way ANOVA was performed using the model weight-loss ~ strain for each individual day separately, and the broad-sense heritability was determined by calculating the interclass correlation : (MSB-MSW)/(MSB+(n-1)*MSW) with MSB = mean square between groups, MSW = mean square within groups, n = mean number of mice per strain. One way ANOVA (model lg.viral.ld ~ suscept, see (Additional file 4: Table S4 and Additional file 5: Table S5 for details) and Tukey HSD post-hoc test was used to determine differences of viral load between strains of different susceptibility (intermediate susceptible, highly susceptible and resistant strains). Correlations for phenotypic traits were calculated and illustrated using R scripts developed by W. Venables (http://goo.gl/nahmV), S. Turner (https://gist.github.com/stephenturner/3492773) and the graphics from the package PerformanceAnalytics (https://cran.r-project.org/).
Gene expression analysis
Array data were analyzed using the R software package . Pre-processing steps included background correction, quantile normalization and annotation using the MmAgilentDesign026655.db , limma , and Agi4 × 44PreProcess  packages. Principal component analysis (PCA) analysis were performed using the affycoretools package . Multi-group comparisons and identification of differentially expressed probe sets (DEPS) were performed with the LIMMA package  using BH correction for multiple testing . DEPS were identified based on an adjusted p-value of < 0.05 and exhibiting more than a two-fold difference in expression levels (abs [log2] > 1). Inflammatory genes that are expressed during influenza infections were selected based on our previous influenza transcriptome studies [35, 104]. GO, KEGG, Reactome enrichment analysis and cluster profiling was performed with the package clusterProfiler . For the identification of genes that exhibited no change in gene expression levels after infection in the lungs of CAST/EiJ mice compared to the other strains, the full expression data set was filtered for probe sets that were not regulated in CAST/EiJ (log-fold change absolute <0.5) but up-regulated in C57BL/6J and 129S1/SvlmJ mice at days 3 and 5 after infection (log-fold change absolute >1). In addition, the data set was filtered for probe sets that were down-regulated in CAST/EiJ (log-fold change < −0.5) but up-regulated in C57BL/6 J and 129S1/SvlmJ mice at days 3 and 5 after infection (log-fold change absolute >1). These filters returned 8 and 2 probe sets, respectively. It should, however, be noted that gene expression differences between strains which are observed in Agilent arrays may not always represent true differences in expression levels between strains but may also result from sequence polymorphisms in transcripts that result in mismatches between probes and transcripts, since the probes on the array are based on the C57BL/6J reference strain. Therefore, we validated true expression changes by analyzing the respective genes in a previously published RNAseq data set . Only genes that were also differentially expressed in the RNAseq data set were further investigated. In this way, eight genes were found to be regulated in CAST/EiJ versus the other strains. The same analysis was performed for blood transcriptomes. Four DEPS were not regulated in CAST/EiJ (log-fold change absolute <0.5) but up-regulated in C57BL/6J and 129S1/SvlmJ mice at days 3 and 5 after infection (log-fold change absolute >1). KEGG pathway analysis was performed using the pathview package in R . DCQ analysis was performed according to the algorithm described in . Gene expression values were pre-processed by subtracting the mean of the mock-infected mice from the infected mice for each mouse strain and by selecting the most variant probe sets for genes with multiple probe sets based on the standard deviation across all groups. DCQ analysis was then performed at the public server provided by the University of Tel Aviv (http://dcq.tau.ac.il/). The results were then visualized by heatmaps using the R software .
Tissues were stored in 10 % neutral-buffered formalin (lungs after careful inflation with formalin), embedded in paraffin, sectioned at 4 μm, and stained with hematoxylin and eosin (HE) for examination by light microscopy. Semi-quantitative assessment of influenza virus-associated inflammation by histopathology in the lung was performed as reported earlier  for the extent of alveolitis and alveolar damage we scored: 0, 0 %; 1, 1-25 %; 2, 25-50 %; 3, > 50 %. For the severity of alveolitis, bronchiolitis and bronchitis we scored: 0, no inflammatory cells; 1, few inflammatory cells; 2, moderate numbers of inflammatory cells; 3, many inflammatory cells. For the presence of alveolar edema, alveolar hemorrhage and type II pneumocyte hyperplasia we scored: 0, no; 1, yes. Finally, for the extent of peribronchial, peribronchiolar and perivascular infiltrates we scored: 0, none; 1, one to two cells thick; 2, three to ten cells thick; 3, more than ten cells thick. Slides were examined without knowledge of the animals.
For detection of influenza A virus antigen, tissues were stained with a primary antibody against the influenza A nucleoprotein as described previously . In each staining procedure, an isotype control was included as a negative control. Semi-quantitative assessment of influenza virus antigen expression in the lungs was performed as reported earlier  with minor modification: for the alveoli, 25 arbitrarily chosen, 20× objective fields of lung parenchyma per animal were examined by light microscopy for the presence of influenza virus nucleoprotein, without the knowledge of the identity of the animals. The score for each animal was presented as the percentage of positive fields. For the bronchi and bronchioles, the percentage of positively staining epithelium was estimated on every slide to provide the score per animal: 0, 0 %; 1, 1-25 %; 2, 25-50 %; 3, > 50 %.
The numbers of neutrophils and intraluminal macrophages were counted in 20 arbitrarily chosen, 40× objective fields per animal. Slides were examined without knowledge of the identity of the animals. Viable neutrophils were identified by immunohistochemistry staining with primary antibody Ly-6G (rat-anti-mouse Ly-6G, 25–5931, eBioscience, San Diego, USA). Binding of the primary antibody was detected using a peroxidase labeled rabbit-anti-rat IgG-HRP (p0450, DAKO, Denmark). Peroxidase activity was revealed using 3-amino-9-ethylcafrbazole (AEC) (Sigma, St Louis, MO, USA). Intraluminal macrophages were identified on the basis of their morphology and location as described previously .
For monitoring hematological parameters, female mice were anesthetized using 160 μl isoflourane at different time points (mock day 3, day 3, 5, 8, 18 and day 30) after infection. Time points were omitted from the analysis when less than 60 % of the infected mice survived infection within a strain. Blood was taken by puncture of the retrobulbar vein plexus, collected in micro tubes containing EDTA (Sarstedt) and immediately measured in the hematologic system VetScan HM5 (Abaxis). The absolute cell numbers of lymphocytes, granulocytes and monocytes were determined and calculated as percentage of white blood cells.
Immune cell populations of single animals were analyzed in cell suspensions derived from homogenized lungs or whole blood. Blood was collected as described for the haematological analysis. Antibody against CD16/32 (clone 2.4G2) was added to blood minimizing unspecific binding to Fcy receptors. Blood cells were directly stained in two different panels with the following antibodies: αCD4-PE (BD Bioscience, Europe), αCD8-APC, αCD11b-APC, αCD19-FITC, αCD115-Alexa Fluor 488, αLy-6G-PE-Cy7 and αNKp46-PerCP-eFluor 710 (eBioscience). After incubation at 4 °C for 30 minutes, in the dark, erythrocytes were lysed and cells fixed in FACS Lysing Solution (BD Bioscience), and washed twice with PBS/2%FCS. Lungs prepared from the same animals were passed through a 100 μm-cell strainer and a density centrifugation using Lympholyte M (Cedarlane) was performed. After washing the cells obtained from the interphase twice with PBS/2%FCS, Fcy receptors were blocked with αCD16/32 and subsequently stained with antibodies as described above. In addition αF4/80-FITC and αCD11c-PE (eBioscience) were used. To exclude dead cells, 1 μg propidium iodide was added to each sample. Samples were measured using an Accuri C6 flow cytometer (BD Bioscience). Data were analyzed using FlowJo version 7.6.5 (Tree Star, Ashland, Oregon). The marker expression and gating strategy of immune cell populations was done as described . Minor subpopulations like CD8 + DC were not considered.
Deposition of data
Array data have been deposited at the gene expression data base GEO (GSE74077), phenotype data will be deposited at the mouse phenome database (MPD, ) upon publication of the manuscript.
Fauci AS. Seasonal and pandemic influenza preparedness: science and countermeasures. J Infect Dis. 2006;194 Suppl 2:S73–6.
Klenk HD, Garten W, Matrosovich M. Molecular mechanisms of interspecies transmission and pathogenicity of influenza viruses: Lessons from the 2009 pandemic. Bioessays. 2011;33(3):180–8.
Kilbourne ED. Influenza pandemics of the 20th century. Emerg Infect Dis. 2006;12(1):9–14.
Russell CJ, Webster RG. The genesis of a pandemic influenza virus. Cell. 2005;123(3):368–71.
Johnson NP, Mueller J. Updating the accounts: global mortality of the 1918–1920 “Spanish” influenza pandemic. Bull Hist Med. 2002;76(1):105–15.
Bautista E, Chotpitayasunondh T, Gao Z, Harper SA, Shaw M, Uyeki TM, et al. Clinical aspects of pandemic 2009 influenza A (H1N1) virus infection. N Engl J Med. 2010;362(18):1708–19.
Faruqui F, Mukundan D. 2009 pandemic influenza: a review. Curr Opin Pediatr. 2010;22(4):530–5.
Albright FS, Orlando P, Pavia AT, Jackson GG, Cannon Albright LA. Evidence for a heritable predisposition to death due to influenza. J Infect Dis. 2008;197(1):18–24.
Gottfredsson M, Halldorsson BV, Jonsson S, Kristjansson M, Kristjansson K, Kristinsson KG, et al. Lessons from the past: familial aggregation analysis of fatal pandemic influenza (Spanish flu) in Iceland in 1918. Proc Natl Acad Sci U S A. 2008;105(4):1303–8.
Horby P, Sudoyo H, Viprakasit V, Fox A, Thai PQ, Yu H, et al. What is the evidence of a role for host genetics in susceptibility to influenza A/H5N1? Epidemiol Infect. 2010;138(11):1–9.
Ferris MT, Heise MT. Quantitative genetics in the study of virus-induced disease. Adv Virus Res. 2014;88:193–225.
Zuniga J, Buendia-Roldan I, Zhao Y, Jimenez L, Torres D, Romo J, et al. Genetic variants associated with severe pneumonia in A/H1N1 influenza infection. Eur Respir J. 2012;39(3):604–10.
Kollmus H, Wilk E, Schughart K. Systems biology and systems genetics-novel innovative approaches to study host-pathogen interactions during influenza infection. Curr Opin Virol. 2014;6C:47–54.
To KK, Zhou J, Chan JF, Yuen KY. Host genes and influenza pathogenesis in humans: an emerging paradigm. Curr Opin Virol. 2015;14:7–15.
Schughart K, Libert C, Kas MJ: Controlling complexity: the clinical relevance of mouse complex genetics. Eur J Hum Genet 2013:10.1038/ejhg.2013.1079.
Roberts A, Pardo-Manuel DeVillena F, Wang W, McMillan L, Threadgill DW. The polymorphism architecture of mouse genetic resources elucidated using genome-wide resequencing data: implications for QTL discovery and systems genetics. Mamm Genome. 2007;18(6–7):473–81.
Peters LL, Robledo RF, Bult CJ, Churchill GA, Paigen BJ, Svenson KL. The mouse as a model for human biology: a resource guide for complex trait analysis. Nat Rev Genet. 2007;8(1):58–69.
Chesler EJ, Lu L, Wang J, Williams RW, Manly KF. WebQTL: rapid exploratory analysis of gene expression and genetic networks for brain and behavior. Nat Neurosci. 2004;7(5):485–6.
Buchner DA, Nadeau JH. Contrasting genetic architectures in different mouse reference populations used for studying complex traits. Genome Res. 2015;25(6):775–91.
Solberg Woods LC. QTL mapping in outbred populations: successes and challenges. Physiol Genomics. 2014;46(3):81–90.
Hunter KW, Crawford NP. The future of mouse QTL mapping to diagnose disease in mice in the age of whole-genome association studies. Annu Rev Genet. 2008;42:131–41.
Iraqi FA, Athamni H, Dorman A, Salymah Y, Tomlinson I, Nashif A, et al. Heritability and coefficient of genetic variation analyses of phenotypic traits provide strong basis for high-resolution QTL mapping in the Collaborative Cross mouse genetic reference population. Mamm Genome. 2014;25(3–4):109–19.
Marquis JF, Gros P. Genetic analysis of resistance to infections in mice: A/J meets C57BL/6 J. Curr Top Microbiol Immunol. 2008;321:27–57.
Vidal SM, Malo D, Marquis JF, Gros P. Forward genetic dissection of immunity to infection in the mouse. Annu Rev Immunol. 2008;26:81–132.
Collaborative Cross Consortium. The genome architecture of the collaborative cross mouse genetic reference population. Genetics. 2012;190(2):389–401.
Keane TM, Goodstadt L, Danecek P, White MA, Wong K, Yalcin B, et al. Mouse genomic variation and its effect on phenotypes and gene regulation. Nature. 2011;477(7364):289–94.
Aylor DL, Valdar W, Foulds-Mathes W, Buus RJ, Verdugo RA, Baric RS, et al. Genetic analysis of complex traits in the emerging Collaborative Cross. Genome Res. 2011;21(8):1213–22.
Mathes WF, Aylor DL, Miller DR, Churchill GA, Chesler EJ, De Villena FP, et al. Architecture of energy balance traits in emerging lines of the Collaborative Cross. Am J Physiol Endocrinol Metab. 2011;300(6):E1124–34.
Vered K, Durrant C, Mott R, Iraqi FA. Susceptibility to klebsiella pneumonaie infection in collaborative cross mice is a complex trait controlled by at least three loci acting at different time points. BMC Genomics. 2014;15:865.
Kelada SN, Aylor DL, Peck BC, Ryan JF, Tavarez U, Buus RJ, et al. Genetic analysis of hematological parameters in incipient lines of the collaborative cross. G3 (Bethesda). 2012;2(2):157–65.
Xiong H, Morrison J, Ferris MT, Gralinski LE, Whitmore AC, Green R, et al. Genomic Profiling of Collaborative Cross Founder Mice Infected with Respiratory Viruses Reveals Novel Transcripts and Infection-Related Strain-Specific Gene and Isoform Expression. G3 (Bethesda). 2014;4(8):1429–44.
Rasmussen AL, Okumura A, Ferris MT, Green R, Feldmann F, Kelly SM, et al. Host genetic diversity enables Ebola hemorrhagic fever pathogenesis and resistance. Science. 2014.
Rogala AR, Morgan AP, Christensen AM, Gooch TJ, Bell TA, Miller DR, et al. The Collaborative Cross as a Resource for Modeling Human Disease: CC011/Unc, a New Mouse Model for Spontaneous Colitis. Mamm Genome. 2014;25(3–4):95–108.
Ferris MT, Aylor DL, Bottomly D, Whitmore AC, Aicher LD, Bell TA, et al. Modeling host genetic regulation of influenza pathogenesis in the collaborative cross. PLoS Pathog. 2013;9(2), e1003196.
Pommerenke C, Wilk E, Srivastava B, Schulze A, Novoselova N, Geffers R, et al. Global transcriptome analysis in influenza-infected mouse lungs reveals the kinetics of innate and adaptive host immune responses. PLoS One. 2012;7(7), e41169.
Altboum Z, Steuerman Y, David E, Barnett-Itzhaki Z, Valadarsky L, Keren-Shaul H, et al. Digital cell quantification identifies global immune cell dynamics during influenza infection. Mol Syst Biol. 2014;10:720.
Bottomly D, Ferris MT, Aicher LD, Rosenzweig E, Whitmore A, Aylor DL, et al. Expression quantitative trait Loci for extreme host response to influenza a in pre-collaborative cross mice. G3 (Bethesda). 2012;2(2):213–21.
Boivin GA, Pothlichet J, Skamene E, Brown EG, Loredo-Osti JC, Sladek R, et al. Mapping of clinical and expression quantitative trait loci in a sex-dependent effect of host susceptibility to mouse-adapted influenza H3N2/HK/1/68. J Immunol. 2012;188(8):3949–60.
Robinson DP, Lorenzo ME, Jian W, Klein SL. Elevated 17beta-estradiol protects females from influenza A virus pathogenesis by suppressing inflammatory responses. PLoS Pathog. 2011;7(7), e1002149.
Gabriel G, Arck PC. Sex, immunity and influenza. J Infect Dis. 2014;209 Suppl 3:S93–9.
Anderson MS, Bluestone JA. The NOD mouse: a model of immune dysregulation. Annu Rev Immunol. 2005;23:447–85.
Makino S, Kunimoto K, Muraoka Y, Mizushima Y, Katagiri K, Tochino Y. Breeding of a non-obese, diabetic strain of mice. Jikken dobutsu Experimental animals. 1980;29(1):1–13.
Lee CN, Lew AM, Shortman K, Wu L. NOD mice are functionally deficient in the capacity of cross-presentation. Immunol Cell Biol. 2015;93(6):548–57.
Ferreira C, Palmer D, Blake K, Garden OA, Dyson J. Reduced regulatory T cell diversity in NOD mice is linked to early events in the thymus. J Immunol. 2014;192(9):4145–52.
Suwanai H, Wilcox MA, Mathis D, Benoist C. A defective Il15 allele underlies the deficiency in natural killer cell activity in nonobese diabetic mice. Proc Natl Acad Sci U S A. 2010;107(20):9305–10.
Fan H, Longacre A, Meng F, Patel V, Hsiao K, Koh JS, et al. Cytokine dysregulation induced by apoptotic cells is a shared characteristic of macrophages from nonobese diabetic and systemic lupus erythematosus-prone mice. J Immunol. 2004;172(8):4834–43.
Allard R, Leclerc P, Tremblay C, Tannenbaum TN. Diabetes and the severity of pandemic influenza A (H1N1) infection. Diabetes Care. 2010;33(7):1491–3.
Lindenmann J, Lance CA, Hobson D. The resistance of A2G mice to myxoviruses. J Immunol. 1963;90:942.
Staeheli P, Grob R, Meier E, Sutcliffe JG, Haller O. Influenza virus-susceptible mice carry Mx genes with a large deletion or a nonsense mutation. Mol Cell Biol. 1988;8(10):4518–23.
Shin DL, Hatesuer B, Bergmann S, Nedelko T, Schughart K. Protection from Severe Influenza Virus Infections in Mice Carrying the Mx1 Influenza Virus Resistance Gene Strongly Depends on Genetic Background. J Virol. 2015;89(19):9998–10009.
Gao S, von der Malsburg A, Dick A, Faelber K, Schroder GF, Haller O, et al. Structure of myxovirus resistance protein a reveals intra- and intermolecular domain interactions required for the antiviral function. Immunity. 2011;35(4):514–25.
Patzina C, Haller O, Kochs G. Structural requirements for the antiviral activity of the human MxA protein against Thogoto and influenza A virus. J Biol Chem. 2014;289(9):6020–7.
Dengler L, Kuhn N, Shin DL, Hatesuer B, Schughart K, Wilk E. Cellular changes in blood indicate severe respiratory disease during influenza infections in mice. PLoS One. 2014;9(7), e103149.
Srivastava B, Blazejewska P, Hessmann M, Bruder D, Geffers R, Mauel S, et al. Host genetic background strongly influences the response to influenza a virus infections. PLoS One. 2009;4(3), e4857.
Boon AC, Finkelstein D, Zheng M, Liao G, Allard J, Klumpp K, Webster R, Peltz G, Webby RJ: H5N1 Influenza Virus Pathogenesis in Genetically Diverse Mice Is Mediated at the Level of Viral Load. MBio 2011, 2 (5):pii: e00171-00111.
Davidson S, Crotta S, McCabe TM, Wack A. Pathogenic potential of interferon alphabeta in acute influenza infection. Nat Commun. 2014;5:3864.
Phillippi J, Xie Y, Miller DR, Bell TA, Zhang Z, Lenarcic AB, et al. Using the emerging Collaborative Cross to probe the immune system. Genes Immun. 2014;15:38–46.
Earl PL, Americo JL, Moss B. Lethal monkeypox virus infection of CAST/EiJ mice is associated with a deficient gamma interferon response. J Virol. 2012;86(17):9105–12.
Uetani K, Hiroi M, Meguro T, Ogawa H, Kamisako T, Ohmori Y, et al. Influenza A virus abrogates IFN-gamma response in respiratory epithelial cells by disruption of the Jak/Stat pathway. Eur J Immunol. 2008;38(6):1559–73.
Crotta S, Davidson S, Mahlakoiv T, Desmet CJ, Buckwalter MR, Albert ML, et al. Type I and type III interferons drive redundant amplification loops to induce a transcriptional signature in influenza-infected airway epithelia. PLoS Pathog. 2013;9(11), e1003773.
Koerner I, Kochs G, Kalinke U, Weiss S, Staeheli P. Protective role of beta interferon in host defense against influenza A virus. J Virol. 2007;81(4):2025–30.
Mordstein M, Neugebauer E, Ditt V, Jessen B, Rieger T, Falcone V, et al. Lambda interferon renders epithelial cells of the respiratory and gastrointestinal tracts resistant to viral infections. J Virol. 2010;84(11):5670–7.
Kolaczkowska E, Kubes P. Neutrophil recruitment and function in health and inflammation. Nat Rev Immunol. 2013;13(3):159–75.
Herbold W, Maus R, Hahn I, Ding N, Srivastava M, Christman JW, et al. Importance of CXC chemokine receptor 2 in alveolar neutrophil and exudate macrophage recruitment in response to pneumococcal lung infection. Infect Immun. 2010;78(6):2620–30.
Wang W, Yang P, Zhong Y, Zhao Z, Xing L, Zhao Y, et al. Monoclonal antibody against CXCL-10/IP-10 ameliorates influenza A (H1N1) virus induced acute lung injury. Cell Res. 2013;23(4):577–80.
Balamayooran G, Batra S, Cai S, Mei J, Worthen GS, Penn AL, et al. Role of CXCL5 in leukocyte recruitment to the lungs during secondhand smoke exposure. Am J Respir Cell Mol Biol. 2012;47(1):104–11.
Huang J, Li Z, Yao X, Li Y, Reng X, Li J, et al. Altered Th1/Th2 commitment contributes to lung senescence in CXCR3-deficient mice. Exp Gerontol. 2013;48(8):717–26.
Arndt PG, Strahan B, Wang Y, Long C, Horiuchi K, Walcheck B. Leukocyte ADAM17 regulates acute pulmonary inflammation. PLoS One. 2011;6(5), e19938.
Harp JR, Onami TM. Naive T cells re-distribute to the lungs of selectin ligand deficient mice. PLoS One. 2010;5(6), e10973.
Pascual DW, Wang X, Kochetkova I, Callis G, Riccardi C. The absence of lymphoid CD8+ dendritic cell maturation in L-selectin−/− respiratory compartment attenuates antiviral immunity. J Immunol. 2008;181(2):1345–56.
Jiao A, Fish SC, Mason LE, Schelling SH, Goldman SJ, Williams CM. A role for endothelial selectins in allergic and nonallergic inflammatory disease. Annals of allergy, asthma & immunology : official publication of the American College of Allergy, Asthma, & Immunology. 2007;98(1):83–8.
Hanna S, Etzioni A. Leukocyte adhesion deficiencies. Ann N Y Acad Sci. 2012;1250:50–5.
Brandes M, Klauschen F, Kuchen S, Germain RN. A systems analysis identifies a feedforward inflammatory circuit leading to lethal influenza infection. Cell. 2013;154(1):197–212.
Tumpey TM, Garcia-Sastre A, Taubenberger JK, Palese P, Swayne DE, Pantin-Jackwood MJ, et al. Pathogenicity of influenza viruses with genes from the 1918 pandemic virus: functional roles of alveolar macrophages and neutrophils in limiting virus replication and mortality in mice. J Virol. 2005;79(23):14933–44.
Goldberg R, Meirovitz A, Hirshoren N, Bulvik R, Binder A, Rubinstein AM, et al. Versatile role of heparanase in inflammation. Matrix biology: journal of the International Society for Matrix Biology. 2013;32(5):234–40.
Vlodavsky I, Iozzo RV, Sanderson RD. Heparanase: multiple functions in inflammation, diabetes and atherosclerosis. Matrix biology: journal of the International Society for Matrix Biology. 2013;32(5):220–2.
Poon IK, Goodall KJ, Phipps S, Chow JD, Pagler EB, Andrews DM, et al. Mice deficient in heparanase exhibit impaired dendritic cell migration and reduced airway inflammation. Eur J Immunol. 2014;44(4):1016–30.
Arikawa-Hirasawa E, Le AH, Nishino I, Nonaka I, Ho NC, Francomano CA, et al. Structural and functional mutations of the perlecan gene cause Schwartz-Jampel syndrome, with myotonic myopathy and chondrodysplasia. Am J Hum Genet. 2002;70(5):1368–75.
Iozzo RV. Perlecan: a gem of a proteoglycan. Matrix biology: journal of the International Society for Matrix Biology. 1994;14(3):203–8.
Kallunki P, Eddy RL, Byers MG, Kestila M, Shows TB, Tryggvason K. Cloning of human heparan sulfate proteoglycan core protein, assignment of the gene (HSPG2) to 1p36.1----p35 and identification of a BamHI restriction fragment length polymorphism. Genomics. 1991;11(2):389–96.
Laguri C, Arenzana-Seisdedos F, Lortat-Jacob H. Relationships between glycosaminoglycan and receptor binding sites in chemokines-the CXCL12 example. Carbohydr Res. 2008;343(12):2018–23.
Schmitz N, Kurrer M, Bachmann MF, Kopf M. Interleukin-1 is responsible for acute lung immunopathology but increases survival of respiratory influenza virus infection. J Virol. 2005;79(10):6441–8.
Karlstetter M, Walczak Y, Weigelt K, Ebert S, Van den Brulle J, Schwer H, et al. The novel activated microglia/macrophage WAP domain protein, AMWAP, acts as a counter-regulator of proinflammatory response. J Immunol. 2010;185(6):3379–90.
Ito T, Carson WF, Cavassani KA, Connett JM, Kunkel SL. CCR6 as a mediator of immunity in the lung and gut. Exp Cell Res. 2011;317(5):613–9.
Rossi DL, Hurst SD, Xu Y, Wang W, Menon S, Coffman RL, et al. Lungkine, a novel CXC chemokine, specifically expressed by lung bronchoepithelial cells. J Immunol. 1999;162(9):5490–7.
Pease JE, Sabroe I. The role of interleukin-8 and its receptors in inflammatory lung disease: implications for therapy. American journal of respiratory medicine: drugs, devices, and other interventions. 2002;1(1):19–25.
Sallusto F, Mackay CR, Lanzavecchia A. The role of chemokine receptors in primary, effector, and memory immune responses. Annu Rev Immunol. 2000;18:593–620.
Gerszten RE, Garcia-Zepeda EA, Lim YC, Yoshida M, Ding HA, Gimbrone Jr MA, et al. MCP-1 and IL-8 trigger firm adhesion of monocytes to vascular endothelium under flow conditions. Nature. 1999;398(6729):718–23.
Patrussi L, Capitani N, Cannizzaro E, Finetti F, Lucherini OM, Pelicci PG, et al. Negative regulation of chemokine receptor signaling and B-cell chemotaxis by p66Shc. Cell death & disease. 2014;5, e1068.
Elsner RA, Ernst DN, Baumgarth N. Single and coexpression of CXCR4 and CXCR5 identifies CD4 T helper cells in distinct lymph node niches during influenza virus infection. J Virol. 2012;86(13):7146–57.
Yazdi AS, Guarda G, D’Ombrain MC, Drexler SK. Inflammatory caspases in innate immunity and inflammation. J Innate Immun. 2010;2(3):228–37.
Acosta-Rodriguez EV, Napolitani G, Lanzavecchia A, Sallusto F. Interleukins 1beta and 6 but not transforming growth factor-beta are essential for the differentiation of interleukin 17-producing human T helper cells. Nat Immunol. 2007;8(9):942–9.
Luft T, Jefford M, Luetjens P, Hochrein H, Masterman KA, Maliszewski C, et al. IL-1 beta enhances CD40 ligand-mediated cytokine secretion by human dendritic cells (DC): a mechanism for T cell-independent DC activation. J Immunol. 2002;168(2):713–22.
Kim KS, Jung H, Shin IK, Choi BR, Kim DH. Induction of interleukin-1 beta (IL-1beta) is a critical component of lung inflammation during influenza A (H1N1) virus infection. J Med Virol. 2015;87(7):1104–12.
Schmidt S, Moser M, Sperandio M. The molecular basis of leukocyte recruitment and its deficiencies. Mol Immunol. 2013;55(1):49–58.
Haller O, Arnheiter H, Lindenmann J. Natural, genetically determined resistance toward influenza virus in hemopoietic mouse chimeras. Role of mononuclear phagocytes. J Exp Med. 1979;150(1):117–26.
R_Core_Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria 2013, http://www.R-project.org/.
Rutledge H, Aylor DL, Carpenter DE, Peck BC, Chines P, Ostrowski LE, et al. Genetic Regulation of Zfp30, CXCL1, and Neutrophilic Inflammation in Mouse Lung. Genetics. 2014;198(2):735–45.
99. Carlson M: MmAgilentDesign026655.db: Agilent Chips that use Agilent design number 026655 annotation data (chip MmAgilentDesign026655).R package version 2.14.0. http://www.bioconductor.org/packages/release/data/annotation/html/MmAgilentDesign026655.db.html 2014.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.
MacDonald J. affycoretools: Functions useful for those doing repetitive analyses with Affymetrix GeneChips. 2008. R package version 1340.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: A practical and powerful apporach to multiple testing. J R Stat Soc. 1995;57:289–300.
Alberts R, Srivastava B, Wu H, Viegas N, Geffers R, Klawonn F, et al. Gene expression changes in the host response between resistant and susceptible inbred mouse strains after influenza A infection. Microbes Infect. 2010;12(4):309–18.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Luo W, Brouwer C. Pathview: an R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics. 2013;29(14):1830–1.
van den Brand JM, Stittelaar KJ, Van Amerongen G, Rimmelzwaan GF, Simon J, De Wit E, et al. Severity of pneumonia due to new H1N1 influenza virus in ferrets is intermediate between that due to seasonal H1N1 virus and highly pathogenic avian influenza H5N1 virus. J Infect Dis. 2010;201(7):993–9.
van den Brand JM, Stittelaar KJ, Leijten LM, Van Amerongen G, Simon JH, Osterhaus AD, et al. Modification of the ferret model for pneumonia from seasonal human influenza A virus infection. Vet Pathol. 2012;49(3):562–8.
van den Brand JM, Stittelaar KJ, Van Amerongen G, Reperant L, De Waal L, Osterhaus AD, et al. Comparison of temporal and spatial dynamics of seasonal H3N2, pandemic H1N1 and highly pathogenic avian influenza H5N1 virus infections in ferrets. PLoS One. 2012;7(8), e42343.
Grubb SC, Bult CJ, Bogue MA. Mouse phenome database. Nucleic Acids Res. 2014;42(1):D825–34.
This work was supported by intra-mural grants from the Helmholtz-Association (Program Infection and Immunity) and a research grant FluResearchNet (No. 01KI07137) from the German Ministry of Education and Research to KS. The PhD work of SRL was partially financed by a scholarship from the “Gesellschaft der Freunde der Tierärztlichen Hochschule Hannover e. V.”. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
We wish to thank Kirsten Marlen Kleemann for maintaining the CC founder colony, Stefanie Edler for excellent technical assistance, Esther Wilk for advice on immune cell phenotyping, Peter van Run, Viroscience lab, Erasmus Medical Center, Rotterdam for support with histological studies, and Irit Gat-Viks and Yael Steuerman, Cell Research and Immunology Department, Tel Aviv University, Tel Aviv, Israel for help with digital cell quantification. Part of this work has been performed as PhD thesis work (SRL) at the University of Veterinary Medicine, Hannover.
The authors declare that they do not have any conflict of interests.
SRL, CP, HK, RB and KS conceived and designed the experiments. SRL, CP infected the mice, RG performed the array analysis. SRL, HK, KS performed the analysis of the data. LD performed the flow cytometry analysis. JMAvdB and TK performed the histology analysis. HK, SRL and KS wrote the manuscript. All authors read and approved the final manuscript.
Heike Kollmus and Klaus Schughart contributed equally as last authors.
Pairwise chi-square tests for death/survival counts between resistant, intermediate susceptible and highly susceptible strains. (DOCX 55 kb)
ANOVA analysis of main effects and interactions on body weight loss. Df: degrees of freedom; Sum Sq: sum of squares, Mean sq: mean sum of squares, Pr: p-value. ANOVA model: body weight loss ~ strain * sex * day; after deleting the non-significant three way interaction strain:day:sex. (DOCX 52 kb)
Results from Tukey post-hoc HSD test for the interaction of strain:day of the ANOVA analysis from Table S2. (XLSX 83 kb)
Post-hoc Tukey-HSD test of differences in viral loads between different groups of susceptible strains at day 3 p.i.. Results from ANOVA analysis (model: lg.viral.ld ~ suscept); lg.viral.ld: log2 viral load on day 3 p.i.; suscept: categories for resistant strains (resist), intermediate susceptible strains (int_susc) and highly susceptible strains (hlg_susc) are shown. Pairwise comparisons: int_susc-hlg_susc: intermediate susceptible versus highly susceptible strains, etc.; diff: difference between means of log2 viral load, p adj: adjusted p-value using Tukey HSD test. n = 36. (DOCX 52 kb)
Post-hoc Tukey-HSD test of differences in viral load between different groups of susceptible strains at day 5 p.i. Results from ANOVA analysis (model: lg.viral.ld ~ suscept); lg.viral.ld: log2 viral load on day 5 p.i.; suscept: categories for resistant strains (resist), intermediate susceptible strains (int_susc) and highly susceptible strains (hlg_susc) are shown. Pairwise comparisons: int_susc-hlg_susc: intermediate susceptible versus highly susceptible strains, etc.; diff: difference between means of log2 viral load, p adj: adjusted p-value using Tukey HSD test. n = 41. (DOCX 52 kb)
Average scoring of numbers of granulocytes and alveolar macrophages in C57BL/6 J, 129S1/SvImJ and CAST/EiJ lungs on day 5 after infection with H3N2 influenza virus from the histopathological study. (DOCX 53 kb)
Gene expression changes of inflammatory genes induced by influenza infection in 129S1/SvImJ, C56BL/6J, PWK/PhJ and CAST/EiJ lungs. Elevated gene expression levels of selected cytokines and chemokines in 129S1/SvImJ, C56BL/6J, PWK/PhJ and CAST/EiJ reflect influenza A infection. Expression values represent normalized log2 transformed signal intensities at different time points p.i. relative to expression levels in mock-infected control mice. (JPG 713 kb)
Gene expression changes of inflammatory genes induced by influenza infection in 129S1/SvImJ, C56BL/6J, PWK/PhJ and CAST/EiJ blood. Elevated gene expression levels of selected cytokines and chemokines in 129S1/SvImJ, C56BL/6J, PWK/PhJ and CAST/EiJ reflect influenza A infection. Expression values represent normalized log2 transformed signal intensities at different time points p.i. relative to expression levels in mock-infected control mice. (JPG 715 kb)
Gene expression changes of genes up-regulated in 129S1/SvImJ and C56BL/6J but not in CAST/EiJ in blood. Changes in the expression levels of probe sets that were up-regulated in 129S1/SvImJ and C56BL/6J but not in CAST/EiJ mice. Expression values represent normalized log2 transformed signal intensities at different time points p.i. relative to expression levels in mock-infected control mice. (JPG 395 kb)
Gene expression of selectins in lungs. Expression values represent normalized log2 transformed signal intensities at different time points p.i. relative to expression levels in mock-infected control mice. Probeset ID Sele: A_51_P455326 (JPG 652 kb)
Marker descriptions of immune cells that were studied in the DCQ analysis. (DOCX 54 kb)
About this article
Cite this article
Leist, S.R., Pilzner, C., van den Brand, J.M. et al. Influenza H3N2 infection of the collaborative cross founder strains reveals highly divergent host responses and identifies a unique phenotype in CAST/EiJ mice. BMC Genomics 17, 143 (2016). https://doi.org/10.1186/s12864-016-2483-y
- Influenza susceptibility
- Collaborative cross
- Mouse model
- Immune response