|Home | About | Journals | Submit | Contact Us | Français|
Microarray expression signature analyses have suggested that extracellular matrix (ECM) gene dysregulation is predictive of metastasis in both mouse mammary tumorigenesis and human breast cancer. We have previously demonstrated that such ECM dysregulation is influenced by hereditary germline-encoded variation. To identify novel metastasis efficiency modifiers, we performed expression QTL (eQTL) mapping in recombinant inbred mice by characterizing genetic loci modulating metastasis-predictive ECM gene expression. Three reproducible eQTLs were observed on chromosomes 7, 17 and 18. Candidate genes were identified by correlation analyses and known associations with metastasis. Seven candidates were identified (Ndn, Pi16, Luc7l, Rrp1b, Brd4, Centd3 and Csf1r). Stable transfection of the highly metastatic Mvt-1 mouse mammary tumor cell line with expression vectors encoding each candidate modulated metastasis-predictive ECM gene expression. Implantation of these cells into mice demonstrated that candidate gene ectopic expression impacts tumor progression. Gene expression analyses facilitated the construction of a transcriptional network that we have termed the ‘Diasporin Pathway’. This pathway contains the seven candidates, as well as metastasis-predictive ECM genes and metastasis suppressors. Brd4 and Rrp1b appear to form a central node within this network, which likely is a consequence of their physical interaction with the metastasis efficiency modifier Sipa1. Furthermore, we demonstrate that the microarray gene expression signatures induced by activation of ECM eQTL genes in the Mvt-1 cell line can be used to accurately predict survival in a human breast cancer cohort. These data imply that the Diasporin Pathway may be an important nexus in tumor progression in both mice and humans.
The clinical burden of solid tumor dissemination is immense, with most cancer-related deaths a consequence of metastasis. For example, those patients who present with non-disseminated breast cancer have a relatively good prognosis with the ten-year survival of these patients being approximately 80% . However, those with disseminated disease have a much poorer prognosis, with the median survival of patients with metastatic breast cancer at diagnosis estimated to be 2–4 years . The therapeutic advances of recent years have done little to improve this, which has led many to surmise that advanced disseminated breast cancer is an incurable condition .
These data misleadingly imply that disease presentation and course in solid cancers are uniform in their natures. This is not the reality in clinical practice however, with substantial variations in outcome frequently observed in patients with seemingly comparable disease patterns at presentation, risk factors and clinicopathological tumor features. There are several explanations as to why this is the case, with one plausible explanation being that current pathological classification systems for primary tumor characteristics do not capture sufficient information for highly accurate prediction of patient outcome. Additional variables may be required to achieve more accurate prognostic indexes. In recent years many investigators have been exploring molecular as well as histopathological markers as additional covariates to consider for prediction of disease outcome. Many reports have now demonstrated that primary tumors with higher propensities to metastasize exhibit characteristic gene expression patterns [3–5], supporting the use of molecular features in clinical prognostic evaluation. There are a number of possible reasons as to why this might be the case, with the most commonly cited being that metastasis-prone tumors develop somatic aberrations early in their evolution that drive these characteristic gene expression patterns .
Another plausible and somewhat complementary explanation for the disparate clinical behaviors of tumors with similar pathological features is the modifying effect of hereditary germline-encoded variation upon the disease process. Work in our laboratory using the highly metastatic polyoma middle-T (PyMT) transgenic mouse model has demonstrated the significant influence of germline variation on tumor progression [7,8]. Tumors induced by the same oncogenic event, activation of the PyMT transgene, on different mouse genetic backgrounds were shown to have significantly different predilections for pulmonary metastasis. These results suggested that metastasis susceptibility was likely to be due not only to somatic events occurring during the evolution of the tumor, but also due to combinations of inherited germline polymorphisms acting in concert to modulate the probability of whether a tumor will progress. Further investigation of these initial observations allowed us to identify two of these polymorphic metastasis susceptibility genes that modulate metastasis in both mouse and human: the Rap-GTPase activating protein (GAP) Sipa1 ; and Rrp1b, a gene that modulates the enzymatic activity of Sipa1 .
Additional studies have allowed us to examine whether there is likely an intimate relationship between hereditary polymorphism and tumor gene expression patterns, i.e. genetic polymorphism may not only modulate metastasis efficiency, but also may be an important factor in the induction of the prognostic signature profiles [11,12]. Examination of previously published metastasis-predictive signatures in both the mouse [11,12] and in humans [3–5] have revealed that dysregulation of extracellular matrix (ECM) gene expression is a common feature of tumors with a higher propensity to metastasize. Indeed, ECM gene dysregulation has been shown to be a very prominent feature of metastatic progression, and may well explain why highly metastatic mouse mammary tumor cell lines are typically more adhesive, invasive, and migratory than the less metastatic lines .
AKXD recombinant inbred (RI) mice  were used to determine whether this ECM dysregulation was at least partially under germline control. This methodology not only allowed us to identify Rrp1b as a dual ECM and metastasis efficiency modifier, but also to demonstrate that the cell growth regulator Bromodomain 4 (Brd4) [15–17], which also binds Sipa1, has similar properties with regards to metastasis and ECM gene expression (Bromodomain 4 activation predicts breast cancer survival; N.Crawford, J.Alsarraj, L.Lukes, R.Walker, H.Yang, M.Lee, K. Ozato, K.Hunter; PNAS, in press). Furthermore, we demonstrated that differential activity of BRD4 in primary breast cancers accurately predicts outcome in five different tumor gene expression datasets.
The current study represents a continuation of this earlier work to characterize novel metastasis efficiency modifier genes. By using techniques similar to those outlined above, we identify five novel candidate dual ECM and metastasis efficiency modifiers. Furthermore, our data demonstrate that these genes are elements of a potential novel transcriptional network, the activation of which is predictive of outcome in human breast cancer.
Microarray hybridization methodology and generation of the microarray expression data from AKXD × PyMT primary tumors has been described previously . Affymetrix CEL files were normalized using the RMA method, averaged for each AKXD RI strain, and loaded into the GeneNetwork web service (www.genenetwork.org) . GeneNetwork databases were then searched for probe sets within the eleven genes from our previously described metastasis signature profile  classified as being an “ECM component” (Col3a1, Fbln2, Fbn1, Mfap5, Mmp2, Mmp16, Nid1, Serpinf1, Serping1, Timp2, and Tnxb). Additionally, we also included probe sets from the two ECM genes (Col1a1 and Col1a2) represented the eleven gene human breast carcinoma metastasis gene signature profile described by Ramaswamy et al .
The Mvt-1 cell line, which is derived from an explant-cell culture of primary mammary tumors from MMTV-c-Myc/MMTV-Vegf bi-transgenic mice , was obtained as a gift from Lalage Wakefield (NCI, Bethesda). These cells were cultured as described . Mammalian expression vectors encoding the full-length cDNAs for all ECM eQTL candidate genes other than Brd4 were obtained from the Mammalian Gene Collection (Supplementary Table 1). A full-length cDNA for Brd4 in the expression vector pFLAG-CMV2 (Sigma-Aldrich, Saint Louis, MO) was obtained as a gift from Keiko Ozato (NICHD, Bethesda). The control cell line was generated using the vector pCMV-SPORT-β-Galactosidase (Invitrogen, Carlsbad, CA). Transfections and clonal isolations were performed as described previously . Colonies were screened by quantitative real-time PCR (qPCR) as described below to identify clones ectopically expressing candidate genes.
Total RNA was isolated from Mvt-1 clonal isolates as described . cDNA was synthesized from RNA isolated from either primary tumor tissues or transfected cell lines using the ThermoScript RT-PCR System (Invitrogen) by following the manufacturer’s protocol. Single RT-PCRs were performed for each Mvt-1 clonal isolate. SYBR Green or TaqMan qPCR was performed to detect the cDNA levels of ECM eQTL candidates and a variety of metastasis-predictive ECM genes (see above) using an ABI PRISM 7500 and/or 7900HT Sequence Detection Systems (Applied Biosystems, Foster City, CA). Primer and fluorogenic probe sequences for chromosome 7, 17 and 18 eQTL candidates and the housekeeping gene Peptidylprolyl Isomerase B (Ppib) are shown in Supplementary Table 2. Primers for metastasis-predictive ECM genes have been described elsewhere . Reactions were performed using QuantiTect SYBR Green Master Mix (Qiagen) or TaqMan Universal PCR Mastermix (Applied Biosystems) per the manufacturer’s protocol. The cDNA level of each eQTL candidate or metastasis-predictive ECM gene was normalized to Ppib cDNA levels using custom-designed primers for SYBR green-quantified target genes or custom-designed primers and fluorogenic probe for TaqMan-quantified target genes.
One hundred thousand cells stably expressing either as single eQTL candidate or pCMV-SPORT-β-Gal were subcutaneously implanted into six week old virgin FVB/NJ mice as described previously . These experiments were performed in compliance with the National Cancer Institute’s Animal Care and Use Committee guidelines. Animals were then aged for 28 days before being euthanized. Total tumor weight was determined by complete dissection of the primary tumor. Metastases counts were determined by surface pulmonary counts.
Transfected Mvt-1 cells were trypsinized and collected, and cell pellets lysed in cell lysis buffer [M-PER® Mammalian Protein Extraction Reagent (PIERCE, Rockford, IL) supplemented with Halt Protease Inhibitor cocktail kit (PIERCE, Rockford, IL)]. SDS-PAGE was performed for 60–90 min at 120 V using the XCell SureLock™ Mini-Cell (Invitrogen, CA) with NuPAGE Novex Bis-Tris Gels (Invitrogen, CA). Proteins were transferred to Immobilon™-P membranes (Millipore, Bedford, MA) and immunoblotted against Kai1 (SC-1087; Santa Cruz) and Nm23-H1 (SC-343; Santa Cruz) to detect the levels of these proteins in response to ectopic expression of ECM eQTL candidate genes. To ensure uniformity in sample loading onto the gel, membranes were re- immunoblotted against Gapdh (Santa Cruz) or β-Actin (Abcam).
RNA extraction and processing for Affymetrix GeneChip analysis have been described elsewhere. High confidence human transcriptional signatures of individual ECM eQTL candidate gene expression were generated in exactly the same manner as described for Rrp1b in our earlier work . Briefly, Affymetrix microarrays were used to compare gene expression in three to four Mvt-1 clonal isolates for ECM eQTLs and three Mvt-1/β-galactosidase clonal isolates. CEL files were analyzed using the Affymetrix GeneChip Probe Level Data RMA option of BRB ArrayTools 3.5.0. Genes with <1.5 fold-change from the gene's median value in 50% of samples, or a log-ratio variation P>0.01 were eliminated from analyses. To identify candidate gene expression signatures, the Class Comparison tool of BRB ArrayTools was performed, using a two-sample t-test with random variance univariate test. P-values for significance were computed based on 10,000 random permutations, at a nominal significance level of each univariate test of 0.0001. Tumor gene expression data from a well described breast cancer cohort  were downloaded from the Rosetta Company website (http://www.rii.com/publications/2002/vantveer.html).
Survival analysis was performed as described . BRB ArrayTools was used to perform unsupervised clustering based on individual ECM eQTL signature gene expression. Clustering was performed using average linkage, the centered correlation metric and center the genes analytical option. Samples were assigned into two groups based on the first bifurcation of the cluster dendogram, and Kaplan-Meier survival analysis was performed using the Survival module of the software package Statistica. Significance of survival analyses was performed using the Cox F-test.
The initial aim of this study was to define ECM expression quantitative trait loci (eQTLs). Traditionally, QTL analysis has involved defining specific genomic regions linked to a physiological trait (e.g. pulmonary metastasis burden). However, the advent of microarray technology has allowed global gene expression patterns to be treated as a quantifiable trait, which in turn has facilitated the characterization of eQTLs. An eQTL is therefore a genetically defined genomic locus associated with variation of gene expression (e.g. ECM gene expression). We utilized AKXD RI mice  for ECM eQTL mapping. RI panels are a specialized genetic mapping resource that has unique properties that make them particularly valuable for the study of complex traits [20,21]. The AKXD RI mice are a particularly useful tool for the study of germline-encoded metastatic propensity since they are derived from AKR/J, a highly metastatic strain, and DBA/2J, a weakly metastatic strain . Detailed methodology for the generation of ECM eQTLs has been described previously .
Briefly, extracellular matrix (ECM) genes are ubiquitous components of metastasis-predictive expression signatures in both human breast tumor tissue [3–5], and in PyMT-induced mouse mammary tumors [11,12]. Genetic mapping experiments were performed to identify genomic intervals associated with ECM component transcriptional control in our mouse metastasis model system using our previously derived microarray expression analyses of the PyMT-induced primary tumors in the AKXD RI panel . eQTLs were mapped using probe sets from thirteen previously identified metastasis-predictive ECM genes. A number of reproducible loci were evident  although the loci on chromosomes 7, 17 and 18 were of most interest due to their highly consistent correlation with expression of metastasis-predictive ECM probe sets. This is best illustrated when ECM eQTLs are mapped with those probe sets most strongly correlated with metastasis susceptibility (11 probe sets from a PyMT-induced metastasis-predictive gene signatures were identified , and five probe sets from the 17 gene human breast carcinoma metastasis gene signature profile described by Ramaswamy et al ). Such analyses demonstrate that expression of most of these sixteen metastasis-predictive ECM probe sets are least suggestively correlated with the loci on chromosomes 7, 17 and 18 (Supplementary Table 3).
Further evidence for the role of these chromosomal regions in modulation of metastatic capacity is derived from the observation that all three ECM eQTLs co-localized with suggestive metastasis QTLs observed after composite interval mapping of the AKXD data. We argue that this observation is consistent with the hypothesis that the ECM eQTLs and the metastasis efficiency QTLs might be causally linked. Chromosomal substitution strains have also been generated for chromosomes 7 and 17, both of which were shown to suppress metastatic capacity . Thus, evidence for a potential causal link between metastasis efficiency and ECM gene expression exists for all three chromosomes.
Correlation analysis was performed to identify potential candidates for the chromosomes 7, 17 and 18 ECM modifiers since the current prevailing hypothesis that most modifiers are likely to result from modest variations in gene expression levels or mRNA stability [23,24]. Whole genome correlation analysis of the microarray data was performed using the Trait Correlation function of the WebQTL database within the GeneNetwork web service  to identify probe sets with high correlation coefficients with each of the ECM genes of interest. Probes sets from the whole genome analysis that reside within the candidate eQTL intervals were then examined to identify genes that were reproducibly present in the ECM eQTL gene lists. Twenty-eight genes located within the region of chromosome 7 spanning the peak likelihood ratio statistic (a derivation of LOD score; LRS) score (physical locations ~ 56–121 Mb; fig. 1a) displayed both a high degree of correlation and a low P value with regard to expression of two or more of the 9 probe sets within metastasis-predictive ECM genes (Supplementary Table 4). Regarding the chromosome 18 ECM eQTL (physical location ~ 37–76 Mb; fig. 1b), expression of probe sets representing sixteen genes located within this locus were correlated with the expression of metastasis-predictive ECM probe sets (Supplementary Table 5). The chromosome 17 ECM eQTL has been described in detail elsewhere , but thirty candidate genes located within the peak LRS score region (physical locations ~ 18–40 Mb) were correlated with metastasis-predictive ECM genes.
Seven genes were observed in the three candidate intervals that demonstrated reproducible associations with ECM gene expression across the AKXD panel (Table 1). Although there are many plausible ECM eQTL candidates within the three eQTL loci (see Supplementary Table 4 & Supplementary Table 5 and Crawford et al ) we chose to further characterize the seven candidates for the following reasons. Three of these genes are of particular interest based on published literatures. Necdin (Ndn) is an imprinted transcription factor on chromosome 7 that has previously been implicated in cellular proliferation, collagen gene expression , and regulation of the metastasis-associated gene Hif1a [26,27]. Brd4, on chromosome 17, is a bromodomain-containing protein that associates with chromatin . Intriguingly, Brd4 has been demonstrated to be binding partner of our recently described metastasis efficiency modifier, Sipa1 . Further characterization of this gene in our laboratory demonstrated that differential expression of Brd4 drives a gene expression signature that predicts outcome in breast cancer (Crawford et al, in submission). Csf1r is a receptor for colony stimulating factor and is important in mammary development, oncogenesis and macrophage recruitment . Studies using knockout mice in the PyMT model have suggested that CSF signaling is an important component for metastatic progression . A fourth gene, Rrp1b, was also considered of great interest since functional genomics studies in our laboratory identified it as potentially binding to the polymorphic PDZ domain of Sipa1 . Furthermore, additional characterization of Rrp1b showed that this gene modulates metastasis efficiency in mice and that the human homolog of Rrp1b contains a polymorphism that is associated with various markers of outcome in breast cancer . The final three genes (Pi16, Luc7l and Centd3) were chosen for further characterization as they possess either a plausible functional role in metastasis (Table 1), a high degree of correlation with metastasis-predictive ECM probe sets (Supplementary Table 5 and Crawford et al ), or a combination of both.
The three reproducible ECM eQTL peaks were then examined to see whether they function in a common transcription regulatory pathway, rather than independently regulating ECM gene expression. Each of the ECM eQTL candidate genes was subjected to eQTL analysis to determine what regions of the genome were associated with steady state mRNA levels (data not shown). Probe sets from each gene in question were examined using the interval mapping function of GeneNetWork to identify genetically defined regions in the AKXD RI panel that were associated with probe intensity on the Affymetrix chips, as a surrogate for steady state gene expression levels. Interestingly, mRNA levels of Ndn, on chromosome 7 appeared to be influenced by the eQTLs on chromosomes 17 and 18, which co-localize with the ECM eQTLs. Analysis of the genes on chromosome 17 revealed the potential influence of the chromosome 18 locus on gene expression of all the candidates in the chromosome 17 locus. The data were therefore consistent with the possibility of a transcriptional cascade, with the chromosome 18 locus modulating expression of candidate genes within the chromosome 17 locus. This locus, in turn, appears to modulate the expression of Ndn on chromosome 7.
To evaluate the role that these genes play in the modulation of metastasis-predictive ECM genes and to gain support for the hypothesized transcriptional pathway, transfection experiments were performed. Mammalian expression vectors were obtained from the MGC collection (Supplementary Table 1) and transfected into the FVB/NJ mouse mammary tumor cell line Mvt-1 . Four independent clonal cell lines stably expressing individual ECM eQTL candidate genes as well as control cell lines stably expressing β-galactosidase were isolated by limiting dilution, and qPCR used to confirm candidate gene ectopic expression. Quantitative real-time PCR was also used to determine the expression of transcripts defining the three ECM eQTLs (Col1a1, Col1a2, Col3a1, Fbn1, Mmp2, Nid1, and Serping1; see Supplementary Table 3) as well as five related ECM genes (Col5a3, Col6a2, Fbln2, Mfap5, and Serpinf1). Ectopic expression of each of the candidates was shown to dysregulate the expression of different metastasis-predictive ECM genes (Table 2). The effects of ectopic expression of Rrp1b and Brd4 upon ECM gene transcription have been described elsewhere  (Crawford et al, in submission).
Mvt-1 cells ectopically expressing eQTL candidate genes were implanted in the vicinity of the mammary fat pad of FVB/NJ mice to assess the effect of those genes upon tumor growth and metastasis. Following a 28 day incubation period, mice were sacrificed and primary tumor weight was measured and pulmonary surface metastases counted (Table 3). The control cell line data represent the tumor growth and metastatic potential of one Mvt-1/β-galactosidase. One control cell line was used to minimize the number of animals in experiments since we have previously demonstrated all of the control cell lines to be indistinguishable for in vitro and in vivo growth and metastasis characteristics (unpublished data). The effects of ectopic expression of Rrp1b and Brd4 on tumor growth and metastasis in this spontaneous metastasis model have been described elsewhere  (Crawford et al, in submission).
Ectopic expression of each of the candidate genes profoundly impacted the tumor growth resulting from subcutaneous implantation of Mvt-1 clones (Table 3). This does not appear to result from a generalized reduction in cellular growth rate since previous growth curve analyses have demonstrated no differences in the in vitro growth characteristics of the ECM eQTL candidate and control cell lines (, Crawford et al in submission and unpublished observations). Furthermore, these data are concurrent with our earlier studies that have demonstrated that ectopic expression of both Rrp1b and Brd4 profoundly reduce tumor growth and metastasis  (Crawford et al, in submission).
A significant reduction in secondary lesion formation was also observed in all but the Mvt-1/Pi16 cell lines. However, given that the ECM eQTL candidate genes may be acting as growth suppressors in vivo, it is unclear at this stage whether this reduction in pulmonary metastatic burden is a primary event or secondary to the reduced tumor growth rate. In contrast to all of the other cell lines, the Mvt-1/Pi16 cell line appears to behave disparately: tumor growth is markedly reduced yet the metastatic potential of these cell lines has been maintained. Therefore, the cell lines ectopically expressing Pi16 display a net increase in the rate of secondary lesion implying that Pi16 might actually be a metastasis enhancer, at least in the highly metastatic Mvt-1 cell line. Further investigation will be required to see whether this is the case in other systems.
Affymetrix microarrays were used to compare gene expression in three to four Mvt-1 clonal isolates for ECM eQTLs and three Mvt-1/β-galactosidase clonal isolates. Probe sets significantly up- and down-regulated in response to ectopic expression of Ndn, Pi16, Luc7l, Centd3 and Csf1r according to these criteria are listed in Supplementary Table 6 to Supplementary Table 10, respectively. Microarray expression signatures indicative of activation of Rrp1b and Brd4 in the Mvt-1 cell line have been described elsewhere  (Crawford et al, in submission).
ECM eQTL mapping has raised the possibility that candidate genes are components of a transcriptional pathway (see above). The feasibility of such a functional relationship is strengthened when one considers that ectopic expression of each candidate in the Mvt-1 cell line has similar phenotypic effects (e.g. modulation of ECM gene expression, tumor growth and metastasis). If this is actually the case then dysregulation (in the form of ectopic expression) of individual candidate genes should impact the expression of other transcriptional pathway components. Therefore, in addition to the microarray analysis described above, the effect of ectopic expression of each candidate gene on the cellular transcriptional profile of Mvt-1 cells was also analyzed by qPCR (Supplementary Table 11). Additionally, we also used qPCR to analyze the expression of a number of other metastasis-related genes including Tgfb1, the metastasis suppressors Cd82/Kai1  and Tfpi  (Supplementary Table 11). A putative transcriptional network was constructed by applying the Occam’s Razor principle, i.e. the simplest explanation is likely correct, by minimizing the number of nodes necessary to explain the results. For example, ectopic expression of Pi16 was shown to up-regulate Brd4 and down regulate Tfpi, which was also down regulated by ectopic expression of Brd4. The down regulation of Tfpi by Pi16 could therefore be explained by action through the Brd4 pathway and the putative pathway drawn to represent this potential cascade. We have subsequently termed this putative network the ‘Diasporin Pathway’ based on the potential role in the diaspora, or dissemination, of tumor cells during the metastatic process. Simultaneously, yeast-two hybrid analysis revealed that our previously identified metastasis efficiency gene, Sipa1, putatively interacted with one of the novel chromosome 17 candidate genes (Rrp1b) . Combining these data suggested that Sipa1 and Rrp1b, plus Brd4 which had previously been identified as a Sipa1-interacting partner , play a key role in the transcriptional regulation of a number of metastasis-associated genes, including several metastasis suppressors (Fig. 2).
As a means of confirming the downstream effects of ECM eQTL candidate dysregulation as predicted by the Diasporin Pathway, we performed immunoblot analysis to quantify the protein levels of the metastasis-associated factors KAI1  and NM23-H1 . KAI1 is a well characterized factor that suppresses tumor metastasis by primarily inhibiting cancer cell motility and invasiveness . Both microarray and qPCR gene expression analyses suggested that ectopic expression of Diasporin Pathway components should enhance the expression of KAI1. Western blot analysis demonstrated, as expected, that Rrp1b expression significantly altered KAI1 levels (Fig. 3a). Ectopic expression of most of the other Diasporin pathway members also had similar effects which may reflect effects of changing the relative ratios of BRD4 and SIPA1 which previously has been shown to have major impact on cell cycle and cell physiology . Similarly, microarray analysis of Mvt-1 clonal isolates suggested that ECM eQTL candidate dysregulation should suppress expression of NM23-H1. Consistent with these observations, we did observe suppression of NM23-H1 in Rrp1b and Brd4 clonal isolates (Fig. 3b). Convincing evidence of protein alteration in the other cell lines was not observed, probably due the lack of sensitivity of this method to subtle changes in protein levels.
Our hypothesis that ECM eQTLs and metastasis susceptibility genes are one and the same leads to the further postulate that the transcriptional signature induced by ectopic expression of these putative dual function genes should induce gene expression signatures that are indicative of both ECM modulation and metastasis. This, in turn, leads to the possibility that such gene expression signatures could be used to predict survival in breast cancer. To evaluate this possibility, the Dutch breast cancer data of van de Vijver  were downloaded and the expression signatures induced by candidate gene ectopic expression mapped to the annotations of the Hu25K custom chip. The resulting gene list was then used to perform unsupervised clustering of human breast cancer microarray datasets. Kaplan-Meier survival analysis was subsequently performed, comparing the survival of the two major clusters formed by the first bifurcation of the dendogram representing high and low levels of candidate gene activation in primary tumor samples (Fig. 4).
Significant survival differences were observed for those gene expression signatures induced by activation of NDN (Fig. 4a), LUC7L (Fig. 4b), PI16 (Fig. 4c), CENTD3 (Fig. 4d), and CSF1R (Fig. 4e) in the Mvt-1 cell line. Additionally, we have previously demonstrated that the Mvt-1/Rrp1b and Mvt-1/Brd4 gene signatures have also been shown to predict outcome in the same breast cancer cohort  (Crawford et al, in submission). This implies that the level of Diasporin Pathway candidate gene activation or candidate gene-associated pathways within a tumor, presumably because of either somatic mutation or germline polymorphism, is an important determinant of the overall likelihood of relapse and/or survival. Further analysis indicated that the differential survival of breast cancer patients was associated with the effects of a distinct subset of genes induced by candidate activation (Supplementary Table 12). The degree of survival difference represented by each of these Diasporin Pathway candidate gene induced gene expression signatures was similar to the original 70-gene signature described by van't Veer and colleagues  (Fig. 4f).
We have used a multi-faceted experimental approach to identify a novel transcriptional nexus that appears to be involved in modulation of both ECM gene expression and tumor progression. We initially sought to map ECM eQTLs in AKXD recombinant inbred mice in order to identify those genomic regions where transcript expression displayed a significant degree of correlation with that of metastasis-predictive ECM genes. Our hypothesis here is that the aberrant ECM gene expression observed in tumors more prone to dissemination in both humans [3–5] and mice  is to some degree influenced by host germline polymorphism. We identified three loci on chromosomes 7, 17 and 18 that appear to be linked to metastasis-predictive ECM gene expression. Subsequently, correlation analysis was performed to facilitate candidate gene identification since the current prevailing hypothesis is that most modifiers are likely to result from modest variations in gene expression levels or mRNA stability [23,24]. This revealed a correlation between expression of a variety of genes within the peak region of linkage of each eQTL and the expression of various metastasis-predictive ECM genes (Supplementary Table 4 & Supplementary Table 5 and Crawford et al ). Seven genes were chosen for further analysis, for the reasons outlined above. We must state, however, that differential functionality of these seven candidate genes within AKXD RI mice is highly unlikely to account for all of the linkage observed with the three ECM eQTLs. Indeed, additional plausible candidate genes are evident within the peak regions of ECM eQTL linkage (Supplementary Table 4 & Supplementary Table 5 and Crawford et al ), and it may well be that these too have an effect upon expression of metastasis-predictive ECM genes as well as tumor progression. The aim of this study was not to perform an exhaustive analysis of all factors controlling metastasis-predictive ECM gene expression in AKXD mice, but to identify high-priority metastasis modifier candidate genes. These genes would then be used as a framework upon which we can build a more comprehensive picture of this ECM and metastasis-related transcriptional network.
This approach facilitated the identification of seven high priority ECM and metastasis modulator candidate genes. Although these genes have diverse cellular functions and localizations, our ECM eQTL data suggested that they may well share some form of functional relationship by virtue of their co-regulation with metastasis-predictive ECM genes. Indeed, ectopic expression of each putative ECM eQTL modifier in the Mvt-1 cell line proved that candidate activation not only modulates the expression of metastasis-predictive ECM genes, but also the expression of other members of the assumed transcriptional network. More significantly, we demonstrate that activation of each of these genes in the Mvt-1 cell line suppresses the inherent tumorigenicity of this highly aggressive cell line, with quite profound suppression of tumor growth in some instances (Table 3). Furthermore, we demonstrate that activation of all candidates other than Pi16 suppress metastasis, although we cannot say at present whether this is independent of the reduction of tumor growth rate. Interestingly, it appears that activation of Pi16 in the Mvt-1 cell line reduces tumor growth yet has no overall effect upon the ability of these cells to form secondary lesions. This implies that Pi16 activation in this cell line has facilitated a net increase in the metastatic capacity. Therefore, it would appear at this point that Pi16, unlike the other ECM eQTL candidates, is a tumor suppressor and a metastasis enhancer. Future experimentation will investigate the effects of Pi16 activation in a less aggressive cell line to see whether this observation can be replicated.
However, in spite of these encouraging initial data, the precise relationships between these genes, especially in the broader physiological sense, remain unclear. What is apparent is that the BRD4-SIPA1-RRP1B relationship seems to be a critical node in this ECM/tumor progression transcriptional network. SIPA1 is a potent modulator of metastatic efficiency in mice, which may relate to its RAPGAP activity . Furthermore, we have demonstrated that SIPA1 is a tumor progression susceptibility gene in humans . Short-hairpin RNA-mediated knockdown of Sipa1 in Mvt-1 cells dysregulates the expression of a variety of metastasis-predictive ECM genes on microarray gene expression analysis (unpublished observations). A plausible mechanism for this ECM dysregulation is that SIPA1 physically interacts with two ECM eQTL candidates identified in this study: Brd4  and Rrp1b . Furthermore, physical interactions between SIPA1 and each of these proteins have opposite effects upon the RapGAP enzymatic activity of SIPA1. Specifically, BRD4 increases the RAPGAP activity of SIPA1 , whereas RRP1B reduces its activity levels . The significance of these divergent effects upon SIPA1 is particularly intriguing, especially as both Brd4 and Rrp1b appear to have similar effects upon tumor progression. RAP1 activity has been shown to play an important role in metastasis [35,36]. Therefore, it appears crucial at this stage that we further investigate the BRD4-SIPA1-RRP1B relationship and the influence that it has upon RAP1 levels within the tumor cell to determine whether the effects of these factors are dependent or independent of this system.
The power and potential clinical relevance of this study becomes apparent when one considers that activation of each of Diasporin Pathway components induces a gene expression signature that can be used to predict outcome in human breast cancer. The approach that we have taken is somewhat similar to that of Bild et al , who showed that gene expression signatures can be identified that reflect the activation status of several oncogenic pathways. Specifically, this group demonstrated that in vitro adenoviral-mediated activation of well characterized oncogenic pathways can induce gene expression signatures that predict outcome in a number of different human cancers. We too demonstrate that observations taken from an in vitro setting can be used with effectiveness in a clinical setting. We do, however acknowledge that this approach does have limitations, the most notable of which is that we have only demonstrated the effectiveness of Diasporin Pathway gene expression signatures in one breast cancer cohort. Future efforts in this regards will concentrate upon determining whether the expression signatures identified here hold predictive value in other breast cancer populations as well as in different cancer sub-types. Nevertheless, if these expression signatures do hold substantial predictive value in multiple patient cohorts then it implies that they could enable clinicians to tailor the treatment of individual patients. Specifically, such tumor gene expression signatures could prove useful in prospectively identifying those patients with a seemingly low risk of relapse at presentation based on traditional clinical variables (e.g. node-negative/ER-positive patients) who are likely to eventually relapse. This in turn could facilitate initiation of adjuvant therapy that would not be administered under current treatment protocols. Indeed, the utility of tumor expression profiling is currently being investigated in clinical trials (e.g. TAILORx, MINDACT ). We believe that further studies like the current work that reveal novel expression signatures not only serve to enhance our understanding of human breast cancer progression, but more importantly hold the promise of improving upon assessment of prognosis at the time of presentation.
We thank Drs. Doug Lowy and Jude Alsarraj for critical comments on this manuscript. This research was supported (in part) by the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research.
Nigel P. S. Crawford, Laboratory of Cancer Biology and Genetics, National Cancer Institute, National Institutes of Health, Bethesda, Maryland 20892, USA.
Renard C. Walker, Laboratory of Cancer Biology and Genetics, National Cancer Institute, National Institutes of Health, Bethesda, Maryland 20892, USA.
Luanne Lukes, Laboratory of Cancer Biology and Genetics, National Cancer Institute, National Institutes of Health, Bethesda, Maryland 20892, USA.
Jennifer S. Officewala, Laboratory of Cancer Biology and Genetics, National Cancer Institute, National Institutes of Health, Bethesda, Maryland 20892, USA.
Robert W. Williams, Department of Anatomy and Neurobiology, University of Tennessee Health Science Center, Memphis, Tennessee, USA.
Kent W. Hunter, Laboratory of Cancer Biology and Genetics, National Cancer Institute, National Institutes of Health, Bethesda, Maryland 20892, USA.