Skip to main content

Context-specific gene regulatory networks subdivide intrinsic subtypes of breast cancer

Abstract

Background

Breast cancer is a highly heterogeneous disease with respect to molecular alterations and cellular composition making therapeutic and clinical outcome unpredictable. This diversity creates a significant challenge in developing tumor classifications that are clinically reliable with respect to prognosis prediction.

Results

This paper describes an unsupervised context analysis to infer context-specific gene regulatory networks from 1,614 samples obtained from publicly available gene expression data, an extension of a previously published methodology. We use the context-specific gene regulatory networks to classify the tumors into clinically relevant subgroups, and provide candidates for a finer sub-grouping of the previously known intrinsic tumors with a focus on Basal-like tumors. Our analysis of pathway enrichment in the key contexts provides an insight into the biological mechanism underlying the identified subtypes of breast cancer.

Conclusions

The use of context-specific gene regulatory networks to identify biological contexts from heterogenous breast cancer data set was able to identify genomic drivers for subgroups within the previously reported intrinsic subtypes. These subgroups (contexts) uphold the clinical relevant features for the intrinsic subtypes and were associated with increased survival differences compared to the intrinsic subtypes. We believe our computational approach led to the generation of novel rationalized hypotheses to explain mechanisms of disease progression within sub-contexts of breast cancer that could be therapeutically exploited once validated.

Background

Complex diseases such as breast tumors frequently have genomic mutations, translocations, and increased or decreased dosage of genes. The complex regulatory arrangements are further permuted, producing extreme heterogeneity in regulation and severe analytic complications. Such heterogeneity prevents existing methods, which often assume a certain level of homogeneity in samples, from learning underlying regulatory mechanisms from molecular measurements of tumor tissues. This inherent heterogeneity also generates a need for specialized therapeutic response, necessitating the development of models of breast cancer that can incorporate such heterogeneity.

Several landmark studies have shown that array-based expression profiling can provide insight into the complexity of breast tumors and can be used to 1) derive a molecular taxonomy for breast cancer, and 2) provide prognostic information better than standard assessment of clinical variables [1]. For example, genomic grade, or proliferation index is a strong predictor of outcome in estrogen receptor alpha (ER) positive disease. Another example is the 21-gene OncotypeDx assay (Genomic Health, Redwood City, CA) used to stratify ER positive patients into risk of recurrence groups following endocrine therapy. From seminal work published by Dr. Charles Perou [2] and others, classification methods have been, and continue to be, used to define “intrinsic” subtypes of breast cancer. These subtypes include Luminal A, Luminal B, Basal-like, HER2-enriched and normal breast-like, and are believed to represent distinct biological entities. Moreover, multiple studies have now confirmed that patient survival significantly differs with respect to intrinsic subtype.

A pathway-based classification of breast cancer shows that intrinsic gene expression signatures can be built using knowledge from pathway activity on previously known subtypes [3]. The aim of the study was to provide a functional interpretation of the gene expression data that can be linked to therapeutic options. The paper by Gatza et al. [3] indicates that the intrinsic subtypes can have further subgroups which may lead to much better understanding of each subtype. Recently, a subgroup of Basal-like tumors associated with poor prognosis has also been reported [4, 5].

Aim of this work

To improve the modeling and inference of regulatory mechanisms from such heterogeneous samples, a biologically based approach to sample and process stratification that models and learns context-specific regulations was proposed and developed [6, 7]. The model hypothesizes that genomic (expression) regulation is comprised of two distinct types: convergent regulation and divergent regulation, the former representing a particular set of genes being modulated by different sets of regulators, and the latter indicating a given set of regulators modulating entirely different sets of genes in different cellular contexts. The model also assumes that when a cell maintains a specific cellular context, (i.e a phenotype) it tightly regulates a battery of genes. It is hypothesized that the set of genes under such tight regulation would show rather deterministic transcriptional activities. When the cell moves away from this cellular context or changes to a different cellular state, the behavior of the same set of genes will not appear as deterministic since their behavior is now under the control of various external agents. In this paper, we will illustrate, using the concepts of conditioning and crosstalk, that systematic inquiry of candidate genes can identify a set of cellular contexts where a set of genes is tightly regulated, and corresponding context-specific gene regulatory networks.

Genomic regulation of breast cancer subtypes may show several common traits, although they have several unique features that make them distinct. The contexts obtained from this approach can be further used to study the underlying biology of the individual subtypes, which can lead to a better understanding of the differences and similarities between the tumors.

In contrast to previous methods, we used an unsupervised method to identify biologically meaningful cellular contexts within breast cancer. Our motivation lies in modeling the heterogeneity of breast cancer with a context-specific approach.

Results and discussion

The results section describes the data collection process, followed by the context analysis, phenotype and functional enrichment analysis and survival analysis.

Breast cancer data collection and processing

Ten breast cancer Affymetrix HG-U133A microarray data sets were downloaded from the NCBI GEO data repository (http://www.ncbi.nlm.nih.gov/geo/). These cohorts contain distinct clinical and molecular features such as ER+/ ER-, PgR+/ PgR-, Grade and LN+ and LN- types. Table 1 lists the data sets along with the number of samples within each cohort. The data from all cohorts were combined and normalized together by RMA normalization. A 2-fold change was used to categorize genes as under-expressed, no change or over-expressed; thus generating a data with ternary values {-1, 0, 1}. The cohorts contain a total of 1,887 samples with some samples repeated in more than one cohort. After removing the duplicates, a total of 1,636 samples were obtained. Additionally, GSE 2603 contains some cell line data that was removed reducing the number of samples to 1,614.

Table 1 Breast cancer cohorts

Many variables in the data sets have low variance and may not contribute to network learning. These variables with low variance across all samples were removed from the data sets. This also reduced the dimensionality of the data and made the network learning process computationally more tractable. Affymetrix probe sets were matched to HUGO gene symbols, probes matching to the same genes were combined by taking the median of the probes with Spearman’s correlation of 0.8. Probe sets with lower correlation values were discarded. After filtering at a variance of 0.14 and combining probes, we reduced the variable size to 5,023 highly variant genes.

Context analysis

A context-specific gene regulatory network was generated for the data using a parallel implementation of the algorithm called ExPattern (available at http://sysbio.fulton.asu.edu/expattern). The steps involved in finding contexts from the breast cancer expression data is illustrated in Figure 1. A graph with context-motifs filtered at a statistical significance of < 0.05 after FDR correction was generated. A total of 1,466 context-motifs generated at this step were clustered using Markov clustering (MCL) [8] to obtain 189 clusters, which are referred to as ”contexts” henceforth in the paper. MCL was performed on the graph with an inflation of 3.0 to keep the granularity high, and connectivity was imposed within clusters, such that each context contained connected context-motifs only. Contexts with less than 80 samples (< 5% of total samples) may not convey meaningful results and thus were discarded, resulting in 41 contexts. Specificity of the contexts was measured by computing pairwise Jaccard distance between the contexts for both samples and genes [9]. The contexts had an average Jaccard distance of 0.96 for genes and 0.85 for samples, indicating that most of the contexts were well separated with little overlap. A summary of context analysis with respect to the number of associated samples and genes is given in Table 2.

Figure 1
figure 1

Context-Mining process flow. The process to analyze heterogeneous biological data to learn context-specific gene regulations is illustrated in this figure. We first identify context-motifs using crosstalk, conditioning and statistical p-value computations. Since some genes can be a driver in a context motif, but a passenger in other context motifs, these context motifs can be chained together to build a interaction graph. In this graph, each edge represents an interaction specific to certain subset of samples (context motif). We now use this property along with graph clustering to identify potential cellular contexts where we should see a set of interactions sharing significant numbers of samples in common. Once cellular contexts are identified, we annotate each context (which includes a subset of samples and a subset of genes) using gene enrichment, subtype enrichment, or survival analysis methods as described in the paper.

Table 2 Contexts summary

Clinical characterization and subtype enrichment

Following clustering, the contexts were analyzed for clinical and molecular marker enrichments. Additionally, intrinsic subtypes were also associated with contexts with statistically significant enriched subtypes. Clinical and molecular markers and intrinsic subtypes associated with each context are listed in Table 2. A reasonably large number of contexts showed enrichment for at least one subtype. The grouping of ER+ intrinsic subtypes (LumA, LumB and Normal) and ER- tumors (Her2 and Basal-like) was clearly evident with the context enrichment. Basal-like tumors associated with low survival, showed high grade consistent with previous studies of Basal-like breast cancer. Additionally, LumA and LumB types were enriched with more than one context and Basal-like tumors were enriched in several contexts. Average Jaccard distance of samples for LumA contexts is 0.75 and LumB context is 0.85. There were no overlapping genes between the LumA and LumB contexts. The average Jaccard distance of samples for Basal-enriched contexts was 0.84, indicating that these groups are highly distinct and may indicate subgroups of Basal-like tumors. Table 2 shows some contexts enriched with multiple intrinsic subtypes, and we studied this further by grouping contexts and intrinsic subtypes based on their co-enrichments, via hierarchical clustering. Enrichments were annotated with ternary values 1, 0, -1, indicating presence, absence and, in the case of some clinical features, presence of negative types. Clinical enrichments ER, PgR, LN status and Grade were encoded as “-1” for ER-, PgR-, LN- and Low grade tumors, respectively, and positive “1” for ER+, PgR+, LN+ and high grade tumors, respectively. Hierarchical clustering was performed using Hamming distance and clusters were chained with complete linkage. The result is shown in Figure 2, which indicates biologically relevant groups for subtypes and clinical features. For example, Basal-like tumors known to be associated with high grade are clustered with grade. Luminal A tumors group with Normal-like tumors and Luminal B group with Her2-like tumors. Additionally, correspondence between ER and PgR states is also observed in the clustering result.

Figure 2
figure 2

Context clusters. Hierarchical clustering of contexts with clinical and subtype enrichments

Functional annotation

Functional annotation on the contexts with gene sets from MSigDB revealed interesting results. The results validate the enrichment of the contexts with ER+ and ER- tumors, and gene sets pertaining to these characteristics were found. Context 16 an ER+ and Luminal-like enriched context showed significant enrichment with Luminal-like breast cancer gene sets (p-values: 6.00E – 12, 1.38E – 10, 1.07E – 08). Context 48, ER-, high grade, Basal-like context was enriched with ER- gene sets and with invasive breast cancer gene sets (p-values: 0.00E + 00). Context 168 (ER-, Basal-like context) showed enrichment with ER- breast cancer gene sets and with Basal-like breast cancer gene sets (p-values: 1.55E – 04, 3.32E – 06). Additional pathways for some selected contexts are included in the Supplement tables 1 - 7 (see Additional file 1 Supplement tables 1-7).

Survival analysis

Survival analysis was performed on the 436 samples out of 1,614 with survival data (see Table 3)The Kaplan-Meier plot in Figure 3 with survival of Basal-like tumors, demonstrates the difference with rest of the tumors (non-Basal) with disease free survival (DFS) as the endpoint. The Kaplan-Meier plot of Basal-like enriched context 130 (genes: GATA3, INPP4B) in Figure 4 not only indicates shorter survival as expected for higher grade, ER- tumors but also a larger separation from the rest of the samples including other Basal-like tumors. Comparison of Figures 3 and 4 clearly indicates a potential sub-grouping within Basal-like tumors. Kaplan-Meier plot of Context 51 (genes: BUB1, DLG7, CENPA, MAD2L1, TTK, MCM10) ER+ tumors also indicates a better survival of ER+ tumors compared to rest of the samples (Figure 5).

Table 3 Context sample survival
Figure 3
figure 3

Survival for all Basal-like samples. Survival plot (in years) for all Basal-like tumors compared to rest of the tumors (all non Basal-like).

Figure 4
figure 4

Survival for Basal-like Context 130. Survival plot (in years) for Context 130 enriched within a subgroup of Basal-like tumors shows poor survival compared to all Basal-like tumors. GATA3 which was under-expressed in this context was correlated with increased tumor size and estrogen and progesterone receptor negativity [20], confirming the poor survival indicated in this context.

Figure 5
figure 5

Survival for ER+ context samples. Context 51 enriched with low grade, ER+ tumors indicates expected increase in survival compared to rest of the tumors.

Discussion

Several contexts of biologic interest and potential translational potential were highlighted by this analysis that appear both expected, and novel. Context 51, indicative of ER-positive and PgR-positive, low grade, Luminal A and normal-like tumors, was significantly enriched for genes associated with cell cycle checkpoint regulation, specifically, the M phase of mitotic cell cycle (BUB1 MAD2L1 TTK). As would be expected for ER+ low grade tumors, which tend to exhibit lower levels of proliferation, this context correlated with an increase in median survival (Figure 5 p = 7.8997e 10–8). Context 89 shared the same enriched subtypes as context 51, and contained just 2 genes from the same family, MAGEA3 and MAGEA6. The potential utility of MAGEA (Melanoma Antigen family A) proteins as a biomarker of the presence of micrometastases and circulating tumor cells has been previously reported [10]. We noted that in this instance, the MAGEA genes were associated with tumors that typically have better outcome. It is interesting to speculate whether analysis of MAGEA proteins in circulating breast tumor cells or micrometastases may enhance prognostication in stage III or IV breast cancer. This has not yet been studied. Contexts 57, 48 and 145 were three of several contexts associated with the Basal-like intrinsic subtype and high grade tumors, each with strikingly different apparent molecular underpinnings. Context 57 contained genes (e.g., TEK) suggestive of highly angiogenic Basal-like breast tumors [11]. This tumor context includes positive lymph node status and a decrease in median survival (5.9 vs 7.3 months). In contrast to context 57, context 48 which contained 27 genes, was significantly associated with cell cycle, with no significant difference in prognosis, perhaps due to low numbers of tumors with survival data within this context. Context 130, a Basal-like context has under-expression of GATA3 which is in concordance with previous studies of Basal-like subgroup, ’claudin-low’ with poor prognosis and more refractory to chemotherapy [5] . Lastly, context 145, again a Basal-like context of high tumor grade and ER negative status contained genes associated with deregulated secretory pathways and mechanisms of docking and fusion of vesicles to target membranes. The gene PSENEN in this context codes for a gamma secretase and is known to play a role in intramembranous processing of proteins such as Notch, a key mediator of cell-fate, tissue patterning and morphogenesis. PSENEN protein is required for Notch pathway signaling [12] and Notch signaling is deregulated in breast cancer [13]. Interestingly, Prat et al have also identified a subtype of Basal-like breast cancer with Notch-associated signaling deregulation [4]. Additional genes in context 145 (such as, MAP3K2) point to deregulated MAPK, NFkB and PKC signaling, all of which are oncogenic in breast cancer and have been reported to be linked to Notch deregulation. As Notch signaling is emerging as an attractive therapeutic target in breast and other cancers [13], this context was of particular interest. There was only one sample with survival data in context 145 for prognostic evaluation, however the trend was an association with poor survival. Context 124 is consistent with the low survival of patients with LumB tumors (p < 1.1897e 10–7). The above summarizes a sampling of contexts which highlight important unanswered questions in translational breast cancer research. Validation of these hypotheses to explain mechanisms of disease progression within sub-contexts of breast cancer have a potential to be therapeutically exploited.

There are a number of well characterized commercially available breast cancer cell lines that mimic various stages of breast cancer progression and biologic characteristics (including luminal A, HER2 enriched, Basal-like, invasive, non-invasive, metastatically competent, etc). Genes of interest identified as part of a specific context can be experimentally manipulated in vitro using breast cancer cell lines that match the phenotypic and/or molecular context of interest. Techniques commonly used to manipulate an individual gene within a viable cell line include RNA interference technology, which specifically eliminates expression of any specified target gene, use of target-selective drugs, or use of exogenous DNA gene expression constructs, which are engineered to introduce and express a specific gene of interest in a cell. The biological and molecular consequences of manipulating expression of a specific gene can then be measured using cell-based and/or molecular techniques to validate a computationally predicted hypothesis. Once verified, this information can be leveraged to develop more accurate prognostic or predictive biomarkers for clinical application.

Conclusions

This papers demonstrates the application of context-specific gene regulatory networks to identify biological contexts within heterogeneous breast cancer data over many samples. This large sample set identifies a finite number of contexts linked with intrinsic subtypes and clinical parameters. Diagnosis of intrinsic subtype is an important step that aids the prognostics for breast cancer. Our analysis of intrinsic subtype gene expression signatures is consistent with previous findings of individual cohort molecular profiling studies. Previously established intrinsic subtypes show different mechanisms indicating a possibility of further grouping of the intrinsic subtypes. Distinct contexts of Basal-like tumors confirm the existence of subgroups within Basal-like tumors as reported in previous studies. The contextual drivers identified for these subgroups can help explain the molecular aspects for the groups. Several new genes were found driving some contexts that have not been previously reported to be associated with known subgroups within these subtypes. Functional annotation of the genes associated with contexts also revealed different characteristics associated with each subgroup that can be biologically validated to define signatures for the groups.

Future work

The results of the experiments in the paper provide a promising approach to finding gene and clinical signatures associated with intrinsic subtypes within breast cancer. Nevertheless, biological validation of the genes involved is necessary and can strengthen the signatures for each context. Future directions include testing the results on a independent data set to group subtypes.

Methods

In this section, we first describe succinctly an approach to infer context-specific gene regulatory networks [7], [14], [15], a metric to associate samples with appropriate context, and then describe statistical tests to identify pathways and clinical phenotypes that are enriched in context.

Inferring context-specific gene regulatory networks

Previously, we developed a method to infer context-specific gene regulatory network from gene expression data [7], [14], [15]. In this section, we describe the method that we have further refined since then, by introducing context-motif mining, followed by graph-clustering of context-motifs to infer contexts and corresponding context-specific gene regulatory networks.

Mining context-motifs

Given a gene g k as a driver gene and a condition defined by a subset of samples M j , the algorithm uses probabilistic measures to identify a set of genes, i.e. passenger genes, that show a coherent molecular pattern within the condition. We define this set of genes, one or more of which function as drivers and the others as passenger genes, context-motif. Formally, a context-motif is represented as C i = (G i , Y i , S i , M i ) where G i represents a set of driver genes, Y i the possible states of the genes (an example would be -1, 0, +1 for a ternary quantized data set), S i a set of passenger genes, and M i the set of samples under which coherent expression is observed.

Coherence of expression pattern and its specificity are measured by two statistics, conditioning (δ k ) and crosstalk (η k ), as given in Eqs. 1 and 2, which determine if a gene k displays a cohesive expression pattern specific to a cellular context regulated by Y=1, where X k is state of driven genes.

δ k = 1 – P(X k = 1 | Y = 1),   (1)

η k = P(X k = 1 | Y ≠ 1)     (2)

Conceptually, conditioning measures the lack of transcriptional coherence in the condition of interest and crosstalk measures the specificity of coherence. This is based on the property that, cell deviates from its regulatory behavior under environmental changes or, in this study, more specifically, the presence of tumor. A change in the cellular context can be used to condition a subset of samples.

Since both crosstalk and conditioning parameters are estimated from observations, the statistical significance (p-value) of these parameters is computed by hypergeometric probability, to determine whether the patterns found in this case are not by chance.

The algorithm to identify all potential context-motifs interrogates every gene in the data set as a potential driver gene (G i ) by being in a specific state (Y i ) across a subset of samples (M i ) and to find all corresponding passenger genes (S i ). As we test every gene in the data set, we also estimate the statistical significance (p-value) of identified context-motif C i via permutation test and multiple testing correction by Storey’s false discovery rate (FDR) [16].

Once the context-motifs are identified with statistical significance, each context is considered to manifest regulatory relationships between the driver genes and corresponding passenger genes, i.e G i g S i , specific to M i with G i (drivers) conditioned on a specific state Y i = y i . A driver g j in context-motif C j could be a passenger in another context-motif C i , conditioned by g i . When such implicit driver-passenger relationships g i g j are added together, a set of context-motifs identified from a given data set can be represented as a graph. The context-motif-specific gene-gene interactions represented in a graph can be further analyzed as described below to reveal context-specific gene regulatory network.

Contexts and context-specific gene regulatory networks

The graph described above consists of several hundreds (or thousands) of context-motifs and thousands of gene interactions, and each interaction is specific to certain subset of samples. Hence, this graph might be sub-divided into sub-networks based on its topological structure, and each sub-network might be associated with subset of samples. We utilize a clustering technique for graph, specifically, Markov clustering, as described in Ramesh et al. [8, 15].

Markov clustering (MCL) is an unsupervised graph clustering algorithm that simulates the flow in a graph using the notion of random walks. If a random walk visits a node in a cluster, it would be likely to visit several other members of the cluster before leaving the cluster [8].

The algorithm consists of two alternating operations; expansion and inflation to simulate the flow. Graph expansion is identical to taking the power of a matrix using matrix multiplication, which homogenizes the flow across different regions of the graph. The second operation, inflation, is mathematically equivalent to taking the Hadamard power of a matrix followed by scaling. Simply, the graph is denoted by a matrix of transition probabilities and expansion computes random walks by assigning probabilities with all pairs of nodes, since there are more paths within a cluster than between clusters the probabilities will be higher within a cluster. To maintain the stochastic property of the matrix, inflation re-scales the columns. Thus, the inflation parameter controls the granularity of the clusters. We use an implementation of Markov clustering based on the algorithm proposed by van Dongen [17].

Sample-Context association

Contexts obtained from clustering consist of quite a few context-motifs each of which is individually represented by a set of variables (genes) and conditions (samples). We developed a method to aggregate all the samples assigned to the context-motifs in a context and to determine if a sample can be specifically associated with the context with statistical significance.

Formally, let N be the number of samples and k i the number of samples in a context motif C i . Now let C be a context made of {C1, C2,…, C m }. In a simple approach, the samples for the context cluster can be assigned by combining all the samples in every context-motifs:

(3)

However, some samples could be present in only one or two context-motifs and may not represent the overall context. Hence, we use a metric to evaluate samples that are consistently present across majority of the context-motifs to systematically associate samples to context. Let C{j} C denote the subset of C in which the sample s j is included. Then, we define a likelihood that sample s j belongs to C, considering the fact that each context motif C i consists of different number of samples, as:

(4)

where s j C indicates s j is assigned to C, and

(5)

to compensate the different sample size associated with each context motif. It’s easy to see 0 ≤ L(s j C) ≤ 1, where L(s j C) = 0 indicates no appearance of the sample in any context motif, while L(s j C) = 1 indicates the presence of the sample in every context motif. K is used to control how favorably one wants to consider context-specificity of sample membership to a given context. The higher the K, the more context-specific the sample membership is.

Enrichment analysis

Intrinsic subtypes of breast cancer

A method, Single Subtype Predictor (SSP), for individual class classification developed by Hu et al. [18] was used to classify tumors from the 1,614 samples into five intrinsic subtypes. The algorithm uses the expression of 306 “intrinsic genes” across 315 samples of known subtypes to define a “centroid” (expression profile) for each subtype (available at https://genome.unc.edu/pubsup/breastTumor/). New tumors are then classified based on the expression profile of these 306 genes, with tumors assigned to the closest subtype centroid using Spearman rank correlation as a measure of distance. Probe sets from the Affymetrix data sets used here were mapped to the 306 genes in the intrinsic gene set, with median log base 2 intensities used when multiple probe sets matched a gene in the “intrinsic” list. The log-transformed expression data for each gene was then mean-centered within each cohort, before comparing them to the subtype centroid for classification.

Phenotype enrichment

Subsequent to clustering of contexts and associating samples to contexts, we study the phenotypic characteristics of each context. We use the intrinsic subtypes, as described above, such as Estrogen receptor (ER) status, Progesterone receptor (PgR) status, lymph node (LN) status and grade of the tumor, as phenotypes. Each of the phenotype determines certain characteristics of the tumor and can reveal therapeutic treatment options. Tumors contexts enriched with these phenotypes can provide interesting biological insights. Enrichment of contexts with a certain phenotype can be performed using hyper-geometric probability with multiple testing correction [16].

Functional annotation: gene set enrichment analysis

In addition to the phenotypic enrichments of a contexts, we also investigate the enrichment of biological functions in each context, using gene set associated with each context. The Molecular Signatures Database (MSigDB) consists of collections of gene sets such as Gene Ontology (GO) gene sets, gene sets for Biological Processes, pathway gene sets, curated sets, and computationally predicted gene expression neighborhoods underlying certain biological characteristics [19]. Genes can be annotated using a method called gene set enrichment analysis, which computes the enrichment of database gene sets with the genes found in the contexts. This method also uses hypergeometric test to measure the significance of the enrichment. A gene annotation tool GATHER was also used in for annotation of contexts (http://gather.genome.duke.edu/). The overall process of mining context-motifs followed by chaining context-motifs to obtain contexts can be illustrated in Figure 1. The process flow diagram also illustrates functional annotation processes for genes within the contexts and phenotype enrichment for samples belonging to each context.

References

  1. Ivshina AV, George J, Senko O, Mow B, Putti TC, Smeds J, Lindahl T, Pawitan Y, Hall P, Nordgren H, Wong JEL, Liu ET, Bergh J, Kuznetsov VA, Miller LD: 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

    Article  CAS  PubMed  Google Scholar 

  2. Fan C, Oh DS, Wessels L, Weigelt B, Nuyten DSA, Nobel AB, van’t Veer LJ, Perou CM: Concordance among gene-expression-based predictors for breast cancer. N Engl J Med 2006, 355(6):560–569. 10.1056/NEJMoa052933

    Article  CAS  PubMed  Google Scholar 

  3. Gatza ML, Lucas JE, Barry WT, Kim JW, Wang Q, Crawford MD, Datto MB, Kelley M, Mathey-Prevot B, Potti A, Nevins JR: A pathway-based classification of human breast cancer. Proc Natl Acad Sci U S A 2010, 107(15):6994–6999. 10.1073/pnas.0912708107

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Prat A, Karginova O, Fan C, Perou CM: Notch-associated expression profiles in basal-like and claudin-low breast cancer molecular subtypes. J Clin Oncol (Meeting Abstracts) 2009, 27(15S):11017. [http://meeting.ascopubs.org/cgi/content/abstract/27/15S/11017]

    Google Scholar 

  5. Prat A, Parker J, Karginova O, Fan C, Livasy C, Herschkowitz J, He X, Perou C: Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer. Breast Cancer Res 2010, 12(5):R68. 10.1186/bcr2635

    Article  PubMed Central  PubMed  Google Scholar 

  6. Dougherty ER, Brun M, Trent JM, Bittner ML: Conditioning-based modeling of contextual genomic regulation. IEEE/ACM Trans Comput Biol Bioinform 2009, 6(2):310–320. 10.1109/TCBB.2007.70247

    Article  CAS  PubMed  Google Scholar 

  7. Sen I, Verdicchio M, Jung S, Trevino R, Bittner M, Kim S: Context-Specific Gene Regulations in Cancer Gene Expression Data. Pacific Symposium on Biocomputing 2009, 14: 75–86.

    Google Scholar 

  8. Dongen SV: Graph Clustering by Flow Simulation. PhD thesis. University of Utrecht; 2000.

    Google Scholar 

  9. Jaccard P: Distribution de la flore alpine dans le bassin des Dranses et dans quelques régions voisines. Bulletin de la Société Vaudoise des Sciences Naturelles 1901, 37: 241–272.

    Google Scholar 

  10. Ooka M, Sakita I, Fujiwara Y, Tamaki Y, Yamamoto H, Aihara T, Miyazaki M, Kadota M, Masuda N, Sugita Y, Iwao K, Monden M: Selection of mRNA markers for detection of lymph node micrometastases in breast cancer patients. Oncol Rep 2000, 7(3):561–566.

    CAS  PubMed  Google Scholar 

  11. Meunier-Carpentier S, Dales JP, Djemli A, Garcia S, Bonnier P, Andrac-Meyer L, Lavaut MN, Allasia C, Charpin C: Comparison of the prognosis indication of VEGFR-1 and VEGFR-2 and Tie2 receptor expression in breast carcinoma. Int J Oncol 2005, 26(4):977–984.

    CAS  PubMed  Google Scholar 

  12. Placanica L, Chien JW, Li YM: Characterization of an atypical gamma-secretase complex from hematopoietic origin. Biochemistry 2010, 49(13):2796–2804. 10.1021/bi901388t

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Yin L, Velazquez OC, Liu ZJ: Notch signaling: emerging molecular targets for cancer therapy. Biochem Pharmacol 2010, 80(5):690–701. 10.1016/j.bcp.2010.03.026

    Article  CAS  PubMed  Google Scholar 

  14. Kim S, Sen I, Bittner ML: Mining molecular contexts of cancer via in-silico conditioning. In Comput Syst Bioinformatics. Volume 6. World Scientific Publishing; 2007:169–179.

    Chapter  Google Scholar 

  15. Ramesh A, Trevino R, VON-Hoff DD, Kim S: Clustering context-specific gene regulatory networks. Pac Symp Biocomput 2010, 444–455.

    Google Scholar 

  16. Storey JD: A direct approach to false discovery rates. Journal of the Royal Statistical Society 2002, Series B(64):479–498.

    Article  Google Scholar 

  17. van Dongen S: A Cluster Algorithm for Graphs. In Tech. Rep. INS-R0010. National Research Institute for Mathematics and Computer Science; 2000.

    Google Scholar 

  18. Hu Z, Fan C, Oh DS, Marron JS, He X, Qaqish BF, Livasy C, Carey LA, Reynolds E, Dressler L, Nobel A, Parker J, Ewend MG, Sawyer LR, Wu J, Liu Y, Nanda R, Tretiakova M, Ruiz Orrico A, Dreher D, Palazzo JP, Perreard L, Nelson E, Mone M, Hansen H, Mullins M, Quackenbush JF, Ellis MJ, Olopade OI, Bernard PS, Perou CM: The molecular portraits of breast tumors are conserved across microarray platforms. BMC Genomics 2006, 7: 96. 10.1186/1471-2164-7-96

    Article  PubMed Central  PubMed  Google Scholar 

  19. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 2005, 102(43):15545–15550. 10.1073/pnas.0506580102

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Yoon NK, Maresh EL, Shen D, Elshimali Y, Apple S, Horvath S, Mah V, Bose S, Chia D, Chang HR, Goodglick L: Higher levels of GATA3 predict better survival in women with breast cancer. Hum Pathol 2010, 41(12):1794–1801. 10.1016/j.humpath.2010.06.010

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Miller LD, Smeds J, George J, Vega VB, Vergara L, Ploner A, Pawitan Y, Hall P, Klaar S, Liu ET, Bergh J: An expression signature for p53 status in human breast cancer predicts mutation status, transcriptional effects, and patient survival. Proc Natl Acad Sci U S A 2005, 102(38):13550–13555. 10.1073/pnas.0506230102

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Sotiriou C, Wirapati P, Loi S, Harris A, Fox S, Smeds J, Nordgren H, Farmer P, Praz V, Haibe-Kains B, Desmedt C, Larsimont D, Cardoso F, Peterse H, Nuyten D, Buyse M, Van de Vijver MJ, Bergh J, Piccart M, Delorenzi M: 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. 10.1093/jnci/djj052

    Article  CAS  PubMed  Google Scholar 

  23. Pawitan Y, Bjohle J, Amler L, Borg AL, Egyhazi S, Hall P, Han X, Holmberg L, Huang F, Klaar S, Liu ET, Miller L, Nordgren H, Ploner A, Sandelin K, Shaw PM, Smeds J, Skoog L, Wedren S, Bergh J: 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–64. 10.1186/bcr1325

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Desmedt C, Piette F, Loi S, Wang Y, Lallemand F, Haibe-Kains B, Viale G, Delorenzi M, Zhang Y, d’Assignies MS, Bergh J, Lidereau R, Ellis P, Harris AL, Klijn JGM, Foekens JA, Cardoso F, Piccart MJ, Buyse M, Sotiriou C: Strong time dependence of the 76-gene prognostic signature for node-negative breast cancer patients in the TRANSBIG multicenter independent validation series. Clin Cancer Res 2007, 13(11):3207–3214. 10.1158/1078-0432.CCR-06-2765

    Article  CAS  PubMed  Google Scholar 

  25. Schmidt M, Bohm D, von Torne C, Steiner E, Puhl A, Pilch H, Lehr HA, Hengstler JG, Kolbl H, Gehrmann M: The humoral immune system has a key prognostic impact in node-negative breast cancer. Cancer Res 2008, 68(13):5405–5413. 10.1158/0008-5472.CAN-07-5206

    Article  CAS  PubMed  Google Scholar 

  26. Zhang Y, Sieuwerts AM, McGreevy M, Casey G, Cufer T, Paradiso A, Harbeck N, Span PN, Hicks DG, Crowe J, Tubbs RR, Budd GT, Lyons J, Sweep FCGJ, Schmitt M, Schittulli F, Golouh R, Talantov D, Wang Y, Foekens JA: The 76-gene signature defines high-risk patients that benefit from adjuvant tamoxifen therapy. Breast Cancer Res Treat 2009, 116(2):303–309. 10.1007/s10549-008-0183-2

    Article  CAS  PubMed  Google Scholar 

  27. 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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Minn AJ, Gupta GP, Padua D, Bos P, Nguyen DX, Nuyten D, Kreike B, Zhang Y, Wang Y, Ishwaran H, Foekens JA, van de Vijver M, Massague J: Lung metastasis genes couple breast tumor size and metastatic spread. Proc Natl Acad Sci U S A 2007, 104(16):6740–6745. 10.1073/pnas.0701138104

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Carroll JS, Meyer CA, Song J, Li W, Geistlinger TR, Eeckhoute J, Brodsky AS, Keeton EK, Fertuck KC, Hall GF, Wang Q, Bekiranov S, Sementchenko V, Fox EA, Silver PA, Gingeras TR, Liu XS, Brown M: Genome-wide analysis of estrogen receptor binding sites. Nat Genet 2006, 38(11):1289–1297. 10.1038/ng1901

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the Computational Biology Division team at TGen and the Computational Systems Biology group at Arizona State University. Our special thanks to Sungwon Jung, who developed some of the tools used in the analysis. SK is partially supported by NIH 1R21LM009706-01, SFAZ CAA 0243-08, and NIH P01 CA109552-01A1. SN is partially supported by SFAZ CAA 0243-08. We also thank the reviewers for DTMBIO for their valuable suggestions that helped us improve our manuscript.

This article has been published as part of BMC Bioinformatics Volume 12 Supplement 2, 2011: Fourth International Workshop on Data and Text Mining in Bioinformatics (DTMBio) 2010. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/12?issue=S2.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Seungchan Kim.

Additional information

Authors' contributions

SN and SK participated in the design of the study. SN performed the data preparation and analysis. MB did the classification of the tumors into intrinsic subtypes. HC did the biological evaluation and wrote discussion section. The draft was initially prepared by SN and SK and was reviewed by HC and MB. All authors reviewed the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material

12859_2011_4451_MOESM1_ESM.pdf

Additional file 1: Functional annotation for contexts Functional annotation for selected contexts is provided as Supplement tables 1-7. Each table lists pathways or gene sets found to be enriched with genes from a context, size of the pathway or gene set, its description, amount of overlap and statistical significance. (PDF 47 KB)

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Nasser, S., Cunliffe, H.E., Black, M.A. et al. Context-specific gene regulatory networks subdivide intrinsic subtypes of breast cancer. BMC Bioinformatics 12 (Suppl 2), S3 (2011). https://doi.org/10.1186/1471-2105-12-S2-S3

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-12-S2-S3

Keywords