Previous Article |
Table of Contents
| Next Article
GENETICS
Susceptibility to tuberculosis: A locus on mouse chromosome 19 (Trl-4) regulates Mycobacterium tuberculosis replication in the lungs





*Department of Biochemistry, Center for Host Resistance, and Cancer Center, McGill University, Montreal, PQ, Canada H3G 1Y6;
Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom OX3 7BN; and
Trudeau Institute, Saranac Lake, NY 12983
Communicated by Louis Siminovitch, Mount Sinai Hospital, Toronto, ON, Canada, March 25, 2003 (received for review February 1, 2003)
| Abstract |
|---|
|
|
|---|
To date, inherited mutations in gp91/phox (14), IL-12 (15), and in the IFN-
receptor genes (16) have been found in a few rare familial cases of infantile TB or of disseminated Mycobacterium bovis (bacillus CalmetteGuérin) infection (4, 17). However, such Mendelian disorders are exceedingly rare. The genetic component of TB susceptibility has been investigated in population- and family-based studies. Case control studies in areas of endemic disease have pointed to several gene variants contributing to TB risk, including the HLA (4), the natural resistance-associated macrophage protein NRAMP1 (1822), the vitamin D receptor (23, 24), and the mannose-binding protein (25). A strong association [logarithm of odds (lod) 3.8] of NRAMP1alleles on 2q35 with susceptibility to TB was independently found by linkage analysis in a large Aboriginal Canadian pedigree in the outbreak situation (26). Major gene effects were recently investigated by whole-genome scan by using 173 affected sib pairs from The Gambia and South Africa; this analysis identified suggestive linkages (lod
2) on chromosomes (Chr.) 15q and Xq (27, 28). These studies suggest that the genetic control of susceptibility to TB in humans is complex.
Such complex genetic traits can be studied in mouse models of disease, where environmental and genetic components can be best controlled, and where single gene effects may have become fixed in inbred, recombinant inbred, and recombinant congenic strains of mice (29). Genetic control of susceptibility to TB is complex in mouse and is influenced by the M. tuberculosis isolate, the route and dose of infection, the mouse strains used, and the phenotypic measure of susceptibility (3035). Using survival time after i.v. injection of 1 x 105 M. tuberculosis, inbred mouse strains are classified (36) into either highly susceptible (CBA, C3H, DBA/2, 129svJ) or highly resistant (C57BL/6, BALB/c). Susceptibility of DBA/2 (D2) and resistance of C57BL/6 (B6) and BALB/c strains, as measured by mean survival time, are also observed after airborne infection with 102 bacilli (37). Using informative backcross and F2 mice issued from C3HeB/Fe and B6 parents and infected i.v. with 106 M. tuberculosis Erdman, Kramnik et al. (38) have mapped a locus on distal Chr.1 (sst1; position 4958 cM) that controls the rate of bacterial replication and granuloma formation in the lung.
Susceptibility of D2 mice to M. tuberculosis infection is characterized by (i) progressive bacterial replication in the lung; (ii) extended neutrophil-dominated lung pathology, including large numbers of acid-fast bacilli and areas of necrosis; and (iii) early death (3639). We have previously mapped three quantitative trait loci (QTL) on distal Chr. 1 (Trl-1; lod 4.8), proximal Chr. 7 (Trl-3; lod 4.7), and proximal Chr. 3 (Trl-2; lod 3.9) that affect survival time of (B6 x D2) F2 mice after i.v. infection with 1 x 105 M. tuberculosis H37Rv (39). In this report, the genetic analysis of D2 susceptibility to pulmonary TB was expanded by using an infection protocol more closely related to the human situation than the previously used i.v. model (35, 38, 39). (B6 x D2) F2 mice were infected by the respiratory route with 2 x 102 M. tuberculosis H37Rv, and pulmonary bacterial load at 90 days was used as a quantitative measure of susceptibility. A whole-genome scan revealed a major locus on Chr. 19 (designated Trl-4) regulating replication of M. tuberculosis in the lung.
| Materials and Methods |
|---|
|
|
|---|
Mycobacteria. M. tuberculosis strain H37Rv (TMC no. 102) was obtained from the Trudeau Mycobacterial Culture Collection as a frozen (-70°C) log phase dispersed culture in Proskauer and Beck medium (Difco) containing 0.01% Tween 80. For each experiment, a vial was thawed, subjected to 5-s ultrasound to break up aggregates, and diluted appropriately in PBS containing 0.01% Tween 80. Mice (810 weeks of age) were inoculated with 2 x 102 colony-forming units (cfu) by aerosol in a Middle-brook airborne infection apparatus (Tri Instruments, Jamaica, NY). Bacilli were enumerated in the lungs of infected mice 90 days postinfection by preparing lung homogenates in PBS containing 0.05% Tween 80 and by plating 10-fold serial dilutions of the homogenates on enriched agar (Middlebrook 7H11; Difco). cfu were enumerated after 34 weeks of incubation at 37°C, and the data are presented as log10 of total cfu count per lung.
Genotyping. Before infection, tail biopsies were obtained, and genomic DNA was prepared (39). A total of 151 microsatellite markers distributed over all Chr. except Chr. Y (
10 cM coverage) were selected (www-genome.wi.mit.edu) and purchased from Research Genetics (Huntsville, AL). Genotyping was performed by standard PCR-based method by using trace amounts of [32P]
-dATP, followed by separation on denaturing polyacrylamide gels, exactly as described (39). Some markers were genotyped by using primer pairs fluorescently labeled, either commercially available or custom-synthesized by Applied Biosystems. In this case, one of the primers was synthesized and labeled with FAM, HEX, or NED phosphoramides. Products were analyzed by capillary separation by using an Applied Biosystems Prism 3700 automated DNA sequencer.
Statistical Analysis. Genome-wide interval mapping analysis between lung cfus (log10) and genetic markers for the identification of QTLs was performed by using MAPMAKER/EXP Ver. 3.0 and MAPMAKER/QTL Ver. 1.1 (40). lod scores were calculated as
2/2ln (10). Permutations of the phenotypes in this sample were conducted by using linear regression in QTL CARTOGRAPHER (41, 42) to obtain empirical significance levels at each locus (10,000 iterates, unless noted otherwise). Genome-wide significance levels were also obtained by using QTL CARTOGRAPHER, providing thresholds appropriate for this particular study. Initial linkage analyses were conducted by using a "free" model involving codominance and dominance effects, yielding
2 statistics and lod scores with 2 df. Tests of specific genetic submodels were conducted by fitting each 1-df model (dominant, recessive, and codominant) and comparing the likelihood ratios [
2(1) = 2[LL(free) - LL(nested)] (LL, log likelihood). For significant linkage regions, similar likelihood-ratio tests were conducted to test whether apparent linkage results were due to mean trait differences between sexes: heterogeneity test for sex-specific effects,
2 = 2[LL(combined)-(LL(male) + LL(female))]. Significant results in this test indicate that the linkage results differ between males and females; nonsignificant results indicate that apparent differences in lod scores are consistent with random variation.
The distribution of lung cfus (log10) in (B6 x D2) F2 mice closely followed a normal distribution (see Fig. 1). Significant differences were noted between mean cfu counts in female vs. male mice, (Table 1; male mean = 7.11, female mean = 6.57; two-tailed t = 6.92, P < 0.001) (43). Linkage to Chr. X was initially detected (lod 4.56), suggesting a possible sex effect. However, further analysis indicated no significant linkage to Chr. X in either males or females when analyzed independently (lod of 0.25 and 0.81, respectively), suggesting that the original linkage simply reflects mean trait differences between genders instead of the influence of any Chr. X loci. To compensate for sex effects, the lung cfu counts (log10) were then adjusted by subtracting the gender-specific mean from each individual to create a "sex-adjusted cfu" value. This measure allows the male and female cfu data to be analyzed jointly on the same scale, free of gender-specific influences.
|
|
| Results |
|---|
|
|
|---|
50- to 75-fold increase in lung bacterial load and some associated mortality. Histological examination shows large diffuse lesions, dominated by neutrophils and containing areas of tissue necrosis (37, 39, 44). This severe lung pathology leads to early and uniform death in this group with a mean survival time of 102 days (range: days 88126) (data not shown). The genetic control of differential lung replication in the late phase of infection was investigated. For these studies, 98 informative male and female (B6 x D2) F2 animals as well as D2 and B6 controls were infected by the respiratory route with 2 x 102 live M. tuberculosis H37Rv, and 90 days later the extent of bacterial replication in the lung was determined (Fig. 1B). The 90-day time point was chosen because (i) parental strains show clear differences in lung cfus at this time, and (ii) D2 mice begin dying beyond 90 days. In this experiment, there was a highly significant 50- to 100-fold difference in cfu counts recovered from susceptible D2 (X = 7.45; range 7.17.9) when compared with resistant B6 controls (X = 5.9; range 5.66.1). cfu counts in (B6 x D2) F2 showed a continuous distribution (log10 cfus 5.558.0) between that of resistant B6 and susceptible D2 parents (Fig. 1B) with minor deviations from normality (Table 1) (43). The mean log10 cfus in the F2 was at 6.89, a value closer to susceptible D2 than to resistant B6 controls, suggesting that susceptibility does not segregate as a recessive trait in this cross. Comparison of cfu counts in male and female (B6 x D2) F2 mice showed a clear gender effect, with females more resistant to M. tuberculosis replication (X = 6.57; range 5.557.25) than males (X = 7.11; range 6.358.0). No significant deviations from normality were observed for the full or gender-specific distributions of lung cfus (Table 1).
The raw cfu data for the total F2 cross behaves as a quantitative trait amenable to study by QTL analysis. QTL mapping was carried out by a whole-genome scan approach, using a total of 151 polymorphic markers providing an average coverage of
10 cM along each chromosome (www-genome.wi.mit.edu/cgi-bin/mouse/index) (Table 3, which is published as supporting information on the PNAS web site, www.pnas.org). The largest gap was estimated at
18.5 cM for Chr. 7 and X. Genome-wide multipoint linkage analysis was performed by using MAPMAKER/EXP Ver. 3.0, MAPMAKER/QTL 1.1 (40), and QTL CARTOGRAPHER (42). Results of this analysis are shown as multipoint lod score traces in Fig. 2 A and B, and numerical data for individual intervals are shown in Table 2. Using raw cfu (log10) data as a quantitative trait, one highly significant linkage was identified on the distal portion of Chr. 19, with maximal lod score obtained for the interval defined by markers D19Mit91 and D19Mit54 (
2 = 25.74; lod 5.59, P = 0.000005). This level of significance was not observed one time in 50,000 permutations of the data. This QTL overlaps
10 cM on Chr. 19, explains 24.1% of the total phenotypic variance in the F2 cross, and was given the designation Trl-4 (for Tuberculosis resistance locus-4). Several weaker linkages were also detected on proximal (D5Mit254;
2 = 17.42; lod 3.78, empirical P = 0.0002) and distal (D5Mit216;
2 = 15.77; lod 3.42, empirical P = 0.04) portions of Chr. 5 and Chr. 7 (D7Mit270;
2 = 14.44; lod 3.14, empirical P = 0.0005), and Chr. 10 (D10Mit194;
2 = 15.83; lod 3.4, empirical P = 0.02). However, none of these linkages reached genome-wide statistical significance at the 0.01 level.
|
|
Possible linkage was initially detected to all Chr. X markers tested (lod 3.34.6). However, further analyses failed to reveal significant linkage to Chr. X in either males or females alone when analyzed independently (lod of 0.25 and 0.81, respectively), suggesting that the original linkage reflects mean differences between genders instead of any Chr. X loci. In this gender-specific analysis, evidence for linkage to Trl-4 was stronger in males (n = 59; lod 3.8, empirical P = 0.0002) than in females (n = 39; lod 1.59, empirical P = 0.02). A likelihood-ratio heterogeneity test (see Materials and Methods) indicated that these apparent differences were not statistically significant [
2(2) = 2(25.7-(17.54 + 7.37) = 1.66, P = 0.44], suggesting that the differences in lod scores between males and females may be due to differences in sample size and random variation. To further explore gender specificity in linkage to Chr. 19 QTL (Trl-4), residual cfu values for F2 mice were examined after controlling for sex effects (Table 1; see Materials and Methods). This sex-adjusted cfu count transformation provided a means to retain all data, whereas it eliminated the effects of baseline sex differences in bacterial replication. Analysis of the sex-adjusted cfus retained solely the Chr.19 hit as significant (
2 = 18.40; lod 4.00; explaining 18% of the phenotypic variance). A result this large was not observed in 50,000 permutations of the data. Conditioning on the genotypes at each of the Chr. X markers genotyped revealed no evidence for improved linkage at any autosomal loci (45). Thus, the linkage to Trl-4 did not appear to be sex-specific.
To visualize the effect of parental alleles (D2, d; B6, b) at Trl-4 on lung cfus, F2 animals were separated according to their genotype at D19Mit54 (Fig. 3B). B6 alleles were associated with reduced bacterial replication and were inherited in an incompletely dominant fashion, with mean cfu counts (log10) of 7.25 ("d/d" homozygotes), 6.85 ("b/d" heterozygotes), and 6.55 ("b/b" homozygotes) (Fig. 3B). We have previously detected a QTL on Chr. 7 (Trl-3) that affects survival of (B6 x D2) F2 mice after i.v. injection of 1 x 105 M. tuberculosis (D7Mit270;
2 = 21.4; lod 4.66) (39). In the present study, D7Mit270 shows suggestive linkage to lung cfus after aerosol infection (Table 2; lod 3.14; 16.2% of the variance), and alignment of lod score plots from both studies (Fig. 2 B vs. C) suggests that Trl-3 may indeed affect both phenotypes. Thus, F2 animals were also separated according to their genotype for D7Mit270 (Trl-3). Results in Fig. 3C also suggest that "b" alleles at the Chr. 7 linkage are associated with reduced M. tuberculosis replication in the lung and are inherited in a codominant fashion as previously noted for their effect on survival (39).
|
Analysis of combined effects of Chr. 19 and 7 loci on lung cfus suggested an additive and very strong effect of the two QTLs. Two-loci linkage analysis yielded a lod of 10.09 for the combined QTLs, explaining 38% of the variation in raw cfus. Regression tests of lung cfus on both D7Mit270 and D19Mit54 did not indicate a significant interaction between the loci (t =-1.46, P = 0.15). However, although few animals were available for analysis, mice homozygous for "b/b" alleles at both loci were completely resistant and showed cfu counts (5.55, 5.80) in the range of B6 controls (range 5.66.1), whereas mice homozygous for "d/d" alleles were completely susceptible and showed cfu counts (7.0, 7.6, and 7.8) similar to those seen in the susceptible D2 controls (range 7.17.9) (Fig. 3 D and E). This effect was specific and was not seen when Chr. 19 haplotypes were analyzed in conjunction with any of the weaker linkages (Chr. 5 and 10) detected in this study (data not shown). Finally, the presence of permissive homozygous "d/d" haplotypes at Trl-4 largely dissipated the protective effect of "b" alleles at Trl-3.
| Discussion |
|---|
|
|
|---|
To examine genetic control of susceptibility of D2 mice, we have carried out independent genome scans in informative (B6 x D2) F2 mice. To sample different aspects of host response to infection that may be under unique or common genetic controls, we used different doses (2 x 102, 1 x 105) and routes of infection (aerosol, i.v.) and monitored different quantitative measures of susceptibility (bacterial replication in the lung, overall survival time). The major conclusions of these experiments are that (i) the genetic control of susceptibility in D2 is complex; (ii) individual QTLs affecting this trait can be mapped in this experimental setting, with four significant loci mapped to date (Trl-1 to Trl-4); and (iii) different infection models reveal different gene effects with little overlap between the different experimental protocols. In a first genome scan (39), we infected (B6 x D2) F2 mice with 1 x 105 M. tuberculosis H37Rv by the i.v. route and used time of survival (log10) as a quantitative measure of susceptibility. The genetic control was found to be complex, with two significant linkages mapping on distal Chr. 1 (Trl-1; lod 4.8) and proximal Chr. 7 (Trl-3; lod 4.7), each accounting for 21% of the phenotypic variance. A third suggestive linkage was mapped to proximal Chr. 3 (Trl-2; lod 3.9) (39). In the present study, M. tuberculosis was introduced by the respiratory route (which closely resembles the mode of infection in humans), and the extent of pulmonary replication (log10 cfus at day 90) was used as a quantitative measure of susceptibility. QTL mapping using either raw or sex adjusted cfu counts revealed a highly significant linkage on Chr. 19 (Trl-4; D19Mit91 and D19Mit54;
2 = 25.67; lod 5.58, P = 0.000005), which explains 24.1% of the total phenotypic variance in the F2 cross. The Trl-4 QTL is distinct from other QTLs previously mapped in murine models of M. tuberculosis infection, including the sst1 locus (susceptibility to tuberculosis 1) (38) originally mapped on Chr. 1 (D1Mit49) in a (C3H x C57BL/6) F2 cross, and which regulates pulmonary growth, inflammatory response, and overall survival after i.v. infection with high dose (1 x 106 i.v.) of M. tuberculosis. Trl-4 is also distinct from the QTLs mapped by Lavebratt et al. (35) on Chrs. 3, 5, 9, and 10 that regulate in a gender-specific fashion differential body weight loss after i.v. infection with high doses of M. tuberculosis (>106 cfu). Importantly, Trl-4 was identified as a locus that regulates replication of M. tuberculosis in the lung after aerosol infection with small numbers of the pathogen.
The interval for Trl-4 on Chr. 19 is
10 cM and contains an estimated 70 transcription units (data not shown). Several of these may be potential candidates by virtue of their established role in host immune and inflammatory responses. Trl-4 maps to a QTL designated Pgia12 (48), previously shown to control onset of arthritis induced by injection of human cartilage-derived proteoglycan (PGIA), a known mouse model of human rheumatoid arthritis. The relationship between Trl-4 and Pgia12 is unknown, but it is interesting to note that both QTLs seem to affect host inflammatory responses. The Trl-4 interval also contains the NF-
B (position 45.8 cM) and I

(chuk; position 45 cM) genes. NF-
B (p52) is a subunit of the NF-
B factors, a group of transcription factors implicated in the induction of numerous genes in response to inflammatory stimuli, as well as pathogen-derived or stress signals (lipopolysaccharide, IL-1, or tumor necrosis factor-
) (49). Mouse mutants lacking functional NF-
B (p52) show absence of B cell follicles in secondary lymphoid organs and cannot produce antibodies to T dependent antigens (50). NF-
B mutant mice become susceptible to Leishmania major infection, which is associated with uncontrolled parasite replication, nonhealing lesions, and failure to develop an IFN-
response (51). NF-
B mutant mice are also susceptible to Toxoplasma gondii (52). I

is one of two known kinases that phosphorylate IkBs (inhibitor of
B kinase) and thus acts as a regulatory subunit of NF-
B factors. I

mutant mice (studied in chimeras) show a phenotype similar to NF-
B mutants, with respect to impaired B cell function (53). Finally, the Trl-4 region also contains the
chain of the cell surface receptor for granulocyte/macrophage colony-stimulating factor (GM-CSF), known as GM-CSFR
(CSF2r
; position 51 cM). GM-CSF, IL-3, and IL-5 are related cytokines that bind to cell surface receptors composed of a cytokine-specific
chain and a
chain common to the three receptors (54). GM-CSF acts as a growth factor for macrophages and granulocytes, and GM-CSF mutant mice show increased susceptibility to pneumonia caused by Pasteurella pneumotropica, group B Streptococcus, Penumocytis carinii, and others (5558). They also display reduced pathogen killing by isolated alveolar macrophages (59), which is concomitant to exaggerated inflammatory response in the lungs (57). Importantly, the only "constitutive" phenotype displayed by GM-CSF and GM-CSFR
mutant mice is "pulmonary alveolar proteinosis" (PAP), a condition characterized by abnormal catabolism of lung surfactant by lung epithelial cells and by alveolar macrophages (54). Together with the recent discovery of a GM-CSFR
mutation in a human PAP patient (60), these findings suggest that GM-CSF plays a major role in lung surfactant homeostasis by these cells. Interestingly, both alveolar type II epithelial cells and alveolar macrophages are invaded in vivo by M. tuberculosis (61) and are key hosts to this pathogen early in infection.
In two genome scans conducted to date in (B6 x D2) F2 mice (scan 1: i.v. infection, survival; scan 2: aerosol, lung replication), Trl-4 constitutes the strongest linkage identified to date with a lod
5.6. Interestingly, the Trl-4 linkage was detected only in scan 2, whereas the Chr. 1 (Trl-1) and Chr. 3 (Trl-2) QTLs were detected only in scan 1. This observation suggests that Trl-4 may affect pulmonary replication per se (after aerosol infection), whereas Trl-1/Trl-2 may influence time of death in the presence of high lung bacterial load. Alternatively, Trl-1/2/4 may be statistical accidents that await validation in larger groups of mice of the same cross or a different cross. This explanation is unlikely for Trl-4, which reaches a very high degree of significance and explains a large proportion of the phenotypic variance in the cross. Thus, we believe that Trl-4 is a major determinant of M. tuberculosis replication in the lungs, after infection by the respiratory route. The Trl-3 QTL was the only QTL detected in both scans 1 (39) and 2 (this study). The observation that in both scans the Trl-3 alleles of B6 are protective and inherited in a codominant fashion suggests that the gene effect is real and thus this QTL affects both the extent of M. tuberculosis replication (scan 2) and the survival to infection (scan 1). Further-more, a major additive effect of Trl-3 and Trl-4 on pulmonary replication of M. tuberculosis (lod 10.09;
40% of variance) was detected in scan 2 (this study). Remarkably, and although only a few animals were available for study, mice homozygous for B6 alleles at Trl-3/Trl-4 were completely resistant to infection and phenotypically undistinguishable from B6 parents. Also, mice homozygous for D2 alleles at both loci were completely susceptible and similar to susceptible D2 controls (Fig. 3 D and E). The independent and combined contribution of Trl-3 and Trl-4 to regulation of M. tuberculosis replication in the lung is currently being investigated in congenic mice.
The Chr. 19 Trl-4 region is syntenic with human 10q, whereas the Chr.7 Trl-3 is syntenic with human 19q13. A possible association of these chromosomal regions with susceptibility to TB in humans can now be tested in population studies from areas where the disease is endemic.
| Acknowledgements |
|---|
| Footnotes |
|---|
To whom correspondence should be addressed. E-mail: philippe.gros{at}mcgill.ca.
| References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
J.-F. Marquis, A. Nantel, R. LaCourse, L. Ryan, R. J. North, and P. Gros Fibrotic Response as a Distinguishing Feature of Resistance and Susceptibility to Pulmonary Infection with Mycobacterium tuberculosis in Mice Infect. Immun., January 1, 2008; 76(1): 78 - 88. [Abstract] [Full Text] [PDF] |
||||
![]() |
B.-S. Yan, A. V. Pichugin, O. Jobe, L. Helming, E. B. Eruslanov, J. A. Gutierrez-Pabello, M. Rojas, Y. V. Shebzukhov, L. Kobzik, and I. Kramnik Progression of Pulmonary Tuberculosis and Efficiency of Bacillus Calmette-Guerin Vaccination Are Genetically Controlled via a Common sst1-Mediated Mechanism of Innate Immunity J. Immunol., November 15, 2007; 179(10): 6919 - 6932. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Keller, R. Hoffmann, R. Lang, S. Brandau, C. Hermann, and S. Ehlers Genetically Determined Susceptibility to Tuberculosis in Mice Causally Involves Accelerated and Enhanced Recruitment of Granulocytes Infect. Immun., July 1, 2006; 74(7): 4295 - 4309. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. V. Pichugin, S. N. Petrovskaya, and A. S. Apt H2 complex controls CD4/CD8 ratio, recurrent responsiveness to repeated stimulations, and resistance to activation-induced apoptosis during T cell response to mycobacterial antigens J. Leukoc. Biol., April 1, 2006; 79(4): 739 - 746. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. B. Kamath, J. Alt, H. Debbabi, C. Taylor, and S. M. Behar The Major Histocompatibility Complex Haplotype Affects T-Cell Recognition of Mycobacterial Antigens but Not Resistance to Mycobacterium tuberculosis in C3H Mice Infect. Immun., December 1, 2004; 72(12): 6790 - 6798. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||