| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology, Pathobiology, and Genetics |
Departments of 1 Human Oncology and 2 Biostatistics and Medical Informatics, School of Medicine and Public Health and 3 Department of Genetics, University of Wisconsin-Madison, Madison, Wisconsin
Requests for reprints: Amy R. Moser, Department of Human Oncology, University of Wisconsin-Madison, K4/310 CSC Box 3684, 600 Highland Avenue, Madison, WI 53792. Phone: 608-265-6520; Fax: 608-263-9947; E-mail: armoser{at}wisc.edu.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
The transition from hyperplastic lesions to tumors is critical for the development of cancers (4). In the human, ductal and lobular hyperplasias are known risk factors for cancer development. For example, in one study, about 20% of women diagnosed with atypical lobular hyperplasia or lobular carcinoma in situ later developed cancer (5). Interestingly, this risk for cancer development applies to both the involved and contralateral breast in these women (5). This transition is poorly understood due to both the difficulty of identifying and tracking early hyperplastic lesions and the lack of good mouse models.
The ApcMin/+ mouse is a well-characterized mouse model that is predisposed to both intestinal and mammary tumors (6). Upon exposure to the alkylating agent ethylnitrosourea (ENU), >90% of B6-ApcMin/+ female mice develop mammary squamous cell carcinomas (SCC) within 65 days, with an average 3.3 tumors per mouse and an average latency of about 56 days. FVBB6 F1 ApcMin/+ mice are resistant to mammary tumors but susceptible to focal alveolar hyperplasias under the same treatment (7). This suggests that FVB carries alleles at modifier loci of ApcMin that act dominantly to affect the progression of hyperplasias to tumors. The identification of such modifiers could provide biomarkers of hyperplasias that could be helpful in the early diagnosis of breast cancer.
The inheritance of a mutation in APC gives predisposition to the autosomal dominant disorder familial adenomatous polyposis in humans (8, 9). Mouse APC shares 90% amino acid identity with human APC, and germline mutations result in similar phenotypes. Min, a dominant allele of Apc in the mouse, was generated by ENU germline mutagenesis and results in a truncation mutation at codon 850 (10). The insight that genetic background affects intestinal tumor development led to the mapping and molecular identification of Mom1 (modifier of Min 1; refs. 11, 12).
Wnt signaling has been observed to be abnormally activated in many cancers including breast cancer in both mice and humans. APC is critically important in the canonical Wnt signaling pathway and, together with GSK-3β and Axin, can participate in the phosphorylation of β-catenin, which is then targeted for degradation (13, 14). Upon Wnt activation or mutations of other components, β-catenin can translocate into the nucleus and bind with transcription factors, such as TCF4, to activate target genes including cylin D1 and c-myc (13, 14). Transgenic expression of β-catenin, cyclin D1, and c-myc in mammary tissue leads to mammary tumors in mice (15–17). Amplification of cyclin D1 was found in >50% human breast cancer, and somatic APC mutations were found in 18% breast cancer (18–20). The aberrant regulation of β-catenin is a prognostic marker in human breast cancer (21, 22). Overall, these data provide strong evidence that activation of canonic Wnt signaling is a common feature in human breast cancer.
In B6 ApcMin/+ female mice, SCC develop rapidly after ENU treatment with nearly all of the mice developing tumors by 65 days after treatment. In contrast, when treated with ENU, FVBxB6 ApcMin/+ female mice develop multiple alveolar hyperplasias and few SCC with long latency (7). Because the initiation step is most likely the same in both of these strains of mice, loss of the wild-type allele of Apc due to ENU-induced mutation, the difference in mammary tumor susceptibility must be due to genes that affect the process of tumor development. The resistance to mammary tumor development in FVBB6 F1 ApcMin/+ females could be due to the effect of modifier loci on any of the steps of tumor development. However, the susceptibility of FVBB6 F1 ApcMin/+ females to hyperplasias suggests the modifier loci do not affect the initiation of hyperplasia. Thus, the modifier loci most likely affect the transition of hyperplasias to mammary tumors or the progression and growth of mammary tumors. As a first step toward understanding the genetic control of mammary tumor development in these mice, we performed backcross analysis, mapping the traits of tumor number and tumor latency. Further molecular identification of these modifier genes would provide more insights into the regulation of the Wnt signaling pathway and mechanisms of breast cancer development.
| Materials and Methods |
|---|
|
|
|---|
For the FVBB6 backcross, wild-type FVB females were mated to B6 ApcMin/+ males (N59-66) to produce F1 females. The F1 Apc+/+ females (n = 72) were backcrossed with B6 ApcMin/+ males (n = 48) to produce N2 females. A total of 307 ApcMin/+ N2 females were produced and treated with 50 mg/kg ENU by i.p. injection when between 35 to 40 days of age and were euthanized
75 days after ENU treatment as described (7). Mice were palpated weekly by the same observer (D. Teske), and palpable tumors were recorded when first noted and then confirmed at necropsy. The total tumor number was determined at necropsy by noting all visible discrete masses. Tumors were collected and fixed in formalin for histologic analysis. We also collected 4-cm regions from the proximal, middle, and distal small intestine and the whole colon for counting intestinal tumors.
For the 129X1B6 backcross, wild-type 129X1 female mice (N = 42) were crossed to B6 ApcMin/+ males (N35-36; n = 7) to produce (129X1B6) F1 mice. Backcross mice were produced by crossing (129X1B6) F1 females to B6 ApcMin/+ males (n = 7), (129X1B6) F1 ApcMin/+ females to B6 male mice (n = 18), and B6 females to (129X1B6) F1 ApcMin/+ male mice (n = 55). A total of 80 129X1B6 backcross ApcMin/+ females were produced. The mice were ENU-treated and followed as described above for the FVB backcross.
Genotyping. For ApcMin/+ genotyping, genomic DNA was obtained after incubating tissues in 100 µL of 50 mmol/L NaOH at 95°C for 1 h and then adding 10 µL of 1 mol/L Tris (pH = 5), modified from the previous report (23). Animals were genotyped for ApcMin/+ by PCR using an allele-specific PCR assay (11).
For both backcross analyses, DNA was isolated from the spleens of backcross females using the DNeasy kit (Qiagen). Simple sequence length polymorphism (SSLP) primers were purchased from Research Genetics or Integrated DNA Technologies. For the FVBxB6 BC, a total of 102 SSLP markers with an average spacing of about 20 cM were genotyped for 180 mice selected from the ends of the tumor number distribution (90 mice with 1 or fewer tumors and 90 with 3 or more). Chromosomes with putative modifier loci were then genotyped with higher density markers for the remaining 127 mice. For the 129 backcross, all 80 mice were genotyped with 79 SSLP markers. Markers are available upon request. The PCR reactions were set up with 20 µL containing 50 to 100 ng genomic DNA, 2.0 µL 10 x buffer, 1.5 µL 25 nmol/L MgCl2, 30 nmol/L forward and reverse primers, and 0.5 µL 10 nmol/L deoxynucleotide triphosphates. The PCR conditions were 2 min at 94°C for 1 cycle; 15 s at 94°C, 45 s at 55°C, and 45 s at 72°C for 35 cycles; and finally, 7 min at 72° for 1 cycle. PCR products were resolved on 4% agarose gels and visualized by ethidium bromide staining.
Statistical analysis. The number of QTL controlling mammary tumor number was estimated as described in Dietrich et al. (11) with the classic formula of Wright, n = (uF1-uP)2/4
2G. Usually this formula underestimates the number of QTL because the assumptions are not exactly satisfied. MSTAT, a program provided by Dr. Norman Drinkwater of the McArdle Laboratory at the University of Wisconsin-Madison (Madison, WI), was used for the Kaplan-Meier analysis of tumor latency.
Linkage analysis was carried out using the parametric and nonparametric methods implemented in R/qtl (ref. 24; referred to hereinafter as R/qtl-P and R/qtl-NP, respectively) and the nonparametric methods implemented in the Q-link program, provided by Dr. Drinkwater. R/qtl-P implements a hidden Markov model to deal with missing genotypes; the expectation-maximization algorithm was used to perform single QTL genome scans (24). For Q-link, the statistic Zw is determined for each of the markers using the Wilcoxon rank-sum test and a logarithm of odds (LOD) score was calculated by the function LODW = 0.5 (log10 e) (ZW) (ref. 25); R/qtl-NP also uses a Wilcoxon rank-sum test but allows for LOD score evaluation between markers (24). For the linkage analyses of tumor latency, nonparametric methods were used because the raw and transformed phenotypes were not Gaussian. Permutation tests were performed (10,000 permutations per phenotype) to determine thresholds for suggestive (P < 0.05) and significant linkage (P < 0.01) at the genome level.
Two methods were used to assess the presence of interactions among QTL. We first considered a model selection procedure as detailed in Lan et al. (26). That procedure identifies potential interactions among QTL and other genome regions that perhaps do not show significant main effects. In short, the procedure first identifies putative QTL using a LOD score profile. In the initial step, we do not require that putative QTL be statistically significant; they are defined as those with LODs >1, where LOD scores are calculated using R/qtl-P. We then consider all possible models allowing for additive effects among the putative QTL and pairwise interactions. The Bayes Information Criterion (BIC; ref. 27) is used to score each model. The BIC balances goodness-of-model fit with the number of model variables. The model with the best (lowest) BIC is then identified. We also considered the two dimensional scans provided by R/qtl. There, each pair of loci is evaluated allowing for main effects and the pairwise interaction between those effects. The likelihood of the full model is then compared with alternatives to test for additive and interaction effects (24).
To combine the FVBB6 and 129X1B6 backcrosses (28), the FVB and 129X1 alleles were assigned the same genotype and the B6 allele was treated as the other genotype. Therefore, a binary allelic pattern was used, and only loci where both FVB and 129X1 carry alleles that affect tumor development in the same direction will be identified as modifiers.
| Results |
|---|
|
|
|---|
By the end of the experiment, 57% of backcross females had developed at least one palpable tumor (meaning that the tumor was detected before necropsy). The classification of mammary tumors as palpable tumors is dependent on both tumor size and location. However, because some tumors did develop in each of the mammary glands in the set of backcross mice, it is expected that tumors in all mice had approximately equal chances of detection. Among a total of 780 visible mammary tumors in the backcross females, about one third (259 of 780) were palpable and two thirds (521 of 780) were detected only at dissection. On average, the backcross ApcMin/+ females developed 0.84 palpable tumors and 1.70 nonpalpable tumors.
Identification of modifiers affecting tumor number in the FVBB6 backcross. A whole genome scan with 102 SSLP markers spaced an average of 20 cM apart was performed with mice at the two extremes of the tumor number distribution: 90 mice with 0 or 1 tumor and 90 mice with <3 tumors (Fig. 1A ). The rest of the mice were then genotyped only on chromosomes with a LOD score of >1. Three loci, on chromosomes 4, 6, and 9, were associated with mammary tumor number. These QTL are designated Mmom1 (mammary modifier of Min; chromosome 9), Mmom2 (chromosome 4), and Mmom3 (chromosome 6) in order of effect. The summary of the QTL analysis is shown in Table 1 , and the LOD profile for each chromosome is shown in Fig. 1B.
|
|
|
50% compared with mice homozygous at both loci (Fig. 2A; 1.7 ± 1.5 versus 3.2 ± 2.0). Mmom3 also has an additive effect in concert with Mmom2 (Fig. 2B; 3.3 ± 2.2 versus 2.0 ± 1.6). Taken together, these analyses indicate that although the locus on chromosome 6 fails to achieve significance in the one-dimensional scans, this locus does influence tumor number in this backcross.
|
A whole genome scan was carried out as with tumor multiplicity (using the mice on the ends of the tumor number distribution) and chromosomes with LOD > 1 were selected for further genotyping and analysis using all of the mice. Modifier loci on chromosomes 4 and 18 were found to significantly affect tumor latency (Fig. 1B and Table 1). The modifiers on chromosomes 6 and 9 did not show a significant effect on tumor latency. The peak position of the QTL affecting latency is the same as the QTL affecting tumor number on chromosome 4 (D4mit82 and Mmom2). This suggests that Mmom2 could affect both tumor number and tumor latency. Less than 50% of mice carrying an FVB allele at Mmom2 developed a palpable mammary tumor by the end of experiment (75 days after ENU), whereas 50% of mice homozygous for B6 alleles had palpable tumors by 57 days after ENU treatment. A locus on chromosome 18 at D18Mit24 (25 cM) specifically affects tumor latency and was designated Mmom4. This locus shows no effect on the tumor multiplicity. This locus maps about 10 cM distal of the Apc locus, which also maps to chromosome 18.
To determine the effect of each modifier on tumor development, we performed Kaplan-Meier analysis stratifying the mice by genotype at D4Mit82 or D18Mit24. The tumor latency for the mice homozygous at each marker is significantly different from the tumor latency for the mice heterozygous at each marker (Fig. 3A and C ). We then stratified for both markers and found additive effects between Mmom2 and Mmom4 (Fig. 3D). When mice were homozygous for B6 alleles at both loci, around 80% of females developed a palpable tumor. When the mice were heterozygous at both loci, only 40% of the mice developed a palpable tumor (P < 0.001). Heterozygosity for either one of the loci has approximately the same effect on tumor latency.
|
As the data from both the FVBB6 and the 129B6 backcrosses were collected in a similar way, we combined these two backcrosses to improve the power of modifier mapping (28). Because different markers are used in these two backcrosses, 1-cM interval mapping was used to calculate the LOD score for each modifier with the software R/qtl. The hypothesis is that if a modifier is shared between different crosses, the combined cross can improve the LOD score and would dilute the effect of cross-specific modifier QTL. Figure 1D shows the LOD profile in the combined analysis and the results are summarized in Table 2. On chromosome 4, the LOD score in the region of Mmom2 was highly significant (P < 0.001) for both tumor number and tumor latency. This suggests that the resistant allele of Mmom2 is shared by the FVB and 129X1 strains. On chromosome 6, the LOD score of Mmom3 reached suggestive significance (P < 0.05) for tumor number with no effect on tumor latency. This provides more evidence that there is a shared modifier on chromosome 6 affecting tumor number. On chromosome 9, the LOD score of Mmom1 decreased compared with FVBB6 backcross and there was no effect on tumor latency. This suggests that the 129X1 and B6 strains share alleles with similar effect at Mmom1. On chromosome 18, the LOD score of Mmom4 was highly significant (P < 0.001) for tumor latency, indicating that FVB and 129X1 strains share alleles at Mmom4. In summary, the combined cross analysis confirmed the results of individual backcross analyses and provided more support for the effect of Mmom2, Mmom3, and Mmom4 on tumor development.
Tissue specificity of modifiers. Because ApcMin/+ mice develop both mammary and intestinal tumors, we could test for the effect of any modifiers on both tumor types. None of the Mmom loci affected intestinal tumor development (Table 3 ). D4Mit13 (71 cM), the most distal marker on chromosome 4 in the FVBB6 backcross, is significantly associated with intestinal tumor number but has no effect on mammary tumor number and latency. D4Mit13 maps to the same position as Mom1 (11). This suggests that the FVB strain carries a resistant allele of Mom1 (11). A database search identified a SNP (rs27560348, T/C, and dbSNP126) between B6 and FVB in the phospholipase gene Pla2g2a. This SNP causes a frameshift mutation and generates a mutant allele for Pla2g2a in B6 mice (12). Taken together, these data provide strong evidence that the modifiers of ApcMin show tissue specificity for mammary tumors and intestinal tumors.
|
| Discussion |
|---|
|
|
|---|
In this analysis, we have identified modifier loci that affect either mammary tumor number or/and tumor latency. Mmom1 on chromosome 9 affects only mammary tumor number, suggesting that it may affect the probability of transition from hyperplasia to tumor. Mmom2 on chromosome 4 affects both tumor number and tumor latency, suggesting that it may affect both the probability of transition from hyperplasia to tumor and the rate of tumor growth. This effect is similar in both backcrosses investigated. Mmom3 on chromosome 6 seems to cooperate with Mmom1 or Mmom2 to affect tumor number. Mmom4 on chromosome 18 affects mammary tumor latency, suggesting that it may affect mammary tumor growth. The additive interaction between Mmom2 and Mmom4 suggests the genes underlying these two loci do not act in the same pathway.
Humans carrying germ-line mutations on APC usually develop colon cancer but very few other cancer types, implying that there are tissue-specific effects of loss of APC function. Our previous data suggest that there is no association between susceptibility to mammary tumors and intestinal tumors in F1 ApcMin/+ mice and raises a question whether there are tissue-specific modifiers of ApcMin (7). Our current data show that the susceptibility to intestinal tumors in this backcross is mostly controlled by Mom1 on chromosome 4, which has no effect on mammary tumor number or tumor latency. Conversely, the mammary modifiers show no association with intestinal tumor multiplicity. Alina et al. (31) reported that recombinant congenic strains (RCS) carrying the ApcMin allele show opposite effects on intestinal and mammary tumors. These opposite effects could be due to one modifier locus with opposite effects in each tissue or the effects of tissue-specific modifiers. Our data favor the hypothesis that the observed phenotype in RCS is caused by different modifiers. The identification of tissue or tumor type–specific modifiers may aid in the identification of the genes underlying the modifier loci based on function. This also indicates the complexity of the phenotype conferred by loss of Apc function on different tissues.
Interaction between modifiers is frequent in colon cancer, lung cancer, and skin cancer (32–35). However, to identify low-penetrance disease-related genes presents a challenge in human and mouse genetics. The powerful genetic manipulation possible in mice makes it feasible to identify these low-penetrance genes. These low-penetrance genes might only be identified when they interact with major QTL or other low-penetrance loci to achieve a statistically significant effect. In our experiments, we used backcrosses to reduce the genetic complexity (36) because of the narrow mammary tumor number range in B6-ApcMin/+, FVBB6 F1-ApcMin/+, and 129X1B6 F1 ApcMin/+ mice. Our results reveal that additive effects between modifiers not only affect mammary tumor number but also regulate mammary tumor latency. These results suggest that multiple loci control each stage of the development of mammary tumors and distinct stages are controlled by different modifier loci.
Only a few QTL affecting mammary tumor susceptibility have been mapped in mice, all of which are in the mouse mammary tumor virus (MMTV)-PyVT transgenic mouse model (37, 38). Mmom1 maps near Apmt2, which interacts with Apmt1 (on chromosome 15) to affect tumor latency but not tumor number (37) in MMTV-PyVT mice. Mmom2 maps near Mmtg1, which affects total mammary tumor weight but not tumor number or tumor latency (38). However, the FVB alleles confer sensitivity in the MMTV-PyVT model, whereas the FVB alleles confer resistance in the ApcMin/+ mouse model. As they affect different phenotypes with different susceptibility/resistant status, these modifier loci most likely represent different genes, although they map to a similar location. Several QTL affecting various aspects of mammary tumor development in the rat have been mapped to the orthologous regions to the loci we have mapped on chromosomes 4. 9, and 18 (39, 40). Some of the rat and mouse loci that map to orthologous regions seem to have similar effects on tumor development. However, given the imprecision of the mapping and the differences in the models, it is premature to conclude that we have mapped the same genes. However, with more precise mapping and further phenotypic analysis, it will be possible to determine whether the same genes have indeed been identified.
Mmom1 maps to the distal part of chromosome 9, which is orthologous to human 3p21-23. This region undergoes frequent loss of heterozygosity (LOH) and hypermethylation in lung cancer and breast cancer (41, 42). In addition, the Mmom2 region is orthologous to human chromosome 9q, which has been shown to undergo LOH in about 40% of breast cancers (43). Therefore, it is possible that the loci we have identified may have a role in the development of human cancer.
B6 and FVB are two strains that are frequently used for mouse models of cancer research. As such, information about loci that alter the course of tumor development or progression can possibly be applied to other models. The identification of QTL is the first important step toward understanding the development of mammary tumors, and the characterization of modifier genes could provide valuable insights into the regulation of mammary tumors. These modifier genes could be biomarkers for early stages of breast cancer and improve the diagnosis of breast cancer.
| 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 Laura Hegge and Jared Finger for assistance with the mouse colony and backcross experiments; Dr. Norman Drinkwater for the Q-Link program and for the advice; Dr. Karl Broman for the advice on using R/qtl; Dr. Ruth Sullivan for the assistance with pathology; Meng Chen for the help with statistical analyses; and the University of Wisconsin Comprehensive Cancer Center Histology Core for assistance with processing and sectioning of tumor samples.
Received 7/ 3/07. Revised 9/27/07. Accepted 10/ 8/07.
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 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 |