The MHC class I Dk molecule supplies vital host resistance during murine cytomegalovirus (MCMV) infection. Natural killer (NK) cells expressing the Ly49G2 inhibitory receptor, which specifically binds Dk, are required to control viral spread. The extent of Dk-dependent host resistance, however, differs significantly amongst related strains of mice, C57L and MA/My. As a result, we predicted that relatively small-effect modifier genetic loci might together shape immune cell features, NK cell reactivity, and the host immune response to MCMV. A robust Dk-dependent genetic effect, however, has so far hindered attempts to identify additional host resistance factors. Thus, we applied genomic mapping strategies and multicolor flow cytometric analysis of immune cells in naive and virus-infected hosts to identify genetic modifiers of the host immune response to MCMV. We discovered and validated many quantitative trait loci (QTL); these were mapped to at least 19 positions on 16 chromosomes. Intriguingly, one newly discovered non-MHC locus (Cmv5) controlled splenic NK cell accrual, secondary lymphoid organ structure, and lymphoid follicle development during MCMV infection. We infer that Cmv5 aids host resistance to MCMV infection by expanding NK cells needed to preserve and protect essential tissue structural elements, to enhance lymphoid remodeling and to increase viral clearance in spleen.
Yellow fever virus (YFV), once a major scourge of humanity, was one of the first viruses studied experimentally in mammalian hosts [1]. In pioneering studies, Nobel Laureate Max Theiler developed an inactivated YFV vaccine [2], and laid the groundwork for investigations into the genetic basis of host resistance to virus infection. Sawyer and Lloyd later observed that different strains of white mice are differently susceptible to YFV [3], and Lynch and Hughes solidified the point that YFV susceptibility is a heritable trait [4]. Many years later, the underlying cause of disease and effect of genetic variance on host resistance to viral infection and pathogenesis is still of vital interest [5], as it promises to reveal yet unknown molecular targets, signaling pathways and cellular networks with relevance to human health and disease. Genetic analysis of host resistance to MCMV has been especially rewarding [6–8]. Recently discovered genes encode virus sensors and ligands, cytokines and receptors, signal transducers, and effector molecules that either increase or decrease host resistance to infection [9–16]. Often, these molecules are related to cellular immunity, including a clutch of polymorphic NK cell receptors that specifically bind ligands on virus-infected cells [17–24]. Yet, our understanding of the genetic influences on NK cells in the response to viral infection remains incomplete. We established several mouse models to explore the effect of MHC class I (MHC I) polymorphism on NK cells in viral immunity [25]. MHC I Dk confers dominant MCMV resistance in MA/My and C57L-derived transgenic Dk mice, while Db-expression in C57L and MA/My-derived congenic M. H2b mice does not [26]. The Dk resistance effect requires NK cells that express Ly49G2 (G2), an inhibitory receptor that binds Dk and protects against viral spread [24]. Thus, decreased Dk expression on infected cells might release the G2-specific brake on NK stimulatory signals, therefore helping to eliminate MCMV targets [25,27]. However, G2' s precise role in host resistance is still under investigation. While analyzing the MHC effect on NK cells, we found that Dk-dependent MCMV resistance is greater in C57L-derived mice, relative to MA/My mice [26]. Thus, C57L genetic modifiers may increase host resistance to infection. However, modifier genetic loci have so far eluded detection, likely due to the prominent role of Dk. Although forward and reverse genetics approaches have uncovered many pathogen resistance genes, neither strategy is ideally suited to resolve smaller genetic effects. Moreover, reverse genetics relies on introducing novel mutations that result in phenotypic abnormalities, so it is not a practical way to identify or characterize natural allele variants with distinct effects on immune function or pathogen resistance. We thus set out to map and characterize variable genetic effects (quantitative trait loci, QTL) that shape immune features in naïve and MCMV-infected mice. We expected that virus-responsive immune cell traits and those that provide critical MCMV resistance would be revealed through common QTL linkages at defined chromosome positions across the genome. To guard against spurious results, we generated an additional cohort of MHC I Dk-disparate offspring that were given lower dose MCMV infection. We found that defined genomic regions controlled either naïve or MCMV-responsive NK cell features, while others affected both sets of traits. While the MHC I D molecule has a central role in virus resistance and NK cell functionality, genetic dissection allowed many additional immune cell and host response traits to be separately mapped and verified. One such newly discovered non-MHC QTL remarkably controlled the integrity of the splenic environment, the prevalence of NK cells, and overall resistance to virus infection.
Previously we assessed the immune response to MCMV with a deep phenomic approach [23]. The original high dose (HD) cohort [23] and a second cohort of MA/My x C57L hybrid mice given lower dose (LD) infection are included in the current analysis (S1 Fig). All mouse genome-wide genotypes were assessed and verified using automated SNP analysis on an Illumina platform. We performed genome scans to detect and map host resistance and immune response modifier QTL. Briefly, pre- and postinfection quantitative traits from both cohorts were separately analyzed using the quantitative genetic mapping package R/qtl in the statistical computing program R [28]. Table 1 consolidates the data, conservatively reporting the most significant genetic intervals (QTL) detected. Results are ordered by chromosome, and then by position with overlapping genetic intervals for the same trait detected in both cohorts, or distinct traits detected in at least one cohort listed. As expected, sex-related naïve mouse body weight (Mbw-X) QTL mapped to overlapping genetic intervals on chr-X in both cohorts, but with slightly different peak logarithm of the odds (LOD) linkage positions, which relates the probability that an observed linkage was not due to chance (Fig 1A and 1B, Table 1). MCMV-induced body weight postinfection (Vbwp-X) QTL also mapped to this interval on chr-X, with near identical positions detected in LD and HD cohorts. As the only Vbwp QTL detected after LD infection, the results further suggested that MCMV-induced weight change was especially sensitive to chr-X weight control. Several smaller-effect mouse body weight QTL (Mbw-3, Mbw-7, Mbw-19) were also found and these coincided with postinfection weight control QTL (Vbwp-3, Vbwp-7, Vbwc-7, Vbwp-19) detected on chromosomes 3,7 and 19, respectively (Table 1). Thus, several QTL were discovered that affected body weight in naive and MCMV-infected mice. Host immune cell and response traits were likewise analyzed. Of the fifteen newly discovered NK cell trait modifier (Nktm) QTL, eight mapped respectively in both cohorts (Table 1). Twenty-four cohort-specific QTL were also identified (Table 1), which suggested that subtle differences in infectious dose, cohort-size, or genetic make-up (see S1 Fig) could have revealed additional genetic associations. In general, broadly profiling genetically diverse mice with genome scans for immune traits provided a sensitive and reliable measure of the host response to viral infection, and led to the discovery of many novel QTL in control of immune cell homeostasis and reactivity. As expected based on prior studies [20,25,29,30], a strong MCMV (log spleen) resistance QTL mapped to H-2D on chr-17 in both cohorts (Fig 1A and 1B, Table 1). No other single QTL was mapped with greater precision. Two critical MCMV-induced weight-control QTL (Vbwp-17, Vbwc-17) also mapped at the H-2D locus, but only in HD-infected mice (Fig 1, Table 1). Thus, these data demonstrated that H-2D controlled viral spread and the severity of weight loss (i. e. morbidity) solely after HD infection. The prominence of Dk in host resistance led us to next assess its effect on NK cells in the response to MCMV. In fact, several NK cell modifier QTL in control of NKp46+ or G2+ NK cells and their response to MCMV were detected on chr-17. Interestingly, both NK cell trait (Nktmg2sp-17, Nktmg2spf-17) and response (Nkrmg2sp-17, Nkrmg2spf-17) modifiers of G2+ NK cells mapped to the H-2D locus. In contrast, QTL regulating Ly49I/U+ (I/U+) NK cells had no association to chr-17, and instead mapped to the NKC on chr-6 (Fig 1A and 1B, Table 1). Importantly, H-2D polymorphism specifically controls baseline G2+ NK cell features in naïve mice in addition to the G2+ NK response to MCMV, virus clearance and morbidity indices. To pursue genetic modifiers of MCMV immunity, we performed two-dimensional genome scans in R/qtl and analyzed each trait in possible two-locus QTL models. A higher LODav1 value is indicative of an additive effect, without evidence of epistasis, while a higher LODint value is suggestive of genetic interaction. Two-locus modeling in R/qtl detected numerous QTL on a range of chromosomes with decidedly significant LODav1 values, including many already discovered in Table 1. Thus, these QTL added to H-2D’s (pos 17. 8 in Table 2) effect on MCMV burden or body weight, without evidence of interaction (Table 2, S2A Fig). Interestingly, seven common QTL positions precisely mapped on chromosomes 1,3, 6 (pos 21. 7, separate from the NKC), 11,17 (pos 12. 8, separate from the MHC), 18 and X added to both MCMV and body weight control, which suggested that a single shared QTL on each chromosome enhanced H-2D' s effect on both traits. Three additional QTL that mapped at similar positions on chromosomes 7,8 and 12 also added to H-2D' s effect on MCMV and body weight control, which suggested these too mediate a shared locus effect. Five additional QTL on chromosomes 2,9, 10,13 and 19 had wider-ranging positions that added to H-2D' s control of either trait. Thus, we could not resolve whether common or distinct QTL on these chromosomes had a role. Nonetheless, both chr-13 (pos 37. 1) and chr-19 (pos 49. 1) QTL intervals that affected MCMV burden also overlapped with postinfection body weight intervals (Table 2). Several Nktm and Nkrm QTL mapped at the same chr-13 interval (Tables 1 and 2), further established strong support of at least one QTL on chr-13. While a chr-2 MCMV control QTL (pos 16. 8) coincided with Nktmg2n-2 and Nktmg2sp-2, a weight-control QTL (pos 31. 8) corresponded better with two MCMV-induced weight control (Vbwp-2, Vbwc-2) QTL locations and the NK response modifier, Nkrmdn-2 (Tables 1 and 2). Thus, two QTL were reliably detected on chr-2. Lastly, a chr-10 QTL for MCMV burden coincided with NK trait modifier QTL (Nktmg2n-10, Nktmg2sp-10), but a chr-10 (pos 13. 4) weight control QTL could not be independently verified. In aggregate, the results implicated at least sixteen new QTL that regulate NK cells, virus clearance and morbidity indices, and most corroborated QTL reported in Table 1. Evidence of epistasis or interactive QTL (iQTL) pairs with significant LODint values that affected MCMV burden or weight control was also uncovered (Table 2, S2B Fig). Of these, all but the QTL on chr-5 and -14 could be independently verified in the genome scans of other traits. Interestingly, a chr-4 weight control iQTL that mapped similarly to Cmv6, coincided with an additive H-2D effect on MCMV burden (Table 1), in addition to its interaction with a chr-12 iQTL that further added to H-2D-dependent control of MCMV clearance and morbidity (Table 2). Thus, two-dimensional genome profiling uncovered and validated at least 18 novel QTL that affected viral spread or weight control during MCMV infection (S2 Fig). We infer that these genetic modifiers of NK cells and body weight together strongly influence overall host resistance to viral infection in both H-2D-dependent and independent ways. NKC-encoded receptor polymorphism has been shown to affect NK cells in the response to MCMV, and MHC-NK gene complex (NKC) epistasis contributes to MCMV control [20,30]. As expected, several modifier-QTL in control of NK cells, and I/U+ or G2+ NK subsets, were mapped to the NKC on chr-6 (Tables 1 and 2; Fig 2A and 2B). Chr-6, however, had no obvious effect on viral clearance or morbidity, so it was unlikely that Ly49 (Klra) polymorphism alone affected the variance in MCMV resistance in this genetic comparison. Nonetheless, genome profiling had uncovered a chr-6 QTL that added to H-2D’s effect on MCMV burden and weight control (Table 2). This result corresponds to a prior study showing that the C57L-derived NKC (NKCl) enhanced NK cell-mediated resistance after HD-infection in C57L mice without Dk [31]. Thus, we tested NK cell responsiveness by stratifying the phenotypic results by MHC and NKC genotypes. Interestingly, LD mice without Dk (i. e. H-2b homozygous mice designated LL) had a lower viral burden if they also had a C57L-derived NKC (Fig 2C). A similar relationship was observed in HD mice, though not significant. While a single amino acid variation distinguishes G2 receptors in MA/My and C57L [24,32], the results suggested that G2 receptor polymorphism could have affected NK cell sensitivity to MCMV. Further analysis of H-2b homozygous offspring revealed that naïve G2-SP NK cells with C57L-derived G2c57l receptors were fewer with lower G2 receptor MFI values, in comparison to mice solely expressing MA/My-derived G2mamy receptors (Fig 2C). Whereas MHC I ligands of inhibitory NK receptors are known to affect NK cell and receptor expression features [33,34], the data suggested that the variant G2 receptors might differently bind H-2b class I or class I-related molecules, which could have further resulted in MCMV resistance differences. To corroborate our analysis of NK cells in H-2b homozygous offspring, we next examined NKC congenic mice. As with the hybrid NK cells, naïve G2+ NK cell percentages and G2 MFI values were lower in NKCc57l than NKCmamy congenic mice (Fig 2D), which suggested that the disparate G2 receptors might also differently license NK cells. However, comparisons of NK cells in NKC congenic mice (S3 Fig), or C57L and M. H2b mice (S4A Fig) revealed similar sensitivity to ex vivo triggering via several different NK stimulatory receptors, and consequently no difference in NK cell licensing. Lower viral loads in NKCc57l than in NKCmamy spleens (Fig 2C and 2D, [31]), nonetheless prompted further analysis of G2+ NK cell reactivity to MCMV. In agreement with results from C57L-derived NKC congenic mice, lower G2+ NK cell percentages and G2 receptor MFI trends were observed in C57L mice, in comparison to MHC-matched, but NKC-disparate M. H2b mice (S4B Fig). More intriguingly, host resistance to MCMV in C57L exceeded M. H2b, which directly corresponded to higher postinfection G2+ NK cell percentages and G2 receptor MFI (S4B–S4D Fig). C57L mice also exhibited significantly increased frequencies and numbers of total NK and G2+ NK cells producing greater amounts of IFN-γ (S4E–S4G Fig). Dk-independent host resistance to MCMV therefore was also magnified in the C57L genetic background. We infer that G2+ NK cell phenotypes and effector functions are distinctly shaped by and sensitive to both NKC and non-NKC genetic effects, which further controls host-specific variations in virus response features. Given the robust Dk effect on NK cells and virus resistance, it was not surprising that the postinfection percentage of NKp46+ NK cells in spleen also mapped to chr-17 (Table 1, Fig 3A and 3B). Precision mapping in the HD cohort, however, resolved a peak QTL location (rs13483002) distal to H-2D (Fig 3B). The data implicated that a distinct locus outside of the MHC on chr-17 affected NK cell percentages after infection. To address the question, we compared allele effect plots for several pre- and postinfection traits. We found that Dk had a profound effect on G2 receptor expression and NK cell subsets in naïve animals (Fig 3C). It also had a dominant effect on spleen MCMV burden after LD- or HD-infection (Fig 3D, compare H-2D genotype classes, LL and LM). More remarkable, Dk had a major impact on postinfection percentages of G2-SP NK cells and the severity of weight loss after infection (Fig 3D). In contrast, a C57L-derived locus (rs13483002) either dominantly prevented or failed to support NKp46+ NK cell accrual after HD-infection (Fig 3D, right panel). This effect was specific to NK cells since Dk affected the overall balance of total lymphocytes in the spleen (Table 1, Fig 3E). A CI for the QTL controlling percentages of NKp46+ NK cells was mapped distal to Dk in the HD-infected cohort (Fig 3E). In support of this location, a very narrow R7-MHC congenic block (~300-kb Dk-spanning segment [25]) was included in the genetic crosses used to generate the HD cohort so that many offspring had recombinant chromosomes between MHC I D and rs13483002. The results demonstrated that the Dk-resistance factor alone was inadequate to support NK accumulation in spleen after HD-infection. Thus, at least two different chr-17 QTL control NK cells responding to MCMV. To corroborate the rs13483002 QTL effect, we compared the NK-cell response to infection in MA/My and two recently generated MA/My-derived strains, M. H2b and M. H2b-TgDk (hereafter referred to as Tg1), which expresses a genomic Dk transgene [25]. MA/My and Tg1 mice display MCMV resistance, whereas M. H2b mice are highly susceptible to infection. Moreover, the M. H2b congenic interval spans the rs13483002 locus (for maps of R7 and M. H2b congenic blocks, see Ref [23]). Overall NK cell percentages and numbers, in addition to spleen sizes, were equivalent in uninfected mice (Fig 4A and 4B). Within days after infection, however, MA/My spleens had grown much larger with significantly increased frequencies of NK1. 1+ (S5 Fig) and NKp46+ NK cells (Fig 4A and 4B), which corresponded to lower viral loads (Fig 4C). In contrast, M. H2b spleens failed to increase in size or frequency of NK cells after infection, regardless of viral dose, which resulted in greatly increased viral burden. Similar differences distinguished MA/My and Tg1 after HD-infection, regardless of enhanced viral control in Tg1 (compared to M. H2b mice). Thus, Dk expression by itself is not sufficient to support accumulation of spleen NK cells in response to MCMV. Higher virus levels in Tg1 than MA/My mice further suggested that a C57L-derived susceptibility allele in the M. H2b congenic interval is sufficient to decrease overall resistance. To test this, we next examined (MA/My x M. H2b) F1 cross mice after MCMV infection. Remarkably, M. H2k/b spleens were smaller with significantly higher MCMV levels, in comparison to MA/My (Fig 4D). Thus, these data demonstrated that a C57L-derived susceptibility allele, provisionally designated Cmv5s, interfered with the NK-mediated resistance effect of Dk. Moreover, the results established that diminished control in M. H2k/b is unlikely related to a Dk gene dose effect, since the same allele also interfered with host resistance in Tg1 mice with multiple integrated Dk gene copies and slightly higher protein expression [32]. Significantly higher percentages of NKp46+ NK cells and larger spleens in infected Cmv5r resistant MA/My mice provide added support of this interpretation. We infer that enhanced host resistance in MA/My due to increased accrual of NK cells is prevented by expression of at least one Cmv5s allele in more susceptible MA/My-derived (e. g. M. H2b) mice. Bekiaris et al. have previously shown that Ly49H+ NK cells are needed to protect spleen white pulp after MCMV infection [35]. Thus, we investigated whether NK cells similarly protected splenic architecture in Cmv5-disparate MA/My mice by comparing secondary lymphoid organ (SLO) structures after MCMV infection. Remarkably, within four days after infection, we observed decidedly ill-defined white pulp (WP) regions that were smaller and had less clearly delineated marginal zones (MZ) in M. H2b and Tg1 mice, in comparison to infected MA/My or uninfected control mice (Fig 5A). Infected MA/My mice had more abundant lymphoid-like cells, in addition to pseudo-nodules of loosely packed mononuclear cells that were heterogeneous in size and shape, appeared to be activated with many mitotic figures, and located at the MZ at the edge of the WP (Fig 5B). In parallel, we observed severe necrosis in the expanded red pulp (RP), which was infiltrated with cells with viral inclusions, granulocytes, and large cells that resemble plasma cells in M. H2b and Tg1 spleens, as compared to infected MA/My or uninfected controls (Fig 5B). To validate the histological findings and Cmv5 genetic regulation of SLO structure after MCMV infection, we next examined spleen sections for NK cells and SLO structural features by IF microscopy. As expected, we observed equivalent numbers of intact NKp46+ cells (ring-like staining patterns) in splenic RP of each strain before infection (Fig 5C). Afterward, NKp46+ NK cells were more closely associated with MAdCAM+ cells at the MZ (Fig 5C). In infected MA/My spleens, we observed intact NKp46+ cells, in addition to fragmented NKp46+ granular material, possibly representing NK cell debris. In contrast, intact NKp46+ cells were rare in infected M. H2b and Tg1 spleens. Instead, extensive NKp46 staining of fragmented granular material, frequently concentrated in foci several times larger than an intact NK cell, was typical of infected M. H2b and Tg1 spleens. These data suggested that necrotic NK cells and debris were much more prevalent in mice without Cmv5 protection. We also examined MAdCAM+ marginal sinus-lining cells and CD169+ macrophages at the MZ. At four days postinfection, MA/My spleens consistently had more discrete (rather than more diffuse) MZ staining patterns, than either M. H2b or Tg1 (Fig 5D). Nonetheless, IF staining for T and B cells verified even more consistent differences and a selective loss of SLO structural features in M. H2b and Tg1 mice after MCMV infection (Fig 5E). Prominently, in M. H2b and Tg1 spleens, WP was marked by a lack of association between B220+ cells and the MZ. In addition, we observed disorganized regions of T and B cells, and a loss of B cell zone integrity. Together, these results establish that the M. H2b congenic segment failed to protect spleen structures and NK cell accrual after infection. Related occurrences in Tg1 mice verified that Dk was inadequate to protect SLO structure after infection, even though it aided viral clearance compared to M. H2b (Fig 4). As the extent of viral control differs in MA/My and Tg1 spleens (Fig 4), we considered that apparent differences in splenic SLO integrity could be due to the variance in viral titers. In line with previous studies [23–25], we reasoned that G2+ NK-mediated resistance in MA/My and Tg1 is equivalent. Thus, we further assessed the effect of G2+ NK cells on SLO structure under both LD- and HD-infection conditions. As expected, G2+ NK cell depletion resulted in significantly diminished spleen weight, splenocyte numbers, and NK cell recovery after infection in both strains (S6A and S6B Fig). Virus levels were correspondingly higher in G2-depleted mice (S6C Fig). Thus, G2+ NK cells are critical. Nonetheless, residual MCMV resistance in G2-depleted MA/My still exceeded that observed in G2-depleted Tg1, which revealed that Cmv5 exerted G2+ NK-independent effects on spleen and MCMV clearance (S6C Fig). Moreover, although spleen structure again waned with increased viral burden, G2-depleted MA/My spleens retained better defined WP and a higher density of small lymphocytes, in comparison to virus level-matched Tg1 (S6C and S6D Fig; compare LD G2-depleted Tg1 to HD G2-depleted MA/My, and HD Tg1 spleen sections). Despite that infected cell inclusion bodies were observed in both strains, Tg1 spleen was remarkable for increased granulocytosis and fibrinoid necrosis (S6C and S6D Fig). These results therefore demonstrate that Cmv5' s effect on SLO structural integrity during infection is independent of MCMV titers, and G2+ NK-mediated resistance. We infer that the combined effect of Dk and Cmv5r homozygosity in MA/My intensifies host resistance by protecting or regenerating spleen SLO structural features, preventing severe necrosis in spleen RP and promoting NK cell expansion needed to mediate specific virus clearance.
Although small-effect QTL were thought to increase the extent of Dk-dependent resistance to MCMV in C57L, compared to MA/My background mice [26], until now they have eluded detection likely due to the robust effect of Dk. A genome-driven integrated approach with multiparametric flow cytometry revealed 56 novel immune and MCMV responsive QTL. Though improbable that any two separately measured traits (e. g. splenic MCMV level and percentage G2+ NK cells postinfection) would map similarly, unless controlled by the same, or another closely linked gene, the combined approach helped to distinguish significant genomic linkages. Importantly, the combined genome scanning approach resulted in the discovery and mapping of multiple QTL, which together shape strain-specific variances associated with MCMV infection. Consistent with its role in licensing G2+ NK cells [25,36], we observed that the H-2D locus regulated both naïve G2+ NK cell percentages and G2 receptor expression without specific effects on features in other NK cell subsets. H-2D also controlled G2+ NK cell percentages, and the severity of weight loss after MCMV exposure. Together the findings confirm that much of the genetic impact on host resistance is through Dk-specific regulation of G2+ NK cells. Intriguingly, an interactive NKC-linked QTL further shaped MHC I D genetic regulation of G2+ NK cells since G2c57l and G2mamy receptors were expressed differently on NK cells in non-Dk mice. Lower G2c57l receptor expression on a smaller percentage of NK cells is suggestive that these NK cells could be licensed, though not confirmed by ex vivo stimulation experiments. As the majority of NK cells express the given G2 receptor allele in C57L-derived NKCm or NKCl congenic mice, further analysis is needed to exclude potential effects due to other common licensing receptors. Moreover, whereas increased host resistance due to NKCl further suggests that G2 receptor polymorphism might affect MCMV clearance differences, further analysis is needed to precisely define the role of G2+ NK cells. In aggregate, these findings suggest that chr-6 holds two or more QTL: one or more in the NKC that influence NK cell traits, and another more proximal locus, reported in Table 2, that regulates postinfection body weight. Therefore, the best interpretation is that an additive chr-6 effect with H-2D on MCMV resistance requires a more proximal locus, possibly in addition to a NKC-linked QTL. Given H-2D’s profound influence over NK cells, it was not surprising that overall postinfection percentages of NKp46+ NK cells were controlled by a chr-17 QTL. The newly discovered Cmv5 locus, however, segregated away from H-2D in both the LD and HD analyses. Two-dimensional genome scans of MCMV resistance and morbidity indices validated that at least one non-MHC chr-17 iQTL added to the genetic effect of H-2D. Though mapped proximally in the LD cohort such that a Cmv5 CI spans the H-2D locus, we favor a more distal map location as obtained in the HD analysis, for several reasons. First, in addition to its larger size, the HD cohort included two different MHC congenic blocks (R7 and M. H2b), which increased both resolution (smaller CIs) and precision (more informative crossovers) in the genetic linkage analysis. Second, less morbidity amongst LD mice also might have limited the map analysis as Cmv5’s effect was most clear after HD-infection. Lastly, analysis of MA/My-derived congenic strains verified that Cmv5 is distinct from and resides distal to H-2D on chr-17. In addition to its effect on NK cells in infected spleen, we found that Cmv5 also regulates spleen size, SLO structure and organization, and red pulp necrosis, which is at least partly via its direct effect in spleen and independent of G2+ NK cells. One possibility is that the Cmv5r resistance allele protects against the loss of splenic organization and other changes by increasing NK-mediated control of viral replication. An intrinsic modifier of cell proliferation or effector function, for example, might facilitate a more rapid NK response and sensing of viral targets. In fact, m157-specific Ly49H+ NK cells were recently shown to protect lymphoid organization in spleen red and white pulp within days after MCMV infection [35]. Extrinsic cues (i. e. cytokines or growth factors), on the other hand, could also aid in accrual of NK cells via increased recruitment, retention or enhanced proliferation in infected spleen tissue. Cmv5 protection of, or enhancement of IL-15 production by, CD8α DC [37–39] or VCAM-1+ CD31- stromal cells in spleen [40] might explain increased NK accumulation. FMS-like tyrosine kinase 3 ligand (Flt3-L), another interesting candidate, is needed to generate pDC and resident DC from Lin-Flt3+ bone marrow precursors [41], which subsequently stimulates DC to increase NK cell activation after MCMV infection [42]. Alternatively, Cmv5’s action may be directed to support specific spleen cell function (s) independent of NK-mediated viral control, especially given that resistance in Tg1 surpasses that in M. H2b mice after LD- or HD-infection, without affecting other Cmv5 features. Thus, Dk-dependent MCMV resistance per se is inadequate to guard against virus-induced loss of splenic architecture. Whether a cell intrinsic effect or not, greater numbers of competent NK cells in the spleen that can limit viral spread would be expected to minimize excessive inflammation and tissue injury [43]. Recently, we performed whole exome sequencing of MA/My, C57L and M. H2b to elucidate potential genetic candidates for Cmv5’s effect (H. Lee, A. G. and M. G. B. , manuscript in preparation). In the critical genetic interval (~17-Mb) flanked by H2D and the M. H2b congenic recombination crossover marked by Sgol1,117 expressed sequences were assessed, including many with non-synonymous or frame shift variants that distinguish MA/My and C57L alleles. Notable among immune-related positional candidates are Tnfrsf21, and Pla2g7 implicated in influenza control [44], though its effect was limited to male mice. The Trem /Trem-like gene cluster [45] also maps nearby Cmv5' s peak position. As these molecules have been shown to regulate innate immune cells [45], inflammation [46,47], dead cell clearance [48,49] and disease [50–53], they too represent interesting positional candidates. During recovery from viral infection, the spleen and other organs can serve as sites of blood formation or extramedullary hematopoiesis (EMH) [54]. Indeed, spleen NK cells were also shown to promote effective EMH by specifically controlling viral spread [55]. It will be important to determine whether Cmv5 regulates EMH, either apart from or in unison with Dk-dependent NK-mediated MCMV resistance. Thus, this genetic model may extend beyond investigating the effect of MHC polymorphism on NK-mediated virus resistance to examine EMH regulation generally, and during viral infection. Genetic regulation of splenic NK cell accumulation might lend additional precision to the host immune response. Whether due to increased splenic recruitment, retention, or enhanced proliferation that result from intrinsic or extrinsic effects, the fewer numbers of NK cells in the spleen of mice that carry a C57L-derived Cmv5s susceptibility allele should likewise decrease NK-mediated virus resistance. Alterations in NK numbers, and possibly activation status, have the potential to influence other immune cells in the response to infection. Specific NK-mediated control of MCMV has been shown to protect CD8α DC, which reciprocally affect the expansion of NK cells after infection [37], and the balance of NK-DC crosstalk has considerable potential to affect T cell immunity [56,57]. During persistent LCMV infection, NK cells were also shown to regulate adaptive immunity through direct lysis of CD8+ and CD4+ T cell effectors [58,59], which is subject to type I IFN regulation of MHC I expression on T cells themselves [60,61]. Likewise in MCMV infection, too few NK cells or inefficient NK activity may result in over exaggerated effector T-cell function, while too many NK cells can also harm otherwise protective T cell immunity [62,63]. Intriguingly, this effect of NK cells extends to humoral immunity, as follicular helper T cells are also a target of NK-mediated lysis after LCMV and other viral infections [64,65]. Thus, Cmv5 regulation of the innate immune response to viral infection is an additional pivotal mechanism governing NK cell numbers and consequently, their impact on other critical immune cells in the host response, inflammation and tissue injury.
All animal experiments conducted in this study were carried out in accordance with the Animal Welfare Act and the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All experiments were approved by the University of Virginia Animal Care and Use Committee (Protocol Number: #3050). MA/My x C57L crossed offspring given HD MCMV were described previously [23]. Congenic blocks were included in the crosses to aid identification of small-effect, and possibly MHC- or NKC-linked QTL. Crossed mice without congenic blocks were also generated and analyzed following LD MCMV infection. A diagram of the breeding strategy is shown in S1 Fig. MA/My-derived MHC congenic M. H2b (i. e. MA/My. L-H2b) and M. H2b-TgDk (Tg1) were described before [25,26]. Tg1 mice have been backcrossed to M. H2b mice for 10+ generations, with marker-assisted genetic selection in the first 4–5 generations. Mice used in the study were managed with a Colony Management System (Jackson Labs, Version 4. 1. 2), maintained in a dedicated animal care facility under SPF conditions and treated in accordance with IACUC regulations and guidelines. Abs and isotype controls for flow cytometry were purchased from BD Biosciences, Biolegend, and R&D Systems. They included anti-CD16/32 (2. 4G2), Ly49G2 (4D11 and AT8), Ly49D (4E5), Ly49A (YE1/48. 10. 6), Ly49I/U (14B11), CD3 (145-2C11), CD19 (6D5) and NKp46 (29A1. 4). Abs were conjugated to FITC, PE, PerCP-Cy5. 5, APC, APC-Cy7 or biotin followed by streptavidin conjugated to APC-Cy7 (Biolegend). Leukocyte stains, dead cell exclusion, compensation and cytometric analyses were performed as described [23,36]. To determine the licensing status and functional cytokine production by NK cells, ex vivo stimulation and intracellular cytokine staining were performed as previously described [66]. For IF staining of leukocytes and MZ, spleens were snap frozen in liquid nitrogen, processed at the UVA Research Histology Core and examined as described [67]. Briefly, frozen sections (6 μm) were fixed with a 1: 1 mix of acetone/ethanol (100%) and stained with Abs against B220 (RA3-632, BD Bioscience), CD3 (145-2C11, Ebiosciences), MAdCAM-1 (MECA-367, BioXcell), CD209b/SIGN-R1 (eBio22D1, Ebiosciences), NKp46 (polyclonal, R&D systems), and CD169 (MOMA-1, AbDSerotec). Non-specific staining and endogenous peroxidase activity were blocked using the Avidin/Biotin blocking kit (Vector Labs), and 3%H2O2 + 10mM NaN3 solution, respectively. The reaction was amplified using the tyramide signal amplification technique (Perkin Elmer, Boston MA). Slides were incubated with primary Abs (1h), followed by secondary biotinylated or FITC-conjugated Abs (Vector Labs). Slides were then incubated with streptavidin-conjugated horseradish peroxidase and with biotinylated or FITC conjugated tyramide. Finally, streptavidin conjugated with Texas Red (Southern Tech) was used to visualize the reaction. Sections in which the primary Ab was omitted served as negative control. Sections were counterstained with DAPI. Spleens fixed with 10% phosphate-buffered formaldehyde were processed and stained at the UVA Research Histology Core and then examined as described [67]. Mice given HD (2 x 105 PFU) MCMV were described previously [23]. Additional crossed mice were infected with LD (1 x 104 PFU) MCMV and analyzed as described before [23]. Spleen and liver fragments were processed for genomic (g) DNA genotyping and quantitative real-time PCR (qPCR) analysis of MCMV genome level as described [23,68]. MA/My, M. H2b and M. H2b-TgDk received HD or LD MCMV and tissues were analyzed 3-7d after infection. Spleen and liver gDNA samples were prepared using a Gentraprep kit. A DNA' footprint' of MHC, NKC and chromosome 19 genotypes was obtained for each mouse using gene-specific PCR strategies as described [23,32,69]. High-resolution melt PCR for H-2D exons 5 and 8 was used to validate cross and congenic mouse genotypes. Genome-wide genotyping was performed at DartMouse (Dartmouth, NH) using an Illumina medium density SNP panel. Mouse genotypes were concatenated and transposed in Excel. Genomic profiles were reviewed for marker alignment and reliability. Those with greater than 20% failure rate (double recombination or non-informative) were discarded. Genomic profiles were then merged with quantitative trait data and stored in a Csv file. LD- and HD-infected mice were studied and analyzed in two separate cohorts. QTL mapping was performed using R/qtl in the statistical computing program R [28]. Csv files were imported into R using read. cross. Calc. genoprob (5 cM step) was used to predict genotypes between markers. Each trait was analyzed using a One QTL scan (Haley-Knott method). QTL significance thresholds were calculated based on 10,000 permutations. QTL effect plots were generated using effectplot after running sim. geno (1 cM step) and n. draws (= 64) to obtain imputations. In some cases, the phenotype x genotype command was run to view phenotype values for individual animals at a given marker. Traits analyzed with two-dimensional genome scans (scantwo) in R/qtl had significance thresholds based on 1000 permutations. Assessment of data integrity for HD mice has been described [23]. Similar methods were used to assess and validate LD results. Statistical analyses included paired and unpaired Student T tests, Pearson correlations and multiple linear regression tests performed in R [70] with select plots drawn using the ggplot2 package [71]. Results obtained for MA/My-derived congenic strains were analyzed using one-way ANOVA in conjunction with Tukey’s test using Prism software (version 6. 0d, 2013).
Uncovering the genetic basis of resistance to viral infection and disease is critical to learning about how immune defenses might be adjusted, how to design better vaccines, and how to elicit effectual immune protection in human populations. Prior studies have shown that both MHC and non-MHC genes support host defenses, or endow specialized immune cells with efficient sensing or responsiveness to infection. Many additional resistance genes remain to be identified, including difficult to detect smaller-effect alleles, which might add to or interact with other genetic factors. Our grasp of the complex interaction involving these genetic elements is thus inadequate. We combined genomic and multiparameter phenotypic analyses to map and identify host genes that control immune cells or sensitivity to viral infection. We reasoned that some might also affect viral clearance. Thus we enumerated a range of immune cell traits in mice before and after infection, which permitted genomic analysis of viral immunity, and mapping of genetic modifiers for each trait. Our study demonstrates that distinct loci collectively regulate both NK cells and host resistance, which provides a framework to understand the genetic interactions, and a variety of potential novel targets to adjust NK cell functionality and host resistance to infection.
blood cells
body weight
medicine and health sciences
immune cells
immune physiology
quantitative trait loci
spleen
immunology
clinical medicine
physiological parameters
mammalian genomics
white blood cells
major histocompatibility complex
animal cells
genetic loci
animal genomics
immune response
cell biology
clinical immunology
nk cells
physiology
genetics
biology and life sciences
cellular types
genomics
Genomic Modifiers of Natural Killer Cells, Immune Responsiveness and Lymphoid Tissue Remodeling Together Increase Host Resistance to Viral Infection