Identification of genes differentially expressed between benign and osteopontin transformed rat mammary epithelial cells

Background Osteopontin is a secreted, integrin-binding and phosphorylated acidic glycoprotein which has an important role in tumor progression. Findings In this study, we have utilized suppressive subtractive hybridization (SSH) to evaluate OPN regulated gene expression, using the Rama 37 benign non-invasive rat mammary cell line and a subclone, Rama 37-OPN. Rama 37-OPN was produced by stably transfecting Rama 37 with an OPN expression vector and it demonstrates increased malignant properties in vitro. Sequence and expression array analysis of the respective cDNA libraries of over 1600 subtracted cDNA fragments revealed 982 ESTs, 45 novel sequences and 659 known genes. The known up-regulated genes in the Rama 37-OPN library code for proteins with a variety of functions including those involved in metabolism, cell adhesion and migration, signal transduction and in apoptosis. Four of the most differentially expressed genes between the benign and in vitro malignant rat mammary cell lines are tumor protein translationally controlled I (TPTI), aryl hydrocarbon receptor nuclear translocator (ARNT), ataxia telangiectasia mutated (ATM) and RAN GTPase (RAN). The largest difference (ca 10,000 fold) between the less aggressively (MCF-7, ZR-75) and more aggressively malignant (MDA MB 231, MDA MB 435S) human breast cancer cell lines is that due to RAN, the next is that due to osteopontin itself. Conclusion The results suggest that enhanced properties associated with the malignant state in vitro induced by osteopontin may be due to, in part, overexpression of RAN GTPase and these biological results are the subject of a subsequent publication [1].


Background
Metastasis is the major cause of treatment failure in breast cancer patients [1]. The extracellular matrix glycophosphoprotein osteopontin (OPN) is normally secreted by osteoblasts and is utilized as an extracellular adhesion molecule [1]. Osteopontin has also been associated with certain aspects of malignant transformation [2] by enhancing malignant cell attachment contributing to anchorage-independent growth and the migration of tumor cells [3,4]. Moreover, transfection of benign, non-metastatic rat mammary cells with the cDNA for OPN in an expression vector endows the transfectants with the ability to overproduce OPN in vitro and to metastasize in vivo [5]. Furthermore, we and others have shown that OPN overexpression is associated with poor prognosis in human primary breast cancer [6,7]. Recently, OPN has been shown to be the single most powerful prognostic factor in a multivariate analysis against outcome, in a large prospective study of breast cancer patients [8]. Circulating plasma levels of OPN are also higher in metastatic breast cancer patients [9]. However, the precise molecular mechanisms of how OPN regulates metastasis remain unclear.
To obtain a better understanding of this phenomenon and, in particular, the downstream target genes in the OPN signaling pathway, the technique of suppression subtractive hybridization (SSH) has been employed. SSH is an efficient method that uses widely available facilities to identify sequences of known and unknown genes [10]. SSH has been used in the present study to determine differential gene expression between the benign rat mammary cell line Rama 37 (R37) and R37 cells stably transfected with an expression vector for OPN, termed R37-OPN cells.
Using the SSH method combined with reverse Northern hybridization, we have identified genes that are differentially expressed between the benign non-invasive rat mammary cell line and the invasive and malignant metastatic Rama 37-OPN. Some of the differentially expressed genes were then tested further by Real Time PCR for the relative level of their expressed mRNAs in a series of cell lines established from human breast cancer and the mRNA species which changes the most has been identified as RAN GTPase. The biological relevance of these changes are the subject of a subsequent paper [1].

Stable transfection of osteopontin gene
To ascertain if OPN expression was related to invasion and other properties associated with the malignant state in vitro, we raised various stable cell lines from R37 cells. These were stably transfected by empty vector pBK-CMV (R37-pBK-CMV) or by the constitutively-active expression vector OPN-pBK-CMV (R37-OPN cells), as described in Methods. Individual clones of the transfectants were combined for subsequent analysis to generate pooled cell lines. Immunoblot analysis using a monoclonal (MAb) antibody to OPN, which recognizes both human and rat OPN (Materials and Methods), showed that the OPN protein was expressed at a low level in R37 and R37-pBK-CMV cells (Fig. 1A). OPN protein expression was increased 10 fold in R37-OPN compared to R37 and R37-pBK-CMV cell lines (Fig. 1A). In R37, R37-pBK-CMV and R37-OPN cell lines, the MAb to OPN recognized a protein of M r 65,000 (Fig. 1A, Lanes 1-3), consistent with the size of OPN from the original rat cell lines [8,11]. The increase in OPN protein in the R37-OPN pooled cell line suggests that the OPN-pBK-CMV vector consistently overexpresses OPN, the empty vector pBK-CMV having no effect on OPN expression in R37 cells.

Effect of OPN on cellular adhesion, anchorageindependent growth and invasion
Adhesion of cells to fibronectin-coated surfaces was assayed using a dye-based system. R37-OPN cells showed a 9.2 and 8 fold increase in cell adhesion to fibronectincoated dishes, in comparison with R37 and R37-pBK-CMV cells, respectively (Student's t-test, p < 0.01) (Fig.  1B). Colony formation was assayed in soft agar, R37-OPN cells induced a 1.8 and 2.3 fold increase in colony number per plate compared to R37 and R37-pBK-CMV cells, respectively (p < 0.01) (Fig. 1C). The ability of cells to migrate through a reconstituted three dimensional collagen gel (Matrigel) and appear on the underside of a polycarbonate membrane was tested as an assay for cell invasion. The cells on the underside of the membrane were stained, scanned and counted using a digital imaging system described in Methods [12]. Migration of R37-OPN transformants was 913 and 1006 fold greater than the parental R37 and R37-pBK-CMV cells, respectively (p ≤ 0.002) (Fig. 1D). These results suggest that high levels of OPN induce cell adhesion, anchorage-independent growth and invasion in vitro, properties consistent with the malignant metastatic state in vivo [5,13].

cDNA library construction by suppression subtractive hybridization
A cDNA library was constructed from the two cell lines. In the forward SSH library, the benign noninvasive R37 was used as the driver and the malignant invasive R37-OPN cell line as the tester. In the reverse suppression subtractive hybridization library, the benign, noninvasive R37 was used as the tester and the malignant invasive R37-OPN cell line as the driver. The forward and reverse subtracted libraries represent genes that are preferentially up and down-regulated in relation to OPN-mediated transformation in mammary cells. To evaluate the efficiency of cDNA subtraction, we compared the transcript levels of the housekeeping gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) by RT-PCR in subtracted and unsubtracted cDNA libraries from R37 and R37-OPN cells, respectively. Detection of GAPDH sequences for both subtractions required 26 PCR cycles with subtracted cDNA as template, whereas only 10 cycles were required to amplify GAPDH from control cDNAs. Thus the commonly expressed gene GAPDH was significantly depleted from the subtracted cDNA libraries. Overall, 90% of the randomly-picked cloned cDNAs yielded sequence information. We have obtained 291 sequences (comprising 108 libraries yielded a complex of previously identified cancer-associated genes, tumor suppressor genes as well as a set of previously uncharacterized genes with no level of redundancy. These data are consistent with a significant enrichment for invasion-associated and suppressor genes from the R37-OPN and R37 libraries, respectively.

Differential screening of the subtracted libraries
The relative levels of mRNAs corresponding to the sequenced cloned cDNAs, isolated from the subtracted benign, non malignant and invasive, malignant libraries, were estimated in R37 and R37-OPN cells by reverse Northern hybridization, using cDNA produced from mRNA from either the R37 or R37-OPN cells as probes.
The hybridization results were normalized using cDNAs corresponding to mRNAs for GAPDH, which showed similar expression between the R37 and R37-OPN cell lines.

Quantitative Real Time PCR
Four of the highest differentially expressed genes that were also associated with invasion and/or malignancy were chosen from the total number of 1685 for validation by quantitative real-time PCR analysis (QPCR) (see Additional file 5) using the relatively non-invasive MCF-7 and the highly invasive MDA-MB-231 and MDA-MB-435S human breast cancer cell lines [25]. The four genes encode for tumor protein translationally-controlled (TPTI), aryl hydrocarbon receptor nuclear translocator (ARNT), ataxia telangiectasia mutated (ATM) and RAN GTPase (RAN). MCF-7 cell: other cell ratios were calculated using the comparative threshold cycle (Ct) method [26] after normalization to a control housekeeping gene for ribosomal RNA S18. Results of the mean ± standard error from three independent experiments are shown. C Soft agar assay was carried out to assess the ability of stably transfected cell lines to grow in an anchorage independent environment. The colony number was assessed after 5 days.
Results of the mean ± standard error from three independent experiments are shown. D R37, R37-pBK-CMV and R37-OPN were plated on ECM-coated filters (500 μg/ml) in Boyden chambers. The number of cells that migrated through the filter after 48 hrs was determined by staining and scanning using a digital imaging system (Methods). Results of the mean ± standard error from three independent experiments are shown.
fold increase, respectively, in RAN expression compared to MCF-7 cells and they also showed a 9785 and 60 fold increase in levels of OPN (see Additional file 5). The relatively non-invasive breast cancer cell line ZR-75 showed almost identical levels of low expression of RAN and OPN to those of the MCF-7 cell line (see Additional file 5).

Discussion
These results are not a simple artefact of a single clone of cells, since they have been obtained with pools of single cell clones of transfectants. R37 cells thus provide a robust system for the study of OPN-induced changes in the parental R37 cells that lead to the invasive, malignant phenotype in vitro and, in a parallel system, to the metastatic state in vivo [5]. Using SSH technology, we have now identified potential downstream effectors of the OPNinduced signaling network that may lead to invasion and ultimately metastasis. Here we report the recovery of gene fragments representing differentially expressed mRNAs between benign, noninvasive R37 [11] and malignant, invasive R37-OPN cells [11] from two cDNA libraries established after SSH, and compare their sequences with those of known genes (see Additional files 1 &2, 3 &4). The SSH procedures used to establish differences in gene expression between the two cell types have permitted the isolation of genes expressed in high and low-abundance classes. It also leads directly to the production of cloned fragments of expressed genes without the necessity to clone subsequently genes of interest.

Cell culture and production of stable transformant cell lines
In vitro tests for malignancy Assays for cell adhesion, colony formation and invasion through Matrigel in Boyden chambers were carried out as previously described ( [13]; see Additional file 6).

Synthesis of SSH cDNA libraries
Two subtracted cDNA libraries, from R37 and R37-OPN cell lines, were synthesized using the PCR-Select™ cDNA subtraction kit (CLONTECH) (see Additional file 6).

Cloning and sequence analysis of OPN-target genes
The forward and reverse subtracted cDNAs were cloned into pCR2.1-TOPO vectors (Invitrogen) and transformed into competent TOPO 10 cells (Invitrogen) (see Additional file 6).

Expression array screening and analysis
A total of 327 combined cDNA inserts from forward subtracted and reverse subtracted libraries were PCR amplified using NP1R and NP2R primers as previously described [27] (see Additional file 6).

Quantitative Real Time RT-PCR (QPCR)
QPCR was used as an independent method to probe the association of identified differentially expressed genes with that of OPN in different human breast cancer cell lines (see Additional file 6).

Statistical treatment of results
All biological experiments were performed at least 3 times. The mean and standard error were calculated and p values less than 0.05 were considered significant as calculated using the Student's t-test.

Authors' contributions
VVK: Carried out the laboratory molecular biology studies and sequence alignment; PGJ: Participated in data analysis and revising the manuscript; PSR: Provided the parental Rama 37 and revising the manuscript and MKE: conceived, coordinated and designed the study, led the research, data analysis and wrote the entire manuscript. All the authors read and approved the final manuscript.

Additional material
Additional File 1