| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology, Pathobiology, and Genetics |
Departments of 1 Internal Medicine, 2 Anatomy and Neurobiology, and 3 Biostatistics, Washington University School of Medicine, St. Louis, Missouri; 4 The Jackson Laboratory, Bar Harbor, Maine; and 5 Department of Pathology, Division of Molecular Diagnostics, University of Pittsburgh, Pittsburgh, Pennsylvania
Requests for reprints: Timothy Graubert, Division of Oncology, Stem Cell Biology Section, Washington University School of Medicine, Campus Box 8007, 660 South Euclid Avenue, St. Louis, MO 63110. Phone: 314-747-4437; Fax: 314-362-9333; E-mail: graubert{at}medicine.wustl.edu.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
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 (68). 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 |
|---|
|
|
|---|
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 pathogenfree 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 |
|---|
|
|
|---|
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).
|
|
40 weeks after ENU exposure (50 weeks of age; Fig. 3A
).
|
|
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 < 109) 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 < 106) 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.
|
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).
|
|
| Discussion |
|---|
|
|
|---|
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.
| Acknowledgments |
|---|
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.
| Footnotes |
|---|
Received 9/22/05. Revised 3/ 1/06. Accepted 3/ 9/06.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. K. Funk, T. J. Maxwell, M. Izumi, D. Edwin, F. Kreisel, T. J. Ley, J. M. Cheverud, and T. A. Graubert Quantitative trait loci associated with susceptibility to therapy-related acute murine promyelocytic leukemia in hCG-PML/RARA transgenic mice Blood, August 15, 2008; 112(4): 1434 - 1442. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Walter Del(5q): gene dosage matters Blood, July 15, 2007; 110(2): 473 - 474. [Full Text] [PDF] |
||||
![]() |
J. M. Joslin, A. A. Fernald, T. R. Tennant, E. M. Davis, S. C. Kogan, J. Anastasi, J. D. Crispino, and M. M. Le Beau Haploinsufficiency of EGR1, a candidate gene in the del(5q), leads to the development of myeloid disorders Blood, July 15, 2007; 110(2): 719 - 726. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Bogue, S. C. Grubb, T. P. Maddatu, and C. J. Bult Mouse Phenome Database (MPD) Nucleic Acids Res., January 12, 2007; 35(suppl_1): D643 - D649. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. A. Larson Etiology and Management of Therapy-Related Myeloid Leukemia Hematology, January 1, 2007; 2007(1): 453 - 459. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |