Heterogeneity of mammary lesions represent molecular differences
© Namba et al; licensee BioMed Central Ltd. 2006
Received: 31 July 2006
Accepted: 05 December 2006
Published: 05 December 2006
Human breast cancer is a heterogeneous disease, histopathologically, molecularly and phenotypically. The molecular basis of this heterogeneity is not well understood. We have used a mouse model of DCIS that consists of unique lines of mammary intraepithelial neoplasia (MIN) outgrowths, the premalignant lesion in the mouse that progress to invasive carcinoma, to understand the molecular changes that are characteristic to certain phenotypes. Each MIN-O line has distinguishable morphologies, metastatic potentials and estrogen dependencies.
We utilized oligonucleotide expression arrays and high resolution array comparative genomic hybridization (aCGH) to investigate whole genome expression patterns and whole genome aberrations in both the MIN-O and tumor from four different MIN-O lines that each have different phenotypes. From the whole genome analysis at 35 kb resolution, we found that chromosome 1, 2, 10, and 11 were frequently associated with whole chromosome gains in the MIN-Os. In particular, two MIN-O lines had the majority of the chromosome gains. Although we did not find any whole chromosome loss, we identified 3 recurring chromosome losses (2F1-2, 3E4, 17E2) and two chromosome copy number gains on chromosome 11. These interstitial deletions and duplications were verified with a custom made array designed to interrogate the specific regions at approximately 550 bp resolution.
We demonstrated that expression and genomic changes are present in the early premalignant lesions and that these molecular profiles can be correlated to phenotype (metastasis and estrogen responsiveness). We also identified expression changes associated with genomic instability. Progression to invasive carcinoma was associated with few additional changes in gene expression and genomic organization. Therefore, in the MIN-O mice, early premalignant lesions have the major molecular and genetic changes required and these changes have important phenotypic significance. In contrast, the changes that occur in the transition to invasive carcinoma are subtle, with few consistent changes and no association with phenotype.
We propose that the early lesions carry the important genetic changes that reflect the major phenotypic information, while additional genetic changes that accumulate in the invasive carcinoma are less associated with the overall phenotype.
The paradigm that cancer progression is a multi-step process, associated with multiple molecular changes as it progresses from preneoplasia to invasive carcinoma , has been challenged by recent molecular data. Gene expression profiling and comparative genomic hybridization (CGH) studies of breast cancer demonstrate that early stages in the human breast cancer such as ductal carcinoma in situ (DCIS), a precursor lesion for invasive carcinoma, has most, if not all, of the molecular attributes of the corresponding invasive carcinoma despite the distinct pathological characteristics [2–5]. This is contrary to the multi-step paradigm that centers on cumulative molecular aberrations with progression. These data suggest an alternative view that the early lesions are already equipped with the molecular changes responsible for tumorigenesis, despite the disparate histological characteristics between the early lesions and the invasive carcinoma.
Breast cancer can be viewed as a heterogeneous disease histopathologically, as well as molecularly. Molecular profiling studies of breast cancer have shown that tumors can be classified into subtypes based on their expression patterns [6–8]. Pathologically, breast lesions are classified by different categories, such as estrogen receptor (ER) status, Her2 status and the degree of differentiation (tumor grade). Both ER and Her2 status are important prognostic factors and portend how a lesion responds to different therapeutic strategies. Both DCIS and invasive ductal carcinoma (IDC) are categorized into three tumor grades . DCIS lesions are also classified into different subtypes by their histological morphology . It is believed that the different classes of lesions have common characteristics that reflect certain clinical outcomes. Gene expression studies of different pathological stages of breast cancer have shown that different tumor grades are associated with distinct expression signatures, confirming the molecular basis for the differences in pathological classification; the same studies show that the profiles of the different stage lesions also have extensive similarities, suggesting that there may not be as many molecular changes associated with tumor progression as was previously thought .
Over the last several years, the MIN-O (mammary intraepithelial neoplasia outgrowth) mouse model has been shown to parallel various aspects of human breast cancer development [11, 12]. MIN (mammary intraepithelial neoplasia) is an early mammary lesion that satisfies the operational definition of premalignancy . The MIN-O mouse was established by transplanting a MIN lesion from a young polyomavirus middle-T (PyVmT) transgenic female mammary fat pad to a host fat pad . The transplanted MIN lesion will grow from the transplanted location to fill the host fat pad. The growth from the transplanted MIN lesion is called MIN-outgrowth (MIN-O). Since the MIN-O mouse eventually develops invasive carcinoma within the MIN-O, tumor in this context refers to invasive carcinoma that emerges within the MIN-O tissue. Different lines of MIN-O mice, which were established by transplanting MIN lesions to separate hosts, have characteristic histopathology, metastatic potential, and tumor latency, suggesting that MIN lesions have different biological outcomes despite the same initiating oncogene (PyVmT) . In addition, in depth gene expression studies in one MIN-O mouse line, demonstrated that the MIN-O and tumor tissues show strikingly similar expression profiles between each early premalignant lesion and its corresponding tumor pair . Thus, this model recapitulates the foundation of human DCIS and the heterogeneous nature of the disease.
In the current study, we assessed the correlation between certain biological outcomes to molecular changes associated with four different lines of MIN-O mice. In addition to the gene expression analysis, we report for the first time the utilization of high resolution array comparative genomic hybridization (aCGH) to investigate whole genome aberrations in both the MIN-O and tumor tissue. In addition to this whole genome analysis (35 kb resolution), interstitial deletions were verified with a custom made array designed to interrogate the specific regions at approximately 550 bp resolution. We demonstrate that expression and genomic changes are present in the early premalignant lesions and that these molecular profiles can be correlated to phenotype (metastasis and estrogen responsiveness). Progression to invasive carcinoma was associated with few additional changes in gene expression and genomic organization. Therefore, in the MIN-O mice, early premalignant lesions have the major molecular and genetic changes required and these changes have important phenotypic significance. In contrast, the changes that occur in the transition to invasive carcinoma are subtle, with few consistent changes and no association with phenotype. Therefore, we propose that the early lesions carry the important genetic changes that reflect the major phenotypic information, while the additional genetic changes accumulated in the invasive carcinoma are less associated with the overall phenotype.
Standard techniques for mammary gland clearing and mammary tissue transplantation, and the establishment of the MIN-O lines, have been described previously [11, 12]. Three-week-old FVB female mice were purchased from Charles River Laboratories (Wilmington, MA). Animals were maintained in a vivarium at the University of California Davis campus under IACUC approved protocols.
To test MIN-O ovarian hormone sensitivity, a 1 mm3 piece of a MIN-O tissue from each line was transplanted in the gland-cleared #4 fat pads of 3-week-old virgin FVB female mice. At the time of transplantation, the experimental group (ovx, n = 8) was ovariectomized, while the control group (intact, n = 4) received mock surgery without ovariectomy. At 5 weeks after transplantation, the MIN-O sizes were measured by visually inspecting the transplantation site under a dissecting microscope. Animals were palpated weekly until a tumor was detected. The experiment was concluded at day 99 post-transplantation. For the tumor ovarian hormone sensitivity experiment, 1 mm3 pieces of tumor tissue that developed from various MIN-O lines were transplanted to intact and ovx hosts (at least 2 animals/4 tumors each for a tumor). Both intact and ovx animals were terminated at the same time and the tumor bearing fat pads were weighed.
Immunohistochemistry of formalin fixed tissues using anti-estrogen receptor was performed as previously described . MIN-O transplanted fat pads from the ovariectomy experiment were used to determine the ER localization in the MIN-Os from the intact and ovariectomized hosts.
Mouse Genome CGH Microarrays
The Agilent Mouse Whole Genome CGH Microarray Kit 44A (G4414A, Agilent Technologies, Santa Clara, CA) consists of roughly 43,000 60-mer oligonucleotide probes that cover both exonic and intronic sequences. These arrays were used to analyze the MIN-Os and tumors from each line. These high resolution arrays have an average spatial resolution of 35 kB, after removal of Repeat Masked (-RM) . Agilent uses SurePrint® technology  to print this entire content on a standard 1" × 3" glass slide.
Regions selected for high resolution custom CGH microarrays
Probe Density (-RM)
DNA was extracted from frozen MIN-O/tumor pairs isolated from the same fat pads of MIN-O lines 4w-4 (n = 5 pairs), 4w-11 (n = 3 pairs), 8w-B (n = 3 pairs) and 8w-D (n = 3 pairs). DNA from FVB/NCrl female livers was used as the reference sample. Samples were prepared using the direct labeling method, according to the Oligonucleotide Array-Based CGH for Genomic DNA Analysis Protocol (Version 2.0, available at ), with minor modifications of the protocol. A total of 2 ug of digested purified DNA was labeled, using the BioPrime Array CGH Genomic Labeling kit (Invitrogen, Carlsbad, CA), for both the experimental and reference channel. During the labeling, cyanine 3-dUTP and cyanine 5-dUTP (Perkin Elmer, Wellesley, MA) were incorporated in either the experimental or the reference sample. For each sample, two arrays were utilized in a dye-swap pair format, in which the experimental sample and reference sample are each labeled separately with cyanine 3 and cyanine 5, respectively, and combined for the first array, and then labeled in the reverse and combined for the second array. After labeling, the corresponding experimental and reference samples were combined and unincorporated nucleotides were removed on a Microcon YM-30 apparatus (Millipore, Billerica, MA). The purified sample was then hybridized onto a single CGH array in a solution containing Mouse Cot-1 DNA (Invitrogen, Carlsbad, CA), Agilent Blocking Agent, and Agilent Hybridization Buffer (Agilent, 5188–5226). Hybridization took place in a rotisserie oven set to 65°C with a rotation speed of 20 rpm. The washing and scanning of the slides was performed in an ozone-free area to prevent the degradation of the cyanine 5 dye. Each slide was washed according to the manufacturer's recommendations with no Stabilization and Drying Solution. The slides were scanned on an Agilent DNA Microarray Scanner B.
TIFF image files generated from the scanner were analyzed using Agilent's Feature Extraction Software 8.1. This software allows a standardized method of extracting data using a standard CGH normalization protocol based on all array features (44k_CGH_0605). The QC report generated was checked to ensure proper hybridization and placement of the grid by the software. Data generated by the Feature Extraction tool were loaded into CGH Analytics 3.1 (Agilent Technologies) to allow visualization of the data. Dye swap pairs were combined to facilitate the visualization of chromosomal aberrations. A smoothing average of 50 MB was used for the full genome arrays. Regions of aberration (deleted or amplified) were determined by the methodology of Barrett et al..
Total RNA extraction and biotin-labeled probe synthesis for gene expression analysis has been described previously [11, 12]. Frozen tissues from MIN-O lines 4w-4 (n = 3 MIN-Os, n = 3 tumors), 4w-11 (n = 3 MIN-Os, n = 3 tumors), 8w-B (n = 4 MIN-Os, n = 4 tumors) and 8w-D (n = 4 MIN-Os, n = 3 tumors) and intact mammary fat pads from 16-day prelactating FVB females (n = 2) were used for hybridization onto Murine Genome U74Av2 GeneChips (Affymetrix, Santa Clara, CA). The scanned images were processed using GCOS 1.4, and then the data was analyzed with GC-RMA in the ArrayAssist software package (Stratagene). GCRMA derived expression values were log-transformed and analyzed in Arrayassist (t-test, unequal variance). Expression comparisons of all MIN-Os to prelactating mammary glands, and all tumors to MIN-Os were carried out in Arrayassist using t-test with FDR correction. 938 probes were differentially expressed with p < 0.01 in the MIN-Os vs prelactating comparison. 55 probes were differentially expressed with p < 0.01 in the tumors vs MIN-Os comparison. For expression comparisons of samples from individual lines, profiles were compared in Arrayassist using t-test without FDR correction and gene lists were generated with criteria of p < 0.01.
Gene lists were subjected to pathway and function analysis in Ingenuity Pathways Analysis (Ingenuity Systems, Redwood City, CA). Annotation and gene ontology of each probe was assigned by NetAffyx (Affymetrix, Inc). The results from the statistical analyses described above were analyzed by ErmineJ software  to determine the number of significant gene ontology (GO) categories using gene score resampling (GSR) . Unfiltered t-test results were imported into Erminej and analyzed using a 100,00 fold iteration re-sampling approach to determine which Biological Processes categories were significant below a Benjamini-Hochberg FDR-corrected p-value of 0.01.
For correlating the genomic changes to the expression differences, chromosome locations of probes and known genes were determined by the UCSC Genome Browser . For the analysis of whole chromosome gain of chromosome 2, expression values of chromosome 2 probes (n = 871) from 8w-B and 4w-4 MIN-Os were compared to those from 8w-D and 4w-11 MIN-Os (t-test, unequal variance). For chromosome 1, 10 and 11 probes (n = 652, 481, and 957, respectively), expression values from 8w-B MIN-Os were compared to those from 8w-D, 4w-4, and 4w-11 (t-test, unequal variance). For the analysis of 2F, 3E, 17E and 11B, the corresponding probes were identified in UCSC Genome Browser. For 2F, 3E and 17E probes, the average expression values of the 8w-D MIN-Os were compared to the average expression values of 8w-B, 4w-4 and 4w-11 MIN-Os. For 11B probes, expression values of the 5189 MIN-O sample were compared to those of the 5207 MIN-O sample.
MINs are heterogeneous
Characteristics of MIN-O lines
Major histological characteristics of MIN-Os*
ER status of MIN-Os
Solid chords & nests
Negative or cytoplasmic
Microacini & organized cysts
Dysplastic cysts. Better differentiated. Microacinar structures.
Prominent and distinct microacinar patterns. Organized cysts. Patches of eosinophilic microacini.
Ovarian hormone sensitivity of MIN-Os
Effect of ovariectomy on MIN-O growth and tumor latency
MIN-O size at 5-weeks posttransplantation
Tumor latency (in weeks)
Control %FPF (n)
% FPF (n)
relative to control
55.00 ± 10.00 (4)
24.29 ± 9.76 (7)
37.50 ± 17.08 (4)
15.00 ± 18.17 (6)
67.50 ± 12.58 (4)
22.86 ± 7.56 (7)
60.00 ± 0.00 (4)
10.71 ± 4.50 (7)
ER expression in the 4w-4, 4w-11, 8w-B and 8w-D MIN-Os were determined by immunohistochemistry of the MIN-O transplanted fat pads (Table 2). In the intact hosts, patches of strong nuclear ER-positive cells were present in the 8w-B MIN-Os, while in the 8w-D MIN-Os, nuclear ER-positive cells were also present, but with less intensity. In 4w-4 and 4w-11 MIN-Os from the intact hosts, nuclear ER staining was rarely found and, if any, weak ER staining was seen in the cytoplasm. In all four lines, MIN-Os from the ovariectomized hosts had no nuclear ER-positive cells, and reduced cytoplasmic staining, if any, compared to those from the intact hosts (data not shown).
To study genomic aberrations associated with MIN-O and tumor samples, we utilized Agilent mouse CGH arrays on 14 MIN-O/tumor pairs. Genomic DNA from at least 3 different MIN-Os and the corresponding tumors from the same fat pad were isolated from 4 different MIN-O mouse lines.
In addition to the whole chromosome gains, the MIN-Os also exhibited chromosomal losses and/or gains of small chromosomal regions. The most common regions deleted were 2F1-2 (21.4% MIN-Os and tumors), 3E4 (7.1% MIN-Os, 21.4% tumors), and 17E2 (7.1% MIN-Os, 21.4% tumors). Other regions that were found to have deletions in one MIN-O/tumor pair were 12A, 16C2, and 3D-E1. Also, the deletion of a small region in 19D was found in a single tumor.
Small interstitial gains (duplication) of a chromosome material were not as common as interstitial chromosome deletions. The only recurring small gain was 11B1.2 (7.1% MIN-Os, 14.2% tumors). Another area of interest was duplication of 11E1 found in one tumor sample. This area has been reported to be amplified in Tg(PyVmT) tumors .
Expression changes associated with tumorigenesis
To understand the expression changes associated with mammary tumorigenesis, as well as the molecular basis of the heterogeneity of mammary lesions, we studied gene expression profiles of MIN-Os and tumor samples. Highly proliferative wild-type mammary gland (prelactating) was used as the control, as it has the same proliferative index as the early preneoplastic tissue (in contrast to the non-proliferative quiescent mammary fat pad). The tissues were analyzed with the Affymetrix GeneChip MG-U74Av2 oligonucleotide arrays.
The gene expression profiles of the MIN-Os were dramatically different from that of the proliferating mammary glands. 937 probes (899 genes – comprised of 19 unknown probes, 103 unknown genes or ESTs, and 777 known genes) were differentially expressed (p < 0.01). 293 out of the 777 known genes were down-regulated in MIN-Os and 484 out of the 777 genes were over-expressed in MIN-Os. The majority of these genes were involved in cell cycle and/or cell death functions [see Additional file 1].
Tumor vs MIN-O expression differences
Probe Set ID
FC in tumors
GO biological function
adiponectin, C1Q and collagen domain containing
Lipid metabolism, Organic acid metabolism
adrenergic receptor, beta 3
amine oxidase, copper containing 3
ankyrin repeat domain 47
anthrax toxin receptor 2
carbonic anhydrase 3
Lipid metabolism, Organic acid metabolism
cell death-inducing DFFA-like effector c
complement component 1, r subcomponent
C-type lectin domain family 3, member b
cysteine dioxygenase 1, cytosolic
Organic acid metabolism
cytochrome P450, family 2, subfamily e, polypeptide 1
early B-cell factor 1
fatty acid binding protein 4, adipocyte
fibroblast activation protein
flavin containing monooxygenase 1
FXYD domain-containing ion transport regulator 1
growth hormone receptor
histidine ammonia lyase
Organic acid metabolism
insulin-like growth factor 1
insulin-like growth factor binding protein 6
laminin, alpha 2
lipase, hormone sensitive
phosphoenolpyruvate carboxykinase 1, cytosolic
Lipid metabolism, Organic acid metabolism
plasma membrane associated protein, S3-12
protein phosphatase 2, regulatory subunit B (B56), alpha isoform
retinoic acid receptor responder (tazarotene induced) 2
retinol binding protein 4, plasma
RIKEN cDNA 4631408O11 gene
serine (or cysteine) peptidase inhibitor, clade A, member 3C
SNF related kinase
solute carrier family 1 (neutral amino acid transporter), member 5
sorbin and SH3 domain containing 1
stearoyl-Coenzyme A desaturase 1
Lipid metabolism, Organic acid metabolism
sulfotransferase family 1A, phenol-preferring, member 1
thioredoxin interacting protein
zinc finger homeobox 1a
Molecular changes associated with each MIN-O line
Recurring genomic aberration by MIN-O lines
When CGH results were analyzed for each MIN-O line, we found recurring genomic changes that were specific to certain lines. MIN-Os from 8w-B and 4w-4 mice carried large chromosome aberrations, frequently including multiple chromosome gains, while large chromosome aberrations were rare in the 4w-11 and 8w-D lines (Figure 1A&E). 100% of 8w-B MIN-Os carried whole chromosome gain of chromosomes 1, 2, 10 and 11. 80% of 4w-4 MIN-Os carried whole chromosome gain of chromosome 2. Tumors arising from both MIN-O lines (8w-B and 4w-4) often acquired additional chromosome gain(s) (6 out of 8 tumors), although the majority of the cases were limited to a gain of one chromosome, and there were no recurring chromosome gains found in tumors, except for chromosome 15.
Chromosome 2 gene expression in MIN-O lines with whole chromosome gains (8w-B & 4w-4 vs. 8w-D & 4w-11)
# of probes
Chromosome 1, 10 and 11 gene expression differences in MIN-O lines with whole chromosome gains (8w-B vs. 8w-D, 4w-4 & 4w-11)
# of probes
In-depth analysis of the data found that all the 8w-D samples carried deletions of 2F1-2, 3E4, and 17E2 regions. A pair of MIN-O and tumor samples from one 8w-D animal (5189) had the same amplification of 11B1 region. We also found a large amplification of 11E in a tumor from one 4w-4 animal. Among the four MIN-O lines, we found the least amount of chromosome abnormality in 4w-11 samples. 4w-11 samples did not carry any whole chromosome aberrations, or large duplications. In one MIN-O/tumor pair of 4w-11, we found a small deletion in the 3D-E1 region.
Expression in the deleted region (2F1-2, chr2: 130,621,986-131,735,939)
U74Av2 probe ID
lower expression in 8w-D MIN-Os*
Expression in the duplicated region (11B, ch11:42,631,278-44,744,905)
U74Av2 probe ID
higher expression in 5189MINO**
Gene expression associated with each MIN-O line
Significant biological functions for MIN-O vs PL expression differences for 4w-4 line
Highly significant GO functions
DNA-dependent DNA replication
organic acid biosynthesis
chromatin assembly or disassembly
regulation of cyclin dependent protein kinase activity
Significant biological functions for MIN-O vs PL expression differences for 4w-11 line
regulation of I-kappaB kinase/NF-kappaB cascade
positive regulation of I-kappaB kinase/NF-kappaB cascade
I-kappaB kinase/NF-kappaB cascade
antigen presentation, exogenous antigen
response to xenobiotic stimulus
response to bacteria
regulation of immune response
defense response to bacteria
Reponse to radiation/UV
response to radiation
detection of light
response to light
Significant biological functions for MIN-O vs PL expression differences for 8w-B line
mitotic spindle organization and biogenesis
spindle organization and biogenesis
microtubule cytoskeleton organization and biogenesis
microtubule polymerization or depolymerization
M phase specific microtubule process
regulation of mitosis
cell cycle checkpoint
interphase of mitotic cell cycle
G1/S transition of mitotic cell cycle
negative regulation of cell cycle
G2/M transition of mitotic cell cycle
regulation of cyclin dependent protein kinase activity
G1 phase of mitotic cell cycle
traversing start control point of mitotic cell cycle
cell cycle arrest
blood vessel morphogenesis
blood vessel development
Significant biological functions for MIN-O vs PL expression differences for 8w-D line
Response to radiation/UV
detection of light
response to radiation
response to light
response to UV
complement activation, classical pathway
positive regulation of immune response
regulation of immune response
cellular defense response (sensu Vertebrata)
antimicrobial humoral response
positive regulation of lymphocyte activation
symbiotic interaction between host and other organism
oxygen and reactive oxygen species metabolism
response to xenobiotic stimulus
T-helper 1 type immune response
detection of biotic stimulus
cytokine and chemokine mediated signaling pathway
antigen presentation, endogenous antigen
blood vessel morphogenesis
blood vessel development
extracellular structure organization and biogenesis
extracellular matrix organization and biogenesis
I-kappaB kinase/NF-kappaB cascade
regulation of I-kappaB kinase/NF-kappaB cascade
positive regulation of I-kappaB kinase/NF-kappaB cascade
Wnt receptor signaling pathway
Ras protein signal transduction
inactivation of MAPK
Although the comparison of all tumors vs. all MIN-Os identified a small number of down-regulated genes, there was no one gene common to the tumor vs. MIN-O gene lists from the four MIN-O lines [see Additional file 2]. Because the expression profiles of tumors were very similar to the corresponding MIN-Os, GO analysis did not result in any significant functional terms. However each line had characteristic expression differences of genes that suggest changes in tumor-associated functions. Tumors developed from 4w-4 MIN-Os were characterized by down-regulation of immunoglobulin and up-regulation of actin and myosin (3- to 4-fold change). Tnfrsf11/RANK, receptor activator of NfκB , and RalA, which is a critical component for Ras-induced tumorigenesis  were also up-regulated in these tumors. In 4w-11 tumors, CEBPA, which is associated with breast cancer progression  was down-regulated and CEBPB, a gene found to be up-regulated in breast cancer as well as mouse mammary tumors , is up-regulated. Moreover, the 4w-11 tumors exhibited down-regulation of actin, myosin, Ig, and MHC class 2 transcripts, which have been shown to be decreased in highly metastatic breast cancer cell lines , as well as decreased during progression to tumor . 8w-B tumors had up-regulation of PLK  and 8w-D tumors had up-regulation of Epha2 .
Molecular changes associated with biological characteristics
Molecular changes in MIN-Os with genetic instability (lines 8w-B and 4w-4)
dysregulation of mitotic spindle and chromosome segregation related genes To further understand the basis of chromosome instability in the MIN-Os, we assessed gene expression profiles of the genetically unstable MIN-O lines – 8w-B and 4w-4. Seven out of eight MIN-O samples from these two lines had whole chromosome gain of chromosome 2. The two lines had 13 genes from chromosome 2 in common in the gene expression comparison between MIN-O and PL [see Additional file 3]. Among these 13 genes, all but one was over-expressed in the MIN-Os. Aurora kinase A (Aurka/STK6/STK15) was the most dramatically over-expressed (9- to 10-fold).
Aurka regulates centrosome duplication and mitotic spindle assembly. Other genes involved in spindle assembly (Ran, Tuba1, Tubb6, Kif22, Nde1, Ndrg1) and centrosome duplication and sister chromatid segregation (Plk1, Gadd45a, Apc11, Cdc20 and Cdkn1a) were also dysregulated. In addition, genes up-regulated in 8w-B and 4w-4 MIN-Os included a large group of E2F4 regulated genes (H2AFZ, RFC2, RFC4, ZNF305, KIF22, RPA3, RRM1, CDK4, CDC2, PLK1, PSMA4, and TK1) .
Expression profiles associated with metastatic MIN-Os
8w-D and 4w-11 are the two MIN-O lines with highest metastatic index (Maglione 2004). From our GO analysis (Table 7), MIN-Os from both lines seem to have dysregulation of significant number of genes that are involved in immune-related responses and apoptotic program. The MIN-Os from these two lines share 209 probes that are dysregulated [see Additional file 4]. Among the 209 probes (186 genes and 18 ESTs), 42 genes are unique to these two lines. Examples of these unique genes include Serpine2, which is overexpressed in metastatic cancer; Ctnnb1, Fliih, Pdlim4/Ril, all of which affect migratory behavior of cells; and Gas6, which is a ligand for the Axl oncogene .
The metastasis signature genes in human cancer  have been reported to be predictive of a high metastatic rate in mouse tumors . 14 out of these 17 signature genes are represented in the MG-U74Av2 array. We compared the expression of these 14 genes in both MIN-Os and tumors, between the two most metastatic MIN-O lines (8w-D and 4w-11) and the two least metastatic MIN-O lines (8w-B and 4w-4). We did not find significant correlation of the expression trends between what has reported with human metastatic breast cancer and our metastatic MIN-O lines (data not shown), yet, further analysis and comparison to the human data is required in order to identify potential correlated GO terms and pathway hubs that may better reflect similarity in the model and the human data.
Expression profiles associated with estrogen sensitivity
8w-B and 8w-D MIN-Os are more sensitive to ovarian ablation than other lines. Transcript levels of ERs in the 8w-B and 8w-D MIN-Os tend to be slightly higher than the 4w-4 and 4w-11 MIN-Os [see Additional file 5], which is consistent with ER expression found in the immunohistochemical analysis. Moreover, the expression of cyp19a1, an aromatase which catalyzes essential reactions for estrogen biosynthesis, is also slightly but statistically significantly higher in the 8w-B and 8w-D lines (fold change = 1.07, p = 0.0068). In these lines, steroid biosynthesis (p = 0.0345 for 8w-B, p = 3 × 10-6 for 8w-D) and sterol metabolism (p = 0.0416 for 8w-D) were represented as highly significant GO functions (Table 7). In addition, angiogenesis related functions were also significant. 271 probes (114 probes unique to only these two lines) were differentially expressed in both 8w-B and 8w-D MIN-Os [see Additional file 6]. Genes dysregulated only in 8w-B and 8w-D MIN-Os include proteins that interact with estrogen or ER, such as Abcg1, Gnai2, Fbln1, Nr4a1, Jun (all up-regulated), and Strn and Cyp1b1 (both down-regulated). Also, NfκB pathway genes – Nfkb2, Relb, and Bcl10 – were down-regulated.
Major chromosome abnormalities associated with tumorigenesis are present at the MIN-O stage
This study is the first molecular analysis of mammary tumor progression in the four MIN-O mouse lines, and it shows that MIN-Os, mouse mammary premalignant lesions, are already genetically and molecularly advanced, with molecular and genomic profiles similar to those of invasive tumors. This study also confirms that these lesions, which originated in the same organ site, are genetically heterogeneous, as reflected from their morphological, genomic, and gene expression characteristics. We hypothesize that even with the same initiating oncogene, the resulting early premalignant lesions and subsequent tumors will have different groups of biologically distinct characteristics that are determined by the genetic changes already present in the early lesions – in this case, in the MIN-Os.
From the whole genome analysis of MIN-O and tumor samples by the array CGH, we found that most lesions already have chromosome aberrations at the MIN-O stage. The results from the expression arrays mirror the CGH results in that major expression changes were seen in the MIN-O stage and expression changes between the MIN-Os and tumors are surprisingly relatively few. The finding that the MIN-Os already have the major genomic and expression changes that have previously been associated with tumorigenesis is consistent with observations in human breast cancer [3, 4]. Human DCIS has been found to have the same chromosome alterations found in invasive ductal carcinoma . In addition, at the gene expression level, ADH, DCIS, and IDC from the same individual are highly similar to each other [4, 37].
Our findings show whole chromosome gains in the MIN-Os and tumors to be fairly common, while no whole chromosome loss was observed. This is consistent with LOH analysis that did not find increased chromosome loss or large deletions in Tg(PyVmT) mammary tumor . A previous CGH study by Hodgson et al. found that chromosomes 2, 11 and 15 most often had recurrent gains in PyVmT mammary tumors . Indeed, the tumor samples we analyzed frequently had whole chromosome gains of these chromosomes, as well as chromosomes 1 and 10. However, there were no genetic changes found consistently during the transition from MIN-O to tumor in the different lines. Moreover, as expected from the CGH analysis, there seem to be no common expression changes that are associated with every MIN-O to tumor transition event.
The whole genome analysis with aCGH also led to the identification of areas of small deletions. Moreover, using a custom-designed array, the aberrations were confirmed and the breakpoints were refined. The smallest area confirmed by the zoom-in array was the 59,704 bp deletion on chromosome 3, which spans only one known gene, 2810410M20Rik. The zoom-in array clearly demonstrated the loss of 56 probes that represented this stretch of chromosome 3 in MIN-Os and the corresponding tumor samples. Moreover, for the other two areas of deletion, the losses of DNA were confirmed to span the same areas in the MIN-Os and the tumors. In general, the same small aberrations as well as large gains were seen in both MIN-Os and tumors, suggesting that these genomic changes were already present in the MIN-Os.
Chromosome aberrations in human breast cancer has been studied with LOH analysis  and array CGH [2, 3, 5, 40–43]. These studies have identified regions of frequent copy number changes in breast cancer. For example, ErbB2 amplicon on 17q21 is amplified in 20–30% of breast cancers . Frequent copy number changes of certain chromosome regions are associated with BRCA1 and BRCA2 tumors, and they can be classified by the genomic copy number alterations . More recently, Fridlyand et al. has shown that breast tumors can be classified into three subtypes based on the copy number alteration phenotypes . The "mixed amplifier" group had low level gains and losses with some recurring copy number changes, while the others two types had either few copy number changes ("1q/16q" group) or significant changes encompassing nearly 60% of the genome with copy number alterations ("complex" group). The patterns of genome alteration in our MIN-O lines resemble those of the first two groups. Interestingly, the majority of human tumors in these two groups were ER positive and that these two groups had better outcome compared to the complex group . Our MIN-O lines may better represent human breast cancers that are low grade and less aggressive than those that are ER negative, high grade lesions with poor outcome.
Major expression changes associated with tumorigenesis
Although there were no common genomic alterations found among all MIN-O samples examined, there seem to be common gene expression changes associated with the transition from normal mammary gland to MIN-O. For example, MIN-Os over-expressed cyclin D1 and Cdkn1a (P21). Cyclin D1 is a well-known oncogene, which is commonly over-expressed in human breast cancer . Moreover, recent reports demonstrated that cyclin D1-associated kinase activity, through its association with CDK4, is required for ErbB2-dependent mammary tumorigenesis [47–49]. Cdkn1a (P21) was also significantly over-expressed in the MIN-Os. Cdkn1a is a cyclin-dependent kinase inhibitor and is commonly known to induce G1 arrest via the p53 pathway . However this is not incongruous with Cdkn1a over-expression in the MIN-Os, as p21 is often over-expressed in human cancer, and more studies have shown p21 as a modulator of cell proliferation, rather than inhibitor [51, 52]. Dysregulation of these genes in the MIN-Os demonstrates that although these tissues are phenotypically premalignant, the MIN-Os already have molecular changes known to be associated with tumorigenesis.
Since major molecular changes had occurred in the MIN-Os, tumors developed from the MIN-Os had surprisingly small genetic differences from the corresponding MIN-O tissues. Many genes coding for extracellular proteins were down-regulated in tumors. Among these genes, Lum and Gsn are two genes whose association with breast cancer has been well studied. Decreased expression of Lum is associated with poor outcome of invasive breast cancer  and down-regulation of Gsn is associated with breast cancer progression .
Changes associated with aneuploidy
Each MIN-O line, established from a distinct MIN lesion in the original transgenic mammary fat pad, has distinct genomic and expression profiles and biological characteristics. This is also similar to human DCIS lesions, because DCIS is subclassified into different histological grades and the pattern of chromosome alterations and gene expression signatures differs between different histological grades [3, 4]. MIN-Os can be classified into subgroups by genomic and molecular profiles that reflect biological behavior (metastatic potential, tumor latency, estrogen sensitivity).
Genetically unstable MIN-O lines, which have multiple whole chromosome gains, almost always had gain of chromosome 2. AurkA was one of the genes on chromosome 2 that was highly over-expressed in these MIN-Os. AurkA is localized to the spindles and its overexpression induces aberrant mitosis that results in centrosome duplication and aneuploidy . Overexpression of AurkA is associated with cancer, including human DCIS [56, 57]. Other genes involved in mitotic spindle assembly and centrosome duplication were also over-expressed in these MIN-Os. Ran, Kif22/Kid, Plk1 and Nde1  are important for mitotic spindle assembly [58–61]. Mitotic spindle components, tublins alpha1 and beta6 were both over-expressed. Plk1 also regulates centrosome assembly and separation, and controls the onset of G2/M transition in the cell cycle by activating cdc25c . Interestingly, beta-tubulin and cdc25c are substrates of Plk1 . Other substrates of Plk1 are subunits of anaphase promoting complex (APC), which regulates progression and exit through mitosis by proteolysis of cell cycle regulators . Two APC related genes, a component of APC, Apc11, and an activator of APC, cdc20, were both up-regulated in these MIN-Os. APC/cdc20 complex is important for metaphase to anaphase transition . These genes highlight dysregulation of G2/M phase, especially during metaphase/anaphase when sister chromatids and spindle poles separate.
Changes associated with metastasis
Highly metastatic MIN-O lines had alteration of genes involved in immune response and apoptotic processes. Recently, studies have shown that the inflammation process is closely linked to cancer development [63, 64]. Innate immune systems promote tumorigenesis by producing factors that activate tumor associated cells and changes stromal microenvironment . It is believed that changes in apoptotic tumor cell death pattern affect the tumor microenvironment and the associated immune response . In particular, for the metastatic processes, tumor cells utilize many of the same factors, such as adhesion molecules, chemokines and receptors, used by inflammatory cells . In the Tg(PyVmT) model, lack of functional macrophages inhibits the development of invasive mammary carcinoma and pulmonary metastasis , and other factors that regulate tumor microenvironment also affect the metastatic potential [69–72]. Gas6 is the ligand for Axl receptor tyrosine kinase, which regulates immune system and inhibits apoptosis [73, 74]. Bcl6 is a transcriptional regulator that has an anti-apoptotic role and also affects humoral immune response . Both of these genes were highly over-expressed in both metastatic MIN-O lines. Moreover, the expression of ICAM-1, an endothelial intercellular adhesion molecule was down-regulated in these MIN-Os. ICAM-1 is known to be down-regulated by tumor derived angiogenic factors in order to inhibit leukocyte infiltration and to escape from anti-tumor immune surveillance [76, 77]. These expression changes suggest the alteration of the genes related to the immune modulation in these MIN-Os, which may contribute to the development of tumor associated environment that promote migration and invasion.
The "Metastatic signature gene set"  was not differentially expressed between the metastatic and non-metastatic MIN-O lines. This seems to be contradictory to the results by Qiu et al  in which the signature gene set was predictive of metastatic behaviour of Tg(PyVmT) tumors. However, this study compares MIN-Os with different metastatic potentials, which express PyVmT as the initiating oncogene in the same genetic background (FVB), and is fundamentally different from the study by Qiu et al in which Tg(PyVmT) tumors from various genetic backgrounds were compared. Thus, while the metastatic signature gene set may be effective in classifying metastatic disease from a group of genetically polymorphic samples, the MIN-O mouse lines will be useful for investigating specific molecular mechanisms responsible for metastasis in the disease tissues.
Changes associated with estrogen dependence
ER status of the MIN-Os, in general, reflected the ovarian hormone sensitivity of the lesions. The 8w-B and 8w-D lines, which were associated with nuclear ER-positive MIN-O cells, were the most sensitive to the ovariectomy. In all four MIN-O lines studied, significantly less ER-positive cells were found in the tumors from the intact hosts, and in the MIN-Os from the ovariectomized hosts.
Two common functions that were affected in the estrogen sensitive MIN-O lines were steroid biosynthesis and angiogenesis. Significance of these two functions as well as the expression differences in the ERs and aromatase may represent the difference in estrogen mediated signaling in the 8w-B and D MIN-Os from the rest of the MIN-O lines, since estrogen modulates endothelial cell behavior, including angiogensis . This concept is further supported by the dysregulation of genes that are known to be regulated by estrogen or ER, such as Cyp1b1 , Nf-kB2 , RelB, Strn , Fbln1  or genes that interact with estrogen or ER, such as Abcg1, Gnai2, and Jun . The down-regulation of two Nf-kB family members, Nf-kB2 and RelB, maybe significant since ER alpha is known to inhibit the Nf-kB activity by various mechanisms, including direct binding to the Nf-kB complex to inhibit transcriptional activation . Jun was up-regulated in estrogen sensitive MIN-Os. ER activates transcription at AP-1 sites via interaction with Jun/Fos/coactivator complex . Dysregulation of these genes seems to draw a picture of increase ER and estrogen activities in these MIN-Os.
In summary, in this mouse model of breast cancer, the majority of the molecular changes responsible for tumorigenesis are already present in early premalignant lesions. The progression to invasive carcinoma is associated with only few additional changes, and these changes are not common in all tumors. Moreover, the molecular heterogeneity of the premalignant lesions represents phenotypic differences of these lesions. Thus, the MIN-O mouse model recapitulates both the preencoding and heterogeneity found in early human breast lesions and may be used as a high fidelity model for the biology of human breast cancer.
List of abbreviations used
- comparative genomic hybridization (CGH):
ductal carcinoma in situ (DCIS), estrogen receptor (ER), invasive ductal carcinoma (IDC), mammary intraepithelial neoplasia outgrowth (MIN-O), MIN (mammary intraepithelial neoplasia), polyomavirus middle-T (PyVmT), array comparative genomic hybridization (aCGH), Gene Ontology (GO)
This work was supported by Grants 11IB-0123 (JPG) and 9FB-0212 (RN) from the California Breast Cancer Research Program and the Investigator Initiated Grant by Department of Pathology at UC Davis. Processing and analysis of gene expression was done in Gene Expression Shared Resource (2 P30 CA93373-04, P30-ES05707). Processing and immunohistochemistry of fixed tissue was done in the mutant mouse pathology lab (CCM, UC Davis).
- Fearon ER, Vogelstein B: A genetic model for colorectal tumorigenesis. Cell. 1990, 61 (5): 759-767. 10.1016/0092-8674(90)90186-I.View ArticlePubMedGoogle Scholar
- Chin K, de Solorzano CO, Knowles D, Jones A, Chou W, Rodriguez EG, Kuo WL, Ljung BM, Chew K, Myambo K, Miranda M, Krig S, Garbe J, Stampfer M, Yaswen P, Gray JW, Lockett SJ: In situ analyses of genome instability in breast cancer. Nat Genet. 2004, 36 (9): 984-988. 10.1038/ng1409.View ArticlePubMedGoogle Scholar
- Hwang ES, DeVries S, Chew KL, Moore DH, Kerlikowske K, Thor A, Ljung BM, Waldman FM: Patterns of chromosomal alterations in breast ductal carcinoma in situ. Clin Cancer Res. 2004, 10 (15): 5160-5167. 10.1158/1078-0432.CCR-04-0165.View ArticlePubMedGoogle Scholar
- Ma XJ, Salunga R, Tuggle JT, Gaudet J, Enright E, McQuary P, Payette T, Pistone M, Stecker K, Zhang BM, Zhou YX, Varnholt H, Smith B, Gadd M, Chatfield E, Kessler J, Baer TM, Erlander MG, Sgroi DC: Gene expression profiles of human breast cancer progression. Proc Natl Acad Sci U S A. 2003, 100 (10): 5974-5979. 10.1073/pnas.0931261100.View ArticlePubMedPubMed CentralGoogle Scholar
- Aubele M, Mattis A, Zitzelsberger H, Walch A, Kremer M, Welzl G, Hofler H, Werner M: Extensive ductal carcinoma In situ with small foci of invasive ductal carcinoma: evidence of genetic resemblance by CGH. Int J Cancer. 2000, 85 (1): 82-86. 10.1002/(SICI)1097-0215(20000101)85:1<82::AID-IJC15>3.0.CO;2-S.View ArticlePubMedGoogle Scholar
- Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, Fluge O, Pergamenschikov A, Williams C, Zhu SX, Lonning PE, Borresen-Dale AL, Brown PO, Botstein D: Molecular portraits of human breast tumours. Nature. 2000, 406 (6797): 747-752. 10.1038/35021093.View ArticlePubMedGoogle Scholar
- Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn M, Jeffrey SS, Thorsen T, Quist H, Matese JC, Brown PO, Botstein D, Eystein Lonning P, Borresen-Dale AL: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci U S A. 2001, 98 (19): 10869-10874. 10.1073/pnas.191367098.View ArticlePubMedPubMed CentralGoogle Scholar
- van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415 (6871): 530-536. 10.1038/415530a.View ArticlePubMedGoogle Scholar
- Mallon E, Osin P, Nasiri N, Blain I, Howard B, Gusterson B: The basic pathology of human breast cancer. J Mammary Gland Biol Neoplasia. 2000, 5 (2): 139-163. 10.1023/A:1026439204849.View ArticlePubMedGoogle Scholar
- Shoker BS, Sloane JP: DCIS grading schemes and clinical implications. Histopathology. 1999, 35 (5): 393-400. 10.1046/j.1365-2559.1999.035005393.x.View ArticlePubMedGoogle Scholar
- Maglione JE, McGoldrick ET, Young LJ, Namba R, Gregg JP, Liu L, Moghanaki D, Ellies LG, Borowsky AD, Cardiff RD, MacLeod CL: Polyomavirus middle T-induced mammary intraepithelial neoplasia outgrowths: single origin, divergent evolution, and multiple outcomes. Mol Cancer Ther. 2004, 3 (8): 941-953.PubMedGoogle Scholar
- Namba R, Maglione JE, Young LJ, Borowsky AD, Cardiff RD, MacLeod CL, Gregg JP: Molecular characterization of the transition to malignancy in a genetically engineered mouse-based model of ductal carcinoma in situ. Mol Cancer Res. 2004, 2 (8): 453-463.PubMedGoogle Scholar
- Cardiff RD, Anver MR, Gusterson BA, Hennighausen L, Jensen RA, Merino MJ, Rehm S, Russo J, Tavassoli FA, Wakefield LM, Ward JM, Green JE: The mammary pathology of genetically engineered mice: the consensus report and recommendations from the Annapolis meeting. Oncogene. 2000, 19 (8): 968-988. 10.1038/sj.onc.1203277.View ArticlePubMedGoogle Scholar
- Maglione JE, Moghanaki D, Young LJ, Manner CK, Ellies LG, Joseph SO, Nicholson B, Cardiff RD, MacLeod CL: Transgenic Polyoma middle-T mice model premalignant mammary disease. Cancer Res. 2001, 61 (22): 8298-8305.PubMedGoogle Scholar
- Namba R, Young LJ, Maglione JE, McGoldrick ET, Liu S, Wurz GT, DeGregorio MW, Borowsky AD, MacLeod CL, Cardiff RD, Gregg JP: Selective estrogen receptor modulators inhibit growth and progression of premalignant lesions in a mouse model of DCIS. Breast Cancer Res. 2005, 7 (6): R881-889. 10.1186/bcr1317.View ArticlePubMedPubMed CentralGoogle Scholar
- RepeatMasker. [http://www.repeatmasker.org]
- Hughes TR, Mao M, Jones AR, Burchard J, Marton MJ, Shannon KW, Lefkowitz SM, Ziman M, Schelter JM, Meyer MR, Kobayashi S, Davis C, Dai H, He YD, Stephaniants SB, Cavet G, Walker WL, West A, Coffey E, Shoemaker DD, Stoughton R, Blanchard AP, Friend SH, Linsley PS: Expression profiling using microarrays fabricated by an ink-jet oligonucleotide synthesizer. Nat Biotechnol. 2001, 19 (4): 342-347. 10.1038/86730.View ArticlePubMedGoogle Scholar
- Agilent Technologies. [http://www.chem.agilent.com]
- Barrett MT, Scheffer A, Ben-Dor A, Sampas N, Lipson D, Kincaid R, Tsang P, Curry B, Baird K, Meltzer PS, Yakhini Z, Bruhn L, Laderman S: Comparative genomic hybridization using oligonucleotide microarrays and total genomic DNA. Proc Natl Acad Sci U S A. 2004, 101 (51): 17765-17770. 10.1073/pnas.0407979101.View ArticlePubMedPubMed CentralGoogle Scholar
- ErmineJ software. [http://microarray.cu-genome.org/ermineJ/]
- Lee HK, Braynen W, Keshav K, Pavlidis P: ErmineJ: tool for functional analysis of gene expressiono data sets. BMC Bioinformatics. 2005, 6: 269-10.1186/1471-2105-6-269.View ArticlePubMedPubMed CentralGoogle Scholar
- UCSC Genome Browser . [http://genome.ucsc.edu/cgi-bin/hgGateway]
- Lin EY, Jones JG, Li P, Zhu L, Whitney KD, Muller WJ, Pollard JW: Progression to malignancy in the polyoma middle T oncoprotein mouse breast cancer model provides a reliable model for human diseases. Am J Pathol. 2003, 163 (5): 2113-2126.View ArticlePubMedPubMed CentralGoogle Scholar
- Hodgson JG, Malek T, Bornstein S, Hariono S, Ginzinger DG, Muller WJ, Gray JW: Copy number aberrations in mouse breast tumors reveal loci and genes important in tumorigenic receptor tyrosine kinase signaling. Cancer Res. 2005, 65 (21): 9695-9704. 10.1158/0008-5472.CAN-05-0755.View ArticlePubMedGoogle Scholar
- Anderson DM, Maraskovsky E, Billingsley WL, Dougall WC, Tometsko ME, Roux ER, Teepe MC, DuBose RF, Cosman D, Galibert L: A homologue of the TNF receptor and its ligand enhance T-cell growth and dendritic-cell function. Nature. 1997, 390 (6656): 175-179. 10.1038/36593.View ArticlePubMedGoogle Scholar
- Lim KH, Baines AT, Fiordalisi JJ, Shipitsin M, Feig LA, Cox AD, Der CJ, Counter CM: Activation of RalA is critical for Ras-induced tumorigenesis of human cells. Cancer Cell. 2005, 7 (6): 533-545. 10.1016/j.ccr.2005.04.030.View ArticlePubMedGoogle Scholar
- Gery S, Tanosaki S, Bose S, Bose N, Vadgama J, Koeffler HP: Down-regulation and growth inhibitory role of C/EBPalpha in breast cancer. Clin Cancer Res. 2005, 11 (9): 3184-3190. 10.1158/1078-0432.CCR-04-2625.View ArticlePubMedGoogle Scholar
- Grimm SL, Rosen JM: The role of C/EBPbeta in mammary gland development and breast cancer. J Mammary Gland Biol Neoplasia. 2003, 8 (2): 191-204. 10.1023/A:1025900908026.View ArticlePubMedGoogle Scholar
- Shi B, Vinyals A, Alia P, Broceno C, Chen F, Adrover M, Gelpi C, Price JE, Fabra A: Differential expression of MHC class II molecules in highly metastatic breast cancer cells is mediated by the regulation of the CIITA transcription Implication of CIITA in tumor and metastasis development. Int J Biochem Cell Biol. 2006, 38 (4): 544-562. 10.1016/j.biocel.2005.07.012.View ArticlePubMedGoogle Scholar
- Tamiolakisl D, Venizelos I, Lambropoulou M, Jivannakis T, Seliniotakis E, Tsikouras P, Limberis V, Tsalkidis A, Papadopoulos N: Gains and losses of HLA class II (DR) and CD4 in atypical hyperplasia, carcinoma in situ and infiltrating ductal carcinoma of the breast. Acta Medica (Hradec Kralove). 2004, 47 (4): 257-262.Google Scholar
- Eckerdt F, Yuan J, Strebhardt K: Polo-like kinases and oncogenesis. Oncogene. 2005, 24 (2): 267-276. 10.1038/sj.onc.1208273.View ArticlePubMedGoogle Scholar
- Kinch MS, Carles-Kinch K: Overexpression and functional alterations of the EphA2 tyrosine kinase in cancer. Clin Exp Metastasis. 2003, 20 (1): 59-68. 10.1023/A:1022546620495.View ArticlePubMedGoogle Scholar
- Ren B, Cam H, Takahashi Y, Volkert T, Terragni J, Young RA, Dynlacht BD: E2F integrates cell cycle progression with DNA repair, replication, and G(2)/M checkpoints. Genes Dev. 2002, 16 (2): 245-256. 10.1101/gad.949802.View ArticlePubMedPubMed CentralGoogle Scholar
- Holland SJ, Powell MJ, Franci C, Chan EW, Friera AM, Atchison RE, McLaughlin J, Swift SE, Pali ES, Yam G, Wong S, Lasaga J, Shen MR, Yu S, Xu W, Hitoshi Y, Bogenberger J, Nor JE, Payan DG, Lorens JB: Multiple roles for the receptor tyrosine kinase axl in tumor formation. Cancer Res. 2005, 65 (20): 9294-9303. 10.1158/0008-5472.CAN-05-0993.View ArticlePubMedGoogle Scholar
- Ramaswamy S, Ross KN, Lander ES, Golub TR: A molecular signature of metastasis in primary solid tumors. Nat Genet. 2003, 33 (1): 49-54. 10.1038/ng1060.View ArticlePubMedGoogle Scholar
- Qiu TH, Chandramouli GV, Hunter KW, Alkharouf NW, Green JE, Liu ET: Global expression profiling identifies signatures of tumor virulence in MMTV-PyMT-transgenic mice: correlation to human disease. Cancer Res. 2004, 64 (17): 5973-5981. 10.1158/0008-5472.CAN-04-0242.View ArticlePubMedGoogle Scholar
- Schuetz CS, Bonin M, Clare SE, Nieselt K, Sotlar K, Walter M, Fehm T, Solomayer E, Riess O, Wallwiener D, Kurek R, Neubauer HJ: Progression-specific genes identified by expression profiling of matched ductal carcinomas in situ and invasive breast tumors, combining laser capture microdissection and oligonucleotide microarray analysis. Cancer Res. 2006, 66 (10): 5278-5286. 10.1158/0008-5472.CAN-05-4610.View ArticlePubMedGoogle Scholar
- Ritland SR, Rowse GJ, Chang Y, Gendler SJ: Loss of heterozygosity analysis in primary mammary tumors and lung metastases of MMTV-MTAg and MMTV-neu transgenic mice. Cancer Res. 1997, 57 (16): 3520-3525.PubMedGoogle Scholar
- Wang ZC, Lin M, Wei LJ, Li C, Miron A, Lodeiro G, Harris L, Ramaswamy S, Tanenbaum DM, Meyerson M, Iglehart JD, Richardson A: Loss of heterozygosity and its correlation with expression profiles in subclasses of invasive breast cancers. Cancer Res. 2004, 64 (1): 64-71. 10.1158/0008-5472.CAN-03-2570.View ArticlePubMedGoogle Scholar
- Jonsson G, Naylor TL, Vallon-Christersson J, Staaf J, Huang J, Ward MR, Greshock JD, Luts L, Olsson H, Rahman N, Stratton M, Ringner M, Borg A, Weber BL: Distinct genomic profiles in hereditary breast tumors identified by array-based comparative genomic hybridization. Cancer Res. 2005, 65 (17): 7612-7621.PubMedGoogle Scholar
- Gunther K, Merkelbach-Bruse S, Amo-Takyi BK, Handt S, Schroder W, Tietze L: Differences in genetic alterations between primary lobular and ductal breast cancers detected by comparative genomic hybridization. J Pathol. 2001, 193 (1): 40-47. 10.1002/1096-9896(2000)9999:9999<::AID-PATH745>3.0.CO;2-N.View ArticlePubMedGoogle Scholar
- Yao J, Weremowicz S, Feng B, Gentleman RC, Marks JR, Gelman R, Brennan C, Polyak K: Combined cDNA array comparative genomic hybridization and serial analysis of gene expression analysis of breast tumor progression. Cancer Res. 2006, 66 (8): 4065-4078. 10.1158/0008-5472.CAN-05-4083.View ArticlePubMedGoogle Scholar
- O'Connell P: Genetic and cytogenetic analyses of breast cancer yield different perspectives of a complex disease. Breast Cancer Res Treat. 2003, 78 (3): 347-357. 10.1023/A:1023037925950.View ArticlePubMedGoogle Scholar
- Slamon DJ, Godolphin W, Jones LA, Holt JA, Wong SG, Keith DE, Levin WJ, Stuart SG, Udove J, Ullrich A, et al: Studies of the HER-2/neu proto-oncogene in human breast and ovarian cancer. Science. 1989, 244 (4905): 707-712. 10.1126/science.2470152.View ArticlePubMedGoogle Scholar
- Fridlyand J, Snijders AM, Ylstra B, Li H, Olshen A, Segraves R, Dairkee S, Tokuyasu T, Ljung BM, Jain AN, McLennan J, Ziegler J, Chin K, Devries S, Feiler H, Gray JW, Waldman F, Pinkel D, Albertson DG: Breast tumor copy number aberration phenotypes and genomic instability. BMC Cancer. 2006, 6: 96-10.1186/1471-2407-6-96.View ArticlePubMedPubMed CentralGoogle Scholar
- Sutherland RL, Musgrove EA: Cyclin D1 and mammary carcinoma: new insights from transgenic mouse models. Breast Cancer Res. 2002, 4 (1): 14-17. 10.1186/bcr411.View ArticlePubMedGoogle Scholar
- Landis MW, Pawlyk BS, Li T, Sicinski P, Hinds PW: Cyclin D1-dependent kinase activity in murine development and mammary tumorigenesis. Cancer Cell. 2006, 9 (1): 13-22. 10.1016/j.ccr.2005.12.019.View ArticlePubMedGoogle Scholar
- Reddy HK, Mettus RV, Rane SG, Grana X, Litvin J, Reddy EP: Cyclin-dependent kinase 4 expression is essential for neu-induced breast tumorigenesis. Cancer Res. 2005, 65 (22): 10174-10178. 10.1158/0008-5472.CAN-05-2639.View ArticlePubMedGoogle Scholar
- Yu Q, Sicinska E, Geng Y, Ahnstrom M, Zagozdzon A, Kong Y, Gardner H, Kiyokawa H, Harris LN, Stal O, Sicinski P: Requirement for CDK4 kinase function in breast cancer. Cancer Cell. 2006, 9 (1): 23-32. 10.1016/j.ccr.2005.12.012.View ArticlePubMedGoogle Scholar
- Liu G, Lozano G: p21 stability: linking chaperones to a cell cycle checkpoint. Cancer Cell. 2005, 7 (2): 113-114. 10.1016/j.ccr.2005.01.019.View ArticlePubMedGoogle Scholar
- Weiss RH: p21Waf1/Cip1 as a therapeutic target in breast and other cancers. Cancer Cell. 2003, 4 (6): 425-429. 10.1016/S1535-6108(03)00308-8.View ArticlePubMedGoogle Scholar
- Yuste L, Montero JC, Esparis-Ogando A, Pandiella A: Activation of ErbB2 by overexpression or by transmembrane neuregulin results in differential signaling and sensitivity to herceptin. Cancer Res. 2005, 65 (15): 6801-6810. 10.1158/0008-5472.CAN-04-4023.View ArticlePubMedGoogle Scholar
- Troup S, Njue C, Kliewer EV, Parisien M, Roskelley C, Chakravarti S, Roughley PJ, Murphy LC, Watson PH: Reduced expression of the small leucine-rich proteoglycans, lumican, and decorin is associated with poor outcome in node-negative invasive breast cancer. Clin Cancer Res. 2003, 9 (1): 207-214.PubMedGoogle Scholar
- Winston JS, Asch HL, Zhang PJ, Edge SB, Hyland A, Asch BB: Downregulation of gelsolin correlates with the progression to breast carcinoma. Breast Cancer Res Treat. 2001, 65 (1): 11-21. 10.1023/A:1006446108411.View ArticlePubMedGoogle Scholar
- Meraldi P, Honda R, Nigg EA: Aurora-A overexpression reveals tetraploidization as a major route to centrosome amplification in p53-/- cells. Embo J. 2002, 21 (4): 483-492. 10.1093/emboj/21.4.483.View ArticlePubMedPubMed CentralGoogle Scholar
- Hoque A, Carter J, Xia W, Hung MC, Sahin AA, Sen S, Lippman SM: Loss of aurora A/STK15/BTAK overexpression correlates with transition of in situ to invasive ductal carcinoma of the breast. Cancer Epidemiol Biomarkers Prev. 2003, 12 (12): 1518-1522.PubMedGoogle Scholar
- Miyoshi Y, Iwao K, Egawa C, Noguchi S: Association of centrosomal kinase STK15/BTAK mRNA expression with chromosomal instability in human breast cancers. Int J Cancer. 2001, 92 (3): 370-373. 10.1002/ijc.1200.View ArticlePubMedGoogle Scholar
- Feng Y, Walsh CA: Mitotic spindle regulation by Nde1 controls cerebral cortical size. Neuron. 2004, 44 (2): 279-293. 10.1016/j.neuron.2004.09.023.View ArticlePubMedGoogle Scholar
- Castro A, Bernis C, Vigneron S, Labbe JC, Lorca T: The anaphase-promoting complex: a key factor in the regulation of cell cycle. Oncogene. 2005, 24 (3): 314-325. 10.1038/sj.onc.1207973.View ArticlePubMedGoogle Scholar
- Peters JM: The anaphase-promoting complex: proteolysis in mitosis and beyond. Mol Cell. 2002, 9 (5): 931-943. 10.1016/S1097-2765(02)00540-3.View ArticlePubMedGoogle Scholar
- van Vugt MA, Medema RH: Getting in and out of mitosis with Polo-like kinase-1. Oncogene. 2005, 24 (17): 2844-2859. 10.1038/sj.onc.1208617.View ArticlePubMedGoogle Scholar
- Myer DL, Bahassi el M, Stambrook PJ: The Plk3-Cdc25 circuit. Oncogene. 2005, 24 (2): 299-305. 10.1038/sj.onc.1208278.View ArticlePubMedGoogle Scholar
- Coussens LM, Werb Z: Inflammation and cancer. Nature. 2002, 420 (6917): 860-867. 10.1038/nature01322.View ArticlePubMedPubMed CentralGoogle Scholar
- Philip M, Rowley DA, Schreiber H: Inflammation as a tumor promoter in cancer induction. Semin Cancer Biol. 2004, 14 (6): 433-439. 10.1016/j.semcancer.2004.06.006.View ArticlePubMedGoogle Scholar
- Bissell MJ, Radisky D: Putting tumours in context. Nat Rev Cancer. 2001, 1 (1): 46-54. 10.1038/35094059.View ArticlePubMedPubMed CentralGoogle Scholar
- Zeh HJ, Lotze MT: Addicted to death: invasive cancer and the immune response to unscheduled cell death. J Immunother. 2005, 28 (1): 1-9. 10.1097/00002371-200501000-00001.View ArticlePubMedGoogle Scholar
- Muller A, Homey B, Soto H, Ge N, Catron D, Buchanan ME, McClanahan T, Murphy E, Yuan W, Wagner SN, Barrera JL, Mohar A, Verastegui E, Zlotnik A: Involvement of chemokine receptors in breast cancer metastasis. Nature. 2001, 410 (6824): 50-56. 10.1038/35065016.View ArticlePubMedGoogle Scholar
- Lin EY, Pollard JW: Macrophages: modulators of breast cancer progression. Novartis Found Symp. 2004, 256: 158-68; discussion 168-72, 259-69.View ArticlePubMedGoogle Scholar
- Forrester E, Chytil A, Bierie B, Aakre M, Gorska AE, Sharif-Afshar AR, Muller WJ, Moses HL: Effect of conditional knockout of the type II TGF-beta receptor gene in mammary epithelia on mammary gland development and polyomavirus middle T antigen induced tumor formation and metastasis. Cancer Res. 2005, 65 (6): 2296-2302. 10.1158/0008-5472.CAN-04-3272.View ArticlePubMedGoogle Scholar
- Muraoka-Cook RS, Kurokawa H, Koh Y, Forbes JT, Roebuck LR, Barcellos-Hoff MH, Moody SE, Chodosh LA, Arteaga CL: Conditional overexpression of active transforming growth factor beta1 in vivo accelerates metastases of transgenic mammary tumors. Cancer Res. 2004, 64 (24): 9002-9011. 10.1158/0008-5472.CAN-04-2111.View ArticlePubMedGoogle Scholar
- Muraoka RS, Dumont N, Ritter CA, Dugger TC, Brantley DM, Chen J, Easterly E, Roebuck LR, Ryan S, Gotwals PJ, Koteliansky V, Arteaga CL: Blockade of TGF-beta inhibits mammary tumor cell viability, migration, and metastases. J Clin Invest. 2002, 109 (12): 1551-1559. 10.1172/JCI200215234.View ArticlePubMedPubMed CentralGoogle Scholar
- Bugge TH, Lund LR, Kombrinck KK, Nielsen BS, Holmback K, Drew AF, Flick MJ, Witte DP, Dano K, Degen JL: Reduced metastasis of Polyoma virus middle T antigen-induced mammary cancer in plasminogen-deficient mice. Oncogene. 1998, 16 (24): 3097-3104. 10.1038/sj.onc.1201869.View ArticlePubMedGoogle Scholar
- Lee WP, Wen Y, Varnum B, Hung MC: Akt is required for Axl-Gas6 signaling to protect cells from E1A-mediated apoptosis. Oncogene. 2002, 21 (3): 329-336. 10.1038/sj.onc.1205066.View ArticlePubMedGoogle Scholar
- Schwartzberg PL: Immunology. Tampering with the immune system. Science. 2001, 293 (5528): 228-229. 10.1126/science.1063291.View ArticlePubMedGoogle Scholar
- Phan RT, Dalla-Favera R: The BCL6 proto-oncogene suppresses p53 expression in germinal-centre B cells. Nature. 2004, 432 (7017): 635-639. 10.1038/nature03147.View ArticlePubMedGoogle Scholar
- Bouma-ter Steege JC, Baeten CI, Thijssen VL, Satijn SA, Verhoeven IC, Hillen HF, Wagstaff J, Griffioen AW: Angiogenic profile of breast carcinoma determines leukocyte infiltration. Clin Cancer Res. 2004, 10 (21): 7171-7178. 10.1158/1078-0432.CCR-04-0742.View ArticlePubMedGoogle Scholar
- Griffioen AW, Damen CA, Martinotti S, Blijham GH, Groenewegen G: Endothelial intercellular adhesion molecule-1 expression is suppressed in human malignancies: the role of angiogenic factors. Cancer Res. 1996, 56 (5): 1111-1117.PubMedGoogle Scholar
- Cid MC, Schnaper HW, Kleinman HK: Estrogens and the vascular endothelium. Ann N Y Acad Sci. 2002, 966: 143-157.View ArticlePubMedGoogle Scholar
- Tsuchiya Y, Nakajima M, Yokoi T: Cytochrome P450-mediated metabolism of estrogens and its regulation in human. Cancer Lett. 2005, 227 (2): 115-124. 10.1016/j.canlet.2004.10.007.View ArticlePubMedGoogle Scholar
- Kalaitzidis D, Gilmore TD: Transcription factor cross-talk: the estrogen receptor and NF-kappaB. Trends Endocrinol Metab. 2005, 16 (2): 46-52. 10.1016/j.tem.2005.01.004.View ArticlePubMedGoogle Scholar
- Lu Q, Pallas DC, Surks HK, Baur WE, Mendelsohn ME, Karas RH: Striatin assembles a membrane signaling complex necessary for rapid, nongenomic activation of endothelial NO synthase by estrogen receptor alpha. Proc Natl Acad Sci U S A. 2004, 101 (49): 17126-17131. 10.1073/pnas.0407492101.View ArticlePubMedPubMed CentralGoogle Scholar
- Moll F, Katsaros D, Lazennec G, Hellio N, Roger P, Giacalone PL, Chalbos D, Maudelonde T, Rochefort H, Pujol P: Estrogen induction and overexpression of fibulin-1C mRNA in ovarian cancer cells. Oncogene. 2002, 21 (7): 1097-1107. 10.1038/sj.onc.1205171.View ArticlePubMedGoogle Scholar
- Kushner PJ, Agard DA, Greene GL, Scanlan TS, Shiau AK, Uht RM, Webb P: Estrogen receptor pathways to AP-1. J Steroid Biochem Mol Biol. 2000, 74 (5): 311-317. 10.1016/S0960-0760(00)00108-4.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/6/275/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.