Abstract
Secondary malignancies are a serious adverse consequence of alkylator chemotherapy. The risk of developing an alkylator-associated malignancy is influenced by genetic background, although the relevant genetic factors are poorly understood. To screen for novel susceptibility factors, we established a mouse model of alkylator-induced malignancy. We exposed mice from 20 inbred strains to the prototypical alkylating agent, N-nitroso-N-ethylurea (ENU). ENU was a potent carcinogen in many of the strains tested, inducing 140 tumors in 240 ENU-treated mice (66% incidence of at least one tumor in evaluable mice), compared with a background incidence of 8% spontaneous tumors in 240 strain-, age-, and sex-matched control mice (relative risk, 8.4; P < 0.0001). A wide variety of tumor histologies were noted, including epithelial carcinomas, soft tissue sarcomas, and hematopoietic tumors. Cancer susceptibility was a heritable trait for the most common tumor types, lung adenocarcinoma (H2 = 0.25), T cell lymphoma (H2 = 0.19), and myeloid malignancies (H2 = 0.10). Quantitative trait locus mapping identified regions on chromosomes 3, 6, 9, and 15 containing candidate genes associated with lung adenoma, lung carcinoma, and lymphoma susceptibility. This novel mouse model recapitulates many features of human alkylator-associated cancer and supports the hypothesis that susceptibility to this syndrome is influenced by inherited polymorphisms that could be used to make informed clinical treatment decisions. (Cancer Res 2006; 66(10): 5029-38)
Introduction
Secondary cancers are a major problem complicating therapy for highly curable malignancies such as breast cancer, Hodgkin's disease, non–Hodgkin's lymphoma, and childhood acute lymphoblastic leukemia (1–3). Hematologic malignancies, particularly acute myeloid leukemia, are the most common class of secondary cancers, but solid tumors such as lung carcinoma also occur at an increased frequency following chemoradiotherapy (4). The prevalence of secondary leukemias is rising for reasons which may include improved patient survival after treatment of the primary tumor and more frequent and dose-intensive use of chemoradiotherapy. Alkylating agents cause a distinct clinical syndrome of treatment-related acute myeloid leukemia (t-AML). Alkylator-associated t-AML has a relatively long latency (5-7 years), frequently evolves from a myelodysplastic (MDS) phase, and is often associated with deletions on the long arm of chromosome 7 and/or chromosome 5 (5). Secondary leukemias have a dismal prognosis and are, therefore, a prime target for prevention and control strategies.
The incidence of secondary malignancies varies widely among studies from 1% to 25%, implying that susceptibility is likely multifactorial. Nongenetic risk factors include patient age, endocrine and metabolic status, primary tumor type, and chemotherapy dose. Inherited t-MDS/AML susceptibility factors include rare familial cancer predisposition syndromes associated with germline mutations in TP53, XPD, or NF1 (6–8). Mice heterozygous for a null Nf1 allele are predisposed to alkylator-induced myeloid malignancies (9), providing strong evidence that NF1 is a susceptibility gene for t-MDS/AML. These mutant alleles may increase the risk of t-MDS/AML in affected individuals, but accounts for disproportionately few cases. More prevalent polymorphisms in these and other genes (each with modest effect) may account for more of the cumulative risk in the general population. Retrospective case-control studies have yielded conflicting findings on the association between t-AML and polymorphisms in cytochrome p450 family members. Increased risk of t-AML was associated with the common allele in the 5′-promoter of CYP3A4 (odds ratio, 0.07 with variant allele) and lower magnitude risks (odds ratios, 2.7-2.8) were associated with poor metabolizer genotypes of CYP2C19 and CYP2D6 (10, 11). These findings were not reproduced in other studies (12, 13). An inactive allele of the detoxifying enzyme NAD(P)H:quinone oxidoreductase (NQ01) was found at higher frequencies in t-AML cases (12), but no association was found in other studies (13, 14). Studies examining variation in conjugating enzymes involved in phase II metabolism have shown mostly null associations. A recent meta-analysis of all published studies of glutathione S-transferase polymorphisms emphasized the modest and variable effect of GSTM1 (odds ratios, 1.05-1.37), GSTP1 (odds ratio, 1.63), and GSTT1 (odds ratios, 1.70-4.62) variants on risk of t-AML (15).
Mouse models are powerful tools to screen for novel genetic factors that underlie complex traits such as susceptibility to drug-induced cancer. We optimized protocols to induce, detect, and analyze tumors in mice using the prototypical alkylating agent, N-nitroso-N-ethylurea (ENU). ENU is a potent mutagen causing point mutations with a bias towards A:T transversions and A:G transitions. The most common reported ENU-induced tumors in mice occur in the lung, thymus, and liver. There is some evidence that genetic variability contributes to susceptibility in these tumors (16). However, the literature on ENU-induced cancer in mice is largely anecdotal, covering a limited number of strains and cancer sites. We did a systematic screen in diverse inbred strains of mice to better define the genetic contribution to alkylator-induced cancer susceptibility.
Materials and Methods
Mutagenesis protocol. To model alkylator-induced neoplasia in mice, we used the prototypical alkylator, ENU. ENU is the most potent mutagen known for mice. At typical doses (100-400 mg/kg), ENU induces point mutations roughly once every 1,000 loci (∼20-30 per mouse genome; ref. 17). The 200 mg/kg dose was chosen for this study because higher doses increase toxicity without significantly increasing the frequency of mutations per genome (18). ENU administration followed minor modifications of established protocols (19). In brief, 1 g Isopac vials (Sigma, St. Louis, MO) were reconstituted by gentle heating at a volume of 10 mL 10% ethanol. The solution was then brought to a final volume of 100 mL with phosphate citrate buffer [0.1 mol/L dibasic sodium phosphate, 0.05 mol/L sodium citrate (pH 5.0)]. Interlot variability was corrected by assaying each reconstituted vial spectrophotometrically immediately before use. We noted <20% lot-to-lot variabilities in ENU concentration. Mice were shipped into our specific pathogen–free barrier facility at 8 weeks of age from the research colony at The Jackson Laboratory (Bar Harbor, ME). Mice (n = 240 per cohort, equal number of each sex) received two doses of ENU (100 mg/kg, i.p.) or no treatment (n = 240 mice) at 9 and 10 weeks of age. Mice were observed daily for signs of illness (e.g., hunched posture, pallor, rapid breathing). All studies were conducted with the approval of the Animal Studies Committee at Washington University School of Medicine. Mice were sacrificed and analyzed when moribund or if they reached 16 months postinjection.
Hematologic analysis. Complete blood counts with automated five-part differentials were determined with a Forcyte veterinary cell counter (Oxford Sciences, Inc., Oxford, CT). For flow cytometric analysis, bone marrow cells were stained with directly conjugated antibodies to Sca1, CD117, CD34, CD11b, Gr1, B220, CD3, CD4, or CD8 (BD PharMingen, San Diego, CA). Scatter-gated events (25,000) were acquired on a FACScan (BD Biosciences) and analyzed with FlowJo software (ver 4.5.1, Treestar, Inc., Ashland, OR).
Histology. Tissues were fixed in 10% formalin and embedded in paraffin. Five-micron sections were stained with H&E. Peripheral blood, bone marrow smears, and spleen touch preps were stained with May-Grünwald/Giemsa. All slides were reviewed in a blinded fashion by a veterinary pathologist and a clinical hematopathologist.
Statistical analysis. Cumulative probability of survival was estimated by the Kaplan-Meier method. Broad-sense heritability (H2) was estimated using an ANOVA model with strain as a random factor and cancer phenotype (absent, present) as the dependent variable (20). The variance in cancer incidence between strains is the total genetic variance, whereas the variance within strains is environmental in origin (21). Broad-sense heritability is the ratio between strain variance to the total phenotypic variance (between-strain variance + within-strain variance). Genetic correlations among cancers were calculated using MANOVA, with strain as the factor and the five cancer phenotypes (lung adenoma, lung adenocarcinoma, lymphoma, myeloid leukemia, and gastrointestinal carcinoma) as the dependent variables. Environmental correlations among cancers are given by the residual correlation matrix from the MANOVA.
Quantitative trait locus mapping. The genotypes of 13,287 single nucleotide polymorphisms (SNP) were obtained for the 20 strains in this study from the Wellcome-CTC Mouse Strain SNP Genotype Set (http://www.well.ox.ac.uk/mouse/INBREDS). This set of markers places an informative SNP approximately every 160 kb across the genome. In a modification of published in silico mapping techniques (22), a nested ANOVA model was used to test for differences in cancer frequency between genotypes at each SNP. The use of dichotomous-dependent variables (cancer absent/present) is in formal violation of the ANOVA model because residuals are not normally distributed, but in practice, dichotomous-dependent variables can be usefully employed in linear models (23). We controlled for multiple comparisons using the false discovery rate approach (24). When a plot of −log10(rank) versus −log10(P) deviated markedly from linearity, the null hypothesis that no quantitative trait locus (QTL) exists was rejected. QTLs were eliminated from further consideration if there was 100% concordance between strain distributions at multiple loci because of the inability to discriminate true QTLs from “shadow” QTLs in these instances (25). Candidate regions were defined by the QTL SNP and adjacent SNPs that retain significant association with the cancer phenotype (P < 0.05). All annotated genes within these intervals were identified using build 35 of the mouse genome.
Results
Hematologic variables vary among inbred mouse strains. We selected 20 inbred strains designated highest priority by the Mouse Phenome Project steering committee (The Jackson Laboratory). These strains were chosen to represent the diversity of inbred mouse strains and include the most commonly used strains. Although some strains in this set are relatively closely related to one another (e.g., SJL/J, SWR/J, and FVB/NJ; C57L/J, C58/J, and C57BL/6J), analyses of the relationship between similarities in cancer experience and genotype indicates no significant similarity in cancer susceptibility among related strains. Therefore, these strains can be considered as independent with respect to cancer susceptibilities. We first asked whether hematologic variables differed among these strains at baseline. Age- and sex-matched mice (n = 24 per strain) were obtained from The Jackson Laboratory research colony. Baseline hematology studies were done at 9 weeks of age. All variables have significant (P < 0.0001) heritable components (Supplementary Table S1). Hematologic variables are not genetically correlated with each other. There are significant sex differences for each of the hematologic variables with males having higher values on all variables, except lymphocyte count for which females have the higher level. There are also significant sex by strain interactions, indicating genetic variation in levels of sexual dimorphism for WBC count, hemoglobin, hematocrit, and platelet count. For these traits, some strains are sexually dimorphic, whereas others are monomorphic. Heritability values for the hematologic variables were also obtained separately by sex. All traits are significantly heritable in both sexes (P < 0.0001), with heritability values ranging from a low of 15% for hematocrit in males to a high of 56% for lymphocyte count in males.
We compared the WBC values obtained in this study to previously reported values for the same strains (26). Overall, there is no significant difference in WBC values between the studies. However, the rank orders of the strain-specific values did vary from study to study. The correlation between strain means, or repeatability, across the two studies is 0.61 in females and 0.56 in males. These are only moderate repeatabilities and show that the different environments and procedures followed in the two studies interact with the strains' genetic backgrounds to modify these phenotypes.
Mice are highly susceptible to alkylator-induced cancer. ENU injections were well-tolerated in all strains. Necropsy data were adequate to characterize tumors in 172 ENU-treated and 191 control mice. The most common tumors were hematopoietic, with the lymphoid lineage affected more often than the myeloid lineage (Table 1). Benign and malignant lung tumors were frequently observed. A variety of other tumors of epithelial origin were noted in the gastrointestinal tract (e.g., stomach, large intestine, rectum) and genitourinary system (e.g., uterus, ovary, kidney). Less common ENU-induced tumors included fibrosarcoma (two cases), schwannoma (one), rhabdomyosarcoma (one), pheochromocytoma (one), breast carcinoma (one), and hemangiosarcoma (one). Overall, ENU-treated mice had a significantly increased risk of developing at least one tumor (relative risk, 8.44; 95% confidence interval, 5.13-13.88; P < 0.0001 by Fisher's exact test). Multiple independent tumors were identified in some mice (n = 30 with two tumors, n = 6 with more than two tumors).
Spectrum of tumors induced by ENU in a survey of 20 inbred strains
Tissue . | Histology . | ENU . | Control . |
---|---|---|---|
Hematopoietic | Lymphoma | 42 | 6 |
Myeloid leukemia | 9 | 1 | |
Myelodysplasia | 4 | 1 | |
Histiocytic sarcoma | 1 | 3 | |
Lung | Adenoma | 13 | 2 |
Adenocarcinoma | 30 | 1 | |
Gastrointestinal | Adenoma, carcinoma | 9 | 0 |
Genitourinary | Adenoma, carcinoma | 12 | 0 |
Breast | Carcinoma | 1 | 1 |
Liver | Adenoma | 7 | 0 |
Hemangioma | 3 | 0 | |
Carcinoma | 3 | 0 | |
Skin | Fibroma | 1 | 0 |
Other | Sarcoma | 10 | 3 |
Total tumors | 145 | 18 | |
Evaluable mice | 172 | 191 | |
Mice with ≥1 tumor | 114 | 15 | |
Cancer incidence | 0.66 | 0.08 |
Tissue . | Histology . | ENU . | Control . |
---|---|---|---|
Hematopoietic | Lymphoma | 42 | 6 |
Myeloid leukemia | 9 | 1 | |
Myelodysplasia | 4 | 1 | |
Histiocytic sarcoma | 1 | 3 | |
Lung | Adenoma | 13 | 2 |
Adenocarcinoma | 30 | 1 | |
Gastrointestinal | Adenoma, carcinoma | 9 | 0 |
Genitourinary | Adenoma, carcinoma | 12 | 0 |
Breast | Carcinoma | 1 | 1 |
Liver | Adenoma | 7 | 0 |
Hemangioma | 3 | 0 | |
Carcinoma | 3 | 0 | |
Skin | Fibroma | 1 | 0 |
Other | Sarcoma | 10 | 3 |
Total tumors | 145 | 18 | |
Evaluable mice | 172 | 191 | |
Mice with ≥1 tumor | 114 | 15 | |
Cancer incidence | 0.66 | 0.08 |
NOTE: Twenty-four mice from each of 20 strains were divided equally into ENU and control groups. The total number of tumors identified in each cohort (and mice bearing at least one tumor / evaluable mice) is shown.
Susceptibility to lung and hematopoietic tumors varies by strain. A wide variation in susceptibility between strains was observed for the three most commonly encountered ENU-induced malignancies (lung carcinoma, lymphoma, and myeloid malignancies). Lung carcinoma occurred in 50% or more of the A/J, PERA/EiJ, and SWR/J mice treated with ENU (Fig. 1A). Lymphoma developed at a rate of 50% or more in AKR/J, C58/J, C57BL/6J, and NOD/LtJ mice (Fig. 1B). Myeloid malignancies (defined here as acute myeloid leukemia, granulocytic sarcoma, or myelodysplasia) were observed in four out of five (80%) evaluable SWR/J mice, three out of six (50%) DBA/2J mice, and three out of nine (33%) BALB/cByJ mice (Fig. 1C). For each of the top three cancer phenotypes, we identified at least five completely resistant strains.
Incidence of three major cancer phenotypes varies by strain. The incidence of ENU-induced (A) lung carcinoma, (B) lymphoma, and (C) myeloid malignancy (acute myeloid leukemia, granulocytic sarcoma, or myelodysplasia) is shown for 20 strains of mice. The denominator is the number of evaluable mice in each strain (n = 6-12).
Incidence of three major cancer phenotypes varies by strain. The incidence of ENU-induced (A) lung carcinoma, (B) lymphoma, and (C) myeloid malignancy (acute myeloid leukemia, granulocytic sarcoma, or myelodysplasia) is shown for 20 strains of mice. The denominator is the number of evaluable mice in each strain (n = 6-12).
ENU-treated animals from many of the strains included in this screen were incidentally found to have small lung adenomas at the time of sacrifice. They have a nonpapillary architecture and are composed predominantly of Clara cells with occasional scattered type II pneumocytes (Fig. 2A and B). In contrast, the ENU-induced lung adenocarcinomas were highly aggressive, symptomatic tumors. Affected mice typically developed large bilateral pulmonary masses which caused rapid, shallow breathing and lethargy. Histologic examination of these tumors revealed typical features of adenocarcinoma with acinar architecture, nuclear atypia, and a propensity to invade neighboring structures (Fig. 2C and D). These tumors developed with a relative long latency of ∼40 weeks after ENU exposure (50 weeks of age; Fig. 3A).
Histology of ENU-induced mouse tumors. A and B, pulmonary alveolar/bronchiolar adenoma. C and D, pulmonary alveolar/bronchiolar carcinoma, acinar type. E-H, T cell lymphoma. I-L, acute myeloid leukemia. A, hyperplasia of alveolar cells exhibiting a lepidic growth pattern (original magnification, ×100). B, densely packed cuboidal to columnar cells with bright eosinophilic cytoplasm, predominantly Clara cells with rare type II pneumocytes interspersed (original magnification, ×400). C, discrete tumor mass invading an adjacent bronchiole (original magnification, ×40). D, well-differentiated adenocarcinoma composed of atypical polygonal cells which have retained the ability to form glandular lumina (original magnification, ×200). E, peripheral blood smear shows marked leukocytosis involving predominantly abnormal lymphoid cells. Two malignant cells are seen undergoing mitosis (original magnification, ×500). F, bone marrow is heavily infiltrated by lymphoid cells. Megakaryocytes are present, but rare mature myeloid forms are seen (original magnification, ×500). G, spleen shows lymphoid cells infiltrating both the red and white pulp (original magnification, ×200). H, malignant lymphoid infiltrate concentrated in periportal regions of the liver (original magnification, ×400). I, peripheral blood smear shows blasts with fine chromatin, irregular nuclear contours, and a modest amount of cytoplasm containing occasional azurophilic granules. Note the presence of large platelets and the absence of mature myeloid cells (original magnification, ×500). J, bone marrow. There is a predominance of myeloblasts with rare mature ring forms (original magnification, ×500). K, spleen. Red and white pulp are diffusely infiltrated by immature, monomorphic cells. Mitotic figures are abundant. Hematopoiesis is markedly reduced (original magnification, ×400). L, liver. Monomorphic blasts infiltrate both the periportal space and the sinuses. In addition, there is destruction of liver parenchyma and extramedullary hematopoiesis (original magnification, ×200).
Histology of ENU-induced mouse tumors. A and B, pulmonary alveolar/bronchiolar adenoma. C and D, pulmonary alveolar/bronchiolar carcinoma, acinar type. E-H, T cell lymphoma. I-L, acute myeloid leukemia. A, hyperplasia of alveolar cells exhibiting a lepidic growth pattern (original magnification, ×100). B, densely packed cuboidal to columnar cells with bright eosinophilic cytoplasm, predominantly Clara cells with rare type II pneumocytes interspersed (original magnification, ×400). C, discrete tumor mass invading an adjacent bronchiole (original magnification, ×40). D, well-differentiated adenocarcinoma composed of atypical polygonal cells which have retained the ability to form glandular lumina (original magnification, ×200). E, peripheral blood smear shows marked leukocytosis involving predominantly abnormal lymphoid cells. Two malignant cells are seen undergoing mitosis (original magnification, ×500). F, bone marrow is heavily infiltrated by lymphoid cells. Megakaryocytes are present, but rare mature myeloid forms are seen (original magnification, ×500). G, spleen shows lymphoid cells infiltrating both the red and white pulp (original magnification, ×200). H, malignant lymphoid infiltrate concentrated in periportal regions of the liver (original magnification, ×400). I, peripheral blood smear shows blasts with fine chromatin, irregular nuclear contours, and a modest amount of cytoplasm containing occasional azurophilic granules. Note the presence of large platelets and the absence of mature myeloid cells (original magnification, ×500). J, bone marrow. There is a predominance of myeloblasts with rare mature ring forms (original magnification, ×500). K, spleen. Red and white pulp are diffusely infiltrated by immature, monomorphic cells. Mitotic figures are abundant. Hematopoiesis is markedly reduced (original magnification, ×400). L, liver. Monomorphic blasts infiltrate both the periportal space and the sinuses. In addition, there is destruction of liver parenchyma and extramedullary hematopoiesis (original magnification, ×200).
Survival analysis. A, lung carcinoma is highly penetrant with a latency of ∼40 weeks after ENU exposure. B, lymphoma has a shorter latency (∼20 weeks) and intermediate penetrance. C, myeloid malignancies have a longer latency (∼30 weeks) and lower penetrance compared with lymphoma. Two independent cohorts of ENU-treated AKR/J mice (n = 12 and n = 40) had similar outcomes. The cumulative probability of survival was computed by the Kaplan-Meier method for each of the three end points with all other outcomes censored. Strains selected for this analysis had the highest incidence of the cancer phenotype and fewer than half of the observations censored.
Survival analysis. A, lung carcinoma is highly penetrant with a latency of ∼40 weeks after ENU exposure. B, lymphoma has a shorter latency (∼20 weeks) and intermediate penetrance. C, myeloid malignancies have a longer latency (∼30 weeks) and lower penetrance compared with lymphoma. Two independent cohorts of ENU-treated AKR/J mice (n = 12 and n = 40) had similar outcomes. The cumulative probability of survival was computed by the Kaplan-Meier method for each of the three end points with all other outcomes censored. Strains selected for this analysis had the highest incidence of the cancer phenotype and fewer than half of the observations censored.
The lymphoid neoplasms observed were also highly aggressive. In most cases, there was symptomatic mediastinal involvement. Enlargement of mesenteric, cervical, or axillary lymph nodes was less commonly noted. Numerous abnormal lymphocytes, including frequent mitotic cells, could be found in peripheral blood smears from these mice (Fig. 2E). Bone marrow and spleen were often extensively involved, so that the megakaryocytes, maturing myeloid elements and erythroid precursors normally present in these hematopoietic tissues were diminished and partly replaced by a lymphoid infiltrate (Fig. 2F and G). In severely affected animals, nonlymphoid organs such as the liver and kidneys were diffusely infiltrated by tumor cells (Fig. 2H). Flow cytometric analysis showed an immature CD4+CD8+ immunophenotype in most cases (n = 22), but single positive CD4+CD8− (n = 4) or CD8+CD4− (n = 11) T cell lymphomas and B cell lymphomas (n = 4) were also seen (Supplementary Fig. S1). The lymphomas had a relatively short latency of ∼20 weeks (Fig. 3B).
The myeloid malignancies were variable in presentation. Many of the animals with acute myeloid leukemia seemed lethargic or tachypneic. Laboratory studies revealed one or more of the following: leukocytosis (20-130,000 cells/μL), anemia (hemoglobin, 1.2-13.1 g/dL), or thrombocytopenia. Peripheral blood and/or bone marrow evaluation revealed excess myeloblasts (Fig. 2I and J); however, there were cases in which the disease was primarily extramedullary (lymph node, mediastinum, or granulocytic sarcoma). Affected mice had splenomegaly, often massive (>20-fold enlargement) with infiltration of the red pulp by myeloblasts that encroach on or replace the megarkaryocytes and erythropoietic elements normally present in the mouse spleen (Fig. 2K). In some cases, the portal triads and the sinuses of the liver were infiltrated by blasts and abnormal extramedullary hematopoiesis could be seen within the liver parenchyma (Fig. 2L). By flow cytometry, cells were positive for CD117 or CD34 and lacked expression of lymphoid markers (Supplementary Fig. S1). A minority of cases were biphenotypic (usually CD117+, CD34+, and CD3lo or CD4lo). Myelodysplasia was usually identified incidentally postmortem in animals which appeared to be grossly normal. These animals had mild to moderate anemia with splenomegaly and no evidence of AML or lymphoma. Histopathologic examination of the spleen revealed extensive extramedullary hematopoiesis with numerous micromegakaryocytes and dysplasia in the myeloid and erythroid lineages. Flow cytometric analysis in these MDS cases typically showed a modest increase in Ter119+ erythroid precursors in the spleen, but was otherwise unremarkable. The myeloid malignancies had an intermediate latency of ∼30 weeks (Fig. 3C).
Given the relatively high incidence of both lymphoid and myeloid malignancies seen in AKR/J mice, a second independent cohort (n = 40 ENU-treated and 6 control mice) was studied. Among 37 evaluable ENU-treated mice, a similar frequency of lymphomas and myeloid malignancies were again observed, with consistent clinical, histologic, and flow cytometric features (Fig. 3C; data not shown).
Susceptibility to lung and hematopoietic tumors is a heritable trait. The interstrain variability in cancer susceptibility was studied to determine which of these phenotypes had a significant heritable component. Most cancers were heritable (Table 2). Lung carcinoma has the highest heritability (H2 = 0.25, P < 10−9) with strains A/J (nine), PERA/EiJ (six), FVB/NJ (four), and SWR/J (three) accounting for most of the cancer cases. Lymphoma had a heritability of 19% (P < 10−6) with AKR/J (seven), C58/J (six), FVB/NJ (six), NOD/LtJ (five), and C57BL/6J (five) accounting for most of the cancer cases. Myeloid malignancies have a heritability of only 10% (P < 0.001) with SWR/J (four), DBA/2J (three), AKR/J (two), NOD/LtJ (two) having multiple cases. Gastrointestinal cancers have negligible heritability (H2 = 0.06, P < 0.05) with strains CAST/EiJ (three), A/J (two), C57L/J (two), and SM/J (two) having more than one case each. Finally, lung adenomas are not significantly heritable (H2 = 0.03) with strains 129S1/SvlmJ (three), C3H/HeJ (2) and SM/J (two) having more than one case each.
Heritability and correlation of cancer phenotypes
![]() |
![]() |
NOTE: Broad-sense Heritability (diagonal, shaded), genetic correlations (below diagonal), and environmental (above diagonal) correlations between cancer types for a set of 20 inbred mouse strains. Values statistically significantly different from zero are in boldface.
The genetic and environmental correlations between cancers are given in Table 2. Interestingly, strains with relatively high rates of lung adenomas or gastrointestinal cancers tended not to have lymphoma or myeloid malignancies. However, there is no correlation between lung adenoma and gastrointestinal cancer rates. This negative genetic association between strains indicates either that the same genes affect both kinds of cancer producing susceptibility to lung adenoma or gastrointestinal cancers and resistance to lymphoma and myeloid malignancies or that closely linked genes have these effects separately. This lack of correlation can also be accounted for by the short latency and high mortality of the hematopoietic tumors obscuring susceptibility to the longer latency lung and gastrointestinal tumors. There are no significant environmental correlations between cancers (Table 2).
QTL mapping. Whole genome scans were done to identify regions associated with susceptibility to the most prevalent alkylator-induced tumors observed in these 20 strains. For each of the major tumor types, loci were found that significantly diverge from a null hypothesis of no association (Fig. 4). A false discovery threshold of 10% was used. For lung adenoma susceptibility, four loci were identified on chromosomes 3, 6, and 9. Thirty-four annotated genes are contained within these four QTLs (Table 3), including several encoding cell surface proteins (Cd276, Cd1d1, and Sdfr1) and transcription factors (Pml and Pax4). Lung adenocarcinoma susceptibility was associated with 10 potential genomic intervals. At two of these (both on chromosome 15), the clustering of sensitive and resistant strains was unique. One of these regions contains the prolactin receptor gene, Prlr (Table 3). The other loci had identical strain distributions, representing either multiple QTLs each with modest effects or the shadow of a single QTL cast on multiple sites with nonsyntenic associations (25). Because we cannot differentiate between these two possibilities, these potential QTLs were eliminated from further consideration. Lymphoma susceptibility mapped to two regions on chromosomes 3 and 15, one of which contains the gene encoding Sprouty1, an antagonist of fibroblast growth factor signaling (Table 3).
Genomic loci associated with susceptibility to alkylator-induced cancers. The log probability ratios (LPR) are shown around the genomic intervals associated with susceptibility to (A) lung carcinoma, (B) lung adenoma, and (C) lymphoma. Horizontal line, a false discovery threshold of 10%. Two QTL on chromosome 3 for lung adenoma susceptibility are depicted on a single plot.
Genomic loci associated with susceptibility to alkylator-induced cancers. The log probability ratios (LPR) are shown around the genomic intervals associated with susceptibility to (A) lung carcinoma, (B) lung adenoma, and (C) lymphoma. Horizontal line, a false discovery threshold of 10%. Two QTL on chromosome 3 for lung adenoma susceptibility are depicted on a single plot.
QTL associated with alkylator-induced cancer susceptibility
Tumor . | SNP . | Log probability ratio . | Interval . | Candidate genes . |
---|---|---|---|---|
Lung adenoma | rs13478675 | 4.44 | 6_27.423.114-6_28.851.991 | Grm8 |
Pax4 | ||||
Fscn3 | ||||
Arf5 | ||||
Gcc1 | ||||
rs13480220 | 3.86 | 9_58.102.941-9_58.995.534 | Stoml1 | |
Pml | ||||
Islr | ||||
Stra6 | ||||
Tbc1d21 | ||||
Cd276 | ||||
Sdfr1 | ||||
Hcn4 | ||||
Loxl1 | ||||
Neo1 | ||||
rs13477230 | 3.62 | 3_81.359.882-3_82.799.868 | Ctso | |
Tdo2 | ||||
Accn5 | ||||
Gucy1b3 | ||||
Gucy1a3 | ||||
Lrat | ||||
Fgg | ||||
Fga | ||||
Fgb | ||||
Plrg1 | ||||
Npy2r | ||||
rs13477242 | 3.62 | 3_85.631.744-3_86.825.116 | Sh3d19 | |
Rps3a | ||||
Mab21l2 | ||||
Cd1d1 | ||||
Dcamkl2 | ||||
Hig1 | ||||
Lrba | ||||
Cd1d2 | ||||
Lung carcinoma | mCV25349597 | 5.28 | 15_9.948.126-15_10.288.804 | Prlr |
mCV24802203 | 4.15 | 15_37.305.956-15_37.631.176 | Ncald | |
Tcfcp2l3 | ||||
Lymphoma | rs3671702 | 3.77 | 15_22.118.644-15_23.644.895 | Cdh18 |
rs3705939 | 3.68 | 3_36.903.931-3_37.866.074 | Spry1 | |
Spata5 |
Tumor . | SNP . | Log probability ratio . | Interval . | Candidate genes . |
---|---|---|---|---|
Lung adenoma | rs13478675 | 4.44 | 6_27.423.114-6_28.851.991 | Grm8 |
Pax4 | ||||
Fscn3 | ||||
Arf5 | ||||
Gcc1 | ||||
rs13480220 | 3.86 | 9_58.102.941-9_58.995.534 | Stoml1 | |
Pml | ||||
Islr | ||||
Stra6 | ||||
Tbc1d21 | ||||
Cd276 | ||||
Sdfr1 | ||||
Hcn4 | ||||
Loxl1 | ||||
Neo1 | ||||
rs13477230 | 3.62 | 3_81.359.882-3_82.799.868 | Ctso | |
Tdo2 | ||||
Accn5 | ||||
Gucy1b3 | ||||
Gucy1a3 | ||||
Lrat | ||||
Fgg | ||||
Fga | ||||
Fgb | ||||
Plrg1 | ||||
Npy2r | ||||
rs13477242 | 3.62 | 3_85.631.744-3_86.825.116 | Sh3d19 | |
Rps3a | ||||
Mab21l2 | ||||
Cd1d1 | ||||
Dcamkl2 | ||||
Hig1 | ||||
Lrba | ||||
Cd1d2 | ||||
Lung carcinoma | mCV25349597 | 5.28 | 15_9.948.126-15_10.288.804 | Prlr |
mCV24802203 | 4.15 | 15_37.305.956-15_37.631.176 | Ncald | |
Tcfcp2l3 | ||||
Lymphoma | rs3671702 | 3.77 | 15_22.118.644-15_23.644.895 | Cdh18 |
rs3705939 | 3.68 | 3_36.903.931-3_37.866.074 | Spry1 | |
Spata5 |
NOTE: SNPs defining each QTL are shown. The log probability ratio at each SNP is equal to log10(1/P). The indicated interval is the genomic region flanking the QTL SNP that retains significant (P < 0.05) association with each cancer phenotype. Annotated genes in each interval are shown. ESTs and predicated genes were omitted.
Discussion
We performed a screen in mice to define the contribution of genetic variability to alkylator-induced cancer susceptibility. We mutagenized mice with ENU and did extensive phenotypic analysis on the resulting tumors. The drug, dose, and schedule of administration were well-tolerated in all 20 strains examined. ENU was a potent carcinogen, inducing a large number of tumors with significant interstrain variability. The most common cancer phenotypes were lung adenocarcinoma and hematopoietic malignancies of both the lymphoid and myeloid lineages. The heritability of these traits is similar to the heritabilities of common human phenotypes studied in partially inbred populations, such as plasma insulin levels (H2 = 0.32), blood pressure (H2 = 0.21-0.45), and measures of lung function (H2 = 0.35-0.41; ref. 27). A large study of cancer incidence in twins found that ∼40% of cancer risk was heritable (28), emphasizing the importance of genetics, environment, and gene/environment interactions in determining cancer risk. Heritability estimates in the current study are, admittedly, for very artificial populations. Nevertheless, they are an accurate measure of the variance in cancer susceptibility attributable to genetic background. The strains chosen are a standard set in mouse genetics. A range of inbred strains and Mus subspecies were included to increase genetic diversity. The high penetrance of specific cancer phenotypes we observed shows the value of this approach for the identification of genetic susceptibility factors for alkylator-induced cancer.
Resting blood counts varied little within strains and were concordant with previously published results (26). There was significant interstrain variability in total leukocyte count, neutrophil count, lymphocyte count, hemoglobin, hematocrit, and platelet count, demonstrating that the size of these hematopoietic compartments is under genetic control. Previous studies have mapped QTLs that affect peripheral blood leukocyte counts, lymphocyte proportions, and hematopoietic stem cell number and function (29, 30). In the present study, baseline hematologic variables were not correlated with cancer incidence after ENU treatment. These results suggest that distinct genetic pathways affect basal hematopoiesis and susceptibility to alkylator-induced hematopoietic malignancies.
We found a high rate of both premalignant and malignant lung tumors in several strains of mice following exposure to ENU. In some strains (e.g., A/J, PERA/EiJ, SWR/J), the incidence of lung cancer exceeded 50%. Several other strains (e.g., the C57-related strains C57L/J, C58/J, and C57BL/6J) were completely resistant. It has been known for many years that certain inbred mouse strains have increased susceptibility to spontaneous lung tumors, and that strains with a high (e.g., A/J, SWR/J) or intermediate (BALB/c) spontaneous rate of lung tumors are also more susceptible to chemical induction of tumors using cigarette smoke, tar, or purified chemical carcinogens, including ENU (16, 31). Our results update these earlier studies using modern strains and extend them to include several commonly used laboratory and wild-derived strains (SPRET/EiJ, MOLF/EiJ, and CAST/EiJ) not previously evaluated for this phenotype.
A large number of pulmonary adenoma susceptibility (Sluc, Pas; refs. 32, 33) and resistance (Par; ref. 34) loci have been mapped, and several candidate genes have been cloned, notably Kras2 (the murine homologue of human KRAS; ref. 35). Kras2 is frequently mutated in ENU-induced pulmonary adenomas (36) and the KrasG12D allele efficiently induces pulmonary adenomas when conditionally expressed in the lung (37). We did association mapping to identify potential cancer susceptibility loci. Several interesting candidate regions were identified, although these results need to be validated and complemented using traditional intercross mapping strategies because of the known statistical power limitations of this technique (38). We found four loci on chromosomes 6, 9, and 3 associated with lung adenoma susceptibility. One of these (on chromosome 6 at ∼28 Mb) is near a susceptibility locus (Sluc7) previously mapped using recombinant congenic strains (30). The remaining three QTLs are not closely linked to previously described loci.
There is increasing awareness that alkylating agents increase the risk of lung adenocarcinoma in humans. This is particularly evident in patients with Hodgkin's disease where the relative risk ranges from 4.2 to 6.8, and is a function of increasing alkylator dose (4). The genetic basis of human susceptibility to lung cancer is not well-understood, but at least one candidate susceptibility locus has been mapped to chromosome 6q23 (39), a region not syntenic to the QTLs we identified in mice. Our analysis has identified two potential lung carcinoma susceptibility loci on chromosome 15. Three susceptibility loci (Sluc 24, 25, and 26) have been mapped to this chromosome previously (32). Sluc24 is tightly linked to the interval we identified containing the prolactin receptor gene. Sluc25 is close to the distal QTL we identified on chromosome 15 at ∼38 Mb.
We found a high frequency of aggressive ENU-induced T cell lymphomas in several strains, with AKR/J and C58/J showing the highest incidence. Most strains showed intermediate susceptibility, and several strains were completely resistant (A/J, SWR/J, and the wild-derived strains). These findings are consistent with data from The Jackson Laboratory in which the frequency of spontaneous lymphomas is reported to be highest in AKR/J and C58/J, and was negligible in A/J and SWR/J (40). Lymphomas arise spontaneously in AKR/J, C58/J, and other strains by somatic recombination of latent murine leukemia viruses that activate oncogenes or inactivate tumor suppressors via insertional mutagenesis (41). Most of the recombinant inbred lines between AKR/J and DBA/2J (AKXD) are also highly susceptible to B and T cell lymphomas. The lymphomas we and others have observed in treated AKR/J and C58/J mice, therefore, represent the cumulative effect of somatic mutations induced by ENU and retroviral integration. Treatment-related lymphoid malignancies are rare in humans, but do occur (42). Therefore, mapping the susceptibility alleles for ENU-induced lymphomas may have potential for improving human health.
Several strains in our survey were susceptible to a spectrum of ENU-induced myeloid malignancies that ranged from myelodysplasia (e.g., SWR/J), to granulocytic sarcoma (e.g., DBA/2J), and acute myeloid leukemia (e.g., AKR/J, BALB/cByJ, and NOD/LtJ). These phenotypes have not been reported in prior ENU screens. None of the commonly used modern inbred strains has a significant spontaneous rate of myeloid leukemias. In the course of generating the inbred AKR/J strain, Jacob Furth described a subline highly susceptible to myeloid leukemia which, unfortunately, was subsequently lost (43). Although most of the AKXD RI lines develop lymphomas exclusively, AKXD-23 mice develop AML with nearly full penetrance by 1 year (44). BXH-2 recombinant inbred mice also develop AML at high penetrance, whereas neither of the parental strains (C57BL/6J and C3H/HeJ) or any of the other BXH RI lines has a significant spontaneous incidence of AML (45). BXH-2 mice acquire infectious retrovirus by maternal-fetal transmission and are viremic throughout life. The identification of common retroviral integration sites in these mice has proved to be a very successful strategy for the identification of novel somatic events that contribute to AML pathogenesis (46), but again does not shed light on germline polymorphisms that affect susceptibility to drug-induced cancer. Numerous transgenic strains engineered to express fusion genes associated with human AML (e.g., AML1/ETO and PML/RARA) fail to develop AML at high penetrance, implying that AML is not a “single hit” disease. ENU mutagenesis in several of these transgenic strains, but not in wild-type littermates, generated cooperating mutations that were sufficient to induce AML (47, 48). We observed AML in some strains after ENU, implying that inherited polymorphisms in these strains provide the first genetic “hit.” This phenotype is the closest correlate to the most prevalent type of alkylator-induced cancer occurring in humans. Therefore, the strains we have identified are of particular interest for mapping susceptibility alleles for t-MDS/AML.
Although we present evidence here that susceptibility of mice to alkylator-induced cancer has a statistically significant and biologically relevant genetic component, a substantial amount of interstrain variability in this phenotype has a nonheritable cause. This is concordant with human twin studies that show up to 60% of cancer risk has an environmental or shared genetic/environmental cause (28). These variables may include nutrition, stress, or epigenetic modification of gene expression.
In addition to the well-established link between alkylating agent exposure and t-MDS/AML in humans, it is now apparent that these drugs also increase the risk of secondary solid tumors such as lung adenocarcinoma (4), bladder carcinoma (49), and sarcoma (50). These observations provide additional validation of our ENU screen in mice and justify further study of these strains to understand the genetic basis of susceptibility to a broad range of alkylator-induced malignancies. The long-term goal of these studies is to create screening tools that will allow for prospective risk stratification of patients so that therapy can be tailored to lower the risk of secondary malignancies.
Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org/).
Acknowledgments
Grant support: NIH PO1 CA101937 (PI Project 6, T. Graubert). Support from AstraZeneca was used to defray Mouse Phenome Project production costs.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank Piia Hanson and Masayo Izumi for expert technical assistance, Dr. Marie LaRegina (Department of Comparative Medicine, Washington University School of Medicine) for analysis and interpretation of solid tumor pathology, and S.C. Grubb (The Jackson Laboratory) for assistance with Mouse Phenome Database data submission. Flow cytometric analysis was done in the Siteman Cancer Center High Speed Cell Sorter Core Facility.