Abstract
The hematopoietic stem and progenitor cell (HSPC) compartment is subject to extensive quantitative genetic variation. We have previously shown that transforming growth factor-beta 2 (TGF-β2) at low concentrations enhances flt3 ligand induced growth of HSPCs, while it is potently antiproliferative at higher concentrations. This in vitro enhancing effect was subject to quantitative genetic variation, for which a quantitative trait locus (QTL) was tentatively mapped to chr.4. Tgfb2+/- mice have a smaller and more slowly cycling HSPC compartment, which has a decreased serial repopulation capacity, and are less susceptible to the lethal effect of high doses of 5-fluorouracil (5-FU). To unequivocally demonstrate that these phenotypes can be attributed to the enhancing effect of TGF-β2 on HSPC proliferation observed in vitro and are therefore subject to mouse strain-dependent variation as well, we generated congenic mice where the telomeric region of chr.4 was introgressed from DBA/2 into C57BL/6 mice. In these mice, the enhancing effect of TGF-β2 on flt3 signaling, but not the generic antiproliferative effect of high concentrations of TGF-β2, was abrogated, confirming the location of this QTL, which we named tb2r1, on chr.4. These mice shared a smaller and more slowly cycling HSPC compartment, increased 5-FU resistance but not a decreased serial repopulation capacity with Tgfb2+/- mice. The concordance of phenotypes between Tgfb2+/- and congenic mice indicates that HSPC frequency and cycling are regulated by tb2r1, while an additional QTL in the telomeric region of chr.4 may regulate the serial repopulation capacity of HSCs.
INTRODUCTION
The hematopoietic stem and progenitor cell (HSPC) compartment, which gives rise to all cells of the blood, is regulated by cell intrinsic and extrinsic mechanisms that endow the system with a capacity to maintain the steady-state production of multiple lineages of blood cells, each with widely differing lifespans, and to provide prompt replacement of these cells after hematopoietic stress (1-5).
A powerful approach to dissect the multiple pathways involved in the regulation of self-renewal, differentiation and proliferation of HSPCs is the investigation of genetic variation in the HSPC compartment by quantitative trait analysis. Quantitative traits vary continuously across genetically different individuals and are inherited in a non-mendelian fashion because of the contribution of multiple loci to the phenotype (6). These loci are called quantitative trait loci, or QTL. Several genes and pathways that show quantitative genetic variation in humans have been successfully modeled in mice (6). Importantly, many disease susceptibility QTL in mice could be translated to humans (7-12). Therefore, the investigation of quantitative variation in the hematopoietic system of mice may not only lead to the identification of novel regulatory mechanisms, but may also provide insight into susceptibility to hematological diseases and toxicity of drugs in the hematopoietic system.
Among inbred mouse strains there is extensive genetically determined variation in the function and kinetics of HSPCs (5, 13-23). Although multiple suggestive QTL have been mapped, for only one has the underlying gene been identified thus far, latexin, which is involved in the regulation of stem cell pool size (24). However, its mechanism of action is yet not understood. We have shown that signaling by the transforming growth factor-beta (TGF-β) isoform TGF-β2 plays a role in mouse strain-dependent variation in the HSPC compartment (25, 26). There are 3 isoforms of TGF-β (TGF-β1, -β2 and -β3), which are encoded on different chromosomes, but bind to the same receptors though with different binding mechanisms and relative affinities (27-30). TGF-βs are potent inhibitors of HSPC proliferation in culture (31-34), although studies in mice with conditional deletion of the type I TGF-β receptor questioned the importance of these observations as these mice had normal steady-state hematopoiesis (35, 36). Our data suggest, however, that one of these isoforms, TGF-β2, is in fact a positive regulator of HSC both in vitro and in vivo. In vitro, TGF-β2 has a biphasic dose response on the proliferation of purified HSPCs, defined as lineage-Sca1+c-kit+ or LSK cells (37). At low concentrations, the effect of this factor is stimulatory and requires factor(s) present in mouse and fetal calf sera (26), while at higher concentrations inhibition of proliferation occurs. This stimulatory effect of TGF-β2 is caused by a specific interaction with flt3 signaling, and not with any other early-acting hematopoietic cytokine (submitted for publication). Studies in mice with a heterozygous deletion of Tgfb2 (Tgfb2-/- mice die at birth (38)) revealed that the frequency of LSK cells as well as their cycling activity and the serial repopulating capacity of HSC were lower than in wildtype (wt) littermates (25). Interestingly, while the stimulatory effect of exogenously added TGF-β2 on the proliferation of HSPCs was relatively subtle and only detected over a limited concentration range, Tgfb2+/- LSK cells showed a profound proliferation defect in response to flt3 ligand (flt3L; submitted for publication) or to growth factor combinations containing flt3L (25). These data strongly suggest that the stimulatory effect of TGF-β2 is predominantly, though not exclusively cell autonomous, and is based on a specific enhancement of flt3 signaling. This effect of TGF-β2 on HSPCs is subject to genetically determined variation in inbred mouse strains. Using BXD recombinant inbred stains of mice a suggestive QTL for this trait was mapped to the telomeric region of chr.4 (25). This region also contains a QTL that contributes to genetic variation in the frequency of LSK cells (20). While it is probable that the phenotypes of Tgfb2+/- mice can be attributed to the enhancing effect of TGF-β2 on HSPC proliferation observed in vitro, and are therefore also likely subject to mouse strain-dependent variation, definitive proof of this contention requires the generation of mouse strains where specifically the proliferative effect of TGF-β2 was absent.
A classical strategy to demonstrate the veracity of a mapped QTL is the construction of congenic mice, where a chromosomal region from one inbred mouse strain is introgressed into the genetic background of a different mouse strain, leading to a change in the phenotype of the acceptor strain for this trait. In this context, we generated appropriate congenic mice, which allowed us to confirm the location of a QTL on chr.4 that regulates TGF-β2 responsiveness of HSPCs. In addition, we show here that these congenic mice shared several phenotypes with Tgfb2+/- mice, including a decreased responsiveness to flt3L, a more slowly cycling and smaller HSPC compartment and decreased lethality from cell cycle-specific cytotoxic drug 5-fluorouracil (5-FU). Together, our findings strongly suggest that genetic variation in flt3 signaling caused by variation in isoform-specific TGF-β2 signaling mechanistically underlies these potentially clinically relevant quantitative traits. Another trait, decreased serial repopulation capacity of HSC, however, was not shared between congenic and Tgfb2+/- mice, suggesting the presence of distinct QTL in this region of chr.4 that are involved in the regulation of the repopulation capacity of HSC.
MATERIALS AND METHODS
Mice
4-6-week-old C57BL/6J mice were purchased from Jackson Laboratories (Bar Harbor, ME). B6.D2-chr.4 mice were generated through the Speed Congenic Service of the Jackson Laboratory (Bar Harbor, ME), using 150 single nucleotide polymorphism (SNP) markers (39, 40). After the initial C57BL/6 and DBA/2 mating, 80-160 animals of each generation after the F1 crossbreeding were scanned using 108 Mit markers, five of which, d4mit155, d4mit37, d4mit251, d4mit234 and d4mit190, were on chr.4. Of these, 20-25 samples with the telomeric chr.4 genotype of DBA/2 were then scanned using 150 SNP markers to select the animals with the highest percent of C57BL/6 in the remaining part of the genome and mated those to C57BL/6 mice to produce the next generation. After six generations (N6) animals were found to be 100% C57BL/6 for all SNP markers except in the telomeric 20cM of chr.4 which was DBA/2. N6 animals were mated to produce homozygous DBA/2 chr.4 telomeric region at generation N6F1. The sex chromosomes were fixed via mating. Animals were housed in a specific pathogen-free facility. Experiments and animal care were performed in accordance with the Mount Sinai Institutional Animal Care and Use Committee (IACUC).
Antibodies and cytokines
Unconjugated anti-CD2, -CD3ε, -CD8α, -CD4, -B220, -Gr1, -Mac1, PE-conjugated anti-CD45.1, biotinylated anti-Thy1 and FITC-conjugated goat anti-rat antibody were purchased from Southern Biotechnologies (Birmingham, AL). FITC-conjugated anti-CD2, -CD3ε, -CD8α, -CD4, -CD19, -B220, -Gr1, and -Mac1, PE-conjugated anti-flt3, PECy7-conjugated streptavidin and APC-AlexaFluor®750-conjugated anti-c-kit were purchased from eBiosciences (San Diego, CA). Unconjugated anti-Ter119, biotinylated anti-Sca1 and -CD34, PE-conjugated anti-Sca1 and -CD34, PerCP-Cy5.5-conjugated streptavidin and anti-CD45.2, APC-conjugated anti-c-kit and goat anti-rat antibody, PerCP-conjugated streptavidin, PECy7-conjugated anti-CD19, APC-Cy7-conjugated streptavidin and anti-CD19 were purchased from Pharmingen (San Diego, CA). Pacific blueconjugated anti-Sca1 were purchased from BioLegend (San Diego, CA). Recombinant cytokines were purchased from R&D Systems (Minneapolis, MN). Recombinant human FLT3 ligand was received from Amgen (Seattle, WA).
Cell sorting and flow cytometry
Bone marrow (BM) cells were prepared by flushing the femora and tibia of mice with cold DMEM (Cellgro, Mediatech, Manassas, VA) containing 2% FBS (Hyclone, Logan, UT) and penicillin/streptomycin (Cellgro). Mononuclear cells were obtained after gradient centrifugation using lymphocyte separation medium (Cellgro). Low-density BM cells were stained with antibodies for lineage antigens (CD2, CD3ε, CD8α, CD4, CD19, B220, Ter119, Gr1, Mac1), Sca1 and c-kit, and isolated by cell sorting, as shown in Figure S1, using FACSVantage SE (Beckon Dickinson, San Diego, CA), Cytopeia (Advanced Cytometry Systems, Seattle, WA) or MoFlo (CytoMation, Fort Collins, CO) sorters, to obtain lineage-Sca1+c-kit+ (LSK) cells. Flow cytometric analysis was performed on a three-laser LSRII or a Special Order five-laser LSRII with DiVa software (Beckon Dickinson). Data were analyzed using FlowJo software. Doublets were excluded by plotting with forward scatter pulse area versus side scatter both for sorting and FACS analysis.
Culture of lineage -Sca1+c-kit+ (LSK) cells
Sorted LSK cells were cultured in triplicate at 20 to 60 cells per well in flat bottom 96-well plates in StemPro34 medium (Invitrogen), 10% FCS (Hyclone, Logan, UT), penicillin/streptomycin, in the presence of various cytokines, as mentioned for each experiment. Within an hour after plating, the exact number of cells per well was determined by visually counting the cells at 200x magnification. After 5 days of liquid culture in a humidified incubator with 5% CO2 at 37°C the cells were again counted.
Cell cycle analysis of LSK cells
BM mononuclear cells were stained with lineage markers (CD2, CD3ε, CD8α, CD4, B220, Gr1, Mac1, Ter119, unlabeled antibodies), followed by goat anti-rat polyclonal IgG-FITC, PE-conjugated anti-Sca1 and APC-conjugated anti-c-kit. The cells were fixed in 1% paraformaldehyde (Sigma, St Louis, MO) and 0.2% Nonidet P40 (Sigma) for 1 hour at 4°C. Fixed cells were washed in PBS and labeled with Hoechst 33342 (Sigma), final concentration of 5μg/mL, for 1 hour at 37°C. Analysis was performed on a triple laser LSRII flow cytometer with DiVA software (Becton Dickinson). Double exclusion using pulse shape was done by plotting area versus height of the UV-excited Hoechst fluorescence. Cell cycle analysis was done using up to 3000 singlet LSK cells.
5-Fluorouracil administration
5-Fluorouracil (5-FU, Sigma) was dissolved in PBS at concentration of 20mg/mL for survival and of 11.25mg/mL for recovery studies, and filtered through 0.2μm SFCA filters (Corning). The effect of sublethal doses of 5-FU (150mg/kg, IP) was followed in 6 mice per experiment of which 3 were bled alternately every other day. Blood counts were performed using a Beckman-Coulter Ac·Tdiff automated hematology machine (Beckman Coulter, Inc., Miami, Florida, USA).
Competitive repopulation assays
LSK cells from donor mouse strain (C57BL/6 or B6.D2-chr.4, CD45.2+) were injected into lethally (700cG followed by 500 cG 3 hours later) irradiated recipient strain (C57BL/6, CD45.1+). Donor cells were mixed with 2·105 bone marrow cell from CD45.1+CD45.2+ C57BL/6 mice in case of primary competitive repopulation transplantations. Peripheral blood cells were analyzed for the expression of CD45.1, CD45.2, and lineage antigens (Thy-1, Gr-1/Mac1, and CD19) 12 weeks after transplantation. Competitive serial transplantations were performed at least 16 weeks after reconstitution by injecting 2·106 bone marrow cells from primary recipients into lethally irradiated secondary recipients (CD45.1+). Three secondary recipients were reconstituted with bone marrow from one primary donor (n=9 primary donors per donor strain). After 3 months, the ratio between the competing cell populations was measured in the peripheral blood of the secondary recipients (n=27 total recipients per donor strain). The data were presented as ratios between the two donor populations determined by the expression CD45 alleles. Changes in reconstitution ratios between primary and secondary were analyzed by evaluating the difference in the logarithm of the reconstitution ratios in the respective recipients. The reason for this strategy was described previously (25). Briefly, the change in the reconstitution ratio between primary and secondary recipients is a better measure of a shift in reconstitution capacity upon serial transplantation than a change in the percentage contribution of CD45.2+ cells between primary and secondary recipients, as the latter method does not weigh the same percent change at high or low contribution level equally (i.e., a 5% change from 20 to 25% is a change in reconstitution ratio from 0.25 to 0.33, whereas a 5% change from 90 to 95% represents a much larger change in reconstitution ratio from 9 to 19). A more accurate measure is the difference of CD45.2+(donor)/CD45.1+CD45.2+(competitor) ratios in primary recipients and secondary recipients. With ratiometric data, the difference between log(CD45.2+/CD45.1+CD45.2+) can be used which added the advantage of log transformation where a ratio smaller than one will give a negative value, and negative ratios will extend over the same numerical ranges as positive ones (e.g., a ratio of 0.01 gives a log ratio of 2, a ratio of 100 gives a log ratio of 2). Thus, serial repopulation data are presented as the difference between log(CD45.2+/CD45.1+CD45.2+) of primary transplantation recipients and the log(CD45.2+/CD45.1+CD45.2+) of secondary transplantation recipients, and will be referred to as Δ log ratio.
Statistical analysis
Student’s 2-tailed t-test for paired samples was used in the calculation of all P values unless indicated otherwise. Data represent mean ±SD.
RESULTS
A locus on chr.4 regulates the stimulatory effect of TGF-β2 on HSPC growth, number and cycling, and 5-fluorouracil lethality
The isoform-specific proliferative effect of TGF-β2 on HSPCs is subject to genetically determined variation in inbred mouse strains, and a suggestive QTL for this trait was mapped to the telomeric region of chr.4 (25). To confirm the location of this QTL, we constructed congenic mice by introgressing 20cM of the telomeric region of chr.4 from DBA/2 into C57BL/6 mice. Although there was wide variation in TGF-β2 responsiveness of HSPCs among various inbred mouse strains, the TGF-β2 dose response was similar in DBA/2 and C57BL/6 mice (26). Constructing congenic mice from these strains may therefore seem ill justified at first sight. However, this QTL was mapped using BXD recombinant inbred mice, which are derived from the allele pool of C57BL/6 and DBA/2 mice. In BXD recombinant inbred strains, where the genome is made up of a patchwork of segments that are homozygously inherited from either progenitor strain, significant variation in TGF-β2 responsiveness was observed (25,26). This phenomenon is explained by the fact that if multiple loci are involved in determining a phenotype, some BXD strains can accumulate predominantly ‘high’ or ‘low’ alleles for this trait and acquire more extreme phenotypes than either of the progenitor strains. In other BXD strains and in the progenitor C57BL/6 and DBA/2 strains, the phenotypic effects of ‘high’ and ‘low’ alleles balance each other. These observations therefore proved that TGF-β2 responsiveness is a quantitative trait, and justified generating congenic mice. As only one or a limited number of QTL are transferred from the DBA/2 to the C57BL/6 background in the congenic mice, we anticipated that if the QTL mapping was correct, the phenotype of the congenic mice and the background C57BL/6 strain would differ. Using C57BL/6 as the acceptor strain also facilitates analysis of in vivo transplantation assays by taking advantage of allelic variation at the CD45 locus to track donor and host contribution to hematopoiesis.
Construction of congenic mice was accomplished by repeated backcrossing of DBA/2 onto C57BL/6 mice and selecting offspring where the telomeric region of chr.4 was heterozygous (see Material and Methods). These congenic mice will be referred to as B6.D2-chr.4 mice hereafter. The transition between the C57BL/6 and DBA/2-derived genome occurred between D4Mit37 (57 cM) and D4Mit251 (66cM) or between chromosome base pairs 04-114064127 and 04-135804867 (Figure 1A). Any difference in the effect of TGF-β2 on the proliferation of HSPCs between B6.D2-chr.4 and C57BL/6 mice must be caused by one or more alleles in the introgressed region of chr.4. If this is the case and if the antiproliferative effect of higher concentration of TGF-β2 is similar in B6.D2-chr.4 and C57BL/6 mice, then any phenotype shared by B6.D2-chr.4 mice and Tgfb2+/- mice can with a very large likelihood be assigned to the stimulatory effect of low concentrations of TGF-β2 on HSPC proliferation. To test this, we cultured LineagenegSca+c-kit+ (LSK) cells (Figure S1) from C57BL/6 and B6.D2-chr.4 mice in serum-containing media for 5 days in the presence of flt3 ligand (flt3L), kit ligand (KL), thrombopoietin (TPO) and increasing concentrations of TGF-β2. While the dose response of exogenously added TGF-β2 on the proliferation of LSK cells from parental C57BL/6 mice was biphasic with a stimulatory effect at low concentrations, this stimulatory effect was absent in B6.D2-chr.4 LSK cells (Figure 1B). Importantly, the inhibitory effect of higher concentrations of TGF-β2 (Figure 1B) and the dose response of TGF-β1 (data not shown) on the proliferation of LSK cells were similar in B6.D2-chr.4 and C57BL/6 mice. We have shown that TGF-β2 specifically enhances flt3 signaling in HSPCs (submitted for publication). Similar to Tgfb2+/- mice, in the absence of any exogenously added TGF-β2, the response of LSK cells to increasing concentrations of flt3L was approximately 50% lower in B6.D2-chr.4 than in C57BL/6 mice (Figure 1C). In contrast, the response to IL-3 was similar in both mice (Figure 1D). Together, these data indicate that only the predominantly, though not exclusively cell autonomous enhancing effect of TGF-β2 on flt3L responsiveness of HSPCs and not the generic antiproliferative effect of TGF-β2 are regulated by a QTL on chr.4. These results in fact suggest that most genetic variation in the stimulatory effect of TGF-β2 is regulated by this QTL. As our observation provides unequivocal confirmation for the location of this QTL, we call this QTL, tb2r1 (for TGF-β2 responsiveness 1).
Figure 1. Tb2r1 is located in the telomeric region of chr.4.
(A) Schematic diagram of generation of B6.D2-chr.4 mice. White portion of chr.4 represents the 20cM telomeric region of DBA/2 introgressed into C57BL/6. (B) TGF-β2 responsiveness of C57BL/6 and B6.D2-chr.4 LSK cells in vitro, in the presence of flt3L, KL and TPO, each at 50ng/mL (n=7, *P < 0.00005). (C) B6.D2-chr.4 and C57BL/6 LSK cell proliferation in response to flt3L only in the presence of 10% serum (n=5, *P < 0.02). (D) B6.D2-chr.4 and C57BL/6 LSK cell proliferation in response to IL-3 only in the presence of 10% serum (n=2, P > 0.3). 100% response is defined as that of C57BL/6 LSK cells at 50ng/mL of flt3L or IL-3.
As Tgfb2+/- mice display a hematopoietic defect, we further analyzed the hematopoietic phenotype of the B6.D2-chr.4 mice. We had previously mapped a QTL regulating the frequency of LSK cells to the telomeric region of chr.4 (20). Consistent with these mapping data, the frequency of LSK cells was significantly lower in B6.D2-chr.4 than in C57BL/6 mice (Figure 2A, Table I), confirming that the same region on chr.4 also contains a QTL regulating LSK frequency. This finding cannot be explained by intrinsic variation in c-kit expression, as the overall c-kit expression level in the bone marrow (BM) was similar in C57BL/6 mice and congenics. The frequency of a subpopulation of LSK cells containing the more primitive stem cells, the CD34-flt3-LSK population, was also decreased in B6.D2-chr.4 mice BM (Figure 2B). BM cellularity, peripheral blood counts and spleen weight were similar in B6.D2-chr.4 and C57BL/6 mice (Table I). Tgfb2+/- mice have a more slowly cycling HSPC compartment than wt littermates (25). Therefore, we examined cell cycle activity of LSK cells in B6.D2-chr.4 mice. Similar to Tgfb2+/- mice, fewer LSK cells in B6.D2-chr.4 mice than in control C57BL/6 mice were in cell cycle, as measured by staining fixed BM cells with DNA-binding dye Hoechst 33342 (Figure 2C), while the fraction of S/G2/M phase cells in the total BM cell population did not differ between B6.D2-chr.4 and C57BL/6 mice (data not shown). We hypothesized that a slower cycling HSPC compartment may render Tgfb2+/- and B6.D2-chr.4 mice more resistant to the effects of cell cycle specific cytotoxic agents, such as 5-FU. The lethal effects of this drug are due to hematopoietic failure as mice are rescued from high doses of 5-FU by BM transplantation (41). Consistent with this notion, the survival of both Tgfb2+/- and B6.D2-chr.4 mice after an intravenous (IV) injection of highly lethal dose of 5-FU (500mg/kg, IV) was better than that of control mice (P < 0.0001 and P=0.0001, respectively; Figure 3A). However, recovery from a sublethal dose of 5-FU was similar in B6.D2-chr.4 and C57BL/6 mice (Figure 3B), and in wt and Tgfb2+/- mice (Figures 3C). We conclude that B6.D2-chr.4 mice and Tgfb2+/- mice share several hematopoietic phenotypes, including the frequency and cycling activity of HSPCs in vivo, and their response to flt3L in vitro. Furthermore, both mouse strains show enhanced resistance to the lethal effects of 5-FU. These observations suggest that these phenotypes of B6.D2-chr.4 mice are likely due to regulation of TGF-β2 responsiveness by tb2r1.
Figure 2. B6.D2-chr.4 mice have a smaller and slower cycling HSPC compartment.
(A) Representative flow cytometric analysis of LSK frequency in C57BL/6 and B6.D2-chr.4 mice. (B) Frequencies of CD34-flt3-LSK, CD34+flt3-LSK and CD34+flt3+LSK cells in BM of C57BL/6 and B6.D2-chr.4 mice (n=8, *P < 0.001). (C) Fraction of LSK cells in S/G2/M phase of cell cycle in C57BL/6 and B6.D2-chr.4 mice as determined by Hoechst 33342 staining of fixed BM cells (n=4, P=0.039). Lines connect data obtained from each individual experiment, where 2 or 3 mice of each genotype were pooled. P value calculated using 2-tailed Student’s t-test for four paired experiments.
Table I.
Hematopoietic parameters in C57BL/6 and B6.D2-chr.4 mice. Data presented as mean±SD. P values calculated using Student’s t-test for paired samples. n≥3 for each group for all parameters.
C57BL/6 | B6.D2 -chr.4 | P value | |
---|---|---|---|
WBC (×103/μL of peripheral blood) | 1.310 ±0.321 | 1.263 ±0.203 | 0.798 |
RBC (×106/μL of peripheral blood) | 8.387 ±0.240 | 7.653 ±0.357 | 0.148 |
Platelet (×105/μL of peripheral blood) | 12.560 ±0.442 | 11.400 ±1.310 | 0.352 |
Spleen weight (mg per g of mouse) | 3.316 ±0.511 | 3.382 ±0.660 | 0.419 |
LineagenegSca+c-kit- cell frequency in BM | 1.295 ±0.454 | 0.962 ±0.415 | 0.054 |
Bone marrow cellularity (×107 per mouse) | 4.34 ±1.15 | 4.90 ±1.89 | 0.242 |
LSK cell frequency in bone marrow | 0.260 ±0.091 | 0.088 ±0.091 | 0.000032 |
Figure 3. Decreased 5-FU lethality in B6.D2-chr.4 and Tgfb2+/- mice.
(A) Analyses of the survival of C57BL/6 and B6.D2-chr.4 mice (n=26 per group, P=0.0001), and of Tgfb2+/- and wt littermate mice (n=20 and n=23, respectively, P < 0.0001) after administration of 5-FU (500mg/kg, IV). P values calculated using Kaplan-Meier statistics. (B) Peripheral leukocyte (left panel) and platelet (right panel) counts after injection of a sublethal dose of 5-FU (150mg/kg, IP) in C57BL/6 and B6.D2-chr.4 mice (n=3 for each data point, P > 0.1). (C) Peripheral leukocyte (left panel) and platelet (right panel) counts after injection of a sublethal dose of 5-FU (150mg/kg, IP) in Tgfb2+/- mice and wt littermates (n=3 for each data point, P > 0.1).
Decreased serial reconstitution capacity is not shared between B6.D2-chr.4 and Tgfb2+/- mice
HSC from Tgfb2+/- mice show a defect in serial repopulation capacity compared to HSC from wt mice (25). Therefore, we tested the serial competitive repopulation capacity of LSK cells from B6.D2-chr.4 mice. 500 purified CD45.2+ LSK cells from C57BL/6 or B6.D2-chr.4 mice were competed with 2·105 CD45.1+CD45.2+ C57BL/6 bone marrow cells in lethally irradiated CD45.1+ C57BL/6 recipients. In the primary transplantation, LSK cells from B6.D2-chr.4 mice competed equally well compared to LSK cells from C57BL/6 mice as determined by the donor CD45.2+ contribution to peripheral blood mononuclear cells after 12 weeks (P > 0.1; Figure 4A). This was not surprising, as the content of most primitive CD34-flt3- stem cell fraction was the same in the LSK population of C57BL/6 and B6.D2-chr.4 mice (Figure S2). For serial transplantation, we injected 2·106 bone marrow cells from competitively repopulated primary recipients into lethally irradiated secondary CD45.1+ recipients, 4 months after the primary transplantation. 12 weeks later, the contribution of CD45.2+ and CD45.1+CD45.2+ cells to hematopoiesis in peripheral blood was assessed in secondary recipients. The data are presented as the logarithm of the ratio of the two donor populations determined by the expression CD45 alleles. The difference in the log reconstitution ratios (Δ log ratio; see Materials and Methods) between primary and secondary recipients was significantly higher for LSK cells from B6.D2-chr.4 mice than for LSK cells from C57BL/6 mice (Figure 4B). These data indicate that HSC from B6.D2-chr.4 mice performed significantly better than HSC from C57BL/6 mice in serial transplantation assays. This phenotype of B6.D2-chr.4 mice is not shared with Tgfb2+/- mice, which show a subtle defect in competitive repopulation capacity that increases with serial transplantation.
Figure 4. Serial competitive reconstitution capacity of HSCs from B6.D2-chr.4 and C57BL/6 mice.
(A) Log ratio of donor/competitor contribution in primary and secondary transplantations of C57BL/6 and B6.D2-chr.4 LSK cells (donor, CD45.2+), and C57BL/6 competitor bone marrow cells (CD45.1+CD45.2+) into lethally irradiated recipients (CD45.1+ C57BL/6) (n=5, P > 0.1, primary transplantation; n=9, P=0.019, secondary transplantation). P values calculated using Mann-Whitney test. (B) The difference in the log ratio, Δ log ratio, of reconstitution ratios between secondary and primary transplantations for recipients of bone marrow from primary recipients repopulated with LSK cells from C57BL/6 or B6.D2-chr.4 mice (n=9, P=0.042). P value calculated using 2-tailed Student’s t-test for unpaired samples.
DISCUSSION
We show here that a strong QTL in the telomeric region of chr.4, tb2r1, regulates the enhancing effect of TGF-β2 on HSPC proliferation in vitro. We have shown that this effect of TGF-β2 is caused by a specific enhancement of flt3 signaling (submitted for publication). As in B6.D2-chr.4 mice the positive regulatory effect of TGF-β2 is absent while the antiproliferative effect of high concentrations of TGF-β2 was similar to C57BL/6 mice, any traits shared between B6.D2-chr.4 and Tgfb2+/- mice are most likely caused by a defect in TGF-β2-mediated enhancement of flt3 signaling. These traits are the frequency and cycling activity of HSPCs in vivo and 5-FU lethality. Alternatively, it is possible, though unlikely, that these hematopoietic phenotypes of B6.D2-chr.4 mice are not caused by a defect in TGF-β2 signaling, though they must be caused by reduced TGF-β2 signaling in Tgfb2+/- mice. In that case, this narrow region of chr.4 contains a remarkable number of QTL that independently affect hematopoiesis.
The lower cycling activity of HSPCs in both Tgfb2+/- and B6.D2-chr.4 mice compared to control mice is likely linked to the decreased lethality from 5-FU in these mouse strains, as it has been shown that mice die of hematopoietic failure after high doses of 5-FU (41). As the cycling activity of total bone marrow was similar in Tgfb2+/- and B6-D2-chr.4 mice compared to their appropriate control groups, our data suggest that the lethality from high doses of 5-FU is due to loss of the more primitive cells in the HSPC compartment. On the other hand, recovery from a sublethal dose, which did not appear to be subject to genetically determined variation in the mouse strains tested here, is likely dependent on more mature progenitors that are less affected by TGF-β2 signaling. The similar rate of recovery from sublethal 5-FU dose in C57BL/6, Tgfb2+/- and B6.D2-chr.4 mice also indicates that 5-FU metabolism is similar in these mouse strains. It cannot be fully excluded that additional QTL in the introgressed region of chr.4 regulated lethality from 5-FU, for example by affecting intestinal toxicity. A QTL regulating lethality from 5-FU is clinically important, however. As genes and pathways that show quantitative genetic variation in mice often also do so in humans (6-12), it is possible that similar variation exists in humans, and that the underlying mechanism is the same.
The only phenotype not shared between Tgfb2+/- and congenic mice was the serial repopulating capacity of HSC, which was enhanced in the latter but decreased in the former. Traits not shared between B6.D2-chr.4 mice, where at least 174 alleles of genes with known protein product might be distinct (Table SI), and Tgfb2+/- mice, where only one allele of one gene is deleted, may be caused by additional QTL in the introgressed region. One explanation for this observation may therefore be that a separate QTL on chr.4 is responsible for the enhanced serial repopulation capacity of HSC from B6.D2-chr.4 mice. In particular, a stronger QTL in this region may override the effect of tb2r1 in regulating serial repopulation capacity. However, as in B6.D2-chr.4 mice only the stimulatory TGF-β2 signaling is compromised, while in Tgfb2+/- mice both stimulatory and inhibitory TGF-β2 signaling are equally affected, another explanation for this discrepancy may be that the defect in serial repopulation capacity of Tgfb2+/- HSC is due to the lower level of antiproliferative TGF-β2 signaling in these mice, which might lead to accelerated exhaustion of HSC. Only gene identification of tb2r1, for which the generation and analysis of the congenic mice described here is a critical step, will shed light on this issue.
TGF-β2 responsiveness of LSK cells was similar in DBA/2 and C57BL/6 mice. It may therefore seem unexpected that in congenic mice, where a region of chr.4 was introgressed from DBA/2 to C57BL/6, responsiveness to the flt3-specific proliferative effect of TGF-β2 was abrogated. However, TGF-β2 responsiveness is a quantitative trait, and, as discussed previously (see Results), there is likely an unbalanced contribution of ‘high’ and ‘low’ responsiveness alleles in B6.D2-chr.4 mice. In addition, the effect of certain alleles can depend on the genetic background, leading to unexpected epistatic interactions (42). Hence, the isolation of a DBA/2 allele on a C57BL/6 background may alter the balance between ‘high’ and ‘low’ alleles and/or change allele-specific epistatic interactions, leading to the observed near complete abrogation of the phenotype in B6.D2-chr.4 mice (42). The phenotype of the B6.D2-chr.4 mice was still somewhat unexpected, however, as according to our quantitative trait analysis DBA/2 mice carried the ‘high’ allele for LSK frequency and TGF-β2 responsiveness on chr.4 (26). Consequently, a higher frequency of LSK cells and a higher responsiveness to TGF-β2 were anticipated in B6.D2-chr.4 mice. Surprisingly, exactly the opposite was observed. Two explanations are possible. One possibility is that the telomeric region of chr.4 contains more than one QTL with opposing effects, and that the introgressed 20cM contains only one of those, the other one(s) lying just outside of the terminal 20cM of chr.4. Another possibility is that the region only contains one QTL, but, because of epistatic interactions with other genes, has a different effect on a C57BL/6 than on a DBA/2 background. The first scenario, implying multiple QTL within the same region, is surprisingly common (43), while evidence in plants supporting the plausibility of the second explanation has been published recently (44).
Identifying the gene underlying tb2r1 is a critical step towards undestanding the mechanism of the effect of TGF-β2 on HSPCs. There are more than 170 known gene products in the introgressed portion of chr.4 in B6.D2-chr.4 mice (Table SI). Some of these, such as Ski (45) and Prdm16 (46), have been described to play a role in TGF-β signaling via their interaction with Smad4 and Smad3, respectively. Others, like Tnfrsf9 (47) and Tnfrsf14 (48), have been implicated in hematopoiesis by knockout mouse models. However, speculating on potential candidates could be misleading as the gene underlying QTL can be very surprising. For example, latexin, regulating HSC pool size in mice (24), is an inhibitor of metalocarboxypeptidases in neurons. The mechanism of its role in the biology of HSC is still unknown. HPE, the hereditary hemochromatosis protein in humans, is highly homologous to major histocompatibility complex (MHC) class I proteins but has no immunological function (49), and is important in iron intake regulation (50). Hence, we have adopted an unbiased approach to achieve gene identification through gene expression profiling in C57BL/6 and B6.D2-chr.4 mice, and through positional cloning using fine-congenic lines with smaller introgressed regions that have been generated from B6.D2-chr.4 mice.
In conclusion, we have unequivocally shown that responsiveness of HSPCs to the proliferative effect of TGF-β2 is regulated by the QTL tb2r1 on chr.4. Given the concordance of phenotypes between Tgfb2+/- and B6.D2-chr.4 mice, our observations indicate that frequency and cycling of HSPCs, in addition to lethality from high doses of 5-FU are also quantitative traits regulated by tb2r1. Finally, an additional QTL in the telomeric region of chr.4 may regulate the serial repopulation capacity of HSCs.
Supplementary Material
ACKNOWLEDGEMENTS
SA performed most of the experiments and wrote the manuscript. LG performed the initial analysis of B6.D2-chr.4 mice. JC performed the 5-FU survival studies on Tgfb2+/- mice. HWS designed and supervised the experiments and co-wrote the manuscript with SA.
Supported by a grant from the NIH RO1 AG016327 and HL073760.
REFERENCES
- 1.Orkin SH, Zon LI. Hematopoiesis: an evolving paradigm for stem cell biology. Cell. 2008;132:631–644. doi: 10.1016/j.cell.2008.01.025. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 2.Harrison DE, Lerner CP. Most primitive hematopoietic stem cells are stimulated to cycle rapidly after treatment with 5-fluorouracil. Blood. 1991;78:1237–1240. [PubMed] [Google Scholar]
- 3.Ogawa M. Differentiation and proliferation of hematopoietic stem cells. Blood. 1993;81:2844–2853. [PubMed] [Google Scholar]
- 4.Szilvassy SJ, Cory S. Phenotypic and functional characterization of competitive long-term repopulating hematopoietic stem cells enriched from 5-fluorouracil-treated murine marrow. Blood. 1993;81:2310–2320. [PubMed] [Google Scholar]
- 5.de Haan G, Van Zant G. Intrinsic and extrinsic control of hemopoietic stem cell numbers: mapping of a stem cell gene. J Exp Med. 1997;186:529–536. doi: 10.1084/jem.186.4.529. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 6.Peters LL, Robledo RF, Bult CJ, Churchill GA, Paigen BJ, Svenson KL. The mouse as a model for human biology: a resource guide for complex trait analysis. Nat Rev Genet. 2007;8:58–69. doi: 10.1038/nrg2025. [DOI] [PubMed] [Google Scholar]
- 7.McIntire JJ, Umetsu SE, Akbari O, Potter M, Kuchroo VK, Barsh GS, Freeman GJ, Umetsu DT, DeKruyff RH. Identification of Tapr (an airway hyperreactivity regulatory locus) and the linked Tim gene family. Nat Immunol. 2001;2:1109–1116. doi: 10.1038/ni739. [DOI] [PubMed] [Google Scholar]
- 8.Ueda H, Howson JM, Esposito L, Heward J, Snook H, Chamberlain G, Rainbow DB, Hunter KM, Smith AN, Di Genova G, Herr MH, Dahlman I, Payne F, Smyth D, Lowe C, Twells RC, Howlett S, Healy B, Nutland S, Rance HE, Everett V, Smink LJ, Lam AC, Cordell HJ, Walker NM, Bordin C, Hulme J, Motzo C, Cucca F, Hess JF, Metzker ML, Rogers J, Gregory S, Allahabadia A, Nithiyananthan R, Tuomilehto-Wolf E, Tuomilehto J, Bingley P, Gillespie KM, Undlien DE, Ronningen KS, Guja C, Ionescu-Tirgoviste C, Savage DA, Maxwell AP, Carson DJ, Patterson CC, Franklyn JA, Clayton DG, Peterson LB, Wicker LS, Todd JA, Gough SC. Association of the T-cell regulatory gene CTLA4 with susceptibility to autoimmune disease. Nature. 2003;423:506–511. doi: 10.1038/nature01621. [DOI] [PubMed] [Google Scholar]
- 9.Korstanje R, Eriksson P, Samnegard A, Olsson PG, Forsman-Semb K, Sen S, Churchill GA, Rollins J, Harris S, Hamsten A, Paigen B. Locating Ath8, a locus for murine atherosclerosis susceptibility and testing several of its candidate genes in mice and humans. Atherosclerosis. 2004;177:443–450. doi: 10.1016/j.atherosclerosis.2004.08.006. [DOI] [PubMed] [Google Scholar]
- 10.Gharani N, Benayed R, Mancuso V, Brzustowicz LM, Millonig JH. Association of the homeobox transcription factor, ENGRAILED 2, 3, with autism spectrum disorder. Mol Psychiatry. 2004;9:474–484. doi: 10.1038/sj.mp.4001498. [DOI] [PubMed] [Google Scholar]
- 11.Wang X, Ria M, Kelmenson PM, Eriksson P, Higgins DC, Samnegard A, Petros C, Rollins J, Bennet AM, Wiman B, de Faire U, Wennberg C, Olsson PG, Ishii N, Sugamura K, Hamsten A, Forsman-Semb K, Lagercrantz J, Paigen B. Positional identification of TNFSF4, encoding OX40 ligand, as a gene that influences atherosclerosis susceptibility. Nat Genet. 2005;37:365–372. doi: 10.1038/ng1524. [DOI] [PubMed] [Google Scholar]
- 12.Hillebrandt S, Wasmuth HE, Weiskirchen R, Hellerbrand C, Keppeler H, Werth A, Schirin-Sokhan R, Wilkens G, Geier A, Lorenzen J, Kohl J, Gressner AM, Matern S, Lammert F. Complement factor 5 is a quantitative trait gene that modifies liver fibrogenesis in mice and humans. Nat Genet. 2005;37:835–843. doi: 10.1038/ng1599. [DOI] [PubMed] [Google Scholar]
- 13.Muller-Sieburg CE, Riblet R. Genetic control of the frequency of hematopoietic stem cells in mice: mapping of a candidate locus to chromosome 1. J Exp Med. 1996;183:1141–1150. doi: 10.1084/jem.183.3.1141. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 14.De Haan G, Van Zant G. Genetic analysis of hemopoietic cell cycling in mice suggests its involvement in organismal life span. Faseb J. 1999;13:707–713. doi: 10.1096/fasebj.13.6.707. [DOI] [PubMed] [Google Scholar]
- 15.Chen J, Astle CM, Harrison DE. Genetic regulation of primitive hematopoietic stem cell senescence. Exp Hematol. 2000;28:442–450. doi: 10.1016/s0301-472x(99)00157-5. [DOI] [PubMed] [Google Scholar]
- 16.Hasegawa M, Baldwin TM, Metcalf D, Foote SJ. Progenitor cell mobilization by granulocyte colony-stimulating factor controlled by loci on chromosomes 2 and 11. Blood. 2000;95:1872–1874. [PubMed] [Google Scholar]
- 17.de Haan G, Bystrykh LV, Weersing E, Dontje B, Geiger H, Ivanova N, Lemischka IR, Vellenga E, Van Zant G. A genetic and genomic analysis identifies a cluster of genes associated with hematopoietic cell turnover. Blood. 2002;100:2056–2062. doi: 10.1182/blood-2002-03-0808. [DOI] [PubMed] [Google Scholar]
- 18.Morrison SJ, Qian D, Jerabek L, Thiel BA, Park IK, Ford PS, Kiel MJ, Schork NJ, Weissman IL, Clarke MF. A genetic determinant that specifically regulates the frequency of hematopoietic stem cells. J Immunol. 2002;168:635–642. doi: 10.4049/jimmunol.168.2.635. [DOI] [PubMed] [Google Scholar]
- 19.Henckaerts E, Geiger H, Langer JC, Rebollo P, Van Zant G, Snoeck HW. Genetically determined variation in the number of phenotypically defined hematopoietic progenitor and stem cells and in their response to early-acting cytokines. Blood. 2002;99:3947–3954. doi: 10.1182/blood.v99.11.3947. [DOI] [PubMed] [Google Scholar]
- 20.Henckaerts E, Langer JC, Snoeck HW. Quantitative genetic variation in the hematopoietic stem cell and progenitor cell compartment and in lifespan are closely linked at multiple loci in BXD recombinant inbred mice. Blood. 2004;104:374–379. doi: 10.1182/blood-2003-12-4304. [DOI] [PubMed] [Google Scholar]
- 21.Geiger H, Szilvassy SJ, Ragland P, Van Zant G. Genetic analysis of progenitor cell mobilization by granulocyte colony-stimulating factor: verification and mechanisms for loci on murine chromosomes 2 and 11. Exp Hematol. 2004;32:60–67. doi: 10.1016/j.exphem.2003.09.015. [DOI] [PubMed] [Google Scholar]
- 22.Bystrykh L, Weersing E, Dontje B, Sutton S, Pletcher MT, Wiltshire T, Su AI, Vellenga E, Wang J, Manly KF, Lu L, Chesler EJ, Alberts R, Jansen RC, Williams RW, Cooke MP, de Haan G. Uncovering regulatory pathways that affect hematopoietic stem cell function using ’genetical genomics’. Nat Genet. 2005;37:225–232. doi: 10.1038/ng1497. [DOI] [PubMed] [Google Scholar]
- 23.van Os R, Ausema A, Noach EJ, van Pelt K, Dontje BJ, Vellenga E, de Haan G. Identification of quantitative trait loci regulating haematopoietic parameters in B6AKRF2 mice. Br J Haematol. 2006;132:80–90. doi: 10.1111/j.1365-2141.2005.05835.x. [DOI] [PubMed] [Google Scholar]
- 24.Liang Y, Jansen M, Aronow B, Geiger H, Van Zant G. The quantitative trait gene latexin influences the size of the hematopoietic stem cell population in mice. Nat Genet. 2007;39:178–188. doi: 10.1038/ng1938. [DOI] [PubMed] [Google Scholar]
- 25.Langer JC, Henckaerts E, Orenstein J, Snoeck HW. Quantitative trait analysis reveals transforming growth factor-beta2 as a positive regulator of early hematopoietic progenitor and stem cell function. J Exp Med. 2004;199:5–14. doi: 10.1084/jem.20030980. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 26.Henckaerts E, Langer JC, Orenstein J, Snoeck HW. The positive regulatory effect of TGF-beta2 on primitive murine hemopoietic stem and progenitor cells is dependent on age, genetic background, and serum factors. J Immunol. 2004;173:2486–2493. doi: 10.4049/jimmunol.173.4.2486. [DOI] [PubMed] [Google Scholar]
- 27.Massague J. TGF-beta signal transduction. Annu Rev Biochem. 1998;67:753–791. doi: 10.1146/annurev.biochem.67.1.753. [DOI] [PubMed] [Google Scholar]
- 28.Massague J, Chen YG. Controlling TGF-beta signaling. Genes Dev. 2000;14:627–644. [PubMed] [Google Scholar]
- 29.Massague J, Wotton D. Transcriptional control by the TGF-beta/Smad signaling system. Embo J. 2000;19:1745–1754. doi: 10.1093/emboj/19.8.1745. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 30.Larsson J, Karlsson S. The role of Smad signaling in hematopoiesis. Oncogene. 2005;24:5676–5692. doi: 10.1038/sj.onc.1208920. [DOI] [PubMed] [Google Scholar]
- 31.Ohta M, Greenberger JS, Anklesaria P, Bassols A, Massague J. Two forms of transforming growth factor-beta distinguished by multipotential haematopoietic progenitor cells. Nature. 1987;329:539–541. doi: 10.1038/329539a0. [DOI] [PubMed] [Google Scholar]
- 32.Goey H, Keller JR, Back T, Longo DL, Ruscetti FW, Wiltrout RH. Inhibition of early murine hemopoietic progenitor cell proliferation after in vivo locoregional administration of transforming growth factor-beta 1. J Immunol. 1989;143:877–880. [PubMed] [Google Scholar]
- 33.Sitnicka E, Ruscetti FW, Priestley GV, Wolf NS, Bartelmez SH. Transforming growth factor beta 1 directly and reversibly inhibits the initial cell divisions of long-term repopulating hematopoietic stem cells. Blood. 1996;88:82–88. [PubMed] [Google Scholar]
- 34.Jacobsen SE, Veiby OP, Myklebust J, Okkenhaug C, Lyman SD. Ability of flt3 ligand to stimulate the in vitro growth of primitive murine hematopoietic progenitors is potently and directly inhibited by transforming growth factor-beta and tumor necrosis factor-alpha. Blood. 1996;87:5016–5026. [PubMed] [Google Scholar]
- 35.Larsson J, Blank U, Helgadottir H, Bjornsson JM, Ehinger M, Goumans MJ, Fan X, Leveen P, Karlsson S. TGF-beta signaling-deficient hematopoietic stem cells have normal self-renewal and regenerative ability in vivo despite increased proliferative capacity in vitro. Blood. 2003;102:3129–3135. doi: 10.1182/blood-2003-04-1300. [DOI] [PubMed] [Google Scholar]
- 36.Larsson J, Blank U, Klintman J, Magnusson M, Karlsson S. Quiescence of hematopoietic stem cells and maintenance of the stem cell pool is not dependent on TGF-beta signaling in vivo. Exp Hematol. 2005;33:592–596. doi: 10.1016/j.exphem.2005.02.003. [DOI] [PubMed] [Google Scholar]
- 37.Spangrude GJ, Heimfeld S, Weissman IL. Purification and characterization of mouse hematopoietic stem cells. Science. 1988;241:58–62. doi: 10.1126/science.2898810. [DOI] [PubMed] [Google Scholar]
- 38.Sanford LP, Ormsby I, Gittenberger-de Groot AC, Sariola H, Friedman R, Boivin GP, Cardell EL, Doetschman T. TGFbeta2 knockout mice have multiple developmental defects that are non-overlapping with other TGFbeta knockout phenotypes. Development. 1997;124:2659–2670. doi: 10.1242/dev.124.13.2659. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 39.Markel P, Shu P, Ebeling C, Carlson GA, Nagle DL, Smutko JS, Moore KJ. Theoretical and empirical issues for marker-assisted breeding of congenic mouse strains. Nat Genet. 1997;17:280–284. doi: 10.1038/ng1197-280. [DOI] [PubMed] [Google Scholar]
- 40.Wakeland E, Morel L, Achey K, Yui M, Longmate J. Speed congenics: a classic technique in the fast lane (relatively speaking) Immunol Today. 1997;18:472–477. doi: 10.1016/s0167-5699(97)01126-2. [DOI] [PubMed] [Google Scholar]
- 41.Floersheim GL. Protection by bone marrow grafts against lethal doses of 5-fluorouracil. Eur J Cancer Clin Oncol. 1984;20:1441–1444. doi: 10.1016/0277-5379(84)90065-8. [DOI] [PubMed] [Google Scholar]
- 42.Mackay TF. The genetic architecture of quantitative traits: lessons from Drosophila. Curr Opin Genet Dev. 2004;14:253–257. doi: 10.1016/j.gde.2004.04.003. [DOI] [PubMed] [Google Scholar]
- 43.Flint J, Valdar W, Shifman S, Mott R. Strategies for mapping and cloning quantitative trait genes in rodents. Nat Rev Genet. 2005;6:271–286. doi: 10.1038/nrg1576. [DOI] [PubMed] [Google Scholar]
- 44.Kroymann J, Mitchell-Olds T. Epistasis and balanced polymorphism influencing complex trait variation. Nature. 2005;435:95–98. doi: 10.1038/nature03480. [DOI] [PubMed] [Google Scholar]
- 45.Luo K. Ski and SnoN: negative regulators of TGF-beta signaling. Curr Opin Genet Dev. 2004;14:65–70. doi: 10.1016/j.gde.2003.11.003. [DOI] [PubMed] [Google Scholar]
- 46.Warner DR, Horn KH, Mudd L, Webb CL, Greene RM, Pisano MM. PRDM16/MEL1: a novel Smad binding protein expressed in murine embryonic orofacial tissue. Biochim Biophys Acta. 2007;1773:814–820. doi: 10.1016/j.bbamcr.2007.03.016. [DOI] [PubMed] [Google Scholar]
- 47.Kwon BS, Hurtado JC, Lee ZH, Kwack KB, Seo SK, Choi BK, Koller BH, Wolisi G, Broxmeyer HE, Vinay DS. Immune responses in 4-1BB (CD137)-deficient mice. J Immunol. 2002;168:5483–5490. doi: 10.4049/jimmunol.168.11.5483. [DOI] [PubMed] [Google Scholar]
- 48.Wang Y, Subudhi SK, Anders RA, Lo J, Sun Y, Blink S, Wang J, Liu X, Mink K, Degrandi D, Pfeffer K, Fu YX. The role of herpesvirus entry mediator as a negative regulator of T cell-mediated responses. J Clin Invest. 2005;115:711–717. doi: 10.1172/JCI200522982. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 49.Drakesmith H, Townsend A. The structure and function of HFE. Bioessays. 2000;22:595–598. doi: 10.1002/1521-1878(200007)22:7<595::AID-BIES1>3.0.CO;2-E. [DOI] [PubMed] [Google Scholar]
- 50.Feder JN, Gnirke A, Thomas W, Tsuchihashi Z, Ruddy DA, Basava A, Dormishian F, Domingo R, Jr., Ellis MC, Fullan A, Hinton LM, Jones NL, Kimmel BE, Kronmal GS, Lauer P, Lee VK, Loeb DB, Mapa FA, McClelland E, Meyer NC, Mintier GA, Moeller N, Moore T, Morikang E, Prass CE, Quintana L, Starnes SM, Schatzman RC, Brunke KJ, Drayna DT, Risch NJ, Bacon BR, Wolff RK. A novel MHC class I-like gene is mutated in patients with hereditary haemochromatosis. Nat Genet. 1996;13:399–408. doi: 10.1038/ng0896-399. [DOI] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.