Cellular processes of v-Src transformation revealed by gene profiling of primary cells - Implications for human cancer
BMC Cancer volume 10, Article number: 41 (2010)
Cell transformation by the Src tyrosine kinase is characterized by extensive changes in gene expression. In this study, we took advantage of several strains of the Rous sarcoma virus (RSV) to characterize the patterns of v-Src-dependent gene expression in two different primary cell types, namely chicken embryo fibroblasts (CEF) and chicken neuroretinal (CNR) cells. We identified a common set of v-Src regulated genes and assessed if their expression is associated with disease-free survival using several independent human tumor data sets.
CEF and CNR cells were infected with transforming, non-transforming, and temperature sensitive mutants of RSV to identify the patterns of gene expression in response to v-Src-transformation. Microarray analysis was used to measure changes in gene expression and to define a common set of v-Src regulated genes (CSR genes) in CEF and CNR cells. A clustering enrichment regime using the CSR genes and two independent breast tumor data-sets was used to identify a 42-gene aggressive tumor gene signature. The aggressive gene signature was tested for its prognostic value by conducting survival analyses on six additional tumor data sets.
The analysis of CEF and CNR cells revealed that cell transformation by v-Src alters the expression of 6% of the protein coding genes of the genome. A common set of 175 v-Src regulated genes (CSR genes) was regulated in both CEF and CNR cells. Within the CSR gene set, a group of 42 v-Src inducible genes was associated with reduced disease- and metastasis-free survival in several independent patient cohorts with breast or lung cancer. Gene classes represented within this group include DNA replication, cell cycle, the DNA damage and stress responses, and blood vessel morphogenesis.
By studying the v-Src-dependent changes in gene expression in two types of primary cells, we identified a set of 42 inducible genes associated with poor prognosis in breast and lung cancer. The identification of these genes provides a set of biomarkers of aggressive tumor behavior and a framework for the study of cancer cells characterized by elevated Src kinase activity.
The v-Src kinase, the product of the Rous sarcoma virus (RSV), has provided a paradigm for the study of signaling pathways and mechanisms of cell transformation by receptor and non-receptor type tyrosine kinases. Its cellular counterpart, c-Src, is a member of a small family of kinases sharing a similar domain organization, overall structure and regulatory mechanism. Members of the Src family of kinases (SFK) contribute to several aspects of the activity of receptor tyrosine kinases including receptor turn-over, reorganization of the cytoskeleton and the initiation of DNA synthesis .
A role for c-Src in the initiation or progression of human cancer has been documented in several studies. Elevated Src kinase activity has been observed in several human cancers and in particular in breast, ovary, lung, bladder, stomach and colon carcinomas . The majority of breast tumors samples (>70%) show elevated Src kinase activity that reflects increased protein levels . While c-Src over-expression is not sufficient to induce cell transformation, c-Src likely cooperates with other tyrosine kinases, such as the EGF receptor, frequently over-expressed in the same tumors. An activating mutation resulting in the deletion of the c-Src C-terminal region adjacent to the negative regulatory tyrosine (Y530) has also been identified in a subset of patients with advanced colon carcinomas . This mutation mimics the oncogenic activation of v-Src, whose C-terminus lacks the C-terminal Src kinase (Csk) phosphorylation site.
Signaling pathways controlling cell proliferation or survival, in particular the Ras and PI3K pathways, have been the subject of intense investigation in v-Src transformed cells [5–7]. More recently, elevated Src kinase activity has been linked to several aspects of tumorigenesis including modification of the tumor micro-environment, vasculogenesis, metastasis and the acquisition of chemoresistance [8, 9]. The mechanisms by which Src controls these properties of tumor cells remain largely unknown. One of the earliest and defining observations of v-Src transformation is the capacity of this oncoprotein to modify the pattern of gene expression. This was revealed by the cloning and characterization of genes aberrantly expressed in v-Src transformed cells, including genes encoding metalloproteinases and chemokines, and the trans-acting factors regulating their expression [10–14]. Investigations based on gene disruption or the use of dominant-negative mutants established the importance of transcription factors such as AP-1, members of the STAT and Ets families, and c-Myc in the proliferation and behavior of v-Src transformed cells [15–22].
Gene profiling studies of transformed cells or tumors characterized by elevated Src kinase activity have documented the changes in gene expression associated with this oncoprotein [23–26]. However, these studies did not assess the transformation-dependence of gene expression . Often they were performed on heterogeneous tumor specimens that limit the sensitivity of the analysis or in immortalized cell lines that may not reveal the full range of Src-dependent changes in gene expression. In this study, we employ a different strategy by comparing the expression profiles of two different primary cell types infected by the Rous sarcoma virus, namely chicken embryo fibroblast (CEF) and chicken neuroretinal cells (CNR). By using temperature sensitive or transformation deficient mutants of v-Src, we define a common set of genes regulated by v-Src. The expression of a group of 42 v-Src inducible genes of the common set was associated with reduced disease-free survival in independent cohorts of patients with breast carcinomas. High expression of this gene set was also observed in primary tumors of patients with reduced bone or lung metastasis-free survival, suggesting that the common program of v-Src transformed cells is enriched for genes associated with a more aggressive tumor phenotype. The definition of this v-Src gene signature provides a set of biomarkers for the identification of tumors with an aggressive phenotype and the framework for the study of cellular properties conferred by high Src kinase activity in human tumors.
Cell culture and RSV strains
Chicken embryo fibroblasts (CEF) were isolated from day 10 embryos while chicken neuroretinal cells (CNR) were isolated from the retina of 7 day old embryos, as described before . CNR cells do not proliferate in the absence of serum or a transforming v-Src kinase and therefore were infected with the wt Schmidt Ruppin-A strain (SR-A) of RSV or with the temperature-sensitive mutant ts NY72-4 RSV. CEF were infected with the RSV strains wt SR-A, ts NY72-4 and td NY315 or with RCASBP-A, a derivative of avian sarcoma viruses lacking a viral Src gene. All viruses belong to the A sub-group of ASV. Characterization of the ts NY72-4 and td NY315 RSV mutants were described before [28, 29]. CNR and CEF were cultured in high glucose Richter's modified MEM medium (HyClone #SH30601) supplemented with 5% cosmic calf serum (HyClone #SH30087.03), 5% tryptose phosphate broth, glutamine, penicillin and streptomycin. All studies were performed with actively dividing cells cultured in medium replenished the day before sample preparation to avoid starvation and acidosis of the transformed cells. Populations of ts NY72-4 RSV-infected CNR cells were expanded at the permissive temperature of 37.5°C while CEF were cultured at the non-permissive temperature of 41.5°C until transferred to the permissive temperature of 37.5°C to induce transformation. For PI3K inhibition, CEF were treated for 8 hours with either a 100 nm dose of wortmannin (Biomol Research Laboratories, Pennsylvania) or a 15 μM dose of LY290042 (Biomol Research Laboratories, Pennsylvania). Animal use was approved by the McMaster University Animal Research Ethics Board (AUP#05-06-26) and conducted in accordance with the guidelines of the Canadian Council on Animal Care.
RNA isolation and northern blotting analysis
Total cell RNA was isolated using Trizol reagent (Invitrogen). Ten μg of total RNA was loaded per well and separated on a 1.2% agarose gel containing 3.75% formaldehyde. RNA was transferred to nylon membrane (Schleicher & Schuell) by capillary transfer, cross-linked, and probed with randomly primed radiolabeled DNA fragments as described before . Signals detected on northern blots were quantified directly by phosphor-imaging using a Storm 820 phosphorimager (Molecular Dynamics) or by scanning autoradiographs on a Umax Astra 1220 U scanner. Images were quantified using ImageJ  and adjusted for loading using GAPDH.
Gene profiling analyses
RNA samples were isolated using Trizol (Invitrogen). RNA quality was assessed by gel electrophoresis and examined with a Bioanalyzer (Agilent). All RNA samples were first analyzed by northern blotting analysis and probed for IL8 and GAPDH expression. RNA samples with a RNA Integrity of less than 9.7 were discarded. Microarray experiments were conducted at the Toronto Centre for Applied Genomics (TCAG) at the Hospital for Sick Children (Toronto, Canada). Biotinylated cRNAs were generated at TCAG and hybridized to Affymetrix Chicken GeneChip arrays using standard Affymetrix protocols (EukGE-WS2v4). GeneChips were scanned using the Affymetrix GeneChip Scanner 3000. Feature intensity was quantified using Command Console software and exported to CEL format.
CEL files were analysed using dChip software version 2007 . Array data were normalized against the median intensity array for each experiment using the invariant set normalization method . Median array intensities for the three experiments were scaled to the same value prior to normalization in order to provide comparable probe intensities for inter-experimental comparison. Expression levels were determined by the model-based expression index method using perfect-match probes only [33, 34]. Two-fold changes in gene expression between experimental conditions within each experiment were determined using log-transformed expression values and statistical significance of expression was determined using one-way ANOVA (RCASBP/NY315/SR-A comparison using RCASBP(A) and NY315 as control groups) or an unpaired two-tailed t-test (NY72-4 CEF and CNR comparisons). Probe-sets whose p-values were greater than 0.05 (pairwise comparisons) or 0.01 (ANOVA) were excluded from further analysis. False discovery rates in all comparisons were estimated to be less than 3% by permutation (10,000 permutations;). For clustering, redundant genes were removed by identifying probe-sets whose Entrez Gene and/or reference sequence IDs were identical. Probe-sets with the lowest p-values were retained. Gene clustering was performed using unsupervised hierarchal clustering by average Euclidean distance . Venn diagrams were drawn using Vennmaster v. 0.35 .
Experimental groups were divided into the following: CEF infected by RCASBP(A), NY315 RSV, or SR-A RSV (experimental group 1); CEF infected by ts NY72-4 RSV and grown either at the permissive or non-permissive temperature (experimental group 2); and CNR cells infected by ts NY72-4 RSV and grown either at the permissive or non-permissive temperature (experimental group 3). A minimum of three biological replicates per experimental condition was used. An additional three RCASBP(A)-infected CEF samples cultured at 37.5°C were used to control for temperature effects (normalized within experimental group 1). Temperature-regulated genes (determined using the same criteria as for the experimental groups) were cross-referenced to all differentially expressed genes from the CEF ts NY72-4 and CNR tsNY72-4 experiments. Genes whose change in magnitude expression was greater than or equal to the same genes found to be differentially expressed in the SR-A/RCASBP(A) experiment were considered to be temperature-regulated and not attributable to v-Src-dependent transformation. This includes 13 genes from the CEF ts NY72-4 and 16 genes from the CNR ts NY72-4 experiments respectively (Additional File 1). Arrays in each experimental group were normalized within said experimental group to prevent skewing of data points from tissue-effects as determined by M-A plots for individual arrays. Array data are accessible through the Gene Expression Omnibus http://www.ncbi.nlm.nih.gov/geo/ using the accession number GSE14489.
Gene ontology and pathway analyses
Gene ontology (GO) analyses were conducted using the DAVID 2.0 bioinformatics resources [38, 39] functional annotation tool using an EASE score of 0.05 and a minimum of two genes per term. Pathway analyses were carried out using Pathway Express software available in the Onto-tools package . Probe-sets associated with significantly differentially expressed genes from each experimental group were converted to the orthologous probe sets from the Affymetrix human U133 Plus 2.0 array. Non-redundant orthologous probe-sets and associated fold change differences (linear-scale) were uploaded to Pathway Express for use in pathway analysis. Pathway analysis was conducted using default settings and only pathways whose corrected gamma-p-value was less than or equal to 0.05 were kept.
Tumor data analyses
Tumor expression data sets were obtained from the Gene Expression Omnibus at the National Center for Biotechnological Information. Breast tumor data from Pawitan et al.  and Ivshina et al.  were used as training data sets (Figure 5, panels Ai and Aii respectively). Breast tumor data from Minn et al. , colon tumor data from the International Genomics Consortium Expression Project for Oncology http://www.intgen.org/, and breast, lung and ovary tumor data from Bild et al.  were used as test sets. All clinical data were obtained from the original publications. Orthologous probe-sets (Hu95av2, U133, U133 Plus 2.0 Affymetrix arrays) for the CSR gene set were obtained from NetAffyx (Affymetrix).
To ascertain if up-regulated genes in the CSR set could be used as markers of poor prognosis in human cancer, orthologous probe-sets corresponding to the up-regulated genes in the CSR set were obtained for the Affymetrix platforms used in the training data sets. Only the up-regulated orthologous CSR probe-sets were used for clustering of the Pawitan and Ivshina training sets. Unsupervised gene and sample clustering (Pearson correlation distance metric ) was performed on the two training data in order to identify clusters of tumors with similar gene expression profiles. Each tumor cluster was subjected to survival analysis in order to identify tumor clusters associated with poor survival. Since the objective was to identify up-regulated CSR genes associated with poor prognosis, the genes whose expression was statistically significantly higher in the poor-prognosis tumor set versus the other three sets were chosen as candidates for the aggressive gene signature (t-test, single-tailed, p < 0.05). The common set of 42 genes identified from the two training sets defined the aggressive gene signature.
For the test data, the sets of tumors expressing high levels of the aggressive signature were defined as the upper quartile of tumors from any given test set with the highest average gene-wise mean-centered expression of the aggressive signature set. Survival analyses were carried out using the Kaplan-Meier product limit estimator method  and statistical significance was evaluated using the log-rank test.
Immunofluorescence and western blotting
CEFs were grown to sub-confluence on glass slides, rinsed with PBS twice and fixed in 1% formaldehyde/PBS for ten minutes. Cells were washed three times for ten minutes in PBS and incubated with anti-CD44-FITC antibody (Abcam, cat. # ab24907) for 2 hours at room temperature and washed three times for ten minutes in PBS. Cells were visualized on a Nikon Eclipse TE2000-U inverted microscope. Western blotting was performed as described before  using the following antibodies: anti-GSK3-β and phospho-S9-GSK3-β (Cell Signaling #9332 and 9336S), anti-heme oxygenase 1 (Stressgen, #SPA-896), anti-PKB (New England Biolabs #9272), anti-phospho-PKB (Cell Signaling #9271), anti-Erk1 (Santa Cruz Biotechnology #sc-94), and anti-Twist1 (Santa Cruz Biotechnology #sc-15393). Autoradiographs were scanned on a Umax Astra 1220 U scanner, quantified using ImageJ and corrected for loading using Erk1 .
Preparation and cloning of DNA fragments
RNA was isolated from 107 CEFs using Trizol reagent (Invitrogen) as per manufacturer's directions. One μg of RNA was reverse transcribed (Invitrogen SuperScript RT) and treated with RNAse A (Invitrogen). Double-stranded DNA probes for northern blotting analyses were PCR amplified using Promega GoTaq and the appropriate forward/reverse primers (Table 1). DNA fragments were gel-purified and cloned into pCR2.1-TA (Invitrogen) vector as per manufacturer's directions. All cloned products were verified by sequencing at the MOBIX sequencing facility, McMaster University. CEF4/IL8 and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) probes were generated as described previously .
Identification of Transformation-Regulated genes in v-Src transformed CEF
To identify genes regulated in a transformation-dependent manner, we characterized the gene expression profiles of CEF transformed by the wt RSV strain Schmidt Ruppin - group A (SR-A) or infected with the transformation-deficient viruses NY 315 RSV or RCASBP(A). NY315 RSV is a group A virus encoding a catalytically active deletion mutant of the v-Src kinase lacking amino acids 2-14 of the SH4 domain. As a result of this deletion, the v-Src kinase of NY315 RSV is not myristoylated, does not associate with the plasma membrane and is non-transforming . RCASBP(A) is a group A replication competent virus of the avian sarcoma virus family lacking an oncogene . Total cellular RNA was isolated and analyzed on the Affymetrix chicken GeneChip, an array comprising 37,703 probe sets representing 19,881 UniGene clusters. A total of 3254 probe sets, corresponding to 2904 unique genes were differentially expressed by two fold or greater in cells infected by either NY315 or SR-A RSV (Additional File 2). Within regulated sequences, 2455, 1691 and 730 transcripts were differentially expressed in pair-wise comparisons between SR-A and RCASBP(A), NY315 and RCASBP(A), and SR-A and NY315 infected CEF, respectively.
Genes regulated in a transformation-dependent manner were identified by one-way ANOVA (p≤0.01) using RCASBP(A) and NY315 as control groups and defined by probe sets with differences in expression greater or equal to two-fold between any two conditions. Unsupervised hierarchal clustering shows that transformation responsive genes cluster to one of two groups, either up-regulated or down-regulated in SR-A RSV transformed CEF but not in the other two conditions (Figure 1). A total of 1095 transcripts corresponding to 418 up-regulated and 535 down-regulated genes were identified in this analysis, defining a set of genes regulated in a transformation-dependent fashion. Approximately 15% of these genes were differentially expressed by five fold or more. We refer to this set of genes as the Transformation-Regulated genes (TR genes) of v-Src transformed CEF (Additional File 3).
Differentially expressed genes were also characterized in CEF infected with the temperature sensitive (ts) mutant NY72-4 RSV. A pair-wise comparison of CEF infected with NY72-4 revealed that 568 transcripts encoded by 477 unique genes (261 up- and 216 down-regulated) were differentially expressed at the permissive versus non-permissive temperatures (Figure 2 and Additional File 4). Cross-referencing of the gene lists provided by the SR-A/NY315/RCASBP(A) and NY72-4 RSV analyses revealed a set of 145 non-redundant genes (or 199 transcripts) corresponding to 84 down- and 61 up-regulated genes associated with cell transformation (Additional File 5).
To control for the effect of temperature change (cold shock), we compared the profile of RCASBP(A) infected CEF at the permissive temperature of 37.5°C and non-permissive temperature of 41.5°C (Additional File 1). Few genes were differentially expressed by a factor of three or more at the two different temperatures. However, some of the highly v-Src responsive genes (HMOX1, UPP1, Aquaporin 1, for instance) were also regulated by temperature change. Genes affected by temperature and wt v-Src, as determined by the analysis of control and SR-A RSV-transformed CEF at 41.5°C, indicate a sub-class of genes that are both stress-responsive and v-Src-regulated.
Identification of v-Src regulated genes in chicken neuroretinal (CNR) cells
In vertebrates, the c-Src kinase is widely expressed and detected in most tissues. In addition, v-Src is capable of transforming a variety of cell types in vitro [48–54]. To compare the biological response of different cell types transformed by v-Src, we characterized the patterns of gene expression of chicken neuroretinal (CNR) cells infected with the temperature-sensitive mutant ts NY72-4 RSV. Unlike CEF, neuroretinal cells do not proliferate in vitro unless transformed by an oncoprotein such as v-Src and cultured in serum-containing medium . Therefore, ts NY72-4 RSV infected CNR are quiescent at the non-permissive temperature of 41.5°C but are actively dividing and transformed at the permissive temperature of 37.5°C. Gene profiling analyses of CNR cells identified 1062 transcripts, consisting of 485 up- and 577 down-regulated RNA species, with a two-fold or greater difference in gene expression between the two temperatures (Additional File 6). A significant fraction of the v-Src regulated genes identified in ts NY72-4 RSV (31%) and SR-A/NY315/RCASBP(A) infected CEF (21%) overlapped with the corresponding set of genes in CNR cells (Figure 2). The total number of transcripts (probe sets) and genes differentially expressed in the three systems of v-Src transformation is provided in Table 2. Together, the total genes regulated in these three sets represent approximately 6% of the protein-coding genes in the chicken genome. When all three systems are compared, a group of 84 common genes are identified as v-Src regulated (Figure 2; Additional File 7).
A second group of 91 genes was regulated in ts NY72-4 RSV infected CEF and CNR cells but not in SR-A RSV transformed CEF (Additional File 8). Up-regulated genes in this list encode several important regulators of cell proliferation or behavior including ornithine decarboxylase (ODC), osteopontin, hyaluran synthase 2 (HAS2), cyclin A and cyclin E2. In contrast, the cyclin kinase inhibitor p27Kip1 was down-regulated in NY72-4 RSV transformed CEF and CNR cells. Since CNR cells are amplified and cultured at the permissive temperature for extensive periods of time, it is unlikely that the regulation of these genes represents a transient effect of ts v-Src activation. A more likely explanation for this discrepancy is that differences in the structure of the v-Src kinase encoded by NY72-4 and SR-A RSV account for this specificity in gene regulation. It is also likely that the list of genes defined above (Additional File 7) is an under-representation of the gene cohort controlled by v-Src in CEF and CNR cells. A more accurate description of this class of genes can be obtained by combining the genes described in Additional File 7 and the genes regulated uniquely in NY72-4 RSV transformed CEF and CNR cells (Additional File 8). This combined list would then consist of a total of 175 genes that we define as the Common set of v-Src Regulated genes or CSR genes.
Validation and characterization of v-Src regulated genes
The expression of v-Src regulated genes was confirmed by northern blotting analyses. In this study, the expression of genes identified by gene profiling was compared to that of previously characterized markers of Src transformation including IL8/CEF-4, CD44 and Nov [10, 14, 56–58]. For all up and down-regulated genes selected, little difference was observed between non-transformed CEF infected with the replication competent virus RCASBP(A) and CEF infected with the transformation-deficient virus NY315 RSV (Figure 3A). As reported before, the kinetics of IL8/CEF-4 RNA accumulation was biphasic and characterized by the rapid accumulation of IL8/CEF-4 (within 1 hr) followed by a temporary reduction of RNA level between 4 and 8 hrs of temperature shift . Other RNA species (VIP, CD44, integrin-α-6) accumulated more slowly i.e. between 2 and 4 hrs following ts v-Src activation and did not show the transient down-regulation observed for the IL8/CEF-4 mRNA. IL8/CEF-4 and CD44 were both induced in cells treated with cycloheximide indicating that the different kinetics of RNA accumulation do not necessarily reflect a different requirement for de novo protein synthesis (data not shown).
Cell heterogeneity was assessed by looking at the surface expression of CD44 by immunofluorescence (Additional File 9). A strong CD44 signal was observed in SR-A RSV transformed CEF but was absent in cells infected with RCASBP(A) and NY315 RSV. All cells expressing the ts mutant of v-Src (NY72-4 RSV) were also positive for CD44 expression at the permissive temperature. Low but detectable signals for CD44 expression were also observed in ts NY72-4 RSV-infected cells at the non-permissive temperature. This may be indicative of a certain degree of leakiness of the ts v-Src kinase, an observation accounting for the lower number of v-Src regulated genes identified in this system.
Characterization of the common set of v-Src regulated genes in CEF and CNR cells
Differences in the profiles of v-Src regulated genes in CEF and CNR cells reflect in part the proliferative and differentiation state of these cells since CNR cells do not proliferate when non-transformed and are partially committed to differentiation [55, 60]. However, a common program of gene expression was also evident in the profiles of v-Src transformed CEF and CNR cells. Genes such as IL8/CEF4, VIP, HMOX1, PLCPI and UPP1, were all activated in v-Src transformed CEF and CNR cells. HMOX1 also provided an example of a gene with partial activation in CEF infected with NY315 RSV, both at the RNA and protein level (Figure 3C &3E).
A comparative analysis of v-Src regulated mRNA species examined by northern blotting and microarray profiling was performed for all mRNA species included in Figure 3. This study revealed a strong correlation (Spearman ρ of 0.83, p < 0.0001) between northern and microarray gene expression estimates. A regression coefficient of 0.91 indicates a nearly 1:1 correspondence between northern blotting and microarray gene expression values (Additional File 10, Additional File 11). Collectively the results of the immunofluorescence, northern, and western blotting analyses validated the expression of a sub-set of genes identified in the gene profiling studies and illustrated the existence of a common program of gene expression controlled by v-Src in different cell types.
Signaling pathways and biological processes regulated during v-Src transformation
Pathways potentially regulated by v-Src-induced changes in gene expression were identified by Pathway Express software . Pathway Express estimates whether genes differentially expressed in a given data set impact a given pathway by incorporating fold-change of gene expression, where those genes are in the hierarchy of a pathway and whether those genes are found to be statistically overrepresented in that pathway. Using this methodology, "ECM-receptor interaction", "focal adhesion" and "phosphatidylinositol signaling" were identified as pathways most significantly altered in all three data sets (corrected γ p-value ≤ 0.05; Table 3). TGFβ signaling was also identified in CEF but not CNR cells. Several pathways were also identified in a single data set and, in particular, in the large cohort of the Transformation-Regulated gene set. Thus, "leukocyte transendothelial migration", "complement and coagulation cascades", "cell adhesion molecules", "MAPK" and "Toll-like receptor signaling" pathways were all identified in this data set (Additional File 12).
Pathway Express stresses the regulatory nature of a gene product. This is best illustrated by the genes encoding the pro-survival regulatory and catalytic sub-units of PI3K, which are part of the "focal adhesion", "PI3K" and "leukocyte transendothelial migration" pathways and are hallmarks of several cancer types (ex. "small cell lung cancer"). The up-regulation of both the catalytic and regulatory subunits of PI3K in SR-A RSV transformed CEF implies that the pathway is strongly activated in these cells (Additional File 3). This was confirmed by looking at the phosphorylation status of PKB/Akt and GSK3β, two downstream targets of the PI3K pathway. High levels of phospho-PKB and phospho-GSK3β were detected in SR-A RSV transformed CEF but not in cells infected with RCASBP(A) or NY315 RSV, even when these cells were actively cycling (Additional File 13). The signal was reduced in cells treated with the specific inhibitor LY294002, indicating that phosphorylation of these proteins was dependent on the PI3K pathway.
Other important regulatory factors are included in the list of Transformation-Regulated genes. An example is provided by Twist1, a bHLH transcription factor promoting metastasis and capable of inducing the epithelial-to-mesenchymal transition (Additional Files 3 &4;). The up-regulation of Twist1, observed as a 26 kDa doublet by western blotting, was confirmed in SR-A RSV transformed CEF (Figure 4).
The Common Set of v-Src Regulated (CSR) genes includes a subset of genes associated with reduced disease-free survival in human cancer
To test whether v-Src-induced changes in gene expression correlate with reduced disease- or metastasis-free survival in human cancers, we undertook a clustering-based enrichment regime to identify up-regulated genes of the CSR set associated with this phenotype. To this end, two independent breast carcinoma datasets with clinical information were selected and used as training sets. Mammary tumor sets were used because of the abundance of good quality array data obtained on more recent Affymetrix platforms and because a large proportion of breast carcinomas (>70%) exhibit high Src activity . It was also reasoned that CSR genes were more likely to be regulated in tumors because their expression is controlled by v-Src in different cell types (i.e. CEF and CNR cells). The aggressive gene signature, defined by clustering analysis, was then tested against additional, independent tumor datasets to assess if it could be used to predict reduced disease-free survival.
Probe-sets corresponding to orthologous genes from the up-regulated CSR set and tumor samples from the two independent breast carcinoma datasets [41, 42] were clustered to identify samples containing unique up-regulated gene clusters associated with poor disease-free survival. Survival analyses of these tumor datasets revealed tumor sample clusters associated with higher disease relapse (Figure 5A). Figure 5Ai shows that sample cluster 3 has poorest estimated mean disease-free survival (5.6 years) compared to the other three clusters (7.6, 7.8, and 6.8 years; clusters 1, 2 and 4 respectively; P < 0.001). In Figure 5A, panel ii, cluster 4 has the poorest estimated mean survival compared to clusters 1, 2 and 3 (7.8, 10.1, 10.7 and 8.4 years respectively; P < 0.01). CSR genes that were up-regulated in the low-survival cluster in relation to the other tumor clusters of each tumor data set were cross-referenced to each other to generate an aggressive tumor gene signature (Table 4, see methods). This set contains 42 unique genes overlapping with an 80-87% concordance with respect to the genes found differentially up-regulated in the two tumor data training sets. This striking overlap represents over 50% of the 80 up-regulated CSR genes (Additional Files 7 &8) indicating that the common program of v-Src gene expression is enriched for genes associated with an aggressive tumor phenotype. Functional classification of these genes showed greatest enrichment for gene ontology biological processes (GOBP) such as "Cell Cycle" ("M-phase", "Mitosis", "Cell Division"), "DNA Metabolic Process" ("Response to DNA damage stimulus"), "Response to Stress", "Cell Proliferation" and "DNA Replication". (Additional File 14)
To assess whether the aggressive tumor gene signature could be used to predict poor disease-free survival, six independent tumor datasets were interrogated to determine if high expression of the 42 genes defining this signature was associated with a poor outcome. By comparing the top 25% of tumors expressing the highest levels of the aggressive gene signature against the bottom 75%, four data sets out of six showed that patients with tumors characterized by high expression of the 42 gene cohort have reduced disease-free or metastasis-free survival (Figure 5B). Reduced relapse-free survival was seen in the breast and lung carcinoma datasets while reduced lung and bone metastasis-free survival was observed in two other breast tumor datasets.
Analysis of tumors expressing the highest levels of the aggressive signature genes in the breast tumor set from Minn et al.,  indicated that the patients from whom the tumors were derived had reduced disease-free survival with respect to lung metastasis (mean estimated survival 5.4 years versus 10.1 years, P < 0.001; Figure 5Bi). Similarly, using bone metastasis data from the same study, the high aggressive signature expressing tumors were associated with reduced metastasis-free survival in patients (mean estimated survival 6.4 years versus 9.3 years, P < 0.05; Figure 5Bii). Analysis of a second breast tumor data set  also shows that patients whose tumors express high levels of the aggressive signature genes have a reduced disease-free survival (mean estimated survival 71 months versus 106 months, P < 0.02; Figure 5Biii). Lung tumor data from the same study indicates a similar trend where patients with tumors expressing high levels of the aggressive signature have reduced disease-free survival (mean estimated survival 29 months versus 54 months, P < 0.001, Figure 5Biv). Analysis of ovarian tumor data from the same study however did not yield any statistically significant difference (data not shown). Since survival data was not available for the colon carcinoma dataset, we asked whether high expression of the aggressive tumor gene signature could identify patients with distal metastases. Analysis of the relative proportions of metastasis in patients exhibiting high expression of the aggressive signature to those that did not, failed to show any statistically significant difference (data not shown). Taken together, these data show that a subset of the v-Src-regulated CSR genes correlates with poor prognosis in patients with breast tumors, exhibiting some specificity for metastasis to lung and bone.
Definition of a gene signature for v-Src transformation of primary cells
Previous gene profiling studies of v-Src transformed cells have relied principally on the use of immortalized rodent cell lines that may not reflect the full spectrum of biological processes altered during v-Src transformation [24, 25]. In this study, we used transformation-deficient (NY315 RSV) and temperature sensitive mutants of the Rous sarcoma virus (NY72-4 RSV) to identify gene products regulated in a transformation-dependent manner in primary cells. Further, all experiments were performed in conditions promoting asynchronous cell division and limiting acidosis. A separate study of RSV transformed CEF was reported previously . However, most of the genes identified in our studies (80%) were not included in the gene set described by these investigators. Moreover, the results of this study differ from ours in some important aspects. First, Masker and co-workers did not assess the transformation-dependence of differentially expressed genes since their analysis was restricted to a pairwise comparison of SR-A RSV transformed CEF and CEF infected with a replication competent virus devoid of any oncogene. Second, several established markers of v-Src transformation such as IL8/CEF-4 and CD44 [10, 14, 56], were poorly regulated in their experimental conditions and did not pass the significance criteria of the analysis. Third, Masker and co-workers identified some transcription factors such as c-Jun and c-Myc, as being up-regulated by v-Src while we did not. We previously established that early immediate genes such as c-Myc are activated by v-Src in quiescent cells but are not differentially expressed when actively dividing cells are compared . Therefore, some of these differences may reflect the experimental conditions of the studies and, in particular, the conditions of cell culture.
A study performed with ts NY72-4 v-Src identified a limited number of v-Src regulated genes in NIH 3T3 fibroblasts . Key transcriptional regulators identified by these authors, such as Hmga2 (high mobility group AT-hook2) and Id4 (inhibitor of DNA binding 4), are also included in the list of v-Src-regulated genes in CNR cells (Additional File 6). In addition, Liu and co-workers described the regulation of important regulators of the cell cycle and angiogenesis by SSeCKS, a metastasis-suppressor gene down-regulated by v-Src . Interestingly, tumor suppressor genes such as DKK3 and gamma FBPα, the avian homolog of HIC-1 (hypermethylated in cancer 1), were repressed markedly by v-Src in CEF. Thus, the abrupt down-regulation of tumor-suppressor genes may be required for v-Src-dependent cell transformation (Figure 3 and Additional File 3; ).
The common set of v-Src regulated genes (CSR genes) is enriched for genes associated with an aggressive tumor phenotype
Since the c-Src kinase is expressed ubiquitously in most vertebrate tissues, we reasoned that genes regulated by v-Src in multiple cell types would better represent targets of this kinase in tumorigenesis. Using CEF and CNR cells, we thus described a common program of gene expression consisting of 175 genes regulated by v-Src in both cell types. This program, referred to as the Common set of v-Src Regulated genes or CSR genes, was then used to interrogate several independent tumor data sets with corresponding clinical information. Beginning with data from two independent cohorts of breast carcinomas as training sets, we performed hierarchal clustering analyses to identify a group of v-Src inducible genes highly expressed in tumors associated with reduced disease-free survival (Table 4). Significantly, the majority (80-87%) of the 42 genes included in this group were identified in both independent tumor data sets. Since the CSR gene cohort consists of a total of 80 up-regulated genes, the 42 genes identified in this analysis suggest that CSR genes are highly enriched for genes associated with a more aggressive tumor phenotype. To confirm this finding, we then interrogated several additional tumor data sets including some used previously by other investigators [43, 44]. This analysis confirmed the expression of the v-Src inducible gene set in primary tumors of patients with reduced disease-free or metastasis-free survival in four of six tumor data sets, including breast and lung carcinomas (Figure 5). In contrast, no correlation was observed for ovarian and colon carcinomas suggesting tissue specificity of the v-Src aggressive tumor gene signature. A recent report by Zhang and co-investigators shed some light on the role of Src in bone metastasis of breast cancer cells . Knockdown of c-Src in the highly metastatic human BoM-1833 cell line impaired the proliferation of these cells in the bone micro-environment but had no effect on lymph node metastasis in a mouse model. Moreover, these authors provided evidence that Src enhanced the survival of the metastatic cancer cell line by mediating the activation of PKB/Akt in response to CXCR4 activation .
Prognostic gene signatures for breast cancer have been described previously by several groups [62–66]. A comparative analysis of the 70-gene, 76-gene and the Gene expression Grade Index (GGI) signatures revealed that they have similar prognostic ability despite limited overlap in gene identity . Recently, Yu et al. showed that a common set of biological processes is enriched among these gene signatures even when defined by different genes. For instance, "Cell Cycle" is enriched in the 70-gene, 76-gene and the 97-gene GGI signatures even though cyclin E2 (CCNE2) is the only gene of this pathway common to all three signatures [63, 65–67]. Interestingly, CCNE2, the only gene found in more than two signatures  was also identified in our aggressive gene signature along with EXO1 and KIF11, found in two other studies [62, 63] (Table 4). Functional classification of our 42-gene aggressive signature shows that three of the 16 enriched GOBP terms ("Cell Cycle", "Mitosis" and "DNA replication") are identical to three of the 36 core GOBP terms described by Yu et al. Twelve remaining terms from our aggressive signature map to parent terms of one or more of the 36 core terms while two, "blood vessel morphogenesis" and "cell proliferation", are unique to our gene set (Additional File 14). This similarity in enriched processes shared with other breast cancer signatures reinforces the hypothesis that there is a common aggressive tumor "pathway-signature" .
Biological processes defined by the CSR gene signature - DNA damage and the stress response
Several genes of the DNA damage and stress responses are included in the 42 genes of the aggressive tumor signature. DNA damage resulting from the production of reactive oxygen species (ROS) or DNA replication stress has been observed in cells transformed by several oncogenes [68–72]. However, evidence of DNA damage, as indicated by the levels of γ H2AX, is presently lacking in v-Src transformed CEF (our unpublished results). Therefore, the pathways responsible for the production of ROS or replicative stress may be attenuated in these cells. Src-transformed cells are known to be more resistant to conditions of oxidative stress and the induction of genes such as HMOX1 or UPP1 enhance cell survival in conditions of hypoxia or glucose depletion i.e. in conditions promoting the production of ROS [73–77]. Therefore, genes of the aggressive tumor signature may limit the effects of oncogenic stress and provide a survival advantage to the cell.
Several genes identified in this study were also regulated in response to temperature change (cold shock). Since temperature change had a modest effect on the activation of these genes with induction levels of 3-fold or less (Additional File 1), and since all v-Src regulated genes affected by temperature change were also identified in the RCASBP/NY315/SR-A RSV analysis at 41.5°C, it is unlikely that these genes represent a class of genes regulated uniquely as a result of temperature change. Genes like HMOX1, Aquaporin 1 and the metallothioneins MT2A and 3 are induced in response to a wide variety of stress conditions and are widely regarded as general stress response genes. Since HMOX1, Aquaporin 1 and UPP1 are activated by temperature change and are included in the Common Set of v-Src Regulated (Additional Files 7 &8), they may be part of a general and previously unrecognized stress-response program controlled by v-Src.
Biological processes affected by v-Src transformation: Motility and Invasiveness
The large cohort of Transformation-Regulated genes (TR genes; Additional File 3) includes several genes that may contribute to a more aggressive tumor phenotype. Src-transformed CEFs are highly motile and characterized by structures mediating cell migration called podosomes [78, 79]. Like the related invadopodia described in other cancer cells, the podosomes are actin-rich structures closely associated with adhesion molecules and ECM-degrading enzymes such as the matrix metalloproteinases MMP-1, MMP-2 or MMP-9 . In CEF, v-Src transformation stimulated the expression of four critical regulators of actin polymerization and podosome formation, namely N-WASP, cortactin, gelsolin and the p41 subunit of the actin related protein complex 2/3 (Arp 2/3 subunit 1b-p41; Additional File 3[81–84]). Two of these factors (Arp 2/3 subunit 1b-p41 and cortactin) were also activated in CNR cells indicating that v-Src controls the expression of these genes in diverse cell types (Additional File 6). The gene encoding the p41 subunit of the Arp2/3 complex is frequently amplified in pancreatic cancer and is a regulator of the motility and invasion of pancreatic cell lines . In contrast, caldesmon and transgelin were repressed by v-Src in CEF. Transgelin is a multi-functional protein with roles independent of podosome formation. In particular, transgelin is a negative regulator of MMP-9 expression and a suspected tumor suppressor [86, 87]. Independent studies have shown that caldesmon, an actin filament cross-linker, antagonizes the action of the Arp 2/3 complex and is a negative regulator of podosome formation and invasion in transformed cells [88–90]. Therefore, v-Src transformation is characterized by changes in gene expression promoting the dynamic remodeling of the actin cytoskeleton and the assembly of podosomes. It is also significant that wt v-Src induces the expression of the p110δ catalytic subunit, a PI3K isoform required for cell migration in macrophages and breast cancer cells (TR genes; Additional File 3; [91, 92]).
Biological processes affected by v-Src transformation: Epithelial-to-mesenchymal transition and drug resistance
Several studies have linked the expression of v-Src and other Src family kinases to resistance of a variety of chemotherapeutic agents including cisplatin, geftinib, paclitaxel, oxaliplatin and tamoxifen [93–97]. Since chemoresistance is often associated with increased cell motility and invasiveness, it has been suggested that v-Src controls these activities by inducing the epithelial-to-mesenchymal transition (EMT) in tumor cells . In support of this notion, Sham and co-workers reported recently that the up-regulation of the basic-helix-loop-helix factor Twist1 by NF-kB increases chemoresistance of PC3 prostate cancer cells treated with daunorubicin or cisplatin . In a separate report, Cheng et al. identified Twist1 and its target Akt2/PKBβ as factors contributing to the metastatic potential of highly invasive breast carcinoma cell lines. In the same study, these investigators showed that Twist1 and Akt2/PKBβ determine the resistance of these cells to paclitaxel . Twist1 is a potent inducer of the EMT and a member of the Transformation-Regulated gene cohort identified in this study (Figure 4 and Additional File 3). Since forced expression of Twist1 is sufficient to induce the EMT in mammary epithelial cells, the regulation of this factor provides a mechanism by which v-Src may induce the EMT and enhance chemoresistance . Twist1 is not the only marker of the EMT identified in this study since N-cadherin mRNA expression was also up-regulated in SR-A transformed CEF (TR genes; Additional File 3).
Elevated Src kinase activity has been described in several unrelated human tumors and in cells derived from these tumors . In human cancers, high Src activity correlates with progression to a more malignant phenotype and the increased metastatic potential of tumor cells [4, 100, 101]. In this study, we define a signature of 42 v-Src inducible genes whose expression is associated with reduced disease-free survival in breast and lung cancer patients from several independent studies. In one dataset, this 42 gene signature was also associated with reduced bone and lung metastasis-free survival. Several genes of the larger cohort of Transformation-Regulated Genes (TR genes) have also been associated with features of aggressive tumors such as invasiveness and chemoresistance. These genes provide a set of biomarkers and candidate therapeutic targets for the treatment of patients with tumors characterized by more aggressive behavior. The functional characterization of these genes represents a roadmap for the study of tumor cells characterized by elevated Src kinase activity.
Thomas SM, Brugge JS: Cellular functions regulated by Src family kinases. Annu Rev Cell Dev Biol. 1997, 13: 513-609. 10.1146/annurev.cellbio.13.1.513.
Irby RB, Yeatman TJ: Role of Src expression and activation in human cancer. Oncogene. 2000, 19 (49): 5636-5642. 10.1038/sj.onc.1203912.
Biscardi JS, Ishizawar RC, Silva CM, Parsons SJ: Tyrosine kinase signalling in breast cancer: epidermal growth factor receptor and c-Src interactions in breast cancer. Breast Cancer Res. 2000, 2 (3): 203-210. 10.1186/bcr55.
Irby RB, Mao W, Coppola D, Kang J, Loubeau JM, Trudeau W, Karl R, Fujita DJ, Jove R, Yeatman TJ: Activating SRC mutation in a subset of advanced human colon cancers [In Process Citation]. Nat Genet. 1999, 21 (2): 187-190. 10.1038/5971.
Sugimoto Y, Erikson RL: Phosphatidylinositol kinase activities in normal and Rous sarcoma virus-transformed cells. Mol Cell Biol. 1985, 5 (11): 3194-3198.
Aftab DT, Kwan J, Martin GS: Ras-independent transformation by v-Src. Proc Natl Acad Sci USA. 1997, 94 (7): 3028-3033. 10.1073/pnas.94.7.3028.
Webb BL, Jimenez E, Martin GS: v-Src generates a p53-independent apoptotic signal. Mol Cell Biol. 2000, 20 (24): 9271-9280. 10.1128/MCB.20.24.9271-9280.2000.
Shah AN, Gallick GE: Src, chemoresistance and epithelial to mesenchymal transition: are they related?. Anticancer Drugs. 2007, 18 (4): 371-375. 10.1097/CAD.0b013e32801265d7.
Zhang XH, Wang Q, Gerald W, Hudis CA, Norton L, Smid M, Foekens JA, Massague J: Latent bone metastasis in breast cancer tied to Src-dependent survival signals. Cancer Cell. 2009, 16 (1): 67-78. 10.1016/j.ccr.2009.05.017.
Sugano S, Stoeckle MY, Hanafusa H: Transformation by Rous sarcoma virus induces a novel gene with homology to a mitogenic platelet protein. Cell. 1987, 49 (3): 321-328. 10.1016/0092-8674(87)90284-4.
Simmons DL, Levy DB, Yannoni Y, Erikson RL: Identification of a phorbol ester-repressible v-src-inducible gene. Proc Natl Acad Sci USA. 1989, 86 (4): 1178-1182. 10.1073/pnas.86.4.1178.
Matrisian LM, Glaichenhaus N, Gesnel MC, Breathnach R: Epidermal growth factor and oncogenes induce transcription of the same cellular mRNA in rat fibroblasts. EMBO J. 1985, 4 (6): 1435-1440.
Gillet G, Guerin M, Trembleau A, Brun G: A Bcl-2-related gene is activated in avian cells transformed by the Rous sarcoma virus. EMBO J. 1995, 14 (7): 1372-1381.
Bedard PA, Alcorta D, Simmons DL, Luk KC, Erikson RL: Constitutive expression of a gene encoding a polypeptide homologous to biologically active human platelet protein in Rous sarcoma virus-transformed fibroblasts. Proc Natl Acad Sci USA. 1987, 84 (19): 6715-6719. 10.1073/pnas.84.19.6715.
Wasylyk C, Maira SM, Sobieszczuk P, Wasylyk B: Reversion of Ras transformed cells by Ets transdominant mutants. Oncogene. 1994, 9 (12): 3665-3673.
Turkson J, Bowman T, Garcia R, Caldenhoven E, De Groot RP, Jove R: Stat3 activation by Src induces specific gene regulation and is required for cell transformation. Mol Cell Biol. 1998, 18 (5): 2545-2552.
Suzuki T, Murakami M, Onai N, Fukuda E, Hashimoto Y, Sonobe MH, Kameda T, Ichinose M, Miki K, Iba H: Analysis of AP-1 function in cellular transformation pathways. J Virol. 1994, 68 (6): 3527-3535.
Lloyd A, Yancheva N, Wasylyk B: Transformation suppressor activity of a Jun transcription factor lacking its activation domain. Nature. 1991, 352 (6336): 635-638. 10.1038/352635a0.
Johnson R, Spiegelman B, Hanahan D, Wisdom R: Cellular transformation and malignancy induced by ras require c-jun. Mol Cell Biol. 1996, 16 (8): 4504-4511.
Granger-Schnarr M, Benusiglio E, Schnarr M, Sassone-Corsi P: Transformation and transactivation suppressor activity of the c-Jun leucine zipper fused to a bacterial repressor. Proc Natl Acad Sci USA. 1992, 89 (10): 4236-4239. 10.1073/pnas.89.10.4236.
Bruder JT, Heidecker G, Rapp UR: Serum-, TPA-, and Ras-induced expression from Ap-1/Ets-driven promoters requires Raf-1 kinase. Genes Dev. 1992, 6 (4): 545-556. 10.1101/gad.6.4.545.
Bromberg JF, Horvath CM, Besser D, Lathem WW, Darnell JE: Stat3 activation is required for cellular transformation by v-src. Mol Cell Biol. 1998, 18 (5): 2553-2558.
Masker K, Golden A, Gaffney CJ, Mazack V, Schwindinger WF, Zhang W, Wang LH, Carey DJ, Sudol M: Transcriptional profile of Rous Sarcoma Virus transformed chicken embryo fibroblasts reveals new signaling targets of viral-src. Virology. 2007, 364 (1): 10-20. 10.1016/j.virol.2007.03.026.
Malek RL, Irby RB, Guo QM, Lee K, Wong S, He M, Tsai J, Frank B, Liu ET, Quackenbush J, et al: Identification of Src transformation fingerprint in human colon cancer. Oncogene. 2002, 21 (47): 7256-7265. 10.1038/sj.onc.1205900.
Paz K, Socci ND, van Nimwegen E, Viale A, Darnell JE: Transformation fingerprint: induced STAT3-C, v-Src and Ha-Ras cause small initial changes but similar established profiles in mRNA. Oncogene. 2004, 23 (52): 8455-8463. 10.1038/sj.onc.1207803.
Liu Y, Gao L, Gelman IH: SSeCKS/Gravin/AKAP12 attenuates expression of proliferative and angiogenic genes during suppression of v-Src-induced oncogenesis. BMC Cancer. 2006, 6: 105-10.1186/1471-2407-6-105.
Crisanti P, Pessac B, Calothy G: Avian neuroretina cells in oncogene studies. Methods Enzymol. 1995, 254: 77-97. full_text.
Garber EA, Mayer BJ, Jove R, Hanafusa H: Analysis of p60v-src mutants carrying lesions involved in temperature sensitivity. J Virol. 1987, 61 (2): 354-360.
Cross FR, Garber EA, Pellman D, Hanafusa H: A short sequence for the p60src n terminus is required for p60src myristylation and membrane association and for cell transformation. Mol Cell Biol. 1984, 4: 1834-1842.
Gagliardi M, Maynard S, Bojovic B, Bedard PA: The constitutive activation of the CEF-4/9E3 chemokine gene depends on C/EBPbeta in v-src transformed chicken embryo fibroblasts. Oncogene. 2001, 20 (18): 2301-2313. 10.1038/sj.onc.1204354.
Abramoff MD, Magelhaes PJ, Ram SJ: Image Processing with ImageJ. Biophotonics International. 2004, 11 (7): 36-42.
Li CW: DNA-Chip Analyzer (dChip). The analysis of gene expression data:methods and software. Edited by: Parmigiani EG, Irizarry R, Zeger SL. 2003, New York: Spinger
Li C, Hung Wong W: Model-based analysis of oligonucleotide arrays: model validation, design issues and standard error application. Genome Biol. 2001, 2 (8): RESEARCH0032-
Li C, Wong WH: Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci USA. 2001, 98 (1): 31-36. 10.1073/pnas.011404098.
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.
Kestler HA, Muller A, Gress TM, Buchholz M: Generalized Venn diagrams: a new method of visualizing complex genetic set relations. Bioinformatics. 2005, 21 (8): 1592-1595. 10.1093/bioinformatics/bti169.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.
Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4 (10): R70-10.1186/gb-2003-4-10-r70.
Draghici S, Khatri P, Tarca AL, Amin K, Done A, Voichita C, Georgescu C, Romero R: A systems biology approach for pathway level analysis. Genome Res. 2007, 17 (10): 1537-1545. 10.1101/gr.6202607.
Pawitan Y, Bjohle J, Amler L, Borg AL, Egyhazi S, Hall P, Han X, Holmberg L, Huang F, Klaar S, et al: Gene expression profiling spares early breast cancer patients from adjuvant therapy: derived and validated in two population-based cohorts. Breast Cancer Res. 2005, 7 (6): R953-964. 10.1186/bcr1325.
Ivshina AV, George J, Senko O, Mow B, Putti TC, Smeds J, Lindahl T, Pawitan Y, Hall P, Nordgren H, et al: Genetic reclassification of histologic grade delineates new clinical subtypes of breast cancer. Cancer Res. 2006, 66 (21): 10292-10301. 10.1158/0008-5472.CAN-05-4414.
Minn AJ, Gupta GP, Siegel PM, Bos PD, Shu W, Giri DD, Viale A, Olshen AB, Gerald WL, Massague J: Genes that mediate breast cancer metastasis to lung. Nature. 2005, 436 (7050): 518-524. 10.1038/nature03799.
Bild AH, Yao G, Chang JT, Wang Q, Potti A, Chasse D, Joshi MB, Harpole D, Lancaster JM, Berchuck A, et al: Oncogenic pathway signatures in human cancers as a guide to targeted therapies. Nature. 2006, 439 (7074): 353-357. 10.1038/nature04296.
Kaplan EL, Meier P: Nonparametric estimation from incomplete observations. Journal of the American Statistical Association. 1958, 53 (282): 457-481. 10.2307/2281868.
Papaconstantinou M, Wu Y, Pretorius HN, Singh N, Gianfelice G, Tanguay RM, Campos AR, Bedard PA: Menin is a regulator of the stress response in Drosophila melanogaster. Mol Cell Biol. 2005, 25 (22): 9960-9972. 10.1128/MCB.25.22.9960-9972.2005.
Petropoulos CJ, Hughes SH: Replication-competent retrovirus vectors for the transfer and expression of gene cassettes in avian cells. J Virol. 1991, 65 (7): 3728-3737.
Boettiger D, Roby K, Brumbaugh J, Biehl J, Holtzer H: Transformation of chicken embryo retinal melanoblasts by a temperature-sensitive mutant of Rous sarcoma virus. Cell. 1977, 11 (4): 881-890. 10.1016/0092-8674(77)90299-9.
Tato F, Alema S, Dlugosz A, Boettiger D, Holtzer H, Cossu G, Pacifici M: Development of 'revertant' myotubes in cultures of Rous sarcoma virus transformed avian myogenic cells. Differentiation. 1983, 24 (2): 131-139. 10.1111/j.1432-0436.1983.tb01312.x.
Pacifici M, Boettiger D, Roby K, Holtzer H: Transformation of chondroblasts by Rous sarcoma virus and synthesis of the sulfated proteoglycan matrix. Cell. 1977, 11 (4): 891-899. 10.1016/0092-8674(77)90300-2.
Menko AS, Boettiger D: Inhibition of chicken embryo lens differentiation and lens junction formation in culture by pp60v-src. Mol Cell Biol. 1988, 8 (4): 1414-1420.
Crisanti-Combes P, Lorinet AM, Girard A, Pessac B, Wasseff M, Calothy G: Effects of Rous sarcoma virus on the differentiation of chick and quail neuroretina cells in culture. Adv Exp Med Biol. 1982, 158: 115-122.
Yoshimura M, Iwasaki Y, Kaji A: In vitro differentiation of chicken embryo skin cells transformed by Rous sarcoma virus. J Cell Physiol. 1981, 109 (3): 373-385. 10.1002/jcp.1041090302.
Kruger A, Anderson SM: The v-src oncogene blocks the differentiation of a murine myeloid progenitor cell line and induces a tumorigenic phenotype. Oncogene. 1991, 6 (2): 245-256.
Gillet G, Michel D, Crisanti P, Guerin M, Herault Y, Pessac B, Calothy G, Brun G, Volovitch M: Serum factors and v-src control two complementary mitogenic pathways in quail neuroretinal cells in culture. Oncogene. 1993, 8 (3): 565-574.
Ladeda V, Frankel P, Feig LA, Foster DA, Bal de Kier Joffe E, Aguirre-Ghiso JA: RalA mediates v-Src, v-Ras, and v-Raf regulation of CD44 and fibronectin expression in NIH3T3 fibroblasts. Biochem Biophys Res Commun. 2001, 283 (4): 854-861. 10.1006/bbrc.2001.4845.
Jamal HH, Cano-Gauci DF, Buick RN, Filmus J: Activated ras and src induce CD44 overexpression in rat intestinal epithelial cells. Oncogene. 1994, 9 (2): 417-423.
Scholz G, Martinerie C, Perbal B, Hanafusa H: Transcriptional down regulation of the nov proto-oncogene in fibroblasts transformed by p60v-src. Mol Cell Biol. 1996, 16 (2): 481-486.
Gonneville L, Martins TJ, Bedard PA: Complex expression pattern of the CEF-4 cytokine in transformed and mitogenically stimulated cells. Oncogene. 1991, 6 (10): 1825-1833.
Calothy G, Poirier F, Dambrine G, Mignatti P, Combes P, Pessac B: Expression of viral oncogenes in differentiating chick embryo neuroretinal cells infected with avian tumor viruses. Cold Spring Harb Symp Quant Biol. 1980, 44 (Pt 2): 983-990.
Yang J, Mani SA, Donaher JL, Ramaswamy S, Itzykson RA, Come C, Savagner P, Gitelman I, Richardson A, Weinberg RA: Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis. Cell. 2004, 117 (7): 927-939. 10.1016/j.cell.2004.06.006.
Yu K, Lee CH, Tan PH, Hong GS, Wee SB, Wong CY, Tan P: A molecular signature of the Nottingham prognostic index in breast cancer. Cancer Res. 2004, 64 (9): 2962-2968. 10.1158/0008-5472.CAN-03-2430.
Sotiriou C, Wirapati P, Loi S, Harris A, Fox S, Smeds J, Nordgren H, Farmer P, Praz V, Haibe-Kains B, et al: Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis. J Natl Cancer Inst. 2006, 98 (4): 262-272.
Haibe-Kains B, Desmedt C, Piette F, Buyse M, Cardoso F, Van't Veer L, Piccart M, Bontempi G, Sotiriou C: Comparison of prognostic gene expression signatures for breast cancer. BMC Genomics. 2008, 9: 394-10.1186/1471-2164-9-394.
Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, Meijer-van Gelder ME, Yu J, et al: Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365 (9460): 671-679.
van't Veer LJ, Dai H, Vijver van de MJ, He YD, Hart AA, Mao M, Peterse HL, Kooy van der K, Marton MJ, Witteveen AT, et al: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415 (6871): 530-536. 10.1038/415530a.
Yu JX, Sieuwerts AM, Zhang Y, Martens JW, Smid M, Klijn JG, Wang Y, Foekens JA: Pathway analysis of gene signatures predicting metastasis of node-negative primary breast cancer. BMC Cancer. 2007, 7: 182-10.1186/1471-2407-7-182.
Vafa O, Wade M, Kern S, Beeche M, Pandita TK, Hampton GM, Wahl GM: c-Myc can induce DNA damage, increase reactive oxygen species, and mitigate p53 function: a mechanism for oncogene-induced genetic instability. Mol Cell. 2002, 9 (5): 1031-1044. 10.1016/S1097-2765(02)00520-8.
Lee AC, Fenster BE, Ito H, Takeda K, Bae NS, Hirai T, Yu ZX, Ferrans VJ, Howard BH, Finkel T: Ras proteins induce senescence by altering the intracellular levels of reactive oxygen species. J Biol Chem. 1999, 274 (12): 7936-7940. 10.1074/jbc.274.12.7936.
Bartkova J, Rezaei N, Liontos M, Karakaidos P, Kletsas D, Issaeva N, Vassiliou LV, Kolettas E, Niforou K, Zoumpourlis VC, et al: Oncogene-induced senescence is part of the tumorigenesis barrier imposed by DNA damage checkpoints. Nature. 2006, 444 (7119): 633-637. 10.1038/nature05268.
Di Micco R, Fumagalli M, Cicalese A, Piccinin S, Gasparini P, Luise C, Schurra C, Garre M, Nuciforo PG, Bensimon A, et al: Oncogene-induced senescence is a DNA damage response triggered by DNA hyper-replication. Nature. 2006, 444 (7119): 638-642. 10.1038/nature05327.
Mallette FA, Gaumont-Leclerc MF, Ferbeyre G: The DNA damage signaling pathway is a critical mediator of oncogene-induced senescence. Genes Dev. 2007, 21 (1): 43-48. 10.1101/gad.1487307.
Deichman GI, Kashkina LM, Mizenina OA, Gorojanskaya EG, Nikiforov MA, Gudkov AV, Dyakova NA, Komelkov AV, Prilutskaya MO, Kushlinsky NE, et al: Mechanisms of unusually high antioxidant activity of RSV-SR-transformed cells and of its suppression by activated p21ras. Int J Cancer. 1996, 66 (6): 747-752. 10.1002/(SICI)1097-0215(19960611)66:6<747::AID-IJC7>3.0.CO;2-#.
Kocanova S, Buytaert E, Matroule JY, Piette J, Golab J, de Witte P, Agostinis P: Induction of heme-oxygenase 1 requires the p38 MAPK and PI3K pathways and suppresses apoptotic cell death following hypericin-mediated photodynamic therapy. Apoptosis. 2007, 12 (4): 731-741. 10.1007/s10495-006-0016-x.
Zuckerbraun BS, Billiar TR: Heme oxygenase-1: a cellular Hercules. Hepatology. 2003, 37 (4): 742-744. 10.1053/jhep.2003.50139.
Choi JW, Yoo BK, Shin CY, Ryu MK, Ryu JH, el Kouni MH, Lee JC, Kim WK, Ko KH: Uridine prevents the glucose deprivation-induced death of immunostimulated astrocytes via the action of uridine phosphorylase. Neurosci Res. 2006, 56 (1): 111-118. 10.1016/j.neures.2006.06.004.
Choi JW, Shin CY, Choi MS, Yoon SY, Ryu JH, Lee JC, Kim WK, El Kouni MH, Ko KH: Uridine protects cortical neurons from glucose deprivation-induced death: possible role of uridine phosphorylase. J Neurotrauma. 2008, 25 (6): 695-707. 10.1089/neu.2007.0409.
Sohara Y, Ishiguro N, Machida K, Kurata H, Thant AA, Senga T, Matsuda S, Kimata K, Iwata H, Hamaguchi M: Hyaluronan activates cell motility of v-Src-transformed cells via Ras-mitogen-activated protein kinase and phosphoinositide 3-kinase-Akt in a tumor-specific manner. Mol Biol Cell. 2001, 12 (6): 1859-1868.
Platek A, Mettlen M, Camby I, Kiss R, Amyere M, Courtoy PJ: v-Src accelerates spontaneous motility via phosphoinositide 3-kinase, phospholipase C and phospholipase D, but abrogates chemotaxis in Rat-1 and MDCK cells. J Cell Sci. 2004, 117 (Pt 20): 4849-4861. 10.1242/jcs.01359.
Gimona M, Buccione R, Courtneidge SA, Linder S: Assembly and biological role of podosomes and invadopodia. Curr Opin Cell Biol. 2008, 20 (2): 235-241. 10.1016/j.ceb.2008.01.005.
Tehrani S, Faccio R, Chandrasekar I, Ross FP, Cooper JA: Cortactin has an essential and specific role in osteoclast actin assembly. Mol Biol Cell. 2006, 17 (7): 2882-2895. 10.1091/mbc.E06-03-0187.
Kempiak SJ, Yamaguchi H, Sarmiento C, Sidani M, Ghosh M, Eddy RJ, Desmarais V, Way M, Condeelis J, Segall JE: A neural Wiskott-Aldrich Syndrome protein-mediated pathway for localized activation of actin polymerization that is regulated by cortactin. J Biol Chem. 2005, 280 (7): 5836-5842. 10.1074/jbc.M410713200.
Yamaguchi H, Lorenz M, Kempiak S, Sarmiento C, Coniglio S, Symons M, Segall J, Eddy R, Miki H, Takenawa T, et al: Molecular mechanisms of invadopodium formation: the role of the N-WASP-Arp2/3 complex pathway and cofilin. J Cell Biol. 2005, 168 (3): 441-452. 10.1083/jcb.200407076.
Webb BA, Eves R, Mak AS: Cortactin regulates podosome formation: roles of the protein interaction domains. Exp Cell Res. 2006, 312 (6): 760-769. 10.1016/j.yexcr.2005.11.032.
Laurila E, Savinainen K, Kuuselo R, Karhu R, Kallioniemi A: Characterization of the 7q21-q22 amplicon identifies ARPC1A, a subunit of the Arp2/3 complex, as a regulator of cell migration and invasion in pancreatic cancer. Genes Chromosomes Cancer. 2009, 48 (4): 330-339. 10.1002/gcc.20643.
Assinder SJ, Stanton JA, Prasad PD: Transgelin: an actin-binding protein and tumour suppressor. Int J Biochem Cell Biol. 2009, 41 (3): 482-486. 10.1016/j.biocel.2008.02.011.
Nair RR, Solway J, Boyd DD: Expression cloning identifies transgelin (SM22) as a novel repressor of 92-kDa type IV collagenase (MMP-9) expression. J Biol Chem. 2006, 281 (36): 26424-26436. 10.1074/jbc.M602703200.
Morita T, Mayanagi T, Sobue K: Reorganization of the actin cytoskeleton via transcriptional regulation of cytoskeletal/focal adhesion genes by myocardin-related transcription factors (MRTFs/MAL/MKLs). Exp Cell Res. 2007, 313 (16): 3432-3445. 10.1016/j.yexcr.2007.07.008.
Yoshio T, Morita T, Kimura Y, Tsujii M, Hayashi N, Sobue K: Caldesmon suppresses cancer cell invasion by regulating podosome/invadopodium formation. FEBS Lett. 2007, 581 (20): 3777-3782. 10.1016/j.febslet.2007.06.073.
Eves R, Webb BA, Zhou S, Mak AS: Caldesmon is an integral component of podosomes in smooth muscle cells. J Cell Sci. 2006, 119 (Pt 9): 1691-1702. 10.1242/jcs.02881.
Vanhaesebroeck B, Jones GE, Allen WE, Zicha D, Hooshmand-Rad R, Sawyer C, Wells C, Waterfield MD, Ridley AJ: Distinct PI(3)Ks mediate mitogenic signalling and cell migration in macrophages. Nat Cell Biol. 1999, 1 (1): 69-71. 10.1038/9045.
Sawyer C, Sturge J, Bennett DC, O'Hare MJ, Allen WE, Bain J, Jones GE, Vanhaesebroeck B: Regulation of breast cancer cell chemotaxis by the phosphoinositide 3-kinase p110delta. Cancer Res. 2003, 63 (7): 1667-1675.
Duxbury MS, Ito H, Zinner MJ, Ashley SW, Whang EE: siRNA directed against c-Src enhances pancreatic adenocarcinoma cell gemcitabine chemosensitivity. J Am Coll Surg. 2004, 198 (6): 953-959. 10.1016/j.jamcollsurg.2004.01.037.
Griffiths GJ, Koh MY, Brunton VG, Cawthorne C, Reeves NA, Greaves M, Tilby MJ, Pearson DG, Ottley CJ, Workman P, et al: Expression of kinase-defective mutants of c-Src in human metastatic colon cancer cells decreases Bcl-xL and increases oxaliplatin- and Fas-induced apoptosis. J Biol Chem. 2004, 279 (44): 46113-46121. 10.1074/jbc.M408550200.
Chen T, Pengetnze Y, Taylor CC: Src inhibition enhances paclitaxel cytotoxicity in ovarian cancer cells by caspase-9-independent activation of caspase-3. Mol Cancer Ther. 2005, 4 (2): 217-224.
Hiscox S, Morgan L, Green TP, Barrow D, Gee J, Nicholson RI: Elevated Src activity promotes cellular invasion and motility in tamoxifen resistant breast cancer cells. Breast Cancer Res Treat. 2006, 97 (3): 263-274. 10.1007/s10549-005-9120-9.
Masumoto N, Nakano S, Fujishima H, Kohno K, Niho Y: v-src induces cisplatin resistance by increasing the repair of cisplatin-DNA interstrand cross-links in human gallbladder adenocarcinoma cells. Int J Cancer. 1999, 80 (5): 731-737. 10.1002/(SICI)1097-0215(19990301)80:5<731::AID-IJC17>3.0.CO;2-H.
Pham CG, Bubici C, Zazzeroni F, Knabb JR, Papa S, Kuntzen C, Franzoso G: Upregulation of Twist-1 by NF-kappaB blocks cytotoxicity induced by chemotherapeutic drugs. Mol Cell Biol. 2007, 27 (11): 3920-3935. 10.1128/MCB.01219-06.
Cheng GZ, Chan J, Wang Q, Zhang W, Sun CD, Wang LH: Twist transcriptionally up-regulates AKT2 in breast cancer cells leading to increased migration, invasion, and resistance to paclitaxel. Cancer Res. 2007, 67 (5): 1979-1987. 10.1158/0008-5472.CAN-06-1479.
Talamonti MS, Roh MS, Curley SA, Gallick GE: Increase in activity and level of pp60c-src in progressive stages of human colorectal cancer. J Clin Invest. 1993, 91 (1): 53-60. 10.1172/JCI116200.
Termuhlen PM, Curley SA, Talamonti MS, Saboorian MH, Gallick GE: Site-specific differences in pp60c-src activity in human colorectal metastases. J Surg Res. 1993, 54 (4): 293-298. 10.1006/jsre.1993.1046.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.
Morey JS, Ryan JC, Van Dolah FM: Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR. Biol Proced Online. 2006, 8: 175-193. 10.1251/bpo126.
Ding Y, Xu L, Jovanovic BD, Helenowski IB, Kelly DL, Catalona WJ, Yang XJ, Pins M, Bergan RC: The methodology used to measure differential gene expression affects the outcome. J Biomol Tech. 2007, 18 (5): 321-330.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/10/41/prepub
BDN is the recipient of scholarships from the Ministère de la Recherche, the Association pour la Recherche sur le Cancer (France) and the European Council of doctorate candidates and young researchers (Eurodoc). This work was made possible by grants from the Association pour la Recherche, the Ligue Nationale contre le Cancer and the Région Rhône Alpes (France) to GG, and the Canadian Institutes of Health Research (MOP#10272) to PAB.
The authors declare that they have no competing interests.
BMM performed microarray and data analyses, statistics, CEF cell culture, immunofluorescence, immuno-blotting, northern blot probe design and cloning, and contributed to experimental design. BDN derived and cultured the CNR cells and participated in experimental design. YW cultured CEF and performed northern blotting analyses. LW and NAR participated in CEF cell culture, treatment and western blotting analysis. GG participated in experimental design and coordination. PAB conceived and coordinated the study, derived the CEF, and contributed to experimental design. BMM and PAB wrote the manuscript. All authors read and approved the manuscript.
Electronic supplementary material
Additional file 9: Expression of CD44 in normal and transformed CEF. Surface expression of CD44 was examined by immunofluorescence in CEF infected with RCASBP(A), NY315 or SR-A RSV. NY72-4 infected CEF were either grown at the non-permissive (41.5°C) or the permissive temperature (37.5°C) for 24 hours prior to fixing and staining for CD44. (PDF 1 MB)
Additional file 10: Comparison of expression data of v-Src regulated genes as determined by gene profiling and northern blotting analyses. Gene expression by northern blotting analysis (Figure 3) was quantified and analyzed against microarray data to confirm a correlation in gene expression as measured by the two methods. Analysis of correlation of log2-transformed gene expression ratios (log2(experimental/baseline)) indicates a strong correlation (Spearman ρ of 0.83; p < 0.0001) between northern blot and microarray gene expression estimates. This ρ value is higher than typically observed in array validations , and consistent with northern blotting as a superior method of gene expression validation . A slope of 0.91 for the regression line indicates a nearly 1:1 ratio between log2 expression ratios of northern and microarray data. (PDF 66 KB)
Additional file 11: Comparison of gene expression data quantified from northern blots (Figure 3) compared to values obtained from microarray analysis.(DOC 54 KB)
Additional file 12: Full Pathway Express output summary. Common pathways found to be dysregulated in the Transformation-Regulated (TR), CEF NY72-4 and CNR NY72-4 gene sets are shown. The number of genes in the pathway refers to the number of genes in the associated KEGG pathway. Input genes refer to the number of differentially expressed genes that were found in that pathway. Corrected γ p-value is a measure of significance as calculated by Pathway Express. N.D. and N.S. indicate not determined and not significant (corrected γ p-value > 0.05) respectively. (DOC 93 KB)
Additional file 13: Activation of the PI3K/PKB-Akt pathway in v-Src-transformed CEF. The activation of PKB-Akt in RCASBP(A), NY315 and SR-A RSV infected CEF was investigated by western blotting analysis. Antibodies for PKB-Akt and the Ser-473 phosphorylated form of PKB-Akt were used to determine the expression and activation of PKB-Akt, respectively. The level of phospho-PKB-Akt was examined in cells treated with 1% DMSO (D; diluent) or the PI3K inhibitors LY290042 (LY) and wortmannin (W). PKB-Akt was hyper-phosphorylated in v-Src transformed CEF but phospho-PKB-Akt levels decreased upon treatment with the PI3K inhibitors. The activation of PKB-Akt coincided with increased Ser-9 phosphorylation of GSK3-β. (PDF 545 KB)
About this article
Cite this article
Maślikowski, B.M., Néel, B.D., Wu, Y. et al. Cellular processes of v-Src transformation revealed by gene profiling of primary cells - Implications for human cancer. BMC Cancer 10, 41 (2010). https://doi.org/10.1186/1471-2407-10-41