|Home | About | Journals | Submit | Contact Us | Français|
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.
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 13days 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.
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.
Influenza A virus represents a major health threat to humans. The 1918 H1N1 pandemic caused at about 30 to 50 million deaths . Seasonal influenza epidemics cause high economic loss, morbidity and deaths every year . 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 .
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–7days post infection (p.i.), whereas C57BL/6J mice loss weight until days 6–8 after infection and regain their initial weight by 14days 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 . 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 ).
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 13days 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.
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 ) and followed body weight and survival over the next 13days. 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 13days post infection (Figure1, 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.
For body weight loss and survival, three different phenotypic response groups can be defined (Figure1). In the first (Figure1A), all infected mice within a strain survived, in the second (Figure1B) a majority but not all individuals within a strain survived, and in the third group (Figure1C), 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 (Figure1A). In contrast, BXD28 belonging to the third group, lost body weight much more rapidly than even the highly susceptible DBA/2J parent (Figure1C). 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. (Figure2). 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.(Figure2), 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.
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 (Figure1B, C). After day 7, surviving mice started to gain weight (Figure1A, 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 . This QTL exhibited genome-wide significance on days 5 (LRS: 19.0, effect size: 30%) and 6 (LRS: 19.4, effect size: 31%, Figure3). However, the effect of Qivr5, although weaker, can also be detected on days 2, 3, and 4 and 7 (Figure3). 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 (Figure3A-G).
Another suggestive QTL peak was found on chromosome 10 on day 1 (LRS: 14.1, effect size: 23%, Figure3A). Its effect is lower at day 2 and not apparent at later days. These observations are in accordance with the PCA (Figure2) 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 . 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%, Figure4A) as well as suggestive QTL peaks on proximal chromosome 16 and on chromosome 17 (Figure4A). 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 (Figure4B) which was not seen in any other trait.
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 (Figure5A-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% Figure5A, B) and its effect was still evident from day 9 until 11 p.i. (Additional file 5: Figure S4 and Figure5C). 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 (Figure5A-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 (Figure5D). In addition, suggestive QTLs were found on chromosomes 2 and 17 for the MTTD trait.
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.
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 (Figure6A, 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 (Figure6C, D).
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 (Figure7A). The inclusion of an interaction term increases the LRS by 23 (Additional file 6: Figure S5). For both Qivr5 - chromosome 9 interaction (Figure7B) and Qivr5 – Qivr19 interaction (Figure7C), the allele combination DBA/2J / DBA/2J showed highest and the combination C57BL/6J / C57BL/6J lowest survival scores. All differences were highly significant.
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 (Figure8). 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  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 , 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.  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 Table1. 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 Table2). These genes were further evaluated based on their known function from the literature and phenotypes in knock-out mouse mutants (see 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 . 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  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  of a total of 26 susceptible strains. Furthermore, five strains that were resistant in our study were susceptible in the study by . 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 . 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 . 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 (Table2). 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/2J 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 (Table2) by using additional information, such as temporal expression after PRM8 infection (Pommerenke et al., PLoS ONE, in press), cis-eQTLs in non-infected lungs , differences in expression between DBA/2J and C57BL/6J , 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 . 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 . 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 . 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 . 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 . 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 . 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) . It is involved in modulating pro-inflammatory responses in macrophages .
Qivr17-2 represents a positive influence of the C57BL/6J 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 . 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/2J 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 . Pla2g7 expression levels in susceptible A/J mice were higher than in resistant C57BL/6J mice after infection with H3N2 virus . 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 . Tnfrsf21 which was identified by  as potential candidate of Qivr17-2 also exhibits a cis-eQTL in non-infected BXD mice  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 .
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 3days 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).
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.
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).
The mouse-adapted influenza strain A/Puerto Rico/8/34 (H1N1; PR8M,  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 (50mg/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 13days 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%.
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 . 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 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 . 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.
The QTL region analysis was initially performed using the QTLminer which has been implemented in GeneNetwork . 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 . 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.
The authors declare that they have no competing interests.
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.
Tatiana Nedelko and Heike Kollmus contributed equally as first authors.
Robert W. Williams and Klaus Schughart contributed equally as senior authors.
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.
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).
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.
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.
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).
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.
Figure S6. Scheme to identify potential candidate genes in the analyzed QTL intervals. The detailed steps are described in the main manuscript.
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.
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.