Skip to main content
  • Research article
  • Open access
  • Published:

Distinct gene loci control the host response to influenza H1N1 virus infection in a time-dependent manner

Abstract

Background

There is strong but mostly circumstantial evidence that genetic factors modulate the severity of influenza infection in humans. Using genetically diverse but fully inbred strains of mice it has been shown that host sequence variants have a strong influence on the severity of influenza A disease progression. In particular, C57BL/6J, the most widely used mouse strain in biomedical research, is comparatively resistant. In contrast, DBA/2J is highly susceptible.

Results

To map regions of the genome responsible for differences in influenza susceptibility, we infected a family of 53 BXD-type lines derived from a cross between C57BL/6J and DBA/2J strains with influenza A virus (PR8, H1N1). We monitored body weight, survival, and mean time to death for 13 days after infection. Qivr5 (quantitative trait for influenza virus resistance on chromosome 5) was the largest and most significant QTL for weight loss. The effect of Qivr5 was detectable on day 2 post infection, but was most pronounced on days 5 and 6. Survival rate mapped to Qivr5, but additionally revealed a second significant locus on chromosome 19 (Qivr19). Analysis of mean time to death affirmed both Qivr5 and Qivr19. In addition, we observed several regions of the genome with suggestive linkage. There are potentially complex combinatorial interactions of the parental alleles among loci. Analysis of multiple gene expression data sets and sequence variants in these strains highlights about 30 strong candidate genes across all loci that may control influenza A susceptibility and resistance.

Conclusions

We have mapped influenza susceptibility loci to chromosomes 2, 5, 16, 17, and 19. Body weight and survival loci have a time-dependent profile that presumably reflects the temporal dynamic of the response to infection. We highlight candidate genes in the respective intervals and review their possible biological function during infection.

Background

Influenza A virus represents a major health threat to humans. The 1918 H1N1 pandemic caused at about 30 to 50 million deaths [1]. Seasonal influenza epidemics cause high economic loss, morbidity and deaths every year [2]. The course and outcome of an influenza A virus infection is influenced by viral and host factors. Host risk factors, like obesity or pregnancy, became evident during the recent swine flu pandemics [3, 4]. Furthermore, genetic factors in humans associated with a higher susceptibility to influenza infections and severe disease outcome have been suspected for the 1918 pandemics, as well as the H5N1 human infections [5–7]. Recently, the importance of IFITM3 as a crucial factor for host susceptibility has been demonstrated in mice and humans [8].

The importance of host factors to host susceptibility and resistance has been demonstrated clearly in animal models. We and others have shown in mouse infection models that the susceptibility of the host to influenza A infection strongly depends on the genetic background [9–17]. In particular, DBA/2J mice are highly susceptible to many influenza A virus subtypes, including those that were directly derived from human isolates without prior adaptation to the mouse [9, 13, 16–18]. In contrast, C57BL/6J mice are more resistant. After infection with mouse-adapted H1N1 (PR8M virus), DBA/2J mice loss weight very rapidly and die within 5–7 days post infection (p.i.), whereas C57BL/6J mice loss weight until days 6–8 after infection and regain their initial weight by 14 days p.i. [9, 18]. Viral load in the lungs of DBA/2J infected mice is much higher and lung pathology is very severe compared to infected C57BL/6J mice. Also, the production of chemokines and cytokines is much higher in DBA/2J mice [9, 18].

However, the genomic regions that are responsible for the differential response after infection with H1N1 have not been determined. Therefore, we used a large family of BXD type recombinant inbred strains generated by crossing C57BL/6J (resistant) to DBA/2J (susceptible) to map genetic loci that modulate disease severity. The BXD genetic reference population (GRP) is made up of a set of progeny strains, each with a defined and fixed genetic architecture. It is one of the largest families of strains, consisting of about 80 fully inbred strains [19, 20] available from the Jackson Laboratory and a new set of 80 additional lines that are still in production at the University of Tennessee. Individuals within each single strain are essentially isogenic (except for the sex chromosomes) and genotypes for the entire family, including most of the new strains, are known [21]. Genetic variation among this family has be exploited extensively in the past to systematically study the genetics of many traits (for examples of phenotypes see the GeneNetwork database [22]).

Here, we infected over 50 of the BXD strains with influenza A H1N1 virus and monitored body weight, survival, and mean time to death for the following 13 days post infection. We identified two significant and several suggestive loci peaks for all three traits. All showed a time-dependent appearance. Data mining of the intervals revealed several candidate genes, several of which may be important for the host response to influenza A virus infection.

Results

Susceptibility to influenza A virus after experimental infection of BXD mouse strains is highly variable

We infected 53 recombinant inbred strains of the BXD population plus the parental strains C57BL/6J, DBA/2J and B6D2F1 mice with mouse-adapted H1N1 virus (PR8M, H1N1 [18]) and followed body weight and survival over the next 13 days. Mice that lost more than 25% of their starting weight were sacrificed and also recorded as dead. After infection, body weight changes were highly variable between the BXD strains over the period of 13 days post infection (Figure 1, Additional file 1: Table S1). Similarly, survival rates (Additional file 2: Figure S1A, B) and mean time to death (MTTD; Additional file 2: Figure S1C) were highly variable between the different BXD strains. Furthermore, the trait percent survival showed a strain-dependent progression over time.

Figure 1
figure 1

BXD strains exhibit variable kinetics of weight loss and survival after infection with Influenza A virus. Mice from 53 BXD and parental strains were infected intra-nasally with 2 × 103 FFU of PR8 virus. Weight loss and survival of infected mice was followed over a period of 13 days. Mortality includes mice that were sacrificed because they had lost more than 25% of body weight. Three phenotypic groups can be distinguished: in the first group (A), all infected mice within a given BXD strain survived, in the second group (B) less than 50% of infected mice within a given BXD strain died, in the third group (C), more than 50% of infected mice within a given BXD strain died. From the weight loss curve of the second and third group it also becomes obvious that non-surviving mice were all approaching the 75% body weight loss endpoint before dying.

For body weight loss and survival, three different phenotypic response groups can be defined (Figure 1). In the first (Figure 1A), all infected mice within a strain survived, in the second (Figure 1B) a majority but not all individuals within a strain survived, and in the third group (Figure 1C), a majority died. Most remarkably, in the first group four strains—BXD9, BXD13, BXD43, BXD97—were highly resistant indicating that the infection may not cause any major pathology (Figure 1A). In contrast, BXD28 belonging to the third group, lost body weight much more rapidly than even the highly susceptible DBA/2J parent (Figure 1C). These results illustrate a large variation of responses within the BXD family.

By day 7 p.i., all infected mice had succumbed to infection in 11 strains, whereas 14 others exhibited more limited mortality (Additional file 2: Figure S1A). The incidence of mortality increased in some strains from day 8 p.i. through day 11 but not thereafter (Additional file 2: Figure S1B). For the MTTD phenotype, 17 strains showed no mortality after infection, similar to the resistant C57BL/6J parent. Three strains—BXD28, BXD18, BXD103—exhibited a MTTD that was even shorter than for the susceptible DBA/2J parent (Additional file 2: Figure S1C). Although this study was conducted over a period of approximately three years, and although mice were received from different sources and different experimenters performed the infection experiments, we did not note any significant influence of these potential confounds and cofactors.

Principal Component Analysis (PCA) of strain mean body weight loss from day 1 until day 7 was carried out to reduce the number and redundancy of related measurements, as well as to evaluate whether this set of measurements can be broken down into statistically and genetically independent processes. This analysis revealed one major component, PC1 that explains most of the weight loss among strains (more than 80%, see Additional file 3: Figure S2). This component corresponds to differences in loss from day 3 to 7 p.i. (Figure 2). In contrast, PC2 accounts for only about 10% of the variance in weight loss (Additional file 3: Figure S2). This component corresponds only to very early weight loss after virus inoculation, on days 1 and 2 p.i.(Figure 2), well before any appreciable virus replication has occurred. Thus, PC2 is highly likely to be a technical, but still interesting effect associated with the stressful procedure of intranasal inoculation and full anesthesia. PC1 most likely reflects the biologically relevant variation among strains in response to viral replication.

Figure 2
figure 2

Principal component analysis of all body weight traits reveals two distinct groups. PCA loading plot of principal components PC1 and PC2 (Supplement Figure 2) for body weight traits form day 1 to day 7 p.i. shows two distinct groups. The variation of body weight loss at day 3 to day 7 (ID 13007 to ID 13011) are mainly represented by factor 1. On the other hand, variation of body weight loss at day 1 (ID 13005) and 2 (ID 13006) are most likely explained by factor 2.

Analysis of body weight loss revealed a significant QTL on chromosome 5 and several suggestive QTLs with time-dependent effects

We mapped body weight loss following infection day by day. For purposes of analysis, body weights of mice that had died or that were euthanized were assigned a weight equal to 75% of their initial weight. It should be noted that mice which died continuously lost weight and were close to 75% body weight loss before they were found dead (Figure 1B, C). After day 7, surviving mice started to gain weight (Figure 1A, B). For this reason, we limited our analyses of body weight traits to the period of day 1 to day 7 p.i. in order to avoid mixing data from dead with recovering mice. Interval mapping of body weight loss detected a significant time-dependent locus on chromosome 5 that we named ‘QTL for influenza virus resistance on chromosome 5’ (Qivr5) adopting the nomenclature proposed by [12]. This QTL exhibited genome-wide significance on days 5 (LRS: 19.0, effect size: 30%) and 6 (LRS: 19.4, effect size: 31%, Figure 3). However, the effect of Qivr5, although weaker, can also be detected on days 2, 3, and 4 and 7 (Figure 3). Most interestingly; the resistance allele at Qivr5 is inherited from the nominally sensitive DBA/2J parental strain, illustrating the genetic complexity of the influenza response. Suggestive loci (LRS values between 10 and 15) map to chromosomes 2, 6, 9, proximal and distal 16, and chromosome 17. The effect of these loci was also time-dependent (Figure 3A-G).

Figure 3
figure 3

Genome-wide linkage analysis for body weight loss after H1N1 infection identifies time-dependent QTLs. Analysis of body weight from day 1 (A), day 2 (B), day 3 (C), day 4 (D), day 5 (E), day 6 (F) and day 7 (G) after PR8 infection revealed time-dependent QTLs. A significant QTL is located on chromosome 5 and suggestive QTLs are observed on chromosomes 2, 6, 9, 10, 16 and 17. Dead or euthanized mice were included for the analysis by manually setting their body weight to 75%. Interval mapping is shown across the entire genome as blue line representing the LRS scores; green line: DBA/2J alleles increase trait values; red line: C57BL/6J alleles increase trait values. Upper x-axis shows mouse chromosomes, lower x-axis shows physical map in mega bases for each chromosome, y-axis represents LRS score. The horizontal lines mark the genome-wide significant thresholds at p<0.05 (red line) and suggestive thresholds at p<0.37 (gray line).

Another suggestive QTL peak was found on chromosome 10 on day 1 (LRS: 14.1, effect size: 23%, Figure 3A). Its effect is lower at day 2 and not apparent at later days. These observations are in accordance with the PCA (Figure 2) that reveals two separate time-dependent influences on body weight variance among strains. The QTL appeared at an early time point after infection, when virus replication has just begun and strong inflammatory host responses are not yet evident [9, 18, 23]. These observations indicate that the effect is most likely related to the experimental protocol, namely the stress to anesthesia and intra-nasal application as well as treatment recovery. Treatment-dependent QTLs were described previously [24]. It is worth to note that mock-infection of the parental DBA/2J and C57BL/6J mice did not lead to a lasting body weight loss over a longer time interval except for a slight drop in body weight on day 1 p.i. (Additional file 4: Figure S3).

We also performed a QTL analysis for the PC1 and PC2 described above. PC1 detects a significant QTL on chromosome 5 (LRS: 19, effect size: 30%, Figure 4A) as well as suggestive QTL peaks on proximal chromosome 16 and on chromosome 17 (Figure 4A). Thus, the PC1 confirmed the significant QTL on chromosome 5 found with the daily body weight loss traits. Interval mapping of PC2 detected no significant QTL and one suggestive QTL on chromosome 13 (Figure 4B) which was not seen in any other trait.

Figure 4
figure 4

Genome-wide linkage scan for PC1 and PC2 body weight loss. Interval mapping of principle component PC1 (A) for body weight loss traits at day 1 to 7 after PR8 infection shows significant QTLs on chromosome 5 and suggestive QTLs on chromosomes 16 and 17. Interval mapping of PC2 (B) for body weight loss traits reveals no significant QTLs and a suggestive QTL on chromosome 13. Interval mapping is shown across the entire genome as blue line representing the LRS scores; green line: DBA/2J alleles increase trait values; red line: C57BL/6J alleles increase trait values. Upper x-axis shows mouse chromosomes, lower x-axis shows physical map in mega bases for each chromosome, y-axis represents LRS score. The horizontal lines mark the genome-wide significant thresholds at p<0.05 (red line) and suggestive thresholds at p<0.37 (gray line).

Survival rate and mean time to death traits confirmed the QTL on chromosome 5 and detected another significant QTL on chromosome 19

The analysis of survival rate traits on days 7, 8 and 11 p.i. revealed significant peaks on chromosomes 5 and 19 and suggestive peaks on chromosomes 1, 2, 3, 10, 16 and 17 (Figure 5A-C). The effects of these QTLs were again time-dependent. Qivr5 was significant at day 7 (LRS: 24.8, effect size: 37%) and day 8 p.i. (LRS: 19.7, effect size: 31% Figure 5A, B) and its effect was still evident from day 9 until 11 p.i. (Additional file 5: Figure S4 and Figure 5C). Similarly, another significant QTL with time-dependent effects was observed on chromosome 19 (Qivr19) at day 8 (LRS: 19.4; effect size: 29%) that could also be detected at day 7 and days 9–11 (Figure 5A-C and Additional file 5: Figure S4). Qivr5 and Qivr19 resulted from a positive influence of DBA/2J on survival, whereas Qivr16 (distal locus) and Qivr17-2 resulted from a positive influence of C57BL/6J alleles to increase survival. Furthermore, MTTD analysis confirmed Qivr5 as a significant QTL (LRS: 20.5, effect size: 32%), and Qivr19 was almost significant for this trait (Figure 5D). In addition, suggestive QTLs were found on chromosomes 2 and 17 for the MTTD trait.

Figure 5
figure 5

Genome-wide linkage scan for survival rate and mean time to death confirms QTLs on chromosomes 5 and 19. Interval mapping of survival rates (in %) after PR8 infection at day 7 (A), day 8 (B) and day 11 (C) show significant QTLs on chromosomes 5 and 19. Suggestive QTLs are observed on chromosomes 2, 3, 10, 16 and 17. Interval mapping using mean time to death confirmed the significant QTL on chromosome 5 and revealed suggestive QTLs on chromosomes 2, 17 and 19 (D). Interval mapping is shown across the entire genome as blue line representing the LRS scores; green line: DBA/2J alleles increase trait values; red line: C57BL/6J alleles increase trait values. Upper x-axis shows mouse chromosomes, lower x-axis shows physical map in mega bases for each chromosome, y-axis represents LRS score. The horizontal lines mark the genome-wide significant thresholds at p<0.05 (red line) and suggestive thresholds at p<0.37 (gray line).

In conclusion, survival and MTTD traits confirmed the significant QTL on chromosome 5, revealed an additional significant QTL on chromosome 19 and several suggestive QTLs. All QTLs showed a time-dependent effect.

Composite interval and pair-scan mapping indicates various interactions of QTLs

The strongest QTLs map to chromosomes 5 and 19. We therefore performed composite interval mapping in which the contributions of these strong QTLs were factored out to reveal possible secondary QTLs. When we controlled for Qivr19, the linkage to distal chromosome 16 increased and became nearly significant for both survival at day 8 and MTTD (Figure 6A, B). The linkage to chromosome 2 increased slightly. However, control for the contribution of Qivr5 did not reveal any additional loci linked to survival or MTTD traits (Figure 6C, D).

Figure 6
figure 6

Composite interval mapping indicates interactions of QTLs on chromosomes 5 and 19. The influence of the markers rs13483633 on chromosome 19 (A, B) and rs13478587 on chromosome 5 (C, D) was factored out and the residual LRS were calculated. Two different traits, survival day 8 (A, C) and mean time of death (B, D) were tested by using the web tool of GeneNetwork. Interval mapping is shown across the entire genome as blue line representing the LRS scores; green line: DBA/2J alleles increase trait values; red line: C57BL/6J alleles increase trait values. Upper x-axis shows mouse chromosomes, lower x-axis shows physical map in mega bases for each chromosome, y-axis represents LRS score. The horizontal lines mark the genome-wide significant thresholds at p<0.05 (red line) and suggestive thresholds at p<0.37 (gray line).

Finally, we analyzed the joint additive effects of QTLs described above as well as possible epistatic interactions among these QTLs. When using a full two-locus model (Trait Variance = Q1 + Q2 + Q1 x Q2), we found two potential interactions between Qivr5 and a new locus on chromosome 9 as well as between Qivr5 and Qivr19, respectively, for the survival trait at day 8 (Figure 7A). The inclusion of an interaction term increases the LRS by 23 (Additional file 6: Figure S5). For both Qivr5 - chromosome 9 interaction (Figure 7B) and Qivr5 – Qivr19 interaction (Figure 7C), the allele combination DBA/2J / DBA/2J showed highest and the combination C57BL/6J / C57BL/6J lowest survival scores. All differences were highly significant.

Figure 7
figure 7

Interactions are observed between several QTLs. The graph displays pair-scan results for the trait percent survival at day 8 p.i. (ID13002) (A). The upper left half of the plot highlights any epistatic interactions, the lower right half provides a summary of LRS of the full model, representing cumulative effects of linear and non-linear terms based on the model ‘Trait Variance = Q1 + Q2 + Q1 x Q2 ’. Subsequently, BXD lines were grouped according to their allele combinations and trait values were compared between groups (B, C). The largest differences are indicated by asterisks. They are highly significant between the respective groups (p < 0.0001). For the various loci, trait values for the following markers were used: Qivr5: rs13478587, chromosome 9: rs6191976, Qivr19: rs13483633. Number of lines and p-values for Qivr5 – chromosome 9: D2/D2: n = 17, D2/B6: n = 17, B6/D2: n = 6, B6/B6 n = 12; D2/D2 vs B6/B6: p < 0.0001, D2/B6 vs B6/B6: p = 0.003. Number of lines and p-values for Qivr5 – Qivr19: D2/D2: n = 18, D2/B6: n = 16, B6/D2: n = 4, B6/B6: n = 14; p-values: D2/D2 vs B6/B6: p < 0.0001, D2/B6 vs B6/B6: p = 0.0008. Abbreviations: D2: DBA/2J, B6: C57BL/6J.

Analysis of QTL regions identified several candidate genes that may contribute to host susceptibility or resistance

In total, the mapping studies revealed five QTLs on chromosomes 2, 5, distal 16, 17 and 19 (Qivr2-2, Qivr5, Qivr16, Qivr17-2 and Qivr19) that did merit further analysis because they were consistently observed in at least two traits and exerted an effect on at least two different days p.i.

For the candidate gene searches, we defined the critical regions of the QTLs manually by considering peak height, its shoulder and bootstraps as guidance (Figure 8). We then investigated these Qivr intervals for potential candidate genes that may be causal for the studied traits (Additional file 7: Figure S6). We first used the QTLminer tool in GeneNetwork [25] to identify all genes in the QTL intervals and to obtain associated GO terms. Next, we select all annotated genes within the QTL intervals that were expressed between day 1 and day 60 after infection of C57BL/6J mice with PR8M virus (the latter data set is derived from a separate study (Pommerenke et al., PLoS ONE, in press). Subsequently, these genes were further characterized for the following attributes: genes that were up- or up-regulated in infected lungs by at least 1.5-fold in C57BL/6J after PR8M infection, genes carrying an insertion or deletion or a non-synonymous nucleotide change in the open reading frame, genes exhibiting a cis-expression QTL in the lung [26], and genes that were differentially expressed by at least 1.5-fold in C57BL/6J and DBA/2J mice in infected lungs between day 1 and day 8 p.i. [23] and Pommerenke et al., PLoS ONE, in press). The strategy of the QTL mining is shown as flow chart in the Additional file 7: Figure S6, and the final results are presented in Table 1. Furthermore, the attributes of all genes located in the QTL regions are listed in detail in the Additional file 8: Table S2). Using these combined attributes, we identified 31 genes as the most likely candidates to regulate the traits controlled by Qivr2-2, Qivr5, Qivr16, Qivr17-2 and Qivr19 (listed in Table 2). These genes were further evaluated based on their known function from the literature and phenotypes in knock-out mouse mutants (see discussion).

Figure 8
figure 8

Detailed maps of QTLs contributing to host susceptibility and resistance to PR8 influenza A infection. Interval mapping was performed for chromosomes 2 (A), 5 (B), 16 (C), 17 (D) and 19 (E). The critical intervals were selected based on the peak shapes and the bootstrap signals: Qivr2-1 (56–68 Mb), Qivr5 (140–153 Mb), Qivr16 (64–78 Mb), Qivr17-2 (30–44 Mb) and Qivr19 (37–45 Mb). The Qivr intervals are delineated by red arrows. The yellow bars represent the frequencies of peak LRS using bootstrap analysis. The multicolored chatters along the top of the graph are hyperlinks to sites with additional genetic and sequence information. The orange chatter along the x-axis indicates the density of SNPs present in the BXD strains. Interval mapping is shown across the entire genome as blue line representing the LRS scores; green line: DBA/2J alleles increase trait values; red line: C57BL/6J alleles increase trait values. Upper x-axis shows mouse chromosomes, lower x-axis shows physical map in mega bases for each chromosome, y-axis represents LRS score. The horizontal lines mark the genome-wide significant thresholds at p<0.05 (red line) and suggestive thresholds at p<0.37 (gray line).

Table 1 Candidate genes in mapped QTL intervals
Table 2 GO-terms and functions observed in knock-out mice of 31 potential candidates from QTL intervals

Discussion

DBA/2J and C57BL/6J mice have been shown previously to differ largely in their susceptibility to H1N1 (PR8M) influenza A virus [9, 18]. Here, we expanded these studies and utilized the BXD recombinant inbred set of mouse strains to map the genomic regions that are responsible for differences in these two mouse strains. We monitored three phenotypic traits, body weight over time, survival over time and mean time to death to identify quantitative trait for influenza resistance. Two significant QTLs, Qivr5 and Qivr19, were found on chromosomes 5 and 19, respectively. Furthermore several suggestive QTLs, Qivr2-2, Qivr16 and Qivr17-2 were observed in at least two traits and at two days on chromosomes 2, 16 and 17, respectively. Composite mapping revealed an additional almost significant QTL at distal chromosome 16, Qivr16.

A similar analysis for host resistance to influenza has been performed previously after infecting 66 BXD strains with H5N1 influenza virus. This study reported three significant QTLs on chromosomes 2, 7, and 17 [12]. Thus, none of these significant QTLs overlaps with the QTLs identified in our analysis. Five of the strains that were resistant (all infected mice survived) in our study were also resistant in the study of [12] where a total of 14 strains were found to be resistant. Five strains that were highly susceptible in our study (100% of infected mice died) were also highly susceptible in the study by [12] of a total of 26 susceptible strains. Furthermore, five strains that were resistant in our study were susceptible in the study by [12]. Thus, there is also not much overlap between the two studies with respect of susceptible and resistant strains. The differences between the two studies are most likely explained by the use of two different influenza virus subtypes. The H1N1 virus from our study represents a subtype with a monobasic hemagglutinin (HA) cleavage site, whereas the H5N1 which was used in the study by Boon et al. is a subtype with a polybasic HA cleavage site. The cellular tropism of these two subtypes for virus replication and processing is quite different, because monobasic viruses are dependent on cell-specific proteases for the processing of the HA whereas polybasic subtypes can be processed by more ubiquitously expressed host proteases, e.g.[27–32]. Therefore, the contribution of host factors to susceptibility may be different between H1 and H5 containing virus subtypes.

Another study described the genetic mapping of susceptibility and resistance factors after infecting a panel of 29 AxB / BxA congenic strains with a mouse-adapted H3N2 influenza virus [33]. The AxB / BxA congenic strains were generated from a cross of susceptible A/J and resistant C57BL/6J parental mouse strains. The authors found three major QTLs on chromosomes 2, 6 and 17. The QTL on chromosome 17 overlaps with the Qivr17-2 locus which we found in our study. Furthermore, the candidate gene Pla2g7 that was identified in their study was also detected as candidate gene in our analysis (see below).

The influence of genetic factors determining the host response to H1N1 influenza virus infections was also examined in mice of the pre-Collaborative Cross collection [34]. In this study, gene expression levels in extreme responders were used to identify expression QTLs (eQTL). One gene that exhibited a cis-eQTL, Sik1 (salt inducible kinase 1), was located in the Qivr17-2 interval from our study, and we also identified it as potential quantitative trait gene (Table 2). This gene is associated with the GO terms ‘negative regulation of transcription from RNA polymerase II promoter, regulation of cell differentiation, and protein kinase cascade’. However, no specific infection-related functions have been yet described for this gene.

One of the most interesting findings in our study was the time-dependent effect of QTLs which we observed in the body weight and survival traits. The peak QTLs for the two significant QTLs, Qivr5 and Qivr19, were found at different times p.i., day 6 and day 8, respectively. In addition, the effects of both QTLs were not only evident at the times p.i. where they exerted the significant peak QTL signals but also several days before and after the peak. Furthermore, for the suggestive QTLs, also time-dependent effects were observed. These results suggest that the causal genes underlying different QTLs act at different time points of the host defense.

Most interestingly, Qivr5 as well as Qivr19 represent a positive influence on body weight, survival and MTTD from the DBA/2J haplotype, the susceptible strain. These findings indicate that genomic regions from the susceptible parent are able to increase resistance when combined with the resistant parental genome. We are now analyzing several BXD strains that were more resistant than the parental C57BL/6J mice in more detail. One possible mechanism to explain such an effect may be that an activator (secreted ligand or transcription factor) is expressed in susceptible DBA/2J mice but the corresponding target (receptor or regulated gene) is mutated. On the other hand, in C57BL/6J mice, the target but not the activator may be mutated. If the wild type alleles are now coming together in a BXD strain, the functional activator finds its functional target and thereby an increased resistance state is achieved.

Both composite and interaction mapping revealed many genetic interactions between C57BL/6J and DBA/2 J alleles. Thus, many genomic regions from the parental strains are able to contribute to the host response and this effect depends strongly on the allele combinations in the respective QTLs. These observations may be studied further in double congenic mouse lines.

We subsequently analyzed the five QTL intervals, Qivr2-2, Qivr5, Qivr16, Qivr17-2 and Qivr19 in more detail to identify genes that may be causal for resistance or susceptibility. In total, 830 genes are located in these intervals. We narrowed down the total number of genes to 31 candidates (Table 2) by using additional information, such as temporal expression after PRM8 infection (Pommerenke et al., PLoS ONE, in press), cis-eQTLs in non-infected lungs [26], differences in expression between DBA/2J and C57BL/6J [23], and sequence variants in the coding regions.

Qivr5 contains the candidate gene Eif2ak1 (eukaryotic translation initiation factor 2 alpha kinase 1) that is a member of eIF2alpha kinases which have been associated with anti-viral host responses [35]. Boon et al. described another eIF2alpha kinase, Eif2ak2 / Pkr (eukaryotic translation initiation factor 2-alpha kinase 2), in the Qirv17 locus after infection with influenza H5N1 [12]. Eif2ak2 plays a critical role in modulating immunoglobulin expression during RSV infection. In addition Eif2ak2 knock-out mouse mutants are more susceptible to influenza infections [36, 37]. We have initiated the generation of a congenic mouse lines for the chromosome 5 interval to verify and further characterize the effect of this region for resistance to influenza infection.

Qivr2-2 contains two candidate genes, Itgb6 (integrin beta 6) and Ifih1 (interferon induced with helicase C domain 1), with known functions in the host defense to viral infections. Itgb6 mouse knock-out mutants exhibit severe pneumonia and an increase in granulocyte recruitment to the lung [38]. The protease-activated receptor 1-mediated enhancement of Itgb6- dependent TGF-beta activation has been proposed to represent one mechanism by which activation of the coagulation cascade contributes to the development of acute lung injury [39]. The Ifih1 gene is also known as MDA5 (Melanoma Differentiation-Associated protein 5). IFIH1 is part of the RIG-I-like receptor (RLR) family, which function as pattern recognition receptors and are activated upon binding of virus dsRNA [40]. IFIH1 functions as cytosolic receptor that leads to the selective activation of type I IFN genes and is indispensable for sustained expression of IFN in response to paramyxovirus infection [41, 42]. Ifih1 mutant knock-out mice exhibit an impaired response to different viral pathogens [43, 44].

Qivr16 contains two potential genes with known functions in the host defense and lung function, Robo1 (roundabout homolog 1 (Drosophila)) and Nrip1 (nuclear receptor interacting protein 1). DBA/2J mice carry a frame shift mutation in the Robo1 gene which might lead to an impaired function of the encoded protein. Robo1 has been described to be involved in guidance and migration of axons, myoblasts, and leukocytes in vertebrates (e.g.[45–47]) but is also expressed in the developing lung [48]. Robo1 knock-out mutants exhibit a delayed lung maturation and bronchial hyperplasia. The latter results suggest that Robo1 may be involved in maintaining proper lung function and it may become essential when lung epithelium is destroyed during an influenza infection. Nrip1/Rip140 functions as a co-activator for cytokine gene promoter activity via direct protein-protein interactions with the NFkappaB subunit RelA and histone acetylase cAMP-responsive element binding protein (CREB)-binding protein (CBP) [49]. It is involved in modulating pro-inflammatory responses in macrophages [50].

Qivr17-2 represents a positive influence of the C57BL/6 J genotype on body weight, survival and MTTD. This QTL is located in a gene-rich region which carries many genes that are involved in the host immune response, in particular the H2 histocompatibility genes which are involved in antigen presentation [51]. Therefore, many candidate genes are found in the Qivr17-2 region. The Lst1 (leukocyte specific transcript 1) gene is of special interest because the DBA/2J allele mice carries a single nucleotide deletion in the first exon resulting in a frame shift of the open reading frame. This mutation most likely results in a non-functional Lst1 protein in DBA/2 J mice. We confirmed the presence of the deletion by sequencing the parental DBA/2J and some BXD strains carrying the DBA/2J allele. The wild type allele was confirmed in C57BL/6J mice and in some BXD strains carrying the C57BL/6J allele. In humans, LST1 plays a role in the regulation of the immune response to inflammatory diseases such as rheumatoid arthritis, microbial infection or Rubella vaccine-induced immunity [52–55]. Also, Lst1 is up-regulated after influenza A infection in C57BL/6J mice starting at day 2 and exhibits a strong peak of expression at day 8 p.i. Pommerenke et al., 2012 (Pommerenke, C., E. Wilk, B. Srivastava, A. Schulze, N. Novoselova, R. Geffers, and K. Schughart. 2012. Global transcriptome analysis in influenza-infected mouse lungs reveals the kinetics of innate and adaptive host immune responses. PLoS ONE. 7:e41169.). Thus, the expression profile and known functions of Lst1 fit well with a possible critical role for the host defense to influenza A virus. We initiated the generation of knock-out mice to evaluate the role of Lst1 in more detail. In addition, a second, most interesting candidate, Pla2g7 (phospholipase A2, group VII (platelet-activating factor acetylhydrolase, plasma)) was identified in the Qivr17-2 interval. In humans, increased activities of certain variants of PLA2G7 were associated with early coronary atherosclerosis and with endothelial dysfunction, but the gene may also exert an anti-inflammatory function [56–60]. The Pla2g7 gene was also identified as a potential candidate gene for susceptibility against infections with H3N2 influenza virus [33]. Pla2g7 expression levels in susceptible A/J mice were higher than in resistant C57BL/6J mice after infection with H3N2 virus [33]. We also showed previously that Pla2g7 exhibits a cis-eQTL between C57BL/6J and DBA/2J in non-infected lungs where the DBA/2J allele shows high levels of expression [26]. Tnfrsf21 which was identified by [33] as potential candidate of Qivr17-2 also exhibits a cis-eQTL in non-infected BXD mice [26] but was not found to be regulated in C57BL/6 mice after infection (data not published). Tapbp (TAP binding protein) plays a major role in the antigen processing and MHC class I presentation by stabilizing the TAP peptide transporter, e.g.[61–65]. Also, Tap2 (transporter 2, ATP-binding cassette, sub-family B (MDR/TAP)) gene is involved in antigen processing and presentation [63, 66]. Gtf2h4 (general transcription factor II H, polypeptide 4) encodes a general transcription factor. Recruitment and activation of Gtf2h4 represents a rate-limiting step for the emergence of HIV from latency and sequence variants have been associated with multiple sclerosis [67–69].

Within the Qivr19 interval, only one gene, Hps1 (Hermansky-Pudlak syndrome 1 homolog (human)), has been associated with the host responses to infection. Mice carrying a natural mutation in the Hps1 gene showed an increased inflammatory response in alveolar macrophages after intranasal challenge with LPS [70].

The GeneNetwork database allows searching for other phenotypic traits that exhibit a genome-wide significant (LRS ≥ 18) within the Qivr intervals identified by our study. Two phenotypic traits, related to neuronal responses (trait ID 11285) and body weight changes (trait ID 12838), are located to the Qivr16 locus. Also, the Qivr17-2 interval contained significant QTLs for other traits. Two traits are related to host infectious diseases, ‘Ectromelia virus survival’ (ID 12672) and ‘Chlamydia psittaci (6BC) infection response’ (ID 11025) and four traits are associated with immune cell responses (ID 10201, 10466, 10238, 10236). In addition two traits described seizure responses (ID 10388, 10507), and one trait has not been disclosed yet (ID 13920). Within the early time chromosome 10 interval, three other traits exhibit their most significant QTLs: ‘3a,5a-THDOC in blood plasma 3 days after cycle 5 of chronic intermittent air vapor’ (ID13027) and two non-disclosed traits. The first trait may relate to stress responses in the central nervous system (ID 13292 and 13846).

Conclusions

The mapping of resistance and susceptibility loci in the BXD population revealed several new QTLs and potential gene candidates that may be critical for the host defense against influenza A virus infection. Body weight and survival QTLs showed a time-dependent profile indicating that the genetic factors in these QTLs are important for the host response in a temporal dynamic fashion. Five QTL regions were examined in detail, and we identified several possible candidate genes that may be critical for the host response to influenza A infections in humans.

Methods

Animals

The mouse inbred strains C57BL/6J, DBA/2J and B6D2F1 were delivered from Janvier, France. Recombinant inbred mouse strains BXD were purchased from three different sources: The Jackson Laboratory, the University of Tennessee Health Science Center (Memphis, TN) and from Harlan, The Netherlands. For the analysis, mice were transferred to the animal facility in Braunschweig and adapted for at least two weeks to the new environment before starting experiments. Animals were maintained under specific pathogen free conditions. All experiments in mice were approved by an external committee and according to the national guidelines of the animal welfare law in Germany (‘Tierschutzgesetz in der Fassung der Bekanntmachung vom 18. Mai 2006 (BGBl. I S. 1206, 1313), das zuletzt durch Artikel 20 des Gesetzes vom 9. Dezember 2010 (BGBl. I S. 1934) geändert worden ist.’). The protocol used in these experiments has been reviewed by an ethics committee and approved by the ‘Niedersächsisches Landesamt für Verbraucherschutz und Lebensmittelsicherheit, Oldenburg, Germany’, according to the German animal welfare law (Permit Numbers: 33.42502/04-108/06, 33.9.42502-04-051/09).

Virus and infection of mice

The mouse-adapted influenza strain A/Puerto Rico/8/34 (H1N1; PR8M, [18] and references therein) was used for all infection studies. Stocks were prepared by infection of 10-day-old embryonated chicken eggs. After mice were anesthetized by intra-peritoneal injection of Ketamin-Xylazine solution in sterile NaCl (50 mg/ml Ketamine, Invesa Arzneimittel GmbH, Freiburg; 2% Xylazine, Bayer Health-Care, Leverkusen) with a dose adjusted to the individual body weight, mice were infected intranasally with 2 × 103 FFU of PR8M in 20 μl of sterile phosphate-buffered saline. Mice were assayed daily for body weight (determined as % of initial weight at day 0) and mortality during 13 days p.i. We used death as the end point for survival. Mice were sacrificed if body weight loss exceeded 25%. It should be noted that for mice that did not show any signs of body weight loss over the entire time period after infection, we do not have additional parameters to verify that they have indeed been infected. However, these cases were very few. In addition, we have ample experience with this infection method and the failure rate, for example with the DBA/2J strain, is less that 5%.

Data handling and statistical analysis

In total, 283 BXD mice and 127 mice from the parental strains or F1 generation were used for the infection experiments. In total 53 BXD strains were infected with an average of 5 mice per strain. We performed all primary calculations using simple features and functions of Microsoft Excel. Three sets of analysis were performed for the following variables: (1) body weight in percentage from starting weight (day 0) using the strain medians to exclude outliers; (2) survival by calculating the survival rate of each strain from day 7 to 13, (3) mean time to death in days. For statistical analyses, tests and visualization we used R, a free software environment for statistical computing and graphics (http://www.R-project.org). In order to test for batch effects or other co-factors, we visualized the data using multidimensional scaling based on the Sammon mapping method [71]. No clusters with respect to any of the co-factors age, weight at day 0, experimenter, time of infection, or source of mice could be found in the visualizations.

QTL mapping

QTL mapping was performed using the web-based complex trait analysis available on the GeneNetwork website (http://www.genenetwork.org) and the mapping module to analyze phenotypes in context of mouse genotypic differences. Interval mapping evaluates potential QTLs at regular intervals and estimates the significance at each location with a graphical representation of the likelihood ratio statistics (LRS) using 2000 permutation tests [19, 22]. LRS values may be converted to LOD scores by dividing by 4.61. For the two locus model the following equation was used: Var = Q1 + Q2 + Q1xQ2 + e, where Var = the between-strain mean variance in the trait, Q1 and Q2 are makers tightly linked to the loci, Q1XQ2 is the ‘additive-by-additive’ epistatic interaction term, and e is the residual error. The original data sets can be obtained at http://www.genenetwork.org with the following identification numbers: body weight: 13005 to 13017; survival: 13000 to 13004 and mean time to death: 12996. We performed full genome scans for epistatic interactions using the pair-scan module that is implemented in GeneNetwork. This module exploits the direct global optimization algorithm developed by [72]. The code compares the fit (as measured by LRS scores) for a purely additive model, a purely epistatic model, and the full model. The code also implements a permutation test (n = 500) and this enabled us to estimate the empirical p value of the alternative models.

Candidate gene discovery

The QTL region analysis was initially performed using the QTLminer which has been implemented in GeneNetwork [25]. By using the automatic function of GeneNetwork we identified significant cis-QTLs with LRS higher than 18 at a genome-wide p-value of < 0.05. Additionally the genes mapped within the analyzed QTLs were surveyed by the National Center for Biotechnology Information (NCBI) Entrez Gene website (http://www.ncbi.nlm.nih.gov/sites/entrez?db=gene) and the Jackson Laboratory's MGI Mouse Genome Database project (http://www.informatics.jax.org/) to identify potential candidate genes. The GeneRIF database (http://www.ncbi.nlm.nih.gov/projects/GeneRIF/GeneRIFhelp.html) was used as a primary source to search for known gene functions and corresponding citations.

All sequence variants between B6 and D2 parental genomes (SNPs, indels) were extracted by using a comparative analysis that relies on approximately 100x whole genome shot gun of DBA/2J [73]. All of these sequence data are available at http://ucscbrowser.genenetwork.org/, the GeneNetwork Variant Browser (http://www.genenetwork.org/webqtl/main.py?FormID=snpBrowser), and the NCBI Short Read Archive (18 files total, e.g., SRX037575, SRX013980, SRX013299, SRX012582, SRX012581, SRX012580); http://www.biomedcentral.com/1471-2105/11/S4/O7.

Authors’ information

Tatiana Nedelko and Heike Kollmus contributed equally as first authors.

Robert W. Williams and Klaus Schughart contributed equally as senior authors.

References

  1. Johnson NP, Mueller J: Updating the accounts: global mortality of the 1918–1920 "Spanish" influenza pandemic. Bull Hist Med. 2002, 76 (1): 105-115. 10.1353/bhm.2002.0022.

    Article  PubMed  Google Scholar 

  2. Fauci AS: Seasonal and pandemic influenza preparedness: science and countermeasures. J Infect Dis. 2006, 194 (Suppl 2): S73-76.

    Article  PubMed  Google Scholar 

  3. Scriven J, McEwen R, Mistry S, Green C, Osman H, Bailey M, Ellis C: Swine flu: a Birmingham experience. Clin Med. 2009, 9 (6): 534-538.

    Article  Google Scholar 

  4. Yates L, Pierce M, Stephens S, Mill A, Spark P, Kurinczuk J, Valappil M, Brocklehurst P, Thomas S, Knight M: Influenza A/H1N1v in pregnancy: an investigation of the characteristics and management of affected women and the relationship to pregnancy outcomes for mother and infant. Health Technol Assess. 2010, 14 (34): 109-182.

    Article  CAS  PubMed  Google Scholar 

  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. 10.1086/524064.

    Article  PubMed  Google Scholar 

  6. Gottfredsson M, Halldorsson BV, Jonsson S, Kristjansson M, Kristjansson K, Kristinsson KG, Love A, Blondal T, Viboud C, Thorvaldsson S, 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-1308. 10.1073/pnas.0707659105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Horby P, Sudoyo H, Viprakasit V, Fox A, Thai PQ, Yu H, Davila S, Hibberd M, Dunstan SJ, Monteerarat Y, 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.

    Article  Google Scholar 

  8. Everitt AR, Clare S, Pertel T, John SP, Wash RS, Smith SE, Chin CR, Feeley EM, Sims JS, Adams DJ, et al: IFITM3 restricts the morbidity and mortality associated with influenza. Nature. 2012, 484 (7395): 519-523. 10.1038/nature10921.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Srivastava B, Blazejewska P, Hessmann M, Bruder D, Geffers R, Mauel S, Gruber AD, Schughart K: Host genetic background strongly influences the response to influenza a virus infections. PLoS ONE. 2009, 4 (3): e4857-10.1371/journal.pone.0004857.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Trammell RA, Toth LA: Genetic susceptibility and resistance to influenza infection and disease in humans and mice. Expert Rev Mol Diagn. 2008, 8 (4): 515-529. 10.1586/14737159.8.4.515.

    Article  CAS  PubMed  Google Scholar 

  11. Ding M, Lu L, Toth LA: Gene expression in lung and basal forebrain during influenza infection in mice. Genes Brain Behav. 2008, 7 (2): 173-183. 10.1111/j.1601-183X.2007.00335.x.

    Article  CAS  PubMed  Google Scholar 

  12. Boon AC, DeBeauchamp J, Hollmann A, Luke J, Kotb M, Rowe S, Finkelstein D, Neale G, Lu L, Williams RW, et al: Host genetic variation affects resistance to infection with a highly pathogenic H5N1 influenza A virus in mice. J Virol. 2009, 83 (20): 10417-10426. 10.1128/JVI.00514-09.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Boon AC, Debeauchamp J, Krauss S, Rubrum A, Webb AD, Webster RG, McElhaney J, Webby RJ: Cross-reactive neutralizing antibodies directed against pandemic H1N1 2009 virus are protective in a highly sensitive DBA/2 influenza mouse model. J Virol. 2010, 84 (15): 7662-7667. 10.1128/JVI.02444-09.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Otte A, Sauter M, Alleva L, Baumgarte S, Klingel K, Gabriel G: Differential host determinants contribute to the pathogenesis of 2009 pandemic H1N1 and human H5N1 influenza A viruses in experimental mouse models. Am J Pathol. 2011, 179 (1): 230-239. 10.1016/j.ajpath.2011.03.041.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. 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-11.

    Article  Google Scholar 

  16. Trammell RA: Liberati TA. 2011, Toth LA: Host genetic background and the innate inflammatory response of lung to influenza virus. Microbes Infect

    Google Scholar 

  17. Pica N, Iyer A, Ramos I, Bouvier NM, Fernandez-Sesma A, Garcia-Sastre A, Lowen AC, Palese P, Steel J: The DBA.2 mouse is susceptible to disease following infection with a broad, but limited, range of influenza A and B viruses. J Virol. 2011, 85 (23): 5-12829.

    Article  Google Scholar 

  18. Blazejewska P, Koscinski L, Viegas N, Anhlan D, Ludwig S, Schughart K: Pathogenicity of different PR8 influenza A virus variants in mice is determined by both viral and host factors. Virology. 2011, 412 (1): 36-45. 10.1016/j.virol.2010.12.047.

    Article  CAS  PubMed  Google Scholar 

  19. Peirce JL, Lu L, Gu J, Silver LM, Williams RW: A new set of BXD recombinant inbred lines from advanced intercross populations in mice. BMC Genet. 2004, 5: 7-

    Article  PubMed Central  PubMed  Google Scholar 

  20. Taylor BA, Wnek C, Kotlus BS, Roeme RN, MacTaggart T, Phillips SJ: Genotyping new BXD recombinant inbred mouse strains and comparison of BXD and consensus maps. Mamm Genome. 1999, 10 (4): 335-348. 10.1007/s003359900998.

    Article  CAS  PubMed  Google Scholar 

  21. Williams RW, Gu J, Qi S, Lu L: The genetic structure of recombinant inbred mice: high-resolution consensus maps for complex trait analysis. Genome Biol. 2001, 2 (11): RESEARCH0046-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. 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-486. 10.1038/nn0504-485.

    Article  CAS  PubMed  Google Scholar 

  23. Alberts R, Srivastava B, Wu H, Viegas N, Geffers R, Klawonn F, Novoselova N, Do Valle TZ, Panthier JJ, Schughart K: Gene expression changes in the host response between resistant and susceptible inbred mouse strains after influenza A infection. Microbes Infect. 2010, 12 (4): 309-318. 10.1016/j.micinf.2010.01.008.

    Article  CAS  PubMed  Google Scholar 

  24. Ziebarth JD, Cook MN, Wang X, Williams RW, Lu L, Cui Y: Treatment- and population-dependent activity patterns of behavioral and expression QTLs. PLoS ONE. 2012, 7 (2): e31805-10.1371/journal.pone.0031805.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Alberts R, Schughart K: QTLminer: identifying genes regulating quantitative traits. BMC Bioinformatics. 2010, 11: 516-10.1186/1471-2105-11-516.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Alberts R, Lu L, Williams RW, Schughart K: Genome-wide analysis of the mouse lung transcriptome reveals novel molecular gene interaction networks and cell-specific expression signatures. Respir Res. 2011, 12: 61-10.1186/1465-9921-12-61.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Zhirnov OP, Matrosovich TY, Matrosovich MN, Klenk HD: Aprotinin, a protease inhibitor, suppresses proteolytic activation of pandemic H1N1v influenza virus. Antivir Chem Chemother. 2011, 21 (4): 169-174. 10.3851/IMP1715.

    Article  CAS  PubMed  Google Scholar 

  28. Okumura Y, Takahashi E, Yano M, Ohuchi M, Daidoji T, Nakaya T, Bottcher E, Garten W, Klenk HD, Kido H: Novel type II transmembrane serine proteases, MSPL and TMPRSS13, proteolytically activate membrane fusion activity of hemagglutinin of highly pathogenic avian influenza viruses and induce their multicycle replication. J Viro. 2010, 84 (10): 5089-5096. 10.1128/JVI.02605-09.

    Article  CAS  Google Scholar 

  29. Bottcher-Friebertshauser E, Freuer C, Sielaff F, Schmidt S, Eickmann M, Uhlendorff J, Steinmetzer T, Klenk HD, Garten W: Cleavage of influenza virus hemagglutinin by airway proteases TMPRSS2 and HAT differs in subcellular localization and susceptibility to protease inhibitors. J Virol. 2010, 84 (11): 5605-5614. 10.1128/JVI.00140-10.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Bertram S, Glowacka I, Blazejewska P, Soilleux E, Allen P, Danisch S, Steffen I, Choi SY, Park Y, Schneider H, et al: TMPRSS2 and TMPRSS4 facilitate trypsin-independent spread of influenza virus in Caco-2 cells. J Virol. 2010, 84 (19): 10016-10025. 10.1128/JVI.00239-10.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Chaipan C, Kobasa D, Bertram S, Glowacka I, Steffen I, Tsegaye TS, Takeda M, Bugge TH, Kim S, Park Y, et al: Proteolytic activation of the 1918 influenza virus hemagglutinin. J Virol. 2009, 83 (7): 3200-3211. 10.1128/JVI.02205-08.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Steinhauer DA: Role of hemagglutinin cleavage for the pathogenicity of influenza virus. Virology. 1999, 258 (1): 1-20. 10.1006/viro.1999.9716.

    Article  CAS  PubMed  Google Scholar 

  33. Boivin GA, Pothlichet J, Skamene E, Brown EG, Loredo-Osti JC, Sladek R, Vidal SM: 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

    Google Scholar 

  34. Bottomly D, Ferris MT, Aicher LD, Rosenzweig E, Whitmore A, Aylor DL, Haagmans BL, Gralinski LE, Bradel-Tretheway BG, Bryan JT, et al: Expression quantitative trait Loci for extreme host response to influenza a in pre-collaborative cross mice. G3 (Bethesda). 2012, 2 (2): 213-221.

    Article  Google Scholar 

  35. Krishnamoorthy J, Mounir Z, Raven JF, Koromilas AE: The eIF2alpha kinases inhibit vesicular stomatitis virus replication independently of eIF2alpha phosphorylation. Cell Cycle. 2008, 7 (15): 2346-2351.

    Article  CAS  PubMed  Google Scholar 

  36. Thakur SA, Zalinger Z, Johnson TR, Imani F: PKR is a novel mediator of CD40 signaling and plays a critical role in modulating immunoglobulin expression during RSV infection. Clin Vaccine Immunol. 2011, 18 (12): 2060-2066. 10.1128/CVI.05476-11.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Bergmann M, Garcia-Sastre A, Carnero E, Pehamberger H, Wolff K, Palese P, Muster T: Influenza virus NS1 protein counteracts PKR-mediated inhibition of replication. J Virol. 2000, 74 (13): 6203-6206. 10.1128/JVI.74.13.6203-6206.2000.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Ludlow A, Yee KO, Lipman R, Bronson R, Weinreb P, Huang X, Sheppard D, Lawler J: Characterization of integrin beta6 and thrombospondin-1 double-null mice. J Cell Mol Med. 2005, 9 (2): 421-437. 10.1111/j.1582-4934.2005.tb00367.x.

    Article  CAS  PubMed  Google Scholar 

  39. Hogmalm A, Sheppard D, Lappalainen U, Bry K: beta6 Integrin subunit deficiency alleviates lung injury in a mouse model of bronchopulmonary dysplasia. Am J Respir Cell Mol Biol. 2010, 43 (1): 88-98. 10.1165/rcmb.2008-0480OC.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Yoneyama M, Kikuchi M, Matsumoto K, Imaizumi T, Miyagishi M, Taira K, Foy E, Loo YM, Gale M, Akira S, et al: Shared and unique functions of the DExD/H-box helicases RIG-I, MDA5, and LGP2 in antiviral innate immunity. J Immunol. 2005, 175 (5): 2851-2858.

    Article  CAS  PubMed  Google Scholar 

  41. Choi MK, Wang Z, Ban T, Yanai H, Lu Y, Koshiba R, Nakaima Y, Hangai S, Savitsky D, Nakasato M, et al: A selective contribution of the RIG-I-like receptor pathway to type I interferon responses activated by cytosolic DNA. Proc Natl Acad Sci U S A. 2009, 106 (42): 17870-17875. 10.1073/pnas.0909545106.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Gitlin L, Benoit L, Song C, Cella M, Gilfillan S, Holtzman MJ, Colonna M: Melanoma differentiation-associated gene 5 (MDA5) is involved in the innate immune response to Paramyxoviridae infection in vivo. PLoS Pathog. 2010, 6 (1): e1000734-10.1371/journal.ppat.1000734.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Gitlin L, Barchet W, Gilfillan S, Cella M, Beutler B, Flavell RA, Diamond MS, Colonna M: Essential role of mda-5 in type I IFN responses to polyriboinosinic:polyribocytidylic acid and encephalomyocarditis picornavirus. Proc Natl Acad Sci U S A. 2006, 103 (22): 8459-8464. 10.1073/pnas.0603082103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Kato H, Takeuchi O, Sato S, Yoneyama M, Yamamoto M, Matsui K, Uematsu S, Jung A, Kawai T, Ishii KJ, et al: Differential roles of MDA5 and RIG-I helicases in the recognition of RNA viruses. Nature. 2006, 441 (7089): 101-105. 10.1038/nature04734.

    Article  CAS  PubMed  Google Scholar 

  45. Andrews W, Barber M, Hernadez-Miranda LR, Xian J, Rakic S, Sundaresan V, Rabbitts TH, Pannell R, Rabbitts P, Thompson H, et al: The role of Slit-Robo signaling in the generation, migration and morphological differentiation of cortical interneurons. Dev Biol. 2008, 313 (2): 648-658. 10.1016/j.ydbio.2007.10.052.

    Article  CAS  PubMed  Google Scholar 

  46. Long H, Sabatier C, Ma L, Plump A, Yuan W, Ornitz DM, Tamada A, Murakami F, Goodman CS, Tessier-Lavigne M: Conserved roles for Slit and Robo proteins in midline commissural axon guidance. Neuron. 2004, 42 (2): 213-223. 10.1016/S0896-6273(04)00179-5.

    Article  CAS  PubMed  Google Scholar 

  47. Xian J, Clark KJ, Fordham R, Pannell R, Rabbitts TH, Rabbitts PH: Inadequate lung development and bronchial hyperplasia in mice with a targeted deletion in the Dutt1/Robo1 gene. Proc Natl Acad Sci U S A. 2001, 98 (26): 15062-15066. 10.1073/pnas.251407098.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Clark K, Hammond E, Rabbitts P: Temporal and spatial expression of two isoforms of the Dutt1/Robo1 gene in mouse development. FEBS Lett. 2002, 523 (1–3): 12-16.

    Article  CAS  PubMed  Google Scholar 

  49. Zschiedrich I, Hardeland U, Krones-Herzig A, Berriel Diaz M, Vegiopoulos A, Muggenburg J, Sombroek D, Hofmann TG, Zawatzky R, Yu X, et al: Coactivator function of RIP140 for NFkappaB/RelA-dependent cytokine gene expression. Blood. 2008, 112 (2): 264-276. 10.1182/blood-2007-11-121699.

    Article  CAS  PubMed  Google Scholar 

  50. Ho PC, Chang KC, Chuang YS, Wei LN: Cholesterol regulation of receptor-interacting protein 140 via microRNA-33 in inflammatory cytokine production. FASEB J. 2011, 25 (5): 1758-1766. 10.1096/fj.10-179267.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Janeway CA, Travers P, Walport MJ, Shlomchik MJ: Immunobiology: The Immune System in Health and Disease. 2005, New York: Garland Science

    Google Scholar 

  52. Ovsyannikova IG, Vierkant RA, Pankratz VS, Jacobson RM, Poland GA: Extended LTA, TNF, LST1 and HLA gene haplotypes and their association with rubella vaccine-induced immunity. PLoS ONE. 2010, 5 (7): e11806-10.1371/journal.pone.0011806.

    Article  PubMed Central  PubMed  Google Scholar 

  53. Schiller C, Nitschke MJ, Seidl A, Kremmer E, Weiss EH: Rat monoclonal antibodies specific for LST1 proteins. Hybridoma (Larchmt). 2009, 28 (4): 281-286. 10.1089/hyb.2009.0021.

    Article  CAS  Google Scholar 

  54. Mulcahy H, O'Rourke KP, Adams C, Molloy MG, O'Gara F: LST1 and NCR3 expression in autoimmune inflammation and in response to IFN-gamma, LPS and microbial infection. Immunogenetics. 2006, 57 (12): 893-903. 10.1007/s00251-005-0057-2.

    Article  CAS  PubMed  Google Scholar 

  55. Kilding R, Iles MM, Timms JM, Worthington J, Wilson AG: Additional genetic susceptibility for rheumatoid arthritis telomeric of the DRB1 locus. Arthritis Rheum. 2004, 50 (3): 763-769. 10.1002/art.20043.

    Article  CAS  PubMed  Google Scholar 

  56. Lavi S, McConnell JP, Rihal CS, Prasad A, Mathew V, Lerman LO, Lerman A: Local production of lipoprotein-associated phospholipase A2 and lysophosphatidylcholine in the coronary circulation: association with early coronary atherosclerosis and endothelial dysfunction in humans. Circulation. 2007, 115 (21): 2715-2721. 10.1161/CIRCULATIONAHA.106.671420.

    Article  CAS  PubMed  Google Scholar 

  57. Kim JY, Hyun YJ, Jang Y, Lee BK, Chae JS, Kim SE, Yeo HY, Jeong TS, Jeon DW, Lee JH: Lipoprotein-associated phospholipase A2 activity is associated with coronary artery disease and markers of oxidative stress: a case–control study. Am J Clin Nutr. 2008, 88 (3): 630-637.

    CAS  PubMed  Google Scholar 

  58. Kostopanagiotou G, Avgerinos E, Costopanagiotou C, Arkadopoulos N, Andreadou I, Diamantopoulou K, Lekka M, Smyrniotis V, Nakos G: Acute lung injury in a rat model of intestinal ischemia-reperfusion: the potential time depended role of phospholipases A(2). J Surg Res. 2008, 147 (1): 108-116. 10.1016/j.jss.2007.07.023.

    Article  CAS  PubMed  Google Scholar 

  59. Sutton BS, Crosslin DR, Shah SH, Nelson SC, Bassil A, Hale AB, Haynes C, Goldschmidt-Clermont PJ, Vance JM, Seo D, et al: Comprehensive genetic analysis of the platelet activating factor acetylhydrolase (PLA2G7) gene and cardiovascular disease in case–control and family datasets. Hum Mol Genet. 2008, 17 (9): 1318-1328. 10.1093/hmg/ddn020.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Wang Q, Hao Y, Mo X, Wang L, Lu X, Huang J, Cao J, Li H, Gu D: PLA2G7 gene polymorphisms and coronary heart disease risk: a meta-analysis. Thromb Res. 2010, 126 (6): 498-503. 10.1016/j.thromres.2010.09.009.

    Article  CAS  PubMed  Google Scholar 

  61. Sadegh-Nasseri S, Chen M, Narayan K, Bouvier M: The convergent roles of tapasin and HLA-DM in antigen presentation. Trends Immunol. 2008, 29 (3): 141-147. 10.1016/j.it.2008.01.001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Turnquist HR, Petersen JL, Vargas SE, McIlhaney MM, Bedows E, Mayer WE, Grandea AG, Van Kaer L, Solheim JC: The Ig-like domain of tapasin influences intermolecular interactions. J Immunol. 2004, 172 (5): 2976-2984.

    Article  CAS  PubMed  Google Scholar 

  63. Garbi N, Tiwari N, Momburg F, Hammerling GJ: A major role for tapasin as a stabilizer of the TAP peptide transporter and consequences for MHC class I expression. Eur J Immunol. 2003, 33 (1): 264-273. 10.1002/immu.200390029.

    Article  CAS  PubMed  Google Scholar 

  64. Williams AP, Peh CA, Purcell AW, McCluskey J, Elliott T: Optimization of the MHC class I peptide cargo is dependent on tapasin. Immunity. 2002, 16 (4): 509-520. 10.1016/S1074-7613(02)00304-7.

    Article  CAS  PubMed  Google Scholar 

  65. Raghuraman G, Lapinski PE, Raghavan M: Tapasin interacts with the membrane-spanning domains of both TAP subunits and enhances the structural stability of TAP1 x TAP2 Complexes. J Biol Chem. 2002, 277 (44): 41786-41794. 10.1074/jbc.M207128200.

    Article  CAS  PubMed  Google Scholar 

  66. Rufer E, Leonhardt RM, Knittler MR: Molecular architecture of the TAP-associated MHC class I peptide-loading complex. J Immunol. 2007, 179 (9): 5717-5727.

    Article  CAS  PubMed  Google Scholar 

  67. Fukuda A, Tokonabe S, Hamada M, Matsumoto M, Tsukui T, Nogi Y, Hisatake K: Alleviation of PC4-mediated transcriptional repression by the ERCC3 helicase activity of general transcription factor TFIIH. J Biol Chem. 2003, 278 (17): 14827-14831. 10.1074/jbc.M213172200.

    Article  CAS  PubMed  Google Scholar 

  68. Kim YK, Bourgeois CF, Pearson R, Tyagi M, West MJ, Wong J, Wu SY, Chiang CM, Karn J: Recruitment of TFIIH to the HIV LTR is a rate-limiting step in the emergence of HIV from latency. EMBO J. 2006, 25 (15): 3596-3604. 10.1038/sj.emboj.7601248.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  69. Briggs FB, Goldstein BA, McCauley JL, Zuvich RL, De Jager PL, Rioux JD, Ivinson AJ, Compston A, Hafler DA, Hauser SL, et al: Variation within DNA repair pathway genes and risk of multiple sclerosis. Am J Epidemiol. 2010, 172 (2): 217-224. 10.1093/aje/kwq086.

    Article  PubMed Central  PubMed  Google Scholar 

  70. Young LR, Borchers MT, Allen HL, Gibbons RS, McCormack FX: Lung-restricted macrophage activation in the pearl mouse model of Hermansky-Pudlak syndrome. J Immunol. 2006, 176 (7): 4361-4368.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  71. Sammon JW: A nonlinear mapping for data structure analysis. IEEE Transactions on Computers. 1969, 18: 401-409.

    Article  Google Scholar 

  72. Ljungberg K, Holmgren S, Carlborg O: Simultaneous search for multiple QTL using the global optimization algorithm DIRECT. Bioinformatics. 2004, 20 (12): 1887-1895. 10.1093/bioinformatics/bth175.

    Article  CAS  PubMed  Google Scholar 

  73. Wang X, Agarwala R, Capra J, Chen Z, Church D, Ciobanu D, Li Z, Lu L, Mozhui K, Mulligan M, et al: High-throughput sequencing of the DBA/2J mouse genome. BMC Bioinformatics. 2010, 11 (Suppl 4): O7-10.1186/1471-2105-1111-S1184-O1187.

    Article  PubMed Central  Google Scholar 

Download references

Acknowledgments

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 and the virtual institute ‘GeNeSys’ funded by the Helmholtz Association. The funders have no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Mice for these experiments were maintained by the animal caretakers of the Central Animal Facilities at the HZI. We wish to thank Christin Fricke for excellent technical assistance.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Klaus Schughart.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

TN, MH and KS conceived and designed the experiments. TN and MH infected the mice and collected the data. HK, RA, RW and KS performed and interpreted the QTL analysis. FK performed the analysis of environmental and experimental factors. SS and LL provided BXD mice for experimental studies. HK, RW and KS wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2012_4293_MOESM1_ESM.xls

Additional file 1: Table S1. Body weight loss from day 1 until day 7 for all strains showing the mean values, SEM per day and number of mice per strain analyzed. (XLS 50 KB)

12864_2012_4293_MOESM2_ESM.jpeg

Additional file 2: Figure S1. Differential susceptibility to influenza A infection among different BXD strains, parental strains and F1 generation. Rank-ordered strain distribution pattern illustrating the percentage of surviving mice per strain after influenza A infection for 53 BXD, parental DBA/2J and C57BL/6J strains and F1 (B6D2F1) mice for day 7 (A, trait ID: 13001) and day 8 (B, trait ID: 13002). Rank-ordered strain distribution pattern showing mean time to death for 53 BXD, parental DBA/2J and C57BL/6J strains and F1 (B6D2F1) mice (C, trait ID: 12996). (JPEG 355 KB)

12864_2012_4293_MOESM3_ESM.jpeg

Additional file 3: Figure S2. Principal component analysis of all body weight traits. Percent contribution of principal components to the total variance for body weight traits from day 1 until day 7 p.i. PC1 contributes about 80% and PC2 about 10% to the total variance. (JPEG 68 KB)

12864_2012_4293_MOESM4_ESM.jpeg

Additional file 4: Figure S3. Body weight changes in mock-infected C57BL/6J and DBA/2J mice. Female DBA/2J (n=4) and C57BL/6J (n=3) mice were intranasally inoculated with 25 μl PBS under anesthesia. Body weight changes for each group of treated mice at various days p.i. is shown with reference to the starting weight (% body weight). Data represent mean values +/- SEM. (JPEG 218 KB)

12864_2012_4293_MOESM5_ESM.jpeg

Additional file 5: Figure S4. Genome-wide linkage analysis for survival after H1N1 infection at days 9 and 10 p.i. Analysis of survival from day 9 (A) and day 10 (B) after PR8 infection revealed suggestive QTLs on chromosomes 5, 16, 17 and 19. Interval mapping is shown across the entire genome as blue line representing the LRS scores; green line: DBA/2J alleles increase trait values; red line: C57BL/6J alleles increase trait values. Upper x-axis shows mouse chromosomes, lower x-axis shows physical map in mega bases for each chromosome, y-axis represents LRS score. The horizontal lines mark the genome-wide significant thresholds at p0.05 (red line) and suggestive thresholds at p0.37 (gray line). (JPEG 616 KB)

12864_2012_4293_MOESM6_ESM.jpeg

Additional file 6: Figure S5. Pair Scan analysis of trait 13002. The table provides the breakdown of trait 13002 as an example for a pair scan analysis study that compares the fit of the alternative models. An analysis of this type can be regenerated rapidly on any of the traits presented in our manuscript. (JPEG 785 KB)

12864_2012_4293_MOESM7_ESM.jpeg

Additional file 7: Figure S6. Scheme to identify potential candidate genes in the analyzed QTL intervals. The detailed steps are described in the main manuscript. (JPEG 758 KB)

12864_2012_4293_MOESM8_ESM.xlsx

Additional file 8: Table S2. List of genes in the QTL intervals that were studied in detail. The QTLminer tool in GeneNetwork was to identify all genes in the QTL intervals and to obtain associated GO terms and KEGG pathways. All annotated genes within the QTL intervals that were expressed between day 1 and day 60 after infection of C57BL/6J mice with PR8M virus were further analyzed for the following attributes: genes that were up- or down-regulated in infected lungs after PR8M infection, genes carrying an insertion or deletion or a non-synonymous nucleotide change in the open reading frame, genes exhibiting a cis-expression QTL in the lung, and genes that were differentially expressed in C57BL/6J and DBA/2J mice in infected lungs between day 1 and day 8 p.i. Abbreviations: FS: frame shift in coding region, SNP: single nucleotide polymorphism, cis-eQTL: cis-regulated expression QTL, B6: C57BL/6J, D2: DBA/2J, ns: non-synonymous codons, stop_G: stop codon gained, stop_L: stop codon lost, DEL: deletion, INS: insertion; ns (…): number of non-synonymous codons. Legends: Expression d0-d60, log2 ≥8 at any day during infection; No: log2 8; ?: no expression data available. Regulation during infection (d0-d60): + - up or down regulation, genes that exhibited at least a 2-fold difference in changes of the expression levels in lungs of infected mice compared to non-infected controls; (+)(-): up-or down-regulation, genes that exhibited at least a 1.5-fold difference in changes of the expression levels in lungs of infected mice compared to the non-infected controls; no: no regulation. SNP: ns (.) number of non-synonymous codons; STOP_L: Stop codon lost, STOP_G: Stop codon gained; no: no SNPs between B6 and D2. Indel: DEL: Deletion; INS: Insertion; no: no INDEL. cisQTL: yes (…): LRS value has to be greater than ≥18; no: no QTL or LRS 18. diff expressed (B6-D2) during infection: Genes were defined as differentially expressed (log2 ≥8) at any day in infected lungs of C57BL/6J and DBA/2J; + -: Genes were defined as differentially expressed exhibiting at least a 2.0 difference in changes of the expression levels in infected lungs of C57BL/6J and DBA/2J; (+)(-) differentially expressed: Genes were defined as differentially expressed exhibiting at least a 1.5 difference in changes of the expression levels in infected lungs of C57BL/6J and DBA/2J; no: no regulation or not expressed; Score: the different attributes were counted. (XLSX 296 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Nedelko, T., Kollmus, H., Klawonn, F. et al. Distinct gene loci control the host response to influenza H1N1 virus infection in a time-dependent manner. BMC Genomics 13, 411 (2012). https://doi.org/10.1186/1471-2164-13-411

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-13-411

Keywords