<figure><img src="image_1.png" /></figure>

HAL
open science

# New insight for pharmacogenomics studies from the transcriptional analysis of two large-scale cancer cell line panels

Benjamin Sadacca, Anne-Sophie Hamy-Petit, Cecile Laurent, Pierre Gestraud, Hélène Bonsang-Kitzis, Alice Pinheiro, Judith Abecassis, Pierre Neuvial, Fabien Reyal

► To cite this version:

Benjamin Sadacca, Anne-Sophie Hamy-Petit, Cecile Laurent, Pierre Gestraud, Hélène Bonsang-Kitzis, et al.. New insight for pharmacogenomics studies from the transcriptional analysis of two large-scale cancer cell line panels. Scientific Reports, 2017, 7 (1), pp.15126. 10.1038/s41598-017-14770-6. hal-01634180

HAL Id: hal-01634180

https://hal.science/hal-01634180

Submitted on 12 Jan 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_3.png" /></figure>

OPEN

# New insight for pharmacogenomics studies from the transcriptional analysis of two large-scale cancer cell line panels

Received: 23 May 2017

Accepted: 12 October 2017

Published online: 09 November 2017

Benjamin Sadacca<sup>1,2,3</sup>, Anne-Sophie Hamy-Petit<sup>1,2</sup>, Cécile Laurent<sup>1,2</sup>, Pierre Gestraud<sup>5</sup>, Hélène Bonsang-Kitzis<sup>1,2,6</sup>, Alice Pinheiro<sup>1,2</sup>, Judith Abecassis<sup>1,2,4,5</sup>, Pierre Neuvial<sup>3,7</sup> & Fabien Reyal<sup>1,2,6</sup>

One of the most challenging problems in the development of new anticancer drugs is the very high attrition rate. The so-called "drug repositioning process" propose to find new therapeutic indications to already approved drugs. For this, new analytic methods are required to optimize the information present in large-scale pharmacogenomics datasets. We analyzed data from the Genomics of Drug Sensitivity in Cancer and Cancer Cell Line Encyclopedia studies. We focused on common cell lines (n = 471), considering the molecular information, and the drug sensitivity for common drugs screened (n = 15). We propose a novel classification based on transcriptomic profiles of cell lines, according to a biological network-driven gene selection process. Our robust molecular classification displays greater homogeneity of drug sensitivity than cancer cell line grouped based on tissue of origin. We then identified significant associations between cell line cluster and drug response robustly found between both datasets. We further demonstrate the relevance of our method using two additional external datasets and distinct sensitivity metrics. Some associations were still found robust, despite cell lines and drug responses' variations. This study defines a robust molecular classification of cancer cell lines that could be used to find new therapeutic indications to known compounds.

One of the most challenging problems in the development of new anticancer drugs is the very high attrition rate. Less than 5% of the drugs entering phase I trials eventually obtain marketing authorization<fnref n="1" />. Clinical trials are the only real way to assess drug efficacy and toxicity, but this approach is inadequate for testing the hundreds of drugs currently being developed². Scientists need to test hundreds of drugs on numerous tumor models therefore frequently make use of tumor-derived cell lines³⁻⁵. Such studies aim to identify genomic biomarkers for predicting the responses of individual patients to the drug and, ultimately, for identifying the best drug for each patient.

In 2012, the first large-scale pharmacogenomics studies provided an unprecedented wealth to the scientific community. The Broad Institute-Cancer Cell Line Encyclopedia (CCLE) provided a collection of 1,036 human cancer cell lines from 36 tumor types, tested for 24 anticancer drugs. The Genomics of Drug Sensitivity in Cancer (GDSC) assessed the sensitivity of 727 cell lines, from 29 tissue types, to 138 drugs. Both datasets contain genome-wide gene expression and sequencing data for a subset of genes. These studies have provided unprecedented amounts of information about molecular profiles and drug sensitivity and have validated several known genetic biomarkers, such as the BRAF-V600E mutation sensitizing melanomas to vemurafenib⁶ or ERBB2 amplification/overexpression conferring sensitivity to lapatinib⁷.

<footnote n="1">¹Residual Tumor & Response to Treatment Laboratory (RT2Lab), PSL Research University, Translational Research Department, F-75248, Paris, France. ²U932 Immunity and Cancer; INSERM; Institut Curie, Paris, France. ³Laboratoire de Mathématiques et Modélisation d'Evry, Université d'Évry Val d'Essonne, UMR CNRS 8071, ENSIIE, USC INRA, Evry Val d'Essonne, France. ⁴Mines Paristech, PSL-Research University, CBIO-Centre for Computational Biology, Mines ParisTech, Fontainebleau, F-77300, France. ⁵Institut Curie, PSL Research University, Mines Paris Tech, Bioinformatics and Computational Systems Biology of Cancer, INSERM U900, F-75005, Paris, France. ⁶Department of Surgery, Institut Curie, Paris, F-75248, France. ⁷Institut de Mathématiques de Toulouse; UMR5219 Université de Toulouse; CNRS UPS IMT, F-31062, Toulouse Cedex 9, France. Correspondence and requests for materials should be addressed to F.R. (email: fabien.reyal@curie.fr)</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Flow chart of the analysis. We apply the same pipeline of analysis independently to CCLE and GDSC. (a) Biologically driven gene selection was performed to build robust clusters of genes. (b) Robust clusters of cell lines were then built using the selected genes. (c) Cell lines clusters have been associated to distinct drug response.</figcaption></figure>

Previous studies assessed drug sensitivity by pooling all the cell lines or by controlling for tissue source. However, with improvements in our knowledge about tumors, it has become clear that genomic, epigenomic, transcriptional, and proteomic analyses of a given cancer can reveal subtypes differing in pathway activity, progression or treatment response<sup>8,9</sup>. Conversely, the recent success of basket studies<sup>10,11</sup> have demonstrated that treatment choices can be based on abnormalities shared by tumors originating from different tissue types.

We present here a comprehensive reanalysis of these two recently published large-scale pharmacogenomics resources. We propose an alternative approach in which cell lines are grouped by transcriptomic profile, based on a biological network-driven gene selection process. This molecular classification of cancer cell lines appeared robust across CCLE and GDSC. We further demonstrated the relevance of this novel classification through the drug response We validate our approach by robustly found in CCLE and GDSC as in two external dataset the significant associations between cell line clusters and drug responses.

## Results

**A biologically driven approach identifies four robust gene modules.** Gene expression profiles were recovered for 471 cell lines, from 24 different tissues, tested in both CCLE and GDSC. Data were curated and annotated with the pipeline of Haibe-Kains *et al.*<sup>12</sup>. We developed a three-step biological network-driven process based on transcriptomic data for identifying robust clusters of genes. This process was applied in parallel for each dataset. We first selected the most variant genes from the set of 12,153 genes common to GDSC and CCLE, by the inflexion point method. We then performed hierarchical consensus clustering<sup>13</sup> to identify robust gene modules. Finally, we used String® database software<sup>14</sup> to analyze our gene selection. The goal was to decrease the heterogeneity of each gene cluster. We retained the genes from our initial selection that had (1) high String® database gene connection indices (greater than 0.7), and (2) similar patterns of expression to other genes within the same biological network (correlation coefficient of at least 0.5) (Fig. 1 step A). This selection process identified four stable clusters in GDSC (n = 183 genes) and five in CCLE (n = 210 genes), including a subset of 170 genes common to the two datasets. Distinct functional gene ontologies were associated to each gene modules based on a gene ontology analysis: (Supplementary Fig. 1) Gene Cluster - Extracellular Matrix (GC-ECM; n_ccle = 48, n_GDSC = 36), Gene Cluster - Migration (GC-Migration; n_ccle = 56, n_GDSC = 75), Gene Cluster - Immunity-Interferon (GC-Immunity; n_ccle = 22, n_GDSC = 14) and Gene Cluster - Epithelial Phenotype (GC-Epithelial; n_ccle = 63, n_GDSC = 58). A set of 21 genes enriched in development processes (GC-Development) was found exclusively in the CCLE dataset.

**Biologically driven gene selection identifies eleven reproducible cell line clusters.** We performed a consensus clustering with the previously selected genes, for each dataset separately, to identify global
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Cell line clustering with CCLE data. (a) Heatmap clustering with 471 cell lines (in columns) and 210 selected genes (in rows) for the CCLE data. (b) EMT status of the cell lines.</figcaption></figure>

differences in gene expression between cancer cell lines (Fig. 1 step B). We obtained eleven stable clusters of cell lines in CCLE and GDSC (Fig. 2a and Supplementary Fig. 2a).

Previous studies reported strong correlations between the expression profiles of identical cell lines¹². We therefore investigated the closeness of the cell line clusters obtained. We defined the similarity between any two cell lines as the number of datasets in which they clustered together (0 = none, 1 = CCLE or GDSC, 2 = CCLE and GDSC). We assessed the consistency between the clustering patterns obtained with CCLE and GDSC data, using a heatmap clustering of the similarity matrix as a visualization tool. The heatmap shows the number of times that two samples are clustered together across datasets (Fig. 3a). Groups of cell lines that frequently cluster with each other are shown in darker shades of blue. The heatmap revealed a well defined 11-block, corresponding to the 11 clusters previously identified. A high degree of consistency between the 11 clusters was observed, with 90% accuracy. As the cell line clusters were highly similar, we use the term “cluster” to denote the same group of cell lines from CCLE and GDSC, unless the dataset is specified.

**Tissue-of-origin or transcriptomic features dominate cell line clusters.** Our eleven clusters can be organized in three major patterns: (i) four clusters of cell lines were derived mostly from tumors from the same tissue of origin. These cell line clusters were named after the organ or cancer subtypes from which most of the cell lines were derived: hematopoietic and lymphoid tissues (HAL), small cell lung cancer (SCLC), skin (SKCM) and breast (BRCA) clusters; (ii) four clusters of cell lines were derived from tissues from the same organ system or had a common embryonic origin: gastrointestinal tract (GI), aerodigestive tract (ADG), glioma and sarcoma (GLSR) and endodermal origin tumors (EDOT) clusters; (iii) three clusters contained cell lines from different tissues of origin. These clusters were named Mixed 1, Mixed 2 and Mixed 3 (Fig. 3b and c). Details provided in Supplementary data 1 and 2).

**Clusters of cell lines with common presumptive tissues of origin.** Four cell line clusters appeared very homogeneous in terms of tissue lineage: HAL, SCLC, SKCM and BRCA. These lineages accounted for 84%, on average, of the cells of their respective clusters. The HAL cluster grouped together all the cell lines originating from hematopoietic and lymphoid tissues. This clear clustering pattern can be accounted for by the hematopoietic phenotype of this type of tumor. The SKCM cluster was the second most homogeneous cell line cluster in terms of tissue type (92% of the cell lines in this group originated from melanomas). Breast cancer is a heterogeneous disease with a growing number of recognized biological subtypes, including ER+ Her2−, Her2+ and triple-negative breast cancer (TNBC), which is the most aggressive subtype. BRCA cluster contained all the breast cancer cell lines defined as ER+ Her2− (7/7) and Her2+ (7/7). However, only about half the cell lines defined as triple-negative
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Clustering similarity. (a) Color-coded heatmap for similarity between CCLE and GDSC clustering; Tag Cloud represents the tissue composition of cell lines cluster, in CCLE (b) and GDSC (c). The importance of each tissue is indicated by font size. The TNBC cell lines belonging to each cluster are indicated by red dots.</figcaption></figure>

belonged to this cluster (11/20 in GDSC, 8/20 in CCLE). The remaining triple-negative breast cancer cell lines were found in six different clusters of cell lines (SCLC, EDOT, Mixed 1, Mixed 2, GLSR and ADG) (Fig. 3b,c). SCLC cluster contained 28% of the lung cancer cell lines and 45% of the small-cell lung carcinoma cell lines. We performed a Gene Set Enrichment Analysis¹⁵ (GSEA) based on our previously defined gene modules to characterize the transcriptomic profile of cell line clusters (Supplementary Fig. 3). The immunity gene module was strongly expressed in the cell lines of the HAL cluster. Leukemia affects both the bone marrow and lymphocytes, potentially accounting for the detection of immunity gene expression in cell lines derived from a tumor system with no stromal environment. In the SKCM cell line cluster, the epithelial phenotype gene module was down-regulated. Furthermore, the activation of the ECM and migration gene modules in this cluster is suggestive of aggressive cancer. In the BRCA and SCLC cell line clusters, the epithelial gene module was expressed, whereas the migration and ECM gene modules were not.

*Clusters of cell lines from tissues of the same organ system or common embryonic origin.* Some clusters could not be defined on the basis of origin from a single tissue type. However, with a more systemic vision, a consistent organization was obtained for four clusters: GI, ADG, GLSR and EDOT. Cell lines derived from tumors of the digestive system belonged to two clusters. The ADG cell line cluster consisted mostly of tumors from the esophagus, upper aerodigestive tract, salivary and also urinary glands, whereas the GI cluster grouped together tumors derived from large intestine, stomach and pancreas cancers. About 70% of the cell lines of the GLSR cluster were derived from tumors of the central nervous system, bone, autonomic ganglia and soft tissue. Finally, the EDOT cell line cluster grouped together cell lines derived from tumors of different tissues (e.g. lung, pancreas, urinary tract) arising from the same germ layer (endoderm). The relevance of the EDOT cluster is supported by studies suggesting that oncogenesis may be initiated by the activation of a common pathway in an endodermal progenitor¹⁶.

The ADG, GI and EDOT clusters all displayed strong expression of the genes of the epithelial phenotype module and weak expression of the ECM gene module. According to GSEA, the migration gene module was less strongly expressed in GI cells. For the EDOT cluster, inconsistencies between the CCLE and GDSC datasets were observed concerning the activation or inhibition of migration gene expression at the transcriptomic level only. The GLSR cluster displayed low levels of expression for the epithelial gene module, and high levels of expression for the ECM and migration modules.

*Clusters of cell lines from tumors with heterogeneous tissues of origin.* Three clusters displayed no particular prevalence of cell lines corresponding to any particular tissue or organ system. They contained cell lines from tumors
<table><caption>**Table 1.** Significant associations found between CCLE, GDSC, GSK and GCSI. In bold associations found significant in at least three datasets. The association between PD0325901 and SKCM had an adjusted p-values of 0.058 (marked with*).</caption><thead><tr><th colspan="3">CCLE vs GDSC</th><th colspan="3">CCLE vs GDSC</th><th colspan="3">CCLE vs GSK</th></tr><tr><th colspan="3">IC50</th><th colspan="3">AUC</th><th colspan="3">IC50</th></tr><tr><th>Drug</th><th>Cluster</th><th>Response</th><th>Drug</th><th>Cluster</th><th>Response</th><th>Drug</th><th>Cluster</th><th>Response</th></tr></thead><tbody><tr><td>Erlotinib</td><td>ADG</td><td>Sensitive</td><td>Erlotinib</td><td>ADG</td><td>Sensitive</td><td>Lapatinib</td><td>Mixed 1</td><td>Resistant</td></tr><tr><td>AZD6244</td><td>SKCM</td><td>Sensitive</td><td>AZD6244</td><td>SKCM</td><td>Sensitive</td><td>Lapatinib</td><td>SKCM</td><td>Resistant</td></tr><tr><td>AZD6244</td><td>BRCA</td><td>Resistant</td><td>AZD6244</td><td>BRCA</td><td>Resistant</td><td></td><td></td><td></td></tr><tr><td>Lapatinib</td><td>SCLC</td><td>Resistant</td><td>Lapatinib</td><td>HAL</td><td>Resistant</td><td></td><td></td><td></td></tr><tr><td>Lapatinib</td><td>ADG</td><td>Sensitive</td><td>Crizotinib</td><td>SKCM</td><td>Resistant</td><td></td><td></td><td></td></tr><tr><td>PD0332991</td><td>GI</td><td>Resistant</td><td>AZD0530</td><td>SKCM</td><td>Resistant</td><td></td><td></td><td></td></tr><tr><td>PD0332991</td><td>HAL</td><td>Sensitive</td><td>PLX4720</td><td>SKCM</td><td>Sensitive</td><td></td><td></td><td></td></tr><tr><td>PLX4720</td><td>SKCM</td><td>Sensitive</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PD0325901</td><td>GI</td><td>Sensitive</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PD0325901</td><td>SKCM</td><td>Sensitive</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th colspan="3">CCLE vs gCSI</th><th colspan="3">GDSC vs gCSI</th><th colspan="3"></th></tr><tr><th colspan="3">IC50</th><th colspan="3">IC50</th><th colspan="3"></th></tr><tr><th>Drug</th><th>Cluster</th><th>Response</th><th>Drug</th><th>Cluster</th><th>Response</th><th>Drug</th><th>Cluster</th><th>Response</th></tr><tr><td>Erlotinib</td><td>ADG</td><td>Sensitive</td><td>PD0325901</td><td>SKCM</td><td>Sensitive*</td><td></td><td></td><td></td></tr><tr><td>Erlotinib</td><td>Mixed 1</td><td>Resistant</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Erlotinib</td><td>GLSR</td><td>Resistant</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Erlotinib</td><td>SKCM</td><td>Resistant</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Lapatinib</td><td>Mixed 1</td><td>Resistant</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Lapatinib</td><td>ADG</td><td>Sensitive</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PD0325901</td><td>BRCA</td><td>Resistant</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PD0325901</td><td>SKCM</td><td>Sensitive</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th colspan="3">CCLE vs gCSI</th><th colspan="3">GDSC vs gCSI</th><th colspan="3"></th></tr><tr><th colspan="3">Mean Viability</th><th colspan="3">Mean Viability</th><th colspan="3"></th></tr><tr><th>Drug</th><th>Cluster</th><th>Response</th><th>Drug</th><th>Cluster</th><th>Response</th><th>Drug</th><th>Cluster</th><th>Response</th></tr><tr><td>Erlotinib</td><td>ADG</td><td>Sensitive</td><td>PD0325901</td><td>SKCM</td><td>Sensitive</td><td></td><td></td><td></td></tr><tr><td>Erlotinib</td><td>Mixed 1</td><td>Resistant</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Erlotinib</td><td>GLSR</td><td>Resistant</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Erlotinib</td><td>SKCM</td><td>Resistant</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Erlotinib</td><td>HAL</td><td>Resistant</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Erlotinib</td><td>SCLC</td><td>Resistant</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PD0325901</td><td>BRCA</td><td>Resistant</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PD0325901</td><td>SKCM</td><td>Sensitive</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

of 11 to 16 different tissues. We named these clusters Mixed 1, Mixed 2 and Mixed 3. All three of these clusters displayed low levels of epithelial phenotype genes, suggesting that the cell lines they contained were probably mesenchymal. These clusters also displayed an upregulation of ECM genes. Mixed 1 and 2 displayed an upregulation of migration gene expression. These results suggest that some of the cell lines may have been metastatic in origin or subject to drift, from the characteristics of the tissue of origin to a less differentiated state. In this case, transcriptomic profile is more relevant than tissue of origin.

**EMT discriminates between cell line clusters.** The identification of an epithelial phenotype gene module led us to investigate the epithelial-mesenchymal status of each cell line. A previous study¹⁷ showed that epithelial/mesenchymal transition (EMT)-associated differences in gene expression were a major determinant of the stratification of cancer cell lines based on transcriptomic profiles. Indeed, we found a significant overlap between our gene selections and a published EMT-derived gene signature consisting of 249 genes¹⁸ (*P* < 0.0001, two-tailed Fisher's exact test). We superimposed epithelial/mesenchymal cell line classifications over our gene expression clusters and found a strong association (Fig. 2b and Supplementary Fig. 3b). According to the EMT signature, five cell line clusters (SCLC, GI, EDOT, ADG and BRCA) contained mostly epithelial cell lines, whereas the Mixed 1, Mixed 3, GLSR and SKCM cell line clusters contained mostly mesenchymal cell lines. The Mixed 2 cell line cluster appeared to contain mostly mesenchymal cell lines in GDSC but almost half the cell lines assigned to this cluster in CCLE were epithelial. The HAL cell lines were not concerned by this stratification. Finally, the epithelial/mesenchymal classification was consistent with that obtained with the epithelial phenotype gene module.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Pseudo *F* value for the 15 drugs common to CCLE and GDSC. The pseudo *F* index have been computed from the IC₅₀ values for each drug. The pseudo *F* statistic is the ratio of between-cluster variance to within-cluster variance. Large values of pseudo *F* indicate well-separated, tight clusters. Drugs are listed in descending order of pseudo *F* values for clustering.</figcaption></figure>

**Cell line clusters are enriched in somatic mutations.** We investigated a common set of 64 genes for the presence of mutations in CCLE and GDSC datasets. However, many inconsistencies between both datasets led us to focus on a set of eight genes (TP53, KRAS, NRAS, APC, PIK3CA, BRAF, PTEN and RB1) for which at least 5% of identical cell lines display mutations in both datasets (Supplementary Information). The mutational profile of cell line clusters was then described based on these genes. Mutation profiles clearly distinguished four clusters (Fig. 2a). The SCLC cluster was enriched in RB1 mutations. The GI cluster was rich in APC and KRAS mutations; NRAS mutations were overrepresented in the HAL cluster and the SKCM cluster was enriched in BRAF mutations. Finally, KRAS mutations were particularly abundant in the EDOT clusters. No significant enrichment in mutations was observed for the GLSR, ADG, BRCA and Mixed 3 cell line clusters (Supplementary Tables 1 and 2). These clusters have fewer mean mutation rates than the other clusters (GDSC: 13% vs. 19%, *t*-test *p*-value = 0.01; CCLE: 17% vs. 22%, *t*-test *p*-value = 0.08).

**Transcriptomic clustering is more consistent than clustering on the basis of tissue of origin in terms of drug responses.** The large-scale drug screening programs of the Broad and Sanger Institutes have provided to the scientific community an unprecedented wealth of publicly available data. Molecular data have been systematically collected for each cell line, but far less information is available for drug screening (Supplementary Information). Moreover, in many cases (25% in CCLE and 45% in GDSC) it was not possible to extract the IC₅₀ from the dose-response curve. In order to overcome these issues, both study also report the AUC (area under the dose response curve) that can always be calculated.

We evaluated whether our clustering was more discriminant than the tissue of origin of the cell lines, in terms of drug response. We calculated a pseudo *F*-statistic separately for IC₅₀ and AUC values for each of the 15 drugs common to CCLE and GDSC. This measurement should capture consistency between the clustering and screening data. It is calculated as the ratio of between-group variance in drug response to the corresponding within-group variance¹⁹. High pseudo *F* values indicate well-separated, compact clusters. We then compared the pseudo *F* values calculated with our clustering method with those obtained for 'tissue partitioning' for a given drug (i.e. each tissue being to correspond to a cluster of cell lines).

Twelve of the fifteen drugs had a higher ratio in CCLE and GDSC for our clustering than for clustering based on tissue of origin with the IC50 (Fig. 4) and ten out of fifteen with the AUC (Supplementary Fig. 4). This trend was confirmed by a *t*-test comparing the pseudo *F* values for our clustering with those for 'tissue partitioning' (IC50: CCLE t.test *p*-value = 0.041, GDSC t.test *p*-value = 0.032, AUC: CCLE t.test *p*-value = 0.011, GDSC t.test *p*-value = 0.043). PLX4720 (Raf kinase B inhibitor) and PD0325901 (MEK1 and MEK2 inhibitors) were drugs with the largest pseudo *F* values in both dataset. Paclitaxel was the only molecule in the panel with a higher pseudo *F* value for tissue partitioning in CCLE and GDSC. As the drug sensitivity results were not used to determine the clustering of the cell lines, these findings provide independent evidence for a major role of mRNA levels in drug sensitivity.

**Robust identification of drug response across datasets.** Subgroups of patients or cell lines defined on the basis of transcriptomic data have been shown to be associated with differences in drug sensitivity⁸,⁹. We sought to identify associations between clusters of cell lines and “sensitive” or “resistant” drug phenotypes, for the
15 drugs tested in both CCLE and GDSC. For each dataset and each drug separately, we investigated whether the mean IC₅₀ of a given cell line cluster differed significantly from those for the other cell line clusters (see Fig. 1 step C and Materials and Methods). Six molecules were found to be significantly associated with six different clusters in both CCLE and GDSC (Table 1 and Supplementary Table 3, Supplementary Fig. 5). The SKCM and GI cell line clusters were both significantly more sensitive than the other cell lines to PD0325901 (MEK 1 and MEK 2 inhibitors) (Fig. 5a). The association of melanoma and PLX4720 (Raf kinase B inhibitor) is already well established and was confirmed by our analysis. Moreover, an inhibitor of MEK 1 and MEK 2, AZD6244, displayed significantly higher levels of activity in cell lines from the SKCM cell line cluster. Both EGFR inhibitors, erlotinib (Fig. 5b) and lapatinib, appeared to be significantly more effective against ADG cell lines than against other cell lines. Hematopoietic and lymphoid tissue cells were sensitive to the CDK4/6 inhibitor PD033991. By contrast, SLCL cell lines appeared to be resistant to lapatinib (EGFR and HER2 inhibitor) and the CDK4/6 inhibitor PD033991 was found inefficient to kill GI cell lines. Finally, AZD6244 (inhibitor of MEK1 and MEK2) appeared ineffective to treat BRCA cell. In addition to variation between drug sensitivity and cell lines, previous studies report variations across the different metrics used to report the drug efficacy¹,²,²⁰. We then performed similar analysis using AUC. More than half of the associations between cell lines clusters and drug sensitivity were found still significant with AUC (Table 1 and Supplementary Information).

We further evaluated the relevance of our clustering regarding the drug sensitivity using two external public datasets. We first study the 118 cell lines tested in common between the CCLE and the GlaxoSmithKline cell line collection (GSK)²¹ on lapatinib and paclitaxel (GDSC was excluded due to small sample size, see Supplementary Information). We found that lapatinib was significantly inactive to kill cells from clusters SKCM and Mixed 1 in both CCLE and GSK (Table 1). Since the set of common cell lines and drugs was small between CCLE, GDSC and GSK (Supplementary Table 4), we consider the Genentech Cell Line Screening Initiative (gCSI)²². A panel of 244 unique cell lines and 5 drugs overlap between CCLE, GDSC and gCSI. Instead of AUC, the gCSI reported the mean viability statistic to measure drug efficacy in addition to the IC₅₀. Eight associations between cell lines clusters and drug sensitivity were found significant using the IC₅₀ and nine with the mean viability statistic. Among them, the sensitivity of ADG to erlotinib and lapatinib as well as the efficacy of PD0325901 to kill cells from SKCM cluster were common to CCLE, GDSC and gCSI (Table 1 and Supplementary Information). Our results suggest that our cell line clustering is able to find significant associations with drugs efficacy robustly in four different dataset, despite the large variations across pharmacological data and drug response measures.

**Distinct drug profiles were associated with the various cell line clusters.** We applied the same procedure to all the drugs tested in the CCLE (24 molecules) and GDSC (129 molecules) studies. For each dataset and each of the 153 drugs separately, we determined whether the mean IC₅₀ of a given cell line cluster was significantly different those of the other cell line clusters (Supplementary Tables 5 and 6). Overall, the most striking result was the very small number of drugs associated with a sensitive profile (88 associations, including 71 unique drugs) compared to drugs associated with a resistant profile (163 associations, including 92 unique drugs) (Supplementary Information). It was particularly interesting to observed that Mixed 2 and Mixed 3 clusters were each sensitive to only one drug: respectively midostaurin and vorinostat. Both drugs are targeted agents (PI3K/mTOR inhibitor and HDAC inhibitor). These clusters are made of several cells from different tissue of origin. However, we were able to identify targeted therapies active to kill those cells. These results provide further evidences that our clustering can identify relevant groups of cell sharing unknown features associated to targeted drugs.

Overall, these results suggest that cancer cell lines can be classified, on the basis of their transcriptomic profile, into 11 clusters that may or may not be specific to the tissue of origin. We demonstrated that transcriptomic clustering was more consistent than clustering on the basis of tissue of origin in terms of drug response whatever the drug sensitivity metric considered. We were also able to find several significant associations between clusters of cell lines and “sensitive” or “resistant” drug phenotypes. Many of these associations were robustly found across four different datasets with three different drug response metrics. As the drug sensitivity results were not used to determine the clustering of the cell lines, these findings provide independent evidence about the relevance of this new classification. Furthermore, we show that when we are trying to associate a group of genes from a consistent biological pathway with a group of cell lines, rather than a single gene with a single drug, robust associations can be established across several pharmacologic datasets.

## Discussion

Despite the progress in the development of *in vivo* models, cancer cell lines remain a key tool in cancer research. Patients are usually treated with combination therapy. However, it is important to better understand the mechanisms involved with monotherapies before moving forward to study combination therapies. Here, we introduce a new cell line classification constructed from 471 cell lines derived from tumors from 24 different tissues. A biological network analysis for the most variant genes identified 11 clusters of cell lines. These clusters appeared robust in two large-scale cell line panels. This biologically driven gene selection process, which is probably less sensitive to sample fluctuations than other methods, made it possible to capture strong biological signals that might be concealed by the noise present in microarray data. Several studies have reported that the incorporation of network information improves the stability of gene selection and the biological interpretability of biomarker signatures for a given prediction accuracy²³⁻²⁵.

In this new classification, a clear distinction was established between non-epithelial cancer cell lines (GLSR, SKCM, Mixed 3) and epithelial cell lines (EDOT, BRCA, GI). This suggests that EMT-associated differences in gene expression are major determinants of the gene expression-based stratification of cancer cell lines. This new molecular clustering system classified more than 65% of the cell lines differently from the currently used
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Distribution of IC₅₀ values for each in CCLE and GDSC. Ordered according to mean IC₅₀ for the cluster. From resistant (left) to sensitive (right).</figcaption></figure>

tissue-of-origin cell line classification system. Only four clusters consisted mostly of cell lines originating from a single tissue. Furthermore, three clusters include cells with expression profiles stronger than that of the original tissue (Mixed clusters). Thus, 25% of the cell lines displayed no link to any tissue of origin or related organ system.

One of the most interesting cases was the triple-negative breast cancer (TNBC). We focused on this subtype, as it is the only subtype of breast cancer without any targeted therapy associated. TNBC were found to be highly heterogeneous, falling into six different clusters. This divergence shows the relevance of studying cell lines from various tumor types. Drug response was dependent on cluster membership, with the EDOT cluster sensitive to chemotherapy, whereas the BRCA cluster was resistant. The widely dispersed TNBC cell lines were mostly mesenchymal, whereas the cell lines of the BRCA cluster were exclusively epithelial. TNBC is increasingly emerging as a heterogeneous disease²⁶,²⁷, with tumors differing in histological features, gene expression profiles, clinical behavior, overall prognosis²⁸ and sensitivity to systemic treatment⁹,²⁹,³⁰. These findings provide strong evidence to suggest that TNBC heterogeneity is reflected at the cell line level. Our results suggest also that particular attention should be paid to the selection of cell lines for studies of particular types or subtypes of cancer.

By analyzing several large-scale public data sets, we demonstrated that drug efficacy is significantly associated to transcriptomic profile. A comparative analysis recently showed that the gene-expression profiles of the 471 cell lines shared by CCLE and GDSC were highly concordant whereas the reported cell-line drug sensitivities for the 15 drugs tested in both studies were highly inconsistent¹². The authors put forward several hypotheses to explain these discrepancies, including differences in experimental protocols, the viability assay and procedures for summarizing dose response and non-observed IC₅₀ (the half maximal inhibitory concentration). Despite discrepancies between the drug sensitivity data retrieved from different databases, we were able to find some robust combinations. Well-known drug associations were found, such as the sensitivity of SKCM lines to vemurafenib. We also found that cancers with BRAF mutations, such as melanoma³¹ and cancers with KRAS/BRAF mutations, such as colorectal cancer³², were more sensitive to MEK inhibitors. Furthermore, CDK4/6 inhibition-induced cell death has been noted in cell lines and xenografts derived from patients with T-cell leukemia³³. SCLC cell lines have been shown to be resistant to lapatinib, but combination with a cytotoxic agent may yield promising results³⁴.

The decline in the number of new treatments approved in recent years is a major challenge for the pharmaceutical industry. One of the reasons for this decline is the lack of systematic evaluation of therapeutic indications for a drug that is either in advanced development phase or has already obtained a marketing authorization. The so-called "drug repositioning process" proposed to find new therapeutic indications to already approved drugs with faster development times and reduced risks. Furthermore, it allows patients to have access to earlier therapeutic advances³⁵. Several robust associations were found. Targeted drugs were found efficient to treat clusters of cell lines constituted of cell from different tissues. These drugs are known to be active in one or several tissues that constitute these clusters. It would be of particular interest to test specifically these drugs on the other tissues
represented in these cell lines clusters. For example, cluster ADG is mostly constituted of upper-aerodigestive, oesophagus and urinary tract cancer cell lines. ADG cluster was particularly sensitivity to the anti EGFR - erlotinib. If EGFR is a validated target for upper-aerodigestive cancer³⁶⁻³⁸ the therapeutic potential of erlotinib has already been highlighted for bladder cancer³⁹ and showed promising results in phase II for oesophagus cancer⁴⁰,⁴¹.

Different types of drugs have been used in the panels. Around 10% of the 153 drugs screened in CCLE and GDSC, and only 1 out of the 15 drugs in common to both studies, are cytotoxic agents. These drugs are expected to be broadly active among the cell line panel since they are not specific molecules. On the contrary, targeted agents are expected to be active only in a subset of cell lines, at least, those carrying the given target. Furthermore, the recent study published by Rees et al.⁴² demonstrated that target's expression and drug sensitivity were correlated in only 31% of the cases. Grouping cell lines on the basis of their transcriptomic profiles makes it possible to identify subsets of cells with common off-target features. It is then more relevant to compare the drug sensitivity between cell lines of these groups rather than examined the correlation of response of each cell line to a particular drug reported by one dataset with the response of the same cell line to the same drug reported by another dataset. These results suggest that when robust clusters of cell lines based on biologically network-driven approach are considered, consistency between drug responses can be achieved.

In conclusion, our cell line classification provides novel insight for pharmacogenomics studies. As cell lines remain the most widely used models for the preclinical evaluation of candidate cancer drugs, further investigation should be made to use this classification in the development of cancer treatments with the aim of reducing the attrition rate.

## Materials and Methods

**Pharmacogenomics data.** We collected data from the Broad and Sanger Institutes. The CCLE profiled 24 anticancer drugs on 1,036 cell lines. The GDSC screened 138 drugs on 727 cell lines. Both datasets contain genome-wide gene expression and massive parallel sequencing data. All data were recovered, curated and annotated with the pipeline developed by Haibe-Kains et al.¹² (the GDSC was referred to the Cancer Genome Project [CGP] in Haibe-Kains et al.). We used this pipeline as described in the original article, but with a different method for the normalization of gene expression. Haibe-Kains et al. normalized gene expression data by frozen robust multiarray analysis, fRMA⁴³. This method was designed to combine several datasets and overcome multiple batch issues. This strategy is relevant when trying to ensure assay reproducibility. Even though this approach would be unlikely to have a major effect on gene expression values, we chose to normalize the gene expression data separately with RMA⁴⁴, to ensure that the two datasets were perfectly independent. Our analysis focused on 471 cell lines and 15 drugs for which we have transcriptomic and drug sensitivity data available in both the CCLE and GDSC studies. We also collected two large datasets to validate our classification. Data from the GlaxoSmithKline cell line collection were retrieved from Haibe-Kains et al.¹². The Genentech Cell Line Screening Initiative data were available from compareDrug Screens R package published by Haverty et al.²².

**Gene expression data.** Transcriptomic data were restricted to the 12,153 genes common to the two technologies used by GDSC and CCLE (Affymetrix GeneChip HG-U133A and HG-U133PLUS2, respectively). The Jetset method⁴⁵ was used to select a unique probe set for given genes. The same probe set was used in both datasets for 83% of the genes.

**Drug sensitivity data.** The micromolar concentration (μM) at which the drug inhibited 50% of maximal cell growth was used to assess drug sensitivity as well as the area under the dose response curve (AUC). We also consider the mean viability statistic when comparing with gCSI. These measurements were converted to a common scale (−log10(M) for IC50, [0,1] for AUC and 1 – mean viability for mean viability), such that high values would be correspond to cell lines sensitive to drugs.

**Gene selection by the inflexion point method.** We selected the most variant genes, based on the inflexion point of the interquartile range (IQR) distribution for gene expression. This method is more data-driven than a fixed threshold for defining the proportion of genes displaying the highest level of variation. The full procedure is described below. For each gene, we: (1) calculated the IQR for all cell lines, (2) sorted the IQR values of the genes in ascending order, to generate an ordered distribution, (3) estimated the major inflection point of the IQR curve as the point on the curve furthest away from a line drawn between the start and end points of the distribution, and (4) retained genes with an IQR higher than the inflection point.

**Gene expression-based identification of cell line clusters.** We developed a biological network-driven process based on transcriptomic data, to identify robust clusters of genes and cell lines. This process can be broken down into two parts: (A) identification of robust clusters of genes, used for (B) identification of robust clusters of cell lines.

(A) The gene selection process is a three-step procedure. (1) We selected the most variant genes from among the 12,153 genes common to GDSC and CCLE, by the inflexion point method. (2) We performed hierarchical consensus clustering (ConsensusClusterPlus R Package) to identify robust gene modules. The consensus-clustering step, based on Pearson distance and Ward linkage, identified robust clusters of genes. It involved hierarchical clustering by resampling (1,000 iterations) randomly selected genes. (3) We identified known biological networks, for each gene cluster separately, using String© database software version 9.1 (http://string-db.org/). We then applied a two-step selection process: (1) we selected strong biological networks by retaining only genes for which connection scores of at least 0.7 were obtained with String© database software, (2) within each biological network, we selected groups of genes for which expression levels were correlated, with a correlation coefficient of at least 0.5. We used the R package clusterProfiler⁴⁶ for comparing and visualizing gene ontologies profiles among
gene modules. (B) We applied a consensus-clustering with hierarchical clustering to the cell line gene expression profiles, using the selected genes to visualize the optimal number of stable cell line clusters.

**Characterization of cell line clusters at the transcriptomic and mutational levels.** Gene set enrichment analysis (GSEA) was performed on genes modules built in step (A) of the biological network-driven process described above. We identified up-regulated or down-regulated gene modules, associated with each cell line cluster. An analysis was first performed to identify genes differentially expressed between a particular cluster and all the other cell lines, based on a linear model. For a given cluster k, cell lines were partitioned into two groups $j = \{Cluster-k, non-Cluster-k\}$. We then performed a differential analysis by comparing the mean gene expression of each group in a linear model (limma R package⁴⁷). The analysis was performed separately for each dataset. The results were used to rank genes in order of significance and to search for overrepresented gene modules, by pre-ranked gene set enrichment analysis (GSEA).

Genes with significantly higher frequencies of mutation in a given cluster were identified by one-tailed Fisher's exact tests. We compared the occurrence of any given mutation in each cell line clusters with that in all the remaining clusters combined.

**Identification of cell line clusters common to different studies.** We studied the likeness between the clusterings for CCLE and GDSC, by clustering the cell lines with a similarity matrix (hierarchical clustering with Pearson's metric and the Ward agglomerative method). The similarity matrix contains the number of times two cell lines are clustered together in each dataset ($0 = never, 1 = only$ in one classification, $2 =$ in both classifications). This similarity matrix constitutes a natural visualization tool for assessing the consistency between two clustering patterns. In particular, if we associate a color gradient to the 0–2 range of real numbers, such that white corresponds to 0, and dark blue corresponds to 2, and if we assume that the matrix is arranged so that items belonging to the same cluster are adjacent to each other (with the same item order used to index both the rows and the columns of the matrix), a matrix corresponding to a perfect consensus will be displayed as a color-coded heatmap characterized by blue blocks along the diagonal, on a white background. The accuracy was calculated as the number of times two cell lines clustered together divided by the number of possible combinations

**EMT cell line classification.** The “epithelial” or “mesenchymal” status of each cell line was defined with the signature identified by Taube *et al.*¹⁸. This epithelial-to-mesenchymal transition signature consists of 159 down-regulated genes and 90 upregulated genes. We performed a hierarchical clustering of cell lines based on these 249 genes and labeled clusters of cell lines according to the overexpression of known epithelial marker genes, known mesenchymal marker genes or neither.

**Definition of breast cancer subtypes.** Breast cancer subtypes were defined with a bimodal mixture of two Gaussian distributions for ESR1, PGR and ERBB2 gene expression. Triple-negative (TN) breast cancer cell lines were defined by an absence of estrogen and progesterone receptor expression and a lack of ERBB2 overexpression/amplification ($n = 31$). We subsequently defined breast cancer cell lines overexpressing ESR1 but with a lower level of ERBB2 expression as the ER+Her2- subtype ($n = 7$), with cell lines overexpressing the ERBB2 gene defined as the Her2+ subtype ($n = 7$).

**Impact of cell line clustering on drug sensitivity.** We investigated the relevance of our clustering for drug sensitivity, by comparing the results obtained for this method with those for ‘tissue partitioning’ (i.e. each tissue of origin being considered to correspond to a cluster of cell lines). We calculated the pseudo F index computed from any drug sensitivity statistic ($IC_{50}$, AUC, mean viability) for each drug. The pseudo F statistic is the ratio of between-cluster variance to within-cluster variance¹⁹. It is defined as $[Between-cluster variance/(N-K)]/[Within-cluster variance/(K-1)]$, where N is the number of observations ($N = 471$) and K is the number of clusters ($K = 11$ or $K = 24$). Large values of pseudo F indicate well-separated, tight clusters.

The sensitivity and resistant phenotypes of each cell line for a given drug were defined by comparing the drug sensitivity measure between cell lines from any given cluster and the cell lines in all remaining clusters combined. We focus on $IC_{50}$ for clarity. For a given cluster k, cell lines were partitioned into two groups $j = \{Cluster-k, non-Cluster-k\}$. We then compared the mean $IC_{50}$ values of the two groups in a t test. The sign of the statistical test was used to define the phenotype as sensitive ($t > 0$) or resistant ($t < 0$). We accounted for multiple testing, by calculating the FDR-adjusted p-value for each drug. An FDR-adjusted p-value < 0.05 was considered significant.

*Supplementary data.* A summary of each cell line cluster, with information regarding tissue composition, molecular profile and drug profile, and other supplementary data for this article can be accessed from the publisher's website.

**Consent for publication.** All authors read and approved the final manuscript.

**Availability of data and materials.** All data analyzed during this study were retrieved and curated based on the pipeline published by Haibe-Kains *et al.*¹². Data are available from the CCLE website (http://www.broadinstitute) and GDSC website (http://www.cancerrxgene.org/downloads/).

**References**

1. Moreno, L. & Pearson, A. D. J. How can attrition rates be reduced in cancer drug discovery? *Expert Opin. Drug Discov.* **8**, 363–8 (2013).

2. Tauzin, B. More than 900 medicines and vaccines in clinical testing offer new hope in the fight against cancer. *DC Med. Dev. Cancer* **1–2** (2009).
3. Weinstein, J. N. et al. An information-intensive approach to the molecular pharmacology of cancer. *Science* **275**, 343–9 (1997).

4. Garnett, M. J. et al. Systematic identification of genomic markers of drug sensitivity in cancer cells. *Nature* **483**, 570–575 (2012).

5. Barretina, J., Caponigro, G. & Stransky, N. The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. *Nature* **483**, 603–607 (2012).

6. Chapman, P. B. et al. Improved survival with vemurafenib in melanoma with BRAF V600E mutation. *N. Engl. J. Med.* **364**, 2507–16 (2011).

7. Konecny, G. E. et al. Activity of the dual kinase inhibitor lapatinib (GW572016) against HER-2-overexpressing and trastuzumab-treated breast cancer cells. *Cancer Res.* **66**, 1630–9 (2006).

8. Heiser, L. M. et al. Subtype and pathway specific responses to anticancer compounds in breast cancer. *Proc. Natl. Acad. Sci.* **109**, 2724–2729 (2012).

9. von Minckwitz, G. et al. Definition and Impact of Pathologic Complete Response on Prognosis After Neoadjuvant Chemotherapy in Various Intrinsic Breast Cancer Subtypes. *J. Clin. Oncol.* **30**, 1796–1804 (2012).

10. Tsimberidou, A.-M. et al. Personalized Medicine in a Phase I Clinical Trials Program: The MD Anderson Cancer Center Initiative. *Clin. Cancer Res.* **18**, 6373–6383 (2012).

11. Hyman, D. M. et al. Vemurafenib in Multiple Nonmelanoma Cancers with BRAF V600 Mutations. *N. Engl. J. Med.* **373**, 726–36 (2015).

12. Haibe-Kains, B. et al. Inconsistency in large pharmacogenomic studies. *Nature* **504**, 389–93 (2013).

13. Monti, S., Tamayo, P., Mesirov, J. & Golub, T. Consensus Clustering - A resampling-based method for class discovery and visualization of gene expression microarray data Stefano. *Mach. Learn.* **52**, 1–34 (2003).

14. Szklarczyk, D. et al. The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. *Nucleic Acids Res.* **39**, D561–D568 (2011).

15. Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. *Proc. Natl. Acad. Sci. USA* **102**, 15545–50 (2005).

16. Quinlan, M. P., Quatela, S. E., Philips, M. R. & Settleman, J. Activated Kras, but Not Hras or Nras, May Initiate Tumors of Endodermal Origin via Stem Cell Expansion. *Mol. Cell. Biol.* **28**, 2659–2674 (2008).

17. Klijn, C. et al. A comprehensive transcriptional portrait of human cancer cell lines. *Nat. Biotechnol.*, doi:https://doi.org/10.1038/ nbt.3080 (2014).

18. Taube, J. H. et al. Core epithelial-to-mesenchymal transition interactome gene-expression signature is associated with claudin-low and metastolic breast cancer subtypes. *Proc. Natl. Acad. Sci. USA* **107**, 15449–15454 (2010).

19. Calinski, T. & Harabasz, J. A dendrite method for cluster analysis. *Commun. Stat. - Theory Methods* **3**, 1–27 (1974).

20. Fallahi-Sichani, M., Honarnejad, S., Heiser, L. M., Gray, J. W. & Sorger, P. K. Metrics other than potency reveal systematic variation in responses to cancer drugs. *Nat. Chem. Biol.* **9**, 708–14 (2013).

21. Greshock, J. et al. Molecular target class is predictive of *in vitro* response profile. *Cancer Res.* **70**, 3677–3686 (2010).

22. Haverty, P. M. et al. Reproducible pharmacogenomic profiling of cancer cell line panels. *Nature* **533**, 333–337 (2016).

23. Sontrop, H. M., Moerland, P. D., van den Ham, R., Reinders, M. J. & Verhaegh, W. F. A comprehensive sensitivity analysis of microarray breast cancer classification under feature variability. *BMC Bioinformatics* **10**, 389 (2009).

24. Cun, Y. & Fröhlich, H. Prognostic gene signatures for patient stratification in breast cancer - accuracy, stability and interpretability of gene selection approaches using prior knowledge on protein-protein interactions. *BMC Bioinformatics* **13**, 69 (2012).

25. Sanavia, T., Aiolli, F., Da San Martino, G., Bisognin, A. & Di Camillo, B. Improving biomarker list stability by integration of biological knowledge in the learning process. *BMC Bioinformatics* **13**, S22 (2012).

26. Bonsang-Kitzis, H. et al. Biological network-driven gene selection identifies a stromal immune module as a key determinant of triple-negative breast carcinoma prognosis. *Oncoimmunology* **37–41** (2015).

27. Lehmann, B. & Bauer, J. Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. *J. Clin. Investig* (2011).

28. Turner, N. C. & Reis-Filho, J. S. Tackling the Diversity of Triple-Negative Breast Cancer. *Clin. Cancer Res.* **19**, 6380–6388 (2013).

29. Liedtke, C. et al. Response to Neoadjuvant Therapy and Long-Term Survival in Patients With Triple-Negative Breast Cancer. *J. Clin. Oncol.* **26**, 1275–1281 (2008).

30. Cortazar, P. et al. Pathological complete response and long-term clinical benefit in breast cancer: the CTNeoBC pooled analysis. *Lancet* **384**, 164–172 (2014).

31. Wang, D., Boerner, S. A., Winkler, J. D. & LoRusso, P. M. Clinical experience of MEK inhibitors in cancer therapy. *Biochim. Biophys. Acta - Mol. Cell Res.* **1773**, 1248–1255 (2007).

32. Yeh, J. J. et al. KRAS/BRAF mutation status and ERK1/2 activation as biomarkers for MEK1/2 inhibitor therapy in colorectal cancer. *Mol. Cancer Ther.* **8**, 834–843 (2009).

33. Sawai, C. M. et al. Therapeutic Targeting of the Cyclin D3:CDK4/6 Complex in T Cell Leukemia. *Cancer Cell* **22**, 452–465 (2012).

34. Minami, T. et al. HER2 As Therapeutic Target for Overcoming ATP-Binding Cassette Transporter-Mediated Chemoresistance in Small Cell Lung Cancer. *Mol. Cancer Ther.* **11**, 830–841 (2012).

35. Ashburn, T. T. & Thor, K. B. Drug repositioning: identifying and developing new uses for existing drugs. *Nat. Rev. Drug Discov.* **3**, 673–683 (2004).

36. Van Allen, E. M. et al. Genomic Correlate of Exceptional Erlotinib Response in Head and Neck Squamous Cell Carcinoma. *JAMA Oncol.* **1**, 238–44 (2015).

37. Vokes, E. E. et al. A phase I study of erlotinib and bevacizumab for recurrent or metastatic squamous cell carcinoma of the head and neck (HNC). *J. Clin. Oncol.* **23**, 5504–5504 (2005).

38. Cohen, E. E. et al. Erlotinib and bevacizumab in patients with recurrent or metastatic squamous-cell carcinoma of the head and neck: a phase I/II study. *Lancet Oncol.* **10**, 247–257 (2009).

39. Pruthi, R. S. et al. A phase II trial of neoadjuvant erlotinib in patients with muscle-invasive bladder cancer undergoing radical cystectomy: Clinical and pathological results. *BJU Int.* **106**, 349–356 (2010).

40. Li, G. et al. Phase II study of concurrent chemoradiation in combination with erlotinib for locally advanced esophageal carcinoma. *Int. J. Radiat. Oncol. Biol. Phys.* **78**, 1407–1412 (2010).

41. Wainberg, Z. A. et al. Phase II trial of modified FOLFOX6 and erlotinib in patients with metastatic or advanced adenocarcinoma of the oesophagus and gastro-oesophageal junction. *Br. J. Cancer* **105**, 760–5 (2011).

42. Rees, M. G. et al. Correlating chemical sensitivity and basal gene expression reveals mechanism of action. *Nat Chem Biol* **12**, 109–116 (2016).

43. McCall, M. N., Bolstad, B. M. & Irizarry, R. A. Frozen robust multiarray analysis (fRMA). *Biostatistics* **11**, 242–253 (2010).

44. Irizarry, R. A. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. *Biostatistics* **4**, 249–264 (2003).

45. Li, Q., Birkbak, N. J., Gyorffy, B., Szallasi, Z. & Eklund, A. C. Jetset: selecting the optimal microarray probe set to represent a gene. *BMC Bioinformatics* **12**, 474 (2011).

46. Yu, G., Wang, L.-G., Han, Y. & He, Q.-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. *OMICS* **16**, 284–7 (2012).

47. Ritchie, M. E. et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. *Nucleic Acids Res.* **43**, e47–e47 (2015).
Acknowledgements

The authors thank Sergio Roman-Roman for reviewing the study and the manuscript. This work was supported by the Department of Translational Research of the Institut Curie and the Fondation ARC (Association pour la Recherche contre le Cancer Grant No. DOC20160604097). Funding was also obtained from the Site de Recherche Integrée en Cancérologie/Institut National du Cancer (Grant No. INCa-DGOS-4654).

Author Contributions

Conceptualization: F.R. P.N. Data curation: B.S. Formal analysis: B.S., J.A. Investigation: B.S., C.L., P.N. Methodology: B.S., P.N., C.L., P.G. Validation: A.P. Project administration: F.R., P.N. Supervision: F.R., P.N. Writing - original draft: B.S., A.S.H.P. All authors reviewed the manuscript.

Additional Information

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-017-14770-6.

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017