Bivona, Trever G. | Hieronymus, Haley | Parker, Joel | Chang, Kenneth | Taron, Miquel | Rosell, Rafael | Moonsamy, Philicia | Dahlman, Kimberly | Miller, Vincent A. | Costa, Carlota | Hannon, Gregory | Sawyers, Charles L.
Nature
2011;471(7339):523-526.
Human lung adenocarcinomas with activating mutations in EGFR (epidermal growth factor receptor) often respond to treatment with EGFRtyrosine kinase inhibitors(TKIs),butthe magnitude of tumour regression is variable and transient1,2. This heterogeneity in treatment response could result from genetic modifiers that regulate the degree to which tumour cells are dependent on mutant EGFR. Through a pooled RNA interference screen, we show that knockdown of FAS and several components of the NF-κB pathway specifically enhanced cell death induced by the EGFR TKI erlotinib in EGFR-mutant lung cancer cells. Activation of NF-κB through overexpression of c-FLIP or IKK (also known as CFLAR and IKBKB, respectively), or silencing of IκB (also known as NFKBIA), rescued EGFR-mutant lung cancer cells from EGFR TKI treatment. Genetic or pharmacologic inhibition of NF-κB enhanced erlotinib-induced apoptosis in erlotinib-sensitive and erlotinib-resistant EGFR-mutant lung cancer models. Increased expression of the NF-κB inhibitor IκB predicted for improved response and survival in EGFR-mutant lung cancer patients treated with EGFR TKI. These data identify NF-κB as a potential companion drug target, together with EGFR, in EGFR-mutant lung cancers and provide insight into the mechanisms by which tumour cells escape from oncogene dependence.
doi:10.1038/nature09870
PMCID: PMC3541675
PMID: 21430781
Introduction
In breast cancer, the basal-like subtype has high levels of genomic instability relative to other breast cancer subtypes with many basal-like-specific regions of aberration. There is evidence that this genomic instability extends to smaller scale genomic aberrations, as shown by a previously described micro-deletion event in the PTEN gene in the Basal-like SUM149 breast cancer cell line.
Methods
We sought to identify if small regions of genomic DNA copy number changes exist by using a high density, gene-centric Comparative Genomic Hybridizations (CGH) array on cell lines and primary tumors. A custom tiling array for CGH (244,000 probes, 200 bp tiling resolution) was created to identify small regions of genomic change, which was focused on previously identified basal-like-specific, and general cancer genes. Tumor genomic DNA from 94 patients and 2 breast cancer cell lines was labeled and hybridized to these arrays. Aberrations were called using SWITCHdna and the smallest 25% of SWITCHdna-defined genomic segments were called micro-aberrations (<64 contiguous probes, ∼ 15 kb).
Results
Our data showed that primary tumor breast cancer genomes frequently contained many small-scale copy number gains and losses, termed micro-aberrations, most of which are undetectable using typical-density genome-wide aCGH arrays. The basal-like subtype exhibited the highest incidence of these events. These micro-aberrations sometimes altered expression of the involved gene. We confirmed the presence of the PTEN micro-amplification in SUM149 and by mRNA-seq showed that this resulted in loss of expression of all exons downstream of this event. Micro-aberrations disproportionately affected the 5′ regions of the affected genes, including the promoter region, and high frequency of micro-aberrations was associated with poor survival.
Conclusion
Using a high-probe-density, gene-centric aCGH microarray, we present evidence of small-scale genomic aberrations that can contribute to gene inactivation. These events may contribute to tumor formation through mechanisms not detected using conventional DNA copy number analyses.
doi:10.1371/journal.pone.0051719
PMCID: PMC3524128
PMID: 23284754
Background
Next-generation sequencing technologies have become important tools for genome-wide studies. However, the quality scores that are assigned to each base have been shown to be inaccurate. If the quality scores are used in downstream analyses, these inaccuracies can have a significant impact on the results.
Results
Here we present ReQON, a tool that recalibrates the base quality scores from an input BAM file of aligned sequencing data using logistic regression. ReQON also generates diagnostic plots showing the effectiveness of the recalibration. We show that ReQON produces quality scores that are both more accurate, in the sense that they more closely correspond to the probability of a sequencing error, and do a better job of discriminating between sequencing errors and non-errors than the original quality scores. We also compare ReQON to other available recalibration tools and show that ReQON is less biased and performs favorably in terms of quality score accuracy.
Conclusion
ReQON is an open source software package, written in R and available through Bioconductor, for recalibrating base quality scores for next-generation sequencing data. ReQON produces a new BAM file with more accurate quality scores, which can improve the results of downstream analysis, and produces several diagnostic plots showing the effectiveness of the recalibration.
doi:10.1186/1471-2105-13-221
PMCID: PMC3447716
PMID: 22946927
Next-generation sequencing; Quality score; Recalibration; Bioinformatics; Bioconductor
Ellis, Matthew J. | Suman, Vera J. | Hoog, Jeremy | Lin, Li | Snider, Jacqueline | Prat, Aleix | Parker, Joel S. | Luo, Jingqin | DeSchryver, Katherine | Allred, D. Craig | Esserman, Laura J. | Unzeitig, Gary W. | Margenthaler, Julie | Babiera, Gildy V. | Marcom, P. Kelly | Guenther, Joseph M. | Watson, Mark A. | Leitch, Marilyn | Hunt, Kelly | Olson, John A.
Purpose
Preoperative aromatase inhibitor (AI) treatment promotes breast-conserving surgery (BCS) for estrogen receptor (ER) –positive breast cancer. To study this treatment option, responses to three AIs were compared in a randomized phase II neoadjuvant trial designed to select agents for phase III investigations.
Patients and Methods
Three hundred seventy-seven postmenopausal women with clinical stage II to III ER-positive (Allred score 6-8) breast cancer were randomly assigned to receive neoadjuvant exemestane, letrozole, or anastrozole. The primary end point was clinical response. Secondary end points included BCS, Ki67 proliferation marker changes, the Preoperative Endocrine Prognostic Index (PEPI), and PAM50-based intrinsic subtype analysis.
Results
On the basis of clinical response rates, letrozole and anastrozole were selected for further investigation; however, no other differences in surgical outcome, PEPI score, or Ki67 suppression were detected. The BCS rate for mastectomy-only patients at presentation was 51%. PAM50 analysis identified AI-unresponsive nonluminal subtypes (human epidermal growth factor receptor 2 enriched or basal-like) in 3.3% of patients. Clinical response and surgical outcomes were similar in luminal A (LumA) versus luminal B tumors; however, a PEPI of 0 (best prognostic group) was highest in the LumA subset (27.1% v 10.7%; P = .004).
Conclusion
Neoadjuvant AI treatment markedly improved surgical outcomes. Ki67 and PEPI data demonstrated that the three agents tested are biologically equivalent and therefore likely to have similar adjuvant activities. LumA tumors were more likely to have favorable biomarker characteristics after treatment; however, occasional paradoxical increases in Ki67 (12% of tumors with > 5% increase after therapy) suggest treatment-resistant cells, present in some LumA tumors, can be detected by post-treatment profiling.
doi:10.1200/JCO.2010.31.6950
PMCID: PMC3107749
PMID: 21555689
Dahlman, Kimberly Brown | Parker, Joel S. | Shamu, Tambudzai | Hieronymus, Haley | Chapinski, Caren | Carver, Brett | Chang, Kenneth | Hannon, Gregory J. | Sawyers, Charles L. | Campbell, Moray
There is significant need to identify novel prostate cancer drug targets because current hormone therapies eventually fail, leading to a drug-resistant and fatal disease termed castration-resistant prostate cancer. To functionally identify genes that, when silenced, decrease prostate cancer cell proliferation or induce cell death in combination with antiandrogens, we employed an RNA interference-based short hairpin RNA barcode screen in LNCaP human prostate cancer cells. We identified and validated four candidate genes (AKT1, PSMC1, STRADA, and TTK) that impaired growth when silenced in androgen receptor positive prostate cancer cells and enhanced the antiproliferative effects of antiandrogens. Inhibition of AKT with a pharmacologic inhibitor also induced apoptosis when combined with antiandrogens, consistent with recent evidence for PI3K and AR pathway crosstalk in prostate cancer cells. Recovery of hairpins targeting a known prostate cancer pathway validates the utility of shRNA library screening in prostate cancer as a broad strategy to identify new candidate drug targets.
doi:10.1371/journal.pone.0034414
PMCID: PMC3324507
PMID: 22509301
Bild, Andrea H | Parker, Joel S | Gustafson, Adam M | Acharya, Chaitanya R | Hoadley, Katherine A | Anders, Carey | Marcom, P Kelly | Carey, Lisa A | Potti, Anil | Nevins, Joseph R | Perou, Charles M
doi:10.1186/bcr2909
PMCID: PMC3236331
doi:10.1200/JCO.2010.28.9199
PMCID: PMC3055864
PMID: 21115855
Background
Recurrent head and neck squamous cell carcinoma (HNSCC) remains a difficult cancer to treat. Here, we describe a patient with HNSCC who had complete response to methotrexate (MTX) after progressing on multiple cytotoxic agents, cetuximab, and AMG-479 (monoclonal antibody against insulin-like growth factor-1 receptor [IGF-1R]).
Methods
The clinical information was collected by a retrospective medical record review under an Institutional Review Board–approved protocol. From 4 tumors and 2 normal mucosal epithelia, global gene expression, and IGF-1R and dihydrofolate reductase (DHFR) protein levels were determined.
Results
Effective target inhibition in the tumor was confirmed by the decreased protein levels of total and phospho-IGF-1R after treatment with AMG-479. Decreased level of DHFR and conversion of a gene expression profile associated with cetuximab-resistance to cetuximab-sensitivity were also observed.
Conclusion
This suggests that the combination of AMG- 479 and MTX or cetuximab may be a promising therapeutic approach in refractory HNSCC.
doi:10.1002/hed.21478
PMCID: PMC3111896
PMID: 20652976
IGF-1R inhibitor; cetuximab; head and neck squamous cell carcinoma; AMG-479; methotrexate
Weigman, Victor J. | Chao, Hann-Hsiang | Shabalin, Andrey A. | He, Xiaping | Parker, Joel S. | Nordgard, Silje H. | Grushko, Tatyana | Huo, Dezheng | Nwachukwu, Chika | Nobel, Andrew | Kristensen, Vessela N. | Børresen-Dale, Anne-Lise | Olopade, Olufunmilayo I. | Perou, Charles M.
Breast cancer is a heterogeneous disease with known expression-defined tumor subtypes. DNA copy number studies have suggested that tumors within gene expression subtypes share similar DNA Copy number aberrations (CNA) and that CNA can be used to further sub-divide expression classes. To gain further insights into the etiologies of the intrinsic subtypes, we classified tumors according to gene expression subtype and next identified subtype-associated CNA using a novel method called SWITCHdna, using a training set of 180 tumors and a validation set of 359 tumors. Fisher’s exact tests, Chi-square approximations, and Wilcoxon rank-sum tests were performed to evaluate differences in CNA by subtype. To assess the functional significance of loss of a specific chromosomal region, individual genes were knocked down by shRNA and drug sensitivity, and DNA repair foci assays performed. Most tumor subtypes exhibited specific CNA. The Basal-like subtype was the most distinct with common losses of the regions containing RB1, BRCA1, INPP4B, and the greatest overall genomic instability. One Basal-like subtype-associated CNA was loss of 5q11–35, which contains at least three genes important for BRCA1-dependent DNA repair (RAD17, RAD50, and RAP80); these genes were predominantly lost as a pair, or all three simultaneously. Loss of two or three of these genes was associated with significantly increased genomic instability and poor patient survival. RNAi knockdown of RAD17, or RAD17/RAD50, in immortalized human mammary epithelial cell lines caused increased sensitivity to a PARP inhibitor and carboplatin, and inhibited BRCA1 foci formation in response to DNA damage. These data suggest a possible genetic cause for genomic instability in Basal-like breast cancers and a biological rationale for the use of DNA repair inhibitor related therapeutics in this breast cancer subtype.
Electronic supplementary material
The online version of this article (doi:10.1007/s10549-011-1846-y) contains supplementary material, which is available to authorized users.
doi:10.1007/s10549-011-1846-y
PMCID: PMC3387500
PMID: 22048815
Basal-like breast cancer; Genome instability; BRCA1 pathway; Copy number aberration; Molecular subtypes; Array CGH
Nielsen, Torsten O. | Parker, Joel S | Leung, Samuel | Voduc, David | Ebbert, Mark | Vickery, Tammi | Davies, Sherri R. | Snider, Jacqueline | Stijleman, Inge J. | Reed, Jerry | Cheang, Maggie C.U. | Mardis, Elaine R. | Perou, Charles M. | Bernard, Philip S. | Ellis, Matthew J.
Purpose
To compare clinical, immunohistochemical and gene expression models of prognosis applicable to formalin-fixed, paraffin-embedded blocks in a large series of estrogen receptor positive breast cancers, from patients uniformly treated with adjuvant tamoxifen.
Methods
qRT-PCR assays for 50 genes identifying intrinsic breast cancer subtypes were completed on 786 specimens linked to clinical (median followup 11.7 years) and immunohistochemical (ER, PR, HER2, Ki67) data. Performance of predefined intrinsic subtype and Risk-Of-Relapse scores was assessed using multivariable Cox models and Kaplan-Meier analysis. Harrell’s C index was used to compare fixed models trained in independent data sets, including proliferation signatures.
Results
Despite clinical ER positivity, 10% of cases were assigned to non-Luminal subtypes. qRT-PCR signatures for proliferation genes gave more prognostic information than clinical assays for hormone receptors or Ki67. In Cox models incorporating standard prognostic variables, hazard ratios for breast cancer disease specific survival over the first 5 years of followup, relative to the most common Luminal A subtype, are 1.99 (95% CI: 1.09–3.64) for Luminal B, 3.65 (1.64–8.16) for HER2-enriched and 17.71 (1.71–183.33) for the basal like subtype. For node-negative disease, PAM50 qRT-PCR based risk assignment weighted for tumor size and proliferation identifies a group with >95% 10 yr survival without chemotherapy. In node positive disease, PAM50-based prognostic models were also superior.
Conclusion
The PAM50 gene expression test for intrinsic biological subtype can be applied to large series of formalin-fixed paraffin-embedded breast cancers, and gives more prognostic information than clinical factors and immunohistochemistry using standard cutpoints.
doi:10.1158/1078-0432.CCR-10-1282
PMCID: PMC2970720
PMID: 20837693
Background
Globally, gastric cancer is the second most common cause of cancer-related death, with the majority of the health burden borne by economically less-developed countries.
Methods
Here, we report a genetic characterization of 50 gastric adenocarcinoma samples, using affymetrix SNP arrays and Illumina mRNA expression arrays as well as Illumina sequencing of the coding regions of 384 genes belonging to various pathways known to be altered in other cancers.
Results
Genetic alterations were observed in the WNT, Hedgehog, cell cycle, DNA damage and epithelial-to-mesenchymal-transition pathways.
Conclusions
The data suggests targeted therapies approved or in clinical development for gastric carcinoma would be of benefit to ~22% of the patients studied. In addition, the novel mutations detected here, are likely to influence clinical response and suggest new targets for drug discovery.
doi:10.1186/1479-5876-9-119
PMCID: PMC3152520
PMID: 21781349
Background
Genomic tests are available to predict breast cancer recurrence and to guide clinical decision making. These predictors provide recurrence risk scores along with a measure of uncertainty, usually a confidence interval. The confidence interval conveys random error and not systematic bias. Standard tumor sampling methods make this problematic, as it is common to have a substantial proportion (typically 30-50%) of a tumor sample comprised of histologically benign tissue. This "normal" tissue could represent a source of non-random error or systematic bias in genomic classification.
Methods
To assess the performance characteristics of genomic classification to systematic error from normal contamination, we collected 55 tumor samples and paired tumor-adjacent normal tissue. Using genomic signatures from the tumor and paired normal, we evaluated how increasing normal contamination altered recurrence risk scores for various genomic predictors.
Results
Simulations of normal tissue contamination caused misclassification of tumors in all predictors evaluated, but different breast cancer predictors showed different types of vulnerability to normal tissue bias. While two predictors had unpredictable direction of bias (either higher or lower risk of relapse resulted from normal contamination), one signature showed predictable direction of normal tissue effects. Due to this predictable direction of effect, this signature (the PAM50) was adjusted for normal tissue contamination and these corrections improved sensitivity and negative predictive value. For all three assays quality control standards and/or appropriate bias adjustment strategies can be used to improve assay reliability.
Conclusions
Normal tissue sampled concurrently with tumor is an important source of bias in breast genomic predictors. All genomic predictors show some sensitivity to normal tissue contamination and ideal strategies for mitigating this bias vary depending upon the particular genes and computational methods used in the predictor.
doi:10.1186/1755-8794-4-54
PMCID: PMC3151208
PMID: 21718502
biomarker validation; genomic assays; breast cancer; normal tissue; bias
The ability to predict metastatic potential could be of great clinical importance, however, it is uncertain if predicting metastasis to specific vital organs is feasible. As a first step in evaluating metastatic predictions, we analyzed multiple primary tumors and metastasis pairs and determined that >90% of 298 gene expression signatures were found to be similarly expressed between matched pairs of tumors and metastases; therefore, primary tumors may be a good predictor of metastatic propensity. Next, using a dataset of >1,000 human breast tumor gene expression microarrays we determined that HER2-enriched subtype tumors aggressively spread to the liver, while basal-like and claudin-low subtypes colonize the brain and lung. Correspondingly, brain and lung metastasis signatures, along with embryonic stem cell, tumor initiating cell, and hypoxia signatures, were also strongly expressed in the basal-like and claudin-low tumors. Interestingly, low “Differentiation Scores,” or high expression of the aforementioned signatures, further predicted for brain and lung metastases. In total, these data identify that depending upon the organ of relapse, a combination of gene expression signatures most accurately predicts metastatic behavior.
Electronic supplementary material
The online version of this article (doi:10.1007/s10549-011-1619-7) contains supplementary material, which is available to authorized users.
doi:10.1007/s10549-011-1619-7
PMCID: PMC3303043
PMID: 21671017
Breast cancer; Gene expression; Intrinsic subtype; Metastasis; Microarray
Dunbier, Anita K. | Anderson, Helen | Ghazoui, Zara | Folkerd, Elizabeth J. | A'Hern, Roger | Crowder, Robert J. | Hoog, Jeremy | Smith, Ian E. | Osin, Peter | Nerurkar, Ashutosh | Parker, Joel S. | Perou, Charles M. | Ellis, Matthew J. | Dowsett, Mitch
Purpose
To determine whether plasma estradiol (E2) levels are related to gene expression in estrogen receptor (ER)–positive breast cancers in postmenopausal women.
Materials and Methods
Genome-wide RNA profiles were obtained from pretreatment core-cut tumor biopsies from 104 postmenopausal patients with primary ER-positive breast cancer treated with neoadjuvant anastrozole. Pretreatment plasma E2 levels were determined by highly sensitive radioimmunoassay. Genes were identified for which expression was correlated with pretreatment plasma E2 levels. Validation was performed in an independent set of 73 ER-positive breast cancers.
Results
The expression of many known estrogen-responsive genes and gene sets was highly significantly associated with plasma E2 levels (eg, TFF1/pS2, GREB1, PDZK1 and PGR; P < .005). Plasma E2 explained 27% of the average expression of these four average estrogen-responsive genes (ie, AvERG; r = 0.51; P < .0001), and a standardized mean of plasma E2 levels and ER transcript levels explained 37% (r, 0.61). These observations were validated in an independent set of 73 ER-positive tumors. Exploratory analysis suggested that addition of the nuclear coregulators in a multivariable analysis with ER and E2 levels might additionally improve the relationship with the AvERG. Plasma E2 and the standardized mean of E2 and ER were both significantly correlated with 2-week Ki67, a surrogate marker of clinical outcome (r = −0.179; P = .05; and r = −0.389; P = .0005, respectively).
Conclusion
Plasma E2 levels are significantly associated with gene expression of ER-positive breast cancers and should be considered in future genomic studies of ER-positive breast cancer. The AvERG is a new experimental tool for the study of putative estrogenic stimuli of breast cancer.
doi:10.1200/JCO.2009.23.9616
PMCID: PMC2834467
PMID: 20124184
Background
Multiple breast cancer gene expression profiles have been developed that appear to provide similar abilities to predict outcome and may outperform clinical-pathologic criteria; however, the extent to which seemingly disparate profiles provide additive prognostic information is not known, nor do we know whether prognostic profiles perform equally across clinically defined breast cancer subtypes. We evaluated whether combining the prognostic powers of standard breast cancer clinical variables with a large set of gene expression signatures could improve on our ability to predict patient outcomes.
Methods
Using clinical-pathological variables and a collection of 323 gene expression "modules", including 115 previously published signatures, we build multivariate Cox proportional hazards models using a dataset of 550 node-negative systemically untreated breast cancer patients. Models predictive of pathological complete response (pCR) to neoadjuvant chemotherapy were also built using this approach.
Results
We identified statistically significant prognostic models for relapse-free survival (RFS) at 7 years for the entire population, and for the subgroups of patients with ER-positive, or Luminal tumors. Furthermore, we found that combined models that included both clinical and genomic parameters improved prognostication compared with models with either clinical or genomic variables alone. Finally, we were able to build statistically significant combined models for pathological complete response (pCR) predictions for the entire population.
Conclusions
Integration of gene expression signatures and clinical-pathological factors is an improved method over either variable type alone. Highly prognostic models could be created when using all patients, and for the subset of patients with lymph node-negative and ER-positive breast cancers. Other variables beyond gene expression and clinical-pathological variables, like gene mutation status or DNA copy number changes, will be needed to build robust prognostic models for ER-negative breast cancer patients. This combined clinical and genomics model approach can also be used to build predictors of therapy responsiveness, and could ultimately be applied to other tumor types.
doi:10.1186/1755-8794-4-3
PMCID: PMC3025826
PMID: 21214954
Fannin, Rick D. | Russo, Mark | O’Connell, Thomas M. | Gerrish, Kevin | Winnike, Jason H. | Macdonald, Jeffrey | Newton, Jack | Malik, Shahid | Sieber, Stella O. | Parker, Joel | Shah, Ruchir | Zhou, Tong | Watkins, Paul B. | Paules, Richard S.
The diagnosis and management of drug-induced liver injury (DILI) is hindered by the limited utility of traditional clinical chemistries. It has recently been shown that hepatotoxicants can produce compound-specific changes in the peripheral blood (PB) transcriptome in rodents, suggesting the blood transcriptome might provide new biomarkers of DILI. To investigate in humans, we used DNA microarrays as well as serum metabolomic methods to characterize changes in the transcriptome and metabolome in serial PB samples obtained from 6 healthy adults treated with a 4 g bolus dose of acetaminophen (APAP) and from 3 receiving placebo. Treatment did not cause liver injury as assessed by traditional liver chemistries. However, 48 hours after exposure, treated subjects showed marked down-regulation of genes involved in oxidative phosphorylation/mitochondrial function that was not observed in the placebos (p <1.66E-19). The magnitude of down-regulation was positively correlated with the percent of APAP converted to the reactive metabolite NAPQI (r = 0.739; p=0.058). In addition, unbiased analysis of the serum metabolome revealed an increase in serum lactate from 24 to 72 hours post dosing in the treated subjects alone (p<0.005). Similar PB transcriptome changes were observed in human overdose patients and rats receiving toxic doses.
Conclusion
The single 4 gm APAP dose produced a transcriptome signature in PB cells characterized by down regulation of oxidative phosphorylation genes accompanied by increased serum lactate. Similar gene expression changes were observed in rats and several patients after consuming hepatotoxic doses of APAP. The timing of the changes and the correlation with NAPQI production are consistent with mechanisms known to underlie APAP hepatoxicity. These studies support the further exploration of the blood transcriptome for biomarkers of DILI.
doi:10.1002/hep.23330
PMCID: PMC2925683
PMID: 19918972
hepatotoxicity; microarray; biomarker; surrogate; mitochondria
Mavrakis, Konstantinos J. | Wolfe, Andrew L. | Oricchio, Elisa | Palomero, Teresa | de Keersmaecker, Kim | McJunkin, Katherine | Zuber, Johannes | James, Taneisha | Chang, Kenneth | Khan, Aly A. | Leslie, Christina S. | Parker, Joel S. | Paddison, Patrick J. | Tam, Wayne | Ferrando, Adolfo | Wendel, Hans-Guido
doi:10.1038/ncb2037
PMCID: PMC2989719
PMID: 20190740
By virtue of their accumulated genetic alterations, tumor cells may acquire vulnerabilities that create opportunities for therapeutic intervention. We have devised a massively parallel strategy for screening short hairpin RNA (shRNA) collections for stable loss-of-function phenotypes. We assayed from 6000 to 20,000 shRNAs simultaneously to identify genes important for the proliferation and survival of five cell lines derived from human mammary tissue. Lethal shRNAs common to these cell lines targeted many known cell-cycle regulatory networks. Cell line–specific sensitivities to suppression of protein complexes and biological pathways also emerged, and these could be validated by RNA interference (RNAi) and pharmacologically. These studies establish a practical platform for genome-scale screening of complex phenotypes in mammalian cells and demonstrate that RNAi can be used to expose genotype-specific sensitivities.
doi:10.1126/science.1149185
PMCID: PMC2981861
PMID: 18239125
Background
The transcription factor c-Myb has been well characterized as an oncogene in several human tumor types, and its expression in the hematopoietic stem/progenitor cell population is essential for proper hematopoiesis. However, the role of c-Myb in mammopoeisis and breast tumorigenesis is poorly understood, despite its high expression in the majority of breast cancer cases (60–80%).
Methodology/Principal Findings
We find that c-Myb high expression in human breast tumors correlates with the luminal/ER+ phenotype and a good prognosis. Stable RNAi knock-down of endogenous c-Myb in the MCF7 luminal breast tumor cell line increased tumorigenesis both in vitro and in vivo, suggesting a possible tumor suppressor role in luminal breast cancer. We created a mammary-derived c-Myb expression signature, comprised of both direct and indirect c-Myb target genes, and found it to be highly correlated with a published mature luminal mammary cell signature and least correlated with a mammary stem/progenitor lineage gene signature.
Conclusions/Significance
These data describe, for the first time, a possible tumor suppressor role for the c-Myb proto-oncogene in breast cancer that has implications for the understanding of luminal tumorigenesis and for guiding treatment.
doi:10.1371/journal.pone.0013073
PMCID: PMC2951337
PMID: 20949095
Introduction
In breast cancer, gene expression analyses have defined five tumor subtypes (luminal A, luminal B, HER2-enriched, basal-like and claudin-low), each of which has unique biologic and prognostic features. Here, we comprehensively characterize the recently identified claudin-low tumor subtype.
Methods
The clinical, pathological and biological features of claudin-low tumors were compared to the other tumor subtypes using an updated human tumor database and multiple independent data sets. These main features of claudin-low tumors were also evaluated in a panel of breast cancer cell lines and genetically engineered mouse models.
Results
Claudin-low tumors are characterized by the low to absent expression of luminal differentiation markers, high enrichment for epithelial-to-mesenchymal transition markers, immune response genes and cancer stem cell-like features. Clinically, the majority of claudin-low tumors are poor prognosis estrogen receptor (ER)-negative, progesterone receptor (PR)-negative, and epidermal growth factor receptor 2 (HER2)-negative (triple negative) invasive ductal carcinomas with a high frequency of metaplastic and medullary differentiation. They also have a response rate to standard preoperative chemotherapy that is intermediate between that of basal-like and luminal tumors. Interestingly, we show that a group of highly utilized breast cancer cell lines, and several genetically engineered mouse models, express the claudin-low phenotype. Finally, we confirm that a prognostically relevant differentiation hierarchy exists across all breast cancers in which the claudin-low subtype most closely resembles the mammary epithelial stem cell.
Conclusions
These results should help to improve our understanding of the biologic heterogeneity of breast cancer and provide tools for the further evaluation of the unique biology of claudin-low tumors and cell lines.
doi:10.1186/bcr2635
PMCID: PMC3096954
PMID: 20813035
Crowder, Robert J. | Phommaly, Chanpheng | Tao, Yu | Hoog, Jeremy | Luo, Jingqin | Perou, Charles M. | Parker, Joel S. | Miller, Melinda A. | Huntsman, David G. | Lin, Li | Snider, Jacqueline | Davies, Sherri R. | Olson, John A. | Watson, Mark A. | Saporita, Anthony | Weber, Jason D. | Ellis, Matthew J.
Several phosphoinositide-3-kinase (PI3K) catalytic subunit inhibitors are currently in clinical trial. We therefore sought to examine relationships between pharmacological inhibition and somatic mutations in PI3K catalytic subunits in ER+ breast cancer, where these mutations are particularly common. RNA interference (RNAi) was used to determine the effect of selective inhibition of PI3K catalytic subunits, p110α and p110β, in ER+ breast cancer cells harboring either mutation (PIK3CA) or gene amplification (PIK3CB). p110α RNAi inhibited growth and promoted apoptosis in all tested ER+ breast cancer cells under estrogen deprived-conditions, whereas p110β RNAi only affected cells harboring PIK3CB amplification. Moreover, dual p110α/p110β inhibition potentiated these effects. In addition, treatment with the clinical grade PI3K catalytic subunit inhibitor BEZ235 also promoted apoptosis in ER+ breast cancer cells. Importantly, estradiol suppressed apoptosis induced by both gene knockdowns and by BEZ235 treatment. Our results suggest that PI3K inhibitors should target both p110α and p110β catalytic subunits, whether wild-type or mutant, and be combined with endocrine therapy for maximal efficacy when treating ER+ breast cancer.
doi:10.1158/0008-5472.CAN-08-4450
PMCID: PMC2811393
PMID: 19366795
breast cancer; estrogen receptor; PI3 kinase; endocrine therapy; synthetic lethality
Background
The goal of this study was to understand gene expression signatures of hepatocellular carcinoma (HCC) recurrence in subjects with hepatitis C virus (HCV) infection. Recurrence-free survival (RFS) following curative resection of HCC in subjects with HCV is highly variable. Traditional clinico-pathological endpoints are recognized as weak predictors of RFS. It has been suggested that gene expression profiling of HCC and nontumoral liver tissue may improve prediction of RFS, aid in understanding of the underlying liver disease, and guide individualized patient management. Frozen samples of the tumors and nontumoral liver were obtained from 47 subjects with HCV-associated HCC. Additional nontumoral liver samples were obtained from HCV-free subjects with metastatic liver tumors. Gene expression profiling data was used to determine the molecular signature of HCV-associated HCC and to develop a predictor of RFS.
Results
The molecular profile of the HCV-associated HCC confirmed central roles for MYC and TGFβ1 in liver tumor development. Gene expression in tumors was found to have poor predictive power with regards to RFS, but analysis of nontumoral tissues yielded a strong predictor for RFS in late-recurring (>1 year) subjects. Importantly, nontumoral tissue-derived gene expression predictor of RFS was highly significant in both univariable and multivariable Cox proportional hazard model analyses.
Conclusions
Microarray analysis of the nontumoral tissues from subjects with HCV-associated HCC delivers novel molecular signatures of RFS, especially among the late-recurrence subjects. The gene expression predictor may hold important insights into the pathobiology of HCC recurrence and de novo tumor formation in cirrhotic patients.
doi:10.1186/1476-4598-9-74
PMCID: PMC2856554
PMID: 20380719
Parker, Joel S. | Mullins, Michael | Cheang, Maggie C.U. | Leung, Samuel | Voduc, David | Vickery, Tammi | Davies, Sherri | Fauron, Christiane | He, Xiaping | Hu, Zhiyuan | Quackenbush, John F. | Stijleman, Inge J. | Palazzo, Juan | Marron, J.S. | Nobel, Andrew B. | Mardis, Elaine | Nielsen, Torsten O. | Ellis, Matthew J. | Perou, Charles M. | Bernard, Philip S.
Purpose
To improve on current standards for breast cancer prognosis and prediction of chemotherapy benefit by developing a risk model that incorporates the gene expression–based “intrinsic” subtypes luminal A, luminal B, HER2-enriched, and basal-like.
Methods
A 50-gene subtype predictor was developed using microarray and quantitative reverse transcriptase polymerase chain reaction data from 189 prototype samples. Test sets from 761 patients (no systemic therapy) were evaluated for prognosis, and 133 patients were evaluated for prediction of pathologic complete response (pCR) to a taxane and anthracycline regimen.
Results
The intrinsic subtypes as discrete entities showed prognostic significance (P = 2.26E-12) and remained significant in multivariable analyses that incorporated standard parameters (estrogen receptor status, histologic grade, tumor size, and node status). A prognostic model for node-negative breast cancer was built using intrinsic subtype and clinical information. The C-index estimate for the combined model (subtype and tumor size) was a significant improvement on either the clinicopathologic model or subtype model alone. The intrinsic subtype model predicted neoadjuvant chemotherapy efficacy with a negative predictive value for pCR of 97%.
Conclusion
Diagnosis by intrinsic subtype adds significant prognostic and predictive information to standard parameters for patients with breast cancer. The prognostic properties of the continuous risk score will be of value for the management of node-negative breast cancers. The subtypes and risk score can also be used to assess the likelihood of efficacy from neoadjuvant chemotherapy.
doi:10.1200/JCO.2008.18.1370
PMCID: PMC2667820
PMID: 19204204
Haugen, Astrid C. | Di Prospero, Nicholas A. | Parker, Joel S. | Fannin, Rick D. | Chou, Jeff | Meyer, Joel N. | Halweg, Christopher | Collins, Jennifer B. | Durr, Alexandra | Fischbeck, Kenneth | Van Houten, Bennett | Pearson, Christopher E.
The neurodegenerative disease Friedreich's ataxia (FRDA) is the most common autosomal-recessively inherited ataxia and is caused by a GAA triplet repeat expansion in the first intron of the frataxin gene. In this disease, transcription of frataxin, a mitochondrial protein involved in iron homeostasis, is impaired, resulting in a significant reduction in mRNA and protein levels. Global gene expression analysis was performed in peripheral blood samples from FRDA patients as compared to controls, which suggested altered expression patterns pertaining to genotoxic stress. We then confirmed the presence of genotoxic DNA damage by using a gene-specific quantitative PCR assay and discovered an increase in both mitochondrial and nuclear DNA damage in the blood of these patients (p<0.0001, respectively). Additionally, frataxin mRNA levels correlated with age of onset of disease and displayed unique sets of gene alterations involved in immune response, oxidative phosphorylation, and protein synthesis. Many of the key pathways observed by transcription profiling were downregulated, and we believe these data suggest that patients with prolonged frataxin deficiency undergo a systemic survival response to chronic genotoxic stress and consequent DNA damage detectable in blood. In conclusion, our results yield insight into the nature and progression of FRDA, as well as possible therapeutic approaches. Furthermore, the identification of potential biomarkers, including the DNA damage found in peripheral blood, may have predictive value in future clinical trials.
Author Summary
Friedreich's ataxia is an inherited disease that causes progressive damage to the nervous system and affects the muscles and heart. The disease is caused by a defect in the frataxin gene, which is involved in iron homeostasis and likely protects against reactive oxygen species. In order to identify mechanisms involved in the nature and progression of the disease, we performed transcriptional profiling and measurements of mitochondrial and nuclear DNA damage on blood cells from FRDA patients. Transcriptional profiling was performed on blood samples from a cohort of 28 children compared to a control group. These data were then validated with a cohort of 14 adults with FRDA compared to a second independent control group. DNA damage was assessed on the blood samples from the 28 FRDA children, plus an additional 19 affected children, by quantitative PCR (QPCR). Transcriptional profiling revealed changes in gene expression consistent with the presence of genotoxic stress in FRDA patients. This finding was supported by the direct evidence that FRDA patients accumulated significantly higher levels of mitochondrial and nuclear DNA damage as compared to controls. The identification of potential biomarkers, including the DNA damage found in peripheral blood, may help identify therapeutic approaches for this devastating disease.
doi:10.1371/journal.pgen.1000812
PMCID: PMC2799513
PMID: 20090835
Cohen, Ezra Eddy Wyssam | Zhu, Hongyan | Lingen, Mark W. | Martin, Leslie E. | Kuo, Wen-Liang | Choi, Eugene A. | Kocherginsky, Masha | Parker, Joel S. | Chung, Christine H. | Rosner, Marsha Rich
Protein Kinase C alpha (PKCα) has been implicated in cancer but the mechanism is largely unknown. Here we show that PKCα promotes head and neck squamous cell carcinoma (SCCHN) by a feed forward network leading to cell cycle deregulation. PKCα inhibitors decrease proliferation in SCCHN cell lines and xenografted tumors. PKCα inhibition or depletion in tumor cells decreases DNA synthesis by suppressing ERK phosphorylation and cyclin E synthesis. Additionally, PKCα down-regulates miR-15a, a microRNA that directly inhibits protein synthesis of cyclin E as well as other cell cycle regulators. Furthermore, both PKCα and cyclin E protein expression are increased in primary tumors, and PKCα inversely correlates with miR15a expression in primary tumors. Finally, PKCα is associated with poor prognosis in SCCHN. These results identify PKCα as a key regulator of HNSCC tumor cell growth by a mechanism involving activation of MAP kinase, an initiator of the cell cycle, and suppression of miR-15a, an inhibitor of DNA synthesis. Although the specific components may be different, this type of feed forward loop network, consisting of a stimulus that activates a positive signal and removes a negative brake, is likely to be a general one that enables induction of DNA synthesis by a variety of growth or oncogenic stimuli.
doi:10.1158/0008-5472.CAN-08-0377
PMCID: PMC2746005
PMID: 19117988