DawnRank: discovering personalized driver genes in cancer
© Hou and Ma; licensee BioMed Central 2014
Received: 20 April 2014
Accepted: 18 July 2014
Published: 31 July 2014
Large-scale cancer genomic studies have revealed that the genetic heterogeneity of the same type of cancer is greater than previously thought. A key question in cancer genomics is the identification of driver genes. Although existing methods have identified many common drivers, it remains challenging to predict personalized drivers to assess rare and even patient-specific mutations. We developed a new algorithm called DawnRank to directly prioritize altered genes on a single patient level. Applications to TCGA datasets demonstrated the effectiveness of our method. We believe DawnRank complements existing driver identification methods and will help us discover personalized causal mutations that would otherwise be obscured by tumor heterogeneity. Source code can be accessed at http://bioen-compbio.bioen.illinois.edu/DawnRank/.
Recent advances in next-generation sequencing (NGS) technologies have provided us with an unprecedented opportunity to better characterize the molecular signatures of human cancers. The critical challenge facing cancer genomics today is to analyze and integrate such information in the most efficient and meaningful way to advance cancer biology, and then to translate that knowledge to clinic ,. A key question in cancer genomics is how to distinguish ‘driver’ mutations, which contribute to tumorigenesis, from functionally neutral ‘passenger’ mutations . The most basic approach is to categorize mutations based on recurrence, that is, the most commonly occurring mutations are more likely to be drivers ,, or by comparing mutation rates in individual genes based on an empirically derived background mutation rate, such as MutSig  and MuSiC . Machine learning based approaches use existing knowledge to help identify drivers. For example, CHASM utilizes random forest to classify driver mutations using alterations trained from known cancer-causing somatic missense mutations . There are several recent methods that use additional information to help predict driver genes and driver pathways. CONEXIC was developed to integrate copy number change and gene expression data to identify potential driver genes located in regions that are amplified or deleted in tumors . Network and pathway-based approaches have become one of the most promising methods to understand drivers due to their ability to model gene-gene interactions by aggregating small effect sizes from individual genes. MEMo and Dendrix rely on predicted mutual exclusivity of driver mutations within pathways or subnetworks ,. MEMo utilizes driver cliques based on known pathways with mutually exclusive mutations in the patient cohort, whereas Dendrix identifies subnetworks (de novo) with potential driver activity as having high coverage and high mutual exclusivity. Another method, PARADIGM-Shift was developed to utilize pathway-level information along with other features (such as expression, methylation, copy number) to infer gain and loss of function for mutations. PARADIGM-Shift works best with small pathways with specific genes of interest, that is, focus genes . DriverNet classifies driver mutations as mutations that propagate outlying downstream differential expression in the transcriptional regulatory network . More recently, TieDIE  was developed to find small cancer driver pathways within a large super-pathway to connect genomic alterations with transcriptomic changes, MAXDRIVER  was proposed to identify driver genes by integrating multiple omics data and heterogeneous networks, and VarWalker  was developed to construct cancer-specific mutation networks to assist driver identification in a patient cohort. However, the scope of these methods is still limited. Most existing methods such as MEMo, Dendrix, and DriverNet require a large number of patient samples to generate reliable results and lack the ability to discover rare and patient-specific drivers. Other methods such as PARADIGM-Shift are designed to determine drivers in small pathways and often require detailed previous knowledge of specific pathways and focus genes to operate effectively. New methods are needed to identify novel and rare drivers when we do not have much prior knowledge of the tumor.
It is now acknowledged that individual tumors of the same type are highly heterogeneous and have diverse genomic alterations ,. This stems from the ‘long-tail phenomenon’ which states that cancer mutations are characterized by a small number of frequently mutated genes and a large number of infrequently mutated genes ,. Discovering rare drivers in the long tail of genetic alterations remains difficult. Therefore, we urgently need methods to assess the impact of patient-specific and rare mutations from individual tumor samples in order to elucidate personalized molecular drivers.
In this work, we introduce a new method called DawnRank that detects driver genes using data from a single patient sample. By only using data from an individual patient sample rather than a large cohort, we identify drivers in a personalized fashion. The single patient approach classifies drivers regardless of mutation frequency, thereby allowing us to focus on rare (infrequent) drivers. DawnRank ranks potential driver genes based on their impact on the overall differential expression of its downstream genes in the molecular interaction network. Mutated genes with higher rankings are more likely to be drivers. Our method builds on the DriverNet rationale that the impact of a potential driver gene can be determined by its effect on the genes that are regulated by it. However, unlike DriverNet, DawnRank can be applied to a single patient at a time. It is also important to note that our approach differs from the aforementioned PARADIGM-Shift. The small-scale network approach of PARADIGM-Shift works well when the user has a target pathway and genes in mind, but is not applicable in determining rare driver genes distributed across multiple pathways with a more comprehensive gene network. DawnRank differs from MEMo and Dendrix as the scope of DawnRank is to determine individual driver genes rather than driver pathways. Although CONEXIC also utilizes pathway information to identify important CNVs, DawnRank has a more generalized framework to assess the impact of both mutations and CNVs at the individual patient level. Such a patient-specific framework also makes DawnRank different from methods such as MAXDRIVER and VarWalker. DawnRank also differs from TieDIE, which predicts subnetworks of interlinking genes that highlight cancer specific subnetworks, whereas DawnRank predicts individual driver genes.
The DawnRank algorithm
In DawnRank, a gene will possess a higher impact score (that is, rank) if the gene is highly connected to differentially expressed downstream genes (directly and indirectly connected). Driver genes tend to display a high-degree of connectivity within the gene network ,. For example, using the number of outgoing edges alone, known driver genes as classified by the Cancer Gene Census (CGC)  have a mean and median of 31.45 and 12 outgoing edges, respectively, whereas genes not typically classified as drivers (not in CGC) have a mean and median of 17.73 and three outgoing edges, respectively. The higher number of outgoing connectivity of known driver genes suggests that the PageRank model would be appropriate to prioritize driver genes based on their impact in the gene interaction network. PageRank has had several adaptations in genomics. GeneRank utilized PageRank to rank the importance of genes in a molecular network . PageRank derivatives (such as SPIA ) have also been used to analyze pathway-level importance. More recently, it was utilized to predict clinical outcome of cancer patients based on gene expression  and to assist subtype identification . Such approaches also show similarity in nature to modeling network impact as a heat diffusion process as used in HotNet  and TieDIE . DawnRank builds on the original PageRank algorithm by providing a way to model a network’s directionality with more stable rankings by utilizing dynamic damping factors (see below).
DawnRank views the gene network as a directed graph. Let N be the number of nodes (in our case, genes) in the directed graph, and A be the adjacency matrix representation of the graph, a 0-1 matrix (if node i links to j, then A ij =1). Note that the current 0-1 adjacency matrix can be naturally extended to consider weighted edges to further distinguish different gene-gene interactions.
The parameter μ follows a Dirichlet prior trained from maximizing the values of μ over 100 random samples. We selected the μ value of 3 because it had the highest average DawnRank scores for known drivers in CGC (see Additional file 1: Supplementary Materials for more details on how μ was trained). Overall, the dynamic damping factor mitigates the large change in the damping factor in nodes with 0 and 1 incoming edges by gradually increasing the damping factor as the gene’s in-degree increases, thereby creating more reliable and more stable rankings. A toy example is in the Additional file 1: Figure S2. We also show that DawnRank performs more reliably with a dynamic damping factor than a static damping factor on the TCGA datasets (see Additional file 1: Figure S1).
DawnRank converges when there is no longer a significant update in the ranks. This is when the magnitude of the difference of the ranks between time t + 1 and the previous time point t falls below a small ε, which we set to 0.001, the same value suggested by . DawnRank also stops when no solution is present after a maximum number of iterations, which we set at 100. In practice, DawnRank always converges for any reasonable μ between 0.01 and 20 within 20 iterations. Nonetheless, there are corner cases at low damping factors (μ < 10− 10) where DawnRank either does not converge or converges very slowly. A detailed discussion of the convergence of DawnRank can be found in the Additional file 1: Supplementary Materials.
Condorcet voting for rank aggregation
We used the output from DawnRank, which we converted to percentile rank format, to represent the ranking of the gene. The penalty heuristic lowers the value of a non-mutated gene when comparing it against a mutated gene. This heuristic serves as both a means to prevent a rare mutation that is impactful in one patient from winning all pairwise comparisons (akin to a candidate winning just because one and only voter that voted for it ranked it higher than any other candidate) and to prevent a low-impact, high-frequency mutation from winning a pairwise comparison against high-impact genes that are not frequently mutated (akin to an unpopular candidate winning just because many voters had a low-preference vote for that candidate). We selected δ by running DawnRank over 100 random patient samples for various instances of δ between 0 and 1 and calculating the precision with respect to CGC genes (see Additional file 1: Supplementary Materials). We found δ to be 0.85.
Results and discussion
We applied DawnRank to TCGA datasets. We first showed that DawnRank outperforms two recent pathway-based methods DriverNet and PARADIGM-Shift. We also showed that DawnRank produces reliable results as compared to cohort-based approaches CHASM and OncodriveFM. We then used the results of DawnRank to determine both potential novel drivers (new genes mutated frequently), and more importantly, potential rare and personalized drivers that previously could not be assessed by other methods. We also applied the method to study breast cancer subtypes and found that the amount of predicted rare drivers has a strong correlation with the degree of genetic heterogeneity in different breast cancer subtypes.
Datasets and network
We applied DawnRank to 512 glioblastoma multiforme (GBM) samples, 504 breast cancer (BRCA) samples, and 572 ovarian cancer (OV) samples in TCGA. The datasets we used in this work include gene expression and coding-region mutation data for three cancer types generated by TCGA –. The data were accessed on 20 May 2013. The mutation data we used included non-synonymous point mutations and insertions and deletions (indels) in coding regions.
We built the gene interaction network used in DawnRank using a variety of sources, including the network used in MEMo  as well as the up-to-date curated information from Reactome , the NCI-Nature Curated PID , and KEGG . The MEMo network consisted of inferred gene-interaction from sources of information such as protein interactions, gene co-expression, protein domain interaction, and text-mined interaction described by . To aggregate all of the networks together, all redundant edges were collapsed to single edges. We included self-loops within the network to account for auto-regulation events . The resulting aggregate network consisted of 11,648 genes and 211,794 edges, and can be downloaded from our supplementary website.
To help evaluate the quality of our results, we obtained a list of 487 known driver genes from the well-studied cancer gene database, CGC . Additionally, we also compared the quality of our results to the Pan-Cancer drivers , a driver gene list built using results from well known cohort-based methods over twelve tumor types. In practice, there is no gold standard of known drivers. However, well-curated cancer gene databases such as CGC and the Pan-Cancer study provide an approximate benchmark of known drivers ,.
DawnRank outperforms pathway-based methods DriverNet and PARADIGM-Shift
DawnRank achieves reliable results as compared to cohort-based approaches CHASM and OncodriveFM
In addition to comparing DawnRank to the pathway-based approaches of DriverNet and PARADIGM-Shift, we compared DawnRank to non-pathway-based approaches as well. We compared DawnRank with CHASM  and OncodriveFM . We ran CHASM version 3.0 using CRAVAT  with the cancer driver analysis mode , and we ran OncodriveFM using the IntOGen Mutation Analysis interface using SIFT, Polyphen2, and Mutation Assessor data . Both CHASM and OncodriveFM can rank genes from most likely drivers to least likely ones. Like in DawnRank, we evaluated all non-synonymous mutations and evaluated each method using precision-recall metrics in comparison to CGC. The results (Additional file 1: Figure S10) show that DawnRank performs comparably to both CHASM and OncodriveFM. These results also demonstrate that the personalized approach of DawnRank can achieve reliable results at the cohort level.
DawnRank can discover both novel and personalized drivers
In this section, we use DawnRank to identify driver genes that may not have been classified as drivers by other methods. We discuss such driver genes determined by DawnRank in two main categories: general novel drivers, and personalized drivers (we are mostly interested in personalized rare drivers).
A general novel driver is a recurrent driver that has not been classified as a driver gene. Here we define novel drivers as genes that satisfy the following requirements: (1) Highly ranked (genes that have an aggregate rank in the top 30 of driver genes based on the patient cohort), which means it is considered as driver in multiple samples; (2) Recurrent (alteration frequency >2% in the patient cohort); (3) Not previously classified as a driver by CGC.
A personalized driver is a gene predicted to be a driver for specific patients. Here we define personalized drivers as drivers classified as significant using DawnRank’s maximally selected rank statistics cutoff (explained later) and as drivers that display higher than expected rankings under the Chauvenet’s criterion for detecting outliers. In other words, personalized driver genes could be recurrently mutated in multiple samples, but they are only ranked high enough to be considered as drivers in specific samples. We are particularly interested in rare drivers, as rare drivers would be obscured by the long tail phenomenon of cancer mutations. Personalized drivers mutated in less than 2% of patients will be considered personalized rare drivers.
Discovering novel drivers with coding mutations
Using the criteria listed above to detect novel drivers, we used the aggregate DawnRank score to select genes with coding mutations that were considered novel drivers. In GBM, we found several consistently high-ranking potential novel driver genes. Of them, TGFBR2, HIF1A, and FOXO3 are the most promising. These genes are highly ranked, frequently mutated, and are involved in several cancer functions and pathways (see Additional file 3). TGFBR2 is a transforming growth factor whose function is to regulate cell division signals and inhibit cell division . TGFBR2 is the third-ranked driver in terms of Condorcet aggregate score, and it is altered in 18.7% of GBM patient samples. TGFBR2 is classified as a potential tumor suppressor (it shows decreased expression in GBM samples), and it is involved in many pathways to cancer including TGFβ signaling, HTLV infection, and Hippo signaling . HIF1A is a hypoxia inducing factor that has been shown to be necessary for the survival of cancer stem cells in gliomas . HIF1A is involved in mTOR pathway, which regulates nutrient sensing for cell growth, which could be responsible for cell growth and metastasis in cancer. HIF1A is ranked as the fourth most impactful driver in GBM, and it is mutated in 13.4% of GBM patient samples . FOXO3 is the eighth-ranked gene. It is a common phosphorylation target of AKT and ERK, and is a key trigger for apoptosis in the PIK3/AKT pathway .
In OV and BRCA, there were fewer candidate novel drivers than GBM, however, we found two genes that show strong potential to become novel drivers: PDPK1 in OV and CENPE in BRCA. PDPK1 is a phosphoinositide dependent protein kinase-1 that interacts with a crucial driver gene AKT1, and together their signaling plays a critical role in activating proliferation and survival pathways within cancer cells such as the PIK3/AKT pathway and the mTOR pathway ,. PDPK1 is ranked second in OV and is mutated in 5.4% of OV samples. The mutation in PDPK1 itself also suggests an important impact in the functionality of the protein binding. For example, in sample TCGA-13-0751, the mutation is in a Glycine to Arginine substitution near a hairpin loop in the middle of the protein (Additional file 1: Figure S6). The amino acid change from G to R changes the dynamic of the protein structure by replacing a non-polar amino acid with a positively charged amino acid, which may change the binding interaction with the negatively charged substrate phosphatidylinositol-3,4,5-trisphosphate, which in turn may affect the phosphorylation and activation of at least 24 kinases .
In BRCA, a potential novel driver we found is CENPE. CENPE is a centromere-associated protein that contributes to mechanics of microtubule-chromosome interactions with mitotic checkpoint. Overexpression of CENPE can lead to excessive cell growth and contribute to tumor proliferation . CENPE is ranked 26th in DawnRank. We observed that it is highly overexpressed in tumor samples. CENPE is interesting as it is ranked differently among different subtypes of breast cancer. All but one of the mutations in CENPE was in Luminal B or Basal subtype. The Basal subtype mutations had the highest average ranking in the 95.2 percentile and the average Luminal B subtype mutations ranked lower in the 90.9 percentile, whereas the Luminal A mutation was much less significant with a ranking in the 78.3 percentile. This shows that even though CENPE has mutations in multiple subtypes, it may harbor different driver potentials in different subtypes.
Discovering novel drivers with copy number changes
We then applied DawnRank to examine the role of CNVs in the three cancer types used in this study. We included genes that displayed at least a two-fold change in copy number change (with respect to the normal baseline) and also had a corresponding change in expression. We treated CNV events in a similar fashion as coding mutations to prioritize the genes with an amplification and deletion and then determined the most impactful CNVs. In addition to known driver genes with CNVs, we also identified novel driver CNVs not in CGC (Additional file 1: Figure S7).
In GBM, known CNV drivers (based on CGC) include the amplification of some well-known driver genes such as EGFR, PDGFRA, and MYC and the deletions of MLLT3 and ANK3. A potential novel driver CNV is the amplification of SEC61G. SEC61G is a proto-oncogene that is required for tumor cell survival . SEC61G is altered in 14.3% of GBM cases and it is ranked 17th among CNVs. ELAVL2 is a potential tumor suppressor gene whose function is to stabilize and control nervous-specific binding proteins . ELAVL2 is altered in 6.7% of GBM cases and is ranked second among CNVs in GBM. In OV, many known drivers are amplifications such as CCNE1 and KRAS and deletions such as MAF or PIK3R1. An example of a novel amplification is FZD5, a ‘frizzled’ gene that has a strong co-expression event with the Wnt signaling pathway in ovarian cancer . FZD5 is altered in 2.9% of OV cases and it is ranked ninth among CNVs in OV. In BRCA, known driver CNVs include amplifications in CCND1, MYC, GATA3, and EGFR and deletions in MAF. An example of a novel amplification event in BRCA is PAK1, which is a breast cancer oncogene that activates MAPK and MET signaling pathways . PAK1 is the third ranked gene with CNVs, and it is altered in 2.4% of the BRCA patient samples.
Discovering personalized drivers
In this section, we discuss the personalized scope of DawnRank by demonstrating its ability to determine personalized novel and rare drivers. The main aspect that distinguishes DawnRank from existing programs is the ability to discover rare and even patient-specific drivers. Even if a gene is altered only in a single patient, DawnRank will still be able to evaluate the impact of that alteration. We determined such personalized drivers using the ‘maximally rank statistics’ method. This method classifies alterations as ‘drivers’ and ‘non-drivers’ by assigning a cutoff that maximizes the number of known driver genes (genes in CGC) ranking higher than the cutoff). Altered genes in the individual patient with a ranking higher than the cutoff would be described as personalized drivers. The average cutoff for GBM, OV, and BRCA was 98.1, 92.8, and 92.2, respectively (measured in percentile rank). We specifically looked for genes that have not been classified as a driver in CGC and genes that have a significantly higher than expected rank in specific patients. A rank is considered significantly high in a patient if the rank is considered as an outlier under Chauvenet’s criterion for outlier detection.
We found 26 potential personalized drivers in GBM. Of these genes, we found that several of the candidate driver genes were involved in important known cancer pathways. Using KEGG to map genes to pathways (see Additional file 3) , mutations in three interferon receptor genes IFNA14, IFNW1, and IFNA17 belong to multiple pathways that have significant impact on cancer, including the Cytokine-Cytokine receptor pathway and the JAK/STAT pathway. Cytokines are important factors in tumor cell control as they are key players in inflammatory and immune response . JAK/STAT is a related pathway activated by the binding of cytokines to its receptor. The pathway is involved in inflammation, proliferation, and invasion/migration . In either case, IFNA14, IFNW1, and IFNA17 are interferon receptor genes that may interact with genes such as TP53 to induce apoptosis of cancer cells . Both the Cytokine-Cytokine receptor pathway and the JAK/STAT pathway are common drug targets in GBM ,, which could lead to the implications that IFNA14, IFNW1, and IFNA17 could be targeted. IFNA14, IFNW1, and IFNA17 are mutated in 4.7%, 5.6%, and 4.5% of GBM patients, respectively. IFNA14 is ranked significantly higher than its average expected ranking in patients TCGA-32-1978 and TCGA-19-2619 while IFNW1 and IFNA17 both have high-ranking mutations in the patient TCGA-28-1749. Although IFNA14, IFNW1, and IFNA17 are examples of personalized drivers, they are not personalized rare drivers. An example of a personalized rare driver is PITX2, which is altered in only two patients (0.4%). PITX2 is in the TGFβ signaling pathway, which can be targeted by cancer drugs such as decanoyl-RVKR chloromethylketone in GBM ,.
We found 77 potential personalized driver genes in OV, and like GBM, many of these genes fall under similar KEGG pathways (see Additional file 3). Figure 5A shows the list of personalized drivers in OV. The most well represented KEGG pathways are the calcium signaling pathway and the MAPK signaling pathway. Calcium signaling is important in the regulation of cancer cell proliferation and/or apoptosis . In OV, four alterations are present in the calcium signaling pathway: mutations in ADORA2B, CACNA1E, and ITPR3 along with a CNV for SLC8A1. Of these alterations, ADORA2B is often associated with cancer cell movement and metastasis, and it has been a drug target in other cancers . The MAPK signaling pathway plays a role in communication with cell growth, and it can be targeted by cancer drugs such as MEK and RAF inhibitors ,. Genes mutated in the MAPK signaling pathway include CACNA1E, TRAF6, and DUSP16. The TRAF6 gene is especially interesting as it upregulates HIF1A and has implications in tumor angiogenesis . TRAF6 is another good example to demonstrate the power of DawnRank because it is mutated in only two patients, TCGA-13-1410 and TCGA-61-1919, but DawnRank was still able to detect it as a rare driver. We looked into the TRAF6 mutation further, including its connectivity in the network and the potential functional impact of the mutation (Figure 5B and C). We found that the mutation in TCGA-13-1410 is considered to have an important functional impact by MutationAssessor, affecting the MATH domain . The mutation is located at a loop defining one wall of the binding site and projects its long side-chain to directly contact the conserved aromatic residue in the peptide substrate, and it likely modulates its interaction with its receptors . TRAF6 is also known to be involved in the NF-κB pathway, which is related to poor outcomes in ovarian cancer .
We found 56 potential rare drivers in BRCA. Unlike in GBM and OV, the distribution of rare drivers in BRCA seemed to spread out among different pathway with no major cancer pathway having more than two candidate rare drivers. Nevertheless, several important cancer pathways are highlighted in BRCA (see Additional file 3). One of them is the KEGG ‘Proteoglycans in Cancer Pathway’. Proteoglycans bind to ligands and receptors that regulate tumor neoplastic growth and angiogenesis . Alterations in genes coding for proteoglycans include ANK3 and FRS2. Of these, FRS2 is especially interesting due to both its scarcity (only one amplification event in all BRCA patients) and its function as a fibroblast growth factor whose amplification activates the FGFR2. It is also highly differentially expressed in the patient, having differential expression 3.21 standard deviations above normal. The calcium signaling pathway, like in ovarian cancer, also has multiple potential drivers in breast cancer. These genes include CACNA1B and GRIN2A. GRIN2A is a potential rare driver as it is mutated in 1.7% of BRCA patients, but it has been known to be recurrent driver mutation in other cancers .
Distribution of personalized rare drivers across breast cancer subtypes
The distribution of known drivers and rare drivers may correspond to the genetic heterogeneity of breast cancer subtypes. Luminal A is an ER+ subtype with less heterogeneity and has been defined by its good prognosis. Patients with Luminal A subtype generally have the most favorable outcome and have a smaller chance of relapse when treated early ,. On the other hand, the Basal-like subtype, which mainly comprises of triple negative breast cancers, is difficult to treat due to tumor heterogeneity. Basal-like breast cancers are known to be more aggressive with poor prognosis ,. Based on our results, we found that Luminal A samples seem to have common drivers more than any other cancer type. The pattern in the Basal-like type is the opposite. There are more potential rare drivers in Basal-like samples than other subtype. This suggests that Luminal A breast cancers are driven more by common drivers while Basal-like breast cancers are driven more by rare drivers.
Potential personalized rare driver alterations for Basal-like breast cancer include LGR5 and BCL2L14. LGR5 is a CNV that is amplified and overexpressed in the Basal-like patient sample TCGA-AN-A0AT. LGR5 is necessary for the stimulation of basal cell growth in breast tissue . BCL2L14 is a key regulator of cell apoptosis at the mitochondria level , and it is altered in the single Basal-like patient TCGA-AN-A0FJ.
Implementation and running time
DawnRank was implemented in R. Using a network of 1,492 nodes and 8,070 edges and a cohort of the BRCA samples we used, we compared the run time of DawnRank, DriverNet, and PARADIGM-Shift on a computer with 16GB RAM and Intel i5-3317U processor. DriverNet was the fastest, with a runtime of 6 min. DawnRank completed the analysis in 9 min for the main script and 12 min for the Condorcet rank aggregator. PARADIGM-Shift took approximately 15 h to run.
It is now acknowledged that individual tumors of the same type are highly heterogeneous and have diverse genomic alterations. Therefore, we urgently need novel methods to identify patient-specific and rare drivers from individual tumor samples in order to elucidate personalized molecular mechanisms in different types of cancer. The goal of DawnRank is to integrate mutation data, gene expression, and network information to discover drivers in a personalized manner. We applied DawnRank to a large number of TCGA datasets. By comparing to previous studies, our results demonstrated the effectiveness of DawnRank: (1) Despite its single-patient scope, DawnRank detects common and known drivers with as much or more precision than existing methods; (2) We can identify rare and novel genes that are potential drivers to specific patients. We believe this method will complement existing driver identification methods and will help us discover potential personalized drivers. The application to breast cancer subtypes further demonstrated that the rare drivers predicted by DawnRank may provide new insights into the molecular explanations of cancer subtypes with higher tumor heterogeneity.
One potential limitation of the current DawnRank method is that we still largely rely on known molecular interaction network. Such network information is still incomplete, which may create false negatives in DawnRank’s prediction. In addition, currently the overall network information we use in DawnRank is not cancer-type specific or patient-specific. Therefore, some of the network level interactions and perturbations specific to certain cancer subtypes or patient samples may be obscured by this approach. Another potential limitation is that DawnRank detects only potential drivers that alter the expression of other genes. However, this may not be the case for all drivers. In addition, some of the predictions from our method could potentially be caused by passenger mutations that coincidentally change the expression of downstream genes. In those cases, it may be useful to additionally utilize cohort and functional-impact based methods such as OncodriveFM. As the community continues to refine our understanding of the interaction networks and cancer driver genes, we expect that DawnRank’s ability to predict driver alterations would also improve.
Our method provides a unique solution to predict potential driver genes in cancer. However, the computationally predicted personalized drivers should not be over-interpreted before additional experimental evidence becomes available. To fully validate such computational predictions, both in vitro and in vivo experimental validations would be needed to comprehensively assess the tumorigenic potentials of predicted drivers in individual patients. Nevertheless, our method provides a reliable solution to prioritize such follow-up experimental validations. Taken together, our personalized approach is promising to discover potential causal generic variants that would be otherwise obscured by tumor heterogeneity. The personalized framework may help determine the optimal treatment strategy for each patient through individualized assessment based on the molecular signatures of their cancers.
Cancer Gene Census
The Cancer Genome Atlas
We thank Benita Katzenellenbogen, Edison Liu, Koichiro Inaki, and ChengXiang Zhai for helpful discussions. We thank Sam Ng (Josh Stuart’s lab at UCSC) for technical assistance to run PARADIGM-Shift. This work was supported by an award from the Interdisciplinary Innovation Initiative (In3) program at the University of Illinois (to J.M.) and an Illinois Distinguished Fellowship (to JPH). We also thank the Cancer Genome Atlas Research Network for making the data publicly available.
- Green ED, Guyer MS: Charting a course for genomic medicine from base pairs to bedside. Nature. 2011, 470: 204-213. 10.1038/nature09764.View ArticlePubMedGoogle Scholar
- Stratton MR: Journeys into the genome of cancer cells. EMBO Mol Med. 2013, 5: 169-172. 10.1002/emmm.201202388.PubMed CentralView ArticlePubMedGoogle Scholar
- Greenman C, Stephens P, Smith R, Dalgliesh GL, Hunter C, Bignell G, Davies H, Teague J, Butler A, Stevens C, Edkins S, O'Meara S, Vastrik I, Schmidt EE, Avis T, Barthorpe S, Bhamra G, Buck G, Choudhury B, Clements J, Cole J, Dicks E, Forbes S, Gray K, Halliday K, Harrison R, Hills K, Hinton J, Jenkinson A, Jones D: Patterns of somatic mutation in human cancer genomes. Nature. 2007, 446: 153-158. 10.1038/nature05610.PubMed CentralView ArticlePubMedGoogle Scholar
- Ding L, Getz G, Wheeler DA, Mardis ER, McLellan MD, Cibulskis K, Sougnez C, Greulich H, Muzny DM, Morgan MB, Fulton L, Fulton RS, Zhang Q, Wendl MC, Lawrence MS, Larson DE, Chen K, Dooling DJ, Sabo A, Hawes AC, Shen H, Jhangiani SN, Lewis LR, Hall O, Zhu Y, Mathew T, Ren Y, Yao J, Scherer SE, Clerc K: Somatic mutations affect key pathways in lung adenocarcinoma. Nature. 2008, 455: 1069-1075. 10.1038/nature07423.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones S, Zhang X, Parsons DW, Lin JC, Leary RJ, Angenendt P, Mankoo P, Carter H, Kamiyama H, Jimeno A, Hong SM, Fu B, Lin MT, Calhoun ES, Kamiyama M, Walter K, Nikolskaya T, Nikolsky Y, Hartigan J, Smith DR, Hidalgo M, Leach SD, Klein AP, Jaffee EM, Goggins M, Maitra A, Iacobuzio-Donahue C, Eshleman JR, Kern SE, Hruban RH: Core signaling pathways in human pancreatic cancers revealed by global genomic analyses. Science. 2008, 321: 1801-1806. 10.1126/science.1164368.PubMed CentralView ArticlePubMedGoogle Scholar
- Banerji S, Cibulskis K, Rangel-Escareno C, Brown KK, Carter SL, Frederick AM, Lawrence MS, Sivachenko AY, Sougnez C, Zou L, Cortes ML, Fernandez-Lopez JC, Peng S, Ardlie KG, Auclair D, Bautista-Pina V, Duke F, Francis J, Jung J, Maffuz-Aziz A, Onofrio RC, Parkin M, Pho NH, Quintanar-Jurado V, Ramos AH, Rebollar-Vega R, Rodriguez-Cuevas S, Romero-Cordoba SL, Schumacher SE, Stransky N: Sequence analysis of mutations and translocations across breast cancer subtypes. Nature. 2012, 486: 405-409. 10.1038/nature11154.PubMed CentralView ArticlePubMedGoogle Scholar
- Dees ND, Zhang Q, Kandoth C, Wendl MC, Schierding W, Koboldt DC, Mooney TB, Callaway MB, Dooling D, Mardis ER, Wilson RK, Ding L: MuSiC: identifying mutational significance in cancer genomes. Genome Res. 2012, 22: 1589-1598. 10.1101/gr.134635.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Carter H, Chen S, Isik L, Tyekucheva S, Velculescu VE, Kinzler KW, Vogelstein B, Karchin R: Cancer-specific high-throughput annotation of somatic mutations: computational prediction of driver missense mutations. Cancer Res. 2009, 69: 6660-6667. 10.1158/0008-5472.CAN-09-1133.PubMed CentralView ArticlePubMedGoogle Scholar
- Akavia UD, Litvin O, Kim J, Sanchez-Garcia F, Kotliar D, Causton HC, Pochanard P, Mozes E, Garraway LA, Pe'er D: An integrated approach to uncover drivers of cancer. Cell. 2010, 143: 1005-1017. 10.1016/j.cell.2010.11.013.PubMed CentralView ArticlePubMedGoogle Scholar
- Ciriello G, Cerami E, Sander C, Schultz N: Mutual exclusivity analysis identifies oncogenic network modules. Genome Res. 2012, 22: 398-406. 10.1101/gr.125567.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Vandin F, Upfal E, Raphael BJ: De novo discovery of mutated driver pathways in cancer. Genome Res. 2012, 22: 375-385. 10.1101/gr.120477.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Ng S, Collisson EA, Sokolov A, Goldstein T, Gonzalez-Perez A, Lopez-Bigas N, Benz C, Haussler D, Stuart JM: PARADIGM-SHIFT predicts the function of mutations in multiple cancers using pathway impact analysis. Bioinformatics. 2012, 28: i640-i646. 10.1093/bioinformatics/bts402.PubMed CentralView ArticlePubMedGoogle Scholar
- Bashashati A, Haffari G, Ding J, Ha G, Lui K, Rosner J, Huntsman DG, Caldas C, Aparicio SA, Shah SP: DriverNet: uncovering the impact of somatic driver mutations on transcriptional networks in cancer. Genome Biol. 2012, 13: R124-10.1186/gb-2012-13-12-r124.PubMed CentralView ArticlePubMedGoogle Scholar
- Paull EO, Carlin DE, Niepel M, Sorger PK, Haussler D, Stuart JM: Discovering causal pathways linking genomic events to transcriptional states using Tied Diffusion Through Interacting Events (TieDIE). Bioinformatics. 2013, 29: 2757-2764. 10.1093/bioinformatics/btt471.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen Y, Hao J, Jiang W, He T, Zhang X, Jiang T, Jiang R: Identifying potential cancer driver genes by genomic data integration. Sci Rep. 2013, 3: 3538-PubMed CentralPubMedGoogle Scholar
- Jia P, Zhao Z: VarWalker: personalized mutation network analysis of putative cancer genes from next-generation sequencing data. PLoS Comput Biol. 2014, 10: e1003460-10.1371/journal.pcbi.1003460.PubMed CentralView ArticlePubMedGoogle Scholar
- Pe’er D, Hacohen N: Principles and strategies for developing network models in cancer. Cell. 2011, 144: 864-873. 10.1016/j.cell.2011.03.001.PubMed CentralView ArticlePubMedGoogle Scholar
- Ding L, Wendl MC, Koboldt DC, Mardis ER: Analysis of next-generation genomic data in cancer: accomplishments and challenges. Hum Mol Genet. 2010, 19: R188-R196. 10.1093/hmg/ddq391.PubMed CentralView ArticlePubMedGoogle Scholar
- Reimand J, Bader GD: Systematic analysis of somatic mutations in phosphorylation signaling predicts novel cancer drivers. Mol Syst Biol. 2013, 9: 637-10.1038/msb.2012.68.PubMed CentralView ArticlePubMedGoogle Scholar
- Page L, Brin S, Motwani R, Winograd T: The PageRank Citation Ranking: Bringing Order to the Web. Technical Report. 1999, Stanford InfoLab, Stanford, CAGoogle Scholar
- Brin S, Page L: The anatomy of a large-scale hypertextual Web search engine. Computer Networks and ISDN Systems. 1998, 30: 107-117. 10.1016/S0169-7552(98)00110-X.View ArticleGoogle Scholar
- Khurana E, Fu Y, Colonna V, Mu XJ, Kang HM, Lappalainen T, Sboner A, Lochovsky L, Chen J, Harmanci A, Das J, Abyzov A, Balasubramanian S, Beal K, Chakravarty D, Challis D, Chen Y, Clarke D, Clarke L, Cunningham F, Evani US, Flicek P, Fragoza R, Garrison E, Gibbs R, Gumus ZH, Herrero J, Kitabayashi N, Kong Y, Lage K: Integrative annotation of variants from 1092 humans: application to cancer genomics. Science. 2013, 342: 1235587-10.1126/science.1235587.PubMed CentralView ArticlePubMedGoogle Scholar
- D’Antonio M, Ciccarelli FD: Integrated analysis of recurrent properties of cancer genes to identify novel drivers. Genome Biol. 2013, 14: R52-10.1186/gb-2013-14-5-r52.PubMed CentralView ArticlePubMedGoogle Scholar
- Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, Rahman N, Stratton MR: A census of human cancer genes. Nat Rev Cancer. 2004, 4: 177-183. 10.1038/nrc1299.PubMed CentralView ArticlePubMedGoogle Scholar
- Morrison JL, Breitling R, Higham DJ, Gilbert DR: GeneRank: using search engine technology for the analysis of microarray experiments. BMC Bioinformatics. 2005, 6: 233-10.1186/1471-2105-6-233.PubMed CentralView ArticlePubMedGoogle Scholar
- Tarca AL, Draghici S, Khatri P, Hassan SS, Mittal P, Kim JS, Kim CJ, Kusanovic JP, Romero R: A novel signaling pathway impact analysis. Bioinformatics. 2009, 25: 75-82. 10.1093/bioinformatics/btn577.PubMed CentralView ArticlePubMedGoogle Scholar
- Winter C, Kristiansen G, Kersting S, Roy J, Aust D, Knosel T, Rummele P, Jahnke B, Hentrich V, Ruckert F, Niedergethmann M, Weichert W, Bahra M, Schlitt HJ, Settmacher U, Friess H, Buchler M, Saeger HD, Schroeder M, Pilarsky C, Grutzmann R: Google goes cancer: improving outcome prediction for cancer patients by network-based ranking of marker genes. PLoS Comput Biol. 2012, 8: e1002511-10.1371/journal.pcbi.1002511.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu Y, Gu Q, Hou JP, Han J, Ma J: A network-assisted co-clustering algorithm to discover cancer subtypes based on gene expression. BMC Bioinformatics. 2014, 15: 37-10.1186/1471-2105-15-37.PubMed CentralView ArticlePubMedGoogle Scholar
- Vandin F, Upfal E, Raphael BJ: Algorithms for detecting significantly mutated pathways in cancer. J Comput Biol. 2011, 18: 507-522. 10.1089/cmb.2010.0265.View ArticlePubMedGoogle Scholar
- Wang X, Tao T, Sun J-T, Shakery A, Zhai C: Dirichletrank: Solving the zero-one gap problem of pagerank. ACM Transactions on Information Systems (TOIS). 2008, 26: 10-Google Scholar
- Vinh NX, Chetty M, Coppel R, Wangikar PP: Issues impacting genetic network reverse engineering algorithm validation using small networks. Biochim Biophys Acta. 1824, 2012: 1434-1441.Google Scholar
- Pihur V, Datta S, Datta S: Finding common genes in multiple cancer types through meta-analysis of microarray experiments: a rank aggregation approach. Genomics. 2008, 92: 400-403. 10.1016/j.ygeno.2008.05.003.View ArticlePubMedGoogle Scholar
- Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature. 2008, 455: 1061-1068. 10.1038/nature07385.View ArticleGoogle Scholar
- Integrated genomic analyses of ovarian carcinoma. Nature. 2011, 474: 609-615. 10.1038/nature10166.View ArticleGoogle Scholar
- Comprehensive molecular portraits of human breast tumours. Nature. 2012, 490: 61-70. 10.1038/nature11412.View ArticleGoogle Scholar
- Croft D, O'Kelly G, Wu G, Haw R, Gillespie M, Matthews L, Caudy M, Garapati P, Gopinath G, Jassal B, Jupe S, Kalatskaya I, Mahajan S, May B, Ndegwa N, Schmidt E, Shamovsky V, Yung C, Birney E, Hermjakob H, D'Eustachio P, Stein L: Reactome: a database of reactions, pathways and biological processes. Nucleic Acids Res. 2011, 39: D691-D697. 10.1093/nar/gkq1018.PubMed CentralView ArticlePubMedGoogle Scholar
- Schaefer CF, Anthony K, Krupa S, Buchoff J, Day M, Hannay T, Buetow KH: PID: the Pathway Interaction Database. Nucleic Acids Res. 2009, 37: D674-D679. 10.1093/nar/gkn653.PubMed CentralView ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M: KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012, 40: D109-D114. 10.1093/nar/gkr988.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu G, Feng X, Stein L: A human functional protein interaction network and its application to cancer data analysis. Genome Biol. 2010, 11: R53-10.1186/gb-2010-11-5-r53.PubMed CentralView ArticlePubMedGoogle Scholar
- Becskei A, Serrano L: Engineering stability in gene networks by autoregulation. Nature. 2000, 405: 590-593. 10.1038/35014651.View ArticlePubMedGoogle Scholar
- Tamborero D, Gonzalez-Perez A, Perez-Llamas C, Deu-Pons J, Kandoth C, Reimand J, Lawrence MS, Getz G, Bader GD, Ding L, Lopez-Bigas N: Comprehensive identification of mutational cancer driver genes across 12 tumor types. Sci Rep. 2013, 3: 2650-PubMed CentralPubMedGoogle Scholar
- Lawrence MS, Stojanov P, Mermel CH, Robinson JT, Garraway LA, Golub TR, Meyerson M, Gabriel SB, Lander ES, Getz G: Discovery and saturation analysis of cancer genes across 21 tumour types. Nature. 2014, 505: 495-501. 10.1038/nature12912.PubMed CentralView ArticlePubMedGoogle Scholar
- DriverNet:, [http://www.bioconductor.org/packages/release/bioc/html/DriverNet.html]
- PARADIGM-Shift: , [http://sysbio.soe.ucsc.edu/paradigm/tutorial/]
- Gonzalez-Perez A, Lopez-Bigas N: Functional impact bias reveals cancer drivers. Nucleic Acids Res. 2012, 40: e169-10.1093/nar/gks743.PubMed CentralView ArticlePubMedGoogle Scholar
- Douville C, Carter H, Kim R, Niknafs N, Diekhans M, Stenson PD, Cooper DN, Ryan M, Karchin R: CRAVAT: cancer-related analysis of variants toolkit. Bioinformatics. 2013, 29: 647-648. 10.1093/bioinformatics/btt017.PubMed CentralView ArticlePubMedGoogle Scholar
- CRAVAT: , [http://www.cravat.us/]
- IntOGen Mutations Analysis: , [http://www.intogen.org/web/mutations/v04/analysis]
- Ye XZ, Xu SL, Xin YH, Yu SC, Ping YF, Chen L, Xiao HL, Wang B, Yi L, Wang QL, Jiang XF, Yang L, Zhang P, Qian C, Cui YH, Zhang X, Bian XW: Tumor-associated microglia/macrophages enhance the invasion of glioma stem-like cells via TGF-beta1 signaling pathway. J Immunol. 2012, 189: 444-453. 10.4049/jimmunol.1103248.View ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M: Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res . 2002, 42: D199-D205. 10.1093/nar/gkt1076.View ArticleGoogle Scholar
- Wang Y, Liu Y, Malek SN, Zheng P, Liu Y: Targeting HIF1alpha eliminates cancer stem cells in hematological malignancies. Cell Stem Cell. 2011, 8: 399-411. 10.1016/j.stem.2011.02.006.PubMed CentralView ArticlePubMedGoogle Scholar
- Hsieh AC, Liu Y, Edlind MP, Ingolia NT, Janes MR, Sher A, Shi EY, Stumpf CR, Christensen C, Bonham MJ, Wang S, Ren P, Martin M, Jessen K, Feldman ME, Weissman JS, Shokat KM, Rommel C, Ruggero D: The translational landscape of mTOR signalling steers cancer initiation and metastasis. Nature. 2012, 485: 55-61. 10.1038/nature10912.PubMed CentralView ArticlePubMedGoogle Scholar
- Sunayama J, Sato A, Matsuda K, Tachibana K, Watanabe E, Seino S, Suzuki K, Narita Y, Shibui S, Sakurada K, Kayama T, Tomiyama A, Kitanaka C: FoxO3a functions as a key integrator of cellular signals that control glioblastoma stem-like cell differentiation and tumorigenicity. Stem Cells. 2011, 29: 1327-1337.PubMedGoogle Scholar
- Meuillet EJ, Zuohe S, Lemos R, Ihle N, Kingston J, Watkins R, Moses SA, Zhang S, Du-Cuny L, Herbst R, Jacoby JJ, Zhou LL, Ahad AM, Mash EA, Kirkpatrick DL, Powis G: Molecular pharmacology and antitumor activity of PHT-427, a novel Akt/phosphatidylinositide-dependent protein kinase 1 pleckstrin homology domain inhibitor. Mol Cancer Ther. 2010, 9: 706-717. 10.1158/1535-7163.MCT-09-0985.PubMed CentralView ArticlePubMedGoogle Scholar
- Ericson K, Gan C, Cheong I, Rago C, Samuels Y, Velculescu VE, Kinzler KW, Huso DL, Vogelstein B, Papadopoulos N: Genetic inactivation of AKT1, AKT2, and PDPK1 in human colorectal cancer cells clarifies their roles in tumor growth regulation. Proc Natl Acad Sci U S A. 2010, 107: 2598-2603. 10.1073/pnas.0914018107.PubMed CentralView ArticlePubMedGoogle Scholar
- Komander D, Fairservice A, Deak M, Kular GS, Prescott AR, Peter Downes C, Safrany ST, Alessi DR, van Aalten DM: Structural insights into the regulation of PDK1 by phosphoinositides and inositol phosphates. EMBO J. 2004, 23: 3918-3928. 10.1038/sj.emboj.7600379.PubMed CentralView ArticlePubMedGoogle Scholar
- Wood KW, Chua P, Sutton D, Jackson JR: Centromere-associated protein E: a motor that puts the brakes on the mitotic checkpoint. Clin Cancer Res. 2008, 14: 7588-7592. 10.1158/1078-0432.CCR-07-4443.View ArticlePubMedGoogle Scholar
- Lu Z, Zhou L, Killela P, Rasheed AB, Di C, Poe WE, McLendon RE, Bigner DD, Nicchitta C, Yan H: Glioblastoma proto-oncogene SEC61gamma is required for tumor cell survival and response to endoplasmic reticulum stress. Cancer Res. 2009, 69: 9105-9111. 10.1158/0008-5472.CAN-09-2775.PubMed CentralView ArticlePubMedGoogle Scholar
- Nord H, Hartmann C, Andersson R, Menzel U, Pfeifer S, Piotrowski A, Bogdan A, Kloc W, Sandgren J, Olofsson T, Hesselager G, Blomquist E, Komorowski J, von Deimling A, Bruder CE, Dumanski JP, Diaz de Stahl T: Characterization of novel and complex genomic aberrations in glioblastoma using a 32K BAC array. Neuro Oncol. 2009, 11: 803-818. 10.1215/15228517-2009-013.PubMed CentralView ArticlePubMedGoogle Scholar
- Yoshioka S, King ML, Ran S, Okuda H, MacLean JA, McAsey ME, Sugino N, Brard L, Watabe K, Hayashi K: WNT7A regulates tumor growth and progression in ovarian cancer through the WNT/beta-catenin pathway. Mol Cancer Res. 2012, 10: 469-482. 10.1158/1541-7786.MCR-11-0177.PubMed CentralView ArticlePubMedGoogle Scholar
- Shrestha Y, Schafer EJ, Boehm JS, Thomas SR, He F, Du J, Wang S, Barretina J, Weir BA, Zhao JJ, Polyak K, Golub TR, Beroukhim R, Hahn WC: PAK1 is a breast cancer oncogene that coordinately activates MAPK and MET signaling. Oncogene. 2012, 31: 3397-3408. 10.1038/onc.2011.515.PubMed CentralView ArticlePubMedGoogle Scholar
- Dranoff G: Cytokines in cancer pathogenesis and cancer therapy. Nat Rev Cancer. 2004, 4: 11-22. 10.1038/nrc1252.View ArticlePubMedGoogle Scholar
- McFarland BC, Ma JY, Langford CP, Gillespie GY, Yu H, Zheng Y, Nozell SE, Huszar D, Benveniste EN: Therapeutic potential of AZD1480 for the treatment of human glioblastoma. Mol Cancer Ther. 2011, 10: 2384-2393. 10.1158/1535-7163.MCT-11-0480.PubMed CentralView ArticlePubMedGoogle Scholar
- Takaoka A, Hayakawa S, Yanai H, Stoiber D, Negishi H, Kikuchi H, Sasaki S, Imai K, Shibue T, Honda K, Taniguchi T: Integration of interferon-alpha/beta signalling to p53 responses in tumour suppression and antiviral defence. Nature. 2003, 424: 516-523. 10.1038/nature01850.View ArticlePubMedGoogle Scholar
- Smith AL, Robin TP, Ford HL: Molecular pathways: targeting the TGF-beta pathway for cancer therapy. Clin Cancer Res. 2012, 18: 4514-4521. 10.1158/1078-0432.CCR-11-3224.View ArticlePubMedGoogle Scholar
- Basque J, Martel M, Leduc R, Cantin AM: Lysosomotropic drugs inhibit maturation of transforming growth factor-beta. Can J Physiol Pharmacol. 2008, 86: 606-612. 10.1139/Y08-063.View ArticlePubMedGoogle Scholar
- Swami M: Signalling: the calcium connection. Nat Rev Cancer. 2010, 10: 738-10.1038/nrc2956.View ArticlePubMedGoogle Scholar
- Desmet CJ, Gallenne T, Prieur A, Reyal F, Visser NL, Wittner BS, Smit MA, Geiger TR, Laoukili J, Iskit S, Rodenko B, Zwart W, Evers B, Horlings H, Ajouaou A, Zevenhoven J, van Vliet M, Ramaswamy S, Wessels LF, Peeper DS: Identification of a pharmacologically tractable Fra-1/ADORA2B axis promoting breast cancer metastasis. Proc Natl Acad Sci U S A. 2013, 110: 5139-5144. 10.1073/pnas.1222085110.PubMed CentralView ArticlePubMedGoogle Scholar
- Boldt S, Weidle UH, Kolch W: The role of MAPK pathways in the action of chemotherapeutic drugs. Carcinogenesis. 2002, 23: 1831-1838. 10.1093/carcin/23.11.1831.View ArticlePubMedGoogle Scholar
- Santarpia L, Lippman SM, El-Naggar AK: Targeting the MAPK-RAS-RAF signaling pathway in cancer therapy. Expert Opin Ther Targets. 2012, 16: 103-119. 10.1517/14728222.2011.645805.PubMed CentralView ArticlePubMedGoogle Scholar
- Sun H, Li XB, Meng Y, Fan L, Li M, Fang J: TRAF6 upregulates expression of HIF-1alpha and promotes tumor angiogenesis. Cancer Res. 2013, 73: 4950-4959. 10.1158/0008-5472.CAN-13-0370.View ArticlePubMedGoogle Scholar
- Reva B, Antipin Y, Sander C: Predicting the functional impact of protein mutations: application to cancer genomics. Nucleic Acids Res. 2011, 39: e118-10.1093/nar/gkr407.PubMed CentralView ArticlePubMedGoogle Scholar
- Ye H, Arron JR, Lamothe B, Cirilli M, Kobayashi T, Shevde NK, Segal D, Dzivenu OK, Vologodskaia M, Yim M, Du K, Singh S, Pike JW, Darnay BG, Choi Y, Wu H: Distinct molecular mechanism for initiating TRAF6 signalling. Nature. 2002, 418: 443-447. 10.1038/nature00888.View ArticlePubMedGoogle Scholar
- Annunziata CM, Stavnes HT, Kleinberg L, Berner A, Hernandez LF, Birrer MJ, Steinberg SM, Davidson B, Kohn EC: Nuclear factor kappaB transcription factors are coexpressed and convey a poor outcome in ovarian cancer. Cancer. 2010, 116: 3276-3284. 10.1002/cncr.25190.PubMed CentralView ArticlePubMedGoogle Scholar
- Iozzo RV, Sanderson RD: Proteoglycans in cancer biology, tumour microenvironment and angiogenesis. J Cell Mol Med. 2011, 15: 1013-1031. 10.1111/j.1582-4934.2010.01236.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang K, Chu K, Wu X, Gao H, Wang J, Yuan YC, Loera S, Ho K, Wang Y, Chow W, Un F, Chu P, Yen Y: Amplification of FRS2 and activation of FGFR/FRS2 signaling pathway in high-grade liposarcoma. Cancer Res. 2013, 73: 1298-1307. 10.1158/0008-5472.CAN-12-2086.View ArticlePubMedGoogle Scholar
- Wei X, Walia V, Lin JC, Teer JK, Prickett TD, Gartner J, Davis S, Program NCS, Stemke-Hale K, Davies MA, Gershenwald JE, Robinson W, Robinson S, Rosenberg SA, Samuels Y: Exome sequencing identifies GRIN2A as frequently mutated in melanoma. Nat Genet. 2011, 43: 442-446. 10.1038/ng.810.PubMed CentralView ArticlePubMedGoogle Scholar
- Guedj M, Marisa L, de Reynies A, Orsetti B, Schiappa R, Bibeau F, MacGrogan G, Lerebours F, Finetti P, Longy M, Bertheau P, Bertrand F, Bonnet F, Martin AL, Feugeas JP, Bieche I, Lehmann-Che J, Lidereau R, Birnbaum D, Bertucci F, de The H, Theillet C: A refined molecular taxonomy of breast cancer. Oncogene. 2012, 31: 1196-1206. 10.1038/onc.2011.301.PubMed CentralView ArticlePubMedGoogle Scholar
- Neve RM, Chin K, Fridlyand J, Yeh J, Baehner FL, Fevr T, Clark L, Bayani N, Coppe JP, Tong F, Speed T, Spellman PT, DeVries S, Lapuk A, Wang NJ, Kuo WL, Stilwell JL, Pinkel D, Albertson DG, Waldman FM, McCormick F, Dickson RB, Johnson MD, Lippman M, Ethier S, Gazdar A, Gray JW: A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer Cell. 2006, 10: 515-527. 10.1016/j.ccr.2006.10.008.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheang MC, Voduc D, Bajdik C, Leung S, McKinney S, Chia SK, Perou CM, Nielsen TO: Basal-like breast cancer defined by five biomarkers has superior prognostic value than triple-negative phenotype. Clin Cancer Res. 2008, 14: 1368-1376. 10.1158/1078-0432.CCR-07-1658.View ArticlePubMedGoogle Scholar
- Plaks V, Brenot A, Lawson DA, Linnemann JR, Van Kappel EC, Wong KC, de Sauvage F, Klein OD, Werb Z: Lgr5-expressing cells are sufficient and necessary for postnatal mammary gland organogenesis. Cell Rep. 2013, 3: 70-78. 10.1016/j.celrep.2012.12.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu X, Pan Z, Zhang L, Sun Q, Wan J, Tian C, Xing G, Yang J, Liu X, Jiang J, He F: JAB1 accelerates mitochondrial apoptosis by interaction with proapoptotic BclGs. Cell Signal. 2008, 20: 230-240. 10.1016/j.cellsig.2007.10.012.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.