- 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
- Sarah R. Leist1,
- Carolin Pilzner1,
- Judith M.A. van den Brand2,
- Leonie Dengler1,
- Robert Geffers3,
- Thijs Kuiken2,
- Rudi Balling4,
- Heike Kollmus†1 and
- Klaus Schughart†1, 5Email author
© Leist et al. 2016
- Received: 15 October 2015
- Accepted: 17 February 2016
- Published: 27 February 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.
- Influenza susceptibility
- Collaborative cross
- Mouse model
- Immune response
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
Infection of CC founder strains leads to highly variable body weight changes
ANOVA-Broad-sense heritability and sex effect per day
P-value sex effect
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
Kinetics of viral load in lungs revealed large differences between resistant and susceptible CC founder strains
Each CC founder strain exhibits a characteristic leukocyte profile in peripheral blood
Correlation analysis of phenotypic traits
Relative lung weight
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.
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
Average scoring of lesions by histopathology and of virus antigen expression by immunohistochemistry in lung tissue of C57BL/6 J, 129S1/SvImJ and CAST/EiJ mice infected with H3N2
Day post infection or mock control
Extent alveolitis/alveolar damage
Severity of alveolitis
Severity of alveolar edema
Presence of interstitial edema
Presence of type II pneumocyte hyperplasia
Severity of bronchiolitis/bronchitis
Degree of peribronchial/perivascular cuffing
Bronchioles and bronchi
Analysis of global gene expression changes in lung and blood of CAST/EiJ mice revealed abnormal immune response after infection
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 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
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
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.
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Fauci AS. Seasonal and pandemic influenza preparedness: science and countermeasures. J Infect Dis. 2006;194 Suppl 2:S73–6.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Kilbourne ED. Influenza pandemics of the 20th century. Emerg Infect Dis. 2006;12(1):9–14.PubMed CentralPubMedView ArticleGoogle Scholar
- Russell CJ, Webster RG. The genesis of a pandemic influenza virus. Cell. 2005;123(3):368–71.PubMedView ArticleGoogle Scholar
- Johnson NP, Mueller J. Updating the accounts: global mortality of the 1918–1920 “Spanish” influenza pandemic. Bull Hist Med. 2002;76(1):105–15.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Faruqui F, Mukundan D. 2009 pandemic influenza: a review. Curr Opin Pediatr. 2010;22(4):530–5.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Ferris MT, Heise MT. Quantitative genetics in the study of virus-induced disease. Adv Virus Res. 2014;88:193–225.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Buchner DA, Nadeau JH. Contrasting genetic architectures in different mouse reference populations used for studying complex traits. Genome Res. 2015;25(6):775–91.PubMed CentralPubMedView ArticleGoogle Scholar
- Solberg Woods LC. QTL mapping in outbred populations: successes and challenges. Physiol Genomics. 2014;46(3):81–90.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Collaborative Cross Consortium. The genome architecture of the collaborative cross mouse genetic reference population. Genetics. 2012;190(2):389–401.PubMed CentralView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Gabriel G, Arck PC. Sex, immunity and influenza. J Infect Dis. 2014;209 Suppl 3:S93–9.PubMedView ArticleGoogle Scholar
- Anderson MS, Bluestone JA. The NOD mouse: a model of immune dysregulation. Annu Rev Immunol. 2005;23:447–85.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Lindenmann J, Lance CA, Hobson D. The resistance of A2G mice to myxoviruses. J Immunol. 1963;90:942.PubMedGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.Google Scholar
- Davidson S, Crotta S, McCabe TM, Wack A. Pathogenic potential of interferon alphabeta in acute influenza infection. Nat Commun. 2014;5:3864.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Kolaczkowska E, Kubes P. Neutrophil recruitment and function in health and inflammation. Nat Rev Immunol. 2013;13(3):159–75.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Arndt PG, Strahan B, Wang Y, Long C, Horiuchi K, Walcheck B. Leukocyte ADAM17 regulates acute pulmonary inflammation. PLoS One. 2011;6(5), e19938.PubMed CentralPubMedView ArticleGoogle Scholar
- Harp JR, Onami TM. Naive T cells re-distribute to the lungs of selectin ligand deficient mice. PLoS One. 2010;5(6), e10973.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Hanna S, Etzioni A. Leukocyte adhesion deficiencies. Ann N Y Acad Sci. 2012;1250:50–5.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Iozzo RV. Perlecan: a gem of a proteoglycan. Matrix biology: journal of the International Society for Matrix Biology. 1994;14(3):203–8.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Yazdi AS, Guarda G, D’Ombrain MC, Drexler SK. Inflammatory caspases in innate immunity and inflammation. J Innate Immun. 2010;2(3):228–37.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Schmidt S, Moser M, Sperandio M. The molecular basis of leukocyte recruitment and its deficiencies. Mol Immunol. 2013;55(1):49–58.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- MacDonald J. affycoretools: Functions useful for those doing repetitive analyses with Affymetrix GeneChips. 2008. R package version 1340.Google Scholar
- 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.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Luo W, Brouwer C. Pathview: an R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics. 2013;29(14):1830–1.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Grubb SC, Bult CJ, Bogue MA. Mouse phenome database. Nucleic Acids Res. 2014;42(1):D825–34.PubMed CentralPubMedView ArticleGoogle Scholar