This Article Abstract Full Text (PDF) Supplemental Material Purchase Article View Shopping Cart Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Montel, V. Articles by Tarin, D. Search for Related Content PubMed PubMed Citation Articles by Montel, V. Articles by Tarin, D.

(American Journal of Pathology. 2005;166:1565-1579.)
© 2005 American Society for Investigative Pathology

Expression Profiling of Primary Tumors and Matched Lymphatic and Lung Metastases in a Xenogeneic Breast Cancer Model

From the Rebecca and John Moores Comprehensive Cancer Center and Department of Pathology, University of California, San Diego, La Jolla, California

	Abstract

Top Abstract Materials and Methods Results Discussion References

 
Using a purpose-designed experimental model, we have defined new, statistically significant, differences in gene expression between heavily and weakly metastatic human breast cancer cell populations, in vivo and in vitro. The differences increased under selection pressures designed to increase metastatic proficiency. Conversely, the expression signatures of primary tumors generated by more aggressive variants, and their matched metastases in the lungs and lymph nodes, all tended to converge. However, the few persisting differences among these selectively enriched malignant growths in the breast, lungs, and lymph nodes were highly statistically significant, implying potential mechanistic involvement of the corresponding genes. The evidence that has emerged from the current work indicates that selective enhancement of metastatic proficiency by serial transplantation co-purifies a subliminal gene expression pattern within the tumor cell population. This signature most likely includes genes participating in metastasis pathogenesis, and we document manageable numbers of candidates for this role. The findings also suggest that metastasis to at least two different organs occurs through closely similar genetic mechanisms.

We therefore used oligonucleotide microarrays to perform high throughput screening of global gene expression in tumors and metastases produced by a unique matched pair of human clonal cell lines of opposite metastatic capabilities, which we have derived from the same breast cancer line, MDA-MB-435,⁶ and confirmed to be isogenic by several methods, including chromosomal analysis⁷ and genetic fingerprinting. There are already some reports of high-density microarray profiling of gene expression patterns in metastatic human primary cancers and metastases,^8-12 but no consensus has yet emerged on any groups of genes that are consistently involved. This may be due to the masking effects resulting from comparisons between samples from individuals of different genetic backgrounds.

The work presented below makes progress from previous approaches by using a tightly controlled, well characterized, xenograft model of breast cancer metastasis.^7,13 This investigative system facilitated direct examination of differences between primary tumors and matched metastases in the lungs and lymph nodes from the same animal and thus eliminated the noise from biological variations between different individuals. To our knowledge, this is the first study to systematically investigate gene expression patterns in matched metastases from both of these organs in the same host. In addition, this investigation provides new data on dynamic gene expression patterns in vivo and in vitro of metastasis-competent and incompetent human cell populations within the same parent tumor, opening a window on the effects of tumor-host interactions on behavior. Such comparison is not possible in samples excised from clinical tumor specimens. Technical advances are also incorporated in this work: tumor cell lineages were labeled with green fluorescent protein (GFP) to enhance accuracy of selection of primary and secondary tumor tissue for analysis. Also, to evaluate the initial screening results, we conducted extensive laboratory studies and computational (training and test) procedures and validated the expression levels of a number of genes of interest.

Our findings indicated that the genes expressed in primary tumors generated by metastasis-competent cell populations differed clearly from those in their metastasis-incompetent counterparts. In contrast, the patterns in metastases were similar to the primary tumors from which they originated, and metastases in the lungs displayed remarkably similar gene expression patterns to those in the lymph nodes of the same animal. Additionally, the patterns observed in tumors and metastases differed from those seen in the parental cell lines in vitro, indicating that the host microenvironment is an active participant in tumor progression and metastasis. These findings have significant implications for defining mechanisms of metastasis and for designing novel effective therapy. They also contribute a manageable list of candidate genes from which to choose targets for interventional studies on mechanisms involved in the process.

	Materials and Methods

Top Abstract Materials and Methods Results Discussion References

 
Cell Lines

The NM-2C5 and M-4A4 lines were isolated in our laboratory from the MDA-MD-435 breast cancer cell line as described by Bao et al⁶ and subsequently transduced with an enhanced green fluorescent protein-expressing vector.¹³ These monoclonal cell lines were routinely cultured in RPMI 1640 medium supplemented with 10% newborn calf serum (Invitrogen, Carlsbad, CA), penicillin, and streptomycin at 37°C in a humidified atmosphere of 5% CO₂-95% air. Cell line LM3 was derived from a lung metastasis produced by M-4A4, and cell line CL16 was obtained from a lung metastasis made by descendants of LM3 after two more similar selection cycles in nude mice. Both are progressively more metastatic variants of the parent line (see below) grown under the same conditions in vitro.

Murine Xenograft Metastasis Model

One million cells in 50 µl of a mixture of RPMI 1640 medium and ECM gel (Sigma Chemical Co., St. Louis, MO) were inoculated into the mammary fat pad of anesthetized mice. Animals were euthanized and autopsied at 3 to 4 months postinoculation when the primary tumors reached 20 mm in diameter. Metastasis formation was assessed by macroscopic observation of all major organs for secondary tumors and confirmed by histological examination. Metastasis was also confirmed by looking for fluorescence of incorporated GFP under blue light ( = 490 nm), which is sufficiently sensitive to detect single cells. Only cell clusters (>1 mm) are regarded as true metastases. Tissues from primary tumors and metastases were snap-frozen and stored at –80°C until used for RNA or protein extraction.

Protein Analyses

Protein detection and quantification were performed either on primary tumors or sera from tumor-bearing mice depending on the protein localization. Tumor homogenates were prepared in 20 mmol/L Tris-HCl, pH 8.0, 5 mmol/L CaCl₂, 1 mmol/L phenylmethylsulfonyl fluoride, 15 µmol/L pepstatin A, and 0.05% (w/v) Brij 35. The Complete EDTA-free protease inhibitor cocktail (Boehringer Ingelheim Chemicals, Petersburg, VA) was added to the extraction buffer. The tissues were homogenized on ice, and the homogenates were centrifuged at 14,000 x g for 40 minutes (4°C). Protein quantitation of the supernatants was performed using the Coomassie Plus Protein Assay kit (Pierce, Rockford, IL). Twenty micrograms of denatured protein samples was separated on a 12.5% sodium dodecyl sulfate-polyacrylamide gel electrophoresis glycine gel, and proteins were transferred onto polyvinylidene difluoride membranes using a semi-dry apparatus (Bio-Rad, Life Science, Hercules, CA) according to the manufacturer’s instructions. Immunodetection was performed as described previously⁷ using the following specific antibodies: monoclonal anti-silver homolog (Pmel-17) (Neomarkers, Inc., Fremont, CA); monoclonal anti-MITF (Neomarkers, Inc.); polyclonal anti-1-antichymotrypsin (DAKO, Carpinteria, CA); polyclonal anti-TRP1 (Santa Cruz Biotechnology, Santa Cruz, CA); monoclonal anti-osteopontin (Chemicon International, Inc., Temecula, CA); monoclonal anti-thrombospondin (BD Transduction Laboratories, San Jose, CA); and polyclonal anti-MMP-8 (Chemicon International, Inc). MMP-8 and OPN protein expressions were quantified using commercially available enzyme-linked immunosorbent assay (ELISA) kits from Amersham Biosciences (San Francisco, CA) and Assay Designs (Ann Arbor, MI), respectively.

Frozen sections of xenograft primary tumors were fixed in 4% buffered formaldehyde. Antigen retrieval was performed using the target retrieval solution at 1:10 dilution (DAKO). Endogenous peroxidase activity was blocked by incubation in 3% H₂O₂. Nonspecific binding of the antibodies to irrelevant proteins was blocked by incubation in 10% goat serum. Proteins of interest were targeted with the same antibodies used in the Western-blot experiments (see above). When monoclonal mouse antibodies were applied, a prior incubation of the section with goat anti-mouse Ig Fab fragments (Jackson Immunoresearch Laboratories, West Grove, PA), which neutralized the Fc domain reactivity of the endogenous host Ig, was performed. The horseradish peroxidase-conjugated secondary antibody was visualized by diaminobenzidine substrate-chromogen 3,3'-diaminobenzidine (DAKO). Sections were counterstained with Hematoxylin-Gills No. 2 solution, dehydrated in alcohol, cleared in xylene, and mounted in Permount (Fisher Chemicals, Lake Forest, CA).

cRNA Preparation and GeneChip Hybridization

Total RNA was extracted from cultured cells and frozen tissue samples with TRIzol reagent (Invitrogen) and cleaned with the DNA-free kit (Ambion, Austin, TX). RNA quality was assessed by running the samples on a native 1% agarose gel and on a Biogem analyzer (Agilent, Palo Alto, CA). cRNA was prepared in the University of California-San Diego Cancer Center Microarray facility as described by the standard Affymetrix microarray protocols. The cRNA was then hybridized to human HG-U133A GeneChip oligonucleotide arrays (Affymetrix, Santa Clara, CA), which interrogate approximately 22,000 transcripts. The arrays were scanned at 560 nm using an argon-ion confocal laser as the excitation source.

Microarray Data Analysis

The DAT files containing the scanned images of each microarray were individually inspected for quality control and digitized by Microarray Analysis Suite 5.0 (Affymetrix). The resultant CEL files containing the raw numerical data for signal intensity at probe level were collectively read and analyzed in dChip software.^14,15 Briefly, each microarray was normalized against a common baseline array using the "invariant probe set" method. After normalization, the model-based expression index of each gene was then calculated according to the PM-MM model. To identify candidate genes that were differentially expressed between any two group of arrays, a screening filter consisting of the following criteria was applied: 1) a fold change (fc) larger than 1.5 or 3; 2) two-tailed P values (paired t-test if applicable) smaller than 0.05; and 3) a minimal difference of 100 between the group mean of model-based expression index. The resultant lists of candidate genes were then sorted according to their corresponding P values with a cut-off of 0.005 to ensure high stringency of the analysis (Tables 1–6) . For subsequent high-level analysis, candidate gene lists from six group-wise comparisons were combined and subjected to hierarchical clustering (centroid-linkage) or classification by linear discriminant analysis (LDA) within the R environment (http://www.R-project.org).

mRNA from the same total RNA samples used for the microarray analyses was reverse transcribed using M-MLV reverse transcriptase and oligo(dT) from the Retroscript cDNA synthesis system (Ambion). The amplification reactions were conducted in 96-well plates in 25-µl reaction volumes containing 12.5 µl of 2X SYBR Green Master Mix (PE Applied Biosystems, Foster City, CA), 50 nmol/L each of forward and reverse primers, and 1 µl of the cDNA and monitored in an ABI Prism 7700 Sequence Detector System (PE Applied Biosystems). The thermal profile for the PCR was 50°C for 2 minutes and 95°C for 10 minutes followed by 40 cycles of 95°C for 10 seconds (denaturation step) and 60°C for 1 minute (annealing and elongation steps). Measurements on each sample were performed in triplicate, and the expression of the tested gene was normalized to a GAPDH standard curve run in duplicates on the same plate.

	Results

Top Abstract Materials and Methods Results Discussion References

 
Our investigative strategy compared the gene expression patterns associated with the metastatic behaviors of three isogenic human breast tumor cell lines and the primary tumors that they generated after orthotopic inoculation in nude or SCID mice. It also examined and compared gene expression profiles of metastases in various organs with each other and with the primary tumors. One of these lines generates tumors that are essentially nonmetastatic (NM-2C5), whereas the other two produce ones that are moderately (LM3) or highly (CL16) metastatic to the lungs and lymph nodes (Figure 1) . The differences are great in magnitude and obvious but not absolute. In some batches of mice inoculated with NM-2C5, occasional metastases are seen, whereas in other batches, there are none. For convenience and brevity, we shall refer to metastatic or nonmetastatic primary tumors. LM3 and CL16 are third and fifth generation descendents of the M-4A4 metastatic line originally isolated by Bao et al⁶ from MDA-MB-435 and were obtained by cyclically culturing and orthotopically re-inoculating the cells of successive generations of metastases as originally described by Fidler.¹⁶ The degree of metastatic aggressiveness attained by the fifth cycle compared with ancestral lines is readily seen from the overwhelming burden of fluorescent cancer cells colonizing the lungs and lymph nodes (Figure 1) .

View larger version (73K): [in this window] [in a new window]   Figure 1. Visualization of GFP-labeled tumors and metastases in our xenogeneic breast cancer metastasis model. Left: Increasing metastatic loads in the five pulmonary lobes of mice bearing 2-cm primary tumors developed after orthotopic inoculation of isogenic clonal cell lines with progressively enriched metastatic capabilities, as indicated. Right: A survey view of the heavy metastatic spread of a CL16 primary tumor in the whole animal. Numerous tumor cells leaving the primary tumor (P) distended the afferent lymphatic vessels (arrowheads). Pulmonary and lymphatic metastases are also indicated. An asterisk indicates the heart.

The tissue samples were obtained from two animal experiments; one in severe combined immunodeficient (SCID) mice and a second in nude mice. In SCIDs, the metastatic line CL16 is overwhelmingly more metastatic than in nude mice, although the metastatic capability of NM-2C5 remains low in both strains. The SCID mouse study was, therefore, the primary experiment, and the one in nude mice was a backup for testing conclusions in another strain. In tissues from SCID mice, the expression levels of 22,000 genes were screened in three NM-2C5 and five CL16 primary tumors as well as dissected and cleaned lung and thoracic lymph node metastases from the same five animals bearing the CL16 tumors. RNA from each primary tumor or metastasis sample was hybridized to a separate individual GeneChip, comprising a total of 18 microarrays on SCID samples (see supplemental material at http://ajp.amjpathol.org). The studies on samples from nude mice were conducted on 4 microarrays hybridized with primary NM-2C5 tumors, 7 arrays with primary LM3 or CL16 tumors, and 10 with lung metastases, totaling 21 arrays. In addition, we screened gene expression patterns in the parent NM-2C5 and M-4A4 cell lines and in CL16 in vitro using RNA from three distinct cultures of each line pooled and hybridized to a separate individual array. The grand total of arrays used in the work described in this communication, therefore, amounts to 42.

To evaluate the degree of consistency of the global gene expression profiles among samples in each biological category, chip-to-chip comparisons of the signal intensities from the complete set of genes present on the array were performed for the data from the SCID animals, which we later used to train the algorithms we used to test the data from the nude animals. The scatter plots obtained by comparing data from each SCID microarray with the others from the same tumor category (ie Mor NM primaries or metastases) showed tight grouping of the points, and the correlation coefficients ranged between 0.89 and 0.98, with the majority between 0.97 and 0.98. These analyses, using thousands of data points from each chip, established that the biological samples in each category showed good consistency and were suitable for more detailed studies.

Quantitative PCR Validation of a Subset of Differentially Expressed Genes Selected from the Microarray Data Analyses

To evaluate the microarray results using a different method, we selected 41 genes of interest to us, from the 22,000 present on the HG-U133A chip, and quantified their expression in vitro (cell lines) and in vivo (primary tumor/metastasis samples) by real-time PCR using human-specific primers. The genes chosen for validation by this independent technology covered the low- and high-signal-intensity range and included 37 differentially expressed genes in one or the other of the pathologically relevant comparisons (nonmetastatic versus metastatic primary tumor or metastases versus primary), as well as four genes that were not differentially expressed (fc <1.5). Numerical data and gene identities are available online on the Web site (http://ajp.amjpathol.org).

The real-time PCR data for the cultured cell lines showed that 76% of the microarray fc were confirmed to be in the same trend and that only 9% had a trend in the opposite direction. The remaining 15% were distributed between 12% false positive (fc >1.5 detected only by microarray quantification) or 3% false negative (fc >1.5 detected only by quantitative PCR (Q-PCR)) categories. It was also found that, because Q-PCR had much greater sensitivity and dynamic range of quantitation, the magnitudes of fold changes of differentially expressed genes were underestimated by microarray analysis. Similar analysis of the datasets for the tumor tissue samples in vivo showed that 49% of the fold changes identified as differentially regulated by microarrays were validated by Q-PCR and that 15% displayed a reverse trend. The frequency of false-positive results was approximately similar to that seen in the pure cell lines, but the false-negative results increased from 3 to 22%. Because the chips used in all our experiments were of the same type, the difference in validated fold changes identified by them in cells versus tissues is attributed to a confounding effect from host (mouse) transcripts in the latter. The value of 76% fc validated on pure human cells with an additional 12% false positive by microarray indicates that the technique is not missing much important information.

Collectively, these results indicate that oligonucleotide microarray technology is a valuable screening tool for selecting potentially interesting candidate genes for further study, but real-time PCR validation is an essential requirement, especially in a xenogeneic system in which human mRNA is co-extracted with mouse (host) transcripts from the intermingled cell populations in the tumors and the metastases. Our validation Q-PCR studies, using confirmed human-specific primers, established that many of the differences in gene expression revealed by HG-U133A microarray analysis between nonmetastatic and metastatic primary tumors and distant metastases were attributable to the human tumor cells within them. It should be noted that, to be very conservative in interpretation, any differential expression of less than fc 1.5 measured by Q-PCR in our validation experiments was recorded as not validated. Therefore, although differential changes seen by microarray were indicated as correct by the Q-PCR readings for many genes, we chose to exclude them in case they were within the margin of experimental error.

Confirmation of differential gene expression at the protein level is also desirable, and we did this by Western blot, ELISA, or immunohistochemistry (Figure 2) for seven of these genes (TYRP-1, MITF, TSP-1, Pmel-17, MMP-8, ACT, and OPN) using antibodies raised against human antigens. Several of these proteins (OPN, ACT, MMP-8, and TSP-1) were differentially detectable in the serum of the tumor-bearing animals, indicating the possibility of detecting potential biomarkers by this approach.

View larger version (92K): [in this window] [in a new window]   Figure 2. Protein validations of candidates selected on microarray-based differential expression between the nonmetastatic and the metastatic tumors in BalbC mice. Changes in mRNA expressions estimated by microarray analyses and quantified by real-time PCR (middle columns) were validated at the protein level by quantitative (ELISA, third column) or semiquantitative methods (Western blotting, left; immunohistochemistry, right) as described in Materials and Methods. Both the "Affy" and the "Q-PCR" fold changes came from the Q-PCR validation table provided in the supplemental material. This explains why, for some genes, the microarray data indicated in this figure may not exactly match those in Table 1 , which are based on SCID comparisons. All of the stained tumor sections were observed at 400x magnification, and a scale bar was inserted in the NM-2C5 tumor section stained with the anti-Pmel-17 antibody. Arrows in the upper left show examples of positive signals of MITF localization in the nuclei of nonmetastatic tumor cells.

Microarray Experiment Design and Comparative Analysis

To enhance discrimination of the metastasis-related signals from nonspecific background expression due to the biological variations of an in vivo system, a "training" data set derived from CL16 tumors and their matched metastases to lymph nodes or lungs from the SCID experiments was first used for preliminary comparisons. The results were then validated with the "test" data set consisting of biological samples from a different strain of host, namely nude mice. This policy enabled us to perform paired t-tests with more statistical power to eliminate gene expression due to individual variation that is irrelevant to metastasis.

All biologically relevant comparisons illustrated in Figure 3 were performed using dChip software as described in Materials and Methods. The comparison of expression levels for each gene between two biological sample groups generated fold changes with attached P values. Instead of focusing on high fold change based on group means, we opted to rank-order the candidate gene lists by high consistency of differential expression (low P values) in specimens from a given category (Tables 1–6) . This approach of emphasizing statistical measures of data consistency over fold change is a defining feature of this work, a direct benefit of the tightly controlled experiment design. In addition, it was chosen because of the limited dynamic range and accuracy of the microarray results, revealed by the Q-PCR validation experiment, and our intention to have a high level of confidence in the candidate genes selected for further study.

View larger version (16K): [in this window] [in a new window]   Figure 3. Schematic diagram of the "horizontal" (black arrows) and "vertical" (white arrows) comparisons performed in the microarray data analyses. dChip software14,15 was used to compare global expression datasets obtained with the HG-U133A array hybridized with the nonmetastatic (NM-2C5) and metastatic (CL16) primary tumors and the spontaneous metastatic deposits developing in the lung (Lung Mets) or in the thoracic lymph nodes (LN Mets) in SCID mice.

Screening among the 22,000 features on the human HG-U133A chip for differentially expressed genes between the NM-2C5 and CL16 tumors in the SCID mice resulted in a list of candidate genes ranging from 72 (55+17) to 1765 (1068+697) genes when different filtering criteria were applied (Table 2) . To evaluate the biological significance of the resultant gene lists, we performed linear discriminant analysis (LDA) to evaluate whether each list would classify additional test samples (from the experiment in nude mice) into the nonmetastatic group or the metastatic group (Figure 4C) correctly. The results suggested that at P < 0.005 and fc >1.5, the classification of the 17 samples from a separate (duplicate) experiment in a different mouse strain was close to optimal, with only two arrays being borderline misclassified (93.5% accuracy). In the corresponding gene list (P < 0.005), 386 candidates were found to be up-regulated more than 1.5 times in the nonmetastatic primary tumor compared with its metastatic counterpart, and 416 were up-regulated in the metastatic primary tumor (Table 2) . Thus, of a total of 802 genes that are consistently differentially expressed in the CL16 tumors, approximately one-half of them have to be down-regulated, and the other one-half, up-regulated. It follows that some genes must be "turned on," but also that others need to be "turned off," for the metastatic phenotype to be triggered. This observation highlights the potential importance of negative regulators (metastasis suppressor) in metastasis.

View larger version (46K): [in this window] [in a new window]   Figure 4. Signature gene clustering and sample classification. Hierarchical cluster analyses of the differentially expressed genes from in vivo and in vitro comparisons. A: Candidate genes differentially expressed (with P < 0.005 and fc >1.5) at least in one of the six in vivo comparisons (illustrated in Figure 3 ) were clustered according to their expression levels in the NM-2C5 and CL16 lines in vitro and in the primary and secondary tumors they generated in SCID mice. B: Candidate genes differentially expressed (with fc >3) at least in one of the three in vitro comparisons of the cell lines (Table 7) were clustered according to their expression levels in NM-2C5, LM3, and CL16 lines. The vertical bars on right of this panel indicate the candidate genes gradually overexpressed (white bar) or down-regulated (gray bar) in regard with the progressive acquirement of high metastatic capabilities. As shown in the color bar, red indicates high expression; green, low expression; and black, intermediate expression. The dendogram on the left indicates the pairing of genes, and the branch length is proportional to the distances between the clusters. C: Linear discriminant analysis using NM-2C5 (blue striped squares) and CL16 (red striped squares) primary tumors in SCID mice as the training set to classify unknown samples (test set) according to the combined candidate gene list used in A for the clustering analysis. Plain blue squares: Samples classified into the nonmetastatic group; plain red squares: samples classified into the metastatic group. The unknown samples include primary tumors (NM-2C5 and LM3), metastases in the lung (Lg) or lymph nodes (LN), and cell lines (underlined). The asterisk shows misclassified samples.

An identical statistical strategy was used in a vertical comparison (Figure 3) between the primary (CL16) tumors and their metastatic deposits in the lungs and lymph nodes. The term "vertical" is used to describe a spatio-temporal comparison between a sample and its biological progeny located in a different organ. The tumor cells constituting the secondary deposits are the direct descendants of the experimentally selected clone of neoplastic cells forming the primary tumor and are, therefore, genetically almost identical to them. Theoretically, one might speculate that the changes that generate metastasis might be evident in both populations and that their expression profiles might be very similar. Indeed, we found (Table 2) that their expression profiles showed little differences in contrast to when the metastatic and nonmetastatic primary tumors were compared (ie, in the horizontal comparison): only 26 (17+49) genes were more than 1.5-fold differentially expressed between lung metastases and their corresponding CL16 primaries at the high level of significance (P < 0.005). For lymph node metastases, the corresponding figure was 114 (56+58) genes.

Because the original breast tumor from which MDA-MB-435 was derived was composed of a mixture of metastatic, less metastatic, and nonmetastatic tumor cell populations, we also compared the expression profiles of metastases with that of sister (isogenic) cells in nonmetastatic NM-2C5 primary tumors, growing in a different set of mice of the same strain (Table 3) . This comparison is clinically and pathogenetically relevant, because it represents the extreme ends of the phenotype but cannot be performed in "wild-type" tumors and is, uniquely, only possible in this experimental system. It would be expected to reveal metastasis-relevant genes that are expressed in both primary and secondary tumors of the metastatic phenotype (CL16) and therefore are not seen as differentially expressed when these two categories are compared. Surprisingly, both lung and lymph node metastases harbor slightly less differentially expressed genes than the parental metastatic primary tumors (CL16) when compared with the NM-2C5 tumors. This result suggests that the CL16 cells residing in a primary tumor have already acquired an enriched metastatic gene expression profile as a result of the cycling and selection procedures and that this profile is inherited by the cells in the metastases.

Combination of Comparisons

Each of the pair-wise comparisons performed above only represents a segment of the "metastatic spectrum." Therefore, we combined all six candidate gene lists to create a working database comprising 1127 differentially expressed genes putatively involved in progression from a basal to an enhanced metastatic potential in the primary tumor, extending to a consummated metastatic profile in lymph node and lung metastases. To visualize changing expression patterns of these genes accompanying transitions between tumors of low, higher, and highest metastatic potency, we performed supervised hierarchical clustering of the biological samples using the "pooled" candidate gene list (Figure 4A) . This procedure graphically demonstrated the similarities between prevailing expression patterns in the cells in metastatic lesions in different organs and in their parent tumors as well as the differences from their phenotypically opposite counterparts. Additionally, we included in this analysis the expression data from cultures of NM-2C5 and CL16 cell lines in vitro. This permitted side-by-side comparison with the in vivo data and indicated that several genes were regulated by the host microenvironment while many others were expressed at similar levels both in vitro and in vivo. Within the NM-2C5 category, the expression levels of these 1127 pooled candidate genes were more stable and conserved between in vivo and in vitro samples. In contrast, the CL16 samples displayed greater magnitudes of changes in expression levels for these genes, and approximately one-half of them were highly inducible by the host microenvironment. Therefore, our in silico reconstruction of tumor progression toward metastasis indicates that the M-4A4 family of cells (ie, M-4A4, LM-3, and CL16) responded more dynamically to the environment in vivo and developed a characteristic signature that segregated with the metastatic phenotype during accentuation by cyclical selection and reinoculation.

Comparisons of Cell Lines in Vitro

To ascertain whether the changing pattern of expression seen in primary tumors and metastases, as the phenotype became enriched, reflected inherited intrinsic changes in the tumor cells as they were cycled, we compared data from cell lines NM-2C5, M-4A4, and CL16 cells growing in vitro. In this situation, there are no intermingled mouse cells, and the results reflect gene expression patterns in "purified" human cells of these tumor lines growing under essentially identical conditions. It was found that 724 genes (Table 2) were differentially expressed by 1.5-fold or more between the first generation M-4A4 and NM-2C5 cell lines. This number increased to 2638 when the fifth generation CL16 was similarly compared with the NM-2C5 line. More importantly, 377 of the total 724 genes generated by the first comparison were also found in the second comparison. In addition, very few of these 377 genes were regulated in opposite directions in M-4A4 and CL16, consistent with the notion that specific gene expression was segregated along with the metastatic phenotype. Hierarchical clustering of the lines demonstrated the inheritance of the changing gene expression pattern as the cell lines evolved more metastatic phenotype (Figure 4B) . Therefore, in the absence of host influence, comparisons of autonomous gene expression profiles reveal a group of human genes that likely represent an essential requirement for metastasis.

	Discussion

Top Abstract Materials and Methods Results Discussion References

 
Although comprising a powerful investigation tool, microarray experiments provide, by design, a summary analysis of (biologically) averaged gene expression profiles. As a result, microarray studies conducted on clinical cancers are often obscured by the uneven amount of cancerous tissues within samples that were harvested from patients of unmatched genetic backgrounds. For the study of metastasis-related gene expression, the difficulties are further compounded by tumor heterogeneity, which translates into uncertainty about the ratio of cell populations in the primary tumor that inherited variable genetic potential for metastasis.

In this study, we therefore designed the biological system to control, if not overcome, the above problems. By using clonal isogenic tumor cell lines of divergent metastatic performance, we focused the enquiry on metastasis-related differential expression and enhanced the potency of the metastatic line by recycling it several times. Also, the use of GFP-labeled human cell lines in immunocompromised unlabeled mice improved the identification of metastases while minimizing contamination by irrelevant tissues during sample dissection, hence enriching gene expression truly associated with the tumor or metastases. On the other hand, the host genetic backgrounds were carefully controlled by performing independent animal experiments in two different stains of mice.

The main findings emerging from this work are 1) that gene expression in clonal metastatic primary breast tumors differed clearly from that seen in isogenic nonmetastatic tumors generated by a different clonal cell population isolated from the same patient; and 2) that the expression patterns in matched pulmonary and lymph nodal metastases closely resembled those in the primary tumor within the same animal. Hierarchical clustering and linear discriminant analysis of the differentially expressed genes from all of the comparisons in two successive experiments in different mouse strains cross-corroborated each other: that metastatic primary tumors and their metastases were very distinct from nonmetastatic primary tumors. Thus the evidence that emerges from this specially designed study of human cancer cells favors the view that cells with a distinct gene expression signature associated with metastasis can be isolated from human neoplasms and that this profile faithfully segregates with the phenotype during cyclical selection procedures designed to enrich metastatic capability.^4,17-22 The tumor cells within the primary tumors and metastases in SCIDs were in the fifth cycle of re-derivation and orthotopic inoculation, and the relatively small numbers of differences between lung and lymph node metastasis signatures suggest that many similar biological mechanisms are being executed in reaching and growing in these two different favored sites. Expression profiles of matched bone, brain, or liver metastases would inform us whether the results obtained in this study apply to other preferred sites for human breast cancer metastasis. Unfortunately, these samples are rarely available in this model. Collectively, this is potentially valuable information for further investigations aiming to find prognostic markers and therapeutic targets.

Although the identities of the candidate genes that may be "driving" the process are of substantial interest, one cannot be sure that they can be extrapolated to naturally occurring "wild-type" human cancers. We believe that, at present, the observations are informative only about the system under study, but they do provide potentially valuable leads about genes and pathways that merit functional validation in the experimental model and further investigation in pure human tumors. Also, the candidates, which we provide in Table 2 and on the study Web site, are composite collections of genes that are differentially expressed in the tumor cells and in the host tissues intermingled with them. Thus, some of the results obtained will be reflecting differences between the mouse stroma in the tumors in different organs. The investigative power of this xenogeneic system lies in the opportunity now available to identify which components of this joint signature are contributed by the tumor cells and which by the host stroma, by the use of species-specific primers and antibodies. Earlier studies^23-25 involving orthotopic versus ectopic inoculation of carcinoma cells have demonstrated that local host tissues influence metastatic behavior by metastasis-competent cells, and opportunities to analyze host gene expression in the process are therefore valuable. In this way, the present work contributes a foundation for the further study of tumor-stromal interactions in metastasis and for the identification of intrinsic tumor cell genes driving the process.

It is appropriate to compare these findings with previous microarray studies describing gene expression profiles in metastases and primary human tumors. In general, such work has been hampered by the limited availability of suitable human tissue samples. However, lately, two major studies investigated gene expression patterns in metastases from carcinomas of the prostate⁸ and from primary adenocarcinomas of different organs⁹ and compared the results with those in corresponding primary tumors from other patients (ie, unmatched primaries). Both concluded that secondary tumors display similar transcriptional profiles to the primary neoplasm, in agreement with Adib et al,²⁶ who described observations on primary ovarian cancers compared with matched omental metastases (although metastasis was hematogenous in the first two studies cited but transcoelomic in the third), and with Weigelt et al,²⁷ , who compared individual hematogenous metastases excised from various organs of eight patients with their primary breast cancers. Conversely, two earlier studies on prostate cancer metastasis reported that expression in metastatic samples differed from that in the primary tumor,^11,12 although, once again, the primaries and metastases were not matched. The discordance of these observations might result either from primary tumor cell heterogeneity, resulting from differences in the time of emergence of the metastatic cell phenotype, or from comparison of metastases and primary tumors obtained from patients with different genotypes. Both issues were bypassed in our model.

Recently Lee and coworkers²⁸ described a study of xenografted tumors profiled with cDNA spotted microarrays comprising 5800 known genes on glass slides. They had derived two MDA-MB-435 variant cell lines, LN435 and Tho435, which preferentially metastasize to the lymph nodes or to the thoracic cavity, respectively, when orthotopically injected in SCID mice. Gene expression profiles of the primary tumors that they generated after inoculation were compared with the profiles of primary tumors generated by the parent MDA-MB-435 line to derive candidate genes potentially involved in organ specificity of breast cancer metastasis. Although the experimental strategy of our work differed from that of Lee et al,²⁸ by profiling metastatic deposits taken directly from the lungs and lymph nodes and comparing their expression patterns with their isogenic primary tumors in the mammary gland, using 22,000 gene Affymetrix arrays, there are several overlaps between the gene lists provided in the Lee et al²⁸ study and in our own. Despite many conceptual differences between these two studies based on the MDA-MB-435 breast cancer line, the cross-corroboration of a number of interesting candidates by separate laboratories strongly supports their involvement in the metastatic process.

The pioneering work of Fidler^1,2 and of many other subsequent investigators^4,22,29-34 has established beyond any reasonable doubt that tumors are composed of heterogeneous populations that differ from each other in many qualities including metastatic capability. This has been confirmed by pathologists who have shown that human tumors show zonal heterogeneity with regard to pigmentation, fibrosis, vascular supply, and many other properties. Our own work, reported here and previously,^6,7 shows conclusively that the parent (MDA-MB-435) cell line from which we obtained NM-2C5 and M-4A4 clones was heterogeneous and that it contained nonmetastatic cell populations. Thus, the deconstruction exercise by which we derived these clones with polar opposite behavior and the subsequent screening process by which we chose the most suitable tumors for study facilitated and clarified our gene expression analysis. This helps to illuminate how investigations on unmatched and unselected fresh wild-type human tumors may give variable results.

In conclusion, this work, based on oligonucleotide microarray gene profiling, adopted a purpose-designed, new strategy to define patterns of gene expression relevant to the dynamic process of tumor metastasis, using a unique isogenic breast cancer model. The gene-sorting strategy was designed to provide very robust candidates (P < 0.005), even with low fold change values (fc >1.5). We recognize that this statistical stringency could potentially exclude some genes of interest, but we concluded that combining it with pair-wise analysis of signatures from other neoplastic lesions in the same animal offered the best chance of serving our intention to identify genes worthy of further work. Moreover, if this strategy was to exclude an important gene related to metastasis, the global microarray screening technique, by nature, would detect related changes in its relevant gene network and thereby could lead to the identification of the missing candidate. Statistical "training and test" validation of the observations made on repeat experiments in different mouse strains together with biochemical validation of expression of 41 genes constitute further novel aspects of this study building upon and extending previous investigations by other laboratories. These results show the power of microarray analysis for initial high-throughput screening while demonstrating the necessity of corroboration with other techniques. The combination of the methods described enables one to have reasonable confidence that some of the genes driving metastasis lie within these lists, although further sifting is needed. The study design permitted previously impossible comparisons to be made (such as between metastases and isogenic nonmetastatic breast cancers) and revealed molecular information pertinent to the emergence and to the ongoing cyclical continuation of the metastatic phenotype. The main purpose of the work has been to design and use a new tool for the identification of genes that are of biological or clinical interest. The candidate genes so obtained can be screened for new information on possible pathways involved and tested in knock-down and overexpression assays, using this model as well as assayed in fresh human tumor samples.

The new information provided here, showing the importance of the inductive effects of the stroma of the host breast on gene expression by mammary carcinoma cells that are about to metastasize and the preservation of this distinctive expression pattern when their descendants make deposits in the lungs and lymph nodes, opens new insights into the mechanisms of the biology of the phenomenon and fresh possibilities for finding prognostic markers and therapeutic targets.

	Acknowledgements

	Footnotes

 
Address reprint requests to David Tarin, Rebecca and John Moores Comprehensive Cancer Center and Department of Pathology, University of California, San Diego, 9500 Gilman Drive, M/C 0912 La Jolla, CA 92093-0912. E-mail: dtarin{at}ucsd.edu

Supported by a grant from the Loppieda Foundation.

V.M. and T.-Y.H. contributed equally to this work.

Accepted for publication February 1, 2005.

	References

Top Abstract Materials and Methods Results Discussion References

 

1. Fidler IJ, Kripke ML: Metastasis results from preexisting variant cells within a malignant tumor. Science 1977, 197:893-895[Abstract/Free Full Text]
2. Fidler IJ: Tumor heterogeneity and the biology of cancer invasion and metastasis. Cancer Res 1978, 38:2651-2660[Abstract/Free Full Text]
3. Fidler IJ: The pathogenesis of cancer metastasis: the ‘seed and soil’ hypothesis revisited. Nat Rev Cancer 2003, 3:453-458[Medline]
4. Price JE, Carr D, Tarin D: Spontaneous and induced metastasis of naturally occurring tumors in mice: analysis of cell shedding into the blood. J Natl Cancer Inst 1984, 73:1319-1326
5. Tarin D, Price JE, Kettlewell MG, Souter RG, Vass AC, Crossley B: Mechanisms of human tumor metastasis studied in patients with peritoneovenous shunts. Cancer Res 1984, 44:3584-3592[Abstract/Free Full Text]
6. Bao L, Pigott R, Matsumura Y, Baban D, Tarin D: Correlation of VLA-4 integrin expression with metastatic potential in various human tumour cell lines. Differentiation 1993, 52:239-246[Medline]
7. Urquidi V, Sloan D, Kawai K, Agarwal D, Woodman AC, Tarin D, Goodison S: Contrasting expression of thrombospondin-1 and osteopontin correlates with absence or presence of metastatic phenotype in an isogenic model of spontaneous human breast cancer metastasis. Clin Cancer Res 2002, 8:61-74[Abstract/Free Full Text]
8. Lapointe J, Li C, Higgins JP, van de Rijn M, Bair E, Montgomery K, Ferrari M, Egevad L, Rayford W, Bergerheim U, Ekman P, DeMarzo AM, Tibshirani R, Botstein D, Brown PO, Brooks JD, Pollack JR: Gene expression profiling identifies clinically relevant subtypes of prostate cancer. Proc Natl Acad Sci USA 2004, 101:811-816[Abstract/Free Full Text]
9. Ramaswamy S, Ross KN, Lander ES, Golub TR: A molecular signature of metastasis in primary solid tumors. Nat Genet 2003, 33:49-54[Medline]
10. van de Vijver MJ, He YD, van’t Veer LJ, Dai H, Hart AA, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ, Parrish M, Atsma D, Witteveen A, Glas A, Delahaye L, van der Velde T, Bartelink H, Rodenhuis S, Rutgers ET, Friend SH, Bernards R: A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med 2002, 347:1999-2009[Abstract/Free Full Text]
11. LaTulippe E, Satagopan J, Smith A, Scher H, Scardino P, Reuter V, Gerald WL: Comprehensive gene expression analysis of prostate cancer reveals distinct transcriptional programs associated with metastatic disease. Cancer Res 2002, 62:4499-4506[Abstract/Free Full Text]
12. Dhanasekaran SM, Barrette TR, Ghosh D, Shah R, Varambally S, Kurachi K, Pienta KJ, Rubin MA, Chinnaiyan AM: Delineation of prognostic biomarkers in prostate cancer. Nature 2001, 412:822-826[Medline]
13. Goodison S, Kawai K, Hihara J, Jiang P, Yang M, Urquidi V, Hoffman RM, Tarin D: Prolonged dormancy and site-specific growth potential of cancer cells spontaneously disseminated from nonmetastatic breast tumors as revealed by labeling with green fluorescent protein. Clin Cancer Res 2003, 9:3808-3814[Abstract/Free Full Text]
14. Li C, Wong WH: Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci USA 2001, 98:31-36[Abstract/Free Full Text]
15. Li C, Wong WH: Model-based analysis of oligonucleotide arrays: model validation, design issues and standard error application. Genome Biol 2001, 2:RESEARCH0032. Epub
16. Fidler IJ: Selection of successive tumour lines for metastasis. Nat New Biol 1973, 242:148-149[Medline]
17. Brunson KW, Beattie G, Nicolsin GL: Selection and altered properties of brain-colonising metastatic melanoma. Nature 1978, 272:543-545[Medline]
18. Tao T, Matter A, Vogel K, Burger MM: Liver-colonizing melanoma cells selected from B-16 melanoma. Int J Cancer 1979, 23:854-857[Medline]
19. Edel G, Grundmann E: Selection of liver-colonizing tumor cells from a murine fibrosarcoma induced by methylcholanthrene. J Cancer Res Clin Oncol 1984, 108:274-280[Medline]
20. Pettaway CA, Pathak S, Greene G, Ramirez E, Wilson MR, Killion JJ, Fidler IJ: Selection of highly metastatic variants of different human prostatic carcinomas using orthotopic implantation in nude mice. Clin Cancer Res 1996, 2:1627-1636[Abstract]
21. Barroga EF, Kadosawa T, Okumura M, Fujinaga T: Establishment and characterization of the growth and pulmonary metastasis of a highly lung metastasizing cell line from canine osteosarcoma in nude mice. J Vet Med Sci 1999, 61:361-367[Medline]
22. Price JE, Carr D, Jones LD, Messer P, Tarin D: Experimental analysis of factors affecting metastatic spread using naturally occurring tumours. Invasion Metastasis 1982, 2:77-112[Medline]
23. Price JE, Polyzos A, Zhang RD, Daniels LM: Tumorigenicity and metastasis of human breast carcinoma cell lines in nude mice. Cancer Res 1990, 50:717-721[Abstract/Free Full Text]
24. Morikawa K, Walker SM, Nakajima M, Pathak S, Jessup JM, Fidler IJ: Influence of organ environment on the growth, selection, and metastasis of human colon carcinoma cells in nude mice. Cancer Res 1988, 48:6863-6871[Abstract/Free Full Text]
25. Bao L, Matsumura Y, Baban D, Sun Y, Tarin D: Effects of inoculation site and Matrigel on growth and metastasis of human breast cancer cells. Br J Cancer 1994, 70:228-232[Medline]
26. Adib TR, Henderson S, Perrett C, Hewitt D, Bourmpoulia D, Ledermann J, Boshoff C: Predicting biomarkers for ovarian cancer using gene-expression microarrays. Br J Cancer 2004, 90:686-692[Medline]
27. Weigelt B, Glas AM, Wessels LF, Witteveen AT, Peterse JL, van’t Veer LJ: Gene expression profiles of primary breast tumors maintained in distant metastases. Proc Natl Acad Sci USA 2003, 100:15901-15905[Abstract/Free Full Text]
28. Lee H, Lin ECK, Liu L, Smith JW: Gene expression profiling of tumor xenografts: in vivo analysis of organ-specific metastasis. Int J Cancer 2003, 107:528-534[Medline]
29. Suzuki N, Withers HR, Koehler MW: Heterogeneity and variability of artificial lung colony-forming ability among clones from mouse fibrosarcoma. Cancer Res 1978, 38:3349-3351[Abstract/Free Full Text]
30. Tarin D, Price JE: Metastatic colonization potential of primary tumour cells in mice. Br J Cancer 1979, 39:740-754[Medline]
31. Neri A, Nicolson GL: Phenotypic drift of metastatic and cell-surface properties of mammary adenocarcinoma cell clones during growth in vitro. Int J Cancer 1981, 28:731-738[Medline]
32. Poste G, Doll J, Fidler IJ: Interactions among clonal subpopulations affect stability of the metastatic phenotype in polyclonal populations of B16 melanoma cells. Proc Natl Acad Sci USA 1981, 78:6226-6230[Abstract/Free Full Text]
33. Poste G: Cellular heterogeneity in malignant neoplasms and the therapy of metastases. Ann NY Acad Sci 1982, 397:34-48[Medline]
34. Wang N, Yu SH, Liener IE, Hebbel RP, Eaton JW, McKhann CF: Characterization of high- and low-metastatic clones derived from a methylcholanthrene-induced murine fibrosarcoma. Cancer Res 1982, 42:1046-1051[Abstract/Free Full Text]

This article has been cited by other articles:

				D. Tarin Comparisons of Metastases in Different Organs: Biological and Clinical Implications Clin. Cancer Res., April 1, 2008; 14(7): 1923 - 1925. [Full Text] [PDF]

				J. DeMasi, M. C. Chao, A. S. Kumar, and P. M. Howley Bovine Papillomavirus E7 Oncoprotein Inhibits Anoikis J. Virol., September 1, 2007; 81(17): 9419 - 9425. [Abstract] [Full Text] [PDF]

				M. Suzuki, E. Mose, C. Galloy, and D. Tarin Osteopontin Gene Expression Determines Spontaneous Metastatic Performance of Orthotopic Human Breast Cancer Xenografts Am. J. Pathol., August 1, 2007; 171(2): 682 - 692. [Abstract] [Full Text] [PDF]

				E. I. Chen, J. Hewel, J. S. Krueger, C. Tiraby, M. R. Weber, A. Kralli, K. Becker, J. R. Yates III, and B. Felding-Habermann Adaptation of Energy Metabolism in Breast Cancer Brain Metastases Cancer Res., February 15, 2007; 67(4): 1472 - 1486. [Abstract] [Full Text] [PDF]

				M. Suzuki, E. S. Mose, V. Montel, and D. Tarin Dormant Cancer Cells Retrieved from Metastasis-Free Organs Regain Tumorigenic and Metastatic Potency Am. J. Pathol., August 1, 2006; 169(2): 673 - 681. [Abstract] [Full Text] [PDF]

				O. Mendez, B. Martin, R. Sanz, R. Aragues, V. Moreno, B. Oliva, V. Stresing, and A. Sierra Underexpression of transcriptional regulators is common in metastatic breast cancer cells overexpressing Bcl-xL Carcinogenesis, June 1, 2006; 27(6): 1169 - 1179. [Abstract] [Full Text] [PDF]

				S. L. Kominsky and N. E. Davidson A "Bone" Fide Predictor of Metastasis? Predicting Breast Cancer Metastasis to Bone J. Clin. Oncol., May 20, 2006; 24(15): 2227 - 2229. [Full Text] [PDF]

				L. Espana, B. Martin, R. Aragues, C. Chiva, B. Oliva, D. Andreu, and A. Sierra Bcl-xL-Mediated Changes in Metabolic Pathways of Breast Cancer Cells: From Survival in the Blood Stream to Organ-Specific Metastasis Am. J. Pathol., October 1, 2005; 167(4): 1125 - 1137. [Abstract] [Full Text] [PDF]

				P. S. Steeg New Insights into the Tumor Metastatic Process Revealed by Gene Expression Profiling Am. J. Pathol., May 1, 2005; 166(5): 1291 - 1294. [Full Text] [PDF]

This Article Abstract Full Text (PDF) Supplemental Material Purchase Article View Shopping Cart Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Montel, V. Articles by Tarin, D.