To identify genes regulated synergistically by cooperating oncogenic mutations at genomic scale, we compared mRNA expression profiles of young adult murine colon (YAMC) cells with those of YAMC cells expressing mutant p53175H
(mp53), activated H-Ras12V
(Ras) or both mutant proteins together (mp53/Ras)17
using Affymetrix microarrays. Using a step-wise procedure, we first identified 538 genes differentially expressed between mp53/Ras and YAMC control cells with a statistical cut off at p < 0.01 (N-test, Westfall-Young adjusted). A further subset of 95 annotated genes that respond synergistically (28 up/67 down) to the combination of mutant p53 and Ras proteins, termed ‘cooperation response genes’ (CRG) was then determined using a synergy score, as described in methods (, Supplementary Table 1, Supplementary File 1
). Expression values and synergy scores for the CRGs derived from TaqMan low-density QPCR array (TLDA) data showed strong positive correlation with the values for the same genes obtained from microarray analysis (Supplementary Figures 1 and 2, Supplementary Table 2
and Supplementary File 2
). Thus CRG identification was confirmed by independent methods, with final CRG selection based on microarray data, due to higher sample replication in this data set.
Identification and characterization of cooperation response genes (CRGs)
CRGs encode proteins involved in the regulation of cell signaling, transcription, apoptosis, metabolism, transport or adhesion ( and Supplementary Table 1
), and in large proportion appear misexpressed in human cancer. For 47 of 75 CRGs tested co-regulation is found in primary human colon cancer and our murine colon cancer cell model (, Supplementary Figure 3
). Furthermore, altered expression of 29 CRGs has been reported in a variety of human cancer types, consistent with the direction of the change in gene expression observed in our experiments (, Supplementary Table 1 and references therein
). Thus, modulation of CRG expression has common features in malignant cell transformation of both murine and human cells.
Differential expression and synergy scores of CRGs in mp53/Ras cells and CRG co-regulation in human colon cancer
The relevance of differentially expressed genes for malignant cell transformation was assessed by genetic perturbation of a series of 24 CRGs and 14 genes responding to p53175H
and/or activated H-Ras12V
in a non-cooperative manner (non-CRGs). Perturbed genes were chosen across a broad range of biological functions, levels of differential expression and synergy scores (, Supplementary Figure 4
, Supplementary File 3
). Gene perturbations were carried out in mp53/Ras cells with the goal to re-establish mRNA expression of the manipulated genes to levels relatively close to those found in YAMC control cells, and to monitor subsequent tumor formation following sub-cutaneous injection of these cells into immuno-compromised mice. Of the perturbed genes, 18 were up- and 20 down-regulated in mp53/Ras cells, relative to YAMC.
Reversal of the changes in CRG expression significantly reduced tumor formation by mp53/Ras cells in 14 out of 24 cases (, left panel; ; Supplementary Figure 5A and Supplementary Table 3
), indicating a critical role in malignant transformation for a surprisingly large fraction of these genes. Perturbation of Plac8, Jag2 and HoxC13 gene expression had the strongest effects. We also combined perturbations of two CRGs, Fas and Rprm, that alone produced significant yet milder changes in tumor formation. This yielded significantly increased efficacy in tumor inhibition as compared with the respective single perturbations (, Supplementary Figure 5B and Supplementary Table 4
). Thus, even genetic perturbations of CRGs with relatively smaller effects when examined on their own show evidence of being essential when analyzed in combination.
Synergistic response of downstream genes to oncogenic mutations is a strong predictor for critical role in malignant transformation
CRG perturbations reduce tumor formation of both mp53/Ras and human cancer cells
In contrast to the multitude of CRG-related effects on tumor inhibition, out of the 14 non-CRG perturbations, only one showed a significant reduction in tumor formation of mp53/Ras cells (, right panel; Supplementary Figure 6, and Supplementary Table 5
). Taken together, our data suggest that among the genes differentially expressed in cancer cells, malignant transformation strongly relies on the class of genes synergistically regulated by cooperating oncogenic mutations (, Supplementary Figure 7
Genetic perturbation experiments were carried out utilizing retrovirus-mediated re-expression of corresponding cDNAs for down-regulated genes (Supplementary Table 6
) and shRNA-dependent stable knock-down using multiple independent targets for over-expressed genes (Supplementary Table 7
). In addition, Plac8 knock down was functionally rescued by expression of shRNA-resistant Plac8 (), confirming specificity of the Plac8 loss-of-function experiments. The extent of all gene perturbations was assessed by quantitative PCR (Supplementary Figure 8
). As expected, the genetic perturbations disrupt tumor formation downstream of the initiating oncogenic mutations. Expression of both mutant p53 and activated Ras proteins remains unaffected by all genetic manipulations that alter the formation of tumors (Supplementary Figure 9
). Moreover, gene perturbations distinguished tumor growth from in vitro cell proliferation, as they generally did not affect cell accumulation in tissue culture (Supplementary Fig 10
Perturbations of CRGs in human cancer cells (, Supplementary Figure 11 and Supplementary Tables 8 and 9
) had similarly strong tumor inhibitory effects to those in the genetically tractable murine mp53/Ras cells, as assessed by xenografts in nude mice. Perturbations of both up- and down-regulated CRGs, i.e. Dffb, Fas, HoxC13, Jag2, Perp, Plac8, Rprm, Zfp385 and Fas + Rprm were performed in human DLD-1 and/or HT-29 colon cancer cell lines using retroviruses (Supplementary Figure 12, Supplementary Tables 6 and 10
) as described above. Similar to mp53/Ras cells, both human cancer cell lines have p53 mutations, whereas with K-Ras (DLD-1) and B-Raf (HT-29) mutations they express activated members of the Ras/Raf signaling pathway distinct from activated H-Ras in mp53/Ras cells. In addition, DLD-1 and HT-29 cells carry further oncogenic lesions such as APC and PIK3CA mutations, with HT-29 cells also exhibiting a mutation in Smad4 (for references, see Supplemental Methods
). The genetic perturbations had no effect on mutant Ras/Raf or p53 protein expression levels in both DLD-1 and HT-29 cells (Supplementary Figure 13
), indicating disruption of the cancer phenotype downstream of oncogenic mutations. Taken together, these experiments indicate the relevance of CRGs to cancer in a variety of backgrounds and genetic contexts.
The data described here indicate that the cooperative nature of malignant cell transformation, to a considerable degree, depends on a class of downstream effector genes regulated synergistically by multiple oncogenic mutations. We show that these cooperation response genes (CRGs) identified here contain a strikingly large fraction of genes (14 out of 24 tested) that are critical to the malignant phenotype, and that their perturbation, singly or in combination, can inhibit formation of tumors containing multiple oncogenic lesions, including p53 deficiency. In contrast, few of the genes differentially expressed in a non-synergistic manner (1 out of 14) significantly reduced tumor growth upon perturbation. Synergistic behavior found in gene expression data thus appears highly informative for identification of genes critically involved in malignant cell transformation (), and provides a rational path to discovery of both cancer cell-specific vulnerabilities and targets for intervention in cancer cells harboring multiple mutations, including p53 loss-of-function.
CRGs represent a set of 95 annotated cellular genes, many of which have been associated with human cancer by virtue of altered gene expression (, Supplementary Table 1
). They are involved in the regulation of cell signaling, transcription, apoptosis and metabolism, and based on our data represent key control points in many facets of cancer cell behavior. We thus consider CRGs as critical nodes in gene networks underlying the malignant phenotype, providing an attractive rationale to explain why several features of cancer cells emerge simultaneously out of the interaction of a few genetic lesions17
Among CRGs and other differentially expressed effector genes we also have identified examples that when perturbed produce significantly larger tumors ( and Supplementary Tables 3 and 5
). This is consistent with the notion that oncogenic mutations can induce strongly anti-proliferative cellular stress responses18-21
. The existence of genes that while responding to oncogenic mutations restrict tumor formation provides direct evidence to support the idea that the state of malignant transformation arises as the result of a finely tuned balance between opposing signals generated by oncogenic mutations15-17, 20, 22, 23
. It is thus reasonable to speculate that tumor suppression via perturbation of CRGs, as shown here, may involve the disruption of this delicate balance. In fact, such targeted disruption downstream of oncogenic mutations may allow selective cancer cell deconstruction yielding intervention strategies with high specificity for cancer cells.
For the 14 CRGs with tumor-inhibitory perturbations, a clear causal role in tumor formation downstream of oncogenic mutations has been shown here for the first time. Moreover, our data indicate that both gene extinctions (eight genes) and gene inductions (six genes) play important roles in this process. For example, we show that re-expression of the down-regulated CRGs Jag2, a Notch ligand, or of HoxC13, a homeobox transcription factor, as well as shRNA-dependent knock down of Plac8 gene expression are each strongly tumor inhibitory in p53 defective murine and human cancer cells. Both Notch signaling24
can play oncogenic roles in haematopoietic malignancies, but are involved in promoting differentiation of epithelial cells26, 27
consistent with the tumor-inhibitory function of Jag2 and HoxC13 in the context of the solid tumor models investigated here. Plac8 is a little investigated gene encoding a cysteine-rich highly conserved peptide expressed in placenta, haematopoietic and epithelial cells that is non-essential for mouse development28
. When over-expressed, Plac8 can suppress p5329
. Its essential role for tumor formation of p53-deficient cancer cells, however, is novel and unexpected. Among the eight down-regulated CRGs is Zfp385, another gene of unknown function. Moreover, there is a considerable number of pro-apoptotic/anti-proliferative genes such as Perp, Rprm, Fas, Dffb and Wnt9a, indicating that Ras activation and p53 deficiency cooperate to extinguish the expression of multiple growth inhibitory genes, each of which contributes significantly to restricting tumor growth in the YAMC model when re-expressed. Out of these genes, Perp, Rprm, and Fas previously have been identified as direct p53 targets, suggesting that their regulation by p53 is highly conditional on Ras activity (Supplementary Table 1 and references therein
). Most of the up-regulated CRGs contributing to tumor growth affect signal transduction. This includes Fgf7, Rgs2, Gpr149, an uncharacterized orphan seven-trans-membrane receptor, and Sod3, which acts on signaling via modulation of metabolites30
. For all of these genes, including Pla2g7, a role in promoting tumor growth is reported here for the first time.
Notably, the efficacy of CRG perturbations performed in human colon cancer cells was comparable to that in the murine colon cell transformation model, suggesting dependence of the malignant state on a similar set of genes in both backgrounds. This is remarkable in light of the fact that these human cancer cells carry oncogenic mutations in genes in addition to Ras or Raf and p53, and suggests that CRGs may play key roles in the generation and maintenance of the cancer cell phenotype in a variety of contexts. CRGs thus may provide a valuable source for identification of much sought ‘Achilles heels’ in human cancer by rational means.